title: ANR: Articulated Neural Rendering for Virtual Avatars

abstract: The combination of traditional rendering with neural networks in Deferred Neural Rendering (DNR) provides a compelling balance between computational complexity and realism of the resulting images. Using skinned meshes for rendering articulating objects is a natural extension for the DNR framework and would open it up to a plethora of applications. However, in this case the neural shading step must account for deformations that are possibly not captured in the mesh, as well as alignment inaccuracies and dynamics -- which can confound the DNR pipeline. We present Articulated Neural Rendering (ANR), a novel framework based on DNR which explicitly addresses its limitations for virtual human avatars. We show the superiority of ANR not only with respect to DNR but also with methods specialized for avatar creation and animation. In two user studies, we observe a clear preference for our avatar model and we demonstrate state-of-the-art performance on quantitative evaluation metrics. Perceptually, we observe better temporal stability, level of detail and plausibility.


title: New plans orthogonal through the block factor

abstract: In the present paper we construct plans orthogonal through the block factor (POTBs). We describe procedures for adding blocks as well as factors to an initial plan and thus generate a bigger plan. Using these procedures we construct POTBs for symmetrical experiments with factors having three or more levels. We also construct a series of plans inter-class orthogonal through the block factor for two-level factors.


title: Analytic critical points of charged Renyi entropies from hyperbolic black holes

abstract: We analytically study phase transitions of holographic charged Renyi entropies in two gravitational systems dual to the $\mathcal{N}=4$ super-Yang-Mills theory at finite density and zero temperature. The first system is the Reissner-Nordstrom-AdS$_5$ black hole, which has finite entropy at zero temperature. The second system is a charged dilatonic black hole in AdS$_5$, which has zero entropy at zero temperature. Hyperbolic black holes are employed to calculate the Renyi entropies with the entangling surface being a sphere. We perturb each system by a charged scalar field, and look for a zero mode signaling the instability of the extremal hyperbolic black hole. Zero modes as well as the leading order of the full retarded Green's function are analytically solved for both systems, in contrast to previous studies in which only the IR (near horizon) instability was analytically treated.


title: The Future of High-Energy Astrophysical Neutrino Flavor Measurements

abstract: We critically examine the ability of future neutrino telescopes, including Baikal-GVD, KM3NeT, P-ONE, TAMBO, and IceCube-Gen2, to determine the flavor composition of high-energy astrophysical neutrinos, ie, the relative number of $\nu_e$, $\nu_\mu$, and $\nu_\tau$, in light of improving measurements of the neutrino mixing parameters. Starting in 2020, we show how measurements by JUNO, DUNE, and Hyper-Kamiokande will affect our ability to determine the regions of flavor composition at Earth that are allowed by neutrino oscillations under different assumptions of the flavor composition that is emitted by the astrophysical sources. From 2020 to 2040, the error on inferring the flavor composition at the source will improve from $> 40\%$ to less than $6\%$. By 2040, under the assumption that pion decay is the principal production mechanism of high-energy astrophysical neutrinos, a sub-dominant mechanism could be constrained to contribute less than 20\% of the flux at 99.7\% credibility. These conclusions are robust in the nonstandard scenario where neutrino mixing is non-unitary, a scenario that is the target of next-generation experiments, in particular the IceCube-Upgrade. Finally, to illustrate the improvement in using flavor composition to test beyond-the-Standard-Model physics, we examine the possibility of neutrino decay and find that, by 2040, combined neutrino telescope measurements will be able to limit the decay rate of the heavier neutrinos to below $1.8\times 10^{-5} (m/\mathrm{eV})$~s$^{-1}$, at 95\% credibility.


title: Quark cluster expansion model for interpreting finite-T lattice QCD thermodynamics

abstract: We present a unified approach to the thermodynamics of hadron-quark-gluon matter at finite temperatures on the basis of a quark cluster expansion in the form of a generalized Beth-Uhlenbeck approach with a generic ansatz for the hadronic phase shifts that fulfills the Levinson theorem. The change in the composition of the system from a hadron resonance gas to a quark-gluon plasma takes place in the narrow temperature interval of $150 - 185$ MeV where the Mott dissociation of hadrons is triggered by the dropping quark mass as a result of the restoration of chiral symmetry. The deconfinement of quark and gluon degrees of freedom is regulated by the Polyakov loop variable that signals the breaking of the $Z(3)$ center symmetry of the color $SU(3)$ group of QCD. We suggest a Polyakov-loop quark-gluon plasma model with $\mathcal{O}(\alpha_s)$ virial correction and solve the stationarity condition of the thermodynamic potential (gap equation) for the Polyakov loop. The resulting pressure is in excellent agreement with lattice QCD simulations up to high temperatures.


title: A Modern Analysis of Hutchinson's Trace Estimator

abstract: The paper establishes the new state-of-art in the accuracy analysis of Hutchinson's trace estimator. Leveraging tools that have not been previously used in this context, particularly hypercontractive inequalities and concentration properties of sub-gamma distributions, we offer an elegant and modular analysis, as well as numerically superior bounds. Besides these improvements, this work aims to better popularize the aforementioned techniques within the CS community.


title: How Does a Neural Network's Architecture Impact Its Robustness to Noisy Labels?

abstract: Noisy labels are inevitable in large real-world datasets. In this work, we explore an area understudied by previous works -- how the network's architecture impacts its robustness to noisy labels. We provide a formal framework connecting the robustness of a network to the alignments between its architecture and target/noise functions. Our framework measures a network's robustness via the predictive power in its representations -- the test performance of a linear model trained on the learned representations using a small set of clean labels. We hypothesize that a network is more robust to noisy labels if its architecture is more aligned with the target function than the noise. To support our hypothesis, we provide both theoretical and empirical evidence across various neural network architectures and different domains. We also find that when the network is well-aligned with the target function, its predictive power in representations could improve upon state-of-the-art (SOTA) noisy-label-training methods in terms of test accuracy and even outperform sophisticated methods that use clean labels.


title: On Polynomial Representations of the DP Color Function: Theta Graphs and Their Generalizations

abstract: DP-coloring (also called correspondence coloring) is a generalization of list coloring that has been widely studied in recent years after its introduction by Dvo\v{r}\'{a}k and Postle in 2015. As the analogue of the chromatic polynomial $P(G,m)$, the DP color function of a graph $G$, denoted $P_{DP}(G,m)$, counts the minimum number of DP-colorings over all possible $m$-fold covers. It is known that, unlike the list color function $P_{\ell}(G,m)$, for any $g \geq 3$ there exists a graph $G$ with girth $g$ such that $P_{DP}(G,m) < P(G,m)$ when $m$ is sufficiently large. Thus, two fundamental open questions regarding the DP color function are: (i) for which $G$ does there exist an $N \in \mathbb{N}$ such that $P_{DP}(G,m) = P(G,m)$ whenever $m \geq N$, (ii) Given a graph $G$ does there always exist an $N \in \mathbb{N}$ and a polynomial $p(m)$ such that $P_{DP}(G,m) = p(m)$ whenever $m \geq N$? In this paper we give exact formulas for the DP color function of a Theta graph based on the parity of its path lengths. This gives an explicit answer, including the formulas for the polynomials that are not the chromatic polynomial, to both the questions above for Theta graphs. We extend this result to Generalized Theta graphs by characterizing the exact parity condition that ensures the DP color function eventually equals the chromatic polynomial. To answer the second question for Generalized Theta graphs, we confirm it for the larger class of graphs with a feedback vertex set of size one.


title: Forcing and anti-forcing polynomials of a polyomino graph

abstract: The forcing number of a perfect matching $M$ in a graph $G$ is the smallest number of edges inside $M$ that can not be contained in other perfect matchings. The anti-forcing number of $M$ is the smallest number of edges outside $M$ whose removal results in a subgraph with a single perfect matching, that is $M$. Recently, in order to investigate the distributions of forcing numbers and anti-forcing numbers, the forcing polynomial and anti-forcing polynomial were proposed, respectively. In this work, the forcing and anti-forcing polynomials of a polyomino graph are obtained. As consequences, the forcing and anti-forcing spectra of this polyomino graph are determined, and the asymptotic behaviors on the degree of freedom and the sum of all anti-forcing numbers are revealed, respectively.


title: Learning by Self-Explanation, with Application to Neural Architecture Search

abstract: Learning by self-explanation is an effective learning technique in human learning, where students explain a learned topic to themselves for deepening their understanding of this topic. It is interesting to investigate whether this explanation-driven learning methodology broadly used by humans is helpful for improving machine learning as well. Based on this inspiration, we propose a novel machine learning method called learning by self-explanation (LeaSE). In our approach, an explainer model improves its learning ability by trying to clearly explain to an audience model regarding how a prediction outcome is made. LeaSE is formulated as a four-level optimization problem involving a sequence of four learning stages which are conducted end-to-end in a unified framework: 1) explainer learns; 2) explainer explains; 3) audience learns; 4) explainer re-learns based on the performance of the audience. We develop an efficient algorithm to solve the LeaSE problem. We apply LeaSE for neural architecture search on CIFAR-100, CIFAR-10, and ImageNet. Experimental results strongly demonstrate the effectiveness of our method.


title: Lithium pollution of a white dwarf records the accretion of an extrasolar planetesimal

abstract: Tidal disruption and subsequent accretion of planetesimals by white dwarfs can reveal the elemental abundances of rocky bodies in exoplanetary systems. Those abundances provide information on the composition of the nebula from which the systems formed, analogous to how meteorite abundances inform our understanding of the early Solar System. We report the detection of Li, Na, K and Ca in the atmosphere of the white dwarf Gaia DR2 4353607450860305024, which we ascribe to accretion of a planetesimal. Using model atmospheres, we determine abundance ratios of these elements, and with the exception of Li, they are consistent with meteoritic values in the Solar System. We compare the measured Li abundance to measurements in old stars and to expectations from Big Bang nucleosynthesis.


title: Lattice gauge equivariant convolutional neural networks

abstract: We propose Lattice gauge equivariant Convolutional Neural Networks (L-CNNs) for generic machine learning applications on lattice gauge theoretical problems. At the heart of this network structure is a novel convolutional layer that preserves gauge equivariance while forming arbitrarily shaped Wilson loops in successive bilinear layers. Together with topological information, for example from Polyakov loops, such a network can in principle approximate any gauge covariant function on the lattice. We demonstrate that L-CNNs can learn and generalize gauge invariant quantities that traditional convolutional neural networks are incapable of finding.


title: BlackHawk: A tool for computing Black Hole evaporation

abstract: We present the public code \texttt{BlackHawk}, a powerful tool to compute the evaporation of any distribution of Schwarzschild and Kerr Black Holes and the emission spectra of Hawking radiation.


title: Decoherence from General Relativity

abstract: It is of great interest to explore matter in nontrivial quantum arrangements, including Schrodinger cat-like states. Such states are sensitive to decoherence from their environment. Recently, in Ref. [1] we computed the rate of decoherence of a piece of superposed matter that primarily only interacts gravitationally, a dark-matter-Schrodinger-cat-state (DMSCS), within the nonrelativistic approximation. In this work we improve this to a general relativistic analysis. We firstly derive a single particle relativistic Schrodinger equation for a probe particle that passes through the DMSCS; the interaction is provided by the weak field metric of general relativity from the source. For a static DMSCS we find a neat generalization of our previous results. We then turn to the interesting new case of a time dependent DMSCS, which can be provided by a coherently oscillating axion field leading to superposed time dependent oscillations in the metric; a truly quantum-general relativistic phenomenon. We use scattering theory to derive the decoherence rate in all these cases. When the DMSCS is in a superposition of distinct density profiles, we find that the decoherence rate can be appreciable. We then consider the novel special case in which the density is not in a superposition, but the phase of its field oscillation is; this is a property that cannot be decohered within the nonrelativistic framework. We find that if the probe particle and/or the DMSCS's velocity dispersion is slow, then the rate of decoherence of the phase is exponentially suppressed. However, if both the probe and the DMSCS's velocity dispersion are relativistic, then the phase can decohere more rapidly. As applications, we find that diffuse galactic axions with superposed phases are robust against decoherence, while dense boson stars and regions near black hole horizons are not, and we discuss implications for experiment.

categories: gr-qc astro-ph.CO hep-ph hep-th quant-ph


title: Disorder effects on Majorana zero modes: Kitaev chain versus semiconductor nanowire

abstract: Majorana zero modes in a superconductor-semiconductor nanowire have been extensively studied during the past decade. Disorder remains a serious problem, preventing the definitive observation of topological Majorana bound states. Thus, it is worthwhile to revisit the simple model, the Kitaev chain, and study the effects of weak and strong disorder on the Kitaev chain. By comparing the role of disorder in a Kitaev chain with that in a nanowire, we find that disorder affects both systems but in a nonuniversal manner. In general, disorder has a much stronger effect on the nanowire than the Kitaev chain, particularly for weak to intermediate disorder. For strong disorder, both the Kitaev chain and nanowire manifest random featureless behavior due to universal Anderson localization. Only the vanishing and strong disorder regimes are thus universal, manifesting respectively topological superconductivity and Anderson localization, but the experimentally relevant intermediate disorder regime is nonuniversal with the details dependent on the disorder realization in the system.


title: Driving galactic winds with magnetic fields at low and high redshift

abstract: Although playing a key role for our understanding of the evolution of galaxies, the exact way how observed galactic outflows are driven is still far from being understood and therefore our understanding of associated feedback mechanisms that control the evolution of galaxies is still plagued by many enigmas. In this work we present a simple toy model that can provide insight on how non-axis-symmetric instabilities in galaxies (bars, spiral-arms, warps) can lead to local exponential magnetic field growth by a radial flows beyond the equipartition value by at least two orders of magnitude on a time-scale of a few $100$ Myr. Our predictions show that the process can lead to galactic outflows in barred spiral galaxies with a mass loading factor $\eta \approx 0.1$, in agreement with our numerical simulations. Moreover, our outflow mechanism could contribute to an understanding of the large fraction of bared spiral galaxies that show signs of galactic outflows in the CHANG-ES survey. Extending our model shows the importance of such processes in high redshift galaxies by assuming equipartition between magnetic energy and turbulent energy. Simple estimates for the star formation rate (SFR) in our model together with cross-correlated masses from the star-forming main-sequence at redshifts $z\sim2$ allow us to estimate the outflow rate and mass loading factors by non-axis-symmetric instabilities and a subsequent radial inflow dynamo, giving mass loading factors of $\eta \approx 0.1$ for galaxies in the range of M$_{\star}=10^9 - 10^{12}$ M$_{\odot}$, in good agreement with recent results of Sinfoni and KMOS$^{3\mathrm{D}}$.


title: Topology of three-dimensional Dirac semimetals and generalized quantum spin Hall systems without gapless edge modes

abstract: Usually the quantum spin Hall states are expected to possess gapless, helical edge modes. Are there clean, non-interacting, quantum spin Hall states without gapless, edge modes? We show the generic, $n$-fold-symmetric, momentum planes of three-dimensional, stable Dirac semi-metals, which are orthogonal to the direction of nodal separation are examples of such generalized quantum spin Hall systems. We demonstrate that the planes lying between two Dirac points and the celebrated Bernevig-Zhang-Hughes model support identical quantized, non-Abelian Berry flux of magnitude $2 \pi$. Consequently, both systems exhibit spin-charge separation in response to electromagnetic, $\pi$-flux vortex. The Dirac points are identified as the unit-strength, monopoles of $SO(5)$ Berry connection, describing topological quantum phase transitions between generalized, quantum spin Hall and trivial insulators. Our work identifies precise bulk invariant and quantized response of Dirac semimetals and shows that many two-dimensional higher-order topological insulators can be understood as generalized quantum spin Hall systems, possessing gapped edge states.

categories: cond-mat.mes-hall cond-mat.mtrl-sci cond-mat.str-el hep-lat hep-th


title: Improving Damped Random Walk parameters for SDSS Stripe 82 Quasars with Pan-STARRS1

abstract: We use the Panoramic Survey Telescope and Rapid Response System 1 Survey (Pan-STARRS1, PS1) data to extend the Sloan Digital Sky Survey (SDSS) Stripe 82 quasar light curves. Combining PS1 and SDSS light curves provides a 15 yr baseline for 9248 quasars - 5 yr longer than prior studies that used only SDSS. We fit the light curves with the damped random walk (DRW) model model - a statistical description of their variability. We correlate the resulting DRW model parameters: asymptotic variability amplitude SF$_{\infty}$, and characteristic timescale $\tau$, with quasar physical properties - black hole mass, bolometric luminosity, and redshift. Using simulated light curves, we find that a longer baseline allows us to better constrain the DRW parameters. After adding PS1 data, the variability amplitude is a stronger function of the black hole mass, and has a weaker dependence on quasar luminosity. In addition, the characteristic timescale $\tau$ dependence on quasar luminosity is marginally weaker. We also make predictions for the fidelity of DRW model parameter retrieval when light curves will be further extended with Zwicky Transient Facility (ZTF) and the Rubin Observatory Legacy Survey of Space and Time (LSST) data. Finally, we show how updated DRW parameters offer an independent method of discovering changing-look quasar candidates (CLQSOs). The candidates are outliers in terms of differences in magnitude and scatter between SDSS and PS1 segments. We identify 40 objects (35 newly reported) with tenfold increase in variability timescale between SDSS and SDSS--PS1 data, which is due to a large change in brightness (over 0.5 mag) - characteristic for CLQSOs.


title: Measuring the topology of reionization with Betti numbers

abstract: The distribution of ionised hydrogen during the epoch of reionization (EoR) has a complex morphology. We propose to measure the three-dimensional topology of ionised regions using the Betti numbers. These quantify the topology using the number of components, tunnels and cavities in any given field. Based on the results for a set of reionization simulations we find that the Betti numbers of the ionisation field show a characteristic evolution during reionization, with peaks in the different Betti numbers characterising different stages of the process. The shapes of their evolutionary curves can be fitted with simple analytical functions. We also observe that the evolution of the Betti numbers shows a clear connection with the percolation of the ionized and neutral regions and differs between different reionization scenarios. Through these properties, the Betti numbers provide a more useful description of the topology than the widely studied Euler characteristic or genus. The morphology of the ionisation field will be imprinted on the redshifted 21-cm signal from the EoR. We construct mock image cubes using the properties of the low-frequency element of the future Square Kilometre Array and show that we can extract the Betti numbers from such datasets if an observation time of 1000 h is used. Even for a much shorter observation time of 100 h, some topological information can be extracted for the middle and later stages of reionization. We also find that the topological information extracted from the mock 21-cm observations can put constraints on reionization models.


title: Amorphous topological phases protected by continuous rotation symmetry

abstract: Protection of topological surface states by reflection symmetry breaks down when the boundary of the sample is misaligned with one of the high symmetry planes of the crystal. We demonstrate that this limitation is removed in amorphous topological materials, where the Hamiltonian is invariant on average under reflection over any axis due to continuous rotation symmetry. While the local disorder caused by the amorphous structure weakens the topological protection, we demonstrate that the edge remains protected from localization. In order to classify such phases we perform a systematic search over all the possible symmetry classes in two dimensions and construct the example models realizing each of the proposed topological phases. Finally, we compute the topological invariant of these phases as an integral along a meridian of the spherical Brillouin zone of an amorphous Hamiltonian.


title: Mid-Infrared Period$-$Luminosity Relations for Miras in the Large Magellanic Cloud

abstract: We present the mid-infrared (mid-IR) period$-$luminosity relations (PLRs) using over 1000 Mira variables in the Large Magellanic Cloud (LMC), for the four Wide-Field Infrared Survey Explorer (WISE) and the four Spitzer bands. These PLRs cover a mid-IR wavelength range from $3.4$ $\mu$m to $22$ $\mu$m and are presented separately for the oxygen-rich (O-rich) and carbon-rich (C-rich) Miras. These relations can be used to measure distances to individual O-rich and/or C-rich Mira stars with the accuracy of $5\%$ and $12\%$, respectively. They are the most accurate Mira PLRs in the mid-IR to date.


title: An Amelioration for the Sign Problem: Adiabatic Quantum Monte Carlo

abstract: We introduce the adiabatic quantum Monte Carlo (AQMC) method, where we gradually crank up the interaction strength, as an amelioration of the sign problem. It is motivated by the adiabatic theorem and will approach the true ground-state if the evolution time is long enough. We demonstrate that the AQMC enhances the average sign exponentially such that low enough temperatures can be accessed and ground-state properties probed. It is a controlled approximation that satisfies the variational theorem and provides an upper bound for the ground-state energy. We first benchmark the AQMC vis-\`a-vis the undoped Hubbard model on the square lattice which is known to be sign-problem-free within the conventional quantum Monte Carlo formalism. Next, we test the AQMC against the density-matrix-renormalization-group approach for the doped four-leg ladder Hubbard model and demonstrate its remarkable accuracy. As a nontrivial example, we apply our method to the Hubbard model at $p=1/8$ doping for a $16\times 8$ system and discuss its ground-state properties. We finally utilize our method and demonstrate the emergence of $U(1)_2\sim SU(2)_1$ topological order in a strongly correlated Chern insulator.


title: Signatures of ultrafast reversal of excitonic order in Ta$_2$NiSe$_5$

abstract: In the presence of electron-phonon coupling, an excitonic insulator harbors two degenerate ground states described by an Ising-type order parameter. Starting from a microscopic Hamiltonian, we derive the equations of motion for the Ising order parameter in the phonon coupled excitonic insulator Ta$_2$NiSe$_5$ and show that it can be controllably reversed on ultrashort timescales using appropriate laser pulse sequences. Using a combination of theory and time-resolved optical reflectivity measurements, we report evidence of such order parameter reversal in Ta$_2$NiSe$_5$ based on the anomalous behavior of its coherently excited order-parameter-coupled phonons. Our work expands the field of ultrafast order parameter control beyond spin and charge ordered materials.


title: Integrable triples in semisimple Lie algebras

abstract: We classify all integrable triples in simple Lie algebras, up to equivalence. The importance of this problem stems from the fact that for each such equivalence class one can construct the corresponding integrable hierarchy of bi-Hamiltonian PDE. The simplest integrable triple $(f,0,e)$ in $\mathfrak{sl}_2$ corresponds to the KdV hierarchy, and the triple $(f,0,e_\theta)$, where $f$ is the sum of negative simple root vectors and $e_\theta$ is the highest root vector of a simple Lie algebra, corresponds to the Drinfeld-Sokolov hierarchy.


title: Shadow analysis for rotating black holes in the presence of plasma for an expanding universe

abstract: We explore the structure of shadow for a Kerr-de Sitter black hole with a non-magnetized, pressureless plasma surrounding it. Specific plasma distributions are considered to separate the Hamilton-Jacobi equation and find the photon regions. An analytic formula describing the boundary curve of the shadow for such a black hole in an expanding universe for an observer at any finite point outside the horizon is derived. We observe deviations which are further explored by calculating the curvature radius at a particular point for such structures in the presence and absence of plasma and calculate the diameter of the shadow.


title: Classification of bad pixels of the Hawaii-2RG detector of the ASTROnomical NearInfraRed CAMera

abstract: ASTRONIRCAM is an infrared camera-spectrograph installed at the 2.5-meter telescope of the CMO SAI. The instrument is equipped with the HAWAII-2RG array. A bad pixels classification of the ASTRONIRCAM detector is proposed. The classification is based on histograms of the difference of consecutive non-destructive readouts of a flat field. Bad pixels are classified into 5 groups: hot (saturated on the first readout), warm (the signal accumulation rate is above the mean value by more than 5 standard deviations), cold (the rate is under the mean value by more than 5 standard deviations), dead (no signal accumulation), and inverse (having a negative signal accumulation in the first readouts). Normal pixels of the ASTRONIRCAM detector account for 99.6% of the total. We investigated the dependence between the amount of bad pixels and the number of cooldown cycles of the instrument. While hot pixels remain the same, the bad pixels of other types may migrate between groups. The number of pixels in each group stays roughly constant. We found that the mean and variance of the bad pixels amount in each group and the transitions between groups do not differ noticeably between normal or slow cooldowns.


title: Toward a Comprehensive Model of Snow Crystal Growth: 10. On the Molecular Dynamics of Structure Dependent Attachment Kinetics

abstract: I examine the molecular dynamics of ice growth from water vapor, focusing on how the attachment kinetics can be augmented by edge-dependent surface diffusion. Although there are significant uncertainties in developing an accurate physical model of this process, it is possible to make some reasonable estimates of surface diffusion rates and admolecule density enhancements, derived from our basic understanding of ice-crystal growth processes. A quantitative model suggests that edge-dependent surface diffusion could substantially enhance terrace nucleation on narrow faceted surfaces, especially at the onset of surface premelting. This result supports our hypothesized mechanism for structure-dependent attachment kinetics, which readily explains the changes in snow crystal growth morphology with temperature depicted in the well-known Nakaya diagram. Many of the model features described here may be amenable to further quantitative investigation using existing computational models of the molecular structure and dynamics of the ice surface.


title: Nonparametric approximation of conditional expectation operators

abstract: Given the joint distribution of two random variables $X,Y$ on some second countable locally compact Hausdorff space, we investigate the statistical approximation of the $L^2$-operator defined by $[Pf](x) := \mathbb{E}[ f(Y) \mid X = x ]$ under minimal assumptions. By modifying its domain, we prove that $P$ can be arbitrarily well approximated in operator norm by Hilbert-Schmidt operators acting on a reproducing kernel Hilbert space. This fact allows to estimate $P$ uniformly by finite-rank operators over a dense subspace even when $P$ is not compact. In terms of modes of convergence, we thereby obtain the superiority of kernel-based techniques over classically used parametric projection approaches such as Galerkin methods. This also provides a novel perspective on which limiting object the nonparametric estimate of $P$ converges to. As an application, we show that these results are particularly important for a large family of spectral analysis techniques for Markov transition operators. Our investigation also gives a new asymptotic perspective on the so-called kernel conditional mean embedding, which is the theoretical foundation of a wide variety of techniques in kernel-based nonparametric inference.


title: Gravitational scattering at the seventh order in $G$: nonlocal contribution at the sixth post-Newtonian accuracy

abstract: A recently introduced approach to the classical gravitational dynamics of binary systems involves intricate integrals (linked to a combination of nonlocal-in-time interactions with iterated $\frac1r$-potential scattering) which have so far resisted attempts at their analytical evaluation. By using computing techniques developed for the evaluation of multi-loop Feynman integrals (notably Harmonic Polylogarithms and Mellin transform) we show how to analytically compute all the integrals entering the nonlocal-in-time contribution to the classical scattering angle at the sixth post-Newtonian accuracy, and at the seventh order in Newton's constant, $G$ (corresponding to six-loop graphs in the diagrammatic representation of the classical scattering angle).


title: Optimal convergence rates in $L^2$ for a first order system least squares finite element method. Part I: homogeneous boundary conditions

abstract: We analyze a divergence based first order system least squares method applied to a second order elliptic model problem with homogeneous boundary conditions. We prove optimal convergence in the $L^2(\Omega)$ norm for the scalar variable. Numerical results confirm our findings.


title: Extensions of dissipative operators with closable imaginary part

abstract: Given a dissipative operator $A$ on a complex Hilbert space $\mathcal{H}$ such that the quadratic form $f\mapsto \mbox{Im}\langle f,Af\rangle$ is closable, we give a necessary and sufficient condition for an extension of $A$ to still be dissipative. As applications, we describe all maximally accretive extensions of strictly positive symmetric operators and all maximally dissipative extensions of a highly singular first-order operator on the interval.


title: Site-Selective Oxygen Vacancy Formation Derived from the Characteristic Crystal Structures of in Sn-Nb complex Oxides

abstract: Divalent tin oxides have attracted considerable attention as novel p-type oxide semiconductors, which are essential for realizing future oxide electronic devices. Recently, p-type Sn2Nb2O7 and SnNb2O6 were developed; however, enhanced hole mobility by reducing defect concentrations is required for practical use. In this work, we investigate the correlation between the formation of oxygen vacancy which may reduce the hole-generation efficiency and hole mobility, and the crystal structure in Sn-Nb complex oxides. Extended X-ray absorption fine structure spectroscopy and Rietveld analysis of x-ray diffraction revealed the preferential formation of oxygen vacancy at the O site bonded to the Sn ions in both the tin niobates. Moreover, a large amount of oxygen vacancy around the Sn ions were found in the p-type Sn2Nb2O7, thereby indicating the effect of oxygen vacancy to the low hole-generation efficiency. The dependence of the formation of oxygen vacancy on the crystal structure can be elucidated from the Sn-O bond strength that is evaluated based on the bond valence sum and Debye temperature. The differences in the bond strengths of the two Sn-Nb complex oxides are correlated through the steric hindrance of Sn2+ with asymmetric electron density distribution. This suggests the importance of the material design with a focus on the local structure around the Sn ions to prevent the formation of oxygen vacancy in p-type Sn2+ oxides.


title: Twist-induced control of near-field heat radiation between magnetic Weyl semimetals

abstract: Due to the large anomalous Hall effect, magnetic Weyl semimetals can support nonreciprocal surface plasmon polariton modes in the absence of an external magnetic field. This implies that magnetic Weyl semimetals can find novel application in (thermal) photonics. In this work, we consider the near-field radiative heat transfer between two magnetic Weyl semimetal slabs and show that the heat transfer can be controlled with a relative rotation of the parallel slabs. Thanks to the intrinsic nonreciprocity of the surface modes, this so-called twisting method does not require surface structuring like periodic gratings. The twist-induced control of heat transfer is due to the mismatch of the surface modes from the two slabs with a relative rotation.


title: ALMA detects molecular gas in the halo of the powerful radio galaxy TXS 0828+193

abstract: Both theoretical and observational results suggest that high-redshift radio galaxies (HzRGs) inhabit overdense regions of the universe and might be the progenitors of local, massive galaxies residing in the centre of galaxy clusters. In this paper we present CO(3-2) line observations of the HzRG TXS 0828+193 (z=2.57) and its environment using the Atacama Large Millimeter/submillimeter Array. In contrast to previous observations, we detect CO emission associated with the HzRG and derive a molecular gas mass of $(0.9\pm0.3)\times10^{10}\,\rm M_{\odot}$. Moreover, we confirm the presence of a previously detected off-source CO emitting region (companion #1), and detect three new potential companions. The molecular gas mass of each companion is comparable to that of the HzRG. Companion #1 is aligned with the axis of the radio jet and has stellar emission detected by Spitzer. Thus this source might be a normal star-forming galaxy or alternatively a result of jet-induced star formation. The newly found CO sources do not have counterparts in any other observing band and could be high-density clouds in the halo of TXS 0828+193 and thus potentially linked to the large-scale filamentary structure of the cosmic web.


title: Topological Weyl magnons and thermal Hall effect in layered honeycomb ferromagnets

abstract: In this work, we study the topological properties and magnon Hall effect of a three-dimensional ferromagnet in the ABC stacking honeycomb lattice, motivated by the recent inelastic neutron scattering study of CrI$_3$. We show that the magnon band structure and Chern numbers of the magnon branches are significantly affected by the interlayer coupling $J_c$, which moreover has a qualitatively different effect in the ABC stacking compared to the AA stacking adopted by other authors. The nontrivial Chern number of the lowest magnon band is stabilized by the next-nearest-neighbour Dzyaloshinsky-Moriya interaction in each honeycomb layer, resulting in the hopping term similar to that in the electronic Haldane model for graphene. However, we also find several gapless Weyl points, separating the non-equivalent Chern insulating phases, tuned by the ratio of the interlayer coupling $J_c$ and the third-neighbour Heisenberg interaction $J_3$. We further show that the topological character of magnon bands results in non-zero thermal Hall conductivity, whose sign and magnitude depend on $J_c$ and the intra-layer couplings. Since the interlayer coupling strength $J_c$ can be easily tuned by applying pressure to the quasi-2D material such as CrI$_3$, this provides a potential route to tuning the magnon thermal Hall effect in an experiment.


title: Inverse analysis of asteroseismic data: a review

abstract: Asteroseismology has emerged as the best way to characterize the global and internal properties of nearby stars. Often, this characterization is achieved by fitting stellar evolution models to asteroseismic observations. The star under investigation is then assumed to have the properties of the best-fitting model, such as its age. However, the models do not fit the observations perfectly. This is due to incorrect or missing physics in stellar evolution calculations, resulting in predicted stellar structures that are discrepant with reality. Through an inverse analysis of the asteroseismic data, it is possible to go further than fitting stellar models, and instead infer details about the actual internal structure of the star at some locations in its interior. Comparing theoretical and observed stellar structures then enables the determination of the locations where the stellar models have discrepant structure, and illuminates a path for improvements to our understanding of stellar evolution. In this invited review, we describe the methods of asteroseismic inversions, and outline the progress that is being made towards measuring the interiors of stars.


title: Diving below the spin-down limit: Constraints on gravitational waves from the energetic young pulsar PSR J0537-6910

abstract: We present a search for continuous gravitational-wave signals from the young, energetic X-ray pulsar PSR J0537-6910 using data from the second and third observing runs of LIGO and Virgo. The search is enabled by a contemporaneous timing ephemeris obtained using NICER data. The NICER ephemeris has also been extended through 2020 October and includes three new glitches. PSR J0537-6910 has the largest spin-down luminosity of any pulsar and is highly active with regards to glitches. Analyses of its long-term and inter-glitch braking indices provided intriguing evidence that its spin-down energy budget may include gravitational-wave emission from a time-varying mass quadrupole moment. Its 62 Hz rotation frequency also puts its possible gravitational-wave emission in the most sensitive band of LIGO/Virgo detectors. Motivated by these considerations, we search for gravitational-wave emission at both once and twice the rotation frequency. We find no signal, however, and report our upper limits. Assuming a rigidly rotating triaxial star, our constraints reach below the gravitational-wave spin-down limit for this star for the first time by more than a factor of two and limit gravitational waves from the $l=m=2$ mode to account for less than 14% of the spin-down energy budget. The fiducial equatorial ellipticity is limited to less than about 3e-5, which is the third best constraint for any young pulsar.


title: Towards a common performance and effectiveness terminology for digital proximity tracing applications

abstract: Digital proximity tracing (DPT) for Sars-CoV-2 pandemic mitigation is a complex intervention with the primary goal to notify app users about possible risk exposures to infected persons. Policymakers and DPT operators need to know whether their system works as expected in terms of speed or yield (performance) and whether DPT is making an effective contribution to pandemic mitigation (also in comparison to and beyond established mitigation measures, particularly manual contact tracing). Thereby, performance and effectiveness are not to be confused. Not only are there conceptual differences but also diverse data requirements. This article describes differences between performance and effectiveness measures and attempts to develop a terminology and classification system for DPT evaluation. We discuss key aspects for critical assessments of whether the integration of additional data measurements into DPT apps - beyond what is required to fulfill its primary notification role - may facilitate an understanding of performance and effectiveness of planned and deployed DPT apps. Therefore, the terminology and a classification matrix may offer some guidance to DPT system operators regarding which measurements to prioritize. DPT developers and operators may also make conscious decisions to integrate measures for epidemic monitoring but should be aware that this introduces a secondary purpose to DPT that is not part of the original DPT design. Ultimately, the integration of further information for epidemic monitoring into DPT involves a trade-off between data granularity and linkage on the one hand, and privacy on the other. Decision-makers should be aware of the trade-off and take it into account when planning and developing DPT notification and monitoring systems or intending to assess the added value of DPT relative to existing contact tracing systems.


title: Vapor-liquid equilibrium predictions of n-alkane/nitrogen mixtures using neural networks

abstract: Understanding fluid phase behavior in high pressure and high temperature conditions is crucial for developing high-fidelity simulations of chemically reacting flows in liquid-fueled combustion systems. The study of vapor-liquid equilibrium (VLE) curves also forms an integral part of the design and modeling of the control processes in chemical and oil-gas industries. The main objective of this study was to develop data-driven models to predict VLE of Type III binary mixtures involving long-chained n-alkanes and nitrogen. Two data-driven models have been proposed in this study, each of which was competent in estimating VLE for the binary systems of C10/N2 and C12/N2, at pressures ranging up to 50-60 MPa. Both the models showed better performance (less average absolute percentage error) in predicting equilibrium pressure of the binary mixtures as compared to the VLE modeled using Peng-Robinson equation of state (PR-EOS). The data-driven models were also able to correctly trace the change in curvature of the vapor phase composition, close to the mixture critical point, at high pressures and temperatures-a feature of the n-alkane/nitrogen system VLE which the PR-EOS model fails to capture.


title: 3-manifolds that bound no definite 4-manifold

abstract: We produce a rational homology 3-sphere that does not smoothly bound either a positive or negative definite 4-manifold. Such a 3-manifold necessarily cannot be rational homology cobordant to a Seifert fibered space or any 3-manifold obtained by Dehn surgery on a knot. The proof requires an analysis of short characteristic covectors in bimodular lattices.


title: Migration of Jupiter mass planets in low viscosity discs

abstract: Type-II migration of giant planets has a speed proportional to the disc's viscosity for values of the alpha viscosity parameter larger than 1.e-4 . At lower viscosities previous studies, based on 2D simulations have shown that migration can be very chaotic and often characterized by phases of fast migration. The reason is that in low-viscosity discs vortices appear due to the Rossby-wave instability at the edges of the gap opened by the planet. Migration is then determined by vortex-planet interactions. Our aim is to study migration in low viscosity 3D discs. We performed numerical simulations using 2D (including self-gravity) and 3D codes. After selecting disc masses for which self-gravity is not important, 3D simulations without self-gravity can be safely used. In our nominal simulation only numerical viscosity is present. We then performed simulations with prescribed viscosity to assess the threshold below which the new migration processes appear. We show that for alpha viscosity <= 1.e-5 two migration modes are possible which differ from classical Type-II migration, in the sense that they are not proportional to the disc's viscosity. The first occurs when the gap opened by the planet is not very deep. This occurs in 3D simulations and/or when a big vortex forms at the outer edge of the planetary gap, diffusing material into the gap. We call this type of migration "vortex-driven migration". This migration is very slow and cannot continue indefinitely, because eventually the vortex dissolves. The second migration mode occurs when the gap is deep so that the planet's eccentricity grows to a value ~0.2 due to inefficient eccentricity damping by corotation resonances. This second, faster migration mode appears to be typical of 2D models in discs with slower damping of temperature's perturbations.


title: On Using Classification Datasets to Evaluate Graph-Level Outlier Detection: Peculiar Observations and New Insights

abstract: It is common practice of the outlier mining community to repurpose classification datasets toward evaluating various detection models. To that end, often a binary classification dataset is used, where samples from one of the classes is designated as the inlier samples, and the other class is substantially down-sampled to create the ground-truth outlier samples. Graph-level outlier detection (GLOD) is rarely studied but has many potentially influential real-world applications. In this study, we identify an intriguing issue with repurposing graph classification datasets for GLOD. We find that ROC-AUC performance of the models changes significantly (flips from high to very low, even worse than random) depending on which class is down-sampled. Interestingly, ROC-AUCs on these two variants approximately sum to 1 and their performance gap is amplified with increasing propagations for a certain family of propagation based outlier detection models. We carefully study the graph embedding space produced by propagation based models and find two driving factors: (1) disparity between within-class densities which is amplified by propagation, and (2)overlapping support (mixing of embeddings) across classes. We also study other graph embedding methods and downstream outlier detectors, and find that the intriguing performance flip issue still widely exists but which version of the downsample achieves higher performance may vary. Thoughtful analysis over comprehensive results further deeper our understanding of the established issue.


title: Synthesis of Supervisors Robust Against Sensor Deception Attacks

abstract: We consider feedback control systems where sensor readings may be compromised by a malicious attacker intending on causing damage to the system. We study this problem at the supervisory layer of the control system, using discrete event systems techniques. We assume that the attacker can edit the outputs from the sensors of the system before they reach the supervisory controller. In this context, we formulate the problem of synthesizing a supervisor that is robust against the class of edit attacks on the sensor readings and present a solution methodology for this problem. This methodology blends techniques from games on automata with imperfect information with results from supervisory control theory of partially-observed discrete event systems. Necessary and sufficient conditions are provided for the investigated problem.


title: First Co-spatial Comparison of Stellar, Neutral-, and Ionized-gas Metallicities in a metal-rich galaxy: M83

abstract: We carry out a comparative analysis of the metallicities from the stellar, neutral-gas, and ionized-gas components in the metal-rich spiral galaxy M83. We analyze spectroscopic observations taken with the Hubble Space Telescope (HST), the Large Binocular Telescope (LBT) and the Very Large Telescope (VLT). We detect a clear depletion of the HI gas, as observed from the HI column densities in the nuclear region of this spiral galaxy. We find column densities of log[$N$(HI) cm$^{-2}$] $<$ 20.0 at galactocentric distances of $<$ 0.18 kpc, in contrast to column densities of log[$N$(HI) cm$^{-2}$] $\sim$ 21.0 in the galactic disk, a trend observed in other nearby spiral galaxies. We measure a metallicity gradient of $-$0.03 $\pm$ 0.01 dex kpc$^{-1}$ for the ionized gas, comparable to the metallicity gradient of a local benchmark of 49 nearby star-forming galaxies of $-$0.026 $\pm$ 0.002 dex kpc$^{-1}$. Our co-spatial metallicity comparison of the multi-phase gas and stellar populations shows excellent agreement outside of the nucleus of the galaxy hinting at a scenario where the mixing of newly synthesized metals from the most massive stars in the star clusters takes longer than their lifetimes ($\sim$10 Myr). Finally, our work shows that caution must be taken when studying the metallicity gradient of the neutral-gas component in star-forming galaxies, since this can be strongly biased, as these environments can be dominated by molecular gas. In these regions the typical metallicity tracers can provide inaccurate abundances as they may trace both the neutral- and molecular-gas components.


title: Topological superconductivity in nanowires proximate to a diffusive superconductor-magnetic insulator bilayer

abstract: We study semiconductor nanowires coupled to a bilayer of a disordered superconductor and a magnetic insulator, motivated by recent experiments reporting possible Majorana-zero-mode signatures in related architectures. Specifically, we pursue a quasiclassical Usadel equation approach that treats superconductivity in the bilayer self-consistently in the presence of spin-orbit scattering, magnetic-impurity scattering, and Zeeman splitting induced by both the magnetic insulator and a supplemental applied field. Within this framework we explore prospects for engineering topological superconductivity in a nanowire proximate to the bilayer. We find that a magnetic-insulator-induced Zeeman splitting, mediated through the superconductor alone, cannot induce a topological phase since the destruction of superconductivity (i.e., Clogston limit) preempts the required regime in which the nanowire's Zeeman energy exceeds the induced pairing strength. However, this Zeeman splitting does reduce the critical applied field needed to access the topological phase transition, with fields antiparallel to the magnetization of the magnetic insulator having an optimal effect. Finally, we show that magnetic-impurity scattering degrades the topological phase, and spin-orbit scattering, if present in the superconductor, pushes the Clogston limit to higher fields yet simultaneously increases the critical applied field strength.


title: Rendezvous Mission for Interstellar Objects Using a Solar Sail-based Statite Concept

abstract: Using the "statite," or static-satelite, concept -- an artificial satellite capable of hovering in place using a solar sail -- this work proposes to create a dynamic orbital slingshot in anticipation of Interstellar Objects (ISOs) passing through our solar system. The existence of these ISOs offers a unique scientific opportunity to answer fundamental scientific questions about the origin of solar system volatiles, the compositions of exo-solar systems, and the transfer rates of material between solar systems. However, due to their high heliocentric velocities and relatively short lead time, it may be extremely difficult to visit ISOs with current satellite propulsion systems. This work investigates the statite concept as applied to ISO missions and demonstrates potential configurations for optimal ISO flyby and rendezvous missions.


title: Newton acceleration on manifolds identified by proximal-gradient methods

abstract: Proximal methods are known to identify the underlying substructure of nonsmooth optimization problems. Even more, in many interesting situations, the output of a proximity operator comes with its structure at no additional cost, and convergence is improved once it matches the structure of a minimizer. However, it is impossible in general to know whether the current structure is final or not; such highly valuable information has to be exploited adaptively. To do so, we place ourselves in the case where a proximal gradient method can identify manifolds of differentiability of the nonsmooth objective. Leveraging this manifold identification, we show that Riemannian Newton-like methods can be intertwined with the proximal gradient steps to drastically boost the convergence. We prove the superlinear convergence of the algorithm when solving some nondegenerated nonsmooth nonconvex optimization problems. We provide numerical illustrations on optimization problems regularized by $\ell_1$-norm or trace-norm.


title: Minimal controllability time for systems with nonlinear drift under a compact convex state constraint

abstract: In this paper we estimate the minimal controllability time for a class of non-linear control systems with a bounded convex state constraint. An explicit expression is given for the controllability time if the image of the control matrix is of co-dimension one. A lower bound for the controllability time is given in the general case. The technique is based on finding a lower dimension system with the similar controllability properties as the original system. The controls corresponding to the minimal time, or time close to the minimal one, are discussed and computed analytically. The effectiveness of the proposed approach is illustrated by a few examples.


title: S-matrices from 4d worldvolume

abstract: We give an example of how conformal field theory methods in worldvolumes of dimension d > 2 could be used to calculate string-like amplitudes. The worldvolume propagator's logarithmic behavior is based on the use of worldvolume superspace (rather than the worldvolume ghost coordinates of a previous paper). Massless states have maximum spin d. Unitarity is not studied. We also touch on some related topics for ordinary spinning strings: calculating closed-string trees in worldsheet superspace, and zeroth-quantization of OSp(1|2) with zero-mode ghosts only.


title: Data Analysis for a Balloon Borne Compton Polarimeter

abstract: The Gamma RAy Polarimeter Experiment (GRAPE), a balloon borne Compton polarimeter for 50-500 keV gamma rays, was successfully flown for the second time in 2014. GRAPE consists of 24 collimated polarimeter modules. Each module is made up of 64 rectangular scintillators elements in a grid of 8 $\times$ 8 (36 plastic scintillators surrounded by 28 CsI scintillators). GRAPE uses photon scattered events between 2 detector elements to measure the polarization of the measured photons. GRAPE was flown from Fort Sumner, New Mexico on the morning of September 26$^{\text{th}}$ 2014. GRAPE was at float altitude for 14.4 hours and the Crab (primary target) was observed for 1.8 hours (initially planned for 8 hours). Background dominates the observations at flight altitudes, which depends on many flight and instrument parameters (altitude, instrument pointing, temperatures, etc). A technique based on the Principle Component Analysis (PCA) was implemented which used the varying parameters to estimate the background for the Crab. Our instrument design allows for an optical crosstalk (a known issue) that arises from light exiting the scintillator element and spilling to neighboring anodes. A model was developed to represent the crosstalk, which was incorporated in our simulations and the instrument response. The analysis was focused on the phase-integrated data (due to limited statistics of the shortened observation period). A power-law spectrum with photon index of 1.70$\pm$0.24 and a normalization of 1.01$\pm$1.35 ph/keV/s/cm$^2$ was measured. A polarization fraction of 0.43$\pm$0.4 and a polarization angle of 56$^\circ \pm$30$^\circ$ was measured for the phase integrated Crab observation in the 70-200 keV energy range. This result was not sufficiently significant enough to further our understanding of the emission mechanism.


title: Interpreting the dual Riccati equation through the LQ reproducing kernel

abstract: In this study, we provide an interpretation of the dual differential Riccati equation of Linear-Quadratic (LQ) optimal control problems. Adopting a novel viewpoint, we show that LQ optimal control can be seen as a regression problem over the space of controlled trajectories, and that the latter has a very natural structure as a reproducing kernel Hilbert space (RKHS). The dual Riccati equation then describes the evolution of the values of the LQ reproducing kernel when the initial time changes. This unveils new connections between control theory and kernel methods, a field widely used in machine learning.


title: BATTPOWER Toolbox: Memory-Efficient and High-Performance Multi-Period AC Optimal Power Flow Solver

abstract: With the introduction of massive renewable energy sources and storage devices, the traditional process of grid operation must be improved in order to be safe, reliable, fast responsive and cost efficient, and in this regard power flow solvers are indispensable. In this paper, we introduce an Interior Point-based (IP) Multi-Period AC Optimal Power Flow (MPOPF) solver for the integration of Stationary Energy Storage Systems (SESS) and Electric Vehicles (EV). The primary methodology is based on: 1) analytic and exact calculation of partial differential equations of the Lagrangian sub-problem, and 2) exploiting the sparse structure and pattern of the coefficient matrix of Newton-Raphson approach in the IP algorithm. Extensive results of the application of proposed methods on several benchmark test systems are presented and elaborated, where the advantages and disadvantages of different existing algorithms for the solution of MPOPF, from the standpoint of computational efficiency, are brought forward. We compare the computational performance of the proposed Schur-Complement algorithm with a direct sparse LU solver. The comparison is performed for two different applicational purposes: SESS and EV. The results suggest the substantial computational performance of Schur-Complement algorithm in comparison with that of a direct LU solver when the number of storage devices and optimisation horizon increase for both cases of SESS and EV. Also, some situations where computational performance is inferior are discussed.


title: Construction of low-energy symmetric Hamiltonians and Hubbard parameters for twisted multilayer systems using ab-initio input

abstract: A computationally efficient workflow for obtaining the low-energy symmetric tight-binding Hamiltonians for twisted multilayer systems is presented in this work. We apply this scheme to twisted bilayer graphene at the first magic angle. As initial step, the full-energy tight-binding Hamiltonian is generated by the Slater-Koster model with parameters fitted to ab-initio data at larger angles. Then, the low-energy symmetric four-band and twelve-band Hamiltonians are constructed using the maximum-localization procedure subjected to crystal and time-reversal-symmetry constraints. Finally, we compute extended Hubbard parameters for both models within the constrained random phase approximation (cRPA) for screening, which again respect the symmetries. The relevant data and results of this work are freely available via an online repository. Our workflow, exemplified in this work on twisted bilayer graphene, is straightforwardly transferable to other twisted multi-layer materials.


title: Physics and astrophysics of ultra-high energy cosmic rays: recent results from the Pierre Auger Observatory

abstract: Ultra-high-energy cosmic rays (UHECRs) are the highest energy messengers in the universe, with energies up to $10^{20}$ eV. Studies of astrophysical particles (nuclei, electrons, neutrinos and photons) at their highest observed energies have implications for fundamental physics as well as astrophysics. The primary particles interact in the atmosphere (or in the Earth) and generate extensive air showers. Analysis of those showers enables one not only to estimate the energy, direction and most probable mass of the primary cosmic particles, but also to obtain information about the properties of their hadronic interactions at energies more than one order of magnitude above that accessible with the current highest energy human-made accelerator. The Pierre Auger Observatory, located in the province of Mendoza, Argentina, is the largest cosmic ray experiment ever built. The observatory was designed as a hybrid detector covering an area of 3000 km$^2$ and has been taking data for almost twenty years. In this paper, a selection of the latest results is presented: the cosmic ray energy spectrum, studies of hadronic physics, searches for a directional anisotropy and studies of mass composition (including the photon and neutrino searches). Finally, the current upgrade ("AugerPrime") of the observatory, which is mostly aimed at improving the sensitivity to the particle type and mass of ultra-high energy cosmic rays, is described.


title: A Deep Chandra X-ray Observatory Study of the Millisecond Pulsar Population in the Globular Cluster Terzan 5

abstract: We present an analysis of 745.6 ks of archival Chandra X-ray Observatory Advanced CCD Imaging Spectrometer data accumulated between 2000 and 2016 of the millisecond pulsar (MSP) population in the rich Galactic globular cluster Terzan 5. Eight of the 37 MSPs with precise positions are found to have plausible X-ray source matches. Despite the deep exposure, the remaining MSPs are either marginally detected or have no obvious X-ray counterparts, which can be attributed to the typically soft thermal spectra of rotation-powered MSPs, which are strongly attenuated by the high intervening absorbing column (~$10^{22}$ cm$^{-2}$) towards the cluster, and in some instances severe source crowding/blending. For the "redback" MSP binaries, PSRs J1748-2446P and J1748-2446ad, and the "black widow" binary PSRs J1748-2446O, we find clear evidence for large-amplitude X-ray variability at the orbital period consistent with an intrabinary shock origin. The third redback MSP in the cluster, PSR J1748-2446A, shows large amplitude variations in flux on time scales of years, possibility due to state transitions or intense flaring episodes from the secondary star.


title: Optimal trading without optimal control

abstract: A hypothetical risk-neutral agent who trades to maximize the expected profit of the next trade will approximately exhibit long-term optimal behavior as long as this agent uses the vector $p = \nabla V (t, x)$ as effective microstructure alphas, where V is the Bellman value function for a smooth relaxation of the problem. Effective microstructure alphas are the steepest-ascent direction of V , equal to the generalized momenta in a dual Hamiltonian formulation. This simple heuristics has wide-ranging practical implications; indeed, most utility-maximization problems that require implementation via discrete limit-order-book markets can be treated by our method.


title: Automatic Recognition of Landmarks on Digital Dental Models

abstract: Fundamental to improving Dental and Orthodontic treatments is the ability to quantitatively assess and cross-compare their outcomes. Such assessments require calculating distances and angles from 3D coordinates of dental landmarks. The costly and repetitive task of hand-labelling dental models impedes studies requiring large sample size to penetrate statistical noise. We have developed techniques and software implementing these techniques to map out automatically, 3D dental scans. This process is divided into consecutive steps - determining a model's orientation, separating and identifying the individual tooth and finding landmarks on each tooth - described in this paper. Examples to demonstrate techniques and the software and discussions on remaining issues are provided as well. The software is originally designed to automate Modified Huddard Bodemham (MHB) landmarking for assessing cleft lip/palate patients. Currently only MHB landmarks are supported, but is extendable to any predetermined landmarks. This software, coupled with intra-oral scanning innovation, should supersede the arduous and error prone plaster model and caliper approach to Dental research and provide a stepping-stone towards automation of routine clinical assessments such as "index of orthodontic treatment need" (IOTN).


title: Hole-Like Fermi Surface in the Overdoped Non-Superconducting Bi$_{1.8}$Pb$_{0.4}$Sr$_2$CuO$_{6+\delta}$

abstract: In high-temperature cuprate superconductors, the anti-ferromagnetic spin fluctuations are thought to have a very important role in naturally producing an attractive interaction between the electrons in the $d$-wave channel. The connection between superconductivity and spin fluctuations is expected to be especially consequential at the overdoped end point of the superconducting dome. In some materials, that point seems to coincide with a Lifshitz transition, where the Fermi surface changes from the hole-like centered at ($\pi, \pi$) to the electron-like, centered at the $\Gamma$ point causing a loss of large momentum anti-ferromagnetic fluctuations. Here, we study the doping dependence of the electronic structure of Bi$_{1.8}$Pb$_{0.4}$Sr$_2$CuO$_{6+\delta}$ in angle-resolved photoemission and find that the superconductivity vanishes at lower doping than at which the Lifshitz transition occurs. This requires a more detailed re-examination of a spin-fluctuation scenario.


title: Bivariate Densities in Bayes Spaces: Orthogonal Decomposition and Spline Representation

abstract: A new orthogonal decomposition for bivariate probability densities embedded in Bayes Hilbert spaces is derived. It allows one to represent a density into independent and interactive parts, the former being built as the product of revised definitions of marginal densities and the latter capturing the dependence between the two random variables being studied. The developed framework opens new perspectives for dependence modelling (which is commonly performed through copulas), and allows for the analysis of dataset of bivariate densities, in a Functional Data Analysis perspective. A spline representation for bivariate densities is also proposed, providing a computational cornerstone for the developed theory.


title: Quantum limited sound attenuation in a dilute atomic Fermi gas

abstract: A new experiment (P. B. Patel et al., Science 370 (2020) 6521, 1222) involving resonantly interacting atoms confined by laser beams sheds light on momentum and energy diffusion in quantum fluids.


title: Weak Convergence of a Collection of Random Functions Defined by the Eigenvectors of Large Dimensional Random Matrices

abstract: For each $n$, let $U_n$ be Haar distributed on the group of $n\times n$ unitary matrices. Let $\bfx_{n,1},\ldots,\bfx_{n,m} $ denote orthogonal nonrandom unit vectors in ${\Bbb C}^n$ and let $\text{\bf u}_{n,k}=(u_k^1,\ldots,u_k^n)^*=U^*\text{\bf x}_{n,k}$, $k=1,\ldots,m$. Define the following functions on [0,1]: $X^{k,k}_n(t)=\sqrt n\sum_{i=1}^{[nt]}(|u_k^i|^2-\tfrac1n)$, $X_n^{k,k'}(t)=\sqrt{2n}\sum_{i=1}^{[nt]}\bar u_k^iu_{k'}^i$, $k<k'$. %("$\bar{\,\,\,\,\,}$" denoting conjugate). Then it is proven that $X_n^{k,k},\Re X_n^{k,k'}$, $\Im X_n^{k,k'}$, considered as random processes in $D[0,1]$, converge weakly, as $n\to\infty$, to $m^2$ independent copies of Brownian bridge. The same result holds for the $m(m+1)/2$ processes in the real case, where $O_n$ is real orthogonal Haar distributed and $\bfx_{n,i}\in{\Bbb R}^n$, with $\sqrt n$ in $X^{k,k}_n$ and $\sqrt{2n}$ in $X_n^{k,k'}$ replaced with $\sqrt{\frac n2}$ and $\sqrt{n}$, respectively. This latter result will be shown to hold for the matrix of eigenvectors of $M_n=(1/s)V_nV_n^T$ where $V_n$ is $n\times s$ consisting of the entries of $\{v_{ij}\},\ i,j=1,2,\ldots$, i.i.d. standardized and symmetrically distributed, with each $\bfx_{n,i}=\{\pm1/\sqrt n,\ldots,\pm1/\sqrt n\}$, and $n/s\to y>0$ as $n\to\infty$. This result extends the result in J.W. Silverstein {\sl Ann. Probab. \bf18} 1174-1194. These results are applied to the detection problem in sampling random vectors mostly made of noise and detecting whether the sample includes a nonrandom vector.


title: If Global or Local Investor Sentiments are Prone to Developing an Impact on Stock Returns, is there an Industry Effect?

abstract: This paper investigates the heterogeneous impacts of either Global or Local Investor Sentiments on stock returns. We study 10 industry sectors through the lens of 6 (so called) emerging countries: China, Brazil, India, Mexico, Indonesia and Turkey, over the 2000 to 2014 period. Using a panel data framework, our study sheds light on a significant effect of Local Investor Sentiments on expected returns for basic materials, consumer goods, industrial, and financial industries. Moreover, our results suggest that from Global Investor Sentiments alone, one cannot predict expected stock returns in these markets.


title: A differential perspective on Gradient Flows on ${\sf CAT}(\kappa)$-spaces and applications

abstract: We review the theory of Gradient Flows in the framework of convex and lower semicontinuous functionals on ${\sf CAT}(\kappa)$-spaces and prove that they can be characterized by the same differential inclusion $y_t'\in-\partial^-{\sf E}(y_t)$ one uses in the smooth setting and more precisely that $y_t'$ selects the element of minimal norm in $-\partial^-{\sf E}(y_t)$. This generalizes previous results in this direction where the energy was also assumed to be Lipschitz. We then apply such result to the Korevaar-Schoen energy functional on the space of $L^2$ and ${\sf CAT}(0)$ valued maps: we define the Laplacian of such $L^2$ map as the element of minimal norm in $-\partial^-{\sf E}(u)$, provided it is not empty. The theory of gradient flows ensures that the set of maps admitting a Laplacian is $L^2$-dense. Basic properties of this Laplacian are then studied.


title: Ranks of Tensor Networks for Eigenspace Projections and the Curse of Dimensionality

abstract: The hierarchical (multi-linear) rank of an order-$d$ tensor is key in determining the cost of representing a tensor as a (tree) Tensor Network (TN). In general, it is known that, for a fixed accuracy, a tensor with random entries cannot be expected to be efficiently approximable without the curse of dimensionality, i.e., a complexity growing exponentially with $d$. In this work, we show that the ground state projection (GSP) of a class of unbounded Hamiltonians can be approximately represented as an operator of low effective dimensionality that is independent of the (high) dimension $d$ of the GSP. This allows to approximate the GSP without the curse of dimensionality.


title: Dissecting a resonance wedge on heteroclinic bifurcations

abstract: This article studies routes to chaos occurring within a resonance wedge for a 3-parametric family of differential equations acting on a 3-sphere. Our starting point is an autonomous vector field whose flow exhibits a weakly attracting heteroclinic network made by two 1-dimensional connections and a 2-dimensional separatrix between two equilibria with different Morse indices. After changing the parameters, while keeping the 1-dimensional connections unaltered, we concentrate our study in the case where the 2-dimensional invariant manifolds of the equilibria do not intersect. We derive the first return map near the ghost of the attractor and we reduce the analysis of the system to a 2-dimensional map on the cylinder. Complex dynamical features arise from a discrete-time Bogdanov-Takens singularity, which may be seen as the organizing center by which one can obtain infinitely many attracting tori, strange attractors, infinitely many sinks and non-trivial contracting wandering domains. These dynamical phenomena occur within a structure that we call resonance wedge. As an application, we may see the "classical" Arnold tongue as a projection of a resonance wedge. The results are general, extend to other contexts and lead to a fine-tuning of the theory.


title: ARGES -- Advanced Renormalisation Group Equation Simplifier

abstract: We present the initial release of ARGES, a toolkit for obtaining renormalisation group equations in perturbation theory. As such, ARGES can handle any perturbatively renormalisable four-dimensional quantum field theory. Notable further features include a symbolic rather than numeric computation, input of unconventional scalar and Yukawa sectors, an interactive evaluation and disentanglement as well as capabilities to inject algebraic simplification rules. We provide a conceptual and practical introduction into ARGES, and highlight similarities and differences with complementary packages.


title: Antiferromagnetism and crystalline-electric field excitations in tetragonal NaCeO2

abstract: We investigate the crystal structure, magnetic properties, and crystalline-electric field of tetragonal, $I4_1/amd$, NaCeO$_2$. In this compound, Ce$^{3+}$ ions form a tetragonally elongated diamond lattice coupled by antiferromagnetic interactions ($\Theta_{CW} = -7.69$ K) that magnetically order below $T_N = 3.18$ K. The Ce$^{3+}$ $J = 5/2$ crystalline-electric field-split multiplet is studied via inelastic neutron scattering to parameterize a $J_{eff} = 1/2$ ground state doublet comprised of states possessing mixed $|m_z \rangle$ character. Neutron powder diffraction data reveal the onset of $A$-type antiferromagnetism with $\mu=0.57(2)$ $\mu_B$ moments aligned along the $c$-axis. The magnetic structure is consistent with the expectations of a frustrated Heisenberg $J_1$-$J_2$ model on the elongated diamond lattice with effective exchange values $J_1 > 4 J_2$ and $J_1 > 0$.


title: GD424 -- a helium-atmosphere white dwarf with a large amount of trace hydrogen in the process of digesting a rocky planetesimal

abstract: The photospheric metal pollution of white dwarfs is now well-established as the signature of the accretion of planetary debris. However, the origin of the trace hydrogen detected in many white dwarfs with helium atmospheres is still debated. Here, we report the analysis of GD424: a metal-polluted, helium-atmosphere white dwarf with a large amount of trace hydrogen. We determined the atmospheric parameters using a hybrid analysis that combines the sensitivity of spectroscopy to the atmospheric composition, $\log(\mathrm{H/He})$, with that of photometry and astrometry to the effective temperature, $T_{\mathrm{eff}}$, and surface gravity, $\log g$. The resulting white dwarf mass, radius, and cooling age are $M_{\mathrm{WD}}=0.77\pm0.01\,\mathrm{M}_{\odot}$, $R_{\mathrm{WD}}=0.0109\pm0.0001\,\mathrm{R}_{\odot}$, and $\tau_\mathrm{cool}=215\pm10$ Myr, respectively. We identified and measured the abundances of 11 photospheric metals and argue that the accretion event is most likely either in the increasing or steady state, and that the disrupted planetesimal resembles either CI chondrites or the bulk Earth in terms of its composition. We suggest that the observed $1.33\times 10^{22}$ g of trace hydrogen in GD424 were at least partly acquired through accretion of water-rich planetary debris in an earlier accretion episode.


title: Privacy Preservation for Wireless Sensor Networks in Healthcare: State of the Art, and Open Research Challenges

abstract: The advent of miniature biosensors has generated numerous opportunities for deploying wireless sensor networks in healthcare. However, an important barrier is that acceptance by healthcare stakeholders is influenced by the effectiveness of privacy safeguards for personal and intimate information which is collected and transmitted over the air, within and beyond these networks. In particular, these networks are progressing beyond traditional sensors, towards also using multimedia sensors, which raise further privacy concerns. Paradoxically, less research has addressed privacy protection, compared to security. Nevertheless, privacy protection has gradually evolved from being assumed an implicit by-product of security measures, and it is maturing into a research concern in its own right. However, further technical and socio-technical advances are needed. As a contribution towards galvanising further research, the hallmarks of this paper include: (i) a literature survey explicitly anchored on privacy preservation, it is underpinned by untangling privacy goals from security goals, to avoid mixing privacy and security concerns, as is often the case in other papers; (ii) a critical survey of privacy preservation services for wireless sensor networks in healthcare, including threat analysis and assessment methodologies; it also offers classification trees for the multifaceted challenge of privacy protection in healthcare, and for privacy threats, attacks and countermeasures; (iii) a discussion of technical advances complemented by reflection over the implications of regulatory frameworks; (iv) a discussion of open research challenges, leading onto offers of directions for future research towards unlocking the door onto privacy protection which is appropriate for healthcare in the twenty-first century.


title: Quantum sensing protocol for motionally chiral Rydberg atoms

abstract: A quantum sensing protocol is proposed for demonstrating the motion-induced chirality of circularly polarised Rydberg atoms. To this end, a cloud of Rydberg atoms is dressed by a bichromatic light field. This allows to exploit the long-lived ground states for implementing a Ramsey interferometer in conjunction with a spin echo pulse sequence for refocussing achiral interactions. Optimal parameters for the dressing lasers are identified. Combining a circularly polarised dipole transition in the Rydberg atom with atomic centre-of-mass motion, the system becomes chiral. The resulting discriminatory chiral energy shifts induced by a chiral mirror are estimated using a macroscopic quantum electrodynamics approach.


title: Active Deep Learning on Entity Resolution by Risk Sampling

abstract: While the state-of-the-art performance on entity resolution (ER) has been achieved by deep learning, its effectiveness depends on large quantities of accurately labeled training data. To alleviate the data labeling burden, Active Learning (AL) presents itself as a feasible solution that focuses on data deemed useful for model training. Building upon the recent advances in risk analysis for ER, which can provide a more refined estimate on label misprediction risk than the simpler classifier outputs, we propose a novel AL approach of risk sampling for ER. Risk sampling leverages misprediction risk estimation for active instance selection. Based on the core-set characterization for AL, we theoretically derive an optimization model which aims to minimize core-set loss with non-uniform Lipschitz continuity. Since the defined weighted K-medoids problem is NP-hard, we then present an efficient heuristic algorithm. Finally, we empirically verify the efficacy of the proposed approach on real data by a comparative study. Our extensive experiments have shown that it outperforms the existing alternatives by considerable margins. Using ER as a test case, we demonstrate that risk sampling is a promising approach potentially applicable to other challenging classification tasks.


title: The necessity and power of random, under-sampled experiments in biology

abstract: A vast array of transformative technologies developed over the past decade has enabled measurement and perturbation at ever increasing scale, yet our understanding of many systems remains limited by experimental capacity. Overcoming this limitation is not simply a matter of reducing costs with existing approaches; for complex biological systems it will likely never be possible to comprehensively measure and perturb every combination of variables of interest. There is, however, a growing body of work - much of it foundational and precedent setting - that extracts a surprising amount of information from highly under sampled data. For a wide array of biological questions, especially the study of genetic interactions, approaches like these will be crucial to obtain a comprehensive understanding. Yet, there is no coherent framework that unifies these methods, provides a rigorous mathematical foundation to understand their limitations and capabilities, allows us to understand through a common lens their surprising successes, and suggests how we might crystalize the key concepts to transform experimental biology. Here, we review prior work on this topic - both the biology and the mathematical foundations of randomization and low dimensional inference - and propose a general framework to make data collection in a wide array of studies vastly more efficient using random experiments and composite experiments.


title: The basal level of gene expression associated with chromatin loosening shapes Waddington landscapes and controls cell differentiation

abstract: The baseline level of transcription, which is variable and difficult to quantify, seriously complicates the normalization of comparative transcriptomic data, but its biological importance remains unappreciated. We show that this currently neglected ingredient is essential for controlling gene network multistability and therefore cellular differentiation. Basal expression is correlated to the degree of chromatin loosening measured by DNA accessibility, and systematically leads to cellular dedifferentiation as assessed by transcriptomic signatures, irrespective of the molecular and cellular tools used. Modeling gene network motifs formally involved in developmental bifurcations, reveals that the epigenetic landscapes of Waddington are restructured by the level of non specific expression, such that the attractors of progenitor and differentiated cells can be mutually exclusive. This mechanism is universal and holds beyond the particular nature of the genes involved, provided the multistable circuits are correctly described with autonomous basal expression. These results explain the relationships long established between gene expression noise, chromatin decondensation and cellular dedifferentiation, and highlight how heterochromatin maintenance is essential for preventing pathological cellular reprogramming, age-related diseases and cancer.


title: Modeling the subjective perspective of consciousness and its role in the control of behaviours

abstract: Consciousness has been hypothesized to operate as a global workspace, which accesses and integrates multimodal information in a unified manner, supports expectation violation monitoring and reduction, and the motivation, programming and control of action. One important yet open issue concerns how the subjective perspective at the core of consciousness, and subjective properties of manifestation of the environment in such perspective as an embodied experience, plays a role in such process. We operationalised the concept of subjective perspective using the principles of the Projective Consciousness Model (PCM), based on the projective geometrical concept of Field of Consciousness. We show how these principles can account for documented relationships between appraisal and distance as an inverse distance law, yield a generative model of affective and epistemic drives based on purely subjective parameters, such as the apparent size of objects, and can be generalised to implement Theory of Mind, in a manner that is consistent with simulation theory. We used simulations of artificial agents, based on psychological rationale, to demonstrate how different model parameters could generate a variety of emergent adaptive and maladaptive behaviours that are relevant to developmental and clinical psychology: the ability to be resilient in the face of obstacles through imaginary projections, the emergence of social approach and joint attention behaviours, the ability to take advantage of false beliefs attributed to others, the emergence of avoidance behaviours as observed in social anxiety disorders, the presence of restricted interests as observed in autism spectrum disorders. The simulation of agents was applied to a specific robotic context, and agents' behaviours were demonstrated by controlling the corresponding robots.


title: Representing Partial Programs with Blended Abstract Semantics

abstract: Synthesizing programs from examples requires searching over a vast, combinatorial space of possible programs. In this search process, a key challenge is representing the behavior of a partially written program before it can be executed, to judge if it is on the right track and predict where to search next. We introduce a general technique for representing partially written programs in a program synthesis engine. We take inspiration from the technique of abstract interpretation, in which an approximate execution model is used to determine if an unfinished program will eventually satisfy a goal specification. Here we learn an approximate execution model implemented as a modular neural network. By constructing compositional program representations that implicitly encode the interpretation semantics of the underlying programming language, we can represent partial programs using a flexible combination of concrete execution state and learned neural representations, using the learned approximate semantics when concrete semantics are not known (in unfinished parts of the program). We show that these hybrid neuro-symbolic representations enable execution-guided synthesizers to use more powerful language constructs, such as loops and higher-order functions, and can be used to synthesize programs more accurately for a given search budget than pure neural approaches in several domains.


title: Gravitational baryogenesis in DGP brane cosmology

abstract: We consider the imbalance of matter and antimatter by using a gravitational baryogenesis mechanism in the background of Dvali-Gabadadze-Porrati (DGP) brane cosmology. By taking into account a flat Friedmann-Lemaitre-Robertson-Walker (FLRW) metric in the DGP brane model, we find that for a radiation dominated universe, $w = 1/3$, the ratio of baryon number density to entropy from the gravitational baryogenesis is not zero, contrary to ordinary general relativity. Also, we study the ratio of baryon number density to entropy against the observational constraints in DGP cosmology.


title: Asymptotic approximation of a modified compressible Navier-Stokes system

abstract: We study the long time asymptotics of a modified compressible Navier-Stokes system (mcNS) inspired by the previous work of Hoff and Zumbrun. We introduce a new decomposition of the momentum field into its irrotational and incompressible parts, and a new method for approximating solutions of the heat equation in terms of Hermite functions in which $n^{th}$ order approximations can be computed for solutions with $n^{th}$ order moments. We then obtain existence of solutions to the mcNS system and show that the approximation in terms of Hermite functions gives the leading order terms in the long-time asymptotics, and under certain assumptions can be evaluated explicitly.


title: Linear-optical dynamics of one-dimensional anyons

abstract: We study the dynamics of bosonic and fermionic anyons defined on a one-dimensional lattice, under the effect of Hamiltonians quadratic in creation and annihilation operators, commonly referred to as linear optics. These anyonic models are obtained from deformations of the standard bosonic or fermionic commutation relations via the introduction of a non-trivial exchange phase between different lattice sites. We study the effects of the anyonic exchange phase on the usual bosonic and fermionic bunching behaviors. We show how to exploit the inherent Aharonov-Bohm effect exhibited by these particles to build a deterministic, entangling two-qubit gate and prove quantum computational universality in these systems. We define coherent states for bosonic anyons and study their behavior under two-mode linear-optical devices. In particular we prove that, for a specific value of the exchange factor, an anyonic mirror can generate cat states, an important resource in quantum information processing with continuous variables.


title: A single-shot measurement of time-dependent diffusion over sub-millisecond timescales using static field gradient NMR

abstract: Time-dependent diffusion behavior is probed over sub-millisecond timescales in a single shot using an NMR static gradient, time-incremented echo train acquisition (SG-TIETA) framework. The method extends the Carr-Purcell-Meiboom-Gill (CPMG) cycle under a static field gradient by discretely incrementing the $\pi$-pulse spacings to simultaneously avoid off-resonance effects and probe a range of timescales ($50 - 500$ microseconds). Pulse spacings are optimized based on a derived ruleset. The remaining effects of pulse inaccuracy are examined and found to be consistent across pure liquids of different diffusivities: water, decane, and octanol-1. A pulse accuracy correction is developed. Instantaneous diffusivity, $D_{\mathrm{inst}}(t)$, curves (i.e., half of the time derivative of the mean-squared displacement in the gradient direction), are recovered from pulse accuracy-corrected SG-TIETA decays using a model-free, log-linear least squares inversion method validated by Monte Carlo simulations. A signal-averaged, 1-minute experiment is described. A flat $D_{\mathrm{inst}}(t)$ is measured on pure dodecamethylcyclohexasiloxane whereas decreasing $D_{\mathrm{inst}}(t)$ are measured on yeast suspensions, consistent with the expected short-time $D_{\mathrm{inst}}(t)$ behavior for confining microstructural barriers on the order of microns.


title: Complete analytical solution to the quantum Yukawa potential

abstract: We present a complete analytical solution to the quantum problem of a particle in the Yukawa potential, using supersymmetry and a systematic expansion of the corresponding super-potentials. Results for the critical screening of the ground state improve in several figures existing results based on both numerical solutions and approximation methods. Our calculation to order $(a_{0}/D)^{2}$ for the squared ground state wavefunction at the origin, which enter in darkonium transitions, yields a correction of $\pi^{4}/216$ to results based on variational techniques.


title: Incommensurate magnetism mediated by Weyl fermions in NdAlSi

abstract: Emergent relativistic quasiparticles in Weyl semimetals are the source of exotic electronic properties such as surface Fermi arcs, the anomalous Hall effect, and negative magnetoresistance, all observed in real materials. Whereas these phenomena highlight the effect of Weyl fermions on the electronic transport properties, less is known about what collective phenomena they may support. Here, we report a new Weyl semimetal, NdAlSi that offers an example. Using neutron diffraction, we report a long-wavelength magnetic order in NdAlSi whose periodicity is linked to the nesting vector between two topologically non-trivial Fermi pockets, which we characterize using density functional theory and quantum oscillation measurements. Our work provides a rare example of Weyl fermions driving collective magnetism.


title: A unifying approach to non-minimal quasi-stationary distributions for one-dimensional diffusions

abstract: Convergence to non-minimal quasi-stationary distributions for one-dimensional diffusions is studied. We give a method of reducing the convergence to the tail behavior of the lifetime via a property which we call the first hitting uniqueness. We apply the results to Kummer diffusions with negative drifts and give a class of initial distributions converging to each non-minimal quasi-stationary distribution.


title: A triangle process on regular graphs

abstract: Switches are operations which make local changes to the edges of a graph, usually with the aim of preserving the vertex degrees. We study a restricted set of switches, called triangle switches. Each triangle switch creates or deletes at least one triangle. Triangle switches can be used to define Markov chains which generate graphs with a given degree sequence and with many more triangles (3-cycles) than is typical in a uniformly random graph with the same degrees. We show that the set of triangle switches connects the set of all $d$-regular graphs on $n$ vertices, for all $d\geq 3$. Hence, any Markov chain which assigns positive probability to all triangle switches is irreducible on these graphs. We also investigate this question for 2-regular graphs.


title: The Nirenberg problem on high dimensional half spheres: The effect of pinching conditions

abstract: In this paper we study the Nirenberg problem on standard half spheres $(\mathbb{S}^n_+,g), \, n \geq 5$, which consists of finding conformal metrics of prescribed scalar curvature and zero boundary mean curvature on the boundary. This problem amounts to solve the following boundary value problem involving the critical Sobolev exponent: \begin{equation*} (\mathcal{P}) \quad \begin{cases} -\D_{g} u \, + \, \frac{n(n-2)}{4} u \, = K \, u^{\frac{n+2}{n-2}},\, u > 0 & \mbox{in } \mathbb{S}^n_+, \frac{\partial u}{\partial \nu }\, =\, 0 & \mbox{on } \partial \mathbb{S}^n_+. \end{cases} \end{equation*} where $K \in C^3(\mathbb{S}^n_+)$ is a positive function. This problem has a variational structure but the related Euler-Lagrange functional $J_K$ lacks compactness. Indeed it admits \emph{critical points at infinity}, which are \emph{limits} of non compact orbits of the (negative) gradient flow. Through the construction of an appropriate \emph{pseudogradient} in the \emph{neighborhood at infinity}, we characterize these \emph{critical points at infinity}, associate to them an index, perform a \emph{Morse type reduction} of the functional $J_K$ in their neighborhood and compute their contribution to the difference of topology between the level sets of $J_K$, hence extending the full Morse theoretical approach to this \emph{non compact variational problem}. Such an approach is used to prove, under various pinching conditions, some existence results for $(\mathcal{P})$ on half spheres of dimension $n \geq 5$.


title: Li-Yau inequalities for general non-local diffusion equations via reduction to the heat kernel

abstract: We establish a reduction principle to derive Li-Yau inequalities for non-local diffusion problems in a very general framework, which covers both the discrete and continuous setting. Our approach is not based on curvature-dimension inequalities but on heat kernel representations of the solutions and consists in reducing the problem to the heat kernel. As an important application we solve a long-standing open problem by obtaining a Li-Yau inequality for positive solutions $u$ to the fractional (in space) heat equation of the form $(-\Delta)^{\beta/2}(\log u)\leq C/t$, where $\beta\in (0,2)$. We also illustrate our general result with an example in the discrete setting by proving a sharp Li-Yau inequality for diffusion on a complete graph.


title: Detecting Hate Speech in Memes Using Multimodal Deep Learning Approaches: Prize-winning solution to Hateful Memes Challenge

abstract: Memes on the Internet are often harmless and sometimes amusing. However, by using certain types of images, text, or combinations of both, the seemingly harmless meme becomes a multimodal type of hate speech -- a hateful meme. The Hateful Memes Challenge is a first-of-its-kind competition which focuses on detecting hate speech in multimodal memes and it proposes a new data set containing 10,000+ new examples of multimodal content. We utilize VisualBERT -- which meant to be the BERT of vision and language -- that was trained multimodally on images and captions and apply Ensemble Learning. Our approach achieves 0.811 AUROC with an accuracy of 0.765 on the challenge test set and placed third out of 3,173 participants in the Hateful Memes Challenge.


title: A Plethora of Polynomials: A Toolbox for Counting Problems

abstract: A wide variety of problems in combinatorics and discrete optimization depend on counting the set $S$ of integer points in a polytope, or in some more general object constructed via discrete geometry and first-order logic. We take a tour through numerous problems of this type. In particular, we consider families of such sets $S_t$ depending on one or more integer parameters $t$, and analyze the behavior of the function $f(t)=|S_t|$. In the examples that we investigate, this function exhibits surprising polynomial-like behavior. We end with two broad theorems detailing settings where this polynomial-like behavior must hold. The plethora of examples illustrates the framework in which this behavior occurs and also gives an intuition for many of the proofs, helping us create a toolbox for counting problems like these.