Person: Hernquist, Lars
Loading...
Email Address
AA Acceptance Date
Birth Date
Research Projects
Organizational Units
Job Title
Last Name
Hernquist
First Name
Lars
Name
Hernquist, Lars
21 results
Search Results
Now showing 1 - 10 of 21
Publication A tilted dark halo origin of the Galactic disk warp and flare(Springer Science and Business Media LLC, 2023-09-14) Han, Jiwon; Conroy, Charles; Hernquist, LarsThe outer disk of the Milky Way Galaxy is warped and flared. Several mechanisms have been proposed to explain these phenomena, but none have quantitatively reproduced both features. Recent work has demonstrated that the Galactic stellar halo is tilted with respect to the disk plane, suggesting that at least some component of the dark matter halo may also be tilted. Here we show that a dark halo tilted in the same direction as the stellar halo can induce a warp and flare in the Galactic disk at the same amplitude and orientation as the data. In our model the warp is visible in both the gas and stars of all ages, which is consistent with the breadth of observational tracers of the warp. These results, in combination with data in the stellar halo, provide compelling evidence that our Galaxy is embedded in a tilted dark matter halo.Publication First Results From the IllustrisTNG Simulations: Matter and Galaxy Clustering(Oxford University Press (OUP), 2017-12-22) Springel, Volker; Pakmor, Rüdiger; Pillepich, Annalisa; Weinberger, Rainer; Nelson, Dylan; Hernquist, Lars; Vogelsberger, Mark; Genel, Shy; Torrey, Paul; Marinacci, Federico; Naiman, JillHydrodynamical simulations of galaxy formation have now reached sufficient volume to make precision predictions for clustering on cosmologically relevant scales. Here, we use our new IllustrisTNG simulations to study the non-linear correlation functions and power spectra of baryons, dark matter, galaxies, and haloes over an exceptionally large range of scales. We find that baryonic effects increase the clustering of dark matter on small scales and damp the total matter power spectrum on scales up to k ∼ 10 h Mpc−1 by 20 per cent. The non-linear two-point correlation function of the stellar mass is close to a power-law over a wide range of scales and approximately invariant in time from very high redshift to the present. The two-point correlation function of the simulated galaxies agrees well with Sloan Digital Sky Survey at its mean redshift z ≃ 0.1, both as a function of stellar mass and when split according to galaxy colour, apart from a mild excess in the clustering of red galaxies in the stellar mass range of109–1010 h−2 M⊙. Given this agreement, the TNG simulations can make valuable theoretical predictions for the clustering bias of different galaxy samples. We find that the clustering length of the galaxy autocorrelation function depends strongly on stellar mass and redshift. Its power-law slope γ is nearly invariant with stellar mass, but declines from γ ∼ 1.8 at redshift z = 0 to γ ∼ 1.6 at redshift z ∼ 1, beyond which the slope steepens again. We detect significant scale dependences in the bias of different observational tracers of large-scale structure, extending well into the range of the baryonic acoustic oscillations and causing nominal (yet fortunately correctable) shifts of the acoustic peaks of around ∼ 5 per cent.Publication Unveiling the Role of the Magnetic Field at the Smallest Scales of Star Formation(American Astronomical Society, 2017) Hull, Charles L. H.; Mocz, Philip; Burkhart, Blakesley; Goodman, Alyssa; Girart, Josep M.; Cortés, Paulo C.; Hernquist, Lars; Springel, Volker; Li, Zhi-Yun; Lai, Shih-PingWe report Atacama Large Millimeter/submillimeter Array (ALMA) observations of polarized dust emission from the protostellar source Ser-emb 8 at a linear resolution of 140 au. Assuming models of dust-grain alignment hold, the observed polarization pattern gives a projected view of the magnetic field structure in this source. Contrary to expectations based on models of strongly magnetized star formation, the magnetic field in Ser-emb 8 does not exhibit an hourglass morphology. Combining the new ALMA data with previous observational studies, we can connect magnetic field structure from protostellar core (̃80,000 au) to disk (̃100 au) scales. We compare our observations with four magnetohydrodynamic gravo-turbulence simulations made with the AREPO code that have initial conditions ranging from super-Alfvénic (weakly magnetized) to sub-Alfvénic (strongly magnetized). These simulations achieve the spatial dynamic range necessary to resolve the collapse of protostars from the parsec scale of star-forming clouds down to the ̃100 au scale probed by ALMA. Only in the very strongly magnetized simulation do we see both the preservation of the field direction from cloud to disk scales and an hourglass-shaped field at <1000 au scales. We conduct an analysis of the relative orientation of the magnetic field and the density structure in both the Ser-emb 8 ALMA observations and the synthetic observations of the four AREPO simulations. We conclude that the Ser-emb 8 data are most similar to the weakly magnetized simulations, which exhibit random alignment, in contrast to the strongly magnetized simulation, where the magnetic field plays a role in shaping the density structure in the source. In the weak-field case, it is turbulence—not the magnetic field—that shapes the material that forms the protostar, highlighting the dominant role that turbulence can play across many orders of magnitude in spatial scale.Publication Zooming in on accretion – I. The structure of halo gas(Oxford University Press (OUP), 2016) Nelson, Dylan; Genel, Shy; Pillepich, Annalisa; Vogelsberger, Mark; Springel, Volker; Hernquist, LarsWe study the properties of gas in and around 1012 M⊙ haloes at z = 2 using a suite of high-resolution cosmological hydrodynamic ‘zoom’ simulations. We quantify the thermal and dynamical structure of these gaseous reservoirs in terms of their mean radial distributions and angular variability along different sightlines. With each halo simulated at three levels of increasing resolution, the highest reaching a baryon mass resolution of ∼10 000 solar masses, we study the interface between filamentary inflow and the quasi-static hot halo atmosphere. We highlight the discrepancy between the spatial resolution available in the halo gas as opposed to within the galaxy itself, and find that stream morphologies become increasingly complex at higher resolution, with large coherent flows revealing density and temperature structure at progressively smaller scales. Moreover, multiple gas components co-exist at the same radius within the halo, making radially averaged analyses misleading. This is particularly true where the hot, quasi-static, high entropy halo atmosphere interacts with cold, rapidly inflowing, low entropy accretion. Haloes at this mass have a well-defined virial shock, associated with a sharp jump in temperature and entropy at ≳ 1.25 rvir. The presence, radius, and radial width of this boundary feature, however, vary not only from halo to halo, but also as a function of angular direction, covering roughly ∼75 per cent of the 4π sphere. We investigate the process of gas virialization as imprinted in the halo structure, and discuss different modes for the accretion of gas from the intergalactic medium.Publication Modelling galactic conformity with the colour–halo age relation in the Illustris simulation(Oxford University Press (OUP), 2015) Bray, Aaron; Pillepich, Annalisa; Sales, Laura V.; Zhu, Emily; Genel, Shy; Rodriguez-Gomez, Vicente; Torrey, P; Nelson, Dylan; Vogelsberger, Mark; Springel, Volker; Eisenstein, Daniel; Hernquist, LarsComparisons between observational surveys and galaxy formation models find that dark matter haloes’ mass can largely explain their galaxies’ stellar mass. However, it remains uncertain whether additional environmental variables, known as assembly bias, are necessary to explain other galaxy properties. We use the Illustris simulation to investigate the role of assembly bias in producing galactic conformity by considering 18 000 galaxies with Mstellar > 2 × 109 M⊙. We find a significant signal of galactic conformity: out to distances of about 10 Mpc, the mean red fraction of galaxies around redder galaxies is higher than around bluer galaxies at fixed stellar mass. Dark matter haloes exhibit an analogous conformity signal, in which the fraction of haloes formed at earlier times (old haloes) is higher around old haloes than around younger ones at fixed halo mass. A plausible interpretation of galactic conformity is the combination of the halo conformity signal with the galaxy colour–halo age relation: at fixed stellar mass, particularly towards the low-mass end, Illustris’ galaxy colours correlate with halo age, with the reddest galaxies (often satellites) preferentially found in the oldest haloes. We explain the galactic conformity effect with a simple semi-empirical model, assigning stellar mass via halo mass (abundance matching) and galaxy colour via halo age (age matching). Regarding comparison to observations, we conclude that the adopted selection/isolation criteria, projection effects, and stacking techniques can have a significant impact on the measured amplitude of the conformity signal.Publication Recoiling black holes: prospects for detection and implications of spin alignment(Oxford University Press (OUP), 2015) Blecha, Laura; Sijacki, Debora; Kelley, Luke; Torrey, P; Vogelsberger, Mark; Nelson, Dylan; Springel, Volker; Snyder, Gregory; Hernquist, LarsSupermassive black hole (BH) mergers produce powerful gravitational wave emission. Asymmetry in this emission imparts a recoil kick to the merged BH, which can eject the BH from its host galaxy altogether. Recoiling BHs could be observed as offset active galactic nuclei (AGN). Several candidates have been identified, but systematic searches have been hampered by large uncertainties regarding their observability. By extracting merging BHs and host galaxy properties from the Illustris cosmological simulations, we have developed a comprehensive model for recoiling AGN. Here, for the first time, we model the effects of BH spin alignment and recoil dynamics based on the gas richness of host galaxies. We predict that if BH spins are not highly aligned, seeing-limited observations could resolve offset AGN, making them promising targets for all-sky surveys. For randomly oriented spins, ≲ 10 spatially offset AGN may be detectable in Hubble Space Telescope-Cosmological Evolution Survey, and >103 could be found with the Panoramic Survey Telescope & Rapid Response System (Pan-STARRS), the Large Synoptic Survey Telescope (LSST), Euclid, and the Wide-Field Infrared Survey Telescope (WFIRST). Nearly a thousand velocity offset AGN are predicted within the Sloan Digital Sky Survey (SDSS) footprint; the rarity of large broad-line offsets among SDSS quasars is likely due in part to selection effects but suggests that spin alignment plays a role in suppressing recoils. None the less, in our most physically motivated model where alignment occurs only in gas-rich mergers, hundreds of offset AGN should be found in all-sky surveys. Our findings strongly motivate a dedicated search for recoiling AGN.Publication The stellar mass assembly of galaxies in the Illustris simulation: growth by mergers and the spatial distribution of accreted stars(Oxford University Press (OUP), 2016) Rodriguez-Gomez, Vicente; Pillepich, Annalisa; Sales, Laura V.; Genel, Shy; Vogelsberger, Mark; Zhu, Qirong; Wellons, Sarah; Nelson, Dylan; Torrey, P; Springel, Volker; Ma, Chung-Pei; Hernquist, LarsWe use the Illustris simulation to study the relative contributions of in situ star formation and stellar accretion to the build-up of galaxies over an unprecedentedly wide range of masses (M* = 109-1012 M⊙), galaxy types, environments, and assembly histories. We find that the ‘two-phase’ picture of galaxy formation predicted by some models is a good approximation only for the most massive galaxies in our simulation – namely, the stellar mass growth of galaxies below a few times 1011 M⊙ is dominated by in situ star formation at all redshifts. The fraction of the total stellar mass of galaxies at z = 0 contributed by accreted stars shows a strong dependence on galaxy stellar mass, ranging from about 10 per cent for Milky Way-sized galaxies to over 80 per cent for M* ≈ 1012 M⊙ objects, yet with a large galaxy-to-galaxy variation. At a fixed stellar mass, elliptical galaxies and those formed at the centres of younger haloes exhibit larger fractions of ex situ stars than disc-like galaxies and those formed in older haloes. On average, ∼50 per cent of the ex situ stellar mass comes from major mergers (stellar mass ratio μ > 1/4), ∼20 per cent from minor mergers (1/10 < μ < 1/4), ∼20 per cent from very minor mergers (μ < 1/10), and ∼10 per cent from stars that were stripped from surviving galaxies (e.g. flybys or ongoing mergers). These components are spatially segregated, with in situ stars dominating the innermost regions of galaxies, and ex situ stars being deposited at larger galactocentric distances in order of decreasing merger mass ratio.Publication Simulated Galaxy Interactions as Probes of Merger Spectral Energy Distributions(IOP Publishing, 2014) Lanz, Lauranne; Hayward, Christopher C.; Zezas, Andreas; Smith, Howard; Ashby, Matthew; Brassington, Nicola; Fazio, Giovanni; Hernquist, LarsWe present the first systematic comparison of ultraviolet-millimeter spectral energy distributions (SEDs) of observed and simulated interacting galaxies. Our sample is drawn from the Spitzer Interacting Galaxy Survey and probes a range of galaxy interaction parameters. We use 31 galaxies in 14 systems which have been observed with Herschel, Spitzer, GALEX, and 2MASS. We create a suite of GADGET-3 hydrodynamic simulations of isolated and interacting galaxies with stellar masses comparable to those in our sample of interacting galaxies. Photometry for the simulated systems is then calculated with the SUNRISE radiative transfer code for comparison with the observed systems. For most of the observed systems, one or more of the simulated SEDs match reasonably well. The best matches recover the infrared luminosity and the star formation rate of the observed systems, and the more massive systems preferentially match SEDs from simulations of more massive galaxies. The most morphologically distorted systems in our sample are best matched to the simulated SEDs that are close to coalescence, while less evolved systems match well with the SEDs over a wide range of interaction stages, suggesting that an SED alone is insufficient for identifying the interaction stage except during the most active phases in strongly interacting systems. This result is supported by our finding that the SEDs calculated for simulated systems vary little over the interaction sequence.Publication Stellar orbit evolution in close circumstellar disc encounters(Oxford University Press (OUP), 2014) Munoz, Diego Jose; Kratter, K.; Vogelsberger, M.; Hernquist, Lars; Springel, V.The formation and early evolution of circumstellar discs often occurs within dense, newborn stellar clusters. For the first time, we apply the moving-mesh code arepo, to circumstellar discs in 3D, focusing on disc–disc interactions that result from stellar flybys. Although a small fraction of stars are expected to undergo close approaches, the outcomes of the most violent encounters might leave an imprint on the discs and host stars that will influence both their orbits and their ability to form planets. We first construct well-behaved 3D models of self-gravitating discs, and then create a suite of numerical experiments of parabolic encounters, exploring the effects of pericentre separation rp, disc orientation and disc–star mass ratio (Md/M*) on the orbital evolution of the host stars. Close encounters (2rp ≲ disc radius) can truncate discs on very short time-scales. If discs are massive, close encounters facilitate enough orbital angular momentum extraction to induce stellar capture. We find that for realistic primordial disc masses Md ≲ 0.1M*, non-colliding encounters induce minor orbital changes, which is consistent with analytic calculations of encounters in the linear regime. The same disc masses produce entirely different results for grazing/colliding encounters. In the latter case, rapidly cooling discs lose orbital energy by radiating away the energy excess of the shock-heated gas, thus causing capture of the host stars into a bound orbit. In rare cases, a tight binary with a circumbinary disc forms as a result of this encounter.Publication A moving mesh unstaggered constrained transport scheme for magnetohydrodynamics(Oxford University Press (OUP), 2016) Mocz, Philip; Pakmor, Rüdiger; Springel, Volker; Vogelsberger, Mark; Marinacci, Federico; Hernquist, LarsWe present a constrained transport (CT) algorithm for solving the 3D ideal magnetohydrodynamic (MHD) equations on a moving mesh, which maintains the divergence-free condition on the magnetic field to machine-precision. Our CT scheme uses an unstructured representation of the magnetic vector potential, making the numerical method simple and computationally efficient. The scheme is implemented in the moving mesh code arepo. We demonstrate the performance of the approach with simulations of driven MHD turbulence, a magnetized disc galaxy, and a cosmological volume with primordial magnetic field. We compare the outcomes of these experiments to those obtained with a previously implemented Powell divergence-cleaning scheme. While CT and the Powell technique yield similar results in idealized test problems, some differences are seen in situations more representative of astrophysical flows. In the turbulence simulations, the Powell cleaning scheme artificially grows the mean magnetic field, while CT maintains this conserved quantity of ideal MHD. In the disc simulation, CT gives slower magnetic field growth rate and saturates to equipartition between the turbulent kinetic energy and magnetic energy, whereas Powell cleaning produces a dynamically dominant magnetic field. Such difference has been observed in adaptive-mesh refinement codes with CT and smoothed-particle hydrodynamics codes with divergence-cleaning. In the cosmological simulation, both approaches give similar magnetic amplification, but Powell exhibits more cell-level noise. CT methods in general are more accurate than divergence-cleaning techniques, and, when coupled to a moving mesh can exploit the advantages of automatic spatial/temporal adaptivity and reduced advection errors, allowing for improved astrophysical MHD simulations.
- «
- 1 (current)
- 2
- 3
- »