The convection that takes place in the innermost shells of massive stars plays an important role in the formation of core-collapse supernova explosions. Upon encountering the supernova shock, additional turbulence is generated, amplifying the explosion. In this work, we study how the convective perturbations evolve during the stellar collapse. Our main aim is to establish their physical properties right before they reach the supernova shock. To this end, we solve the linearized hydrodynamics equations perturbed on a stationary background flow. The latter is approximated by the spherical transonic Bondi accretion, while the convective perturbations are modeled as a combination of entropy and vorticity waves. We follow their evolution from large radii, where convective shells are initially located, down to small radii, where they are expected to encounter the accretion shock above the proto-neutron star. Considering typical vorticity perturbations with a Mach number $\sim 0.1$ and entropy perturbations with magnitude $\sim 0.05 k_\mathrm{b}/\mathrm{baryon}$, we find that the advection of these perturbations down to the shock generates acoustic waves with a relative amplitude $\delta p/\gamma p \lesssim 10\%$, in agreement with published numerical simulations. The velocity perturbations consist of contributions from acoustic and vorticity waves with values reaching $\sim 10\%$ of the sound speed ahead of the shock. The perturbation amplitudes decrease with increasing $\ell$ and initial radii of the convective shells.

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.

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.

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.

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.

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 success of the neutrino mechanism of core-collapse supernovae relies on the supporting action of two hydrodynamic instabilities: neutrino-driven convection and the Standing Accretion Shock Instability (SASI). Depending on the structure of the stellar progenitor, each of these instabilities can dominate the evolution of the gain region prior to the onset of explosion, with implications for the ensuing asymmetries. Here we examine the flow dynamics in the neighborhood of explosion by means of parametric two-dimensional, time-dependent hydrodynamic simulations for which the linear stability properties are well understood. We find that systems for which the convection parameter is sub-critical (SASI-dominated) develop explosions once large-scale, high-entropy bubbles are able to survive for several SASI oscillation cycles. These long-lived structures are seeded by the SASI during shock expansions. Finite-amplitude initial perturbations do not alter this outcome qualitatively, though they can lead to significant differences in explosion times. Supercritical systems (convection-dominated) also explode by developing large-scale bubbles, though the formation of these structures is due to buoyant activity. Non-exploding systems achieve a quasi-steady state in which the time-averaged flow adjusts itself to be convectively sub-critical. We characterize the turbulent flow using a spherical Fourier-Bessel decomposition, identifying the relevant scalings and connecting temporal and spatial components. Finally, we verify the applicability of these principles on the general relativistic, radiation-hydrodynamic simulations of Mueller, Janka, & Heger (2012), and discuss implications for the three-dimensional case.

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.

Some general properties of the advective-acoustic instability are described and understood using a toy model which is simple enough to allow for analytical estimates of the eigenfrequencies. The essential ingredients of this model, in the unperturbed regime, are a stationary shock and a subsonic region of deceleration. For the sake of analytical simplicity, the 2D unperturbed flow is parallel and the deceleration is produced adiabatically by an external potential. The instability mechanism is determined unambiguously as the consequence of a cycle between advected and acoustic perturbations. The purely acoustic cycle, considered alone, is proven to be stable in this flow. Its contribution to the instability can be either constructive or destructive. A frequency cut-off is associated to the advection time through the region of deceleration. This cut-off frequency explains why the instability favours eigenmodes with a low frequency and a large horizontal wavelength. The relation between the instability occurring in this highly simplified toy model and the properties of SASI observed in the numerical simulations of stellar core-collapse is discussed. This simple set up is proposed as a benchmark test to evaluate the accuracy, in the linear regime, of numerical simulations involving this instability. We illustrate such benchmark simulations in a companion paper.

The physical processes involved in the advective-acoustic instability are investigated with 2D numerical simulations. Simple toy models, developped in a companion paper, are used to describe the coupling between acoustic and entropy/vorticity waves, produced either by a stationary shock or by the deceleration of the flow. Using two Eulerian codes based on different second order upwind schemes, we confirm the results of the perturbative analysis. The numerical convergence with respect to the computation mesh size is studied with 1D simulations. We demonstrate that the numerical accuracy of the quantities which depend on the physics of the shock is limited to a linear convergence. We argue that this property is likely to be true for most current numerical schemes dealing with SASI in the core-collapse problem, and could be solved by the use of advanced techniques for the numerical treatment of the shock. We propose a strategy to choose the mesh size for an accurate treatment of the advective-acoustic coupling in future numerical simulations.

A perturbative analysis is used to investigate the effect of rotation on the instability of a steady accretion shock (SASI) in a simple toy-model, in view of better understanding supernova explosions in which the collapsing core contains angular momentum. A cylindrical geometry is chosen for the sake of simplicity. Even when the centrifugal force is very small, rotation can have a strong effect on the non axisymmetric modes of SASI by increasing the growth rate of the spiral modes rotating in the same direction as the steady flow. Counter-rotating spiral modes are significantly damped, while axisymmetric modes are hardly affected by rotation. The growth rates of spiral modes have a nearly linear dependence on the specific angular momentum of the flow. The fundamental one-armed spiral mode (m=1) is favoured for small rotation rates, whereas stronger rotation rates favour the mode m=2. A WKB analysis of higher harmonics indicates that the efficiency of the advective-acoustic cycles associated to spiral modes is strongly affected by rotation in the same manner as low frequency modes, whereas the purely acoustic cycles are stable. These results suggest that the linear phase of SASI in rotating core-collapse supernovae naturally selects a spiral mode rotating in the same direction of the flow, as observed in the 3D numerical simulations of Blondin & Mezzacappa (2007). This emphasizes the need for a 3D approach of rotating core-collapse, before conclusions on the explosion mechanisms and pulsar kicks can be drawn.

By 2D hydrodynamic simulations including a detailed equation of state and neutrino transport, we investigate the interplay between different non-radial hydrodynamic instabilities that play a role during the postbounce accretion phase of collapsing stellar cores. The convective mode of instability, which is driven by negative entropy gradients caused by neutrino heating or by time variations of the shock strength, can be identified clearly by the development of typical Rayleigh-Taylor mushrooms. However, in cases where the gas in the postshock region is rapidly advected towards the gain radius, the growth of such a buoyancy instability can be suppressed. In such a situation the shocked flow nevertheless can develop non-radial asymmetry with an oscillatory growth of the amplitude. This phenomenon has been termed ``standing accretion shock instability'' (SASI). It is shown here that the SASI oscillations can trigger convective instability and like the latter they lead to an increase of the average shock radius and of the mass in the gain layer. Both hydrodynamic instabilities in combination stretch the advection time of matter through the neutrino-heating layer and thus enhance the neutrino energy deposition in support of the neutrino-driven explosion mechanism. A rapidly contracting and more compact nascent NS turns out to be favorable for explosions, because the accretion luminosity and neutrino heating are larger and the growth rate of the SASI is higher. Moreover, we show that the oscillation period of the SASI and a variety of other features in our simulations agree with estimates for the advective-acoustic cycle (AAC), in which perturbations are carried by the accretion flow from the shock to the neutron star and pressure waves close an amplifying global feedback loop. (abridged)

We analyze the linear stability of a stalled accretion shock in a perfect gas with a parametrized cooling function L ~ rho^{beta-alpha} P^alpha. The instability is dominated by the l=1 mode if the shock radius exceeds 2-3 times the accretor radius, depending on the parameters of the cooling function. The growth rate and oscillation period are comparable to those observed in the numerical simulations of Blondin & Mezzacappa (2006). The instability mechanism is analyzed by separately measuring the efficiencies of the purely acoustic cycle and the advective-acoustic cycle. These efficiencies are estimated directly from the eigenspectrum, and also through a WKB analysis in the high frequency limit. Both methods prove that the advective-acoustic cycle is unstable, and that the purely acoustic cycle is stable. Extrapolating these results to low frequency leads us to interpret the dominant mode as an advective-acoustic instability, different from the purely acoustic interpretation of Blondin & Mezzacappa (2006). A simplified characterization of the instability is proposed, based on an advective-acoustic cycle between the shock and the radius r_nabla where the velocity gradients of the stationary flow are strongest. The importance of the coupling region in this mechanism calls for a better understanding of the conditions for an efficient advective-acoustic coupling in a decelerated, nonadiabatic flow, in order to extend these results to core-collapse supernovae.

A toy model is analyzed in order to evaluate the linear stability of the gain region immediately behind a stalled accretion shock, after core bounce. This model demonstrates that a negative entropy gradient is not sufficient to warrant linear instability. The stability criterion is governed by the ratio \chi of the advection time through the gain region divided by the local timescale of buoyancy. The gain region is linearly stable if \chi< 3. The classical convective instability is recovered in the limit \chi\gg 3. For \chi>3, perturbations are unstable in a limited range of horizontal wavelengths centered around twice the vertical size H of the gain region. The threshold horizontal wavenumbers k_{min} and k_{max} follow simple scaling laws such that Hk_{min}\propto 1/{\chi} and Hk_{max}\propto \chi. The convective stability of the l=1 mode in spherical accretion is discussed, in relation with the asymmetric explosion of core collapse supernovae. The advective stabilization of long wavelength perturbations weakens the possible influence of convection alone on a global l=1 mode.

The puzzle of birth velocities of pulsars (pulsar kicks) could be solved by an asymmetric explosion of type II Supernovae. We propose a simple hydrodynamical mechanism in order to explain this asymmetry, through the advective-acoustic cycle (Foglizzo 2002) : during the phase of stalled shock, an instability based on the cycle between advected perturbations (entropy / vorticity) and acoustic perturbations can develop between the shock and the surface of the nascent neutron star. Eigenfrequencies are computed numerically, improving the calculation of Houck & Chevalier (1992). The linear instability is dominated by a mode l=1, as observed in the numerical simulations of Blondin et al. (2003) and Scheck et al. (2004). The frequency dependence of the growth rate reveals the presence of the advective-acoustic cycle.

The instability of Bondi-Hoyle-Lyttleton accretion, observed in numerical simulations, is analyzed through known physical mechanisms and possible numerical artefacts. The mechanisms of the longitudinal and transverse instabilities, established within the accretion line model, are clarified. They cannot account for the instability of BHL accretion at moderate Mach number when the pressure forces within the shock cone are taken into account. The advective-acoustic instability is considered in the context of BHL accretion when the shock is detached from the accretor. This mechanism naturally explains the stability of the flow when the shock is weak, and the instability when the accretor is small. In particular, it is a robust proof of the instability of 3D accretion when gamma=5/3 if the accretor is small enough, even for moderate shock strength (M sim 3). The numerical artefacts that may be present in existing numerical simulations are reviewed, with particular attention paid to the advection of entropy/vorticity perturbations and the artificial acoustic feedback from the accretor boundary condition. Several numerical tests are proposed to test these mechanisms.

The linear stability of stalled accretion shocks is investigated in the context of core collapse of type II supernovae. We focus on a particular instability mechanism based on the coupling of acoustic perturbations with advected ones (vorticity, entropy). This advective-acoustic cycle takes place between the shock and the nascent neutron star. Both adiabatic and non-adiabatic processes may contribute to this coupling, but only adiabatic ones are considered in this first approach. The growth time of the adiabatic instability scales like the advection time, and is dominated by low degree modes l=0,1,2. Non radial modes (l=1,2) found unstable by Blondin et al. (2003) can be related to this mechanism.

The linear stability of a shocked isothermal accretion flow onto a black hole is investigated in the inviscid limit. The outer shock solution, which was previously found to be stable with respect to axisymmetric perturbations, is, however, generally unstable to non-axisymmetric ones. Eigenmodes and growth rates are obtained by numerical integration of the linearized equations. The mechanism of this instability is based on the cycle of acoustic waves between their corotation radius and the shock. It is a form of the Papaloizou-Pringle instability, modified by advection and the presence of the shock. As such it can be generalized to non isothermal shocked accretion flows. Blobs and vortices are generated by the shock as a by-product of the instability.

The linear stability of isothermal Bondi accretion with a shock is studied analytically in the asymptotic limit of high incident Mach number M_1. The flow is unstable with respect to radial perturbations as expected by Nakayama (1993), due to post-shock acceleration. Its growth time scales like the advection time from the shock r_sh to the sonic point r_son. The growth rate of non-radial perturbations l=1 is higher by a factor M_1^{2/3}, and is therefore intermediate between the advection and acoustic frequencies. Besides these instabilities based on post-shock acceleration, our study revealed another generic mechanism based on the cycle of acoustic and vortical perturbations between the shock and the sonic radius, independently of the sign of post-shock acceleration. The vortical-acoustic instability is fundamentally non-radial. It is fed by the efficient excitation of vorticity waves by the isothermal shock perturbed by acoustic waves. The growth rate exceeds the advection rate by a factor log M_1. Unstable modes cover a wide range of frequencies from the fundamental acoustic frequency ~c/r_sh up to a cut-off ~c/r_son associated with the sonic radius. The highest growth rate is reached for l=1 modes near the cut-off. The additional cycle of acoustic waves between the shock and the sonic radius is responsible for variations of the growth rate by a factor up to 3 depending on its phase relative to the vortical-acoustic cycle. The instability also exists, with a similar growth rate, below the fundamental acoustic frequency down to the advection frequency, as vorticity waves are efficiently coupled to the region of pseudosound. These results open new perspectives to address the stability of shocked accretion flows.

In the radial flow of gas into a black hole (i.e. Bondi accretion), the infall of any entropy or vorticity perturbation produces acoustic waves propagating outward. The dependence of this acoustic flux on the shape of the perturbation is investigated in detail. This is the key process in the mechanism of the entropic-acoustic instability proposed by Foglizzo & Tagger (2000) to explain the instability of Bondi-Hoyle-Lyttleton accretion. These acoustic waves create new entropy and vorticity perturbations when they reach the shock, thus closing the entropic-acoustic cycle. With an adiabatic index 1<gamma<=5/3, the linearized equations describing the perturbations of the Bondi flow are studied analytically and solved numerically. The fundamental frequency of this problem is the cut-off frequency of acoustic refraction, below which ingoing acoustic waves are refracted out. This cut-off is significantly smaller than the Keplerian frequency at the sonic radius and depends on the latitudinal number l of the perturbations. When advected adiabatically inward, entropy and vorticity perturbations trigger acoustic waves propagating outward, with an efficiency which is highest for non radial perturbations l=1. The outgoing acoustic flux produced by the advection of vorticity perturbations is always moderate and peaks at rather low frequency. By contrast, the acoustic flux produced by an entropy wave is highest close to the refraction cut-off. It can be very large if gamma is close to 5/3. These results suggest that the shocked Bondi flow with gamma=5/3 is strongly unstable with respect to the entropic-acoustic mechanism.

A new instability mechanism is described in accretion flows where the gas is accelerated from a stationary shock to a sonic surface. The instability is based on a cycle of acoustic and entropic waves in this subsonic region of the flow. When advected adiabatically inward, entropy perturbations trigger acoustic waves propagating outward. If a shock is present at the outer boundary, acoustic waves reaching the shock produce new entropy perturbations, thus creating an entropic-acoustic cycle between the shock and the sonic surface. The interplay of acoustic and entropy perturbations is estimated analytically using a simplified model based on the compact nozzle approximation. According to this model, the entropic-acoustic cycle is unstable if the sound speed at the sonic surface significantly exceeds the sound speed immediately after the shock. The growth rate scales like the inverse of the advection time from the outer shock to the sonic point. The frequency of the most unstable perturbations is comparable to the refraction cutoff, defined as the frequency below which acoustic waves propagating inward are significantly refracted outward. This generic mechanism should occur in Bondi-Hoyle-Lyttleton accretion, and also in shocked accretion discs.

The adiabatic shock produced by a compact object moving supersonically relative to a gas with uniform entropy and no vorticity is a source of entropy gradients and vorticity. We investigate these analytically. The non-axisymmetric Rayleigh-Taylor and axisymmetric Kelvin-Helmholtz linear instabilities are potential sources of destabilization of the subsonic accretion flow after the shock. A local Lagrangian approach is used in order to evaluate the efficiency of these linear instabilities. However, the conditions required for such a WKB type approximation are fulfilled only marginally: a quantitative estimate of their local growth rate integrated along a flow line shows that their growth time is at best comparable to the time needed for advection onto the accretor, even at high Mach number and for a small accretor size. Despite this apparently low efficiency, several features of these mechanisms qualitatively match those observed in numerical simulations: in a gas with uniform entropy, the instability occurs only for supersonic accretors. It is nonaxisymmetric, and begins close to the accretor in the equatorial region perpendicular to the symmetry axis. The mechanism is more efficient for a small, highly supersonic accretor, and also if the shock is detached. We also show by a 3-D numerical simulation an example of unstable accretion of a subsonic flow with non-uniform entropy at infinity. This instability is qualitatively similar to the one observed in 3-D simulations of the Bondi-Hoyle-Lyttleton flow, although it involves neither a bow shock nor an accretion line.

The excitation mechanism of low degree acoustic modes is investigated through the analysis of the stochastic time variations of their energy. The correlation between the energies of two different modes is interpreted as the signature of the occurrence rate of their excitation source. The different correlations determined by Foglizzo et al. (1998) constrain the physical properties of an hypothetical source of excitation, which would act in addition to the classical excitation by the turbulent convection. Particular attention is drawn to the effect of coronal mass ejections. The variability of their occurrence rate with the solar cycle could account for the variation of the correlation between IPHIR and GOLF data. Such an interpretation would suggest that the mean correlation between low degree acoustic modes is at least 0.2 % at solar minimum.

We have studied the statistical properties of the energy of individual acoustic modes, extracted from 310 days of GOLF data near the solar minimum. The exponential distribution of the energy of each mode is clearly seen. The modes are found to be uncorrelated with a +-0.6 % accuracy, thus supporting the hypothesis of stochastic excitation by the solar convection. Nevertheless, the same analysis performed on the same modes just before the solar maximum, using IPHIR data, rejects the hypothesis of no correlation at a 99.3 % confidence level. A simple model suggests that 31.3 +-9.4 % of the energy of each mode is coherent among the modes studied in IPHIR data, correponding to a mean correlation of 10.7 +-5.9 %.

We prove that the sonic surface of axisymmetric meridional stationary flows is always attached to the accretor, however small, if the adiabatic index of the gas is gamma=5/3. Using local expansions near a point-like accretor, we extend Bondi's classification of spherically symmetric flows to axisymmetric flows, introducing the possibility of angular sectors reached by no flow lines, and singular directions of infinite mass flux, in addition to the angular regions of subsonic and supersonic accretion. For gamma<5/3, we show the impossibility of subsonic accretion onto a point-like accretor when the entropy of the flow is not uniform. The special case gamma=5/3 is treated separately. We analyse the influence of the adiabatic index and Mach number of the flow at infinity on the mass accretion rate of shocked spherical flows. We propose an interpolation formula for the mass accretion rate of axisymmetric flows as a function of the Mach number and the adiabatic index, in the range 9/7<gamma<5/3.

We describe the effects of both magnetic buoyancy and differential rotation on a disc of isothermal gas embedded in a purely azimuthal magnetic field, in order to study the evolution and interplay of Parker and shearing instabilities. Both instabilities occur on the slow MHD branch of the dispersion relation, and can affect the same waves. We put a stress on the natural polarization properties of the slow MHD waves to get a better understanding of the physics involved. The mechanism of the shearing instability is described in details. Waves with a vertical wavelength shorter than the scale height of the disc are subject to both the Parker and the transient shearing instabilities. They occur in different ranges of radial wavenumbers, \ie at different times in the shearing evolution ; these ranges can overlap or, on the contrary, be separated by a phase of wave-like oscillations, depending on the strength of differential rotation. Our results can be applied to both galactic and accretion discs.

We present a detailed study of the growth of the Parker instability in a differentially rotating disk embedded in an azimuthal equilibrium magnetic field, such as the interstellar gas or an accretion disk. Basic properties of the instability without shear are first recalled. Differential rotation is modeled in the shearing sheet approximation, classical in the theory of spiral density waves. The action of differential rotation is reduced to two different effects, (i) a linear time-dependence of the radial wavenumber, and (ii) a radial differential force. We present both exact numerical solutions, and approximate analytical ones based on the WKB approximation in the limit of weak differential rotation. Most important are (i) a transient natural stabilization of the Parker mode due to the radial differential force (ii) the generation of magnetosonic and Alfvenic waves, and (iii) in a certain parameter range a possible ``turn-over'' of the perturbation whereby, quite surprisingly, matter which had started being elevated by the instability may end up dropping towards the disk midplane. A simplified model shows the possible observable effects of this turn-over.