Source author record

J. D. Emberson

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

ResearcherUnclaimed source record

astro-ph.CO astro-ph.GA gr-qc

Catalog footprint

What is connected

6works

3topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2026arXiv

Emulator-Based Inference of Cosmological Subgrid Models

The formation of structure in the Universe at large scales is dominated by gravity, with baryonic physics becoming significant at $\sim{\rm Mpc}$ scales. To capture the impact of baryonic physics, cosmological simulations must model gas dynamics and a host of relevant astrophysical processes. A recent extension of the Hardware/Hybrid Accelerated Cosmology Code (HACC) couples its gravity solver with a modern smoothed particle hydrodynamics method. This extension incorporates sub-resolution models for chemical enrichment, black hole and star formation, AGN kinetic and thermal feedback, supernova-driven feedback, galactic winds, and metal-line cooling. We present an inference framework based on high-fidelity emulators to aid in model calibration against observational targets, e.g., the galaxy stellar mass function, radial gas density profiles, and the cluster gas fraction. The emulators are trained on simulation suites comprising 64 boxes with side-length $128\,h^{-1}$Mpc and 16 boxes with side-length $256\,h^{-1}$Mpc with $2\times 512^3$ and $2\times 1024^3$ particles, respectively. Our analysis reveals two distinct AGN kinetic feedback modes -- a low-feedback mode yielding strong agreement with the observed radial gas density profiles of massive X-ray clusters, and a high-feedback mode providing a better fit to cluster gas fraction data, but systematically underestimating gas densities in inner regions.

preprint2015arXiv

Evolution of Low Mass Galactic Subhalos and Dependence on Concentration

We carry out a detailed study of the orbital dynamics and structural evolution of over 6000 subhalos in the Via Lactea II simulation, from infall to present. By analyzing subhalos with masses down to m = 4e5 Msun, we find that lower mass subhalos, which are not strongly affected by dynamical friction, exhibit behaviors qualitatively different from those found previously for more massive ones. Furthermore, there is a clear trend of subhalos that fell into the host earlier being less concentrated. We show that the concentration at infall characterizes various aspects of subhalo evolution. In particular, tidal effects truncate the growth of less concentrated subhalos at larger distances from the host; subhalos with smaller concentrations have larger infall radii. The concentration at infall is further shown to be a determining factor for the subsequent mass loss of subhalos within the host, and also for the evolution of their internal structure in the v_max-r_max plane. Our findings raise the prospects of using the concentration to predict the tidal evolution of subhalos, which will be useful for obtaining analytic models of galaxy formation, as well as for near field cosmology.

preprint2015arXiv

Precision reconstruction of the dark matter-neutrino relative velocity from N-body simulations

Discovering the mass of neutrinos is a principle goal in high energy physics and cosmology. In addition to cosmological measurements based on two-point statistics, the neutrino mass can also be estimated by observations of neutrino wakes resulting from the relative motion between dark matter and neutrinos. Such a detection relies on an accurate reconstruction of the dark matter-neutrino relative velocity which is affected by non-linear structure growth and galaxy bias. We investigate our ability to reconstruct this relative velocity using large N-body simulations where we evolve neutrinos as distinct particles alongside the dark matter. We find that the dark matter velocity power spectrum is overpredicted by linear theory whereas the neutrino velocity power spectrum is underpredicted. The magnitude of the relative velocity observed in the simulations is found to be lower than what is predicted in linear theory. Since neither the dark matter nor the neutrino velocity fields are directly observable from galaxy or 21 cm surveys, we test the accuracy of a reconstruction algorithm based on halo density fields and linear theory. Assuming prior knowledge of the halo bias, we find that the reconstructed relative velocities are highly correlated with the simulated ones with correlation coefficients of 0.94, 0.93, 0.91 and 0.88 for neutrinos of mass 0.05, 0.1, 0.2 and 0.4 eV. We confirm that the relative velocity field reconstructed from large scale structure observations such as galaxy or 21 cm surveys can be accurate in direction and, with appropriate scaling, magnitude.

preprint2013arXiv

High Performance P3M N-body code: CUBEP3M

This paper presents CUBEP3M, a publicly-available high performance cosmological N-body code and describes many utilities and extensions that have been added to the standard package. These include a memory-light runtime SO halo finder, a non-Gaussian initial conditions generator, and a system of unique particle identification. CUBEP3M is fast, its accuracy is tuneable to optimize speed or memory, and has been run on more than 27,000 cores, achieving within a factor of two of ideal weak scaling even at this problem size. The code can be run in an extra-lean mode where the peak memory imprint for large runs is as low as 37 bytes per particles, which is almost two times leaner than other widely used N-body codes. However, load imbalances can increase this requirement by a factor of two, such that fast configurations with all the utilities enabled and load imbalances factored in require between 70 and 120 bytes per particles. CUBEP3M is well designed to study large scales cosmological systems, where imbalances are not too large and adaptive time-stepping not essential. It has already been used for a broad number of science applications that require either large samples of non-linear realizations or very large dark matter N-body simulations, including cosmological reionization, halo formation, baryonic acoustic oscillations, weak lensing or non-Gaussian statistics. We discuss the structure, the accuracy, known systematic effects and the scaling performance of the code and its utilities, when applicable.

preprint2013arXiv

The Opacity of the Intergalactic Medium During Reionization: Resolving Small-Scale Structure

Early in the reionization process, the intergalactic medium (IGM) would have been quite inhomogeneous on small scales, due to the low Jeans mass in the neutral IGM and the hierarchical growth of structure in a cold dark matter Universe. This small-scale structure acted as an important sink during the epoch of reionization, impeding the progress of the ionization fronts that swept out from the first sources of ionizing radiation. Here we present results of high-resolution cosmological hydrodynamics simulations that resolve the cosmological Jeans mass of the neutral IGM in representative volumes several Mpc across. The adiabatic hydrodynamics we follow are appropriate in an unheated IGM, before the gas has had a chance to respond to the photoionization heating. Our focus is determination of the resolution required in cosmological simulations in order to sufficiently sample and resolve small-scale structure regulating the opacity of an unheated IGM. We find that a dark matter particle mass of m_dm < 50 Msun and box size of L > 1 Mpc are required. With our converged results we show how the mean free path of ionizing radiation and clumping factor of ionized hydrogen depends upon the ultraviolet background (UVB) flux and redshift. We find, for example at z = 10, clumping factors typically of 10 to 20 for an ionization rate of Gamma ~ 0.3 - 3 x 1e-12 s^-1, with corresponding mean free paths of ~ 3 - 15 Mpc, extending previous work on the evolving mean free path to considerably smaller scales and earlier times.

preprint2012arXiv

Interpolation in waveform space: enhancing the accuracy of gravitational waveform families using numerical relativity

Matched-filtering for the identification of compact object mergers in gravitational-wave antenna data involves the comparison of the data stream to a bank of template gravitational waveforms. Typically the template bank is constructed from phenomenological waveform models since these can be evaluated for an arbitrary choice of physical parameters. Recently it has been proposed that singular value decomposition (SVD) can be used to reduce the number of templates required for detection. As we show here, another benefit of SVD is its removal of biases from the phenomenological templates along with a corresponding improvement in their ability to represent waveform signals obtained from numerical relativity (NR) simulations. Using these ideas, we present a method that calibrates a reduced SVD basis of phenomenological waveforms against NR waveforms in order to construct a new waveform approximant with improved accuracy and faithfulness compared to the original phenomenological model. The new waveform family is given numerically through the interpolation of the projection coefficients of NR waveforms expanded onto the reduced basis and provides a generalized scheme for enhancing phenomenological models.