Showing 1–50 of 218 results for author: Blackman, E G

Second-generation planet formation after tidal disruption from common envelope evolution

Authors: Luke Chamandy, Jason Nordhaus, Eric G. Blackman, Emily Wilson

Abstract: We propose that certain white dwarf (WD) planets, such as WD 1856+534 b, may form out of material from a stellar companion that tidally disrupts from common envelope (CE) evolution with the WD progenitor star. The disrupted companion shreds into an accretion disc, out of which a gas giant protoplanet forms due to gravitational instability. To explore this scenario, we make use of detailed stellar… ▽ More We propose that certain white dwarf (WD) planets, such as WD 1856+534 b, may form out of material from a stellar companion that tidally disrupts from common envelope (CE) evolution with the WD progenitor star. The disrupted companion shreds into an accretion disc, out of which a gas giant protoplanet forms due to gravitational instability. To explore this scenario, we make use of detailed stellar evolution models consistent with WD 1856+534. The minimum mass companion that produces a gravitationally-unstable disk after tidal disruption is $\sim$$0.15\,\mathrm{M}_\odot$. Planet formation from tidal disruption is a new channel for producing second-generation planets around WD. △ Less

Submitted 19 July, 2024; originally announced July 2024.

Comments: 8 pages including appendix, 2 figures, 1 table. Submitted to A & A. Comments welcome

Structure and Dynamics of Magneto-Inertial, Differentially Rotating Laboratory Plasmas

Authors: V. Valenzuela-Villaseca, L. G. Suttle, F. Suzuki-Vidal, J. W. D. Halliday, D. R. Russell, S. Merlini, E. R. Tubman, J. D. Hare, J. P. Chittenden, M. E. Koepke, E. G. Blackman, S. V. Lebedev

Abstract: We present a detailed characterization of the structure and evolution of differentially rotating plasmas driven on the MAGPIE pulsed-power generator (1.4 MA peak current, 240 ns rise-time). The experiments were designed to simulate physics relevant to accretion discs and jets on laboratory scales. A cylindrical aluminium wire array Z pinch enclosed by return posts with an overall azimuthal off-set… ▽ More We present a detailed characterization of the structure and evolution of differentially rotating plasmas driven on the MAGPIE pulsed-power generator (1.4 MA peak current, 240 ns rise-time). The experiments were designed to simulate physics relevant to accretion discs and jets on laboratory scales. A cylindrical aluminium wire array Z pinch enclosed by return posts with an overall azimuthal off-set angle was driven to produce ablation plasma flows that propagate inwards in a slightly off-radial trajectory, injecting mass, angular momentum, and confining ram pressure to a rotating plasma column on the axis. However, the plasma is free to expand axially, forming a collimated, differentially rotating axial jet that propagates at $\approx 100$ km/s. The density profile of the jet corresponds to a dense shell surrounding a low-density core, which is consistent with the centrifugal barrier effect being sustained along the jet's propagation. We show analytically that, as the rotating plasma accretes mass, conservation of mass and momentum implies plasma radial growth scaling as $r \propto t^{1/3}$. As the characteristic moment of inertia increases, the rotation velocity is predicted to decrease and settle on a characteristic value $\approx 20$ km/s. We find that both predictions are in agreement with Thomson scattering and optical self-emission imaging measurements. △ Less

Submitted 29 March, 2024; originally announced March 2024.

Comments: Submitted to the Journal of Plasma Physics. 22 pages, 9 figures

Relativistic Jet Motion in the Radio-quiet LINER Galaxy KISSR872

Authors: Preeti Kharb, Eric G. Blackman, Eric Clausen-Brown, Mousumi Das, Daniel A. Schwartz, Aneta Siemiginowska, Smitha Subramanian, Sravani Vaddi

Abstract: We report superluminal jet motion with an apparent speed of $β_\mathrm{app}=1.65\pm0.57$ in the radio-quiet (RQ) low ionisation nuclear emission line region (LINER) galaxy, KISSR872. This result comes from two epoch phase-referenced very long baseline interferometry (VLBI) observations at 5 GHz. The detection of bulk relativistic motion in the jet of this extremely radio faint AGN, with a total 1.… ▽ More We report superluminal jet motion with an apparent speed of $β_\mathrm{app}=1.65\pm0.57$ in the radio-quiet (RQ) low ionisation nuclear emission line region (LINER) galaxy, KISSR872. This result comes from two epoch phase-referenced very long baseline interferometry (VLBI) observations at 5 GHz. The detection of bulk relativistic motion in the jet of this extremely radio faint AGN, with a total 1.4 GHz flux density of 5 mJy in the 5.4 arcsec resolution Very Large Array (VLA) FIRST survey image and 1.5 mJy in the $\sim5$ milli-arcsec resolution Very Long Baseline Array (VLBA) image, is the first of its kind in a RQ LINER galaxy. The presence of relativistic jets in lower accretion rate objects like KISSR872, with an Eddington ratio of 0.04, reveals that even RQ AGN can harbor relativistic jets, and evidentiates their universality over a wide range of accretion powers. △ Less

Submitted 18 December, 2023; originally announced December 2023.

Comments: 9 pages, 4 figures, accepted for publication in ApJ

On The Structure of Plasma Jets in the Rotating Plasma Experiment

Authors: V. Valenzuela-Villaseca, L. G. Suttle, F. Suzuki-Vidal, J. W. D. Halliday, D. R. Russell, S. Merlini, E. R. Tubman, J. D. Hare, J. P. Chittenden, M. E. Koepke, E. G. Blackman, S. V. Lebedev

Abstract: Recent pulsed-power experiments have demonstrated the formation of astrophysically-relevant, differentially rotating plasmas [1]. Key features of the plasma flows are the discovery of a quasi-Keplerian rotation curve, the launching of highly-collimated angular-momentum-transporting axial jets, and a hollow density structure sustained by the centrifugal barrier effect. In this communication we disc… ▽ More Recent pulsed-power experiments have demonstrated the formation of astrophysically-relevant, differentially rotating plasmas [1]. Key features of the plasma flows are the discovery of a quasi-Keplerian rotation curve, the launching of highly-collimated angular-momentum-transporting axial jets, and a hollow density structure sustained by the centrifugal barrier effect. In this communication we discuss several features of the plasma structure in these experiments through order-of-magnitude models. First, we show that the observed rotation velocity would produce a centrifugal force strong enough to support the hollow density profile. Second, we show that the axial jet should diverge much faster than what was observed, were it not for a magnetized halo with 3T which surrounds the jet and exerts pressure on the interface.Finally, we discuss the temperature structure in the axial jet and plasma halo.We show that a 3T magnetic field would also suppress electron heat conduction,leading to the flat profile observed experimentally. We also find that the axial jet is efficiently radiatively cooled,whereas the halo is not, which would explain the thermal decoupling between the two regions. △ Less

Submitted 15 February, 2024; v1 submitted 4 December, 2023; originally announced December 2023.

Comments: Second version with clarifications on text and figures, and explicit calculations regarding the omega-effect. Manuscript re-submitted to IEEE Transactions On Plasma Science. 4 Figures

Cooling and Instabilities in Colliding Radiative Flows with Toroidal Magnetic Fields

Authors: R. N. Markwick, A. Frank, E. G. Blackman, J. Carroll-Nellenback, S. V. Lebedev, D. R. Russell, J. W. D. Halliday, L. G. Suttle, P. M. Hartigan

Abstract: We report on the results of a simulation based study of colliding magnetized plasma flows. Our set-up mimics pulsed power laboratory astrophysical experiments but, with an appropriate frame change, are relevant to astrophysical jets with internal velocity variations. We track the evolution of the interaction region where the two flows collide. Cooling via radiative loses are included in the calcul… ▽ More We report on the results of a simulation based study of colliding magnetized plasma flows. Our set-up mimics pulsed power laboratory astrophysical experiments but, with an appropriate frame change, are relevant to astrophysical jets with internal velocity variations. We track the evolution of the interaction region where the two flows collide. Cooling via radiative loses are included in the calculation. We systematically vary plasma beta ($β_m$) in the flows, the strength of the cooling ($Λ_0$) and the exponent ($α$) of temperature-dependence of the cooling function. We find that for strong magnetic fields a counter-propagating jet called a "spine" is driven by pressure from shocked toroidal fields. The spines eventually become unstable and break apart. We demonstrate how formation and evolution of the spines depends on initial flow parameters and provide a simple analytic model that captures the basic features of the flow. △ Less

Submitted 13 November, 2023; originally announced November 2023.

Comments: 14 pages, 16 figures. Submitted to MNRAS

How negative feedback and the ambient environment limit the influence of recombination in common envelope evolution

Authors: Luke Chamandy, Jonathan Carroll-Nellenback, Eric G. Blackman, Adam Frank, Yisheng Tu, Baowei Liu, Yangyuxin Zou, Jason Nordhaus

Abstract: We perform 3D hydrodynamical simulations to study recombination and ionization during the common envelope (CE) phase of binary evolution, and develop techniques to track the ionic transitions in time and space. We simulate the interaction of a $2\,M_\odot$ red giant branch primary and a $1\,M_\odot$ companion modeled as a particle. We compare a run employing a tabulated equation of state (EOS) tha… ▽ More We perform 3D hydrodynamical simulations to study recombination and ionization during the common envelope (CE) phase of binary evolution, and develop techniques to track the ionic transitions in time and space. We simulate the interaction of a $2\,M_\odot$ red giant branch primary and a $1\,M_\odot$ companion modeled as a particle. We compare a run employing a tabulated equation of state (EOS) that accounts for ionization and recombination, with a run employing an ideal gas EOS. During the first half of the simulations, $\sim15$ per cent more mass is unbound in the tabulated EOS run due to the release of recombination energy, but by simulation end the difference has become negligible. We explain this as being a consequence of (i) the tabulated EOS run experiences a shallower inspiral and hence smaller orbital energy release at late times because recombination energy release expands the envelope and reduces drag, and (ii) collision and mixing between expanding envelope gas, ejecta and circumstellar ambient gas assists in unbinding the envelope, but does so less efficiently in the tabulated EOS run where some of the energy transferred to bound envelope gas is used for ionization. The rate of mass unbinding is approximately constant in the last half of the simulations and the orbital separation steadily decreases at late times. A simple linear extrapolation predicts a CE phase duration of $\sim2\,\mathrm{yr}$, after which the envelope would be unbound. △ Less

Submitted 4 January, 2024; v1 submitted 28 April, 2023; originally announced April 2023.

Comments: 20 pages, 18 figures, 2 tables. Accepted for publication in MNRAS

Why the observed spin evolution of older-than-solar like stars might not require a dynamo mode change

Authors: Ketevan Kotorashvili, Eric G. Blackman, James E. Owen

Abstract: The spin evolution of main sequence stars has long been of interest for basic stellar evolution, stellar aging, stellar activity, and consequent influence on companion planets. Observations of older than solar late-type main-sequence stars have been interpreted to imply that a change from a dipole-dominated magnetic field to one with more prominent higher multipoles might be necessary to account f… ▽ More The spin evolution of main sequence stars has long been of interest for basic stellar evolution, stellar aging, stellar activity, and consequent influence on companion planets. Observations of older than solar late-type main-sequence stars have been interpreted to imply that a change from a dipole-dominated magnetic field to one with more prominent higher multipoles might be necessary to account for the data. The spin-down models that lead to this inference are essentially tuned to the sun. Here we take a different approach which considers individual stars as fixed points rather than just the Sun. We use a time-dependent theoretical model to solve for the spin evolution of low-mass main-sequence stars that includes a Parker-type wind and a time-evolving magnetic field coupled to the spin. Because the wind is exponentially sensitive to the stellar mass over radius and the coronal base temperature, the use of each observed star as a separate fixed point is more appropriate and, in turn, produces a set of solution curves that produces a solution envelope rather than a simple line. This envelope of solution curves, unlike a single line fit, is consistent with the data and does not unambiguously require a modal transition in the magnetic field to explain it. △ Less

Submitted 4 April, 2023; v1 submitted 11 January, 2023; originally announced January 2023.

Comments: 8 pages, 4 figures; Accepted for publication in MNRAS

The messy death of a multiple star system and the resulting planetary nebula as observed by JWST

Authors: Orsola De Marco, Muhammad Akashi, Stavros Akras, Javier Alcolea, Isabel Aleman, Philippe Amram, Bruce Balick, Elvire De Beck, Eric G. Blackman, Henri M. J. Boffin, Panos Boumis, Jesse Bublitz, Beatrice Bucciarelli, Valentin Bujarrabal, Jan Cami, Nicholas Chornay, You-Hua Chu, Romano L. M. Corradi, Adam Frank, Guillermo Garcia-Segura, D. A. Garcia-Hernandez, Jorge Garcia-Rojas, Veronica Gomez-Llanos, Denise R. Goncalves, Martin A. Guerrero , et al. (44 additional authors not shown)

Abstract: Planetary nebulae (PNe), the ejected envelopes of red giant stars, provide us with a history of the last, mass-losing phases of 90 percent of stars initially more massive than the Sun. Here, we analyse James Webb Space Telescope (JWST) Early Release Observation (ERO) images of the PN NGC3132. A structured, extended H2 halo surrounding an ionised central bubble is imprinted with spiral structures,… ▽ More Planetary nebulae (PNe), the ejected envelopes of red giant stars, provide us with a history of the last, mass-losing phases of 90 percent of stars initially more massive than the Sun. Here, we analyse James Webb Space Telescope (JWST) Early Release Observation (ERO) images of the PN NGC3132. A structured, extended H2 halo surrounding an ionised central bubble is imprinted with spiral structures, likely shaped by a low-mass companion orbiting the central star at 40-60 AU. The images also reveal a mid-IR excess at the central star interpreted as a dusty disk, indicative of an interaction with another, closer companion. Including the previously known, A-type visual companion, the progenitor of the NGC3132 PN must have been at least a stellar quartet. The JWST images allow us to generate a model of the illumination, ionisation and hydrodynamics of the molecular halo, demonstrating the power of JWST to investigate complex stellar outflows. Further, new measurements of the A-type visual companion allow us to derive the value for the mass of the progenitor of a central star to date with excellent precision: 2.86+/-0.06 Mo. These results serve as path finders for future JWST observations of PNe providing unique insight into fundamental astrophysical processes including colliding winds, and binary star interactions, with implications for supernovae and gravitational wave systems. △ Less

Submitted 6 January, 2023; originally announced January 2023.

Comments: 32 pages, 5 figures for the main article. 12 pages 8 figures for the supplementary material

Journal ref: Nature Astronomy, 2022, Vol. 6, p. 1421

Morphology of Shocked Lateral Outflows in Colliding Hydrodynamic Flows

Authors: R. N. Markwick, A. Frank, J. Carroll-Nellenback, E. G. Blackman, P. M. Hartigan, S. V. Lebedev, D. R. Russel, J. W. D. Halliday, L. G. Suttle

Abstract: Supersonic interacting flows occurring in phenomena such as protostellar jets give rise to strong shocks, and have been demonstrated in several laboratory experiments. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in three dimensions. We introduce variations in the flow parameters of density, velocity, and cross sectional radius of the colliding f… ▽ More Supersonic interacting flows occurring in phenomena such as protostellar jets give rise to strong shocks, and have been demonstrated in several laboratory experiments. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in three dimensions. We introduce variations in the flow parameters of density, velocity, and cross sectional radius of the colliding flows %radius in order to study the propagation and conical shape of the bow shock formed by collisions between two, not necessarily symmetric, hypersonic flows. We find that the motion of the interaction region is driven by imbalances in ram pressure between the two flows, while the conical structure of the bow shock is a result of shocked lateral outflows (SLOs) being deflected from the horizontal when the flows are of differing cross-section. △ Less

Submitted 11 December, 2022; originally announced December 2022.

Journal ref: Physics of Plasmas, Volume 29, Issue 10, id.102901, 14 pp (October 2022)

Ion and Electron Acoustic Bursts during Anti-Parallel Magnetic Reconnection Driven by Lasers

Authors: Shu Zhang, Abraham Chien, Lan Gao, Hantao Ji, Eric G. Blackman, Russ Follett, Dustin H. Froula, Joseph Katz, Chikang Li, Andrew Birkel, Richard Petrasso, John Moody, Hui Chen

Abstract: Magnetic reconnection converts magnetic energy into thermal and kinetic energy in plasma. Among numerous candidate mechanisms, ion acoustic instabilities driven by the relative drift between ions and electrons, or equivalently electric current, have been suggested to play a critical role in dissipating magnetic energy in collisionless plasmas. However, their existence and effectiveness during reco… ▽ More Magnetic reconnection converts magnetic energy into thermal and kinetic energy in plasma. Among numerous candidate mechanisms, ion acoustic instabilities driven by the relative drift between ions and electrons, or equivalently electric current, have been suggested to play a critical role in dissipating magnetic energy in collisionless plasmas. However, their existence and effectiveness during reconnection have not been well understood due to ion Landau damping and difficulties in resolving the Debye length scale in the laboratory. Here we report a sudden onset of ion acoustic bursts measured by collective Thomson scattering in the exhaust of anti-parallel magnetically driven reconnection using high-power lasers. The ion acoustic bursts are followed by electron acoustic bursts with electron heating and bulk acceleration. We reproduce these observations with 1D and 2D particle-in-cell simulations in which electron outflow jet drives ion-acoustic instabilities, forming double layers. These layers induce electron two-stream instabilities that generate electron acoustic bursts and energize electrons. Our results demonstrate the importance of ion and electron acoustic dynamics during reconnection when ion Landau damping is ineffective, a condition applicable to a range of astrophysical plasmas including near-Earth space, stellar flares, and black hole accretion engines. △ Less

Submitted 29 March, 2023; v1 submitted 26 September, 2022; originally announced September 2022.

Pre-planetary nebulae: a context for principles, progress, and questions on how binaries and magnetic fields produce jets

Authors: Eric G. Blackman

Abstract: Astrophysical outflows treated initially as spherically symmetric often show evidence for asymmetry once seen at higher resolution. The preponderance of aspherical and multipolar planetary nebulae (PN) and pre-planetary nebulae (PPN) was evident after many observations from the Hubble Space Telescope. Binary interactions have long been thought to be essential for shaping asymmetric PN/PPN, but how… ▽ More Astrophysical outflows treated initially as spherically symmetric often show evidence for asymmetry once seen at higher resolution. The preponderance of aspherical and multipolar planetary nebulae (PN) and pre-planetary nebulae (PPN) was evident after many observations from the Hubble Space Telescope. Binary interactions have long been thought to be essential for shaping asymmetric PN/PPN, but how? PPN are the more kinematically demanding of the two, and warrant particular focus. I address how progress from observation and theory suggests two broad classes of accretion driven PPN jets: one for wider binaries (PPN-W) where the companion is outside the outer radius of the giant and accretes via Roche lobe overflow, and the other which occurs in the later stages of common envelope evolution (CEE) for close binaries (PPN-C). The physics within these scenarios connects to progress and open questions about the role and origin of magnetic fields in the engines and in astrophysical jets more generally. △ Less

Submitted 28 February, 2022; v1 submitted 15 February, 2022; originally announced February 2022.

Comments: 20 Pages, 4 figures; invited review, accepted for the Proceedings of IAU Symposium 366: "The Origin of Outflows in Evolved Stars"; L. Decin, A.A. Zÿlstra and C. Gielen, editors

Jets from main sequence and white dwarf companions during common envelope evolution

Authors: Yangyuxin Zou, Luke Chamandy, Jonathan Carroll-Nellenback, Eric G. Blackman, Adam Frank

Abstract: It has long been speculated that jet feedback from accretion onto the companion during a common envelope (CE) event could affect the orbital evolution and envelope unbinding process, but this conjecture has heretofore remained largely untested. We present global 3D hydrodynamical simulations of CE evolution (CEE) that include a jet subgrid model and compare them with an otherwise identical model w… ▽ More It has long been speculated that jet feedback from accretion onto the companion during a common envelope (CE) event could affect the orbital evolution and envelope unbinding process, but this conjecture has heretofore remained largely untested. We present global 3D hydrodynamical simulations of CE evolution (CEE) that include a jet subgrid model and compare them with an otherwise identical model without a jet. Our binary consists of a $2M_\odot$ red giant branch primary and a $1M_\odot$ or $0.5M_\odot$ main sequence or white dwarf secondary companion modeled as a point particle. We run the simulations for 10 orbits (40 days). Our jet model adds mass at a constant rate $\dot{M}_\mathrm{j}$ of order the Eddington rate, with maximum velocity $v_\mathrm{j}$ of order the escape speed, to two spherical sectors with the jet axis perpendicular to the orbital plane, and supplies kinetic energy at the rate $\sim\dot{M}_\mathrm{j} v_\mathrm{j}^2/40$. We explore the influence of the jet on orbital evolution, envelope morphology and envelope unbinding, and assess the dependence of the results on jet mass-loss rate, launch speed, companion mass, opening angle, and whether or not subgrid accretion is turned on. In line with our theoretical estimates, we find that in all cases the jet becomes choked around the time of first periastron passage. We also find that jets lead to increases in unbound mass of up to $\sim10\%$, as compared to simulations which do not include a jet. △ Less

Submitted 7 September, 2022; v1 submitted 11 February, 2022; originally announced February 2022.

Comments: 16 pages, 8 figures

Characterization of Quasi-Keplerian, Differentially Rotating, Free-Boundary Laboratory Plasmas

Authors: V. Valenzuela-Villaseca, L. G. Suttle, F. Suzuki-Vidal, J. W. D. Halliday, S. Merlini, D. R. Russell, E. R. Tubman, J. D. Hare, J. P. Chittenden, M. E. Koepke, E. G. Blackman, S. V. Lebedev

Abstract: We present results from pulsed-power driven differentially rotating plasma experiments designed to simulate physics relevant to astrophysical disks and jets. In these experiments, angular momentum is injected by the ram pressure of the ablation flows from a wire array Z pinch. In contrast to previous liquid metal and plasma experiments, rotation is not driven by boundary forces. Axial pressure gra… ▽ More We present results from pulsed-power driven differentially rotating plasma experiments designed to simulate physics relevant to astrophysical disks and jets. In these experiments, angular momentum is injected by the ram pressure of the ablation flows from a wire array Z pinch. In contrast to previous liquid metal and plasma experiments, rotation is not driven by boundary forces. Axial pressure gradients launch a rotating plasma jet upwards, which is confined by a combination of ram, thermal, and magnetic pressure of a surrounding plasma halo. The jet has subsonic rotation, with a maximum rotation velocity $23 \pm 3$ km/s. The rotational velocity profile is quasi-Keplerian with a positive Rayleigh discriminant $κ^2 \propto r^{-2.8\pm0.8}$ rad$^2$/s$^2$. The plasma completes $0.5 - 2$ full rotations in the experimental time frame ($\sim 150$ ns). △ Less

Submitted 25 January, 2022; originally announced January 2022.

Comments: 4 figures

Direct measurement of non-thermal electron acceleration from magnetically driven reconnection in a laboratory plasma

Authors: Abraham Chien, Lan Gao, Shu Zhang, Hantao Ji, Eric G. Blackman, William Daughton, Adam Stanier, Ari Le, Fan Guo, Russ Follett, Hui Chen, Gennady Fiksel, Gabriel Bleotu, Robert C. Cauble, Sophia N. Chen, Alice Fazzini, Kirk Flippo, Omar French, Dustin H. Froula, Julien Fuchs, Shinsuke Fujioka, Kenneth Hill, Sallee Klein, Carolyn Kuranz, Philip Nilson , et al. (2 additional authors not shown)

Abstract: Magnetic reconnection is a ubiquitous astrophysical process that rapidly converts magnetic energy into some combination of plasma flow energy, thermal energy, and non-thermal energetic particles, including energetic electrons. Various reconnection acceleration mechanisms in different low-$β$ (plasma-to-magnetic pressure ratio) and collisionless environments have been proposed theoretically and stu… ▽ More Magnetic reconnection is a ubiquitous astrophysical process that rapidly converts magnetic energy into some combination of plasma flow energy, thermal energy, and non-thermal energetic particles, including energetic electrons. Various reconnection acceleration mechanisms in different low-$β$ (plasma-to-magnetic pressure ratio) and collisionless environments have been proposed theoretically and studied numerically, including first- and second-order Fermi acceleration, betatron acceleration, parallel electric field acceleration along magnetic fields, and direct acceleration by the reconnection electric field. However, none of them have been heretofore confirmed experimentally, as the direct observation of non-thermal particle acceleration in laboratory experiments has been difficult due to short Debye lengths for \textit{in-situ} measurements and short mean free paths for \textit{ex-situ} measurements. Here we report the direct measurement of accelerated non-thermal electrons from low-$β$ magnetically driven reconnection in experiments using a laser-powered capacitor coil platform. We use kiloJoule lasers to drive parallel currents to reconnect MegaGauss-level magnetic fields in a quasi-axisymmetric geometry. The angular dependence of the measured electron energy spectrum and the resulting accelerated energies, supported by particle-in-cell simulations, indicate that the mechanism of direct electric field acceleration by the out-of-plane reconnection electric field is at work. Scaled energies using this mechanism show direct relevance to astrophysical observations. Our results therefore validate one of the proposed acceleration mechanisms by reconnection, and establish a new approach to study reconnection particle acceleration with laboratory experiments in relevant regimes. △ Less

Submitted 24 January, 2022; originally announced January 2022.

Comments: 13 pages, 6 figures

Effects of Charge Exchange on the Evaporative Wind of HD 209458b

Authors: Alex Debrecht, Jonathan Carroll-Nellenback, Adam Frank, Eric G. Blackman, Luca Fossati, Ruth Murray-Clay, John McCann

Abstract: The role of charge exchange in shaping exoplanet photoevaporation remains a topic of contention. Exchange of electrons between stellar wind protons from the exoplanet's host star and neutral hydrogen from the planet's wind has been proposed as a mechanism to create "energetic neutral atoms" (ENAs), which could explain the high absorption line velocities observed in systems where mass loss is occur… ▽ More The role of charge exchange in shaping exoplanet photoevaporation remains a topic of contention. Exchange of electrons between stellar wind protons from the exoplanet's host star and neutral hydrogen from the planet's wind has been proposed as a mechanism to create "energetic neutral atoms" (ENAs), which could explain the high absorption line velocities observed in systems where mass loss is occurring. In this paper we present results from 3D hydrodynamic simulations of the mass loss of a planet similar to HD 209458b. We self-consistently launch a planetary wind by calculating the ionization and heating resulting from incident high-energy radiation, inject a stellar wind into the simulation, and allow electron exchange between the stellar and planetary winds. We predict the potential production of ENAs by the wind-wind interaction analytically, then present the results of our simulations, which confirm the analytic limits. Within the limits of our hydrodynamic simulation, we find that charge exchange with the stellar wind properties examined here is unable to explain the absorption observed at high Doppler velocities. △ Less

Submitted 18 January, 2022; originally announced January 2022.

Comments: arXiv admin note: text overlap with arXiv:1906.00075

Cooling and Instabilities in Colliding Flows

Authors: R. N. Markwick, A. Frank, J. Carroll-Nellenback, B. Liu, E. G. Blackman, S. V. Lebedev, P. M. Hartigan

Abstract: Collisional self-interactions occurring in protostellar jets give rise to strong shocks, the structure of which can be affected by radiative cooling within the flow. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in both one and three dimensions with a power law cooling function. The characteristic length and time scales for cooling are temperature… ▽ More Collisional self-interactions occurring in protostellar jets give rise to strong shocks, the structure of which can be affected by radiative cooling within the flow. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in both one and three dimensions with a power law cooling function. The characteristic length and time scales for cooling are temperature dependent and thus may vary as shocked gas cools. When the cooling length decreases sufficiently rapidly the system becomes unstable to the radiative shock instability, which produces oscillations in the position of the shock front; these oscillations can be seen in both the one and three dimensional cases. Our simulations show no evidence of the density clumping characteristic of a thermal instability, even when the cooling function meets the expected criteria. In the three-dimensional case, the nonlinear thin shell instability (NTSI) is found to dominate when the cooling length is sufficiently small. When the flows are subjected to the radiative shock instability, oscillations in the size of the cooling region allow NTSI to occur at larger cooling lengths, though larger cooling lengths delay the onset of NTSI by increasing the oscillation period. △ Less

Submitted 7 September, 2021; originally announced September 2021.

Comments: 12 pages, 21 figures. Accepted for publication in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 508, Issue 2, pp.2266-2278 (Dec 2021)

Scaling of Turbulent Viscosity and Resistivity: Extracting a Scale-dependent Turbulent Magnetic Prandtl Number

Authors: Xin Bian, Jessica K. Shang, Eric G. Blackman, Gilbert W. Collins, Hussein Aluie

Abstract: Turbulent viscosity $ν_t$ and resistivity $η_t$ are perhaps the simplest models for turbulent transport of angular momentum and magnetic fields, respectively. The associated turbulent magnetic Prandtl number $Pr_t\equiv ν_t/η_t$ has been well recognized to determine the final magnetic configuration of accretion disks. Here, we present an approach to determining these ''effective transport'' coeffi… ▽ More Turbulent viscosity $ν_t$ and resistivity $η_t$ are perhaps the simplest models for turbulent transport of angular momentum and magnetic fields, respectively. The associated turbulent magnetic Prandtl number $Pr_t\equiv ν_t/η_t$ has been well recognized to determine the final magnetic configuration of accretion disks. Here, we present an approach to determining these ''effective transport'' coefficients acting at different length-scales using coarse-graining and recent results on decoupled kinetic and magnetic energy cascades [Bian & Aluie 2019]. By analyzing the kinetic and magnetic energy cascades from a suite of high-resolution simulations, we show that our definitions of $ν_t$, $η_t$, and $Pr_t$ have power-law scalings in the ''decoupled range.'' We observe that $Pr_t\approx1 \text{~to~}2$ at the smallest inertial-inductive scales, increasing to $\approx 5$ at the largest scales. However, based on physical considerations, our analysis suggests that $Pr_t$ has to become scale-independent and of order unity in the decoupled range at sufficiently high Reynolds numbers (or grid-resolution), and that the power-law scaling exponents of velocity and magnetic spectra become equal. In addition to implications to astrophysical systems, the scale-dependent turbulent transport coefficients offer a guide for large eddy simulation modeling. △ Less

Submitted 2 July, 2021; originally announced July 2021.

Asteroid Magnetization from the Early Solar Wind

Authors: Atma Anand, Jonathan Carroll-Nellenback, Eric G. Blackman, John A. Tarduno

Abstract: Magnetic fields provide an important probe of the thermal, material, and structural history of planetary and sub-planetary bodies. Core dynamos are a potential source of magnetic fields for differentiated bodies, but evidence of magnetization in undifferentiated bodies requires a different mechanism. Here we study the amplified field provided by the stellar wind to an initially unmagnetized body u… ▽ More Magnetic fields provide an important probe of the thermal, material, and structural history of planetary and sub-planetary bodies. Core dynamos are a potential source of magnetic fields for differentiated bodies, but evidence of magnetization in undifferentiated bodies requires a different mechanism. Here we study the amplified field provided by the stellar wind to an initially unmagnetized body using analytic theory and numerical simulations, employing the resistive MHD AstroBEAR adaptive mesh refinement (AMR) multiphysics code. We obtain a broadly applicable scaling relation for the peak magnetization achieved once a wind advects, piles-up, and drapes a body with magnetic field, reaching a quasi-steady state. We find that the dayside magnetic field for a sufficiently conductive body saturates when it balances the sum of incoming solar wind ram, magnetic, and thermal pressures. Stronger amplification results from pileup by denser and faster winds. Careful quantification of numerical diffusivity is required for accurately interpreting the peak magnetic field strength from simulations, and corroborating with theory. As specifically applied to the Solar System, we find that early solar wind-induced field amplification is a viable source of magnetization for observed paleointensities in meteorites from some undifferentiated bodies. This mechanism may also be applicable to other Solar System bodies, including metal-rich bodies to be visited in future space missions such as the asteroid (16) Psyche. △ Less

Submitted 21 October, 2021; v1 submitted 18 May, 2021; originally announced May 2021.

Comments: 12 pages, 6 figures, accepted by MNRAS

The Formation of Discs in the Interior of AGB Stars from the Tidal Disruption of Planets and Brown Dwarfs

Authors: Gabriel Guidarelli, Jason Nordhaus, Jonathan Carroll-Nellenback, Luke Chamandy, Eric G. Blackman, Adam Frank

Abstract: A significant fraction of isolated white dwarfs host magnetic fields in excess of a MegaGauss. Observations suggest that these fields originate in interacting binary systems where the companion is destroyed thus leaving a singular, highly-magnetized white dwarf. In post-main-sequence evolution, radial expansion of the parent star may cause orbiting companions to become engulfed. During the common… ▽ More A significant fraction of isolated white dwarfs host magnetic fields in excess of a MegaGauss. Observations suggest that these fields originate in interacting binary systems where the companion is destroyed thus leaving a singular, highly-magnetized white dwarf. In post-main-sequence evolution, radial expansion of the parent star may cause orbiting companions to become engulfed. During the common envelope phase, as the orbital separation rapidly decreases, low-mass companions will tidally disrupt as they approach the giant's core. We hydrodynamically simulate the tidal disruption of planets and brown dwarfs, and the subsequent accretion disc formation, in the interior of an asymptotic giant branch star. These dynamically formed discs are commensurate with previous estimates, suggesting strong magnetic fields may originate from these tidal disruption events. △ Less

Submitted 30 April, 2021; originally announced May 2021.

Comments: Submitted to MNRAS

On the shear-current effect: toward understanding why theories and simulations have mutually and separately conflicted

Authors: Hongzhe Zhou, Eric G. Blackman

Abstract: The shear-current effect (SCE) of mean-field dynamo theory refers to the combination of a shear flow and a turbulent coefficient $β_{21}$ with a favorable negative sign for exponential mean-field growth, rather than positive for diffusion. There have been long standing disagreements among theoretical calculations and comparisons of theory with numerical experiments as to the sign of kinetic (… ▽ More The shear-current effect (SCE) of mean-field dynamo theory refers to the combination of a shear flow and a turbulent coefficient $β_{21}$ with a favorable negative sign for exponential mean-field growth, rather than positive for diffusion. There have been long standing disagreements among theoretical calculations and comparisons of theory with numerical experiments as to the sign of kinetic ($β^u_{21}$) and magnetic ($β^b_{21}$) contributions. To resolve these discrepancies, we combine an analytical approach with simulations, and show that unlike $β^b_{21}$, the kinetic SCE $β^u_{21}$ has a strong dependence on the kinetic energy spectral index and can transit from positive to negative values at $\mathcal{O}(10)$ Reynolds numbers if the spectrum is not too steep. Conversely, $β^b_{21}$ is always negative regardless of the spectral index and Reynolds numbers. For very steep energy spectra, the positive $β^u_{21}$ can dominate even at energy equipartition $u_\text{rms}\simeq b_\text{rms}$, resulting in a positive total $β_{21}$ even though $β^b_{21}<0$. Our findings bridge the gap between the seemingly contradictory results from the second-order-correlation approximation (SOCA) versus the spectral-$τ$ closure (STC), for which opposite signs for $β^u_{21}$ have been reported, with the same sign for $β^b_{21}<0$. The results also offer an explanation for the simulations that find $β^u_{21}>0$ and an inconclusive overall sign of $β_{21}$ for $\mathcal{O}(10)$ Reynolds numbers. The transient behavior of $β^u_{21}$ is demonstrated using the kinematic test-field method. We compute dynamo growth rates for cases with or without rotation, and discuss opportunities for further work. △ Less

Submitted 13 September, 2021; v1 submitted 22 April, 2021; originally announced April 2021.

Comments: 17 pages, 12 figures; accepted by MNRAS

Influence of inhomogeneous stochasticity on the falsifiability of mean-field theories and examples from accretion disc modeling

Authors: Hongzhe Zhou, Eric G. Blackman

Abstract: Despite spatial and temporal fluctuations in turbulent astrophysical systems, mean-field theories can be used to describe their secular evolution. However, observations taken over time scales much shorter than dynamical time scales capture a system in a single state of its turbulence ensemble. Comparing with mean-field theory can falsify the latter only if the theory is additionally supplied with… ▽ More Despite spatial and temporal fluctuations in turbulent astrophysical systems, mean-field theories can be used to describe their secular evolution. However, observations taken over time scales much shorter than dynamical time scales capture a system in a single state of its turbulence ensemble. Comparing with mean-field theory can falsify the latter only if the theory is additionally supplied with a quantified precision. The central limit theorem provides appropriate estimates to the precision only when fluctuations contribute linearly to an observable and with constant coherent scales. Here we introduce an error propagation formula that relaxes both limitations, allowing for nonlinear functional forms of observables and inhomogeneous coherent scales and amplitudes of fluctuations. The method is exemplified in the context of accretion disc theories, where inhomogeneous fluctuations in the surface temperature are propagated to the disc emission spectrum--the latter being a nonlinear and non-local function of the former. The derived precision depends non-monotonically on emission frequency. Using the same method, we investigate how binned spectral fluctuations in telescope data change with the spectral resolving power. We discuss the broader implications for falsifiability of a mean-field theory. △ Less

Submitted 30 August, 2021; v1 submitted 12 December, 2020; originally announced December 2020.

Comments: 10 pages, 4 figures; accepted by MNRAS

Successive common envelope events from multiple planets

Authors: Luke Chamandy, Eric G. Blackman, Jason Nordhaus, Emily Wilson

Abstract: Many stars harbour multi-planet systems. As these stars expand late in their evolutions, the innermost planet may be engulfed, leading to a common envelope (CE) event. Even if this is insufficient to eject the envelope, it may expand the star further, causing additional CE events, with the last one unbinding what remains of the envelope. This multi-planet CE scenario may have broad implications fo… ▽ More Many stars harbour multi-planet systems. As these stars expand late in their evolutions, the innermost planet may be engulfed, leading to a common envelope (CE) event. Even if this is insufficient to eject the envelope, it may expand the star further, causing additional CE events, with the last one unbinding what remains of the envelope. This multi-planet CE scenario may have broad implications for stellar and planetary evolution across a range of systems. We develop a simplified version and show that it may be able to explain the recently observed planet WD 1856 b. △ Less

Submitted 18 February, 2021; v1 submitted 22 November, 2020; originally announced November 2020.

Comments: 5 pages, 2 figures, accepted for publication in MNRAS Letters

Major Scientific Challenges and Opportunities in Understanding Magnetic Reconnection and Related Explosive Phenomena in Solar and Heliospheric Plasmas

Authors: H. Ji, J. Karpen, A. Alt, S. Antiochos, S. Baalrud, S. Bale, P. M. Bellan, M. Begelman, A. Beresnyak, A. Bhattacharjee, E. G. Blackman, D. Brennan, M. Brown, J. Buechner, J. Burch, P. Cassak, B. Chen, L. -J. Chen, Y. Chen, A. Chien, L. Comisso, D. Craig, J. Dahlin, W. Daughton, E. DeLuca , et al. (83 additional authors not shown)

Abstract: Magnetic reconnection underlies many explosive phenomena in the heliosphere and in laboratory plasmas. The new research capabilities in theory/simulations, observations, and laboratory experiments provide the opportunity to solve the grand scientific challenges summarized in this whitepaper. Success will require enhanced and sustained investments from relevant funding agencies, increased interagen… ▽ More Magnetic reconnection underlies many explosive phenomena in the heliosphere and in laboratory plasmas. The new research capabilities in theory/simulations, observations, and laboratory experiments provide the opportunity to solve the grand scientific challenges summarized in this whitepaper. Success will require enhanced and sustained investments from relevant funding agencies, increased interagency/international partnerships, and close collaborations of the solar, heliospheric, and laboratory plasma communities. These investments will deliver transformative progress in understanding magnetic reconnection and related explosive phenomena including space weather events. △ Less

Submitted 16 September, 2020; originally announced September 2020.

Comments: 4 pages (including a title page), white paper submitted to Helio2050 workshop at https://www.hou.usra.edu/meetings/helio2050/. arXiv admin note: substantial text overlap with arXiv:2004.00079

Persistent mysteries of jet engines, formation, propagation, and particle acceleration: have they been addressed experimentally?

Authors: Eric G. Blackman, Sergey V. Lebedev

Abstract: The physics of astrophysical jets can be divided into three regimes: (i) engine and launch (ii) propagation and collimation, (iii) dissipation and particle acceleration. Since astrophysical jets comprise a huge range of scales and phenomena, practicality dictates that most studies of jets intentionally or inadvertently focus on one of these regimes, and even therein, one body of work may be simply… ▽ More The physics of astrophysical jets can be divided into three regimes: (i) engine and launch (ii) propagation and collimation, (iii) dissipation and particle acceleration. Since astrophysical jets comprise a huge range of scales and phenomena, practicality dictates that most studies of jets intentionally or inadvertently focus on one of these regimes, and even therein, one body of work may be simply boundary condition for another. We first discuss long standing persistent mysteries that pertain the physics of each of these regimes, independent of the method used to study them. This discussion makes contact with frontiers of plasma astrophysics more generally. While observations theory, and simulations, and have long been the main tools of the trade, what about laboratory experiments? Jet related experiments have offered controlled studies of specific principles, physical processes, and benchmarks for numerical and theoretical calculations. We discuss what has been done to date on these fronts. Although experiments have indeed helped us to understand certain processes, proof of principle concepts, and benchmarked codes, they have yet to solved an astrophysical jet mystery on their own. A challenge is that experimental tools used for jet-related experiments so far, are typically not machines originally designed for that purpose, or designed with specific astrophysical mysteries in mind. This presents an opportunity for a different way of thinking about the development of future platforms: start with the astrophysical mystery and build an experiment to address it. △ Less

Submitted 17 September, 2020; originally announced September 2020.

Comments: 37 pages, including 11 figures at the end; for New Astronomy Reviews special volume '100 Years of Jets' (eds. Rob Fender and Ralph Wijers)

Necessity of ventilation for mitigating virus transmission quantified simply

Authors: Eric G. Blackman, Gourab Ghoshal

Abstract: To mitigate the SARS-CoV-2 pandemic, officials have employed social distancing and stay-at-home measures, with increased attention to room ventilation emerging only more recently. Effective distancing practices for open spaces can be ineffective for poorly ventilated spaces, both of which are commonly filled with turbulent air. This is typical for indoor spaces that use mixing ventilation. While t… ▽ More To mitigate the SARS-CoV-2 pandemic, officials have employed social distancing and stay-at-home measures, with increased attention to room ventilation emerging only more recently. Effective distancing practices for open spaces can be ineffective for poorly ventilated spaces, both of which are commonly filled with turbulent air. This is typical for indoor spaces that use mixing ventilation. While turbulence initially reduces the risk of infection near a virion-source, it eventually increases the exposure risk for all occupants in a space without ventilation. To complement detailed models aimed at precision, minimalist frameworks are useful to facilitate order of magnitude estimates for how much ventilation provides safety, particularly when circumstances require practical decisions with limited options. Applying basic principles of transport and diffusion, we estimate the time-scale for virions injected into a room of turbulent air to infect an occupant, distinguishing cases of low vs. high initial virion mass loads and virion-destroying vs. virion-reflecting walls. We consider the effect of an open window as a proxy for ventilation. When the airflow is dominated by isotropic turbulence, the minimum area needed to ensure safety depends only on the ratio of total viral load to threshold load for infection. The minimalist estimates here convey simply that the equivalent of ventilation by modest sized open window in classrooms and workplaces significantly improves safety. △ Less

Submitted 18 August, 2020; v1 submitted 20 June, 2020; originally announced June 2020.

Comments: 8 pages, 1 table, 1 figure, submitted

Common Envelope Evolution on the Asymptotic Giant Branch: Unbinding within a Decade?

Authors: Luke Chamandy, Eric G. Blackman, Adam Frank, Jonathan Carroll-Nellenback, Yisheng Tu

Abstract: Common envelope (CE) evolution is a critical but still poorly understood progenitor phase of many high-energy astrophysical phenomena. Although 3D global hydrodynamic CE simulations have become more common in recent years, those involving an asymptotic giant branch (AGB) primary are scarce, due to the high computational cost from the larger dynamical range compared to red giant branch (RGB) primar… ▽ More Common envelope (CE) evolution is a critical but still poorly understood progenitor phase of many high-energy astrophysical phenomena. Although 3D global hydrodynamic CE simulations have become more common in recent years, those involving an asymptotic giant branch (AGB) primary are scarce, due to the high computational cost from the larger dynamical range compared to red giant branch (RGB) primaries. But CE evolution with AGB progenitors is desirable to simulate because such events are the likely progenitors of most bi-polar planetary nebulae (PNe), and prominent observational testing grounds for CE physics. Here we present a high resolution global simulation of CE evolution involving an AGB primary and $1\,\mathrm{M}_\odot$ secondary, evolved for $20$ orbital revolutions. During the last $16$ of these orbits, the envelope unbinds at an almost constant rate of about $0.1$-$0.2\,\mathrm{M}_\odot\,\mathrm{yr}^{-1}$. If this rate were maintained, the envelope would be unbound in less than $10\,\mathrm{yr}$. The dominant source of this unbinding is consistent with inspiral; we assess the influence of the ambient medium to be subdominant. We compare this run with a previous run that used an RGB phase primary evolved from the same $2\,\mathrm{M}_\odot$ main sequence star to assess the influence of the evolutionary state of the primary. When scaled appropriately, the two runs are quite similar, but with some important differences. △ Less

Submitted 27 May, 2020; v1 submitted 14 April, 2020; originally announced April 2020.

Comments: 12 pages, 12 figures, published in MNRAS

Major Scientific Challenges and Opportunities in Understanding Magnetic Reconnection and Related Explosive Phenomena throughout the Universe

Authors: H. Ji, A. Alt, S. Antiochos, S. Baalrud, S. Bale, P. M. Bellan, M. Begelman, A. Beresnyak, E. G. Blackman, D. Brennan, M. Brown, J. Buechner, J. Burch, P. Cassak, L. -J. Chen, Y. Chen, A. Chien, D. Craig, J. Dahlin, W. Daughton, E. DeLuca, C. F. Dong, S. Dorfman, J. Drake, F. Ebrahimi , et al. (75 additional authors not shown)

Abstract: This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process. This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process. △ Less

Submitted 31 March, 2020; originally announced April 2020.

Comments: 6 pages, 1 figure, white paper submitted to both Plasma 2020 and Astro 2020 Decadal Surveys

Improving TBI Protection Measures and Standards for Combat Helmets

Authors: Eric G. Blackman

Abstract: How well do current combat helmets protect against Traumatic Brain Injury? To answer this question it is necessary to evaluate both the measures and standards of protection. We define measure as a physical test that the helmet must be subject to and the standard as the quantitative threshold of performance in this test that the helmet must satisfy to be acceptable. We find that both the measures a… ▽ More How well do current combat helmets protect against Traumatic Brain Injury? To answer this question it is necessary to evaluate both the measures and standards of protection. We define measure as a physical test that the helmet must be subject to and the standard as the quantitative threshold of performance in this test that the helmet must satisfy to be acceptable. We find that both the measures are inadequate and the standards are too low. Current combat helmets simply do not protect against TBI effectively as they have not been designed to do so. We explain this inadequacy of combat helmet protection against TBI, and how the situation can be dramatically improved both in the short and long term. △ Less

Submitted 14 March, 2020; originally announced March 2020.

Comments: Report of archival value from 2006-2007 (small update 2011) to the Defense Science Study Group (institute for Defense Analyses, Alexandria VA); 32 pages

A Turbulent-Entropic Instability and the Fragmentation of Star-Forming Clouds

Authors: Eric Keto, George B. Field, Eric G. Blackman

Abstract: The kinetic energy of supersonic turbulence within interstellar clouds is subject to cooling by dissipation in shocks and subsequent line radiation. The clouds are therefore susceptible to a condensation process controlled by the specific entropy. In a form analogous to the thermodynamic entropy, the entropy for supersonic turbulence is proportional to the log of the product of the mean turbulent… ▽ More The kinetic energy of supersonic turbulence within interstellar clouds is subject to cooling by dissipation in shocks and subsequent line radiation. The clouds are therefore susceptible to a condensation process controlled by the specific entropy. In a form analogous to the thermodynamic entropy, the entropy for supersonic turbulence is proportional to the log of the product of the mean turbulent velocity and the size scale. We derive a dispersion relation for the growth of entropic instabilities in a spherical self-gravitating cloud and find that there is a critical maximum dissipation time scale, about equal to the crossing time, that allows for fragmentation and subsequent star formation. However, the time scale for the loss of turbulent energy may be shorter or longer, for example with rapid thermal cooling or the injection of mechanical energy. Differences in the time scale for energy loss in different star-forming regions may result in differences in the outcome, for example, in the initial mass function. △ Less

Submitted 8 January, 2020; originally announced January 2020.

Comments: To be published in MNRAS

Synthesizing Observations and Theory to Understand Galactic Magnetic Fields: Progress and Challenges

Authors: Rainer Beck, Luke Chamandy, Ed Elson, Eric G. Blackman

Abstract: Constraining dynamo theories of magnetic field origin by observation is indispensable but challenging, in part because the basic quantities measured by observers and predicted by modelers are different. We clarify these differences and sketch out ways to bridge the divide. Based on archival and previously unpublished data, we then compile various important properties of galactic magnetic fields fo… ▽ More Constraining dynamo theories of magnetic field origin by observation is indispensable but challenging, in part because the basic quantities measured by observers and predicted by modelers are different. We clarify these differences and sketch out ways to bridge the divide. Based on archival and previously unpublished data, we then compile various important properties of galactic magnetic fields for nearby spiral galaxies. We consistently compute strengths of total, ordered, and regular fields, pitch angles of ordered and regular fields, and we summarize the present knowledge on azimuthal modes, field parities, and the properties of non-axisymmetric spiral features called magnetic arms. We review related aspects of dynamo theory, with a focus on mean-field models and their predictions for large-scale magnetic fields in galactic discs and halos. Further, we measure the velocity dispersion of HI gas in arm and inter-arm regions in three galaxies, M 51, M 74, and NGC 6496, since spiral modulation of the root-mean-square turbulent speed has been proposed as a driver of non-axisymmetry in large-scale dynamos. We find no evidence for such a modulation and place upper limits on its strength, helping to narrow down the list of mechanisms to explain magnetic arms. Successes and remaining challenges of dynamo models with respect to explaining observations are briefly summarized, and possible strategies are suggested. With new instruments like the Square Kilometre Array (SKA), large data sets of magnetic and non-magnetic properties from thousands of galaxies will become available, to be compared with theory. △ Less

Submitted 21 December, 2019; v1 submitted 18 December, 2019; originally announced December 2019.

Comments: Submitted on 10 October 2019 to "Galaxies", accepted on 16 December 2019

MSC Class: J.2 ACM Class: J.2

Bipolar Planetary Nebulae from Outflow Collimation by Common Envelope Evolution

Authors: Yangyuxin Zou, Adam Frank, Zhuo Chen, Thomas Reichardt, Orsola De Marco, Eric G. Blackman, Jason Nordhaus, Bruce Balick, Jonathan Carroll-Nellenback, Luke Chamandy, Baowei Liu

Abstract: The morphology of bipolar planetary nebulae (PNe) can be attributed to interactions between a fast wind from the central engine and dense toroidal shaped ejecta left over from common envelope (CE) evolution. Here we use the 3-D hydrodynamic AMR code AstroBEAR to study the possibility that bipolar PN outflows can emerge collimated even from an uncollimated spherical wind in the aftermath of a CE ev… ▽ More The morphology of bipolar planetary nebulae (PNe) can be attributed to interactions between a fast wind from the central engine and dense toroidal shaped ejecta left over from common envelope (CE) evolution. Here we use the 3-D hydrodynamic AMR code AstroBEAR to study the possibility that bipolar PN outflows can emerge collimated even from an uncollimated spherical wind in the aftermath of a CE event. The output of a single CE simulation via the SPH code PHANTOM serves as the initial conditions. Four cases of winds, all with high enough momenta to account for observed high momenta preplanetary nebula outflows, are injected spherically from the region of the CE binary remnant into the ejecta. We compare cases with two different momenta and cases with no radiative cooling versus application of optically thin emission via a cooling curve to the outflow. Our simulations show that in all cases highly collimated bipolar outflows result from deflection of the spherical wind via the interaction with the CE ejecta. Significant asymmetries between the top and bottom lobes are seen in all cases. The asymmetry is strongest for the lower momentum case with radiative cooling. While real post CE winds may be aspherical, our models show that collimation via "inertial confinement" will be strong enough to create jet-like outflows even beginning with maximally uncollimated drivers. Our simulations reveal detailed shock structures in the shock focused inertial confinement (SFIC) model and develop a lens-shaped inner shock that is a new feature of SFIC driven bipolar lobes. △ Less

Submitted 3 June, 2020; v1 submitted 3 December, 2019; originally announced December 2019.

Comments: 17 pages, 16 figures

Cloud fragmentation cascades and feedback: on reconciling an unfettered inertial range with a low star formation rate

Authors: Eric G. Blackman

Abstract: Molecular cloud complexes exhibit both (i) an unfettered Larson-type spectrum over much of their dynamic range, whilst (ii) still producing a much lower star-formation rate than were this cascade to remain unfettered all the way down to star-forming scales. Here we explain the compatibility of these attributes with minimalist considerations of a mass-conserving fragmentation cascade, combined with… ▽ More Molecular cloud complexes exhibit both (i) an unfettered Larson-type spectrum over much of their dynamic range, whilst (ii) still producing a much lower star-formation rate than were this cascade to remain unfettered all the way down to star-forming scales. Here we explain the compatibility of these attributes with minimalist considerations of a mass-conserving fragmentation cascade, combined with estimates of stellar feedback. Of importance is that the amount of feedback needed to abate fragmentation and truncate the complex decreases with decreasing scale. The scale at which the feedback momentum matches the free-fall momentum marks a transition scale below most of the cascade is truncated and the molecular cloud complex dissipated. For a $10^6M_\odot$ GMC complex starting with radius of $\sim 50$pc, the combined feedback from young stellar objects, supernovae, radiation, and stellar winds for a GMC cloud complex can truncate the cascade within an outer free-fall time but only after the cascade reaches parsec scales. △ Less

Submitted 4 February, 2020; v1 submitted 18 October, 2019; originally announced October 2019.

Comments: 6 pages, accepted by MNRAS

Determination of Starspot Covering Fraction as a function of Stellar Age from Observational Data

Authors: Fiona Nichols-Fleming, Eric G. Blackman

Abstract: The association of starspots with magnetic fields leads to an expectation that quantities which correlate with magnetic field strength may also correlate with {starspot} coverage. Since younger stars spin faster and are more magnetically active, assessing whether {starspot} coverage correlates with shorter rotation periods and stellar youth tests these principles. Here we analyze the {starspot} co… ▽ More The association of starspots with magnetic fields leads to an expectation that quantities which correlate with magnetic field strength may also correlate with {starspot} coverage. Since younger stars spin faster and are more magnetically active, assessing whether {starspot} coverage correlates with shorter rotation periods and stellar youth tests these principles. Here we analyze the {starspot} covering fraction versus stellar age for M{-}, G{-}, K{-}, and F{-}type stars based on previously determined variability and rotation periods of over 30,000 {\textit{Kepler}} main-sequence stars. We determine the correlation between age and variability using single and dual power law best fits. We find that {starspot} coverage does indeed decrease with age. Only when the data {are} binned in an effort to remove the effects of activity cycles of individual stars, do statistically significant power law fits emerge for each stellar type. {Using bin averages,} we then find that the {starspot} covering fraction scales with the {X}-ray to bolometric ratio to the power $λ$ with {$0.22\pm 0.03 < λ< 0.32\pm 0.09$} for {G-type} stars of rotation period below 15 days and for the full range of F{-} and M{-type} stars. For K{-}type stars, we find two branches of $λ$ separated by variability bins, with the lower branch showing nearly constant starspot coverage and the upper branch {$λ\sim 0.35\pm 0.04$.} G{-}type stars with periods longer than $15$ days exhibit a transition to steeper power law of {$ λ\sim 2.4 \pm 1.0$.} The potential connection to previous rotation-age measurements suggesting a magnetic breaking transition at the solar age, corresponding to period of $24.5$ is also of interest. △ Less

Submitted 13 November, 2019; v1 submitted 6 September, 2019; originally announced September 2019.

Comments: 10 Pages, 7 figures, accepted by MNRAS

How Drag Force Evolves in Global Common Envelope Simulations

Authors: Luke Chamandy, Eric G. Blackman, Adam Frank, Jonathan Carroll-Nellenback, Yangyuxin Zou, Yisheng Tu

Abstract: We compute the forces, torque and rate of work on the companion-core binary due to drag in global simulations of common envelope (CE) evolution for three different companion masses. Our simulations help to delineate regimes when conventional analytic drag force approximations are applicable. During and just prior to the first periastron passage of the in-spiral phase, the drag force is reasonably… ▽ More We compute the forces, torque and rate of work on the companion-core binary due to drag in global simulations of common envelope (CE) evolution for three different companion masses. Our simulations help to delineate regimes when conventional analytic drag force approximations are applicable. During and just prior to the first periastron passage of the in-spiral phase, the drag force is reasonably approximated by conventional analytic theory and peaks at values proportional to the companion mass. Good agreement between global and local 3D "wind tunnel" simulations, including similar net drag force and flow pattern, is obtained for comparable regions of parameter space. However, subsequent to the first periastron passage, the drag force is up to an order of magnitude smaller than theoretical predictions, quasi-steady, and depends only weakly on companion mass. The discrepancy is exacerbated for larger companion mass and when the inter-particle separation reduces to the Bondi-Hoyle-Lyttleton accretion radius, creating a turbulent thermalized region. Greater flow symmetry during this phase leads to near balance of opposing gravitational forces in front of and behind the companion, hence a small net drag. The reduced drag force at late times helps explain why companion-core separations necessary for envelope ejection are not reached by the end of limited duration CE simulations. △ Less

Submitted 3 October, 2019; v1 submitted 16 August, 2019; originally announced August 2019.

Comments: 14 pages, 8 figures, MNRAS in press

Hydrodynamic Simulations of Disrupted Planetary Accretion Discs Inside the Core of an AGB Star

Authors: Gabriel Guidarelli, Jason Nordhaus, Luke Chamandy, Zhuo Chen, Eric G. Blackman, Adam Frank, Jonathan Carroll-Nellenback, Baowei Liu

Abstract: Volume complete sky surveys provide evidence for a binary origin for the formation of isolated white dwarfs with magnetic fields in excess of a MegaGauss. Interestingly, not a single high-field magnetic white dwarf has been found in a detached system suggesting that if the progenitors are indeed binaries, the companion must be removed or merge during formation. An origin scenario consistent with o… ▽ More Volume complete sky surveys provide evidence for a binary origin for the formation of isolated white dwarfs with magnetic fields in excess of a MegaGauss. Interestingly, not a single high-field magnetic white dwarf has been found in a detached system suggesting that if the progenitors are indeed binaries, the companion must be removed or merge during formation. An origin scenario consistent with observations involves the engulfment, inspiral, and subsequent tidal disruption of a low-mass companion in the interior of a giant star during a common envelope phase. Material from the shredded companion forms a cold accretion disc embedded in the hot ambient around the proto-white dwarf. Entrainment of hot material may evaporate the disc before it can sufficiently amplify the magnetic field, which typically requires at least a few orbits of the disc.Using three-dimensional hydrodynamic simulations of accretion discs with masses between 1 and 10 times the mass of Jupiter inside the core of an Asymptotic Giant Branch star, we find that the discs survive for at least 10 orbits (and likely for 100 orbits), sufficient for strong magnetic fields to develop. △ Less

Submitted 31 July, 2019; originally announced August 2019.

Comments: Submitted to MNRAS

Effects of Radiation Pressure on the Evaporative Wind of HD 209458b

Authors: Alex Debrecht, Jonathan Carroll-Nellenback, Adam Frank, Eric G. Blackman, Luca Fossati, John McCann, Ruth Murray-Clay

Abstract: The role of radiation pressure in shaping exoplanet photoevaporation remains a topic of contention. Radiation pressure from the exoplanet's host star has been proposed as a mechanism to drive the escaping atmosphere into a "cometary" tail and explain the high velocities observed in systems where mass loss is occurring. In this paper we present results from high-resolution 3-D hydrodynamic simulati… ▽ More The role of radiation pressure in shaping exoplanet photoevaporation remains a topic of contention. Radiation pressure from the exoplanet's host star has been proposed as a mechanism to drive the escaping atmosphere into a "cometary" tail and explain the high velocities observed in systems where mass loss is occurring. In this paper we present results from high-resolution 3-D hydrodynamic simulations of a planet similar to HD 209458b. We self-consistently launch a wind flowing outward from the planet by calculating the ionization and heating resulting from incident high-energy radiation, and account for radiation pressure. We first present a simplified calculation, setting a limit on the Lyman-$α$ flux required to drive the photo-evaporated planetary material to larger radii and line-of-sight velocities. We then present the results of our simulations, which confirm the limits determined by our analytic calculation. We thus demonstrate that, within the limits of our hydrodynamic simulation and for the Lyman-$α$ fluxes expected for HD 209458, radiation pressure is unlikely to significantly affect photoevaporative winds or to explain the high velocities at which wind material is observed, though further possibilities remain to be investigated. △ Less

Submitted 3 February, 2020; v1 submitted 31 May, 2019; originally announced June 2019.

Generalized quenching of large-scale magnetic dynamos in anisotropic flows

Authors: Hongzhe Zhou, Eric G. Blackman

Abstract: The buildup of small-scale magnetic helicity which accompanies the oppositely signed growth on large scales is central to conventional dynamical quenching theories of mean-field dynamos. However, the conventional formalism presumes isotropic turbulence and thereby excludes part of the magnetic Lorentz back-reaction. This renders it insufficient to predict the full quenching for general anisotropic… ▽ More The buildup of small-scale magnetic helicity which accompanies the oppositely signed growth on large scales is central to conventional dynamical quenching theories of mean-field dynamos. However, the conventional formalism presumes isotropic turbulence and thereby excludes part of the magnetic Lorentz back-reaction. This renders it insufficient to predict the full quenching for general anisotropic flows. To overcome this deficiency, we derive a new generalized quenching formalism that includes the full back-reacting Lorentz force, and a new "selective-damping-$τ$" closure which conserves magnetic helicity. We apply the formalism to examples of $\bmα^2$ dynamos and show its predicted quenching for different cases of turbulence---isotropic helical, anisotropic helical, and anisotropic non-helical. It predicts stronger-than-conventional quenching in general, but reduces to the conventional case in the helical isotropic limit. △ Less

Submitted 9 June, 2020; v1 submitted 3 May, 2019; originally announced May 2019.

Comments: 16 pages, 2 figures; submitted to Journal of Plasma Physics

Solving the Riemann Problem for Realistic Astrophysical Fluids

Authors: Zhuo Chen, Matthew S. B. Coleman, Eric G. Blackman, Adam Frank

Abstract: We present new methods to solve the Riemann problem both exactly and approximately for general equations of state (EoS) to facilitate realistic modeling and understanding of astrophysical flows. The existence and uniqueness of the new exact general EoS Riemann solution can be guaranteed if the EoS is monotone regardless of the physical validity of the EoS. We confirm that: (1) the solution of the… ▽ More We present new methods to solve the Riemann problem both exactly and approximately for general equations of state (EoS) to facilitate realistic modeling and understanding of astrophysical flows. The existence and uniqueness of the new exact general EoS Riemann solution can be guaranteed if the EoS is monotone regardless of the physical validity of the EoS. We confirm that: (1) the solution of the new exact general EoS Riemann solver and the solution of the original exact Riemann solver match when calculating perfect gas Euler equations; (2) the solution of the new Harten-Lax-van Leer-Contact (HLLC) general EoS Riemann solver and the solution of the original HLLC Riemann solver match when working with perfect gas EoS; and (3) the solution of the new HLLC general EoS Riemann solver approaches the new exact solution. We solve the EoS with two methods, one is to interpolate 2D EoS tables by the bi-linear interpolation method, and the other is to analytically calculate thermodynamic variables at run-time. The interpolation method is more general as it can work with other monotone and realistic EoS while the analytic EoS solver introduced here works with a relatively idealized EoS. Numerical results confirm that the accuracy of the two EoS solvers is similar. We study the efficiency of these two methods with the HLLC general EoS Riemann solver and find that analytic EoS solver is faster in the test problems. However, we point out that a combination of the two EoS solvers may become favorable in some specific problems. Throughout this research, we assume local thermal equilibrium. △ Less

Submitted 2 November, 2019; v1 submitted 11 March, 2019; originally announced March 2019.

Comments: 42 pages, 19 figures, accepted by the Journal of Computational Physics. Corrected some grammatical errors and a typo in the latest version

Minimalist large scale dynamo from shear-driven inhomogeneity

Authors: F. Ebrahimi, E. G. Blackman

Abstract: We show that non-axisymmetric, non-helical perturbations in an unstratified shear flow produce a shear-plane averaged electromotive force (EMF) proportional to a spatially dependent kinetic helicity. This new "shear-driven $α$-effect" can amplify shear-plane averaged magnetic fields, but is undiscoverable when invoking the homogeneity assumption for dynamo coefficients. It is a candidate mechanism… ▽ More We show that non-axisymmetric, non-helical perturbations in an unstratified shear flow produce a shear-plane averaged electromotive force (EMF) proportional to a spatially dependent kinetic helicity. This new "shear-driven $α$-effect" can amplify shear-plane averaged magnetic fields, but is undiscoverable when invoking the homogeneity assumption for dynamo coefficients. It is a candidate mechanism underlying large scale field generation in astrophysical, laboratory, and simulated shear flows, and highlights the essentiality of correctly averaging dynamo coefficients and fields. △ Less

Submitted 5 May, 2019; v1 submitted 12 February, 2019; originally announced February 2019.

Comments: 4 pages

Energy Budget and Core-Envelope Motion in Common Envelope Evolution

Authors: Luke Chamandy, Yisheng Tu, Eric G. Blackman, Jonathan Carroll-Nellenback, Adam Frank, Baowei Liu, Jason Nordhaus

Abstract: We analyze a 3D hydrodynamic simulation of common envelope evolution to understand how energy is transferred between various forms and whether theory and simulation are mutually consistent given the setup. Virtually all of the envelope unbinding in the simulation occurs before the end of the rapid plunge-in phase, here defined to coincide with the first periastron passage. In contrast, the total e… ▽ More We analyze a 3D hydrodynamic simulation of common envelope evolution to understand how energy is transferred between various forms and whether theory and simulation are mutually consistent given the setup. Virtually all of the envelope unbinding in the simulation occurs before the end of the rapid plunge-in phase, here defined to coincide with the first periastron passage. In contrast, the total envelope energy is nearly constant during this time because positive energy transferred to the gas from the core particles is counterbalanced by the negative binding energy from the closer proximity of the inner layers to the plunged-in secondary. During the subsequent slow spiral-in phase, energy continues to transfer to the envelope from the red giant core and secondary core particles. We also propose that relative motion between the centre of mass of the envelope and the centre of mass of the particles could account for the offsets of planetary nebula central stars from the nebula's geometric centre. △ Less

Submitted 25 March, 2019; v1 submitted 28 December, 2018; originally announced December 2018.

Comments: 17 pages, 10 figures, 5 tables, accepted for publication in MNRAS

Photoevaporative Flows From Exoplanet Atmospheres: A 3-D Radiative Hydrodynamic Parameter Study

Authors: Alex Debrecht, Jonathan Carroll-Nellenback, Adam Frank, John McCann, Ruth Murray-Clay, Eric G. Blackman

Abstract: The photoionization-driven evaporation of planetary atmospheres has emerged as a potentially fundamental process for planets on short period orbits. While 1-D studies have proven the effectiveness of stellar fluxes at altering the atmospheric mass and composition for sub-Jupiter mass planets, there remains much that is uncertain with regard to the larger-scale, multidimensional nature of such "pla… ▽ More The photoionization-driven evaporation of planetary atmospheres has emerged as a potentially fundamental process for planets on short period orbits. While 1-D studies have proven the effectiveness of stellar fluxes at altering the atmospheric mass and composition for sub-Jupiter mass planets, there remains much that is uncertain with regard to the larger-scale, multidimensional nature of such "planetary wind" flows. In this paper we use a new radiation-hydrodynamic platform to simulate atmospheric evaporative flows. Using the AstroBEAR AMR multiphysics code in a co-rotating frame centered on the planet, we model the transfer of ionizing photons into the atmosphere, the subsequent launch of the wind and the wind's large scale evolution subject to tidal and non-inertial forces. We run simulations for planets of 0.263 and 0.07 Jupiter masses and stellar fluxes of $2 \times 10^{13}$ and $2 \times 10^{14}$ photons/cm^2/s. Our results reveal new, potentially observable planetary wind flow patterns, including the development, in some cases, of an extended neutral tail lagging behind the planet in its orbit. △ Less

Submitted 22 November, 2018; originally announced November 2018.

Calculating turbulent transport tensors by averaging single plume dynamics and application to dynamos

Authors: Hongzhe Zhou, Eric G. Blackman

Abstract: Transport coefficients in turbulence are comprised of correlation functions between turbulent fluctuations and efficient methods to calculate them are desirable. For example, in mean field dynamo theories used to model the growth of large scale magnetic fields of stars and galaxies, the turbulent electromotive force is commonly approximated by a series of tensor products of turbulent transport coe… ▽ More Transport coefficients in turbulence are comprised of correlation functions between turbulent fluctuations and efficient methods to calculate them are desirable. For example, in mean field dynamo theories used to model the growth of large scale magnetic fields of stars and galaxies, the turbulent electromotive force is commonly approximated by a series of tensor products of turbulent transport coefficients with successively higher order spatial derivatives of the mean magnetic field. One ingredient of standard models is the kinematic coefficient of the zeroth order term, namely the averaged kinetic pseudotensor $\bmα$, that converts toroidal to poloidal fields. Here we demonstrate an efficient way to calculate this quantity for rotating stratified turbulence, whereby the pre-averaged quantity is calculated for the motion of a single plume, and the average is then taken over an ensemble of plumes of different orientations. We calculate the plume dynamics in the most convenient frame, before transforming back to the lab frame and averaging. Our concise configuration space calculation gives essentially identical results to previous lengthier approaches. The present application exemplifies what is a broadly applicable method. △ Less

Submitted 28 November, 2018; v1 submitted 16 October, 2018; originally announced October 2018.

Comments: 6 pages; 1 figure; accepted by MNRAS Letters

Accretion in Common Envelope Evolution

Authors: Luke Chamandy, Adam Frank, Eric G. Blackman, Jonathan Carroll-Nellenback, Baowei Liu, Yisheng Tu, Jason Nordhaus, Zhuo Chen, Bo Peng

Abstract: Common envelope evolution (CEE) occurs in some binary systems involving asymptotic giant branch (AGB) or red giant branch (RGB) stars, and understanding this process is crucial for understanding the origins of various transient phenomena. CEE has been shown to be highly asymmetrical and global 3D simulations are needed to help understand the dynamics. We perform and analyze hydrodynamic CEE simula… ▽ More Common envelope evolution (CEE) occurs in some binary systems involving asymptotic giant branch (AGB) or red giant branch (RGB) stars, and understanding this process is crucial for understanding the origins of various transient phenomena. CEE has been shown to be highly asymmetrical and global 3D simulations are needed to help understand the dynamics. We perform and analyze hydrodynamic CEE simulations with the adaptive mesh refinement (AMR) code AstroBEAR, and focus on the role of accretion onto the companion star. We bracket the range of accretion rates by comparing a model that removes mass and pressure using a subgrid accretion prescription with one that does not. Provided a pressure-release valve, such as a bipolar jet, is available, super-Eddington accretion could be common. Finally, we summarize new results pertaining to the energy budget, and discuss the overall implications relating to the feasibility of unbinding the envelope in CEE simulations. △ Less

Submitted 10 October, 2018; originally announced October 2018.

Comments: 4 pages, 3 figures, to appear in the proceedings of IAU Symposium 343: Why Galaxies Care About AGB Stars, (eds.) F. Kerschbaum, M. Groenewegen and H. Olofsson

Journal ref: Proc. IAU 14 (2018) 235-238

Accretion in common envelope evolution

Authors: Luke Chamandy, Adam Frank, Eric G. Blackman, Jonathan Carroll-Nellenback, Baowei Liu, Yisheng Tu, Jason Nordhaus, Zhuo Chen, Bo Peng

Abstract: Common envelope evolution (CEE) is presently a poorly understood, yet critical, process in binary stellar evolution. Characterizing the full 3D dynamics of CEE is difficult in part because simulating CEE is so computationally demanding. Numerical studies have yet to conclusively determine how the envelope ejects and a tight binary results, if only the binary potential energy is used to propel the… ▽ More Common envelope evolution (CEE) is presently a poorly understood, yet critical, process in binary stellar evolution. Characterizing the full 3D dynamics of CEE is difficult in part because simulating CEE is so computationally demanding. Numerical studies have yet to conclusively determine how the envelope ejects and a tight binary results, if only the binary potential energy is used to propel the envelope. Additional power sources might be necessary and accretion onto the inspiraling companion is one such source. Accretion is likely common in post-asymptotic giant branch (AGB) binary interactions but how it operates and how its consequences depend on binary separation remain open questions. Here we use high resolution global 3D hydrodynamic simulations of CEE with the adaptive mesh refinement (AMR) code AstroBEAR, to bracket the range of CEE companion accretion rates by comparing runs that remove mass and pressure via a subgrid accretion model with those that do not. The results show that if a pressure release valve is available, super-Eddington accretion may be common. Jets are a plausible release valve in these environments, and they could also help unbind and shape the envelopes. △ Less

Submitted 9 May, 2018; originally announced May 2018.

Comments: 14 pages, submitted to MNRAS

Generation of a Circumstellar Gas Disk by Hot Jupiter WASP-12b

Authors: Alex Debrecht, Jonathan Carroll-Nellenback, Adam Frank, Luca Fossati, Eric G. Blackman, Ian Dobbs-Dixon

Abstract: Observations of transiting extra-solar planets provide rich sources of data for probing the in-system environment. In the WASP-12 system, a broad depression in the usually-bright MgII h&k lines has been observed, in addition to atmospheric escape from the extremely hot Jupiter WASP-12b. It has been hypothesized that a translucent circumstellar cloud is formed by the outflow from the planet, causin… ▽ More Observations of transiting extra-solar planets provide rich sources of data for probing the in-system environment. In the WASP-12 system, a broad depression in the usually-bright MgII h&k lines has been observed, in addition to atmospheric escape from the extremely hot Jupiter WASP-12b. It has been hypothesized that a translucent circumstellar cloud is formed by the outflow from the planet, causing the observed signatures. We perform 3D hydrodynamic simulations of the full system environment of WASP-12, injecting a planetary wind and stellar wind from their respective surfaces. We find that a torus of density high enough to account for the lack of MgII h&k line core emission in WASP-12 can be formed in approximately 13 years. We also perform synthetic observations of the Lyman-alpha spectrum at different points in the planet's orbit, which demonstrate that significant absorption occurs at all points in the orbit, not just during transits, as suggested by the observations. △ Less

Submitted 1 May, 2018; originally announced May 2018.

Experimental Confirmation of the Standard Magnetorotational Instability Mechanism with a Spring-Mass Analogue

Authors: Derek M. H. Hung, Eric G. Blackman, Kyle J. Caspary, Erik P. Gilson, Hantao Ji

Abstract: The Magnetorotational Instability (MRI) has long been considered a plausibly ubiquitous mechanism to destabilize otherwise stable Keplerian flows to support radially outward transport of angular momentum. Such an efficient transport process would allow fast accretion in astrophysical objects such as stars and black holes to release copious kinetic energy that powers many of the most luminous sourc… ▽ More The Magnetorotational Instability (MRI) has long been considered a plausibly ubiquitous mechanism to destabilize otherwise stable Keplerian flows to support radially outward transport of angular momentum. Such an efficient transport process would allow fast accretion in astrophysical objects such as stars and black holes to release copious kinetic energy that powers many of the most luminous sources in the universe. But the standard MRI under a purely vertical magnetic field has heretofore never been directly measured despite numerous efforts over more than a decade. Here we report an unambiguous laboratory demonstration of the spring-mass analogue to the standard MRI by comparing motion of a spring-tethered ball within different rotating flows. The experiment corroborates the theory: efficient outward angular momentum transport manifests only for cases with a weak spring in quasi-Keperian flow. Our experimental method accomplishes this in a new way, thereby connecting solid and fluid mechanics to plasma astrophysics. △ Less

Submitted 11 December, 2018; v1 submitted 10 January, 2018; originally announced January 2018.

Comments: 7 pages; 4 figures; accepted by Communications Physics

Mass, energy, and momentum capture from stellar winds by magnetized and unmagnetized planets: implications for atmospheric erosion and habitability

Authors: Eric G. Blackman, John A. Tarduno

Abstract: The extent to which a magnetosphere protects its planetary atmosphere from stellar wind ablation depends upon how well it prevents energy and momentum exchange with the atmosphere and how well it traps otherwise escaping plasma. We focus on the former, and provide a formalism for estimating approximate upper limits on mass, energy, and momentum capture, and use them to constrain loss rates. Our ap… ▽ More The extent to which a magnetosphere protects its planetary atmosphere from stellar wind ablation depends upon how well it prevents energy and momentum exchange with the atmosphere and how well it traps otherwise escaping plasma. We focus on the former, and provide a formalism for estimating approximate upper limits on mass, energy, and momentum capture, and use them to constrain loss rates. Our approach quantifies a competition between the local deflection of incoming plasma by a planetary magnetic field and the increase in area for solar wind mass capture provided by this magnetosphere when the wind and planetary fields incur magnetic reconnection. Solar wind capture can be larger for a magnetized planet versus an unmagnetized planet with small ionosphere, even if the rate of energy transfer is less. We find this to be the case throughout Earth's history since the lunar forming impact-- the likely start of the geodynamo. Focusing of incoming charged particles near the magnetic poles can increase the local energy flux for a strong enough planetary field and weak enough wind, and we find that this may be the case for Earth's future. The future protective role of Earth's magnetic field would then depend on its ability to trap outgoing plasma. This competition between increased collection area and reduced inflow speed from a magnetosphere is likely essential in determining the net protective role of planetary fields and its importance in sustaining habitability, and can help explain the current ion loss rates from Mars, Venus and Earth. △ Less

Submitted 6 October, 2018; v1 submitted 2 January, 2018; originally announced January 2018.

Comments: 10 pages; 2 figures; accepted and in press, MNRAS

Derivation and precision of mean field electrodynamics with mesoscale fluctuations

Authors: Hongzhe Zhou, Eric G. Blackman, Luke Chamandy

Abstract: Mean field electrodynamics (MFE) facilitates practical modeling of secular, large scale properties of astrophysical or laboratory systems with fluctuations.Practitioners commonly assume wide scale separation between mean and fluctuating quantities, to justify equality of ensemble and spatial or temporal averages.Often however, real systems do not exhibit such scale separation. This raises two ques… ▽ More Mean field electrodynamics (MFE) facilitates practical modeling of secular, large scale properties of astrophysical or laboratory systems with fluctuations.Practitioners commonly assume wide scale separation between mean and fluctuating quantities, to justify equality of ensemble and spatial or temporal averages.Often however, real systems do not exhibit such scale separation. This raises two questions: (I) what are the appropriate generalized equations of MFE in the presence of mesoscale fluctuations? (II) how precise are theoretical predictions from MFE? We address both by first deriving the equations of MFE for different types of averaging, along with mesoscale correction terms that depend on the ratio of averaging scale to variation scale of the mean. We then show that even if these terms are small, predictions of MFE can still have a significant precision error. This error has an intrinsic contribution from the dynamo input parameters and a filtering contribution from differences in the way observations and theory are projected through the measurement kernel.Minimizing the sum of these contributions can produce an optimal scale of averaging that makes the theory maximally precise.The precision error is important to quantify when comparing to observations because it quantifies the resolution of predictive power. We exemplify these principles for galactic dynamos, comment on broader implications, and identify possibilities for further work. △ Less

Submitted 21 March, 2019; v1 submitted 6 October, 2017; originally announced October 2017.

Comments: 35 pages, 10 figures. Published in Journal of Plasma Physics; Mar. 21 2019: Typos of a sign in (3.6) and indices of $ξ$ in the paragraphs containing (B7)-(B9) corrected

Journal ref: Journal of Plasma Physics, 84(3) (2018)

Wind-accelerated orbital evolution in binary systems with giant stars

Authors: Zhuo Chen, Eric G. Blackman, Jason Nordhaus, Adam Frank, Jonathan Carroll-Nellenback

Abstract: Using 3D radiation-hydrodynamic simulations and analytic theory, we study the orbital evolution of asymptotic-giant-branch (AGB) binary systems for various initial orbital separations and mass ratios, and thus different initial accretion modes. The time evolution of binary separations and orbital periods are calculated directly from the averaged mass loss rate, accretion rate and angular momentum… ▽ More Using 3D radiation-hydrodynamic simulations and analytic theory, we study the orbital evolution of asymptotic-giant-branch (AGB) binary systems for various initial orbital separations and mass ratios, and thus different initial accretion modes. The time evolution of binary separations and orbital periods are calculated directly from the averaged mass loss rate, accretion rate and angular momentum loss rate. We separately consider spin-orbit synchronized and zero spin AGB cases. We find that the the angular momentum carried away by the mass loss together with the mass transfer can effectively shrink the orbit when accretion occurs via wind-Roche-lobe overflow. In contrast, the larger fraction of mass lost in Bondi-Hoyle-Lyttleton accreting systems acts to enlarge the orbit. Synchronized binaries tend to experience stronger orbital period decay in close binaries. We also find that orbital period decay is faster when we account for the nonlinear evolution of the accretion mode as the binary starts to tighten. This can increase the fraction of binaries that result in common envelope, luminous red novae, Type Ia supernovae and planetary nebulae with tight central binaries. The results also imply that planets in the the habitable zone around white dwarfs are unlikely to be found. △ Less

Submitted 31 August, 2017; v1 submitted 4 May, 2017; originally announced May 2017.

Comments: 10 pages, 10 figures and 6 tables submitted to MNRAS

Mass Transfer and Disc Formation in AGB Binary Systems

Authors: Zhuo Chen, Adam Frank, Eric G. Blackman, Jason Nordhaus, Jonathan Carroll-Nellenback

Abstract: We investigate mass transfer and the formation of disc in binary systems using a combination of numerical simulations and theory. We consider six models distinguished by binary separation, secondary mass and outflow mechanisms. Each system consists of an asymptotic-giant-branch (AGB) star and an accreting secondary. The AGB star loses its mass via a wind. In one of our six models, the AGB star inc… ▽ More We investigate mass transfer and the formation of disc in binary systems using a combination of numerical simulations and theory. We consider six models distinguished by binary separation, secondary mass and outflow mechanisms. Each system consists of an asymptotic-giant-branch (AGB) star and an accreting secondary. The AGB star loses its mass via a wind. In one of our six models, the AGB star incurs a short period of outburst. In all cases, the secondary accretes part of the ejected mass and also influences the mass-loss rate of the AGB star. The ejected mass may remain gravitationally bound to the binary system and form a circumbinary disk, or contribute to an accretion disk around the secondary. In other cases, the ejecta will escape the binary system. The accretion rate on to the secondary changes non-linearly with binary separation. In our closest binary simulations, our models exemplify the wind Roche lobe overflow while in our wide binary cases, the mass transfer exhibits Bondi-Hoyle accretion. The morphologies of the outflows in the binary systems are varied. The variety may provide clues to how the late AGB phase influences planetary nebulae shaping. We employ the adaptive-mesh-refinement code ASTROBEAR for our simulations and include ray-tracing, radiation transfer, cooling and dust formation. To attain the highest computational efficiency and the most stable results, all simulations are run in the corotating frame. △ Less

Submitted 13 April, 2017; v1 submitted 20 February, 2017; originally announced February 2017.

Comments: 14 pages, 31 figures and 4 tables. accepted by MNRAS