Person:

Hernquist, Lars

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

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

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

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

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

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. showed that luminous quasars at zgsim 6 can form through hierarchical mergers of gas-rich galaxies, and that these systems are expected to evolve from starburst through quasar phases. Here, we calculate the dust properties of simulated quasars and their progenitors using a three-dimensional Monte Carlo radiative transfer code, ART2 (All-wavelength Radiative Transfer with Adaptive Refinement Tree). ART2 incorporates a radiative equilibrium algorithm which treats dust emission self-consistently, an adaptive grid method which can efficiently cover a large dynamic range in both spatial and density scales, a multiphase model of the interstellar medium which accounts for the observed scaling relations of molecular clouds, and a supernova-origin model for dust which can explain the existence of dust in cosmologically young objects. By applying ART2 to the hydrodynamic simulations of Li et al., we reproduce the observed spectral energy distribution (SED) and inferred dust properties of SDSS J1148+5251, the most distant Sloan quasar. We find that the dust and infrared emission are closely associated with the formation and evolution of the quasar host. The system evolves from a cold to a warm ultraluminous infrared galaxy (ULIRG) owing to heating and feedback from stars and the active galactic nucleus (AGN). Furthermore, the AGN activity 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.

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

We present a constrained transport (CT) algorithm for solving the 3D ideal magnetohydrodynamic (MHD) equations on a moving mesh, which maintains the divergence-free condition on the magnetic field to machine-precision. Our CT scheme uses an unstructured representation of the magnetic vector potential, making the numerical method simple and computationally efficient. The scheme is implemented in the moving mesh code arepo. We demonstrate the performance of the approach with simulations of driven MHD turbulence, a magnetized disc galaxy, and a cosmological volume with primordial magnetic field. We compare the outcomes of these experiments to those obtained with a previously implemented Powell divergence-cleaning scheme. While CT and the Powell technique yield similar results in idealized test problems, some differences are seen in situations more representative of astrophysical flows. In the turbulence simulations, the Powell cleaning scheme artificially grows the mean magnetic field, while CT maintains this conserved quantity of ideal MHD. In the disc simulation, CT gives slower magnetic field growth rate and saturates to equipartition between the turbulent kinetic energy and magnetic energy, whereas Powell cleaning produces a dynamically dominant magnetic field. Such difference has been observed in adaptive-mesh refinement codes with CT and smoothed-particle hydrodynamics codes with divergence-cleaning. In the cosmological simulation, both approaches give similar magnetic amplification, but Powell exhibits more cell-level noise. CT methods in general are more accurate than divergence-cleaning techniques, and, when coupled to a moving mesh can exploit the advantages of automatic spatial/temporal adaptivity and reduced advection errors, allowing for improved astrophysical MHD simulations.

We present the first systematic comparison of ultraviolet-millimeter spectral energy distributions (SEDs) of observed and simulated interacting galaxies. Our sample is drawn from the Spitzer Interacting Galaxy Survey and probes a range of galaxy interaction parameters. We use 31 galaxies in 14 systems which have been observed with Herschel, Spitzer, GALEX, and 2MASS. We create a suite of GADGET-3 hydrodynamic simulations of isolated and interacting galaxies with stellar masses comparable to those in our sample of interacting galaxies. Photometry for the simulated systems is then calculated with the SUNRISE radiative transfer code for comparison with the observed systems. For most of the observed systems, one or more of the simulated SEDs match reasonably well. The best matches recover the infrared luminosity and the star formation rate of the observed systems, and the more massive systems preferentially match SEDs from simulations of more massive galaxies. The most morphologically distorted systems in our sample are best matched to the simulated SEDs that are close to coalescence, while less evolved systems match well with the SEDs over a wide range of interaction stages, suggesting that an SED alone is insufficient for identifying the interaction stage except during the most active phases in strongly interacting systems. This result is supported by our finding that the SEDs calculated for simulated systems vary little over the interaction sequence.

The formation and early evolution of circumstellar discs often occurs within dense, newborn stellar clusters. For the first time, we apply the moving-mesh code arepo, to circumstellar discs in 3D, focusing on disc–disc interactions that result from stellar flybys. Although a small fraction of stars are expected to undergo close approaches, the outcomes of the most violent encounters might leave an imprint on the discs and host stars that will influence both their orbits and their ability to form planets. We first construct well-behaved 3D models of self-gravitating discs, and then create a suite of numerical experiments of parabolic encounters, exploring the effects of pericentre separation rp, disc orientation and disc–star mass ratio (Md/M*) on the orbital evolution of the host stars. Close encounters (2rp ≲ disc radius) can truncate discs on very short time-scales. If discs are massive, close encounters facilitate enough orbital angular momentum extraction to induce stellar capture. We find that for realistic primordial disc masses Md ≲ 0.1M*, non-colliding encounters induce minor orbital changes, which is consistent with analytic calculations of encounters in the linear regime. The same disc masses produce entirely different results for grazing/colliding encounters. In the latter case, rapidly cooling discs lose orbital energy by radiating away the energy excess of the shock-heated gas, thus causing capture of the host stars into a bound orbit. In rare cases, a tight binary with a circumbinary disc forms as a result of this encounter.