Jérôme Guilet

Transport of magnetic flux and the vertical structure of accretion discs – I. Uniform diffusion coefficients

Standard models of accretion discs study the transport of mass on a viscous time-scale but do not consider the transport of magnetic flux. The evolution of a large-scale poloidal magnetic field is, however, an important problem because of its role in the launching of jets and winds and in determining the intensity of turbulence. As a consequence, the transport of poloidal magnetic flux should be considered on an equal basis to the transport of mass. In this paper, we develop a formalism to study such a transport of mass and magnetic flux in a thin accretion disc. The governing equations are derived by performing an asymptotic expansion in the limit of a thin disc, in the regime where the magnetic field is dominated by its vertical component. Turbulent viscosity and resistivity are included, with an arbitrary vertical profile that can be adjusted to mimic the vertical structure of the turbulence. At a given radius and time, the rates of transport of mass and magnetic flux are determined by a one-dimensional problem in the vertical direction, in which the radial gradients of various quantities appear as source terms. We solve this problem to obtain the transport rates and the vertical structure of the disc. This paper is then restricted to the idealized case of uniform diffusion coefficients, while a companion paper will study more realistic vertical profiles of these coefficients. We show the advection of weak magnetic fields to be significantly faster than the advection of mass, contrary to what a crude vertical averaging might suggest. This results from the larger radial velocities away from the mid-plane, which barely affect the mass accretion owing to the low density in these regions but do affect the advection of magnetic flux. Possible consequences of this larger accretion velocity include a potentially interesting time dependence with the magnetic flux distribution evolving faster than the mass distribution. If the disc is not too thin, this fast advection may also partially solve the long-standing problem of too efficient diffusion of an inclined magnetic field.

On the linear growth mechanism driving the standing accretion shock instability

During stellar core collapse, which eventually leads to a supernovae explosion, the stalled shock is unstable due to the standing accretion shock instability (SASI). This instability induces large-scale non spherical oscillations of the shock, which have crucial consequences on the dynamics and the geometry of the explosion. While the existence of this instability has been firmly established, its physical origin remains somewhat uncertain. Two mechanisms have indeed been proposed to explain its linear growth. The first is an advective-acoustic cycle, where the instability results from the interplay between advected perturbations (entropy and vorticity) and an acoustic wave. The second mechanism is purely acoustic and assumes that the shock is able to amplify trapped acoustic waves. Several arguments favouring the advective-acoustic cycle have already been proposed, however none was entirely conclusive for realistic flow parameters. In this article we give two new arguments which unambiguously show that the instability is not purely acoustic, and should be attributed to the advective-acoustic cycle. First, we extract a radial propagation timescale by comparing the frequencies of several unstable harmonics that differ only by their radial structure. The extracted time matches the advective-acoustic time but strongly disagrees with a purely acoustic interpretation. Second, we present a method to compute purely acoustic modes, by artificially removing advected perturbations below the shock. All these purely acoustic modes are found to be stable, showing that the advected wave is essential to the instability mechanism.

THE SATURATION OF SASI BY PARASITIC INSTABILITIES

The standing accretion shock instability (SASI) is commonly believed to be responsible for large amplitude dipolar oscillations of the stalled shock during core collapse, potentially leading to an asymmetric supernovae explosion. The degree of asymmetry depends on the amplitude of SASI, but the nonlinear saturation mechanism has never been elucidated. We investigate the role of parasitic instabilities as a possible cause of nonlinear SASI saturation. As the shock oscillations create both vorticity and entropy gradients, we show that both Kelvin-Helmholtz and Rayleigh-Taylor types of instabilities are able to grow on a SASI mode if its amplitude is large enough. We obtain simple estimates of their growth rates, taking into account the effects of advection and entropy stratification. In the context of the advective-acoustic cycle, we use numerical simulations to demonstrate how the acoustic feedback can be decreased if a parasitic instability distorts the advected structure. The amplitude of the shock deformation is estimated analytically in this scenario. When applied to the set up of Fernandez & Thompson, this saturation mechanism is able to explain the dramatic decrease of the SASI power when both the nuclear dissociation energy and the cooling rate are varied. Our results open new perspectives for anticipating the effect, on the SASI amplitude, of the physical ingredients involved in the modeling of the collapsing star.

Transport of magnetic flux and the vertical structure of accretion discs – II. Vertical profile of the diffusion coefficients

We investigate the radial transport of magnetic flux in a thin accretion disc, the turbulence being modelled by effective diffusion coefficients (viscosity and resistivity). Both turbulent diffusion and advection by the accretion flow contribute to flux transport, and they are likely to act in opposition. We study the consequences of the vertical variation of the diffusion coefficients, due to a varying strength of the turbulence. For this purpose, we consider three different vertical profiles of these coefficients. The first one is aimed at mimicking the turbulent stress profile observed in numerical simulations of MHD turbulence in stratified discs. This enables us to confirm the robustness of the main result of Paper I obtained for uniform diffusion coefficients that, for weak magnetic fields, the contribution of the accretion flow to the transport velocity of magnetic flux is much larger than the transport velocity of mass. We then consider the presence of a dead zone around the equatorial plane, where the physical resistivity is high while the turbulent viscosity is low. We find that it amplifies the previous effect: weak magnetic fields can be advected orders of magnitude faster than mass, for dead zones with a large vertical extension. The ratio of advection to diffusion, determining the maximum inclination of the field at the surface of the disc, is however not much affected. Finally, we study the effect of a non-turbulent layer at the surface of the disc, which has been suggested as a way to reduce the diffusion of the magnetic flux. We find that the reduction of the diffusion requires the conducting layer to extend below the height at which the magnetic pressure equals the thermal pressure. As a consequence, if the absence of turbulence is caused by the large-scale magnetic field, the highly conducting layer is inefficient at reducing the diffusion.

Shallow water analogue of the standing accretion shock instability: experimental demonstration and a two-dimensional model.

Despite the sphericity of the collapsing stellar core, the birth conditions of neutron stars can be highly nonspherical due to a hydrodynamical instability of the shocked accretion flow. Here we report the first laboratory experiment of a shallow water analogue, based on the physics of hydraulic jumps. Both the experiment and its shallow water modeling demonstrate a robust linear instability and nonlinear properties of symmetry breaking, in a system which is one million times smaller and about one hundred times slower than its astrophysical analogue.

The Explosion Mechanism of Core-Collapse Supernovae: Progress in Supernova Theory and Experiments

The explosion of core-collapse supernova depends on a sequence of events taking place in less than a second in a region of a few hundred kilometers at the center of a supergiant star, after the stellar core approaches the Chandrasekhar mass and collapses into a proto-neutron star, and before a shock wave is launched across the stellar envelope. Theoretical efforts to understand stellar death focus on the mechanism which transforms the collapse into an explosion. Progress in understanding this mechanism is reviewed with particular attention to its asymmetric character. We highlight a series of successful studies connecting observations of supernova remnants and pulsars properties to the theory of core-collapse using numerical simulations. The encouraging results from first principles models in axisymmetric simulations is tempered by new puzzles in 3D. The diversity of explosion paths and the dependence on the pre-collapse stellar structure is stressed, as well as the need to gain a better understanding of hydrodynamical and MHD instabilities such as SASI and neutrino-driven convection. The shallow water analogy of shock dynamics is presented as a comparative system where buoyancy effects are absent. This dynamical system can be studied numerically and also experimentally with a water fountain. The potential of this complementary research tool for supernova theory is analyzed. We also review its potential for public outreach in science museums.

Global evolution of the magnetic field in a thin disc and its consequences for protoplanetary systems

The strength and structure of the large-scale magnetic field in protoplanetary discs are still unknown, although they could have important consequences for the dynamics and evolution of the disc. Using a mean-field approach in which we model the effects of turbulence through enhanced diffusion coefficients, we study the time-evolution of the large-scale poloidal magnetic field in a global model of a thin accretion disc, with particular attention to protoplanetary discs. With the transport coefficients usually assumed, the magnetic field strength does not significantly increase radially inwards, leading to a relatively weak magnetic field in the inner part of the disc. We show that with more realistic transport coefficients that take into account the vertical structure of the disc and the back-reaction of the magnetic field on the flow as obtained by Guilet & Ogilvie (2012), the magnetic field can significantly increase radially inwards.

Numerical simulations of the magnetorotational instability in protoneutron stars – I. Influence of buoyancy

The magneto-rotational instability (MRI) is considered to be a promising mechanism to amplify the magnetic field in fast rotating protoneutron stars. In contrast to accretion disks, radial buoyancy driven by entropy and lepton fraction gradients is expected to have a dynamical role as important as rotation and shear. We investigate the poorly known impact of buoyancy on the non-linear phase of the MRI, by means of three dimensional numerical simulations of a local model in the equatorial plane of a protoneutron star. The use of the Boussinesq approximation allows us to utilise a shearing box model with clean shearing periodic boundary conditions, while taking into account the buoyancy driven by radial entropy and composition gradients. We find significantly stronger turbulence and magnetic fields in buoyantly unstable flows. On the other hand, buoyancy has only a limited impact on the strength of turbulence and magnetic field amplification for buoyantly stable flows in the presence of a realistic thermal diffusion. The properties of the turbulence are, however, significantly affected in the latter case. In particular, the toroidal components of the magnetic field and of the velocity become even more dominant with respect to the poloidal ones. Furthermore, we observed in the regime of stable buoyancy the formation of long lived coherent structures such as channel flows and zonal flows. Overall, our results support the ability of the MRI to amplify the magnetic field significantly even in stably stratified regions of protoneutron stars.

Angular momentum redistribution by SASI spiral modes and consequences for neutron star spins

In the collapsing core of massive stars, the standing accretion shock instability (SASI) can drive spiral modes that efficiently redistribute angular momentum. This process can impart a spin to the forming neutron star even when the progenitor star is nonrotating. Here we develop the first analytical description of the angular momentum redistribution driven by a spiral mode of the SASI. Our analysis, valid in the limit of small mode amplitude, shows that the angular momentum separation is driven by the Reynolds stress generated by the spiral mode. The resulting solutions compare favorably with previous three-dimensional hydrodynamic simulations of the SASI in the linear and weakly non-linear phases. Reasonable agreement is also found when extrapolating the solutions into the fully non-linear phase. A Reynolds-decomposition of the flow is performed in the saturated state of these simulations, showing that outward angular momentum transport by the Reynolds stress and the fluctuating component of the mass flux balance inward transport by advection. We derive an approximate analytic expression for the maximum angular momentum deposited in the neutron star as a function of the mass accretion rate, shock radius, shock compression ratio, and amplitude of the spiral mode at the time of explosion. Implications for the birth spin periods of neutron stars are discussed.

New insights on the spin-up of a neutron star during core collapse

The spin of a neutron star at birth may be impacted by the asymmetric character of the explosion of its massive progenitor. During the first second after bounce, the spiral mode of the Standing Accretion Shock Instability (SASI) is able to redistribute angular momentum and spin-up a neutron star born from a non-rotating progenitor. Our aim is to assess the robustness of this process. We perform 2D numerical simulations of a simplified setup in cylindrical geometry to investigate the timescale over which the dynamics is dominated by a spiral or a sloshing mode. We observe that the spiral mode prevails only if the ratio of the initial shock radius to the neutron star radius exceeds a critical value. In that regime, both the degree of asymmetry and the average expansion of the shock induced by the spiral mode increase monotonously with this ratio, exceeding the values obtained when a sloshing mode is artificially imposed. With a timescale of 2-3 SASI oscillations, the dynamics of SASI takes place fast enough to affect the spin of the neutron star before the explosion. The spin periods deduced from the simulations are compared favorably to analytical estimates evaluated from the measured saturation amplitude of the SASI wave. Despite the simplicity of our setup, numerical simulations revealed unexpected stochastic variations, including a reversal of the direction of rotation of the shock. Our results show that the spin up of neutron stars by SASI spiral modes is a viable mechanism even though it is not systematic.