Robert Feldmann

THE X-FACTOR IN GALAXIES. I. DEPENDENCE ON ENVIRONMENT AND SCALE

Characterizing the conversion factor between CO emission and column density of molecular hydrogen, X CO, is crucial in studying the gaseous content of galaxies, its evolution, and relation to star formation. In most cases the conversion factor is assumed to be close to that of giant molecular clouds (GMCs) in the Milky Way, except possibly for mergers and star-bursting galaxies. However, there are physical grounds to expect that it should also depend on the gas metallicity, surface density, and strength of the interstellar radiation field. The X CO factor may also depend on the scale on which CO emission is averaged due to effects of limited resolution. We study the dependence of X CO on gas properties and averaging scale using a model that is based on a combination of results of sub-parsec scale magnetohydrodynamic simulations and on the gas distribution from self-consistent cosmological simulations of galaxy formation. Our model predicts X CO (2-4) × 1020 K–1 cm–2 km–1 s, consistent with the Galactic value, for interstellar medium conditions typical for the Milky Way. For such conditions the predicted X CO varies by only a factor of two for gas surface densities in the range . However, the model also predicts that more generally on the scale of GMCs, X CO is a strong function of metallicity and depends on the column density and the interstellar UV flux. We show explicitly that neglecting these dependencies in observational estimates can strongly bias the inferred distribution of H2 column densities of molecular clouds to have a narrower and offset range compared to the true distribution. We find that when averaged on ~kiloparsec scales the X-factor depends only weakly on radiation field and column density, but is still a strong function of metallicity. The predicted metallicity dependence can be approximated as X COZ –γ with γ 0.5-0.8.

The formation of submillimetre-bright galaxies from gas infall over a billion years

Submillimetre-bright galaxies at high redshift are the most luminous, heavily star-forming galaxies in the Universe and are characterized by prodigious emission in the far-infrared, with a flux of at least five millijanskys at a wavelength of 850 micrometres. They reside in haloes with masses about 1013 times that of the Sun, have low gas fractions compared to main-sequence disks at a comparable redshift, trace complex environments and are not easily observable at optical wavelengths. Their physical origin remains unclear. Simulations have been able to form galaxies with the requisite luminosities, but have otherwise been unable to simultaneously match the stellar masses, star formation rates, gas fractions and environments. Here we report a cosmological hydrodynamic galaxy formation simulation that is able to form a submillimetre galaxy that simultaneously satisfies the broad range of observed physical constraints. We find that groups of galaxies residing in massive dark matter haloes have increasing rates of star formation that peak at collective rates of about 500–1,000 solar masses per year at redshifts of two to three, by which time the interstellar medium is sufficiently enriched with metals that the region may be observed as a submillimetre-selected system. The intense star formation rates are fuelled in part by the infall of a reservoir gas supply enabled by stellar feedback at earlier times, not through major mergers. With a lifetime of nearly a billion years, our simulations show that the submillimetre-bright phase of high-redshift galaxies is prolonged and associated with significant mass buildup in early-Universe proto-clusters, and that many submillimetre-bright galaxies are composed of numerous unresolved components (for which there is some observational evidence).

THE HUBBLE SEQUENCE IN GROUPS: THE BIRTH OF THE EARLY-TYPE GALAXIES

The physical mechanisms and timescales that determine the morphological signatures and the quenching of star formation of typical (~L*) elliptical galaxies are not well understood. To address this issue, we have simulated the formation of a group of galaxies with sufficient resolution to track the evolution of gas and stars inside about a dozen galaxy group members over cosmic history. Galaxy groups, which harbor many elliptical galaxies in the universe, are a particularly promising environment to investigate morphological transformation and star formation quenching, due to their high galaxy density, their relatively low velocity dispersion, and the presence of a hot intragroup medium. Our simulation reproduces galaxies with different Hubble morphologies and, consequently, enables us to study when and where the morphological transformation of galaxies takes place. The simulation does not include feedback from active galactic nuclei showing that it is not an essential ingredient for producing quiescent, red elliptical galaxies in galaxy groups. Ellipticals form, as suspected, through galaxy mergers. In contrast with what has often been speculated, however, these mergers occur at z>1, before the merging progenitors enter the virial radius of the group and before the group is fully assembled. The simulation also shows that quenching of star formation in the still star-forming elliptical galaxies lags behind their morphological transformation, but, once started, is taking less than a billion years to complete. As long envisaged the star formation quenching happens as the galaxies approach and enter the finally assembled group, due to quenching of gas accretion and (to a lesser degree) stripping. A similar sort is followed by unmerged, disk galaxies, which, as they join the group, are turned into the red-and-dead disks that abound in these environments.

FIRE-2 Simulations: Physics versus Numerics in Galaxy Formation

The Feedback In Realistic Environments (FIRE) project explores feedback in cosmological galaxy formation simulations. Previous FIRE simulations used an identical source code (“FIRE-1”) for consistency. Motivated by the development of more accurate numerics – including hydrodynamic solvers, gravitational softening, and supernova coupling algorithms – and exploration of new physics (e.g. magnetic fields), we introduce “FIRE-2”, an updated numerical implementation of FIRE physics for the GIZMO code. We run a suite of simulations and compare against FIRE-1: overall, FIRE-2 improvements do not qualitatively change galaxy-scale properties. We pursue an extensive study of numerics versus physics. Details of the star-formation algorithm, cooling physics, and chemistry have weak effects, provided that we include metal-line cooling and star formation occurs at higher-than-mean densities. We present new resolution criteria for high-resolution galaxy simulations. Most galaxy-scale properties are robust to numerics we test, provided: (1) Toomre masses are resolved; (2) feedback coupling ensures conservation, and (3) individual supernovae are time-resolved. Stellar masses and profiles are most robust to resolution, followed by metal abundances and morphologies, followed by properties of winds and circum-galactic media (CGM). Central (∼kpc) mass concentrations in massive (>L*) galaxies are sensitive to numerics (via trapping/recycling of winds in hot halos). Multiple feedback mechanisms play key roles: supernovae regulate stellar masses/winds; stellar mass-loss fuels late star formation; radiative feedback suppresses accretion onto dwarfs and instantaneous star formation in disks. We provide all initial conditions and numerical algorithms used.

The Argo simulation – I. Quenching of massive galaxies at high redshift as a result of cosmological starvation

Observations show a prevalence of high redshift galaxies with large stellar masses and predominantly passive stellar populations. A variety of processes have been suggested that could reduce the star formation in such galaxies to observed levels, including quasar mode feedback, virial shock heating, or galactic winds driven by stellar feedback. However, the main quenching mechanisms have yet to be identified. Here we study the origin of star formation quenching using Argo, a cosmological zoom-in simulation that follows the evolution of a massive galaxy at z > 2. This simulation adopts the same sub-grid recipes of the Eris simulations, which have been shown to form realistic disk galaxies, and, in one version, adopts also a mass and spatial resolution identical to Eris. The resulting galaxy has properties consistent with those of observed, massive (M 10 11 M ) galaxies at z 2 and with abundance matching predictions. Our models do not include AGN feedback indicating that supermassive black holes likely play a subordinate role in determining masses and sizes of massive galaxies at high z. The specific star formation rate (sSFR) of the simulated galaxy matches the observed M - sSFR relation at early times. This period of smooth stellar mass growth comes to a sudden halt at z = 3:5 when the sSFR drops by almost an order of magnitude within a few hundred Myr. The suppression is initiated by a leveling off and a subsequent reduction of the cool gas accretion rate onto the galaxy, and not by feedback processes. This “cosmological starvation” occurs as the parent dark matter halo switches from a fast collapsing mode to a slow accretion mode. Additional mechanisms, such as perhaps radio mode feedback from an AGN, are needed to quench any residual star formation of the galaxy and to maintain a low sSFR until the present time.

(Star)bursts of FIRE: observational signatures of bursty star formation in galaxies

Galaxy formation models are now able to reproduce observed relations such as the relation between galaxies’ star formation rates (SFRs) and stellar masses (M^*) and the stellar-mass–halo-mass relation. We demonstrate that comparisons of the short-time-scale variability in galaxy SFRs with observational data provide an additional useful constraint on the physics of galaxy formation feedback. We apply SFR indicators with different sensitivity time-scales to galaxies from the Feedback in Realistic Environments (FIRE) simulations. We find that the SFR–M* relation has a significantly greater scatter when the Hα-derived SFR is considered compared with when the far-ultraviolet (FUV)-based SFR is used. This difference is a direct consequence of bursty star formation because the FIRE galaxies exhibit order-of-magnitude SFR variations over time-scales of a few Myr. We show that the difference in the scatter between the simulated Hα- and FUV-derived SFR–M^* relations at z = 2 is consistent with observational constraints. We also find that the Hα/FUV ratios predicted by the simulations at z = 0 are similar to those observed for local galaxies except for a population of low-mass (M* ≲ 10^(9.5) M_⊙) simulated galaxies with lower Hα/FUV ratios than observed. We suggest that future cosmological simulations should compare the Hα/FUV ratios of their galaxies with observations to constrain the feedback models employed.

INTRINSIC ALIGNMENT OF CLUSTER GALAXIES: THE REDSHIFT EVOLUTION

We present measurements of two types of cluster galaxy alignments based on a volume limited and highly pure (≥90%) sample of clusters from the GMBCG catalog derived from Data Release 7 of the Sloan Digital Sky Survey (SDSS DR7). We detect a clear brightest cluster galaxy (BCG) alignment (the alignment of major axis of the BCG toward the distribution of cluster satellite galaxies). We find that the BCG alignment signal becomes stronger as the redshift and BCG absolute magnitude decrease and becomes weaker as BCG stellar mass decreases. No dependence of the BCG alignment on cluster richness is found. We can detect a statistically significant (≥3σ) satellite alignment (the alignment of the major axes of the cluster satellite galaxies toward the BCG) only when we use the isophotal fit position angles (P.A.s), and the satellite alignment depends on the apparent magnitudes rather than the absolute magnitudes of the BCGs. This suggests that the detected satellite alignment based on isophotal P.A.s from the SDSS pipeline is possibly due to the contamination from the diffuse light of nearby BCGs. We caution that this should not be simply interpreted as non-existence of the satellite alignment, but rather that we cannot detect them with our current photometric SDSS data. We perform our measurements on both SDSS r-band and i-band data, but do not observe a passband dependence of the alignments.

A stellar feedback origin for neutral hydrogen in high-redshift quasar-mass haloes

Observations of quasar pairs reveal that quasar host halos at z ~ 2 have large covering fractions of cool dense gas (≳ 60% for Lyman limit systems within a projected virial radius). Most simulations have so far failed to explain these large observed covering fractions. We analyze a new set of 15 simulated massive halos with explicit stellar feedback from the FIRE project, covering the halo mass range M_h ≈ 2 x 10^(12) - 10^(13) M_☉ at z = 2. This extends our previous analysis of the circum-galactic medium of high-redshift galaxies to more massive halos. Feedback from active galactic nuclei (AGN) is not included in these simulations. We find covering fractions consistent with those observed around z ~ 2 quasars. The large HI covering fractions arise from star formation-driven galactic winds, including winds from low-mass satellite galaxies that interact with the cosmological infalling filaments in which they are typically embedded. The simulated covering fractions increase with both halo mass and redshift over the ranges covered, as well as with resolution. Our simulations predict that galaxies occupying dark matter halos of mass similar to quasars but without a luminous AGN should have Lyman limit system covering fractions comparable to quasars. This prediction can be tested by measuring covering fractions transverse to sub-millimeter galaxies or to more quiescent galaxies selected based on their high stellar mass.

The Argo Simulation: II. The Early Build-up of the Hubble Sequence

The Hubble sequence is a common classification scheme for the structure of galaxies. Despite the tremendous usefulness of this diagnostic, we still do not fully understand when, where, and how this morphological ordering was put in place. Here, we investigate the morphological evolution of a sample of 22 high redshift (z > 3) galaxies extracted from the Argo simulation. Argo is a cosmological zoom-in simulation of a group-sized halo and its environment. It adopts the same high resolution ( 10 4 M , 100 pc) and sub-grid physical model that was used in the Eris simulation but probes a sub-volume almost ten times bigger with as many as 45 million gas and star particles in the zoom-in region. Argo follows the early assembly of galaxies with a broad range of stellar masses (logM?=M 8 11 at z’ 3), while resolving properly their structural properties. We recover a diversity of morphologies, including latetype/irregular disc galaxies with flat rotation curves, spheroid dominated early-type discs, and a massive elliptical galaxy, already established at z 3. We identify major mergers as the main trigger for the formation of bulges and the steepening of the circular velocity curves. Minor mergers and non-axisymmetric perturbations (stellar bars) drive the bulge growth in some cases. The specific angular momenta of the simulated disc components fairly match the values inferred from nearby galaxies of similar M? once the expected redshift evolution of disc sizes is accounted for. We conclude that morphological transformations of high redshift galaxies of intermediate mass are likely triggered by processes similar to those at low redshift and result in an early build-up of the Hubble sequence.

The equilibrium view on dust and metals in galaxies: Galactic outflows drive low dust-to-metal ratios in dwarf galaxies

Most galaxy evolution simulations as well as a variety of observational methods assume a linear scaling between the (galaxy-averaged) dust-to-gas ratio D and metallicity Z of the interstellar medium (ISM). Indeed, nearby galaxies with solar or moderately sub-solar metallicities clearly follow this trend albeit with significant scatter. However, a growing number of observations show that the linear scaling breaks down for metal-poor galaxies (Z<0.2 Z_sun), highlighting the need for a more sophisticated modeling of the dust-to-metal ratio of galaxies. Here we study the co-evolution of dust and metal abundances in galaxies with the help of a dynamical, one-zone model that incorporates dust formation and destruction processes in addition to gas inflows, outflows, and metal enrichment. The dynamical model is consistent with various observational constraints, including the stellar mass -- metallicity relation, the stellar mass -- halo mass relation, and the observed Z -- D relation for both metal-poor and metal-rich galaxies. The functional form of the Z -- D relation follows from a basic equilibrium ansatz, similar to the ideas used previously to model the stellar mass -- metallicity relation. Galactic outflows regulate the inflow rate of gas from the cosmic web for galaxies of a given star formation rate. The mass loading factor of outflows thus dictates the rate at which the dust and metal content of the ISM is diluted. The stellar mass dependence of the mass loading factor drives the evolution of metallicities, dust-to-gas ratios, and dust-to-metal ratios in galaxies.