BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Andrew Pontzen

Andrew Pontzen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.COastro-ph.GAastro-ph.IMgr-qccond-mat.quant-gashep-thArtificial Intelligenceastro-ph.SRhep-phMachine Learning

Catalog footprint

What is connected

42works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

42 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.

0 citations0 reviewsNo code linkedOpen access
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', 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.

0 citations0 reviewsNo code linkedOpen access
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' 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' dynamical history.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

The Eridanus II (EriII) 'ultra-faint' 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'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'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'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'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's SC, and it predicts substantial tidal tails around EriII'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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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' optical sizes. Time variability introduces further scatter in the $M_{\star}-M_{HI}$ relation and affects a galaxy'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.

0 citations0 reviewsNo code linkedOpen access
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'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'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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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'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.

0 citations0 reviewsNo code linkedOpen access
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 "genetic modification'' 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.

0 citations0 reviewsNo code linkedOpen access
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 "genetic modifications" 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A framework for testing isotropy with the cosmic microwave background

We present a new framework for testing the isotropy of the Universe using cosmic microwave background data, building on the nested-sampling ANICOSMO code. Uniquely, we are able to constrain the scalar, vector and tensor degrees of freedom alike; previous studies only considered the vector mode (linked to vorticity). We employ Bianchi type VII$_h$ cosmologies to model the anisotropic Universe, from which other types may be obtained by taking suitable limits. In a separate development, we improve the statistical analysis by including the effect of Bianchi power in the high-$\ell$, as well as the low-$\ell$, likelihood. To understand the effect of all these changes, we apply our new techniques to WMAP data. We find no evidence for anisotropy, constraining shear in the vector mode to $(σ_V/H)_0 < 1.7 \times 10^{-10}$ (95% CL). For the first time, we place limits on the tensor mode; unlike other modes, the tensor shear can grow from a near-isotropic early Universe. The limit on this type of shear is $(σ_{T,\rm reg}/H)_0 < 2.4 \times 10^{-7}$ (95% CL).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Bursty star formation feedback and cooling outflows

We study how outflows of gas launched from a central galaxy undergoing repeated starbursts propagate through the circumgalactic medium (CGM), using the simulation code RAMSES. We assume that the outflow from the disk can be modelled as a rapidly moving bubble of hot gas at $\mathrm{\sim1\;kpc}$ above disk, then ask what happens as it moves out further into the halo around the galaxy on $\mathrm{\sim 100\;kpc}$ scales. To do this we run 60 two-dimensional simulations scanning over parameters of the outflow. Each of these is repeated with and without radiative cooling, assuming a primordial gas composition to give a lower bound on the importance of cooling. In a large fraction of radiative-cooling cases we are able to form rapidly outflowing cool gas from in situ cooling of the flow. We show that the amount of cool gas formed depends strongly on the 'burstiness' of energy injection; sharper, stronger bursts typically lead to a larger fraction of cool gas forming in the outflow. The abundance ratio of ions in the CGM may therefore change in response to the detailed historical pattern of star formation. For instance, outflows generated by star formation with short, intense bursts contain up to 60 per cent of their gas mass at temperatures $<5 \times 10^4\,\mathrm{K}$; for near-continuous star formation the figure is $\lesssim$ 5 per cent. Further study of cosmological simulations, and of idealised simulations with e.g., metal-cooling, magnetic fields and/or thermal conduction, will help to understand the precise signature of bursty outflows on observed ion abundances.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Cosmological $N$-body simulations with suppressed variance

We present and test a method that dramatically reduces variance arising from the sparse sampling of wavemodes in cosmological simulations. The method uses two simulations which are fixed (the initial Fourier mode amplitudes are fixed to the ensemble average power spectrum) and paired (with initial modes exactly out of phase). We measure the power spectrum, monopole and quadrupole redshift-space correlation functions, halo mass function and reduced bispectrum at $z=1$. By these measures, predictions from a fixed pair can be as precise on non-linear scales as an average over 50 traditional simulations. The fixing procedure introduces a non-Gaussian correction to the initial conditions; we give an analytic argument showing why the simulations are still able to predict the mean properties of the Gaussian ensemble. We anticipate that the method will drive down the computational time requirements for accurate large-scale explorations of galaxy bias and clustering statistics, enabling more precise comparisons with theoretical models, and facilitating the use of numerical simulations in cosmological data interpretation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

How isotropic is the Universe?

A fundamental assumption in the standard model of cosmology is that the Universe is isotropic on large scales. Breaking this assumption leads to a set of solutions to Einstein's field equations, known as Bianchi cosmologies, only a subset of which have ever been tested against data. For the first time, we consider all degrees of freedom in these solutions to conduct a general test of isotropy using cosmic microwave background temperature and polarization data from Planck. For the vector mode (associated with vorticity), we obtain a limit on the anisotropic expansion of $(σ_V/H)_0 < 4.7 \times 10^{-11}$ (95% CI), which is an order of magnitude tighter than previous Planck results that used CMB temperature only. We also place upper limits on other modes of anisotropic expansion, with the weakest limit arising from the regular tensor mode, $(σ_{T,\rm reg}/H)_0<1.0 \times 10^{-6}$ (95% CI). Including all degrees of freedom simultaneously for the first time, anisotropic expansion of the Universe is strongly disfavoured, with odds of 121,000:1 against.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

How to quench a galaxy

We show how the interplay between active galactic nuclei (AGN) and merger history determines whether a galaxy quenches star formation at high redshift. We first simulate, in a full cosmological context, a galaxy of total dynamical mass $10^{12}\,M_{\odot}$ at $z=2$. Then we systematically alter the accretion history of the galaxy by minimally changing the linear overdensity in the initial conditions. This "genetic modification" approach allows the generation of three sets of $Λ$CDM initial conditions leading to maximum merger ratios of 1:10, 1:5 and 2:3 respectively. The changes leave the final halo mass, large scale structure and local environment unchanged, providing a controlled numerical experiment. Interaction between the AGN physics and mergers in the three cases lead respectively to a star-forming, temporarily-quenched and permanently-quenched galaxy. However the differences do not primarily lie in the black hole accretion rates, but in the kinetic effects of the merger: the galaxy is resilient against AGN feedback unless its gaseous disk is first disrupted. Typical accretion rates are comparable in the three cases, falling below $0.1\,M_{\odot}$ yr$^{-1}$, equivalent to around $2\%$ of the Eddington rate or $10^{-3}$ times the pre-quenching star formation rate, in agreement with observations. This low level of black hole accretion can be sustained even when there is insufficient dense cold gas for star formation. Conversely, supernova feedback is too distributed to generate outflows in high-mass systems, and cannot maintain quenching over periods longer than the halo gas cooling time.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

SILC: a new Planck Internal Linear Combination CMB temperature map using directional wavelets

We present new clean maps of the CMB temperature anisotropies (as measured by Planck) constructed with a novel internal linear combination (ILC) algorithm using directional, scale-discretised wavelets --- Scale-discretised, directional wavelet ILC or SILC. Directional wavelets, when convolved with signals on the sphere, can separate the anisotropic filamentary structures which are characteristic of both the CMB and foregrounds. Extending previous component separation methods, which use the frequency, spatial and harmonic signatures of foregrounds to separate them from the cosmological background signal, SILC can additionally use morphological information in the foregrounds and CMB to better localise the cleaning algorithm. We test the method on Planck data and simulations, demonstrating consistency with existing component separation algorithms, and discuss how to optimise the use of morphological information by varying the number of directional wavelets as a function of spatial scale. We find that combining the use of directional and axisymmetric wavelets depending on scale could yield higher quality CMB temperature maps. Our results set the stage for the application of SILC to polarisation anisotropies through an extension to spin wavelets.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Spin-SILC: CMB polarisation component separation with spin wavelets

We present Spin-SILC, a new foreground component separation method that accurately extracts the cosmic microwave background (CMB) polarisation $E$ and $B$ modes from raw multifrequency Stokes $Q$ and $U$ measurements of the microwave sky. Spin-SILC is an internal linear combination method that uses spin wavelets to analyse the spin-2 polarisation signal $P = Q + iU$. The wavelets are additionally directional (non-axisymmetric). This allows different morphologies of signals to be separated and therefore the cleaning algorithm is localised using an additional domain of information. The advantage of spin wavelets over standard scalar wavelets is to simultaneously and self-consistently probe scales and directions in the polarisation signal $P = Q + iU$ and in the underlying $E$ and $B$ modes, therefore providing the ability to perform component separation and $E$-$B$ decomposition concurrently for the first time. We test Spin-SILC on full-mission Planck simulations and data and show the capacity to correctly recover the underlying cosmological $E$ and $B$ modes. We also demonstrate a strong consistency of our CMB maps with those derived from existing component separation methods. Spin-SILC can be combined with the pseudo- and pure $E$-$B$ spin wavelet estimators presented in a companion paper to reliably extract the cosmological signal in the presence of complicated sky cuts and noise. Therefore, it will provide a computationally-efficient method to accurately extract the CMB $E$ and $B$ modes for future polarisation experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Genetically modified halos: towards controlled experiments in $Λ$CDM galaxy formation

We propose a method to generate `genetically-modified' (GM) initial conditions for high-resolution simulations of galaxy formation in a cosmological context. Building on the Hoffman-Ribak algorithm, we start from a reference simulation with fully random initial conditions, then make controlled changes to specific properties of a single halo (such as its mass and merger history). The algorithm demonstrably makes minimal changes to other properties of the halo and its environment, allowing us to isolate the impact of a given modification. As a significant improvement over previous work, we are able to calculate the abundance of the resulting objects relative to the reference simulation. Our approach can be applied to a wide range of cosmic structures and epochs; here we study two problems as a proof-of-concept. First, we investigate the change in density profile and concentration as the collapse time of three individual halos are varied at fixed final mass, showing good agreement with previous statistical studies using large simulation suites. Second, we modify the $z=0$ mass of halos to show that our theoretical abundance calculations correctly recover the halo mass function. The results demonstrate that the technique is robust, opening the way to controlled experiments in galaxy formation using hydrodynamic zoom simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Inverted initial conditions: exploring the growth of cosmic structure and voids

We introduce and explore "paired" cosmological simulations. A pair consists of an A and B simulation with initial conditions related by the inversion $δ_A(x, t_{initial})=-δ_B(x,t_{initial})$ (underdensities substituted for overdensities and vice versa). We argue that the technique is valuable for improving our understanding of cosmic structure formation. The A and B fields are by definition equally likely draws from ΛCDM initial conditions, and in the linear regime evolve identically up to the overall sign. As non-linear evolution takes hold, a region that collapses to form a halo in simulation A will tend to expand to create a void in simulation B. Applications include (i) contrasting the growth of A-halos and B-voids to test excursion-set theories of structure formation; (ii) cross-correlating the density field of the A and B universes as a novel test for perturbation theory; and (iii) canceling error terms by averaging power spectra between the two boxes. Generalizations of the method to more elaborate field transformations are suggested.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Milking the spherical cow: on aspherical dynamics in spherical coordinates

Galaxies and the dark matter halos that host them are not spherically symmetric, yet spherical symmetry is a helpful simplifying approximation for idealised calculations and analysis of observational data. The assumption leads to an exact conservation of angular momentum for every particle, making the dynamics unrealistic. But how much does that inaccuracy matter in practice for analyses of stellar distribution functions, collisionless relaxation, or dark matter core-creation? We provide a general answer to this question for a wide class of aspherical systems; specifically, we consider distribution functions that are "maximally stable", i.e. that do not evolve at first order when external potentials (which arise from baryons, large scale tidal fields or infalling substructure) are applied. We show that a spherically-symmetric analysis of such systems gives rise to the false conclusion that the density of particles in phase space is ergodic (a function of energy alone). Using this idea we are able to demonstrate that: (a) observational analyses that falsely assume spherical symmetry are made more accurate by imposing a strong prior preference for near-isotropic velocity dispersions in the centre of spheroids; (b) numerical simulations that use an idealised spherically-symmetric setup can yield misleading results and should be avoided where possible; and (c) triaxial dark matter halos (formed in collisionless cosmological simulations) nearly attain our maximally-stable limit, but their evolution freezes out before reaching it.

0 citations0 reviewsNo code linkedOpen access
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'. Tagging is a fast way to obtain stellar population dynamics: instead of tracking gas and star formation, it `paints' 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' tagging schemes, in which particles are regularly tagged throughout the evolution of a galaxy.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Cold dark matter heats up

One of the principal discoveries in modern cosmology is that standard model particles (including baryons, leptons and photons) together comprise only 5% of the mass-energy budget of the Universe. The remaining 95% consists of dark energy and dark matter (DM). Consequently our picture of the universe is known as ΛCDM, with Λ denoting dark energy and CDM cold dark matter. ΛCDM is being challenged by its apparent inability to explain the low density of DM measured at the centre of cosmological systems, ranging from faint dwarf galaxies to massive clusters containing tens of galaxies the size of the Milky Way. But before making conclusions one should carefully include the effect of gas and stars, which were historically seen as merely a passive component during the assembly of galaxies. We now understand that these can in fact significantly alter the DM component, through a coupling based on rapid gravitational potential fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Constraints on ionising photon production from the large-scale Lyman-alpha forest

Recent work has shown that the z~2.5 Lyman-alpha forest on large scales encodes information about the galaxy and quasar populations that keep the intergalactic medium photoionized. We present the first forecasts for constraining the populations with data from current and next-generation surveys. At a minimum the forest should tell us whether galaxies or, conversely, quasars dominate the photon production. The number density and clustering strength of the ionising sources might be estimated to sub-10% precision with a DESI-like survey if degeneracies (e.g., with the photon mean-free-path, small-scale clustering power normalization and potentially other astrophysical effects) can be broken by prior information. We demonstrate that, when inhomogeneous ionisation is correctly handled, constraints on dark energy do not degrade.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Faint dwarfs as a test of DM models: WDM vs. CDM

We use high resolution Hydro$+$N-Body cosmological simulations to compare the assembly and evolution of a small field dwarf (stellar mass ~ 10$^{6-7}$ M$\odot$, total mass 10$^{10}$ M$\odot$ in $Λ$ dominated CDM and 2keV WDM cosmologies. We find that star formation (SF) in the WDM model is reduced and delayed by 1-2 Gyr relative to the CDM model, independently of the details of SF and feedback. Independent of the DM model, but proportionally to the SF efficiency, gas outflows lower the central mass density through `dynamical heating', such that all realizations have circular velocities $<$ 20kms at 500$~$pc, in agreement with local kinematic constraints. As a result of dynamical heating, older stars are less centrally concentrated than younger stars, similar to stellar population gradients observed in nearby dwarf galaxies. Introducing an important diagnostic of SF and feedback models, we translate our simulations into artificial color-magnitude diagrams and star formation histories in order to directly compare to available observations. The simulated galaxies formed most of their stars in many $\sim$10 Myr long bursts. The CDM galaxy has a global SFH, HI abundance and Fe/H and alpha-elements distribution well matched to current observations of dwarf galaxies. These results highlight the importance of directly including `baryon physics' in simulations when 1) comparing predictions of galaxy formation models with the kinematics and number density of local dwarf galaxies and 2) differentiating between CDM and non-standard models with different DM or power spectra.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Quintessence in a quandary: prior dependence in dark energy models

The archetypal theory of dark energy is quintessence: a minimally coupled scalar field with a canonical kinetic energy and potential. By studying random potentials we show that quintessence imposes a restricted set of priors on the equation of state of dark energy. Focusing on the commonly-used parametrisation, $w(a)\approx w_0+w_a(1-a)$, we show that there is a natural scale and direction in the $(w_0, w_a)$ plane that distinguishes quintessence as a general framework. We calculate the expected information gain for a given survey and show that, because of the non-trivial prior information, it is a function of more than just the figure of merit. This allows us to make a quantitative case for novel survey strategies. We show that the scale of the prior sets target observational requirements for gaining significant information. This corresponds to a figure of merit FOM$\gtrsim 200$, a requirement that future galaxy redshift surveys will meet.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Scale-dependent bias in the BAO-scale intergalactic neutral hydrogen

I discuss fluctuations in the neutral hydrogen (HI) density of the z~2.3 intergalactic medium and show that their relation to cosmic overdensity is strongly scale-dependent. This behaviour arises from a linearized version of the well-known "proximity effect", in which bright sources suppress atomic hydrogen density. Using a novel, systematic and detailed linear-theory radiative transfer calculation, I demonstrate how HI density consequently anti-correlates with total matter density when averaged on scales exceeding the Lyman-limit mean-free-path. The radiative transfer thumbprint is highly distinctive and should be measurable in the Lyman-alpha forest. Effects extend to sufficiently small scales to generate significant distortion of the correlation function shape around the baryon acoustic oscillation peak, although the peak location shifts only by 1.2 percent for a mean source bias of b_j=3. The distortion changes significantly with b_j and other astrophysical parameters; measuring it should provide a helpful observational constraint on the nature of ionizing photon sources in the near future.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Astropy: A Community Python Package for Astronomy

We present the first public version (v0.2) of the open-source and community-developed Python package, Astropy. This package provides core astronomy-related functionality to the community, including support for domain-specific file formats such as Flexible Image Transport System (FITS) files, Virtual Observatory (VO) tables, and common ASCII table formats, unit and physical quantity conversions, physical constants specific to astronomy, celestial coordinate and time transformations, world coordinate system (WCS) support, generalized containers for representing gridded as well as tabular data, and a framework for cosmological transformations and conversions. Significant functionality is under active development, such as a model fitting framework, VO client and server tools, and aperture and point spread function (PSF) photometry tools. The core development team is actively making additions and enhancements to the current code base, and we encourage anyone interested to participate in the development of future Astropy versions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Estimating the large-scale angular power spectrum in the presence of systematics: a case study of Sloan Digital Sky Survey quasars

The angular power spectrum is a powerful statistic for analysing cosmological signals imprinted in the clustering of matter. However, current galaxy and quasar surveys cover limited portions of the sky, and are contaminated by systematics that can mimic cosmological signatures and jeopardise the interpretation of the measured power spectra. We provide a framework for obtaining unbiased estimates of the angular power spectra of large-scale structure surveys at the largest scales using quadratic estimators. The method is tested by analysing the 600 CMASS mock catalogues constructed by Manera et al. (2013) for the Baryon Oscillation Spectroscopic Survey (BOSS). We then consider the Richards et al. (2009) catalogue of photometric quasars from the Sixth Data Release (DR6) of the Sloan Digital Sky Survey (SDSS), which is known to include significant stellar contamination and systematic uncertainties. Focusing on the sample of ultraviolet-excess (UVX) sources, we show that the excess clustering power present on the largest-scales can be largely mitigated by making use of improved sky masks and projecting out the modes corresponding to the principal systematics. In particular, we find that the sample of objects with photometric redshift $1.3 < z_p < 2.2$ exhibits no evidence of contamination when using our most conservative mask and mode projection. This indicates that any residual systematics are well within the statistical uncertainties. We conclude that, using our approach, this sample can be used for cosmological studies.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Baryons Matter: Why Luminous Satellite Galaxies Have Reduced Central Masses

Using high resolution cosmological hydrodynamical simulations of Milky Way-massed disk galaxies, we demonstrate that supernovae feedback and tidal stripping lower the central masses of bright (-15 < M_V < -8) satellite galaxies. These simulations resolve high density regions, comparable to giant molecular clouds, where stars form. This resolution allows us to adopt a prescription for H_2 formation and destruction that ties star formation to the presence of shielded, molecular gas. Before infall, supernova feedback from the clumpy, bursty star formation captured by this physically motivated model leads to reduced dark matter (DM) densities and shallower inner density profiles in the massive satellite progenitors (Mvir > 10^9 Msun, Mstar > 10^7 Msun) compared to DM-only simulations. The progenitors of the lower mass satellites are unable to maintain bursty star formation histories, due to both heating at reionization and gas loss from initial star forming events, preserving the steep inner density profile predicted by DM-only simulations. After infall, tidal stripping acts to further reduce the central densities of the luminous satellites, particularly those that enter with cored dark matter halos, increasing the discrepancy in the central masses predicted by baryon+DM and DM-only simulations. We show that DM-only simulations, which neglect the baryonic effects described in this work, produce denser satellites with larger central velocities. We provide a simple correction to the central DM mass predicted for satellites by DM-only simulations. We conclude that DM-only simulations should be used with great caution when interpreting kinematic observations of the Milky Way's dwarf satellites.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Conserved actions, maximum entropy and dark matter halos

We use maximum entropy arguments to derive the phase space distribution of a virialized dark matter halo. Our distribution function gives an improved representation of the end product of violent relaxation. This is achieved by incorporating physically motivated dynamical constraints (specifically on orbital actions) which prevent arbitrary redistribution of energy. We compare the predictions with three high-resolution dark matter simulations of widely varying mass. The numerical distribution function is accurately predicted by our argument, producing an excellent match for the vast majority of particles. The remaining particles constitute the central cusp of the halo (<4% of the dark matter). They can be accounted for within the presented framework once the short dynamical timescales of the centre are taken into account.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Cusp-core transformations in dwarf galaxies: observational predictions

The presence of a dark matter core in the central kiloparsec of many dwarf galaxies has been a long standing problem in galaxy formation theories based on the standard cold dark matter paradigm. Recent cosmological simulations, based on Smooth Particle Hydrodynamics and rather strong feedback recipes have shown that it was indeed possible to form extended dark matter cores using baryonic processes related to a more realistic treatment of the interstellar medium. Using adaptive mesh refinement, together with a new, stronger supernovae feedback scheme that we have recently implemented in the RAMSES code, we show that it is also possible to form a prominent dark matter core within the well-controlled framework of an isolated, initially cuspy, 10 billion solar masses dark matter halo. Although our numerical experiment is idealized, it allows a clean and unambiguous identification of the dark matter core formation process. Our dark matter inner profile is well fitted by a pseudo-isothermal profile with a core radius of 800 pc. The core formation mechanism is consistent with the one proposed recently by Pontzen & Governato. We highlight two key observational predictions of all simulations that find cusp-core transformations: (i) a bursty star formation history with peak to trough ratio of 5 to 10 and a duty cycle comparable to the local dynamical time; and (ii) a stellar distribution that is hot with v/sigma=1. We compare the observational properties of our model galaxy with recent measurements of the isolated dwarf WLM. We show that the spatial and kinematical distribution of stars and HI gas are in striking agreement with observations, supporting the fundamental role played by stellar feedback in shaping both the stellar and dark matter distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

How supernova feedback turns dark matter cusps into cores

We propose and successfully test against new cosmological simulations a novel analytical description of the physical processes associated with the origin of cored dark matter density profiles. In the simulations, the potential in the central kiloparsec changes on sub-dynamical timescales over the redshift interval 4 > z > 2 as repeated, energetic feedback generates large underdense bubbles of expanding gas from centrally-concentrated bursts of star formation. The model demonstrates how fluctuations in the central potential irreversibly transfer energy into collisionless particles, thus generating a dark matter core. A supply of gas undergoing collapse and rapid expansion is therefore the essential ingredient. The framework, based on a novel impulsive approximation, breaks with the reliance on adiabatic approximations which are inappropriate in the rapidly-changing limit. It shows that both outflows and galactic fountains can give rise to cusp-flattening, even when only a few per cent of the baryons form stars. Dwarf galaxies maintain their core to the present time. The model suggests that constant density dark matter cores will be generated in systems of a wide mass range if central starbursts or AGN phases are sufficiently frequent and energetic.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

MAGICC haloes: confronting simulations with observations of the circumgalactic medium at z=0

We explore the circumgalactic medium (CGM) of two simulated star-forming galaxies with luminosities L ~ 0.1 and 1 L* generated using the smooth particle hydrodynamic code GASOLINE. These simulations are part of the Making Galaxies In a Cosmological Context (MAGICC) program in which the stellar feedback is tuned to match the stellar mass-halo mass relationship. For comparison, each galaxy was also simulated using a 'lower feedback' (LF) model which has strength comparable to other implementations in the literature. The 'MAGICC feedback' (MF) model has a higher incidence of massive stars and an approximately two times higher energy input per supernova. Apart from the low-mass halo using LF, each galaxy exhibits a metal-enriched CGM that extends to approximately the virial radius. A significant fraction of this gas has been heated in supernova explosions in the disc and subsequently ejected into the CGM where it is predicted to give rise to substantial O VI absorption. The simulations do not yet address the question of what happens to the O VI when the galaxies stop forming stars. Our models also predict a reservoir of cool H I clouds that show strong Lyαabsorption to several hundred kpc. Comparing these models to recent surveys with the Hubble Space Telescope, we find that only the MF models have sufficient O VI and H I gas in the CGM to reproduce the observed distributions. In separate analyses, these same MF models also show better agreement with other galaxy observables (e.g. rotation curves, surface brightness profiles and H I gas distribution). We infer that the CGM is the dominant reservoir of baryons for galaxy haloes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The coupling between the core/cusp and missing satellite problems

We calculate the energy that baryons must inject in cold dark matter (CDM) haloes in order to remove centrally-divergent DM cusps on scales relevant to observations of dwarf spheroidal galaxies (dSphs). We estimate that the CDM haloes often associated with the Milky Way's dSphs (M_vir/M_\odot \sim 10^{9-10}) require ΔE/erg \sim 10^{53-55} in order to form cores on scales comparable to the luminous size of these galaxies. While supernova type II (SNeII) explosions can in principle generate this energy, the actual contribution is limited by the low star formation efficiency implied by the abundance of luminous satellites. Considering that CDM's well-known `core/cusp' and `missing satellite' problems place opposing demands on star formation efficiencies, existing observational evidences for large cores in the most luminous dSphs require that CDM models invoke some combination of the following: (i) efficient (of order unity) coupling of SNeII energy into dark matter particles, (ii) star formation histories peaking at unexpectedly high redshifts (z>6), (iii) a top-heavy stellar IMF, and/or (iv) substantial satellite disruption or other stochastic effects to ease the substructure abundance constraints. Our models show that the tension between CDM problems on small scales would increase if cored DM profiles were to be found in fainter dwarves.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Avoiding bias in reconstructing the largest observable scales from partial-sky data

Obscuration due to Galactic emission complicates the extraction of information from cosmological surveys, and requires some combination of the (typically imperfect) modeling and subtraction of foregrounds, or the removal of part of the sky. This particularly affects the extraction of information from the largest observable scales. Maximum-likelihood estimators for reconstructing the full-sky spherical harmonic coefficients from partial-sky maps have recently been shown to be susceptible to contamination from within the sky cut, arising due to the necessity to band-limit the data by smoothing prior to reconstruction. Using the WMAP 7-year data, we investigate modified implementations of such estimators which are robust to the leakage of contaminants from within masked regions. We provide a measure, based on the expected amplitude of residual foregrounds, for selecting the most appropriate estimator for the task at hand. We explain why the related quadratic maximum-likelihood estimator of the angular power spectrum does not suffer from smoothing-induced bias.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Linearization of homogeneous, nearly-isotropic cosmological models

Homogeneous, nearly-isotropic Bianchi cosmological models are considered. Their time evolution is expressed as a complete set of non-interacting linear modes on top of a Friedmann-Robertson-Walker background model. This connects the extensive literature on Bianchi models with the more commonly-adopted perturbation approach to general relativistic cosmological evolution. Expressions for the relevant metric perturbations in familiar coordinate systems can be extracted straightforwardly. Amongst other possibilities, this allows for future analysis of anisotropic matter sources in a more general geometry than usually attempted. We discuss the geometric mechanisms by which maximal symmetry is broken in the context of these models, shedding light on the origin of different Bianchi types. When all relevant length-scales are super-horizon, the simplest Bianchi I models emerge (in which anisotropic quantities appear parallel transported). Finally we highlight the existence of arbitrarily long near-isotropic epochs in models of general Bianchi type (including those without an exact isotropic limit).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The cut-sky cosmic microwave background is not anomalous

The observed angular correlation function of the cosmic microwave background has previously been reported to be anomalous, particularly when measured in regions of the sky uncontaminated by Galactic emission. Recent work by Efstathiou et al. presents a Bayesian comparison of isotropic theories, casting doubt on the significance of the purported anomaly. We extend this analysis to all anisotropic Gaussian theories with vanishing mean (<delta T> = 0), using the much wider class of models to confirm that the anomaly is not likely to point to new physics. On the other hand if there is any new physics to be gleaned, it results from low-l alignments which will be better quantified by a full-sky statistic. We also consider quadratic maximum likelihood power spectrum estimators that are constructed assuming isotropy. The underlying assumptions are therefore false if the ensemble is anisotropic. Nonetheless we demonstrate that, for theories compatible with the observed sky, these estimators (while no longer optimal) remain statistically superior to pseudo-C_l power spectrum estimators.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The nature of HI absorbers in GRB afterglows: clues from hydrodynamic simulations

In recent work, we have shown that it is possible to link quantitatively many aspects of damped Lyman alpha (DLA) absorbers in the spectra of quasars to high resolution simulations of galaxy formation. Using runs from the same series of hydrodynamic numerical studies, we consider the expected properties of Lyman alpha absorbers seen in the spectra of high redshift (z>2) gamma ray burst afterglows (GRB-DLAs). If GRBs are associated with the death of massive stars, their afterglows provide insights into otherwise unprobed regions of protogalactic objects, but detailed physical interpretations are currently embryonic. We find that median impact parameters (measured from the potential minimum) are approximately 1 kpc for GRBs compared with 4 kpc for QSO-DLAs. However, an equally important difference is that GRB-DLAs are predominantly associated with halos of mass 10^10<M_vir/M_sol<10^12, an order of magnitude larger than the hosts of QSO-DLAs. Accordingly, there are differences in the stellar properties of hosts. Our simulations accurately predict the form of the GRB-DLA HI column density distribution, producing quantitative agreement for N_HI>10^19 cm^-2, but they somewhat underpredict the incidence of low column densities N_HI<10^19 cm^-2. Line-of-sight neutral gas metallicities predicted by our simulations (10^-2 < Z/Z_sol < 1) are consistent with the modest observational constraints. Because of large internal dispersions in gas metallicities, this agreement is not significantly compromised by imposing a cut-off on the metallicity of stars able to launch GRBs (Z_star<Z_sol/3), confounding claims that the observed metallicity of GRB-DLAs poses a challenge to current GRB models. (Abridged.)

0 citations0 reviewsNo code linkedOpen access