Search

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., Comment: 18 pages, 14 figures. Accepted for publication in ApJ

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., Comment: 15 pages, 22 postscript figures. Accepted for publication by the Astrophysical Journal. High resolution figures can be found at http://flash.uchicago.edu/~zhang/research_paper.html

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., Comment: Accepted to PRD, typos corrected, added references, extension of the analysis

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., Comment: Accepted to Journal of Computational Physics; several additional numerical tests included. High resolution fig. 5 and additional animations of simulations can be found at http://flash.uchicago.edu/~apol/JCP/

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., Comment: 19 pages, 5 figures. Submitted to Astrophysical Journal

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., Comment: Submitted to ApJ

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., Comment: 4 pages, 2 figures

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., Comment: ApJ, in press

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., Comment: 17 pages, 5 figures. To appear in Science, January 2003

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., Comment: 7 pages + 5 figures, IAU Symposium 212 on Massive Stars, D. Reidel Conf. Series, ed. E. van den Hucht

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., Comment: To Appear in ApJ

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., Comment: 9 pages + 6 figures (jpeg/separate), Jan van Paradijs Symposium, ed. Ed. van Heuvel et al., Kluver Press, Acknowledment added

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., Comment: 13 pages, 5 figures (BW=ps)/color=gif), latex, aipproc.sty. In 20th Texas Conference on Relativistic Astrophysics, AIP-Publ

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., Comment: 18 pages, 11 figures, 8 additional colour figures available from http://star-www.dur.ac.uk/~csf/SantaBarbara/ Accepted in Astrophysical Journal

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., Comment: Contributions to the Proceedings of the MPA/ESO Conference "Evolution of LSS: from Recombination to Garching", Garching, August 2-7, 1998. LaTeX (uses tp.sty, psfig.sty), 5 pages, 4 figs. Full-resolution version of Fig.2 is available at ftp://charon.nmsu.edu/pub/kravtsov/PAPERS/Garching98/

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), Comment: submitted to the Astrophys.Journal; 21 pages, LaTeX (uses emulateapj.sty); full resolution versions of figs.1 and 2 are available at http://astro.nmsu.edu/~akravtso/GROUP/group_publications.html or at ftp://charon.nmsu.edu/pub/kravtsov/PAPERS/Bias/

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., Comment: 26 pages, including 8 figures; LaTeX2e, uses aaspp4 and graphicx packages. To be published in Ap. J., Vol. 503

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., Comment: 32 pages, Tex, ApJ, in press

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), Comment: Accepted for publication in ApJ, substantional revisions after the first version, LaTeX 23 pages, 18 figs. (uses emulateapj.sty), Full-resolution version of Fig.9 is available upon request

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., Comment: revised in accord with published version, LaTeX (uses aaspp4.sty), 42 pages, 12 figures, all PostScript figures are available at http://astro.nmsu.edu/~akravtso/GROUP/group_publications.html or at ftp://charon.nmsu.edu/pub/kravtsov/PAPERS/ART

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., Comment: 7 pages, 3 figures, macros ''aaspp.sty'. LaTeX Style. Astrophysical Journal Letters, submitted Jul. 1995, revised Aug. 1995, resubmitted Sep. 1996

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., Comment: 17 pages LATEX (uses mn.sty). 12 PostScript figures and the text are available at ftp://charon.nmsu.edu/pub/aklypin/HYDRO Corrected ftp adderess

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., Comment: 20 pages, 14 figures, macros 'crckapb.sty'.For a postscript version with high resolution figures, send E-Mail to pah@cfa141.harvard.edu; To appear in ""Thermonuclear Supernovae" eds. R.Canal, P. Ruiz-Lapuente, J. Isern, (Dordrecht, Kluwer Academic Press)