title: A Multiplicative Property for Zero-Sums II

abstract: Let $G=C_n\oplus C_{mn}$ with $n\geq 2$ and $m\geq 1$, and let $k\in [0,n-1]$. It is known that any sequence of $mn+n-1+k$ terms from $G$ must contain a nontrivial zero-sum of length at most $mn+n-1-k$. The associated inverse question is to characterize those sequences with maximal length $mn+n-2+k$ that fail to contain a nontrivial zero-sum subsequence of length at most $mn+n-1-k$. For $k\leq 1$, this is the inverse question for the Davenport Constant. For $k=n-1$, this is the inverse question for the $\eta(G)$ invariant concerning short zero-sum subsequences. The structure in both these cases is known, and the structure for $k\in [2,n-2]$ when $m=1$ was studied previously with it conjectured that they must have the form $S=e_1^{[n-1]}\boldsymbol{\cdot} e_2^{[n-1]}\boldsymbol{\cdot} (e_1+e_2)^{[k]}$ for some basis $(e_1,e_2)$, with the conjecture established in many cases. We focus on $m\geq 2$. Assuming the conjectured structure holds for $k\in [2,n-2]$ in $C_n\oplus C_n$, we characterize the structure of all sequences of maximal length $mn+n-2+k$ in $C_n\oplus C_{mn}$ that fail to contain a nontrivial zero-sum of length at most $mn+n-1-k$, showing they must have either have the form $S=e_1^{[n-1]}\boldsymbol{\cdot} e_2^{[sn-1]}\boldsymbol{\cdot} (e_1+e_2)^{[(m-s)n+k]}$ for some $s\in [1,m]$ and basis $(e_1,e_2)$ with $\mathsf{ord}(e_2)=mn$, or else have the form $S=g_1^{[n-1]}\boldsymbol{\cdot} g_2^{[n-1]}\boldsymbol{\cdot} (g_1+g_2)^{[(m-1)n+k]}$ for some generating set $\{g_1,g_2\}$ with $\mathsf{ord}(g_1+g_2)=mn$. Additionally, we give a new proof of the precise structure in the case $k=n-1$ for $m=1$. Combined with known results, our results unconditionally establish the structure of extremal sequences in $G=C_n\oplus C_{mn}$ in many cases.


title: Entanglement witnessing by arbitrarily many independent observers recycling a local quantum shared state

abstract: We investigate the scenario where an observer, Alice, shares a two-qubit state with an arbitrary number of observers, Bobs, via sequentially and independently recycling the qubit in possession of the first Bob. It is known that there exist entangled states which can be used to have an arbitrarily long sequence of Bobs who can violate the Clauser-Horne-Shimony-Holt (CHSH) Bell inequality with the single Alice. We show that there exist entangled states that do not violate the Bell inequality and whose entanglement can be detected by an arbitrary number of Bobs by suitably choosing the entanglement witness operator and the unsharp measurement settings by the Bobs. This proves that the set of states that can be used to witness entanglement sequentially is larger than those that can witness sequential violation of local realism. There exist, therefore, two-party quantum correlations that are Bell "classical", but whose entanglement "nonclassicality" can be witnessed sequentially and independently by an arbitrarily large number of observers at one end of the shared state with the single observer at the other end.


title: Limit Cycles Bifurcating from a Periodic Annulus in Discontinuous Planar Piecewise Linear Hamiltonian differential System with Three Zones

abstract: In this paper, we study the number of limit cycles that can bifurcating from a periodic annulus in discontinuous planar piecewise linear Hamiltonian differential system with three zones separated by two parallel straight lines. We prove that if the central subsystem, i.e. the system defined between the two parallel lines, has a real center and the others subsystems have centers or saddles, then we have at least three limit cycles that appear after perturbations of periodic annulus. For this, we study the number of zeros of a Melnikov function for piecewise Hamiltonian system and present a normal form for this system in order to simplify the computations.


title: Example-Driven Model-Based Reinforcement Learning for Solving Long-Horizon Visuomotor Tasks

abstract: In this paper, we study the problem of learning a repertoire of low-level skills from raw images that can be sequenced to complete long-horizon visuomotor tasks. Reinforcement learning (RL) is a promising approach for acquiring short-horizon skills autonomously. However, the focus of RL algorithms has largely been on the success of those individual skills, more so than learning and grounding a large repertoire of skills that can be sequenced to complete extended multi-stage tasks. The latter demands robustness and persistence, as errors in skills can compound over time, and may require the robot to have a number of primitive skills in its repertoire, rather than just one. To this end, we introduce EMBER, a model-based RL method for learning primitive skills that are suitable for completing long-horizon visuomotor tasks. EMBER learns and plans using a learned model, critic, and success classifier, where the success classifier serves both as a reward function for RL and as a grounding mechanism to continuously detect if the robot should retry a skill when unsuccessful or under perturbations. Further, the learned model is task-agnostic and trained using data from all skills, enabling the robot to efficiently learn a number of distinct primitives. These visuomotor primitive skills and their associated pre- and post-conditions can then be directly combined with off-the-shelf symbolic planners to complete long-horizon tasks. On a Franka Emika robot arm, we find that EMBER enables the robot to complete three long-horizon visuomotor tasks at 85% success rate, such as organizing an office desk, a file cabinet, and drawers, which require sequencing up to 12 skills, involve 14 unique learned primitives, and demand generalization to novel objects.


title: Gambler ruin random walks and brownian motions in reserves modeling, application to pensions funds sustainability

abstract: We used the random walk to model the problem of reserves. The classic case of a stochastic process is the example of random walks, which are used to study a set of phenomena and, particularly, as in this article, models of reserves evolution. Random walks also allow the construction of significant complex systems and are also used as an instrument of analysis, being used in the sense of giving a theoretical characteristic to other types of systems. Our goal is primarily to study reserves to see how to ensure that pension funds are sustainable. This classic approach to the study of pension funds makes possible to draw interesting conclusions about the problem of reserves.


title: Chiral Surface Wave propagation with Anomalous Spin-momentum Locking

abstract: The ability to control the directionality of surface waves by manipulating its polarization has been of great significance for applications in spintronics and polarization-based optics. Surface waves with evanescent tails are found to possess an inherent in-plane transverse spin which is dependent on the propagation direction while an out-of-plane transverse spin does not naturally occur for surface waves and requires a specific surface design. Here, we introduce a new type of surface waves called chiral surface wave which has two transverse spins, an in-plane one that is inherent to any surface wave and an out-of-plane spin which is enforced by the design due to strong x-to-y coupling and broken rotational symmetry. The two transverse spins are locked to the momentum. Our study opens a new direction for metasurface designs with enhanced and controlled spin-orbit interaction by adding an extra degree of freedom to control the propagation direction as well as the transverse spin of surface waves.


title: Critical Tori for Mean Curvature Energies in Killing Submersions

abstract: We study surface energies depending on the mean curvature in total spaces of Killing submersions, which extend the classical notion of Willmore energy. Based on a symmetry reduction procedure, we construct vertical tori critical for these mean curvature energies. These vertical tori are based on closed curves critical for curvature energy functionals in Riemannian 2-space forms. The binormal evolution of these critical curves in Riemannian 3-space forms generates rotational tori solutions for an ample family of Weingarten surfaces. Therefore, we also introduce some correspondence results between these two types of tori and illustrate their relation.


title: Direct and clean loading of nanoparticles into optical traps at millibar pressures

abstract: Nanoparticles levitated by optical fields under vacuum conditions have applications in quantum science, the study of nanothermodynamics and precision sensing. Existing techniques for loading optical traps require ambient conditions, and often involve dispersion in liquids, which can contaminate delicate optics and lead to enhanced optical absorption and heating. Here we present a clean, dry and generic mechanism for directly loading optical traps at pressures down to 1\,mbar, exploiting Laser Induced Acoustic Desorption. Our method allows rapid and efficient trapping, and is also suitable for site-selective loading of nanofabricated particles grown on a silicon substrate.


title: Introduction to Neural Network Verification

abstract: Deep learning has transformed the way we think of software and what it can do. But deep neural networks are fragile and their behaviors are often surprising. In many settings, we need to provide formal guarantees on the safety, security, correctness, or robustness of neural networks. This book covers foundational ideas from formal verification and their adaptation to reasoning about neural networks and deep learning.


title: Nematic liquid crystals in a rectangular confinement: solution landscape and bifurcation

abstract: We study the solution landscape and bifurcation diagrams of nematic liquid crystals confined on a rectangle, using a reduced two-dimensional Landau--de Gennes framework in terms of two geometry-dependent variables: half short edge length $\lambda$ and aspect ratio $b$. First, we analytically prove that, for any $b$ with a small enough $\lambda$ or for a large enough $b$ with a fixed domain size, there is a unique stable solution that has two line defects on the opposite short edges. Second, we numerically construct solution landscapes by varying $\lambda$ and $b$, and report a novel X state, which emerges from saddle-node bifurcation and serves as the parent state in such a solution landscape. Various new classes are then found among these solution landscapes, including the X class, the S class, and the L class. By tracking the Morse indices of individual solutions, we present bifurcation diagrams for nematic equilibria, thus illustrating the emergence mechanism of critical points and several effects of geometrical anisotropy on confined defect patterns.


title: Community detection for weighted bipartite networks

abstract: The bipartite network appears in various areas, such as biology, sociology, physiology, and computer science. \cite{rohe2016co} proposed Stochastic co-Blockmodel (ScBM) as a tool for detecting community structure of binary bipartite graph data in network studies. However, ScBM completely ignores edge weight and is unable to explain the block structure of a weighted bipartite network. Here, to model a weighted bipartite network, we introduce a Bipartite Distribution-Free model by releasing ScBM's distribution restriction. We also build an extension of the proposed model by considering the variation of node degree. Our models do not require a specific distribution on generating elements of the adjacency matrix but only a block structure on the expected adjacency matrix. Spectral algorithms with theoretical guarantees on the consistent estimation of node labels are presented to identify communities. Our proposed methods are illustrated by simulated and empirical examples.


title: Operatorial formulation of the ghost rotationally-invariant slave-Boson theory

abstract: We propose an operatorially-exact formalism for describing the equilibrium and quantum-dynamical properties of many-electron systems interacting locally on a lattice, called "ghost rotationally-invariant slave-Boson theory" (g-RISB). We demonstrate that our theoretical framework reduces to the recently-developed ghost Gutzwiller approximation (g-GA) at the mean-field level. Furthermore, we introduce the time-dependent mean-field g-RISB action, generalizing the time-dependent GA theory. Since the g-RISB is based on exact reformulation of the many-body problem, it may pave the way to the development of practical implementations for adding systematically quantum-fluctuation corrections towards the exact solution, in arbitrary dimension.


title: Detectability of Population III stellar remnants as X-ray binaries from tidal captures in the local Universe

abstract: We assess the feasibility of detecting the compact object remnants from Population III (Pop III) stars in nearby dense star clusters, where they become luminous again as X-ray binaries (XRBs) and tidal disruption events (TDEs) via strong tidal encounters. Analytically modelling the formation of Pop III stars, coupled with a top-heavy initial mass function predicted by numerical simulations, we derive the number of (active) Pop III XRBs and TDEs in the present-day Milky Way (MW) nuclear star cluster as $\sim 0.06-0.3$ and $\lesssim 4\times 10^{-6}$, rendering any detection unlikely. The detection probability, however, can be significantly boosted when surveying all massive star clusters from the MW and neighboring galaxy clusters. Specifically, we predict $\sim 1.5-6.5$ and $\sim 40-2800$ active Pop III XRBs in the MW and the Virgo cluster, respectively. Our Pop III XRBs are dominated ($\sim 99\%$) by black holes with a typical mass and luminosity of $\sim 45$ $\rm M_{\odot}$ and $\sim 10^{36}$ $\rm erg\ s^{-1}$. Deep surveys of nearby ($\lesssim 30-300$ $\rm Mpc$) galaxy clusters for such Pop III XRBs are well within reach of next-generation X-ray telescopes, such as ATHENA and LYNX.


title: CondNet: Conditional Classifier for Scene Segmentation

abstract: The fully convolutional network (FCN) has achieved tremendous success in dense visual recognition tasks, such as scene segmentation. The last layer of FCN is typically a global classifier (1x1 convolution) to recognize each pixel to a semantic label. We empirically show that this global classifier, ignoring the intra-class distinction, may lead to sub-optimal results. In this work, we present a conditional classifier to replace the traditional global classifier, where the kernels of the classifier are generated dynamically conditioned on the input. The main advantages of the new classifier consist of: (i) it attends on the intra-class distinction, leading to stronger dense recognition capability; (ii) the conditional classifier is simple and flexible to be integrated into almost arbitrary FCN architectures to improve the prediction. Extensive experiments demonstrate that the proposed classifier performs favourably against the traditional classifier on the FCN architecture. The framework equipped with the conditional classifier (called CondNet) achieves new state-of-the-art performances on two datasets. The code and models are available at https://git.io/CondNet.


title: Simultaneous dilation and translation tilings of $\mathbb R^n$

abstract: We solve the wavelet set existence problem. That is, we characterize the full-rank lattices $\Gamma\subset \mathbb R^n$ and invertible $n \times n$ matrices $A$ for which there exists a measurable set $W$ such that $\{W + \gamma: \gamma \in \Gamma\}$ and $\{A^j(W): j\in \mathbb Z\}$ are tilings of $\mathbb R^n$. The characterization is a non-obvious generalization of the one found by Ionascu and Wang, which solved the problem in the case $n = 2$. As an application of our condition and a theorem of Margulis, we also strengthen a result of Dai, Larson, and the second author on the existence of wavelet sets by showing that wavelet sets exist for matrix dilations, all of whose eigenvalues $\lambda$ satisfy $|\lambda| \ge 1$. As another application, we show that the Ionascu-Wang characterization characterizes those dilations whose product of two smallest eigenvalues in absolute value is $\ge 1$.


title: Optical and spin-coherence properties of rubidium atoms trapped in solid neon

abstract: In this work, we measure the properties of ensembles of rubidium atoms trapped in solid neon that are relevant for use as quantum sensors of magnetic fields: the spin coherence of the trapped atoms and the ability to optically control and measure their spin state. We use the rubidium atoms as an AC magnetometer (by employing an appropriate dynamical decoupling sequence) and demonstrate NMR detection of Ne-21 atoms co-trapped in the neon matrix.


title: Hierarchy of Ideal Flatbands in Chiral Twisted Multilayer Graphene Models

abstract: We propose models of twisted multilayer graphene that exhibit exactly flat Bloch bands with arbitrary Chern numbers and ideal band geometries. The models are constructed by twisting two sheets of Bernal-stacked multiple graphene layers with only inter-sublattice couplings. Analytically we show that flatband wavefunctions in these models exhibit a momentum space holomorphic character, leading to ideal band geometries. We also explicitly demonstrate a generic "wavefunction exchange" mechanism that generates the high Chern numbers of these ideal flatbands. The ideal band geometries and high Chern numbers of the flatbands imply the possibility of hosting exotic fractional Chern insulators which do not have analogues in continuum Landau levels. We numerically verify that these exotic fractional Chern insulators are model states for short-range interactions, characterized by exact ground-state degeneracies at zero energy and infinite particle-cut entanglement gaps.


title: Structured transmittance illumination coherence holography

abstract: The coherence holography offers an unconventional way to reconstruct the hologram where an incoherent light illumination is used for reconstruction purposes, and object encoded into the hologram is reconstructed as the distribution of the complex coherence function. Measurement of the coherence function usually requires an interferometric setup and array detectors. This paper presents an entirely new idea of reconstruction of the complex coherence function in the coherence holography without an interferometric setup. This is realized by structured pattern projections on the incoherent source structure and implementing measurement of the cross-covariance of the intensities by a single-pixel detector. This technique, named structured transmittance illumination coherence holography (STICH), helps to reconstruct the complex coherence from the intensity measurement in a single-pixel detector without an interferometric setup and also keeps advantages of the intensity correlations. A simple experimental setup is presented as a first step to realize the technique, and results based on the computer modeling of the experimental setup are presented to show validation of the idea.


title: Data-driven controllers and the need for perception systems in underwater manipulation

abstract: The underwater environment poses a complex problem for developing autonomous capabilities for Underwater Vehicle Manipulator Systems (UVMSs). The modeling of UVMSs is a complicated and costly process due to the highly nonlinear dynamics and the presence of unknown hydrodynamical effects. This is aggravated in tasks where the manipulation of objects is necessary, as this may not only introduce external disturbances that can lead to a fast degradation of the control system performance, but also requires the coordinating with a vision system for the correct grasping and operation of the object. In this article, we introduce a control strategy for UVMSs working with unknown payloads. The proposed control strategy is based on a data-driven optimal controller. We present a number of experimental results showing the benefits of the proposed strategy. Furthermore, we include a discussion regarding the visual perception requirements for the UVMS in order to achieve full autonomy in underwater manipulation tasks of unknown payloads.


title: Gorenstein points in $\mathbb{P}^3$ via Hadamard product of projective varieties

abstract: We show how to construct a stick figure of lines in $\mathbb{P}^3$ using the Hadamard product of projective varieties. Then, applying the results of Migliore and Nagel, we use such stick figure to build a Gorenstein set of points with given $h-$vector ${\mathbf h}$. Since the Hadamard product is a coordinate-wise product, we show, at the end, how the coordinates of the points, in the Gorenstein set, can be directly determined.


title: Homography augumented momentum constrastive learning for SAR image retrieval

abstract: Deep learning-based image retrieval has been emphasized in computer vision. Representation embedding extracted by deep neural networks (DNNs) not only aims at containing semantic information of the image, but also can manage large-scale image retrieval tasks. In this work, we propose a deep learning-based image retrieval approach using homography transformation augmented contrastive learning to perform large-scale synthetic aperture radar (SAR) image search tasks. Moreover, we propose a training method for the DNNs induced by contrastive learning that does not require any labeling procedure. This may enable tractability of large-scale datasets with relative ease. Finally, we verify the performance of the proposed method by conducting experiments on the polarimetric SAR image datasets.


title: A Bayesian hierarchical model for disease mapping that accounts for scaling and heavy-tailed latent effects

abstract: In disease mapping, the relative risk of a disease is commonly estimated across different areas within a region of interest. The number of cases in an area is often assumed to follow a Poisson distribution whose mean is decomposed as the product between an offset and the logarithm of the disease's relative risk. The log risk may be written as the sum of fixed effects and latent random effects. The BYM2 model decomposes each latent effect into a weighted sum of independent and spatial effects. We build on the BYM2 model to allow for heavy-tailed latent effects and accommodate potentially outlying risks, after accounting for the fixed effects. We assume a scale mixture structure wherein the variance of the latent process changes across areas and allows for outlier identification. We propose two prior specifications for this scale mixture parameter. These are compared through simulation studies and in the analysis of Zika cases from the first (2015-2016) epidemic in Rio de Janeiro city, Brazil. The simulation studies show that, in terms of the model assessment criterion WAIC and outlier detection, the two proposed parametrisations perform better than the model proposed by Congdon (2017) to capture outliers. In particular, the proposed parametrisations are more efficient, in terms of outlier detection, than Congdon's when outliers are neighbours. Our analysis of Zika cases finds 19 out of 160 districts of Rio as potential outliers, after accounting for the socio-development index. Our proposed model may help prioritise interventions and identify potential issues in the recording of cases.


title: Characteristic polynomials of random truncations: moments, duality and asymptotics

abstract: We study moments of characteristic polynomials of truncated Haar distributed matrices from the three classical compact groups O(N), U(N) and Sp(2N). For finite matrix size we calculate the moments in terms of hypergeometric functions of matrix argument and give explicit integral representations highlighting the duality between the moment and the matrix size as well as the duality between the orthogonal and symplectic cases. Asymptotic expansions in strong and weak non-unitarity regimes are obtained. Using the connection to matrix hypergeometric functions, we establish limit theorems for the log-modulus of the characteristic polynomial evaluated on the unit circle.


title: On the purity conjecture of Nisnevich for torsors under reductive group schemes

abstract: Let $R$ be a regular semilocal integral domain containing an infinite field $k$. Let $f\in R$ be an element such that for all maximal ideals $\mathfrak m$ of $R$ we have $f\notin\mathfrak m^2$. Let $\mathbf G$ be a reductive group scheme over $R$. Under an isotropy assumption on $\mathbf G$ we show that a $\mathbf G$-torsor over the localization $R_f$ is trivial, provided it is rationally trivial. We show that it is not true without the isotropy assumption. Finally, if $\mathbf G$ is a commutative group scheme of multiplicative type and the regular semilocal ring contains a field of characteristic zero, we prove an analogue of Nisnevich purity conjecture for higher \'etale cohomology groups. The first statement is derived from its abstract version concerning presheaves of pointed sets satisfying some properties. The counterexample is constructed by providing a torsor over a local family of affine lines that cannot be extended to the family of projective lines. The latter is accomplished using the technique of affine Grassmannians.


title: Fine-grained Meta-Theorems for Vertex Integrity

abstract: Vertex Integrity is a graph measure which sits squarely between two more well-studied notions, namely vertex cover and tree-depth, and that has recently gained attention as a structural graph parameter. In this paper we investigate the algorithmic trade-offs involved with this parameter from the point of view of algorithmic meta-theorems for First-Order (FO) and Monadic Second Order (MSO) logic. Our positive results are the following: (i) given a graph $G$ of vertex integrity $k$ and an FO formula $\phi$ with $q$ quantifiers, deciding if $G$ satisfies $\phi$ can be done in time $2^{O(k^2q+q\log q)}+n^{O(1)}$; (ii) for MSO formulas with $q$ quantifiers, the same can be done in time $2^{2^{O(k^2+kq)}}+n^{O(1)}$. Both results are obtained using kernelization arguments, which pre-process the input to sizes $2^{O(k^2)}q$ and $2^{O(k^2+kq)}$ respectively. The complexities of our meta-theorems are significantly better than the corresponding meta-theorems for tree-depth, which involve towers of exponentials. However, they are worse than the roughly $2^{O(kq)}$ and $2^{2^{O(k+q)}}$ complexities known for corresponding meta-theorems for vertex cover. To explain this deterioration we present two formula constructions which lead to fine-grained complexity lower bounds and establish that the dependence of our meta-theorems on $k$ is the best possible. More precisely, we show that it is not possible to decide FO formulas with $q$ quantifiers in time $2^{o(k^2q)}$, and that there exists a constant-size MSO formula which cannot be decided in time $2^{2^{o(k^2)}}$, both under the ETH. Hence, the quadratic blow-up in the dependence on $k$ is unavoidable and vertex integrity has a complexity for FO and MSO logic which is truly intermediate between vertex cover and tree-depth.


title: Probing gas kinematics and PDR structure around O-type stars in Sh 2-305 HII region

abstract: We report an observational study of the Galactic HII region Sh 2-305/S305 using the [CII] 158 $\mu$m line data, which are used to examine the gas dynamics and structure of photodissociation regions. The integrated [CII] emission map at [39.4, 49.5] km s$^{-1}$ spatially traces two shell-like structures (i.e., inner and outer neutral shells) having a total mass of $\sim$565 M$_\odot$. The inner neutral shell encompasses an O9.5V star at its centre and has a compact ring-like appearance. However, the outer shell is seen with more extended and diffuse [CII] emission, hosting an O8.5V star at its centre and surrounds the inner neutral shell. The velocity channel maps and position-velocity diagrams confirm the presence of a compact [CII] shell embedded in the diffuse outer shell, and both the shells seem to expand with $v_{\rm exp}\sim$1.3 km s$^{-1}$. The outer shell appears to be older than the inner shell, hinting that these shells are formed sequentially. The [CII] profiles are examined toward S305, which are either double-peaked or blue-skewed and have the brighter redshifted component. The redshifted and blueshifted components spatially trace the inner and outer neutral shell geometry, respectively. The ionized, neutral, and molecular zones in S305 are seen adjacent to one another around the O-type stars. The regularly spaced dense molecular and dust clumps (mass $\sim$10-10$^{3}$ M$_{\odot}$) are investigated around the neutral shells, which might have originated due to gravitational instability in the shell of collected materials.


title: Effective field theory of fluctuating wall in open systems: from a kink in Josephson junction to general domain wall

abstract: We investigate macroscopic behaviors of fluctuating domain walls in nonequilibrium open systems with the help of the effective field theory based on symmetry. Since the domain wall in open systems breaks the translational symmetry, there appears a gapless excitation identified as the Nambu-Goldstone (NG) mode, which shows the non-propagating diffusive behavior in contrast to those in closed systems. After demonstrating the presence of the diffusive NG mode in the $(2+1)$-dimensional dissipative Josephson junction, we provide a symmetry-based general analysis for open systems breaking the one-dimensional translational symmetry. A general effective Lagrangian is constructed based on the Schwinger-Keldysh formalism, which supports the presence of the gapless diffusion mode in the fluctuation spectrum in the thin wall regime. Besides, we also identify a term peculiar to the open system, which possibly leads to the instability in the thick-wall regime or the nonlinear Kardar-Parisi-Zhang coupling in the thin-wall regime although it is absent in the Josephson junction.


title: Is backreaction in cosmology a relativistic effect? On the need for an extension of Newton's theory to non-Euclidean topologies

abstract: Cosmological backreaction corresponds to the effect of inhomogeneities of structure on the global expansion of the Universe. The main question surrounding this phenomenon is whether or not it is important enough to lead to measurable effects on the scale factor evolution eventually explaining its acceleration or the Hubble tension. One of the most important result on this subject is the Buchert-Ehlers theorem (Buchert \& Ehlers, 1997) stating that backreaction is exactly zero when calculated using Newton's theory of gravitation, which may not be the case in general relativity. It is generally said that this result implies that backreaction is a purely relativistic effect. We will show that this is not necessarily the case, in the sense that this implication does not apply to a universe which is still well described by Newton's theory on small scales but has a non-Euclidean topology. The theorem should therefore be generalised to account for such a scenario. In a heuristic calculation where we construct a theory which is locally Newtonian but defined on a non-Euclidean topology, we show that backreaction is non-zero, meaning that it might be non-relativistic depending on the topological class of our Universe. However, that construction is not unique and remains to be justified from a non-relativistic limit of general relativity.


title: Funny business from the large $N_c$ finite temperature crossover

abstract: It is well known that the deconfinement transition temperature for $SU(N_c)$ gauge theory is almost independent of $N_c$, and the transition is first order for $N_c \ge 3$. In the real world ($N_c=3$, light quarks) it is a crossover located far away from the pure gauge value. What happens to the transition temperature at fixed fermion mass if the number of fermion flavors is held constant ($N_f=2$) and $N_c$ is varied? There are multiple plausible stories, only one of which appears to be true when the systems are simulated on the lattice. I describe the physics issues which surround the question and my lattice - based answer to it.


title: Wave Function Identity: A New Symmetry for 2-electron Systems in an Electromagnetic Field

abstract: Stationary-state Schr{\"o}dinger-Pauli theory is a description of electrons with a spin moment in an external electromagnetic field. For 2-electron systems as described by the Schr{\"o}dinger-Pauli theory Hamiltonian with a symmetrical binding potential, we report a new symmetry operation of the electronic coordinates. The symmetry operation is such that it leads to the equality of the transformed wave function to the wave function. This equality is referred to as the Wave Function Identity. The symmetry operation is a two-step process: an interchange of the spatial coordinates of the electrons whilst keeping their spin moments unchanged, followed by an inversion. The Identity is valid for arbitrary structure of the binding potential, arbitrary electron interaction of the form $w(|{\bf{r}} - {\bf{r}}'|)$, all bound electronic states, and arbitrary dimensionality. It is proved that the exact wave functions satisfy the Identity. On application of the permutation operation for fermions to the identity, it is shown that the parity of the singlet states is even and that of triplet states odd. As a consequence, it follows that at electron-electron coalescence, the singlet state wave functions satisfy the cusp coalescence constraint, and triplet state wave functions the node coalescence condition. Further, we show that the parity of the singlet state wave functions about all points of electron-electron coalescence is even, and that of the triplet state wave functions odd. The Wave Function Identity and the properties on parity, together with the Pauli principle, are then elucidated by application to the $2$-dimensional $2$-electron `artificial atoms' or semiconductor quantum dots in a magnetic field in their first excited singlet $2^{1}S$ and triplet $2^{3} S$ states.


title: Lack of Value Definiteness in Quantum Gravity

abstract: In this paper we consider the observables describing fundamental spatiotemporal properties and relations in the context of Quantum Gravity (QG). As we will show, in both Loop Quantum Gravity and in String Theory, these observables are non-commuting ones. By analogy with some recent arguments put forward in the context of non-relativistic quantum mechanics (QM), we shall propose to interpret these physical quantities as ontologically indeterminate -- i.e., indeterminate in a non-epistemic, mind-independent way. This result has two important consequences for current debates in philosophy of physics. First, it shows that ontological indeterminacy may extend well beyond the non-relativistic case, thereby also suggesting a conceptual continuity between QM and QG. Second, when applied to QG, the notion of indeterminacy may provide a unified philosophical framework for understanding very distinct approaches that are usually considered incompatible.


title: Torque-free manipulation of nanoparticle rotations via embedded spins

abstract: Spin angular momentum and mechanical rotation both contribute to the total angular momentum of rigid bodies, leading to spin-rotational coupling via the Einstein-de Haas and Barnett effects. Here we show that the revolutions of symmetric nanorotors can be strongly affected by a small number of intrinsic spins. The resulting dynamics are observable with freely rotating nanodiamonds with embedded nitrogen-vacancy centers and persist for realistically-shaped near-symmetric particles, opening the door to torque-free schemes to control their rotations at the quantum level.


title: Multilingual Document-Level Translation Enables Zero-Shot Transfer From Sentences to Documents

abstract: Document-level neural machine translation (DocNMT) achieves coherent translations by incorporating cross-sentence context. However, for most language pairs there's a shortage of parallel documents, although parallel sentences are readily available. In this paper, we study whether and how contextual modeling in DocNMT is transferable via multilingual modeling. We focus on the scenario of zero-shot transfer from teacher languages with document level data to student languages with no documents but sentence level data, and for the first time treat document-level translation as a transfer learning problem. Using simple concatenation-based DocNMT, we explore the effect of 3 factors on the transfer: the number of teacher languages with document level data, the balance between document and sentence level data at training, and the data condition of parallel documents (genuine vs. backtranslated). Our experiments on Europarl-7 and IWSLT-10 show the feasibility of multilingual transfer for DocNMT, particularly on document-specific metrics. We observe that more teacher languages and adequate data balance both contribute to better transfer quality. Surprisingly, the transfer is less sensitive to the data condition, where multilingual DocNMT delivers decent performance with either backtranslated or genuine document pairs.


title: Infraparticle quantum fields and the formation of photon clouds

abstract: A non-perturbative and exactly solvable quantum field theoretical model for a "dressed Dirac field" is presented, that exhibits all the kinematical features of QED: an appropriate delocalization of the charged field as a prerequisite for the global Gauss Law, superselected photon clouds (asymptotic expectation values of the Maxwell field), infraparticle nature of charged particles that cannot be separated from their photon clouds, broken Lorentz symmetry. The model serves as an intermediate leg on a new roadmap towards full QED, formulated as an off-shell theory, i.e., including a perturbative construction of its interacting charged fields. It also fills a gap in recent discussions of the "Infrared Triangle", and points the way towards a new scattering theory for theories with massless particles of helicity $\geq 1$, in which infraparticles can be described at the level of charged fields, rather than just states.


title: Generalized Larcombe-Fenessey invariants of matrix powers

abstract: In this article, we have found significant generalization of the invariance properties of powers of matrices discovered by Larcombe, Fenessey and further explored by Zeilberger. Moreover, we found interesting new results exhibiting similar phenomena in a more general setup.


title: Distribution of certain $\ell$-regular partitions and triangular numbers

abstract: Let $pod_{\ell}(n)$ be the number of $\ell$-regular partitions of $n$ with distinct odd parts. In this article, prove that for any positive integer $k$, the set of non-negative integers $n$ for which $pod_{\ell}(n)\equiv 0 \pmod{p^{k}}$ has density one under certain conditions on $\ell$ and $p$. For $p \in \{3,5,7\}$, we also exhibit multiplicative identities for $pod_{p}(n)$ modulo $p.$


title: Classical Double Copy of Worldline Quantum Field Theory

abstract: The recently developed worldline quantum field theory (WQFT) formalism for the classical gravitational scattering of massive bodies is extended to massive, charged point particles coupling to bi-adjoint scalar field theory, Yang-Mills theory, and dilaton-gravity. We establish a classical double copy relation in these WQFTs for classical observables (deflection, radiation). The bi-adjoint scalar field theory fixes the locality structure of the double copy from Yang-Mills to dilaton-gravity. Using this the eikonal scattering phase (or free energy of the WQFT) is computed to next-to-leading order (NLO) in coupling constants using the double copy as well as directly finding full agreement. We clarify the relation of our approach to previous studies in the effective field theory formalism. Finally, the equivalence of the WQFT double copy to the double copy relation of the classical limit of quantum scattering amplitudes is shown explicitly up to NLO.


title: Relation-Guided Pre-Training for Open-Domain Question Answering

abstract: Answering complex open-domain questions requires understanding the latent relations between involving entities. However, we found that the existing QA datasets are extremely imbalanced in some types of relations, which hurts the generalization performance over questions with long-tail relations. To remedy this problem, in this paper, we propose a Relation-Guided Pre-Training (RGPT-QA) framework. We first generate a relational QA dataset covering a wide range of relations from both the Wikidata triplets and Wikipedia hyperlinks. We then pre-train a QA model to infer the latent relations from the question, and then conduct extractive QA to get the target answer entity. We demonstrate that by pretraining with propoed RGPT-QA techique, the popular open-domain QA model, Dense Passage Retriever (DPR), achieves 2.2%, 2.4%, and 6.3% absolute improvement in Exact Match accuracy on Natural Questions, TriviaQA, and WebQuestions. Particularly, we show that RGPT-QA improves significantly on questions with long-tail relations


title: Security-Hardening Software Libraries with Ada and SPARK -- A TCP Stack Use Case

abstract: This white paper demonstrates how the assurance, reliability, and security of an existing professional-grade, open-source embedded TCP/IP stack implementation written in the C programming language is significantly enhanced by adopting the SPARK technology. A multifaceted approach achieves this. Firstly, the TCP layer's C code is being replaced with formally verified SPARK, a subset of the Ada programming language supported by formal verification tools. Then the lower layers, still written in C and on which the TCP layer depends, are modeled using SPARK contracts and validated using symbolic execution with KLEE. Finally, formal contracts for the upper layers are defined to call the TCP layer. The work allowed the detection and correction of two bugs in the TCP layer. In an increasingly connected world, where Cyber Security is of paramount importance, the powerful approach detailed in this work can be applied to any existing critical C library to harden their reliability and security significantly.


title: All that Glitters is not Gold: Relational Events Models with Spurious Events

abstract: As relational event models are an increasingly popular model for studying relational structures, the reliability of large-scale event data collection becomes more and more important. Automated or human-coded events often suffer from non-negligible false-discovery rates in event identification. And most sensor data is primarily based on actors' spatial proximity for predefined time windows; hence, the observed events could relate either to a social relationship or random co-location. Both examples imply spurious events that may bias estimates and inference. We propose the Relational Event Model for Spurious Events (REMSE), an extension to existing approaches for interaction data. The model provides a flexible solution for modeling data while controlling for spurious events. Estimation of our model is carried out in an empirical Bayesian approach via data augmentation. Based on a simulation study, we investigate the properties of the estimation procedure. To demonstrate its usefulness in two distinct applications, we employ this model to combat events from the Syrian civil war and student co-location data. Results from the simulation and the applications identify the REMSE as a suitable approach to modeling relational event data in the presence of spurious events.


title: Enabling variable high spatial resolution retrieval from a long pulse BOTDA sensor

abstract: In the field of Internet of Things, there is an urgent need for sensors with large-scale sensing capability for scenarios such as intelligent monitoring of production lines and urban infrastructure. Brillouin optical time domain analysis (BOTDA) sensors, which can monitor thousands of continuous points simultaneously, show great advantages in these applications. We propose a convolutional neural network (CNN) to process the data of conventional Brillouin optical time domain analysis (BOTDA) sensors, which achieves unprecedented performance improvement that allows to directly retrieve higher spatial resolution (SR) from the sensing system that use long pump pulses. By using the simulated Brillouin gain spectrums (BGSs) as the CNN input and the corresponding high SR BFS as the output target, the trained CNN is able to obtain a SR higher than the theoretical value determined by the pump pulse width. In the experiment, the CNN accurately retrieves 0.5-m hotspots from the measured BGS with pump pulses from 20 to 50 ns, and the acquired BFS is in great agreement with 45/40 ns differential pulse-width pair (DPP) measurement results. Compared with the DPP technique, the proposed CNN demonstrates a 2-fold improvement in BFS uncertainty with only half the measurement time. In addition, by changing the training datasets, the proposed CNN can obtain tunable high SR retrieval based on conventional BOTDA sensors that use long pulses without any requirement of hardware modifications. The proposed data post-processing approach paves the way to enable novel high spatial resolution BOTDA sensors, which brings substantial improvement over the state-of-the-art techniques in terms of system complexity, measurement time and reliability, etc.


title: Competitive X-ray and Optical Cooling in the Collisionless Shocks of WR 140

abstract: WR 140 is a long-period, highly eccentric Wolf-Rayet star binary system with exceptionally well-determined orbital and stellar parameters. Bright, variable X-ray emission is generated in shocks produced by the collision of the winds of the WC7pd+O5.5fc component stars. We discuss the variations in the context of the colliding-wind model using broad-band spectrometry from the RXTE, SWIFT, and NICER observatories obtained over 20 years and nearly 1000 observations through 3 consecutive 7.94-year orbits including 3 periastron passages. The X-ray luminosity varies as expected with the inverse of the stellar separation over most of the orbit: departures near periastron are produced when cooling shifts to excess optical emission in CIII $\lambda5696$ in particular. We use X-ray absorption to estimate mass-loss rates for both stars and to constrain the system morphology. The absorption maximum coincides closely with inferior conjunction of the WC star and provides evidence of the ion-reflection mechanism that underlies the formation of collisionless shocks governed by magnetic fields probably generated by the Weibel instability. Comparisons with K-band emission and HeI $\lambda$10830 absorption show that both are correlated after periastron with the asymmetric X-ray absorption. Dust appears within a few days of periastron suggesting formation within shocked gas near the stagnation point. X-ray flares seen in $\eta$ Carinae have not occurred in WR 140, suggesting the absence of large-scale wind inhomogeneities. Relatively constant soft emission revealed during the X-ray minimum is probably not from recombining plasma entrained in outflowing shocked gas.


title: STAR HE-Linac Complete Detailed Design Report

abstract: This Document contains a complete technical description of the system devoted to the upgrade of the STAR Linear Accelerator (STAR Linac). According to the Contract signed between Universita della Calabria (UniCal) and Istituto Nazionale di Fisica Nucleare (INFN) on May 7th, 2021, INFN is committed to install, test and commission the upgrade of the STAR Linac denominated STAR-HE-Linac (STAR High Energy Linac), hereafter STAR-HEL. The technical components as well as the installation / test procedures and the ancillary equipment involved in such an upgrade are the object of this Document, named Complete Detailed Design Report. A technical offer was submitted by INFN in the frame of its participation to the tender issued by UniCal, describing a possible energy upgrade of the STAR Linac, with an electron beam energy boosted from 65 MeV up to 150 MeV by means of Radiofrequency (RF) accelerating sections and power stations based on S-band technology (ie, 2856 MHz RF frequency). Following the Contract signature, INFN conceived and conceptually designed a technology change that offers several advantages both on performances and on operational reliability of STAR-HEL, based on adopting C-band technology (ie, 5712 MHz RF frequency) for accelerating sections and RF power stations. Such a technology change was illustrated in a dedicated document named "improvement option for STAR-2 HE-Linac", addressed to UniCal STAR Management Board and tender R.U.P. for approval on June 9th, 2021. INFN received a formal letter of approval on June 16th, 2021.


title: Rapid stellar and binary population synthesis with COMPAS

abstract: Compact Object Mergers: Population Astrophysics and Statistics (COMPAS; https://compas.science) is a public rapid binary population synthesis code. COMPAS generates populations of isolated stellar binaries under a set of parametrized assumptions in order to allow comparisons against observational data sets, such as those coming from gravitational-wave observations of merging compact remnants. It includes a number of tools for population processing in addition to the core binary evolution components. COMPAS is publicly available via the github repository https://github.com/TeamCOMPAS/COMPAS/, and is designed to allow for flexible modifications as evolutionary models improve. This paper describes the methodology and implementation of COMPAS. It is a living document which will be updated as new features are added to COMPAS; the current document describes COMPAS v02.21.00.


title: An Ultra-Fast Method for Simulation of Realistic Ultrasound Images

abstract: Convolutional neural networks (CNNs) have attracted a rapidly growing interest in a variety of different processing tasks in the medical ultrasound community. However, the performance of CNNs is highly reliant on both the amount and fidelity of the training data. Therefore, scarce data is almost always a concern, particularly in the medical field, where clinical data is not easily accessible. The utilization of synthetic data is a popular approach to address this challenge. However, but simulating a large number of images using packages such as Field II is time-consuming, and the distribution of simulated images is far from that of the real images. Herein, we introduce a novel ultra-fast ultrasound image simulation method based on the Fourier transform and evaluate its performance in a lesion segmentation task. We demonstrate that data augmentation using the images generated by the proposed method substantially outperforms Field II in terms of Dice similarity coefficient, while the simulation is almost 36000 times faster (both on CPU).


title: A Bernstein-type Inequality for High Dimensional Linear Processes with Applications to Robust Estimation of Time Series Regressions

abstract: Time series regression models are commonly used in time series analysis. However, in modern real-world applications, serially correlated data with an ultra-high dimension and fat tails are prevalent. This presents a challenge in developing new statistical tools for time series analysis. In this paper, we propose a novel Bernstein-type inequality for high-dimensional linear processes and apply it to investigate two high-dimensional robust estimation problems: (1) time series regression with fat-tailed and correlated covariates and errors, and (2) fat-tailed vector autoregression. Our proposed approach allows for exponential increases in dimension with sample size under mild moment and dependence conditions, while ensuring consistency in the estimation process.


title: Holographic Evolution with Dynamical Boundary Gravity

abstract: Holography has provided valuable insights into the time evolution of strongly coupled gauge theories in a fixed spacetime. However, this framework is insufficient if this spacetime is dynamical. We present a scheme to evolve a four-dimensional, strongly interacting gauge theory coupled to four-dimensional dynamical gravity in the semiclassical regime. As in previous work, we use holography to evolve the quantum gauge theory stress tensor, whereas the four-dimensional metric evolves according to Einstein's equations coupled to the expectation value of the stress tensor. The novelty of our approach is that both the boundary and the bulk spacetimes are constructed dynamically, one time step at a time. We focus on Friedmann-Lema\^itre-Robertson-Walker geometries and evolve far-from-equilibrium initial states that lead to asymptotically expanding, flat or collapsing Universes


title: Formation and evolution of young massive clusters in galaxy mergers: the SMUGGLE view

abstract: Galaxy mergers are known to host abundant young massive cluster (YMC) populations, whose formation mechanism is still not well-understood. Here, we present a high-resolution galaxy merger simulation with explicit star formation and stellar feedback prescriptions to investigate how mergers affect the properties of the interstellar medium and YMCs. Compared with a controlled simulation of an isolated galaxy, the mass fraction of dense and high-pressure gas is much higher in mergers. Consequently, the mass function of both molecular clouds and YMCs becomes shallower and extends to higher masses. Moreover, cluster formation efficiency is significantly enhanced and correlates positively with the star formation rate surface density and gas pressure. We track the orbits of YMCs and investigate the time evolution of tidal fields during the course of the merger. At an early stage of the merger, the tidal field strength correlates positively with YMC mass, $\lambda_{\rm tid}\propto M^{0.71}$, which systematically affects the shape of the mass function and age distribution of the YMCs. At later times, most YMCs closely follow the orbits of their host galaxies, gradually sinking into the center of the merger remnant due to dynamical friction, and are quickly dissolved via efficient tidal disruption. Interestingly, YMCs formed during the first passage, mostly in tidal tails and bridges, are distributed over a wide range of galactocentric radii, greatly increasing their survivability because of the much weaker tidal field in the outskirts of the merger system. These YMCs are promising candidates for globular clusters that survive to the present day.


title: Tasting Flavoured Majorana Dark Matter

abstract: We study a simplified model of flavoured Majorana dark matter in the Dark Minimal Flavour Violation framework. The model extends the Standard Model by a dark matter flavour triplet and a scalar mediator, through which the new dark fermions couple to right-handed up-type quarks. This interaction is governed by a new coupling matrix $\lambda$ which is assumed to constitute the only new source of flavour and CP violation. We analyse the parameter space of this model by using constraints from collider searches, $D^0-\bar{D}^0$ mixing, cosmology and direct dark matter searches. Throughout our study, we point out crucial differences between the Majorana and Dirac dark matter cases. After performing a combined analysis within the context of all the experimental constraints mentioned above, we analyse which flavour for the dark matter particle is preferred by experimental data. We further investigate if this model is capable of explaining the large measured value of the direct CP asymmetry $\Delta A^\text{dir}_{CP}$ in charm decays. We find that significant enhancements with respect to the Standard Model expectation are compatible with all constraints, and even the central value of the measurement can be reached. We also advertise the flavour-violating final state with two same-sign top quarks produced in association with missing transverse energy as a smoking-gun signature for flavoured Majorana dark matter at the LHC.


title: Characterizing Heavy Neutral Fermions via their Decays

abstract: Many extensions of the Standard Model of particle physics contain new electrically-neutral fermions. Should one of these particles be discovered, questions will naturally arise regarding its nature. For instance: is it a self-conjugate particle (i.e., is it a Dirac or a Majorana fermion)?, does it interact via the Standard Model force carriers or something else? One set of well-motivated particles in this class are Heavy Neutral Leptons (HNLs), Standard Model gauge-singlet fermions that mix with the neutrinos and may be produced in meson decays. We demonstrate that measuring the three body decays of the HNL (or phenomenologically similar heavy fermions) can help determine whether they are Majorana or Dirac fermions. We also investigate the ability to distinguish among different models for the physics responsible for the HNL decay. We compare the reach assuming full and partial event reconstruction, and propose experimental analyses. Should a new fermion be discovered, studying its three body decays provides a powerful diagnostic tool of its nature.


title: Hunting for tetraquarks in ultra-pheripheral heavy ion collisions

abstract: Ultra-peripheral heavy ion collisions constitute an ideal setup to look for exotic hadrons because of their low event multiplicity and the possibility of an efficient background rejection. We propose to look for four-quark states produced by photon-photon fusion in these collisions at center-of-mass energy per nucleon pair $\sqrt{s_{\rm NN}}=5.5~{\rm TeV}$. In particular, we focus on those states that would represent a definite smoking gun for the compact tetraquark model. We show that the $X(6900)$, a likely $cc\bar{c}\bar{c}$ compact state, is a perfect candidate for this search, and estimate a production cross section ranging from around $250$ nb to $1150$ nb, depending on its quantum numbers. Furthermore, we discuss the importance of ultra-peripheral collisions to the search for the scalar and tensor partners of the $X(3872)$ predicted by the diquarkonium model, and not yet observed. The completion of such a flavor-spin multiplet would strongly speak in favor of the compact tetraquark model.


title: Robust marginalization of baryonic effects for cosmological inference at the field level

abstract: We train neural networks to perform likelihood-free inference from $(25\,h^{-1}{\rm Mpc})^2$ 2D maps containing the total mass surface density from thousands of hydrodynamic simulations of the CAMELS project. We show that the networks can extract information beyond one-point functions and power spectra from all resolved scales ($\gtrsim 100\,h^{-1}{\rm kpc}$) while performing a robust marginalization over baryonic physics at the field level: the model can infer the value of $\Omega_{\rm m} (\pm 4\%)$ and $\sigma_8 (\pm 2.5\%)$ from simulations completely different to the ones used to train it.

categories: astro-ph.CO astro-ph.GA astro-ph.IM cs.CV cs.LG


title: Eccentric millisecond pulsars by resonant convection

abstract: Eccentric millisecond pulsars (eMSPs) with white dwarf companions exhibit orbital eccentricities orders of magnitude larger than predicted by turbulent convection in the white dwarfs' red giant progenitors. The orbital periods of eMSPs cluster around $P=20-30$ d, remarkably close to the red giant convective eddy turnover time $t_{\rm eddy}$. We propose that the anomalously large eccentricities are resonantly driven by convective flows somehow made coherent when the turnover time matches the tidally locked red giant's spin period, which is also the tidal forcing period. Numerical simulations of rotating red giants and magnetic field studies of stars show some evidence for especially ordered flow patterns when the convective Rossby number $P/t_{\rm eddy}$ is of order unity. We show that resonant convection boosts eccentricities by a factor of $(t_{\rm nuc}/P)^{1/2}\approx 3\times 10^3$ over the random-walk values that characterize conventional MSPs, in good agreement with observations ($t_{\rm nuc}$ is the giant's nuclear burning time-scale). We also show how variations in the eddy turnover time arising from red giant metallicity variations can reproduce the observed effective width of the resonance, $\Delta P/P\approx 0.4$.


title: Combined analysis of neutron star natal kicks using proper motions and parallax measurements for radio pulsars and Be X-ray binaries

abstract: Supernova explosion and the associated neutron star natal kicks are important events on a pathway of a binary to become a gravitational wave source, an X-ray binary or a millisecond radio pulsar. Weak natal kicks often lead to binary survival, while strong kicks frequently disrupt the binary. In this article, we aim to further constrain neutron star natal kicks in binaries. We explore binary population synthesis models by varying prescription for natal kick, remnant mass and mass accretion efficiency. We introduce a robust statistical technique to analyse combined observations of different nature. Using this technique, we further test different models using parallax and proper motion measurements for young isolated radio pulsars and similar measurements for Galactic Be X-ray binaries. Our best model for natal kicks is consistent with both measurements and contains a fraction of $w=0.2\pm 0.1$ weak natal kicks with $\sigma_1 = 45^{+25}_{-15}$ km/s, the remaining natal kicks are drawn from the high-velocity component, same as in previous works: $\sigma_2 = 336$ km/s. We found that currently used models for natal kicks of neutron stars produced by electron capture supernova (combination of maxwellian $\sigma=265$ km/s and $\sigma = 30$ km/s for electron capture) are inconsistent or marginally consistent with parallaxes and proper motions measured for isolated radio pulsars. We suggest a new model for natal kicks of ecSN, which satisfy both observations of isolated radio pulsars and Be X-ray binaries.


title: Self-consistent population spectral synthesis with FADO: mean stellar metallicity of galaxies in spectral synthesis methods

abstract: In this work, we investigate the reliability of spectral synthesis methods in the estimation of the mean stellar age and metallicity, addressing the question of which signal-to-noise ratios (S/N) are needed to determine these quantities. To address this problem we used simulated spectra containing stellar and nebular emission, reproducing the evolution of a galaxy for a constant and exponentially declining star formation law. The spectra have been degraded to different S/N and analysed with three different spectral synthesis codes: FADO, STARLIGHT, and STECKMAP assuming similar fitting set-ups and the same spectral bases. For S/N > 5 all tools considered show a large diversity in the results. FADO and STARLIGHT find median differences in light-weighted mean stellar ages of ~0.1 dex, while STECKMAP shows a higher value of ~0.2 dex. For S/N > 50 the median differences in FADO are ~0.03 dex (~7%), a factor 3 and 4 lower than the 0.08 dex (~20%) and 0.11 dex (~30%) obtained from STARLIGHT and STECKMAP, respectively. Our results indicate that phases of high specific star formation rate (sSFR) in galaxies require analysis tools that do not neglect the nebular continuum emission in the fitting process, since purely stellar models would have strong problems in the estimation of star formation history, even in presence of high S/N spectra. The median values of these differences are of the order of 7% (FADO), 20% (STARLIGHT), and 30% (STECKMAP) for light-weighted quantities, and 20% (FADO), 60% (STARLIGHT), and 20% (STECKMAP) for mass-weighted quantities. That implies a severe overestimation of the mass-to-light ratio and stellar mass, even in the presence of a mild contribution from the nebular continuum. Our work underlines the importance of a self-consistent treatment of nebular emission, which is the only route towards a reliable determination of the assembly of any high-sSFR galaxy.


title: Axions from Neutron Star Mergers

abstract: Axion-like particles (ALP) can in principle be produced in very hot and dense astrophysical environments, escape from the extreme object where such conditions are met, and then be converted in gamma--rays in the magnetic fields intervening between the event and the Earth. This process potentially offers a new window on both the physics of the axions, and the inner working of the astrophysical objects where they are produced. Interestingly, while this process has been studied for core--collapse supernovae and other extreme astrophysical events, no estimate exists for Neutron Star Mergers, objects recently identified through the detection of gravitational waves. In this work we study the production of ALPs in neutron star mergers, finding that for a large region of the ALP parameter space its magnitude at the source is such to produce a sizable gamma-ray signal at Earth. We show detection forecasts for such events placed in nearby galaxies, finding that they are potentially observable with the Fermi--LAT, thus opening a new window into both the astrophysics of these cataclysmic events, and of new particles beyond the standard model.


title: Study of $\mathrm{^{26}Al}$ in the COSI 2016 superpressure balloon flight

abstract: The Compton Spectrometer and Imager (COSI) is a balloon-borne compact Compton telescope designed to survey the $\gamma$-ray sky from 0.2 to 5 MeV. COSI's wide field-of-view (FOV) and excellent energy resolution from high-purity germanium detectors make it uniquely capable of probing this under-explored energy regime. In particular, it can facilitate understanding of stellar nucleosynthesis through studies of diffuse emission from the radioisotope $\mathrm{^{26}Al}$ at 1.809 MeV. In 2016, COSI was launched from Wanaka, New Zealand on a NASA superpressure balloon and flew for 46 days. The flight was a technologic and scientific success, boasting live detection and polarization studies of GRB160530A, spectral analysis of the Crab Nebula and the 511-keV positron annihilation emission at the Galactic Center, and detection of Cygnus X-1. This article details the first maximum-likelihood search for the 1.809 MeV signature of Galactic $\mathrm{^{26}Al}$ in the 2016 data. The analysis reveals a promising excess around the expected energies of an $\mathrm{^{26}Al}$ signature with 3.7$\sigma$ significance and a measured flux of (17.0 $\pm$ 4.9) $\times$ 10$^{-4}$ ph cm$^{-2}$ s$^{-1}$. Further exploration is currently underway to solidify the measurement.


title: Testing Exotic Scalars with HiggsBounds

abstract: The program HiggsBounds is a well-established tool for testing Beyond-the-Standard Model (BSM) theories with an extended Higgs sector against experimental limits from collider searches at LEP, Tevatron and LHC. Thus far, it could be applied to any neutral or charged Higgs bosons originating from the modified Higgs sector. Implicitly, these particles were assumed to exhibit a somewhat hierarchical Yukawa structure as present in the Standard Model, where in particular the couplings to first generation fermions could be neglected. In this work, we extend the HiggsBounds functionalities to go beyond these restrictions, thus making the code applicable to any neutral or charged BSM scalars. Moreover, we develop a new approach to implement experimental searches whose kinematic acceptance depends significantly on the values of the involved couplings. We achieve this by recasting the searches to general scalar models. Using this approach we incorporate relevant current experimental limits from LHC searches for exotic scalars, and present the implications of these limits for a dark matter scalar mediator model, a flipped Two-Higgs-Doublet Model and a supersymmetric model with R-parity violation.


title: Pole inflation from non-minimal coupling to gravity

abstract: Transforming canonical scalars to the Einstein frame can give a multi-field generalization of pole inflation (namely, a scalar with a divergent kinetic term) at vanishing field-dependent Planck mass. However, to obtain an attractor, the scalar potential must obey certain non-generic conditions. These are automatically satisfied in Quantum Field Theories with dimension-less couplings. The resulting models of pole inflation have special inflationary predictions determined by the full RG running of couplings. Acceptable predictions for the tensor/scalar ratio arise for perturbative but moderately large couplings, so we explore the possible QFT runnings: to confinement, to an IR fixed point, and to a UV fixed point.


title: Stellar Evolution confronts Axion Models

abstract: Axion production from astrophysical bodies is a topic in continuous development, because of theoretical progress in the estimate of stellar emission rates and, especially, because of improved stellar observations. We carry out a comprehensive analysis of the most informative astrophysics data, revisiting the bounds on axion couplings to photons, nucleons and electrons, and reassessing the significance of various hints of anomalous stellar energy losses. We confront the performance of various theoretical constructions in accounting for these hints, while complying with the observational limits on axion couplings. We identify the most favorable models, and the regions in the mass/couplings parameter space which are preferred by the global fit. Finally, we scrutinize the discovery potential for such models at upcoming helioscopes, namely IAXO and its scaled versions.


title: Void Galaxy Distribution: A Challenge for $\Lambda$CDM

abstract: We extract void catalogs from the Sloan Digital Sky Survey Data Release 16 (SDSS DR16) survey and also from the Millennium simulation. We focus our comparison on distribution of galaxies brighter than $M_r < -18$ inside voids and study the mean separation of void galaxies, distance from the void center, and the radial density profile. We find that mean separation of void galaxies depends on void size, as bigger voids have lower mean separation in both samples. However, void galaxies in the observation sample seem to have generally larger mean-distance than simulated ones at any given void size. In addition, observed void galaxies tend to reside closer to the void center than those in the simulation. This discrepancy is also shown in the density profile of voids. Regardless of the void size, the central densities of real void profiles are higher than the ones in the predicted simulated catalog.


title: The role of the chiral anomaly in polarized deeply inelastic scattering II: Topological screening and transitions from emergent axion-like dynamics

abstract: In [1], we demonstrated that the structure function $g_1(x_B,Q^2)$ measured in polarized deeply inelastic scattering (DIS) is dominated by the triangle anomaly in both the Bjorken limit of large $Q^2$ and the Regge limit of small $x_B$. In the worldline formulation of quantum field theory, the triangle anomaly arises from the imaginary part of the worldline effective action. We show explicitly how a Wess-Zumino-Witten term coupling the topological charge density to a primordial isosinglet ${\bar \eta}$ arises in this framework. We demonstrate the fundamental role played by this contribution both in topological mass generation of the $\eta^\prime$ and in the cancellation of the off-forward pole arising from the triangle anomaly in the proton's helicity $\Sigma(Q^2)$. We recover the striking result by Shore and Veneziano that $\Sigma\propto \sqrt{\chi'(0)}$, where $\chi'$ is the slope of the QCD topological susceptibility in the forward limit. We construct an axion-like effective action for $g_1$ at small $x_B$ that describes the interplay between gluon saturation and the topology of the QCD vacuum. In particular, we outline the role of "over-the-barrier" sphaleron-like transitions in spin diffusion at small $x_B$. Such topological transitions can be measured in polarized DIS at a future Electron-Ion Collider.


title: Polarization signatures of the head-on collision model for Type Ia supernovae: How much asymmetry is too much?

abstract: In a previous paper, we showed that the asymmetric ejecta produced by (zero impact parameter) head-on collisions of carbon-oxygen white dwarfs allow these progenitor models for Type Ia supernovae (SNe Ia) to cover the observed two-dimensional (2D) distribution of Si II line depths (Branch plot). In this paper, we study the polarization signature associated with the 2D asymmetric ejecta of the collision model and a double-detonation model using similar TARDIS radiative transfer simulations along different lines of sight with a spherical photosphere, combined with a new 3D Monte Carlo polarization code. We show that the polarization $Q$ can be parametrized as a product $Q=Q_{\max}Q_{\rm{x}}$ of a radial structure component $Q_{\max}$ which is insensitive to the model specifics and is shown to be universally around $Q_{\max}\sim 5\%$, and a cancellation component $Q_{\rm{x}}$ which depends on the asymmetry details. The continuum polarization is found to be low for both the collision and double-detonation models with $Q\sim 0.5\%$. However, the irregular Si distribution in the 2D head-on collision model results in Si II line polarization reaching $Q\sim 3\%$ ($Q_{\rm{x}} \lesssim 50\%$) in tension with observations (mostly $\lesssim 1.2\%$). In contrast, we show that the double-detonation model also covers the Branch plot, and yet results in low line polarization $Q\lesssim 0.7\%$ ($Q_{\rm{x}} \sim 10\%$) consistent with previous results and most SNe Ia. These results strengthen the case for asymmetric explosions as progenitors of SNe Ia, emphasizing an additional requirement for large polarization cancellations to account for the low observed line polarizations.


title: A discrete Darboux-Lax scheme for integrable difference equations

abstract: We propose a discrete Darboux-Lax scheme for deriving auto-B\"acklund transformations and constructing solutions to quad-graph equations that do not necessarily possess the 3D consistency property. As an illustrative example we use the Adler-Yamilov type system which is related to the nonlinear Schr\"odinger (NLS) equation [19]. In particular, we construct an auto-B\"acklund transformation for this discrete system, its superposition principle, and we employ them in the construction of the one- and two-soliton solutions of the Adler-Yamilov system.


title: Higher-Order Kinematical Effects in Deeply Virtual Compton Scattering

abstract: We study the deeply virtual Compton scattering cross-section in twist-two generalized parton distribution (GPD) approximation, and show that different choices of light-cone vectors and gauges for the final photon polarization will lead to different higher-order kinematical corrections to the cross-section formula. The choice of light-cone vectors affects kinematic corrections at the twist-three level, accounting for the differences between the cross-section formulas in the literature. On the other hand, kinematical corrections from higher-twist GPDs should eliminate the light-cone dependence at twist three. Those light-cone dependencies are studied systematically at JLab 12 GeV and future EIC kinematics. They serve as the intrinsic systematic uncertainties in extracting the Compton form factors through the cross-section formula. More importantly, they are also necessary for understanding cross-section measurements with higher-twist precision and to reconstruct higher-order Compton form factors.


title: Shear-driven solidification and nonlinear elasticity in epithelial tissues

abstract: Biological processes, from morphogenesis to tumor invasion, spontaneously generate shear stresses inside living tissue. The mechanisms that govern the transmission of mechanical forces in epithelia and the collective response of the tissue to bulk shear deformations remain, however, poorly understood. Using a minimal cell-based computational model, we investigate the constitutive relation of confluent tissues under simple shear deformation. We show that an initially undeformed fluidlike tissue acquires finite rigidity above a critical applied strain. This is akin to the shear-driven rigidity observed in other soft matter systems. Interestingly, shear-driven rigidity can be understood by a critical scaling analysis in the vicinity of the second order critical point that governs the liquid-solid transition of the undeformed system. We further show that a solidlike tissue responds linearly only to small strains and but then switches to a nonlinear response at larger stains, with substantial stiffening. Finally, we propose a mean-field formulation for cells under shear that offers a simple physical explanation of shear-driven rigidity and nonlinear response in a tissue.

categories: cond-mat.soft cond-mat.dis-nn physics.bio-ph q-bio.TO


title: AMBER: A Semi-numerical Abundance Matching Box for the Epoch of Reionization

abstract: The Abundance Matching Box for the Epoch of Reionization (AMBER) is a semi-numerical code for modeling the cosmic dawn. The new algorithm is not based on the excursion set formalism for reionization, but takes the novel approach of calculating the reionization-redshift field $z_\mathrm{re}(\boldsymbol{x})$ assuming that hydrogen gas encountering higher radiation intensity are photoionized earlier. Redshift values are assigned while matching the abundance of ionized mass according to a given mass-weighted ionization fraction $\bar{x}_\mathrm{i}(z)$. The code has the unique advantage of allowing users to directly specify the reionization history through the redshift midpoint $z_\mathrm{mid}$, duration $\Delta_\mathrm{z}$, and asymmetry $A_\mathrm{z}$ input parameters. The reionization process is further controlled through the minimum halo mass $M_\mathrm{min}$ for galaxy formation and the radiation mean free path $l_\mathrm{mfp}$ for radiative transfer. We implement improved methods for constructing density, velocity, halo, and radiation fields, which are essential components for modeling reionization observables. We compare AMBER with two other semi-numerical methods and find that our code more accurately reproduces the results from radiation-hydrodynamic simulations. The parallelized code is over four orders of magnitude faster than radiative transfer simulations and will efficiently enable large-volume models, full-sky mock observations, and parameter-space studies. AMBER will be made publicly available to facilitate and transform studies of the EoR.


title: Learning through structure: towards deep neuromorphic knowledge graph embeddings

abstract: Computing latent representations for graph-structured data is an ubiquitous learning task in many industrial and academic applications ranging from molecule synthetization to social network analysis and recommender systems. Knowledge graphs are among the most popular and widely used data representations related to the Semantic Web. Next to structuring factual knowledge in a machine-readable format, knowledge graphs serve as the backbone of many artificial intelligence applications and allow the ingestion of context information into various learning algorithms. Graph neural networks attempt to encode graph structures in low-dimensional vector spaces via a message passing heuristic between neighboring nodes. Over the recent years, a multitude of different graph neural network architectures demonstrated ground-breaking performances in many learning tasks. In this work, we propose a strategy to map deep graph learning architectures for knowledge graph reasoning to neuromorphic architectures. Based on the insight that randomly initialized and untrained (i.e., frozen) graph neural networks are able to preserve local graph structures, we compose a frozen neural network with shallow knowledge graph embedding models. We experimentally show that already on conventional computing hardware, this leads to a significant speedup and memory reduction while maintaining a competitive performance level. Moreover, we extend the frozen architecture to spiking neural networks, introducing a novel, event-based and highly sparse knowledge graph embedding algorithm that is suitable for implementation in neuromorphic hardware.


title: The colour matrix at next-to-leading-colour accuracy for tree-level multi-parton processes

abstract: We investigate the next-to-leading-colour (NLC) contributions to the colour matrix in the fundamental and the colour-flow decompositions for tree-level processes with all gluons, one quark pair and two quark pairs. By analytical examination of the colour factors, we find the non-zero elements in the colour matrix at NLC. At this colour order, together with the symmetry of the phase-space, it is reduced from factorial to polynomial the scaling of the contributing dual amplitudes as the number of partons participating in the scattering process is increased. This opens a path to an accurate tree-level matrix element generator of which all factorial complexity is removed, without resulting to Monte Carlo sampling over colour.


title: Normal, Dust-Obscured Galaxies in the Epoch of Reionization

abstract: Over the past decades, rest-frame ultraviolet (UV) observations have provided large samples of UV luminous galaxies at redshift (z) greater than 6, during the so-called epoch of reionization. While a few of these UV identified galaxies revealed significant dust reservoirs, very heavily dust-obscured sources at these early times have remained elusive. They are limited to a rare population of extreme starburst galaxies, and companions of rare quasars. These studies conclude that the contribution of dust-obscured galaxies to the cosmic star formation rate density at $z>6$ is sub-dominant. Recent ALMA and Spitzer observations have identified a more abundant, less extreme population of obscured galaxies at $z=3-6$. However, this population has not been confirmed in the reionization epoch so far. Here, we report the discovery of two dust-obscured star forming galaxies at $z=6.6813\pm0.0005$ and $z=7.3521\pm0.0005$. These objects are not detected in existing rest-frame UV data, and were only discovered through their far-infrared [CII] lines and dust continuum emission as companions to typical UV-luminous galaxies at the same redshift. The two galaxies exhibit lower infrared luminosities and star-formation rates than extreme starbursts, in line with typical star-forming galaxies at $z\sim7$. This population of heavily dust-obscured galaxies appears to contribute 10-25 per cent to the $z>6$ cosmic star formation rate density.


title: Scaling laws and a general theory for the growth of public companies

abstract: Publicly traded companies are fundamental units of contemporary economies and markets and are important mechanisms through which humans interact with their environments. Understanding the general properties that underlie the processes of their growth has long been of interest, yet fundamental debates about the effects of firm size on growth have persisted. Here we develop a scaling framework that focuses on company size as the critical feature determining a variety of tradeoffs, and use this to reveal novel systematic behavior across the diversity of publicly-traded companies. Using a large database of 31,553 US companies over nearly 70 years, and 3,160 Chinese companies over 24 year, we show how the dynamics of companies expressed as scaling relationships leads to a quantitative, analytic theory for their growth. This theory produces several predictions that are in good agreement with data for both the US and China, whose markets have strikingly different histories and underlying structures. In both cases sales scale sublinearly with assets and exhibit nearly identical exponents leading, surprisingly and nontrivially, to assets that grow as a power law in time rather than exponentially, as often assumed. On the other hand, liabilities scale linearly in the US (exponent of $1.0$) but superlinearly in China (exponent of $1.09$). We show that such small differences in scaling exponents can have a significant impact on the character and long-term evolution of growth trajectories. These results illustrate that while companies are part of a larger class of growth phenomena driven by incomes and costs that scale with size, they are unique in that they grow following a temporal power-function which sets them apart from organisms, cities, nations, and markets, whose growth over time is often exponential.


title: Deep Policies for Online Bipartite Matching: A Reinforcement Learning Approach

abstract: The challenge in the widely applicable online matching problem lies in making irrevocable assignments while there is uncertainty about future inputs. Most theoretically-grounded policies are myopic or greedy in nature. In real-world applications where the matching process is repeated on a regular basis, the underlying data distribution can be leveraged for better decision-making. We present an end-to-end Reinforcement Learning framework for deriving better matching policies based on trial-and-error on historical data. We devise a set of neural network architectures, design feature representations, and empirically evaluate them across two online matching problems: Edge-Weighted Online Bipartite Matching and Online Submodular Bipartite Matching. We show that most of the learning approaches perform consistently better than classical baseline algorithms on four synthetic and real-world datasets. On average, our proposed models improve the matching quality by 3--10\% on a variety of synthetic and real-world datasets. Our code is publicly available at https://github.com/lyeskhalil/CORL.


title: Improved multispecies Dougherty collisions

abstract: The Dougherty model Fokker-Planck operator is extended to describe nonlinear full-f collisions between multiple species in plasmas. Simple relations for cross-species interactions are developed which obey conservation laws, and reproduce familiar velocity and temperature relaxation rates. This treatment of multispecies Dougherty collisions, valid for arbitrary mass ratios, satisfies the H-Theorem unlike analogous Bhatnagar-Gross-Krook operators.


title: Weighed average over Narain moduli space as $T \bar T$ deformation of CFT target space

abstract: We consider the weighted average of a two dimensional CFT, whose target space is $T^2$, over its Narain moduli space. We take as the weighing function the integral kernel which gives rise to $T \bar T$ deformation when applied to the world sheet moduli data of the partition function viewed as vacuum amplitude when the world sheet is a torus. We compute the smeared partition function where this kernel is applied to the target space moduli. Smearing the partition function over the parameter space of a field theory generally leads to the breakdown in the ability to write the partition function as a sum over Boltzmann factor with unit coefficients. The weight function inspired by the $T \bar T$ deformation appears to be an exception to this general expectation. We show that this smearing leads to a marginal deformation corresponding to the overall rescaling of the target space $T^2$.


title: Generalized Orbital Angular Momentum Symmetry in Parametric Amplification

abstract: We investigate interesting symmetry properties verified by the down-converted beams produced in optical parametric amplification with structured light. We show that the Poincar\'e sphere symmetry, previously demonstrated for first-order spatial modes, translates to a multiple Poincar\'e sphere structure for higher orders. Each one of these multiple spheres is associated with a two-dimensional subspace defined by a different value of the orbital angular momentum. Therefore, the symmetry verified by first order modes is reproduced independently in each subspace. This effect can be useful for parallel control of independently correlated beams.


title: Quenching of star formation from a lack of inflowing gas to galaxies

abstract: Star formation in half of massive galaxies was quenched by the time the Universe was three billion years old. Very low amounts of molecular gas appear responsible for this, at least in some cases, though morphological gas stabilization, shock heating, or activity associated with accretion onto a central supermassive black hole is invoked in other cases. Recent studies of quenching by gas depletion have been based upon upper limits that are insufficiently sensitive to determine this robustly, or stacked emission with its problems of averaging. Here we report 1.3mm observations of dust emission from six strongly lensed galaxies where star formation has been quenched, with magnifications of up to a factor of 30. Four of the six galaxies are undetected in dust emission, with an estimated upper limit on the dust mass of 0.0001 times the stellar mass, and by proxy (assuming a Milky Way molecular gas-to-dust ratio) 0.01 times the stellar mass in molecular gas. This is two orders of magnitude less molecular gas per unit stellar mass than seen in star forming galaxies at similar redshifts. It remains difficult to extrapolate from these small samples, but these observations establish that gas depletion is responsible for a cessation of star formation in some fraction of high-redshift galaxies.


title: Learning to Guide Human Attention on Mobile Telepresence Robots with 360 Vision

abstract: Mobile telepresence robots (MTRs) allow people to navigate and interact with a remote environment that is in a place other than the person's true location. Thanks to the recent advances in 360 degree vision, many MTRs are now equipped with an all-degree visual perception capability. However, people's visual field horizontally spans only about 120 degree of the visual field captured by the robot. To bridge this observability gap toward human-MTR shared autonomy, we have developed a framework, called GHAL360, to enable the MTR to learn a goal-oriented policy from reinforcements for guiding human attention using visual indicators. Three telepresence environments were constructed using datasets that are extracted from Matterport3D and collected from a real robot respectively. Experimental results show that GHAL360 outperformed the baselines from the literature in the efficiency of a human-MTR team completing target search tasks.


title: Monotonicity for continuous-time random walks

abstract: Consider continuous-time random walks on Cayley graphs where the rate assigned to each edge depends only on the corresponding generator. We show that the limiting speed is monotone increasing in the rates for infinite Cayley graphs that arise from Coxeter systems, but not for all Cayley graphs. On finite Cayley graphs, we show that the distance -- in various senses -- to stationarity is monotone decreasing in the rates for Coxeter systems and for abelian groups, but not for all Cayley graphs. We also find several examples of surprising behaviour in the dependence of the distance to stationarity on the rates. This includes a counterexample to a conjecture on entropy of Benjamini, Lyons, and Schramm. We also show that the expected distance at any fixed time for random walks on $\mathbb{Z}^+$ is monotone increasing in the rates for arbitrary rate functions, which is not true on all of $\mathbb{Z}$. Various intermediate results are also of interest.


title: Toward Reusable Science with Readable Code and Reproducibility

abstract: An essential part of research and scientific communication is researchers' ability to reproduce the results of others. While there have been increasing standards for authors to make data and code available, many of these files are hard to re-execute in practice, leading to a lack of research reproducibility. This poses a major problem for students and researchers in the same field who cannot leverage the previously published findings for study or further inquiry. To address this, we propose an open-source platform named RE3 that helps improve the reproducibility and readability of research projects involving R code. Our platform incorporates assessing code readability with a machine learning model trained on a code readability survey and an automatic containerization service that executes code files and warns users of reproducibility errors. This process helps ensure the reproducibility and readability of projects and therefore fast-track their verification and reuse.


title: 4D Einstein-Gauss-Bonnet Gravity From Non-Einsteinian Phase

abstract: We set up an Einstein-Gauss-Bonnet theory in four dimensions, based on the recent formulation of pure gravity with extra dimensions of vanishing metrical length [1]. In absence of torsion, the effective field equations depend only on the four-metric, reflecting a quadratic curvature nonlinearity but no higher than second order derivatives. In contrast with recent proposals to obtain 4dEGB theories through a singular rescaling of the Gauss-Bonnet coupling, this formalism requires no (classical) regularization of divergences and is inequivalent to Horndeski gravity, while being generally covariant and independent of compactification. Notably, the vacuum field equations admit FLRW cosmologies containing nonsingular bounce and self-accelerating Universe, and spherically symmetric black holes more general than Schwarzschild. The quadratic theory of gravity emerging here is shown to be preserved against the inclusion of higher Lovelock densities in 5+2n dimensions.


title: Monopole-like configurations in the O(3) spin model at the upper critical dimension

abstract: We present a high-precision Monte Carlo study of the O(3) spin theory on the lattice in four dimensions. This model exhibits interesting dynamical features, in particular in the broken-symmetry phase, where suitable boundary conditions can be used to enforce monopole-like topological excitations. We investigate the Euclidean time propagation and the features of these excitations close to the critical point, where our numerical results show an excellent quantitative agreement with analytic predictions derived from purely quantum-field-theoretical tools by G. Delfino. We conclude by commenting on the implications of our findings for a conjectured violation of Derrick's theorem at the quantum level and on the consequences in various areas of physics, ranging from condensed matter to astro-particle physics.


title: Coast Sargassum Level Estimation from Smartphone Pictures

abstract: Since 2011, significant and atypical arrival of two species of surface dwelling algae, Sargassum natans and Sargassum Fluitans, have been detected in the Mexican Caribbean. This massive accumulation of algae has had a great environmental and economic impact. Therefore, for the government, ecologists, and local businesses, it is important to keep track of the amount of sargassum that arrives on the Caribbean coast. High-resolution satellite imagery is expensive or may be time delayed. Therefore, we propose to estimate the amount of sargassum based on ground-level smartphone photographs. From the computer vision perspective, the problem is quite difficult since no information about the 3D world is provided, in consequence, we have to model it as a classification problem, where a set of five labels define the amount. For this purpose, we have built a dataset with more than one thousand examples from public forums such as Facebook or Instagram and we have tested several state-of-the-art convolutional networks. As a result, the VGG network trained under fine-tuning showed the best performance. Even though the reached accuracy could be improved with more examples, the current prediction distribution is narrow, so the predictions are adequate for keeping a record and taking quick ecological actions.


title: Coalgebra symmetry for discrete systems

abstract: In this paper we introduce the notion of coalgebra symmetry for discrete systems. With this concept we prove that all discrete radially symmetric systems in standard form are quasi-integrable and that all variational discrete quasi-radially symmetric systems in standard form are Poincar\'e--Lyapunov--Nekhoroshev maps of order $N-2$, where $N$ are the degrees of freedom of the system. We also discuss the integrability properties of several vector systems which are generalisations of well-known one degree of freedom discrete integrable systems, including two $N$ degrees of freedom autonomous discrete Painlev\'e I equations and an $N$ degrees of freedom McMillan map.


title: Multi-Modal Model Predictive Control through Batch Non-Holonomic Trajectory Optimization: Application to Highway Driving

abstract: Standard Model Predictive Control (MPC) or trajectory optimization approaches perform only a local search to solve a complex non-convex optimization problem. As a result, they cannot capture the multi-modal characteristic of human driving. A global optimizer can be a potential solution but is computationally intractable in a real-time setting. In this paper, we present a real-time MPC capable of searching over different driving modalities. Our basic idea is simple: we run several goal-directed parallel trajectory optimizations and score the resulting trajectories based on user-defined meta cost functions. This allows us to perform a global search over several locally optimal motion plans. Although conceptually straightforward, realizing this idea in real-time with existing optimizers is highly challenging from technical and computational standpoints. With this motivation, we present a novel batch non-holonomic trajectory optimization whose underlying matrix algebra is easily parallelizable across problem instances and reduces to computing large batch matrix-vector products. This structure, in turn, is achieved by deriving a linearization-free multi-convex reformulation of the non-holonomic kinematics and collision avoidance constraints. We extensively validate our approach using both synthetic and real data sets (NGSIM) of traffic scenarios. We highlight how our algorithm automatically takes lane-change and overtaking decisions based on the defined meta cost function. Our batch optimizer achieves trajectories with lower meta cost, up to 6x faster than competing baselines.