Source author record

Eliot Quataert

Eliot Quataert appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

162works
15topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

162 published item(s)

preprint2022arXiv

Cosmic-Ray Driven Outflows to Mpc Scales from $L_{\ast}$ Galaxies

We study the effects of cosmic rays (CRs) on outflows from star-forming galaxies in the circum and inter-galactic medium (CGM/IGM), in high-resolution, fully-cosmological FIRE-2 simulations (accounting for mechanical and radiative stellar feedback, magnetic fields, anisotropic conduction/viscosity/CR diffusion and streaming, and CR losses). We showed previously that massive ($M_{\rm halo}\gtrsim 10^{11}\,M_{\odot}$), low-redshift ($z\lesssim 1-2$) halos can have CR pressure dominate over thermal CGM pressure and balance gravity, giving rise to a cooler CGM with an equilibrium density profile. This dramatically alters outflows. Absent CRs, high gas thermal pressure in massive halos "traps" galactic outflows near the disk, so they recycle. With CRs injected in supernovae as modeled here, the low-pressure halo allows "escape" and CR pressure gradients continuously accelerate this material well into the IGM in "fast" outflows, while lower-density gas at large radii is accelerated in-situ into "slow" outflows that extend to $>$Mpc scales. CGM/IGM outflow morphologies are radically altered: they become mostly volume-filling (with inflow in a thin mid-plane layer) and coherently biconical from the disk to $>$Mpc. The CR-driven outflows are primarily cool ($T\sim10^{5}\,$K) and low-velocity. All of these effects weaken and eventually vanish at lower halo masses ($\lesssim 10^{11}\,M_{\odot}$) or higher redshifts ($z\gtrsim 1-2$), reflecting the ratio of CR to thermal+gravitational pressure in the outer halo. We present a simple analytic model which explains all of the above phenomena.

preprint2022arXiv

Effects of Different Cosmic Ray Transport Models on Galaxy Formation

Cosmic rays (CRs) with ~GeV energies can contribute significantly to the energy and pressure budget in the interstellar, circumgalactic, and intergalactic medium (ISM, CGM, IGM). Recent cosmological simulations have begun to explore these effects, but almost all studies have been restricted to simplified models with constant CR diffusivity and/or streaming speeds. Physical models of CR propagation/scattering via extrinsic turbulence and self-excited waves predict transport coefficients which are complicated functions of local plasma properties. In a companion paper, we consider a wide range of observational constraints to identify proposed physically-motivated cosmic-ray propagation scalings which satisfy both detailed Milky Way (MW) and extra-galactic $γ$-ray constraints. Here, we compare the effects of these models relative to simpler 'diffusion+streaming' models on galaxy and CGM properties at dwarf through MW mass scales. The physical models predict large local variations in CR diffusivity, with median diffusivity increasing with galacto-centric radii and decreasing with galaxy mass and redshift. These effects lead to a more rapid dropoff of CR energy density in the CGM (compared to simpler models), in turn producing weaker effects of CRs on galaxy star formation rates (SFRs), CGM absorption profiles and galactic outflows. The predictions of the more physical CR models tend to lie 'in between' models which ignore CRs entirely and models which treat CRs with constant diffusivity.

preprint2022arXiv

Galaxies lacking dark matter produced by close encounters in a cosmological simulation

The standard cold dark matter plus cosmological constant model predicts that galaxies form within dark-matter haloes, and that low-mass galaxies are more dark-matter dominated than massive ones. The unexpected discovery of two low-mass galaxies lacking dark matter immediately provoked concerns about the standard cosmology and ignited explorations of alternatives, including self-interacting dark matter and modified gravity. Apprehension grew after several cosmological simulations using the conventional model failed to form adequate numerical analogues with comparable internal characteristics (stellar masses, sizes, velocity dispersions and morphologies). Here we show that the standard paradigm naturally produces galaxies lacking dark matter with internal characteristics in agreement with observations. Using a state-of-the-art cosmological simulation and a meticulous galaxy-identification technique, we find that extreme close encounters with massive neighbours can be responsible for this. We predict that approximately 30 percent of massive central galaxies (with at least 1e11 solar masses in stars) harbour at least one dark-matter-deficient satellite (with 1e8 - 1e9 solar masses in stars). This distinctive class of galaxies provides an additional layer in our understanding of the role of interactions in shaping galactic properties. Future observations surveying galaxies in the aforementioned regime will provide a crucial test of this scenario.

preprint2022arXiv

Gas infall and radial transport in cosmological simulations of Milky Way-mass disks

Observations indicate that a continuous supply of gas is needed to maintain observed star formation rates in large, disky galaxies. To fuel star formation, gas must reach the inner regions of such galaxies. Despite its crucial importance for galaxy evolution, how and where gas joins galaxies is poorly constrained observationally and is rarely explored in fully cosmological simulations. To investigate gas accretion in the vicinity of galaxies, we analyze the FIRE-2 cosmological zoom-in simulations for 4 Milky Way mass galaxies (M_halo ~ 10E12 solar masses), focusing on simulations with cosmic ray physics. We find that at z~0, gas approaches the disk with angular momentum similar to the gaseous disk edge and low radial velocities, piling-up near the edge and settling into full rotational support. Accreting gas moves predominantly parallel to the disk with small but nonzero vertical velocity components, and joins the disk largely in the outskirts as opposed to "raining" down onto the disk. Once in the disk, gas trajectories are complex, being dominated by spiral arm induced oscillations and feedback. However, time and azimuthal averages show clear but slow net radial infall with transport speeds of 1-3 km/s and net mass fluxes through the disk on the order of one solar mass per year, comparable to the star formation rates of the galaxies and decreasing towards galactic center as gas is sunk into star formation. These rates are slightly higher in simulations without cosmic rays (1-7 km/s, ~4-5 solar masses per year). We find overall consistency of our results with observational constraints and discuss prospects of future observations of gas flows in and around galaxies.

preprint2022arXiv

High-frequency heating of the solar wind triggered by low-frequency turbulence

The fast solar wind's high speeds and nonthermal features require that significant heating occurs well above the Sun's surface. Two leading theories seem incompatible: low-frequency "Alfvénic" turbulence, which transports energy outwards and is observed ubiquitously by spacecraft but struggles to explain the observed dominance of ion over electron heating; and high-frequency ion-cyclotron waves (ICWs), which explain the nonthermal heating of ions but lack an obvious source. Here, we argue that the recently proposed "helicity barrier" effect, which limits electron heating by inhibiting the turbulent cascade of energy to the smallest scales, can unify these two paradigms. Our six-dimensional simulations show how the helicity barrier causes the large-scale energy to grow in time, generating small parallel scales and high-frequency ICW heating from low-frequency turbulence. The resulting turbulence and ion distribution function also closely match in-situ measurements from Parker Solar Probe and other spacecraft, explaining, among other features, the decades-long puzzle of the steep "transition range" observed in magnetic fluctuation spectra. The theory predicts a causal link between plasma expansion and the ion-to-electron heating ratio. Given the observational association between wind speed and expansion, we argue that the helicity barrier could play a key role in regulating the bimodal speed distribution of the solar wind.

preprint2022arXiv

Hot-mode accretion and the physics of thin-disk galaxy formation

We use FIRE simulations to study disk formation in z~0, Milky Way-mass galaxies, and conclude that a key ingredient for the formation of thin stellar disks is the ability for accreting gas to develop an aligned angular momentum distribution via internal cancellation *prior* to joining the galaxy. Among galaxies with a high fraction (>70%) of their young stars in a thin disk (h/R~0.1) we find that: (i) hot, virial-temperature gas dominates the inflowing gas mass on halo scales (>~20 kpc), with radiative losses offset by compression heating; (ii) this hot accretion proceeds until angular momentum support slows inward motion, at which point the gas cools to T~10^4 K or less; (iii) prior to cooling, the accreting gas develops an angular momentum distribution that is aligned with the galaxy disk, and while cooling transitions from a quasi-spherical spatial configuration to a more flattened, disk-like configuration. We show that the existence of this "rotating cooling flow" accretion mode is strongly correlated with the fraction of stars forming in a thin disk among a sample of 17 z~0 galaxies spanning a halo mass range of 10^10.5 solar masses to 10^12 solar masses, or a stellar mass range 10^8 solar masses to 10^11 solar masses. Notably, galaxies with a thick disk or irregular morphology do not undergo significant angular momentum alignment of gas prior to accretion and show no correspondence between halo gas cooling and flattening. Our results suggest that rotating cooling flows (or, more generally, rotating subsonic flows) that become coherent and angular momentum-supported prior to accretion onto the galaxy are likely a necessary condition for the formation of thin, star-forming disk galaxies in a LambdaCDM universe.

preprint2022arXiv

Long-term evolution of supercritical black hole accretion with outflows: a subgrid feedback model for cosmological simulations

We study the long-term evolution of the global structure of axisymmetric accretion flows onto a black hole (BH) at rates substantially higher than the Eddington value ($\dot{M}_{\rm Edd}$), performing two-dimensional hydrodynamical simulations with and without radiative diffusion. In the high-accretion optically-thick limit, where the radiation energy is efficiently trapped within the inflow, the accretion flow becomes adiabatic and comprises of turbulent gas in the equatorial region and strong bipolar outflows. As a result, the mass inflow rate decreases toward the center as $\dot{M}_{\rm in}\propto r^{p}$ with $p\sim 0.5-0.7$ and a small fraction of the inflowing gas feeds the nuclear BH. Thus, super-Eddington accretion is sustained only when a larger amount of gas is supplied from larger radii at $> 100-1000~\dot{M}_{\rm Edd}$. The global structure of the flow settles down to a quasi-steady state in millions of the orbital timescale at the BH event horizon, which is $> 10-100$ times longer than that addressed in previous (magneto-)RHD simulation studies. Energy transport via radiative diffusion accelerates the outflow near the poles in the inner region but does not change the overall properties of the accretion flow compared to the cases without diffusion. Based on our simulation results, we provide a mechanical feedback model for super-Eddington accreting BHs. This can be applied as a sub-grid model in large-scale cosmological simulations that do not sufficiently resolve galactic nuclei, and to the formation of the heaviest gravitational-wave sources via accretion in dense environments.

preprint2022arXiv

Numerical Simulations of the Random Angular Momentum in Convection: Implications for Supergiant Collapse to Form Black Holes

During the core collapse of massive stars that do not undergo a canonical energetic explosion, some of the hydrogen envelope of a red supergiant (RSG) progenitor may infall onto the newborn black hole (BH). Within the Athena++ framework, we perform three-dimensional, hydrodynamical simulations of idealized models of supergiant convection and collapse in order to assess whether the infall of the convective envelope can give rise to rotationally-supported material, even if the star has zero angular momentum overall. Our dimensionless, polytropic models are applicable to the optically-thick hydrogen envelope of non-rotating RSGs and cover a factor of 20 in stellar radius. At all radii, the specific angular momentum due to random convective flows implies associated circularization radii of 10 - 1500 times the innermost stable circular orbit of the BH. During collapse, the angular momentum vector of the convective flows is approximately conserved and is slowly varying on the timescale relevant to forming disks at small radii. Our results indicate that otherwise failed explosions of RSGs lead to the formation of rotationally-supported flows that are capable of driving outflows to large radii and powering observable transients. When the BH is able to accrete most of the hydrogen envelope, the final BH spin parameter is $\sim$ 0.5, even though the star is non-rotating. For fractional accretion of the envelope, the spin parameter is generally lower and never exceeds 0.8. We discuss the implications of our results for transients produced by RSG collapse to a black hole.

preprint2022arXiv

Observational Signatures of Black Hole Accretion: Rotating vs. Spherical Flows with Tilted Magnetic Fields

We study the observational signatures of magnetically arrested black hole accretion with non-rotating inflow onto a rotating black hole; we consider a range of angles between the black hole spin and the initial magnetic field orientation. We compare the results of our General Relativistic Magneto-Hydrodynamic simulations to more commonly used rotating initial conditions and to the Event Horizon Telescope (EHT) observations of M87. We find that the mm intensity images, polarization images, and synchrotron emission spectra are very similar among the different simulations when post-processed with the same electron temperature model; observational differences due to different electron temperature models are significantly larger than those due to the different realizations of magnetically arrested accretion. The orientation of the mm synchrotron polarization is particularly insensitive to the initial magnetic field orientation, the electron temperature model, and the rotation of the inflowing plasma. The largest difference among the simulations with different initial rotation and magnetic tilt is in the strength and stability of the jet; spherical inflow leads to kink-unstable jets. We discuss the implications of our results for current and future EHT observations and for theoretical models of event-horizon-scale black hole accretion.

preprint2022arXiv

Optical to X-ray Signatures of Dense Circumstellar Interaction in Core-Collapse Supernovae

Progenitors of core-collapse supernovae (SNe) can shed significant mass to circumstellar material (CSM) in the months--years preceding core-collapse. The ensuing SN explosion launches ejecta that may subsequently collide with this CSM, producing shocks that can power emission across the electromagnetic spectrum. In this work we explore the thermal signatures of dense CSM interaction, when the CSM density profile is truncated at some outer radius. CSM with optical depth $>c/v$ (where $v$ is the shock velocity) will produce primarily $\sim$blackbody optical/UV emission whereas lower optical-depth CSM will power bremsstrahlung X-ray emission. Focusing on the latter, we derive light-curves and spectra of the resulting X-ray transients, that include a detailed treatment of Comptonization. Due to strong photoelectric absorption, the X-ray light-curve is dominated by the `post-interaction' phase that occurs after the shock reaches the CSM truncation radius. We treat this regime here for the first time. Using these results, we present the phase-space of optical, UV, and X-ray transients as a function of CSM properties, and discuss detectability prospects. We find that ROSAT would not have been sensitive to CSM X-ray transients but that eROSITA is expected to detect many such events. Future wide-field UV missions such as ULTRASAT will dramatically enhance sensitivity to large optical-depth CSM configurations. Finally, we present a framework within which CSM properties may be directly inferred from observable features of X-ray transients. This can serve as an important tool for studying stellar mass loss using SN X-ray detections.

preprint2022arXiv

Reconciling cosmic-ray transport theory with phenomenological models motivated by Milky-Way data

Phenomenological models of cosmic-ray (CR) transport in the Milky Way (MW) can reproduce a wide range of observations assuming that CRs scatter off of magnetic-field fluctuations with spectrum $\propto k^{-δ}$ and $δ\sim [1.4,1.67]$. We study the extent to which such models can be reconciled with current microphysical theories of CR transport, specifically self-confinement due to the streaming instability and/or extrinsic turbulence due to a cascade of MHD fast modes. We first review why it is that on their own neither theory is compatible with observations. We then highlight that CR transport is a strong function of local plasma conditions in the multi-phase interstellar medium (ISM), and may be diffusive due to turbulence in some regions and streaming due to self-confinement in others. A multi-phase combination of scattering mechanisms can in principle reproduce the main trends in the proton spectrum and the boron-to-carbon ratio (B/C). However, models with a combination of scattering by self-excited waves and fast-mode turbulence require significant fine-tuning due to fast-mode damping, unlike phenomenological models that assume undamped Kolmogorov turbulence. The assumption that fast modes follow a weak cascade is also not well justified theoretically, as the weak cascade is suppressed by wave steepening and weak-shock dissipation even in subsonic turbulence. These issues suggest that there may be a significant theoretical gap in our understanding of MHD turbulence. We discuss a few topics at the frontier of MHD turbulence theory that bear on this (possible) gap and that may be relevant for CR scattering.

preprint2022arXiv

Research Opportunities in Plasma Astrophysics

Major scientific questions and research opportunities are described on 10 unprioritized plasma astrophysics topics: (1) magnetic reconnection, (2) collisionless shocks and particle acceleration, (3) waves and turbulence, (4) magnetic dynamos, (5) interface and shear instabilities, (6) angular momentum transport, (7) dusty plasmas, (8) radiative hydrodynamics, (9) relativistic, pair-dominated and strongly magnetized plasmas, (10) jets and outflows. Note that this is a conference report from a Workshop on Opportunities in Plasma Astrophysics (WOPA, https://w3.pppl.gov/conferences/2010/WOPA/) in January 2010, that attracted broad representation from the community and was supported by the U.S. Department of Energy, National Aeronautics and Space Administration, National Science Foundation, American Physical Society's Topical Group for Plasma Astrophysics and Division of Plasma Physics, and Center for Magnetic Self-Organization in Laboratory and Astrophysical Plasmas. Although there has been much planning and many developments in both science and infrastructure since the report was written, most of the motivation, priorities, problems and technical challenges discussed therein remain unaddressed and are relevant at the time of posting.

preprint2022arXiv

Supercritical growth pathway to overmassive black holes at cosmic dawn: coevolution with massive quasar hosts

Observations of the most luminous quasars at high redshifts ($z > 6$) have revealed that the largest supermassive black holes (SMBHs) at those epochs tend to be substantially overmassive relative to their host galaxies compared to the local relations, suggesting they experienced rapid early growth phases. We propose an assembly model for the SMBHs that end up in rare massive $\sim10^{12}~M_{\odot}$ host halos at $z \sim 6-7$, applying a kinetic feedback prescription for BHs accreting above the Eddington rate, provided by radiation hydrodynamic simulations for the long-term evolution of the accretion-flow structure. The large inflow rates into these halos during their assembly enable the formation of $>10^9~M_{\odot}$ SMBHs by $z \sim 6$, even starting from stellar-mass seeds at $z \sim 30$, and even in the presence of outflows that reduce the BH feeding rate, especially at early times. This mechanism also naturally yields a high BH-to-galaxy mass ratio of $> 0.01$ before the SMBH mass reaches $M_{\rm BH} > 10^9~M_{\odot}$ by $z \sim 6$. These fast-growing SMBH progenitors are bright enough to be detected by upcoming observations with the James Webb Space Telescope over a wide range of redshift ($7 < z < 15$), regardless of how they were seeded.

preprint2022arXiv

Testing Physical Models for Cosmic Ray Transport Coefficients on Galactic Scales: Self-Confinement and Extrinsic Turbulence at GeV Energies

The microphysics of ~GeV cosmic ray (CR) transport on galactic scales remain deeply uncertain, with almost all studies adopting simple prescriptions (e.g. constant-diffusivity). We explore different physically-motivated, anisotropic, dynamical CR transport scalings in high-resolution cosmological FIRE simulations of dwarf and ~$L_{\ast}$ galaxies where scattering rates vary with local plasma properties motivated by extrinsic turbulence (ET) or self-confinement (SC) scenarios, with varying assumptions about e.g. turbulent power spectra on un-resolved scales, Alfven-wave damping, etc. We self-consistently predict observables including $γ$-rays ($L_γ$), grammage, residence times, and CR energy densities to constrain the models. We demonstrate many non-linear dynamical effects (not captured in simpler models) tend to enhance confinement. For example, in multi-phase media, even allowing arbitrary fast transport in neutral gas does not substantially reduce CR residence times (or $L_γ$), as transport is rate-limited by the ionized WIM and 'inner CGM' gaseous halo ($10^{4}-10^{6}$ K gas within 10-30 kpc), and $L_γ$ can be dominated by trapping in small 'patches.' Most physical ET models contribute negligible scattering of ~1-10 GeV CRs, but it is crucial to account for anisotropy and damping (especially of fast modes) or else scattering rates would violate observations. We show that the most widely-assumed scalings for SC models produce excessive confinement by factors >100 in the WIM and inner CGM, where turbulent and Landau damping dominate. This suggests either a breakdown of quasi-linear theory used to derive the CR transport parameters in SC, or that other novel damping mechanisms dominate in intermediate-density ionized gas.

preprint2022arXiv

The In-situ Origins of Dwarf Stellar Outskirts in FIRE-2

Extended, old, and round stellar halos appear to be ubiquitous around high-mass dwarf galaxies ($10^{8.5}<M_\star/M_\odot<10^{9.6}$) in the observed universe. However, it is unlikely that these dwarfs have undergone a sufficient number of minor mergers to form stellar halos that are composed of predominantly accreted stars. Here, we demonstrate that FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulations are capable of producing dwarf galaxies with realistic structure, including both a thick disk and round stellar halo. Crucially, these stellar halos are formed in-situ, largely via the outward migration of disk stars. However, there also exists a large population of "non-disky" dwarfs in FIRE-2 that lack a well-defined disk/halo and do not resemble the observed dwarf population. These non-disky dwarfs tend to be either more gas poor or to have burstier recent star formation histories than the disky dwarfs, suggesting that star formation feedback may be preventing disk formation. Both classes of dwarfs underscore the power of a galaxy's intrinsic shape -- which is a direct quantification of the distribution of the galaxy's stellar content -- to interrogate the feedback implementation in simulated galaxies.

preprint2021arXiv

A stripped-companion origin for Be stars: clues from the putative black holes HR 6819 and LB-1

HR 6819 is a bright ($V=5.36$), blue star recently proposed to be a triple containing a detached black hole (BH). We show that the system is a binary and does not contain a BH. Using spectral decomposition, we disentangle the observed composite spectra into two components: a rapidly rotating Be star and a slowly rotating B star with low surface gravity $(\log g \approx 2.75)$. Both stars show periodic radial velocity (RV) variability, but the RV semi-amplitude of the B star's orbit is $K_{\rm B}= (62.7 \pm 1)\,\rm km\,s^{-1}$, while that of the Be star is only $K_{\rm Be} = (4.5\pm 2)\,\rm km\,s^{-1}$. This implies that the B star is less massive by at least a factor of 10. The surface abundances of the B star bear imprints of CNO burning. We argue that the B star is a bloated, recently stripped helium star with mass $\approx 0.5\,M_{\odot}$ that is currently contracting to become a hot subdwarf. The orbital motion of the Be star obviates the need for a BH to explain the B star's motion. A stripped-star model reproduces the observed luminosity of the system, while a normal star with the B star's temperature and gravity would be more than 10 times too luminous. HR 6819 and the binary LB-1 probably formed through similar channels. We use MESA models to investigate their evolutionary history, finding that they likely formed from intermediate-mass ($3-7\,M_{\odot}$) primaries stripped by slightly lower-mass secondaries and are progenitors to Be + sdOB binaries such as $ϕ$ Persei. The lifetime of their current evolutionary phase is on average $2\times 10^5$ years, of order half a percent of the total lifetime of the Be phase. This implies that many Be stars have hot subdwarf and white dwarf companions, and that a substantial fraction ($20-100\%$) of field Be stars form through accretion of material from a binary companion.

preprint2021arXiv

Exploring the epoch of hydrogen reionization using FRBs

We describe three different methods for exploring the hydrogen reionization epoch using fast radio bursts (FRBs) and provide arguments for the existence of FRBs at high redshift (z). The simplest way, observationally, is to determine the maximum dispersion measure (DM$_{\rm max}$) of FRBs for an ensemble that includes bursts during the reionization. The DM$_{\rm max}$ provides information regarding reionization much like the optical depth of the CMB to Thomson scattering does, and it has the potential to be more accurate than constraints from Planck, if DM$_{\rm max}$ can be measured to a precision better than 500 $\mbox{pc cm}^{-3}$. Another method is to measure redshifts of about 40 FRBs between z of 6-10 with$\sim10\%$ accuracy to obtain the average electron density in 4 different z-bins with $\sim4\%$ accuracy. These two methods don't require knowledge of the FRB luminosity function and its possible redshift evolution. Finally, we show that the reionization history is reflected in the number of FRBs per unit DM, given a fluence limited survey of FRBs that includes bursts during the reionization epoch; we show using FIRE simulations that the contributions to DM from the FRB host galaxy $\&$ CGM during the reionization era is a small fraction of the observed DM. This third method requires no redshift information but does require knowledge of the FRB luminosity function.

preprint2021arXiv

Magnetically Modified Spherical Accretion in GRMHD: Reconnection-Driven Convection and Jet Propagation

We present 3D general relativistic magnetohydrodynamic(GRMHD) simulations of zero angular momentum accretion around a rapidly rotating black hole, modified by the presence of initially uniform magnetic fields. We consider serveral angles between the magnetic field direction and the black hole spin. In the resulting flows, the midplane dynamics are governed by magnetic reconnection-driven turbulence in a magnetically arrested (or a nearly arrested) state. Electromagnetic jets with outflow efficiencies ~10-200% occupy the polar regions, reaching several hundred gravitational radii before they dissipate due to the kink instability. The jet directions fluctuate in time and can be tilted by as much as ~30 degrees with respect to black hole spin, but this tilt does not depend strongly on the tilt of the initial magnetic field. A jet forms even when there is no initial net vertical magnetic flux since turbulent, horizon-scale fluctuations can generate a net vertical field locally. Peak jet power is obtained for an initial magnetic field tilted by 40-80 degrees with respect to the black hole spin because this maximizes the amount of magnetic flux that can reach the black hole. These simulations may be a reasonable model for low luminosity black hole accretion flows such as Sgr A* or M87.

preprint2021arXiv

The Effects of Tilt on the Time Variability of Millimeter and Infrared Emission from Sagittarius A*

Using a combination of general-relativistic magnetohydrodynamics simulations and ray tracing of synchrotron emission, we study the effect of modest (24 degrees) misalignment between the black hole spin and plasma angular momentum, focusing on the variability of total flux, image centroids, and image sizes. We consider both millimeter and infrared (IR) observables motivated by Sagittarius A* (Sgr A*), though our results apply more generally to optically thin flows. For most quantities, tilted accretion is more variable, primarily due to a significantly hotter and denser "coronal" region well off the disk midplane. We find (1) a 150% increase in millimeter light curve variability when adding tilt to the flow; (2) the tilted image centroid in the millimeter shifts on a scale of 3.7 microarcseconds over 28 hours (5000 gravitational times) for some electron temperature models; (3) tilted disk image diameters in the millimeter can be 10% larger (52 versus 47 microarcseconds) than those of aligned disks at certain viewing angles; (4) the tilted models produce significant IR flux, similar to that seen in Sgr A*, with comparable or even greater variability than observed; (5) for some electron models, the tilted IR centroid moves by more than 50 microarcseconds over several hours, in a similar fashion to the centroid motion detected by the GRAVITY interferometer.

preprint2021arXiv

The impact of astrophysical dust grains on the confinement of cosmic rays

We argue that charged dust grains could significantly impact the confinement and transport of galactic cosmic rays. For sub-GeV to ~1000GeV cosmic rays, small-scale parallel Alfvén waves, which isotropize cosmic rays through gyro-resonant interactions, are also gyro-resonant with charged grains. If the dust is nearly stationary, as in the bulk of the interstellar medium, Alfvén waves are damped by dust. This will reduce the amplitude of Alfvén waves produced by the cosmic rays through the streaming instability, thus enhancing cosmic-ray transport. In well-ionized regions, the dust damping rate is larger by a factor of ~10 than other mechanisms that damp parallel Alfvén waves at the scales relevant for ~GeV cosmic rays, suggesting that dust could play a key role in regulating cosmic-ray transport. In astrophysical situations in which the dust moves through the gas with super-Alfvénic velocities, Alfvén waves are rendered unstable, which could directly scatter cosmic rays. This interaction has the potential to create a strong feedback mechanism where dust, driven through the gas by radiation pressure, then strongly enhances the confinement of cosmic rays, increasing their capacity to drive outflows. This mechanism may act in the circumgalactic medium around star-forming galaxies and active galactic nuclei.

preprint2020arXiv

Ab Initio Horizon-Scale Simulations of Magnetically Arrested Accretion in Sagittarius A* Fed by Stellar Winds

We present 3D general relativistic magnetohydrodynamic (GRMHD) simulations of the accretion flow surrounding Sagittarius A* that are initialized using larger-scale MHD simulations of the $\sim$ 30 Wolf--Rayet (WR) stellar winds in the Galactic center. The properties of the resulting accretion flow on horizon scales are set not by ad hoc initial conditions but by the observationally constrained properties of the WR winds with limited free parameters. For this initial study we assume a non-spinning black hole. Our simulations naturally produce a $\sim 10^{-8} M_\odot$ yr$^{-1}$ accretion rate, consistent with previous phenomenological estimates. We find that a magnetically arrested flow is formed by the continuous accretion of coherent magnetic field being fed from large radii. Near the event horizon, the magnetic field is so strong that it tilts the gas with respect to the initial angular momentum and concentrates the originally quasi-spherical flow to a narrow disk-like structure. We also present 230 GHz images calculated from our simulations where the inclination angle and physical accretion rate are not free parameters but are determined by the properties of the WR stellar winds. The image morphology is highly time variable. Linear polarization on horizon scales is coherent with weak internal Faraday rotation.

preprint2020arXiv

Black widow evolution: magnetic braking by an ablated wind

Black widows are close binary systems in which a millisecond pulsar is orbited by a companion a few per cent the mass of the sun. It has been suggested that the pulsar's rotationally powered $γ$-ray luminosity gradually evaporates the companion, eventually leaving behind an isolated millisecond pulsar. The evaporation efficiency is determined by the temperature $T_{\rm ch}\propto F^{2/3}$ to which the outflow is heated by the flux $F$ on a dynamical time-scale. Evaporation is most efficient for companions that fill their Roche lobes. In this case, the outflow is dominated by a cap around the L1 point with an angle $θ_g\sim (T_{\rm ch}/T_g)^{1/2}$, and the evaporation time is $t_{\rm evap}=0.46(T_{\rm ch}/T_g)^{-2}\textrm{ Gyr}$, where $T_g>T_{\rm ch}$ is the companion's virial temperature. We apply our model to the observed black widow population, which has increased substantially over the last decade, considering each system's orbital period, companion mass, and pulsar spin-down power. While the original (Fruchter et al. 1988) black widow evaporates its companion on a few Gyr time-scale, direct evaporation on its own is too weak to explain the overall population. We propose instead that the evaporative wind couples to the companion's magnetic field, removes angular momentum from the binary, and maintains stable Roche-lobe overflow. While a stronger wind carries more mass, it also reduces the Alfvén radius, making this indirect magnetic braking mechanism less dependent on the flux $t_{\rm mag}\propto t_{\rm evap}^{1/3}$. This reduces the scatter in evolution times of observed systems, thus better explaining the combined black widow and isolated millisecond pulsar populations.

preprint2020arXiv

But What About... Cosmic Rays, Magnetic Fields, Conduction, & Viscosity in Galaxy Formation

We present a suite of high-resolution cosmological simulations, using the FIRE-2 feedback physics together with explicit treatment of magnetic fields, anisotropic conduction and viscosity, and cosmic rays (CRs) injected by supernovae (including anisotropic diffusion, streaming, adiabatic, hadronic and Coulomb losses). We survey systems from ultra-faint dwarf ($M_{\ast}\sim 10^{4}\,M_{\odot}$, $M_{\rm halo}\sim 10^{9}\,M_{\odot}$) through Milky Way masses, systematically vary CR parameters (e.g. the diffusion coefficient $κ$ and streaming velocity), and study an ensemble of galaxy properties (masses, star formation histories, mass profiles, phase structure, morphologies). We confirm previous conclusions that magnetic fields, conduction, and viscosity on resolved ($\gtrsim 1\,$pc) scales have small effects on bulk galaxy properties. CRs have relatively weak effects on all galaxy properties studied in dwarfs ($M_{\ast} \ll 10^{10}\,M_{\odot}$, $M_{\rm halo} \lesssim 10^{11}\,M_{\odot}$), or at high redshifts ($z\gtrsim 1-2$), for any physically-reasonable parameters. However at higher masses ($M_{\rm halo} \gtrsim 10^{11}\,M_{\odot}$) and $z\lesssim 1-2$, CRs can suppress star formation by factors $\sim 2-4$, given relatively high effective diffusion coefficients $κ\gtrsim 3\times10^{29}\,{\rm cm^{2}\,s^{-1}}$. At lower $κ$, CRs take too long to escape dense star-forming gas and lose energy to hadronic collisions, producing negligible effects on galaxies and violating empirical constraints from $γ$-ray emission. But around $κ\sim 3\times10^{29}\,{\rm cm^{2}\,s^{-1}}$, CRs escape the galaxy and build up a CR-pressure-dominated halo which supports dense, cool ($T\ll 10^{6}$ K) gas that would otherwise rain onto the galaxy. CR heating (from collisional and streaming losses) is never dominant.

preprint2020arXiv

No missing photons for reionization: moderate ionizing photon escape fractions from the FIRE-2 simulations

We present the escape fraction of hydrogen ionizing photons (f_esc) from a sample of 34 high-resolution cosmological zoom-in simulations of galaxies at z>5 in the Feedback in Realistic Environments project, post-processed with a Monte Carlo radiative transfer code for ionizing radiation. Our sample consists of 8500 halos in M_vir~10^8--10^{12} M_sun (M_star~10^4--10^{10} M_sun) at z=5--12. We find the sample average <f_esc> increases with halo mass for M_vir~10^8--10^{9.5} M_sun, becomes nearly constant for M_vir~10^{9.5}--10^{11} M_sun, and decreases at M_vir>10^{11} M_sun. Equivalently, <f_esc> increases with stellar mass up to M_star~10^8 M_sun and decreases at higher masses. Even applying single-star stellar population synthesis models, we find a moderate <f_esc>~0.2 for galaxies at M_star~10^8 M_sun. Nearly half of the escaped ionizing photons come from stars 1--3 Myr old and the rest from stars 3--10 Myr old. Binaries only have a modest effect, boosting <f_esc> by ~25--35% and the number of escaped photons by 60--80%. Most leaked ionizing photons are from vigorously star-forming regions that usually contain a feedback-driven kpc-scale superbubble surrounded by a dense shell. The shell is forming stars while accelerated, so new stars formed earlier in the shell are already inside the shell. Young stars in the bubble and near the edge of the shell can fully ionize some low-column-density paths pre-cleared by feedback, allowing a large fraction of their ionizing photons to escape. The decrease of <f_esc> at the high-mass end is due to dust attenuation, while at the low-mass end, <f_esc> decreases owing to inefficient star formation (and hence feedback). At fixed mass, <f_esc> tends to increase with redshift. Our simulations produce sufficient ionizing photons for cosmic reionization.

preprint2020arXiv

Not so fast: LB-1 is unlikely to contain a 70 $M_{\odot}$ black hole

The recently discovered binary LB-1 has been reported to contain a $\sim$\,$70\,M_{\odot}$ black hole (BH). The evidence for the unprecedentedly high mass of the unseen companion comes from reported radial velocity (RV) variability of the H$α$ emission line, which has been proposed to originate from an accretion disk around a BH. We show that there is in fact no evidence for RV variability of the H$α$ emission line, and that its apparent shifts instead originate from shifts in the luminous star's H$α$ absorption line. If not accounted for, such shifts will cause a stationary emission line to appear to shift in anti-phase with the luminous star. We show that once the template spectrum of a B star is subtracted from the observed Keck/HIRES spectra of LB-1, evidence for RV variability vanishes. Indeed, the data rule out periodic variability of the line with velocity semi-amplitude $K_{\rm Hα} > 1.3\,\rm km\,s^{-1}$. This strongly suggests that the observed H$α$ emission does not originate primarily from an accretion disk around a BH, and thus that the mass ratio cannot be constrained from the relative velocity amplitudes of the emission and absorption lines. The nature of the unseen companion remains uncertain, but a "normal" stellar-mass BH with mass $5\lesssim M/M_{\odot}\lesssim 20 $ seems most plausible. The H$α$ emission likely originates primarily from circumbinary material, not from either component of the binary.

preprint2020arXiv

On the Comparison of AGN with GRMHD Simulations: I. Sgr A*

We present models of Galactic Center emission in the vicinity of Sagittarius A* that use parametrizations of the electron temperature or energy density. These models include those inspired by two-temperature general relativistic magnetohydrodynamic (GRMHD) simulations as well as jet-motivated prescriptions generalizing equipartition of particle and magnetic energies. From these models, we calculate spectra and images and classify them according to their distinct observational features. Some models produce morphological and spectral features, e.g., image sizes, the sub-mm bump and low frequency spectral slope compatible with observations. Models with spectra consistent with observations produce the most compact images, with the most prominent, asymmetric photon rings. Limb brightened outflows are also visible in many models. Of all the models we consider, that which represents the current data the best is one in which electrons are relativistically hot when magnetic pressure is larger than the thermal pressure, but cold (i.e., negligibly contributing to the emission) otherwise. This work is part of a series also applying the "observing" simulations methodology to near-horizon regions of supermassive black holes in M87 and 3C 279.

preprint2020arXiv

Properties of the Circumgalactic Medium in Cosmic Ray-Dominated Galaxy Halos

We investigate the impact of cosmic rays (CRs) on the circumgalactic medium (CGM) in FIRE-2 simulations, for ultra-faint dwarf through Milky Way (MW)-mass halos hosting star-forming (SF) galaxies. Our CR treatment includes injection by supernovae, anisotropic streaming and diffusion along magnetic field lines, collisional and streaming losses, with constant parallel diffusivity $κ\sim3\times10^{29}\,\mathrm{cm^2\ s^{-1}}$ chosen to match $γ$-ray observations. With this, CRs become more important at larger halo masses and lower redshifts, and dominate the pressure in the CGM in MW-mass halos at $z\lesssim 1-2$. The gas in these "CR-dominated" halos differs significantly from runs without CRs: the gas is primarily cool (a few $\sim10^{4}\,$K), and the cool phase is volume-filling and has a thermal pressure below that needed for virial or local thermal pressure balance. Ionization of the "low" and "mid" ions in this diffuse cool gas is dominated by photo-ionization, with O VI columns $\gtrsim 10^{14.5}\,\mathrm{cm^{-2}}$ at distances $\gtrsim 150\,\mathrm{kpc}$. CR and thermal gas pressure are locally anti-correlated, maintaining total pressure balance, and the CGM gas density profile is determined by the balance of CR pressure gradients and gravity. Neglecting CRs, the same halos are primarily warm/hot ($T\gtrsim 10^{5}\,$K) with thermal pressure balancing gravity, collisional ionization dominates, O VI columns are lower and Ne VIII higher, and the cool phase is confined to dense filaments in local thermal pressure equilibrium with the hot phase.

preprint2020arXiv

Self-consistent proto-globular cluster formation in cosmological simulations of high-redshift galaxies

We report the formation of bound star clusters in a sample of high-resolution cosmological zoom-in simulations of z>5 galaxies from the FIRE project. We find that bound clusters preferentially form in high-pressure clouds with gas surface densities over 10^4 Msun pc^-2, where the cloud-scale star formation efficiency is near unity and young stars born in these regions are gravitationally bound at birth. These high-pressure clouds are compressed by feedback-driven winds and/or collisions of smaller clouds/gas streams in highly gas-rich, turbulent environments. The newly formed clusters follow a power-law mass function of dN/dM~M^-2. The cluster formation efficiency is similar across galaxies with stellar masses of ~10^7-10^10 Msun at z>5. The age spread of cluster stars is typically a few Myrs and increases with cluster mass. The metallicity dispersion of cluster members is ~0.08 dex in [Z/H] and does not depend on cluster mass significantly. Our findings support the scenario that present-day old globular clusters (GCs) were formed during relatively normal star formation in high-redshift galaxies. Simulations with a stricter/looser star formation model form a factor of a few more/fewer bound clusters per stellar mass formed, while the shape of the mass function is unchanged. Simulations with a lower local star formation efficiency form more stars in bound clusters. The simulated clusters are larger than observed GCs due to finite resolution. Our simulations are among the first cosmological simulations that form bound clusters self-consistently in a wide range of high-redshift galaxies.

preprint2020arXiv

Sound-Wave Instabilities in Dilute Plasmas with Cosmic Rays: Implications for Cosmic-Ray Confinement and the Perseus X-ray Ripples

Weakly collisional, magnetised plasmas characterised by anisotropic viscosity and conduction are ubiquitous in galaxies, halos and the intracluster medium (ICM). Cosmic rays (CRs) play an important role in these environments as well, by providing additional pressure and heating to the thermal plasma. We carry out a linear stability analysis of weakly collisional plasmas with cosmic rays using Braginskii MHD for the thermal gas. We assume that the CRs stream at the Alfvén speed, which in a weakly collisional plasma depends on the pressure anisotropy ($Δp$) of the thermal plasma. We find that this $Δp$-dependence introduces a phase shift between the CR-pressure and gas-density fluctuations. This drives a fast-growing acoustic instability: CRs offset the damping of acoustic waves by anisotropic viscosity and give rise to wave growth when the ratio of CR pressure to gas pressure is $\gtrsim αβ^{-1/2}$, where $β$ is the ratio of thermal to magnetic pressure, and $α$, typically $\lesssim 1$, depends on other dimensionless parameters. In high-$β$ environments like the ICM, this condition is satisfied for small CR pressures. We speculate that the instability studied here may contribute to the scattering of high-energy CRs and to the excitation of sound waves in galaxy-halo, group and cluster plasmas, including the long-wavelength X-ray fluctuations in \textit{Chandra} observations of the Perseus cluster. It may also be important in the vicinity of shocks in dilute plasmas (e.g., cluster virial shocks or galactic wind termination shocks), where the CR pressure is locally enhanced.

preprint2020arXiv

Synthetic Gaia surveys from the FIRE cosmological simulations of Milky Way-mass galaxies

With Gaia Data Release 2, the astronomical community is entering a new era of multidimensional surveys of the Milky Way. This new phase-space view of our Galaxy demands new tools for comparing observations to simulations of Milky-Way-mass galaxies in a cosmological context, to test the physics of both dark matter and galaxy formation. We present ananke, a framework for generating synthetic phase-space surveys from high-resolution baryonic simulations, and use it to generate a suite of synthetic surveys resembling Gaia DR2 in data structure, magnitude limits, and observational errors. We use three cosmological simulations of Milky-Way-mass galaxies from the Latte suite of the Feedback In Realistic Environments (FIRE) project, which feature self-consistent clustering of star formation in dense molecular clouds and thin stellar/gaseous disks in live cosmological halos with satellite dwarf galaxies and stellar halos. We select three solar viewpoints from each simulation to generate nine synthetic Gaia-like surveys. We sample synthetic stars by assuming each star particle (of mass 7070 $M_{\odot}$) represents a single stellar population. At each viewpoint, we compute dust extinction from the simulated gas metallicity distribution and apply a simple error model to produce a synthetic Gaia-like survey that includes both observational properties and a pointer to the generating star particle. We provide the complete simulation snapshot at $z = 0$ for each simulated galaxy. We describe data access points, the data model, and plans for future upgrades. These synthetic surveys provide a tool for the scientific community to test analysis methods and interpret Gaia data.

preprint2020arXiv

The Effect of Jet-Ejecta Interaction on the Viewing Angle Dependence of Kilonova Light Curves

The merger of two neutron stars produces an outflow of radioactive heavy nuclei. Within a second of merger, the central remnant is expected to also launch a relativistic jet, which shock-heats and disrupts a portion of the radioactive ejecta. Within a few hours, emission from the radioactive material gives rise to an ultraviolet, optical, and infrared transient (a kilonova). We use the endstates of a suite of 2D relativistic hydrodynamic simulations of jet-ejecta interaction as initial conditions for multi-dimensional Monte Carlo radiation transport simulations of the resulting viewing angle-dependent light curves and spectra starting at 1.5h after merger. We find that on this timescale, jet shock heating does not affect the kilonova emission. However, the jet disruption to the density structure of the ejecta does change the light curves. The jet carves a channel into the otherwise spheroidal ejecta, revealing the hot, inner regions. As seen from near ($\lesssim 30 °$) the jet axis, the kilonova is brighter by a factor of a few and bluer. The strength of this effect depends on the jet parameters, since the light curves of more heavily disrupted ejecta are more strongly affected. The light curves and spectra are also more heavily modified in the ultraviolet than in the optical.

preprint2020arXiv

The Effects of Tilt on the Images of Black Hole Accretion Flows

We analyze two 3D general-relativistic magnetohydrodynamic accretion simulations in the context of how they would manifest in Event Horizon Telescope (EHT) observations of supermassive black holes. The two simulations differ only in whether the initial angular momentum of the plasma is aligned with the rapid (a = 0.9) spin of the black hole. Both have low net magnetic flux. Ray tracing is employed to generate resolved images of the synchrotron emission. When using parameters appropriate for Sgr A* and assuming a viewing angle aligned with the black hole spin, we find the most prominent difference is that the central shadow in the image is noticeably eccentric in tilted models, with the ring of emission relatively unchanged. Applying this procedure to M87 with a viewing angle based on the large-scale jet, we find that adding tilt increases the angular size of the ring for fixed black hole mass and distance, while at the same time increasing the number of bright spots in the image. Our findings illustrate observable features that can distinguish tilted from aligned flows. They also show that tilted models can be viable for M87, and that not accounting for tilt can bias inferences of physical parameters. Future modeling of horizon-scale observations should account for potential angular momentum misalignment, which is likely generic at the low accretion rates appropriate for EHT targets.

preprint2020arXiv

The Impact of Type Ia Supernovae in Quiescent Galaxies: I. Formation of the Multiphase Interstellar medium

A cool phase of the interstellar medium has been observed in many giant elliptical galaxies, but its origin remains unclear. We propose that uneven heating from Type Ia supernovae (SNe Ia), together with radiative cooling, can lead to the formation of the cool phase. The basic idea is that since SNe Ia explode randomly, gas parcels which are not directly heated by SN shocks will cool, forming multiphase gas. We run a series of idealized high-resolution numerical simulations, and find that cool gas develops even when the overall SNe heating rate $H$ exceeds the cooling rate $C$ by a factor as large as 1.4. We also find that the time for multiphase gas development depends on the gas temperature. When the medium has a temperature $T = 3\times 10^6$ K, the cool phase forms within one cooling time \tc; however, the cool phase formation is delayed to a few times \tc\ for higher temperatures. The main reason for the delay is turbulent mixing. Cool gas formed this way would naturally have a metallicity lower than that of the hot medium. For constant $H/C$, there is more turbulent mixing for higher temperature gas. We note that this mechanism of producing cool gas cannot be captured in cosmological simulations, which usually fail to resolve individual SN remnants.

preprint2020arXiv

The Impact of Type Ia Supernovae in Quiescent Galaxies: II. Energetics and Turbulence

Type Ia supernovae (SNe Ia) provide unique and important feedback in quiescent galaxies, but their impact has been underappreciated. In this paper, we analyze a series of high-resolution simulations to examine the energetics and turbulence of the medium under SNe Ia. We find that when SN remnants are resolved, their effects differ distinctly from a volumetric heating term, as is commonly assumed in unresolved simulations. First, the net heating is significantly higher than expected, by 30$\pm$10\% per cooling time. This is because a large fraction of the medium is pushed into lower densities which cool inefficiently. Second, the medium is turbulent; the root-mean-squared (RMS) velocity of the gas to 20-50 km s$^{-1}$ on a driving scale of tens of parsec. The velocity field of the medium is dominated by compressional modes, which are larger than the solenoidal components by a factor of 3-7. Third, the hot gas has a very broad density distribution. The ratio between the density fluctuations and the RMS Mach number, parameterized as $b$, is 2-20. This is in contrast to previous simulations of turbulent media, which have found $b\lesssim$ 1. The reason for the difference is mainly caused by the \textit{localized} heating of SNe Ia, which creates a large density contrast. Last, the typical length scale of a density fluctuation grows with time, forming increasingly larger bubbles and filamentary ridges. These underlying density fluctuations need to be included when X-ray observations are interpreted.

preprint2020arXiv

The maximum accretion rate of hot gas in dark matter halos

We revisit the question of 'hot mode' versus 'cold mode' accretion onto galaxies using steady-state cooling flow solutions and idealized 3D hydrodynamic simulations. We demonstrate that for the hot accretion mode to exist, the cooling time is required to be longer than the free-fall time near the radius where the gas is rotationally-supported, R_circ, i.e. the existence of the hot mode depends on physical conditions at the galaxy scale rather than on physical conditions at the halo scale. When allowing for the depletion of the halo baryon fraction relative to the cosmic mean, the longer cooling times imply that a virialized gaseous halo may form in halo masses below the threshold of ~10^12 M_sun derived for baryon-complete halos. We show that for any halo mass there is a maximum accretion rate for which the gas is virialized throughout the halo and can accrete via the hot mode of Mdot_crit ~ 0.7(v_c/100 km/s)^5.4 (R_circ / 10 kpc) (Z/Z_sun)^-0.9 M_sun / yr, where Z and v_c are the metallicity and circular velocity measured at R_circ. For accretion rates >~Mdot_crit the volume-filling gas phase can in principle be `transonic' -- virialized in the outer halo but cool and free-falling near the galaxy. We compare Mdot_crit to the average star formation rate (SFR) in halos at 0<z<10 implied by the stellar-mass halo-mass relation. For a plausible metallicity evolution with redshift, we find that SFR <~ Mdot_crit at most masses and redshifts, suggesting that the SFR of galaxies could be primarily sustained by the hot mode in halo masses well below the classic threshold of ~10^12 M_sun.

preprint2020arXiv

The Structure of Radiatively Inefficient Black Hole Accretion Flows

We run three long-timescale general-relativistic magnetohydrodynamic simulations of radiatively inefficient accretion flows onto non-rotating black holes. Our aim is to achieve steady-state behavior out to large radii and understand the resulting flow structure. A simulation with adiabatic index Gamma = 4/3 and small initial alternating poloidal magnetic field loops is run to a time of 440,000 GM/c^3, reaching inflow equilibrium inside a radius of 370 GM/c^2. Variations with larger alternating field loops and with Gamma = 5/3 are run to 220,000 GM/c^3, attaining equilibrium out to 170 GM/c^2 and 440 GM/c^2. There is no universal self-similar behavior obtained at radii in inflow equilibrium: the Gamma = 5/3 simulation shows a radial density profile with power law index ranging from -1 in the inner regions to -1/2 in the outer regions, while the others have a power-law slope ranging from -1/2 to close to -2. Both simulations with small field loops reach a state with polar inflow of matter, while the more ordered initial field has polar outflows. However, unbound outflows remove only a factor of order unity of the inflowing material over a factor of ~300 in radius. Our results suggest that the dynamics of radiatively inefficient accretion flows are sensitive to how the flow is fed from larger radii, and may differ appreciably in different astrophysical systems. Millimeter images appropriate for Sgr A* are qualitatively (but not quantitatively) similar in all simulations, with a prominent asymmetric image due to Doppler boosting.

preprint2020arXiv

The Surprisingly Small Impact of Magnetic Fields On The Inner Accretion Flow of Sagittarius A* Fueled By Stellar Winds

We study the flow structure in 3D magnetohydrodynamic (MHD) simulations of accretion onto Sagittarius A* via the magnetized winds of the orbiting Wolf-Rayet stars. These simulations cover over 3 orders of magnitude in radius to reach $\approx$ 300 gravitational radii, with only one poorly constrained parameter (the magnetic field in the stellar winds). Even for winds with relatively weak magnetic fields (e.g., plasma $β$ $\sim$ $10^6$), flux freezing/compression in the inflowing gas amplifies the field to $β$ $\sim$ few well before it reaches the event horizon. Overall, the dynamics, accretion rate, and spherically averaged flow profiles (e.g., density, velocity) in our MHD simulations are remarkably similar to analogous hydrodynamic simulations. We attribute this to the broad distribution of angular momentum provided by the stellar winds, which sources accretion even absent much angular momentum transport. We find that the magneto-rotational instability is not important because of i) strong magnetic fields that are amplified by flux freezing/compression, and ii) the rapid inflow/outflow times of the gas and inefficient radiative cooling preclude circularization. The primary effect of magnetic fields is that they drive a polar outflow that is absent in hydrodynamics. The dynamical state of the accretion flow found in our simulations is unlike the rotationally supported tori used as initial conditions in horizon scale simulations, which could have implications for models being used to interpret Event Horizon Telescope and GRAVITY observations of Sgr A*.

preprint2020arXiv

Thermal Instability of Halo Gas Heated by Streaming Cosmic Rays

Heating of virialized gas by streaming cosmic rays (CRs) may be energetically important in galaxy halos, groups and clusters. We present a linear thermal stability analysis of plasmas heated by streaming CRs. We separately treat equilibria with and without background gradients, and with and without gravity. We include both CR streaming and diffusion along the magnetic-field direction. Thermal stability depends strongly on the ratio of CR pressure to gas pressure, which determines whether modes are isobaric or isochoric. Modes with $\mathbf{k \cdot B }\neq 0$ are strongly affected by CR diffusion. When the streaming time is shorter than the CR diffusion time, thermally unstable modes (with $\mathbf{k \cdot B }\neq 0$) are waves propagating at a speed $\propto$ the Alfvén speed. Halo gas in photoionization equilibrium is thermally stable independent of CR pressure, while gas in collisional ionization equilibrium is unstable for physically realistic parameters. In gravitationally stratified plasmas, the oscillation frequency of thermally overstable modes can be higher in the presence of CR streaming than the buoyancy/free-fall frequency. This may modify the critical $t_{\rm cool}/t_{\rm ff}$ at which multiphase gas is present. The criterion for convective instability of a stratified, CR-heated medium can be written in the familiar Schwarzschild form $d s_{\rm eff} / d z < 0$, where $s_{\rm eff}$ is an effective entropy involving the gas and CR pressures. We discuss the implications of our results for the thermal evolution and multiphase structure of galaxy halos, groups and clusters.

preprint2019arXiv

Direct Detection of Black Hole-Driven Turbulence in the Centers of Galaxy Clusters

Supermassive black holes (SMBHs) are thought to provide energy that prevents catastrophic cooling in the centers of massive galaxies and galaxy clusters. However, it remains unclear how this "feedback" process operates. We use high-resolution optical data to study the kinematics of multi-phase filamentary structures by measuring the velocity structure function (VSF) of the filaments over a wide range of scales in the centers of three nearby galaxy clusters: Perseus, Abell 2597 and Virgo. We find that the motions of the filaments are turbulent in all three clusters studied. There is a clear correlation between features of the VSFs and the sizes of bubbles inflated by SMBH driven jets. Our study demonstrates that SMBHs are the main driver of turbulent gas motions in the centers of galaxy clusters and suggests that this turbulence is an important channel for coupling feedback to the environment. Our measured amplitude of turbulence is in good agreement with Hitomi Doppler line broadening measurement and X-ray surface brightness fluctuation analysis, suggesting that the motion of the cold filaments is well-coupled to that of the hot gas. The smallest scales we probe are comparable to the mean free path in the intracluster medium (ICM). Our direct detection of turbulence on these scales provides the clearest evidence to date that isotropic viscosity is suppressed in the weakly-collisional, magnetized intracluster plasma.

preprint2019arXiv

Low-frequency variability in massive stars: Core generation or surface phenomenon?

Bowman et al. (2019) reported low-frequency photometric variability in 164 O- and B-type stars observed with K2 and TESS. They interpret these motions as internal gravity waves, which could be excited stochastically by convection in the cores of these stars. The detection of internal gravity waves in massive stars would help distinguish between massive stars with convective or radiative cores, determine core size, and would provide important constraints on massive star structure and evolution. In this work, we study the observational signature of internal gravity waves generated by core convection. We calculate the \textit{wave transfer function}, which links the internal gravity wave amplitude at the base of the radiative zone to the surface luminosity variation. This transfer function varies by many orders of magnitude for frequencies $\lesssim 1 \, {\rm d}^{-1}$, and has regularly-spaced peaks near $1 \, {\rm d}^{-1}$ due to standing modes. This is inconsistent with the observed spectra which have smooth ``red noise'' profiles, without the predicted regularly-spaced peaks. The wave transfer function is only meaningful if the waves stay predominately linear. We next show that this is the case: low frequency traveling waves do not break unless their luminosity exceeds the radiative luminosity of the star; and, the observed luminosity fluctuations at high frequencies are so small that standing modes would be stable to nonlinear instability. These simple calculations suggest that the observed low-frequency photometric variability in massive stars is not due to internal gravity waves generated in the core of these stars. We finish with a discussion of (sub)surface convection, which produces low-frequency variability in low-mass stars, very similar to that observed in Bowman et al. (2019) in higher mass stars.

preprint2019arXiv

The Fate of AGB Wind in Massive Galaxies and the ICM

Asymptotic Giant Branch (AGB) winds from evolved stars not only provide a non-trivial amount of mass and energy return, but also produce dust grains in massive elliptical galaxies. Due to the fast stellar velocity and the high ambient temperature, the wind is thought to form a comet-like tail, similar to Mira in the Local Bubble. Many massive elliptical galaxies and cluster central galaxies host extended dusty cold filaments. The fate of the cold dusty stellar wind and its relation to cold filaments are not well understood. In this work, we carry out both analytical and numerical studies of the interaction between an AGB wind and the surrounding hot gas. We find that the cooling time of the tail is inversely proportional to the ambient pressure. In the absence of cooling, or in low pressure environments (e.g., the outskirts of elliptical galaxies), AGB winds are quickly mixed into the hot gas, and all the AGB winds have similar appearance and head-to-tail ratio. In high pressure environments, such as the Local Bubble and the central regions of massive elliptical galaxies, some of the gas in the mixing layer between the stellar wind and the surrounding hot gas can cool efficiently and cause the tail to become longer. Our simulated tail of Mira itself has similar length and velocity to that observed, and appears similar to the simulated AGB tail in the central regions of massive galaxies. We speculate that instead of thermal instability, the induced condensation at the mixing layer of AGB winds may be the origin of cold filaments in massive galaxies and galaxy clusters. This naturally explains the existence of dust and PAH in the filaments.

preprint2019arXiv

The Progenitors of Calcium-Strong Transients

A new class of faint, spectroscopically peculiar transients has emerged in the last decade. We term these events "calcium-strong transients" (CaSTs) because of their atypically high calcium-to-oxygen nebular line ratios. Previous studies have struggled to deduce the identity of their progenitors due to a combination of their extremely extended radial distributions with respect to their host galaxies and their relatively high rate of occurrence. In this work, we find that the CaST radial distribution is consistent with the radial distribution of two populations of stars: old (ages > 5 Gyr), low-metallicity (Z/Zsol < 0.3) stars and globular clusters. While no obvious progenitor scenario arises from considering old, metal-poor stars, the alternative production site of globular clusters leads us to narrow down the list of possible candidates to three binary scenarios: mergers of helium and oxygen/neon white dwarfs; tidal disruptions of helium white dwarfs by neutron stars; and stable accretion from low-mass helium-burning stars onto white dwarfs. While rare in the field, these binary systems can be formed dynamically at much higher rates in globular clusters. Subsequent binary hardening both increases their interaction rate and ejects them from their parent globular clusters prior to mass transfer contact. Their production in, and ejection from, globular clusters may explain their radial distribution and the absence of globular clusters at their explosion site. This model predicts a currently undiscovered high rate of CaSTs in nuclear star clusters. Alternatively, an undetermined progenitor scenario involving old, low-metallicity stars may instead hold the key to understanding CaSTs.

preprint2016arXiv

A Diagnostic for Localizing Red Giant Differential Rotation

We present a simple diagnostic that can be used to constrain the location of the differential rotation in red giants with measured mixed mode rotational splittings. Specifically, in red giants with radii $\sim 4R_\odot$, the splittings of p-dominated modes (sound wave dominated) relative to those of g-dominated modes (internal gravity wave dominated) are sensitive to how much of the differential rotation resides in the outer convection zone versus the radiative interior of the red giant. An independently measured surface rotation rate significantly aids breaking degeneracies in interpreting the measured splittings. We apply our results to existing observations of red giants, particularly those of Kepler-56, and find that most of the differential rotation resides in the radiative region rather than in the convection zone. This conclusion is consistent with results in the literature from rotational inversions, but our results are insensitive to some of the uncertainties in the inversion process and can be readily applied to large samples of red giants with even a modest number of measured rotational splittings. We argue that differential rotation in the radiative interior strongly suggests that angular momentum transport in red giants is dominated by local fluid instabilities rather than large-scale magnetic stresses.

preprint2016arXiv

Binary Stars Can Provide the "Missing Photons" Needed for Reionization

Empirical constraints on reionization require galactic ionizing photon escape fractions fesc>20%, but recent high-resolution radiation-hydrodynamic calculations have consistently found much lower values ~1-5%. While these models include strong stellar feedback and additional processes such as runaway stars, they almost exclusively consider stellar evolution models based on single (isolated) stars, despite the fact that most massive stars are in binaries. We re-visit these calculations, combining radiative transfer and high-resolution cosmological simulations from the Feedback in Realistic Environments (FIRE) project. For the first time, we use a stellar evolution model that includes a physically and observationally motivated treatment of binaries (the BPASS model). Binary mass transfer and mergers enhance the population of massive stars at late times (>3 Myr) after star formation, which in turn strongly enhances the late-time ionizing photon production (especially at low metallicities). These photons are produced after feedback from massive stars has carved escape channels in the ISM, and so efficiently leak out of galaxies. As a result, the time-averaged "effective" escape fraction (ratio of escaped ionizing photons to observed 1500 A photons) increases by factors 4-10, sufficient to explain reionization. While important uncertainties remain, we conclude that binary evolution may be critical for understanding the ionization of the Universe.

preprint2016arXiv

Collisionless Isotropization of the Solar-Wind Protons by Compressive Fluctuations and Plasma Instabilities

Compressive fluctuations are a minor yet significant component of astrophysical plasma turbulence. In the solar wind, long-wavelength compressive slow-mode fluctuations lead to changes in $β_{\parallel \mathrm p}\equiv 8πn_{\mathrm p}k_{\mathrm B}T_{\parallel \mathrm p}/B^2$ and in $R_{\mathrm p}\equiv T_{\perp \mathrm p}/T_{\parallel \mathrm p}$, where $T_{\perp \mathrm p}$ and $T_{\parallel \mathrm p}$ are the perpendicular and parallel temperatures of the protons, $B$ is the magnetic field strength, and $n_{\mathrm p}$ is the proton density. If the amplitude of the compressive fluctuations is large enough, $R_{\mathrm p}$ crosses one or more instability thresholds for anisotropy-driven microinstabilities. The enhanced field fluctuations from these microinstabilities scatter the protons so as to reduce the anisotropy of the pressure tensor. We propose that this scattering drives the average value of $R_{\mathrm p}$ away from the marginal stability boundary until the fluctuating value of $R_{\mathrm p}$ stops crossing the boundary. We model this "fluctuating-anisotropy effect" using linear Vlasov--Maxwell theory to describe the large-scale compressive fluctuations. We argue that this effect can explain why, in the nearly collisionless solar wind, the average value of $R_{\mathrm p}$ is close to unity.

preprint2016arXiv

Magnetorotational Turbulence and Dynamo in a Collisionless Plasma

We present results from the first 3D kinetic numerical simulation of magnetorotational turbulence and dynamo, using the local shearing-box model of a collisionless accretion disc. The kinetic magnetorotational instability grows from a subthermal magnetic field having zero net flux over the computational domain to generate self-sustained turbulence and outward angular-momentum transport. Significant Maxwell and Reynolds stresses are accompanied by comparable viscous stresses produced by field-aligned ion pressure anisotropy, which is regulated primarily by the mirror and ion-cyclotron instabilities through particle trapping and pitch-angle scattering. The latter endow the plasma with an effective viscosity that is biased with respect to the magnetic-field direction and spatio-temporally variable. Energy spectra suggest an Alfvén-wave cascade at large scales and a kinetic-Alfvén-wave cascade at small scales, with strong small-scale density fluctuations and weak non-axisymmetric density waves. Ions undergo non-thermal particle acceleration, their distribution accurately described by a kappa distribution. These results have implications for the properties of low-collisionality accretion flows, such as that near the black hole at the Galactic center.

preprint2016arXiv

PIC Simulations of the Effect of Velocity Space Instabilities on Electron Viscosity and Thermal Conduction

In low-collisionality plasmas, velocity-space instabilities are a key mechanism providing an effective collisionality for the plasma. We use particle-in-cell (PIC) simulations to study the interplay between electron and ion-scale velocity-space instabilities and their effect on electron pressure anisotropy, viscous heating, and thermal conduction. The adiabatic invariance of the magnetic moment in low-collisionality plasmas leads to pressure anisotropy, $p_{\perp,j} > p_{||,j}$, if the magnetic field $\vec{B}$ is amplified ($p_{\perp,j}$ and $p_{||,j}$ denote the pressure of species $j$ [electron, ion] perpendicular and parallel to $\vec{B}$). If the resulting anisotropy is large enough, it can in turn trigger small-scale plasma instabilities. Our PIC simulations explore the nonlinear regime of the mirror, ion-cyclotron, and electron whistler instabilities, through continuous amplification of the magnetic field $|\vec{B}|$ by an imposed shear in the plasma. In the regime $1 \lesssim β_j \lesssim 20$ ($β_j \equiv 8πp_j/|\vec{B}|^2$), the saturated electron pressure anisotropy, $Δp_e/p_{||,e}$, is determined mainly by the (electron-lengthscale) whistler marginal stability condition, with a modest factor of $\sim 1.5-2$ decrease due to the trapping of electrons by the mirrors. We explicitly calculate the mean free path of the electrons and ions along the mean magnetic field and provide a simple physical prescription for the mean free path and thermal conductivity in low-collisionality $β_j \gtrsim 1$ plasmas. Our results imply that velocity-space instabilities likely decrease the thermal conductivity of plasma in the outer parts of massive, hot, galaxy clusters. We also discuss the implications of our results for electron heating and thermal conduction in low-collisionality accretion flows onto black holes, including Sgr A* in the Galactic Center.

preprint2016arXiv

Reconciling dwarf galaxies with LCDM cosmology: Simulating a realistic population of satellites around a Milky Way-mass galaxy

Low-mass "dwarf" galaxies represent the most significant challenges to the cold dark matter (CDM) model of cosmological structure formation. Because these faint galaxies are (best) observed within the Local Group (LG) of the Milky Way (MW) and Andromeda (M31), understanding their formation in such an environment is critical. We present first results from the Latte Project: the Milky Way on FIRE (Feedback in Realistic Environments). This simulation models the formation of a MW-mass galaxy to z = 0 within LCDM cosmology, including dark matter, gas, and stars at unprecedented resolution: baryon particle mass of 7070 Msun with gas kernel/softening that adapts down to 1 pc (with a median of 25 - 60 pc at z = 0). Latte was simulated using the GIZMO code with a mesh-free method for accurate hydrodynamics and the FIRE-2 model for star formation and explicit feedback within a multi-phase interstellar medium. For the first time, Latte self-consistently resolves the spatial scales corresponding to half-light radii of dwarf galaxies that form around a MW-mass host down to Mstar > 10^5 Msun. Latte's population of dwarf galaxies agrees with the LG across a broad range of properties: (1) distributions of stellar masses and stellar velocity dispersions (dynamical masses), including their joint relation; (2) the mass-metallicity relation; and (3) a diverse range of star-formation histories, including their mass dependence. Thus, Latte produces a realistic population of dwarf galaxies at Mstar > 10^5 Msun that does not suffer from the "missing satellites" or "too big to fail" problems of small-scale structure formation. We conclude that baryonic physics can reconcile observed dwarf galaxies with standard LCDM cosmology.

preprint2016arXiv

Resonance locking as the source of rapid tidal migration in the Jupiter and Saturn moon systems

The inner moons of Jupiter and Saturn migrate outwards due to tidal energy dissipation within the planets, the details of which remain poorly understood. We demonstrate that resonance locking between moons and internal oscillation modes of the planet can produce rapid tidal migration. Resonance locking arises due to the internal structural evolution of the planet and typically produces an outward migration rate comparable to the age of the solar system. Resonance locking predicts a similar migration timescale but a different effective tidal quality factor $Q$ governing the migration of each moon. It also predicts nearly constant migration timescales a function of semi-major axis, such that effective $Q$ values were larger in the past. Recent measurements of Jupiter and Saturn's moon systems find effective $Q$ values that are smaller than expected (and are different between moons), and which correspond to migration timescales of $\sim$10 Gyr. If confirmed, the measurements are broadly consistent with resonance locking as the dominant source of tidal dissipation in Jupiter and Saturn. Resonance locking also provides solutions to several problems posed by current measurements: it naturally explains the exceptionally small $Q$ governing Rhea's migration, it allows the large heating rate of Enceladus to be achieved in an equilibrium eccentricity configuration, and it resolves evolutionary problems arising from present-day migration/heating rates.

preprint2016arXiv

Stellar and Quasar Feedback in Concert: Effects on AGN Accretion, Obscuration, and Outflows

We study the interaction of feedback from active galactic nuclei (AGN) and a multi-phase interstellar medium (ISM), in simulations including explicit stellar feedback, multi-phase cooling, accretion-disk winds, and Compton heating. We examine radii ~0.1-100 pc around a black hole (BH), where the accretion rate onto the BH is determined and where AGN-powered winds and radiation couple to the ISM. We conclude: (1) The BH accretion rate is determined by exchange of angular momentum between gas and stars in gravitational instabilities. This produces accretion rates ~0.03-1 Msun/yr, sufficient to power luminous AGN. (2) The gas disk in the galactic nucleus undergoes an initial burst of star formation followed by several Myrs where stellar feedback suppresses the star formation rate (SFR). (3) AGN winds injected at small radii with momentum fluxes ~L/c couple efficiently to the ISM and have dramatic effects on ISM properties within ~100 pc. AGN winds suppress the nuclear SFR by factors ~10-30 and BH accretion rate by factors ~3-30. They increase the outflow rate from the nucleus by factors ~10, consistent with observational evidence for galaxy-scale AGN-driven outflows. (4) With AGN feedback, the predicted column density distribution to the BH is consistent with observations. Absent AGN feedback, the BH is isotropically obscured and there are not enough optically-thin sightlines to explain Type-I AGN. A 'torus-like' geometry arises self-consistently as AGN feedback evacuates gas in polar regions.

preprint2016arXiv

Supernova feedback in a local vertically stratified medium: interstellar turbulence and galactic winds

We use local Cartesian simulations with a vertical gravitational potential to study how supernova (SN) feedback in stratified galactic discs drives turbulence and launches galactic winds. Our analysis includes three disc models with gas surface densities ranging from Milky Way-like galaxies to gas-rich ultra-luminous infrared galaxies (ULIRGs), and two different SN driving schemes (random and correlated with local gas density). In order to isolate the physics of SN feedback, we do not include additional feedback processes. We find that, in these local box calculations, SN feedback excites relatively low mass-weighted gas turbulent velocity dispersions ~3-7 km/s and low wind mass loading factors < 1 in all the cases we study. The low turbulent velocities and wind mass loading factors predicted by our local box calculations are significantly below those suggested by observations of gas-rich and rapidly star-forming galaxies; they are also in tension with global simulations of disc galaxies regulated by stellar feedback. Using a combination of numerical tests and analytic arguments, we argue that local Cartesian boxes cannot predict the properties of galactic winds because they do not capture the correct global geometry and gravitational potential of galaxies. The wind mass loading factors are in fact not well-defined in local simulations because they decline significantly with increasing box height. More physically realistic calculations (e.g., including a global galactic potential and disc rotation) will likely be needed to fully understand disc turbulence and galactic outflows, even for the idealized case of feedback by SNe alone.

preprint2016arXiv

Supernova Feedback in an Inhomogeneous Interstellar Medium

Supernova (SN) feedback is one of the key processes shaping the interstellar medium (ISM) of galaxies. SNe contribute to (and in some cases may dominate) driving turbulence in the ISM and accelerating galactic winds. Modern cosmological simulations have sufficient resolution to capture the main structures in the ISM of galaxies, but are typically still not capable of explicitly resolving all of the small-scale stellar feedback processes, including the expansion of supernova remnants (SNRs). We perform a series of controlled three-dimensional hydrodynamic (adaptive mesh refinement) simulations of single SNRs expanding in an inhomogeneous density field with statistics motivated by those of the turbulent ISM. We use these to quantify the momentum and thermal energy injection from SNe as a function of spatial scale and the density, metallicity, and structure of the ambient medium. We develop a series of analytic formulae that we fit to the simulations. These formulae can be used as a basis for a more predictive sub-resolution model for SN feedback for galaxy formation simulations. We then use simulations of multiple, stochastically driven SNe that resolve the key phases of SNRs to test the sub-resolution model, and show that it accurately captures the turbulent kinetic energy and thermal energy in the ISM. By contrast, proposed SN feedback models in the literature based on `delayed cooling' significantly overpredict the late-time thermal energy and momentum in SNRs.

preprint2016arXiv

The Evolution and Fate of Super-Chandrasekhar Mass White Dwarf Merger Remnants

We present stellar evolution calculations of the remnant of the merger of two carbon-oxygen white dwarfs (CO WDs). We focus on cases that have a total mass in excess of the Chandrasekhar mass. After the merger, the remnant manifests as an $L \sim 3 \times 10^4 L_\odot$ source for $\sim 10^4$ yr. A dusty wind may develop, leading these sources to be self-obscured and to appear similar to extreme AGB stars. Roughly $\sim 10$ such objects should exist in the Milky Way and M31 at any time. As found in previous work, off-center carbon fusion is ignited within the merger remnant and propagates inward via a carbon flame, converting the WD to an oxygen-neon (ONe) composition. By following the evolution for longer than previous calculations, we demonstrate that after carbon-burning reaches the center, neutrino-cooled Kelvin-Helmholtz contraction leads to off-center neon ignition in remnants with masses $\ge 1.35 M_\odot$. The resulting neon-oxygen flame converts the core to a silicon WD. Thus, super-Chandrasekhar WD merger remnants do not undergo electron-capture induced collapse as traditionally assumed. Instead, if the remnant mass remains above the Chandrasekhar mass, we expect that it will form a low-mass iron core and collapse to form a neutron star. Remnants that lose sufficient mass will end up as massive, isolated ONe or Si WDs.

preprint2016arXiv

The formation of massive, quiescent galaxies at cosmic noon

The cosmic noon (z~1.5-3) marked a period of vigorous star formation for most galaxies. However, about a third of the more massive galaxies at those times were quiescent in the sense that their observed stellar populations are inconsistent with rapid star formation. The reduced star formation activity is often attributed to gaseous outflows driven by feedback from supermassive black holes, but the impact of black hole feedback on galaxies in the young Universe is not yet definitively established. We analyze the origin of quiescent galaxies with the help of ultra-high resolution, cosmological simulations that include feedback from stars but do not model the uncertain consequences of black hole feedback. We show that dark matter halos with specific accretion rates below ~0.25-0.4 per Gyr preferentially host galaxies with reduced star formation rates and red broad-band colors. The fraction of such halos in large dark matter only simulations matches the observed fraction of massive quiescent galaxies (~10^10-10^11 Msun). This strongly suggests that halo accretion rate is the key parameter determining which massive galaxies at z~1.5-3 become quiescent. Empirical models that connect galaxy and halo evolution, such as halo occupation distribution or abundance matching models, assume a tight link between galaxy properties and the masses of their parent halos. These models will benefit from adding the specific accretion rate of halos as a second model parameter.

preprint2016arXiv

The impact of stellar feedback on hot gas in galaxy haloes: the Sunyaev-Zel'dovich effect and soft X-ray emission

The thermal Sunyaev-Zel'dovich (SZ) effect and soft X-ray emission are routinely observed around massive galaxies and in galaxy groups and clusters. We study these observational diagnostics of galaxy haloes for a suite of cosmological `zoom-in' simulations from the `Feedback In Realistic Environments' project, which spans a large range in halo mass 10^10-10^13 Msun). We explore the effect of stellar feedback on the hot gas observables. The properties of our simulated groups, such as baryon fractions, SZ flux, and X-ray luminosities (L_X), are broadly consistent with existing observations, even though feedback from active galactic nuclei is not included. We make predictions for future observations of lower-mass objects for both SZ and diffuse X-ray measurements, finding that they are not just scaled-down versions of massive galaxies, but more strongly affected by galactic winds driven by star formation. Low-mass haloes (<~10^11 Msun) retain a low fraction of their baryons, which results in a strong suppression of the SZ signal. Our simulations therefore predict a scaling with halo mass that is steeper than self-similar for haloes less massive than 10^13 Msun. For halo masses <~10^12 Msun, L_X is time-variable and correlated primarily with the star formation rate (SFR). For these objects, the diffuse X-ray emission is powered mostly by galactic winds and the gas dominating the X-ray emission is flowing out with radial velocities close to the halo's circular velocity. For halo masses >~10^13 Msun, on the other hand, L_X is much less variable and not correlated with the SFR, because the emission originates from the quasi-hydrostatic, virialized halo gas.

preprint2016arXiv

Thermal Runaway During the Evolution of ONeMg Cores towards Accretion-Induced Collapse

We study the evolution of degenerate electron cores primarily composed of the carbon burning products oxygen, neon, and magnesium (hereafter ONeMg cores) that are undergoing compression. Electron capture reactions on A=20 and A=24 isotopes reduce the electron fraction and heat the core. We develop and use a new capability of the Modules for Experiments in Stellar Astrophysics (MESA) stellar evolution code that provides a highly accurate implementation of these key reactions. These new accurate rates and the ability of MESA to perform extremely small spatial zoning demonstrates a thermal runaway in the core triggered by the temperature and density sensitivity of the Ne-20 electron capture reactions. Both analytics and numerics show that this thermal runaway does not trigger core convection, but rather leads to a centrally concentrated (r < km) thermal runaway that will subsequently launch an oxygen deflagration wave from the center of the star. We use MESA to perform a parameter study that quantifies the influence of the magnesium mass fraction, the central temperature, the compression rate, and uncertainties in the electron capture reaction rates on the ONeMg core evolution. This allows us to establish a lower limit on the central density at which the oxygen deflagration wave initiates of $ρ_c > 8.5 \times 10^9\, \textrm{g cm}^{-3}$. Based on previous work and order-of-magnitude calculations, we expect objects which ignite oxygen at or above these densities to collapse and form a neutron star. Calculations such as these are an important step in producing more realistic progenitor models for studies of the signature of accretion-induced collapse.

preprint2016arXiv

Turbulent Chemical Diffusion in Convectively Bounded Carbon Flames

It has been proposed that mixing induced by convective overshoot can disrupt the inward propagation of carbon deflagrations in super-asymptotic giant branch stars. To test this theory, we study an idealized model of convectively bounded carbon flames with 3D hydrodynamic simulations of the Boussinesq equations using the pseudospectral code Dedalus. Because the flame propagation timescale is much longer than the convection timescale, we approximate the flame as fixed in space, and only consider its effects on the buoyancy of the fluid. By evolving a passive scalar field, we derive a {\it turbulent} chemical diffusivity produced by the convection as a function of height, $D_{\rm t}(z)$. Convection can stall a flame if the chemical mixing timescale, set by the turbulent chemical diffusivity, $D_{\rm t}$, is shorter than the flame propagation timescale, set by the thermal diffusivity, $κ$, i.e., when $D_{\rm t}>κ$. However, we find $D_{\rm t}<κ$ for most of the flame because convective plumes are not dense enough to penetrate into the flame. Extrapolating to realistic stellar conditions, this implies that convective mixing cannot stall a carbon flame and that "hybrid carbon-oxygen-neon" white dwarfs are not a typical product of stellar evolution.

preprint2015arXiv

A Narrow Short-Duration GRB Jet from a Wide Central Engine

We use two-dimensional relativistic hydrodynamic numerical calculations to show that highly collimated relativistic jets can be produced in neutron star merger models of short-duration gamma ray bursts without the need for a highly directed engine or a large net magnetic flux. Even a hydrodynamic engine generating a very wide sustained outflow on small scales can in principle produce a highly collimated relativistic jet, facilitated by a dense surrounding medium which provides a cocoon surrounding the jet core. An oblate geometry to the surrounding gas significantly enhances the collimation process. Previous numerical simulations have shown that the merger of two neutron stars produces an oblate, expanding cloud of dynamical ejecta. We show that this gas can efficiently collimate the central engine power much as the surrounding star does in long-duration GRB models. For typical short-duration GRB central engine parameters, we find jets with opening angles of order 10 degrees in which a large fraction of the total outflow power of the central engine resides in highly relativistic material. These results predict large differences in the opening angles of outflows from binary neutron star mergers versus neutron star-black hole mergers.

preprint2015arXiv

An Extended Magnetohydrodynamics Model for Relativistic Weakly Collisional Plasmas

Black holes that accrete far below the Eddington limit are believed to accrete through a geometrically thick, optically thin, rotationally supported plasma that we will refer to as a radiatively inefficient accretion flow (RIAF). RIAFs are typically collisionless in the sense that the Coulomb mean free path is large compared to $GM/c^2$, and relativistically hot near the event horizon. In this paper we develop a phenomenological model for the plasma in RIAFs, motivated by the application to sources such as Sgr A* and M87. The model is derived using Israel-Stewart theory, which considers deviations up to second order from thermal equilibrium, but modified for a magnetized plasma. This leads to thermal conduction along magnetic field lines and a difference in pressure, parallel and perpendicular to the field lines (which is equivalent to anisotrotropic viscosity). In the non-relativistic limit, our model reduces to the widely used Braginskii theory of magnetized, weakly collisional plasmas. We compare our model to the existing literature on dissipative relativistic fluids, describe the linear theory of the plasma, and elucidate the physical meaning of the free parameters in the model. We also describe limits of the model when the conduction is saturated and when the viscosity implies a large pressure anisotropy. In future work, the formalism developed in this paper will be used in numerical models of RIAFs to assess the importance of non-ideal processes for the dynamics and radiative properties of slowly accreting black holes.

preprint2015arXiv

An Origin for Multi-Phase Gas in Galactic Winds and Halos

The physical origin of high velocity cool gas seen in galactic winds remains unknown. Following Wang (1995), we argue that radiative cooling in initially hot thermally-driven outflows can produce fast neutral atomic and photoionized cool gas. The inevitability of adiabatic cooling from the flow's initial 10^7-10^8K temperature and the shape of the cooling function for T<10^7K imply that outflows with hot gas mass-loss rate relative to star formation rate of beta=Mdot_hot/Mdot_star > 0.5 cool radiatively on scales ranging from the size of the energy injection region to tens of kpc. We highlight the beta and star formation rate surface density dependence of the column density, emission measure, radiative efficiency, and velocity. At r_cool, the gas produces X-ray and then UV/optical line emission with a total power bounded by 10^{-2} L_star if the flow is powered by steady-state star formation with luminosity L_star. The wind is thermally unstable at r_cool, potentially leading to a multi-phase medium. Cooled winds decelerate significantly in the extended gravitational potential of galaxies. The cool gas precipitated from hot outflows may explain its prevalence in galactic halos. We forward a picture of winds whereby cool clouds are initially accelerated by the ram pressure of the hot flow, but are rapidly shredded by hydrodynamical instabilities, thereby increasing beta, seeding radiative and thermal instability, and cool gas rebirth. If the cooled wind shocks as it sweeps up the circumgalactic medium, its cooling time is short, thus depositing cool gas far out into the halo. Finally, conduction can dominate energy transport in low-beta hot winds, leading to flatter temperature profiles than otherwise expected, potentially consistent with X-ray observations of some starbursts.

preprint2015arXiv

Dynamics of Dusty Radiation Pressure Driven Shells and Clouds: Fast Outflows from Galaxies, Star Clusters, Massive Stars, and AGN

It is typically assumed that radiation pressure driven winds are accelerated to an asymptotic velocity of V ~ v_esc, where v_esc is the escape velocity from the central source. We note that this is not the case for dusty shells and clouds. Instead, if the shell or cloud is initially optically-thick to the UV emission from the source of luminosity L, then there is a significant boost in V that reflects the integral of the momentum absorbed as it is accelerated. For shells reaching a generalized Eddington limit, we show that V ~ (4R_UV L/M_sh c)^1/2, in both point-mass and isothermal-sphere potentials, where R_UV is the radius where the shell becomes optically-thin to UV photons, and M_sh is the mass of the shell. The asymptotic velocity significantly exceeds v_esc for typical parameters, and can explain the ~1000-2000km/s outflows observed from rapidly star-forming galaxies and active galactic nuclei if the surrounding halo has low gas density. Similarly fast outflows from massive stars can be accelerated on few - 10^3 yr timescales. These results carry over to clouds that subtend only a small fraction of the solid angle from the source of radiation and that expand as a consequence of their internal sound speed. We further consider the dynamics of shells that sweep up a dense circumstellar or circumgalactic medium. We calculate the "momentum ratio" Mdot v/(L/c) in the shell limit and show that it can only significantly exceed ~2 if the effective optical depth of the shell to re-radiated FIR photons is much larger than unity. We discuss simple prescriptions for the properties of galactic outflows for use in large-scale cosmological simulations. We also briefly discuss applications to the dusty ejection episodes of massive stars, the disruption of giant molecular clouds, and AGN.

preprint2015arXiv

Electron Thermodynamics in GRMHD Simulations of Low-Luminosity Black Hole Accretion

Simple assumptions made regarding electron thermodynamics often limit the extent to which general relativistic magnetohydrodynamic (GRMHD) simulations can be applied to observations of low-luminosity accreting black holes. We present, implement, and test a model that self-consistently evolves an electron entropy equation and takes into account the effects of spatially varying electron heating and relativistic anisotropic thermal conduction along magnetic field lines. We neglect the back-reaction of electron pressure on the dynamics of the accretion flow. Our model is appropriate for systems accreting at $\ll 10^{-5}$ of the Eddington rate, so radiative cooling by electrons can be neglected. It can be extended to higher accretion rates in the future by including electron cooling and proton-electron Coulomb collisions. We present a suite of tests showing that our method recovers the correct solution for electron heating under a range of circumstances, including strong shocks and driven turbulence. Our initial applications to axisymmetric simulations of accreting black holes show that (1)~physically-motivated electron heating rates yield electron temperature distributions significantly different from the constant electron to proton temperature ratios assumed in previous work, with higher electron temperatures concentrated in the coronal region between the disc and the jet; (2)~electron thermal conduction significantly modifies the electron temperature in the inner regions of black hole accretion flows if the effective electron mean free path is larger than the local scale-height of the disc (at least for the initial conditions and magnetic field configurations we study). The methods developed in this work are important for producing more realistic predictions for the emission from accreting black holes such as Sagittarius A* and M87; these applications will be explored in future work.

preprint2015arXiv

Evolution of Accretion Discs around a Kerr Black Hole using Extended Magnetohydrodynamics

Black holes accreting well below the Eddington rate are believed to have geometrically thick, optically thin, rotationally supported accretion discs in which the Coulomb mean free path is large compared to $GM/c^2$. In such an environment, the disc evolution may differ significantly from ideal magnetohydrodynamic predictions. We present non-ideal global axisymmetric simulations of geometrically thick discs around a rotating black hole. The simulations are carried out using a new code ${\rm\it grim}$, which evolves a covariant extended magnetohydrodynamics model derived by treating non-ideal effects as a perturbation of ideal magnetohydrodynamics. Non-ideal effects are modeled through heat conduction along magnetic field lines, and a difference between the pressure parallel and perpendicular to the field lines. The model relies on an effective collisionality in the disc from wave-particle scattering and velocity-space (mirror and firehose) instabilities. We find that the pressure anisotropy grows to match the magnetic pressure, at which point it saturates due to the mirror instability. The pressure anisotropy produces outward angular momentum transport with a magnitude comparable to that of MHD turbulence in the disc, and a significant increase in the temperature in the wall of the jet. We also find that, at least in our axisymmetric simulations, conduction has a small effect on the disc evolution because (1) the heat flux is constrained to be parallel to the field and the field is close to perpendicular to temperature gradients, and (2) the heat flux is choked by an increase in effective collisionality associated with the mirror instability.

preprint2015arXiv

Fast Luminous Blue Transients from Newborn Black Holes

Newborn black holes in collapsing massive stars can be accompanied by a fallback disk. The accretion rate is typically super-Eddington and strong disk outflows are expected. Such outflows could be directly observed in some failed explosions of compact (blue supergiants or Wolf-Rayet stars) progenitors, and may be more common than long-duration gamma-ray bursts. Using an analytical model, we show that the fallback disk outflows produce blue UV-optical transients with a peak bolometric luminosity of ~10^(42-43) erg s^-1 (peak R-band absolute AB magnitudes of -16 to -18) and an emission duration of ~ a few to ~ 10 days. The spectra are likely dominated intermediate mass elements, but will lack much radioactive nuclei and iron-group elements. The above properties are broadly consistent with some of the rapid blue transients detected by Pan-STARRS and PTF. This scenario can be distinguished from alternative models using radio observations within a few years after the optical peak.

preprint2015arXiv

Forged in FIRE: cusps, cores, and baryons in low-mass dwarf galaxies

We present ultra-high resolution cosmological hydrodynamic simulations of $M_*\simeq10^{4-6.3}M_{\odot}$ dwarf galaxies that form within $M_{v}=10^{9.5-10}M_{\odot}$ dark matter halos. Our simulations rely on the FIRE implementation of star formation feedback and were run with high enough force and mass resolution to directly resolve stellar and dark matter structure on the ~200 pc scales of interest for classical and ultra-faint dwarfs in the Local Group. The resultant galaxies sit on the $M_*$ vs. $M_{v}$ relation required to match the Local Group stellar mass function. They have bursty star formation histories and also form with half-light radii and metallicities that broadly match those observed for local dwarfs at the same stellar mass. We demonstrate that it is possible to create a large (~1 kpc) dark matter core in a cosmological simulation of an $M_*\simeq10^{6.5}M_{\odot}$ dwarf galaxy that resides within an $M_{v}=10^{10}M_{\odot}$ halo -- precisely the scale of interest for resolving the Too Big to Fail problem. However, these large cores are not ubiquitous and appear to correlate closely with the star formation histories of the dwarfs: dark matter cores are largest in systems that form their stars late ($z\lesssim2$), after the early epoch of cusp building mergers has ended. Our $M_*\simeq10^4M_{\odot}$ dwarf retains a cuspy dark matter halo density profile that matches almost identically that of a dark-matter only run of the same system. Despite forming in a field environment, this very low mass dwarf has observable properties that match closely to those of ultra-faint satellite galaxies of the Milky Way, including a uniformly old stellar population (>10 Gyr). Though ancient, most of the stars in our ultra-faint form after reionization; the UV field acts mainly to suppress fresh gas accretion, not to boil away gas that is already present in the proto-dwarf.

preprint2015arXiv

Gusty, gaseous flows of FIRE: galactic winds in cosmological simulations with explicit stellar feedback

We present an analysis of the galaxy-scale gaseous outflows from the FIRE (Feedback in Realistic Environments) simulations. This suite of hydrodynamic cosmological zoom simulations resolves formation of star-forming giant molecular clouds to $z=0$, and features an explicit stellar feedback model on small scales. Our simulations reveal that high redshift galaxies undergo bursts of star formation followed by powerful gusts of galactic outflows that eject much of the ISM and temporarily suppress star formation. At low redshift, however, sufficiently massive galaxies corresponding to L*-progenitors develop stable disks and switch into a continuous and quiescent mode of star formation that does not drive outflows far into the halo. Mass-loading factors for winds in L*-progenitors are $η\approx 10$ at high redshift, but decrease to $η\ll 1$ at low redshift. Although lower values of $η$ are expected as halos grow in mass over time, we show that the strong suppression of outflows with decreasing redshift cannot be explained by mass evolution alone. Circumgalactic outflow velocities are variable and broadly distributed, but typically range between one and three times the circular velocity of the halo. Much of the ejected material builds a reservoir of enriched gas within the circumgalactic medium, some of which could be later recycled to fuel further star formation. However, a fraction of the gas that leaves the virial radius through galactic winds is never regained, causing most halos with mass $M_h \le 10^{12} M_{\odot}$ to be deficient in baryons compared to the cosmic mean by $z=0$.

preprint2015arXiv

Linear Vlasov Theory in the Shearing Sheet Approximation with Application to the Magneto-Rotational Instability

We derive the conductivity tensor for axisymmetric perturbations of a hot, collisionless, and charge-neutral plasma in the shearing sheet approximation. Our results generalize the well-known linear Vlasov theory for uniform plasmas to differentially rotating plasmas and can be used for wide range of kinetic stability calculations. We apply these results to the linear theory of the magneto-rotational instability (MRI) in collisionless plasmas. We show analytically and numerically how the general kinetic theory results derived here reduce in appropriate limits to previous results in the literature, including the low frequency guiding center (or "kinetic MHD") approximation, Hall MHD, and the gyro-viscous approximation. We revisit the cold plasma model of the MRI and show that, contrary to previous results, an initially unmagnetized collisionless plasma is linearly stable to axisymmetric perturbations in the cold plasma approximation. In addition to their application to astrophysical plasmas, our results provide a useful framework for assessing the linear stability of differentially rotating plasmas in laboratory experiments.

preprint2015arXiv

Local Radiation Hydrodynamic Simulations of Massive Star Envelopes at the Iron Opacity Peak

We perform three dimensional radiation hydrodynamic simulations of the structure and dynamics of radiation dominated envelopes of massive stars at the location of the iron opacity peak. One dimensional hydrostatic calculations predict an unstable density inversion at this location, whereas our simulations reveal a complex interplay of convective and radiative transport whose behavior depends on the ratio of the photon diffusion time to the dynamical time. The latter is set by the ratio of the optical depth per pressure scale height, $τ_0$, to $τ_c=c/c_g$, where $c_g \approx$ 50 km/s is the isothermal sound speed in the gas alone. When $τ_0 \gg τ_c$, convection reduces the radiation acceleration and removes the density inversion. The turbulent energy transport in the simulations agrees with mixing length theory and provides its first numerical calibration in the radiation dominated regime. When $τ_0 \ll τ_c$, convection becomes inefficient and the turbulent energy transport is negligible. The turbulent velocities exceed $c_g$, driving shocks and large density fluctuations that allow photons to preferentially diffuse out through low-density regions. However, the effective radiation acceleration is still larger than the gravitational acceleration so that the time average density profile contains a modest density inversion. In addition, the simulated envelope undergoes large-scale oscillations with periods of a few hours. The turbulent velocity field may affect the broadening of spectral lines and therefore stellar rotation measurements in massive stars, while the time variable outer atmosphere could lead to variations in their mass loss and stellar radius.

preprint2015arXiv

Numerical Simulations of Internal Wave Generation by Convection in Water

Water's density maximum at 4C makes it well suited to study internal gravity wave excitation by convection: an increasing temperature profile is unstable to convection below 4C, but stably stratified above 4C. We present numerical simulations of a water-like fluid near its density maximum in a two dimensional domain. We successfully model the damping of waves in the simulations using linear theory, provided we do not take the weak damping limit typically used in the literature. In order to isolate the physical mechanism exciting internal waves, we use the novel spectral code Dedalus to run several simplified model simulations of our more detailed simulation. We use data from the full simulation as source terms in two simplified models of internal wave excitation by convection: bulk excitation by convective Reynolds stresses, and interface forcing via the mechanical oscillator effect. We find excellent agreement between the waves generated in the full simulation and the simplified simulation implementing the bulk excitation mechanism. The interface forcing simulations over excite high frequency waves because they assume the excitation is by the "impulsive" penetration of plumes, which spreads energy to high frequencies. However, we find the real excitation is instead by the "sweeping" motion of plumes parallel to the interface. Our results imply that the bulk excitation mechanism is a very accurate heuristic for internal wave generation by convection.

preprint2015arXiv

PIC Simulations of Continuously Driven Mirror and Ion Cyclotron Instabilities in High Beta Astrophysical and Heliospheric Plasmas

We use particle-in-cell (PIC) simulations to study the nonlinear evolution of ion velocity space instabilities in an idealized problem in which a background velocity shear continuously amplifies the magnetic field. We simulate the astrophysically relevant regime where the shear timescale is long compared to the ion cyclotron period, and the plasma beta is ~ 1-100. The background field amplification in our calculation is meant to mimic processes such as turbulent fluctuations or MHD-scale instabilities. The field amplification continuously drives a pressure anisotropy with the perpendicular pressure larger than the parallel pressure, and the plasma becomes unstable to the mirror and ion cyclotron instabilities. In all cases, the nonlinear state is dominated by the mirror instability, not the ion cyclotron instability, and the plasma pressure anisotropy saturates near the threshold for the linear mirror instability. The magnetic field fluctuations initially undergo exponential growth but saturate in a secular phase in which the fluctuations grow on the same timescale as the background magnetic field (with δB ~ 0.3 <B> in the secular phase). At early times, the ion magnetic moment is well-conserved but once the fluctuation amplitudes exceed δB ~ 0.1 <B>, the magnetic moment is no longer conserved but instead changes on a timescale comparable to that of the mean magnetic field. We discuss the implications of our results for low-collisionality astrophysical plasmas, including the near-Earth solar wind and low-luminosity accretion disks around black holes.

preprint2015arXiv

Super-Eddington Stellar Winds Driven by Near-Surface Energy Deposition

We develop analytic and numerical models of the properties of super-Eddington stellar winds, motivated by phases in stellar evolution when super-Eddington energy deposition (via, e.g., unstable fusion, wave heating, or a binary companion) heats a region near the stellar surface. This appears to occur in luminous blue variables (LBVs), Type IIn supernovae progenitors, classical novae, and X-ray bursts. We show that when the wind kinetic power exceeds Eddington, the photons are trapped and behave like a fluid. Convection does not play a significant role in the wind energy transport. The wind properties depend on the ratio of a characteristic speed in the problem vc ~ (Edot G)^{1/5} (where Edot is the heating rate) to the stellar escape speed near the heating region vesc(r_h). For vc > vesc(r_h) the wind kinetic power at large radii Edot_w ~ Edot. For vc < vesc(r_h), most of the energy is used to unbind the wind material and thus Edot_w < Edot. Multidimensional hydrodynamic simulations without radiation diffusion using FLASH and one-dimensional hydrodynamic simulations with radiation diffusion using MESA are in good agreement with the analytic predictions. The photon luminosity from the wind is itself super-Eddington but in many cases the photon luminosity is likely dominated by `internal shocks' in the wind. We discuss the application of our models to eruptive mass loss from massive stars and argue that the wind models described here can account for the broad properties of LBV outflows and the enhanced mass loss in the years prior to Type IIn core-collapse supernovae.

preprint2015arXiv

The creation and persistence of a misaligned gas disc in a simulated early-type galaxy

Massive early-type galaxies commonly have gas discs which are kinematically misaligned with the stellar component. These discs feel a torque from the stars and the angular momentum vectors are expected to align quickly. We present results on the evolution of a misaligned gas disc in a cosmological simulation of a massive early-type galaxy from the Feedback In Realistic Environments project. This galaxy experiences a merger which, together with a strong galactic wind, removes most of the original gas disc. The galaxy subsequently reforms a gas disc through accretion of cold gas, but it is initially 120 degrees misaligned with the stellar rotation axis. This misalignment persists for about 2 Gyr before the gas-star misalignment angle drops below 20 degrees. The time it takes for the gaseous and stellar components to align is much longer than previously thought, because the gas disc is accreting a significant amount of mass for about 1.5 Gyr after the merger, during which the angular momentum change induced by accreted gas dominates over that induced by stellar torques. Once the gas accretion rate has decreased sufficiently, the gas disc decouples from the surrounding halo gas and realigns with the stellar component in about 6 dynamical times. During the late evolution of the misaligned gas disc, the centre aligns faster than the outskirts, resulting in a warped disc. We discuss the observational consequences of the long survival of our misaligned gas disc and how our results can be used to calibrate merger rate estimates from observed gas misalignments.

preprint2015arXiv

The Difficulty of Getting High Escape Fractions of Ionizing Photons from High-redshift Galaxies: a View from the FIRE Cosmological Simulations

We present a series of high-resolution (20-2000 Msun, 0.1-4 pc) cosmological zoom-in simulations at z~6 from the Feedback In Realistic Environment (FIRE) project. These simulations cover halo masses 10^9-10^11 Msun and rest-frame ultraviolet magnitude Muv = -9 to -19. These simulations include explicit models of the multi-phase ISM, star formation, and stellar feedback, which produce reasonable galaxy properties at z = 0-6. We post-process the snapshots with a radiative transfer code to evaluate the escape fraction (fesc) of hydrogen ionizing photons. We find that the instantaneous fesc has large time variability (0.01%-20%), while the time-averaged fesc over long time-scales generally remains ~5%, considerably lower than the estimate in many reionization models. We find no strong dependence of fesc on galaxy mass or redshift. In our simulations, the intrinsic ionizing photon budgets are dominated by stellar populations younger than 3 Myr, which tend to be buried in dense birth clouds. The escaping photons mostly come from populations between 3-10 Myr, whose birth clouds have been largely cleared by stellar feedback. However, these populations only contribute a small fraction of intrinsic ionizing photon budgets according to standard stellar population models. We show that fesc can be boosted to high values, if stellar populations older than 3 Myr produce more ionizing photons than standard stellar population models (as motivated by, e.g., models including binaries). By contrast, runaway stars with velocities suggested by observations can enhance fesc by only a small fraction. We show that "sub-grid" star formation models, which do not explicitly resolve star formation in dense clouds with n >> 1 cm^-3, will dramatically over-predict fesc.

preprint2015arXiv

The Diversity of Transients from Magnetar Birth

Strongly-magnetized, rapidly-rotating neutron stars are contenders for the central engines of both long-duration gamma-ray bursts (LGRBs) and hydrogen-poor super-luminous supernovae (SLSNe-I). Models for typical (~minute long) LGRBs invoke magnetars with high dipole magnetic fields (Bd > 1e15 G) and short spin-down times, while models for SLSNe-I invoke neutron stars with weaker fields and longer spin-down times of weeks. Here we identify a transition region in the space of Bd and birth period for which a magnetar can power both a long GRB and a luminous SN. In particular, we show that a 2 ms period magnetar with a spin-down time of ~1e4 s can explain the observations of both the ultra-long GRB 111209 and its associated luminous SN2011kl. For magnetars with longer spin down times, we predict even longer duration (~1e6 s) GRBs and brighter supernovae, a correlation that extends to Swift J2058+05 (commonly interpreted as a tidal disruption event). We further show that previous estimates of the maximum rotational energy of a proto-magnetar were too conservative and energies up to Emax ~1-2e53 erg are possible. The magnetar model can therefore comfortably accommodate the extreme energy requirements recently posed by the most luminous supernova ASASSN-15lh. The high ionization flux from a pulsar wind nebula powering ASASSN-15lh may lead to an "ionization break-out" X-ray burst over the coming months, which would be accompanied by an abrupt change in the optical spectrum. We conclude by briefly contrasting millisecond magnetar and black hole models for SLSNe and ultra-long GRBs.

preprint2015arXiv

The interplay of disk wind and dynamical ejecta in the aftermath of neutron star - black hole mergers

We explore the evolution of the different ejecta components generated during the merger of a neutron star (NS) and a black hole (BH). Our focus is the interplay between material ejected dynamically during the merger, and the wind launched on a viscous timescale by the remnant accretion disk. These components are expected to contribute to an electromagnetic transient and to produce r-process elements, each with a different signature when considered separately. Here we introduce a two-step approach to investigate their combined evolution, using two- and three-dimensional hydrodynamic simulations. Starting from the output of a merger simulation, we identify each component in the initial condition based on its phase space distribution, and evolve the accretion disk in axisymmetry. The wind blown from this disk is injected into a three-dimensional computational domain where the dynamical ejecta is evolved. We find that the wind can suppress fallback accretion on timescales longer than ~100 ms. Due to self-similar viscous evolution, the disk accretion at late times nevertheless approaches a power-law time dependence $\propto t^{-2.2}$. This can power some late-time GRB engine activity, although the available energy is significantly less than in traditional fallback models. Inclusion of radioactive heating due to the r-process does not significantly affect the fallback accretion rate or the disk wind. We do not find any significant modification to the wind properties at large radius due to interaction with the dynamical ejecta. This is a consequence of the different expansion velocities of the two components.

preprint2015arXiv

The Origin and Evolution of the Galaxy Mass-Metallicity Relation

We use high-resolution cosmological zoom-in simulations from the Feedback in Realistic Environment (FIRE) project to study the galaxy mass-metallicity relations (MZR) from z=0-6. These simulations include explicit models of the multi-phase ISM, star formation, and stellar feedback. The simulations cover halo masses Mhalo=10^9-10^13 Msun and stellar mass Mstar=10^4-10^11 Msun at z=0 and have been shown to produce many observed galaxy properties from z=0-6. For the first time, our simulations agree reasonably well with the observed mass-metallicity relations at z=0-3 for a broad range of galaxy masses. We predict the evolution of the MZR from z=0-6 as log(Zgas/Zsun)=12+log(O/H)-9.0=0.35[log(Mstar/Msun)-10]+0.93 exp(-0.43 z)-1.05 and log(Zstar/Zsun)=[Fe/H]-0.2=0.40[log(Mstar/Msun)-10]+0.67 exp(-0.50 z)-1.04, for gas-phase and stellar metallicity, respectively. Our simulations suggest that the evolution of MZR is associated with the evolution of stellar/gas mass fractions at different redshifts, indicating the existence of a universal metallicity relation between stellar mass, gas mass, and metallicities. In our simulations, galaxies above Mstar=10^6 Msun are able to retain a large fraction of their metals inside the halo, because metal-rich winds fail to escape completely and are recycled into the galaxy. This resolves a long-standing discrepancy between "sub-grid" wind models (and semi-analytic models) and observations, where common sub-grid models cannot simultaneously reproduce the MZR and the stellar mass functions.

preprint2015arXiv

The Properties of G-modes in Layered Semi-Convection

We study low frequency waves that propagate in a region of layered semi-convection. Layered semi-convection is predicted to be present in stellar and planetary interiors and can significantly modify the rate of thermal and compositional mixing. We derive a series of analytical dispersion relations for plane-parallel layered semi-convection in the Boussinesq approximation using a matrix transfer formalism. We find that like a continuously stratified medium, a semi-convective staircase -- in which small convective regions are separated by sharp density jumps -- supports internal gravity waves (g-modes). When the wavelength is much longer than the distance between semi-convective steps, these behave nearly like g-modes in a continuously stratified medium. However, the g-mode period spacing in a semi-convective region is systematically {\em smaller} than in a continuously stratified medium, and it decreases with decreasing mode frequency. When the g-mode wavelength becomes comparable to the distance between semi-convective steps, the g-mode frequencies deviate significantly from those of a continuously stratified medium (the frequencies are higher). G-modes with vertical wavelengths smaller than the distance between semi-convective steps are evanescent and do not propagate in the staircase. Thus, there is a lower cutoff frequency for a given horizontal wavenumber. We generalize our results to gravito-inertial waves relevant for rapidly rotating stars and planets. Finally, we assess the prospects for detecting layered semi-convection using astero/planetary seismology.

preprint2015arXiv

The spin rate of pre-collapse stellar cores: wave-driven angular momentum transport in massive stars

The core rotation rates of massive stars have a substantial impact on the nature of core-collapse supernovae and their compact remnants. We demonstrate that internal gravity waves (IGW), excited via envelope convection during a red supergiant phase or during vigorous late time burning phases, can have a significant impact on the rotation rate of the pre-SN core. In typical ($10 \, M_\odot \lesssim M \lesssim 20 \, M_\odot$) supernova progenitors, IGW may substantially spin down the core, leading to iron core rotation periods $P_{\rm min,Fe} \gtrsim 30 \, {\rm s}$. Angular momentum (AM) conservation during the supernova would entail minimum NS rotation periods of $P_{\rm min,NS} \gtrsim 3 \, {\rm ms}$. In most cases, the combined effects of magnetic torques and IGW AM transport likely lead to substantially longer rotation periods. However, the stochastic influx of AM delivered by IGW during shell burning phases inevitably spin up a slowly rotating stellar core, leading to a maximum possible core rotation period. We estimate maximum iron core rotation periods of $P_{\rm max,Fe} \lesssim 5 \times 10^3 \, {\rm s}$ in typical core-collapse supernova progenitors, and a corresponding spin period of $P_{\rm max, NS} \lesssim 500 \, {\rm ms}$ for newborn neutron stars. This is comparable to the typical birth spin periods of most radio pulsars. Stochastic spin-up via IGW during shell O/Si burning may thus determine the initial rotation rate of most neutron stars. For a given progenitor, this theory predicts a Maxwellian distribution in pre-collapse core rotation frequency that is uncorrelated with the spin of the overlying envelope.

preprint2014arXiv

Acceleration of Relativistic Electrons by MHD Turbulence: Implications for Non-thermal Emission from Black Hole Accretion Disks

We use analytic estimates and numerical simulations of test particles interacting with magnetohydrodynamic (MHD) turbulence to show that subsonic MHD turbulence produces efficient second-order Fermi acceleration of relativistic particles. This acceleration is not well-described by standard quasi-linear theory but is a consequence of resonance broadening of wave-particle interactions in MHD turbulence. We provide momentum diffusion coefficients that can be used for astrophysical and heliospheric applications and discuss the implications of our results for accretion flows onto black holes. In particular, we show that particle acceleration by subsonic turbulence in radiatively inefficient accretion flows can produce a non-thermal tail in the electron distribution function that is likely important for modeling and interpreting the emission from low luminosity systems such as Sgr A* and M87.

preprint2014arXiv

Dynamical resonance locking in tidally interacting binary systems

We examine the dynamics of resonance locking in detached, tidally interacting binary systems. In a resonance lock, a given stellar or planetary mode is trapped in a highly resonant state for an extended period of time, during which the spin and orbital frequencies vary in concert to maintain the resonance. This phenomenon is qualitatively similar to resonance capture in planetary dynamics. We show that resonance locks can accelerate the course of tidal evolution in eccentric systems and also efficiently couple spin and orbital evolution in circular binaries. Previous analyses of resonance locking have not treated the mode amplitude as a fully dynamical variable, but rather assumed the adiabatic (i.e. Lorentzian) approximation valid only in the limit of relatively strong mode damping. We relax this approximation, analytically derive conditions under which the fixed point associated with resonance locking is stable, and further check these analytic results using numerical integrations of the coupled mode, spin, and orbital evolution equations. These show that resonance locking can sometimes take the form of complex limit cycles or even chaotic trajectories. We provide simple analytic formulae that define the binary and mode parameter regimes in which resonance locks of some kind occur (stable, limit cycle, or chaotic). We briefly discuss the astrophysical implications of our results for white dwarf and neutron star binaries as well as eccentric stellar binaries.

preprint2014arXiv

Galactic r-process enrichment by neutron star mergers in cosmological simulations of a Milky Way-mass galaxy

We quantify the stellar abundances of neutron-rich r-process nuclei in cosmological zoom-in simulations of a Milky Way-mass galaxy from the Feedback In Realistic Environments project. The galaxy is enriched with r-process elements by binary neutron star (NS) mergers and with iron and other metals by supernovae. These calculations include key hydrodynamic mixing processes not present in standard semi-analytic chemical evolution models, such as galactic winds and hydrodynamic flows associated with structure formation. We explore a range of models for the rate and delay time of NS mergers, intended to roughly bracket the wide range of models consistent with current observational constraints. We show that NS mergers can produce [r-process/Fe] abundance ratios and scatter that appear reasonably consistent with observational constraints. At low metallicity, [Fe/H]<-2, we predict there is a wide range of stellar r-process abundance ratios, with both supersolar and subsolar abundances. Low-metallicity stars or stars that are outliers in their r-process abundance ratios are, on average, formed at high redshift and located at large galactocentric radius. Because NS mergers are rare, our results are not fully converged with respect to resolution, particularly at low metallicity. However, the uncertain rate and delay time distribution of NS mergers introduces an uncertainty in the r-process abundances comparable to that due to finite numerical resolution. Overall, our results are consistent with NS mergers being the source of most of the r-process nuclei in the Universe.

preprint2014arXiv

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_sun, and stellar masses ~10^4-10^11. Our simulations include fully explicit treatment of both the multi-phase ISM (molecular through hot) and stellar feedback. The stellar feedback inputs (energy, momentum, mass, and metal fluxes) are taken directly from stellar population models. These sources of stellar feedback, with zero adjusted parameters, reproduce the observed relation between stellar and halo mass up to M_halo~10^12 M_sun (including dwarfs, satellites, MW-mass disks, and small groups). By extension, this leads to reasonable agreement with the stellar mass function for M_star<10^11 M_sun. We predict weak redshift evolution in the M_star-M_halo relation, consistent with current constraints to z>6. We find that the M_star-M_halo relation is insensitive to numerical details, but is sensitive to the feedback physics. Simulations with only supernova feedback fail to reproduce the observed stellar masses, particularly in dwarf and high-redshift galaxies: radiative feedback (photo-heating and radiation pressure) is necessary to disrupt GMCs and enable efficient coupling of later supernovae to the gas. Star formation rates 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 in the simulation, and is instead determined by self-regulation via stellar feedback. Feedback reduces star formation rates considerably and produces a reservoir of gas that leads to rising late-time star formation histories significantly different from the halo accretion history. Feedback also produces large short-timescale variability in galactic SFRs, especially in dwarfs. Many of these properties are not captured by common 'sub-grid' galactic wind models.

preprint2014arXiv

Hot Galactic Winds Constrained by the X-Ray Luminosities of Galaxies

Galactic superwinds may be driven by very hot outflows generated by overlapping supernovae within the host galaxy. We use the Chevalier & Clegg (CC85) wind model and the observed correlation between X-ray luminosities of galaxies and their SFRs to constrain the mass loss rates (\dot{M}_hot) across a wide range of star formation rates (SFRs), from dwarf starbursts to ultra-luminous infrared galaxies. We show that for fixed thermalization efficiency and mass loading rate, the X-ray luminosity of the hot wind scales as L_X ~ SFR^2, significantly steeper than is observed for star-forming galaxies: L_X ~ SFR. Using this difference we constrain the mass-loading and thermalization efficiency of hot galactic winds. For reasonable values of the thermalization efficiency (<~ 1) and for SFR >~ 10 M_sun/yr we find that \dot{M}_hot/SFR <~ 1, significantly lower than required by integrated constraints on the efficiency of stellar feedback in galaxies, and potentially too low to explain observations of winds from rapidly star-forming galaxies. In addition, we highlight the fact that heavily mass-loaded winds cannot be described by the adiabatic CC85 model because they become strongly radiative.

preprint2014arXiv

Linear Instabilities Driven by Differential Rotation in Very Weakly Magnetized Plasmas

We study the linear stability of weakly magnetized differentially rotating plasmas in both collisionless kinetic theory and Braginskii's theory of collisional, magnetized plasmas. We focus on the very weakly magnetized limit that is important for understanding how astrophysical magnetic fields originate and are amplified at high redshift. We show that the single instability of fluid theory - the magnetorotational instability mediated by magnetic tension - is replaced by two distinct instabilities, one associated with ions and one with electrons. Each of these has a different way of tapping into the free energy of differential rotation. The ion instability is driven by viscous transport of momentum across magnetic field lines due to a finite ion cyclotron frequency (gyroviscosity); the fastest growing modes have wavelengths significantly longer than MHD and Hall MHD predictions. The electron instability is a whistler mode driven unstable by the temperature anisotropy generated by differential rotation; the growth time can be orders of magnitude shorter than the rotation period. The electron instability is an example of a broader class of instabilities that tap into the free energy of differential rotation or shear via the temperature anisotropy they generate. We briefly discuss the application of our results to the stability of planar shear flows and show that such flows are linearly overstable in the presence of fluid gyroviscosity. We also briefly describe the implications of our results for magnetic field amplification in the virialized halos of high redshift galaxies.

preprint2014arXiv

Magnetized Gas Clouds can Survive Acceleration by a Hot Wind

We present three-dimensional magnetohydrodynamic simulations of magnetized gas clouds accelerated by hot winds. We initialize gas clouds with tangled internal magnetic fields and show that this field suppresses the disruption of the cloud: rather than mixing into the hot wind as found in hydrodynamic simulations, cloud fragments end up co-moving and in pressure equilibrium with their surroundings. We also show that a magnetic field in the hot wind enhances the drag force on the cloud by a factor ~(1+v_A^2/v_wind^2)$, where v_A is the Alfven speed in the wind and v_wind measures the relative speed between the cloud and the wind. We apply this result to gas clouds in several astrophysical contexts, including galaxy clusters, galactic winds, the Galactic center, and the outskirts of the Galactic halo. Our results can explain the prevalence of cool gas in galactic winds and galactic halos and how such cool gas survives in spite of its interaction with hot wind/halo gas. We also predict that drag forces can lead to a deviation from Keplerian orbits for the G2 cloud in the galactic center.

preprint2014arXiv

Observational Signatures of Galactic Winds Powered by Active Galactic Nuclei

We predict the observational signatures of galaxy scale outflows powered by active galactic nuclei (AGN). Most of the emission is produced by the forward shock driven into the ambient interstellar medium (ISM) rather than by the reverse shock. AGN powered galactic winds with energetics suggested by phenomenological feedback arguments should produce spatially extended 1-10 keV X-ray emission of 10^(41-44) erg/s, significantly in excess of the spatially extended X-ray emission associated with normal star forming galaxies. The presence of such emission is a direct test of whether AGN outflows significantly interact with the ISM of their host galaxy. We further show that even radio quiet quasars should have a radio luminosity comparable to or in excess of the far infrared-radio correlation of normal star forming galaxies. This radio emission directly constrains the total kinetic energy flux in AGN-powered galactic winds. Radio emission from AGN wind shocks can also explain the recently highlighted correlations between radio luminosity and the kinematics of AGN narrow-line regions in radio quiet quasars.

preprint2014arXiv

Outflows from accretion disks formed in neutron star mergers: effect of black hole spin

The accretion disk that forms after a neutron star merger is a source of neutron-rich ejecta. The ejected material contributes to a radioactively-powered electromagnetic transient, with properties that depend sensitively on the composition of the outflow. Here we investigate how the spin of the black hole remnant influences mass ejection on the thermal and viscous timescales. We carry out two-dimensional, time-dependent hydrodynamic simulations of merger remnant accretion disks including viscous angular momentum transport and approximate neutrino self-irradiation. The gravity of the spinning black hole is included via a pseudo-Newtonian potential. We find that a disk around a spinning black hole ejects more mass, up to a factor of several, relative to the non-spinning case. The enhanced mass loss is due to energy release by accretion occurring deeper in the gravitational potential, raising the disk temperature and hence the rate of viscous heating in regions where neutrino cooling is ineffective. The mean electron fraction of the outflow increases moderately with BH spin due to a highly-irradiated (though not neutrino-driven) wind component. While the bulk of the ejecta is still very neutron-rich, thus generating heavy r-process elements, the leading edge of the wind contains a small amount of Lanthanide-free material. This component can give rise to a ~1 day blue optical `bump' in a kilonova light curve, even in the case of prompt BH formation, which may facilitate its detection.

preprint2014arXiv

The Role of Magnetic Fields in Starburst Galaxies as Revealed by OH Megamasers

We present estimates of magnetic field strengths in the interstellar media of starburst galaxies derived from measurements of Zeeman splitting associated with OH megamasers. The results for eight galaxies with Zeeman detections suggest that the magnetic energy density in the interstellar medium of starburst galaxies is comparable to their hydrostatic gas pressure, as in the Milky Way. We discuss the significant uncertainties in this conclusion, and possible measurements that could reduce these uncertainties. We also compare the Zeeman splitting derived magnetic field estimates to magnetic field strengths estimated using synchrotron fluxes and assuming that the magnetic field and cosmic rays have comparable energy densities, known as the "minimum energy" argument. We find that the minimum energy argument systematically underestimates magnetic fields in starburst galaxies, and that the conditions that would be required to produce agreement between the minimum energy estimate and the Zeeman derived estimate of interstellar medium magnetic fields are implausible. The conclusion that magnetic fields in starburst galaxies exceed the minimum energy magnetic fields is consistent with starburst galaxies adhering to the linearity of the FIR-radio correlation.

preprint2013arXiv

Dense Molecular Gas: A Sensitive Probe of Stellar Feedback Models

We show that the mass fraction of GMC gas (n>100 cm^-3) in dense (n>>10^4 cm^-3) star-forming clumps, observable in dense molecular tracers (L_HCN/L_CO(1-0)), is a sensitive probe of the strength and mechanism(s) of stellar feedback. Using high-resolution galaxy-scale simulations with pc-scale resolution and explicit models for feedback from radiation pressure, photoionization heating, stellar winds, and supernovae (SNe), we make predictions for the dense molecular gas tracers as a function of GMC and galaxy properties and the efficiency of stellar feedback. In models with weak/no feedback, much of the mass in GMCs collapses into dense sub-units, predicting L_HCN/L_CO(1-0) ratios order-of-magnitude larger than observed. By contrast, models with feedback properties taken directly from stellar evolution calculations predict dense gas tracers in good agreement with observations. Changing the strength or timing of SNe tends to move systems along, rather than off, the L_HCN-L_CO relation (because SNe heat lower-density material, not the high-density gas). Changing the strength of radiation pressure (which acts efficiently in the highest density gas), however, has a much stronger effect on L_HCN than on L_CO. We predict that the fraction of dense gas (L_HCN/L_CO(1-0)) increases with increasing GMC surface density; this drives a trend in L_HCN/L_CO(1-0) with SFR and luminosity which has tentatively been observed. Our results make specific predictions for enhancements in the dense gas tracers in unusually dense environments such as ULIRGs and galactic nuclei (including the galactic center).

preprint2013arXiv

Internal Gravity Wave Excitation by Turbulent Convection

We calculate the flux of internal gravity waves (IGWs) generated by turbulent convection in stars. We solve for the IGW eigenfunctions analytically near the radiative-convective interface in a local, Boussinesq, and cartesian domain. We consider both discontinuous and smooth transitions between the radiative and convective regions and derive Green's functions to solve for the IGWs in the radiative region. We find that if the radiative-convective transition is smooth, the IGW flux depends on the exact form of the buoyancy frequency near the interface. IGW excitation is most efficient for very smooth interfaces, which gives an upper bound on the IGW flux of ~ F_conv (d/H), where F_conv is the flux carried by the convective motions, d is the width of the transition region, and H is the pressure scale height. This can be much larger than the standard result in the literature for a discontinuous radiative-convective transition, which gives a wave flux ~ F_conv M$, where M is the convective Mach number. However, in the smooth transition case, the most efficiently excited perturbations will break in the radiative zone. The flux of IGWs which do not break and are able to propagate in the radiative region is at most F_conv M^(5/8) (d/H)^(3/8), larger than the discontinuous transition result by (MH/d)^(-3/8). The transition region in the Sun is smooth for the energy-bearing waves; as a result, we predict that the IGW flux is a few to five times larger than previous estimates. We discuss the implications of our results for several astrophysical applications, including IGW driven mass loss and the detectability of convectively excited IGWs in main sequence stars.

preprint2013arXiv

Perpendicular Ion Heating by Reduced Magnetohydrodynamic Turbulence

Recent theoretical studies argue that the rate of stochastic ion heating in low-frequency Alfvén-wave turbulence is given by $Q_\perp = c_1 [(δu)^3 /ρ] \exp(-c_2/ε)$, where $δu$ is the rms turbulent velocity at the scale of the ion gyroradius $ρ$, $ε= δu/v_{\perp \rm i}$, $v_{\perp \rm i}$ is the perpendicular ion thermal speed, and $c_1$ and $c_2$ are dimensionless constants. We test this theoretical result by numerically simulating test particles interacting with strong reduced magnetohydrodynamic (RMHD) turbulence. The heating rates in our simulations are well fit by this formula. The best-fit values of $c_1$ are $\sim 1$. The best-fit values of $c_2$ decrease (i.e., stochastic heating becomes more effective) as the grid size and Reynolds number of the RMHD simulations increase. As an example, in a $1024^2 \times 256$ RMHD simulation with a dissipation wavenumber of order the inverse ion gyroradius, we find $c_2 = 0.21$. We show that stochastic heating is significantly stronger in strong RMHD turbulence than in a field of randomly phased Alfvén waves with the same power spectrum, because coherent structures in strong RMHD turbulence increase orbit stochasticity in the regions where ions are heated most strongly. We find that $c_1$ increases by a factor of $\sim 3$ while $c_2$ changes very little as the ion thermal speed increases from values $\ll v_{\rm A}$ to values $\sim v_{\rm A}$, where $v_{\rm A}$ is the Alfvén speed. We discuss the importance of these results for perpendicular ion heating in the solar wind.

preprint2013arXiv

Setting the Stage for Circumstellar Interaction in Core-Collapse Supernovae II: Wave-Driven Mass Loss in Supernova Progenitors

Supernovae (SNe) powered by interaction with circumstellar material provide evidence for intense stellar mass loss during the final years leading up to core collapse. We have argued that during and after core neon burning, internal gravity waves excited by core convection can tap into the core fusion power and transport a super-Eddington energy flux out to the stellar envelope, potentially unbinding up to ~ 1 solar mass of material. In this work, we explore the internal conditions of SN progenitors using the MESA 1-D stellar evolution code, in search of those most susceptible to wave-driven mass loss. We focus on simple, order of magnitude considerations applicable to a wide range of progenitors. Wave-driven mass loss during core neon and oxygen fusion happens preferentially in either lower mass (<~ 20 solar mass ZAMS) stars or massive, sub-solar metallicity stars. Roughly 20 per cent of the SN progenitors we survey can excite ~ 10^46 - 10^48 erg of energy in waves that can potentially drive mass loss within a few months to a decade of core collapse. This energy can generate a circumstellar environment with 10^-3 - 1 solar mass reaching ~ 100 AU before explosion. We predict a correlation between the energy associated with pre-SN mass ejection and the time to core collapse, with the most intense mass loss preferentially happening closer to core collapse. During silicon burning, a ~ 5 day long phase for our progenitor models, wave energy may inflate ~ 10^-3 - 1 solar mass of the stellar envelope to ~ 10 - 100s of solar radii. This suggests that some nominally compact SN progenitors (Type Ibc progenitors) will have a significantly different SN shock breakout signature than traditionally assumed. We discuss the implications of our results for the core-collapse SN mechanism, Type IIn SNe, Type IIb SNe from extended progenitors (e.g., SNe 1993j and 2011dh), and observed pre-SN outbursts.

preprint2013arXiv

The Efficiency of Second-Order Fermi Acceleration by Weakly Compressible MHD Turbulence

We investigate the effects of pitch-angle scattering on the efficiency of particle heating and acceleration by MHD turbulence using phenomenological estimates and simulations of non-relativistic test particles interacting with strong, subsonic MHD turbulence. We include an imposed pitch-angle scattering rate, which is meant to approximate the effects of high frequency plasma waves and/or velocity space instabilities. We focus on plasma parameters similar to those found in the near-Earth solar wind, though most of our results are more broadly applicable. An important control parameter is the size of the particle mean free path lambda_{mfp} relative to the scale of the turbulent fluctuations L. For small scattering rates, particles interact quasi-resonantly with turbulent fluctuations in magnetic field strength. Scattering increases the long-term efficiency of this resonant heating by factors of a few-10, but the distribution function does not develop a significant non-thermal power-law tail. For higher scattering rates, the interaction between particles and turbulent fluctuations becomes non-resonant, governed by particles heating and cooling adiabatically as they encounter turbulent density fluctuations. Rapid pitch-angle scattering can produce a power-law tail in the proton distribution function but this requires fine-tuning of parameters. Moreover, in the near-Earth solar wind, a significant power-law tail cannot develop by this mechanism because the particle acceleration timescales are longer than the adiabatic cooling timescale set by the expansion of the solar wind. Our results thus imply that MHD-scale turbulent fluctuations are unlikely to be the origin of the v^{-5} tail in the proton distribution function observed in the solar wind.

preprint2013arXiv

Tidal resonance locks in inspiraling white dwarf binaries

We calculate the tidal response of helium and carbon/oxygen (C/O) white dwarf (WD) binaries inspiraling due to gravitational wave emission. We show that resonance locks, previously considered in binaries with an early-type star, occur universally in WD binaries. In a resonance lock, the orbital and spin frequencies evolve in lockstep, so that the tidal forcing frequency is approximately constant and a particular normal mode remains resonant, producing efficient tidal dissipation and nearly synchronous rotation. We show that analogous locks between the spin and orbital frequencies can occur not only with global standing modes, but even when damping is so efficient that the resonant tidal response becomes a traveling wave. We derive simple analytic formulas for the tidal quality factor Q and tidal heating rate during a g-mode resonance lock, and verify our results numerically. We find that Q ~ 10^7 for orbital periods ~ 1 - 2 hr in C/O WDs, and Q ~ 10^9 for P_orb ~ 3 - 10 hr in helium WDs. Typically tidal heating occurs sufficiently close to the surface that the energy should be observable as surface emission. Moreover, near an orbital period of ~ 10 min, the tidal heating rate reaches ~ 10^{-2} L_\sun, rivaling the luminosities of our fiducial WD models. Recent observations of the 13-minute double-WD binary J0651 are roughly consistent with our theoretical predictions. Tides naturally tend to generate differential rotation; however, we show that the fossil magnetic field strength of a typical WD can maintain solid-body rotation down to at least P_orb ~ 10 min even in the presence of a tidal torque concentrated near the WD surface.

preprint2012arXiv

3-D Radiative Transfer Calculations of Radiation Feedback from Massive Black Holes: Outflow of Mass from the Dusty "Torus"

Observational and theoretical arguments suggest that the momentum carried in mass outflows from AGN can reach several times L / c, corresponding to outflow rates of hundreds of solar masses per year. Radiation pressure on lines alone may not be sufficient to provide this momentum deposition, and the transfer of reprocessed IR radiation in dusty nuclear gas has been postulated to provide the extra enhancement. The efficacy of this mechanism, however, will be sensitive to multi-dimensional effects such as the tendency for the reprocessed radiation to preferentially escape along sight-lines of lower column density. We use Monte Carlo radiative transfer calculations to determine the radiation force on dusty gas residing within approximately 10 parsecs from an accreting super-massive black hole. We calculate the net rate of momentum deposition in the surrounding gas and estimate the mass-loss rate in the resulting outflow as a function of solid angle for different black hole luminosities, sightline-averaged column densities, clumping parameters, and opening angles of the dusty gas. We find that these dust-driven winds carry momentum fluxes of 1-5 times L / c and correspond to mass-loss rates of 10-100 solar masses per year for a 10^8 solar mass black hole radiating at or near its Eddington limit. These results help to explain the origin of high velocity molecular and atomic outflows in local ULIRGs, and can inform numerical simulations of galaxy evolution including AGN feedback.

preprint2012arXiv

An Analytic Model of Angular Momentum Transport by Gravitational Torques: From Galaxies to Massive Black Holes

We present analytic calculations of angular momentum transport and gas inflow in galaxies, from scales of ~kpc to deep in the potential of a central black hole (BH). We compare these analytic calculations to numerical simulations and use them to develop a sub-grid model of BH growth that can be incorporated into semi-analytic models or cosmological simulations. Both analytic calculations and simulations argue that the strongest torque on gas arises when non-axisymmetric perturbations to the stellar gravitational potential produces orbit crossings and shocks in the gas. This is true both at large radii ~0.01-1 kpc, where bar-like modes dominate the non-axisymmetric potential, and at smaller radii <10 pc, where a lopsided/eccentric disk dominates. The traditional orbit crossing criterion is not always adequate to predict the locations of, and inflow due to, shocks in gas+stellar disks with finite sound speeds. We derive a modified criterion that predicts the presence of shocks in stellar dominated systems even absent formal orbit crossing. We then derive analytic expressions for the loss of angular momentum and the resulting gas inflow rates in the presence of such shocks. We test our analytic predictions using hydrodynamic simulations at a range of galactic scales, and show that they successfully predict the mass inflow rates and quasi-steady gas surface densities with small scatter (0.3 dex). We use our analytic results to construct a new estimate of the BH accretion rate given galaxy properties at larger radii. This captures the key scalings in the numerical simulations. Alternate estimates such as the local viscous accretion rate or the spherical Bondi rate fail systematically to reproduce the simulations.

preprint2012arXiv

An Explanation for the Slopes of Stellar Cusps in Galaxy Spheroids

The stellar surface mass density profiles at the centers of typical ~L* and lower-mass spheroids exhibit power law 'cusps' with $Σ\propto R^(-n)$, where 0.5<n<1 for radii ~1-100 pc. Observations and theory support models in which these cusps are formed by dissipative gas inflows and nuclear starbursts in gas-rich mergers. At these comparatively large radii, stellar relaxation is unlikely to account for or strongly modify the cuspy stellar profiles. We argue that the power-law surface density profiles observed are a natural consequence of the gravitational instabilities that dominate angular momentum transport in the gravitational potential of a central massive black hole. The dominant mode at these radii is an m=1 lopsided/eccentric disk instability, in which stars torquing the gas can drive rapid inflow and accretion. Such a mode first generically appears at large radii and propagates inwards by exciting eccentricities at smaller and smaller radii, where M*(<R)<<M_BH. When the stellar surface density profile is comparatively shallow with n<1/2, the modes cannot efficiently propagate to R=0 and so gas piles up and star formation steepens the profile. But if the profile is steeper than n=1, the inwards propagation of eccentricity is strongly damped, suppressing inflow and bringing n down again. Together these results produce an equilibrium slope of 1/2 < n < 1 in the potential of the central black hole. These physical arguments are supported by nonlinear numerical simulations of gas inflow in galactic nuclei. Together, these results naturally explain the observed stellar density profiles of 'cusp' elliptical galaxies.

preprint2012arXiv

Inhomogeneous accretion discs and the soft states of black hole X-ray binaries

Observations of black hole binaries (BHBs) have established a rich phenomenology of X-ray states. The soft states range from the low variability, accretion disc dominated thermal state (TD) to the higher variability, non-thermal steep power law state (SPL). The disc component in all states is typically modeled with standard thin disc accretion theory. However, this theory is inconsistent with optical/UV spectral, variability, and gravitational microlensing observations of active galactic nuclei (AGNs), the supermassive analogs of BHBs. An inhomogeneous disc (ID) model with large (~0.4 dex) temperature fluctuations in each radial annulus can qualitatively explain all of these AGN observations. The inhomogeneity may be a consequence of instabilities in radiation dominated discs, and therefore may be present in BHBs as well. We show that ID models can explain many features of the TD and SPL states of BHBs. The observed relationships between spectral hardness, disc fraction, and rms variability amplitude in BHBs are reproduced with temperature fluctuations similar to those inferred in AGNs, suggesting a unified picture of luminous accretion discs across orders of magnitude in black hole mass. This picture can be tested with spectral fitting of ID models, X-ray polarization observations, and radiation MHD simulations. If BHB accretion discs are indeed inhomogeneous, only the most disc dominated states (disc fraction > 0.95) can be used to robustly infer black hole spin using current continuum fitting methods.

preprint2012arXiv

Local 2D Particle-in-cell simulations of the collisionless MRI

The magnetorotational instability (MRI) is a crucial mechanism of angular momentum transport in a variety of astrophysical accretion disks. In systems accreting at well below the Eddington rate, such as the central black hole in the Milky Way (Sgr A*), the rate of Coulomb collisions between particles is very small, making the disk evolve essentially as a collisionless plasma. We present a nonlinear study of the collisionless MRI using first-principles particle-in-cell (PIC) plasma simulations. In this initial study we focus on local two-dimensional (axisymmetric) simulations, deferring more realistic three-dimensional simulations to future work. For simulations with net vertical magnetic flux, the MRI continuously amplifies the magnetic field until the Alfvén velocity, v_A, is comparable to the speed of light, c (independent of the initial value of v_A/c). This is consistent with the lack of saturation of MRI channel modes in analogous axisymmetric MHD simulations. The amplification of the magnetic field by the MRI generates a significant pressure anisotropy in the plasma (with the perpendicular pressure being larger than the parallel pressure). We find that this pressure anisotropy in turn excites mirror modes and that the volume averaged pressure anisotropy remains near the threshold for mirror mode excitation. Particle energization is due to both reconnection and viscous heating associated with the pressure anisotropy. Reconnection produces a distinctive power-law component in the energy distribution function of the particles, indicating the likelihood of non-thermal ion and electron acceleration in collisionless accretion disks. This has important implications for interpreting the observed emission -- from the radio to the gamma-rays -- of systems such as Sgr A*.

preprint2012arXiv

Nonlinear Tides in Close Binary Systems

We study the excitation and damping of tides in close binary systems, accounting for the leading order nonlinear corrections to linear tidal theory. These nonlinear corrections include two distinct effects: three-mode nonlinear interactions and nonlinear excitation of modes by the time-varying gravitational potential of the companion. This paper presents the formalism for studying nonlinear tides and studies the nonlinear stability of the linear tidal flow. Although the formalism is applicable to binaries containing stars, planets, or compact objects, we focus on solar type stars with stellar or planetary companions. Our primary results include: (1) The linear tidal solution often used in studies of binary evolution is unstable over much of the parameter space in which it is employed. More specifically, resonantly excited gravity waves are unstable to parametric resonance for companion masses M' > 10-100 M_Earth at orbital periods P = 1-10 days. The nearly static equilibrium tide is, however, parametrically stable except for solar binaries with P < 2-5 days. (2) For companion masses larger than a few Jupiter masses, the dynamical tide causes waves to grow so rapidly that they must be treated as traveling waves rather than standing waves. (3) We find a novel form of parametric instability in which a single parent wave excites a very large number of daughter waves (N = 10^3[P / 10 days]) and drives them as a single coherent unit with growth rates that are ~N times faster than the standard three wave parametric instability. (4) Independent of the parametric instability, tides excite a wide range of stellar p-modes and g-modes by nonlinear inhomogeneous forcing; this coupling appears particularly efficient at draining energy out of the dynamical tide and may be more important than either wave breaking or parametric resonance at determining the nonlinear dissipation of the dynamical tide.

preprint2012arXiv

On the Structure of Hot Gas in Halos: Implications for the Lx-Tx Relation & Missing Baryons

We present one-dimensional models of the hot gas in dark-matter halos, which both predict the existence of cool cores and explain their structure. Our models are directly applicable to semi-analytic models (SAMs) of galaxy formation. We have previously argued that filaments of cold (~10^4 K) gas condense out of the intracluster medium (ICM) in hydrostatic and thermal equilibrium when the ratio of the thermal instability timescale to the free-fall time $t_{TI}/t_{ff}$ falls below 5-10. This criterion corresponds to an upper limit on the density of the ICM and motivates a model in which a density core forms wherever $t_{TI}/t_{ff} \lesssim 10$. Consistent with observations and numerical simulations, this model predicts larger and more tenuous cores for lower-mass halos---while the core density in a cluster may be as large as ~ 0.1 cm^{-3}, the core density in the Galactic halo should not exceed ~ 10^{-4} cm^{-3}. We can also explain the large densities in smaller mass halos (galactic `coronae') if we include the contribution of the central galaxy to the gravitational potential. Our models produce a favorable match to the observational X-ray luminosity-temperature (Lx-Tx) relation. For halo masses $\lesssim 10^{13}$ Msun the core size approaches the virial radius. Thus, most of the baryons in such halos cannot be in the hot ICM, but either in the form of stars or in the form of hot gas beyond the virial radius. Because of the smaller mass in the ICM and much larger mass available for star formation, the majority of the baryons in low mass halos ($\lesssim 10^{13}$ Msun) can be expelled beyond the virial radius due to supernova feedback. This can account for the baryons `missing' from low mass halos, such as the Galactic halo.

preprint2012arXiv

Resonance Broadening and Heating of Charged Particles in Magnetohydrodynamic Turbulence

The heating, acceleration, and pitch-angle scattering of charged particles by MHD turbulence are important in a wide range of astrophysical environments, including the solar wind, accreting black holes, and galaxy clusters. We simulate the interaction of high-gyrofrequency test particles with fully dynamical simulations of subsonic MHD turbulence, focusing on the parameter regime with beta ~ 1, where beta is the ratio of gas to magnetic pressure. We use the simulation results to calibrate analytical expressions for test particle velocity-space diffusion coefficients and provide simple fits that can be used in other work. The test particle velocity diffusion in our simulations is due to a combination of two processes: interactions between particles and magnetic compressions in the turbulence (as in linear transit-time damping; TTD) and what we refer to as Fermi Type-B (FTB) interactions, in which charged particles moving on field lines may be thought of as beads spiralling around moving wires. We show that test particle heating rates are consistent with a TTD resonance which is broadened according to a decorrelation prescription that is Gaussian in time. TTD dominates the heating for v_s >> v_A (e.g. electrons), where v_s is the thermal speed of species s and v_A is the Alfven speed, while FTB dominates for v_s << v_A (e.g. minor ions). Proton heating rates for beta ~ 1 are comparable to the turbulent cascade rate. Finally, we show that velocity diffusion of collisionless, large gyrofrequency particles due to large-scale MHD turbulence does not produce a power-law distribution function.

preprint2012arXiv

Stellar Feedback & Bulge Formation in Clumpy Disks

We use numerical simulations of isolated galaxies to study the effects of stellar feedback on the formation and evolution of giant star-forming gas 'clumps' in high-redshift, gas-rich galaxies. Such galactic disks are unstable to the formation of bound gas-rich clumps whose properties initially depend only on global disk properties, not the microphysics of feedback. In simulations without stellar feedback, clumps turn an order-unity fraction of their mass into stars and sink to the center, forming a large bulge and kicking most of the stars out into a much more extended stellar envelope. By contrast, strong radiative stellar feedback disrupts even the most massive clumps after they turn ~10-20% of their mass into stars, in a timescale of ~10-100 Myr, ejecting some material into a super-wind and recycling the rest of the gas into the diffuse ISM. This suppresses the bulge formation rate by direct 'clump coalescence' by a factor of several. However, the galactic disks do undergo significant internal evolution in the absence of mergers: clumps form and disrupt continuously and torque gas to the galactic center. The resulting evolution is qualitatively similar to bar/spiral evolution in simulations with a more homogeneous ISM.

preprint2012arXiv

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 super-winds that enrich the IGM and shape the galaxy mass function and mass-metallicity relation. In previous papers, we introduced new numerical methods for implementing stellar feedback on sub-GMC through galactic scales in galaxy simulations. This includes radiation pressure (UV through IR), SNe (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 SFR. The mass-loading efficiency (wind mass loss rate divided by SFR) scales inversely with circular velocity, consistent with momentum-conservation expectations. We study the contributions of each feedback mechanism to galactic winds in a range of galaxy models, from SMC-like dwarfs & MW-analogues to z~2 clumpy disks. In massive, gas-rich systems (local starbursts and high-z galaxies), radiation pressure dominates the wind generation. For MW-like spirals and dwarf galaxies the gas densities are much lower, and shock-heated gas from SNe and stellar winds dominates production of large-scale outflows. In all models, however, winds have a multi-phase structure that depends on interactions between multiple feedback mechanisms operating on different spatial & time scales: any single mechanism fails to reproduce the winds observed. We provide fitting functions for 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 explicit dependence on gas surface density that can be very important in both low-density dwarf galaxies and high-density gas-rich galaxies.

preprint2012arXiv

The Dynamics of Rayleigh-Taylor Stable and Unstable Contact Discontinuities with Anisotropic Thermal Conduction

We study the effects of anisotropic thermal conduction along magnetic field lines on an accelerated contact discontinuity in a weakly collisional plasma. We first perform a linear stability analysis similar to that used to derive the Rayleigh-Taylor instability (RTI) dispersion relation. We find that anisotropic conduction is only important for compressible modes, as incompressible modes are isothermal. Modes grow faster in the presence of anisotropic conduction, but growth rates do not change by more than a factor of order unity. We next run fully non-linear numerical simulations of a contact discontinuity with anisotropic conduction. The non-linear evolution can be thought of as a superposition of three physical effects: temperature diffusion due to vertical conduction, the RTI, and the heat flux driven buoyancy instability (HBI). In simulations with RTI-stable contact discontinuities, the temperature discontinuity spreads due to vertical heat conduction. This occurs even for initially horizontal magnetic fields due to the initial vertical velocity perturbation and numerical mixing across the interface. The HBI slows this temperature diffusion by reorienting initially vertical magnetic field lines to a more horizontal geometry. In simulations with RTI-unstable contact discontinuities, the dynamics are initially governed by temperature diffusion, but the RTI becomes increasingly important at late times. We discuss the possible application of these results to supernova remnants, solar prominences, and cold fronts in galaxy clusters.

preprint2012arXiv

The Effects of Anisotropic Viscosity on Turbulence and Heat Transport in the Intracluster Medium

In the intracluster medium (ICM) of galaxy clusters, heat and momentum are transported almost entirely along (but not across) magnetic field lines. We perform the first fully self-consistent Braginskii-MHD simulations of galaxy clusters including both of these effects. Specifically, we perform local and global simulations of the magnetothermal instability (MTI) and the heat-flux-driven buoyancy instability (HBI) and assess the effects of viscosity on their saturation and astrophysical implications. We find that viscosity has only a modest effect on the saturation of the MTI. As in previous calculations, we find that the MTI can generate nearly sonic turbulent velocities in the outer parts of galaxy clusters, although viscosity somewhat suppresses the magnetic field amplification. At smaller radii in cool-core clusters, viscosity can decrease the linear growth rates of the HBI. However, it has less of an effect on the HBI's nonlinear saturation, in part because three-dimensional interchange motions (magnetic flux tubes slipping past each other) are not damped by anisotropic viscosity. In global simulations of cool core clusters, we show that the HBI robustly inhibits radial thermal conduction and thus precipitates a cooling catastrophe. The effects of viscosity are, however, more important for higher entropy clusters. We argue that viscosity can contribute to the global transition of cluster cores from cool-core to non cool-core states: additional sources of intracluster turbulence, such as can be produced by AGN feedback or galactic wakes, suppress the HBI, heating the cluster core by thermal conduction; this makes the ICM more viscous, which slows the growth of the HBI, allowing further conductive heating of the cluster core and a transition to a non cool-core state.

preprint2012arXiv

The Long-Term Evolution of Double White Dwarf Mergers

In this paper, we present a model for the long-term evolution of the merger of two unequal mass C/O white dwarfs (WDs). After the dynamical phase of the merger, magnetic stresses rapidly redistribute angular momentum, leading to nearly solid-body rotation on a viscous timescale of 1e-4 to 1 yr, long before significant cooling can occur. Due to heating during the dynamical and viscous phases, the less massive WD is transformed into a hot, slowly rotating, and radially extended envelope supported by thermal pressure. Following the viscous phase of evolution, the maximum temperature near the envelope base may already be high enough to begin off-center convective carbon-burning. If not, Kelvin-Helmholtz contraction of the inner region of the envelope on a thermal timescale of 1e3-1e4 yr compresses the base of the envelope, again yielding off-center burning. As a result, the long-term evolution of the merger remnant is similar to that seen in previous calculations: the burning shell diffuses inwards over ~1e4 yr, eventually yielding a high-mass O/Ne WD or a collapse to a neutron star. During the cooling and shell-burning phases, the merger remnant radiates near the Eddington limit. Given the double WD merger rate of a few per 1000 yr, tens of these ~1e38 erg/s sources should exist in a Milky Way-type galaxy. While the end result is similar to that of previous studies, the physical picture and the dynamical state of the matter in our model differ from previous work. Furthermore, remaining uncertainties related to the convective structure near the photosphere and mass loss during the thermal evolution may significantly affect our conclusions. Thus, future work within the context of the physical model presented here is required to better address the eventual fate of double WD mergers, including those for which one or both of the components is a He WD.

preprint2012arXiv

The Observational Signatures of Convectively Excited Internal Gravity Modes in Main Sequence Stars

We predict the flux and surface velocity perturbations produced by convectively excited gravity modes (g-modes) in main sequence stars. Core convection in massive stars can excite g-modes to sufficient amplitudes to be detectable with high precision photometry by Kepler and CoRoT, if the thickness of the convective overshoot region is < 30 per cent of a pressure scale height. The g-modes manifest as excess photometric variability, with amplitudes of ~ 10 micromagnitudes at frequencies ~ 10 microHz (0.8 1/d) near the solar metallicity zero-age main sequence. The flux variations are largest for stars with M > 5 solar masses, but are potentially detectable down to M ~ 2 - 3 solar masses. During the main sequence evolution, radiative damping decreases such that ever lower frequency modes reach the stellar surface and flux perturbations reach up to ~ 100 micromagnitudes at the terminal-age main sequence. Using the same convective excitation model, we confirm previous predictions that solar g-modes produce surface velocity perturbations of < 0.3 mm/s. This implies that stochastically excited g-modes are more easily detectable in the photometry of massive main sequence stars than in the Sun.

preprint2012arXiv

The Stability of Massive Main Sequence Stars as a Function of Metallicity

We investigate the pulsational stability of massive (M >~ 120 Msun) main sequence stars of a range of metallicities, including primordial, Population III stars. We include a formulation of convective damping motivated by numerical simulations of the interaction between convection and periodic shear flows. We find that convective viscosity is likely strong enough to stabilize radial pulsations whenever nuclear-burning (the epsilon-mechanism) is the dominant source of driving. This suggests that massive main sequence stars with Z <~ 2 x 10^-3 are pulsationally stable and are unlikely to experience pulsation-driven mass loss on the main sequence. These conclusions are, however, sensitive to the form of the convective viscosity and highlight the need for further high-resolution simulations of the convection-oscillation interaction. For more metal-rich stars (Z >~ 2 x 10^-3), the dominant pulsational driving arises due to the kappa-mechanism arising from the iron-bump in opacity and is strong enough to overcome convective damping. Our results highlight that even for oscillations with periods a few orders of magnitude shorter than the outer convective turnover time, the "frozen-in" approximation for the convection-oscillation interaction is inappropriate, and convective damping should be taken into account when assessing mode stability.

preprint2012arXiv

The Structure of the Interstellar Medium of Star Forming Galaxies

We present numerical methods for including stellar feedback in galaxy-scale simulations. We include heating by SNe (I & II), gas recycling and shock-heating from O-star & AGB winds, HII photoionization, and radiation pressure from stellar photons. The energetics and time-dependence are taken directly from stellar evolution models. We implement these in simulations with pc-scale resolution, modeling galaxies from SMC-like dwarfs and MW analogues to massive z~2 starburst disks. Absent feedback, gas cools and collapses without limit. With feedback, the ISM reaches a multi-phase steady state in which GMCs continuously form, disperse, and re-form. Our primary results include: (1) Star forming galaxies generically self-regulate at Toomre Q~1. Most of the volume is in diffuse hot gas with most of the mass in dense GMC complexes. The phase structure and gas mass at high densities are much more sensitive probes of stellar feedback physics than integrated quantities (Toomre Q or gas velocity dispersion). (2) Different feedback mechanisms act on different scales: radiation & HII pressure are critical to prevent runaway collapse of dense gas in GMCs. SNe and stellar winds dominate the dynamics of volume-filling hot gas; however this primarily vents out of the disk. (3) The galaxy-averaged SFR is determined by feedback. For given feedback efficiency, restricting star formation to molecular gas or modifying the cooling function has little effect; but changing feedback mechanisms directly translates to shifts off the Kennicutt-Schmidt relation. (4) Self-gravity leads to marginally-bound GMCs with an ~M^-2 mass function with a cutoff at the Jeans mass; they live a few dynamical times before being disrupted by stellar feedback and turn ~1-10% of their mass into stars (increasing from dwarfs through starburst galaxies). Low-mass GMCs are preferentially unbound.

preprint2012arXiv

The Viscous Evolution of White Dwarf Merger Remnants

The merger of two white dwarfs (WDs) creates a differentially rotating remnant which is unstable to magnetohydrodynamic instabilities. These instabilities can lead to viscous evolution on a time-scale short compared to the thermal evolution of the remnant. We present multi-dimensional hydrodynamic simulations of the evolution of WD merger remnants under the action of an $α$-viscosity. We initialize our calculations using the output of eight WD merger simulations from Dan et al. (2011), which span a range of mass ratios and total masses. We generically find that the merger remnants evolve towards spherical states on time-scales of hours, even though a significant fraction of the mass is initially rotationally supported. The viscous evolution unbinds only a very small amount of mass $(< 10^{-5} M_\odot)$. Viscous heating causes some of the systems we study with He WD secondaries to reach conditions of nearly dynamical burning. It is thus possible that the post-merger viscous phase triggers detonation of the He envelope in some WD mergers, potentially producing a Type Ia supernova via a double detonation scenario. Our calculations provide the proper initial conditions for studying the long-term thermal evolution of WD merger remnants. This is important for understanding WD mergers as progenitors of Type Ia supernovae, neutron stars, R Coronae Borealis stars and other phenomena.

preprint2012arXiv

Thermal Instability in Gravitationally-Stratified Plasmas: Implications for Multi-Phase Structure in Clusters and Galaxy Halos

We study the interplay among cooling, heating, conduction, and magnetic fields in gravitationally stratified plasmas using simplified, plane-parallel numerical simulations. Since the physical heating mechanism remains uncertain in massive halos such as groups or clusters, we adopt a simple, observationally-motivated prescription which enforces global thermal equilibrium when averaged over large scales. The plasma remains susceptible to local thermal instability, however, and cooling drives an inward flow of material. In contrast to previous results, we argue that the thermal stability of the plasma is independent of its convective stability. We find that the ratio of the cooling timescale to the dynamical timescale t_cool/t_ff controls the saturation of the thermal instability: when t_cool/t_ff < 1, the plasma develops extended multi-phase structure, whereas when t_cool / t_ff > 1 it does not. (In a companion paper, we show that the criterion for thermal instability in a spherical potential is somewhat less stringent, t_cool / t_ff < 10.) When thermal conduction is anisotropic with respect to the magnetic field, the criterion for multi-phase structure is essentially independent of the thermal conductivity of the plasma. Our criterion for local thermal instability to produce multi-phase structure is an extension of the cold vs. hot accretion modes in galaxy formation that applies at all radii in hot halos, not just to the virial shock. We show that this criterion is consistent with data on multi-phase gas in the ACCEPT sample of clusters; in addition, when t_cool / t_ff > 1, the net cooling rate to low temperatures and the mass flux to small radii are suppressed enough relative to models without heating to be qualitatively consistent with star formation rates and x-ray line emission in groups and clusters.

preprint2012arXiv

Wave-Driven Mass Loss in the Last Year of Stellar Evolution: Setting the Stage for the Most Luminous Core-Collapse Supernovae

During the late stages of stellar evolution in massive stars (C fusion and later), the fusion luminosity in the core of the star exceeds the star's Eddington luminosity. This can drive vigorous convective motions which in turn excite internal gravity waves. The local wave energy flux excited by convection is itself well above Eddington during the last few years in the life of the star. We suggest that an interesting fraction of the energy in gravity waves can, in some cases, convert into sound waves as the gravity waves propagate (tunnel) towards the stellar surface. The subsequent dissipation of the sound waves can unbind up to several $M_\odot$ of the stellar envelope. This wave-driven mass loss can explain the existence of extremely large stellar mass loss rates just prior to core-collapse, which are inferred via circumstellar interaction in some core-collapse supernovae (e.g., SNe 2006gy and PTF 09uj, and even Type IIn supernovae more generally). An outstanding question is understanding what stellar parameters (mass, rotation, metallicity, age) are the most susceptible to wave-driven mass loss. This depends on the precise internal structure of massive stars and the power-spectrum of internal gravity waves excited by stellar convection.

preprint2011arXiv

A relativistic jetted outburst from a massive black hole fed by a tidally disrupted star

While gas accretion onto some massive black holes (MBHs) at the centers of galaxies actively powers luminous emission, the vast majority of MBHs are considered dormant. Occasionally, a star passing too near a MBH is torn apart by gravitational forces, leading to a bright panchromatic tidal disruption flare (TDF). While the high-energy transient Swift J164449.3+573451 ("Sw 1644+57") initially displayed none of the theoretically anticipated (nor previously observed) TDF characteristics, we show that the observations (Levan et al. 2011) suggest a sudden accretion event onto a central MBH of mass ~10^6-10^7 solar masses. We find 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 with a smaller-scale blazar. The phenomenologically novel Sw 1644+57 thus connects the study of TDFs and active galaxies, opening a new vista on disk-jet interactions in BHs and magnetic field generation and transport in accretion systems.

preprint2011arXiv

Galaxy-Scale Outflows Driven by Active Galactic Nuclei

We present hydrodynamical simulations of major mergers of galaxies and study the effects of winds produced by active galactic nuclei (AGN) on interstellar gas in the AGN's host galaxy. We consider winds with initial velocities ~ 10,000 km/s and an initial momentum (energy) flux of ~ tau_w L/c (~ 0.01 tau_w L), with tau_w ~ 1-10. The AGN wind sweeps up and shock heats the surrounding interstellar gas, leading to a galaxy-scale outflow with velocities ~ 1000 km/s, peak mass outflow rates comparable to the star formation rate, and a total ejected gas mass ~ 3 x 10^9 M_sun. Large momentum fluxes, tau_w > 3, are required for the AGN-driven galactic outflow to suppress star formation and accretion in the black hole's host galaxy. Less powerful AGN winds (tau_w < 3) still produce a modest galaxy-scale outflow, but the outflow has little global effect on the ambient interstellar gas. We argue that this mechanism of AGN feedback can plausibly produce the high velocity outflows observed in post-starburst galaxies and the massive molecular and atomic outflows observed in local ultra-luminous infrared galaxies. Moreover, the outflows from local ultra-luminous infrared galaxies are inferred to have tau_w ~ 10, comparable to what we find is required for AGN winds to regulate the growth of black holes and set the M_BH-sigma relation. We conclude by discussing theoretical mechanisms that can lead to AGN wind mass-loading and momentum/energy fluxes large enough to have a significant impact on galaxy formation.

preprint2011arXiv

Gyrokinetic Simulations of Solar Wind Turbulence from Ion to Electron Scales

The first three-dimensional, nonlinear gyrokinetic simulation of plasma turbulence resolving scales from the ion to electron gyroradius with a realistic mass ratio is presented, where all damping is provided by resolved physical mechanisms. The resulting energy spectra are quantitatively consistent with a magnetic power spectrum scaling of $k^{-2.8}$ as observed in \emph{in situ} spacecraft measurements of the "dissipation range" of solar wind turbulence. Despite the strongly nonlinear nature of the turbulence, the linear kinetic \Alfven wave mode quantitatively describes the polarization of the turbulent fluctuations. The collisional ion heating is measured at sub-ion-Larmor radius scales, which provides the first evidence of the ion entropy cascade in an electromagnetic turbulence simulation.

preprint2011arXiv

Incorporating Kinetic Physics into a Two-Fluid Solar-Wind Model with Temperature Anisotropy and Low-Frequency Alfven-Wave Turbulence

We develop a 1D solar-wind model that includes separate energy equations for the electrons and protons, proton temperature anisotropy, collisional and collisionless heat flux, and an analytical treatment of low-frequency, reflection-driven, Alfven-wave turbulence. To partition the turbulent heating between electron heating, parallel proton heating, and perpendicular proton heating, we employ results from the theories of linear wave damping and nonlinear stochastic heating. We account for mirror and oblique firehose instabilities by increasing the proton pitch-angle scattering rate when the proton temperature anisotropy exceeds the threshold for either instability. We numerically integrate the equations of the model forward in time until a steady state is reached, focusing on two fast-solar-wind-like solutions. These solutions are consistent with a number of observations, supporting the idea that Alfven-wave turbulence plays an important role in the origin of the solar wind.

preprint2011arXiv

On The GeV & TeV Detections of the Starburst Galaxies M82 & NGC 253

The GeV and TeV emission from M82 and NGC 253 observed by Fermi, HESS, and VERITAS constrains the physics of cosmic rays (CRs) in these dense starbursts. We argue that the gamma rays are predominantly hadronic in origin, as expected by previous studies. The measured fluxes imply that pionic losses are efficient for CR protons in both galaxies: we show that a fraction F_cal ~ 0.2 - 0.4 of the energy injected in high energy primary CR protons is lost to inelastic proton-proton collisions (pion production) before escape, producing gamma rays, neutrinos, and secondary electrons and positrons. We discuss the factor ~2 uncertainties in this estimate, including supernova rate and leptonic contributions to the GeV-TeV emission. We argue that gamma-ray data on ULIRGs like Arp 220 can test whether M82 and NGC 253 are truly calorimetric, and we present upper limits on Arp 220 from the Fermi data. We show that the observed ratio of the GeV to GHz fluxes of the starbursts suggests that non-synchrotron cooling processes are important for cooling the CR electron/positron population. We briefly reconsider previous predictions in light of the gamma-ray detections, including the starburst contribution to the gamma-ray background and CR energy densities. Finally, as a guide for future studies, we list the brightest star-forming galaxies on the sky and present updated predictions for their gamma-ray and neutrino fluxes.

preprint2011arXiv

PTF10iya: A short-lived, luminous flare from the nuclear region of a star-forming galaxy

We present the discovery and characterisation of PTF10iya, a short-lived (dt ~ 10 d, with an optical decay rate of ~ 0.3 mag per d), luminous (M_g ~ -21 mag) transient source found by the Palomar Transient Factory. The ultraviolet/optical spectral energy distribution is reasonably well fit by a blackbody with T ~ 1-2 x 10^4 K and peak bolometric luminosity L_BB ~ 1-5 x 10^44 erg per s (depending on the details of the extinction correction). A comparable amount of energy is radiated in the X-ray band that appears to result from a distinct physical process. The location of PTF10iya is consistent with the nucleus of a star-forming galaxy (z = 0.22405 +/- 0.00006) to within 350 mas (99.7 per cent confidence radius), or a projected distance of less than 1.2 kpc. At first glance, these properties appear reminiscent of the characteristic "big blue bump" seen in the near-ultraviolet spectra of many active galactic nuclei (AGNs). However, emission-line diagnostics of the host galaxy, along with a historical light curve extending back to 2007, show no evidence for AGN-like activity. We therefore consider whether the tidal disruption of a star by an otherwise quiescent supermassive black hole may account for our observations. Though with limited temporal information, PTF10iya appears broadly consistent with the predictions for the early "super-Eddington" phase of a solar-type star disrupted by a ~ 10^7 M_sun black hole. Regardless of the precise physical origin of the accreting material, the large luminosity and short duration suggest that otherwise quiescent galaxies can transition extremely rapidly to radiate near the Eddington limit; many such outbursts may have been missed by previous surveys lacking sufficient cadence.

preprint2011arXiv

Self-Regulated Star Formation in Galaxies via Momentum Input from Massive Stars

Feedback from massive stars is believed to play a critical role in shaping the galaxy mass function, the structure of the interstellar medium (ISM), and the low efficiency of star formation, but the exact form of the feedback is uncertain. In this paper, the first in a series, we present and test a novel numerical implementation of stellar feedback resulting from momentum imparted to the ISM by radiation, supernovae, and stellar winds. We employ a realistic cooling function, and find that a large fraction of the gas cools to <100K, so that the ISM becomes highly inhomogeneous. Despite this, our simulated galaxies reach an approximate steady state, in which gas gravitationally collapses to form giant molecular clouds (GMCs), dense clumps, and stars; subsequently, stellar feedback disperses the GMCs, repopulating the diffuse ISM. This collapse and dispersal cycle is seen in models of SMC-like dwarfs, the Milky-Way, and z~2 clumpy disk analogues. The simulated global star formation efficiencies are consistent with the observed Kennicutt-Schmidt relation. Moreover, the star formation rates are nearly independent of the numerically imposed high-density star formation efficiency, density threshold, and density scaling. This is a consequence of the fact that, in our simulations, star formation is regulated by stellar feedback limiting the amount of very dense gas available for forming stars. In contrast, in simulations without stellar feedback, i.e. under the action of only gravity and gravitationally-induced turbulence, the ISM experiences runaway collapse to very high densities. In these simulations without feedback, the global star formation rates exceed observed galactic star formation rates by 1-2 orders of magnitude, demonstrating that stellar feedback is crucial to the regulation of star formation in galaxies.

preprint2011arXiv

Spectroscopic Signatures of the Tidal Disruption of Stars by Massive Black Holes

During the tidal disruption of a star by a massive black hole (BH) of mass MBH <~ 10^7 Msun, stellar debris falls back to the BH at a rate well above the Eddington rate. A fraction of this gas is subsequently blown away from the BH, producing an optically bright flare of radiation. We predict the spectra and spectral evolution of tidal disruption events, focusing on the photoionized gas outside this outflow's photosphere. The spectrum will show absorption lines that are strongly blueshifted relative to the host galaxy, very broad (0.01-0.1c), and strongest at UV wavelengths (e.g., C IV, Ly alpha, O VI), lasting ~ 1 month for a 10^6 Msun BH. Meanwhile, supernovae in galactic nuclei are a significant source of confusion in optical surveys for tidal disruption events: we estimate that nuclear Type Ia supernovae are two orders of magnitude more common than tidal disruption events at z ~ 0.1 for ground-based surveys. Nuclear Type II supernovae occur at a comparable rate but can be excluded by pre-selecting red galaxies. Supernova contamination can be reduced to a manageable level using high-resolution follow-up imaging with adaptive optics or the Hubble Space Telescope. Our predictions should help optical transient surveys capitalize on their potential for discovering tidal disruption events.

preprint2011arXiv

Swift 1644+57: The Longest Gamma-ray Burst?

Swift recently discovered an unusual gamma-ray and x-ray transient (Sw 1644+57) that was initially identified as a long-duration gamma-ray burst (GRB). However, the ~ 10 keV x-ray emission has persisted for over a month with a luminosity comparable to its peak value. The astrometric coincidence of the source with the center of its host galaxy, together with other considerations, motivated the interpretation that Sw 1644+57 was produced by an outburst from a 10^{6-7} M_sun black hole at the center of the galaxy. Here we consider the alternate possibility that Sw 1644+57 is indeed a long-duration GRB, albeit a particularly long one. We discuss the general properties of very long-duration, low-power GRB-like transients associated with the core-collapse of a massive star. Both neutron star (magnetar) spindown and black hole accretion can power such events. The requirements for producing low-power, very long-duration GRBs by magnetar spindown are similar to those for powering extremely luminous supernovae by magnetar spindown, suggesting a possible connection between these two unusual types of transients. Alternatively, Sw 1644+57 be associated with the faintest core-collapse explosions: the collapse of a rotating red supergiant in a nominally failed supernova can power accretion onto a solar-mass black hole for up to ~ 100 days; the jet produced by black hole accretion inevitably unbinds the outer envelope of the progenitor, leading to a weak ~ 10^{49} erg explosion. In both neutron star and black hole models, a jet can burrow through the host star in a few days, with a kinetic luminosity ~ 10^{45-46} ergs/s, sufficient to power the observed emission of Sw 1644+57.

preprint2011arXiv

The Radial Velocity Signature of Tides Raised in Stars Hosting Exoplanets

Close-in, massive exoplanets raise significant tides in their stellar hosts. We compute the radial velocity (RV) signal due to this fluid motion in the equilibrium tide approximation. The predicted radial velocities in the observed sample of exoplanets exceed 1 m/s for 17 systems, with the largest predicted signal being ~30 m/s for WASP-18 b. Tidally-induced RV's are thus detectable with present methods. Both tidal fluid flow and the epicyclic motion of a slightly eccentric orbit produce an RV signal at twice the orbital frequency. If care is not taken, the tidally induced RV may, in some cases, be confused with a finite orbital eccentricity. Indeed, WASP-18 b is reported to have an eccentric orbit with small e=0.009 and pericenter longitude ω=-π/2. Whereas such a close alignment of the orbit and line of sight to the observer requires fine tuning, this phase in the RV signal is naturally explained by the tidal velocity signature of an e=0 orbit. Additionally, the equilibrium tide estimate for the amplitude is in rough agreement with the data. Thus the reported eccentricity for WASP-18 b is instead likely a signature of the tidally-induced RV in the stellar host. Measurement of both the orbital and tidal velocity for non-transiting planets may allow planet mass and inclination to be separately determined solely from radial velocity data. We suggest that high precision fitting of RV data should include the tidal velocity signal in those cases where it may affect the determination of orbital parameters.

preprint2011arXiv

Thermal Instability & the Feedback Regulation of Hot Halos in Clusters, Groups, and Galaxies

Observations of clusters and groups imply that such halos are roughly in global thermal equilibrium, with heating balancing cooling when averaged over sufficiently long time- and length-scales; the ICM is, however, very likely to be locally thermally unstable. Using simple observationally-motivated heating prescriptions, we show that local thermal instability (TI) can produce a multi-phase medium---with ~ 10000 K cold filaments condensing out of the hot ICM---only when the ratio of the TI timescale in the hot plasma (t_{TI}) to the free-fall timescale (t_{ff}) satisfies t_{TI}/t_{ff} <~ 10. This criterion quantitatively explains why cold gas and star formation are preferentially observed in low-entropy clusters and groups. In addition, the interplay among heating, cooling, and TI reduces the net cooling rate and the mass accretion rate at small radii by factors of ~ 100 relative to cooling-flow models. This dramatic reduction is in line with observations. The feedback efficiency required to prevent a cooling-flow is ~ 0.001 for clusters and decreases for lower mass halos; supernova heating may be energetically sufficient to balance cooling in galactic halos. We further argue that the ICM self-adjusts so that t_{TI}/t_{ff} >~ 10 at all radii. When this criterion is not satisfied, cold filaments condense out of the hot phase and reduce the density of the ICM. These cold filaments can power the black hole and/or stellar feedback required for global thermal balance, which drives t_{TI}/t_{ff} >~ 10. In comparison to clusters, groups have central cores with lower densities and larger radii. This can account for the deviations from self-similarity in the X-ray luminosity-temperature (L_X-T_X) relation. The high-velocity clouds observed in the Galactic halo can also be due to local TI producing multi-phase gas close to the virial radius.

preprint2011arXiv

Tidal asteroseismology: Kepler's KOI-54

We develop a general framework for interpreting and analyzing high-precision lightcurves from eccentric stellar binaries. Although our methods are general, we focus on the recently discovered Kepler system KOI-54, a face-on binary of two A stars with $e=0.83$ and an orbital period of 42 days. KOI-54 exhibits strong ellipsoidal variability during its periastron passage; its lightcurve also contains ~20 pulsations at perfect harmonics of the orbital frequency, and another ~10 nonharmonic pulsations. Analysis of such data is a new form of asteroseismology in which oscillation amplitudes and phases rather than frequencies contain information that can be mined to constrain stellar properties. We qualitatively explain the physics of mode excitation and the range of harmonics expected to be observed. To quantitatively model observed pulsation spectra, we develop and apply a linear, tidally forced, nonadiabatic stellar oscillation formalism including the Coriolis force. We produce temporal power spectra for KOI-54 that are semi-quantitatively consistent with the observations. Both stars in the KOI-54 system are expected to be rotating pseudosynchronously, with resonant nonaxisymmetric modes providing a key contribution to the total torque; such resonances provide a possible explanation for the two largest-amplitude harmonic pulsations observed in KOI-54, although we find quantitative problems with this interpretation. We show in detail that the nonharmonic pulsations observed in KOI-54 can be produced by nonlinear three-mode coupling. The methods developed in this paper can be generalized in the future to determine the best-fit stellar parameters given pulsation data. We also derive an analytic model of KOI-54's ellipsoidal variability, including both tidal distortion and stellar irradiation, which can be used to model other similar systems.

preprint2011arXiv

Turbulent Pressure Support in the Outer Parts of Galaxy Clusters

We use three-dimensional MHD simulations with anisotropic thermal conduction to study turbulence due to the magnetothermal instability (MTI) in the intracluster medium (ICM) of galaxy clusters. The MTI grows on timescales of ~1 Gyr and is capable of driving vigorous, sustained turbulence in the outer parts of galaxy clusters if the temperature gradient is maintained in spite of the rapid thermal conduction. If this is the case, turbulence due to the MTI can provide up to 5-30% of the pressure support beyond r_500 in galaxy clusters, an effect that is strongest for hot, massive clusters. The turbulence driven by the MTI is generally additive to other sources of turbulence in the ICM, such as that produced by structure formation. This new source of non-thermal pressure support reduces the observed Sunyaev-Zel'dovich (SZ) signal and X-ray pressure gradient for a given cluster mass and introduces a cluster mass and temperature gradient-dependent bias in SZ and X-ray mass estimates of clusters. This additional physics may also need to be taken into account when estimating the matter power spectrum normalization, sigma-8, through simulation templates from the observed amplitude of the SZ power spectrum.

preprint2010arXiv

Can conduction induce convection? The non-linear saturation of buoyancy instabilities in dilute plasmas

We study the effects of anisotropic thermal conduction on low-collisionality, astrophysical plasmas using two and three-dimensional magnetohydrodynamic simulations. For weak magnetic fields, dilute plasmas are buoyantly unstable for either sign of the temperature gradient: the heat-flux-driven buoyancy instability (HBI) operates when the temperature increases with radius while the magnetothermal instability (MTI) operates in the opposite limit. In contrast to previous results, we show that, in the presence of a sustained temperature gradient, the MTI drives strong turbulence and operates as an efficient magnetic dynamo (akin to standard, adiabatic convection). Together, the turbulent and magnetic energies contribute up to ~10% of the pressure support in the plasma. In addition, the MTI drives a large convective heat flux, ~1.5% of rho c_s^3. These findings are robust even in the presence of an external source of strong turbulence. Our results on the nonlinear saturation of the HBI are consistent with previous studies but we explain physically why the HBI saturates quiescently by re-orienting the magnetic field (suppressing the conductive heat flux through the plasma), while the MTI saturates by generating sustained turbulence. We also systematically study how an external source of turbulence affects the saturation of the HBI: such turbulence can disrupt the HBI only on scales where the shearing rate of the turbulence is faster than the growth rate of the HBI. In particular, our results provide a simple mapping between the level of turbulence in a plasma and the effective isotropic thermal conductivity. We discuss the astrophysical implications of these findings, with a particular focus on the intracluster medium of galaxy clusters.

preprint2010arXiv

How Do Massive Black Holes Get Their Gas?

We use multi-scale SPH simulations to follow the inflow of gas from galactic scales to <0.1pc, where the gas begins to resemble a traditional Keplerian accretion disk. The key ingredients are gas, stars, black holes (BHs), self-gravity, star formation, and stellar feedback. We use ~100 simulations to survey a large parameter space of galaxy properties and subgrid models for the ISM physics. We generate initial conditions for our simulations of galactic nuclei (<~300pc) using galaxy scale simulations, including both major mergers and isolated bar-(un)stable disk galaxies. For sufficiently gas-rich, disk-dominated systems, a series of gravitational instabilities generates large accretion rates of up to 1-10 M_sun/yr onto the BH (at <<0.1pc); sufficient to fuel the most luminous quasars. The BH accretion rate is highly time variable, given fixed conditions at ~kpc. At >~10pc, our simulations resemble the 'bars within bars' model, but the gas exhibits diverse morphologies, including spirals, rings, clumps, and bars; their duty cycle is modest, complicating attempts to correlate BH accretion with nuclear morphology. At ~1-10pc, the gravitational potential becomes dominated by the BH and bar-like modes are no longer present. However, the gas becomes unstable to a standing, eccentric disk or a single-armed spiral mode (m=1), driving the gas to sub-pc scales. Proper treatment of this mode requires including star formation and the self-gravity of both the stars and gas. We predict correlations between BHAR and SFR at different galactic nuclei: nuclear SF is more tightly coupled to AGN activity, but correlations exist at all scales.

preprint2010arXiv

Perpendicular Ion Heating by Low-Frequency Alfven-Wave Turbulence in the Solar Wind

We consider ion heating by turbulent Alfven waves (AWs) and kinetic Alfven waves (KAWs) with perpendicular wavelengths comparable to the ion gyroradius and frequencies smaller than the ion cyclotron frequency. When the turbulence amplitude exceeds a certain threshold, an ion's orbit becomes chaotic. The ion then interacts stochastically with the time-varying electrostatic potential, and the ion's energy undergoes a random walk. Using phenomenological arguments, we derive an analytic expression for the rates at which different ion species are heated, which we test by simulating test particles interacting with a spectrum of randomly phased AWs and KAWs. We find that the stochastic heating rate depends sensitively on the quantity epsilon = dv/vperp, where vperp is the component of the ion velocity perpendicular to the background magnetic field B0, and dv (dB) is the rms amplitude of the velocity (magnetic-field) fluctuations at the gyroradius scale. In the case of thermal protons, when epsilon << eps1, where eps1 is a constant, a proton's magnetic moment is nearly conserved and stochastic heating is extremely weak. However, when epsilon > eps1, the proton heating rate exceeds the cascade power that would be present in strong balanced KAW turbulence with the same value of dv, and magnetic-moment conservation is violated. For the random-phase waves in our test-particle simulations, eps1 is approximately 0.2. For protons in low-beta plasmas, epsilon is approximately dB/B0 divided by the square root of beta, and epsilon can exceed eps1 even when dB/B0 << eps1. At comparable temperatures, alpha particles and minor ions have larger values of epsilon than protons and are heated more efficiently as a result. We discuss the implications of our results for ion heating in coronal holes and the solar wind.

preprint2010arXiv

Self-Regulated Black Hole Growth via Momentum Deposition in Galaxy Merger Simulations

We perform hydrodynamical simulations of major galaxy mergers using new methods for calculating the growth of massive black holes (BH) in galactic nuclei and their impact on the surrounding galaxy. We model BH growth by including a subgrid model for accretion produced by angular momentum transport on unresolved scales. The impact of the BHs radiation on surrounding gas is approximated by depositing momentum into the ambient gas, which produces an outward force away from the BH. We argue that these phenomenological models for BH growth and feedback better approximate the interaction between the BH and dense gas in galaxies than previous models. We show that this physics leads to self-regulated black hole growth: during the peak of activity, the accretion rate onto the BH is largely determined by the physics of BH feedback, not the subgrid accretion model. The BH significantly modifies the gas dynamics in the galactic nucleus (< 300 pc), but does not generate large-scale galactic outflows. Integrated over an entire galaxy merger, BH feedback has little effect on the total number of stars formed, but is crucial for setting the BHs mass.

preprint2010arXiv

The Growth of Massive Black Holes in Galaxy Merger Simulations with Feedback by Radiation Pressure

We study the growth of massive black holes (BH) in galaxies using smoothed particle hydrodynamic simulations of major galaxy mergers with new implementations of BH accretion and feedback. The effect of BH accretion on gas in its host galaxy is modeled by depositing momentum at a rate ~ tau L/c into the ambient gas, where L is the luminosity produced by accretion onto the BH and tau is the wavelength-averaged optical depth of the galactic nucleus to the AGN's radiation (a free parameter of our model). The accretion rate onto the BH is relatively independent of our subgrid accretion model and is instead determined by the BH's dynamical impact on its host galaxy: BH accretion is thus self-regulated rather than `supply limited.' We show that the final BH mass and total stellar mass formed during a merger are more robust predictions of the simulations than the time dependence of the star formation rate or BH accretion rate. In particular, the latter depend on the assumed interstellar medium physics, which determines when and where the gas fragments to form star clusters; this in turn affects the fuel available for further star formation and BH growth. Simulations over a factor of ~ 30 in galaxy mass are consistent with the observed M_BH-sigma relation for a mean optical depth of tau ~ 25. This requires that most BH growth occur when the galactic nucleus is optically thick to far-infrared radiation, consistent with the hypothesized connection between ultra-luminous infrared galaxies and quasars. We find tentative evidence for a shallower M_BH-sigma relation in the lowest mass galaxies, sigma < 100 km/s. Our results demonstrate that feedback-regulated BH growth and consistency with the observed M_BH-sigma relation do not require that BH feedback terminate star formation in massive galaxies or unbind large quantities of cold gas.

preprint2010arXiv

The Nuclear Stellar Disk in Andromeda: A Fossil from the Era of Black Hole Growth

The physics of angular momentum transport from galactic scales (~10-100 pc) to much smaller radii is one of the oustanding problems in our understanding of the formation and evolution of super-massive black holes (BHs). Seemingly unrelated observations have discovered that there is a lopsided stellar disk of unknown origin orbiting the BH in M31, and possibly many other systems. We show that these nominally independent puzzles are in fact closely related. Multi-scale simulations of gas inflow from galactic to BH scales show that when sufficient gas is driven towards a BH, gravitational instabilities form a lopsided, eccentric disk that propagates inwards from larger radii. The lopsided stellar disk exerts a strong torque on the remaining gas, driving inflows that fuel the growth of the BH and produce quasar-level luminosities. The same disk can produce significant obscuration along many sightlines and thus may be the putative 'torus' invoked to explain obscured active galactic nuclei and the cosmic X-ray background. The stellar relic of this disk is long lived and retains the eccentric pattern. Simulations that yield quasar-level accretion rates produce relic stellar disks with kinematics, eccentric patterns, precession rates, and surface density profiles in reasonable agreement with observations of M31. The observed properties of nuclear stellar disks can thus be used to constrain the formation history of super-massive BHs.

preprint2010arXiv

The Physics of the FIR-Radio Correlation: I. Calorimetry, Conspiracy, and Implications

(Abridged) The far-infrared (FIR) and radio luminosities of star-forming galaxies are linearly correlated over a very wide range in star formation rate, from normal spirals like the Milky Way to the most intense starbursts. Using one-zone models of cosmic ray (CR) injection, cooling, and escape in star-forming galaxies, we attempt to reproduce the observed FIR-radio correlation over its entire span. We show that ~2% of the kinetic energy from supernova explosions must go into primary CR electrons and that ~10 - 20% must go into primary CR protons. Secondary electrons and positrons are likely comparable to or dominate primary electrons in dense starburst galaxies. We discuss the implications of our models for the magnetic field strengths of starbursts, the detectability of starbursts by Fermi, and cosmic ray feedback. Overall, our models indicate that both CR protons and electrons escape from low surface density galaxies, but lose most of their energy before escaping dense starbursts. The FIR-radio correlation is caused by a combination of the efficient cooling of CR electrons (calorimetry) in starbursts and a conspiracy of several factors. For lower surface density galaxies, the decreasing radio emission caused by CR escape is balanced by the decreasing FIR emission caused by the low effective UV dust opacity. In starbursts, bremsstrahlung, ionization, and Inverse Compton cooling decrease the radio emission, but they are countered by secondary electrons/positrons and the decreasing critical synchrotron frequency, which both increase the radio emission. Our conclusions hold for a broad range of variations on our fiducial model.

preprint2010arXiv

Thermal Instability with Anisotropic Thermal Conduction and Adiabatic Cosmic Rays: Implications for Cold Filaments in Galaxy Clusters

Observations of the cores of nearby galaxy clusters show H$α$ and molecular emission line filaments. We argue that these are the result of {\em local} thermal instability in a {\em globally} stable galaxy cluster core. We present local, high resolution, two-dimensional magnetohydrodynamic simulations of thermal instability for conditions appropriate to the intracluster medium (ICM); the simulations include thermal conduction along magnetic field lines and adiabatic cosmic rays. Thermal conduction suppresses thermal instability along magnetic field lines on scales smaller than the Field length ($\gtrsim$10 kpc for the hot, diffuse ICM). We show that the Field length in the cold medium must be resolved both along and perpendicular to the magnetic field in order to obtain numerically converged results. Because of negligible conduction perpendicular to the magnetic field, thermal instability leads to fine scale structure in the perpendicular direction. Filaments of cold gas along magnetic field lines are thus a natural consequence of thermal instability with anisotropic thermal conduction. Nonlinearly, filaments of cold ($\sim 10^4$ K) gas should have lengths (along the magnetic field) comparable to the Field length in the cold medium $\sim 10^{-4}$ pc! Observations show, however, that the atomic filaments in clusters are far more extended, $\sim 10$ kpc. Cosmic ray pressure support (or a small scale turbulent magnetic pressure) may resolve this discrepancy: even a small cosmic ray pressure in the diffuse ICM, $\sim 10^{-4}$ of the thermal pressure, can be adiabatically compressed to provide significant pressure support in cold filaments. This is qualitatively consistent with the large population of cosmic rays invoked to explain the atomic and molecular line ratios observed in filaments.

preprint2010arXiv

Time Dependent Models of Flares from Sagittarius A*

The emission from Sgr A*, the supermassive black hole in the Galactic Center, shows order of magnitude variability ("flares") a few times a day that is particularly prominent in the near-infrared (NIR) and X-rays. We present a time-dependent model for these flares motivated by the hypothesis that dissipation of magnetic energy powers the flares. We show that episodic magnetic reconnection can occur near the last stable circular orbit in time-dependent magnetohydrodynamic simulations of black hole accretion - the timescales and energetics of these events are broadly consistent with the flares from Sgr A*. Motivated by these results, we present a spatially one-zone time-dependent model for the electron distribution function in flares, including energy loss due to synchrotron cooling and adiabatic expansion. Synchrotron emission from transiently accelerated particles can explain the NIR/X-ray lightcurves and spectra of a luminous flare observed 4 April 2007. A significant decrease in the magnetic field strength during the flare (coincident with the electron acceleration) is required to explain the simultaneity and symmetry of the simultaneous lightcurves. Our models predict that the NIR and X-ray spectral indices differ by 0.5 and that there is only modest variation in the spectral index during flares. We also explore implications of this model for longer wavelength (radio-submm) emission seemingly associated with X-ray and NIR flares; we argue that a few hour decrease in the submm emission is a more generic consequence of large-scale magnetic reconnection than delayed radio emission from adiabatic expansion.

preprint2010arXiv

Turbulence in Galaxy Cluster Cores: a Key to Cluster Bimodality?

We study the effects of externally imposed turbulence on the thermal properties of galaxy cluster cores, using three-dimensional numerical simulations including magnetic fields, anisotropic thermal conduction, and radiative cooling. The imposed "stirring" crudely approximates the effects of galactic wakes, waves generated by galaxies moving through the intracluster medium (ICM), and/or turbulence produced by a central active galactic nucleus. The simulated clusters exhibit a strong bimodality. Modest levels of turbulence, ~100 km/s (~10% of the sound speed), suppress the heat-flux-driven buoyancy instability (HBI), resulting in an isotropically tangled magnetic field and a quasi-stable, high entropy, thermal equilibrium with no cooling catastrophe. Thermal conduction dominates the heating of the cluster core, but turbulent mixing is critical because it suppresses the HBI and (to a lesser extent) the thermal instability. Lower levels of turbulent mixing (approximately less than 100 km/s) are insufficient to suppress the HBI, rapidly leading to a thermal runaway and a cool-core cluster. Remarkably, then, small fluctuations in the level of turbulence in galaxy cluster cores can initiate transitions between cool-core (low entropy) and non cool-core (high entropy) states.

preprint2009arXiv

A Maximum Stellar Surface Density in Dense Stellar Systems

We compile observations of the surface mass density profiles of dense stellar systems, including globular clusters in the Milky Way and nearby galaxies, massive star clusters in nearby starbursts, nuclear star clusters in dwarf spheroidals and late-type disks, ultra-compact dwarfs, and galaxy spheroids spanning the range from low-mass cusp bulges and ellipticals to massive core ellipticals. We show that in all cases the maximum stellar surface density attained in the central regions of these systems is similar, Sigma_max ~ 10^11 M_sun/kpc^2 (~20 g/cm^2), despite the fact that the systems span 7 orders of magnitude in total stellar mass M_star, 5 in effective radius R_e, and have a wide range in effective surface density M_star/R_e^2. The surface density limit is reached on a wide variety of physical scales in different systems and is thus not a limit on three-dimensional stellar density. Given the very different formation mechanisms involved in these different classes of objects, we argue that a single piece of physics likely determines Sigma_max. The radiation fields and winds produced by massive stars can have a significant influence on the formation of both star clusters and galaxies, while neither supernovae nor black hole accretion are important in star cluster formation. We thus conclude that feedback from massive stars likely accounts for the observed Sigma_max, plausibly because star formation reaches an Eddington-like flux that regulates the growth of these diverse systems. This suggests that current models of galaxy formation, which focus on feedback from supernovae and active galactic nuclei, are missing a crucial ingredient.

preprint2009arXiv

Are Most Low-Luminosity AGN Really Obscured?

At low Eddington ratio (mdot), two effects make it harder to detect AGN given some selection criteria. First, even with fixed accretion physics, AGN are diluted/less luminous relative to their hosts; the magnitude of this depends on host properties and so on luminosity and redshift. Second, they may transition to a radiatively inefficient state, changing SED shape and dramatically decreasing in optical/IR luminosity. These effects lead to differences in observed AGN samples, even at fixed bolometric luminosity and after correction for obscuration. The true Eddington ratio distribution may depend strongly on luminosity, but this will be seen only in surveys robust to dilution and radiative inefficiency (X-ray or narrow-line samples); selection effects imply that AGN in optical samples will have uniformly high mdot. This also implies that different selection methods yield systems with different hosts: the clustering of faint optical/IR sources will be weaker than that of X-ray sources, and optical/IR Seyferts will reside in more disk-dominated galaxies while X-ray selected Seyferts will preferentially occupy early-type systems. If observed mdot distributions are correct, a large fraction of low-luminosity AGN currently classified as 'obscured' are in fact diluted and/or radiatively inefficient, not obscured by gas or dust. This is equally true if X-ray hardness is used as a proxy for obscuration, since radiatively inefficient SEDs near mdot~0.01 are X-ray hard. These effects can explain most of the claimed luminosity/redshift dependence in the 'obscured' AGN population, with the true obscured fraction as low as 20%.

preprint2009arXiv

Compact High-Redshift Galaxies Are the Cores of the Most Massive Present-Day Spheroids

Observations suggest that effective radii of high-z massive spheroids are as much as a factor ~6 smaller than low-z galaxies of comparable mass. Given the apparent absence of low-z counterparts, this has often been interpreted as indicating that the high density, compact red galaxies must be 'puffed up' by some mechanism. We compare the ensemble of high-z observations with large samples of well-observed low-z ellipticals. At the same physical radii, the stellar surface mass densities of low and high-z systems are comparable. Moreover, the abundance of high surface density material at low redshift is comparable to or larger than that observed at z>1-2, consistent with the continuous buildup of spheroids over this time. The entire population of compact, high-z red galaxies may be the progenitors of the high-density cores of present-day ellipticals, with no need for a decrease in stellar density from z=2 to z=0. The primary difference between low and high-z systems is thus the observed low-density material at large radii in low-z spheroids (rather than the high-density material in high-z spheroids). Such low-density material may either (1) assemble at z<2 or (2) be present, but not yet detected, at z>2. Mock observations of low-z massive systems show that the high-z observations do not yet probe sufficiently low surface brightness material to detect the low surface density 'wings' (if present). Thus, if the high-z galaxies resemble the most massive systems today, their inferred effective radii could be under-estimated by factors ~2-4. This difference arises because massive systems at low redshift are not well-fit by single Sersic profiles. We discuss implications of our results for physical models of galaxy evolution.

preprint2009arXiv

Constraining Low-Frequency Alfvenic Turbulence in the Solar Wind Using Density Fluctuation Measurements

One proposed mechanism for heating the solar wind, from close to the sun to beyond 10 AU, invokes low-frequency, oblique, Alfven-wave turbulence. Because small-scale oblique Alfven waves (kinetic Alfven waves) are compressive, the measured density fluctuations in the solar wind place an upper limit on the amplitude of kinetic Alfven waves and hence an upper limit on the rate at which the solar wind can be heated by low-frequency, Alfvenic turbulence. We evaluate this upper limit for both coronal holes at 5 solar radii and in the near-Earth solar wind. At both radii, the upper limit we find is consistent with models in which the solar wind is heated by low-frequency Alfvenic turbulence. At 1 AU, the upper limit on the turbulent heating rate derived from the measured density fluctuations is within a factor of 2 of the measured solar wind heating rate. Thus if low-frequency Alfvenic turbulence contributes to heating the near-Earth solar wind, kinetic Alfven waves must be one of the dominant sources of solar wind density fluctuations at frequencies of order 1 Hz. We also present a simple argument for why density fluctuation measurements do appear to rule out models in which the solar wind is heated by non-turbulent high-frequency waves ``sweeping'' through the ion-cyclotron resonance, but are compatible with heating by low-frequency Alfvenic turbulence.

preprint2009arXiv

On the Interpretation of Magnetic Helicity Signatures in the Dissipation Range of Solar Wind Turbulence

Measurements of small-scale turbulent fluctuations in the solar wind find a non-zero right-handed magnetic helicity. This has been interpreted as evidence for ion cyclotron damping. However, theoretical and empirical evidence suggests that the majority of the energy in solar wind turbulence resides in low frequency anisotropic kinetic Alfven wave fluctuations that are not subject to ion cyclotron damping. We demonstrate that a dissipation range comprised of kinetic Alfven waves also produces a net right-handed fluctuating magnetic helicity signature consistent with observations. Thus, the observed magnetic helicity signature does not necessarily imply that ion cyclotron damping is energetically important in the solar wind.

preprint2009arXiv

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, 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 in the optical and near-IR. Some 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 strongest for BHs with MBH ~ 10^5 - 10^6 Msun 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 observational challenges in the optical. Pan-STARRS should detect at least several events per year--many more if current theoretical models of super-Eddington outflows are correct. 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.

preprint2009arXiv

Radio Emission from Supernova Remnants: Implications for Post-Shock Magnetic Field Amplification and the Magnetic Fields of Galaxies

Using observations from the literature, we show that the non-thermal radio luminosity (L) of supernova remnants (SNRs) is a strong function of the average gas surface density (Sigma) of the galaxy in which the remnants reside, from normal spirals to luminous starbursts. We combine a simple theory for electron cooling in SNRs with the observed radio luminosities to estimate the remnant magnetic field strength (B_SNR): the correlation between L and Sigma implies that B_SNR also increases with Sigma. We explore two interpretations of this correlation: (1) B_SNR is generated by post-shock magnetic field amplification, with B_SNR^2 proportional to Sigma and (2) B_SNR results from shock-compression of the ambient ISM magnetic field (B_ISM), with B_ISM being larger in denser galaxies. We find that shock compression is, on average, sufficient to produce the observed radio emission from SNRs in the densest starbursts; amplification of post-shock magnetic fields is not required. By contrast, in normal spirals post-shock field amplification (by a factor of a few - 10) is consistent with the data; we find tentative evidence that both the Alfven speed and the ratio of B_SNR^2 to the post-shock pressure ("epsilon_B") are constant in SNRs from galaxy to galaxy. We discuss observational tests that can be used to distinguish between these two interpretations of the radio luminosities of SNRs. Regardless of which is correct, the radio emission from SNRs provides an upper limit to B_ISM that is independent of the minimum energy assumption. For the densest starbursts, the ISM magnetic energy density is below that required for hydrostatic equilibrium; thus magnetic fields are not dynamically important on the largest scales in starbursts, in contrast with spiral galaxies like our own. This dichotomy may have implications for galactic dynamo theory.

preprint2009arXiv

The Disruption of Giant Molecular Clouds by Radiation Pressure and the Efficiency of Star Formation in Galaxies

Star formation is slow, in the sense that the gas consumption time is much longer than the dynamical time. It is also inefficient; essentially all star formation in local galaxies takes place in giant molecular clouds (GMCs), but the fraction of a GMC converted to stars is very small, ~5%. In the most luminous starbursts, the GMC lifetime is shorter than the main sequence lifetime of even the most massive stars, so that supernovae can play no role in GMC disruption. We investigate the disruption of GMCs across a wide range of galaxies, from normal spirals to the densest starbursts; we take into account the effects of HII gas pressure, shocked stellar winds, protostellar jets, and radiation pressure produced by the absorption and scattering of starlight on dust grains. In the Milky Way, we find that a combination of three mechanisms, jets, HII gas pressure, and radiation pressure, disrupts the clouds. In more rapidly star forming galaxies such as ``clump'' galaxies at high-redshift, ultra-luminous infrared galaxies (ULIRGs) and submillimeter galaxies, radiation pressure dominates natal cloud distribution. We predict the presence of 10-20 clusters with masses ~10^7 Msun in local ULIRGs such as Arp 220 and a similar number of clusters with M_* ~ 10^8 Msun in high redshift clump galaxies; submillimeter galaxies will have even more massive clusters. We find that the mass fraction of a GMC that ends up in stars is an increasing function of the gas surface density of a galaxy, reaching ~35% in the most luminous starbursts. Furthermore, the disruption of bubbles by radiation pressure stirs the interstellar medium to velocities of ~10 km/s in normal galaxies and to ~100 km/s in ULIRGs like Arp 220, consistent with observations. Thus, radiation pressure may play a dominant role in the ISM of star-forming galaxies.

preprint2009arXiv

The Heating of Test Particles in Numerical Simulations of Alfvenic Turbulence

We study the heating of charged test particles in three-dimensional numerical simulations of weakly compressible magnetohydrodynamic (MHD) turbulence (``Alfvenic turbulence''); these results are relevant to particle heating and acceleration in the solar wind, solar flares, accretion disks onto black holes, and other astrophysics and heliospheric environments. The physics of particle heating depends on whether the gyrofrequency of a particle is comparable to the frequency of a turbulent fluctuation that is resolved on the computational domain. Particles with these frequencies nearly equal undergo strong perpendicular heating (relative to the local magnetic field) and pitch angle scattering. By contrast, particles with large gyrofrequency undergo strong parallel heating. Simulations with a finite resistivity produce additional parallel heating due to parallel electric fields in small-scale current sheets. Many of our results are consistent with linear theory predictions for the particle heating produced by the Alfven and slow magnetosonic waves that make up Alfvenic turbulence. However, in contrast to linear theory predictions, energy exchange is not dominated by discrete resonances between particles and waves; instead, the resonances are substantially ``broadened.'' We discuss the implications of our results for solar and astrophysics problems, in particular the thermodynamics of the near-Earth solar wind. We conclude that Alfvenic turbulence produces significant parallel heating via the interaction between particles and magnetic field compressions (``slow waves''). However, on scales above the proton Larmor radius, Alfvenic turbulence does not produce significant perpendicular heating of protons or minor ions.

preprint2009arXiv

When Should We Treat Galaxies as Isolated?

Traditionally, secular evolution is defined as evolution of systems where the internal growth of structure and instabilities dominates the growth via external drivers (e.g. accretion/mergers). Most study has focused on 'isolated' galaxies, where seed asymmetries may represent realistic cosmological substructure, but subsequent evolution ignores galaxy growth. Large-scale modes in the disk then grow on a timescale of order a disk rotation period (0.1-1 Gyr). If, however, galaxies evolve cosmologically on a shorter timescale, then it may not be appropriate to consider them 'isolated.' We outline simple scalings to ask whether the timescale for secular evolution is shorter than the timescale for cosmological accretion and mergers. This is the case in a narrow, but important range of perturbation amplitudes corresponding to substructure or mode/bar fractional amplitudes 0.01-0.1, a range of interest for observed strong bars and pseudobulges. At smaller amplitudes <<0.1, systems are not isolated: typical disks will grow by accretion at a comparable level over even a single dynamical time. At larger amplitudes >>0.1, the evolution is no longer secular; direct gravitational evolution of the seed swamps the internal disk response. We derive criteria for when disks can be well-approximated as 'isolated' as a function of mass, redshift, and disk stability. The relevant parameter space shrinks at higher mass, higher disk stability, and higher-z as accretion rates increase. Cosmological rates of galaxy evolution also define a maximum bar/mode lifetime of practical interest, of ~0.1/H(z). Longer-lived modes will de-couple from their drivers (if driven) and encounter cosmological effects.

preprint2008arXiv

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 collisional 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 Larmor 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 cyclotron frequency due to the strong anisotropy of the turbulent fluctuations, k_parallel << k_perp (implied by critical balance). In this limit, the turbulence is optimally described by gyrokinetics; 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. Cascade 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 dissipation range spectrum is predicted to be an exponential fall off; the power-law behavior 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.

preprint2008arXiv

Astrometry with the Keck-Interferometer: the ASTRA project and its science

The sensitivity and astrometry upgrade ASTRA of the Keck Interferometer is introduced. After a brief overview of the underlying interferometric principles, the technology and concepts of the upgrade are presented. The interferometric dual-field technology of ASTRA will provide the KI with the means to observe two objects simultaneously, and measure the distance between them with a precision eventually better than 100 uas. This astrometric functionality of ASTRA will add a unique observing tool to fields of astrophysical research as diverse as exo-planetary kinematics, binary astrometry, and the investigation of stars accelerated by the massive black hole in the center of the Milky Way as discussed in this contribution.

preprint2008arXiv

Spherical Accretion with Anisotropic Thermal Conduction

We study the effects of anisotropic thermal conduction on magnetized spherical accretion flows using global axisymmetric MHD simulations. In low collisionality plasmas, the Bondi spherical accretion solution is unstable to the magnetothermal instability (MTI). The MTI grows rapidly at large radii where the inflow is subsonic. For a weak initial field, the MTI saturates by creating a primarily radial magnetic field, i.e., by aligning the field lines with the background temperature gradient. The saturation is quasilinear in the sense that the magnetic field is amplified by a factor of $\sim 10-30$ independent of the initial field strength (for weak fields). In the saturated state, the conductive heat flux is much larger than the convective heat flux, and is comparable to the field-free (Spitzer) value (since the field lines are largely radial). The MTI by itself does not appreciably change the accretion rate $\dot M$ relative to the Bondi rate $\dot M_B$. However, the radial field lines created by the MTI are amplified by flux freezing as the plasma flows in to small radii. Oppositely directed field lines are brought together by the converging inflow, leading to significant resistive heating. When the magnetic energy density is comparable to the gravitational potential energy density, the plasma is heated to roughly the virial temperature; the mean inflow is highly subsonic; most of the energy released by accretion is transported to large radii by thermal conduction; and the accretion rate $\dot M \ll \dot M_B$. The predominantly radial magnetic field created by the MTI at large radii in spherical accretion flows may account for the stable Faraday rotation measure towards Sgr A* in the Galactic Center.

preprint2007arXiv

Electron Heating in Hot Accretion Flows

Local (shearing box) simulations of the nonlinear evolution of the magnetorotational instability in a collisionless plasma show that angular momentum transport by pressure anisotropy ($p_\perp \ne p_\parallel$, where the directions are defined with respect to the local magnetic field) is comparable to that due to the Maxwell and Reynolds stresses. Pressure anisotropy, which is effectively a large-scale viscosity, arises because of adiabatic invariants related to $p_\perp$ and $p_\parallel$ in a fluctuating magnetic field. In a collisionless plasma, the magnitude of the pressure anisotropy, and thus the viscosity, is determined by kinetic instabilities at the cyclotron frequency. Our simulations show that $\sim 50$ % of the gravitational potential energy is directly converted into heat at large scales by the viscous stress (the remaining energy is lost to grid-scale numerical dissipation of kinetic and magnetic energy). We show that electrons receive a significant fraction ($\sim [T_e/T_i]^{1/2}$) of this dissipated energy. Employing this heating by an anisotropic viscous stress in one dimensional models of radiatively inefficient accretion flows, we find that the radiative efficiency of the flow is greater than 0.5% for $\dot{M} \gtrsim 10^{-4} \dot{M}_{Edd}$. Thus a low accretion rate, rather than just a low radiative efficiency, is necessary to explain the low luminosity of many accreting black holes. For Sgr A* in the Galactic Center, our predicted radiative efficiencies imply an accretion rate of $\approx 3 \times 10^{-8} M_\odot {\rm yr^{-1}}$ and an electron temperature of $\approx 3 \times 10^{10}$ K at $\approx 10$ Schwarzschild radii; the latter is consistent with the brightness temperature inferred from VLBI observations.

preprint2007arXiv

Faraday Rotation in Global Accretion Disk Simulations: Implications for Sgr A*

These Faraday rotation calculations of hot, thick accretion flows are motivated by the measured steady rotation measure (RM) of $\approx -6 \times 10^5$ rad m$^{-2}$ from Sgr A*. In our numerical simulations, the quasi-steady state structure of the accretion flow, and the RM it produces, depends on the initial magnetic field. In spite of this dependence, we can draw several robust conclusions about Faraday rotation produced by geometrically thick accretion disks: i) the time averaged RM does not depend that sensitively on the viewing angle, but the stability of the RM can. Equatorial viewing angles show significant variability in RM (including sign reversals), while polar viewing angles are relatively stable if there is a large scale magnetic field threading the disk at large radii. ii) Most of the RM is produced at small radii for polar viewing angles while all radii contribute significantly near the midplane of the disk. Our simulations confirm previous analytic arguments that the accretion rate onto Sgr A* must satisfy $\dot M_{\rm in} \ll \dot M_{\rm Bondi} \sim 10^{-5} \mpy$ in order to not over-produce the measured RM. We argue that the steady RM $\approx -6 \times 10^5$ rad m$^{-2}$ from Sgr A* has two plausible explanations: 1) it is produced at $\sim 100$ Schwarzschild radii, requires $\dot{M}_{\rm in} \approx 3 \times 10^{-8} M_\odot$ yr$^{-1}$, and we view the flow at an angle of $\sim 30^\circ$ relative to the rotation axis of the disk; in our simulations, the variation in RM across a finite-sized source is sufficient to depolarize the emission below $\approx$ 100 GHz, consistent with observations. 2) Alternatively, the RM may be produced in the relatively spherical inflowing plasma near the circularization radius at $\sim 10^3-10^4$ Schwarzschild radii.

preprint2007arXiv

Luminosity Dependence in the Fundamental Plane Projections of Elliptical Galaxies

We analyze the fundamental plane projections of elliptical galaxies as a function of luminosity, using a sample of approximately 80,000 galaxies drawn from Data Release 4 (DR4) of the Sloan Digital Sky Survey (SDSS). We separate brightest cluster galaxies (BCGs) from our main sample and reanalyze their photometry due to a problem with the default pipeline sky subtraction for BCGs. The observables we consider are effective radius (R_e), velocity dispersion (sigma), dynamical mass (M_dyn ~ R_e sigma2), effective density (sigma2/R_e2), and effective surface brightness (mu_e). With the exception of the L-M_dyn correlation, we find evidence of variations in the slope (i.e. the power-law index) of the fundamental plane projections with luminosity for our normal elliptical galaxy population. In particular, the radius-luminosity and Faber-Jackson relations are steeper at high luminosity relative to low luminosity, and the more luminous ellipticals become progressively less dense and have lower surface brightnesses than lower luminosity ellipticals. These variations can be understood as arising from differing formation histories, with more luminous galaxies having less dissipation. Data from the literature and our reanalysis of BCGs show that BCGs have radius-luminosity and Faber-Jackson relations steeper than the brightest non-BCG ellipticals in our sample, consistent with significant growth of BCGs via dissipationless mergers. The variations in slope we find in the Faber-Jackson relation of non-BCGs are qualitatively similar to that reported in the black hole mass-velocity dispersion (M_BH-sigma) correlation. This similarity is consistent with a roughly constant value of M_BH/M_star over a wide range of early type galaxies, where M_star is the stellar mass.

preprint2007arXiv

Proto-Neutron Star Winds with Magnetic Fields and Rotation

We solve the 1D neutrino-heated non-relativistic MHD wind problem for conditions that range from slowly rotating (spin period P > 10 ms) protoneutron stars (PNSs) with surface field strengths typical of radio pulsars (B < 10^13 G), to "proto-magnetars" with B ~ 10^14-10^15 G in their hypothesized rapidly rotating initial states (P ~ 1 ms). We use the simulations of Bucciantini et al. (2006) to map our monopole results onto a more physical dipole geometry and to estimate the spindown of PNSs when their winds are relativistic. We then quantify the effects of rotation and magnetic fields on the mass loss, energy loss, and r-process nucleosynthesis in PNS winds. We describe the evolution of PNS winds through the Kelvin-Helmholtz cooling epoch, emphasizing the transition between (1) thermal neutrino-driven, (2) non-relativistic magnetically-dominated, and (3) relativistic magnetically-dominated outflows. We find that proto-magnetars with P ~ 1 ms and B > 10^15 G drive relativistic winds with luminosities, energies, and Lorentz factors (magnetization sigma ~ 0.1-1000) consistent with those required to produce long duration gamma-ray bursts and hyper-energetic supernovae (SNe). A significant fraction of the rotational energy may be extracted in only a few seconds, sufficiently rapidly to alter the asymptotic energy of the SN remnant, its morphology, and, potentially, its nucleosynthetic yield. Winds from PNSs with more modest rotation periods (2 - 10 ms) and with magnetar-strength fields produce conditions significantly more favorable for the r-process than winds from slowly rotating PNSs. Lastly, we show that energy and momentum deposition by convectively-excited waves further increase the likelihood of successful r-process in PNS winds.

preprint2007arXiv

The Effects of Thermal Conduction on Radiatively-Inefficient Accretion Flows

We quantify the effects of electron thermal conduction on the properties of hot accretion flows, under the assumption of spherical symmetry. Electron heat conduction is important for low accretion rate systems where the electron cooling time is longer than the conduction time of the plasma, such as Sgr A* in the Galactic Center. For accretion flows with density profiles similar to the Bondi solution (n[r] ~ r^[-3/2]), we show that heat conduction leads to super-virial temperatures, implying that conduction significantly modifies the dynamics of the accretion flow. We then self-consistently solve for the dynamics of spherical accretion in the presence of saturated conduction and electron heating. We find that the accretion rate onto the central object can be reduced by ~1-3 orders of magnitude relative to the canonical Bondi rate. Electron conduction may thus be an important ingredient in explaining the low radiative efficiencies and low accretion rates inferred from observations of low-luminosity galactic nuclei. The solutions presented in this paper may also describe the nonlinear saturation of the magnetothermal instability in hot accretion flows.

preprint2006arXiv

The Ionization State of Sodium in Galactic Winds

Roughly 80% of Ultraluminous Infrared Galaxies (ULIRGs) show blue shifted absorption in the resonance lines of neutral sodium, indicating that cool winds are common in such objects, as shown by Rupke et al and by Martin. The neutral sodium (NaI) columns indicated by these absorption lines are ~ 10^{13}-3x10^{14}/cm^2, while the bolometric luminosity varies by a factor of only four. We show that the gas in ULIRG outflows is likely to be in photoionization equilibrium. The very small ULIRG sample of Goldader et al. demonstrates that the ratio of ultraviolet flux to far infrared flux varies by a factor $\sim100$ from object to object. While the Goldader sample does not overlap with those of Rupke et al. and Martin, we show that such a large variation in ultraviolet flux will produce a similar variation in the column of neutral sodium for a fixed mass flux and density. However, if the cold gas is in pressure equilibrium with a hot outflow with a mass loss rate similar to the star formation rate, the range of ionization state is significantly smaller. Measurements of the UV flux for objects in the Martin and Rupke et al. catalogs will definitively determine if photoionization effects are responsible for the wide variation seen in the sodium columns. If they are, a determination of the gas density and mass loss rate in the cool winds will follow, with attendant improvements in our understanding of wind driving mechanisms and of the effects of galaxies on their surroundings.

preprint2005arXiv

Radiation Pressure Supported Starburst Disks and AGN Fueling

We consider the structure of marginally Toomre-stable starburst disks under the assumption that radiation pressure on dust grains provides the dominant vertical support against gravity. This is particularly appropriate when the disk is optically thick to its own IR radiation, as in the central regions of ULIRGs. Because the disk radiates at its Eddington limit, the Schmidt-law for star formation changes in the optically-thick limit, with the star formation rate per unit area scaling as Sigma_g/kappa, where Sigma_g is the gas surface density and kappa is the mean opacity. We show that optically thick starburst disks have a characteristic flux and dust effective temperature of F ~ 10^{13} L_sun/kpc^2 and T_eff ~ 90K, respectively. We compare our predictions with observations and find good agreement. We extend our model from many-hundred parsec scales to sub-parsec scales and address the problem of fueling AGN. We assume that angular momentum transport proceeds via global torques rather than a local viscosity. We account for the radial depletion of gas due to star formation and find a strong bifurcation between two classes of disk models: (1) solutions with a starburst on large scales that consumes all of the gas with little fueling of a central AGN and (2) models with an outer large-scale starburst accompanied by a more compact starburst on 1-10 pc scales and a bright central AGN. The luminosity of the latter models is in many cases dominated by the AGN. We show that the vertical thickness of the starburst disk on pc scales can approach h ~ r, perhaps accounting for the nuclear obscuration in some Type 2 AGN. We also argue that the disk of young stars in the Galactic Center may be the remnant of such a compact nuclear starburst.

preprint2004arXiv

Magnetar Spindown, Hyper-Energetic Supernovae, and Gamma Ray Bursts

The Kelvin-Helmholtz cooling epoch, lasting tens of seconds after the birth of a neutron star in a successful core-collapse supernova, is accompanied by a neutrino-driven wind. For magnetar-strength ($\sim10^{15}$ G) large scale surface magnetic fields, this outflow is magnetically-dominated during the entire cooling epoch.Because the strong magnetic field forces the wind to co-rotate with the protoneutron star,this outflow can significantly effect the neutron star's early angular momentum evolution, as in analogous models of stellar winds (e.g. Weber & Davis 1967). If the rotational energy is large in comparison with the supernova energy and the spindown timescale is short with respect to the time required for the supernova shockwave to traverse the stellar progenitor, the energy extracted may modify the supernova shock dynamics significantly. This effect is capable of producing hyper-energetic supernovae and, in some cases, provides conditions favorable for gamma ray bursts. We estimate spindown timescales for magnetized, rotating protoneutron stars and construct steady-state models of neutrino-magnetocentrifugally driven winds. We find that if magnetars are born rapidly rotating, with initial spin periods ($P$) of $\sim1$ millisecond, that of order $10^{51}-10^{52}$ erg of rotational energy can be extracted in $\sim10$ seconds. If magnetars are born slowly rotating ($P\gtrsim10$ ms) they can spin down to periods of $\sim1$ second on the Kelvin-Helmholtz timescale.

preprint2004arXiv

The distribution and cosmic evolution of massive black hole spins

We study the expected distribution of massive black hole (MBH) spins and its evolution with cosmic time in the context of hierarchical galaxy formation theories. Our model uses Monte Carlo realizations of the merger hierarchy in a LCDM cosmology, coupled to semi-analytical recipes, to follow the merger history of dark matter halos, the dynamics of the MBHs they host, and their growth via gas accretion and binary coalescences. The coalescence of comparable mass holes increases the spin of MBHs, while the capture of smaller companions in randomly-oriented orbits acts to spin holes down. We find that, given the distribution of MBH binary mass ratios in hierarchical models, binary coalescences alone do not lead to a systematic spin-up or spin-down of MBHs with time: the spin distribution retains memory of its initial conditions. By contrast, because of the Bardeen-Petterson effect, gas accretion via a thin disk tends to spin holes up even if the direction of the spin axis changes randomly in time. In our models, accretion dominates over black hole captures and efficiently spins holes up. The spin distribution is heavily skewed towards fast-rotating Kerr holes, is already in place at early epochs, and does not change much below redshift 5. If accretion is via a thin disk, about 70% of all MBHs are maximally rotating and have radiative efficiencies approaching 30% (assuming a "standard'' spin-efficiency conversion). Even in the conservative case where accretion is via a geometrically thick disk, about 80% of all MBHs have spin parameters a/m > 0.8 and accretion efficiencies > 12%. Rapidly spinning holes with high radiative efficiencies may satisfy constraints based on comparing the local MBH mass density with the mass density inferred from luminous quasars (Soltan's argument).

preprint2004arXiv

The Formation and Evolution of the First Massive Black Holes

The first massive astrophysical black holes likely formed at high redshifts (z>10) at the centers of low mass (~10^6 Msun) dark matter concentrations. These black holes grow by mergers and gas accretion, evolve into the population of bright quasars observed at lower redshifts, and eventually leave the supermassive black hole remnants that are ubiquitous at the centers of galaxies in the nearby universe. The astrophysical processes responsible for the formation of the earliest seed black holes are poorly understood. The purpose of this review is threefold: (1) to describe theoretical expectations for the formation and growth of the earliest black holes within the general paradigm of hierarchical cold dark matter cosmologies, (2) to summarize several relevant recent observations that have implications for the formation of the earliest black holes, and (3) to look into the future and assess the power of forthcoming observations to probe the physics of the first active galactic nuclei.

preprint1997arXiv

Are Particles in Advection-Dominated Accretion Flows Thermal?

We investigate the form of the momentum distribution function for protons and electrons in an advection-dominated accretion flow (ADAF). We show that for all accretion rates, Coulomb collisions are too inefficient to thermalize the protons. The proton distribution function is therefore determined by the viscous heating mechanism, which is unknown. The electrons, however, can exchange energy quite efficiently through Coulomb collisions and the emission and absorption of synchrotron photons. We find that for accretion rates greater than \sim 10^{-3} of the Eddington accretion rate, the electrons have a thermal distribution throughout the accretion flow. For lower accretion rates, the electron distribution function is determined by the electron's source of heating, which is primarily adiabatic compression. Using the principle of adiabatic invariance, we show that an adiabatically compressed collisionless gas maintains a thermal distribution until the particle energies become relativistic. We derive a new, non-thermal, distribution function which arises for relativistic energies and provide analytic formulae for the synchrotron radiation from this distribution. Finally, we discuss its implications for the emission spectra from ADAFs.

preprint1995arXiv

Observational Searches for Solar g-modes: Some Theoretical Considerations

We argue that the solar g-modes are unlikely to have caused the discrete peaks in the power spectrum of the solar wind flux observed by Thomson et al. (1995). The lower limit to the energy of individual g-modes, using the amplitudes given by Thomson et al., is estimated to be at least 10$^{36}$ erg for low order g-modes; the resulting surface velocity amplitude is at least 50 cm s$^{-1}$, larger than the observational upper limit (5 cm s$^{-1}$). We suggest that the most likely source for the excitation of solar g-modes is turbulent stresses in the convection zone. The surface velocity amplitude of low degree and low order g-modes resulting from this process is estimated to be of order 10$^{-2}$ cm s$^{-1}$. This amplitude is interestingly close to the detection threshold of the SOHO satellite. The long lifetime of g-modes ($\sim 10^6$ years for low order modes) should be helpful in detecting these small amplitude pulsations.

preprint1995arXiv

On the Validity of the Classical Apsidal Motion Formula for Tidal Distortion

We check the validity of the widely used classical apsidal motion formula as a function of orbital parameters, stellar structure, and stellar rotation rate by comparing dynamical calculations of the periastron advance with the static tidal formula. We find that the classical formula gives very accurate results when the periods of the low order quadrupole g, f and p modes are smaller than the periastron passage time by a factor of about 7 or more. However, when this condition is not satisfied, the difference between the classical formula and the exact result can be quite large, and even periastron recession can result. The largest difference arises when one of the low order modes of the star is nearly resonant with an integer multiple of the orbital frequency minus twice the rotation rate of the star. The resonance of higher order g-modes (number of radial nodes $\gta 4$) with the orbit is very unlikely to cause significant deviation from the classical result because of their weak coupling to the tidal force and thus their small contribution to the apsidal motion. Resonances involving rotational modes of the star are also unlikely to make much contribution to the apsidal motion because of their small overlap with the tidal force, even though they have periods comparable to the periastron passage time. We apply our work to two famous binary systems (AS Cam and DI Her) which show abnormally small apsidal motion, and conclude that dynamical effects are unimportant for these systems, i.e. the static tide assumption is an excellent approximation.