Showing papers on "Star formation published in 2020"

Evolutionary roads leading to low effective spins, high black hole masses, and O1/O2 rates for LIGO/Virgo binary black holes

Abstract: All ten LIGO/Virgo binary black hole (BH-BH) coalescences reported following the O1/O2 runs have near-zero effective spins. There are only three potential explanations for this. If the BH spin magnitudes are large, then: (i) either both BH spin vectors must be nearly in the orbital plane or (ii) the spin angular momenta of the BHs must be oppositely directed and similar in magnitude. Then there is also the possibility that (iii) the BH spin magnitudes are small. We consider the third hypothesis within the framework of the classical isolated binary evolution scenario of the BH-BH merger formation. We test three models of angular momentum transport in massive stars: A mildly efficient transport by meridional currents (as employed in the Geneva code), an efficient transport by the Tayler-Spruit magnetic dynamo (as implemented in the MESA code), and a very-efficient transport (as proposed by Fuller et al.) to calculate natal BH spins. We allow for binary evolution to increase the BH spins through accretion and account for the potential spin-up of stars through tidal interactions. Additionally, we update the calculations of the stellar-origin BH masses, including revisions to the history of star formation and to the chemical evolution across cosmic time. We find that we can simultaneously match the observed BH-BH merger rate density and BH masses and BH-BH effective spins. Models with efficient angular momentum transport are favored. The updated stellar-mass weighted gas-phase metallicity evolution now used in our models appears to be key for obtaining an improved reproduction of the LIGO/Virgo merger rate estimate. Mass losses during the pair-instability pulsation supernova phase are likely to be overestimated if the merger GW170729 hosts a BH more massive than 50âMâS. We also estimate rates of black hole-neutron star (BH-NS) mergers from recent LIGO/Virgo observations. If, in fact. angular momentum transport in massive stars is efficient, then any (electromagnetic or gravitational wave) observation of a rapidly spinning BH would indicate either a very effective tidal spin up of the progenitor star (homogeneous evolution, high-mass X-ray binary formation through case A mass transfer, or a spin-up of a Wolf-Rayet star in a close binary by a close companion), significant mass accretion by the hole, or a BH formation through the merger of two or more BHs (in a dense stellar cluster). (Less)

The Evolution of the Star-Forming Interstellar Medium Across Cosmic Time

Abstract: Over the past decade, increasingly robust estimates of the dense molecular gas content in galaxy populations between redshift z = 0 and the peak of cosmic galaxy/star formation (z ∼ 1–3) have becom

Nuclear star clusters

Abstract: We review the current knowledge about nuclear star clusters (NSCs), the spectacularly dense and massive assemblies of stars found at the centers of most galaxies Recent observational and theoretical works suggest that many NSC properties, including their masses, densities, and stellar populations, vary with the properties of their host galaxies Understanding the formation, growth, and ultimate fate of NSCs, therefore, is crucial for a complete picture of galaxy evolution Throughout the review, we attempt to combine and distill the available evidence into a coherent picture of NSC evolution Combined, this evidence points to a clear transition mass in galaxies of $$\sim 10^9\,M_\odot$$ where the characteristics of nuclear star clusters change We argue that at lower masses, NSCs are formed primarily from globular clusters that inspiral into the center of the galaxy, while at higher masses, star formation within the nucleus forms the bulk of the NSC We also discuss the co-existence of NSCs and central black holes, and how their growth may be linked The extreme densities of NSCs and their interaction with massive black holes lead to a wide range of unique phenomena including tidal disruption and gravitational-wave events Finally, we review the evidence that many NSCs end up in the halos of massive galaxies stripped of the stars that surrounded them, thus providing valuable tracers of the galaxies’ accretion histories

The VLA/ALMA Nascent Disk and Multiplicity (VANDAM) Survey of Orion Protostars. A Statistical Characterization of Class 0 and I Protostellar Disks

Abstract: We have conducted a survey of 328 protostars in the Orion molecular clouds with ALMA at 0.87 mm at a resolution of $\sim$0.1" (40 au), including observations with the VLA at 9 mm toward 148 protostars at a resolution of $\sim$0.08" (32 au). This is the largest multi-wavelength survey of protostars at this resolution by an order of magnitude. We use the dust continuum emission at 0.87 mm and 9 mm to measure the dust disk radii and masses toward the Class 0, Class I, and Flat Spectrum protostars, characterizing the evolution of these disk properties in the protostellar phase. The mean dust disk radii for the Class 0, Class I, and Flat Spectrum protostars are 44.9$^{+5.8}_{-3.4}$, 37.0$^{+4.9}_{-3.0}$, and 28.5$^{+3.7}_{-2.3}$ au, respectively, and the mean protostellar dust disk masses are 25.9$^{+7.7}_{-4.0}$, 14.9$^{+3.8}_{-2.2}$, 11.6$^{+3.5}_{-1.9}$ Earth masses, respectively. The decrease in dust disk masses is expected from disk evolution and accretion, but the decrease in disk radii may point to the initial conditions of star formation not leading to the systematic growth of disk radii or that radial drift is keeping the dust disk sizes small. At least 146 protostellar disks (35% out of 379 detected 0.87 mm continuum sources plus 42 non-detections) have disk radii greater than 50 au in our sample. These properties are not found to vary significantly between different regions within Orion. The protostellar dust disk mass distributions are systematically larger than that of Class II disks by a factor of $>$4, providing evidence that the cores of giant planets may need to at least begin their formation during the protostellar phase.

The AGORA high-resolution galaxy simulations comparison project: Public data release

Abstract: As part of the AGORA High-resolution Galaxy Simulations Comparison Project (Kim et al. 2014, 2016) we have generated a suite of isolated Milky Way-mass galaxy simulations using 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical galaxy formation community. In these simulations we adopted identical galactic disk initial conditions, and common physics models (e.g., radiative cooling and ultraviolet background by a standardized package). Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production were carefully constrained. Here we release the simulation data to be freely used by the community. In this release we include the disk snapshots at 0 and 500Myr of evolution per each code as used in Kim et al. (2016), from simulations with and without star formation and feedback. We encourage any member of the numerical galaxy formation community to make use of these resources for their research - for example, compare their own simulations with the AGORA galaxies, with the common analysis yt scripts used to obtain the plots shown in our papers, also available in this release.

A compendium of distances to molecular clouds in the Star Formation Handbook

Abstract: Accurate distances to local molecular clouds are critical for understanding the star and planet formation process, yet distance measurements are often obtained inhomogeneously on a cloud-by-cloud basis. We have recently developed a method that combines stellar photometric data with Gaia DR2 parallax measurements in a Bayesian framework to infer the distances of nearby dust clouds to a typical accuracy of ∼5%. After refining the technique to target lower latitudes and incorporating deep optical data from DECam in the southern Galactic plane, we have derived a catalog of distances to molecular clouds in Reipurth (2008, Star Formation Handbook, Vols. I and II) which contains a large fraction of the molecular material in the solar neighborhood. Comparison with distances derived from maser parallax measurements towards the same clouds shows our method produces consistent distances with ≲10% scatter for clouds across our entire distance spectrum (150 pc−2.5 kpc). We hope this catalog of homogeneous distances will serve as a baseline for future work.

Kraken reveals itself – the merger history of the Milky Way reconstructed with the E-MOSAICS simulations

Abstract: Globular clusters (GCs) formed when the Milky Way experienced a phase of rapid assembly. We use the wealth of information contained in the Galactic GC population to quantify the properties of the satellite galaxies from which the Milky Way assembled. To achieve this, we train an artificial neural network on the E-MOSAICS cosmological simulations of the co-formation and co-evolution of GCs and their host galaxies. The network uses the ages, metallicities, and orbital properties of GCs that formed in the same progenitor galaxies to predict the stellar masses and accretion redshifts of these progenitors. We apply the network to Galactic GCs associated with five progenitors: {\it Gaia}-Enceladus, the Helmi streams, Sequoia, Sagittarius, and the recently discovered, `low-energy' GCs, which provide an excellent match to the predicted properties of the enigmatic galaxy `Kraken'. The five galaxies cover a narrow stellar mass range [$M_\star=(0.6{-}4.6)\times10^8~{\rm M}_\odot$], but have widely different accretion redshifts ($z_{\rm acc}=0.57{-}2.65$). All accretion events represent minor mergers, but Kraken likely represents the most major merger ever experienced by the Milky Way, with stellar and virial mass ratios of $r_{M_\star}=1$:$31^{+34}_{-16}$ and $r_{M_{\rm h}}=1$:$7^{+4}_{-2}$, respectively. The progenitors match the $z=0$ relation between GC number and halo virial mass, but have elevated specific frequencies, suggesting an evolution with redshift. Even though these progenitors likely were the Milky Way's most massive accretion events, they contributed a total mass of only $\log{(M_{\rm \star,tot}/{\rm M}_\odot)}=9.0\pm0.1$, similar to the stellar halo. This implies that the Milky Way grew its stellar mass mostly by in-situ star formation. We conclude by organising these accretion events into the most detailed reconstruction to date of the Milky Way's merger tree.

Streams, substructures and the early history of the Milky Way

Abstract: The advent of Gaia's 2nd data release in combination with large spectroscopic surveys are revolutionizing our understanding of the Galaxy. Thanks to these and the knowledge accumulated thus far, a more mature picture of the evolution of the early Milky Way is emerging: * Two of the traditional Galactic components, i.e. the stellar halo and the thick disk, appear to be intimately linked: stars with halo-like kinematics originate in similar proportions, from a "heated" (thick) disk and from debris from a system named Gaia-Enceladus. Gaia-Enceladus was the last big merger event experienced by the Milky Way and probably completed around 10 Gyr ago. The puffed-up stars now present in the halo as a consequence of the merger have thus exposed the existence of a disk component at z ~ 1.8. * The Helmi streams, Sequoia, and Thamnos are amongst the newly uncovered or better characterized merger events. Knowledge of their progenitor's properties, star formation and chemical histories is still incomplete. * Debris' from different objects often overlap in phase-space. A task for the next years will be to use spectroscopic surveys for chemical labelling and to disentangle events from one another using dimensions other than only phase-space, metallicity or [alpha/Fe]. * These surveys will also provide line-of-sight velocities missing for faint stars and more accurate distance determinations for distant objects. The resulting samples of stars will cover a much wider volume of the Galaxy allowing, for example, linking kinematic substructures in the inner halo to spatial overdensities in the outer halo. * All the results obtained so far are in-line with expectations of current cosmological models. Yet, tailored hydrodynamical simulations as well as "constrained" cosmological simulations are needed to push our knowledge of the assembly of the Milky Way back to the earliest times. [abridged]

High-redshift star formation in the Atacama large millimetre/submillimetre array era.

Abstract: The Atacama Large Millimetre/submillimetre Array (ALMA) is currently in the process of transforming our view of star-forming galaxies in the distant (z≳1) universe. Before ALMA, most of what we kne...

Astrochemistry During the Formation of Stars

Abstract: Star-forming regions show a rich and varied chemistry, including the presence of complex organic molecules—in both the cold gas distributed on large scales and the hot regions close to young stars ...

Binary black holes in the pair instability mass gap

Abstract: Pair instability (PI) and pulsational PI prevent the formation of black holes (BHs) with mass $\gtrsim{}60$ M$_\odot$ from single star evolution. Here, we investigate the possibility that BHs with mass in the PI gap form via stellar mergers and multiple stellar mergers, facilitated by dynamical encounters in young star clusters. We analyze $10^4$ simulations, run with the direct N-body code nbody6++gpu coupled with the population synthesis code MOBSE. We find that up to $\sim{}6$~% of all simulated BHs have mass in the PI gap, depending on progenitor's metallicity. This formation channel is strongly suppressed in metal-rich (Z = 0.02) star clusters, because of stellar winds. BHs with mass in the PI gap are initially single BHs but can efficiently acquire companions through dynamical exchanges. We find that $\sim{}$21%, 10% and 0.5% of all binary BHs have at least one component in the PI mass gap at metallicity Z = 0.0002, 0.002 and 0.02, respectively. Based on the evolution of the cosmic star formation rate and metallicity, and under the assumption that all stars form in young star clusters, we predict that $\sim{}5$~% of all binary BH mergers detectable by advanced LIGO and Virgo at their design sensitivity have at least one component in the PI mass gap.

Strong dependence of Type Ia supernova standardization on the local specific star formation rate

Abstract: As part of an on-going effort to identify, understand and correct for astrophysics biases in the standardization of Type Ia supernovae (SN Ia) for cosmology, we have statistically classified a large sample of nearby SNe Ia into those that are located in predominantly younger or older environments. This classification is based on the specific star formation rate measured within a projected distance of 1 kpc from each SN location (LsSFR). This is an important refinement compared to using the local star formation rate directly, as it provides a normalization for relative numbers of available SN progenitors and is more robust against extinction by dust. We find that the SNe Ia in predominantly younger environments are ΔY = 0.163 ± 0.029 mag (5.7σ) fainter than those in predominantly older environments after conventional light-curve standardization. This is the strongest standardized SN Ia brightness systematic connected to the host-galaxy environment measured to date. The well-established step in standardized brightnesses between SNe Ia in hosts with lower or higher total stellar masses is smaller, at ΔM = 0.119 ± 0.032 mag (4.5σ), for the same set of SNe Ia. When fit simultaneously, the environment-age offset remains very significant, with ΔY = 0.129 ± 0.032 mag (4.0σ), while the global stellar mass step is reduced to ΔM = 0.064 ± 0.029 mag (2.2σ). Thus, approximately 70% of the variance from the stellar mass step is due to an underlying dependence on environment-based progenitor age. Also, we verify that using the local star formation rate alone is not as powerful as LsSFR at sorting SNe Ia into brighter and fainter subsets. Standardization that only uses the SNe Ia in younger environments reduces the total dispersion from 0.142 ± 0.008 mag to 0.120 ± 0.010 mag. We show that as environment-ages evolve with redshift, a strong bias, especially on the measurement of the derivative of the dark energy equation of state, can develop. Fortunately, data that measure and correct for this effect using our local specific star formation rate indicator, are likely to be available for many next-generation SN Ia cosmology experiments.

The VLA/ALMA Nascent Disk and Multiplicity (VANDAM) Survey of Orion Protostars. II. A Statistical Characterization of Class 0 and Class I Protostellar Disks

Abstract: We have conducted a survey of 328 protostars in the Orion molecular clouds with the Atacama Large Millimeter/submillimeter Array at 0.87 mm at a resolution of ∼0.″1 (40 au), including observations with the Very Large Array at 9 mm toward 148 protostars at a resolution of ∼0.″08 (32 au). This is the largest multiwavelength survey of protostars at this resolution by an order of magnitude. We use the dust continuum emission at 0.87 and 9 mm to measure the dust disk radii and masses toward the Class 0, Class I, and flat-spectrum protostars, characterizing the evolution of these disk properties in the protostellar phase. The mean dust disk radii for the Class 0, Class I, and flat-spectrum protostars are 44.9-3.4+5.8, 37.0-3.0+4.9, and 28.5-2.3+3.7 au, respectively, and the mean protostellar dust disk masses are 25.9-4.0+7.7, 14.9-2.2+3.8, 11.6-1.9+3.5 M⊙, respectively. The decrease in dust disk masses is expected from disk evolution and accretion, but the decrease in disk radii may point to the initial conditions of star formation not leading to the systematic growth of disk radii or that radial drift is keeping the dust disk sizes small. At least 146 protostellar disks (35% of 379 detected 0.87 mm continuum sources plus 42 nondetections) have disk radii greater than 50 au in our sample. These properties are not found to vary significantly between different regions within Orion. The protostellar dust disk mass distributions are systematically larger than those of Class II disks by a factor of >4, providing evidence that the cores of giant planets may need to at least begin their formation during the protostellar phase.

But what about.: Cosmic rays, magnetic fields, conduction, and viscosity in galaxy formation

Abstract: We present and study a large suite of high-resolution cosmological zoom-in simulations, using the FIRE-2 treatment of mechanical and radiative feedback from massive stars, together with explicit treatment of magnetic fields, anisotropic conduction and viscosity (accounting for saturation and limitation by plasma instabilities at high β), and cosmic rays (CRs) injected in supernovae shocks (including anisotropic diffusion, streaming, adiabatic, hadronic and Coulomb losses). We survey systems from ultrafaint dwarf (⁠M∗∼10⁴M⊙⁠, M_(halo)∼10⁹M⊙⁠) through Milky Way/Local Group (MW/LG) masses, systematically vary uncertain CR parameters (e.g. the diffusion coefficient κ and streaming velocity), and study a broad ensemble of galaxy properties [masses, star formation (SF) histories, mass profiles, phase structure, morphologies, etc.]. We confirm previous conclusions that magnetic fields, conduction, and viscosity on resolved (⁠≳1 pc) scales have only small effects on bulk galaxy properties. CRs have relatively weak effects on all galaxy properties studied in dwarfs (⁠M∗≪10¹⁰M⊙, M_(halo)≲10¹¹M⊙⁠), or at high redshifts (z ≳ 1–2), for any physically reasonable parameters. However, at higher masses (⁠M_(halo)≳10¹¹M⊙) and z ≲ 1–2, CRs can suppress SF and stellar masses by factors ∼2–4, given reasonable injection efficiencies and relatively high effective diffusion coefficients κ≳3×10²⁹cm²s⁻¹⁠. At lower κ, CRs take too long to escape dense star-forming gas and lose their energy to collisional hadronic losses, producing negligible effects on galaxies and violating empirical constraints from spallation and γ-ray emission. At much higher κ CRs escape too efficiently to have appreciable effects even in the CGM. But around κ∼3×10²⁹cm²s⁻¹⁠, CRs escape the galaxy and build up a CR-pressure-dominated halo which maintains approximate virial equilibrium and supports relatively dense, cool (T ≪ 10⁶ K) gas that would otherwise rain on to the galaxy. CR ‘heating’ (from collisional and streaming losses) is never dominant.

Molecular outflows in local galaxies: Method comparison and a role of intermittent AGN driving

Abstract: We report new detections and limits from a NOEMA and ALMA CO(1-0) search for molecular outflows in 13 local galaxies with high far-infrared surface brightness, and combine these with local universe CO outflow results from the literature. The CO line ratios and spatial outflow structure of our targets provide some constraints on the conversion steps from observables to physical quantities such as molecular mass outflow rates. Where available, ratios between outflow emission in higher J CO transitions and in CO(1-0) are typically consistent with excitation R-i1 less than or similar to 1. However, for IRAS 13120 5453, R-31 = 2.10 +/- 0.29 indicates optically thin CO in the outflow. Like much of the outflow literature, we use ff CO(1 0) = 0.8, and we present arguments for using C = 1 in deriving molecular mass outflow rates. (M)over dot(out) = CM(out)v(out)/R-out. We compare the two main methods for molecular outflow detection: CO millimeter interferometry and Herschel OH-based spectroscopic outflow searches. For 26 sources studied with both methods, we find an 80% agreement in detecting vout & 150 km s 1 outflows, and non-matches can be plausibly ascribed to outflow geometry and signal-to-noise ratio. For a published sample of 12 bright ultraluminous infrared galaxies with detailed OH-based outflow modeling, CO outflows are detected in all but one. Outflow masses, velocities, and sizes for these 11 sources agree well between the two methods, and modest remaining di fferences may relate to the di fferent but overlapping regions sampled by CO emission and OH absorption. Outflow properties correlate better with active galactic nucleus (AGN) luminosity and with bolometric luminosity than with far-infrared surface brightness. The most massive outflows are found for systems with current AGN activity, but significant outflows in nonAGN systems must relate to star formation or to AGN activity in the recent past. We report scaling relations for the increase of outflow mass, rate, momentum rate, and kinetic power with bolometric luminosity. Short flow times of similar to 10(6) yr and some sources with resolved multiple outflow episodes support a role of intermittent driving, likely by AGNs.

The ALPINE–ALMA [C ii] Survey: Multiwavelength Ancillary Data and Basic Physical Measurements

Abstract: We present the ancillary data and basic physical measurements for the galaxies in the ALMA Large Program to Investigate C⁺ at Early Times (ALPINE) survey—the first large multiwavelength survey that aims at characterizing the gas and dust properties of 118 main-sequence galaxies at redshifts 4.4 3.5σ) and the surrounding far-infrared continuum in conjunction with a wealth of optical and near-infrared data. We outline in detail the spectroscopic data and selection of the galaxies as well as the ground- and space-based imaging products. In addition, we provide several basic measurements including stellar masses, star formation rates (SFR), rest-frame ultra-violet (UV) luminosities, UV continuum slopes (β), and absorption line redshifts, as well as Hα emission derived from Spitzer colors. We find that the ALPINE sample is representative of the 4 < z < 6 galaxy population selected by photometric methods and only slightly biased toward bluer colors (Δβ ~ 0.2). Using [C II] as tracer of the systemic redshift (confirmed for one galaxy at z = 4.5 out of 118 for which we obtained optical [O III]λ3727A emission), we confirm redshifted Lyα emission and blueshifted absorption lines similar to findings at lower redshifts. By stacking the rest-frame UV spectra in the [C II] rest frame, we find that the absorption lines in galaxies with high specific SFR are more blueshifted, which could be indicative of stronger winds and outflows.

Dust masses of young disks: constraining the initial solid reservoir for planet formation

Abstract: Context. Recent years have seen building evidence that planet formation starts early, in the first ~0.5 Myr. Studying the dust masses available in young disks enables us to understand the origin of planetary systems given that mature disks are lacking the solid material necessary to reproduce the observed exoplanetary systems, especially the massive ones.Aims. We aim to determine if disks in the embedded stage of star formation contain enough dust to explain the solid content of the most massive exoplanets.Methods. We use Atacama Large Millimeter/submillimeter Array (ALMA) Band 6 (1.1–1.3 mm) continuum observations of embedded disks in the Perseus star-forming region together with Very Large Array (VLA) Ka -band (9 mm) data to provide a robust estimate of dust disk masses from the flux densities measured in the image plane.Results. We find a strong linear correlation between the ALMA and VLA fluxes, demonstrating that emission at both wavelengths is dominated by dust emission. For a subsample of optically thin sources, we find a median spectral index of 2.5 from which we derive the dust opacity index β = 0.5, suggesting significant dust growth. Comparison with ALMA surveys of Orion shows that the Class I dust disk mass distribution between the two regions is similar, but that the Class 0 disks are more massive in Perseus than those in Orion. Using the DIANA opacity model including large grains, with a dust opacity value of κ 9 mm = 0.28 cm2 g−1 , the median dust masses of the embedded disks in Perseus are 158 M ⊕ for Class 0 and 52 M ⊕ for Class I from the VLA fluxes. The lower limits on the median masses from ALMA fluxes are 47 M ⊕ and 12 M ⊕ for Class 0 and Class I, respectively, obtained using the maximum dust opacity value κ 1.3 mm = 2.3 cm2 g−1 . The dust masses of young Class 0 and I disks are larger by at least a factor of ten and three, respectively, compared with dust masses inferred for Class II disks in Lupus and other regions.Conclusions. The dust masses of Class 0 and I disks in Perseus derived from the VLA data are high enough to produce the observed exoplanet systems with efficiencies acceptable by planet formation models: the solid content in observed giant exoplanets can be explained if planet formation starts in Class 0 phase with an efficiency of ~15%. A higher efficiency of ~30% is necessary if the planet formation is set to start in Class I disks.

Cosmological inference using gravitational wave standard sirens : a mock data analysis

Abstract: The observation of binary neutron star merger GW170817, along with its optical counterpart, provided the first constraint on the Hubble constant H-0 using gravitational wave standard sirens. When no counterpart is identified, a galaxy catalog can be used to provide the necessary redshift information. However, the true host might not be contained in a catalog which is not complete out to the limit of gravitational-wave detectability. These electromagnetic and gravitational-wave selection effects must be accounted for. We describe and implement a method to estimate H-0 using both the counterpart and the galaxy catalog standard siren methods. We perform a series of mock data analyses using binary neutron star mergers to confirm our ability to recover an unbiased estimate of H-0. Our simulations used a simplified universe with no redshift uncertainties or galaxy clustering, but with different magnitude-limited catalogs and assumed host galaxy properties, to test our treatment of both selection effects. We explore how the incompleteness of catalogs affects the final measurement of H-0, as well as the effect of weighting each galaxy's likelihood of being a host by its luminosity. In our most realistic simulation, where the simulated catalog is about three times denser than the density of galaxies in the local universe, we find that a 4.4% measurement precision can be reached using galaxy catalogs with 50% completeness and similar to 250 binary neutron star detections with sensitivity similar to that of Advanced LIGO's second observing run.

The ALPINE-ALMA [C II] survey: Star-formation-driven outflows and circumgalactic enrichment in the early Universe

Abstract: We study the efficiency of galactic feedback in the early Universe by stacking the [C II] 158 μ m emission in a large sample of normal star-forming galaxies at 4 |≲500 km s−1 . The significance of these features increases when stacking the subset of galaxies with star formation rates (SFRs) higher than the median (SFRmed = 25 M ⊙ yr−1 ), thus confirming their star-formation-driven nature. The estimated mass outflow rates are comparable to the SFRs, yielding mass-loading factors of the order of unity (similarly to local star-forming galaxies), suggesting that star-formation-driven feedback may play a lesser role in quenching galaxies at z > 4. From the stacking analysis of the datacubes, we find that the combined [C II] core emission (|v | ) of the higher-SFR galaxies is extended on physical sizes of ∼30 kpc (diameter scale), well beyond the analogous [C II] core emission of lower-SFR galaxies and the stacked far-infrared continuum. The detection of such extended metal-enriched gas, likely tracing circumgalactic gas enriched by past outflows, corroborates previous similar studies, confirming that baryon cycle and gas exchanges with the circumgalactic medium are at work in normal star-forming galaxies already at early epochs.

EDGE: the mass–metallicity relation as a critical test of galaxy formation physics

Abstract: We introduce the 'Engineering Dwarfs at Galaxy Formation's Edge' (EDGE) project to study the cosmological formation and evolution of the smallest galaxies in the Universe. In this first paper, we explore the effects of resolution and sub-grid physics on a single low-mass halo (M_halo=109{ M}_☉), simulated to redshift z = 0 at a mass and spatial resolution of 20{ M}_☉ and ∼3 pc. We consider different star formation prescriptions, supernova feedback strengths, and on-the-fly radiative transfer (RT). We show that RT changes the mode of galactic self-regulation at this halo mass, suppressing star formation by causing the interstellar and circumgalactic gas to remain predominantly warm (∼104 K) even before cosmic reionization. By contrast, without RT, star formation regulation occurs only through starbursts and their associated vigorous galactic outflows. In spite of this difference, the entire simulation suite (with the exception of models without any feedback) matches observed dwarf galaxy sizes, velocity dispersions, V-band magnitudes, and dynamical mass-to-light-ratios. This is because such structural scaling relations are predominantly set by the host dark matter halo, with the remaining model-to-model variation being smaller than the observational scatter. We find that only the stellar mass-metallicity relation differentiates the galaxy formation models. Explosive feedback ejects more metals from the dwarf, leading to a lower metallicity at a fixed stellar mass. We conclude that the stellar mass-metallicity relation of the very smallest galaxies provides a unique constraint on galaxy formation physics.

The cosmic merger rate density evolution of compact binaries formed in young star clusters and in isolated binaries

Abstract: Next generation ground-based gravitational-wave detectors will observe binary black hole (BBH) mergers up to redshift $z\gtrsim{}10$, probing the evolution of compact binary (CB) mergers across cosmic time. Here, we present a new data-driven model to estimate the cosmic merger rate density (MRD) evolution of CBs, by coupling catalogs of CB mergers with observational constraints on the cosmic star formation rate density and on the metallicity evolution of the Universe. We adopt catalogs of CB mergers derived from recent $N-$body and population-synthesis simulations, to describe the MRD of CBs formed in young star clusters (hereafter, dynamical CBs) and in the field (hereafter, isolated CBs). The local MRD of dynamical BBHs is $\mathcal{R}_{\rm BBH}=64^{+34}_{-20}$ Gpc$^{-3}$ yr$^{-1}$, consistent with the 90% credible interval from the first and second observing run (O1 and O2) of the LIGO-Virgo collaboration, and with the local MRD of isolated BBHs ($\mathcal{R}_{\rm BBH}=50^{+71}_{-37}$ Gpc$^{-3}$ yr$^{-1}$). The local MRD of dynamical and isolated black hole - neutron star binaries is $\mathcal{R}_{\rm BHNS}=41^{+33}_{-23}$ and $49^{+48}_{-34}$~Gpc$^{-3}$ yr$^{-1}$, respectively. Both values are consistent with the upper limit inferred from O1 and O2. Finally, the local MRD of dynamical binary neutron stars (BNSs, $\mathcal{R}_{\rm BNS}=151^{+59}_{-38}$ Gpc$^{-3}$ yr$^{-1}$) is a factor of two lower than the local MRD of isolated BNSs ($\mathcal{R}_{\rm BNS}=283^{+97}_{-75}$ Gpc$^{-3}$ yr$^{-1}$). The MRD for all CB classes grows with redshift, reaching its maximum at $z \in [1.5,2.5]$, and then decreases. This trend springs from the interplay between cosmic star formation rate, metallicity evolution and delay time of binary compact objects.

Ejective and preventative: the IllustrisTNG black hole feedback and its effects on the thermodynamics of the gas within and around galaxies

Abstract: Supermassive black holes (SMBHs) that reside at the centres of galaxies can inject vast amounts of energy into the surrounding gas and are thought to be a viable mechanism to quench star formation in massive galaxies. Here, we study the $10^{9-12.5}\, \mathrm{M_\odot }$ stellar mass central galaxy population of the IllustrisTNG simulation, specifically the TNG100 and TNG300 volumes at z = 0, and show how the three components – SMBH, galaxy, and circumgalactic medium (CGM) – are interconnected in their evolution. We find that gas entropy is a sensitive diagnostic of feedback injection. In particular, we demonstrate how the onset of the low-accretion black hole (BH) feedback mode, realized in the IllustrisTNG model as a kinetic, BH-driven wind, leads not only to star formation quenching at stellar masses $\gtrsim 10^{10.5}\, \mathrm{M_\odot }$ but also to a change in thermodynamic properties of the (non-star-forming) gas, both within the galaxy and beyond. The IllustrisTNG kinetic feedback from SMBHs increases the average gas entropy, within the galaxy and in the CGM, lengthening typical gas cooling times from $10\!-\!100\, \mathrm{Myr}$ to $1\!-\!10\, \mathrm{Gyr}$, effectively ceasing ongoing star formation and inhibiting radiative cooling and future gas accretion. In practice, the same active galactic nucleus (AGN) feedback channel is simultaneously ‘ejective’ and ‘preventative’ and leaves an imprint on the temperature, density, entropy, and cooling times also in the outer reaches of the gas halo, up to distances of several hundred kiloparsecs. In the IllustrisTNG model, a long-lasting quenching state can occur for a heterogeneous CGM, whereby the hot and dilute CGM gas of quiescent galaxies contains regions of low-entropy gas with short cooling times.

Blast Waves from Magnetar Flares and Fast Radio Bursts

Abstract: Magnetars younger than one century are expected to be hyper active. Besides winds powered by rotation they generate frequent magnetic flares, which launch powerful blast waves into the wind. These internal shocks act as masers producing fast (millisecond) radio bursts (FRBs) with the following properties. (1) GHz radio emission occurs at radii $r\sim 10^{14}$ cm and lasts $\lesssim 1$ ms in observer's time. (2) Induced scattering in the surrounding wind does not suppress the radio burst. (3) The emission has linear polarization set by the magnetar rotation axis. (4) The emission drifts to lower frequencies during the burst, and its duration broadens at lower frequencies. (5) Blast waves in inhomogeneous winds may emit variable bursts; periodicity might appear on sub-ms timescales if the magnetar rotates with $\sim 1$ s period. However, the observed FRB structure is likely changed by lensing effects during propagation through the host galaxy. (6) The FRBs from magnetars are expected to repeat, with rare strong bursts (up to $\sim 10^{43}$ erg) or more frequent weak bursts. (7) When a repeating flare strikes the wind bubble in the tail of a previous flare, the FRB turns into a bright optical flash. Its luminosity may approach that of a supernova Ia and last seconds. The rate of these optical flashes in the universe is much lower than the FRB rate, however it may exceed the supernova rate. Locations of hyper-active magnetars in their host galaxies depend on how they form: magnetars created in supernovae explosions will trace star formation regions, and magnetars formed in mergers of compact objects will be offset. The merger magnetars are expected to be most energetic and particularly hyper-active.

The 1.28 GHz MeerKAT DEEP2 Image

Abstract: We present the confusion-limited 1.28 GHz MeerKAT DEEP2 image covering one $\approx 68'$ FWHM primary beam area with $7.6''$ FWHM resolution and $0.55 \pm 0.01$ $\mu$Jy/beam rms noise. Its J2000 center position $\alpha=04^h 13^m 26.4^s$, $\delta=-80^\circ 00' 00''$ was selected to minimize artifacts caused by bright sources. We introduce the new 64-element MeerKAT array and describe commissioning observations to measure the primary beam attenuation pattern, estimate telescope pointing errors, and pinpoint $(u,v)$ coordinate errors caused by offsets in frequency or time. We constructed a 1.4 GHz differential source count by combining a power-law count fit to the DEEP2 confusion $P(D)$ distribution from $0.25$ to $10$ $\mu$Jy with counts of individual DEEP2 sources between $10$ $\mu$Jy and $2.5$ mJy. Most sources fainter than $S \sim 100$ $\mu$Jy are distant star-forming galaxies obeying the FIR/radio correlation, and sources stronger than $0.25$ $\mu$Jy account for $\sim93\%$ of the radio background produced by star-forming galaxies. For the first time, the DEEP2 source count has reached the depth needed to reveal the majority of the star formation history of the universe. A pure luminosity evolution of the 1.4 GHz local luminosity function consistent with the Madau & Dickinson (2014) model for the evolution of star-forming galaxies based on UV and infrared data underpredicts our 1.4 GHz source count in the range $-5 \lesssim \log[S(\mathrm{Jy})] \lesssim -4$.

Detection of large-scale X-ray bubbles in the Milky Way halo

Abstract: The halo of the Milky Way provides a laboratory to study the properties of the shocked hot gas that is predicted by models of galaxy formation. There is observational evidence of energy injection into the halo from past activity in the nucleus of the Milky Way1–4; however, the origin of this energy (star formation or supermassive-black-hole activity) is uncertain, and the causal connection between nuclear structures and large-scale features has not been established unequivocally. Here we report soft-X-ray-emitting bubbles that extend approximately 14 kiloparsecs above and below the Galactic centre and include a structure in the southern sky analogous to the North Polar Spur. The sharp boundaries of these bubbles trace collisionless and non-radiative shocks, and corroborate the idea that the bubbles are not a remnant of a local supernova5 but part of a vast Galaxy-scale structure closely related to features seen in γ-rays6. Large energy injections from the Galactic centre7 are the most likely cause of both the γ-ray and X-ray bubbles. The latter have an estimated energy of around 1056 erg, which is sufficient to perturb the structure, energy content and chemical enrichment of the circumgalactic medium of the Milky Way. Observations from the eROSITA telescope reveal soft-X-ray-emitting bubbles extending above and below the Galactic plane, which arose from energy injected into the Galactic halo from past activity in the Galactic centre.

Both starvation and outflows drive galaxy quenching

Abstract: Star-forming galaxies can in principle be transformed into passive systems by a multitude of processes that quench star formation, such as the halting of gas accretion (starvation) or the rapid removal of gas in AGN-driven outflows. However, it remains unclear which processes are the most significant, primary drivers of the SF-passive bimodality. We address this key issue in galaxy evolution by studying the chemical properties of 80,000 local galaxies in SDSS DR7. In order to distinguish between different quenching mechanisms, we analyse the stellar metallicities of star-forming, green valley and passive galaxies. We find that the significant difference in stellar metallicity between passive galaxies and their star-forming progenitors implies that for galaxies at all masses, quenching must have involved an extended phase of starvation. However, some form of gas ejection also has to be introduced into our models to best match the observed properties of local passive galaxies, indicating that, while starvation is likely to be the prerequisite for quenching, it is the combination of starvation and outflows that is responsible for quenching the majority of galaxies. Closed-box models indicate that the duration of the quenching phase is 2-3 Gyr, with an $e$-folding time of 2-4 Gyr, after which further star formation is prevented by an ejective/heating mode. Alternatively, leaky-box models find a longer duration for the quenching phase of 5-7 Gyr and an $e$-folding time of $\sim$1 Gyr, with outflows becoming increasingly important with decreasing stellar mass. Finally, our analysis of local green valley galaxies indicates that quenching is slower in the local Universe than at high-redshift.

The ALMA spectroscopic survey large program: the infrared excess of z = 1.5-10 UV-selected galaxies and the implied high-redshift star formation history

Abstract: We make use of sensitive (9.3 μJy beam−1 rms) 1.2 mm continuum observations from the Atacama Large Millimeter/submillimeter Array (ALMA) Spectroscopic Survey in the Hubble Ultra-Deep Field (ASPECS) large program to probe dust-enshrouded star formation from 1362 Lyman-break galaxies spanning the redshift range z = 1.5–10 (to ~7–28 M ⊙ yr−1 at 4σ over the entire range). We find that the fraction of ALMA-detected galaxies in our z = 1.5–10 samples increases steeply with stellar mass, with the detection fraction rising from 0% at 109.0 M ⊙ to ${85}_{-18}^{+9}$% at >1010 M ⊙. Moreover, on stacking all 1253 low-mass ( ${10}^{9.5}\,{M}_{\odot }$ and an SMC-like relation at lower masses. Using stellar mass and β measurements for z ~ 2 galaxies over the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey, we derive a new empirical relation between β and stellar mass and then use this correlation to show that our IRX–β and IRX–stellar mass relations are consistent with each other. We then use these constraints to express the IRX as a bivariate function of β and stellar mass. Finally, we present updated estimates of star formation rate density determinations at z > 3, leveraging present improvements in the measured IRX and recent probes of ultraluminous far-IR galaxies at z > 2.

The Molecular Cloud Lifecycle.

Abstract: Giant molecular clouds (GMCs) and their stellar offspring are the building blocks of galaxies. The physical characteristics of GMCs and their evolution are tightly connected to galaxy evolution. The macroscopic properties of the interstellar medium propagate into the properties of GMCs condensing out of it, with correlations between e.g. the galactic and GMC scale gas pressures, surface densities and volume densities. That way, the galactic environment sets the initial conditions for star formation within GMCs. After the onset of massive star formation, stellar feedback from e.g. photoionisation, stellar winds, and supernovae eventually contributes to dispersing the parent cloud, depositing energy, momentum and metals into the surrounding medium, thereby changing the properties of galaxies. This cycling of matter between gas and stars, governed by star formation and feedback, is therefore a major driver of galaxy evolution. Much of the recent debate has focused on the durations of the various evolutionary phases that constitute this cycle in galaxies, and what these can teach us about the physical mechanisms driving the cycle. We review results from observational, theoretical, and numerical work to build a dynamical picture of the evolutionary lifecycle of GMC evolution, star formation, and feedback in galaxies.

Quiescent galaxies 1.5 billion years after the Big Bang and their progenitors

Abstract: We report two secure (z = 3.775; 4.012) and one tentative (z ≈ 3.767) spectroscopic confirmations of massive and quiescent galaxies through K-band observations with Keck/MOSFIRE and Very Large Telescope/X-Shooter. The stellar continuum emission, absence of strong nebular emission lines, and lack of significant far-infrared detections confirm the passive nature of these objects, disfavoring the alternative solution of low-redshift dusty star-forming interlopers. We derive stellar masses of log(M⋆/M⊙) ~ 11 and ongoing star formation rates placing these galaxies ≳ 1–2 dex below the main sequence at their redshifts. The adopted parameterization of the star formation history suggests that these sources experienced a strong (〈SFR〉 ~1200 - 3500 M⊙ yr⁻¹) and short (~50 Myr) burst of star formation, peaking ~150–500 Myr before the time of observation, all properties reminiscent of the characteristics of submillimeter galaxies (SMGs) at z > 4. We investigate this connection by comparing the comoving number densities and the properties of these two populations. We find a fair agreement only with the deepest submillimeter surveys detecting not only the most extreme starbursts but also more normal galaxies. We support these findings by further exploring the Illustris TNG cosmological simulation, retrieving populations of both fully quenched massive galaxies at z ~ 3–4 and SMGs at z ~ 4−5, with number densities and properties in agreement with the observations at z ~ 3 but in increasing tension at higher redshift. Nevertheless, as suggested by the observations, not all of the progenitors of quiescent galaxies at these redshifts shine as bright SMGs in their past, and, similarly, not all bright SMGs quench by z ~ 3, both fractions depending on the threshold assumed to define the SMGs themselves.