Search

The origin of the proton-rich trans-iron isotopes in the solar system is still uncertain. Single-degenerate thermonuclear supernovae (SNIa) with n-capture nucleosynthesis seeds assembled in the external layers of the progenitor's rapidly accreting white dwarf phase may produce these isotopes. We calculate the stellar structure of the accretion phase of five white dwarf models with initial masses >~ 0.85Msun using the stellar code MESA. The near-surface layers of the 1, 1.26, 1.32 and 1.38Msun models are most representative of the regions in which the bulk of the p nuclei are produced during SNIa explosions, and for these models we also calculate the neutron-capture nucleosynthesis in the external layers. Contrary to previous rapidly-accreting white dwarf models at lower mass, we find that the H-shell ashes are the main site of n-capture nucleosynthesis. We find high neutron densities up to several 10^15 cm^-3 in the most massive WDs. Through the recurrence of the H-shell ashes these intermediate neutron densities can be sustained effectively for a long time leading to high neutron exposures with a strong production up to Pb. Both the neutron density and the neutron exposure increase with increasing the mass of the accreting WD. Finally, the SNIa nucleosynthesis is calculated using the obtained abundances as seeds. We obtain solar to super-solar abundances for p-nuclei with A>96. Our models show that SNIa are a viable p-process production site., Accepted for publication in MNRAS. 49 pages, 17 figures, 5 tables

Atmospheric parameters and chemical compositions for ten stars with metallicities in the region of -2.2< [Fe/H], Comment: 25 pages, 7 figures

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, Comment: 17 pages, 5 figures, revised version, accepted for publication in MNRAS

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., 13 pages, 11 figures, accepted to MNRAS

The direct impact of white dwarfs has been suggested as a plausible channel for type Ia supernovae. In spite of their (a priori) rareness, in highly populated globular clusters and in galactic centers, where the amount of white dwarfs is considerable, the rate of violent collisions between two of them might be non-negligible. Even more, there are indications that binary white dwarf systems orbited by a third stellar-mass body have an important chance to induce a clean head-on collision. Therefore, this scenario represents a source of contamination for the supernova light-curves sample that it is used as standard candles in cosmology, and it deserves further investigation. Some groups have conducted numerical simulations of this scenario, but their results show several differences. In this paper we address some of the possible sources of these differences, presenting the results of high resolution hydrodynamical simulations jointly with a detailed nuclear post-processing of the nuclear abundances, to check the viability of white dwarf collisions to produce significant amounts of 56Ni. To that purpose, we use a 2D-axial symmetric smoothed particle hydrodynamic code to obtain a resolution considerably higher than in previous studies. In this work, we also study how the initial mass and nuclear composition affect the results. The gravitational wave emission is also calculated, as this is a unique signature of this kind of events. All calculated models produce a significant amount of 56Ni, ranging from 0.1 Msun to 1.1 Msun, compatible not only with normal-Branch type Ia supernova but also with the subluminous and super-Chandrasekhar subset. Nevertheless, the distribution mass-function of white dwarfs favors collisions among 0.6-0.7 Msun objects, leading to subluminous events., Comment: 24 pages, 12 figures, accepted for publication in MNRAS

Internal flows inside gravitationally stable astrophysical objects, such as the Sun, normal and compact stars, are rotating, highly compressed and extremely subsonic. Such low Mach number flows are usually encountered when studying, for example, the dynamo action in stars and planets or the nuclear burst on neutron stars and white dwarfs. Handling of such flows numerically on time-scales longer than the dynamical one is complicated and challenging. The aim of this paper is to address the numerical problems associated with the modelling of internal quasi-stationary, rotating low Mach number flows in stars and to discuss possible solution scenarios. It is shown that the quasi-symmetric approximate factorization method (AFM) as a pre-conditioner within a non-linear Newton-type defect-correction solution procedure is best suited for modelling quasi-stationary weakly compressible flows with moderate low Mach numbers. This method is robust as it can be applied to model time-dependent compressible flows without further modifications. The AFM-pre-conditioning techniques are shown to be extendable into three dimensions with an arbitrary equation of state. Classical dimensional splitting techniques, however, such as the alternating direction implicit or line-Gauss–Seidel methods are not suited for modelling compressible low Mach number flows. It is also argued that hot and low Mach number astrophysical flows cannot be considered as an asymptotic limit of incompressible flows, but rather as highly compressed flows with extremely stiff pressure terms. We show that, unlike the pseudo-pressure in incompressible fluids, a Poisson-like treatment for the pressure would smooth unnecessarily physically induced acoustic perturbations, thereby violating the conservation of the total energy. Results of several hydrodynamical calculations are presented, which demonstrate the capability of the solver to search for solutions, that correspond to stationary, viscous and rotating flows with a Mach number as small as as well as to fluid flows that are subject to ultra-strong Newtonian and general relativistic gravitational fields.

We consider the influence of potential quark matter existing at high densities in neutron star interiors on gravitational waves (GW) emitted in a binary neutron star merger event. Two types of equations of state (EoS) at zero temperatures are used, one describing pure nuclear matter, the other nuclear matter with a phase transition to quark matter at very high densities. Binary equilibrium sequences close to the innermost stable circular orbit (ISCO) are calculated to determine the GW frequencies just before merger. It is found that EoS effects begin to play a role for gravitational masses larger than $M_\infty\simeq1.5M_\odot$. The difference in the gravitational wave frequency at the ISCO grows to up to $\simeq 10%$ for the maximal allowed mass given by the EoSs used. Then, we perform 3D hydrodynamic simulations for each EoS varying the initial mass and determine the characteristic GW frequencies of the merger remnants. The implications of quark matter show up mainly in a different collapse behaviour of the merger remnant. If the collapse does not take place immediately after merger, we find a phase difference between two EoS's in the post-merger GW signal. We also compare the GW frequencies emitted by the merger remnant to values from simulations using a polytropic EoS and find an imprint of the non-constant adiabatic index of our EoSs. All calculations are based on the conformally flat (CF) approximation to general relativity and the GW signal from the merger simulation is extracted up to quadrupole order., 13 pages, 8 figures