Context. Classical novae are thermonuclear explosions hosted by accreting white dwarfs in stellar binary systems. Material piles up on top of the white dwarf star under mildly degenerate conditions, driving a thermonuclear runaway. The energy released by the suite of nuclear processes operating at the envelope (mostly proton-capture reactions and beta-decays) heats the material up to peak temperatures ranging from 100 to 400 MK. In these events, about 10-3 - 10-7 Msun, enriched in CNO and, sometimes, other intermediate-mass elements (e.g., Ne, Na, Mg, Al), are ejected into the interstellar medium. Aims. To date, most of the efforts undertaken in the modeling of classical nova outbursts have focused on the early stages of the explosion and ejection, ignoring the interaction of the ejecta, first with the accretion disk orbiting the white dwarf, and ultimately with the secondary star. Methods. A suite of three-dimensional, SPH simulations of the interaction between the nova ejecta, the accretion disk, and the stellar companion have been performed to fill this gap, aimed at testing the influence of the different parameters (i.e., mass and velocity of the ejecta, mass and geometry of the accretion disk) on the dynamical and chemical properties of the system. Results. We discuss the conditions that lead to the disruption of the accretion disk and to mass loss from the binary system. In addition, we discuss the likelihood of chemical contamination of the stellar secondary induced by the impact with the nova ejecta and its potential effect on the next nova cycle.

Hydrodynamical instabilities and shocks are ubiquitous in astrophysical scenarios. Therefore, an accurate numerical simulation of these phenomena is mandatory to correctly model and understand many astrophysical events, such as Supernovas, stellar collisions, or planetary formation. In this work, we attempt to address many of the problems that the smoothed particle hydrodynamics (SPH) technique has when dealing with subsonic hydrodynamical instabilities or shocks. To that aim we built a new SPH code named SPHYNX, that includes many of the recent advances in the SPH technique and some other new ones, which we present here. SPHYNX is of Newtonian type and grounded in the Euler-Lagrange formulation of the SPH technique. Its distinctive features are: the use of an integral approach to estimating the gradients; the use of a flexible family of interpolators called sinc kernels, which suppress pairing instability; and the incorporation of a new type of volume element which provides a better partition of the unity. Unlike other modern formulations, which consider volume elements linked to pressure, our volume element choice relies on density. SPHYNX conserves mass, linear and angular momentum, energy, entropy, and preserves kernel normalization even in strong shocks. The coupling between the integral approach to calculate gradients and the new family of volume elements reduces the so-called tensile instability. Both features help to suppress the damp which often prevents the growth of hydrodynamic instabilities in regular SPH codes. On the whole, SPHYNX has passed the verification tests described below. For identical particle setting and initial conditions the results were similar (or better in some particular cases) than those obtained with other SPH schemes such as GADGET-2, PSPH or with the recent density-independent formulation (DISPH) and conservative reproducing kernel (CRKSPH) techniques.

SN2011fe was detected by the Palomar Transient Factory on August 24th 2011 in M101 few hours after the explosion. From the early spectra it was immediately realized that it was a Type Ia supernova thus making this event the brightest one discovered in the last twenty years. In this paper the observations performed with the instruments on board of INTEGRAL (SPI, IBIS/ISGRI, JEM-X and OMC) before and after the maximum of the optical light as well as the interpretation in terms of the existing models of $\gamma$--ray emission from such kind of supernovae are reported. All INTEGRAL high-energy have only been able to provide upper limits to the expected emission due to the decay of $^{56}$Ni. These bounds allow to reject explosions involving a massive white dwarf in the sub--Chandrasekhar scenario. On the other hand, the optical light curve obtained with the OMC camera suggests that the event was produced by a delayed detonation of a CO white dwarf that produced $\sim 0.5$ M$\odot$ of $^{56}$Ni. In this particular case, INTEGRAL would have only been able to detect the early $\gamma$--ray emission if the supernova had occurred at a distance of 2 -3 Mpc, although the brightest event could be visible up to distances larger by a factor two.

SN2011fe was detected by the Palomar Transient Factory on August 24th 2011 in M101 a few hours after the explosion. From the early optical spectra it was immediately realized that it was a Type Ia supernova thus making this event the brightest one discovered in the last twenty years. The distance of the event offered the rare opportunity to perform a detailed observation with the instruments on board of INTEGRAL to detect the gamma-ray emission expected from the decay chains of $^{56}$Ni. The observations were performed in two runs, one before and around the optical maximum, aimed to detect the early emission from the decay of $^{56}$Ni and another after this maximum aimed to detect the emission of $^{56}$Co. The observations performed with the instruments on board of INTEGRAL (SPI, IBIS/ISGRI, JEMX and OMC) have been analyzed and compared with the existing models of gamma-ray emission from such kind of supernovae. In this paper, the analysis of the gamma-ray emission has been restricted to the first epoch. Both, SPI and IBIS/ISGRI, only provide upper-limits to the expected emission due to the decay of $^{56}$Ni. These upper-limits on the gamma-ray flux are of 7.1 $\times$ 10$^{-5}$ ph/s/cm$^2$ for the 158 keV line and of 2.3 $\times$ 10$^{-4}$ ph/s/cm$^2$ for the 812 keV line. These bounds allow to reject at the $2\sigma$ level explosions involving a massive white dwarf, $\sim 1$ M$\odot$ in the sub--Chandrasekhar scenario and specifically all models that would have substantial amounts of radioactive $^{56}$Ni in the outer layers of the exploding star responsible of the SN2011fe event. The optical light curve obtained with the OMC camera also suggests that SN2011fe was the outcome of the explosion, possibly a delayed detonation although other models are possible, of a CO white dwarf that synthesized $\sim 0.55$ M$_\odot$ of $^{56}$Ni. For this specific model.

The behavior of IAD_0 scheme, a fully conservative SPH scheme based on a tensor formulation, is analyzed in connection with several astrophysical scenarios, and compared to the same simulations carried out with the standard SPH technique. The proposed hydrodynamic scheme is validated using a variety of numerical tests that cover important topics in astrophysics, such as the evolution of supernova remnants, the stability of self-gravitating bodies and the coalescence of compact objects. The results suggest that the SPH scheme built with the integral approach to the derivatives premise improves the results of the standard SPH technique. In particular, it is observed a better development of hydrodynamic instabilities, an improved description of self-gravitant structures in equilibrium and a reasonable description of the process of coalescence of two white dwarfs. A good energy, and linear and angular momentum conservation, generally better than that of standard SPH, was also obtained. In addition the new scheme is less susceptible to suffer pairing instability.

In this paper we develop and test a fully conservative SPH scheme based on a tensor formulation that can be applied to simulate astrophysical systems. In the proposed scheme, derivatives are calculated from an integral expression that leads to a tensor (instead of a vectorial) estimation of gradients and reduces to the standard formulation in the continuum limit. The new formulation improves the interpolation of physical magnitudes, leading to a set of conservative equations that resembles those of standard SPH. The resulting scheme is verified using a variety of well-known tests, all of them simulated in two dimensions. We also discuss an application of the proposed tensor method to astrophysics by simulating the stability of a Sun-like polytrope calculated in three dimensions.

Observational evidences point to a common explosion mechanism of Type Ia supernovae based on a delayed detonation of a white dwarf. However, all attempts to find a convincing ignition mechanism based on a delayed detonation in a destabilized, expanding, white dwarf have been elusive so far. One of the possibilities that has been invoked is that an inefficient deflagration leads to pulsation of a Chandrasekhar-mass white dwarf, followed by formation of an accretion shock that confines a carbon-oxygen rich core, while transforming the kinetic energy of the collapsing halo into thermal energy of the core, until an inward moving detonation is formed. This chain of events has been termed Pulsating Reverse Detonation (PRD). In this work we present three dimensional numerical simulations of PRD models from the time of detonation initiation up to homologous expansion. Different models characterized by the amount of mass burned during the deflagration phase, M_defl, give explosions spanning a range of kinetic energies, K ~ (1.0-1.2) foes, and 56Ni masses, M(56Ni) ~ 0.6-0.8 M_sun, which are compatible with what is expected for typical Type Ia supernovae. Spectra and light curves of angle-averaged spherically symmetric versions of the PRD models are discussed. Type Ia supernova spectra pose the most stringent requirements on PRD models.

Observational evidences point to a common explosion mechanism of Type Ia supernovae based on a delayed detonation of a white dwarf. Although several scenarios have been proposed and explored by means of one, two, and three-dimensional simulations, the key point still is the understanding of the conditions under which a stable detonation can form in a destabilized white dwarf. One of the possibilities that have been invoked is that an inefficient deflagration leads to the pulsation of a Chandrasekhar-mass white dwarf, followed by formation of an accretion shock around a carbon-oxygen rich core. The accretion shock confines the core and transforms kinetic energy from the collapsing halo into thermal energy of the core, until an inward moving detonation is formed. This chain of events has been termed Pulsating Reverse Detonation (PRD). In this work we explore the robustness of the detonation ignition for different PRD models characterized by the amount of mass burned during the deflagration phase, M_defl. The evolution of the white dwarf up to the formation of the accretion shock has been followed with a three-dimensional hydrodynamical code with nuclear reactions turned off. We found that detonation conditions are achieved for a wide range of M_defl. However, if the nuclear energy released during the deflagration phase is close to the white dwarf binding energy (~ 0.46 foes -> M_defl ~ 0.30 M_sun) the accretion shock cannot heat and confine efficiently the core and detonation conditions are not robustly achieved.

The axisymmetric form of the hydrodynamic equations within the smoothed particle hydrodynamics (SPH) formalism is presented and checked using idealized scenarios taken from astrophysics (free fall collapse, implosion and further pulsation of a sun-like star), gas dynamics (wall heating problem, collision of two streams of gas) and inertial confinement fusion (ICF, -ablative implosion of a small capsule-). New material concerning the standard SPH formalism is given. That includes the numerical handling of those mass points which move close to the singularity axis, more accurate expressions for the artificial viscosity and the heat conduction term and an easy way to incorporate self-gravity in the simulations. The algorithm developed to compute gravity does not rely in any sort of grid, leading to a numerical scheme totally compatible with the lagrangian nature of the SPH equations.

A set of interpolating functions of the type f(v)={(sin[v pi/2])/(v pi/2)}^n is analyzed in the context of the smoothed-particle hydrodynamics (SPH) technique. The behaviour of these kernels for several values of the parameter n has been studied either analytically as well as numerically in connection with several tests carried out in two dimensions. The main advantage of this kernel relies in its flexibility because for n=3 it is similar to the standard widely used cubic-spline, whereas for n>3 the interpolating function becomes more centrally condensed, being well suited to track discontinuities such as shock fronts and thermal waves.

Deflagration models poorly explain the observed diversity of SNIa. Current multidimensional simulations of SNIa predict a significant amount of, so far unobserved, carbon and oxygen moving at low velocities. It has been proposed that these drawbacks can be resolved if there is a sudden jump to a detonation (delayed detonation), but this kind of models has been explored mainly in one dimension. Here we present new three-dimensional delayed detonation models in which the deflagraton-to-detonation transition (DDT) takes place in conditions like those favored by one-dimensional models. We have used a SPH code adapted to SNIa with algorithms devised to handle subsonic as well as supersonic combustion fronts. The starting point was a C-O white dwarf of 1.38 solar masses. When the average density on the flame surface reached 2-3x10^7 g/cm^3 a detonation was launched. The detonation wave processed more than 0.3 solar masses of carbon and oxygen, emptying the central regions of the ejecta of unburned fuel and raising its kinetic energy close to the fiducial 10^51 ergs expected from a healthy Type Ia supernova. The final amount of 56Ni synthesized also was in the correct range. However, the mass of carbon and oxygen ejected is still too high. The three-dimensional delayed detonation models explored here show an improvement over pure deflagration models, but they still fail to coincide with basic observational constraints. However, there are many aspects of the model that are still poorly known (geometry of flame ignition, mechanism of DDT, properties of detonation waves traversing a mixture of fuel and ashes). Therefore, it will be worth pursuing its exploration to see if a good SNIa model based on the three-dimensional delayed detonation scenario can be obtained.

We calculate detailed NLTE synthetic spectra of a Pulsating Reverse Detonation (PRD) model, a novel explosion mechanism for Type Ia supernovae. While the hydro models are calculated in 3-D, the spectra use an angle averaged hydro model and thus some of the 3-D details are lost, but the overall average should be a good representation of the average observed spectra. We study the model at 3 epochs: maximum light, seven days prior to maximum light, and 5 days after maximum light. At maximum the defining Si II feature is prominent, but there is also a prominent C II feature, not usually observed in normal SNe Ia near maximum. We compare to the early spectrum of SN 2006D which did show a prominent C II feature, but the fit to the observations is not compelling. Finally we compare to the post-maximum UV+optical spectrum of SN 1992A. With the broad spectral coverage it is clear that the iron-peak elements on the outside of the model push too much flux to the red and thus the particular PRD realizations studied would be intrinsically far redder than observed SNe Ia. We briefly discuss variations that could improve future PRD models.

The physical structure of a nuclear flame is a basic ingredient of the theory of Type Ia supernovae (SNIa). Assuming an exponential density reduction with several characteristic times we have followed the evolution of a planar nuclear flame in an expanding background from an initial density 6.6 10^7 g/cm3 down to 2 10^6 g/cm3. The total amount of synthesized intermediate-mass elements (IME), from silicon to calcium, was monitored during the calculation. We have made use of the computed mass fractions, X_IME, of these elements to give an estimation of the total amount of IME synthesized during the deflagration of a massive white dwarf. Using X_IME and adopting the usual hypothesis that turbulence decouples the effective burning velocity from the laminar flame speed, so that the relevant flame speed is actually the turbulent speed on the integral length-scale, we have built a simple geometrical approach to model the region where IME are thought to be produced. It turns out that a healthy production of IME involves the combination of not too short expansion times, t_c > 0.2 s, and high turbulent intensities. According to our results it could be difficult to produce much more than 0.2 solar masses of intermediate-mass elements within the deflagrative paradigma. The calculations also suggest that the mass of IME scales with the mass of Fe-peak elements, making it difficult to conciliate energetic explosions with low ejected nickel masses, as in the well observed SN1991bg or in SN1998de. Thus a large production of Si-peak elements, especially in combination with a low or a moderate production of iron, could be better addressed by either the delayed detonation route in standard Chandrasekhar-mass models or, perhaps, by the off-center helium detonation in the sub Chandrasekhar-mass scenario.

We describe a mechanism by which a failed deflagration of a Chandrasekhar-mass carbon-oxygen white dwarf can turn into a successful thermonuclear supernova explosion, without invoking an ad hoc high-density deflagration-detonation transition. Following a pulsating phase, an accretion shock develops above a core of 1 M_sun composed of carbon and oxygen, inducing a converging detonation. A three-dimensional simulation of the explosion produced a kinetic energy of 1.05E51 ergs and 0.70 M_sun of 56Ni, ejecting scarcely 0.01 M_sun of C-O moving at low velocities. The mechanism works under quite general conditions and is flexible enough to account for the diversity of normal Type Ia supernovae. In given conditions the detonation might not occur, which would reflect in peculiar signatures in the gamma and UV-wavelengths

In this paper, we review the present state of theoretical models of thermonuclear supernovae, and compare their predicitions with the constraints derived from observations of Type Ia supernovae. The diversity of explosion mechanisms usually found in one-dimensional simulations is a direct consequence of the impossibility to resolve the flame structure under the assumption of spherical symmetry. Spherically symmetric models have been successful in explaining many of the observational features of Type Ia supernovae, but they rely on two kinds of empirical models: one that describes the behaviour of the flame on the scales unresolved by the code, and another that takes account of the evolution of the flame shape. In contrast, three-dimensional simulations are able to compute the flame shape in a self-consistent way, but they still need a model for the propagation of the flame in the scales unresolved by the code. Furthermore, in three dimensions the number of degrees of freedom of the initial configuration of the white dwarf at runaway is much larger than in one dimension. Recent simulations have shown that the sensitivity of the explosion output to the initial conditions can be extremely large. New paradigms of thermonuclear supernovae have emerged from this situation, as the Pulsating Reverse Detonation. The resolution of all these issues must rely on the predictions of observational properties of the models, and their comparison with current Type Ia supernova data, including X-ray spectra of Type Ia supernova remnants.

In this paper we address the theory of Type Ia supernovae from the moment of carbon runaway up to several hours after the explosion. We have concentrated on the boiling-pot model: a deflagration characterized by the (nearly-) simultaneous ignition of a number of bubbles that pervade the core of the white dwarf. Thermal fluctuations larger than >1% of the background temperature (7x10^8 K) on lengthscales of < 1m could be the seeds of the bubbles. Variations of the homogeneity of the temperature perturbations can lead to two alternative configurations at carbon runaway: if the thermal gradient is small, all the bubbles grow to a common characteristic size related to the value of the thermal gradient, but if the thermal gradient is large enough, the size spectrum of the bubbles extends over several orders of magnitude. The explosion phase has been studied with the aid of a smoothed particle hydrodynamics code suited to simulate thermonuclear supernovae. In spite of important procedural differences and different physical assumptions, our results converge with the most recent calculations of 3D deflagrations in white dwarfs carried out in supernova studies by different groups. For large initial numbers of bubbles (>3-4 per octant), the explosion produces about 0.45 solar masses of 56Ni, and the kinetic energy of the ejecta is 0.45x10^{51} ergs. However, all three-dimensional deflagration models share three main drawbacks: 1) the scarce synthesis of intermediate-mass elements, 2) the loss of chemical stratification of the ejecta due to mixing by Rayleigh-Taylor instabilities during the first second of the explosion, and 3) the presence of big clumps of 56Ni at the photosphere at the time of maximum brightness.

We have explored in three dimensions the fate of a massive white dwarf as a function of different initial locations of carbon ignition, with the aid of a SPH code. The calculated models cover a variety of possibilities ranging from the simultaneous ignition of the central volume of the star to the off-center ignition in multiple scattered spots. In the former case, there are discussed the possibility of a transition to a detonation when the mean density of the nuclear flame decreases below 2x10**7 g cm**-3, and its consequences. In the last case, the dependence of the results on the number of initial igniting spots and the chance of some of these models to evolve to the pulsating delayed detonation scenario are also outlined.

A method for integrating the chemical equations associated with nuclear combustion at high temperature is presented and extensively checked. Following the idea of E. M\"uller, the feedback between nuclear rates and temperature was taken into account by simultaneously computing molar fraction changes and temperature response in the same matrix. The resulting algorithm is very stable and efficient at calculating nuclear combustion in explosive scenarios, especially in those situations where the reacting material manages to climb to the nuclear statistical equilibrium regime. The numerical scheme may be useful not only for those who carry out hydrodynamical simulations of explosive events, but also as a tool to investigate the properties of a nuclear system approaching equilibrium through a variety of thermodynamical trajectories.

In recent years, it has become clear from multidimensional simulations that the outcome of deflagrations depends strongly on the initial configuration of the flame. We have studied under which conditions this configuration could consist of a number of scattered, isolated, hot bubbles. Afterwards, we have calculated the evolution of deflagrations starting from different numbers of bubbles. We have found that starting from 30 bubbles a mild explosion is produced M(Ni56)=0.56 solar masses, while starting from 10 bubbles the star becomes only marginally unbound (K = 0.05 foes).

We present several hydrodynamical simulations of thermonuclear supernovae dealing with multiple delayed detonations. The calculations were carried out in three dimensions, making possible to study the influence of geometry of the flame front in two aspects. First, the evolution of its fractal dimension during the deflagration phase has been followed until a critical value is reached such that the deflagration may turn into a detonation. Second, as the resulting detonation could probably be scattered through the flame, the effect of its initial location on the detonation propagation, final energetics and nucleosynthesis has been explored.

The explosion of a white dwarf of mass 1.36 M$_\odot$ has been simulated in three dimensions with the aid of a SPH code. The explosion follows the delayed detonation paradigma. In this case the deflagration-detonation transition is induced by the large corrugation of the flame front resulting from Rayleigh-Taylor instability and turbulence. The nucleosynthetic yields have been calculated, showing that some neutronized isotopes such as 54^Fe or 58^Ni are not overproduced with respect to the solar system ratios. The distribution of intermediate-mass elements is also compatible with the spectra of normal SNIa. The excepcion is, however, the abundance of carbon and oxygen, which are overproduced.

We present a straightforward integration method to compute the abundance and temperature evolution in explosive scenarios. In this approach the thermal equation is implicitely coupled with chemical equations in order to avoid instabilities and ensure a gentle transition from the normal combustion regime to the quasi (QSE) and complete nuclear statistical equilibrium (NSE). Two nuclear networks, with 14 nuclei (alpha-network) and 86 nuclei (including protons and neutrons) respectively, have been considered. The scheme is suitable to cope with a variety of explosive burning regimes.