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

Renyue Cen

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

ResearcherUnclaimed source record
astro-ph.COastro-ph.GAastro-phastro-ph.HEastro-ph.IMastro-ph.SRgr-qchep-phhep-th

Catalog footprint

What is connected

90works
9topics
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

90 published item(s)

preprint2025arXiv

Rapid formation of a very massive star >50000 $M_\odot$ and subsequently an IMBH from runaway collisions. Direct N-body and Monte Carlo simulations of dense star clusters

Context. We present simulations of a massive young star cluster using \textsc{Nbody6++GPU} and \textsc{MOCCA}. The cluster is initially more compact than previously published models, with one million stars, a total mass of $5.86 \times 10^5~\mathrm{M}_{\odot}$, and a half-mass radius of $0.1~\mathrm{pc}$. Aims. We analyse the formation and growth of a very massive star (VMS) through successive stellar collisions and investigate the subsequent formation of an intermediate-mass black hole (IMBH) in the core of a dense star cluster. Methods. We use both direct \textit{N}-body and Monte Carlo simulations, incorporating updated stellar evolution prescriptions (SSE/BSE) tailored to massive stars and VMSs. These include revised treatments of stellar radii, rejuvenation, and mass loss during collisions. While the prescriptions represent reasonable extrapolations into the VMS regime, the internal structure and thermal state of VMSs formed through stellar collisions remain uncertain, and future work may require further refinement. Results. We find that runaway stellar collisions in the cluster core produce a VMS exceeding $5 \times 10^4~\mathrm{M}_{\odot}$ within 5 Myr, which subsequently collapses into an IMBH. Conclusions. Our model suggests that dense stellar environments may enable the formation of very massive stars and massive black hole seeds through runaway stellar collisions. These results provide a potential pathway for early black hole growth in star clusters and offer theoretical context for interpreting recent JWST observations of young, compact clusters at high redshift.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

Resolving the Origins and Pathways of Ionizing Radiation Escape with UV Integral Field Spectroscopy

The Epoch of Reionization marks the last major phase transition in the early Universe, during which the majority of neutral hydrogen once filling the intergalactic medium was ionized by the first galaxies. The James Webb Space Telescope (JWST) is now identifying promising galaxy candidates capable of producing sufficient ionizing photons to drive this transformation. However, the fraction of these photons that escape into intergalactic space--the escape fraction--remains highly uncertain. Stellar feedback is thought to play a critical role in carving low-density channels that allow ionizing radiation to escape, but the dominant mechanisms, their operation, and their connection to observable signatures are not well understood. Local analogs of high-redshift galaxies offer a powerful alternative for studying these processes, since ionizing radiation is unobservable at high redshift due to intergalactic absorption. However, current UV space-based instrumentation lacks the spatial resolution and sensitivity required to fully address this problem. The core challenge lies in the multiscale nature of LyC escape: ionizing photons are generated on scales of 1--100 pc in super star clusters but must traverse the circumgalactic medium which can extend beyond 100 kpc. The proposed Habitable Worlds Observatory (HWO) will provide a platform for future UV instruments capable of resolving these scales. In this article, we present a science case for understanding how LyC photons escape from star-forming galaxies and define the observational requirements for future instruments aboard HWO, including a UV integral field spectrograph (IFS).

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

AMBER: A Semi-numerical Abundance Matching Box for the Epoch of Reionization

The Abundance Matching Box for the Epoch of Reionization (AMBER) is a semi-numerical code for modeling the cosmic dawn. The new algorithm is not based on the excursion set formalism for reionization, but takes the novel approach of calculating the reionization-redshift field $z_\mathrm{re}(\boldsymbol{x})$ assuming that hydrogen gas encountering higher radiation intensity are photoionized earlier. Redshift values are assigned while matching the abundance of ionized mass according to a given mass-weighted ionization fraction $\bar{x}_\mathrm{i}(z)$. The code has the unique advantage of allowing users to directly specify the reionization history through the redshift midpoint $z_\mathrm{mid}$, duration $Δ_\mathrm{z}$, and asymmetry $A_\mathrm{z}$ input parameters. The reionization process is further controlled through the minimum halo mass $M_\mathrm{min}$ for galaxy formation and the radiation mean free path $l_\mathrm{mfp}$ for radiative transfer. We implement improved methods for constructing density, velocity, halo, and radiation fields, which are essential components for modeling reionization observables. We compare AMBER with two other semi-numerical methods and find that our code more accurately reproduces the results from radiation-hydrodynamic simulations. The parallelized code is over four orders of magnitude faster than radiative transfer simulations and will efficiently enable large-volume models, full-sky mock observations, and parameter-space studies. AMBER will be made publicly available to facilitate and transform studies of the EoR.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Detecting preheating in proto-clusters with Lyman-$α$ Forest Tomography

Studies of low redshift galaxy clusters suggest the intra-cluster medium (ICM) has experienced non-gravitational heating during the formation phase of the clusters. Using simple phenomenological heating prescriptions, we simulate the effect of this preheating of the nascent ICM in galaxy proto-clusters and examine its effect on Lyman-$α$ (Ly$α$) forest tomographic maps. We analyse a series of cosmological zoom-in simulations of proto-clusters within the framework of the Ly$α$ transmission-dark matter (DM) density distribution. We find that the more energy is injected into the proto-ICM at $z$ = 3, the more the distribution at high DM density tilts towards higher Ly$α$ transmission. This effect has been confirmed in both low-resolution simulations adopting a preheating scheme based on entropy floors, as well as in higher-resolution simulations with another scheme based on energy floors. The evolution of the slope of this distribution is shown to vary with redshift. The methodology developed here can be applied to current and upcoming Ly$α$ forest tomographic survey data to help constrain feedback models in galaxy proto-clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Osaka Feedback Model II: Modeling Supernova Feedback Based on High-Resolution Simulations

Feedback from supernovae (SNe) is an essential mechanism that self-regulates the growth of galaxies, and a better model of SN feedback is still needed in galaxy formation simulations. In the first part of this paper, using an Eulerian hydrodynamic code Athena++, we find universal scaling relations for the time evolution of momentum and radius for a superbubble, when the momentum and time are scaled by those at the shell-formation time. In the second part of this paper, we develop an SN feedback model based on the Athena++ simulation results utilizing Voronoi tessellation around each star particle, and implement it into the GADGET3-Osaka smoothed particle hydrodynamic code. Our feedback model was demonstrated to be isotropic and conservative in terms of energy and momentum. We examined the mass/energy/metal loading factors and find that our stochastic thermal feedback model produced galactic outflow that carries metals high above the galactic plane but with weak suppression of star formation. Additional mechanical feedback further suppressed star formation and brought the simulation results in better agreement with the observations of the Kennicutt--Schmidt relation, with all the results being within the uncertainties of observed data. We argue that both thermal and mechanical feedback are necessary for the SN feedback model of galaxy evolution when an individual SN bubble is unresolved.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Pilot-WINGS: An extended MUSE view of the structure of Abell 370

We investigate the strong-lensing cluster Abell 370 (A370) using a wide Integral Field Unit (IFU) spectroscopic mosaic from the Multi-Unit Spectroscopic Explorer (MUSE). IFU spectroscopy provides significant insight into the structure and mass content of galaxy clusters, yet IFU-based cluster studies focus almost exclusively on the central Einstein-radius region. Covering over 14 arcmin$^2$, the new MUSE mosaic extends significantly beyond the A370 Einstein radius, providing, for the first time, a detailed look at the cluster outskirts. Combining these data with wide-field, multi-band Hubble Space Telescope (HST) imaging from the BUFFALO project, we analyse the distribution of objects within the cluster and along the line of sight. Identifying 416 cluster galaxies, we use kinematics to trace the radial mass profile of the halo, providing a mass estimate independent from the lens model. We also measure radially-averaged properties of the cluster members, tracking their evolution as a function of infall. Thanks to the high spatial resolution of our data, we identify six cluster members acting as galaxy-galaxy lenses, which constrain localized mass distributions beyond the Einstein radius. Finally, taking advantage of MUSE's 3D capabilities, we detect and analyse multiple spatially extended overdensities outside of the cluster that influence lensing-derived halo mass estimates. We stress that much of this work is only possible thanks to the robust, extended IFU coverage, highlighting its importance even in less optically dense cluster regions. Overall, this work showcases the power of combining HST+MUSE, and serves as the initial step towards a larger and wider program targeting several clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

SCORCH. III. Analytical Models of Reionization with Varying Clumping Factors

In the Simulations and Constructions of the Reionization of Cosmic Hydrogen (SCORCH) project, we compare analytical models of the hydrogen ionization fraction with radiation-hydrodynamic simulations. We derive analytical models of the mass-weighted hydrogen ionization fraction from the local ionization balance equations as a more accurate alternative to the widely adopted model based on the volume filling factor. In particular, our model has a recombination term quadratic in the ionization fraction, which is consistent with the two-body interaction nature of recombination. Then, we use the radiation-hydrodynamic simulations to study the clumping factors needed to solve the analytical equations, and provide accurate fitting functions. We find that the ionized hydrogen clumping factors from our radiative transfer simulations are significantly different than those from other simulations that use a uniform photoionization background. In addition to redshift dependence, we also see the dependence of ionized hydrogen clumping factor on ionization fraction, and we incorporate this into our fits. We calculate the reionization histories using our analytical models and clumping factors and compare with widely adopted models, and all of our models achieve $<7\%$ difference from simulation results while the other models have $>20\%$ deviations. The Thomson optical depths from reionization calculated from our analytical models result in $<5\%$ deviation from simulations, while the previous analytical models have $>20\%$ difference in and could result in biased conclusions of the IGM reionization.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Absorption line abundances in the SMC-like galaxy UGC 5282: evidence for ISM dilution from inflows on kiloparsec scales

We present an HST Cosmic Origins Spectrograph (COS) spectrum of the QSO SDSS J095109.12+330745.8 ($z_{em} = 0.645$) whose sightline passes through the SMC-like dwarf galaxy UGC 5282 ($M_B = -16.0$, $cz = 1577$ km s$^{-1}$ ), 1.2 kpc in projection from the central HII region of the galaxy. Damped Ly$α$ (DLA) absorption is detected at the redshift of UGC 5282 with $\log$[$N$(HI) cm$^{-2}\:$] = 20.89$^{+0.12}_{-0.21}\:\:$. Analysis of the accompanying SII, PII and OI metal lines yields a neutral gas metallicity, $Z_{HI}$, of [S/H] $\simeq$ [P/H] $= -0.80\pm0.24$. The metallicity of ionized gas from the central HII region, $Z_{HII}\:$, measured from its emission lines is [O/H]$= -0.37\pm0.10$, a difference of $+0.43\pm 0.26$ from $Z_{HI}\:$. This difference $δ$ is consistent with that seen towards HII regions in other star-forming galaxies and supports the idea that ionized gas near star forming regions shows systematically higher metallicities than exist in the rest of a galaxy's neutral interstellar medium (ISM). The positive values of $δ$ found in UGC 5282 (and the other star forming galaxies) is likely due to infalling low metallicity gas from the intergalactic medium that mixes with the galaxy's ISM on kpc scales. This model is also consistent with broad Ly$α$ emission detected at the bottom of the DLA absorption, offset by $\sim 125$ km s$^{-1}$ from the absorption velocity. Models of galaxy evolution that attempt to replicate population characteristics, such as the mass-metallicity relation, may need to start with a galaxy metallicity represented by $Z_{HI}$ rather than that measured traditionally from $Z_{HII}\:$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Assembly Conformity of Structure Growth: Fossil versus Normal Groups of Galaxies

Using a semi-analytic method calibrated to the global star formation history and the stellar mass function at $z=0$, we attempt to understand the most stellar deficient galaxy groups. We argue such groups are a kind of fossil group (FGs) -- in comparison to the normal groups of galaxies, they assemble both halo and stellar mass earlier. We find there is a central galaxy and satellite conformity between these FGs and normal groups: centrals and satellites in the former form earlier and more stellar deficient than their counterparts of the latter. We term this effect "Assembly Conformity" of dark matter halos. This effect accounts for about 70 percent of the difference in stellar content between FGs and normal groups. When split by the peak redshift for the star formation rate of a group, the mass functions of satellite halos on either side of the peak redshift are found to be indistinguishable between FGs and normal groups, indicating a self-similarity of halo assembly with respect to the peak. The "baryonic environmental" effect due to ram-pressure and gas heating accounts for about 30 percent of the difference in stellar content. While the total stellar mass of FGs is lower than that of normal groups, we predict that the mass of the brightest central galaxy of FGs is, on average, higher than that of normal groups. We also predict that in the central galaxies of FGs, there is a negative stellar age gradient from the center outward, where the opposite is expected for those in normal groups.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Circumnuclear Molecular Gas in Low-redshift Quasars and Matched Star-forming Galaxies

A series of gravitational instabilities in a circumnuclear gas disk (CND) are required to trigger gas transport to a central supermassive black hole (SMBH) and ignite Active Galactic Nuclei (AGNs). A test of this scenario is to investigate whether an enhanced molecular gas mass surface density ($Σ_{\rm mol}$) is found in the CND-scale of quasars relative to a comparison sample of inactive galaxies. Here we performed sub-kpc resolution CO(2-1) observations with ALMA of four low-redshift ($z \sim 0.06$), luminous ($\sim 10^{45}$ erg s$^{-1}$) quasars with each matched to a different star-forming galaxy, having similar redshift, stellar mass, and star-formation rate. We detected CO(2-1) emission from all quasars, which show diverse morphologies. Contrary to expectations, $Σ_{\rm mol}$ of the quasar sample, computed from the CO(2-1) luminosity, tends to be smaller than the comparison sample at $r < 500$ pc; there is no systematic enhancement of $Σ_{\rm mol}$ in our quasars. We discuss four possible scenarios that would explain the lower molecular gas content (or CO(2-1) luminosity as an actual observable) at the CND-scale of quasars, i.e., AGN-driven outflows, gas-rich minor mergers, time-delay between the onsets of a starburst-phase and a quasar-phase, and X-ray-dominated region (XDR) effects on the gas chemical abundance and excitation. While not extensively discussed in the literature, XDR effects can have an impact on molecular mass measurements particularly in the vicinity of luminous quasar nuclei; therefore higher resolution molecular gas observations, which are now viable using ALMA, need to be considered.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Evidence for Infalling Gas in a Lyman-$α$ Blob

Lyman-$α$ blobs (LABs) are spatially extended nebulae of emission in the Ly$α$ line of hydrogen, seen at high redshifts$^{1,2}$, and most commonly found in the dense environment of star-forming galaxies$^{3,4}$. The origin of Ly$α$ emission in the LABs is still unclear and under debate$^{5}$. Proposed powering sources generally fall into two categories: (1) photoionization, galactic super-winds/outflows, resonant scattering of Ly$α$ photons from starbursts or active galactic nuclei (AGNs)$^{6,7,8,9,10}$ and (2) cooling radiation from cold streams of gas accreting onto galaxies$^{12}$. Here we analyze the gas kinematics within a LAB providing rare observational evidence for infalling gas. This is consistent with the release of gravitational accretion energy as cold streams radiate Ly$α$ photons. It also provides direct evidence for possible cold streams feeding the central galaxies. The infalling gas is not important by mass but hints at more than one mechanism to explain the origin of the extended Ly$α$ emission around young galaxies. It is also possible that the infalling gas may represent material falling back to the galaxy from where it originated, forming a galactic fountain.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Non-steady heating of cool cores of galaxy clusters by ubiquitous turbulence and AGN

Recent cosmological simulations have shown that turbulence should be generally prevailing in clusters because clusters are continuously growing through matter accretion. Using one-dimensional hydrodynamic simulations, we study the heating of cool-core clusters by the ubiquitous turbulence as well as feedback from the central active galactic nuclei (AGNs) for a wide range of cluster and turbulence parameters, focusing on the global stability of the core. We find that the AGN shows intermittent activities in the presence of moderate turbulence similar to the one observed with Hitomi. The cluster core maintains a quasi-equilibrium state for most of the time because the heating through turbulent diffusion is nearly balanced with radiative cooling. The balance is gradually lost because of slight dominance of the radiative cooling, and the AGN is ignited by increased gas inflow. Finally, when the AGN bursts, the core is heated almost instantaneously. Thanks to the pre-existing turbulence, the heated gas is distributed throughout the core without becoming globally unstable and causing catastrophic cooling, and the core recovers the quasi-equilibrium state. The AGN bursts can be stronger in lower-mass clusters. Predictions of our model can be easily checked with future X-ray missions like XRISM and Athena.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Physics of Non-Universal Larson's Relation

From a new perspective, we re-examine self-gravity and turbulence jointly, in hopes of understanding the physical basis for one of the most important empirical relations governing clouds in the interstellar medium (ISM), the Larson's Relation relating velocity dispersion ($σ_R$) to cloud size ($R$). We report on two key new findings. First, the correct form of the Larson's Relation is $σ_R=α_v^{1/5}σ_{pc}(R/1pc)^{3/5}$, where $α_v$ is the virial parameter of clouds and $σ_{pc}$ is the strength of the turbulence, if the turbulence has the Kolmogorov spectrum. Second, the amplitude of the Larson's Relation, $σ_{pc}$, is not universal, differing by a factor of about two between clouds on the Galactic disk and those at the Galactic center, evidenced by observational data.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Physics of Prodigious Lyman Continuum Leakers

An analysis of the dynamics of a star formation event is performed. It is shown that galaxies able to drive leftover gas to sufficient altitudes in a few million years are characterized by two basic properties: small sizes (<1kpc) and high star formation rate surface densities (Sigma_SFR > 10 Msun/yr/kpc2). For the parameter space of relevance, the outflow is primarily driven by supernovae with radiation pressure being significant but subdominant. Our analysis provides the unifying physical origin for a diverse set of observed LyC leakers, including the green-peas galaxies, [SII]-weak galaxies, Lyman-alpha emitters, with these two characteristics as the common denominator. Among verifiable physical properties of LyC leakers, we predict that (1) the newly formed stellar masses are are typically in the range of 1e8-1e10 Msun, except perhaps ULIRGs, (2) the outflow velocities are typically in the range typically of 100-600km/s, but may exceed 1e3 km/s in ULIRGs, with a strong positive correlation between the stellar masses formed and the outflow velocities, (3) the overall escape fraction of galaxies is expected to increase with increasing redshift, given the cosmological trend that galaxies become denser and more compact with increasing redshift. In addition, two interesting by-product predictions are also borne out. First, ULIRGs appear to be in a parameter region where they should be prodigious LyC leakers, unless there is a large ram-pressure. Second, Lyman break galaxies (LBGs) are not supposed to be prodigious LyC leakers in our model, given their claimed effective radii exceeding 1kpc.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The BUFFALO HST Survey

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Where do quasar hosts lie with respect to the size-mass relation of galaxies?

The evolution of the galaxy size - stellar mass (Mstellar) relation has been a puzzle for over a decade. High redshift galaxies are significantly more compact than galaxies observed today, at an equivalent mass, but how much of this apparent growth is driven by progenitor bias, minor mergers, secular processes, or feedback from AGN is unclear. To help disentangle the physical mechanisms at work by addressing the latter, we study the galaxy size - Mstellar relation of 32 carefully-selected broad-line AGN hosts at 1.2 < z < 1.7 (7.5 < log M_BH < 8.5; L_bol/L_Edd > 0.1). Using HST with multi-band photometry and state-of-the-art modeling techniques, we measure half-light radii while accounting for uncertainties from subtracting bright central point sources. We find AGN hosts to have sizes ranging from 1 to 6 kpc at Mstellar ~ 0.3 - 1 x 10^11 Msun. Thus, many hosts have intermediate sizes as compared to equal-mass star-forming and quiescent galaxies. While inconsistent with the idea that AGN feedback may induce an increase in galaxy sizes, this finding is consistent with hypotheses in which AGNs preferentially occur in systems with prior concentrated gas reservoirs, or are involved in secular compaction processes perhaps responsible for simultaneously building bulges and shutting down star formation. If driven by minor mergers, which do not grow central black holes as fast as they do bulge-like stellar structures, such a process would explain both the galaxy size - mass relation observed here and the evolution in the black hole, bulge mass relation described in a companion paper.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

C IV and He II Line Emission of Lyman Alpha Blobs: Powered by Shock Heated Gas

Utilizing {\it ab initio} ultra-high resolution hydrodynamical simulations, we investigate the properties of the interstellar and circum-galactic medium of Ly$α$ Blobs (LABs) at $z=3$, focusing on three important emission lines: Ly$α$ 1216Å, \heii 1640Å and \civ 1449Å. Their relative strengths provide a powerful probe of the thermodynamic properties of the gas when confronted with observations. By adjusting the dust attenuation effect using one parameter and matching the observed size-luminosity relation of LABs using another parameter, we show that our simulations can reproduce the observed \civ/\lya\ and \heii/\lya\ ratios adequately. This analysis provides the first successful physical model to account for simultaneously the LAB luminosity function, luminosity-size relation, and the \civ/Ly$α$ and \heii/Ly$α$ ratios, with only two parameters. The physical underpinning for this model is that, in addition to the stellar component for the \lya\ emission, the \lya\ and \civ\ emission lines due to shock heated gas are primarily collisional excitation driven and the \heii\ emission line collisional ionization driven. We find that the density, temperature and metallicity of the gas responsible for each emission line is significantly distinct, in a multi-phase interstellar and circumgalactic medium that is shock-heated primarily by supernovae and secondarily by gravitational accretion of gas.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Formation of globular clusters in atomic-cooling halos via rapid gas condensation and fragmentation during the epoch of reionization

We investigate the formation of metal-poor globular clusters (GCs) at the center of two dark matter halos with $M_{\rm halo} \sim4\times10^7 M_\odot$ at $z>10$ using cosmological radiation-hydrodynamics simulations. We find that very compact ($\lesssim$ 1 pc) and massive ($\sim6\times10^5 M_\odot$) clusters form rapidly when pristine gas collapses isothermally with the aid of efficient Ly$α$ emission during the transition from molecular-cooling halos to atomic-cooling halos. Because the local free-fall time of dense star-forming gas is very short ($\ll 1\,{\rm Myr}$), a large fraction of the collapsed gas is turned into stars before stellar feedback processes blow out the gas and shut down star formation. Although the early stage of star formation is limited to a small region of the central star-forming disk, we find that the disk quickly fragments due to metal enrichment from supernovae. Sub-clusters formed in the fragmented clouds eventually merge with the main cluster at the center. The simulated clusters closely resemble the local GCs in mass and size but show a metallicity spread that is much wider than found in the local GCs. We discuss a role of pre-enrichment by Pop III and II stars as a potential solution to the latter issue. Although not without shortcomings, it is encouraging that a naive blind (not tuned) cosmological simulation presents a possible channel for the formation of at least some massive GCs.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Reconstructing the redshift evolution of escaped ionizing flux from early galaxies with Planck and HST observations

While galaxies at $6 \lesssim z \lesssim 10$ are believed to dominate the epoch of cosmic reionization, the escape fraction of ionizing flux $f_\mathrm{esc}$ and the photon production rate $\dot n_γ$ from these galaxies must vary with redshift to simultaneously match CMB and low-redshift observations. We constrain $f_\mathrm{esc}(z)$ and $\dot n_γ(z)$ with Planck 2015 measurements of the Thomson optical depth $τ$, recent low multipole E-mode polarization measurements from Planck 2016, SDSS BAO data, and $3 \lesssim z \lesssim 10$ galaxy observations. We compare different galaxy luminosity functions that are calibrated to HST observations, using both parametric and non-parametric statistical methods that marginalize over the effective clumping factor $C_\mathrm{HII}$, the LyC production efficiency $ξ_\mathrm{ion}$, and the time-evolution of the UV limiting magnitude $dM_\mathrm{SF}/dz$. Using a power-law model, we find $f_\mathrm{esc} \lesssim 0.5$ at $z=8$ with slope $β\gtrsim 2.0$ at $68\%$ confidence with little dependence on the galaxy luminosity function or data, although there is non-negligible probability for no redshift evolution $β\sim 0$ or small escape fraction $f_\mathrm{esc} \sim 10^{-2}$. A non-parametric form for $f_\mathrm{esc}(z)$ evolves significantly with redshift, yielding $f_\mathrm{esc} \sim 0.2, 0.3, 0.6$ at $z=6,9,12$, respectively. However, a model-independent reconstruction of $\dot n_γ(z)$ predicts a suppressed escaped photon production rate at $z=9$ for the latest Planck data compared to the other models, implying a quicker period of reionization. We find evidence for redshift evolution in the limiting magnitude of the galaxy luminosity function for empirical models of the galaxy luminosity function.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Physical Nature of the Most Metal-Poor Damped Lyman Alpha Systems

Utilizing the high-resolution, large-scale LAOZI cosmological simulations we investigate the nature of the metal-poor (${\rm [Z/H]<-2}$) damped Lyman alpha systems (mpDLA) at $z=3$. The following physical picture of mpDLAs emerges. The majority of mpDLAs inhabit regions $\ge 20$~kpc from the host galaxy center on infalling cold gas streams originating from the intergalactic medium, with infall velocity of $\sim 100$ km/s and temperature of $\sim 10^{4}$ K. For each host galaxy, on average, about $1\%$ of the area within a radius $150$~kpc is covered by mpDLAs. The mpDLAs are relatively diffuse ($n_{\rm{gas}} \sim 10^{-2}$ cm$^{-3}$), Jeans quasi-stable, and have very low star formation rate ($\dotΣ \le 10^{-4} M_{\odot} \rm{\ yr}^{-1} \rm{\ kpc}^{-2}$). As mpDLAs migrate inward to the galaxy center, they mix with high metallicity gas and stellar outflows in the process, removing themselves from the metal-poor category and rendering the central ($\le 5$ kpc) regions of galaxies devoid of mpDLAs. Thus, the central regions of the host galaxies are populated by mostly metal-rich DLAs instead of mpDLAs. All observables of the simulated mpDLAs are in excellent agreement with observations, except the gas density, which is about a factor of ten lower than the value inferred observationally. However, the observationally inferred value is based on simplified assumptions that are not borne out in the simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Upper Limit on Star Formation and Metal Enrichment in Minihalos

An analysis of negative radiative feedback from resident stars in minihalos is performed. It is found that the most effective mechanism to suppress star formation is provided by infrared photons from resident stars via photo-detachment of ${\rm H^-}$. It is shown that a stringent upper bound on (total stellar mass, metallicity) of ($\sim 1000{\rm M_\odot}$, $-3.3\pm 0.2$) in any newly minted atomic cooling halo can be placed, with the actual values possibly significantly lower. This has both important physical ramifications on formation of stars and supermassive black seeds in atomic cooling halos at high redshift, pertaining to processes of low temperature metal cooling, dust formation and fragmentation, and direct consequences on the faint end galaxy luminosity function at high redshift and cosmological reionization. The luminosity function of galaxies at the epoch of reionization may be substantially affected due to the combined effect of a diminished role of minihalos and an enhanced contribution from Pop III stars in atomic cooling halos. Upcoming results on reionization optical depth from Planck High-Frequency Instrument data may provide a significant constraint on and a unique probe of this star formation physical process in minihalos. As a numerical example, in the absence of significant contributions from minihalos with virial masses below $1.5\times 10^{8}{\rm M_\odot}$ the reionization optical depth is expected to be no greater than $0.065$, whereas allowing for minihalos of masses as low as ($10^7{\rm M_\odot}$, $10^{6.5}{\rm M_\odot}$) to form stars unconstrained by this self-regulation physical process, the reionization optical depth is expected to exceed $(0.075,0.085)$, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

UV Absorption Line Ratios in Circumgalactic Medium at Low Redshift in Realistic Cosmological Hydrodynamic Simulations

Utilizing high-resolution cosmological hydrodynamic simulations we investigate various ultra-violet absorption lines in the circumgalactic medium of star forming galaxies at low redshift, in hopes of checking and alleviating the claimed observational conundrum of the ratio of NV to OVI absorbers, among others. We find a satisfactory agreement between simulations and extant observational data with respect to the ratios of the following four line pairs examined, NV/OVI, SiIV/OVI, NIII/OVI and NII/OVI. For the pairs involving nitrogen lines, we examine two cases of nitrogen abundance, one with constant N/O ratio and the other with varying N/O ratio, with the latter motivated by theoretical considerations of two different synthetic sources of nitrogen that is empirically verified independently. Along a separate vector, for all line pairs, we examine two cases of radiation field, one with the Haardt-Madau background radiation field and the other with an additional local radiation field sourced by hot gas in the host galaxy. In all cases, two-sample Kolmogorov-Smirnov tests indicate excellent agreements. We find that the apparent agreements between simulations and observations will be strongly tested, if the bulk of current upper limits of various line ratios are turned into actual detections. We show that an increase in observational sensitivity by 0.2 dex will already start to significantly constrain the models.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Coevolution Between Supermassive Black Holes and Bulges Is Not Via Internal Feedback Regulation But By Rationed Gas Supply Due To Angular Momentum Distribution

We reason that, without physical fine-tuning, neither the supermassive black holes (SMBHs) nor the stellar bulges can self-regulate or inter-regulate by driving away already fallen cold gas to produce the observed correlation between them. We suggest an alternative scenario where the observed mass ratios of the SMBHs to bulges reflect the angular momentum distribution of infallen gas such that the mass reaching the stable accretion disc is a small fraction of that reaching the bulge region, averaged over the cosmological time scales. We test this scenario using high resolution, large-scale cosmological hydrodynamic simulations (without AGN feedback), assuming the angular momentum distribution of gas landing in the bulge region to yield a Mestel disc that is supported by independent simulations resolving the Bondi radii of SMBHs. A mass ratio of $0.1-0.3\%$ between the very low angular momentum gas that free-falls to the sub-parsec region to accrete to the SMBH and the overall star formation rate is found. This ratio is found to increase with increasing redshift to within a factor of $\sim 2$, suggesting that the SMBH to bulge ratio is nearly redshift independent, with a modest increase with redshift, a testable prediction. Furthermore, the duty cycle of active galactic nuclei (AGN) with high Eddington ratios is expected to increase significantly with redshift. Finally, while SMBHs and bulges are found to coevolve on $\sim 30-150$Myr time scales or longer, there is indication that, on shorer time scales, the SMBH accretion rate and star formation may be less correlated.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Don't Forget the Forest for the Trees: The Stellar-Mass Halo-Mass Relation in Different Environments

The connection between dark matter halos and galactic baryons is often not well-constrained nor well-resolved in cosmological hydrodynamical simulations. Thus, Halo Occupation Distribution (HOD) models that assign galaxies to halos based on halo mass are frequently used to interpret clustering observations, even though it is well-known that the assembly history of dark matter halos is related to their clustering. In this paper we use high-resolution hydrodynamical cosmological simulations to compare the halo and stellar mass growth of galaxies in a large-scale overdensity to those in a large-scale underdensity (on scales of about 20 Mpc). The simulation reproduces assembly bias, that halos have earlier formation times in overdense environments than in underdense regions. We find that the stellar mass to halo mass ratio is larger in overdense regions in central galaxies residing in halos with masses between 10$^{11}$-10$^{12.9}$ M$_{\odot}$. When we force the local density (within 2 Mpc) at z=0 to be the same for galaxies in the large-scale over- and underdensities, we find the same results. We posit that this difference can be explained by a combination of earlier formation times, more interactions at early times with neighbors, and more filaments feeding galaxies in overdense regions. This result puts the standard practice of assigning stellar mass to halos based only on their mass, rather than considering their larger environment, into question.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Inflow Generated X-ray Corona Around Supermassive Black Holes and Unified Model for X-ray Emission

Three-dimensional hydrodynamic simulations, covering the spatial domain from hundreds of Schwarzschild radii to $2\ \mathrm{pc}$ around the central supermassive black hole of mass $10^8 M_\odot$, with detailed radiative cooling processes, are performed. Generically found is the existence of a significant amount of shock heated, high temperature ($\geq 10^8\ \mathrm{K}$) coronal gas in the inner ($\leq 10^4 r_\mathrm{sch}$) region. It is shown that the composite bremsstrahlung emission spectrum due to coronal gas of various temperatures are in reasonable agreement with the overall ensemble spectrum of AGNs and hard X-ray background. Taking into account inverse Compton processes, in the context of the simulation-produced coronal gas, our model can readily account for the wide variety of AGN spectral shape, which can now be understood physically. The distinguishing feature of our model is that X-ray coronal gas is, for the first time, an integral part of the inflow gas and its observable characteristics are physically coupled to the concomitant inflow gas. One natural prediction of our model is the anti-correlation between accretion disk luminosity and spectral hardness: as the luminosity of SMBH accretion disk decreases, the hard X-ray luminosity increases relative to the UV/optical luminosity.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the Diffuse Lyman-alpha Halo Around Lyman-alpha Emitting Galaxies

Ly$α$ photons scattered by neutral hydrogen atoms in the circumgalactic media or produced in the halos of star-forming galaxies are expected to lead to extended Ly$α$ emission around galaxies. Such low surface brightness Ly$α$ halos (LAHs) have been detected by stacking Ly$α$ images of high-redshift star-forming galaxies. We study the origin of LAHs by performing radiative transfer modeling of nine $z=3.1$ Lyman-Alpha Emitters (LAEs) in a high resolution hydrodynamic cosmological galaxy formation simulation. We develop a method of computing the mean Ly$α$ surface brightness profile of each LAE by effectively integrating over many different observing directions. Without adjusting any parameters, our model yields an average Ly$α$ surface brightness profile in remarkable agreement with observations. We find that observed LAHs cannot be accounted for solely by photons originating from the central LAE and scattered to large radii by hydrogen atoms in the circumgalactic gas. Instead, Ly$α$ emission from regions in the outer halo is primarily responsible for producing the extended LAHs seen in observations, which potentially includes both star-forming and cooling radiation. With the limit on the star formation contribution set by the ultra-violet (UV) halo measurement, we find that cooling radiation can play an important role in forming the extended LAHs. We discuss the implications and caveats of such a picture.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Quantifying Distributions of Lyman Continuum Escape Fraction

Simulations have indicated that most of the escaped Lyman continuum photons escape through a minority of solid angles with near complete transparency, with the remaining majority of the solid angles largely opaque, resulting in a very broad and skewed probability distribution function (PDF) of the escape fraction when viewed at different angles. Thus, the escape fraction of Lyman continuum photons of a galaxy observed along a line of sight merely represents the properties of the interstellar medium along that line of sight, which may be an ill-representation of true escape fraction of the galaxy averaged over its full sky. Here we study how Lyman continuum photons escape from galaxies at $z=4-6$, utilizing high-resolution large-scale cosmological radiation-hydrodynamic simulations. We compute the PDF of the mean escape fraction ($\left<f_{\rm esc,1D}\right>$) averaged over mock observational samples, as a function of the sample size, compared to the true mean (had you an infinite sample size). We find that, when the sample size is small, the apparent mean skews to the low end. For example, for a true mean of 6.7%, an observational sample of (2,10,50) galaxies at $z=4$ would have have 2.5% probability of obtaining the sample mean lower than $\left<f_{\rm esc,1D}\right>=$(0.007%, 1.8%, 4.1%) and 2.5% probability of obtaining the sample mean being greater than (43%, 18%, 11%). Our simulations suggest that at least $\sim$ 100 galaxies should be stacked in order to constrain the true escape fraction within 20% uncertainty.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

SCORCH I: The Galaxy-Halo Connection in the First Billion Years

SCORCH (Simulations and Constructions of the Reionization of Cosmic Hydrogen) is a new project to study the Epoch of Reionization (EoR). In this first paper, we probe the connection between observed high-redshift galaxies and simulated dark matter halos to better understand the abundance and evolution of the primary source of ionizing radiation. High-resolution N-body simulations are run to quantify the abundance of dark matter halos as a function of mass $M$, accretion rate $\dot{M}$, and redshift $z$. A new fit for the halo mass function $dn/dM$ is $\approx 20\%$ more accurate at the high-mass end where bright galaxies are expected to reside. A novel approach is used to fit the halo accretion rate function $dn/d\dot{M}$ in terms of the halo mass function. Abundance matching against the observed galaxy luminosity function is used to estimate the luminosity-mass relation and the luminosity-accretion-rate relation. The inferred star formation efficiency is not monotonic with $M$ nor $\dot{M}$, but reaches a maximum value at a characteristic mass $\sim 2 \times 10^{11}\ M_\odot$ and a characteristic accretion rate $\sim 6 \times 10^2\ M_\odot/{\rm yr}$ at $z \approx 6$. We find a universal EoR luminosity-accretion-rate relation and construct a fiducial model for the galaxy luminosity function. The Schechter parameters evolve such that $ϕ_\star$ decreases, $M_\star$ is more positive (fainter), and $α$ is more negative (steeper) at higher redshifts. We forecast for the upcoming James Webb Space Telescope and show that with apparent magnitude limit $m_{\rm AB} \approx 31\ (32)$, it can observe $\gtrsim 11\ (24)$ unlensed galaxies per square degree per unit redshift at least down to $M_\star$ at $z \lesssim 13\ (14)$.

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

Testing Dark Matter Halo Models of Quasars With Thermal Sunyaev-Zeldovich Effect

A statistical analysis of stacked Compton$-y$ maps of quasar hosts with a median redshift of $1.5$ using Millennium Simulation is performed to address two issues, one on the feedback energy from quasars and the other on testing dark matter halo models for quasar hosts. On the first, we find that, at the resolution of FWHM=$10$ arcmin obtained by Planck data, the observed thermal Sunyaev-Zeldovich (tSZ) effect can be entirely accounted for and explained by the thermal energy of halos sourced by gravitational collapse of halos, without a need to invoke additional, large energy sources, such as quasar or stellar feedback. Allowing for uncertainties of dust temperature in the calibration of observed Comton$-y$ maps, the maximum additional feedback energy is $\sim 25\%$ of that previously suggested. Second, we show that, with FWHM=$1$ arcmin beam, tSZ measurements will provide a potentially powerful test of quasar-hosting dark matter halo models, limited only by possible observational systematic uncertainties, not by statistical ones, even in the presence of possible quasar feedback.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Non-Linear Power Spectrum of the Lyman Alpha Forest

The Lyman alpha forest power spectrum has been measured on large scales by the BOSS survey in SDSS-III at $z\sim 2.3$, has been shown to agree well with linear theory predictions, and has provided the first measurement of Baryon Acoustic Oscillations at this redshift. However, the power at small scales, affected by non-linearities, has not been well examined so far. We present results from a variety of hydrodynamic simulations to predict the redshift space non-linear power spectrum of the Lyman Alpha transmission for several models, testing the dependence on resolution and box size. A new fitting formula is introduced to facilitate the comparison of our simulation results with observations and other simulations. The non-linear power spectrum has a generic shape determined by a transition scale from linear to non-linear anisotropy, and a Jeans scale below which the power drops rapidly. In addition, we predict the two linear bias factors of the Lyman Alpha forest and provide a better physical interpretation of their values and redshift evolution. The dependence of these bias factors and the non-linear power on the amplitude and slope of the primordial fluctuations power spectrum, the temperature-density relation of the intergalactic medium, and the mean Lyman Alpha transmission, as well as the redshift evolution, is investigated and discussed in detail. A preliminary comparison to the observations shows that the predicted redshift distortion parameter is in good agreement with the recent determination of Blomqvist et al., but the density bias factor is lower than observed. We make all our results publicly available in the form of tables of the non-linear power spectrum that is directly obtained from all our simulations, and parameters of our fitting formula.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Topology of neutral hydrogen distribution with the Square Kilometer Array

Morphology of the complex HI gas distribution can be quantified by statistics like the Minkowski functionals, and can provide a way to statistically study the large scale structure in the HI maps both at low redshifts, and during the epoch of reionization (EoR). At low redshifts, the 21cm emission traces the underlying matter distribution. Topology of the HI gas distribution, as measured by the genus, could be used as a "standard ruler". This enables the determination of distance-redshift relation and also the discrimination of various models of dark energy and of modified gravity. The topological analysis is also sensitive to certain primordial non-Gaussian features. Compared with two-point statistics, the topological statistics are more robust against the nonlinear gravitational evolution, bias, and redshift-space distortion. The HI intensity map observation naturally avoids the sparse sampling distortion, which is an important systematic in optical galaxy survey. The large cosmic volume accessible to SKA would provide unprecedented accuracy using such a measurement... [abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Towards simulating star formation in turbulent high-z galaxies with mechanical supernova feedback

Feedback from supernovae is essential to understanding the self-regulation of star formation in galaxies. However, the efficacy of the process in a cosmological context remains unclear due to excessive radiative losses during the shock propagation. To better understand the impact of SN explosions on the evolution of galaxies, we perform a suite of high-resolution (12 pc), zoom-in cosmological simulations of a Milky Way-like galaxy at z=3 with adaptive mesh refinement. We find that SN explosions can efficiently regulate star formation, leading to the stellar mass and metallicity consistent with the observed mass-metallicity relation and stellar mass-halo mass relation at z~3. This is achieved by making three important changes to the classical feedback scheme: i) the different phases of SN blast waves are modelled directly by injecting radial momentum expected at each stage, ii) the realistic time delay of SNe, commencing at as early as 3 Myr, is required to disperse very dense gas before a runaway collapse sets in at the galaxy centre via mergers of gas clumps, and iii) a non-uniform density distribution of the ISM is taken into account below the computational grid scale for the cell in which SN explodes. The last condition is motivated by the fact that our simulations still do not resolve the detailed structure of a turbulent ISM in which the fast outflows can propagate along low-density channels. The simulated galaxy with the SN feedback model shows strong outflows, which carry approximately ten times larger mass than star formation rate, as well as smoothly rising circular velocity. Other feedback models that do not meet the three conditions form too many stars, producing a peaked rotation curve. Our results suggest that understanding the structure of the turbulent ISM may be crucial to assess the role of SN and other feedback processes in galaxy formation theory.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A New Model for Dark Matter Halos Hosting Quasars

A new model for quasar-hosting dark matter halos, meeting two physical conditions, is put forth. First, significant interactions are taken into consideration to trigger quasar activities. Second, satellites in very massive halos at low redshift are removed from consideration, due to their deficiency of cold gas. We analyze the {\em Millennium Simulation} to find halos that meet these two conditions and simultaneously match two-point auto-correlation functions of quasars and cross-correlation functions between quasars and galaxies at $z=0.5-3.2$. %The found halos have some distinct properties worth noting. The masses of found quasar hosts decrease with decreasing redshift, with the mass thresholds being $[(2-5)\times 10^{12}, (2-5)\times 10^{11}, (1-3)\times 10^{11}]\msun$ for median luminosities of $\sim[10^{46}, 10^{46}, 10^{45}]$erg/s at $z=(3.2, 1.4, 0.53)$, respectively, an order of magnitude lower than those inferred based on halo occupation distribution modeling. In this model quasar hosts are primarily massive central halos at $z\ge 2-3$ but increasingly dominated by lower mass satellite halos experiencing major interactions towards lower redshift. But below $z=1$ satellite halos in groups more massive than $\sim 2\times 10^{13}\msun$ do not host quasars. Whether for central or satellite halos, imposing the condition of significant interactions substantially boosts the clustering strength compared to the total population with the same mass cut. The inferred lifetimes of quasars at $z=0.5-3.2$ of $3-30$Myr are in agreement with observations. Quasars at $z\sim 2$ would be hosted by halos of mass $\sim 5\times 10^{11}\msun$ in this model, compared to $\sim 3\times 10^{12}\msun$ previously thought, which would help reconcile with the observed, otherwise puzzling high covering fractions for Lyman limit systems around quasars.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Diverse Properties of Interstellar Medium Embedding Gamma-Ray Bursts at the Epoch of Reionization

Analysis is performed on ultra-high resolution large-scale cosmological radiation-hydrodynamic simulations to, for the first time, quantify the physical environment of long-duration gamma-ray bursts (GRBs) at the epoch of reionization. We find that, on parsec scales, 13% of GRBs remain in high density ($\ge 10^4$cm$^{-3}$) low-temperature star-forming regions, whereas 87% of GRBs occur in low-density ($\sim 10^{-2.5}$cm$^{-3}$) high temperature regions heated by supernovae. More importantly, the spectral properties of GRB afterglows, such as the neutral hydrogen column density, total hydrogen column density, dust column density, gas temperature and metallicity of intervening absorbers, vary strongly from sightline to sightline. Although our model explains extant limited observationally inferred values with respect to circumburst density, metallicity, column density and dust properties, a substantially larger sample of high-z GRB afterglows would be required to facilitate a statistically solid test of the model. Our findings indicate that any attempt to infer the physical properties (such as metallicity) of the interstellar medium of the host galaxy based on a very small number of (usually one) sightlines would be precarious. Utilizing high-z GRBs to probe interstellar medium and intergalactic medium should be undertaken properly taking into consideration the physical diversities of the interstellar medium.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Escape fraction of ionizing photons during reionization: effects due to supernova feedback and runaway OB stars

The fraction of hydrogen ionizing photons escaping from galaxies into the intergalactic medium is a critical ingredient in the theory of reionization. We use two zoomed-in, high-resolution (4 pc), cosmological radiation hydrodynamic simulations with adaptive mesh refinement to investigate the impact of two physical mechanisms (supernova feedback and runaway OB stars) on the escape fraction (f_esc) at the epoch of reionization (z>7). We implement a new, physically motivated supernova feedback model that can approximate the Sedov solutions at all (from the free expansion to snowplow) stages. We find that there is a significant time delay of about ten million years between the peak of star formation and that of escape fraction, due to the time required for the build-up and subsequent destruction of the star-forming cloud by supernova feedback. Consequently, the photon number-weighted mean escape fraction for dwarf galaxies in halos of mass 10^8-10^10.5 Msun is found to be <fesc>~11%, although instantaneous values of f_esc>20% are common when star formation is strongly modulated by the supernova explosions. We find that the inclusion of runaway OB stars increases the mean escape fraction by 22% to <fesc>~14%. As supernovae resulting from runaway OB stars tend to occur in less dense environments, the feedback effect is enhanced and star formation is further suppressed in halos with Mvir>10^9 Msun in the simulation with runaway OB stars compared with the model without them. While both our models produce enough ionizing photons to maintain a fully ionized universe at z>7 as observed, a still higher amount of ionizing photons at z>9 appears necessary to accommodate the high observed electron optical depth inferred from cosmic microwave background observations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Evolution of Cold Streams and Emergence of the Hubble Sequence

A new physical framework for the emergence of the Hubble sequence is outlined, based on novel analyses performed to quantify the evolution of cold streams of a large sample of galaxies from a state-of-the-art ultra-high resolution, large-scale adaptive mesh-refinement hydrodynamic simulation in a fully cosmological setting. It is found that the following three key physical variables of galactic cold inflows crossing the virial sphere substantially decrease with decreasing redshift: the number of streams N_{90} that make up 90% of concurrent inflow mass flux, average inflow rate per stream dot M_{90} and mean (mass flux weighted) gas density in the streams n_{gas}. Another key variable, the stream dimensionless angular momentum parameter lambda, instead is found to increase with decreasing redshift. Assimilating these trends and others leads naturally to a physically coherent scenario for the emergence of the Hubble sequence, including the following expectations: (1) the predominance of a mixture of disproportionately small irregular and complex disk galaxies at z>2 when most galaxies have multiple concurrent streams, (2) the beginning of the appearance of flocculent spirals at z~1-2 when the number of concurrent streams are about 2-3, (3) the grand-design spiral galaxies appear at z<1 when galaxies with only one major cold stream significantly emerge. These expected general trends are in good accord with observations. Early type galaxies are those that have entered a perennial state of zero cold gas stream, with their abundance increasing with decreasing redshift.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Frequent Spin Reorientation of Galaxies due to Local Interactions

We study the evolution of angular momenta of ($M_*=10^{10}-10^{12}\msun$) galaxies utilizing large-scale ultra-high resolution cosmological hydrodynamic simulations and find that spin of the stellar component changes direction frequently, caused by major mergers, minor mergers, significant gas inflows and torques by nearby systems. The rate and nature of change of spin direction can not be accounted for by large-scale tidal torques, because the latter fall short in rates by orders of magnitude and because the apparent random swings of the spin direction are inconsistent with alignment by linear density field. The implications for galaxy formation as well as intrinsic alignment of galaxies are profound. Assuming the large-scale tidal field is the sole alignment agent, a new picture emerging is that intrinsic alignment of galaxies would be a balance between slow large-scale coherent torquing and fast spin reorientation by local interactions. What is still open is whether other processes, such as feeding galaxies with gas and stars along filaments or sheets, introduce coherence for spin directions of galaxies along the respective structures.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Gas Loss in Simulated Galaxies as They Fall into Clusters

We use high-resolution cosmological hydrodynamic galaxy formation simulations to gain insights on how galaxies lose their cold gas at low redshift as they migrate from the field to the high density regions of clusters of galaxies. We find that beyond three cluster virial radii, the fraction of gas-rich galaxies is constant, representing the field. Within three cluster-centric radii, the fraction of gas-rich galaxies declines steadily with decreasing radius, reaching $\mathbf{<10%}$ near the cluster center. Our results suggest that galaxies start to feel the impact of the cluster environment on their gas content well beyond the cluster virial radius. We show that almost all gas-rich galaxies at the cluster virial radius are falling in for the first time at nearly radial orbits. Furthermore, we find that almost no galaxy moving outward at the cluster virial radius is gas-rich (with gas to baryon ratio greater than $\mathbf{1%}$). These results suggest that galaxies that fall into clusters lose their cold gas within a single radial round-trip.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Gaussian Random Field: Physical Origin of Sersic Profiles

While the Sersic profile family provide adequate fits for the surface brightness profiles of observed galaxies, the physical origin is unknown. We show that, if the cosmological density field are seeded by random gaussian fluctuations, as in the standard cold dark matter model, galaxies with steep central profiles have simultaneously extended envelopes of shallow profiles in the outskirts, whereas galaxies with shallow central profiles are accompanied by steep density profiles in the outskirts. These properties are in accord with those of the Sersic profile family. Moreover, galaxies with steep central profiles form their central regions in smaller denser subunits that possibly merge subsequently, which naturally leads to formation of bulges. In contrast, galaxies with shallow central profiles form their central regions in a coherent fashion without significant substructure, a necessary condition for disk galaxy formation. Thus, the scenario is self-consistent with respect to the correlation between observed galaxy morphology and Sersic index. We predict further that clusters of galaxies should display a similar trend, which should be verifiable observationally.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the Origin of the Hubble Sequence: I. Insights on Galaxy Color Migration from Cosmological Simulations

An analysis of more than 3000 galaxies resolved at better than 114 pc/h at z=0.62 in a LAOZI cosmological adaptive mesh refinement hydrodynamic simulation is performed and insights gained on star formation quenching and color migration. The vast majority of red galaxies are found to be within three virial radii of a larger galaxy, at the onset of quenching when the specific star formation rate experiences the sharpest decline to fall below ~10^{-2}-10^{-1}/Gyr (depending on the redshift). We shall thus call this mechanism "environment quenching", which encompasses satellite quenching. Two physical processes are largely responsible: ram-pressure stripping first disconnects the galaxy from the cold gas supply on large scales, followed by a longer period of cold gas starvation taking place in high velocity dispersion environment, during the early part of which the existing dense cold gas in the central region (<10kpc) is consumed by in situ star formation. Quenching is found to be more efficient, but not faster, on average, in denser environment. Throughout this quenching period and the ensuing one in the red sequence galaxies follow nearly vertical tracks in the color-stellar-mass diagram. In contrast, individual galaxies of all masses grow most of their stellar masses in the blue cloud, prior to the onset of quenching, and progressively more massive blue galaxies with already relatively older mean stellar ages continue to enter the red sequence. Consequently, correlations among observables of red galaxies - such as the age-mass relation - are largely inherited from their blue progenitors at the onset of quenching. While the color makeup of the entire galaxy population strongly depends on environment, which is a direct result of environment quenching, physical properties of blue galaxies as a sub-population show little dependence on environment.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the Reversal of SFR-Density Relation at z=1: Insights from Simulations

Recent large surveys have found a reversal of the star formation rate (SFR)-density relation at z=1 from that at z=0 (e.g. Elbaz et al.; Cooper et al.), while the sign of the slope of the color-density relation remains unchanged (e.g. Cucciati et al.; Quadri et al.). We use state-of-the-art adaptive mesh refinement cosmological hydrodynamic simulations of a 21x24x20 (Mpc/h)$^3$ region centered on a cluster to examine the SFR-density and color-density relations of galaxies at z=0 and z=1. The local environmental density is defined by the dark matter mass in spheres of radius 1 Mpc/h, and we probe two decades of environmental densities. Our simulations produce a large increase of SFR with density at z=1, as in the observations of Elbaz et al. We also find a significant evolution to z=0, where the SFR-density relation is much flatter. The color-density relation in our simulations is consistent from z=1 to z=0, in agreement with observations. We find that the increase in the median SFR with local density at z=1 is due to a growing population of star-forming galaxies in higher-density environments. At z=0 and z=1 both the SFR and cold gas mass are tightly correlated with the galaxy halo mass, and therefore the correlation between median halo mass and local density is an important cause of the SFR-density relation at both redshifts. We also show that the local density on 1 Mpc/h scales affects galaxy SFRs as much as halo mass at z=0. Finally, we find indications that the role of the 1 Mpc/h scale environment reverses from z=0 to z=1: at z=0 high-density environments depress galaxy SFRs, while at z=1 high-density environments tend to increase SFRs.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Penrose At Work

If a supermassive black hole has some material orbiting around it at close to its innermost stable circular orbit (ISCO), then, when it plunges into a second supermassive black hole, the orbiting material has a velocity dispersion of order of speed of light about the orbital velocity of its host black hole. It becomes plausible that some of the orbiting material will be "catapulted" to the negative-energy ergosphere orbits of the second black hole at the plunge. This may provide an astrophysically plausible way to extract energy from the black hole, originally suggested by Penrose.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Reionization on Large Scales IV: Predictions for the 21 cm signal incorporating the light cone effect

We present predictions for the 21 cm brightness temperature power spectrum during the Epoch of Reionization (EoR). We discuss the implications of the "light cone" effect, which incorporates evolution of the neutral hydrogen fraction and 21 cm brightness temperature along the line of sight. Using a novel method calibrated against radiation-hydrodynamic simulations, we model the neutral hydrogen density field and 21 cm signal in large volumes ($L = 2$ Gpc/$h$). The inclusion of the light cone effect leads to a relative decrease of about 50% in the 21 cm power spectrum on all scales. We also find that the effect is more prominent at the midpoint of reionization and later. The light cone effect also can introduce an anisotropy along the line of sight. By decomposing the 3D power spectrum into components perpendicular to and along the line of sight, we find that in our fiducial reionization model, there is no significant anisotropy. However, parallel modes can contribute up to 40% more power for shorter reionization scenarios. The scales on which the light cone effect is relevant are comparable to scales where one measures the baryon acoustic oscillation. We argue that due to its large comoving scale and introduction of anisotropy, the light cone effect is important when considering redshift space distortions and future application to the Alcock-Paczynski test for the determination of cosmological parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Shock Waves and Cosmic Ray Acceleration in the Outskirts of Galaxy Clusters

The outskirts of galaxy clusters are continuously disturbed by mergers and gas infall along filaments, which in turn induce turbulent flow motions and shock waves. We examine the properties of shocks that form within $r_{200}$ in sample galaxy clusters from structure formation simulations. While most of these shocks are weak and inefficient accelerators of cosmic rays (CRs), there are a number of strong, energetic shocks which can produce large amounts of CR protons via diffusive shock acceleration. We show that the energetic shocks reside mostly in the outskirts and a substantial fraction of them are induced by infall of the warm-hot intergalactic medium from filaments. As a result, the radial profile of the CR pressure in the intracluster medium is expected to be broad, dropping off more slowly than that of the gas pressure, and might be even temporarily inverted, peaking in the outskirts. The volume-integrated momentum spectrum of CR protons inside $r_{200}$ has the power-law slope of $4.25 - 4.5$, indicating that the average Mach number of the shocks of main CR production is in the range of $\left< M_s \right>_{\rm CR} \approx 3 - 4$. We suggest that some radio relics with relatively flat radio spectrum could be explained by primary electrons accelerated by energetic infall shocks with $M_s >~ 3$ induced in the cluster outskirts.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Temporal Self-Organization in Galaxy Formation

We report on the discovery of a relation between the number of star formation (SF) peaks per unit time, $ν_{\rm peak}$, and the size of the temporal smoothing window function, $Δt$, used to define the peaks: $ν_{\rm peak}\proptoΔt^{1-ϕ}$ ($ϕ\sim 1.618$). This relation holds over the range of $Δt=10$ to $1000$Myr that can be reliably computed, using a large sample of galaxies obtained from a state-of-the-art cosmological hydrodynamic simulation. This means that the temporal distribution of SF peaks in galaxies as a population is fractal with a Hausdorff fractal dimension equal to $ϕ-1$. This finding reveals, for the first time, that the superficially chaotic process of galaxy formation is underlined by a temporal self-organization up to at least one gigayear. It is tempting to suggest that, given the known existence of spatial fractals (such as the power-law two-point function of galaxies), there is a joint spatio-temporal self-organization in galaxy formation. From an observational perspective, it will be urgent to devise diagnostics to probe SF histories of galaxies with good temporal resolution to facilitate a test of this prediction. If confirmed, it would provide unambiguous evidence for a new picture of galaxy formation that is interaction driven, cooperative and coherent in and between time and space. Unravelling its origin may hold the key to understanding galaxy formation.

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

Composition of Low Redshift Halo Gas

Halo gas in low-z (z<0.5) >0.1L* galaxies in high-resolution, large-scale cosmological hydrodynamic simulations is examined with respect to three components: (cold, warm, hot) with temperatures equal to (<10^5, 10^{5-6}, >10^6)K, respectively. The warm component is compared, utilizing O VI λλ1032, 1038 absorption lines, to observations and agreement is found with respect to the galaxy-O VI line correlation, the ratio of O VI line incidence rate in blue to red galaxies and the amount of O VI mass in star-forming galaxies. A detailed account of the sources of warm halo gas (stellar feedback heating, gravitational shock heating and accretion from the intergalactic medium), inflowing and outflowing warm halo gas metallicity disparities and their dependencies on galaxy types and environment is also presented. Having the warm component securely anchored, our simulations make the following additional predictions. First, cold gas is the primary component in inner regions, with its mass comprising 50% of all gas within galacto-centric radius r=(30,150)kpc in (red, blue) galaxies. Second, at r>(30,200)kpc in (red, blue) galaxies the hot component becomes the majority. Third, the warm component is a perpetual minority, with its contribution peaking at ~30% at r=100-300kpc in blue galaxies and never exceeding 5% in red galaxies. The significant amount of cold gas in low-z early-type galaxies found in simulations, in agreement with recent observations (Thom et al.), is intriguing, so is the dominance of hot gas at large radii in blue galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Enzo: An Adaptive Mesh Refinement Code for Astrophysics

This paper describes the open-source code Enzo, which uses block-structured adaptive mesh refinement to provide high spatial and temporal resolution for modeling astrophysical fluid flows. The code is Cartesian, can be run in 1, 2, and 3 dimensions, and supports a wide variety of physics including hydrodynamics, ideal and non-ideal magnetohydrodynamics, N-body dynamics (and, more broadly, self-gravity of fluids and particles), primordial gas chemistry, optically-thin radiative cooling of primordial and metal-enriched plasmas (as well as some optically-thick cooling models), radiation transport, cosmological expansion, and models for star formation and feedback in a cosmological context. In addition to explaining the algorithms implemented, we present solutions for a wide range of test problems, demonstrate the code's parallel performance, and discuss the Enzo collaboration's code development methodology.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Heavy dust obscuration of z=7 galaxies in a cosmological hydrodynamic simulation

Hubble Space Telescope observations with Wide Field Camera 3/IR reveal that galaxies at z~7 have very blue ultraviolet (UV) colors, consistent with these systems being dominated by young stellar populations with moderate or little attenuation by dust. We investigate UV and optical properties of the high-z galaxies in the standard cold dark matter model using a high-resolution adaptive mesh refinement cosmological hydrodynamic simulation. For this purpose, we perform panchromatic three-dimensional dust radiative transfer calculations on 198 galaxies of stellar mass 5x10^8-3x10^{10} Msun with three parameters, the dust-to-metal ratio, the extinction curve, and the fraction of directly escaped light from stars (\fesc). Our stellar mass function is found to be in broad agreement with Gonzalez et al., independent of these parameters. We find that our heavily dust-attenuated galaxies (A_V~1.8) can also reasonably match modest UV-optical colors, blue UV slopes, as well as UV luminosity functions, provided that a significant fraction (~10%) of light directly escapes from them. The observed UV slope and scatter are better explained with an Small Magellanic Cloud-type extinction curve, whereas Milky Way-type dust predicts too blue UV colors due to the 2175A bump. We expect that upcoming observations by ALMA will be able to test this heavily obscured model.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Infrared Properties of z=7 Galaxies from Cosmological Simulations

Three-dimensional panchromatic dust radiative transfer calculations are performed on a set of 198 galaxies of stellar masses in the range 5x10^8-3x10^10 Msun from a cosmological hydrodynamic simulation (resolved at 29pc/h) at z=7. In a companion paper (Kimm & Cen), the stellar mass and UV luminosity functions, and UV-optical and FUV-NUV colors are shown to be in good agreement with observations, if an SMC-type dust extinction curve is adopted. Here we make useful predictions, self-consistently, of the infrared properties of these z=7 simulated galaxies that can be confronted with upcoming ALMA data. Our findings are as follows. (1) The effective radius in the rest-frame MIPS 70 micron band is in the range of 80-400pc proper for z=7 galaxies with L_FIR=10^{11.3-12}Lsun. (2) The median of the peak wavelength of the far-infrared (FIR) spectral energy distribution is in the range of 45-60 micron, depending on the dust-to-metal ratio. (3) For star formation rate in the range 3-100 Msun/yr the median FIR to bolometric luminosity ratio is 60-90%. (4) The FIR luminosity function displays a power law in the high end with a slope of -3.1 +- 0.4, instead of the usual exponential decline.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Validity of Hydrostatic Equilibrium in Galaxy Clusters from Cosmological Hydrodynamical Simulations

We examine the validity of the hydrostatic equilibrium (HSE) assumption for galaxy clusters using one of the highest-resolution cosmological hydrodynamical simulations. We define and evaluate several effective mass terms corresponding to the Euler equations of the gas dynamics, and quantify the degree of the validity of HSE in terms of the mass estimate. We find that the mass estimated under the HSE assumption (the HSE mass) deviates from the true mass by up to ~ 30 %. This level of departure from HSE is consistent with the previous claims, but our physical interpretation is rather different. We demonstrate that the inertial term in the Euler equations makes a negligible contribution to the total mass, and the overall gravity of the cluster is balanced by the thermal gas pressure gradient and the gas acceleration term. Indeed the deviation from the HSE mass is well explained by the acceleration term at almost all radii. We also clarify the confusion of previous work due to the inappropriate application of the Jeans equations in considering the validity of HSE from the gas dynamics extracted from cosmological hydrodynamical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Calibrating an updated SPH scheme within GCD+

We adapt a modern scheme of smoothed particle hydrodynamics (SPH) to our tree N-body/SPH galactic chemodynamics code GCD+. The applied scheme includes imple- mentations of the artificial viscosity switch and artificial thermal conductivity pro- posed by Morris & Monaghan (1997), Rosswog & Price (2007) and Price (2008), to model discontinuities and Kelvin-Helmholtz instabilities more accurately. We first present hydrodynamics test simulations and contrast the results to runs undertaken without artificial viscosity switch or thermal conduction. In addition, we also explore the different levels of smoothing by adopting larger or smaller smoothing lengths, i.e. a larger or smaller number of neighbour particles, Nnb. We demonstrate that the new version of GCD+ is capable of modelling Kelvin-Helmholtz instabilities to a simi- lar level as the mesh code, Athena. From the Gresho vortex, point-like explosion and self-similar collapse tests, we conclude that setting the smoothing length to keep the number of neighbour particles as high as Nnb~58 is preferable to adopting smaller smoothing lengths. We present our optimised parameter sets from the hydrodynamics tests.

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

Coincidences between OVI and OVII Lines: Insights from High Resolution Simulations of the Warm-Hot Intergalactic Medium

With high resolution (0.46kpc/h), adaptive mesh-refinement Eulerian cosmological hydrodynamic simulations we compute properties of O VI and O VII absorbers from the warm-hot intergalactic medium (WHIM). Our new simulations are in broad agreement with previous simulations, with ~40% of the intergalactic medium being in the WHIM at z=0. It is found (1) The amount of gas in the WHIM at temperature below and above 10^6K is about equal within uncertainties. (1) Our simulations are in excellent agreement with observed properties of O VI absorbers, with respect to the line incidence rate and Doppler width-column density relation. (2) Velocity structures within absorbing regions are a significant, and for large Doppler width clouds, a dominant contributor to the Doppler widths of both O VI and O VII absorbers. A non-negligible fraction (in number and mass) of O VI and O VII clouds can arise from gas of temperature lower than 10^5, until the Doppler width is well in excess of 100km/s. (3) Strong O VI absorbers are predominantly collisionally ionized. About (61%, 57%, 39%) of O VI absorbers in the column density ranges of log N(OVI) cm^2=(12.5-13,13-14,>14) have temperature lower than 10^5K. (4) Quantitative prediction is made for the presence of broad and shallow O VI lines, which current observations may have largely missed. Upcoming observations by COS may be able to provide a test. (5) The reported 3 sigma upper limit on the mean column density of coincidental O VII lines at the location of detected O VI lines by Yao et al is above the predicted value by a factor of 2.5-4. (6) The claimed observational detection of O VII lines by Nicastro et al, if true, is 2 sigma above what our simulations predict.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Comparisons of Cosmological MHD Galaxy Cluster Simulations to Radio Observations

Radio observations of galaxy clusters show that there are $μ$G magnetic fields permeating the intra-cluster medium (ICM), but it is hard to accurately constrain the strength and structure of the magnetic fields without the help of advanced computer simulations. We present qualitative comparisons of synthetic VLA observations of simulated galaxy clusters to radio observations of Faraday Rotation Measure (RM) and radio halos. The cluster formation is modeled using adaptive mesh refinement (AMR) magneto-hydrodynamic (MHD) simulations with the assumption that the initial magnetic fields are injected into the ICM by active galactic nuclei (AGNs) at high redshift. In addition to simulated clusters in Xu et al. (2010, 2011), we present a new simulation with magnetic field injections from multiple AGNs. We find that the cluster with multiple injection sources is magnetized to a similar level as in previous simulations with a single AGN. The RM profiles from simulated clusters, both $|RM|$ and the dispersion of RM ($σ_{RM}$), are consistent at a first-order with the radial distribution from observations. The correlations between the $σ_{RM}$ and X-ray surface brightness from simulations are in a broad agreement with the observations, although there is an indication that the simulated clusters could be slightly over-dense and less magnetized with respect to those in the observed sample. In addition, the simulated radio halos agree with the observed correlations between the radio power versus the cluster X-ray luminosity and between the radio power versus the radio halo size. These studies show that the cluster wide magnetic fields that originate from AGNs and are then amplified by the ICM turbulence (Xu et al. 2010) match observations of magnetic fields in galaxy clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Effects on Galaxy Evolution: Pair Interactions versus Environment

In a hierarchical universe, mergers may be an important mechanism not only in increasing the mass of galaxies but also in driving the color and morphological evolution of galaxies. We use a large sample of ~1000 simulated galaxies of stellar mass greater than 10^9.6 solar masses (for ~4800 observations at multiple redshifts) from a high-res (0.46 h^{-1} kpc) cosmological simulation to determine under what circumstances being a member of a pair influences galaxy properties at z <= 0.2. We identify gravitationally bound pairs, and find a relative fraction of blue-blue, red-red, and blue-red pairs that agrees with observations (Lin et al. 2010). Pairs tend to avoid the extreme environments of clusters and void centres. While pairs in groups can include galaxies that are both blue, both red, or one of each color, in the field it is rare for pair galaxies to both be red. We find that physically bound pairs closer than 250 h^{-1} kpc tend to have higher sSFRs than the galaxy population as a whole. However, the sSFR of a bound galaxy relative to galaxies in a comparable local density environment (determined by the distance to the fifth nearest neighbor, rho_5), differs depending on the local density. In regions of high rho_5 the bound population has a higher fraction of star-forming (bluer) galaxies, whereas there is very little difference between bound and unbound galaxies in low rho_5 regions. This effect on the star-forming fraction may be driven by the higher fraction of bound HI-rich galaxies compared to unbound galaxies, particularly at high local densities. It appears that being in a pair has an incremental, but not overwhelming, effect on the star formation rate of the paired galaxies, compared to the more pronounced trend where galaxies overall have low sSFR (are red) in clusters and higher sSFR (blue) at the centre of voids. This trend depends most strongly on rho_5.(abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Nature of Lyman Alpha Blobs: Powered by Extreme Starbursts

We present a new model for the observed Lyman alpha blobs (LABs) within the context of the standard cold dark matter model. In this model, LABs are the most massive halos with the strongest clustering (proto-clusters) undergoing extreme starbursts in the high-z universe. Aided by calculations of detailed radiative transfer of Lya photons through ultra-high resolution (159pc) large-scale (>30Mpc) adaptive mesh-refinement cosmological hydrodynamic simulations with galaxy formation, this model is shown to be able to, for the first time, reproduce simultaneously the global Lya luminosity function and luminosity-size relation of the observed LABs. Physically, a combination of dust attenuation of Lya photons within galaxies, clustering of galaxies, and complex propagation of Lya photons through circumgalactic and intergalactic medium gives rise to the large sizes and frequently irregular isophotal shapes of LABs that are observed. A generic and unique prediction of this model is that there should be strong far-infrared (FIR) sources within each LAB, with the most luminous FIR source likely representing the gravitational center of the proto-cluster, not necessarily the apparent center of the Lya emission of the LAB or the most luminous optical source. Upcoming ALMA observations should unambiguously test this prediction. If verified, LABs will provide very valuable laboratories for studying formation of galaxies in the most overdense regions of the universe at a time when global star formation is most vigorous.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Physics of Coevolution of Galaxies and Supermassive Black Holes

A new model for coevolution of galaxies and supermassive black holes (SMBH) is presented that is physically based. The evolutionary track starts with an event that triggers a significant starburst in the central region of a galaxy. In this model, the main SMBH growth takes place in post-starburst phase fueled by recycled gas from inner bulge stars in a self-regulated fashion on a time scale that is substantially longer than 100Myrs and at a diminishing Eddington ratio with time. We argue that the SMBH cannot gorge itself during the starburst phase, despite the abundant supply of cold gas, because star formation is a preferred mode of gas consumption in such an environment than accretion to the central SMBH. We also show that feedback from star formation is at least as strong as that from AGN and thus, if star formation is in need of being quenched, AGN feedback generally does not play the primary role. The predicted relation between SMBH mass and bulge mass/velocity dispersion is consistent with observations. A clear prediction is that early-type galaxy hosts of high Eddingtion rate AGNs are expected to be light-blue to green in optical color, gradually evolving to the red sequences with decreasing AGN luminosity. A suite of falsifiable predictions and implications with respect to relationships between various types of galaxies and AGN, and others, are made. For those where comparisons to extant observations are possible, the model appears to be in good standing.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Reionization on Large Scales I: A Parametric Model Constructed from Radiation-Hydrodynamic Simulations

We present a new method for modeling inhomogeneous cosmic reionization on large scales. Utilizing high-resolution radiation-hydrodynamic simulations with 2048^3 dark matter particles, 2048^3 gas cells, and 17 billion adaptive rays in a L = 100 Mpc/h box, we show that the density and reionization-redshift fields are highly correlated on large scales (>~ 1 Mpc/h). This correlation can be statistically represented by a scale-dependent linear bias. We construct a parametric function for the bias, which is then used to filter any large-scale density field to derive the corresponding spatially varying reionization-redshift field. The parametric model has three free parameters which can be reduced to one free parameter when we fit the two bias parameters to simulations results. We can differentiate degenerate combinations of the bias parameters by combining results for the global ionization histories and correlation length between ionized regions. Unlike previous semi-analytic models, the evolution of the reionization-redshift field in our model is directly compared cell by cell against simulations and preforms well in all tests. Our model maps the high-resolution, intermediate-volume radiation-hydrodynamic simulations onto lower-resolution, larger-volume N-body simulations (>~ 2 Gpc/h) in order to make mock observations and theoretical predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Reionization on Large Scales III: Predictions for Low-ell Cosmic Microwave Background Polarization and High-ell Kinetic Sunyaev-Zel'dovich Observables

We present new predictions for temperature (on small angular scales) and polarization (on large angular scales) CMB anisotropies induced during the epoch of reionization (EoR). Using a novel method calibrated from Radiation-Hydrodynamic simulations we model the EoR in large volumes (L >~ 2 Gpc/h) in the context of galactic reionization. We find that the EoR contribution to the kinetic Sunyaev- Zel'dovich power spectrum (patchy kSZ) ranges between ~0.6 - 2.8 muK^2 at ell = 3000, for the parameter space we explored. These patchy kSZ power spectra are calculated from large 15 Deg x 15 Deg maps that are found to be necessary. Decreasing the size of these maps biases the overall patchy kSZ power to higher values. We find that the amplitude of the patchy kSZ power spectrum at ell = 3000 follows simple scalings of D_ell=3000^kSZ propto <z> and D_ell=3000^kSZ propto Delz^0.47 for the mean redshift (<z>) of reionization and the duration (dz). Using the constraints on <z> from WMAP 7-year results and the lower limit on dz from EDGES we find a lower limit of ~ 0.4 muK^2 on the kSZ at ell = 3000. Planck will constrain the mean redshift and the Thomson optical depth from the low-ell polarization power spectrum. Future measurements of the high-ell CMB power spectrum from the South Pole Telescope (SPT) and the Atacama Cosmology Telescope (ACT) should detect the patchy kSZ signal if the cross correlation between the cosmic infrared background and the thermal Sunyaev Zel'dovich effect is constrained. We show that the combination of temperature and polarization measurements constrains both <z> and dz. The patchy kSZ maps, power spectra templates and the polarization power spectra will be publicly available.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Environmentally Driven Global Evolution of Galaxies

Utilizing high-resolution large-scale galaxy formation simulations of the standard cold dark matter model, we examine global trends in the evolution of galaxies due to gravitational shock heating by collapse of large halos and large-scale structure. We find two major global trends. (1) The mean specific star formation rate (sSFR) at a given galaxy mass is a monotonically increasing function with increasing redshift. (2) The mean sSFR at a given redshift is a monotonically increasing function of decreasing galaxy mass that steepens with decreasing redshift. The general dimming trend with time merely reflects the general decline of gas inflow rate with increasing time. The differential evolution of galaxies of different masses with redshift is a result of gravitational shock heating of gas due to formation of large halos (groups and clusters) and large-scale structure that move a progressively larger fraction of galaxies and their satellites into environments where gas has too high an entropy to cool to continue feeding resident galaxies. Overdense regions where larger halos are preferentially located begin to be heated earlier and have higher temperatures than lower density regions at any given time, causing sSFR of larger galaxies to fall below the general dimming trend at higher redshift than less massive galaxies and galaxies with high sSFR to gradually shift to lower density environments at lower redshift. We find that several noted cosmic downsizing phenomena are different manifestations of these general trends. We also find that the great migration of galaxies from blue cloud to red sequence as well as color-density relation, among others, may arise naturally in this picture.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Extended Lyman-Alpha Emission around Star-forming Galaxies

Lyman-alpha (Lya) photons that escape the interstellar medium of star-forming galaxies may be resonantly scattered by neutral hydrogen atoms in the circumgalactic and intergalactic media, thereby increasing the angular extent of the galaxy's Lya emission. We present predictions of this extended, low surface brightness Lya emission based on radiative transfer modeling in a cosmological reionization simulation. The extended emission can be detected from stacked narrowband images of Lya emitters (LAEs) or of Lyman break galaxies (LBGs). Its average surface brightness profile has a central cusp, then flattens to an approximate plateau beginning at an inner characteristic scale below ~0.2 Mpc (comoving), then steepens again beyond an outer characteristic scale of ~1 Mpc. The inner scale marks the transition from scattered light of the central source to emission from clustered sources, while the outer scale marks the spatial extent of scattered emission from these clustered sources. Both scales tend to increase with halo mass, UV luminosity, and observed Lya luminosity. The extended emission predicted by our simulation is already within reach of deep narrowband photometry using large ground-based telescopes. Such observations would test radiative transfer models of emission from LAEs and LBGs, and they would open a new window on the circumgalactic environment of high-redshift star-forming galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Far-Infrared Properties of Lyman Break Galaxies from Cosmological Simulations

Utilizing state-of-the-art, adaptive mesh-refinement cosmological hydrodynamic simulations with ultra-high resolution (114h-1pc) and large sample size (>3300 galaxies of stellar mass >10^9Msun), we show how the stellar light of Lyman Break Galaxies at z=2 is distributed between optical/ultra-violet (UV) and far-infrared (FIR) bands. With a single scalar parameter for dust obscuration we can simultaneously reproduce the observed UV luminosity function for the entire range (3-100 Msun/yr) and extant FIR luminosity function at the bright end (>20Msun/yr). We quantify that galaxies more massive or having higher SFR tend to have larger amounts of dust obscuration mostly due to a trend in column density and in a minor part due to a mass (or SFR)-metallicity relation. It is predicted that the FIR luminosity function in the range SFR=1-100Msun/yr is a powerlaw with a slope about -1.7. We further predict that there is a "galaxy desert" at SFR(FIR) < 0.02 (SFR(UV)/10Msun/yr)^2.1 Msun/yr in the SFR(UV)-SFR(FIR) plane. Detailed distributions of SFR(FIR) at a fixed SFR(UV) are presented. Upcoming observations by ALMA should test this model. If confirmed, it validates the predictions of the standard cold dark matter model and has important implications on the intrinsic SFR function of galaxies at high redshift.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Grand Unification II: Hot Accretion and AGN Jets

We show that a quasi-spherical (QS) hot accretion flow is expected to operate in all SMBHs, with its rate being capped at dot m_QS,max~0.001 (M/1E8Msun) in units of the Eddington rate. It is then proposed that AGN jet power is proportional to the product of the hot accretion rate and a SMBH spin-dependent efficiency for energy extraction in the form of jets. Predictions from this model include (1) while radio jets should emerge from all SMBHs, the maximum jet power goes with SMBH mass approximately as 1E43.6(M/1E8Msun)^{2}erg/s, (2) even although there are two separate underlying populations of radio-quiet (RQ) and radio-loud (RL) AGNs, any bimodality in the observed AGN radio loudness distribution is likely due to a selection effect, such as some imposed optical magnitude limits, (3) the RL fraction of quasars is expected to decrease with redshift in the cold dark matter model, (4) host galaxies of RQ and RL quasars should be drawn from the same underlying elliptical galaxy population, although RL quasars may have SMBHs that are somewhat more massive than their RQ counterparts and RL quasars may reside predominantly in core elliptical galaxies, (5) RL low-luminosity AGNs and LINERs may represent the long, declining "trailing" phase following the initial, more luminous AGN phase, (6) a broad anti-correlation between radio- loudness and disk accretion rate is expected, (7) RL AGNs may be expected to be more abundant in Type Ia supernovae than RQ AGNs, (8) among RL AGNs a correlation between radio power and clustering strength is predicted, and (9) RQ and RL AGNs should have, on average, similar IR-optical-UV properties. ~

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Inconsequence of Galaxy Major Mergers in Driving Star Formation at z>1: Insights from Cosmological Simulations

Utilizing a high-resolution (114 pc/h) adaptive mesh-refinement cosmological galaxy formation simulation of the standard cold dark matter model with a large (2000-3000 galaxies with stellar mass greater than 1e9 Msun) statistical sample, we examine the role of major mergers in driving star formation at z>1 in a cosmological setting, after validating that some of the key properties of simulated galaxies are in reasonable agreement with observations, including luminosity functions, SF history, effective sizes and damped Lyman alpha systems. We find that major mergers have a relatively modest effect on star formation, in marked contrast to previous idealized merger simulations of disk galaxies that show up to two orders of magnitude increase in star formation rate. At z=2.4-3.7, major mergers tend to increase the specific star formation rate by 10-25% for galaxies in the entire stellar mass range 10^9-10^12 Msun probed. Their effect appears to increase with decreasing redshift, but is capped at 60% at z=1.4-2.4. Two factors may account for this modest effect. First, SFR of galaxies not in major mergers are much higher at z>1 than local disk galaxy counterparts. Second, most galaxies at z>1 have small sizes and contain massive dense bulges, which suppress the merger induced structural effects and gas inflow enhancement. Various other predictions are also made that will provide verifiable tests of the model.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Simulated void galaxies in the standard cold dark matter model

We analyze a (120 h^{-1} Mpc)^3 adaptive mesh refinement hydrodynamic simulation that contains a higher-resolution 31 x 31 x 35 h^{-3} Mpc subvolume centered on a ~30 Mpc diameter void. Our detailed ~1 kpc resolution allows us to identify 1300 galaxies within this void to a limiting halo mass of ~10^{10} M_sun. Nearly 1000 galaxies are found to be in underdense regions, with 300 galaxies residing in regions less than half the mean density of the simulation volume. We construct mock observations of the stellar and gas properties of these systems, and reproduce the range of colors and luminosities observed in the SDSS for nearby (z < 0.03) galaxies. We find no trends with density for the most luminous (M_r < -18) galaxies, however our dwarf void galaxies (M_r > -16), though they are less reliably resolved, typically appear bluer, with higher rates of star formation and specific star formation and lower mean stellar ages than galaxies in average density environments. We find a significant population of low luminosity (M_r ~ -14) dwarf galaxies that is preferentially located in low density regions and specifically in the void center. This population may help to reduce, but not remove, the discrepancy between the predicted and observed number of void galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Star Formation Feedback and Metal Enrichment History Of The Intergalactic Medium

Using hydrodynamic simulations we compute the metal enrichment history of the intergalactic medium (IGM). We show that galactic superwind (GSW) feedback can transport metals to the IGM and that the properties of simulated metal absorbers match observations. The distance of influence of GSW is typically limited to >0.5Mpc and within regions of overdensity >10. Most CIV and OVI absorbers are located within shocked regions of elevated temperature (T>2x10^4K), overdensity (>10), and metallicity ([-2.5,-0.5]). OVI absorbers have typically higher metallicity, lower density and higher temperature than CIV absorbers. For OVI absorbers collisional ionization dominates over the entire redshift range z=0-6, whereas for CIV absorbers the transition occurs at moderate redshift z~3 from collisionally dominated to photoionization dominated. We find that the observed column density distributions for CIV and OVI in the range log N cm^2=12-15 are reasonably reproduced by the simulations. The evolution of mass densities contained in CIV and OVI lines, Omega_CIV and Omega_OVI, is also in good agreement with observations, which shows a near constancy at low redshifts and an exponential drop beyond redshift z=3-4. For both CIV and OVI, most absorbers are transient and the amount of metals probed by CIV and OVI lines of column log N cm^2=12-15 is only ~2% of total metal density at any epoch. While gravitational shocks from large-scale structure formation dominate the energy budget (80-90%) for turning about 50% of IGM to the warm-hot intergalactic medium (WHIM) by z=0, GSW feedback shocks are energetically dominant over gravitational shocks at z > 1-2. Most of the so-called "missing metals" at z=2-3 are hidden in a warm-hot (T=10^{4.5-7}K) gaseous phase, heated up by GSW feedback shocks. Their mass distribution is broadly peaked at $δ=1-10$ in the IGM, outside virialized halos.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

21cmFAST: A Fast, Semi-Numerical Simulation of the High-Redshift 21-cm Signal

We introduce a powerful semi-numeric modeling tool, 21cmFAST, designed to efficiently simulate the cosmological 21-cm signal. Our code generates 3D realizations of evolved density, ionization, peculiar velocity, and spin temperature fields, which it then combines to compute the 21-cm brightness temperature. Although the physical processes are treated with approximate methods, we compare our results to a state-of-the-art large-scale hydrodynamic simulation, and find good agreement on scales pertinent to the upcoming observations (>~ 1 Mpc). The power spectra from 21cmFAST agree with those generated from the numerical simulation to within 10s of percent, down to the Nyquist frequency. We show results from a 1 Gpc simulation which tracks the cosmic 21-cm signal down from z=250, highlighting the various interesting epochs. Depending on the desired resolution, 21cmFAST can compute a redshift realization on a single processor in just a few minutes. Our code is fast, efficient, customizable and publicly available, making it a useful tool for 21-cm parameter studies.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Comparison Of Reionization Models: Radiative Transfer Simulations And Approximate, Semi-Numeric Models

We compare the predictions of four different algorithms for the distribution of ionized gas during the Epoch of Reionization. These algorithms are all used to run a 100 Mpc/h simulation of reionization with the same initial conditions. Two of the algorithms are state-of-the-art ray-tracing radiative transfer codes that use disparate methods to calculate the ionization history. The other two algorithms are fast but more approximate schemes based on iterative application of a smoothing filter to the underlying source and density fields. We compare these algorithms' resulting ionization and 21 cm fields using several different statistical measures. The two radiative transfer schemes are in excellent agreement with each other (with the cross-correlation coefficient of the ionization fields >0.8 for k < 10 h/Mpc and in good agreement with the analytic schemes (>0.6 for k < 1 h/Mpc). When used to predict the 21cm power spectrum at different times during reionization, all ionization algorithms agree with one another at the 10s of percent level. This agreement suggests that the different approximations involved in the ray tracing algorithms are sensible and that semi-numerical schemes provide a numerically-inexpensive, yet fairly accurate, description of the reionization process.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Crawling the Cosmic Network: Exploring the Morphology of Structure in the Galaxy Distribution

Although coherent large-scale structures such as filaments and walls are apparent to the eye in galaxy redshift surveys, they have so far proven difficult to characterize with computer algorithms. This paper presents a procedure that uses the eigenvalues and eigenvectors of the Hessian matrix of the galaxy density field to characterize the morphology of large-scale structure. By analysing the smoothed density field and its Hessian matrix, we can determine the types of structure - walls, filaments, or clumps - that dominate the large-scale distribution of galaxies as a function of scale. We have run the algorithm on mock galaxy distributions in a LCDM cosmological N-body simulation and the observed galaxy distributions in the Sloan Digital Sky Survey. The morphology of structure is similar between the two catalogues, both being filament-dominated on 10-20 h^{-1} Mpc smoothing scales and clump-dominated on 5 h^{-1} Mpc scales. There is evidence for walls in both distributions, but walls are not the dominant structures on scales smaller than ~25 h^{-1} Mpc. Analysis of the simulation suggests that, on a given comoving smoothing scale, structures evolve with time from walls to filaments to clumps, where those found on smaller smoothing scales are further in this progression at a given time.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Crawling the Cosmic Network: Identifying and Quantifying Filamentary Structure

We present the Smoothed Hessian Major Axis Filament Finder (SHMAFF), an algorithm that uses the eigenvectors of the Hessian matrix of the smoothed galaxy distribution to identify individual filamentary structures. Filaments are traced along the Hessian eigenvector corresponding to the largest eigenvalue, and are stopped when the axis orientation changes more rapidly than a preset threshold. In both N-body simulations and the Sloan Digital Sky Survey (SDSS) main galaxy redshift survey data, the resulting filament length distributions are approximately exponential. In the SDSS galaxy distribution, using smoothing lengths of 10 h^{-1} Mpc and 15 h^{-1} Mpc, we find filament lengths per unit volume of 1.9x10^{-3} h^2 Mpc^{-2} and 7.6x10^{-4} h^2 Mpc^{-2}, respectively. The filament width distributions, which are much more sensitive to non-linear growth, are also consistent between the real and mock galaxy distributions using a standard cosmology. In SDSS, we find mean filament widths of 5.5 h^{-1} Mpc and 8.4 h^{-1} Mpc on 10 h^{-1} Mpc and 15 h^{-1} Mpc smoothing scales, with standard deviations of 1.1 h^{-1} Mpc and 1.4 h^{-1} Mpc, respectively. Finally, the spatial distribution of filamentary structure in simulations is very similar between z=3 and z=0 on smoothing scales as large as 15 h^{-1} Mpc, suggesting that the outline of filamentary structure is already in place at high redshift.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Radiative Transfer Modeling of Lyman Alpha Emitters. I. Statistics of Spectra and Luminosity

We combine a cosmological reionization simulation with box size of 100Mpc/h on a side and a Monte Carlo Lyman-alpha (Lya) radiative transfer code to model Lyman Alpha Emitters (LAEs) at z~5.7. The model introduces Lya radiative transfer as the single factor for transforming the intrinsic Lya emission properties into the observed ones. Spatial diffusion of Lya photons from radiative transfer results in extended Lya emission and only the central part with high surface brightness can be observed. Because of radiative transfer, the appearance of LAEs depends on density and velocity structures in circumgalactic and intergalactic media as well as the viewing angle, which leads to a broad distribution of apparent (observed) Lya luminosity for a given intrinsic Lya luminosity. Radiative transfer also causes frequency diffusion of Lya photons. The resultant Lya line is asymmetric with a red tail. The peak of the Lya line shifts towards longer wavelength and the shift is anti-correlated with the apparent to intrinsic Lya luminosity ratio. The simple radiative transfer model provides a new framework for studying LAEs. It is able to explain an array of observed properties of z~5.7 LAEs in Ouchi et al. (2008), producing Lya spectra, morphology, and apparent Lya luminosity function (LF) similar to those seen in observation. The broad distribution of apparent Lya luminosity at fixed UV luminosity provides a natural explanation for the observed UV LF, especially the turnover towards the low luminosity end. The model also reproduces the observed distribution of Lya equivalent width (EW) and explains the deficit of UV bright, high EW sources. Because of the broad distribution of the apparent to intrinsic Lya luminosity ratio, the model predicts effective duty cycles and Lya escape fractions for LAEs.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Radiative Transfer Modeling of Lyman Alpha Emitters. II. New Effects in Galaxy Clustering

We study the clustering properties of z~5.7 Lyman-alpha emitters (LAEs) in a cosmological reionization simulation with a full Lya radiative transfer calculation. Lya radiative transfer substantially modifies the intrinsic Lya emission properties, compared to observed ones, depending on the density and velocity structure environment around the Lya emitting galaxy. This environment-dependent Lya selection introduces new features in LAE clustering, suppressing (enhancing) the line-of-sight (transverse) density fluctuations and giving rise to scale-dependent galaxy bias. In real space, the contours of the three-dimensional two-point correlation function of LAEs appear to be prominently elongated along the line of sight on large scales, an effect that is opposite to and much stronger than the linear redshift-space distortion effect. The projected two-point correlation function is greatly enhanced in amplitude by a factor of up to a few, compared to the case without the environment dependent selection effect. The new features in LAE clustering can be understood with a simple, physically motivated model, where Lya selection depends on matter density, velocity, and their gradients. We discuss the implications and consequences of the effects on galaxy clustering from Lya selection in interpreting clustering measurements and in constraining cosmology and reionization from LAEs.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Nature of Damped Lyman Alpha Systems and Their Hosts in the Standard Cold Dark Matter Universe

Using adaptive mesh-refinement cosmological hydrodynamic simulations with a physically motivated supernova feedback prescription we show that the standard cold dark matter model can account for extant observed properties of damped Lyman alpha systems (DLAs). We then examine the properties of DLA host galaxies. We find: (1) While DLA hosts roughly trace the overall population of galaxies at all redshifts, they are always gas rich. (2) The history of DLA evolution reflects primarily the evolution of the underlying cosmic density, galaxy size and galaxy interactions. With higher density and more interactions at high redshift DLAs are larger in both absolute terms and in relative terms with respect to virial radii of halos. (3) The variety of DLAs at high redshift is richer with a large contribution coming from galactic filaments, created through close galaxy interactions. The portion of gaseous disks of galaxies where most stars reside makes relatively small contribution to DLA incidence at z=3-4. (4) The vast majority of DLAs arise in halos of mass M_h=10^10-10^12 Msun at z=1.6-4. At z=3-4, 20-30% of DLA hosts are Lyman Break Galaxies (LBGs). (5) Galactic winds play an indispensable role in shaping the kinematic properties of DLAs. Specifically, the high velocity width DLAs are a mixture of those arising in high mass, high velocity dispersion halos and those arising in smaller mass systems where cold gas clouds are entrained to high velocities by galactic winds. (6) In agreement with observations, we see a weak but noticeable evolution in DLA metallicity. The metallicity distribution centers at [Z/H]=-1.5 to -1 at z=3-4, with the peak moving to [Z/H]=-0.75 at z=1.6 and [Z/H]=-0.5 by z=0. (7) The star formation rate of DLA hosts is concentrated in the range 0.3-30Msun/yr at z=3-4, gradually shifting lower to peak at ~0.5-1 Msun/yr by z=0.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The State of the Universe at z~6

In the context stellar reionization in the standard cold dark matter model, we analyze observations at z~6 and are able to draw three significant conclusions with respect to star formation and the state of the intergalactic medium (IGM) at z~6. (1) An initial stellar mass function (IMF) more efficient, by a factor of 10-20, in producing ionizing photons than the standard Salpeter IMF is required at z~6. This may be achieved by having either (A) a metal-enriched IMF with and a lower mass cutoff of >= 30Msun or (B) 2-4% of stellar mass being Population III massive metal-free stars at z~6. While there is no compelling physical reason or observational evidence to support (A), (B) could be fulfilled plausibly by continued existence of some pockets of uncontaminated, metal-free gas for star formation. (2) The volume-weighted neutral fraction of the IGM of <f_HI>_V~ 10^-4 at z=5.8 inferred from the SDSS observations of QSO absorption spectra provides enough information to ascertain that reionization is basically complete with at most ~0.1-1% of IGM that is un-ionized at z=5.8. (3) Barring some extreme evolution of the IMF, the neutral fraction of the IGM is expected to rise quickly toward high redshift from the point of HII bubble percolation, with the mean neutral fraction of the IGM expected to reach 6-12% at z=6.5, 13-27% at z=7.7 and 22-38% at z=8.8.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The 21 cm Forest as a Probe of the Reionization and the Temperature of the Intergalactic Medium

Using high redshift radio sources as background, the 21cm forest observations probe the neutral hydrogen absorption signatures of early structures along the lines of sight. Directly sensitive to the spin temperature, it complements the 21cm tomography observations, and provides information on the temperature as well as the ionization state of the intergalactic medium (IGM). We use a radiative transfer simulation to investigate the 21cm forest signals during the epoch of reionization. We first confirmed that the optical depth and equivalent width (EW) are good representations of the ionization and thermal state of the IGM. The features selected by their relative optical depth are excellent tracers of ionization fields, while those selected by absolute optical depth are very sensitive to the IGM temperature, so the IGM temperature information could potentially be extracted, thus breaking a degeneracy in 21cm tomographic observation. From the distributions of EWs and the number evolution of absorbers and leakers with different EWs, we see clearly the cosmological evolution of ionization state of the IGM. The number density of potentially observable features decreases rapidly with increasing gas temperature. The sensitivity of the proposed EW statistic to the IGM temperature makes it a unique and potentially powerful probe of reionization. Missing small-scale structures, such as small filaments and minihalos that are unresolved in our current simulation, and lack of an accurate calculation of the IGM temperature, however, likely have rendered the presented signals quantitatively inaccurate. Finally we discuss the requirements of the background radio sources for such observations, and find that signals with equivalent widths larger than 1kHz are hopeful to be detected.(Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Ionizing Photon Escape Fractions from High Redshift Dwarf Galaxies

It has been argued that low-luminosity dwarf galaxies are the dominant source of ionizing radiation during cosmological reionization. The fraction of ionizing radiation that escapes into the intergalactic medium from dwarf galaxies with masses less than ~10^9.5 solar masses plays a critical role during this epoch. Using an extensive suite of very high resolution (0.1 pc), adaptive mesh refinement, radiation hydrodynamical simulations of idealized and cosmological dwarf galaxies, we characterize the behavior of the escape fraction in galaxies between 3 x 10^6 and 3 x 10^9 solar masses with different spin parameters, amounts of turbulence, and baryon mass fractions. For a given halo mass, escape fractions can vary up to a factor of two, depending on the initial setup of the idealized halo. In a cosmological setting, we find that the time-averaged photon escape fraction always exceeds 25% and reaches up to 80% in halos with masses above 10^8 solar masses with a top-heavy IMF. The instantaneous escape fraction can vary up to an order of magnitude in a few million years and tend to be positively correlated with star formation rate. We find that the mean of the star formation efficiency times ionizing photon escape fraction, averaged over all atomic cooling (T_vir > 8000 K) galaxies, ranges from 0.02 for a normal IMF to 0.03 for a top-heavy IMF, whereas smaller, molecular cooling galaxies in minihalos do not make a significant contribution to reionizing the universe due to a much lower star formation efficiency. These results provide the physical basis for cosmological reionization by stellar sources, predominately atomic cooling dwarf galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

CosmoMHD: A Cosmological Magnetohydrodynamics Code

In this era of precision cosmology, a detailed physical understanding on the evolution of cosmic baryons is required. Cosmic magnetic fields, though still poorly understood, may represent an important component in the global cosmic energy flow that affects the baryon dynamics. We have developed an Eulerian-based cosmological magnetohydrodynamics code (CosmoMHD) with modern shock capturing schemes to study the formation and evolution of cosmic structures in the presence of magnetic fields. The code solves the ideal MHD equations as well as the non-equilibrium rate equations for multiple species, the Vlasov equation for dynamics of collisionless particles, the Poisson's equation for the gravitational potential field and the equation for the evolution of the intergalactic ionizing radiation field. In addition, a detailed star formation prescription and feedback processes are implemented. Several methods for solving the MHD by high-resolution schemes with finite-volume and finite-difference methods are implemented. The divergence-free condition of the magnetic fields is preserved at a level of computer roundoff error via the constraint transport method. We have also implemented a high-resolution method via dual-equation formulations to track the thermal energy accurately in very high Mach number or high Alfven-Mach number regions. Several numerical tests have demonstrated the efficacy of the proposed schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

The Extended Starformation History of the First Generation of Stars, and the Reionization of Cosmic Hydrogen

Population-III (Pop-III) starformation (SF) is thought to be quenched when the metallicity of the star-forming gas reaches a critical level. At high z, when the general intergalactic medium (IGM) was enriched with metals, the fraction of primordial gas already collapsed in minihalos was significantly larger than the fraction of primordial gas that had already been involved in Pop-III SF. We argue that this minihalo gas remained largely in a metal-free state, until these minihalos merged into large systems and formed stars. As a result, the era of Pop-III SF was significantly prolonged, leading to an integrated Pop-III SF an order of magnitude larger than expected for an abrupt transition redshift. The contribution of Pop-III SF to the reionization of hydrogen could have been significant until z~10 and may have extended to z~6. Our modeling allows for gradual enrichment of the IGM, feedback from photo-ionization and screening of reionization by minihalos. Nevertheless, extended Pop-III SF can result in complex, multi-peaked reionization histories. The contribution of Pop-III stars to reionization will be tested by the three-year WMAP results: (1) if Pop-III stars do not contribute to reionization, tau_es<0.05-0.06 and a rapid reionization at z~6 is expected; (2) if the product of star formation efficiency and escape fraction for Pop-III stars is significantly larger than for Pop-II stars, then a maximum tau_es=0.21 is achievable; (3) in a scenario where the product of star formation efficiency and escape fraction for Pop-III stars is comparable to that for Pop-II stars, tau_es=0.09-0.12 would be observed, with reionization histories characterized by an extended ionization plateau from z=7-12. This result holds regardless of the redshift where the IGM becomes enriched with metals.

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
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
preprint1997arXiv

Sizes, Shapes, and Correlations of Lyman Alpha Clouds and Their Evolution in the CDM$+Λ$ Universe

This study analyzes the sizes, shapes and correlations of $\lya$ clouds produced by a hydrodynamic simulation of a spatially flat CDM universe with a non-zero cosmological constant ($Ω_0=0.4$, $Λ_0=0.6$, $σ_8 =0.79$), over the redshift range $2\le z \le 4$. The $\lya$ clouds range in size from several kiloparsecs to about a hundred kiloparsecs in proper units, and they range in shape from roundish, high column density regions with $\nhi\ge 10^{15} cm^{-2}$ to low column density sheet-like structures with $\nhi \le 10^{13} cm^{-2}$ at z=3. The most common shape found in the simulation resembles that of a flattened cigar. The physical size of a typical cloud grows with time roughly as $(1+z)^{-3/2}$ while its shape hardly evolves (except for the most dense regions $ρ_{cut}>30$). Our result indicates that any simple model with a population of spheres (or other shapes) of a uniform size is oversimplified; if such a model agrees with observational evidence, it is probably only by coincidence. We also illustrate why the use of double quasar sightlines to set lower limits on cloud sizes is useful only when the perpendicular sightline separation is small ($Δr \le 50h^{-1}$ kpc). Finally, we conjecture that high column density $\lya$ clouds ($\nhi\ge 10^{15} cm^{-2}$) may be the progenitors of the lower redshift faint blue galaxies. This seems plausible because their correlation length, number density (extrapolated to lower redshift) and their masses are in fair agreement with those observed.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Testing Cosmological Models With A $\lya$ Forest Statistic: The High End Of The Optical Depth Distribution

We pay particular attention to the high end of the $\lya$ optical depth distribution of a quasar spectrum. Based on the flux distribution (Miralda-Escudé et al 1996), a simple yet seemingly cosmological model -differentiating statistic, $Δ_{τ_0}$ -- the cumulative probability of a quasar spectrum with $\lya$ optical depth greater than a high value $τ_0$ -- is emphasized. It is shown that two different models -- the cold dark matter model with a cosmological constant and the mixed hot and cold dark matter model, both normalized to COBE and local galaxy cluster abundance -- yield quite different values of $Δ_{τ_0}$: 0.13 of the former versus 0.058 of the latter for $τ_0=3.0$ at $z=3$. Moreover, it is argued that $Δ_{τ_0}$ may be fairly robust to compute theoretically because it does not seem to depend sensitively on small variations of simulations parameters such as radiation field, cooling, feedback process, radiative transfer, resolution and simulation volume within the plausible ranges of the concerned quantities. Furthermore, it is illustrated that $Δ_{τ_0}$ can be obtained sufficiently accurately from currently available observed quasar spectra for $τ_0\sim 3.0-4.0$, when observational noise is properly taken into account. We anticipate that analyses of observations of quasar $\lya$ absorption spectra over a range of redshift may be able to constrain the redshift evolution of the amplitude of the density fluctuations on small-to-intermediate scales, therefore providing an independent constraint on $Ω_0$, $Ω_{0,HDM}$ and $Λ_0$.

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