Showing 1–50 of 119 results for author: Ludlow, A

Milky Way-mass disc galaxies with low-mass stellar haloes have diverse merger histories

Authors: Katy L. Proctor, Aaron D. Ludlow, Claudia del P. Lagos, Aaron S. G. Robotham

Abstract: We use the Eagle simulation to study the relationship between the stellar haloes of Milky-Way-mass galaxies and their recent merger histories. The stellar mass ratio of the most massive merger that occurred since $z=1$ is strongly correlated with the $z=0$ fraction of ex situ stars in the galaxy, $f_{\rm ex\ situ}$, but is weakly correlated with stellar halo mass fraction, $f_{\rm SH}$, particular… ▽ More We use the Eagle simulation to study the relationship between the stellar haloes of Milky-Way-mass galaxies and their recent merger histories. The stellar mass ratio of the most massive merger that occurred since $z=1$ is strongly correlated with the $z=0$ fraction of ex situ stars in the galaxy, $f_{\rm ex\ situ}$, but is weakly correlated with stellar halo mass fraction, $f_{\rm SH}$, particularly for disc galaxies. Contrary to common belief, our results suggest that disc galaxies with low mass stellar haloes do not necessarily have quiescent merger histories; in fact, roughly one quarter have experienced a merger with stellar mass ratio $> 0.1$ since $z=1$. We demonstrate that the population of disc galaxies with low $f_{\rm SH}$ and active merger histories undergo mergers with satellites whose orbits are more circular than average and are approximately co-planar with the disc; instead of contributing significantly to the stellar halo, these mergers lead to discs that contain a substantial fraction of ex situ stellar mass and are thicker and more extended than those of quiescent galaxies. Such mergers also supply fuel that often incites a significant episode of in situ star formation in the disc. Our results suggest that seemingly quiescent disc galaxies with low-mass stellar haloes actually have diverse merger histories, which limits the extent to which stellar haloes alone can act as observable tracers of the hierarchical assembly of galaxies. △ Less

Submitted 16 July, 2024; originally announced July 2024.

Comments: 11 pages, 8 figures. Submitted to MNRAS

Clock-line-mediated Sisyphus Cooling

Authors: Chun-Chia Chen, Jacob L. Siegel, Benjamin D. Hunt, Tanner Grogan, Youssef S. Hassan, Kyle Beloy, Kurt Gibble, Roger C. Brown, Andrew D. Ludlow

Abstract: We demonstrate sub-recoil Sisyphus cooling using the long-lived $^{3}\mathrm{P}_{0}$ clock state in alkaline-earth-like ytterbium. A 1388 nm optical standing wave nearly resonant with the $^{3}\textrm{P}_{0}$$\,\rightarrow$$\,^{3}\textrm{D}_{1}$ transition creates a spatially periodic light shift of the $^{3}\textrm{P}_{0}$ clock state. Following excitation on the ultranarrow clock transition, we… ▽ More We demonstrate sub-recoil Sisyphus cooling using the long-lived $^{3}\mathrm{P}_{0}$ clock state in alkaline-earth-like ytterbium. A 1388 nm optical standing wave nearly resonant with the $^{3}\textrm{P}_{0}$$\,\rightarrow$$\,^{3}\textrm{D}_{1}$ transition creates a spatially periodic light shift of the $^{3}\textrm{P}_{0}$ clock state. Following excitation on the ultranarrow clock transition, we observe Sisyphus cooling in this potential, as the light shift is correlated with excitation to $^{3}\textrm{D}_{1}$ and subsequent spontaneous decay to the $^{1}\textrm{S}_{0}$ ground state. We observe that cooling enhances the loading efficiency of atoms into a 759 nm magic-wavelength one-dimensional (1D) optical lattice, as compared to standard Doppler cooling on the $^{1}\textrm{S}_{0}$$\,\rightarrow\,$$^{3}\textrm{P}_{1}$ transition. Sisyphus cooling yields temperatures below 200 nK in the weakly confined, transverse dimensions of the 1D optical lattice. These lower temperatures improve optical lattice clocks by facilitating the use of shallow lattices with reduced light shifts, while retaining large atom numbers to reduce the quantum projection noise. This Sisyphus cooling can be pulsed or continuous and is applicable to a range of quantum metrology applications. △ Less

Submitted 19 June, 2024; originally announced June 2024.

Comments: 8 pages, 6 figures

Ratchet Loading and Multi-Ensemble Operation in an Optical Lattice Clock

Authors: Youssef S. Hassan, Takumi Kobayashi, Tobias Bothwell, Jacob L. Seigel, Benjamin D. Hunt, Kyle Beloy, Kurt Gibble, Tanner Grogan, Andrew D. Ludlow

Abstract: We demonstrate programmable control over the spatial distribution of ultra-cold atoms confined in an optical lattice. The control is facilitated through a combination of spatial manipulation of the magneto-optical trap and atomic population shelving to a metastable state. We first employ the technique to load an extended (5 mm) atomic sample with uniform density in an optical lattice clock, reduci… ▽ More We demonstrate programmable control over the spatial distribution of ultra-cold atoms confined in an optical lattice. The control is facilitated through a combination of spatial manipulation of the magneto-optical trap and atomic population shelving to a metastable state. We first employ the technique to load an extended (5 mm) atomic sample with uniform density in an optical lattice clock, reducing atomic interactions and realizing remarkable frequency homogeneity across the atomic cloud. We also prepare multiple spatially separated atomic ensembles and realize multi-ensemble clock operation within the standard one-dimensional (1D) optical lattice clock architecture. Leveraging this technique, we prepare two oppositely spin-polarized ensembles that are independently addressable, offering a platform for implementing spectroscopic protocols for enhanced tracking of local oscillator phase. Finally, we demonstrate a relative fractional frequency instability at one second of $2.4(1) \times10^{-17}$ between two ensembles, useful for characterization of intra-lattice differential systematics. △ Less

Submitted 30 May, 2024; originally announced May 2024.

Comments: 10 pages, 6 figures

Excited-Band Coherent Delocalization for Improved Optical Lattice Clock Performance

Authors: Jacob Siegel, William McGrew, Youssef Hassan, Chun-Chia Chen, Kyle Beloy, Tanner Grogan, Xiaogang Zhang, Andrew Ludlow

Abstract: We implement coherent delocalization as a tool for improving the two primary metrics of atomic clock performance: systematic uncertainty and instability. By decreasing atomic density with coherent delocalization, we suppress cold-collision shifts and two-body losses. Atom loss attributed to Landau-Zener tunneling in the ground lattice band would compromise coherent delocalization at low trap depth… ▽ More We implement coherent delocalization as a tool for improving the two primary metrics of atomic clock performance: systematic uncertainty and instability. By decreasing atomic density with coherent delocalization, we suppress cold-collision shifts and two-body losses. Atom loss attributed to Landau-Zener tunneling in the ground lattice band would compromise coherent delocalization at low trap depths for our $^{171}$Yb atoms; hence, we implement for the first time delocalization in excited lattice bands. Doing so increases the spatial distribution of atoms trapped in the vertically-oriented optical lattice by $\sim7$ times. At the same time we observe a reduction of the cold-collision shift by 6.5(8) times, while also making inelastic two-body loss negligible. With these advantages, we measure the trap-light-induced quenching rate and natural lifetime of the ${}^3$P${}_0$ excited-state as $5.7(7)\times10^{-4}$ $E_r^{-1}s^{-1}$ and 19(2) s, respectively. △ Less

Submitted 12 February, 2024; v1 submitted 7 February, 2024; originally announced February 2024.

The distribution of stellar orbits in Eagle galaxies -- the effect of mergers, gas accretion, and secular evolution

Authors: Giulia Santucci, Claudia Del P. Lagos, Katherine E. Harborne, Aaron Ludlow, Caro Foster, Richard McDermid, Adriano Poci, Katy L. Proctor, Sabine Thater, Glenn van de Ven, Ling Zhu, Daniel Walo Martin

Abstract: The merger history of a galaxy is thought to be one of the major factors determining its internal dynamics, with galaxies having undergone different types or mergers (e.g. dry, minor or major mergers) predicted to show different dynamical properties. We study the instantaneous orbital distribution of galaxies in the Eagle simulation, colouring the orbits of the stellar particles by their stellar a… ▽ More The merger history of a galaxy is thought to be one of the major factors determining its internal dynamics, with galaxies having undergone different types or mergers (e.g. dry, minor or major mergers) predicted to show different dynamical properties. We study the instantaneous orbital distribution of galaxies in the Eagle simulation, colouring the orbits of the stellar particles by their stellar age, in order to understand whether stars form in particular orbits (e.g. in a thin or thick disc). We first show that Eagle reproduces well the observed stellar mass fractions in different stellar orbital families as a function of stellar mass and spin parameter at z = 0. We find that the youngest stars reside in a thin disc component that can extend to the very inner regions of galaxies, and that older stars have warmer orbits, with the oldest ones showing orbits consistent with both hot and counter-rotating classifications, which is consistent with the trend found in the Milky-Way and other disc galaxies. We also show that counter-rotating orbits trace galaxy mergers - in particular dry mergers, and that in the absence of mergers, counter-rotating orbits can also be born from highly misaligned gas accretion that leads to star formation. △ Less

Submitted 10 October, 2023; originally announced October 2023.

Comments: Submitted to MNRAS, 21 pages

Identifying the discs, bulges, and intra-halo light of simulated galaxies through structural decomposition

Authors: Katy L. Proctor, Claudia del P. Lagos, Aaron D. Ludlow, Aaron S. G. Robotham

Abstract: We perform a structural decomposition of galaxies identified in three cosmological hydrodynamical simulations by applying Gaussian Mixture Models (GMMs) to the kinematics of their stellar particles. We study the resulting disc, bulge, and intra-halo light (IHL) components of galaxies whose host dark matter haloes have virial masses in the range $M_{200}=10^{11}$-- $10^{15}\,{\rm M_\odot}$. Our dec… ▽ More We perform a structural decomposition of galaxies identified in three cosmological hydrodynamical simulations by applying Gaussian Mixture Models (GMMs) to the kinematics of their stellar particles. We study the resulting disc, bulge, and intra-halo light (IHL) components of galaxies whose host dark matter haloes have virial masses in the range $M_{200}=10^{11}$-- $10^{15}\,{\rm M_\odot}$. Our decomposition technique isolates galactic discs whose mass fractions, $f_{\rm disc}$, correlate strongly with common alternative morphology indicators; for example, $f_{\rm disc}$ is approximately equal to $κ_{\rm co}$, the fraction of stellar kinetic energy in co-rotation. The primary aim of our study, however, is to characterise the IHL of galaxies in a consistent manner and over a broad mass range, and to analyse its properties from the scale of galactic stellar haloes up to the intra-cluster light. Our results imply that the IHL fraction, $f_{\rm IHL}$, has appreciable scatter and is strongly correlated with galaxy morphology: at fixed stellar mass, the IHL of disc galaxies is typically older and less massive than that of spheroids. Above $M_{200}\approx 10^{13}\,{\rm M_\odot}$, we find, on average, $f_{\rm IHL}\approx 0.45$, albeit with considerable scatter. The transition radius beyond which the IHL dominates the stellar mass of a galaxy is roughly $30\,{\rm kpc}$ for $M_{200}\lesssim 10^{12.8}\,{\rm M_\odot}$, but increases strongly towards higher masses. However, we find that no alternative IHL definitions -- whether based on the ex-situ stellar mass, or the stellar mass outside a spherical aperture -- reproduce our dynamically-defined IHL masses. △ Less

Submitted 6 December, 2023; v1 submitted 24 July, 2023; originally announced July 2023.

Comments: 16 pages, 11 figures. MNRAS published

Optimal binary gratings for multi-wavelength magneto-optical traps

Authors: Oliver S. Burrow, Robert J. Fasano, Wesley Brand, Michael W. Wright, Wenbo Li, Andrew D. Ludlow, Erling Riis, Paul F. Griffin, Aidan S. Arnold

Abstract: Grating magneto-optical traps are an enabling quantum technology for portable metrological devices with ultracold atoms. However, beam diffraction efficiency and angle are affected by wavelength, creating a single-optic design challenge for laser cooling in two stages at two distinct wavelengths - as commonly used for loading e.g. Sr or Yb atoms into optical lattice or tweezer clocks. Here, we opt… ▽ More Grating magneto-optical traps are an enabling quantum technology for portable metrological devices with ultracold atoms. However, beam diffraction efficiency and angle are affected by wavelength, creating a single-optic design challenge for laser cooling in two stages at two distinct wavelengths - as commonly used for loading e.g. Sr or Yb atoms into optical lattice or tweezer clocks. Here, we optically characterize a wide variety of binary gratings at different wavelengths to find a simple empirical fit to experimental grating diffraction efficiency data in terms of dimensionless etch depth and period for various duty cycles. The model avoids complex 3D light-grating surface calculations, yet still yields results accurate to a few percent across a broad range of parameters. Gratings optimized for two (or more) wavelengths can now be designed in an informed manner suitable for a wide class of atomic species enabling advanced quantum technologies. △ Less

Submitted 18 November, 2023; v1 submitted 29 June, 2023; originally announced June 2023.

Comments: 11 pages, 3 figures, 2 tables. See also 4-page supplementary document

Spurious heating of stellar motions by dark matter particles in cosmological simulations of galaxy formation

Authors: Aaron D. Ludlow, S. Michael Fall, Matthew J. Wilkinson, Joop Schaye, Danail Obreschkow

Abstract: We use two cosmological simulations to study the impact of spurious heating of stellar motions within simulated galaxies by dark matter (DM) particles. The simulations share the same numerical and subgrid parameters, but one used a factor of 7 more DM particles. Many galaxy properties are unaffected by spurious heating, including their masses, star formation histories, and the spatial distribution… ▽ More We use two cosmological simulations to study the impact of spurious heating of stellar motions within simulated galaxies by dark matter (DM) particles. The simulations share the same numerical and subgrid parameters, but one used a factor of 7 more DM particles. Many galaxy properties are unaffected by spurious heating, including their masses, star formation histories, and the spatial distribution of their gaseous baryons. The distribution and kinematics of stellar and DM particles, however, are affected. Below a resolution-dependent virial mass, $M_{200}^{\rm spur}$, galaxies have higher characteristic velocities, larger sizes, and more angular momentum in the simulation with lower DM mass resolution; haloes have higher central densities and lower velocity dispersions. Above $M_{200}^{\rm spur}$, galaxies and haloes have similar properties in both runs. The differences arise due to spurious heating, which transfers energy from DM to stellar particles, causing galaxies to heat up and haloes to cool down. The value of $M_{200}^{\rm spur}$ can be derived from an empirical disc heating model, and coincides with the mass below which the predicted {\em spurious} velocity dispersion exceeds the {\em measured} velocity dispersion of simulated galaxies. We predict that galaxies in the $100^3\, {\rm Mpc}^3$ \eagle\, run and IllustrisTNG-100 are robust to spurious collisional effects at their half-mass radii provided $M_{200}^{\rm spur}\approx 10^{11.7}{\rm M_\odot}$; for the $25^3\, {\rm Mpc}^3$ \eagle\, run and IllustrisTNG-50, we predict $M_{200}^{\rm spur}\approx 10^{11}{\rm M_\odot}$. Suppressing spurious heating at smaller/larger radii, or for older/younger stellar populations, requires haloes to be resolved with more/fewer DM particles. △ Less

Submitted 4 September, 2023; v1 submitted 9 June, 2023; originally announced June 2023.

Comments: 17 pages, 15 figures. MNRAS published. A basic python script for calculating spurious disc heating rates is available at https://github.com/AaronDLudlow/hot-disc.git

The impact of stochastic modeling on the predictive power of galaxy formation simulations

Authors: Josh Borrow, Matthieu Schaller, Yannick M. Bahe, Joop Schaye, Aaron D. Ludlow, Sylvia Ploeckinger, Folkert S. J. Nobels, Edoardo Altamura

Abstract: All modern galaxy formation models employ stochastic elements in their sub-grid prescriptions to discretise continuous equations across the time domain. In this paper, we investigate how the stochastic nature of these models, notably star formation, black hole accretion, and their associated feedback, that act on small ($<$ kpc) scales, can back-react on macroscopic galaxy properties (e.g. stellar… ▽ More All modern galaxy formation models employ stochastic elements in their sub-grid prescriptions to discretise continuous equations across the time domain. In this paper, we investigate how the stochastic nature of these models, notably star formation, black hole accretion, and their associated feedback, that act on small ($<$ kpc) scales, can back-react on macroscopic galaxy properties (e.g. stellar mass and size) across long ($>$ Gyr) timescales. We find that the scatter in scaling relations predicted by the EAGLE model implemented in the SWIFT code can be significantly impacted by random variability between re-simulations of the same object, even when galaxies are resolved by tens of thousands of particles. We then illustrate how re-simulations of the same object can be used to better understand the underlying model, by showing how correlations between galaxy stellar mass and black hole mass disappear at the highest black hole masses ($M_{\rm BH} > 10^8$ M$_\odot$), indicating that the feedback cycle may be interrupted by external processes. We find that although properties that are collected cumulatively over many objects are relatively robust against random variability (e.g. the median of a scaling relation), the properties of individual galaxies (such as galaxy stellar mass) can vary by up to 25\%, even far into the well-resolved regime, driven by bursty physics (black hole feedback) and mergers between galaxies. We suggest that studies of individual objects within cosmological simulations be treated with caution, and that any studies aiming to closely investigate such objects must account for random variability within their results. △ Less

Submitted 22 September, 2023; v1 submitted 15 November, 2022; originally announced November 2022.

Comments: Accepted for publication in MNRAS

The physical drivers of gas turbulence in simulated disc galaxies

Authors: Esteban Jiménez, Claudia del P. Lagos, Aaron D. Ludlow, Emily Wisnioski

Abstract: We use the EAGLE cosmological simulations to study the evolution of the vertical velocity dispersion of cold gas, $σ_{z}$, in central disc galaxies and its connection to stellar feedback, gravitational instabilities, cosmological gas accretion and galaxy mergers. To isolate the impact of feedback, we analyse runs that turn off stellar and (or) AGN feedback in addition to a run that includes both.… ▽ More We use the EAGLE cosmological simulations to study the evolution of the vertical velocity dispersion of cold gas, $σ_{z}$, in central disc galaxies and its connection to stellar feedback, gravitational instabilities, cosmological gas accretion and galaxy mergers. To isolate the impact of feedback, we analyse runs that turn off stellar and (or) AGN feedback in addition to a run that includes both. The evolution of $σ_z$ and its dependence on stellar mass and star formation rate in EAGLE are in good agreement with observations. Galaxies hosted by haloes of similar virial mass, $\rm M_{200}$, have similar $σ_z$ values even in runs where feedback is absent. The prevalence of local instabilities in discs is uncorrelated with $σ_z$ at low redshift and becomes only weakly correlated at high redshifts and in galaxies hosted by massive haloes. $σ_z$ correlates most strongly with the specific gas accretion rate onto the disc as well as with the degree of misalignment between the inflowing gas and the disc's rotation axis. These correlations are significant across all redshifts and halo masses, with misaligned accretion being the primary driver of high gas turbulence at redshifts $z \lesssim 1$ and for halo masses $\rm M_{200} \lesssim 10^{11.5} M_{\odot}$. Galaxy mergers increase $σ_z$, but because they are rare in our sample, they play only a minor role in its evolution. Our results suggest that the turbulence of cold gas in EAGLE discs results from a complex interplay of different physical processes whose relative importance depends on halo mass and redshift. △ Less

Submitted 13 July, 2023; v1 submitted 18 October, 2022; originally announced October 2022.

Comments: 22 pages, 12 figures. Accepted for publication in MNRAS

The impact of spurious collisional heating on the morphological evolution of simulated galactic discs

Authors: Matthew J. Wilkinson, Aaron D. Ludlow, Claudia del P. Lagos, S. Michael Fall, Joop Schaye, Danail Obreschkow

Abstract: We use a suite of idealised N-body simulations to study the impact of spurious heating of star particles by dark matter particles on the kinematics and morphology of simulated galactic discs. We find that spurious collisional heating leads to a systematic increase of the azimuthal velocity dispersion ($σ_φ$) of stellar particles and a corresponding decrease in their mean azimuthal velocities (… ▽ More We use a suite of idealised N-body simulations to study the impact of spurious heating of star particles by dark matter particles on the kinematics and morphology of simulated galactic discs. We find that spurious collisional heating leads to a systematic increase of the azimuthal velocity dispersion ($σ_φ$) of stellar particles and a corresponding decrease in their mean azimuthal velocities ($\overline{v}_φ$). The rate of heating is dictated primarily by the number of dark matter halo particles (or equivalently, by the dark matter particle mass at fixed halo mass) and by radial gradients in the local dark matter density along the disc; it is largely insensitive to the stellar particle mass. Galaxies within haloes resolved with fewer than $\approx 10^6$ dark matter particles are particularly susceptible to spurious morphological evolution, irrespective of the total halo mass (with even more particles required to prevent heating of the galactic centre). Collisional heating transforms galactic discs from flattened structures into rounder spheroidal systems, causing them to lose rotational support in the process. It also affects the locations of galaxies in standard scaling relations that link their various properties: at fixed stellar mass, it increases the sizes of galaxies, and reduces their mean stellar rotation velocities and specific angular momenta. Our results urge caution when extrapolating simulated galaxy scaling relations to low masses where spurious collisional effects can bias their normalisation, slope and scatter. △ Less

Submitted 30 January, 2023; v1 submitted 16 August, 2022; originally announced August 2022.

Comments: Minor changes to match MNRAS published version. 20 pages, 13 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 519, Issue 4, March 2023, Pages 5942-5961

The cosmology dependence of the concentration-mass-redshift relation

Authors: Daniel López-Cano, Raúl E. Angulo, Aaron D. Ludlow, M. Zennaro, S. Contreras, Jonás Chaves-Montero, G. Aricò

Abstract: The concentrations of dark matter haloes provide crucial information about their internal structure and how it depends on mass and redshift -- the so-called concentration-mass-redshift relation, denoted $c(M,z)$. We present here an extensive study of the cosmology-dependence of $c(M,z)$ that is based on a suite of 72 gravity-only, full N-body simulations in which the following cosmological paramet… ▽ More The concentrations of dark matter haloes provide crucial information about their internal structure and how it depends on mass and redshift -- the so-called concentration-mass-redshift relation, denoted $c(M,z)$. We present here an extensive study of the cosmology-dependence of $c(M,z)$ that is based on a suite of 72 gravity-only, full N-body simulations in which the following cosmological parameters were varied: $σ_{8}$, $Ω_{\mathrm{M}}$, $Ω_{\mathrm{b}}$, $n_{\mathrm{s}}$, $h$, $M_ν$, $w_{0}$ and $w_{\mathrm{a}}$. We characterize the impact of these parameters on concentrations for different halo masses and redshifts. In agreement with previous works, and for all cosmologies studied, we find that there exists a tight correlation between the characteristic densities of dark matter haloes within their scale radii, $r_{-2}$, and the critical density of the Universe at a suitably defined formation time. This finding, when combined with excursion set modelling of halo formation histories, allows us to accurately predict the concentrations of dark matter haloes as a function of mass, redshift, and cosmology. We use our simulations to test the reliability of a number of published models for predicting halo concentration and highlight when they succeed or fail to reproduce the cosmological $c(M,z)$ relation. △ Less

Submitted 27 July, 2022; originally announced July 2022.

Comments: 11 pages, 9 figures

Sub-recoil clock-transition laser cooling enabling shallow optical lattice clocks

Authors: X. Zhang, K. Beloy, Y. S. Hassan, W. F. McGrew, C-C Chen, J. L. Siegel, T. Grogan, A. D. Ludlow

Abstract: Laser cooling is a key ingredient for quantum control of atomic systems in a variety of settings. In divalent atoms, two-stage Doppler cooling is typically used to bring atoms to the uK regime. Here, we implement a pulsed radial cooling scheme using the ultranarrow 1S0-3P0 clock transition in ytterbium to realize sub-recoil temperatures, down to tens of nK. Together with sideband cooling along the… ▽ More Laser cooling is a key ingredient for quantum control of atomic systems in a variety of settings. In divalent atoms, two-stage Doppler cooling is typically used to bring atoms to the uK regime. Here, we implement a pulsed radial cooling scheme using the ultranarrow 1S0-3P0 clock transition in ytterbium to realize sub-recoil temperatures, down to tens of nK. Together with sideband cooling along the one-dimensional lattice axis, we efficiently prepare atoms in shallow lattices at an energy of 6 lattice recoils. Under these conditions key limits on lattice clock accuracy and instability are reduced, opening the door to dramatic improvements. Furthermore, tunneling shifts in the shallow lattice do not compromise clock accuracy at the 10-19 level. △ Less

Submitted 17 June, 2022; originally announced June 2022.

Optical Atomic Clock aboard an Earth-orbiting Space Station (OACESS): Enhancing searches for physics beyond the standard model in space

Authors: Vladimir Schkolnik, Dmitry Budker, Oliver Fartmann, Victor Flambaum, Leo Hollberg, Tigran Kalaydzhyan, Shimon Kolkowitz, Markus Krutzik, Andrew Ludlow, Nathan Newbury, Christoph Pyrlik, Laura Sinclair, Yevgeny Stadnik, Ingmari Tietje, Jun Ye, Jason Williams

Abstract: We present a concept for a high-precision optical atomic clock (OAC) operating on an Earth-orbiting space station. This pathfinder science mission will compare the space-based OAC with one or more ultra-stable terrestrial OACs to search for space-time-dependent signatures of dark scalar fields that manifest as anomalies in the relative frequencies of station-based and ground-based clocks. This ope… ▽ More We present a concept for a high-precision optical atomic clock (OAC) operating on an Earth-orbiting space station. This pathfinder science mission will compare the space-based OAC with one or more ultra-stable terrestrial OACs to search for space-time-dependent signatures of dark scalar fields that manifest as anomalies in the relative frequencies of station-based and ground-based clocks. This opens the possibility of probing models of new physics that are inaccessible to purely ground-based OAC experiments where a dark scalar field may potentially be strongly screened near Earth's surface. This unique enhancement of sensitivity to potential dark matter candidates harnesses the potential of space-based OACs. △ Less

Submitted 8 August, 2022; v1 submitted 20 April, 2022; originally announced April 2022.

Comments: Submitted to Quantum Sci. Technol. special issue: Cold atom in space

Journal ref: Quantum Sci. Technol. 8, 014003 (2023)

Drivers of asymmetry in synthetic H I emission-line profiles of galaxies in the EAGLE simulation

Authors: Aditya Manuwal, Aaron D. Ludlow, Adam R. H. Stevens, Ruby J. Wright, Aaron S. G. Robotham

Abstract: We study the shapes of spatially integrated H I emission-line profiles of galaxies in the EAGLE simulation using three separate measures of the profile's asymmetry. We show that the subset of EAGLE galaxies whose gas fractions and stellar masses are consistent with those in the xGASS survey also have similar H I line asymmetries. Central galaxies with symmetric H I line profiles typically correspo… ▽ More We study the shapes of spatially integrated H I emission-line profiles of galaxies in the EAGLE simulation using three separate measures of the profile's asymmetry. We show that the subset of EAGLE galaxies whose gas fractions and stellar masses are consistent with those in the xGASS survey also have similar H I line asymmetries. Central galaxies with symmetric H I line profiles typically correspond to rotationally supported H I and stellar disks, but those with asymmetric line profiles may or may not correspond to dispersion-dominated systems. Galaxies with symmetric H I emission lines are, on average, more gas rich than those with asymmetric lines, and also exhibit systematic differences in their specific star formation rates, suggesting that turbulence generated by stellar or AGN feedback may be one factor contributing to H I line asymmetry. The line asymmetry also correlates strongly with the dynamical state of a galaxy's host dark matter halo: older, more relaxed haloes host more-symmetric galaxies than those hosted by unrelaxed ones. At fixed halo mass, asymmetric centrals tend to be surrounded by a larger number of massive subhaloes than their symmetric counterparts, and also experience higher rates of gas accretion and outflow. At fixed stellar mass, central galaxies have, on average, more symmetric H I emission lines than satellites; for the latter, ram pressure and tidal stripping are significant sources of asymmetry. △ Less

Submitted 29 December, 2021; v1 submitted 23 September, 2021; originally announced September 2021.

Comments: Accepted for publication in MNRAS

Optical coherence between atomic species at the second scale: improved clock comparisons via differential spectroscopy

Authors: May E. Kim, William F. McGrew, Nicholas V. Nardelli, Ethan R. Clements, Youssef S. Hassan, Xiaogang Zhang, Jose L. Valencia, Holly Leopardi, David B. Hume, Tara M. Fortier, Andrw D. Ludlow, David R. Leibrandt

Abstract: Comparisons of high-accuracy optical atomic clocks \cite{Ludlow2015} are essential for precision tests of fundamental physics \cite{Safronova2018}, relativistic geodesy \cite{McGrew2018, Grotti2018, Delva2019}, and the anticipated redefinition of the SI second \cite{Riehle2018}. The scientific reach of these applications is restricted by the statistical precision of interspecies comparison measure… ▽ More Comparisons of high-accuracy optical atomic clocks \cite{Ludlow2015} are essential for precision tests of fundamental physics \cite{Safronova2018}, relativistic geodesy \cite{McGrew2018, Grotti2018, Delva2019}, and the anticipated redefinition of the SI second \cite{Riehle2018}. The scientific reach of these applications is restricted by the statistical precision of interspecies comparison measurements. The instability of individual clocks is limited by the finite coherence time of the optical local oscillator (OLO), which bounds the maximum atomic interrogation time. In this letter, we experimentally demonstrate differential spectroscopy \cite{Hume2016}, a comparison protocol that enables interrogating beyond the OLO coherence time. By phase-coherently linking a zero-dead-time (ZDT) \cite{Schioppo2017} Yb optical lattice clock with an Al$^+$ single-ion clock via an optical frequency comb and performing synchronised Ramsey spectroscopy, we show an improvement in comparison instability relative to our previous result \cite{network2020frequency} of nearly an order of magnitude. To our knowledge, this result represents the most stable interspecies clock comparison to date. △ Less

Submitted 24 November, 2021; v1 submitted 20 September, 2021; originally announced September 2021.

Spurious heating of stellar motions in simulated galactic disks by dark matter halo particles

Authors: Aaron D. Ludlow, S. Michael Fall, Joop Schaye, Danail Obreschkow

Abstract: We use idealized N-body simulations of equilibrium stellar disks embedded within course-grained dark matter haloes to study the effects of spurious collisional heating on disk structure and kinematics. Collisional heating artificially increases the vertical and radial velocity dispersions of disk stars, as well as the thickness and size of disks; the effects are felt at all galacto-centric radii.… ▽ More We use idealized N-body simulations of equilibrium stellar disks embedded within course-grained dark matter haloes to study the effects of spurious collisional heating on disk structure and kinematics. Collisional heating artificially increases the vertical and radial velocity dispersions of disk stars, as well as the thickness and size of disks; the effects are felt at all galacto-centric radii. The integrated effects of collisional heating are determined by the mass of dark matter halo particles (or equivalently, by the number of particles at fixed halo mass), their local density and characteristic velocity dispersion, but are largely insensitive to the stellar particle mass. The effects can therefore be reduced by increasing the mass resolution of dark matter in cosmological simulations, with limited benefits from increasing the baryonic (or stellar) mass resolution. We provide a simple empirical model that accurately captures the effects of spurious collisional heating on the structure and kinematics of simulated disks, and use it to assess the importance of disk heating for simulations of galaxy formation. We find that the majority of state-of-the-art zoom simulations, and a few of the highest-resolution, smallest-volume cosmological runs, are in principle able to resolve thin stellar disks in Milky Way-mass haloes, but most large-volume cosmological simulations cannot. For example, dark matter haloes resolved with fewer than $\approx 10^6$ particles will collisionally heat stars near the stellar half-mass radius such that their vertical velocity dispersion increases by $\gtrsim 10$ per cent of the halo's virial velocity in approximately one Hubble time. △ Less

Submitted 23 September, 2021; v1 submitted 7 May, 2021; originally announced May 2021.

Comments: Minor changes to match MNRAS accepted version. Main text: 19 pages, 10 Figure (Appendix: 5 pages, 6 figures)

Characterization and suppression of background light shifts in an optical lattice clock

Authors: R. J. Fasano, Y. J. Chen, W. F. McGrew, W. J. Brand, R. W. Fox, A. D. Ludlow

Abstract: Experiments involving optical traps often require careful control of the ac Stark shifts induced by strong confining light fields. By carefully balancing light shifts between two atomic states of interest, optical traps at the magic wavelength have been especially effective at suppressing deleterious effects stemming from such shifts. Highlighting the power of this technique, optical clocks today… ▽ More Experiments involving optical traps often require careful control of the ac Stark shifts induced by strong confining light fields. By carefully balancing light shifts between two atomic states of interest, optical traps at the magic wavelength have been especially effective at suppressing deleterious effects stemming from such shifts. Highlighting the power of this technique, optical clocks today exploit Lamb-Dicke confinement in magic-wavelength optical traps, in some cases realizing shift cancellation at the ten parts per billion level. Theory and empirical measurements can be used at varying levels of precision to determine the magic wavelength where shift cancellation occurs. However, lasers exhibit background spectra from amplified spontaneous emission or other lasing modes which can easily contaminate measurement of the magic wavelength and its reproducibility in other experiments or conditions. Indeed, residual light shifts from laser background have plagued optical lattice clock measurements for years. In this work, we develop a simple theoretical model allowing prediction of light shifts from measured background spectra. We demonstrate good agreement between this model and measurements of the background light shift from an amplified diode laser in an Yb optical lattice clock. Additionally, we model and experimentally characterize the filtering effect of a volume Bragg grating bandpass filter, demonstrating that application of the filter can reduce background light shifts from amplified spontaneous emission well below the $10^{-18}$ fractional clock frequency level. This demonstration is corroborated by direct clock comparisons between a filtered amplified diode laser and a filtered titanium:sapphire laser. △ Less

Submitted 22 March, 2021; originally announced March 2021.

Comments: 11 pages, 7 figures

The orbital evolution of UFDs and GCs in an evolving Galactic potential

Authors: Benjamin M. Armstrong, Kenji Bekki, Aaron D. Ludlow

Abstract: We use the second Gaia data release to investigate the kinematics of 17 ultra-faint dwarf galaxies (UFDs) and 154 globular clusters (GCs) in the Milky Way, focusing on the differences between static and evolving models of the Galactic potential. An evolving potential modifies a satellite's orbit relative to its static equivalent, though the difference is small compared to existing uncertainties on… ▽ More We use the second Gaia data release to investigate the kinematics of 17 ultra-faint dwarf galaxies (UFDs) and 154 globular clusters (GCs) in the Milky Way, focusing on the differences between static and evolving models of the Galactic potential. An evolving potential modifies a satellite's orbit relative to its static equivalent, though the difference is small compared to existing uncertainties on orbital parameters. We find that the UFD Boötes II is likely on its first passage around the Milky Way. Depending on the assumed mass of the Milky Way, the UFDs Triangulum II, Hydrus I, Coma Berenices, Draco II, and Ursa Major II, as well as the GC Pyxis, may also be on first infall so may be useful for constraining the mass of the Galaxy. We identify a clear kinematic distinction between metal-rich (${\rm [Fe/H]}>-1.1$) and metal-poor GCs (${\rm [Fe/H]}\leq-1.1$). Although most metal-rich clusters occupy predominately prograde orbits, with low eccentricities ($e\approx 0.35$) and similar specific angular momenta and orbital planes as the Galactic disc, 7 show potentially retrograde orbits, the origin of which is unclear. Metal-poor clusters have more diverse orbits, higher eccentricities ($e\approx 0.65$), and half have orbital planes offset from the disc by 60 to 120 degrees. The UFDs have similar $θ$ and $φ$ to the metal-poor GCs, suggesting a similar origin. We provide a catalogue of orbital parameters for UFDs and GCs for two different Galaxy masses and their observational uncertainties. △ Less

Submitted 6 December, 2020; v1 submitted 25 November, 2020; originally announced November 2020.

Comments: 23 pages, 12 figures, accepted for publication in MNRAS

Pre-processing, group accretion and the orbital trajectories of associated subhaloes

Authors: Lucie Bakels, Aaron D. Ludlow, Chris Power

Abstract: We use a high-resolution cosmological dark matter-only simulation to study the orbital trajectories of haloes and subhaloes in the environs of isolated hosts. We carefully tally all apsis points and use them to distinguish haloes that are infalling for the first time from those that occupy more evolved orbits. We find that roughly 21 per cent of subhaloes within a host's virial radius are currentl… ▽ More We use a high-resolution cosmological dark matter-only simulation to study the orbital trajectories of haloes and subhaloes in the environs of isolated hosts. We carefully tally all apsis points and use them to distinguish haloes that are infalling for the first time from those that occupy more evolved orbits. We find that roughly 21 per cent of subhaloes within a host's virial radius are currently on first infall, and have not yet reached their first orbital pericentre; roughly 44 per cent are still approaching their first apocentre after infall. For the range of host masses studied, roughly half of all accreted systems were pre-processed prior to infall, and about 20 per cent were accreted in groups. We confirm that the entire population of accreted subhaloes -- often referred to as "associated" subhaloes -- extend far beyond the virial radii of their hosts, with roughly half currently residing at distances that exceed $\approx 1.2\times r_{200}$. Many of these backsplash haloes have gained orbital energy since infall, and occupy extreme orbits that carry them well past their initial turnaround radii. Such extreme orbits are created during the initial accretion and dissolution of loosely bound groups, but also through penetrating encounters between subhaloes on subsequent orbits. The same processes may also give rise to unexpectedly abrupt losses of orbital energy. These effects combine, giving rise to a large variation in the ratio of sequent apocentres for accreted systems. We find that, within 2 virial radii from host centres, the concentrations of first-infall halos are remarkably similar those of isolated field halos, whereas backsplash haloes, as well as systems that were pre-processed, are considerably more concentrated. △ Less

Submitted 7 January, 2021; v1 submitted 12 August, 2020; originally announced August 2020.

Comments: 17 pages, 12 figures, submitted to MNRAS

Observational constraints on the slope of the radial acceleration relation at low accelerations

Authors: Kyle A. Oman, Margot M. Brouwer, Aaron D. Ludlow, Julio F. Navarro

Abstract: The radial acceleration relation (RAR) locally relates the `observed' acceleration inferred from the dynamics of a system to the acceleration implied by its baryonic matter distribution. The relation as traced by galaxy rotation curves is one-to-one with remarkably little scatter, implying that the dynamics of a system can be predicted simply by measuring its density profile as traced by e.g. stel… ▽ More The radial acceleration relation (RAR) locally relates the `observed' acceleration inferred from the dynamics of a system to the acceleration implied by its baryonic matter distribution. The relation as traced by galaxy rotation curves is one-to-one with remarkably little scatter, implying that the dynamics of a system can be predicted simply by measuring its density profile as traced by e.g. stellar light or gas emission lines. Extending the relation to accelerations below those usually probed by practically observable kinematic tracers is challenging, especially once accounting for faintly emitting baryons, such as the putative warm-hot intergalactic medium, becomes important. We show that in the low-acceleration regime, the (inverted) RAR predicts an unphysical, declining enclosed baryonic mass profile for systems with `observed' acceleration profiles steeper than $g_{\rm obs}\propto r^{-1}$ (corresponding to density profiles steeper than isothermal - $ρ(r)\propto r^{-2}$). If the RAR is tantamount to a natural law, such acceleration profiles cannot exist. We apply this argument to test the compatibility of an extrapolation of the rotation curve-derived RAR to low accelerations with data from galaxy-galaxy weak lensing, dwarf spheroidal galaxy stellar kinematic, and outer Milky~Way dynamical measurements, fully independent of the uncertainties inherent in direct measurements of the baryonic matter distribution. In all cases we find that the data weakly favour a break to a steeper low-acceleration slope. Improvements in measurements and modelling of the outer Milky~Way, and weak lensing, seem like the most promising path toward stronger constraints on the low-acceleration behaviour of the RAR. △ Less

Submitted 16 December, 2020; v1 submitted 11 June, 2020; originally announced June 2020.

Comments: 11 pages, 8 figures (+appendices), v2 corrects a factual error in Sec. 2

Optical Atomic Clock Comparison through Turbulent Air

Authors: Martha I. Bodine, Jean-Daniel Deschênes, Isaac H. Khader, William C. Swann, Holly Leopardi, Kyle Beloy, Tobias Bothwell, Samuel M. Brewer, Sarah L. Bromley, Jwo-Sy Chen, Scott A. Diddams, Robert J. Fasano, Tara M. Fortier, Youssef S. Hassan, David B. Hume, Dhruv Kedar, Colin J. Kennedy, Amanda Koepke, David R. Leibrandt, Andrew D. Ludlow, William F. McGrew, William R. Milner, Daniele Nicolodi, Eric Oelker, Thomas E. Parker , et al. (10 additional authors not shown)

Abstract: We use frequency comb-based optical two-way time-frequency transfer (O-TWTFT) to measure the optical frequency ratio of state-of-the-art ytterbium and strontium optical atomic clocks separated by a 1.5 km open-air link. Our free-space measurement is compared to a simultaneous measurement acquired via a noise-cancelled fiber link. Despite non-stationary, ps-level time-of-flight variations in the fr… ▽ More We use frequency comb-based optical two-way time-frequency transfer (O-TWTFT) to measure the optical frequency ratio of state-of-the-art ytterbium and strontium optical atomic clocks separated by a 1.5 km open-air link. Our free-space measurement is compared to a simultaneous measurement acquired via a noise-cancelled fiber link. Despite non-stationary, ps-level time-of-flight variations in the free-space link, ratio measurements obtained from the two links, averaged over 30.5 hours across six days, agree to $6\times10^{-19}$, showing that O-TWTFT can support free-space atomic clock comparisons below the $10^{-18}$ level. △ Less

Submitted 11 September, 2020; v1 submitted 1 June, 2020; originally announced June 2020.

Journal ref: Phys. Rev. Research 2, 033395 (2020)

Frequency Ratio Measurements with 18-digit Accuracy Using a Network of Optical Clocks

Authors: Boulder Atomic Clock Optical Network, Collaboration, :, Kyle Beloy, Martha I. Bodine, Tobias Bothwell, Samuel M. Brewer, Sarah L. Bromley, Jwo-Sy Chen, Jean-Daniel Deschênes, Scott A. Diddams, Robert J. Fasano, Tara M. Fortier, Youssef S. Hassan, David B. Hume, Dhruv Kedar, Colin J. Kennedy, Isaac Khader, Amanda Koepke, David R. Leibrandt, Holly Leopardi, Andrew D. Ludlow, William F. McGrew, William R. Milner, Nathan R. Newbury , et al. (13 additional authors not shown)

Abstract: Atomic clocks occupy a unique position in measurement science, exhibiting higher accuracy than any other measurement standard and underpinning six out of seven base units in the SI system. By exploiting higher resonance frequencies, optical atomic clocks now achieve greater stability and lower frequency uncertainty than existing primary standards. Here, we report frequency ratios of the $^{27}$Al… ▽ More Atomic clocks occupy a unique position in measurement science, exhibiting higher accuracy than any other measurement standard and underpinning six out of seven base units in the SI system. By exploiting higher resonance frequencies, optical atomic clocks now achieve greater stability and lower frequency uncertainty than existing primary standards. Here, we report frequency ratios of the $^{27}$Al$^+$, $^{171}$Yb and $^{87}$Sr optical clocks in Boulder, Colorado, measured across an optical network spanned by both fiber and free-space links. These ratios have been evaluated with measurement uncertainties between $6\times10^{-18}$ and $8\times10^{-18}$, making them the most accurate reported measurements of frequency ratios to date. This represents a critical step towards redefinition of the SI second and future applications such as relativistic geodesy and tests of fundamental physics. △ Less

Submitted 29 May, 2020; originally announced May 2020.

Comments: 51 pages, 12 figures, 6 tables

Modeling motional energy spectra and lattice light shifts in optical lattice clocks

Authors: K. Beloy, W. F. McGrew, X. Zhang, D. Nicolodi, R. J. Fasano, Y. S. Hassan, R. C. Brown, A. D. Ludlow

Abstract: We develop a model to describe the motional (i.e., external degree of freedom) energy spectra of atoms trapped in a one-dimensional optical lattice, taking into account both axial and radial confinement relative to the lattice axis. Our model respects the coupling between axial and radial degrees of freedom, as well as other anharmonicities inherent in the confining potential. We further demonstra… ▽ More We develop a model to describe the motional (i.e., external degree of freedom) energy spectra of atoms trapped in a one-dimensional optical lattice, taking into account both axial and radial confinement relative to the lattice axis. Our model respects the coupling between axial and radial degrees of freedom, as well as other anharmonicities inherent in the confining potential. We further demonstrate how our model can be used to characterize lattice light shifts in optical lattice clocks, including shifts due to higher multipolar (magnetic dipole and electric quadrupole) and higher order (hyperpolarizability) coupling to the lattice field. We compare results for our model with results from other lattice light shift models in the literature under similar conditions. △ Less

Submitted 13 April, 2020; originally announced April 2020.

Comments: 19 pages; accepted to PRA

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

A comparison of $\text{H}_2$ formation models at high redshift

Authors: Alexander Schäbe, Emilio Romano-Díaz, Cristiano Porciani, Aaron D. Ludlow, Matteo Tomassetti

Abstract: Modelling the molecular gas that is routinely detected through CO observations of high-redshift galaxies constitutes a major challenge for ab initio simulations of galaxy formation. We carry out a suite of cosmological hydrodynamic simulations to compare three approximate methods that have been used in the literature to track the formation and evolution of the simplest and most abundant molecule,… ▽ More Modelling the molecular gas that is routinely detected through CO observations of high-redshift galaxies constitutes a major challenge for ab initio simulations of galaxy formation. We carry out a suite of cosmological hydrodynamic simulations to compare three approximate methods that have been used in the literature to track the formation and evolution of the simplest and most abundant molecule, H$_2$. Namely, we consider: i) a semi-empirical procedure that associates H$_2$ to dark-matter haloes based on a series of scaling relations inferred from observations, ii) a model that assumes chemical equilibrium between the H$_2$ formation and destruction rates, and iii) a model that fully solves the out-of-equilibrium rate equations and accounts for the unresolved structure of molecular clouds. We study the impact of finite spatial resolution and show that robust H$_2$ masses at redshift $z\approx 4$ can only be obtained for galaxies that are sufficiently metal enriched in which H$_2$ formation is fast. This corresponds to H$_2$ reservoirs with masses $M_{\mathrm{H_2}}\gtrsim 6\times 10^9 \mathrm{M}_\odot$. In this range, equilibrium and non-equilibrium models predict similar molecular masses (but different galaxy morphologies) while the semi-empirical method produces less H$_2$. The star formation rates as well as the stellar and H$_2$ masses of the simulated galaxies are in line with those observed in actual galaxies at similar redshifts that are not massive starbursts. The H$_2$ mass functions extracted from the simulations at $z\approx 4$ agree well with recent observations that only sample the high-mass end. However, our results indicate that most molecular material at high $z$ lies yet undetected in reservoirs with $10^9<M_{\mathrm H_2}<10^{10} \mathrm{M}_\odot$. △ Less

Submitted 3 September, 2020; v1 submitted 9 March, 2020; originally announced March 2020.

Comments: 19 pages, 12 figures, 1 table, accepted for publication in MNRAS, comments are welcome

Coherent Optical Clock Down-Conversion for Microwave Frequencies with 10-18 Instability

Authors: Takuma Nakamura, Josue Davila-Rodriguez, Holly Leopardi, Jeff A. Sherman, Tara M. Fortier, Xiaojun Xie, Joe C. Campbell, William F. McGrew, Xiaogang Zhang, Youssef S. Hassan, Daniele Nicolodi, Kyle Beloy, Andrew D. Ludlow, Scott A. Diddams, Franklyn Quinlan

Abstract: Optical atomic clocks are poised to redefine the SI second, thanks to stability and accuracy more than one hundred times better than the current microwave atomic clock standard. However, the best optical clocks have not seen their performance transferred to the electronic domain, where radar, navigation, communications, and fundamental research rely on less stable microwave sources. By comparing t… ▽ More Optical atomic clocks are poised to redefine the SI second, thanks to stability and accuracy more than one hundred times better than the current microwave atomic clock standard. However, the best optical clocks have not seen their performance transferred to the electronic domain, where radar, navigation, communications, and fundamental research rely on less stable microwave sources. By comparing two independent optical-to-electronic signal generators, we demonstrate a 10 GHz microwave signal with phase that exactly tracks that of the optical clock phase from which it is derived, yielding an absolute fractional frequency instability of 1*10-18 in the electronic domain. Such faithful reproduction of the optical clock phase expands the opportunities for optical clocks both technologically and scientifically for time-dissemination, navigation, and long-baseline interferometric imaging. △ Less

Submitted 9 March, 2020; v1 submitted 5 March, 2020; originally announced March 2020.

Comments: 19 page, 10 figures (including Supplementary Text)

Local Group star formation in warm and self-interacting dark matter cosmologies

Authors: Mark R. Lovell, Wojciech Hellwing, Aaron Ludlow, Jesús Zavala, Andrew Robertson, Azadeh Fattahi, Carlos S. Frenk, Jennifer Hardwick

Abstract: The nature of the dark matter can affect the collapse time of dark matter haloes, and can therefore be imprinted in observables such as the stellar population ages and star formation histories of dwarf galaxies. In this paper we use high resolution hydrodynamical simulations of Local Group-analogue (LG) volumes in cold dark matter (CDM), sterile neutrino warm dark matter (WDM) and self-interacting… ▽ More The nature of the dark matter can affect the collapse time of dark matter haloes, and can therefore be imprinted in observables such as the stellar population ages and star formation histories of dwarf galaxies. In this paper we use high resolution hydrodynamical simulations of Local Group-analogue (LG) volumes in cold dark matter (CDM), sterile neutrino warm dark matter (WDM) and self-interacting dark matter (SIDM) models with the EAGLE galaxy formation code to study how galaxy formation times change with dark matter model. We are able to identify the same haloes in different simulations, since they share the same initial density field phases. We find that the stellar mass of galaxies depends systematically on resolution, and can differ by as much as a factor of two in haloes of a given dark matter mass. The evolution of the stellar populations in SIDM is largely identical to that of CDM, but in WDM early star formation is instead suppressed. The time at which LG haloes can begin to form stars through atomic cooling is delayed by $\sim$200~Myr in WDM models compared to CDM. It will be necessary to measure stellar ages of old populations to a precision of better than 100~Myr, and to address degeneracies with the redshift of reionization -- and potentially other baryonic processes -- in order to use these observables to distinguish between dark matter models. △ Less

Submitted 25 September, 2020; v1 submitted 25 February, 2020; originally announced February 2020.

Comments: 17 pages, 13 figures. Version accepted by MNRAS. Contact: lovell@hi.is

Characterising the Structure of Halo Merger Trees Using a Single Parameter: The Tree Entropy

Authors: Danail Obreschkow, Pascal J. Elahi, Claudia del P. Lagos, Rhys J. J. Poulton, Aaron D. Ludlow

Abstract: Linking the properties of galaxies to the assembly history of their dark matter haloes is a central aim of galaxy evolution theory. This paper introduces a dimensionless parameter $s\in[0,1]$, the "tree entropy", to parametrise the geometry of a halo's entire mass assembly hierarchy, building on a generalisation of Shannon's information entropy. By construction, the minimum entropy ($s=0$) corresp… ▽ More Linking the properties of galaxies to the assembly history of their dark matter haloes is a central aim of galaxy evolution theory. This paper introduces a dimensionless parameter $s\in[0,1]$, the "tree entropy", to parametrise the geometry of a halo's entire mass assembly hierarchy, building on a generalisation of Shannon's information entropy. By construction, the minimum entropy ($s=0$) corresponds to smoothly assembled haloes without any mergers. In contrast, the highest entropy ($s=1$) represents haloes grown purely by equal-mass binary mergers. Using simulated merger trees extracted from the cosmological $N$-body simulation SURFS, we compute the natural distribution of $s$, a skewed bell curve peaking near $s=0.4$. This distribution exhibits weak dependences on halo mass $M$ and redshift $z$, which can be reduced to a single dependence on the relative peak height $δ_{\rm c}/σ(M,z)$ in the matter perturbation field. By exploring the correlations between $s$ and global galaxy properties generated by the SHARK semi-analytic model, we find that $s$ contains a significant amount of information on the morphology of galaxies $-$ in fact more information than the spin, concentration and assembly time of the halo. Therefore, the tree entropy provides an information-rich link between galaxies and their dark matter haloes. △ Less

Submitted 11 February, 2020; v1 submitted 27 November, 2019; originally announced November 2019.

Comments: 20 pages, 18 figures

To beta or not to beta: can higher-order Jeans analysis break the mass-anisotropy degeneracy in simulated dwarfs?

Authors: Anna Genina, Justin I. Read, Carlos S. Frenk, Shaun Cole, Alejandro Benitez-Llambay, Aaron D. Ludlow, Julio F. Navarro, Kyle A. Oman, Andrew Robertson

Abstract: We test a non-parametric higher-order Jeans analysis method, GravSphere, on 32 simulated dwarf galaxies comparable to classical Local Group dwarfs like Fornax. The galaxies are selected from the APOSTLE suite of cosmological hydrodynamics simulations with Cold Dark Matter (CDM) and Self-Interacting Dark Matter (SIDM) models, allowing us to investigate cusps and cores in density distributions. We f… ▽ More We test a non-parametric higher-order Jeans analysis method, GravSphere, on 32 simulated dwarf galaxies comparable to classical Local Group dwarfs like Fornax. The galaxies are selected from the APOSTLE suite of cosmological hydrodynamics simulations with Cold Dark Matter (CDM) and Self-Interacting Dark Matter (SIDM) models, allowing us to investigate cusps and cores in density distributions. We find that, for CDM dwarfs, the recovered enclosed mass profiles have a bias of no more than 10 per cent, with a 50 per cent scatter in the inner regions and a 20 per cent scatter near the half-light radius, consistent with standard mass estimators. The density profiles are also recovered with a bias of no more than 10 per cent and a scatter of 30 per cent in the inner regions. For SIDM dwarfs, the mass and density profiles are recovered within our 95 per cent confidence intervals, but are biased towards cuspy dark matter distributions. This is mainly due to a lack of sufficient constraints from the data. We explore the sources of scatter in the accuracy of the recovered profiles and suggest a $χ^2$ statistic to separate successful models from biased ones. Finally, we show that the uncertainties on the mass profiles obtained with GravSphere are smaller than those for comparable Jeans methods, and that they can be further improved if stronger priors, motivated by cosmological simulations, are placed on the velocity anisotropy. We conclude that GravSphere is a promising Jeans-based approach for modelling dark matter distributions in dwarf galaxies. △ Less

Submitted 5 August, 2020; v1 submitted 20 November, 2019; originally announced November 2019.

Comments: 22 pages, 13 figures. Accepted in MNRAS

Baryonic clues to the puzzling diversity of dwarf galaxy rotation curves

Authors: Isabel M. E. Santos-Santos, Julio F. Navarro, Andrew Robertson, Alejandro Benítez-Llambay, Kyle A. Oman, Mark R. Lovell, Carlos S. Frenk, Aaron D. Ludlow, Azadeh Fattahi, Adam Ritz

Abstract: We use a compilation of disc galaxy rotation curves to assess the role of the luminous component ("baryons") in the rotation curve diversity problem. As in earlier work, we find that rotation curve shape correlates with baryonic surface density: high surface density galaxies have rapidly-rising rotation curves consistent with cuspy cold dark matter halos; slowly-rising rotation curves (characteris… ▽ More We use a compilation of disc galaxy rotation curves to assess the role of the luminous component ("baryons") in the rotation curve diversity problem. As in earlier work, we find that rotation curve shape correlates with baryonic surface density: high surface density galaxies have rapidly-rising rotation curves consistent with cuspy cold dark matter halos; slowly-rising rotation curves (characteristic of galaxies with inner mass deficits or "cores") occur only in low surface density galaxies. The correlation, however, seems too weak to be the main driver of the diversity. In addition, dwarf galaxies exhibit a clear trend, from "cuspy" systems where baryons are unimportant in the inner mass budget to "cored" galaxies where baryons actually dominate. This trend constrains the various scenarios proposed to explain the diversity, such as (i) baryonic inflows and outflows during galaxy formation; (ii) dark matter self-interactions; (iii) variations in the baryonic mass structure coupled to rotation velocities through the "mass discrepancy-acceleration relation" (MDAR); or (iv) non-circular motions in gaseous discs. Together with analytical modeling and cosmological hydrodynamical simulations, our analysis shows that each of these scenarios has promising features, but none seems to fully account for the observed diversity. The MDAR, in particular, is inconsistent with the observed trend between rotation curve shape and baryonic importance; either the trend is caused by systematic errors in the data or the MDAR does not apply. The origin of the dwarf galaxy rotation curve diversity and its relation to the structure of cold dark matter halos remains an open issue. △ Less

Submitted 28 April, 2020; v1 submitted 20 November, 2019; originally announced November 2019.

Comments: 20 pages, 16 figures, 1 table. Accepted for publication in MNRAS. Updated acknowledgements

Numerical convergence of hydrodynamical simulations of galaxy formation: the abundance and internal structure of galaxies and their cold dark matter haloes

Authors: Aaron D. Ludlow, Joop Schaye, Matthieu Schaller, Richard Bower

Abstract: We address the issue of numerical convergence in cosmological smoothed particle hydrodynamics simulations using a suite of runs drawn from the EAGLE project. Our simulations adopt subgrid models that produce realistic galaxy populations at a fiducial mass and force resolution, but systematically vary the latter in order to study their impact on galaxy properties. We provide several analytic criter… ▽ More We address the issue of numerical convergence in cosmological smoothed particle hydrodynamics simulations using a suite of runs drawn from the EAGLE project. Our simulations adopt subgrid models that produce realistic galaxy populations at a fiducial mass and force resolution, but systematically vary the latter in order to study their impact on galaxy properties. We provide several analytic criteria that help guide the selection of gravitational softening for hydrodynamical simulations, and present results from runs that both adhere to and deviate from them. Unlike dark matter-only simulations, hydrodynamical simulations exhibit a strong sensitivity to gravitational softening, and care must be taken when selecting numerical parameters. Our results--which focus mainly on star formation histories, galaxy stellar mass functions and sizes--illuminate three main considerations. First, softening imposes a minimum resolved escape speed, $v_ε$, due to the binding energy between gas particles. Runs that adopt such small softening lengths that $v_ε\gt 10\,{\rm km s^{-1}}$ (the sound speed in ionised $\sim 10^4\,{\rm K}$ gas) suffer from reduced effects of photo-heating. Second, feedback from stars or active galactic nuclei may suffer from numerical over-cooling if the gravitational softening length is chosen below a critical value, $ε_{\rm eFB}$. Third, we note that small softening lengths exacerbate the segregation of stars and dark matter particles in halo centres, often leading to the counter-intuitive result that galaxy sizes {\em increase} as softening is reduced. The structure of dark matter haloes in hydrodynamical runs respond to softening in a way that reflects the sensitivity of their galaxy populations to numerical parameters. △ Less

Submitted 7 May, 2020; v1 submitted 14 August, 2019; originally announced August 2019.

Comments: 28 pages, 15 figures; Accepted 2020 January 29 (Figures with reduced resolution; hi-res figures available on request)

Ramsey-Bordé matter-wave interferometry for laser frequency stabilization at $10^{-16}$ frequency instability and below

Authors: Judith Olson, Richard W. Fox, Tara M. Fortier, Todd F. Sheerin, Roger C. Brown, Holly Leopardi, Richard E. Stoner, Chris W. Oates, Andrew D. Ludlow

Abstract: We demonstrate Ramsey-Bordé (RB) atom interferometry for high performance laser stabilization with fractional frequency instability $<2 \times 10^{-16}$ for timescales between 10 and 1000s. The RB spectroscopy laser interrogates two counterpropagating $^{40}$Ca beams on the $^1$S$_0$ -- $^3$P$_1$ transition at 657 nm, yielding 1.6 kHz linewidth interference fringes. Fluorescence detection of the e… ▽ More We demonstrate Ramsey-Bordé (RB) atom interferometry for high performance laser stabilization with fractional frequency instability $<2 \times 10^{-16}$ for timescales between 10 and 1000s. The RB spectroscopy laser interrogates two counterpropagating $^{40}$Ca beams on the $^1$S$_0$ -- $^3$P$_1$ transition at 657 nm, yielding 1.6 kHz linewidth interference fringes. Fluorescence detection of the excited state population is performed on the (4s4p) $^3$P$_1$ -- (4p$^2$) $^3$P$_0$ transition at 431 nm. Minimal thermal shielding and no vibration isolation are used. These stability results surpass performance from other thermal atomic or molecular systems by one to two orders of magnitude, and further improvements look feasible. △ Less

Submitted 15 July, 2019; originally announced July 2019.

Energy equipartition between stellar and dark matter particles in cosmological simulations results in spurious growth of galaxy sizes

Authors: Aaron D. Ludlow, Joop Schaye, Matthieu Schaller, Jack Richings

Abstract: The impact of 2-body scattering on the innermost density profiles of dark matter haloes is well established. We use a suite of cosmological simulations and idealised numerical experiments to show that 2-body scattering is exacerbated in situations where there are two species of unequal mass. This is a consequence of mass segregation and reflects a flow of kinetic energy from the more to less massi… ▽ More The impact of 2-body scattering on the innermost density profiles of dark matter haloes is well established. We use a suite of cosmological simulations and idealised numerical experiments to show that 2-body scattering is exacerbated in situations where there are two species of unequal mass. This is a consequence of mass segregation and reflects a flow of kinetic energy from the more to less massive particles. This has important implications for the interpretation of galaxy sizes in cosmological hydrodynamic simulations, which nearly always model stars with less massive particles than are used for the dark matter. We compare idealised models as well as simulations from the EAGLE project that differ only in the mass resolution of the dark matter component, but keep sub-grid physics, baryonic mass resolution and gravitational force softening fixed. If the dark matter particle mass exceeds the mass of stellar particles, then galaxy sizes--quantified by their projected half-mass radii, ${\rm R_{50}}$--increase systematically with time until ${\rm R_{50}}$ exceeds a small fraction of the redshift-dependent mean inter-particle separation, $l$ (${\rm R_{50}}\geq 0.05\times l$). Our conclusions should also apply to simulations that adopt different hydrodynamic solvers, subgrid physics or adaptive softening, but in that case may need quantitative revision. Any simulation employing a stellar-to-dark matter particle mass ratio greater than unity will escalate spurious energy transfer from dark matter to baryons on small scales. △ Less

Submitted 10 July, 2019; v1 submitted 24 March, 2019; originally announced March 2019.

Comments: 6 pages, 3 figures; replaced by MNRAS accepted version

The connection between halo concentrations and assembly histories: a probe of gravity?

Authors: Piotr Oleśkiewicz, Carlton M. Baugh, Aaron Ludlow

Abstract: We use two high resolution N-body simulations, one assuming general relativity and the other the Hu-Sawicki form of $f\left(R\right)$ gravity with $\vert \bar{f}_{R} \vert = 10^{-6}$, to investigate the concentration--formation time relation of dark matter haloes. We stack haloes in logarithmically spaced mass bins to fit median density profiles and extract median formation times. At fixed mass, h… ▽ More We use two high resolution N-body simulations, one assuming general relativity and the other the Hu-Sawicki form of $f\left(R\right)$ gravity with $\vert \bar{f}_{R} \vert = 10^{-6}$, to investigate the concentration--formation time relation of dark matter haloes. We stack haloes in logarithmically spaced mass bins to fit median density profiles and extract median formation times. At fixed mass, haloes in modified gravity are more concentrated than those in GR, especially at low masses and at low redshift, and do not follow the concentration--formation time relation seen in GR. We assess the sensitivity of the relation to how concentration and formation time are defined, as well as to the segregation of the halo population by the amount of gravitational screening. We find a clear difference between halo concentrations and assembly histories displayed in modified gravity and those in GR. Existing models for the mass--concentration--redshift relation that have gained success in cold and warm dark matter models require revision in $f\left(R\right)$ gravity. △ Less

Submitted 18 March, 2019; originally announced March 2019.

Comments: 12 pages, 10 figures

Measurements of $^{27}$Al$^{+}$ and $^{25}$Mg$^{+}$ magnetic constants for improved ion clock accuracy

Authors: S. M. Brewer, J. -S. Chen, K. Beloy, A. M. Hankin, E. R. Clements, C. W. Chou, W. F. McGrew, X. Zhang, R. J. Fasano, D. Nicolodi, H. Leopardi, T. M. Fortier, S. A. Diddams, A. D. Ludlow, D. J. Wineland, D. R. Leibrandt, D. B. Hume

Abstract: We have measured the quadratic Zeeman coefficient for the ${^{1}S_{0} \leftrightarrow {^{3}P_{0}}}$ optical clock transition in $^{27}$Al$^{+}$, $C_{2}=-71.944(24)$~MHz/T$^{2}$, and the unperturbed hyperfine splitting of the $^{25}$Mg$^{+}$ $^{2}S_{1/2}$ ground electronic state, $ΔW / h = 1~788~762~752.85(13)$~Hz, with improved uncertainties. Both constants are relevant to the evaluation of the… ▽ More We have measured the quadratic Zeeman coefficient for the ${^{1}S_{0} \leftrightarrow {^{3}P_{0}}}$ optical clock transition in $^{27}$Al$^{+}$, $C_{2}=-71.944(24)$~MHz/T$^{2}$, and the unperturbed hyperfine splitting of the $^{25}$Mg$^{+}$ $^{2}S_{1/2}$ ground electronic state, $ΔW / h = 1~788~762~752.85(13)$~Hz, with improved uncertainties. Both constants are relevant to the evaluation of the $^{27}$Al$^{+}$ quantum-logic clock systematic uncertainty. The measurement of $C_{2}$ is in agreement with a previous measurement and a new calculation at the $1~σ$ level. The measurement of $ΔW$ is in good agreement with a recent measurement and differs from a previously published result by approximately $2σ$. With the improved value for $ΔW$, we deduce an improved value for the nuclear-to-electronic g-factor ratio $g_{I}/g_{J} = 9.299 ~308 ~313(60) \times 10^{-5}$ and the nuclear g-factor for the $^{25}$Mg nucleus $g_{I} = 1.861 ~957 ~82(28) \times 10^{-4}$. Using the values of $C_{2}$ and $ΔW$ presented here, we derive a quadratic Zeeman shift of the $^{27}$Al$^{+}$ quantum-logic clock of $Δν/ ν= -(9241.8 \pm 3.7) \times 10^{-19}$, for a bias magnetic field of $B \approx 0.12$~mT. △ Less

Submitted 22 May, 2019; v1 submitted 11 March, 2019; originally announced March 2019.

Comments: 7 pages, 5 figures, updated

Journal ref: Phys. Rev. A 100, 013409 (2019)

Optical-Clock-Based Time Scale

Authors: Jian Yao, Jeff A. Sherman, Tara Fortier, Holly Leopardi, Thomas Parker, William McGrew, Xiaogang Zhang, Daniele Nicolodi, Robert Fasano, Stefan Schäffer, Kyle Beloy, Joshua Savory, Stefania Romisch, Chris Oates, Scott Diddams, Andrew Ludlow, Judah Levine

Abstract: A time scale is a procedure for accurately and continuously marking the passage of time. It is exemplified by Coordinated Universal Time (UTC), and provides the backbone for critical navigation tools such as the Global Positioning System (GPS). Present time scales employ microwave atomic clocks, whose attributes can be combined and averaged in a manner such that the composite is more stable, accur… ▽ More A time scale is a procedure for accurately and continuously marking the passage of time. It is exemplified by Coordinated Universal Time (UTC), and provides the backbone for critical navigation tools such as the Global Positioning System (GPS). Present time scales employ microwave atomic clocks, whose attributes can be combined and averaged in a manner such that the composite is more stable, accurate, and reliable than the output of any individual clock. Over the past decade, clocks operating at optical frequencies have been introduced which are orders of magnitude more stable than any microwave clock. However, in spite of their great potential, these optical clocks cannot be operated continuously, which makes their use in a time scale problematic. In this paper, we report the development of a hybrid microwave-optical time scale, which only requires the optical clock to run intermittently while relying upon the ensemble of microwave clocks to serve as the flywheel oscillator. The benefit of using clock ensemble as the flywheel oscillator, instead of a single clock, can be understood by the Dick-effect limit. This time scale demonstrates for the first time sub-nanosecond accuracy for a few months, attaining a fractional frequency uncertainty of 1.45*10-16 at 30 days and reaching the 10-17 decade at 50 days, with respect to UTC. This time scale significantly improves the accuracy in timekeeping and could change the existing time-scale architectures. △ Less

Submitted 10 April, 2019; v1 submitted 18 February, 2019; originally announced February 2019.

Journal ref: Phys. Rev. Applied 12, 044069 (2019)

Galaxy And Mass Assembly (GAMA): Environmental Quenching of Centrals and Satellites in Groups

Authors: L. J. M. Davies, A. S. G. Robotham, C. del P. Lagos, S. P. Driver, A. R. H. Stevens, Y. M. Bahé, M. Alpaslan, M. N. Bremer, M. J. I. Brown, S. Brough, J. Bland-Hawthorn, L. Cortese, P. Elahi, M. W. Grootes, B. W. Holwerda, A. D. Ludlow, S. McGee, M. Owers, S. Phillipps

Abstract: Recently a number of studies have found a similarity between the passive fraction of central and satellite galaxies when controlled for both stellar and halo mass. These results suggest that the quenching processes that affect galaxies are largely agnostic to central/satellite status, which contradicts the traditional picture of increased satellite quenching via environmental processes such as str… ▽ More Recently a number of studies have found a similarity between the passive fraction of central and satellite galaxies when controlled for both stellar and halo mass. These results suggest that the quenching processes that affect galaxies are largely agnostic to central/satellite status, which contradicts the traditional picture of increased satellite quenching via environmental processes such as stripping, strangulation and starvation. Here we explore this further using the Galaxy And Mass Assembly (GAMA) survey which extends to ~2dex lower in stellar mass than SDSS, is more complete for closely-separated galaxies (>95% compared to >70%), and identifies lower-halo-mass groups outside of the very local Universe (M$_{\mathrm{halo}}\sim10^{12}$M$_{\odot}$ at $0.1<z<0.2$). As far as possible we aim to replicate the selections, completeness corrections and central/satellite division of one of the previous studies but find clear differences between passive fractions of centrals and satellites. We also find that our passive fractions increase with both halo-to-satellite mass ratio and central-to-second rank mass ratio. This suggests that quenching is more efficient in satellites that are low-mass for their halo ($i.e$ at high halo-to-satellite mass ratio in comparison to low halo-to-satellite mass ratio) and are more likely to be passive in older groups - forming a consistent picture of environmental quenching of satellites. We then discuss potential explanations for the previously observed similarity, such as dependence on the group-finding method. △ Less

Submitted 6 January, 2019; originally announced January 2019.

Comments: 15 pages, 8 figures, MNRAS accepted

Halo concentrations from extended Press-Schechter merger histories

Authors: Andrew J. Benson, Aaron Ludlow, Shaun Cole

Abstract: We apply the model relating halo concentration to formation history proposed by Ludlow et al. to merger trees generated using an algorithm based on excursion set theory. We find that while the model correctly predicts the median relation between halo concentration and mass, it underpredicts the scatter in concentration at fixed mass. Since the same model applied to N-body merger trees predicts the… ▽ More We apply the model relating halo concentration to formation history proposed by Ludlow et al. to merger trees generated using an algorithm based on excursion set theory. We find that while the model correctly predicts the median relation between halo concentration and mass, it underpredicts the scatter in concentration at fixed mass. Since the same model applied to N-body merger trees predicts the correct scatter, we postulate that the missing scatter is due to the lack of any environmental dependence in merger trees derived from excursion set theory. We show that a simple modification to the merger tree construction algorithm, which makes merger rates dependent on environment, can increase the scatter by the required amount, and simultaneously provide a qualitatively correct correlation between environment and formation epoch in the excursion set merger trees. △ Less

Submitted 14 December, 2018; originally announced December 2018.

Comments: 10 pages, 7 figures, submitted to MNRAS

Numerical convergence of simulations of galaxy formation: the abundance and internal structure of cold dark matter haloes

Authors: Aaron D. Ludlow, Joop Schaye, Richard Bower

Abstract: We study the impact of numerical parameters on the properties of cold dark matter haloes formed in collisionless cosmological simulations. We quantify convergence in the median spherically-averaged circular velocity profiles for haloes of widely varying particle number, as well as in the statistics of their structural scaling relations and mass functions. In agreement with prior work focused on si… ▽ More We study the impact of numerical parameters on the properties of cold dark matter haloes formed in collisionless cosmological simulations. We quantify convergence in the median spherically-averaged circular velocity profiles for haloes of widely varying particle number, as well as in the statistics of their structural scaling relations and mass functions. In agreement with prior work focused on single haloes, our results suggest that cosmological simulations yield robust halo properties for a wide range of gravitational softening parameters, $ε$, provided: 1) $ε$ is not larger than a "convergence radius", $r_{\rm conv}$, which is dictated by 2-body relaxation and determined by particle number, and 2) a sufficient number of timesteps are taken to accurately resolve particle orbits with short dynamical times. Provided these conditions are met, median circular velocity profiles converge to within $\approx 10$ per cent for radii beyond which the local 2-body relaxation timescale exceeds the Hubble time by a factor $κ\equiv t_{\rm relax}/t_{\rm H}\gt 0.177$, with better convergence attained for higher $κ$. We provide analytic estimates of $r_{\rm conv}$ that build on previous attempts in two ways: first, by highlighting its explicit (but weak) softening-dependence and, second, by providing a simpler criterion in which $r_{\rm conv}$ is determined entirely by the mean inter-particle spacing, $l$; for example, better than $10$ per cent convergence in circular velocity for $r\gt 0.05\,l$. We show how these analytic criteria can be used to assess convergence in structural scaling relations for dark matter haloes as a function of their mass or maximum circular speed. △ Less

Submitted 10 July, 2019; v1 submitted 13 December, 2018; originally announced December 2018.

Comments: 24 pages, 15 figures; MNRAS-accepted version

Towards Adoption of an Optical Second: Verifying Optical Clocks at the SI Limit

Authors: W. F. McGrew, X. Zhang, H. Leopardi, R. J. Fasano, D. Nicolodi, K. Beloy, J. Yao, J. A. Sherman, S. A. Schäffer, J. Savory, R. C. Brown, S. Römisch, C. W. Oates, T. E. Parker, T. M. Fortier, A. D. Ludlow

Abstract: The pursuit of ever more precise measures of time and frequency is likely to lead to the eventual redefinition of the second in terms of an optical atomic transition. To ensure continuity with the current definition, based on a microwave transition between hyperfine levels in ground-state $^{133}$Cs, it is necessary to measure the absolute frequency of candidate standards, which is done by compari… ▽ More The pursuit of ever more precise measures of time and frequency is likely to lead to the eventual redefinition of the second in terms of an optical atomic transition. To ensure continuity with the current definition, based on a microwave transition between hyperfine levels in ground-state $^{133}$Cs, it is necessary to measure the absolute frequency of candidate standards, which is done by comparing against a primary cesium reference. A key verification of this process can be achieved by performing a loop closure$-$comparing frequency ratios derived from absolute frequency measurements against ratios determined from direct optical comparisons. We measure the $^1$S$_0\!\rightarrow^3$P$_0$ transition of $^{171}$Yb by comparing the clock frequency to an international frequency standard with the aid of a maser ensemble serving as a flywheel oscillator. Our measurements consist of 79 separate runs spanning eight months, and we determine the absolute frequency to be 518 295 836 590 863.71(11) Hz, the uncertainty of which is equivalent to a fractional frequency of $2.1\times10^{-16}$. This absolute frequency measurement, the most accurate reported for any transition, allows us to close the Cs-Yb-Sr-Cs frequency measurement loop at an uncertainty of $<$3$\times10^{-16}$, limited by the current realization of the SI second. We use these measurements to tighten the constraints on variation of the electron-to-proton mass ratio, $μ=m_e/m_p$. Incorporating our measurements with the entire record of Yb and Sr absolute frequency measurements, we infer a coupling coefficient to gravitational potential of $k_\mathrmμ=(-1.9\pm 9.4)\times10^{-7}$ and a drift with respect to time of $\frac{\dotμ}μ=(5.3 \pm 6.5)\times10^{-17}/$yr. △ Less

Submitted 14 November, 2018; originally announced November 2018.

The signal of decaying dark matter with hydrodynamical simulations

Authors: Mark R. Lovell, David Barnes, Yannick Bahé, Joop Schaye, Matthieu Schaller, Tom Theuns, Sownak Bose, Robert A. Crain, Claudio dalla Vecchia, Carlos S. Frenk, Wojciech Hellwing, Scott T. Kay, Aaron D. Ludlow, Richard G. Bower

Abstract: Dark matter particles may decay, emitting photons. Drawing on the EAGLE family of hydrodynamic simulations of galaxy formation -- including the APOSTLE and C-EAGLE simulations -- we assess the systematic uncertainties and scatter on the decay flux from different galaxy classes, from Milky Way satellites to galaxy clusters, and compare our results to studies of the 3.55~keV line. We demonstrate tha… ▽ More Dark matter particles may decay, emitting photons. Drawing on the EAGLE family of hydrodynamic simulations of galaxy formation -- including the APOSTLE and C-EAGLE simulations -- we assess the systematic uncertainties and scatter on the decay flux from different galaxy classes, from Milky Way satellites to galaxy clusters, and compare our results to studies of the 3.55~keV line. We demonstrate that previous detections and non-detections of this line are consistent with a dark matter interpretation. For example, in our simulations the width of the the dark matter decay line for Perseus-analogue galaxy clusters lies in the range 1300-1700~\kms. Therefore, the non-detection of the 3.55~keV line in the centre of the Perseus cluster by the {\it Hitomi} collaboration is consistent with detections by other instruments. We also consider trends with stellar and halo mass and evaluate the scatter in the expected fluxes arising from the anisotropic halo mass distribution and from object-to-object variations. We provide specific predictions for observations with {\it XMM-Newton} and with the planned X-ray telescopes {\it XRISM} and {\it ATHENA}. If future detections of unexplained X-ray lines match our predictions, including line widths, we will have strong evidence that we have discovered the dark matter. △ Less

Submitted 11 October, 2018; originally announced October 2018.

Comments: 20 pages, 16 Figures. Highlights: Figs 4, 13, 15. To be submitted to MNRAS, comments welcome. Contact: lovell@hi.is

Baryon-induced dark matter cores in the EAGLE simulations

Authors: Alejandro Benitez-Llambay, Carlos S. Frenk, Aaron D. Ludlow, Julio F. Navarro

Abstract: We examine the formation of dark matter (DM) cores in dwarf galaxies simulated with the EAGLE model of galaxy formation. As in earlier work, we find that the star formation (SF) gas density threshold ($ρ_{\rm th}$) plays a critical role. At low thresholds (LT), gas is unable to reach densities high enough to dominate the gravitational potential before being dispersed by feedback from supernovae. L… ▽ More We examine the formation of dark matter (DM) cores in dwarf galaxies simulated with the EAGLE model of galaxy formation. As in earlier work, we find that the star formation (SF) gas density threshold ($ρ_{\rm th}$) plays a critical role. At low thresholds (LT), gas is unable to reach densities high enough to dominate the gravitational potential before being dispersed by feedback from supernovae. LT runs show little effect on the inner DM profile, even in systems with extended and bursty SF, two ingredients often cited as critical for core formation. For higher thresholds, gas is able to dominate the gravitational potential before being ejected by feedback. This can lead to a substantial reduction in the inner DM content, but only if the gas is gravitationally important over an extended period of time, allowing the halo to contract before gas removal. Rapid assembly and removal of gas in short SF bursts is less effective at altering the inner DM content. Subsequent gas accretion may draw DM back in and reform a cusp, unless SF is bursty enough to prevent it, preserving the core. Thus, for the EAGLE SF+feedback model, there is no simple relation between core formation and SF history, contrary to recent claims. The dependence of the inner DM content of dwarfs on $ρ_{\rm th}$ hinders robust predictions and the interpretation of observations. A simulation of a $(12 \rm \ Mpc)^3$ volume with high $ρ_{\rm th}$ results in dwarfs with sizeable cores over a limited halo mass range, but with insufficient variety in mass profiles to explain the observed diversity of dwarf galaxy rotation curves. △ Less

Submitted 23 August, 2019; v1 submitted 9 October, 2018; originally announced October 2018.

Comments: Updated to match the published version

Atomic clock performance beyond the geodetic limit

Authors: W. F. McGrew, X. Zhang, R. J. Fasano, S. A. Schäffer, K. Beloy, D. Nicolodi, R. C. Brown, N. Hinkley, G. Milani, M. Schioppo, T. H. Yoon, A. D. Ludlow

Abstract: The passage of time is tracked by counting oscillations of a frequency reference, such as Earth's revolutions or swings of a pendulum. By referencing atomic transitions, frequency (and thus time) can be measured more precisely than any other physical quantity, with the current generation of optical atomic clocks reporting fractional performance below the $10^{-17}$ level. However, the theory of re… ▽ More The passage of time is tracked by counting oscillations of a frequency reference, such as Earth's revolutions or swings of a pendulum. By referencing atomic transitions, frequency (and thus time) can be measured more precisely than any other physical quantity, with the current generation of optical atomic clocks reporting fractional performance below the $10^{-17}$ level. However, the theory of relativity prescribes that the passage of time is not absolute, but impacted by an observer's reference frame. Consequently, clock measurements exhibit sensitivity to relative velocity, acceleration and gravity potential. Here we demonstrate optical clock measurements surpassing the present-day ability to account for the gravitational distortion of space-time across the surface of Earth. In two independent ytterbium optical lattice clocks, we demonstrate unprecedented levels in three fundamental benchmarks of clock performance. In units of the clock frequency, we report systematic uncertainty of $1.4 \times 10^{-18}$, measurement instability of $3.2 \times 10^{-19}$ and reproducibility characterised by ten blinded frequency comparisons, yielding a frequency difference of $[-7 \pm (5)_{stat} \pm (8)_{sys}] \times 10^{-19}$. While differential sensitivity to gravity could degrade the performance of these optical clocks as terrestrial standards of time, this same sensitivity can be used as an exquisite probe of geopotential. Near the surface of Earth, clock comparisons at the $1 \times 10^{-18}$ level provide 1 cm resolution along gravity, outperforming state-of-the-art geodetic techniques. These optical clocks can further be used to explore geophysical phenomena, detect gravitational waves, test general relativity and search for dark matter. △ Less

Submitted 30 July, 2018; originally announced July 2018.

Comments: 32 pages, 5 figures

Journal ref: Nature 564, 87-90 (2018)

First observation with global network of optical atomic clocks aimed for a dark matter detection

Authors: P. Wcisło, P. Ablewski, K. Beloy, S. Bilicki, M. Bober, R. Brown, R. Fasano, R. Ciuryło, H. Hachisu, T. Ido, J. Lodewyck, A. Ludlow, W. McGrew, P. Morzyński, D. Nicolodi, M. Schioppo, M. Sekido, R. Le Targat, P. Wolf, X. Zhang, B. Zjawin, M. Zawada

Abstract: We report on the first earth-scale quantum sensor network based on optical atomic clocks aimed at dark matter (DM) detection. Exploiting differences in the susceptibilities to the fine-structure constant of essential parts of an optical atomic clock, i.e. the cold atoms and the optical reference cavity, we can perform sensitive searches for dark matter signatures without the need of real-time comp… ▽ More We report on the first earth-scale quantum sensor network based on optical atomic clocks aimed at dark matter (DM) detection. Exploiting differences in the susceptibilities to the fine-structure constant of essential parts of an optical atomic clock, i.e. the cold atoms and the optical reference cavity, we can perform sensitive searches for dark matter signatures without the need of real-time comparisons of the clocks. We report a two orders of magnitude improvement in constraints on transient variations of the fine-structure constant, which considerably improves the detection limit for the standard model (SM) - DM coupling. We use Yb and Sr optical atomic clocks at four laboratories on three continents to search for both topological defect (TD) and massive scalar field candidates. No signal consistent with a dark-matter coupling is identified, leading to significantly improved constraints on the DM-SM couplings. △ Less

Submitted 12 June, 2018; originally announced June 2018.

Journal ref: Science Advances, 4, eaau4869 (2018)

Faraday-shielded, DC Stark-free optical lattice clock

Authors: K. Beloy, X. Zhang, W. F. McGrew, N. Hinkley, T. H. Yoon, D. Nicolodi, R. J. Fasano, S. A. Schäffer, R. C. Brown, A. D. Ludlow

Abstract: We demonstrate the absence of a DC Stark shift in an ytterbium optical lattice clock. Stray electric fields are suppressed through the introduction of an in-vacuum Faraday shield. Still, the effectiveness of the shielding must be experimentally assessed. Such diagnostics are accomplished by applying high voltage to six electrodes, which are grounded in normal operation to form part of the Faraday… ▽ More We demonstrate the absence of a DC Stark shift in an ytterbium optical lattice clock. Stray electric fields are suppressed through the introduction of an in-vacuum Faraday shield. Still, the effectiveness of the shielding must be experimentally assessed. Such diagnostics are accomplished by applying high voltage to six electrodes, which are grounded in normal operation to form part of the Faraday shield. Our measurements place a constraint on the DC Stark shift at the $10^{-20}$ level, in units of the clock frequency. Moreover, we discuss a potential source of error in strategies to precisely measure or cancel non-zero DC Stark shifts, attributed to field gradients coupled with the finite spatial extent of the lattice-trapped atoms. With this consideration, we find that Faraday shielding, complemented with experimental validation, provides both a practically appealing and effective solution to the problem of DC Stark shifts in optical lattice clocks. △ Less

Submitted 28 March, 2018; originally announced March 2018.

Comments: 5 pages + supplemental material; accepted to PRL

Journal ref: Phys. Rev. Lett. 120, 183201 (2018)

Hyperpolarizability and operational magic wavelength in an optical lattice clock

Authors: R. C. Brown, N. B. Phillips, K. Beloy, W. F. McGrew, M. Schioppo, R. J. Fasano, G. Milani, X. Zhang, N. Hinkley, H. Leopardi, T. H. Yoon, D. Nicolodi, T. M. Fortier, A. D. Ludlow

Abstract: Optical clocks benefit from tight atomic confinement enabling extended interrogation times as well as Doppler- and recoil-free operation. However, these benefits come at the cost of frequency shifts that, if not properly controlled, may degrade clock accuracy. Numerous theoretical studies have predicted optical lattice clock frequency shifts that scale nonlinearly with trap depth. To experimentall… ▽ More Optical clocks benefit from tight atomic confinement enabling extended interrogation times as well as Doppler- and recoil-free operation. However, these benefits come at the cost of frequency shifts that, if not properly controlled, may degrade clock accuracy. Numerous theoretical studies have predicted optical lattice clock frequency shifts that scale nonlinearly with trap depth. To experimentally observe and constrain these shifts in an $^{171}$Yb optical lattice clock, we construct a lattice enhancement cavity that exaggerates the light shifts. We observe an atomic temperature that is proportional to the optical trap depth, fundamentally altering the scaling of trap-induced light shifts and simplifying their parametrization. We identify an "operational" magic wavelength where frequency shifts are insensitive to changes in trap depth. These measurements and scaling analysis constitute an essential systematic characterization for clock operation at the $10^{-18}$ level and beyond. △ Less

Submitted 24 October, 2017; v1 submitted 29 August, 2017; originally announced August 2017.

Comments: 5 + 2 pages, 3 figures, added supplemental

Journal ref: Phys. Rev. Lett. 119, 253001 (2017)

The vertical structure of gaseous galaxy discs in cold dark matter halos

Authors: Alejandro Benitez-Llambay, Julio F. Navarro, Carlos S. Frenk, Aaron D. Ludlow

Abstract: We study the vertical structure of polytropic, $P\propto ρ^Γ$, centrifugally-supported gaseous discs embedded in cold dark matter (CDM) halos. At fixed radius $R$, the shape of the vertical density profile depends only weakly on whether the disc is self-gravitating (SG) or not (NSG). The disc thickness, set by the midplane sound speed and circular velocity, $(c_s/V_c)R$, in the NSG case, and by th… ▽ More We study the vertical structure of polytropic, $P\propto ρ^Γ$, centrifugally-supported gaseous discs embedded in cold dark matter (CDM) halos. At fixed radius $R$, the shape of the vertical density profile depends only weakly on whether the disc is self-gravitating (SG) or not (NSG). The disc thickness, set by the midplane sound speed and circular velocity, $(c_s/V_c)R$, in the NSG case, and by the sound speed and surface density, $c_s^2/GΣ$, in SG discs, is smaller than either of these scales. SG discs are typically Toomre unstable, NSG discs are stable. Exponential discs in CDM halos with roughly flat circular velocity curves generally "flare" outwards. For the polytropic equation of state of the EAGLE simulations, discs whose mass and size match observational constraints are stable (NSG) for $M_d< 3\times 10^9\, M_\odot$ and unstable (SG) at higher masses, if fully gaseous. We test these analytic results using a set of idealized SPH simulations and find excellent agreement. Our results clarify the role of the gravitational softening on the thickness of simulated discs, and on the onset of radial instabilities. EAGLE low-mass discs are non-self-gravitating so the softening plays no role in their vertical structure. High-mass discs, on the other hand, are expected to be self-gravitating and unstable, and may be artificially thickened and stabilized unless gravity is well resolved. Simulations with spatial resolution high enough to not compromise the vertical structure of a disc also resolve the onset of their instabilities, but the converse is not true: resolving instabilities does not guarantee that the vertical structure is resolved. △ Less

Submitted 25 July, 2017; originally announced July 2017.

Comments: Submitted to MNRAS

Photonic-chip supercontinuum with tailored spectra for precision frequency metrology

Authors: David Carlson, Daniel Hickstein, Alexander Lind, Judith Olson, Richard Fox, Roger Brown, Andrew Ludlow, Qing Li, Daron Westly, Holly Leopardi, Tara Fortier, Kartik Srinivasan, Scott Diddams, Scott Papp

Abstract: Supercontinuum generation using chip-integrated photonic waveguides is a powerful approach for spectrally broadening pulsed laser sources with very low pulse energies and compact form factors. When pumped with a mode-locked laser frequency comb, these waveguides can coherently expand the comb spectrum to more than an octave in bandwidth to enable self-referenced stabilization. However, for applica… ▽ More Supercontinuum generation using chip-integrated photonic waveguides is a powerful approach for spectrally broadening pulsed laser sources with very low pulse energies and compact form factors. When pumped with a mode-locked laser frequency comb, these waveguides can coherently expand the comb spectrum to more than an octave in bandwidth to enable self-referenced stabilization. However, for applications in frequency metrology and precision spectroscopy, it is desirable to not only support self-referencing, but also to generate low-noise combs with customizable broadband spectra. In this work, we demonstrate dispersion-engineered waveguides based on silicon nitride that are designed to meet these goals and enable precision optical metrology experiments across large wavelength spans. We perform a clock comparison measurement and report a clock-limited relative frequency instability of $3.8\times10^{-15}$ at $τ= 2$ seconds between a 1550 nm cavity-stabilized reference laser and NIST's calcium atomic clock laser at 657 nm using a two-octave waveguide-supercontinuum comb. △ Less

Submitted 13 February, 2017; v1 submitted 10 February, 2017; originally announced February 2017.

Journal ref: Phys. Rev. Applied 8, 014027 (2017)

A compact, thermal noise limited reference cavity for ultra-low noise microwave generation

Authors: J. Davila-Rodriguez, F. N. Baynes, A. D. Ludlow, T. M. Fortier, H. Leopardi, S. A. Diddams, F. Quinlan

Abstract: We demonstrate an easy to manufacture, 25 mm long ultra-stable optical reference cavity for transportable photonic microwave generation systems. Employing a rigid holding geometry that is first-order insensitive to the squeezing force and a cavity geometry that improves the thermal noise limit at room temperature, we observe a laser phase noise that is nearly thermal noise limited for three freque… ▽ More We demonstrate an easy to manufacture, 25 mm long ultra-stable optical reference cavity for transportable photonic microwave generation systems. Employing a rigid holding geometry that is first-order insensitive to the squeezing force and a cavity geometry that improves the thermal noise limit at room temperature, we observe a laser phase noise that is nearly thermal noise limited for three frequency decades (1 Hz to 1 kHz offset) and supports 10 GHz generation with phase noise near -100 dBc/Hz at 1 Hz offset and <-173 dBc/Hz for all offsets >600 Hz. The fractional frequency stability reaches $2\times10^{-15}$ at 0.1 s of averaging. △ Less

Submitted 6 January, 2017; originally announced January 2017.

Comments: Letter

The origin of the mass discrepancy-acceleration relation in $Λ$CDM

Authors: Julio F. Navarro, Alejandro Benítez-Llambay, Azadeh Fattahi, Carlos S. Frenk, Aaron D. Ludlow, Kyle A. Oman, Matthieu Schaller, Tom Theuns

Abstract: We examine the origin of the mass discrepancy--radial acceleration relation (MDAR) of disk galaxies. This is a tight empirical correlation between the disk centripetal acceleration and that expected from the baryonic component. The MDAR holds for most radii probed by disk kinematic tracers, regardless of galaxy mass or surface brightness. The relation has two characteristic accelerations; $a_0$, a… ▽ More We examine the origin of the mass discrepancy--radial acceleration relation (MDAR) of disk galaxies. This is a tight empirical correlation between the disk centripetal acceleration and that expected from the baryonic component. The MDAR holds for most radii probed by disk kinematic tracers, regardless of galaxy mass or surface brightness. The relation has two characteristic accelerations; $a_0$, above which all galaxies are baryon-dominated; and $a_{\rm min}$, an effective minimum aceleration probed by kinematic tracers in isolated galaxies. We use a simple model to show that these trends arise naturally in $Λ$CDM. This is because: (i) disk galaxies in $Λ$CDM form at the centre of dark matter haloes spanning a relatively narrow range of virial mass; (ii) cold dark matter halo acceleration profiles are self-similar and have a broad maximum at the centre, reaching values bracketed precisely by $a_{\rm min}$ and $a_0$ in that mass range; and (iii) halo mass and galaxy size scale relatively tightly with the baryonic mass of a galaxy in any successful $Λ$CDM galaxy formation model. Explaining the MDAR in $Λ$CDM does not require modifications to the cuspy inner mass profiles of dark haloes, although these may help to understand the detailed rotation curves of some dwarf galaxies and the origin of extreme outliers from the main relation. The MDAR is just a reflection of the self-similar nature of cold dark matter haloes and of the physical scales introduced by the galaxy formation process. △ Less

Submitted 26 June, 2017; v1 submitted 19 December, 2016; originally announced December 2016.

Comments: 8 pages, 3 figures, MNRAS submitted