Showing papers by "Eliot Quataert published in 2017"

Convection Destroys the Core/Mantle Structure in Hybrid C/O/Ne White Dwarfs

Abstract: A hybrid C/O/Ne white dwarf (WD)—an unburned C/O core surrounded by an O/Ne/Na mantle—can be formed if the carbon flame is quenched in a super-AGB star or white dwarf merger remnant. We show that this segregated hybrid structure becomes unstable to rapid mixing within 2000 years of the onset of WD cooling. Carbon burning includes a weak reaction that removes electrons, resulting in a lower electron-to-baryon ratio () in the regions processed by carbon burning compared to the unburned C/O core, making the O/Ne mantle denser than the C/O core as the WD cools. This is unstable to efficient mixing. We use the results of models with different size C/O cores to quantify the rate at which the cores mix with the mantle as they cool. In all cases, we find that the WDs undergo significant core/mantle mixing on timescales shorter than the time available to grow the WD to the Chandrasekhar mass (MCh) by accretion. As a result, hybrid WDs that reach MCh due to later accretion will have lower central carbon fractions than assumed thus far. We briefly discuss the implications of these results for the possibility of SNe Ia from hybrid WDs.

Accretion-induced Collapse from Helium Star + White Dwarf Binaries

Abstract: Accretion-induced collapse (AIC) occurs when an O/Ne white dwarf (WD) grows to nearly the Chandrasekhar mass ($M_{\rm Ch}$), reaching central densities that trigger electron captures in the core. Using Modules for Experiments in Stellar Astrophysics ($\texttt{MESA}$), we present the first true binary simulations of He star + O/Ne WD binaries, focusing on a $1.5 M_\odot$ He star in a 3 hour orbital period with $1.1-1.3 M_\odot$ O/Ne WDs. The helium star fills its Roche lobe after core helium burning is completed and donates helium on its thermal timescale to the WD, $\dot{M}\approx3\times10^{-6} M_\odot$/yr, a rate high enough that the accreting helium burns stably on the WD. The accumulated carbon/oxygen ashes from the helium burning undergo an unstable shell flash that initiates an inwardly moving carbon burning flame. This flame is only quenched when it runs out of carbon at the surface of the original O/Ne core. Subsequent accumulation of fresh carbon/oxygen layers also undergo thermal instabilities, but no mass loss is triggered, allowing $M_{\rm WD}\rightarrow M_{\rm Ch}$, triggering the onset of AIC. We also discuss the scenario of accreting C/O WDs that experience shell carbon ignitions to become O/Ne WDs, and then, under continuing mass transfer, lead to AIC. Studies of the AIC event rate using binary population synthesis should include all of these channels, especially this latter channel, which has been previously neglected but might dominate the rate.

Testing the Recovery of Intrinsic Galaxy Sizes and Masses of z ̃ 2 Massive Galaxies Using Cosmological Simulations

Abstract: Accurate measurements of galaxy masses and sizes are key to tracing galaxy evolution over time. Cosmological zoom-in simulations provide an ideal test bed for assessing the recovery of galaxy properties from observations. Here, we utilize galaxies with M_* ~ 10^10 - 10^(11.5) M_☉ at z ~ 1.7–2 from the MassiveFIRE cosmological simulation suite, part of the Feedback in Realistic Environments (FIRE) project. Using mock multi-band images, we compare intrinsic galaxy masses and sizes to observational estimates. We find that observations accurately recover stellar masses, with a slight average underestimate of ~ 0.06 dex and a ~ 0.15 dex scatter. Recovered half-light radii agree well with intrinsic half-mass radii when averaged over all viewing angles, with a systematic offset of ~ 0.1 dex (with the half-light radii being larger) and a scatter of ~ 0.2 dex. When using color gradients to account for mass-to-light variations, recovered half-mass radii also exceed the intrinsic half-mass radii by ~ 1 dex. However, if not properly accounted for, aperture effects can bias size estimates by ~ 0.1 dex. No differences are found between the mass and size offsets for star-forming and quiescent galaxies. Variations in viewing angle are responsible for ~25% of the scatter in the recovered masses and sizes. Our results thus suggest that the intrinsic scatter in the mass–size relation may have previously been overestimated by ~25%. Moreover, orientation-driven scatter causes the number density of very massive galaxies to be overestimated by ~ 0.5 dex at M_* ~ 10^(11.5) M_☉_.

Testing the Recovery of Intrinsic Galaxy Sizes and Masses of z~2 Massive Galaxies Using Cosmological Simulations

Abstract: Accurate measurements of galaxy masses and sizes are key to tracing galaxy evolution over time. Cosmological zoom-in simulations provide an ideal test bed for assessing the recovery of galaxy properties from observations. Here, we utilize galaxies with $M_*\sim10^{10}-10^{11.5}M_{\odot}$ at z~1.7-2 from the MassiveFIRE cosmological simulation suite, part of the Feedback in Realistic Environments (FIRE) project. Using mock multi-band images, we compare intrinsic galaxy masses and sizes to observational estimates. We find that observations accurately recover stellar masses, with a slight average underestimate of ~0.06 dex and a ~0.15 dex scatter. Recovered half-light radii agree well with intrinsic half-mass radii when averaged over all viewing angles, with a systematic offset of ~0.1 dex (with the half-light radii being larger) and a scatter of ~0.2 dex. When using color gradients to account for mass-to-light variations, recovered half-mass radii also exceed the intrinsic half-mass radii by ~0.1 dex. However, if not properly accounted for, aperture effects can bias size estimates by ~0.1 dex. No differences are found between the mass and size offsets for star-forming and quiescent galaxies. Variations in viewing angle are responsible for ~25% of the scatter in the recovered masses and sizes. Our results thus suggest that the intrinsic scatter in the mass-size relation may have previously been overestimated by ~25%. Moreover, orientation-driven scatter causes the number density of very massive galaxies to be overestimated by ~0.5 dex at $M_*\sim10^{11.5}M_{\odot}$.

PIC Simulations of Velocity-Space Instabilities in a Decreasing Magnetic Field: Viscosity and Thermal Conduction

Abstract: We use particle-in-cell (PIC) simulations of a collisionless, electron-ion plasma with a decreasing background magnetic field, $B$, to study the effect of velocity-space instabilities on the viscous heating and thermal conduction of the plasma. If $B$ decreases, the adiabatic invariance of the magnetic moment gives rise to pressure anisotropies with $p_{||,j} > p_{\perp,j}$ ($p_{||,j}$ and $p_{\perp,j}$ represent the pressure of species $j$ ($=i$ or $e$) parallel and perpendicular to the magnetic field). Linear theory indicates that, for sufficiently large anisotropies, different velocity-space instabilities can be triggered. These instabilities, which grow on scales comparable to the electron and ion Larmor radii, in principle have the ability to pitch-angle scatter the particles, limiting the growth of the anisotropies. Our PIC simulations focus on the nonlinear, saturated regime of the instabilities. This is done through the permanent decrease of the magnetic field by an imposed shear in the plasma. Our results show that, in the regime $2 \lesssim \beta_j \lesssim 20$ ($\beta_j \equiv 8\pi p_j/B^2$), the saturated ion and electron pressure anisotropies are controlled by the combined effect of the oblique ion firehose (OIF) and the fast magnetosonic/whistler (FM/W) instabilities. These instabilities grow preferentially on the ion Larmor radius scale, and make the ion and electron pressure anisotropies nearly equal: $\Delta p_e/p_{||,e} \approx \Delta p_i/p_{||,i}$ (where $\Delta p_j=p_{\perp,j} - p_{||,j}$). We also quantify the thermal conduction of the plasma by directly calculating the mean free path of electrons along the mean magnetic field, which we find strongly depends on whether $B$ decreases or increases. Our results can be applied in studies of low collisionality plasmas such as the solar wind, the intracluster medium, and some accretion disks around black holes.

The Effects of Magnetic Fields on the Dynamics of Radiation Pressure Dominated Massive Star Envelopes

Abstract: We use three dimensional radiation magneto-hydrodynamic simulations to study the effects of magnetic fields on the energy transport and structure of radiation pressure dominated main sequence massive star envelopes at the region of the iron opacity peak. We focus on the regime where the local thermal timescale is shorter than the dynamical timescale, corresponding to inefficient convective energy transport. We begin with initially weak magnetic fields relative to the thermal pressure, from 100-1000G in differing geometries. The unstable density inversion amplifies the magnetic field, increasing the magnetic energy density to values close to equipartition with the turbulent kinetic energy density. By providing pressure support, the magnetic field's presence significantly increases the density fluctuations in the turbulent envelope, thereby enhancing the radiative energy transport by allowing photons to diffuse out through low density regions. Magnetic buoyancy brings small scale magnetic fields to the photosphere and increases the vertical energy transport with the energy advection velocity proportional to the Alfv\'en velocity, although in all cases we study photon diffusion still dominates the energy transport. The increased radiative and advective energy transport causes the stellar envelope to shrink by several scale heights. We also find larger turbulent velocity fluctuations compared to the purely hydrodynamic case, reaching $\approx$ 100 km/s at the stellar photosphere. The photosphere also shows vertical oscillations with similar averaged velocities and periods of a few hours. The increased turbulent velocity and oscillations will have strong impacts on the line broadening and periodic signals in massive stars.

Fast and Luminous Transients from the Explosions of Long-lived Massive White Dwarf Merger Remnants

Abstract: We study the evolution and final outcome of long-lived (${\approx}10^5$ years) remnants from the merger of a He white dwarf (WD) with a more massive C/O or O/Ne WD. Using Modules for Experiments in Stellar Astrophysics ($\texttt{MESA}$), we show that these remnants have a red giant configuration supported by steady helium burning, adding mass to the WD core until it reaches $M_{\rm core}\approx 1.12-1.20 M_\odot$. At that point, the base of the surface convection zone extends into the burning layer, mixing the helium burning products (primarily carbon and magnesium) throughout the convective envelope. Further evolution depletes the convective envelope of helium, and dramatically slows the mass increase of the underlying WD core. The WD core mass growth re-initiates after helium depletion, as then an uncoupled carbon burning shell is ignited and proceeds to burn the fuel from the remaining metal-rich extended envelope. For large enough initial total merger masses, O/Ne WD cores would experience electron-capture triggered collapse to neutron stars (NSs) after growing to near Chandrasekhar mass ($M_{\rm Ch}$). Massive C/O WD cores could suffer the same fate after a carbon-burning flame converts them to O/Ne. The NS formation would release ${\approx}10^{50}$ ergs into the remaining extended low mass envelope. Using the STELLA radiative transfer code, we predict the resulting optical light curves from these exploded envelopes. Reaching absolute magnitudes of $M_V\approx -17$, these transients are bright for about one week, and have many features of the class of luminous, rapidly evolving transients studied by Drout and collaborators.

Pressure-anisotropy-induced nonlinearities in the kinetic magnetorotational instability.

Abstract: In collisionless and weakly collisional plasmas, such as hot accretion flows onto compact objects, the magnetorotational instability (MRI) can differ significantly from the standard (collisional) MRI. In particular, pressure anisotropy with respect to the local magnetic-field direction can both change the linear MRI dispersion relation and cause nonlinear modifications to the mode structure and growth rate, even when the field and flow perturbations are very small. This work studies these pressure-anisotropy-induced nonlinearities in the weakly nonlinear, high-ion-beta regime, before the MRI saturates into strong turbulence. Our goal is to better understand how the saturation of the MRI in a low-collisionality plasma might differ from that in the collisional regime. We focus on two key effects: (i) the direct impact of self-induced pressure-anisotropy nonlinearities on the evolution of an MRI mode, and (ii) the influence of pressure anisotropy on the ‘parasitic instabilities’ that are suspected to cause the mode to break up into turbulence. Our main conclusions are: (i) The mirror instability regulates the pressure anisotropy in such a way that the linear MRI in a collisionless plasma is an approximate nonlinear solution once the mode amplitude becomes larger than the background field (just as in magnetohyrodynamics). This implies that differences between the collisionless and collisional MRI become unimportant at large amplitudes. (ii) The break up of large-amplitude MRI modes into turbulence via parasitic instabilities is similar in collisionless and collisional plasmas. Together, these conclusions suggest that the route to magnetorotational turbulence in a collisionless plasma may well be similar to that in a collisional plasma, as suggested by recent kinetic simulations. As a supplement to these findings, we offer guidance for the design of future kinetic simulations of magnetorotational turbulence.

Accretion-Induced Collapse From Helium Star + White Dwarf Binaries

Abstract: Accretion-induced collapse (AIC) occurs when an O/Ne white dwarf (WD) grows to nearly the Chandrasekhar mass ($M_{\rm Ch}$), reaching central densities that trigger electron captures in the core. Using Modules for Experiments in Stellar Astrophysics ($\texttt{MESA}$), we present the first true binary simulations of He star + O/Ne WD binaries, focusing on a $1.5 M_\odot$ He star in a 3 hour orbital period with $1.1-1.3 M_\odot$ O/Ne WDs. The helium star fills its Roche lobe after core helium burning is completed and donates helium on its thermal timescale to the WD, $\dot{M}\approx3\times10^{-6} M_\odot$/yr, a rate high enough that the accreting helium burns stably on the WD. The accumulated carbon/oxygen ashes from the helium burning undergo an unstable shell flash that initiates an inwardly moving carbon burning flame. This flame is only quenched when it runs out of carbon at the surface of the original O/Ne core. Subsequent accumulation of fresh carbon/oxygen layers also undergo thermal instabilities, but no mass loss is triggered, allowing $M_{\rm WD}\rightarrow M_{\rm Ch}$, triggering the onset of AIC. We also discuss the scenario of accreting C/O WDs that experience shell carbon ignitions to become O/Ne WDs, and then, under continuing mass transfer, lead to AIC. Studies of the AIC event rate using binary population synthesis should include all of these channels, especially this latter channel, which has been previously neglected but might dominate the rate.

Fast and Luminous Transients from the Explosions of Long Lived Massive White Dwarf Merger Remnants

Abstract: We study the evolution and final outcome of long-lived (${\approx}10^5$ years) remnants from the merger of a He white dwarf (WD) with a more massive C/O or O/Ne WD. Using Modules for Experiments in Stellar Astrophysics ($\texttt{MESA}$), we show that these remnants have a red giant configuration supported by steady helium burning, adding mass to the WD core until it reaches $M_{\rm core}\approx 1.12-1.20 M_\odot$. At that point, the base of the surface convection zone extends into the burning layer, mixing the helium burning products (primarily carbon and magnesium) throughout the convective envelope. Further evolution depletes the convective envelope of helium, and dramatically slows the mass increase of the underlying WD core. The WD core mass growth re-initiates after helium depletion, as then an uncoupled carbon burning shell is ignited and proceeds to burn the fuel from the remaining metal-rich extended envelope. For large enough initial total merger masses, O/Ne WD cores would experience electron-capture triggered collapse to neutron stars (NSs) after growing to near Chandrasekhar mass ($M_{\rm Ch}$). Massive C/O WD cores could suffer the same fate after a carbon-burning flame converts them to O/Ne. The NS formation would release ${\approx}10^{50}$ ergs into the remaining extended low mass envelope. Using the STELLA radiative transfer code, we predict the resulting optical light curves from these exploded envelopes. Reaching absolute magnitudes of $M_V\approx -17$, these transients are bright for about one week, and have many features of the class of luminous, rapidly evolving transients studied by Drout and collaborators.

Stochastic Electron Acceleration by the Whistler Instability in a Growing Magnetic Field

Abstract: We use 2D particle-in-cell (PIC) simulations to study the effect of the saturated whistler instability on the viscous heating and nonthermal acceleration of electrons in a shearing, collisionless plasma with a growing magnetic field, \textbf{B}. In this setup, an electron pressure anisotropy with $p_{\perp,e} > p_{||,e}$ naturally arises due to the adiabatic invariance of the electron magnetic moment ($p_{||,e}$ and $p_{\perp,e}$ are the pressures parallel and perpendicular to \textbf{B}). If the anisotropy is large enough, the whistler instability arises, efficiently scattering the electrons and limiting $\Delta p_e$ ($\equiv p_{\perp,e}-p_{||,e}$). In this context, $\Delta p_e$ taps into the plasma velocity shear, producing electron heating by the so called anisotropic viscosity. In our simulations, we permanently drive the growth of $|\textbf{B}|$ by externally imposing a plasma shear, allowing us to self-consistently capture the long-term, saturated whistler instability evolution. We find that besides the viscous heating, the scattering by whistler modes can stochastically accelerate electrons to nonthermal energies. This acceleration is most prominent when initially $\beta_e\sim 1$, gradually decreasing its efficiency for larger values of $\beta_e$ ($\equiv 8\pi p_e/|\textbf{B}|^2$). If initially $\beta_e \sim 1$, the final electron energy distribution can be approximately described by a thermal component, plus a power-law tail with spectral index $\sim 3.7$. In these cases, the nonthermal tail accounts for $\sim 5\%$ of the electrons, and for $\sim 15\%$ of their kinetic energy. We discuss the implications of our results for electron heating and acceleration in low-collisionality astrophysical environments, such as low-luminosity accretion flows.

Stochastic Electron Acceleration by the Whistler Instability in a Growing Magnetic Field

Abstract: UC Berkeley-Chile Fund / NSF, AST 13-33612, 17-15054 / Simons Foundation / David and Lucile Packard Foundation / National Science Foundation, ACI-1053575 / supercomputing infrastructure of the NLHPC at the Center for Mathematical Modeling of University of Chile ECM-02

Pressure-anisotropy-induced nonlinearities in the kinetic magnetorotational instability

Abstract: In collisionless and weakly collisional plasmas, such as hot accretion flows onto compact objects, the magnetorotational instability (MRI) can differ significantly from the standard (collisional) MRI. In particular, pressure anisotropy with respect to the local magnetic-field direction can both change the linear MRI dispersion relation and cause nonlinear modifications to the mode structure and growth rate, even when the field and flow perturbations are small. This work studies these pressure-anisotropy-induced nonlinearities in the weakly nonlinear, high-ion-beta regime, before the MRI saturates into strong turbulence. Our goal is to better understand how the saturation of the MRI in a low collisionality plasma might differ from that in the collisional regime. We focus on two key effects: (i) the direct impact of self-induced pressure-anisotropy nonlinearities on the evolution of an MRI mode, and (ii) the influence of pressure anisotropy on the "parasitic instabilities" that are suspected to cause the mode to break up into turbulence. Our main conclusions are: (i) The mirror instability regulates the pressure anisotropy in such a way that the linear MRI in a collisionless plasma is an approximate nonlinear solution once the mode amplitude becomes larger than the background field (just as in MHD). This implies that differences between the collisionless and collisional MRI become unimportant at large amplitudes. (ii) The break up of large amplitude MRI modes into turbulence via parasitic instabilities is similar in collisionless and collisional plasmas. Together, these conclusions suggest that the route to magnetorotational turbulence in a collisionless plasma may well be similar to that in a collisional plasma, as suggested by recent kinetic simulations. As a supplement to these findings, we offer guidance for the design of future kinetic simulations of magnetorotational turbulence.