Surface detonation in type Ia supernova explosions?

We explore the evolution of thermonuclear supernova explosions when the progenitor white dwarf star ignites asymmetrically off-center. Several numerical simulations are carried out in two and three dimensions to test the consequences of different initial flame configurations such as spherical bubbles displaced from the center, more complex deformed configurations, and teardrop-shaped ignitions. The burning bubbles float towards the surface while releasing energy due to the nuclear reactions. If the energy release is too small to gravitationally unbind the star, the ash sweeps around it, once the burning bubble approaches the surface. Collisions in the fuel on the opposite side increase its temperature and density and may -- in some cases -- initiate a detonation wave which will then propagate inward burning the core of the star and leading to a strong explosion. However, for initial setups in two dimensions that seem realistic from pre-ignition evolution, as well as for all three-dimensional simulations the collimation of the surface material is found to be too weak to trigger a detonation.Comment: 5 pages, 3 figures, in: Proceedings of the SciDAC 2006 Meeting, Denver June 25-26 2006, also available at http://herald.iop.org/jpcs46/m51/gbr//link/40

Band structure approach to the resonant x-ray scattering

We study the resonance behaviour of the forbidden 600 and 222 x-ray Bragg peaks in Ge using LDA band structure methods. These Bragg peaks remain forbidden in the resonant dipole scattering approximation even taking into account the non local nature of the band states. However they become allowed at resonance if the eigenstates of the unoccupied conduction band involve a hybridization of p like and d like atomic states. We show that the energy dependence of the resonant behaviour, including the phase of the scattering, is a direct measure of this p-d hybridization.and obtain quantitative agreement with experiment. A simple physical picture involving a product of dipole and quadrupolar transition matrix elements explains this behaviour and shows that it should be generally true for cases where the resonating atom is not at an inversion center. This has strong implications for the description of the resonance behavior of x-ray scattering in materials where the resonant atom is not at an inversion center such as V2O3 and in ferro and antiferro electric and piezo electric materials in general.Comment: 4 pages, 5figure

The Possibility of Emersion of the Outer Layers in a Massive Star Simultaneously with Iron-Core Collapse: A Hydrodynamic Model

We analyze the behavior of the outer envelope in a massive star during and after the collapse of its iron core into a protoneutron star (PNS) in terms of the equations of one-dimensional spherically symmetric ideal hydrodynamics. The profiles obtained in the studies of the evolution of massive stars up to the final stages of their existence, immediately before a supernova explosion (Boyes et al. 1999), are used as the initial data for the distribution of thermodynamic quantities in the envelope.We use a complex equation of state for matter with allowances made for arbitrary electron degeneracy and relativity, the appearance of electron-positron pairs, the presence of radiation, and the possibility of iron nuclei dissociating into free nucleons and helium nuclei. We performed calculations with the help of a numerical scheme based on Godunov's method. These calculations allowed us to ascertain whether the emersion of the outer envelope in a massive star is possible through the following two mechanisms: first, the decrease in the gravitational mass of the central PNS through neutrino-signal emission and, second, the effect of hot nucleon bubbles, which are most likely formed in the PNS corona, on the envelope emersion. We show that the second mechanism is highly efficient in the range of acceptable masses of the nucleon bubbles ($\leq 0.01M_\odot$) simulated in our hydrodynamic calculations in a rough, spherically symmetric approximation.Comment: 14 pages, 11 figure

New Relativistic Effects in the Dynamics of Nonlinear Hydrodynamical Waves

In Newtonian and relativistic hydrodynamics the Riemann problem consists of calculating the evolution of a fluid which is initially characterized by two states having different values of uniform rest-mass density, pressure and velocity. When the fluid is allowed to relax, one of three possible wave-patterns is produced, corresponding to the propagation in opposite directions of two nonlinear hydrodynamical waves. New effects emerge in a special relativistic Riemann problem when velocities tangential to the initial discontinuity surface are present. We show that a smooth transition from one wave-pattern to another can be produced by varying the initial tangential velocities while otherwise maintaining the initial states unmodified. These special relativistic effects are produced by the coupling through the relativistic Lorentz factors and do not have a Newtonian counterpart.Comment: 4 pages, 5 figure

Mathisson-Papapetrou equations in metric and gauge theories of gravity in a Lagrangian formulation

We present a simple method to derive the semiclassical equations of motion for a spinning particle in a gravitational field. We investigate the cases of classical, rotating particles (pole-dipole particles), as well as particles with intrinsic spin. We show that, starting with a simple Lagrangian, one can derive equations for the spin evolution and momentum propagation in the framework of metric theories of gravity and in theories based on a Riemann-Cartan geometry (Poincare gauge theory), without explicitly referring to matter current densities (spin and energy-momentum). Our results agree with those derived from the multipole expansion of the current densities by the conventional Papapetrou method and from the WKB analysis for elementary particles.Comment: 28 page

Double-detonation supernovae of sub-Chandrasekhar mass white dwarfs

In the "double-detonation sub-Chandrasekhar" model for type Ia supernovae, a carbon-oxygen (C + O) white dwarf accumulates sufficient amounts of helium such that a detonation ignites in that layer before the Chandrasekhar mass is reached. This detonation is thought to trigger a secondary detonation in the C + O core. By means of one- and two-dimensional hydrodynamic simulations, we investigate the robustness of this explosion mechanism for generic 1-M_sun models and analyze its observable predictions. Also a resolution dependence in numerical simulations is analyzed. The propagation of thermonuclear detonation fronts, both in helium and in the carbon-oxygen mixture, is computed by means of both a level-set function and a simplified description for nuclear reactions. The decision whether a secondary detonation is triggered in the white dwarf's core or not is made based on criteria given in the literature. In a parameter study involving different initial flame geometries for He-shell masses of 0.2 and 0.1 M_sun, we find that a secondary detonation ignition is a very robust process. Converging shock waves originating from the detonation in the He shell generate the conditions for a detonation near the center of the white dwarf in most of the cases considered. Finally, we follow the complete evolution of three selected models with 0.2 M_sun of He through the C/O-detonation phase and obtain nickel-masses of about 0.40 to 0.45 M_sun. Although we have not done a complete scan of the possible parameter space, our results show that sub-Chandrasekhar models are not good candidates for normal or sub-luminous type Ia supernovae. The chemical composition of the ejecta features significant amounts of nickel in the outer layers at high expansion velocities, which is inconsistent with near-maximum spectra. (abbreviated)Comment: 11 pages, 10 figures, PDFLaTeX, accepted for publication in A&

Gravitation: Global Formulation and Quantum Effects

A nonintegrable phase-factor global approach to gravitation is developed by using the similarity of teleparallel gravity with electromagnetism. The phase shifts of both the COW and the gravitational Aharonov-Bohm effects are obtained. It is then shown, by considering a simple slit experiment, that in the classical limit the global approach yields the same result as the gravitational Lorentz force equation of teleparallel gravity. It represents, therefore, the quantum mechanical version of the classical description provided by the gravitational Lorentz force equation. As teleparallel gravity can be formulated independently of the equivalence principle, it will consequently require no generalization of this principle at the quantum level.Comment: Latex (IOP style), 14 pages, 3 figures. To appear in Classical and Quantum Gravit

Three-Dimensional Simulations of a Starburst-Driven Galactic Wind

We have performed a series of three-dimensional simulations of a starburst-driven wind in an inhomogeneous interstellar medium. The introduction of an inhomogeneous disk leads to differences in the formation of a wind, most noticeably the absence of the ``blow-out'' effect seen in homogeneous models. A wind forms from a series of small bubbles that propagate into the tenuous gas between dense clouds in the disk. These bubbles merge and follow the path of least resistance out of the disk, before flowing freely into the halo. Filaments are formed from disk gas that is broken up and accelerated into the outflow. These filaments are distributed throughout a biconical structure within a more spherically distributed hot wind. The distribution of the inhomogeneous interstellar medium in the disk is important in determining the morphology of this wind, as well as the distribution of the filaments. While higher resolution simulations are required in order to ascertain the importance of mixing processes, we find that soft X-ray emission arises from gas that has been mass-loaded from clouds in the disk, as well as from bow shocks upstream of clouds, driven into the flow by the ram pressure of the wind, and the interaction between these shocks.Comment: 37 pages, 16 figures, mpg movie can be obtained at http://www.mso.anu.edu.au/~jcooper/movie/video16.mpg, accepted for publication in Ap

Experimental feasibility of measuring the gravitational redshift of light using dispersion in optical fibers

This paper describes a new class of experiments that use dispersion in optical fibers to convert the gravitational frequency shift of light into a measurable phase shift or time delay. Two conceptual models are explored. In the first model, long counter-propagating pulses are used in a vertical fiber optic Sagnac interferometer. The second model uses optical solitons in vertically separated fiber optic storage rings. We discuss the feasibility of using such an instrument to make a high precision measurement of the gravitational frequency shift of light.Comment: 11 pages, 12 figure

Numerical Methods for the Stochastic Landau-Lifshitz Navier-Stokes Equations

The Landau-Lifshitz Navier-Stokes (LLNS) equations incorporate thermal fluctuations into macroscopic hydrodynamics by using stochastic fluxes. This paper examines explicit Eulerian discretizations of the full LLNS equations. Several CFD approaches are considered (including MacCormack's two-step Lax-Wendroff scheme and the Piecewise Parabolic Method) and are found to give good results (about 10% error) for the variances of momentum and energy fluctuations. However, neither of these schemes accurately reproduces the density fluctuations. We introduce a conservative centered scheme with a third-order Runge-Kutta temporal integrator that does accurately produce density fluctuations. A variety of numerical tests, including the random walk of a standing shock wave, are considered and results from the stochastic LLNS PDE solver are compared with theory, when available, and with molecular simulations using a Direct Simulation Monte Carlo (DSMC) algorithm