title: Impact of thermal fluctuations on transport in antiferromagnetic semimetals

abstract: Recent demonstrations on manipulating antiferromagnetic (AF) order have triggered a growing interest in antiferromagnetic metal (AFM), and potential high-density spintronic applications demand further improvements in the anisotropic magnetoresistance (AMR). The antiferromagnetic semimetals (AFS) are newly discovered materials that possess massless Dirac fermions that are protected by the crystalline symmetries. In this material, a reorientation of the AF order may break the underlying symmetries and induce a finite energy gap. As such, the possible phase transition from the semimetallic to insulating phase gives us a choice for a wide range of resistance ensuring a large AMR. To further understand the robustness of the phase transition, we study thermal fluctuations of the AF order in AFS at a finite temperature. For macroscopic samples, we find that the thermal fluctuations effectively decrease the magnitude of the AF order by renormalizing the effective Hamiltonian. Our finding suggests that the insulating phase exhibits a gap narrowing at elevated temperatures, which leads to a substantial decrease in AMR. We also examine spatially correlated thermal fluctuations for microscopic samples by solving the microscopic Landau-Lifshitz-Gilbert equation finding a qualitative difference of the gap narrowing in the insulating phase. For both cases, the semimetallic phase shows a minimal change in its transmission spectrum illustrating the robustness of the symmetry protected states in AFS. Our finding may serve as a guideline for estimating and maximizing AMR of the AFS samples at elevated temperatures.


title: Berry curvature and Hall viscosities in an anisotropic Dirac semi-metal

abstract: We investigate parity-odd non-dissipative transport in an anisotropic Dirac semi-metal in two spatial dimensions. The analysis is relevant for interacting electronic systems with merging Dirac points at charge neutrality. For such systems the dispersion relation is relativistic in one direction and non-relativistic in the other. We give a proposal how to calculate the Berry curvature for this system and use it to derive more than one odd viscosities, in contrast to rotationally invariant systems. We observe that in such a model the odd part of stress tensor is parameterised by two independent transport coefficients and one that is identically zero.


title: Introduction to the Dicke model: from equilibrium to nonequilibrium, and vice versa

abstract: The Dicke model describes the coupling between a quantized cavity field and a large ensemble of two-level atoms. When the number of atoms tends to infinity, this model can undergo a transition to a superradiant phase, belonging to the mean-field Ising universality class. The superradiant transition was first predicted for atoms in thermal equilibrium and was recently realized with a quantum simulator made of atoms in an optical cavity, subject to both dissipation and driving. In this Progress Report, we offer an introduction to some theoretical concepts relevant to the Dicke model, reviewing the critical properties of the superradiant phase transition, and the distinction between equilibrium and nonequilibrium conditions. In addition, we explain the fundamental difference between the superradiant phase transition and the more common lasing transition. Our report mostly focuses on the steady states of atoms in single-mode optical cavities, but we also mention some aspects of real-time dynamics, as well as other quantum simulators, including superconducting qubits, trapped ions, and using spin-orbit coupling for cold atoms. These realizations differ in regard to whether they describe equilibrium or non-equilibrium systems.


title: Redshifted broad absorption line quasars found via machine-learned spectral similarity

abstract: We report the discovery of 31 new redshifted broad absorption line quasars (RSBALs) from the Sloan Digital Sky Survey (SDSS). The number of previously known such objects is 19. The identification of the new objects was enabled by calculating similarities between quasar spectra in the SDSS. Using these similarities we look for the objects that are similar to the ones in the original sample, visually inspecting only hundreds, out of over 160,000 spectra considered. We compare the performance of several similarity measures, as well as different methods of employing them, in finding the RSBALs. We find that decision tree based similarities recover the most objects, and that an ensemble of methods performs better than any single one. As the similarities are not tailored for the specific problem of finding RSBALs, they could be used for searching for other types of quasars. The similarities and the code for their calculation are available online.


title: Implications of a Dark Matter-Neutrino Coupling at Hyper-Kamiokande

abstract: Dark matter and neutrinos provide the two most compelling pieces of evidence for new physics beyond the Standard Model of Particle Physics but they are often treated as two different sectors. In this paper, we consider how neutrino observables can be used to constrain the parameter space of different models where active neutrinos interact with DM via mediators of different spins. We study for the first time the sensitivity of the Hyper-Kamiokande detector to neutrinos produced from MeV Dark Matter annihilation. In particular, we find that thermally produced DM candidates with masses between 15-30 MeV could be fully excluded.


title: Classical phase diagram of the stuffed honeycomb lattice

abstract: We investigate the classical phase diagram of the stuffed honeycomb Heisenberg lattice, which consists of a honeycomb lattice with a superimposed triangular lattice formed by sites at the center of each hexagon. This lattice encompasses and interpolates between the honeycomb, triangular and dice lattices, preserving the hexagonal symmetry while expanding the phase space for potential spin liquids. We use a combination of iterative minimization, classical Monte Carlo and analytical techniques to determine the complete ground state phase diagram. It is quite rich, with a variety of non-coplanar and non-collinear phases not found in the previously studied limits. In particular, our analysis reveals the triangular lattice critical point to be a multicritical point with two new phases vanishing via second order transitions at the critical point. We analyze these phases within linear spin wave theory and discuss consequences for the S = 1/2 spin liquid.


title: $\mathcal{N}=(3,3)$ holography on ${\rm AdS}_3 \times ({\rm S}^3 \times {\rm S}^3 \times {\rm S}^1)/\mathbb Z_2$

abstract: We consider string theory on AdS$_3$ $\times$ (S$^3$ $\times$ S$^3$ $\times$ S$^1)/\mathbb Z_2$, a background supporting $\mathcal N=(3,3)$ spacetime supersymmetry. We propose that string theory on this background is dual to the symmetric product orbifold of $\mathcal S_0/\mathbb Z_2$ where $\mathcal S_0$ is a theory of four free fermions and one free boson. We show that the BPS spectra of the two sides of the duality match precisely. Furthermore, we compute the elliptic genus of the dual CFT and that of the supergravity limit of string theory and demonstrate that they match, hence providing non-trivial support for the holographic proposal.


title: The Six Cylinders Problem: $\mathbb{D}_{3}$-symmetry Approach

abstract: Motivated by a question of W. Kuperberg, we study the 18-dimensional manifold of configurations of 6 non-intersecting infinite cylinders of radius $r,$ all touching the unit ball in $\mathbb{R}^{3}.$ We find a configuration with \[ r=\frac{1}{8}\left( 3+\sqrt{33}\right) \approx1.093070331\ .\] We believe that this value is the maximal possible.


title: RELICS: A Strong Lens Model for SPT-CLJ0615-5746, a z=0.972 Cluster

abstract: We present a lens model for the cluster SPT-CLJ0615$-$5746, which is the highest redshift ($z=0.972$) system in the Reionization of Lensing Clusters Survey (RELICS), making it the highest redshift cluster for which a full strong lens model is published. We identify three systems of multiply-imaged lensed galaxies, two of which we spectroscopically confirm at $z=1.358$ and $z=4.013$, which we use as constraints for the model. We find a foreground structure at $z\sim0.4$, which we include as a second cluster-sized halo in one of our models; however two different statistical tests find the best-fit model consists of one cluster-sized halo combined with three individually optimized galaxy-sized halos, as well as contributions from the cluster galaxies themselves. We find the total projected mass density within $r=26.7"$ (the region where the strong lensing constraints exist) to be $M=2.51^{+0.15}_{-0.09}\times 10^{14}$~M$_{\odot}$. If we extrapolate out to $r_{500}$, our projected mass density is consistent with the mass inferred from weak lensing and from the Sunyaev-Zel'dovich effect ($M\sim10^{15}$~M$_{\odot}$). This cluster is lensing a previously reported $z\sim10$ galaxy, which, if spectroscopically confirmed, will be the highest-redshift strongly lensed galaxy known.


title: A Family of Pb-based Superconductors with Variable Cubic to Hexagonal Packing

abstract: We describe three previously unreported superconductors, BaPb3, Ba0.89Sr0.11Pb3 and Ba0.5Sr0.5Pb3. These three materials, together with SrPb3, form a distinctive isoelectronic family of intermetallic superconductors based on the stacking of Pb planes, with crystal structures that display a hexagonal to cubic perovskite-like progression, as rarely seen in metals. The superconducting transition temperatures (Tc) are similar for all - 2.2 K for BaPb3, 2.7 K for Ba0.89Sr0.11Pb3 and 2.6 K for Ba0.5Sr0.5Pb3, and the previously reported Tc of SrPb3, ~ 2 K, is confirmed. The materials are moderate coupling superconductors, and calculations show that the electronic densities of states at the Fermi energy are primarily contributed by Pb. The observations suggest that the Pb-stacking variation has only a minor effect on the superconductivity.


title: Finiteness spaces and generalized power series

abstract: We consider Ribenboim's construction of rings of generalized power series. Ribenboim's construction makes use of a special class of partially ordered monoids and a special class of their subsets. While the restrictions he imposes might seem conceptually unclear, we demonstrate that they are precisely the appropriate conditions to represent such monoids as internal monoids in an appropriate category of Ehrhard's finiteness spaces. Ehrhard introduced finiteness spaces as the objects of a categorical model of classical linear logic, where a set is equipped with a class of subsets to be thought of as finitary. Morphisms are relations preserving the finitary structure. The notion of finitary subset allows for a sharper analysis of computational structure than is available in the relational model. For example, fixed point operators fail to be finitary. In the present work, we take morphisms to be partial functions preserving the finitary structure rather than relations. The resulting category is symmetric monoidal closed, complete and cocomplete. Any pair of an internal monoid in this category and a ring induces a ring of generalized power series by an extension of the Ribenboim construction based on Ehrhard's notion of linearization of a finiteness space. We thus further generalize Ribenboim's constructions. We give several examples of rings which arise from this construction, including the ring of Puiseux series and the ring of formal power series generated by a free monoid.


title: Solitary Waves of the Two-Dimensional Camassa-Holm--Nonlinear Schr\"odinger Equation

abstract: In this work, we study solitary waves in a (2+1)-dimensional variant of the defocusing nonlinear Schr\"odinger (NLS) equation, the so-called Camassa-Holm NLS (CH-NLS) equation. We use asymptotic multiscale expansion methods to reduce this model to a Kadomtsev--Petviashvili (KP) equation. The KP model includes both the KP-I and KP-II versions, which possess line and lump soliton solutions. Using KP solitons, we construct approximate solitary wave solutions on top of the stable continuous-wave solution of the original CH-NLS model, which are found to be of both the dark and anti-dark type. We also use direct numerical simulations to investigate the validity of the approximate solutions, study their evolution, as well as their head-on collisions.


title: Strategic Monte Carlo Methods for State and Parameter Estimation in High Dimensional Nonlinear Problems

abstract: In statistical data assimilation one seeks the largest maximum of the conditional probability distribution $P(\mathbf{X},\mathbf{p}|\mathbf{Y})$ of model states, $\mathbf{X}$, and parameters,$\mathbf{p}$, conditioned on observations $\mathbf{Y}$ through minimizing the `action', $A(\mathbf{X}) = -\log P(\mathbf{X},\mathbf{p}|\mathbf{Y})$. This determines the dominant contribution to the expected values of functions of $\mathbf{X}$ but does not give information about the structure of $P(\mathbf{X},\mathbf{p}|\mathbf{Y})$ away from the maximum. We introduce a Monte Carlo sampling method, called Strategic Monte Carlo (SMC) sampling, for estimating $P(\mathbf{X}, \mathbf{p}|\mathbf{Y})$ in the neighborhood of its largest maximum to remedy this limitation. SMC begins with a systematic variational annealing (VA) procedure for finding the smallest minimum of $A(\mathbf{X})$. SMC generates accurate estimates for the mean, standard deviation and other higher moments of $P(\mathbf{X},\mathbf{p}|\mathbf{Y})$. Additionally, the random search allows for an understanding of any multimodal structure that may underly the dynamics of the problem. SMC generates a gaussian probability control term based on the paths determined by VA to minimize a cost function $A(\mathbf{X},\mathbf{p})$. This probability is sampled during the Monte Carlo search of the cost function to constrain the search to high probability regions of the surface thus substantially reducing the time necessary to sufficiently explore the space.


title: Whither Tone Ring Ring?

abstract: An extremely simple model captures the essence of the interaction of a banjo tone ring with the wood rim. The large scale, low frequency resonances of the assembled system are related to the weights and resonant frequencies of the tone ring and rim separately. Very crude measurements satisfy the derived relations within about 5% for the lowest frequency modes and give qualitative agreement for the next ones on a particular, heavy-tone-ring resonator banjo. The two combined sub-systems become increasingly independent for higher frequency, shorter-lived modes. Nevertheless, the ringing sounds of the struck individual parts, which dominate the perception of their pitch and sustain, are related by the simple model to the sound when the parts are struck when combined into one.


title: Dynamic Chain Graph Models for Ordinal Time Series Data

abstract: This paper introduces sparse dynamic chain graph models for network inference in high dimensional non-Gaussian time series data. The proposed method parametrized by a precision matrix that encodes the intra time-slice conditional independence among variables at a fixed time point, and an autoregressive coefficient that contains dynamic conditional independences interactions among time series components across consecutive time steps. The proposed model is a Gaussian copula vector autoregressive model, which is used to model sparse interactions in a high-dimensional setting. Estimation is achieved via a penalized EM algorithm. In this paper, we use an efficient coordinate descent algorithm to optimize the penalized log-likelihood with the smoothly clipped absolute deviation penalty. We demonstrate our approach on simulated and genomic datasets. The method is implemented in an R package tsnetwork.


title: Time pressure and honesty in a deception game

abstract: Previous experiments have found mixed results on whether honesty is intuitive or requires deliberation. Here we add to this literature by building on prior work of Capraro (2017). We report a large study (N=1,389) manipulating time pressure vs time delay in a deception game. We find that, in this setting, people are more honest under time pressure, and that this result is not driven by confounds present in earlier work.


title: Discovery of a sub-Keplerian disk with jet around a 20Msun young star. ALMA observations of G023.01-00.41

abstract: It is well established that Solar-mass stars gain mass via disk accretion, until the mass reservoir of the disk is exhausted and dispersed, or condenses into planetesimals. Accretion disks are intimately coupled with mass ejection via polar cavities, in the form of jets and less collimated winds, which allow mass accretion through the disk by removing a substantial fraction of its angular momentum. Whether disk accretion is the mechanism leading to the formation of stars with much higher masses is still unclear. Here, we are able to build a comprehensive picture for the formation of an O-type star, by directly imaging a molecular disk which rotates and undergoes infall around the central star, and drives a molecular jet which arises from the inner disk regions. The accretion disk is truncated between 2000-3000au, it has a mass of about a tenth of the central star mass, and is infalling towards the central star at a high rate (6x10^-4 Msun/yr), as to build up a very massive object. These findings, obtained with the Atacama Large Millimeter/submillimeter Array at 700au resolution, provide observational proof that young massive stars can form via disk accretion much like Solar-mass stars.


title: Baseline Needs More Love: On Simple Word-Embedding-Based Models and Associated Pooling Mechanisms

abstract: Many deep learning architectures have been proposed to model the compositionality in text sequences, requiring a substantial number of parameters and expensive computations. However, there has not been a rigorous evaluation regarding the added value of sophisticated compositional functions. In this paper, we conduct a point-by-point comparative study between Simple Word-Embedding-based Models (SWEMs), consisting of parameter-free pooling operations, relative to word-embedding-based RNN/CNN models. Surprisingly, SWEMs exhibit comparable or even superior performance in the majority of cases considered. Based upon this understanding, we propose two additional pooling strategies over learned word embeddings: (i) a max-pooling operation for improved interpretability; and (ii) a hierarchical pooling operation, which preserves spatial (n-gram) information within text sequences. We present experiments on 17 datasets encompassing three tasks: (i) (long) document classification; (ii) text sequence matching; and (iii) short text tasks, including classification and tagging. The source code and datasets can be obtained from https:// github.com/dinghanshen/SWEM.


title: Neutrino and Collider Implications of a Left-Right Extended Zee Model

abstract: We study a simple left-right symmetric (LRS) extension of the Zee model for neutrino mass generation. An extra $SU(2)_{L/R}$ singlet charged scalar helps in generating a loop-induced Majorana mass for neutrinos in this model. The right-handed neutrinos in this case are very light of the order of a few eV to a few MeV which makes this scenario quite different from other LRS models. We have analyzed the scalar potential and Higgs spectrum in detail, which also play an important role for the neutrino phenomenology. We identified the parameter regions in the model which satisfy the experimentally observed neutrino masses and mixings along with other experimental constraints. We have then studied the collider signatures of the charged scalar at $e^+e^-$ colliders with different benchmark points. It is possible to get a huge enhancement in the production cross-section of the charged scalar at lepton collider compared to the hadron colliders, resulting in a much stronger signal which can be easily observed at the upcoming ILC or CLIC experiments.


title: A cosmic dance at z~3: Detecting the host galaxies of the dual AGN system LBQS0302-0019 and Jil with HAWK-I+GRAAL

abstract: We recently discovered that the luminous radio-quiet QSO LBQS0302-0019 at z=3.286 is likely accompanied by an obscured AGN at 20kpc projected distance, which we dubbed Jil. It represents the tightest candidate obscured/unobscured dual AGN system at z>3. To verify the dual AGN scenario we obtained deep $K_s$ band (rest-frame $V$ band) imaging with the VLT/HAWK-I+GRAAL instrument at 0.4" resolution during science verification in January 2018. Indeed, we detect the individual host galaxies of the QSO and Jil with estimated stellar masses of $\log(M_\star/M_{\odot})=11.4\pm0.5$ and $\log(M_\star/M_{\odot})=10.9\pm0.5$, respectively. Near-IR spectra obtained with VLT-KMOS reveal a clear [OIII] $\lambda$5007 line detection at the location of Jil which does not contribute significantly to the $K_s$ band flux. Both observations therefore corroborate the dual AGN scenario. A comparison to Illustris cosmological simulations suggests a parent halo mass of $\log(M_\mathrm{halo}/M_{\odot})=13.2\pm0.5$ for this interacting galaxy system, corresponding to a very massive dark matter halo at that epoch.


title: X-ray tomography of extended objects: a comparison of data acquisition approaches

abstract: The penetration power of x-rays allows one to image large objects. For example, centimeter-sized specimens can be imaged with micron-level resolution using synchrotron sources. In this case, however, the limited beam diameter and detector size preclude the acquisition of the full sample in a single take, necessitating strategies for combining data from multiple regions. Object stitching involves the combination of local tomography data from overlapping regions, while projection stitching involves the collection of projections at multiple offset positions from the rotation axis followed by data merging and reconstruction. We compare these two approaches in terms of radiation dose applied to the specimen, and reconstructed image quality. Object stitching involves an easier data alignment problem, and immediate viewing of subregions before the entire dataset has been acquired. Projection stitching is more dose-efficient, and avoids certain artifacts of local tomography; however, it also involves a more difficult data assembly and alignment procedure, in that it is more sensitive to accumulative registration error.


title: Magnetic suppression of zonal flows on a beta plane

abstract: Zonal flows in rotating systems have been previously shown to be suppressed by the imposition of a background magnetic field aligned with the direction of rotation. Understanding the physics behind the suppression may be important in systems found in astrophysical fluid dynamics, such as stellar interiors. However, the mechanism of suppression has not yet been explained. In the idealized setting of a magnetized beta plane, we provide a theoretical explanation that shows how magnetic fluctuations directly counteract the growth of weak zonal flows. Two distinct calculations yield consistent conclusions. The first, which is simpler and more physically transparent, extends the Kelvin-Orr shearing wave to include magnetic fields and shows that weak, long-wavelength shear flow organizes magnetic fluctuations to absorb energy from the mean flow. The second calculation, based on the quasilinear, statistical CE2 framework, is valid for arbitrary wavelength zonal flow and predicts a self-consistent growth rate of the zonal flow. We find that a background magnetic field suppresses zonal flow if the bare Alfven frequency is comparable to or larger than the bare Rossby frequency. However, suppression can occur for even smaller magnetic fields if the resistivity is sufficiently small enough to allow sizable magnetic fluctuations. Our calculations reproduce the $\eta/B_0^2 = \text{const.}$ scaling that describes the boundary of zonation, as found in previous work, and we explicitly link this scaling to the amplitude of magnetic fluctuations.


title: Magnetic Field Induced Weyl Semimetal from Wannier-Function-based Tight-Binding Model

abstract: Weyl semimetals (WSMs) have Weyl nodes where conduction and valence bands meet in the absence of inversion or time-reversal symmetry (TRS), or both. Interesting phenomena are expected in WSMs such as the chiral magnetic effect, anomalous Hall conductivity or Nernst effect, and unique quantum oscillations. The TRS-broken WSM phase can be driven from a topological Dirac semimetal (DSM) by magnetic field or magnetic dopants, considering that DSMs have degenerate Weyl nodes stabilized by rotational symmetry, i.e. Dirac nodes. Here we develop a Wannier-function-based tight-binding (WF-TB) model to investigate the formation of Weyl nodes and nodal rings induced by B field in the topological DSM Na$_3$Bi. The field is applied along the rotational axis. Remarkably, our study based on the WF-TB model shows that upon B field each Dirac node is split into four separate Weyl nodes along the rotational axis near the Fermi level; two nodes with single chiral charge (single Weyl nodes) and two with double chiral charge (double Weyl nodes). This result is in contrast to the common belief that each Dirac node splits into two Weyl nodes. In the effective models, the existence of double Weyl nodes ensures nonzero cubic terms in momentum. We further examine the evolution of Fermi arcs at a side surface as a function of chemical potential. The number of Fermi arcs is consistent with the corresponding Fermi surface Chern numbers, corroborating our finding of the double Weyl nodes. Furthermore, our study reveals the existence of nodal rings in the mirror plane below the Fermi level upon B field. These nodal rings persist with spin-orbit coupling, in contrast to typical nodal ring semimetals. Our WF-TB model can be used to compute interesting features arising from Berry curvature such as anomalous Hall and thermal conductivities, and our findings can be applied to other topological DSMs like Cd$_3$As$_2$.


title: Closed-Form Projection Method for Regularizing a Function Defined by a Discrete Set of Noisy Data and for Estimating its Derivative and Fractional Derivative

abstract: We present a closed-form finite-dimensional projection method for regularizing a function defined by a discrete set of measurement data, which have been contaminated by random, zero mean errors, and for estimating the derivative and fractional derivative of this function by linear combinations of a few low degree trigonometric or Jacobi polynomials. Our method takes advantage of the fact that there are known infinite-dimensional singular value decompositions of the operators of integration and fractional integration.


title: A Bug Bounty Perspective on the Disclosure of Web Vulnerabilities

abstract: Bug bounties have become increasingly popular in recent years. This paper discusses bug bounties by framing these theoretically against so-called platform economy. Empirically the interest is on the disclosure of web vulnerabilities through the Open Bug Bounty (OBB) platform between 2015 and late 2017. According to the empirical results based on a dataset covering nearly 160 thousand web vulnerabilities, (i) OBB has been successful as a community-based platform for the dissemination of web vulnerabilities. The platform has also attracted many productive hackers, (ii) but there exists a large productivity gap, which likely relates to (iii) a knowledge gap and the use of automated tools for web vulnerability discovery. While the platform (iv) has been exceptionally fast to evaluate new vulnerability submissions, (v) the patching times of the web vulnerabilities disseminated have been long. With these empirical results and the accompanying theoretical discussion, the paper contributes to the small but rapidly growing amount of research on bug bounties. In addition, the paper makes a practical contribution by discussing the business models behind bug bounties from the viewpoints of platforms, ecosystems, and vulnerability markets.


title: Detection of the Galactic Warm Neutral Medium in HI 21cm absorption

abstract: We report a deep Giant Metrewave Radio Telescope (GMRT) search for Galactic H{\sc i} 21cm absorption towards the quasar B0438$-$436, yielding the detection of wide, weak H{\sc i} 21cm absorption, with a velocity-integrated H{\sc i} 21cm optical depth of $0.0188 \pm 0.0036$~km~s$^{-1}$. Comparing this with the H{\sc i} column density measured in the Parkes Galactic All-Sky Survey gives a column density-weighted harmonic mean spin temperature of $3760 \pm 365$~K, one of the highest measured in the Galaxy. This is consistent with most of the H{\sc i} along the sightline arising in the stable warm neutral medium (WNM). The low peak H{\sc i} 21cm optical depth towards B0438$-$436 implies negligible self-absorption, allowing a multi-Gaussian joint decomposition of the H{\sc i} 21cm absorption and emission spectra. This yields a gas kinetic temperature of $\rm T_k \leq (4910 \pm 1900)$~K, and a spin temperature of $\rm T_s = (1000 \pm 345)$~K for the gas that gives rise to the H{\sc i} 21cm absorption. Our data are consistent with the H{\sc i} 21cm absorption arising from either the stable WNM, with $\rm T_s \ll T_k$, $\rm T_k \approx 5000$~K, and little penetration of the background Lyman-$\alpha$ radiation field into the neutral hydrogen, or from the unstable neutral medium, with $\rm T_s \approx T_k \approx 1000\;K$.


title: Can the universe be described by a wavefunction?

abstract: Suppose we assume that in gently curved spacetime (a) causality is not violated to leading order (b) the Birkoff theorem holds to leading order and (c) CPT invariance holds. Then we argue that the `mostly empty' universe we observe around us cannot be described by an exact wavefunction $\Psi$. Rather, the weakly coupled particles we see are approximate quasiparticles arising as excitations of a `fuzz'. The `fuzz' {\it does} have an exact wavefunction $\Psi_{fuzz}$, but this exact wavefunction does not directly describe local particles. The argument proceeds by relating the cosmological setting to the black hole information paradox, and then using the small corrections theorem to show the impossibility of an exact wavefunction describing the visible universe.


title: Modular Decomposition of Graphs and the Distance Preserving Property

abstract: Given a graph $G$, a subgraph $H$ is isometric if $d_H(u,v) = d_G(u,v)$ for every pair $u,v\in V(H)$, where $d$ is the distance function. A graph $G$ is distance preserving (dp) if it has an isometric subgraph of every possible order. A graph is sequentially distance preserving (sdp) if its vertices can be ordered such that deleting the first $i$ vertices results in an isometric subgraph, for all $i\ge1$. We introduce a generalisation of the lexicographic product of graphs, which can be used to non-trivially describe graphs. This generalisation is the inverse of the modular decomposition of graphs, which divides the graph into disjoint clusters called modules. Using these operations, we give a necessary and sufficient condition for graphs to be dp. Finally, we show that the Cartesian product of a dp graph and an sdp graph is dp.


title: The fractional angular momentum realized by a cold neutral atom

abstract: Inspired by the electromagnetic duality, we propose an approach to realize the fractional angular momentum by using a cold atom which possesses a permanent magnetic dipole momentum. This atom interacts with two electric fields and is trapped by a harmonic potential which enable the motion of the atom to be planar and rotationally symmetric. We show that eigenvalues of the canonical angular momentum of the cold atom can take fractional values when the atom is cooled down to its lowest kinetic energy level. The fractional part of canonical angular momentum is dual to that of the fractional angular momenta realized by using a charged particle. Another approach of getting the fractional angular momentum is also presented. The differences between these two approaches are investigated.


title: MINLO t-channel single-top plus jet

abstract: We present a next-to-leading order accurate simulation of t-channel single-top plus jet production matched to parton showers via the POWHEG method. The calculation underlying the simulation is enhanced with a process-specific implementation of the multi-scale improved NLO (MINLO) method, such that it gives physical predictions all through phase space, including regions where the jet additional to the t-channel single-top process is unresolved. We further describe a tuning procedure for the MINLO Sudakov form factor, fitting the coefficient of the first subleading term in its exponent using an artificial neural-network. The latter tuning, implemented as a straightforward event-by-event reweighting, renders the MINLO simulation NLO accurate for t-channel single-top observables, in addition to those of the analogous single-top plus jet process.


title: Calibrating the HISA temperature: Measuring the temperature of the Riegel-Crutcher cloud

abstract: HI self absorption (HISA) clouds are clumps of cold neutral hydrogen (HI) visible in front of warm background gas, which makes them ideal places to study the properties of the cold atomic component of the interstellar medium (ISM). The Riegel-Crutcher (R-C) cloud is the most striking HISA feature in the Galaxy. It is one of the closest HISA clouds to us and is located in the direction of the Galactic Centre, which provides a bright background. High-resolution interferometric measurements have revealed the filamentary structure of this cloud, however it is difficult to accurately determine the temperature and the density of the gas without optical depth measurements. In this paper we present new HI absorption observations with the Australia Telescope Compact Array (ATCA) against 46 continuum sources behind the Riegel-Crutcher cloud to directly measure the optical depth of the cloud. We decompose the complex HI absorption spectra into Gaussian components using an automated machine learning algorithm. We find 300 Gaussian components, from which 67 are associated with the R-C cloud (0 < v_{LSR} < 10 kms, FWHM < 10 kms). Combining the new HI absorption data with HI emission data from previous surveys we calculate the spin temperature and find it to be between 20 and 80 K. Our measurements uncover a temperature gradient across the cloud with spin temperatures decreasing towards positive Galactic latitudes. We also find three new OH absorption lines associated with the cloud, which support the presence of molecular gas.


title: NNLOPS accurate predictions for $W^+W^-$ production

abstract: We present novel predictions for the production of $W^+W^-$ pairs in hadron collisions that are next-to-next-to-leading order accurate and consistently matched to a parton shower (NNLOPS). All diagrams that lead to the process $pp\to e^- \bar \nu_e\;\mu^+\nu_\mu+X$ are taken into account, thereby including spin correlations and off-shell effects. For the first time full NNLOPS accuracy is achieved for a $2\to 4$ process. We find good agreement, at the 1$\sigma$ level, with the $W^+W^-$ rates measured by ATLAS and CMS. The importance of NNLOPS predictions is evident from differential distributions sensitive to soft-gluon effects and from the large impact ($10$% and more) of including next-to-next-to-leading order corrections on top of MiNLO. We define a charge asymmetry for the $W$ bosons and the leptons in $W^+W^-$ production at the LHC, which is sensitive to the $W$ polarizations and hence can be used as a probe of new physics.


title: Product type potential on the $X\,Y$ model: selection of maximizing probability and a large deviation principle

abstract: Given an interval $[a,b]$ the associated $X\,Y$ model is the space $\Omega=[a,b]^\mathbb{N}$ with an a priori probability $\nu $ on the state space $[a,b]$. We will present here the case of the product type potential on the $X\,Y$ model and in this setting we can show the explicit expression of the equilibrium probability. We will also consider questions about Ergodic Optimization, maximizing probabilities, subactions and we will show selection of a maximizing probability, when temperature goes to zero. Finally we show a large deviation principle when temperature goes to zero and we present an explicit expression for the deviation function.


title: Measure of gap and inequalities in basic education students proficiencies

abstract: This study uses students performance on standardized tests as evidence of the quality of education and introduces a methodology based on the comparison of performance distributions to produce indicators for both the level achieved by the students and the learning gap between social groups, two inseparable dimensions of quality of education. In the first case, the study compares the distribution of the group observed with a reference distribution, which represents an ideal situation of where students should be. In the second, it compares the performance distribution of students belonging to social groups defined by socioeconomic characteristics. This article uses the Kullback-Leibler divergence to characterize the differences between the distributions. This measure takes into account types of diferences not considered by other measures and have solid conceptual justifications. The proposed methodology is used to describe the quality of Brazilian basic education using the test results applied biannually to all Brazilian students of basic education.


title: Single-electron $G^{(2)}$ function at nonzero temperatures

abstract: The single-particle state is not expected to demonstrate second-order coherence. This proposition, correct in the case of a pure quantum state, is not verified in the case of a mixed state. Here I analyze the consequences of this fact for the second-order correlation function, $G ^{(2)}$, of electrons injected on top of the Fermi sea with nonzero temperature. At zero temperature, the function $G ^{(2)}$ unambiguously demonstrates whether the injected state is a single- or a multi-particle state: $G^{(2)}_{}$ vanishes in the former case, while it does not vanish in the latter case. However, at nonzero temperatures, when the quantum state of injected electrons is a mixed state, the purely single-particle contribution makes the function $G ^{(2)}_{}$ to be non vanishing even in the case of a single-electron injection. The single-particle contribution puts the lower limit to the second-order correlation function of electrons injected into conductors at nonzero temperatures. The existence of a single-particle contribution to $G ^{(2)}_{}$ can be verified experimentally by measuring the cross-correlation electrical noise.


title: Distributed Symmetry-Breaking with Improved Vertex-Averaged Complexity

abstract: We study the distributed message-passing model in which a communication network is represented by a graph G=(V,E). Usually, the measure of complexity that is considered in this model is the worst-case complexity, which is the largest number of rounds performed by a vertex v\in V. While often this is a reasonable measure, in some occasions it does not express sufficiently well the actual performance of the algorithm. For example, an execution in which one processor performs r rounds, and all the rest perform significantly less rounds than r, has the same running time as an execution in which all processors perform the same number of rounds r. On the other hand, the latter execution is less efficient in several respects, such as energy efficiency, task execution efficiency, local-neighborhood efficiency and simulation efficiency. Consequently, a more appropriate measure is required in these cases. Recently, the vertex-averaged complexity was proposed by \cite{Feuilloley2017}. In this measure, the running time is the worst-case sum of rounds of communication performed by all of the graph's vertices, averaged over the number of vertices. Feuilloley \cite{Feuilloley2017} showed that leader-election admits an algorithm with a vertex-averaged complexity significantly better than its worst-case complexity. On the other hand, for O(1)-coloring of rings, the worst-case and vertex-averaged complexities are the same. This complexity is O\left(\log^{*}n\right) [12]. It remained open whether the vertex-averaged complexity of symmetry-breaking in general graphs can be better than the worst-case complexity. We answer this question in the affirmative, by showing a number of results with improved vertex-averaged complexity.


title: Light-induced atomic desorption of lithium

abstract: We demonstrate loading of a Li magneto-optical trap using light-induced atomic desorption. The magneto-optical trap confines up to approximately $4\times10^{4}$ $^{7}\text{Li}$ atoms with loading rates up to approximately $4\times10^{3}$ atoms per second. We study the Li desorption rate as a function of the desorption wavelength and power. The extracted wavelength threshold for desorption of Li from fused silica is approximately $470$ nm. In addition to desorption of lithium, we observe light-induced desorption of background gas molecules. The vacuum pressure increase due to the desorbed background molecules is $\lesssim50$ % and the vacuum pressure decreases back to its base value with characteristic timescales on the order of seconds when we extinguish the desorption light. By examining both the loading and decay curves of the magneto-optical trap, we are able to disentangle the trap decay rates due to background gases and desorbed lithium. Our results show that light-induced atomic desorption can be a viable Li vapor source for compact devices and sensors.


title: Fast Neural Machine Translation Implementation

abstract: This paper describes the submissions to the efficiency track for GPUs at the Workshop for Neural Machine Translation and Generation by members of the University of Edinburgh, Adam Mickiewicz University, Tilde and University of Alicante. We focus on efficient implementation of the recurrent deep-learning model as implemented in Amun, the fast inference engine for neural machine translation. We improve the performance with an efficient mini-batching algorithm, and by fusing the softmax operation with the k-best extraction algorithm. Submissions using Amun were first, second and third fastest in the GPU efficiency track.


title: Inverse Rational Control: Inferring What You Think from How You Forage

abstract: Complex behaviors are often driven by an internal model, which integrates sensory information over time and facilitates long-term planning. Inferring an agent's internal model is a crucial ingredient in social interactions (theory of mind), for imitation learning, and for interpreting neural activities of behaving agents. Here we describe a generic method to model an agent's behavior under an environment with uncertainty, and infer the agent's internal model, reward function, and dynamic beliefs. We apply our method to a simulated agent performing a naturalistic foraging task. We assume the agent behaves rationally --- that is, they take actions that optimize their subjective utility according to their understanding of the task and its relevant causal variables. We model this rational solution as a Partially Observable Markov Decision Process (POMDP) where the agent may make wrong assumptions about the task parameters. Given the agent's sensory observations and actions, we learn its internal model and reward function by maximum likelihood estimation over a set of task-relevant parameters. The Markov property of the POMDP enables us to characterize the transition probabilities between internal belief states and iteratively estimate the agent's policy using a constrained Expectation-Maximization (EM) algorithm. We validate our method on simulated agents performing suboptimally on a foraging task currently used in many neuroscience experiments, and successfully recover their internal model and reward function. Our work lays a critical foundation to discover how the brain represents and computes with dynamic beliefs.


title: Low-redshift Lyman continuum leaking galaxies with high [OIII]/[OII] ratios

abstract: We present observations with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope of five star-forming galaxies at redshifts z in the range 0.2993-0.4317 and with high emission-line flux ratios O32=[OIII]5007/[OII]3727 ~ 8-27 aiming to detect the Lyman continuum (LyC) emission. We detect LyC emission in all galaxies with the escape fractions fesc(LyC) in a range of 2-72 per cent. A narrow Ly-alpha emission line with two peaks in four galaxies and with three peaks in one object is seen in medium-resolution COS spectra with a velocity separation between the peaks Vsep varying from ~153 km/s to ~345 km/s. We find a general increase of the LyC escape fraction with increasing O32 and decreasing stellar mass M*, but with a large scatter of fesc(LyC). A tight anti-correlation is found between fesc(LyC) and Vsep making Vsep a good parameter for the indirect determination of the LyC escape fraction. We argue that one possible source driving the escape of ionizing radiation is stellar winds and radiation from hot massive stars.


title: Pooling of Causal Models under Counterfactual Fairness via Causal Judgement Aggregation

abstract: In this paper we consider the problem of combining multiple probabilistic causal models, provided by different experts, under the requirement that the aggregated model satisfy the criterion of counterfactual fairness. We build upon the work on causal models and fairness in machine learning, and we express the problem of combining multiple models within the framework of opinion pooling. We propose two simple algorithms, grounded in the theory of counterfactual fairness and causal judgment aggregation, that are guaranteed to generate aggregated probabilistic causal models respecting the criterion of fairness, and we compare their behaviors on a toy case study.


title: Vector boson star solutions with a quartic order self-interaction

abstract: We investigate boson star (BS) solutions in the Einstein-Proca theory with the quartic order self-interaction of the vector field $\lambda (A^\mu\bar{A}_\mu)^2/4$ and the mass term $\mu \bar{A}^\mu A_\mu/2$, where $A_\mu$ is the complex vector field and ${\bar A}_\mu$ is the complex conjugate of $A_\mu$, and $\lambda$ and $\mu$ are the coupling constant and the mass of the vector field, respectively. The vector BSs are characterized by the two conserved quantities, the Arnowitt-Deser-Misner (ADM) mass and the Noether charge associated with the global $U(1)$ symmetry. We show that in comparison with the case without the self-interaction $\lambda=0$, the maximal ADM mass and Noether charge increase for $\lambda>0$ and decrease for $\lambda<0$. We also show that there exists the critical central amplitude of the temporal component of the vector field above which there is no vector BS solution, and for $\lambda>0$ it can be expressed by the simple analytic expression. For a sufficiently large positive coupling $\Lambda:=M_{pl}^2\lambda /(8\pi\mu^2) \gg 1$, the maximal ADM mass and Noether charge of the vector BSs are obtained from the critical central amplitude and of ${\cal O}[\sqrt{\lambda}M_{pl}^3/\mu^2 \ln (\lambda M_{pl}^2/\mu^2)]$, which is different from that of the scalar BSs, ${\cal O}(\sqrt{\lambda_\phi}M_{pl}^3/\mu_\phi^2)$, where $\lambda_\phi$ and $\mu_\phi$ are the coupling constant and the mass of the complex scalar field.


title: All macroscopic quantum states are fragile and hard to prepare

abstract: We study the effect of local decoherence on arbitrary quantum states. Adapting techniques developed in quantum metrology, we show that the action of generic local noise processes -- though arbitrarily small -- always yields a state whose Quantum Fisher Information (QFI) with respect to local observables is linear in system size N, independent of the initial state. This implies that all macroscopic quantum states, which are characterized by a QFI that is quadratic in N, are fragile under decoherence, and cannot be maintained if the system is not perfectly isolated. We also provide analytical bounds on the effective system size, and show that the effective system size scales as the inverse of the noise parameter p for small p for all the noise channels considered, making it increasingly difficult to generate macroscopic or even mesoscopic quantum states. In turn, we also show that the preparation of a macroscopic quantum state, with respect to a conserved quantity, requires a device whose QFI is already at least as large as the one of the desired state. Given that the preparation device itself is classical and not a perfectly isolated macroscopic quantum state, the preparation device needs to be quadratically bigger than the macroscopic target state.


title: Measurement of $D_s^{\pm}$ production asymmetry in $pp$ collisions at $\sqrt{s} =7$ and 8 TeV

abstract: The inclusive $D_s^{\pm}$ production asymmetry is measured in $pp$ collisions collected by the LHCb experiment at centre-of-mass energies of $\sqrt{s} =7$ and 8 TeV. Promptly produced $D_s^{\pm}$ mesons are used, which decay as $D_s^{\pm}\to\phi\pi^{\pm}$, with $\phi\to K^+K^-$. The measurement is performed in bins of transverse momentum, $p_{\rm T}$, and rapidity, $y$, covering the range $2.5<p_{\rm T}<25.0$ GeV$/c$ and $2.0<y<4.5$. No kinematic dependence is observed. Evidence of nonzero $D_s^{\pm}$ production asymmetry is found with a significance of 3.3 standard deviations.


title: Overcoming dispersive spreading of quantum wave packets via periodic nonlinear kicking

abstract: We propose the suppression of dispersive spreading of wave packets governed by the free-space Schr\"odinger equation with a periodically pulsed nonlinear term. Using asymptotic analysis, we construct stroboscopically-dispersionless quantum states that are physically reminiscent of, but mathematically different from, the well-known one-soliton solutions of the nonlinear Schr\"odinger equation with a constant (time-independent) nonlinearity. Our analytics are strongly supported by full numerical simulations. The predicted dispersionless wave packets can move with arbitrary velocity and can be realized in experiments involving ultracold atomic gases with temporally controlled interactions.


title: Confidence Region of Singular Subspaces for Low-rank Matrix Regression

abstract: Low-rank matrix regression refers to the instances of recovering a low-rank matrix based on specially designed measurements and the corresponding noisy outcomes. In the last decade, numerous statistical methodologies have been developed for efficiently recovering the unknown low-rank matrices. However, in some applications, the unknown singular subspace is scientifically more important than the low-rank matrix itself. In this article, we revisit the low-rank matrix regression model and introduce a two-step procedure to construct confidence regions of the singular subspace. The procedure involves the de-biasing for the typical low-rank estimators after which we calculate the empirical singular vectors. We investigate the distribution of the joint projection distance between the empirical singular subspace and the unknown true singular subspace. We specifically prove the asymptotical normality of the joint projection distance with data-dependent centering and normalization when $r^{3/2}(m_1+m_2)^{3/2}=o(n/\log n)$ where $m_1, m_2$ denote the matrix row and column sizes, $r$ is the rank and $n$ is the number of independent random measurements. Consequently, we propose data-dependent confidence regions of the true singular subspace which attains any pre-determined confidence level asymptotically. In addition, non-asymptotical convergence rates are also established. Numerical results are presented to demonstrate the merits of our methods.

categories: math.ST cs.IT math.IT stat.ME stat.ML stat.TH


title: An electron acceptor molecule in a nanomesh: F4TCNQ on h-BN/Rh(111)

abstract: The adsorption of molecules on surfaces affects the surface dipole and thus changes in the work function may be expected. The effect in change of work function is particularly strong if charge between substrate and adsorbate is involved. Here we report the deposition of a strong electron acceptor molecule, tetrafluorotetracyanoquinodimethane C$_{12}$F$_4$N$_4$ (F$_{4}$TCNQ) on a monolayer of hexagonal boron nitride nanomesh ($h$-BN on Rh(111)). The work function of the F$_{4}$TCNQ/$h$-BN/Rh system increases upon increasing molecular coverage. The magnitude of the effect indicates electron transfer from the substrate to the F$_{4}$TCNQ molecules. Density functional theory calculations confirm the work function shift and predict doubly charged F$_{4}$TCNQ$^{2-}$ in the nanomesh pores, where the $h$-BN is closest to the Rh substrate, and to have the largest binding energy there. The preferred adsorption in the pores is conjectured from a series of ultraviolet photoelectron spectroscopy data, where the $\sigma$ bands in the pores are first attenuated. Scanning tunneling microscopy measurements indicate that F$_{4}$TCNQ molecules on the nanomesh are mobile at room temperature, as "hopping" between neighboring pores is observed.


title: Concave regression: value-constrained estimation and likelihood ratio-based inference

abstract: We propose a likelihood ratio statistic for forming hypothesis tests and confidence intervals for a nonparametrically estimated univariate regression function, based on the shape restriction of concavity (alternatively, convexity). Dealing with the likelihood ratio statistic requires studying an estimator satisfying a null hypothesis, that is, studying a concave least-squares estimator satisfying a further equality constraint. We study this null hypothesis least-squares estimator (NLSE) here, and use it to study our likelihood ratio statistic. The NLSE is the solution to a convex program, and we find a set of inequality and equality constraints that characterize the solution. We also study a corresponding limiting version of the convex program based on observing a Brownian motion with drift. The solution to the limit problem is a stochastic process. We study the optimality conditions for the solution to the limit problem and find that they match those we derived for the solution to the finite sample problem. This allows us to show the limit stochastic process yields the limit distribution of the (finite sample) NLSE. We conjecture that the likelihood ratio statistic is asymptotically pivotal, meaning that it has a limit distribution with no nuisance parameters to be estimated, which makes it a very effective tool for this difficult inference problem. We provide a partial proof of this conjecture, and we also provide simulation evidence strongly supporting this conjecture.


title: Learning Nonlinear Brain Dynamics: van der Pol Meets LSTM

abstract: Many real-world data sets, especially in biology, are produced by complex nonlinear dynamical systems. In this paper, we focus on brain calcium imaging (CaI) of different organisms (zebrafish and rat), aiming to build a model of joint activation dynamics in large neuronal populations, including the whole brain of zebrafish. We propose a new approach for capturing dynamics of temporal SVD components that uses the coupled (multivariate) van der Pol (VDP) oscillator, a nonlinear ordinary differential equation (ODE) model describing neural activity, with a new parameter estimation technique that combines variable projection optimization and stochastic search. We show that the approach successfully handles nonlinearities and hidden state variables in the coupled VDP. The approach is accurate, achieving 0.82 to 0.94 correlation between the actual and model-generated components, and interpretable, as VDP's coupling matrix reveals anatomically meaningful positive (excitatory) and negative (inhibitory) interactions across different brain subsystems corresponding to spatial SVD components. Moreover, VDP is comparable to (or sometimes better than) recurrent neural networks (LSTM) for (short-term) prediction of future brain activity; VDP needs less parameters to train, which was a plus on our small training data. Finally, the overall best predictive method, greatly outperforming both VDP and LSTM in short- and long-term predictive settings on both datasets, was the new hybrid VDP-LSTM approach that used VDP to simulate large domain-specific dataset for LSTM pretraining; note that simple LSTM data-augmentation via noisy versions of training data was much less effective.


title: Autonomous Thermalling as a Partially Observable Markov Decision Process (Extended Version)

abstract: Small uninhabited aerial vehicles (sUAVs) commonly rely on active propulsion to stay airborne, which limits flight time and range. To address this, autonomous soaring seeks to utilize free atmospheric energy in the form of updrafts (thermals). However, their irregular nature at low altitudes makes them hard to exploit for existing methods. We model autonomous thermalling as a POMDP and present a receding-horizon controller based on it. We implement it as part of ArduPlane, a popular open-source autopilot, and compare it to an existing alternative in a series of live flight tests involving two sUAVs thermalling simultaneously, with our POMDP-based controller showing a significant advantage.


title: Testing small study effects in multivariate meta-analysis

abstract: Small study effects occur when smaller studies show different, often larger, treatment effects than large ones, which may threaten the validity of systematic reviews and meta-analyses. The most well-known reasons for small study effects include publication bias, outcome reporting bias and clinical heterogeneity. Methods to account for small study effects in univariate meta-analysis have been extensively studied. However, detecting small study effects in a multivariate meta-analysis setting remains an untouched research area. One of the complications is that different types of selection processes can be involved in the reporting of multivariate outcomes. For example, some studies may be completely unpublished while others may selectively report multiple outcomes. In this paper, we propose a score test as an overall test of small study effects in multivariate meta-analysis. Two detailed case studies are given to demonstrate the advantage of the proposed test over various naive applications of univariate tests in practice. Through simulation studies, the proposed test is found to retain nominal Type I error with considerable power in moderate sample size settings. Finally, we also evaluate the concordance between the proposed test with the naive application of univariate tests by evaluating 44 systematic reviews with multiple outcomes from the Cochrane Database.


title: Online Optimization as a Feedback Controller: Stability and Tracking

abstract: This paper develops and analyzes feedback-based online optimization methods to regulate the output of a linear time-invariant (LTI) dynamical system to the optimal solution of a time-varying convex optimization problem. The design of the algorithm is based on continuous-time primal-dual dynamics, properly modified to incorporate feedback from the LTI dynamical system, applied to a proximal augmented Lagrangian function. The resultant closed-loop algorithm tracks the solution of the time-varying optimization problem without requiring knowledge of (time-varying) disturbances in the dynamical system. The analysis leverages integral quadratic constraints to provide linear matrix inequality (LMI) conditions that guarantee global exponential stability and bounded tracking error. Analytical results show that, under a sufficient time-scale separation between the dynamics of the LTI dynamical system and the algorithm, the LMI conditions can be always satisfied. The paper further proposes a modified algorithm that can track an approximate solution trajectory of the constrained optimization problem under less restrictive assumptions. As an illustrative example, the proposed algorithms are showcased for power transmission systems, to compress the time scales between secondary and tertiary control, and allow to simultaneously power re-balancing and tracking of DC optimal power flow points.


title: Hydrogen bonding in acrylamide and its role in the scattering behavior of acrylamide-based block copolymers

abstract: Hydrogen bonding plays a role in the microphase separation behavior of many block copolymers, such as those used in lithography, where the stronger interactions due to H-bonding can lead to a smaller period for the self-assembled structures, allowing the production of higher resolution templates. However, current statistical thermodynamic models used in descriptions of microphase separation, such as the Flory-Huggins approach, do not take into account some important properties of hydrogen bonding, such as site specificity and cooperativity. In this combined theoretical and experimental study, a step is taken toward the development of a more complete theory of hydrogen bonding in polymers, using polyacrylamide as a model system. We begin by developing a set of association models to describe hydrogen bonding in amides. Both models with one association constant and two association constants are considered. This theory is used to fit IR spectroscopy data from acrylamide solutions in chloroform, thereby determining the model parameters. These parameters are then employed to calculate the scattering function of the disordered state of a diblock copolymer with one polyacrylamide block and one non-hydrogen-bonding block in the random phase approximation. It is then shown that the expression for the inverse scattering function with hydrogen bonding is the same as that without hydrogen bonding, but with the Flory-Huggins parameter $\chi$ replaced by an effective value $\chi_\text{eff}=\chi+\delta\chi_\text{HB}\left(f\right)$, where the hydrogen-bonding contribution $\delta\chi_\text{HB}$ depends on the volume fraction $f$ of the hydrogen-bonding block. We find that models with two constants give better predictions of bond energy in the acrylamide dimer and more realistic asymptotic behavior of the association constants and $\delta\chi_\text{HB}$ in the limit of high temperatures.


title: Topologically protected states in $\delta$-doped junctions with band inversion

abstract: A topological boundary can be formed at the interface between a trivial and a topological insulator. The difference in the topological index across the junction leads to robust gapless surface states. Optical studies of these states are scarce in the literature, the reason being the difficulty to isolate their response from that of the bulk. In this work, we propose to deposit a $\delta$ layer of donor impurities in close proximity to a topological boundary to help detecting gapless surface states. As we will show, gapless surface states are robust against this perturbation and they enhance intraband optical transitions as measured by the oscillator strength. These results allow to understand the interplay of surface and bulk states in topological insulators.


title: On the Computational Complexity of Model Checking for Dynamic Epistemic Logic with S5 Models

abstract: Dynamic epistemic logic (DEL) is a logical framework for representing and reasoning about knowledge change for multiple agents. An important computational task in this framework is the model checking problem, which has been shown to be PSPACE-hard even for S5 models and two agents---in the presence of other features, such as multi-pointed models. We answer open questions in the literature about the complexity of this problem in more restricted settings. We provide a detailed complexity analysis of the model checking problem for DEL, where we consider various combinations of restrictions, such as the number of agents, whether the models are single-pointed or multi-pointed, and whether postconditions are allowed in the updates. In particular, we show that the problem is already PSPACE-hard in (1) the case of one agent, multi-pointed S5 models, and no postconditions, and (2) the case of two agents, only single-pointed S5 models, and no postconditions. In addition, we study the setting where only semi-private announcements are allowed as updates. We show that for this case the problem is already PSPACE-hard when restricted to two agents and three propositional variables. The results that we obtain in this paper help outline the exact boundaries of the restricted settings for which the model checking problem for DEL is computationally tractable.


title: Formation of multiple low mass stars, brown dwarfs and planemos via gravitational collapse

abstract: The origin of very low-mass stars (VLMS) and brown dwarfs (BDs) is still an unresolved topic of star formation. We here present numerical simulations of the formation of VLMS, BDs, and planet mass objects (planemos) resulting from the gravitational collapse and fragmentation of solar mass molecular cores with varying rotation rates and initial density perturbations. Our simulations yield various types of binary systems including the combinations VLMS-VLMS, BD-BD, planemo-planemo, VLMS-BD, VLMS-planemos, BD-planemo. Our scheme successfully addresses the formation of wide VLMS and BD binaries with semi-major axis up to 441 AU and produces a spectrum of mass ratios closer to the observed mass ratio distribution (q > 0.5). Molecular cores with moderate values of the ratio of kinetic to gravitational potential energy (0.16 <= beta <= 0.21) produce planemos. Solar mass cores with rotational parameters beta outside of this range yield either VLMS/BDs or a combination of both. With regard to the mass ratios we find that for both types of binary systems the mass ratio distribution varies in the range 0.31 <= q <= 0.74. We note that in the presence of radiative feedback, the length scale of fragmentation would increase by approximately two orders of magnitude, implying that the formation of binaries may be efficient for wide orbits, while being suppressed for short-orbit systems.


title: Interface characteristics in an {\alpha}+{\beta} titanium alloy

abstract: The alpha/beta interface in Ti-6Al-2Sn-4Zr-6Mo (Ti-6246) is investigated via centre of symmetry analysis, both as-grown and after 10% cold work. Semi-coherent interface steps are observed at a spacing of 4.5 +/-1.13 atoms in the as-grown condition, in good agreement with theory prediction (4.37 atoms). Lattice accommodation is observed, with elongation along [-1 2 -1 0]alpha and contraction along [1 0 -1 0]alpha . Deformed alpha exhibited larger, less coherent steps with slip bands lying in {110}beta. This indicates dislocation pile-up at the grain boundary, a precursor to globularisation, offering insight into the effect of deformation processing on the interface, which is important for titanium alloy processing route design.


title: Covering numbers for bounded variation functions

abstract: In this paper, we provide upper and lower estimates for the minimal number of functions needed to represent a bounded variation function with an accuracy of epsilon with respect to ${\bf L}^1$-distance.


title: Ridge Regression Estimated Linear Probability Model Predictions of O-glycosylation in Proteins with Structural and Sequence Data

abstract: The likelihood of O-GlcNAc glycosylation in human proteins is predicted using the ridge regression estimated linear probability model (LPM). To achieve this, sequences from three similar post-translational modifications (PTMs) of proteins occurring at, or very near, the S or T site are analyzed: N-glycosylation, O-mucin type (O-GalNAc) glycosylation, and phosphorylation. Results found include: 1) The consensus composite sequon for O-glycosylation does NOT have W on either side of the glycosylation site. 2) The same holds for the consensus sequon for phosphorylation. 3) For LPM estimation, N-glycosylated sequences are found to be good approximations to non-O-glycosylatable sequences. 4) The selective positioning of an amino acid along the sequence, differentiates the PTMs of proteins. 5) Some N-glycosylated sequences are also phosphorylated at the S or T site. 6) ASA values for N-glycosylated sequences are stochastically larger than those for O-GlcNAc glycosylated sequences. 7) Structural attributes (beta turn II, II', helix, beta bridges, beta hairpin, and the phi angle) are significant LPM predictors of O-GlcNAc glycosylation. The LPM with sequence and structural data as explanatory variables yields a Kolmogorov-Smirnov (KS) statistic value of 99%. 8) With only sequence data, the KS statistic erodes to 80%, underscoring the germaneness of structural information, which is sparse on O-glycosylated sequences. With 50% as the cutoff probability for predicting O-GlcNAc glycosylation, this LPM mispredicts 21% of out-of-sample O-GlcNAc glycosylated sequences as not being glycosylated. The 95% confidence interval around this mispredictions rate is 16% to 26%


title: The Kinetics of Primary Alpha Plate Growth in Titanium Alloys

abstract: The kinetics of primary alpha-Ti colony/Widmanstatten plate growth from the beta are examined, comparing model to experiment. The plate growth velocity depends sensitively both on the diffusivity D(T) of the rate-limiting species and on the supersaturation around the growing plate. These result in a maxima in growth velocity around 40 K below the transus, once sufficient supersaturation is available to drive plate growth. In Ti-6246, the plate growth velocity was found to be around 0.32 um min-1 at 850 oC, which was in good agreement with the model prediction of 0.36 um min-1 . The solute field around the growing plates, and the plate thickness, was found to be quite variable, due to the intergrowth of plates and soft impingement. This solute field was found to extend to up to 30 nm, and the interface concentration in the beta was found to be around 6.4 at.% Mo. It was found that increasing O content will have minimal effect on the plate lengths expected during continuous cooling; in contrast, Mo approximately doubles the plate lengths obtained for every 2 wt.% Mo reduction. Alloys using V as the beta stabiliser instead of Mo are expected to have much faster plate growth kinetics at nominally equivalent V contents. These findings will provide a useful tool for the integrated design of alloys and process routes to achieve tailored microstructures.


title: Error-Controlled Exploration of Chemical Reaction Networks with Gaussian Processes

abstract: For a theoretical understanding of the reactivity of complex chemical systems, relative energies of stationary points on potential energy hypersurfaces need to be calculated to high accuracy. Due to the large number of intermediates present in all but the simplest chemical processes, approximate quantum chemical methods are required that allow for fast evaluations of the relative energies, but at the expense of accuracy. Despite the plethora of benchmark studies, the accuracy of a quantum chemical method is often difficult to assess. Moreover, a significant improvement of a method's accuracy (e.g., through reparameterization or systematic model extension) is rarely possible. Here, we present a new approach that allows for the systematic, problem-oriented, and rolling improvement of quantum chemical results through the application of Gaussian processes. Due to its Bayesian nature, reliable error estimates are provided for each prediction. A reference method of high accuracy can be employed, if the uncertainty associated with a particular calculation is above a given threshold. The new data point is then added to a growing data set in order to continuously improve the model, and as a result, all subsequent predictions. Previous predictions are validated by the updated model to ensure that uncertainties remain within the given confidence bound, which we call backtracking. We demonstrate our approach at the example of a complex chemical reaction network.

categories: physics.chem-ph cond-mat.mtrl-sci physics.comp-ph physics.data-an


title: Super-stability in the Student-Project Allocation Problem with Ties

abstract: The Student-Project Allocation problem with lecturer preferences over Students (SPA-S) involves assigning students to projects based on student preferences over projects, lecturer preferences over students, and the maximum number of students that each project and lecturer can accommodate. This classical model assumes that each project is offered by one lecturer and that preference lists are strictly ordered. Here, we study a generalisation of SPA-S where ties are allowed in the preference lists of students and lecturers, which we refer to as the Student-Project Allocation problem with lecturer preferences over Students with Ties (SPA-ST). We investigate stable matchings under the most robust definition of stability in this context, namely super-stability. We describe the first polynomial-time algorithm to find a super-stable matching or to report that no such matching exists, given an instance of SPA-ST. Our algorithm runs in $O(L)$ time, where $L$ is the total length of all the preference lists. Finally, we present results obtained from an empirical evaluation of the linear-time algorithm based on randomly-generated SPA-ST instances. Our main finding is that, whilst super-stable matchings can be elusive when ties are present in the students' and lecturers' preference lists, the probability of such a matching existing is significantly higher if ties are restricted to the lecturers' preference lists.


title: On crossing families of complete geometric graphs

abstract: A crossing family is a collection of pairwise crossing segments, this concept was introduced by Aronov et. al. (1994). They prove that any set of $n$ points (in general position) in the plain contains a crossing family of size $\sqrt{n/12}$. In this paper we present a generalization of the concept and give several results regarding this generalization.


title: Non-homogeneous persistent random walks and averaged environment for the L\'evy-Lorentz gas

abstract: We consider transport properties for a non-homogeneous persistent random walk, that may be viewed as a mean-field version of the L\'evy-Lorentz gas, namely a 1-d model characterized by a fat polynomial tail of the distribution of scatterers' distance, with parameter $\alpha$. By varying the value of $\alpha$ we have a transition from normal transport to superdiffusion, which we characterize by appropriate continuum limits.


title: Truth, Disjunction, and Induction

abstract: By a well-known result of Kotlarski, Krajewski, and Lachlan (1981), first-order Peano arithmetic $PA$ can be conservatively extended to the theory $CT^{-}[PA]$ of a truth predicate satisfying compositional axioms, i.e., axioms stating that the truth predicate is correct on atomic formulae and commutes with all the propositional connectives and quantifiers. This results motivates the general question of determining natural axioms concerning the truth predicate that can be added to $CT^{-}[PA]$ while maintaining conservativity over $PA$. Our main result shows that conservativity fails even for the extension of $CT^{-}[PA]$ obtained by the seemingly weak axiom of disjunctive correctness $DC$ that asserts that the truth predicate commutes with disjunctions of arbitrary finite size. In particular, $CT^{-}[PA] + DC$ implies $Con(PA)$. Our main result states that the theory $CT^{-}[PA] + DC$ coincides with the theory $CT_0[PA]$ obtained by adding $\Delta_0$-induction in the language with the truth predicate. This result strengthens earlier work by Kotlarski (1986) and Cie\'sli\'nski (2010). For our proof we develop a new general form of Visser's theorem on non-existence of infinite descending chains of truth definitions and prove it by reduction to (L\"ob's version of) G\"odel's second incompleteness theorem, rather than by using the Visser-Yablo paradox, as in Visser's original proof (1989).


title: Coded FFT and Its Communication Overhead

abstract: We propose a coded computing strategy and examine communication costs of coded computing algorithms to make distributed Fast Fourier Transform (FFT) resilient to errors during the computation. We apply maximum distance separable (MDS) codes to a widely used "Transpose" algorithm for parallel FFT. In the uncoded distributed FFT algorithm, the most expensive step is a single "all-to-all" communication. We compare this with communication overhead of coding in our coded FFT algorithm. We show that by using a systematic MDS code, the communication overhead of coding is negligible in comparison with the communication costs inherent in an uncoded FFT implementation if the number of parity nodes is at most $o(\log K)$, where $K$ is the number of systematic nodes.


title: Weakly almost periodic topologies, idempotents and ideals

abstract: Let (G,tau_G) be a topological group. We establish relationships between weakly almost periodic topologies on G coarser than tau_G, central idempotents in the weakly almost periodic compactification G^W, and certain ideals in the algebra of weakly almost periodic functions W(G). We gain decompositions of weakly almost periodic representations, generalizing many from the literature. We look at the role of pre-locally compact topologies, unitarizable topologies, and extend or decompositions to Fourier-Stieltjes algebras B(G).


title: The length and depth of compact Lie groups

abstract: Let $G$ be a connected Lie group. An unrefinable chain of $G$ is a chain of subgroups $G = G_0 > G_1 > \cdots > G_t = 1$, where each $G_i$ is a maximal connected subgroup of $G_{i-1}$. In this paper, we introduce the notion of the length (respectively, depth) of $G$, defined as the maximal (respectively, minimal) length of such a chain, and we establish several new results for compact groups. In particular, we compute the exact length and depth of every compact simple Lie group, and draw conclusions for arbitrary connected compact Lie groups $G$. We obtain best possible bounds on the length of $G$ in terms of its dimension, and characterize the connected compact Lie groups that have equal length and depth. The latter result generalizes a well known theorem of Iwasawa for finite groups. More generally, we establish a best possible upper bound on $\dim G'$ in terms of the chain difference of $G$, which is its length minus its depth.


title: Engineering Structural Robustness in Power Grid Networks Susceptible to Coherent Swing Instability

abstract: Networked power grid systems are susceptible to a phenomenon known as Coherent Swing Instability (CSI), in which a subset of machines in the grid lose synchrony with the rest of the network. We develop network level evaluation metrics to (i) identify community substructures in the power grid network, (ii) determine weak points in the network that are particularly sensitive to CSI, and (iii) produce an engineering approach for the addition of transmission lines to reduce the incidences of CSI in existing networks, or design new power grid networks that are robust to CSI by their network design. For simulations on a reduced model for the American Northeast power grid, where a block of buses representing the New England region exhibit a strong propensity for CSI, we show that modifying the network's connectivity structure can markedly improve the grid's resilience to CSI. Our analysis provides a versatile diagnostic tool for evaluating the efficacy of adding lines to a power grid which is known to be prone to CSI. This is a particularly relevant problem in large-scale power systems, where improving stability and robustness to interruptions by increasing overall network connectivity is not feasible due to financial and infrastructural constraints.


title: A Comparison of Galactic Cosmic Ray Proton and Helium Nuclei Spectra From ~10 MeV/nuc to 1 TeV/nuc Using New Voyager and Higher Energy Magnetic Spectrometer Measurements - Are There Differences In the Source Spectra of The Two Nuclei

abstract: This paper determines the relative source spectra of cosmic ray H and He nuclei using a Leaky Box model for galactic propagation and the observed spectra of these nuclei from ~10 MeV/nuc to ~1 TeV/nuc. The observations consist of Voyager 1 measurements up to several hundred MeV/nuc in local interstellar space and measurements above ~10 GeV/nuc where solar modulation effects are small by experiments on BESS, PAMELA and AMS-2. Using BESS and PAMELA measurements which agree with each other, the observed spectra for H and He nuclei and the H/He ratio are well fit by source rigidity spectra for both nuclei which are ~P-2.24 over the entire range of rigidities corresponding to energies between 10 MeV/nuc and several hundred GeV/nuc. In this case, the H/He rigidity source ratio is 5.0 + 1. The recent and presumably more accurate measurements of these spectra above 10 GeV/nuc made by AMS-2 do not entirely agree with the earlier measurements, however. In particular the H spectrum is found to be steeper than that of He by about 0.10 in the spectral exponent. Using the same model for galactic propagation the AMS-2 data leads to source spectra of H and He which are ~P-2.24 up to a break rigidity ~6-8 GV. At higher rigidities the He source spectrum continues to be ~P-2.24 but the required source spectrum for H steepens to an index ~P-2.36 above ~8 GV and, as a result, the H/He source ratio decreases with increasing rigidity using the AMS-2 data.


title: Instability in a Vlasov-Fokker-Planck binary mixture

abstract: This paper is concerned with a kinetic model of a Vlasov-Fokker-Planck system used to describe the evolution of two species of particles interacting through a potential and a thermal reservoir at given temperature. We prove that at low temperature, the homogeneous equilibrium is dynamically unstable under certain perturbations. Our work is motivated by a problem arising in \cite{EGM1}.


title: Strong Thermo-mechanical Squeezing in a far detuned Membrane-in-the-middle System

abstract: We demonstrate 8.5 dB thermal squeezing of a membrane oscillator using the dynamical backaction effect and electrostatic feedback in an optomechanical membrane-in-the-middle setup. We show that strong squeezing can be obtained even in the far detuning regime of a sideband-resolved system. By using the dielectrophoretic force of a metallic needle kept in close proximity to the membrane, we implement the one-quadrature active feedback scheme to prevent the divergence of the amplified quadrature and surpass the 3 dB limit of mechanical squeezing. We also discuss different regions of the sideband spectrum where strong squeezing can be obtained. Although the demonstration here is classical, this technique is equally applicable to prepare the mechanical oscillator in a quantum squeezed state.


title: Performing Co-Membership Attacks Against Deep Generative Models

abstract: In this paper we propose a new membership attack method called co-membership attacks against deep generative models including Variational Autoencoders (VAEs) and Generative Adversarial Networks (GANs). Specifically, membership attack aims to check whether a given instance x was used in the training data or not. A co-membership attack checks whether the given bundle of n instances were in the training, with the prior knowledge that the bundle was either entirely used in the training or none at all. Successful membership attacks can compromise the privacy of training data when the generative model is published. Our main idea is to cast membership inference of target data x as the optimization of another neural network (called the attacker network) to search for the latent encoding to reproduce x. The final reconstruction error is used directly to conclude whether x was in the training data or not. We conduct extensive experiments on a variety of datasets and generative models showing that: our attacker network outperforms prior membership attacks; co-membership attacks can be substantially more powerful than single attacks; and VAEs are more susceptible to membership attacks compared to GANs.


title: Mapping the Calogero model on the Anyon model

abstract: We explicitly map the N-body one dimensional Calogero eigenstates in a harmonic well to the lowest Landau level sector of N-body eigenstates of the two dimensional anyon model in a harmonic well. The mapping is achieved in terms of a convolution kernel that uses as input the scattering eigenstates of the free Calogero model on the infinite line, which are obtained in an operator formulation.

categories: math-ph cond-mat.stat-mech hep-th math.MP


title: The impact of the cosmic variance on $H_0$ on cosmological analyses

abstract: The current $3.8 \sigma$ tension between local (arXiv:1804.10655) and global (arXiv:1605.02985) measurements of $H_0$ cannot be fully explained by the concordance $\Lambda$CDM model. It could be produced by unknown systematics or by physics beyond the standard model. On the other hand, it is well known that linear perturbation theory predicts a cosmic variance on the Hubble parameter $H_0$, which leads to systematic errors on its local determination. Here, we study how including in the likelihood the cosmic variance on $H_0$ affects statistical inference. In particular we consider the $\gamma$CDM, $w$CDM and $\gamma w$CDM parametric extensions of the standard model, which we constrain with the latest CMB, BAO, SNe Ia, RSD and $H_0$ data. We learn two important lessons. First, the systematic error from cosmic variance is - independently of the model - approximately $\sigma_{\text{cv}}\approx 0.88$ km s$^{-1}$ Mpc$^{-1}$ (1.2\% $H_0^{\text{loc}}$) when considering the redshift range $0.0233 \le z \le 0.15$, which is relative to the main analysis of (arXiv:1804.10655), and $\sigma_{\text{cv}}\approx 1.5$ km s$^{-1}$ Mpc$^{-1}$ (2.1\% $H_0^{\text{loc}}$) when considering the wider redshift range $0.01 \le z \le 0.15$. Although $\sigma_{\text{cv}}$ affects the total error budget on local $H_0$, it does not significantly alleviate the tension which remains at $\approx 3 \sigma$. Second, cosmic variance, besides shifting the constraints, can change the results of model selection: much of the statistical advantage of non-standard models is to alleviate the now-reduced tension. We conclude that, when constraining non-standard models it is important to include the cosmic variance on $H_0$ if one wants to use the local determination of the Hubble constant by Riess et al. (arXiv:1804.10655). Doing the contrary could potentially bias the conclusions.


title: Boolean Decision Rules via Column Generation

abstract: This paper considers the learning of Boolean rules in either disjunctive normal form (DNF, OR-of-ANDs, equivalent to decision rule sets) or conjunctive normal form (CNF, AND-of-ORs) as an interpretable model for classification. An integer program is formulated to optimally trade classification accuracy for rule simplicity. Column generation (CG) is used to efficiently search over an exponential number of candidate clauses (conjunctions or disjunctions) without the need for heuristic rule mining. This approach also bounds the gap between the selected rule set and the best possible rule set on the training data. To handle large datasets, we propose an approximate CG algorithm using randomization. Compared to three recently proposed alternatives, the CG algorithm dominates the accuracy-simplicity trade-off in 7 out of 15 datasets. When maximized for accuracy, CG is competitive with rule learners designed for this purpose, sometimes finding significantly simpler solutions that are no less accurate.


title: Generic Conditions for Forecast Dominance

abstract: Recent studies have analyzed whether one forecast method dominates another under a class of consistent scoring functions. While the existing literature focuses on empirical tests of forecast dominance, little is known about the theoretical conditions under which one forecast dominates another. To address this question, we derive a new characterization of dominance among forecasts of the mean functional. We present various scenarios under which dominance occurs. Unlike existing results, our results allow for the case that the forecasts' underlying information sets are not nested, and allow for uncalibrated forecasts that suffer, e.g., from model misspecification or parameter estimation error. We illustrate the empirical relevance of our results via data examples from finance and economics.


title: Impact of a primordial magnetic field on cosmic microwave background $B$ modes with weak lensing

abstract: We discuss the manner in which the primordial magnetic field (PMF) suppresses the cosmic microwave background (CMB) $B$ mode due to the weak-lensing (WL) effect. The WL effect depends on the lensing potential (LP) caused by matter perturbations, the distribution of which at cosmological scales is given by the matter power spectrum (MPS). Therefore, the WL effect on the CMB $B$ mode is affected by the MPS. Considering the effect of the ensemble average energy density of the PMF, which we call "the background PMF," on the MPS, the amplitude of MPS is suppressed in the wave number range of $k>0.01~h$ Mpc$^{-1}$.The MPS affects the LP and the WL effect in the CMB $B$ mode; however, the PMF can damp this effect. Previous studies of the CMB $B$ mode with the PMF have only considered the vector and tensor modes. These modes boost the CMB $B$ mode in the multipole range of $\ell > 1000$, whereas the background PMF damps the CMB $B$ mode owing to the WL effect in the entire multipole range. The matter density in the Universe controls the WL effect. Therefore, when we constrain the PMF and the matter density parameters from cosmological observational data sets, including the CMB $B$ mode, we expect degeneracy between these parameters. The CMB $B$ mode also provides important information on the background gravitational waves, inflation theory, matter density fluctuations, and the structure formations at the cosmological scale through the cosmological parameter search. If we study these topics and correctly constrain the cosmological parameters from cosmological observations including the CMB $B$ mode, we need to correctly consider the background PMF.


title: The VIMOS Public Extragalactic Redshift Survey (VIPERS). The complexity of galaxy populations at 0.4< z<1.3 revealed with unsupervised machine-learning algorithms

abstract: Various galaxy classification schemes have been developed so far to constrain the main physical processes regulating evolution of different galaxy types. In the era of a deluge of astrophysical information and recent progress in machine learning, a new approach to galaxy classification becomes imperative. We employ a Fisher Expectation-Maximization unsupervised algorithm working in a parameter space of 12 rest-frame magnitudes and spectroscopic redshift. The model (DBk) and the number of classes (12) were established based on the joint analysis of standard statistical criteria and confirmed by the analysis of the galaxy distribution with respect to a number of classes and their properties. This new approach allows us to classify galaxies based just on their redshifts and UV-NIR spectral energy distributions. The FEM unsupervised algorithm has automatically distinguished 12 classes: 11 classes of VIPERS galaxies and an additional class of broad-line AGNs. After a first broad division into blue, green and red categories we obtained a further sub-division into three red, three green, and five blue galaxy classes. The FEM classes follow the galaxy sequence from the earliest to the latest types that is reflected in their colours (which are constructed from rest-frame magnitudes used in classification procedure) but also their morphological, physical, and spectroscopic properties (not included in the classification scheme). We demonstrate that the members of each class share similar physical and spectral properties. In particular, we are able to find three different classes of red passive galaxy populations. Thus, we demonstrate the potential of an unsupervised approach to galaxy classification and we retrieve the complexity of galaxy populations at z~0.7, a task that usual simpler colour-based approaches cannot fulfil.


title: The VIMOS Public Extragalactic Redshift Survey (VIPERS). Unsupervised classification with photometric redshifts: a method to accurately classify large galaxy samples without spectroscopic information

abstract: Techniques to classify galaxies solely based on photometry will be necessary for future large cosmology missions, such as Euclid or LSST. However, the precision of classification is always lower in photometric surveys and can be systematically biased with respect to classifications based upon spectroscopic data. We verified how precisely the detailed classification scheme introduced by Siudek et al. (2018, hereafter: S1) for galaxies at z~0.7 could be reproduced if only photometric data are available. We applied the Fisher Expectation-Maximization (FEM) unsupervised clustering algorithm to 54,293 VIPERS galaxies working in a parameter space of reliable photometric redshifts and 12 corresponding rest-frame magnitudes. The FEM algorithm distinguishes four main groups: (1) red, (2) green, (3) blue, and (4) outliers. Each group is further divided into 3, 3, 4, and 2 subclasses, respectively. The accuracy of reproducing galaxy classes using spectroscopic data is high: 92%, 84%, 96% for red, green, and blue classes, respectively, except for dusty star-forming galaxies. The presented verification of the photometric classification demonstrates that large photometric samples can be used to distinguish different galaxy classes at z > 0.5 with an accuracy provided so far only by spectroscopic data except for particular galaxy classes.


title: Diffusion Maps for Textual Network Embedding

abstract: Textual network embedding leverages rich text information associated with the network to learn low-dimensional vectorial representations of vertices. Rather than using typical natural language processing (NLP) approaches, recent research exploits the relationship of texts on the same edge to graphically embed text. However, these models neglect to measure the complete level of connectivity between any two texts in the graph. We present diffusion maps for textual network embedding (DMTE), integrating global structural information of the graph to capture the semantic relatedness between texts, with a diffusion-convolution operation applied on the text inputs. In addition, a new objective function is designed to efficiently preserve the high-order proximity using the graph diffusion. Experimental results show that the proposed approach outperforms state-of-the-art methods on the vertex-classification and link-prediction tasks.


title: On the failure of the H\"ormander multiplier theorem in a limiting case

abstract: We discuss the H\"ormander multiplier theorem for $L^p$ boundedness of Fourier multipliers in which the multiplier belongs to a fractional Sobolev space with smoothness $s$. We show that this theorem does not hold in the limiting case $|1/p-1/2|=s/n$.


title: Current-voltage characteristic and shot noise of shift current photovoltaics

abstract: We theoretically study the current-voltage relation, $I-V$ characteristic, of the photovoltaics due to the shift current, i.e., the photocurrent generated ${\it without}$ the external dc electric field in noncentrosymmetric crystals through the Berry connection of the Bloch wavefunctions. We find that the $I-V$ characteristic and shot noise are controlled by the difference of group velocities between conduction and valence bands, i.e., $v_{11}-v_{22}$, and the relaxation time $\tau$. Since the shift current itself is independent of these quantities, there are wide possibilities to design it to maximize the energy conversion rate and also to suppress the noise. We propose that the Landau levels in noncentrosymmetric two-dimensional systems are the promising candidate for energy conversion.


title: Structure Learning from Time Series with False Discovery Control

abstract: We consider the Granger causal structure learning problem from time series data. Granger causal algorithms predict a 'Granger causal effect' between two variables by testing if prediction error of one decreases significantly in the absence of the other variable among the predictor covariates. Almost all existing Granger causal algorithms condition on a large number of variables (all but two variables) to test for effects between a pair of variables. We propose a new structure learning algorithm called MMPC-p inspired by the well known MMHC algorithm for non-time series data. We show that under some assumptions, the algorithm provides false discovery rate control. The algorithm is sound and complete when given access to perfect directed information testing oracles. We also outline a novel tester for the linear Gaussian case. We show through our extensive experiments that the MMPC-p algorithm scales to larger problems and has improved statistical power compared to existing state of the art for large sparse graphs. We also apply our algorithm on a global development dataset and validate our findings with subject matter experts.


title: Fairness GAN

abstract: In this paper, we introduce the Fairness GAN, an approach for generating a dataset that is plausibly similar to a given multimedia dataset, but is more fair with respect to protected attributes in allocative decision making. We propose a novel auxiliary classifier GAN that strives for demographic parity or equality of opportunity and show empirical results on several datasets, including the CelebFaces Attributes (CelebA) dataset, the Quick, Draw!\ dataset, and a dataset of soccer player images and the offenses they were called for. The proposed formulation is well-suited to absorbing unlabeled data; we leverage this to augment the soccer dataset with the much larger CelebA dataset. The methodology tends to improve demographic parity and equality of opportunity while generating plausible images.