Search

Combining cosmological hydrodynamic simulations and radiative transfer (RT) calculations, we present predictions of multi-wavelength radiative properties of the first galaxies at z ∼ 6–5. We find that intermittent star formation due to supernova (SN) feedback causes the escape fraction of UV photons to fluctuate rapidly, which then produces the observed diversity of SEDs for high-z galaxies. The simulated galaxies make rapid transition between UV-bright and IR-bright phase, and our RT calculations suggest that dust temperatures in the first galaxies are higher than z < 3 galaxies with ∼ 60 K.

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., 5 pages, 4 figures, Proceedings IAU Symposium No. 352 "Uncovering early galaxy evolution in the ALMA and JWST era", Viana do Castelo, Portugal, June 2019

Over the past decades, rest-frame ultraviolet (UV) observations have provided large samples of UV luminous galaxies at redshift (z) greater than 6, during the so-called epoch of reionization. While a few of these UV identified galaxies revealed significant dust reservoirs, very heavily dust-obscured sources at these early times have remained elusive. They are limited to a rare population of extreme starburst galaxies, and companions of rare quasars. These studies conclude that the contribution of dust-obscured galaxies to the cosmic star formation rate density at $z>6$ is sub-dominant. Recent ALMA and Spitzer observations have identified a more abundant, less extreme population of obscured galaxies at $z=3-6$. However, this population has not been confirmed in the reionization epoch so far. Here, we report the discovery of two dust-obscured star forming galaxies at $z=6.6813\pm0.0005$ and $z=7.3521\pm0.0005$. These objects are not detected in existing rest-frame UV data, and were only discovered through their far-infrared [CII] lines and dust continuum emission as companions to typical UV-luminous galaxies at the same redshift. The two galaxies exhibit lower infrared luminosities and star-formation rates than extreme starbursts, in line with typical star-forming galaxies at $z\sim7$. This population of heavily dust-obscured galaxies appears to contribute 10-25 per cent to the $z>6$ cosmic star formation rate density., 16 pages, 7 figures, 4 tables, published in Nature

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$\alpha$ 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$\alpha$, 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., Comment: Accepted for publication in MNRAS; 20 pages, 26 figures

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., Comment: 24 pages, 19 figures, accepted for publication in MNRAS

The galaxy Malin 1 contains the largest stellar disk known but the formation mechanism of this structure has been elusive. In this paper, we report a Malin 1 analogue in the 100 Mpc IllustrisTNG simulation and describe its formation history. At redshift zero, this massive galaxy, having a maximum circular velocity $V_{\rm max}$ of 430 ${\rm km\ s^{-1}}$, contains a 100 kpc gas/stellar disk with morphology similar to Malin 1. The simulated galaxy reproduces well many observed features of Malin 1's vast disk, including its stellar ages, metallicities, and gas rotation curve. We trace the extended disk back in time and find that a large fraction of the cold gas at redshift zero originated from the cooling of hot halo gas, triggered by the merger of a pair of intruding galaxies. Our finding provides a novel way to form large galaxy disks as extreme as Malin 1 within the current galaxy formation framework., 5 page, 5 figures; resubmitted to MNRAS Letters after first referee report; Comments welcome. Animations of gas/stars can be found here at https://youtu.be/Ibl4Cyybw2Q and https://youtu.be/3aIpCu6e9dQ, or downloaded at http://www.tng-project.org/movies/tng/TNG_Malin_1_gas.mpeg and http://www.tng-project.org/movies/tng/TNG_Malin_1_star.mpeg

Recent observations have successfully detected UV-bright and infrared-bright galaxies in the epoch of reionization. However, the origin of their radiative properties has not been understood yet. Combining cosmological hydrodynamic simulations and radiative transfer calculations, we present predictions of multiwavelength radiative properties of the first galaxies at z ∼ 6–15. Using zoom-in initial conditions, we investigate three massive galaxies and their satellites in different environment and halo masses at z = 6: $M_{\rm h}= 2.4\times 10^{10}\,$, $1.6\times 10^{11}\, $, and $0.7\times 10^{12}\, {\rm M_{\odot }}$. We find that most of the gas and dust are ejected from star-forming regions by supernova feedback, which allows the UV photons to escape. We show that the peak of the spectral energy distribution (SED) rapidly changes between UV and infrared wavelengths on a time-scale of ∼ 100 Myr due to intermittent star formation and feedback, and the escape fraction of UV photons fluctuates in the range of 0.2–0.8 at z < 10 with a time-averaged value of 0.3. When dusty gas covers the star-forming regions, the galaxies become bright in the observed-frame sub-millimeter wavelengths. We predict the detectability of high-z galaxies with the Atacama Large Millimeter Array (ALMA). For a sensitivity limit of $0.1\, {\rm mJy}$ at $850\, {\rm \mu m}$, the detection probability of galaxies in haloes $M_{\rm h}\gtrsim 10^{11}\, \, {\rm M_{\odot }}$ at z ≲ 7 exceeds fifty per cent. We argue that supernova feedback can produce the observed diversity of SEDs for high-z galaxies.

The "Disk of satellites" (DoS) around Milky Way is a highly debated topic with conflicting interpretations of observations and their theoretical models. We perform a comprehensive analysis of all dwarfs detected in the Milky Way and find that the DoS structure depends strongly on the plane identification method and the sample size. In particular, we demonstrate that a smaller sample size produces a higher anisotropy of the spatial distribution and a stronger clustering of the angular momentum of the satellites. Moreover, we calculate the evolution of the 11 classical satellites with proper motion measurements and find that the thin DoS they currently reside in is transient. Furthermore, we analyze two cosmological simulations using the same initial conditions of a Milky Way-sized galaxy, an N-body run with dark matter only and a hydrodynamic one with both baryonic and dark matter, and find that the hydrodynamic simulation produces more anisotropic distribution of satellites than the N-body one. Our results suggest that an anisotropic distribution of satellites in galaxies can originate from baryonic processes in the hierarchical structure formation model, but the claimed highly-flattened, coherently-rotating DoS of the Milky Way may be a small-number selection effect. These findings may help resolve the contradictory claims of DoS in other galaxies and the discrepancy among numerical simulations., Comment: Published in ApJ. 7 pages, 5 figures

A large number of high-redshift galaxies have been discovered via their narrow-band Lyα line or broad-band continuum colours in recent years. The nature of the escaping process of photons from these early galaxies is crucial to understand galaxy evolution and the cosmic reionization. Here, we investigate the escape of Lyα, non-ionizing UV-continuum (λ = 1300-1600 A in rest frame), and ionizing photons (λ < 912 A) from galaxies by combining a cosmological hydrodynamic simulation with three-dimensional multiwavelength 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 (f_esc) of these different photons shows a complex dependence on redshift and galaxy properties: $f_{\rm esc}^{\rm \rm {Ly{\alpha }}}$ and $f_{\rm esc}^{\rm 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 $f_{\rm esc}^{\rm Ion}$ remains roughly constant at ∼0.2 from z ∼ 0 to 10, and it does not show clear dependence on galaxy properties. $f_{\rm esc}^{\rm \rm {Ly{\alpha }}}$ correlates more strongly with $f_{\rm esc}^{\rm UV}$ than with $f_{\rm esc}^{\rm Ion}$. In addition, we find a relation between the emergent Lyα 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.

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., 7 pages, 7 figures, accepted for publication in MNRAS

Observations of globular clusters show that they have universal lognormal mass functions with a characteristic peak at $\sim 2\times 10^{5}\, {\rm{M_{\odot}}}$, but the origin of this peaked distribution is highly debated. Here we investigate the formation and evolution of star clusters in interacting galaxies using high-resolution hydrodynamical simulations performed with two different codes in order to mitigate numerical artifacts. We find that massive star clusters in the range of $\sim 10^{5.5} - 10^{7.5}\, {\rm{M_{\odot}}}$ form preferentially in the highly-shocked regions produced by galaxy interactions. The nascent cluster-forming clouds have high gas pressures in the range of $P/k \sim 10^8 - 10^{12}\, \rm{K}\,\rm{cm^{-3}}$, which is $\sim 10^4 - 10^8$ times higher than the typical pressure of the interstellar medium but consistent with recent observations of a pre-super star cluster cloud in the Antennae Galaxies. Furthermore, these massive star clusters have quasi-lognormal initial mass functions with a peak around $\sim 10^{6}\, {\rm{M_{\odot}}}$. The number of clusters declines with time due to destructive processes, but the shape and the peak of the mass functions do not change significantly during the course of galaxy collisions. Our results suggest that gas-rich galaxy mergers may provide a favorable environment for the formation of massive star clusters such as globular clusters, and that the lognormal mass functions and the unique peak may originate from the extreme high-pressure conditions of the birth clouds and may survive the dynamical evolution., Comment: Accepted in ApJ. 16 pages, 9 figures

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.

We investigate the nature of CO emission from z~6 quasars by combining non-LTE radiative transfer calculations with merger-driven models of z~6 quasar formation that arise naturally in LCDM cosmological simulations. We consider four model quasars formed in 10^12-10^13 M_sun halos from different merging histories. Our main results follow. Owing to massive starbursts and funneling of dense gas into the nuclear regions of merging galaxies, the CO is highly excited and the flux density peaks between J=5-8. The CO morphology of z~6 quasars often exhibits multiple emission peaks which arise from H2 concentrations which have not yet fully coalesced. Quasars at z~6 display a large range of sightline dependent line widths such that the lines are narrowest when the rotating H2 gas associated with the quasar is viewed face-on (when L_B is largest), and broadest when the gas is seen edge-on (when L_B is lowest). Thus for all models selection effects exist such that quasars selected for optical luminosity are preferentially face-on which may result in detected CO line widths narrower than the median. The sightline averaged line width is reflective of the circular velocity (V_c) of the host halo, and ranges from sigma~300-650 km/s. For optically selected QSOs, 10-25% (halo-mass dependant) of sightlines have narrow line widths compatible with the sole CO detection at z>6, J1148+5251. When accounting for both the temporal evolution of CO line widths, as well as the redshift evolution of halo V_c, these models self-consistently account for the CO line widths of both z~2 sub-mm galaxies and QSO's. Finally, the dynamical mass derived from the sightline averaged line widths provides a good estimate of the total mass, and allows for a stellar bulge and SMBH consistent with the local M_BH-M_bulge relation. [abridged]

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., Version accepted by ApJ (5/21/2007), for high-resolution figures see http://kicp.uchicago.edu/~brant/astro-ph/massive_high_redshift_galaxies/

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.

We study gravitational instability and consequent star formation in a wide range of isolated disk galaxies, using three-dimensional, smoothed particle hydrodynamics simulations at resolution sufficient to fully resolve gravitational collapse. Stellar feedback is represented by an isothermal equation of state. Absorbing sink particles are inserted in dynamically bound, converging regions with number density $n > 10^3$ cm$^{-3}$ to directly measure the mass of gravitationally collapsing gas available for star formation. Our models quantitatively reproduce not only the observed Schmidt law, but also the observed star formation threshold in disk galaxies. Our results suggest that the dominant physical mechanism determining the star formation rate is just the strength of gravitational instability, with feedback primarily functioning to maintain a roughly constant effective sound speed., ApJL accepted, version with high quality color images can be found in http://research.amnh.org/~yuexing/astro-ph/0407247.pdf

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., 12 pages, 11 figures, accepted for publication in MNRAS

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, Comment: 10 pages, 11 figures, accepted for publication in ApJ

New far-infrared and sub-millimeter photometry from the Herschel Space Observatory is presented for 61 nearby galaxies from the Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel (KINGFISH) sample. The spatially-integrated fluxes are largely consistent with expectations based on Spitzer far-infrared photometry and extrapolations to longer wavelengths using popular dust emission models. Dwarf irregular galaxies are notable exceptions, as already noted by other authors, as their 500um emission shows evidence for a sub-millimeter excess. In addition, the fraction of dust heating attributed to intense radiation fields associated with photo-dissociation regions is found to be (21+/-4)% larger when Herschel data are included in the analysis. Dust masses obtained from the dust emission models of Draine & Li are found to be on average nearly a factor of two higher than those based on single-temperature modified blackbodies, as single blackbody curves do not capture the full range of dust temperatures inherent to any galaxy. The discrepancy is largest for galaxies exhibiting the coolest far-infrared colors., Accepted for publication in ApJ

The detection of large quantities of dust in z ~ 6 quasars by infrared and radio surveys presents puzzles for the formation and evolution of dust in these early systems. Previously (Li et al. 2007), we showed that luminous quasars at z > 6 can form through hierarchical mergers of gas-rich galaxies. Here, we calculate the dust properties of simulated quasars and their progenitors using a three-dimensional Monte Carlo radiative transfer code, ART^2 -- All-wavelength Radiative Transfer with Adaptive Refinement Tree. ART^2 incorporates a radiative equilibrium algorithm for dust emission, an adaptive grid for inhomogeneous density, a multiphase model for the ISM, and a supernova-origin dust model. We reproduce the SED and dust properties of SDSS J1148+5251, and find that the infrared emission are closely associated with the formation and evolution of the quasar host. The system evolves from a cold to a warm ULIRG owing to heating and feedback from stars and AGN. Furthermore, the AGN has significant implications for the interpretation of observation of the hosts. Our results suggest that vigorous star formation in merging progenitors is necessary to reproduce the observed dust properties of z~6 quasars, supporting a merger-driven origin for luminous quasars at high redshifts and the starburst-to-quasar evolutionary hypothesis. (Abridged), 26 pages, 22 figures, accepted by ApJ. Version with full resolution images is available at http://www.cfa.harvard.edu/~yxli/ARTDUST/astroph0706.3706.pdf

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), 11 pages, 8 figures, accepted to ApJ

We model gravitational collapse leading to star formation in a wide range of isolated disk galaxies using a three-dimensional, smoothed particle hydrodynamics code. The model galaxies include a dark matter halo and a disk of stars and isothermal gas. Absorbing sink particles are used to directly measure the mass of gravitationally collapsing gas. They reach masses characteristic of stellar clusters. In this paper, we describe our galaxy models and numerical methods, followed by an investigation of the gravitational instability in these galaxies. Gravitational collapse forms star clusters with correlated positions and ages, as observed, for example, in the Large Magellanic Cloud. Gravitational instability alone acting in unperturbed galaxies appears sufficient to produce flocculent spiral arms, though not more organized patterns. Unstable galaxies show collapse in thin layers in the galactic plane; associated dust will form thin dust lanes in those galaxies, in agreement with observations. (abridged), Comment: 49 pages, 22 figures, to appear in ApJ (July, 2005), version with high quality color images can be fond in http://research.amnh.org/~yuexing/astro-ph/0501022.pdf

(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., Comment: 21 pages, 15 figures, to appear on ApJ March 2006. (Paper I is ApJ 626, 823, 2005, or astro-ph/0501022)

We present a high-resolution simulation of globular cluster formation in a galaxy merger. For the first time in such a simulation, individual star clusters are directly identified and followed on their orbits. We quantitatively compare star formation in the merger to that in the unperturbed galaxies. The merging galaxies show a strong starburst, in sharp contrast to their isolated progenitors. Most star clusters form in the tidal features. With a mass range of $5\times10^{5}$--$5\times 10^{6} M_{\odot}$, they are identified as globular clusters. The merger remnant is an elliptical galaxy. Clusters with different mass or age have different radial distributions in the galaxy. Our results show that the high specific frequency and bimodal distribution of metallicity observed in elliptical galaxies are natural products of gas-rich mergers, supporting a merger origin for the ellipticals and their globular cluster systems., Comment: ApJL accepted, version with high quality color images can be found in http://research.amnh.org/~yuexing/astro-ph/0407248.pdf

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.

The Ly{alpha} emission has been observed from galaxies over a redshift span z {approx} 0-8.6. However, the evolution of high-redshift Ly{alpha} emitters (LAEs), and the link between these populations and local galaxies, remains poorly understood. Here, we investigate the Ly{alpha} properties of progenitors of a local L* galaxy by combining cosmological hydrodynamic simulations with three-dimensional radiative transfer calculations using the new ART{sup 2} code. We find that the main progenitor (the most massive one) of a Milky-Way-like galaxy has a number of Ly{alpha} properties close to those of observed LAEs at z {approx} 2-6, but most of the fainter ones appear to fall below the detection limits of current surveys. The Ly{alpha} 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 blueshifted Ly{alpha} line profiles characteristic of gas inflow, and that the Ly{alpha} emission by excitation cooling increases with redshift, and becomes dominant at z {approx}> 6. Our results suggest that some observed LAEs at z {approx} 2-6 with luminosity of L{sub Ly{alpha}} {approx} 10{sup 42}-10{sup 43} erg s{sup -1} may be similarmore » to the main progenitor of the Milky Way at high redshift, and that they may evolve into present-day L* galaxies.« less