BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Matthew W. Kunz

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

ResearcherUnclaimed source record
astro-ph.SRastro-ph.HEphysics.plasm-phastro-ph.EPastro-ph.GAastro-ph.COphysics.space-phastro-phphysics.comp-phphysics.flu-dyn

Catalog footprint

What is connected

20works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

20 published item(s)

preprint2022arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The vertical shear instability in poorly ionised, magnetized protoplanetary discs

Protoplanetary discs should exhibit a weak vertical variation in their rotation profiles. Typically this `vertical shear' issues from a baroclinic effect driven by the central star's radiation field, but it might also arise during the launching of a magnetocentrifugal wind. As a consequence, protoplanetary discs are subject to a hydrodynamical instability, the `vertical shear instability' (VSI), whose breakdown into turbulence could transport a moderate amount of angular momentum and facilitate, or interfere with, the process of planet formation. Magnetic fields may suppress the VSI, however, either directly via magnetic tension or indirectly through magnetorotational turbulence. On the other hand, protoplanetary discs exhibit notoriously low ionisation fractions, and non-ideal effects, if sufficiently dominant, may come to the VSI's rescue. In this paper we develop a local linear theory that explores how non-ideal MHD influences the VSI, while also launching additional diffusive shear instabilities. We derive a set of analytical criteria that establish when the VSI prevails, and then show how it can be applied to a realistic global model of a protoplanetary disc. Our calculations suggest that within ~10au the VSI should have little trouble emerging in the main body of the disk, but beyond that, and in the upper regions of the disc, its onset depends sensitively on the size of the preponderant dust grains.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Triggering tearing in a forming current sheet with the mirror instability

We study the time-dependent formation and evolution of a current sheet (CS) in a magnetized, collisionless, high-beta plasma using hybrid-kinetic particle-in-cell simulations. An initially tearing-stable Harris sheet is frozen into a persistently driven incompressible flow so that its characteristic thickness gradually decreases in time. As the CS thins, the strength of the reconnecting field increases, and adiabatic invariance in the inflowing fluid elements produces a field-biased pressure anisotropy with excess perpendicular pressure. At large plasma beta, this anisotropy excites the mirror instability, which deforms the reconnecting field on ion-Larmor scales and dramatically reduces the effective thickness of the CS. Tearing modes whose wavelengths are comparable to that of the mirrors then become unstable, triggering reconnection on smaller scales and at earlier times than would have occurred if the thinning CS were to have retained its Harris profile. A novel method for identifying and tracking X-points is introduced, yielding X-point separations that are initially intermediate between the perpendicular and parallel mirror wavelengths in the upstream plasma. These mirror-stimulated tearing modes ultimately grow and merge to produce island widths comparable to the CS thickness, an outcome we verify across a range of CS formation timescales and initial CS widths. Our results may find their most immediate application in the tearing disruption of magnetic folds generated by turbulent dynamo in weakly collisional, high-beta, astrophysical plasmas.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Formation and evolution of protostellar accretion discs. I. Angular-momentum budget, gravitational self-regulation, and numerical convergence

We investigate the formation and early evolution of a protostellar disc from a magnetized pre-stellar core using non-ideal magnetohydrodynamic (MHD) simulations including ambipolar diffusion and Ohmic dissipation. The dynamical contraction of the pre-stellar core ultimately leads to the formation of a first hydrostatic core, after ambipolar diffusion decouples the magnetic field from the predominantly neutral gas. The hydrostatic core accumulates angular momentum from the infalling material, evolving into a rotationally supported torus; this `first hydrostatic torus' then forms an accreting protostar and a rotationally supported disc. The disc spreads out by gravitational instability, reaching $\sim$30 au in diameter at $\sim$3 kyr after protostar formation. The total mass and angular momentum of the protostar-disc system are determined mainly by accretion of gas from an infalling pseudo-disc, which has low specific angular momentum because of magnetic braking; their removal from the protostar-disc system by outflow and disc magnetic braking are negligible, in part because the magnetic field is poorly coupled there. The redistribution of angular momentum within the protostar-disc system is facilitated mainly by gravitational instability; this allows formation of relatively large discs even when the specific angular momentum of infalling material is low. We argue that such discs should remain marginally unstable as they grow (with Toomre $Q\sim 1$-$2$), an idea that is broadly consistent with recent observational estimates for Class 0/I discs. We discuss the numerical convergence of our results, and show that properly treating the inner boundary condition is crucial for achieving convergence at an acceptable computational cost.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Production and persistence of extreme two-temperature plasmas in radiative relativistic turbulence

Turbulence is a predominant process for energizing electrons and ions in collisionless astrophysical plasmas, and thus is responsible for shaping their radiative signatures (luminosity, spectra, and variability). To better understand the kinetic properties of a collisionless radiative plasma subject to externally driven turbulence, we investigate particle-in-cell simulations of relativistic plasma turbulence with external inverse Compton cooling acting on the electrons. We find that ions continuously heat up while electrons gradually cool down (due to the net effect of radiation), and hence the ion-to-electron temperature ratio $T_i/T_e$ grows in time. We show that $T_i/T_e$ is limited only by the size and duration of the simulations (reaching $T_i/T_e \sim 10^3$), indicating that there are no efficient collisionless mechanisms of electron-ion thermal coupling. This result has implications for models of radiatively inefficient accretion flows, such as observed in the Galactic Center and in M87, for which so-called two-temperature plasmas with $T_i/T_e \gg 1$ have been invoked to explain their low luminosity. Additionally, we find that electrons acquire a quasi-thermal distribution (dictated by the competition of turbulent particle energization and radiative cooling), while ions undergo efficient nonthermal acceleration (acquiring a harder distribution than in equivalent non-radiative simulations). There is a modest nonthermal population of high-energy electrons that are beamed intermittently in space, time, and direction; these beamed electrons may explain rapid flares in certain high-energy astrophysical systems (e.g., in the Galactic Center). These numerical results demonstrate that extreme two-temperature plasmas can be produced and maintained by relativistic radiative turbulence.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Linear Vlasov theory of a magnetised, thermally stratified atmosphere

The stability of a collisionless, magnetised plasma to local convective disturbances is examined, with a focus on kinetic and finite-Larmor-radius effects. Specific application is made to the outskirts of galaxy clusters, which contain hot and tenuous plasma whose temperature increases in the direction of gravity. At long wavelengths (the "drift-kinetic" limit), we obtain the kinetic version of the magnetothermal instability (MTI) and its Alfvénic counterpart (Alfvénic MTI), which were previously discovered and analysed using a magnetofluid (i.e. Braginskii) description. At sub-ion-Larmor scales, we discover an overstability driven by the electron temperature gradient of kinetic-Alfvén drift waves -- the electron MTI (eMTI) -- whose growth rate is even larger than the standard MTI. At intermediate scales, we find that ion finite-Larmor-radius effects tend to stabilise the plasma. We discuss the physical interpretation of these instabilities in detail, and compare them both with previous work on magnetised convection in a collisional plasma and with temperature-gradient-driven drift-wave instabilities well-known to the magnetic-confinement-fusion community. The implications of having both fluid and kinetic scales simultaneously driven unstable by the same temperature gradient are briefly discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Magnetorotational Turbulence and Dynamo in a Collisionless Plasma

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Magnetically driven accretion in protoplanetary discs

We characterize magnetically driven accretion at radii between 1 au and 100 au in protoplanetary discs, using a series of local non-ideal magnetohydrodynamic (MHD) simulations. The simulations assume a Minimum Mass Solar Nebula (MMSN) disc that is threaded by a net vertical magnetic field of specified strength. Confirming previous results, we find that the Hall effect has only a modest impact on accretion at 30 au, and essentially none at 100 au. At 1-10 au the Hall effect introduces a pronounced bi-modality in the accretion process, with vertical magnetic fields aligned to the disc rotation supporting a strong laminar Maxwell stress that is absent if the field is anti-aligned. In the anti-aligned case, we instead find evidence for bursts of turbulent stress at 5-10 au, which we tentatively identify with the non-axisymmetric Hall-shear instability. The presence or absence of these bursts depends upon the details of the adopted chemical model, which suggests that appreciable regions of actual protoplanetary discs might lie close to the borderline between laminar and turbulent behaviour. Given the number of important control parameters that have already been identified in MHD models, quantitative predictions for disc structure in terms of only radius and accretion rate appear to be difficult. Instead, we identify robust qualitative tests of magnetically driven accretion. These include the presence of turbulence in the outer disc, independent of the orientation of the vertical magnetic fields, and a Hall-mediated bi-modality in turbulent properties extending from the region of thermal ionization to 10 au.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Firehose and Mirror Instabilities in a Collisionless Shearing Plasma

Hybrid-kinetic numerical simulations of firehose and mirror instabilities in a collisionless plasma are performed in which pressure anisotropy is driven as the magnetic field is changed by a persistent linear shear $S$. For a decreasing field, it is found that mostly oblique firehose fluctuations grow at ion Larmor scales and saturate with energies $\sim$$S^{1/2}$; the pressure anisotropy is pinned at the stability threshold by particle scattering off microscale fluctuations. In contrast, nonlinear mirror fluctuations are large compared to the ion Larmor scale and grow secularly in time; marginality is maintained by an increasing population of resonant particles trapped in magnetic mirrors. After one shear time, saturated order-unity magnetic mirrors are formed and particles scatter off their sharp edges. Both instabilities drive sub-ion-Larmor--scale fluctuations, which appear to be kinetic-Alfvén-wave turbulence. Our results impact theories of momentum and heat transport in astrophysical and space plasmas, in which the stretching of a magnetic field by shear is a generic process.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Thanatology in Protoplanetary Discs: the combined influence of Ohmic, Hall, and ambipolar diffusion on dead zones

Protoplanetary discs are poorly ionised due to their low temperatures and high column densities, and are therefore subject to three "non-ideal" magnetohydrodynamic effects: Ohmic dissipation, ambipolar diffusion, and the Hall effect. The existence of magnetically driven turbulence in these discs has been a central question since the discovery of the magnetorotational instability. Early models considered Ohmic diffusion only and led to a scenario of layered accretion, in which a magnetically "dead" zone in the disc midplane is embedded within magnetically "active" surface layers at distances ~1-10 au from the central protostellar object. Recent work has suggested that a combination of Ohmic dissipation and ambipolar diffusion can render both the midplane and surface layers of the disc inactive and that torques due to magnetically driven outflows are required to explain the observed accretion rates. We reassess this picture by performing three-dimensional numerical simulations that include, for the first time, all three non-ideal MHD effects. We find that the Hall effect can generically "revive" dead zones by producing a dominant azimuthal magnetic field and a large-scale Maxwell stress throughout the midplane, provided the angular velocity and magnetic field satisfy Omega.B > 0. The attendant large magnetic pressure modifies the vertical density profile and substantially increases the disc scale height beyond its hydrostatic value. Outflows are produced, but are not necessary to explain accretion rates <10^{-7} Msun/yr. The flow in the disc midplane is essentially laminar, suggesting that dust sedimentation may be efficient. These results demonstrate that, if the MRI is relevant for driving mass accretion in protoplanetary discs, one must include the Hall effect to obtain even qualitatively correct results.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Magnetic self-organisation in Hall-dominated magnetorotational turbulence

The magnetorotational instability (MRI) is the most promising mechanism by which angular momentum is efficiently transported outwards in astrophysical discs. However, its application to protoplanetary discs remains problematic. These discs are so poorly ionised that they may not support magnetorotational turbulence in regions referred to as `dead zones'. It has recently been suggested that the Hall effect, a non-ideal magnetohydrodynamic (MHD) effect, could revive these dead zones by enhancing the magnetically active column density by an order of magnitude or more. We investigate this idea by performing local, three-dimensional, resistive Hall-MHD simulations of the MRI in situations where the Hall effect dominates over Ohmic dissipation. As expected from linear stability analysis, we find an exponentially growing instability in regimes otherwise linearly stable in resistive MHD. However, instead of vigorous and sustained magnetorotational turbulence, we find that the MRI saturates by producing large-scale, long-lived, axisymmetric structures in the magnetic and velocity fields. We refer to these structures as zonal fields and zonal flows, respectively. Their emergence causes a steep reduction in turbulent transport by at least two orders of magnitude from extrapolations based upon resistive MHD, a result that calls into question contemporary models of layered accretion. We construct a rigorous mean-field theory to explain this new behaviour and to predict when it should occur. Implications for protoplanetary disc structure and evolution, as well as for theories of planet formation, are briefly discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Pegasus: A New Hybrid-Kinetic Particle-in-Cell Code for Astrophysical Plasma Dynamics

We describe Pegasus, a new hybrid-kinetic particle-in-cell code tailored for the study of astrophysical plasma dynamics. The code incorporates an energy-conserving particle integrator into a stable, second-order--accurate, three-stage predictor-predictor-corrector integration algorithm. The constrained transport method is used to enforce the divergence-free constraint on the magnetic field. A delta-f scheme is included to facilitate a reduced-noise study of systems in which only small departures from an initial distribution function are anticipated. The effects of rotation and shear are implemented through the shearing-sheet formalism with orbital advection. These algorithms are embedded within an architecture similar to that used in the popular astrophysical magnetohydrodynamics code Athena, one that is modular, well-documented, easy to use, and efficiently parallelized for use on thousands of processors. We present a series of tests in one, two, and three spatial dimensions that demonstrate the fidelity and versatility of the code.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Bridging the gap: disk formation in the Class 0 phase with ambipolar diffusion and Ohmic dissipation

Context: Ideal MHD simulations have revealed catastrophic magnetic braking (MB) in the protostellar phase, which prevents the formation of a centrifugal disk around a nascent protostar. Aims: We determine if non-ideal MHD, including the effects of ambipolar diffusion and Ohmic dissipation determined from a detailed chemical network model, allows for disk formation at the earliest stages of star formation (SF). Methods: We employ the axisymmetric thin-disk approximation in order to resolve a dynamic range of 9 orders of magnitude in length and 16 in density, while also calculating partial ionization using up to 19 species in a detailed chemical equilibrium model. MB is applied using a steady-state approximation, and a barotropic relation is used to capture the thermal evolution. Results: We resolve the formation of the first and second cores, with expansion waves at the periphery of each, a magnetic diffusion shock, and prestellar infall profiles at larger radii. Power-law profiles in each region can be understood analytically. After the formation of the second core, centrifugal support rises rapidly and a low-mass disk of radius ~10 R_Sun is formed, when the second core has mass ~0.001 M_Sun. The mass-to-flux ratio is ~10,000 times the critical value in the central region. Conclusions: A small centrifugal disk can form in the earliest stage of SF, due to a shut-off of MB caused by magnetic field dissipation in the first core region. There is enough angular momentum loss to allow the second collapse to occur directly, and a low-mass stellar core to form with a surrounding disk. The disk mass and size will depend upon how the angular momentum transport mechanisms within the disk can keep up with mass infall onto the disk. We estimate that the disk will remain <~10 AU, undetectable even by ALMA, in the early Class 0 phase.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Buoyancy Instabilities in a Weakly Collisional Intracluster Medium

The intracluster medium of galaxy clusters is a weakly collisional, high-beta plasma in which the transport of heat and momentum occurs primarily along magnetic-field lines. Anisotropic heat conduction allows convective instabilities to be driven by temperature gradients of either sign, the magnetothermal instability (MTI) in the outskirts of non-isothermal clusters and the heat-flux buoyancy-driven instability (HBI) in their cooling cores. We employ the Athena MHD code to investigate the nonlinear evolution of these instabilities, self-consistently including the effects of anisotropic viscosity (i.e. Braginskii pressure anisotropy), anisotropic conduction, and radiative cooling. We highlight the importance of the microscale instabilities that inevitably accompany and regulate the pressure anisotropies generated by the HBI and MTI. We find that, in all but the innermost regions of cool-core clusters, anisotropic viscosity significantly impairs the ability of the HBI to reorient magnetic-field lines orthogonal to the temperature gradient. Thus, while radio-mode feedback appears necessary in the central few tens of kpc, conduction may be capable of offsetting radiative losses throughout most of a cool core over a significant fraction of the Hubble time. Magnetically-aligned cold filaments are then able to form by local thermal instability. Viscous dissipation during the formation of a cold filament produces accompanying hot filaments, which can be searched for in deep Chandra observations of nearby cool-core clusters. In the case of the MTI, anisotropic viscosity maintains the coherence of magnetic-field lines over larger distances than in the inviscid case, providing a natural lower limit for the scale on which the field can fluctuate freely. In the nonlinear state, the magnetic field exhibits a folded structure in which the field-line curvature and field strength are anti-correlated.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The HBI in a quasi-global model of the intracluster medium

In this paper we investigate how convective instabilities influence heat conduction in the intracluster medium (ICM) of cool-core galaxy clusters. The ICM is a high-beta, weakly collisional plasma in which the transport of momentum and heat is aligned with the magnetic field. The anisotropy of heat conduction, in particular, gives rise to instabilities that can access energy stored in a temperature gradient of either sign. We focus on the heat-flux buoyancy-driven instability (HBI), which feeds on the outwardly increasing temperature profile of cluster cool cores. Our aim is to elucidate how the global structure of a cluster impacts on the growth and morphology of the linear HBI modes when in the presence of Braginskii viscosity, and ultimately on the ability of the HBI to thermally insulate cores. We employ an idealised quasi-global model, the plane-parallel atmosphere, which captures the essential physics -- e.g. the global radial profile of the cluster -- while letting the problem remain analytically tractable. Our main result is that the dominant HBI modes are localised to the the innermost (~<20%) regions of cool cores. It is then probable that, in the nonlinear regime, appreciable field-line insulation will be similarly localised. Thus, while radio-mode feedback appears necessary in the central few tens of kpc, heat conduction may be capable of offsetting radiative losses throughout most of a cool core over a significant fraction of the Hubble time. Finally, our linear solutions provide a convenient numerical test for the nonlinear codes that tackle the saturation of such convective instabilities in the presence of anisotropic transport.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dynamical stability of a thermally stratified intracluster medium with anisotropic momentum and heat transport

In weakly-collisional plasmas such as the intracluster medium (ICM), heat and momentum transport become anisotropic with respect to the local magnetic field direction. Anisotropic heat conduction causes the slow magnetosonic wave to become buoyantly unstable to the magnetothermal instability (MTI) when the temperature increases in the direction of gravity and to the heat-flux--driven buoyancy instability (HBI) when the temperature decreases in the direction of gravity. The local changes in magnetic field strength that attend these instabilities cause pressure anisotropies that viscously damp motions parallel to the magnetic field. In this paper we employ a linear stability analysis to elucidate the effects of anisotropic viscosity (i.e. Braginskii pressure anisotropy) on the MTI and HBI. By stifling the convergence/divergence of magnetic field lines, pressure anisotropy significantly affects how the ICM interacts with the temperature gradient. Instabilities which depend upon the convergence/divergence of magnetic field lines to generate unstable buoyant motions (the HBI) are suppressed over much of the wavenumber space, whereas those which are otherwise impeded by field-line convergence/divergence (the MTI) are strengthened. As a result, the wavenumbers at which the HBI survives largely unsuppressed in the ICM have parallel components too small to rigorously be considered local. This is particularly true as the magnetic field becomes more and more orthogonal to the temperature gradient. In contrast, the fastest-growing MTI modes are unaffected by anisotropic viscosity. However, we find that anisotropic viscosity couples slow and Alfven waves in such a way as to buoyantly destabilise Alfvenic fluctuations when the temperature increases in the direction of gravity. Consequently, many wavenumbers previously considered MTI-stable or slow-growing are in fact maximally unstable. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Magnetic Field Diffusion and the Formation of Circumstellar Disks

A non-ideal MHD collapse calculation employing the axisymmetric thin-disk approximation is used to resolve cloud core collapse down to the scales of the second (stellar) core. Rotation and a magnetic braking torque are included in the model, and the partial ionization resulting in ambipolar diffusion and Ohmic dissipation is calculated from a detailed chemical network. We find that a centrifugal disk can indeed form in the earliest stage of star formation, due to a shut-off of magnetic braking caused by magnetic field diffusion in the first core region. Thus, there is no catastrophic magnetic braking in a model with realistic non-ideal MHD.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Formation of Interstellar Clouds: Parker Instability with Phase Transitions

We follow numerically the nonlinear evolution of the Parker instability in the presence of phase transitions from a warm to a cold HI interstellar medium in two spatial dimensions. The nonlinear evolution of the system favors modes that allow the magnetic field lines to cross the galactic plane. Cold HI clouds form with typical masses ~= 10^5 M_sun, mean densities ~= 20 cm^-3, mean magnetic field strengths ~= 4.3 muG (rms field strengths ~= 6.4 muG), mass-to-flux ratios ~= 0.1 - 0.3 relative to critical, temperatures ~= 50 K, (two-dimensional) turbulent velocity dispersions ~= 1.6 km s^-1, and separations ~= 500 pc, in agreement with observations. The maximum density and magnetic field strength are ~= 10^3 cm^-3 and ~= 20 muG, respectively. Approximately 60% of all HI mass is in the warm neutral medium. The cold neutral medium is arranged into sheet-like structures both perpendicular and parallel to the galactic plane, but it is also found almost everywhere in the galactic plane, with the density being highest in valleys of the magnetic field lines. `Cloudlets' also form whose physical properties are in quantitative agreement with those observed for such objects by Heiles (1967). The nonlinear phase of the evolution takes ~< 30 Myr, so that, if the instability is triggered by a nonlinear perturbation such as a spiral density shock wave, interstellar clouds can form within a time suggested by observations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Initial Core Mass Function due to Ambipolar Diffusion in Molecular Clouds

We show that the ambipolar-diffusion--initiated fragmentation of molecular clouds leads simply and naturally to an initial core mass function (CMF) which is very similar to the initial stellar mass function (IMF) and in excellent agreement with existing observations. This agreement is robust provided that the three (input) free parameters remain within their range of values suggested by observations. Other, observationally testable, predictions are made.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

The Nonisothermal Stage of Magnetic Star Formation. I. Formulation of the Problem and Method of Solution

We formulate the problem of the formation and subsequent evolution of fragments (or cores) in magnetically-supported, self-gravitating molecular clouds in two spatial dimensions. The six-fluid (neutrals, electrons, molecular and atomic ions, positively-charged, negatively-charged, and neutral grains) physical system is governed by the radiative, nonideal magnetohydrodynamic (RMHD) equations. The magnetic flux is not assumed to be frozen in any of the charged species. Its evolution is determined by a newly-derived generalized Ohm's law, which accounts for the contributions of both elastic and inelastic collisions to ambipolar diffusion and Ohmic dissipation. The species abundances are calculated using an extensive chemical-equilibrium network. Both MRN and uniform grain size distributions are considered. The thermal evolution of the protostellar core and its effect on the dynamics are followed by employing the grey flux-limited diffusion approximation. Realistic temperature-dependent grain opacities are used that account for a variety of grain compositions. We have augmented the publicly-available Zeus-MP code to take into consideration all these effects and have modified several of its algorithms to improve convergence, accuracy and efficiency. Results of magnetic star formation simulations that accurately track the evolution of a protostellar fragment from a density ~10^3 cm^-3 to a density ~10^15 cm^-3, while rigorously accounting for both nonideal MHD processes and radiative transfer, are presented in a separate paper.

0 citations0 reviewsNo code linkedOpen access