A simple approach to the supernova progenitor–explosion connection

We present a new approach to understand the landscape of supernova explosion energies, ejected nickel masses, and neutron star birth masses. In contrast to other recent parametric approaches, our model predicts the properties of neutrino-driven explosions based on the pre-collapse stellar structure without the need for hydrodynamic simulations. The model is based on physically motivated scaling laws and simple differential equations describing the shock propagation, the contraction of the neutron star, the neutrino emission, the heating conditions, and the explosion energetics. Using model parameters compatible with multi-D simulations and a fine grid of thousands of supernova progenitors, we obtain a variegated landscape of neutron star and black hole formation similar to other parametrized approaches and find good agreement with semi-empirical measures for the ‘explodability’ of massive stars. Our predicted explosion properties largely conform to observed correlations between the nickel mass and explosion energy. Accounting for the coexistence of outflows and downflows during the explosion phase, we naturally obtain a positive correlation between explosion energy and ejecta mass. These correlations are relatively robust against parameter variations, but our results suggest that there is considerable leeway in parametric models to widen or narrow the mass ranges for black hole and neutron star formation and to scale explosion energies up or down. Our model is currently limited to an all-or-nothing treatment of fallback and there remain some minor discrepancies between model predictions and observational constraints.

Progenitor-Dependent Explosion Dynamics in Self-Consistent, Axisymmetric Simulations of Neutrino-Driven Core-Collapse Supernovae

We present self-consistent, axisymmetric core-collapse supernova simulations performed with the Prometheus-Vertex code for 18 pre-supernova models in the range of 11–28 M ⊙, including progenitors recently investigated by other groups. All models develop explosions, but depending on the progenitor structure, they can be divided into two classes. With a steep density decline at the Si/Si–O interface, the arrival of this interface at the shock front leads to a sudden drop of the mass-accretion rate, triggering a rapid approach to explosion. With a more gradually decreasing accretion rate, it takes longer for the neutrino heating to overcome the accretion ram pressure and explosions set in later. Early explosions are facilitated by high mass-accretion rates after bounce and correspondingly high neutrino luminosities combined with a pronounced drop of the accretion rate and ram pressure at the Si/Si–O interface. Because of rapidly shrinking neutron star radii and receding shock fronts after the passage through their maxima, our models exhibit short advection timescales, which favor the efficient growth of the standing accretion-shock instability. The latter plays a supportive role at least for the initiation of the re-expansion of the stalled shock before runaway. Taking into account the effects of turbulent pressure in the gain layer, we derive a generalized condition for the critical neutrino luminosity that captures the explosion behavior of all models very well. We validate the robustness of our findings by testing the influence of stochasticity, numerical resolution, and approximations in some aspects of the microphysics.

The role of nuclear physics in supernovae and the evolution of neutron stars, Neutrino Opacities, Equation of State, Transport Coefficients, and Dark Matter Production

Thesis (Ph.D.)--University of Washington, 2016-08

Outer crust of a cold non-accreting magnetar

The outer-crust structure and composition of a cold, non-accreting magnetar are studied. We model the outer crust to be made of fully equilibrated matter where ionized nuclei form a Coulomb crystal embedded in an electron gas. The main effects of the strong magnetic field are those of quantizing the electron motion in Landau levels and of modifying the nuclear single-particle levels producing, on average, an increased binding of nucleons in nuclei present in the Coulomb lattice. The effect of a homogeneous and constant magnetic field on nuclear masses has been predicted by using a covariant density functional in which induced currents and axial deformation due to the presence of a magnetic field that breaks time-reversal symmetry have been included self-consistently in the nucleon and meson equations of motion. Although not yet observed, for Ba 1016 G both effects contribute to produce different compositions - odd-mass nuclei are frequently predicted - and to increase the neutron-drip pressure as compared to a typical neutron star. Specifically, in such a regime, the magnetic-field effects on nuclei favor the appearance of heavier nuclei at low pressures. As B increases, such heavier nuclei are also preferred up to larger pressures. For the most extreme magnetic field considered, B=1018 G, and for the models studied, almost the whole outer crust is made of 4092Zr52.

On neutron stars in ƒ (R) theories: Small radii, large masses and large energy emitted in a merger

In the context of ƒ (R) gravity theories, we show that the apparent mass of a neutron star as seen from an observer at infinity is numerically calculable but requires careful matching, first at the star’s edge, between interior and exterior solutions, none of them being totally Schwarzschild-like but presenting instead small oscillations of the curvature scalar R; and second at large radii, where the Newtonian potential is used to identify the mass of the neutron star. We find that for the same equation of state, this mass definition is always larger than its general relativistic counterpart. We exemplify this with quadratic R^2 and Hu-Sawicki-like modifications of the standard General Relativity action. Therefore, the finding of two-solar mass neutron stars basically imposes no constraint on stable ƒ (R) theories. However, star radii are in general smaller than in General Relativity, which can give an observational handle on such classes of models at the astrophysical level. Both larger masses and smaller matter radii are due to much of the apparent effective energy residing in the outer metric for scalar-tensor theories. Finally, because the ƒ (R) neutron star masses can be much larger than General Relativity counterparts, the total energy available for radiating gravitational waves could be of order several solar masses, and thus a merger of these stars constitutes an interesting wave source.

Self-bound models of compact stars and recent mass-radius measurements

The exact composition of a specific class of compact stars, historically referred to as ""neutron stars,'' is still quite unknown. Possibilities ranging from hadronic to quark degrees of freedom, including self-bound versions of the latter, have been proposed. We specifically address the suitability of strange star models (including pairing interactions) in this work, in the light of new measurements available for four compact stars. The analysis shows that these data might be explained by such an exotic equation of state, actually selecting a small window in parameter space, but still new precise measurements and also further theoretical developments are needed to settle the subject.

Simple analytic model for astrophysical S factors

We propose a physically transparent analytic model of astrophysical S factors as a function of a center-of-mass energy E of colliding nuclei (below and above the Coulomb barrier) for nonresonant fusion reactions. For any given reaction, the S(E) model contains four parameters [two of which approximate the barrier potential, U(r)]. They are easily interpolated along many reactions involving isotopes of the same elements; they give accurate practical expressions for S(E) with only several input parameters for many reactions. The model reproduces the suppression of S(E) at low energies (of astrophysical importance) due to the shape of the low-r wing of U(r). The model can be used to reconstruct U(r) from computed or measured S(E). For illustration, we parametrize our recent calculations of S(E) (using the Sao Paulo potential and the barrier penetration formalism) for 946 reactions involving stable and unstable isotopes of C, O, Ne, and Mg (with nine parameters for all reactions involving many isotopes of the same elements, e. g., C+O). In addition, we analyze astrophysically important (12)C+(12)C reaction, compare theoretical models with experimental data, and discuss the problem of interpolating reliably known S(E) values to low energies (E less than or similar to 2-3 MeV).

Paramagnetism of caesium titanium alum

Caesium titanium alum, CsTi(SO4)(2) . 12H(2)O, is a beta alum and exhibits a large trigonal field and a dynamic Jahn-Teller effect. Exact calculations of the linear (2)T(2)xe Jahn-Teller coupling show that in the strict S-6 Site symmetry the ground multiplet consists of a Kramers doublet 2 Gamma(6) with magnetic splitting factors g(parallel to)=1.1 and g perpendicular to=0, a Gamma(4) Gamma(5) doublet at similar to 60 cm(-1) with g(parallel to)=2.51 and g(perpendicular to)=0.06 and another Gamma(4) Gamma(5) doublet at similar to 270 cm(-1) with g(parallel to)=1.67 and g(perpendicular to)=1.83. The controversial g values observed below 4.2 K, g(parallel to)=1.25 and g(perpendicular to)=1.14, are shown to arise from low symmetry distortions. These distortions couple the vibronic levels and induce into the ground state the off-diagonal axial Zeeman interaction that exists between the first excited and the ground vibronic levels. (C) 1997 American Institute of Physics.

Compact Stars and Magnetized CFL Matter

The stability of the color flavor locked phase in the presence of a strong magnetic field is investigated within the phenomenological MIT bag model. It is found that the minimum value of the energy per baryon in a color flavor locked state at vanishing pressure is lower than the corresponding one for unpaired magnetized strange quark matter and, as the magnetic field increases, the energy per baryon decreases. This implies that magnetized color flavor locked matter is more stable and could become the ground state inside neutron stars. The anisotropy of the pressures is discussed. The mass-radius relation for such stars is also studied.

Magnetized color flavor locked state and compact stars

The stability of the color flavor locked phase in the presence of a strong magnetic field is investigated within the phenomenological MIT bag model, taking into account the variation of the strange quark mass, the baryon density, the magnetic field, as well as the bag and gap parameters. It is found that the minimum value of the energy per baryon in a color flavor locked state at vanishing pressure is lower than the corresponding one for unpaired magnetized strange quark matter and, as the magnetic field increases, the energy per baryon decreases. This implies that magnetized colorflavor locked matter is more stable and could become the ground state inside neutron stars. The mass-radius relation for such stars is also studied.

Multiwavelength Observations of EXO 0748-676. II. Emission-Line Behavior

We present optical and ultraviolet spectra, light curves, and Doppler tomograms of the low-mass X-ray binary EXO 0748-676. Using an extensive set of 15 emission-line tomograms, we show that, along with the usual emission from the stream and ``hot spot,'' there is extended nonaxisymmetric emission from the disk rim. Some of the emission and Hα and Hβ absorption features lend weight to the hypothesis that part of the stream overflows the disk rim and forms a two phase medium. The data are consistent with a 1.35 Msolar neutron star with a main-sequence companion and hence a mass ratio q~0.34.

Are the magnetic fields of millisecond pulsars similar to 10(8) G?

It is generally assumed that the magnetic fields of millisecond pulsars (MSPs) are similar to 10(8) G. We argue that this may not be true and the fields may be appreciably greater. We present six evidences for this: (1) The similar to 10(8)G field estimate is based on magnetic dipole emission losses which is shown to be questionable; (2) The MSPs in low mass X-ray binaries (LMXBs) are claimed to have < 10(11) G on the basis of a Rayleygh-Taylor instability accretion argument. We show that the accretion argument is questionable and the upper limit 10(11) G may be much higher; (3) Low magnetic field neutron stars have difficulty being produced in LMXBs; (4) MSPs may still be accreting indicating a much higher magnetic field; (5) The data that predict similar to 10(8) G for MSPs also predict ages on the order of, and greater than, ten billion years, which is much greater than normal pulsars. If the predicted ages are wrong, most likely the predicted similar to 10(8) G fields of MSPs are wrong; (6) When magnetic fields are measured directly with cyclotron lines in X-ray binaries, fields a parts per thousand << 10(8) G are indicated. Other scenarios should be investigated. One such scenario is the following. Over 85% of MSPs are confirmed members of a binary. It is possible that all MSPs are in large separation binaries having magnetic fields > 10(8) G with their magnetic dipole emission being balanced by low level accretion from their companions.

Origin of 10(15)-10(16) G magnetic fields in the central engine of gamma ray bursts

Various authors have suggested that the gamma-ray burst (GRB) central engine is a rapidly rotating, strongly magnetized, (similar to 10(15)-10(16) G) compact object. The strong magnetic field can accelerate and collimate the relativistic flow and the rotation of the compact object can be the energy source of the GRB. The major problem in this scenario is the difficulty of finding an astrophysical mechanism for obtaining such intense fields. Whereas, in principle, a neutron star could maintain such strong fields, it is difficult to justify a scenario for their creation. If the compact object is a black hole, the problem is more difficult since, according to general relativity it has ""no hair"" (i.e., no magnetic field). Schuster, Blackett, Pauli, and others have suggested that a rotating neutral body can create a magnetic field by non-minimal gravitational-electromagnetic coupling (NMGEC). The Schuster-Blackett form of NMGEC was obtained from the Mikhail and Wanas`s tetrad theory of gravitation (MW). We call the general theory NMGEC-MW. We investigate here the possible origin of the intense magnetic fields similar to 10(15)-10(16) G in GRBs by NMGEC-MW. Whereas these fields are difficult to explain astrophysically, we find that they are easily explained by NMGEC-MW. It not only explains the origin of the similar to 10(15)-10(16) G fields when the compact object is a neutron star, but also when it is a black hole.

Astrophysical S factors for fusion reactions involving C, O, Ne, and Mg isotopes

Using the Sao Paulo potential and the barrier penetration formalism we have calculated the astrophysical factor S(E) for 946 fusion reactions involving stable and neutron-rich isotopes of C, O, Ne, and Mg for center-of-mass energies E varying from 2 to approximate to 18-30 MeV (covering the range below and above the Coulomb barrier). We have parameterized the energy dependence, S(E), by an accurate universal 9-parameter analytic expression and present tables of fit parameters for all the reactions. We also discuss the reduced 3-parameter version of our fit which is highly accurate at energies below the Coulomb barrier, and outline the procedure for calculating the reaction rates. The results can be easily converted to thermonuclear or pycnonuclear reaction rates to simulate various nuclear burning phenomena, in particular, stellar burning at high temperatures and nucleosynthesis in high density environments. (C) 2010 Elsevier Inc. All rights reserved

Low-temperature structural effects in the (TMTSF)2PF6 and AsF6 Bechgaard salts

We present a detailed low-temperature investigation of the statics and dynamics of the anions and methyl groups in the organic conductors (TMTSF) 2PF6 and (TMTSF)2AsF6 (TMTSF: tetramethyl-tetraselenafulvalene). The 4 K neutron-scattering structure refinement of the fully deuterated (TMTSF)2PF6-D12 salt allows locating precisely the methyl groups at 4 K. This structure is compared to the one of the fully hydrogenated (TMTSF)2PF6-H12 salt previously determined at the same temperature. Surprisingly, it is found that deuteration corresponds to the application of a negative pressure of 5×102 MPa to the H12 salt. Accurate measurements of the Bragg intensity show anomalous thermal variations at low temperature both in the deuterated PF 6 and AsF6 salts. Two different thermal behaviors have been distinguished. Small Bragg-angle measurements reflect the presence of low-frequency modes at characteristic energies θE = 8.3 K and θE = 6.7 K for the PF6-D12 and AsF6-D12 salts, respectively. These modes correspond to the low-temperature methyl group motion. Large Bragg-angle measurements evidence an unexpected structural change around 55 K, which probably corresponds to the linkage of the anions to the methyl groups via the formation of F...D-CD2 bonds observed in the 4 K structural refinement. Finally we show that the thermal expansion coefficient of (TMTSF)2PF6 is dominated by the librational motion of the PF6 units. We quantitatively analyze the low-temperature variation of the lattice expansion via the contribution of Einstein oscillators, which allows us to determine for the first time the characteristic frequency of the PF6 librations: θE ≈ 50 K and θE = 76 K for the PF6-D12 and PF6-H12 salts, respectively. © 2013 American Physical Society.