No Arabic abstract
We use deep Chandra imaging to measure the distribution of X-ray luminosities (L_X) for samples of star-forming galaxies as a function of stellar mass and redshift, using a Bayesian method to push below the nominal X-ray detection limits. Our luminosity distributions all show narrow peaks at L_X < 10^{42} erg/s that we associate with star formation, as opposed to AGN that are traced by a broad tail to higher L_X. Tracking the luminosity of these peaks as a function of stellar mass reveals an X-ray main sequence with a constant slope ~0.63 +/- 0.03 over 8.5 < log M*/Msun < 11.5 and 0.1 < z < 4, with a normalization that increases with redshift as (1+z)^{3.79+/-0.12}. We also compare the peak X-ray luminosities with UV-to-IR tracers of star formation rates (SFRs) to calibrate the scaling between L_X and SFR. We find that L_X propto SFR^{0.83} x (1+z)^{1.3}, where the redshift evolution and non-linearity likely reflect changes in high-mass X-ray binary populations of star-forming galaxies. Using galaxies with a broader range of SFR, we also constrain a stellar-mass-dependent contribution to L_X, likely related to low-mass X-ray binaries. Using this calibration, we convert our X-ray main sequence to SFRs and measure a star-forming main sequence with a constant slope ~0.76+/-0.06 and a normalization that evolves with redshift as (1+z)^{2.95+/-0.33}. Based on the X-ray emission, there is no evidence for a break in the main sequence at high stellar masses, although we cannot rule out a turnover given the uncertainties in the scaling of L_X to SFR.
In the Milky Way, cosmic rays (CRs) are dynamically important in the interstellar medium, contribute to hydrostatic balance, and may help regulate star formation. However, we know far less about the importance of CRs in galaxies whose gas content or star formation rate differ significantly from those of the Milky Way. Here we construct self-consistent models for hadronic CR transport, losses, and contribution to pressure balance as a function of galaxy properties, covering a broad range of parameters from dwarfs to extreme starbursts. While the CR energy density increases from $sim 1$ eV cm$^{-3}$ to $sim 1$ keV cm$^{-3}$ over the range from sub-Milky Way dwarfs to bright starbursts, strong hadronic losses render CRs increasingly unimportant dynamically as the star formation rate surface density increases. In Milky Way-like systems, CR pressure is typically comparable to turbulent gas and magnetic pressure at the galactic midplane, but the ratio of CR pressure to gas pressure drops to $sim 10^{-3}$ in dense starbursts. Galaxies also become increasingly CR calorimetric and gamma-ray bright in this limit. The degree of calorimetry at fixed galaxy properties is sensitive to the assumed model for CR transport, and in particular to the time CRs spend interacting with neutral ISM, where they undergo strong streaming losses. We also find that in some regimes of parameter space hydrostatic equilibrium discs cannot exist, and in Paper II of this series we use this result to derive a critical surface in the plane of star formation surface density and gas surface density beyond which CRs may drive large-scale galactic winds.
We map the co-eval growth of galaxies and their central supermassive black holes in detail by measuring the incidence of Active Galactic Nuclei (AGN) in galaxies as a function of star formation rate (SFR) and redshift (to z~4). We combine large galaxy samples with deep Chandra X-ray imaging to measure the probability distribution of specific black hole accretion rates (LX relative to stellar mass) and derive robust AGN fractions and average specific accretion rates. First, we consider galaxies along the main sequence of star formation. We find a linear correlation between the average SFR and both the AGN fraction and average specific accretion rate across a wide range in stellar mass ($M_* sim 10^{8.5-11.5}M_odot$) and to at least z~2.5, indicating that AGN in main-sequence galaxies are driven by the stochastic accretion of cold gas. We also consider quiescent galaxies and find significantly higher AGN fractions than predicted, given their low SFRs, indicating that AGN in quiescent galaxies are fuelled by additional mechanisms (e.g. stellar winds). Next, we bin galaxies according to their SFRs relative to the main sequence. We find that the AGN fraction is significantly elevated for galaxies that are still star-forming but with SFRs below the main sequence, indicating further triggering mechanisms enhance AGN activity within these sub-main-sequence galaxies. We also find that the incidence of high-accretion-rate AGN is enhanced in starburst galaxies and evolves more mildly with redshift than within the rest of the galaxy population, suggesting mergers play a role in driving AGN activity in such high-SFR galaxies.
Cosmic rays (CRs) are a plausible mechanism for launching winds of cool material from the discs of star-forming galaxies. However, there is no consensus on what types of galaxies likely host CR-driven winds, or what role these winds might play in regulating galaxies star formation rates. Using a detailed treatment of the transport and losses of hadronic CRs developed in the previous paper in this series, here we develop a semi-analytic model that allows us to assess the viability of using CRs to launch cool winds from galactic discs. In particular, we determine the critical CR fluxes -- and corresponding star formation rate surface densities -- above which hydrostatic equilibrium within a given galaxy is precluded because CRs drive the gas off in a wind or otherwise render it unstable. We show that, for star-forming galaxies with lower gas surface densities typical of the Galaxy and local dwarfs, the locus of this CR stability curve patrols the high side of the observed distribution of galaxies in the Kennicutt-Schmidt parameter space of star formation rate versus gas surface density. However, hadronic losses render CRs unable to drive winds in galaxies with higher surface densities. Our results show that quiescent, low surface density galaxies like the Milky Way are poised on the cusp of instability, such that small changes to ISM parameters can lead to the launching of CR-driven outflows, and we suggest that, as a result, CR feedback sets an ultimate limit to the star formation efficiency of most modern galaxies.
We present the ~800 star formation rate maps for the SAMI Galaxy Survey based on H{alpha} emission maps, corrected for dust attenuation via the Balmer decrement, that are included in the SAMI Public Data Release 1. We mask out spaxels contaminated by non-stellar emission using the [O III]/H{beta}, [N II]/H{alpha}, [S II]/H{alpha}, and [O I]/H{alpha} line ratios. Using these maps, we examine the global and resolved star-forming main sequences of SAMI galaxies as a function of morphology, environmental density, and stellar mass. Galaxies further below the star-forming main sequence are more likely to have flatter star formation profiles. Early-type galaxies split into two populations with similar stellar masses and central stellar mass surface densities. The main sequence population has centrally-concentrated star formation similar to late-type galaxies, while galaxies >3{sigma} below the main sequence show significantly reduced star formation most strikingly in the nuclear regions. The split populations support a two-step quenching mechanism, wherein halo mass first cuts off the gas supply and remaining gas continues to form stars until the local stellar mass surface density can stabilize the reduced remaining fuel against further star formation. Across all morphologies, galaxies in denser environments show a decreased specific star formation rate from the outside in, supporting an environmental cause for quenching, such as ram-pressure stripping or galaxy interactions.
We present a framework for modelling the star-formation histories of galaxies as a stochastic process. We define this stochastic process through a power spectrum density with a functional form of a broken power-law. Star-formation histories are correlated on short timescales, the strength of this correlation described by a power-law slope, $alpha$, and they decorrelate to resemble white noise over a timescale that is proportional to the timescale of the break in the power spectrum density, $tau_{rm break}$. We use this framework to explore the properties of the stochastic process that, we assume, gives rise to the log-normal scatter about the relationship between star-formation rate and stellar mass, the so-called galaxy star-forming main sequence. Specifically, we show how the measurements of the normalisation and width ($sigma_{rm MS}$) of the main sequence, measured in several passbands that probe different timescales, give a constraint on the parameters of the underlying power spectrum density. We first derive these results analytically for a simplified case where we model observations by averaging over the recent star-formation history. We then run numerical simulations to find results for more realistic observational cases. As a proof of concept, we use observational estimates of the main sequence scatter at $zsim0$ and $M_{star}approx10^{10}~M_{odot}$ measured in H$alpha$, UV+IR and the u-band, and show that combination of these point to $tau_{rm break}=178^{+104}_{-66}$ Myr, when assuming $alpha=2$. This implies that star-formation histories of galaxies lose memory of their previous activity on a timescale of $sim200$ Myr, highlighting the importance of baryonic effects that act over the dynamical timescales of galaxies.