Source author record

Eric G. Blackman

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

ResearcherUnclaimed source record

astro-ph.SR astro-ph.HE physics.plasm-ph astro-ph.GA astro-ph astro-ph.EP physics.flu-dyn astro-ph.IM physics.space-ph astro-ph.CO nlin.CD physics.ao-ph physics.ed-ph physics.hist-ph physics.med-ph physics.soc-ph

Catalog footprint

What is connected

59works

16topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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

Building this map preview

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

preprint2024arXiv

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

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.

preprint2023arXiv

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

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.

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2021arXiv

Successive common envelope events from multiple planets

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.

preprint2021arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

Bipolar Planetary Nebulae from Outflow Collimation by Common Envelope Evolution

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

Generalized quenching of large-scale magnetic dynamos in anisotropic flows

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.

preprint2020arXiv

Improving TBI Protection Measures and Standards for Combat Helmets

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.

preprint2020arXiv

Necessity of ventilation for mitigating virus transmission quantified simply

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.

preprint2020arXiv

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

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.

preprint2019arXiv

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

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.

preprint2018arXiv

Accretion in Common Envelope Evolution

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.

preprint2016arXiv

Large scale dynamo action precedes turbulence in shearing box simulations of the magnetorotational instability

We study the dynamo generation (exponential growth) of large scale (planar averaged) fields in unstratified shearing box simulations of the magnetorotational instability (MRI). In contrast to previous studies restricted to horizontal ($x$-$y$) averaging, we demonstrate the presence of large scale fields when either horizontal or vertical ($y$-$z$) averaging is employed. By computing planar averaged fields and power spectra, we find large scale dynamo action in the early MRI growth phase---a previously unidentified feature. Fast growing horizontal low modes and fiducial vertical modes over a narrow range of wave numbers amplify these planar averaged fields in the MRI growth phase, before turbulence sets in. The large scale field growth requires linear fluctuations but not nonlinear turbulence (as defined by mode-mode coupling) and grows as a direct global mode of the MRI. Only by vertical averaging, can it be shown that the growth of horizontal low wavenumber MRI modes directly feed-back to the initial vertical field providing a clue as to why the large scale vertical field sustains against turbulent diffusion in the saturation regime. We compute the terms in the planar averaged mean field equations to identify the individual contributions to large scale field growth for both vertical and horizontal averaging. The large scale fields obtained from such vertical averaging are found to compare well with global cylindrical simulations and quasilinear analytical analysis from a previous study by Ebrahimi \& Blackman. We discuss the potential implications of these new results for understanding large scale MRI dynamo saturation and turbulence.

preprint2016arXiv

Vector cylindrical harmonics for low-dimensional convection models

Approximate empirical models of thermal convection can allow us to identify the essential properties of the flow in simplified form, and to produce empirical estimates using only a few parameters. Such "low-dimensional" empirical models can be constructed systematically by writing numerical or experimental measurements as superpositions of a set of appropriate basis modes, a process known as Galerkin projection. For three-dimensional convection in a cylinder, those basis modes should be vector-valued, mutually orthogonal, and defined in cylindrical coordinates. Here we construct such a basis set and demonstrate that it has these desired properties and boundary conditions when the exact constraint of incompressibility is relaxed. We show its use for representing sample simulation data and point out its potential for low-dimensional convection models.

preprint2015arXiv

Detecting the oldest geodynamo and attendant shielding from the solar wind: Implications for habitability

The onset and nature of the earliest geomagnetic field is important for understanding the evolution of the core, atmosphere and life on Earth. A record of the early geodynamo is preserved in ancient silicate crystals containing minute magnetic inclusions. These data indicate the presence of a geodynamo during the Paleoarchean, between 3.4 and 3.45 billion years ago. While the magnetic field sheltered Earth's atmosphere from erosion at this time, standoff of the solar wind was greatly reduced, and similar to that during modern extreme solar storms. These conditions suggest that intense radiation from the young Sun may have modified the atmosphere of the young Earth by promoting loss of volatiles, including water. Such effects would have been more pronounced if the field were absent or very weak prior to 3.45 billion years ago, as suggested by some models of lower mantle evolution. The frontier is thus trying to obtain geomagnetic field records that are >>3.45 billion-years-old, as well as constraining solar wind pressure for these times. In this review we suggest pathways for constraining these parameters and the attendant history of Earth's deep interior, hydrosphere and atmosphere. In particular, we discuss new estimates for solar wind pressure for the first 700 million years of Earth history, the competing effects of magnetic shielding versus solar ion collection, and bounds on the detection level of a geodynamo imposed by the presence of external fields. We also discuss the prospects for constraining Hadean-Paleoarchean magnetic field strength using paleointensity analyses of zircons.

preprint2015arXiv

On the resilience of helical magnetic fields to turbulent diffusion and the astrophysical implications

The extent to which large scale magnetic fields are susceptible to turbulent diffusion is important for interpreting the need for in situ large scale dynamos in astrophysics and for observationally inferring field strengths compared to kinetic energy. By solving coupled equations for magnetic energy and magnetic helicity in a system initiated with isotropic turbulence and an arbitrarily helical large scale field, we quantify the decay rate of the latter for a bounded or periodic system. The energy associated with the non-helical magnetic field rapidly decays by turbulent diffusion, but the decay rate of the helical component depends on whether the ratio of its magnetic energy to the turbulent kinetic energy exceeds a critical value given by M_{1,c} =(k_1/k_2)^2, where k_1 and k_2 are the wave numbers of the large and forcing scales. Turbulently diffusing helical fields to small scales while conserving magnetic helicity requires a rapid increase in total magnetic energy. As such, only when the helical fields are sub-critical can they so diffuse. When super-critical, the large scale helical field decays slowly, at a rate determined by microphysical dissipation even when macroscopic turbulence is present. Amplification of small scale magnetic helicity abates the turbulent diffusion. Two implications are that: (1) Standard arguments supporting the need for in situ large scale dynamos based on the otherwise rapid turbulent diffusion of large scale fields require re-thinking since only the non-helical field is so diffused in a closed system. Boundary terms could however provide potential pathways for rapid change of the large scale helical field. (2) Since M_{1,c} <<1 for k_1 << k_2, the presence of long-lived ordered large scale helical fields, as in extragalactic jets, does not guarantee that the magnetic field dominates the kinetic energy.

preprint2014arXiv

Explaining the Observed Relation Between Stellar Activity and Rotation

Observations of late-type main-sequence stars have revealed empirical scalings of coronal activity versus rotation period or Rossby number $Ro$ (a ratio of rotation period to convective turnover time) which has hitherto lacked explanation. For $Ro >> 1$, the activity observed as X-ray to bolometric flux varies as $Ro^{-q}$ with $2\le q \le 3$, whilst $|q| < 0.12$ for $Ro << 1$. Here we explain the transition between these two regimes and the power law in the $Ro >> 1$ regime by constructing an expression for the coronal luminosity based on dynamo magnetic field generation and magnetic buoyancy. We explain the $Ro<<1$ behavior from the inference that observed rotation is correlated with internal differential rotation and argue that once the shear time scale is shorter than the convective turnover time, eddies will be shredded on the shear time scale and so the eddy correlation time actually becomes the shear time and the convection time drops out of the equations. We explain the $Ro >> 1$ behavior using a dynamo saturation theory based on magnetic helicity buildup and buoyant loss.

preprint2014arXiv

Magnetic Helicity and Large Scale Magnetic Fields: A Primer

Magnetic fields of laboratory, planetary, stellar, and galactic plasmas commonly exhibit significant order on large temporal or spatial scales compared to the otherwise random motions within the hosting system. Such ordered fields can be measured in the case of planets, stars, and galaxies, or inferred indirectly by the action of their dynamical influence, such as jets. Whether large scale fields are amplified in situ or a remnant from previous stages of an object's history is often debated for objects without a definitive magnetic activity cycle. Magnetic helicity, a measure of twist and linkage of magnetic field lines, is a unifying tool for understanding large scale field evolution for both mechanisms of origin. Its importance stems from its two basic properties: (1) magnetic helicity is typically better conserved than magnetic energy; and (2) the magnetic energy associated with a fixed amount of magnetic helicity is minimized when the system relaxes this helical structure to the largest scale available. Here I discuss how magnetic helicity has come to help us understand the saturation of and sustenance of large scale dynamos, the need for either local or global helicity fluxes to avoid dynamo quenching, and the associated observational consequences. I also discuss how magnetic helicity acts as a hindrance to turbulent diffusion of large scale fields, and thus a helper for fossil remnant large scale field origin models in some contexts. I briefly discuss the connection between large scale fields and accretion disk theory as well. The goal here is to provide a conceptual primer to help the reader efficiently penetrate the literature.

preprint2014arXiv

On characterizing nonlocality and anisotropy for the magnetorotational instability

The extent to which angular momentum transport in accretion discs is primarily local or non-local and what determines this is an important avenue of study for understanding accretion engines. Taking a step along this path, we analyze simulations of the magnetorotational instability (MRI) by calculating energy and stress power spectra in stratified isothermal shearing box simulations in several new ways. We divide our boxes in two regions, disc and corona where the disc is the MRI unstable region and corona is the magnetically dominated region. We calculate the fractional power in different quantities, including magnetic energy and Maxwell stresses and find that they are dominated by contributions from the lowest wave numbers. This is even more dramatic for the corona than the disc, suggesting that transport in the corona region is dominated by larger structures than the disc. By calculating averaged power spectra in one direction of $k$ space at a time, we also show that the MRI turbulence is strongly anisotropic on large scales when analyzed by this method, but isotropic on small scales. Although the shearing box itself is meant to represent a local section of an accretion disc, the fact that the stress and energy are dominated by the largest scales highlights that the locality is not captured within the box. This helps to quantify the intuitive importance of global simulations for addressing the question of locality of transport, for which similar analyses can be performed.

preprint2014arXiv

Ribbons characterize magnetohydrodynamic magnetic fields better than lines: a lesson from dynamo theory

Blackman & Brandenburg argued that magnetic helicity conservation in dynamo theory can in principle be captured by diagrams of mean field dynamos when the magnetic fields are represented by ribbons or tubes, but not by lines. Here we present such a schematic ribbon diagram for the $α^2$ dynamo that tracks magnetic helicity and provides distinct scales of large scale magnetic helicity, small scale magnetic helicity, and kinetic helicity involved in the process. This also motivates our construction of a new "2.5 scale" minimalist generalization of the helicity-evolving equations for the α^2 dynamo that separately allows for these three distinct length scales while keeping only two dynamical equations. We solve these equations and, as in previous studies, find that the large scale field first grows at a rate independent of the magnetic Reynolds number R_M before quenching to an R_M dependent regime. But we also show that the larger the ratio of the wavenumber where the small scale current helicity resides to that of the forcing scale, the earlier the non-linear dynamo quenching occurs, and the weaker the large scale field is at the turnoff from linear growth. The harmony between the theory and the schematic diagram exemplifies a general lesson that magnetic fields in MHD are better visualized as two-dimensional ribbons (or pairs of lines) rather than single lines.

preprint2014arXiv

Sensitivity of the Magnetorotational Instability to the shear parameter in stratified simulations

The magnetorotational instability (MRI) is a shear instability and thus its sensitivity to the shear parameter $q = - d\lnΩ/d\ln r $ is of interest to investigate. Motivated by astrophysical disks, most (but not all) previous MRI studies have focused on the Keplerian value of $ q=1.5$. Using simulation with 8 vertical density scale heights, we contribute to the subset of studies addressing the the effect of varying $q$ in stratified numerical simulations. We discuss why shearing boxes cannot easily be used to study $q>2$ and thus focus on $q<2$. As per previous simulations, which were either unstratified or stratified with a smaller vertical domain, we find that the $q$ dependence of stress for the stratified case is not linear, contrary to the Shakura-Sunyaev model. We find that the scaling agrees with \cite{1996MNRAS.281L..21A} who found it to be proportional to the shear to vorticity ratio $q/(2-q)$. We also find however, that the shape of the magnetic and kinetic energy spectra are relatively insensitive to $q$ and that the ratio of Maxwell stress to magnetic energy ratio also remains nearly independent of $q$. This is consistent with a theoretical argument in which the rate of amplification of the azimuthal field depends linearly on $q$ and the turbulent correlation time $τ$ depends inversely on $q$. As such, we measure the correlation time of the turbulence and find that indeed it is inversely proportional to $q$.

preprint2014arXiv

Using Kinematic Properties of Pre-Planetary Nebulae to Constrain Engine Paradigms

Some combination of binary interactions and accretion plausibly conspire to produce the ubiquitous collimated outflows from planetary nebulae (PN) and their presumed pre-planetary nebulae (PPN) precursors. But which accretion engines are viable? The difficulty in observationally resolving the engines warrants the pursuit of indirect constraints. We show how kinematic outflow data for 19 PPN can be used to determine the minimum required accretion rates. We consider main sequence (MS) and white dwarf (WD) accretors and five example accretion rates inferred from published models to compare with the minima derived from outflow momentum conservation. While our primary goal is to show the method in anticipation of more data and better theoretical constraints, taking the present results at face value already rule out modes of accretion: Bondi-Hoyle Lyttleton (BHL) wind accretion and wind Roche lobe overflow (M-WRLOF, based on Mira parameters) are too feeble for all 19/19 objects for a MS accretor. For a WD accretor, BHL is ruled out for 18/19 objects and M-WRLOF for 15/19 objects. Roche lobe overflow (RLOF) from the primary at the Red Rectangle level can accommodate 7/19 objects, though RLOF modes with higher accretion rates are not yet ruled out. Accretion modes operating from within common envelope evolution can accommodate all 19 objects, if jet collimation can be maintained. Overall, sub-Eddington rates for a MS accretor are acceptable but 8/19 would require super-Eddington rates for a WD.

preprint2013arXiv

A study of large scale dynamo growth rates from numerical simulations and implications for mean field theories

Understanding large scale magnetic field growth in turbulent plasmas in the magnetohydrodynamic limit is a goal of magnetic dynamo theory. In particular, assessing how well large scale helical field growth and saturation in simulations matches that predicted by existing theories is important for progress. Using numerical simulations of isotropically forced turbulence without large scale shear with the implications, we focus on several aspects of this comparison that have not been previously tested: (1) Leading mean field dynamo theories which break the field into large and small scales predict that large scale helical field growth rates are determined by the difference between kinetic helicity and current helicity with no dependence on the non-helical energy in small scale magnetic fields. Our simulations show that the growth rate of the large scale field from fully helical forcing is indeed unaffected by the presence or absence of small scale magnetic fields amplified in a precursor non-helical dynamo. However, because the precursor non helical dynamo in our simulations produced fields that were strongly sub-equipartition with respect to the kinetic energy, we cannot yet rule out the potential influence of stronger non- helical small scale fields. (2) We have identified two features in our simulations which cannot be explained by the most minimalist versions of two-scale mean field theory: (i) fully helical small scale forcing produces significant non-helical large scale magnetic energy and (ii) the saturation of the large scale field growth is time-delayed with respect to what minimalist theory predicts. We comment on desirable generalizations to the theory in this context and future desired work.

preprint2013arXiv

Mass loading and knot formation in AGN jets by stellar winds

Jets from active galaxies propagate from the central black hole out to the radio lobes on scales of hundreds of kiloparsecs. The jets may encounter giant stars with strong stellar winds and produce observable signatures. For strong winds and weak jets, the interaction may truncate the jet flow during its transit via the mass loading. For weaker jets, the interaction can produce knots in the jet. We present recent 3DMHD numerical simulations to model the evolution of this jet-wind interaction and its observational consequences. We explore (i) the relative mechanical luminosity of the radio jets and the stellar winds (ii) the impact parameter between the jets' axis and the stellar orbital path (iii) the relative magnetic field strength of the jets and the stellar winds.

preprint2013arXiv

MHD Shock-Clump Evolution with Self-Contained Magnetic Fields

We study the interaction of strong shock waves with magnetized clumps. Previous numerical work focused on the simplified scenario in which shocked clumps are immersed in a globally uniform magnetic field that extends through both the clump and the ambient medium. Here we consider the complementary circumstance in which the field is completely self-contained within the clumps. This could arise naturally during clump formation via dynamical or thermal instabilities for example as magnetic field pinches off from the ambient medium. Using our AMR MHD code AstroBEAR, we carry out a series of simulations with magnetized clumps that have different self-contained magnetic field configurations. We find that the clump and magnetic evolution are sensitive to the fraction of magnetic field aligned with versus perpendicular to the shock normal. The relative strength of magnetic pressure and tension in the different field configurations allows us to analytically understand the different cases of post-shock evolution. We also show how turbulence and the mixing it implies depends of the initial field configuration and suggest ways in which observed shock-clump morphology may be used as a proxy for identifying internal field topologies a posteriori.

preprint2013arXiv

On deriving flux freezing in magnetohydrodynamics by direct differentiation

The magnetic flux freezing theorem or Alfvén's theorem, is a basic principle of ideal magnetohydrodynamics (MHD), a commonly used approximation to describe the aspects of astrophysical and laboratory plasmas. The theorem states that the magnetic flux---the integral of magnetic field penetrating a surface---is conserved in time as that surface is distorted in in time by fluid motions. Pedagogues of MHD commonly derive flux freezing without showing how to take the material derivative of a general flux integral and/or assuming a vanishing field divergence from the outset. Here I avoid these shortcomings and derive flux freezing by direct differentiation, explicitly using a Jacobian to transform between the evolving field-penetrating surface at different times. The approach is instructive for its generality and helps elucidate the role of magnetic monopoles in breaking flux freezing. The paucity of appearances of this derivation in standard MHD texts suggests that its pedagogic value is underappreciated.

preprint2013arXiv

Particle-in-cell simulations of particle energization via shock drift acceleration from low Mach number quasi-perpendicular shocks in solar flares

Low Mach number, high beta fast mode shocks can occur in the magnetic reconnection outflows of solar flares. These shocks, which occur above flare loop tops, may provide the electron energization responsible for some of the observed hard X-rays and contemporaneous radio emission. Here we present new 2D particle-in-cell simulations of low Mach number/high beta quasi-perpendicular shocks. The simulations show that electrons above a certain energy threshold experience shock-drift-acceleration. The transition energy between the thermal and non-thermal spectrum and the spectral index from the simulations are consistent with some of the X-ray spectra from RHESSI in the energy regime of $E\lesssim 40\sim 100$ keV. Plasma instabilities associated with the shock structure such as the modified-two-stream and the electron whistler instabilities are identified using numerical solutions of the kinetic dispersion relations. We also show that the results from PIC simulations with reduced ion/electron mass ratio can be scaled to those with the realistic mass ratio.

preprint2013arXiv

Resilience of helical fields to turbulent diffusion II: direct numerical simulations

The recent study of Blackman and Subramanian (Paper I) indicates that large scale helical magnetic fields are resilient to turbulent diffusion in the sense that helical fields stronger than a critical value, decay on slow (~resistively mediated), rather than fast ($\sim$ turbulent) time scales. This gives more credence to potential fossil field origin models of the magnetic fields in stars, galaxies and compact objects. We analyze a suite of direct numerical simulations (DNS) of decaying large scale helical magnetic fields in the presence of non-helical turbulence to further study the physics of helical field decay. We study two separate cases: (1) the initial field is large enough to decay resistively, is tracked until it transitions to decay fast, and the critical large scale helical field at that transition is sought; (2) the case of Paper I, wherein there is a critical initial helical field strength below which the field undergoes fast decay right from the beginning. In case (1), both DNS and solutions of the two scale model (from Paper 1), reveal that the transition energy, $E_{c1}$, is independent of the scale of the turbulent forcing, within a small range of $\Rm$. We also find that the kinetic alpha, $α_K$, is subdominant to magnetic alpha, $α_M$, in the DNS, justifying an assumption in the two scale model. For case (2), we show exact solutions of two scale model in the limit of $η\rightarrow 0$ in fully helical case, leading to the transition energy, $E_{c2} = (k_1/\kf)^2 M_{eq}$, where $k_1$ and $\kf$ are the large scale and small turbulent forcing scale respectively and $M_{eq}$ is the equipartition magnetic energy. The DNS in this case agree qualitatively with the two scale model but the $R_M$ currently achievable, is too small to satisfy a condition $3/R_M << (k_1/k_f)^2$, necessary to robustly reveal the transition, $E_{c2}$ (Abridged).

preprint2013arXiv

Size of discs formed by wind accretion in binaries can be underestimated if the role of wind-driving force is ignored

Binary systems consisting of a secondary accreting form a wind-emitting primary are ubiquitous in astrophysics. The phenomenology of such Bondi-Hoyle-Lyttleton (BHL) accretors is particularly rich when an accretion disc forms around the secondary. The outer radius of such discs is commonly estimated from the net angular momentum produced by a density variation of material across the BHL or Bondi accretion cylinder, as the latter is tilted with respect to the direction to the primary due to orbital motion. But this approach has ignored the fact that the wind experiences an outward driving force that the secondary does not. In actuality, the accretion stream falls toward a retarded point in the secondary's orbit as the secondary is pulled toward the primary relative to the stream. The result is a finite separation or "accretion stream impact parameter" (ASIP) separating the secondary and stream. When the orbital radius a_o exceeds the BHL radius r_b, the ratio of outer disc radius estimated as the ASIP to the conventional estimate a_o^{1/2}/r_b^{1/2}>1. We therefore predict that discs will form at larger radii from the secondary than traditional estimates. This agrees with the importance of the ASIP emphasized by Huarte-Espinosa et al. (2013) and the practical consequence that resolving the initial outer radius of such an accretion disc in numerical simulations can be less demanding than what earlier estimates would suggest.

preprint2013arXiv

The Formation and Evolution of Wind-Capture Disks In Binary Systems

We study the formation, evolution and physical properties of accretion disks formed via wind capture in binary systems. Using the AMR code AstroBEAR, we have carried out high resolution 3D simulations that follow a stellar mass secondary in the co-rotating frame as it orbits a wind producing AGB primary. We first derive a resolution criteria, based on considerations of Bondi-Hoyle flows, that must be met in order to properly resolve the formation of accretion disks around the secondary. We then compare simulations of binaries with three different orbital radii (10, 15, 20 AU). Disks are formed in all three cases, however the size of the disk and, most importantly, its accretion rate decreases with orbital radii. In addition, the shape of the orbital motions of material within the disk becomes increasingly elliptical with increasing binary separation. The flow is mildly unsteady with "fluttering" around the bow shock observed. The disks are generally well aligned with the orbital plane after a few binary orbits. We do not observe the presence of any large scale, violent instabilities (such as the flip-flop mode). For the first time, moreover, it is observed that the wind component that is accreted towards the secondary has a vortex tube-like structure, rather than a column-like one as it was previously thought. In the context of AGB binary systems that might be precursors to Pre-Planetary and Planetary Nebula, we find that the wind accretion rates at the chosen orbital separations are generally too small to produce the most powerful outflows observed in these systems if the companions are main sequence stars but marginally capable if the companions are white dwarfs. It is likely that many of the more powerful PPN and PN involve closer binaries than the ones considered here. The results also demonstrate principles of broad relevance to all wind-capture binary systems.

preprint2012arXiv

Accretion Disks and Dynamos: Toward a Unified Mean Field Theory

Conversion of gravitational energy into radiation in accretion discs and the origin of large scale magnetic fields in astrophysical rotators have often been distinct topics of research. In semi-analytic work on both problems it has been useful to presume large scale symmetries, necessarily resulting in mean field theories. MHD turbulence makes the underlying systems locally asymmetric and nonlinear. Synergy between theory and simulations should aim for the development of practical mean field models that capture essential physics and can be used for observational modeling. Mean field dynamo (MFD) theory and alpha-viscosity accretion theory exemplify such ongoing pursuits. 21st century MFD theory has more nonlinear predictive power compared to 20th century MFD theory, whereas accretion theory is still in a 20th century state. In fact, insights from MFD theory are applicable to accretion theory and the two are artificially separated pieces of what should be a single theory. I discuss pieces of progress that provide clues toward a unified theory. A key concept is that large scale magnetic fields can be sustained via local or global magnetic helicity fluxes or via relaxation of small scale magnetic fluctuations, without the kinetic helicity driver of 20th century textbooks. These concepts may help explain the formation of large scale fields that supply non-local angular momentum transport via coronae and jets in a unified theory of accretion and dynamos. In diagnosing the role of helicities and helicity fluxes in disk simulations, each disk hemisphere should be studied separately to avoid being misled by cancelation that occurs as a result of reflection asymmetry. The fraction of helical field energy in disks is expected to be small compared to the total field in each hemisphere as a result of shear, but can still be essential for large scale dynamo action.

preprint2012arXiv

Comparing Poynting flux dominated magnetic tower jets with kinetic-energy dominated jets

Magnetic Towers represent one of two fundamental forms of MHD outflows. Driven by magnetic pressure gradients, these flows have been less well studied than magneto-centrifugally launched jets even though magnetic towers may well be as common. Here we present new results exploring the behavior and evolution of magnetic tower outflows and demonstrate their connection with pulsed power experimental studies and purely hydrodynamic jets which might represent the asymptotic propagation regimes of magneto-centrifugally launched jets. High-resolution AMR MHD simulations (using the AstroBEAR code) provide insights into the underlying physics of magnetic towers and help us constrain models of their propagation. Our simulations have been designed to explore the effects of thermal energy losses and rotation on both tower flows and their hydro counterparts. We find these parameters have significant effects on the stability of magnetic towers, but mild effects on the stability of hydro jets. Current-driven perturbations in the Poynting Flux Dominated (PDF) towers are shown to be amplified in both the cooling and rotating cases. Our studies of the long term evolution of the towers show that the formation of weakly magnetized central jets within the tower are broken up by these instabilities becoming a series of collimated clumps which magnetization properties vary over time. In addition to discussing these results in light of laboratory experiments, we address their relevance to astrophysical observations of young star jets and outflow from highly evolved solar type stars.

preprint2012arXiv

From Bipolar to Elliptical: Simulating the Morphological Evolution of Planetary Nebulae

The majority of Proto-planetary nebulae (PPN) are observed to have bipolar morphologies. The majority of mature PN are observed to have elliptical shapes. In this paper we address the evolution of PPN/PN morphologies attempting to understand if a transition from strongly bipolar to elliptical shape can be driven by changes in the parameters of the mass loss process. To this end we present 2.5D hydrodynamical simulations of mass loss at the end stages of stellar evolution for intermediate mass stars. We track changes in wind velocity, mass loss rate and mass loss geometry. In particular we focus on the transition from mass loss dominated by a short duration jet flow (driven during the PPN phase) to mass loss driven by a spherical fast wind (produced by the central star of the PN). We address how such changes in outflow characteristics can change the nebula from a bipolar to an elliptical morphology. Our results show that including a period of jet formation in the temporal sequence of PPN to PN produces realistic nebular synthetic emission geometries. More importantly such a sequence provides insight, in principle, into the apparent difference in morphology statistics characterizing PPN and PN systems. In particular we find that while jet driven PPN can be expected to be dominated by bipolar morphologies, systems that begin with a jet but are followed by a spherical fast wind will evolve into elliptical nebulae. Furthermore, we find that spherical nebulae are highly unlikely to ever derive from either bipolar PPN or elliptical PN.

preprint2012arXiv

Not Much Helicity is Needed to Drive Large Scale Dynamos

Understanding the in situ amplification of large scale magnetic fields in turbulent astrophysical rotators has been a core subject of dynamo theory. When turbulent velocities are helical, large scale dynamos that substantially amplify fields on scales that exceed the turbulent forcing scale arise, but the minimum sufficient fractional kinetic helicity f_h,C has not been previously well quantified. Using direct numerical simulations for a simple helical dynamo, we show that f_h,C decreases as the ratio of forcing to large scale wave numbers k_F/k_min increases. From the condition that a large scale helical dynamo must overcome the backreaction from any non-helical field on the large scales, we develop a theory that can explain the simulations. For k_F/k_min>8 we find f_h,C< 3%, implying that very small helicity fractions strongly influence magnetic spectra for even moderate scale separation.

preprint2012arXiv

On the structure and stability of magnetic tower jets

Modern theoretical models of astrophysical jets combine accretion, rotation, and magnetic fields to launch and collimate supersonic flows from a central source. Near the source, magnetic field strengths must be large enough to collimate the jet requiring that the Poynting flux exceeds the kinetic-energy flux. The extent to which the Poynting flux dominates kinetic energy flux at large distances from the engine distinguishes two classes of models. In magneto-centrifugal launch (MCL) models, magnetic fields dominate only at scales $\lesssim 100$ engine radii, after which the jets become hydrodynamically dominated (HD). By contrast, in Poynting flux dominated (PFD) magnetic tower models, the field dominates even out to much larger scales. To compare the large distance propagation differences of these two paradigms, we perform 3-D ideal MHD AMR simulations of both HD and PFD stellar jets formed via the same energy flux. We also compare how thermal energy losses and rotation of the jet base affects the stability in these jets. For the conditions described, we show that PFD and HD exhibit observationally distinguishable features: PFD jets are lighter, slower, and less stable than HD jets. Unlike HD jets, PFD jets develop current-driven instabilities that are exacerbated as cooling and rotation increase, resulting in jets that are clumpier than those in the HD limit. Our PFD jet simulations also resemble the magnetic towers that have been recently created in laboratory astrophysical jet experiments.

preprint2012arXiv

Particle-in-cell simulations of particle energization from low Mach number fast mode shocks

Astrophysical shocks are often studied in the high Mach number limit but weakly compressive fast shocks can occur in magnetic reconnection outflows and are considered to be a site of particle energization in solar flares. Here we study the microphysics of such perpendicular, low Mach number collisionless shocks using two-dimensional particle-in-cell (PIC) simulations with a reduced ion/electron mass ratio and employ a moving wall boundary method for initial generation of the shock. This moving wall method allows for more control of the shock speed, smaller simulation box sizes, and longer simulation times than the commonly used fixed wall, reflection method of shock formation. Our results, which are independent of the shock formation method, reveal the prevalence shock drift acceleration (SDA) of both electron and ions in a purely perpendicular shock with Alfvén Mach number $M_A=6.8$ and ratio of thermal to magnetic pressure $β=8$. We determine the respective minimum energies required for electrons and ions to incur SDA. We derive a theoretical electron distribution via SDA that compares to the simulation results. We also show that a modified two-stream instability due to the incoming and reflecting ions in the shock transition region acts as the mechanism to generate collisionless plasma turbulence that sustains the shock.

preprint2011arXiv

Anderson and Stoner Published White Dwarf Mass Limits Before Chandrasekhar

In their engaging recountals of Chandrasekhar's extraordinary career (Physics Today, vol 63, Issue 12, Dec 2010), neither Dyson nor Wali mention that Chandrasekhar was the third person not the first, to publish a white dwarf mass limit incorporating a relativistic treatment of degenerate electrons. As it has become a common misconception that Chandrasekhar was the first, a clarifying reminder on this historical point is warranted. In short, the white dwarf mass limit widely attributed to Chandrasekhar (1931) should be the specific white dwarf mass limit calculated for a polytrope. The insight that a relativistic treatment of degeneracy leads to the existence of a white dwarf mass limit first appeared in papers of W. Anderson (1929) and E.C. Stoner (1930) for a uniform density star. Accordingly, Chandrasekhar (1931) cites Stoner (1930) and points out that the polytrope white dwarf mass limit is less than Stoner's uniform density white dwarf mass limit by about 20%.

preprint2011arXiv

Fast Shocks From Magnetic Reconnection Outflows

Magnetic reconnection is commonly perceived to drive flow and particle acceleration in flares of solar, stellar, and astrophysical disk coronae but the relative roles of different acceleration mecha- nisms in a given reconnection environment are not well understood. We show via direct numerical simulations that reconnection outflows produce weak fast shocks, when conditions for fast recon- nection are met and the outflows encounter an obstacle. The associated compression ratios lead to a Fermi acceleration particle spectrum that is significantly steeper than the strong fast shocks commonly studied, but consistent with the demands of solar flares. While this is not the only acceleration mechanism operating in a reconnection environment, it is plausibly a ubiquitous one.

preprint2011arXiv

Modeling Poynting flux vs. kinetic-energy dominated jets

We present 3D-MHD AMR simulations of Poynting flux dominated (PFD) jets formed by injection of magnetic energy. We compare their evolution with a hydrodynamic jet which is formed by injecting kinetic energy with the same energy flux than the PFD jets. We predict characteristic emission distributions for each of these jets. Current-driven perturbations in PFD jets are amplified by both cooling and rotation for the regimes studied: Shocks and thermal pressure support are weakened by cooling, making the jets more susceptible to kinking. Rotation amplifies the toroidal magnetic field which also exacerbates the kink instability.

preprint2010arXiv

From Bipolar to Elliptical: Morphological Changes in the Temporal Evolution of PN

Proto-planetary nebulae (pPN) and planetary nebulae (PN) seem to be formed by interacting winds from asymptotic giant branch (AGB) stars. The observational issue that most pPN are bipolar but most older PN are elliptical is addressed. We present 2.5D hydrodynamical numerical simulations of episodic cooling interacting winds to investigate the long term evolution of PN morphologies. We track wind acceleration, decrease in mass-loss and episodic change in wind geometry from spherical (AGB) to collimated (pPN) and back to spherical again (PN). This outflow sequence is found to produce realistic PN dynamics and morphological histories. Effects from different AGB distributions and jet duty cycles are also investigated.

preprint2010arXiv

Quantifying the Imprecision of Accretion Theory and Implications for Multi-Epoch Observations of Protoplanetary Discs

If accretion disc emission results from turbulent dissipation, then axisymmetric accretion theory must be used as a mean field theory: turbulent flows are at most axisymmetric only when suitably averaged. Spectral predictions therefore have an intrinsic imprecision that must be quantified to interpret the variability exhibited by a source observed at different epochs. We quantify contributions to the stochastic imprecision that come from azimuthal and radial averaging and show that the imprecision is minimized for a particular choice of radial averaging, which in turn, corresponds to an optimal spectral resolution of a telescope for a spatially unresolved source. If the optimal spectral resolution is less than that of the telescope then the data can be binned to compare to the theoretical prediction of minimum imprecision. Little stochastic variability is predicted at radii much larger than that at which the dominant eddy turnover time ($\sim$ orbit time) exceeds the time interval between observations; the epochs would then be sampling the same member of the stochastic ensemble. We discuss the application of these principles to protoplanetary discs for which there is presently a paucity of multi-epoch data but for which such data acquisition projects are underway.

preprint2009arXiv

Comparisons and Connections between Mean Field Dynamo Theory and Accretion Disc Theory

The origin of large scale magnetic fields in astrophysical rotators, and the conversion of gravitational energy into radiation near stars and compact objects via accretion have been subjects of active research for a half century. Magnetohydrodynamic turbulence makes both problems highly nonlinear, so both subjects have benefitted from numerical simulations.However, understanding the key principles and practical modeling of observations warrants testable semi-analytic mean field theories that distill the essential physics. Mean field dynamo (MFD) theory and alpha-viscosity accretion disc theory exemplify this pursuit. That the latter is a mean field theory is not always made explicit but the combination of turbulence and global symmetry imply such. The more commonly explicit presentation of assumptions in 20th century textbook MFDT has exposed it to arguably more widespread criticism than incurred by 20th century alpha-accretion theory despite complementary weaknesses. In the 21st century however, MFDT has experienced a breakthrough with a dynamical saturation theory that consistently agrees with simulations. Such has not yet occurred in accretion disc theory, though progress is emerging. Ironically however, for accretion engines, MFDT and accretion theory are presently two artificially uncoupled pieces of what should be a single coupled theory. Large scale fields and accretion flows are dynamically intertwined because large scale fields likely play a key role in angular momentum transport. I discuss and synthesize aspects of recent progress in MFDT and accretion disc theory to suggest why the two likely conspire in a unified theory.

preprint2009arXiv

Coronae as Consequence of Large Scale Magnetic Fields in Turbulent Accretion Disks

Non-thermal X-ray emission in compact accretion engines can be interpreted to result from magnetic dissipation in an optically thin magnetized corona above an optically thick accretion disk. If coronal magnetic field originates in the disk and the disk is turbulent, then only magnetic structures large enough for their turbulent shredding time to exceed their buoyant rise time survive the journey to the corona. We use this concept and a physical model to constrain the minimum fraction of magnetic energy above the critical scale for buoyancy as a function of the observed coronal to bolometric emission. Our results suggest that a significant fraction of the magnetic energy in accretion disks resides in large scale fields, which in turn provides circumstantial evidence for significant non-local transport phenomena and the need for large scale magnetic field generation. For the example of Seyfert AGN, for which of order 30 per cent of the bolometric flux is in the X-ray band, we find that more than 20 per cent of the magnetic energy must be of large enough scale to rise and dissipate in the corona.

preprint2009arXiv

Magnetic Nested-wind Scenarios for Bipolar Outflows: Pre-planetary and YSO nebular shaping

We present results of a series of magnetohydrodynamic (MHD) and hydro- dynamic (HD) 2.5D simulations of the morphology of outflows driven by nested wide-angle winds - i.e. winds which eminate from a central star as well as from an orbiting accretion disk. While our results are broadly relevent to nested wind systems we have tuned the parameters of the simulations to touch on issues in both Young Stellar Objects and Planetary Nebula studies. In particular our studies connect to open issues in the early evolution of Planetary Nebulae. We find that nested MHD winds exhibit marked morphological differences from the single MHD wind case along both dimensions of the flow. Nested HD winds on the other hand give rise mainly to geometric distortions of an outflow that is topologically similar to the flow arising from a single stellar HD wind. Our MHD results are insensitive to changes in ambient temperature between ionized and un-ionized circumstellar environments. The results are sensitive to the relative mass-loss rates, and to the relative speeds of the stellar and disk winds. We also present synthetic emission maps of both nested MHD and HD simulations. We find that nested MHD winds show knots of emission appearing on-axis that do not appear in the HD case.

preprint2009arXiv

Stochastic Wobble of Accretion Discs and Jets from Turbulent Rocket Torques

Models of accretion discs and their associated outflows often incorporate assumptions of axisymmetry and symmetry across the disc plane. However, for turbulent discs these symmetries only apply to averaged quantities and do not apply locally. The local asymmetries can induce local imbalances in outflow power across the disc mid-plane, which can in turn induce local tilting torques. Here we calculate the effect of the resulting stochastic torques on disc annuli that are a consequence of standard mean field accretion disc models. The torques induce a random walk of the vector perpendicular to the plane of each averaged annulus. This random walk is characterized by a radially dependent diffusion coefficient which we calculate for small angle tilt. We use the coefficient to calculate a radially dependent time scale for annular tilt and associated jet wobble. The wobble time depends on the square of the wander angle so the age of a given system determines the maximum wobble angle. We apply this to examples of blazars, young stellar objects and binary engines of pre-planetary nebulae and microquasars. It is noteworthy that for an averaging time $t_w\sim 3$ days, we estimate a wobble angle for jets in SS433 of $θ\sim 0.8$ degrees, not inconsistent with observational data. In general the non-periodic nature of the stochastic wobble could distinguish it from faster periodic jet precession.

preprint2008arXiv

A Search for Molecular Gas in the Nucleus of M87 and Implications for the Fueling of Supermassive Black Holes

Supermassive black holes in giant elliptical galaxies are remarkably faint given their expected accretion rates. This motivates models of radiatively inefficient accretion, due to either ion-electron thermal decoupling, generation of outflows that inhibit accretion, or settling of gas to a gravitationally unstable disk that forms stars in preference to feeding the black hole. The latter model predicts the presence of cold molecular gas in a thin disk around the black hole. Here we report Submillimeter Array observations of the nucleus of the giant elliptical galaxy M87 that probe 230 GHz continuum and CO(J=2--1) line emission. Continuum emission is detected from the nucleus and several knots in the jet, including one that has been undergoing flaring behavior. We estimate a conservative upper limit on the mass of molecular gas within ~100pc and +-400km/s line of sight velocity of the central black hole of ~8x10^6Msun, which includes an allowance for possible systematic errors associated with subtraction of the continuum. Ignoring such errors, we have a 3 sigma sensitivity to about 3x10^6Msun. In fact, the continuum-subtracted spectrum shows weak emission features extending up to 4 sigma above the RMS dispersion of the line-free channels. These may be artifacts of the continuum subtraction process. Alternatively, if they are interpreted as CO emission, then the implied molecular gas mass is ~5x10^6Msun spread out over a velocity range of 700km/s. These constraints on molecular gas mass are close to the predictions of the model of self-gravitating, star-forming accretion disks fed by Bondi accretion (Tan & Blackman 2005).

preprint2008arXiv

Identifying Deficiencies of Standard Accretion Disk Theory: Lessons from a Mean-Field Approach

Turbulent viscosity is frequently used in accretion disk theory to replace the microphysical viscosity in order to accomodate the observational need for in- stabilities in disks that lead to enhanced transport. However, simply replacing the microphysical transport coefficient by a single turbulent transport coeffi- cient hides the fact that the procedure should formally arise as part of a closure in which the hydrodynamic or magnetohydrodynamic equations are averaged, and correlations of turbulent fluctuations are replaced by transport coefficients. Here we show how a mean field approach leads quite naturally two transport coefficients, not one, that govern mass and angular momentum transport. In particular, we highlight that the conventional approach suffers from a seemingly inconsistent neglect of turbulent diffusion in the surface density equation. We constrain these new transport coefficients for specific cases of inward, outward, and zero net mass transport. In addition, we find that one of the new transport terms can lead to oscillations in the mean surface density which then requires a constant or small inverse Rossby number for disks to maintain a monotonic power-law surface density.

preprint2008arXiv

Magnetic Fields in Paradigms of Planetary Nebulae and Related MHD Frontiers

Many, if not all, post AGB stellar systems swiftly transition from a spherical to a powerful aspherical pre-planetary nebula (pPNE) outflow phase before waning into a PNe. The pPNe outflows require engine rotational energy and a mechanism to extract this energy into collimated outflows. Just radiation and rotation are insufficient but a symbiosis between rotation, differential rotation and large scale magnetic fields remains promising. Present observational evidence for magnetic fields in evolved stars is suggestive of dynamically important magnetic fields, but both theory and observation are rife with research opportunity. I discuss how magnetohydrodynamic outflows might arise in pPNe and PNe and distinguish different between approaches that address shaping vs. those that address both launch and shaping. Scenarios involving dynamos in single stars, binary driven dynamos, or accretion engines cannot be ruled out. One appealing paradigm involves accretion onto the primary post-AGB white dwarf core from a low mass companion whose decaying accretion supply rate owers first the pPNe and then the lower luminosity PNe. Determining observational signatures of different MHD engines is a work in progress. Accretion disk theory and large scale dynamos pose many of their own fundamental challenges, some of which I discuss in a broader context.

preprint2008arXiv

Outflow Driven Turbulence in Molecular Clouds

In this paper we explore the relationship between protostellar outflows and turbulence in molecular clouds. Using 3-D numerical simulations we focus on the hydrodynamics of multiple outflows interacting within a parsec scale volume. We explore the extent to which transient outflows injecting directed energy and momentum into a sub-volume of a molecular cloud can be converted into random turbulent motions. We show that turbulence can readily be sustained by these interactions and show that it is possible to broadly characterize an effective driving scale of the outflows. We compare the velocity spectrum obtained in our studies to that of isotropically forced hydrodynamic turbulence finding that in outflow driven turbulence a power law is indeed achieved. However we find a steeper spectrum (beta ~ 3) is obtained in outflow driven turbulence models than in isotropically forced simulations (beta ~ 2). We discuss possible physical mechanisms responsible for these results as well and their implications for turbulence in molecular clouds where outflows will act in concert with other processes such as gravitational collapse.

preprint2005arXiv

Extracting Rotational Energy in Supernova Progenitors: Transient Poynting Flux Growth vs. Turbulent Dissipation

Observational evidence for anisotropy in supernovae (SN) may signal the importance of angular momentum and differential rotation in the progenitors. Free energy in differential rotation and rotation can be extracted magnetically or via turbulent dissipation. The importance that magnetohydrodyamic jets and coronae may play in driving SN motivates understanding large scale dynamos in SN progenitors. We develop a dynamical large scale interface dynamo model in which the differential rotation and rotation deplete both through Poynting flux and turbulent diffusion. We apply the model to a differentially rotating core surrounded by a convection zone of a SN progenitor from a initial 15$M_\odot$ star. Unlike the Sun, the dynamo is transient because the differential rotation is primarily due to the initial collapse. Up to $\sim 10^{51}$erg can be drained into time-integrated Poynting flux and heat, the relative fraction of which depends on the relative amount of turbulence in the shear layer vs. convection zone and the fraction of the shear layer into which the magnetic field penetrates. Both sinks can help facilitate explosions and could lead to different levels of anisotropy and pulsar kicks. In all cases, the poloidal magnetic field is much weaker than the toroidal field, and the Poynting flux is lower than previous estimates which invoke the magnitude of the total magnetic energy. A signature of a large scale dynamo is that the oscillation of the associated Poynting flux on $\sim 1$ sec time scales, implying the same for the energy delivery to a SN.

preprint2004arXiv

How Spectral Shapes of Magnetic Energy and Magnetic Helicity Influence their Respective Decay Time Scales

It is shown that in magnetohydrodynamics, the shapes of the magnetic energy and magnetic helicity spectra can influence whether the magnetic helicity resistive dissipation time scale is longer, equal to, or shorter than that of magnetic energy. The calculations highlight that magnetic helicity need not always dissipate significantly more slowly than magnetic energy. While this may be implicitly understood in detailed studies of magnetic helicity in plasma devices, it is not sufficiently emphasized in elementary discussions comparing magnetic helicity vs. magnetic energy evolution. Measuring magnetic energy and helicity decay times may provide a posteriori insight into the initial spectra.

preprint2001arXiv

Dynamos in Asymptotic-Giant-Branch Stars As the Origin of Magnetic Fields Shaping Planetary Nebulale

Planetary nebulae are thought to be formed when a slow wind from the progenitor giant star is overtaken by a subsequent fast wind generated as the star enters its white dwarf stage$^{1}$. A shock forms near the boundary between the winds, which creates a relatively dense shell that provides the characteristic appearance of a planetary nebula. A spherically symmetric wind will produce a spherically symmetric shell, yet over half of known planetary nebulae are not spherical; rather, they are elliptical or bipolar in shape$^{2}$. While a magnetic field could launch and collimate a bipolar outflow, the origin of such a field has hitherto been unclear, as previous work suggested that a field could not be generated${^3}$. Here we show that an asymptotic-giant-branch (AGB) star can indeed generate a strong magnetic field, in a dynamo at the interface between a rapidly rotating core and the more slowly rotating envelope of the star. The field is strong enough to shape the bipolar outflows that produce the observed bipolar planetary nebulae. Magnetic braking of the stellar core during this process may also explain the puzzlingly$^{4}$ slow rotation of most white dwarf stars.

preprint2000arXiv

Reprocessed emission line profiles from dense clouds in geometrically thick accretion engines

The central engines of active galactic nuclei (AGN) contain cold, dense material as well as hot X-ray emitting gas. The standard paradigm for the engine geometry is a cold thin disc sandwiched between hot X-ray coronae. Strong support for this geometry in Seyferts comes from the study of fluorescent iron line profiles, although the evidence is not ubiquitously air tight. The thin disc model of line profiles in AGN and in X-ray binaries should be bench marked against other plausible possibilities. One proposed alternative is an engine consisting of dense clouds embedded in an optically thin, geometrically thick X-ray emitting engine. This model is further motivated by studies of geometrically thick engines such as advection dominated accretion flows (ADAFs). Here we compute the reprocessed iron line profiles from dense clouds embedded in geometrically thick, optically thin X-ray emitting discs near a Schwarzchild black hole. We consider a range of cloud distributions and disc solutions, including ADAFs, pure radial infall, and bipolar outflows. We find that such models can reproduce line profiles similar to those from geometrically thin, optically thick discs and might help alleviate some of the problems encountered from the latter.

preprint1999arXiv

Reprocessed emission from warped accretion discs with application to X-ray iron line profiles

Flourescent iron line profiles currently provide the best diagnostic for active galactic nuclei (AGN) engine geometries. Here we construct a method for calculating the relativistic iron line profile from an arbitrarily warped accretion disc, illuminated from above and below by hard X-ray sources. This substantially generalises previous calculations of reprocessing by accretion discs by including non-axisymmetric effects. We include a relativistic treatment of shadowing by ray-tracing photon paths along Schwarzchild geodesics. We apply this method to two classes of warped discs, and generate a selection of resulting line profiles. New profile features include the possibility of sharper red, and softer blue fall-offs, a time varying line profile if the warp precesses about the disc, and some differences between `twisted' and `twist-free' warps. We discuss some qualitative implications of the line profiles in the context of Type I and II Seyfert AGN.

Eric G. Blackman

What is connected

Connect this record

See the researcher in context

Building this map preview

59 published item(s)

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

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

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

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

Successive common envelope events from multiple planets

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

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

Bipolar Planetary Nebulae from Outflow Collimation by Common Envelope Evolution

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

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

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

Generalized quenching of large-scale magnetic dynamos in anisotropic flows

Improving TBI Protection Measures and Standards for Combat Helmets

Necessity of ventilation for mitigating virus transmission quantified simply

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

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

Accretion in Common Envelope Evolution

Large scale dynamo action precedes turbulence in shearing box simulations of the magnetorotational instability

Vector cylindrical harmonics for low-dimensional convection models

Detecting the oldest geodynamo and attendant shielding from the solar wind: Implications for habitability

On the resilience of helical magnetic fields to turbulent diffusion and the astrophysical implications

Explaining the Observed Relation Between Stellar Activity and Rotation

Magnetic Helicity and Large Scale Magnetic Fields: A Primer

On characterizing nonlocality and anisotropy for the magnetorotational instability

Ribbons characterize magnetohydrodynamic magnetic fields better than lines: a lesson from dynamo theory

Sensitivity of the Magnetorotational Instability to the shear parameter in stratified simulations

Using Kinematic Properties of Pre-Planetary Nebulae to Constrain Engine Paradigms

A study of large scale dynamo growth rates from numerical simulations and implications for mean field theories

Mass loading and knot formation in AGN jets by stellar winds

MHD Shock-Clump Evolution with Self-Contained Magnetic Fields

On deriving flux freezing in magnetohydrodynamics by direct differentiation

Particle-in-cell simulations of particle energization via shock drift acceleration from low Mach number quasi-perpendicular shocks in solar flares

Resilience of helical fields to turbulent diffusion II: direct numerical simulations

Size of discs formed by wind accretion in binaries can be underestimated if the role of wind-driving force is ignored

The Formation and Evolution of Wind-Capture Disks In Binary Systems

Accretion Disks and Dynamos: Toward a Unified Mean Field Theory

Comparing Poynting flux dominated magnetic tower jets with kinetic-energy dominated jets

From Bipolar to Elliptical: Simulating the Morphological Evolution of Planetary Nebulae

Not Much Helicity is Needed to Drive Large Scale Dynamos

On the structure and stability of magnetic tower jets

Particle-in-cell simulations of particle energization from low Mach number fast mode shocks

Anderson and Stoner Published White Dwarf Mass Limits Before Chandrasekhar

Fast Shocks From Magnetic Reconnection Outflows

Modeling Poynting flux vs. kinetic-energy dominated jets

From Bipolar to Elliptical: Morphological Changes in the Temporal Evolution of PN

Quantifying the Imprecision of Accretion Theory and Implications for Multi-Epoch Observations of Protoplanetary Discs

Comparisons and Connections between Mean Field Dynamo Theory and Accretion Disc Theory

Coronae as Consequence of Large Scale Magnetic Fields in Turbulent Accretion Disks

Magnetic Nested-wind Scenarios for Bipolar Outflows: Pre-planetary and YSO nebular shaping

Stochastic Wobble of Accretion Discs and Jets from Turbulent Rocket Torques

A Search for Molecular Gas in the Nucleus of M87 and Implications for the Fueling of Supermassive Black Holes

Identifying Deficiencies of Standard Accretion Disk Theory: Lessons from a Mean-Field Approach

Magnetic Fields in Paradigms of Planetary Nebulae and Related MHD Frontiers

Outflow Driven Turbulence in Molecular Clouds

Extracting Rotational Energy in Supernova Progenitors: Transient Poynting Flux Growth vs. Turbulent Dissipation

How Spectral Shapes of Magnetic Energy and Magnetic Helicity Influence their Respective Decay Time Scales

Dynamos in Asymptotic-Giant-Branch Stars As the Origin of Magnetic Fields Shaping Planetary Nebulale

Reprocessed emission line profiles from dense clouds in geometrically thick accretion engines

Reprocessed emission from warped accretion discs with application to X-ray iron line profiles