title: Predictive power of grand unification from quantum gravity

abstract: If a grand-unified extension of the asymptotically safe Reuter fixed-point for quantum gravity exists, it determines free parameters of the grand-unified scalar potential. All quartic couplings take their fixed-point values in the trans-Planckian regime. They are irrelevant parameters that are, in principle, computable for a given particle content of the grand unified model. In turn, the direction of spontaneous breaking of the grand-unified gauge symmetry becomes predictable. For the flow of the couplings below the Planck mass, gauge and Yukawa interactions compete for the determination of the minimum of the effective potential.


title: An SFP--FCC Method for Pricing and Hedging Early-exercise Options under L\'evy Processes

abstract: This paper extends the Singular Fourier--Pad\'e (SFP) method proposed by Chan (2018) to pricing/hedging early-exercise options--Bermudan, American and discrete-monitored barrier options--under a L\'evy process. The current SFP method is incorporated with the Filon--Clenshaw--Curtis (FCC) rules invented by Dom\'inguez et al. (2011), and we call the new method SFP--FCC. The main purpose of using the SFP--FCC method is to require a small number of terms to yield fast error convergence and to formulate option pricing and option Greek curves rather than individual prices/Greek values. We also numerically show that the SFP--FCC method can retain a global spectral convergence rate in option pricing and hedging when the risk-free probability density function is piecewise smooth. Moreover, the computational complexity of the method is $\mathcal{O}((L-1)(N+1)(\tilde{N} \log \tilde{N}) )$ with $N$ a (small) number of complex Fourier series terms, $\tilde{N}$ a number of Chebyshev series terms and $L$, the number of early-exercise/monitoring dates. Finally, we show that our method is more favourable than existing techniques in numerical experiments.


title: Computing spatially resolved rotational hydration entropies from atomistic simulations

abstract: For a first principles understanding of macromolecular processes, a quantitative understanding of the underlying free energy landscape and in particular its entropy contribution is crucial. The stability of biomolecules, such as proteins, is governed by the hydrophobic effect, which arises from competing enthalpic and entropic contributions to the free energy of the solvent shell. While the statistical mechanics of liquids, as well as molecular dynamics simulations have provided much insight, solvation shell entropies remain notoriously difficult to calculate, especially when spatial resolution is required. Here, we present a method that allows for the computation of spatially resolved rotational solvent entropies via a non-parametric k-nearest-neighbor density estimator. We validated our method using analytic test distributions and applied it to atomistic simulations of a water box. With an accuracy of better than 9.6%, the obtained spatial resolution should shed new light on the hydrophobic effect and the thermodynamics of solvation in general.

categories: physics.bio-ph physics.chem-ph physics.comp-ph stat.CO


title: PSTGF : time-independent R-Matrix atomic electron-impact code

abstract: STGF is a community code employed for outer-region R-matrix calculations, describing electron-impact collisional processes. It is widely recognised that the original version of STGF was written by M. J. Seaton in 1983, but through constant refinement over the next decades by worldwide contributors has evolved into its current form that more reflects modern coding practice and current computer architectures. Despite its current wide acceptance, it was never formally published. Therefore, we present an updated high-performance parallel version of PSTGF, that balances the requirements of small university clusters, yet can exploit the computational power of cutting edge supercomputers. There are many improvements over the original STGF, but most noticeably, the full introduction of MQDT options that provide subsequent integration with ICFT (Intermediate Coupling Frame Transformation) codes, and for either Breit-Pauli / DARC (Dirac Atomic R-matrix Codes), better load balancing, high levels of vectorisation and simplified output. Semantically, the program is full fortran 90 in conjunction with MPI (Message Passing Interface) though has CUDA fortran options for the most numerically intensive code sections.


title: Subdiffusion in one-dimensional Hamiltonian chains with sparse interactions

abstract: We establish rigorously that transport is slower than diffusive for a class of disordered one-dimensional Hamiltonian chains. This is done by deriving quantitative bounds on the variance in equilibrium of the energy or particle current, as a function of time. The slow transport stems from the presence of rare insulating regions (Griffiths regions). In many-body disordered quantum chains, they correspond to regions of anomalously high disorder, where the system is in a localized phase. In contrast, we deal with quantum and classical disordered chains where the interactions, usually referred to as anharmonic couplings in classical systems, are sparse. The system hosts thus rare regions with no interactions and, since the chain is Anderson localized in the absence of interactions, the non-interacting rare regions are insulating. Part of the mathematical interest of our model is that it is one of the few non-integrable models where the diffusion constant can be rigorously proven not to be infinite.


title: Surprisingly Strong K-band Emission Found in Low Luminosity Active Galactic Nuclei

abstract: We examine the near-infrared (NIR) emission from low-luminosity AGNs (LLAGNs). Our galaxy sample includes 15 objects with detected 2-10 keV X-ray emission, dynamical black hole mass estimates from the literature, and available Gemini/NIFS integral field spectroscopy (IFU) data. We find evidence for red continuum components at the center of most galaxies, consistent with the hot dust emission seen in higher luminosity AGN. We decompose the spectral data cubes into a stellar and continuum component, assuming the continuum component comes from thermal emission from hot dust. We detect nuclear thermal emission in 14 out of 15 objects. This emission causes weaker CO absorption lines and redder continuum ($2.05-2.28\:\mu$m) in our $K$-band data, as expected from hot dust around an AGN. The NIR emission is clearly correlated with the 2-10 keV X-ray flux, with a Spearman coefficient of $r_{spearman}=0.69$ suggesting a $>99\%$ significance of correlation, providing further evidence of an AGN origin. Our sample has typical X-ray and NIR fluxes $3-4$ orders of magnitude less luminous than previous work studying the NIR emission from AGN. We find that the ratio of NIR to X-ray emission increases towards lower Eddington ratios. The NIR emission in our sample is often brighter than the X-ray emission, with our $K$-band AGN luminosities comparable to or greater than the 2-10 keV X-ray luminosities in all objects with Eddington ratios below $0.01\%$. The nature of this LLAGN NIR emission remains unclear, with one possibility being an increased contribution from jet emission at these low luminosities. These observations suggest JWST will be a useful tool for detecting the lowest luminosity AGN.


title: Dirichlet Depths for Point Process

abstract: Statistical depths have been well studied for multivariate and functional data over the past few decades, but remain under-explored for point processes. A first attempt on the notion of point process depth was conducted recently where the depth was defined as a weighted product of two terms: (1) the probability of the number of events in each process and (2) the depth of the event times conditioned on the number of events by using a Mahalanobis depth. We point out that multivariate depths such as the Mahalanobis depth cannot be directly used because they often neglect the important ordered property in the point process events. To deal with this problem, we propose a model-based approach for point processes systematically. In particular, we develop a Dirichlet-distribution-based framework on the conditional depth term, where the new methods are referred to as Dirichlet depths. We examine the mathematical properties of the new depths and conduct the asymptotic analysis. In addition, we illustrate the new methods using various simulated and real experiment data. It is found that the proposed framework provides a proper center-outward rank and the new methods have superior decoding performance to previous methods in two neural spike train datasets.


title: Magnetic topology in fluids

abstract: We study the evolution of turbulent magnetic fields from a topological point of view, invoking commonplace mathematical tools from general topology and dynamical systems theory which connect magnetic field evolution to time reversal invariance, entropy increase and the second law of thermodynamics. We show that in fact magnetic topology is well-defined only in the phase space corresponding to a dynamical system governed by the induction equation. Hence the field's topology and stochasticity can be studied in terms of the corresponding phase space trajectories rather than the field lines in real Euclidean space. In fact, our results suggest that magnetic field lines should not be taken too literally because their existence and uniqueness and more importantly continuity in time require strong mathematical conditions, hardly satisfied in astrophysical systems. As for magnetic topology change, it is shown that the phase space topology is preserved in time for a magnetic field which, besides satisfying few continuity conditions, solves a time reversal invariant induction equation. What breaks the time symmetry in the induction equation is the presence of non-ideal plasma effects at small scales such as resistivity, which results from random collisions between diffusing electrons and other particles. The small scale, stochastic disturbances produced thereby are super-linearly amplified by Richardson diffusion in the turbulent cascade, which are eventually manifested as large scale reconnection events, somehow similar to stretching quantum fluctuations during inflation to seed large scale cosmological structures. This suggests that reconnection is rooted in the second law of thermodynamics that dictates entropy increase which in turn breaks the time symmetry.


title: Multitype Integer Monoid Optimization and Applications

abstract: Configuration integer programs (IP) have been key in the design of algorithms for NP-hard high-multiplicity problems since the pioneering work of Gilmore and Gomory [Oper. Res., 1961]. Configuration IPs have a variable for each possible configuration, which describes a placement of items into a location, and whose value corresponds to the number of locations with that placement. In high multiplicity problems items come in types, and are represented succinctly by a vector of multiplicities; solving the configuration IP then amounts to deciding whether the input vector of multiplicities of items of each type can be decomposed into a given number of configurations. We make this implicit notion explicit by observing that the set of all input vectors decomposable into configurations forms a monoid, and solving the configuration IP is the Monoid Decomposition problem. Motivated by applications, we enrich this problem in two ways. First, sometimes each configuration additionally has an objective value, yielding an optimization problem of finding a "best" decomposition under the given objective. Second, there are often different types of configurations for different types of locations. The resulting problem is to optimize over decompositions of the input multiplicity vector into configurations of several types, and we call it Multitype Integer Monoid Optimization, or MIMO. We develop fast exact algorithms for various MIMO with few or many location types and with various objectives. Our algorithms build on a novel proximity theorem connecting the solutions of a certain configuration IP to those of its continuous relaxation. We then cast several fundamental scheduling and bin packing problems as MIMOs, and thereby obtain new or substantially faster algorithms for them. We complement our positive algorithmic results by hardness results.


title: Resonant $a_0(980)$ state in triangle rescattering $D_s^+\to \pi^+\pi^0\eta$ decays

abstract: We study the $D_s^+\to \pi^+(a_0(980)^0\to)\pi^0\eta$, $\pi^0(a_0(980)^+\to)\pi^+\eta$ decays, which have been recently measured by the BESIII collaboration. We propose that $D_s^+\to \pi^{+(0)}(a_0(980)^{0(+)}\to)\pi^{0(+)}\eta$ receives the contributions from the triangle rescattering processes, where $M^0$ and $\rho^+$ in $D_s^+\to M^0 \rho^+$, by exchanging $\pi^{0(+)}$, are formed as $a_0(980)^{0(+)}$ and $\pi^{+(0)}$, respectively, with $M^0=(\eta,\eta')$. Accordingly, we calculate that ${\cal B}(D_s^+\to a_0(980)^{0(+)}\pi^{+(0)})=(1.7\pm 0.2\pm 0.1)\times 10^{-2}$ and ${\cal B}(D_s^+\to \pi^{+(0)}(a_0(980)^{0(+)}\to)\pi^{0(+)}\eta)=(1.4\pm 0.1\pm 0.1)\times 10^{-2}$, being consistent with the data.


title: Learning Invariants through Soft Unification

abstract: Human reasoning involves recognising common underlying principles across many examples. The by-products of such reasoning are invariants that capture patterns such as "if someone went somewhere then they are there", expressed using variables "someone" and "somewhere" instead of mentioning specific people or places. Humans learn what variables are and how to use them at a young age. This paper explores whether machines can also learn and use variables solely from examples without requiring human pre-engineering. We propose Unification Networks, an end-to-end differentiable neural network approach capable of lifting examples into invariants and using those invariants to solve a given task. The core characteristic of our architecture is soft unification between examples that enables the network to generalise parts of the input into variables, thereby learning invariants. We evaluate our approach on five datasets to demonstrate that learning invariants captures patterns in the data and can improve performance over baselines.


title: Nonlinear response of wrinkled premixed flames to time- and space-dependent forcing and stretch

abstract: Premixed-flame wrinkling is studied via a Michelson-Sivashinsky (MS) type of evolution equation retaining the Darrieus-Landau (DL) instability, a curvature effect and a geometric nonlinearity. Here it also keeps forcing by longitudinal shearflow and wrinkle stretch by transverse flow; both imposed stimuli vary in time and space as to make the front slope comprise a given fluctuating spatial harmonics and unknown pole-decomposed pieces. A DL-free Burgers version is examined in parallel, also with Neumann conditions and symmetry. As is shown for both models, solving Ntot equations of motion for the poles in principle yields the front dynamics, the arclength increment V(t) and its time-average Vav. Yet this could be worked out analytically (or nearly so) only in high-frequency HF or low-frequency LF limits. These tackle one or two pairs of poles per cell, then a large number of pairs Ntot forming two piles viewed as continua, one per crest. Despite ample pole motions that make some commute between crests, Vav grows in a nearly parabolic way with the combined intensity of forcing and stretch. LF stimuli and DL instability can induce multiple branches and relaxation phenomena. Numerical t-averages are needed even if V(t) is analytically known. For Ntot=1,2 and short wrinkles, or Burgers fronts, Vav transitions from quadratic to sublinear as the forcing grows ; for longer wrinkles Vav keeps its MS value at moderate forcing, then bifurcates to an ultimately sublinear growth that depends on the stimulus phases. For very long wrinkles, coupled integral equations give analytical slope and pole-density profiles, but pile heights/contents need a t-dependent numerical search of up to two roots to get V(t). A summary, a discussion and hints of generalizations are provided, and open problems are evoked.


title: Knowledge Discovery In Nanophotonics Using Geometric Deep Learning

abstract: We present here a new approach for using the intelligence aspects of artificial intelligence for knowledge discovery rather than device optimization in electromagnetic (EM) nanostructures. This approach uses training data obtained through full-wave EM simulations of a series of nanostructures to train geometric deep learning algorithms to assess the range of feasible responses as well as the feasibility of a desired response from a class of EM nanostructures. To facilitate the knowledge discovery and reduce the computation complexity, our approach combines the dimensionality reduction technique (using an autoencoder) with convex-hull and one-class support-vector-machine (SVM) algorithms to find the range of the feasible responses in the latent (or the reduced) response space of the EM nanostructure. We show that by using a small set of training instances (compared to all possible structures), our approach can provide better than 95% accuracy in assessing the feasibility of a given response. More importantly, the one-class SVM algorithm can be trained to provide the degree of feasibility (or unfeasibility) of a response from a given nanostructure. This important information can be used to modify the initial structure to an alternative one that can enable an initially unfeasible response. To show the applicability of our approach, we apply it to two important classes of binary metasurfaces (MSs), formed by array of plasmonic nanostructures, and periodic MSs formed by an array of dielectric nanopillars. In addition to theoretical results, we show the experimental results obtained by fabricating several MSs of the second class. Our theoretical and experimental results confirm the unique features of this approach for knowledge discovery in EM nanostructures.


title: ReduKtor: How We Stopped Worrying About Bugs in Kotlin Compiler

abstract: Bug localization is well-known to be a difficult problem in software engineering, and specifically in compiler development, where it is beneficial to reduce the input program to a minimal reproducing example; this technique is more commonly known as delta debugging. What additionally contributes to the problem is that every new programming language has its own unique quirks and foibles, making it near impossible to reuse existing tools and approaches with full efficiency. In this experience paper we tackle the delta debugging problem w.r.t. Kotlin, a relatively new programming language from JetBrains. Our approach is based on a novel combination of program slicing, hierarchical delta debugging and Kotlin-specific transformations, which are synergistic to each other. We implemented it in a prototype called ReduKtor and did extensive evaluation on both synthetic and real Kotlin programs; we also compared its performance with classic delta debugging techniques. The evaluation results support the practical usability of our approach to Kotlin delta debugging and also shows the importance of using both language-agnostic and language-specific techniques to achieve best reduction efficiency and performance.


title: Simple model for efficient search of high-mobility organic semiconductors

abstract: High charge mobility in active layers of organic electronic devices is often necessary for their efficient operation. As a result, search for high-mobility materials among the plethora of synthesizable organic semiconductors is of paramount importance for organic electronics. However, a model for rapid but reliable prediction of charge mobility in various organic semiconductors is still lacking. To solve this issue, we propose a simple analytical model that considers the most essential factors governing the charge transport in these materials: intermolecular electronic coupling, charge delocalization, electron-phonon interaction, and static and dynamic disorder. The suggested model efficiently predicts charge mobility for organic semiconductors of various chemical structures and sizes, and significantly outperforms another approach usually used for screening of organic semiconductors, - Marcus model, - which overlooks high-mobility materials possessing small conjugated cores. We thus anticipate that the suggested model is an efficient tool for the search of the high-mobility organic semiconductors, and its possible integration with crystal structure prediction can boost the development of organic electronics.


title: Spin-orbit coupling in the hydrogen atom, the Thomas precession, and the exact solution of Dirac's equation

abstract: Bohr's model of the hydrogen atom can be extended to account for the observed spin-orbit interaction, either with the introduction of the Thomas precession, or with the stipulation that, during a spin-flip transition, the orbital radius remains intact. In other words, if there is a desire to extend Bohr's model to accommodate the spin of the electron, then experimental observations mandate the existence of the Thomas precession, which is a questionable hypothesis, or the existence of artificially robust orbits during spin-flip transitions. This is tantamount to admitting that Bohr's model, which is a poor man's way of understanding the hydrogen atom, is of limited value, and that one should really rely on Dirac's equation for the physical meaning of spin, for the mechanism that gives rise to the gyromagnetic coefficient g=2, for Zeeman splitting, for relativistic corrections to Schr\"odinger's equation, for Darwin's term, and for the correct 1/2 factor in the spin-orbit coupling energy.


title: Quantum optics with single molecules in a three-dimensional polymeric platform

abstract: The successful development of future photonic quantum technologies heavily depends on the possibility of realizing robust, reliable and, crucially, scalable nanophotonic devices. In integrated networks, quantum emitters can be deployed as single-photon sources or non-linear optical elements, provided their transition linewidth is broadened only by spontaneous emission. However, conventional fabrication approaches are hardly scalable, typically detrimental for the emitter coherence properties and bear limitations in terms of geometries and materials. Here we introduce an alternative platform, based on molecules embedded in polymeric photonic structures. Three-dimensional patterns are achieved via direct laser writing around selected molecular emitters, which preserve near-Fourier-limited fluorescence. By using an integrated polymeric design, record-high photon fluxes from a single cold molecule are reported. The proposed technology allows to conceive a novel class of quantum devices, including integrated multi-photon interferometers, arrays of indistinguishable single photon sources and hybrid electro-optical nanophotonic devices.


title: Numerically exact treatment of many body self-organization in a cavity

abstract: We investigate the full quantum evolution of ultracold interacting bosonic atoms on a chain and coupled to an optical cavity. Extending the time-dependent matrix product state techniques and the many-body adiabatic elimination technique to capture the global coupling to the cavity mode and the open nature of the cavity, we examine the long time behavior of the system beyond the mean-field elimination of the cavity field. We investigate the many body steady states and the self-organization transition for a wide range of parameters. We show that in the self-organized phase the steady state consists in a mixture of the mean-field predicted density wave states and excited states with additional defects. In particular, for large dissipation strengths a steady state with a fully mixed atomic sector is obtained crucially different from the predicted mean-field state.


title: Hyper-Differential Sensitivity Analysis of Uncertain Parameters in PDE-Constrained Optimization

abstract: Many problems in engineering and sciences require the solution of large scale optimization constrained by partial differential equations (PDEs). Though PDE-constrained optimization is itself challenging, most applications pose additional complexity, namely, uncertain parameters in the PDEs. Uncertainty quantification (UQ) is necessary to characterize, prioritize, and study the influence of these uncertain parameters. Sensitivity analysis, a classical tool in UQ, is frequently used to study the sensitivity of a model to uncertain parameters. In this article, we introduce "hyper-differential sensitivity analysis" which considers the sensitivity of the solution of a PDE-constrained optimization problem to uncertain parameters. Our approach is a goal-oriented analysis which may be viewed as a tool to complement other UQ methods in the service of decision making and robust design. We formally define hyper-differential sensitivity indices and highlight their relationship to the existing optimization and sensitivity analysis literatures. Assuming the presence of low rank structure in the parameter space, computational efficiency is achieved by leveraging a generalized singular value decomposition in conjunction with a randomized solver which converts the computational bottleneck of the algorithm into an embarrassingly parallel loop. Two multi-physics examples, consisting of nonlinear steady state control and transient linear inversion, demonstrate efficient identification of the uncertain parameters which have the greatest influence on the optimal solution.


title: Advantages of $q$-logarithm representation over $q$-exponential representation from the sense of scale and shift on nonlinear systems

abstract: Addition and subtraction of observed values can be computed under the obvious and implicit assumption that the scale unit of measurement should be the same for all arguments, which is valid even for any nonlinear systems. This paper starts with the distinction between exponential and non-exponential family in the sense of the scale unit of measurement. In the simplest nonlinear model ${dy}/{dx}=y^{q}$, it is shown how typical effects such as rescaling and shift emerge in the nonlinear systems and affect observed data. Based on the present results, the two representations, namely the $q$-exponential and the $q$-logarithm ones, are proposed. The former is for rescaling, the latter for unified understanding with a fixed scale unit. As applications of these representations, the corresponding entropy and the general probability expression for unified understanding with a fixed scale unit are presented. For the theoretical study of nonlinear systems, $q$-logarithm representation is shown to have significant advantages over $q$-exponential representation.


title: An evolving broad iron line from the first Galactic ultraluminous X-ray pulsar Swift J0243.6+6124

abstract: We present a spectral study of the ultraluminous Be/X-ray transient pulsar Swift J0243.6+6124 using Neutron Star Interior Composition Explorer (NICER) observations during the system's 2017--2018 giant outburst. The 1.2--10~keV energy spectrum of the source can be approximated with an absorbed cut-off power law model. We detect strong, luminosity-dependent emission lines in the 6--7 keV energy range. A narrow 6.42 keV line, observed in the sub-Eddington regime, is seen to evolve into a broad Fe-line profile in the super-Eddington regime. Other features are found at 6.67 and 6.97 keV in the Fe-line complex. An asymmetric broad line profile, peaking at 6.67 keV, is possibly due to Doppler effects and gravitational redshift. The 1.2--79 keV broadband spectrum from NuSTAR and NICER observations at the outburst peak is well described by an absorbed cut-off power law plus multiple Gaussian lines and a blackbody component. Physical reflection models are also tested to probe the broad iron line feature. Depending on the mass accretion rate, we found emission sites that are evolving from ~5000 km to a range closer to the surface of the neutron star. Our findings are discussed in the framework of the accretion disk and its implication on the magnetic field, the presence of optically thick accretion curtain in the magnetosphere, jet emission, and the massive, ultra-fast outflow expected at super-Eddington accretion rates. We do not detect any signatures of a cyclotron absorption line in the NICER or NuSTAR data.


title: Interactive Martingale Tests for the Global Null

abstract: Global null testing is a classical problem going back about a century to Fisher's and Stouffer's combination tests. In this work, we present simple martingale analogs of these classical tests, which are applicable in two distinct settings: (a) the online setting in which there is a possibly infinite sequence of $p$-values, and (b) the batch setting, where one uses prior knowledge to preorder the hypotheses. Through theory and simulations, we demonstrate that our martingale variants have higher power than their classical counterparts even when the preordering is only weakly informative. Finally, using a recent idea of "masking" $p$-values, we develop a novel interactive test for the global null that can take advantage of covariates and repeated user guidance to create a data-adaptive ordering that achieves higher detection power against structured alternatives.


title: Generating Synthetic Cosmological Data with GalSampler

abstract: As part of the effort to meet the needs of the Large Synoptic Survey Telescope Dark Energy Science Collaboration (LSST DESC) for accurate, realistically complex mock galaxy catalogs, we have developed GalSampler, an open-source python package that assists in generating large volumes of synthetic cosmological data. The key idea behind GalSampler is to recast hydrodynamical simulations and semi-analytic models as physically-motivated galaxy libraries. GalSampler populates a new, larger-volume halo catalog with galaxies drawn from the baseline library; by using weighted sampling guided by empirical modeling techniques, GalSampler inherits statistical accuracy from the empirical model and physically-motivated complexity from the baseline library. We have recently used GalSampler to produce the cosmoDC2 extragalactic catalog made for the LSST DESC Data Challenge 2. Using cosmoDC2 as a guiding example, we outline how GalSampler can continue to support ongoing and near-future galaxy surveys such as the Dark Energy Survey (DES), the Dark Energy Spectroscopic Instrument (DESI), WFIRST, and Euclid.


title: Entanglement front generated by an impurity travelling in an isolated many-body quantum system

abstract: We investigate the effect on the entanglement dynamics of an impurity moving at constant velocity in a closed quantum system. We focus on one-dimensional strongly-correlated lattice models, both in the presence of integrable and chaotic dynamics. In the former, the slow impurity is preceded by fast quasiparticles carrying an "endogenous" entanglement front which decays in time as a power-law; on the contrary, a fast impurity drags itself an "exogenous" entanglement front which never fades. We argue that these effects are valid for generic systems whose correlations propagate inside a light-cone. To assess the fully chaotic regime, we formulate a random circuit model which supports a moving impurity and a sharp lightcone. Although the qualitative behavior is similar to the integrable case, the endogenous regime is only visible at short times due to the onset of diffusive energy transport. Our predictions are supported by numerical simulations in the different regimes.


title: Multilingual Neural Machine Translation for Zero-Resource Languages

abstract: In recent years, Neural Machine Translation (NMT) has been shown to be more effective than phrase-based statistical methods, thus quickly becoming the state of the art in machine translation (MT). However, NMT systems are limited in translating low-resourced languages, due to the significant amount of parallel data that is required to learn useful mappings between languages. In this work, we show how the so-called multilingual NMT can help to tackle the challenges associated with low-resourced language translation. The underlying principle of multilingual NMT is to force the creation of hidden representations of words in a shared semantic space across multiple languages, thus enabling a positive parameter transfer across languages. Along this direction, we present multilingual translation experiments with three languages (English, Italian, Romanian) covering six translation directions, utilizing both recurrent neural networks and transformer (or self-attentive) neural networks. We then focus on the zero-shot translation problem, that is how to leverage multi-lingual data in order to learn translation directions that are not covered by the available training material. To this aim, we introduce our recently proposed iterative self-training method, which incrementally improves a multilingual NMT on a zero-shot direction by just relying on monolingual data. Our results on TED talks data show that multilingual NMT outperforms conventional bilingual NMT, that the transformer NMT outperforms recurrent NMT, and that zero-shot NMT outperforms conventional pivoting methods and even matches the performance of a fully-trained bilingual system.


title: Cosmological Correlators, In-In Formalism and Double Copy

abstract: We are investigating if the double copy structure as product of scattering amplitudes of gauge theories applies to cosmological correlators computed in a class of theories for inflation by the operatorial version of the In-In formalism of Schwinger-Keldysh. We consider tree level momentum-space correlators involving primordial gravitational waves with different polarizations and the scalar curvature fluctuations on a three dimensional fixed spatial slice. The correlators are sum of terms factorized in a time dependent scalar factor, which takes into account the curved background where energy is not conserved, and in a so-called tensor factor, constructed by polarization tensors. In the latter we recognize scattering amplitudes in four dimensional Minkowski space spanned by three points gravitational amplitudes related by double copy to those of gauge theories. Our study indicates that gravitational waves are double copy of gluons and the primordial scalar curvature is double copy of a scalar with Higgs-like interactions.


title: Near-term quantum algorithms for linear systems of equations

abstract: Solving linear systems of equations is essential for many problems in science and technology, including problems in machine learning. Existing quantum algorithms have demonstrated the potential for large speedups, but the required quantum resources are not immediately available on near-term quantum devices. In this work, we study near-term quantum algorithms for linear systems of equations of the form $Ax = b$. We investigate the use of variational algorithms and analyze their optimization landscapes. There exist types of linear systems for which variational algorithms designed to avoid barren plateaus, such as properly-initialized imaginary time evolution and adiabatic-inspired optimization, suffer from a different plateau problem. To circumvent this issue, we design near-term algorithms based on a core idea: the classical combination of variational quantum states (CQS). We exhibit several provable guarantees for these algorithms, supported by the representation of the linear system on a so-called Ansatz tree. The CQS approach and the Ansatz tree also admit the systematic application of heuristic approaches, including a gradient-based search. We have conducted numerical experiments solving linear systems as large as $2^{300} \times 2^{300}$ by considering cases where we can simulate the quantum algorithm efficiently on a classical computer. These experiments demonstrate the algorithms' ability to scale to system sizes within reach in near-term quantum devices of about $100$-$300$ qubits.


title: The canonical foliation on null hypersurfaces in low regularity

abstract: Let $\mathcal{H}$ denote the future outgoing null hypersurface emanating from a spacelike 2-sphere $S$ in a vacuum spacetime $(\mathcal{M},\mathbf{g})$. In this paper we study the so-called canonical foliation on $\mathcal{H}$ introduced by Klainerman and Nicol\`o and show that the corresponding geometry is controlled locally only in terms of the initial geometry on $S$ and the $L^2$ curvature flux through $\mathcal{H}$. In particular, we show that the ingoing and outgoing null expansions $\mathrm{tr} \chi$ and $\mathrm{tr} \underline{\chi}$ are both locally uniformly bounded. The proof of our estimates relies on a generalisation of the methods of Klainerman and Rodnianski, and Alexakis, Shao and Wang where the geodesic foliation on null hypersurfaces $\mathcal{H}$ is studied. The results of this paper, while of independent interest, are essential for the proof of the spacelike-characteristic bounded $L^2$ curvature theorem by Czimek and Graf.


title: Deep Learning the Morphology of Dark Matter Substructure

abstract: Strong gravitational lensing is a promising probe of the substructure of dark matter halos. Deep learning methods have the potential to accurately identify images containing substructure, and differentiate WIMP dark matter from other well motivated models, including vortex substructure of dark matter condensates and superfluids. This is crucial in future efforts to identify the true nature of dark matter. We implement, for the first time, a classification approach to identifying dark matter substructure based on simulated strong lensing images with different substructure. Utilizing convolutional neural networks trained on sets of simulated images, we demonstrate the feasibility of deep neural networks to reliably distinguish among different types of dark matter substructure. With thousands of strong lensing images anticipated with the coming launch of LSST, we expect that supervised and unsupervised deep learning models will play a crucial role in determining the nature of dark matter.


title: Inserting one edge into a simple drawing is hard

abstract: A {\em simple drawing} $D(G)$ of a graph $G$ is one where each pair of edges share at most one point: either a common endpoint or a proper crossing. An edge $e$ in the complement of $G$ can be {\em inserted} into $D(G)$ if there exists a simple drawing of $G+e$ extending $D(G)$. As a result of Levi's Enlargement Lemma, if a drawing is rectilinear (pseudolinear), that is, the edges can be extended into an arrangement of lines (pseudolines), then any edge in the complement of $G$ can be inserted. In contrast, we show that it is NP -complete to decide whether one edge can be inserted into a simple drawing. This remains true even if we assume that the drawing is pseudocircular, that is, the edges can be extended to an arrangement of pseudocircles. On the positive side, we show that, given an arrangement of pseudocircles $\mathcal{A}$ and a pseudosegment $\sigma$, it can be decided in polynomial time whether there exists a pseudocircle $\Phi_\sigma$ extending $\sigma$ for which $\mathcal{A}\cup\{\Phi_\sigma\}$ is again an arrangement of pseudocircles.


title: A gravitational collapse singularity theorem consistent with black hole evaporation

abstract: The global hyperbolicity assumption present in gravitational collapse singularity theorems is in tension with the quantum mechanical phenomenon of black hole evaporation. In this work I show that the causality conditions in Penrose's theorem can be almost completely removed. As a result, it is possible to infer the formation of spacetime singularities even in the absence of predictability and hence compatibly with quantum field theory and black hole evaporation.


title: Onset of ferromagnetism for strongly correlated electrons in one-dimensional chains

abstract: The existence of the Nagaoka ferromagnetism is examined in the context of the one-dimensional $U=\infty$ Hubbard model. We construct the exact quantum partition function to describe the physics of such a regime. Our calculation reveals that, while the ground state in an open chain is always spin-degenerate, in a finite size closed chain with at least one vacancy, the ground state can only be ferromagnetic when the number of electrons is less or equal to three. Our results shed more light on a very recent experimental verification of Nagaoka ferromagnetism in a quantum dot set up.


title: Chudnovsky-Ramanujan Type Formulae for non-Compact arithmetic triangle groups

abstract: We develop a uniform method to derive Chudnovsky-Ramanujan type formulae for triangle groups based on a generalization of a method of Chudnovsky and Chudnovsky; in particular, we carry out the method systematically for non-compact arithmetic triangle groups and one non-Fuchsian covering. As a result, we derive all rational Ramanujan type series given by Chan-Cooper for levels 1-4, as well as two additional rational series of a similar form prescribed by Chan-Cooper for these levels, but not found in the paper of Chan-Cooper. These two additional series were first found by Z.-W. Sun in a slightly different form. We also derive additional rational series of a similar form, but not found in the papers of Chan-Cooper nor Z.-W. Sun. As an ingredient in the method, we give an algorithm to rigorously confirm the singular values of normalized Eisenstein series of weight 2, which may be of independent interest.


title: Spins, percolation and height functions

abstract: To highlight certain similarities in combinatorial representations of several well known two-dimensional models of statistical mechanics, we introduce and study a new family of models which specializes to these cases after a proper tuning of the parameters. To be precise, our model consists of two independent standard Potts models, with possibly different numbers of spins and different coupling constants (the four parameters of the model), defined jointly on a graph embedded in a surface and its dual graph, and conditioned on the event that the primal and dual interfaces between spins of different value do not intersect. We also introduce naturally related height function and bond percolation models, and we discuss their basic properties and mutual relationship. As special cases we recover the standard Potts and random cluster model, the six-vertex model and loop $O(n)$ model, the random current, double random current and XOR-Ising model.


title: Audio-Visual Speech Separation and Dereverberation with a Two-Stage Multimodal Network

abstract: Background noise, interfering speech and room reverberation frequently distort target speech in real listening environments. In this study, we address joint speech separation and dereverberation, which aims to separate target speech from background noise, interfering speech and room reverberation. In order to tackle this fundamentally difficult problem, we propose a novel multimodal network that exploits both audio and visual signals. The proposed network architecture adopts a two-stage strategy, where a separation module is employed to attenuate background noise and interfering speech in the first stage and a dereverberation module to suppress room reverberation in the second stage. The two modules are first trained separately, and then integrated for joint training, which is based on a new multi-objective loss function. Our experimental results show that the proposed multimodal network yields consistently better objective intelligibility and perceptual quality than several one-stage and two-stage baselines. We find that our network achieves a 21.10% improvement in ESTOI and a 0.79 improvement in PESQ over the unprocessed mixtures. Moreover, our network architecture does not require the knowledge of the number of speakers.


abstract: The implications of sufficiently large quantum computers for widely used public-key cryptography is well-documented and increasingly discussed by the security community. An April 2016 report by the National Institute of Standards and Technology (NIST), notably, calls out the need for new standards to replace cryptosystems based on integer factorization and discrete logarithm problems, which have been shown to be vulnerable to Shor's quantum algorithm for prime factorization. Specifically, widely used RSA, ECDSA, ECDH, and DSA cryptosystems will need to be replaced by post-quantum cryptography (PQC) alternatives (also known as quantum-resistant or quantum-safe cryptography). Failure to transition before sufficiently powerful quantum computers are realized will jeopardize the security of public key cryptosystems which are widely deployed within communication protocols, digital signing mechanisms, authentication frameworks, and more. To avoid this, NIST has actively led a PQC standardization effort since 2016, leveraging a large and international research community. On January 31-February 1, 2019, the Computing Community Consortium (CCC) held a workshop in Washington, D.C. to discuss research challenges associated with PQC migration. Entitled, "Identifying Research Challenges in Post Quantum Cryptography Migration and Cryptographic Agility", participants came from three distinct yet related communities: cryptographers contributing to the NIST PQC standards effort, applied cryptographers with expertise in creating cryptographic solutions and implementing cryptography in real-world settings, and industry practitioners with expertise in deploying cryptographic standards within products and compute infrastructures. Discussion centered around two key themes: identifying constituent challenges in PQC migration and imagining a new science of "cryptographic agility".


title: The Science of Murray Gell-Mann

abstract: This article summarizes some of the most important scientific contributions of Murray Gell-Mann (1929-2019). (Invited article for Current Science, Indian Academy of Sciences.)


title: The spacelike-characteristic Cauchy problem of general relativity in low regularity

abstract: In this paper we study the spacelike-characteristic Cauchy problem for the Einstein vacuum equations. We prove that given initial data on a maximal compact spacelike hypersurface $\Sigma \simeq \overline{B(0,1)} \subset \mathbb{R}^3$ and the outgoing null hypersurface $\mathcal{H}$ emanating from $\partial \Sigma$, the time of existence of a solution to the Einstein vacuum equations is controlled by low regularity bounds on the initial data at the level of curvature in $L^2$. The proof uses the bounded $L^2$ curvature theorem by Klainerman, Szeftel and Rodnianski, the extension procedure for the constraint equations by Czimek, Cheeger-Gromov theory in low regularity developed by Czimek, the canonical foliation on null hypersurfaces in low regularity by Czimek and Graf, and global elliptic estimates for spacelike maximal hypersurfaces.


title: Spectral diagnostics of cool flare loops observed by SST: I. Inversion of the Ca II 8542 \r{A} and H$\beta$ lines

abstract: Flare loops form an integral part of eruptive events, being detected in the range of temperatures from X-rays down to cool chromospheric-like plasmas. While the hot loops are routinely observed by the Solar Dynamics Observatory's Atmospheric Imaging Assembly (SDO/AIA), cool loops seen off-limb are rare. In this paper we employ unique observations of the SOL2017-09-10T16:06 X8.2-class flare which produced an extended arcade of loops. The Swedish 1-m Solar Telescope (SST) made a series of spectral images of the cool off-limb loops in the Ca II 8542 \r{A} and the hydrogen H$\beta$ lines. Our focus is on the loop apices. Non-LTE spectral inversion is achieved through the construction of extended grids of models covering a realistic range of plasma parameters. The Multilevel Accelerated Lambda Iterations (MALI) code solves the non-LTE radiative-transfer problem in a 1D externally-illuminated slab, approximating the studied loop segment. Inversion of the Ca II 8542 \r{A} and H$\beta$ lines yields two similar solutions, both indicating high electron densities around $2 \times 10^{12}$ cm$^{-3}$ and relatively large microturbulence around 25 kms$^{-1}$. These are in reasonable agreement with other independent studies of the same or similar events. In particular, the high electron densities in the range $10^{12} - 10^{13}$ cm$^{-3}$ are consistent with those derived from the SDO's Helioseismic and Magnetic Imager white-light observations. The presence of such high densities in solar eruptive flares supports the loop interpretation of the optical continuum emission of stars which manifest superflares.


title: Casimir effect in free-fall towards a Schwarzschild black hole

abstract: In this paper we discuss the Casimir effect in a small cavity, freely falling from spatial infinity in spacetime geometry outside of a Schwarzschild black hole. Our main goal is to search for possible changes in the vacuum energy, as well as particle creation inside the falling cavity, with respect to a comoving observer. Working in the Lema\^itre chart and assuming a cavity size $L$ much smaller than the Schwarzschild radius ($L/r_g\ll 1$), we solve the Klein-Gordon equation for a massless scalar field confined within the cavity in the reference frame of the comoving observer. We follow Schwinger's proper time approach, evaluating the one-loop effective action for the field in the falling cavity hence evaluating the corrections to the vacuum energy. We find a small reduction in the absolute value of Casimir energy as the cavity approaches the black hole horizon due to the changing spacetime geometry. Since the spacetime geometry for the cavity changes dynamically, we further find the energy density of the created particles due to the dynamical Casimir effect. These dynamical contributions exactly match the deficit to the static Casimir energy. Combined, the observer measures a net increase in energy within the cavity as she falls.


title: Zero temperature momentum distribution of an impurity in a polaron state of one-dimensional Fermi and Tonks-Girardeau gases

abstract: We investigate the momentum distribution function of a single distinguishable impurity particle which formed a polaron state in a gas of either free fermions or Tonks-Girardeau bosons in one spatial dimension. We obtain a Fredholm determinant representation of the distribution function for the Bethe ansatz solvable model of an impurity-gas $\delta$-function interaction potential at zero temperature, in both repulsive and attractive regimes. We deduce from this representation the fourth power decay at a large momentum, and a weakly divergent (quasi-condensate) peak at a finite momentum. We also demonstrate that the momentum distribution function in the limiting case of infinitely strong interaction can be expressed through a correlation function of the one-dimensional impenetrable anyons.


title: Unipotent Representations Attached to the Principal Nilpotent Orbit

abstract: In this paper, we construct and classify the special unipotent representations of a real reductive group attached to the principal nilpotent orbit. We give formulas for the $\mathbf{K}$-types, associated varieties, and Langlands parameters of all such representations.


title: Subgroups of the mapping class group of the torus generated by powers of Dehn twists

abstract: We study subgroups of the mapping class group of the torus generated by powers generated by powers of Dehn twists. We give a criterion to show when a collection of powers Dehn twists generates a free group using the ping pong lemma. We show that the subgroup generated by three uniform powers of Dehn twists can be either the whole mapping class group, a direct product of a free group of rank two with the cyclic group of order two, or a free group of rank at most three. We characterize subgroups generated by a collection of (respectively squares of) Dehn twists, if there is a pair among so that the geometric intersection number of the corresponding simple closed curves is two (respectively one). We also determine the subgroup generated by uniform powers of all Dehn twists.


title: Test black holes, scattering amplitudes and perturbations of Kerr spacetime

abstract: It has been suggested that amplitudes for quantum higher-spin massive particles exchanging gravitons lead, via a classical limit, to results for scattering of spinning black holes in general relativity, when the massive particles are in a certain way minimally coupled to gravity. Such limits of such amplitudes suggest, at least at lower orders in spin, up to second order in the gravitational constant $G$, that the classical aligned-spin scattering function for an arbitrary-mass-ratio two-spinning-black-hole system can be obtained by a simple kinematical mapping from that for a spinning test black hole scattering off a stationary background Kerr black hole. Here we test these suggestions, at orders beyond the reach of the post-Newtonian and post-Minkowskian results used in their initial partial verifications, by confronting them with results from general-relativistic "self-force" calculations of the linear perturbations of a Kerr spacetime sourced by a small orbiting body, here considering only results for circular orbits in the equatorial plane. We translate between scattering and circular-orbit results by assuming the existence of a local-in-time canonical Hamiltonian governing the conservative dynamics of generic (bound and unbound) aligned-spin orbits, while employing the associated first law of spinning binary mechanics. We confirm, through linear order in the mass ratio, some previous conjectures which would begin to fill in the spin-dependent parts of the conservative dynamics for arbitrary-mass-ratio aligned-spin binary black holes at the fourth-and-a-half and fifth post-Newtonian orders.


title: Uniform asymptotics of Toeplitz determinants with Fisher-Hartwig singularities

abstract: We obtain an asymptotic formula for $n\times n$ Toeplitz determinants as $n\to \infty$, for real valued symbols with any fixed number of Fisher-Hartwig singularities, which is uniform with respect to the location of the singularities. As an application, we prove a conjecture by Fyodorov and Keating regarding moments of averages of the characteristic polynomial of the Circular Unitary Ensemble. In addition, we obtain an asymptotic formula regarding the momentum of impenetrable bosons in one dimension with periodic boundary conditions.

categories: math-ph math.CA math.CV math.FA math.MP math.PR


title: On an irreducibility type condition for the ergodicity of nonconservative semigroups

abstract: We propose a simple criterion, inspired from the irreducible aperiodic Markov chains, to derive the exponential convergence of general positive semi-groups. When not checkable on the whole state space, it can be combined to the use of Lyapunov functions. It differs from the usual generalization of irreducibility and is based on the accessibility of the trajectories of the underlying dynamics. It allows to obtain new existence results of principal eigenelements, and their exponential attractiveness, for a nonlocal selection-mutation population dynamics model defined in a space-time varying environment.


title: Thermodynamic resources in continuous-variable quantum systems

abstract: Thermodynamic resources, beyond their well-known usefulness in work extraction and other thermodynamic tasks, are often important also in tasks that are not evidently thermodynamic. Here we develop a framework for identifying such resources in diverse applications of bosonic continuous-variable systems. Introducing the class of bosonic linear thermal operations to model operationally-feasible processes, we apply this model to identify uniquely quantum properties of bosonic states that refine classical notions of thermodynamic resourcefulness. Among these are (1) a suite of temperature-like quantities generalizing the equilibrium temperature to quantum, non-equilibrium scenarios; (2) signal-to-noise ratios quantifying a system's capacity to carry information in phase-space displacement; and (3) well-established non-classicality measures quantifying the resolution in sensing and parameter estimation tasks.


title: Ramanujan graphs and exponential sums over function fields

abstract: We prove that $q+1$-regular Morgenstern Ramanujan graphs $X^{q,g}$ (depending on $g\in\mathbb{F}_q[t]$) have diameter at most $\left(\frac{4}{3}+\varepsilon\right)\log_{q}|X^{q,g}|+O_{\varepsilon}(1)$ (at least for odd $q$ and irreducible $g$) provided that a twisted Linnik-Selberg conjecture over $\mathbb{F}_q(t)$ is true. This would break the 30 year-old upper bound of $2\log_{q}|X^{q,g}|+O(1)$, a consequence of a well-known upper bound on the diameter of regular Ramanujan graphs proved by Lubotzky, Phillips, and Sarnak using the Ramanujan bound on Fourier coefficients of modular forms. We also unconditionally construct infinite families of Ramanujan graphs that prove that $\frac{4}{3}$ cannot be improved.


title: Consistency relations for large-scale structure in modified gravity and the matter bispectrum

abstract: We study perturbation theory for large-scale structure in the most general scalar-tensor theories propagating a single scalar degree of freedom, which include Horndeski theories and beyond. We model the parameter space using the effective field theory of dark energy. For Horndeski theories, the gravitational field and fluid equations are invariant under a combination of time-dependent transformations of the coordinates and fields. This symmetry allows one to construct a physical adiabatic mode which fixes the perturbation-theory kernels in the squeezed limit and ensures that the well-known consistency relations for large-scale structure, originally derived in general relativity, hold in modified gravity as well. For theories beyond Horndeski, instead, one generally cannot construct such an adiabatic mode. Because of this, the perturbation-theory kernels are modified in the squeezed limit and the consistency relations for large-scale structure do not hold. We show, however, that the modification of the squeezed limit depends only on the linear theory. We investigate the observational consequences of this violation by computing the matter bispectrum. In the squeezed limit, the largest effect is expected when considering the cross-correlation between different tracers. Moreover, the individual contributions to the 1-loop matter power spectrum do not cancel in the infrared limit of the momentum integral, modifying the power spectrum on non-linear scales.


title: Hipster random walks

abstract: We introduce and study a family of random processes on trees we call hipster random walks, special instances of which we heuristically connect to the min-plus binary trees introduced by Robin Pemantle and studied by Auffinger and Cable (2017; arXiv:1709.07849), and to the critical random hierarchical lattice studied by Hambly and Jordan (2004). We prove distributional convergence for the processes by showing that their evolutions can be understood as a discrete analogues of certain convection-diffusion equations, then using a combination of coupling arguments and results from the numerical analysis literature on convergence of numerical approximations of PDEs.


title: Document classification methods

abstract: Information on different fields which are collected by users requires appropriate management and organization to be structured in a standard way and retrieved fast and more easily. Document classification is a conventional method to separate text based on their subjects among scientific text, web pages and digital library. Different methods and techniques are proposed for document classifications that have advantages and deficiencies. In this paper, several unsupervised and supervised document classification methods are studied and compared.


title: Spatiotemporal Attention Networks for Wind Power Forecasting

abstract: Wind power is one of the most important renewable energy sources and accurate wind power forecasting is very significant for reliable and economic power system operation and control strategies. This paper proposes a novel framework with spatiotemporal attention networks (STAN) for wind power forecasting. This model captures spatial correlations among wind farms and temporal dependencies of wind power time series. First of all, we employ a multi-head self-attention mechanism to extract spatial correlations among wind farms. Then, temporal dependencies are captured by the Sequence-to-Sequence (Seq2Seq) model with a global attention mechanism. Finally, experimental results demonstrate that our model achieves better performance than other baseline approaches. Our work provides useful insights to capture non-Euclidean spatial correlations.


title: Machine learning in healthcare -- a system's perspective

abstract: A consequence of the fragmented and siloed healthcare landscape is that patient care (and data) is split along multitude of different facilities and computer systems and enabling interoperability between these systems is hard. The lack interoperability not only hinders continuity of care and burdens providers, but also hinders effective application of Machine Learning (ML) algorithms. Thus, most current ML algorithms, designed to understand patient care and facilitate clinical decision-support, are trained on limited datasets. This approach is analogous to the Newtonian paradigm of Reductionism in which a system is broken down into elementary components and a description of the whole is formed by understanding those components individually. A key limitation of the reductionist approach is that it ignores the component-component interactions and dynamics within the system which are often of prime significance in understanding the overall behaviour of complex adaptive systems (CAS). Healthcare is a CAS. Though the application of ML on health data have shown incremental improvements for clinical decision support, ML has a much a broader potential to restructure care delivery as a whole and maximize care value. However, this ML potential remains largely untapped: primarily due to functional limitations of Electronic Health Records (EHR) and the inability to see the healthcare system as a whole. This viewpoint (i) articulates the healthcare as a complex system which has a biological and an organizational perspective, (ii) motivates with examples, the need of a system's approach when addressing healthcare challenges via ML and, (iii) emphasizes to unleash EHR functionality - while duly respecting all ethical and legal concerns - to reap full benefits of ML.


title: A prototype for a serious digital game to teach linguistic ontologies

abstract: The objective of ontologies is to increase the compression of a given domain by eliminating interpretation problems. Among kinds of ontologies are linguistics ontologies which are ontologies used to simplify the interface between domain knowledge and linguistic components. Digital games have received increasing interest from educators in recent years for their potential to enhance the language learning and linguistic learning experience. Within the literature are games to teach ontologies of a specific domain, and games that use ontologies to facilitate the understanding of a given domain. Other educational games teach linguistics or vocabulary in contexts in which language is useful and meaningful. Although games help to understand difficult topics, the use of games that seek to meet the learning objectives of linguistics is not very popular and those focused on teaching linguistic ontologies are scarce. To solve the lack of the recreational resource for teaching linguistics in this document a prototype of a digital game called onto-ling is proposed. The goal is for the player to learn the relationship between concepts according to semantics, types of concepts and relationships through a game of levels.


title: Modeling Traffic Congestion with Spatiotemporal Big Data for An Intelligent Freeway Monitoring System

abstract: Traffic congestion is a complex, nonlinear spatiotemporal modeling problem. By collecting and analyzing a vast quantity and different categories of information, traffic flow, and road congestion can be predicted and controlled on an intelligent transportation system. This report provides an analysis of traveling time across Taiwan from North to South, vice versa. We analyze traffic in a national freeway between Tainan and Kaohsiung section, which represents the common trip of the population in Southern Taiwan. The data is recorded using the Electronic Toll Collection System (ETC) provided by Ministry of Transportation in Taiwan. We use MapReduce framework to process data into a smaller task which can be distributed on several computer clusters to speed up the process. The results show that the spatiotemporal model of traffic flow is strongly influenced by direction, working hour, and holidays with a recurring pattern for each week. The distinctive pattern inside the spatiotemporal dataset can be used on an AI-powered decision-making system for future development.


title: Policy Prediction Network: Model-Free Behavior Policy with Model-Based Learning in Continuous Action Space

abstract: This paper proposes a novel deep reinforcement learning architecture that was inspired by previous tree structured architectures which were only useable in discrete action spaces. Policy Prediction Network offers a way to improve sample complexity and performance on continuous control problems in exchange for extra computation at training time but at no cost in computation at rollout time. Our approach integrates a mix between model-free and model-based reinforcement learning. Policy Prediction Network is the first to introduce implicit model-based learning to Policy Gradient algorithms for continuous action space and is made possible via the empirically justified clipping scheme. Our experiments are focused on the MuJoCo environments so that they can be compared with similar work done in this area.


title: A Linearly Constrained Nonparametric Framework for Imitation Learning

abstract: In recent years, a myriad of advanced results have been reported in the community of imitation learning, ranging from parametric to non-parametric, probabilistic to non-probabilistic and Bayesian to frequentist approaches. Meanwhile, ample applications (e.g., grasping tasks and human-robot collaborations) further show the applicability of imitation learning in a wide range of domains. While numerous literature is dedicated to the learning of human skills in unconstrained environment, the problem of learning constrained motor skills, however, has not received equal attention yet. In fact, constrained skills exist widely in robotic systems. For instance, when a robot is demanded to write letters on a board, its end-effector trajectory must comply with the plane constraint from the board. In this paper, we aim to tackle the problem of imitation learning with linear constraints. Specifically, we propose to exploit the probabilistic properties of multiple demonstrations, and subsequently incorporate them into a linearly constrained optimization problem, which finally leads to a non-parametric solution. In addition, a connection between our framework and the classical model predictive control is provided. Several examples including simulated writing and locomotion tasks are presented to show the effectiveness of our framework.


title: Extending and Automating Basic Probability Theory with Propositional Computability Logic

abstract: Classical probability theory is formulated using sets. In this paper, we extend classical probability theory with propositional computability logic. Unlike other formalisms, computability logic is built on the notion of events/games, which is central to probability theory. The probability theory based on CoL is therefore useful for {\it automating} uncertainty reasoning. We describe some basic properties of this new probability theory. We also discuss a novel isomorphism between the set operations and computability logic operations.


title: Where are the Keys? -- Learning Object-Centric Navigation Policies on Semantic Maps with Graph Convolutional Networks

abstract: Emerging object-based SLAM algorithms can build a graph representation of an environment comprising nodes for robot poses and object landmarks. However, while this map will contain static objects such as furniture or appliances, many moveable objects (e.g. the car keys, the glasses, or a magazine), are not suitable as landmarks and will not be part of the map due to their non-static nature. We show that Graph Convolutional Networks can learn navigation policies to find such unmapped objects by learning to exploit the hidden probabilistic model that governs where these objects appear in the environment. The learned policies can generalise to object classes unseen during training by using word vectors that express semantic similarity as representations for object nodes in the graph. Furthermore, we show that the policies generalise to unseen environments with only minimal loss of performance. We demonstrate that pre-training the policy network with a proxy task can significantly speed up learning, improving sample efficiency.


title: A Karhunen-Loeve expansion for one-mode open quantum harmonic oscillators using the eigenbasis of the two-point commutator kernel

abstract: This paper considers one-mode open quantum harmonic oscillators with a pair of conjugate position and momentum variables driven by vacuum bosonic fields according to a linear quantum stochastic differential equation. Such systems model cavity resonators in quantum optical experiments. Assuming that the quadratic Hamiltonian of the oscillator is specified by a positive definite energy matrix, we consider a modified version of the quantum Karhunen-Loeve expansion of the system variables proposed recently. The expansion employs eigenvalues and eigenfunctions of the two-point commutator kernel for linearly transformed system variables. We take advantage of the specific structure of this eigenbasis in the one-mode case (including its connection with the classical Ornstein-Uhlenbeck process). These results are applied to computing quadratic-exponential cost functionals which provide robust performance criteria for risk-sensitive control of open quantum systems.

categories: quant-ph cs.SY eess.SY math-ph math.MP math.OC math.PR


title: Searching for Accurate Binary Neural Architectures

abstract: Binary neural networks have attracted tremendous attention due to the efficiency for deploying them on mobile devices. Since the weak expression ability of binary weights and features, their accuracy is usually much lower than that of full-precision (i.e. 32-bit) models. Here we present a new frame work for automatically searching for compact but accurate binary neural networks. In practice, number of channels in each layer will be encoded into the search space and optimized using the evolutionary algorithm. Experiments conducted on benchmark datasets and neural architectures demonstrate that our searched binary networks can achieve the performance of full-precision models with acceptable increments on model sizes and calculations.


title: On a sufficient condition for infinite horizon optimal control problems

abstract: New form of sufficient optimality condition is obtained in comparison with the Mangasarian sufficiency theorem. Both finite and infinite values of objective functional are allowed since concepts of overtaking and weakly overtaking optimality are implied. Examples, where new conditions can be applied are presented. The conditions are shown to be both necessary and sufficient, when Hamiltonian is linear with respect to state and control.


title: How far can we push deconvolution? A SCUBA-2 test case

abstract: How far can we use multi-wavelength cross-identifications to deconvolve far-infrared images? In this short research note I explore a test case of CLEAN deconvolutions of simulated confused 850 micron SCUBA-2 data, and explore the possible scientific applications of combining this data with ostensibly deeper TolTEC Large Scale Structure (LSS) survey 1.1mm-2mm data. I show that the SCUBA-2 can be reconstructed to the 1.1mm LMT resolution and achieve an 850 micron deconvolved sensitivity of 0.7 mJy RMS, an improvement of at least ~1:5x over naive point source filtered images. The TolTEC/SCUBA-2 combination can constrain cold (<10K) observed-frame colour temperatures, where TolTEC alone cannot.


title: Operator content of entanglement spectra after global quenches in the transverse field Ising chain

abstract: We consider the time evolution of the gaps of the entanglement spectrum for a block of consecutive sites in finite transverse field Ising chains after sudden quenches of the magnetic field. We provide numerical evidence that, whenever we quench at or across the quantum critical point, the time evolution of the ratios of these gaps allows to obtain universal information. They encode the low-lying gaps of the conformal spectrum of the Ising boundary conformal field theory describing the spatial bipartition within the imaginary time path integral approach to global quenches at the quantum critical point.

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


title: Activity-controlled Annealing of Colloidal Monolayers

abstract: Molecular motors are essential to the living, they generate additional fluctuations that boost transport and assist assembly. Self-propelled colloids, that consume energy to move, hold similar potential for the man-made assembly of microparticles. Yet, experiments showing their use as a powerhouse in materials science lack. Our work explores the design of man-made materials controlled by fluctuations, arising from the internal forces generated by active colloids. Here we show a massive acceleration of the annealing of a monolayer of passive beads by moderate addition of self-propelled microparticles. We rationalize our observations with a model of collisions that drive active fluctuations to overcome kinetic barriers and activate the annealing. The experiment is quantitatively compared with Brownian dynamic simulations that further unveil a dynamical transition in the mechanism of annealing. Active dopants travel uniformly in the system or co-localize at the grain boundaries as a result of the persistence of their motion. Our findings uncover the potential of man-made materials controlled by internal activity and lay the groundwork for the rise of materials science beyond equilibrium.


title: Lorentz Symmetry Fractionalization and Dualities in (2+1)d

abstract: We discuss symmetry fractionalization of the Lorentz group in (2+1)$d$ non-spin quantum field theory (QFT), and its implications for dualities. We prove that two inequivalent non-spin QFTs are dual as spin QFTs if and only if they are related by a Lorentz symmetry fractionalization with respect to an anomalous $\mathbb{Z}_2$ one-form symmetry. Moreover, if the framing anomalies of two non-spin QFTs differ by a multiple of 8, then they are dual as spin QFTs if and only if they are also dual as non-spin QFTs. Applications to summing over the spin structures, time-reversal symmetry, and level/rank dualities are explored. The Lorentz symmetry fractionalization naturally arises in Chern-Simons matter dualities that obey certain spin/charge relations, and is instrumental for the dualities to hold when viewed as non-spin theories.


title: Gravitational Waves from Supernova Mass Loss and Natal Kicks in Close Binaries

abstract: Some fraction of compact binaries that merge within a Hubble time may have formed from two massive stars in isolation. For this isolated-binary formation channel, binaries need to survive two supernova (SN) explosions in addition to surviving common-envelope evolution. For the SN explosions, both the mass loss and natal kicks change the orbital characteristics, producing either a bound or unbound binary. We show that gravitational waves (GWs) may be produced not only from the core-collapse SN process, but also from the SN mass loss and SN natal kick during the pre-SN to post-SN binary transition. We model the dynamical evolution of a binary at the time of the second SN explosion with an equation of motion that accounts for the finite timescales of the SN mass loss and the SN natal kick. From the dynamical evolution of the binary, we calculate the GW burst signals associated with the SN natal kicks. We find that such GW bursts may be of interest to future mid-band GW detectors like DECIGO. We also find that the energy radiated away from the GWs emitted due to the SN mass loss and natal kick may be a significant fraction, ${\gtrsim}10\%$, of the post-SN binary's orbital energy. For unbound post-SN binaries, the energy radiated away in GWs tends to be higher than that of bound binaries.


title: Synthetic mean-field interactions in photonic lattices

abstract: Photonic lattices are usually considered to be limited by their lack of methods to include interactions. We address this issue by introducing mean-field interactions through optical components which are external to the photonic lattice. The proposed technique to realise mean-field interacting photonic lattices relies on a Suzuki-Trotter decomposition of the unitary evolution for the full Hamiltonian. The technique realises the dynamics in an analogous way to that of a step-wise numerical implementation of quantum dynamics, in the spirit of digital quantum simulation. It is a very versatile technique which allows for the emulation of interactions that do not only depend on inter-particle separations or do not decay with particle separation. We detail the proposed experimental scheme and consider two examples of interacting phenomena, self-trapping and the decay of Bloch oscillations, that are observable with the proposed technique.


title: Constraints on the production and escape of ionizing radiation from the emission-line spectra of metal-poor star-forming galaxies

abstract: We explore the production and escape of ionizing photons in young galaxies by investigating the ultraviolet and optical emission-line properties of models of ionization-bounded and density-bounded HII regions, active-galactic-nucleus (AGN) narrow-line regions and radiative shocks computed all using the same physically-consistent description of element abundances and depletion on to dust grains down to very low metallicities. We compare these models with a reference sample of metal-poor star-forming galaxies and Lyman-continuum (LyC) leakers at various redshifts, which allows the simultaneous exploration of more spectral diagnostics than typically available at once for individual subsamples. We confirm that current single- and binary-star population synthesis models do not produce hard-enough radiation to account for the high-ionization emission of the most metal-poor galaxies. Introducing either an AGN or radiative-shock component brings models into agreement with observations. A published model including X-ray binaries is an attractive alternative to reproduce the observed rise in HeII4686/Hbeta ratio with decreasing oxygen abundance in metal-poor star-forming galaxies, but not the high observed HeII4686/Hbeta ratios of galaxies with large EW(Hbeta). A source of harder ionizing radiation appears to be required in these extreme objects, such as an AGN or radiative-shock component, perhaps linked to an initial-mass-function bias toward massive stars at low metallicity. This would also account for the surprisingly high [OI]/[OIII] ratios of confirmed LyC leakers relative to ionization-bounded models. We find no simple by-eye diagnostic of the nature of ionizing sources and the escape of LyC photon, which require proper simultaneous fits of several lines to be discriminated against.


title: Superradiant Searches for Dark Photons in Two Stage Atomic Transitions

abstract: We study a new mechanism to discover dark photon fields, by resonantly triggering two photon transitions in cold gas preparations. Using coherently prepared cold parahydrogen, coupling sensitivity for sub-meV mass dark photon fields can be advanced by orders of magnitude, with a modified light-shining-through-wall setup. We calculate the effect of a background dark photon field on the dipole moment and corresponding transition rate of cold parahydrogen pumped into its first vibrational excited state by counter-propagating laser beams. The nonlinear amplification of two photon emission triggered by dark photons in a cold parahydrogen sample is numerically simulated to obtain the expected dark photon coupling sensitivity.


title: The assembly of dusty galaxies at $z \geq 4$: statistical properties

abstract: The recent discovery of high redshift dusty galaxies implies a rapid dust enrichment of their interstellar medium (ISM). To interpret these observations, we run a cosmological simulation in a 30$h^{-1}$ cMpc/size volume down to $z \approx 4$. We use the hydrodynamical code dustyGadget, which accounts for the production of dust by stellar populations and its evolution in the ISM. We find that the cosmic dust density parameter ($\Omega_{\rm d}$) is mainly driven by stellar dust at $z \gtrsim 10$, so that mass- and metallicity-dependent yields are required to assess the dust content in the first galaxies. At $z \lesssim 9$ the growth of grains in the ISM of evolved systems (Log$(M_{\star}/M_{\odot})>8.5$) significantly increases their dust mass, in agreement with observations in the redshift range $4 \lesssim z < 8$. Our simulation shows that the variety of high redshift galaxies observed with ALMA can naturally be accounted for by modeling the grain-growth timescale as a function of the physical conditions in the gas cold phase. In addition, the trends of dust-to-metal (DTM) and dust-to-gas (${\cal D}$) ratios are compatible with the available data. A qualitative investigation of the inhomogeneous dust distribution in a representative massive halo at $z \approx 4$ shows that dust is found from the central galaxy up to the closest satellites along polluted filaments with $\rm Log({\cal D}) \leq -2.4$, but sharply declines at distances $d \gtrsim 30$ kpc along many lines of sight, where $\rm Log({\cal D}) \lesssim -4.0$.


title: Wide-Field Survey of Dwarf Satellite Systems Around 10 Hosts in the Local Volume

abstract: We present the results of an extensive search for dwarf satellite galaxies around 10 primary host galaxies in the Local Volume (D$<$12 Mpc) using archival CFHT/MegaCam imaging data. The hosts span a wide range in properties, with stellar masses ranging from that of the LMC to ${\sim}3$ times that of the Milky Way (MW). The surveyed hosts are: NGC 1023, NGC 1156, NGC 2903, NGC 4258, NGC 4565, NGC 4631, NGC 5023, M51, M64, and M104. We detect satellite candidates using a consistent semi-automated detection algorithm that is optimized for the detection of low surface brightness objects. Depending on the host, our completeness limit is $M_g{\sim}-8$ to $-10$ (assuming the distance of the host). We detect objects with surface brightness down to $\mu_{0,g}{\sim}26$ mag arcsec$^{-2}$ at $\gtrsim90\%$ completeness. The survey areas of the six best-surveyed hosts cover most of the inner projected $R<150$ kpc area, which roughly doubles the number of MW-mass hosts surveyed at this level of area and luminosity completeness. The number of detected candidates range from 1 around M64 to 33 around NGC 4258. In total, 153 candidates are found, of which 93 are new. While we defer an analysis of the satellite luminosity functions of the hosts until distance information is available for the candidates, we do show that the candidates are primarily red, spheroid systems with properties roughly consistent with known satellites in the Local Group.


title: Discovering the $h\to Z \gamma$ Decay in $t \bar t$ Associated Production

abstract: We explore the prospects to discover the $h \to Z \gamma$ decay in $t\bar t$-associated production, featuring a signal-to-background ratio of ${\cal O}(1)$. Performing a detailed analysis of the semi-leptonic $t \bar t $-decay channel, we demonstrate that the production mode could lead to a $\sim5\,\sigma$ discovery at the high-luminosity LHC, while the effective $h Z \gamma$ coupling could be extracted with a $\sim15\,\%$ accuracy. Extending the analysis to potential future $pp$ colliders with 27 TeV and 100 TeV center-of-mass energies, we also show that the latter would allow precision measurements at the few percent level, rendering possible precise extractions of the spin and CP properties of the Higgs boson.


title: Dictionary for the type II nongeometric flux compactifications

abstract: We study the $T$-dual completion of the four-dimensional ${\cal N}=1$ type II effective potentials in the presence of (non-)geometric fluxes. First, we invoke a cohomology version of the $T$-dual transformations among the various moduli, axions and the fluxes appearing in the type IIA and type IIB effective supergravities. This leads to some useful observations about a significant mixing of the standard NS-NS fluxes with the (non-)geometric fluxes on the mirror side. Further, using our $T$-duality rules, we establish an explicit mapping among the $F$-terms, $D$-terms, tadpole conditions as well as the Bianchi identities of the two theories. Secondly, we propose what we call a set of "axionic flux polynomials", which depend on all the axionic moduli and the fluxes. This subsequently helps in presenting the two scalar potentials in a concise and manifestly $T$-dual form, which can be directly utilized for various phenomenological purposes as we illustrate in a couple of examples.


title: Beyond the exoplanet mass-radius relation

abstract: Mass and radius are two fundamental properties for characterising exoplanets, but only for a relatively small fraction of exoplanets are they both available. Mass is often derived from radial velocity measurements, while the radius is almost always measured using the transit method. For a large number of exoplanets, either the radius or the mass is unknown, while the host star has been characterised. Several mass-radius relations that are dependent on the planet's type have been published that often allow us to predict the radius. Our goal is to derive the radius of exoplanets using only observables extracted from spectra used primarily to determine radial velocities and spectral parameters. Our objective is to obtain a mass-radius relation independent of the planet's type. We worked with a database of confirmed exoplanets with known radii and masses, as well as the planets from our Solar System. Using random forests, a machine learning algorithm, we computed the radius of exoplanets and compared the results to the published radii. Our code, Bem, is available online. The estimated radii reproduces the spread in radius found for high mass planets better than previous mass-radius relations. The average radius error is 1.8 $R_{\oplus}$ across the whole range of radii from 1-22 $R_{\oplus}$. We find that a random forest algorithm is able to derive reliable radii, especially for planets between 4 $R_{\oplus}$ and 20 $R_{\oplus}$ for which the error is under 25$\%$. The algorithm has a low bias yet a high variance, which could be reduced by limiting the growth of the forest, or adding more data. The random forest algorithm is a promising method for deriving exoplanet properties. We show that the exoplanet's mass and equilibrium temperature are the relevant properties that constrain the radius, and do so with higher accuracy than the previous methods.


title: Utilizing Shell Galaxies

abstract: Stellar shells are low surface brightness features in the form of open, concentric arcs, formed in close-to-radial collisions of galaxies. They occur in at least 10% of early-type galaxies and a small portion of spirals and their unique kinematics carry valuable information about the host galaxies. We discuss a method using measurements of the number and distribution of shells to estimate the mass distribution of the galaxies and the time since the merger. The method is applied on the shells of NGC 4993 - a galaxy hosting the electromagnetic counterpart of the gravitational wave event GW170817, to estimate the probable time since the galactic merger. We used analytical calculations and particle simulations to show that, in special cases, when kinematic data are available, further constraints on mass distribution and merger time can be derived. Applying the methods to the rapidly growing sample of known shell galaxies will constrain the dark-matter content in the galaxies and reveal detailed information on the recent merger history of the Universe.


title: An Imprint of the Galactic Magnetic Field in the Diffuse Unpolarized Dust Emission

abstract: It is well known that aligned, aspherical dust grains emit polarized radiation and that the degree of polarization depends on the angle $\psi$ between the interstellar magnetic field and the line of sight. However, anisotropy of the dust absorption cross sections also modulates the $total\ intensity$ of the radiation as the viewing geometry changes. We report a detection of this effect in the high Galactic latitude $Planck$ data, finding that the 353 GHz dust intensity per $N_{\rm HI}$ is smaller when the Galactic magnetic field is mostly in the plane of the sky and larger when the field is mostly along the line of sight. These variations are of opposite sign and roughly equal magnitude as the changes in polarized intensity per $N_{\rm HI}$ with $\psi$, as predicted. In principle, the variation in intensity can be used in conjunction with the dust polarization angle to constrain the full 3D orientation of the Galactic magnetic field.


title: Cusp-to-core transition in low-mass dwarf galaxies induced by dynamical heating of cold dark matter by primordial black holes

abstract: We performed a series of high-resolution $N$-body simulations to examine whether dark matter candidates in the form of primordial black holes (PBHs) can solve the cusp-core problem in low-mass dwarf galaxies. If some fraction of the dark matter in low-mass dwarf galaxies consists of PBHs and the rest is cold dark matter, dynamical heating of the cold dark matter by the PBHs induces a cusp-to-core transition in the total dark matter profile. The mechanism works for PBHs in the 25-100 M$_{\sun}$ mass window, consistent with the LIGO detections, but requires a lower limit on the PBH mass fraction of 1$\%$ of the total dwarf galaxy dark matter content. The cusp-to-core transition time-scale is between 1 and 8 Gyr. This time-scale is also a constant multiple of the relaxation time between cold dark matter particles and PBHs, which depends on the mass, the mass fraction and the scale radius of the initial density profile of PBHs. We conclude that dark matter cores occur naturally in halos comprised of cold dark matter and PBHs, without the need to invoke baryonic processes.


title: The Lyman Continuum Escape Survey -- II: Ionizing Radiation as a Function of the [OIII]/[OII] Line Ratio

abstract: We discuss the rest-frame optical emission line spectra of a large (~50) sample of z=3.1 Lyman alpha emitters (LAEs) whose physical properties suggest such sources are promising analogs of galaxies in the reionization era. Reliable Lyman continuum escape fractions have now been determined for a large sample of such LAEs from the Lyman Continuum Escape Survey (LACES) undertaken via deep HST imaging in the SSA22 survey area reported in Fletcher et al. (2019). Using new measures of [OII] emission secured from Keck MOSFIRE spectra we re-examine, for a larger sample, earlier claims that Lyman continuum leakages may correlate with the nebular emission line ratio [OIII]/[OII] as expected for density-bound HII regions. We find that a large [OIII]/[OII] line ratio is indeed a necessary condition for Lyman continuum leakage, strengthening earlier claims made using smaller samples at various redshifts. However, not all LAEs with large [OIII]/[OII] line ratios are leakers and leaking radiation appears not to be associated with differences in other spectral diagnostics. This suggests the detection of leaking radiation is modulated by an additional property, possibly the viewing angle for porous HII regions. We discuss our new results in the context of the striking bimodality of LAE leakers and non-leakers found in the LACES program and the implications for the sources of cosmic reionization.


title: Multiple bursts of star formation in young star clusters: The case of the Orion Nebula Cluster

abstract: Young star clusters (YSCs) with resolved stellar populations are well suited for studying star-cluster formation. In most cases, the (pre-main-sequence) stellar populations found in the YSCs are coeval with an intrinsic age spread of up to 1Myr. Such observations can be understood as the YSCs having formed in one burst, which star formation was truncated by stellar feedback. The recent discovery that the colour-magnitude diagram of the Orion Nebula Clusters (ONC) contains three well defined age-separated populations appears to shatter this model. The implication is that the ONC formed in three bursts, with star formation still on-going in the last burst. We present new observational results focusing on the three populations in the ONC using OmegaCAM photometry and Gaia DR2 measurements. We also describe a theoretical model which may explain these observations by an interplay between stellar feedback and cluster dynamics.


title: Explaining GRB prompt emission with sub-photospheric dissipation and Comptonization

abstract: Even though the observed spectra for GRB prompt emission is well constrained, no single radiation mechanism can robustly explain its distinct non-thermal nature. Here we explore the radiation mechanism with the photospheric emission model using our Monte Carlo Radiative Transfer (MCRaT) code. We study the sub-photospheric Comptonization of fast cooled synchrotron photons while the Maxwellian electrons and mono-energetic protons are accelerated to relativistic energies by repeated dissipation events. Unlike previous simulations, we implement a realistic photon to electron number ratio $N_{\gamma}/N_{e} \sim 10^5$ consistent with the observed radiative efficiency of a few percent. We show that it is necessary to have a critical number of episodic energy injection events $N_{rh,cr} \sim {\rm few}\ 10{\rm s}-100$ in the jet in addition to the electron-proton Coulomb coupling in order to inject sufficient energy $E_{inj,cr} \sim 2500-4000\ m_e c^2$ per electron and produce an output photon spectrum consistent with observations. The observed GRB spectrum can be generated when the electrons are repeatedly accelerated to highly relativistic energies $\gamma_{e,in} \sim {\rm few}\ 10{\rm s}-100$ in a jet with bulk Lorentz factor $\Gamma \sim 30-100$, starting out from moderate optical depths $\tau_{in} \sim 20-40$. The shape of the photon spectrum is independent of the initial photon energy distribution and baryonic energy content of the jet and hence independent of the emission mechanism, as expected for photospheric emission.


title: Confined phases of one-dimensional spinless fermions coupled to $Z_2$ gauge theory

abstract: We investigate a quantum many-body lattice system of one-dimensional spinless fermions interacting with a dynamical $Z_2$ gauge field. The gauge field mediates long-range attraction between fermions resulting in their confinement into bosonic dimers. At strong coupling we develop an exactly solvable effective theory of such dimers with emergent constraints. Even at generic coupling and fermion density, the model can be rewritten as a local spin chain. Using the Density Matrix Renormalization Group the system is shown to form a Luttinger liquid, indicating the emergence of fractionalized excitations despite the confinement of lattice fermions. In a finite chain we observe the doubling of the period of Friedel oscillations which paves the way towards experimental detection of confinement in this system. We discuss the possibility of a Mott phase at the commensurate filling $2/3$.

categories: cond-mat.str-el cond-mat.quant-gas hep-lat hep-th


title: Towards sub-kpc scale kinematics of molecular and ionized gas of star-forming galaxies at $z\sim1$

abstract: We compare the molecular and ionized gas kinematics of two strongly lensed galaxies at $z\sim1$ based on observations from ALMA and MUSE. We derive the CO and [OII] rotation curves and dispersion profiles of these two galaxies. We find a difference between the observed molecular and ionized gas rotation curves for one of the galaxies, the Cosmic Snake, for which we obtain a spatial resolution of few hundred parsecs along the major axis. The rotation curve of the molecular gas is steeper than the rotation curve of the ionized gas. In the second galaxy, A521, the molecular and ionized gas rotation curves are consistent, but the spatial resolution is only of few kpc on the major axis. Using simulations, we investigate the effect of the thickness of the gas disk and effective radius on the observed rotation curves and find that a more extended and thicker disk smooths the curve. We also find that the presence of a strongly inclined thick disk (>1 kpc) can smooth the rotation curve because it degrades the spatial resolution along the line of sight. By building a model using a stellar disk and two gas disks, we reproduce the rotation curves of the Cosmic Snake with a molecular gas disk that is more massive and more radially and vertically concentrated than the ionized gas disk. Finally, we also obtain an intrinsic velocity dispersion in the Cosmic Snake of 18.5+-7 km/s and 19.5+-6 km/s for the molecular and ionized gas, respectively, which is consistent with a molecular disk with a smaller and thinner disk. For A521, the intrinsic velocity dispersion values are 11+-8 km/s and 54+-11 km/s, with a higher value for the ionized gas. This could indicate that the ionized gas disk is thicker and more turbulent in this galaxy. These results highlight the different spatial distribution of the molecular and ionized gas disks at $z\sim1$ and suggest the presence of thick ionized gas disks at this epoch.