title: Near-Field Scanning Microwave Microscopy in the Single Photon Regime

abstract: The microwave properties of nano-scale structures are important in a wide variety of applications in quantum technology. Here we describe a low-power cryogenic near-field scanning microwave microscope (NSMM) which maintains nano-scale dielectric contrast down to the single microwave photon regime, up to $10^{9}$ times lower power than in typical NSMMs. We discuss the remaining challenges towards developing nano-scale NSMM for quantum coherent interaction with two-level systems as an enabling tool for the development of quantum technologies in the microwave regime.


title: Negative Friction Coefficients in Superlubric Graphite-Hexagonal Boron Nitride Heterojunctions

abstract: Negative friction coefficients, where friction is reduced upon increasing normal load, are predicted for superlubric graphite-hexagonal boron nitride heterojunctions. The origin of this counterintuitive behavior lies in the load-induced suppression of the moir\'e superstructure out-of-plane distortions leading to a less dissipative interfacial dynamics. Thermally induced enhancement of the out-of-plane fluctuations leads to an unusual increase of friction with temperature. The highlighted frictional mechanism is of a general nature and is expected to appear in many layered material heterojunctions.


title: Towards Reliable, Automated General Movement Assessment for Perinatal Stroke Screening in Infants Using Wearable Accelerometers

abstract: Perinatal stroke (PS) is a serious condition that, if undetected and thus untreated, often leads to life-long disability, in particular Cerebral Palsy (CP). In clinical settings, Prechtl's General Movement Assessment (GMA) can be used to classify infant movements using a Gestalt approach, identifying infants at high risk of developing PS. Training and maintenance of assessment skills are essential and expensive for the correct use of GMA, yet many practitioners lack these skills, preventing larger-scale screening and leading to significant risks of missing opportunities for early detection and intervention for affected infants. We present an automated approach to GMA, based on body-worn accelerometers and a novel sensor data analysis method-Discriminative Pattern Discovery (DPD)-that is designed to cope with scenarios where only coarse annotations of data are available for model training. We demonstrate the effectiveness of our approach in a study with 34 newborns (21 typically developing infants and 13 PS infants with abnormal movements). Our method is able to correctly recognise the trials with abnormal movements with at least the accuracy that is required by newly trained human annotators (75%), which is encouraging towards our ultimate goal of an automated PS screening system that can be used population-wide.


title: With Great Speed Come Small Buffers: Space-Bandwidth Tradeoffs for Routing

abstract: We consider the Adversarial Queuing Theory (AQT) model, where packet arrivals are subject to a maximum average rate $0\le\rho\le1$ and burstiness $\sigma\ge0$. In this model, we analyze the size of buffers required to avoid overflows in the basic case of a path. Our main results characterize the space required by the average rate and the number of distinct destinations: we show that $O(k d^{1/k})$ space suffice, where $d$ is the number of distinct destinations and $k=\lfloor 1/\rho \rfloor$; and we show that $\Omega(\frac 1 k d^{1/k})$ space is necessary. For directed trees, we describe an algorithm whose buffer space requirement is at most $1 + d' + \sigma$ where $d'$ is the maximum number of destinations on any root-leaf path.


title: Strategyproof Facility Location for Three Agents on a Circle

abstract: We consider the facility location problem in a metric space, focusing on the case of three agents. We show that selecting the reported location of each agent with probability proportional to the distance between the other two agents results in a mechanism that is strategyproof in expectation, and dominates the random dictator mechanism in terms of utilitarian social welfare. We further improve the upper bound for three agents on a circle to 7/6 (whereas random dictator obtains 4/3); and provide the first lower bounds for randomized strategyproof facility location in any metric space, using linear programming.


title: Reconfiguration Graph for Vertex Colourings of Weakly Chordal Graphs

abstract: The reconfiguration graph $R_k(G)$ of the $k$-colourings of a graph $G$ contains as its vertex set the $k$-colourings of $G$ and two colourings are joined by an edge if they differ in colour on just one vertex of $G$. We show that for each $k \geq 3$ there is a $k$-colourable weakly chordal graph $G$ such that $R_{k+1}(G)$ is disconnected. We also introduce a subclass of $k$-colourable weakly chordal graphs which we call $k$-colourable compact graphs and show that for each $k$-colourable compact graph $G$ on $n$ vertices, $R_{k+1}(G)$ has diameter $O(n^2)$. We show that this class contains all $k$-colourable co-chordal graphs and when $k = 3$ all $3$-colourable $(P_5, \overline{P_5}, C_5)$-free graphs. We also mention some open problems.


title: Radiation Efficiency Aware High-Mobility Massive MIMO with Antenna Selection

abstract: High-mobility wireless communications have received a lot of attentions in the past few years. In this paper, we consider angle-domain Doppler shifts compensation scheme to reduce the channel time variation for the high-mobility uplink from a high-speed terminal (HST) to a base station (BS). We propose to minimize Doppler spread by antenna weighting under the constraint of maintaining radiation efficiency. The sequential parametric convex approximation (SPCA) algorithm is exploited to solve the above non-convex problem. Moreover, in order to save the RF chains, we further exploit the idea of antenna selection for high-mobility wireless communications. Simulations verify the effectiveness of the proposed studies.


title: Ultralow-field Nuclear Magnetic Resonance Asymmetric Spectroscopy

abstract: Ultralow-field nuclear magnetic resonance (NMR) provides a new regime for many applications ranging from materials science to fundamental physics. However, the experimentally observed spectra show asymmetric amplitudes, differing greatly from those predicted by the standard theory. Its physical origin remains unclear, as well as how to suppress it. Here we provide a comprehensive model to explain the asymmetric spectral amplitudes, further observe more unprecedented asymmetric spectroscopy and find a way to eliminate it. Moreover, contrary to the traditional idea that asymmetric phenomena were considered as a nuisance, we show that more information can be gained from the asymmetric spectroscopy, e.g., the light shift of atomic vapors and the sign of Land$\acute{\textrm{e}}$ $g$ factor of NMR systems.


title: A set-oriented path following method for the approximation of parameter dependent attractors

abstract: In this work we present a set-oriented path following method for the computation of relative global attractors of parameter-dependent dynamical systems. We start with an initial approximation of the relative global attractor for a fixed parameter $\lambda_0$ computed by a set-oriented subdivision method. By using previously obtained approximations of the parameter-dependent relative global attractor we can track it with respect to a one-dimensional parameter $\lambda > \lambda_0$ without restarting the whole subdivision procedure. We illustrate the feasibility of the set-oriented path following method by exploring the dynamics in models for shear flows during the transition to turbulence.


title: Nonlocal Cosmology II --- Cosmic acceleration without fine tuning or dark energy

abstract: We present an improved version of our original cosmological model to explain the current phase of cosmological acceleration without resorting to a cosmological constant or any other mass scale. Like the original, this phenomenological approach is based on an effective quantum gravitational action, but now depends on the original nonlocal dimensionless scalar $X = \square^{-1} R$ only through $Y = \square^{-1} g^{\mu\nu} X_{,\mu} X_{,\nu}$. Both $X$ and $Y$ are quiescent during the radiation-dominated ($R=0$) era, both only grow logarithmically during matter dominance, and neither affects the propagation of gravitational radiation. However, while $X$ has the same sign for gravitationally bound systems as for cosmology, we show that the sign of $Y$ differs for the two cases: it is positive for cosmology and negative for strongly gravitationally bound systems. We can therefore enforce the $\Lambda$CDM expansion history by making a suitable choice of the nonlocal distortion function $f(Y)$ for $Y > 0$, while ensuring that there is no change in the heavily constrained solar system phenomenology simply by making $f$ vanish for $Y < 0$ without discontinuity. The required $f(Y>0)$ is determined numerically to have a strikingly simple exponential form.


title: Calibrating gravitational-wave detectors with GW170817

abstract: The waveform of a compact binary coalescence is predicted by general relativity. It is therefore possible to directly constrain the response of a gravitational-wave (GW) detector by analyzing a signal's observed amplitude and phase evolution as a function of frequency. GW signals alone constrain the relative amplitude and phase between different frequencies within the same detector and between different detectors. We analyze GW170817's ability to calibrate the LIGO/Virgo detectors, finding a relative amplitude calibration precision of approximately $\pm20\%$ and relative phase precision of $\pm15^\circ$ (1-$\sigma$ uncertainty) between the LIGO Hanford and Livingston detectors. Incorporating additional information about the distance and inclination of the source from electromagnetic observations, the relative amplitude of the LIGO detectors can be tightened to $\sim\pm15\%$. We investigate the ability of future events to improve astronomical calibration. By simulating the cumulative uncertainties from an ensemble of detections, we find that with several hundred events with electromagnetic counterparts, or several thousand events without counterparts, we reach percent-level astronomical calibration. This corresponds to $\sim$5-10 years of operation at advanced LIGO and Virgo design sensitivity. It is to be emphasized that direct {\em in-situ}\/ measurements of detector calibration provide significantly higher precision than astronomical sources, and already constrain the calibration to a few percent in amplitude and a few degrees in phase. In this sense, our astronomical calibrators only corroborate existing calibration measurements. Nonetheless, astrophysical calibration may become an important corroboration of existing calibration methods, providing a completely independent constraint of potential systematics.


title: Breaking the Softmax Bottleneck via Learnable Monotonic Pointwise Non-linearities

abstract: The Softmax function on top of a final linear layer is the de facto method to output probability distributions in neural networks. In many applications such as language models or text generation, this model has to produce distributions over large output vocabularies. Recently, this has been shown to have limited representational capacity due to its connection with the rank bottleneck in matrix factorization. However, little is known about the limitations of Linear-Softmax for quantities of practical interest such as cross entropy or mode estimation, a direction that we explore here. As an efficient and effective solution to alleviate this issue, we propose to learn parametric monotonic functions on top of the logits. We theoretically investigate the rank increasing capabilities of such monotonic functions. Empirically, our method improves in two different quality metrics over the traditional Linear-Softmax layer in synthetic and real language model experiments, adding little time or memory overhead, while being comparable to the more computationally expensive mixture of Softmaxes.


title: A fast linearized finite difference method for the nonlinear multi-term time-fractional wave equation

abstract: In this paper, we study a fast and linearized finite difference method to solve the nonlinear time-fractional wave equation with multi fractional orders. We first propose a discretization to the multi-term Caputo derivative based on the recently established fast L2-1{\sigma} formula and a weighted approach. Then we apply the discretization to construct a fully fast linearized discrete scheme for the nonlinear problem under consideration. The nonlinear term, which just fulfills the Lipschitz condition, will be evaluated on the previous time level. Therefore only linear systems are needed to be solved for obtaining numerical solutions. The proposed scheme is shown to have second-order unconditional convergence with respect to the discrete H1-norm. Numerical examples are provided to justify the efficiency.


title: Curve-shortening flow of open, elastic curves in $\mathbb{R}^2$ with repelling endpoints: A minimizing movement approach

abstract: We study an $L^{2}$-type gradient flow of an immersed elastic curve in $\mathbb{R}^{2}$ whose endpoints repel each other via a Coulomb potential. By De Giorgi's minimizing movements scheme we prove long-time existence of the flow. The work is complemented by several numerical experiments.


title: Surface acoustic wave modulation of single photon emission from GaN/InGaN nanowire quantum dots

abstract: On-chip quantum information processing requires controllable quantum light sources that can be operated on-demand at high-speeds and with the possibility of in-situ control of the photon emission wavelength and its optical polarization properties. Here, we report on the dynamic control of the optical emission from core-shell GaN/InGaN nanowire (NW) heterostructures using radio frequency surface acoustic waves (SAWs). The SAWs are excited on the surface of a piezoelectric lithium niobate crystal equipped with a SAW delay line onto which the NWs were mechanically transferred. Luminescent quantum dot (QD)-like exciton localization centers induced by compositional fluctuations within the InGaN nanoshell were identified using stroboscopic micro-photoluminescence (micro-PL) spectroscopy. They exhibit narrow and almost fully linearly polarized emission lines in the micro-PL spectra and a pronounced anti-bunching signature of single photon emission in the photon correlation experiments. When the nanowire is perturbed by the propagating SAW, the embedded QD is periodically strained and its excitonic transitions are modulated by the acousto-mechanical coupling, giving rise to a spectral fine-tuning within a ~1.5 meV bandwidth at the acoustic frequency of ~330 MHz. This outcome can be further combined with spectral detection filtering for temporal control of the emitted photons. The effect of the SAW piezoelectric field on the QD charge population and on the optical polarization degree is also observed. The advantage of the acousto-optoelectric over other control schemes is that it allows in-situ manipulation of the optical emission properties over a wide frequency range (up to GHz frequencies).


title: DeepHoops: Evaluating Micro-Actions in Basketball Using Deep Feature Representations of Spatio-Temporal Data

abstract: How much is an on-ball screen worth? How much is a backdoor cut away from the ball worth? Basketball is one of a number of sports which, within the past decade, have seen an explosion in quantitative metrics and methods for evaluating players and teams. However, it is still challenging to evaluate individual off-ball events in terms of how they contribute to the success of a possession. In this study, we develop an end-to-end deep learning architecture DeepHoops to process a unique dataset composed of spatio-temporal tracking data from NBA games in order to generate a running stream of predictions on the expected points to be scored as a possession progresses. We frame the problem as a multi-class sequence classification problem in which our model estimates probabilities of terminal actions taken by players (e.g. take field goal, turnover, foul etc.) at each moment of a possession based on a sequence of ball and player court locations preceding the said moment. Each of these terminal actions is associated with an expected point value, which is used to estimate the expected points to be scored. One of the challenges associated with this problem is the high imbalance in the action classes. To solve this problem, we parameterize a downsampling scheme for the training phase. We demonstrate that DeepHoops is well-calibrated, estimating accurately the probabilities of each terminal action and we further showcase the model's capability to evaluate individual actions (potentially off-ball) within a possession that are not captured by boxscore statistics.


title: Rejuvenation and Memory Effects in a Structural Glass

abstract: We show numerically that a three-dimensional model for structural glass displays aging, rejuvenation and memory effects when submitted to a temperature cycle. These effects indicate that the free energy landscape of structural glasses may possess the complex hierarchical structure that characterize materials such as spin and polymer glasses. We use the theoretical concept of marginal stability to interpret our results, and explain in which physical conditions a complex aging dynamics can emerge in dense supercooled liquids, paving the way for future experimental studies of complex aging dynamics in colloidal and granular glasses.


title: TREVR: A general $N log^2 N$ radiative transfer algorithm

abstract: We present TREVR (Tree-based REVerse Ray Tracing), a general algorithm for computing the radiation field, including absorption, in astrophysical simulations. TREVR is designed to handle large numbers of sources and absorbers; it is based on a tree data structure and is thus suited to codes that use trees for their gravity or hydrodynamics solvers (e.g. Adaptive Mesh Refinement). It achieves computational speed while maintaining a specified accuracy via controlled lowering of the resolution of both sources and rays from each source. TREVR computes the radiation field in order ${N\log N_{source}}$ time without absorption and order ${N \log N_{source} \log{N}}$ time with absorption. These scalings arise from merging sources of radiation according to an opening angle criterion and walking the tree structure to trace a ray to a depth that gives the chosen accuracy for absorption. The absorption-depth refinement criterion is unique to TREVR. We provide a suite of tests demonstrating the algorithm's ability to accurately compute fluxes, ionization fronts and shadows.


title: Smooth approximation is not a selection principle for the transport equation with rough vector field

abstract: In this paper we analyse the selection problem for weak solutions of the transport equation with rough vector field. We answer in the negative the question whether solutions of the equation with a regularized vector field converge to a unique limit, which would be the selected solution of the limit problem. To this aim, we give a new example of a vector field which admits infinitely many flows. Then we construct a smooth approximating sequence of the vector field for which the corresponding solutions have subsequences converging to different solutions of the limit equation.


title: Variable stars in Palomar 13; an evaporating globular cluster

abstract: We present new CCD $VI$ photometry of the distant globular cluster Pal 13. Fourier decomposition of the light curves of the three cluster member RRab stars lead to estimations of [Fe/H]=-1.65, and a distance of 23.67$\pm$0.57 kpc. Light and colour near minimum phases for RRab stars leads to an estimate of $E(B-V)$=0.104 $\pm$ 0.001. A $V/(V-I)$ colour-magnitude diagram, built exclusively with likely star members, shows consistency with the above parameters and an age of 12 Gyrs. A search of variable stars in the field of view of our images revealed the variability of a red giant cluster member and of three probably non-member stars; two RRab stars and one W Virginis star or CW. The GAIA proper motions of member stars in Pal 13 show a significant scatter, consistent with the scenario of the cluster being tidally stripped.


title: The Arboricity Captures the Complexity of Sampling Edges

abstract: In this paper, we revisit the problem of sampling edges in an unknown graph $G = (V, E)$ from a distribution that is (pointwise) almost uniform over $E$. We consider the case where there is some a priori upper bound on the arboriciy of $G$. Given query access to a graph $G$ over $n$ vertices and of average degree $d$ and arboricity at most $\alpha$, we design an algorithm that performs $O\!\left(\frac{\alpha}{d} \cdot \frac{\log^3 n}{\varepsilon}\right)$ queries in expectation and returns an edge in the graph such that every edge $e \in E$ is sampled with probability $(1 \pm \varepsilon)/m$. The algorithm performs two types of queries: degree queries and neighbor queries. We show that the upper bound is tight (up to poly-logarithmic factors and the dependence in $\varepsilon$), as $\Omega\!\left(\frac{\alpha}{d} \right)$ queries are necessary for the easier task of sampling edges from any distribution over $E$ that is close to uniform in total variational distance. We also prove that even if $G$ is a tree (i.e., $\alpha = 1$ so that $\frac{\alpha}{d}=\Theta(1)$), $\Omega\left(\frac{\log n}{\log\log n}\right)$ queries are necessary to sample an edge from any distribution that is pointwise close to uniform, thus establishing that a $\mathrm{poly}(\log n)$ factor is necessary for constant $\alpha$. Finally we show how our algorithm can be applied to obtain a new result on approximately counting subgraphs, based on the recent work of Assadi, Kapralov, and Khanna (ITCS, 2019).


title: Assessing Olfaction Using Ultrasonic Vocalization Recordings in Mouse Pups with a Sono-olfactometer

abstract: Olfaction is the first sensory modality to develop during fetal life in mammals, and plays a key role in the various behaviors of neonates such as feeding and social interaction. Odorant cues (i.e., mother or predator scents) can trigger potentiation or inhibition of ultrasonic vocalizations (USV) emitted by pups following their isolation. Here, we report how USV are inhibited by olfactory cues using a sono-olfactometer that has been designed to quantify precisely olfaction in pups congenitally infected by cytomegalovirus. This olfactory-driven behavioral test assesses the USV emitted in presence of unfamiliar odorants such as citral scent or adult male mouse scent. We measure the number of USV emitted as an index of odorant detection during the three periods of the 5-min isolation time of the pup into the sono-olfactometer: first period without any odorant, second period with odorant exposure and last period with exhaust odorant. This protocol can be easily used to reveal olfactory deficits in pups with altered olfactory system due to toxic lesions or infectious diseases.


title: Cyclic annealing as an iterated random map

abstract: Disordered magnets, martensitic mixed crystals, and glassy solids can be irreversibly deformed by subjecting them to external deformation. The deformation produces a smooth, reversible response punctuated by abrupt relaxation "glitches". Under appropriate repeated forward and reverse deformation producing multiple glitches, a strict repetition of a single sequence of microscopic configurations often emerges. We exhibit these features by describing the evolution of the system configuration from glitch to glitch as a mapping of $\mathcal{N}$ states into one-another. A map $\mathbf{U}$ controls forward deformation; a second map $\mathbf{D}$ controls reverse deformation. Iteration of a given sequence of forward and reverse maps, e.g. $\mathbf{DDDDUUU}$ necessarily produces a convergence to a fixed cyclic repetition of states covering multiple glitches. The repetition may have a period of more than one strain cycle, as recently observed in simulations. Using numerical sampling, we characterize the convergence properties of four types of random maps implementing successive physical restrictions. The most restrictive is the much-studied Preisach model. These maps show only the most qualitative resemblance to annealing simulations. However, they suggest further properties needed for a realistic mapping scheme.


title: A free-energy stable nodal discontinuous Galerkin approximation with summation-by-parts property for the Cahn-Hilliard equation

abstract: We present a nodal Discontinuous Galerkin (DG) scheme for the Cahn-Hilliard equation that satisfies the summation-by-parts simultaneous-approximation-term (SBP-SAT) property. The latter permits us to show that the discrete free-energy is bounded, and as a result, the scheme is provably stable. The scheme and the stability proof are presented for general curvilinear three-dimensional hexahedral meshes. We use the Bassi-Rebay 1 (BR1) scheme to compute interface fluxes, and an IMplicit-EXplicit (IMEX) scheme to integrate in time. Lastly, we test the theoretical findings numerically and present examples for two and three-dimensional problems.


title: Clifford-based spectral action and renormalization group analysis of the gauge couplings

abstract: The Spectral Action Principle in noncommutative geometry derives the actions of the Standard Model and General Relativity (along with several other gravitational terms) by reconciling them in a geometric setting, and hence offers an explanation for their common origin. However, one of the requirements in the minimal formalism, unification of the gauge coupling constants, is not satisfied, since the basic construction does not introduce anything new that can change the renormalization group (RG) running of the Standard Model. On the other hand, it has been recently argued that incorporating structure of the Clifford algebra into the finite part of the spectral triple, the main object that encodes the complete information of a noncommutative space, gives rise to five additional scalar fields in the basic framework. We investigate whether these scalars can help to achieve unification. We perform a RG analysis at the one-loop level, allowing possible mass values of these scalars to float from the electroweak scale to the putative unification scale. We show that out of twenty configurations of mass hierarchy in total, there does not exist even a single case that can lead to unification. In consequence, we confirm that the spectral action formalism requires a model-construction scheme beyond the (modified) minimal framework.


title: On the qubit routing problem

abstract: We introduce a new architecture-agnostic methodology for mapping abstract quantum circuits to realistic quantum computing devices with restricted qubit connectivity, as implemented by Cambridge Quantum Computing's tket compiler. We present empirical results showing the effectiveness of this method in terms of reducing two-qubit gate depth and two-qubit gate count, compared to other implementations.


title: A Fast Approximate Approach to Microlensing Survey Analysis

abstract: Microlensing can be used to discover exoplanets of a wide range of masses with orbits beyond ~ 1 AU, and even free-floating planets. The WFIRST mission will use microlensing to discover approximately 1600 planets by monitoring ~100 million stars to find ~50000 microlensing events. Modelling each microlensing event, especially the ones involving two or more lenses, is typically complicated and time-consuming, and analyzing thousands of WFIRST microlensing events is possibly infeasible using current methods. Here, we present an algorithm that is able to rapidly evaluate thousands of simulated WFIRST binary-lens microlensing light curves, returning an estimate for the physical parameters of the lens systems. We find that this algorithm can recover projected separations between the planet and the star very well for low-mass-ratio events, and can also estimate mass ratios within an order of magnitude for events with wide and close caustic topologies.


title: Unsupervised Grounding of Plannable First-Order Logic Representation from Images

abstract: Recently, there is an increasing interest in obtaining the relational structures of the environment in the Reinforcement Learning community. However, the resulting "relations" are not the discrete, logical predicates compatible to the symbolic reasoning such as classical planning or goal recognition. Meanwhile, Latplan (Asai and Fukunaga 2018) bridged the gap between deep-learning perceptual systems and symbolic classical planners. One key component of the system is a Neural Network called State AutoEncoder (SAE), which encodes an image-based input into a propositional representation compatible to classical planning. To get the best of both worlds, we propose First-Order State AutoEncoder, an unsupervised architecture for grounding the first-order logic predicates and facts. Each predicate models a relationship between objects by taking the interpretable arguments and returning a propositional value. In the experiment using 8-Puzzle and a photo-realistic Blocksworld environment, we show that (1) the resulting predicates capture the interpretable relations (e.g. spatial), (2) they help obtaining the compact, abstract model of the environment, and finally, (3) the resulting model is compatible to symbolic classical planning.


title: Interrupted Motility Induced Phase Separation in Aligning Active Colloids

abstract: Switching on high activity in a relatively dense system of active Janus colloids, we observe fast clustering, followed by cluster aggregation towards full phase separation. The phase separation process is however interrupted when large enough clusters start breaking apart. Following the cluster size distribution as a function of time, we identify three successive dynamical regimes. Tracking both the particle positions and orientations, we characterize the structural ordering and alignment in the growing clusters and thereby unveil the mechanisms at play in these regimes. In particular, we identify how alignment between the neighboring particles is responsible for the interruption of the full phase separation. This experimental study, which provides the first large scale observation of the phase separation kinetics in active colloids, combined with single particle analysis of the local mechanisms, points towards the new physics observed when both alignment and short-range repulsions are present.


title: Fractional Conductance in Strongly Interacting 1D Systems

abstract: We study one dimensional clean systems with few channels and strong electron-electron interactions. We find that in several circumstances, even when time reversal symmetry holds, they may lead to two terminal fractional quantized conductance and fractional shot noise. The condition on the commensurability of the Fermi momenta of the different channels and the strength of interactions resulting in such remarkable phenomena are explored using abelian bosonization. Finite temperature and length effects are accounted for by a generalization of the Luther-Emery re-fermionization at specific values of the interaction strength. We discuss the connection of our model to recent experiments in confined 2DEG, featuring possible fractional conductance plateaus. One of the most dominant observed fractions, with two terminal conductance equals to $\frac{2}{5}\frac{e^{2}}{h}$, is found in several scenarios of our model. Finally, we discuss how at very small energy scales the conductance returns to an integer value and the role of disorder.


title: Boosted Top quark polarization

abstract: In top quark production, the polarization of top quarks, decided by the chiral structure of couplings, is likely to be modified in the presence of any new physics contribution to the production. Hence the same is a good discriminator for those new physics models wherein the couplings have a chiral structure different than that in the Standard Model (SM). In this note we construct probes of the polarization of a top quark decaying hadronically, using easily accessible kinematic variables such as the energy fraction or angular correlations of the decay products. Tagging the boosted top quark using the usual jet sub structure technique we study robustness of these observables for a benchmark process, $W^{\prime} \to tb$. We demonstrate that the energy fraction of b-jet in the laboratory frame and a new angular variable, constructed by us in the top rest frame, are both very powerful tools to discriminate between the left and right polarized top quarks. Based on the polarization sensitive angular variables, we construct asymmetries which reflect the polarization. We study the efficiency of these variables for two new physics processes where which give rise to boosted top quarks: (i) decay of the top squark in the context of supersymmetry searches, and (ii) decays of the Kaluza-Klein(KK) graviton and KK gluon, in Randall Sundrum(RS) model. Remarkably, it is found that the asymmetry can vary over a wide range about +20\% to -20\%. The dependence of asymmetry on top quark couplings of the new particles present in these models beyond the SM (BSM) is also investigated in detail.


title: Structure and equation-of-state of a disordered system of shape anisotropic patchy particles

abstract: We present a new two dimensional model for elliptic (i.e., shape anisotropic) patchy colloids, where the impenetrable core of the particles is decorated on its co-)vertices by Kern-Frenkel type of patches. Using (i) well-documented criteria for the overlap of the undecorated ellipses and (ii) proposing new criteria that give evidence if the patchy regions of two particles interact we perform extensive Monte Carlo simulations in the NPT ensemble. Considering elliptic particles of aspect ratios \{kappa} = 2, 4, and 6 with two patches located at the co-vertices and choosing representative values for the temperature we study on a semi-quantitative level the emerging disordered phases in terms of snapshots and radial distribution functions and present results for the equation-of-state.


title: Spatial And Temporal Changes Of The Geomagnetic Field: Insights From Forward And Inverse Core Field Models

abstract: Observational constraints on geomagnetic field changes from interannual to millenial periods are reviewed, and the current resolution of field models (covering archeological to satellite eras) is discussed. With the perspective of data assimilation, emphasis is put on uncertainties entaching Gauss coefficients, and on the statistical properties of ground-based records. These latter potentially call for leaving behind the notion of geomagnetic jerks. The accuracy at which we recover interannual changes also requires considering with caution the apparent periodicity seen in the secular acceleration from satellite data. I then address the interpretation of recorded magnetic fluctuations in terms of core dynamics, highlighting the need for models that allow (or pre-suppose) a magnetic energy orders of magnitudes larger than the kinetic energy at large length-scales, a target for future numerical simulations of the geodynamo. I finally recall the first attempts at implementing geomagnetic data assimilation algorithms.


title: Monodromy of rational curves on toric surfaces

abstract: For an ample line bundle $\mathcal{L}$ on a complete toric surface $X$, we consider the subset $V_{\mathcal{L}} \subset \vert \mathcal{L} \vert$ of irreducible, nodal, rational curves contained in the smooth locus of $X$. We study the monodromy map from the fundamental group of $V_{\mathcal{L}}$ to the permutation group on the set of nodes of a reference curve $C \in V_{\mathcal{L}}$. We identify a certain obstruction map $\varPsi_{X}$ defined on the set of nodes of $C$ and show that the image of the monodromy is exactly the group of deck transformations of $\varPsi_{X}$, provided that $\mathcal{L}$ is sufficiently big (in a sense we precise below). Along the way, we provide a handy tool to compute the image of the monodromy for any pair $(X, \mathcal{L})$. Eventually, we present a family of pairs $(X, \mathcal{L})$ with small $\mathcal{L}$ and for which the image of the monodromy is strictly smaller than expected.


title: Ab-initio calculations of carbon and boron nitride allotropes and their structural phase transitions using periodic coupled cluster theory

abstract: We present an ab-initio study of boron nitride as well as carbon allotropes. Their relative thermodynamic stabilities and structural phase transitions from low- to high-density phases are investigated. Pressure-temperature phase diagrams are calculated and compared to experimental findings. The calculations are performed using quantum chemical wavefunction based as well as density functional theories. Our findings reveal that predicted energy differences often depend significantly on the choice of the employed method. Comparison between calculated and experimental results allows for benchmarking the accuracy of various levels of theory. The produced results show that quantum chemical wavefunction based theories allow for achieving systematically improvable estimates. We find that on the level of coupled cluster theories the low- and high-density phases of boron nitride become thermodynamically degenerate at 0 K. This is in agreement with recent experimental findings, indicating that cubic boron nitride is not the thermodynamically stable allotrope at ambient conditions. Furthermore we employ the calculated results to assess transition probabilities from graphitic low-density to diamond-like high-density phases in an approximate manner. We conclude that the stacking order of the parent graphitic material is crucial for the possible formation of meta-stable wurtzite boron nitride and hexagonal carbon diamond also known as lonsdaleite.


title: The inviscid limit of Navier-Stokes equations for vortex-wave data on $\mathbb{R}^2$

abstract: We establish the inviscid limit of the incompressible Navier-Stokes equations on the whole plane $\mathbb{R}^2$ for initial data having vorticity as a superposition of point vortices and a regular component. In particular, this rigorously justifies the vortex-wave system from the physical Navier-Stokes flows in the vanishing viscosity limit, a model that was introduced by Marchioro and Pulvirenti in the early 90s to describe the dynamics of point vortices in a regular ambient vorticity background. The proof rests on the previous analysis of Gallay in his derivation of the vortex-point system.


title: Statistics and Samples in Distributional Reinforcement Learning

abstract: We present a unifying framework for designing and analysing distributional reinforcement learning (DRL) algorithms in terms of recursively estimating statistics of the return distribution. Our key insight is that DRL algorithms can be decomposed as the combination of some statistical estimator and a method for imputing a return distribution consistent with that set of statistics. With this new understanding, we are able to provide improved analyses of existing DRL algorithms as well as construct a new algorithm (EDRL) based upon estimation of the expectiles of the return distribution. We compare EDRL with existing methods on a variety of MDPs to illustrate concrete aspects of our analysis, and develop a deep RL variant of the algorithm, ER-DQN, which we evaluate on the Atari-57 suite of games.


title: Searches for third-generation scalar leptoquarks in $\sqrt{s}$ = 13 TeV pp collisions with the ATLAS detector

abstract: Limits are set on the pair production of scalar leptoquarks, where all possible decays of the leptoquark into a quark ($t$, $b$) and a lepton ($\tau$, $\nu$) of the third generation are considered. The limits are presented as a function of the leptoquark mass and the branching ratio into charged leptons for up-type ($\mathrm{LQ}_{3}^{u}$ $\rightarrow t \nu / b \tau$) and down-type ($\mathrm{LQ}_{3}^{d}$ $\rightarrow b \nu / t \tau$) leptoquarks. Many results are reinterpretations of previously published ATLAS searches. In all cases, LHC proton-proton collision data at a centre-of-mass energy of $\sqrt{s}$ = 13 TeV recorded by the ATLAS detector in 2015 and 2016 are used, corresponding to an integrated luminosity of 36.1 fb$^{-1}$. Masses below 800 GeV are excluded for both $\mathrm{LQ}_{3}^{u}$ and $\mathrm{LQ}_{3}^{d}$ independently of the branching ratio, with masses below about 1 TeV being excluded for the limiting cases of branching ratios equal to zero or unity.


title: Really Small Shoe Boxes - On Realistic Quantum Resource Estimation

abstract: Reliable resource estimation and benchmarking of quantum algorithms is a critical component of the development cycle of viable quantum applications for quantum computers of all sizes. Determining critical resource bottlenecks in algorithms, especially when resource intensive error correction protocols are required, will be crucial to reduce the cost of implementing viable algorithms on actual quantum hardware.


title: Formation and eruption of sigmoidal structure from a weak field region of NOAA 11942

abstract: Using observations from Solar Dynamics Observatory, we studied an interesting example of a sigmoid formation and eruption from small-scale flux canceling regions of active region (AR) 11942. Analysis of HMI and AIA observations infer that initially the AR is compact and bipolar in nature, evolved to sheared configuration consisting of inverse J-shaped loops hosting a filament channel over a couple of days. By tracking the photospheric magnetic features, shearing and converging motions are observed to play a prime role in the development of S-shaped loops and further flux cancellation leads to tether-cutting reconnection of J-loops. This phase is co-temporal with the filament rise motion followed by sigmoid eruption at 21:32 UT on January 6. The flux rope rises in phases of slow (v$_{avg}$ = 26 km~s$^{-1}$) and fast (a$_{avg}$= 55 ms$^{-2}$) rise motion categorizing the CME as slow with an associated weak C1.0 class X-ray flare. The flare ribbon separation velocity peaks at around peak time of the flare at which maximum reconnection rate (2.14 Vcm$^{-1}$) occurs. Further, the EUV light-curves of 131, 171\AA~have delayed peaks of 130 minutes compared to 94\AA~and is explained by differential emission measure. Our analysis suggests that the energy release is proceeded in a much long time duration, manifesting the onset of filament rise and eventual eruption driven by converging and canceling flux in the photosphere. Unlike strong eruption events, the observed slow CME and weak flare are indications of slow runway tether-cutting reconnection where most of the sheared arcade is relaxed during the extended post phase of the eruption.


title: Existence of densities for stochastic evolution equations driven by fractional Brownian motion

abstract: In this work, we prove a version of H\"{o}rmander's theorem for a stochastic evolution equation driven by a trace-class fractional Brownian motion with Hurst exponent $\frac{1}{2} < H < 1$ and an analytic semigroup on a given separable Hilbert space. In contrast to the classical finite-dimensional case, the Jacobian operator in typical solutions of parabolic stochastic PDEs is not invertible which causes a severe difficulty in expressing the Malliavin matrix in terms of an adapted process. Under a H\"{o}rmander's bracket condition and some algebraic constraints on the vector fields combined with the range of the semigroup, we prove the law of finite-dimensional projections of such solutions has a density w.r.t Lebesgue measure. The argument is based on rough path techniques and a suitable analysis on the Gaussian space of the fractional Brownian motion.


title: A Fully Lagrangian Meshfree Framework for PDEs on Evolving Surfaces

abstract: We propose a novel framework to solve PDEs on moving manifolds, where the evolving surface is represented by a moving point cloud. This has the advantage of avoiding the need to discretize the bulk volume around the surface, while also avoiding the need to have a global mesh. Distortions in the point cloud as a result of the movement are fixed by local adaptation. We first establish a comprehensive Lagrangian framework for arbitrary movement of curves and surfaces given by point clouds. Collision detection algorithms between point cloud surfaces are introduced, which also allow the handling of evolving manifolds with topological changes. We then couple this Lagrangian framework with a meshfree Generalized Finite Difference Method (GFDM) to approximate surface differential operators, which together give a method to solve PDEs on evolving manifolds. The applicability of this method is illustrated with a range of numerical examples, which include advection-diffusion equations with large deformations of the surface, curvature dependent geometric motion, and wave equations on evolving surfaces.


title: Addressing the ${\mathbf R_{D^{(*)}}}$ anomalies with an ${\mathbf S_1}$ leptoquark from $\mathbf{SO(10)}$ grand unification

abstract: Motivated by the $R_{D^{(*)}}$ anomalies, we investigate an $\mathrm{SO}(10)$ grand unification scenario where a charge $-1/3$ scalar leptoquark ($S_1$) remains as the only new physics candidate at the TeV scale. This leptoquark along with the Standard Model (SM) Higgs doublet originates from the same 10-dimensional real scalar multiplet in the $\mathrm{SO}(10)$ grand unification framework taking its mass close to the electroweak scale. We explicitly show how the gauge coupling unification is achieved with only one intermediate symmetry-breaking scale at which the Pati-Salam gauge group is broken into the SM group. We investigate the phenomenological implications of our scenario and show that an $S_1$ with a specific Yukawa texture can explain the $R_{D^{(*)}}$ anomalies. We perform a multiparameter scan considering the relevant flavour constraints on $R_{D^{(*)}}$, $F_L{(D^*)}$, $P_\tau(D^*)$ and $R_{K^{(*)}}^{\nu\nu}$ as well as the constraint coming from the $Z\to\tau\tau$ decay and the latest $\tau\tau$ resonance search data at the LHC. Our analysis shows that a single leptoquark solution to the observed $R_{D^{(*)}}$ anomalies with $S_1$ is still a viable solution.


title: The fluctuations of the giant cluster for percolation on random split trees

abstract: A split tree of cardinality $n$ is constructed by distributing $n$ "balls" in a subset of vertices of an infinite tree which encompasses many types of random trees such as $m$-ary search trees, quad trees, median-of-$(2k+1)$ trees, fringe-balanced trees, digital search trees and random simplex trees. In this work, we study Bernoulli bond percolation on arbitrary split trees of large but finite cardinality $n$. We show for appropriate percolation regimes that depend on the cardinality $n$ of the split tree that there exists a unique giant cluster, the fluctuations of the size of the giant cluster as $n \rightarrow \infty$ are described by an infinitely divisible distribution that belongs to the class of stable Cauchy laws. This work generalizes the results for the random $m$-ary recursive trees in Berzunza (2015). Our approach is based on a remarkable decomposition of the size of the giant percolation cluster as a sum of essentially independent random variables which may be useful for studying percolation on other trees with logarithmic height; for instance in this work we study also the case of regular trees.


title: Time-averaged height distribution of the Kardar-Parisi-Zhang interface

abstract: We study the complete probability distribution $\mathcal{P}\left(\bar{H},t\right)$ of the time-averaged height $\bar{H}=(1/t)\int_0^t h(x=0,t')\,dt'$ at point $x=0$ of an evolving 1+1 dimensional Kardar-Parisi-Zhang (KPZ) interface $h\left(x,t\right)$. We focus on short times and flat initial condition and employ the optimal fluctuation method to determine the variance and the third cumulant of the distribution, as well as the asymmetric stretched-exponential tails. The tails scale as $-\ln\mathcal{P}\sim\left|\bar{H}\right|^{3/2} \! /\sqrt{t}$ and $-\ln\mathcal{P}\sim\left|\bar{H}\right|^{5/2} \! /\sqrt{t}$, similarly to the previously determined tails of the one-point KPZ height statistics at specified time $t'=t$. The optimal interface histories, dominating these tails, are markedly different. Remarkably, the optimal history, $h\left(x=0,t\right)$, of the interface height at $x=0$ is a non-monotonic function of time: the maximum (or minimum) interface height is achieved at an intermediate time. We also address a more general problem of determining the probability density of observing a given height history of the KPZ interface at point $x=0$.


title: Dispersionless Multi-Dimensional Integrable Systems and Related Conformal Structure Generating Equations of Mathematical Physics

abstract: Using diffeomorphism group vector fields on $\mathbb{C}$-multiplied tori and the related Lie-algebraic structures, we study multi-dimensional dispersionless integrable systems that describe conformal structure generating equations of mathematical physics. An interesting modification of the devised Lie-algebraic approach subject to spatial-dimensional invariance and meromorphicity of the related differential-geometric structures is described and applied in proving complete integrability of some conformal structure generating equations. As examples, we analyze the Einstein-Weyl metric equation, the modified Einstein-Weyl metric equation, the Dunajski heavenly equation system, the first and second conformal structure generating equations and the inverse first Shabat reduction heavenly equation. We also analyze the modified Pleba\'nski heavenly equations, the Husain heavenly equation and the general Monge equation along with their multi-dimensional generalizations. In addition, we construct superconformal analogs of the Whitham heavenly equation.


title: Matrix-free multigrid solvers for phase-field fracture problems

abstract: In this work, we present a framework for the matrix-free solution to a monolithic quasi-static phase-field fracture model with geometric multigrid methods. Using a standard matrix based approach within the Finite Element Method requires lots of memory, which eventually becomes a serious bottleneck. A matrix-free approach overcomes this problems and greatly reduces the amount of required memory, allowing to solve larger problems on available hardware. One key challenge is concerned with the crack irreversibility for which a primal-dual active set method is employed. Here, the Active-Set values of fine meshes must be available on coarser levels of the multigrid algorithm. The developed multigrid method provides a preconditioner for a generalized minimal residual solver. This method is used for solving the linear equations inside Newton's method for treating the overall nonlinear-monolithic discrete displacement/phase-field formulation. Several numerical examples demonstrate the performance and robustness of our solution technology. Mesh refinement studies, variations in the phase-field regularization parameter, iterations numbers of the linear and nonlinear solvers, and some parallel performances are conducted to substantiate the efficiency of the proposed solver for single fractures, multiple pressurized fractures, and a L-shaped panel test in three dimensions.


title: Global H\"older estimates for 2D linearized Monge-Amp\`ere equations with right-hand side in divergence form

abstract: We establish global H\"older estimates for solutions to inhomogeneous linearized Monge-Amp\`ere equations in two dimensions with the right hand side being the divergence of a bounded vector field. These equations arise in the semi-geostrophic equations in meteorology and in the approximation of convex functionals subject to a convexity constraint using fourth order Abreu type equations. Our estimates hold under natural assumptions on the domain, boundary data and Monge-Amp\`ere measure being bounded away from zero and infinity. They are an affine invariant and degenerate version of global H\"older estimates by Murthy-Stampacchia and Trudinger for second order elliptic equations in divergence form.


title: Disentangling fission fragment evaporation models using TALYS

abstract: Several computer codes based on phenomenological models are being developed with the aim of obtaining fission observables, such as neutron and gamma multiplicities and product yields. Key points in these calculations, which are handled differently by the various codes, are the sharing of the total excitation energy between the fragments and the generation of angular momenta. After the initial states of the fragments is set, the de-excitation through the emission of neutrons and photons can be model. However, many models also include tunable parameters that are used to make the calculations conform with literature data. As a result, it could be possible to obtain good agreement with experimental results even with questionable assumptions on the initial conditions. To test the assumptions made by different fission models the code DE$\ell$FIN has been developed. DE$\ell$FIN starts with the fission fragments as generated by the models and calculates the average neutron emission as a function of mass in a transparent and coherent way using the nuclear reaction code TALYS. Hence, the probabilities of neutron emission in competition with $\gamma$ de-excitation comes from parameters of the TALYS code that have not been optimised. The results are then compared looking at general trends and at the difference to what is obtained with the stand-alone versions of the codes. In this study, the output of DE$\ell$FIN for five of the most commonly used fission codes is compared to the results of the stand-alone versions, and to experimental data, for the reactions $^{235}$U(n$_{\mathrm{th}}$,f), $^{239}$Pu(n$_{\mathrm{th}}$,f) and $^{252}$Cf(sf).


title: The BLAS API of BLASFEO: optimizing performance for small matrices

abstract: BLASFEO is a dense linear algebra library providing high-performance implementations of BLAS- and LAPACK-like routines for use in embedded optimization and other applications targeting relatively small matrices. BLASFEO defines an API which uses a packed matrix format as its native format. This format is analogous to the internal memory buffers of optimized BLAS, but it is exposed to the user and it removes the packing cost from the routine call. For matrices fitting in cache, BLASFEO outperforms optimized BLAS implementations, both open-source and proprietary. This paper investigates the addition of a standard BLAS API to the BLASFEO framework, and proposes an implementation switching between two or more algorithms optimized for different matrix sizes. Thanks to the modular assembly framework in BLASFEO, tailored linear algebra kernels with mixed column- and panel-major arguments are easily developed. This BLAS API has lower performance than the BLASFEO API, but it nonetheless outperforms optimized BLAS and especially LAPACK libraries for matrices fitting in cache. Therefore, it can boost a wide range of applications, where standard BLAS and LAPACK libraries are employed and the matrix size is moderate. In particular, this paper investigates the benefits in scientific programming languages such as Octave, SciPy and Julia.


title: On the image of polynomials evaluated on incidence algebras: a counter-example and a solution

abstract: In this paper, we investigate the subset obtained by evaluations of a fixed multilinear polynomial on a given algebra. We provide an example of a multilinear polynomial, whose image is not a vector subspace; namely, the product of two commutators need not to be a subspace whenever evaluated on certain subalgebras of upper triangular matrices (the so-called incidence algebras). In the last part of the paper, given that the field is infinite, we reduce the problem of the description of the image of a polynomial evaluated on an incidence algebra to the study of evaluations of a certain family of polynomials on its Jacobson radical. In particular, we are able to describe the image of multilinear polynomials evaluated on the algebra of upper triangular matrices.


title: Time versus Hardware: Reducing Qubit Counts with a (Surface Code) Data Bus

abstract: We introduce a data bus, for reducing the qubit counts within quantum computations (protected by surface codes). For general computations, an automated trade-off analysis (software tool and source code are open sourced and available online) is performed to determine to what degree qubit counts are reduced by the data bus: is the time penalty worth the qubit count reductions? We provide two examples where the qubit counts are convincingly reduced: 1) interaction of two surface code patches on NISQ machines with 28 and 68 qubits, and 2) very large-scale circuits with a structure similar to state-of-the-art quantum chemistry circuits. The data bus has the potential to transform all layers of the quantum computing stack (e.g., as envisioned by Google, IBM, Riggeti, Intel), because it simplifies quantum computation layouts, hardware architectures and introduces lower qubits counts at the expense of a reasonable time penalty.


title: Dynamics of weighted composition operators on spaces of continuous functions

abstract: Our study is focused on the dynamics of weighted composition operators defined on a locally convex space $E\hookrightarrow (C(X),\tau_p)$ with $X$ being a topological Hausdorff space containing at least two different points and such that the evaluations $\{\delta_x:\ x\in X\}$ are linearly independent in $E'$. We prove, when $X$ is compact and $E$ is a Banach space containing a nowhere vanishing function, that a weighted composition operator $C_{\varphi,\omega}$ is never weakly supercyclic on $E$. We also prove that if the symbol $\varphi$ lies in the unit ball of $A(\mathbb{D})$, then every weighted composition operator can never be $\tau_p$-supercyclic neither on $C(\mathbb{D})$ nor on the disc algebra $A(\mathbb{D})$. Finally, we obtain Ansari-Bourdon type results and conditions on the spectrum for arbitrary weakly supercyclic operators, and we provide necessary conditions for a composition operator to be weakly supercyclic on the space of holomorphic functions defined in non necessarily simply connected planar domains. As a consequence, we show that no composition operator can be weakly supercyclic neither on the space of holomorphic functions on the punctured disc nor in the punctured plane.


title: First-principle study of the melting temperature of MgO

abstract: Using first-principles only, we calculate the melting point of MgO, also called periclase or magnesia. The random phase approximation (RPA) is used to include the exact exchange as well as local and non-local many-body correlation terms, in order to provide high accuracy. Using the free energy method, we obtain the melting temperature directly from the internal energies calculated with DFT. The free energy differences between the ensembles generated by the molecular dynamics simulations are calulated with thermodynamic integration or thermodynamic perturbation theory. The predicted melting temperature is $T_m^\text{RPA} =3043\pm 86$ K, and the values obtained with the PBE and SCAN functionals are $T_m^\text{PBE} = 2747 \pm 59$ K and $T_m^\text{SCAN} = 3032\pm 53 $ K.


title: First-order perturbation theory for material changes in the surrounding of open optical resonators

abstract: The single-mode approximation of the resonant state expansion has proven to give accurate first-order approximations of resonance shifts and linewidth changes when modifying the material properties inside open optical resonators. Here, we extend this first-order perturbation theory to modifications of the material properties in the surrounding medium. As a side product of our derivations, we retrieve the already known analytical normalization condition for resonant states. We apply our theory to two example systems: A metallic nanosphere and a one-dimensional photonic crystal slab.


title: Time and Energy-Optimal Lane Change Maneuvers for Cooperating Connected Automated Vehicles

abstract: We derive optimal control policies for a Connected and Automated Vehicle (CAV) cooperating with neighboring CAVs to implement a highway lane change maneuver. We optimize the maneuver time and subsequently minimize the associated energy consumption of all cooperating vehicles in this maneuver. We prove structural properties of the optimal policies which simplify the solution derivations and lead to analytical optimal control expressions. The solutions, when they exist, are guaranteed to satisfy safety constraints for all vehicles involved in the maneuver. Simulation results show the effectiveness of the proposed solution and significant performance improvements compared to maneuvers performed by human-driven vehicles.


title: Convergence of fully discrete implicit and semi-implicit approximations of nonlinear parabolic equations

abstract: The article addresses the convergence of implicit and semi-implicit, fully discrete approximations of a class of nonlinear parabolic evolution problems. Such schemes are popular in the numerical solution of evolutions defined with the $p$-Laplace operator since the latter lead to linear systems of equations in the time steps. The semi-implicit treatment of the operator requires introducing a regularization parameter that has to be suitably related to other discretization parameters. To avoid restrictive, unpractical conditions, a careful convergence analysis has to be carried out. The arguments presented in this article show that convergence holds under a moderate condition that relates the step size to the regularization parameter but which is independent of the spatial resolution.


title: Cross-Sensor Periocular Biometrics in a Global Pandemic: Comparative Benchmark and Novel Multialgorithmic Approach

abstract: The massive availability of cameras results in a wide variability of imaging conditions, producing large intra-class variations and a significant performance drop if heterogeneous images are compared for person recognition. However, as biometrics is deployed, it is common to replace damaged or obsolete hardware, or to exchange information between heterogeneous applications. Variations in spectral bands can also occur. For example, surveillance face images (typically acquired in the visible spectrum, VIS) may need to be compared against a legacy iris database (typically acquired in near-infrared, NIR). Here, we propose a multialgorithmic approach to cope with periocular images from different sensors. With face masks in the front line against COVID-19, periocular recognition is regaining popularity since it is the only face region that remains visible. We integrate different comparators with a fusion scheme based on linear logistic regression, in which scores are represented by log-likelihood ratios. This allows easy interpretation of scores and the use of Bayes thresholds for optimal decision-making since scores from different comparators are in the same probabilistic range. We evaluate our approach in the context of the Cross-Eyed Competition, whose aim was to compare recognition approaches when NIR and VIS periocular images are matched. Our approach achieves EER=0.2% and FRR of just 0.47% at FAR=0.01%, representing the best overall approach of the competition. Experiments are also reported with a database of VIS images from different smartphones. We also discuss the impact of template size and computation times, with the most computationally heavy comparator playing an important role in the results. Lastly, the proposed method is shown to outperform other popular fusion approaches, such as the average of scores, SVMs or Random Forest.


title: Report on Tests and Measurements of Hadronic Interaction Properties with Air Showers

abstract: We present a summary of recent tests and measurements of hadronic interaction properties with air showers. This report has a special focus on muon density measurements. Several experiments reported deviations between simulated and recorded muon densities in extensive air showers, while others reported no discrepancies. We combine data from eight leading air shower experiments to cover shower energies from PeV to tens of EeV. Data are combined using the z-scale, a unified reference scale based on simulated air showers. Energy-scales of experiments are cross-calibrated. Above 10 PeV, we find a muon deficit in simulated air showers for each of the six considered hadronic interaction models. The deficit is increasing with shower energy. For the models EPOS-LHC and QGSJet-II.04, the slope is found significant at 8 sigma.


title: Self-induced Magnetic Flux Structure in the Magnetic Superconductor RbEuFe$_4$As$_4$

abstract: We report an unusual enhancement of the magnetic induction in single crystals of the magnetic superconductor RbEuFe$_4$As$_4$ , highlighting the interplay between superconducting and magnetic subsystems in this material. Contrary to the conventional Meissner expulsion of magnetic flux below the superconducting transition temperature, we observe a substantial boost of the magnetic flux density upon approaching the magnetic transition temperature, Tm. Direct imaging of the flux evolution with a magneto-optical technique, shows that the magnetic subsystem serves as an internal magnetic flux pump, drawing Abrikosov vortices from the surface, while the superconducting subsystem controls their conveyance into the bulk of the magnetic superconductor via a peculiar self-organized critical state.


title: Hybrid Direct-Indirect Adaptive Control of Nonlinear System with Unmatched Uncertainty

abstract: In this paper, we present a hybrid direct-indirect model reference adaptive controller (MRAC), to address a class of problems with matched and unmatched uncertainties. In the proposed architecture, the unmatched uncertainty is estimated online through a companion observer model. Upon convergence of the observer, the unmatched uncertainty estimate is remodeled into a state dependent linear form to augment the nominal system dynamics. Meanwhile, a direct adaptive controller designed for a switching system cancels the effect of matched uncertainty in the system and achieves reference model tracking. We demonstrate that the proposed hybrid controller can handle a broad class of nonlinear systems with both matched and unmatched uncertainties


title: Modifying the Chi-square and the CMH test for population genetic inference: adapting to over-dispersion

abstract: Evolve and resequence studies provide a popular approach to simulate evolution in the lab and explore its genetic basis. In this context, the chi-square test, Fishers exact test, as well as the Cochran-Mantel-Haenszel test are commonly used to infer genomic positions affected by selection from temporal changes in allele frequency. However, the null model associated with these tests does not match the null hypothesis of actual interest. Indeed due to genetic drift and possibly other additional noise components such as pool sequencing, the null variance in the data can be substantially larger than accounted forby these common test statistics. This leads to p-values that are systematically too small and therefore a huge number of false positive results. Even, if the ranking rather than the actual p-values is of interest, a naive application of the mentioned tests will give misleading results, as the amount of over-dispersion varies from locus to locus. We therefore propose adjusted statistics that take the over-dispersion into account while keeping the formulas simple. This is particularly useful in genome-wide applications, where millions of SNPs can be handled with little computational effort. We then apply the adapted test statistics to real data fromDrosophila, and investigate how in-formation from intermediate generations can be included when avail-able. The obtained formulas may also be useful in other situations, provided that the null variance either is known or can be estimated.


title: Boundary-weighted Domain Adaptive Neural Network for Prostate MR Image Segmentation

abstract: Accurate segmentation of the prostate from magnetic resonance (MR) images provides useful information for prostate cancer diagnosis and treatment. However, automated prostate segmentation from 3D MR images still faces several challenges. For instance, a lack of clear edge between the prostate and other anatomical structures makes it challenging to accurately extract the boundaries. The complex background texture and large variation in size, shape and intensity distribution of the prostate itself make segmentation even further complicated. With deep learning, especially convolutional neural networks (CNNs), emerging as commonly used methods for medical image segmentation, the difficulty in obtaining large number of annotated medical images for training CNNs has become much more pronounced that ever before. Since large-scale dataset is one of the critical components for the success of deep learning, lack of sufficient training data makes it difficult to fully train complex CNNs. To tackle the above challenges, in this paper, we propose a boundary-weighted domain adaptive neural network (BOWDA-Net). To make the network more sensitive to the boundaries during segmentation, a boundary-weighted segmentation loss (BWL) is proposed. Furthermore, an advanced boundary-weighted transfer leaning approach is introduced to address the problem of small medical imaging datasets. We evaluate our proposed model on the publicly available MICCAI 2012 Prostate MR Image Segmentation (PROMISE12) challenge dataset. Our experimental results demonstrate that the proposed model is more sensitive to boundary information and outperformed other state-of-the-art methods.


title: A Mean Field Theory of Batch Normalization

abstract: We develop a mean field theory for batch normalization in fully-connected feedforward neural networks. In so doing, we provide a precise characterization of signal propagation and gradient backpropagation in wide batch-normalized networks at initialization. Our theory shows that gradient signals grow exponentially in depth and that these exploding gradients cannot be eliminated by tuning the initial weight variances or by adjusting the nonlinear activation function. Indeed, batch normalization itself is the cause of gradient explosion. As a result, vanilla batch-normalized networks without skip connections are not trainable at large depths for common initialization schemes, a prediction that we verify with a variety of empirical simulations. While gradient explosion cannot be eliminated, it can be reduced by tuning the network close to the linear regime, which improves the trainability of deep batch-normalized networks without residual connections. Finally, we investigate the learning dynamics of batch-normalized networks and observe that after a single step of optimization the networks achieve a relatively stable equilibrium in which gradients have dramatically smaller dynamic range. Our theory leverages Laplace, Fourier, and Gegenbauer transforms and we derive new identities that may be of independent interest.


title: Hyperbolicity of asymmetric lemon billiards

abstract: Asymmetric lemon billiards was introduced in [CMZZ], where the billiard table $Q(r,b,R)$ is the intersection of two round disks with radii $r\le R$, respectively, and $b$ measures the distance between the two centers. It is conjectured [BZZ] that the asymmetric lemon billiards is hyperbolic when the arc $\Gamma_r$ is a major arc and $R$ is large. In this paper we prove this conjecture for sufficiently large $R$.


title: Two-dimensional electronic transport in rubrene: the impact of inter-chain coupling

abstract: Organic semi-conductors have unique electronic properties and are important systems both at the fundamental level and also for their applications in electronic devices. In this article we focus on the particular case of rubrene which has one of the best electronic transport properties for application purposes. We show that this system can be well simulated by simple tight-binding systems representing one-dimensional (1D) chains that are weakly coupled to their neighboring chains in the same plane. This makes in principle this rubrene system somehow intermediate between 1D and isotropic 2D models. We analyse in detail the dc-transport and terahertz conductivity in the 1D and in the anisotropic 2D models. The transient localisation scenario allows us to reproduce satisfactorily some basics results such as mobility anisotropy and orders of magnitude as well as ac-conductivity in the terahertz range. This model shows in particular that even a weak inter-chain coupling is able to improve notably the propagation along the chains. This suggest also that a strong inter-chain coupling is important to get organic semi-conductors with the best possible transport properties for applicative purposes.


title: Economic MPC using a Cyclic Horizon with Application to Networked Control Systems

abstract: In this paper, we analyze an economic model predictive control scheme with terminal region and cost, where the system is optimally operated in a certain subset of the state space. The predictive controller operates with a cyclic horizon, which means that starting from an initial length, the horizon is reduced by one at each time step before it is restored to its maximum length again after one cycle. We give performance guarantees for the closed loop, and under a suitable dissipativity condition, establish convergence to the optimal subset. Moreover, we present conditions under which asymptotic stability of the optimal subset can be guaranteed. The results are illustrated in a practical example from the context of Networked Control Systems, which initially motivated the development of the theory presented in this paper.


title: Maximising the Number of Cycles in Graphs with Forbidden Subgraphs

abstract: Fix $k \ge 2$ and let $H$ be a graph with $\chi(H) = k+1$ containing a critical edge. We show that for sufficiently large $n$, the unique $n$-vertex $H$-free graph containing the maximum number of cycles is $T_k(n)$. This resolves both a question and a conjecture of Arman, Gunderson and Tsaturian.


title: Domain Partitioning Network

abstract: Standard adversarial training involves two agents, namely a generator and a discriminator, playing a mini-max game. However, even if the players converge to an equilibrium, the generator may only recover a part of the target data distribution, in a situation commonly referred to as mode collapse. In this work, we present the Domain Partitioning Network (DoPaNet), a new approach to deal with mode collapse in generative adversarial learning. We employ multiple discriminators, each encouraging the generator to cover a different part of the target distribution. To ensure these parts do not overlap and collapse into the same mode, we add a classifier as a third agent in the game. The classifier decides which discriminator the generator is trained against for each sample. Through experiments on toy examples and real images, we show the merits of DoPaNet in covering the real distribution and its superiority with respect to the competing methods. Besides, we also show that we can control the modes from which samples are generated using DoPaNet.


title: Coloring squares of graphs with mad constraints

abstract: A proper vertex $k$-coloring of a graph $G=(V,E)$ is an assignment $c:V\to \{1,2,\ldots,k\}$ of colors to the vertices of the graph such that no two adjacent vertices are associated with the same color. The square $G^2$ of a graph $G$ is the graph defined by $V(G)=V(G^2)$ and $uv \in E(G^2)$ if and only if the distance between $u$ and $v$ is at most two. We denote by $\chi(G^2)$ the chromatic number of $G^2$, which is the least integer $k$ such that a $k$-coloring of $G^2$ exists. By definition, at least $\Delta(G)+1$ colors are needed for this goal, where $\Delta(G)$ denotes the maximum degree of the graph $G$. In this paper, we prove that the square of every graph $G$ with $\text{mad}(G)<4$ and $\Delta(G) \geqslant 8$ is $(3\Delta(G)+1)$-choosable and even correspondence-colorable. Furthermore, we show a family of $2$-degenerate graphs $G$ with $\text{mad}(G)<4$, arbitrarily large maximum degree, and $\chi(G^2)\geqslant \frac{5\Delta(G)}{2}$, improving the result of Kim and Park.


title: Exclusive vector meson production in heavy ion collisions

abstract: We discuss the salient features of exclusive vector meson production in heavy ion collisions at LHC energies. Special attention is paid to the space-time picture of the process. We account for both coherent and incoherent contributions. The explicit quantitative predictions are given for the $\rho$-meson differential cross section in lead-lead collisions in different kinematical configurations relevant for the LHCb and ALICE experiments.


title: Deep CNN-based Speech Balloon Detection and Segmentation for Comic Books

abstract: We develop a method for the automated detection and segmentation of speech balloons in comic books, including their carrier and tails. Our method is based on a deep convolutional neural network that was trained on annotated pages of the Graphic Narrative Corpus. More precisely, we are using a fully convolutional network approach inspired by the U-Net architecture, combined with a VGG-16 based encoder. The trained model delivers state-of-the-art performance with an F1-score of over 0.94. Qualitative results suggest that wiggly tails, curved corners, and even illusory contours do not pose a major problem. Furthermore, the model has learned to distinguish speech balloons from captions. We compare our model to earlier results and discuss some possible applications.


title: A Nonparametric Multi-view Model for Estimating Cell Type-Specific Gene Regulatory Networks

abstract: We present a Bayesian hierarchical multi-view mixture model termed Symphony that simultaneously learns clusters of cells representing cell types and their underlying gene regulatory networks by integrating data from two views: single-cell gene expression data and paired epigenetic data, which is informative of gene-gene interactions. This model improves interpretation of clusters as cell types with similar expression patterns as well as regulatory networks driving expression, by explaining gene-gene covariances with the biological machinery regulating gene expression. We show the theoretical advantages of the multi-view learning approach and present a Variational EM inference procedure. We demonstrate superior performance on both synthetic data and real genomic data with subtypes of peripheral blood cells compared to other methods.

categories: stat.ML cs.LG q-bio.GN q-bio.MN q-bio.QM


title: Unitary time-evolution in stochastic time-dependent Hilbert spaces

abstract: In this work we study the unitary time-evolutions of quantum systems defined on infinite-dimensional separable time-dependent Hilbert spaces. Two possible cases are considered: a quantum system defined on a stochastic interval and another one defined on a Hilbert space with stochastic integration measure (stochastic time-dependent scalar product). The formulations of the two problems and a comparison with the general theory of open quantum systems are discussed. Possible physical applications of the situations considered are analyzed.


title: Relative cubulations and groups with a 2-sphere boundary

abstract: We introduce a new kind of action of a relatively hyperbolic group on a CAT(0) cube complex, called a relatively geometric action. We provide an application to characterize finite-volume Kleinian groups in terms of action on cube complexes, analogous to the results of Markovic and Ha\"issinsky in the closed case.


title: The reflection principle in the control problem of the heat equation

abstract: We consider the control problem for the generalized heat equation for a Schroedinger operator on a domain with a reflection symmetry with respect to a hyperplane. We show that if this system is null-controllable, then so is the system on its respective parts and the corresponding control cost does not exceed the one on the whole domain. As an application, we obtain null-controllability results for the heat equation on half-spaces, orthants, and sectors of angle $\pi/2$. As a byproduct, we also obtain explicit control cost bounds for the heat equation on certain triangles and corresponding prisms in terms of geometric parameters of the control set.


title: Evaluating the Search Phase of Neural Architecture Search

abstract: Neural Architecture Search (NAS) aims to facilitate the design of deep networks for new tasks. Existing techniques rely on two stages: searching over the architecture space and validating the best architecture. NAS algorithms are currently compared solely based on their results on the downstream task. While intuitive, this fails to explicitly evaluate the effectiveness of their search strategies. In this paper, we propose to evaluate the NAS search phase. To this end, we compare the quality of the solutions obtained by NAS search policies with that of random architecture selection. We find that: (i) On average, the state-of-the-art NAS algorithms perform similarly to the random policy; (ii) the widely-used weight sharing strategy degrades the ranking of the NAS candidates to the point of not reflecting their true performance, thus reducing the effectiveness of the search process. We believe that our evaluation framework will be key to designing NAS strategies that consistently discover architectures superior to random ones.


title: Effective-field theory analysis of the $\tau^{-}\rightarrow (K \pi)^{-}\nu_{\tau}$ decays

abstract: We analyze the $\tau^-\to(K\pi)^-\nu_\tau$ decays within an effective field theory description of heavy new physics (NP) modifying the SM left-handed weak charged current and include refined SM input for the participant meson form factors exploiting chiral symmetry, dispersion relations and (lattice) data. We include the leading dimension six operators and work at linear order in the effective couplings. Within this setting we: i) follow the derivation in Phys. Rev. Lett. 120 (2018) no.14, 141803 where it was proved unambiguously that it is impossible to understand the BaBar anomaly in the CP asymmetry measurement within this framework. We allow for reasonable variations of the hadronic input involved and study the associated uncertainty; ii) first show that the anomalous bump present in the published Belle data for the $K_S\pi^-$ invariant mass distribution close to threshold cannot be due to heavy NP; iii) first bind the heavy NP effective couplings using $\tau^-\to(K\pi)^-\nu_\tau$ decays and show that they are competitive with those found in hyperon semileptonic decays (but clearly not with those obtained in Kaon (semi)leptonic decays for NP scalar currents). We put forward that the comparison of the considered tau decays with (semi)leptonic kaon and hyperon decays provide with meaningful tests of lepton universality for (NP) tensor interactions. We also compare the SM predictions with the possible deviations caused by NP in a couple of Dalitz plot distributions, in the forward-backward asymmetry and in the di-meson invariant mass distribution and discuss the most interesting measurements to be performed at Belle-II using these decays data.


title: Integrally Closed Residuated Lattices

abstract: A residuated lattice is defined to be integrally closed if it satisfies the equations x\x = e and x/x = e. Every integral, cancellative, or divisible residuated lattice is integrally closed, and, conversely, every bounded integrally closed residuated lattice is integral. It is proved that the mapping a -> (a\e)\e on any integrally closed residuated lattice is a homomorphism onto a lattice-ordered group. A Glivenko-style property is then established for varieties of integrally closed residuated lattices with respect to varieties of lattice-ordered groups, showing in particular that integrally closed residuated lattices form the largest variety of residuated lattices admitting this property with respect to lattice-ordered groups. The Glivenko property is used to obtain a sequent calculus admitting cut-elimination for the variety of integrally closed residuated lattices and to establish the decidability, indeed PSPACE-completenes, of its equational theory. Finally, these results are related to previous work on (pseudo) BCI-algebras, semi-integral residuated partially ordered monoids, and algebras for Casari's comparative logic.


title: Total Variation and Mean Curvature PDEs on $\mathbb{R}^d \rtimes S^{d-1}$

abstract: Total variation regularization and total variation flows (TVF) have been widely applied for image enhancement and denoising. To include a generic preservation of crossing curvilinear structures in TVF we lift images to the homogeneous space $M = \mathbb{R}^d \rtimes S^{d-1}$ of positions and orientations as a Lie group quotient in SE(d). For d = 2 this is called 'total roto-translation variation' by Chambolle & Pock. We extend this to d = 3, by a PDE-approach with a limiting procedure for which we prove convergence. We also include a Mean Curvature Flow (MCF) in our PDE model on M. This was first proposed for d = 2 by Citti et al. and we extend this to d = 3. Furthermore, for d = 2 we take advantage of locally optimal differential frames in invertible orientation scores (OS). We apply our TVF and MCF in the denoising/enhancement of crossing fiber bundles in DW-MRI. In comparison to data-driven diffusions, we see a better preservation of bundle boundaries and angular sharpness in fiber orientation densities at crossings. We support this by error comparisons on a noisy DW-MRI phantom. We also apply our TVF and MCF in enhancement of crossing elongated structures in 2D images via OS, and compare the results to nonlinear diffusions (CED-OS) via OS.


title: Limiting magnetic field for minimal deformation of a magnetised neutron star

abstract: In this work we study the structure of neutron stars under the effect of a poloidal magnetic field and determine the limiting highest magnetic field intensity which still allows a satisfactory description of magnetic neutron stars in the spherical symmetry regime. We describe different compositions of stars (nucleonic, hyperonic, and hybrid), using three state-of-the-art relativistic mean field models for the microscopic description of matter, which are in agreement with experimental and observational data. The structure of stars is described by the general relativistic solution of both Einstein's field equations assuming a spherical symmetry, and Einstein-Maxwell's field equations assuming an axi-symmetric deformation. We find a limiting magnetic moment of the order of $2\times 10^{31}$Am$^2$, which corresponds to magnetic fields of the order of 10$^{16}$ G at the surface, and $ \sim 10^{17}$ G at the centre of the star, above which the deformation due to the magnetic field is not negligible. We show that the intensity of the magnetic field developed in the star depends on the EoS, and, for a given baryonic mass and fixed magnetic moment, larger fields are attained with softer EoS. We also show that the appearance of exotic degrees of freedom, such as hyperons or a quark core, is disfavored in the presence of a very strong magnetic field. As a consequence, a highly magnetized nucleonic star may suffer an internal conversion due to the decay of the magnetic field, which could be accompanied by a sudden cooling of the star or a gamma ray burst.


title: Visualization of Local Conductance in MoS2/WSe2 Heterostructure Transistors

abstract: The vertical stacking of van der Waals (vdW) materials introduces a new degree of freedom to the research of two-dimensional (2D) systems. The interlayer coupling strongly influences the band structure of the heterostructures, resulting in novel properties that can be utilized for electronic and optoelectronic applications. Based on microwave microscopy studies, we report quantitative electrical imaging on gated molybdenum disulfide (MoS2)/tungsten diselenide (WSe2) heterostructure devices, which exhibit an intriguing anti-ambipolar effect in the transfer characteristics. Interestingly, in the region with significant source-drain current, electrons in the n-type MoS2 and holes in the p-type WSe2 segments are nearly balanced, whereas the heterostructure area is depleted of mobile charges. The configuration is analogous to the p-i-n diode, where the injected carriers dominate in the recombination current. The spatial evolution of local conductance can be ascribed to the lateral band bending and formation of depletion regions along the line of MoS2-heterostructure-WSe2. Our work vividly demonstrates the microscopic origin of novel transport behaviors, which is important for the vibrant field of vdW heterojunction research.


title: The strong running coupling from the gauge sector of Domain Wall lattice QCD with physical quark masses

abstract: We report on the first computation of the strong running coupling at the physical point (physical pion mass) from the ghost-gluon vertex, computed from lattice simulations with three flavors of Domain Wall fermions. We find $\alpha_{\overline{\rm MS}}(m_Z^2)=0.1172(11)$, in remarkably good agreement with the world-wide average. Our computational bridge to this value is the Taylor-scheme strong coupling, which has been revealed of great interest by itself because it can be directly related to the quark-gluon interaction kernel in continuum approaches to the QCD bound-state problem.


title: Aperiodic Weighted Automata and Weighted First-Order Logic

abstract: By fundamental results of Sch\"utzenberger, McNaughton and Papert from the 1970s, the classes of first-order definable and aperiodic languages coincide. Here, we extend this equivalence to a quantitative setting. For this, weighted automata form a general and widely studied model. We define a suitable notion of a weighted first-order logic. Then we show that this weighted first-order logic and aperiodic polynomially ambiguous weighted automata have the same expressive power. Moreover, we obtain such equivalence results for suitable weighted sublogics and finitely ambiguous or unambiguous aperiodic weighted automata. Our results hold for general weight structures, including all semirings, average computations of costs, bounded lattices, and others.


title: From Stress Chains to Acoustic Emission

abstract: A numerical scheme using the combined finite-discrete element methods (FDEM) is employed to study a model of an earthquake system comprising a granular layer embedded in a formation. When the formation is driven so as to shear the granular layer, a system of stress chains emerges. The stress chains endow the layer with resistance to shear and on failure launch broadcasts into the formation. These broadcasts, received as acoustic emission, provide a remote monitor of the state of the granular layer, of the earthquake system.


title: BFKL Pomeron and the survival factor

abstract: We consider the absorptive corrections and the rapidity gap survival factor which are necessary to provide the unitarization of the BFKL Pomeron. In particular we discuss the role of the enhanced screening diagrams.


title: Deep Long Asymmetric Occultation in EPIC 204376071

abstract: We have discovered a young M star of mass $0.16\,M_\odot$ and radius $0.63\,R_\odot$, likely in the Upper Sco Association, that exhibits only a single $80\%$ deep occultation of 1-day duration. The star has frequent flares and a low-amplitude rotational modulation, but is otherwise quiet over 160 days of cumulative observation during K2 Campaigns C2 and C15. We discuss how such a deep eclipse is not possible by one star crossing another in any binary or higher-order stellar system in which no mass transfer has occurred. The two possible explanations we are left with are (1) orbiting dust or small particles (e.g., a disk bound to a smaller orbiting body, or unbound dust that emanates from such a body); or (2) a transient accretion event of dusty material near the corotation radius of the star. In either case, the time between such occultation events must be longer than $\sim$80 days. We model a possible orbiting occulter both as a uniform elliptically shaped surface (e.g., an inclined circular disk) and as a `dust sheet' with a gradient of optical depth behind its leading edge. The required masses in such dust features are then $\gtrsim 3 \times 10^{19}$ g and $\gtrsim 10^{19}$ g, for the two cases, respectively.


title: Learned Step Size Quantization

abstract: Deep networks run with low precision operations at inference time offer power and space advantages over high precision alternatives, but need to overcome the challenge of maintaining high accuracy as precision decreases. Here, we present a method for training such networks, Learned Step Size Quantization, that achieves the highest accuracy to date on the ImageNet dataset when using models, from a variety of architectures, with weights and activations quantized to 2-, 3- or 4-bits of precision, and that can train 3-bit models that reach full precision baseline accuracy. Our approach builds upon existing methods for learning weights in quantized networks by improving how the quantizer itself is configured. Specifically, we introduce a novel means to estimate and scale the task loss gradient at each weight and activation layer's quantizer step size, such that it can be learned in conjunction with other network parameters. This approach works using different levels of precision as needed for a given system and requires only a simple modification of existing training code.


title: Quantum Package 2.0: An Open-Source Determinant-Driven Suite of Programs

abstract: \textsc{Quantum Package} is an open-source programming environment for quantum chemistry specially designed for wave function methods. Its main goal is the development of determinant-driven selected configuration interaction (sCI) methods and multi-reference second-order perturbation theory (PT2). The determinant-driven framework allows the programmer to include any arbitrary set of determinants in the reference space, hence providing greater methodological freedoms. The sCI method implemented in \textsc{Quantum Package} is based on the CIPSI (Configuration Interaction using a Perturbative Selection made Iteratively) algorithm which complements the variational sCI energy with a PT2 correction. Additional external plugins have been recently added to perform calculations with multireference coupled cluster theory and range-separated density-functional theory. All the programs are developed with the IRPF90 code generator, which simplifies collaborative work and the development of new features. \textsc{Quantum Package} strives to allow easy implementation and experimentation of new methods, while making parallel computation as simple and efficient as possible on modern supercomputer architectures. Currently, the code enables, routinely, to realize runs on roughly 2\,000 CPU cores, with tens of millions of determinants in the reference space. Moreover, we have been able to push up to 12\,288 cores in order to test its parallel efficiency. In the present manuscript, we also introduce some key new developments: i) a renormalized second-order perturbative correction for efficient extrapolation to the full CI limit, and ii) a stochastic version of the CIPSI selection performed simultaneously to the PT2 calculation at no extra cost.

categories: physics.comp-ph cond-mat.str-el physics.atom-ph physics.chem-ph


title: The Schinzel Hypothesis for Polynomials

abstract: The Schinzel hypothesis is a famous conjectural statement about primes in value sets of polynomials, which generalizes the Dirichlet theorem about primes in an arithmetic progression. We consider the situation that the ring of integers is replaced by a polynomial ring and prove the Schinzel hypothesis for a wide class of them: polynomials in at least one variable over the integers, polynomials in several variables over an arbitrary field, etc. We achieve this goal by developing a version over rings of the Hilbert specialization property. A polynomial Goldbach conjecture is deduced, along with a result on spectra of rational functions.


title: Optical analysis of a CMB cosmic string candidate

abstract: The complexity of the cosmological scenario regarding cosmic strings (CSs) stands still in the way of a complete understanding. We describe here a promising strategy for the possible detection of these elusive physical entities. It is based on the search of strong gravitational lensing events in the location area of the CS candidate (CSc-1), which was declared in a previous work by CMB analysis. Using photometric and geometric criteria, we identified pairs of candidates of lensed galaxies (LGCs) in the "string field" (SF), which were then compared with the average density of background galaxy pairs in a set of "control fields" (CFs). We found an excess of $22\%$ (per sq. deg.) of the LGCs in SF, which exceeds the estimated cosmic dispersion. We also found that the number of LGCs is in excess of $29.2\%$ in the angular separation bin $[8'', 9'']$. Finally, we analysed the possibility of a preferred orientation of the line connecting the centres of the LGCs. The orientation is statistically significant for an angular separation bin $[4'',6'']$. Therefore, we found two "windows" for the preferred angular separation for LGCs along the possible CS. However, the confirmation of the gravitational lensing origin of our LGCs requires spectroscopic observations which seem to be justified by the present results. We plan to acquire their spectra as well as to continue the study of the spectral and morphological features of the LGCs in the CSc-1 field and to analyse the other CS-candidates using the same strategy.


title: Assembly of 2N entangled fermions into multipartite composite bosons

abstract: An even number of fermions can behave in a bosonic way. The simplest scenario involves two fermions which can form a single boson. But four fermions can either behave as two bipartite bosons or further assemble into a single four-partite bosonic molecule. In general, for 2N fermions there are many possible arrangements into composite bosons. The question is: what determines which fermionic arrangement is going to be realized in a given situation and can such arrangement be considered truly bosonic? This work aims to find the answer to the above question. We propose an entanglement-based method to assess bosonic quality of fermionic arrangements and apply it to study how the ground state of the extended one-dimensional Hubbard model changes as the strength of intra-particle interactions increases.