Showing papers by "Eliot Quataert published in 2019"

Low-frequency Variability in Massive Stars: Core Generation or Surface Phenomenon?

Abstract: 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.

A Resolution Study of Magnetically Arrested Disks

Abstract: We investigate numerical convergence in simulations of magnetically arrested disks around spinning black holes. Using the general-relativistic magnetohydrodynamics code Athena++, we study the same system at four resolutions (up to an effective 512 × 256 × 512 cells) and with two different spatial reconstruction algorithms. The accretion rate and general large-scale structure of the flow agree across the simulations. This includes the amount of magnetic flux accumulated in the saturated state and the ensuing suppression of the magnetorotational instability from the strong field. The energy of the jet and the efficiency with which spin energy is extracted via the Blandford–Znajek process also show convergence. However the spatial structure of the jet shows variation across the set of grids employed, as do the Lorentz factors. Small-scale features of the turbulence, as measured by correlation lengths, are not fully converged. Despite convergence of a number of aspects of the flow, modeling of synchrotron emission shows that variability is not converged and decreases with increasing resolution even at our highest resolutions.

Tilted Disks around Black Holes: A Numerical Parameter Survey for Spin and Inclination Angle

Abstract: We conduct a systematic study of the properties of tilted accretion flows around spinning black holes, covering a range of tilt angles and black hole spins, using the general-relativistic magnetohydrodynamics code Athena++. The same initial magnetized torus is evolved around black holes with spins ranging from 0 to 0.9, with inclinations ranging from 0 degrees to 24 degrees. The tilted disks quickly reach a warped and twisted shape that rigidly precesses about the black hole spin axis with deformations in shape large enough to hinder the application of linear bending wave theory. Magnetized polar outflows form, oriented along the disk rotation axes. At sufficiently high inclinations a pair of standing shocks develops in the disks. These shocks dramatically affect the flow at small radii, driving angular momentum transport. At high spins they redirect material more effectively than they heat it, reducing the dissipation rate relative to the mass accretion rate and lowering the radiative efficiency of the flow.

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

Abstract: 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.

Tilted Disks around Black Holes: A Numerical Parameter Survey for Spin and Inclination Angle

Abstract: We conduct a systematic study of the properties of tilted accretion flows around spinning black holes, covering a range of tilt angles and black hole spins, using the general-relativistic magnetohydrodynamics code Athena++. The same initial magnetized torus is evolved around black holes with spins ranging from 0 to 0.9, with inclinations ranging from 0 degrees to 24 degrees. The tilted disks quickly reach a warped and twisted shape that rigidly precesses about the black hole spin axis with deformations in shape large enough to hinder the application of linear bending wave theory. Magnetized polar outflows form, oriented along the disk rotation axes. At sufficiently high inclinations a pair of standing shocks develops in the disks. These shocks dramatically affect the flow at small radii, driving angular momentum transport. At high spins they redirect material more effectively than they heat it, reducing the dissipation rate relative to the mass accretion rate and lowering the radiative efficiency of the flow.

Weak Shock Propagation with Accretion. II. Stability of Self-similar Solutions to Radial Perturbations

Abstract: Coughlin et al. (2018) (Paper I) derived and analyzed a new regime of self-similarity that describes weak shocks (Mach number of order unity) in the gravitational field of a point mass. These solutions are relevant to low energy explosions, including failed supernovae. In this paper, we develop a formalism for analyzing the stability of shocks to radial perturbations, and we demonstrate that the self-similar solutions of Paper I are extremely weakly unstable to such radial perturbations. Specifically, we show that perturbations to the shock velocity and post-shock fluid quantities (the velocity, density, and pressure) grow with time as $t^{\alpha}$, where $\alpha \le 0.12$, implying that the ten-folding timescale of such perturbations is roughly ten orders of magnitude in time. We confirm these predictions by performing high-resolution, time-dependent numerical simulations. Using the same formalism, we also show that the Sedov-Taylor blastwave is trivially stable to radial perturbations provided that the self-similar, Sedov-Taylor solutions extend to the origin, and we derive simple expressions for the perturbations to the post-shock velocity, density, and pressure. Finally, we show that there is a third, self-similar solution (in addition to the the solutions in Paper I and the Sedov-Taylor solution) to the fluid equations that describes a rarefaction wave, i.e., an outward-propagating sound wave of infinitesimal amplitude. We interpret the stability of shock propagation in light of these three distinct self-similar solutions.

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

Abstract: 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.

Weak Shock Propagation with Accretion. II. Stability of Self-similar Solutions to Radial Perturbations.

Abstract: Coughlin et al. (2018) (Paper I) derived and analyzed a new regime of self-similarity that describes weak shocks (Mach number of order unity) in the gravitational field of a point mass. These solutions are relevant to low energy explosions, including failed supernovae. In this paper, we develop a formalism for analyzing the stability of shocks to radial perturbations, and we demonstrate that the self-similar solutions of Paper I are extremely weakly unstable to such radial perturbations. Specifically, we show that perturbations to the shock velocity and post-shock fluid quantities (the velocity, density, and pressure) grow with time as $t^{\alpha}$, where $\alpha \le 0.12$, implying that the ten-folding timescale of such perturbations is roughly ten orders of magnitude in time. We confirm these predictions by performing high-resolution, time-dependent numerical simulations. Using the same formalism, we also show that the Sedov-Taylor blastwave is trivially stable to radial perturbations provided that the self-similar, Sedov-Taylor solutions extend to the origin, and we derive simple expressions for the perturbations to the post-shock velocity, density, and pressure. Finally, we show that there is a third, self-similar solution (in addition to the the solutions in Paper I and the Sedov-Taylor solution) to the fluid equations that describes a rarefaction wave, i.e., an outward-propagating sound wave of infinitesimal amplitude. We interpret the stability of shock propagation in light of these three distinct self-similar solutions.

A Resolution Study of Magnetically Arrested Disks

Abstract: We investigate numerical convergence in simulations of magnetically arrested disks around spinning black holes. Using the general-relativistic magnetohydrodynamics code Athena++, we study the same system at four resolutions (up to an effective 512 by 256 by 512 cells) and with two different spatial reconstruction algorithms. The accretion rate and general large-scale structure of the flow agree across the simulations. This includes the amount of magnetic flux accumulated in the saturated state and the ensuing suppression of the magnetorotational instability from the strong field. The energy of the jet and the efficiency with which spin energy is extracted via the Blandford-Znajek process also show convergence. However the spatial structure of the jet shows variation across the set of grids employed, as do the Lorentz factors. Small-scale features of the turbulence, as measured by correlation lengths, are not fully converged. Despite convergence of a number of aspects of the flow, modeling of synchrotron emission shows that variability is not converged and decreases with increasing resolution even at our highest resolutions.

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

Abstract: 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.

Gravity and Light: Combining Gravitational Wave and Electromagnetic Observations in the 2020s

Abstract: As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly something completely unexpected. As we build this sample, we will explore exotic astrophysical topics ranging from nucleosynthesis, stellar evolution, general relativity, high-energy astrophysics, nuclear matter, to cosmology. The discovery potential is extraordinary, and investments in this area will yield major scientific breakthroughs. Here we outline some of the most exciting scientific questions that can be answered by combining GW and EM observations.

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

Abstract: 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.