The SPT 0311–58 system at
Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher.
Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?
Some links on this page may take you to non-federal websites. Their policies may differ from this site.
-
Abstract z = 6.900 is an extremely massive structure within the reionization epoch and offers a chance to understand the formation of galaxies at an extreme peak in the primordial density field. We present 70 mas Atacama Large Millimeter/submillimeter Array observations of the dust continuum and [Cii ] 158μ m emission in the central pair of galaxies and reach physical resolutions of ∼100–350 pc, among the most detailed views of any reionization-era system to date. The observations resolve the source into at least a dozen kiloparsec-size clumps. The global kinematics and high turbulent velocity dispersion within the galaxies present a striking contrast to recent claims of dynamically cold thin-disk kinematics in some dusty galaxies just 800 Myr later atz ∼ 4. We speculate that both gravitational interactions and fragmentation from massive parent disks have likely played a role in the overall dynamics and formation of clumps in the system. Each clump individually is comparable in mass to other 6 <z < 8 galaxies identified in rest-UV/optical deep field surveys, but with star formation rates elevated by a factor of ~3-5. Internally, the clumps themselves bear close resemblance to greatly scaled-up versions of virialized cloud-scale structures identified in low-redshift galaxies. Our observationsmore » -
Abstract The Baldwin, Philips, & Terlevich diagram of [O
iii ]/Hβ versus [Nii ]/Hα (hereafter N2-BPT) has long been used as a tool for classifying galaxies based on the dominant source of ionizing radiation. Recent observations have demonstrated that galaxies atz ∼ 2 reside offset from local galaxies in the N2-BPT space. In this paper, we conduct a series of controlled numerical experiments to understand the potential physical processes driving this offset. We model nebular line emission in a large sample of galaxies, taken from thesimba cosmological hydrodynamic galaxy formation simulation, using thecloudy photoionization code to compute the nebular line luminosities from Hii regions. We find that the observed shift toward higher [Oiii ]/Hβ and [Nii ]/Hα values at high redshift arises from sample selection: when we consider only the most massive galaxiesM *∼ 1010–11M ⊙, the offset naturally appears, due to their high metallicities. We predict that deeper observations that probe lower-mass galaxies will reveal galaxies that lie on a locus comparable toz ∼ 0 observations. Even when accounting for samples-selection effects, we find that there is a subtle mismatch between simulations and observations. To resolve this discrepancy, we investigate the impact of varying ionization parameters, Hii region densities, gas-phase abundance patterns, and increasing radiation field hardness on N2-BPT diagrams. We find that either decreasing themore » -
Abstract One of the most common methods for inferring galaxy attenuation curves is via spectral energy distribution (SED) modeling, where the dust attenuation properties are modeled simultaneously with other galaxy physical properties. In this paper, we assess the ability of SED modeling to infer these dust attenuation curves from broadband photometry, and suggest a new flexible model that greatly improves the accuracy of attenuation curve derivations. To do this, we fit mock SEDs generated from the simba cosmological simulation with the prospector SED fitting code. We consider the impact of the commonly assumed uniform screen model and introduce a new nonuniform screen model parameterized by the fraction of unobscured stellar light. This nonuniform screen model allows for a nonzero fraction of stellar light to remain unattenuated, resulting in a more flexible attenuation curve shape by decoupling the shape of the UV attenuation curve from the optical attenuation curve. The ability to constrain the dust attenuation curve is significantly improved with the use of a nonuniform screen model, with the median offset in UV attenuation decreasing from −0.30 dex with a uniform screen model to −0.17 dex with the nonuniform screen model. With this increase in dust attenuation modeling accuracy, wemore »Free, publicly-accessible full text available May 1, 2023
-
Abstract We investigate the fine-structure [C ii ] line at 158 μ m as a molecular gas tracer by analyzing the relationship between molecular gas mass ( M mol ) and [C ii ] line luminosity ( L [C II ] ) in 11,125 z ≃ 6 star-forming, main-sequence galaxies from the simba simulations, with line emission modeled by the Simulator of Galaxy Millimeter/Submillimeter Emission. Though most (∼50%–100%) of the gas mass in our simulations is ionized, the bulk (>50%) of the [C ii ] emission comes from the molecular phase. We find a sublinear (slope 0.78 ± 0.01) log L [ C II ] – log M mol relation, in contrast with the linear relation derived from observational samples of more massive, metal-rich galaxies at z ≲ 6. We derive a median [C ii ]-to- M mol conversion factor of α [C II ] ≃ 18 M ⊙ / L ⊙ . This is lower than the average value of ≃30 M ⊙ / L ⊙ derived from observations, which we attribute to lower gas-phase metallicities in our simulations. Thus, a lower, luminosity-dependent conversion factor must be applied when inferring molecular gas masses from [C ii ] observations ofmore »Free, publicly-accessible full text available April 1, 2023
-
Abstract The H i gas content is a key ingredient in galaxy evolution, the study of which has been limited to moderate cosmological distances for individual galaxies due to the weakness of the hyperfine H i 21 cm transition. Here we present a new approach that allows us to infer the H i gas mass M HI of individual galaxies up to z ≈ 6, based on a direct measurement of the [C ii ]-to-H i conversion factor in star-forming galaxies at z ≳ 2 using γ -ray burst afterglows. By compiling recent [C ii ]-158 μ m emission line measurements we quantify the evolution of the H i content in galaxies through cosmic time. We find that M HI starts to exceed the stellar mass M ⋆ at z ≳ 1, and increases as a function of redshift. The H i fraction of the total baryonic mass increases from around 20% at z = 0 to about 60% at z ∼ 6. We further uncover a universal relation between the H i gas fraction M HI / M ⋆ and the gas-phase metallicity, which seems to hold from z ≈ 6 to z = 0. The majority of galaxiesmore »
-
Abstract Observations of cold molecular gas reservoirs are critical for understanding the shutdown of star formation in massive galaxies. While dust continuum is an efficient and affordable tracer, this method relies upon the assumption of a “normal” molecular-gas to dust mass ratio, δ GDR , typically of order 100. Recent null detections of quiescent galaxies in deep dust continuum observations support a picture where the cold gas and dust have been rapidly depleted or expelled. In this work, we present another viable explanation: a significant fraction of galaxies with low star formation per unit stellar mass are predicted to have extreme δ GDR ratios. We show that simulated massive quiescent galaxies at 0 < z < 3 in the simba cosmological simulations have δ GDR values that extend >4 orders of magnitude. The dust in most simulated quiescent galaxies is destroyed significantly more rapidly than the molecular gas depletes, and cannot be replenished. The transition from star-forming to quiescent halts dust formation via star formation processes, with dust subsequently destroyed by supernova shocks and thermal sputtering of dust grains embedded in hot plasma. After this point, the dust growth rate in the models is not sufficient to overcome the lossmore »
-
ABSTRACT We develop a cosmological model for the evolution of dust grains in galaxies with a distribution of sizes in order to understand the origin of the Milky Way dust extinction curve. Our model considers the formation of active dust in evolved stars, growth by accretion and coagulation, and destruction processes via shattering, sputtering, and astration in the ISM of galaxies over cosmic time. Our main results follow. Galaxies in our cosmological model with masses comparable to the Milky Way’s at z ∼ 0 exhibit a diverse range of extinction laws, though with slopes and bump strengths comparable to the range observed in the Galaxy. The progenitors of the Milky Way have steeper slopes, and only flatten to slopes comparable to the Galaxy at z ∼ 1. This owes to increased grain growth rates at late times/in high-metallicity environments driving up the ratio of large to small grains, with a secondary dependence on the graphite-to-silicate ratio evolution. The UV bump strengths depend primarily on the graphite-to-silicate ratio, and remain broadly constant in MW-like galaxies between z = 3 and z = 0, though show slight variability. Our models span comparable regions of bump-slope space as sightlines in the Galaxy do, though theremore »
-
null (Ed.)ABSTRACT The relation between infrared excess (IRX) and UV spectral slope (βUV) is an empirical probe of dust properties of galaxies. The shape, scatter, and redshift evolution of this relation are not well understood, however, leading to uncertainties in estimating the dust content and star formation rates (SFRs) of galaxies at high redshift. In this study, we explore the nature and properties of the IRX–βUV relation with a sample of z = 2–6 galaxies ($M_*\approx 10^9\!-\!10^{12}\, \mathrm{M}_\odot$) extracted from high-resolution cosmological simulations (MassiveFIRE) of the Feedback in Realistic Environments (FIRE) project. The galaxies in our sample show an IRX–βUV relation that is in good agreement with the observed relation in nearby galaxies. IRX is tightly coupled to the UV optical depth, and is mainly determined by the dust-to-star geometry instead of total dust mass, while βUV is set both by stellar properties, UV optical depth, and the dust extinction law. Overall, much of the scatter in the IRX–βUV relation of our sample is found to be driven by variations of the intrinsic UV spectral slope. We further assess how the IRX–βUV relation depends on viewing direction, dust-to-metal ratio, birth-cloud structures, and the dust extinction law and we present a simplemore »
-
ABSTRACT Matching the number counts of high-z submillimetre-selected galaxies (SMGs) has been a long-standing problem for galaxy formation models. In this paper, we use 3D dust radiative transfer to model the submm emission from galaxies in the simba cosmological hydrodynamic simulations, and compare predictions to the latest single-dish observational constraints on the abundance of 850 μm-selected sources. We find good agreement with the shape of the integrated 850 μm luminosity function, and the normalization is within 0.25 dex at >3 mJy, unprecedented for a fully cosmological hydrodynamic simulation, along with good agreement in the redshift distribution of bright SMGs. The agreement is driven primarily by simba’s good match to infrared measures of the star formation rate (SFR) function between z = 2 and 4 at high SFRs. Also important is the self-consistent on-the-fly dust model in simba, which predicts, on average, higher dust masses (by up to a factor of 2.5) compared to using a fixed dust-to-metals ratio of 0.3. We construct a light-cone to investigate the effect of far-field blending, and find that 52 per cent of sources are blends of multiple components, which makes a small contribution to the normalization of the bright end of the number counts. We provide new fits to the 850 μm luminosity asmore »