Eliot Quataert

Galaxies on FIRE (Feedback In Realistic Environments): Stellar feedback explains cosmologically inefficient star formation

We present a series of high-resolution cosmological simulations of galaxy formation to z = 0, spanning halo masses ∼10^8–10^(13) M⊙, and stellar masses ∼10^4–10^(11) M⊙. Our simulations include fully explicit treatment of the multiphase interstellar medium and stellar feedback. The stellar feedback inputs (energy, momentum, mass, and metal fluxes) are taken directly from stellar population models. These sources of feedback, with zero adjusted parameters, reproduce the observed relation between stellar and halo mass up to M_(halo) ∼ 10^(12) M⊙. We predict weak redshift evolution in the M*–M_(halo) relation, consistent with current constraints to z > 6. We find that the M*–M_(halo) relation is insensitive to numerical details, but is sensitive to feedback physics. Simulations with only supernova feedback fail to reproduce observed stellar masses, particularly in dwarf and high-redshift galaxies: radiative feedback (photoheating and radiation pressure) is necessary to destroy giant molecular clouds and enable efficient coupling of later supernovae to the gas. Star formation rates (SFRs) agree well with the observed Kennicutt relation at all redshifts. The galaxy-averaged Kennicutt relation is very different from the numerically imposed law for converting gas into stars, and is determined by self-regulation via stellar feedback. Feedback reduces SFRs and produces reservoirs of gas that lead to rising late-time star formation histories, significantly different from halo accretion histories. Feedback also produces large short-time-scale variability in galactic SFRs, especially in dwarfs. These properties are not captured by common ‘sub-grid’ wind models.

Nonthermal Electrons in Radiatively Inefficient Accretion Flow Models of Sagittarius A

We investigate radiatively inefficient accretion flow models for Sgr A*, the supermassive black hole in our Galactic center, in light of new observational constraints. Confirmation of linear polarization in the submillimeter emission argues for accretion rates much less than the canonical Bondi rate. We consider models with low accretion rates and calculate the spectra produced by a hybrid electron population consisting of both thermal and nonthermal particles. The thermal electrons produce the submillimeter emission and can account for its linear polarization properties. As noted in previous work, the observed low-frequency radio spectrum can be explained if a small fraction (≈1.5%) of the electron thermal energy resides in a soft power-law tail. In the innermost region of the accretion flow, turbulence and/or magnetic reconnection events may occasionally accelerate a fraction of the electrons into a harder power-law tail. We show that the synchrotron emission from these electrons, or the Compton upscattering of synchrotron photons by the same electrons, may account for the X-ray flares observed by Chandra.

ASTROPHYSICAL GYROKINETICS: KINETIC AND FLUID TURBULENT CASCADES IN MAGNETIZED WEAKLY COLLISIONAL PLASMAS

This paper presents a theoretical framework for understanding plasma turbulence in astrophysical plasmas. It is motivated by observations of electromagnetic and density fluctuations in the solar wind, interstellar medium and galaxy clusters, as well as by models of particle heating in accretion disks. All of these plasmas and many others have turbulent motions at weakly collisional and collisionless scales. The paper focuses on turbulence in a strong mean magnetic field. The key assumptions are that the turbulent fluctuations are small compared to the mean field, spatially anisotropic with respect to it and that their frequency is low compared to the ion cyclotron frequency. The turbulence is assumed to be forced at some system-specific outer scale. The energy injected at this scale has to be dissipated into heat, which ultimately cannot be accomplished without collisions. A kinetic cascade develops that brings the energy to collisional scales both in space and velocity. The nature of the kinetic cascade in various scale ranges depends on the physics of plasma fluctuations that exist there. There are four special scales that separate physically distinct regimes: the electron and ion gyroscales, the mean free path and the electron diffusion scale. In each of the scale ranges separated by these scales, the fully kinetic problem is systematically reduced to a more physically transparent and computationally tractable system of equations, which are derived in a rigorous way. In the inertial range above the ion gyroscale, the kinetic cascade separates into two parts: a cascade of Alfvenic fluctuations and a passive cascade of density and magnetic-field-strength fluctuations. The former are governed by the reduced magnetohydrodynamic (RMHD) equations at both the collisional and collisionless scales; the latter obey a linear kinetic equation along the (moving) field lines associated with the Alfvenic component (in the collisional limit, these compressive fluctuations become the slow and entropy modes of the conventional MHD). In the dissipation range below ion gyroscale, there are again two cascades: the kinetic-Alfven-wave (KAW) cascade governed by two fluid-like electron reduced magnetohydrodynamic (ERMHD) equations and a passive cascade of ion entropy fluctuations both in space and velocity. The latter cascade brings the energy of the inertial-range fluctuations that was Landau-damped at the ion gyroscale to collisional scales in the phase space and leads to ion heating. The KAW energy is similarly damped at the electron gyroscale and converted into electron heat. Kolmogorov-style scaling relations are derived for all of these cascades. The relationship between the theoretical models proposed in this paper and astrophysical applications and observations is discussed in detail.This paper presents a theoretical framework for understand ing plasma turbulence in astrophysical plasmas. It is motivated by observations of electromagnetic and densit y fluctuations in the solar wind, interstellar medium and galaxy clusters, as well as by models of particle heating in accretion disks. All of these plasmas and many others have turbulent motions at weakly collisional an d collisionless scales. This paper focuses on turbulence in a strong mean magnetic field (the guide field). T he key assumptions behind the theory developed here are that the turbulent fluctuations are anisotropic wit h respect to the mean field and that their frequency is low compared to the ion cyclotron frequency. The turbulen ce is assumed to be stirred (forced) at some system-specific outer scale L. The energy injected at this scale has to be dissipated into h eat, which ultimately cannot be accomplished without collisions. A kinetic cascadedevelops that brings the energy to collisional scales both in space and velocity. The nature of the kinetic c as ade in various scale ranges depends on the physics of plasma fluctuations that can exist there. There ar four special scales that separate physically distinct regimes: the electron gyroscale ρe, the ion gyroscaleρi , the mean free path λmfpi and the electron heat diffusion scale (mi/me)λmfpi (me andmi are electron and ion masses). In each of the scale ranges sepa rated by these scales, a number of physically meaningful and rigorously ju stifiable simplifications of the fully kinetic plasma description are possible. These are derived systematicall y vi a hierarchy of asymptotic expansions. The result is that, in each scale range, the fully kinetic proble m is reduced to a more physically transparent and computationally tractable system of equations, which are d erived in a rigorous way. In the “inertial range” above the ion gyroscale, the kinetic cascade separates into two parts: a cascade of Alfvénic fluctuations and a passive cascade of density and magnetic-field-strength fluc tuations. The former are governed by two fluid-like Reduced Magnetohydrodynamic (RMHD) equations at both the c ollisional and collisionless scales; the latter obey a linear kinetic equation along the (moving) field lines as ociated with the Alfvénic component (in the collisional limit, these passive fluctuations become the sl ow and entropy modes of the conventional MHD). In the “dissipation range” between the ion and electron gyroscales, there are again two cascades: the kineticAlfvén-wave (KAW) cascade governed by two fluid-like Electr on Reduced Magnetohydrodynamic (ERMHD) equations and a passive cascade of ion entropy fluctuations b oth in space and velocity. The latter cascade brings the energy of the inertial-range fluctuations that was dampe d by collisionless wave-particle interaction at the ion gyroscle to collisional scales and leads to ion heating. The KAW energy is similarly damped at the electron gyroscale and is converted into electron heat. The relation ship between the theoretical models proposed in this paper and astrophysical applications and observations is d i cussed in detail. Subject headings: magnetic fields—methods: analytical—MHD—plasmas—turbul ence

A Possible Relativistic Jetted Outburst from a Massive Black Hole Fed by a Tidally Disrupted Star

A recent bright emission observed by the Swift satellite is due to the sudden accretion of a star onto a massive black hole. Gas accretion onto some massive black holes (MBHs) at the centers of galaxies actively powers luminous emission, but most MBHs are considered dormant. Occasionally, a star passing too near an MBH is torn apart by gravitational forces, leading to a bright tidal disruption flare (TDF). Although the high-energy transient Sw 1644+57 initially displayed none of the theoretically anticipated (nor previously observed) TDF characteristics, we show that observations suggest a sudden accretion event onto a central MBH of mass about 106 to 107 solar masses. There is evidence for a mildly relativistic outflow, jet collimation, and a spectrum characterized by synchrotron and inverse Compton processes; this leads to a natural analogy of Sw 1644+57 to a temporary smaller-scale blazar.

CONVECTION-DOMINATED ACCRETION FLOWS

Nonradiating advection-dominated accretion flows are convectively unstable in the radial direction. We calculate the two-dimensional (r-θ) structure of such flows assuming that (1) convection transports angular momentum inward, opposite to normal viscosity, and (2) viscous transport by other mechanisms (e.g., magnetic fields) is weak (α 1). Under such conditions convection dominates the dynamics of the accretion flow and leads to a steady state structure that is marginally stable to convection. We show that the marginally stable flow has a constant temperature and rotational velocity on spherical shells, a net flux of energy from small to large radii, zero net accretion rate, and a radial density profile of ρ ∝ r-1/2, flatter than the ρ ∝ r-3/2 profile characteristic of spherical accretion flows. This solution accurately describes the full two-dimensional structure of recent axisymmetric numerical simulations of advection-dominated accretion flows.

Stellar feedback in galaxies and the origin of galaxy‐scale winds

Feedback from massive stars is believed to play a critical role in driving galactic superwinds that enrich the intergalactic medium and shape the galaxy mass function, massmetallicity relation, and other global galaxy properties. In previous papers, we have introduced new numerical methods for implementing stellar feedback on sub-GMC through galactic scales in numerical simulations of galaxies; the key physical processes include radiation pressure in the UV through IR, supernovae (Type-I & II), stellar winds (“fast” O star through “slow” AGB winds), and HII photoionization. Here, we show that these feedback mechanisms drive galactic winds with outflow rates as high as 10 20 times the galaxy star formation rate. The mass-loading efficiency (wind mass loss rate divided by the star formation rate) scales roughly as _ Mwind= _ M / V 1 c (where Vc is the galaxy circular velocity), consistent with simple momentum-conservation expectations. We use our suite of simulations to study the relative contribution of each feedback mechanism to the generation of galactic winds in a range of galaxy models, from SMC-like dwarfs and Milky-way analogues to z 2 clumpy disks. In massive, gas-rich systems (local starbursts and high-z galaxies), radiation pressure dominates the wind generation. By contrast, for MW-like spirals and dwarf galaxies the gas densities are much lower and sources of shock-heated gas such as supernovae and stellar winds dominate the production of large-scale outflows. In all of our models, however, the winds have a complex multi-phase structure that depends on the interaction between multiple feedback mechanisms operating on different spatial and time scales: any single feedback mechanism fails to reproduce the winds observed. We use our simulations to provide fitting functions to the wind mass-loading and velocities as a function of galaxy properties, for use in cosmological simulations and semi-analytic models. These differ from typically-adopted formulae with an explicit dependence on the gas surface density that can be very important in both low-density dwarf galaxies and high-density gas-rich galaxies.

Dynamical friction and galaxy merging time‐scales

The time-scale for galaxies within merging dark matter haloes to merge with each other is an important ingredient in galaxy formation models. Accurate estimates of merging time-scales are required for predictions of astrophysical quantities such as black hole binary merger rates, the build-up of stellar mass in central galaxies and the statistical properties of satellite galaxies within dark matter haloes. In this paper, we study the merging time-scales of extended dark matter haloes using N-body simulations. We compare these results to standard estimates based on the Chandrasekhar theory of dynamical friction. We find that these standard predictions for merging time-scales, which are often used in semi-analytic galaxy formation models, are systematically shorter than those found in simulations. The discrepancy is approximately a factor of 1.7 for M sat /M host ≈ 0.1 and becomes larger for more disparate satellite-to-host mass ratios, reaching a factor of ∼3.3 for M sat /M host ≈ 0.01. Based on our simulations, we propose a new, easily implementable fitting formula that accurately predicts the time-scale for an extended satellite to sink from the virial radius of a host halo down to the halos centre for a wide range of M sat /M host and orbits. Including a central bulge in each galaxy changes the merging time-scale by ≤10 per cent. To highlight one concrete application of our results, we show that merging time-scales often used in the literature overestimate the growth of stellar mass by satellite accretion by ≈40 per cent, with the extra mass gained in low mass ratio mergers.

Magnetic Fluctuation Power Near Proton Temperature Anisotropy Instability Thresholds in the Solar Wind

The proton temperature anisotropy in the solar wind is known to be constrained by the theoretical thresholds for pressure-anisotropy-driven instabilities. Here, we use approximately 1x10;{6} independent measurements of gyroscale magnetic fluctuations in the solar wind to show for the first time that these fluctuations are enhanced along the temperature anisotropy thresholds of the mirror, proton oblique firehose, and ion cyclotron instabilities. In addition, the measured magnetic compressibility is enhanced at high plasma beta (beta_{ parallel} greater, similar1) along the mirror instability threshold but small elsewhere, consistent with expectations of the mirror mode. We also show that the short wavelength magnetic fluctuation power is a strong function of collisionality, which relaxes the temperature anisotropy away from the instability conditions and reduces correspondingly the fluctuation power.

Optical Flares from the Tidal Disruption of Stars by Massive Black Holes

A star that wanders too close to a massive black hole (BH) is shredded by the BH’s tidal gravity. Stellar gas falls back to the BH at a rate initially exceeding the Eddington rate, releasing a flare of energy. In anticipation of upcoming transient surveys, we predict the light curves and spectra of tidal flares as a function of time, highlighting the unique signatures of tidal flares at optical and near-infrared wavelengths. A reasonable fraction of the gas initially bound to the BH is likely blown away when the fallback rate is super-Eddington at early times. This outflow produces an optical luminosity comparable to that of a supernova; such events have durations of � 10 days and may have been missed in supernova searches that exclude the nuclear regions of galaxies. When the fallback rate subsides below Eddington, the gas accretes onto the BH via a thin disk whose emission peaks in the UV to soft X-rays. Some of this emission is reprocessed by the unbound stellar debris, producing a spectrum of very broad emission lines (with no corresponding narrow forbidden lines). These lines are the strongest for BHs with MBH � 10 5 10 6 M⊙ and thus optical surveys are particularly sensitive to the lowest mass BHs in galactic nuclei. Calibrating our models to ROSAT and GALEX observations, we predict detection rates for Pan-STARRS, PTF, and LSST and highlight some of the observational challenges associated with studying tidal disruption events in the optical. Upcoming surveys such as Pan-STARRS should detect at least tens of events per year, and may detect many more if current models of outflows during super-Eddington accretion are reasonably accurate. These surveys will significantly improve our knowledge of stellar dynamics in galactic nuclei, the physics of super-Eddington accretion, the demography of intermediate mass BHs, and the role of tidal disruption in the growth of massive BHs.

A Model of Turbulence in Magnetized Plasmas: Implications for the Dissipation Range in the Solar Wind

This paper studies the turbulent cascade of magnetic energy in weakly col- lisional magnetized plasmas. A cascade model is presented, based on the assumptions of local nonlinear energy transfer in wavenumber space, critical balance between linear propagation and nonlinear interaction times, and the applicability of linear dissipation rates for the nonlinearly turbulent plasma. The model follows the nonlinear cascade of energy from the driving scale in the MHD regime, through the transition at the ion Lar- mor radius into the kinetic Alfven wave regime, in which the turbulence is dissipated by kinetic processes. The turbulent fluctuations remain at frequencies below the ion cy- clotron frequency due to the strong anisotropy of the turbulent fluctuations, kk ≪ k⊥ (implied by critical balance). In this limit, the turbulence is optimally described by gy- rokinetics; it is shown that the gyrokinetic approximation is well satisfied for typical slow solar wind parameters. Wave phase velocity measurements are consistent with a kinetic Alfven wave cascade and not the onset of ion cyclotron damping. The conditions under which the gyrokinetic cascade reaches the ion cyclotron frequency are established. Cas- cade model solutions imply that collisionless damping provides a natural explanation for the observed range of spectral indices in the dissipation range of the solar wind. The dis- sipation range spectrum is predicted to be an exponential fall off; the power-law behav- ior apparent in observations may be an artifact of limited instrumental sensitivity. The cascade model is motivated by a programme of gyrokinetic simulations of turbulence and particle heating in the solar wind.