Simulation of the Magnetothermal Instability

In many magnetized, dilute astrophysical plasmas, thermal conduction occurs almost exclusively parallel to magnetic field lines. In this case, the usual stability criterion for convective stability, the Schwarzschild criterion, which depends on entropy gradients, is modified. In the magnetized long mean free path regime, instability occurs for small wavenumbers when (dP/dz)(dln T/dz) > 0, which we refer to as the Balbus criterion. We refer to the convective-type instability that results as the magnetothermal instability (MTI). We use the equations of MHD with anisotropic electron heat conduction to numerically simulate the linear growth and nonlinear saturation of the MTI in plane-parallel atmospheres that are unstable according to the Balbus criterion. The linear growth rates measured from the simulations are in excellent agreement with the weak field dispersion relation. The addition of isotropic conduction, e.g. radiation, or strong magnetic fields can damp the growth of the MTI and affect the nonlinear regime. The instability saturates when the atmosphere becomes isothermal as the source of free energy is exhausted. By maintaining a fixed temperature difference between the top and bottom boundaries of the simulation domain, sustained convective turbulence can be driven. MTI-stable layers introduced by isotropic conduction are used to prevent the formation of unresolved, thermal boundary layers. We find that the largest component of the time-averaged heat flux is due to advective motions as opposed to the actual thermal conduction itself. Finally, we explore the implications of this instability for a variety of astrophysical systems, such as neutron stars, the hot intracluster medium of galaxy clusters, and the structure of radiatively inefficient accretion flows.Comment: Accepted for publication in Astrophysics and Space Science as proceedings of the 6th High Energy Density Laboratory Astrophysics (HEDLA) Conferenc

Buoyancy driven rotating boundary currents

The structure of boundary currents formed from intermediately dense water introduced into a rotating, stably stratified, two-layer environment is investigated in a series of laboratory experiments, performed for Froude numbers ranging from 0.01 to 1. The thickness and streamwise velocity profiles in quasi-steady currents are measured using a pH activated tracer (thymol blue) and found to compare favorably to simplified analytic solutions and numerical models. Currents flowing along sloping boundaries in a stratified background exhibit robust stability at all experimental Froude numbers. Such stability is in sharp contrast to the unequivocal instability of such currents flowing against vertical boundaries, or of currents flowing along slopes in a uniform background. The presence of a variety of wave mechanisms in the ambient medium might account for the slower and wider observed structures and the stability of the currents, by effecting the damping of disturbances through wave radiation.Comment: 9 pages with 2 figures to appear in Ann NYAS "Long range effects in physics and astrophysics

The role of Hall diffusion in the magnetically threaded thin accretion discs

We study role of the Hall diffusion in the magnetic star-disc interaction. In a simplified steady state configuration, the total torque is calculated in terms of the fastness parameter and a new term because of the Hall diffusion. We show the total torque reduces as the Hall term becomes more significant. Also, the critical fastness parameter (at which the total torque is zero) reduces because of the Hall diffusion.Comment: Accepted for publication in Astrophysics and Space Scienc

Magnetic fields in protoplanetary disks

Magnetic fields likely play a key role in the dynamics and evolution of protoplanetary discs. They have the potential to efficiently transport angular momentum by MHD turbulence or via the magnetocentrifugal acceleration of outflows from the disk surface, and magnetically-driven mixing has implications for disk chemistry and evolution of the grain population. However, the weak ionisation of protoplanetary discs means that magnetic fields may not be able to effectively couple to the matter. I present calculations of the ionisation equilibrium and magnetic diffusivity as a function of height from the disk midplane at radii of 1 and 5 AU. Dust grains tend to suppress magnetic coupling by soaking up electrons and ions from the gas phase and reducing the conductivity of the gas by many orders of magnitude. However, once grains have grown to a few microns in size their effect starts to wane and magnetic fields can begin to couple to the gas even at the disk midplane. Because ions are generally decoupled from the magnetic field by neutral collisions while electrons are not, the Hall effect tends to dominate the diffusion of the magnetic field when it is able to partially couple to the gas. For a standard population of 0.1 micron grains the active surface layers have a combined column of about 2 g/cm^2 at 1 AU; by the time grains have aggregated to 3 microns the active surface density is 80 g/cm^2. In the absence of grains, x-rays maintain magnetic coupling to 10% of the disk material at 1 AU (150 g/cm^2). At 5 AU the entire disk thickness becomes active once grains have aggregated to 1 micron in size.Comment: 11 pages, 11 figs, aastex.cls. Accepted for publication in Astrophysics & Space Science. v3 corrects bibliograph

Density probability distribution in one-dimensional polytropic gas dynamics

We discuss the generation and statistics of the density fluctuations in highly compressible polytropic turbulence, based on a simple model and one-dimensional numerical simulations. Observing that density structures tend to form in a hierarchical manner, we assume that density fluctuations follow a random multiplicative process. When the polytropic exponent $\gamma$ is equal to unity, the local Mach number is independent of the density, and our assumption leads us to expect that the probability density function (PDF) of the density field is a lognormal. This isothermal case is found to be singular, with a dispersion $\sigma_s^2$ which scales like the square turbulent Mach number $\tilde M^2$, where $s\equiv \ln \rho$ and $\rho$ is the fluid density. This leads to much higher fluctuations than those due to shock jump relations. Extrapolating the model to the case $\gamma \not =1$, we find that, as the Mach number becomes large, the density PDF is expected to asymptotically approach a power-law regime, at high densities when $\gamma<1$, and at low densities when $\gamma>1$. This effect can be traced back to the fact that the pressure term in the momentum equation varies exponentially with $s$, thus opposing the growth of fluctuations on one side of the PDF, while being negligible on the other side. This also causes the dispersion $\sigma_s^2$ to grow more slowly than $\tilde M^2$ when $\gamma\not=1$. In view of these results, we suggest that Burgers flow is a singular case not approached by the high-$\tilde M$ limit, with a PDF that develops power laws on both sides.Comment: 9 pages + 12 postscript figures. Submitted to Phys. Rev.

Plane shearing waves of arbitrary form: exact solutions of the Navier-Stokes equations

We present exact solutions of the incompressible Navier-Stokes equations in a background linear shear flow. The method of construction is based on Kelvin's investigations into linearized disturbances in an unbounded Couette flow. We obtain explicit formulae for all three components of a Kelvin mode in terms of elementary functions. We then prove that Kelvin modes with parallel (though time-dependent) wave vectors can be superposed to construct the most general plane transverse shearing wave. An explicit solution is given, with any specified initial orientation, profile and polarization structure, with either unbounded or shear-periodic boundary conditions.Comment: 6 pages, 2 figures; version published in the European Physical Journal Plu

Multicomponent theory of buoyancy instabilities in magnetized plasmas: The case of magnetic field parallel to gravity

We investigate electromagnetic buoyancy instabilities of the electron-ion plasma with the heat flux based on not the magnetohydrodynamic (MHD) equations, but using the multicomponent plasma approach when the momentum equations are solved for each species. We consider a geometry in which the background magnetic field, gravity, and stratification are directed along one axis. The nonzero background electron thermal flux is taken into account. Collisions between electrons and ions are included in the momentum equations. No simplifications usual for the one-fluid MHD-approach in studying these instabilities are used. We derive a simple dispersion relation, which shows that the thermal flux perturbation generally stabilizes an instability for the geometry under consideration. This result contradicts to conclusion obtained in the MHD-approach. We show that the reason of this contradiction is the simplified assumptions used in the MHD analysis of buoyancy instabilities and the role of the longitudinal electric field perturbation which is not captured by the ideal MHD equations. Our dispersion relation also shows that the medium with the electron thermal flux can be unstable, if the temperature gradients of ions and electrons have the opposite signs. The results obtained can be applied to the weakly collisional magnetized plasma objects in laboratory and astrophysics.Comment: Accepted for publication in Astrophysics & Space Scienc

Are the magnetic fields of millisecond pulsars ~ 10^8 G?

It is generally assumed that the magnetic fields of millisecond pulsars (MSPs) are $\sim 10^{8}$G. We argue that this may not be true and the fields may be appreciably greater. We present six evidences for this: (1) The $\sim 10^{8}$ G field estimate is based on magnetic dipole emission losses which is shown to be questionable; (2) The MSPs in low mass X-ray binaries (LMXBs) are claimed to have $< 10^{11}$ G on the basis of a Rayleygh-Taylor instability accretion argument. We show that the accretion argument is questionable and the upper limit $10^{11}$ G may be much higher; (3) Low magnetic field neutron stars have difficulty being produced in LMXBs; (4) MSPs may still be accreting indicating a much higher magnetic field; (5) The data that predict $\sim 10^{8}$ G for MSPs also predict ages on the order of, and greater than, ten billion years, which is much greater than normal pulsars. If the predicted ages are wrong, most likely the predicted $\sim 10^{8}$ G fields of MSPs are wrong; (6) When magnetic fields are measured directly with cyclotron lines in X-ray binaries, fields $\gg 10^{8}$ G are indicated. Other scenarios should be investigated. One such scenario is the following. Over 85% of MSPs are confirmed members of a binary. It is possible that all MSPs are in large separation binaries having magnetic fields $> 10^{8}$ G with their magnetic dipole emission being balanced by low level accretion from their companions.Comment: 16 pages, accept for publication in Astrophysics and Space Scienc

Resistive and magnetized accretion flows with convection

We considered the effects of convection on the radiatively inefficient accretion flows (RIAF) in the presence of resistivity and toroidal magnetic field. We discussed the effects of convection on transports of angular momentum and energy. We established two cases for the resistive and magnetized RIAFs with convection: assuming the convection parameter as a free parameter and using mixing-length theory to calculate convection parameter. A self-similar method was used to solve the integrated equations that govern the behavior of the presented model. The solutions showed that the accretion and rotational velocities decrease by adding the convection parameter, while the sound speed increases. Moreover, by using mixing-length theory to calculate convection parameter, we found that the convection can be important in RIAFs with magnetic field and resistivity.Comment: 7 pages, 3 figures, accepted by Ap&S

Radial and vertical angular momentum transport in protostellar discs

Angular momentum in protostellar discs can be transported either radially, through turbulence induced by the magnetorotational instability (MRI), or vertically, through the torque exerted by a large-scale magnetic field. We present a model of steady-state discs where these two mechanisms operate at the same radius and derive approximate criteria for their occurrence in an ambipolar diffusion dominated disc. We obtain "weak field'' solutions - which we associate with the MRI channel modes in a stratified disc - and transform them into accretion solutions with predominantly radial angular-momentum transport by implementing a turbulent-stress prescription based on published results of numerical simulations. We also analyze "intermediate field strength'' solutions in which both radial and vertical transport operate at the same radial location. Our results suggest, however, that this overlap is unlikely to occur in real discs.Comment: 5 pages, 2 figures, 1 table, aastex.cls. Accepted for publication in Astrophysics & Space Scienc