Engineering topics

Hernquist, Lars

Publications and source records attributed to Hernquist, Lars.

At least 19 records

DOE OSTI2023

Interpreting Sunyaev–Zel’dovich observations with MillenniumTNG: mass and environment scaling relations

ABSTRACT Sunyaev–Zel’dovich (SZ) measurements can dramatically improve our understanding of the intergalactic medium and the role of feedback processes in galaxy formation, allowing us to calibrate important astrophysical systematics in cosmological constraints from weak lensing galaxy clustering surveys. However, the signal is only measured in a two-dimensional projection, and its correct interpretation relies on understanding the connection between observable quantities and the underlying intrinsic properties of the gas, in addition to the relation between the gas and the underlying matter distribution. One way to address these challenges is through the use of hydrodynamical simulations such as the high-resolution, large-volume MillenniumTNG suite. We find that measurements of the optical depth, τ, and the Compton-y parameter, Y, receive large line-of-sight contributions that can be removed effectively by applying a compensated aperture photometry filter. In contrast with other τ probes (e.g. X-rays and fast radio bursts), the kinematic SZ-inferred τ receives most of its signal from a confined cylindrical region around the halo due to the velocity decorrelation along the line of sight. Additionally, we perform fits to the Y–M and τ–M scaling relations and report best-fitting parameters adopting the smoothly broken power law formalism. We note that subgrid physics modelling can broaden the error bar on these by 30 per cent for intermediate-mass haloes (${\sim }10^{13} \, {\rm M}_{\odot }$). The scatter of the scaling relations can be captured by an intrinsic dependence on concentration and an extrinsic dependence on tidal shear. Finally, we comment on the effect of using galaxies rather than haloes in observations, which can bias the inferred profiles by ∼20 per cent for L* galaxies.

Hadzhiyska, Boryana (ORCID:0000000223123121)↗

DOE OSTI2023

NANCY: Next-generation All-sky Near-infrared Community surveY

The Nancy Grace Roman Space Telescope is capable of delivering an unprecedented all-sky, high-spatial resolution, multi-epoch infrared map to the astronomical community. This opportunity arises in the midst of numerous ground- and space-based surveys that will provide extensive spectroscopy and imaging together covering the entire sky (such as Rubin/LSST, Euclid, UNIONS, SPHEREx, DESI, SDSS-V, GALAH, 4MOST, WEAVE, MOONS, PFS, UVEX, NEO Surveyor, etc.). Roman can uniquely provide uniform high-spatial-resolution (~0.1 arcsec) imaging over the entire sky, vastly expanding the science reach and precision of all of these near-term and future surveys. This imaging will not only enhance other surveys, but also facilitate completely new science. By imaging the full sky over two epochs, Roman can measure the proper motions for stars across the entire Milky Way, probing 100 times fainter than Gaia out to the very edge of the Galaxy. Here, we propose NANCY: a completely public, all-sky survey that will create a high-value legacy dataset benefiting innumerable ongoing and forthcoming studies of the universe. NANCY is a pure expression of Roman's potential: it images the entire sky, at high spatial resolution, in a broad infrared bandpass that collects as many photons as possible. The majority of all ongoing astronomical surveys would benefit from incorporating observations of NANCY into their analyses, whether these surveys focus on nearby stars, the Milky Way, near-field cosmology, or the broader universe.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2023

The first quiescent galaxies in TNG300

We identify the first quiescent galaxies in TNG300, the largest volume of the IllustrisTNG cosmological simulation suite, and explore their quenching processes and time evolution to z = 0. We find that the first quiescent galaxies with stellar masses M_* > 3 × 10^10 M_⊙ and specific star formation rates sSFR < 10^−11 yr^−1 emerge at z ∼ 4.2 in TNG300. Suppression of star formation in these galaxies begins with a thermal mode of active galactic nucleus feedback at z ∼ 6, and a kinetic feedback mode acts in each galaxy by z ∼ 4.7 to complete the quenching process, which occurs on a time-scale of ∼0.35 Gyr. Surprisingly, we find that the majority of these galaxies are not the main progenitors of their z = 0 descendants; instead, four of the five galaxies fall into more massive galaxies in subsequent mergers at a range of redshifts 2.5 < z < 0.2. By z = 0, these descendants are the centres of galaxy clusters with average stellar masses of 8 × 10^11 M_⊙. We make predictions for the first quenched galaxies to be located by the JWST.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Cosmological structure formation and soliton phase transition in fuzzy dark matter with axion self-interactions

Herein we investigate cosmological structure formation in fuzzy dark matter (FDM) with the attractive self-interaction (SI) with numerical simulations. Such a SI would arise if the FDM boson were an ultra-light axion, which has a strong CP symmetry-breaking scale (decay constant). Although weak, the attractive SI may be strong enough to counteract the quantum ‘pressure’ and alter structure formation. We find in our simulations that the SI can enhance small-scale structure formation, and soliton cores above a critical mass undergo a phase transition, transforming from dilute to dense solitons.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

H α emission in local galaxies: star formation, time variability, and the diffuse ionized gas

ABSTRACT The nebular recombination line H α is widely used as a star formation rate (SFR) indicator in the local and high-redshift Universe. We present a detailed H α radiative transfer study of high-resolution isolated Milky-Way and Large Magellanic Cloud simulations that include radiative transfer, non-equilibrium thermochemistry, and dust evolution. We focus on the spatial morphology and temporal variability of the H α emission, and its connection to the underlying gas and star formation properties. The H α and H β radial and vertical surface brightness profiles are in excellent agreement with observations of nearby galaxies. We find that the fraction of H α emission from collisional excitation amounts to fcol ∼ 5–$10{{\ \rm per\ cent}}$, only weakly dependent on radius and vertical height, and that scattering boosts the H α luminosity by $\sim 40{{\ \rm per\ cent}}$. The dust correction via the Balmer decrement works well (intrinsic H α emission recoverable within 25 per cent), though the dust attenuation law depends on the amount of attenuation itself both on spatially resolved and integrated scales. Important for the understanding of the H α–SFR connection is the dust and helium absorption of ionizing radiation (Lyman continuum [LyC] photons), which are about $f_{\rm abs}\approx 28{{\ \rm per\ cent}}$ and $f_{\rm He}\approx 9{{\ \rm per\ cent}}$, respectively. Together with an escape fraction of $f_{\rm esc}\approx 6{{\ \rm per\ cent}}$, this reduces the available budget for hydrogen line emission by nearly half ($f_{\rm H}\approx 57{{\ \rm per\ cent}}$). We discuss the impact of the diffuse ionized gas, showing – among other things – that the extraplanar H α emission is powered by LyC photons escaping the disc. Future applications of this framework to cosmological (zoom-in) simulations will assist in the interpretation of spectroscopy of high-redshift galaxies with the upcoming James Webb Space Telescope.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Galaxy assembly bias and large-scale distribution: a comparison between IllustrisTNG and a semi-analytic model

In this work, we compare large scale structure observables for stellar mass selected samples at z = 0, as predicted by two galaxy models, the hydrodynamical simulation IllustrisTNG and the Santa-Cruz semi-analytic model (SC-SAM). Although both models have been independently calibrated to match observations, rather than each other, we find good agreement between the two models for two-point clustering and galaxy assembly bias signatures. The models also show a qualitatively similar response of occupancy and clustering to secondary halo parameters other than mass, such as formation history and concentration, although with some quantitative differences. Thus, our results demonstrate that the galaxy–halo relationships in SC-SAM and TNG are quite similar to first order. However, we also find areas in which the models differ. For example, we note a strong correlation between halo gas content and environment in TNG, which is lacking in the SC-SAM, as well as differences in the occupancy predictions for low-mass haloes. Moreover, we show that higher order statistics, such as cumulants of the density field, help us to accurately describe the galaxy distribution and discriminate between models that show degenerate behaviour for two-point statistics. Our results suggest that SAMs are a promising cost-effective and intuitive method for generating mock catalogues for next generation cosmological surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Extensions to models of the galaxy–halo connection

ABSTRACT We explore two widely used empirical models for the galaxy–halo connection, subhalo abundance matching (SHAM) and the halo occupation distribution (HOD), and compare them with the hydrodynamical simulation IllustrisTNG (TNG) for multiple statistics quantifying the galaxy distribution at $n_{\rm gal}\approx 1.3\times 10^{-3}\, ({\rm Mpc}\,h^{-1})^{-3}$. We observe that in their most straightforward implementations, both models fail to reproduce the two-point clustering measured in TNG. We find that SHAM models that use the relaxation velocity, Vrelax, and the peak velocity, Vpeak, perform best, and match the clustering reasonably well, although neither captures adequately the one-halo clustering. Splitting the total sample into sub-populations, we discover that SHAM overpredicts the clustering of high-mass, blue, star-forming, and late-forming galaxies and underpredicts that of low-mass, red, quiescent, and early-forming galaxies. We also study various baryonic effects, finding that subhaloes in the dark-matter-only simulation have consistently higher values of their SHAM-proxy properties than their full-physics counterparts. We then consider a 2D implementation of the HOD model augmented with a secondary parameter (environment, velocity anisotropy, σ2Rhalf-mass, and total potential) tuned so as to match the two-point clustering of the IllustrisTNG galaxies on large scales. We analyse these galaxy populations adopting alternative statistical tools such as galaxy–galaxy lensing, void–galaxy cross-correlations, and cumulants of the density field, finding that the hydrodynamical galaxy distribution disfavours σ2Rhalf-mass and the total potential as secondary parameters, while the environment and velocity anisotropy samples are consistent with full physics across all statistical probes examined. Our results demonstrate the power of examining multiple statistics for determining the secondary parameters that are vital for understanding the galaxy–halo connection.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Stellar and weak lensing profiles of massive galaxies in the Hyper-Suprime Cam survey and in hydrodynamic simulations

ABSTRACT We perform a consistent comparison of the mass and mass profiles of massive (M⋆ > 1011.4 M⊙) central galaxies at z ∼ 0.4 from deep Hyper Suprime-Cam (HSC) observations and from the Illustris, TNG100, and Ponos simulations. Weak lensing measurements from HSC enable measurements at fixed halo mass and provide constraints on the strength and impact of feedback at different halo mass scales. We compare the stellar mass function (SMF) and the Stellar-to-Halo Mass Relation (SHMR) at various radii and show that the radius at which the comparison is performed is important. In general, Illustris and TNG100 display steeper values of α where $M_{\star } \propto M_{\rm vir}^{\alpha }$. These differences are more pronounced for Illustris than for TNG100 and in the inner rather than outer regions of galaxies. Differences in the inner regions may suggest that TNG100 is too efficient at quenching in situ star formation at Mvir ≃ 1013 M⊙ but not efficient enough at Mvir ≃ 1014 M⊙. The outer stellar masses are in excellent agreement with our observations at Mvir ≃ 1013 M⊙, but both Illustris and TNG100 display excess outer mass as Mvir ≃ 1014 M⊙ (by ∼0.25 and ∼0.12 dex, respectively). We argue that reducing stellar growth at early times in $M_\star \sim 10^{9-10} \, \mathrm{M}_{\odot }$ galaxies would help to prevent excess ex-situ growth at this mass scale. The Ponos simulations do not implement AGN feedback and display an excess mass of ∼0.5 dex at r < 30 kpc compared to HSC which is indicative of overcooling and excess star formation in the central regions. The comparison of the inner profiles of Ponos and HSC suggests that the physical scale over which the central AGN limits star formation is r ≲ 20 kpc. Joint comparisons between weak lensing and galaxy stellar profiles are a direct test of whether simulations build and deposit galaxy mass in the correct dark matter haloes and thereby provide powerful constraints on the physics of feedback and galaxy growth. Our galaxy and weak lensing profiles are publicly available to facilitate comparisons with other simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Joint galaxy–galaxy lensing and clustering constraints on galaxy formation

ABSTRACT We compare predictions for galaxy–galaxy lensing profiles and clustering from the Henriques et al. public version of the Munich semi-analytical model (SAM) of galaxy formation and the IllustrisTNG suite, primarily TNG300, with observations from KiDS + GAMA and SDSS-DR7 using four different selection functions for the lenses (stellar mass, stellar mass and group membership, stellar mass and isolation criteria, and stellar mass and colour). We find that this version of the SAM does not agree well with the current data for stellar mass-only lenses with $M_\ast \gt 10^{11}\, \mathrm{ M}_\odot$. By decreasing the merger time for satellite galaxies as well as reducing the radio-mode active galactic nucleus accretion efficiency in the SAM, we obtain better agreement, both for the lensing and the clustering, at the high-mass end. We show that the new model is consistent with the signals for central galaxies presented in Velliscig et al. Turning to the hydrodynamical simulation, TNG300 produces good lensing predictions, both for stellar mass-only (χ2 = 1.81 compared to χ2 = 7.79 for the SAM) and locally brightest galaxy samples (χ2 = 3.80 compared to χ2 = 5.01). With added dust corrections to the colours it matches the SDSS clustering signal well for red low-mass galaxies. We find that both the SAMs and TNG300 predict $\sim 50\, {{\ \rm per\ cent}}$ excessive lensing signals for intermediate-mass red galaxies with 10.2 < log10M*[M⊙] < 11.2 at $r \approx 0.6\, h^{-1}\, \text{Mpc}$, which require further theoretical development.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Limitations to the ‘basic’ HOD model and beyond

ABSTRACT We use the IllustrisTNG cosmological, hydrodynamical simulations to test fundamental assumptions of the mass-based halo occupation distribution (HOD) approach to modelling the galaxy–halo connection. By comparing the clustering of galaxies measured in the 300 Mpc TNG box (TNG300) with that predicted by the standard (basic) HOD model, we find that, on average, the ‘basic’ HOD model underpredicts the real-space correlation function in the TNG300 box by ∼15 per cent on scales of $1 \,\,\lt\,\, r \,\,\lt\,\, 20 \ {\rm Mpc}\, h^{-1}$, which is well beyond the target precision demanded of next-generation galaxy redshift surveys. We perform several tests to establish the robustness of our findings to systematic effects, including the effect of finite box size and the choice of halo finder. In our exploration of ‘secondary’ parameters with which to augment the ‘basic’ HOD, we find that the local environment of the halo, the velocity dispersion anisotropy, β, and the product of the half-mass radius and the velocity dispersion, σ2Rhalfmass, are the three most effective measures of assembly bias that help reconcile the ‘basic’ HOD-predicted clustering with that in TNG300. In addition, we test other halo properties such as halo spin, formation epoch, and halo concentration. We also find that at fixed halo mass, galaxies in one type of environment cluster differently from galaxies in another. We demonstrate that a more complete model of the galaxy–halo connection can be constructed if we combine both mass and local environment information about the halo.

42 ENGINEERING↗

DOE OSTI2019

First Star-Forming Structures in Fuzzy Cosmic Filaments

In hierarchical models of structure formation, the first galaxies form in low-mass dark matter potential wells, probing the behavior of dark matter on kiloparsec scales. Even though these objects are below the detection threshold of current telescopes, future missions will open an observational window into this emergent world. In this Letter, we investigate how the first galaxies are assembled in a “fuzzy” dark matter (FDM) cosmology where dark matter is an ultralight ~10 –22 eV boson and the primordial stars are expected to form along dense dark matter filaments. Furthermore, using a first-of-its-kind cosmological hydrodynamical simulation, we explore the interplay between baryonic physics and unique wavelike features inherent to FDM. In our simulation, the dark matter filaments show coherent interference patterns on the boson de Broglie scale and develop cylindrical solitonlike cores, which are unstable under gravity and collapse into kiloparsec-scale spherical solitons. Features of the dark matter distribution are largely unaffected by the baryonic feedback. On the contrary, the distributions of gas and stars, which do form along the entire filament, exhibit central cores imprinted by dark matter—a smoking gun signature of FDM.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1995

Double core evolution. 7: The infall of a neutron star through the envelope of its massive star companion

Binary systems with properties similar to those of high-mass X-ray binaries are evolved through the common envelope phase. Three-dimensional simulations show that the timescale of the infall phase of the neutron star depends upon the evolutionary state of its massive companion. We find that tidal torques more effectively accelerate common envelope evolution for companions in their late core helium-burning stage and that the infall phase is rapid (approximately several initial orbital periods). For less evolved companions the decay of the orbit is longer; however, once the neutron star is deeply embedded within the companion's envelope the timescale for orbital decay decreases rapidly. As the neutron star encounters the high-density region surrounding the helium core of its massive companion, the rate of energy loss from the orbit increases dramatically leading to either partial or nearly total envelope ejection. The outcome of the common envelope phase depends upon the structure of the evolved companion. In particular, it is found that the entire common envelope can be ejected by the interaction of the neutron star with a red supergiant companion in binaries with orbital periods similar to those of long-period Be X-ray binaries. For orbital periods greater than or approximately equal to 0.8-2 yr (for companions of mass 12-24 solar mass) it is likely that a binary will survive the common envelope phase. For these systems, the structure of the progenitor star is characterized by a steep density gradient above the helium core, and the common envelope phase ends with a spin up of the envelope to within 50%-60% of corotation and with a slow mass outflow. The efficiency of mass ejection is found to be approximately 30%-40%. For less evolved companions, there is insufficient energy in the orbit to unbind the common envelope and only a fraction of it is ejected. Since the timescale for orbital decay is always shorter than the mass-loss timescale from the common envelope, the two cores will likely merge to form a Thorne-Zytkow object. Implications for the origin of Cyg X-3, an X-ray source consisting of a Wolf-Rayet star and a compact companion, and for the fate of the remnant binary consisting of a helium star and a neutron star are briefly discussed.

Terman, James L.↗

NASA NTRS1995

Using X-rays to determine which compact groups are illusory

If the large-scale galaxy distribution is filamentary, as suggested by some observations and recent hydrodynamical simulations, then lengthwise views of filaments will apparently produce compact groups (CGs) that are in reality stretched out along the line of sight. This possibility has been advocated recently by Hernquist, Katz, & Weinberg (1995). Here we propose a test for this hypothesis using X-ray emission from CGs. The observable quantity Q identical with L(sub x)a(sup 3)(sup p)/(L(sup 2)(sup g))T(sup 1/2)(sub x) should be proportional to the axis ratio of the group, a/c, where a and c are the short and long axes of a prolate distribution, a(sub p) is the radius of the group projected onto the sky, L(sub x) is the bolometric X-ray luminosity, L(sub g) is the group blue luminosity, and T(sub x) is the gas temperature. We find that the distribution of Q is consistent with the notion that many CGs with unusually small values of a/c are frauds, i.e., that the values of Q are anomalously small. Absent other information, it is equally possible that CGs are very gas-poor relative to rich clusters; however, this can be tested using the Sunyaev-Zeldovich effect. If the groups have a close to normal ratio of gas to total mass, but are simply stretched out along the line of sight, a Sunyaev-Zeldovich signal should be detectable.

Ostriker, Jeremiah P.↗

NASA NTRS1995

Physically detached 'compact groups'

A small fraction of galaxies appear to reside in dense compact groups, whose inferred crossing times are much shorter than a Hubble time. These short crossing times have led to considerable disagreement among researchers attempting to deduce the dynamical state of these systems. In this paper, we suggest that many of the observed groups are not physically bound but are chance projections of galaxies well separated along the line of sight. Unlike earlier similar proposals, ours does not require that the galaxies in the compact group be members of a more diffuse, but physically bound entity. The probability of physically separated galaxies projecting into an apparent compact group is nonnegligible if most galaxies are distributed in thin filaments. We illustrate this general point with a specific example: a simulation of a cold dark matter universe, in which hydrodynamic effects are included to identify galaxies. The simulated galaxy distribution is filamentary and end-on views of these filaments produce apparent galaxy associations that have sizes and velocity dispersions similar to those of observed compact groups. The frequency of such projections is sufficient, in principle, to explain the observed space density of groups in the Hickson catalog. We discuss the implications of our proposal for the formation and evolution of groups and elliptical galaxies. The proposal can be tested by using redshift-independent distance estimators to measure the line-of-sight spatial extent of nearby compact groups.

Hernquist, Lars↗

NASA NTRS1995

A comparison of two algorithms for simulating collisionless systems

Two completely different simulation algorithms are compared by applying them to the same stellar dynamical problems: one is a self-consistent field (SCF) method for solving Poisson's equation and the other is a phase-space method for integrating the collisionless Boltzmann equation. We consider simulations of spherical stellar systems which are initially far from equilibrium and relax to their final states by gravitational collapse. The initial conditions consist of either uniform-density spheres or nonequilibrium models having Plummer density profiles, in which velocity dispersions are assigned according to given virial ratios. If a few tens of radial expansion terms with hundreds of thousands of particles are used in the SCF code, excellent agreement is found between the results it generates and those obtained with the phase-space solver, provided that a sufficiently large number of grid cells are employed with the latter. These findings imply that for simulating collisionless systmes over many dynamical times, the SCF approach based on sampling phase space is competitive with the approach treating phase space as a continuous fluid. The results of our tests make it possible to estimate the number of particles and basis functions required in situations like those modeled. Limitations of the SCF method and the choice of an optimal set of basis functions are also discussed.

Hozumi, Shunsuke↗

NASA NTRS1994

Star-forming galaxy models: Blending star formation into TREESPH

We have incorporated star-formation algorithms into a hybrid N-body/smoothed particle hydrodynamics code (TREESPH) in order to describe the star forming properties of disk galaxies over timescales of a few billion years. The models employ a Schmidt law of index n approximately 1.5 to calculate star-formation rates, and explicitly include the energy and metallicity feedback into the Interstellar Medium (ISM). Modeling the newly formed stellar population is achieved through the use of hybrid SPH/young star particles which gradually convert from gaseous to collisionless particles, avoiding the computational difficulties involved in creating new particles. The models are shown to reproduce well the star-forming properties of disk galaxies, such as the morphology, rate of star formation, and evolution of the global star-formation rate and disk gas content. As an example of the technique, we model an encounter between a disk galaxy and a small companion which gives rise to a ring galaxy reminiscent of the Cartwheel (AM 0035-35). The primary galaxy in this encounter experiences two phases of star forming activity: an initial period during the expansion of the ring, and a delayed phase as shocked material in the ring falls back into the central regions.

Mihos, J. Christopher↗

NASA NTRS1994

Dense stellar cores in merger remnants

We use numerical models which include star formation to analyze the mass profiles of remnants formed by mergers of disk galaxies. During a merger, dissipation in gas and ensuing star formation leave behind a dense stellar core in the remnant. Rather than joining smoothly onto a de Vaucouleurs profile, the starburst population leads to a sharp break in the surface density profile at a few percent of the effective radius. While our results are preliminary, the lack of such signatures in most elliptical galaxies suggests that mergers of gas-rich disk galaxies may not have contributed greatly to the populations of present-day ellipticals.

Mihos, J. Christopher↗

NASA NTRS1994

Kinematic misalignments in remnants of multiple mergers

Observations indicate that the rotation and minor axes in elliptical galaxies are usually well aligned. However, computer models of dissipationless mergers of pairs of disk galaxies have failed to establish such a correlation, implying that most ellipticals probably did not originate in this manner. Here we consider remnants produced from repeated merging in small groups of a half-dozen or so disk galaxies and find that the internal kinematics of such objects are more representative of actual elliptical galaxies. Although our coverage of parameter space is not complete, the remnants produced in our simulations invariably possess small kinematic misalignments. These results suggest that the majority of elliptical galaxies have had evolutionary histories more complex than those envisaged by simple merger hypotheses.

Weil, Melinda L.↗