Source author record

Lars Hernquist

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

215works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

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

215 published item(s)

preprint2025arXiv

The THESAN-ZOOM project: Star-formation efficiencies in high-redshift galaxies

Recent JWST observations hint at unexpectedly intense cosmic star-formation in the early Universe, often attributed to enhanced star-formation efficiencies (SFEs). Here, we analyze the SFE in THESAN-ZOOM, a novel zoom-in radiation-hydrodynamic simulation campaign of high-redshift ($z \gtrsim 3$) galaxies employing a state-of-the-art galaxy formation model resolving the multiphase interstellar medium (ISM). The halo-scale SFE ($ε^{\ast}_{\rm halo}$) - the fraction of baryons accreted by a halo that are converted to stars - follows a double power-law dependence on halo mass, with a mild redshift evolution above $M_{\rm halo} \gtrsim 10^{9.5}\,{\rm M}_{\odot}$. The power-law slope is roughly $1/3$ at large halo masses, consistent with expectations when gas outflows are momentum-driven. At lower masses, the slope is roughly $2/3$ and is more aligned with the energy-driven outflow scenario. $ε^{\ast}_{\rm halo}$ is a factor of $2-3$ larger than commonly assumed in empirical galaxy-formation models at $M_{\rm halo} \lesssim 10^{11}\,{\rm M}_{\odot}$. On galactic (kpc) scales, the Kennicutt-Schmidt (KS) relation of neutral gas is universal in THESAN-ZOOM, following $Σ_{\rm SFR} \propto Σ_{\rm gas}^2$, indicative of a turbulent energy balance in the ISM maintained by stellar feedback. The rise of $ε^{\ast}_{\rm halo}$ with halo mass can be traced primarily to increasing gas surface densities in massive galaxies, while the underlying KS relation and neutral, star-forming gas fraction remain unchanged. Although the increase in $ε^{\ast}_{\rm halo}$ with redshift is relatively modest, it is sufficient to explain the large observed number density of UV-bright galaxies at $z \gtrsim 12$. However, reproducing the brightest sources at $M_{\rm UV} \lesssim -21$ may require extrapolating the SFE beyond the halo mass range directly covered by THESAN-ZOOM.

preprint2024arXiv

Gas Morphology of Milky Way-like Galaxies in the TNG50 Simulation: Signals of Twisting and Stretching

We present an in-depth analysis of gas morphologies for a sample of 25 Milky Way-like galaxies from the IllustrisTNG TNG50 simulation. We constrain the morphology of cold, warm, hot gas, and gas particles as a whole using a Local Shell Iterative Method (LSIM) and explore its observational implications by computing the hard-to-soft X-ray ratio, which ranges between $10^{-3}$-$10^{-2}$ in the inner $\sim 50 \rm kpc$ of the distribution and $10^{-5}$-$10^{-4}$ at the outer portion of the hot gas distribution. We group galaxies into three main categories: simple, stretched, and twisted. These categories are based on the radial reorientation of the principal axes of the reduced inertia tensor. We find that a vast majority ($77\%$) of the galaxies in our sample exhibit twisting patterns in their radial profiles. Additionally, we present detailed comparisons between 1) the gaseous distributions belonging to individual temperature regimes, 2) the cold gas distributions and stellar distributions, and 3) the gaseous distributions and dark matter (DM) halos. We find a strong correlation between the morphological properties of the cold gas and stellar distributions. Furthermore, we find a correlation between gaseous distributions with DM halo that increases with gas temperature, implying that we may use the warm-hot gaseous morphology as a tracer to probe the DM morphology. Finally, we show gaseous distributions exhibit significantly more prolate morphologies than the stellar distributions and DM halos, which we hypothesize is due to stellar and AGN feedback.

preprint2022arXiv

Breaking baryon-cosmology degeneracy with the electron density power spectrum

Uncertain feedback processes in galaxies affect the distribution of matter, currently limiting the power of weak lensing surveys. If we can identify cosmological statistics that are robust against these uncertainties, or constrain these effects by other means, then we can enhance the power of current and upcoming observations from weak lensing surveys such as DES, Euclid, the Rubin Observatory, and the Roman Space Telescope. In this work, we investigate the potential of the electron density auto-power spectrum as a robust probe of cosmology and baryonic feedback. We use a suite of (magneto-)hydrodynamic simulations from the CAMELS project and perform an idealized analysis to forecast statistical uncertainties on a limited set of cosmological and physically-motivated astrophysical parameters. We find that the electron number density auto-correlation, measurable through either kinematic Sunyaev-Zel'dovich observations or through Fast Radio Burst dispersion measures, provides tight constraints on $Ω_{m}$ and the mean baryon fraction in intermediate-mass halos, $\bar{f}_{\mathrm{bar}}$. By obtaining an empirical measure for the associated systematic uncertainties, we find these constraints to be largely robust to differences in baryonic feedback models implemented in hydrodynamic simulations. We further discuss the main caveats associated with our analysis, and point out possible directions for future work.

preprint2022arXiv

Dynamics of Intermediate-Mass Black Holes Wandering in the Milky Way Galaxy Using the Illustris TNG50 Simulation

The detection of Intermediate-Mass Black Holes (IMBHs) in dwarf galaxies is crucial to closing the gap in the wide mass distribution of black holes ($\sim 3 \, \rm M_{\odot}$ to $\sim 5 \times 10^{10} \, \rm M_{\odot}$). IMBHs originally located at the center of dwarfs that later collide with the Milky Way (MW) could be wandering, undetected, in our Galaxy. We used TNG50, the highest-resolution run of the IllustrisTNG project, to study the kinematics and dynamics of star clusters, in the appropriate mass range, acting as IMBH proxies in a MW analog galaxy. We showed that $\sim 87\%$ of our studied IMBHs drift inward. The radial velocity of these sinking IMBHs has a median magnitude of $\sim 0.44 \, \rm ckpc \, h^{-1} \, Gyr^{-1}$ and no dependence on the black hole mass. The central $1 \, \rm ckpc \, h^{-1}$ has the highest number density of IMBHs in the galaxy. A physical toy model with linear drag forces was developed to explain the orbital circularization with time. These findings constrain the spatial distribution of IMBHs, suggesting that future searches should focus on the central regions of the Galaxy. Additionally, we found that the 3D velocity distribution of IMBHs with respect to the galactic center has a mean of $\sim 180 \, \rm km \, s^{-1}$ and larger variance with decreasing radius. Remarkably, the velocity distribution relative to the local gas shows significantly lower values, with a mean of $\sim 88 \, \rm km \, s^{-1}$. These results are instrumental for predicting the accretion and radiation properties of IMBHs, facilitating their detection with future surveys.

preprint2022arXiv

Early-type galaxy density profiles from IllustrisTNG: III. Effects on outer kinematic structure

Early-type galaxies (ETGs) possess total density profiles close to isothermal, which can lead to non-Gaussian line-of-sight velocity dispersion (LOSVD) under anisotropic stellar orbits. However, recent observations of local ETGs in the MASSIVE Survey reveal outer kinematic structures at $1.5 R_{\mathrm{eff}}$ (effective radius) that are inconsistent with fixed isothermal density profiles; the authors proposed varying density profiles as an explanation. We aim to verify this conjecture and understand the influence of stellar assembly on these kinematic features through mock ETGs in IllustrisTNG. We create mock Integral-Field-Unit observations to extract projected stellar kinematic features for 207 ETGs with stellar mass $M_{\ast}\geqslant 10^{11} \mathrm{M_{\odot}}$ in TNG100-1. The mock observations reproduce the key outer ($1.5R_{\mathrm{eff}}$) kinematic structures in the MASSIVE ETGs, including the puzzling positive correlation between velocity dispersion profile outer slope $γ_{\mathrm{outer}}$ and the kurtosis $h_{4}$'s gradient. We find that $h_{4}$ is uncorrelated with stellar orbital anisotropy beyond $R_{\mathrm{eff}}$; instead we find that the variations in $γ_{\mathrm{outer}}$ and outer $h_{4}$ (a good proxy for $h_{4}$ gradient) are both driven by variations of the density profile at the outskirts across different ETGs. These findings corroborate the proposed conjecture and rule out velocity anisotropy as the origin of non-Gaussian outer kinematic structure in ETGs. We also find that the outer kurtosis and anisotropy correlate with different stellar assembly components, with the former related to minor mergers or flyby interactions while the latter is mainly driven by major mergers, suggesting distinct stellar assembly origins that decorrelates the two quantities.

preprint2022arXiv

Fast, Slow, Early, Late: Quenching Massive Galaxies at z~0.8

We investigate the stellar populations for a sample of 161 massive, mainly quiescent galaxies at $\langle z_{\rm obs} \rangle=0.8$ with deep Keck/DEIMOS rest-frame optical spectroscopy (HALO7D survey). With the fully Bayesian framework Prospector, we simultaneously fit the spectroscopic and photometric data with an advanced physical model (including non-parametric star-formation histories, emission lines, variable dust attenuation law, and dust and AGN emission) together with an uncertainty and outlier model. We show that both spectroscopy and photometry are needed to break the dust-age-metallicity degeneracy. We find a large diversity of star-formation histories: although the most massive ($M_{\star}>2\times10^{11}~M_{\odot}$) galaxies formed the earliest (formation redshift of $z_{\rm f}\approx5-10$ with a short star-formation timescale of $τ_{\rm SF}\lesssim1~\mathrm{Gyr}$), lower-mass galaxies have a wide range of formation redshifts, leading to only a weak trend of $z_{\rm f}$ with $M_{\star}$. Interestingly, several low-mass galaxies with have formation redshifts of $z_{\rm f}\approx5-8$. Star-forming galaxies evolve about the star-forming main sequence, crossing the ridgeline several times in their past. Quiescent galaxies show a wide range and continuous distribution of quenching timescales ($τ_{\rm quench}\approx0-5~\mathrm{Gyr}$) with a median of $\langleτ_{\rm quench}\rangle=1.0_{-0.9}^{+0.8}~\mathrm{Gyr}$ and of quenching epochs of $z_{\rm quench}\approx0.8-5.0$ ($\langle z_{\rm quench}\rangle=1.3_{-0.4}^{+0.7}$). This large diversity of quenching timescales and epochs points toward a combination of internal and external quenching mechanisms. In our sample, rejuvenation and "late bloomers" are uncommon. In summary, our analysis supports the "grow & quench" framework and is consistent with a wide and continuously-populated diversity of quenching timescales.

preprint2022arXiv

Formation and fate of low metallicity stars in IllustrisTNG50

Low metallicity stars give rise to unique spectacular transients and are of immense interest for understanding stellar evolution. Their importance has only grown further with the recent detections of mergers of stellar mass black holes that likely originate mainly from low metallicity progenitor systems. Moreover, the formation of low metallicity stars is intricately linked to galaxy evolution, in particular to early enrichment and to later accretion and mixing of lower metallicity gas. Because low metallicity stars are difficult to observe directly, cosmological simulations are crucial for understanding their formation. Here we quantify the rates and locations of low metallicity star formation using the high-resolution TNG50 magnetohydrodynamical cosmological simulation, and we examine where low metallicity stars end up at $z=0$. We find that $20\%$ of stars with $Z_*<0.1\,\mathrm{Z_\odot}$ form after $z=2$, and that such stars are still forming in galaxies of all masses at $z=0$ today. Moreover, most low-metallicity stars at $z=0$ reside in massive galaxies. We analyse the radial distribution of low metallicity star formation, and discuss the curious case of seven galaxies in TNG50 that form stars from primordial gas even at $z=0$.

preprint2022arXiv

Galactic angular momentum in the IllustrisTNG simulation -- I. Connection to morphology, halo spin, and black hole mass

We use the TNG100 simulation of the IllustrisTNG project to investigate the stellar specific angular momenta ($j_{\ast}$) of $\sim$12,000 central galaxies at $z=0$ in a full cosmological context, with stellar masses ($M_{\ast}$) ranging from $10^{9}$ to $10^{12} \, {\rm M}_{\odot}$. We find that the $j_{\ast}$-$M_{\ast}$ relations for early-type and late-type galaxies in IllustrisTNG are in good overall agreement with observations, and that these galaxy types typically `retain' $\sim$10-20 and $\sim$50-60 per cent of their host haloes' specific angular momenta, respectively, with some dependence on the methodology used to measure galaxy morphology. We present results for kinematic as well as visual-like morphological measurements of the simulated galaxies. Next, we explore the scatter in the $j_{\ast}$-$M_{\ast}$ relation with respect to the spin of the dark matter halo and the mass of the supermassive black hole (BH) at the galactic centre. We find that galaxies residing in faster spinning haloes, as well as those hosting less massive BHs, tend to have a higher specific angular momentum. We also find that, at fixed galaxy or halo mass, halo spin and BH mass are anticorrelated with each other, probably as a consequence of more efficient gas flow toward the galactic centre in slowly rotating systems. Finally, we show that halo spin plays an important role in determining galaxy sizes - larger discs form at the centres of faster-rotating haloes - although the trend breaks down for massive galaxies with $M_{\ast} \gtrsim 10^{11} \, {\rm M}_{\odot}$, roughly the mass scale at which a galaxy's stellar mass becomes dominated by accreted stars.

preprint2022arXiv

Galaxy Formation in the Santa Cruz semi-analytic model compared with IllustrisTNG -- I. Galaxy scaling relations, dispersions, and residuals at z=0

We present the first results from applying the Santa Cruz semi-analytic model (SAM) for galaxy formation on merger trees extracted from a dark matter only version of the IllustrisTNG (TNG) simulations. We carry out a statistical comparison between the predictions of the Santa Cruz SAM and TNG for a subset of central galaxy properties at z = 0, with a focus on stellar mass, cold and hot gas mass, star formation rate (SFR), and black hole (BH) mass. We find fairly good agreement between the mean predictions of the two methods for stellar mass functions and the stellar mass vs. halo mass (SMHM) relation, and qualitatively good agreement between the SFR or cold gas mass vs. stellar mass relation and quenched fraction as a function of stellar mass. There are greater differences between the predictions for hot (circumgalactic) gas mass and BH mass as a function of halo mass. Going beyond the mean relations, we also compare the dispersion in the predicted scaling relations, and the correlation in residuals on a halo-by-halo basis between halo mass and galaxy property scaling relations. Intriguingly, we find similar correlations between residuals in SMHM in the SAM and in TNG, suggesting that these relations may be shaped by similar physical processes. Other scaling relations do not show significant correlations in the residuals, indicating that the physics implementations in the SAM and TNG are significantly different.

preprint2022arXiv

Gas-phase metallicity break radii of star-forming galaxies in IllustrisTNG

We present radial gas-phase metallicity profiles, gradients, and break radii at redshift $z = 0 - 3$ from the TNG50-1 star-forming galaxy population. These metallicity profiles are characterized by an emphasis on identifying the steep inner gradient and flat outer gradient. From this, the break radius, $r_{\rm Break}$, is defined as the region where the transition occurs. We observe the break radius having a positive trend with mass that weakens with redshift. When normalized by the stellar half-mass radius, the break radius has a weaker relation with both mass and redshift. To test if our results are dependent on the resolution or adopted physics of TNG50-1, the same analysis is performed in TNG50-2 and Illustris-1. We find general agreement between each of the simulations in their qualitative trends; however, the adopted physics between TNG and Illustris differ and therefore the breaks, normalized by galaxy size, deviate by a factor of $\sim$2. In order to understand where the break comes from, we define two relevant time-scales: an enrichment time-scale and a radial gas mixing time-scale. We find that $r_{\rm Break}$ occurs where the gas mixing time-scale is $\sim$10 times as long as the enrichment time-scale in all three simulation runs, with some weak mass and redshift dependence. This implies that galactic disks can be thought of in two-parts: a star-forming inner disk with a steep gradient and a mixing-dominated outer disk with a flat gradient, with the break radius marking the region of transition between them.

preprint2022arXiv

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

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 $f_{\rm col}\sim5-10\%$, only weakly dependent on radius and vertical height, and that scattering boosts the H$α$ luminosity by $\sim40\%$. The dust correction via the Balmer decrement works well (intrinsic H$α$ emission recoverable within $25\%$), 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}\approx28\%$ and $f_{\rm He}\approx9\%$, respectively. Together with an escape fraction of $f_{\rm esc}\approx6\%$, this reduces the available budget for hydrogen line emission by nearly half ($f_{\rm H}\approx57\%$). 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.

preprint2022arXiv

Mass of the dynamically hot inner stellar halo predicts the ancient accreted stellar mass

Galactic dynamical structures are fossil records of the assembly histories of galaxies. By analyzing the cosmological hydrodynamical simulation TNG50, we find that a dynamical structure that we call the "hot inner stellar halo," defined by stars on dynamically hot orbits with circularity $λ_z < 0.5$ at $3.5\,{\rm kpc}<r \lesssim 2\,R_e$, is a strong indicator of the mass of accreted satellite galaxies. We find a strong correlation between the mass of this hot inner stellar halo and the total ex situ stellar mass. There is a similarly strong correlation with the stellar mass of the most massive secondary galaxy ever merged. These TNG50 correlations are compatible with those predicted by other simulations, for example by TNG100 across the whole mass range under study (galaxy stellar masses, $M_*$, in the $10^{10.3-11.6}$\,\Msun\, range) and by EAGLE for $M_* \gtrsim 10^{10.6} $\,\Msun\, galaxies.\ This shows that our predictions are robust across different galaxy formation and feedback models and hold across a wide range of numerical resolution. The hot inner stellar halo is a product of massive and typically ancient mergers, with inner-halo stars exhibiting three main physical origins: accreted and stripped from massive satellites, dynamically heated by mergers from the bulge and/or disk in the main progenitor, and formed from star formation triggered during mergers. The mass of the hot inner stellar halo defined in this paper is a quantity that can be robustly obtained for real galaxies by applying a population-orbit superposition method to integral-field-unit spectroscopy data, out to a distance of $\sim2\,R_e$, which is possible with current observations. Hence, this paper shows that integral-field-unit observations and dynamical models of the inner regions of galaxies provide a way to quantitatively determine the mass of ancient accreted satellites.

preprint2022arXiv

Morphological decomposition of TNG50 galaxies: methodology and catalogue

We present MORDOR (MORphological DecOmposeR, a new algorithm for structural decomposition of simulated galaxies based on stellar kinematics. The code measures the properties of up to five structural components (a thin/cold and a thick/warm disc, a classical and a secular bulge, and a spherical stellar halo), and determines the properties of a stellar bar (if present). A comparison with other algorithms presented in the literature yields overall good agreement, with MORDOR displaying a higher flexibility in correctly decomposing systems and identifying bars in crowded environments (e.g. with ongoing fly-bys, often observable in cosmological simulations). We use MORDOR to analyse galaxies in the TNG50 simulation and find the following: ($i$) the thick disc component undergoes the strongest evolution in the binding energy-circularity plane, as expected when disc galaxies decrease their turbulent-rotational support with cosmic time; ($ii$) smaller galaxies (with stellar mass, $10^{9} \lesssim M_{*} / {\rm M_{\odot}} \leq 5 \times 10^{9}$) undergo a major growth in their disc components after $z\sim 1$, whereas ($iii$) the most massive galaxies ($5 \times 10^{10} < M_{*} / {\rm M_{\odot}} \leq 5\times10^{11}$) evolve toward more spheroidal dominated objects down to $z=0$ due to frequent gravitational interactions with satellites; ($iv$) the fraction of barred galaxies grows rapidly at high redshift and stabilizes below $z\sim 2$, except for the most massive galaxies that show a decrease in the bar occupation fraction at low redshift; ($v$) galaxies with $M_{*} \sim 10^{11}~{\rm M_{\odot}}$ exhibit the highest relative occurrence of bars at $z=0$, in agreement with observational studies. We publicly release MORDOR and the morphological catalogue of TNG50 galaxies.

preprint2022arXiv

Percent-level constraints on baryonic feedback with spectral distortion measurements

High-significance measurements of the monopole thermal Sunyaev-Zel'dovich CMB spectral distortions have the potential to tightly constrain poorly understood baryonic feedback processes. The sky-averaged Compton-y distortion and its relativistic correction are measures of the total thermal energy in electrons in the observable universe and their mean temperature. We use the CAMELS suite of hydrodynamic simulations to explore possible constraints on parameters describing the subgrid implementation of feedback from active galactic nuclei and supernovae, assuming a PIXIE-like measurement. The small 25 Mpc/h CAMELS boxes present challenges due to the significant cosmic variance. We utilize machine learning to construct interpolators through the noisy simulation data. Using the halo model, we translate the simulation halo mass functions into correction factors to reduce cosmic variance where required. Our results depend on the subgrid model. In the case of IllustrisTNG, we find that the best-determined parameter combination can be measured to ~2% and corresponds to a product of AGN and SN feedback. In the case of SIMBA, the tightest constraint is ~0.2% on a ratio between AGN and SN feedback. A second orthogonal parameter combination can be measured to ~8%. Our results demonstrate the significant constraining power a measurement of the late-time spectral distortion monopoles would have for baryonic feedback models.

preprint2022arXiv

Probing the $z\gtrsim6$ quasars in a universe with IllustrisTNG physics: Impact of gas-based black hole seeding models

We explore implications of a range of black hole (BH) seeding prescriptions on the formation of the brightest $z\gtrsim6$ quasars in cosmological hydrodynamic simulations. The underlying galaxy formation model is the same as in IllustrisTNG. Using constrained initial conditions, we study the growth of BHs in rare overdense regions (forming $\gtrsim10^{12}M_{\odot}/h$ halos by $z=7$) using a $(9~\mathrm{Mpc}/h)^3$ simulated volume. BH growth is maximal within halos that are compact and have a low tidal field. For these halos, we consider an array of gas-based seeding prescriptions wherein $M_{\mathrm{seed}}=10^4-10^6~M_{\odot}/h$ seeds are inserted in halos above critical thresholds for halo mass and dense, metal-poor gas mass (defined as $\tilde{M}_{\mathrm{h}}$ and $\tilde{M}_{\mathrm{sf,mp}}$, respectively, in units of $M_{\mathrm{seed}}$). We find that a seed model with $\tilde{M}_{\mathrm{sf,mp}}=5$ and $\tilde{M}_{\mathrm{h}}=3000$ successfully produces a $z\sim6$ quasar with $\sim10^9~M_{\odot}$ mass and $\sim10^{47}~\mathrm{ergs~s^ {-1}}$ luminosity. BH mergers play a crucial role at $z\gtrsim9$, causing an early boost in BH mass at a time when accretion-driven BH growth is negligible. When more stringent seeding conditions are applied (for e.g., $\tilde{M}_{\mathrm{sf,mp}}=1000$), the relative paucity of BH seeds results in a much lower merger rate. In this case, $z\gtrsim6$ quasars can only be formed if we enhance the maximum allowed BH accretion rates (by factors $\gtrsim10$) compared to the accretion model used in IllustrisTNG. This can be achieved either by allowing for super-Eddington accretion, or by reducing the radiative efficiency. Our results show that progenitors of $z\sim6$ quasars have distinct BH merger histories for different seeding models, which will be distinguishable with LISA observations.

preprint2022arXiv

Sunyaev-Zel'dovich effect and X-ray scaling relations of galaxies, groups and clusters in the IllustrisTNG simulations

Observable thermodynamical properties of the intracluster medium (ICM) reflect the complex interplay between AGN feedback and the gravitational collapse of haloes. Using the large volume TNG300 simulation of the IllustrisTNG project we provide predictions for X-ray and Sunyaev-Zel'dovich (SZ) scaling relations for a sample of over 30,000 haloes that cover a wide mass range from galaxies to massive galaxy clusters ($M_{\rm 500crit}$ $\in [10^{12}$ M$_{\odot} - 2\times 10^{15}$ M$_{\odot}$]). We produce mock X-ray observations of simulated haloes using methods that are consistent with observational techniques. Thus, we investigate the scaling relations between the soft-band X-ray luminosity, spectroscopic temperature, gas mass fraction, $Y_{\rm X}$ and $Y_{\rm SZ}$ as a function of halo mass, and we find broad agreement between IllustrisTNG and the observed relations. Our results highlight the scatter and bias introduced by estimated masses, and thus the importance of converting simulated ICM properties to the observable space when comparing simulations to current X-ray observations. The wide range of halo masses in our sample provides new insights into the shape of the X-ray and SZ scaling relations across three orders of magnitude in mass. Our findings show strong evidence for a break in $z=0$ scaling relations. We introduce a smoothly broken power law model which robustly captures the location of this break, the width of the transition region around the break, as well as the slope dependence on halo mass. Our results inform the next generation of subgrid black hole feedback models and provide predictions for ongoing and future observational surveys.

preprint2022arXiv

The buildup of strongly barred galaxies in the TNG100 simulation

We analyse the properties of strongly barred disc galaxies using the TNG100 simulation, a cosmological hydrodynamical realisation of the IllustrisTNG suite. We identify 270 disc galaxies at $z=0$ in the stellar mass range $M_{*}=10^{10.4-11}M_{\odot}$, of which 40 per cent are barred. Of the detected bars, more than half are strong. We find that the fraction of barred galaxies increases with stellar mass, in agreement with observational results. Strongly barred galaxies exhibit, overall, lower gas-to-stellar mass ratios compared to unbarred galaxies. The majority of barred galaxies are quenched (sSFR $\sim10^{-11.7} $yr$^{-1}$), whereas unbarred galaxies continue to be active (sSFR $\sim10^{-10.3}$yr$^{-1}$) on the main sequence of star-forming galaxies. We explore the evolution of strongly barred and unbarred galaxies to investigate their formation and quenching histories. We find that strong bars form between $0.5< z< 1.5$, with more massive galaxies hosting older bars. Strong bars form in galaxies with an early-established prominent disc component, undergoing periods of enhanced star formation and black hole accretion, possibly assisted by cosmological inflows. Unbarred galaxies, on the other hand, assemble most of their mass and disc component at late times. The nuclear region of strongly barred galaxies quenches shortly after bar formation, while unbarred galaxies remain active across time. Our findings are indicative of bar quenching, possibly assisted by nuclear feedback processes. We conclude that the cosmological environment, together with small scale feedback processes, determine the chances of a galaxy to form a bar and to rapidly quench its central region.

preprint2022arXiv

The Circumgalactic Medium from the CAMELS Simulations: Forecasting Constraints on Feedback Processes from Future Sunyaev-Zeldovich Observations

The cycle of baryons through the circumgalactic medium (CGM) is important to understand in the context of galaxy formation and evolution. In this study we forecast constraints on the feedback processes heating the CGM with current and future Sunyaev-Zeldovich (SZ) observations. To constrain these processes, we use a suite of cosmological simulations, the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS), that varies four different feedback parameters of two previously existing hydrodynamical simulations, IllustrisTNG and SIMBA. We capture the dependencies of SZ radial profiles on these feedback parameters with an emulator, calculate their derivatives, and forecast future constraints on these feedback parameters from upcoming experiments. We find that for a DESI-like (Dark Energy Spectroscopic Instrument) galaxy sample observed by the Simons Observatory all four feedback parameters are able to be constrained (some within the $10\%$ level), indicating that future observations will be able to further restrict the parameter space for these sub-grid models. Given the modeled galaxy sample and forecasted errors in this work, we find that the inner SZ profiles contribute more to the constraining power than the outer profiles. Finally, we find that, despite the wide range of AGN feedback parameter variation in the CAMELS simulation suite, we cannot reproduce the tSZ signal of galaxies selected by the Baryon Oscillation Spectroscopic Survey as measured by the Atacama Cosmology Telescope.

preprint2022arXiv

The distinct stellar-to-halo mass relations of satellite and central galaxies: insights from the IllustrisTNG simulations

We study the stellar-to-halo mass relation (SHMR) for central and satellite galaxies with total dynamical masses above 10^10.5 Msun using the suite of cosmological magneto-hydrodynamical simulations IllustrisTNG. In particular, we quantify environmental effects on satellite populations from TNG50, TNG100, and TNG300 located within the virial radius of group- and cluster-like hosts with total masses of 10^12-15.2 Msun. At fixed stellar mass, the satellite SHMR exhibits a distinct shift towards lower dynamical mass compared to the SHMR of centrals. Conversely, at fixed dynamical mass, satellite galaxies appear to have larger stellar-to-total mass fractions than centrals by up to a factor of a few. The systematic deviation from the central SHMR is larger for satellites in more massive hosts, at smaller cluster-centric distances, with earlier infall times, and that inhabit higher local density environments; moreover, it is in place already at early times (z < 2). Systematic environmental effects might contribute to the perceived galaxy-to-galaxy variation in the measured SHMR when galaxies cannot be separated into satellites and centrals. The SHMR of satellites exhibits a larger scatter than centrals, over the whole range of dynamical mass (by up to 0.8 dex). The shift of the satellite SHMR results mostly from tidal stripping of their dark matter, which affects satellites in an outside-in fashion: the departure of the satellite SHMR from the centrals' relation diminishes for measurements of dynamical mass in progressively smaller apertures. Finally, we provide a family of fitting functions for the SHMR predicted by IllustrisTNG.

preprint2022arXiv

The evolution of the barred galaxy population in the TNG50 simulation

We use the magnetic-hydrodynamical simulation TNG50 to study the evolution of barred massive disc galaxies. Massive spiral galaxies are already present as early as $z=4$, and bar formation takes place already at those early times. The bars grow longer and stronger as the host galaxies evolve, with the bar sizes increasing at a pace similar to that of the disc scale lengths. The bar fraction mildly evolves with redshift for galaxies with $M_{*}\geq10^{10}M\odot$, being greater than $\sim40\%$ at $0.5<z<3$ and $\sim30\%$ at $z=0$. When bars larger than a given physical size ($\geq 2\,\rm kpc$) or the angular resolution limit of twice the I-band angular PSF FWHM of the HST are considered, the bar fraction dramatically decreases with increasing redshift, reconciling the theoretical predictions with observational data. We find that barred galaxies have an older stellar population, lower gas fractions and star formation rates than unbarred galaxies. In most cases, the discs of barred galaxies assembled earlier and faster than the discs of unbarred galaxies. We also find that barred galaxies are typical in haloes with larger concentrations and smaller spin parameters than unbarred galaxies. Furthermore, the inner regions of barred galaxies are more baryon-dominated than those of unbarred galaxies but have comparable global stellar mass fractions. Our findings suggest that the bar population could be used as a potential tracer of the buildup of disc galaxies and their host haloes. With this paper, we release a catalogue of barred galaxies in TNG50 at $6$ redshifts between $z=4$ and $z=0$.

preprint2022arXiv

The formation of low surface brightness galaxies in the IllustrisTNG simulation

We explore the nature of low surface brightness galaxies (LSBGs) in the hydrodynamic cosmological simulation TNG100 of the IllustrisTNG project, selecting a sample of LSBGs ($r$-band effective surface brightness $μ_r > 22.0$ mag arcsec$^{-2}$) at $z=0$ over a wide range of stellar masses ($M_{\ast} = 10^{9}$-$10^{12}$ M$_\odot$). We find LSBGs of all stellar masses, although they are particularly prevalent at $M_{\ast} < 10^{10}$ M$_\odot$. We show that the specific star formation rates of LSBGs are not significantly different from those of high surface brightness galaxies (HSBGs) but, as a population, LSBGs are systematically less massive and more extended than HSBGs, and tend to display late-type morphologies according to a kinematic criterion. At fixed stellar mass, we find that haloes hosting LSBGs are systematically more massive and have a higher baryonic fraction than those hosting HSBGs. We find that LSBGs have higher stellar specific angular momentum and halo spin parameter values compared to HSBGs, as suggested by previous works. We track the evolution of these quantities back in time, finding that the spin parameters of the haloes hosting LSBGs and HSBGs exhibit a clear bifurcation at $z \sim 2$, which causes a similar separation in the evolutionary tracks of other properties such as galactic angular momentum and effective radius, ultimately resulting in the values observed at $z =$ 0. The higher values of specific stellar angular momentum and halo spin in LSBGs seem to be responsible for their extended nature, preventing material from collapsing into the central regions of the galaxies, also causing LSBGs to host less massive black holes at their centres.

preprint2022arXiv

The observability of galaxy merger signatures in nearby gas-rich spirals

Galaxy mergers are crucial to understanding galaxy evolution, therefore we must determine their observational signatures to select them from large IFU galaxy samples such as MUSE and SAMI. We employ 24 high-resolution idealised hydrodynamical galaxy merger simulations based on the "Feedback In Realistic Environment" (FIRE-2) model to determine the observability of mergers to various configurations and stages using synthetic images and velocity maps. Our mergers cover a range of orbital configurations at fixed 1:2.5 stellar mass ratio for two gas rich spirals at low redshift. Morphological and kinematic asymmetries are computed for synthetic images and velocity maps spanning each interaction. We divide the interaction sequence into three: (1) the pair phase; (2) the merging phase; and (3) the post-coalescence phase. We correctly identify mergers between first pericentre passage and 500 Myr after coalescence using kinematic asymmetry with 66% completeness, depending upon merger phase and the field-of-view of the observation. We detect fewer mergers in the pair phase (40%) and many more in the merging and post-coalescence phases (97%). We find that merger detectability decreases with field-of-view, except in retrograde mergers, where centrally concentrated asymmetric kinematic features enhances their detectability. Using a cut-off derived from a combination of photometric and kinematic asymmetry, we increase these detections to 89% overall, 79% in pairs, and close to 100% in the merging and post-coalescent phases. By using this combined asymmetry cut-off we mitigate some of the effects caused by smaller fields-of-view subtended by massively multiplexed integral field spectroscopy programmes.

preprint2022arXiv

The THESAN project: predictions for multi-tracer line intensity mapping in the epoch of reionization

Line intensity mapping (LIM) is rapidly emerging as a powerful technique to study galaxy formation and cosmology in the high-redshift Universe. We present LIM estimates of select spectral lines originating from the interstellar medium (ISM) of galaxies and 21 cm emission from neutral hydrogen gas in the Universe using the large volume, high resolution THESAN reionization simulations. A combination of sub-resolution photo-ionization modelling for HII regions and Monte Carlo radiative transfer calculations is employed to estimate the dust-attenuated spectral energy distributions (SEDs) of high-redshift galaxies ($z\gtrsim5.5$). We show that the derived photometric properties such as the ultraviolet (UV) luminosity function and the UV continuum slopes match observationally inferred values, demonstrating the accuracy of the SED modelling. We provide fits to the luminosity--star formation rate relation (L-SFR) for the brightest emission lines and find that important differences exist between the derived scaling relations and the widely used low-$z$ ones because the interstellar medium of reionization era galaxies is generally less metal-enriched than in their low redshift counterparts. We use these relations to construct line intensity maps of nebular emission lines and cross correlate with the 21 cm emission. Interestingly, the wavenumber at which the correlation switches sign ($k_\mathrm{transition}$) depends heavily on the reionization model and to a lesser extent on the targeted emission line, which is consistent with the picture that $k_\mathrm{transition}$ probes the typical sizes of ionized regions. The derived scaling relations and intensity maps represent a timely state-of-the-art framework for forecasting and interpreting results from current and upcoming LIM experiments.

preprint2022arXiv

The THESAN project: properties of the intergalactic medium and its connection to Reionization-era galaxies

The high-redshift intergalactic medium (IGM) and the primeval galaxy population are rapidly becoming the new frontier of extra-galactic astronomy. We investigate the IGM properties and their connection to galaxies at $z\geq5.5$ under different assumptions for the ionizing photon escape and the nature of dark matter, employing our novel THESAN radiation-hydrodynamical simulation suite, designed to provide a comprehensive picture of the emergence of galaxies in a full reionization context. Our simulations have realistic `late' reionization histories, match available constraints on global IGM properties and reproduce the recently-observed rapid evolution of the mean free path of ionizing photons. We additionally examine high-z Lyman-$α$ transmission. The optical depth evolution is consistent with data, and its distribution suggests an even-later reionization than simulated, although with a strong sensitivity to the source model. We show that the effects of these two unknowns can be disentangled by characterising the spectral shape and separation of Lyman-$α$ transmission regions, opening up the possibility to observationally constrain both. For the first time in simulations, THESAN reproduces the modulation of the Lyman-$α$ flux as a function of galaxy distance, demonstrating the power of coupling a realistic galaxy formation model with proper radiation-hydrodynamics. We find this feature to be extremely sensitive on the timing of reionization, while being relatively insensitive to the source model. Overall, THESAN produces a realistic IGM and galaxy population, providing a robust framework for future analysis of the high-z Universe.

preprint2022arXiv

Unifying Sunyaev-Zel'dovich and X-ray predictions from clusters to galaxy groups: the impact of X-ray mass estimates on the $Y-M$ scaling relation

One of the main limitations in precision cluster cosmology arises from systematic errors and uncertainties in estimating cluster masses. Using the Mock-X pipeline, we produce synthetic X-ray images and derive cluster and galaxy group X-ray properties for a sample of over 30,000 simulated galaxy groups and clusters with $M_{\rm 500crit}$ between $10^{12}$ and $2\times 10^{15}$ M$_{\odot}$ in IllustrisTNG. We explore the similarities and differences between IllustrisTNG predictions of the Sunyaev-Zel'dovich and X-ray scaling relations with mass. We find a median hydrostatic mass bias $b = 0.125 \pm 0.003$ for $M_{\rm 500crit}$ $>10^{13}$ M$_{\odot}$. The bias increases to $b = 0.17 \pm 0.004$ when masses are derived from synthetic X-ray observations. We model how different underlying assumptions about the dependence of $Y_{\rm X}$ on halo mass can generate biases in the observed $Y_{\rm SZ} - M_{Y_{\rm X}}$ scaling relation. In particular, the simplifying assumption that $Y_{\rm X} - M_{\rm tot}$ is self-similar at all mass scales largely hides the break in $Y_{\rm SZ} - M_{\rm tot}$ and overestimates $Y_{\rm SZ}$ at galaxy and groups scales. We show that calibrating the $Y_{\rm X}-$mass proxy using a new model for a smoothly broken power law reproduces the true underlying $Y_{\rm SZ} - M_{\rm tot}$ scaling relation with high accuracy. Moreover, $M_{Y_{\rm X}}$ estimates calibrated with this method lead to $Y_{\rm SZ} - M_{Y_{\rm X}}$ predictions that are not biased by the presence of lower mass clusters or galaxy groups in the sample. Finally, we show that our smoothly broken power law model provides a robust way to derive the $Y_{\rm X}-$mass proxy, significantly reducing the level of mass bias for clusters, groups, and galaxies.

preprint2022arXiv

UV to submillimetre luminosity functions of TNG50 galaxies

We apply the radiative transfer (RT) code SKIRT on a sample of ~14000 low-redshift (z<= 0.1) galaxies extracted from the TNG50 simulation to enable an apples-to-apples comparison with observations. The RT procedure is calibrated via comparison of a subsample of TNG50 galaxies with the DustPedia observational sample: we compare several luminosity and colour scaling relations and spectral energy distributions in different specific SFR bins. We consistently derive galaxy luminosity functions for the TNG50 simulation in 14 broadband filters from UV to submillimetre wavelengths and investigate the effects of the aperture, orientation, radiative transfer recipe, and numerical resolution. We find that, while our TNG50+RT fiducial model agrees well with the observed luminosity functions at the knee (+/- 0.04 dex typical agreement), the TNG50+RT luminosity functions evaluated within 5R_1/2 are generally higher than observed at both the faint and bright ends, by 0.004 (total IR)-0.27 (UKIDSS H) dex and 0.12 (SPIRE250)-0.8 (GALEX FUV) dex, respectively. A change in the aperture does affect the bright end of the luminosity function, easily by up to 1 dex depending on the choice. However, we also find that the galaxy luminosity functions of a worse-resolution run of TNG50 (TNG50-2, with 8 times worse mass resolution than TNG50, similar to TNG100) are in better quantitative agreement with observational constraints. Finally, we publicly release the photometry for the TNG50 sample in 53 broadbands from FUV to submillimetre, in three orientations and four apertures, as well as galaxy spectral energy distributions.

preprint2021arXiv

Finding universal relations in subhalo properties with artificial intelligence

We use a generic formalism designed to search for relations in high-dimensional spaces to determine if the total mass of a subhalo can be predicted from other internal properties such as velocity dispersion, radius, or star-formation rate. We train neural networks using data from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project and show that the model can predict the total mass of a subhalo with high accuracy: more than 99% of the subhalos have a predicted mass within 0.2 dex of their true value. The networks exhibit surprising extrapolation properties, being able to accurately predict the total mass of any type of subhalo containing any kind of galaxy at any redshift from simulations with different cosmologies, astrophysics models, subgrid physics, volumes, and resolutions, indicating that the network may have found a universal relation. We then use different methods to find equations that approximate the relation found by the networks and derive new analytic expressions that predict the total mass of a subhalo from its radius, velocity dispersion, and maximum circular velocity. We show that in some regimes, the analytic expressions are more accurate than the neural networks. We interpret the relation found by the neural network and approximated by the analytic equation as being connected to the virial theorem.

preprint2021arXiv

High and low Sérsic index bulges in Milky Way- and M31-like galaxies: origin and connection to the bar with TNG50

We study bulge formation in MW/M31-like galaxies in a $Λ$-cold dark matter scenario, focusing on the origin of high- and low-Sersic index bulges. For this purpose we use TNG50, a simulation of the IllustrisTNG project that combines a resolution of $\sim 8 \times 10^4 M_{\odot}$ in stellar particles with a cosmological volume 52 cMpc in extent. We parametrize bulge surface brightness profiles by the Sérsic index and the bulge-to-total (B/T) ratio obtained from two-component photometric decompositions. In our sample of 287 MW/M31-like simulated galaxies, $17.1\%$ of photometric bulges exhibit high-Sérsic indices and $82.9\%$ show low-Sérsic indices. We study the impact that the environment, mergers and bars have in shaping the surface brightness profiles. We explore two different definitions for local environment and find no correlation between bulge properties and the environment where they reside. Simulated galaxies with higher Sérsic indices show, on average, a higher fraction of ex-situ stars in their kinematically selected bulges. For this bulge population the last significant merger (total mass ratio $m_{\rm sat}/m_{\rm host} > 0.1$) occurs, on average, at later times. However, a substantial fraction of low-Sérsic index bulges also experience a late significant merger. We find that bars play an important role in the development of the different types of photometric bulges. We show that the fraction of simulated galaxies with strong bars is smaller for the high- than for the low-Sérsic index population, reaching differences of $20\%$ at $z > 1$. Simulated galaxies with high fractions of ex-situ stars in the bulge do not develop strong bars. Conversely, simulated galaxies with long-lived strong bars have bulges with ex-situ fractions, $f_{\rm ex-situ} < 0.2$.

preprint2021arXiv

Hot and counter-rotating star-forming disk galaxies in IllustrisTNG and their real-world counterparts

A key feature of a large population of low-mass, late-type disk galaxies are star-forming disks with exponential light distributions. They are typically also associated with thin and flat morphologies, blue colors, and dynamically cold stars moving along circular orbits within co-planar thin gas disks. However, the latter features do not necessarily always imply the former, in fact, a variety of different kinematic configurations do exist. In this work, we use the cosmological hydrodynamical IllustrisTNG Simulation to study the nature and origin of dynamically hot, sometimes even counter-rotating, star-forming disk galaxies in the lower stellar mass range (between $5\times 10^9\,\mathrm{M_{\odot}}$ and $2\times 10^{10}\,\mathrm{M_{\odot}}$). We find that being dynamically hot arises in most cases as an induced transient state, for example due to galaxy interactions and merger activities, rather than as an age-dependent evolutionary phase of star-forming disk galaxies. The dynamically hot but still actively star-forming disks show a common feature of hosting kinematically misaligned gas and stellar disks, and centrally concentrated on-going star formation. The former is often accompanied by disturbed gas morphologies, while the latter is reflected in low gas and stellar spins in comparison to their dynamically cold, normal disk counterparts. Interestingly, observed galaxies from MaNGA with kinematic misalignment between gas and stars show remarkably similar general properties as the IllustrisTNG galaxies, and therefore are plausible real-world counterparts. In turn, this allows us to make predictions for the stellar orbits and gas properties of these misaligned galaxies.

preprint2021arXiv

On the formation of massive quiescent galaxies with diverse morphologies in the TNG50 simulation

Observations have shown that the star-formation activity and the morphology of galaxies are closely related, but the underlying physical connection is not well understood. Using the TNG50 simulation, we explore the quenching and the morphological evolution of the 102 massive quiescent galaxies in the mass range of $10.5<\log(M_{\rm stellar}/M_{\odot})<11.5$ selected at $z=0$. The morphology of galaxies is quantified based on their kinematics, and we measure the quenching timescale of individual galaxies directly from star formation history. We show that galaxies tend to be quenched more rapidly if they: (i) are satellites in massive halos, (ii) have lower star-forming gas fractions, or (iii) inject a larger amount of black hole kinetic feedback energy. By following the global evolutionary pathways, we conclude that quiescent discs are mainly disc galaxies that are recently and slowly quenched. Approximately half of the quiescent ellipticals at $z=0$ are rapidly quenched at higher redshifts while still disc-like. While being quiescent, they gradually become more elliptical mostly by disc heating, yet these ellipticals still retain some degree of rotation. The other half of quiescent ellipticals with the most random motion-dominated kinematics build up large spheroidal components before quenching primarily by mergers, or in some cases, misaligned gas accretion. However, the mergers that contribute to morphological transformation do not immediately quench galaxies in many cases. In summary, we find that quenching and morphological transformation are decoupled. We conclude that the TNG black hole feedback -- in combination with the stochastic merger history of galaxies -- leads to a large diversity of quenching timescales and a rich morphological landscape.

preprint2021arXiv

Submillimetre galaxies in cosmological hydrodynamical simulations -- an opportunity for constraining feedback models

Submillimetre galaxies (SMGs) have long posed a challenge for theorists, and self-consistently reproducing the properties of the SMG population in a large-volume cosmological hydrodynamical simulation has not yet been achieved. We use a scaling relation derived from previous simulations plus radiative transfer calculations to predict the submm flux densities of simulated SMGs drawn from cosmological simulations from the Illustris and IllustrisTNG projects based on the simulated galaxies' star formation rates (SFRs) and dust masses and compare the predicted number counts with observations. We find that the predicted SMG number counts based on IllustrisTNG are significantly less than observed (more than 1 dex at $S_{850} \gtrsim 4$ mJy). The simulation from the original Illustris project yields more SMGs than IllustrisTNG: the predicted counts are consistent with those observed at both $S_{850} \lesssim 5$ mJy and $S_{850} \gtrsim 9$ mJy and only a factor of $\sim 2$ lower than observed at intermediate flux densities. The redshift distribution of SMGs with $S_{850} > 3$ mJy in IllustrisTNG is consistent with the observed distribution, whereas the Illustris redshift distribution peaks at significantly lower redshift (1.5 vs. 2.8). We demonstrate that IllustrisTNG hosts fewer SMGs than Illustris because in the former, high-mass ($M_{\star} \sim 10^{11} \, \text{M}_{\odot}$) $z \sim 2-3$ galaxies have lower dust masses and SFRs than in Illustris owing to differences in the sub-grid models for stellar or/and active galactic nucleus (AGN) feedback between the two simulations (we unfortunately cannot isolate the specific cause(s) post hoc). Our results demonstrate that because our method enables predicting SMG number counts in post-processing with a negligible computational expense, SMGs can provide useful constraints for tuning sub-grid models in future large-volume cosmological simulations.

preprint2021arXiv

The IllustrisTNG Simulations: Public Data Release

We present the full public release of all data from the TNG50, TNG100 and TNG300 simulations of the IllustrisTNG project. IllustrisTNG is a suite of large volume, cosmological, gravo-magnetohydrodynamical simulations run with the moving-mesh code Arepo. TNG includes a comprehensive model for galaxy formation physics, and each TNG simulation self-consistently solves for the coupled evolution of dark matter, cosmic gas, luminous stars, and supermassive blackholes from early time to the present day, z=0. Each of the flagship runs -- TNG50, TNG100, and TNG300 -- are accompanied by lower-resolution and dark-matter only counterparts, and we discuss scientific and numerical cautions and caveats relevant when using TNG. Full volume snapshots are available at 100 redshifts; halo and subhalo catalogs at each snapshot and merger trees are also released. The data volume now directly accessible online is ~1.1 PB, including 2,000 full volume snapshots and ~110,000 high time-resolution subbox snapshots. Data access and analysis examples are available in IDL, Python, and Matlab. We describe improvements and new functionality in the web-based API, including on-demand visualization and analysis of galaxies and halos, exploratory plotting of scaling relations and other relationships between galactic and halo properties, and a new JupyterLab interface. This provides an online, browser-based, near-native data analysis platform which supports user computation with fully local access to TNG data, alleviating the need to download large simulated datasets.

preprint2021arXiv

The impact of galactic feedback on the shapes of dark-matter haloes

We quantify the impact of galaxy formation on dark matter halo shapes using cosmological simulations at redshift $z=0$. The haloes are drawn from the IllustrisTNG project, a suite of magneto-hydrodynamic simulations of galaxies. We focus on haloes of mass $10^{10-14} M_\odot$ from the 50-Mpc (TNG50) and 100-Mpc (TNG100) boxes, and compare them to dark matter-only (DMO) analogues and other simulations e.g. NIHAO and Eagle. We further quantify the prediction uncertainty by varying the baryonic feedback models in a series of smaller 25 Mpc $h^{-1}$ boxes. We find that: (i) galaxy formation results in rounder haloes compared to the DMO simulations, in qualitative agreement with past hydrodynamic models. Haloes of mass $\approx 2\times 10^{12} M_\odot$ are most spherical, with an average minor-to-major axis ratio of $\left< s \right> \approx 0.75$ in the inner halo, an increase of 40 per cent compared to their DMO counterparts. No significant change in halo shape is found for low-mass $10^{10} M_\odot$ haloes; (ii) stronger feedback, e.g. increasing galactic wind speed, reduces the impact of baryons; (iii) the inner halo shape correlates with the stellar mass fraction, which can explain the dependence of halo shapes on different feedback models; (iv) the fiducial and weaker feedback models are most consistent with observational estimates of the Milky Way halo shape. Yet, at fixed halo mass, very diverse and possibly unrealistic feedback models all predict inner halo shapes that are closer to one another than to the DMO results. This implies that a larger observational sample would be required to statistically distinguish between different baryonic prescriptions due to large halo-to-halo variation in halo shapes.

preprint2020arXiv

A missing outskirts problem? Comparisons between stellar halos in the Dragonfly Nearby Galaxies Survey and the TNG100 simulation

Low surface brightness galactic stellar halos provide a challenging but promising path towards unraveling the past assembly histories of individual galaxies. Here, we present detailed comparisons between the stellar halos of Milky Way-mass disk galaxies observed as part of the Dragonfly Nearby Galaxies Survey (DNGS) and stellar mass-matched galaxies in the TNG100 run of the IllustrisTNG project. We produce stellar mass maps as well as mock $g$ and $r$-band images for randomly oriented simulated galaxies, convolving the latter with the Dragonfly PSF and taking care to match the background noise, surface brightness limits and spatial resolution of DNGS. We measure azimuthally averaged stellar mass density and surface brightness profiles, and find that the DNGS galaxies generally have less stellar mass (or light) at large radii (>20 kpc) compared to their mass-matched TNG100 counterparts, and that simulated galaxies with similar surface density profiles tend to have low accreted mass fractions for their stellar mass. We explore potential solutions to this apparent "missing outskirts problem" by implementing several ad-hoc adjustments within TNG100 at the stellar particle level. Although we are unable to identify any single adjustment that fully reconciles the differences between the observed and simulated galaxy outskirts, we find that artificially delaying the disruption of satellite galaxies and reducing the spatial extent of in-situ stellar populations result in improved matches between the outer profile shapes and stellar halo masses, respectively. Further insight can be achieved with higher resolution simulations that are able to better resolve satellite accretion, and with larger samples of observed galaxies.

preprint2020arXiv

A redshift-dependent IRX-$β$ dust attenuation relation for TNG50 galaxies

We study the relation between the UV-slope, $β$, and the ratio between the infrared- and UV-luminosities (IRX) of galaxies from TNG50, the latest installment of the IllustrisTNG galaxy formation simulations. We select 7280 star-forming main-sequence (SFMS) galaxies with stellar mass $10^9M_\odot$ at redshifts $0 \leq z \leq 4$ and perform radiative transfer with SKIRT to model effects of interstellar medium dust on the emitted stellar light. Assuming a Milky Way (MW) dust type and a dust-to-metal ratio of 0.3, we find that TNG50 SFMS galaxies generally agree with observationally-derived IRX-$β$ relations at $z \lesssim 1$. However, we find a redshift-dependent systematic offset with respect to empirically-derived local relations, with the TNG50 IRX-$β$ relation shifting towards lower $β$ and steppening at higher redshifts. This is partially driven by variations in the dust-uncorrected UV-slope of galaxies, due to different star-formation histories of galaxies selected at different cosmic epochs; we suggest the remainder of the effect is caused by differences in the effective dust attenuation curves (EDACs) of galaxies as a function of redshift. We find a typical galaxy-to-galaxy variation of 0.3 dex in IRX at fixed $β$, correlated with intrinsic galaxy properties: galaxies with higher star-formation rates, star-formation efficiencies, gas metallicities and stellar masses exhibit larger IRX values. We demonstrate a degeneracy between stellar age, dust geometry and dust composition: $z=4$ galaxies with a Small Magellanic Cloud dust type follow the same IRX-$β$ relation as low-redshift galaxies with MW dust. We provide a redshift-dependent fitting function for the IRX-$β$ relation for MW dust based on our models.

preprint2020arXiv

Atomic hydrogen in IllustrisTNG galaxies: the impact of environment parallelled with local 21-cm surveys

We investigate the influence of environment on the cold-gas properties of galaxies at z=0 within the TNG100 cosmological, magnetohydrodynamic simulation, part of the IllustrisTNG suite. We extend previous post-processing methods for breaking gas cells into their atomic and molecular phases, and build detailed mocks to comprehensively compare to the latest surveys of atomic hydrogen (HI) in nearby galaxies, namely ALFALFA and xGASS. We use TNG100 to explore the HI content, star formation activity, and angular momentum of satellite galaxies, each as a function of environment, and find that satellites are typically a factor of ~3 poorer in HI than centrals of the same stellar mass, with the exact offset depending sensitively on parent halo mass. Due to the large physical scales on which HI measurements are made (~45--245 kpc), contributions from gas not bound to the galaxy of interest but in the same line of sight crucially lead to larger HI mass measurements in the mocks in many cases, ultimately aligning with observations. This effect is mass-dependent and naturally greater for satellites than centrals, as satellites are never isolated by definition. We also show that HI stripping in TNG100 satellites is closely accompanied by quenching, in tension with observational data that instead favour that HI is preferentially stripped before star formation is reduced.

preprint2020arXiv

Cosmological insights into the assembly of the radial and compact stellar halo of the Milky Way

Recent studies using Gaia DR2 have identified a massive merger in the history of the Milky Way (MW) whose debris is markedly radial and counterrotating. This event, known as the Gaia-Enceladus/Gaia-Sausage (GE/GS), is also hypothesized to have built the majority of the inner stellar halo. We use the cosmological hydrodynamic simulation Illustris to place this merger in the context of galaxy assembly within $Λ$CDM. From $\sim$150 MW analogs, $\sim 80 \%$ have experienced at least one merger of similar mass and infall time as GE/GS. Within this sample, 37 have debris as radial as that of the GE/GS, which we dub the Ancient Radial Mergers (ARMs). Counterrotation is not rare among ARMs, with $43 \%$ having $>40 \%$ of their debris in counterrotating orbits. However, the compactness inferred for the GE/GS debris, given its large $β$ and its substantial contribution to the stellar halo, is more difficult to reproduce. The median radius of ARM debris is r$_{*,deb}\simeq 45$kpc, while GE/GS is thought to be mostly contained within $r\sim 30$ kpc. For most MW analogs, a few mergers are required to build the inner stellar halo, and ARM debris only accounts for $\sim 12 \%$ of inner accreted stars. Encouragingly, we find one ARM that is both compact and dominates the inner halo of its central, making it our best GE/GS analog. Interestingly, this merger deposits a significant number of stars (M$_*\simeq1.5 \times 10^9 M_\odot$) in the outer halo, suggesting that an undiscovered section of GE/GS may await detection.

preprint2020arXiv

Dust entrainment in galactic winds

Winds driven by stellar feedback are an essential part of the galactic ecosystem and are the main mechanism through which low-mass galaxies regulate their star formation. These winds are generally observed to be multi-phase with detections of entrained neutral and molecular gas. They are also thought to enrich the circum-galactic medium around galaxies with metals and dust. This ejected dust encodes information about the integrated star formation and outflow history of the galaxy. It is therefore, important to understand how much dust is entrained and driven out of the disc by galactic winds. Here we demonstrate that stellar feedback is efficient in driving dust-enriched winds and eject enough material to account for the amount of extraplanar dust observed in nearby galaxies. The amount of dust in the wind depends on the sites from where they are launched, with dustier galaxies launching more dust enriched outflows. Moreover, the outflowing cold-dense gas is significantly more dust-enriched than the volume filling hot tenuous material, naturally reproducing the complex multiphase structure of the outflowing wind observed in nearby galaxies. These results provide an important new insight into the dynamics, structure, and composition of galactic winds and their role in determining the dust content of the extragalactic gas in galaxies.

preprint2020arXiv

Early-type galaxy density profiles from IllustrisTNG: I. Galaxy correlations and the impact of baryons

We explore the isothermal total density profiles of early-type galaxies (ETGs) in the IllustrisTNG simulation. For the selected 559 ETGs at $z = 0$ with stellar mass $10^{10.7}\mathrm{M}_{\odot} \leqslant M_{\ast} \leqslant 10^{11.9}\mathrm{M}_{\odot}$, the total power-law slope has a mean of $\langleγ^{\prime}\rangle = 2.011 \pm 0.007$ and a scatter of $σ_{γ^{\prime}} = 0.171$ over the radial range 0.4 to 4 times the stellar half mass radius. Several correlations between $γ^{\prime}$ and galactic properties including stellar mass, effective radius, stellar surface density, central velocity dispersion, central dark matter fraction and in-situ-formed stellar mass ratio are compared to observations and other simulations, revealing that IllustrisTNG reproduces many correlation trends, and in particular, $γ^{\prime}$ is almost constant with redshift below $z = 2$. Through analyzing IllustrisTNG model variations we show that black hole kinetic winds are crucial to lowering $γ^{\prime}$ and matching observed galaxy correlations. The effects of stellar winds on $γ^{\prime}$ are subdominant compared to AGN feedback, and differ due to the presence of AGN feedback from previous works. The density profiles of the ETG dark matter halos are well-described by steeper-than-NFW profiles, and they are steeper in the full physics (FP) run than their counterparts in the dark matter only (DMO) run. Their inner density slopes anti-correlates (remain constant) with the halo mass in the FP (DMO) run, and anti-correlates with the halo concentration parameter $c_{200}$ in both types of runs. The dark matter halos of low-mass ETGs are contracted whereas high-mass ETGs are expanded, suggesting that variations in the total density profile occur through the different halo responses to baryons.

preprint2020arXiv

Ejective and preventative: the IllustrisTNG black hole feedback and its effects on the thermodynamics of the gas within and around galaxies

Supermassive black holes (SMBHs) which reside at the centres of galaxies can inject vast amounts of energy into the surrounding gas and are thought to be a viable mechanism to quench star-formation in massive galaxies. Here we study the $10^{9\textrm{--}12.5}\,\mathrm{M_\odot}$ stellar mass central galaxy population of the IllustrisTNG simulation, specifically the TNG100 and TNG300 volumes at \zeq{0}, and show how the three components -- SMBH, galaxy, and circumgalactic medium (CGM) -- are interconnected in their evolution. We find that gas entropy is a sensitive diagnostic of feedback injection. In particular, we demonstrate how the onset of the low-accretion BH feedback mode, realised in the IllustrisTNG model as a kinetic, BH-driven wind, leads not only to star-formation quenching at stellar masses $\gtrsim10^{10.5}\mathrm{M_\odot}$ but also to a change in thermodynamic properties of the (\emph{non}-star-forming) gas, both within the galaxy and beyond. The IllustrisTNG kinetic feedback from SMBHs increases the average gas entropy, within the galaxy and in the CGM, lengthening typical gas cooling times from $10\textrm{--}100\,\mathrm{Myr}$ to $1\textrm{--}10\,\mathrm{Gyr}$, effectively ceasing ongoing star-formation and inhibiting radiative cooling and future gas accretion. In practice, the same AGN feedback channel is simultaneously `ejective' and `preventative' and leaves an imprint on the temperature, density, entropy, and cooling times also in the outer reaches of the gas halo, up to distances of several hundred kiloparsecs. In the IllustrisTNG model, a long-lasting quenching state can occur for a heterogeneous CGM, whereby the hot and dilute CGM gas of quiescent galaxies contains regions of low-entropy gas with short cooling times.

preprint2020arXiv

First results from SMAUG: Characterization of Multiphase Galactic Outflows from a Suite of Local Star-Forming Galactic Disk Simulations

Large scale outflows in star-forming galaxies are observed to be ubiquitous, and are a key aspect of theoretical modeling of galactic evolution in a cosmological context, the focus of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project. Gas blown out from galactic disks, similar to gas within galaxies, consists of multiple phases with large contrasts of density, temperature, and other properties. To study multiphase outflows as emergent phenomena, we run a suite of ~pc-resolution local galactic disk simulations using the TIGRESS framework. Explicit modeling of the interstellar medium (ISM), including star formation and self-consistent radiative heating plus supernova feedback, regulates ISM properties and drives the outflow. We investigate the scaling of outflow mass, momentum, energy, and metal loading factors with galactic disk properties, including star formation rate (SFR) surface density (Σ_SFR~10^{-4}-1 M_sun/kpc^2/yr), gas surface density (~1-100 M_sun/pc^2), and total midplane pressure (or weight) (~10^3-10^6 k_B cm^{-3} K). The main components of outflowing gas are mass-delivering cool gas (T~10^4 K) and energy/metal-delivering hot gas (T~10^6 K). Cool mass outflow rates measured at outflow launch points (one or two scale heights) are 1-100 times the SFR (decreasing with Σ_SFR), although in massive galaxies most mass falls back due to insufficient outflow velocity. The hot galactic outflow carries mass comparable to 10% of the SFR, together with 10-20% of the energy and 30-60% of the metal mass injected by SN feedback. The characteristic outflow velocities of both phases scale very weakly with SFR, as v_out \propto Σ_SFR^{0.1~0.2}, consistent with observations. Importantly, our analysis demonstrates that in any physically-motivated cosmological wind model, it is crucial to include at least two distinct thermal wind components.

preprint2020arXiv

First results from SMAUG: Insights into star formation conditions from spatially-resolved ISM properties in TNG50

Physical and chemical properties of the interstellar medium (ISM) at sub-galactic ($\sim$kpc) scales play an indispensable role in controlling the ability of gas to form stars. As part of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project, in this paper, we use the TNG50 cosmological simulation to explore the physical parameter space of 8 resolved ISM properties in star-forming regions to constrain the areas of this hyperspace over which most star-forming environments exist. We deconstruct our simulated galaxies spanning a wide range of mass (M$_\star = 10^{7-11}$ M$_\odot$) and redshift ($0 \leq z \leq 3$) into kpc-sized regions, and statistically analyze the gas/stellar surface densities, gas metallicity, vertical stellar velocity dispersion, epicyclic frequency and dark-matter volumetric density representative of each region in the context of their star formation activity and galactic environment (radial galactocentric location). By examining the star formation rate (SFR) weighted distributions of these properties, we show that stars primarily form in two spatially distinct environmental regimes, which are brought about by an underlying bi-component radial SFR surface density profile in galaxies. We examine how the relative prominence of these two regimes depends on host galaxy mass and cosmic time. We also compare our findings with those from integral field spectroscopy observations and achieve a good overall agreement. Further, using dimensionality reduction, we characterise the aforementioned hyperspace to reveal a high-degree of multicollinearity in relationships amongst ISM properties that drive the distribution of star formation at kpc-scales. Based on this, we show that a reduced 3D representation underpinned by a multi-variate radius relationship is sufficient to capture most of the variance in the original 8D space.

preprint2020arXiv

High redshift JWST predictions from IllustrisTNG: II. Galaxy line and continuum spectral indices and dust attenuation curves

We present predictions for high redshift ($z=2-10$) galaxy populations based on the IllustrisTNG simulation suite and a full Monte Carlo dust radiative transfer post-processing. Specifically, we discuss the ${\rm H}_α$ and ${\rm H}_β$ + $[\rm O \,III]$ luminosity functions up to $z=8$. The predicted ${\rm H}_β$ + $[\rm O \,III]$ luminosity functions are consistent with present observations at $z\lesssim 3$ with $\lesssim 0.1\,{\rm dex}$ differences in luminosities. However, the predicted ${\rm H}_α$ luminosity function is $\sim 0.3\,{\rm dex}$ dimmer than the observed one at $z\simeq 2$. Furthermore, we explore continuum spectral indices, the Balmer break at $4000$Å(D4000) and the UV continuum slope $β$. The median D4000 versus sSFR relation predicted at $z=2$ is in agreement with the local calibration despite a different distribution pattern of galaxies in this plane. In addition, we reproduce the observed $A_{\rm UV}$ versus $β$ relation and explore its dependence on galaxy stellar mass, providing an explanation for the observed complexity of this relation. We also find a deficiency in heavily attenuated, UV red galaxies in the simulations. Finally, we provide predictions for the dust attenuation curves of galaxies at $z=2-6$ and investigate their dependence on galaxy colors and stellar masses. The attenuation curves are steeper in galaxies at higher redshifts, with bluer colors, or with lower stellar masses. We attribute these predicted trends to dust geometry. Overall, our results are consistent with present observations of high redshift galaxies. Future JWST observations will further test these predictions.

preprint2020arXiv

Kinematic decomposition of IllustrisTNG disk galaxies: morphology and relation with morphological structures

We recently developed an automated method, auto-GMM to decompose simulated galaxies. It extracts kinematic structures in an accurate, efficient, and unsupervised way. We use auto-GMM to study the stellar kinematic structures of disk galaxies from the TNG100 run of IllustrisTNG. We identify four to five structures that are commonly present among the diverse galaxy population. Structures having strong to moderate rotation are defined as cold and warm disks, respectively. Spheroidal structures dominated by random motions are classified as bulges or stellar halos, depending on how tightly bound they are. Disky bulges are structures that have moderate rotation but compact morphology. Across all disky galaxies and accounting for the stellar mass within 3 half-mass radii, the kinematic spheroidal structures, obtained by summing up stars of bulges and halos, contribute ~45% of the total stellar mass, while the disky structures constitute 55%. This study also provides important insights about the relationship between kinematically and morphologically derived galactic structures. Comparing the morphology of kinematic structures with that of traditional bulge+disk decomposition, we conclude: (1) the morphologically decomposed bulges are composite structures comprised of a slowly rotating bulge, an inner halo, and a disky bulge; (2) kinematically disky bulges, akin to what are commonly called pseudo bulges in observations, are compact disk-like components that have rotation similar to warm disks; (3) halos contribute almost 30% of the surface density of the outer part of morphological disks when viewed face-on; and (4) both cold and warm disks are often truncated in central regions.

preprint2020arXiv

Limitations to the "basic'' HOD model and beyond

We make use of 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 $\sim$ 15\% on scales of $1 \ {\rm Mpc}/h < r < 20 \ {\rm Mpc}/h$, 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, $σ^2 R_{\rm halfmass}$, 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.

preprint2020arXiv

Redshift Evolution of the Fundamental Plane Relation in the IllustrisTNG Simulation

We investigate the fundamental plane (FP) evolution of early-type galaxies in the IllustrisTNG-100 simulation (TNG100) from redshift $z=0$ to $z=2$. We find that a tight plane relation already exists as early as $z=2$. Its scatter stays as low as $\sim 0.08$ dex across this redshift range. Both slope parameters $b$ and $c$ (where $R \propto σ^b I^c$ with $R$, $σ$, and $I$ being the typical size, velocity dispersion, and surface brightness) of the plane evolve mildly since $z=2$, roughly consistent with observations. The FP residual $\rm Res$ ($\equiv\,a\,+\,b\log σ\,+\,c\log I\,-\,\log R$, where $a$ is the zero point of the FP) is found to strongly correlate with stellar age, indicating that stellar age can be used as a crucial fourth parameter of the FP. However, we find that $4c+b+2=δ$, where $δ\sim 0.8$ for FPs in TNG, rather than zero as is typically inferred from observations. This implies that a tight power-law relation between the dynamical mass-to-light ratio $M_{\rm dyn}/L$ and the dynamical mass $M_{\rm dyn}$ (where $M_{\rm dyn}\equiv 5σ^2R/G$, with $G$ being the gravitational constant) is not present in the TNG100 simulation. Recovering such a relation requires proper mixing between dark matter and baryons, as well as star formation occurring with correct efficiencies at the right mass scales. This represents a powerful constraint on the numerical models, which has to be satisfied in future hydrodynamical simulations.

preprint2020arXiv

Resolving small-scale cold circumgalactic gas in TNG50

We use the high-resolution TNG50 cosmological magnetohydrodynamical simulation to explore the properties and origin of cold circumgalactic medium (CGM) gas around massive galaxies (M* > 10^11 Msun) at intermediate redshift (z~0.5). We discover a significant abundance of small-scale, cold gas structure in the CGM of 'red and dead' elliptical systems, as traced by neutral HI and MgII. Halos can host tens of thousands of discrete absorbing cloudlets, with sizes of order a kpc or smaller. With a Lagrangian tracer analysis, we show that cold clouds form due to strong drho/rho >> 1 gas density perturbations which stimulate thermal instability. These local overdensities trigger rapid cooling from the hot virialized background medium at ~10^7 K to radiatively inefficient ~10^4 K clouds, which act as cosmologically long-lived, 'stimulated cooling' seeds in a regime where the global halo does not satisfy the classic tcool/tff < 10 criterion. Furthermore, these small clouds are dominated by magnetic rather than thermal pressure, with plasma beta << 1, suggesting that magnetic fields may play an important role. The number and total mass of cold clouds both increase with resolution, and the ~8x10^4 Msun cell mass of TNG50 enables the ~few hundred pc, small-scale CGM structure we observe to form. Finally, we make a preliminary comparison against observations from the COS-LRG, LRG-RDR, COS-Halos, and SDSS LRG surveys. We broadly find that our recent, high-resolution cosmological simulations produce sufficiently high covering fractions of extended, cold gas as observed to surround massive galaxies.

preprint2020arXiv

Separate Universe Simulations with IllustrisTNG: baryonic effects on power spectrum responses and higher-order statistics

We measure power spectrum response functions in the presence of baryonic physical processes using separate universe simulations with the IllustrisTNG galaxy formation model. The response functions describe how the small-scale power spectrum reacts to long-wavelength perturbations and they can be efficiently measured with the separate universe technique by absorbing the effects of the long modes into a modified cosmology. Specifically, we focus on the total first-order matter power spectrum response to an isotropic density fluctuation $R_1(k,z)$, which is fully determined by the logarithmic derivative of the nonlinear matter power spectrum ${\rm dln}P_m(k,z)/{\rm dln}k$ and the growth-only response function $G_1(k,z)$. We find that $G_1(k,z)$ is not affected by the baryonic physical processes in the simulations at redshifts $z < 3$ and on all scales probed ($k \lesssim 15h/{\rm Mpc}$, i.e. length scales $\gtrsim 0.4 {\rm Mpc}/h$). In practice, this implies that the power spectrum fully specifies the baryonic dependence of its response function. Assuming an idealized lensing survey setup, we evaluate numerically the baryonic impact on the squeezed-lensing bispectrum and the lensing super-sample power spectrum covariance, which are given in terms of responses. Our results show that these higher-order lensing statistics can display varying levels of sensitivity to baryonic effects compared to the power spectrum, with the squeezed-bispectrum being the least sensitive. We also show that ignoring baryonic effects on lensing covariances slightly overestimates the error budget (and is therefore conservative from the point of view of parameter error bars) and likely has negligible impact on parameter biases in inference analyses.

preprint2020arXiv

Simulating the "hidden giant" in cold and self-interacting dark matter models

We perform a series of controlled N-body simulations to study realizations of the recently discovered Antlia 2 galaxy in cold dark matter (CDM) and self-interacting dark matter (SIDM) scenarios. Our simulations contain six benchmark models, where we vary the initial halo concentration and the self-scattering cross section. We adopt well-motivated initial stellar and halo masses, and our fiducial orbit has a small pericenter. After evolving in the Milky Way's tidal field, the simulated galaxies experience significant mass loss and their stellar distributions expand accordingly. These tidal effects are more prominent if the initial halo concentration is lower and if the self-scattering cross section is larger. Our results show that Antlia 2-like galaxies could be realized in CDM if the halo concentration is low and the stellar distribution is diffuse at the infall time, while these conditions could be relaxed in SIDM. We also find all the simulated galaxies predict approximately the same stellar velocity dispersion after imposing selection criteria for stellar particles. This has important implications for testing dark matter models using tidally disturbed systems.

preprint2020arXiv

The Angular Momentum of the Circumgalactic Medium in the TNG100 Simulation

We present an analysis of the angular momentum content of the circumgalactic medium (CGM) using TNG100, one of the flagship runs of the IllustrisTNG project. We focus on Milky Way-mass halos ($\sim 10^{12} \; M_{\odot}$) at $z=0$ but also analyze other masses and redshifts up to $z=5$. We find that the CGM angular momentum properties are strongly correlated with the stellar angular momentum of the corresponding galaxy: the CGM surrounding high-angular momentum galaxies has a systematically higher angular momentum and is better aligned to the rotational axis of the galaxy itself than the CGM surrounding low-angular momentum galaxies. Both the hot and cold phases of the CGM show this dichotomy, though it is stronger for colder gas. The CGM of high-angular momentum galaxies is characterized by a large wedge of cold gas with rotational velocities at least $\sim1/2$ of the halo's virial velocity, extending out to $\sim 1/2$ of the virial radius, and by biconical polar regions dominated by radial velocities suggestive of galactic fountains; both of these features are absent from the CGM of low-angular momentum galaxies. These conclusions are general to halo masses $\lesssim 10^{12} \; M_{\odot}$ and for $z \lesssim 2$, but they do not apply for more massive halos or at the highest redshift studied. By comparing simulations run with alterations to the fiducial feedback model, we identify the better alignment of the CGM to high-angular momentum galaxies as a feedback-independent effect and the galactic winds as a dominant influence on the CGM's angular momentum.

preprint2020arXiv

The Diversity and Variability of Star Formation Histories in Models of Galaxy Evolution

Understanding the variability of galaxy star formation histories (SFHs) across a range of timescales provides insight into the underlying physical processes that regulate star formation within galaxies. We compile the SFHs of galaxies at $z=0$ from an extensive set of models, ranging from cosmological hydrodynamical simulations (Illustris, IllustrisTNG, Mufasa, Simba, EAGLE), zoom simulations (FIRE-2, g14, and Marvel/Justice League), semi-analytic models (Santa Cruz SAM) and empirical models (UniverseMachine), and quantify the variability of these SFHs on different timescales using the power spectral density (PSD) formalism. We find that the PSDs are well described by broken power-laws, and variability on long timescales ($\gtrsim1$ Gyr) accounts for most of the power in galaxy SFHs. Most hydrodynamical models show increased variability on shorter timescales ($\lesssim300$ Myr) with decreasing stellar mass. Quenching can induce $\sim0.4-1$ dex of additional power on timescales $>1$ Gyr. The dark matter accretion histories of galaxies have remarkably self-similar PSDs and are coherent with the in-situ star formation on timescales $>3$ Gyr. There is considerable diversity among the different models in their (i) power due to SFR variability at a given timescale, (ii) amount of correlation with adjacent timescales (PSD slope), (iii) evolution of median PSDs with stellar mass, and (iv) presence and locations of breaks in the PSDs. The PSD framework is a useful space to study the SFHs of galaxies since model predictions vary widely. Observational constraints in this space will help constrain the relative strengths of the physical processes responsible for this variability.

preprint2020arXiv

The effect of differential accretion on the Gravitational Wave Background and the present day MBH Binary population

Massive black hole binaries (MBHBs) form as a consequence of galaxy mergers. However, it is still unclear whether they typically merge within a Hubble time, and how accretion may affect their evolution. These questions will be addressed by pulsar timing arrays (PTAs), which aim to detect the GW background (GWB) emitted by MBHBs during the last Myrs of inspiral. Here we investigate the influence of differential accretion on MBHB merger rates, chirp masses and the resulting GWB spectrum. We evolve a MBHB sample from the Illustris hydrodynamic cosmological simulation using semi-analytic models and for the first time self-consistently evolve their masses with binary accretion models. In all models, MBHBs coalesce with median total masses up to $1.5 \times 10^8 M_{\odot}$, up to $3-4$ times larger than in models neglecting accretion. In our model with the largest plausible impact, the median mass ratio of coalescing MBHBs increases by a factor $3.6$, the coalescence rate by $52.3\%$, and the GWB amplitude by a factor $4.0$, yielding a dimensionless GWB strain $A_{yr^{-1}} = 1 \times 10^{-15}$. Our model that favours accretion onto the primary MBH reduces the median mass ratio of coalescing MBHBs by a factor of $2.9$, and yields a GWB amplitude $A_{yr^{-1}} = 3.1 \times 10^{-16}$. This is nearly indistinguishable from our model neglecting accretion, despite higher MBHB masses at coalescence. \textbf{We further predict binary separation and mass ratio distributions of stalled MBHBs in the low-redshift universe, and find that these depend sensitively on binary accretion models. This presents the potential for combined EM and GW observational constraints on merger rates and accretion models of MBHB populations.}

preprint2020arXiv

The fate of disk galaxies in IllustrisTNG clusters

We study the stellar morphological evolution of disc galaxies within clusters in the TNG50 and TNG100 runs from the IllustrisTNG simulation suite. We select satellites of masses $10^{9.7} \leq M_{*,z=0}/\text{M}_{\odot} \leq 10^{11.6}$ residing in clusters of masses $10^{14} \lesssim M_{\text{200c,z=0}}/\text{M}_{\odot} \leq 10^{14.6}$ at $z=0$ and that were discs at accretion according to a kinematic morphology indicator (the circularity fraction). These are traced from the time of accretion to $z=0$ and compared to a control sample of central galaxies mass-matched at accretion. Most cluster discs become non-discy by $z=0$, in stark contrast with the control discs, of which a significant fraction remains discy over the same timescales. Cluster discs become non-discy accompanied by gas removal and star formation quenching, loss of dark matter and little growth or a loss of stellar mass. In contrast, control discs transform while also losing gas mass and quenching, but growing significantly in dark matter and stellar mass. Most cluster satellites change morphologies on similar timescales regardless of stellar mass, in $\sim0.5-4$ Gyr after accretion. Cluster discs that experienced more numerous and closer pericentric passages show the largest change in morphology. Morphological change in all cases requires the presence of a gravitational perturbation to drive stellar orbits to non-discy configurations, along with gas removal/heating to prevent replenishment of the disc through continued star-formation. For cluster discs, the perturbation is impulsive tidal shocking at pericentres and not tidal stripping of outer disc stellar material, whereas for control discs, a combination of mergers and AGN feedback appears to be the key driving force behind morphological transformations.

preprint2020arXiv

The Formation History of Subhalos and the Evolution of Satellite Galaxies

Satellites constitute an important fraction of the overall galaxy population and are believed to form in dark matter subhalos. Here we use the cosmological hydrodynamic simulation TNG100 to investigate how the formation histories of subhalos affect the properties and evolution of their host galaxies. We use a scaled formation time ($a_{\rm nf}$) to characterize the mass assembly histories of the subhalos before they are accreted by massive host halos. We find that satellite galaxies in young subhalos (low $a_{\rm nf}$) are less massive and more gas rich, and have stronger star formation and a higher fraction of ex situ stellar mass than satellites in old subhalos (high $a_{\rm nf}$). Furthermore, these low $a_{\rm nf}$ satellites require longer timescales to be quenched as a population than the high $a_{\rm nf}$ counterparts. We find very different merger histories between satellites in fast accretion (FA, $a_{\rm nf}<1.3$) and slow accretion (SA, $a_{\rm nf}>1.3$) subhalos. For FA satellites, the galaxy merger frequency dramatically increases just after accretion, which enhances the star formation at accretion. While, for SA satellites, the mergers occur smoothly and continuously across the accretion time. Moreover, mergers with FA satellites happen mainly after accretion, while a contrary trend is found for SA satellites. Our results provide insight into the evolution and star formation quenching of the satellite population.

preprint2020arXiv

The Formation of Ultra-Diffuse Galaxies in Clusters

We use the IllustrisTNG cosmological hydrodynamical simulation to study the formation of ultra-diffuse galaxies (UDGs) in galaxy clusters. We supplement the simulations with a realistic mass-size relation for galaxies at the time of infall into the cluster, as well as an analytical model to describe the tidally-induced evolution of their stellar mass, velocity dispersion and size. The model assumes "cuspy" NFW halos and, contrary to recent claims, has no difficulty reproducing the observed number of UDGs in clusters. Our results further suggest that the UDG population consists of a mixture of "normal" low surface brightness galaxies such as those found in the field ("born" UDGs, or B-UDGs), as well as a distinct population that owe their large size and low surface brightness to the effects of cluster tides ("tidal", or T-UDGs). The simulations indicate that T-UDGs entered the cluster earlier and should be more prevalent than B-UDGs near the cluster centres. T-UDGs should also have, at given stellar mass, lower velocity dispersion, higher metallicities, and lower dark matter content than B-UDGs. Our results suggest that systems like DF-44 are consistent with having been born as UDGs, while others such as DF2, DF4 and VLSB-D are possibly extreme T-UDG examples.

preprint2020arXiv

The relationship between fine galaxy stellar morphology and star formation activity in cosmological simulations: a deep learning view

Hydrodynamical simulations of galaxy formation and evolution attempt to fully model the physics that shapes galaxies. The agreement between the morphology of simulated and real galaxies, and the way the morphological types are distributed across galaxy scaling relations are important probes of our knowledge of galaxy formation physics. Here we propose an unsupervised deep learning approach to perform a stringent test of the fine morphological structure of galaxies coming from the Illustris and IllustrisTNG (TNG100 and TNG50) simulations against observations from a subsample of the Sloan Digital Sky Survey. Our framework is based on PixelCNN, an autoregressive model for image generation with an explicit likelihood. We adopt a strategy that combines the output of two PixelCNN networks in a metric that isolates the fine morphological details of galaxies from the sky background. We are able to \emph{quantitatively} identify the improvements of IllustrisTNG, particularly in the high-resolution TNG50 run, over the original Illustris. However, we find that the fine details of galaxy structure are still different between observed and simulated galaxies. This difference is driven by small, more spheroidal, and quenched galaxies which are globally less accurate regardless of resolution and which have experienced little improvement between the three simulations explored. We speculate that this disagreement, that is less severe for quenched disky galaxies, may stem from a still too coarse numerical resolution, which struggles to properly capture the inner, dense regions of quenched spheroidal galaxies.

preprint2020arXiv

Tidally induced warps of spiral galaxies in IllustrisTNG

Warps are common features in both stellar and gaseous disks of nearby spiral galaxies with the latter usually easier to detect. Several theories have been proposed in the literature to explain their formation and prevalence, including tidal interactions with external galaxies. Observational correlations also suggest the importance of tides for warp formation. Here, we use the TNG100 run from the magnetohydrodynamical cosmological simulation suite IllustrisTNG to investigate the connection between interactions and the formation of gas warps. We find that in the sample of well-resolved gas-rich spiral galaxies ($10^{10}\lesssim\mathrm{M_{*}/M_{\odot}}\lesssim10^{11}$ at $z=0$) from the simulation TNG100-1, about $16\%$ possess the characteristic S-shaped warp. Around one third of these objects have their vertical morphology induced by interactions with other galaxies. Half of these interactions end with the perturber absorbed by the host. Warps formed in interactions are more asymmetrical than the remaining sample, however after the interaction the asymmetry decreases with time. We find that warps induced by interactions survive on average for $<1$ Gyr. The angle between the orbital angular momentum of the perturber and the angular momentum of the host's disk that most likely leads to warp formation is around 45 degrees. While our main goal is to investigate tidally induced warps, we find that during interactions in addition to tides, new gas that is accreted from infalling satellites also can contribute to warp formation.

preprint2020arXiv

X-ray Signatures of Black Hole Feedback: Hot Galactic Atmospheres in IllustrisTNG and X-ray Observations

Hot gaseous atmospheres that permeate galaxies and extend far beyond their stellar distribution, where they are commonly referred to as the circumgalactic medium (CGM), imprint important information about feedback processes powered by the stellar populations of galaxies and their central supermassive black holes (SMBH). In this work we study the properties of this hot X-ray emitting medium using the IllustrisTNG cosmological simulations. We analyse their mock X-ray spectra, obtained from the diffuse and metal-enriched gas in TNG100 and TNG50, and compare the results with X-ray observations of nearby early-type galaxies. The simulations reproduce the observed X-ray luminosities ($L_{\rm X}$) and temperature ($T_{\rm X})$ at small ($<R_{\rm e}$) and intermediate ($<5R_{\rm e}$) radii reasonably well. We find that the X-ray properties of lower mass galaxies depend on their star formation rates. In particular, in the magnitude range where the star-forming and quenched populations overlap, $M_{\rm K}\sim-24$ $ (M_*\sim10^{10.7}M_\odot)$, we find that the X-ray luminosities of star-forming galaxies are on average about an order of magnitude higher than those of their quenched counterparts. We show that this diversity in $L_{\rm X}$ is a direct manifestation of the quenching mechanism in the simulations, where the galaxies are quenched due to gas expulsion driven by SMBH kinetic feedback. The observed dichotomy in $L_{\rm X}$ is thus an important observable prediction for the SMBH feedback-based quenching mechanisms implemented in state-of-the-art cosmological simulations. While the current X-ray observations of star forming galaxies are broadly consistent with the predictions of the simulations, the observed samples are small and more decisive tests are expected from the sensitive all-sky X-ray survey with eROSITA.

preprint2019arXiv

Antlia2's role in driving the ripples in the outer gas disk of the Galaxy

We employ the earlier published proper motions of the newly discovered Antlia 2 dwarf galaxy derived from Gaia data to calculate its orbital distribution in the cosmologically recent past. Using these observationally motivated orbits, we calculate the effect of the Antlia 2 dwarf galaxy on the outer HI disk of the Milky Way, using both test particle and Smoothed Particle Hydrodynamics simulations. We find that orbits with low pericenters, $\sim$ 10 kpc, produce disturbances that match the observed outer HI disk perturbations. We have independently recalculated the proper motion of the Antlia 2 dwarf from Gaia data and found a proper motion of $(μ_αcosδ, μ_δ) = (-0.068,0.032) \pm (0.023,-0.031)~\rm mas/yr$, which agrees with results from Torrealba et al. (2019) within the errors, but gives lower mean pericenters, e.g., $\sim$ 15 kpc for our fiducial model of the Milky Way. We also show that the Sagittarius dwarf galaxy interaction does not match the observed perturbations in the outer gas disk. Thus, Antlia 2 may be the driver of the observed large perturbations in the outer gas disk of the Galaxy. The current location of the Antlia 2 dwarf galaxy closely matches that predicted by an earlier dynamical analysis (Chakrabarti \& Blitz 2009) of the dwarf that drove ripples in the outer Galaxy, and, in particular, its orbit is nearly coplanar to the Galactic disk. If the Antlia 2 dwarf galaxy is responsible for the perturbations in the outer Galactic disk, it would have a specific range of proper motions that we predict here; this can be tested soon with Gaia DR-3 and Gaia DR-4 data.

preprint2019arXiv

Baryons in the CosmicWeb of IllustrisTNG -- II: the Connection among Galaxies, Halos, their Formation Time and their Location in the Cosmic Web

The connections among galaxies, the dark matter halos where they form and the properties of the large-scale Cosmic Web still need to be completely disentangled. We use the cosmological hydrodynamical simulation TNG100 of the IllustrisTNG suite to quantify the effects played by the large-scale density field and the Cosmic Web morphology on the relation between halo mass and galaxy stellar mass. We select objects with total dynamical mass in the range $\geq 6.3\times 10^{10} h ^{-1}\, M_{\odot}$ up to a few $10^{14} h^{-1} \, M_{\odot}$ between redshift $z=4$ and redshift $z=0$. A Cosmic Web class (knot, filament, sheet, void) is assigned to each region of the volume using a density field deformation tensor-based method. We find that galaxy stellar mass strongly correlates with total dynamical mass and formation time, and more weakly with large-scale overdensity and Cosmic Web class. The latter two quantities correlate with each other, but are not entirely degenerate. Furthermore, we find that at fixed halo mass, galaxies with stellar mass lower than the median value are more likely to be found in voids and sheets, whereas galaxies with stellar mass higher than the median are more likely to be found in filaments and knots. Finally, we find that the dependence on environment is stronger for satellites than for centrals, and discuss the physical implications of these results.

preprint2019arXiv

Correlations Between Black Holes and Host Galaxies in the Illustris and IllustrisTNG Simulations

We study black hole - host galaxy correlations, and the relation between the over-massiveness (the distance from the average $M_{BH}-σ$ relation) of super-massive black holes (SMBHs) and star formation histories of their host galaxies in the Illustris and TNG100 simulations. We find that both simulations are able to produce black hole scaling relations in general agreement with observations at $z=0$, but with noticeable discrepancies. Both simulations show an offset from the observations for the $M_{BH}-σ$ relation, and the relation between $M_{BH}$ and the Sersic index. The relation between $M_{BH}$ and stellar mass $M_*$ is tighter than the observations, especially for TNG100. For massive galaxies in both simulations, the hosts of over-massive SMBHs (those above the mean $M_{BH}-σ$ relation) tend to have larger Sersic indices and lower baryon conversion efficiency, suggesting a multidimensional link between SMBHs and properties of their hosts. In Illustris, the hosts of over-massive SMBHs have formed earlier and have lower present-day star formation rates, in qualitative agreement with the observations for massive galaxies with $σ>100 \rm km/s$. For low-mass galaxies, such a correlation still holds in Illustris but does not exist in the observed data. For TNG100, the correlation between SMBH over-massiveness and star formation history is much weaker. The hosts of over-massive SMBHs generally have consistently larger star formation rates throughout history. These galaxies have higher stellar mass as well, due to the strong $M_{BH}-M_*$ correlation. Our findings show that simulated SMBH scaling relations and correlations are sensitive to features in the modeling of SMBHs.

preprint2019arXiv

Galaxy Formation with BECDM -- II. Cosmic Filaments and First Galaxies

Bose-Einstein Condensate Dark Matter (BECDM; also known as Fuzzy Dark Matter) is motivated by fundamental physics and has recently received significant attention as a serious alternative to the established Cold Dark Matter (CDM) model. We perform cosmological simulations of BECDM gravitationally coupled to baryons and investigate structure formation at high redshifts ($z \gtrsim 5$) for a boson mass $m=2.5\cdot 10^{-22}~{\rm eV}$, exploring the dynamical effects of its wavelike nature on the cosmic web and the formation of first galaxies. Our BECDM simulations are directly compared to CDM as well as to simulations where the dynamical quantum potential is ignored and only the initial suppression of the power spectrum is considered -- a Warm Dark Matter-like ("WDM") model often used as a proxy for BECDM. Our simulations confirm that "WDM" is a good approximation to BECDM on large cosmological scales even in the presence of the baryonic feedback. Similarities also exist on small scales, with primordial star formation happening both in isolated haloes and continuously along cosmic filaments; the latter effect is not present in CDM. Global star formation and metal enrichment in these first galaxies are delayed in BECDM/"WDM" compared to the CDM case: in BECDM/"WDM" first stars form at $z\sim 13$/$13.5$ while in CDM star formation starts at $z\sim 35$. The signature of BECDM interference, not present in "WDM", is seen in the evolved dark matter power spectrum: although the small scale structure is initially suppressed, power on kpc scales is added at lower redshifts. Our simulations lay the groundwork for realistic simulations of galaxy formation in BECDM.

preprint2019arXiv

Galaxy interactions in IllustrisTNG-100, I: The power and limitations of visual identification

We present a sample of 446 galaxy pairs constructed using the cosmological simulation IllustrisTNG-100 at z = 0, with M$_{FoF, dm}$ = 10$^{11}$-10$^{13.5}$ M$_{\odot}$. We produce ideal mock SDSS g-band images of all pairs to test the reliability of visual classification schema employed to produce samples of interacting galaxies. We visually classify each image as interacting or not based on the presence of a close neighbour, the presence of stellar debris fields, disturbed discs, and/or tidal features. By inspecting the trajectories of the pairs, we determine that these indicators correctly identify interacting galaxies $\sim45\%$ of the time. We subsequently split the sample into the visually identified interacting pairs (VIP; 38 pairs) and those which are interacting but are not visually identified (nonVIP; 47 pairs). We find that VIP have undergone a close passage nearly twice as recently as the nonVIP, and typically have higher stellar masses. Further, the VIP sit in dark matter haloes that are approximately 2.5 times as massive, in environments nearly 2 times as dense, and are almost a factor of 10 more affected by the tidal forces of their surroundings than the nonVIP. These factors conspire to increase the observability of tidal features and disturbed morphologies, making the VIP more likely to be identified. Thus, merger rate calculations which rely on stellar morphologies are likely to be significantly biased toward massive galaxy pairs which have recently undergone a close passage.

preprint2019arXiv

High redshift JWST predictions from IllustrisTNG: Dust modelling and galaxy luminosity functions

The James Webb Space Telescop (JWST) promises to revolutionise our understanding of the early Universe, and contrasting its upcoming observations with predictions of the $Λ$CDM model requires detailed theoretical forecasts. Here, we exploit the large dynamic range of the IllustrisTNG simulation suite, TNG50, TNG100, and TNG300, to derive multi-band galaxy luminosity functions from $z=2$ to $z=10$. We put particular emphasis on the exploration of different dust attenuation models to determine galaxy luminosity functions for the rest-frame ultraviolet (UV), and apparent wide NIRCam bands. Our most detailed dust model is based on continuum Monte Carlo radiative transfer calculations employing observationally calibrated dust properties. This calibration results in constraints on the redshift evolution of the dust attenuation normalisation and dust-to-metal ratios yielding a stronger redshift evolution of the attenuation normalisation compared to most previous theoretical studies. Overall we find good agreement between the rest-frame UV luminosity functions and observational data for all redshifts, also beyond the regimes used for the dust-model calibrations. Furthermore, we also recover the observed high redshift ($z=4-6$) UV luminosity versus stellar mass relation, the H$α$ versus star formation rate relation, and the H$α$ luminosity function at $z=2$. The bright end ($M_{\rm UV}>-19.5$) cumulative galaxy number densities are consistent with observational data. For the F200W NIRCam band, we predict that JWST will detect $\sim 80$ ($\sim 200$) galaxies with a signal-to-noise ratio of $10$ ($\sim 5$) within the NIRCam field of view, $2.2\times2.2 \,{\rm arcmin}^{2}$, for a total exposure time of $10^5{\rm s}$ in the redshift range $z=8 \pm 0.5$. These numbers drop to $\sim 10$ ($\sim 40$) for an exposure time of $10^4{\rm s}$.

preprint2019arXiv

The relationship between black hole mass and galaxy properties: Examining the black hole feedback model in IllustrisTNG

Supermassive black hole feedback is thought to be responsible for the lack of star formation, or quiescence, in a significant fraction of galaxies. We explore how observable correlations between the specific star formation rate (sSFR), stellar mass (M$_{\rm{star}}$), and black hole mass (M$_{\rm{BH}}$) are sensitive to the physics of black hole feedback in a galaxy formation model. We use the IllustrisTNG simulation suite, specifically the TNG100 simulation and ten model variations that alter the parameters of the black hole model. Focusing on central galaxies at $z = 0$ with M$_{\rm{star}} > 10^{10}$ M$_{\odot}$, we find that the sSFR of galaxies in IllustrisTNG decreases once the energy from black hole kinetic winds at low accretion rates becomes larger than the gravitational binding energy of gas within the galaxy stellar radius. This occurs at a particular M$_{\rm{BH}}$ threshold above which galaxies are found to sharply transition from being mostly star-forming to mostly quiescent. As a result of this behavior, the fraction of quiescent galaxies as a function of M$_{\rm{star}}$ is sensitive to both the normalization of the M$_{\rm{BH}}$-M$_{\rm{star}}$ relation and the M$_{\rm{BH}}$ threshold for quiescence in IllustrisTNG. Finally, we compare these model results to observations of 91 central galaxies with dynamical M$_{\rm{BH}}$ measurements with the caveat that this sample is not representative of the whole galaxy population. While IllustrisTNG reproduces the observed trend that quiescent galaxies host more massive black holes, the observations exhibit a broader scatter in M$_{\rm{BH}}$ at a given M$_{\rm{star}}$ and show a smoother decline in sSFR with M$_{\rm{BH}}$.

preprint2018arXiv

Efficacy of early stellar feedback in low gas surface density environments

We present a suite of high resolution radiation hydrodynamic simulations of a small patch ($1 \ {\rm kpc}^2$) of the inter-stellar medium (ISM) performed with Arepo-RT, with the aim to quantify the efficacy of various feedback processes like supernovae explosions (SNe), photoheating and radiation pressure in low gas surface density galaxies ($Σ_{\rm gas} \simeq 10 \ {\rm M}_\odot \ {\rm pc}^{-2}$). We show that radiation fields decrease the star formation rate and therefore the total stellar mass formed by a factor of $\sim 2$. This increases the gas depletion timescale and brings the simulated Kennicutt-Schmidt relation closer to the observational estimates. Radiation feedback coupled with SNe is more efficient at driving outflows with the mass and energy loading increasing by a factor of $\sim 10$. This increase is mainly driven by the additional entrainment of medium density ($10^{-2} \leq n< 1 \ {\rm cm}^{-3}$), warm ($300 \ {\rm K}\leq T<8000 \ {\rm K}$) material. Therefore including radiation fields tends to launch colder, denser and higher mass and energy loaded outflows. This is because photoheating of the high density gas around a newly formed star over-pressurises the region, causing it to expand. This reduces the ambient density in which the SNe explode by a factor of $10-100$ which in turn increases their momentum output by a factor of $\sim 1.5-2.5$. Finally, we note that in these low gas surface density environments, radiation fields primarily impact the ISM via photoheating and radiation pressure has only a minimal role in regulating star formation.

preprint2016arXiv

A moving mesh unstaggered constrained transport scheme for magnetohydrodynamics

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

preprint2016arXiv

Discriminating Topology in Galaxy Distributions using Network Analysis

(abridged) The large-scale distribution of galaxies is generally analyzed using the two-point correlation function. However, this statistic does not capture the topology of the distribution, and it is necessary to resort to higher order correlations to break degeneracies. We demonstrate that an alternate approach using network analysis can discriminate between topologically different distributions that have similar two-point correlations. We investigate two galaxy point distributions, one produced by a cosmological simulation and the other by a Lévy walk. For the cosmological simulation, we adopt the redshift $z = 0.58$ slice from Illustris (Vogelsberger et al. 2014A) and select galaxies with stellar masses greater than $10^8$$M_\odot$. The two point correlation function of these simulated galaxies follows a single power-law, $ξ(r) \sim r^{-1.5}$. Then, we generate Lévy walks matching the correlation function and abundance with the simulated galaxies. We find that, while the two simulated galaxy point distributions have the same abundance and two point correlation function, their spatial distributions are very different; most prominently, \emph{filamentary structures}, absent in Lévy fractals. To quantify these missing topologies, we adopt network analysis tools and measure diameter, giant component, and transitivity from networks built by a conventional friends-of-friends recipe with various linking lengths. Unlike the abundance and two point correlation function, these network quantities reveal a clear separation between the two simulated distributions; therefore, the galaxy distribution simulated by Illustris is not a Lévy fractal quantitatively. We find that the described network quantities offer an efficient tool for discriminating topologies and for comparing observed and theoretical distributions.

preprint2016arXiv

Large-scale mass distribution in the Illustris simulation

Observations at low redshifts thus far fail to account for all of the baryons expected in the Universe according to cosmological constraints. A large fraction of the baryons presumably resides in a thin and warm-hot medium between the galaxies, where they are difficult to observe due to their low densities and high temperatures. Cosmological simulations of structure formation can be used to verify this picture and provide quantitative predictions for the distribution of mass in different large-scale structure components. Here we study the distribution of baryons and dark matter at different epochs using data from the Illustris simulation. We identify regions of different dark matter density with the primary constituents of large-scale structure, allowing us to measure mass and volume of haloes, filaments and voids. At redshift zero, we find that 49 % of the dark matter and 23 % of the baryons are within haloes more massive than the resolution limit of $2\times 10^8$ M$_\odot$. The filaments of the cosmic web host a further 45 % of the dark matter and 46 % of the baryons. The remaining 31 % of the baryons reside in voids. The majority of these baryons have been transported there through active galactic nuclei feedback. We note that the feedback model of Illustris is too strong for heavy haloes, therefore it is likely that we are overestimating this amount. Categorizing the baryons according to their density and temperature, we find that 17.8 % of them are in a condensed state, 21.6 % are present as cold, diffuse gas, and 53.9 % are found in the state of a warm-hot intergalactic medium.

preprint2016arXiv

Massive Black Hole Binary Mergers in Dynamical Galactic Environments

Gravitational Waves (GW) have now been detected from stellar-mass black hole binaries, and the first observations of GW from Massive Black Hole (MBH) Binaries are expected within the next decade. Pulsar Timing Arrays (PTA), which can measure the years long periods of GW from MBHB, have excluded many standard predictions for the amplitude of a stochastic GW Background (GWB). We use coevolved populations of MBH and galaxies from hydrodynamic, cosmological simulations ('Illustris') to calculate a predicted GWB. The most advanced predictions so far have included binary hardening mechanisms from individual environmental processes. We present the first calculation including all of the environmental mechanisms expected to be involved: dynamical friction, stellar 'loss-cone' scattering, and viscous drag from a circumbinary disk. We find that MBH binary lifetimes are generally multiple gigayears, and only a fraction coalesce by redshift zero. For a variety of parameters, we find all GWB amplitudes to be below the most stringent PTA upper limit of $A_{\textrm{yr}^{-1}} \approx 10^{-15}$. Our fairly conservative fiducial model predicts an amplitude of $A_{\textrm{yr}^{-1}} \approx 0.4\times 10^{-15}$---less than a factor of three below the current limit. At lower frequencies, we find $A_{0.1\,\textrm{yr}^{-1}} \approx 1.5\times 10^{-15}$ with spectral indices between $-0.4$ and $-0.6$---significantly flatter than the canonical value of $-2/3$ due to purely GW-driven evolution. Typical MBHB driving the GWB signal come from redshifts around $0.3$, with total masses of a few times $10^9\,M_\odot$, and in host galaxies with very large stellar masses. Even without GWB detections, our results can be connected to observations of dual AGN to constrain binary evolution.

preprint2016arXiv

Shock finding on a moving-mesh: II. Hydrodynamic shocks in the Illustris universe

Hydrodynamical shocks are a manifestation of the non-linearity of the Euler equations and play a fundamental role in cosmological gas dynamics. In this work, we identify and analyse shocks in the Illustris simulation, and contrast the results with those of non-radiative runs. We show that simulations with more comprehensive physical models of galaxy formation pose new challenges for shock finding algorithms due to radiative cooling and star-forming processes, prompting us to develop a number of methodology improvements. We find in Illustris a total shock surface area which is about 1.4 times larger at the present epoch compared to non-radiative runs, and an energy dissipation rate at shocks which is higher by a factor of around 7. Remarkably, shocks with Mach numbers above and below $\mathcal{M}\approx10$ contribute about equally to the total dissipation across cosmic time. This is in sharp contrast to non-radiative simulations, and we demonstrate that a large part of the difference arises due to strong black hole radio-mode feedback in Illustris. We also provide an overview of the large diversity of shock morphologies, which includes complex networks of halo-internal shocks, shocks on to cosmic sheets, feedback shocks due to black holes and galactic winds, as well as ubiquitous accretion shocks. In high redshift systems more massive than $10^{12}\,\mathrm{M}_\odot$ we discover the existence of a double accretion shock pattern in haloes. They are created when gas streams along filaments without being shocked at the outer accretion shock, but then forms a second, roughly spherical accretion shock further inside.

preprint2016arXiv

Simulating the Carbon Footprint of Galactic Halos

We compare simulations, including the Illustris simulations, to observations of CIV and CII absorption at z=2-4. These are the CIV column density distribution function in the column density range $10^{12} - 10^{15}$ cm$^{-2}$, the CIV equivalent width distribution at 0.1 - 2 Å, and the covering fractions and equivalent widths of CIV 1548 and CII 1337 around DLAs. In the context of the feedback models we investigate, all CIV observations favour the use of more energetic wind models, which are better able to enrich the gas surrounding halos. We propose two ways to achieve this; an increased wind velocity and an increase in wind thermal energy. However, even our most energetic wind models do not produce enough absorbers with CIV equivalent width > 0.6 Å, which in our simulations are associated with the most massive haloes. All simulations are in reasonable agreement with the CII covering fraction and equivalent widths around Damped Lyman-$α$ absorbers, although there is a moderate deficit in one bin 10 - 100 kpc from the DLA. Finally, we show that the CIV in our simulations is predominantly photoionized.

preprint2016arXiv

The Mass Profile of the Milky Way to the Virial Radius from the Illustris Simulation

We use particle data from the Illustris simulation, combined with individual kinematic constraints on the mass of the Milky Way (MW) at specific distances from the Galactic center, to infer the radial distribution of the MW's dark matter halo mass. Our method allows us to convert any constraint on the mass of the MW within a fixed distance to a full circular velocity profile to the MW's virial radius. As primary examples, we take two recent (and discrepant) measurements of the total mass within 50 kpc of the Galaxy and find they imply very different mass profiles and stellar masses for the Galaxy. The dark-matter-only version of the Illustris simulation enables us to compute the effects of galaxy formation on such constraints on a halo-by-halo basis; on small scales, galaxy formation enhances the density relative to dark-matter-only runs, while the total mass density is approximately 20% lower at large Galactocentric distances. We are also able to quantify how current and future constraints on the mass of the MW within specific radii will be reflected in uncertainties on its virial mass: even a measurement of M(<50 kpc) with essentially perfect precision still results in a 20% uncertainty on the virial mass of the Galaxy, while a future measurement of M(<100 kpc) with 10% errors would result in the same level of uncertainty. We expect that our technique will become even more useful as (1) better kinematic constraints become available at larger distances and (2) cosmological simulations provide even more faithful representations of the observable Universe.

preprint2016arXiv

The Network Behind the Cosmic Web

The concept of the cosmic web, viewing the Universe as a set of discrete galaxies held together by gravity, is deeply engrained in cosmology. Yet, little is known about the most effective construction and the characteristics of the underlying network. Here we explore seven network construction algorithms that use various galaxy properties, from their location, to their size and relative velocity, to assign a network to galaxy distributions provided by both simulations and observations. We find that a model relying only on spatial proximity offers the best correlations between the physical characteristics of the connected galaxies. We show that the properties of the networks generated from simulations and observations are identical, unveiling a deep universality of the cosmic web.

preprint2016arXiv

The Role of Star-Formation and AGN in Dust Heating of z=0.3-2.8 Galaxies - II. Informing IR AGN fraction estimates through simulations

A key question in extragalactic studies is the determination of the relative roles of stars and AGN in powering dusty galaxies at $z\sim$1-3 where the bulk of star-formation and AGN activity took place. In Paper I, we present a sample of $336$ 24$μ$m-selected (Ultra)Luminous Infrared Galaxies, (U)LIRGs, at $z \sim 0.3$-$2.8$, where we focus on determining the AGN contribution to the IR luminosity. Here, we use hydrodynamic simulations with dust radiative transfer of isolated and merging galaxies, to investigate how well the simulations reproduce our empirical IR AGN fraction estimates and determine how IR AGN fractions relate to the UV-mm AGN fraction. We find that: 1) IR AGN fraction estimates based on simulations are in qualitative agreement with the empirical values when host reprocessing of the AGN light is considered; 2) for star-forming galaxy-AGN composites our empirical methods may be underestimating the role of AGN, as our simulations imply $>$50% AGN fractions, $\sim$3$\times$ higher than previous estimates; 3) 6% of our empirically classified "SFG" have AGN fractions $\gtrsim$50%. While this is a small percentage of SFGs, if confirmed, would imply the true number density of AGN may be underestimated; 4) this comparison depends on the adopted AGN template -- those that neglect the contribution of warm dust lower the empirical fractions by up to 2$\times$; and 5) the IR AGN fraction is only a good proxy for the intrinsic UV-mm AGN fraction when the extinction is high ($A_V\gtrsim 1$ or up to and including coalescence in a merger).

preprint2015arXiv

An analysis of the evolving comoving number density of galaxies in hydrodynamical simulations

The cumulative comoving number-density of galaxies as a function of stellar mass or central velocity dispersion is commonly used to link galaxy populations across different epochs. By assuming that galaxies preserve their number-density in time, one can infer the evolution of their properties, such as masses, sizes, and morphologies. However, this assumption does not hold in the presence of galaxy mergers or when rank ordering is broken owing to variable stellar growth rates. We present an analysis of the evolving comoving number density of galaxy populations found in the Illustris cosmological hydrodynamical simulation focused on the redshift range $0\leq z \leq 3$. Our primary results are as follows: 1) The inferred average stellar mass evolution obtained via a constant comoving number density assumption is systematically biased compared to the merger tree results at the factor of $\sim$2(4) level when tracking galaxies from redshift $z=0$ out to redshift $z=2(3)$; 2) The median number density evolution for galaxy populations tracked forward in time is shallower than for galaxy populations tracked backward in time; 3) A similar evolution in the median number density of tracked galaxy populations is found regardless of whether number density is assigned via stellar mass, stellar velocity dispersion, or dark matter halo mass; 4) Explicit tracking reveals a large diversity in galaxies' assembly histories that cannot be captured by constant number-density analyses; 5) The significant scatter in galaxy linking methods is only marginally reduced by considering a number of additional physical and observable galaxy properties as realized in our simulation. We provide fits for the forward and backward median evolution in stellar mass and number density and discuss implications of our analysis for interpreting multi-epoch galaxy property observations.

preprint2015arXiv

Galactic Angular Momentum in the Illustris Simulation: Feedback and the Hubble Sequence

We study the stellar angular momentum of thousands of galaxies in the Illustris cosmological simulation, which captures gravitational and gas dynamics within galaxies, as well as feedback from stars and black holes. We find that the angular momentum of the simulated galaxies matches observations well, and in particular two distinct relations are found for late-type versus early-type galaxies. The relation for late-type galaxies corresponds to the value expected from full conservation of the specific angular momentum generated by cosmological tidal torques. The relation for early-type galaxies corresponds to retention of only ~30% of that, but we find that those early-type galaxies with low angular momentum at z=0 nevertheless reside at high redshift on the late-type relation. Some of them abruptly lose angular momentum during major mergers. To gain further insight, we explore the scaling relations in simulations where the galaxy formation physics is modified with respect to the fiducial model. We find that galactic winds with high mass-loading factors are essential for obtaining the high angular momentum relation typical for late-type galaxies, while AGN feedback largely operates in the opposite direction. Hence, feedback controls the stellar angular momentum of galaxies, and appears to be instrumental for establishing the Hubble sequence.

preprint2015arXiv

Galaxy Morphology and Star Formation in the Illustris Simulation at z=0

We study how optical galaxy morphology depends on mass and star formation rate (SFR) in the Illustris Simulation. To do so, we measure automated galaxy structures in 10808 simulated galaxies at z=0 with stellar masses 10^9.7 < M_*/M_sun < 10^12.3. We add observational realism to idealized synthetic images and measure non-parametric statistics in rest-frame optical and near-IR images from four directions. We find that Illustris creates a morphologically diverse galaxy population, occupying the observed bulge strength locus and reproducing median morphology trends versus stellar mass, SFR, and compactness. Morphology correlates realistically with rotation, following classification schemes put forth by kinematic surveys. Type fractions as a function of environment agree roughly with data. These results imply that connections among mass, star formation, and galaxy structure arise naturally from models matching global star formation and halo occupation functions when simulated with accurate methods. This raises a question of how to construct experiments on galaxy surveys to better distinguish between models. We predict that at fixed halo mass near 10^12 M_sun, disc-dominated galaxies have higher stellar mass than bulge-dominated ones, a possible consequence of the Illustris feedback model. While Illustris galaxies at M_* ~ 10^11 M_sun have a reasonable size distribution, those at M_* ~ 10^10 M_sun have half-light radii larger than observed by a factor of two. Furthermore, at M_* ~ 10^10.5-10^11 M_sun, a relevant fraction of Illustris galaxies have distinct "ring-like" features, such that the bright pixels have an unusually wide spatial extent.

preprint2015arXiv

Hot Gaseous Coronae around Spiral Galaxies: Probing the Illustris Simulation

The presence of hot gaseous coronae around present-day massive spiral galaxies is a fundamental prediction of galaxy formation models. However, our observational knowledge remains scarce, since to date only four gaseous coronae were detected around spirals with massive stellar bodies ($\gtrsim2\times10^{11} \ \rm{M_{\odot}}$). To explore the hot coronae around lower mass spiral galaxies, we utilized Chandra X-ray observations of a sample of eight normal spiral galaxies with stellar masses of $(0.7-2.0)\times10^{11} \ \rm{M_{\odot}}$. Although statistically significant diffuse X-ray emission is not detected beyond the optical radii ($\sim20$ kpc) of the galaxies, we derive $3σ$ limits on the characteristics of the coronae. These limits, complemented with previous detections of NGC 1961 and NGC 6753, are used to probe the Illustris Simulation. The observed $3σ$ upper limits on the X-ray luminosities and gas masses exceed or are at the upper end of the model predictions. For NGC 1961 and NGC 6753 the observed gas temperatures, metal abundances, and electron density profiles broadly agree with those predicted by Illustris. These results hint that the physics modules of Illustris are broadly consistent with the observed properties of hot coronae around spiral galaxies. However, a shortcoming of Illustris is that massive black holes, mostly residing in giant ellipticals, give rise to powerful radio-mode AGN feedback, which results in under luminous coronae for ellipticals.

preprint2015arXiv

Hydrogen Reionization in the Illustris universe

Hydrodynamical simulations of galaxy formation such as the Illustris simulations have progressed to a state where they approximately reproduce the observed stellar mass function from high to low redshift. This in principle allows self-consistent models of reionization that exploit the accurate representation of the diffuse gas distribution together with the realistic growth of galaxies provided by these simulations, within a representative cosmological volume. In this work, we apply and compare two radiative transfer algorithms implemented in a GPU-accelerated code to the $106.5\,{\rm Mpc}$ wide volume of Illustris in postprocessing in order to investigate the reionization transition predicted by this model. We find that the first generation of galaxies formed by Illustris is just about able to reionize the universe by redshift $z\sim 7$, provided quite optimistic assumptions about the escape fraction and the resolution limitations are made. Our most optimistic model finds an optical depth of $τ\simeq 0.065$, which is in very good agreement with recent Planck 2015 determinations. Furthermore, we show that moment-based approaches for radiative transfer with the M1 closure give broadly consistent results with our angular-resolved radiative transfer scheme. In our favoured fiducial model, 20\% of the hydrogen is reionized by redshift $z=9.20$, and this rapidly climbs to 80\% by redshift $z=6.92$. It then takes until $z=6.24$ before 99\% of the hydrogen is ionized. On average, reionization proceeds `inside-out' in our models, with a size distribution of reionized bubbles that progressively features regions of ever larger size while the abundance of small bubbles stays fairly constant.

preprint2015arXiv

Mapping galaxy encounters in numerical simulations: The spatial extent of induced star formation

We employ a suite of 75 simulations of galaxies in idealised major mergers (stellar mass ratio ~2.5:1), with a wide range of orbital parameters, to investigate the spatial extent of interaction-induced star formation. Although the total star formation in galaxy encounters is generally elevated relative to isolated galaxies, we find that this elevation is a combination of intense enhancements within the central kpc and moderately suppressed activity at large galacto-centric radii. The radial dependence of the star formation enhancement is stronger in the less massive galaxy than in the primary, and is also more pronounced in mergers of more closely aligned disc spin orientations. Conversely, these trends are almost entirely independent of the encounter's impact parameter and orbital eccentricity. Our predictions of the radial dependence of triggered star formation, and specifically the suppression of star formation beyond kph-scales, will be testable with the next generation of integral-field spectroscopic surveys.

preprint2015arXiv

Modeling Galactic Conformity with the Color-Halo Age Relation in the Illustris Simulation

Comparisons between observational surveys and galaxy formation models find that the mass of dark matter haloes can largely explain galaxies' stellar mass. However, it remains uncertain whether additional environmental variables, generally referred to 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 $M_{stellar}$ > $2 \times 10^9$ $M_{\odot}$. 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 can be given as a combination of the halo conformity signal with the galaxy color-halo age relation: at fixed stellar mass, particularly toward the low-mass end, Illustris' galaxy colors correlate with halo age, with the reddest galaxies (often satellites) being preferentially found in the oldest haloes. In fact, we can explain the galactic conformity effect with a simple semi-empirical model, by assigning stellar mass based on halo mass (abundance matching) and by assigning galaxy color based on halo age (age matching). We investigate other interpretations for the galactic conformity, particularly its dependence on the isolation criterion and on the central-satellite information. 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.

preprint2015arXiv

On the Assembly of Dwarf Galaxies in Clusters and their Efficient Formation of Globular Clusters

Galaxy clusters contain a large population of low mass dwarf elliptical galaxies whose exact origin is unclear: their colors, structural properties and kinematics differ substantially from those of dwarf irregulars in the field. We use the Illustris cosmological simulation to study differences in the assembly paths of dwarf galaxies (3e8 < M_*/M_sun < 1e10) according to their environment. We find that cluster dwarfs achieve their maximum total and stellar mass on average ~ 8 and ~ 4.5 Gyr ago (or redshifts z = 1.0 and z = 0.4, respectively), around the time of infall into the clusters. In contrast, field dwarfs not subjected to environmental stripping, reach their maximum mass at redshift z = 0. This different assembly history naturally produces a color bimodality, with blue isolated dwarfs and redder cluster dwarfs exhibiting negligible star-formation today. The cessation of star formation happens over median times 3.5-5 Gyr depending on stellar mass, and shows a large scatter (~ 1-8 Gyr), with the lower values associated with starburst events that occur at infall through the virial radius or pericentric passages. We argue that such starbursts together with the early assembly of cluster dwarfs can provide a natural explanation for the higher specific frequency of globular clusters (GCs) in cluster dwarfs, as found observationally. We present a simple model for the formation and stripping of GCs that supports this interpretation. The origin of dwarf ellipticals in clusters is, therefore, consistent with an environmentally-driven evolution of field dwarf irregulars. However the z = 0 field analogs of cluster dwarf progenitors have today stellar masses a factor ~ 3 larger --a difference arising from the early truncation of star formation in cluster dwarfs.

preprint2015arXiv

Planet-disc interaction on a freely moving mesh

General-purpose, moving-mesh schemes for hydrodynamics have opened the possibility of combining the accuracy of grid-based numerical methods with the flexibility and automatic resolution adaptivity of particle-based methods. Due to their supersonic nature, Keplerian accretion discs are in principle a very attractive system for applying such freely moving mesh techniques. However, the high degree of symmetry of simple accretion disc models can be difficult to capture accurately by these methods, due to the generation of geometric grid noise and associated numerical diffusion, which is absent in polar grids. To explore these and other issues, in this work we study the idealized problem of two-dimensional planet-disc interaction with the moving-mesh code AREPO. We explore the hydrodynamic evolution of discs with planets through a series of numerical experiments that vary the planet mass, the disc viscosity and the mesh resolution, and compare the resulting surface density, vortensity field and tidal torque with results from the literature. We find that the performance of the moving-mesh code in this problem is in accordance with published results, showing good consistency with grid codes written in polar coordinates. We also conclude that grid noise and mesh distortions do not introduce excessive numerical diffusion. Finally, we show how the moving-mesh approach can naturally increase resolution in regions of high densityaround planets and planetary wakes, while retaining the background flow at low resolution. This provides an alternative to the difficult task of implementing adaptive mesh refinement in conventional polar-coordinate codes.

preprint2015arXiv

Recoiling black holes: prospects for detection and implications of spin alignment

Supermassive black hole (BH) mergers produce powerful gravitational wave (GW) 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, less than about 10 spatially-offset AGN may be detectable in HST-COSMOS, and > 10^3 could be found with Pan-STARRS, LSST, Euclid, and WFIRST. Nearly a thousand velocity-offset AGN are predicted within the 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. Nonetheless, 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.

preprint2015arXiv

Reducing noise in moving-grid codes with strongly-centroidal Lloyd mesh regularization

A method for improving the accuracy of hydrodynamical codes that use a moving Voronoi mesh is described. Our scheme is based on a new regularization scheme that constrains the mesh to be centroidal to high precision while still allowing the cells to move approximately with the local fluid velocity, thereby retaining the quasi-Lagrangian nature of the approach. Our regularization technique significantly reduces mesh noise that is attributed to changes in mesh topology and deviations from mesh regularity. We demonstrate the advantages of our method on various test problems, and note in particular improvements obtained in handling shear instabilities, mixing, and in angular momentum conservation. Calculations of adiabatic jets in which shear excites Kelvin Helmholtz instability show reduction of mesh noise and entropy generation. In contrast, simulations of the collapse and formation of an isolated disc galaxy are nearly unaffected, showing that numerical errors due to the choice of regularization do not impact the outcome in this case.

preprint2015arXiv

Reproducing the Kinematics of Damped Lyman-alpha Systems

We examine the kinematic structure of Damped Lyman-alpha Systems (DLAs) in a series of cosmological hydrodynamic simulations using the AREPO code. We are able to match the distribution of velocity widths of associated low ionisation metal absorbers substantially better than earlier work. Our simulations produce a population of DLAs dominated by halos with virial velocities around 70 km/s, consistent with a picture of relatively small, faint objects. In addition, we reproduce the observed correlation between velocity width and metallicity and the equivalent width distribution of SiII. Some discrepancies of moderate statistical significance remain; too many of our spectra show absorption concentrated at the edge of the profile and there are slight differences in the exact shape of the velocity width distribution. We show that the improvement over previous work is mostly due to our strong feedback from star formation and our detailed modelling of the metal ionisation state.

preprint2015arXiv

Stellar orbit evolution in close circumstellar disc encounters

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 3-D, focusing on disc-disc interactions that result from stellar fly-bys. 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 3-D models of self-gravitating discs, and then create a suite of numerical experiments of parabolic encounters, exploring the effects of pericenter separation r_p, disc orientation and disc-star mass ratio (M_d/M_*) on the orbital evolution of the host stars. Close encounters (2r_p<~ 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 M_d<~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.

preprint2015arXiv

The diverse evolutionary paths of simulated high-z massive, compact galaxies to z=0

Massive quiescent galaxies have much smaller physical sizes at high redshift than today. The strong evolution of galaxy size may be caused by progenitor bias, major and minor mergers, adiabatic expansion, and/or renewed star formation, but it is difficult to test these theories observationally. Herein, we select a sample of 35 massive, compact galaxies ($M_* = 1-3 \times 10^{11}$ M$_\odot$, $M_*/R^{1.5} > 10^{10.5}$ M$_\odot$/kpc$^{1.5}$) at $z=2$ in the cosmological hydrodynamical simulation Illustris and trace them forward to $z=0$ to uncover their evolution and identify their descendants. By $z=0$, the original factor of 3 difference in stellar mass spreads to a factor of 20. The dark matter halo masses similarly spread from a factor of 5 to 40. The galaxies' evolutionary paths are diverse: about half acquire an ex-situ envelope and are the core of a more massive descendant, a third survive undisturbed and gain very little mass, 15% are consumed in a merger with a more massive galaxy, and a small remainder are thoroughly mixed by major mergers. The galaxies grow in size as well as mass, and only $\sim$10% remain compact by $z=0$. The majority of the size growth is driven by the acquisition of ex-situ mass. The most massive galaxies at $z=0$ are the most likely to have compact progenitors, but this trend possesses significant dispersion which precludes a direct linkage to compact galaxies at $z=2$. The compact galaxies' merger rates are influenced by their $z=2$ environments, so that isolated or satellite compact galaxies (which are protected from mergers) are the most likely to survive to the present day.

preprint2015arXiv

The Formation of Massive, Compact Galaxies at z=2 in the Illustris Simulation

Massive, quiescent galaxies at high redshift have been found to be considerably more compact than galaxies of similar mass in the local universe. How these compact galaxies formed has yet to be determined, though several progenitor populations have been proposed. Here we investigate the formation processes and quantify the assembly histories of such galaxies in Illustris, a suite of hydrodynamical cosmological simulations encompassing a sufficiently large volume to include rare objects, while simultaneously resolving the internal structure of galaxies. We select massive (~10^11 solar masses) and compact (stellar half-mass radius < 2 kpc) galaxies from the simulation at z=2. Within the Illustris suite, we find that these quantities are not perfectly converged, but are reasonably reliable for our purposes. The resulting population is composed primarily of quiescent galaxies, but we also find several star-forming compact galaxies. The simulated compact galaxies are similar to observed galaxies in star formation activity and appearance. We follow their evolution at high redshift in the simulation and find that there are multiple pathways to form these compact galaxies, dominated by two mechanisms: (i) intense, centrally concentrated starbursts generally triggered by gas-rich major mergers between z~2-4, reducing the galaxies' half-mass radii by a factor of a few to below 2 kpc, and (ii) assembly at very early times when the universe was much denser; the galaxies formed compact and remained so until z~2.

preprint2015arXiv

The Illustris Simulation: Public Data Release

We present the full public release of all data from the Illustris simulation project. Illustris is a suite of large volume, cosmological hydrodynamical simulations run with the moving-mesh code Arepo and including a comprehensive set of physical models critical for following the formation and evolution of galaxies across cosmic time. Each simulates a volume of (106.5 Mpc)^3 and self-consistently evolves five different types of resolution elements from a starting redshift of z=127 to the present day, z=0. These components are: dark matter particles, gas cells, passive gas tracers, stars and stellar wind particles, and supermassive black holes. This data release includes the snapshots at all 136 available redshifts, halo and subhalo catalogs at each snapshot, and two distinct merger trees. Six primary realizations of the Illustris volume are released, including the flagship Illustris-1 run. These include three resolution levels with the fiducial "full" baryonic physics model, and a dark matter only analog for each. In addition, we provide four distinct, high time resolution, smaller volume "subboxes". The total data volume is ~265 TB, including ~800 full volume snapshots and ~30,000 subbox snapshots. We describe the released data products as well as tools we have developed for their analysis. All data may be directly downloaded in its native HDF5 format. Additionally, we release a comprehensive, web-based API which allows programmatic access to search and data processing tasks. In both cases we provide example scripts and a getting-started guide in several languages: currently, IDL, Python, and Matlab. This paper addresses scientific issues relevant for the interpretation of the simulations, serves as a pointer to published and on-line documentation of the project, describes planned future additional data releases, and discusses technical aspects of the release.

preprint2015arXiv

The Illustris simulation: the evolving population of black holes across cosmic time

We study the properties of black holes and their host galaxies across cosmic time in the Illustris simulation. Illustris is a large scale cosmological hydrodynamical simulation which resolves a (106.5 Mpc)^3 volume with more than 12 billion resolution elements and includes state-of-the-art physical models relevant for galaxy formation. We find that the black hole mass density for redshifts z = 0 - 5 and the black hole mass function at z = 0 predicted by Illustris are in very good agreement with the most recent observational constraints. We show that the bolometric and hard X-ray luminosity functions of AGN at z = 0 and 1 reproduce observational data very well over the full dynamic range probed. Unless the bolometric corrections are largely underestimated, this requires radiative efficiencies to be on average low, epsilon_r <= 0.1, noting however that in our model radiative efficiencies are degenerate with black hole feedback efficiencies. Cosmic downsizing of the AGN population is in broad agreement with the findings from X-ray surveys, but we predict a larger number density of faint AGN at high redshifts than currently inferred. We also study black hole -- host galaxy scaling relations as a function of galaxy morphology, colour and specific star formation rate. We find that black holes and galaxies co-evolve at the massive end, but for low mass, blue and star-forming galaxies there is no tight relation with either their central black hole masses or the nuclear AGN activity.

preprint2015arXiv

The impact of feedback on cosmological gas accretion

We investigate how the way galaxies acquire their gas across cosmic time in cosmological hydrodynamic simulations is modified by a comprehensive physical model for baryonic feedback processes. To do so, we compare two simulations -- with and without feedback -- both evolved with the moving mesh code AREPO. The feedback runs implement the full physics model of the Illustris simulation project, including star formation driven galactic winds and energetic feedback from supermassive blackholes. We explore: (a) the accretion rate of material contributing to the net growth of galaxies and originating directly from the intergalactic medium, finding that feedback strongly suppresses the raw, as well as the net, inflow of this "smooth mode" gas at all redshifts, regardless of the temperature history of newly acquired gas. (b) At the virial radius the temperature and radial flux of inflowing gas is largely unaffected at z=2. However, the spherical covering fraction of inflowing gas at 0.25 rvir decreases substantially, from more than 80% to less than 50%, while the rates of both inflow and outflow increase, indicative of recycling across this boundary. (c) The fractional contribution of smooth accretion to the total accretion rate is lower in the simulation with feedback, by roughly a factor of two across all redshifts. Moreover, the smooth component of gas with a cold temperature history, is entirely suppressed in the feedback run at z<1. (d) The amount of time taken by gas to cross from the virial radius to the galaxy -- the "halo transit time" -- increases in the presence of feedback by a factor of ~2-3, and is notably independent of halo mass. We discuss the possible implications of this invariance for theoretical models of hot halo gas cooling.

preprint2015arXiv

The Impact of Galactic Feedback on the Circumgalactic Medium

Galactic feedback strongly affects the way galactic environments are enriched. We examine this connection by performing a suite of cosmological hydrodynamic simulations, exploring a range of parameters based on the galaxy formation model developed in Vogelsberger et al. 2013 (henceforth V13). We examine the effects of AGN feedback, wind mass loading, wind specific energy, and wind metal-loading on the properties of the circumgalactic medium (CGM) of galaxies with $M_\text{halo} > 10^{11} M_\odot$. Note that while the V13 model was tuned to match observations including the stellar mass function, no explicit tuning was done for the CGM. The wind energy per unit outflow mass has the most significant effect on the CGM enrichment. High energy winds launch metals far beyond the virial radius. AGN feedback also has a significant effect, but only at $z < 3$. We compare to high redshift HI and CIV observations. All our simulations produce the observed number of Damped Lyman-$α$ Absorbers. At lower column density, several of our simulations produce enough Lyman Limit Systems (LLS) $100$ kpc from the galaxy, but in all cases the LLS abundance falls off with distance faster than observations, with too few LLS at $200$ kpc. Further, in all models the CIV abundance drops off too sharply with distance, with too little CIV $100$-$200$ kpc from the galaxy. Higher energy wind models produce more extended CIV but also produce less stars, in tension with star-formation rate density observations. This highlights the fact that circumgalactic observations are a strong constraint on galactic feedback models.

preprint2015arXiv

The merger rate of galaxies in the Illustris Simulation: a comparison with observations and semi-empirical models

We have constructed merger trees for galaxies in the Illustris Simulation by directly tracking the baryonic content of subhalos. These merger trees are used to calculate the galaxy-galaxy merger rate as a function of descendant stellar mass, progenitor stellar mass ratio, and redshift. We demonstrate that the most appropriate definition for the mass ratio of a galaxy-galaxy merger consists in taking both progenitor masses at the time when the secondary progenitor reaches its maximum stellar mass. Additionally, we avoid effects from `orphaned' galaxies by allowing some objects to `skip' a snapshot when finding a descendant, and by only considering mergers which show a well-defined `infall' moment. Adopting these definitions, we obtain well-converged predictions for the galaxy-galaxy merger rate with the following main features, which are qualitatively similar to the halo-halo merger rate except for the last one: a strong correlation with redshift that evolves as $\sim (1+z)^{2.4-2.8}$, a power law with respect to mass ratio, and an increasing dependence on descendant stellar mass, which steepens significantly for descendant stellar masses greater than $\sim 2 \times 10^{11} \, {\rm M_{\odot}}$. These trends are consistent with observational constraints for medium-sized galaxies ($M_{\ast} \gtrsim 10^{10} \, {\rm M_{\odot}}$), but in tension with some recent observations of the close pair fraction for massive galaxies ($M_{\ast} \gtrsim 10^{11} \, {\rm M_{\odot}}$), which report a nearly constant or decreasing evolution with redshift. Finally, we provide a fitting function for the galaxy-galaxy merger rate which is accurate over a wide range of stellar masses, progenitor mass ratios, and redshifts.

preprint2015arXiv

Zooming in on accretion - I. The structure of halo gas

We study the properties of gas in and around 10^12 solar mass halos 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 interaction of 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. We investigate the process of gas virialization and identify different regimes for the heating of gas as it accretes from the intergalactic medium. Haloes at this mass have a well-defined virial shock, associated with a sharp jump in temperature and entropy at ~1.25 r_vir. 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 ~85% of the 4pi sphere. Our findings are relevant for the proper interpretation of observations pertaining to the circumgalactic medium, including evidence for large amounts of cold gas surrounding massive haloes at intermediate redshifts.

preprint2014arXiv

A Constrained Transport Scheme for MHD on Unstructured Static and Moving Meshes

Magnetic fields play an important role in many astrophysical systems and a detailed understanding of their impact on the gas dynamics requires robust numerical simulations. Here we present a new method to evolve the ideal magnetohydrodynamic (MHD) equations on unstructured static and moving meshes that preserves the magnetic field divergence-free constraint to machine precision. The method overcomes the major problems of using a cleaning scheme on the magnetic fields instead, which is non-conservative, not fully Galilean invariant, does not eliminate divergence errors completely, and may produce incorrect jumps across shocks. Our new method is a generalization of the constrained transport (CT) algorithm used to enforce the $\nabla\cdot \mathbf{B}=0$ condition on fixed Cartesian grids. Preserving $\nabla\cdot \mathbf{B}=0$ at the discretized level is necessary to maintain the orthogonality between the Lorentz force and $\mathbf{B}$. The possibility of performing CT on a moving mesh provides several advantages over static mesh methods due to the quasi-Lagrangian nature of the former (i.e., the mesh generating points move with the flow), such as making the simulation automatically adaptive and significantly reducing advection errors. Our method preserves magnetic fields and fluid quantities in pure advection exactly.

preprint2014arXiv

ALMA Imaging of HCN, CS and dust in Arp 220 and NGC 6240

We report ALMA Band 7 (350 GHz) imaging at 0.4 - 0.6arcsec resolution and Band 9 (696 GHz) at ~0.25arcsec resolution of the luminous IR galaxies Arp 220 and NGC 6240. The long wavelength dust continuum is used to estimate ISM masses for Arp 220 East, West and NGC 6240 of 1.9, 4.2 and 1.6x10^9 msun within radii of 69, 65 and 190 pc. The HCN emission was modeled to derive the emissivity distribution as a function of radius and the kinematics of each nuclear disk, yielding dynamical masses consistent with the masses and sizes derived from the dust emission. In Arp 220, the major dust and gas concentrations are at radii less than 50 pc in both counter-rotating nuclear disks. The thickness of the disks in Arp 220estimated from the velocity dispersion and rotation velocities are 10-20 pc and the mean gas densities are n_H2 ~10^5 cm^-3 at R < 50 pc. We develop an analytic treatment for the molecular excitation (including photon trapping), yielding volume densities for both the HCN and CS emission with n_H2 ~2x10^5 cm^-3. The agreement of the mean density from the total mass and size with that required for excitation suggests that the volume is essentially filled with dense gas, i.e. it is not cloudy or like swiss cheese.

preprint2014arXiv

Complexity Phenomena and ROMA of the Magnetospheric Cusp, Hydrodynamic Turbulence, and the Cosmic Web

Dynamic Complexity is a phenomenon exhibited by a nonlinearly interacting system within which multitudes of different sizes of large scale coherent structures emerge, resulting in a globally nonlinear stochastic behavior vastly different from that could be surmised from the underlying equations of interaction. The hallmark of such nonlinear, complex phenomena is the appearance of intermittent fluctuating events with the mixing and distributions of correlated structures at all scales. We briefly review here a relatively recent method, ROMA (rank-ordered multifractal analysis), explicitly constructed to analyze the intricate details of the distribution and scaling of such types of intermittent structures. This method is then applied to the analyses of selected examples related to the dynamical plasmas of the cusp region of the magnetosphere, velocity fluctuations of classical hydrodynamic turbulence, and the distribution of the structures of the cosmic gas obtained through large scale, moving mesh simulations. Differences and similarities of the analyzed results among these complex systems will be contrasted and highlighted. The first two examples have direct relevance to the geospace environment and are summaries of previously reported findings. The third example on the cosmic gas, though involving phenomena much larger in spatiotemporal scales, with its highly compressible turbulent behavior and the unique simulation technique employed in generating the data, provides direct motivations of applying such analysis to studies of similar multifractal processes in various extreme environments. These new results are both exciting and intriguing.

preprint2014arXiv

Damped Lyman-alpha absorbers as a probe of stellar feedback

We examine the abundance, clustering and metallicity of Damped Lyman-alpha Absorbers (DLAs) in a suite of hydrodynamic cosmological simulations using the moving mesh code AREPO. We incorporate models of supernova and AGN feedback, as well as molecular hydrogen formation. We compare our simulations to the column density distribution function at $z=3$, the total DLA abundance at $z=2-4$, the measured DLA bias at $z=2.3$ and the DLA metallicity distribution at $z=2-4$. Our preferred models produce populations of DLAs in good agreement with most of these observations. The exception is the DLA abundance at $z < 3$, which we show requires stronger feedback in $10^{11-12} \, h^{-1} M_\odot$ mass halos. While the DLA population probes a wide range of halo masses, we find the cross-section is dominated by halos of mass $10^{10} - 10^{11} \, h^{-1} M_\odot$ and virial velocities $50 - 100 \;\mathrm{km/s}$. The simulated DLA population has a linear theory bias of $1.7$, whereas the observations require $2.17 \pm 0.2$. We show that non-linear growth increases the bias in our simulations to $2.3$ at $k=1\; \mathrm{Mpc/}h$, the smallest scale observed. The scale-dependence of the bias is, however, very different in the simulations compared against the observations. We show that, of the observations we consider, the DLA abundance and column density function provide the strongest constraints on the feedback model.

preprint2014arXiv

Determining the stellar masses of submillimetre galaxies: the critical importance of star formation histories

Submillimetre (submm) galaxies are among the most rapidly star-forming and most massive high-redshift galaxies; thus, their properties provide important constraints on galaxy evolution models. However, there is still a debate about their stellar masses and their nature in the context of the general galaxy population. To test the reliability of their stellar mass determinations, we used a sample of simulated submm galaxies for which we derived stellar masses via spectral energy distribution (SED) modelling (with Grasil, Magphys, Hyperz and LePhare) adopting various star formation histories (SFHs). We found that the assumption of SFHs with two independent components leads to the most accurate stellar masses. Exponentially declining SFHs (tau) lead to lower masses (albeit still consistent with the true values), while the assumption of single-burst SFHs results in a significant mass underestimation. Thus, we conclude that studies based on the higher masses inferred from fitting the SEDs of real submm galaxies with double SFHs are most likely to be correct, implying that submm galaxies lie on the high-mass end of the main sequence of star-forming galaxies. This conclusion appears robust to assumptions of whether or not submm galaxies are driven by major mergers, since the suite of simulated galaxies modelled here contains examples of both merging and isolated galaxies. We identified discrepancies between the true and inferred stellar ages (rather than the dust attenuation) as the primary determinant of the success/failure of the mass recovery. Regardless of the choice of SFH, the SED-derived stellar masses exhibit a factor of ~2 scatter around the true value; this scatter is an inherent limitation of the SED modelling due to simplified assumptions. Finally, we found that the contribution of active galactic nuclei does not have any significant impact on the derived stellar masses.

preprint2014arXiv

Formation of massive protostars in atomic cooling haloes

We present the highest-resolution three-dimensional simulation to date of the collapse of an atomic cooling halo in the early Universe. We use the moving-mesh code arepo with the primordial chemistry module introduced in Greif (2014), which evolves the chemical and thermal rate equations for over more than 20 orders of magnitude in density. Molecular hydrogen cooling is suppressed by a strong Lyman-Werner background, which facilitates the near-isothermal collapse of the gas at a temperature of about $10^4\,$K. Once the central gas cloud becomes optically thick to continuum emission, it settles into a Keplerian disc around the primary protostar. The initial mass of the protostar is about $0.1\,{\rm M}_\odot$, which is an order of magnitude higher than in minihaloes that cool via molecular hydrogen. The high accretion rate and efficient cooling of the gas catalyse the fragmentation of the disc into a small protostellar system with 5-10 members. After about 12 yr, strong gravitational interactions disrupt the disc and temporarily eject the primary protostar from the centre of the cloud. By the end of the simulation, a secondary clump has collapsed at a distance of $\simeq 150\,$au from the primary clump. If this clump undergoes a similar evolution as the first, the central gas cloud may evolve into a wide binary system. High accretion rates of both the primary and secondary clumps suggest that fragmentation is not a significant barrier for forming at least one massive black hole seed.

preprint2014arXiv

Galaxy mergers on a moving mesh: a comparison with smoothed-particle hydrodynamics

Galaxy mergers have been investigated for decades using smoothed particle hydrodynamics (SPH), but recent work highlighting inaccuracies inherent in the traditional SPH technique calls into question the reliability of previous studies. We explore this issue by comparing a suite of Gadget-3 SPH simulations of idealised (i.e., non-cosmological) isolated discs and galaxy mergers with otherwise identical calculations performed using the moving-mesh code Arepo. When black hole (BH) accretion and active galactic nucleus (AGN) feedback are not included, the star formation histories (SFHs) obtained from the two codes agree well. When BHs are included, the code- and resolution-dependent variations in the SFHs are more significant, but the agreement is still good, and the stellar mass formed over the course of a simulation is robust to variations in the numerical method. During a merger, the gas morphology and phase structure are initially similar prior to the starburst phase. However, once a hot gaseous halo has formed from shock heating and AGN feedback (when included), the agreement is less good. In particular, during the post-starburst phase, the SPH simulations feature more prominent hot gaseous haloes and spurious clumps, whereas with Arepo, gas clumps and filaments are less apparent and the hot halo gas can cool more efficiently. We discuss the origin of these differences and explain why the SPH technique yields trustworthy results for some applications (such as the idealised isolated disc and galaxy merger simulations presented here) but not others (e.g., gas flows onto galaxies in cosmological hydrodynamical simulations).

preprint2014arXiv

Halo Mass and Assembly History Exposed in the Faint Outskirts: the Stellar and Dark Matter Haloes of Illustris Galaxies

We use the Illustris Simulations to gain insight into the build-up of the outer, low-surface brightness regions which surround galaxies. We characterize the stellar haloes by means of the logarithmic slope of the spherically-averaged stellar density profiles, alphaSTARS at z=0, and we relate these slopes to the properties of the underlying Dark-Matter (DM) haloes, their central galaxies, and their assembly histories. We analyze a sample of ~5,000 galaxies resolved with more than 5x10^4 particles each, and spanning a variety of morphologies and halo masses (3x10^11 < Mvir < 10^14 Msun). We find a strong trend between stellar halo slope and total halo mass, where more massive objects have shallower stellar haloes than the less massive ones (-5.5 \pm 0.5 < alphaSTARS <-3.5 \pm 0.2 in the studied mass range). At fixed halo mass, we show that disk-like, blue, young, and more massive galaxies are surrounded by significantly steeper stellar haloes than elliptical, red, older, and less massive galaxies. Overall, the stellar density profiles fall off much more steeply than the underlying DM, and no clear trend holds between stellar slope and DM halo concentration. However, DM haloes which formed more recently, or which accreted larger fractions of stellar mass from infalling satellites, exhibit shallower stellar haloes than their older analogs with similar masses, by up to Delta(alphaSTARS) ~ 0.5-0.7. Our findings, combined with the most recent measurements of the strikingly different stellar power-law indexes for M31 and the Milky Way, appear to favour a massive M31, and a Milky Way characterized by a much quieter accretion history over the past 10 Gyrs than its companion.

preprint2014arXiv

Introducing the Illustris Project: Simulating the coevolution of dark and visible matter in the Universe

We introduce the Illustris Project, a series of large-scale hydrodynamical simulations of galaxy formation. The highest resolution simulation, Illustris-1, covers a volume of $(106.5\,{\rm Mpc})^3$, has a dark mass resolution of ${6.26 \times 10^{6}\,{\rm M}_\odot}$, and an initial baryonic matter mass resolution of ${1.26 \times 10^{6}\,{\rm M}_\odot}$. At $z=0$ gravitational forces are softened on scales of $710\,{\rm pc}$, and the smallest hydrodynamical gas cells have an extent of $48\,{\rm pc}$. We follow the dynamical evolution of $2\times 1820^3$ resolution elements and in addition passively evolve $1820^3$ Monte Carlo tracer particles reaching a total particle count of more than $18$ billion. The galaxy formation model includes: primordial and metal-line cooling with self-shielding corrections, stellar evolution, stellar feedback, gas recycling, chemical enrichment, supermassive black hole growth, and feedback from active galactic nuclei. At $z=0$ our simulation volume contains about $40,000$ well-resolved galaxies covering a diverse range of morphologies and colours including early-type, late-type and irregular galaxies. The simulation reproduces reasonably well the cosmic star formation rate density, the galaxy luminosity function, and baryon conversion efficiency at $z=0$. It also qualitatively captures the impact of galaxy environment on the red fractions of galaxies. The internal velocity structure of selected well-resolved disk galaxies obeys the stellar and baryonic Tully-Fisher relation together with flat circular velocity curves. In the well-resolved regime the simulation reproduces the observed mix of early-type and late-type galaxies. Our model predicts a halo mass dependent impact of baryonic effects on the halo mass function and the masses of haloes caused by feedback from supernova and active galactic nuclei.

preprint2014arXiv

Introducing the Illustris Project: the evolution of galaxy populations across cosmic time

We present an overview of galaxy evolution across cosmic time in the Illustris Simulation. Illustris is an N-body/hydrodynamical simulation that evolves 2*1820^3 resolution elements in a (106.5Mpc)^3 box from cosmological initial conditions down to z=0 using the AREPO moving-mesh code. The simulation uses a state-of-the-art set of physical models for galaxy formation that was tuned to reproduce the z=0 stellar mass function and the history of the cosmic star-formation rate density. We find that Illustris successfully reproduces a plethora of observations of galaxy populations at various redshifts, for which no tuning was performed, and provide predictions for future observations. In particular, we discuss (a) the buildup of galactic mass, showing stellar mass functions and the relations between stellar mass and halo mass from z=7 to z=0, (b) galaxy number density profiles around massive central galaxies out to z=4, (c) the gas and total baryon content of both galaxies and their halos for different redshifts, and as a function of mass and radius, and (d) the evolution of galaxy specific star-formation rates up to z=8. In addition, we (i) present a qualitative analysis of galaxy morphologies from z=5 to z=0, for the stellar as well as the gaseous components, and their appearance in HST mock observations, (ii) follow galaxies selected at z=2 to their z=0 descendants, and quantify their growth and merger histories, and (iii) track massive z=0 galaxies to high redshift and study their joint evolution in star-formation activity and compactness. We conclude with a discussion of several disagreements with observations, and lay out possible directions for future research.

preprint2014arXiv

Numerical Convergence in Smoothed Particle Hydrodynamics

We study the convergence properties of smoothed particle hydrodynamics (SPH) using numerical tests and simple analytic considerations. Our analysis shows that formal numerical convergence is possible in SPH only in the joint limit $N \rightarrow \infty$, $h \rightarrow 0$, and $N_{nb} \rightarrow \infty$, where $N$ is the total number of particles, $h$ is the smoothing length, and $N_{nb}$ is the number of neighbor particles within the smoothing volume used to compute smoothed estimates. Previous work has generally assumed that the conditions $N \rightarrow \infty$ and $h \rightarrow 0$ are sufficient to achieve convergence, while holding $N_{nb}$ fixed. We demonstrate that if $N_{nb}$ is held fixed as the resolution is increased, there will be a residual source of error that does not vanish as $N \rightarrow \infty$ and $h \rightarrow 0$. Formal numerical convergence in SPH is possible only if $N_{nb}$ is increased systematically as the resolution is improved. Using analytic arguments, we derive an optimal compromise scaling for $N_{nb}$ by requiring that this source of error balance that present in the smoothing procedure. For typical choices of the smoothing kernel, we find $N_{nb} \propto N^{1/2}$. This means that if SPH is to be used as a numerically convergent method, the required computational cost does not scale with particle number as $O(N)$, but rather as $O(N^{1+δ})$, where $δ\approx 1/2$, with a weak dependence on the form of the smoothing kernel.

preprint2014arXiv

Properties of galaxies reproduced by a hydrodynamic simulation

Previous simulations of the growth of cosmic structures have broadly reproduced the 'cosmic web' of galaxies that we see in the Universe, but failed to create a mixed population of elliptical and spiral galaxies due to numerical inaccuracies and incomplete physical models. Moreover, because of computational constraints, they were unable to track the small scale evolution of gas and stars to the present epoch within a representative portion of the Universe. Here we report a simulation that starts 12 million years after the Big Bang, and traces 13 billion years of cosmic evolution with 12 billion resolution elements in a volume of $(106.5\,{\rm Mpc})^3$. It yields a reasonable population of ellipticals and spirals, reproduces the distribution of galaxies in clusters and statistics of hydrogen on large scales, and at the same time the metal and hydrogen content of galaxies on small scales.

preprint2014arXiv

Simulated Galaxy Interactions as Probes of Merger Spectral Energy Distributions

We 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 simulated SEDs close to coalescence, while less evolved systems match well with SEDs over a wide range of interaction stages, suggesting that an SED alone is insufficient to identify 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.

preprint2014arXiv

Synthetic Galaxy Images and Spectra from the Illustris Simulation

We present our methods for generating a catalog of 7,000 synthetic images and 40,000 integrated spectra of redshift z = 0 galaxies from the Illustris Simulation. The mock data products are produced by using stellar population synthesis models to assign spectral energy distributions (SED) to each star particle in the galaxies. The resulting synthetic images and integrated SEDs therefore properly reflect the spatial distribution, stellar metallicity distribution, and star formation history of the galaxies. From the synthetic data products it is possible to produce monochromatic or color-composite images, perform SED fitting, classify morphology, determine galaxy structural properties, and evaluate the impacts of galaxy viewing angle. The main contribution of this paper is to describe the production, format, and composition of the image catalog that makes up the Illustris Simulation Obsevatory. As a demonstration of this resource, we derive galactic stellar mass estimates by applying the SED fitting code FAST to the synthetic galaxy products, and compare the derived stellar masses against the true stellar masses from the simulation. We find from this idealized experiment that systematic biases exist in the photometrically derived stellar mass values that can be reduced by using a fixed metallicity in conjunction with a minimum galaxy age restriction.

preprint2014arXiv

The colors of satellite galaxies in the Illustris Simulation

Observationally, the fraction of blue satellite galaxies decreases steeply with host halo mass, and their radial distribution around central galaxies is significantly shallower in massive (M_* >10e11M_sun) than in Milky Way like systems. Theoretical models, based primarily on semi-analytical techniques, have had a long-standing problem with reproducing these trends, instead predicting too few blue satellites in general but also estimating a radial distribution that is too shallow, regardless of primary mass. In this Letter, we use the Illustris cosmological simulation to study the properties of satellite galaxies around isolated primaries. For the first time, we find good agreement between theory and observations. We identify the main source of this success relative to earlier work to be a consequence of the large gas contents of satellites at infall, a factor ~5-10 times larger than in semi-analytical models. Because of their relatively large gas reservoirs, satellites can continue to form stars long after infall, with a typical timescale for star-formation to be quenched ~2 Gyr in groups but more than ~5 Gyr for satellites around Milky Way like primaries. The gas contents we infer are consistent with z=0 observations of HI gas in galaxies, although we find large discrepancies among reported values in the literature. A testable prediction of our model is that the gas-to-stellar mass ratio of satellite progenitors should vary only weakly with cosmic time.

preprint2014arXiv

The star formation main sequence and stellar mass assembly of galaxies in the Illustris simulation

Understanding the physical processes that drive star formation is a key challenge for galaxy formation models. In this article we study the tight correlation between the star formation rate (SFR) and stellar mass of galaxies at a given redshift, how halo growth influences star formation, and star formation histories of individual galaxies. We study these topics using Illustris, a state-of-the-art cosmological hydrodynamical simulation of galaxy formation. Illustris reproduces the observed relation (the star formation main sequence; SFMS) between SFR and stellar mass at redshifts z=0 and z=4, but at intermediate redshifts of z~2, the simulated SFMS has a significantly lower normalisation than reported by observations. The scatter in the relation is consistent with the observed scatter. However, the fraction of outliers above the SFR-stellar mass relation in Illustris is less than that observed. Galaxies with halo masses of ~10^{12} solar masses dominate the SFR density of the Universe, in agreement with the results of abundance matching. Furthermore, more-massive galaxies tend to form the bulk of their stars at high redshift, which indicates that `downsizing' occurs in Illustris. We also studied the star formation histories of individual galaxies, including the use of a principal component analysis decomposition. We find that for fixed stellar mass, galaxies that form earlier have more-massive black holes at z=0, indicating that star formation and black hole growth are tightly linked processes in Illustris. While many of the properties of normal star-forming galaxies are well-reproduced in the Illustris simulation, forming a realistic population of starbursts will likely require higher resolution and probably a more sophisticated treatment of star formation and feedback from stars and black holes.

preprint2014arXiv

The total infrared luminosity may significantly overestimate the star formation rate of recently quenched galaxies

The total infrared (IR) luminosity is very useful for estimating the star formation rate (SFR) of galaxies, but converting the IR luminosity into an SFR relies on assumptions that do not hold for all galaxies. We test the effectiveness of the IR luminosity as an SFR indicator by applying it to synthetic spectral energy distributions generated from three-dimensional hydrodynamical simulations of isolated disc galaxies and galaxy mergers. In general, the SFR inferred from the IR luminosity agrees well with the true instantaneous SFR of the simulated galaxies. However, for the major mergers in which a strong starburst is induced, the SFR inferred from the IR luminosity can overestimate the instantaneous SFR during the post-starburst phase by greater than two orders of magnitude. Even though the instantaneous SFR decreases rapidly after the starburst, the stars that were formed in the starburst remain dust-obscured and thus produce significant IR luminosity. Consequently, use of the IR luminosity as an SFR indicator may cause one to conclude that post-starburst galaxies are still star-forming, whereas in reality, star formation was recently quenched.

preprint2013arXiv

A discontinuous Galerkin method for solving the fluid and MHD equations in astrophysical simulations

A discontinuous Galerkin (DG) method suitable for large-scale astrophysical simulations on Cartesian meshes as well as arbitrary static and moving Voronoi meshes is presented. Most major astrophysical fluid dynamics codes use a finite volume (FV) approach. We demonstrate that the DG technique offers distinct advantages over FV formulations on both static and moving meshes. The DG method is also easily generalized to higher than second-order accuracy without requiring the use of extended stencils to estimate derivatives (thereby making the scheme highly parallelizable). We implement the technique in the AREPO code for solving the fluid and the magnetohydrodynamic (MHD) equations. By examining various test problems, we show that our new formulation provides improved accuracy over FV approaches of the same order, and reduces post-shock oscillations and artificial diffusion of angular momentum. In addition, the DG method makes it possible to represent magnetic fields in a locally divergence-free way, improving the stability of MHD simulations and moderating global divergence errors, and is a viable alternative for solving the MHD equations on meshes where Constrained-Transport (CT) cannot be applied. We find that the DG procedure on a moving mesh is more sensitive to the choice of slope limiter than is its FV method counterpart. Therefore, future work to improve the performance of the DG scheme even further will likely involve the design of optimal slope limiters. As presently constructed, our technique offers the potential of improved accuracy in astrophysical simulations using the moving mesh AREPO code as well as those employing adaptive mesh refinement (AMR).

preprint2013arXiv

A model for cosmological simulations of galaxy formation physics

We present a new comprehensive model of the physics of galaxy formation designed for large-scale hydrodynamical simulations of structure formation using the moving mesh code AREPO. Our model includes primordial and metal line cooling with self-shielding corrections, stellar evolution and feedback processes, gas recycling, chemical enrichment, a novel subgrid model for the metal loading of outflows, black hole (BH) seeding, BH growth and merging procedures, quasar- and radio-mode feedback, and a prescription for radiative electro-magnetic (EM) feedback from active galactic nuclei (AGN). The metal mass loading of outflows can be adjusted independently of the wind mass loading. This is required to simultaneously reproduce the stellar mass content of low mass haloes and their gas oxygen abundances. Radiative EM AGN feedback is implemented assuming an average spectral energy distribution and a luminosity-dependent scaling of obscuration effects. This form of feedback suppresses star formation more efficiently than continuous thermal quasar-mode feedback alone, but is less efficient than mechanical radio-mode feedback in regulating star formation in massive haloes. We contrast simulation predictions for different variants of our galaxy formation model with key observations. Our best match model reproduces, among other things, the cosmic star formation history, the stellar mass function, the stellar mass - halo mass relation, g-, r-, i-, z-band SDSS galaxy luminosity functions, and the Tully-Fisher relation. We can achieve this success only if we invoke very strong forms of stellar and AGN feedback such that star formation is adequately reduced in both low and high mass systems. In particular, the strength of radio-mode feedback needs to be increased significantly compared to previous studies to suppress efficient cooling in massive, metal-enriched haloes.

preprint2013arXiv

A physical model for cosmological simulations of galaxy formation: multi-epoch validation

We present a multi-epoch analysis of the galaxy populations formed within the cosmological hydrodynamical simulations presented in Vogelsberger et al. (2013). These simulations explore the performance of a recently implemented feedback model which includes primordial and metal line radiative cooling with self-shielding corrections; stellar evolution with associated mass loss and chemical enrichment; feedback by stellar winds; black hole seeding, growth and merging; and AGN quasar- and radio-mode heating with a phenomenological prescription for AGN electro-magnetic feedback. We illustrate the impact of the model parameter choices on the resulting simulated galaxy population properties at high and intermediate redshifts. We demonstrate that our scheme is capable of producing galaxy populations that broadly reproduce the observed galaxy stellar mass function extending from redshift z=0 to z=3. We also characterise the evolving galactic B-band luminosity function, stellar mass to halo mass ratio, star formation main sequence, Tully-Fisher relation, and gas-phase mass-metallicity relation and confront them against recent observational estimates. This detailed comparison allows us to validate elements of our feedback model, while also identifying areas of tension that will be addressed in future work.

preprint2013arXiv

Empirical Constraints for the Magnitude and Composition of Galactic Winds

Galactic winds are a key physical mechanism for understanding galaxy formation and evolution, yet empirical and theoretical constraints for the character of winds are limited and discrepant. Recent empirical models find that local star-forming galaxies have a deficit of oxygen that scales with galaxy stellar mass. The oxygen deficit provides unique empirical constraints on the magnitude of mass loss, composition of outflowing material and metal reaccretion onto galaxies. We formulate the oxygen deficit constraints so they may be easily implemented into theoretical models of galaxy evolution. We parameterize an effective metal loading factor which combines the uncertainties of metal outflows and metal reaccretion into a single function of galaxy virial velocity. We determine the effective metal loading factor by forward-fitting the oxygen deficit. The effective metal loading factor we derive has important implications for the implementation of mass loss in models of galaxy evolution.

preprint2013arXiv

Following the flow: tracer particles in astrophysical fluid simulations

We present two numerical schemes for passive tracer particles in the hydrodynamical moving-mesh code AREPO, and compare their performance for various problems, from simple setups to cosmological simulations. The purpose of tracer particles is to allow the flow to be followed in a Lagrangian way, tracing the evolution of the fluid with time, and allowing the thermodynamical history of individual fluid parcels to be recorded. We find that the commonly-used `velocity field tracers', which are advected using the fluid velocity field, do not in general follow the mass flow correctly, and explain why this is the case. This method can result in orders-of-magnitude biases in simulations of driven turbulence and in cosmological simulations, rendering the velocity field tracers inappropriate for following these flows. We then discuss a novel implementation of `Monte Carlo tracers', which are moved along with fluid cells, and are exchanged probabilistically between them following the mass flux. This method reproduces the mass distribution of the fluid correctly. The main limitation of this approach is that it is more diffusive than the fluid itself. Nonetheless, we show that this novel approach is more reliable than what has been employed previously and demonstrate that it is appropriate for following hydrodynamical flows in mesh-based codes. The Monte Carlo tracers can also naturally be transferred between fluid cells and other types of particles, such as stellar particles, so that the mass flow in cosmological simulations can be followed in its entirety.

preprint2013arXiv

Global Star Formation Rates and Dust Emission Over the Galaxy Interaction Sequence

We measured and modeled the spectral energy distributions (SEDs) in 28 bands from the ultraviolet to the far-infrared (FIR) for 31 interacting galaxies in 14 systems. The sample is drawn from the Spitzer Interacting Galaxy Survey, which probes a range of galaxy interaction parameters at multiple wavelengths with an emphasis on the infrared bands. The subset presented in this paper consists of all galaxies for which FIR Herschel SPIRE observations are publicly available. Our SEDs combine the Herschel photometry with multi-wavelength data from Spitzer, GALEX, Swift UVOT, and 2MASS. While the shapes of the SEDs are broadly similar across our sample, strongly interacting galaxies typically have more mid-infrared emission relative to their near-infrared and FIR emission than weakly or moderately interacting galaxies. We modeled the full SEDs to derive host galaxy star formation rates (SFR), specific star formation rates (sSFR), stellar masses, dust temperatures, dust luminosities, and dust masses. We find increases in the dust luminosity and mass, SFR, and cold (15-25 K) dust temperature as the interaction progresses from moderately to strongly interacting and between non- interacting and strongly interacting galaxies. We also find increases in the SFR between weakly and strongly interacting galaxies. In contrast, the sSFR remains unchanged across all the interaction stages. The SFR derived from the SED modeling agrees well with rates estimated by proportionality relations that depend on infrared emission. (abridged)

preprint2013arXiv

Modeling Mid-Infrared Diagnostics of Obscured Quasars and Starbursts

We analyze the link between active galactic nuclei (AGN) and mid-infrared flux using dust radiative transfer calculations of starbursts realized in hydrodynamical simulations. Focusing on the effects of galaxy dust, we evaluate diagnostics commonly used to disentangle AGN and star formation in ultraluminous infrared galaxies (ULIRGs). We examine these quantities as a function of time, viewing angle, dust model, AGN spectrum, and AGN strength in merger simulations representing two possible extremes of the ULIRG population: one is a typical gas-rich merger at z ~ 0, and the other is characteristic of extremely obscured starbursts at z ~ 2 to 4. This highly obscured burst begins star-formation-dominated with significant PAH emission, and ends with a ~10^9 yr period of red near-IR colors. At coalescence, when the AGN is most luminous, dust obscures the near-infrared AGN signature, reduces the relative emission from polycyclic aromatic hydrocarbons (PAHs), and enhances the 9.7 micron absorption by silicate grains. Although generally consistent with previous interpretations, our results imply none of these indicators can unambiguously estimate the AGN luminosity fraction in all cases. Motivated by the simulations, we show that a combination of the extinction feature at 9.7 micron, the PAH strength, and a near-infrared slope can simultaneously constrain the AGN fraction and dust grain distribution for a wide range of obscuration. We find that this indicator, accessible to the James Webb Space Telescope, may estimate the AGN power as tightly as the hard X-ray flux alone, thereby providing a valuable future cross-check and constraint for large samples of distant ULIRGs.

preprint2013arXiv

Moving mesh cosmology: tracing cosmological gas accretion

We investigate the nature of gas accretion onto haloes and galaxies at z=2 using cosmological hydrodynamic simulations run with the moving mesh code AREPO. Implementing a Monte Carlo tracer particle scheme to determine the origin and thermodynamic history of accreting gas, we make quantitative comparisons to an otherwise identical simulation run with the smoothed particle hydrodynamics (SPH) code GADGET-3. Contrasting these two numerical approaches, we find significant physical differences in the thermodynamic history of accreted gas in haloes above 10^10.5 solar masses. In agreement with previous work, GADGET simulations show a cold fraction near unity for galaxies forming in massive haloes, implying that only a small percentage of accreted gas heats to an appreciable fraction of the virial temperature during accretion. The same galaxies in AREPO show a much lower cold fraction, <20% in haloes above 10^11 solar masses. This results from a hot gas accretion rate which, at this same halo mass, is an order of magnitude larger than with GADGET, while the cold accretion rate is also lower. These discrepancies increase for more massive systems, and we explain both as due to numerical inaccuracies in the standard formulation of SPH. We also observe that the relatively sharp transition from cold to hot mode dominated accretion, at a halo mass of ~10^11, is a consequence of comparing past gas temperatures to a constant threshold value independent of virial temperature. Examining the spatial distribution of accreting gas, we find that gas filaments in GADGET tend to remain collimated and flow coherently to small radii, or artificially fragment and form a large number of purely numerical "blobs". Similar gas streams in AREPO show increased heating and disruption at 0.25-0.5 virial radii and contribute to the hot gas accretion rate in a manner distinct from classical cooling flows.

preprint2013arXiv

Resolving The Generation of Starburst Winds in Galaxy Mergers

We study galaxy super-winds driven in major mergers, using pc-resolution simulations with detailed models for stellar feedback that can self-consistently follow the formation/destruction of GMCs and generation of winds. The models include molecular cooling, star formation at high densities in GMCs, and gas recycling and feedback from SNe (I&II), stellar winds, and radiation pressure. We study mergers of systems from SMC-like dwarfs and Milky Way analogues to z~2 starburst disks. Multi-phase super-winds are generated in all passages, with outflow rates up to ~1000 M_sun/yr. However, the wind mass-loading efficiency (outflow rate divided by SFR) is similar to that in isolated galaxy counterparts of each merger: it depends more on global galaxy properties (mass, size, escape velocity) than on the dynamical state of the merger. Winds tend to be bi- or uni-polar, but multiple 'events' build up complex morphologies with overlapping, differently-oriented bubbles/shells at a range of radii. The winds have complex velocity and phase structure, with material at a range of speeds up to ~1000 km/s, and a mix of molecular, ionized, and hot gas that depends on galaxy properties and different feedback mechanisms. These simulations resolve a problem in some 'sub-grid' models, where simple wind prescriptions can dramatically suppress merger-induced starbursts. But despite large mass-loading factors (>~10) in the winds, the peak SFRs are comparable to those in 'no wind' simulations. Wind acceleration does not act equally, so cold dense gas can still lose angular momentum and form stars, while blowing out gas that would not have participated in the starburst in the first place. Considerable wind material is not unbound, and falls back on the disk at later times post-merger, leading to higher post-starburst SFRs in the presence of stellar feedback. This may require AGN feedback to explain galaxy quenching.

preprint2013arXiv

Self-Perpetuating Spiral Arms in Disk Galaxies

The causes of spiral structure in galaxies remain uncertain. Leaving aside the grand bisymmetric spirals with their own well-known complications, here we consider the possibility that multi-armed spiral features originate from density inhomogeneities orbiting within disks. Using high-resolution N-body simulations, we follow the motions of stars under the influence of gravity, and show that mass concentrations with properties similar to those of giant molecular clouds can induce the development of spiral arms through a process termed swing amplification. However, unlike in earlier work, we demonstrate that the eventual response of the disk can be highly non-linear, significantly modifying the formation and longevity of the resulting patterns. Contrary to expectations, ragged spiral structures can thus survive at least in a statistical sense long after the original perturbing influence has been removed.

preprint2013arXiv

Star Formation in Galaxy Mergers with Realistic Models of Stellar Feedback & the Interstellar Medium

We use simulations with realistic models for stellar feedback to study galaxy mergers. These high resolution (1 pc) simulations follow formation and destruction of individual GMCs and star clusters. The final starburst is dominated by in situ star formation, fueled by gas which flows inwards due to global torques. The resulting high gas density results in rapid star formation. The gas is self gravitating, and forms massive (~10^10 M_sun) GMCs and subsequent super-starclusters (masses up to 10^8 M_sun). However, in contrast to some recent simulations, the bulk of new stars which eventually form the central bulge are not born in superclusters which then sink to the center of the galaxy, because feedback efficiently disperses GMCs after they turn several percent of their mass into stars. Most of the mass that reaches the nucleus does so in the form of gas. The Kennicutt-Schmidt law emerges naturally as a consequence of feedback balancing gravitational collapse, independent of the small-scale star formation microphysics. The same mechanisms that drive this relation in isolated galaxies, in particular radiation pressure from IR photons, extend over seven decades in SFR to regulate star formation in the most extreme starbursts (densities >10^4 M_sun/pc^2). Feedback also drives super-winds with large mass loss rates; but a significant fraction of the wind material falls back onto the disks at later times, leading to higher post-starburst SFRs in the presence of stellar feedback. Strong AGN feedback is required to explain sharp cutoffs in star formation rate. We compare the predicted relic structure, mass profile, morphology, and efficiency of disk survival to simulations which do not explicitly resolve GMCs or feedback. Global galaxy properties are similar, but sub-galactic properties and star formation rates can differ significantly.

preprint2012arXiv

A General Model for the CO-H2 Conversion Factor in Galaxies with Applications to the Star Formation Law

The most common means of converting an observed CO line intensity into a molecular gas mass requires the use of a conversion factor (Xco). While in the Milky Way this quantity does not appear to vary significantly, there is good reason to believe that Xco will depend on the larger-scale galactic environment. Utilising numerical models, we investigate how varying metallicities, gas temperatures and velocity dispersions in galaxies impact the way CO line emission traces the underlying H2 gas mass, and under what circumstances Xco may differ from the Galactic mean value. We find that, due to the combined effects of increased gas temperature and velocity dispersion, Xco is depressed below the Galactic mean in high surface density environments such as ULIRGs. In contrast, in low metallicity environments, Xco tends to be higher than in the Milky Way, due to photodissociation of CO in metal-poor clouds. At higher redshifts, gas-rich discs may have gravitationally unstable clumps which are warm (due to increased star formation) and have elevated velocity dispersions. These discs tend to have Xco values ranging between present-epoch gas-rich mergers and quiescent discs at low-z. This model shows that on average, mergers do have lower Xco values than disc galaxies, though there is significant overlap. Xco varies smoothly with the local conditions within a galaxy, and is not a function of global galaxy morphology. We combine our results to provide a general fitting formula for Xco as a function of CO line intensity and metallicity. We show that replacing the traditional approach of using one constant Xco for starbursts and another for discs with our best-fit function produces star formation laws that are continuous rather than bimodal, and that have significantly reduced scatter.

preprint2012arXiv

A Simple Model for the Density Profiles of Isolated Dark Matter Halos

We explore the possibility that the density profiles of elliptical galaxies and cold dark matter (CDM) halos found in cosmological simulations can be understood in terms of the same physical process, collisionless gravitational collapse. To investigate this, we study a simplified model, the collapse of a perfectly cold Plummer sphere. First, we examine an N-body simulation of this model with particles constrained to purely radial orbits. This results in a final state characterized by a profile slightly steeper than ρ\propto r^{-2} at small radii and behaving as ρ\propto r^{-4} at large radii, which can be understood in terms of simple analytic arguments. Next, we repeat our simulation without the restriction of radial orbits. This results in a shallower inner density profile, like those found in elliptical galaxies and CDM halos. We attribute this change to the radial orbit instability (ROI) and propose a form of the distribution function (DF) motivated by a physical picture of collapse. As evidence of the link between our model and CDM halos, we find that our collapse simulation has a final state with pseudo-phase-space density which scales roughly as ρ/σ^3 \propto r^{-1.875}, like that observed in CDM halos from cosmological simulations (Navarro et al. 2010). The velocity anisotropy profile is also qualitatively similar to that found near the centers of these halos. We argue that the discrepancy at large radii (where CDM halos scale as ρ\propto r^{-3}) is due to the presence of the cosmological background or continued infall. This leads us to predict that the outer CDM halo density profile is not "universal," but instead depends on cosmological environment (be it an underdense void or overdense region).

preprint2012arXiv

How to distinguish starbursts and quiescently star-forming galaxies: The `bimodal' submillimetre galaxy population as a case study

In recent work (arXiv:1101.0002) we have suggested that the high-redshift (z ~ 2-4) bright submillimetre galaxy (SMG) population is heterogeneous, with major mergers contributing both at early stages, where quiescently star-forming discs are blended into one submm source (`galaxy-pair SMGs'), and late stages, where mutual tidal torques drive gas inflows and cause strong starbursts. Here we combine hydrodynamic simulations of major mergers with 3-D dust radiative transfer calculations to determine observational diagnostics that can distinguish between quiescently star-forming SMGs and starburst SMGs via integrated data alone. We fit the far-IR SEDs of the simulated galaxies with the optically thin single-temperature modified blackbody, the full form of the single-temperature modified blackbody, and a power-law temperature-distribution model. The effective dust temperature, T_dust, and power-law index of the dust emissivity in the far-IR, β, derived can significantly depend on the fitting form used, and the intrinsic β of the dust is not recovered. However, for all forms used here, there is a T_dust above which almost all simulated galaxies are starbursts, so a T_dust cut is very effective at selecting starbursts. Simulated merger-induced starbursts also have higher L_IR/M_gas and L_IR/L_FUV than quiescently star-forming galaxies and lie above the star formation rate-stellar mass relation. These diagnostics can be used to test our claim that the SMG population is heterogeneous and to observationally determine what star formation mode dominates a given galaxy population. We comment on applicability of these diagnostics to ULIRGs that would not be selected as SMGs. These `hot-dust ULIRGs' are typically starburst galaxies lower in mass than SMGs, but they can also simply be SMGs observed from a different viewing angle.

preprint2012arXiv

Moving mesh cosmology: characteristics of galaxies and haloes

We discuss cosmological hydrodynamic simulations of galaxy formation performed with the new moving-mesh code AREPO, which promises higher accuracy compared with the traditional SPH technique that has been widely employed for this problem. We use an identical set of physics in corresponding simulations carried out with the well-tested SPH code GADGET, adopting also the same high-resolution gravity solver. We are thus able to compare both simulation sets on an object-by-object basis, allowing us to cleanly isolate the impact of different hydrodynamical methods on galaxy and halo properties. In accompanying papers, we focus on an analysis of the global baryonic statistics predicted by the simulation codes, (Vogelsberger et al. 2011) and complementary idealized simulations that highlight the differences between the hydrodynamical schemes (Sijacki et al. 2011). Here we investigate their influence on the baryonic properties of simulated galaxies and their surrounding haloes. We find that AREPO leads to significantly higher star formation rates for galaxies in massive haloes and to more extended gaseous disks in galaxies, which also feature a thinner and smoother morphology than their GADGET counterparts. Consequently, galaxies formed in AREPO have larger sizes and higher specific angular momentum than their SPH correspondents. The more efficient cooling flows in AREPO yield higher densities and lower entropies in halo centers (and the opposite trend in halo outskirts) leading to higher star formation rates of massive galaxies. While both codes agree to acceptable accuracy on a number of baryonic properties of cosmic structures, our results clearly demonstrate that galaxy formation simulations greatly benefit from the use of more accurate hydrodynamical techniques such as AREPO.

preprint2012arXiv

Moving mesh cosmology: numerical techniques and global statistics

We present the first hydrodynamical simulations of structure formation using the new moving mesh code AREPO and compare the results with GADGET simulations based on a traditional smoothed particle hydrodynamics (SPH) technique. The two codes share the same Tree-PM gravity solver and include identical sub-resolution physics, but employ different methods to solve the equations of hydrodynamics. This allows us to assess the impact of hydro-solver uncertainties on the results of cosmological studies of galaxy formation. We focus on predictions for global baryon statistics, such as the cosmic star formation rate density, after we introduce our simulation suite and numerical methods. Properties of individual galaxies and haloes are examined by Keres et al. (2011), while a third paper by Sijacki et al. (2011) uses idealised simulations to analyse the differences between the hydrodynamical schemes. We find that the global baryon statistics differ significantly between the two simulation approaches. AREPO shows higher star formation rates at late times, lower mean temperatures, and different gas mass fractions in characteristic phases of the intergalactic medium, in particular a reduced amount of hot gas. Although both codes use the same implementation of cooling, more gas cools out of haloes in AREPO compared with GADGET towards low redshifts. We show that this is caused by a higher heating rate with SPH in the outer parts of haloes, owing to viscous dissipation of SPH's inherent sonic velocity noise and SPH's efficient damping of subsonic turbulence injected in the halo infall region, and because of a higher efficiency of gas stripping in AREPO. As a result of such differences, AREPO leads also to more disk-like morphologies compared to GADGET. Our results indicate that inaccuracies in hydrodynamic solvers can lead to comparatively large systematic differences.

preprint2012arXiv

Moving Mesh Cosmology: Properties of Gas Disks

We compare the structural properties of galaxies formed in cosmological simulations using the smoothed particle hydrodynamics (SPH) code GADGET with those using the moving-mesh code AREPO. Both codes employ identical gravity solvers and the same sub-resolution physics but use very different methods to track the hydrodynamic evolution of gas. This permits us to isolate the effects of the hydro solver on the formation and evolution of galactic gas disks in GADGET and AREPO haloes with comparable numerical resolution. In a matching sample of GADGET and AREPO haloes we fit simulated gas disks with exponential profiles. We find that the cold gas disks formed using the moving mesh approach have systematically larger disk scale lengths and higher specific angular momenta than their GADGET counterparts across a wide range in halo masses. For low mass galaxies differences between the properties of the simulated galaxy disks are caused by an insufficient number of resolution elements which lead to the artificial angular momentum transfer in our SPH calculation. We however find that galactic disks formed in massive halos, resolved with 10^6 particles/cells, are still systematically smaller in the GADGET run by a factor of ~2. The reasons for this are: 1) The excessive heating of haloes close to the cooling radius due to spurious dissipation of the subsonic turbulence in GADGET; and 2) The efficient delivery of low angular momentum gaseous blobs to the bottom of the potential well. While this large population of gaseous blobs in GADGET originates from the filaments which are pressure confined and fragment due to the SPH surface tension while infalling into hot halo atmospheres, it is essentially absent in the moving mesh calculation, clearly indicating numerical rather than physical origin of the blob material.

preprint2012arXiv

Moving mesh cosmology: properties of neutral hydrogen in absorption

We examine the distribution of neutral hydrogen in cosmological simulations carried out with the new moving-mesh code AREPO and compare it with the corresponding GADGET simulations based on the smoothed particle hydrodynamics (SPH) technique. The two codes use identical gravity solvers and baryonic physics implementations, but very different methods for solving the Euler equations, allowing us to assess how numerical effects associated with the hydro-solver impact the results of simulations. Here we focus on an analysis of the neutral gas, as detected in quasar absorption lines. We find that the high column density regime probed by Damped Lyman-alpha (DLA) and Lyman Limit Systems (LLS) exhibits significant differences between the codes. GADGET produces spurious artefacts in large halos in the form of gaseous clumps, boosting the LLS cross-section. Furthermore, it forms halos with denser central baryonic cores than AREPO, which leads to a substantially greater DLA cross-section from smaller halos. AREPO thus produces a significantly lower cumulative abundance of DLAs, which is intriguingly in much closer agreement with observations. The column density function, however, is not altered enough to significantly reduce the discrepancy with the observed value. For the low column density gas probed by the Lyman-alpha forest, the codes differ only at the level of a few percent, suggesting that this regime is quite well described by both methods, a fact that is reassuring for the many Lyman-alpha studies carried out with SPH thus far. While the residual differences are smaller than the errors on current Lyman-alpha forest data, we note that this will likely change for future precision experiments.

preprint2012arXiv

Moving mesh cosmology: the hydrodynamics of galaxy formation

We present a detailed comparison between the well-known SPH code GADGET and the new moving-mesh code AREPO on a number of hydrodynamical test problems. Through a variety of numerical experiments we establish a clear link between test problems and systematic numerical effects seen in cosmological simulations of galaxy formation. Our tests demonstrate deficiencies of the SPH method in several sectors. These accuracy problems not only manifest themselves in idealized hydrodynamical tests, but also propagate to more realistic simulation setups of galaxy formation, ultimately affecting gas properties in the full cosmological framework, as highlighted in papers by Vogelsberger et al. (2011) and Keres et al. (2011). We find that an inadequate treatment of fluid instabilities in GADGET suppresses entropy generation by mixing, underestimates vorticity generation in curved shocks and prevents efficient gas stripping from infalling substructures. In idealized tests of inside-out disk formation, the convergence rate of gas disk sizes is much slower in GADGET due to spurious angular momentum transport. In simulations where we follow the interaction between a forming central disk and orbiting substructures in a halo, the final disk morphology is strikingly different. In AREPO, gas from infalling substructures is readily depleted and incorporated into the host halo atmosphere, facilitating the formation of an extended central disk. Conversely, gaseous sub-clumps are more coherent in GADGET simulations, morphologically transforming the disk as they impact it. The numerical artefacts of the SPH solver are particularly severe for poorly resolved flows, and thus inevitably affect cosmological simulations due to their hierarchical nature. Our numerical experiments clearly demonstrate that AREPO delivers a physically more reliable solution.

preprint2012arXiv

Multi-Dimensional, Compressible Viscous Flow on a Moving Voronoi Mesh

Numerous formulations of finite volume schemes for the Euler and Navier-Stokes equations exist, but in the majority of cases they have been developed for structured and stationary meshes. In many applications, more flexible mesh geometries that can dynamically adjust to the problem at hand and move with the flow in a (quasi) Lagrangian fashion would, however, be highly desirable, as this can allow a significant reduction of advection errors and an accurate realization of curved and moving boundary conditions. Here we describe a novel formulation of viscous continuum hydrodynamics that solves the equations of motion on a Voronoi mesh created by a set of mesh-generating points. The points can move in an arbitrary manner, but the most natural motion is that given by the fluid velocity itself, such that the mesh dynamically adjusts to the flow. Owing to the mathematical properties of the Voronoi tessellation, pathological mesh-twisting effects are avoided. Our implementation considers the full Navier-Stokes equations and has been realized in the AREPO code both in 2D and 3D. We propose a new approach to compute accurate viscous fluxes for a dynamic Voronoi mesh, and use this to formulate a finite volume solver of the Navier-Stokes equations. Through a number of test problems, including circular Couette flow and flow past a cylindrical obstacle, we show that our new scheme combines good accuracy with geometric flexibility, and hence promises to be competitive with other highly refined Eulerian methods. This will in particular allow astrophysical applications of the AREPO code where physical viscosity is important, such as in the hot plasma in galaxy clusters, or for viscous accretion disk models.

preprint2012arXiv

Stellar Feedback & Bulge Formation in Clumpy Disks

We use numerical simulations of isolated galaxies to study the effects of stellar feedback on the formation and evolution of giant star-forming gas 'clumps' in high-redshift, gas-rich galaxies. Such galactic disks are unstable to the formation of bound gas-rich clumps whose properties initially depend only on global disk properties, not the microphysics of feedback. In simulations without stellar feedback, clumps turn an order-unity fraction of their mass into stars and sink to the center, forming a large bulge and kicking most of the stars out into a much more extended stellar envelope. By contrast, strong radiative stellar feedback disrupts even the most massive clumps after they turn ~10-20% of their mass into stars, in a timescale of ~10-100 Myr, ejecting some material into a super-wind and recycling the rest of the gas into the diffuse ISM. This suppresses the bulge formation rate by direct 'clump coalescence' by a factor of several. However, the galactic disks do undergo significant internal evolution in the absence of mergers: clumps form and disrupt continuously and torque gas to the galactic center. The resulting evolution is qualitatively similar to bar/spiral evolution in simulations with a more homogeneous ISM.

preprint2012arXiv

Submillimetre galaxies in a hierarchical universe: number counts, redshift distribution, and implications for the IMF

High-redshift submillimetre galaxies (SMGs) are some of the most rapidly star-forming galaxies in the Universe. Historically, galaxy formation models have had difficulty explaining the observed number counts of SMGs. We combine a semi-empirical model with 3-D hydrodynamical simulations and 3-D dust radiative transfer to predict the number counts of unlensed SMGs. Because the stellar mass functions, gas and dust masses, and sizes of our galaxies are constrained to match observations, we can isolate uncertainties related to the dynamical evolution of galaxy mergers and the dust radiative transfer. The number counts and redshift distributions predicted by our model agree well with observations. Isolated disc galaxies dominate the faint (S_1.1 < ~1 mJy, or S_850 < ~2 mJy) population. The brighter sources are a mix of merger-induced starbursts and galaxy-pair SMGs; the latter subpopulation accounts for ~30 - 50 per cent of all SMGs at all S_1.1 > ~0.5 mJy (S_850 > ~1 mJy). The mean redshifts are ~3.0 - 3.5, depending on the flux cut, and the brightest sources tend to be at higher redshifts. Because the galaxy-pair SMGs will be resolved into multiple fainter sources by ALMA, the bright ALMA counts should be as much as 2 times less than those observed using single-dish telescopes. The agreement between our model, which uses a Kroupa IMF, and observations suggests that the IMF in high-redshifts starbursts need not be top-heavy; if the IMF were top-heavy, our model would over-predict the number counts. We conclude that the difficulty some models have reproducing the observed SMG counts is likely indicative of more general problems -- such as an under-prediction of the abundance of massive galaxies or a star formation rate--stellar mass relation normalisation lower than that observed -- rather than a problem specific to the SMG population.

preprint2012arXiv

The Effects of Patchy Reionization on Satellite Galaxies of the Milky Way

We combine the high-resolution Aquarius simulations with three-dimensional models of reionization based on the initial density field of the Aquarius parent simulation, Millennium-II, to study the impact of patchy reionization on the faint satellite population of Milky Way halos. Because the Aquarius suite consists of zoom-in simulations of halos in the Millennium-II volume, we follow the formation of substructure and the growth of reionization bubbles due to the larger environment simultaneously, and thereby determine the reionization redshifts of satellite candidates. We do this for four different reionization models, and also compare results to instantaneous reionization. Using a simple procedure for selecting satellites and assigning luminosities in the simulations, we compare the resulting satellite populations. We find that the overall number of satellites depends sensitively on the reionization model, with a factor of 3-4 variation between the four models for a given host halo, although the difference is entirely in the population of faint satellites (M_V > -10). In addition, we find that for a given reionization model the total number of satellites differs by 10%-20% between the patchy and homogeneous scenarios, provided that the redshift is chosen appropriately for the instantaneous case. However, the halo-halo scatter from the six Aquarius halos is large, up to a factor of 2-3, and so is comparable to the difference between reionization scenarios. In order to use the population of faint dwarf galaxies around the Milky Way as a probe of the local reionization history, then, it is necessary to first better understand the general distribution of substructure around Milky Way-mass halos.

preprint2012arXiv

The Origin of the Microlensing Events Observed Towards the LMC and the Stellar Counterpart of the Magellanic Stream

We introduce a novel theoretical model to explain the long-standing puzzle of the nature of the microlensing events reported towards the Large Magellanic Cloud (LMC) by the MACHO and OGLE collaborations. We propose that a population of tidally stripped stars from the Small Magellanic Cloud (SMC) located ~4-10 kpc behind a lensing population of LMC disk stars can naturally explain the observed event durations (17-71 days), event frequencies and spatial distribution of the reported events. Differences in the event frequencies reported by the OGLE (~0.33 /yr) and MACHO (~1.75 /yr) surveys appear to be naturally accounted for by their different detection efficiencies and sensitivity to faint sources. The presented models of the Magellanic System were constructed without prior consideration of the microlensing implications. These results favor a scenario for the interaction history of the Magellanic Clouds, wherein the Clouds are on their first infall towards the Milky Way and the SMC has recently collided with the LMC 100-300 Myr ago, leading to a large number of faint sources distributed non-uniformly behind the LMC disk. In contrast to self-lensing models, microlensing events are also expected to occur in fields off the LMC's stellar bar since the stellar debris is not expected to be concentrated in the bar region. This scenario leads to a number of observational tests: the sources are low-metallicity SMC stars, they exhibit high velocities relative to LMC disk stars that may be detectable via proper motion studies, and, most notably, there should exist a stellar counterpart to the gaseous Magellanic Stream and Bridge with a V-band surface brightness of > 32 mag/arcsec^2. In particular, the stellar Bridge should contain enough RR Lyrae stars to be detected by the ongoing OGLE survey of this region.

preprint2012arXiv

The Role of Dwarf Galaxy Interactions in Shaping the Magellanic System and Implications for Magellanic Irregulars

We present a novel pair of numerical models of the interaction history between the Large and Small Magellanic Clouds (LMC and SMC, respectively) and our Milky Way (MW) in light of recent high precision proper motions (Kallivayalil et al. 2006a,b). Given the new velocities, cosmological simulations of structure formation favor a scenario where the Magellanic Clouds (MCs) are currently on their first infall towards our Galaxy (Boylan-Kolchin et al. 2011, Busha et al. 2011). We illustrate here that the observed irregular morphology and internal kinematics of the MCs (in gas and stars) are naturally explained by interactions between the LMC and SMC, rather than gravitational interactions with the MW. This picture further supports a first infall scenario (Besla et a. 2007). In particular, we demonstrate that the Magellanic Stream, a band of HI gas trailing behind the MCs 150 degrees across the sky, can be accounted for by the action of LMC tides on the SMC before the system was accreted by the MW. We further demonstrate that the off-center, warped stellar bar of the LMC and its one-armed spiral, can be naturally explained by a recent direct collision with the SMC. Such structures are key morphological characteristics of a class of galaxies referred to as Magellanic Irregulars (de Vaucouleurs & Freeman 1972), the majority of which are not associated with massive spiral galaxies. We infer that dwarf-dwarf galaxy interactions are important drivers for the morphological evolution of Magellanic Irregulars and can dramatically affect the efficiency of baryon removal from dwarf galaxies via the formation of extended tidal bridges and tails. Such interactions are important not only for the evolution of dwarf galaxies but also have direct consequences for the buildup of baryons in our own MW, as LMC-mass systems are believed to be the dominant building blocks of MW-type halos.

preprint2012arXiv

Why Are AGN and Host Galaxies Misaligned?

It is well-established observationally that the characteristic angular momentum axis on small scales around AGN, traced by radio jets and the putative torus, is not well-correlated with the large-scale angular momentum axis of the host galaxy. In this paper, we show that such misalignments arise naturally in high-resolution simulations in which we follow angular momentum transport and inflows from galaxy to sub-pc scales near AGN, triggered either during galaxy mergers or by instabilities in isolated disks. Sudden misalignments can sometimes be caused by single massive clumps falling into the center slightly off-axis, but more generally, they arise even when the gas inflows are smooth and trace only global gravitational instabilities. When several nested, self-gravitating modes are present, the inner ones can precess and tumble in the potential of the outer modes. Resonant angular momentum exchange can flip or re-align the spin of an inner mode on a short timescale, even without the presence of massive clumps. We therefore do not expect that AGN and their host galaxies will be preferentially aligned, nor should the relative alignment be an indicator of the AGN fueling mechanism. We discuss implications of this conclusion for AGN feedback and BH spin evolution. The misalignments may mean that even BHs accreting from smooth large-scale disks will not be spun up to maximal rotation, and so have more modest radiative efficiencies and inefficient jet formation. Even more random orientations are possible if there is further, un-resolved clumpiness in the gas, and more ordered accretion may occur if the inflow is slower and not self-gravitating.

preprint2011arXiv

Dynamics of the Magellanic Clouds in a LCDM Universe

We examine Milky Way-Magellanic Cloud systems selected from the Millennium-II Simulation in order to place the orbits of the Magellanic Clouds in a cosmological context. Our analysis shows that satellites massive enough to be LMC analogs are typically accreted at late times. Moreover, those that are accreted at early times and survive to the present have orbital properties that are discrepant with those observed for the LMC. The high velocity of the LMC, coupled with the dearth of unbound orbits seen in the simulation, argues that the mass of the MW's halo is unlikely to be less than 2 x 10^12 Msun. This conclusion is further supported by statistics of halos hosting satellites with masses, velocities, and separations comparable to those of the LMC. We further show that: (1) LMC and SMC-mass objects are not particularly uncommon in MW-mass halos; (2) the apparently high angular momentum of the LMC is not cosmologically unusual; and (3) it is rare for a MW halo to host a LMC-SMC binary system at z=0, but high speed binary pairs accreted at late times are possible. Based on these results, we conclude that the LMC was accreted within the past four Gyr and is currently making its first pericentric passage about the MW.

preprint2011arXiv

K+A Galaxies as the Aftermath of Gas-Rich Mergers: Simulating the Evolution of Galaxies as Seen by Spectroscopic Surveys

Models of poststarburst (or "K+A") galaxies are constructed by combining fully three-dimensional hydrodynamic simulations of galaxy mergers with radiative transfer calculations of dust attenuation. Spectral line catalogs are generated automatically from moderate-resolution optical spectra calculated as a function of merger progress in each of a large suite of simulations. The mass, gas fraction, orbital parameters, and mass ratio of the merging galaxies are varied systematically, showing that the lifetime and properties of the K+A phase are strong functions of merger scenario. K+A durations are generally less than ~0.1-0.3 Gyr, significantly shorter than the commonly assumed 1 Gyr, which is obtained only in rare cases, owing to a wide variation in star formation histories resulting from different orbital and progenitor configurations. Combined with empirical merger rates, the model lifetimes predict rapidly-rising K+A fractions as a function of redshift that are consistent with results of large spectroscopic surveys, resolving tension between the observed K+A abundance and that predicted when one assumes the K+A duration is the lifetime of A stars (~1 Gyr). The effects of dust attenuation, viewing angle, and aperture bias on our models are analyzed. In some cases, the K+A features are longer-lived and more pronounced when AGN feedback removes dust from the center, uncovering the young stars formed during the burst. In this picture, the K+A phase begins during or shortly after the bright starburst/AGN phase in violent mergers, and thus offers a unique opportunity to study the effects of quasar and star formation feedback on the gas reservoir and evolution of the remnant. Analytic fitting formulae are provided for the estimates of K+A incidence as a function of merger scenario.

preprint2011arXiv

Recoiling Black Holes in Merging Galaxies

Gravitational-wave (GW) recoil of merging supermassive black holes (SMBHs) may influence the co-evolution of SMBHs and their host galaxies. We examine this possibility using SPH/N-body simulations of gaseous galaxy mergers in which the merged BH receives a recoil kick. With our suite of over 200 merger simulations, we identify systematic trends in the behavior of recoiling BHs. Our main results are as follows. (1) While BHs kicked at nearly the central escape speed (vesc) are essentially "lost" to the galaxy, in gas rich mergers, BHs kicked with up to about 0.7 vesc may be confined to the central few kpc of the galaxy. (2) The inflow of cold gas during a gas-rich major merger may cause a rapid increase in central escape speed; in such cases recoil trajectories will depend on the timing of the BH merger relative to the change in vesc. (3) Recoil events generally reduce the lifetimes of bright active galactic nuclei (AGN) but may actually extend AGN lifetimes at lower luminosities. (4) Recoiling AGN may be observable via kinematic offsets (v > 500 km s^-1) or spatial offsets (R > 1 kpc) for lifetimes of up to about 10 - 100 Myr. (5) Rapidly-recoiling BHs may be up to about 5 times less massive than their stationary counterparts. These mass deficits lower the normalization of the M - sigma relation and contribute to both intrinsic and overall scatter. (6) Finally, the displacement of AGN feedback by a recoil event causes higher central star formation rates in the merger remnant, thereby extending the starburst phase of the merger and creating a denser, more massive stellar cusp.

preprint2011arXiv

Recoiling Black Holes in Merging Galaxies: Relationship to AGN Lifetimes and Merger Remnant Properties

Central supermassive black holes (SMBHs) are a ubiquitous feature of locally-observed galaxies, and ample evidence suggests that the growth of SMBHs and their host galaxies is closely linked. However, in the event of a merger, gravitational-wave (GW) recoil may displace a SMBH from its galactic center, or eject it entirely. To explore the consequences of this phenomenon, we use hydrodynamic simulations of gaseous galaxy mergers that include a range of BH recoil velocities. We have generated a suite of over 200 simulations with more than 60 merger models, enabling us to identify systematic trends in the behavior of recoiling BHs -- specifically (i) their dynamics, (ii) their observable signatures, and (iii) their effects on BH/galaxy co-evolution. (i) Recoiling BH trajectories depend heavily on the gas content of the host galaxy; maximal BH displacements from the center may vary by up to an order of magnitude between gas-rich and gas-poor mergers. In some cases, recoil trajectories also depend on the timing of the BH merger relative to the formation of the galaxy merger remnant. (ii) Recoiling BHs may be observable as offset active galactic nuclei (AGN) via either kinematic offsets (v > 800 km s^-1) or spatial offsets (R > 1 kpc) for lifetimes of about 1 - 100 Myr. In addition, recoil events affect the total AGN lifetime. GW recoil generally reduces the lifetimes of bright AGN, but may extend lower-luminosity AGN lifetimes. (iii) Rapidly-recoiling BHs may be up to about 5 times less massive than their stationary counterparts. These mass deficits lower the normalization of the M - sigma relation and contribute to both intrinsic and overall scatter. Furthermore, recoil events displace AGN feedback from the galactic center, which enhances central star formation rates. This results in longer starburst phases and higher central stellar densities in merger remnants.

preprint2011arXiv

Relation Between Globular Clusters and Supermassive Black Holes in Ellipticals as a Manifestation of the Black Hole Fundamental Plane

We analyze the relation between the mass of the central supermassive black hole (Mbh) and the number of globular clusters (Ngc) in elliptical galaxies and bulges as a ramification of the black hole fundamental plane, the theoretically predicted and observed multi-variable correlation between Mbh and bulge binding energy. Although the tightness of the Mbh-Ngc correlation suggests an unlikely causal link between supermassive black holes and globular clusters, such a correspondence can exhibit small scatter even if the physical relationship is indirect. We show that the relatively small scatter of the Mbh-Ngc relation owes to the mutual residual correlation of Mbh and Ngc with stellar mass when the velocity dispersion is held fixed. Thus, present observations lend evidence for feedback-regulated models in which the bulge binding energy is most important; they do not necessarily imply any `special' connection between globular clusters and Mbh. This raises the question of why Ngc traces the formation of ellipticals and bulges sufficiently well to be correlated with binding energy.

preprint2011arXiv

The CO-H2 Conversion Factor in Disc Galaxies and Mergers

Relating the observed CO emission from giant molecular clouds (GMCs) to the underlying H2 column density is a long-standing problem in astrophysics. While the Galactic CO-H2 conversion factor (Xco) appears to be reasonably constant, observations indicate that Xco may be depressed in high-surface density starburst environments. Using a multi-scale approach, we investigate the dependence of Xco on the galactic environment in numerical simulations of disc galaxies and galaxy mergers. Xco is proportional to the GMC surface density divided by the integrated CO intensity, Wco, and Wco is related to the kinetic temperature and velocity dispersion in the cloud. In disc galaxies (except within the central ~ kpc), the galactic environment is largely unimportant in setting the physical properties of GMCs provided they are gravitationally bound. The temperatures are roughly constant at ~10 K due to the balance of CO cooling and cosmic ray heating, giving a nearly constant CO-H2 conversion factor in discs. In mergers, the velocity dispersion of the gas rises dramatically during coalescence. The gas temperature also rises as it couples well to the warm (~50 K) dust at high densities (n > 10^4 cm^-3). The rise in velocity dispersion and temperature combine to offset the rise in surface density in mergers, causing Xco to drop by a factor of ~2-10 compared to the disc simulation. This model predicts that high-resolution ALMA observations of nearby ULIRGs should show velocity dispersions of ~10-100 km/s, and brightness temperatures comparable to the dust temperatures.

preprint2011arXiv

The Kennicutt-Schmidt Star Formation Relation at z~2

Recent observations of excited CO emission lines from z~2 disc galaxies have shed light on the Kennicutt-Schmidt relation at high-z via observed SFR-CO (J=2-1) and (J=3-2) relations. Here, we describe a novel methodology for utilising these observations of high-excitation CO to derive the underlying Schmidt (SFR-rho^N) relationship. To do this requires an understanding of the potential effects of differential CO excitation with SFR. If the most heavily star-forming galaxies have a larger fraction of their gas in highly excited CO states than the lower SFR galaxies, then the observed molecular SFR-CO^alpha index, alpha, will be less than the underlying (volumetric) Schmidt index, N. Utilising a combination of SPH models of galaxy evolution and molecular line radiative transfer, we present the first calculations of CO excitation in z~2 disc galaxies with the aim of developing a mapping between various observed SFR-CO relationships and the underlying Schmidt relation. We find that even in relatively luminous z~2 discs, differential excitation does indeed exist, resulting in alpha < N for highly excited CO lines. This means that an observed (e.g.) SFR-CO (J=3-2) relation does not map linearly to SFR-H2 relation. We utilise our model results to provide a mapping from alpha to N for the range of Schmidt indices N=1-2. By comparing to recent observational surveys, we find that the observed SFR-CO (J=2-1) and SFR-CO (J=3-2) relations suggest that an underlying SFR rho^1.5 relation describes z~2 disc galaxies.

preprint2011arXiv

The Metallicity Evolution of Interacting Galaxies

Nuclear inflows of metal-poor interstellar gas triggered by galaxy interactions can account for the systematically lower central oxygen abundances observed in local interacting galaxies. Here, we investigate the metallicity evolution of a large set of simulations of colliding galaxies. Our models include cooling, star formation, feedback, and a new stochastic method for tracking the mass recycled back to the interstellar medium from stellar winds and supernovae. We study the influence of merger-induced inflows, enrichment, gas consumption, and galactic winds in determining the nuclear metallicity. The central metallicity is primarily a competition between the inflow of low-metallicity gas and enrichment from star formation. An average depression in the nuclear metallicity of ~0.07 is found for gas-poor disk-disk interactions. Gas-rich disk-disk interactions, on the other hand, typically have an enhancement in the central metallicity that is positively correlated with the gas content. The simulations fare reasonably well when compared to the observed mass-metallicity and separation-metallicity relationships, but further study is warranted.

preprint2011arXiv

The Origins of AGN Obscuration: The 'Torus' as a Dynamical, Unstable Driver of Accretion

Multi-scale simulations have made it possible to follow gas inflows onto massive black holes (BHs) from galactic scales to the accretion disk. When sufficient gas is driven towards the BH, gravitational instabilities generically form lopsided, eccentric disks that propagate inwards. The lopsided stellar disk exerts a strong torque on the gas disk, driving inflows that fuel rapid BH growth. Here, we investigate whether the same gas disk is the 'torus' invoked to explain obscured AGN. The disk is generically thick and has characteristic ~1-10 pc sizes and masses resembling those required of the torus. The scale heights and obscured fractions of the predicted torii are substantial even in the absence of strong stellar feedback providing the vertical support. Rather, they can be maintained by strong bending modes and warps excited by the inflow-generating instabilities. Other properties commonly attributed to feedback processes may be explained by dynamical effects: misalignment between torus and host galaxy, correlations between local SFR and turbulent gas velocities, and dependence of obscured fractions on AGN luminosity or SFR. We compare the predicted torus properties with observations of gas surface density profiles, kinematics, scale heights, and SFR densities in AGN nuclei, and find that they are consistent. We argue that it is not possible to reproduce these observations and the observed column density (N_H) distribution without a clumpy gas distribution, but allowing for clumping on small scales the predicted N_H distribution is in good agreement with observations from 10^20-27 cm^-2. We examine how N_H scales with galaxy and AGN properties, and find that AGN feedback may be necessary to explain some trends with luminosity and/or redshift. The torus is not merely a bystander or passive fuel source for accretion, but is itself the mechanism driving accretion.

preprint2011arXiv

The SEDs of Interacting Galaxies

The evolution of galaxies is greatly influenced by their interactions. As part of a program to study interacting galaxies, we have measured and modeled the spectral energy distri- butions (SEDs) from the ultraviolet (UV) to the far-infrared (FIR). We describe the constraints imposed on star formation histories by these SEDs, and the variations therein seen across the interaction sequence, and we compare the results of different star formation rate prescriptions applied to the data. The sample itself is based on the Spitzer Interacting Galaxy Survey (SIGS) of 111 galaxies in 50 systems, a project designed to probe a range of galaxy interaction parameters in the infrared. Our SEDs combine the Spitzer results with multiwavelength data from other missions, in particular GALEX and Herschel. The subset presented here is the sample for which FIR Herschel observations are currently publicly available.

preprint2011arXiv

The Signatures of Large-scale Temperature and Intensity Fluctuations in the Lyman-alpha Forest

It appears inevitable that reionization processes would have produced large-scale temperature fluctuations in the intergalactic medium. Using toy temperature models and detailed heating histories from cosmological simulations of HeII reionization, we study the consequences of inhomogeneous heating for the Ly-alpha forest. The impact of temperature fluctuations in physically well-motivated models can be surprisingly subtle. In fact, we show that temperature fluctuations at the level predicted by our reionization simulations do not give rise to detectable signatures in the types of statistics that have been employed previously. However, because of the aliasing of small-scale density power to larger scale modes in the line-of-sight Ly-alpha forest power spectrum, earlier analyses were not sensitive to 3D modes with >~ 30 comoving Mpc wavelengths -- scales where temperature fluctuations are likely to be relatively largest. The ongoing Baryon Oscillation Spectroscopic Survey (BOSS) aims to measure the 3D power spectrum of the Ly-alpha forest, P_F, from a large sample of quasars in order to avoid this aliasing. We find that physically motivated temperature models can alter P_F at an order unity level at k <~ 0.1 comoving Mpc^{-1}, a magnitude that should be easily detectable with BOSS. Fluctuations in the intensity of the ultraviolet background can also alter P_F significantly. These signatures will make it possible for BOSS to study the thermal impact of HeII reionization at 2 < z < 3 and to constrain models for the sources of the ionizing background. Future spectroscopic surveys could extend this measurement to even higher redshifts, potentially detecting the thermal imprint of hydrogen reionization.

preprint2011arXiv

What Does a Submillimeter Galaxy Selection Actually Select? The Dependence of Submillimeter Flux Density on Star Formation Rate and Dust Mass

We perform 3-D dust radiative transfer (RT) calculations on hydrodynamic simulations of isolated and merging disk galaxies in order to quantitatively study the dependence of observed-frame submillimeter (submm) flux density on galaxy properties. We find that submm flux density and star formation rate (SFR) are related in dramatically different ways for quiescently star-forming galaxies and starbursts. Because the stars formed in the merger-induced starburst do not dominate the bolometric luminosity and the rapid drop in dust mass and more compact geometry cause a sharp increase in dust temperature during the burst, starbursts are very inefficient at boosting submm flux density (e.g., a $\ga16$x boost in SFR yields a $\la 2$x boost in submm flux density). Moreover, the ratio of submm flux density to SFR differs significantly between the two modes; thus one cannot assume that the galaxies with highest submm flux density are necessarily those with the highest bolometric luminosity or SFR. These results have important consequences for the bright submillimeter-selected galaxy (SMG) population. Among them are: 1. The SMG population is heterogeneous. In addition to merger-driven starbursts, there is a subpopulation of galaxy pairs, where two disks undergoing a major merger but not yet strongly interacting are blended into one submm source because of the large ($\ga 15$", or $\sim 130$ kpc at $z = 2$) beam of single-dish submm telescopes. 2. SMGs must be very massive ($M_{\star} \ga 6 \times 10^{10} \msun$). 3. The infall phase makes the SMG duty cycle a factor of a few greater than what is expected for a merger-driven starburst. (Abridged.)

preprint2010arXiv

A Non-Parametric Estimate of Mass 'Scoured' in Galaxy Cores

We present a simple estimate of the mass 'deficits' in cored spheroids, as a function of galaxy mass and radius within the galaxy. Previous attempts to measure such deficits depended on fitting some functional form to the profile at large radii and extrapolating inwards; this is sensitive to the assumed functional form and does not allow for variation in nuclear profile shapes. We take advantage of larger data sets to directly construct stellar mass profiles of observed systems and measure the stellar mass enclosed in a series of physical radii (M(<R)), for samples of cusp and core spheroids at the same stellar mass. There is a significant bimodality in this distribution at small radii, and we non-parametrically measure the median offset between core and cusp populations (the deficit Delta_M(<R)). We construct the scoured mass profile as a function of radius, without reference to any assumed functional form. The mass deficit rises in power-law fashion (Delta_M(<R) R^{1.3-1.8}) from a significant but small mass at R<10pc, to asymptote to a maximum ~0.5-2 M_BH at ~100pc. At larger radii there is no statistically significant separation between populations; the upper limit to the cumulative scoured mass at ~kpc is ~2-4 M_BH. This does not depend strongly on stellar mass. The dispersion in M(<R) appears larger in the core population, possibly reflecting the fact that scouring increases the scatter in profile shapes. These results are in good agreement with models of scouring from BH binary systems.

preprint2010arXiv

A Physical Model for z~2 Dust Obscured Galaxies

We present a physical model for the origin of z~2 Dust-Obscured Galaxies (DOGs), a class of high-redshift ULIRGs selected at 24 micron which are particularly optically faint (24/R>1000). By combining N-body/SPH simulations of high redshift galaxy evolution with 3D polychromatic dust radiative transfer models, we find that luminous DOGs (with F24 > 0.3 mJy at z~2 are well-modeled as extreme gas-rich mergers in massive (~5x10^12-10^13 Msun) halos, with elevated star formation rates (~500-1000 Msun/yr) and/or significant AGN growth (Mdot > 0.5 Msun/yr), whereas less luminous DOGs are more diverse in nature. At final coalescence, merger-driven DOGs transition from being starburst dominated to AGN dominated, evolving from a "bump" to a power-law shaped mid-IR (IRAC) spectral energy distribution (SED). After the DOG phase, the galaxy settles back to exhibiting a "bump" SED with bluer colors and lower star formation rates. While canonically power-law galaxies are associated with being AGN-dominated, we find that the power-law mid-IR SED can owe both to direct AGN contribution, as well as to a heavily dust obscured stellar bump at times that the galaxy is starburst dominated. Thus power-law galaxies can be either starburst or AGN dominated. Less luminous DOGs can be well-represented either by mergers, or by massive ($M_{\rm baryon} ~5x10^11 Msun) secularly evolving gas-rich disc galaxies (with SFR > 50 Msun/yr). By utilising similar models as those employed in the SMG formation study of Narayanan et al. (2010), we investigate the connection between DOGs and SMGs. We find that the most heavily star-forming merger driven DOGs can be selected as Submillimetre Galaxies (SMGs), while both merger-driven and secularly evolving DOGs typically satisfy the BzK selection criteria.

preprint2010arXiv

Cold Fronts by Merging of Shocks

Cold fronts (CFs) are found in most galaxy clusters, as well as in some galaxies and groups of galaxies. We propose that some CFs are relics of merging between two shocks propagating in the same direction. Such shock mergers typically result in a quasi-spherical, factor ~1.4-2.7 discontinuity in density and in temperature. These CFs may be found as far out as the virial shock, unlike what is expected in other CF formation models. As a demonstration of this effect, we use one dimensional simulations of clusters and show that shock induced cold fronts form when perturbations such as explosions or mergers occur near the cluster's centre. Perturbations at a cluster's core induce periodic merging between the virial shock and outgoing secondary shocks. These collisions yield a distinctive, concentric, geometric sequence of CFs which trace the expansion of the virial shock.

preprint2010arXiv

Cold Fronts from Shock Collisions

Cold fronts (CFs) are found in most galaxy clusters, as well as in some galaxies and groups of galaxies. We propose that some CFs are relics of collisions between trailing shocks. Such a collision typically results in a spherical, factor ~1.4-2.7 density/temperature discontinuity. These CFs may be found as far as the virial shock, unlike in other CF formation models. As a demonstration of this effect, we use one dimensional simulations where halo reverberations involving periodic collisions between the virial shock and outgoing secondary shocks exist. These collisions yield a distinctive, concentric geometric sequence of CFs which trace the expansion of the virial shock.

preprint2010arXiv

Constraints on Black Hole Growth, Quasar Lifetimes, and Eddington Ratio Distributions from the SDSS Broad Line Quasar Black Hole Mass Function

We present an estimate of the black hole mass function (BHMF) of broad line quasars (BLQSOs) that self-consistently corrects for incompleteness and the statistical uncertainty in the mass estimates, based on a sample of 9886 quasars at 1 < z < 4.5 drawn from the Sloan Digital Sky Survey. We find evidence for `cosmic downsizing' of black holes in BLQSOs, where the peak in their number density shifts to higher redshift with increasing black hole mass. The cosmic mass density for black holes seen as BLQSOs peaks at z ~ 2. We estimate the completeness of the SDSS as a function of black hole mass and Eddington ratio, and find that at z > 1 it is highly incomplete at M_BH < 10^9 M_Sun and L / L_Edd < 0.5. We also estimate a lower limit on the lifetime of a single BLQSO phase and we place constraints on the maximum mass of a black hole in a BLQSO. Our estimated distribution of BLQSO Eddington ratios peaks at L / L_Edd ~ 0.05 and has a dispersion of ~ 0.4 dex, implying that most BLQSOs are not radiating at or near the Eddington limit; however the location of the peak is subject to considerable uncertainty. The steep increase in number density of BLQSOs toward lower Eddington ratios is expected if the BLQSO accretion rate monotonically decays with time. Furthermore, our estimated lifetime and Eddington ratio distributions imply that the majority of the most massive black holes spend a significant amount of time growing in an earlier obscured phase, a conclusion which is independent of the unknown obscured fraction. These results are consistent with models for self-regulated black hole growth, at least for massive systems at z > 1, where the BLQSO phase occurs at the end of a fueling event when black hole feedback unbinds the accreting gas, halting the accretion flow.

preprint2010arXiv

Cosmic Rays from Pulsars and Magnetars

We compare the expected abundance of cosmic-ray electrons and positrons from pulsars and magnetars. We assume that the distribution of infant pulsars and magnetars follows that of high-mass stars in the Milky Way and that the production rate of cosmic rays is proportional to the spin-down and magnetic-decay power of pulsars and magnetars, respectively. In combination with primary and secondary cosmic-ray leptons from other sources (especially supernova remnants), we find that both magnetars and pulsars can easily account for the observed cosmic-ray spectrum, in particular the dip seen by HESS at several TeV and the increase in positron fraction found by PAMELA.

preprint2010arXiv

Cross-Correlation between Damped Ly-alpha Systems and Ly-break Galaxies in Cosmological SPH Simulations

We calculate the cross-correlation function (CCF) between damped Ly-a systems (DLAs) and Lyman break galaxies (LBGs) using cosmological hydrodynamic simulations at z=3. We compute the CCF with two different methods. First, we assume that there is one DLA in each dark matter halo if its DLA cross section is non-zero. In our second approach we weight the pair-count by the DLA cross section of each halo, yielding a cross-section-weighted CCF. We also compute the angular CCF for direct comparison with observations. Finally, we calculate the auto-correlation functions of LBGs and DLAs, and their bias against the dark matter distribution. For these different approaches, we consistently find that there is good agreement between our simulations and observational measurements by Cooke et al. and Adelberger et al. Our results thus confirm that the spatial distribution of LBGs and DLAs can be well described within the framework of the concordance Lambda CDM model. We find that the correlation strengths of LBGs and DLAs are consistent with the actual observations, and in the case of LBGs it is higher than would be predicted by low-mass galaxy merger models.

preprint2010arXiv

Do sub-millimeter galaxy number counts provide evidence for a top-heavy IMF?

Theoretical models have had difficulty matching the observed number density of sub-millimeter galaxies (SMGs), causing some authors (e.g., Baugh et al. 2005) to suggest that SMGs provide evidence for a top-heavy initial mass function (IMF). To test this claim, we have, for the first time, combined high-resolution 3-D hydrodynamic simulations of isolated and merging massive, gas-rich galaxies, radiative transfer, and a semi-empirical merger rate model to predict the number density of SMGs. Our model can reproduce the observed SMG number density even when using a standard (Kroupa) IMF. The agreement is due to a combination of relatively long sub-mm duty cycles for mergers (a few times 10^8 years for our most massive models), which owe to our combination of high-resolution 3-D hydrodynamic simulations and dust radiative transfer; sufficient number densities of massive, gas-rich mergers; and the decrease in sub-mm counts observed by recent deep/wide surveys (e.g., Austermann et al. 2010) relative to previous surveys. Our results suggest that the observed SMG number counts do not provide evidence for a top-heavy IMF at high redshift.

preprint2010arXiv

Lensing reconstruction of cluster-mass cross-correlation with cosmic microwave background polarization

We extend our maximum likelihood method for reconstructing the cluster-mass cross-correlation from cosmic microwave background (CMB) temperature anisotropies and develop new estimators that utilize six different quadratic combinations of CMB temperature and polarization fields. Our maximum likelihood estimators are constructed with delensed CMB temperature and polarization fields by using an assumed model of the convergence field and they can be iteratively applied to a set of clusters, approaching to the optimal condition for the lensing reconstruction as the assumed initial model is refined. Using smoothed particle hydrodynamics simulations, we create a catalog of realistic clusters obtainable from the current Sunyaev-Zel'dovich (SZ) surveys, and we demonstrate the ability of the maximum likelihood estimators to reconstruct the cluster-mass cross-correlation from the massive clusters. The iTT temperature estimator provides a signal-to-noise ratio of a factor 3 larger than the iEB polarization estimator, unless the detector noise for measuring polarization anisotropies is controlled under 3 microK.

preprint2010arXiv

Lyman-alpha Emitters and Lyman-break Galaxies at z=3-6 in Cosmological SPH Simulations

We study the properties of Lyman-alpha emitters (LAEs) and Lyman-break galaxies (LBGs) at z=3-6 using cosmological SPH simulations. We investigate two simple scenarios for explaining the observed Ly-a and rest-frame UV luminosity functions (LFs) of LAEs: (i) the "escape fraction" scenario, in which the "effective" escape fraction (including the IGM attenuation) of Ly-a photons is f_Lya ~0.1 (0.15) at z=3 (6), and (ii) the "stochastic" scenario, in which the fraction of LAEs that are turned on at z=3 (6) is \Cstoc ~0.07 (0.2) after correcting for the IGM attenuation. Our comparisons with a number of different observations suggest that the stochastic scenario is preferred over the escape fraction scenario. We find that the mean values of stellar mass, metallicity and black hole mass hosted by LAEs are all smaller in the stochastic scenario than in the escape fraction scenario. In our simulations, the galaxy stellar mass function evolves rapidly, as expected in hierarchical structure formation. However, its evolution is largely compensated by a beginning decline in the specific star formation rate, resulting in little evolution of the rest-frame UV LF from z=6 to 3. The rest-frame UV LF of both LAEs and LBGs at z=3 & 6 can be described well by the stochastic scenario provided the extinction is moderate, E(B-V) ~0.15, for both populations, although our simulation might be overpredicting the number of bright LBGs at z=6. We also discuss the correlation function and bias of LAEs. The Ly-a LFs at z=6 in a field-of-view of 0.2 deg^2 show a significantly larger scatter owing to cosmic variance relative to that in a 1 deg^2 field, and the scatter seen in the current observational estimates of the Ly-a LF can be accounted for by cosmic variance.

preprint2010arXiv

Mergers and Bulge Formation in Lambda-CDM: Which Mergers Matter?

We use a suite of semi-empirical models to predict galaxy merger rates and contributions to bulge growth as functions of merger mass, redshift, and mass ratio. The models use empirical halo occupation constraints to identify mergers, together with high-resolution simulations to quantify how mergers with different properties contribute to the bulge population. We find good agreement with a variety of observational constraints, and provide fitting functions for merger rates and contributions to bulge growth. We identify several robust conclusions. (1) Major mergers dominate formation and assembly of L* bulges and the spheroid mass density, minor mergers contribute ~30%. (2) This is mass-dependent: bulge formation is dominated by more minor mergers in lower-mass systems. At higher masses, bulges form in major mergers near L*, but subsequently assemble in minor mergers. (3) The minor/major contribution is also morphology-dependent: higher B/T systems form in more major mergers, lower B/T systems form in situ from minor mergers. (4) Low-mass galaxies, being gas-rich, require more major mergers to reach the same B/T as high-mass systems. (5) Absolute merger rates increase with galaxy mass. (6) Predicted rates agree well with observations, but suggest that some morphology-selected samples include contamination from minor mergers. (7) Predicted rates agree with integrated growth in bulge mass with cosmic time, but with factor ~2 uncertainty - half the bulge mass density could come from non-mergers. We consider ~1000 model variations and quantify resulting uncertainties. Conclusions regarding the major/minor contribution to bulge growth are very robust, absolute merger rates have systematic factor ~2 uncertainties.

preprint2010arXiv

Mergers in Lambda-CDM: Uncertainties in Theoretical Predictions and Interpretations of the Merger Rate

Different methodologies lead to order-of-magnitude variations in predicted galaxy merger rates. We examine and quantify the dominant uncertainties. Different halo merger rates and subhalo 'destruction' rates agree to within a factor ~2 given proper care in definitions. If however (sub)halo masses are not appropriately defined or are under-resolved, the major merger rate can be dramatically suppressed. The dominant differences in galaxy merger rates owe to baryonic physics. Hydrodynamic simulations without feedback and older models that do not agree with the observed galaxy mass function propagate factor ~5 bias in the resulting merger rates. However, if the model matches the galaxy mass function, properties of central galaxies are sufficiently converged to give small differences in merger rates. But variations in baryonic physics of satellites also have dramatic effects. The known problem of satellite 'over-quenching' in most semi-analytic models (SAMs), whereby SAM satellites are too efficiently stripped of gas, could lead to order-of-magnitude under-estimates of merger rates for low-mass, gas-rich galaxies. Fixing the satellite properties to observations tends to predict higher merger rates, but with factor ~2 empirical uncertainties. Choice of mass ratio definition matters: at low masses, most true major mergers (in baryonic/dynamical galaxy mass) will appear to be minor mergers in their stellar or luminosity mass ratio. Observations and models using these criteria may underestimate major merger rates by factors ~5. Orbital parameters and gas fractions also introduce factor ~3 differences in amount of bulge formed by mergers, even for fixed mass ratio encounters.

preprint2010arXiv

Minimum Radii of Super-Earths: Constraints from Giant Impacts

The detailed interior structure models of super-Earth planets show that there is degeneracy in the possible bulk compositions of a super-Earth at a given mass and radius, determined via radial velocity and transit measurements, respectively. In addition, the upper and lower envelopes in the mass--radius relationship, corresponding to pure ice planets and pure iron planets, respectively, are not astrophysically well motivated with regard to the physical processes involved in planet formation. Here we apply the results of numerical simulations of giant impacts to constrain the lower bound in the mass--radius diagram that could arise from collisional mantle stripping of differentiated rocky/iron planets. We provide a very conservative estimate for the minimum radius boundary for the entire mass range of large terrestrial planets. This envelope is a readily testable prediction for the population of planets to be discovered by the Kepler mission.

preprint2010arXiv

On Sizes, Kinematics, M/L Gradients, and Light Profiles of Massive Compact Galaxies at z~2

We present a detailed analysis of the structure and resolved stellar populations of simulated merger remnants, and compare them to observations of compact quiescent galaxies at z ~ 2. We find that major merging is a viable mechanism to produce systems of ~ 10^11 Msun and ~ 1 kpc size, provided the gas fraction at the time of final coalescence is high (~ 40%), and provided that the progenitors are compact star-forming galaxies, as expected at high redshift. Their integrated spectral energy distributions and velocity dispersions are in good agreement with the observations, and their position in the (v_{maj}/sigma, ellipticity) diagram traces the upper envelope of the distribution of lower redshift early-type galaxies. The simulated merger remnants show time- and sightline-dependent M/L ratio gradients that result from a superposition of radially dependent stellar age, stellar metallicity, and extinction. The median ratio of effective radius in rest-frame V-band light to that in mass surface density is ~ 2 during the quiescent remnant phase. This is typically expressed by a negative color gradient (i.e., red core), which we expect to correlate with the integrated color of the system. Finally, the simulations differ from the observations in their surface brightness profile shape. The simulated remnants are typically best fit by high (n >> 4) Sersic indices, whereas observed quiescent galaxies at z ~ 2 tend to be less cuspy (median n ~ 2.3). Limiting early star formation in the progenitors may be required to prevent the simulated merger remnants from having extended wings.

preprint2010arXiv

Orbital structure of merger remnants: Trends with gas fraction in 1:1 mergers

Since the violent relaxation in hierarchical merging is incomplete, elliptical galaxies retain a wealth of information about their formation pathways in their present-day orbital structure. A variety of evidence indicates that gas-rich major mergers play an important role in the formation of elliptical galaxies. We simulate 1:1 disk mergers at seven different initial gas fractions ranging from 0 to 40%, using the TreeSPH code Gadget-2. We classify the stellar orbits in each remnant and construct radial profiles of the orbital content, intrinsic shape, and orientation. The dissipationless remnants are typically prolate-triaxial, dominated by box orbits within r_c ~ 1.5Reff, and by tube orbits in their outer parts. As the gas fraction increases, the box orbits within r_c are increasingly replaced by a population of short axis tubes (z-tubes) with near zero net rotation, and the remnants become progressively more oblate and round. The long axis tube (x-tube) orbits are highly streaming and relatively insensitive to the gas fraction, implying that their angular momentum is retained from the dynamically cold initial conditions. Outside r_c, the orbital structure is essentially unchanged by the gas. The 15-20% gas remnants often display disk-like kinematically distinct cores (KDCs). These remnants show an interesting resemblance, in both their velocity maps and intrinsic orbital structure, to the KDC galaxy NGC4365 (van den Bosch et al. 2008). At 30-40% gas, the remnants are rapidly rotating, with sharp embedded disks on ~ 1Reff scales. We predict a characteristic, physically intuitive orbital structure for 1:1 disk merger remnants, with a distinct transition between 1 and 3Reff that will be readily observable with combined data from the 2D kinematics surveys SAURON and SMEAGOL.

preprint2010arXiv

Quasi-Resonant Theory of Tidal Interactions

When a spinning system experiences a transient gravitational encounter with an external perturber, a quasi-resonance occurs if the spin frequency of the victim matches the peak orbital frequency of the perturber. Such encounters are responsible for the formation of long tails and bridges of stars during galaxy collisions. For high-speed encounters, the resulting velocity perturbations can be described within the impulse approximation. The traditional impulse approximation, however, does not distinguish between prograde and retrograde encounters, and therefore completely misses the resonant response. Here, using perturbation theory, we compute the effects of quasi-resonant phenomena on stars orbiting within a disk. Explicit expressions are derived for the velocity and energy change to the stars induced by tidal forces from an external gravitational perturber passing either on a straight line or parabolic orbit. Comparisons with numerical restricted three-body calculations illustrate the applicability of our analysis.

preprint2010arXiv

Recoiling Black Holes in Merging Galaxies: Relationship to AGN Lifetimes, Starbursts, and the M-sigma Relation

Gravitational-wave (GW) recoil of merging supermassive black holes (SMBHs) may influence the co-evolution of SMBHs and their host galaxies. We examine this possibility using SPH/N-body simulations of gaseous galaxy mergers in which the merged BH receives a recoil kick. This enables us to follow recoiling BHs in self-consistent, evolving merger remnants. In contrast to recent studies on similar topics, we conduct a large parameter study, generating a suite of over 200 simulations with more than 60 merger models and a range of recoil velocities (vk). Our main results are as follows. (1) BHs kicked at nearly the central escape speed (vesc) may oscillate on large orbits for up to a Hubble time, but in gas-rich mergers, BHs kicked with up to ~ 0.7 vesc may be confined to the central few kpc of the galaxy, owing to gas drag and steep central potentials. (2) vesc in gas-rich mergers may increase rapidly during final coalescence, in which case trajectories may depend on the timing of the BH merger relative to the formation of the potential well. (3) Recoil events generally reduce the lifetimes of bright active galactic nuclei (AGN), but may actually extend AGN lifetimes at lower luminosities. (4) Kinematically-offset AGN (v > 800 km s^-1) may be observable for up to ~ 10 Myr either immediately after the recoil or during pericentric passages through a gas-rich remnant. (5) Spatially-offset AGN (R > 1 kpc) generally have low luminosities and lifetimes of ~ 1 - 100 Myr. (6) Rapidly-recoiling BHs may be up to ~ 5 times less massive than their stationary counterparts. This lowers the normalization of the M-sigma relation and contributes to both intrinsic and overall scatter. (7) Finally, the displacement of AGN feedback after a recoil event enhances central star formation rates, thereby extending the starburst phase of the merger and creating a denser stellar cusp. [Abridged.]

preprint2010arXiv

Simulations and Analytic Calculations of Bubble Growth During Hydrogen Reionization

We present results from a large volume simulation of Hydrogen reionization. We combine 3d radiative transfer calculations and an N-body simulation, describing structure formation in the intergalactic medium, to detail the growth of HII regions around high redshift galaxies. Our simulation tracks 1024^3 dark matter particles, in a box of co-moving side length 65.6 Mpc/h. This large volume allows us to accurately characterize the size distribution of HII regions throughout most of the reionization process. At the same time, our simulation resolves many of the small galaxies likely responsible for reionization. It confirms a picture anticipated by analytic models: HII regions grow collectively around highly-clustered sources, and have a well-defined characteristic size, which evolves from a sub-Mpc scale at the beginning of reionization to R > 10 Mpc towards the end. We present a detailed statistical description of our results, and compare them with a numerical scheme based on the analytic model by Furlanetto, Zaldarriaga, and Hernquist. We find that the analytic calculation reproduces the size distribution of HII regions and the 21 cm power spectrum of the radiative transfer simulation remarkably well. The ionization field from the simulation, however, has more small scale structure than the analytic calculation, owing to Poisson scatter in the simulated abundance of galaxies on small scales. We propose and validate a simple scheme to incorporate this scatter into our calculations. Our results suggest that analytic calculations are sufficiently accurate to aid in predicting and interpreting the results of future 21 cm surveys. In particular, our fast numerical scheme is useful for forecasting constraints from future 21 cm surveys, and in constructing mock surveys to test data analysis procedures.

preprint2010arXiv

Simulations of the Magellanic Stream in a First Infall Scenario

Recent high precision proper motions from the Hubble Space Telescope (HST) suggest that the Large and Small Magellanic Clouds (LMC and SMC, respectively) are either on their first passage or on an eccentric long period (>6 Gyr) orbit about the Milky Way (MW). This differs markedly from the canonical picture in which the Clouds travel on a quasi-periodic orbit about the MW (period of ~2 Gyr). Without a short period orbit about the MW, the origin of the Magellanic Stream, a young (1-2 Gyr old) coherent stream of HI gas that trails the Clouds ~150 degrees across the sky, can no longer be attributed to stripping by MW tides and/or ram pressure stripping by MW halo gas. We propose an alternative formation mechanism in which material is removed by LMC tides acting on the SMC before the system is accreted by the MW. We demonstrate the feasibility and generality of this scenario using an N-body/SPH simulation with cosmologically motivated initial conditions constrained by the observations. Under these conditions we demonstrate that it is possible to explain the origin of the Magellanic Stream in a first infall scenario. This picture is generically applicable to any gas-rich dwarf galaxy pair infalling towards a massive host or interacting in isolation.

preprint2010arXiv

Spiral-Induced Star Formation in the Outer Disks of Galaxies

The outer regions of galactic disks have received increased attention since ultraviolet observations with GALEX demonstrated that nearly 30% of galaxies have UV emission beyond their optical extents, indicating star formation activity. These galaxies have been termed extended UV (XUV) disks. Here, we address whether these observations contradict the gas surface density threshold for star formation inferred from Halpha radial profiles of galaxies. We run smoothed particle hydrodynamics simulations of isolated disk galaxies with fiducial star formation prescriptions and show that over-densities owing to the presence of spiral structure can induce star formation in extended gas disks. For direct comparison with observations, we use the 3-D radiative transfer code Sunrise to create simulated FUV and K_s band images. We find that galaxies classified as Type I XUV disks are a natural consequence of spiral patterns, but we are unable to reproduce Type II XUV disks. We also compare our results to studies of the Kennicutt-Schmidt relation in outer disks.

preprint2009arXiv

A Measurement of Small Scale Structure in the 2.2 < z < 4.2 Lyman-alpha Forest

The amplitude of fluctuations in the Ly-a forest on small spatial scales is sensitive to the temperature of the IGM and its spatial fluctuations. The temperature of the IGM and its spatial variations contain important information about hydrogen and helium reionization. We present a new measurement of the small-scale structure in the Ly-a forest from 40 high resolution, high signal-to-noise, VLT spectra at z=2.2-4.2. We convolve each Ly-a forest spectrum with a suitably chosen wavelet filter, which allows us to extract the amount of small-scale structure in the forest as a function of position across each spectrum. We compare these measurements with high resolution hydrodynamic simulations of the Ly-a forest which track more than 2 billion particles. This comparison suggests that the IGM temperature close to the cosmic mean density (T_0) peaks near z=3.4, at which point it is greater than 20,000 K at 2-sigma confidence. The temperature at lower redshift is consistent with the fall-off expected from adiabatic cooling ($T_0 \propto (1+z)^2$), after the peak temperature is reached near z=3.4. At z=4.2 our results favor a temperature of T_0 = 15-20,000 K. However, owing mostly to uncertainties in the mean transmitted flux at this redshift, a cooler IGM model with T_0 = 10,000 K is only disfavored at the 2-sigma level here, although such cool IGM models are strongly discrepant with the z ~ 3-3.4 measurement. We do not detect large spatial fluctuations in the IGM temperature at any redshift covered by our data set. The simplest interpretation of our measurements is that HeII reionization completes sometime near z ~ 3.4, although statistical uncertainties are still large [Abridged].

preprint2009arXiv

A New Empirical Method to Infer the Starburst History of the Universe from Local Galaxy Properties

The centers of bulges are formed dissipationally via gas inflows over short timescales: the 'starburst' mode of star formation (SF). Recent work has shown that detailed observations can be used to separate the stellar mass profile of these 'burst relic' components in local systems. Together with the assumption that some Kennicutt-Schmidt law holds, and that the burst was indeed a dissipational gas-rich event, we show that the observed profiles can be inverted to obtain the time and space-dependent SF history of each burst. Performing this with a large sample of well-studied spheroids, we show that the implied bursts scale in magnitude, mass, peak SFR, and spatial extent with galaxy mass in simple manner, and provide fits to these correlations. Burst masses are ~10% the total spheroid mass; timescales a mass-independent ~10^8 yr; peak SFR ~M_burst/t_burst; and they decay in power-law fashion ~t^-2.4. Sizes are ~0.1 R_e(spheroid), but grow with time. Combined with measurements of the nuclear stellar population ages of these systems (i.e. burst times), it is possible to re-construct the burst contribution to the distribution of SFRs and IR luminosity functions at all redshifts. The burst LF agrees well with observed IR LFs at the brightest luminosities, at z=0-2. At low-L, bursts are always unimportant; the transition to their importance increases from ULIRG luminosities at z~0 to HyLIRG luminosities at z~2. At all redshifts, bursts are a small fraction (~5-10%) of the total SFR density. We discuss possible implications of tension between maximum relic stellar mass densities in massive systems, and estimated number counts of brightest sub-millimeter galaxies.

preprint2009arXiv

Are Most Low-Luminosity AGN Really Obscured?

At low Eddington ratio (mdot), two effects make it harder to detect AGN given some selection criteria. First, even with fixed accretion physics, AGN are diluted/less luminous relative to their hosts; the magnitude of this depends on host properties and so on luminosity and redshift. Second, they may transition to a radiatively inefficient state, changing SED shape and dramatically decreasing in optical/IR luminosity. These effects lead to differences in observed AGN samples, even at fixed bolometric luminosity and after correction for obscuration. The true Eddington ratio distribution may depend strongly on luminosity, but this will be seen only in surveys robust to dilution and radiative inefficiency (X-ray or narrow-line samples); selection effects imply that AGN in optical samples will have uniformly high mdot. This also implies that different selection methods yield systems with different hosts: the clustering of faint optical/IR sources will be weaker than that of X-ray sources, and optical/IR Seyferts will reside in more disk-dominated galaxies while X-ray selected Seyferts will preferentially occupy early-type systems. If observed mdot distributions are correct, a large fraction of low-luminosity AGN currently classified as 'obscured' are in fact diluted and/or radiatively inefficient, not obscured by gas or dust. This is equally true if X-ray hardness is used as a proxy for obscuration, since radiatively inefficient SEDs near mdot~0.01 are X-ray hard. These effects can explain most of the claimed luminosity/redshift dependence in the 'obscured' AGN population, with the true obscured fraction as low as 20%.

preprint2009arXiv

Discriminating Between the Physical Processes that Drive Spheroid Size Evolution

Massive galaxies at high-z have smaller effective radii than those today, but similar central densities. Their size growth therefore relates primarily to the evolving abundance of low-density material. Various models have been proposed to explain this evolution, which have different implications for galaxy, star, and BH formation. We compile observations of spheroid properties as a function of redshift and use them to test proposed models. Evolution in progenitor gas-richness with redshift gives rise to initial formation of smaller spheroids at high-z. These systems can then evolve in apparent or physical size via several channels: (1) equal-density 'dry' mergers, (2) later major or minor 'dry' mergers with less-dense galaxies, (3) adiabatic expansion, (4) evolution in stellar populations & mass-to-light-ratio gradients, (5) age-dependent bias in stellar mass estimators, (6) observational fitting/selection effects. If any one of these is tuned to explain observed size evolution, they make distinct predictions for evolution in other galaxy properties. Only model (2) is consistent with observations as a dominant effect. It is the only model which allows for an increase in M_BH/M_bulge with redshift. Still, the amount of merging needed is larger than that observed or predicted. We therefore compare cosmologically motivated simulations, in which all these effects occur, & show they are consistent with all the observational constraints. Effect (2), which builds up an extended low-density envelope, does dominate the evolution, but effects 1,3,4, & 6 each contribute ~20% to the size evolution (a net factor ~2). This naturally also predicts evolution in M_BH-sigma similar to that observed.

preprint2009arXiv

Mergers, AGN, and 'Normal' Galaxies: Contributions to the Distribution of Star Formation Rates and Infrared Luminosity Functions

We use a novel method to predict the contribution of normal star-forming galaxies, merger-induced bursts, and obscured AGN, to IR luminosity functions (LFs) and global SFR densities. We use empirical halo occupation constraints to populate halos with galaxies and determine the distribution of normal and merging galaxies. Each system can then be associated with high-resolution hydrodynamic simulations. We predict the distribution of observed luminosities and SFRs, from different galaxy classes, as a function of redshift from z=0-6. We provide fitting functions for the predicted LFs, quantify the uncertainties, and compare with observations. At all redshifts, 'normal' galaxies dominate the LF at moderate luminosities ~L* (the 'knee'). Merger-induced bursts increasingly dominate at L>>L*; at the most extreme luminosities, AGN are important. However, all populations increase in luminosity at higher redshifts, owing to increasing gas fractions. Thus the 'transition' between normal and merger-dominated sources increases from the LIRG-ULIRG threshold at z~0 to bright Hyper-LIRG thresholds at z~2. The transition to dominance by obscured AGN evolves similarly, at factor of several higher L_IR. At all redshifts, non-merging systems dominate the total luminosity/SFR density, with merger-induced bursts constituting ~5-10% and AGN ~1-5%. Bursts contribute little to scatter in the SFR-stellar mass relation. In fact, many systems identified as 'ongoing' mergers will be forming stars in their 'normal' (non-burst) mode. Counting this as 'merger-induced' star formation leads to a stronger apparent redshift evolution in the contribution of mergers to the SFR density.

preprint2009arXiv

Recovering Stellar Population Properties and Redshifts from Broad-Band Photometry of Simulated Galaxies: Lessons for SED Modeling

We present a detailed analysis of our ability to determine stellar masses, ages, reddening and extinction values, and star formation rates of high-redshift galaxies by modeling broad-band SEDs with stellar population synthesis. In order to do so, we computed synthetic optical-to-NIR SEDs for model galaxies taken from hydrodynamical merger simulations placed at redshifts 1.5 < z < 3. Viewed under different angles and during different evolutionary phases, the simulations represent a wide variety of galaxy types (disks, mergers, spheroids). We show that simulated galaxies span a wide range in SEDs and color, comparable to these of observed galaxies. In all star-forming phases, dust attenuation has a large effect on colors, SEDs, and fluxes. The broad-band SEDs were then fed to a standard SED modeling procedure and resulting stellar population parameters were compared to their true values. Disk galaxies generally show a decent median correspondence between the true and estimated mass and age, but suffer from large uncertainties. During the merger itself, we find larger offsets (e.g., log M_recovered - log M_true = -0.13^{+0.10}_{-0.14}). E(B-V) values are generally recovered well, but the estimated total visual absorption Av is consistently too low, increasingly so for larger optical depths. Since the largest optical depths occur during the phases of most intense star formation, it is for the highest SFRs that we find the largest underestimates. The masses, ages, E(B-V), Av, and SFR of merger remnants (spheroids) are very well reproduced. We discuss possible biases in SED modeling results caused by mismatch between the true and template star formation history, dust distribution, metallicity variations and AGN contribution.

preprint2009arXiv

Resonant stripping as the origin of dwarf spheroidal galaxies

Dwarf spheroidal galaxies are the most dark matter dominated systems in the nearby Universe and their origin is one of the outstanding puzzles of how galaxies form. Dwarf spheroidals are poor in gas and stars, making them unusually faint, and those known as ultra-faint dwarfs have by far the lowest measured stellar content of any galaxy. Previous theories require that dwarf spheroidals orbit near giant galaxies like the Milky Way, but some dwarfs have been observed in the outskirts of the Local Group. Here we report simulations of encounters between dwarf disk galaxies and somewhat larger objects. We find that the encounters excite a process, which we term ``resonant stripping'', that can transform them into dwarf spheroidals. This effect is distinct from other mechanisms proposed to form dwarf spheroidals, including mergers, galaxy-galaxy harassment, or tidal and ram pressure stripping, because it is driven by gravitational resonances. It may account for the observed properties of dwarf spheroidals in the Local Group, including their morphologies and kinematics. Resonant stripping predicts that dwarf spheroidals should form through encounters, leaving detectable long stellar streams and tails.

preprint2009arXiv

Substructure depletion in the Milky Way halo by the disk

We employ numerical simulations and simple analytical estimates to argue that dark matter substructures orbiting in the inner regions of the Galaxy can be efficiently destroyed by disk shocking, a dynamical process known to affect globular star clusters. We carry out a set of fiducial high-resolution collisionless simulations in which we adiabatically grow a disk, allowing us to examine the impact of the disk on the substructure abundance. We also track the orbits of dark matter satellites in the high-resolution Aquarius simulations and analytically estimate the cumulative halo and disk shocking effect. Our calculations indicate that the presence of a disk with only 10% of the total Milky Way mass can significantly alter the mass function of substructures in the inner parts of halos. This has important implications especially for the relatively small number of satellites seen within ~30 kpc of the Milky Way center, where disk shocking is expected to reduce the substructure abundance by a factor of ~2 at 10^9 M$_{\odot}$ and ~3 at 10^7 M$_{\odot}$. The most massive subhalos with 10^10 M$_{\odot}$ survive even in the presence of the disk. This suggests that there is no inner missing satellite problem, and calls into question whether these substructures can produce transient features in disks, like multi-armed spiral patterns. Also, the depletion of dark matter substructures through shocking on the baryonic structures of the disk and central bulge may aggravate the problem to fully account for the observed flux anomalies in gravitational lens systems, and significantly reduces the dark matter annihilation signal expected from nearby substructures in the inner halo.

preprint2009arXiv

The Effects of Gas on Morphological Transformation in Mergers: Implications for Bulge and Disk Demographics

Transformation of disks into spheroids via mergers is a well-accepted element of galaxy formation models. However, recent simulations have shown that bulge formation is suppressed in increasingly gas-rich mergers. We investigate the global implications of these results in a cosmological framework, using independent approaches: empirical halo-occupation models (where galaxies are populated in halos according to observations) and semi-analytic models. In both, ignoring the effects of gas in mergers leads to the over-production of spheroids: low and intermediate-mass galaxies are predicted to be bulge-dominated (B/T~0.5 at <10^10 M_sun), with almost no bulgeless systems), even if they have avoided major mergers. Including the different physical behavior of gas in mergers immediately leads to a dramatic change: bulge formation is suppressed in low-mass galaxies, observed to be gas-rich (giving B/T~0.1 at <10^10 M_sun, with a number of bulgeless galaxies in good agreement with observations). Simulations and analytic models which neglect the similarity-breaking behavior of gas have difficulty reproducing the strong observed morphology-mass relation. However, the observed dependence of gas fractions on mass, combined with suppression of bulge formation in gas-rich mergers, naturally leads to the observed trends. Discrepancies between observations and models that ignore the role of gas increase with redshift; in models that treat gas properly, galaxies are predicted to be less bulge-dominated at high redshifts, in agreement with the observations. We discuss implications for the global bulge mass density and future observational tests.

preprint2009arXiv

The Formation of High Redshift Submillimeter Galaxies

We describe a model for the formation of \zsim 2 Submillimeter Galaxies (SMGs) which simultaneously accounts for both average and bright SMGs while providing a reasonable match to their mean observed spectral energy distributions (SEDs). By coupling hydrodynamic simulations of galaxy mergers with the high resolution 3D polychromatic radiative transfer code Sunrise, we find that a mass sequence of merger models which use observational constraints as physical input naturally yield objects which exhibit black hole, bulge, and H2 gas masses similar to those observed in SMGs. The dominant drivers behind the 850 micron flux are the masses of the merging galaxies and the stellar birthcloud covering fraction. The most luminous (S850 ~ 15 mJy) sources are recovered by ~10^13 Msun 1:1 major mergers with a birthcloud covering fraction close to unity, whereas more average SMGs ~5-7 mJy) may be formed in lower mass halos ~5x10^12 Msun. These models demonstrate the need for high spatial resolution hydrodynamic and radiative transfer simulations in matching both the most luminous sources as well as the full SEDs of SMGs. While these models suggest a natural formation mechanism for SMGs, they do not attempt to match cosmological statistics of galaxy populations; future efforts along this line will help ascertain the robustness of these models.

preprint2009arXiv

The formation of the first stars and galaxies

Observations made using large ground-based and space-borne telescopes have probed cosmic history all the way from the present-day to a time when the Universe was less than a tenth of its present age. Earlier on lies the remaining frontier, where the first stars, galaxies, and massive black holes formed. They fundamentally transformed the early Universe by endowing it with the first sources of light and chemical elements beyond the primordial hydrogen and helium produced in the Big Bang. The interplay of theory and upcoming observations promises to answer the key open questions in this emerging field.

preprint2009arXiv

The Star-Forming Molecular Gas in High Redshift Submillimeter Galaxies

We present a model for the CO molecular line emission from high redshift Submillimeter Galaxies (SMGs). By combining hydrodynamic simulations of gas rich galaxy mergers with the polychromatic radiative transfer code, Sunrise, and the 3D non-LTE molecular line radiative transfer code, Turtlebeach, we show that if SMGs are typically a transient phase of major mergers, their observed compact CO spatial extents, broad line widths, and high excitation conditions (CO SED) are naturally explained. In this sense, SMGs can be understood as scaled-up analogs to local ULIRGs. We utilize these models to investigate the usage of CO as an indicator of physical conditions. We find that care must be taken when applying standard techniques. The usage of CO line widths as a dynamical mass estimator from SMGs can possibly overestimate the true enclosed mass by a factor ~1.5-2. At the same time, assumptions of line ratios of unity from CO J=3-2 (and higher lying lines) to CO (J=1-0) will oftentimes lead to underestimates of the inferred gas mass. We provide tests for these models by outlining predictions for experiments which are imminently feasible with the current generation of bolometer arrays and radio-wave spectrometers.

preprint2008arXiv

A Characteristic Division Between the Fueling of Quasars and Seyferts: Five Simple Tests

Given the existence of the M_BH-sigma relation, models of self-regulated black hole (BH) growth require both a fuel supply and growth of the host bulge to deepen the potential, or else the system will either starve or self-regulate without sustained activity. This suggests that bright quasars must be triggered in major mergers: a large fraction of the galaxy must be converted to new bulge mass in a dynamical time or less. Low-luminosity AGN, in contrast, require little bulge growth and small gas supplies, and could be triggered in more common non-merger events. This predicts a transition to merger-induced fueling around the traditional quasar-Seyfert luminosity divide (growth of BH masses above/below 10^7 M_sun). We compile observations to test several predictions of such a division, including: (1) A transition to bulge-dominated hosts. (2) A transition between 'pseudobulges' and 'classical' bulges hosting the remnant BHs: pseudobulges are formed in secular processes and minor mergers, whereas classical bulges are relics of major mergers. (3) An increase in the amplitude of small-scale clustering where mergers are more efficient. (4) Different redshift evolution, with gas-rich merger rates rising to redshifts z>2 while secular processes are relatively constant in time. (5) An increasing prominence of post-starburst features in more luminous systems. Our compilation of observations provides tentative evidence for the predicted division around the Seyfert-quasar threshold. We discuss how future observations can improve these constraints and break degeneracies between different fueling models.

preprint2008arXiv

Dissipation and Extra Light in Galactic Nuclei: IV. Evolution in the Scaling Relations of Spheroids

We develop a model for the origins and redshift evolution of spheroid scaling relations. We consider spheroid sizes, velocity dispersions, masses, profile shapes (Sersic indices), and black hole (BH) masses, and their related scalings. Our approach combines advantages of observational constraints in halo occupation models and hydrodynamic merger simulations. This allows us to separate the relative roles of dissipation, dry mergers, formation time, and progenitor evolution, and identify their effects on scalings at each redshift. Dissipation is the most important factor determining spheroid sizes and fundamental plane (FP) scalings, and can account for the FP tilt and differences between disk and spheroid scalings. Because disks at high-z have higher gas fractions, mergers are more gas-rich, yielding more compact spheroids. This predicts mass-dependent evolution in spheroid sizes, in agreement with observations. This relates to subtle evolution in the FP, important to studies that assume a fixed intrinsic FP. This also predicts mild evolution in BH-host correlations, towards larger BHs at higher z. Dry mergers are significant, but only for massive systems which form early: they form compact, but undergo dry mergers (consistent with observations) such that their sizes at later times are similar to spheroids of similar mass formed more recently. We model descendants of observed compact high-z spheroids: most will become cores of BCGs, with sizes, velocity dispersions, and BH masses consistent with observations, but we identify a fraction that might survive to z=0 intact.

preprint2008arXiv

Galaxy morphology, kinematics and clustering in a hydrodynamic simulation of a LambdaCDM universe

We explore galaxy properties and their link with environment and clustering using a population of ~1000 galaxies formed in a high resolution hydrodynamic simulation of the Lambda CDM cosmology. At the redshift we concentrate on, z=1, the spatial resolution is 1.4 proper kpc/h and Milky-way sized disk galaxies contain ~10^5 particles within their virial radii. We include supermassive black hole accretion and feedback as well as a multiphase model for star formation. We find that a number of familiar qualitative relationships hold approximately between galaxy properties, for example, galaxies lie between two broad extremes of type, where ``late'' types tend to be smaller in size, have lower circular velocities, younger stars, higher star formation rates, larger disk to bulge ratios and lower Sersic indices than ``early types''. As in previous studies the stellar component of disk galaxies is not as rotationally supported as in observations. Bulges contain too much of the stellar mass, although disks do have scale lengths compatible with observations. The addition of black hole physics to the simulations does not appear to have an impact on the angular momentum results, nor do we find that it is affected in an identical simulation with significantly lower mass resolution. Despite this, we can profitably use the rank order of either disk to total ratio, Sersic index, or galaxy age to separate galaxies into morphological classes and examine the density-morphology relation and morphology dependence of clustering. We find that while at redshift z=0, the well known preponderance of early types in dense environments is seen, at z=1 the density-morphology relation becomes flatter and late type galaxies are even seen to have a higher clustering amplitude than early types (abridged).

preprint2008arXiv

How Do Disks Survive Mergers?

We develop a physical model for how galactic disks survive and/or are destroyed in interactions. Based on dynamical arguments, we show gas primarily loses angular momentum to internal torques in a merger. Gas within some characteristic radius (a function of the orbital parameters, mass ratio, and gas fraction of the merging galaxies), will quickly lose angular momentum to the stars sharing the perturbed disk, fall to the center and be consumed in a starburst. A similar analysis predicts where violent relaxation of the stellar disks is efficient. Our model allows us to predict the stellar and gas content that will survive to re-form a disk in the remnant, versus being violently relaxed or contributing to a starburst. We test this in hydrodynamic simulations and find good agreement as a function of mass ratio, orbital parameters, and gas fraction, in simulations spanning a wide range in these properties and others, including different prescriptions for gas physics and feedback. In an immediate sense, the amount of disk that re-forms can be understood in terms of well-understood gravitational physics, independent of details of ISM gas physics or feedback. This allows us to explicitly quantify the requirements for such feedback to (indirectly) enable disk survival, by changing the pre-merger gas content and distribution. The efficiency of disk destruction is a strong function of gas content: we show how and why sufficiently gas-rich major mergers can, under general conditions, yield systems with small bulges (B/T<0.2). We provide prescriptions for inclusion of our results in semi-analytic models.

preprint2008arXiv

Incidence Rate of GRB-host-DLAs at High Redshift

We study the incidence rate of damped Ly-a systems associated with the host galaxies of gamma-ray bursts (GRB-host-DLAs) as functions of neutral hydrogen column density (N_HI) and projected star formation rate (SFR) using cosmological SPH simulations. Assuming that the occurrence of GRBs is correlated with the local SFR, we find that the median N_HI of GRB-host-DLAs progressively shifts to lower N_HI values with increasing redshift, and the incidence rate of GRB-host-DLAs with log N_HI > 21.0 decreases rapidly at z>=6. Our results suggest that the likelihood of observing the signature of IGM attenuation in GRB afterglows increases towards higher redshift, because it will not be blocked by the red damping wing of DLAs in the GRB host galaxies. This enhances the prospects of using high-redshift GRBs to probe the reionization history of the Universe. The overall incidence rate of GRB-host-DLAs decreases monotonically with increasing redshift, whereas that of QSO-DLAs increases up to z=6. A measurement of the difference between the two incidence rates would enable an estimation of the value of η_grb, which is the mass fraction of stars that become GRBs for a given amount of star formation. Our predictions can be tested by upcoming high-z GRB missions, including JANUS (Joint Astrophysics Nascent Universe Scout) and SVOM (Space multi-band Variable Object Monitor).

preprint2007arXiv

Observational Evidence for the Co-evolution of Galaxy Mergers, Quasars, and the Blue/Red Galaxy Transition

We compile a number of observations to estimate the time-averaged rate of formation or buildup of red sequence galaxies, as a function of mass and redshift. Comparing this with the mass functions of mergers and quasar hosts, and independently comparing their clustering properties as a function of redshift, we find that these populations trace the same mass distribution, with similar evolution, at redshifts 0<z<~1.5. Knowing one of the quasar, merger, or elliptical mass/luminosity functions, it is possible to predict the others. Allowing for greater model dependence, we compare the rate of early-type buildup with the implied merger and quasar triggering rates as a function of mass and redshift and find agreement. Over this redshift range, observed merger fractions can account for the entire bright quasar luminosity function and buildup of the red sequence at all but the highest masses at low redshift (>~10^11 M_solar at z<~0.3) where 'dry' mergers appear to dominate. This supports a necessary prediction of theories where mergers between gas-rich galaxies produce ellipticals with an associated phase of quasar activity, after which the remnant becomes red. These populations trace a similar characteristic transition mass, possibly reflecting the mass above which the elliptical population is mostly (>~50%) assembled at a given redshift, which increases with redshift over the observed range in a manner consistent with suggestions that cosmic downsizing may apply to red galaxy assembly as well as star formation. These mass distributions as a function of redshift do not uniformly trace the all/red/blue galaxy population, ruling out models in which quasar activity is generically associated with star formation or is long lived in 'old' systems.

preprint2007arXiv

Photometric Properties of the Most Massive High-Redshift Galaxies

We calculate the observable properties of the most massive high-redshift galaxies in the hierarchical formation scenario where stellar spheroid and supermassive black hole growth are fueled by gas-rich mergers. Combining high-resolution hydrodynamical simulations of the hierarchical formation of a z~6 quasar, stellar population synthesis models, template AGN spectra, prescriptions for interstellar and intergalactic absorption, and the response of modern telescopes, the photometric evolution of galaxies destined to host z~6 quasars are modeled at redshifts z~4-14. These massive galaxies, with enormous stellar masses of M_star ~10^11.5-10^12 M_sun. and star formation rates of SFR~10^3-10^4 M_sun yr^-1 at z>~7, satisfy a variety of photometric selection criteria based on Lyman-break techniques including V-band dropouts at z>~5, i-band dropouts at z>~6, and z-band dropouts at z>~7. The observability of the most massive high-redshift galaxies is assessed and compared with a wide range of existing and future photometric surveys including SDSS, GOODS/HUDF, NOAO WDFS, UKIDSS, the IRAC Shallow Survey, Pan-STARRS, LSST, and SNAP. Massive stellar spheroids descended from z~6 quasars will likely be detected at z~4 by existing surveys, but owing to their low number densities the discovery of quasar progenitor galaxies at z>7 will likely require future surveys of large portions of the sky (>~0.5%) at wavelengths lambda>1 micron. The detection of rare, star-bursting, massive galaxies at redshifts z>~6 would provide support for the hierarchical formation of the earliest quasars and characterize the primitive star-formation histories of the most luminous elliptical galaxies.

preprint2006arXiv

An Upper Limit to the Degree of Evolution Between Supermassive Black Holes and their Host Galaxies

We describe a model-independent integral constraint which defines an upper limit to the allowed degree of evolution in the ratio of black hole (BH) mass to host galaxy luminosity or mass, as a function of redshift. Essentially, if the BH/host ratio is excessive at redshift z, then it would imply that the total mass density in BHs above some M_BH(min) is larger at that redshift than at z=0, which is impossible. This argument requires no knowledge of host or BH properties, only a lower limit to the observed luminosity density in the brightest galaxies at some z. We calculate this upper limit from a variety of luminosity and mass functions in different bands from redshifts z=0-2. We show that it is consistent with passive evolution of spheroid populations (with a fixed M_BH/M_host relation) in all cases, and provides tighter constraints than have generally been obtained previously, ruling out at >6 sigma observational and theoretical estimates suggesting that M_BH/M_host was significantly larger at high redshifts than locally, although relatively weak (factor ~2 by z=2) evolution is still allowed. We discuss a variety of possible 'loopholes' or changes in the BH/host populations and correlations, and show that they typically lower the upper limits and strengthen our conclusions.

preprint2006arXiv

Formation of Primordial Stars in a LCDM Universe

We study the formation of the first generation of stars in the standard cold dark matter model, using a very high-resolution hydordynamic simulations. Our simulation achieves a dynamic range of 10^{10} in length scale. With accurate treatment of atomic and molecular physics, it allows us to study the chemo-thermal evolution of primordial gas clouds to densities up to n = 10^{16}/cc without assuming any a priori equation of state; a six orders of magnitudes improvement over previous three-dimensional calculations. All the relevant atomic and molecular cooling and heating processes, including cooling by collision-induced continuum emission, are implemented. For calculating optically thick H2 cooling at high densities, we use the Sobolev method. To examine possible gas fragmentation owing to thermal instability, we compute explicitly the growth rate of isobaric perturbations. We show that the cloud core does not fragment in either the low-density or high-density regimes. We also show that the core remains stable against gravitational deformation and fragmentation. We obtain an accurate gas mass accretion rate within a 10 Msun innermost region around the protostar. The protostar is accreting the surrounding hot gas at a rate of 0.001-0.01 Msun/yr. From these findings we conclude that primordial stars formed in early minihalos are massive. We carry out proto-stellar evolution calculations using the obtained accretion rate. The resulting mass of the first star is M_ZAMS = 60-100 Msun, with the exact mass dependent on the actual accretion rate.

preprint2006arXiv

The Relation Between Quasar and Merging Galaxy Luminosity Functions and the Merger-Induced Star Formation Rate of the Universe

Using a model for self-regulated growth of black holes (BHs) in mergers involving gas-rich galaxies, we study the relationship between quasars and the population of merging galaxies and predict the merger-induced star formation rate density of the Universe. Mergers drive nuclear gas inflows, fueling starbursts and 'buried quasars' until accretion feedback expels the gas, rendering a briefly visible optical quasar. Star formation is shut down and accretion declines, leaving a passively evolving remnant with properties typical of red, elliptical galaxies. Based on evolution of these events in our simulations, we demonstrate that the observed statistics of merger rates, luminosity functions (LFs) and mass functions, SFR distributions, specific SFRs, quasar and quasar host galaxy LFs, and elliptical/red galaxy LFs are self-consistent and follow from one another as predicted by the merger hypothesis. We use our simulations to de-convolve both quasar and merging galaxy LFs to determine the birthrate of black holes of a given final mass and merger rates as a function of stellar mass. We use this to predict the merging galaxy LF in several observed wavebands, color-magnitude relations, mass functions, absolute and specific SFR distributions and SFR density, and quasar host galaxy LFs, as a function of redshift from z=0-6. We invert this and predict e.g. quasar LFs from observed merger LFs or SFR distributions. Our results agree well with observations, but idealized models of quasar lightcurves are ruled out by comparison of merger and quasar observations at >99.9% confidence. Using only observations of quasars, we estimate the contribution of mergers to the SFR density of the Universe even to high redshifts z~4.

preprint2005arXiv

The Evolution of the M-sigma Relation

(Abridged) We examine the evolution of the black hole mass - stellar velocity dispersion (M-sigma) relation over cosmic time using simulations of galaxy mergers that include feedback from supermassive black hole growth. We consider mergers of galaxies varying the properties of the progenitors to match those expected at redshifts z=0-6. We find that the slope of the resulting M-sigma relation is the same at all redshifts considered. For the same feedback efficiency that reproduces the observed amplitude of the M-sigma relation at z=0, there is a weak redshift-dependence to the normalization that results from an increasing velocity dispersion for a given galactic stellar mass. We develop a formalism to connect redshift evolution in the M-sigma relation to the scatter in the local relation at z=0. We show that the scatter in the local relation places severe constraints on the redshift evolution of both the normalization and slope of the M-sigma relation. Furthermore, we demonstrate that cosmic downsizing introduces a black hole mass-dependent dispersion in the M-sigma relation and that the skewness of the distribution about the locally observed M-sigma relation is sensitive to redshift evolution in the normalization and slope. In principle, these various diagnostics provide a method for differentiating between theories for producing the M-sigma relation. In agreement with existing constraints, our simulations imply that hierarchical structure formation should produce the relation with small intrinsic scatter.

preprint2004arXiv

A QED Model for the Origin of Bursts from SGRs and AXPs

We propose a model to account for the bursts from soft gamma repeaters (SGRs) and anomalous X-ray pulsars (AXPs) in which quantum electrodynamics (QED) plays a vital role. In our theory, that we term ``fast-mode breakdown,'' magnetohydrodynamic (MHD) waves that are generated near the surface of a neutron star and propagate outward through the magnetosphere will be modified by the polarization of the vacuum. For neutron star magnetic fields $B_\rmscr{NS} \gtrsim B_\rmscr{QED} \approx 4.4 \times 10^{13}$ G, the interaction of the wave fields with the vacuum produces non-linearities in fast MHD waves that can steepen in a manner akin to the growth of hydrodynamic shocks. Under certain conditions, fast modes can develop field discontinuities on scales comparable to an electron Compton wavelength, at which point the wave energy will be dissipated through electron-positron pair production. We show that this process operates if the magnetic field of the neutron star is sufficiently strong and the ratio of the wavelength of the fast mode to its amplitude is sufficiently small, in which case the wave energy will be efficiently converted into an extended pair-plasma fireball. The radiative output from this fireball will consist of hard X-rays and soft $γ$-rays, with a spectrum similar to those seen in bursts from SGRs and AXPs (abridged).

preprint2004arXiv

Detecting Sunyaev-Zel'dovich clusters with PLANCK: I. Construction of all-sky thermal and kinetic SZ-maps

All-sky thermal and kinetic Sunyaev-Zel'dovich (SZ) maps are presented for assessing how well the PLANCK-mission can find and characterise clusters of galaxies, especially in the presence of primary anisotropies of the cosmic microwave background (CMB) and various galactic and ecliptic foregrounds. The maps have been constructed from numerical simulations of structure formation in a standard LCDM cosmology and contain all clusters out to redshifts of z = 1.46 with masses exceeding 5e13 M_solar/h. By construction, the maps properly account for the evolution of cosmic structure, the halo-halo correlation function, the evolving mass function, halo substructure and adiabatic gas physics. The velocities in the kinetic map correspond to the actual density environment at the cluster positions. We characterise the SZ-cluster sample by measuring the distribution of angular sizes, the integrated thermal and kinetic Comptonisations, the source counts in the three relevant PLANCK-channels, and give the angular power-spectra of the SZ-sky. While our results are broadly consistent with simple estimates based on scaling relations and spherically symmetric cluster models, some significant differences are seen which may affect the number of cluster detectable by PLANCK.

preprint2004arXiv

Massive galaxies at redshift 2 in cosmological hydrodynamic simulations

We study the properties of galaxies at z=2 in a Lambda cold dark matter universe, using two different types of hydrodynamic simulation methods -- Eulerian TVD and smoothed particle hydrodynamics (SPH) -- and a spectrophotometric analysis in the U_n, G, R filter set. The simulated galaxies at z=2 satisfy the color-selection criteria proposed by Adelberger et al. (2004) and Steidel et al. (2004) when we assume Calzetti extinction with E(B-V)=0.15. We find that the number density of simulated galaxies brighter than R<25.5 at z=2 is about 1e-2 h^3 Mpc^-3 for E(B-V)=0.15, which is roughly twice that of the number density found by Erb et al. (2004) for the UV bright sample. This suggests that roughly half of the massive galaxies with M*>10^{10} Msun/h at z=2 are UV bright population, and the other half is bright in the infra-red wavelengths. The most massive galaxies at z=2 have stellar masses >= 10^{11-12} Msun. They typically have been continuously forming stars with a rate exceeding 30 Msun/yr over a few Gyrs from z=10 to z=2, together with significant contribution by starbursts reaching up to 1000 Msun/yr which lie on top of the continuous component. TVD simulations indicate a more sporadic star formation history than the SPH simulations. Our results do not imply that hierarchical galaxy formation fails to account for the observed massive galaxies at z>=1. The global star formation rate density in our simulations peaks at z>=5, a much higher redshift than predicted by the semianalytic models. This star formation history suggests early build-up of the stellar mass density, and predicts that 70 (50, 30)% of the total stellar mass at z=0 had already been formed by z=1 (2, 3). Upcoming observations by Spitzer and Swift might help to better constrain the star formation history at high redshift.

preprint2004arXiv

Massive galaxies in cosmological simulations: UV-selected sample at redshift z=2

We study the properties of galaxies at z=2 in a Lambda CDM universe, using two different types of hydrodynamic simulation methods (Eulerian TVD and SPH) and a spectrophotometric analysis in the Un, G, R filter set. The simulated galaxies at z=2 satisfy the color-selection criteria proposed by Adelberger et al. (2004) when we assume Calzetti extinction with E(B-V)=0.15. We find that the number density of simulated galaxies brighter than R<25.5 at z=2 is about 2e-2 h^3/Mpc^3, roughly one order of magnitude larger than that of Lyman break galaxies at z=3. The most massive galaxies at z=2 have stellar masses >~1e11 Msun, and their observed-frame G-R colors lie in the range 0.0<G-R<1.0. They typically have been continuously forming stars with a rate exceeding 30 Msun/yr over a few Gyrs from z=10 to z=2, although the TVD simulation indicates a more sporadic star formation history than the SPH simulations. Of order half of their stellar mass was already assembled by z~4. The reddest massive galaxies at z=2 with G-R >= 1.0 and Mstar>1e10 Msun/h finished the build-up of their stellar mass by z~3. Interestingly, our study suggests that the majority of the most massive galaxies at z=2 should be detectable at rest-frame UV wavelengths, contrary to some recent claims made on the basis of near-IR studies of galaxies at the same epoch, provided the median extinction is less than E(B-V)<0.3. However, our results also suggest that the fraction of stellar mass contained in galaxies that pass the color-selection criteria could be as low as 50% of the total stellar mass in the Universe at z=2. Our simulations suggest that the missing stellar mass is contained in fainter (R>25.5) and intrinsically redder galaxies. Our results do not suggest that hierarchical galaxy formation fails to account for the massive galaxies at z>=1. (abridged)

preprint2003arXiv

Early Structure Formation and Reionization in a Warm Dark Matter Cosmology

We study first structure formation in Lambda-dominated universes using large cosmological N-body/SPH simulations. We consider a standard LCDM model and a LWDM model in which the mass of the dark matter particles is taken to be m_X=10 keV. The linear power spectrum for the LWDM model has a characteristic cut-off at a wavenumber k=200 /Mpc, suppressing the formation of low mass (< 10^6 Msun) nonlinear objects early on. The absence of low mass halos in the WDM model makes the formation of primordial gas clouds with molecular hydrogen very inefficient at high redshifts. The first star-forming gas clouds form at z~21 in the WDM model, considerably later than in the CDM counterpart, and the abundance of these gas clouds differs by an order of magnitude between the two models. We carry out radiative transfer calculations by embedding massive Population III stars in the gas clouds. We show that the volume fraction of ionized gas rises up close to 100% by z=18 in the CDM case, whereas that of the WDM model remains extremely small at a level of a few percent. Thus the WDM model with m_X=10 keV is strongly inconsistent with the observed high optical depth by the WMAP satellite.

preprint2002arXiv

Cosmological SPH simulations: A hybrid multi-phase model for star formation

We present a model for star formation and supernova feedback that describes the multi-phase structure of star forming gas on scales that are typically not resolved in cosmological simulations. Our approach includes radiative heating and cooling, the growth of cold clouds embedded in an ambient hot medium, star formation in these clouds, feedback from supernovae in the form of thermal heating and cloud evaporation, galactic winds and outflows, and metal enrichment. Implemented using SPH, our scheme is a significantly modified and extended version of the grid-based method of Yepes et al. (1997), and enables us to achieve high dynamic range in simulations of structure formation. We discuss properties of the feedback model in detail and show that it predicts a self-regulated, quiescent mode of star formation, which, in particular, stabilises the star forming gaseous layers of disk galaxies. The parameterisation of this mode can be reduced to a single free quantity which determines the overall timescale for star formation. We fix this parameter to match the observed rates of star formation in local disk galaxies. When normalised in this manner, cosmological simulations nevertheless overproduce the observed cosmic abundance of stellar material. We are thus motivated to extend our feedback model to include galactic winds associated with star formation. Using small-scale simulations of individual star-forming disk galaxies, we show that these winds produce either galactic fountains or outflows, depending on the depth of the gravitational potential. Moreover, outflows from galaxies in these simulations drive chemical enrichment of the intergalactic medium, in principle accounting for the presence of metals in the Lyman alpha forest. (abridged)

preprint2002arXiv

High redshift galaxies and the Lyman-alpha forest in a CDM universe

We use a cosmological hydrodynamic simulation of a cold dark matter universe to investigate theoretically the relationship between high redshift galaxies and the Lyman=alpha forest at redshift z=3. Galaxies in the simulation are surrounded by halos of hot gas, which nevertheless contain enough neutral hydrogen to cause a Ly-alpha flux decrement, its strength increasing with galaxy mass. A comparison with recent observational data by Adelberger et. al on the Ly-alpha forest around galaxies reveals that actual galaxies may have systematically less Ly-alpha absorption within 1 Mpc of them than our simulated galaxies. In order to investigate this possibility, we add several simple prescriptions for galaxy feedback on the IGM to the evolved simulation. These include the effect of photoionizing background radiation coming from galactic sources, galactic winds whose only effect is to deposit thermal energy into the IGM, and another, kinetic model for winds, which are assumed to evacuate cavities in the IGM around galaxies. We find that only the latter is able to produce a large effect, enough to match the tentative observational data, given the energy available from star formation in the simulated galaxies. Another intriguing possibility is that a selection effect is responsible, so that galaxies with low Ly-alpha absorption are preferentially included in the sample. This is also viable, but predicts very different galaxy properties (including clustering) than the other scenarios.

preprint2002arXiv

Towards a Precise Measurement of Matter Clustering: Lyman-alpha Forest Data at Redshifts 2-4

We measure the filling factor, correlation function, and power spectrum of transmitted flux in a large sample of Lya forest spectra, comprised of 30 Keck HIRES spectra and 23 Keck LRIS spectra. We infer the linear matter power spectrum P(k) from the flux power spectrum P_F(k), using an improved version of the method of Croft et al. (1998) that accounts for the influence of z-space distortions, non- linearity, and thermal broadening on P_F(k). The evolution of the shape and amplitude of P(k) over the range z= 2-4 is consistent with gravitational instability, implying that non-gravitational fluctuations do not make a large contribution. Our fiducial measurement of P(k) comes from data with <z> = 2.72. It has amplitude Delta^2(k_p)=0.74^0.20_-0.16 at wavenumber k_p=0.03 (km/s)^-1 and is well described by a power-law of index -2.43 +/- 0.06 or by a CDM-like power spectrum with shape parameter Gamma'=1.3^+0.7_-0.5*10^-3 (km/s) at z=2.72. For Omega_m=0.4, Omega_Lam=0.6, the best-fit Gamma =0.16 (h^-1mpc)^-1, in good agreement with the 2dF Galaxy Redshift Survey, and the best-fit sigma_8=0.82 (Gamma/0.15)^-0.44. Matching the observed cluster mass function and our Delta^2(k_p) in spatially flat models requires Omega_m=0.38^+0.10_-0.08 + 2.2 (Gamma-0.15). Matching Delta^2(k_p) in COBE-normalized, flat CDM models with no tensor fluctuations requires Omega_m = (0.29 +/-0.04) n^-2.89 h_65^-1.9. The Lya forest complements other probes of P(k) by constraining a regime of redshift and lengthscale not accessible by other means, and the consistency of these inferred parameters with independent estimates provides further support for inflation, cold dark matter, and vacuum energy (abridged).

preprint2001arXiv

Simulations of Damped Lyman-Alpha and Lyman Limit Absorbers in Different Cosmologies: Implications for Structure Formation at High Redshift

We use hydrodynamic cosmological simulations to study damped Lyman-alpha (DLA) and Lyman limit (LL) absorption at redshifts z=2-4 in five variants of the cold dark matter scenario. Our standard simulations resolve the formation of dense concentrations of neutral gas in halos with circular velocity v_c roughly 140 km/s for Omega_m=1 and 90 km/s for Omega_m=0.4, at z=2; an additional LCDM simulation resolves halos down to v_c approximately 50 km/s at z=3. We find a clear relation between HI column density and projected distance to the center of the nearest galaxy, with DLA absorption usually confined to galactocentric radii less than 10-15 kpc and LL absorption arising out to projected separations of 30 kpc or more. Detailed examination provides evidence of non-equilibrium effects on absorption cross-section. If we consider only absorption in the halos resolved by our standard simulations, then all five models fall short of reproducing the observed abundance of DLA and LL systems at these redshifts. If we extrapolate to lower halo masses, we find all four models are consistent with the observed abundance of DLA systems if the the extrapolated behavior extends to circular velocities roughly 50-80 km/s, and they may produce too much absorption if the relation continues to 40 km/s. Our results suggest that LL absorption is closely akin to DLA absorption, arising in less massive halos or at larger galactocentric radii but not caused by processes acting on a radically different mass scale.

preprint2000arXiv

Baryons in the Warm-Hot Intergalactic Medium

Approximately 30-40% of all baryons in the present day universe reside in a warm-hot intergalactic medium (WHIM), with temperatures between 10^5<T<10^7 K. This is a generic prediction from six hydrodynamic simulations of currently favored structure formation models having a wide variety of numerical methods, input physics, volumes, and spatial resolutions. Most of these warm-hot baryons reside in diffuse large-scale structures with a median overdensity around 10-30, not in virialized objects such as galaxy groups or galactic halos. The evolution of the WHIM is primarily driven by shock heating from gravitational perturbations breaking on mildly nonlinear, non-equilibrium structures such as filaments. Supernova feedback energy and radiative cooling play lesser roles in its evolution. WHIM gas is consistent with observations of the 0.25 keV X-ray background without being significantly heated by non-gravitational processes because the emitting gas is very diffuse. Our results confirm and extend previous work by Cen & Ostriker and Dave' et al.

preprint2000arXiv

Cooling Radiation and the Lyman-alpha Luminosity of Forming Galaxies

We examine the cooling radiation from forming galaxies in hydrodynamic simulations of the LCDM model (cold dark matter with a cosmological constant), focusing on the Ly-alpha line luminosities of high-redshift systems. Primordial composition gas condenses within dark matter potential wells, forming objects with masses and sizes comparable to the luminous regions of observed galaxies. As expected, the energy radiated in this process is comparable to the gravitational binding energy of the baryons, and the total cooling luminosity of the galaxy population peaks at z ~= 2. However, in contrast to the classical picture of gas cooling from the \sim 10^6 K virial temperature of a typical dark matter halo, we find that most of the cooling radiation is emitted by gas with T < 20,000 K. As a consequence, roughly 50% of this cooling radiation emerges in the Ly-alpha line. While a galaxy's cooling luminosity is usually smaller than the ionizing continuum luminosity of its young stars, the two are comparable in the most massive systems, and the cooling radiation is produced at larger radii, where the Ly-alpha photons are less likely to be extinguished by dust. We suggest, in particular, that cooling radiation could explain the two large (\sim 100 kpc), luminous (L_{Ly-alpha} \sim 10^{44} erg s^{-1}) ``blobs'' of Ly-alpha emission found in Steidel et al.'s (1999) narrow band survey of a z = 3 proto-cluster. Our simulations predict objects of the observed luminosity at about the right space density, and radiative transfer effects can account for the observed sizes and line widths. We discuss observable tests of this hypothesis for the nature of the Ly-alpha blobs, and we present predictions for the contribution of cooling radiation to the Ly-alpha luminosity function of galaxies as a function of redshift.

preprint2000arXiv

Simulating the effects of intergalactic grey dust

Using a high-resolution cosmological hydrodynamic simulation, we present a method to constrain extinction due to intergalactic grey dust based on the observed magnitudes of distant Type IA supernovae. We apply several simple prescriptions to relate the intergalactic dust density to the gas density in the simulation, thereby obtaining dust extinctions that may be directly compared to the observed distribution of supernova magnitudes. Our analysis is sensitive to the spatial distribution of grey dust, but is not dependent on its intrinsic properties such as its opacity or grain size. We present an application of our technique to the supernova data of Perlmutter et al., who find that their high redshift sample is ~0.2 magnitudes fainter than the expectation for a non-accelerating, low-density universe. We find that for grey dust to be responsible, it must be distributed quite smoothly, e.g., tracing intergalactic gas. More realistic dust distributions, such as dust tracing the metal density, are inconsistent with observations at the 1.5-2 sigma level. Upcoming observations and improved modelling of the dust distribution should lead to stronger constraints on intergalactic grey dust extinction.

preprint1997arXiv

Metal Lines in Cosmological Models of Lyman-Alpha Absorbers

The metal absorption lines found in association with \lya absorbers of moderate to high HI column density contain valuable information about the metallicity and ionization conditions within the absorbers and offer a stronger test of models of the intergalactic medium at \sim 3$ than HI absorption lines alone. We have developed a method to predict the strengths of metal absorption lines within the framework of cosmological models for the \lya forest. The method consists of evaluating a quantity, the Line Observability Index, for a database of hundreds of candidate metal lines, allowing a comprehensive identification of the lines the model predicts to be detectable associated with a \lya absorber of a given HI column density and metallicity. Applying this technique to a particular class of models at \sim 2-4$, we predict that the OVI(1032 Å, 1038 Å) doublet is the only practical probe of the metallicity of low column density absorbers (HI \simlt 10^{14.5} cm^{-2}$), that CIV (1548 Å) is the strongest line with rest wavelength $λ_r > 1216$ Å regardless of HI$, and that the strongest metal lines should be CIII(977 Å) and SiIII(1206.5 Å), which peak at HI \sim 10^{17} cm^{-2}$.

preprint1997arXiv

Parallel TreeSPH

We describe PTreeSPH, a gravity treecode combined with an SPH hydrodynamics code designed for massively parallel supercomputers having distributed memory. Our computational algorithm is based on the popular TreeSPH code of Hernquist & Katz (1989). PTreeSPH utilizes a domain decomposition procedure and a synchronous hypercube communication paradigm to build self-contained subvolumes of the simulation on each processor at every timestep. Computations then proceed in a manner analogous to a serial code. We use the Message Passing Interface (MPI) communications package, making our code easily portable to a variety of parallel systems. PTreeSPH uses individual smoothing lengths and timesteps, with a communication algorithm designed to minimize exchange of information while still providing all information required to accurately perform SPH computations. We have additionally incorporated cosmology, periodic boundary conditions with forces calculated using a quadrupole Ewald summation method, and radiative cooling and heating from a parameterized ionizing background following Katz, Weinberg & Hernquist (1996). The addition of other physical processes, such as star formation, is straightforward. A cosmological simulation from z=49 to z=2 with 64^3 gas particles and 64^3 dark matter particles requires ~6000 node-hours on a Cray T3D, with a communications overhead of ~10% and is load balanced to a ~90% level. When used on the new Cray T3E, this code will be capable of performing cosmological hydrodynamical simulations down to z=0 with ~2x10^6 particles, or to z=2 with ~10^7 particles, in a reasonable amount of time. Even larger simulations will be practical in situations where the matter is not highly clustered or when periodic boundaries are not required.

preprint1996arXiv

Voigt-Profile Analysis of the Lyman-alpha Forest in a Cold Dark Matter Universe

We use an automated Voigt-profile fitting procedure to extract statistical properties of the Ly$α$ forest in a numerical simulation of an $Ω=1$, cold dark matter (CDM) universe. Our analysis method is similar to that used in most observational studies of the forest, and we compare the simulations to recently published results derived from Keck HIRES spectra. With the Voigt-profile decomposition analysis, the simulation reproduces the large number of weak lines ($N_{\rm HI}\la 10^{13}\cdunits$) found in the HIRES spectra. The column density distribution evolves significantly between $z=3$ and $z=2$, with the number of lines at fixed column density dropping by a factor $\sim 1.6$ in the range where line blending is not severe. At $z=3$, the $b$-parameter distribution has a median of $35 \kms$ and a dispersion of $20 \kms$, in reasonable agreement with the observed values. The comparison between our new analysis and recent data strengthens earlier claims that the \lya forest arises naturally in hierarchical structure formation as photoionized gas falls into dark matter potential wells. However, there are two statistically signficant discrepancies between the simulated forest and the HIRES results: the model produces too many lines at $z=3$ by a factor $\sim 1.5-2$, and it produces more narrow lines ($b<20 \kms$) than are seen in the data. The first result is sensitive to our adopted normalization of the mean \lya optical depth, and the second is sensitive to our assumption that helium reionization has not significantly raised gas temperatures at $z=3$. It is therefore too early to say whether these discrepancies indicate a fundamental problem with the high-redshift structure of the $Ω=1$ CDM model or reflect errors of detail in our modeling of the gas distribution or the observational procedure.

preprint1995arXiv

A Parallel Self--Consistent Field Code

We describe a version of an algorithm for evolving self-gravitating collections of particles that should be nearly ideal for parallel architectures. Our method is derived from the ``self-consistent field'' (SCF) approach suggested previously by Clutton-Brock and others. Owing to the use of a global description of the gravitational field, the particles in an SCF simulation do not interact with one another directly, minimizing communications overhead between nodes in a parallel implementation. Ideal load-balancing is achieved since precisely the same number of operations are needed to compute the acceleration for each particle. Consequently, the SCF technique is perfectly scalable and the size of feasible applications will grow in simple proportion to advances in computational hardware. We describe an SCF code developed for and tested on a Connection Machine 5. Empirical tests demonstrate the efficient and scalable nature of the algorithm. Depending on the application, simulations with particle numbers in the range $N\sim 10^7 - 10^{8.5}$ are now possible. Larger platforms should make simulations with billions of particles feasible in the near future. Specific astrophysical applications are discussed in the context of collisionless dynamics.

preprint1995arXiv

Cosmological Simulations with TreeSPH

We describe numerical methods for incorporating gas dynamics into cosmological simulations and present illustrative applications to the cold dark matter (CDM) scenario. Our evolution code, a version of TreeSPH (Hernquist \& Katz 1989) generalized to handle comoving coordinates and periodic boundary conditions, combines smoothed--particle hydrodynamics (SPH) with the hierarchical tree method for computing gravitational forces. The Lagrangian hydrodynamics approach and individual time steps for gas particles give the algorithm a large dynamic range, which is essential for studies of galaxy formation in a cosmological context. The code incorporates radiative cooling for an optically thin, primordial composition gas in ionization equilibrium with a user-specified ultraviolet background. We adopt a phenomenological prescription for star formation that gradually turns cold, dense, Jeans-unstable gas into collisionless stars, returning supernova feedback energy to the surrounding medium. In CDM simulations, some of the baryons that fall into dark matter potential wells dissipate their acquired thermal energy and condense into clumps with roughly galactic masses. The resulting galaxy population is insensitive to assumptions about star formation; we obtain similar baryonic mass functions and galaxy correlation functions from simulations with star formation and from simulations without star formation in which we identify galaxies directly from the cold, dense gas.

preprint1995arXiv

Excitation of Activity in Galaxies by Minor Mergers

Mergers between gas--rich disks and less--massive dwarf galaxies are studied using numerical simulation. As the orbit of a satellite decays through dynamical friction, the primary disk develops large-amplitude spirals in response to its tidal forcing. While these features arise in both the stars and the gas in the disk, the non--axisymmetric structures in the gas differ slightly from those in the stars. In particular, as a consequence of the formation of strong shocks in the gas and the effects of radiative cooling, the gas response tends to lead the stellar response, enabling the stars to strongly torque the gas. These torques deprive the gas of its angular momentum, forcing a significant fraction of it into the inner regions of the disk. The radial inflows induced by these mergers accumulate large quantities of interstellar gas in the nuclear regions of the host disks. In some cases, nearly half of all the gas initially distributed throughout the disk winds up in a dense ``cloud'' several hundred parsecs in extent. The models reported here do not include star formation and, so, we cannot determine the ultimate fate of the gas. Nevertheless, given the high densities in the nuclear gas, it is plausible to identify these concentrations of dense gas in the remnants with those accompanying intense starbursts in some active galaxies. Therefore, the calculations here provide a framework for interpreting the origin of nuclear activity in otherwise quiescent disk galaxies. To the extent that galaxy formation is a chaotic process in which large structures are built up by the accretion of smaller fragments, our models may also be relevant to starbursts and the onset of nuclear activity in proto--galaxies at high redshifts.

preprint1995arXiv

Small Scale Structure and High Redshift HI

Cosmological simulations with gas dynamics suggest that the Lyman-alpha forest is produced mainly by "small scale structure" --- filaments and sheets that are the high redshift analog of today's galaxy superclusters. There is no sharp distinction between Lyman-alpha clouds and "Gunn-Peterson" absorption produced by the fluctuating IGM -- the Lyman-alpha forest {\it is} the Gunn-Peterson effect. Lyman limit and damped Lyman-alpha absorption arises in the radiatively cooled gas of forming galaxies. At $z~2-3$, most of the gas is in the photoionized, diffuse medium associated with the Lyman-alpha forest, but most of the {\it neutral} gas is in damped Lyman-alpha systems. We discuss generic evolution of cosmic gas in a hierarchical scenario of structure formation, with particular attention to the prospects for detecting 21cm emission from high redshift HI. A scaling argument based on the present-day cluster mass function suggests that objects with $M_{HI} >~ 5e11 h^{-1} \msun$ should be extremely rare at $z~3$, so detections with existing instruments will be difficult. An instrument like the proposed Square Kilometer Array could detect individual damped Lyman-alpha systems at high redshift, making it possible to map structure in the high redshift universe in much the same way that today's galaxy redshift surveys map the local large scale structure.

preprint1995arXiv

The Disruption of the Sagittarius Dwarf Galaxy.

Numerical simulations of dwarf spheroidal galaxies undergoing several close encounters with the Milky Way are described. By comparing our models to observed properties of the recently discovered dwarf galaxy in Sagittarius (Sgr), we discuss implications of our results for the formation and evolution of the Milky Way system. We find that existing observations are not sufficient to allow us to place precise limits on either the orbit or the initial state of the dwarf. Debris from the ongoing tidal stripping of the Sagittarius galaxy are expected to form moving groups in the halo of the Galaxy and the discovery of such stars would strongly constrain the history and dynamical state of the dwarf. Furthermore, if Sgr is presently being disrupted, we predict that its remains will be detectable as a moving group in the halo for more than 1 Gyr. Thus, if similar accretion events have occurred in the recent history of the Galaxy, their aftereffects may still be observable.

preprint1994arXiv

Models of Galaxies with Central Black Holes: Adiabatic Growth in Spherical Galaxies

We add black holes to nonrotating, spherical galaxy models, with the assumption that the black-hole growth is slow compared with the dynamical time but fast compared with the relaxation time. The outcome differs depending on whether the core of the initial galaxy does or does not resemble that of an isothermal sphere. For the isothermal case the previously-known results are confirmed and sharpened: the black hole induces cusps in the density ($ρ\sim r^{-3/2}$) and velocity dispersion ($v^2\sim r^{-1}$), and a tangential anisotropy in the velocity distribution away from the center. For the non-isothermal case the induced density cusp is steeper, and the induced anisotropy is larger and penetrates right to the center. The cusp around the black hole is insensitive to anisotropy in the initial velocity distribution, and also to the origin of the black hole, unless its mass comes exclusively from the stars of lowest angular momentum, in which case the cusp is suppressed. We discuss the implications for the interpretation of evidence for massive black holes in galactic nuclei.