title: Enhanced Sensitivity to Ultralight Bosonic Dark Matter in the Spectra of the Linear Radical SrOH

abstract: Coupling between Standard Model particles and theoretically well-motivated ultralight dark matter (UDM) candidates can lead to time variation of fundamental constants including the proton-to-electron mass ratio $\mu\equiv m_{p}/m_{e}\approx1836$. The presence of nearly-degenerate vibrational energy levels of different character in polyatomic molecules can result in significantly enhanced relative energy shifts in molecular spectra originating from $\partial_{t}\mu$, relaxing experimental complexity required for high-sensitivity measurements. We analyze the amplification of UDM effects in the spectrum of laser-cooled strontium monohydroxide (SrOH). SrOH is the first and, so far, the only polyatomic molecule to be directly laser cooled to sub-millikelvin temperatures, opening the possibility of long experimental coherence times. Because of the high enhancement factors ($\left|Q_{\mu}\right|\approx10^{3}$), measurements of the $\tilde{X}\left(200\right)\leftrightarrow\tilde{X}\left(03^{1}0\right)$ rovibrational transitions of SrOH in the microwave regime can result in $\sim10^{-17}$ fractional uncertainty in $\delta\mu/\mu$ with one day of integration. Furthermore, ultracold SrOH provides a promising platform for suppressing systematic errors. More complex SrOR radicals with additional vibrational modes arising from larger ligands R could lead to even greater enhancement factors.


title: On A Polytime Factorization Algorithm for Multilinear Polynomials Over F2

abstract: In 2010, A. Shpilka and I. Volkovich established a prominent result on the equivalence of polynomial factorization and identity testing. It follows from their result that a multilinear polynomial over the finite field of order 2 can be factored in time cubic in the size of the polynomial given as a string. Later, we have rediscovered this result and provided a simple factorization algorithm based on computations over derivatives of multilinear polynomials. The algorithm has been applied to solve problems of compact representation of various combinatorial structures, including Boolean functions and relational data tables. In this paper, we describe an improvement of this factorization algorithm and report on preliminary experimental analysis.


title: Finding forbidden minors in sublinear time: a $n^{1/2+o(1)}$-query one-sided tester for minor closed properties on bounded degree graphs

abstract: Let $G$ be an undirected, bounded degree graph with $n$ vertices. Fix a finite graph $H$, and suppose one must remove $\varepsilon n$ edges from $G$ to make it $H$-minor free (for some small constant $\varepsilon > 0$). We give an $n^{1/2+o(1)}$-time randomized procedure that, with high probability, finds an $H$-minor in such a graph. As an application, suppose one must remove $\varepsilon n$ edges from a bounded degree graph $G$ to make it planar. This result implies an algorithm, with the same running time, that produces a $K_{3,3}$ or $K_5$ minor in $G$. No prior sublinear time bound was known for this problem. By the graph minor theorem, we get an analogous result for any minor-closed property. Up to $n^{o(1)}$ factors, this resolves a conjecture of Benjamini-Schramm-Shapira (STOC 2008) on the existence of one-sided property testers for minor-closed properties. Furthermore, our algorithm is nearly optimal, by an $\Omega(\sqrt{n})$ lower bound of Czumaj et al (RSA 2014). Prior to this work, the only graphs $H$ for which non-trivial one-sided property testers were known for $H$-minor freeness are the following: $H$ being a forest or a cycle (Czumaj et al, RSA 2014), $K_{2,k}$, $(k\times 2)$-grid, and the $k$-circus (Fichtenberger et al, Arxiv 2017).


title: Diagonals of operators and Blaschke's enigma

abstract: We introduce new techniques allowing one to construct diagonals of bounded Hilbert space operators and operator tuples under "Blaschke-type" assumptions. This provides a new framework for a number of results in the literature and identifies, often large, subsets in the set of diagonals of arbitrary bounded operators (and their tuples). Moreover, our approach leads to substantial generalizations of the results due to Bourin, Herrero and Stout having assumptions of a similar nature.


title: Towards a viable scalar interpretation of $R_{D^{(*)}}$

abstract: Recent measurements of semileptonic B-meson decays seemingly imply violations of lepton flavor universality beyond the Standard Model predictions. With three-level explanations based on extended Higgs sectors being strongly challenged by the measurements of the $B_c^-$ lifetime, new theories invoking leptoquark or vector fields appear as the only feasible answer. However, in this work we show that simple scalar extensions of the Standard Model still offer a possible solution to the $B$ physics puzzle, owing to sizeable loop-level corrections which mimic the effects of new vector contributions. Considering a simplified model characterised by a charged and a neutral scalar particle, we verify the compatibility of the observed $R_{D^{(*)}}$ signal with the relevant collider bounds. We also study an embedding of the simplified model into a three-Higgs-doublet framework, and investigate its main phenomenological constraints.


title: Spiral bandwidth of four-wave mixing in Rb vapour

abstract: Laguerre-Gauss beams, and more generally the orbital angular momentum of light (OAM) provide valuable research tools for optical manipulation, processing, imaging and communication. Here we explore the high-efficiency frequency conversion of OAM in a four-wave mixing process in rubidium vapour. Conservation of the OAM in the two pump beams determines the total OAM shared by the generated light fields at 420 nm and 5.2 um - but not its distribution between them. We experimentally investigate the spiral bandwidth of the generated light modes as a function of pump OAM. A small pump OAM is transferred almost completely to the 420 nm beam. Increasing the total pump OAM broadens the OAM spectrum of the generated light, indicating OAM entanglement between the generated light fields. This clears the path to high-efficiency OAM entanglement between widely disparate wavelengths.


title: Hierarchically Structured Reinforcement Learning for Topically Coherent Visual Story Generation

abstract: We propose a hierarchically structured reinforcement learning approach to address the challenges of planning for generating coherent multi-sentence stories for the visual storytelling task. Within our framework, the task of generating a story given a sequence of images is divided across a two-level hierarchical decoder. The high-level decoder constructs a plan by generating a semantic concept (i.e., topic) for each image in sequence. The low-level decoder generates a sentence for each image using a semantic compositional network, which effectively grounds the sentence generation conditioned on the topic. The two decoders are jointly trained end-to-end using reinforcement learning. We evaluate our model on the visual storytelling (VIST) dataset. Empirical results from both automatic and human evaluations demonstrate that the proposed hierarchically structured reinforced training achieves significantly better performance compared to a strong flat deep reinforcement learning baseline.


title: Spitzer Space Telescope Infrared Observations of the Binary Neutron Star Merger GW170817

abstract: We present Spitzer Space Telescope 3.6 and 4.5 micron observations of the binary neutron star merger GW170817 at 43, 74, and 264 days post-merger. Using the final observation as a template, we uncover a source at the position of GW170817 at 4.5 micron with a brightness of 22.9+/-0.3 AB mag at 43 days and 23.8+/-0.3 AB mag at 74 days (the uncertainty is dominated by systematics from the image subtraction); no obvious source is detected at 3.6 micron to a 3-sigma limit of >23.3 AB mag in both epochs. The measured brightness is dimmer by a factor of about 2-3 times compared to our previously published kilonova model, which is based on UV, optical, and near-IR data at <30 days. However, the observed fading rate and color (m_{3.6}-m_{4.5}> 0 AB mag) are consistent with our model. We suggest that the discrepancy is likely due to a transition to the nebular phase, or a reduced thermalization efficiency at such late time. Using the Spitzer data as a guide, we briefly discuss the prospects of observing future binary neutron star mergers with Spitzer (in LIGO/Virgo Observing Run 3) and the James Webb Space Telescope (in LIGO/Virgo Observing Run 4 and beyond).


title: Masking: A New Perspective of Noisy Supervision

abstract: It is important to learn various types of classifiers given training data with noisy labels. Noisy labels, in the most popular noise model hitherto, are corrupted from ground-truth labels by an unknown noise transition matrix. Thus, by estimating this matrix, classifiers can escape from overfitting those noisy labels. However, such estimation is practically difficult, due to either the indirect nature of two-step approaches, or not big enough data to afford end-to-end approaches. In this paper, we propose a human-assisted approach called Masking that conveys human cognition of invalid class transitions and naturally speculates the structure of the noise transition matrix. To this end, we derive a structure-aware probabilistic model incorporating a structure prior, and solve the challenges from structure extraction and structure alignment. Thanks to Masking, we only estimate unmasked noise transition probabilities and the burden of estimation is tremendously reduced. We conduct extensive experiments on CIFAR-10 and CIFAR-100 with three noise structures as well as the industrial-level Clothing1M with agnostic noise structure, and the results show that Masking can improve the robustness of classifiers significantly.


title: KvN mechanics approach to the time-dependent frequency harmonic oscillator

abstract: Using the Ermakov-Lewis invariants appearing in KvN mechanics, the time-dependent frequency harmonic oscillator is studied. The analysis builds upon the operational dynamical model, from which it is possible to infer quantum or classical dynamics; thus, the mathematical structure governing the evolution will be the same in both cases. The Liouville operator associated with the time-dependent frequency harmonic oscillator can be transformed using an Ermakov-Lewis invariant, which is also time dependent and commutes with itself at any time. Finally, because the solution of the Ermakov equation is involved in the evolution of the classical state vector, we explore some analytical and numerical solutions.


title: Depth-Limited Solving for Imperfect-Information Games

abstract: A fundamental challenge in imperfect-information games is that states do not have well-defined values. As a result, depth-limited search algorithms used in single-agent settings and perfect-information games do not apply. This paper introduces a principled way to conduct depth-limited solving in imperfect-information games by allowing the opponent to choose among a number of strategies for the remainder of the game at the depth limit. Each one of these strategies results in a different set of values for leaf nodes. This forces an agent to be robust to the different strategies an opponent may employ. We demonstrate the effectiveness of this approach by building a master-level heads-up no-limit Texas hold'em poker AI that defeats two prior top agents using only a 4-core CPU and 16 GB of memory. Developing such a powerful agent would have previously required a supercomputer.


title: Learning Maximum-A-Posteriori Perturbation Models for Structured Prediction in Polynomial Time

abstract: MAP perturbation models have emerged as a powerful framework for inference in structured prediction. Such models provide a way to efficiently sample from the Gibbs distribution and facilitate predictions that are robust to random noise. In this paper, we propose a provably polynomial time randomized algorithm for learning the parameters of perturbed MAP predictors. Our approach is based on minimizing a novel Rademacher-based generalization bound on the expected loss of a perturbed MAP predictor, which can be computed in polynomial time. We obtain conditions under which our randomized learning algorithm can guarantee generalization to unseen examples.


title: Deformations of pairs of Kleinian singularities

abstract: Kleinian singularities, i.e., the varieties corresponding to the algebras of invariants of Kleinian groups are of fundamental importance for Algebraic geometry, Representation theory and Singularity theory. The filtered deformations of these algebras of invariants were classified by Slodowy (the commutative case) and Losev (the general case). To an inclusion of Kleinian groups, there is the corresponding inclusion of algebras of invariants. We classify deformations of these inclusions when a smaller subgroup is normal in the larger.


title: Thermo-magnetic evolution of the QCD strong coupling

abstract: We study the one-loop gluon polarization tensor at zero and finite temperature in the presence of a magnetic field, to extract the thermo-magnetic evolution of the strong coupling $\alpha_s$. We analyze four distinct regimes, to wit, the small and large field cases, both at zero and at high temperature. From a renormalization group analysis we show that at zero temperature, either for small or large magnetic fields, and for a fixed transferred momentum $Q^2$, $\alpha_s$ grows with the field strength with respect to its vacuum value. However, at high temperature and also for a fixed value of $Q^2$ we find two different cases: When the magnetic field is even larger than the squared temperature, $\alpha_s$ also grows with the field strength. On the contrary, when the squared temperature is larger than the magnetic field, a turnover behavior occurs and $\alpha_s$ decreases with the field strength. This thermo-magnetic behavior of $\alpha_s$ can help explain the inverse magnetic catalysis phenomenon found by lattice QCD calculations.


title: On the physical basis of biological signaling by interface pulses

abstract: Currently, biological signaling is envisaged as a combination of activation and movement, triggered by local molecular interactions and molecular diffusion, respectively. However, we here suggest, that other fundamental physical mechanisms might play an at least equally important role. We have recently shown that lipid interfaces permit the excitation and propagation of sound pulses. Here we demonstrate that these reversible perturbations can control the activity of membrane-embedded enzymes without a requirement for molecular transport. They can thus facilitate rapid communication between distant biological entities at the speed of sound, which is here of the order of 1 m/s within the membrane. The mechanism described provides a new physical framework for biological signaling that is fundamentally different from the molecular approach that currently dominates the textbooks.


title: Topological states of non-Hermitian systems

abstract: Recently, the search for topological states of matter has turned to non-Hermitian systems, which exhibit a rich variety of unique properties without Hermitian counterparts. Lattices modeled through non-Hermitian Hamiltonians appear in the context of photonic systems, where one needs to account for gain and loss, circuits of resonators, and also when modeling the lifetime due to interactions in condensed matter systems. Here we provide a brief overview of this rapidly growing subject, the search for topological states and a bulk-boundary correspondence in non-Hermitian systems.


title: The elastic and directed percolation backbone

abstract: We argue that the elastic backbone (EB) (union of shortest paths) on a cylindrical system, recently studied by Sampaio Filho et al. [Phys. Rev. Lett. 120, 175701 (2018)], is in fact the backbone of two-dimensional directed percolation (DP). We simulate the EB on the same system as considered by these authors, and also study the DP backbone directly using an algorithm that allows backbones to be generated in a completely periodic manner. We find that both the EB in the bulk and the DP backbone have a fractal dimension of $d_{b} = d_{B,\rm DP} = 1.681\,02(15)$ at the identical critical point $p_{c,\rm{DP}} \approx 0.705\,485\,22$. We also measure the fractal dimension at the edge of the EB system and for the full DP clusters, and find $d_e = d_{\rm DP} = 1.840 \, 54 (4)$. We argue that those two fractal dimensions follow from the DP exponents as $d_{B,\rm DP} = 2-2\beta/\nu_\parallel = 1.681 \, 07 2 (12) $ and $d_{\rm DP} = 2-\beta/\nu_\parallel = 1.840\, 536 (6)$. Our fractal dimensions differ from the value 1.750(3) found by Sampaio Filho et al., whose value may represent a crossover effect.


title: Concentration-Dependent Swelling and Structure of Ionic Microgels: Simulation and Theory of a Coarse-Grained Model

abstract: We study swelling and structural properties of ionic microgel suspensions within a comprehensive coarse-grained model that combines the polymeric and colloidal natures of microgels as permeable, compressible, charged spheres governed by effective interparticle interactions. The model synthesizes the Flory-Rehner theory of cross-linked polymer gels, the Hertz continuum theory of effective elastic interactions, and a theory of density-dependent effective electrostatic interactions. Implementing the model using Monte Carlo simulation and thermodynamic perturbation theory, we compute equilibrium particle size distributions, swelling ratios, volume fractions, net valences, radial distribution functions, and static structure factors as functions of concentration. Trial Monte Carlo moves comprising particle displacements and size variations are accepted or rejected based on the total change in elastic and electrostatic energies. The theory combines first-order thermodynamic perturbation and variational free energy approximations. For illustrative system parameters, theory and simulation agree closely at concentrations ranging from dilute to beyond particle overlap. With increasing concentration, as microgels deswell, we predict a decrease in the net valence and an unusual saturation of pair correlations. Comparison with experimental data for deionized, aqueous suspensions of PNIPAM particles demonstrates the capacity of the coarse-grained model to predict and interpret measured swelling behavior.


title: Anisotropic magnetocaloric effect in single crystals of CrI$_3$

abstract: We report a systematic investigation of dc magnetization and ac susceptibility, as well as anisotropic magnetocaloric effect in bulk CrI$_3$ single crystals. A second-stage magnetic transition was observed just below the Curie temperature $T_c$, indicating a two-step magnetic ordering. The low temperature thermal demagnetization could be well fitted by the spin-wave model rather than the single-particle model, confirming its localized magnetism. The maximum magnetic entropy change $-\Delta S_M^{max} \sim 5.65$ J kg$^{-1}$ K$^{-1}$ and the corresponding adiabatic temperature change $\Delta T_{ad} \sim 2.34$ K are achieved from heat capacity analysis with the magnetic field up to 9 T. Anisotropy of $\Delta S_M(T,H)$ was further investigated by isothermal magnetization, showing that the difference of $-\Delta S_M^{max}$ between the $ab$ plane and the $c$ axis reaches a maximum value $\sim$ 1.56 J kg$^{-1}$ K$^{-1}$ with the field change of 5 T. With the scaling analysis of $\Delta S_M$, the rescaled $\Delta S_M(T,H)$ curves collapse onto a universal curve, indicating a second-order type of the magnetic transition. Furthermore, the $-\Delta S_M^{max}$ follows the power law of $H^n$ with $n = 0.64(1)$, and the relative cooling power RCP depends on $H^m$ with $m = 1.12(1)$.


title: A theory on the absence of spurious solutions for nonconvex and nonsmooth optimization

abstract: We study the set of continuous functions that admit no spurious local optima (i.e. local minima that are not global minima) which we term \textit{global functions}. They satisfy various powerful properties for analyzing nonconvex and nonsmooth optimization problems. For instance, they satisfy a theorem akin to the fundamental uniform limit theorem in the analysis regarding continuous functions. Global functions are also endowed with useful properties regarding the composition of functions and change of variables. Using these new results, we show that a class of nonconvex and nonsmooth optimization problems arising in tensor decomposition applications are global functions. This is the first result concerning nonconvex methods for nonsmooth objective functions. Our result provides a theoretical guarantee for the widely-used $\ell_1$ norm to avoid outliers in nonconvex optimization.


title: STREAMFINDER II: A possible fanning structure parallel to the GD-1 stream in Pan-STARRS1

abstract: STREAMFINDER is a new algorithm that we have built to detect stellar streams in an automated and systematic way in astrophysical datasets that possess any combination of positional and kinematic information. In Paper I, we introduced the methodology and the workings of our algorithm and showed that it is capable of detecting ultra-faint and distant halo stream structures containing as few as $\sim 15$ members ($\Sigma_{\rm G} \sim 33.6\, {\rm mag \, arcsec^{-2}}$) in the Gaia dataset. Here, we test the method with real proper motion data from the Pan-STARRS1 survey, and by selecting targets down to $r_{0}=18.5$ mag we show that it is able to detect the GD-1 stellar stream, whereas the structure remains below a useful detection limit when using a Matched Filter technique. The radial velocity solutions provided by STREAMFINDER for GD-1 candidate members are found to be in good agreement with observations. Furthermore, our algorithm detects a $\sim 40^{\deg}$ long structure approximately parallel to GD-1, and which fans out from it, possibly a sign of stream-fanning due to the triaxiality of the Galactic potential. This analysis shows the promise of this method for detecting and analysing stellar streams in the upcoming Gaia DR2 catalogue.


title: Bias-Reduced Uncertainty Estimation for Deep Neural Classifiers

abstract: We consider the problem of uncertainty estimation in the context of (non-Bayesian) deep neural classification. In this context, all known methods are based on extracting uncertainty signals from a trained network optimized to solve the classification problem at hand. We demonstrate that such techniques tend to introduce biased estimates for instances whose predictions are supposed to be highly confident. We argue that this deficiency is an artifact of the dynamics of training with SGD-like optimizers, and it has some properties similar to overfitting. Based on this observation, we develop an uncertainty estimation algorithm that selectively estimates the uncertainty of highly confident points, using earlier snapshots of the trained model, before their estimates are jittered (and way before they are ready for actual classification). We present extensive experiments indicating that the proposed algorithm provides uncertainty estimates that are consistently better than all known methods.


title: Gradient descent in hyperbolic space

abstract: Gradient descent generalises naturally to Riemannian manifolds, and to hyperbolic $n$-space, in particular. Namely, having calculated the gradient at the point on the manifold representing the model parameters, the updated point is obtained by travelling along the geodesic passing in the direction of the gradient. Some recent works employing optimisation in hyperbolic space have not attempted this procedure, however, employing instead various approximations to avoid a calculation that was considered to be too complicated. In this tutorial, we demonstrate that in the hyperboloid model of hyperbolic space, the necessary calculations to perform gradient descent are in fact straight-forward. The advantages of the approach are then both illustrated and quantified for the optimisation problem of computing the Fr\'echet mean (i.e. barycentre) of points in hyperbolic space.


title: Measuring Congressional District Meandering

abstract: In recent decades, state legislatures have often drawn U.S. Congressional voting districts that look---to the human eye---to be rather twisted. In this paper, we propose a method to measure how much districts "meander" via a computation of the medial axis of the region. We then compare this to the medial axis of the convex hull of the district to obtain the {\em medial-hull ratio}: a dimensionless quantity that captures the district's irregularity. We compute this quantity for many example Congressional districts.


title: Possible Worlds, Incompleteness and Undefinability

abstract: This short squib looks at how using a broader definition of G\"odel numbering to mimic the accessibility relation between possible worlds results in two-world systems that sidestep undecidable sentences as well as the Liar paradox.


title: Spontaneous and Stimulated Emissions of Quantum Free-Electron Wavepackets - QED Analysis

abstract: Do the wavepacket-size of free-electron wavefunction and its history have physical effect in its interaction with light? Here we answer this problem by analyzing a QED model, considering both spontaneous and stimulated emission of quantized radiation field. For coherent radiation (Glauber state), we confirm that stimulated emission/absorption of photons has a dependence on wavepacket size that decays when it exceeds the interacting radiation wavelength, consistently and complementarily with Schrodinger equation analysis of wavepacket acceleration in classical electromagnetic field. Furthermore, the stimulated emission of modulated electron wavepacket with coherently-bunched profiles has characteristic harmonic emission spectrum that is also wavepacket size dependent but beyond the frequency cut-off. In either case, there is no wavepacket dependent emission of Fock state radiation, and particularly the vacuum state spontaneous emission is wavepacket-independent. The transition of radiation emission from the classical point-particle limit to the quantum electron wavefunction limit is demonstrated in electron wavepacket representation. It indicates a way for measuring the wavepacket size of single electron wavefunction, and suggests a new direction for exploring light-matter interaction fundamentally.

categories: quant-ph physics.acc-ph physics.optics physics.plasm-ph


title: The Exoplanet Population Observation Simulator. I - The Inner Edges of Planetary Systems

abstract: The Kepler survey provides a statistical census of planetary systems out to the habitable zone. Because most planets are non-transiting, orbital architectures are best estimated using simulated observations of ensemble populations. Here, we introduce EPOS, the Exoplanet Population Observation Simulator, to estimate the prevalence and orbital architectures of multi-planet systems based on the latest Kepler data release, DR25. We estimate that at least 42% of sun-like stars have nearly coplanar planetary systems with 7 or more exoplanets. The fraction of stars with at least one planet within 1 au could be as high as 100% depending on assumptions about the distribution of single transiting planets. We estimate an occurrence rate of planets in the habitable zone around sun-like stars of eta_earth=36+-14%. The innermost planets in multi-planet systems are clustered around an orbital period of 10 days (0.1 au), reminiscent of the protoplanetary disk inner edge or could be explained by a planet trap at that location. Only a small fraction of planetary systems have the innermost planet at long orbital periods, with fewer than ~8% and ~3% having no planet interior to the orbit of Mercury and Venus, respectively. These results reinforce the view that the solar system is not a typical planetary system, but an outlier among the distribution of known exoplanetary systems. We predict that at least half of the habitable zone exoplanets are accompanied by (non-transiting) planets at shorter orbital periods, hence knowledge of a close-in exoplanet could be used as a way to optimize the search for Earth-size planets in the Habitable Zone with future direct imaging missions.


title: A Triple Origin for the Heavy and Low-Spin Binary Black Holes Detected by LIGO/Virgo

abstract: We explore the masses, merger rates, eccentricities, and spins for field binary black holes driven to merger by a third companion through the Lidov-Kozai mechanism. Using a population synthesis approach, we model the creation of stellar-mass black hole triples across a range of different initial conditions and stellar metallicities. We find that the production of triple-mediated mergers is enhanced at low metallicities by a factor of ~100 due to the lower black hole natal kicks and reduced stellar mass loss. These triples naturally yield heavy binary black holes with near-zero effective spins, consistent with most of the mergers observed to date. This process produces a merger rate of between 2 and 25 Gpc^-3 yr^-1 in the local universe, suggesting that the Lidov-Kozai mechanism can potentially explain all of the low-spin, heavy black hole mergers observed by Advanced LIGO/Virgo. Finally, we show that triples admit a unique eccentricity and spin distribution that will allow this model to be tested in the near future.


title: Electric control of the heat flux through electrophononic effects

abstract: We demonstrate a fully electric control of the heat flux, which can be continuously modulated by an externally applied electric field in PbTiO$_3$, a prototypical ferroelectric perovskite, revealing the mechanisms by which experimentally accessible fields can be used to tune the thermal conductivity by as much as 50% at room temperature.


title: Stellar Mass Profiles of Quiescent Galaxies in Different Environments at $z\sim0$

abstract: We present the stellar mass profiles of 147 isolated quiescent galaxies in very low-density environments (i.e., void regions) in the local Universe ($0.01<z<0.06$) from the Sloan Digital Sky Survey. These galaxies have stellar masses between $ 9.8\lesssim \log(M_{\ast}/M_{\odot})\lesssim11.2$ and they represent $\sim15\%$ of the whole galaxy population in the void regions down to $M_{r} = -19$. We do not find any isolated quiescent galaxies with $\log(M_{\ast}/M_{\odot})\gtrsim11.2$. We compare the stellar mass profiles of these isolated quiescent galaxies with the profiles of stellar mass-matched samples of the quiescent galaxies in group and cluster environments. We find that, at fixed mass, quiescent galaxies in voids have similar central ($1$ kpc) mass densities ($\Sigma_1$) and central velocity dispersions ($\sigma_1$) compared to their counterparts in groups and clusters. We show that quiescent galaxies in voids have at most $10-25\%$ smaller half-mass (and half-light) sizes compared to quiescent galaxies in groups and clusters. We conclude that for the intermediate stellar mass range of $10^{10}-10^{11}M_{\odot}$ in the local Universe, environmental mechanisms have no significant additional effect on the mass profiles of the quiescent galaxies.


title: Fast scrambling on sparse graphs

abstract: Given a quantum many-body system with few-body interactions, how rapidly can quantum information be hidden during time evolution? The fast scrambling conjecture is that the time to thoroughly mix information among N degrees of freedom grows at least logarithmically in N. We derive this inequality for generic quantum systems at infinite temperature, bounding the scrambling time by a finite decay time of local quantum correlations at late times. Using Lieb-Robinson bounds, generalized Sachdev-Ye-Kitaev models, and random unitary circuits, we propose that a logarithmic scrambling time can be achieved in most quantum systems with sparse connectivity. These models also elucidate how quantum chaos is not universally related to scrambling: we construct random few-body circuits with infinite Lyapunov exponent but logarithmic scrambling time. We discuss analogies between quantum models on graphs and quantum black holes, and suggest methods to experimentally study scrambling with as many as 100 sparsely-connected quantum degrees of freedom.

categories: cond-mat.str-el cond-mat.stat-mech hep-th quant-ph


title: Spin-orbit-driven electron pairing in two dimensions

abstract: We show that the spin-orbit interaction (SOI) arising due to the in-plane electric field of the Coulomb repulsion between electrons in a two-dimensional quantum well produces an attractive component in the pair interaction Hamiltonian that depends on the spins and momenta of electrons. If the Rashba SOI constant of the material is high enough the attractive component overcomes the Coulomb repulsion and the centrifugal barrier, which leads to the formation of the two-electron bound states. There are two distinct types of two-electron bound states. The relative bound states are formed by the electrons orbiting around their common barycenter. They have the triplet spin structure and are independent of the center-of-mass momentum. In contrast, the convective bound states are formed because of the center-of-mass motion, which couples the electrons with opposite spins. The binding energy in the meV range is attainable for realistic conditions.


title: Emission from the Ionized Gaseous Halos of Low Redshift Galaxies and Their Neighbors

abstract: Using a sample of nearly half a million galaxies, intersected by over 8 million lines of sight from the Sloan Digital Sky Survey Data Release 12, we extend our previous study of the recombination radiation emitted by the gaseous halos of nearby galaxies. We identify an inflection in the radial profile of the H$\alpha$+N[{\small II}] radial emission profile at a projected radius of $\sim 50$ kpc and suggest that beyond this radius the emission from ionized gas in spatially correlated halos dominates the profile. We confirm that this is a viable hypothesis using results from a highly simplified theoretical treatment in which the dark matter halo distribution from cosmological simulations is straightforwardly populated with gas. Whether we fit the fraction of halo gas in a cooler (T $= 12,000$ K), smooth ($c = 1$) component (0.26 for galaxies with M$_* = 10^{10.88}$ M$_\odot$ and 0.34 for those with M$_* = 10^{10.18}$ M$_\odot$) or take independent values of this fraction from published hydrodynamical simulations (0.19 and 0.38, respectively), this model successfully reproduces the radial location and amplitude of the observed inflection. We also observe that the physical nature of the gaseous halo connects to primary galaxy morphology beyond any relationship to the galaxy's stellar mass and star formation rate. We explore whether the model reproduces behavior related to the central galaxy's stellar mass, star formation rate, and morphology. We find that it is unsuccessful in reproducing the observations at this level of detail and discuss various shortcomings of our simple model that may be responsible.


title: Cosmological constraints with self-interacting sterile neutrinos

abstract: In this work we revisit the question of whether Cosmology can be made compatible with scenarios with light sterile neutrinos, as invoked to explain the SBL anomalies, in the presence of self-interaction among sterile neutrinos mediated by massive gauge bosons. We examine this proposal by deriving the cosmological predictions of the model in a wide range of the model parameters including the effective interaction strength $G_X$, sterile neutrino mass $m_{\rm st}$ and active-sterile mixings. With those we perform a statistical analysis of the cosmological data from BBN, CMB, and BAO data to infer the posterior probabilities of the sterile self-interaction model parameters. BBN mostly provides information about the effective interaction strength and we find that $\log_{10}(G_X/G_F)\geq 3.6$ can describe the primordial abundances at 95% CL. Our analysis of CMB and BAO data show that when allowing a wide prior for the sterile neutrino mass its posterior is bounded to $m_{\rm st}\leq 0.91~{\rm eV}$ (95% CL) considering CMB data only and $m_{\rm st}\leq 0.37~{\rm eV}$ (95% CL) when adding the BAO information. So the mass bounds are slighly relaxed compared with that of a non-interacting sterile neutrino model but a sterile neutrino mass of 1 eV is still excluded at more than $2\sigma$ CL. Conversely if fixing the sterile neutrino mass and mixing to the values prefered by short baseline data we find that while that CMB data alone favors the self-interacting scenario, including the BAO information severly degrades the agreement with the model. Altogether we conclude then that adding the self-interaction can alleviate the tension between eV sterile neutrinos and CMB data, but when including also the BAO results the self-interacting sterile neutrino model cannot lead to a satisfactory description of the data.


title: Sign-Problem-Free Fermionic Quantum Monte Carlo: Developments and Applications

abstract: Reliable simulations of correlated quantum systems, including high-temperature superconductors and frustrated magnets, are increasingly desired nowadays to further understanding of essential features in such systems. Quantum Monte Carlo (QMC) is a unique numerically-exact and intrinsically-unbiased method to simulate interacting quantum many-body systems. More importantly, when QMC simulations are free from the notorious fermion-sign problem, they can reliably simulate interacting quantum models with large system size and low temperature to reveal low-energy physics such as spontaneously-broken symmetries and universal quantum critical behaviors. Here, we concisely review recent progresses made in developing new sign-problem-free QMC algorithms, including those employing Majorana representation and those utilizing hot-spot physics. We also discuss applications of these novel sign-problem-free QMC algorithms in simulations of various interesting quantum many-body models. Finally, we discuss possible future directions of designing sign-problem-free QMC methods.

categories: cond-mat.str-el cond-mat.quant-gas cond-mat.stat-mech cond-mat.supr-con


title: Searching for All-Scale Anisotropies in the Arrival Directions of Cosmic Rays above the Ankle

abstract: The Pierre Auger Observatory has recently reported the detection of a dipole anisotropy in the arrival directions of cosmic rays above 8 EeV with a post-trial significance of more than 5.2$\sigma$. This observation has profound consequences for the distribution and composition of candidate sources of cosmic rays above the ankle (3-5 EeV). In this paper we search for the presence of anisotropies on all angular scales in public Auger data. The analysis follows a likelihood-based reconstruction method that automatically accounts for variations in the observatory's angular acceptance and background rate. Our best-fit dipole anisotropy in the equatorial plane has an amplitude of 5.3 $\pm$ 1.3 percent and right ascension angle of 103 $\pm$ 15 degrees, consistent with the results of the Pierre Auger Collaboration. We do not find evidence for the presence of medium- or small-scale anisotropies. The method outlined in this paper is well-suited for the future analyses of cosmic ray anisotropies below the ankle, where cosmic ray detection in surface arrays is not fully efficient and dominated by systematic uncertainties.


title: A catalogue of faint local radio AGN and the properties of their host galaxies

abstract: We present a catalogue of local ($z<0.1$) galaxies that contain faint AGN. We select these objects by identifying galaxies that exhibit a significant excess in their radio luminosities, compared to what is expected from the observed levels of star-formation activity in these systems. This is achieved by comparing the optical (spectroscopic) star formation rate (SFR) to the 1.4~GHz luminosity measured from the FIRST survey. The majority of the AGN identified in this study are fainter than those in previous work, such as in the Best and Heckman (2012) catalogue. We show that these faint AGN make a non-negligible contribution to the radio luminosity function at low luminosities (below $10^{22.5}$ W Hz$^{-1}$), and host $\sim$13 per cent of the local radio luminosity budget. Their host galaxies are predominantly high stellar-mass systems (with a median stellar mass of $10^{11}\textrm{M}_{\odot}$), are found across a range of environments (but typically in denser environments than star-forming galaxies) and have early-type morphologies. This study demonstrates a general technique to identify AGN in galaxy populations where reliable optical SFRs can be extracted using spectro-photometry and where radio data are also available so that a radio excess can be measured. Our results also demonstrate that it is unsafe to infer SFRs from radio emission \textit{alone}, even if bright AGN have been excluded from a sample, since there is a significant population of faint radio AGN which may contaminate the radio-derived SFRs.


title: $b \to c \tau \nu_{\tau}$ Decays: A Catalogue to Compare, Constrain, and Correlate New Physics Effects

abstract: In this article, we have predicted the standard model (SM) values of the asymmetric and angular observables in $B\to D^{(\ast)}\tau\nu_{\tau}$ decays, using the results of the new up-to-date analysis in $B\to D^{(*)}\ell\nu_{\ell}$. We have also revisited the SM prediction of the inclusive ratio $\mathcal{R}_{X_c}$, and have given its values in different schemes of the charm quark mass. This is the first analysis which includes all the known corrections in the SM. In addition, we have analysed the $b\to c\tau\nu_\tau$ decay modes in a model-independent framework of effective field theory beyond the standard model. Considering all the possible combinations of the effective operators in $b \to c \tau\nu_{\tau}$ decays and using the Akaike Information Criterion, we find out the scenarios which can best explain the available data on these channels. In the selected scenarios, best-fit values and correlations of the new parameters are extracted. Using these results, predictions are made on various observables in the exclusive and inclusive semitaunic $b \to c $ decays. The graphical correlations between these observables are shown, which are found to be useful in discriminating various new physics scenarios.


title: Spin Response and Collective Modes in Simple Metal Dichalcogenides

abstract: Transition metal dichalcogenide (TMD) monolayers are interesting materials in part because of their strong spin-orbit coupling. This leads to intrinsic spin-splitting of opposite signs in opposite valleys, so the valleys are intrinsically spin-polarized when hole-doped. We study spin response in a simple model of these materials, with an eye to identifying sharp collective modes (i.e, spin-waves) that are more commonly characteristic of ferromagnets. We demonstrate that such modes exist for arbitrarily weak repulsive interactions, even when they are too weak to induce spontaneous ferromagnetism. The behavior of the spin response is explored for a range of hole dopings and interaction strengths.


title: "Direct" Gas-Phase Metallicity in Local Analogs of High-Redshift Galaxies: Empirical Metallicity Calibrations for High-Redshift Star-Forming Galaxies

abstract: We study the direct gas-phase oxygen abundance using the well-detected auroral [OIII]4363 line in the stacked spectra of a sample of local analogs of high-redshift galaxies. These local analogs share the same location as $z\sim2$ star-forming galaxies on the Baldwin-Phillips-Terlevich diagram. This type of analog has the same ionized interstellar medium (ISM) properties as high-redshift galaxies. We establish empirical metallicity calibrations between the direct gas-phase oxygen abundances ($7.8<12+\log(\rm{O/H})<8.4$) and the N2/O3N2 indices in our local analogs. We find significant systematic offsets between the metallicity calibrations for our local analogs of high-redshift galaxies and those derived from the local HII regions and a sample of local reference galaxies selected from the Sloan Digital Sky Survey (SDSS). The N2 and O3N2 metallicities will be underestimated by 0.05-0.1~dex relative to our calibration if one simply applies the local metallicity calibration in previous studies to high-redshift galaxies. Local metallicity calibrations also cause discrepancies of metallicity measurements in high-redshift galaxies using the N2 and O3N2 indicators. In contrast, our new calibrations produce consistent metallicities between these two indicators. We also derive metallicity calibrations for R23, O32, log([OIII]/H\beta), and log([NeIII]/[OII]) indices in our local analogs, which show significant offset compared to those in the SDSS reference galaxies. By comparing with MAPPINGS photoionization models, the different empirical metallicity calibration relations in the local analogs and the SDSS reference galaxies can be shown to be primarily due to the change of ionized ISM conditions.


title: Supercurrent anomaly and gauge invariance in N=1 supersymmetric Yang-Mills theory

abstract: We analyse Feynman diagram calculational issues related to the quantum breaking of supercurrent conservation in a supersymmetric non-abelian Yang-Mills theory. For the sake of simplicity, we take a zero mass gauge field multiplet interacting with a massless Majorana spin-$1/2$ field in the adjoint representation of $SU(2)$. We shed light on a long-standing controversy regarding the perturbative evaluation of the supercurrent anomaly in connection with gauge and superconformal symmetry in different frameworks. We find that only superconformal symmetry is unambiguously broken using an invariant four dimensional regularization and compare with the triangle AVV anomaly. Subtleties related to momentum routing invariance in the loops of diagrams and Clifford algebra evaluation inside divergent integrals are also discussed in connection with finite and undetermined quantities in Feynman amplitudes.


title: Nonreciprocal superposition state in antiferromagnetic optospintronics

abstract: The absence of net magnetization, which forbids any stray magnetic fields, is one of the greatest advantages of antiferromagnets in device applications. In conventional antiferromagnets, however, spin current cannot be extracted without the aid of a static magnetic field. Here, we develop a theory of antiferromagnetic opto-spintronics to resolve this fundamental dilemma. By coupling a linearly polarized photon and nonreciprocal magnon bands, we construct a superposition state of left- and right-handed magnon states with opposite group velocities. We numerically demonstrate that by using this superposition state, an antiferromagnetic spin current can be efficiently generated without a net magnetic field including net magnetization. We also find that the breakdown of the superposition state induces the stripe superfluid phase of a two-component Bose-Einstein condensate. Our results lay the foundation for manipulating the superposition states of emergent particles in devices.


title: Performance of quantum error correction with coherent errors

abstract: We compare the performance of quantum error correcting codes when memory errors are unitary with the more familiar case of dephasing noise. For a wide range of codes we analytically compute the effective logical channel that results when the error correction steps are performed noiselessly. Our examples include the entire family of repetition codes, the 5-qubit, Steane, Shor, and surface codes. When errors are measured in terms of the diamond norm, we find that the error correction is typically much more effective for unitary errors than for dephasing. We observe this behavior for a wide range of codes after a single level of encoding, and in the thresholds of concatenated codes using hard decoders. We show that this holds with great generality by proving a bound on the performance of any stabilizer code when the noise at the physical level is unitary. By comparing the diamond norm error $D'_\diamond$ of the logical qubit with the same quantity at the physical level $D_\diamond$, we show that $D'_\diamond \le c D^d_\diamond $ where $d$ is the distance of the code and $c$ is constant that depends on the code but not on the error. This bound compares very favorably to the performance of error correction for dephasing noise and other Pauli channels, where an error correcting code of odd distance $d$ will exhibit a scaling $D'_\diamond \sim D_\diamond^{(d+1)/2}$.


title: Trans-Neptunian Objects Transiently Stuck in Neptune's Mean Motion Resonances: Numerical simulations of the current population

abstract: A substantial fraction of our solar system's trans-Neptunian objects (TNOs) are in mean motion resonance with Neptune. Many of these objects were likely caught into resonances by planetary migration---either smooth or stochastic---approximately 4 Gyr ago. Some, however, gravitationally scattered off of Neptune and became transiently stuck in more recent events. Here, we use numerical simulations to predict the number of transiently-stuck objects, captured from the current actively scattering population, that occupy 111 resonances at semimajor axes $a=$30--100 au. Our source population is an observationally constrained model of the currently-scattering TNOs. We predict that, integrated across all resonances at these distances, the current transient sticking population comprises 40\% of total transiently-stuck+scattering TNOs, suggesting that these objects should be treated as a single population. We compute the relative distribution of transiently-stuck objects across all $p$:$q$ resonances with $1/6 \le q/p < 1$, $p<40$, and $q<20$, providing predictions for the population of transient objects with $H_r < 8.66$ in each resonance. We find that the relative populations are approximately proportional to each resonance's libration period and confirm that the importance of transient sticking increases with semimajor axis in the studied range. We calculate the expected distribution of libration amplitudes for stuck objects and demonstrate that observational constraints indicate that both the total number and the amplitude-distribution of 5:2 resonant TNOs are inconsistent with a population dominated by transient sticking from the current scattering disk. The 5:2 resonance hence poses a challenge for leading theories of Kuiper belt sculpting.


title: The stochastic gravitational-wave background in the absence of horizons

abstract: Gravitational-wave astronomy has the potential to explore one of the deepest and most puzzling aspects of Einstein's theory: the existence of black holes. A plethora of ultracompact, horizonless objects have been proposed to arise in models inspired by quantum gravity. These objects may solve Hawking's information-loss paradox and the singularity problem associated with black holes, while mimicking almost all of their classical properties. They are, however, generically unstable on relatively short timescales. Here, we show that this "ergoregion instability" leads to a strong stochastic background of gravitational waves, at a level detectable by current and future gravitational-wave detectors. The absence of such background in the first observation run of Advanced LIGO already imposes the most stringent limits to date on black-hole alternatives, showing that certain models of "quantum-dressed" stellar black holes can be at most a small percentage of the total population. The future LISA mission will allow for similar constraints on supermassive black-hole mimickers.


title: On the consistency of the expansion with the perturbations

abstract: Assuming a simple form for the growth index gamma(z) depending on two parameters gamma_0 = gamma(z=0) and gamma_1 = gamma'(z=0), we show that these parameters can be constrained using background expansion data. We explore systematically the preferred region in this parameter space. Inside General Relativity we obtain that models with a quasi-static growth index and gamma_1 = -0.02 are favoured. We find further the lower bounds gamma_0 > 0.53 and gamma_1 > -0.15 for models inside GR. Models outside GR having the same background expansion as LCDM and arbitrary gamma(z) with gamma_0 = gamma_0^{LCDM}, satisfy G_{eff,0}>G for gamma_1 > gamma_1^{LCDM}, and G_{eff,0}<G for gamma_1 < gamma_1^{LCDM}. The first models will cross downwards the value G_{eff}=G on very low redshifts z<0.3, while the second models will cross upwards G_{eff}=G in the same redshift range. This makes the realization of such modified gravity models even more problematic.


title: Collider Bounds on 2-Higgs Doublet Models with $U(1)_X$ Gauge Symmetries

abstract: 2-Higgs Doublet Models (2HDMs) typically need to invoke an ad-hoc discrete symmetry to avoid severe flavor bounds and in addition feature massless neutrinos, thus falling short of naturally complying with existing data. However, when augmented by an Abelian gauge symmetry naturally incorporating neutrino masses via a type-I seesaw mechanism while at the same time escaping flavor changing interactions, such enlarged 2HDMs become very attractive phenomenologically. In such frameworks, the distinctive element is the $Z'$ gauge boson generated by the spontaneous breaking of the Abelian group $U(1)_X$. In this work, we derive updated collider bounds on it. Several theoretical setups are possible, each with different and sometimes suppressed couplings to quarks and leptons. Thus, complementary data from dijet and dilepton resonance searches need to be considered to fully probe these objects. We employ the corresponding datasets as obtained at the Large Hadron Collider (LHC) at the 13 TeV CMs energy for $\mathcal{L}=12,36$ and $300$ fb$^{-1}$ of luminosity. Moreover, we present the potential sensitivity to such $Z'$s of the High Luminosity LHC (HL-LHC) and High Energy LHC (HE-LHC).


title: Nearly Flat Chern Bands in Moir\'e Superlattices

abstract: Topology and electron interactions are two central themes in modern condensed matter physics. Here we propose graphene based systems where both the band topology and interaction effects can be simply controlled with electric fields. We study a number of systems of twisted double layers with small twist angle where a moir\'e super-lattice is formed. Each layer is chosen to be either AB stacked bilayer graphene (BG), ABC stacked trilayer graphene (TG), or hexagonal boron nitride (h-BN). In these systems a vertical applied electric field enables control of the bandwidth, and interestingly also the Chern number. We find that the Chern numbers of the bands associated with each of the two microscopic valleys can be $\pm 0,\pm 1,\pm 2, \pm 3$ depending on the specific system and vertical electrical field. We show that these graphene moir\'e super-lattices are promising platforms to realize a number of fascinating many-body phenomena, including (Fractional) Quantum Anomalous Hall Effects. We also discuss conceptual similarities and implications for modeling twisted bilayer graphene systems.


title: Multilevel Models Allow Modular Specification of What and Where to Regularize, Especially in Small Area Estimation

abstract: Through the lense of multilevel model (MLM) specification and regularization, this is a connect-the-dots introductory summary of Small Area Estimation, e.g. small group prediction informed by a complex sampling design. While a comprehensive book is (Rao and Molina 2015), the goal of this paper is to get interested researchers up to speed with some current developments. We first provide historical context of two kinds of regularization: 1) the regularization 'within' the components of a predictor and 2) the regularization 'between' outcome and predictor. We focus on the MLM framework as it allows the analyst to flexibly control the targets of the regularization. The flexible control is useful when analysts want to overcome shortcomings in design-based estimates. We'll describe the precision deficiencies (high variance) typical of design-based estimates of small groups. We then highlight an interesting MLM example from (Chaudhuri and Ghosh 2011) that integrates both kinds of regularization (between and within). The key idea is to use the design-based variance to control the amount of 'between' regularization and prior information to regularize the components 'within' a predictor. The goal is to let the design-based estimate have authority (when precise) but defer to a model-based prediction when imprecise. We conclude by discussing optional criteria to incorporate into a MLM prediction and possible entrypoints for extensions.


title: Magnetization current and anomalous Hall effect for massive Dirac electrons

abstract: Existing investigations of the anomalous Hall effect i.e. a current flowing transverse to the electric field in the absence of an external magnetic field) are concerned with the transport current. However, for many applications one needs to know the total current, including its pure magnetization part. In this paper, we employ the two-dimensional massive Dirac equation to find the exact universal total current flowing along a potential step of arbitrary shape. For a spatially slowly varying potential we find the current density $\mathbf{j}(\vec r)$ and the energy distribution of the current density $\mathbf{j}^\varepsilon(\vec r)$. The latter turns out to be unexpectedly nonuniform, behaving like a $\delta$-function at the border of the classically accessible area at energy~$\varepsilon$. To demonstrate explicitly the difference between the magnetization and transport currents we consider the transverse shift of an electron ray in an electric field.


title: Improving CNN classifiers by estimating test-time priors

abstract: The problem of different training and test set class priors is addressed in the context of CNN classifiers. We compare two different approaches to estimating the new priors: an existing Maximum Likelihood Estimation approach (optimized by an EM algorithm or by projected gradient descend) and a proposed Maximum a Posteriori approach, which increases the stability of the estimate by introducing a Dirichlet hyper-prior on the class prior probabilities. Experimental results show a significant improvement on the fine-grained classification tasks using known evaluation-time priors, increasing the top-1 accuracy by 4.0% on the FGVC iNaturalist 2018 validation set and by 3.9% on the FGVCx Fungi 2018 validation set. Estimation of the unknown test set priors noticeably increases the accuracy on the PlantCLEF dataset, allowing a single CNN model to achieve state-of-the-art results and outperform the competition-winning ensemble of 12 CNNs. The proposed Maximum a Posteriori estimation increases the prediction accuracy by 2.8% on PlantCLEF 2017 and by 1.8% on FGVCx Fungi, where the existing MLE method would lead to a decrease accuracy.


title: Optimal control of thin liquid films and transverse mode effects

abstract: We consider the control of a three-dimensional thin liquid film on a flat substrate, inclined at a non-zero angle to the horizontal. Controls are applied via same-fluid blowing and suction through the substrate surface. We consider both overlying and hanging films, where the liquid lies above or below the substrate, respectively. We study the weakly nonlinear evolution of the system, which is governed by a forced Kuramoto--Sivashinsky equation in two space dimensions. The uncontrolled problem exhibits three ranges of dynamics depending on the incline of the substrate: stable flat film solution, bounded chaotic dynamics, or unbounded exponential growth of unstable transverse modes. We proceed with the assumption that we may actuate at every point on the substrate. The main focus is the optimal control problem, which we first study in the special case that the forcing may only vary in the spanwise direction. The structure of the Kuramoto--Sivashinsky equation allows the explicit construction of optimal controls in this case using the classical theory of linear quadratic regulators. Such controls are employed to prevent the exponential growth of transverse waves in the case of a hanging film, revealing complex dynamics for the streamwise and mixed modes. We then consider the optimal control problem in full generality, and prove the existence of an optimal control. For numerical simulations, we employ an iterative gradient descent algorithm. In the final section, we consider the effects of transverse mode forcing on the chaotic dynamics present in the streamwise and mixed modes for the case of a vertical film flow. Coupling through nonlinearity allows us to reduce the average energy in solutions without directly forcing the linearly unstable dominant modes.


title: Morphosyntactic Tagging with a Meta-BiLSTM Model over Context Sensitive Token Encodings

abstract: The rise of neural networks, and particularly recurrent neural networks, has produced significant advances in part-of-speech tagging accuracy. One characteristic common among these models is the presence of rich initial word encodings. These encodings typically are composed of a recurrent character-based representation with learned and pre-trained word embeddings. However, these encodings do not consider a context wider than a single word and it is only through subsequent recurrent layers that word or sub-word information interacts. In this paper, we investigate models that use recurrent neural networks with sentence-level context for initial character and word-based representations. In particular we show that optimal results are obtained by integrating these context sensitive representations through synchronized training with a meta-model that learns to combine their states. We present results on part-of-speech and morphological tagging with state-of-the-art performance on a number of languages.


abstract: This paper addresses an R(p,q)-deformed conformal Virasoro algebra with an arbitrary conformal dimension Delta. Wellknown deformations constructed in the literature are deduced as particular cases. Then, the special case of the conformal dimension Delta=1 is elucidated for its interesting properties. The R(p,q)-KdV equation, associated with the deformed Virasoro algebra, is also derived and discussed. Finally, the (p,q)-deformed energy-momentum tensor, consistent with the central extension term, is computed and analyzed.


title: The Roles of Supervised Machine Learning in Systems Neuroscience

abstract: Over the last several years, the use of machine learning (ML) in neuroscience has been rapidly increasing. Here, we review ML's contributions, both realized and potential, across several areas of systems neuroscience. We describe four primary roles of ML within neuroscience: 1) creating solutions to engineering problems, 2) identifying predictive variables, 3) setting benchmarks for simple models of the brain, and 4) serving itself as a model for the brain. The breadth and ease of its applicability suggests that machine learning should be in the toolbox of most systems neuroscientists.


title: Geodesically Equivalent Metrics on Homogenous Spaces

abstract: Two metrics on a manifold are geodesically equivalent if sets of their unparameterized geodesics coincide. In this paper we show that if two left $G$-invariant metrics of arbitrary signature on homogenous space $G/H$ are geodesically equivalent, they are affinely equivalent, i.e. they have the same Levi-Civita connection. We also prove that existence of non-proportional, geodesically equivalent, $G$-invariant metrics on homogenous space $G/H$ implies that their holonomy algebra cannot be full. We give an algorithm for finding all left invariant metrics geodesically equivalent to a given left invariant metric on a Lie group. Using that algorithm we prove that no two left invariant metric, of any signature, on sphere $S^3$ are geodesically equivalent. However, we present examples of Lie groups that admit geodesically equivalent, non-proportional, left-invariant metrics.


title: Sparse and Constrained Attention for Neural Machine Translation

abstract: In NMT, words are sometimes dropped from the source or generated repeatedly in the translation. We explore novel strategies to address the coverage problem that change only the attention transformation. Our approach allocates fertilities to source words, used to bound the attention each word can receive. We experiment with various sparse and constrained attention transformations and propose a new one, constrained sparsemax, shown to be differentiable and sparse. Empirical evaluation is provided in three languages pairs.


title: Petrov-Galerkin Method for Fully Distributed-Order Fractional Partial Differential Equations

abstract: Distributed-order PDEs are tractable mathematical models for complex multiscaling anomalous transport, where derivative orders are distributed over a range of values. We develop a fast and stable Petrov-Galerkin spectral method for such models by employing Jacobi \textit{poly-fractonomial}s and Legendre polynomials as temporal and spatial basis/test functions, respectively. By defining the proper underlying \textit{ distributed Sobolev} spaces and their equivalent norms, we prove the well-posedness of the weak formulation, and thereby carry out the corresponding stability and error analysis. We finally provide several numerical simulations to study the performance and convergence of proposed scheme.


title: Operational Dynamical Modeling of spin 1/2 relativistic particles: the Dirac equation and its classical limit

abstract: The formalism of Operational Dynamical Modeling [Phys. Rev. Lett. {\bf 109}, 190403 (2012)] is employed to analyze dynamics of spin half relativistic particles. We arrive at the Dirac equation from specially constructed relativistic Ehrenfest theorems by assuming that the coordinates and momenta do not commute. Forbidding creation of antiparticles and requiring the commutativity of the coordinates and momenta lead to classical Spohn's equation [Ann. Phys. {\bf 282}, 420 (2000)]. Moreover, Spohn's equation turns out to be the classical Koopman-von Neumann theory underlying the Dirac equation.


title: Stochastic modified equations for the asynchronous stochastic gradient descent

abstract: We propose a stochastic modified equations (SME) for modeling the asynchronous stochastic gradient descent (ASGD) algorithms. The resulting SME of Langevin type extracts more information about the ASGD dynamics and elucidates the relationship between different types of stochastic gradient algorithms. We show the convergence of ASGD to the SME in the continuous time limit, as well as the SME's precise prediction to the trajectories of ASGD with various forcing terms. As an application of the SME, we propose an optimal mini-batching strategy for ASGD via solving the optimal control problem of the associated SME.


title: Baryons under Strong Magnetic Fields or in Theories with Space-dependent $\theta$-term

abstract: Baryonic states are sufficiently complex to reveal physics that is hidden in the mesonic bound states. Using gauge/gravity correspondence we study analytically and numerically baryons in theories with space-dependent $\theta$-term, or theories under strong magnetic fields. Such holographic studies on baryons are accommodated in a generic analytic framework we develop for anisotropic theories, where their qualitative features are common irrespective of the source that triggers the anisotropy. We find that the distribution of the quarks forming the state, depends on the angle between the baryon and the anisotropic direction. Its shape is increasingly elliptic with respect to the strength of the field sourcing the anisotropy, counterbalancing the broken rotational invariance on the gluonic degrees of freedom. Strikingly, the baryons dissociate in stages with a process that depends on the proximity of the quarks to the anisotropic direction, where certain quark pairs abandon the bound state first, followed by the closest pairs to them as the temperature increases. This observation may also serve as a way to identify the nature of certain exotic states. Finally, we investigate holographic baryons with decreased number of quarks and explain why in theories under consideration the presence of anisotropy does not modify the universal stability condition in contrast to the usual trend.


title: Quadratic Jahn-Teller effect of fullerene anions

abstract: The quadratic Jahn-Teller effect of C$_{60}^{n-}$ ($n=$ 1-5) is investigated from the first principles. Employing the density functional theory calculations with hybrid functional, the quadratic vibronic coupling constants of C$_{60}^-$ were derived. The warping of the adiabatic potential energy surface of C$_{60}^-$ by the quadratic vibronic coupling is estimated about 2 meV, which is much smaller than the Jahn-Teller stabilization energy ($\approx$ 50 meV). Because of the selection rule and the vibronic reduction, the quadratic coupling slightly modifies the vibronic states of C$_{60}$ anions. Particularly, in the case of C$_{60}^{3-}$, parity and symmetry selection rule significantly reduces the effect of quadratic coupling on vibronic states. The present results confirm that the low-energy vibronic dynamics of C$_{60}^{n-}$ is of pseudorotational type.


title: The implementation of a Fast-Folding pipeline for long-period pulsar searching in the PALFA survey

abstract: The Pulsar Arecibo L-Band Feed Array (PALFA) survey, the most sensitive blind search for radio pulsars yet conducted, is ongoing at the Arecibo Observatory in Puerto Rico. The vast majority of the 180 pulsars discovered by PALFA have spin periods shorter than 2 seconds. Pulsar surveys may miss long-period radio pulsars due to the summing of a finite number of harmonic components in conventional Fourier analyses (typically $\sim$16), or due to the strong effect of red noise at low modulation frequencies. We address this reduction in sensitivity by using a time-domain search technique: the Fast-Folding Algorithm (FFA). We designed a program that implements a FFA-based search in the PALFA processing pipeline, and tested the efficiency of the algorithm by performing tests under both ideal, white noise conditions, as well as with real PALFA observational data. In the two scenarios, we show that the time-domain algorithm has the ability to outperform the FFT-based periodicity search implemented in the survey. We perform simulations to compare the previously reported PALFA sensitivity with that obtained using our new FFA implementation. These simulations show that for a pulsar having a pulse duty cycle of roughly 3%, the performance of our FFA pipeline exceeds that of our FFT pipeline for pulses with DM $\lesssim$ 40 pc cm$^{-3}$ and for periods as short as $\sim$500 ms, and that the survey sensitivity is improved by at least a factor of two for periods $\gtrsim$ 6 sec. Discoveries from the implementation of the algorithm in PALFA are also presented in this paper.


title: Measurement and Significance of Wilson Loops in Synthetic Gauges Fields

abstract: We study Wilson loops as a necessary tool for unambiguous identification of non-Abelian synthetic gauge fields, with attention to certain crucial but often overlooked features, such as the requirement of at least three distinct loops. We devise a method to determine the complete Wilson loop matrix from the time evolved amplitudes of the internal atomic states of laser-coupled ultracold atoms that does not require lattice confinement. The analysis is done in the context of a new cyclic model that can realize both Abelian and non-Abelian structures within a single configuration with continuous variation possible between U(1) and U(2) gauge groups by varying the detuning of the laser fields.


title: Classifier-agnostic saliency map extraction

abstract: Currently available methods for extracting saliency maps identify parts of the input which are the most important to a specific fixed classifier. We show that this strong dependence on a given classifier hinders their performance. To address this problem, we propose classifier-agnostic saliency map extraction, which finds all parts of the image that any classifier could use, not just one given in advance. We observe that the proposed approach extracts higher quality saliency maps than prior work while being conceptually simple and easy to implement. The method sets the new state of the art result for localization task on the ImageNet data, outperforming all existing weakly-supervised localization techniques, despite not using the ground truth labels at the inference time. The code reproducing the results is available at https://github.com/kondiz/casme . The final version of this manuscript is published in Computer Vision and Image Understanding and is available online at https://doi.org/10.1016/j.cviu.2020.102969 .


title: The Geometric Blueprint of Perovskites

abstract: Perovskite minerals form an essential component of the Earth's mantle, and synthetic crystals are ubiquitous in electronics, photonics, and energy technology. The extraordinary chemical diversity of these crystals raises the question on how many and which perovskites are yet to be discovered. Here we show that the "no-rattling" principle postulated by Goldschmidt in 1926, describing the geometric conditions under which a perovskite can form, is much more effective than previously thought, and allows us to predict new perovskites with a fidelity of 80%. By supplementing this principle with inferential statistics and internet data mining we establish that currently known perovskites are only the tip of the iceberg, and we enumerate ninety thousand hitherto-unknown compounds awaiting to be studied. Our results suggest that geometric blueprints may enable the systematic screening of millions of compounds, and offer untapped opportunities in structure prediction and materials design.


title: Asymmetric Dark Matter, Inflation and Leptogenesis from B-L Symmetry Breaking

abstract: We propose a unified setup for dark matter, inflation and baryon asymmetry generation through the neutrino mass seesaw mechanism. Our scenario emerges naturally from an extended gauge group containing $B-L$ as a non-commutative symmetry, broken by a singlet scalar that also drives inflation. Its decays reheat the universe, producing the lightest right-handed neutrino. Automatic matter parity conservation leads to the stability of an asymmetric dark matter candidate, directly linked to the matter-antimatter asymmetry in the universe.


title: Tale of stable interacting dark energy, observational signatures, and the $H_0$ tension

abstract: We investigate the observational consequences of a novel class of stable interacting dark energy (IDE) models, featuring interactions between dark matter (DM) and dark energy (DE). In the first part of our work, we start by considering two IDE models which are known to present early-time linear perturbation instabilities. Applying a transformation depending on the dark energy equation of state (EoS) to the DM-DE coupling, we then obtain two novel stable IDE models. Subsequently, we derive robust and accurate constraints on the parameters of these models, assuming a constant EoS $w_x$ for the DE fluid, in light of some of the most recent publicly available cosmological data. These include Cosmic Microwave Background (CMB) temperature and polarization anisotropy measurements from the \textit{Planck} satellite, a selection of Baryon Acoustic Oscillation measurements, Supernovae Type-Ia luminosity distance measurements from the JLA sample, and measurements of the Hubble parameter up to redshift $2$ from cosmic chronometers. Our analysis displays a mild preference for the DE fluid residing in the phantom region ($w_x<-1$), with significance up to 95\% confidence level, while we obtain new upper limits on the coupling parameter between the dark components. The preference for a phantom DE suggests a coupling function $Q<0$, thus a scenario where energy flows from the DE to the DM. We also examine the possibility of addressing the $H_0$ and $\sigma_8$ tensions, finding that only the former can be partially alleviated. Finally, we perform a Bayesian model comparison analysis to quantify the possible preference for the two IDE models against the standard concordance $\Lambda$CDM model, finding that the latter is always preferred with the strength of the evidence ranging from positive to very strong.


title: NOEMA Observations of a Molecular Cloud in the low-metallicity Galaxy Kiso 5639

abstract: A giant star-forming region in a metal-poor dwarf galaxy has been observed in optical lines with the 10-m Gran Telescopio Canarias and in the emission line of CO(1-0) with the NOEMA mm-wave interferometer. The metallicity was determined to be 12+log(O/H)=7.83+-0.09, from which we estimate a conversion factor of alpha_CO~100 Msun/pc2/(K km/s) and a molecular cloud mass of ~2.9x10^7 Msun. This is an enormous concentration of molecular mass at one end of a small galaxy, suggesting a recent accretion. The molecular cloud properties seem normal: the surface density, 120 Msun/pc2, is comparable to that of a standard giant molecular cloud, the cloud's virial ratio of ~1.8 is in the star-formation range, and the gas consumption time, 0.5 Gyr, at the present star formation rate is typical for molecular regions. The low metallicity implies that the cloud has an average visual extinction of only 0.8 mag, which is close to the threshold for molecule formation. With such an extinction threshold, molecular clouds in metal-poor regions should have high surface densities and high internal pressures. If high pressure is associated with the formation of massive clusters, then metal-poor galaxies such as dwarfs in the early universe could have been the hosts of metal-poor globular clusters.


title: Sample Compression for Real-Valued Learners

abstract: We give an algorithmically efficient version of the learner-to-compression scheme conversion in Moran and Yehudayoff (2016). In extending this technique to real-valued hypotheses, we also obtain an efficient regression-to-bounded sample compression converter. To our knowledge, this is the first general compressed regression result (regardless of efficiency or boundedness) guaranteeing uniform approximate reconstruction. Along the way, we develop a generic procedure for constructing weak real-valued learners out of abstract regressors; this may be of independent interest. In particular, this result sheds new light on an open question of H. Simon (1997). We show applications to two regression problems: learning Lipschitz and bounded-variation functions.


title: Algorithmic and algebraic aspects of unshuffling permutations

abstract: A permutation is said to be a square if it can be obtained by shuffling two order-isomorphic patterns. The definition is intended to be the natural counterpart to the ordinary shuffle of words and languages. In this paper, we tackle the problem of recognizing square permutations from both the point of view of algebra and algorithms. On the one hand, we present some algebraic and combinatorial properties of the shuffle product of permutations. We follow an unusual line consisting in defining the shuffle of permutations by means of an unshuffling operator, known as a coproduct. This strategy allows to obtain easy proofs for algebraic and combinatorial properties of our shuffle product. We besides exhibit a bijection between square $(213,231)$-avoiding permutations and square binary words. On the other hand, by using a pattern avoidance criterion on directed perfect matchings, we prove that recognizing square permutations is {\bf NP}-complete.


title: Evolving Real-Time Heuristics Search Algorithms with Building Blocks

abstract: The research area of real-time heuristics search has produced quite many algorithms. In the landscape of real-time heuristics search research, it is not rare to find that an algorithm X that appears to perform better than algorithm Y on a group of problems, performed worse than Y for another group of problems. If these published algorithms are combined to generate a more powerful space of algorithms, then that novel space of algorithms may solve a distribution of problems more efficiently. Based on this intuition, a recent work Bulitko 2016 has defined the task of finding a combination of heuristics search algorithms as a survival task. In this evolutionary approach, a space of algorithms is defined over a set of building blocks published algorithms and a simulated evolution is used to recombine these building blocks to find out the best algorithm from that space of algorithms. In this paper, we extend the set of building blocks by adding one published algorithm, namely lookahead based A-star shaped local search space generation method from LSSLRTA-star, plus an unpublished novel strategy to generate local search space with Greedy Best First Search. Then we perform experiments in the new space of algorithms, which show that the best algorithms selected by the evolutionary process have the following property: the deeper is the lookahead depth of an algorithm, the lower is its suboptimality and scrubbing complexity.


title: Probing the Shape of Quantum Surfaces: the Quadrupole Moment Operator

abstract: The standard toolkit of operators to probe quanta of geometry in loop quantum gravity consists in area and volume operators as well as holonomy operators. New operators have been defined, in the U(N) framework for intertwiners, which allow to explore the finer structure of quanta of geometry. However these operators do not carry information on the global shape of the intertwiners. Here we introduce dual multipole moments for continuous and discrete surfaces, defined through the normal vector to the surface, taking special care to maintain parametrization invariance. These are raised to multipole operators probing the shape of quantum surfaces. Further focusing on the quadrupole moment, we show that it appears as the Hessian matrix of the large spin Gaussian approximation of coherent intertwiners, which is the standard method for extracting the semi-classical regime of spinfoam transition amplitudes. This offers an improvement on the usual loop quantum gravity techniques, which mostly focus on the volume operator, in the perspective of modeling (quantum) gravitational waves as shape fluctuations waves propagating on spin network states.


title: Does a distinct quasi many-body localized phase exist? A numerical study of a translationally invariant system in the thermodynamic limit

abstract: We consider a quench in an infinite spin ladder describing a system with two species of bosons in the limit of strong interactions. If the heavy bosonic species has infinite mass the model becomes a spin chain with quenched binary disorder which shows true Anderson localization (AL) or many-body localization (MBL). For finite hopping amplitude $J'$ of the heavy particles, on the other hand, we find an exponential polarization decay with a relaxation rate which depends monotonically on $J'$. Furthermore, the entanglement entropy changes from a constant (AL) or logarithmic (MBL) scaling in time $t$ for $J'=0$ to a sub-ballistic power-law, $S_{\textrm{ent}}\sim t^\alpha$ with $\alpha< 1$, for finite $J'$. We do not find a distinct regime in time where the dynamics for $J'\neq 0$ shows the characteristics of an MBL phase. Instead, we discover a time regime with distinct dephasing and entanglement times, different both from a localized and a fully ergodic phase.

categories: cond-mat.dis-nn cond-mat.quant-gas cond-mat.stat-mech cond-mat.str-el


title: An imaging algorithm for a lunar orbit interferometer array

abstract: Radio astronomical observation below 30 MHz is hampered by the refraction and absorption of the ionosphere, and the radio frequency interference (RFI), so far high angular resolution sky intensity map is not available. An interferometer array on lunar orbit provides a perfect observatory in this frequency band: it is out of the ionosphere and the Moon helps to block the RFIs from the Earth. The satellites can make observations on the far side of the Moon and then send back the data on the near side part of the orbit. However, for such array the traditional imaging algorithm is not applicable: the field of view is very wide (almost whole sky), and for baselines distributed on a plane, there is a mirror symmetry between the two sides of the plane. A further complication is that for each baseline, the Moon blocks part of the sky, but as the satellites orbit the Moon, both the direction of the baseline and the blocked sky change, so even imaging algorithms which can deal with non-coplanar baseline may not work in this case. Here we present an imaging algorithm based on solving the linear mapping equations relating the sky intensity to the visibilities. We show that the mirror symmetry can be broken by the three dimensional baseline distribution generated naturally by the precession of the orbital plane of the satellites. The algorithm is applicable and good maps could be reconstructed, even though for each baseline the sky blocking by the Moon is different. We also investigate how the map-making is affected by inhomogeneous baseline distributions.


title: E6 inspired SUSY models with Custodial Symmetry

abstract: The breakdown of E_6 within the supersymmetric (SUSY) Grand Unified Theories (GUTs) can result in SUSY extensions of the standard model (SM) based on the SM gauge group together with extra U(1) gauge symmetry under which right-handed neutrinos have zero charge. In these U(1)_N extensions of the minimal supersymmetric standard model (MSSM) a single discrete \tilde{Z}^H_2 symmetry may be used to suppress the most dangerous operators, that give rise to proton decay as well as non-diagonal flavour transitions at low energies. The SUSY models under consideration involves Z' and extra exotic matter beyond the MSSM. We discuss leptogenesis within this SUSY model and argue that the extra exotic states may lead to the non--standard Higgs decays.


title: Mathematics of Smoothed Particle Hydrodynamics: a Study via Nonlocal Stokes Equations

abstract: Smoothed Particle Hydrodynamics (SPH) is a popular numerical technique developed for simulating complex fluid flows. Among its key ingredients is the use of nonlocal integral relaxations to local differentiations. Mathematical analysis of the corresponding nonlocal models on the continuum level can provide further theoretical understanding of SPH. We present, in this part of a series of works on the mathematics of SPH, a nonlocal relaxation to the conventional linear steady state Stokes system for incompressible viscous flows. The nonlocal continuum model is characterized by a smoothing length $\delta$ which measures the range of nonlocal interactions. It serves as a bridge between the discrete approximation schemes that involve a nonlocal integral relaxation and the local continuum models. We show that for a class of carefully chosen nonlocal operators, the resulting nonlocal Stokes equation is well-posed and recovers the original Stokes equation in the local limit when $\delta$ approaches zero. We also discuss the implications of our finding on the design of numerical methods.


title: Magnetostatic problems in fractal domains

abstract: We consider a magnetostatic problem in a 3D "cylindrical" domain of Koch type. We prove existence and uniqueness results for both the fractal and pre-fractal problems and we investigate the convergence of the pre-fractal solutions to the limit fractal one. We consider the numerical approximation of the pre-fractal problems via FEM and we prove a priori error estimates. Some numerical simulations are also shown. Our long term motivation includes studying problems that appear in quantum physics in fractal domains.

categories: math-ph math.AP math.FA math.MG math.MP math.NA


title: Learning What Information to Give in Partially Observed Domains

abstract: In many robotic applications, an autonomous agent must act within and explore a partially observed environment that is unobserved by its human teammate. We consider such a setting in which the agent can, while acting, transmit declarative information to the human that helps them understand aspects of this unseen environment. In this work, we address the algorithmic question of how the agent should plan out what actions to take and what information to transmit. Naturally, one would expect the human to have preferences, which we model information-theoretically by scoring transmitted information based on the change it induces in weighted entropy of the human's belief state. We formulate this setting as a belief MDP and give a tractable algorithm for solving it approximately. Then, we give an algorithm that allows the agent to learn the human's preferences online, through exploration. We validate our approach experimentally in simulated discrete and continuous partially observed search-and-recover domains. Visit http://tinyurl.com/chitnis-corl-18 for a supplementary video.


title: The Superconducting and Pseudogap Phase Diagram of High-Tc Cuprates

abstract: We derive analytic expressions for the critical temperatures of the superconducting (SC) and pseudogap (PG) phases of the high-Tc cuprates, which are in excellent agreement with the experimental data for single-layered materials such as LSCO, Bi2201 and Hg1201. Our effective Hamiltonian, defined in the oxygen square sub-lattices formed by the alternate hybridization of $p_x$ and $p_y$ orbitals with the $3d$ copper orbitals, provides an unified explanation for the $d_{x^2-y^2}$ symmetry of both the SC and PG order parameters. Attractive and repulsive interactions involve holes of the two different sublattices and can be derived from the spin-fermion model. Optimal doping occurs when the chemical potential vanishes. For $N$-layered cuprates, the growth of the optimal temperature with $N$, as well as the trend of the SC and AF domes to superimpose, can be simply understood. Our results for the optimal SC transition temperature are in excellent agreement with the experiments for $N=2$ materials of the $Bi$ and $Hg$ families. For $N=3$ the agreement is still satisfactory, while for $N>3$, it becomes poor. The explanation for these facts allows us to suggest a method for increasing the critical SC temperature in cuprates.


title: Extracting 21cm signal by frequency and angular filtering

abstract: Extracting the neutral hydrogen (HI) signal is a great challenge for cosmological 21cm experiments, both the astrophysical foregrounds and the receiver noise are typically several orders of magnitude greater than the 21cm signal. However, the different properties of the 21cm signal, foreground, and noise can be exploited to separate these components. The foregrounds are generally smooth or correlated over the frequency space along a line of sight (l.o.s.), while both the 21cm signal and the noise varies stochastically along the same l.o.s. The foreground can be removed by filtering out the smooth component in the frequency space. The receiver noise is basically uncorrelated for observations at different times, hence for surveys they are also uncorrelated in the different directions, while the 21cm signal which traces the large scale structure are correlated up to certain scales. In this exercise, we apply Wiener filters in frequency and angular space, to extract the 21cm signals. We found that the method works well. The inaccurate knowledge about beam could degrade the reconstruction, but the overall result is still good, showing that the method is fairly robust.


title: On the Selection of Initialization and Activation Function for Deep Neural Networks

abstract: The weight initialization and the activation function of deep neural networks have a crucial impact on the performance of the training procedure. An inappropriate selection can lead to the loss of information of the input during forward propagation and the exponential vanishing/exploding of gradients during back-propagation. Understanding the theoretical properties of untrained random networks is key to identifying which deep networks may be trained successfully as recently demonstrated by Schoenholz et al. (2017) who showed that for deep feedforward neural networks only a specific choice of hyperparameters known as the `edge of chaos' can lead to good performance. We complete this analysis by providing quantitative results showing that, for a class of ReLU-like activation functions, the information propagates indeed deeper for an initialization at the edge of chaos. By further extending this analysis, we identify a class of activation functions that improve the information propagation over ReLU-like functions. This class includes the Swish activation, $\phi_{swish}(x) = x \cdot \text{sigmoid}(x)$, used in Hendrycks & Gimpel (2016), Elfwing et al. (2017) and Ramachandran et al. (2017). This provides a theoretical grounding for the excellent empirical performance of $\phi_{swish}$ observed in these contributions. We complement those previous results by illustrating the benefit of using a random initialization on the edge of chaos in this context.


title: Quantitative spectral classification of Galactic O stars

abstract: Our goal is to provide a quantification of several spectral classification criteria for O stars. We collected high-spectral resolution spectra of 105 Galactic O-type stars from various archives. We measured equivalent widths of classification lines. We defined average values of classification criteria for given spectral types and luminosity classes. We find that the ratio HeI 4471 to HeII 4542 well matches the published ratios for spectral types. We have quantified equivalent width ratios of helium and silicon lines among O8-O9.7 stars to refine spectral class typing in this spectral range. We present quantitative criteria to separate between luminosity class V, IV-III-II (grouped), and I among O3-O8.5 stars, mainly based on the strength of HeII 4686. We find that these criteria also define very well the f, (f), and ((f)) classes for O3-O7.5 stars. Among O9-O9.7 stars we quantify the ratios of HeII 4686 to HeI 4713 and SiIV 4089 to HeI 4026 for all luminosity classes. The tabulated values of the classification criteria should help classify any new O-type stars. The final step of the classification process should rely on a direct comparison to standard stars of the assigned spectral type and luminosity class.


title: Optimal Sketching Bounds for Exp-concave Stochastic Minimization

abstract: We derive optimal statistical and computational complexity bounds for exp-concave stochastic minimization in terms of the effective dimension. For common eigendecay patterns of the population covariance matrix, this quantity is significantly smaller than the ambient dimension. Our results reveal interesting connections to sketching results in numerical linear algebra. In particular, our statistical analysis highlights a novel and natural relationship between algorithmic stability of empirical risk minimization and ridge leverage scores, which play significant role in sketching-based methods. Our main computational result is a fast implementation of a sketch-to-precondition approach in the context of exp-concave empirical risk minimization.


title: Characterization of Argon Plasma in a variable Multi-pole line Cusp Magnetic Field Configuration

abstract: This paper demonstrates a detailed characterization of argon plasma in a variable multi-pole line cusp magnetic field (VMMF). The VMMF has been produced by placing six electromagnets (with embedded profiled vacoflux-50 core) over a large cylindrical volume (1 m axial length and 40 cm diameter). The magnetic field have been measured by hall probe method and compared with simulated magnetic field by performing simulation using FEMM tools. Results from magnetic field simulation indicate that the rate of change of pole magnetic field (maximum magnetic field) with respect to magnet current for vacoflux-50 core is high (7.53 G/A) as compared to the simple air core electromagnet (2.15 G/A). The area of the nearly field free region (null region) in the chamber volume can be controlled without changing a number of pole magnets. From the experimental results, it has been observed that in this field configuration the confinement of the primary electrons increases and leak width of plasma decreases with increasing the magnetic field. Thus the mean density, particle confinement time and the stability of the plasma increase with increasing magnetic field. In addition to this, it has been found that the radial uniformity of the plasma density explicitly depends on the VMMF. It is also shown that the VMMF controls the scavenging of confined primary electrons and confinement of primary electron increased with magnetic field which helps to boost up the plasma density.


title: Colouring Square-Free Graphs without Long Induced Paths

abstract: The complexity of {\sc Colouring} is fully understood for $H$-free graphs, but there are still major complexity gaps if two induced subgraphs $H_1$ and $H_2$ are forbidden. Let $H_1$ be the $s$-vertex cycle $C_s$ and $H_2$ be the $t$-vertex path $P_t$. We show that {\sc Colouring} is polynomial-time solvable for $s=4$ and $t\leq 6$, strengthening several known results. Our main approach is to initiate a study into the boundedness of the clique-width of atoms (graphs with no clique cutset) of a hereditary graph class. We first show that the classifications of boundedness of clique-width of $H$-free graphs and $H$-free atoms coincide. We then show that this is not the case if two graphs are forbidden: we prove that $(C_4,P_6)$-free atoms have clique-width at most~18. Our key proof ingredients are a divide-and-conquer approach for bounding the clique-width of a subclass of $C_4$-free graphs and the construction of a new bound on the clique-width for (general) graphs in terms of the clique-width of recursively defined subgraphs induced by homogeneous pairs and triples of sets. As a complementary result we prove that {\sc Colouring} is \NP-complete for $s=4$ and $t\geq 9$, which is the first hardness result on {\sc Colouring} for $(C_4,P_t)$-free graphs. Combining our new results with known results leads to an almost complete dichotomy for \cn restricted to $(C_s,P_t)$-free graphs.