Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userThe Effect of Cosmic Ray Diffusion on the Parker Instability
57
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The Parker instability, which has been considered as a process governing the structure of the interstellar medium, is induced by the buoyancy of magnetic field and cosmic rays. In previous studies, while the magnetic field has been fully incorporated in the context of isothermal magnetohydrodynamics, cosmic rays have been normally treated with the simplifying assumption of infinite diffusion along magnetic field lines but no diffusion across them. The cosmic ray diffusion is, however, finite. In this work, we take into account fully the diffusion process of cosmic rays in a linear stability analysis of the Parker instability. Cosmic rays are described with the diffusion-convection equation. With realistic values of cosmic ray diffusion coefficients expected in the interstellar medium, we show that the result of previous studies with the simplifying assumption on cosmic ray diffusion applies well. Finiteness of parallel diffusion decreases the growth rate of the Parker instability, while the relatively smaller perpendicular diffusion has no significant effect. We discuss the implication of our result on the role of the Parker instability in the interstellar medium.
rate research
Read More
We present the results obtained from linear stability analysis and 2.5-dimensional magnetohydrodynamic (MHD) simulations of the magnetorotational instability (MRI), including the effects of cosmic rays (CRs). We took into account of the CR diffusion along the magnetic field but neglect the cross-field-line diffusion. Two models are considered in this paper: shearing box model and differentially rotating cylinder model. We studied how MRI is affected by the initial CR pressure (i.e., energy) distribution. In the shearing box model, the initial state is uniform distribution. Linear analysis shows that the growth rate of MRI does not depend on the value of CR diffusion coefficient. In the differentially rotating cylinder model, the initial state is a constant angular momentum polytropic disk threaded by weak uniform vertical magnetic field. Linear analysis shows that the growth rate of MRI becomes larger if the CR diffusion coefficient is larger. Both results are confirmed by MHD simulations. The MHD simulation results show that the outward movement of matter by the growth of MRI is not impeded by the CR pressure gradient, and the centrifugal force which acts to the concentrated matter becomes larger. Consequently, the growth rate of MRI is increased. On the other hand, if the initial CR pressure is uniform, then the growth rate of the MRI barely depends on the value of the CR diffusion coefficient.
We examine the evolution of the Parker instability in galactic disks using 3D numerical simulations. We consider a local Cartesian box section of a galactic disk, where gas, magnetic fields and cosmic rays are all initially in a magnetohydrostatic equilibrium. This is done for different choices of initial cosmic ray density and magnetic field. The growth rates and characteristic scales obtained from the models, as well as their dependences on the density of cosmic rays and magnetic fields, are in broad agreement with previous (linearized, ideal) analytical work. However, this non-ideal instability develops a multi-modal 3D structure, which cannot be quantitatively predicted from the earlier linearized studies. This 3D signature of the instability will be of importance in interpreting observations. As a preliminary step towards such interpretations, we calculate synthetic polarized intensity and Faraday rotation measure maps, and the associated structure functions of the latter, from our simulations; these suggest that the correlation scales inferred from rotation measure maps are a possible probe for the cosmic ray content of a given galaxy. Our calculations highlight the importance of cosmic rays in these measures, making them an essential ingredient of realistic models of the interstellar medium.
A linear stability analysis has been done to a magnetized disk under a linear gravity. We have reduced the linearized perturbation equations to a second-order differential equation which resembles the Schr{o}dinger equation with the potential of a harmonic oscillator. Depending on the signs of energy and potential terms, eigensolutions can be classified into ``continuum and ``discrete families. When magnetic field is ignored, the continuum family is identified as the convective mode, while the discrete family as acoustic-gravity waves. If the effective adiabatic index $gamma$ is less than unity, the former develops into the convective instability. When a magnetic field is included, the continuum and discrete families further branch into several solutions with different characters. The continuum family is divided into two modes: one is the original Parker mode, which is a slow MHD mode modulated by the gravity, and the other is a stable Alfven mode. The Parker modes can be either stable or unstable depending on $gamma$. When $gamma$ is smaller than a critical value $gamma_{cr}$, the Parker mode becomes unstable. The discrete family is divided into three modes: a stable fast MHD mode modulated by the gravity, a stable slow MHD mode modulated by the gravity, and an unstable mode which is also attributed to a slow MHD mode. The unstable discrete mode does not always exist. Even though the unstable discrete mode exists, the Parker mode dominates it if the Parker mode is unstable. However, if $gamma ge gamma_{cr}$, the discrete mode could be the only unstable one. When $gamma$ is equal $gamma_{cr}$, the minimum growth time of the unstable discrete mode is $1.3 times 10^8$ years with a corresponding length scale of 2.4 kpc. It is suggestive that the corrugatory features seen in the Galaxy and external galaxies are related to the unstable discrete mode.
Stellar winds are an integral part of the underlying dynamo, the motor of stellar activity. The wind controls the stars angular momentum loss, which depends on the magnetic field geometry which varies significantly in time and latitude. Here we study basic properties of a self-consistent model that includes simple representations of both the global stellar dynamo in a spherical shell and the exterior in which the wind accelerates and becomes supersonic. We numerically solve an axisymmetric mean-field model for the induction, momentum, and continuity equations using an isothermal equation of state. The model allows for the simultaneous generation of a mean magnetic field and the development of a Parker wind. The resulting flow is transonic at the critical point, which we arrange to be between the inner and outer radii of the model. The boundary conditions are assumed to be such that the magnetic field is antisymmetric about the equator, i.e., dipolar. At the solar rotation rate, the dynamo is oscillatory and of $alpha^2$ type. In most of the domain, the magnetic field corresponds to that of a split monopole. The magnetic energy flux is largest between the stellar surface and the critical point. The angular momentum flux is highly variable in time and can reach negative values, especially at midlatitudes. At rapid rotation of up to 50 times the solar value, most of the magnetic field is lost along the axis within the inner tangential cylinder of the model. The model reveals unexpected features that are not generally anticipated from models that are designed to reproduce the solar wind: highly variable angular momentum fluxes even from just an $alpha^2$ dynamo in the star. A major caveat of our isothermal models with a magnetic field produced by a dynamo is the difficulty to reach small enough plasma betas without the dynamo itself becoming unrealistically strong inside the star.
We study the role of ambipolar diffusion (AD) on the non-linear evolution of the MRI in protoplanetary disks using the strong coupling limit, which applies when the electron recombination time is much shorter than the orbital time. The effect of AD in this limit is characterized by the dimensionless number Am, the frequency of which neutral particles collide with ions normalized to the orbital frequency. We perform three-dimensional unstratified shearing-box simulations of the MRI over a wide range of Am as well as different magnetic field strengths and geometries. The saturation level of the MRI turbulence depends on the magnetic geometry and increases with the net magnetic flux. There is an upper limit to the net flux for sustained turbulence, corresponding to the requirement that the most unstable vertical wavelength be less than the disk scale height. Correspondingly, at a given Am, there exists a maximum value of the turbulent stress alpha_max. For Am<1, the largest stress is associated with a field geometry that has both net vertical and toroidal flux. In this case, we confirm the results of linear analyses that show the fastest growing mode has a non-zero radial wave number with growth rate exceeding the pure vertical field case. We find there is a very tight correlation between the turbulent stress (alpha) and the plasma beta=P_gas/P_mag~1/(2alpha) at the saturated state of the MRI turbulence regardless of field geometry, and alpha_max rapidly decreases with decreasing Am. In particular, we quote alpha_max~0.007 for Am=1 and alpha_max~0.0006 for Am=0.1.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
