BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Jeremiah P. Ostriker

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

ResearcherUnclaimed source record
astro-ph.GAastro-ph.COastro-phastro-ph.HEastro-ph.IMastro-ph.SRgr-qc

Catalog footprint

What is connected

56works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

56 published item(s)

preprint2022arXiv

A Parameter Space Exploration of High Resolution Numerically Evolved Early Type Galaxies Including AGN Feedback and Accurate Dynamical Treatment of Stellar Orbits

An extensive exploration of the model parameter space of axisymmetric Early-Type Galaxies (ETGs) hosting a central supermassive Black Hole (SMBH) is conducted by means of high resolution hydrodynamical simulations performed with our code MACER. Global properties such as 1) total SMBH accreted mass, 2) final X-ray luminosity and temperature of the X-ray emitting halos, 3) total amount of new stars formed from the cooling gas, 4) total ejected mass in form of supernovae and AGN feedback induced galactic winds, are obtained as a function of galaxy structure and internal dynamics. In addition to the galactic dark matter halo, the model galaxies are also embedded in a group/cluster dark matter halo; finally cosmological accretion is also included, with amount and time dependence derived from cosmological simulations. Angular momentum conservation leads to the formation of cold HI disks; these disks further evolve under the action of star formation induced by disk instabilities, of the associated mass discharge onto the central SMBH, and of the consequent AGN feedback. At the end of the simulations, the hot (metal enriched) gas mass is roughly $10\%$ the mass in the old stars, with twice as much having been ejected into the intergalactic medium. The cold gas disks are a $\approx$ kpc in size, and the metal rich new stars are in $0.1$ kpc disks. The masses of cold gas and new stars are roughly $0.1\%$ the mass of the old stars. Overall, the final systems appear to reproduce quite successfully the main global properties of real ETGs.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Massive black hole formation in dense stellar environments: Enhanced X-ray detection rates in high velocity dispersion nuclear star clusters

We analyze Chandra X-ray Observatory imaging of 108 galaxies hosting nuclear star clusters (NSCs) to search for signatures of massive black holes (BHs). NSCs are extremely dense stellar environments with conditions that can theoretically facilitate massive BH formation. Recent work by Stone et al. (2017) finds that sufficiently dense NSCs should be unstable to the runaway growth of a stellar mass BH into a massive BH via tidal captures. Furthermore, there is a velocity dispersion threshold ($40\;\rm{km\;s^{-1}}$) above which NSCs should inevitably form a massive BH. To provide an observational test of these theories, we measure X-ray emission from NSCs and compare to the measured velocity dispersion and tidal capture runaway timescale. We find that NSCs above the $40\;\rm{km\;s^{-1}}$ threshold are X-ray detected at roughly twice the rate of those below (after accounting for contamination from X-ray binaries). These results are consistent with a scenario in which dense, high-velocity NSCs can form massive BHs, providing a formation pathway that does not rely on conditions found only at high redshift.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Metallicity Distribution Function in Outer Halo Fields of Simulated Elliptical Galaxies Compared to Observations of NGC 5128

Stellar metallicity distribution functions (MDF) have been measured for resolved stellar populations in the outer halos of many galaxies in nearby groups. Among them, the MDF of NGC 5128, the central giant elliptical in the Centaurus group, provides essential constraints for theories of massive galaxy formation and hierarchical assembly. To investigate the formation and chemical evolution history of the outer halo of giant elliptical galaxies, we examine the chemical properties of three zoom-in high resolution cosmological hydrodynamical simulations of an NGC 5128-like giant elliptical galaxy and compare their outer halo MDFs to the observed one of NGC 5128. Even though the simulated galaxies have different merging histories and age distributions, all predicted MDFs are in good qualitative agreement with the observed one. The median metallicity of the simulated galaxies is on average $\rm [M/H]=-0.41 \pm 0.06$ compared to the observed value of $\rm [M/H]=-0.38 \pm 0.02$ for NGC 5128, and the dispersion in metallicity is $\sim 0.77$ dex for both observed and simulated galaxies. We investigate the origin of the stars ending up in the outer halo field of simulated galaxies and show that most have an `accreted' origin, formed in other small galaxies and later accreted in mergers. Only $\sim 15$ percent of the stars are formed `in situ' within the main progenitor of galaxy and radially migrate outwards. We show that the contribution of metal-rich in situ stars is sub-dominant in the outer halos of our simulated galaxies, but can be prominent in the inner regions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Infra-Red Emission from Cold Gas Dusty Disks in Massive Ellipticals

What is the expected infrared output of elliptical galaxies? Here we report the latest findings obtained in this high time resolution (~10 years) and high spatial resolution (2.5 parsec at center) study. We add a set of grain physics to the MACER code, including (a) dust grains made in passive stellar evolution; (b) dust grain growth due to collision and sticking; (c) grain destruction due to thermal sputtering; (d) dust cooling of hot gas via inelastic collisions; and (e) radiation pressure on dust grains. The code improvements enable us to analyze the effects of dust on metal depletion and AGN obscuration, and also to assess the infrared output of the modeled galaxies. We simulate a representative massive elliptical galaxy of a central stellar velocity dispersion ~ 260 km/s and modest rotation. We find that: (1) the circumnuclear disk (of a size <~ 1 kpc) is dusty in its outer region where most of the metals are in dust grains, while in the inner disk most of the dust grains are destroyed by the AGN irradiation; (2) the dusty disk is optically thick to both the starlight within the disk and the radiation from the central AGN. Thus the AGN will be obscured behind the disk, and the latter is of a covering factor ~ 0.2; (3) the dust infrared emission is mainly due to the AGN irradiation. The median infrared luminosity is ~ 2e44 erg/s, and it can reach >~ 1e46 erg/s during outbursts; (4) the duty cycles of the AGN activities, star formation, and the dust infrared luminosity roughly match observations, e.g., in most of its lifetime, the simulated galaxy is a stereotypical "quiescent" elliptical galaxy with L_{IR} ~ 1e11*L_{solar}, while the star formation rate can exceed 250 M_{solar}/yr during central outbursts.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

NIHAO-UHD: The properties of MW-like stellar disks in high resolution cosmological simulations

Simulating thin and extended galactic disks has long been a challenge in computational astrophysics. We introduce the NIHAO-UHD suite of cosmological hydrodynamical simulations of Milky Way mass galaxies and study stellar disk properties such as stellar mass, size and rotation velocity which agree well with observations of the Milky Way and local galaxies. In particular, the simulations reproduce the age-velocity dispersion relation and a multi-component stellar disk as observed for the Milky Way. Half of our galaxies show a double exponential vertical profile, while the others are well described by a single exponential model which we link to the disk merger history. In all cases, mono-age populations follow a single exponential whose scale height varies monotonically with stellar age and radius. The scale length decreases with stellar age while the scale height increases. The general structure of the stellar disks is already set at time of birth as a result of the inside-out and upside-down formation. Subsequent evolution modifies this structure by increasing both the scale length and height of all mono-age populations. Thus, our results put tight constraints on how much dynamical memory stellar disks can retain over cosmological timescales. Our simulations demonstrate that it is possible to form thin galactic disks in cosmological simulations provided there are no significant stellar mergers at low redshifts. Most of the stellar mass is formed in-situ with only a few percent ($\lesssim5\%$) brought in by merging satellites at early times. Redshift zero snapshots and halo catalogues are publicly available.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

MIND THE GAP: Is The Too Big To Fail Problem Resolved?

The faintness of satellite systems in galaxy groups has contributed to the widely discussed "missing satellite" and "too big to fail" issues. Using techniques based on Tremaine & Richstone (1977), we show that there is no problem with the luminosity function computed from modern codes per se, but that the gap between first and second brightest systems is too big {\it given} the luminosity function, that the same large gap is found in modern, large scale baryonic $Λ$CDM simulations such as EAGLE and IllustrisTNG, is even greater in dark matter only simulations, and finally, that this is most likely due to gravitationally induced merging caused by classical dynamical friction. Quantitatively the gap is larger in the computed simulations than in the randomized ones by $1.79 \pm 1.04$, $1.51 \pm 0.93$, $3.43 \pm 1.44$ and $3.33 \pm 1.35$ magnitudes in the EAGLE, IllustrisTNG, and dark matter only simulations of EAGLE and IllustrisTNG respectively. Furthermore the anomalous gaps in the simulated systems are even larger than in the real data by over half a magnitude and are still larger in the dark matter only simulations. Briefly stated, $Λ$CDM does not have a problem with an absence of "too big to fail" galaxies. Statistically significant large gaps between first and second brightest galaxies are to be expected.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

On the Role of Hot Feedback Mode in Active Galactic Nuclei Feedback in an Elliptical Galaxy

Depending on the value of the accretion rate, black hole accretion is divided into cold and hot modes. The two modes have distinctly different physics and correspond to two feedback modes. Most previous feedback works either focus only on one mode, or the accretion physics is not always properly adopted. Here, by performing hydrodynamical numerical simulations of AGN feedback in an elliptical galaxy, we show that including both is important, and gives different results from including just one or the other. We specifically focus on the wind and radiation (but neglecting the jet) feedback in the hot mode, to explore whether this particular mode of feedback can play any important feedback role. For this aim, we have run two test models. In one model, we always adopt the cold mode no matter what value the accretion rate is; in another model, we turn off the AGN once it enters into the hot mode. We have calculated the AGN light curves, black hole growth, star formation, and AGN energy duty-cycle; and compared the results with the model in which the two modes are correctly included. Important differences are found. For example, if we were to adopt only the cold mode, the total mass of newly formed stars would become two orders of magnitude smaller, and the fraction of energy ejected within the high accretion regime (i.e., $L_{\rm BH} > 2\% L_{\rm Edd}$) would be too small to be consistent with observations. We have also investigated the respective roles of wind and radiation in the hot mode.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Hyper-Eddington accretion flows onto massive black holes

We study very-high rate spherically symmetric accretion flows onto a massive black hole (BH; 10^2 < M_BH < 10^6 Msun) embedded in a dense gas cloud with a low abundance of metals, performing one-dimensional hydrodynamical simulations which include multi-frequency radiation transfer and non-equilibrium primordial chemistry. We find that rapid gas supply from the Bondi radius at a hyper-Eddington rate can occur without being impeded by radiation feedback when (n/10^5 cm^-3) > (M_BH/10^4Msun)^{-1}(T/10^4 K)^{3/2}, where n and T are the density and temperature of ambient gas outside of the Bondi radius. The resulting accretion rate in this regime is steady, and larger than 3000 times the Eddington rate. At lower Bondi rates, the accretion is episodic due to radiative feedback and the average rate is limited below the Eddington rate. For the hyper-Eddington case, the steady solution consists of two parts: a radiation-dominated central core, where photon trapping due to electron scattering is important, and an accreting envelope which follows a Bondi profile with T~8000 K. When the emergent luminosity is limited below the Eddington luminosity because of photon trapping, radiation from the central region does not affect the gas dynamics at larger scales. We apply our result to the rapid formation of massive BHs in protogalaxies with a virial temperature of T_vir> 10^4 K. Once a seed BH forms at the center of the galaxy, it can grow up to a maximum ~10^5 (T_vir/10^4 K) Msun via gas accretion independent of the initial BH mass. Finally, we discuss possible observational signatures of rapidly accreting BHs with/without allowance for dust. We suggest that these systems could explain Lya emitters without X-rays and luminous infrared sources with hot dust emission, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Launching cosmic-ray-driven outflows from the magnetized interstellar medium

We present a hydrodynamical simulation of the turbulent, magnetized, supernova (SN)-driven interstellar medium (ISM) in a stratified box that dynamically couples the injection and evolution of cosmic rays (CRs) and a self-consistent evolution of the chemical composition. CRs are treated as a relativistic fluid in the advection-diffusion approximation. The thermodynamic evolution of the gas is computed using a chemical network that follows the abundances of H+, H, H2, CO, C+, and free electrons and includes (self-)shielding of the gas and dust. We find that CRs perceptibly thicken the disk with the heights of 90% (70%) enclosed mass reaching ~1.5 kpc (~0.2 kpc). The simulations indicate that CRs alone can launch and sustain strong outflows of atomic and ionized gas with mass loading factors of order unity, even in solar neighborhood conditions and with a CR energy injection per SN of 10^50 erg, 10% of the fiducial thermal energy of an SN. The CR-driven outflows have moderate launching velocities close to the midplane (~100 km/s) and are denser (ρ~1e-24 - 1e-26 g/cm^3), smoother, and colder than the (thermal) SN-driven winds. The simulations support the importance of CRs for setting the vertical structure of the disk as well as the driving of winds.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Towards Self-Consistent Modelling of the Sgr A* Accretion Flow: Linking Theory and Observation

The interplay between supermassive black holes (SMBHs) and their environments is believed to command an essential role in galaxy evolution. The majority of these SMBHs are in the radiative inefficient accretion phase where this interplay remains elusive, but suggestively important, due to few observational constraints. To remedy this, we directly fit 2-D hydrodynamic simulations to Chandra observations of Sgr A* with Markov Chain Monte Carlo sampling, self-consistently modelling the 2-D inflow-outflow solution for the first time. We find the temperature and density at flow onset are consistent with the origin of the gas in the stellar winds of massive stars in the vicinity of Sgr A*. We place the first observational constraints on the angular momentum of the gas and estimate the centrifugal radius, r$_c$ $\approx$ 0.056 r$_b$ $\approx8\times10^{-3}$ pc, where r$_b$ is the Bondi radius. Less than 1\% of the inflowing gas accretes onto the SMBH, the remainder being ejected in a polar outflow. For the first time, we decouple the quiescent point-like emission from the spatially extended flow. We find this point-like emission, accounting for $\sim4$\% of the quiescent flux, is spectrally too steep to be explained by unresolved flares, nor Bremsstrahlung, but is likely a combination of a relatively steep synchrotron powerlaw and the high energy tail of inverse-Compton emission. With this self-consistent model of the accretion flow structure, we make predictions for the flow dynamics and discuss how future X-ray spectroscopic observations can further our understanding of the Sgr A* accretion flow.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Dynamical Potential-Density Pair for Star Clusters With Nearly Isothermal Interiors

We present a potential-density pair designed to model nearly isothermal star clusters (and similar self-gravitating systems) with a central core and an outer turnover radius, beyond which density falls off as $r^{-4}$. In the intermediate zone, the profile is similar to that of an isothermal sphere (density $ρ\propto r^{-2}$), somewhat less steep than the King 62 profile, and with the advantage that many dynamical quantities can be written in a simple closed form. We derive new analytic expressions for the cluster binding energy and velocity dispersion, and apply these to create toy models for cluster core collapse and evaporation. We fit our projected surface brightness profiles to observed globular and open clusters, and find that the quality of the fit is generally at least as good as that for the surface brightness profiles of King 62. This model can be used for convenient computation of the dynamics and evolution of globular and nuclear star clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Two-Parameter Matching Scheme for Massive Galaxies and Dark Matter Haloes

Halo Abundance Matching has been used to construct a one-parameter mapping between galaxies and dark matter haloes by assuming that halo mass and galaxy luminosity (or stellar mass) are monotonically related. While this approach has been reasonably successful, it is known that galaxies must be described by at least two parameters, as can be seen from the two-parameter Fundamental Plane on which massive early-type galaxies lie. In this paper, we derive a connection between initial dark matter density perturbations in the early universe and present-day virialized dark matter haloes by assuming simple spherical collapse combined with conservation of mass and energy. We find that $z = 0$ halo concentration, or alternatively the inner slope of the halo density profile $α$, is monotonically and positively correlated with the collapse redshift of the halo. This is qualitatively similar to the findings of some previous works based on numerical simulations, with which we compare our results. We then describe how the halo mass and concentration (or inner slope $α$) can be used as two halo parameters in combination with two parameters of early-type galaxies to create an improved abundance matching scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Kinetic energy from supernova feedback in high-resolution galaxy simulations

We describe a new method for adding a prescribed amount of kinetic energy to simulated gas modeled on a cartesian grid by directly altering grid cells' mass and velocity in a distributed fashion. The method is explored in the context of supernova feedback in high-resolution ($\sim 10$ pc) hydrodynamic simulations of galaxy formation. Resolution-dependence is a primary consideration in our application of the method and simulations of isolated explosions (performed at different resolutions) motivate a resolution-dependent scaling for the injected fraction of kinetic energy that we apply in cosmological simulations of a $10^9$ Msun dwarf halo. We find that in high density media ($\gtrsim$ 50 cm$^{-3}$) with coarse resolution ($\gtrsim 4$ pc per cell), results are sensitive to the initial kinetic energy fraction due to early and rapid cooling. In our galaxy simulations, the deposition of small amounts of supernova energy in kinetic form (as little as 1%) has a dramatic impact on the evolution of the system, resulting in an order of magnitude suppression of stellar mass. The overall behavior of the galaxy in the two highest resolution simulations we perform appears to converge. We discuss the resulting distribution of stellar metallicities, an observable sensitive to galactic wind properties, and find that while the new method demonstrates increased agreement with observed systems, significant discrepancies remain, likely due to simplistic assumptions that neglect contributions from Type Ia supernovae and stellar winds.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Supernova Feedback and the Hot Gas Filling Fraction of the Interstellar Medium

Supernovae (SN), the most energetic stellar feedback mechanism, are crucial for regulating the interstellar medium (ISM) and launching galactic winds. We explore how supernova remnants (SNRs) create a multiphase medium by performing 3D hydrodynamical simulations at various SN rates, $S$, and ISM average densities, $\bar{n}$. The evolution of a SNR in a self-consistently generated three-phase ISM is qualitatively different from that in a uniform or a two-phase warm/cold medium. By travelling faster and further in the low-density hot phase, the domain of a SNR increases by $>10^{2.5}$. Varying $\bar{n}$ and $S$, we find that a steady state can only be achieved when the hot gas volume fraction $f_{\rm{V,hot}}\lesssim 0.6 \pm 0.1 $. Above that level, overlapping SNRs render connecting topology of the hot gas, and the ISM is subjected to thermal runaway. Photoelectric heating (PEH) has a surprisingly strong impact on $f_{\rm{V,hot}}$. For $\bar{n}\gtrsim 3 \cm-3 $, a reasonable PEH rate is able to suppress the thermal runaway. Overall, we determine the critical SN rate for the onset of thermal runaway to be $S_{\rm{crit}} = 200 (\bar{n}/1\cm-3)^k (E_{\rm{SN}}/10^{51}\erg)^{-1} \kpc^{-3} \myr-1$, where $k = (1.2,2.7)$ for $\bar{n} \leq 1$ and $> 1\cm-3 $, respectively. We present a fitting formula of the ISM pressure $P(\bar{n}$, $S$), which can be used as an effective equation of state in cosmological simulations. Despite the 5 orders of magnitude span of $(\bar{n},S)$, the average Mach number varies little: $\mathcal{M} \approx \ 0.5\pm 0.2, \ 1.2\pm 0.3,\ 2.3\pm 0.9$ for the hot, warm and cold phases, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The impact of mechanical AGN feedback on the formation of massive early-type galaxies

We employ cosmological hydrodynamical simulations to investigate the effects of AGN feedback on the formation of massive galaxies with present-day stellar masses of $M_{stel} = 8.8 \times 10^{10} - 6.0 \times 10^{11} M_{sun}$. Using smoothed particle hydrodynamics simulations with a pressure-entropy formulation that allows an improved treatment of contact discontinuities and fluid mixing, we run three sets of simulations of 20 halos with different AGN feedback models: (1) no feedback, (2) thermal feedback, and (3) mechanical and radiation feedback. We assume that seed black holes are present at early cosmic epochs at the centre of emerging dark matter halos and trace their mass growth via gas accretion and mergers with other black holes. Both feedback models successfully recover the observed M_BH - sigma relation and black hole-to-stellar mass ratio for simulated central early-type galaxies. The baryonic conversion efficiencies are reduced by a factor of two compared to models without any AGN feedback at all halo masses. However, massive galaxies simulated with thermal AGN feedback show a factor of ~10-100 higher X-ray luminosities than observed. The mechanical/radiation feedback model reproduces the observed correlation between X-ray luminosities and velocity dispersion, e.g. for galaxies with sigma = 200 km/s, the X-ray luminosity is reduced from $10^{42}$ erg/s to $10^{40}$ erg/s. It also efficiently suppresses late time star formation, reducing the specific star formation rate from $10^{-10.5}$ $yr^{-1}$ to $10^{-14}$ $yr^{-1}$ on average and resulting in quiescent galaxies since z=2, whereas the thermal feedback model shows higher late time in-situ star formation rates than observed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The stellar accretion origin of stellar population gradients in massive galaxies at large radii

We investigate the evolution of stellar population gradients from $z=2$ to $z=0$ in massive galaxies at large radii ($r > 2R_{\mathrm{eff}}$) using ten cosmological zoom simulations of halos with $6 \times 10^{12} M_{\odot} < M_{\mathrm{halo}} < 2 \times 10^{13}M_{\odot}$. The simulations follow metal cooling and enrichment from SNII, SNIa and AGB winds. We explore the differential impact of an empirical model for galactic winds that reproduces the mass-metallicity relation and its evolution with redshift. At larger radii the galaxies, for both models, become more dominated by stars accreted from satellite galaxies in major and minor mergers. In the wind model, fewer stars are accreted, but they are significantly more metal poor resulting in steep global metallicity ($\langle \nabla Z_{\mathrm{stars}} \rangle= -0.35$ dex/dex) and color (e.g. $\langle \nabla g-r \rangle = -0.13$ dex/dex) gradients in agreement with observations. In contrast, colour and metallicity gradients of the models without winds are inconsistent with observations. Age gradients are in general mildly positive at $z=0$ ($\langle \nabla Age_{\mathrm{stars}} \rangle= 0.04$ dex/dex) with significant differences between the models at higher redshift. We demonstrate that for the wind model, stellar accretion is steepening existing in-situ metallicity gradients by about 0.2 dex by the present day and helps to match observed gradients of massive early-type galaxies at large radii. Colour and metallicity gradients are significantly steeper for systems which have accreted stars in minor mergers, while galaxies with major mergers have relatively flat gradients, confirming previous results. This study highlights the importance of stellar accretion for stellar population properties of massive galaxies at large radii, which can provide important constraints for formation models.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

AGN feedback in an isolated elliptical galaxy: the effect of strong radiative feedback in the kinetic mode

Based on two-dimensional high resolution hydrodynamic numerical simulation, we study the mechanical and radiative feedback effects from the central AGN on the cosmological evolution of an isolated elliptical galaxy. Physical processes such as star formation and supernovae are considered. The inner boundary of the simulation domain is carefully chosen so that the fiducial Bondi radius is resolved and the accretion rate of the black hole is determined self-consistently. In analogy to previous works, we assume that the specific angular momentum of the galaxy is low. It is well-known that when the accretion rates are high and low, the central AGNs will be in cold and hot accretion modes, which correspond to the radiative and kinetic feedback modes, respectively. The emitted spectrum from the hot accretion flows is harder than that from the cold accretion flows, which could result in a higher Compton temperature accompanied by a more efficient radiative heating, according to previous theoretical works. Such a difference of the Compton temperature between the two feedback modes, the focus of this study, has been neglected in previous works. Significant differences in the kinetic feedback mode are found as a result of the stronger Compton heating and accretion becomes more chaotic. More importantly, if we constrain models to correctly predict black hole growth and AGN duty cycle after cosmological evolution, we find that the favored model parameters are constrained: mechanical feedback efficiency diminishes with decreasing luminosity (the maximum efficiency being $\simeq 10^{-3.5}$) and X-ray Compton temperature increases with decreasing luminosity, although models with fixed mechanical efficiency and Compton temperature can be found that are satisfactory as well. We conclude that radiative feedback in the kinetic mode is much more important than previously thought.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Co-Evolution of Galactic Nuclei and Globular Cluster Systems

We revisit the hypothesis that dense galactic nuclei are formed from inspiraling globular clusters. Recent advances in understanding of the continuous formation of globular clusters over cosmic time and the concurrent evolution of the galaxy stellar distribution allow us to construct a simple model that matches the observed spatial and mass distributions of clusters in the Galaxy and the giant elliptical galaxy M87. In order to compare with observations, we model the effects of dynamical friction and dynamical evolution, including stellar mass loss, tidal stripping of stars, and tidal disruption of clusters by the growing galactic nucleus. We find that inspiraling globular clusters form a dense central structure, with mass and radius comparable to the typical values in observed nuclear star clusters (NSCs) in late-type and low-mass early-type galaxies. The density contrast associated with the NSC is less pronounced in giant elliptical galaxies. Our results indicate that the NSC mass as a fraction of mass of the galaxy stellar spheroid scales as M_NSC/M_* = 0.0025 M_{*,11}^{-0.5}. Thus disrupted globular clusters could contribute most of the mass of NSCs in galaxies with stellar mass below 10^{11} Msun. The inner part of the accumulated cluster may seed the growth of a central black hole via stellar dynamical core collapse, thereby relieving the problem of how to form luminous quasars at high redshift. The seed black hole may reach about 10^5 Msun within 1 Gyr of the beginning of globular cluster formation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Consequences of Mechanical and Radiative Feedback from Black Holes in Disc Galaxy Mergers

We study the effect of AGN mechanical and radiation feedback on the formation of bulge dominated galaxies via mergers of disc galaxies. The merging galaxies have mass-ratios of 1:1 to 6:1 and include pre-existing hot gaseous halos to properly account for the global impact of AGN feedback. Using smoothed particle hydrodynamics simulation code (GADGET-3) we compare three models with different AGN feedback models: (1) no black hole and no AGN feedback; (2) thermal AGN feedback; and (3) mechanical and radiative AGN feedback. The last model is motivated by observations of broad line quasars which show winds with initial velocities of $v_w \ge$ 10,000 km/s and also heating associated with the central AGN X-ray radiation. The primary changes in gas properties due to mechanical AGN feedback are lower thermal X-ray luminosity from the final galaxy - in better agreement with observations - and galactic outflows with higher velocity $\sim 1000$ km/s similar to recent direct observations of nearby merger remnants. The kinetic energy of the outflowing gas is a factor of $\sim$ 20 higher than in the thermal feedback case. All merger remnants with momentum-based AGN feedback with $v_w \sim 10,000$ km/s and $ε_w=2 \times 10^{-3}$, independent of their progenitor mass-ratios, reproduce the observed relations between stellar velocity dispersion and black hole mass ($M_{\rm bh} - σ$) as well as X-ray luminosity ($L_X - σ$) with $10^{37.5}$ erg/s $\lesssim L_X (0.3-8~{\rm keV}) \lesssim 10^{39.2}$ erg/s for velocity dispersions in the range of 120 km/s $\lesssim σ\lesssim$ 190 km/s. In addition, the mechanical feedback produces a much greater AGN variability. We also show that gas is more rapidly and impulsively stripped from the galactic centres driving a moderate increase in galaxy size and decrease in central density with the mechanical AGN feedback model.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Grain Physics and IR Dust Emission in AGN Environments

We study the effects of a detailed dust treatment on the properties and evolution of early-type galaxies containing central black holes, as determined by AGN feedback. We find that during cooling flow episodes, radiation pressure on the dust in and interior to infalling shells of cold gas can greatly impact the amount of gas able to be accreted and therefore the frequency of AGN bursts. However, the overall hydrodynamic evolution of all models, including mass budget, is relatively robust to the assumptions on dust. We find that IR re-emission from hot dust can dominate the bolometric luminosity of the galaxy during the early stages of an AGN burst, reaching values in excess of $10^{46}$ erg/s. The AGN-emitted UV is largely absorbed, but the optical depth in the IR does not exceed unity, so the radiation momentum input never exceeds $L_{\rm BH}/c$. We constrain the viability of our models by comparing the AGN duty cycle, broadband luminosities, dust mass, black hole mass, and other model predictions to current observations. These constraints force us to models wherein the dust to metals ratios are $\simeq 1%$ of the Galactic value, and only models with a dynamic dust to gas ratio are able to produce both quiescent galaxies consistent with observations and high obscured fractions during AGN "on" phases. During AGN outbursts, we predict that a large fraction of the FIR luminosity can be attributed to warm dust emission ($\simeq100$ K) from dense dusty gas within $\leq 1$ kpc reradiating the AGN UV emission.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Astronomical Reach of Fundamental Physics

Using basic physical arguments, we derive by dimensional and physical analysis the characteristic masses and sizes of important objects in the Universe in terms of just a few fundamental constants. This exercise illustrates the unifying power of physics and the profound connections between the small and the large in the Cosmos we inhabit. We focus on the minimum and maximum masses of normal stars, the corresponding quantities for neutron stars, the maximum mass of a rocky planet, the maximum mass of a white dwarf, and the mass of a typical galaxy. To zeroth order, we show that all these masses can be expressed in terms of either the Planck mass or the Chandrasekar mass, in combination with various dimensionless quantities. With these examples we expose the deep interrelationships imposed by Nature between disparate realms of the Universe and the amazing consequences of the unifying character of physical law.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Understanding black hole mass assembly via accretion and mergers at late times in cosmological simulations

Accretion is thought to primarily contribute to the mass accumulation history of supermassive black holes throughout cosmic time. While this may be true at high redshifts, at lower redshifts and for the most massive black holes mergers themselves might add significantly to the mass budget. We evolve SMBHs from $4 > z > 0$ using merger trees derived from hydrodynamical cosmological simulations of a cluster and void region, scaled to the observed value of the stellar mass fraction to account for overcooling. Mass gains from gas accretion proportional to bulge growth and BH-BH mergers are tracked, as are black holes that remain "orbiting" due to insufficient dynamical friction in a merger remnant, as well as those that are ejected due to gravitational recoil. We find that gas accretion remains the dominant source of mass accumulation in almost all SMBHs; mergers contribute $2.5\pm0.1\%$ for all SMBHs in the cluster and $1.0\pm0.1\%$ in the void since $z = 4$. However, mergers are significant for massive SMBHs. The fraction of mass accumulated from mergers for central BHs generally increases for larger values of the host bulge mass: in the void, the fraction is $2\%$ at $M_{*, bul} = 10^{10} M_{\odot}$, increasing to $4\%$ at $M_{*, bul} \gtrsim 10^{11} M_{\odot}$, and in the cluster it is $4\%$ at $M_{*, bul} = 10^{10} M_{\odot}$ and $23\%$ at $10^{12} M_{\odot}$. We find that $40\%$ of SMBHs and $\approx 8\%$ of the total SMBH mass is found orbiting in the cluster region at $z = 0$. The existence of orbiting and ejected SMBHs requires modification of the Soltan argument. We estimate this correction to the integrated accreted mass density of SMBHs to be in the range $6-21\%$, with a mean value of $11\pm3\%$. We also calculate the total energy output and strain from gravitational waves emitted by merging SMBHs, and obtain a signal potentially detectable by pulsar timing arrays.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Radiation driven outflow in active galactic nuclei: the feedback effects of scattered and reprocessed photons

We perform time-dependent, 2DHD numerical simulations to study the dynamics of a slowly rotating accretion flow from sub-pc to pc scales under the irradiation from the central AGN. Compared to previous work, we improve the calculation of the radiative force due to X-rays. More importantly, in addition to radiative pressure and radiative heating/cooling directly from the central AGN, in the momentum equation we also include the force due to the scattered and reprocessed photons. We find that the accretion flow properties change significantly due to this "re-radiation" effect. The inflow rate at the inner boundary is reduced, while the outflow rate at the outer boundary is enhanced by about one order of magnitude. This effect is more significant when the density at the outer boundary is higher. The properties of outflows such as velocity, momentum and energy fluxes, and the ratio of outflow rate and the accretion rate, are calculated. We find that the efficiency of transferring the radiation power into the kinetic power of outflow is typically $10^{-3}$, far below the value of $\sim 0.05$ which is assumed in some cosmological simulations. The effect of the temperature of the gas at the outer boundary ($T_0$) is investigated. When $T_0$ is high, the emitted luminosity of the accretion flow oscillates. This is because in this case the gas around the Bondi radius can be more easily heated to be above the virial temperature due to its high internal energy. Another question we hope to address is the so-called "sub-Eddington" puzzle. Observationally, the luminosity of almost all AGNs are sub-Eddington, while theoretically the luminosity of an accretion flow can easily be super-Eddington. We find that even when the re-radiation effect is included and outflow does become much stronger, the luminosity, while reduced, can still be super-Eddington.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The effect of metal enrichment and galactic winds on galaxy formation in cosmological zoom simulations

We investigate the differential effects of metal cooling and galactic stellar winds on the cosmological formation of individual galaxies with three sets of cosmological, hydrodynamical zoom simulations of 45 halos in the mass range 10^11<M_halo<10^13M_sun. Models including both galactic winds and metal cooling (i) suppress early star formation at z>1 and predict reasonable star formation histories, (ii) produce galaxies with high cold gas fractions (30-60 per cent) at high redshift, (iii) significantly reduce the galaxy formation efficiencies for halos (M_halo<10^12M_sun) at all redshifts in agreement with observational and abundance matching constraints, (iv) result in high-redshift galaxies with reduced circular velocities matching the observed Tully-Fisher relation at z~2, and (v) significantly increase the sizes of low-mass galaxies (M_stellar<3x10^10M_sun) at high redshift resulting in a weak size evolution - a trend in agreement with observations. However, the low redshift (z<0.5) star formation rates of massive galaxies are higher than observed (up to ten times). No tested model predicts the observed size evolution for low-mass and high-mass galaxies simultaneously. Due to the delayed onset of star formation in the wind models, the metal enrichment of gas and stars is delayed and agrees well with observational constraints. Metal cooling and stellar winds are both found to increase the ratio of in situ formed to accreted stars - the relative importance of dissipative vs. dissipationless assembly. For halo masses below ~10^12M_sun, this is mainly caused by less stellar accretion and compares well to predictions from semi-analytical models but still differs from abundance matching models. For higher masses, the fraction of in situ stars is over-predicted due to the unrealistically high star formation rates at low redshifts.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Calibration of Nonthermal Pressure in Global Dark Matter Simulations of Clusters of Galaxies

We present a new method for incorporating nonthermal pressure from bulk motions of gas into an analytic model of the intracluster medium in clusters of galaxies, which is based on a polytropic equation of state and hydrostatic equilibrium inside gravitational potential wells drawn from cosmological dark matter simulations. The pressure is allowed to have thermal and nonthermal components with different radial distributions; the overall level of nonthermal support is based on the dynamical state of the halo, such that it is lower in more relaxed clusters. This level is normalized by comparison to pressure profiles derived from X-ray observations, and to a high resolution hydrodynamical simulation. The nonthermal pressure fraction measured at r_500 is typically in the range 10-20%, increasing with cluster mass and with redshift. The resulting model cluster properties are in accord with Sunyaev-Zel'dovich (SZ) effect observations of clusters. Inclusion of nonthermal pressure reduces the expected angular power spectrum of SZ fluctuations in the microwave sky by 24%.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Cosmic Mach Number: A Sensitive Probe for the Growth of Structure

In this Letter, we investigate the potential power of the Cosmic Mach Number (CMN), which is the ratio between the mean velocity and the velocity dispersion of galaxies as a function of cosmic scales, to constrain cosmologies. We first measure the CMN from 5 catalogues of galaxy peculiar velocity surveys at low redshift (0.002<z<0.03), and use them to contrast cosmological models. Overall, current data is consistent with the WMAP7 LCDM model. We find that the CMN is highly sensitive to the growth of structure on scales 0.01<k<0.1 h/Mpc in Fourier space. Therefore, modified gravity models, and models with massive neutrinos, in which the structure growth generally deviates from that in the LCDM model in a scale-dependent way, can be well differentiated from the LCDM model using future CMN data.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Forming Early-Type Galaxies in LambdaCDM Simulations -I. Assembly histories

We present a sample of nine high resolution cosmological simulations in the mass range of M_vir=7x10^11-4x10^12 M_sun starting from LambdaCDM initial conditions. Our simulations include primordial radiative cooling, photoionization, star formation, supernova II feedback, but exclude supernova driven winds and AGN feedback. The simulated galaxies assemble in two phases, with the initial growth dominated by compact (r<r_eff) in situ star formation fueled by cold, low entropy gas streams resulting in a very similar mean assembly redshift of z_{f,ins}~2.5 for the in situ stellar component in all galaxies. The late growth is dominated by accretion of old stars formed in subunits outside the main galaxy (r>r_eff) resulting in an assembly redshift of z_{f,acc}~0.5-1.5 with much larger scatter. We find a positive correlation between the fraction of accreted stars and the final mass of our galaxies. We show that gravitational feedback strongly suppresses late star formation in massive galaxies contributing to the observed galaxy color bimodality. The accretion of stellar material is also responsible for the observed size growth of early-type galaxies. In addition, we find that the dark matter fractions within the stellar half-mass radii continuously increase towards lower redshift from about f_DM~0.05 at z~3 to f_DM~0.1-0.3 at z=0. Furthermore, the logarithmic slope of the total density profile is nearly isothermal at the present-day (gamma'~1.9-2.2). Finally, the input of gravitational heating lowers the central dark matter densities in the galaxies, with the effect being smaller compared to simulations without supernova feedback.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

How do minor mergers promote inside-out growth of ellipticals, transforming the size, density profile and dark matter fraction?

There is observational evidence for inside-out growth of elliptical galaxies since $z \gtrsim 2-3$, which is not driven by in-situ star formation. Many systems at high redshift have small sizes $\sim 1kpc$ and surface brightness profiles with low Sersic indices n. The most likely descendants have, on average, grown by a factor of two in mass and a factor of four in size, indicating $r \propto M^α$ with $α\gtrsim 2$. They also have surface brightness profiles with $n \gtrsim 5$. This evolution can be qualitatively explained on the basis of two assumptions: compact ellipticals predominantly grow by collisionless minor or intermediate 'dry' mergers, and they are embedded in massive dark matter halos. We draw these conclusions from idealized collisionless mergers spheroidal galaxies - with and without dark matter - with mass ratios of 1:1, 1:5, and 1:10. The sizes evolve as $r \propto M^α$ with $α< 2$ for mass-ratios of 1:1. For minor mergers of galaxies embedded in dark matter halos, the sizes grow significantly faster and the profile shapes change more rapidly. Mergers with moderate mass-ratios of 1:5 give $α\sim 2.3$ and a final Sersic index of $n = 9.5$ after doubling the stellar mass. This is accompanied by a significant increase of the dark matter fraction within the stellar half-mass radius, driven by the strong size increase probing larger, dark matter dominated regions. Only a few intermediate mass-ratio mergers of galaxies embedded in massive dark matter halos can result in the observed concurrent inside-out growth and the rapid evolution in profile shapes. Apart from negative stellar metallicity gradients such a 'minor' merger scenario also predicts significantly lower dark matter fractions for $z \sim 2$ compact quiescent galaxies and their rare present day analogues (abbreviated).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Radiative and Momentum Based Mechanical AGN Feedback in a 3-Dimensional Galaxy Evolution Code

We study the growth of black holes (BHs) in galaxies using three-dimensional smoothed particle hydrodynamic simulations with new implementations of the momentum mechanical feedback, and restriction of accreted elements to those that are gravitationally bound to the BH. We also include the feedback from the X-ray radiation emitted by the BH, which heats the surrounding gas in the host galaxies, and adds radial momentum to the fluid. We perform simulations of isolated galaxies and merging galaxies and test various feedback models with the new treatment of the Bondi radius criterion. We find that overall the BH growth is similar to what has been obtained by earlier workers using the Springel, Di Matteo, & Hernquist algorithms. However, the outflowing wind velocities and mechanical energy emitted by winds are considerably higher (v_w ~ 1000-3000 km/s) compared to the standard thermal feedback model (v_w ~ 50-100 km/s). While the thermal feedback model emits only 0.1 % of BH released energy in winds, the momentum feedback model emits more than 30 % of the total energy released by the BH in winds. In the momentum feedback model, the degree of fluctuation in both radiant and wind output is considerably larger than in the standard treatments. We check that the new model of the BH mass accretion agrees with analytic results for the standard Bondi problem.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Relaxation and Stripping: The Evolution of Sizes, Dispersions and Dark Matter Fractions in Major and Minor Mergers of Elliptical Galaxies

We revisit collisionless major and minor mergers of spheroidal galaxies in the context of the size evolution of elliptical galaxies. The simulations are performed as a series of mergers with mass-ratios of 1:1 and 1:10 for models representing pure bulges as well as bulges embedded in dark matter halos. For major and minor mergers, respectively, we identify and analyze two different processes, violent relaxation and stripping, leading to size evolution and a change of the dark matter fraction within the observable effective radius. Violent relaxation - which is the dominant mixing process for major mergers but less important for minor mergers - scatters relatively more dark matter particles than bulge particles to small. Stripping in minor mergers assembles stellar satellite particles at large radii in halo dominated regions of the massive host. This strongly increases the size of the bulge into regions with higher dark matter fractions leaving the inner host structure almost unchanged. A factor of two mass increase by minor mergers increases the dark matter fraction by 20 per cent. We present analytic corrections to simple one-component virial estimates for the evolution of the gravitational radii. If such a two-component system grows by minor mergers alone its size growth, $r_{\mathrm{e}} \propto M^α$, reaches values of $α\approx 2.4$, significantly exceeding the simple theoretical limit of $α= 2$. For major mergers the sizes grow with $α\lesssim 1$. Our results indicate that minor mergers of galaxies embedded in massive dark matter halos provide a potential mechanism for explaining the rapid size growth and the build-up of massive elliptical systems predicting significant dark matter fractions and radially biased velocity dispersions at large radii (abbreviated)

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The imminent detection of gravitational waves from massive black-hole binaries with pulsar timing arrays

Recent observations of massive galaxies indicate that they double in mass and quintuple in size between redshift z = 1 and the present, despite undergoing very little star formation, suggesting that galaxy mergers drive the evolution. Since these galaxies will contain supermassive black holes, this suggests a larger black hole merger rate, and therefore a larger gravitational-wave signal, than previously expected. We calculate the merger-driven evolution of the mass function, and find that merger rates are 10 to 30 times higher and gravitational waves are 3 to 5 times stronger than previously estimated, so that the gravitational-wave signal may already be detectable with existing data from pulsar timing arrays. We also provide an explanation for the disagreement with past estimates that were based on dark matter halo simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Environmental effects on the growth of super massive black holes and AGN feedback

We investigate how environmental effects by gas stripping alter the growth of a super massive black hole (SMBH) and its host galaxy evolution, by means of 1D hydrodynamical simulations that include both mechanical and radiative AGN feedback effects. By changing the truncation radius of the gas distribution (R_t), beyond which gas stripping is assumed to be effective, we simulate possible environments for satellite and central galaxies in galaxy clusters and groups. The continuous escape of gas outside the truncation radius strongly suppresses star formation, while the growth of the SMBH is less affected by gas stripping because the SMBH accretion is primarily ruled by the density of the central region. As we allow for increasing environmental effects - the truncation radius decreasing from about 410 to 50 kpc - we find that the final SMBH mass declines from about 10^9 to 8 x 10^8 Msol, but the outflowing mass is roughly constant at about 2 x 10^10 Msol. There are larger change in the mass of stars formed, which declines from about 2 x 10^10 to 2 x 10^9 Msol, and the final thermal X-ray gas, which declines from about 10^9 to 5 x 10^8 Msol, with increasing environmental stripping. Most dramatic is the decline in the total time that the objects would be seen as quasars, which declines from 52 Myr (for R_t = 377 kpc) to 7.9 Myr (for R_t = 51 kpc). The typical case might be interpreted as a red and dead galaxy having episodic cooling flows followed by AGN feedback effects resulting in temporary transitions of the overall galaxy color from red to green or to blue, with (cluster) central galaxies spending a much larger fraction of their time in the elevated state than do satellite galaxies.(Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The cosmological size and velocity dispersion evolution of massive early-type galaxies

We analyze 40 cosmological re-simulations of individual massive galaxies with present-day stellar masses of $M_{*} > 6.3 \times 10^{10} M_{\odot}$ in order to investigate the physical origin of the observed strong increase in galaxy sizes and the decrease of the stellar velocity dispersions since redshift $z \approx 2$. At present 25 out of 40 galaxies are quiescent with structural parameters (sizes and velocity dispersions) in agreement with local early type galaxies. At z=2 all simulated galaxies with $M_* \gtrsim 10^{11}M_{\odot}$ (11 out of 40) at z=2 are compact with projected half-mass radii of $\approx$ 0.77 ($\pm$0.24) kpc and line-of-sight velocity dispersions within the projected half-mass radius of $\approx$ 262 ($\pm$28) kms$^{-1}$ (3 out of 11 are already quiescent). Similar to observed compact early-type galaxies at high redshift the simulated galaxies are clearly offset from the local mass-size and mass-velocity dispersion relations. Towards redshift zero the sizes increase by a factor of $\sim 5-6$, following $R_{1/2} \propto (1+z)^α$ with $α= -1.44$ for quiescent galaxies ($α= -1.12$ for all galaxies). The velocity dispersions drop by about one-third since $z \approx 2$, following $σ_{1/2} \propto (1+z)^β$ with $β= 0.44$ for the quiescent galaxies ($β= 0.37$ for all galaxies). The simulated size and dispersion evolution is in good agreement with observations and results from the subsequent accretion and merging of stellar systems at $z\lesssim 2$ which is a natural consequence of the hierarchical structure formation. A significant number of the simulated massive galaxies (7 out of 40) experience no merger more massive than 1:4 (usually considered as major mergers). On average, the dominant accretion mode is stellar minor mergers with a mass-weighted mass-ratio of 1:5. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Effects of X-Ray Feedback from AGN on Host Galaxy Evolution

Hydrodynamic simulations of galaxies with active galactic nuclei (AGN) have typically employed feedback that is purely local: i.e., an injection of energy to the immediate neighborhood of the black hole. We perform GADGET-2 simulations of massive elliptical galaxies with an additional feedback component: an observationally calibrated X-ray radiation field which emanates from the black hole and heats gas out to large radii from the galaxy center. We find that including the heating and radiation pressure associated with this X-ray flux in our simulations enhances the effects which are commonly reported from AGN feedback. This new feedback model is twice as effective as traditional feedback at suppressing star formation, produces 3 times less star formation in the last 6 Gyr, and modestly lowers the final BH mass (30%). It is also significantly more effective than an X-ray background in reducing the number of satellite galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A New Test of the Statistical Nature of the Brightest Cluster Galaxies

A novel statistic is proposed to examine the hypothesis that all cluster galaxies are drawn from the same luminosity distribution (LD). In such a "statistical model" of galaxy LD, the brightest cluster galaxies (BCGs) are simply the statistical extreme of the galaxy population. Using a large sample of nearby clusters, we show that BCGs in high luminosity clusters (e.g., L_tot > 4x10^11 L_sun) are unlikely (probability <3x10^-4) to be drawn from the LD defined by all red cluster galaxies more luminous than M_r=-20. On the other hand, BCGs in less luminous clusters are consistent with being the statistical extreme. Applying our method to the second brightest galaxies, we show that they are consistent with being the statistical extreme, which implies that the BCGs are also distinct from non-BCG luminous, red, cluster galaxies. We point out some issues with the interpretation of the classical tests proposed by Tremaine & Richstone (1977) that are designed to examine the statistical nature of BCGs, investigate the robustness of both our statistical test and those of TR against difficulties in photometry of galaxies of large angular size, and discuss the implication of our findings on surveys that use the luminous red galaxies to measure the baryon acoustic oscillation features in the galaxy power spectrum.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Dissipational versus Dissipationless Galaxy Formation and the Dark Matter Content of Galaxies

We examine two extreme models for the build-up of the stellar component of luminous elliptical galaxies. In one case, we assume the build-up of stars is dissipational, with centrally accreted gas radiating away its orbital and thermal energy; the dark matter halo will undergo adiabatic contraction and the central dark matter density profile will steepen. For the second model, we assume the central galaxy is assembled by a series of dissipationless mergers of stellar clumps that have formed far from the nascent galaxy. In order to be accreted, these clumps lose their orbital energy to the dark matter halo via dynamical friction, thereby heating the central dark matter and smoothing the dark matter density cusp. The central dark matter density profiles differ drastically between these models. For the isolated elliptical galaxy, NGC 4494, the central dark matter densities follow the power-laws r^(-0.2) and r^(-1.7) for the dissipational and dissipationless models, respectively. By matching the dissipational and dissipationless models to observations of the stellar component of elliptical galaxies, we examine the relative contributions of dissipational and dissipationless mergers to the formation of elliptical galaxies and look for observational tests that will distinguish between these models. Comparisons to strong lensing brightest cluster galaxies yield median M*/L_B ratios of 2.1+/-0.8 and 5.2+/-1.7 at z=0.39 for the dissipational and dissipationless models, respectively. For NGC 4494, the best-fit dissipational and dissipationless models have M*/L_B=2.97 and 3.96. Comparisons to expected stellar mass-to-light ratios from passive evolution and population syntheses appear to rule out a purely dissipational formation mechanism for the central stellar regions of giant elliptical galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Feedback from central black holes in elliptical galaxies. II: Can purely mechanical energy feedback models work?

By using high-resolution 1D hydrodynamical simulations, we investigate the effects of purely mechanical feedback from super massive black holes (SMBHs) in the evolution of elliptical galaxies for a broad range of feedback efficiencies and compare the results to four major observational constraints. In particular, we focus on 1) the central black hole to stellar mass ratio of the host galaxy, 2) the lifetime of the luminous quasar phase, 3) the mass of stars formed in the host galaxy within the last Gyr, and 4) the X-ray luminosity of the hot diffuse gas. As a result, we try to pin down the most successful range of mechanical feedback efficiencies. We find that while low feedback efficiencies result in too much growth of the SMBH, high efficiencies totally blow out the hot interstellar gas, and the models are characterized by very low thermal X-ray luminosity well below the observed range. The net lifetime of the quasar phase is strongly coupled to the mass ratio between SMBH and its host galaxy, while the X-ray luminosity is generally correlated to the recent star formation within the last Gyr. When considering the popularly adopted model of the constant feedback efficiency, the feedback energy deposited into the ambient medium should be more than 0.01% of the SMBH accretion energy to be consistent with the SMBH mass to stellar mass ratio in the local universe. Yet, the X-ray luminosity of the hot gas favors about 0.005% of the accretion energy as the mechanical AGN feedback energy. We conclude that the purely mechanical feedback mode is unlikely to be simultaneously compatible with all four observable tests, even allowing a broad range of feedback efficiencies, and that including both radiative and mechanical feedback together may be a solution to comply the observational constraints. [abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Momentum Driving: which physical processes dominate AGN feedback?

The deposition of mechanical feedback from a supermassive black hole (SMBH) in an active galactic nucleus (AGN) into the surrounding galaxy occurs via broad-line winds which must carry mass and radial momentum as well as energy. The effect can be summarized by the dimensionless parameter $η=dot{M_outflow}/dot{M_accretion}= (2 ε_w c^2)/v_w^2$ where ($\epslion_w \equiv dot{E}_w/(dot{M_accretion} c^2)$) is the efficiency by which accreted matter is turned into wind energy in the disc surrounding the central SMBH. The outflowing mass and omentum are proportional to $η$, and many prior treatments have essentially assumed that $η=0$. We perform one- and two-dimensional simulations and find that the growth of the central SMBH is very sensitive to the inclusion of the mass and momentum driving but is insensitive to the assumed mechanical efficiency. For example in representative calculations, the omission of momentum and mass feedback leads to an hundred fold increase in the mass of the SMBH to over $10^{10} \Msun$. When allowance is made for momentum driving, the final SMBH mass is much lower and the wind efficiencies which lead to the most observationally acceptable results are relatively low with $ε_w \lesssim 10^{-4}$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Templates for the Sunyaev-Zel'dovich Angular Power Spectrum

We present templates for the Sunyaev-Zel'dovich (SZ) angular power spectrum based on four models for the nonlinear gas distribution. The frequency-dependent SZ temperature fluctuations, with thermal (TSZ) and kinetic (KSZ) contributions, are calculated by tracing through a dark matter simulation, processed to include gas in dark matter halos and in the filamentary intergalactic medium. Different halo gas models are compared to study how star formation, energetic feedback, and nonthermal pressure support influence the angular power spectrum. The standard model has been calibrated to reproduce the stellar and gas fractions and X-ray scaling relations measured from low redshift clusters and groups. The other models illustrate the current theoretical and empirical uncertainties relating to properties of the intracluster medium. Relative to the standard model, their angular power spectra differ by approximately 50% (TSZ), 20% (KSZ), and 40% (SZ at 148 GHz) for l=3000, sigma_8=0.8, and homogeneous reionization at z=10. The angular power spectrum decreases in amplitude as gas mass and binding energy is removed through star formation, and as gas is pushed out to larger radii by energetic feedback. With nonthermal pressure support, less pressure is required to maintain hydrostatic equilibrium, thus reducing the thermal contribution to the SZ power. We also calculate the SZ templates as a function of sigma_8 and quantify this dependence. Assuming C_l is proprotional to (sigma_8/0.8)^alpha, the effective scaling index ranges from 7<alpha_tsz<9, 4.5<alpha_ksz<5.5, and 6.5<alpha_sz(148 GHz)<8 at l=3000 for 0.6<sigma_8<1. The template spectra are publicly available and can be used when fitting for the SZ contribution to the cosmic microwave background on arcminute scales.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Two Phases of Galaxy Formation

Cosmological simulations of galaxy formation appear to show a two-phase character with a rapid early phase at z>2 during which in-situ stars are formed within the galaxy from infalling cold gas followed by an extended phase since z<3 during which ex-situ stars are primarily accreted. In the latter phase massive systems grow considerably in mass and radius by accretion of smaller satellite stellar systems formed at quite early times (z>3) outside of the virial radius of the forming central galaxy. These tentative conclusions are obtained from high resolution re-simulations of 39 individual galaxies in a full cosmological context with present-day virial halo masses ranging from 7e11 M_sun h^-1 < M_vir < 2.7e13 M_sun h^-1 and central galaxy masses between 4.5e10 M_sun h^-1 < M_* < 3.6e11 M_sun h^-1. The simulations include the effects of a uniform UV background, radiative cooling, star formation and energetic feedback from SNII. The importance of stellar accretion increases with galaxy mass and towards lower redshift. In our simulations lower mass galaxies ($M_* < 9e10 M_sun h^-1) accrete about 60 per cent of their present-day stellar mass. High mass galaxy ($M_* > 1.7e11 M_sun h^-1) assembly is dominated by accretion and merging with about 80 per cent of the stars added by the present-day. In general the simulated galaxies approximately double their mass since z=1. For massive systems this mass growth is not accompanied by significant star formation. The majority of the in-situ created stars is formed at z>2, primarily out of cold gas flows. We recover the observational result of archaeological downsizing, where the most massive galaxies harbor the oldest stars. We find that this is not in contradiction with hierarchical structure formation. Most stars in the massive galaxies are formed early on in smaller structures, the galaxies themselves are assembled late.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Gravitational Heating Helps Make Massive Galaxies Red and Dead

We study the thermal formation history of four simulated galaxies that were shown in Naab et al. (2007) to reproduce a number of observed properties of elliptical galaxies. The temperature of the gas in the galaxies is steadily increasing with decreasing redshift, although much of the gas has a cooling time shorter than the Hubble time. The gas is being heated and kept hot by gravitational heating processes through the release of potential energy from infalling stellar clumps. The energy is dissipated in supersonic collisions of infalling gas lumps with the ambient gas and through the dynamical capturing of satellite systems causing gravitational wakes that transfer energy to the surrounding gas. Furthermore dynamical friction from the infalling clumps pushes out dark matter, lowering the central dark matter density by up to a factor of two from z=3 to z=0. In galaxies in which the late formation history (z<2) is dominated by minor merging and accretion the energy released (E~5x10^{59} ergs) from gravitational feedback is sufficient to form red and dead elliptical galaxies by z~1 even in the absence of supernova and AGN feedback.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Minor mergers and the size evolution of elliptical galaxies

Using a high resolution hydrodynamical cosmological simulation of the formation of a massive spheroidal galaxy we show that elliptical galaxies can be very compact and massive at high redshift in agreement with recent observations. Accretion of stripped in-falling stellar material increases the size of the system with time and the central concentration is reduced by dynamical friction of the surviving stellar cores. In a specific case of a spheroidal galaxy with a final stellar mass of $1.5 \times 10^{11} M_{\odot}$ we find that the effective radius $r_e$ increases from $0.7 \pm 0.2 \rm kpc$ at z = 3 to $r_e = 2.4 \pm 0.4 \rm kpc$ at z = 0 with a concomitant decrease in the effective density of an order of magnitude and a decrease of the central velocity dispersion by approximately 20% over this time interval. A simple argument based on the virial theorem shows that during the accretion of weakly bound material (minor mergers) the radius can increase as the square of the mass in contrast to the usual linear rate of increase for major mergers. By undergoing minor mergers compact high redshift spheroids can evolve into present-day systems with sizes and concentrations similar to observed local ellipticals. This indicates that minor mergers may be the main driver for the late evolution of sizes and densities of early-type galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Simulations of the Microwave Sky

We create realistic, full-sky, half-arcminute resolution simulations of the microwave sky matched to the most recent astrophysical observations. The primary purpose of these simulations is to test the data reduction pipeline for the Atacama Cosmology Telescope (ACT) experiment; however, we have widened the frequency coverage beyond the ACT bands to make these simulations applicable to other microwave background experiments. Some of the novel features of these simulations are that the radio and infrared galaxy populations are correlated with the galaxy cluster populations, the CMB is lensed by the dark matter structure in the simulation via a ray-tracing code, the contribution to the thermal and kinetic Sunyaev-Zel'dovich (SZ) signals from galaxy clusters, groups, and the IGM has been included, and the gas prescription to model the SZ signals matches the most recent X-ray observations. Regarding the contamination of cluster SZ flux by radio galaxies, we find for 148 GHz (90 GHz) only 3% (4%) of halos have their SZ decrements contaminated at a level of 20% or more. We find the contamination levels higher for infrared galaxies. However, at 90 GHz, less than 20% of clusters with M_{200} > 2.5 x 10^{14} Msun and z<1.2 have their SZ decrements filled in at a level of 20% or more. At 148 GHz, less than 20% of clusters with M_{200} > 2.5 x 10^{14} Msun and z<0.8 have their SZ decrements filled in at a level of 50% or larger. Our models also suggest that a population of very high flux infrared galaxies, which are likely lensed sources, contribute most to the SZ contamination of very massive clusters at 90 and 148 GHz. These simulations are publicly available and should serve as a useful tool for microwave surveys to cross-check SZ cluster detection, power spectrum, and cross-correlation analyses.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Effects of the Ionizing Radiation Background on Galaxy Evolution

We find that the amount and nature of the assumed ionizing background can strongly affect galaxy formation and evolution. Galaxy evolution simulations typically incorporate an ultraviolet background which falls off rapidly above z=3; e.g., that of Haardt & Madau (1996). However, this decline may be too steep to fit the WMAP constraints on electron scattering optical depth or observations of intermediate redshift (z ~ 2-4) Ly-alpha forest transmission. As an alternative, we present simulations of the cosmological formation of individual galaxies with UV backgrounds that decline more slowly at high redshift: both a simple intensity rescaling and the background recently derived by Faucher-Giguere (2009), which softens the spectrum at higher redshifts. We also test an approximation of the X-ray background with a similar z-dependence. We find for the test galaxies that an increase in either the intensity or hardness of ionizing radiation generically pushes star formation towards lower redshifts: although overall star formation in the simulation boxes is reduced by 10-25%, the galaxies show a factor of ~2 increase in the fraction of stars within a 30 kpc radius that are formed after z=1. Other effects include late gas inflows enhanced up to 30 times, stellar half-mass radii decreased by up to 30%, central velocity dispersions increased up to 40%, and a strong reduction in substructure. The magnitude of the effects depends on the environmental/accretion properties of the particular galaxy.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Accurate Realizations of the Ionized Gas in Galaxy Clusters: Calibrating Feedback

Using the full, three-dimensional potential of galaxy cluster halos (drawn from an N-body simulation of the current, most favored cosmology), the distribution of the X-ray emitting gas is found by assuming a polytropic equation of state and hydrostatic equilibrium, with constraints from conservation of energy and pressure balance at the cluster boundary. The resulting properties of the gas for these simulated redshift zero clusters (the temperature distribution, mass-temperature and luminosity-temperature relations, and the gas fraction) are compared with observations in the X-ray of nearby clusters. The observed properties are reproduced only under the assumption that substantial energy injection from non-gravitational sources has occurred. Our model does not specify the source, but star formation and AGN may be capable of providing this energy, which amounts to 3 to 5 x10^{-5} of the rest mass in stars (assuming ten percent of the gas initially in the cluster forms stars). With the method described here it is possible to generate realistic X-ray and Sunyaev-Zel'dovich cluster maps and catalogs from N-body simulations, with the distributions of internal halo properties (and their trends with mass, location, and time) taken into account.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Growth of structure seeded by primordial black holes

We discuss the possibilities for primordial black holes (PBHs) to grow via the accretion of dark matter. In agreement with previous works, we find that accretion during the radiation-dominated era does not lead to a significant mass increase. However, during matter-domination, PBHs may grow by up to two orders of magnitude in mass through the acquisition of large dark matter halos. We discuss the possibility of PBHs being an important component in dark matter halos of galaxies as well as their potential to explain the ultra-luminous x-ray sources (ULXs) observed in nearby galactic disks. We point out that although PBHs are ruled out as the dominant component of dark matter, there is still a great deal of parameter space open to them playing a role in the modern-day universe. For example, a primordial halo population of PBHs each at $10^{2.5} M_\odot$ making up 0.1% of the dark matter grow to $10^{4.5} M_\odot$ via the accumulation of dark matter halos to account for $\sim 10%$ of the dark matter mass by a redshift of $z \approx 30$. These intermediate mass black holes may then ``light up'' when passing through molecular clouds, becoming visible as ULXs at the present day, or they may form the seeds for supermassive black holes at the centers of galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Testing a new analytic model for gravitational lensing probabilities

We study gravitational lensing with a multiple lens plane approach, proposing a simple analytical model for the probability distribution function (PDF) of the dark matter convergence, kappa, for the different lens planes in a given cosmology as a function of redshift and smoothing angle, theta. The model is fixed solely by the variance of kappa, which in turn is fixed by the amplitude of the power spectrum, sigma_8. We test the PDF against a high resolution Tree-Particle-Mesh simulation and find that it is far superior to the Gaussian or the lognormal, especially for small values of theta << 1 arcmin and at large values of kappa relevant to strong lensing. With this model, we predict the probabilities of strong lensing by a single plane or by multiple planes. We find that for theta ~ 10 arcsec, a single plane accounts for almost all (~ 98%) of the strong lensing cases for source redshift unity. However, for a more typical source redshift of 4, about 12% of the strong lensing cases will result from the contribution of a secondary clump of matter along the line of sight, introducing a systematic error in the determination of the surface density of clusters, typically overestimating it by about 2-5%. We also find that matter inhomogenieties introduce a dispersion in the value of the angular diameter distance about its cosmological mean. The probable error relative to the mean increases with redshift to a value of about 8% for z ~ 6 and theta ~ 10 arcsec.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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)

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

What does the local black hole mass distribution tell us about the evolution of the quasar luminosity function?

We present a robust method to derive the duty cycle of QSO activity based on the empirical QSO luminosity function and on the present-day linear relation between the masses of supermassive black holes and those of their spheroidal host stellar systems. It is found that the duty cycle is substantially less than unity, with characteristic values in the range $3-6\times 10^{-3}$. Finally, we tested the expectation that the QSO luminosity evolution and the star formation history should be roughly parallel, as a consequence of the above--mentioned relation between BH and galaxy masses.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Physical Bias of Galaxies From Large-Scale Hydrodynamic Simulations

We analyze a new large-scale ($100h^{-1}$Mpc) numerical hydrodynamic simulation of the popular $Λ$CDM cosmological model, including in our treatment dark matter, gas and star-formation, on the basis of standard physical processes. The method, applied with a numerical resolution of $<200h^{-1}$kpc (which is still quite coarse for following individual galaxies, especially in dense regions), attempts to estimate where and when galaxies form. We then compare the smoothed galaxy distribution with the smoothed mass distribution to determine the "bias" defined as $b\equiv (δM/M)_{gal}/(δM/M)_{total}$ on scales large compared with the code numerical resolution (on the basis of resolution tests given in the appendix of this paper). We find that (holding all variables constant except the quoted one) bias increases with decreasing scale, with increasing galactic age or metallicity and with increasing redshift of observations. At the $8h^{-1}$Mpc fiducial comoving scale bias (for bright regions) is 1.35 at $z=0$ reaching to 3.6 at $z=3$, both numbers being consistent with extant observations. We also find that $(10-20)h^{-1}$Mpc voids in the distribution of luminous objects are as observed (i.e., observed voids are not an argument against CDM-like models) and finally that the younger systems should show a colder Hubble flow than do the early type galaxies (a testable proposition). Surprisingly, little evolution is found in the amplitude of the smoothed galaxy-galaxy correlation function (as a function of {\it comoving} separation). Testing this prediction vs observations will allow a comparison between this work and that of Kauffmann et al which is based on a different physical modelingmethod.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Where are the Baryons?

New, high resolution, large-scale, cosmological hydrodynamic galaxy formation simulations of a standard cold dark matter model (with a cosmological constant) are utilized to predict the distribution of baryons at the present and at moderate redshift. It is found that the average temperature of baryons is an increasing function of time, with most of the baryons at the present time having a temperature in the range 10^{5-7} K. Thus, not only is the universe dominated by dark matter, but more than one half of the normal matter is yet to be detected. Detection of this warm/hot gas poses an observational challenge, requiring sensitive EUV and X-ray satellites. Signatures include a soft, cosmic X-ray background, apparent warm components in hot clusters due to both intrinsic warm intra-cluster gas and warm inter-cluster gas projected onto clusters along the line of sight, absorption lines in X-ray and UV quasar spectra [e.g., O VI (1032,1038)A lines, OVII 574 eV line], strong emission lines (e.g., O VIII 653 eV line) and low redshift, broad, low column density $\lya$ absorption lines. We estimate that approximately 1/4 of the extragalactic soft X-ray background (SXRB) (at 0.7 keV) arises from the warm/hot gas, half of it coming from $z<0.65$ and three-quarters from $z<1.00$, so the source regions should be identifiable on deep optical images.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Background X-ray Emission from Hot Gas in CDM and CDM+Lambda Universes: Spectral Signatures

We present a new treatment of two popular models for the growth of structure, examining the X-ray emission from hot gas with allowance for spectral line emission from various atomic species, primarily ``metals". The X-ray emission from the bright cluster sources is not significantly changed from prior work and, as noted earlier, shows the CDM$+Λ$ model (LCDM) to be consistent but the standard, COBE normalized model (SCDM) to be inconsistent with existing observations --- after allowance for still the considerable numerical modelling uncertainties. But we find one important new result. Radiation in the softer band 0.5-1.0keV is predominantly emitted by gas far from cluster centers (hence ``background"). This background emission dominates over the cluster emission below 1keV and observations of it should show clear spectral signatures indicating its origin. In particular the ``iron blend" should be seen prominantly in this spectral bin from cosmic background hot gas at high galactic latitudes and should show shadowing against the SMC indicating its extragalactic origin. Certain OVII lines also provide a signature of this gas which emits a spectrum characteristic of $10^{6.6\pm 0.6}$K gas. Recent ASCA observations of the X-ray background tentatively indicate the presence of component with exactly the spectral features we predict here.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Topology of Large-Scale Structure by Galaxy Type: Hydrodynamic Simulations

The topology of large scale structure is studied as a function of galaxy type using the genus statistic. In hydrodynamical cosmological CDM simulations, galaxies form on caustic surfaces (Zeldovich pancakes) then slowly drain onto filaments and clusters. The earliest forming galaxies in the simulations (defined as ``ellipticals") are thus seen at the present epoch preferentially in clusters (tending toward a meatball topology), while the latest forming galaxies (defined as ``spirals") are seen currently in a spongelike topology. The topology is measured by the genus (= number of ``donut" holes - number of isolated regions) of the smoothed density-contour surfaces. The measured genus curve for all galaxies as a function of density obeys approximately the theoretical curve expected for random-phase initial conditions, but the early forming elliptical galaxies show a shift toward a meatball topology relative to the late forming spirals. Simulations using standard biasing schemes fail to show such an effect. Large observational samples separated by galaxy type could be used to test for this effect.

0 citations0 reviewsNo code linkedOpen access
preprint1994arXiv

Testing Cosmogonic Models with Gravitational Lensing

Gravitational lensing provides a strict test of cosmogonic models. Detailed numerical propagation of light rays through a universe having a distribution of inhomogeneities derived from the standard CDM (cold dark matter) scenario with the aid of massive, fully nonlinear computer simulations is used to test the model. It predicts that we should have seen far more widely split quasar images than have been found. These and other inconsistencies rule out the COBE normalized CDM model with $Ω=1$ and $H_0=50$km/s/Mpc; but variants of this model might be constructed, which could pass the stringent tests provided by strong gravitational lensing.

0 citations0 reviewsNo code linkedOpen access