Magnetic activity in accretion disc boundary layers

We use three dimensional magnetohydrodynamic simulations to study the structure of the boundary layer between an accretion disc and a non-rotating, unmagnetized star. Under the assumption that cooling is efficient, we obtain a narrow but highly variable transition region in which the radial velocity is only a small fraction of the sound speed. A large fraction of the energy dissipation occurs in high density gas adjacent to the hydrostatic stellar envelope, and may therefore be reprocessed and largely hidden from view of the observer. As suggested by Pringle (1989), the magnetic field energy in the boundary layer is strongly amplified by shear, and exceeds that in the disc by an order of magnitude. These fields may play a role in generating the magnetic activity, X-ray emission, and outflows in disc systems where the accretion rate is high enough to overwhelm the stellar magnetosphere.Comment: MNRAS, in press. Movies from the simulations are available at http://star-www.st-and.ac.uk/~pja3/movies.htm

Magnetorotational instability in protoplanetary discs: The effect of dust grains

We investigate the linear growth and vertical structure of the MRI in protoplanetary discs when dust grains are well mixed with the gas over the entire disc thickness. All the grains have the same radius (a = 0.1, 1 or 3 micron) and constitute 1 % of the total mass of the gas. Solutions are obtained at R = 5 and 10 AU for a minimum-mass solar nebula model and different choices of the initially vertical magnetic field strength (B), configuration of the diffusivity tensor and grain sizes. We find that when no grains are present, or they are > 1 micron, the midplane remains magnetically coupled for B up to a few gauss at both radii. In contrast, when a population of small grains (a = 0.1 micron) is present, the disc is magnetically inactive for z/H < 2 and only B < 50 mG couple to the fluid. At 5 AU, Ohmic diffusion dominates for z/H < 1 when B < a few mG, irrespective of the properties of the grain population. Conversely, at 10 AU this diffusion term is unimportant in all the scenarios studied here. For z/H > 5, ambipolar diffusion is severe and prevents the field from coupling to the gas for all B. Hall diffusion is dominant for a wide range of field strengths at both radii when dust grains are present. The growth rate, wavenumber and range of magnetic field strengths for which MRI-unstable modes exist are all drastically diminished when dust grains are present, particularly when they are small (a ~ 0.1 micron). We conclude that in protoplanetary discs, the magnetic field is able to couple to the gas and shear over a wide range of fluid conditions even when small dust grains are well mixed with the gas. Despite the low magnetic coupling, MRI modes grow for an extended range of magnetic field strengths and Hall diffusion largely determines the properties of the perturbations in the inner regions of the disc (abridged).Comment: 17 pages, 11 figures. Submitted to MNRA

Oscillations of MHD shock waves on the surfaces of T Tauri stars

This work treats the matter deceleration in a magnetohydrodynamics radiative shock wave at the surface of a star. The problem is relevant to classical T Tauri stars where infalling matter is channeled along the star's magnetic field and stopped in the dense layers of photosphere. A significant new aspect of the present work is that the magnetic field has an arbitrary angle with respect to the normal to the star's surface. We consider the limit where the magnetic field at the surface of the star is not very strong in the sense that the inflow is super Alfv\'enic. In this limit the initial deceleration and heating of plasma (at the entrance to the cooling zone) occurs in a fast magnetohydrodynamic shock wave. To calculate the intensity of radiative losses we use "real" and "power-law" radiative functions. We determine the stability/instability of the radiative shock wave as a function of parameters of the incoming flow: velocity, strength of the magnetic field, and its inclination to the surface of the star. In a number of simulation runs with the "real" radiative function, we find a simple criterion for stability of the radiative shock wave. For a wide range of parameters, the periods of oscillation of the shock wave are of the order 0.02-0.2 sec.Comment: 13 pages, 11 figures, submitted to the MNRAS after revision, see animations at http://www.astro.cornell.edu/us-rus/shock.ht

Linear Polarization in Gamma-Ray Bursts: The Case for an Ordered Magnetic Field

Peer reviewe