Magnetars are highly magnetized neutron stars whose magnetic dipole ranges from $10^{14}$ to $10^{15}$ G. The MRI is considered to be a promising mechanism to amplify the magnetic field in fast-rotating protoneutron stars and form magnetars. This scenario is supported by many local studies showing that magnetic fields could be amplified by the MRI on small scales. However, the efficiency of the MRI at generating a dipole field is still unknown. To answer this question, we study the MRI dynamo in an idealized global model of a fast rotating protoneutron star with differential rotation. We perform 3D incompressible MHD simulations in spherical geometry with explicit diffusivities where the differential rotation is forced at the outer boundary. We vary the initial magnetic field and investigated different magnetic boundary conditions. These simulations were compared to local shearing box simulations. We obtain a self-sustained turbulent MRI-driven dynamo, whose saturated state is independent of the initial magnetic field. The MRI generates a strong turbulent magnetic field of $B \geq 2\times 10^{15}$ G and a non-dominant magnetic dipole, which represents systematically about $5\%$ of the averaged magnetic field strength. Interestingly, this dipole is tilted towards the equatorial plane. We find that local shearing box models can reproduce fairly well several characteristics of global MRI turbulence such as the kinetic and magnetic spectra. The turbulence is nonetheless more vigorous in the local models than in the global ones. Overall, our results support the ability of the MRI to form magnetar-like large-scale magnetic fields. They furthermore predict the presence of a stronger small-scale magnetic field. The resulting magnetic field could be important to power outstanding stellar explosions, such as superluminous supernovae and GRBs.

The release of spin-down energy by a magnetar is a promising scenario to power several classes of extreme explosive transients. However, it lacks a firm basis because magnetar formation still represents a theoretical challenge. Using the first three-dimensional simulations of a convective dynamo based on a protoneutron star interior model, we demonstrate that the required dipolar magnetic field can be consistently generated for sufficiently fast rotation rates. The dynamo instability saturates in the magnetostrophic regime with the magnetic energy exceeding the kinetic energy by a factor of up to 10. Our results are compatible with the observational constraints on galactic magnetar field strength and provide strong theoretical support for millisecond protomagnetar models of gamma-ray burst and superluminous supernova central engines.

The magnetic field is believed to play an important role in at least some core-collapse supernovae if its magnitude reaches $10^{15}\,\rm{G}$, which is a typical value for a magnetar. In the presence of fast rotation, such a strong magnetic field can drive powerful jet-like explosions if it has the large-scale coherence of a dipole. The topology of the magnetic field is, however, probably much more complex with strong multipolar and small-scale components and the consequences for the explosion are so far unclear. We investigate the effects of the magnetic field topology on the dynamics of core-collapse supernovae and the properties of forming proto-neutron star (PNS) by comparing pre-collapse fields of different multipolar orders and radial profiles. Using axisymmetric special relativistic MHD simulations and a two-moment neutrino transport, we find that higher multipolar magnetic configurations lead to generally less energetic explosions, slower expanding shocks and less collimated outflows. Models with a low order multipolar configuration tend to produce more oblate PNS, surrounded in some cases by a rotationally supported toroidal structure of neutron-rich material. Moreover, magnetic fields which are distributed on smaller angular scales produce more massive and faster rotating central PNS, suggesting that higher-order multipolar configurations tend to decrease the efficiency of the magnetorotational launching mechanism. Even if our dipolar models systematically display a far more efficient extraction of the rotational energy of the PNS, fields distributed on smaller angular scales are still capable of powering magnetorotational explosions and shape the evolution of the central compact object.

A toy model of the post-shock region of core-collapse supernovae is used to study the non-linear development of turbulent motions driven by convection in the presence of advection. Our numerical simulations indicate that buoyant perturbations of density are able to trigger self-sustained convection only when the instability is not linearly stabilized by advection. Large amplitude perturbations produced by strong shock oscillations or combustion inhomogeneities before the collapse of the progenitor are efficiently shredded through phase mixing and generate a turbulent cascade. Our model enables us to investigate several physical arguments that had been proposed to explain the impact of the dimensionality on the onset of explosions in global simulations of core-collapse supernovae. Three-dimensional (3D) simulations are found to lead to higher entropy values than two-dimensional (2D) ones. We attribute this to greater turbulent mixing and dissipation of the kinetic energy into heat in 3D. Our results show that the increase of entropy is enhanced with finer numerical resolution and larger perturbation amplitude.

Geometrically thick tori with constant specific angular momentum have been widely used in the last decades to construct numerical models of accretion flows onto black holes. Such discs are prone to a global non-axisymmetric hydrodynamic instability, known as Papaloizou-Pringle instability (PPI), which can redistribute angular momentum and also lead to an emission of gravitational waves. It is, however, not clear yet how the development of the PPI is affected by the presence of a magnetic field and by the concurrent development of the magnetorotational instability (MRI). We present a numerical analysis using three-dimensional GRMHD simulations of the interplay between the PPI and the MRI considering, for the first time, an analytical magnetized equilibrium solution as initial condition. In the purely hydrodynamic case, the PPI selects as expected the large-scale $m=1$ azimuthal mode as the fastest growing and non-linearly dominant mode. However, when the torus is threaded by a weak toroidal magnetic field, the development of the MRI leads to the suppression of large-scale modes and redistributes power across smaller scales. If the system starts with a significantly excited $m=1$ mode, the PPI can be dominant in a transient phase, before being ultimately quenched by the MRI. Such dynamics may well be important in compact star mergers and tidal disruption events.

The merger of two neutron stars may give birth to a long-lived hypermassive neutron star. If it harbours a strong magnetic field of magnetar strength, its spin-down could explain several features of short gamma-ray burst afterglows. The magnetorotational instability (MRI) has been proposed as a mechanism to amplify the magnetic field to the required strength. Previous studies have, however, neglected neutrinos, which may have an important impact on the MRI by inducing a viscosity and drag. We investigate the impact of these neutrino effects on the linear growth of the MRI by applying a local stability analysis to snapshots of a neutron star merger simulation. We find that neutrinos have a significant impact inside the hypermassive neutron star, but have at most a marginal effect in the torus surrounding it. Inside the hypermassive neutron star, the MRI grows in different regimes depending on the radius and on the initial magnetic field strength. For magnetic fields weaker than $10^{13}-10^{14}\,{\rm G}$, the growth rate of the MRI is significantly reduced due to the presence of neutrinos. We conclude that neutrinos should be taken into account when studying the growth of the MRI from realistic initial magnetic fields. Current numerical simulations, which neglect neutrino viscosity, are only consistent, i.e. in the adopted ideal regime, if they start from artificially strong initial magnetic fields above $\sim10^{14}\,{\rm G}$. One should be careful when extrapolating these results to lower initial magnetic fields, where the MRI growth is strongly affected by neutrino viscosity or drag.

We discuss the relevance of the magnetorotational instability (MRI) in core-collapse supernovae (CCSNe). Our recent numerical studies show that in CCSNe, the MRI is terminated by parasitic instabilities of the Kelvin-Helmholtz type. To determine whether the MRI can amplify initially weak magnetic fields to dynamically relevant strengths in CCSNe, we performed three-dimensional simulations of a region close to the surface of a differentially rotating proto-neutron star in non-ideal magnetohydrodynamics with two different numerical codes. We find that under the conditions prevailing in proto-neutron stars, the MRI can amplify the magnetic field by (only) one order of magnitude. This severely limits the role of MRI channel modes as an agent amplifying the magnetic field in proto-neutron stars starting from small seed fields.

Pulsars may either be spun up or down by hydrodynamic instabilities during the supernova explosion of massive stars. Besides rapidly-rotating cases related to bipolar explosions, stellar rotation may affect the explosion of massive stars in the more common situations where the centrifugal force is minor. Using 2D simulations of a simplified setup in cylindrical geometry, we examine the impact of rotation on the Standing Accretion Shock Instability (SASI) and the corotation instability, also known as low-T/|W|. The influence of rotation on the saturation amplitude of these instabilities depends on the specific angular momentum in the accretion flow and the ratio of the shock to the neutron star radii. The spiral mode of SASI becomes more vigorous with faster rotation only if this ratio is large enough. A corotation instability develops at large rotation rates and impacts the dynamics more dramatically, leading to a strong one-armed spiral wave. Non-axisymmetric instabilities are able to redistribute angular momentum radially and affect the pulsar spin at birth. A systematic study of the relationship between the core rotation period of the progenitor and the initial pulsar spin is performed. Stellar rotation rates for which pulsars are spun up or down by SASI are estimated. Rapidly spinning progenitors are modestly spun down by spiral modes, less than $\sim 30\%$, when a corotation instability develops. Given the observational constraints on pulsar spin periods at birth, this suggests that rapid rotation might not play a significant hydrodynamic role in most core-collapse supernovae.

Extremely strong magnetic fields of the order of $10^{15}\,{\rm G}$ are required to explain the properties of magnetars, the most magnetic neutron stars. Such a strong magnetic field is expected to play an important role for the dynamics of core-collapse supernovae, and in the presence of rapid rotation may power superluminous supernovae and hypernovae associated to long gamma-ray bursts. The origin of these strong magnetic fields remains, however, obscure and most likely requires an amplification over many orders of magnitude in the protoneutron star. One of the most promising agents is the magnetorotational instability (MRI), which can in principle amplify exponentially fast a weak initial magnetic field to a dynamically relevant strength. We describe our current understanding of the MRI in protoneutron stars and show recent results on its dependence on physical conditions specific to protoneutron stars such as neutrino radiation, strong buoyancy effects and large magnetic Prandtl number.

To take advantage of the astrophysical potential of Gamma-Ray Bursts (GRBs), Chinese and French astrophysicists have engaged the SVOM mission (Space-based multi-band astronomical Variable Objects Monitor). Major advances in GRB studies resulting from the synergy between space and ground observations, the SVOM mission implements space and ground instrumentation. The scientific objectives of the mission put a special emphasis on two categories of GRBs: very distant GRBs at z$>$5 which constitute exceptional cosmological probes, and faint/soft nearby GRBs which allow probing the nature of the progenitors and the physics at work in the explosion. These goals have a major impact on the design of the mission: the on-board hard X-ray imager is sensitive down to 4 keV and computes on line image and rate triggers, and the follow-up telescopes on the ground are sensitive in the NIR. At the beginning of the next decade, SVOM will be the main provider of GRB positions and spectral parameters on very short time scale. The SVOM instruments will operate simultaneously with a wide range of powerful astronomical devices. This rare instrumental conjunction, combined with the relevance of the scientific topics connected with GRB studies, warrants a remarkable scientific return for SVOM. In addition, the SVOM instrumentation, primarily designed for GRB studies, composes a unique multi-wavelength observatory with rapid slew capability that will find multiple applications for the whole astronomy community beyond the specific objectives linked to GRBs. This report lists the scientific themes that will benefit from observations made with SVOM, whether they are specific GRB topics, or more generally all the issues that can take advantage of the multi-wavelength capabilities of SVOM.

Recent numerical simulations suggest that the model by Zahn (1992, A&A, 265, 115) for the turbulent mixing of chemical elements due to differential rotation in stellar radiative zones is valid. We investigate the robustness of this result with respect to the numerical configuration and Reynolds number of the flow. We compare results from simulations performed with two different numerical codes, including one that uses the shearing-box formalism. We also extensively study the dependence of the turbulent diffusion coefficient on the turbulent Reynolds number. The two numerical codes used in this study give consistent results. The turbulent diffusion coefficient is independent of the size of the numerical domain if at least three large turbulent structures fit in the box. Generally, the turbulent diffusion coefficient depends on the turbulent Reynolds number. However, our simulations suggest that an asymptotic regime is obtained when the turbulent Reynolds number is larger than $10^3$. Shear mixing in the regime of small P\'eclet numbers can be investigated numerically both with shearing-box simulations and simulations using explicit forcing. Our results suggest that Zahn's model is valid at large turbulent Reynolds numbers.

Whether the magnetorotational instability (MRI) can amplify initially weak magnetic fields to dynamically relevant strengths in core collapse supernovae is still a matter of active scientific debate. Recent numerical studies have shown that the first phase of MRI growth dominated by channel flows is terminated by parasitic instabilities of the Kelvin-Helmholtz type that disrupt MRI channel flows and quench further magnetic field growth. However, it remains to be prop- erly assessed by what factor the initial magnetic field can be amplified and how it depends on the initial field strength and the amplitude of the perturbations. Different termination criteria leading to different estimates of the amplification factor were proposed within the parasitic model. To determine the amplification factor and test which criterion is a better predictor of the MRI termination, we perform three-dimensional shearing-disc and shearing-box simula- tions of a region close to the surface of a differentially rotating proto-neutron star in non-ideal MHD with two different numerical codes. We find that independently of the initial magnetic field strength, the MRI channel modes can amplify the magnetic field by, at most, a factor of 100. Under the conditions found in proto-neutron stars a more realistic value for the mag- netic field amplification is of the order of 10. This severely limits the role of the MRI channel modes as an agent amplifying the magnetic field in proto-neutron stars starting from small seed fields. A further amplification should therefore rely on other physical processes, such as for example an MRI-driven turbulent dynamo.

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.

We determine how MRI-turbulent stresses depend on gas pressure via a suite of unstratified shearing box simulations. Earlier numerical work reported only a very weak dependence at best, results that call into question the canonical alpha-disk model and the thermal stability results that follow from it. Our simulations, in contrast, exhibit a stronger relationship, and show that previous work was box-size limited: turbulent `eddies' were artificially restricted by the numerical domain rather than by the scale height. Zero-net-flux runs without physical diffusion coefficients yield a stress proportional to $P^{0.5}$, where P is pressure. The stresses are also proportional to the grid length and hence remain numerically unconverged. The same runs with physical diffusivities, however, give a result closer to an alpha-disk: the stress is proportional to $P^{0.9}$. Net-flux simulations without explicit diffusion exhibit stresses proportional to $P^{0.5}$, but stronger imposed fields weaken this correlation. In summary, compressibility is important for the saturation of the MRI, but the exact stress-pressure relationship is difficult to ascertain in local simulations because of numerical convergence issues and the influence of any imposed flux. As a consequence, the interpretation of thermal stability behaviour in local simulations is a problematic enterprise.

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.

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.

The magnetorotational instability (MRI) is a promising mechanism to amplify the magnetic field in fast rotating protoneutron stars. The diffusion of neutrinos trapped in the PNS induces a transport of momentum, which can be modelled as a viscosity on length-scales longer than the neutrino mean free path. This neutrino-viscosity can slow down the growth of MRI modes to such an extent that a minimum initial magnetic field strength of $\gtrsim 10^{12}\, {\rm G}$ is needed for the MRI to grow on a sufficiently short time-scale to potentially affect the explosion. It is uncertain whether the magnetic field of fast rotating progenitor cores is strong enough to yield such an initial magnetic field in PNS. At MRI wavelengths shorter than the neutrino mean free path, on the other hand, neutrino radiation does not act as a viscosity but rather induces a drag on the velocity with a damping rate independent of the wavelength. We perform a linear analysis of the MRI in this regime, and apply our analytical results to the PNS structure from a one-dimensional numerical simulation. We show that in the outer layers of the PNS, the MRI can grow from weak magnetic fields at wavelengths shorter than the neutrino mean free path, while deeper in the PNS MRI growth takes place in the viscous regime and requires a minimum magnetic field strength.

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 non-rotating. 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.

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. The magnetic-field profile adjusts to reach an equilibrium value of the plasma $\beta$ parameter (the ratio of midplane thermal pressure to magnetic pressure) in the inner part of the disc. This value of $\beta$ depends strongly on the aspect ratio of the disc and on the turbulent magnetic Prandtl number, and lies in the range $10^4-10^7$ for protoplanetary discs. Such a magnetic field is expected to affect significantly the dynamics of protoplanetary discs by increasing the strength of MHD turbulence and launching an outflow. We discuss the implications of our results for the evolution of protoplanetary discs and for the formation of powerful jets as observed in T-Tauri star systems.

The great diversity of extrasolar planetary systems has challenged our understanding of how planets form, and how their orbits evolve as they form. Among the various processes that may account for this diversity, the gravitational interaction between planets and their parent protoplanetary disc plays a prominent role in shaping young planetary systems. Planet-disc forces are large, and the characteristic times for the evolution of planets orbital elements are much shorter than the lifetime of protoplanetary discs. The determination of such forces is challenging, because it involves many physical mechanisms and it requires a detailed knowledge of the disc structure. Yet, the intense research of the past few years, with the exploration of many new avenues, represents a very significant improvement on the state of the discipline. This chapter reviews current understanding of planet-disc interactions, and highlights their role in setting the properties and architecture of observed planetary systems.

[Abridged] The migration of low mass planets has been studied in hydrodynamical disc models for more than three decades, but the impact of a magnetic field in the protoplanetary disc is less known. When the disc's magnetic field is strong enough to prevent horseshoe motion, the corotation torque is replaced by a torque arising from magnetic resonances. For weak enough magnetic fields, horseshoe motion and a corotation torque exist, and recent turbulent MHD simulations have reported the existence of a new component of the corotation torque in the presence of a mean toroidal field. The aim of this paper is to investigate the physical origin and the properties of this new corotation torque. We performed MHD simulations of a low mass planet embedded in a 2D laminar disc threaded by a weak toroidal magnetic field, with the effects of turbulence modelled by a viscosity and a resistivity. We confirm that the interaction between the magnetic field and the horseshoe motion results in an additional corotation torque on the planet, which we dub the MHD torque excess. It is caused by the accumulation of the magnetic field along the downstream separatrices of the planet's horseshoe region, which gives rise to an azimuthally asymmetric underdense region at that location. The properties of the MHD torque excess are characterised by varying the slope of the density, temperature and magnetic field profiles, as well as the diffusion coefficients and the strength of the magnetic field. The sign of the MHD torque excess depends on the density and temperature gradients only, and is positive for profiles expected in protoplanetary discs. Its magnitude is in turn mainly determined by the strength of the magnetic field and the turbulent resistivity. The MHD torque excess can be strong enough to reverse migration, even when the magnetic pressure is less than one percent of the thermal pressure.

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.

The evolution of a large-scale poloidal magnetic field in accretion discs is an important problem because of its role in the launching of jets and winds and in determining the intensity of turbulence. In this paper, we develop a formalism to calculate the transport magnetic flux in a thin accretion disc, thus determining its evolution on a viscous/resistive timescale. 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. The present paper is then restricted to the idealised 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.

Despite the sphericity of the collapsing stellar core, the birth conditions of neutron stars can be highly non spherical 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 hundred times slower than its astrophysical analogue.

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.

We investigate the dynamics of an Alfven surface (where the Alfven speed equals the advection velocity) in the context of core collapse supernovae during the phase of accretion on the proto-neutron star. Such a surface should exist even for weak magnetic fields because the advection velocity decreases to zero at the center of the collapsing core. In this decelerated flow, Alfven waves created by the standing accretion shock instability (SASI) or convection accumulate and amplify while approaching the Alfven surface. We study this amplification using one dimensional MHD simulations with explicit physical dissipation. In the linear regime, the amplification continues until the Alfven wavelength becomes as small as the dissipative scale. A pressure feedback that increases the pressure in the upstream flow is created via a non linear coupling. We derive analytic formulae for the maximum amplification and the non linear coupling and check them with numerical simulations to a very good accuracy. We also characterize the non linear saturation of this amplification when compression effects become important, leading to either a change of the velocity gradient, or a steepening of the Alfven wave. Applying these results to core collapse supernovae shows that the amplification can be fast enough to affect the dynamics, if the magnetic field is strong enough for the Alfven surface to lie in the region of strong velocity gradient just above the neutrinosphere. This requires the presence of a strong magnetic field in the progenitor star, which would correspond to the formation of a magnetar under the assumption of magnetic flux conservation. An extrapolation of our analytic formula (taking into account the nonlinear saturation) suggests that the Alfven wave could reach an amplitude of B ~ 10^15 G, and that the pressure feedback could significantly contribute to the pressure below the shock.

The effect of a magnetic field on the linear phase of the advective-acoustic instability is investigated, as a first step toward a magnetohydrodynamic (MHD) theory of the stationary accretion shock instability taking place during stellar core collapse. We study a toy model where the flow behind a planar stationary accretion shock is adiabatically decelerated by an external potential. Two magnetic field geometries are considered: parallel or perpendicular to the shock. The entropy-vorticity wave, which is simply advected in the unmagnetized limit, separates into five different waves: the entropy perturbations are advected, while the vorticity can propagate along the field lines through two Alfven waves and two slow magnetosonic waves. The two cycles existing in the unmagnetized limit, advective-acoustic and purely acoustic, are replaced by up to six distinct MHD cycles. The phase differences among the cycles play an important role in determining the total cycle efficiency and hence the growth rate. Oscillations in the growth rate as a function of the magnetic field strength are due to this varying phase shift. A vertical magnetic field hardly affects the cycle efficiency in the regime of super-Alfvenic accretion that is considered. In contrast, we find that a horizontal magnetic field strongly increases the efficiencies of the vorticity cycles that bend the field lines, resulting in a significant increase of the growth rate if the different cycles are in phase. These magnetic effects are significant for large-scale modes if the Alfven velocity is a sizable fraction of the flow velocity.

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, which 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 (2009a), 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.