-
Dynamical Friction and Black Holes in Ultralight Dark Matter Solitons
Authors:
Russell Boey,
Yourong Wang,
Emily Kendall,
Richard Easther
Abstract:
We numerically simulate the motion of a black hole as it plunges radially through an ultralight dark matter soliton. We investigate the timescale in which dynamical friction reduces the kinetic energy of the black hole to a minimum, and consider the sensitivity of this timescale to changes in the ULDM particle mass, the total soliton mass, and the mass of the black hole. We contrast our numerical…
▽ More
We numerically simulate the motion of a black hole as it plunges radially through an ultralight dark matter soliton. We investigate the timescale in which dynamical friction reduces the kinetic energy of the black hole to a minimum, and consider the sensitivity of this timescale to changes in the ULDM particle mass, the total soliton mass, and the mass of the black hole. We contrast our numerical results with a semi-analytic treatment of dynamical friction, and find that the latter is poorly suited to this scenario. In particular, we find that the back-reaction of the soliton to the presence of the black hole is significant, resulting in oscillations in the coefficient of dynamical friction which cannot be described in the simple semi-analytical framework. Furthermore, we observe a late-time reheating effect, in which a significant amount of kinetic energy is transferred back to the black hole after an initial damping phase. This complicates the discussion of ULDM dynamical friction on the scales relevant to the final parsec problem.
△ Less
Submitted 13 March, 2024;
originally announced March 2024.
-
Aspherical ULDM Collapse: Variation in the Core-Halo Mass Relation
Authors:
Emily Kendall,
Mateja Gosenca,
Richard Easther
Abstract:
Ultralight dark matter (ULDM) is an interesting alternative to the cold dark matter (CDM) paradigm. Due to the extremely low mass of the constituent particle ($\sim 10^{-22}$ eV), ULDM can exhibit quantum effects up to kiloparsec scales. In particular, runaway collapse in the centres of ULDM halos is prevented by quantum pressure, providing a possible resolution to the 'core-cusp problem' of CDM.…
▽ More
Ultralight dark matter (ULDM) is an interesting alternative to the cold dark matter (CDM) paradigm. Due to the extremely low mass of the constituent particle ($\sim 10^{-22}$ eV), ULDM can exhibit quantum effects up to kiloparsec scales. In particular, runaway collapse in the centres of ULDM halos is prevented by quantum pressure, providing a possible resolution to the 'core-cusp problem' of CDM. However, the the detailed relationship between the ULDM core mass and that of the overall halo is poorly understood. We simulate the collapse of both spherical and aspherical isolated ULDM overdensities using AxioNyx, finding that the central cores of collapsed halos undergo sustained oscillatory behaviour which affects both their peak density and overall morphology. The variability in core morphology increases with the asphericity of the initial overdensity and remnants of initial asphericity persist long after collapse. Furthermore, the peak central densities are higher in spherical configurations. Consequently, astrophysically realistic halos may exhibit substantial departures from theoretical core-halo profiles and we would expect a significant variance of the properties of halos with the same mass.
△ Less
Submitted 17 May, 2023;
originally announced May 2023.
-
Commuter Count: Inferring Travel Patterns from Location Data
Authors:
Nathan Musoke,
Emily Kendall,
Mateja Gosenca,
Lillian Guo,
Lerh Feng Low,
Angela Xue,
Richard Easther
Abstract:
In this Working Paper we analyse computational strategies for using aggregated spatio-temporal population data acquired from telecommunications networks to infer travel and movement patterns between geographical regions. Specifically, we focus on hour-by-hour cellphone counts for the SA-2 geographical regions covering the whole of New Zealand. This Working Paper describes the implementation of the…
▽ More
In this Working Paper we analyse computational strategies for using aggregated spatio-temporal population data acquired from telecommunications networks to infer travel and movement patterns between geographical regions. Specifically, we focus on hour-by-hour cellphone counts for the SA-2 geographical regions covering the whole of New Zealand. This Working Paper describes the implementation of the inference algorithms, their ability to produce models of travel patterns during the day, and lays out opportunities for future development.
△ Less
Submitted 30 March, 2023;
originally announced March 2023.
-
Multifield Ultralight Dark Matter
Authors:
Mateja Gosenca,
Andrew Eberhardt,
Yourong Wang,
Benedikt Eggemeier,
Emily Kendall,
J. Luna Zagorac,
Richard Easther
Abstract:
Ultralight dark matter (ULDM) is usually taken to be a single scalar field. Here we explore the possibility that ULDM consists of $N$ light scalar fields with only gravitational interactions. This configuration is more consistent with the underlying particle physics motivations for these scenarios than a single ultralight field. ULDM halos have a characteristic granular structure that increases st…
▽ More
Ultralight dark matter (ULDM) is usually taken to be a single scalar field. Here we explore the possibility that ULDM consists of $N$ light scalar fields with only gravitational interactions. This configuration is more consistent with the underlying particle physics motivations for these scenarios than a single ultralight field. ULDM halos have a characteristic granular structure that increases stellar velocity dispersion and can be used as observational constraints on ULDM models. In multifield simulations, we find that inside a halo the amplitude of the total density fluctuations decreases as $1/\sqrt{N}$ and that the fields do not become significantly correlated over cosmological timescales. Smoother halos heat stellar orbits less efficiently, reducing the velocity dispersion relative to the single field case and thus weakening the observational constraints on the field mass. Analytically, we show that for $N$ equal-mass fields with mass $m$ the ULDM contribution to the stellar velocity dispersion scales as $1/(N m^3)$. Lighter fields heat the most efficiently and if the smallest mass $m_L$ is significantly below the other field masses the dispersion scales as $1/(N^2 m_L^3)$.
△ Less
Submitted 17 January, 2023;
originally announced January 2023.
-
Soliton Formation and the Core-Halo Mass Relation: An Eigenstate Perspective
Authors:
J. Luna Zagorac,
Emily Kendall,
Nikhil Padmanabhan,
Richard Easther
Abstract:
UltraLight Dark Matter (ULDM) is an axion-like dark matter candidate with an extremely small particle mass. ULDM halos consist of a spherically symmetric solitonic core and an NFW-like skirt. We simulate halo creation via soliton mergers and use these results to explore the core-halo mass relation. We calculate the eigenstates of the merged halos and use these to isolate the solitonic core and cal…
▽ More
UltraLight Dark Matter (ULDM) is an axion-like dark matter candidate with an extremely small particle mass. ULDM halos consist of a spherically symmetric solitonic core and an NFW-like skirt. We simulate halo creation via soliton mergers and use these results to explore the core-halo mass relation. We calculate the eigenstates of the merged halos and use these to isolate the solitonic core and calculate its relative contribution to the halo mass. We compare this approach to using a fitting function to isolate the core and find a difference in masses up to 30%. We analyze three families of simulations: equal-mass mergers, unequal-mass mergers, and halos with a two-step merger history. Setting the halo mass to the initial mass in the simulation does not yield a consistent core-halo relationship. Excluding material "ejected" by the collision yields a core-halo relationship with a slope of 1/3 for simultaneous mergers and roughly 0.4 for two-step mergers. Our findings suggest there is no universal core-halo mass relationship for ULDM and shed light on the differing results for the core-halo relationship previously reported in the literature.
△ Less
Submitted 28 April, 2023; v1 submitted 19 December, 2022;
originally announced December 2022.
-
Beyond time-homogeneity for continuous-time multistate Markov models
Authors:
Emmett B. Kendall,
Jonathan P. Williams,
Gudmund H. Hermansen,
Frederic Bois,
Vo Hong Thanh
Abstract:
Multistate Markov models are a canonical parametric approach for data modeling of observed or latent stochastic processes supported on a finite state space. Continuous-time Markov processes describe data that are observed irregularly over time, as is often the case in longitudinal medical data, for example. Assuming that a continuous-time Markov process is time-homogeneous, a closed-form likelihoo…
▽ More
Multistate Markov models are a canonical parametric approach for data modeling of observed or latent stochastic processes supported on a finite state space. Continuous-time Markov processes describe data that are observed irregularly over time, as is often the case in longitudinal medical data, for example. Assuming that a continuous-time Markov process is time-homogeneous, a closed-form likelihood function can be derived from the Kolmogorov forward equations -- a system of differential equations with a well-known matrix-exponential solution. Unfortunately, however, the forward equations do not admit an analytical solution for continuous-time, time-inhomogeneous Markov processes, and so researchers and practitioners often make the simplifying assumption that the process is piecewise time-homogeneous. In this paper, we provide intuitions and illustrations of the potential biases for parameter estimation that may ensue in the more realistic scenario that the piecewise-homogeneous assumption is violated, and we advocate for a solution for likelihood computation in a truly time-inhomogeneous fashion. Particular focus is afforded to the context of multistate Markov models that allow for state label misclassifications, which applies more broadly to hidden Markov models (HMMs), and Bayesian computations bypass the necessity for computationally demanding numerical gradient approximations for obtaining maximum likelihood estimates (MLEs). Supplemental materials are available online.
△ Less
Submitted 3 December, 2023; v1 submitted 6 November, 2022;
originally announced November 2022.
-
Robust Coordinate Ascent Variational Inference with Markov chain Monte Carlo simulations
Authors:
Neil Dey,
Emmett B. Kendall
Abstract:
Variational Inference (VI) is a method that approximates a difficult-to-compute posterior density using better behaved distributional families. VI is an alternative to the already well-studied Markov chain Monte Carlo (MCMC) method of approximating densities. With each algorithm, there are of course benefits and drawbacks; does there exist a combination of the two that mitigates the flaws of both?…
▽ More
Variational Inference (VI) is a method that approximates a difficult-to-compute posterior density using better behaved distributional families. VI is an alternative to the already well-studied Markov chain Monte Carlo (MCMC) method of approximating densities. With each algorithm, there are of course benefits and drawbacks; does there exist a combination of the two that mitigates the flaws of both? We propose a method to combine Coordinate Ascent Variational Inference (CAVI) with MCMC. This new methodology, termed Hybrid CAVI, seeks to improve the sensitivity to initialization and convergence problems of CAVI by proposing an initialization using method of moments estimates obtained from a short MCMC burn-in period. Unlike CAVI, Hybrid CAVI proves to also be effective when the posterior is not from a conditionally conjugate exponential family.
△ Less
Submitted 23 March, 2022;
originally announced March 2022.
-
Transmission of coherent information at the onset of interactions
Authors:
Emily Kendall,
Barbara Šoda,
Achim Kempf
Abstract:
In this work, we investigate the parameters governing the rate at which a quantum channel arises at the onset of an interaction between two systems, $A$ and $B$. In particular, when system $A$ is pre-entangled with an ancilla, $\tilde{A}$, we quantify the early-time transmission of pre-existing entanglement by calculating the leading order change in coherent information of the complementary channe…
▽ More
In this work, we investigate the parameters governing the rate at which a quantum channel arises at the onset of an interaction between two systems, $A$ and $B$. In particular, when system $A$ is pre-entangled with an ancilla, $\tilde{A}$, we quantify the early-time transmission of pre-existing entanglement by calculating the leading order change in coherent information of the complementary channel ($A\rightarrow B'$). We show that, when $A$ and $B$ are initially unentangled and $B$ is pure, there is no change in coherent information to first order, while the leading (second) order change is divergent. However, this divergence may be regulated by embedding the conventional notion of coherent information into what we call the family of $n$-coherent informations, defined using $n$-Rényi entropies. We find that the rate of change of the $n$-coherent information at the onset of the interaction is governed by a quantity, which we call the $n$-exposure, which captures the extent to which the initial coherent information of $A$ with $\tilde{A}$ is exposed to or `seen by' the interaction Hamiltonian between $A$ and $B$. We give examples in qubit systems and in the light-matter interaction.
△ Less
Submitted 5 May, 2022; v1 submitted 13 October, 2021;
originally announced October 2021.
-
Robust inference for geographic regression discontinuity designs: assessing the impact of police precincts
Authors:
Emmett B. Kendall,
Brenden Beck,
Joseph Antonelli
Abstract:
We study variation in policing outcomes attributable to differential policing practices in New York City (NYC) using geographic regression discontinuity designs (GeoRDDs). By focusing on small geographic windows near police precinct boundaries we can estimate local average treatment effects of police precincts on arrest rates. We propose estimands and develop estimators for the GeoRDD when the dat…
▽ More
We study variation in policing outcomes attributable to differential policing practices in New York City (NYC) using geographic regression discontinuity designs (GeoRDDs). By focusing on small geographic windows near police precinct boundaries we can estimate local average treatment effects of police precincts on arrest rates. We propose estimands and develop estimators for the GeoRDD when the data come from a spatial point process. Additionally, standard GeoRDDs rely on continuity assumptions of the potential outcome surface or a local randomization assumption within a window around the boundary. These assumptions, however, can easily be violated in realistic applications. We develop a novel and robust approach to testing whether there are differences in policing outcomes that are caused by differences in police precincts across NYC. Importantly, this approach is applicable to standard regression discontinuity designs with both numeric and point process data. This approach is robust to violations of traditional assumptions made, and is valid under weaker assumptions. We use a unique form of resampling to provide a valid estimate of our test statistic's null distribution even under violations of standard assumptions. This procedure gives substantially different results in the analysis of NYC arrest rates than those that rely on standard assumptions.
△ Less
Submitted 3 February, 2024; v1 submitted 30 June, 2021;
originally announced June 2021.
-
Simple numerical algorithm for generating Hamiltonian cycles and edge labels on planar cubic maps
Authors:
Emily Kendall
Abstract:
In this work we present an algorithm with which any arbitrary cubic planar map may be constructed through successive edge insertion while simultaneously constructing a set of proper edge labels and Hamiltonian cycles for each configuration. We present a publicly available Python implementation of this algorithm, and discuss both theoretical and numerical support for its validity, with reference to…
▽ More
In this work we present an algorithm with which any arbitrary cubic planar map may be constructed through successive edge insertion while simultaneously constructing a set of proper edge labels and Hamiltonian cycles for each configuration. We present a publicly available Python implementation of this algorithm, and discuss both theoretical and numerical support for its validity, with reference to the well-known Four Colour Theorem.
△ Less
Submitted 12 April, 2021;
originally announced April 2021.
-
The Dynamics of Entropies at the Onset of Interactions
Authors:
Emily Kendall,
Achim Kempf
Abstract:
At the onset of an interaction between two initially independent systems, each system tends to experience an increase in its n-Renyi entropies, such as its von Neumann entropy (n = 1) and its mixedness (n = 2). We here ask which properties of a system determine how quickly its Renyi entropies increase and, therefore, how sensitive the system is to becoming entangled. We find that the rate at which…
▽ More
At the onset of an interaction between two initially independent systems, each system tends to experience an increase in its n-Renyi entropies, such as its von Neumann entropy (n = 1) and its mixedness (n = 2). We here ask which properties of a system determine how quickly its Renyi entropies increase and, therefore, how sensitive the system is to becoming entangled. We find that the rate at which the n-Renyi entropy increases in an interaction is determined by a quantity which we term the n-fragility of the system. The 2-fragility is closely related to the notion of 2-norm coherence, in that it too quantifies the extent to which a density matrix is off-diagonal with respect to the eigenbasis of a reference operator. Nevertheless, the 2-fragility is not a coherence monotone in the resource theoretic sense since it depends also on the eigenvalues of the reference operator. It is this additional sensitivity to the eigenvalues of the reference operator, here the interaction Hamiltonian, which enables the 2-fragility to quantify the rate of entropy production in interactions. We give an example using the light-matter interaction and we anticipate applications to the study of the rates at which two systems exchange classical and quantum information when starting to interact.
△ Less
Submitted 6 September, 2020; v1 submitted 6 April, 2020;
originally announced April 2020.
-
The Core-Cusp Problem Revisited: ULDM vs. CDM
Authors:
Emily Kendall,
Richard Easther
Abstract:
The core-cusp problem is a widely cited motivation for the exploration of dark matter models beyond standard CDM. One such alternative is ULDM; extremely light scalar particles exhibiting wavelike properties on kiloparsec scales. Astrophysically realistic ULDM halos are expected to consist of inner solitonic cores embedded in NFW-like outer halos. The presence of the solitonic core suggests that U…
▽ More
The core-cusp problem is a widely cited motivation for the exploration of dark matter models beyond standard CDM. One such alternative is ULDM; extremely light scalar particles exhibiting wavelike properties on kiloparsec scales. Astrophysically realistic ULDM halos are expected to consist of inner solitonic cores embedded in NFW-like outer halos. The presence of the solitonic core suggests that ULDM may resolve the core-cusp discrepancy associated with pure NFW halos without recourse to baryonic physics. However, it has been demonstrated that the density of ULDM halos can exceed those of comparable NFW configurations at some radii and halo masses, apparently exacerbating the problem rather than solving it. This situation arises because, although solitonic cores are flat at their centres, they obey an inverse mass-radius scaling relationship. Meanwhile, the mass of the inner soliton increases with the total halo mass, and therefore the inner core becomes more peaked at large halo masses. We describe a parameterisation of the radial density profiles of ULDM halos that allows for environmental variability of the core-halo mass relation in order to investigate this issue in more detail. For halos up to 10^12 solar masses we find feasible ULDM profiles for which the central density is lower than their NFW counterparts at astrophysically accessible radii. However, comparisons to observed profiles do not strongly favour either option; both give reasonable fits to subsets of the data for some parameter choices. Consequently, we find that robust tests of the core-cusp problem in ULDM will require more comprehensive observational data and simulations that include baryonic feedback.
△ Less
Submitted 27 January, 2020; v1 submitted 7 August, 2019;
originally announced August 2019.
-
PyUltraLight: A Pseudo-Spectral Solver for Ultralight Dark Matter Dynamics
Authors:
Faber Edwards,
Emily Kendall,
Shaun Hotchkiss,
Richard Easther
Abstract:
PyUltraLight simulates the dynamics of ultralight dark matter in a non-expanding background. PyUltraLight can describe the evolution of several interacting ultralight dark matter halos or one or more halos orbiting a central, fixed Newtonian potential, the latter scenario corresponding to dwarf galaxies orbiting a massive central galaxy. We verify PyUltraLight by showing that it reproduces qualita…
▽ More
PyUltraLight simulates the dynamics of ultralight dark matter in a non-expanding background. PyUltraLight can describe the evolution of several interacting ultralight dark matter halos or one or more halos orbiting a central, fixed Newtonian potential, the latter scenario corresponding to dwarf galaxies orbiting a massive central galaxy. We verify PyUltraLight by showing that it reproduces qualitative dynamical features of previously published simulations and demonstrate that it has excellent energy-conservation properties. PyUltraLight is implemented in a Python-based Jupyter notebook, solving the Schrodinger-Poisson equation governing ultralight scalar field dark matter dynamics in the non-relativistic regime using a symmetrised split-step pseudospectral algorithm. The notebook interface makes it simple to specify simulation parameters and visualise the resulting output but performance-critical routines are managed via calls to computationally efficient compiled libraries. PyUltraLight runs on standard desktop hardware with support for shared memory mutlithreading and is available on GitHub.
△ Less
Submitted 16 October, 2018; v1 submitted 11 July, 2018;
originally announced July 2018.