title: Thermoelectric property of a one dimensional channel in the presence of a transverse magnetic field

abstract: We studied the thermal conduction through a quantum point contact (QPC), defined in GaAs-AlGaAs heterostructure, in the presence of a transverse magnetic field. A shift in the position of thermo-voltage peak is observed with increasing field. The position of the thermo-voltage peak follows the Cutler-Mott relation in the small field regime (B < 0.5 T); it starts diverging from the Cutler-Mott relation in the moderate field regime, where a cubic magnetic field term dominates over the trivial quadratic term; eventually the shift saturates in the large field regime (B > 3.0 T). Our results suggest that additional calibration is necessary when using QPC as thermometry, especially when the transverse magnetic field is applied.


title: Shear induced collective diffusivity in an emulsion of viscous drops using dynamic structure factor: effects of viscosity ratio

abstract: The shear induced collective diffusivity in an emulsion of viscous drops, specifically as a function of viscosity ratio, was numerically computed. An initially randomly packed layer of viscous drops spreading due to drop-drop interactions in an imposed shear has been simulated. The shear induced collective diffusivity coefficient was computed using a self-similar solution of the drop concentration profile. We also obtained the collective diffusivity computing the dynamic structure factor from the simulated drop positions--an analysis typically applied only to homogeneous systems. The two quantities computed using different methods are in agreement including their predictions of nonmonotonic variations with increasing capillary number and viscosity ratio. The computed values were also found to match with past measurements. The gradient diffusivity coefficient computed here was expectedly one order of magnitude larger than the self-diffusivity coefficient for a dilute emulsion previously computed using pair-wise simulation of viscous drops. Although self-diffusivity computed previously showed nonmonotonic variation with capillary number, its variation with viscosity ratio is in contrast to nonmonotonic variation of gradient diffusivity found here. The difference in variation could arise from drops not reaching equilibrium deformation between interactions--an effect absent in the pair-wise simulation used for computation of self-diffusivity--or from an intrinsic difference in physics underlying the two diffusivities. We offer a qualitative explanation of the nonmonotonic variation by relating it to average nonmonotonic drop deformation. We also provide empirical correlations of the collective diffusivity as a function of viscosity ratio and capillary number.


title: Direct evaluation of dynamical large-deviation rate functions using a variational ansatz

abstract: We describe a simple form of importance sampling designed to bound and compute large-deviation rate functions for time-extensive dynamical observables in continuous-time Markov chains. We start with a model, defined by a set of rates, and a time-extensive dynamical observable. We construct a reference model, a variational ansatz for the behavior of the original model conditioned on atypical values of the observable. Direct simulation of the reference model provides an upper bound on the large-deviation rate function associated with the original model, an estimate of the tightness of the bound, and, if the ansatz is chosen well, the exact rate function. The exact rare behavior of the original model does not need to be known in advance. We use this method to calculate rate functions for currents and counting observables in a set of network- and lattice models taken from the literature. Straightforward ansatze yield bounds that are tighter than bounds obtained from Level 2.5 of large deviations via approximations that involve uniform scalings of rates. We show how to correct these bounds in order to recover the rate functions exactly. Our approach is complementary to more specialized methods, and offers a physically transparent framework for approximating and calculating the likelihood of dynamical large deviations.


title: Disentangling genuine from matter-induced CP violation in neutrino oscillations

abstract: We prove that, in any flavor transition, neutrino oscillation CP violating asymmetries in matter have two disentangled components: (a) a CPT-odd T-invariant term, non-vanishing iff there are interactions with matter; (b) a T-odd CPT-invariant term, non-vanishing iff there is genuine CP violation. As function of the baseline, these two terms are distinct L-even and L-odd observables, respectively. In the experimental region of terrestrial accelerator neutrinos, we calculate their approximate expressions from which we prove that, at medium baselines, the CPT-odd component is small and nearly $\delta$-independent, so it can be subtracted from the experimental CP asymmetry as a theoretical background, provided the hierarchy is known. At long baselines, on the other hand, we find that (i) a Hierarchy-odd term in the CPT-odd component dominates the CP asymmetry for energies above the first oscillation node, and (ii) the CPT-odd term vanishes, independent of the CP phase $\delta$, at E = 0.92 GeV(L/1300 km) near the second oscillation maximum, where the T-odd term is almost maximal and proportional to $\sin\delta$. A measurement of the CP asymmetry in these energy regions would thus provide separate information on (i) the neutrino mass ordering, and (ii) direct evidence of genuine CP violation in the lepton sector.


title: Heavy Neutrinos in displaced vertex searches at the LHC and HL-LHC

abstract: We study the sensitivity of displaced vertex searches for heavy neutrinos produced in W boson decays in the LHC detectors ATLAS, CMS and LHCb. We also propose a new search that uses the muon chambers to detect muons from heavy neutrino decays outside the tracker. The sensitivity estimates are based on benchmark models in which the heavy neutrinos mix exclusively with one of the three Standard Model generations. In the most sensitive mass regime the displaced vertex searches can improve existing constraints on the mixing with the first two Standard Model generations by more than four orders of magnitude and by three orders of magnitude for the mixing with the third generation.


title: Exact Solution to Haldane-BCS-Hubbard Model Along the Symmetric Lines: Interaction Induced Topological Phase Transition

abstract: We propose a Haldane-BCS-Hubbard model on a honeycomb lattice, which is composed of two copies of the Haldane model of the quantum anomalous Hall effect, an equal-spin pairing term and an onsite Hubbard interaction term. For any interaction strength, this model is exactly solvable along the symmetric line where the hopping and pairing amplitudes are equal to each other. The ground state of the Haldane-BCS-Hubbard model is a topological superconducting state at weak interaction with two chiral Majorana edge states. A strong interaction drives the system across a topological quantum phase transition to a topologically trivial superconductor. A $\mathbb{Z}_{2}$ symmetry of the Hamiltonian, which is a composition of the bond-centered inversion and a gauge transformation, is spontaneously broken by the interaction, resulting a finite antiferromagnetic order in the $y$-direction.


title: The C*-algebra of compact perturbations of diagonal operators

abstract: We study the elementary C*-algebra whose elements are the sum of a diagonal plus a compact operator. We describe the structure of the unitary group, the sets of ideals, automorhisms and projections.


title: Data-driven identification of vehicle dynamics using Koopman operator

abstract: This paper presents the results of identification of vehicle dynamics using the Koopman operator. The basic idea is to transform the state space of a nonlinear system (a car in our case) to a higher-dimensional space, using so-called basis functions, where the system dynamics is linear. The selection of basis functions is crucial and there is no general approach on how to select them, this paper gives some discussion on this topic. Two distinct approaches for selecting the basis functions are presented. The first approach, based on Extended Dynamic Mode Decomposition, relies heavily on expert basis selection and is completely data-driven. The second approach utilizes the knowledge of the nonlinear dynamics, which is used to construct eigenfunctions of the Koopman operator which are known by definition to evolve linearly along the nonlinear system trajectory. The eigenfunctions are then used as basis functions for prediction. Each approach is presented with a numerical example and discussion on the feasibility of the approach for a nonlinear


title: Density shocks in the relativistic expansion of highly charged one component plasmas

abstract: In a previous paper we showed that dynamical density shocks occur in the non-relativistic expansion of dense single component plasmas relevant to ultrafast electron microscopy; and we showed that fluid models capture these effects accurately. We show that the non-relativistic decoupling of the relative and center of mass motions ceases to apply and this coupling leads to novel behavior in the relativistic dynamics under planar, cylindrical, and spherical symmetries. In cases where the relative motion of the bunch is relativistic, we show that a dynamical shock emerges even in the case of a uniform bunch with cold initial conditions; and that density shocks are in general enhanced when the relative motion becomes relativistic. Furthermore, we examine the effect of an extraction field on the relativistic dynamics of a planar symmetric bunch.


title: Multi-Robot Routing for Persistent Monitoring with Latency Constraints

abstract: In this paper we study a multi-robot path planning problem for persistent monitoring of an environment. We represent the areas to be monitored as the vertices of a weighted graph. For each vertex, there is a constraint on the maximum time spent by the robots between visits to that vertex, called the latency, and the objective is to find the minimum number of robots that can satisfy these latency constraints. The decision version of this problem is known to be PSPACE-complete. We present a $O(\log \rho)$ approximation algorithm for the problem where $\rho$ is the ratio of the maximum and the minimum latency constraints. We also present an orienteering based heuristic to solve the problem and show through simulations that in most of the cases the heuristic algorithm gives better solutions than the approximation algorithm. We evaluate our algorithms on large problem instances in a patrolling scenario and in a persistent scene reconstruction application. We also compare the algorithms with an existing solver on benchmark instances.


title: Supermassive black holes at high redshifts

abstract: MeV blazars are the most luminous persistent sources in the Universe and emit most of their energy in the MeV band. These objects display very large jet powers and accretion luminosities and are known to host black holes with a mass often exceeding $10^9 M_{\odot}$. An MeV survey, performed by a new generation MeV telescope which will bridge the entire energy and sensitivity gap between the current generation of hard X-ray and gamma-ray instruments, will detect $>$1000 MeV blazars up to a redshift of $z=5-6$. Here we show that this would allow us: 1) to probe the formation and growth mechanisms of supermassive black holes at high redshifts, 2) to pinpoint the location of the emission region in powerful blazars, 3) to determine how accretion and black hole spin interplay to power the jet.


title: A super-Earth and two sub-Neptunes transiting the nearby and quiet M dwarf TOI-270

abstract: We report the Transiting Exoplanet Survey Satellite discovery of three small planets transiting one of the nearest and brightest M dwarf hosts to date, TOI-270 (TIC 259377017; K-mag 8.3; 22.5 parsec). The M3V-type star is transited by the super-Earth-sized TOI-270 b (1.247+0.089-0.083 R_earth) and the sub-Neptune-sized exoplanets TOI-270 c (2.42+-0.13 R_earth) and TOI-270 d (2.13+-0.12 R_earth). The planets orbit close to a mean-motion resonant chain, with periods (3.36, 5.66, and 11.38 days) near ratios of small integers (5:3 and 2:1). TOI-270 is a prime target for future studies since: 1) its near-resonance allows detecting transit timing variations for precise mass measurements and dynamical studies; 2) its brightness enables independent radial velocity mass measurements; 3) the outer planets are ideal for atmospheric characterisation via transmission spectroscopy; and 4) the quiet star enables future searches for habitable zone planets. Altogether, very few systems with small, temperate exoplanets are as suitable for such complementary and detailed characterisation as TOI-270.


title: Quantum phases and topological properties of interacting fermions in one-dimensional superlattices

abstract: The realization of artificial gauge fields in ultracold atomic gases has opened up a path towards experimental studies of topological insulators and, as an ultimate goal, topological quantum matter in many-body systems. As an alternative to the direct implementation of two-dimensional lattice Hamiltonians that host the quantum Hall effect and its variants, topological charge-pumping experiments provide an additional avenue towards studying many-body systems. Here, we consider an interacting two-component gas of fermions realizing a family of one-dimensional superlattice Hamiltonians with onsite interactions and a unit cell of three sites, whose groundstates would be visited in an appropriately defined charge pump. First, we investigate the grandcanonical quantum phase diagram of individual Hamiltonians, focusing on insulating phases. For a certain commensurate filling, there is a sequence of phase transitions from a band insulator to other insulating phases (related to the physics of ionic Hubbard models) for some members of the manifold of Hamiltonians. Second, we compute the Chern numbers for the whole manifold in a many-body formulation and show that, related to the aforementioned quantum phase transitions, a topological transition results in a change of the value and sign of the Chern number. We provide both an intuitive and conceptual explanation and argue that these properties could be observed in quantum-gas experiments.


title: Necessary conditions involving Lie brackets for impulsive optimal control problems

abstract: We obtain higher order necessary conditions for a minimum of a Mayer optimal control problem connected with a nonlinear, control-affine system, where the controls range on an m-dimensional Euclidean space. Since the allowed velocities are unbounded and the absence of coercivity assumptions makes big speeds quite likely, minimizing sequences happen to converge toward "impulsive", namely discontinuous, trajectories. As is known, a distributional approach does not make sense in such a nonlinear setting, where instead a suitable embedding in the graph space is needed. We will illustrate how the chance of using impulse perturbations makes it possible to derive a Higher Order Maximum Principle which includes both the usual needle variations (in space-time) and conditions involving iterated Lie brackets. An example, where a third order necessary condition rules out the optimality of a given extremal, concludes the paper.


title: Discrete Statistical Models with Rational Maximum Likelihood Estimator

abstract: A discrete statistical model is a subset of a probability simplex. Its maximum likelihood estimator (MLE) is a retraction from that simplex onto the model. We characterize all models for which this retraction is a rational function. This is a contribution via real algebraic geometry which rests on results due to Huh and Kapranov on Horn uniformization. We present an algorithm for constructing models with rational MLE, and we demonstrate it on a range of instances. Our focus lies on models familiar to statisticians, like Bayesian networks, decomposable graphical models, and staged trees.


title: Nuclear Environments Inspection with Micro Aerial Vehicles: Algorithms and Experiments

abstract: In this work, we address the estimation, planning, control and mapping problems to allow a small quadrotor to autonomously inspect the interior of hazardous damaged nuclear sites. These algorithms run onboard on a computationally limited CPU. We investigate the effect of varying illumination on the system performance. To the best of our knowledge, this is the first fully autonomous system of this size and scale applied to inspect the interior of a full scale mock-up of a Primary Containment Vessel (PCV). The proposed solution opens up new ways to inspect nuclear reactors and to support nuclear decommissioning, which is well known to be a dangerous, long and tedious process. Experimental results with varying illumination conditions show the ability to navigate a full scale mock-up PCV pedestal and create a map of the environment, while concurrently avoiding obstacles.


title: Laser Driven Ultra-compact Undulator for Synchrotron Radiation

abstract: Laser wakefield accelerators have emerged as a promising candidate for compact synchrotron radiation and even x-ray free electron lasers. Today, to make the electrons emit electromagnetic radiation, the trajectories of laser wakefield accelerated electrons are deflected by transverse wakefield, counter-propagating laser field or external permanent magnet insertion device. Here, we propose a novel type of undulator which has a few hundred microns of period and tens of Tesla of magnetic field. The undulator consists of a bifilar capacitor-coil target which sustains strong discharge current that generates helical magnetic field around the coil axis when irradiated by a high energy laser. Coupling this undulator with state-of-the-art laser wakefield accelerators can, simultaneously, produce ultra-bright quasi-monochromatic x-rays with tunable energy ranging 5-250 keV and optimize the free electron laser parameter and gain length compared with permanent magnet based undulator. This concept may pave a path toward ultra-compact synchrotron radiation and even x-ray free electron lasers.


title: Who and When to Screen: Multi-Round Active Screening for Recurrent Infectious Diseases Under Uncertainty

abstract: Controlling recurrent infectious diseases is a vital yet complicated problem. In this paper, we propose a novel active screening model (ACTS) and algorithms to facilitate active screening for recurrent diseases (no permanent immunity) under infection uncertainty. Our contributions are: (1) A new approach to modeling multi-round network-based screening/contact tracing under uncertainty, which is a common real-life practice in a variety of diseases; (2) Two novel algorithms, Full- and Fast-REMEDY. Full-REMEDY considers the effect of future actions and finds a policy that provides high solution quality, where Fast-REMEDY scales linearly in the size of the network; (3) We evaluate Full- and Fast-REMEDY on several real-world datasets which emulate human contact and find that they control diseases better than the baselines. To the best of our knowledge, this is the first work on multi-round active screening with uncertainty for diseases with no permanent immunity.


title: State complexity of the multiples of the Thue-Morse set

abstract: The Thue-Morse set is the set of those nonnegative integers whose binary expansions have an even number of $1$. We obtain an exact formula for the state complexity of the multiplication by a constant of the Thue-Morse set $\mathcal{T}$ with respect with any base $b$ which is a power of $2$. Our proof is constructive and we are able to explicitly provide the minimal automaton of the language of all $2^p$-expansions of the set of integers $m\mathcal{T}$ for any positive integers $m$ and $p$.


title: Overdamped antiferromagnetic strange metal state in Sr$_3$IrRuO$_7$

abstract: The unconventional electronic ground state of Sr$_3$IrRuO$_7$ is explored via resonant x-ray scattering techniques and angle-resolved photoemission measurements. As the Ru content approaches $x=0.5$ in Sr$_3$(Ir$_{1-x}$Ru$_x$)$_2$O$_7$, intermediate to the $J_{eff}=1/2$ Mott state in Sr$_3$Ir$_2$O$_7$ and the quantum critical metal in Sr$_3$Ru$_2$O$_7$, a thermodynamically distinct metallic state emerges. The electronic structure of this intermediate phase lacks coherent quasiparticles, and charge transport exhibits a linear temperature dependence over a wide range of temperatures. Spin dynamics associated with the long-range antiferromagnetism of this phase show nearly local, overdamped magnetic excitations and an anomalously large energy scale of 200 meV---an energy far in excess of exchange energies present within either the Sr$_3$Ir$_2$O$_7$ or Sr$_3$Ru$_2$O$_7$ solid-solution endpoints. Overdamped quasiparticle dynamics driven by strong spin-charge coupling are proposed to explain the incoherent spectral features of the strange metal state in Sr$_3$IrRuO$_7$.


title: Scaling attractors in multi-field inflation

abstract: Multi-field inflation with a curved scalar geometry has been found to support background trajectories that violate the slow-roll, slow-turn conditions and thus have the potential to evade the swampland constraints. In order to understand how generic this novel behaviour is and what conditions lead to it, we perform a classification of dynamical attractors of two-field inflation that are of the scaling type. Scaling solutions form a one-parameter generalization of De Sitter solutions with a constant value of the first Hubble flow parameter $\epsilon$ and, as we argue and demonstrate, form a natural starting point for the study of non-slow-roll slow-turn behaviour. All scaling solutions can be classified as critical points of a specific dynamical system. We recover known multi-field inflationary attractors as approximate scaling solutions and classify their stability using dynamical system techniques. In particular, we discover that dynamical bifurcations play an integral role in the transition between geodesic and non-geodesic motion and discuss the ability of scaling solutions to describe realistic multi-field models. We revisit the criteria for background stability and show cases where the usual criteria found in the literature do not capture the background evolution of the system.


title: Deep Residual Autoencoder for quality independent JPEG restoration

abstract: In this paper we propose a deep residual autoencoder exploiting Residual-in-Residual Dense Blocks (RRDB) to remove artifacts in JPEG compressed images that is independent from the Quality Factor (QF) used. The proposed approach leverages both the learning capacity of deep residual networks and prior knowledge of the JPEG compression pipeline. The proposed model operates in the YCbCr color space and performs JPEG artifact restoration in two phases using two different autoencoders: the first one restores the luma channel exploiting 2D convolutions; the second one, using the restored luma channel as a guide, restores the chroma channels explotining 3D convolutions. Extensive experimental results on three widely used benchmark datasets (i.e. LIVE1, BDS500, and CLASSIC-5) show that our model is able to outperform the state of the art with respect to all the evaluation metrics considered (i.e. PSNR, PSNR-B, and SSIM). This results is remarkable since the approaches in the state of the art use a different set of weights for each compression quality, while the proposed model uses the same weights for all of them, making it applicable to images in the wild where the QF used for compression is unkwnown. Furthermore, the proposed model shows a greater robustness than state-of-the-art methods when applied to compression qualities not seen during training.


title: Hadron transverse momentum distributions of the Tsallis normalized and unnormalized statistics

abstract: The exact analytical formulas for the transverse momentum distributions of the Bose-Einstein, Fermi-Dirac and Maxwell-Boltzmann statistics of particles with nonzero mass in the framework of the Tsallis normalized and Tsallis unnormalized (also known as Tsallis-1 and Tsallis-2) statistics have been consistently derived. The final exact results were expressed in terms of the series expansions in the integral representation. The zeroth term approximation to both quantum and classical statistics of particles has been introduced. We have revealed that the phenomenological classical Tsallis distribution (widely used in high energy physics) is equal to the distribution of the Tsallis unnormalized statistics in the zeroth term approximation, but the phenomenological quantum Tsallis distributions (introduced by definition on the basis of the generalized entropy of the ideal gas) do not correspond to the distributions of the Tsallis statistics. We have found that in the ranges of the entropic parameter relevant to the processes of high-energy physics ($q<1$ for Tsallis-1 and $q>1$ for Tsallis-2) the Tsallis statistics is divergent. Therefore, to obtain physical results, we have regularized the Tsallis statistics by introducing an upper cut-off in the series expansion. The exact numerical results for the Bose-Einstein, Fermi-Dirac and Maxwell-Boltzmann statistics of particles in the Tsallis normalized and unnormalized statistics have been obtained. We observed that the exact results of the Tsallis statistics strongly enhanced the production of high-$p_{T}$ hadrons in comparison with the usual phenomenological Tsallis distribution function at the same values of $q$. The $q$-duality of the Tsallis normalized and unnormalized statistics for the massive particles was studied.


title: Correct Approximation of IEEE 754 Floating-Point Arithmetic for Program Verification

abstract: Verification of programs using floating-point arithmetic is challenging on several accounts. One of the difficulties of reasoning about such programs is due to the peculiarities of floating-point arithmetic: rounding errors, infinities, non-numeric objects (NaNs), signed zeroes, denormal numbers, different rounding modes, etc. One possibility to reason about floating-point arithmetic is to model a program computation path by means of a set of ternary constraints of the form z = x op y and use constraint propagation techniques to infer new information on the variables' possible values. In this setting, we define and prove the correctness of algorithms to precisely bound the value of one of the variables x, y or z, starting from the bounds known for the other two. We do this for each of the operations and for each rounding mode defined by the IEEE 754 binary floating-point standard, even in the case the rounding mode in effect is only partially known. This is the first time that such so-called filtering algorithms are defined and their correctness is formally proved. This is an important slab for paving the way to formal verification of programs that use floating-point arithmetics.


title: Counterintuitive VSM Behavior under CVR Incorporating Distribution System

abstract: This paper analyses the impact of conservation by voltage reduction (CVR) on voltage stability margin (VSM) considering transmission and distribution (T&D) systems. VSM is determined by P-V curve analysis using PSSE and GridLAB-D solvers to co-simulate the T&D systems under CVR and No CVR conditions. ZIP loads with profile [ZIP] = [0.4 0.3 0.3] are used to model the load. The paper discusses the counterintuitive result: under CVR, the VSM is reduced. Theoretical justification for the reduced VSM under CVR is the increase in the effective impedance between generation and load and this is proved using an extended 2-bus system. The paper shares T&D co-simulation results with IEEE 9-bus transmission system and a larger 123-bus distribution system and with distributed generation (DG) in unity power factor (UPF) and volt-VAR control (VVC) mode.


title: An EEG-based Stereoscopic Research to Reveal the Brain's Response to What Happens Before and After Watching 2D and 3D Movies

abstract: Despite knowing the reality of three-dimensional (3D) technology in the form of eye fatigue, this technology continues to be retained by people (especially the young community). To check what happens before and after watching a 2D and 3D movie and how this condition influences the human brain's power spectrum density (PSD), a five-member test group was arranged. In this study, electroencephalogram (EEG) was used as a neuroimaging method. EEG recordings of five individuals were taken both before and after watching 2D and 3D movies. After 2D/3D EEG recording, this record was divided into three stages for analysis. These stages consisted of Relax, Watching, and Rest. This benchmarking analysis included I) before and after watching the 2D movie (R2b and R2a), II) before and after watching the 3D movie (R3b and R3a), and III) after watching the 2D/3D movie (R2a and R3a). In the Relax and Rest stages, the 2D/3D EEG power differences in all channels of brain regions for the five EEG bands, including delta ({\delta}), theta ({\theta}), alpha ({\alpha}), beta (\b{eta}), and gamma ({\gamma}), were analyzed and compared.


title: The efficiency of quantum engines in the Poschl-Teller oscillator model

abstract: In this work, the quantum-mechanical properties of the strongly non-linear quantum oscillator described by the Poschl-Teller PT model is examined. This model has a relation with two well-known models, the free particle FP in a box and the harmonic oscillator HO, as described in our previous works. Using the PT model, a quantum-mechanical analog of the Joule-Brayton JB cycle and Otto cycle have been constructed through changes of both, the width of the well and its the quantum state. The efficiency of quantum engines based on the Poschl-Teller-like potential is derived, which is analogous to classical thermodynamic engine


title: Analyzing Occupancy-Driven Thermal Dynamics in Smart Buildings

abstract: The fact that a proper HVAC control strategy can reduce the energy consumption of a building by up to 45% has driven significant research in demand-based HVAC control. This paper presents a novel framework for modeling and analysis of thermal dynamics in smart buildings that incorporates building's thermal properties, a stochastic occupancy model and heating strategies. Each zone of a building is modeled with the help of discrete time Markov rewards formalism where the states represent the occupancy of that zone (either occupied or empty), and the state rewards incorporate the thermal dynamics and heating strategy. To demonstrate the applicability of our proposed framework, we evaluate and compare six different heating strategies for the two zone scenario of a university building. The obtained quantitative results from the PRISM probabilistic model checker show that one of the evaluated control strategies (viz. selective strategy) satisfies our requirement in terms of maintaining the occupants' comfort while being up to 13.5 times more cost effective when compared to the other evaluated strategies. Such evaluations demonstrate the framework's ability to assist in selecting the control strategy tailored around the occupancy pattern and building's thermal property.


title: On generalization of theorems of Pestryakov

abstract: In 1987, A.V. Pestryakov proved a series of theorems for cardinal functions of the space $B_{\alpha}(X)$ of all real-valued functions of Baire class $\alpha$ $(\alpha> 0)$, and he conjectured that most of these theorems are true for spaces containing all finite linear combinations of characteristic functions of zero-sets in $X$. In this paper we investigate for which theorems of Pestriakov generalizations are valid. Also we prove some additional propositions for function spaces applying the theory of selection principles.


title: Inverse wave scattering in the Laplace domain: a factorization method approach

abstract: Let $\Delta_{\Lambda}\le \lambda_{\Lambda}$ be a semi-bounded self-adjoint realization of the Laplace operator with boundary conditions (Dirichlet, Neumann, semi-transparent) assigned on the Lipschitz boundary of a bounded obstacle $\Omega$. Let $u^{\Lambda}_{f}$ and $u^{0}_{f}$ denote the solutions of the wave equations corresponding to $\Delta_{\Lambda}$ and to the free Laplacian $\Delta$ respectively, with a source term $f$ concentrated at time $t=0$ (a pulse). We show that for any fixed $\lambda>\lambda_{\Lambda}\ge 0$ and any fixed $B\subset\subset{\mathbb R}^{n}\backslash\bar\Omega$, the obstacle $\Omega$ can be reconstructed by the data $$ F^{\Lambda}_{\lambda}f(x):=\int_{0}^{\infty}e^{-\sqrt\lambda\,t}\big(u^{\Lambda}_{f}(t,x)-u^{0}_{f}(t,x)\big)\,dt\,,\qquad x\in B\,,\ f\in L^{2}({\mathbb R}^{n})\,,\ \mbox{supp}(f)\subset B\,. $$ A similar result holds in the case of screens reconstruction, when the boundary conditions are assigned only on a part of the boundary. Our method exploits the factorized form of the resolvent difference $(-\Delta_{\Lambda}+\lambda)^{-1}-(-\Delta+\lambda)^{-1}$.


title: Real monodromy action

abstract: The monodromy group is an invariant for parameterized systems of polynomial equations that encodes structure of the solutions over the parameter space. Since the structure of real solutions over real parameter spaces are of interest in many applications, real monodromy action is investigated here. A naive extension of monodromy action from the complex numbers to the real numbers is shown to be very restrictive. Therefore, we define a real monodromy structure which need not be a group but contains tiered characteristics about the real solutions. This real monodromy structure is applied to an example in kinematics which summarizes all the ways performing loops parameterized by leg lengths can cause a mechanism to change poses.


title: Water Distribution System Design Using Multi-Objective Particle Swarm Optimisation

abstract: Application of the multi-objective particle swarm optimisation (MOPSO) algorithm to design of water distribution systems is described. An earlier MOPSO algorithm is augmented with (a) local search, (b) a modified strategy for assigning the leader, and (c) a modified mutation scheme. For one of the benchmark problems described in the literature, the effect of each of the above features on the algorithm performance is demonstrated. The augmented MOPSO algorithm (called MOPSO+) is applied to five benchmark problems, and in each case, it finds non-dominated solutions not reported earlier. In addition, for the purpose of comparing Pareto fronts (sets of non-dominated solutions) obtained by different algorithms, a new criterion is suggested, and its usefulness is pointed out with an example. Finally, some suggestions regarding future research directions are made.


title: Microscopic modelling of general time-dependent quantum Markov processes

abstract: Master equations are typically adopted to describe the dynamics of open quantum systems. Such equations are either in integro-differential or in time-local form, with the latter class more frequently adopted due to the simpler numerical methods developed to obtain the corresponding solution. Here we show that any time-local master equation with positive rates in the generator, i.e. any CP-divisible quantum process, admits a microscopic model whose reduced dynamics is well described by the given equation.


title: Disjoint dominating and 2-dominating sets in graphs

abstract: A graph $G$ is a $D\!D_2$-graph if it has a pair $(D,D_2)$ of disjoint sets of vertices of $G$ such that $D$ is a dominating set and $D_2$ is a 2-dominating set of $G$. We provide several characterizations and hardness results concerning $D\!D_2$-graphs.


title: Hot exozodiacal dust: an exocometary origin?

abstract: We aim to explore two exozodiacal dust production mechanisms, first re-investigating the Poynting-Robertson drag pile-up scenario, and then elaborating on the less explored, but promising exocometary dust delivery scenario. We developped a new versatile, numerical model that calculates the dust dynamics, with non orbit-averaged equations for the grains close to the star. The model includes dust sublimation and incorporates a radiative transfer code for direct comparison to the observations. We consider in this study four stellar types, three dust compositions, and we assume a parent belt at 50 au. We find that, in the case of the Poynting-Robertson drag pile-up scenario, it is impossible to produce long-lived submicron-sized grains close to the star. The inward drifting grains fill in the region between the parent belt and the sublimation distance, producing an unrealistically strong mid-infrared excess compared to the near-infrared excess. The dust pile-up at the sublimation radius is by far insufficient to boost the near-IR flux of the exozodi to the point where it dominates over the mid-infrared excess. In the case of the exocometary dust delivery scenario, we find that a narrow ring can form close to the sublimation zone, populated with large grains several tens to several hundred of micrometers in radius. Although not perfect, this scenario provides a better match to the observations, especially if the grains are carbon-rich. We also find that the required number of active exocomets to sustain the observed dust level is reasonable. We conclude that the hot exozodiacal dust detected by near-infrared interferometry is unlikely to result from inwards grains migration by Poynting-Robertson drag from a distant parent belt, but could instead have an exocometary origin. [Abridged]


title: A differential form approach to the genus of Open Riemann surfaces

abstract: We will show that any open Riemann surface $M$ of finite genus is biholomorphic to an open set of a compact Riemann surface. Moreover, we will introduce a quotient space of forms in $M$ that determines if $M$ has finite genus and also the minimal genus where $M$ can be holomorphically embedded.


title: Arbitrage with Power Factor Correction using Energy Storage

abstract: The importance of reactive power compensation for power factor (PF) correction will significantly increase with the large-scale integration of distributed generation interfaced via inverters producing only active power. In this work, we focus on co-optimizing energy storage for performing energy arbitrage as well as local power factor correction. The joint optimization problem is non-convex, but can be solved efficiently using a McCormick relaxation along with penalty-based schemes. Using numerical simulations on real data and realistic storage profiles, we show that energy storage can correct PF locally without reducing arbitrage profit. It is observed that active and reactive power control is largely decoupled in nature for performing arbitrage and PF correction (PFC). Furthermore, we consider a real-time implementation of the problem with uncertain load, renewable and pricing profiles. We develop a model predictive control based storage control policy using auto-regressive forecast for the uncertainty. We observe that PFC is primarily governed by the size of the converter and therefore, look-ahead in time in the online setting does not affect PFC noticeably. However, arbitrage profit are more sensitive to uncertainty for batteries with faster ramp rates compared to slow ramping batteries.


title: Compression and Interpretability of Deep Neural Networks via Tucker Tensor Layer: From First Principles to Tensor Valued Back-Propagation

abstract: This work aims to help resolve the two main stumbling blocks in the application of Deep Neural Networks (DNNs), that is, the exceedingly large number of trainable parameters and their physical interpretability. This is achieved through a tensor valued approach, based on the proposed Tucker Tensor Layer (TTL), as an alternative to the dense weight-matrices of DNNs. This allows us to treat the weight-matrices of general DNNs as a matrix unfolding of a higher order weight-tensor. By virtue of the compression properties of tensor decompositions, this enables us to introduce a novel and efficient framework for exploiting the multi-way nature of the weight-tensor in order to dramatically reduce the number of DNN parameters. We also derive the tensor valued back-propagation algorithm within the TTL framework, by extending the notion of matrix derivatives to tensors. In this way, the physical interpretability of the Tucker decomposition is exploited to gain physical insights into the NN training, through the process of computing gradients with respect to each factor matrix. The proposed framework is validated on both synthetic data, and the benchmark datasets MNIST, Fashion-MNIST, and CIFAR-10. Overall, through the ability to provide the relative importance of each data feature in training, the TTL back-propagation is shown to help mitigate the "black-box" nature inherent to NNs. Experiments also illustrate that the TTL achieves a 66.63-fold compression on MNIST and Fashion-MNIST, while, by simplifying the VGG-16 network, it achieves a 10\% speed up in training time, at a comparable performance.


title: Wiretap Secret Key Capacity of Tree-PIN

abstract: We consider the problem of multiterminal secret key agreement (SKA) in wiretapped source model where terminals have access to samples of correlated random variables from a publicly known joint probability distribution. The adversary has access to a side information variable, that is correlated with terminals' variables. We focus on a special type of terminal variables in this model, known as Tree-PIN, where the relation between variables of the terminals can be represented by a tree. The study of Tree-PIN source model is of practical importance as it can be realized in wireless network environments. We derive the wiretap secret key capacity of Tree-PIN with independent leakage, and give lower and upper bounds on the maximum achievable secret key length in finite-length regime. We then prove an upper bound and a lower bound for the wiretap secret key capacity of a wiretapped PIN and give two conditions for which these bounds are tight. We also extend our main result to two other related models and prove their corresponding capacities. At the end, we argue how our analysis suggests that public interaction is required for achieving the multiterminal WSK capacity.


title: Deep Switch Networks for Generating Discrete Data and Language

abstract: Multilayer switch networks are proposed as artificial generators of high-dimensional discrete data (e.g., binary vectors, categorical data, natural language, network log files, and discrete-valued time series). Unlike deconvolution networks which generate continuous-valued data and which consist of upsampling filters and reverse pooling layers, multilayer switch networks are composed of adaptive switches which model conditional distributions of discrete random variables. An interpretable, statistical framework is introduced for training these nonlinear networks based on a maximum-likelihood objective function. To learn network parameters, stochastic gradient descent is applied to the objective. This direct optimization is stable until convergence, and does not involve back-propagation over separate encoder and decoder networks, or adversarial training of dueling networks. While training remains tractable for moderately sized networks, Markov-chain Monte Carlo (MCMC) approximations of gradients are derived for deep networks which contain latent variables. The statistical framework is evaluated on synthetic data, high-dimensional binary data of handwritten digits, and web-crawled natural language data. Aspects of the model's framework such as interpretability, computational complexity, and generalization ability are discussed.


title: Astro2020 Science White Paper: Dynamical Processes in the Planet-Forming Environment

abstract: The transfer of circumstellar disk mass and momentum onto the protostar and out into the environment occurs via a variety of mechanisms including magnetospheric accretion, jets, outflows, and disk winds. The interplay of these processes determine both the conditions under which planet formation occurs and the lifetime of the disk. Metallic emission lines, along with the Balmer series of hydrogen, probe the kinematics of gas within the planet-forming and central regions of circumstellar disks. High-spectral resolution study of these emission lines provides critical information on mass and momentum loss, turbulence, and disk wind origins.


title: Markov partitions for toral $\mathbb{Z}^2$-rotations featuring Jeandel-Rao Wang shift and model sets

abstract: We define a partition $\mathcal{P}_0$ and a $\mathbb{Z}^2$-rotation ($\mathbb{Z}^2$-action defined by rotations) on a 2-dimensional torus whose associated symbolic dynamical system is a minimal proper subshift of the Jeandel-Rao aperiodic Wang shift defined by 11 Wang tiles. We define another partition $\mathcal{P}_\mathcal{U}$ and a $\mathbb{Z}^2$-rotation on $\mathbb{T}^2$ whose associated symbolic dynamical system is equal to a minimal and aperiodic Wang shift defined by 19 Wang tiles. This proves that $\mathcal{P}_\mathcal{U}$ is a Markov partition for the $\mathbb{Z}^2$-rotation on $\mathbb{T}^2$. We prove in both cases that the toral $\mathbb{Z}^2$-rotation is the maximal equicontinuous factor of the minimal subshifts and that the set of fiber cardinalities of the factor map is $\{1,2,8\}$. The two minimal subshifts are uniquely ergodic and are isomorphic as measure-preserving dynamical systems to the toral $\mathbb{Z}^2$-rotations. It provides a construction of these Wang shifts as model sets of 4-to-2 cut and project schemes. A do-it-yourself puzzle is available in the appendix to illustrate the results.


title: On scaling limits of random trees and maps with a prescribed degree sequence

abstract: We study a configuration model on bipartite planar maps in which, given $n$ even integers, one samples a planar map with $n$ faces uniformly at random with these face degrees. We prove that when suitably rescaled, such maps always admit nontrivial subsequential limits as $n \to \infty$ in the Gromov-Hausdorff-Prokhorov topology. Further, we show that they converge in distribution towards the celebrated Brownian sphere, and more generally a Brownian disk for maps with a boundary, if and only if there is no inner face with a macroscopic degree, or, if the perimeter is too big, the maps degenerate and converge to the Brownian tree. By first sampling the degrees at random with an appropriate distribution, this model recovers that of size-conditioned Boltzmann maps associated with critical weights in the domain of attraction of a stable law with index $\alpha\in [1,2]$. The Brownian tree and disks then appear respectively in the case $\alpha=1$ and $\alpha=2$, whereas in the case $\alpha \in (1,2)$ our results partially recover previous known ones. Our proofs rely on known bijections with labelled plane trees, which are similarly sampled uniformly at random given $n$ outdegrees. Along the way, we obtain some results on the geometry of such trees, such as a convergence to the Brownian tree but only in the weaker sense of subtrees spanned by random vertices, which are of independent interest.


title: Demonstration of nonstoquastic Hamiltonian in coupled superconducting flux qubits

abstract: Quantum annealing (QA) is a heuristic algorithm for finding low-energy configurations of a system, with applications in optimization, machine learning, and quantum simulation. Up to now, all implementations of QA have been limited to qubits coupled via a single degree of freedom. This gives rise to a stoquastic Hamiltonian that has no sign problem in quantum Monte Carlo (QMC) simulations. In this paper, we report implementation and measurements of two superconducting flux qubits coupled via two canonically conjugate degrees of freedom (charge and flux) to achieve a nonstoquastic Hamiltonian. Such coupling can enhance performance of QA processors, extend the range of quantum simulations. We perform microwave spectroscopy to extract circuit parameters and show that the charge coupling manifests itself as a YY interaction in the computational basis. We observe destructive interference in quantum coherent oscillations between the computational basis states of the two-qubit system. Finally, we show that the extracted Hamiltonian is nonstoquastic over a wide range of parameters.


title: Singular Hamiltonians in models with spontaneous Lorentz symmetry breaking

abstract: Many current models which "violate Lorentz symmetry" do so via a vector or tensor field which takes on a vacuum expectation value, thereby spontaneously breaking the underlying Lorentz symmetry of the Lagrangian. One common way to construct such a model is to posit a smooth potential for this field; the natural low-energy solution of such a model would then be excepted to have the tensor field near the minimum of its potential. It is shown in this work that some such models, while appearing well-posed at the level of the Lagrangian, have a Hamiltonian which is singular on the vacuum manifold and are therefore ill-posed. I illustrate this pathology for an antisymmetric rank-2 tensor field, and find sufficient conditions under which this pathology occurs for more general field theories.


title: LiDAR-Camera Calibration under Arbitrary Configurations: Observability and Methods

abstract: LiDAR-camera calibration is a precondition for many heterogeneous systems that fuse data from LiDAR and camera. However, the constraint from common field of view and the requirement for strict time synchronization make the calibration a challenging problem. In this paper, we propose a novel LiDAR-camera calibration method aiming to eliminate these two constraints. Specifically, we capture a scan of 3D LiDAR when both the environment and the sensors are stationary, then move the camera to reconstruct the 3D environment using the sequentially obtained images. Finally, we align 3D visual points to the laser scan based on tightly couple graph optimization method to calculate the extrinsic parameters between LiDAR and camera. Under this design, the configuration of these two sensors are free from the common field of view constraint owing to the extended view from the moving camera. And we also eliminate the requirement for strict time synchronization as we only use the single scan of laser data when the sensors are stationary. We theoretically derive the conditions of minimal observability for our method and prove that the accuracy of calibration is improved by collecting more observations from multiple scattered calibration targets. We validate our method on both simulation platform and real-world datasets. Experiments show that our method achieves higher accuracy than other comparable methods, which is in accordance with our theoretical analysis. In addition, the proposed method is beneficial to not only plane measurement error based chessboard, but also other point measurement error based calibration targets, such as boxes and polygonal boards.


title: The practice of self-citations: a longitudinal study

abstract: In this article, we discuss the outcomes of an experiment where we analysed whether and to what extent the introduction, in 2012, of the new research assessment exercise in Italy (a.k.a. Italian Scientific Habilitation) affected self-citation behaviours in the Italian research community. The Italian Scientific Habilitation attests to the scientific maturity of researchers and in Italy, as in many other countries, is a requirement for accessing to a professorship. To this end, we obtained from ScienceDirect 35,673 articles published from 1957 and 2016 by the participants to the 2012 Italian Scientific Habilitation, that resulted in the extraction of 1,379,050 citations retrieved through Semantic Publishing technologies. Our analysis showed an overall increment in author self-citations (i.e. where the citing article and the cited article share at least one author) in several of the 24 academic disciplines considered. However, we depicted a stronger causal relation between such increment and the rules introduced by the 2012 Italian Scientific Habilitation in 10 out of 24 disciplines analysed.


title: Variation of stable birational types in positive characteristic

abstract: Let k be an uncountable algebraically closed field and let Y be a smooth projective k-variety which does not admit a decomposition of the diagonal. We prove that Y is not stably birational to a very general hypersurface of any given degree and dimension. We use this to study the variation of the stable birational types of Fano hypersurfaces over fields of arbitrary characteristic. This had been initiated by Shinder, whose method works in characteristic zero.


title: Finite frequency noise in chiral Luttinger liquid coupled to phonons

abstract: We study transport between Quantum Hall (QH) edge states at filling factor $\nu = 1$ in the presence of electron-acoustic-phonon coupling. Performing a Bogoliubov-Valatin (BV) trasformation the low-energy spectrum of interacting electron-phonon system is presented. The electron-phonon interaction splits the spectrum into charged and neutral "downstream" and neutral "upstream" modes with different velocities. In the regimes of dc and periodic ac biases the tunelling current and non-equilibrium finite frequency non-symmetrized noise are calculated perturbatively in tunneling coupling of quantum point contact (QPC). We show that the presence of electron-phonon interaction strongly modifies noise and current relations compared to free-fermion case.


title: On the linear twist of degree 1 functions in the extended Selberg class

abstract: Given a degree 1 function $F\in\mathcal{S}^{\sharp}$ and a real number $\alpha$, we consider the linear twist $F(s,\alpha)$, proving that it satisfies a functional equation reflecting $s$ into $1-s$, which can be seen as a Hurwitz-Lerch type of functional equation. We also derive some results on the distribution of the zeros of the linear twist.


title: Softened Goldstone-Symmetry Breaking

abstract: We propose a new way of breaking the Goldstone symmetry in composite Higgs models, restoring the global symmetry in the mixings between the elementary and composite fermions by completing the former to full representations of this symmetry. The Goldstone symmetry is in turn broken softly by vector-like mass terms in the {\it elementary} sector only. The resulting softened explicit breaking allows for a light Higgs boson, as found at the LHC, and a heavy top quark, without the need of light top partners around the Goldstone scale $f \sim {\rm TeV} \ll m_{\rm comp.}$, which remain elusive at the LHC, while we recover the standard scenario in the limit of infinite vector-like masses.


title: Dead-time correction for spectroscopic photon counting pixel detectors

abstract: Modern photon counting pixel detectors enabled a revolution in applications at synchrotron light sources and beyond in the last decade. One of the limitations of current detectors is reduced counting linearity or even paralysis at high counting rates, due to dead-time which results in photon pile-up. Existing dead-time and pile-up models fail to reproduce the complexity of dead-time effects on photon counting, resulting in empirical calibrations for particular detectors at best, imprecise linearization methods, or no linearization. This problem will increase in the future as many synchrotron light sources plan significant brilliance upgrades and free-electron lasers plan moving to a quasi-continuous operation mode. We present here the first models that use the actual behavior of the analog pre-amplifiers in spectroscopic photon counting pixel detectors with constant current discharge (e.g., Medipix family of detectors) to deduce more accurate analytical models and optimal linearization methods. In particular, for detectors with at least two counters per pixel, we completely eliminate the need of calibration, or previous knowledge of the detector and beam parameters (dead-time, integration time, large sets of synchrotron filling patterns). This is summarized in several models with increasing complexity and accuracy. Finally, we present a general empirical approach applicable to any particular cases where the analytical approach is not sufficiently precise.


title: Lifting images of standard representations of symmetric groups

abstract: We investigate closed subgroups $G \subseteq \mathrm{Sp}_{2g}(\mathbb{Z}_2)$ whose modulo-$2$ images coincide with the image $\mathfrak{S}_{2g + 1} \subseteq \mathrm{Sp}_{2g}(\mathbb{F}_2)$ of $S_{2g + 1}$ or the image $\mathfrak{S}_{2g + 2} \subseteq \mathrm{Sp}_{2g}(\mathbb{F}_2)$ of $S_{2g + 2}$ under the standard representation. We show that when $g \geq 2$, the only closed subgroup $G \subseteq \mathrm{Sp}_{2g}(\mathbb{Z}_2)$ surjecting onto $\mathfrak{S}_{2g + 2}$ is its full inverse image in $\mathrm{Sp}_{2g}(\mathbb{Z}_2)$, while all subgroups $G \subseteq \mathrm{Sp}_{2g}(\mathbb{Z}_2)$ surjecting onto $\mathfrak{S}_{2g + 1}$ are open and contain the level-$8$ principal congruence subgroup of $\mathrm{Sp}_{2g}(\mathbb{Z}_2)$. As an immediate application, we are able to strengthen a result of Zarhin on $2$-adic Galois representations associated to hyperelliptic curves. We also prove an elementary corollary concerning even-degree polynomials with full Galois group.


title: Measurement of strange baryon-antibaryon interactions with femtoscopic correlations

abstract: Two-particle correlation functions were measured for $\rm p\overline{p}$, $\rm p\overline{\Lambda}$, $\rm \overline{p}\Lambda$, and $\Lambda\overline{\Lambda}$ pairs in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=2.76$ TeV and $\sqrt{s_{\rm NN}}=5.02$ TeV recorded by the ALICE detector. From a simultaneous fit to all obtained correlation functions, real and imaginary components of the scattering lengths, as well as the effective ranges, were extracted for combined $\rm p\overline{\Lambda}$ and $\rm \overline{p}\Lambda$ pairs and, for the first time, for $\Lambda\overline{\Lambda}$ pairs. Effective averaged scattering parameters for heavier baryon-antibaryon pairs, not measured directly, are also provided. The results reveal similarly strong interaction between measured baryon-antibaryon pairs, suggesting that they all annihilate in the same manner at the same pair relative momentum $k^{*}$. Moreover, the reported significant non-zero imaginary part and negative real part of the scattering length provide motivation for future baryon-antibaryon bound state searches.


title: Looking for the Devil in the Details: Learning Trilinear Attention Sampling Network for Fine-grained Image Recognition

abstract: Learning subtle yet discriminative features (e.g., beak and eyes for a bird) plays a significant role in fine-grained image recognition. Existing attention-based approaches localize and amplify significant parts to learn fine-grained details, which often suffer from a limited number of parts and heavy computational cost. In this paper, we propose to learn such fine-grained features from hundreds of part proposals by Trilinear Attention Sampling Network (TASN) in an efficient teacher-student manner. Specifically, TASN consists of 1) a trilinear attention module, which generates attention maps by modeling the inter-channel relationships, 2) an attention-based sampler which highlights attended parts with high resolution, and 3) a feature distiller, which distills part features into a global one by weight sharing and feature preserving strategies. Extensive experiments verify that TASN yields the best performance under the same settings with the most competitive approaches, in iNaturalist-2017, CUB-Bird, and Stanford-Cars datasets.


title: Deep Reinforcement Learning with Feedback-based Exploration

abstract: Deep Reinforcement Learning has enabled the control of increasingly complex and high-dimensional problems. However, the need of vast amounts of data before reasonable performance is attained prevents its widespread application. We employ binary corrective feedback as a general and intuitive manner to incorporate human intuition and domain knowledge in model-free machine learning. The uncertainty in the policy and the corrective feedback is combined directly in the action space as probabilistic conditional exploration. As a result, the greatest part of the otherwise ignorant learning process can be avoided. We demonstrate the proposed method, Predictive Probabilistic Merging of Policies (PPMP), in combination with DDPG. In experiments on continuous control problems of the OpenAI Gym, we achieve drastic improvements in sample efficiency, final performance, and robustness to erroneous feedback, both for human and synthetic feedback. Additionally, we show solutions beyond the demonstrated knowledge.


title: Constraining Stellar Photospheres as an Essential Step for Transmission Spectroscopy of Small Exoplanets

abstract: Transmission spectra probe the atmospheres of transiting exoplanets, but these observations are also subject to signals introduced by magnetic active regions on host stars. Here we outline scientific opportunities in the next decade for providing useful constraints on stellar photospheres and inform interpretations of transmission spectra of the smallest ($R<4\,R_{\odot}$) exoplanets. We identify and discuss four primary opportunities: (1) refining stellar magnetic active region properties through exoplanet crossing events; (2) spectral decomposition of active exoplanet host stars; (3) joint retrievals of stellar photospheric and planetary atmospheric properties with studies of transmission spectra; and (4) continued visual transmission spectroscopy studies to complement longer-wavelength studies from $\textit{JWST}$. We make five recommendations to the Astro2020 Decadal Survey Committee: (1) identify the transit light source (TLS) effect as a challenge to precise exoplanet transmission spectroscopy and an opportunity ripe for scientific advancement in the coming decade; (2) include characterization of host star photospheric heterogeneity as part of a comprehensive research strategy for studying transiting exoplanets; (3) support the construction of ground-based extremely large telescopes (ELTs); (4) support multi-disciplinary research teams that bring together the heliophysics, stellar physics, and exoplanet communities to further exploit transiting exoplanets as spatial probes of stellar photospheres; and (5) support visual transmission spectroscopy efforts as complements to longer-wavelength observational campaigns with $\textit{JWST}$.


title: Cohomological representations of parahoric subgroups

abstract: We generalize a cohomological construction of representations due to Lusztig from the hyperspecial case to arbitrary parahoric subgroups of a reductive group over a local field, which splits over an unramified extension. We compute the character of these representations on certain very regular elements.


title: An Ultra Deep Field survey with WFIRST

abstract: Studying the formation and evolution of galaxies at the earliest cosmic times, and their role in reionization, requires the deepest imaging possible. Ultra-deep surveys like the HUDF and HFF have pushed to mag \mAB$\,\sim\,$30, revealing galaxies at the faint end of the LF to $z$$\,\sim\,$9$\,-\,$11 and constraining their role in reionization. However, a key limitation of these fields is their size, only a few arcminutes (less than a Mpc at these redshifts), too small to probe large-scale environments or clustering properties of these galaxies, crucial for advancing our understanding of reionization. Achieving HUDF-quality depth over areas $\sim$100 times larger becomes possible with a mission like the Wide Field Infrared Survey Telescope (WFIRST), a 2.4-m telescope with similar optical properties to HST, with a field of view of $\sim$1000 arcmin$^2$, $\sim$100$\times$ the area of the HST/ACS HUDF. This whitepaper motivates an Ultra-Deep Field survey with WFIRST, covering $\sim$100$\,-\,$300$\times$ the area of the HUDF, or up to $\sim$1 deg$^2$, to \mAB$\,\sim\,$30, potentially revealing thousands of galaxies and AGN at the faint end of the LF, at or beyond $z$\,$\sim$\,9$\,-\,$10 in the epoch of reionization, and tracing their LSS environments, dramatically increasing the discovery potential at these redshifts. (Note: This paper is a somewhat expanded version of one that was submitted as input to the Astro2020 Decadal Survey, with this version including an Appendix (which exceeded the Astro2020 page limits), describing how the science drivers for a WFIRST Ultra Deep Field might map into a notional observing program, including the filters used and exposure times needed to achieve these depths.)


title: The unpolarized two-loop massive pure singlet Wilson coefficients for deep-inelastic scattering

abstract: We calculate the massive two--loop pure singlet Wilson coefficients for heavy quark production in the unpolarized case analytically in the whole kinematic region and derive the threshold and asymptotic expansions. We also recalculate the corresponding massless two--loop Wilson coefficients. The complete expressions contain iterated integrals with elliptic letters. The contributing alphabets enlarge the Kummer-Poincar\'e letters by a series of square-root valued letters. A new class of iterated integrals, the Kummer-elliptic integrals, are introduced. For the structure functions $F_2$ and $F_L$ we also derive improved asymptotic representations adding power corrections. Numerical results are presented.


title: Weighting gates in circuit complexity and holography

abstract: Motivated by recent studies of quantum computational complexity in quantum field theory and holography, we discuss how weighting certain classes of gates building up a quantum circuit more heavily than others does affect the complexity. Utilizing Nielsen's geometric approach to circuit complexity, we investigate the effects for a regulated field theory for which the optimal circuit is a representation of $GL(N,\mathbb{R})$. More precisely, we work out how a uniformly chosen weighting factor acting on the entangling gates affects the complexity and, particularly, its divergent behavior. We show that assigning a higher cost to the entangling gates increases the complexity. Employing the penalized and the unpenalized complexities for the $\mathcal{F}_{\kappa=2}$ cost, we further find an interesting relation between the latter and the one based on the unpenalized $\mathcal{F}_{\kappa=1}$ cost. In addition, we exhibit how imposing such penalties modifies the leading order UV divergence in the complexity. We show that appropriately tuning the gate weighting eliminates the additional logarithmic factor, thus, resulting in a simple power law scaling. We also compare the circuit complexity with holographic predictions, specifically, based on the complexity=action conjecture, and relate the weighting factor to certain bulk quantities. Finally, we comment on certain expectations concerning the role of gate penalties in defining complexity in field theory and also speculate on possible implications for holography.


title: Fecundity regulation in a spatial birth-and-death process

abstract: We study a Markov birth-and-death process on a space of locally finite configurations, which describes an ecological model with a density dependent fecundity regulation mechanism. We establish existence and uniqueness of this process and analyze its properties. In particular, we show global time-space boundedness of the population density and, using a constructed Foster-Lyapunov-type function, we study return times to certain level sets of tempered configurations. We find also sufficient conditions that the degenerate invariant distribution is unique for the considered process.


title: Young massive star cluster formation in the Galactic Centre is driven by global gravitational collapse of high-mass molecular clouds

abstract: Young massive clusters (YMCs) are the most compact, high-mass stellar systems still forming at the present day. The precursor clouds to such systems are, however, rare due to their large initial gas mass reservoirs and rapid dispersal timescales due to stellar feedback. Nonetheless, unlike their high-z counterparts, these precursors are resolvable down to the sites of individually forming stars, and hence represent the ideal environments in which to test the current theories of star and cluster formation. Using high angular resolution (1$^{\prime\prime}$ / 0.05pc) and sensitivity ALMA observations of two YMC progenitor clouds in the Galactic Centre, we have identified a suite of molecular line transitions -- e.g. c-C$_{3}$H$_{2} $($7-6$) -- that are believed to be optically thin, and reliably trace the gas structure in the highest density gas on star-forming core scales. We conduct a virial analysis of the identified core and proto-cluster regions, and show that half of the cores (5/10) and both proto-clusters are unstable to gravitational collapse. This is the first kinematic evidence of global gravitational collapse in YMC precursor clouds at such an early evolutionary stage. The implications are that if these clouds are to form YMCs, then they likely do so via the "conveyor-belt" mode, whereby stars continually form within dispersed dense gas cores as the cloud undergoes global gravitational collapse. The concurrent contraction of both the cluster-scale gas and embedded (proto)stars ultimately leads to the high (proto)stellar density in YMCs.


title: $q$-Racah ensemble and $q$-P$\left(E_7^{(1)}/A_{1}^{(1)}\right)$ discrete Painlev\'e equation

abstract: The goal of this paper is to investigate the missing part of the story about the relationship between the orthogonal polynomial ensembles and Painlev\'e equations. Namely, we consider the $q$-Racah polynomial ensemble and show that the one-interval gap probabilities in this case can be expressed through a solution of the discrete $q$-P$\left(E_7^{(1)}/A_{1}^{(1)}\right)$ equation. Our approach also gives a new Lax pair for this equation. This Lax pair has an interesting additional involutive symmetry structure.


title: Gravitational dressing, soft charges, and perturbative gravitational splitting

abstract: In gauge theories and gravity, field variables are generally not gauge-invariant observables, but such observables may be constructed by "dressing" these or more general operators. Dressed operators create particles, together with their gauge or gravitational fields which typically extend to infinity. This raises an important question of how well quantum information can be localized; one version of this is the question of whether soft charges fully characterize a given localized charge or matter distribution. This paper finds expressions for the non-trivial soft charges of such dressed operators. However, a large amount of flexibility in the dressing indicates that the soft charges, and other asymptotic observables, are not inherently correlated with details of the charge or matter distribution. Instead, these asymptotic observables can be changed by adding a general radiative (source-free) field configuration to the original one. A dressing can be chosen, perturbatively, so that the asymptotic observables are independent of details of the distribution, besides its total electric or Poincar\'e charges. This provides an approach to describing localization of information in gauge theories or gravity, and thus subsystems, that avoids problems associated with nonlocality of operator subalgebras. Specifically, this construction suggests the notions of electromagnetic or gravitational splittings, which involve networks of Hilbert space embeddings in which the charges play an important role.


title: Supersolidity around a critical point in dipolar Bose Einstein condensates

abstract: We explore spatial symmetry breaking of a dipolar Bose Einstein condensate in the thermodynamic limit and reveal a critical point in the phase diagram at which crystallization occurs via a second order phase transition. This behavior is traced back to the significant effects of quantum fluctuations in dipolar condensates, which moreover stabilize a new supersolid phase, namely a regular honeycomb pattern with maximal modulational contrast and near-perfect superfluidity.


title: An analog simulation experiment to study free oscillations of a damped simple pendulum

abstract: The characteristics of drive-free oscillations of a damped simple pendulum under sinusoidal potential force field differ from those of the damped harmonic oscillations. The frequency of oscillation of a large amplitude simple pendulum decreases with increasing amplitude. Many prototype mechanical simple pendulum have been fabricated with precision and studied earlier in view of introducing them in undergraduate physics laboratories. However, fabrication and maintenance of such mechanical pendulum require special skill. In this work, we set up an analog electronic simulation experiment to serve the purpose of studying the force-free oscillations of a damped simple pendulum. We present the details of the setup and some typical results of our experiment. The experiment is simple enough to implement in undergraduate physics laboratories.


title: The soft X-ray spectrometer polarimeter SolpeX

abstract: We present a novel X-ray assembly of functionally related instrument blocks intended to measure solar flare and active region (AR) spectra from within the Russian instrument complex KORTES, to be mounted aboard the International Space Station (ISS). SolpeX consists of three blocks: fast-rotating multiple flat crystal Bragg spectrometer, pin-hole X-ray spectral imager and Bragg polarimeter. This combination of measuring blocks will offer an opportunity to detect/measure possible X-ray polarization in soft X-ray emission lines/continuum and record spectra of solar flares, in particular during their impulsive phases. Polarized Bremsstrahlung and line emission may arise from presence of directed particle beams colliding with denser regions of flares. As a result of evaporation, the X-ray spectral-components are expected to be Doppler shifted, which will also be measured. In this paper, we present details of the construction of three SolpeX blocks and discuss their functionality. Delivery of KORTES with SolpeX to ISS is expected in 2020/2021.


title: Episodic Memory Reader: Learning What to Remember for Question Answering from Streaming Data

abstract: We consider a novel question answering (QA) task where the machine needs to read from large streaming data (long documents or videos) without knowing when the questions will be given, which is difficult to solve with existing QA methods due to their lack of scalability. To tackle this problem, we propose a novel end-to-end deep network model for reading comprehension, which we refer to as Episodic Memory Reader (EMR) that sequentially reads the input contexts into an external memory, while replacing memories that are less important for answering \emph{unseen} questions. Specifically, we train an RL agent to replace a memory entry when the memory is full, in order to maximize its QA accuracy at a future timepoint, while encoding the external memory using either the GRU or the Transformer architecture to learn representations that considers relative importance between the memory entries. We validate our model on a synthetic dataset (bAbI) as well as real-world large-scale textual QA (TriviaQA) and video QA (TVQA) datasets, on which it achieves significant improvements over rule-based memory scheduling policies or an RL-based baseline that independently learns the query-specific importance of each memory.


title: Markov-chain-inspired search for MH370

abstract: Markov-chain models are constructed for the probabilistic description of the drift of marine debris from Malaysian Airlines flight MH370. En route from Kuala Lumpur to Beijing, the MH370 mysteriously disappeared in the southeastern Indian Ocean on 8 March 2014, somewhere along the arc of the 7th ping ring around the Inmarsat-3F1 satellite position when the airplane lost contact. The models are obtained by discretizing the motion of undrogued satellite-tracked surface drifting buoys from the global historical data bank. A spectral analysis, Bayesian estimation, and the computation of most probable paths between the Inmarsat arc and confirmed airplane debris beaching sites are shown to constrain the crash site, near 25$^{\circ}$S on the Inmarsat arc.

categories: stat.AP math.ST physics.ao-ph physics.pop-ph stat.TH


title: A solution to a Partitioned Balanced Tournament Design of side 11

abstract: We provide a solution to the Partitioned Balanced Tournament Design $(PBTD)$ of side $11$. The solution has been generated by exploiting combinatorial optimization methods.


title: Simulating the effect of photoheating feedback during reionization

abstract: We present self-consistent radiation hydrodynamic simulations of hydrogen reionization performed with Arepo-RT complemented by a state-of-the-art galaxy formation model. We examine how photoheating feedback, due to reionization, shapes the galaxies properties. Our fiducial model completes reionization by $z\approx6$ and matches observations of the Ly$\alpha$ forest, the CMB electron scattering optical depth, the high-redshift UV luminosity function, and stellar mass function. Contrary to previous works, photoheating suppresses star formation rates by more than $50\%$ only in halos less massive than $\sim10^{8.4}\ M_\odot$ ($\sim10^{8.8}\ M_\odot$) at $z=6$ $(z=5)$, suggesting inefficient photoheating feedback from photons within galaxies. The use of a uniform UV background that heats up the gas at $z\approx10.7$ generates an earlier onset of suppression of star formation compared to our fiducial model. This discrepancy can be mitigated by adopting a UV background model with a more realistic reionization history. In the absence of stellar feedback, photoheating alone is only able to quench halos less massive than $\sim10^9\ M_\odot$ at $z\gtrsim5$, implying that photoheating feedback is sub-dominant in regulating star formation. In addition, stellar feedback, implemented as a non-local galactic wind scheme in the simulations, weakens the strength of photoheating feedback by reducing the amount of stellar sources. Most importantly, photoheating does not leave observable imprints in the UV luminosity function, stellar mass function, or the cosmic star formation rate density. The feasibility of using these observables to detect imprints of reionization therefore requires further investigation.


title: Enhanced proximity effect in zigzag-shaped Majorana Josephson junctions

abstract: High density superconductor-semiconductor-superconductor junctions have a small induced superconducting gap due to the quasiparticle trajectories with a large momentum parallel to the junction having a very long flight time. Because a large induced gap protects Majorana modes, these long trajectories constrain Majorana devices to a low electron density. We show that a zigzag-shaped geometry eliminates these trajectories, allowing the robust creation of Majorana states with both the induced gap $E_\textrm{gap}$ and the Majorana size $\xi_\textrm{M}$ improved by more than an order of magnitude for realistic parameters. In addition to the improved robustness of Majoranas, this new zigzag geometry is insensitive to the geometric details and the device tuning.


title: Cool baryon and quark matter in holographic QCD

abstract: We establish a holographic bottom-up model which covers both the baryonic and quark matter phases in cold and dense QCD. This is obtained by including the baryons using simple approximation schemes in the V-QCD model, which also includes the backreaction of the quark matter to the dynamics of pure Yang-Mills. We examine two approaches for homogeneous baryon matter: baryons as a thin layer of noninteracting matter in the holographic bulk, and baryons with a homogeneous bulk gauge field. We find that the second approach exhibits phenomenologically reasonable features. At zero temperature, the vacuum, baryon, and quark matter phases are separated by strongly first order transitions as the chemical potential varies. The equation of state in the baryonic phase is found to be stiff, i.e., the speed of sound clearly exceeds the value $c_s^2=1/3$ of conformal plasmas at high baryon densities.


title: A complete set of in-medium splitting functions to any order in opacity

abstract: In this Letter we report the first calculation of all ${\cal O}(\alpha_s)$ medium-induced branching processes to any order in opacity. Our splitting functions results are presented as iterative solutions to matrix equations with initial conditions set by the leading order branchings in the vacuum. The flavor and quark mass dependence of the in-medium $q \rightarrow qg$, $g\rightarrow gg$, $q \rightarrow g q$, $g \rightarrow q\bar{q}$ processes is fully captured by the light-front wavefunction formalism and the color representation of the parent and daughter partons. We include the explicit solutions to second order in opacity as supplementary material and present numerical results in a realistic strongly-interacting medium produced in high center-of-mass energy heavy ion collisions at the Large Hadron Collider. Our numerical simulations show that the second order in opacity corrections can change the energy dependence of the in-medium shower intensity. We further find corrections to the longitudinal and angular distributions of the in-medium splitting kernels that may have important implications for jet substructure phenomenology.


title: Planar master integrals for four-loop form factors

abstract: We present the complete set of planar master integrals relevant to the calculation of three-point functions in four-loop massless Quantum Chromodynamics. Employing direct parametric integrations for a basis of finite integrals, we give analytic results for the Laurent expansion of conventional integrals in the parameter of dimensional regularization through to terms of weight eight.


title: The effect of rotation on the thermal instability of stratified galactic atmospheres - I. Local analysis

abstract: Observations show that (i) multiple gas phases can coexist in the atmospheres of galaxies and clusters; (ii) these atmospheres may be significantly rotating in the inner parts, with typical velocities that approach or even exceed the local sound speed. The thermal instability is a natural candidate to explain the formation of cold structures via condensation of a hotter gas phase. Here we systematically study the effect of rotation on the thermal stability of stratified plane-parallel atmospheres, using both analytical arguments and numerical simulations. We find that the formation of cold structures starting from small isobaric perturbations is enhanced in the regions where the rotation of the system is dynamically important (i.e. when the rotational velocity becomes comparable to the sound speed). In particular, the threshold value of the ratio between the cooling and dynamical time $t_{\rm cool}/t_{\rm dyn}$ below which condensations can form is increased by a factor up to $\sim 10$ in the presence of significant rotation. We briefly discuss the implications of our results for galaxies and clusters.


title: Dynamical scar states in driven fracton systems

abstract: One-dimensional fracton systems can exhibit perfect localization, failing to reach thermal equilibrium under arbitrary local unitary time evolution. We investigate how this nonergodic behavior manifests in the dynamics of a driven fracton system, specifically a one-dimensional Floquet quantum circuit model featuring conservation of a U(1) charge and its dipole moment. For a typical basis of initial conditions, a majority of states heat up to a thermal state at near-infinite temperature. In contrast, a small number of states flow to a localized steady state under the Floquet time evolution. We refer to these athermal steady states as ``dynamical scars," in analogy with the scar states observed in the spectra of certain many-body Hamiltonians. Despite their small number, these dynamical scars are experimentally relevant due to their high overlap with easily-prepared product states. Each scar state displays a single agglomerated fracton peak, in agreement with the steady-state configurations of fractonic random circuits. The details of these scars are insensitive to the precise form of the Floquet operator, which is constructed from random unitary matrices. Rather, dynamical scar states arise directly from fracton conservation laws, providing a concrete mechanism for the appearance of scars in systems with constrained quantum dynamics.

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


title: Thermal stability of winds driven by radiation pressure in super-Eddington accretion disks

abstract: Ultraluminous X-ray sources (ULXs) are mainly powered by accretion in neutron stars or stellar-mass black holes. Accreting at rates exceeding the Eddington limit by factors of a few up to hundreds, radiation pressure is expected to inflate the accretion disc, and drive fast winds that have in fact been observed at significant fractions of the speed of light. Given the super-Eddington luminosity, the accretion disc will be thicker than in sub-Eddington accretors such as common active galactic nuclei and X-ray binaries, leading to a different spectral energy distribution and, possibly, a different thermal status of the wind. Here we show the first attempt to calculate the photoionization balance of the winds driven by strong radiation pressure in thick discs with a focus on ULXs hosting black holes or non-magnetic neutron stars. We find that the winds are generally in thermally stable equilibrium, but long-term variations in the accretion rate and the inclination due to precession may have significant effects on the wind appearance and stability. Our model trends can explain the observed correlation between the spectral residuals around 1 keV and the ULX spectral state. We also find a possible correlation between the spectral hardness of the ULX, the wind velocity and the ionization parameter in support of the general scenario.


title: Band twisting and resilience to disorder in long-range topological superconductors

abstract: Planar topological superconductors with power-law-decaying pairing display different kinds of topological phase transitions where quasiparticles dubbed nonlocal-massive Dirac fermions emerge. These exotic particles form through long-range interactions between distant Majorana modes at the boundary of the system. We show how these propagating-massive Dirac fermions neither mix with bulk states nor Anderson-localize up to large amounts of static disorder despite being finite energy. Analyzing the density of states (DOS) and the band spectrum of the long-range topological superconductor, we identify the formation of an edge gap and a surprising double peak structure in the DOS which can be linked to a twisting of energy bands with nontrivial topology. Our findings are amenable to experimental verification in the near future using atom arrays on conventional superconductors, planar Josephson junctions on two-dimensional electron gases, and Floquet driving of topological superconductors.


title: Convection physics and tidal synchronization of the subdwarf binary NY Virginis

abstract: Asteroseismological analysis of NY Vir suggests that at least the outer 55 per cent of the star (in radius) rotates as a solid body and is tidally synchronized to the orbit. Detailed calculation of tidal dissipation rates in NY Vir fails to account for this synchronization. Recent observations of He core burning stars suggest that the extent of the convective core may be substantially larger than that predicted with theoretical models. We conduct a parametric investigation of sdB models generated with the Cambridge STARS code to artificially extend the radial extent of the convective core. These models with extended cores still fail to account for the synchronization. Tidal synchronization may be achievable with a non-MLT treatment of convection.


title: Electronic Ground State in Bilayer Graphene with Realistic Coulomb Interactions

abstract: Both insulating and conducting electronic behaviors have been experimentally seen in clean bilayer graphene samples at low temperature, and there is still no consensus on the nature of the interacting ground state at half-filling and in the absence of a magnetic field. Theoretically, several possibilities for the insulating ground states have been predicted for weak interaction strength. However, a recent renormalization-group calculation on a Hubbard model for charge-neutral bilayer graphene with short-range interactions suggests the emergence of low-energy Dirac fermions that would stabilize the metallic phase for weak interactions. Using a non-perturbative projective quantum Monte Carlo, we calculate the ground state for bilayer graphene using a realistic model for the Coulomb interaction that includes both short-range and long-range contributions. We find that a finite critical onsite interaction is needed to gap bilayer graphene, thereby confirming the Hubbard model expectations even in the presence of a long-range Coulomb potential, in agreement with our theoretical renormalization group analysis. In addition, we also find that the critical onsite interactions necessary to destabilize the metallic ground state decreases with increasing interlayer coupling.


title: The effect of rotation on the thermal instability of stratified galactic atmospheres - II. The formation of High Velocity Clouds

abstract: Whether High Velocity Clouds (HVCs) can form by condensation of the hot ($T \sim 10^6 \, {\rm K}$) Galactic corona as a consequence of thermal instabilities has been controversial. Here we re-examine this problem and we suggest that rotation of the corona might be a missing key ingredient. We do this by studying the evolution of the models of rotating galactic coronae presented in Sormani et al. (2018) under the presence of cooling and thermal conduction. We combine a linear stability analysis with the results of local and global hydrodynamical simulations. We find that condensations are likely to occur in regions where the corona has substantial rotational support. Under reasonably general assumptions on the rotation profile of the corona, the locations where condensations are expected are in remarkable agreement with the observed location of the major non-magellanic HVCs complexes in our Galaxy (namely, at distances $< 15 \, \rm kpc$ from the Sun and within $30^\circ$ from the disc plane). We conclude that HVCs can form by thermal instabilities provided that (i) the corona is rotating substantially in the inner ($R < 50 \, \rm kpc$) parts, as suggested by current observational data and predicted by cosmological simulations of galaxy formation; (ii) close to the disc the corona is well-represented by a nearly-equilibrium stratified rotating structure (as opposed to a fast cooling flow). Our results also suggest that a better understanding of the disc-halo interface, including supernova feedback, is critical to understand the origin of HVCs.


title: Cosmic Reionization On Computers: Reionization Histories of Present-day Galaxies

abstract: We examine the reionization history of present-day galaxies by explicitly tracing the building blocks of halos from the Cosmic Reionization On Computers project. We track dark matter particles that belong to $z=0$ halos to trace the neutral fractions at corresponding positions during rapid global reionization. The resulting particle reionization histories allow us to explore different definitions of a halo's reionization redshift and to account for the neutral content of the interstellar medium. Consistent with previous work, we find a systematic trend of reionization redshift with mass - present day halos with higher masses have earlier reionization times. Finally, we quantify the spread of reionization times within each halo, which also has a mass dependence.


title: Quantum superpositions of causal orders as an operational resource

abstract: Causal nonseparability refers to processes where events take place in a coherent superposition of different causal orders. These may be the key resource for experimental violations of causal inequalities and have been recently identified as resources for concrete information-theoretic tasks. Here, we take a step forward by deriving a complete operational framework for causal nonseparability as a resource. Our first contribution is a formal definition of quantum control of causal orders, a stronger form of causal nonseparability (with the celebrated quantum switch as best-known example) where the causal orders of events for a target system are coherently controlled by a control system. We then build a resource theory -- for both generic causal nonseparability and quantum control of causal orders -- with a physically-motivated class of free operations, based on process-matrix concatenations. We present the framework explicitly in the mindset with a control register. However, our machinery is versatile, being applicable also to scenarios with a target register alone. Moreover, an important subclass of our operations not only is free with respect to causal nonseparability and quantum control of causal orders but also preserves the very causal structure of causal processes. Hence, our treatment contains, as a built-in feature, the basis of a resource theory of quantum causal networks too. As applications, first, we establish a sufficient condition for pure-process free convertibility. This imposes a hierarchy of quantum control of causal orders with the quantum switch at the top. Second, we prove that causal-nonseparability distillation exists, i.e. we show how to convert multiple copies of a process with arbitrarily little causal nonseparability into fewer copies of a quantum switch. Our findings reveal conceptually new, unexpected phenomena, with both fundamental and practical implications.


title: Binary Red Supergiants II: Discovering and Characterizing B-type Companions

abstract: The percentage of massive main sequence OB stars in binary systems is thought to be as high as 100%. However, very few Galactic binary red supergiants (RSGs) have been identified, despite the fact that these stars are the evolved descendants of OB stars. As shown in our recent paper, binary RSGs will likely have B-type companions, as dictated by stellar evolution considerations. Such a system will have a very unique photometric signature due to the shape of the spectral energy distribution. Using photometric cutoffs it should therefore be possible to detect candidate RSG+B star binary systems. Here we present our spectoscopic follow-up observations of such candidates. Out of our initial list of 280 candidates in M31 and M33, we observed 149 and confirmed 63 as newly discovered RSG+B star binary systems. Additional spectra of four candidate systems in the Small Magellanic Cloud confirmed all of them as new RSG+B star binaries including the first known RSG+Be star system. By fitting BSTAR06 and MARCS model atmospheres to the newly-obtained spectra we place estimates on the temperatures and subtypes of both the B stars and RSGs. Overall, we have found 87 new RSG+B star binary systems in M31, M33 and the Small and Large Magellanic Clouds. Our future studies are aimed at determining the binary fraction of RSGs.


title: Cool circumgalactic gas of passive galaxies from cosmological inflow

abstract: The circumgalactic medium (CGM) of galaxies consists of a multiphase gas with components at very different temperatures, from $10^ {4}$ K to $10^ {7}$ K. One of the greatest puzzle about this medium is the presence of a large amount of low-temperature ($T\sim10^4$ K) gas around quiescent early-type galaxies (ETGs). Using semi-analytical parametric models, we describe the cool CGM around massive, low-redshift ETGs as the cosmological accretion of gas into their dark matter halos, resulting in an inflow of clouds from the external parts of the halos to the central galaxies. We compare our predictions with the observations of the COS-LRG collaboration. We find that inflow models can successfully reproduce the observed kinematics, the number of absorbers and the column densities of the cool gas. Our MCMC fit returns masses of the cool clouds of about $10^5\ \rm{M}_{\odot}$ and shows that they must evaporate during their journey due to hydrodynamic interactions with the hot gas. We conclude that the cool gas present in the halos of ETGs likely cannot reach the central regions and feed the galaxy star formation, thus explaining why these passive objects are no longer forming stars.


title: Non-equilibrium Ionization States Within Galactic Outflows: Explaining Their O VI and N V Column Densities

abstract: We present a suite of one-dimensional spherically-symmetric hydrodynamic simulations that study the atomic ionization structure of galactic outflows. We track the ionization state of the outflowing gas with a non-equilibrium atomic chemistry network that includes photoionization, photo-heating, and ion-by-ion cooling. Each simulation describes a steady-state outflow that is defined by its mass and energy input rates, sonic radius, metallicity, and UV flux from both the host galaxy and meta-galactic background. We find that for a large range of parameter choices, the ionization state of the material departs strongly from what it would be in photo-ionization equilibrium, in conflict with what is commonly assumed in the analysis of observations. In addition, nearly all the models reproduce the low N V to O VI column density ratios and the relatively high O VI column densities that are observed.