title: Global strong solutions in $\mathbb{R}^3$ for ionic Vlasov-Poisson systems

abstract: Systems of Vlasov-Poisson type are kinetic models describing dilute plasma. The structure of the model differs according to whether it describes the electrons or positively charged ions in the plasma. In contrast to the electron case, where the well-posedness theory for Vlasov-Poisson systems is well established, the well-posedness theory for ion models has been investigated more recently. In this article, we prove global well-posedness for two Vlasov-Poisson systems for ions, posed on the whole three-dimensional Euclidean space $\mathbb{R}^3$, under minimal assumptions on the initial data and the confining potential.


title: Anisotropy-induced soliton excitation in magnetized strong-rung spin ladders

abstract: We report low temperature electron spin resonance experimental and theoretical studies of an archetype $S=1/2$ strong-rung spin ladder material (C$_{5}$H$_{12}$N)$_{2}$CuBr$_{4}$. Unexpected dynamics is detected deep in the Tomonaga-Luttinger spin liquid regime. Close to the point where the system is half-magnetized (and believed to be equivalent to a gapless easy plane chain in zero field) we observed orientation-dependent spin gap and anomalous $g$-factor values. Field theoretical analysis demonstrates that the observed low-energy excitation modes in magnetized (C$_{5}$H$_{12}$N)$_{2}$CuBr$_{4}$ are solitonic excitations caused by Dzyaloshinskii-Moriya interaction presence.


title: Interpretation of bulk viscosity as the generalized Chaplygin gas

abstract: The cosmological observations suggest that the presently accelerating universe should be filled by an exotic form of matter, violating the strong energy condition, of unknown nature and origin. We propose the viscous dark matter of a source of acceleration in the form of Chaplygin gas which is characterized by equation of state in the phenomenological form $p=-\frac{A}{\rho^{\alpha}}$, where $p$ and $\rho$ are pressure and energy density respectively ($A$ and $\alpha$ are constants). Chaplygin gas is interpreted in terms of viscous matter and without the cosmological constant. The acceleration effect is caused only by viscosity in this class of cosmological models. We show that bulk viscosity effects introduced to the standard FRW cosmology give rise to the natural unification of both dark matter and dark energy. We show that dust viscous cosmological models are structurally stable if $m < 1/2$ ($1+\alpha=1/2-m$).


title: Understanding Notions of Stationarity in Non-Smooth Optimization

abstract: Many contemporary applications in signal processing and machine learning give rise to structured non-convex non-smooth optimization problems that can often be tackled by simple iterative methods quite effectively. One of the keys to understanding such a phenomenon---and, in fact, one of the very difficult conundrums even for experts---lie in the study of "stationary points" of the problem in question. Unlike smooth optimization, for which the definition of a stationary point is rather standard, there is a myriad of definitions of stationarity in non-smooth optimization. In this article, we give an introduction to different stationarity concepts for several important classes of non-convex non-smooth functions and discuss the geometric interpretations and further clarify the relationship among these different concepts. We then demonstrate the relevance of these constructions in some representative applications and how they could affect the performance of iterative methods for tackling these applications.


title: Holographic Projection of Electromagnetic Maxwell Theory

abstract: The 4D Maxwell theory with single-sided planar boundary is considered. As a consequence of the presence of the boundary, two broken Ward identities are recovered, which, on-shell, give rise to two conserved currents living on the edge. A Ka\c{c}-Moody algebra formed by a subset of the bulk fields is obtained with central charge proportional to the inverse of the Maxwell coupling constant, and the degrees of freedom of the boundary theory are identified as two vector fields, also suggesting that the 3D theory should be a gauge theory. Finally the holographic contact between bulk and boundary theory is reached in two inequivalent ways, both leading to a unique 3D action describing a new gauge theory of two coupled vector fields with a topological Chern-Simons term with massive coefficient. In order to check that the 3D projection of 4D Maxwell theory is well defined, we computed the energy-momentum tensor and the propagators. The role of discrete symmetries is briefly discussed.


title: The path to instability in compact multi-planetary systems

abstract: The dynamical stability of tightly packed exoplanetary systems remains poorly understood. While for a two-planet system a sharp stability boundary exists, numerical simulations of three and more planet systems show that they can experience instability on timescales up to billions of years. Moreover, an exponential trend between the planet orbital separation measured in units of Hill radii and the survival time has been reported. While these findings have been observed in numerous numerical simulations, little is known of the actual mechanism leading to instability. Contrary to a constant diffusion process, planetary systems seem to remain dynamically quiescent for most of their lifetime before a very short unstable phase. In this work, we show how the slow chaotic diffusion due to the overlap of three-body resonances dominates the timescale leading to the instability for initially coplanar and circular orbits. While the last instability phase is related to scattering due to two-planet mean motion resonances (MMR), for circular orbits the two-planets MMR are too far separated to destabilize systems initially away from them. We develop an analytical model to generalize the empirical trend obtained for equal mass and equally-spaced planets to general systems. We obtain an analytical estimate of the survival time consistent with simulations over four orders of magnitude for the planet to star mass ratio $\epsilon$, and 6 to 8 orders of magnitude for the instability time. We also confirm that measuring the orbital spacing in terms of Hill radii is not adapted and that the right spacing unit scales as $\epsilon^{1/4}$. We predict that beyond a certain spacing, the three-planet resonances are not overlapped, which results in an increase of the survival time. We finally discuss the extension of our result to more general systems, containing more planets on initially non circular orbits.


title: Layerwise learning for quantum neural networks

abstract: With the increased focus on quantum circuit learning for near-term applications on quantum devices, in conjunction with unique challenges presented by cost function landscapes of parametrized quantum circuits, strategies for effective training are becoming increasingly important. In order to ameliorate some of these challenges, we investigate a layerwise learning strategy for parametrized quantum circuits. The circuit depth is incrementally grown during optimization, and only subsets of parameters are updated in each training step. We show that when considering sampling noise, this strategy can help avoid the problem of barren plateaus of the error surface due to the low depth of circuits, low number of parameters trained in one step, and larger magnitude of gradients compared to training the full circuit. These properties make our algorithm preferable for execution on noisy intermediate-scale quantum devices. We demonstrate our approach on an image-classification task on handwritten digits, and show that layerwise learning attains an 8% lower generalization error on average in comparison to standard learning schemes for training quantum circuits of the same size. Additionally, the percentage of runs that reach lower test errors is up to 40% larger compared to training the full circuit, which is susceptible to creeping onto a plateau during training.


title: Water Dynamics Around Proteins: T- and R-States of Hemoglobin and Melittin

abstract: The water dynamics, as characterized by the local hydrophobicity (LH), is investigated for tetrameric hemoglobin and dimeric melittin. For the T0 to R0 transition in Hb it is found that LH provides additional molecular-level insight into the Perutz mechanism, i.e., the breaking and formation of salt bridges at the alpha1 / beta2 and alpha2 / beta1 interface is accompanied by changes in LH. For Hb in cubic water boxes with 90 Aengstroem and 120 Aengstroem edge length it is observed that following a decrease in LH as a consequence of reduced water density or change of water orientation at the protein/water interface the alpha / beta interfaces are destabilized; this is a hallmark of the Perutz stereochemical model for the T to R transition in Hb. The present work thus provides a dynamical view of the classical structural model relevant to the molecular foundations of Hb function. For dimeric melittin, earlier results by Cheng and Rossky (Nature, 1998, 392, 696-699) are confirmed and interpreted on the basis of LH from simulations in which the protein structure is frozen. For the flexible melittin dimer the changes in the local hydration can be as much as 30 % than for the rigid dimer, reflecting the fact that protein and water dynamics are coupled.


title: A search for fast radio burst-like emission from Fermi gamma-ray bursts

abstract: We report the results of the rapid follow-up observations of gamma-ray bursts (GRBs) detected by the Fermi satellite to search for associated fast radio bursts. The observations were conducted with the Australian Square Kilometre Array Pathfinder at frequencies from 1.2-1.4 GHz. A set of 20 bursts, of which four were short GRBs, were followed up with a typical latency of about one minute, for a duration of up to 11 hours after the burst. The data was searched using 4096 dispersion measure trials up to a maximum dispersion measure of 3763 pc cm$^{-3}$, and for pulse widths $w$ over a range of duration from 1.256 to 40.48 ms. No associated pulsed radio emission was observed above $26 {\rm Jy ms} (w/1 {\rm ms})^{-1/2}$ for any of the 20 GRBs.


title: Explicit uniform bounds for Brauer groups of singular K3 surfaces

abstract: Let $k$ be a number field. We give an explicit bound, depending only on $[k:\mathbf{Q}]$ and the discriminant of the N\'{e}ron--Severi lattice, on the size of the Brauer group of a K3 surface $X/k$ that is geometrically isomorphic to the Kummer surface attached to a product of isogenous CM elliptic curves. As an application, we show that the Brauer--Manin set for such a variety is effectively computable. Conditional on GRH, we can also make the explicit bound depend only on $[k:\mathbf{Q}]$ and remove the condition that the elliptic curves be isogenous. In addition, we show how to obtain a bound, depending only on $[k:\mathbf{Q}]$, on the number of $\mathbf{C}$-isomorphism classes of singular K3 surfaces defined over $k$, thus proving an effective version of the strong Shafarevich conjecture for singular K3 surfaces.


title: Crossings between non-homotopic edges

abstract: We call a multigraph {\em non-homotopic} if it can be drawn in the plane in such a way that no two edges connecting the same pair of vertices can be continuously transformed into each other without passing through a vertex, and no loop can be shrunk to its end-vertex in the same way. It is easy to see that a non-homotopic multigraph on $n>1$ vertices can have arbitrarily many edges. We prove that the number of crossings between the edges of a non-homotopic multigraph with $n$ vertices and $m>4n$ edges is larger than $c\frac{m^2}{n}$ for some constant $c>0$, and that this bound is tight up to a polylogarithmic factor. We also show that the lower bound is not asymptotically sharp as $n$ is fixed and $m$ tends to infinity.


title: Equipartition of Energy Defines the Size-Thickness Relationship in Liquid-Exfoliated Nanosheets

abstract: Liquid phase exfoliation is a commonly used method to produce 2D nanosheets from a range of layered crystals. However, such nanosheets display broad size and thickness distributions and correlations between area and thickness, issues that limit nanosheet application potential. To understand the factors controlling the exfoliation process, we have liquid-exfoliated 11 different layered materials, size-selecting each into fractions before using AFM to measure the nanosheet length, width, and thickness distributions for each fraction. The resultant data show a clear power-law scaling of nanosheet area with thickness for each material. We have developed a simple nonequilibrium thermodynamics-based model predicting that the power-law prefactor is proportional to both the ratios of in-plane-tearing/out-of-plane-peeling energies and in-plane/out-of-plane moduli. By comparing the experimental data with the modulus ratio calculated from first-principles, we find close agreement between experiment and theory. This supports our hypothesis that energy equipartition holds between nanosheet tearing and peeling during sonication-assisted exfoliation.

categories: physics.chem-ph cond-mat.mes-hall cond-mat.stat-mech physics.app-ph


title: Clustering and Halo Abundances in Early Dark Energy Cosmological Models

abstract: LCDM cosmological models with Early Dark Energy (EDE) have been proposed to resolve tensions between the Hubble constant H0 = 100h km/s/Mpc measured locally, giving h ~ 0.73, and H0 deduced from Planck cosmic microwave background (CMB) and other early universe measurements plus LCDM, giving h ~ 0.67. EDE models do this by adding a scalar field that temporarily adds dark energy equal to about 10% of the cosmological energy density at the end of the radiation-dominated era at redshift z ~ 3500. Here we compare linear and nonlinear predictions of a Planck-normalized LCDM model including EDE giving h = 0.728 with those of standard Planck-normalized LCDM with h = 0.678. We find that nonlinear evolution reduces the differences between power spectra of fluctuations at low redshifts. As a result, at z = 0 the halo mass functions on galactic scales are nearly the same, with differences only 1-2%. However, the differences dramatically increase at high redshifts. The EDE model predicts 50% more massive clusters at z = 1 and twice more galaxy-mass halos at z = 4. Even greater increases in abundances of galaxy-mass halos at higher redshifts may make it easier to reionize the universe with EDE. Predicted galaxy abundances and clustering will soon be tested by JWST observations. Positions of baryonic acoustic oscillations (BAOs) and correlation functions differ by about 2% between the models -- an effect that is not washed out by nonlinearities. Both standard LCDM and the EDE model studied here agree well with presently available acoustic-scale observations, but DESI and Euclid measurements will provide stringent new tests.


title: Can Autonomous Vehicles Identify, Recover From, and Adapt to Distribution Shifts?

abstract: Out-of-training-distribution (OOD) scenarios are a common challenge of learning agents at deployment, typically leading to arbitrary deductions and poorly-informed decisions. In principle, detection of and adaptation to OOD scenes can mitigate their adverse effects. In this paper, we highlight the limitations of current approaches to novel driving scenes and propose an epistemic uncertainty-aware planning method, called \emph{robust imitative planning} (RIP). Our method can detect and recover from some distribution shifts, reducing the overconfident and catastrophic extrapolations in OOD scenes. If the model's uncertainty is too great to suggest a safe course of action, the model can instead query the expert driver for feedback, enabling sample-efficient online adaptation, a variant of our method we term \emph{adaptive robust imitative planning} (AdaRIP). Our methods outperform current state-of-the-art approaches in the nuScenes \emph{prediction} challenge, but since no benchmark evaluating OOD detection and adaption currently exists to assess \emph{control}, we introduce an autonomous car novel-scene benchmark, \texttt{CARNOVEL}, to evaluate the robustness of driving agents to a suite of tasks with distribution shifts.


title: Multi-orbital charge density wave excitations and concomitant phonon anomalies in Bi$_2$Sr$_2$LaCuO$_{6+\delta}$

abstract: Charge density waves (CDWs) are ubiquitous in under-doped cuprate superconductors. As a modulation of the valence electron density, CDWs in hole-doped cuprates possess both Cu-3d and O-2p orbital character owing to the strong hybridization of these orbitals near the Fermi level. Here, we investigate under-doped Bi$_2$Sr$_{1.4}$La$_{0.6}$CuO$_{6+\delta}$ using resonant inelastic X-ray scattering (RIXS) and find that a short-range CDW exists at both Cu and O sublattices in the copper-oxide (CuO2) planes with a comparable periodicity and correlation length. Furthermore, we uncover bond-stretching and bond-buckling phonon anomalies concomitant to the CDWs. Comparing to slightly over-doped Bi$_2$Sr$_{1.8}$La$_{0.2}$CuO$_{6+\delta}$, where neither CDWs nor phonon anomalies appear, we highlight that a sharp intensity anomaly is induced in the proximity of the CDW wavevector (QCDW) for the bond-buckling phonon, in concert with the diffused intensity enhancement of the bond-stretching phonon at wavevectors much greater than QCDW. Our results provide a comprehensive picture of the quasi-static CDWs, their dispersive excitations, and associated electron-phonon anomalies, which are key for understanding the competing electronic instabilities in cuprates.


title: Two-Way Source-Channel Coding

abstract: We propose an adaptive lossy joint source-channel coding (JSCC) scheme for sending correlated sources over two-terminal discrete-memoryless two-way channels (DM-TWCs). The main idea is to couple the independent operations of the terminals via an adaptive coding mechanism, which can mitigate cross-interference resulting from simultaneous channel transmissions and concurrently exploit the sources' correlation to reduce the end-to-end reconstruction distortions. Our adaptive JSCC scheme not only subsumes existing lossy coding methods for two-way simultaneous communication but also improves their performance. Furthermore, we derive outer bounds for our two-way lossy transmission problem and establish complete JSCC theorems in some special settings. In these special cases, a non-adaptive separate source-channel coding (SSCC) scheme achieves the optimal performance, thus simplifying the design of the source-channel communication system.


title: Interfacing topological insulators and ferrimagnets: Bi$_2$Te$_3$ and Fe$_3$O$_4$ heterostructures grown by molecular beam epitaxy

abstract: Relying on the magnetism induced by the proximity effect in heterostructures of topological insulators and magnetic insulators is one of the promising routes to achieve the quantum anomalous Hall effect. Here we investigate heterostructures of Bi$_2$Te$_3$ and Fe$_3$O$_4$. By growing two different types of heterostructures by molecular beam epitaxy, Fe$_3$O$_4$ on Bi$_2$Te$_3$ and Bi$_2$Te$_3$ on Fe$_3$O$_4$, we explore differences in chemical stability, crystalline quality, electronic structure, and transport properties. We find the heterostructure Bi$_2$Te$_3$ on Fe$_3$O$_4$ to be a more viable approach, with transport signatures in agreement with a gap opening in the topological surface states.


title: Limit theory of combinatorial optimization for random geometric graphs

abstract: In the random geometric graph $G(n,r_n)$, $n$ vertices are placed randomly in Euclidean $d$-space and edges are added between any pair of vertices distant at most $r_n$ from each other. We establish strong laws of large numbers (LLNs) for a large class of graph parameters, evaluated for $G(n,r_n)$ in the thermodynamic limit with $nr_n^d =$ const., and also in the dense limit with $n r_n^d \to \infty$, $r_n \to 0$. Examples include domination number, independence number, clique-covering number, eternal domination number and triangle packing number. The general theory is based on certain subadditivity and superadditivity properties, and also yields LLNs for other functionals such as the minimum weight for the travelling salesman, spanning tree, matching, bipartite matching and bipartite travelling salesman problems, for a general class of weight functions with at most polynomial growth of order $d-\varepsilon$, under thermodynamic scaling of the distance parameter.


title: The calculation of the Mittag-Leffler function

abstract: The problem of calculating the Mittag-Leffler function $E_{\rho,\mu} (z)$ is considered in the paper. To solve this problem integral representations for the function $E_{\rho,\mu}(z)$ are transformed in such a way that they could not contain complex variables and parameters. Integral representations written in this form allow one to use standard methods of numerical integration to calculate integrals contained in them. To verify the correctness of the integral representations obtained the function $E_{\rho,\mu}(z)$ was calculated both with the use of obtained formulas and with the use of known representations of the Mittag-Leffler function. The calculation results demonstrate their exact matching. This fact is indicative of the correctness of new integral representations of the function $E_{\rho,\mu}(z)$ that were obtained.


title: Modeling liquid migration in active swollen gel spheres

abstract: Liquid migration in active soft solids is a very common phenomenon in Nature at different scales: from cells to leaves. It can be caused by mechanical as well as chemical actions. The work focuses on the migration of liquid provoked by remodeling processes in an active impermeable gel sphere. Within this context, we present a consistent mathematical theory capable to gain a deep understanding of the phenomenon in both steady and transient conditions.


title: First Search for the $K_L \to \pi^0 \gamma$ Decay

abstract: We report the first search for the $K_L \to \pi^0 \gamma$ decay, which is forbidden by Lorentz invariance, using the data from 2016 to 2018 at the J-PARC KOTO experiment. With a single event sensitivity of $(7.1\pm 0.3_{\rm stat.} \pm 1.6_{\rm syst.})\times 10^{-8}$, no candidate event was observed in the signal region. The upper limit on the branching fraction was set to be $1.7\times 10^{-7}$ at the 90\% confidence level.


title: Seeding Method for Ice Nucleation under Shear

abstract: Hydrodynamic flow can have complex and far-reaching consequences on the rate of homogenous nucleation. We present a general formalism for calculating the nucleation rates of simply sheared systems. We have derived an extension to the conventional Classical Nucleation Theory, explicitly embodying the shear rate. Seeded Molecular Dynamics simulations form the backbone of our approach. The framework can be used for moderate supercoolings, at which temperatures brute-force methods are practically infeasible. The competing energetic and kinetic effects of shear arise naturally from the equations. We show how the theory can be used to identify shear regimes of ice nucleation behaviour for the mW water model, unifying disparate trends reported in the literature. At each temperature, we define a crossover shear rate in the limit of $1000-10,000 \ s^{-1}$, beyond which the nucleation rate increases steadily upto a maximum, at the optimal shear rate. For $235$, $240$, $255$ and $260 \ K$, the optimal shear rates are in the range of $\approx 10^6-10^7 \ s^{-1}$. For very high shear rates beyond $10^8 \ s^{-1}$, nucleation is strongly inhibited. Our results indicate that the shear-dependent nucleation rate curves have a non-monotonic dependence on temperature.


title: Interstellar extinction toward symbiotic stars

abstract: Using diffuse interstellar bands (DIBs) at 5780 A, 5797 A and 6613 A, visible in the high resolution spectra, and measuring their equivalent widths, we estimate the interstellar extinction toward seven symbiotic stars. We find $E_{B-V}$= $1.28 \pm 0.10$ for AS 289, $E_{B-V}$= $1.55 \pm 0.10$ for BI Cru, $E_{B-V}$= $0.63 \pm 0.10$ for HD330036, $E_{B-V}$= $0.33 \pm 0.05$ for V2756 Sgr, $E_{B-V}$= $0.30 \pm 0.05$ for V2905 Sgr, $E_{B-V}$= $1.52 \pm 0.11$ for V417 Cen, $E_{B-V}$= $0.81 \pm 0.10$ for PN Sa 3-22. The derived values are in agreement with the extinction through the Galaxy.


title: Scalable Method for Linear Optimization of Industrial Processes

abstract: In the development of industrial digital twins, the optimization problem of technological and business processes often arises. In many cases, this problem can be reduced to a large-scale linear programming (LP) problem. The article is devoted to the new method for solving large-scale LP problems. This method is called the "apex-method". The apex-method uses the predictor-corrector framework. The predictor step calculates a point belonging to the feasible region of LP problem. The corrector step calculates a sequence of points converging to the exact solution of the LP problem. The article gives a formal description of the apex-method and provides information about its parallel implementation in C++ language by using the MPI library. The results of large-scale computational experiments on a cluster computing system to study the scalability of the apex method are presented.


title: The $Z_5$ model of two-component dark matter

abstract: Scenarios for multi-component scalar dark matter based on a single $Z_N$ ($N\geq 4$) symmetry are simple and well-motivated. In this paper we investigate, for the first time, the phenomenology of the $Z_5$ model for two-component dark matter. This model, which can be seen as an extension of the well-known singlet scalar model, features two complex scalar fields--the dark matter particles--that are Standard Model singlets but have different charges under a $Z_5$ symmetry. The interactions allowed by the $Z_5$ give rise to novel processes between the dark matter particles that affect their relic densities and their detection prospects, which we study in detail. The key parameters of the model are identified and its viable regions are characterized by means of random scans. We show that, unlike the singlet scalar model, dark matter masses below the TeV are still compatible with present data. Even though the dark matter density turns out to be dominated by the lighter component, we find that current and future direct detection experiments may be sensitive to signals from both dark matter particles.


title: Approximating Euclidean by Imprecise Markov Decision Processes

abstract: Euclidean Markov decision processes are a powerful tool for modeling control problems under uncertainty over continuous domains. Finite state imprecise, Markov decision processes can be used to approximate the behavior of these infinite models. In this paper we address two questions: first, we investigate what kind of approximation guarantees are obtained when the Euclidean process is approximated by finite state approximations induced by increasingly fine partitions of the continuous state space. We show that for cost functions over finite time horizons the approximations become arbitrarily precise. Second, we use imprecise Markov decision process approximations as a tool to analyse and validate cost functions and strategies obtained by reinforcement learning. We find that, on the one hand, our new theoretical results validate basic design choices of a previously proposed reinforcement learning approach. On the other hand, the imprecise Markov decision process approximations reveal some inaccuracies in the learned cost functions.


title: From Newton's second law to Euler's equations of perfect fluids

abstract: Vlasov equations can be formally derived from N-body dynamics in the mean-field limit. In some suitable singular limits, they may themselves converge to fluid dynamics equations. Motivated by this heuristic, we introduce natural scalings under which the incompressible Euler equations can be rigorously derived from N-body dynamics with repulsive Coulomb interaction. Our analysis is based on the modulated energy methods of Brenier and Serfaty.


title: Does the $\ell_1$-norm Learn a Sparse Graph under Laplacian Constrained Graphical Models?

abstract: We consider the problem of learning a sparse graph under the Laplacian constrained Gaussian graphical models. This problem can be formulated as a penalized maximum likelihood estimation of the Laplacian constrained precision matrix. Like in the classical graphical lasso problem, recent works made use of the $\ell_1$-norm regularization with the goal of promoting sparsity in Laplacian constrained precision matrix estimation. However, we find that the widely used $\ell_1$-norm is not effective in imposing a sparse solution in this problem. Through empirical evidence, we observe that the number of nonzero graph weights grows with the increase of the regularization parameter. From a theoretical perspective, we prove that a large regularization parameter will surprisingly lead to a complete graph, i.e., every pair of vertices is connected by an edge. To address this issue, we introduce the nonconvex sparsity penalty, and propose a new estimator by solving a sequence of weighted $\ell_1$-norm penalized sub-problems. We establish the non-asymptotic optimization performance guarantees on both optimization error and statistical error, and prove that the proposed estimator can recover the edges correctly with a high probability. To solve each sub-problem, we develop a projected gradient descent algorithm which enjoys a linear convergence rate. Finally, an extension to learn disconnected graphs is proposed by imposing additional rank constraint. We propose a numerical algorithm based on based on the alternating direction method of multipliers, and establish its theoretical sequence convergence. Numerical experiments involving synthetic and real-world data sets demonstrate the effectiveness of the proposed method.


title: Classification of one dimensional dynamical systems by countable structures

abstract: We study the complexity of the classification problem of conjugacy on dynamical systems on some compact metrizable spaces. Especially we prove that the conjugacy equivalence relation of interval dynamical systems is Borel bireducible to isomorphism equivalence relation of countable graphs. This solves a special case of the Hjorth's conjecture which states that every orbit equivalence relation induced by a continuous action of the group of all homeomorphisms of the closed unit interval is classifiable by countable structures. We also prove that conjugacy equivalence relation of Hilbert cube homeomorphisms is Borel bireducible to the universal orbit equivalence relation.


title: Optimizing Cherenkov photons generation and propagation in CORSIKA for CTA Monte-Carlo simulations

abstract: COsmic Ray SImulations for KAscade) is a program for detailed simulation of extensive air showers initiated by high energy cosmic ray particles in the atmosphere, and is used today by almost all the major instruments that aim at measuring primary and secondary cosmic rays on the ground. The Cherenkov Telescope Array (CTA), currently under construction, is the next-generation instrument in the field of very-high-energy gamma-ray astronomy. Detailed CORSIKA Monte Carlo simulations will be regularly performed in parallel to CTA operations to estimate the instrument response functions, necessary to extract the physical properties of the cosmic sources from the measurements during data analysis. The estimated CPU time associated with these simulations is very high, of the order of 200 million HS06 hours per year. Code optimization becomes a necessity towards fast productions and limited costs. We propose in this paper multiple code transformations that aim to facilitate automatic vectorization done by the compiler, ensuring minimal external libraries requirement and high hardware portability.


title: Interacting tachyonic scalar field

abstract: We discuss the coupling between dark energy and matter by considering a homogeneous tachyonic scalar field as a candidate for dark energy. We obtained the functional form of scale factor by assuming that the coupling strength depends linearly on the Hubble parameter and energy density. We also estimated the cosmic age of the universe for different values of coupling constant.


title: On clusters and the multi-isoperimetric profile in Riemannian manifolds with bounded geometry

abstract: For a complete Riemannian manifold with bounded geometry, we prove the existence of isoperimetric clusters and also the compactness theorem for sequence of clusters in a larger space obtained by adding finitely many limit manifolds at infinity. Moreover, we show that isoperimetric clusters are bounded. We introduce and prove the Holder continuity of the multi-isoperimetric profile which has been explored by Emanuel Milman and Joe Neeman with a Gaussian-weighted notion of perimeter. We yield a proof of classical existence theorem, e.g. in space forms, for isoperimetric cluster using the results presented here. The results in this work generalize previous works of Stefano Nardulli, Andrea Mondino, Frank Morgan, Matteo Galli and Manuel Ritor\'e from the classical Riemannian and sub-Riemannian isoperimetric problem to the context of Riemannian isoperimetric clusters and also Frank Morgan and Francesco Maggi works on the clusters theory in the Euclidean setting.


title: A metabolomic measure of energy metabolism moderates how an inflammatory miRNA relates to rs-fMRI network and motor control in football athletes

abstract: Collision sports athletes experience many head acceleration events (HAEs) per season. The effects of these subconcussive events are largely understudied since HAEs may produce no overt symptoms, and are likely to diffusely manifest across multiple scales of study (e.g., molecular, cellular network, and behavior). This study integrated resting-state fMRI with metabolome, transcriptome and computational virtual reality (VR) behavior measures to assess the effects of exposure to HAEs on players in a collegiate American football team. Permutation-based mediation and moderation analysis was used to investigate relationships between network fingerprint, changes in omic measures and VR metrics over the season. Change in an energy cycle fatty acid, tridecenedioate, moderated the relationship between 1) miR-505 and DMN fingerprint and 2) the relationship between DMN fingerprint and worsening VR Balance measures (all p less than or equal to 0.05). In addition, the similarity in DMN over the season was negatively related to cumulative number of HAEs above 80G, and DMN fingerprint was less similar across the season in athletes relative to age-matched non-athletes. miR-505 was also positively related to average number of HAEs above 25G per session. It is important to note that tridecenedioate has a double bond making it a candidate for ROS scavenging. These findings between a candidate ROS-related metabolite, inflammatory miRNA, altered brain imaging and diminished behavioral performance suggests that impact athletes may experience chronic neuroinflammation. The rigorous permutation-based mediation/moderation may provide a methodology for investigating complex multi-scale biological data within humans alone and thus assist study of other functional brain problems.


title: An operator-theoretical study of the specific heat and the critical magnetic field in the BCS-Bogoliubov model of superconductivity

abstract: In the preceding paper, introducing a cutoff, the present author gave a proof of the statement that the transition to a superconducting state is a second-order phase transition in the BCS-Bogoliubov model of superconductivity on the basis of fixed-point theorems, and solved the long-standing problem of the second-order phase transition from the viewpoint of operator theory. In this paper we study the temperature dependence of the specific heat and the critical magnetic field in the model from the viewpoint of operator theory. We first show some properties of the solution to the BCS-Bogoliubov gap equation with respect to the temperature, and give the exact and explicit expression for the gap in the specific heat divided by the specific heat. We then show that it does not depend on superconductors and is a universal constant. Moreover, we show that the critical magnetic field is smooth with respect to the temperature, and point out the behavior of both the critical magnetic field and its derivative.

categories: math-ph cond-mat.stat-mech cond-mat.supr-con math.AP math.MP


title: Symmetries and Differential Invariants for Viscid Flows on a Curve

abstract: In this paper, flows of a viscid fluids on curves are considered. Symmetry algebras and the corresponding fields of differential invariants are found. We study their dependence on thermodynamic states of media, and provide classification of thermodynamic states.


title: Influence of Various Temporal Recoding on Pavlovian Eyeblink Conditioning in The Cerebellum

abstract: We consider the Pavlovian eyeblink conditioning (EBC) via repeated presentation of paired conditioned stimulus (tone) and unconditioned stimulus (airpuff). The influence of various temporal recoding of granule cells on the EBC is investigated in a cerebellar network where the connection probability $p_c$ from Golgi to granule cells is changed. In an optimal case of $p_c^*~(=0.029)$, individual granule cells show various well- and ill-matched firing patterns relative to the unconditioned stimulus. Then, these variously-recoded signals are fed into the Purkinje cells (PCs) through parallel-fibers (PFs). In the case of well-matched PF-PC synapses, their synaptic weights are strongly depressed through strong long-term depression (LTD). On the other hand, practically no LTD occurs for the ill-matched PF-PC synapses. This type of "effective" depression at the PF-PC synapses coordinates firings of PCs effectively, which then make effective inhibitory coordination on cerebellar nucleus neuron [which elicits conditioned response (CR; eyeblink)]. When the learning trial passes a threshold, acquisition of CR begins. In this case, the timing degree ${\cal T}_d$ of CR becomes good due to presence of the ill-matched firing group which plays a role of protection barrier for the timing. With further increase in the trial, strength $\cal S$ of CR (corresponding to the amplitude of eyelid closure) increases due to strong LTD in the well-matched firing group. Thus, with increasing the learning trial, the (overall) learning efficiency degree ${\cal L}_e$ (taking into consideration both timing and strength of CR) for the CR is increased, and eventually it becomes saturated. By changing $p_c$ from $p_c^*$, we also investigate the influence of various temporal recoding on the EBC. It is thus found that, the more various in temporal recoding, the more effective in learning for the Pavlovian EBC.


title: Cancellation theorem for motivic spaces with finite flat transfers

abstract: We show that the category of motivic spaces with transfers along finite flat morphisms, over a perfect field, satisfies all the properties we have come to expect of good categories of motives. In particular we establish the analog of Voevodsky's cancellation theorem.


title: Quantum ergodicity for Eisenstein series on hyperbolic surfaces of large genus

abstract: We give a quantitative estimate for the quantum mean absolute deviation on hyperbolic surfaces of finite area in terms of geometric parameters such as the genus, number of cusps and injectivity radius. It implies a delocalisation result of quantum ergodicity type for eigenfunctions of the Laplacian on hyperbolic surfaces of finite area that Benjamini-Schramm converge to the hyperbolic plane. We show that this is generic for Mirzakhani's model of random surfaces chosen uniformly with respect to the Weil-Petersson volume. Depending on the particular sequence of surfaces considered this gives a result of delocalisation of most cusp forms or of Eisenstein series.


title: Parameter estimation for strong phase transitions in supranuclear matter using gravitational-wave astronomy

abstract: At supranuclear densities, explored in the core of neutron stars, a strong phase transition from hadronic matter to more exotic forms of matter might be present. To test this hypothesis, binary neutron-star mergers offer a unique possibility to probe matter at densities that we can not create in any existing terrestrial experiment. In this work, we show that, if present, strong phase transitions can have a measurable imprint on the binary neutron-star coalescence and the emitted gravitational-wave signal. We construct a new parameterization of the supranuclear equation of state that allows us to test for the existence of a strong phase transition and extract its characteristic properties purely from the gravitational-wave signal of the inspiraling neutron stars. We test our approach using a Bayesian inference study simulating 600 signals with three different equations of state and find that for current gravitational-wave detector networks already twelve events might be sufficient to verify the presence of a strong phase transition. Finally, we use our methodology to analyze GW170817 and GW190425, but do not find any indication that a strong phase transition is present at densities probed during the inspiral.


title: Joints in Random Forests

abstract: Decision Trees (DTs) and Random Forests (RFs) are powerful discriminative learners and tools of central importance to the everyday machine learning practitioner and data scientist. Due to their discriminative nature, however, they lack principled methods to process inputs with missing features or to detect outliers, which requires pairing them with imputation techniques or a separate generative model. In this paper, we demonstrate that DTs and RFs can naturally be interpreted as generative models, by drawing a connection to Probabilistic Circuits, a prominent class of tractable probabilistic models. This reinterpretation equips them with a full joint distribution over the feature space and leads to Generative Decision Trees (GeDTs) and Generative Forests (GeFs), a family of novel hybrid generative-discriminative models. This family of models retains the overall characteristics of DTs and RFs while additionally being able to handle missing features by means of marginalisation. Under certain assumptions, frequently made for Bayes consistency results, we show that consistency in GeDTs and GeFs extend to any pattern of missing input features, if missing at random. Empirically, we show that our models often outperform common routines to treat missing data, such as K-nearest neighbour imputation, and moreover, that our models can naturally detect outliers by monitoring the marginal probability of input features.


title: Photon Collider for energies 1-2 TeV

abstract: We discuss a photon collider based on the $e^+e^-$ \ linear collider with energies of $2E = 1\div 2$ TeV in cms (ILC, CLIC, ...). Previously, this energy range was considered hopeless for the experiment in the foreseeable future. We discuss the realization of the TeV PLC based on modern lasers. A small modification of the laser-optical system prepared for the photon collider with energy of $E \approx 250$ GeV will be sufficient if the parameters are chosen optimally. The high-energy part of the photon spectrum does not depend on design details and is well separated from the low-energy part. That is a narrow band near the upper boundary, about 5\% wide. The high-energy integrated \ggam luminosity is about 1/5 and the maximum differential luminosity is about 1/4 of the corresponding values for the photon collider with $E \approx 250$ GeV.


title: LSBert: A Simple Framework for Lexical Simplification

abstract: Lexical simplification (LS) aims to replace complex words in a given sentence with their simpler alternatives of equivalent meaning, to simplify the sentence. Recently unsupervised lexical simplification approaches only rely on the complex word itself regardless of the given sentence to generate candidate substitutions, which will inevitably produce a large number of spurious candidates. In this paper, we propose a lexical simplification framework LSBert based on pretrained representation model Bert, that is capable of (1) making use of the wider context when both detecting the words in need of simplification and generating substitue candidates, and (2) taking five high-quality features into account for ranking candidates, including Bert prediction order, Bert-based language model, and the paraphrase database PPDB, in addition to the word frequency and word similarity commonly used in other LS methods. We show that our system outputs lexical simplifications that are grammatically correct and semantically appropriate, and obtains obvious improvement compared with these baselines, outperforming the state-of-the-art by 29.8 Accuracy points on three well-known benchmarks.


title: A new perturbative solution to the motion around triangular Lagrangian points in the elliptic restricted three-body problem

abstract: The equations of motion of planar elliptic restricted three body problem are transformed to four decoupled Hill's equations. By using the Floquet theorem analytic solution to the oscillator equations with time dependent periodic coefficients are presented. We show that the new analytic approach is valid for system parameters $0 < e \leq 0.05$ and $0 < \mu \leq 0.01$ where $e$ denotes the eccentricity of primaries while $\mu$ is the mass parameter, respectively. We also clarify the transformation details that provide the applicability of the method.


title: Streaming Transformer ASR with Blockwise Synchronous Beam Search

abstract: The Transformer self-attention network has shown promising performance as an alternative to recurrent neural networks in end-to-end (E2E) automatic speech recognition (ASR) systems. However, Transformer has a drawback in that the entire input sequence is required to compute both self-attention and source--target attention. In this paper, we propose a novel blockwise synchronous beam search algorithm based on blockwise processing of encoder to perform streaming E2E Transformer ASR. In the beam search, encoded feature blocks are synchronously aligned using a block boundary detection technique, where a reliability score of each predicted hypothesis is evaluated based on the end-of-sequence and repeated tokens in the hypothesis. Evaluations of the HKUST and AISHELL-1 Mandarin, LibriSpeech English, and CSJ Japanese tasks show that the proposed streaming Transformer algorithm outperforms conventional online approaches, including monotonic chunkwise attention (MoChA), especially when using the knowledge distillation technique. An ablation study indicates that our streaming approach contributes to reducing the response time, and the repetition criterion contributes significantly in certain tasks. Our streaming ASR models achieve comparable or superior performance to batch models and other streaming-based Transformer methods in all tasks considered.


title: ELMV: an Ensemble-Learning Approach for Analyzing Electrical Health Records with Significant Missing Values

abstract: Many real-world Electronic Health Record (EHR) data contains a large proportion of missing values. Leaving substantial portion of missing information unaddressed usually causes significant bias, which leads to invalid conclusion to be drawn. On the other hand, training a machine learning model with a much smaller nearly-complete subset can drastically impact the reliability and accuracy of model inference. Data imputation algorithms that attempt to replace missing data with meaningful values inevitably increase the variability of effect estimates with increased missingness, making it unreliable for hypothesis validation. We propose a novel Ensemble-Learning for Missing Value (ELMV) framework, which introduces an effective approach to construct multiple subsets of the original EHR data with a much lower missing rate, as well as mobilizing a dedicated support set for the ensemble learning in the purpose of reducing the bias caused by substantial missing values. ELMV has been evaluated on a real-world healthcare data for critical feature identification as well as a batch of simulation data with different missing rates for outcome prediction. On both experiments, ELMV clearly outperforms conventional missing value imputation methods and ensemble learning models.


title: Dynamics of an imprecise stochastic Holling II one-predator two-prey system with jumps

abstract: Groups in ecology are often affected by sudden environmental perturbations. Parameters of stochastic models are often imprecise due to various uncertainties. In this paper, we formulate a stochastic Holling II one-predator two-prey system with jumps and interval parameters. Firstly, we prove the existence and uniqueness of the positive solution. Moreover, the sufficient conditions for the extinction and persistence in the mean of the solution are obtained.


title: Enhanced four-wave-mixing and 3rd order optical nonlinearity in SiN nanowires integrated with graphene oxide films

abstract: Layered 2D graphene oxide (GO) films are integrated with silicon nitride (SiN) waveguides to experimentally demonstrate an enhanced Kerr nonlinearity via four wave mixing (FWM). Owing to the strong light matter interaction between the SiN waveguides and the highly nonlinear GO films, the FWM performance of the hybrid waveguides is significantly improved. SiN waveguides with both uniformly coated and patterned GO films are fabricated based on a transfer free, layer by layer GO coating method together with standard photolithography and lift off processes, yielding precise control of the film thickness, placement and coating length. Detailed FWM measurements are carried out for the fabricated devices with different numbers of GO layers and at different pump powers. By optimizing the trade off between the nonlinearity and loss, we obtain a significant improvement in the FWM conversion efficiency of 7.3 dB for a uniformly coated device with 1 layer of GO and 9.1 dB for a patterned device with 5 layers of GO. We also obtain a significant increase in FWM bandwidth for the patterned devices. A detailed analysis of the influence of pattern length and position on the FWM performance is performed. Based on the FWM measurements, the dependence of GOs third-order nonlinearity on layer number and pump power is also extracted, revealing interesting physical insights about the 2D layered GO films. Finally, we obtain an enhancement in the effective nonlinear parameter of the hybrid waveguides by over a factor of 100. These results verify the enhanced nonlinear optical performance of SiN waveguides achievable by incorporating 2D layered GO films.


title: An even better Density Increment Theorem and its application to Hadwiger's Conjecture

abstract: In 1943, Hadwiger conjectured that every graph with no $K_t$ minor is $(t-1)$-colorable for every $t\ge 1$. In the 1980s, Kostochka and Thomason independently proved that every graph with no $K_t$ minor has average degree $O(t\sqrt{\log t})$ and hence is $O(t\sqrt{\log t})$-colorable. Recently, Norin, Song and the author showed that every graph with no $K_t$ minor is $O(t(\log t)^{\beta})$-colorable for every $\beta > 1/4$, making the first improvement on the order of magnitude of the $O(t\sqrt{\log t})$ bound. More recently, the author showed that every graph with no $K_t$ minor is $O(t (\log t)^{\beta})$-colorable for every $\beta > 0$; more specifically, they are $t \cdot 2^{ O((\log \log t)^{2/3}) }$-colorable. In combination with that work, we show in this paper that every graph with no $K_t$ minor is $O(t (\log \log t)^{6})$-colorable.


title: Hadronic vs leptonic models for $\gamma$-ray emission from VER J2227+608

abstract: Recent observations of VER J2227+608 reveal a broken power $\gamma$-ray spectrum with the spectral index increasing from $\sim 1.8$ in the GeV energy range to $\sim 2.3$ in the TeV range. Such a spectral break can be attributed to radiative energy loss of energetic electrons in the leptonic scenario for the $\gamma$-ray emission, which, in combination with characteristic age of the nearby pulsar, can be used to constrain magnetic field in the emission region. We show that the radio and X-ray observations can also be explained in such a scenario. In the hadronic scenario, the spectral break can be attributed to diffusion of energetic ions in a turbulent medium and detailed spectral measurement can be used to constrain the diffusion coefficient. These two models, however, predict drastically different spectra above 100 TeV, which will be uncovered with future high-resolution observations, such as LHAASO and/or CTA.


title: Distributed Policy Synthesis of Multi-Agent Systems With Graph Temporal Logic Specifications

abstract: We study the distributed synthesis of policies for multi-agent systems to perform \emph{spatial-temporal} tasks. We formalize the synthesis problem as a \emph{factored} Markov decision process subject to \emph{graph temporal logic} specifications. The transition function and task of each agent are functions of the agent itself and its neighboring agents. In this work, we develop another distributed synthesis method, which improves the scalability and runtime by two orders of magnitude compared to our prior work. The synthesis method decomposes the problem into a set of smaller problems, one for each agent by leveraging the structure in the model, and the specifications. We show that the running time of the method is linear in the number of agents. The size of the problem for each agent is exponential only in the number of neighboring agents, which is typically much smaller than the number of agents. We demonstrate the applicability of the method in case studies on disease control, urban security, and search and rescue. The numerical examples show that the method scales to hundreds of agents with hundreds of states per agent and can also handle significantly larger state spaces than our prior work.


title: On properties of composites and monoid domains

abstract: In this paper I consider all possible properties from commutative algebra for polynomial composites and monoid domains. The aim is full characterization of these structures. I start with the examination of group, ring, modules properties, graded, but also the study of invertible elements, irreducible elements, ideals, etc. in these structures. In the second part of the work I give examples of the use of composites and monoid domains in cryptology. Each such polynomial is the sum of the products of the variable and the coefficient. And what if subsequent coefficient sets are appropriate cryptographic systems? Similarly, monoid domains can be a very good tool between encrypting and decrypting messages.


title: Stability of a star-shaped network with local Kelvin-Voigt damping and non-smooth coefficient at interface

abstract: In this paper, we study the stability problem of a star-shaped network of elastic strings with a local Kelvin-Voigt damping. Under the assumption that the damping coefficients have some singularities near the transmission point, we prove that the semigroup corresponding to the system is polynomially stable and the decay rates depends on the speed of the degeneracy. This result improves the decay rate of the semigroup associated to the system on an earlier result of Z.~Liu and Q.~Zhang in \cite{LZ} involving the wave equation with local Kelvin-Voigt damping and non-smooth coefficient at interface.


title: Relative Deviation Margin Bounds

abstract: We present a series of new and more favorable margin-based learning guarantees that depend on the empirical margin loss of a predictor. We give two types of learning bounds, both distribution-dependent and valid for general families, in terms of the Rademacher complexity or the empirical $\ell_\infty$ covering number of the hypothesis set used. Furthermore, using our relative deviation margin bounds, we derive distribution-dependent generalization bounds for unbounded loss functions under the assumption of a finite moment. We also briefly highlight several applications of these bounds and discuss their connection with existing results.


title: A Liouville theorem for an integral equation of the Ginzburg-Landau type

abstract: In this paper, we are concerned with a Liouville-type result of the nonlinear integral equation \begin{equation*} u(x)=\overrightarrow{l}+C_*\int_{\mathbb{R}^{n}}\frac{u(1-|u|^{2})}{|x-y|^{n-\alpha}}dy. \end{equation*} Here $u: \mathbb{R}^{n} \to \mathbb{R}^{k}$ is a bounded, uniformly continuous and differentiable function with $k \geq 1$ and $1<\alpha<n$, $\overrightarrow{l} \in \mathbb{R}^{k}$ is a constant vector, and $C_*$ is a real constant. If $u$ is the finite energy solution, we prove that $|\overrightarrow{l}| \in \{0,1\}$. Furthermore, we also give a Liouville type theorem (i.e., $u \equiv \overrightarrow{l}$).


title: A Hybrid Fast Multipole Method for Cosmological N-body Simulations

abstract: We investigate a hybrid numerical algorithm aimed at the large-scale cosmological N-body simulation for the on-going and the future high precious sky surveys. It makes use of a truncated Fast Multiple Method (FMM) for short-range gravity, incorporating with a Particle Mesh (PM) method for long-range potential, which is applied to deal with extremely large particle number. In this work, we present a specific strategy to modify a conventional FMM by a Gaussian shaped factor and provide quantitative expressions for the interaction kernels between multipole expansions. Moreover, a proper multipole acceptance criteria for the hybrid method is introduced to solve potential precision loss induced by the truncation. Such procedures reduce the mount of computation than an original FMM and decouple the global communication. A simplified version of code is introduced to verify the hybrid algorithm, accuracy and parallel implementation.


title: What they do when in doubt: a study of inductive biases in seq2seq learners

abstract: Sequence-to-sequence (seq2seq) learners are widely used, but we still have only limited knowledge about what inductive biases shape the way they generalize. We address that by investigating how popular seq2seq learners generalize in tasks that have high ambiguity in the training data. We use SCAN and three new tasks to study learners' preferences for memorization, arithmetic, hierarchical, and compositional reasoning. Further, we connect to Solomonoff's theory of induction and propose to use description length as a principled and sensitive measure of inductive biases. In our experimental study, we find that LSTM-based learners can learn to perform counting, addition, and multiplication by a constant from a single training example. Furthermore, Transformer and LSTM-based learners show a bias toward the hierarchical induction over the linear one, while CNN-based learners prefer the opposite. On the SCAN dataset, we find that CNN-based, and, to a lesser degree, Transformer- and LSTM-based learners have a preference for compositional generalization over memorization. Finally, across all our experiments, description length proved to be a sensitive measure of inductive biases.


title: Regular orbits of quasisimple linear groups II

abstract: Let $V$ be a finite-dimensional vector space over a finite field, and suppose $G \leq \Gamma \mathrm{L}(V)$ is a group with a unique subnormal quasisimple subgroup $E(G)$ that is absolutely irreducible on $V$. A base for $G$ is a set of vectors $B\subseteq V$ with pointwise stabiliser $G_B=1$. If $G$ has a base of size 1, we say that it has a regular orbit on $V$. In this paper we investigate the minimal base size of groups $G$ with $E(G)/Z(E(G)) \cong \mathrm{PSL}_n(q)$ in defining characteristic, with an aim of classifying those with a regular orbit on $V$.


title: Quantum oscillations with magnetic hysteresis observed in CeTe$_{3}$ thin films

abstract: We have performed magnetotransport measurements in CeTe$_{3}$ thin films down to $0.2~{\rm K}$. It is known that CeTe$_{3}$ has two magnetic transitions at $T_{\rm N1} \approx 3~{\rm K}$ and $T_{\rm N2} \approx 1~{\rm K}$. A clear Shubnikov-de-Haas (SdH) oscillation was observed at $4~{\rm K}$, demonstrating the strong two-dimensional nature in this material. Below $T_{\rm N2}$, the SdH oscillation has two frequencies, indicating that the Fermi surface could be slightly modulated due to the second magnetic transition. We also observed a magnetic hysteresis in the SdH oscillation below $T_{\rm N1}$. Especially, there is a unique spike in the magnetoresistance at $B \approx 0.6~{\rm T}$ only when the magnetic field is swept from a high enough field (more than $2~{\rm T}$) to zero field.


title: NECI: N-Electron Configuration Interaction with emphasis on state-of-the-art stochastic methods

abstract: We present NECI, a state-of-the-art implementation of the Full Configuration Interaction Quantum Monte Carlo algorithm, a method based on a stochastic application of the Hamiltonian matrix on a sparse sampling of the wave function. The program utilizes a very powerful parallelization and scales efficiently to more than 24000 CPU cores. In this paper, we describe the core functionalities of NECI and recent developments. This includes the capabilities to calculate ground and excited state energies, properties via the one- and two-body reduced density matrices, as well as spectral and Green's functions for ab initio and model systems. A number of enhancements of the bare FCIQMC algorithm are available within NECI, allowing to use a partially deterministic formulation of the algorithm, working in a spin-adapted basis or supporting transcorrelated Hamiltonians. NECI supports the FCIDUMP file format for integrals, supplying a convenient interface to numerous quantum chemistry programs and it is licensed under GPL-3.0.


title: Elimination of Redundant Polynomial Constraints and Its Use in Constrained Control

abstract: The reduction of constraints to obtain minimal representations of sets is a very common problem in many engineering applications. While well-established methodologies exist for the case of linear constraints, the problem of how to detect redundant non-linear constraints is an open problem. In this paper we present a novel methodology based on Sum of Squares for the elimination of redundant polynomial constraints. The paper also presents some relevant applications of the presented method to constrained control problems. In particular, we show how the proposed method can be used in the Model Predictive Control and in the Reference Governor frameworks to reduce the computational burden of the online algorithms. Furthermore, this method can also be used to eliminate the terminal constraints in MPC in a simple way that is independent from the cost function.


title: Spin and orbital Edelstein effect in an oxide two-dimensional electron gas: theory and application to AlO$_x$/SrTiO$_{3}$

abstract: The Edelstein effect provides the purely electrical generation and control of a homogeneous magnetization in primarily nonmagnetic materials with broken inversion symmetry. Usually, only the spin density response to an external electric field is discussed. Here, we report on the electrically induced magnetization containing spin as well as orbital contributions in the topological oxide two-dimensional electron gas at the interface between SrTiO$_3$ and AlO. We find that in this particular system the orbital Edelstein effect exceeds the spin Edelstein effect by more than one order of magnitude. The main reason are orbital moments of different magnitude in the Rashba-like-split band pairs.


title: Anisotropic conformal invariance of lightlike geodesics in pseudo-Finsler manifolds

abstract: In this paper, we prove that lightlike geodesics of a pseudo-Finsler manifold and its focal points are preserved up to reparametrization by anisotropic conformal changes, using the Chern connection and the anisotropic calculus and the fact that geodesics are critical points of the energy functional and Jacobi fields, the kernel of its index form. This result has applications to the study of Finsler spacetimes.


title: On $p$-Laplacian reaction-diffusion problems with dynamical boundary conditions in perforated media

abstract: This paper deals with the homogenization of the $p$-Laplacian reaction-diffusion problems in a domain containing periodically distributed holes of size $\varepsilon$, with a dynamical boundary condition of pure-reactive type. We generalize our previous results established in the case where the diffusion is modeled by the Laplacian operator, i.e., with $p=2$. We prove the convergence of the homogenization process to a nonlinear $p$-Laplacian reaction-diffusion equation defined on a unified domain without holes with zero Dirichlet boundary condition and with extra terms coming from the influence of the nonlinear dynamical boundary conditions.


title: A statistical standard siren measurement of the Hubble constant from the LIGO/Virgo gravitational wave compact object merger GW190814 and Dark Energy Survey galaxies

abstract: We present a measurement of the Hubble constant $H_0$ using the gravitational wave (GW) event GW190814, which resulted from the coalescence of a 23 $M_\odot$ black hole with a 2.6 $M_\odot$ compact object, as a standard siren. No compelling electromagnetic counterpart has been identified for this event, thus our analysis accounts for thousands of potential host galaxies within a statistical framework. The redshift information is obtained from the photometric redshift (photo-$z$) catalog from the Dark Energy Survey. The luminosity distance is provided by the LIGO/Virgo gravitational wave sky map. Since this GW event has the second-smallest localization volume after GW170817, GW190814 is likely to provide the best constraint on cosmology from a single standard siren without identifying an electromagnetic counterpart. Our analysis uses photo-$z$ probability distribution functions and corrects for photo-$z$ biases. We also reanalyze the binary-black hole GW170814 within this updated framework. We explore how our findings impact the $H_0$ constraints from GW170817, the only GW merger associated with a unique host galaxy. From a combination of GW190814, GW170814 and GW170817, our analysis yields $H_0 = 72.0^{+ 12}_{- 8.2 }~{\rm km~s^{-1}~Mpc^{-1}}$ (68\% Highest Density Interval, HDI) for a prior in $H_0$ uniform between $[20,140]~{\rm km~s^{-1}~Mpc^{-1}}$. The addition of GW190814 and GW170814 to GW170817 improves the 68\% HDI from GW170817 alone by $\sim 18\%$, showing how well-localized mergers without counterparts can provide a significant contribution to standard siren measurements, provided that a complete galaxy catalog is available at the location of the event.


title: A Multiresolution Mesoscale Approach for Microscale Hydrodynamics

abstract: A new class of multiscale scheme is presented for micro-hydrodynamic problems based on a dual representation of the fluid observables. The hybrid model is first tested against the classical flow between two parallel plates and then applied to a plug flow within a micrometer-sized striction and a shear flow within a microcavity. Both cases demonstrate the capability of the multiscale approach to reproduce the correct macroscopic hydrodynamics also in the presence of refined grids (one and two levels), while retaining the correct thermal fluctuations, embedded in the multiparticle collision method. This provides the first step toward a novel class of fully mesoscale hybrid approaches able to capture the physics of fluids at the micro- and nanoscales whenever a continuum representation of the fluid falls short of providing the complete physical information, due to a lack of resolution and thermal fluctuations.


title: An approximation to zeros of the Riemann zeta function using fractional calculus

abstract: In this document, as far as the authors know, an approximation to the zeros of the Riemann zeta function has been obtained for the first time using only derivatives of constant functions, which was possible only because a fractional iterative method was used. This iterative method, valid for one and several variables, uses the properties of fractional calculus, in particular the fact that the fractional derivatives of constants are not always zero, to find multiple zeros of a function using a single initial condition. This partly solves the intrinsic problem of iterative methods that if we want to find N zeros it is necessary to give N initial conditions. Consequently, the method is suitable for approximating nontrivial zeros of the Riemann zeta function when the absolute value of its imaginary part tends to infinity. The deduction of the iterative method is presented, some examples of its implementation, and finally 53 different values near to the zeros of the Riemann zeta function are shown.


title: Flow-Based Network Creation Games

abstract: Network Creation Games(NCGs) model the creation of decentralized communication networks like the Internet. In such games strategic agents corresponding to network nodes selfishly decide with whom to connect to optimize some objective function. Past research intensively analyzed models where the agents strive for a central position in the network. This models agents optimizing the network for low-latency applications like VoIP. However, with today's abundance of streaming services it is important to ensure that the created network can satisfy the increased bandwidth demand. To the best of our knowledge, this natural problem of the decentralized strategic creation of networks with sufficient bandwidth has not yet been studied. We introduce Flow-Based NCGs where the selfish agents focus on bandwidth instead of latency. In essence, budget-constrained agents create network links to maximize their minimum or average network flow value to all other network nodes. Equivalently, this can also be understood as agents who create links to increase their connectivity and thus also the robustness of the network. For this novel type of NCG we prove that pure Nash equilibria exist, we give a simple algorithm for computing optimal networks, we show that the Price of Stability is 1 and we prove an (almost) tight bound of 2 on the Price of Anarchy. Last but not least, we show that our models do not admit a potential function.


title: Fast Green's function method for ultrafast electron-boson dynamics

abstract: The interaction of electrons with quantized phonons and photons underlies the ultrafast dynamics of systems ranging from molecules to solids, and it gives rise to a plethora of physical phenomena experimentally accessible using time-resolved techniques. Green's function methods offer an invaluable interpretation tool since scattering mechanisms of growing complexity can be selectively incorporated in the theory. Currently, however, real-time Green's function simulations are either prohibitively expensive due to the cubic scaling with the propagation time or do neglect the feedback of electrons on the bosons, thus violating energy conservation. We put forward a computationally efficient Green's function scheme which overcomes both limitations. The numerical effort scales linearly with the propagation time while the simultaneous dressing of electrons and bosons guarantees the fulfillment of all fundamental conservation laws. We present a real-time study of the phonon-driven relaxation dynamics in an optically excited narrow band-gap insulator, highlighting the nonthermal behavior of the phononic degrees of freedom. Our formulation paves the way to first-principles simulations of electron-boson systems with unprecedented long propagation times.


title: Odd degree isolated points on $X_1(N)$ with rational $j$-invariant

abstract: Let $C$ be a curve defined over a number field $k$. We say a closed point $x\in C$ of degree $d$ is isolated if it does not belong to an infinite family of degree $d$ points parametrized by the projective line or a positive rank abelian subvariety of the curve's Jacobian. Building on work of Bourdon, Ejder, Liu, Odumodu, and Viray, we characterize elliptic curves with rational $j$-invariant which give rise to an isolated point of odd degree on $X_1(N)/\mathbb{Q}$ for some positive integer $N$.


title: UV astronomy with small satellites

abstract: Small satellite platforms with high performance avionics are becoming more affordable. So far, with a few exceptions, small satellites have been mainly dedicated to earth observation. However, astronomy is a fascinating field with a history of large missions and a future of promising large mission candidates. This prompts many questions; can the recent affordability of small satellites change the landscape of space astronomy? What are the potential applications and scientific topics of interest, where small satellites could be instrumental for astronomy? What are the requirements and objectives that need to be fulfilled to successfully address the astronomical investigations of interest? Which kind of instrumentation suits the small platforms and the scientific use cases best? This paper discusses possible scientific use cases that can be achievable with a relatively small telescope aperture of 36 cm, as an example. The result of this survey points to a specific niche market -astronomy observation in the UV spectral range. UV astronomy is a research field which has had valuable scientific impact. It is, however, not the focus of many current or past astronomical investigations. UV astronomy measurements cannot be made from earth, due to atmospheric absorption in this spectral range. The research field is currently sparsely addressed but of scientific interest for a large community. Small satellites offer the opportunity to provide more means of research for UV astronomy. Therefore, this paper also presents an instrument design with a modest telescope aperture, a spectrometer and a detector that is suitable for observations in the UV. The observatory design can be accommodated on small platforms for the selected scientific use cases. It fulfills the scientific objectives and requirements of those use cases.


title: The word problem of the Brin-Higman-Thompson groups

abstract: We show that the word problem of the Brin-Higman-Thompson group $n G_{k,1}$ is {\sf coNP}-complete for all $n \ge 2$ and all $k \ge 2$. For this we prove that $n G_{k,1}$ is finitely generated, and that $n G_{k,1}$ contains a subgroup of $2 G_{2,1}$ that can represent bijective circuits. We also show that for all $n \ge 1$ and $k \ge 2$: \ If $\,K = 1 + (k-1)\,N\,$ for some $N \ge 1$, then $n G_{K,1} \le n G_{k,1}$. In particular, $n G_{K,1} \le n G_{2,1}$ for all $K \ge 2$.


title: Fully Abstract and Robust Compilation and How to Reconcile the Two, Abstractly

abstract: The most prominent formal criterion for secure compilation is full abstraction, the preservation and reflection of contextual equivalence. Recent work introduced robust compilation, defined as the preservation of robust satisfaction of hyperproperties, i.e., their satisfaction against arbitrary attackers. In this paper, we initially set out to compare these two approaches to secure compilation. To that end, we provide an exact description of the hyperproperties that are robustly satisfied by programs compiled with a fully abstract compiler, and show that they can be meaningless or trivial. We then propose a novel criterion for secure compilation formulated in the framework of Mathematical Operational Semantics (MOS), guaranteeing both full abstraction and the preservation of robust satisfaction of hyperproperties in a more sensible manner.


title: Fast Multi-Level Foreground Estimation

abstract: Alpha matting aims to estimate the translucency of an object in a given image. The resulting alpha matte describes pixel-wise to what amount foreground and background colors contribute to the color of the composite image. While most methods in literature focus on estimating the alpha matte, the process of estimating the foreground colors given the input image and its alpha matte is often neglected, although foreground estimation is an essential part of many image editing workflows. In this work, we propose a novel method for foreground estimation given the alpha matte. We demonstrate that our fast multi-level approach yields results that are comparable with the state-of-the-art while outperforming those methods in computational runtime and memory usage.


title: Positivity--preserving numerical scheme for hyperbolic systems with $\delta\,-$ shock solutions and its convergence analysis

abstract: Godunov type numerical schemes for the class of hyperbolic systems, admitting non-classical $\delta-$ shocks are proposed. It is shown that the numerical approximations converge to the solution and preserve the physical properties of the system such as positive density and bounded velocity. The scheme has been extended to positivity preserving and velocity bound preserving second-order accurate scheme by using appropriate slope limiters. The numerical results are compared with the existing the literature and the scheme is shown to capture the solution efficiently. The paper presents a hyperbolic system, for which an entropy satisfying scheme is constructed through an appropriate decoupling of the system into two scalar conservation laws with discontinuous flux.


title: On 2-Clubs in Graph-Based Data Clustering: Theory and Algorithm Engineering

abstract: Editing a graph into a disjoint union of clusters is a standard optimization task in graph-based data clustering. Here, complementing classic work where the clusters shall be cliques, we focus on clusters that shall be 2-clubs, that is, subgraphs of diameter two. This naturally leads to the two NP-hard problems 2-Club Cluster Editing (the allowed editing operations are edge insertion and edge deletion) and 2-Club Cluster Vertex Deletion (the allowed editing operations are vertex deletions). Answering an open question from the literature, we show that 2-Club Cluster Editing is W[2]-hard with respect to the number of edge modifications, thus contrasting the fixed-parameter tractability result for the classic Cluster Editing problem (considering cliques instead of 2-clubs). Then focusing on 2-Club Cluster Vertex Deletion, which is easily seen to be fixed-parameter tractable, we show that under standard complexity-theoretic assumptions it does not have a polynomial-size problem kernel when parameterized by the number of vertex deletions. Nevertheless, we develop several effective data reduction and pruning rules, resulting in a competitive solver, clearly outperforming a standard CPLEX solver in most instances of an established biological test data set.


title: Fluid drag reduction by magnetic confinement

abstract: The frictional forces of a viscous liquid flow are a major energy loss issue and severely limit microfluidics practical use. Reducing this drag by more than a few tens of percent remain illusive. Here, we show how cylindrical liquid-in-liquid flow leads to drag reduction of 60-99% for sub mm and mm sized channels, irrespective of whether the viscosity of the transported liquid is larger or smaller than that of the encapsulating one. In contrast to lubrication or sheath flow, we do not require the continuous flow of the encapsulating lubricant, here made up of a ferrofluid held in place by magnetic forces. In a laminar flow model with appropriate boundary conditions, we introduce a modified Reynolds number with a scaling that depends on geometrical factors and viscosity ratio of the two liquids. It explains our whole range of data and reveal the key design parameters for optimizing the drag reduction values. Our results therefore open the route to microfluidics designs with pressure gradients possibly reduced by orders of magnitudes.


title: Extracting Non-Gaussian Governing Laws from Data on Mean Exit Time

abstract: Motivated by the existing difficulties in establishing mathematical models and in observing the system state time series for some complex systems, especially for those driven by non-Gaussian Levy motion, we devise a method for extracting non-Gaussian governing laws with observations only on mean exit time. It is feasible to observe mean exit time for certain complex systems. With the observations, a sparse regression technique in the least squares sense is utilized to obtain the approximated function expression of mean exit time. Then, we learn the generator and further identify the stochastic differential equations through solving an inverse problem for a nonlocal partial differential equation and minimizing an error objective function. Finally, we verify the efficacy of the proposed method by three examples with the aid of the simulated data from the original systems. Results show that the method can apply to not only the stochastic dynamical systems driven by Gaussian Brownian motion but also those driven by non-Gaussian Levy motion, including those systems with complex rational drift.


title: Signatures of ubiquitous magnetic reconnection in the lower solar atmosphere

abstract: Ellerman Bomb-like brightenings of the hydrogen Balmer line wings in the quiet Sun (QSEBs) are a signature of the fundamental process of magnetic reconnection at the smallest observable scale in the solar lower atmosphere. We analyze high spatial resolution observations (0.1 arcsec) obtained with the Swedish 1-m Solar Telescope to explore signatures of QSEBs in the H$\beta$ line. We find that QSEBs are ubiquitous and uniformly distributed throughout the quiet Sun, predominantly occurring in intergranular lanes. We find up to 120 QSEBs in the FOV for a single moment in time; this is more than an order of magnitude higher than the number of QSEBs found in earlier H$\alpha$ observations. This suggests that about half a million QSEBs could be present in the lower solar atmosphere at any given time. The QSEB brightening found in the H$\beta$ line wings also persist in the line core with a temporal delay and spatial offset towards the nearest solar limb. Our results suggest that QSEBs emanate through magnetic reconnection along vertically extended current sheets in the solar lower atmosphere. The apparent omnipresence of small-scale magnetic reconnection may play an important role in the energy balance of the solar chromosphere.


title: Kink-driven magnetic reconnection in relativistic jets: consequences for X-ray polarimetry of BL Lacs

abstract: We investigate with relativistic MHD simulations the dissipation physics of BL Lac jets, by studying the synchrotron polarization signatures of particles accelerated by the kink instability in a magnetically-dominated plasma column. The nonlinear stage of the kink instability generates current sheets, where particles can be efficiently accelerated via magnetic reconnection. We identify current sheets as regions where s = J d/B is above some predefined threshold (where B is the field strength, J the current density and d the grid scale), and assume that the particle injection efficiency scales as proportional to the square of the current. X-ray emitting particles have short cooling times, so they only probe the field geometry of their injection sites. In contrast, particles emitting in the optical band, which we follow self-consistently as they propagate away from their injection sites while cooling, sample a larger volume, and so they may be expected to produce different polarimetric signatures. We find that the degree of polarization is roughly the same between X-ray and optical bands, because even the optical-emitting particles do not travel far from the current sheet where they were injected, due to lack of sufficient kink-generated turbulence. The polarization angle shows a different temporal evolution between the two bands, due to the different regions probed by X-ray and optical emitting particles. In view of the upcoming IXPE satellite, our results can help constrain whether kink-induced reconnection (as opposed to shocks) can be the source of multi-wavelength emission from BL Lacs.


title: Non-exponential decay of a collective excitation in an atomic ensemble coupled to a one-dimensional waveguide

abstract: We study the dynamics of a single excitation coherently shared amongst an ensemble of atoms and coupled to a one-dimensional wave guide. The coupling between the matter and the light field gives rise to collective phenomena such as superradiant states with an enhanced initial decay rate, but also to the coherent exchange of the excitation between the atoms. We find that the competition between the two phenomena provides a characteristic dynamics for the decay of the excitations, and remarkably exhibits an algebraic behavior, instead of the expected standard exponential one, for a large number of atoms. The analysis is first performed for a chiral waveguide, where the problem can be solved analytically, and then is extended to the bidirectional waveguide.


title: Online 3D Bin Packing with Constrained Deep Reinforcement Learning

abstract: We solve a challenging yet practically useful variant of 3D Bin Packing Problem (3D-BPP). In our problem, the agent has limited information about the items to be packed into the bin, and an item must be packed immediately after its arrival without buffering or readjusting. The item's placement also subjects to the constraints of collision avoidance and physical stability. We formulate this online 3D-BPP as a constrained Markov decision process. To solve the problem, we propose an effective and easy-to-implement constrained deep reinforcement learning (DRL) method under the actor-critic framework. In particular, we introduce a feasibility predictor to predict the feasibility mask for the placement actions and use it to modulate the action probabilities output by the actor during training. Such supervisions and transformations to DRL facilitate the agent to learn feasible policies efficiently. Our method can also be generalized e.g., with the ability to handle lookahead or items with different orientations. We have conducted extensive evaluation showing that the learned policy significantly outperforms the state-of-the-art methods. A user study suggests that our method attains a human-level performance.


title: The Sci-hub Effect: Sci-hub downloads lead to more article citations

abstract: Citations are often used as a metric of the impact of scientific publications. Here, we examine how the number of downloads from Sci-hub as well as various characteristics of publications and their authors predicts future citations. Using data from 12 leading journals in economics, consumer research, neuroscience, and multidisciplinary research, we found that articles downloaded from Sci-hub were cited 1.72 times more than papers not downloaded from Sci-hub and that the number of downloads from Sci-hub was a robust predictor of future citations. Among other characteristics of publications, the number of figures in a manuscript consistently predicts its future citations. The results suggest that limited access to publications may limit some scientific research from achieving its full impact.


title: Adaptive regularisation for ensemble Kalman inversion

abstract: We propose a new regularisation strategy for the classical ensemble Kalman inversion (EKI) framework. The strategy consists of: (i) an adaptive choice for the regularisation parameter in the update formula in EKI, and (ii) criteria for the early stopping of the scheme. In contrast to existing approaches, our parameter choice does not rely on additional tuning parameters which often have severe effects on the efficiency of EKI. We motivate our approach using the interpretation of EKI as a Gaussian approximation in the Bayesian tempering setting for inverse problems. We show that our parameter choice controls the symmetrised Kulback-Leibler divergence between consecutive tempering measures. We further motivate our choice using a heuristic statistical discrepancy principle. We test our framework using electrical impedance tomography with the complete electrode model. Parameterisations of the unknown conductivity are employed which enable us to characterise both smooth or a discontinuous (piecewise-constant) fields. We show numerically that the proposed regularisation of EKI can produce efficient, robust and accurate estimates, even for the discontinuous case which tends to require larger ensembles and more iterations to converge. We compare the proposed technique with a standard method of choice and demonstrate that the proposed method is a viable choice to address computational efficiency of EKI in practical/operational settings.


title: Crepant resolutions of double covers: On the Cynk-Hulek criterion for crepant resolutions of double cover

abstract: A collection $S = \{D_1,\ldots, D_n\}$ of divisors in a smooth variety $X$ is an {\em arrangement} if intersections of all subsets of $S$ are smooth. We show that a double cover of $X$ ramified on an arrangement has a crepant resolution under additional hypotheses. Namely, we assume that all intersection components that change the canonical divisor when blown up satisfy are {\em splayed}, a property of the tangent spaces of the components first studied by Faber. This strengthens a result of Cynk and Hulek, which requires a stronger hypothesis on the intersection components. Further, we study the singular subscheme of the union of the divisors in $S$ and prove that it has a primary decomposition where the primary components are supported on exactly the subvarieties which are blown up in the course of constructing the crepant resolution of the double cover.


title: Performance of ALICE AD modules in the CERN PS test beam

abstract: Two modules of the AD detector have been studied with the test beam at the T10 facility at CERN. The AD detector is made of scintillator pads read out by wave-length shifters (WLS)coupled to clean fibres that carry the produced light to photo-multiplier tubes (PMTs). In ALICE the AD is used to trigger and study the physics of diffractive and ultra-peripheral collisions as well as for a variety of technical tasks like beam-gas background monitoring or as a luminometer. The position dependence of the modules' efficiency has been measured and the effect of hits on the WLS or PMTs has been evaluated. The charge deposited by pions and protons has been measured at different momenta of the test beam. The time resolution is determined as a function of the deposited charge. These results are important ingredients to better understand the AD detector, to benchmark the corresponding simulations, and very importantly they served as a baseline for a similar device, the Forward Diffractive Detector (FDD), being currently built and that will be in operation in ALICE during the LHC Runs 3 and 4.


title: Solution of matching equations of IDA-PBC by Pfaffian differential equations

abstract: Finding the general solution of partial differential equations (PDEs) is essential for controller design in newly developed methods. Interconnection and damping assignment passivity based control (IDA-PBC) is one of such methods in which the solution to corresponding PDEs which are called matching equations, is needed to apply it in practice. In this paper, these matching equations are transformed to corresponding Pfaffian differential equations. Furthermore, it is shown that upon satisfaction of the integrability condition, the solution to the corresponding third-order Pfaffian differential equation may be obtained quite easily. The method is applied to the PDEs of IDA-PBC in some benchmark systems such as Magnetic levitation system, Pendubot, and underactuated cable driven robot to verify its applicability.


title: Suggestive Annotation of Brain Tumour Images with Gradient-guided Sampling

abstract: Machine learning has been widely adopted for medical image analysis in recent years given its promising performance in image segmentation and classification tasks. As a data-driven science, the success of machine learning, in particular supervised learning, largely depends on the availability of manually annotated datasets. For medical imaging applications, such annotated datasets are not easy to acquire. It takes a substantial amount of time and resource to curate an annotated medical image set. In this paper, we propose an efficient annotation framework for brain tumour images that is able to suggest informative sample images for human experts to annotate. Our experiments show that training a segmentation model with only 19% suggestively annotated patient scans from BraTS 2019 dataset can achieve a comparable performance to training a model on the full dataset for whole tumour segmentation task. It demonstrates a promising way to save manual annotation cost and improve data efficiency in medical imaging applications.


title: A Fokker-Planck approach to the study of robustness in gene expression

abstract: We study several Fokker-Planck equations arising from a stochastic chemical kinetic system modeling a gene regulatory network in biology. The densities solving the Fokker-Planck equations describe the joint distribution of the messenger RNA and micro RNA content in a cell. We provide theoretical and numerical evidences that the robustness of the gene expression is increased in the presence of micro RNA. At the mathematical level, increased robustness shows in a smaller coefficient of variation of the marginal density of the messenger RNA in the presence of micro RNA. These results follow from explicit formulas for solutions. Moreover, thanks to dimensional analyses and numerical simulations we provide qualitative insight into the role of each parameter in the model. As the increase of gene expression level comes from the underlying stochasticity in the models, we eventually discuss the choice of noise in our models and its influence on our results.


title: Classification of torus bundles that bound rational homology circles

abstract: In this article, we completely classify torus bundles over the circle that bound 4-manifolds with the rational homology of the circle. Along the way, we classify certain integral surgeries along chain links that bound rational homology balls and explore a connection to 3-braid closures whose double branched covers bound rational homology 4-balls.


title: Relaxed regularization for linear inverse problems

abstract: We consider regularized least-squares problems of the form $\min_{x} \frac{1}{2}\Vert Ax - b\Vert_2^2 + \mathcal{R}(Lx)$. Recently, Zheng et al., 2019, proposed an algorithm called Sparse Relaxed Regularized Regression (SR3) that employs a splitting strategy by introducing an auxiliary variable $y$ and solves $\min_{x,y} \frac{1}{2}\Vert Ax - b\Vert_2^2 + \frac{\kappa}{2}\Vert Lx - y\Vert_2^2 + \mathcal{R}(x)$. By minimizing out the variable $x$ we obtain an equivalent system $\min_{y} \frac{1}{2} \Vert F_{\kappa}y - g_{\kappa}\Vert_2^2+\mathcal{R}(y)$. In our work we view the SR3 method as a way to approximately solve the regularized problem. We analyze the conditioning of the relaxed problem in general and give an expression for the SVD of $F_{\kappa}$ as a function of $\kappa$. Furthermore, we relate the Pareto curve of the original problem to the relaxed problem and we quantify the error incurred by relaxation in terms of $\kappa$. Finally, we propose an efficient iterative method for solving the relaxed problem with inexact inner iterations. Numerical examples illustrate the approach.


title: Unlabelled Data Improves Bayesian Uncertainty Calibration under Covariate Shift

abstract: Modern neural networks have proven to be powerful function approximators, providing state-of-the-art performance in a multitude of applications. They however fall short in their ability to quantify confidence in their predictions - this is crucial in high-stakes applications that involve critical decision-making. Bayesian neural networks (BNNs) aim at solving this problem by placing a prior distribution over the network's parameters, thereby inducing a posterior distribution that encapsulates predictive uncertainty. While existing variants of BNNs based on Monte Carlo dropout produce reliable (albeit approximate) uncertainty estimates over in-distribution data, they tend to exhibit over-confidence in predictions made on target data whose feature distribution differs from the training data, i.e., the covariate shift setup. In this paper, we develop an approximate Bayesian inference scheme based on posterior regularisation, wherein unlabelled target data are used as "pseudo-labels" of model confidence that are used to regularise the model's loss on labelled source data. We show that this approach significantly improves the accuracy of uncertainty quantification on covariate-shifted data sets, with minimal modification to the underlying model architecture. We demonstrate the utility of our method in the context of transferring prognostic models of prostate cancer across globally diverse populations.


title: Tensor estimation with structured priors

abstract: We consider rank-one symmetric tensor estimation when the tensor is corrupted by Gaussian noise and the spike forming the tensor is a structured signal coming from a generalized linear model. The latter is a mathematically tractable model of a non-trivial hidden lower-dimensional latent structure in a signal. We work in a large dimensional regime with fixed ratio of signal-to-latent space dimensions. Remarkably, in this asymptotic regime, the mutual information between the spike and the observations can be expressed as a finite-dimensional variational problem, and it is possible to deduce the minimum-mean-square-error from its solution. We discuss, on examples, properties of the phase transitions as a function of the signal-to-noise ratio. Typically, the critical signal-to-noise ratio decreases with increasing signal-to-latent space dimensions. We discuss the limit of vanishing ratio of signal-to-latent space dimensions and determine the limiting tensor estimation problem. We also point out similarities and differences with the case of matrices.