The faster light-curve evolution of low-luminosity Type Ia supernovae (SNe Ia) suggests that they could result from the explosion of white dwarf (WD) progenitors below the Chandrasekhar mass ($M_{\rm Ch}$). Here we present 1D non-LTE time-dependent radiative transfer simulations of pure central detonations of carbon-oxygen WDs with a mass ($M_\rm{tot}$) between 0.88 M$_{\odot}$ and 1.15 M$_{\odot}$, and a $^{56}\rm{Ni}$ yield between 0.08 M$_{\odot}$ and 0.84 M$_{\odot}$. Their lower ejecta density compared to $M_{\rm Ch}$ models results in a more rapid increase of the luminosity at early times and an enhanced $\gamma$-ray escape fraction past maximum light. Consequently, their bolometric light curves display shorter rise times and larger post-maximum decline rates. Moreover, the higher $M(^{56}\rm{Ni})/M_\rm{tot}$ ratio at a given $^{56}\rm{Ni}$ mass enhances the temperature and ionization level in the spectrum-formation region for the less luminous models, giving rise to bluer colours at maximum light and a faster post-maximum evolution of the $B-V$ colour. For sub-$M_{\rm Ch}$ models fainter than $M_B\approx -18.5$ mag at peak, the greater bolometric decline and faster colour evolution lead to a larger $B$-band post-maximum decline rate, $\Delta M_{15}(B)$. In particular, all of our previously-published $M_{\rm Ch}$ models (standard and pulsational delayed detonations) are confined to $\Delta M_{15}(B) < 1.4$ mag, while the sub-$M_{\rm Ch}$ models with $M_\rm{tot}\lesssim 1$ M$_{\odot}$ extend beyond this limit to $\Delta M_{15}(B)\approx 1.65$ mag for a peak $M_B\approx -17$ mag, in better agreement with the observed width-luminosity relation (WLR). Regardless of the precise ignition mechanism, these simulations suggest that fast-declining SNe Ia at the faint end of the WLR could result from the explosion of WDs whose mass is significantly below the Chandrasekhar limit.

We explore the physics of SN Ia light curves and spectra using the 1-D non-LTE time-dependent radiative-transfer code CMFGEN. Rather than adjusting ejecta properties to match observations, we select as input one "standard" 1-D Chandrasekhar-mass delayed-detonation hydrodynamical model, and then explore the sensitivity of radiation and gas properties on radiative-transfer modeling assumptions. The correct computation of SN Ia radiation is not exclusively a solution to an "opacity problem", characterized by the treatment of a large number of lines. It is also key to treat important atomic processes consistently. Besides handling line blanketing in non-LTE, we show that including forbidden line transitions of metals is increasingly important for the temperature and ionization of the gas beyond maximum light. Non-thermal ionization and excitation are also critical since they affect the color evolution and the Delta-M15 of our model. While impacting little the bolometric luminosity, a more complete treatment of decay routes leads to enhanced line blanketing, e.g., associated with 48Ti in the U and B bands. Overall, we find that SN Ia radiation properties are influenced in a complicated way by the atomic data we employ, so that obtaining converged results is a challenge. We nonetheless obtain a good match to the golden standard type Ia SN 2005cf in the optical and near-IR, from 5 to 60d after explosion, suggesting that assuming spherical symmetry is not detrimental to SN Ia radiative-transfer modeling at these times. Multi-D effects no doubt matter, but they are perhaps less important than accurately treating non-LTE processes [abridged].

Observations of SN 2011fe at early times reveal an evolution analogous to a fireball model of constant color. In contrast, our unmixed delayed detonations of Chandrasekhar-mass white dwarfs (DDC series) exhibit a faster brightening concomitant with a shift in color to the blue. In this paper, we study the origin of these discrepancies. We find that strong chemical mixing largely resolves the photometric mismatch at early times, but it leads to an enhanced line broadening that contrasts, for example, with the markedly narrow SiII6355A line of SN 2011fe. We also explore an alternative configuration with pulsational-delayed detonations (PDDEL model series). Because of the pulsation, PDDEL models retain more unburnt carbon, have little mass at high velocity, and have a much hotter outer ejecta after the explosion. The pulsation does not influence the inner ejecta, so PDDEL and DDC models exhibit similar radiative properties beyond maximum. However, at early times, PDDEL models show bluer optical colors and a higher luminosity, even for weak mixing. Their early-time radiation is derived primarily from the initial shock-deposited energy in the outer ejecta rather than radioactive decay heating. Furthermore, PDDEL models show short-lived CII lines, reminiscent of SN 2013dy. They typically exhibit lines that are weaker, narrower, and of near-constant width, reminiscent of SN 2011fe. In addition to multi-dimensional effects, varying configurations for such ``pulsations" offer a source of spectral diversity amongst SNe Ia. PDDEL and DDC models also provide one explanation for low- and high-velocity gradient SNe Ia.

The high metal content and fast expansion of supernova (SN) Ia ejecta lead to considerable line overlap in their optical spectra. Uncertainties in composition and ionization further complicate the process of line identification. In this paper, we focus on the 5900A emission feature seen in SN Ia spectra after bolometric maximum, a line which in the last two decades has been associated with [CoIII]5888A or NaID. Using non-LTE time-dependent radiative-transfer calculations based on Chandrasekhar-mass delayed-detonation models, we find that NaID line emission is extremely weak at all post-maximum epochs. Instead, we predict the presence of [CoIII]5888A after maximum in all our SN Ia models, which cover a range from 0.12 to 0.87Msun of 56Ni. We also find that the [CoIII]5888A forbidden line is present within days of bolometric maximum, and strengthens steadily for weeks thereafter. Both predictions are confirmed by observations. Rather than trivial taxonomy, these findings confirm that it is necessary to include forbidden-line transitions in radiative-transfer simulations of SNe Ia, both to obtain the correct ejecta cooling rate and to match observed optical spectra.

Aims. The Sun shows strong variability in its magnetic activity, from Grand minima to Grand maxima, but the nature of the variability is not fully understood, mostly because of the insufficient length of the directly observed solar activity records and of uncertainties related to long-term reconstructions. Here we present a new adjustment-free reconstruction of solar activity over three millennia and study its different modes. Methods. We present a new adjustment-free, physical reconstruction of solar activity over the past three millennia, using the latest verified carbon cycle, 14C production, and archeomagnetic field models. This great improvement allowed us to study different modes of solar activity at an unprecedented level of details. Results. The distribution of solar activity is clearly bi-modal, implying the existence of distinct modes of activity. The main regular activity mode corresponds to moderate activity that varies in a relatively narrow band between sunspot numbers about 20 and 67. The existence of a separate Grand minimum mode with reduced solar activity, which cannot be explained by random fluctuations of the regular mode, is confirmed at a high confidence level. The possible existence of a separate Grand maximum mode is also suggested, but the statistics is too low to reach a confident conclusion. Conclusions. The Sun is shown to operate in distinct modes - a main general mode, a Grand minimum mode corresponding to an inactive Sun, and a possible Grand maximum mode corresponding to an unusually active Sun. These results provide important constraints for both dynamo models of Sun-like stars and investigations of possible solar influence on Earth's climate.

The delayed-detonation explosion mechanism applied to a Chandrasekhar-mass white dwarf offers a very attractive model to explain the inferred characteristics of Type Ia supernovae (SNe Ia). The resulting ejecta are chemically stratified, have the same mass and roughly the same asymptotic kinetic energy, but exhibit a range in 56Ni mass. We investigate the contemporaneous photometric and spectroscopic properties of a sequence of delayed-detonation models, characterized by 56Ni masses between 0.18 and 0.81 Msun. Starting at 1d after explosion, we perform the full non-LTE, time-dependent radiative transfer with the code CMFGEN, with an accurate treatment of line blanketing, and compare our results to SNe Ia at bolometric maximum. Despite the 1D treatment, our approach delivers an excellent agreement to observations. We recover the range of SN Ia luminosities, colours, and spectral characteristics from the near-UV to 1 micron, for standard as well as low-luminosity 91bg-like SNe Ia. Our models predict an increase in rise time to peak with increasing 56Ni mass, from ~15 to ~21d, yield peak bolometric luminosities that match Arnett's rule to within 10%, and reproduce the much smaller scatter in near-IR magnitudes compared to the optical. We reproduce the morphology of individual spectral features, the stiff dependence of the R(Si) spectroscopic ratio on 56Ni mass, and the onset of blanketing from TiII/ScII in low-luminosity SNe Ia with a 56Ni mass <0.3 Msun. We find that ionization effects, which often dominate over abundance variations, can produce high-velocity features in CaII lines, even in 1D. Distinguishing between different SN Ia explosion mechanisms is a considerable challenge but the results presented here provide additional support to the viability of the delayed-detonation model.

We present a study of the peculiar Type Ia supernova 2001ay (SN 2001ay). The defining features of its peculiarity are: high velocity, broad lines, and a fast rising light curve, combined with the slowest known rate of decline. It is one magnitude dimmer than would be predicted from its observed value of Delta-m15, and shows broad spectral features. We base our analysis on detailed calculations for the explosion, light curves, and spectra. We demonstrate that consistency is key for both validating the models and probing the underlying physics. We show that this SN can be understood within the physics underlying the Delta-m15 relation, and in the framework of pulsating delayed detonation models originating from a Chandrasekhar mass, white dwarf, but with a progenitor core composed of 80% carbon. We suggest a possible scenario for stellar evolution which leads to such a progenitor. We show that the unusual light curve decline can be understood with the same physics as has been used to understand the Delta-m15 relation for normal SNe Ia. The decline relation can be explained by a combination of the temperature dependence of the opacity and excess or deficit of the peak luminosity, alpha, measured relative to the instantaneous rate of radiative decay energy generation. What differentiates SN 2001ay from normal SNe Ia is a higher explosion energy which leads to a shift of the Ni56 distribution towards higher velocity and alpha < 1. This result is responsible for the fast rise and slow decline. We define a class of SN 2001ay-like SNe Ia, which will show an anti-Phillips relation.

We present ultraviolet (UV) spectroscopy and photometry of four Type Ia supernovae (SNe 2004dt, 2004ef, 2005M, and 2005cf) obtained with the UV prism of the Advanced Camera for Surveys on the Hubble Space Telescope. This dataset provides unique spectral time series down to 2000 Angstrom. Significant diversity is seen in the near maximum-light spectra (~ 2000--3500 Angstrom) for this small sample. The corresponding photometric data, together with archival data from Swift Ultraviolet/Optical Telescope observations, provide further evidence of increased dispersion in the UV emission with respect to the optical. The peak luminosities measured in uvw1/F250W are found to correlate with the B-band light-curve shape parameter dm15(B), but with much larger scatter relative to the correlation in the broad-band B band (e.g., ~0.4 mag versus ~0.2 mag for those with 0.8 < dm15 < 1.7 mag). SN 2004dt is found as an outlier of this correlation (at > 3 sigma), being brighter than normal SNe Ia such as SN 2005cf by ~0.9 mag and ~2.0 mag in the uvw1/F250W and uvm2/F220W filters, respectively. We show that different progenitor metallicity or line-expansion velocities alone cannot explain such a large discrepancy. Viewing-angle effects, such as due to an asymmetric explosion, may have a significant influence on the flux emitted in the UV region. Detailed modeling is needed to disentangle and quantify the above effects.

Incomplete carbon-oxygen detonation with reactions terminating after burning of 12C in the leading 12C + 12C reaction (C-detonation) may occur in the low density outer layers of white dwarfs exploding as Type Ia supernovae (SNIa). Previous studies of carbon-oxygen detonation structure and stability at low densities were performed under the assumption that the "velocity" of a detonation wave derives from complete burning of carbon and oxygen to iron. In fact, at densities <10^6 g/cm3 the detonation in SNIa may release less than a half of the available nuclear energy. In this paper we study basic properties of such detonations. We find that the length of an unsupported steady-state C-detonation is 30-100 times greater than previously estimated, and that the decreased energy has a drastic effect on the detonation stability. In contrast to complete detonations which are one-dimensionally stable, C-detonations may be one-dimensionally unstable and propagate by periodically re-igniting themselves via spontaneous burning. The re-ignition period at rho < 10^6 g/cm3 is estimated to be greater than the time-scale of a SNIa explosion. This suggests that propagation and quenching of C-detonations at these densities could be affected by the instability. Potential observational implications of this effect are discussed.

High-quality observations of $B$ and $V$ light curves obtained at Las Campanas Observatory for local Type Ia Supernovae (SNe Ia) show clear evidence that SNe Ia with the same brightness decline or stretch may have systematic and independent deviations at times < 5 days before and at times > 30 days after maximum light. This suggests the existence of two independent secondary parameters which control the shape of SN Ia light curves in addition to the brightness decline relation. stretch. The differences are consistent in morphology of the time dependence and size with predictions by models within the delayed detonation scenario. The secondary parameters may reflect two independent physical effects caused by variations in the progenitor and accretion rates, and link the LC variations in shape with the intrinsic, absolute brightness.

We present an analysis of well-posedness of constrained evolution of 3+1 formulations of GR. In this analysis we explicitly take into account the energy and momentum constraints as well as possible algebraic constraints on the evolution of high-frequency perturbations of solutions of Einstein's equations. In this respect, our approach is principally different from standard analyses of well-posedness of free evolution in general relativity. Our study reveals the existence of subsets of the linearized Einstein's equations that control the well-posedness of constrained evolution. It is demonstrated that the well-posedness of ADM, BSSN and other 3+1 formulations derived from ADM by adding combinations of constraints to the right-hand-side of ADM and/or by linear transformation of the dynamical ADM variables depends entirely on the properties of the gauge. For certain classes of gauges we formulate conditions for well-posedness of constrained evolution. This provides a new basis for constructing stable numerical integration schemes for a classical Arnowitt--Deser--Misner (ADM) and many other 3+1 formulations of general relativity.

We present HST images of the remnant of SN 1885 seen in absorption against M31's bulge via resonance lines of Ca I, Ca II, Fe I, and Fe II. Viewed in CaII H & K line absorption, the remnant appears as a nearly black circular spot with an outermost angular radius of 0.40" +/- 0.025" implying r = 1.52 pc and a 120 yr average expansion velocity of 12400 +/-1400 km/s. The strongest Ca II absorption is organized in a broken ring structure with a radius of 0.20" (=6000 km/s) with several apparent absorption `clumps' of an angular size near the pixel scale of 0.05" (= 1500 km/s). The detection of Ca II clumps is the first direct evidence for some instabilities and the existence of a deflagration phase in SNe Ia or, alternatively, mixing induced by radioactive decay of 56^Ni over time scales of seconds or days. However, the degree of mixing allowed by the observed images is much smaller than current 3D calculations for Rayleigh-Taylor dominated deflagration fronts. The images also require a central region of no or little Ca but iron group elements indicative of burning under sufficiently high densities for electron capture taking place, i.e., burning prior to a significant pre-expansion of the WD.

The thermonuclear explosion of a massive white dwarf in a Type Ia supernova explosion is characterized by vastly disparate spatial and temporal scales. The extreme dynamic range inherent to the problem prevents the use of direct numerical simulation and forces modelers to resort to subgrid models to describe physical processes taking place on unresolved scales. We consider the evolution of a model thermonuclear flame in a constant gravitational field on a periodic domain. The gravitational acceleration is aligned with the overall direction of the flame propagation, making the flame surface subject to the Rayleigh-Taylor instability. The flame evolution is followed through an extended initial transient phase well into the steady-state regime. The properties of the evolution of flame surface are examined. We confirm the form of the governing equation of the evolution suggested by Khokhlov (1995). The mechanism of vorticity production and the interaction between vortices and the flame surface are discussed. The results of our investigation provide the bases for revising and extending previous subgrid-scale model.

We present the method for computation of fluid flows that are characterized by the large degree of expansion/contraction and in which the fluid velocity is dominated by the bulk component associated with the expansion/contraction and/or rotation of the flow. We consider the formulation of Euler equations of fluid dynamics in a homologously expanding/contracting and/or rotating reference frame. The frame motion is adjusted to minimize local fluid velocities. Such approach allows to accommodate very efficiently large degrees of change in the flow extent. Moreover, by excluding the contribution of the bulk flow to the total energy the method eliminates the high Mach number problem in the flows of interest. An important practical advantage of the method is that it can be easily implemented with virtually any implicit or explicit Eulerian hydrodynamic scheme and adaptive mesh refinement (AMR) strategy. We also consider in detail equation invariance and existence of conservative formulation of equations for special classes of expanding/contracting reference frames. Special emphasis is placed on extensive numerical testing of the method for a variety of reference frame motions, which are representative of the realistic applications of the method. We study accuracy, conservativity, and convergence properties of the method both in problems which are not its optimal applications as well as in systems in which the use of this method is maximally beneficial. Such detailed investigation of the numerical solution behavior is used to define the requirements that need to be considered in devising problem-specific fluid motion feedback mechanisms.

We report observations of SN 2004dt obtained with the Very Large Telescope of the European Southern Observatory on August 13.30, 2004 when the supernova was more than a week before optical maximum. SN 2004dt showed strong lines of \ion{O}{1}, \ion{Mg}{2}, \ion{Si}{2}, and \ion{Ca}{2} with typical velocities of absorption minimum around 17,000 \kms. The line profiles show material moving at velocities as high as 25,000 \kms in these lines. The observations also reveal absorption lines from \ion{S}{2} and \ion{Si}{3} with a velocity of only 11,000 \kms. The highest velocity in the \ion{S}{2} features can be traced no higher than 15,000 \kms, much lower than those of O, Mg, Si, and Ca. SN 2004dt has a polarization spectrum unlike any previously observed. The variation of the polarization across some \ion{Si}{2} lines approaches 2%, making SN 2004dt the most highly polarized SN Ia ever observed. In contrast, the strong line of O I at 777.4 nm shows little or no polarization signature. The degree of polarization points to a richly-structured partially burned silicon layer with substantial departure from spherical symmetry. A geometry that would account for the observations is one in which the distribution of oxygen is essentially spherically symmetric, but with bubbles of intermediate-mass elements with significant opacity within the oxygen substrate.

We study a Type Ia supernova explosion using large-scale three-dimensional numerical simulations based on reactive fluid dynamics with a simplified mechanism for nuclear reactions and energy release. The initial deflagration stage of the explosion involves a subsonic turbulent thermonuclear flame propagating in the gravitational field of an expanding white dwarf. The deflagration produces an inhomogeneous mixture of unburned carbon and oxygen with intermediate-mass and iron-group elements in central parts of the star. During the subsequent detonation stage, a supersonic detonation wave propagates through the material unburned by the deflagration. The total energy released in this delayed-detonation process, (1.3-1.6)x10^51 ergs, is consistent with a typical range of kinetic energies obtained from observations. In contrast to the deflagration model that releases only about 0.6x10^51 ergs, the delayed-detonation model does not leave carbon, oxygen, and intermediate-mass elements in central parts of a white dwarf. This removes the key disagreement between three-dimensional simulations and observations, and makes a delayed detonation the mostly likely mechanism for Type Ia supernova explosions.

We study a type Ia supernova explosion using three-dimensional numerical simulations based on reactive fluid dynamics. We consider a delayed-detonation model that assumes a deflagration-to-detonation transition. In contrast to the pure deflagration model, the delayed-detonation model releases enough energy to account for a healthy explosion, and does not leave carbon, oxygen, and intermediate-mass elements in central parts of a white dwarf. This removes the key disagreement between simulations and observations, and makes a delayed detonation the mostly likely mechanism for type Ia supernovae.

High-quality spectropolarimetry (range 417-860 nm; spectral resolution 1.27 nm and 0.265 nm/pixel) of the SN Ia 2001el were obtained with the ESO Very Large Telescope Melipal (+ FORS1) at 5 epochs. The spectra a week before maximum and around maximum indicate photospheric expansion velocities of about 10,000 km s$^{-1}$. Prior to optical maximum, the linear polarization of the continuum was $\approx 0.2 - 0.3 %$ with a constant position angle, showing that SN 2001el has a well-defined axis of symmetry. The polarization was nearly undetectable a week after optical maximum. The spectra are similar to those of the normally-bright SN 1994D with the exception of a strong double-troughed absorption feature seen around 800 nm (FWHM about 22 nm). This feature could be an important clue to the binary nature of SN Ia, perhaps associated with an accretion disk, or to the nature of the thermonuclear burning. If modeled in terms of an oblate spheroid, the continuum polarization implies a minor to major axis ratio of around 0.9 if seen equator-on; this level of asymmetry would produce an absolute luminosity dispersion of about 0.1 mag when viewed at different viewing angles. If typical for SNe Ia, this would create an RMS scatter of several hundredths of a magnitude around the mean brightness-decline relation.

Large-scale three-dimensional numerical simulations of the deflagration stage of a thermonuclear supernova explosion show the formation and evolution of a highly convoluted turbulent flame in a gravitational field of an expanding carbon-oxygen white dwarf. The flame dynamics is dominated by the gravity-induced Rayleigh-Taylor instability that controls the burning rate. The thermonuclear deflagration releases enough energy to produce a healthy explosion. The turbulent flame, however, leaves large amounts of unburnt and partially burnt material near the star center, whereas observations imply these materials only in outer layers. This disagreement could be resolved if the deflagration triggers a detonation.

Core collapse supernovae(SN) are the final stages of evolution in massive stars during which the central region collapses. Recent explosion scenarios assumed that the ejection is due to energy deposition by neutrinos into the envelope but detailed models do not produce powerful explosions. There is new and mounting evidence for an asphericity and, in particular, for axial symmetry in several SNe which may be hard to reconcile within the spherical picture. This evidence includes the observed high polarization and its variation with time, pulsar kicks, high velocity iron-group and intermediate-mass elements material observed in remnants, direct observations of the debris of SN1987A etc. Some of the new evidence is discussed in more detail. To be in agreement with the observations, any successful mechanism must invoke some sort of axial symmetry for the explosion. We consider jet-induced/dominated explosions of core collapse SNe. Our study is based on detailed 3-D hydrodynamical and radiation transport models We find that the observations can be explained by low velocity, massive jets which stall well within the SN envelope. Such outflows may be produced by MHD- mechanisms, convective dominated accretion disks on the central object or asymmetric neutrino emissions. Asymmetric density/chemical distributions and, for SN2002ap, off-center energy depositions have been identified as crucial for the interpretation of the polarization.

Extensive early observations proved that the ejecta of supernova 1987A (SN 1987A) are aspherical. Fifteen years after the supernova explosion, the Hubble Space Telescope has resolved the rapidly expanding ejecta. The late-time images and spectroscopy provide a geometrical picture that is consistent with early observations and suggests a highly structured, axially symmetric geometry. We present here a new synthesis of the old and new data. We show that the Bochum event, presumably a clump of $^{56}$Ni, and the late-time image, the locus of excitation by $^{44}$Ti, are most naturally accounted for by sharing a common position angle of about 14\degree, the same as the mystery spot and early speckle data on the ejecta, and that they are both oriented along the axis of the inner circumstellar ring at 45\degree to the plane of the sky. We also demonstrate that the polarization represents a prolate geometry with the same position angle and axis as the early speckle data and the late-time image and hence that the geometry has been fixed in time and throughout the ejecta. The Bochum event and the Doppler kinematics of the [Ca II]/[O II] emission in spatially resolved HST spectra of the ejecta can be consistently integrated into this geometry. The radioactive clump is deduced to fall approximately along the axis of the inner circumstellar ring and therefore to be redshifted in the North whereas the [Ca II]/[O II] 7300 \AA emission is redshifted in the South. We present a jet-induced model for the explosion and argue that such a model can account for many of the observed asymmetries. In the jet models, the oxygen and calcium are not expected to be distributed along the jet, but primarily in an expanding torus that shares the plane and northern blue shift of the inner circumstellar ring.

We discuss the possible connection between supernova explosions (SN) and gamma-ray bursters (GRB) from the perspective of our current understanding of SN physics. Recent evidence strongly suggests that the explosion mechanism of core collapse SN is intrinsically aspherical. Typically, a neutron star is formed. However, the observed properties of the expanding SN envelopes remnants make these objects very unlikely candidates for GRBs. Most candidates for a GRB/SN connection seem to require the prompt or delayed formation of a black hole. These include the collapse of very massive stars (e.g. hypernovae) and 'classical' SNe with a significant fallback of material over time scales of hours to days, resulting in the collapse of the neutron star to a black hole. We suggest the merger of a neutron star with a white dwarf as a subclass of thermonuclear SNe and a potential candidate for a SN/GRB connection.

Core collapse supernovae (SN) are the final stages of stellar evolution in massive stars during which the central region collapses, forms a neutron star (NS), and the outer layers are ejected. Recent explosion scenarios assumed that the ejection is due to energy deposition by neutrinos into the envelope but detailed models do not produce powerful explosions. There is new and mounting evidence for an asphericity and, in particular, for axial symmetry in several SNe which may be hard to reconcile within the spherical picture. This evidence includes the observed high polarization and its variation with time, pulsar kicks, high velocity iron-group and intermediate-mass elements material in remnants & direct observations of the debris of SN87A.Any successful mechanism must invoke some sort of axial symmetry for the explosion. As limiting case, we consider jet-induced/dominated explosions of "classical" core collapse SNe. Our study is based on detailed 3-D hydro and radiation transport models. The jet properties and of the progenitor structure influence on the final density and chemical structure. Our calculations show that low velocity, massive jets can explain the observations. Both asymmetric ionization and density/chemical distributions have been identified as crucial for the formation of asymmetric photospheres. Even within the picture of jet-induced explosion, the latter effect alone fails to explain early polarization in core collapse SNe with a massive, hydrogen-rich envelopes such as SN1999em.

Core collapse supernovae (SN) are the final stages of stellar evolution in massive stars during which the central region collapses, forms a neutron star (NS), and the outer layers are ejected. Recent explosion scenarios assumed that the ejection is due to energy deposition by neutrinos into the envelope but detailed models do not produce powerful explosions. There is mounting evidence for an asphericity in the SN which is difficult to explain within this picture. This evidence includes the observed high polarization, pulsar kicks, high velocity iron-group and intermediate-mass elements material observed in remnants, etc. The discovery of highly magnetars revived the idea that the basic mechanism for the ejection of the envelope is related to a highly focused MHD-jet formed at the NS. Our 3-D hydro simulations of the jet propagation through the star confirmed that the mechanism can explain the asphericities. In this paper, detailed 3-D models for jet induced explosions of "classical" core collapse supernovae are presented. We demonstrate the influence of the jet properties and of the underlaying progenitor structure on the final density and chemical structure. Finally, we discuss the observational consequences, predictions and tests of this scenario.

The paper describes a physical model and numerical algorithm for modeling Type Ia supernova (SNIa) explosions in three dimensions and presents first results of modeling a deflagration explosion in a nonrotating, Chandrasekhar-mass carbon-oxygen (CO) white dwarf. Simulations show that the turbulent flame speed grows exponentially, reaches approximately 30% of the speed of sound, and then declines as the large-scale turbulence is frozen by expansion. The freezing of turbulent motions appears to be a crucial physical mechanism regulating the rate of deflagration in SNIa. The energy of the explosion is comparable to that of a typical SNIa. However, the presence of the outer layer of unburned CO and the formation of intermediate mass elements and pockets of unburned CO near the center pose problems for the modeling of SNIa spectra. Delayed detonation is a way to alleviate these problems and to produce consistent spectra.

Thermonuclear explosions of Type Ia supernovae (SNIa) involve turbulent deflagrations, detonations, and possibly a deflagration-to-detonation transition. A phenomenological delayed detonation model of SNIa successfully explains many observational properties of SNIa including monochromatic light curves, spectra, brightness - decline and color - decline relations. Observed variations among SNia are explained as a result of varying nickel mass synthesised in an explosion of a Chandrasekhar mass C/O white dwarf. Based on theoretical models of SNIa, the value of the Hubble constant H_o \simeq 67km/s/Mpc was determined without the use of secondary distance indicators. The cause for the nickel mass variations in SNIa is still debated. It may be a variation of the initial C/O ratio in a supernova progenitor, rotation, or other effects.

We have simulated the formation of an X-ray cluster in a cold dark matter universe using 12 different codes. The codes span the range of numerical techniques and implementations currently in use, including SPH and grid methods with fixed, deformable or multilevel meshes. The goal of this comparison is to assess the reliability of cosmological gas dynamical simulations of clusters in the simplest astrophysically relevant case, that in which the gas is assumed to be non-radiative. We compare images of the cluster at different epochs, global properties such as mass, temperature and X-ray luminosity, and radial profiles of various dynamical and thermodynamical quantities. On the whole, the agreement among the various simulations is gratifying although a number of discrepancies exist. Agreement is best for properties of the dark matter and worst for the total X-ray luminosity. Even in this case, simulations that adequately resolve the core radius of the gas distribution predict total X-ray luminosities that agree to within a factor of two. Other quantities are reproduced to much higher accuracy. For example, the temperature and gas mass fraction within the virial radius agree to about 10%, and the ratio of specific kinetic to thermal energies of the gas agree to about 5%. Various factors contribute to the spread in calculated cluster properties, including differences in the internal timing of the simulations. Based on the overall consistency of results, we discuss a number of general properties of the cluster we have modelled.

We numerically studied the explosion of a supernova caused by supersonic jets present in its center. The jets are assumed to be generated by a magneto-rotational mechanism when a stellar core collapses into a neutron star. We simulated the process of the jet propagation through the star, jet breakthrough, and the ejection of the supernova envelope by the lateral shocks generated during jet propagation. The end result of the interaction is a highly nonspherical supernova explosion with two high-velocity jets of material moving in polar directions, and a slower moving, oblate, highly distorted ejecta containing most of the supernova material. The jet-induced explosion is entirely due to the action of the jets on the surrounding star and does not depend on neutrino transport or re-acceleration of a stalled shock. The jet mechanism can explain the observed high polarization of Type Ib,c and Type II supernovae, pulsar kicks, very high velocity material observed in supernova remnants, indications that radioactive material was carried to the hydrogen-rich layers in SN1987A, and some others observations that are very difficult or impossible to explain by the neutrino energy deposition mechanism. The breakout of the jet from a compact, hydrogen- deficient core may account for the gamma-ray bursts and radio outburst associated with SN1998bw/GRB980425.

We present analysis of the evolution of dark matter halos in dense environments of groups and clusters in dissipationless cosmological simulations. The premature destruction of halos in such environments, known as the overmerging, reduces the predictive power of N-body simulations and makes difficult any comparison between models and observations. We analyze the possible processes that cause the overmerging and assess the extent to which this problem can be cured with current computer resources and codes. Using both analytic estimates and high resolution numerical simulations, we argue that the overmerging is mainly due to the lack of numerical resolution. We find that the force and mass resolution required for a simulated halo to survive in galaxy groups and clusters is extremely high and was almost never reached before: ~1-3 kpc and 10^8-10^9 Msun, respectively. We use the high-resolution Adaptive Refinement Tree (ART) N-body code to run cosmological simulations with the particle mass of \approx 2x10^8/h Msun} and the spatial resolution of \approx 1-2/h kpc, and show that in these simulations the halos do survive in regions that would appear overmerged with lower force resolution. Nevertheless, the halo identification in very dense environments remains a challenge even with the resolution this high. We present two new halo finding algorithms developed to identify both isolated and satellite halos that are stable (existed at previous moments) and gravitationally bound. To illustrate the use of the satellite halos that survive the overmerging, we present a series of halo statistics, that can be compared with those of observed galaxies. (Abridged)

We study the evolution of the halo-halo correlation function and bias in a LCDM model using very high-resolution N-body simulations. The high force and mass resolution allows dark matter (DM) halos to survive in the tidal fields in high-density regions and thus prevents the ambiguities related with the ``overmerging problem.'' Numbers of galaxy-size halos in cluster-like objects in our simulation are similar to the numbers of galaxies observed in real clusters. This allows us to estimate for the first time the evolution of the correlation function and bias at small (down to ~100/h kpc) scales. We compare particle distribution, dark matter correlation function, density profiles, and halo mass function produced with our N-body code and corresponding results of the AP3M simulations. We find that at all epochs the 2-point correlation function of galaxy-size halos xihh is well approximated by a power-law with slope ~1.6-1.8. The difference between the shape of xihh and the shape of the correlation function of matter results in the scale-dependent bias at scales <7/h Mpc, which we find to be a generic prediction of the hierarchical models, independent of the epoch and of the model details. We find that our results agree well with existing clustering data at different redshifts, indicating the general success of the picture of structure formation in which galaxies form inside the host DM halos. Particularly, we find an excellent agreement in both slope and the amplitude between xihh(z=0) in our simulation and the galaxy correlation function measured using the APM galaxy survey. At high redshifts, the observed clustering of the Lyman-break galaxies is also reasonably well reproduced by the models.

We study the evolution of the halo-halo correlation function and bias in four cosmological models (LCDM, OCDM, tauCDM, and SCDM) using very high-resolution N-body simulations. The high force and mass resolution allows dark matter (DM) halos to survive in the tidal fields of high-density regions and thus prevents the ambiguities related with the ``overmerging problem.'' This allows us to estimate for the first time the evolution of the correlation function and bias at small (down to ~100/h kpc) scales. We find that at all epochs the 2-point correlation function of galaxy-size halos xi_hh is well approximated by a power-law with slope ~1.6-1.8. The difference between the shape of xi_hh and the shape of the correlation function of matter results in the scale-dependent bias at scales <7/h Mpc, which we find to be a generic prediction of the hierarchical models. The bias evolves rapidly from a high value of ~2-5 at z~3-7 to the anti-bias of b~0.5-1 at small <5/h Mpc scales at z=0. We find that our results agree well with existing clustering data at different redshifts. Particularly, we find an excellent agreement in both slope and the amplitude between xi_hh(z=0) in our LCDM simulation and the galaxy correlation function measured using the APM galaxy survey. At high redshifts, the observed clustering of the Lyman-break galaxies is also well reproduced by the models. The agreement with the data at high and low z indicates the general success of the hierarchical models of structure formation in which galaxies form inside the host DM halos. (Abridged)

Examination of the process of silicon burning, the burning stage that leads to the production of the iron peak nuclei, reveals that the nuclear evolution is dominated by large groups of nuclei in mutual equilibrium. These quasi-equilibrium (QSE) groups form well in advance of the global Nuclear Statistical Equilibrium (NSE). We present an improved ``minimal'' nuclear network, which takes advantage of quasi-equilibrium in order to calculate the nuclear evolution and energy generation while further reducing the computational cost compared to a conventional \alpha-chain network. During silicon burning, the resultant \emph{QSE-reduced} \alpha network is twice as fast as the full \alpha network it replaces and requires the tracking of only half as many abundance variables, without significant loss of accuracy. When the QSE-reduced \alpha network is used in combination with a conventional \alpha network stretching from He to Si, the combined \alpha7 network provides an accurate approximation for all of the burning stages from He burning to NSE, while tracking only 7 abundances. These reductions in computational cost and the number of species evolved make the \alpha7 network well suited for inclusion within hydrodynamic simulations, particularly those in multi-dimension.

The gamma-ray light curves and spectra are presented for a set of theoretical Type Ia supernova models including deflagration, detonation, delayed detonation, and pulsating delayed detonations of Chandrasekhar mass white dwarfs as well as merger scenarios that may involve more than the Chandrasekhar mass and helium detonations of sub-Chandrasekhar mass white dwarfs. The results have been obtained with a Monte Carlo radiation transport scheme which takes into account all relevant gamma-transitions and interaction processes. The result is a set of accurate line profiles which are characteristic of the initial Ni-mass distribution of the supernova models. The gamma-rays probe the isotopic rather than just the elemental distribution of the radioactive elements in the ejecta. Details of the line profiles including the line width, shift with respect to the rest frame, and line ratios are discussed. With sufficient energy and temporal resolution, different model scenarios can clearly be distinguished. Observational strategies are discussed for current and immediately upcoming generations of satellites (CGRO and INTEGRAL) as well as projected future missions including concepts such as Laue telescopes. With CGRO, it is currently possible with sufficiently early observations (near optical maximum) to distinguish helium detonations from explosions of Chandrasekhar mass progenitors and of those involving mergers up to a distance of about 15 Mpc. This translates into one target of opportunity every eight years. SNe Ia up to about 10 Mpc would allow detailed CGRO studies of line ratios of Co lines.

We present a new high-resolution N-body algorithm for cosmological simulations. The algorithm employs a traditional particle-mesh technique on a cubic grid and successive multilevel relaxations on the finer meshes, introduced recursively in a fully adaptive manner in the regions where the density exceeds a predefined threshold. The mesh is generated to effectively match an arbitrary geometry of the underlying density field -- a property particularly important for cosmological simulations. In a simulation the mesh structure is not created at every time step but is properly adjusted to the evolving particle distribution. The algorithm is fast and effectively parallel. We present a detailed description of the methodology, implementation, and tests of the code. We further use the code to study the structure of dark matter halos in high-resolution (2/h kpc) simulations of standard CDM (Omega=1, h=0.5, sigma_8=0.63) and LCDM (Omega_Lambda=1-Omega_0=0.7, h=0.7, sigma_8=1.0) models. We find that halo density profiles in both CDM and LCDM models are well fitted by the analytical model presented recently by Navarro et al., which predicts a singular [$\rho(r)\propto r^{-1}$] behavior of the halo density profiles at small radii. We therefore conclude that halos formed in the $\Lambda$CDM model have structure similar to CDM halos and thus cannot explain the dynamics of the central parts of dwarf spiral galaxies, as inferred from the galaxies' rotation curves.

A fully threaded tree (FTT) for adaptive refinement of regular meshes is described. By using a tree threaded at all levels, tree traversals for finding nearest neighbors are avoided. All operations on a tree including tree modifications are O(N), where N is a number of cells, and are performed in parallel. An efficient implementation of the tree is described that requires 2N words of memory. A filtering algorithm for removing high-frequency noise during mesh refinement is described. A FTT can be used in various numerical applications. In this paper, it is applied to the integration of the Euler equations of fluid dynamics. An adaptive-mesh time stepping algorithm is described in which different time steps are used at different levels of the tree. Time stepping and mesh refinement are interleaved to avoid extensive buffer layers of fine mesh which were otherwise required ahead of moving shocks. Test examples are presented, and the FTT performance is evaluated. The three-dimensional simulation of the interaction of a shock wave and a spherical bubble is carried out that shows the development of azimuthal perturbations on the bubble surface.

We derive the criteria for deflagration to detonation transition (DDT) in a Type Ia supernova. The theory is based on the two major assumptions: (i) detonation is triggered via the Zeldovich gradient mechanism inside a region of mixed fuel and products, (ii) the mixed region is produced by a turbulent mixing of fuel and products either inside an active deflagration front or during the global expansion and subsequent contraction of an exploding white dwarf. We determine the critical size of the mixed region required to initiate a detonation in a degenerate carbon-oxygen mixture. This critical length is much larger than the width of the reaction front of a Chapman-Jouguet detonation. However, at densities greater than simeq 5 x 10^6 g cm^-3, it is much smaller than the size of a white dwarf. We derive the critical turbulent intensity required to create the mixed region inside an active deflagration front in which a detonation can form. We conclude that the density rho_tr at which a detonation can form in a carbon-oxygem white dwarf is low, less than 2 - 5 x 10^7 g cm^-3, but greater than 5 x 10^6 g cm^-3.

We compare the observed correlations between the maximum brightness, postmaximum decline rate and color at maximum light of Type Ia supernovae (SN Ia) with model predictions. The observations are based on a total of 40 SN Ia with 29 SN of the Calan Tololo Supernova Search and 11 local SN which cover a range of 2 mag in the absolute visual brightness. The observed correlations are not tight, one dimensional relations. Supernovae with the same postmaximum decline or the same color have a spread in visual magnitude of about 0.7 mag. The dispersion in the color-magnitude relation may result from uncertainties in the distance determinations or the interstellar reddening within the host galaxy. The dispersion in the decline rate-magnitude relation suggests that an intrinsic spread in the supernova properties exists that cannot be accounted for by any single relation between visual brightness and postmaximum decline. Theoretical correlations are derived from a grid of models which encompasses delayed detonations, pulsating delayed detonations, the merging scenario and helium detonations. We find that the observed correlations can be understood in terms of explosions of Chandrasekhar mass white dwarfs. Our models show an intrinsic spread in the relations of about 0.5 mag in the maximum brightness and about 0.1 mag in the B-V color. Our study provides strong evidence against the mechanism of helium detonation for subluminous, red SN Ia.

We numerically simulate some of the most critical physical processes in galaxy formation: The supernova feedback, in conjunction with gasdynamics and gravity, plays a crucial role in determining how galaxies arise within the context of a model for large-scale structure. Our treatment incorporates a multi-phase model of the interstellar medium and includes the effects of cooling, heating and metal enrichment by supernovae, and evaporation of cold clouds. The star formation happens inside the clouds of cold gas, which are produced via thermal instability. We simulate the galaxy formation in standard biased CDM model for a variety of parameters and for several resolutions in the range 2--20$h^{-1}$kpc. In our picture, supernova feedback regulates the evolution of the gas components and star formation. The efficiency of cloud evaporation by supernova strongly influences star formation rates. This feedback results in a steady rate of star formation in large galaxies (mass larger than $2-3x10^{11}\Msun$) at a level of $(1-10)\Msun\yr$ for $z<3$. Supernova feedback has an even stronger effect on the evolution of dwarf galaxies, most of which have a small fraction of stars and extremely low luminosities: $M_R>-15$. In the case of both large and small galaxies, the distribution of luminous matter (stars) is strongly BIASED with respect to the dark matter. We find an approximate biasing measure of the form $\rho_{lum}= (\rho_{dm}/133)^{1.7}$ for z=0 and overdensities exceeding 1000. Deviations from this relation (a factor 2-3) depend on the environment. For halo masses exceeding $2x10^{10}\Msun$, the dependence of the absolute magnitude on the total mass can be approximated as $M_V=-18.5-4\log(M_{tot}/10^{11}\Msun)$, with a scatter of less than 0.5mag.

This paper outlines a theoretical approach for predicting the onset of detonation in unconfined turbulent flames which is relevant both to problems of terrestrial combustion and to thermonuclear burning in Type Ia supernovae. Two basic assumuptions are made: 1) the gradient mechanism is the inherent mechanism that leads to DDT in unconfined conditions, and 2) the sole mechanism for preparing the gradient in induction time is by turbulent mixing and local flame quenching. The criterion for DDT is derived in terms of the one-dimensional detonation wave thickness, the laminar flame speed, and the laminar flame thickness in the reactive gas. This approach gives a lower-bound criterion for DDT for conditions where shock preheating, wall effects, and interactions with obstacles are absent. Regions in parameter space where unconfined DDT can and cannot occur are determined. A subsequent paper will address these issues specifically in the astrophysical context.

From the spectra and light curves it is clear that SNIa are thermonuclear explosions of white dwarfs. However, details of the explosion are highly under debate. Here, we present detailed models which are consistent with respect to the explosion mechanism, the optical and infrared light curves (LC), and the spectral evolution. This leaves the description of the burning front and the structure of the white dwarf as the only free parameters. The explosions are calculated using one-dimensional Lagrangian codes including nuclear networks. Subsequently, optical and IR-LCs are constructed. Detailed NLTE-spectra are computed for several instants of time using the density, chemical and luminosity structure resulting from the LCs. The general methods and critical tests are presented (sect. 2). Different models for the thermonuclear explosion are discussed including detonations, deflagrations, delayed detonations, pulsating delayed detonations (PDD) and helium detonations (sect.3). Comparisons between theoretical and observed LCs and spectra provide an insight into details of the explosion and nature of the progenitor stars (sect. 4/5). We try to answer several related questions. Are subluminous SNe Ia a group different from `normal' SN Ia (sect. 5)? Can we understand observed properties of the LCs and spectra (sect. 4)? What do we learn about the progenitor evolution and its metallicity (sect. 3, Figs. 4,5)? Do successful SN~Ia models depend on the type of the host galaxy (Table 2)? Using both the spectral and LC information, theoretical models allow for a determination of the Hubble constant independent from `local' distance indicators such as delta-Cephei stars. Ho is found to be 67 +- 9km/s/Mpc and, from SN1988U, qo equals 0.7 +- 1. within 95 percent confidence levels.

Turbulent combustion is three-dimensional. Turbulence in a Type Ia supernova is driven on large scales by the buoyancy of burning products. The turbulent cascade penetrates down to very small scales, and makes the rate of deflagration independent of the microphysics. The competition between the turbulent cascade and the freeze-out of turbulent motions due to stellar expansion determines the largest scale participating in the cascade. This sets the bulk rate of a deflagration in a supernova. The freeze-out limits the bulk rate of deflagration to a value that makes a powerful explosion impossible. Two-dimensional simulations cannot capture these essential elements of turbulent combustion, even in principle. A powerful delayed detonation explosion can take place if the burning makes a transition to a detonation. A deflagration to detonation transition (DDT) can occur in a layer of mixed cold fuel and hot burning products created either inside an active turbulent burning region by a high intensity of turbulence, or by mixing cold fuel with ashes of a dead deflagration front during the global pulsation of a star, or by both.