title: Audio Watermarking over the Air With Modulated Self-Correlation

abstract: We propose a novel audio watermarking system that is robust to the distortion due to the indoor acoustic propagation channel between the loudspeaker and the receiving microphone. The system utilizes a set of new algorithms that effectively mitigate the impact of room reverberation and interfering sound sources without using dereverberation procedures. The decoder has low-latency and it operates asynchronously, which alleviates the need for explicit synchronization with the encoder. It is also robust to standard audio processing operations in legacy watermarking systems, e.g., compression and volume change. The effectiveness of the system is established with a real-time system under general room conditions.


title: Fully resonant scalars on asymptotically AdS wormholes

abstract: In this work we show the existence of asymptotically AdS wormhole geometries where the scalar probe has an equispaced, fully resonant spectrum, as that of a scalar on AdS spacetime, and explore its dynamics when non-linearities are included. The spacetime is a solution of Einstein-Gauss-Bonnet theory with a single maximally symmetric vacuum. Introducing a non-minimal coupling between the scalar probe and the Ricci scalar remarkably leads to a fully resonant spectrum for a scalar field fulfilling reflective boundary conditions at both infinities. Applying perturbative methods, which are particularly useful for unveiling the dynamics at time scales of order $\varepsilon^{-2}$ (where $\varepsilon$ characterizes the amplitude of the initial perturbation), we observe both direct and inverse energy cascades between modes. This motivates us to explore the energy returns in the case in which the dynamics is dominated by a single mode. We find numerical and perturbative evidence that near exact returns do exist in this regime. We also provide some comments on the fully backreracting case and provide a proof of the universality of the weakly non-linear dynamics around AdS, in the context of Lovelock theories with generic couplings, up to times of order $\varepsilon^{-2}$.


title: Dynamics of interacting dark soliton stripes

abstract: In the present work we examine the statics and dynamics of multiple parallel dark soliton stripes in a two-dimensional Bose-Einstein condensate. Our principal goal is to study the effect of the interaction between the stripes on the transverse instability of the individual stripes. We use a recently developed adiabatic invariant formulation to derive a quasi-analytical prediction for the stripe equilibrium position and for the Bogoliubov-de Gennes spectrum of excitations of stationary stripes. The cases of two-, three- and four-stripe states are studied in detail. We subsequently test our predictions against numerical simulations of the full two-dimensional Gross-Pitaevskii equation. We find that the number of unstable eigenmodes increases as the number of stripes increases due to (unstable) relative motions between the stripes. Their corresponding growth rates do not significantly change, although for large chemical potentials, the larger the stripe number, the larger the maximal instability growth rate. The instability induced dynamics of multiple stripe states and their decay into vortices are also investigated.


title: Fits to Non-Supersymmetric SO(10) Models with Type I and II Seesaw Mechanisms Using Renormalization Group Evolution

abstract: We consider numerical fits to non-supersymmetric $\mathrm{SO}(10)$-based models in which neutrino mass is generated by the type-I or type-II seesaw mechanism or a combination of both. The fits are performed with a sophisticated top-down procedure, taking into account the renormalization group equations of the gauge and Yukawa couplings, integrating out relevant degrees of freedom at their corresponding mass scales, and using recent data for the Standard Model observables. We find acceptable fits for normal neutrino mass ordering only and with neutrino mass generated by either type-I seesaw only or a combination of types I and II seesaw in which type-I seesaw is dominant. Furthermore, we find predictions from the best fit regarding the small neutrino masses, the effective neutrinoless double beta decay mass, and the leptonic CP-violating phase. Finally, we show that the fits are rather insensitive to the chosen value of the unification scale.


title: An assessment of multicomponent flow models and interface capturing schemes for spherical bubble dynamics

abstract: Numerical simulation of bubble dynamics and cavitation is challenging; even the seemingly simple problem of a collapsing spherical bubble is difficult to compute accurately with a general, three-dimensional, compressible, multicomponent flow solver. Difficulties arise due to both the physical model and the numerical method chosen for its solution. We consider the 5-equation model of Allaire et al. [1], the 5-equation model of Kapila et al. [2], and the 6-equation model of Saurel et al. [3] as candidate approaches for spherical bubble dynamics, and both MUSCL and WENO interface-capturing methods are implemented and compared. We demonstrate the inadequacy of the traditional 5-equation model of Allaire et al. [1] for spherical bubble collapse problems and explain the corresponding advantages of the augmented model of Kapila et al. [2] for representing this phenomenon. Quantitative comparisons between the augmented 5-equation and 6-equation models for three-dimensional bubble collapse problems demonstrate the versatility of pressure-disequilibrium models. Lastly, the performance of pressure disequilibrium model for representing a three-dimensional spherical bubble collapse for different bubble interior/exterior pressure ratios is evaluated for different numerical methods. Pathologies associated with each factor and their origins are identified and discussed.


title: A study of vectorization for matrix-free finite element methods

abstract: Vectorization is increasingly important to achieve high performance on modern hardware with SIMD instructions. Assembly of matrices and vectors in the finite element method, which is characterized by iterating a local assembly kernel over unstructured meshes, poses difficulties to effective vectorization. Maintaining a user-friendly high-level interface with a suitable degree of abstraction while generating efficient, vectorized code for the finite element method is a challenge for numerical software systems and libraries. In this work, we study cross-element vectorization in the finite element framework Firedrake via code transformation and demonstrate the efficacy of such an approach by evaluating a wide range of matrix-free operators spanning different polynomial degrees and discretizations on two recent CPUs using three mainstream compilers. Our experiments show that our approaches for cross-element vectorization achieve 30\% of theoretical peak performance for many examples of practical significance, and exceed 50\% for cases with high arithmetic intensities, with consistent speed-up over (intra-element) vectorization restricted to the local assembly kernels.


title: Singularities and unsteady separation for the inviscid two-dimensional Prandtl's system

abstract: We consider the inviscid unsteady Prandtl system in two dimensions, motivated by the fact that it should model to leading order separation and singularity formation for the original viscous system. We give a sharp expression for the maximal time of existence of regular solutions, showing that singularities only happen at the boundary or on the set of zero vorticity, and that they correspond to boundary layer separation. We then exhibit new Lagrangian formulae for backward self-similar profiles, and study them also with a different approach that was initiated by Elliott-Smith-Cowley and Cassel-Smith-Walker. One particular profile is at the heart of the so-called Van-Dommelen and Shen singularity, and we prove its generic appearance (that is, for an open and dense set of blow-up solutions) for any prescribed Eulerian outer flow. We comment on the connexion between these results and the full viscous Prandtl system. This paper combines ideas for transport equations, such as Lagrangian coordinates and incompressibility, and for singularity formation, such as self-similarity and renormalisation, in a novel manner, and designs a new way to study singularities for quasilinear transport equations.


title: Structural stability of shock waves in 2D compressible elastodynamics

abstract: We study the two-dimensional structural stability of shock waves in a compressible isentropic inviscid elastic fluid in the sense of the local-in-time existence and uniqueness of discontinuous shock front solutions of the equations of compressible elastodynamics in two space dimensions. By the energy method based on a symmetrization of the wave equation and giving an a priori estimate without loss of derivatives for solutions of the constant coefficients linearized problem we find a condition sufficient for the uniform stability of rectilinear shock waves. Comparing this condition with that for the uniform stability of shock waves in isentropic gas dynamics, we make the conclusion that the elastic force plays stabilizing role. In particular, we show that, as in isentropic gas dynamics, all compressive shock waves are uniformly stable for convex equations of state. Moreover, for some particular deformations (and general equations of state), by the direct test of the uniform Kreiss-Lopatinski condition we show that the stability condition found by the energy method is not only sufficient but also necessary for uniform stability. As is known, uniform stability implies structural stability of corresponding curved shock waves.


title: Equivariant Steinberg Summands

abstract: We construct Steinberg summands of $G$-equivariant spectra with $\mathrm{GL}_n(\mathbb{F}_p)$-action. We prove a lemma about their fixed points when $G$ is a $p$-group, and then use this lemma to compute the fixed points of the Steinberg summand of the equivariant classifying space of $(\mathbb{Z}/p)^n$. These results will be used in a companion paper to study the layers in the mod $p$ symmetric power filtration for $H\underline{\mathbb{F}}_p$.


title: The Average-Case Complexity of Counting Cliques in Erdos-Renyi Hypergraphs

abstract: We consider the problem of counting $k$-cliques in $s$-uniform Erdos-Renyi hypergraphs $G(n,c,s)$ with edge density $c$, and show that its fine-grained average-case complexity can be based on its worst-case complexity. We prove the following: 1. Dense Erdos-Renyi graphs and hypergraphs: Counting $k$-cliques on $G(n,c,s)$ with $k$ and $c$ constant matches its worst-case time complexity up to a $\mathrm{polylog}(n)$ factor. Assuming randomized ETH, it takes $n^{\Omega(k)}$ time to count $k$-cliques in $G(n,c,s)$ if $k$ and $c$ are constant. 2. Sparse Erdos-Renyi graphs and hypergraphs: When $c = \Theta(n^{-\alpha})$, we give several algorithms exploiting the sparsity of $G(n, c, s)$ that are faster than the best known worst-case algorithms. Complementing this, based on a fine-grained worst-case assumption, our results imply a different average-case phase diagram for each fixed $\alpha$ depicting a tradeoff between a runtime lower bound and $k$. Surprisingly, in the hypergraph case ($s \ge 3$), these lower bounds are tight against our algorithms exactly when $c$ is above the Erd\H{o}s-R\'{e}nyi $k$-clique percolation threshold. This is the first worst-case-to-average-case hardness reduction for a problem on Erd\H{o}s-R\'{e}nyi hypergraphs that we are aware of. We also give a variant of our result for computing the parity of the $k$-clique count that tolerates higher error probability.


title: Computational Tactile Flow for Anthropomorphic Grippers

abstract: Grasping objects requires tight integration between visual and tactile feedback. However, there is an inherent difference in the scale at which both these input modalities operate. It is thus necessary to be able to analyze tactile feedback in isolation in order to gain information about the surface the end-effector is operating on, such that more fine-grained features may be extracted from the surroundings. For tactile perception of the robot, inspired by the concept of the tactile flow in humans, we present the computational tactile flow to improve the analysis of the tactile feedback in robots using a Shadow Dexterous Hand. In the computational tactile flow model, given a sequence of pressure values from the tactile sensors, we define a virtual surface for the pressure values and define the tactile flow as the optical flow of this surface. We provide case studies that demonstrate how the computational tactile flow maps reveal information on the direction of motion and 3D structure of the surface, and feedback regarding the action being performed by the robot.


title: Etude d'un dispositif d'aide a l'intention d'{\'e}l{\`e}ves en difficult{\'e} dans la r{\'e}solution d'une situation-probl{\`e}me math{\'e}matique

abstract: In this paper, we analyze an aid session tested by an elementary school teachers. This aid session has been set up by a teacher for some students with difficulties after the work in the whole class. We first show how this aid session can help the pupils, especially through its chrono-topo-meso-genetic functions. Then we analyze a second situation in which the aid session is implemented, looking for clues of these functions in some episodes involving students who participated in the aid session.


title: The vacuum as a Lagrangian subspace

abstract: We unify and generalize the notions of vacuum and amplitude in linear quantum field theory in curved spacetime. Crucially, the generalized notion admits a localization in spacetime regions and on hypersurfaces. The underlying concept is that of a Lagrangian subspace of the space of complexified germs of solutions of the equations of motion on hypersurfaces. Traditional vacua and traditional amplitudes correspond to the special cases of definite and real Lagrangian subspaces respectively. Further, we introduce both infinitesimal and asymptotic methods for vacuum selection that involve a localized version of Wick rotation. We provide examples from Klein-Gordon theory in settings involving different types of regions and hypersurfaces to showcase generalized vacua and the application of the proposed vacuum selection methods. A recurrent theme is the occurrence of mixed vacua, where propagating solutions yield definite Lagrangian subspaces and evanescent solutions yield real Lagrangian subspaces. The examples cover Minkowski space, Rindler space, Euclidean space and de Sitter space. A simple formula allows for the calculation of expectation values for observables in the generalized vacua.


title: Asymptotically flat wormhole solutions with variable equation-of-state parameter

abstract: In this paper, we study exact wormhole solutions in the framework of general relativity with a general equation of state that reduced to a linear equation of state asymptotically. By considering a special shape function, we find classes of solutions which are asymptotically flat. We study the violation of NEC as the main ingredient in the wormhole physics. We investigate the possibility of finding wormhole solutions with asymptotically different state parameter. We show that in principle, wormhole with vanishing redshift function and the selected shape function, cannot satisfy NEC at large infinity. We present solutions which have the positive total amount of mater in the "volume integral quantifier" method. For this class of solutions, fluid near the wormhole throat is in the phantom regime and at some $r=r_{2}$, the phantom regime is connected to a dark energy regime. Thus, we need small amount of exotic matter to construct wormhole solutions.


title: MP net as Abstract Model of Communication for Message-passing Applications

abstract: MP net is a formal model specifically designed for the field of parallel applications that use a message passing interface. The main idea is to use MP net as a comprehensible way of presenting the actual structure of communication within MPI applications. The goal is to provide users with the kind of feedback that can help them to check quickly whether or not the actual communication within their application corresponds to the intended one. This paper introduces MP net that focuses on the communication part of parallel applications and emphasizes its spatial character, which is rather hidden in sequential (textual) form.


title: Fluidic Fabric Muscle Sheets for Wearable and Soft Robotics

abstract: Conformable robotic systems are attractive for applications in which they can be used to actuate structures with large surface areas, to provide forces through wearable garments, or to realize autonomous robotic systems. We present a new family of soft actuators that we refer to as Fluidic Fabric Muscle Sheets (FFMS). They are composite fabric structures that integrate fluidic transmissions based on arrays of elastic tubes. These sheet-like actuators can strain, squeeze, bend, and conform to hard or soft objects of arbitrary shapes or sizes, including the human body. We show how to design and fabricate FFMS actuators via facile apparel engineering methods, including computerized sewing techniques. Together, these determine the distributions of stresses and strains that can be generated by the FFMS. We present a simple mathematical model that proves effective for predicting their performance. FFMS can operate at frequencies of 5 Hertz or more, achieve engineering strains exceeding 100%, and exert forces greater than 115 times their own weight. They can be safely used in intimate contact with the human body even when delivering stresses exceeding 10$^\text{6}$ Pascals. We demonstrate their versatility for actuating a variety of bodies or structures, and in configurations that perform multi-axis actuation, including bending and shape change. As we also show, FFMS can be used to exert forces on body tissues for wearable and biomedical applications. We demonstrate several potential use cases, including a miniature steerable robot, a glove for grasp assistance, garments for applying compression to the extremities, and devices for actuating small body regions or tissues via localized skin stretch.


title: Efficient Off-Policy Meta-Reinforcement Learning via Probabilistic Context Variables

abstract: Deep reinforcement learning algorithms require large amounts of experience to learn an individual task. While in principle meta-reinforcement learning (meta-RL) algorithms enable agents to learn new skills from small amounts of experience, several major challenges preclude their practicality. Current methods rely heavily on on-policy experience, limiting their sample efficiency. The also lack mechanisms to reason about task uncertainty when adapting to new tasks, limiting their effectiveness in sparse reward problems. In this paper, we address these challenges by developing an off-policy meta-RL algorithm that disentangles task inference and control. In our approach, we perform online probabilistic filtering of latent task variables to infer how to solve a new task from small amounts of experience. This probabilistic interpretation enables posterior sampling for structured and efficient exploration. We demonstrate how to integrate these task variables with off-policy RL algorithms to achieve both meta-training and adaptation efficiency. Our method outperforms prior algorithms in sample efficiency by 20-100X as well as in asymptotic performance on several meta-RL benchmarks.


title: Markov Chain Models of Refugee Migration Data

abstract: The application of Markov chains to modelling refugee crises is explored, focusing on local migration of individuals at the level of cities and days. As an explicit example we apply the Markov chains migration model developed here to UNHCR data on the Burundi refugee crisis. We compare our method to a state-of-the-art `agent-based' model of Burundi refugee movements, and highlight that Markov chain approaches presented here can improve the match to data while simultaneously being more algorithmically efficient.


title: A Quantum Annealing-Based Approach to Extreme Clustering

abstract: Clustering, or grouping, dataset elements based on similarity can be used not only to classify a dataset into a few categories, but also to approximate it by a relatively large number of representative elements. In the latter scenario, referred to as extreme clustering, datasets are enormous and the number of representative clusters is large. We have devised a distributed method that can efficiently solve extreme clustering problems using quantum annealing. We prove that this method yields optimal clustering assignments under a separability assumption, and show that the generated clustering assignments are of comparable quality to those of assignments generated by common clustering algorithms, yet can be obtained a full order of magnitude faster.


title: Quantum Circuits for Scale Transformations

abstract: In this note we develop quantum circuits for exactly simulating the thermal properties of the quantum XY/Ising chain. These circuits are applicable to the simplest integrable lattice models for which the exact momentum-space scattering-matrix can be constructed in terms of local two-body gates. We discuss a circuit analogy for Wilsonian renormalization group flow and show how it can be used to simulate a global coarse graining transformation. We then construct a circuit that transforms the ground state of these models into the thermo-field double state, and show how it can be used to implement a local scale transformation.


title: Pulsar timing in extreme mass ratio binaries: a general relativistic approach

abstract: The detection of a pulsar (PSR) in a tight, relativistic orbit around a supermassive or intermediate mass black hole - such as those in the Galactic centre or in the centre of Globular clusters - would allow for precision tests of general relativity (GR) in the strong-field, non-linear regime. We present a framework for calculating the theoretical time-frequency signal from a PSR in such an Extreme Mass Ratio Binary (EMRB). This framework is entirely relativistic with no weak-field approximations and so able to account for all higher-order strong-field gravitational effects, relativistic spin dynamics, the convolution with astrophysical effects and the combined impact on the PSR timing signal. Specifically we calculate both the spacetime path of the pulsar radio signal and the complex orbital and spin dynamics of a spinning pulsar around a Kerr black hole, accounting for spacetime curvature and frame dragging, relativistic and gravitational time delay, gravitational light bending, temporal and spatial dispersion induced by the presence of plasma along the line of sight and relativistic aberration. This then allows for a consistent time-frequency solution to be generated. Such a framework is key for assessing the use of PSR as probes of strong field GR, helping to inform the detection of an EMRB system hosting a PSR and, most essentially, for providing an accurate theoretical basis to then compare with observations to test fundamental physics.


title: Spectral Properties of Laplacians on Snowflake domains and filled Julia sets

abstract: We present eigenvalue data and pictures of eigenfunctions of the classic and quadratic snowflake fractal and of quadratic filled Julia sets. Furthermore, we approximate the area and box-counting dimension of selected Julia sets to compare the eigenvalue counting function with the Weyl term.


title: Massive MIMO Optimization with Compatible Sets

abstract: Massive multiple-input multiple-output (MIMO) is expected to be a vital component in future 5G systems. As such, there is a need for new modeling in order to investigate the performance of massive MIMO not only at the physical layer, but also higher up the networking stack. In this paper, we present general optimization models for massive MIMO, based on mixed-integer programming and compatible sets, with both maximum ratio combing and zero forcing precoding schemes. We then apply our models to the case of joint device scheduling and power control for heterogeneous devices and traffic demands, in contrast to existing power control schemes that consider only homogeneous users and saturated scenarios. Our results show substantial benefits in terms of energy usage can be achieved without sacrificing throughput, and that both signalling overhead and the complexity of end devices can be reduced by abrogating the need for uplink power control through efficient scheduling.


title: Next-to-Leading Order Ab Initio Electron-Phonon Scattering

abstract: Electron-phonon (e-ph) interactions are usually treated in the lowest order of perturbation theory. Here we derive next-to-leading order e-ph interactions, and compute from first principles the associated two-phonon e-ph scattering rates. The derivation involves Matsubara sums of the relevant two-loop Feynman diagrams, and the numerical calculations are challenging since they involve Brillouin zone integrals over two crystal momenta and depend critically on the intermediate state lifetimes. Using random grids and Monte Carlo integration, together with a self-consistent update of the intermediate state lifetimes, we compute and converge the two-phonon scattering rates, using GaAs as a case study. For the longitudinal optical phonon in GaAs, we find that the two-phonon scattering rates are as large as nearly half the value of the leading-order rates. The energy and temperature dependence of the two-phonon processes are analyzed. We show that including the two-phonon processes is important to accurately predicting the electron mobility in GaAs.


title: Precursors to Exciton Condensation in Quantum Hall Bilayers

abstract: Tunneling spectroscopy reveals evidence for interlayer electron-hole correlations in quantum Hall bilayer two-dimensional electron systems at layer separations near, but above, the transition to the incompressible exciton condensate at total Landau level filling $\nu_T=1$. These correlations are manifested by a non-linear suppression of the Coulomb pseudogap which inhibits low energy interlayer tunneling in weakly-coupled bilayers. The pseudogap suppression is strongest at $\nu_T=1$ and grows rapidly as the critical layer separation for exciton condensation is approached from above.


title: Faster Algorithms for the Geometric Transportation Problem

abstract: Let $R$ and $B$ be two point sets in $\mathbb{R}^d$, with $|R|+ |B| = n$ and where $d$ is a constant. Next, let $\lambda : R \cup B \to \mathbb{N}$ such that $\sum_{r \in R } \lambda(r) = \sum_{b \in B} \lambda(b)$ be demand functions over $R$ and $B$. Let $\|\cdot\|$ be a suitable distance function such as the $L_p$ distance. The transportation problem asks to find a map $\tau : R \times B \to \mathbb{N}$ such that $\sum_{b \in B}\tau(r,b) = \lambda(r)$, $\sum_{r \in R}\tau(r,b) = \lambda(b)$, and $\sum_{r \in R, b \in B} \tau(r,b) \|r-b\|$ is minimized. We present three new results for the transportation problem when $\|r-b\|$ is any $L_p$ metric: - For any constant $\varepsilon > 0$, an $O(n^{1+\varepsilon})$ expected time randomized algorithm that returns a transportation map with expected cost $O(\log^2(1/\varepsilon))$ times the optimal cost. - For any $\varepsilon > 0$, a $(1+\varepsilon)$-approximation in $O(n^{3/2}\varepsilon^{-d} \operatorname{polylog}(U) \operatorname{polylog}(n))$ time, where $U = \max_{p\in R\cup B} \lambda(p)$. - An exact strongly polynomial $O(n^2 \operatorname{polylog}n)$ time algorithm, for $d = 2$.


title: Sketching a Proof of the Maldacena Conjecture at Small Radius

abstract: At small AdS radius, the superstring on $AdS_5\times S^5$ was conjectured by Maldacena to be equivalent to ${\cal N}=4$ super-Yang-Mills at small `t Hooft coupling where thickened Feynman diagrams can be used to compute scattering amplitudes. It was previously shown that the pure spinor worldsheet action of the $AdS_5\times S^5$ superstring can be expressed as the sum of a BRST-trivial term and a `B-term' which is antisymmetric in worldsheet derivatives. Using the explicit form of the pure spinor vertex operators, it will be argued here that the free super-Yang-Mills Feynman diagrams are described by the BRST-trivial term where the thickened propagators are the regions of the string worldsheet near the AdS boundary and the holes are the regions near the AdS horizon. Evidence will then be presented that the antisymmetric B-term generates the super-Yang-Mills vertex so that, at small radius and arbitrary genus, the superstring amplitudes correctly reproduce the super-Yang-Mills Feynman diagram expansion.


title: Quadratic Gorenstein rings and the Koszul property I

abstract: Let $R$ be a standard graded Gorenstein algebra over a field presented by quadrics. Conca, Rossi, and Valla have shown that such a ring is Koszul if $\mathrm{reg}\, R \leq 2$ or if $\mathrm{reg}\, = 3$ and $c= \mathrm{codim}\, R \leq 4$, and they ask whether this is true for $\mathrm{reg}\, R = 3$ in general. We determine sufficient conditions on a non-Koszul quadratic Cohen-Macaulay ring $R$ that guarantee the Nagata idealization $\tilde{R} = R \ltimes \omega_R(-a-1)$ is a non-Koszul quadratic Gorenstein ring. We use this to negatively answer the question of Conca-Rossi-Valla, constructing non-Koszul quadratic Gorenstein rings of regularity 3 for all $c \geq 9$.


title: Caps and progression-free sets in $\mathbb{Z}_m^n$

abstract: We study progression-free sets in the abelian groups $G=(\mathbb{Z}_m^n,+)$. Let $r_k(\mathbb{Z}_m^n)$ denote the maximal size of a set $S \subset \mathbb{Z}_m^n$ that does not contain a proper arithmetic progression of length $k$. We give lower bound constructions, which e.g. include that $r_3(\mathbb{Z}_m^n) \geq C_m \frac{((m+2)/2)^n}{\sqrt{n}}$, when $m$ is even. When $m=4$ this is of order at least $3^n/\sqrt{n}\gg \vert G \vert^{0.7924}$. Moreover, if the progression-free set $S\subset \mathbb{Z}_4^n$ satisfies a technical condition, which dominates the problem at least in low dimension, then $|S|\leq 3^n$ holds. We present a number of new methods which cover lower bounds for several infinite families of parameters $m,k,n$, which includes for example: $r_6(\mathbb{Z}_{125}^n) \geq (85-o(1))^n$. For $r_3(\mathbb{Z}_4^n)$ we determine the exact values, when $n \leq 5$, e.g. $r_3(\mathbb{Z}_4^5)=124$, and for $r_4(\mathbb{Z}_4^n)$ we determine the exact values, when $n \leq 4$, e.g. $r_4(\mathbb{Z}_4^4)=128$.


title: The spectral function of Mott-insulating Hubbard ladders: From fractionalized excitations to coherent quasi-particles

abstract: We study the spectral function of two-leg Hubbard ladders with the time-dependent density matrix renormalization group method (tDMRG). The high-resolution spectrum displays features of spin-charge separation and a scattering continuum of excitations with coherent bands of bound states `leaking' from it. As the inter-leg hopping is increased, the continuum in the bonding channel moves to higher energies and spinon and holon branches merge into a single coherent quasi-particle band. Simultaneously, the spectrum undergoes a crossover from a regime with two minima at incommensurate values of $k_x$ (a Mott insulator), to one with a single minimum at $k_x=\pi$ (a band insulator). We identify the presence of a continuum of scattering states consisting of a triplon and a polaron. We analyze the processes leading to quasiparticle formation by studying the time evolution of charge and spin degrees of freedom in real space after the hole is created. At short times, incoherent holons and spinons are emitted but after a characteristic time $\tau$ charge and spin form polarons that propagate coherently.


title: Simple, Fast, Accurate Intent Classification and Slot Labeling for Goal-Oriented Dialogue Systems

abstract: With the advent of conversational assistants, like Amazon Alexa, Google Now, etc., dialogue systems are gaining a lot of traction, especially in industrial setting. These systems typically consist of Spoken Language understanding component which, in turn, consists of two tasks - Intent Classification (IC) and Slot Labeling (SL). Generally, these two tasks are modeled together jointly to achieve best performance. However, this joint modeling adds to model obfuscation. In this work, we first design framework for a modularization of joint IC-SL task to enhance architecture transparency. Then, we explore a number of self-attention, convolutional, and recurrent models, contributing a large-scale analysis of modeling paradigms for IC+SL across two datasets. Finally, using this framework, we propose a class of 'label-recurrent' models that otherwise non-recurrent, with a 10-dimensional representation of the label history, and show that our proposed systems are easy to interpret, highly accurate (achieving over 30% error reduction in SL over the state-of-the-art on the Snips dataset), as well as fast, at 2x the inference and 2/3 to 1/2 the training time of comparable recurrent models, thus giving an edge in critical real-world systems.


title: Nonlinear Inverse Spin Galvanic Effect in Anisotropic Disorder-free Systems

abstract: Spin transport phenomena in solid materials suffer limitations from spin relaxation associated to disorder or lack of translational invariance. Ultracold atoms, free of that disorder, can provide a platform to observe phenomena beyond the usual two-dimensional electron gas. By generalizing the approach used for isotropic two-dimensional electron gases, we theoretically investigate the inverse spin galvanic effect in the two-level atomic system in the presence of anisotropic Rashba-Dresselhaus spin-orbit couplings (SOC) and external magnetic field. We show that the combination of the SOC results in an asymmetric case: the total spin polarization considered for a small momentum has a longer spin state than in a two-dimensional electron gas when the SOC field prevails over the external electric field. Our results can be relevant for advancing experimental and theoretical investigations in spin dynamics as a basic approach for studying spin state control.


title: Minimizing polynomial functions on quantum computers

abstract: This expository paper reviews some of the recent uses of computational algebraic geometry in classical and quantum optimization. The paper assumes an elementary background in algebraic geometry and adiabatic quantum computing (AQC), and concentrates on presenting concrete examples (with Python codes tested on a quantum computer) of applying algebraic geometry constructs: solving binary optimization, factoring, and compiling. Reversing the direction, we also briefly describe a novel use of quantum computers to compute Groebner bases for toric ideals. We also show how Groebner bases play a role in studying AQC at a fundamental level within a Morse theory framework. We close by placing our work in perspective, by situating this leg of the journey, as part of a marvelous intellectual expedition that began with our ancients over 4000 years ago.


title: Exact Distributions of Currents and Frenesy for Markov Bridges

abstract: We consider discrete-time Markov bridges, chains whose initial and final states coincide. We derive exact finite-time formulae for the joint probability distributions of additive functionals of trajectories. We apply our theory to time-integrated currents and frenesy of enzymatic reactions, which may include absolutely irreversible transitions. We discuss the information that frenesy carries about the currents and show that bridges may violate known uncertainty relations in certain cases. Numerical simulations are in perfect agreement with our theory.

categories: cond-mat.stat-mech cond-mat.dis-nn physics.bio-ph physics.data-an


title: Design and Evaluation of Self-Assembled Actin-Based Nano-Communication

abstract: The tremendous progress in nanotechnology over the last century, makes it possible to engineer tiny nanodevices, which they need a nano-communication network to interact. Two solutions are proposed in literature to create a nano-communication system, either by using the classical electromagnetic paradigm with Terahertz band, or using the bio-inspired molecular communication. However, Terahertz is suffering from molecular absorption and scattering losses at nano level, and the achievable throughput of molecular communication is very low. In this paper, we propose a new solution to establish a wired nano-communication. Self-assembled actin-based is a new method that takes advantage of actin filaments self-assembly to create a nano wire between a transmitter and a receiver, and use electrons as information carriers. VPython framework is used in this paper to perform stochastic simulations of the nano wire formation. The algorithms used for the simulations are presented. The stability of the constructed nano wire is analyzed, and the error probability is calculated. Self-assembled actin-based method promises a fast and stable nano-communication system with a very high achievable throughput.


title: Quadratic Gorenstein rings and the Koszul property II

abstract: A question of Conca, Rossi, and Valla asks whether every quadratic Gorenstein ring $R$ of regularity three is Koszul. In a previous paper, we use idealization to answer their question, proving that in nine or more variables there exist quadratic Gorenstein rings of regularity three which are not Koszul. In this paper, we study the analog of the Conca-Rossi-Valla question when the regularity of $R$ is four or more. Let $R$ be a quadratic Gorenstein ring having $\mathrm{codim}\, R = c$ and $\mathrm{reg}\, R = r \ge 4$. We prove that if $c = r+1$ then $R$ is always Koszul, and for every $c \geq r+2$, we construct quadratic Gorenstein rings that are not Koszul, answering questions of Matsuda and Migliore-Nagel concerning the $h$-vectors of quadratic Gorenstein rings.


title: Slow Fibonacci Walks

abstract: For a positive integer $n$, we study the number of steps to reach $n$ by a {\it Fibonacci walk} for some starting pair $a_1$ and $a_2$ satisfying the recurrence of $a_{k+2}=a_{k+1}+a_k$. The problem of slow Fibonacci walks, first suggested by Richard Stanley, is to determine the maximum number $s(n)$ of steps for such a Fibonacci walk ending at $n$. Stanley conjectured that for most $n$, there is a slow Fibonacci walk reaching $n = a_s$ with the property that $a_{s+1}$ is the integer closest to $\phi n$ where $\phi=(1+\sqrt{5})/2$. We prove that this is true for only a positive fraction of $n$. We give explicit formulas for the choice of the starting pairs and the determination of $s(n)$ by giving a characterization theorem. We also derive a number of density results concerning the distribution of down and up cases (that is, those $n$ with $a_{s+1}=\lfloor \phi n\rfloor$ or $\lceil \phi n \rceil$, respectively), as well as for more general `paradoxical' cases.


title: Gelfand-Tsetlin polytopes: a story of flow and order polytopes

abstract: Gelfand-Tsetlin polytopes are prominent objects in algebraic combinatorics. The number of integer points of the Gelfand-Tsetlin polytope $\mathrm{GT}(\lambda)$ is equal to the dimension of the corresponding irreducible representation of $GL(n)$. It is well-known that the Gelfand-Tsetlin polytope is a marked order polytope; the authors have recently shown it to be a flow polytope. In this paper, we draw corollaries from this result and establish a general theory connecting marked order polytopes and flow polytopes.


title: Switching to nonhyperbolic cycles from codimension two bifurcations of equilibria of delay differential equations

abstract: In this paper we perform the parameter-dependent center manifold reduction near the generalized Hopf (Bautin), fold-Hopf, Hopf-Hopf and transcritical-Hopf bifurcations in delay differential equations (DDEs). This allows us to initialize the continuation of codimension one equilibria and cycle bifurcations emanating from these codimension two bifurcation points. The normal form coefficients are derived in the functional analytic perturbation framework for dual semigroups (sun-star calculus) using a normalization technique based on the Fredholm alternative. The obtained expressions give explicit formulas which have been implemented in the freely available numerical software package DDE-BifTool. While our theoretical results are proven to apply more generally, the software implementation and examples focus on DDEs with finitely many discrete delays. Together with the continuation capabilities of DDE-BifTool, this provides a powerful tool to study the dynamics near equilibria of such DDEs. The effectiveness is demonstrated on various models.


title: Almost dominant generalized slices and convolution diagrams over them

abstract: Let $G$ be a connected reductive complex algebraic group with a maximal torus $T$. We denote by $\Lambda$ the cocharacter lattice of $(T,G)$. Let $\Lambda^+ \subset \Lambda$ be the submonoid of dominant coweights. For $\lambda \in \Lambda^+,\,\mu \in \Lambda,\,\mu \leqslant \lambda$, in arXiv:1604.03625, authors defined a generalized transversal slice $\overline{\mathcal{W}}^\lambda_\mu$. This is an algebraic variety of the dimension $\langle 2\rho^{\vee}, \lambda-\mu \rangle$, where $2\rho^{\vee}$ is the sum of positive roots of $G$. In this paper, we construct an isomorphism $\overline{\mathcal{W}}^\lambda_\mu \simeq \overline{\mathcal{W}}^\lambda_{\mu^+} \times {\mathbb{A}}^{\langle 2\rho^{\vee},\, \mu^+-\mu\rangle}$ for $\mu \in \Lambda$ such that $\langle \alpha^{\vee},\mu\rangle \geqslant -1$ for any positive root $\alpha^{\vee}$, here $\mu^+ \in W\mu$ is the dominant representative in the Weyl group orbit of $\mu$. We consider the example when $\lambda$ is minuscule, $\mu \in W\lambda$ and describe natural coordinates, Poisson structure on $\overline{\mathcal{W}}^\lambda_\mu \simeq {\mathbb{A}}^{\langle 2\rho^\vee,\,\lambda-\mu \rangle}$ and its $T\times {\mathbb{C}}^\times$-character. We apply these results to compute $T \times {\mathbb{C}}^\times$-characters of tangent spaces at fixed points of convolution diagrams $\widetilde{\mathcal{W}}^{\underline{\lambda}}_\mu$ with minuscule $\lambda_i$. We also apply our results to construct open coverings by affine spaces of convolution diagrams $\widetilde{\mathcal{W}}^{\underline{\lambda}}_\mu$ over slices with $\mu$ such that $\langle \alpha^{\vee},\mu\rangle \geqslant -1$ for any positive root $\alpha^{\vee}$ and minuscule $\lambda_i$ and to compute Poincar\'e polynomials of such convolution diagrams $\widetilde{\mathcal{W}}^{\underline{\lambda}}_{\mu}$.


title: The icosahedra of edge length 1

abstract: Retaining the combinatorial Euclidean structure of a regular icosahedron, namely the 20 equiangular (planar) triangles, the 30 edges of length 1, and the 12 different vertices together with the incidence structure, we investigate variations of the regular icosahedron admitting self-intersections of faces. We determine all rigid equivalence classes of these icosahedra with non-trivial automorphism group and find one curve of flexible icosahedra. Visualisations and explicit data for this paper are available under http://algebra.data.rwth-aachen.de/Icosahedra/visualplusdata.html.


title: Quantum Key Distribution with no Shared Reference Frame

abstract: Any quantum communication task requires a common reference frame (i.e. phase, coordinate system). In particular, Quantum Key Distribution requires different bases for preparation and measurements of states which are obviously based on the existence of a common frame of reference. Here we show how QKD can be achieved in the absence of any common frame of reference. We study the coordinate reference frame, where the two parties do not even share a single direction, but the method can be generalized to other general frames of reference, pertaining to other groups of transformations.


title: Preprocessing Ambiguous Imprecise Points

abstract: Let ${R} = \{R_1, R_2, ..., R_n\}$ be a set of regions and let $ X = \{x_1, x_2, ..., x_n\}$ be an (unknown) point set with $x_i \in R_i$. Region $R_i$ represents the uncertainty region of $x_i$. We consider the following question: how fast can we establish order if we are allowed to preprocess the regions in $R$? The preprocessing model of uncertainty uses two consecutive phases: a preprocessing phase which has access only to ${R}$ followed by a reconstruction phase during which a desired structure on $X$ is computed. Recent results in this model parametrize the reconstruction time by the ply of ${R}$, which is the maximum overlap between the regions in ${R}$. We introduce the ambiguity $A({R})$ as a more fine-grained measure of the degree of overlap in ${R}$. We show how to preprocess a set of $d$-dimensional disks in $O(n \log n)$ time such that we can sort $X$ (if $d=1$) and reconstruct a quadtree on $X$ (if $d\geq 1$ but constant) in $O(A({R}))$ time. If $A({R})$ is sub-linear, then reporting the result dominates the running time of the reconstruction phase. However, we can still return a suitable data structure representing the result in $O(A({R}))$ time. In one dimension, ${R}$ is a set of intervals and the ambiguity is linked to interval entropy, which in turn relates to the well-studied problem of sorting under partial information. The number of comparisons necessary to find the linear order underlying a poset $P$ is lower-bounded by the graph entropy of $P$. We show that if $P$ is an interval order, then the ambiguity provides a constant-factor approximation of the graph entropy. This gives a lower bound of $\Omega(A({R}))$ in all dimensions for the reconstruction phase (sorting or any proximity structure), independent of any preprocessing; hence our result is tight.


title: Shrinking the Sample Covariance Matrix using Convex Penalties on the Matrix-Log Transformation

abstract: For $q$-dimensional data, penalized versions of the sample covariance matrix are important when the sample size is small or modest relative to $q$. Since the negative log-likelihood under multivariate normal sampling is convex in $\Sigma^{-1}$, the inverse of its covariance matrix, it is common to add to it a penalty which is also convex in $\Sigma^{-1}$. More recently, Deng-Tsui (2013) and Yu et al.(2017) have proposed penalties which are functions of the eigenvalues of $\Sigma$, and are convex in $\log \Sigma$, but not in $\Sigma^{-1}$. The resulting penalized optimization problem is not convex in either $\log \Sigma$ or $\Sigma^{-1}$. In this paper, we note that this optimization problem is geodesically convex in $\Sigma$, which allows us to establish the existence and uniqueness of the corresponding penalized covariance matrices. More generally, we show the equivalence of convexity in $\log \Sigma$ and geodesic convexity for penalties on $\Sigma$ which are strictly functions of their eigenvalues. In addition, when using such penalties, we show that the resulting optimization problem reduces to to a $q$-dimensional convex optimization problem on the eigenvalues of $\Sigma$, which can then be readily solved via Newton-Raphson. Finally, we argue that it is better to apply these penalties to the shape matrix $\Sigma/(\det \Sigma)^{1/q}$ rather than to $\Sigma$ itself. A simulation study and an example illustrate the advantages of applying the penalty to the shape matrix.


title: Attack-resilient Estimation for Linear Discrete-time Stochastic Systems with Input and State Constraints

abstract: In this paper, an attack-resilient estimation algorithm is presented for linear discrete-time stochastic systems with state and input constraints. It is shown that the state estimation errors of the proposed estimation algorithm are practically exponentially stable.


title: Protoplanetary Disk Masses from Radiative Transfer Modeling: A Case Study in Taurus

abstract: Measuring the masses of protoplanetary disks is crucial for understanding their planet-forming potential. Typically, dust masses are derived from (sub-)millimeter flux density measurements plus assumptions for the opacity, temperature, and optical depth of the dust. Here we use radiative transfer models to quantify the validity of these assumptions with the aim of improving the accuracy of disk dust mass measurements. We first carry out a controlled exploration of disk parameter space. We find that the disk temperature is a strong function of disk size, while the optical depth depends on both disk size and dust mass. The millimeter-wavelength spectral index can be significantly shallower than the naive expectation due to a combination of optical depth and deviations from the Rayleigh-Jeans regime. We fit radiative transfer models to the spectral energy distributions (SEDs) of 132 disks in the Taurus-Auriga region using a Markov chain Monte Carlo approach. We used all available data to produce the most complete SEDs used in any extant modeling study. We perform the fitting twice: first with unconstrained disk sizes and again imposing the disk size--brightness relation inferred for sources in Taurus. This constraint generally forces the disks to be smaller, warmer, and more optically thick. From both sets of fits, we find disks to be $\sim$1--5 times more massive than when derived using (sub-)millimeter measurements and common assumptions. With the uncertainties derived from our model fitting, the previously measured dust mass--stellar mass correlation is present in our study but only significant at the 2$\sigma$ level.


title: Black hole ringdown: the importance of overtones

abstract: It is possible to infer the mass and spin of the remnant black hole from binary black hole mergers by comparing the ringdown gravitational wave signal to results from studies of perturbed Kerr spacetimes. Typically these studies are based on the fundamental quasinormal mode of the dominant $\ell=m=2$ harmonic. By modeling the ringdown of accurate numerical relativity simulations, we find that the fundamental mode alone is insufficient to recover the true underlying mass and spin, unless the analysis is started very late in the ringdown. Including higher overtones associated with this $\ell=m=2$ harmonic resolves this issue, and provides an unbiased estimate of the true remnant parameters. Further, including overtones allows for the modeling of the ringdown signal for all times beyond the peak strain amplitude, indicating that the linear quasinormal regime starts much sooner than previously expected. A model for the ringdown beginning at the peak strain amplitude can exploit the higher signal-to-noise ratio in detectors, reducing uncertainties in the extracted remnant quantities. Tests of the no-hair theorem should consider incorporating overtones in the analysis.


title: Magnetic properties and domain structure of ultrathin yttrium iron garnet/Pt bilayers

abstract: We report on the structure, magnetization, magnetic anisotropy, and domain morphology of ultrathin yttrium iron garnet (YIG)/Pt films with thickness ranging from 3 to 90 nm. We find that the saturation magnetization is close to the bulk value in the thickest films and decreases towards low thickness with a strong reduction below 10 nm. We characterize the magnetic anisotropy by measuring the transverse spin Hall magnetoresistance as a function of applied field. Our results reveal strong easy plane anisotropy fields of the order of 50-100 mT, which add to the demagnetizing field, as well as weaker in-plane uniaxial anisotropy ranging from 10 to 100 $\mu$T. The in-plane easy axis direction changes with thickness, but presents also significant fluctuations among samples with the same thickness grown on the same substrate. X-ray photoelectron emission microscopy reveals the formation of zigzag magnetic domains in YIG films thicker than 10 nm, which have dimensions larger than several 100 $\mu$m and are separated by achiral N\'{e}el-type domain walls. Smaller domains characterized by interspersed elongated features are found in YIG films thinner than 10 nm.


title: An extension of the Glauberman ZJ-Theorem

abstract: Let $p$ be an odd prime and let $J_o(X)$, $J_r(X)$ and $J_e(X)$ denote the three different versions of Thompson subgroups for a $p$-group $X$. In this article, we first prove an extension of Glauberman's replacement theorem. Secondly, we prove the following: Let $G$ be a $p$-stable group and $P\in Syl_p(G)$. Suppose that $C_G(O_{p}(G))\leq O_{p}(G)$. If $D$ is a strongly closed subgroup in $P$, then $Z(J_o(D))$, $\Omega(Z(J_r(D)))$ and $\Omega(Z(J_e(D)))$ are normal subgroups of $G$. Thirdly, we show the following: Let $G$ be a $\text{Qd}(p)$-free group and $P\in Syl_p(G)$. If $D$ is a strongly closed subgroup in $P$, then the normalizers of the subgroups $Z(J_o(D))$, $\Omega(Z(J_r(D)))$ and $\Omega(Z(J_e(D)))$ control strong $G$-fusion in $P$. We also prove a similar result for a $p$-stable and $p$-constrained group. Lastly, we give a $p$-nilpotency criteria, which is an extension of Glauberman-Thompson $p$-nilpotency theorem.


title: Distinguishing Local and non-Local Demagnetization in Ferromagnetic FePt Nanoparticles

abstract: Time-resolved coherent X-ray diffraction is used to measure the spatially resolved magnetization structure within FePt nanoparticles during laser-induced ultrafast demagnetization. The momentum-dependent X-ray magnetic diffraction shows that demagnetization proceeds at different rates at different X-ray momentum transfer. We show that the observed momentum-dependent scattering has the signature of inhomogeneous demagnetization within the nanoparticles, with the demagnetization proceeding more rapidly at the boundary of the nanoparticle. A shell region of reduced magnetization forms and moves inwards at a supermagnonic velocity. Spin-transport calculations show that the shell formation is driven by superdiffusive spin flux mainly leaving the nanoparticle into the surrounding carbon. Quantifying this non-local contribution to the demagnetization allows us to separate it from the local demagnetization.


title: Presentations of Transversal Valuated Matroids

abstract: Given $d$ row vectors of $n$ tropical numbers, $d<n$, the tropical Stiefel map constructs a version of their row space, whose Pl\"ucker coordinates are tropical determinants. We explicitly describe the fibers of this map. From the viewpoint of matroid theory, the tropical Stiefel map defines a generalization of transversal matroids in the valuated context, and our results are the valuated generalizations of theorems of Brualdi and Dinolt, Mason and others on the set of all set families that present a given transversal matroid. We show that a connected valuated matroid is transversal if and only if all of its connected initial matroids are. The duals of our results describe complete stable intersections via valuated strict gammoids.


title: GANs for Semi-Supervised Opinion Spam Detection

abstract: Online reviews have become a vital source of information in purchasing a service (product). Opinion spammers manipulate reviews, affecting the overall perception of the service. A key challenge in detecting opinion spam is obtaining ground truth. Though there exists a large set of reviews online, only a few of them have been labeled spam or non-spam. In this paper, we propose spamGAN, a generative adversarial network which relies on limited set of labeled data as well as unlabeled data for opinion spam detection. spamGAN improves the state-of-the-art GAN based techniques for text classification. Experiments on TripAdvisor dataset show that spamGAN outperforms existing spam detection techniques when limited labeled data is used. Apart from detecting spam reviews, spamGAN can also generate reviews with reasonable perplexity.


title: Derivation of a homogenized bending--torsion theory for rods with micro-heterogeneous prestrain

abstract: In this paper we investigate rods made of nonlinearly elastic, composite--materials that feature a micro-heterogeneous prestrain that oscillates (locally periodic) on a scale that is small compared to the length of the rod. As a main result we derive a homogenized bending-torsion theory for rods as $\Gamma$-limit from 3D nonlinear elasticity by simultaneous homogenization and dimension reduction under the assumption that the prestrain is of the order of the diameter of the rod. The limit model features a spontaneous curvature-torsion tensor that captures the macroscopic effect of the micro-heterogeneous prestrain. We device a formula that allows to compute the spontaneous curvature-torsion tensor by means of a weighted average of the given prestrain. The weight in the average depends on the geometry of the composite and invokes correctors that are defined with help of boundary value problems for the system of linear elasticity. The definition of the correctors depends on a relative scaling parameter $\gamma$, which monitors the ratio between the diameter of the rod and the period of the composite's microstructure. We observe an interesting size-effect: For the same prestrain a transition from flat minimizers to curved minimizers occurs by just changing the value of $\gamma$. Moreover, in the paper we analytically investigate the microstructure-properties relation in the case of isotropic, layered composites, and consider applications to nematic liquid-crystal-elastomer rods and shape programming.


title: A simple model of Small-World Quantum Networks

abstract: A simple model of small world quantum networks, in which a central node plays essential role, is introduced for sharing entanglement over long distances. In view of the challenges in setting up advanced quantum labs which allows only few nodes in a network to play a central role, this kind of small world network may be more relevant for entanglement distribution. It is shown that by only adding a small number of short-cuts, it is possible to produce maximally entangled pairs between arbitrary nodes in the network. Besides this, the threshold values for the initial amount of entanglement and the distance between nodes, for obtaining a highly entangled states shared between remote points are also investigated.


title: Integral forms in vertex operator algebras, a survey

abstract: We give a brief survey of recent work on integral forms in vertex operator algebras (VOAs).


title: Circular orbit of a test particle and phase transition of a black hole

abstract: The radius of the circular orbit for the time-like or light-like test particle in a background of general spherically symmetric spacetime is viewed as a characterized quantity for the thermodynamic phase transition of the corresponding black hole. We generally show that the phase transition information of a black hole can be reflected by its surrounding particle's circular orbit.


title: Sundials in the Shade: An Internet-wide Perspective on ICMP Timestamps

abstract: ICMP timestamp request and response packets have been standardized for nearly 40 years, but have no modern practical application, having been superseded by NTP. However, ICMP timestamps are not deprecated, suggesting that while hosts must support them, little attention is paid to their implementation and use. In this work, we perform active measurements and find 2.2 million hosts on the Internet responding to ICMP timestamp requests from over 42,500 unique autonomous systems. We develop a methodology to classify timestamp responses, and find 13 distinct classes of behavior. Not only do these behaviors enable a new fingerprinting vector, some behaviors leak important information about the host e.g., OS, kernel version, and local timezone.


title: Algebraic $K$-theory of planar cuspidal curves

abstract: In this paper, we evaluate the algebraic $K$-groups of a planar cuspidal curve over a perfect $\mathbb{F}_p$-algebra relative to the cusp point. A conditional calculation of these groups was given earlier by Hesselholt, assuming a conjecture on the structure of certain polytopes. Our calculation here, however, is unconditional and illustrates the advantage of the new setup for topological cyclic homology by Nikolaus-Scholze, which is used throughout. The only input necessary for our calculation is the evaluation by the Buenos Aires Cyclic Homology group and by Larsen of the structure of Hochschild complex of the coordinate ring as a mixed complex, that is, as an object of the infinity category of chain complexes with circle action.


title: Self-interacting sterile neutrino dark matter: the heavy-mediator case

abstract: For active-sterile mixing to be responsible for the full relic abundance of dark matter additional new physics is needed beyond the keV-scale sterile neutrino itself. The extra ingredient we consider here is the presence of self-interactions among the sterile neutrinos. We examine whether active-to-sterile conversion is amplified enough in this scenario that the observed abundance of dark matter can be obtained with a subconstraint mixing angle. This turns out never to be the case in the region we explore: either self-interactions have too small an impact and cannot escape bounds on the mass and mixing angle, or they have too great an impact and cause dark matter to be overproduced. The sharp transition from marginal to excessive effectiveness occurs because a resonance criterion is met in the effective in-medium mixing angle. Once the system goes resonant the game is as good as over, as nonlinearity in the Boltzmann equation leads to runaway production of sterile neutrinos, beginning at a plasma temperature of a few hundred MeV and typically ending at a few tens of MeV. The scenario is therefore ruled out largely by its own dynamics. In this study we focus exclusively on mediators heavier than $\sim 1$ GeV; future work will extend the analysis to lighter mediators, allowing for contact to be made with the kinds of scenarios motivated by issues of small-scale structure.


title: Deep Neural Networks Improve Radiologists' Performance in Breast Cancer Screening

abstract: We present a deep convolutional neural network for breast cancer screening exam classification, trained and evaluated on over 200,000 exams (over 1,000,000 images). Our network achieves an AUC of 0.895 in predicting whether there is a cancer in the breast, when tested on the screening population. We attribute the high accuracy of our model to a two-stage training procedure, which allows us to use a very high-capacity patch-level network to learn from pixel-level labels alongside a network learning from macroscopic breast-level labels. To validate our model, we conducted a reader study with 14 readers, each reading 720 screening mammogram exams, and find our model to be as accurate as experienced radiologists when presented with the same data. Finally, we show that a hybrid model, averaging probability of malignancy predicted by a radiologist with a prediction of our neural network, is more accurate than either of the two separately. To better understand our results, we conduct a thorough analysis of our network's performance on different subpopulations of the screening population, model design, training procedure, errors, and properties of its internal representations.


title: Local Versus Global Distances for Zigzag Persistence Modules

abstract: This short note establishes explicit and broadly applicable relationships between persistence-based distances computed locally and globally. In particular, we show that the bottleneck distance between two zigzag persistence modules restricted to an interval is always bounded above by the distance between the unrestricted versions. While this result is not surprising, it could have different practical implications. We give two related applications for metric graph distances, as well as an extension for the matching distance between multiparameter persistence modules.


title: The X-ray variable sky as seen by MAXI: the future of dust echo tomography with bright Galactic X-ray bursts

abstract: Bright, short duration X-ray flares from accreting compact objects produce thin, dust scattering rings that enable dust echo tomography: high precision distance measurements and mapping of the line-of-sight distribution of dust. This work looks to the past activity of X-ray transient outbursts in order to predict the number of sight lines available for dust echo tomography. We search for and measure the properties of 3$\sigma$ significant flares in the 2-4 keV light curves of all objects available in the public MAXI archive. We derive a fluence sensitivity limit of $10^{-3}$ erg cm$^{-2}$ for the techniques used to analyze the light curves. This limits the study mainly to flares from Galactic X-ray sources. We obtain the number density of flares and estimate the total fluence of the corresponding dust echoes. However, the sharpness of a dust echo ring depends on the duration of a flare relative to quiescence. We select flares that are shorter than their corresponding quiescent period to calculate a number density distribution for dust echo rings as a function of fluence. The results are fit with a power law of slope $-2.3 \pm 0.1$. Extrapolating this to dimmer flares, we estimate that the next generation of X-ray telescopes will be 30 times more sensitive than current observatories, resulting in 10-30 dust ring echoes per year. The new telescopes will also be 10-100 times more sensitive than Chandra to dust ring echoes from the intergalactic medium.


title: Intermittent stickiness synchronization

abstract: This work uses the statistical properties of Finite-Time Lyapunov Exponents (FTLEs) to investigate the Intermittent Stickiness Synchronization (ISS) observed in the mixed phase space of high-dimensional Hamiltonian systems. Full Stickiness Synchronization (SS) occurs when all FTLEs from a chaotic trajectory tend to zero for arbitrary long time-windows. This behavior is a consequence of the sticky motion close to regular structures which live in the high-dimensional phase space and affects all unstable directions proportionally by the same amount, generating a kind of collective motion. Partial SS occurs when at least one FTLE approaches to zero. Thus, distinct degrees of partial SS may occur, depending on the values of nonlinearity and coupling parameters, on the dimension of the phase space and on the number of positive FTLEs. Through filtering procedures used to precisely characterize the sticky motion, we are able to compute the algebraic decay exponents of the ISS and to obtain remarkable evidence about the existence of a universal behavior related to the decay of time correlations encoded in such exponents. In addition, we show that even though the probability to find full SS is small compared to partial SSs, the full SS may appear for very long times due to the slow algebraic decay of time correlations in mixed phase space. In this sense, observations of very late intermittence between chaotic motion and full SSs become rare events.


title: Super-Resolution Fluorescence Imaging of Carbon Nanotubes Using a Nonlinear Excitonic Process

abstract: Highly efficient exciton-exciton annihilation process unique to one-dimensional systems is utilized for super-resolution imaging of air-suspended carbon nanotubes. Through the comparison of fluorescence signals in linear and sublinear regimes at different excitation powers, we extract the efficiency of the annihilation processes using conventional confocal microscopy. Spatial images of the annihilation rate of the excitons have resolution beyond the diffraction limit. We investigate excitation power dependence of the annihilation processes by experiment and Monte Carlo simulation, and the resolution improvement of the annihilation images can be quantitatively explained by the superlinearity of the annihilation process. We have also developed another method in which the cubic dependence of the annihilation rate on exciton density is utilized to achieve further sharpening of single nanotube images.


title: Standing waves of fixed period for $n+1$ vortex filaments

abstract: The $n+1$ vortex filament problem has explicit solutions consisting of $n$ parallel filaments of equal circulation in the form of nested polygons uniformly rotating around a central filament which has circulation of opposite sign. We show that when the relation between temporal and spatial periods is fixed at certain rational numbers, these configurations have an infinite number of homographic time dependent standing wave patterns that bifurcate from these uniformly rotating central configurations.


title: Experimental demonstration of switching entangled photons based on the Rydberg blockade effect

abstract: The long-range interaction between Rydberg-excited atoms endows a medium with large optical nonlinearity. Here, we demonstrate an optical switch to operate on a single photon from an entangled photon pair under a Rydberg electromagnetically induced transparency configuration. With the presence of the Rydberg blockade effect, we switch on a gate field to make the atomic medium nontransparent thereby absorbing the single photon emitted from another atomic ensemble via the spontaneous four-wave mixing process. In contrast to the case without a gate field, more than 50% of the photons sent to the switch are blocked, and finally achieve an effective single-photon switch. There are on average 1-2 gate photons per effective blockade sphere in one gate pulse. This switching effect on a single entangled photon depends on the principal quantum number and the photon number of the gate field. Our experimental progress is significant in the quantum information process especially in controlling the interaction between Rydberg atoms and entangled photon pairs.


title: Riemann--Hilbert Problems

abstract: These lectures introduce the method of nonlinear steepest descent for Riemann-Hilbert problems. This method finds use in studying asymptotics associated to a variety of special functions such as the Painlev\'{e} equations and orthogonal polynomials, in solving the inverse scattering problem for certain integrable systems, and in proving universality for certain classes of random matrix ensembles. These lectures highlight a few such applications.


title: On Extracting Data from Tables that are Encoded using HTML

abstract: Tables are a common means to display data in human-friendly formats. Many authors have worked on proposals to extract those data back since this has many interesting applications. In this article, we summarise and compare many of the proposals to extract data from tables that are encoded using HTML and have been published between $2000$ and $2018$. We first present a vocabulary that homogenises the terminology used in this field; next, we use it to summarise the proposals; finally, we compare them side by side. Our analysis highlights several challenges to which no proposal provides a conclusive solution and a few more that have not been addressed sufficiently; simply put, no proposal provides a complete solution to the problem, which seems to suggest that this research field shall keep active in the near future. We have also realised that there is no consensus regarding the datasets and the methods used to evaluate the proposals, which hampers comparing the experimental results.


title: Order-disorder transition in a two-dimensional associating lattice gas

abstract: We study an associating lattice gas (ALG) using Monte Carlo simulation and solutions on Husimi lattices. In this model, the molecules have an orientational degree of freedom and the interactions depend on the relative orientations of nearest-neighbor molecules. We focus on the transition between the high-density liquid (HDL) phase and the isotropic gas phase in the limit of full occupancy ($\mu \to \infty$). Simulation results show a continuous phase transition at $\tau_c=k_BT_c/\gamma=0.4763(1)$ (where $-\gamma$ is the bond energy) between the low-temperature HDL phase, with a non-vanishing mean orientation of the molecules, and the high-temperature isotropic phase. Results for critical exponents and the Binder cumulant indicate that the transition belongs to the three-state Potts model universality class, even though the ALG Hamiltonian does not have the full permutation symmetry of the Potts model. In contrast with simulation, the Husimi lattice results furnish a discontinuous phase transition, characterized by a discontinuity of the nematic order parameter. The transition temperatures ($\tau_c=0.51403$ and $0.51207$ for trees built with triangles and hexagons, respectively) are slightly higher than the one found via simulation. Since the Husimi lattice studies show that the ALG phase diagram features a discontinuous gas-HDL line for finite $\mu$, three possible scenarios arise. The first is that in the limit $\mu \to \infty$ the first-order line ends in a critical point; the second is a change in the nature of the transition at some finite chemical potential; the third is that the entire line is one of continuous phase transitions. Results from other ALG models and the fact that mean-field approximations show a discontinuous phase transition for the three-state Potts model (known to possess a continuous transition) lends some weight to the third alternative.


title: The Impact of Sex Education on Sexual Activity, Pregnancy, and Abortion

abstract: The purpose of this study is to find a relation between sex education and abortion in the United States. Accordingly, multivariate logistic regression is employed to study the relation between abortion and frequency of sex, pre-marriage sex, and pregnancy by rape. The finding shows the odds of abortion among those who have had premarital sex, more frequent sex before marriage, and been the victim of rape is higher than those who have not experienced any of these incidents. The output identified with one unit increase in pre-marriage sex the log-odds of abortion increases by 0.47. Similarly, it shows by one unit increase in the frequency of sex, the log-odds of abortion increases by 0.39. Also, for every additional pregnancy by rape, there is an expectation of a 3.17 increase in the log-odds of abortion. The findings of this study also suggests abortion is associated with sex education. Despite previous findings, this study shows the factors of age, having children, and social standing is not considered a burden to parents and thereby do not have a causal relation to abortion.


title: Publicising Lattice Field Theory through Visualisation

abstract: The gluon field configurations that form the foundation of every lattice QCD calculation contain a rich diversity of emergent nonperturbative structure. Visualisations of these phenomena not only serve to explain the concept of a nontrivial vacuum but also entertain a diverse audience from research funding panels to the next generation of science enthusiasts. In this brief review, a collection of QCD-vacuum visualisations is presented including the structure of chromo-electromagnetic fields, centre-cluster evolution at finite temperature, the structure of projected centre vortices, and novel correlations between the electromagnetic fields of QED and the chromo-electromagnetic fields of QCD in QED+QCD dynamical-fermion simulations from the QCDSF collaboration.


title: Prospection: Interpretable Plans From Language By Predicting the Future

abstract: High-level human instructions often correspond to behaviors with multiple implicit steps. In order for robots to be useful in the real world, they must be able to to reason over both motions and intermediate goals implied by human instructions. In this work, we propose a framework for learning representations that convert from a natural-language command to a sequence of intermediate goals for execution on a robot. A key feature of this framework is prospection, training an agent not just to correctly execute the prescribed command, but to predict a horizon of consequences of an action before taking it. We demonstrate the fidelity of plans generated by our framework when interpreting real, crowd-sourced natural language commands for a robot in simulated scenes.


title: Exact correlators in the Gaussian Hermitian matrix model

abstract: We present the $W_{1+\infty}$ constraints for the Gaussian Hermitian matrix model, where the constructed constraint operators yield the $W_{1+\infty}$ $n$-algebra. For the Virasoro constraints, we note that the constraint operators give the null 3-algebra. With the help of our Virasoro constraints, we derive a new effective formula for correlators in the Gaussian Hermitian matrix model.


title: Note on bounds for symmetric divergence measures

abstract: I. Sason obtained the tight bounds for symmetric divergence measures are derived by applying the results established by G. L. Gilardoni. In this article, we are going to report two kinds of extensions for the above results, namely classical q-extension and non-commutative extension.


title: Black holes in Gauss-Bonnet and Chern-Simons-scalar theory

abstract: We carry out the stability analysis of the Schwarzschild black hole in Gauss-Bonnet and Chern-Simons-scalar theory. Here, we introduce two quadratic scalar couplings ($\phi_1^2,\phi_2^2$) to Gauss-Bonnet and Chern-Simons terms, where the former term is parity-even, while the latter one is parity-odd. The perturbation equation for the scalar $\phi_1$ is the Klein-Gordon equation with an effective mass, while the perturbation equation for $\phi_2$ is coupled to the parity-odd metric perturbation, providing a system of two coupled equations. It turns out that the Schwarzschild black hole is unstable against $\phi_1$ perturbation, leading to scalarized black holes, while the black hole is stable against $\phi_2$ and metric perturbations, implying no scalarized black holes.


title: LookUP: Vision-Only Real-Time Precise Underground Localisation for Autonomous Mining Vehicles

abstract: A key capability for autonomous underground mining vehicles is real-time accurate localisation. While significant progress has been made, currently deployed systems have several limitations ranging from dependence on costly additional infrastructure to failure of both visual and range sensor-based techniques in highly aliased or visually challenging environments. In our previous work, we presented a lightweight coarse vision-based localisation system that could map and then localise to within a few metres in an underground mining environment. However, this level of precision is insufficient for providing a cheaper, more reliable vision-based automation alternative to current range sensor-based systems. Here we present a new precision localisation system dubbed "LookUP", which learns a neural-network-based pixel sampling strategy for estimating homographies based on ceiling-facing cameras without requiring any manual labelling. This new system runs in real time on limited computation resource and is demonstrated on two different underground mine sites, achieving real time performance at ~5 frames per second and a much improved average localisation error of ~1.2 metre.


title: Inequalities related to some types of entropies and divergences

abstract: The aim of this paper is to discuss new results concerning some kinds of parametric extended entropies and divergences. As a result of our studies for mathematical properties on entropy and divergence, we give new bounds for the Tsallis quasilinear entropy and divergence by applying the Hermite-Hadamard inequality. We also give bounds for biparametrical extended entropies and divergences which have been given in \cite{7}. In addition, we study $(r,q)$-quasilinear entropies and divergences as alternative biparametrical extended entropy and divergence, and then we give bounds for them. Finally we obtain inequalities for an extended Lin's divergence and some characterizations of Fermi-Dirac entropy and Bose-Einstein entropy.


title: High Landau levels of 2D electrons near the topological transition caused by interplay of spin-orbit and Zeeman energy shifts

abstract: In the presence of spin-orbit coupling two branches of the energy spectrum of 2D electrons get shifted in the momentum space. Application of in-plane magnetic field causes the splitting of the branches in energy. When both, spin-orbit coupling and Zeeman splitting are present, the branches of energy spectrum cross at certain energy. Near this energy, the Landau quantization becomes peculiar since semiclassical trajectories, corresponding to individual branches, get coupled. We study this coupling as a function of proximity to the topological transition. Remarkably, the dependence on the proximity is strongly asymmetric reflecting the specifics of the behavior of the trajectories near the crossing. Equally remarkable, on one side of the transition, the magnitude of coupling is an oscillating function of this proximity. These oscillations can be interpreted in terms of the St{\"u}ckelberg interference. Scaling of characteristic detuning with magnetic length is also unusual. This unusual behavior cannot be captured by simply linearizing the Fermi contours near the crossing point.


title: Global rough solution for $L^2$-critical semilinear heat equation in the negative Sobolev space

abstract: In this paper, we consider the Cauchy global problem for the $L^2$-critical semilinear heat equations $\partial_t h=\Delta h\pm |h|^{\frac4d}h, $ with $h(0,x)=h_0$, where $h$ is an unknown real function defined on $ \R^+\times\R^d$. In most of the studies on this subject, the initial data $h_0$ belongs to Lebesgue spaces $L^p(\R^d)$ for some $p\ge 2$ or to subcritical Sobolev space $H^{s}(\R^d)$ with $s>0$. We here prove that there exists some positive constant $\varepsilon_0$ depending on $d$, such that the Cauchy problem is locally and globally well-posed for any initial data $h_0$ which is radial, supported away from origin and in the negative Sobolev space $\dot H^{-\varepsilon_0}(\R^d)$ including $L^p(\R^d)$ with certain $p<2$ as subspace. Furthermore, unconditional uniqueness, and $L^2$-estimate both as time $t\to0$ and $t\to +\infty$ were considered.


title: Castelnuovo-Mumford regularity of representations of certain product categories

abstract: We show in this paper that representations of a finite product of categories satisfying certain combinatorial conditions have finite Castelnuovo-Mumford regularity if and only if they are presented in finite degrees, and hence the category consisting of them is abelian. These results apply to examples such as the categories $\mathrm{FI}^m$ and $\mathrm{FI}_G^m$.


title: An Exact Method for Constrained Maximization of the Conditional Value-at-Risk of a Class of Stochastic Submodular Functions

abstract: We consider a class of risk-averse submodular maximization problems (RASM) where the objective is the conditional value-at-risk (CVaR) of a random nondecreasing submodular function at a given risk level. We propose valid inequalities and an exact general method for solving RASM under the assumption that we have an efficient oracle that computes the CVaR of the random function. We demonstrate the proposed method on a stochastic set covering problem that admits an efficient CVaR oracle for the random coverage function.


title: Well-posedness for the Navier-Stokes equations in critical mixed-norm Lebesgue spaces

abstract: We study the Cauchy problem in $n$-dimensional space for the system of Navier-Stokes equations in critical mixed-norm Lebesgue spaces. Local well-posedness and global well-posedness of solutions are established in the class of critical mixed-norm Lebesgue spaces. Being in the mixed-norm Lebesgue spaces, both of the initial data and the solutions could be singular at certain points or decaying to zero at infinity with different rates in different spatial variable directions. Some of these singular rates could be very strong and some of the decaying rates could be significantly slow. Besides other interests, the results of the paper particularly show an interesting phenomena on the persistence of the anisotropic behavior of the initial data under the evolution. To achieve the goals, fundamental analysis theory such as Young's inequality, time decaying of solutions for heat equations, the boundedness of the Helmholtz-Leray projection, and the boundedness of the Riesz tranfroms are developed in mixed-norm Lebesgue spaces. These fundamental analysis results are independently topics of great interests and they are potentially useful in other problems.


title: First Results from the rapid-response spectrophotometric characterization of Near-Earth Objects

abstract: As part of our multi-observatory, multi-filter campaign, we present \rmi color observations of 82 Near-Earth Objects (NEOs) obtained with the RATIR instrument on the 1.5m robotic telescope at the San Pedro Martir's National Observatory in Mexico. Our project is particularly focused on rapid response observations of small ($\lesssim 850$ m) NEOs. The rapid response and the use of spectrophotometry allows us to constrain the taxonomic classification of NEOs with high efficiency. Here we present the methodology of our observations and our result, suggesting that the ratio of C-type to S-type asteroids in a size range of $\sim$30-850m is 1.1, which is in accordance with our previous results. We also find that 10$\%$ of all NEOs in our sample are neither C- nor S-type asteroids


title: Numerical simulation of liquid film formation and its heat transfer through vapor bubble expansion in a microchannel

abstract: The evaporation of vapor bubbles inside a microchannel is important to realize a device with high cooling performance. The liquid film formed on the solid surface is essential for evaporative heat transfer from solid to fluid; its formation process and heat transfer characteristics need to be investigated. The expansion process of a single vapor bubble via evaporative heat transfer in microchannels was evaluated via a numerical simulation in this study. In the calculation model, the working fluid used was saturated FC-72 at 0.1013 MPa and the channel diameter was 200 $\mu$m. The superheat of the initial temperature field and wall were considered as parameters. To evaluate the heat transfer characteristics, the time variation of liquid film thickness was evaluated. The averaged liquid film thickness had a correlation with the capillary number. Additionally, the dominant heat transfer mode was estimated by decomposing the heat transfer rate into the heat-transfer rate through the liquid film, rear edge, and wake. When the superheat was low, the heat transfer mostly occurred via liquid film evaporation; the heat flux through the liquid film could be predicted using the liquid film thickness. On the other hand, in cases of higher superheat, owing to rapid expansion of the vapor bubble, no evaporative heat transfer occurred through the liquid film around the bubble head. It could be inferred from this study that the relationship between the thickness of the thermal boundary layer of the bubble and liquid film thickness is important for predicting the cooling effect of this phenomena. When the vapor bubble grows in the high superheat liquid, the rapid growth makes the liquid film thick, and the thick liquid film prevents the heat transfer between the liquid-vapor interface and heated wall.


title: A Learning Framework for Distribution-Based Game-Theoretic Solution Concepts

abstract: The past few years have seen several works on learning economic solutions from data; these include optimal auction design, function optimization, stable payoffs in cooperative games and more. In this work, we provide a unified learning-theoretic methodology for modeling such problems, and establish tools for determining whether a given economic solution concept can be learned from data. Our learning theoretic framework generalizes a notion of function space dimension -- the graph dimension -- adapting it to the solution concept learning domain. We identify sufficient conditions for the PAC learnability of solution concepts, and show that results in existing works can be immediately derived using our methodology. Finally, we apply our methods in other economic domains, yielding a novel notion of PAC competitive equilibrium and PAC Condorcet winners.


title: Initial Mass Function Variation in two Elliptical Galaxies using Near-Infrared Tracers

abstract: Using integral field spectroscopy, we demonstrate that gravity-sensitive absorption features in the zJ-band (0.9--1.35 \micron) can constrain the low-mass stellar initial mass function (IMF) in the cores of two elliptical galaxies, M85 and M87. Compared to the visible bands, the near-infrared (NIR) is more sensitive to light from low-mass dwarf stars, whose relative importance is the primary subject of the debate over IMF variations in nearby galaxies. Our analysis compares the observed spectra to the latest stellar population synthesis models by employing two different methods: equivalent widths and spectral fitting. We find that the IMF slopes in M85 are similar to the canonical Milky Way IMF with a median IMF-mismatch parameter $\alpha_{K} = 1.26$. In contrast, we find that the IMF in M87 is steeper than a Salpeter IMF with $\alpha_{K} = 2.77$. The derived stellar population parameters, including the IMF slopes, are consistent with those from recent results in the visible bands based on spectroscopic and kinematic techniques. Certain elemental abundances, e.g. Na and Fe, have dramatic effects on the IMF-sensitive features and therefore the derived IMF slopes. We show evidence for a high [Na/H] $\sim$ 0.65 dex in the core of M85 from two independent \ion{Na}{1} absorption features. The high Na abundance may be the result of a recent galactic merger involving M85. This suggests that including [Na/H] in the stellar population model parameters is critical for constraining the IMF slopes in M85. These results confirm the viability of using NIR absorption features to investigate IMF variation in nearby galaxies.


title: The hydrogen molecule $\rm{H}_{2}$ in inclined configuration in a weak magnetic field

abstract: Highly accurate variational calculations, based on a few-parameter, physically adequate trial function, are carried out for the hydrogen molecule \hh in inclined configuration, where the molecular axis forms an angle $\theta$ with respect to the direction of a uniform constant magnetic field ${\bf B}$, for $B=0,\, 0.1,\, 0.175$ and $0.2\,$a.u. Three inclinations $\theta=0^\circ,\,45^\circ,\,90^\circ$ are studied in detail with emphasis to the ground state $1_g$. Diamagnetic and paramagnetic susceptibilities are calculated (for $\theta=45^\circ$ for the first time), they are in agreement with the experimental data and with other calculations. For $B=0,\, 0.1$ and $0.2\,$a.u. potential energy curves $E$ vs $R$ are built for each inclination, they are interpolated by simple, two-point Pad\'e approximant $Pade[2/6](R)$ with accuracy of not less than 4 significant digits. Spectra of rovibrational states are calculated for the first time. It was found that the optimal configuration of the ground state for $B \leq B_{cr}=0.178\,$a.u. corresponds always to the parallel configuration, $\theta=0$, thus, it is a $^1\Sigma_g$ state. The state $1_g$ remains bound for any magnetic field, becoming metastable for $B > B_{cr}$, while for $B_{cr} < B < 12$\,a.u. the ground state corresponds to two isolated hydrogen atoms with parallel spins.


title: The moving groups as the origin of the vertical phase space spiral

abstract: Using the Gaia data release 2 (DR2), we analyzed the distribution of stars in the close vicinity of the Sun in the full 3D position-velocity space. We have found no evidence of incomplete phase mixing in the vertical direction of the disk, which could be originated by some external events. We show that the vertical phase space spiral $Z$-$V_z$ is produced by the well-known moving groups (MGs), mainly by Coma-Berenices, Pleiades-Hyades and Sirius, when the statistical characteristics (mean, median, or mode) of the azimuthal velocity $V_\varphi$ are used to analyze the distribution in the vertical position-velocity plane. This result does not invoke external perturbations and is independent on the internal dynamical mechanisms that originate the MGs. Our conclusions counterbalance current arguments in favor of short-lived (between 300 and 900 Myr) structures in the solar neighborhood. Contrarily, they support the hypothesis of a longer formation time scale (around a few Gyr) for the MGs.


title: Numerically Stable Polynomially Coded Computing

abstract: We study the numerical stability of polynomial based encoding methods, which has emerged to be a powerful class of techniques for providing straggler and fault tolerance in the area of coded computing. Our contributions are as follows: 1) We construct new codes for matrix multiplication that achieve the same fault/straggler tolerance as the previously constructed MatDot Codes and Polynomial Codes. Unlike previous codes that use polynomials expanded in a monomial basis, our codes uses a basis of orthogonal polynomials. 2) We show that the condition number of every $m \times m$ sub-matrix of an $m \times n, n \geq m$ Chebyshev-Vandermonde matrix, evaluated on the $n$-point Chebyshev grid, grows as $O(n^{2(n-m)})$ for $n > m$. An implication of this result is that, when Chebyshev-Vandermonde matrices are used for coded computing, for a fixed number of redundant nodes $s=n-m,$ the condition number grows at most polynomially in the number of nodes $n$. 3) By specializing our orthogonal polynomial based constructions to Chebyshev polynomials, and using our condition number bound for Chebyshev-Vandermonde matrices, we construct new numerically stable techniques for coded matrix multiplication. We empirically demonstrate that our constructions have significantly lower numerical errors compared to previous approaches which involve inversion of Vandermonde matrices. We generalize our constructions to explore the trade-off between computation/communication and fault-tolerance. 4) We propose a numerically stable specialization of Lagrange coded computing. Motivated by our condition number bound, our approach involves the choice of evaluation points and a suitable decoding procedure that involves inversion of an appropriate Chebyshev-Vandermonde matrix. Our approach is demonstrated empirically to have lower numerical errors as compared to standard methods.


title: Walking with Confidence: Safety Regulation for Full Order Biped Models

abstract: Safety guarantees are valuable in the control of walking robots, as falling can be both dangerous and costly. Unfortunately, set-based tools for generating safety guarantees (such as sums-of-squares optimization) are typically restricted to simplified, low-dimensional models of walking robots. For more complex models, methods based on hybrid zero dynamics can ensure the local stability of a pre-specified limit cycle, but provide limited guarantees. This paper combines the benefits of both approaches by using sums-of-squares optimization on a hybrid zero dynamics manifold to generate a guaranteed safe set for a 10-dimensional walking robot model. Along with this set, this paper describes how to generate a controller that maintains safety by modifying the manifold parameters when on the edge of the safe set. The proposed approach, which is applied to a bipedal Rabbit model, provides a roadmap for applying sums-of-squares verification techniques to high dimensional systems. This opens the door for a broad set of tools that can generate safety guarantees and regulating controllers for complex walking robot models.


title: Traffic flow models with looking ahead-behind dynamics

abstract: Motivated by the traffic flow model with Arrhenius look-ahead relaxation dynamics introduced in [A. Sopasakis and M.A. Katsoulakis, SIAM J. Appl. Math., 66 (2006), p. 921--944], this paper proposes a traffic flow model with look ahead relaxation-behind intensification by inserting look behind intensification dynamics to the flux. Finite time shock formation conditions in the proposed model with various types of interaction potentials are identified. Several numerical experiments are performed in order to demonstrate the performance of the modified model. It is observed that, compare to other well-known traffic flow models, the model equipped with look ahead relaxation-behind intensification has both enhanced dispersive and smoothing effects


title: On Sampling Random Features From Empirical Leverage Scores: Implementation and Theoretical Guarantees

abstract: Random features provide a practical framework for large-scale kernel approximation and supervised learning. It has been shown that data-dependent sampling of random features using leverage scores can significantly reduce the number of features required to achieve optimal learning bounds. Leverage scores introduce an optimized distribution for features based on an infinite-dimensional integral operator (depending on input distribution), which is impractical to sample from. Focusing on empirical leverage scores in this paper, we establish an out-of-sample performance bound, revealing an interesting trade-off between the approximated kernel and the eigenvalue decay of another kernel in the domain of random features defined based on data distribution. Our experiments verify that the empirical algorithm consistently outperforms vanilla Monte Carlo sampling, and with a minor modification the method is even competitive to supervised data-dependent kernel learning, without using the output (label) information.


title: Generalised vector products in three-dimensional geometry

abstract: In three-dimensional Euclidean geometry, the scalar product produces a number associated to two vectors, while the vector product computes a vector perpendicular to them. These are key tools of physics, chemistry and engineering and supported by a rich vector calculus of 18th and 19th century results. This paper extends this calculus to arbitrary metrical geometries on three-dimensional space, generalising key results of Lagrange, Jacobi, Binet and Cauchy in a purely algebraic setting which applies also to general fields, including finite fields. We will then apply these vector theorems to set up the basic framework of rational trigonometry in the three-dimensional affine space and the related two-dimensional projective plane, and show an example of its applications to relativistic geometry.