Researcher profile

Andrew Pontzen

Andrew Pontzen contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
15works
0followers
10topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2024arXiv

Analog vacuum decay from vacuum initial conditions

Ultracold atomic gases can undergo phase transitions that mimic relativistic vacuum decay, allowing us to empirically test early-Universe physics in tabletop experiments. We investigate the physics of these analog systems, going beyond previous analyses of the classical equations of motion to study quantum fluctuations in the cold-atom false vacuum. We show that the fluctuation spectrum of this vacuum state agrees with the usual relativistic result in the regime where the classical analogy holds, providing further evidence for the suitability of these systems for studying vacuum decay. Using a suite of semiclassical lattice simulations, we simulate bubble nucleation from this analog vacuum state in a 1D homonuclear potassium-41 mixture, finding qualitative agreement with instanton predictions. We identify realistic parameters for this system that will allow us to study vacuum decay with current experimental capabilities, including a prescription for efficiently scanning over decay rates, and show that this setup will probe the quantum (rather than thermal) decay regime at temperatures $T\lesssim10\,\mathrm{nK}$. Our results help lay the groundwork for using upcoming cold-atom experiments as a new probe of nonperturbative early-Universe physics.

preprint2024arXiv

Towards Accurate Field-Level Inference of Massive Cosmic Structures

We investigate the accuracy requirements for field-level inference of cluster and void masses using data from galaxy surveys. We introduce a two-step framework that takes advantage of the fact that cluster masses are determined by flows on larger scales than the clusters themselves. First, we determine the integration accuracy required to perform field-level inference of cosmic initial conditions on these large scales, by fitting to late-time galaxy counts using the Bayesian Origin Reconstruction from Galaxies (BORG) algorithm. A 20-step COLA integrator is able to accurately describe the density field surrounding the most massive clusters in the Local Super-Volume ($<135\,h^{-1}\,\mathrm{Mpc}$), but does not by itself lead to converged virial mass estimates. Therefore we carry out `posterior resimulations&#39;, using full $N$-body dynamics while sampling from the inferred initial conditions, and thereby obtain estimates of masses for nearby massive clusters. We show that these are in broad agreement with existing estimates, and find that mass functions in the Local Super-Volume are compatible with $Λ$CDM.

preprint2022arXiv

Discovering the building blocks of dark matter halo density profiles with neural networks

The density profiles of dark matter halos are typically modeled using empirical formulae fitted to the density profiles of relaxed halo populations. We present a neural network model that is trained to learn the mapping from the raw density field containing each halo to the dark matter density profile. We show that the model recovers the widely-used Navarro-Frenk-White (NFW) profile out to the virial radius, and can additionally describe the variability in the outer profile of the halos. The neural network architecture consists of a supervised encoder-decoder framework, which first compresses the density inputs into a low-dimensional latent representation, and then outputs $ρ(r)$ for any desired value of radius $r$. The latent representation contains all the information used by the model to predict the density profiles. This allows us to interpret the latent representation by quantifying the mutual information between the representation and the halos&#39; ground-truth density profiles. A two-dimensional representation is sufficient to accurately model the density profiles up to the virial radius; however, a three-dimensional representation is required to describe the outer profiles beyond the virial radius. The additional dimension in the representation contains information about the infalling material in the outer profiles of dark matter halos, thus discovering the splashback boundary of halos without prior knowledge of the halos&#39; dynamical history.

preprint2022arXiv

EDGE: the puzzling ellipticity of Eridanus II&#39;s star cluster and its implications for dark matter at the heart of an ultra-faint dwarf

The Eridanus II (EriII) &#39;ultra-faint&#39; dwarf has a large ($15\,\text{pc}$) and low mass ($4.3\times10^3\,\text{M}_\odot$) star cluster (SC) offset from its centre by $23\pm3\,\text{pc}$ in projection. Its size and offset are naturally explained if EriII has a central dark matter core, but such a core may be challenging to explain in a $Λ$CDM cosmology. In this paper, we revisit the survival and evolution of EriII&#39;s SC, focussing for the first time on its puzzlingly large ellipticity ($0.31^{+0.05}_{-0.06}$). We perform a suite of 960 direct $N$-body simulations of SCs, orbiting within a range of spherical background potentials fit to ultra-faint dwarf (UFD) galaxy simulations. We find only two scenarios that come close to explaining EriII&#39;s SC. In the first, EriII has a low density dark matter core (of size $\sim70\,\text{pc}$ and density $\lesssim2\times10^8\,\text{M}_{\odot}\,\text{kpc}^{-3}$). In this model, the high ellipticity of EriII&#39;s SC is set at birth, with the lack of tidal forces in the core allowing its ellipticity to remain frozen in for long times. In the second, EriII&#39;s SC orbits in a partial core, with its high ellipticity owing to its imminent tidal destruction. However, this latter model struggles to reproduce the large size of EriII&#39;s SC, and it predicts substantial tidal tails around EriII&#39;s SC that should have already been seen in the data. This leads us to favour the cored model. We discuss potential caveats to these findings, and the implications of the cored model for galaxy formation and the nature of dark matter.

preprint2022arXiv

EDGE: The sensitivity of ultra-faint dwarfs to a metallicity-dependent initial mass function

Motivated by the observed bottom-light initial mass function (IMF) in faint dwarfs, we study how a metallicity-dependent IMF affects the feedback budget and observables of an ultra-faint dwarf galaxy. We model the evolution of a low-mass ($\approx 8 \, \times \, 10^{8} \, \rm M_{\odot}$) dark matter halo with cosmological, zoomed hydrodynamical simulations capable of resolving individual supernovae explosions, which we complement with an empirically motivated subgrid prescription for systematic IMF variations. In this framework, at the low gas metallicities typical of faint dwarfs, the IMF of newborn stellar populations becomes top-heavy, increasing the efficiency of supernova and photoionization feedback in regulating star formation. This results in a 100-fold reduction of the final stellar mass of the dwarf compared to a canonical IMF, at fixed dynamical mass. The increase in the feedback budget is nonetheless met by increased metal production from more numerous massive stars, leading to nearly constant iron content at $z=0$. A metallicity-dependent IMF therefore provides a mechanism to produce low-mass ($\rm M_{\star}\sim 10^3 \rm M_{\odot}$), yet enriched ($\rm [Fe/H]\approx -2$) field dwarf galaxies, thus opening a self-consistent avenue to populate the plateau in $\rm [Fe/H]$ at the faintest end of the mass--metallicity relation.

preprint2022arXiv

EDGE: What shapes the relationship between HI and stellar observables in faint dwarf galaxies?

We show how the interplay between feedback and mass-growth histories introduces scatter in the relationship between stellar and neutral gas properties of field faint dwarf galaxies ($M_{\star} \lessapprox 10^{6} M_{\odot}$). Across a suite of cosmological, high-resolution zoomed simulations, we find that dwarf galaxies of stellar masses $10^5 \leq M_{\star} \leq 10^{6} M_{\odot}$ are bimodal in their cold gas content, being either HI-rich or HI-deficient. This bimodality is generated through the coupling between (i) the modulation of HI contents by the background of ultraviolet radiation (UVB) at late times and (ii) the significant scatter in the stellar-mass-halo-mass relationship induced by reionization. Furthermore, our HI-rich dwarfs exhibit disturbed and time-variable neutral gas distributions primarily due to stellar feedback. Over the last four billion years, we observe order-of-magnitude changes around the median $M_{HI}$, factor-of-a-few variations in HI spatial extents, and spatial offsets between HI and stellar components regularly exceeding the galaxies&#39; optical sizes. Time variability introduces further scatter in the $M_{\star}-M_{HI}$ relation and affects a galaxy&#39;s detectability in HI at any given time. These effects will need to be accounted for when interpreting observations of the population of faint, HI-bearing dwarfs by the combination of optical and radio wide, deep surveys.

preprint2022arXiv

Galaxy mergers can initiate quenching by unlocking an AGN-driven transformation of the baryon cycle

We use zoom simulations to show how merger-driven disruption of the gas disc in a galaxy provides its central active galactic nucleus (AGN) with fuel to drive outflows that entrain and expel a significant fraction of the circumgalactic medium (CGM). This in turn suppresses replenishment of the interstellar medium, causing the galaxy to quench up to several Gyr after the merger. We start by performing a zoom simulation of a present-day star-forming disc galaxy with the EAGLE galaxy formation model. Then, we re-simulate the galaxy with controlled changes to its initial conditions, using the genetic modification technique. These modifications either increase or decrease the stellar mass ratio of the galaxy&#39;s last significant merger, which occurs at $z\approx 0.74$. The halo reaches the same present-day mass in all cases, but changing the mass ratio of the merger yields markedly different galaxy and CGM properties. We find that a merger can unlock rapid growth of the central supermassive black hole if it disrupts the co-rotational motion of gas in the black hole&#39;s vicinity. Conversely, if a less disruptive merger occurs and gas close to the black hole is not disturbed, the AGN does not strongly affect the CGM, and consequently the galaxy continues to form stars. Our result illustrates how a unified view of AGN feedback, the baryon cycle and the interstellar medium is required to understand how mergers and quenching are connected over long timescales.

preprint2021arXiv

Angular momentum evolution can be predicted from cosmological initial conditions

The angular momentum of dark matter haloes controls their spin magnitude and orientation, which in turn influences the galaxies therein. However, the process by which dark matter haloes acquire angular momentum is not fully understood; in particular, it is unclear whether angular momentum growth is stochastic. To address this question, we extend the genetic modification technique to allow control over the angular momentum of any region in the initial conditions. Using this technique to produce a sequence of modified simulations, we can then investigate whether changes to the angular momentum of a specified region in the evolved universe can be accurately predicted from changes in the initial conditions alone. We find that the angular momentum in regions with modified initial conditions can be predicted between 2 and 4 times more accurately than expected from applying tidal torque theory. This result is masked when analysing the angular momentum of haloes, because particles in the outskirts of haloes dominate the angular momentum budget. We conclude that the angular momentum of Lagrangian patches is highly predictable from the initial conditions, with apparent chaotic behaviour being driven by stochastic changes to the arbitrary boundary defining the halo.

preprint2021arXiv

Quenching and morphological evolution due to circumgalactic gas expulsion in a simulated galaxy with a controlled assembly history

We examine the influence of dark matter halo assembly on the evolution of a simulated $\sim L^\star$ galaxy. Starting from a zoom-in simulation of a star-forming galaxy evolved with the EAGLE galaxy formation model, we use the genetic modification technique to create a pair of complementary assembly histories: one in which the halo assembles later than in the unmodified case, and one in which it assembles earlier. Delayed assembly leads to the galaxy exhibiting a greater present-day star formation rate than its unmodified counterpart, whilst in the accelerated case the galaxy quenches at $z\simeq 1$, and becomes spheroidal. We simulate each assembly history nine times, adopting different seeds for the random number generator used by EAGLE&#39;s stochastic subgrid implementations of star formation and feedback. The systematic changes driven by differences in assembly history are significantly stronger than the random scatter induced by this stochasticity. The sensitivity of $\sim L^\star$ galaxy evolution to dark matter halo assembly follows from the close coupling of the growth histories of the central black hole (BH) and the halo, such that earlier assembly fosters the formation of a more massive BH, and more efficient expulsion of circumgalactic gas. In response to this expulsion, the circumgalactic medium reconfigures at a lower density, extending its cooling time and thus inhibiting the replenishment of the interstellar medium. Our results indicate that halo assembly history significantly influences the evolution of $\sim L^\star$ central galaxies, and that the expulsion of circumgalactic gas is a crucial step in quenching them.

preprint2020arXiv

EDGE: From quiescent to gas-rich to star-forming low-mass dwarf galaxies

We study how star formation is regulated in low-mass field dwarf galaxies ($10^5 \leq M_{\star} \leq 10^6 \, \text{M}_{\odot}$), using cosmological high-resolution ($3 \, \text{pc}$) hydrodynamical simulations. Cosmic reionization quenches star formation in all our simulated dwarfs, but three galaxies with final dynamical masses of $3 \times 10^{9} \,\text{M}_{\odot}$ are subsequently able to replenish their interstellar medium by slowly accreting gas. Two of these galaxies re-ignite and sustain star formation until the present day at an average rate of $10^{-5} \, \text{M}_{\odot} \, \text{yr}^{-1}$, highly reminiscent of observed low-mass star-forming dwarf irregulars such as Leo T. The resumption of star formation is delayed by several billion years due to residual feedback from stellar winds and Type Ia supernovae; even at $z=0$, the third galaxy remains in a temporary equilibrium with a large gas content but without any ongoing star formation. Using the &#34;genetic modification&#39;&#39; approach, we create an alternative mass growth history for this gas-rich quiescent dwarf and show how a small $(0.2\,\mathrm{dex})$ increase in dynamical mass can overcome residual stellar feedback, re-igniting star formation. The interaction between feedback and mass build-up produces a diversity in the stellar ages and gas content of low-mass dwarfs, which will be probed by combining next-generation HI and imaging surveys.

preprint2020arXiv

GenetIC -- a new initial conditions generator to support genetically modified zoom simulations

We present genetIC, a new code for generating initial conditions for cosmological N-body simulations. The code allows precise, user-specified alterations to be made to arbitrary regions of the simulation (while maintaining consistency with the statistical ensemble). These &#34;genetic modifications&#34; allow, for example, the history, mass, or environment of a target halo to be altered in order to study the effect on their evolution. The code natively supports initial conditions with nested zoom regions at progressively increasing resolution. Modifications in the high-resolution region must propagate self-consistently onto the lower resolution grids; to enable this while maintaining a small memory footprint, we introduce a Fourier-space filtering approach to generating fields at variable resolution. Due to a close correspondence with modifications, constrained initial conditions can also be produced by genetIC (for example with the aim of matching structures in the local Universe). We test the accuracy of modifications performed within zoom initial conditions. The code achieves sub-percent precision, which is easily sufficient for current applications in galaxy formation.

preprint2020arXiv

Massive neutrinos and degeneracies in Lyman-alpha forest simulations

Using a suite of hydrodynamical simulations with cold dark matter, baryons, and neutrinos, we present a detailed study of the effect of massive neutrinos on the 1-D and 3-D flux power spectra of the Lyman-$α$ (Ly$α$) forest. The presence of massive neutrinos in cosmology induces a scale- and time-dependent suppression of structure formation that is strongest on small scales. Measuring this suppression is a key method for inferring neutrino masses from cosmological data, and is one of the main goals of ongoing and future surveys like eBOSS, DES, LSST, Euclid or DESI. The clustering in the Ly$α$ forest traces the quasi-linear power at late times and on small scales. In combination with observations of the cosmic microwave background, the forest therefore provides some of the tightest constraints on the sum of the neutrino masses. However there is a well-known degeneracy between $Σm_ν$ and the amplitude of perturbations in the linear matter power spectrum. We study the corresponding degeneracy in the 1-D flux power spectrum of the Ly$α$ forest, and for the first time also study this degeneracy in the 3-D flux power spectrum. We show that the non-linear effects of massive neutrinos on the Ly$α$ forest, beyond the effect of linear power amplitude suppression, are negligible, and this degeneracy persists in the Ly$α$ forest observables to a high precision. We discuss the implications of this degeneracy for choosing parametrisations of the Ly$α$ forest for cosmological analysis.

preprint2019arXiv

Nonlinear Dynamics of the Cold Atom Analog False Vacuum

We investigate the nonlinear dynamics of cold atom systems that can in principle serve as quantum simulators of false vacuum decay. The analog false vacuum manifests as a metastable vacuum state for the relative phase in a two-species Bose-Einstein condensate (BEC), induced by a driven periodic coupling between the two species. In the appropriate low energy limit, the evolution of the relative phase is approximately governed by a relativistic wave equation exhibiting true and false vacuum configurations. In previous work, a linear stability analysis identified exponentially growing short-wavelength modes driven by the time-dependent coupling. These modes threaten to destabilize the analog false vacuum. Here, we employ numerical simulations of the coupled Gross-Pitaevski equations (GPEs) to determine the non-linear evolution of these linearly unstable modes. We find that unless a physical mechanism modifies the GPE on short length scales, the analog false vacuum is indeed destabilized. We briefly discuss various physically expected corrections to the GPEs that may act to remove the exponentially unstable modes. To investigate the resulting dynamics in cases where such a removal mechanism exists, we implement a hard UV cutoff that excludes the unstable modes as a simple model for these corrections. We use this to study the range of phenomena arising from such a system. In particular, we show that by modulating the strength of the time-dependent coupling, it is possible to observe the crossover between a second and first order phase transition out of the false vacuum.

preprint2017arXiv

Large-scale fluctuations in the cosmic ionising background: the impact of beamed source emission

When modelling the ionisation of gas in the intergalactic medium after reionisation, it is standard practice to assume a uniform radiation background. This assumption is not always appropriate; models with radiative transfer show that large-scale ionisation rate fluctuations can have an observable impact on statistics of the Lyman-alpha forest. We extend such calculations to include beaming of sources, which has previously been neglected but which is expected to be important if quasars dominate the ionising photon budget. Beaming has two effects: first, the physical number density of ionising sources is enhanced relative to that directly observed; and second, the radiative transfer itself is altered. We calculate both effects in a hard-edged beaming model where each source has a random orientation, using an equilibrium Boltzmann hierarchy in terms of spherical harmonics. By studying the statistical properties of the resulting ionisation rate and HI density fields at redshift $z\sim 2.3$, we find that the two effects partially cancel each other; combined, they constitute a maximum $5\%$ correction to the power spectrum $P_{\mathrm{HI}}(k)$ at $k=0.04 \, h/\mathrm{Mpc}$. On very large scales ($k<0.01\, h/\mathrm{Mpc}$) the source density renormalisation dominates; it can reduce, by an order of magnitude, the contribution of ionising shot-noise to the intergalactic HI power spectrum. The effects of beaming should be considered when interpreting future observational datasets.

preprint2015arXiv

Particle tagging and its implications for stellar population dynamics

We establish a controlled comparison between the properties of galactic stellar halos obtained with hydrodynamical simulations and with `particle tagging&#39;. Tagging is a fast way to obtain stellar population dynamics: instead of tracking gas and star formation, it `paints&#39; stars directly onto a suitably defined subset of dark matter particles in a collisionless, dark-matter-only simulation.Our study shows that there are conditions under which particle tagging generates good fits to the hydrodynamical stellar density profiles of a central Milky-Way-like galaxy and its most prominent substructure. Phase-space diffusion processes are crucial to reshaping the distribution of stars in infalling spheroidal systems and hence the final stellar halo. We conclude that the success of any particular tagging scheme hinges on this diffusion being taken into account, at a minimum by making use of `live&#39; tagging schemes, in which particles are regularly tagged throughout the evolution of a galaxy.