Source author record

Yuexing Li

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

22works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

22 published item(s)

preprint2022arXiv

FOREVER22: galaxy formation in protocluster regions

We present results from a new cosmological hydrodynamics simulation campaign of protocluster (PC) regions, FOREVER22: FORmation and EVolution of galaxies in Extremely-overdense Regions motivated by SSA22. The simulations cover a wide range of cosmological scales using three different zoom set-ups in a parent volume of $(714.2~\rm cMpc)^{3}$: PCR (Proto-Cluster Region; $V= (28.6~{\rm cMpc})^{3} $, SPH particle mass, $m_{\rm{SPH}} = 4.1 \times 10^{6}~\rm M_{\odot}$ and final redshift, $z_{\rm end}=2.0$), BCG (Brightest proto-Cluster Galaxy; $V \sim (10~{\rm cMpc})^{3} $, $m_{\rm SPH} = 5.0\times10^{5}~\rm M_{\odot}$ and $z_{\rm end}=4.0$ ), and First ( $V \sim (3~{\rm cMpc})^{3} $, $m_{\rm SPH} = 7.9 \times 10^{3}~\rm M_{\odot}$ and $z_{\rm end}=9.5$) runs, that allow to focus on different aspects of galaxy formation. In the PCR runs, we follow 10 PCs, each harbouring 1 - 4 SMBHs with $M_{\rm BH} \ge 10^{9}~\rm M_{\odot}$. One of the PC cores shows a spatially close arrangement of seven starburst galaxies with ${\rm SFR} \gtrsim 100~\rm M_{\odot}~yr^{-1}$ each, that are dust-obscured and would appear as submillimeter galaxies with flux $\gtrsim 1~$ mJy at $1.1~ \rm mm$ in observations. The BCG runs show that the total SFRs of haloes hosting BCGs are affected by AGN feedback, but exceed $1000~\rm M_{\odot}~yr^{-1}$ at $z \lesssim 6$. The First runs resolve mini-haloes hosting population (Pop) III stars and we show that, in PC regions, the dominant stellar population changes from Pop III to Pop II at $z \gtrsim 20$, and the first galaxies with ${\rm SFR} \gtrsim 18~\rm M_{\odot}~yr^{-1}$ form at $z \sim 10$. These can be prime targets for future observations with the James Webb Space Telescope. Our simulations successfully reproduce the global star formation activities in observed PCs and suggest that PCs can kickstart cosmic reionization.

preprint2020arXiv

ART$^2$: A 3D Parallel Multi-wavelength Radiative Transfer Code for Continuum and Atomic and Molecular Lines

ART$^2$ is a 3D multi-wavelength Monte Carlo radiative transfer (RT) code that couples continuum and emission lines to track the propagation of photons and their interactions with the interstellar medium (ISM). The previous version of ART$^2$, which included continuum and Ly$α$ line, has been extensively applied to hydrodynamics simulations to study multi-band properties of galaxies and ISM. Here, we describe new implementations of non-local thermodynamic equilibrium RT of molecular and atomic fine structure emission lines, and the parallelization of the code using a number of novel methods. The new ART$^2$ can efficiently and self-consistently produce a full spectrum that includes both continuum and lines such as [CII], [NII], [OIII], Ly$α$, and CO. These essential features, together with the multi-phase ISM model and the adaptive grid, make ART$^2$ a multi-purpose code to study multi-wavelength properties of a wide range of astrophysical systems from planetary disks to large-scale structures. To demonstrate the capability of the new ART$^2$, we applied it to two hydrodynamics simulations: the zoom-in Milky Way Simulation to obtain panchromatic properties of individual galaxies, and the large-scale IllustrisTNG100 Simulation to obtain global properties such as the line intensity mappings. These products are essential for a broad array of studies such as the correlations between physical and panchromatic properties and their evolution. By enabling direct comparison between numerical simulations and multi-band observations, ART$^2$ provides a crucial theoretical framework for the interpretations of existing observations, the plan for future surveys, and the synergy between multi-band galaxy surveys and line intensity mappings. Therefore, ART$^2$ is a powerful and versatile tool to bridge the gap between theories and observations of cosmic structures.

preprint2019arXiv

Galaxy evolution and radiative properties in the early Universe: multi-wavelength analysis in cosmological simulations

Recent observations have successfully detected UV or infrared flux from galaxies at the epoch of reionization. However, the origin of their radiative properties has not been fully understood yet. Combining cosmological hydrodynamic simulations and radiative transfer calculations, we present theoretical predictions of multi-wavelength radiative properties of the first galaxies at z=6-15. We find that most of the gas and dust are ejected from star-forming regions due to supernova (SN) feedback, which allows UV photons to escape. We show that the peak of SED rapidly shifts between UV and infrared wavelengths on a timescale of 100 Myr due to intermittent star formation and feedback. When dusty gas covers the star-forming regions, the galaxies become bright in the observed-frame sub-millimeter wavelengths. In addition, we find that the escape fraction of ionizing photons also changes between 1-40% at z>10. The mass fraction of HII region changes with the star formation history, resulting in the fluctuations of metal lines and Lyman-alpha line luminosities. In the starbursting phase of galaxies with the halo mass $\sim 10^{11}\,{\rm M_{\odot}}$ ($10^{12}\,{\rm M_{\odot}}$), the simulated galaxy has $L_{\rm [O\,III]} \sim 10^{42}\, (10^{43})\, {\rm erg\,s^{-1}}$, which is consistent with the observed star-forming galaxies at z>7. Our simulations suggest that deep [C II] observation with ALMA can trace the distribution of neutral gas extending over $\sim 20$ physical kpc. We also find that the luminosity ratio $L_{\rm [O\,III]}/L_{\rm [C\,II]}$ decreases with bolometric luminosity due to metal enrichment. Our simulations show that the combination of multi-wavelength observations by ALMA and JWST will be able to reveal the multi-phase ISM structure and the transition from starbursting to outflowing phases of high-z galaxies.

preprint2015arXiv

Cold accretion in early galaxy formation and its Lyman-alpha signatures

The Lyman-alpha (Lya) emission has played an important role in detecting high-redshift galaxies, including recently distant ones at redshift z > 7. It may also contain important information on the origin of these galaxies. Here, we investigate the formation of a typical L* galaxy and its observational signatures at the earliest stage, by combining a cosmological hydrodynamic simulation with three-dimensional radiative transfer calculations using the newly improved ART^2 code. Our cosmological simulation uses the Aquila initial condition which zooms in onto a Milky Way-like halo with high resolutions, and our radiative transfer couples multi-wavelength continuum, Lya line, and ionization of hydrogen. We find that the modeled galaxy starts to form at redshift z ~ 24 through efficient accretion of cold gas, which produces a strong Lya line with a luminosity of L(Lya) ~ 10^42 erg/s as early as z ~ 14. The Lya emission appears to trace the cold, dense gas. The lines exhibit asymmetric, single-peak profiles, and are shifted to the blue wing, a characteristic feature of gas inflow. Moreover, the contribution to the total Lya luminosity by excitation cooling increases with redshift, and it becomes dominant at z >~ 6. We predict that L* galaxies such as the modeled one may be detected at z <~ 8 by JWST and ALMA with a reasonable integration time. Beyond redshift 12, however, only Lya line may be observable by spectroscopic surveys. Our results suggest that Lya line is one of the most powerful tools to detect the first generation of galaxies, and to decipher their formation mechanism.

preprint2014arXiv

Distinctive 21 cm structures of the first stars, galaxies, and quasars

Observations of the redshifted 21 cm line with upcoming radio telescopes promise to transform our understanding of the cosmic reionization. To unravel the underlying physical process, we investigate the 21 cm structures of three different ionizing sources, Population (Pop) III stars, the first galaxies, and the first quasars, by using radiative transfer simulations that include both ionization of neutral hydrogen and resonant scattering of Lya photons. We find that Pop III stars and quasars produce a smooth transition from an ionized and hot state to a neutral and cold one, owing to their hard spectral energy distribution with abundant ionizing photons, in contrast to the sharp transition in galaxies. Furthermore, Lya scattering plays a dominant role in producing the 21 cm signal as it determines the relation between hydrogen spin temperature and gas kinetic temperature. This effect, also called Wouthuysen-Field coupling, depends strongly on the ionizing source. It is the strongest around galaxies, where the spin temperature is highly coupled to that of the gas, resulting in extended absorption troughs in the 21 cm brightness temperature. On the other hand, in the case of Pop III stars, the 21 cm signal shows both emission and absorption regions around a small HII bubble. For quasars, a large emission region in the 21 cm signal is produced, and the absorption region decreases as the size of the HII bubble becomes large due to the limited traveling time of photons. We predict that future surveys from large radio arrays such as MWA, LOFAR and SKA may be able to detect the 21 cm signals of primordial galaxies and quasars, but not likely Pop III stars due to its small angular diameter.

preprint2014arXiv

Escape of Lyman-alpha and continuum photons from star-forming galaxies

A large number of high-redshift galaxies have been discovered via their narrow-band Lya line or broad-band continuum colors in recent years. The nature of the escaping process of photons from these early galaxies is crucial to understanding galaxy evolution and the cosmic reionization. Here, we investigate the escape of Lya, non-ionizing UV-continuum (l = 1300 - 1600 angstrom in rest frame), and ionizing photons (l < 912 angstrom) from galaxies by combining a cosmological hydrodynamic simulation with three-dimensional multi-wavelength radiative transfer calculations. The galaxies are simulated in a box of 5^3 h^-3 Mpc^3 with high resolutions using the Aquila initial condition which reproduces a Milky Way-like galaxy at redshift z=0. We find that the escape fraction (fesc) of these different photons shows a complex dependence on redshift and galaxy properties: fesc(Lya) and fesc(UV) appear to evolve with redshift, and they show similar, weak correlations with galaxy properties such as mass, star formation, metallicity, and dust content, while fesc(Ion) remains roughly constant at ~ 0.2 from z ~ 0 - 10, and it does not show clear dependence on galaxy properties. fesc(Lya) correlates more strongly with fesc(UV) than with fesc(Ion). In addition, we find a relation between the emergent Lya luminosity and the ionizing photon emissivity of Lyman Alpha Emitters (LAEs). By combining this relation with the observed luminosity functions of LAEs at different redshift, we estimate the contribution from LAEs to the reionization of intergalactic medium (IGM). Our result suggests that ionizing photons from LAEs alone are not sufficient to ionize IGM at z > 6, but they can maintain the ionization of IGM at z ~ 0 - 5.

preprint2014arXiv

Numerical Convergence in Smoothed Particle Hydrodynamics

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

preprint2014arXiv

Origin of the Correlations Between Supermassive Black Holes and Their Host Galaxies

Observations have shown that supermassive black holes in nearby elliptical galaxies correlate tightly with the stellar velocity dispersion (the $\MBH - σ$ relation) and the stellar mass (the $\MBH - \Mhost$ relation) of their host spheroids. However, the origin of these correlations remains ambiguous. In a previous paper by Zhu et al., we proposed a model which links the M-$σ$ relation to the the dynamical state of the system and the $\MBH - \Mhost$ relation to the self-regulation of galaxy growth. To test this model, we compile a sample of observed galaxies with different properties and examine the dependence of the above correlations on these parameters. We find that galaxies that satisfy the the $\MBH - σ$ correlation appear to have reached virial equilibrium, as indicated by the ratio between kinetic energy and gravitational potential, 2K/U $\sim$ 1. Furthermore, the ratio of black hole accretion rate to star formation rate remains nearly constant, BHAR /SFR $\sim$ $10^{-3}$, in active galaxies over a wide range of mass in the redshift range z=0 - 3. These results confirm our theoretical model that the observed correlations have different origins: the $\MBH - σ$ relation may result from galaxy relaxation, while the $\MBH$ - $\Mhost$ relation may be due to self-regulated black hole accretion and star formation in galaxies.

preprint2013arXiv

Effects of photon trapping on the Lyman-alpha properties of star-forming galaxies

Recent observations show that a large number of Lyman-alpha emitters (LAEs) at high redshift z >= 3 have unusually high Lya equivalent widths (EW > 400 angstrom). However, the origin of these high EWs is an open question. Here, we investigate the impacts of photon trapping on the Lya EW and other properties by tracking the Lya radiative transfer in spherical galactic clouds. We find that the delayed escape of the Lya photons can change the Lya properties significantly. During the transition phase from optically thick to optically thin where the Lya photons can escape simultaneously, the EW can be boosted to ~ 1000 angstrom, the Lya luminosity can be increased by a factor of a few, and the line profile can be significantly broadened. The boost factor appears to depend on the galaxy properties such as mass and star formation rate and timescale, therefore future investigation combing 3D Lya RT calculations with cosmological simulations of galaxy formation and evolution is needed to fully understand the Lya properties of early star-forming galaxies.

preprint2013arXiv

Extended Lyman-alpha emission from interacting galaxies at high redshifts

Recent observations have discovered a population of extended Lya sources, dubbed Lyman-alpha blobs (LABs), at high redshift z ~ 3 - 6.6. These LABs typically have a luminosity of L ~ 10^42-10^44 erg/s, and a size of tens of kiloparsecs, with some giant ones reaching up to D ~ 100 kpc. However, the origin of these LABs is not well understood. In this paper, we investigate a merger model for the formation of LABs by studying Lya emission from interacting galaxies at high redshifts by means of a combination of hydrodynamics simulations with three dimensional radiative transfer calculations. Our galaxy simulations focus on a set of binary major mergers of galaxies with a mass range of 3-7 *10^12 Msun in the redshift range of z ~ 3 -7, and we use the newly improved ART^2 code to perform the radiative transfer calculations which couple multi-wavelength continuum, ionization of hydrogen, and Lya line emission. We find that intense star formation and enhanced cooling induced by gravitational interaction produce strong Lya emission from these merging galaxies. The Lya emission appears to be extended due to the extended distribution of sources and gas. During the close encounter of galaxy progenitors when the star formation rate peaks at ~ 10^3 Msun/yr, our model produces Lya blobs with luminosity of L ~ 10^42-10^44 erg/s, and size of D ~ 10-20 kpc at z>6 and D ~ 20-50 kpc at z ~ 3, in broad agreement with observations in the same redshift range. Our results suggest that merging galaxies may produce some typical LABs as observed, but the giant ones may be produced by mergers more massive than those in our model, or a combination of mergers and cold accretion from filaments on a large scale.

preprint2013arXiv

The Effects of Local Primordial Non-Gaussianity on the Formation and Evolution of Galaxies

Thanks to the rapid progress in precision cosmology in the last few years, we now have access to physical observables that may constrain the theory of inflation through the non-Gaussianity (NG) signatures in the cosmic microwave background radiation and the distribution of large-scale structure. Numerical modeling of the NG signals from different inflation models is essential to correctly interpret current and near future data from large-scale structure surveys. In this study, we use high-resolution cosmological hydrodynamical simulations to investigate the effects of primordial NG on the formation and evolution of galaxies from the cosmic dawn to the present day. Focusing on the local type primordial NG, we find that it may affect the formation history of stars and black holes in galaxies, and their distribution. Compared to the Gaussian case, large non-Gaussian potential with $f_{NL} \gtrsim 10^3$ leads to earlier collapse of the first structures, more massive galaxies especially at high redshifts, stronger clustering of galaxies, and higher halo bias. However, for smaller NG with $f_{NL} \lesssim 10^2$, the effect is significantly weaker. Observations of the distribution and properties of high-redshift, rare objects such as the first galaxies and quasars may provide further constraints on the primordial NG.

preprint2012arXiv

ART^2 : Coupling Lyman-alpha Line and Multi-wavelength Continuum Radiative Transfer

Narrow-band Lya line and broad-band continuum have played important roles in the discovery of high-redshift galaxies in recent years. Hence, it is crucial to study the radiative transfer of both Lya and continuum photons in the context of galaxy formation and evolution in order to understand the nature of distant galaxies. Here, we present a three-dimensional Monte Carlo radiative transfer code, All-wavelength Radiative Transfer with Adaptive Refinement Tree (ART^2), which couples Lya line and multi-wavelength continuum, for the study of panchromatic properties of galaxies and interstellar medium. This code is based on the original version of Li et al., and features three essential modules: continuum emission from X-ray to radio, Lya emission from both recombination and collisional excitation, and ionization of neutral hydrogen. The coupling of these three modules, together with an adaptive refinement grid, enables a self-consistent and accurate calculation of the Lya properties. As an example, we apply ART^2 to a cosmological simulation that includes both star formation and black hole growth, and study in detail a sample of massive galaxies at redshifts z=3.1 - 10.2. We find that these galaxies are Lya emitters (LAEs), whose Lya emission traces the dense gas region, and that their Lya lines show a shape characteristic of gas inflow. Furthermore, the Lya properties, including photon escape fraction, emergent luminosity, and equivalent width, change with time and environment. Our results suggest that LAEs evolve with redshift, and that early LAEs such as the most distant one detected at z ~ 8.6 may be dwarf galaxies with a high star formation rate fueled by infall of cold gas, and a low Lya escape fraction.

preprint2012arXiv

Scaling Relations Between Low-mass Black Holes and Their Host Galaxies

It is well established that supermassive black holes in nearby elliptical galaxies correlate tightly with the kinematic property ($\mbhsigma$ correlation) and stellar mass ($\mbhhost$ correlation) of their host spheroids. However, it is not clear what the relations would be at the low-mass end, and how they evolve. Here, we investigate these relations in low-mass systems ($\MBH \sim \rm{10^{6}- 10^{8}}\, \Msun$) using the Aquila Simulation, a high-resolution cosmological hydrodynamic simulation which follows the formation and evolution of stars and black holes in a Milky Way-size galaxy and its substructures. We find a number of interesting results on the origin and evolution of the scaling relations in these systems: (1) there is a strong redshift evolution in the $\mbhsigma$ relation, but a much weaker one in the $\mbhhost$ relation; (2) there is a close link between the $\mbhsigma$ relation and the dynamical state of the system -- the galaxies that fall on the observed correlation appear to have reached virial equilibrium. (3) the star formation and black hole growth are self-regulated in galaxies -- the ratio between black hole accretion rate and star formation rate remains nearly constant in a wide redshift span $z = 0-6$. These findings suggest that the observed correlations have different origins: the $\mbhsigma$ relation may be the result of virial equilibrium, while the $\mbhhost$ relation may the result of self-regulated star formation and black hole growth in galaxies.

preprint2012arXiv

The Co-evolution of Cosmic Entropy and Structures in the Universe

According to the second law of thermodynamics, the arrow of time points to an ever increasing entropy of the Universe. However, exactly how the entropy evolves with time and what drives the growth remain largely unknown. Here, for the first time, we quantify the evolving entropy of cosmic structures using a large-scale cosmological hydrodynamical simulation. Our simulation starts from initial conditions predicted by the leading LambdaCDM cosmology, self-consistently evolves the dynamics of both dark and baryonic matter, star formation, black hole growth and feedback processes, from the cosmic dawn to the present day. Tracing the entropy contributions of these distinct components in the simulation, we find a strong link between entropy growth and structure formation. The entropy is dominated by that of the black holes in all epochs, and its evolution follows the same path as that of galaxies: it increases rapidly from a low-entropy state at high redshift until z~2, then transits to a slower growth. Our results suggest that cosmic entropy may co-evolve with cosmic structure, and that its growth may be driven mainly by the formation of black holes in galaxies. We predict that the entropy will continue to increase in the near future, but likely at a constant rate.

preprint2012arXiv

Were progenitors of local L* galaxies Lyman-alpha emitters at high redshift?

The Lya emission has been observed from galaxies over a redshift span z ~ 0 - 8.6. However, the evolution of high-redshift Lya emitters (LAEs), and the link between these populations and local galaxies, remain poorly understood. Here, we investigate the Lya properties of progenitors of a local L* galaxy by combining cosmological hydrodynamic simulations with three-dimensional radiative transfer calculations using the new ART^2 code. We find that the main progenitor (the most massive one) of a Milky Way-like galaxy has a number of Lya properties close to those of observed LAEs at z ~ 2 - 6, but most of the fainter ones appear to fall below the detection limits of current surveys. The Lya photon escape fraction depends sensitively on a number of physical properties of the galaxy, such as mass, star formation rate, and metallicity, as well as galaxy morphology and orientation. Moreover, we find that high-redshift LAEs show blue-shifted Lya line profiles characteristic of gas inflow, and that the Lya emission by excitation cooling increases with redshift, and becomes dominant at z > 6. Our results suggest that some observed LAEs at z ~ 2-6 with luminosity of L_Lya ~ 10^{42-43} ergs/s may be similar to the main progenitor of the Milky Way at high redshift, and that they may evolve into present-day L* galaxies.

preprint2011arXiv

Accretion onto Intermediate-mass Seed Black Holes in Primordial Galaxies

The origin of the supermassive black holes that power the most distant quasars observed is largely unknown. One hypothesis is that they grew rapidly from intermediate-mass seeds (~100 M_sun) left by the first stars. However, some previous studies argued that accretion onto these black holes was too low to build up the mass due to strong suppression by radiative feedback. Here, we re-exam the accretion process of such a black hole embedded in a primordial gas cloud, by considering a wide range of physical and numerical parameters not explored before. We find that, while radiative heating and pressure indeed suppress accretion effectively, self-gravity of the gas eventually overcomes the feedback effects and boosts the accretion to the Eddington rate after one free-fall timescale of the cloud. Moreover, for a given black hole mass, there exists a critical density above which the accretion can reach Eddington limit. Furthermore, we find a universal correlation between black hole accretion rate and ambient gas density, which may serve as a realistic recipe for black hole growth in simulations.

preprint2007arXiv

Correlations Between Central Massive Objects And Their Host Galaxies: From Bulgeless Spirals to Ellipticals

Recent observations by Ferrarese et al. (2006) and Wehner et al. (2006) reveal that a majority of galaxies contain a central massive object (CMO), either a supermassive black hole (SMBH) or a compact stellar nucleus, regardless of the galaxy mass or morphological type, and that there is a tight relation between the masses of CMOs and those of the host galaxies. Several recent studies show that feedback from black holes can successfully explain the $\msigma$ correlation in massive elliptical galaxies that contain SMBHs. However, puzzles remain in spirals or dwarf spheroids that do not appear to have black holes but instead harbor a compact central stellar cluster. Here we use three-dimensional, smoothed particle hydrodynamics simulations of isolated galaxies to study the formation and evolution of CMOs in bulgeless disk galaxies, and simulations of merging galaxies to study the transition of the CMO--host mass relation from late-type bulgeless spirals to early-type ellipticals. Our results suggest that the observed correlations may be established primarily by the depletion of gas in the central region by accretion and star-formation, and may hold for all galaxy types. A systematic search for CMOs in the nuclei of bulgeless disk galaxies would offer a test of this conclusion. (Abridged)

preprint2007arXiv

Photometric Properties of the Most Massive High-Redshift Galaxies

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

preprint2006arXiv

Star Formation in Isolated Disk Galaxies. II. Schmidt Laws and Efficiency of Gravitational Collapse

(Abridged). We model gravitational instability in a wide range of isolated disk galaxies, using GADGET, a three-dimensional, smoothed particle hydrodynamics code. The model galaxies include a dark matter halo and a disk of stars and isothermal gas. The global Schmidt law observed in disk galaxies is quantitatively reproduced by our models. We find that the surface density of star formation rate directly correlates with the strength of local gravitational instability. The local Schmidt laws of individual galaxies in our models show clear evidence of star formation thresholds. Our results suggest that the non-linear development of gravitational instability determines the local and global Schmidt laws, and the star formation thresholds.

preprint2006arXiv

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

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

preprint2005arXiv

The stellar mass spectrum from non-isothermal gravoturbulent fragmentation

Identifying the processes that determine the initial mass function of stars (IMF) is a fundamental problem in star formation theory. One of the major uncertainties is the exact chemical state of the star forming gas and its influence on the dynamical evolution. Most simulations of star forming clusters use an isothermal equation of state (EOS). However, theoretical predictions and observations suggest that the effective polytropic exponent gamma in the EOS varies with density. We address these issues and study the effect of a piecewise polytropic EOS on the formation of stellar clusters in turbulent, self-gravitating molecular clouds using three-dimensional, smoothed particle hydrodynamics simulations. To approximate the results of published predictions of the thermal behavior of collapsing clouds, we increase the polytropic exponent gamma from 0.7 to 1.1 at some chosen density n_c, which we vary. The change of thermodynamic state at n_c selects a characteristic mass scale for fragmentation M_ch, which we relate to the peak of the observed IMF. Our investigation generally supports the idea that the distribution of stellar masses depends mainly on the thermodynamic state of the star-forming gas. The thermodynamic state of interstellar gas is a result of the balance between heating and cooling processes, which in turn are determined by fundamental atomic and molecular physics and by chemical abundances. Given the abundances, the derivation of a characteristic stellar mass can thus be based on universal quantities and constants.

preprint2004arXiv

Effects of the Equation of State on the Formation of Star Clusters

We study the effect of varying the equation of state on the formation of stellar clusters in turbulent molecular clouds, using three-dimensional, smoothed particle hydrodynamics simulations. Our results show that the equation of state helps determine how strongly self-gravitating gas fragments. The degree of fragmentation decreases with increasing \new{polytropic exponent} $γ$ in the range $0.2 < γ< 1.4$, although the total amount of \new{mass accreted onto collapsed fragments} appears to remain roughly constant through that range. Low values of $γ$ are expected to lead to the formation of dense clusters of low-mass stars, while $γ>1$ probably results in the formation of isolated and massive stars. Fragmentation and collapse ceases entirely for $γ> 1.4$ as expected from analytic arguments. The mass spectrum of overdense gas clumps is roughly log-normal for {\em non}-self-gravitating turbulent gas, but changes to a power-law under the action of gravity. The spectrum of collapsed cores, on the other hand, remains log-normal for $γ\le 1$, but flattens markedly for $γ>1$. The density PDFs approach log-normal, with widths that decrease with increasing $γ$. Primordial gas may have effective $γ> 1$, in which case these results could help explain why models of the formation of the first stars tend to produce isolated, massive objects.