BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

F. -K. Thielemann

F. -K. Thielemann appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.SRastro-phastro-ph.HEnucl-thnucl-exastro-ph.GAastro-ph.IMastro-ph.COastro-ph.EPhep-phphysics.hist-ph

Catalog footprint

What is connected

40works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

40 published item(s)

preprint2022arXiv

Enrichment of the Galactic disc with neutron-capture elements: Gd, Dy, and Th

The study of the origin of heavy elements is one of the main goals of nuclear astrophysics. In this paper, we present new observational data for the heavy $r$-process elements gadolinium (Gd, Z=64), dysprosium (Dy, Z=66) and thorium (Th, Z=90) in a sample of 276 Galactic disc stars ( --1.0$<$[Fe/H]$<$+0.3). The stellar spectra have a high resolution of 42,000 and 75,000, and the signal-to-noise ratio higher than 100. The LTE abundances of Gd, Dy and Th have been determined by comparing the observed and synthetic spectra for three Gd lines (149 stars), four Dy lines (152 stars) and the Th line at 4019.13 A (170 stars). For about 70% of the stars in our sample Gd and Dy are measured for the first time, and Th for 95% of the stars. Typical errors vary from 0.07 to 0.16 dex. This paper provides the first extended set of Th observations in the Milky Way disc. Together with europium (Eu, Z = 63) data from our previous studies, we have compared these new observations with nucleosynthesis predictions and Galactic Chemical Evolution simulations. We confirm that [Gd/Fe] and [Dy/Fe] show the same behavior of Eu. We study with GCE simulations the evolution of [Th/Fe] in comparison with [Eu/Fe], showing that unlike Eu either the Th production is metallicity dependent in case of a unique source of the r-process in the Galaxy, or the frequency of the Th-rich r-process source is decreasing with the increasing of [Fe/H].

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Women Scientists Who Made Nuclear Astrophysics

Female role models reduce the impact on women of stereotype threat, i.e., of being at risk of conforming to a negative stereotype about one's social, gender, or racial group. This can lead women scientists to underperform or to leave their scientific career because of negative stereotypes such as, not being as talented or as interested in science as men. Sadly, history rarely provides role models for women scientists; instead, it often renders these women invisible. In response to this situation, we present a selection of twelve outstanding women who helped to develop nuclear astrophysics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Galactic evolution of rapid neutron capture process abundances: the inhomogeneous approach

For the origin of heavy r-process elements, different sources have been proposed, e.g., core-collapse supernovae or neutron star mergers. Old metal-poor stars carry the signature of the astrophysical source(s). Among the elements dominantly made by the r-process, europium (Eu) is relatively easy to observe. In this work we simulate the evolution of europium in our galaxy with the inhomogeneous chemical evolution model 'ICE', and compare our results with spectroscopic observations. We test the most important parameters affecting the chemical evolution of Eu: (a) for neutron star mergers the coalescence time scale of the merger ($t_{\mathrm{coal}}$) and the probability to experience a neutron star merger event after two supernova explosions occurred and formed a double neutron star system ($P_{\mathrm{NSM}}$) and (b) for the sub-class of magneto-rotationally driven supernovae ("Jet-SNe"), their occurrence rate compared to standard supernovae ($P_{\mathrm{Jet-SN}}$). We find that the observed [Eu/Fe] pattern in the galaxy can be reproduced by a combination of neutron star mergers and magneto-rotationally driven supernovae as r-process sources. While neutron star mergers alone seem to set in at too high metallicities, Jet-SNe provide a cure for this deficiency at low metallicities. Furthermore, we confirm that local inhomogeneities can explain the observed large spread in the europium abundances at low metallicities. We also predict the evolution of [O/Fe] to test whether the spread in $α$-elements for inhomogeneous models agrees with observations and whether this provides constraints on supernova explosion models and their nucleosynthesis.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Mn abundances in the stars of the Galactic disc with metallicities -1.0 < [Fe/H] < 0.3

In this work we present and discuss the observations of the Mn abundances for 247 FGK dwarfs, located in the Galactic disc with metallicity -1<Fe/H]<+0.3. The observed stars belong to the substructures of the Galaxy thick and thin discs, and to the Hercules stream. The observations were conducted using the 1.93 m telescope at Observatoire de Haute-Provence (OHP, France) equipped with the echelle type spectrographs ELODIE and SOPHIE. The abundances were derived under the LTE approximation, with an average error for the [Mn/Fe] ratio of 0.10 dex. For most of the stars in the sample Mn abundances are not available in the literature. We obtain an evolution of [Mn/Fe] ratio with the metallicity [Fe/H] consistent with previous data compilations. In particular, within the metallicity range covered by our stellar sample the [Mn/Fe] ratio is increasing with the increase of metallicity. This due to the contribution to the Galactic chemical evolution of Mn and Fe from thermonuclear supernovae. We confirm the baseline scenario where most of the Mn in the Galactic disc and in the Sun is made by thermonuclear supernovae. In particular, the effective contribution from core-collapse supernovae to the Mn in the Solar system is about 10-20%. However, present uncertainties affecting the production of Mn and Fe in thermonuclear supernovae are limiting the constraining power of the observed [Mn/Fe] trend in the Galactic discs on, e.g., the frequency of different thermonuclear supernovae populations. The different production of these two elements in different types of thermonuclear supernovae needs to be disentangled by the dependence of their relative production on the metallicity of the supernova progenitor.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Primordial heavy element production

A number of possible mechanisms have been suggested to generate density inhomogeneities in the early Universe which could survive until the onset of primordial nucleosynthesis and generate neutron-rich regions. In this work we are not concerned with how the inhomogeneities were generated but we want to focus on the effect of such inhomogeneities on primordial nucleosynthesis. One of the proposed signatures of inhomogeneity, the synthesis of very heavy elements by neutron capture in a primordial r-process, was analyzed for varying baryon to photon ratios $η$ and fluctuation length scales $L$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Pushing 1D CCSNe to explosions: model and SN 1987A

We report on a method, PUSH, for triggering core-collapse supernova explosions of massive stars in spherical symmetry. We explore basic explosion properties and calibrate PUSH such that the observables of SN1987A are reproduced. Our simulations are based on the general relativistic hydrodynamics code AGILE combined with the detailed neutrino transport scheme IDSA for electron neutrinos and ALS for the muon and tau neutrinos. To trigger explosions in the otherwise non-exploding simulations, we rely on the neutrino-driven mechanism. The PUSH method locally increases the energy deposition in the gain region through energy deposition by the heavy neutrino flavors. Our setup allows us to model the explosion for several seconds after core bounce. We explore the progenitor range 18-21M$_{\odot}$. Our studies reveal a distinction between high compactness (HC) and low compactness (LC) progenitor models, where LC models tend to explore earlier, with a lower explosion energy, and with a lower remnant mass. HC models are needed to obtain explosion energies around 1 Bethe, as observed for SN1987A. However, all the models with sufficiently high explosion energy overproduce $^{56}$Ni. We conclude that fallback is needed to reproduce the observed nucleosynthesis yields. The nucleosynthesis yields of $^{57-58}$Ni depend sensitively on the electron fraction and on the location of the mass cut with respect to the initial shell structure of the progenitor star. We identify a progenitor and a suitable set of PUSH parameters that fit the explosion properties of SN1987A when assuming 0.1M$_{\odot}$ of fallback. We predict a neutron star with a gravitational mass of 1.50M$_{\odot}$. We find correlations between explosion properties and the compactness of the progenitor model in the explored progenitors. However, a more complete analysis will require the exploration of a larger set of progenitors with PUSH.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The long-term evolution of neutron star merger remnants - I. The impact of r-process nucleosynthesis

We follow the longterm evolution of the dynamic ejecta of neutron star mergers for up to 100 years and over a density range of roughly 40 orders of magnitude. We include the nuclear energy input from the freshly synthesized, radioactively decaying nuclei in our simulations and study its effects on the remnant dynamics. Although the nuclear heating substantially alters the longterm evolution, we find that running nuclear networks over purely hydrodynamic simulations (i.e. without heating) yields actually acceptable nucleosynthesis results. The main dynamic effect of the radioactive heating is to quickly smooth out inhomogeneities in the initial mass distribution, subsequently the evolution proceeds self-similarly and after 100 years the remnant still carries the memory of the initial binary mass ratio. We also explore the nucleosynthetic yields for two mass ejection channels. The dynamic ejecta very robustly produce "strong" r-process elements with $A > 130$ with a pattern that is essentially independent of the details of the merging system. From a simple model we find that neutrino-driven winds yield "weak" r-process contributions with $50 < A < 130$ whose abundance patterns vary substantially between different merger cases. This is because their electron fraction, set by the ratio of neutrino luminosities, varies considerably from case to case. Such winds do not produce any $^{56}{\rm Ni}$, but a range of radioactive isotopes that are long-lived enough to produce a second, radioactively powered electromagnetic transient in addition to the "macronova" from the dynamic ejecta. While our wind model is very simple, it nevertheless demonstrates the potential of such neutrino-driven winds for electromagnetic transients and it motivates further, more detailed neutrino-hydrodynamic studies. The properties of the mentioned transients are discussed in more detail in a companion paper.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A comparative study of statistical models for nuclear equation of state of stellar matter

We compare three different statistical models for the equation of state (EOS) of stellar matter at subnuclear densities and temperatures (0.5-10 MeV) expected to occur during the collapse of massive stars and supernova explosions. The models introduce the distributions of various nuclear species in nuclear statistical equilibrium, but use somewhat different nuclear physics inputs. It is demonstrated that the basic thermodynamical quantities of stellar matter under these conditions are similar, except in the region of high densities and low temperatures. We demonstrate that mass and isotopic distributions have considerable differences related to the different assumptions of the models on properties of nuclei at these stellar conditions. Overall, the three models give similar trends, but the details reflect the uncertainties related to the modeling of medium effects, such as the temperature and density dependence of surface and bulk energies of heavy nuclei, and the nuclear shell structure effects. We discuss importance of new physics inputs for astrophysical calculations from experimental data obtained in intermediate energy heavy-ion collisions, in particular, the similarities of the conditions reached during supernova explosions and multifragmentation reactions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Abundances of neutron-capture elements in stars of the galactic disk substructures

The aim of this work is to present and discuss the observations of the iron peak (Fe, Ni) and neutron-capture element (Y, Zr, Ba, La, Ce, Nd, Sm, and Eu) abundances for 276 FGK dwarfs, located in the galactic disk with metallicity -1 < [Fe/H] < +0.3. Atmospheric parameters and chemical composition of the studied stars were determined from an high resolution, high signal-to-noise echelle spectra obtained with the echelle spectrograph ELODIE at the Observatoire de Haute-Provence (France). Effective temperatures were estimated by the line depth ratio method and from the H_α line-wing fitting. Surface gravities (log g) were determined by parallaxes and the ionization balance of iron. Abundance determinations were carried out using the LTE approach, taking the hyperfine structure for Eu into account, and the abundance of Ba was computed under the NLTE approximation. We are able to assign most of the stars in our sample to the substructures of the Galaxy thick disk, thin disk, or Hercules stream according to their kinematics. The classification of 27 stars is uncertain. For most of the stars in the sample, the abundances of neutron-capture elements have not been measured earlier. For all of them, we provide the chemical composition and discuss the contribution from different nucleosynthesis processes.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Production of carbon-rich presolar grains from massive stars

About a year after core collapse supernova, dust starts to condense in the ejecta. In meteorites, a fraction of C-rich presolar grains (e.g., silicon carbide (SiC) grains of Type-X and low density graphites) are identified as relics of these events, according to the anomalous isotopic abundances. Several features of these abundances remain unexplained and challenge the understanding of core-collapse supernovae explosions and nucleosynthesis. We show, for the first time, that most of the measured C-rich grain abundances can be accounted for in the C-rich material from explosive He burning in core-collapse supernovae with high shock velocities and consequent high temperatures. The inefficiency of the $^{12}$C($α$,$γ$)$^{16}$O reaction relative to the rest of the $α$-capture chain at $T > 3.5\times10^8 \mathrm{K}$ causes the deepest He-shell material to be carbon rich and silicon rich, and depleted in oxygen. The isotopic ratio predictions in part of this material, defined here as the C/Si zone, are in agreement with the grain data. The high-temperature explosive conditions that our models reach at the bottom of the He shell, can also be representative of the nucleosynthesis in hypernovae or in the high-temperature tail of a distribution of conditions in asymmetric supernovae. Finally, our predictions are consistent with the observation of large $^{44}$Ca/$^{40}$Ca observed in the grains. This is due to the production of $^{44}$Ti together with $^{40}$Ca in the C/Si zone, and/or to the strong depletion of $^{40}$Ca by neutron captures.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Silicon carbide grains of type C provide evidence for the production of the unstable isotope $^{32}$Si in supernovae

Carbon-rich grains are observed to condense in the ejecta of recent core-collapse supernovae, within a year after the explosion. Silicon carbide grains of type X are C-rich grains with isotpic signatures of explosive supernova nucleosynthesis have been found in primitive meteorites. Much rarer silicon carbide grains of type C are a special sub-group of SiC grains from supernovae. They show peculiar abundance signatures for Si and S, isotopically heavy Si and isotopically light S, which appear to to be in disagreement with model predictions. We propose that C grains are formed mostly from C-rich stellar material exposed to lower SN shock temperatures than the more common type X grains. In this scenario, extreme $^{32}$S enrichments observed in C grains may be explained by the presence of short-lived $^{32}$Si ($τ$$_{1/2}$ = 153 years) in the ejecta, produced by neutron capture processes starting from the stable Si isotopes. No mixing from deeper Si-rich material and/or fractionation of Si from S due to molecular chemistry is needed to explain the $^{32}$S enrichments. The abundance of $^{32}$Si in the grains can provide constraints on the neutron density reached during the supernova explosion in the C-rich He shell material. The impact of the large uncertainty of the neutron capture cross sections in the $^{32}$Si region is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Have Superheavy Elements been Produced in Nature?

We discuss the possibility whether superheavy elements can be produced in Nature by the astrophysical rapid neutron capture process. To this end we have performed fully dynamical network r-process calculations assuming an environment with neutron-to-seed ratio large enough to produce superheavy nuclei. Our calculations include two sets of nuclear masses and fission barriers and include all possible fission channels and the associated fission yield distributions. Our calculations produce superheavy nuclei with A ~ 300 that however decay on timescales of days.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Long-term evolution of massive star explosions

We examine simulations of core-collapse supernovae in spherical symmetry. Our model is based on general relativistic radiation hydrodynamics with three-flavor Boltzmann neutrino transport. We discuss the different supernova phases, including the long-term evolution up to 20 seconds after the onset of explosion during which the neutrino fluxes and mean energies decrease continuously. In addition, the spectra of all flavors become increasingly similar, indicating the change from charged- to neutral-current dominance. Furthermore, it has been shown recently by several groups independently, based on sophisticated supernova models, that collective neutrino flavor oscillations are suppressed during the early mass-accretion dominated post-bounce evolution. Here we focus on the possibility of collective flavor flips between electron and non-electron flavors during the later, on the order of seconds, evolution after the onset of an explosion with possible application for the nucleosynthesis of heavy elements.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Neutrino-driven wind simulations and nucleosynthesis of heavy elements

Neutrino-driven winds, which follow core-collapse supernova explosions, present a fascinating nuclear astrophysics problem that requires understanding advanced astrophysics simulations, the properties of matter and neutrino interactions under extreme conditions, the structure and reactions of exotic nuclei, and comparisons against forefront astronomical observations. The neutrino-driven wind has attracted vast attention over the last 20 years as it was suggested to be a candidate for the astrophysics site where half of the heavy elements are produced via the r-process. In this review, we summarize our present understanding of neutrino-driven winds from the dynamical and nucleosynthesis perspectives. Rapid progress has been made during recent years in understanding the wind with improved simulations and better micro physics. The current status of the fields is that hydrodynamical simulations do not reach the extreme conditions necessary for the r-process and the proton or neutron richness of the wind remains to be investigated in more detail. However, nucleosynthesis studies and observations point already to neutrino-driven winds to explain the origin of lighter heavy elements, such as Sr, Y, Zr.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Core-collapse supernova explosions triggered by a quark-hadron phase transition during the early post-bounce phase

We explore explosions of massive stars, which are triggered via the quark-hadron phase transition during the early post bounce phase of core-collapse supernovae. We construct a quark equation of state, based on the bag model for strange quark matter. The transition between the hadronic and the quark phases is constructed applying Gibbs conditions. The resulting quark-hadron hybrid equations of state are used in core-collapse supernova simulations, based on general relativistic radiation hydrodynamics and three flavor Boltzmann neutrino transport in spherical symmetry. The formation of a mixed phase reduces the adiabatic index, which induces the gravitational collapse of the central protoneutron star. The collapse halts in the pure quark phase, where the adiabatic index increases. A strong accretion shock forms, which propagates towards the protoneutron star surface. Due to the density decrease of several orders of magnitude, the accretion shock turns into a dynamic shock with matter outflow. This moment defines the onset of the explosion in supernova models that allow for a quark-hadron phase transition, where otherwise no explosions could be obtained. The shock propagation across the neutrinospheres releases a burst of neutrinos. This serves as a strong observable identification for the structural reconfiguration of the stellar core. The ejected matter expands on a short timescale and remains neutron-rich. These conditions might be suitable for the production of heavy elements via the r-process. The neutron-rich material is followed by proton-rich neutrino-driven ejecta in the later cooling phase of the protoneutron star where the vp-process might occur.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Cross sections for proton-induced reactions on Pd isotopes at energies relevant for the gamma process

Proton-activation reactions on natural and enriched palladium samples were investigated via the activation technique in the energy range of E_p=2.75 MeV to 9 MeV, close to the upper end of the respective Gamow window of the gamma process. We have determined cross sections for 102Pd(p,gamma)103Ag, 104Pd(p,gamma)105Ag, and 105Pd(p,n)105Ag, as well as partial cross sections of 104Pd(p,n)104Ag^g, 105Pd(p,gamma)106Ag^m, 106Pd(p,n)106Ag^m, and 110Pd(p,n)110Ag^m with uncertainties between 3% and 15% for constraining theoretical Hauser-Feshbach rates and for direct use in gamma-process calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

GRIPS - Gamma-Ray Imaging, Polarimetry and Spectroscopy

We propose to perform a continuously scanning all-sky survey from 200 keV to 80 MeV achieving a sensitivity which is better by a factor of 40 or more compared to the previous missions in this energy range. The Gamma-Ray Imaging, Polarimetry and Spectroscopy (GRIPS) mission addresses fundamental questions in ESA's Cosmic Vision plan. Among the major themes of the strategic plan, GRIPS has its focus on the evolving, violent Universe, exploring a unique energy window. We propose to investigate $γ$-ray bursts and blazars, the mechanisms behind supernova explosions, nucleosynthesis and spallation, the enigmatic origin of positrons in our Galaxy, and the nature of radiation processes and particle acceleration in extreme cosmic sources including pulsars and magnetars. The natural energy scale for these non-thermal processes is of the order of MeV. Although they can be partially and indirectly studied using other methods, only the proposed GRIPS measurements will provide direct access to their primary photons. GRIPS will be a driver for the study of transient sources in the era of neutrino and gravitational wave observatories such as IceCUBE and LISA, establishing a new type of diagnostics in relativistic and nuclear astrophysics. This will support extrapolations to investigate star formation, galaxy evolution, and black hole formation at high redshifts.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Non-standard s process in low metallicity massive rotating stars

Context. Rotation is known to affect the nucleosynthesis of light elements in massive stars, mainly by rotation-induced mixing. In particular, rotation boosts the primary nitrogen production. Models of rotating stars are able to reproduce the nitrogen observed in low-Z halo stars. Aims. Here we present the first grid of stellar models for rotating massive stars at low Z, where a full s-process network is used to study the impact of rotation-induced mixing on the nucleosynthesis of heavy elements. Methods. We used the Geneva stellar evolution code that includes an enlarged reaction network with nuclear species up to bismuth to calculate 25 M$_\odot$ models at three different Z and with different initial rotation rates. Results. First, we confirm that rotation-induced mixing leads to a production of primary $^{22}$Ne, which is the main neutron source for the s process in massive stars. Therefore rotation boosts the s process in massive stars at all Z. Second, the neutron-to-seed ratio increases with decreasing Z in models including rotation, which leads to the complete consumption of all iron seeds at Z < 1e-3 by the end of core He-burning. Thus at low Z, the iron seeds are the main limitation for this boosted s process. Third, as Z decreases, the production of elements up to the Ba peak increases at the expense of the elements of the Sr peak. We studied the impact of the initial rotation rate and of the uncertain $^{17}$O$(α,γ)$ rate (which strongly affects the neutron poison strength of $^{16}$O) on our results. This study shows that rotating models can produce significant amounts of elements up to Ba over a wide range of Z. Fourth, compared to the He-core, the primary $^{22}$Ne production in the He-shell is even higher (> 1% in mass fraction at all Z), which could open the door for an explosive neutron capture nucleosynthesis in the He-shell, with a primary neutron source.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Strange matter in core-collapse supernovae

We discuss the possible impact of strange quark matter on the evolution of core-collapse supernovae with emphasis on low critical densities for the quark-hadron phase transition. For such cases the hot proto-neutron star can collapse to a more compact hybrid star configuration hundreds of milliseconds after core-bounce. The collapse triggers the formation of a second shock wave. The latter leads to a successful supernova explosion and leaves an imprint on the neutrino signal. These dynamical features are discussed with respect to their compatibility with recent neutron star mass measurements which indicate a stiff high density nuclear matter equation of state.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Exploring the QCD phase transition in core collapse supernova simulations in spherical symmetry

For finite chemical potential effective models of QCD predict a first order phase transition. In favour for the search of such a phase transition in nature, we construct an equation of state for strange quark matter based on the MIT bag model. We apply this equation of state to highly asymmetric core collapse supernova matter with finite temperatures and large baryon densities. The phase transition is constructed using the general Gibbs conditions, which results in an extended coexistence region between the pure hadronic and pure quark phases in the phase diagram, i.e. the mixed phase. The supernovae are simulated via general relativistic radiation hydrodynamics based on three flavor Boltzmann neutrino transport in spherical symmetry. During the dynamical evolution temperatures above 10 MeV, baryon densities above nuclear saturation density and a proton-to-baryon ratio below 0.2 are obtained. At these conditions the phase transition is triggered which leads to a significant softening of the EoS for matter in the mixed phase. As a direct consequence of the stiffening of the EoS again for matter in the pure quark phase, a shock wave forms at the boundary between the mixed and the pure hadronic phases. This shock is accelerated and propagates outward which releases a burst of neutrinos dominated by electron anti-neutrinos due to the lifted degeneracy of the shock-heated hadronic material. We discuss the radiation-hydrodynamic evolution of the phase transition at the example of several low and intermediate mass Fe-core progenitor stars and illustrate the expected neutrino signal from the phase transition.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Neutron-induced astrophysical reaction rates for translead nuclei

Neutron-induced reaction rates, including fission, are calculated in the temperature range 1.d8 <T (K) < 1.d10 within the framework of the statistical model for targets with atomic number 83 < Z < 119 (from Po to Uuo) from the neutron to the proton drip-line. Four sets of rates have been calculated, utilizing - where possible - consistent nuclear data for neutron separation energies and fission barriers from Thomas-Fermi (TF), Extended Thomas-Fermi plus Strutinsky Integral (ETFSI), Finite-Range Droplet Model (FRDM) and Hartree-Fock-Bogolyubov (HFB) predictions. Tables of calculated values as well as analytic seven parameter fits in the standard REACLIB format are supplied. We also discuss the sensitivity of the rates to the input, aiming at a better understanding of the uncertainties introduced by the nuclear input.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Nucleosynthesis in Two-Dimensional Delayed Detonation Models of Type Ia Supernova Explosions

The nucleosynthetic characteristics of various explosion mechanisms of Type Ia supernovae (SNe Ia) is explored based on three two-dimensional explosion simulations representing extreme cases: a pure turbulent deflagration, a delayed detonation following an approximately spherical ignition of the initial deflagration, and a delayed detonation arising from a highly asymmetric deflagration ignition. Apart from this initial condition, the deflagration stage is treated in a parameter-free approach. The detonation is initiated when the turbulent burning enters the distributed burning regime. This occurs at densities around $10^{7}$ g cm$^{-3}$ -- relatively low as compared to existing nucleosynthesis studies for one-dimensional spherically symmetric models. The burning in these multidimensional models is different from that in one-dimensional simulations as the detonation wave propagates both into unburned material in the high density region near the center of a white dwarf and into the low density region near the surface. Thus, the resulting yield is a mixture of different explosive burning products, from carbon-burning products at low densities to complete silicon-burning products at the highest densities, as well as electron-capture products synthesized at the deflagration stage. In contrast to the deflagration model, the delayed detonations produce a characteristic layered structure and the yields largely satisfy constraints from Galactic chemical evolution. In the asymmetric delayed detonation model, the region filled with electron capture species (e.g., $^{58}$Ni, $^{54}$Fe) is within a shell, showing a large off-set, above the bulk of $^{56}$Ni distribution, while species produced by the detonation are distributed more spherically (abridged).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Protoneutron star evolution and the neutrino driven wind in general relativistic neutrino radiation hydrodynamics simulations

Massive stars end their life in an explosion event with kinetic energies of the order 1 Bethe. Immediately after the explosion has been launched, a region of low density and high entropy forms behind the ejecta which is continuously subject to neutrino heating. The neutrinos emitted from the remnant at the center, the protoneutron star (PNS), heat the material above the PNS surface. This heat is partly converted into kinetic energy and the material accelerates to an outflow that is known as the neutrino driven wind. For the first time, we simulate the collapse, bounce, explosion and the neutrino driven wind phases consistently over more than 20 seconds. Our numerical model is based on spherically symmetric general relativistic radiation hydrodynamics using spectral three flavor Boltzmann neutrino transport. In simulations where no explosions are obtained naturally, we model neutrino driven explosions for low and intermediate mass Fe-core progenitor stars by enhancing the charged current reaction rates. In the case of a special progenitor star, the O-Ne-Mg-core, the explosion in spherical symmetry was obtained without enhanced opacities. The post explosion evolution is in qualitative agreement with static steady-state and parametrized dynamic models of the neutrino driven wind. On the other hand, we find generally smaller neutrino luminosities and mean neutrino energies as well as a different evolutionary behavior of the neutrino luminosities and mean neutrino energies. The neutrino driven wind is proton-rich for more than 10 seconds and the contraction of the PNS differs from the assumptions made for the conditions at the inner boundary in previous neutrino driven wind studies. Despite the moderately large entropies per baryon of about 100 and the fast expansion timescale, the conditions found in our model are unlikely to favor...

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Stellar (n,gamma) cross sections of p-process isotopes PartI: 102Pd, 120Te, 130,132Ba,and 156Dy

We have investigated the (n,gamma) cross sections of p-process isotopes with the activation technique. The measurements were carried out at the Karlsruhe Van de Graaff accelerator using the 7Li(p,n)7Be source for simulating a Maxwellian neutron distribution of kT = 25 keV. Stellar cross section measurements are reported for the light p-process isotopes 102Pd, 120Te, 130,132Ba, and 156Dy. In a following paper the cross sections of 168Yb, 180W, 184Os, 190Pt, and 196Hg will be discussed. The data are extrapolated to p-process energies by including information from evaluated nuclear data libraries. The results are compared to standard Hauser-Feshbach models frequently used in astrophysics.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Strange quark matter in explosive astrophysical systems

Explosive astrophysical systems, such as supernovae or compact star binary mergers, provide conditions where strange quark matter can appear. The high degree of isospin asymmetry and temperatures of several MeV in such systems may cause a transition to the quark phase already around saturation density. Observable signals from the appearance of quark matter can be predicted and studied in astrophysical simulations. As input in such simulations, an equation of state with an integrated quark matter phase transition for a large temperature, density and proton fraction range is required. Additionally, restrictions from heavy ion data and pulsar observation must be considered. In this work we present such an approach. We implement a quark matter phase transition in a hadronic equation of state widely used for astrophysical simulations and discuss its compatibility with heavy ion collisions and pulsar data. Furthermore, we review the recently studied implications of the QCD phase transition during the early post-bounce evolution of core-collapse supernovae and introduce the effects from strong interactions to increase the maximum mass of hybrid stars. In the MIT bag model, together with the strange quark mass and the bag constant, the strong coupling constant $α_s$ provides a parameter to set the beginning and extension of the quark phase and with this the mass and radius of hybrid stars.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The neutrino signal from protoneutron star accretion and black hole formation

We discuss the formation of stellar mass black holes via protoneutron star (PNS) collapse. In the absence of an earlier explosion, the PNS collapses to a black hole due to the continued mass accretion onto the PNS. We present an analysis of the emitted neutrino spectra of all three flavors during the PNS contraction. Special attention is given to the physical conditions which depend on the input physics, e.g. the equation of state (EoS) and the progenitor model. The PNSs are modeled as the central object in core collapse simulations using general relativistic three-flavor Boltzmann neutrino transport in spherical symmetry. The simulations are launched from several massive progenitors of 40 and 50 solar mass. We analyze the electron-neutrino luminosity dependencies and construct a simple approximation for the electron-neutrino luminosity, which depends only on the physical conditions at the electron-neutrinosphere. In addition, we analyze different mu/tau-neutrino pair-reactions separately and compare the differences during the post-bounce phase of failed core collapse supernova explosions of massive progenitors. We also investigate the connection between the increasing mu/tau-neutrino luminosity and the PNS contraction during the accretion phase before black hole formation. Comparing the different post bounce phase of the progenitor models under investigation, we find large differences in the emitted neutrino spectra. These differences and the analysis of the electron-neutrino luminosity indicate a strong progenitor model dependency of the emitted neutrino signal.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

(n,$γ$) Cross Sections of Light p Nuclei -- Towards an Updated Database for the p Process

The nucleosynthesis of elements beyond iron is dominated by the s and r processes. However, a small amount of stable isotopes on the proton-rich side cannot be made by neutron capture and are thought to be produced by photodisintegration reactions on existing seed nuclei in the so-called "p process". So far most of the p-process reactions are not yet accessible by experimental techniques and have to be inferred from statistical Hauser-Feshbach model calculations. The parametrization of these models has to be constrained by measurements on stable proton-rich nuclei. A series of (n,$γ$) activation measurements on p nuclei, related by detailed balance to the respective photodisintegrations, were carried out at the Karlsruhe Van de Graaff accelerator using the $^7$Li(p,n)$^7$Be source for simulating a Maxwellian neutron distribution of kT= 25 keV. We present here preliminary results of our extended measuring program in the mass range between A=74 and A=132, including first experimental (n,$γ$) cross sections of $^{74}$Se, $^{84}$Sr, $^{120}$Te and $^{132}$Ba, and an improved value for $^{130}$Ba. In all cases we find perfect agreement with the recommended MACS predictions from the Bao et al. compilation.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

$p$-Process simulations with a modified reaction library

We have performed $p$-process simulations with the most recent stellar $(n,γ)$ cross sections from the "Karlsruhe Astrophysical Database of Nucleosynthesis in Stars" project (version v0.2, http://nuclear-astrophysics.fzk.de/kadonis). The simulations were carried out with a parametrized supernova type II shock front model (``$γ$ process'') of a 25 solar mass star and compared to recently published results. A decrease in the normalized overproduction factor could be attributed to lower cross sections of a significant fraction of seed nuclei located in the Bi and Pb region around the $N$=126 shell closure.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Angular momentum transport during X-ray bursts on neutron stars: a numerical general relativistic hydrodynamical study

The distribution of angular momentum of the matter during X-ray bursts on neutron stars is studied by means of 3D axi-symmetric general relativistic hydrodynamics. The set of fully general relativistic Navier-Stokes equations is solved implicitly using the implicit solver GR-I-RMHD in combination with a third order spatial and second order temporal advection scheme. The viscous operators are formulated using a Kerr-like metric in the fixed background of a slowly rotating neutron star whose radius coincides with the corresponding last stable orbit. The importance of these operators and their possible simplifications are discussed as well. In the rotating case and depending on the viscosity parameter, $α_\mm{tur}$, it is found that the viscously-initiated fronts at the center of bursts propagate at much faster speed than the fluid motion. These fast fronts act to decouple angular momentum from matter: angular momentum is transported outwards while matter sinks inwards into the deep gravitational well of the neutron star, thereby enhancing the compression of matter necessary for initiating ignition, that subsequently spreads over the whole surface of the neutron star on the viscous time scale. Based on the numerical simulations, we find that a viscosity parameter $α_\mm{tur} = \mathcal{O}(0.1)$ is most suitable for fitting observations of neutron stars during X-ray bursts. It is argued that the spin up observed in the cooling tails of X-ray bursts is a transient phase, which eventually should be followed by a spin down phase. This delay can be attributed to a significant lengthening of the viscous time scale due to rapid cooling of matter in the outer layers.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

R-process nucleosynthesis calculations with complete nuclear physics input

The r-process constitutes one of the major challenges in nuclear astrophysics. Its astrophysical site has not yet been identified but there is observational evidence suggesting that at least two possible sites should contribute to the solar system abundance of r-process elements and that the r-process responsible for the production of elements heavier than Z=56 operates quite robustly producing always the same relative abundances. From the nuclear-physics point of view the r-process requires the knowledge of a large number of reaction rates involving exotic nuclei. These include neutron capture rates, beta-decays and fission rates, the latter for the heavier nuclei produced in the r-process. We have developed for the first time a complete database of reaction rates that in addition to neutron-capture rates and beta-decay half-lives includes all possible reactions that can induce fission (neutron-capture, beta-decay and spontaneous fission) and the corresponding fission yields. In addition, we have implemented these reaction rates in a fully implicit reaction network. We have performed r-process calculations for the neutrino-driven wind scenario to explore whether or not fission can contribute to provide a robust r-process pattern.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Stellar evolution of massive stars at very low metallicities

Recently, measurements of abundances in extremely metal poor (EMP) stars have brought new constraints on stellar evolution models. In an attempt to explain the origin of the abundances observed, we computed pre--supernova evolution models, explosion models and the related nucleosynthesis. In this paper, we start by presenting the pre-SN models of rotating single stars with metallicities ranging from solar metallicity down to almost metal free. We then review key processes in core-collapse and bounce, before we integrate them in a simplistic parameterization for 3D MHD models, which are well underway and allow one to follow the evolution of the magnetic fields during collapse and bounce. Finally, we present explosive nucleosynthesis results including neutrino interactions with matter, which are calculated using the outputs of the explosion models. The main results of the pre-SN models are the following. First, primary nitrogen is produced in large amount in models with an initial metallicity $Z=10^{-8}$. Second, at the same metallicity of $Z=10^{-8}$ and for models with an initial mass larger than about 60 Mo, rotating models may experience heavy mass loss (up to more than half of the initial mass of the star). The chemical composition of these winds can qualitatively reproduce the abundance patterns observed at the surface of carbon-rich EMP stars. Explosive nucleosynthesis including neutrino-matter interactions produce improved abundances for iron group elements, in particular for scandium and zinc. It also opens the way to a new neutrino and proton rich process ($ν$p-process) able to contribute to the nucleosynthesis of elements with A > 64. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Composition of the Innermost Core Collapse Supernova Ejecta

With presently known input physics and computer simulations in 1D, a self-consistent treatment of core collapse supernovae does not yet lead to successful explosions, while 2D models show some promise. Thus, there are strong indications that the delayed neutrino mechanism works combined with a multi-D convection treatment for unstable layers. On the other hand there is a need to provide correct nucleosynthesis abundances for the progressing field of galactic evolution and observations of low metallicity stars. The innermost ejecta is directly affected by the explosion mechanism, i.e. most strongly the yields of Fe-group nuclei for which an induced piston or thermal bomb treatment will not provide the correct yields because the effect of neutrino interactions is not included. We apply parameterized variations to the neutrino scattering cross sections and alternatively, parameterized variations are applied to the neutrino absorption cross sections on nucleons in the ``gain region''. We find that both measures lead to similar results, causing explosions and a Ye>0.5 in the innermost ejected layers, due to the combined effect of a short weak interaction time scale and a negligible electron degeneracy, unveiling the proton-neutron mass difference. We include all weak interactions (electron and positron capture, beta-decay, neutrino and antineutrino capture on nuclei, and neutrino and antineutrino capture on nucleons) and present first nucleosynthesis results for these innermost ejected layers to discuss how they improve predictions for Fe-group nuclei. The proton-rich environment results in enhanced abundances of 45Sc, 49Ti, and 64Zn as requested by chemical evolution studies and observations of low metallicity stars as well as appreciable production of nuclei in the mass range up to A=80.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

A Finite Difference Representation of Neutrino Radiation Hydrodynamics in Spherically Symmetric General Relativistic Space-Time

We present an implicit finite difference representation for general relativistic radiation hydrodynamics in spherical symmetry. Our code, Agile-Boltztran, solves the Boltzmann transport equation for the angular and spectral neutrino distribution functions in self-consistent simulations of stellar core collapse and postbounce evolution. It implements a dynamically adaptive grid in comoving coordinates. Most macroscopically interesting physical quantities are defined by expectation values of the distribution function. We optimize the finite differencing of the microscopic transport equation for a consistent evolution of important expectation values. We test our code in simulations launched from progenitor stars with 13 solar masses and 40 solar masses. ~0.5 s after core collapse and bounce, the protoneutron star in the latter case reaches its maximum mass and collapses further to form a black hole. When the hydrostatic gravitational contraction sets in, we find a transient increase in electron flavor neutrino luminosities due to a change in the accretion rate. The muon- and tauon-neutrino luminosities and rms energies, however, continue to rise because previously shock-heated material with a non-degenerate electron gas starts to replace the cool degenerate material at their production site. We demonstrate this by supplementing the concept of neutrinospheres with a more detailed statistical description of the origin of escaping neutrinos. We compare the evolution of the 13 solar mass progenitor star to simulations with the MGFLD approximation, based on a recently developed flux limiter. We find similar results in the postbounce phase and validate this MGFLD approach for the spherically symmetric case with standard input physics.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Prediction of Astrophysical Reaction Rates: Methods, Data Needs, and Consequences for Nucleosynthesis Studies

The majority of nuclear reactions in astrophysics involve unstable nuclei which are not fully accessible by experiments yet. Therefore, there is high demand for reliable predictions of cross sections and reaction rates by theoretical means. The majority of reactions can be treated in the framework of the statistical model (Hauser-Feshbach). The global parametrizations of the nuclear properties needed for predictions far off stability probe our understanding of the strong force and take it to its limit. The sensitivity of astrophysical scenarios to nuclear inputs is illustrated in the framework of a detailed nucleosynthesis study in type II supernovae. Abundances resulting from calculations in the same explosion model with two different sets of reaction rates are compared. Key reactions and required nuclear information are identified.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Coalescing Neutron Stars: A Solution to the R-Process Problem?

Recent calculations place questions on the ability of supernovae to produce r-process nuclei in the correct amounts. We present results from 3D Newtonian SPH calculations of mergers of equal mass neutron star binaries. We find the amounts of ejected mass to be very promising for an explanation of the observed r-process abundances. Preliminary abundance calculations show that practically all the material is subject to r-process conditions. The calculated abundance patterns reproduce the basic features of the solar r-process abundances very well. Especially the peak around A= 195 is easily reproduced without any tuning of the initial entropies. Thus, eventually all the observed r-process material could be explained by mass ejection during neutron star mergers.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Mass ejection in neutron star mergers

We present the results of 3D Newtonian SPH simulations of the merger of a neutron star binary. The microscopic properties of matter are described by the physical equation of state of Lattimer and Swesty (LS-EOS). To test for the robustness of our results we check the sensitivity to the approximations of our model as well as to the binary system parameters. The main and new result is that for the realistic LS-EOS, depending on the initial spin, between 4e-3 and 4e-2 solar masses of material become unbound. If, as suggested, large parts of this matter consist of r-process nuclei, neutron star mergers could account for the whole observed r-process material in the Galaxy.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

An Approximation for the rp-Process

Hot (explosive) hydrogen burning or the Rapid Proton Capture Process (rp-process) occurs in a number of astrophysical environments. Novae and X-ray bursts are the most prominent ones, but accretion disks around black holes and other sites are candidates as well. The expensive and often multidimensional hydro calculations for such events require an accurate prediction of the thermonuclear energy generation, while avoiding full nucleosynthesis network calculations. In the present investigation we present an approximation scheme applicable in a temperature range which covers the whole range of all presently known astrophysical sites. It is based on the concept of slowly varying hydrogen and helium abundances and assumes a kind of local steady flow by requiring that all reactions entering and leaving a nucleus add up to a zero flux. This scheme can adapt itself automatically and covers situations at low temperatures, characterized by a steady flow of reactions, as well as high temperature regimes where a $(p,γ)-(γ,p)$-equilibrium is established. In addition to a gain of a factor of 15 in computational speed over a full network calculation, and an energy generation accurate to more than 15 %, this scheme also allows to predict correctly individual isotopic abundances. Thus, it delivers all features of a full network at a highly reduced cost and can easily be implemented in hydro calculations.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

The Nuclear Level Density and the Determination of Thermonuclear Rates for Astrophysics

The prediction of cross sections for nuclei far off stability is crucial in the field of nuclear astrophysics. In recent calculations the nuclear level density -- as an important ingredient to the statistical model (Hauser-Feshbach) -- has shown the highest uncertainties. We present a global parametrization of nuclear level densities based on the back-shifted Fermi-Gas formalism. Employment of an energy-dependent level density parameter $a$ and microscopic corrections from a recent FRDM mass formula by Möller et al.\ leads to a highly improved fit of level densities at the neutron-separation energy in the mass range $20\le A \le 245$. The importance of using proper microscopic corrections from mass formulae is emphasized. The resulting level description is well suited for astrophysical applications. The level density can also provide clues to the applicability of the statistical model which is only correct for a high density of excited states. Using the above description one can derive a ``map'' for the applicability of the model for reactions of stable and unstable nuclei with neutral and charged particles.

0 citations0 reviewsNo code linkedOpen access