Showing 1–50 of 244 results for author: Christensen, R

Upgrading the Submillimeter Array: wSMA and beyond

Authors: Paul K. Grimes, Garrett K. Keating, Raymond Blundell, Robert D. Christensen, Mark Gurwell, Attila Kovacs, Timothy Norton, Scott N. Paine, Ramprasad Rao, Edward C. -Y. Tong, Jonathan Weintroub, David Wilner, Robert W. Wilson, Lingzhen Zeng, Qizhou Zhang

Abstract: The Submillimeter Array (SMA) is an array of 8 antennas operating at millimeter and submillimeter wavelengths on Maunakea, Hawaii, operated by the Smithsonian Astrophysical Observatory and Academia Sinica Institute of Astronomy and Astrophysics, Taiwan. Over the past several years, we have been preparing a major upgrade to the SMA that will replace the aging original receiver cryostats and receive… ▽ More The Submillimeter Array (SMA) is an array of 8 antennas operating at millimeter and submillimeter wavelengths on Maunakea, Hawaii, operated by the Smithsonian Astrophysical Observatory and Academia Sinica Institute of Astronomy and Astrophysics, Taiwan. Over the past several years, we have been preparing a major upgrade to the SMA that will replace the aging original receiver cryostats and receiver cartridges with all new cryostats and new 230 and 345 GHz receiver designs. This wideband upgrade (wSMA) will also include significantly increased instantaneous bandwidth, improved sensitivity, and greater capabilities for dual frequency observations. In this paper, we will describe the wSMA receiver upgrade and status, as well as the future upgrades that will be enabled by the deployment of the wSMA receivers. △ Less

Submitted 24 June, 2024; originally announced June 2024.

Comments: To be published in the proceedings of SPIE Astronomical Telescopes + Instrumentation 2024, paper number 13096-122

The Effects of a Stably Stratified Region with radially varying Electrical Conductivity on the Formation of Zonal Winds on Gas Planets

Authors: Paula N. Wulff, Ulrich R. Christensen, Wieland Dietrich, Johannes Wicht

Abstract: The outer areas of Jupiter and Saturn have multiple zonal winds, reaching the high latitudes, that penetrate deep into the planets' interiors, as suggested by gravity measurements. These characteristics are replicable in numerical simulations by including both a shallow stably stratified layer, below a convecting envelope, and increasing electrical conductivity. A dipolar magnetic field, assumed t… ▽ More The outer areas of Jupiter and Saturn have multiple zonal winds, reaching the high latitudes, that penetrate deep into the planets' interiors, as suggested by gravity measurements. These characteristics are replicable in numerical simulations by including both a shallow stably stratified layer, below a convecting envelope, and increasing electrical conductivity. A dipolar magnetic field, assumed to be generated by a dynamo below our model, is imposed. We find that the winds' depth into the stratified layer depends on the local product of the squared magnetic field strength and electrical conductivity. The key for the drop-off of the zonal winds is a meridional circulation which perturbs the density structure in the stable layer. In the stable region its dynamics is governed by a balance between Coriolis and electromagnetic forces. Our models suggest that a stable layer extending into weakly conducting regions could account for the observed deep zonal wind structures. △ Less

Submitted 30 January, 2024; originally announced February 2024.

Comments: Accepted by JGR - Planets

Bursting with Feedback: The Relationship between Feedback Model and Bursty Star Formation Histories in Dwarf Galaxies

Authors: Bianca Azartash-Namin, Anna Engelhardt, Ferah Munshi, B. W. Keller, Alyson M. Brooks, Jordan Van Nest, Charlotte R. Christensen, Tom Quinn, James Wadsley

Abstract: We use high-resolution cosmological simulations to compare the effect of bursty star formation histories on dwarf galaxy structure for two different subgrid supernovae (SNe) feedback models in dwarf galaxies with stellar masses from $5000 <$ M$_*$/M$_\odot$ $< 10^{9}$. Our simulations are run using two distinct supernova feedback models: superbubble and blastwave. We show that both models are capa… ▽ More We use high-resolution cosmological simulations to compare the effect of bursty star formation histories on dwarf galaxy structure for two different subgrid supernovae (SNe) feedback models in dwarf galaxies with stellar masses from $5000 <$ M$_*$/M$_\odot$ $< 10^{9}$. Our simulations are run using two distinct supernova feedback models: superbubble and blastwave. We show that both models are capable of producing galaxies that are cored and reproduce observed scaling relations for metallicity, luminosity, mass, and size. We show that continuous bursty star formation and the resulting stellar feedback are able to sustain dark matter cores in the higher dwarf galaxy mass regime, while the majority of ultra-faint and classical dwarfs retain cuspy central dark matter density profiles. We find that both subgrid SN models are able to create bursty star formation histories. We find that effective core formation peaks at M$_*$/M$_\odot$ $\simeq 5 \times 10^{-3}$ for both feedback models. Galaxies simulated with superbubble feedback peak at lower mean burstiness values relative to blastwave feedback, indicating that core formation in the superbubble sample may be less motivated by the burstiness of star formation. △ Less

Submitted 11 January, 2024; originally announced January 2024.

Comments: 18 pages, 7 figures

Precise Coil Alignment for Dynamic Wireless Charging of Electric Vehicles with RFID Sensing

Authors: Haijian Sun, Xiang Ma, Rose Qingyang Hu, Randy Christensen

Abstract: Electric vehicle (EV) has emerged as a transformative force for the sustainable and environmentally friendly future. To alleviate range anxiety caused by battery and charging facility, dynamic wireless power transfer (DWPT) is increasingly recognized as a key enabler for widespread EV adoption, yet it faces significant technical challenges, primarily in precise coil alignment. This article begins… ▽ More Electric vehicle (EV) has emerged as a transformative force for the sustainable and environmentally friendly future. To alleviate range anxiety caused by battery and charging facility, dynamic wireless power transfer (DWPT) is increasingly recognized as a key enabler for widespread EV adoption, yet it faces significant technical challenges, primarily in precise coil alignment. This article begins by reviewing current alignment methodologies and evaluates their advantages and limitations. We observe that achieving the necessary alignment precision is challenging with these existing methods. To address this, we present an innovative RFID-based DWPT coil alignment system, utilizing coherent phase detection and a maximum likelihood estimation algorithm, capable of achieving sub-10 cm accuracy. This system's efficacy in providing both lateral and vertical misalignment estimates has been verified through laboratory and experimental tests. We also discuss potential challenges in broader system implementation and propose corresponding solutions. This research offers a viable and promising solution for enhancing DWPT efficiency. △ Less

Submitted 19 December, 2023; originally announced December 2023.

Comments: submitted to IEEE magazine for potential publication. 5 figure

Environment Matters: Predicted Differences in the Stellar Mass--Halo Mass Relation and History of Star Formation for Dwarf Galaxies

Authors: Charlotte R. Christensen, Alyson Brooks, Ferah Munshi, Claire Riggs, Jordan Van Nest, Hollis Akins, Thomas R Quinn, Lucas Chamberland

Abstract: We are entering an era in which we will be able to detect and characterize hundreds of dwarf galaxies within the Local Volume. It is already known that a strong dichotomy exists in the gas content and star formation properties of field dwarf galaxies versus satellite dwarfs of larger galaxies. In this work, we study the more subtle differences that may be detectable in galaxies as a function of di… ▽ More We are entering an era in which we will be able to detect and characterize hundreds of dwarf galaxies within the Local Volume. It is already known that a strong dichotomy exists in the gas content and star formation properties of field dwarf galaxies versus satellite dwarfs of larger galaxies. In this work, we study the more subtle differences that may be detectable in galaxies as a function of distance from a massive galaxy, such as the Milky Way. We compare smoothed particle hydrodynamic simulations of dwarf galaxies formed in a Local Volume-like environment (several Mpc away from a massive galaxy) to those formed nearer to Milky Way-mass halos. We find that the impact of environment on dwarf galaxies extends even beyond the immediate region surrounding Milky Way-mass halos. Even before being accreted as satellites, dwarf galaxies near a Milky Way-mass halo tend to have higher stellar masses for their halo mass than more isolated galaxies. Dwarf galaxies in high-density environments also tend to grow faster and form their stars earlier. We show observational predictions that demonstrate how these trends manifest in lower quenching rates, higher HI fractions, and bluer colors for more isolated dwarf galaxies. △ Less

Submitted 8 November, 2023; originally announced November 2023.

Comments: 26 pages, 14 figures, submitted to ApJ

Quantum Two-Way Communication Protocol Beyond Superdense Coding

Authors: Lorenzo Valentini, Kristian Skafte Jensen, René Bødker Christensen, Marco Chiani, Petar Popovski

Abstract: We introduce a generalization of one-way superdense coding to two-way communication protocols for transmitting classical bits by using entangled quantum pairs. The proposed protocol caters for provision of entangled pairs, which is assumed to be given in one-way superdense coding. The proposed protocol gives a 50% increase in both data rate and energy efficiency compared to the classical protocol.… ▽ More We introduce a generalization of one-way superdense coding to two-way communication protocols for transmitting classical bits by using entangled quantum pairs. The proposed protocol caters for provision of entangled pairs, which is assumed to be given in one-way superdense coding. The proposed protocol gives a 50% increase in both data rate and energy efficiency compared to the classical protocol. Even when decoherence is taken into consideration, the quantum protocol performs better as long as the decoherence time is not extremely short. △ Less

Submitted 6 September, 2023; originally announced September 2023.

Comments: 3 pages

MSC Class: 81P45; 94A40

Oxide layer formation prevents deteriorating ion migration in thermoelectric Cu$_2$Se during operation in air

Authors: Rasmus S. Christensen, Peter S. Thorup, Lasse R. Jørgensen, Martin Roelsgaard, Karl F. F. Fischer, Ann-Christin Dippel, Bo Brummerstedt Iversen

Abstract: Cu$_2$Se is a mixed ionic-electronic conductor with outstanding thermoelectric performance originally envisioned for space missions. Applications were discontinued due to material instability, where elemental Cu grows at the electrode interfaces during operation in vacuum. Here, we show that when Cu$_2$Se is operating in air, formation of an oxide surface layer suppresses Cu$^+$ migration along th… ▽ More Cu$_2$Se is a mixed ionic-electronic conductor with outstanding thermoelectric performance originally envisioned for space missions. Applications were discontinued due to material instability, where elemental Cu grows at the electrode interfaces during operation in vacuum. Here, we show that when Cu$_2$Se is operating in air, formation of an oxide surface layer suppresses Cu$^+$ migration along the current direction. In operando X-ray scattering and electrical resistivity measurements quantify Cu$^+$ migration through refinement of atomic occupancies and phase composition analysis. Cu deposition can be prevented during operation in air, irrespective of a critical voltage, if the thermal gradient is applied along the current direction. Maximum entropy electron density analysis provides experimental evidence that Cu$^+$ migration pathways under thermal and electrical gradients differ substantially from equilibrium diffusion. The study establishes new promise for inexpensive sustainable Cu$_2$Se in thermoelectric applications, and it underscores the importance of atomistic insight into materials during thermoelectric operating conditions. △ Less

Submitted 7 August, 2023; originally announced August 2023.

B-mode polarization forecasts for GreenPol

Authors: U. Fuskeland, A. Kaplan, I. K. Wehus, H. K. Eriksen, P. R. Christensen, S. von Hausegger, H. Liu, P. M. Lubin, P. R. Meinhold, P. Naselsky, H. Thommesen, A. Zonca

Abstract: We present tensor-to-scalar ratio forecasts for GreenPol, a hypothetical ground-based B-mode experiment aiming to survey the cleanest regions of the Northern Galactic Hemisphere at five frequencies between 10 and 44 GHz. Its primary science goal would be to measure large-scale CMB polarization fluctuations at multipoles $\ell \lesssim 500$, and thereby constrain the primordial tensor-to-scalar rat… ▽ More We present tensor-to-scalar ratio forecasts for GreenPol, a hypothetical ground-based B-mode experiment aiming to survey the cleanest regions of the Northern Galactic Hemisphere at five frequencies between 10 and 44 GHz. Its primary science goal would be to measure large-scale CMB polarization fluctuations at multipoles $\ell \lesssim 500$, and thereby constrain the primordial tensor-to-scalar ratio. The observations for the suggested experiment would take place at the Summit Station (72deg N, 38deg W) on Greenland, at an altitude of 3216 m above sea level. In this paper we simulate various experimental setups, and derive limits on the tensor-to-scalar ratio after CMB component separation using a Bayesian component separation implementation called Commander. When combining the proposed experiment with Planck HFI observations for constraining polarized thermal dust emission, we find a projected limit of r<0.02 at 95 % confidence for the baseline configuration. This limit is very robust with respect to a range of important experimental parameters, including sky coverage, detector weighting, foreground priors etc. Overall, GreenPol would have the possibility to provide deep CMB polarization measurements of the Northern Galactic Hemisphere at low frequencies. △ Less

Submitted 18 April, 2024; v1 submitted 22 June, 2023; originally announced June 2023.

Comments: 10 pages, 8 figures. Published in A&A

Journal ref: A&A, 684, A128 (2024)

Private Product Computation using Quantum Entanglement

Authors: René Bødker Christensen, Petar Popovski

Abstract: In this work, we show that a pair of entangled qubits can be used to compute a product privately. More precisely, two participants with a private input from a finite field can perform local operations on a shared, Bell-like quantum state, and when these qubits are later sent to a third participant, the third participant can determine the product of the inputs, but without learning more about the i… ▽ More In this work, we show that a pair of entangled qubits can be used to compute a product privately. More precisely, two participants with a private input from a finite field can perform local operations on a shared, Bell-like quantum state, and when these qubits are later sent to a third participant, the third participant can determine the product of the inputs, but without learning more about the individual inputs. We give a concrete way to realize this product computation for arbitrary finite fields of prime order. △ Less

Submitted 27 October, 2023; v1 submitted 10 May, 2023; originally announced May 2023.

Comments: 9 pages. Identical to published journal paper

MSC Class: 81P40; 68P27; 94A60

Journal ref: IEEE Trans. Quant. Eng., Vol. 4 (2023)

Reliable Quantum Communications based on Asymmetry in Distillation and Coding

Authors: Lorenzo Valentini, René Bødker Christensen, Petar Popovski, Marco Chiani

Abstract: The reliable provision of entangled qubits is an essential precondition in a variety of schemes for distributed quantum computing. This is challenged by multiple nuisances, such as errors during the transmission over quantum links, but also due to degradation of the entanglement over time due to decoherence. The latter can be seen as a constraint on the latency of the quantum protocol, which bring… ▽ More The reliable provision of entangled qubits is an essential precondition in a variety of schemes for distributed quantum computing. This is challenged by multiple nuisances, such as errors during the transmission over quantum links, but also due to degradation of the entanglement over time due to decoherence. The latter can be seen as a constraint on the latency of the quantum protocol, which brings the problem of quantum protocol design into the context of latency-reliability constraints. We address the problem through hybrid schemes that combine: (1) indirect transmission based on teleportation and distillation; (2) direct transmission, based on quantum error correction (QEC). The intuition is that, at present, the quantum hardware offers low fidelity, which demands distillation; on the other hand, low latency can be obtained by QEC techniques. It is shown that, in the proposed framework, the distillation protocol gives rise to asymmetries that can be exploited by asymmetric quantum error correcting code (QECC), which sets the basis for unique hybrid distillation and coding design. Our results show that ad-hoc asymmetric codes give, compared to conventional QEC, a performance boost and codeword size reduction both in a single link and in a quantum network scenario. △ Less

Submitted 3 June, 2024; v1 submitted 1 May, 2023; originally announced May 2023.

Comments: 29 pages, 12 figures

Journal ref: IEEE Transactions on Quantum Engineering, vol. 5, pp. 1-13, 2024

Diurnal and Seasonal Mapping of Martian Ices With EMIRS

Authors: Aurélien Stcherbinine, Christopher S. Edwards, Michael D. Smith, Michael J. Wolff, Christopher Haberle, Eman Al Tunaiji, Nathan M. Smith, Kezman Saboi, Saadat Anwar, Lucas Lange, Philip R. Christensen

Abstract: Condensation and sublimation of ices at the surface of the planet is a key part of both the Martian H$_2$O and CO$_2$ cycles, either from a seasonal or diurnal aspect. While most of the ice is located within the polar caps, surface frost is known to be formed during nighttime down to equatorial latitudes. Here, we use data from the Emirates Mars Infrared Spectrometer (EMIRS) onboard the Emirates M… ▽ More Condensation and sublimation of ices at the surface of the planet is a key part of both the Martian H$_2$O and CO$_2$ cycles, either from a seasonal or diurnal aspect. While most of the ice is located within the polar caps, surface frost is known to be formed during nighttime down to equatorial latitudes. Here, we use data from the Emirates Mars Infrared Spectrometer (EMIRS) onboard the Emirates Mars Mission (EMM) to monitor the diurnal and seasonal evolution of the ices at the surface of Mars over almost one Martian year. The unique local time coverage provided by the instrument allows us to observe the apparition of equatorial CO$_2$ frost in the second half of the Martian night around the equinoxes, to its sublimation at sunrise. △ Less

Submitted 27 June, 2023; v1 submitted 14 March, 2023; originally announced March 2023.

Journal ref: Geophysical Research Letters, 50, e2023GL103629

Multilevel Skeletonization Using Local Separators

Authors: J. Andreas Bærentzen, Rasmus Emil Christensen, Emil Toftegaard Gæde, Eva Rotenberg

Abstract: In this paper we give a new, efficient algorithm for computing curve skeletons, based on local separators. Our efficiency stems from a multilevel approach, where we solve small problems across levels of detail and combine these in order to quickly obtain a skeleton. We do this in a highly modular fashion, ensuring complete flexibility in adapting the algorithm for specific types of input or for ot… ▽ More In this paper we give a new, efficient algorithm for computing curve skeletons, based on local separators. Our efficiency stems from a multilevel approach, where we solve small problems across levels of detail and combine these in order to quickly obtain a skeleton. We do this in a highly modular fashion, ensuring complete flexibility in adapting the algorithm for specific types of input or for otherwise targeting specific applications. Separator based skeletonization was first proposed by Bærentzen and Rotenberg in [ACM Tran. Graphics'21], showing high quality output at the cost of running times which become prohibitive for large inputs. Our new approach retains the high quality output, and applicability to any spatially embedded graph, while being orders of magnitude faster for all practical purposes. We test our skeletonization algorithm for efficiency and quality in practice, comparing it to local separator skeletonization on the University of Groningen Skeletonization Benchmark [Telea'16]. △ Less

Submitted 13 March, 2023; originally announced March 2023.

Fillet-based RRT*: A Rapid Convergence Implementation of RRT* for Curvature Constrained Vehicles

Authors: James Swedeen, Greg Droge, Randall Christensen

Abstract: Rapidly exploring random trees (RRTs) have proven effective in quickly finding feasible solutions to complex motion planning problems. RRT* is an extension of the RRT algorithm that provides probabilistic asymptotic optimality guarantees when using straight-line motion primitives. This work provides extensions to RRT and RRT* that employ fillets as motion primitives, allowing path curvature constr… ▽ More Rapidly exploring random trees (RRTs) have proven effective in quickly finding feasible solutions to complex motion planning problems. RRT* is an extension of the RRT algorithm that provides probabilistic asymptotic optimality guarantees when using straight-line motion primitives. This work provides extensions to RRT and RRT* that employ fillets as motion primitives, allowing path curvature constraints to be considered when planning. Two fillets are developed, an arc-based fillet that uses circular arcs to generate paths that respect maximum curvature constraints and a spline-based fillet that uses Bezier curves to additionally respect curvature continuity requirements. Planning with these fillets is shown to far exceed the performance of RRT* using Dubin's path motion primitives, approaching the performance of planning with straight-line path primitives. Path sampling heuristics are also introduced to accelerate convergence for nonholonomic motion planning. Comparisons to established RRT* approaches are made using the Open Motion Planning Library (OMPL). △ Less

Submitted 22 February, 2023; originally announced February 2023.

Comments: 53 pages, 23 figures, submitted to Journal of Intelligent and Robotic Systems

Transforming First-Year Calculus Teaching for Engineering Students -- Blocks with Field Specific Examples, Problems, and Exams

Authors: René Bødker Christensen, Bettina Dahl, Lisbeth Fajstrup

Abstract: Contribution: We demonstrate that it is feasible to include field specific problems in introductory mathematics courses to motivate engineering students. This is done in a way that still allows large parts of the course to be common to all students, ensuring economic viability. Background: Many first-year engineering students perceive mathematics courses as being too abstract and far from their… ▽ More Contribution: We demonstrate that it is feasible to include field specific problems in introductory mathematics courses to motivate engineering students. This is done in a way that still allows large parts of the course to be common to all students, ensuring economic viability. Background: Many first-year engineering students perceive mathematics courses as being too abstract and far from their chosen study programme. This may lead to a lack of motivation and effort, thus decreasing the learning outcomes. Intended outcomes: That engineering students recognize that the calculus and linear algebra courses are relevant for their future work within their specific field. This is intended to improve their learning outcomes. Application design: The courses have been restructured into smaller subunits, each of which has a corresponding workshop treating a real-world problem from the specific field of a given group of students. These workshops are developed in consultation with the relevant fields of study, and they are given a prominent role in the course for instance by forming the basis for the oral exams. Findings: Based on the feedback from students, we find that inclusion of field specific problems does help to highlight the applicability and importance of mathematics in engineering. When implementing such a solution, however, there are a number of challenges to keep in mind. △ Less

Submitted 12 February, 2023; originally announced February 2023.

Comments: 6 pages

MSC Class: 97D40; 97H60; 97I99

Saturn's Magnetic Field at Unprecedented Detail Achieved by Cassini's Close Encounters

Authors: Hao Cao, Michele K. Dougherty, Gregory J. Hunt, Emma J. Bunce, Ulrich R. Christensen, Krishan K. Khurana, Margaret G. Kivelson

Abstract: The last 22.5 orbits of the Cassini mission brought the spacecraft to less than 3000 km from Saturn's 1-bar surface. These close encounters offered an unprecedented view of Saturn's magnetic field, including contributions from the internal dynamo, the ionosphere, and the magnetosphere. In this chapter, we highlight the new picture of Saturn's magnetic field from the Cassini mission including the p… ▽ More The last 22.5 orbits of the Cassini mission brought the spacecraft to less than 3000 km from Saturn's 1-bar surface. These close encounters offered an unprecedented view of Saturn's magnetic field, including contributions from the internal dynamo, the ionosphere, and the magnetosphere. In this chapter, we highlight the new picture of Saturn's magnetic field from the Cassini mission including the persistent yet time-varying low-latitude field-aligned currents, Alfvén waves planet-ward of the D-ring, extreme axisymmetry, and high-degree magnetic moments. We then discuss the implications and new questions raised for Saturn's innermost magnetosphere, equatorial ionosphere, and interior. We conclude this chapter with an outlook for the future exploration of Saturn and other giant planets. △ Less

Submitted 6 January, 2023; originally announced January 2023.

Comments: Peer reviewed and accepted for publication as Chapter 5 of a multi-volume work edited by Kevin Baines, Michael Flasar, Norbert Krupp, and Thomas Stallard, entitled Cassini at Saturn: The Grand Finale, to be published by Cambridge University Press. 15 pages, 10 figures

Direct driving of simulated planetary jets by upscale energy transfer

Authors: Vincent G. A. Böning, Paula Wulff, Wieland Dietrich, Johannes Wicht, Ulrich R. Christensen

Abstract: The precise mechanism that forms jets and large-scale vortices on the giant planets is unknown. An inverse cascade has been suggested. Alternatively, energy may be directly injected by small-scale convection. Our aim is to clarify whether an inverse cascade feeds zonal jets and large-scale eddies in a system of rapidly rotating, deep, geostrophic spherical-shell convection. We analyze the nonlinea… ▽ More The precise mechanism that forms jets and large-scale vortices on the giant planets is unknown. An inverse cascade has been suggested. Alternatively, energy may be directly injected by small-scale convection. Our aim is to clarify whether an inverse cascade feeds zonal jets and large-scale eddies in a system of rapidly rotating, deep, geostrophic spherical-shell convection. We analyze the nonlinear scale-to-scale transfer of kinetic energy in such simulations as a function of the azimuthal wave number, m. We find that the main driving of the jets is associated with upscale transfer directly from the small convective scales to the jets. This transfer is very nonlocal in spectral space, bypassing large-scale structures. The jet formation is thus not driven by an inverse cascade. Instead, it is due to a direct driving by Reynolds stresses from small-scale convective flows. Initial correlations are caused by the effect of uniform background rotation and shell geometry on the flows. While the jet growth suppresses convection, it increases the correlation of the convective flows, which further amplifies the jet growth until it is balanced by viscous dissipation. To a much smaller extent, energy is transferred upscale to large-scale vortices directly from the convective scales, mostly outside the tangent cylinder. There, large-scale vortices are not driven by an inverse cascade either. Inside the tangent cylinder, the transfer to large-scale vortices is weaker, but more local in spectral space, leaving open the possibility of an inverse cascade as a driver of large-scale vortices. In addition, large-scale vortices receive kinetic energy from the jets via forward transfer. We therefore suggest a jet instability as an alternative formation mechanism of largescale vortices. Finally, we find that the jet kinetic energy scales as $\ell^{-5}$, the same as for the zonostrophic regime. △ Less

Submitted 19 December, 2022; originally announced December 2022.

Comments: 15 pages, 14 figures, accepted for publication in A&A

Journal ref: A&A 670, A15 (2023)

Microscopic theory of cavity-enhanced interactions of dipolaritons

Authors: Esben R. Christensen, Arturo Camacho-Guardian, Ovidiu Cotlet, Atac Imamoglu, Michiel Wouters, Georg M. Bruun, Iacopo Carusotto

Abstract: We develop a theoretical model for the interaction strength between a pair of exciton-polaritons with finite and parallel static electric-dipole moments. In stark contrast to simple perturbative arguments predicting the interaction to be reduced due to the hybridisation with cavity photons, we find that the strong coupling to a cavity mode leads to a drastic enhancement of the polariton interactio… ▽ More We develop a theoretical model for the interaction strength between a pair of exciton-polaritons with finite and parallel static electric-dipole moments. In stark contrast to simple perturbative arguments predicting the interaction to be reduced due to the hybridisation with cavity photons, we find that the strong coupling to a cavity mode leads to a drastic enhancement of the polariton interactions as compared to bare excitons. This enhancement is caused by the small polariton mass due to its photonic component, which increases the kinetic cost associated with minimising the dipole-dipole repulsion. Microscopic numerical calculations for one-dimensional geometries show that the interaction strength for the lower polariton branch is maximal when the exciton and the cavity mode are close to resonance, and can easily exceed state-of-the-art cavity decay rates for realistic confinement lengths, thus entering the polariton blockade regime. Analytical extension of our results to two-dimensional polariton systems highlight dipolar polaritons as a most promising solid-state platform to realize strongly interacting photonic systems. △ Less

Submitted 5 December, 2022; originally announced December 2022.

Comments: Main manuscript: 8 pages, 5 figures. Appendix: 4 pages, 3 figures

Private Randomness Agreement and its Application in Quantum Key Distribution Networks

Authors: René Bødker Christensen, Petar Popovski

Abstract: We define a variation on the well-known problem of private message transmission. This new problem called private randomness agreement (PRA) gives two participants access to a public, authenticated channel alongside the main channels, and the 'message' is not fixed a priori. Instead, the participants aim to agree on a random string completely unknown to a computationally unbounded adversary. We… ▽ More We define a variation on the well-known problem of private message transmission. This new problem called private randomness agreement (PRA) gives two participants access to a public, authenticated channel alongside the main channels, and the 'message' is not fixed a priori. Instead, the participants aim to agree on a random string completely unknown to a computationally unbounded adversary. We define privacy and reliability, and show that PRA cannot be solved in a single round. We then show that it can be solved in three rounds, albeit with exponential cost, and give an efficient four-round protocol based on polynomial evaluation. △ Less

Submitted 14 February, 2023; v1 submitted 11 October, 2022; originally announced October 2022.

Comments: 6 pages

MSC Class: 94A60; 94A64; 81P94

Journal ref: IEEE Communications Letters; vol. 27, no. 2, February 2023. pp. 477-481

Migrating Thermal Tides in the Martian Atmosphere during Aphelion Season Observed by EMM/EMIRS

Authors: Siteng Fan, François Forget, Michael D. Smith, Sandrine Guerlet, Khalid M. Badri, Samuel A. Atwood, Roland M. B. Young, Christopher S. Edwards, Philip R. Christensen, Justin Deighan, Hessa R. Al Matroushi, Antoine Bierjon, Jiandong Liu, Ehouarn Millour

Abstract: Temperature profiles retrieved using the first set of data of the Emirates Mars InfraRed Spectrometer (EMIRS) obtained during the science phase of the Emirates Mars Mission (EMM) are used for the analysis of migrating thermal tides in the Martian atmosphere. The selected data cover a solar longitude (LS) range of 60°-90° of Martian Year (MY) 36. The novel orbit design of the Hope Probe leads to a… ▽ More Temperature profiles retrieved using the first set of data of the Emirates Mars InfraRed Spectrometer (EMIRS) obtained during the science phase of the Emirates Mars Mission (EMM) are used for the analysis of migrating thermal tides in the Martian atmosphere. The selected data cover a solar longitude (LS) range of 60°-90° of Martian Year (MY) 36. The novel orbit design of the Hope Probe leads to a good geographic and local time coverage that significantly improves the analysis. Wave mode decomposition suggests dominant diurnal tide and important semi-diurnal tide with maximal amplitudes of 6K and 2K, respectively, as well as the existence of ~0.5K ter-diurnal tide. The results agree well with predictions by the Mars Planetary Climate Model (PCM), but the observed diurnal tide has an earlier phase (3h), and the semi-diurnal tide has an unexpectedly large wavelength (~200km). △ Less

Submitted 4 September, 2022; originally announced September 2022.

Comments: 11 pages, 5 figures, accepted by GRL special issue "The First Results from the Emirates Mars Mission (EMM)"

Smartwatch-based ecological momentary assessments for occupant wellness and privacy in buildings

Authors: Clayton Miller, Renee Christensen, Jin Kai Leong, Mahmoud Abdelrahman, Federico Tartarini, Matias Quintana, Andre Matthias Müller, Mario Frei

Abstract: This paper describes the adaptation of an open-source ecological momentary assessment smart-watch platform with three sets of micro-survey wellness-related questions focused on i) infectious disease (COVID-19) risk perception, ii) privacy and distraction in an office context, and iii) triggers of various movement-related behaviors in buildings. This platform was previously used to collect data for… ▽ More This paper describes the adaptation of an open-source ecological momentary assessment smart-watch platform with three sets of micro-survey wellness-related questions focused on i) infectious disease (COVID-19) risk perception, ii) privacy and distraction in an office context, and iii) triggers of various movement-related behaviors in buildings. This platform was previously used to collect data for thermal comfort, and this work extends its use to other domains. Several research participants took part in a proof-of-concept experiment by wearing a smartwatch to collect their micro-survey question preferences and perception responses for two of the question sets. Participants were also asked to install an indoor localization app on their phone to detect where precisely in the building they completed the survey. The experiment identified occupant information such as the tendencies for the research participants to prefer privacy in certain spaces and the difference between infectious disease risk perception in naturally versus mechanically ventilated spaces. △ Less

Submitted 11 August, 2022; originally announced August 2022.

Journal ref: 17th International Conference on Indoor Air Quality and Climate, INDOOR AIR 2022

A Semi-automatic Cell Tracking Process Towards Completing the 4D Atlas of C. elegans Development

Authors: Andrew Lauziere, Ryan Christensen, Hari Shroff

Abstract: The nematode Caenorhabditis elegans (C. elegans) is used as a model organism to better understand developmental biology and neurobiology. C. elegans features an invariant cell lineage, which has been catalogued and observed using fluorescence microscopy images. However, established methods to track cells in late-stage development fail to generalize once sporadic muscular twitching has begun. We bu… ▽ More The nematode Caenorhabditis elegans (C. elegans) is used as a model organism to better understand developmental biology and neurobiology. C. elegans features an invariant cell lineage, which has been catalogued and observed using fluorescence microscopy images. However, established methods to track cells in late-stage development fail to generalize once sporadic muscular twitching has begun. We build upon methodology which uses skin cells as fiducial markers to carry out cell tracking despite random twitching. In particular, we present a cell nucleus segmentation and tracking procedure which was integrated into a 3D rendering GUI to improve efficiency in tracking cells across late-stage development. Results on images depicting aforementioned muscle cell nuclei across three test embryos suggest the fiducial markers in conjunction with a classic tracking paradigm overcome sporadic twitching. △ Less

Submitted 2 August, 2022; v1 submitted 27 July, 2022; originally announced July 2022.

Path Planning with Uncertainty for Aircraft Under Threat of Detection from Ground-Based Radar

Authors: Austin Costley, Greg Droge, Randall Christensen, Robert C. Leishman, James Swedeen

Abstract: Mission planners for aircraft operating under threat of detection by ground-based radar systems are concerned with the probability of detection. Current path planning methods for such scenarios consider the aircraft pose, radar position, and radar parameters to be deterministic and known. This paper presents a framework for incorporating uncertainty in these quantities into a radar detection model… ▽ More Mission planners for aircraft operating under threat of detection by ground-based radar systems are concerned with the probability of detection. Current path planning methods for such scenarios consider the aircraft pose, radar position, and radar parameters to be deterministic and known. This paper presents a framework for incorporating uncertainty in these quantities into a radar detection model that is used by a path planner. The developed path planner evaluates the radar detection risk in the presence of uncertainties and uses linear covariance analysis to efficiently generate error budgets. The error budgets identify the contribution of each source of uncertainty (e.g., sensor measurement noise, radar position uncertainty) to the overall variability in the probability of detection. The framework is applied to a modified visibility graph path planner that uses the detection risk and its variability to calculate path adjustments, which maintain the detection risk below a specified threshold. The results show that the framework is effective at providing actionable information to the mission planner that improves the final planned path and reduces the detection risk. △ Less

Submitted 8 July, 2022; originally announced July 2022.

Three-body continuum states and Efimov physics in non-integer geometry

Authors: E. Garrido, E. R. Christensen, A. S. Jensen

Abstract: Continuum structures of three short-range interacting particles in a deformed external one-body field are investigated. We use the equivalent $d$-method employing non-integer dimension, $d$, in a spherical calculation with a dimension-dependent angular momentum barrier. We focus on dimensions close to the critical dimension, $d=d_E$, between two and three, defined by zero two-body energies, where… ▽ More Continuum structures of three short-range interacting particles in a deformed external one-body field are investigated. We use the equivalent $d$-method employing non-integer dimension, $d$, in a spherical calculation with a dimension-dependent angular momentum barrier. We focus on dimensions close to the critical dimension, $d=d_E$, between two and three, defined by zero two-body energies, where the Efimov effect can occur. We design for this dimension region a schematic, long-distance realistic, square-well based, three-body spherical model, which is used to derive analytic expressions for the wave functions, scattering lengths, phase shifts, and elastic scattering cross sections. The procedure and the results are universal, valid for all short-range potentials, and for large scattering lengths. We discuss the properties and validity of the derived expressions by means of the simplest system of three identical bosons. The derived expressions are particularly useful for very small energies, where full numerical calculations are often not feasible. For energies where the numerical calculations can be performed, a good agreement with the analytic results is found. These model results may be tested by scattering experiments for three particles in an equivalent external deformed oscillator potential. The cross sections all vanish in the zero-energy limit for $d<3$ with definite $d$-dependent power of energy. △ Less

Submitted 27 June, 2022; originally announced June 2022.

Comments: To be published in Physical Review A

Semi-Private Computation of Data Similarity with Applications to Data Valuation and Pricing

Authors: René Bødker Christensen, Shashi Raj Pandey, Petar Popovski

Abstract: Consider two data providers that want to contribute data to a certain learning model. Recent works have shown that the value of the data of one of the providers is dependent on the similarity with the data owned by the other provider. It would thus be beneficial if the two providers can calculate the similarity of their data, while keeping the actual data private. In this work, we devise multipart… ▽ More Consider two data providers that want to contribute data to a certain learning model. Recent works have shown that the value of the data of one of the providers is dependent on the similarity with the data owned by the other provider. It would thus be beneficial if the two providers can calculate the similarity of their data, while keeping the actual data private. In this work, we devise multiparty computation-protocols to compute similarity of two data sets based on correlation, while offering controllable privacy guarantees. We consider a simple model with two participating providers and develop methods to compute exact and approximate correlation, respectively, with controlled information leakage. Both protocols have computational and communication complexities that are linear in the number of data samples. We also provide general bounds on the maximal error in the approximation case, and analyse the resulting errors for practical parameter choices. △ Less

Submitted 11 April, 2023; v1 submitted 14 June, 2022; originally announced June 2022.

Comments: 11 pages

MSC Class: 94A60; 68P27

Journal ref: IEEE Transactions on Information Forensics and Security (2023). Vol 18, pp. 1978-1988

Sensitivity of Single-Pulse Radar Detection to Radar State Uncertainty

Authors: Austin Costley, Randall Christensen, Robert C. Leishman, Greg Droge

Abstract: Mission planners for aircraft operating under threat of detection from ground-based radar systems are often concerned with the probability of detection. Current approaches to path planning in such environments consider the radar state (i.e. radar position and parameters) to be deterministic and known. In practice, there is uncertainty in the radar state which induces uncertainty in the probability… ▽ More Mission planners for aircraft operating under threat of detection from ground-based radar systems are often concerned with the probability of detection. Current approaches to path planning in such environments consider the radar state (i.e. radar position and parameters) to be deterministic and known. In practice, there is uncertainty in the radar state which induces uncertainty in the probability of detection. This paper presents a method to incorporate the uncertainty of the radar state in a single-pulse radar detection model. The method linearizes the radar detection model with respect to the the radar state and uses the linearized models to estimate, to the first order, the variance of the probability of detection. The results in this paper validate the linearization using Monte Carlo analysis and illustrate the sensitivity of the probability of detection to radar state uncertainty. △ Less

Submitted 21 March, 2022; originally announced March 2022.

Two-body continuum states in non-integer geometry

Authors: Esben Rohan Christensen, Eduardo Garrido, Aksel Stenholm Jensen

Abstract: Wave functions, phase shifts and corresponding elastic cross sections are investigated for two short-range interacting particles in a deformed external oscillator field. For this we use the equivalent $d$-method employing a non-integer dimension $d$. Using a square-well potential, we derive analytic expressions for scattering lengths and phase shifts. In particular, we consider the dimension,… ▽ More Wave functions, phase shifts and corresponding elastic cross sections are investigated for two short-range interacting particles in a deformed external oscillator field. For this we use the equivalent $d$-method employing a non-integer dimension $d$. Using a square-well potential, we derive analytic expressions for scattering lengths and phase shifts. In particular, we consider the dimension, $d_E$, for infinite scattering length, where the Efimov effect occurs by addition of a third particle. We give explicitly the equivalent continuum wave functions in $d$ and ordinary three dimensional (3D) space, and show that the phase shifts are the same in both methods. Consequently the $d$-method can be used to obtain low-energy two-body elastic cross sections in an external field. △ Less

Submitted 2 March, 2022; originally announced March 2022.

Comments: 12 pages, 8 figures. Submitted to PRA for publication

Journal ref: Phys. Rev. A 105, 033308 (2022)

Sensitivity of Single-Pulse Radar Detection to Aircraft Pose Uncertainties

Authors: Austin Costley, Randall Christensen, Greg Droge, Robert C. Leishman

Abstract: Mission planners for aircraft that operate in radar detection environments are often concerned the probability of detection. The probability of detection is a nonlinear function of the aircraft pose and radar position. Current path planning techniques for this application assume that the aircraft pose is deterministic. In practice, however, the aircraft pose is estimated using a navigation filter… ▽ More Mission planners for aircraft that operate in radar detection environments are often concerned the probability of detection. The probability of detection is a nonlinear function of the aircraft pose and radar position. Current path planning techniques for this application assume that the aircraft pose is deterministic. In practice, however, the aircraft pose is estimated using a navigation filter and therefore contains uncertainty. The uncertainty in the aircraft pose induces uncertainty in the probability of detection, but this phenomenon is generally not considered when path planning. This paper provides a method for combining aircraft pose uncertainty with single-pulse radar detection models to aid mission planning efforts. The method linearizes the expression for the probability of detection and three radar cross section models. The linearized models are then used to determine the variability of the probability of detection induced by uncertainty in the aircraft pose. The results of this paper validate the linearization using Monte Carlo analysis and explore the sensitivity of the probability of detection to aircraft pose uncertainty. △ Less

Submitted 17 January, 2022; originally announced January 2022.

A mobile ion in a Fermi sea

Authors: Esben Rohan Christensen, Arturo Camacho-Guardian, Georg M. Bruun

Abstract: The remarkable single particle control of individual ions combined with the versatility of ultracold atomic gases makes hybrid ion-atom system an exciting new platform for quantum simulation of few- and many-body quantum physics. Here, we study theoretically the properties of a mobile ion immersed in a quantum degenerate gas of fermionic atoms. Using an effective low-energy atom-ion interaction to… ▽ More The remarkable single particle control of individual ions combined with the versatility of ultracold atomic gases makes hybrid ion-atom system an exciting new platform for quantum simulation of few- and many-body quantum physics. Here, we study theoretically the properties of a mobile ion immersed in a quantum degenerate gas of fermionic atoms. Using an effective low-energy atom-ion interaction together with a well established approach that includes exactly two-body correlations, we calculate the full spectral response of the ion and demonstrate the existence of several quasiparticle branches, which are charged analogues of the Fermi polaron observed in neutral atomic gases. Due to the long-range nature of the atom-ion interaction, these ionic Fermi polarons have several properties distinct from their neutral counterparts such as the simultaneous presence of several stable states and smooth transitions from repulsive to attractive polarons with increasing interaction strength. Surprisingly, the residue of the ionic polaron is shown to increase with the Fermi density for fixed interaction strength, which is in marked contrast to the neutral polaron. The properties of the ionic polaron approach that of the neutral polaron only in the low density limit where the average interparticle spacing is larger than the characteristic length of the atom-ion interaction. We finally analyse the effects of the Fermi gas on the molecular ions, which are bound atom-dimer states. △ Less

Submitted 29 October, 2021; originally announced October 2021.

Comments: Manuscript: 5 pages, 4 figures. Supplemental: 1 page

Sensitivity Study of Fiducial-Aided Navigation of Unmanned Aerial Vehicles

Authors: Amanda J. Strate, Randall Christensen

Abstract: The possible applications and benefits of autonomous Unmanned Aerial Vehicle (UAV) use in urban areas are gaining considerable attention. Before these possibilities can be realized, it is essential that UAVs be able to navigate reliably and precisely in urban environments. The most common means of determining the location of a UAV is to utilize position measurements from Global Navigation Satellit… ▽ More The possible applications and benefits of autonomous Unmanned Aerial Vehicle (UAV) use in urban areas are gaining considerable attention. Before these possibilities can be realized, it is essential that UAVs be able to navigate reliably and precisely in urban environments. The most common means of determining the location of a UAV is to utilize position measurements from Global Navigation Satellite Systems (GNSS). In urban environments, however, GNSS measurements are significantly degraded due to occlusions and multipath. This research analyzes the use of camera Line-of-Sight (LOS) measurements to self-describing fiducials as a replacement for conventional GNSS measurements. An extended Kalman filter (EKF) is developed and validated for the purpose of combining continuous measurements from an Inertial Measurement Unit (IMU) with the discrete LOS measurements to accurately estimate the states of a UAV. The sensitivity of the estimation error covariance to various system parameters is assessed, including IMU grade, fiducial placement, vehicle altitude, and image processing frequency. △ Less

Submitted 4 October, 2021; originally announced October 2021.

Comments: 9 pages, 13 figures

Analytical Aircraft State and IMU Signal Generator from Smoothed Reference Trajectory

Authors: Austin Costley, Randall Christensen, Robert C. Leishman, Greg Droge

Abstract: This work presents a method for generating position, attitude, and velocity states for an aircraft following a smoothed reference trajectory. The method also generates accelerometer and gyro measurements consistent with the aircraft states. This work describes three corner smoothing algorithms for generating a smoothed reference trajectory from a series waypoints that accounts for limitations in t… ▽ More This work presents a method for generating position, attitude, and velocity states for an aircraft following a smoothed reference trajectory. The method also generates accelerometer and gyro measurements consistent with the aircraft states. This work describes three corner smoothing algorithms for generating a smoothed reference trajectory from a series waypoints that accounts for limitations in the physical system. The smoothed reference trajectory, and curvilinear motion theory are used to generate the lower-order states. A coordinated turn maneuver is then applied to generate accelerometer and gyro measurement estimates. △ Less

Submitted 25 May, 2021; originally announced May 2021.

A Closed-Loop Linear Covariance Framework for Vehicle Path Planning in a Static Uncertain Obstacle Fiel

Authors: Randall Christensen, Greg Droge, Robert Leishman

Abstract: Path planning in an uncertain environment is a key enabler of true vehicle autonomy. Over the past two decades, numerous approaches have been developed to account for errors in the vehicle path while navigating complex and often uncertain environments. An important capability of such planning is the prediction of vehicle dispersion covariances about a candidate path. This work develops a new close… ▽ More Path planning in an uncertain environment is a key enabler of true vehicle autonomy. Over the past two decades, numerous approaches have been developed to account for errors in the vehicle path while navigating complex and often uncertain environments. An important capability of such planning is the prediction of vehicle dispersion covariances about a candidate path. This work develops a new closed-loop linear covariance (CL-LinCov) framework applicable to a wide range of autonomous system architectures. Important features of the developed framework include the (1) separation of high-level guidance from low-level control, (2) support for output-feedback controllers with internal states, dynamics, and output, and (3) multi-use continuous sensors for navigation state propagation, guidance, and feedback control. The closed-loop nature of the framework preserves the important coupling between the system dynamics, exogenous disturbances, and the guidance, navigation, and control algorithms. The developed framework is applied to a simplified model of an unmanned aerial vehicle and validated by comparison via Monte Carlo analysis. The utility of the CL-LinCov information is illustrated by its application to path planning in a static, uncertain obstacle field via a modified version of the Rapidly Exploring Random Tree algorithm. △ Less

Submitted 5 October, 2021; v1 submitted 25 May, 2021; originally announced May 2021.

Comments: 41 pages, 21 figures

An Algorithmic Approach to Entanglement-Assisted Quantum Error-Correcting Codes from the Hermitian Curve

Authors: René Bødker Christensen, Carlos Munuera, Francisco Revson F. Pereira, Diego Ruano

Abstract: We study entanglement-assisted quantum error-correcting codes (EAQECCs) arising from classical one-point algebraic geometry codes from the Hermitian curve with respect to the Hermitian inner product. Their only unknown parameter is $c$, the number of required maximally entangled quantum states since the Hermitian dual of an AG code is unknown. In this article, we present an efficient algorithmic a… ▽ More We study entanglement-assisted quantum error-correcting codes (EAQECCs) arising from classical one-point algebraic geometry codes from the Hermitian curve with respect to the Hermitian inner product. Their only unknown parameter is $c$, the number of required maximally entangled quantum states since the Hermitian dual of an AG code is unknown. In this article, we present an efficient algorithmic approach for computing $c$ for this family of EAQECCs. As a result, this algorithm allows us to provide EAQECCs with excellent parameters over any field size. △ Less

Submitted 11 January, 2022; v1 submitted 5 May, 2021; originally announced May 2021.

Comments: 23 pages

MSC Class: 94B65; 81P70; 94B05

Cavity-immune features in the spectra of superradiant crossover laser pulses

Authors: Mikkel Tang, Stefan A. Schäffer, Asbjørn A. Jørgensen, Martin R. Henriksen, Bjarke T. R. Christensen, Jörg H. Müller, Jan W. Thomsen

Abstract: Lasing in the bad cavity regime has promising applications in precision metrology due to the reduced sensitivity to cavity noise. Here we investigate the spectral properties and phase behavior of pulsed lasing on the $^1$S$_0 - ^3$P$_1$ line of $^{88}$Sr in a mK thermal ensemble, as first described in arxiv:1903.12593. The system operates in a regime where the Doppler-broadened atomic transition l… ▽ More Lasing in the bad cavity regime has promising applications in precision metrology due to the reduced sensitivity to cavity noise. Here we investigate the spectral properties and phase behavior of pulsed lasing on the $^1$S$_0 - ^3$P$_1$ line of $^{88}$Sr in a mK thermal ensemble, as first described in arxiv:1903.12593. The system operates in a regime where the Doppler-broadened atomic transition linewidth is several times larger than the cavity linewidth. We find that by detuning the cavity resonance, the influence of the cavity noise on the peak lasing frequency can be eliminated to first order despite the system not being deep in the bad cavity regime. Experimental results are compared to a model based on a Tavis-Cummings Hamiltonian, which enables us to investigate the interplay between different thermal velocity classes as the underlying mechanism for the reduction in cavity noise. These velocity-dependent dynamics can occur in pulsed lasing and during the turn-on behavior of lasers in the superradiant crossover regime. △ Less

Submitted 19 May, 2021; v1 submitted 27 April, 2021; originally announced April 2021.

Comments: 9 pages, 8 figures

Journal ref: Phys. Rev. Research 3, 033258 (2021)

An Exact Hypergraph Matching Algorithm for Nuclear Identification in Embryonic Caenorhabditis elegans

Authors: Andrew Lauziere, Ryan Christensen, Hari Shroff, Radu Balan

Abstract: Finding an optimal correspondence between point sets is a common task in computer vision. Existing techniques assume relatively simple relationships among points and do not guarantee an optimal match. We introduce an algorithm capable of exactly solving point set matching by modeling the task as hypergraph matching. The algorithm extends the classical branch and bound paradigm to select and aggreg… ▽ More Finding an optimal correspondence between point sets is a common task in computer vision. Existing techniques assume relatively simple relationships among points and do not guarantee an optimal match. We introduce an algorithm capable of exactly solving point set matching by modeling the task as hypergraph matching. The algorithm extends the classical branch and bound paradigm to select and aggregate vertices under a proposed decomposition of the multilinear objective function. The methodology is motivated by Caenorhabditis elegans, a model organism used frequently in developmental biology and neurobiology. The embryonic C. elegans contains seam cells that can act as fiducial markers allowing the identification of other nuclei during embryo development. The proposed algorithm identifies seam cells more accurately than established point-set matching methods, while providing a framework to approach other similarly complex point set matching tasks. △ Less

Submitted 9 July, 2022; v1 submitted 20 April, 2021; originally announced April 2021.

Comments: 20 pages, 11 figures

Stellar Migration and Chemical Enrichment in the Milky Way Disc: A Hybrid Model

Authors: James W. Johnson, David H. Weinberg, Fiorenzo Vincenzo, Jonathan C. Bird, Sarah R. Loebman, Alyson M. Brooks, Thomas R. Quinn, Charlotte R. Christensen, Emily J. Griffith

Abstract: We develop a hybrid model of galactic chemical evolution that combines a multi-ring computation of chemical enrichment with a prescription for stellar migration and the vertical distribution of stellar populations informed by a cosmological hydrodynamic disc galaxy simulation. Our fiducial model adopts empirically motivated forms of the star formation law and star formation history, with a gradien… ▽ More We develop a hybrid model of galactic chemical evolution that combines a multi-ring computation of chemical enrichment with a prescription for stellar migration and the vertical distribution of stellar populations informed by a cosmological hydrodynamic disc galaxy simulation. Our fiducial model adopts empirically motivated forms of the star formation law and star formation history, with a gradient in outflow mass loading tuned to reproduce the observed metallicity gradient. With this approach, the model reproduces many of the striking qualitative features of the Milky Way disc's abundance structure: (i) the dependence of the [O/Fe]-[Fe/H] distribution on radius $R_\text{gal}$ and midplane distance $|z|$; (ii) the changing shapes of the [O/H] and [Fe/H] distributions with $R_\text{gal}$ and $|z|$; (iii) a broad distribution of [O/Fe] at sub-solar metallicity and changes in the [O/Fe] distribution with $R_\text{gal}$, $|z|$, and [Fe/H]; (iv) a tight correlation between [O/Fe] and stellar age for [O/Fe] $>$ 0.1; (v) a population of young and intermediate-age $α$-enhanced stars caused by migration-induced variability in the Type Ia supernova rate; (vi) non-monotonic age-[O/H] and age-[Fe/H] relations, with large scatter and a median age of $\sim$4 Gyr near solar metallicity. Observationally motivated models with an enhanced star formation rate $\sim$2 Gyr ago improve agreement with the observed age-[Fe/H] and age-[O/H] relations, but worsen agreement with the observed age-[O/Fe] relation. None of our models predict an [O/Fe] distribution with the distinct bimodality seen in the observations, suggesting that more dramatic evolutionary pathways are required. All code and tables used for our models are publicly available through the Versatile Integrator for Chemical Evolution (VICE; https://pypi.org/project/vice). △ Less

Submitted 17 March, 2021; originally announced March 2021.

Comments: 24 pages; 18 figures; to be submitted to MNRAS; comments welcome

Linking Zonal Winds and Gravity II: explaining the equatorially antisymmetric gravity moments of Jupiter

Authors: Wieland Dietrich, Paula Wulff, Johannes Wicht, Ulrich R. Christensen

Abstract: The recent gravity field measurements of Jupiter (Juno) and Saturn (Cassini) confirm the existence of deep zonal flows reaching to a depth of 5\% and 15\% of the respective radius. Relating the zonal wind induced density perturbations to the gravity moments has become a major tool to characterise the interior dynamics of gas giants. Previous studies differ with respect to the assumptions made on h… ▽ More The recent gravity field measurements of Jupiter (Juno) and Saturn (Cassini) confirm the existence of deep zonal flows reaching to a depth of 5\% and 15\% of the respective radius. Relating the zonal wind induced density perturbations to the gravity moments has become a major tool to characterise the interior dynamics of gas giants. Previous studies differ with respect to the assumptions made on how the wind velocity relates to density anomalies, on the functional form of its decay with depth, and on the continuity of antisymmetric winds across the equatorial plane. Most of the suggested vertical structures exhibit a rather smooth radial decay of the zonal wind, which seems at odds with the observed secular variation of the magnetic field and the prevailing geostrophy of the zonal winds. Moreover, the results relied on an artificial equatorial regularisation or ignored the equatorial discontinuity altogether. We favour an alternative structure, where the equatorially antisymmetric zonal wind in an equatorial latitude belt between $\pm 21^\circ$ remains so shallow that it does not contribute to the gravity signal. The winds at higher latitudes suffice to convincingly explain the measured gravity moments. Our results indicate that the winds are geostrophic, i.e. constant along cylinders, in the outer $3000\,$ km and decay rapidly below. The preferred wind structure is 50\% deeper than previously thought, agrees with the measured gravity moment, is compliant with the magnetic constraints and the requirement of an adiabatic atmosphere and unbiased by the treatment of the equatorial discontinuity. △ Less

Submitted 25 February, 2021; originally announced February 2021.

Charged polarons and molecules in a Bose-Einstein Condensate

Authors: Esben Rohan Christensen, Arturo Camacho-Guardian, Georg M. Bruun

Abstract: We investigate the properties a mobile ion immersed in a Bose-Einstein condensate (BEC) using different theoretical approaches. A coherent state variational ansatz predicts that the ion spectral function exhibits several branches in addition to polaronic quasiparticle states, and we employ a diagrammatic analysis of the ion-atom scattering in the BEC to identify them as arising from the binding of… ▽ More We investigate the properties a mobile ion immersed in a Bose-Einstein condensate (BEC) using different theoretical approaches. A coherent state variational ansatz predicts that the ion spectral function exhibits several branches in addition to polaronic quasiparticle states, and we employ a diagrammatic analysis of the ion-atom scattering in the BEC to identify them as arising from the binding of an increasing number of bosons to the ion. We develop a simplified model describing the formation of these molecular ions showing that their spectral weight scales with the number of bound atoms. The number of atoms in the dressing cloud around the ion are calculated from thermodynamic arguments, and we finally show that the dynamics ensuing the injection of an ion into the BEC exhibits various regimes governed by coherent quasiparticle propagation and decay. △ Less

Submitted 22 December, 2020; v1 submitted 21 December, 2020; originally announced December 2020.

Comments: Manuscript: 4 pages, 3 figures. Supplemental: 3 pages, 2 figures

Journal ref: Phys. Rev. Lett. 126, 243001 (2021)

GPS-Denied Navigation Using SAR Images and Neural Networks

Authors: Teresa White, Jesse Wheeler, Colton Lindstrom, Randall Christensen, Kevin R. Moon

Abstract: Unmanned aerial vehicles (UAV) often rely on GPS for navigation. GPS signals, however, are very low in power and easily jammed or otherwise disrupted. This paper presents a method for determining the navigation errors present at the beginning of a GPS-denied period utilizing data from a synthetic aperture radar (SAR) system. This is accomplished by comparing an online-generated SAR image with a re… ▽ More Unmanned aerial vehicles (UAV) often rely on GPS for navigation. GPS signals, however, are very low in power and easily jammed or otherwise disrupted. This paper presents a method for determining the navigation errors present at the beginning of a GPS-denied period utilizing data from a synthetic aperture radar (SAR) system. This is accomplished by comparing an online-generated SAR image with a reference image obtained a priori. The distortions relative to the reference image are learned and exploited with a convolutional neural network to recover the initial navigational errors, which can be used to recover the true flight trajectory throughout the synthetic aperture. The proposed neural network approach is able to learn to predict the initial errors on both simulated and real SAR image data. △ Less

Submitted 22 October, 2020; originally announced October 2020.

Comments: 5 pages, 5 figures

Sensitivity of BPA SAR Image Formation to Initial Position, Velocity, and Attitude Navigation Errors

Authors: Colton Lindstrom, Randall Christensen, Jacob Gunther

Abstract: The Back-Projection Algorithm (BPA) is a time domain matched filtering technique to form synthetic aperture radar (SAR) images. To produce high quality BPA images, precise navigation data for the radar platform must be known. Any error in position, velocity, or attitude results in improperly formed images corrupted by shifting, blurring, and distortion. This paper develops analytical expressions t… ▽ More The Back-Projection Algorithm (BPA) is a time domain matched filtering technique to form synthetic aperture radar (SAR) images. To produce high quality BPA images, precise navigation data for the radar platform must be known. Any error in position, velocity, or attitude results in improperly formed images corrupted by shifting, blurring, and distortion. This paper develops analytical expressions that characterize the relationship between navigation errors and image formation errors. These analytical expressions are verified via simulated image formation and real data image formation. △ Less

Submitted 21 September, 2020; originally announced September 2020.

Comments: This work has been submitted to the IEEE for possible publication. Copyright may be transferred without notice, after which this version may no longer be accessible

Ultra-faint dwarfs in a Milky Way context: Introducing the Mint Condition DC Justice League Simulations

Authors: Elaad Applebaum, Alyson M. Brooks, Charlotte R. Christensen, Ferah Munshi, Thomas R. Quinn, Sijing Shen, Michael Tremmel

Abstract: We present results from the "Mint" resolution DC Justice League suite of Milky Way-like zoom-in cosmological simulations, which extend our study of nearby galaxies down into the ultra-faint dwarf (UFD) regime for the first time. The mass resolution of these simulations is the highest ever published for cosmological Milky Way zoom-in simulations run to $z=0$, with initial star (dark matter) particl… ▽ More We present results from the "Mint" resolution DC Justice League suite of Milky Way-like zoom-in cosmological simulations, which extend our study of nearby galaxies down into the ultra-faint dwarf (UFD) regime for the first time. The mass resolution of these simulations is the highest ever published for cosmological Milky Way zoom-in simulations run to $z=0$, with initial star (dark matter) particle masses of 994 (17900) M$_\odot$, and a force resolution of 87 pc. We study the surrounding dwarfs and UFDs, and find the simulations match the observed dynamical properties of galaxies with $-3 < M_V < -19$, and reproduce the scatter seen in the size-luminosity plane for $r_h\gtrsim200$ pc. We predict the vast majority of nearby galaxies will be observable by the Vera Rubin Observatory's co-added Legacy Survey of Space and Time (LSST). We additionally show that faint dwarfs with velocity dispersions $\lesssim5$ km/s result from severe tidal stripping of the host halo. We investigate the quenching of UFDs in a hydrodynamical Milky Way context, and find that the majority of UFDs are quenched prior to interactions with the Milky Way, though some of the quenched UFDs retain their gas until infall. Additionally these simulations yield some unique dwarfs that are the first of their kind to be simulated, e.g., an HI-rich field UFD, a late-forming UFD that has structural properties similar to Crater 2, as well as a compact dwarf satellite that has no dark matter at $z=0$. △ Less

Submitted 17 November, 2020; v1 submitted 25 August, 2020; originally announced August 2020.

Comments: 33 pages, 17 figures, ApJ accepted version

Quenching timescales of dwarf satellites around Milky Way-mass hosts

Authors: Hollis B. Akins, Charlotte R. Christensen, Alyson M. Brooks, Ferah Munshi, Elaad Applebaum, Anna Engelhardt, Lucas Chamberland

Abstract: Observations of the low-mass satellites in the Local Group have shown high fractions of gas-poor, quiescent galaxies relative to isolated dwarfs, implying that the host halo environment plays an important role in the quenching of dwarf galaxies. In this work, we present measurements of the quenched fractions and quenching timescales of dwarf satellite galaxies in the DC Justice League suite of 4 h… ▽ More Observations of the low-mass satellites in the Local Group have shown high fractions of gas-poor, quiescent galaxies relative to isolated dwarfs, implying that the host halo environment plays an important role in the quenching of dwarf galaxies. In this work, we present measurements of the quenched fractions and quenching timescales of dwarf satellite galaxies in the DC Justice League suite of 4 high-resolution cosmological zoom-in simulations of Milky Way-mass halos. We show that these simulations accurately reproduce the satellite luminosity functions of observed nearby galaxies, as well as the variation in satellite quenched fractions from $M_* \sim 10^5$ solar masses to $10^{10}$ solar masses. We then trace the histories of satellite galaxies back to $z \sim 15$, and find that many satellites with $M_* \sim 10^{6-8}$ solar masses quench within 2 Gyr of infall into the host halo, while others in the same mass range remain star-forming for as long as 5 Gyr. We show that this scatter can be explained by the satellite's gas mass and the ram pressure it feels at infall. Finally, we identify a characteristic stellar mass scale of $10^8$ solar masses above which infalling satellites are largely resistant to rapid environmental quenching. △ Less

Submitted 2 February, 2021; v1 submitted 6 August, 2020; originally announced August 2020.

Comments: 16 pages, 7 figures, accepted to ApJ

Zero-Error Tracking for Autonomous Vehicles through Epsilon-Trajectory Generation

Authors: Clint Ferrin, Greg Droge, Randall Christensen

Abstract: This paper presents a control method and trajectory planner for vehicles with first-order nonholonomic constraints that guarantee asymptotic convergence to a time-indexed trajectory. To overcome the nonholonomic constraint, a fixed point in front of the vehicle can be controlled to track a desired trajectory, albeit with a steady-state error. To eliminate steady state error, a sufficiently smooth… ▽ More This paper presents a control method and trajectory planner for vehicles with first-order nonholonomic constraints that guarantee asymptotic convergence to a time-indexed trajectory. To overcome the nonholonomic constraint, a fixed point in front of the vehicle can be controlled to track a desired trajectory, albeit with a steady-state error. To eliminate steady state error, a sufficiently smooth trajectory is reformulated for the new reference point such that, when tracking the new trajectory, the vehicle asymptotically converges to the original trajectory. The resulting zero-error tracking law is demonstrated through a novel framework for creating time-indexed Clothoids. The Clothoids can be planned to pass through arbitrary waypoints using traditional methods yet result in trajectories that can be followed with zero steady-state error. The results of the control method and planner are illustrated in simulation wherein zero-error tracking is demonstrated. △ Less

Submitted 20 July, 2020; originally announced July 2020.

Inside Out and Upside-Down: The Roles of Gas Cooling and Dynamical Heating in Shaping the Stellar Age-Velocity Relation

Authors: Jonathan C. Bird, Sarah R. Loebman, David H. Weinberg, Alyson Brooks, Thomas R. Quinn, Charlotte R. Christensen

Abstract: Kinematic studies of disk galaxies, using individual stars in the Milky Way or statistical studies of global disk kinematics over time, provide insight into how disks form and evolve. We use a high-resolution, cosmological zoom-simulation of a Milky Way-mass disk galaxy h277 to tie together local disk kinematics and the evolution of the disk over time. The present-day stellar age-velocity relation… ▽ More Kinematic studies of disk galaxies, using individual stars in the Milky Way or statistical studies of global disk kinematics over time, provide insight into how disks form and evolve. We use a high-resolution, cosmological zoom-simulation of a Milky Way-mass disk galaxy h277 to tie together local disk kinematics and the evolution of the disk over time. The present-day stellar age-velocity relationship (AVR) of h277 is nearly identical to that of the analogous solar-neighborhood measurement in the Milky Way. A crucial element of this success is the simulation's dynamically cold multi-phase ISM, which allows young stars to form with a low velocity dispersion ($σ_{\mathrm{birth}}$$\sim 6 - 8\ \mathrm{km}\, \mathrm{s}^{-1}$) at late times. Older stars are born kinematically hotter (i.e., the disk settles over time in an "upside-down" formation scenario), and are subsequently heated after birth. The disk also grows "inside-out", and many of the older stars in the solar neighborhood at $z=0$ are present because of radial mixing. We demonstrate that the evolution of $σ_{\mathrm{birth}}$ in h277 can be explained by the same model used to describe the general decrease in velocity dispersion observed in disk galaxies from $z\sim 2-3$ to the present-day, in which the disk evolves in quasi-stable equilibrium and the ISM velocity dispersion decreases over time due to a decreasing gas fraction. Thus, our results tie together local observations of the Milky Way's AVR with observed kinematics of high $z$ disk galaxies. △ Less

Submitted 26 May, 2020; originally announced May 2020.

Comments: Submitting to MNRAS. Comments welcomed. This simulation is publicly available. Contact Alyson Brooks for details

Statistics of CMB polarization angles

Authors: James Creswell, Nadia Dachlythra, Hao Liu, Pavel Naselsky, Per Rex Christensen

Abstract: We study the distribution functions of the CMB polarization angle $ψ$, focusing on the Planck 2018 CMB maps. We extend the model of Preece & Battye (2014) of Gaussian correlated $Q$ and $U$ Stokes parameters to allow nonzero means. When the variances of $Q$ and $U$ are equal and their covariance and means are zero, the polarization angle is uniformly distributed. Otherwise the uniform distribution… ▽ More We study the distribution functions of the CMB polarization angle $ψ$, focusing on the Planck 2018 CMB maps. We extend the model of Preece & Battye (2014) of Gaussian correlated $Q$ and $U$ Stokes parameters to allow nonzero means. When the variances of $Q$ and $U$ are equal and their covariance and means are zero, the polarization angle is uniformly distributed. Otherwise the uniform distribution will be modulated by harmonics with $2ψ$ and $4ψ$ phases. These modulations are visible in the Planck 2018 CMB maps. Furthermore, the mean value of $U$ is peculiar compared to the power spectrum. △ Less

Submitted 6 January, 2020; originally announced January 2020.

Comments: 14 pages, 6 figures

Linking Zonal Winds and Gravity: The Relative Importance of Dynamic Self Gravity

Authors: Johannes Wicht, Wieland Dietrich, Paula Wulff, Ulrich R. Christensen

Abstract: Recent precise measurements at Jupiter's and Saturn's gravity fields constrain the properties of the zonal flows in the outer envelopes of these planets. A simplified dynamic equation, sometimes called the thermal wind or thermo-gravitational wind equation, establishes a link between zonal flows and the related buoyancy perturbation, which in turn can be exploited to yield the dynamic gravity pert… ▽ More Recent precise measurements at Jupiter's and Saturn's gravity fields constrain the properties of the zonal flows in the outer envelopes of these planets. A simplified dynamic equation, sometimes called the thermal wind or thermo-gravitational wind equation, establishes a link between zonal flows and the related buoyancy perturbation, which in turn can be exploited to yield the dynamic gravity perturbation. Whether or not the action of the dynamic gravity perturbation needs to be explicitly included in this equation, an effect we call the Dynamic Self Gravity (DSG), has been a matter of intense debate. We show that, under reasonable assumptions, the equation can be solved (semi) analytically. This allows us to quantify the impact of the DSG on each gravity harmonic, practically independent of the zonal flow or the details of the planetary interior model. The impact decreases with growing spherical harmonic degree l. For degrees l=2 to about l=4, the DSG is a first order effect and should be taken into account in any attempt of inverting gravity measurements for zonal flow properties. For degrees of about l=5 to roughly l=10, the relative impact of DSG is about 10% and thus seems worthwhile to include, in particular since this comes at little extra costs with the method presented here. For yet higher degrees, is seems questionable whether gravity measurements or interior models will ever reach the required precision equivalent of the DSG impact of only a few percent of less. △ Less

Submitted 5 September, 2019; originally announced September 2019.

Comments: 28 pages, 6 figures

On Steane-Enlargement of Quantum Codes from Cartesian Product Point Sets

Authors: René Bødker Christensen, Olav Geil

Abstract: In this work, we study quantum error-correcting codes obtained by using Steane-enlargement. We apply this technique to certain codes defined from Cartesian products previously considered by Galindo et al. in [4]. We give bounds on the dimension increase obtained via enlargement, and additionally give an algorithm to compute the true increase. A number of examples of codes are provided, and their p… ▽ More In this work, we study quantum error-correcting codes obtained by using Steane-enlargement. We apply this technique to certain codes defined from Cartesian products previously considered by Galindo et al. in [4]. We give bounds on the dimension increase obtained via enlargement, and additionally give an algorithm to compute the true increase. A number of examples of codes are provided, and their parameters are compared to relevant codes in the literature, which shows that the parameters of the enlarged codes are advantageous. Furthermore, comparison with the Gilbert-Varshamov bound for stabilizer quantum codes shows that several of the enlarged codes match or exceed the parameters promised by the bound. △ Less

Submitted 13 August, 2019; originally announced August 2019.

Comments: 12 pages

MSC Class: 94B27; 81Q99

Branching ratios in the $β$ decay of $^{16}$N

Authors: O. S. Kirsebom, E. R. Christensen

Abstract: In this brief note, we present the results of an experiment performed at the ISOLDE Decay Station at CERN in which several of the branching ratios in the $β$ decay of $^{16}$N were determined with high precision and accuracy. In this brief note, we present the results of an experiment performed at the ISOLDE Decay Station at CERN in which several of the branching ratios in the $β$ decay of $^{16}$N were determined with high precision and accuracy. △ Less

Submitted 23 May, 2019; originally announced May 2019.

Comments: 4 pages, 3 tables

Steane-Enlargement of Quantum Codes from the Hermitian Curve

Authors: René Bødker Christensen, Olav Geil

Abstract: In this paper, we study the construction of quantum codes by applying Steane-enlargement to codes from the Hermitian curve. We cover Steane-enlargement of both usual one-point Hermitian codes and of order bound improved Hermitian codes. In particular, the paper contains two constructions of quantum codes whose parameters are described by explicit formulae, and we show that these codes compare favo… ▽ More In this paper, we study the construction of quantum codes by applying Steane-enlargement to codes from the Hermitian curve. We cover Steane-enlargement of both usual one-point Hermitian codes and of order bound improved Hermitian codes. In particular, the paper contains two constructions of quantum codes whose parameters are described by explicit formulae, and we show that these codes compare favourably to existing, comparable constructions in the literature. △ Less

Submitted 22 April, 2019; originally announced April 2019.

Comments: 11 pages

MSC Class: 94B27; 81Q99

Lasing on a narrow transition in a cold thermal strontium ensemble

Authors: Stefan A. Schäffer, Mikkel Tang, Martin R. Henriksen, Asbjørn A. Jørgensen, Bjarke T. R. Christensen, Jan W. Thomsen

Abstract: Highly stable laser sources based on narrow atomic transitions provide a promising platform for direct generation of stable and accurate optical frequencies. Here we investigate a simple system operating in the high-temperature regime of cold atoms. The interaction between a thermal ensemble of $^{88}$Sr at mK temperatures and a medium-finesse cavity produces strong collective coupling and facilit… ▽ More Highly stable laser sources based on narrow atomic transitions provide a promising platform for direct generation of stable and accurate optical frequencies. Here we investigate a simple system operating in the high-temperature regime of cold atoms. The interaction between a thermal ensemble of $^{88}$Sr at mK temperatures and a medium-finesse cavity produces strong collective coupling and facilitates high atomic coherence which causes lasing on the dipole forbidden $^1$S$_0 \leftrightarrow ^3$P$_1$ transition. We experimentally and theoretically characterize the lasing threshold and evolution of such a system, and investigate decoherence effects in an unconfined ensemble. We model the system using a Tavis-Cummings model, and characterize velocity-dependent dynamics of the atoms as well as the dependency on the cavity-detuning. △ Less

Submitted 19 November, 2019; v1 submitted 29 March, 2019; originally announced March 2019.

Comments: 9 pages, 7 figures

Journal ref: Phys. Rev. A 101, 013819 (2020)

Roll of the Dice: A Stochastically Sampled IMF Alters the Stellar Content of Simulated Dwarf Galaxies

Authors: Elaad Applebaum, Alyson M. Brooks, Thomas R. Quinn, Charlotte R. Christensen

Abstract: Cosmological simulations are reaching the resolution necessary to study ultra-faint dwarf galaxies. Observations indicate that in small populations, the stellar initial mass function (IMF) is not fully populated; rather, stars are sampled in a way that can be approximated as coming from an underlying probability density function. To ensure the accuracy of cosmological simulations in the ultra-fain… ▽ More Cosmological simulations are reaching the resolution necessary to study ultra-faint dwarf galaxies. Observations indicate that in small populations, the stellar initial mass function (IMF) is not fully populated; rather, stars are sampled in a way that can be approximated as coming from an underlying probability density function. To ensure the accuracy of cosmological simulations in the ultra-faint regime, we present an improved treatment of the IMF. We implement a self-consistent, stochastically populated IMF in cosmological hydrodynamic simulations. We test our method using high-resolution simulations of a Milky Way halo, run to $z=6$, yielding a sample of nearly 100 galaxies. We also use an isolated dwarf galaxy to investigate the resulting systematic differences in galaxy properties. We find that a stochastic IMF in simulations makes feedback burstier, strengthening feedback, and quenching star formation earlier in small dwarf galaxies. For galaxies in halos with mass $\lesssim10^{8.5}$ M$_\odot$, a stochastic IMF typically leads to lower stellar mass compared to a continuous IMF, sometimes by more than an order of magnitude. We show that existing methods of ensuring discrete supernovae incorrectly determine the mass of the star particle and its associated feedback. This leads to overcooling of surrounding gas, with at least ${\sim}$10 per cent higher star formation and ${\sim}$30 per cent higher cold gas content. Going forward, to accurately model dwarf galaxies and compare to observations, it will be necessary to incorporate a stochastically populated IMF that samples the full spectrum of stellar masses. △ Less

Submitted 20 February, 2020; v1 submitted 31 October, 2018; originally announced November 2018.

Comments: 14 pages, 11 figures, Updated to reflect accepted version, MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, 2020, 492, 8-21