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

L. Hernquist

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

ResearcherUnclaimed source record
astro-ph.GAastro-ph.COastro-phastro-ph.IM

Catalog footprint

What is connected

12works
4topics
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

12 published item(s)

preprint2025arXiv

Infalling ultra-faint dwarfs as emissaries of the Axiverse

Recent discoveries of ultra-faint dwarf galaxies (UFDs) infalling onto the Milky Way, namely Leo K \& M at $r \simeq 450$kpc, considerably strengthens the case that UFDs constitute a distinct galaxy class that is inherently smaller and fainter, and metal-poorer than the classical dwarf spheroidals (dSph). This distinction is at odds with the inherent continuity of galaxy halo masses formed under scale-free gravity for any standard dark-matter (DM) model. Here, we show that distinct galaxy classes do evolve in cosmological simulations of multiple light bosons representing the ``Axiverse'' proposal of string theory, where a discrete mass spectrum of axions is generically predicted to span many decades in mass. In this context, the observed UFD class we show corresponds to a relatively heavy boson of $3\times 10^{-21}$ eV, including Leo K \& M, whereas a lighter axion of $10^{-22}$ eV comprises the bulk of DM in all larger galaxies including the dSphs. Although Leo M is larger in size than Leo K, we predict its velocity dispersion to be smaller $(\simeq 1.7$km/s) than that of Leo K $(\simeq 4.5$km/s) because of the inverse de Broglie scale dependence on momentum. This scenario can be definitively tested using millisecond pulsars close to the Galactic center, where the Compton frequencies of the heavy and light bosons imprint monotone timing residuals that may be detected by the Square Kilometre Array (SKA) on timescales of approximately one week and four months, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Introducing the THESAN project: radiation-magnetohydrodynamic simulations of the epoch of reionization

We introduce the THESAN project, a suite of large volume (L = 95.5 cMpc) radiation-magnetohydrodynamic simulations that simultaneously model the large-scale statistical properties of the IGM during reionization and the resolved characteristics of the galaxies responsible for it. The flagship simulation has dark matter (DM) and baryonic mass resolutions of $3.1 \times 10^6$ $M_\odot$ and $5.8 \times 10^5 $ $M_\odot$, respectively. The gravitational forces are softened on scales of 2.2 ckpc with the smallest cell sizes reaching 10 pc at z=5.5, enabling predictions down to the atomic cooling limit. The simulations use an efficient radiation hydrodynamics solver (AREPO-RT) that precisely captures the interaction between ionizing photons and gas, coupled to well-tested galaxy formation (IllustrisTNG) and dust models to accurately predict the properties of galaxies. Through a complementary set of medium resolution simulations we investigate the changes to reionization introduced by different assumptions for ionizing escape fractions, varying DM models, and numerical convergence. The simulations produce realistic reionization histories that match the observed evolution of the global neutral hydrogen fraction and electron scattering optical depth to reionization. They also match a wealth of high-redshift observationally inferred data, including the stellar-to-halo-mass relation, stellar mass function, star formation rate density, and the mass-metallicity relation, despite the galaxy formation model being mainly calibrated at z=0. We demonstrate that different reionization models give rise to varied bubble size distributions that imprint unique signatures on the 21cm emission, especially on the slope of the power spectrum at large spatial scales, enabling current and upcoming 21cm experiments to accurately characterise the sources that dominate the ionizing photon budget. [abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The THESAN project: Lyman-alpha emission and transmission during the Epoch of Reionization

The visibility of high-redshift Lyman-alpha emitting galaxies (LAEs) provides important constraints on galaxy formation processes and the Epoch of Reionization (EoR). However, predicting realistic and representative statistics for comparison with observations represents a significant challenge in the context of large-volume cosmological simulations. The THESAN project offers a unique framework for addressing such limitations by combining state-of-the-art galaxy formation (IllustrisTNG) and dust models with the Arepo-RT radiation-magneto-hydrodynamics solver. In this initial study we present Lyman-alpha centric analysis for the flagship simulation that resolves atomic cooling haloes throughout a (95.5 cMpc)^3 region of the Universe. To avoid numerical artifacts we devise a novel method for accurate frequency-dependent line radiative transfer in the presence of continuous Hubble flow, transferable to broader astrophysical applications as well. Our scalable approach highlights the utility of LAEs and red damping-wing transmission as probes of reionization, which reveal nontrivial trends across different galaxies, sightlines, and frequency bands that can be modelled in the framework of covering fractions. In fact, after accounting for environmental factors influencing large-scale ionized bubble formation such as redshift and UV magnitude, the variation across galaxies and sightlines mainly depends on random processes including peculiar velocities and self-shielded systems that strongly impact unfortunate rays more than others. Throughout the EoR local and cosmological optical depths are often greater than or less than unity such that the exp(-tau) behavior leads to anisotropic and bimodal transmissivity. Future surveys will benefit by targeting both rare bright objects and Goldilocks zone LAEs to infer the presence of these (un)predictable (dis)advantages.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The stellar halos of ETGs in the IllustrisTNG simulations: the photometric and kinematic diversity of galaxies at large radii

We characterize the photometric and kinematic properties of simulated early-type galaxy (ETG) stellar halos, and compare them to observations. We select a sample of ~1200 ETGs in the TNG100 and TNG50 simulations, spanning a stellar mass range of $10^{10.3}-10^{12}M_{\odot}$ and within the range of (g-r) colour and lambda-ellipticity diagram populated by observed ETGs. We determine photometric parameters, intrinsic shapes, and kinematic observables in their extended stellar halos. We study the variation in kinematics from center to halo and connect it to a change in the intrinsic shape of the galaxies. We find that the simulated galaxy sample reproduces the diversity of kinematic properties observed in ETG halos. Simulated fast rotators (FRs) divide almost evenly in one third having flat lambda profiles and high halo rotational support, a third with gently decreasing profiles, and another third with low halo rotation. Slow rotators (SRs) tend to have increased rotation in the outskirts, with half of them exceeding lambda=0.2. For $M_{*}>10^{11.5}M_{\odot}$ halo rotation is unimportant. A similar variety of properties is found for the stellar halo intrinsic shapes. Rotational support and shape are deeply related: the kinematic transition to lower rotational support is accompanied by a change towards rounder intrinsic shape. Triaxiality in the halos of FRs increases outwards and with stellar mass. Simulated SRs have relatively constant triaxiality profiles. Simulated stellar halos show a large variety of structural properties, with quantitative but no clear qualitative differences between FRs and SRs. At the same stellar mass, stellar halo properties show a gradual transition and significant overlap between the two families, despite the clear bimodality in the central regions. This is in agreement with observations of extended photometry and kinematics. [abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Excitation of coupled stellar motions in the Galactic Disk by orbiting satellites

We use a set of high-resolution N-body simulations of the Galactic disk to study its interactions with the population of satellites predicted cosmologically. One simulation illustrates that multiple passages of massive satellites with different velocities through the disk generate a wobble, having the appearance of rings in face-on projections of the stellar disk. They also produce flares in the disk outer parts and gradually heat the disk through bending waves. A different numerical experiment shows that an individual satellite as massive as the Sagittarius dwarf galaxy passing through the disk will drive coupled horizontal and vertical oscillations of stars in underdense regions, with small significant associated heating. This experiment shows that vertical excursions of stars in these low-density regions can exceed 1 kpc in the Solar neighborhood, resembling the coherent vertical oscillations recently detected locally. They can also induce non-zero vertical streaming motions as large as 10-20 km s$^{-1}$, consistent with recent observations in the Galactic disk. This phenomenon appears as a local ring, with modest associated disk heating.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The information content of stellar halos: Stellar population gradients and accretion histories in early-type Illustris galaxies

Long dynamical timescales in the outskirts of galaxies preserve the information content of their accretion histories, for example in the form of stellar population gradients. We present a detailed analysis of the stellar halo properties of a statistically representative sample of early-type galaxies from the Illustris simulation and show that stellar population gradients at large radii can indeed be used to infer basic properties of galactic accretion histories. We measure metallicity, age, and surface-brightness profiles in quiescent Illustris galaxies ranging from $\mathrm{M}_\star = 10^{10} - 2\times 10^{12}\;\mathrm{M}_\odot$ and show that they are in reasonable agreement with observations. At fixed mass, galaxies that accreted little of their stellar halo material tend to have steeper metallicity and surface-brightness profiles between 2 - 4 effective radii (R$_e$) than those with larger accreted fractions. Profiles of metallicity and surface-brightness in the stellar halo typically flatten from z = 1 to the present. This suggests that the accretion of stars into the stellar halo tends to flatten metallicity and surface-brightness profiles, a picture which is supported by the tight correlation between the two gradients in the stellar halo. We find no statistical evidence of additional information content related to accretion histories in stellar halo metallicity profiles beyond what is contained in surface-brightness profiles. Age gradients in the stellar halo do not appear to be sensitive to galactic accretion histories, and none of the stellar population gradients studied are strongly correlated with the mean merger mass-ratio. Future observations that reach large radii outside galaxies will have the best potential to constrain galactic accretion histories.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Confirmation of Wide-Field Signatures in Redshifted 21 cm Power Spectra

We confirm our recent prediction of the "pitchfork" foreground signature in power spectra of high-redshift 21 cm measurements where the interferometer is sensitive to large-scale structure on all baselines. This is due to the inherent response of a wide-field instrument and is characterized by enhanced power from foreground emission in Fourier modes adjacent to those considered to be the most sensitive to the cosmological H I signal. In our recent paper, many signatures from the simulation that predicted this feature were validated against Murchison Widefield Array (MWA) data, but this key pitchfork signature was close to the noise level. In this paper, we improve the data sensitivity through the coherent averaging of 12 independent snapshots with identical instrument settings and provide the first confirmation of the prediction with a signal-to-noise ratio > 10. This wide-field effect can be mitigated by careful antenna designs that suppress sensitivity near the horizon. Simple models for antenna apertures that have been proposed for future instruments such as the Hydrogen Epoch of Reionization Array and the Square Kilometre Array indicate they should suppress foreground leakage from the pitchfork by ~40 dB relative to the MWA and significantly increase the likelihood of cosmological signal detection in these critical Fourier modes in the three-dimensional power spectrum.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Empirical Covariance Modeling for 21 cm Power Spectrum Estimation: A Method Demonstration and New Limits from Early Murchison Widefield Array 128-Tile Data

The separation of the faint cosmological background signal from bright astrophysical foregrounds remains one of the most daunting challenges of mapping the high-redshift intergalactic medium with the redshifted 21 cm line of neutral hydrogen. Advances in mapping and modeling of diffuse and point source foregrounds have improved subtraction accuracy, but no subtraction scheme is perfect. Precisely quantifying the errors and error correlations due to missubtracted foregrounds allows for both the rigorous analysis of the 21 cm power spectrum and for the maximal isolation of the "EoR window" from foreground contamination. We present a method to infer the covariance of foreground residuals from the data itself in contrast to previous attempts at a priori modeling. We demonstrate our method by setting limits on the power spectrum using a 3 h integration from the 128-tile Murchison Widefield Array. Observing between 167 and 198 MHz, we find at 95% confidence a best limit of Delta^2(k) < 3.7 x 10^4 mK^2 at comoving scale k = 0.18 hMpc^-1 and at z = 6.8, consistent with existing limits.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Stellar populations of stellar halos: Results from the Illustris simulation

The influence of both major and minor mergers is expected to significantly affect gradients of stellar ages and metallicities in the outskirts of galaxies. Measurements of observed gradients are beginning to reach large radii in galaxies, but a theoretical framework for connecting the findings to a picture of galactic build-up is still in its infancy. We analyze stellar populations of a statistically representative sample of quiescent galaxies over a wide mass range from the Illustris simulation. We measure metallicity and age profiles in the stellar halos of quiescent Illustris galaxies ranging in stellar mass from $10^{10}$ to $10^{12} M_\odot$, accounting for observational projection and luminosity-weighting effects. We find wide variance in stellar population gradients between galaxies of similar mass, with typical gradients agreeing with observed galaxies. We show that, at fixed mass, the fraction of stars born in-situ within galaxies is correlated with the metallicity gradient in the halo, confirming that stellar halos contain unique information about the build-up and merger histories of galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Lyman-alpha Cooling Emission from Galaxy Formation

Recent studies have shown that galaxies accrete most of their baryons via the cold mode, from streams with temperatures T~10^4-10^5 K. At these temperatures, the streams should radiate primarily in the Lya line and have therefore been proposed as a model to power the extended, high-redshift objects known as Lya blobs and other high-redshift Lya sources. We introduce a new Lya radiative transfer code, aRT, and apply it to cosmological hydrodynamical simulations. We address physical and numerical issues that are critical to making accurate predictions for the cooling luminosity, but that have been mostly neglected or treated simplistically so far. We highlight the importance of self-shielding and of properly treating sub-resolution models in simulations. Most existing simulations do not self-consistently incorporate these effects, which can lead to order-of-magnitude errors in the predicted cooling luminosity. Using a combination of post-processing ionizing radiative transfer and re-simulation techniques, we develop an approximation to the consistent evolution of the self-shielded gas. We quantify the dependence of the Lya cooling luminosity on halo mass at z=3 for the simplified problem of pure gas accretion. While cooling in massive halos (without additional energy input from star formation and AGN) is in principle sufficient to produce L_alpha~10^43-10^44 erg s^-1 blobs, this requires including energy released in gas of density sufficient to form stars, but which is kept 100% gaseous in our optimistic estimates. Excluding emission from such dense gas yields lower luminosities by up to one to two orders of magnitude at high masses, making it difficult to explain the observed Lya blobs with pure cooling. Resonant scattering produces diffuse Lya halos, even for centrally concentrated emission, and broad double peaked line profiles. [Abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

A Century of Cosmology: A Direct Precision Measurement of the Intergalactic Lyman-alpha Opacity at 2<z<4.2

We directly measure the evolution of the intergalactic Lyman-alpha effective optical depth, tau_eff, over the redshift range 2<z<4.2 from a sample of 86 high-resolution, high-signal-to-noise quasar spectra obtained with Keck/ESI, Keck/HIRES, and Magellan/MIKE. We find that our estimates of the quasar continuum levels in the Ly-alpha forest obtained by spline fitting are systematically biased low, but that this bias can be accounted for using mock spectra. The mean fractional error <Delta C/C_true> is <1% at z=2, 4% at z=3, and 12% at z=4. We provide estimates of the level of absorption arising from metals in the Ly-alpha forest based on both direct and statistical metal removal results in the literature, finding that this contribution is ~6-9% at z=3 and decreases monotonically with redshift. The high precision of our measurement indicates significant departures from the best-fit power-law redshift evolution, particularly near z=3.2.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Extended Perturbation Theory for the Local Density Distribution Function

Perturbation theory makes it possible to calculate the probability distribution function (PDF) of the large scale density field in the small variance limit. For top hat smoothing and scale-free Gaussian initial fluctuations, the result depends only on the linear variance, sigma_linear, and its logarithmic derivative with respect to the filtering scale -(n_linear+3)=dlog sigma_linear^2/dlog L (Bernardeau 1994). In this paper, we measure the PDF and its low-order moments in scale-free simulations evolved well into the nonlinear regime and compare the results with the above predictions, assuming that the spectral index and the variance are adjustable parameters, n_eff and sigma_eff=sigma, where sigma is the true, nonlinear variance. With these additional degrees of freedom, results from perturbation theory provide a good fit of the PDFs, even in the highly nonlinear regime. The value of n_eff is of course equal to n_linear when sigma << 1, and it decreases with increasing sigma. A nearly flat plateau is reached when sigma >> 1. In this regime, the difference between n_eff and n_linear increases when n_linear decreases. For initial power-spectra with n_linear=-2,-1,0,+1, we find n_eff ~ -9,-3,-1,-0.5 when sigma^2 ~ 100.

0 citations0 reviewsNo code linkedOpen access