We analyze the different fates of the circumgalactic medium (CGM) in FIRE-2 cosmological simulations, focusing on the redshifts z=0.25 and z=2 representative of recent surveys. Our analysis includes 21 zoom-in simulations covering the halo mass range Mh(z=0) ~ 10^10 - 10^12 Msun. We analyze both where the gas ends up after first leaving the CGM (its "proximate" fate), as well as its location at z=0 (its "ultimate" fate). Of the CGM at z=2, about half is found in the ISM or stars of the central galaxy by z=0 in Mh(z=2) ~ 5e11 Msun halos, but most of the CGM in lower-mass halos ends up in the IGM. This is so even though most of the CGM in M_h(z=2) ~ 5e10 Msun halos first accretes onto the central galaxy before being ejected into the IGM. On the other hand, most of the CGM mass at z=0.25 remains in the CGM by z=0 at all halo masses analyzed. Of the CGM gas that subsequently accretes onto the central galaxy in the progenitors of Mh(z=0) ~10^12 Msun halos, most of it is cool (T~10^4 K) at z=2 but hot (~Tvir) at z=0.25, consistent with the expected transition from cold mode to hot mode accretion. Despite the transition in accretion mode, at both z=0.25 and z=2 >~80% of the cool gas in Mh >~ 10^11 Msun halos will accrete onto a galaxy. We find that the metallicity of CGM gas is typically a poor predictor of both its proximate and ultimate fates. This is because there is in general little correlation between the origin of CGM gas and its fate owing to substantial mixing while in the CGM.

We present a framework for characterizing the large scale movement of baryons relative to dark matter in cosmological simulations, requiring only the initial conditions and final state of the simulation. This is performed using the spread metric which quantifies the distance in the final conditions between initially neighbouring particles, and by analysing the baryonic content of final haloes relative to that of the initial Lagrangian regions defined by their dark matter component. Applying this framework to the SIMBA cosmological simulations, we show that 40% (10%) of cosmological baryons have moved $> 1h^{-1}~ {\rm Mpc}^{-1}$ ($3h^{-1}~ {\rm Mpc}^{-1}$) by $z=0$, due primarily to entrainment of gas by jets powered by AGN, with baryons moving up to $12h^{-1}~ {\rm Mpc}^{-1}$ away in extreme cases. Baryons decouple from the dynamics of the dark matter component due to hydrodynamic forces, radiative cooling, and feedback processes. As a result, only 60% of the gas content in a given halo at $z=0$ originates from its Lagrangian region, roughly independent of halo mass. A typical halo in the mass range $M_{\rm vir} = 10^{12}$--$10^{13}{\rm M}_\odot$ only retains 20% of the gas originally contained in its Lagrangian region. We show that up to 20% of the gas content in a typical Milky Way mass halo may originate in the region defined by the dark matter of another halo. This inter-Lagrangian baryon transfer may have important implications for the origin of gas and metals in the circumgalactic medium of galaxies, as well as for semi-analytic models of galaxy formation and "zoom-in" simulations.

We use a particle tracking analysis to study the origins of the circumgalactic medium (CGM), separating it into (1) accretion from the intergalactic medium (IGM), (2) wind from the central galaxy, and (3) gas ejected from other galaxies. Our sample consists of 21 FIRE-2 simulations, spanning the halo mass range log(Mh/Msun) ~ 10-12 , and we focus on z=0.25 and z=2. Owing to strong stellar feedback, only ~L* halos retain a baryon mass >~50% of their cosmic budget. Metals are more efficiently retained by halos, with a retention fraction >~50%. Across all masses and redshifts analyzed >~60% of the CGM mass originates as IGM accretion (some of which is associated with infalling halos). Overall, the second most important contribution is wind from the central galaxy, though gas ejected or stripped from satellites can contribute a comparable mass in ~L* halos. Gas can persist in the CGM for billions of years, resulting in well-mixed halo gas. Sight lines through the CGM are therefore likely to intersect gas of multiple origins. For low-redshift ~L* halos, cool gas (T<10^4.7 K) is distributed on average preferentially along the galaxy plane, however with strong halo-to-halo variability. The metallicity of IGM accretion is systematically lower than the metallicity of winds (typically by >~1 dex), although CGM and IGM metallicities depend significantly on the treatment of subgrid metal diffusion. Our results highlight the multiple physical mechanisms that contribute to the CGM and will inform observational efforts to develop a cohesive picture.

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.

We use cosmological hydrodynamic simulations with stellar feedback from the FIRE project to study the physical nature of Lyman limit systems (LLSs) at z<1. At these low redshifts, LLSs are closely associated with dense gas structures surrounding galaxies, such as galactic winds, dwarf satellites, and cool inflows from the intergalactic medium. Our analysis is based on 14 zoom-in simulations covering the halo mass range M_h~10^9-10^13 Msun at z=0, which we convolve with the dark matter halo mass function to produce cosmological statistics. We find that the majority of cosmologically-selected LLSs are associated with halos in the mass range 10^10 < M_h < 10^12 Msun. The incidence and HI column density distribution of simulated absorbers with columns 10^16.2 < N_HI < 2x10^20 cm^-2 are consistent with observations. High-velocity outflows (with radial velocity exceeding the halo circular velocity by a factor >~2) tend to have higher metallicities ([X/H] ~ -0.5) while very low metallicity ([X/H] < -2) LLSs are typically associated with gas infalling from the intergalactic medium. However, most LLSs occupy an intermediate region in metallicity-radial velocity space, for which there is no clear trend between metallicity and radial kinematics. Metal-enriched inflows arise in the FIRE simulations as a result of galactic winds that fall back onto galaxies at low redshift. The overall simulated LLS metallicity distribution has a mean (standard deviation) [X/H] = -0.9 (0.4) and does not show significant evidence for bimodality, in contrast to recent observational studies but consistent with LLSs arising from halos with a broad range of masses and metallicities.

We present an analysis of the flow of metals through the circumgalactic medium (CGM) in the Feedback in Realistic Environments (FIRE) simulations of galaxy formation, ranging from isolated dwarfs to $L*$ galaxies. We find that nearly all metals produced in high-redshift galaxies are carried out in winds that reach $0.25 R_{\rm vir}$. When measured at $0.25 R_{\rm vir}$ the metallicity of outflows is slightly higher than the interstellar medium (ISM) metallicity. Many metals thus reside in the CGM. Cooling and recycling from this reservoir determine the metal budget in the ISM. The outflowing metal flux decreases by a factor of $\sim2-5$ between $0.25 R_{\rm vir}$ and $R_{\rm vir}$. Furthermore, outflow metallicity is typically lower at $R_{\rm vir}$ owing to dilution of the remaining outflow by metal-poor material swept up from the CGM. The inflow metallicity at $R_{\rm vir}$ is generally low, but outflow and inflow metallicities are similar in the inner halo. At low redshift, massive galaxies no longer generate outflows that reach the CGM, causing a divergence in CGM and ISM metallicity. Dwarf galaxies continue to generate outflows, although they preferentially retain metal ejecta. In all but the least massive galaxy considered, a majority of the metals are within the halo at $z=0$. We measure the fraction of metals in CGM, ISM and stars, and quantify the thermal state of CGM metals in each halo. The total amount of metals in the low-redshift CGM of two simulated $L*$ galaxies is consistent with estimates from the COS halos survey, while for the other two it appears to be lower.

We study the structure, age and metallicity gradients, and dynamical evolution using a cosmological zoom-in simulation of a Milky Way-mass galaxy from the Feedback in Realistic Environments project. In the simulation, stars older than 6 Gyr were formed in a chaotic, bursty mode and have the largest vertical scale heights (1.5-2.5 kpc) by z=0, while stars younger than 6 Gyr were formed in a relatively calm, stable disk. The vertical scale height increases with stellar age at all radii, because (1) stars that formed earlier were thicker "at birth", and (2) stars were kinematically heated to an even thicker distribution after formation. Stars of the same age are thicker in the outer disk than in the inner disk (flaring). These lead to positive vertical age gradients and negative radial age gradients. The radial metallicity gradient is neg- ative at the mid-plane, flattens at larger disk height |Z|, and turns positive above |Z|~1.5kpc. The vertical metallicity gradient is negative at all radii, but is steeper at smaller radii. These trends broadly agree with observations in the Milky Way and can be naturally understood from the age gradients. The vertical stellar density profile can be well-described by two components, with scale heights 200-500 pc and 1-1.5 kpc, respectively. The thick component is a mix of stars older than 4 Gyr which formed through a combination of several mechanisms. Our results also demonstrate that it is possible to form a thin disk in cosmological simulations even with strong stellar feedback.

Context. Stars form in dense, dusty clumps of molecular clouds, but little is known about their origin, their evolution and their detailed physical properties. In particular, the relationship between the mass distribution of these clumps (also known as the "clump mass function", or CMF) and the stellar initial mass function (IMF), is still poorly understood. Aims. In order to better understand how the CMF evolve toward the IMF, and to discern the "true" shape of the CMF, large samples of bona-fide pre- and proto-stellar clumps are required. Two such datasets obtained from the Herschel infrared GALactic Plane Survey (Hi-GAL) have been described in paper I. Robust statistical methods are needed in order to infer the parameters describing the models used to fit the CMF, and to compare the competing models themselves. Methods. In this paper we apply Bayesian inference to the analysis of the CMF of the two regions discussed in Paper I. First, we determine the Bayesian posterior probability distribution for each of the fitted parameters. Then, we carry out a quantitative comparison of the models used to fit the CMF. Results. We have compared the results from several methods implementing Bayesian inference, and we have also analyzed the impact of the choice of priors and the influence of various constraints on the statistical conclusions for the preferred values of the parameters. We find that both parameter estimation and model comparison depend on the choice of parameter priors. Conclusions. Our results confirm our earlier conclusion that the CMFs of the two Hi-GAL regions studied here have very similar shapes but different mass scales. Furthermore, the lognormal model appears to better describe the CMF measured in the two Hi-GAL regions studied here. However, this preliminary conclusion is dependent on the choice of parameters priors.

Stars form in dense, dusty clumps of molecular clouds, but little is known about their origin and evolution. In particular, the relationship between the mass distribution of these clumps (also known as the "clump mass function", or CMF) and the stellar initial mass function (IMF), is still poorly understood. In order to discern the "true" shape of the CMF and to better understand how the CMF may evolve toward the IMF, large samples of bona-fide pre- and proto-stellar clumps are required. The sensitive observations of the Herschel Space Observatory (HSO) are now allowing us to look at large clump populations in various clouds with different physical conditions. We analyse two fields in the Galactic plane mapped by HSO during its science demonstration phase, as part of the more complete and unbiased Herschel infrared GALactic Plane Survey (Hi-GAL). These fields undergo a source-extraction and flux-estimation pipeline, which allows us to obtain a sample with thousands of clumps. Starless and proto-stellar clumps are separated using both color and positional criteria to find those coincident with MIPS 24 micron sources. We describe the probability density functions of the power-law and lognormal models that are used to fit the CMFs, and we then find their best-fit parameters. For the lognormal model we apply several statistical techniques to the data and compare their results. The CMFs of the two SDP fields show very similar shapes, but very different mass scales. This similarity is confirmed by the values of the best-fit parameters of either the power-law or lognormal model. The power-law model leads to almost identical CMF slopes, whereas the lognormal model shows that the CMFs have similar widths. The similar CMF shape but different mass scale represents an evidence that the overall process of star formation in the two regions is very different.

The Vela-D region, according to the nomenclature given by Murphy & May (1991), of the star forming complex known as the Vela Molecular Ridge (VMR), has been recently analyzed in details by Olmi et al. (2009), who studied the physical properties of 141 pre- and proto-stellar cold dust cores, detected by the ``Balloon-borne Large-Aperture Submillimeter Telescope'' (BLAST) during a much larger (55 sq. degree) Galactic Plane survey encompassing the whole VMR. This survey's primary goal was to identify the coldest, dense dust cores possibly associated with the earliest phases of star formation. In this work, the dynamical state of the Vela-D cores is analyzed. Comparison to dynamical masses of a sub-sample of the Vela-D cores estimated from the 13CO survey of Elia et al. (2007), is complicated by the fact that the 13CO linewidths are likely to trace the lower density intercore material, in addition to the dense gas associated with the compact cores observed by BLAST. In fact, the total internal pressure of these cores, if estimated using the 13CO linewidths, appears to be higher than the cloud ambient pressure. If this were the case, then self-gravity and surface pressure would be insufficient to bind these cores and an additional source of external confinement (e.g., magnetic field pressure) would be required. However, if one attempts to scale down the 13CO linewidths, according to the observations of high-density tracers in a small sample of sources, then most proto-stellar cores would result effectively gravitationally bound.

The Balloon-borne Large-Aperture Submillimeter Telescope (BLAST) carried out a 250, 350 and 500 micron survey of the galactic plane encompassing the Vela Molecular Ridge, with the primary goal of identifying the coldest dense cores possibly associated with the earliest stages of star formation. Here we present the results from observations of the Vela-D region, covering about 4 square degrees, in which we find 141 BLAST cores. We exploit existing data taken with the Spitzer MIPS, IRAC and SEST-SIMBA instruments to constrain their (single-temperature) spectral energy distributions, assuming a dust emissivity index beta = 2.0. This combination of data allows us to determine the temperature, luminosity and mass of each BLAST core, and also enables us to separate starless from proto-stellar sources. We also analyze the effects that the uncertainties on the derived physical parameters of the individual sources have on the overall physical properties of starless and proto-stellar cores, and we find that there appear to be a smooth transition from the pre- to the proto-stellar phase. In particular, for proto-stellar cores we find a correlation between the MIPS24 flux, associated with the central protostar, and the temperature of the dust envelope. We also find that the core mass function of the Vela-D cores has a slope consistent with other similar (sub)millimeter surveys.