title: Evolving Evaluation Functions for Collectible Card Game AI

abstract: In this work, we presented a study regarding two important aspects of evolving feature-based game evaluation functions: the choice of genome representation and the choice of opponent used to test the model. We compared three representations. One simpler and more limited, based on a vector of weights that are used in a linear combination of predefined game features. And two more complex, based on binary and n-ary trees. On top of this test, we also investigated the influence of fitness defined as a simulation-based function that: plays against a fixed weak opponent, plays against a fixed strong opponent, and plays against the best individual from the previous population. For a testbed, we have chosen a recently popular domain of digital collectible card games. We encoded our experiments in a programming game, Legends of Code and Magic, used in Strategy Card Game AI Competition. However, as the problems stated are of general nature we are convinced that our observations are applicable in the other domains as well.


title: Ultraviolet-Infrared Bounds and Minimum Coupling in Effective Field Theories

abstract: We provide a simple new argument for a lower bound on the coupling of a $U(1)$ gauge interaction in an effective field theory (EFT), originally obtained from the Weak Gravity Conjecture. Our argument employs basic principles of quantum mechanical energy-time uncertainty and Lorentz invariance, plus Bekenstein's entropy bound on the ultraviolet (UV) and infrared (IR) scales of an EFT. We show that using an alternative UV-IR relation based on the Cohen-Kaplan-Nelson (CKN) bound results in a stronger lower bound on the $U(1)$ coupling, consistent with the more stringent nature of the CKN relation. Applicability of our reasoning to other interactions is briefly discussed. We also slightly extend the CKN bound, by accounting for the effective degrees of freedom, and examine some of its phenomenological implications.


title: Archives and AI: An Overview of Current Debates and Future Perspectives

abstract: The digital transformation is turning archives, both old and new, into data. As a consequence, automation in the form of artificial intelligence techniques is increasingly applied both to scale traditional recordkeeping activities, and to experiment with novel ways to capture, organise and access records. We survey recent developments at the intersection of Artificial Intelligence and archival thinking and practice. Our overview of this growing body of literature is organised through the lenses of the Records Continuum model. We find four broad themes in the literature on archives and artificial intelligence: theoretical and professional considerations, the automation of recordkeeping processes, organising and accessing archives, and novel forms of digital archives. We conclude by underlining emerging trends and directions for future work, which include the application of recordkeeping principles to the very data and processes which power modern artificial intelligence, and a more structural, yet critically-aware, integration of artificial intelligence into archival systems and practice.


title: An $X(3872)$-like peak in $B\to (J/\psi\pi^+\pi^-) K\pi$ due to triangle singularity

abstract: Triangle mechanisms for $B^0\to (J/\psi\pi^+\pi^-) K^+\pi^-$ are studied. Experimentally, an $X(3872)$ peak has been observed in this process. When the final $(J/\psi\pi^+\pi^-)\pi$ invariant mass is around the $D^*\bar D^*$ threshold, one of the triangle mechanisms causes a triangle singularity and generates a sharp $X(3872)$-like peak in the $J/\psi\pi^+\pi^-$ invariant mass distribution. The Breit-Wigner mass and width fitted to the spectrum are 3871.68 MeV (a few keV above the $D^{*0}\bar{D}^0$ threshold) and $\sim$0.4 MeV, respectively.These Breit-Wigner parameters hardly depends on a choice of the model parameters. Comparing with the precisely measured $X(3872)$ mass and width, $3871.69\pm 0.17$ MeV and $< 1.2$ MeV, the agreement is remarkable. When studying the $X(3872)$ signal from this process, this non-resonant contribution has to be understood in advance. We also study a charge analogous process $B^0\to (J/\psi\pi^0\pi^-) K^+\pi^0$. A similar triangle singularity exists and generates an $X^-(3876)$-like peak.


title: Iterated learning for emergent systematicity in VQA

abstract: Although neural module networks have an architectural bias towards compositionality, they require gold standard layouts to generalize systematically in practice. When instead learning layouts and modules jointly, compositionality does not arise automatically and an explicit pressure is necessary for the emergence of layouts exhibiting the right structure. We propose to address this problem using iterated learning, a cognitive science theory of the emergence of compositional languages in nature that has primarily been applied to simple referential games in machine learning. Considering the layouts of module networks as samples from an emergent language, we use iterated learning to encourage the development of structure within this language. We show that the resulting layouts support systematic generalization in neural agents solving the more complex task of visual question-answering. Our regularized iterated learning method can outperform baselines without iterated learning on SHAPES-SyGeT (SHAPES Systematic Generalization Test), a new split of the SHAPES dataset we introduce to evaluate systematic generalization, and on CLOSURE, an extension of CLEVR also designed to test systematic generalization. We demonstrate superior performance in recovering ground-truth compositional program structure with limited supervision on both SHAPES-SyGeT and CLEVR.


title: Upper bounds on the average number of colors in the non-equivalent colorings of a graph

abstract: A coloring of a graph is an assignment of colors to its vertices such that adjacent vertices have different colors. Two colorings are equivalent if they induce the same partition of the vertex set into color classes. Let $\mathcal{A}(G)$ be the average number of colors in the non-equivalent colorings of a graph $G$. We give a general upper bound on $\mathcal{A}(G)$ that is valid for all graphs $G$ and a more precise one for graphs $G$ of order $n$ and maximum degree $\Delta(G)\in \{1,2,n-2\}$.


title: The color dipole model bounds with the gluon-gluon recombination correction

abstract: We present nonlinear (NL) and higher twist (HT) corrections to the color dipole model (CDM) bounds at low values of $x$ and $Q^{2}$ using the parameterization method. Consistency between the bounds at this region describe that a transition from the linear to the nonlinear behavior is dependence on the behavior of the gluon distribution function. The parameters in the color dipole model are comparable with the color dipole bounds at low values of $Q^{2}$. Consequently, the obtained reduced cross sections at low and moderate $Q^{2}$ values due to the NL+HT effects show a good agreement with the H1 data.


title: Complex pressure-temperature structural phase diagram of honeycomb iridate Cu$_2$IrO$_3$

abstract: $\mathrm{Cu_2IrO_3}$ is among the newest layered honeycomb iridates and a promising candidate to harbor a Kitaev quantum spin liquid state. Here, we investigate the pressure and temperature dependence of its structure through a combination of powder x-ray diffraction and x-ray absorption fine structure measurements, as well as $ab$-$initio$ evolutionary structure search. At ambient pressure, we revise the previously proposed $C2/c$ solution with a related but notably more stable $P2_1/c$ structure. Pressures below 8 GPa drive the formation of Ir-Ir dimers at both ambient and low temperatures, similar to the case of $\mathrm{Li_2IrO_3}$. At higher pressures, the structural evolution dramatically depends on temperature. A large discontinuous reduction of the Ir honeycomb interplanar distance is observed around 15 GPa at room temperature, likely driven by a collapse of the O-Cu-O dumbbells. At 15 K, pressures beyond 20 GPa first lead to an intermediate phase featuring a continuous reduction of the interplanar distance, which then collapses at 30 GPa across yet another phase transition. However, the resulting structure around 40 GPa is not the same at room and low temperatures. Remarkably, the reduction in interplanar distance leads to an apparent healing of the stacking faults at room temperature, but not at 15 K. Possible implications on the evolution of electronic structure of $\mathrm{Cu_2IrO_3}$ with pressure are discussed.


title: The Resolved Sunyaev--Zel'dovich Profiles of Nearby Galaxy Groups

abstract: Much of the baryons in galaxy groups are thought to have been driven out to large distances ($\gtrsim$$R_{500}$) by feedback, but there are few constraining observations of this extended gas. This work presents the resolved Sunyaev--Zel'dovich (SZ) profiles for a stacked sample of 10 nearby galaxy groups within the mass range log$_{10}(M_{500}[M_{\odot}]) = 13.6 -13.9$. We measured the SZ profiles using the publicly available $y$-map from the Planck Collaboration as well as our own $y$-maps constructed from more recent versions of $Planck$ data. The $y$-map extracted from the latest data release yielded a significant SZ detection out to 3 $R_{500}$. In addition, the stacked profile from these data were consistent with simulations that included AGN feedback. Our best-fit model using the latest $Planck$ data suggested a baryon fraction $\approx 5.6\%$ within $R_{500}$. This is significantly lower than the cosmic value of $\approx 16\%$, supporting the idea that baryons have been driven to large radii by AGN feedback. Lastly, we discovered a significant ($\sim 3\sigma$) "bump" feature near $\sim 2$ $R_{500}$ that is most likely the signature of internal accretion shocks.


title: Combining Broad and Narrow Case Definitions in Matched Case-Control Studies: Firearms in the Home and Suicide Risk

abstract: Does having firearms in the home increase suicide risk? To test this hypothesis, a matched case-control study can be performed, in which suicide case subjects are compared to living controls who are similar in observed covariates in terms of their retrospective exposure to firearms at home. In this application, cases can be defined using a broad case definition (suicide) or a narrow case definition (suicide occurred at home). The broad case definition offers a larger number of cases but the narrow case definition may offer a larger effect size. Moreover, restricting to the narrow case definition may introduce selection bias (i.e., bias due to selecting samples based on characteristics affected by the treatment) because exposure to firearms in the home may affect the location of suicide and thus the type of a case a subject is. We propose a new sensitivity analysis framework for combining broad and narrow case definitions in matched case-control studies, that considers the unmeasured confounding bias and selection bias simultaneously. We develop a valid randomization-based testing procedure using only the narrow case matched sets when the effect of the unmeasured confounder on receiving treatment and the effect of the treatment on case definition among the always-cases are controlled by sensitivity parameters. We then use the Bonferroni method to combine the testing procedures using the broad and narrow case definitions. With the proposed methods, we find robust evidence that having firearms at home increases suicide risk.


title: Context-aware demand prediction in bike sharing systems: incorporating spatial, meteorological and calendrical context

abstract: Bike sharing demand is increasing in large cities worldwide. The proper functioning of bike-sharing systems is, nevertheless, dependent on a balanced geographical distribution of bicycles throughout a day. In this context, understanding the spatiotemporal distribution of check-ins and check-outs is key for station balancing and bike relocation initiatives. Still, recent contributions from deep learning and distance-based predictors show limited success on forecasting bike sharing demand. This consistent observation is hypothesized to be driven by: i) the strong dependence between demand and the meteorological and situational context of stations; and ii) the absence of spatial awareness as most predictors are unable to model the effects of high-low station load on nearby stations. This work proposes a comprehensive set of new principles to incorporate both historical and prospective sources of spatial, meteorological, situational and calendrical context in predictive models of station demand. To this end, a new recurrent neural network layering composed by serial long-short term memory (LSTM) components is proposed with two major contributions: i) the feeding of multivariate time series masks produced from historical context data at the input layer, and ii) the time-dependent regularization of the forecasted time series using prospective context data. This work further assesses the impact of incorporating different sources of context, showing the relevance of the proposed principles for the community even though not all improvements from the context-aware predictors yield statistical significance.


title: Anisotropy-Exchange Resonance as a Mechanism for Entangled State Switching

abstract: We explore the three-particle spin model of an $S_{1}=\frac{1}{2}$ particle (e.g. a stationary electron) interacting with two spin-coupled $S_{\text{2,3}}$ particles with exchange coupling and magnetic anisotropy. We find that in the case of $S_{2,3}=1$ particles, the coupled particle entanglement states can be prepared, controlled, and read by the $S_{1}$ particle. We also find that for particular resonance conditions of the magnetic anisotropy strength $D$ and exchange coupling strength $J$, the entanglement state switching behavior is maximized and is robust against a range of anisotropic application of the exchange coupling.


title: On Wholesale Electricity Prices and Market Values in a Carbon-Neutral Energy System

abstract: Climate and energy policy targets of the European Commission aim to make Europe the first climate-neutral continent by 2050. For low-carbon and net-neutral energy systems primarily based on variable renewable power generation, issues related to the market integration, cannibalisation of revenues, and cost recovery of wind and solar photovoltaics have become major concerns. The traditional discussion of the merit-order effect expects wholesale power prices in a system with 100 % renewable energy sources to alternate between very high and very low values. Unlike previous work, we present a structured and technology-specific analysis of the cross-sectoral demand bidding effect for the price formation in low-carbon power markets. Starting from a stylised market arrangement and by successively augmenting it with all relevant technologies, we construct and quantify the cross-sectoral demand bidding effects in future European power markets with the cross-sectoral market modelling framework SCOPE SD. As the main contribution, we explain and substantiate the market clearing effects of new market participants in detail. Hereby, we put a special focus on hybrid heat supply systems consisting of combined heat and power plant, fuel boiler, thermal storage and electrical back up and derive the opportunity costs of these systems. Furthermore, we show the effects of cross-border integration for a large-scale European net-neutral energy scenario. Finally, the detailed information on market clearing effects allows us to evaluate the resulting revenues of all major technology categories on future electricity markets.


title: Fusing multimodal neuroimaging data with a variational autoencoder

abstract: Neuroimaging studies often involve the collection of multiple data modalities. These modalities contain both shared and mutually exclusive information about the brain. This work aims at finding a scalable and interpretable method to fuse the information of multiple neuroimaging modalities using a variational autoencoder (VAE). To provide an initial assessment, this work evaluates the representations that are learned using a schizophrenia classification task. A support vector machine trained on the representations achieves an area under the curve for the classifier's receiver operating characteristic (ROC-AUC) of 0.8610.


title: Towards A Multi-agent System for Online Hate Speech Detection

abstract: This paper envisions a multi-agent system for detecting the presence of hate speech in online social media platforms such as Twitter and Facebook. We introduce a novel framework employing deep learning techniques to coordinate the channels of textual and im-age processing. Our experimental results aim to demonstrate the effectiveness of our methods for classifying online content, training the proposed neural network model to effectively detect hateful instances in the input. We conclude with a discussion of how our system may be of use to provide recommendations to users who are managing online social networks, showcasing the immense potential of intelligent multi-agent systems towards delivering social good.


title: Homology Conditions for RT Surfaces in Double Holography

abstract: Recently, a novel formula for computing entropy in theories coupled to semi-classical gravity has been devised. Using this so-called island formula the entropy of semi-classical black holes follows a Page curve. Here, we study the relation between this novel entropy and semi-classical entropy in the context of doubly-holographic models. Double holography allows for two different $d$-dimensional descriptions of a black hole coupled to a non-gravitational bath, both of which allow a holographic computation of von Neumann entropy in bath subregions. We argue that the correct homology constraint for Ryu-Takayanagi surfaces depends on which of those $d$-dimensional perspectives is taken. As a consequence the von Neumann entropies of a fixed subregion in both descriptions can disagree. We discuss how the von Neumann entropies in both descriptions are related to the entropy computed by the island formula and coarse grained entropy. Moreover, we argue that the way operators transform between the two descriptions depends on their complexity. A simple toy model is introduced to demonstrate that a sufficiently complicated map between two descriptions of the system can give rise to an island formula and wormholes. Lastly, we speculate about the relation between double-holography and black hole complementarity.


title: Dirac-Born-Infeld warm inflation realization in the strong dissipation regime

abstract: We consider warm inflation with a Dirac-Born-Infeld (DBI) kinetic term in which both the nonequilibrium dissipative particle production and the sound speed parameter slow the motion of the inflaton field. We find that a low sound speed parameter removes, or at least strongly suppresses, the growing function appearing in the scalar of curvature power spectrum of warm inflation, which appears due to the temperature dependence in the dissipation coefficient. As a consequence of that, a low sound speed helps to push warm inflation into the strong dissipation regime, which is an attractive regime from a model building and phenomenological perspective. In turn, the strong dissipation regime of warm inflation softens the microscopic theoretical constraints on cold DBI inflation. The present findings, along with the recent results from swampland criteria, give a strong hint that warm inflation may consistently be embedded into string theory.


title: Systematic study and uncertainty evaluation of $P,T$-odd molecular enhancement factors in BaF

abstract: A measurement of the magnitude of the electric dipole moment of the electron (eEDM) larger than that predicted by the Standard Model (SM) of particle physics is expected to have a huge impact on the search for physics beyond the SM. Polar diatomic molecules containing heavy elements experience enhanced sensitivity to parity ($P$) and time-reversal ($T$)-violating phenomena, such as the eEDM and the scalar-pseudoscalar (S-PS) interaction between the nucleons and the electrons, and are thus promising candidates for measurements. The NL-\textit{e}EDM collaboration is preparing an experiment to measure the eEDM and S-PS interaction in a slow beam of cold BaF molecules [Eur. Phys. J. D, 72, 197 (2018)]. Accurate knowledge of the electronic structure parameters, $W_d$ and $W_s$, connecting the eEDM and the S-PS interaction to the measurable energy shifts is crucial for the interpretation of these measurements. In this work we use the finite field relativistic coupled cluster approach to calculate the $W_d$ and $W_s$ parameters in the ground state of the BaF molecule. Special attention was paid to providing a reliable theoretical uncertainty estimate based on investigations of the basis set, electron correlation, relativistic effects and geometry. Our recommended values of the two parameters, including conservative uncertainty estimates, are 3.13 $\pm$ $0.12 \times 10^{24}\frac{\text{Hz}}{e\cdot \text{cm}}$ for $W_d$ and 8.29 $\pm$ 0.12 kHz for $W_s$.


title: Prediction of clinical tremor severity using Rank Consistent Ordinal Regression

abstract: Tremor is a key diagnostic feature of Parkinson's Disease (PD), Essential Tremor (ET), and other central nervous system (CNS) disorders. Clinicians or trained raters assess tremor severity with TETRAS scores by observing patients. Lacking quantitative measures, inter- or intra- observer variabilities are almost inevitable as the distinction between adjacent tremor scores is subtle. Moreover, clinician assessments also require patient visits, which limits the frequency of disease progress evaluation. Therefore it is beneficial to develop an automated assessment that can be performed remotely and repeatably at patients' convenience for continuous monitoring. In this work, we proposed to train a deep neural network (DNN) with rank-consistent ordinal regression using 276 clinical videos from 36 essential tremor patients. The videos are coupled with clinician assessed TETRAS scores, which are used as ground truth labels to train the DNN. To tackle the challenge of limited training data, optical flows are used to eliminate irrelevant background and statistic objects from RGB frames. In addition to optical flows, transfer learning is also applied to leverage pre-trained network weights from a related task of tremor frequency estimate. The approach was evaluated by splitting the clinical videos into training (67%) and testing sets (0.33%). The mean absolute error on TETRAS score of the testing results is 0.45, indicating that most of the errors were from the mismatch of adjacent labels, which is expected and acceptable. The model predications also agree well with clinical ratings. This model is further applied to smart phone videos collected from a PD patient who has an implanted device to turn "On" or "Off" tremor. The model outputs were consistent with the patient tremor states. The results demonstrate that our trained model can be used as a means to assess and track tremor severity.


title: Quantifying and Maximizing the Benefits of Back-End Noise Adaption on Attention-Based Speech Recognition Models

abstract: This work analyzes how attention-based Bidirectional Long Short-Term Memory (BLSTM) models adapt to noise-augmented speech. We identify crucial components for noise adaptation in BLSTM models by freezing model components during fine-tuning. We first freeze larger model subnetworks and then pursue a fine-grained freezing approach in the encoder after identifying its importance for noise adaptation. The first encoder layer is shown to be crucial for noise adaptation, and the weights are shown to be more important than the other layers. Appreciable accuracy benefits are identified when fine-tuning on a target noisy environment from a model pretrained with noisy speech relative to fine-tuning from a model pretrained with only clean speech when tested on the target noisy environment. For this analysis, we produce our own dataset augmentation tool and it is open-sourced to encourage future efforts in exploring noise adaptation in ASR.


title: Partial Order in the Normal Category Arising from Normal Bands

abstract: The notion of normal category was introduced by KSS Nambooripad in connection with the study of the structure of regular semigroups using cross connections\cite{nambooripad1994theory}. It is an abstraction of the category of principal left ideals of a regular semigroup. A normal band is a semigroup $B$ satisfying $a^2=a$ and $abca=acba$ for all $a,b,c \in B$. Since the normal bands are regular semigroups, the category $\mathcal{L}(B)$ of principal left ideals of a normal band $B$ is a normal category. One of the special properties of this category is that the morphism sets admit a partial order compatible with the composition of morphisms. In this article we derive several properties of this partial order and obtain a new characterization of this partial order.


title: Learning Good State and Action Representations via Tensor Decomposition

abstract: The transition kernel of a continuous-state-action Markov decision process (MDP) admits a natural tensor structure. This paper proposes a tensor-inspired unsupervised learning method to identify meaningful low-dimensional state and action representations from empirical trajectories. The method exploits the MDP's tensor structure by kernelization, importance sampling and low-Tucker-rank approximation. This method can be further used to cluster states and actions respectively and find the best discrete MDP abstraction. We provide sharp statistical error bounds for tensor concentration and the preservation of diffusion distance after embedding. We further prove that the learned state/action abstractions provide accurate approximations to latent block structures if they exist, enabling function approximation in downstream tasks such as policy evaluation.


title: Deterministic and game separability for regular languages of infinite trees

abstract: We show that it is decidable whether two regular languages of infinite trees are separable by a deterministic language, resp., a game language. We consider two variants of separability, depending on whether the set of priorities of the separator is fixed, or not. In each case, we show that separability can be decided in EXPTIME, and that separating automata of exponential size suffice. We obtain our results by reducing to infinite duration games with {\omega}-regular winning conditions and applying the finite-memory determinacy theorem of B\"uchi and Landweber.


title: Fault Tolerant Max-Cut

abstract: In this work, we initiate the study of fault tolerant Max Cut, where given an edge-weighted undirected graph $G=(V,E)$, the goal is to find a cut $S\subseteq V$ that maximizes the total weight of edges that cross $S$ even after an adversary removes $k$ vertices from $G$. We consider two types of adversaries: an adaptive adversary that sees the outcome of the random coin tosses used by the algorithm, and an oblivious adversary that does not. For any constant number of failures $k$ we present an approximation of $(0.878-\epsilon)$ against an adaptive adversary and of $\alpha_{GW}\approx 0.8786$ against an oblivious adversary (here $\alpha_{GW}$ is the approximation achieved by the random hyperplane algorithm of [Goemans-Williamson J. ACM `95]). Additionally, we present a hardness of approximation of $\alpha_{GW}$ against both types of adversaries, rendering our results (virtually) tight. The non-linear nature of the fault tolerant objective makes the design and analysis of algorithms harder when compared to the classic Max Cut. Hence, we employ approaches ranging from multi-objective optimization to LP duality and the ellipsoid algorithm to obtain our results.


title: Fighting off field dependence in MSSM Higgs-mass corrections of order $\alpha_t\,\alpha_s$ and $\alpha_t^2$

abstract: The connection between gauge and Higgs sectors makes supersymmetric extensions of the Standard Model predictive frameworks for the derivation of Higgs masses. In this paper, we study the contamination of such predictions by field-renormalization constants, in the MSSM with two-loop gaugeless corrections of $\mathcal{O}{\left(\alpha_{t,b}\,\alpha_s,\,\alpha_{t,b}^2\right)}$ and full momentum dependence, and demonstrate how strict perturbative expansions allow to systematically neutralize the dependence on such unphysical objects. On the other hand, the popular procedure consisting in an iterative pole search remains explicitly dependent on field counterterms. We then analyze the magnitude of the intrinsic uncertainty that this feature implies for the iterative method, both in non-degenerate and near-degenerate regimes, and conclude that this strategy does not improve on the predictions of the more straightforward expansion. We also discuss several features related to the inclusion of the orders $\alpha_{t,b}\,\alpha_s$ and $\alpha_{t,b}^2$ in the so-called 'fixed-order' approach, such as the resummation of UV-logarithms for heavy supersymmetric spectra.


title: Generalized Adiabatic Theorems: Quantum Systems Driven by Modulated Time-Varying Fields

abstract: We present generalized adiabatic theorems for closed and open quantum systems that can be applied to slow modulations of rapidly varying fields, such as oscillatory fields that occur in optical experiments and light induced processes. The generalized adiabatic theorems show that a sufficiently slow modulation conserves the dynamical modes of time dependent reference Hamiltonians. In the limiting case of modulations of static fields, the standard adiabatic theorems are recovered. Applying these results to periodic fields shows that they remain in Floquet states rather than in energy eigenstates. More generally, these adiabatic theorems can be applied to transformations of arbitrary time-dependent fields, by accounting for the rapidly varying part of the field through the dynamical normal modes, and treating the slow modulation adiabatically. As examples, we apply the generalized theorem to (a) predict the dynamics of a two level system driven by a frequency modulated resonant oscillation, a pathological situation beyond the applicability of earlier results, and (b) to show that open quantum systems driven by slowly turned-on incoherent light, such as biomolecules under natural illumination conditions, can only display coherences that survive in the steady state.


title: Tangent Vector Variational Quantum Eigensolver: A Robust Variational Quantum Eigensolver against the inaccuracy of derivative

abstract: Observing rapid developments of both the number of qubits and quantum volume, especially with recent advances in ion-trap quantum computers, it is no doubt that Fault-Tolerant-Quantum-Computer (FTQC) will be realized in the near future. Since FTQC requires 10,000 physical qubits for every 100 logical ones, it will be used as the first large-scale Noisy-Intermediate-Scale-Quantum (NISQ) . The Variational Quantum Eigensolver (VQE) method will be used until large-scale FTQC with more than 100 logical qubits are realized. Therefore, the VQE method must be improved with respect to both accuracy and time to solution using large resource of the near FTQC . In this paper, we propose Tangent-Vector VQE (TVVQE) method to manage these issues. The method optimizes the norm of tangent vector of trial energy. We demonstrate the calculation of energy levels on Hydrogen molecule, Hubbard model, and Lithium Hydride molecule and reveal that TVVQE has a potential to calculate ground and excited energy levels more accurately than other VQE methods.


title: Relationship among state reopening policies, health outcomes and economic recovery through first wave of the COVID-19 pandemic in the U.S

abstract: State governments in the U.S. have been facing difficult decisions involving tradeoffs between economic and health-related outcomes during the COVID-19 pandemic. Despite evidence of the effectiveness of government-mandated restrictions mitigating the spread of contagion, these orders are stigmatized due to undesirable economic consequences. This tradeoff resulted in state governments employing mandates in widely different ways. We compare the different policies states implemented during periods of restriction (lockdown) and reopening with indicators of COVID-19 spread and consumer card spending at each state during the first wave of the pandemic in the U.S. between March and August 2020. We find that while some states enacted reopening decisions when the incidence rate of COVID-19 was minimal or sustained in its relative decline, other states relaxed socioeconomic restrictions near their highest incidence and prevalence rates experienced so far. Nevertheless, all states experienced similar trends in consumer card spending recovery, which was strongly correlated with reopening policies following the lockdowns and relatively independent from COVID-19 incidence rates at the time. Our findings suggest that consumer card spending patterns can be attributed to government mandates rather than COVID-19 incidence in the states. We estimate the recovery in states that reopened in late April was more than the recovery in states that did not reopen in the same period - 15% for consumer card spending and 18% for spending by high income households. This result highlights the important role of state policies in minimizing health impacts while promoting economic recovery and helps planning effective interventions in subsequent waves and immunization efforts.


title: Traces for factorization homology in dimension 1

abstract: We construct a circle-invariant trace from the factorization homology of the circle $ {\sf trace} \colon \int^\alpha_{{\mathbb S}^1} \\underline{\sf End}(V) \longrightarrow \uno $ associated to a dualizable object $V\in {\boldsymbol{\mathfrak X}}$ in a symmetric monoidal $\infty$-category. This proves a conjecture of To\"en--Vezzosi on existence of circle-invariant traces. Underlying our construction is a calculation of the factorization homology over the circle of the walking adjunction in terms of the paracyclic category of Getzler--Jones: $ \int_{{\mathbb S}^1} {\sf Adj} ~\simeq~ {\bDelta_{\circlearrowleft}}^{\triangleleft\!\triangleright} ~. $ This calculation exhibits a form of Poincar\'e duality for 1-dimensional factorization homology.


title: Euclidean and Hyperbolic Asymmetric Topological Quantum Codes

abstract: In the last three decades, several constructions of quantum error-correcting codes were presented in the literature. Among these codes, there are the asymmetric ones, i.e., quantum codes whose $Z$-distance $d_z$ is different from its $X$-distance $d_x$. The topological quantum codes form an important class of quantum codes, where the toric code, introduced by Kitaev, was the first family of this type. After Kitaev's toric code, several authors focused attention on investigating its structure and the constructions of new families of topological quantum codes over Euclidean and hyperbolic surfaces. As a consequence of establishing the existence and the construction of asymmetric topological quantum codes in Theorem \ref{main}, the main result of this paper, we introduce the class of hyperbolic asymmetric codes. Hence, families of Euclidean and hyperbolic asymmetric topological quantum codes are presented. An analysis regarding the asymptotic behavior of their distances $d_x$ and $d_z$ and encoding rates $k/n$ versus the compact orientable surface's genus is provided due to the significant difference between the asymmetric distances $d_x$ and $d_z$ when compared with the corresponding parameters of topological codes generated by other tessellations. This inherent unequal error-protection is associated with the nontrivial homological cycle of the $\{p,q\}$ tessellation and its dual, which may be appropriately explored depending on the application, where $p\neq q$ and $(p-2)(q-2)\ge 4$. Three families of codes derived from the $\{7,3\}$, $\{5,4\}$, and $\{10,5\}$ tessellations are highlighted.


title: Weak Alfv\'enic turbulence in relativistic plasmas II: Current sheets and dissipation

abstract: Alfv\'{e}n waves as excited in black hole accretion disks and neutron star magnetospheres are the building blocks of turbulence in relativistic, magnetized plasmas. A large reservoir of magnetic energy is available in these systems, such that the plasma can be heated significantly even in the weak turbulence regime. We perform high-resolution three-dimensional simulations of counter-propagating Alfv\'{e}n waves, showing that an $E_{B_{\perp}}(k_{\perp}) \propto k_{\perp}^{-2}$ energy spectrum develops as a result of the weak turbulence cascade in relativistic magnetohydrodynamics and its infinitely magnetized (force-free) limit. The plasma turbulence ubiquitously generates current sheets, which act as locations where magnetic energy dissipates. We show that current sheets form as a natural result of nonlinear interactions between counter-propagating Alfv\'{e}n waves. These current sheets form due to the compression of elongated eddies, driven by the shear induced by growing higher order modes, and undergo a thinning process until they break-up into small-scale turbulent structures. We explore the formation of {current sheets} both in overlapping waves and in localized wave packet collisions. The relativistic interaction of localized Alfv\'{e}n waves induces both Alfv\'{e}n waves and fast waves and efficiently mediates the conversion and dissipation of electromagnetic energy in astrophysical systems. Plasma energization through reconnection in current sheets emerging during the interaction of Alfv\'{e}n waves can potentially explain X-ray emission in black hole accretion coronae and neutron star magnetospheres.

categories: astro-ph.HE astro-ph.SR physics.plasm-ph physics.space-ph


title: Weak Alfv\'{e}nic turbulence in relativistic plasmas. Part 1. Dynamical equations and basic dynamics of interacting resonant triads

abstract: Alfv\'{e}n wave collisions are the primary building blocks of the non-relativistic turbulence that permeates the heliosphere and low-to-moderate energy astrophysical systems. However, many astrophysical systems such as gamma-ray bursts, pulsar and magnetar magnetospheres, and active galactic nuclei have relativistic flows or energy densities. To better understand these high energy systems, we derive reduced relativistic MHD equations and employ them to examine weak Alfv\'{e}nic turbulence, dominated by three-wave interactions, in reduced relativistic magnetohydrodynamics, including the force-free, infinitely magnetized limit. We compare both numerical and analytical solutions to demonstrate that many of the findings from non-relativistic weak turbulence are retained in the relativistic system. But, an important distinction in the relativistic limit is the inapplicability of a formally incompressible limit, i.e, there exists finite coupling to the compressible fast mode regardless of the strength of the magnetic field. Since fast modes can propagate across field lines, this mechanism provides a route for energy to escape strongly magnetized systems, e.g., magnetar magnetospheres. However, we find that the fast-Alfv\'{e}n coupling is diminished in the limit of oblique propagation.


title: Sketches image analysis: Web image search engine usingLSH index and DNN InceptionV3

abstract: The adoption of an appropriate approximate similarity search method is an essential prereq-uisite for developing a fast and efficient CBIR system, especially when dealing with large amount ofdata. In this study we implement a web image search engine on top of a Locality Sensitive Hashing(LSH) Index to allow fast similarity search on deep features. Specifically, we exploit transfer learningfor deep features extraction from images. Firstly, we adopt InceptionV3 pretrained on ImageNet asfeatures extractor, secondly, we try out several CNNs built on top of InceptionV3 as convolutionalbase fine-tuned on our dataset. In both of the previous cases we index the features extracted within ourLSH index implementation so as to compare the retrieval performances with and without fine-tuning.In our approach we try out two different LSH implementations: the first one working with real numberfeature vectors and the second one with the binary transposed version of those vectors. Interestingly,we obtain the best performances when using the binary LSH, reaching almost the same result, in termsof mean average precision, obtained by performing sequential scan of the features, thus avoiding thebias introduced by the LSH index. Lastly, we carry out a performance analysis class by class in terms ofrecall againstmAPhighlighting, as expected, a strong positive correlation between the two.


title: Ultrabroadband nanocavity of hyperbolic phonon polaritons in 1D-like {\alpha}-MoO3

abstract: The exploitation of phonon-polaritons in nanostructured materials offers a pathway to manipulate infrared (IR) light for nanophotonic applications. Notably, hyperbolic phonons polaritons (HP2) in polar bidimensional crystals have been used to demonstrate strong electromagnetic field confinement, ultraslow group velocities, and long lifetimes (~ up to 8 ps). Here we present nanobelts of {\alpha}-phase molybdenum trioxide ({\alpha}-MoO3) as a low-dimensional medium supporting HP2 modes in the mid- and far-IR ranges. By real-space nanoimaging, with IR illuminations provided by synchrotron and tunable lasers, we observe that such HP2 response happens via formation of Fabry-Perot resonances. We remark an anisotropic propagation which critically depends on the frequency range. Our findings are supported by the convergence of experiment, theory, and numerical simulations. Our work shows that the low dimensionality of natural nanostructured crystals, like {\alpha}-MoO3 nanobelts, provides an attractive platform to study polaritonic light-matter interactions and offer appealing cavity properties that could be harnessed in future designs of compact nanophotonic devices.


title: Near-Optimal Cayley Expanders for Abelian Groups

abstract: We give an efficient deterministic algorithm that outputs an expanding generating set for any finite abelian group. The size of the generating set is close to the randomized construction of Alon and Roichman (1994), improving upon various deterministic constructions in both the dependence on the dimension and the spectral gap. By obtaining optimal dependence on the dimension we resolve a conjecture of Azar, Motwani, and Naor (1998) in the affirmative. Our technique is an extension of the bias amplification technique of Ta-Shma (2017), who used random walks on expanders to obtain expanding generating sets over the additive group of n-bit strings. As a consequence, we obtain (i) randomness-efficient constructions of almost k-wise independent variables, (ii) a faster deterministic algorithm for the Remote Point Problem, (iii) randomness-efficient low-degree tests, and (iv) randomness-efficient verification of matrix multiplication.


title: Modeling Social Readers: Novel Tools for Addressing Reception from Online Book Reviews

abstract: Readers' responses to literature have received scant attention in computational literary studies. The rise of social media offers an opportunity to capture a segment of these responses while data-driven analysis of these responses can provide new critical insight into how people "read". Posts discussing an individual book on Goodreads, a social media platform that hosts user discussions of popular literature, are referred to as "reviews", and consist of plot summaries, opinions, quotes, or some mixture of these. Since these reviews are written by readers, computationally modeling them allows one to discover the overall non-professional discussion space about a work, including an aggregated summary of the work's plot, an implicit ranking of the importance of events, and the readers' impressions of main characters. We develop a pipeline of interlocking computational tools to extract a representation of this reader generated shared narrative model. Using a corpus of reviews of five popular novels, we discover the readers' distillation of the main storylines in a novel, their understanding of the relative importance of characters, as well as the readers' varying impressions of these characters. In so doing, we make three important contributions to the study of infinite vocabulary networks: (i) an automatically derived narrative network that includes meta-actants; (ii) a new sequencing algorithm, REV2SEQ, that generates a consensus sequence of events based on partial trajectories aggregated from the reviews; and (iii) a new "impressions" algorithm, SENT2IMP, that provides finer, non-trivial and multi-modal insight into readers' opinions of characters.


title: Pedestrian Detection in 3D Point Clouds using Deep Neural Networks

abstract: Detecting pedestrians is a crucial task in autonomous driving systems to ensure the safety of drivers and pedestrians. The technologies involved in these algorithms must be precise and reliable, regardless of environment conditions. Relying solely on RGB cameras may not be enough to recognize road environments in situations where cameras cannot capture scenes properly. Some approaches aim to compensate for these limitations by combining RGB cameras with TOF sensors, such as LIDARs. However, there are few works that address this problem using exclusively the 3D geometric information provided by LIDARs. In this paper, we propose a PointNet++ based architecture to detect pedestrians in dense 3D point clouds. The aim is to explore the potential contribution of geometric information alone in pedestrian detection systems. We also present a semi-automatic labeling system that transfers pedestrian and non-pedestrian labels from RGB images onto the 3D domain. The fact that our datasets have RGB registered with point clouds enables label transferring by back projection from 2D bounding boxes to point clouds, with only a light manual supervision to validate results. We train PointNet++ with the geometry of the resulting 3D labelled clusters. The evaluation confirms the effectiveness of the proposed method, yielding precision and recall values around 98%.


title: What can the millions of random treatments in nonexperimental data reveal about causes?

abstract: We propose a new method to estimate causal effects from nonexperimental data. Each pair of sample units is first associated with a stochastic 'treatment' - differences in factors between units - and an effect - a resultant outcome difference. It is then proposed that all such pairs can be combined to provide more accurate estimates of causal effects in observational data, provided a statistical model connecting combinatorial properties of treatments to the accuracy and unbiasedness of their effects. The article introduces one such model and a Bayesian approach to combine the $O(n^2)$ pairwise observations typically available in nonexperimnetal data. This also leads to an interpretation of nonexperimental datasets as incomplete, or noisy, versions of ideal factorial experimental designs. This approach to causal effect estimation has several advantages: (1) it expands the number of observations, converting thousands of individuals into millions of observational treatments; (2) starting with treatments closest to the experimental ideal, it identifies noncausal variables that can be ignored in the future, making estimation easier in each subsequent iteration while departing minimally from experiment-like conditions; (3) it recovers individual causal effects in heterogeneous populations. We evaluate the method in simulations and the National Supported Work (NSW) program, an intensively studied program whose effects are known from randomized field experiments. We demonstrate that the proposed approach recovers causal effects in common NSW samples, as well as in arbitrary subpopulations and an order-of-magnitude larger supersample with the entire national program data, outperforming Statistical, Econometrics and Machine Learning estimators in all cases...


title: The Gaia-ESO survey: A lithium depletion boundary age for NGC 2232

abstract: Astrometry and photometry from {\it Gaia} and spectroscopic data from the {\it Gaia}-ESO Survey (GES) are used to identify the lithium depletion boundary (LDB) in the young cluster NGC 2232. A specialised spectral line analysis procedure was used to recover the signature of undepleted lithium in very low luminosity cluster members. An age of $38\pm 3$ Myr is inferred by comparing the LDB location in absolute colour-magnitude diagrams (CMDs) with the predictions of standard models. This is more than twice the age derived from fitting isochrones to low-mass stars in the CMD with the same models. Much closer agreement between LDB and CMD ages is obtained from models that incorporate magnetically suppressed convection or flux-blocking by dark, magnetic starspots. The best agreement is found at ages of $45-50$\,Myr for models with high levels of magnetic activity and starspot coverage fractions $>50$ per cent, although a uniformly high spot coverage does not match the CMD well across the full luminosity range considered.


title: How the 'Auction Cube' Supports the Selection of Auction Designs in Industrial Procurement

abstract: It is well known that rightly applied reverse auctions offer big commercial potential to procurement departments. However, the sheer number of auction types often overwhelms users in practice. And since the implications of a wrongly chosen auction type are equally well known, the overall usage of reverse auctions lacks its potential significantly. In this paper, a novel method is being proposed that guides the user in selecting the right combination of basic auction forms for single lot events, considering both market-, as well as supplier-related, bijective criteria.


title: Trapping Ions and Atoms Optically

abstract: Isolating neutral and charged particles from the environment is essential in precision experiments. For decades, this has been achieved by trapping ions with radio-frequency (rf) fields and neutral particles with optical fields. Recently, trapping of ions by interaction with light has been demonstrated. This might permit combining the advantages of optical trapping and ions. For example, by superimposing optical traps to investigate ensembles of ions and atoms in absence of any radiofrequency fields, as well as to benefit from the versatile and scalable trapping geometries featured by optical lattices. In particular, ions provide individual addressability, electronic and motional degrees of freedom that can be coherently controlled and detected via high fidelity, state-dependent operations. Their long-range Coulomb interaction is significantly larger compared to those of neutral atoms and molecules. This qualifies to study ultra-cold interaction and chemistry of trapped ions and atoms, as well as to provide a novel platform for higher-dimensional experimental quantum simulations. The aim of this topical review is to present the current state of the art and to discuss current challenges and the prospects of the emerging field.


title: Thermomagnetic renormalons in a scalar self-interacting $\lambda \phi^4$ theory

abstract: In this article we extend a previous discussion about the influence of an external magnetic field on renormalons in a self interacting scalar theory by including now temperature effects, in the imaginary formalism, together with an external weak external magnetic field. We show that the location of poles in the Borel plane does not change, getting their residues, however, a dependence on temperature and on the magnetic field.The effects of temperature and the magnetic field strength on the residues turn out to be opposite. We present a detailed discussion about the evolution of these residues, showing technical details involved in the calculation.


title: Efficient Integration of Aggregate Data and Individual Patient Data in One-Way Mixed Models

abstract: Often both Aggregate Data (AD) studies and Individual Patient Data (IPD) studies are available for specific treatments. Combining these two sources of data could improve the overall meta-analytic estimates of treatment effects. Moreover, often for some studies with AD, the associated IPD maybe available, albeit at some extra effort or cost to the analyst. We propose a method for combining treatment effects across trials when the response is from the exponential family of distribution and hence a generalized linear model structure can be used. We consider the case when treatment effects are fixed and common across studies. Using the proposed combination method, we evaluate the wisdom of choosing AD when IPD is available by studying the relative efficiency of analyzing all IPD studies versus combining various percentages of AD and IPD studies. For many different models design constraints under which the AD estimators are the IPD estimators, and hence fully efficient, are known. For such models we advocate a selection procedure that chooses AD studies over IPD studies in a manner that force least departure from design constraints and hence ensures a fully efficient combined AD and IPD estimator.


title: Variance Optimization and Control Regularity for Mean-Field Dynamics

abstract: We study a family of optimal control problems in which one aims at minimizing a cost that mixes a quadratic control penalization and the variance of the system, both for finitely many agents and for the mean-field dynamics as their number goes to infinity. While solutions of the discrete problem always exist in a unique and explicit form, the behavior of their macroscopic counterparts is very sensitive to the magnitude of the time horizon and penalization parameter. When one minimizes the final variance, there always exists a Lipschitz-in-space optimal controls for the infinite dimensional problem, which can be obtained as a suitable extension of the optimal controls for the finite-dimensional problems. The same holds true for variance maximizations whenever the time horizon is sufficiently small. On the contrary, for large final times (or equivalently for small penalizations of the control cost), it can be proven that there does not exist Lipschitz-regular optimal controls for the macroscopic problem.


title: Synthesizing time-series wound prognosis factors from electronic medical records using generative adversarial networks

abstract: Wound prognostic models not only provide an estimate of wound healing time to motivate patients to follow up their treatments but also can help clinicians to decide whether to use a standard care or adjuvant therapies and to assist them with designing clinical trials. However, collecting prognosis factors from Electronic Medical Records (EMR) of patients is challenging due to privacy, sensitivity, and confidentiality. In this study, we developed time series medical generative adversarial networks (GANs) to generate synthetic wound prognosis factors using very limited information collected during routine care in a specialized wound care facility. The generated prognosis variables are used in developing a predictive model for chronic wound healing trajectory. Our novel medical GAN can produce both continuous and categorical features from EMR. Moreover, we applied temporal information to our model by considering data collected from the weekly follow-ups of patients. Conditional training strategies were utilized to enhance training and generate classified data in terms of healing or non-healing. The ability of the proposed model to generate realistic EMR data was evaluated by TSTR (test on the synthetic, train on the real), discriminative accuracy, and visualization. We utilized samples generated by our proposed GAN in training a prognosis model to demonstrate its real-life application. Using the generated samples in training predictive models improved the classification accuracy by 6.66-10.01% compared to the previous EMR-GAN. Additionally, the suggested prognosis classifier has achieved the area under the curve (AUC) of 0.975, 0.968, and 0.849 when training the network using data from the first three visits, first two visits, and first visit, respectively. These results indicate a significant improvement in wound healing prediction compared to the previous prognosis models.


title: The Tracking Machine Learning challenge : Throughput phase

abstract: This paper reports on the second "Throughput" phase of the Tracking Machine Learning (TrackML) challenge on the Codalab platform. As in the first "Accuracy" phase, the participants had to solve a difficult experimental problem linked to tracking accurately the trajectory of particles as e.g. created at the Large Hadron Collider (LHC): given O($10^5$) points, the participants had to connect them into O($10^4$) individual groups that represent the particle trajectories which are approximated helical. While in the first phase only the accuracy mattered, the goal of this second phase was a compromise between the accuracy and the speed of inference. Both were measured on the Codalab platform where the participants had to upload their software. The best three participants had solutions with good accuracy and speed an order of magnitude faster than the state of the art when the challenge was designed. Although the core algorithms were less diverse than in the first phase, a diversity of techniques have been used and are described in this paper. The performance of the algorithms are analysed in depth and lessons derived.


title: Sketching approximability of all finite CSPs

abstract: A constraint satisfaction problem (CSP), $\textsf{Max-CSP}(\mathcal{F})$, is specified by a finite set of constraints $\mathcal{F} \subseteq \{[q]^k \to \{0,1\}\}$ for positive integers $q$ and $k$. An instance of the problem on $n$ variables is given by $m$ applications of constraints from $\mathcal{F}$ to subsequences of the $n$ variables, and the goal is to find an assignment to the variables that satisfies the maximum number of constraints. In the $(\gamma,\beta)$-approximation version of the problem for parameters $0 \leq \beta < \gamma \leq 1$, the goal is to distinguish instances where at least $\gamma$ fraction of the constraints can be satisfied from instances where at most $\beta$ fraction of the constraints can be satisfied. In this work we consider the approximability of this problem in the context of sketching algorithms and give a dichotomy result. Specifically, for every family $\mathcal{F}$ and every $\beta < \gamma$, we show that either a linear sketching algorithm solves the problem in polylogarithmic space, or the problem is not solvable by any sketching algorithm in $o(\sqrt{n})$ space. In particular, we give non-trivial approximation algorithms using polylogarithmic space for infinitely many constraint satisfaction problems. We also extend previously known lower bounds for general streaming algorithms to a wide variety of problems, and in particular the case of $q=k=2$, where we get a dichotomy, and the case when the satisfying assignments of the constraints of $\mathcal{F}$ support a distribution on $[q]^k$ with uniform marginals. Prior to this work, other than sporadic examples, the only systematic classes of CSPs that were analyzed considered the setting of Boolean variables $q=2$, binary constraints $k=2$, singleton families $|\mathcal{F}|=1$ and only considered the setting where constraints are placed on literals rather than variables.


title: Polarization signatures of a high-velocity scatterer in nebular-phase spectra of Type II supernovae

abstract: Type II supernovae (SNe) often exhibit a linear polarization, arising from free-electron scattering, with complicated optical signatures, both in the continuum and in lines. Focusing on the early nebular phase, at a SN age of 200d, we conduct a systematic study of the polarization signatures associated with a 56Ni `blob' that breaks spherical symmetry. Our ansatz, supported by nonLTE radiative transfer calculations, is that the primary role of such a 56Ni blob is to boost the local density of free electrons, which is otherwise reduced following recombination in SNe II. Using 2D polarized radiation transfer modeling, we explore the influence of such an electron-density enhancement, varying its magnitude N_e_fac, its velocity location V_blob, and its spatial extent. For plausible N_e_fac values of a few tens, a high-velocity blob can deliver a continuum polarization P_cont of 0.5-1.0% at 200d. Our simulations reproduce the analytic scalings for P_cont, and in particular the linear growth with the blob radial optical depth. The most constraining information is, however, carried by polarized line photons. For a high V_blob, the polarized spectrum appears as a replica of the full spectrum, scaled down by a factor 100 to 1000 (i.e., 1/P_cont), and redshifted by an amount V_blob(1-cos(alpha_los)), where alpha_los is the line of sight angle. As V_blob is reduced, the redshift decreases and the replication deteriorates. Lines whose formation region overlap with the blob appear weaker and narrower in the polarized flux. Because of its dependence on inclination (~ sin^2 alpha_los), the polarization preferentially reveals asymmetries in the plane perpendicular to the line of sight. With the adequate choice of electron-density enhancement, some of these results may apply to asymmetric explosions in general, or to the polarization signatures from newly-formed dust in the outer ejecta. [abridged]


title: High-order space-time finite element methods for the Poisson-Nernst-Planck equations: Positivity and unconditional energy stability

abstract: We present a novel class of high-order space-time finite element schemes for the Poisson-Nernst-Planck (PNP) equations. We prove that our schemes are mass conservative, positivity preserving, and unconditionally energy stable for any order of approximation. To the best of our knowledge, this is the first class of (arbitrarily) high-order accurate schemes for the PNP equations that simultaneously achieve all these three properties. This is accomplished via (1) using finite elements to directly approximate the so-called entropy variable instead of the density variable, and (2) using a discontinuous Galerkin (DG) discretization in time. The entropy variable formulation, which was originally developed by Metti et al. [17] under the name of a log-density formulation, guarantees both positivity of densities and a continuous-in-time energy stability result. The DG in time discretization further ensures an unconditional energy stability in the fully discrete level for any approximation order, where the lowest order case is exactly the backward Euler discretization and in this case we recover the method of Metti et al. [17].


title: The science case for spacecraft exploration of the Uranian satellites: Candidate ocean worlds in an ice giant system

abstract: The 27 satellites of Uranus are enigmatic, with dark surfaces coated by material that could be rich in organics. Voyager 2 imaged the southern hemispheres of Uranus' five largest 'classical' moons Miranda, Ariel, Umbriel, Titania, and Oberon, as well as the largest ring moon Puck, but their northern hemispheres were largely unobservable at the time of the flyby and were not imaged. Additionally, no spatially resolved datasets exist for the other 21 known moons, and their surface properties are essentially unknown. Because Voyager 2 was not equipped with a near-infrared mapping spectrometer, our knowledge of the Uranian moons' surface compositions, and the processes that modify them, is limited to disk-integrated datasets collected by ground- and space-based telescopes. Nevertheless, images collected by the Imaging Science System on Voyager 2 and reflectance spectra collected by telescope facilities indicate that the five classical moons are candidate ocean worlds that might currently have, or had, liquid subsurface layers beneath their icy surfaces. To determine whether these moons are ocean worlds, and investigate Uranus' ring moons and irregular satellites, close-up observations and measurements made by instruments onboard a Uranus orbiter are needed.


title: Explicit formulas for the inverses of Toeplitz matrices, with applications

abstract: We derive novel explicit formulas for the inverses of truncated block Toeplitz matrices that correspond to a multivariate minimal stationary process. The main ingredients of the formulas are the Fourier coefficients of the phase function attached to the spectral density of the process. The derivation of the formulas is based on a recently developed finite prediction theory applied to the dual process of the stationary process. We illustrate the usefulness of the formulas by two applications. The first one is a strong convergence result for solutions of general block Toeplitz systems for a multivariate short-memory process. The second application is closed-form formulas for the inverses of truncated block Toeplitz matrices corresponding to a multivariate ARMA process. The significance of the latter is that they provide us with a linear-time algorithm to compute the solutions of corresponding block Toeplitz systems.


title: Global 21cm Absorption Signal from Superconducting Cosmic Strings

abstract: Superconducting cosmic strings emit electromagnetic waves between the times of recombination and reionization. Hence, they have an effect on the global 21cm signal. We compute the resulting absorption features, focusing on strings with critical current, study their dependence on the string tension $\mu$, and compare with observational results. For string tensions in the range of $G \mu = 10^{-10}$, where $G$ is Newton's gravitational constant, there is an interesting amplification of the two characteristic absorption features, one during the cosmic dawn, $z \lesssim 30$, and the other during the cosmic dark age, $z \sim 80$, the former being comparable in amplitude to what was observed by the EDGES experiment.


title: Evolution of cooperation with asymmetric social interactions

abstract: How cooperation emerges in human societies is both an evolutionary enigma, and a practical problem with tangible implications for societal health. Population structure has long been recognized as a catalyst for cooperation because local interactions enable reciprocity. Analysis of this phenomenon typically assumes bi-directional social interactions, even though real-world interactions are often uni-directional. Uni-directional interactions -- where one individual has the opportunity to contribute altruistically to another, but not conversely -- arise in real-world populations as the result of organizational hierarchies, social stratification, popularity effects, and endogenous mechanisms of network growth. Here we expand the theory of cooperation in structured populations to account for both uni- and bi-directional social interactions. Even though directed interactions remove the opportunity for reciprocity, we find that cooperation can nonetheless be favored in directed social networks and that cooperation is provably maximized for networks with an intermediate proportion of directed interactions, as observed in many empirical settings. We also identify two simple structural motifs that allow efficient modification of interaction directionality to promote cooperation by orders of magnitude. We discuss how our results relate to the concepts of generalized and indirect reciprocity.


title: Neural network surrogate of QuaLiKiz using JET experimental data to populate training space

abstract: Within integrated tokamak plasma modelling, turbulent transport codes are typically the computational bottleneck limiting their routine use outside of post-discharge analysis. Neural network (NN) surrogates have been used to accelerate these calculations while retaining the desired accuracy of the physics-based models. This paper extends a previous NN model, known as QLKNN-hyper-10D, by incorporating the impact of impurities, plasma rotation and magnetic equilibrium effects. This is achieved by adding a light impurity fractional density ($n_{imp,light} / n_e$) and its normalized gradient, the normalized pressure gradient ($\alpha$), the toroidal Mach number ($M_{tor}$) and the normalized toroidal flow velocity gradient. The input space was sampled based on experimental data from the JET tokamak to avoid the curse of dimensionality. The resulting networks, named QLKNN-jetexp-15D, show good agreement with the original QuaLiKiz model, both by comparing individual transport quantity predictions as well as comparing its impact within the integrated model, JINTRAC. The profile-averaged RMS of the integrated modelling simulations is <10% for each of the 5 scenarios tested. This is non-trivial given the potential numerical instabilities present within the highly nonlinear system of equations governing plasma transport, especially considering the novel addition of momentum flux predictions to the model proposed here. An evaluation of all 25 NN output quantities at one radial location takes $\sim$0.1 ms, $10^4$ times faster than the original QuaLiKiz model. Within the JINTRAC integrated modelling tests performed in this study, using QLKNN-jetexp-15D resulted in a speed increase of only 60 - 100 as other physics modules outside of turbulent transport become the bottleneck.


title: First M87 Event Horizon Telescope Results. VII. Polarization of the Ring

abstract: In 2017 April, the Event Horizon Telescope (EHT) observed the near-horizon region around the supermassive black hole at the core of the M87 galaxy. These 1.3 mm wavelength observations revealed a compact asymmetric ring-like source morphology. This structure originates from synchrotron emission produced by relativistic plasma located in the immediate vicinity of the black hole. Here we present the corresponding linear-polarimetric EHT images of the center of M87. We find that only a part of the ring is significantly polarized. The resolved fractional linear polarization has a maximum located in the southwest part of the ring, where it rises to the level of about 15%. The polarization position angles are arranged in a nearly azimuthal pattern. We perform quantitative measurements of relevant polarimetric properties of the compact emission and find evidence for the temporal evolution of the polarized source structure over one week of EHT observations. The details of the polarimetric data reduction and calibration methodology are provided. We carry out the data analysis using multiple independent imaging and modeling techniques, each of which is validated against a suite of synthetic data sets. The gross polarimetric structure and its apparent evolution with time are insensitive to the method used to reconstruct the image. These polarimetric images carry information about the structure of the magnetic fields responsible for the synchrotron emission. Their physical interpretation is discussed in an accompanying publication.


title: Structured Matrix Approximations via Tensor Decompositions

abstract: We provide a computational framework for approximating a class of structured matrices; here, the term structure is very general, and may refer to a regular sparsity pattern (e.g., block-banded), or be more highly structured (e.g., symmetric block Toeplitz). The goal is to uncover {\it additional latent structure} that will in turn lead to computationally efficient algorithms when the new structured matrix approximations are employed in the place of the original operator. Our approach has three steps: map the structured matrix to tensors, use tensor compression algorithms, and map the compressed tensors back to obtain two different matrix representations -- sum of Kronecker products and block low-rank format. The use of tensor decompositions enables us to uncover latent structure in the problem and leads to compressed representations of the original matrix that can be used efficiently in applications. The resulting matrix approximations are memory efficient, easy to compute with, and preserve the error that is due to the tensor compression in the Frobenius norm. Our framework is quite general. We illustrate the ability of our method to uncover block-low-rank format on structured matrices from two applications: system identification, space-time covariance matrices. In addition, we demonstrate that our approach can uncover sum of structured Kronecker products structure on several matrices from the SuiteSparse collection. Finally, we show that our framework is broad enough to encompass and improve on other related results from the literature, as we illustrate with the approximation of a three-dimensional blurring operator.


title: Machine Learning Applications for Therapeutic Tasks with Genomics Data

abstract: Thanks to the increasing availability of genomics and other biomedical data, many machine learning approaches have been proposed for a wide range of therapeutic discovery and development tasks. In this survey, we review the literature on machine learning applications for genomics through the lens of therapeutic development. We investigate the interplay among genomics, compounds, proteins, electronic health records (EHR), cellular images, and clinical texts. We identify twenty-two machine learning in genomics applications across the entire therapeutics pipeline, from discovering novel targets, personalized medicine, developing gene-editing tools all the way to clinical trials and post-market studies. We also pinpoint seven important challenges in this field with opportunities for expansion and impact. This survey overviews recent research at the intersection of machine learning, genomics, and therapeutic development.


title: Crystals of gauged solitons, force free plasma and resurgence

abstract: We show that the (3+1)-dimensional gauged non-linear sigma model minimally coupled to a U(1) gauge field possesses analytic solutions representing gauged solitons at finite Baryon density whose electromagnetic field is a Force Free Plasma. These gauged solitons manifest a crystalline structure and generate in a very natural way persistent currents able to support Force Free Plasma electromagnetic fields. The trajectories of charged test particles moving within these configurations can be characterized. Quite surprisingly, despite the non-integrable nature of the theory, some of the perturbations of these gauged solitons allow to identify a proper resurgent parameter. In particular, the perturbations of the solitons profile are related to the Lam\'e operator. On the other hand, the electromagnetic perturbations on the configurations satisfy a two-dimensional effective Schrodinger equation, where the soliton background interacts with the electromagnetic perturbations through an effective two-dimensional periodic potential. We studied numerically the band energy spectrum for different values of the free parameters of the theory and we found that bands-gaps are modulated by the potential strength. Finally we compare our crystal solutions with those of the (1+1)- dimesional Gross-Neveu model.


title: First M87 Event Horizon Telescope Results. VIII. Magnetic Field Structure near The Event Horizon

abstract: Event Horizon Telescope (EHT) observations at 230 GHz have now imaged polarized emission around the supermassive black hole in M87 on event-horizon scales. This polarized synchrotron radiation probes the structure of magnetic fields and the plasma properties near the black hole. Here we compare the resolved polarization structure observed by the EHT, along with simultaneous unresolved observations with the Atacama Large Millimeter/submillimeter Array, to expectations from theoretical models. The low fractional linear polarization in the resolved image suggests that the polarization is scrambled on scales smaller than the EHT beam, which we attribute to Faraday rotation internal to the emission region. We estimate the average density n_e of order 10^4-7 cm-3, magnetic field strength B of order 1-30 G, and electron temperature Te of order (1-12) x 10^10 K of the radiating plasma in a simple one-zone emission model. We show that the net azimuthal linear polarization pattern may result from organized, poloidal magnetic fields in the emission region. In a quantitative comparison with a large library of simulated polarimetric images from general relativistic magnetohydrodynamic (GRMHD) simulations, we identify a subset of physical models that can explain critical features of the polarimetric EHT observations while producing a relativistic jet of sufficient power. The consistent GRMHD models are all of magnetically arrested accretion disks, where near-horizon magnetic fields are dynamically important. We use the models to infer a mass accretion rate onto the black hole in M87 of (3-20) x 10^-4 Msun yr-1.


title: On deformation spaces of toric singularities and on singularities of K-moduli of Fano varieties

abstract: Firstly, we see that the bases of the miniversal deformations of isolated $\mathbb{Q}$-Gorenstein toric singularities are quite restricted. In particular, we classify the analytic germs of embedding dimension $\leq 2$ which are the bases of the miniversal deformations of isolated $\mathbb{Q}$-Gorenstein toric singularities. Secondly, we show that the deformation spaces of isolated Gorenstein toric $3$-fold singularities appear, in a weak sense, as singularities of the K-moduli stack of K-semistable Fano varieties of every dimension $\geq 3$. As a consequence, we prove that the number of local branches of the K-moduli stack of K-semistable Fano varieties and of the K-moduli space of K-polystable Fano varieties is unbounded in each dimension $\geq 3$.


title: Towards Maximum Optical Efficiency of Ensembles of Colloidal Nanorods

abstract: Experimental and theoretical studies of colloidal nanoparticles have primarily focused on accurate characterization and simulation of observable characteristics, such as resonant wavelength. In this Letter, we tackle the optimal design of colloidal-nanoparticle ensembles: what is the largest possible optical response, which designs might achieve them, and can such response be experimentally demonstrated? We combine theory and experiment to answer each of these questions. We derive general bounds on the maximum cross-sections per volume, and we apply an analytical antenna model to show that resonant nanorods should nearly achieve such bounds. We use a modified seed-mediated synthesis approach to synthesize ensembles of gold nanorods with small polydispersity, i.e., small variations in size and aspect ratio. Polydispersity is the key determinant of how closely such ensembles can approach their respective bounds yet is difficult to characterize experimentally without near-field measurements. We show that a certain "extinction metric," connecting extinction cross-section per volume with the radiative efficiencies of the nanoparticles, offers a quantitative prediction of polydispersity via quantities that can be rapidly measured with far-field characterization tools. Our predictions apply generally across all plasmonic materials and offers a roadmap to the largest possible optical response of nanoparticle ensembles.


title: Moir\'e tuning of spin excitations: Individual Fe atoms on MoS$_2$/Au(111)

abstract: Magnetic adatoms on surfaces are exchange coupled to their environment, inducing fast relaxation of excited spin states and decoherence. These interactions can be suppressed by inserting decoupling layers between the magnetic adsorbate and the metallic substrate. Using low-temperature scanning tunneling microscopy and spectroscopy, we investigate the magnetic properties of single Fe atoms on Au(111) covered by a monolayer of MoS$_2$ and explore the influence of the moir\'e structure on the spin excitation spectra. Fe adatoms absorbed at minima of the moir\'e structure exhibit almost pure inelastic spin excitations. Kondo correlations develop for adsorption sites off the moir\'e minimum, culminating in fully developed Kondo resonances on moir\'e maxima. The increase in Kondo correlations is accompanied by a decrease in the single-ion anisotropy, which we rationalize by poor-person's scaling. Local variations in the tunneling spectra indicate pronounced interference between tunneling paths involving the spin-carrying orbitals and a hybrid Fe-S orbital, which forms despite the nominally inert nature of the terminating S layer.


title: Application of Gaussian process regression to plasma turbulent transport model validation via integrated modelling

abstract: This paper outlines an approach towards improved rigour in tokamak turbulence transport model validation within integrated modelling. Gaussian process regression (GPR) techniques were applied for profile fitting during the preparation of integrated modelling simulations allowing for rigourous sensitivity tests of prescribed initial and boundary conditions as both fit and derivative uncertainties are provided. This was demonstrated by a JETTO integrated modelling simulation of the JET ITER-like-wall H-mode baseline discharge #92436 with the QuaLiKiz quasilinear turbulent transport model, which is the subject of extrapolation towards a deuterium-tritium plasma. The simulation simultaneously evaluates the time evolution of heat, particle, and momentum fluxes over $\sim10$ confinement times, with a simulation boundary condition at $\rho_{tor} = 0.85$. Routine inclusion of momentum transport prediction in multi-channel flux-driven transport modelling is not standard and is facilitated here by recent developments within the QuaLiKiz model. Excellent agreement was achieved between the fitted and simulated profiles for $n_e$, $T_e$, $T_i$, and $\Omega_{tor}$ within $2\sigma$, but the simulation underpredicts the mid-radius $T_i$ and overpredicts the core $n_e$ and $T_e$ profiles for this discharge. Despite this, it was shown that this approach is capable of deriving reasonable inputs, including derivative quantities, to tokamak models from experimental data. Furthermore, multiple figures-of-merit were defined to quantitatively assess the agreement of integrated modelling predictions to experimental data within the GPR profile fitting framework.


title: Single eigenvalue fluctuations of general Wigner-type matrices

abstract: We consider the single eigenvalue fluctuations of random matrices of general Wigner-type, under a one-cut assumption on the density of states. For eigenvalues in the bulk, we prove that the asymptotic fluctuations of a single eigenvalue around its classical location are Gaussian with a universal variance. Our method is based on a dynamical approach to mesoscopic linear spectral statistics which reduces their behavior on short scales to that on larger scales. We prove a central limit theorem for linear spectral statistics on larger scales via resolvent techniques and show that for certain classes of test functions, the leading-order contribution to the variance agrees with the GOE/GUE cases.


title: Degrees of Restriction for Two-Dimensional Automata

abstract: A three-way (resp., two-way) two-dimensional automaton has a read-only input head that moves in three (resp., two) directions on a finite array of cells labelled by symbols of the input alphabet. Restricting the input head movement of a two-dimensional automaton results in a model that is weaker in terms of recognition power. In this paper, we introduce the notion of "degrees of restriction" for two-dimensional automata, and we develop sets of extended two-dimensional automaton models that allow for some bounded number of restricted moves. We establish recognition hierarchies for both deterministic and nondeterministic extended three-way two-dimensional automata, and we find similar hierarchies for both deterministic and nondeterministic extended two-way two-dimensional automata. We also prove incomparability results between nondeterministic and deterministic extended three-way two-dimensional automata. Lastly, we consider closure properties for some operations on languages recognized by extended three-way two-dimensional automata.


title: Scalar Adjective Identification and Multilingual Ranking

abstract: The intensity relationship that holds between scalar adjectives (e.g., nice < great < wonderful) is highly relevant for natural language inference and common-sense reasoning. Previous research on scalar adjective ranking has focused on English, mainly due to the availability of datasets for evaluation. We introduce a new multilingual dataset in order to promote research on scalar adjectives in new languages. We perform a series of experiments and set performance baselines on this dataset, using monolingual and multilingual contextual language models. Additionally, we introduce a new binary classification task for English scalar adjective identification which examines the models' ability to distinguish scalar from relational adjectives. We probe contextualised representations and report baseline results for future comparison on this task.


title: Automated Estimation of Total Lung Volume using Chest Radiographs and Deep Learning

abstract: Total lung volume is an important quantitative biomarker and is used for the assessment of restrictive lung diseases. In this study, we investigate the performance of several deep-learning approaches for automated measurement of total lung volume from chest radiographs. 7621 posteroanterior and lateral view chest radiographs (CXR) were collected from patients with chest CT available. Similarly, 928 CXR studies were chosen from patients with pulmonary function test (PFT) results. The reference total lung volume was calculated from lung segmentation on CT or PFT data, respectively. This dataset was used to train deep-learning architectures to predict total lung volume from chest radiographs. The experiments were constructed in a step-wise fashion with increasing complexity to demonstrate the effect of training with CT-derived labels only and the sources of error. The optimal models were tested on 291 CXR studies with reference lung volume obtained from PFT. The optimal deep-learning regression model showed an MAE of 408 ml and a MAPE of 8.1\% and Pearson's r = 0.92 using both frontal and lateral chest radiographs as input. CT-derived labels were useful for pre-training but the optimal performance was obtained by fine-tuning the network with PFT-derived labels. We demonstrate, for the first time, that state-of-the-art deep learning solutions can accurately measure total lung volume from plain chest radiographs. The proposed model can be used to obtain total lung volume from routinely acquired chest radiographs at no additional cost and could be a useful tool to identify trends over time in patients referred regularly for chest x-rays.


title: Quadrangular ${\mathbb Z}_{p}^{l}$-actions on Riemann surfaces

abstract: Let $p \geq 3$ be a prime integer and, for $l \geq 1$, let $G \cong {\mathbb Z}_{p}^{l}$ be a group of conformal automorphisms of some closed Riemann surface $S$ of genus $g \geq 2$. By the Riemann-Hurwitz formula, either $p \leq g+1$ or $p=2g+1$. If $l=1$ and $p=2g+1$, then $S/G$ is the sphere with exactly three cone points and, if moreover $p \geq 7$, then $G$ is the unique $p$-Sylow subgroup of ${\rm Aut}(S)$. If $l=1$ and $p=g+1$, then $S/G$ is the sphere with exactly four cone points and, if moreover $p \geq 13$, then $G$ is again the unique $p$-Sylow subgroup. The above unique facts permited many authors to obtain algebraic models and the corresponding groups ${\rm Aut}(S)$ in these situations. Now, let us assume $l \geq 2$. If $p \geq 5$, then either (i) $p^{l} \leq g-1$ or (ii) $S/G$ has genus zero, $p^{l-1}(p-3) \leq 2(g-1)$ and $2 \leq l \leq r-1$, where $r \geq 3$ is the number of cone points of $S/G$. Let us assume we are in case (ii). If $r=3$, then $l=2$ and $S$ happens to be the classical Fermat curve of degree $p$, whose group of automorphisms is well known. The next case, $r=4$, is studied in this paper. We provide an algebraic curve representation for $S$, a description of its group of conformal automorphisms, a discussion of its field of moduli and an isogenous decomposition of its jacobian variety.


title: Octopus: A Zero-Cost Architecture for Stream Network Monitoring

abstract: Considering the growing demand and popularity of Do-It-Yourself (DIY) networks, low-cost devices managed by the people, for the people and the ease of deployment of localised/decentralised Internet services, it is convenient for such networks to have an efficient low-cost monitoring platform to ensure availability, responsiveness and users' Quality of Experience (QoE). This is especially relevant for the developing world where Community Networks (CN) are increasing in popularity and complexity. In this letter, we discuss Octopus: an architecture for stream network monitoring.


title: Reconciling design-based and model-based causal inferences for split-plot experiments

abstract: The split-plot design assigns different interventions at the whole-plot and sub-plot levels, respectively, and induces a group structure on the final treatment assignments. A common strategy is to use the OLS fit of the outcome on the treatment indicators coupled with the robust standard errors clustered at the whole-plot level. It does not give consistent estimator for the causal effects of interest when the whole-plot sizes vary. Another common strategy is to fit the linear mixed-effects model of the outcome with Normal random effects and errors. It is a purely model-based approach and can be sensitive to violations of parametric assumptions. In contrast, the design-based inference assumes no outcome models and relies solely on the controllable randomization mechanism determined by the physical experiment. We first extend the existing design-based inference based on the {\htf} estimator to the Hajek estimator, and establish the finite-population central limit theorem for both under split-plot randomization. We then reconcile the results with those under the model-based approach, and propose two regression strategies, namely (i) the WLS fit of the unit-level data based on the inverse probability weighting and (ii) the OLS fit of the aggregate data based on whole-plot total outcomes, to reproduce the Hajek and {\htf} estimators from least squares, respectively. This, together with the asymptotic conservativeness of the corresponding cluster-robust covariances for estimating the true design-based covariances as we establish in the process, justifies the validity of regression-based estimators for design-based inference. In light of the flexibility of regression formulation with covariate adjustment, we further extend the theory to the case with covariates and demonstrate the efficiency gain by regression-based covariate adjustment via both asymptotic theory and simulation.


title: Effective temperatures for single particle system under dichotomous noise

abstract: Three different definitions of effective temperature -- $\mathcal{T}_{\rm k}$, $\mathcal{T}_{\rm i}$ and $\mathcal{T}_{\rm r}$ related to kinetic theory, system entropy and response theory, respectively -- are applied in the description of a non-equilibrium generalised massive Langevin model in contact with dichotomous noise. The differences between the definitions of $\mathcal{T}$ naturally wade out as the reservoir reaches its white-noise limit, approaching Gaussian features. The same framework is employed in its overdamped version as well, showing the loss of inertial contributions to the dynamics of the system also makes the three mentioned approaches for effective temperature equivalent.


title: Approach to nonlinear magnetohydrodynamic simulations in stellarator geometry

abstract: The capability to model the nonlinear magnetohydrodynamic (MHD) evolution of stellarator plasmas is developed by extending the M3D-$C^1$ code to allow non-axisymmetric domain geometry. We introduce a set of logical coordinates, in which the computational domain is axisymmetric, to utilize the existing finite-element framework of M3D-$C^1$. A $C^1$ coordinate mapping connects the logical domain to the non-axisymmetric physical domain, where we use the M3D-$C^1$ extended MHD models essentially without modifications. We present several numerical verifications on the implementation of this approach, including simulations of the heating, destabilization, and equilibration of stellarator plasmas with strongly anisotropic thermal conductivity, and of the relaxation of stellarator equilibria to integrable and non-integrable magnetic field configurations in realistic geometries.


title: Proximal Learning for Individualized Treatment Regimes Under Unmeasured Confounding

abstract: Data-driven individualized decision making has recently received increasing research interests. Most existing methods rely on the assumption of no unmeasured confounding, which unfortunately cannot be ensured in practice especially in observational studies. Motivated by the recent proposed proximal causal inference, we develop several proximal learning approaches to estimating optimal individualized treatment regimes (ITRs) in the presence of unmeasured confounding. In particular, we establish several identification results for different classes of ITRs, exhibiting the trade-off between the risk of making untestable assumptions and the value function improvement in decision making. Based on these results, we propose several classification-based approaches to finding a variety of restricted in-class optimal ITRs and develop their theoretical properties. The appealing numerical performance of our proposed methods is demonstrated via an extensive simulation study and one real data application.


title: A Finite-temperature Phase Transition for the Ising Spin-glass in $d\geq 2$

abstract: It is believed that the $\pm J$ Ising spin-glass does not order at finite temperatures in dimension $d=2$. However, using a graphical representation and a contour argument, we prove rigorously the existence of a finite-temperature phase transition in $d\geq 2$ with $T_c \geq 0.4$. In the graphical representation, the low-temperature phase allows for the coexistence of multiple infinite clusters each with a rigidly aligned spin-overlap state. These clusters correlate negatively with each other, and are entropically stable without breaking any global symmetry. They can emerge in most graph structures and disorder measures.

categories: cond-mat.dis-nn cond-mat.stat-mech math-ph math.MP


title: Textual Analysis of Communications in COVID-19 Infected Community on Social Media

abstract: During the COVID-19 pandemic, people started to discuss about pandemic-related topics on social media. On subreddit \textit{r/COVID19positive}, a number of topics are discussed or being shared, including experience of those who got a positive test result, stories of those who presumably got infected, and questions asked regarding the pandemic and the disease. In this study, we try to understand, from a linguistic perspective, the nature of discussions on the subreddit. We found differences in linguistic characteristics (e.g. psychological, emotional and reasoning) across three different categories of topics. We also classified posts into the different categories using SOTA pre-trained language models. Such classification model can be used for pandemic-related research on social media.


title: All-optical fluorescence blinking control in quantum dots with ultrafast mid-infrared pulses

abstract: Photoluminescence (PL) intermittency is a ubiquitous phenomenon detrimentally reducing the temporal emission intensity stability of single colloidal quantum dots (CQDs) and the emission quantum yield of their ensembles. Despite efforts for blinking reduction via chemical engineering of the QD architecture and its environment, blinking still poses barriers to the application of QDs, particularly in single-particle tracking in biology or in single-photon sources. Here, we demonstrate the first deterministic all-optical suppression of quantum dot blinking using a compound technique of visible and mid-infrared (MIR) excitation. We show that moderate-field ultrafast MIR pulses (5.5 $\mu$m, 150 fs) can switch the emission from a charged, low quantum yield 'grey' trion state to the 'bright' exciton state in CdSe/CdS core-shell quantum dots resulting in a significant reduction of the QD intensity flicker. Quantum-tunneling simulations suggest that the MIR fields remove the excess charge from trions with reduced emission quantum yield to restore higher brightness exciton emission. Our approach can be integrated with existing single-particle tracking or super-resolution microscopy techniques without any modification to the sample and translates to other emitters presenting charging-induced PL intermittencies, such as single-photon emissive defects in diamond and two-dimensional materials.

categories: cond-mat.mes-hall cond-mat.mtrl-sci physics.chem-ph physics.optics


title: Makespan Trade-offs for Visiting Triangle Edges

abstract: We study a primitive vehicle routing-type problem in which a fleet of $n$unit speed robots start from a point within a non-obtuse triangle $\Delta$, where $n \in \{1,2,3\}$. The goal is to design robots' trajectories so as to visit all edges of the triangle with the smallest visitation time makespan. We begin our study by introducing a framework for subdividing $\Delta$into regions with respect to the type of optimal trajectory that each point $P$ admits, pertaining to the order that edges are visited and to how the cost of the minimum makespan $R_n(P)$ is determined, for $n\in \{1,2,3\}$. These subdivisions are the starting points for our main result, which is to study makespan trade-offs with respect to the size of the fleet. In particular, we define $ R_{n,m} (\Delta)= \max_{P \in \Delta} R_n(P)/R_m(P)$, and we prove that, over all non-obtuse triangles $\Delta$: (i) $R_{1,3}(\Delta)$ ranges from $\sqrt{10}$ to $4$, (ii) $R_{2,3}(\Delta)$ ranges from $\sqrt{2}$ to $2$, and (iii) $R_{1,2}(\Delta)$ ranges from $5/2$ to $3$. In every case, we pinpoint the starting points within every triangle $\Delta$ that maximize $R_{n,m} (\Delta)$, as well as we identify the triangles that determine all $\inf_\Delta R_{n,m}(\Delta)$ and $\sup_\Delta R_{n,m}(\Delta)$ over the set of non-obtuse triangles.