General relativistic effects on neutrino-driven wind from young, hot neutron star and the r-process nucleosynthesis

Neutrino-driven wind from young hot neutron star, which is formed by supernova explosion, is the most promising candidate site for r-process nucleosynthesis. We study general relativistic effects on this wind in Schwarzschild geometry in order to look for suitable conditions for a successful r-process nucleosynthesis. It is quantitatively discussed that the general relativistic effects play a significant role in increasing entropy and decreasing dynamic time scale of the neutrino-driven wind. Exploring wide parameter region which determines the expansion dynamics of the wind, we find interesting physical conditions which lead to successful r-process nucleosynthesis. The conditions which we found realize in the neutrino-driven wind with very short dynamic time scale $\tau_{\rm dyn} \sim 6$ ms and relatively low entropy $S \sim 140$. We carry out the $\alpha$-process and r-process nucleosynthesis calculation on these conditions by the use of our single network code including over 3000 isotopes, and confirm quantitatively that the second and third r-process abundance peaks are produced in the neutrino-driven wind.Comment: Accepted for publication in Ap

Supernovae versus Neutron Star Mergers as the Major r-Process Sources

I show that recent observations of r-process abundances in metal-poor stars are difficult to explain if neutron star mergers (NSMs) are the major r-process sources. In contrast, such observations and meteoritic data on Hf182 and I129 in the early solar system support a self-consistent picture of r-process enrichment by supernovae (SNe). While further theoretical studies of r-process production and enrichment are needed for both SNe and NSMs, I emphasize two possible direct observational tests of the SN r-process model: gamma rays from decay of r-process nuclei in SN remnants and surface contamination of the companion by SN r-process ejecta in binaries.Comment: 5 pages, to appear in ApJ

Accretion dynamics in neutron star black hole binaries

We perform three-dimensional, Newtonian hydrodynamic simulations with a nuclear equation of state to investigate the accretion dynamics in neutron star black hole systems. We find as a general result that non-spinning donor stars yield larger circularization radii than corotating donors. Therefore, the matter from a neutron star without spin will more likely settle into an accretion disk outside the Schwarzschild radius. With the used stiff equation of state we find it hard to form an accretion disk that is promising to launch a gamma-ray burst. In all relevant cases the core of the neutron star survives and keeps orbiting the black hole as a mini neutron star for the rest of the simulation time (up to several hundred dynamical neutron star times scales). The existence of this mini neutron star leaves a clear imprint on the gravitational wave signal which thus can be used to probe the physics at supra-nuclear densities.Comment: submitted to MNRAS, 23 pages, 16 figure

On the Dynamics of Proto-Neutron Star Winds and r-Process Nucleosynthesis

We study here the formation of heavy r-process nuclei in the high-entropy environment of rapidly expanding neutrino-driven winds from compact objects. In particular, we explore the sensitivity of the element creation in the A>130 region to the low-temperature behavior of the outflows. For this purpose we employ a simplified model of the dynamics and thermodynamical evolution for radiation dominated, adiabatic outflows. It consists of a first stage of fast, exponential cooling, followed by a second phase of slower evolution, either assuming constant density and temperature or a power-law decay of these quantities. These cases are supposed to capture the most relevant effects of a strong deceleration or decreasing acceleration of the transsonic outflows, respectively, e.g. in a wind termination shock caused by the collision with the slower, preceding supernova ejecta. We find that not only the transition temperature between the two expansion phases can make a big difference in the formation of the platinum peak, but also the detailed cooling law during the later phase. Unless the transition temperature and corresponding (free neutron) density become too small (T < 2*10^8 K), a lower temperature or faster temperature decline during this phase allow for a stronger appearance of the third abundance peak. Since the nuclear photodisintegration rates between ~2*10^8 K and ~10^9 K are more sensitive to the temperature than the n-capture rates are to the free neutron density, a faster cooling in this temperature regime shifts the r-process path closer to the n-drip line. With low (gamma,n)- but high beta-decay rates, the r-processing then does not proceed through a (gamma,n)-(n,gamma) equilibrium but through a quasi-equilibrium of (n,gamma)-reactions and beta-decays, as recently also pointed out by Wanajo.Comment: 18 pages, 14 figures with 25 eps plots; referee comments included; accepted by Astronomy & Astrophysic

Neutrino signatures and the neutrino-driven wind in Binary Neutron Star Mergers

We present VULCAN/2D multi-group flux-limited-diffusion radiation hydrodynamics simulations of binary neutron star (BNS) mergers, using the Shen equation of state, covering ~100 ms, and starting from azimuthal-averaged 2D slices obtained from 3D SPH simulations of Rosswog & Price for 1.4 Msun (baryonic) neutron stars with no initial spins, co-rotating spins, and counter-rotating spins. Snapshots are post-processed at 10 ms intervals with a multi-angle neutrino-transport solver. We find polar-enhanced neutrino luminosities, dominated by $\bar{\nu}_e$ and ``$\nu_\mu$'' neutrinos at peak, although $\nu_e$ emission may be stronger at late times. We obtain typical peak neutrino energies for $\nu_e$, $\bar{\nu}_e$, and ``$\nu_\mu$'' of ~12, ~16, and ~22 MeV. The super-massive neutron star (SMNS) formed from the merger has a cooling timescale of ~1 s. Charge-current neutrino reactions lead to the formation of a thermally-driven bipolar wind with ~10$^{-3}$ Msun/s, baryon-loading the polar regions, and preventing any production of a GRB prior to black-hole formation. The large budget of rotational free energy suggests magneto-rotational effects could produce a much greater polar mass loss. We estimate that ~10$^{-4}$ Msun of material with electron fraction in the range 0.1-0.2 become unbound during this SMNS phase as a result of neutrino heating. We present a new formalism to compute the $\nu_i\bar{\nu}_i$ annihilation rate based on moments of the neutrino specific intensity computed with our multi-angle solver. Cumulative annihilation rates, which decay as $t^{-1.8}$, decrease over our 100 ms window from a few 10$^{50}$ to ~10$^{49}$ erg/s, equivalent to a few 10$^{54}$ to ~10$^{53}$ $e^-e^+$ pairs per second.Comment: 23 pages, 20 figures, 2 tables, submitted to ApJ, high resolution version of the paper available at http://hermes.as.arizona.edu/~luc/ms.pd

Determining the main-sequence mass of Type II supernova progenitors

We present radiation-hydrodynamics simulations of core-collapse supernova (SN) explosions, artificially generated by driving a piston at the base of the envelope of a rotating or non-rotating red-supergiant progenitor star. We search for trends in ejecta kinematics in the resulting Type II-Plateau (II-P) SN, exploring dependencies with explosion energy and pre-SN stellar-evolution model. We recover the trivial result that larger explosion energies yield larger ejecta velocities in a given progenitor. However, we emphasise that for a given explosion energy, the increasing helium-core mass with main-sequence mass of such Type II-P SN progenitors leads to ejection of core-embedded oxygen-rich material at larger velocities. We find that the photospheric velocity at 15d after shock breakout is a good and simple indicator of the explosion energy in our selected set of pre-SN models. This measurement, combined with the width of the nebular-phase OI6303-6363A line, can be used to place an upper-limit on the progenitor main-sequence mass. Using the results from our simulations, we find that the current, but remarkably scant, late-time spectra of Type II-P SNe support progenitor main-sequence masses inferior to ~20Msun and thus, corroborate the inferences based on the direct, but difficult, progenitor identification in pre-explosion images. The narrow width of OI6303-6363A in Type II-P SNe with nebular spectra does not support high-mass progenitors in the range 25-30Msun. Combined with quantitative spectroscopic modelling, such diagnostics offer a means to constrain the main-sequence mass of the progenitor, the mass fraction of the core ejected, and thus, the mass of the compact remnant formed.Comment: accepted to MNRA

Coalescing neutron stars - a step towards physical models III. Improved numerics and different neutron star masses and spins

(Abridged) In this paper we present a compilation of results from our most advanced neutron star merger simulations, including a description of the employed numerical procedures and a more complete overview over a large number of computed models. The three-dimensional hydrodynamic simulations were done with a code based on the Piecewise Parabolic Method with up to five levels of nested Cartesian grids. The simulations are basically Newtonian, but gravitational-wave emission and the corresponding back-reaction are taken into account. The use of a physical nuclear equation of state allows us to follow the thermodynamic history of the stellar medium and to compute the energy and lepton number loss due to the emission of neutrinos. The computed models differ concerning the neutron star masses and mass ratios, the neutron star spins, the numerical resolution expressed by the cell size of the finest grid and the number of grid levels, and the calculation of the temperature from the solution of the entropy equation instead of the energy equation. Our simulations show that the details of the gravitational-wave emission are still sensitive to the numerical resolution, even in our highest-quality calculations. The amount of mass which can be ejected from neutron star mergers depends strongly on the angular momentum of the system. Our results do not support the initial conditions of temperature and proton-to-nucleon ratio assumed in recent work for producing a solar r-process pattern for nuclei around and above the A approx 130 peak. The improved models confirm our previous conclusion that gamma-ray bursts are not powered by neutrino emission during the dynamical phase of the merging of two neutron stars.Comment: accpeted by A&A; some clarifying text changes due to referee comment

Conformally Flat Smoothed Particle Hydrodynamics: Application to Neutron Star Mergers

We present a new 3D SPH code which solves the general relativistic field + hydrodynamics equations in the conformally flat approximation. Several test cases are considered to test different aspects of the code. We finally apply then the code to the coalescence of a neutron star binary system. The neutron stars are modeled by a polytropic equation of state (EoS) with adiabatic indices $\Gamma=2.0$, $\Gamma=2.6$ and $\Gamma=3.0$. We calculate the gravitational wave signals, luminosities and frequency spectra by employing the quadrupole approximation for emission and back reaction in the slow motion limit. In addition, we consider the amount of ejected mass.Comment: 23 pages, 12 figures. Accepted for publication in Phys. Rev. D. v3: Final Versio

General relativistic simulations of pasive-magneto-rotational core collapse with microphysics

This paper presents results from axisymmetric simulations of magneto-rotational stellar core collapse to neutron stars in general relativity using the passive field approximation for the magnetic field. These simulations are performed using a new general relativistic numerical code specifically designed to study this astrophysical scenario. The code is based on the conformally-flat approximation of Einstein's field equations and conservative formulations of the magneto-hydrodynamics equations. The code has been recently upgraded to incorporate a tabulated, microphysical equation of state and an approximate deleptonization scheme. This allows us to perform the most realistic simulations of magneto-rotational core collapse to date, which are compared with simulations employing a simplified (hybrid) equation of state, widely used in the relativistic core collapse community. Furthermore, state-of-the-art (unmagnetized) initial models from stellar evolution are used. In general, stellar evolution models predict weak magnetic fields in the progenitors, which justifies our simplification of performing the computations under the approach that we call the passive field approximation for the magnetic field. Our results show that for the core collapse models with microphysics the saturation of the magnetic field cannot be reached within dynamical time scales by winding up the poloidal magnetic field into a toroidal one. We estimate the effect of other amplification mechanisms including the magneto-rotational instability (MRI) and several types of dynamos.Comment: 25 pages, 15 figures, accepted for publication in Astronomy & Astrophysics July 31, 2007. Added 1 figure and a new subsectio

Relativistic electron beams above thunderclouds

Non-luminous relativistic electron beams above thunderclouds have been detected by the radio signals of low frequency &amp;sim;40–400 kHz which they radiate. The electron beams occur &amp;sim;2–9 ms after positive cloud-to-ground lightning discharges at heights between &amp;sim;22–72 km above thunderclouds. Intense positive lightning discharges can also cause sprites which occur either above or prior to the electron beam. One electron beam was detected without any luminous sprite which suggests that electron beams may also occur independently of sprites. Numerical simulations show that beams of electrons partially discharge the lightning electric field above thunderclouds and thereby gain a mean energy of &amp;sim;7 MeV to transport a total charge of &amp;sim;−10 mC upwards. The impulsive current &amp;sim;3 &amp;times; 10&lt;sup&gt;&amp;minus;3&lt;/sup&gt; Am&lt;sup&gt;−2&lt;/sup&gt; associated with relativistic electron beams above thunderclouds is directed downwards and needs to be considered as a novel element of the global atmospheric electric circuit