Astrophysical Sites Research Articles

The Relation between Escape and Scattering Times of Energetic Particles in a Turbulent Magnetized Plasma: Application to Solar Flares

Abstract A knowledge of the particle escape time from the acceleration regions of many space and astrophysical sources is of critical importance in the analysis of emission signatures produced by these particles and in the determination of the acceleration and transport mechanisms at work. This Letter addresses this general problem, in particular in solar flares, where in addition to scattering by turbulence, the magnetic field convergence from the acceleration region toward its boundaries also influences the particle escape. We test an (approximate) analytic relation between escape and scattering times, and the field convergence rate, based on the work of Malyshkin & Kulsrud, valid for both strong and weak diffusion limits and isotropic pitch-angle distributions of the injected particles, with a numerical model of particle transport. To this end, a kinetic Fokker–Planck transport model of particles is solved with a stochastic differential equation scheme, assuming different initial pitch-angle distributions. This approach enables further insights into the phase-space dynamics of the transport process, which would otherwise not be accessible. We find that in general the numerical results agree well with the analytic equation for the isotropic case; however, there are significant differences in the weak diffusion regime for non-isotopic cases, especially for distributions beamed along the magnetic field lines. The results are important in the interpretation of observations of energetic particles in solar flares and other similar space and astrophysical acceleration sites, and for the determination of acceleration-transport coefficients, commonly used in Fokker–Planck–type kinetic equations.

Large-scale Compression Acceleration during Magnetic Reconnection in a Low-β Plasma

Abstract In solar flares and other astrophysical systems, a major challenge for solving the particle acceleration problem associated with magnetic reconnection is the enormous scale separation between kinetic scales and the observed reconnection scale. Because of this, it has been difficult to draw any definite conclusions by just using kinetic simulations. A particle acceleration model that solves the energetic particle transport equation can capture the main acceleration physics found in kinetic simulations and thus provide a practical way to make observable predictions and directly compare model results with observations. Here we study compression particle acceleration in magnetic reconnection by solving the Parker (diffusion–advection) transport equation using velocity and magnetic fields from two-dimensional magnetohydrodynamics (MHD) simulations of a low-β high-Lundquist-number reconnection layer. We show that the compressible reconnection layer can give significant particle acceleration, leading to the formation of power-law particle energy distributions. We analyze the acceleration rate and find that the acceleration in the reconnection layer is a mixture of first- and second-order Fermi processes. When including a guide field, we find that the spectrum becomes steeper and both the power-law cutoff energy and maximum particle energy decrease as plasma becomes less compressible. This model produces a 2D particle distribution that one can use to generate a radiation map and directly compare with solar flare observations. This provides a framework to explain particle acceleration at large-scale astrophysical reconnection sites, such as solar flares.

Neutrino spectral split in the exact many-body formalism

We consider the many-body system of neutrinos interacting with each other through neutral current weak force. Emerging many-body effects in such a system could play important roles in some astrophysical sites such as the core collapse supernovae. In the literature this many-body system is usually treated within the mean field approximation which is an effective one-body description based on omitting entangled neutrino states. In this paper, we consider the original many-body system in an effective two flavor mixing scenario under the single angle approximation and present a solution without using the mean field approximation. Our solution is formulated around a special class of many-body eigenstates which do not undergo any level crossings as the neutrino self interaction rate decreases while the neutrinos radiate from the supernova. In particular, an initial state which consists of electron neutrinos and antineutrinos of an orthogonal flavor can be entirely decomposed in terms of those eigenstates. Assuming that the conditions are perfectly adiabatic so that the evolution of these eigenstates follow their variation with the interaction rate, we show that this initial state develops a spectral split at exactly the same energy predicted by the mean field formulation.

Stochastic Chemical Evolution of Galactic Subhalos and the Origin of r-process Elements

Abstract Mergers of compact binaries (of a neutron star and another neutron star or a black hole, NSMs) are suggested to be the promising astrophysical site of the r-process. While the average coalescence timescale of NSMs appears to be , most of previous chemical evolution models indicate that the observed early appearance and large dispersion of in Galactic halo stars at favors shorter coalescence times of 1–10 Myr. We argue that this is not the case for the models assuming the formation of the Galactic halo from clustering of subhalos with different star formation histories as suggested by Ishimaru et al. We present a stochastic chemical evolution model of the subhalos, in which the site of the r-process is assumed to be mainly NSMs with a coalescence timescale of . In view of the scarcity of NSMs, their occurrence in each subhalo is computed with a Monte Carlo method. Our results show that the less massive subhalos evolve at lower metallicities and generate highly r-process-enhanced stars. An assembly of these subhalos leaves behind the large star-to-star scatters of in the Galactic halo as observed. However, the observed scatters of [Sr/Ba] at low metallicities indicate the presence of an additional site that partially contributes to the enrichment of light neutron-capture elements such as Sr. The high enhancements of at low metallicities found in our low-mass subhalo models also qualitatively reproduce the abundance signatures of the stars in the recently discovered ultra-faint dwarf galaxy Reticulum II. Therefore, our results suggest NSMs as the dominant sources of r-process elements in the Galactic halo.

Indirect methods in nuclear astrophysics

It is well known that measuring cross-sections of thermonuclear reactions at the low energies typical of astrophysical sites is very difficult. This is due to the presence of the Coulomb barrier between the interacting nuclei. For non-explosive scenarios at astrophysical sites, the relevant energies typically span from few tens to few hundreds of keV while the Coulomb barrier is in the order of MeV. The fusion processes then proceed via tunnel effect and their cross-sections are strongly depending on the probability of penetration through the barrier. In a first approximation, this probability is given by the Gamow factor. Owing to the exponential decrease in this factor with energy, the cross-section values of thermonuclear fusion processes in stellar systems often reach values as small as micro- and nanobarn and even lower ones. Neutron-induced reactions, in spite of the absence of Coulomb barrier, are also difficult to measure. Indeed, it is the possible presence of a centrifugal barrier that can hinder the measurement of the values of the cross-sections of these processes. In either cases the cross-sections of astrophysical nuclear processes result in experimental difficulties, due to the low signal-to-noise ratio, that have been a challenge for scientists since the setting of this scientific field: Nuclear Astrophysics. In the last two to three decades, experimental improvements, including the construction of underground laboratories, allowed for the first time the measurement of cross-sections of astrophysical nuclear processes in the relevant energy region for astrophysics. Also, indirect methods were developed. As a general and common feature, using these methods it is possible to relate the features --typically the cross-section-- of a process that is experimentally simpler to measure, although not directly linked to astrophysics, to those of another process that is of interest for this latter field. This didactic paper will briefly describe some of these indirect methods with a special emphasis on Trojan Horse and on its application also to reactions that involve the use of radioactive ion beams and to neutron-induced reactions.

Probing the Early Universe through nuclear physics

Big Bang Nucleosynthesis (BBN) requires several nuclear physics inputs and nuclear reaction rates. An up-to-date compilation of direct cross sections of d(d,p)t, d(d,n)3 He and 3He(d,p)4He reactions is given, being these ones among the most uncertain bare-nucleus cross sections. An intense experimental effort has been carried on in the last decade to apply the Trojan Horse Method (THM) to study reactions of relevance for the BBN and measure their astrophysical S(E)-factor. The reaction rates and the relative error for the four reactions of interest are then numerically calculated in the temperature ranges of relevance for BBN (0.01 < T9 < 10). These value were then used as input physics for primordial nucleosynthesis calculations in order to evaluate their impact on the calculated primordial abundances and isotopical composition for H, He and Li. New results on the 7Be(n, α)4 He reaction rate were also taken into account. These were compared with the observational primordial abundance estimates in different astrophysical sites. Reactions to be studied in perspective will also be discussed

The importance of observing astrophysical tau neutrinos

The evidence of a new population of diffuse high-energy neutrinos, obtained by IceCube, has opened a new era in the field of neutrino physics. While most of the detected events are without any source counterpart, they are compatible with the standard picture of cosmic neutrinos undergoing 3-flavor neutrino oscillations. We analyze the implications of neutrino oscillations for the present and future experiments, focusing particularly on tau neutrinos. In fact tau neutrinos are very important: even if they are not produced in astrophysical sites, they have to exist due to oscillations, and their observation should be regarded as a basic proof in support of this scenario. Moreover, IceCube's measurement of the flux of muon neutrinos implies that the flux of tau neutrinos is measured within 20%, just assuming standard neutrino oscillations. On this basis, after discussing the experimental signatures of tau neutrinos, we predict the rates for ντ detection in the present and future neutrino telescopes. We show that the present IceCube detector is close to observe the first tau neutrinos, with a probability of about 90%. Moreover the next generation of IceCube can identify about 2 neutrinos per year, reaching an evidence of 5σ in about 10 years, despite the present uncertainty on the spectrum and on the production mechanism. The non observation of these neutrino events would have dramatic implications, such as the questioning of cosmic neutrino observations or the violation of neutrino oscillations over cosmological scales.

The first direct double neutron star merger detection: Implications for cosmic nucleosynthesis

Context. The astrophysical r-process site where about half of the elements, heavier than iron are produced, has been a puzzle for several decades. Here we discuss the role of one of the leading ideas – neutron star mergers (NSMs) – in the light of the first direct detection of such an event in both gravitational (GW) and electromagnetic (EM) waves. Aims. Our aim is to understand the implications of the first GW/EM observations of a NSM for cosmic nucleosynthesis. Methods. We analyse bolometric and NIR lightcurves of the first detected double NSM and compare them to nuclear reaction network-based macronova models. Results. The slope of the bolometric lightcurve is consistent with the radioactive decay of neutron star ejecta with Ye ≲ 0.3 (but not larger), which provides strong evidence for an r-process origin of the electromagnetic emission. This rules out in particular “nickel winds” as major source of the emission. We find that the NIR lightcurves can be well fitted either with or without lanthanide-rich ejecta. Our limits on the ejecta mass together with estimated rates directly confirm earlier purely theoretical or indirect observational conclusions that double neutron star mergers are indeed a major site of cosmic nucleosynthesis. If the ejecta mass was typical, NSMs can easily produce all of the estimated Galactic r-process matter, and – depending on the real rate – potentially even more. This could be a hint that the event ejected a particularly large amount of mass, maybe due to a substantial difference between the component masses. This would be compatible with the mass limits obtained from the GW-observation. Conclusions. The recent observations suggests that NSMs are responsible for a broad range of r-process nuclei and that they are at least a major, but likely the dominant r-process site in the Universe.

From Nuclei to the Cosmos: Tracing Heavy-Element Production with the Oldest Stars

Understanding the origin of the elements has been a decades-long pursuit, with many open questions remaining. Old stars found in the Milky Way and its dwarf satellite galaxies can provide answers because they preserve clean element abundance patterns of the nucleosynthesis processes that operated some 13 billion years ago, enabling reconstruction of the chemical evolution of the elements. This review focuses on the astrophysical signatures of heavy neutron-capture elements made in the s-, i-, and r-processes found in old stars. A highlight is the recently discovered r-process galaxy Reticulum II, which was enriched by a neutron star merger. These results show that old stars in dwarf galaxies provide a novel means to constrain the astrophysical site of the r-process, ushering in much-needed progress on this major outstanding question. This nuclear astrophysics research complements the many experimental and theoretical nuclear physics efforts into heavy-element formation, and also aligns with results on the gravitational-wave signature of neutron star mergers.

A Stringent Limit on the Mass Production Rate of r-process Elements in the Milky Way

Abstract We analyze data from several studies of metal-poor stars in the Milky Way, focusing individually on the main r-process elements (Eu) as well as the lighter neutron-capture element Sr, at the neutron-magic peak N = 50. Because these elements were injected in an explosion, we calculate the mass swept up when the blast wave first becomes radiative, yielding a lower limit for the dilution of such elements and hence a lower limit on the ejecta mass that is incorporated into the next generation of stars. Our study demonstrates that in order to explain the largest enhancements in [Eu/Fe] observed in stars at low [Fe/H] metallicities, individual r-process production events must synthesize a minimum of roughly 10−3 M ⊙ of r-process material. This provides a critical constraint on galactic chemical evolution models. We also show independently that if the site of Mg production is the same as that of Eu, individual injection events must synthesize up to ∼10−3 M ⊙ of r-process material. On the other hand, demanding that Sr traces Mg production results in r-process masses per event of ∼10−5 M ⊙ . This suggests that the astrophysical sites responsible for the genesis of the main r-process elements need to operate at a drastically reduced rate when compared to standard core-collapse supernovae.

The R-process Alliance: First Release from the Southern Search for R-process-enhanced Stars in the Galactic Halo*

Abstract The recent detection of a binary neutron star merger and the clear evidence of the decay of radioactive material observed in this event have, after 60 years of effort, provided an astrophysical site for the rapid neutron-capture (r-) process which is responsible for the production of the heaviest elements in our universe. However, observations of metal-poor stars with highly enhanced r-process elements have revealed abundance patterns suggesting that multiple sites may be involved. To address this issue, and to advance our understanding of the r-process, we have initiated an extensive search for bright (V &lt; 13.5), very metal-poor ([Fe/H] &lt; −2) stars in the Milky Way halo exhibiting strongly enhanced r-process signatures. This paper presents the first sample collected in the southern hemisphere using the echelle spectrograph on du Pont 2.5 m telescope at Las Campanas Observatory. We have observed and analyzed 107 stars with −3.13 &lt; [Fe/H] &lt; −0.79. Of those, 12 stars are strongly enhanced in heavy r-process elements (r-II), 42 stars show moderate enhancements of heavy r-process material (r-I), and 20 stars exhibit low abundances of the heavy r-process elements and higher abundances of the light r-process elements relative to the heavy ones (limited-r). This search is more successful at finding r-process-enhanced stars compared to previous searches, primarily due to a refined target selection procedure that focuses on red giants.

Neutron Skins and Neutron Stars in the Multimessenger Era.

The historical first detection of a binary neutron star merger by the LIGO-Virgo Collaboration [B. P. Abbott etal., Phys. Rev. Lett. 119, 161101 (2017)PRLTAO0031-900710.1103/PhysRevLett.119.161101] is providing fundamental new insights into the astrophysical site for the r process and on the nature of dense matter. A set of realistic models of the equation of state (EOS) that yield an accurate description of the properties of finite nuclei, support neutron stars of two solar masses, and provide a Lorentz covariant extrapolation to dense matter are used to confront its predictions against tidal polarizabilities extracted from the gravitational-wave data. Given the sensitivity of the gravitational-wave signal to the underlying EOS, limits on the tidal polarizability inferred from the observation translate into constraints on the neutron-star radius. Based on these constraints, models that predict a stiff symmetry energy, and thus large stellar radii, can be ruled out. Indeed, we deduce an upper limit on the radius of a 1.4M_{⊙} neutron star of R_{⋆}^{1.4}<13.76 km. Given the sensitivity of the neutron-skin thickness of ^{208}Pb to the symmetry energy, albeit at a lower density, we infer a corresponding upper limit of about R_{skin}^{208}≲0.25 fm. However, if the upcoming PREX-II experiment measures a significantly thicker skin, this may be evidence of a softening of the symmetry energy at high densities-likely indicative of a phase transition in the interior of neutron stars.

The Origin of r-process Elements in the Milky Way

Abstract Some of the heavy elements, such as gold and europium (Eu), are almost exclusively formed by the rapid neutron capture process (r-process). However, it is still unclear which astrophysical site between core-collapse supernovae and neutron star–neutron star (NS–NS) mergers produced most of the r-process elements in the universe. Galactic chemical evolution (GCE) models can test these scenarios by quantifying the frequency and yields required to reproduce the amount of europium (Eu) observed in galaxies. Although NS–NS mergers have become popular candidates, their required frequency (or rate) needs to be consistent with that obtained from gravitational wave measurements. Here, we address the first NS–NS merger detected by LIGO/Virgo (GW170817) and its associated gamma-ray burst and analyze their implication for the origin of r-process elements. The range of NS–NS merger rate densities of 320–4740 Gpc−3 yr−1 provided by LIGO/Virgo is remarkably consistent with the range required by GCE to explain the Eu abundances in the Milky Way with NS–NS mergers, assuming the solar r-process abundance pattern for the ejecta. Under the same assumption, this event has produced about 1–5 Earth masses of Eu, and 3–13 Earth masses of gold. When using theoretical calculations to derive Eu yields, constraining the role of NS–NS mergers becomes more challenging because of nuclear astrophysics uncertainties. This is the first study that directly combines nuclear physics uncertainties with GCE calculations. If GW170817 is a representative event, NS–NS mergers can produce Eu in sufficient amounts and are likely to be the main r-process site.

The Roles of Fluid Compression and Shear in Electron Energization during Magnetic Reconnection

Abstract Particle acceleration in space and astrophysical reconnection sites is an important unsolved problem in studies of magnetic reconnection. Earlier kinetic simulations have identified several acceleration mechanisms that are associated with particle drift motions. Here, we show that, for sufficiently large systems, the energization processes due to particle drift motions can be described as fluid compression and shear, and that the shear energization is proportional to the pressure anisotropy of energetic particles. By analyzing results from fully kinetic simulations, we show that the compression energization dominates the acceleration of high-energy particles in reconnection with a weak guide field, and the compression and shear effects are comparable when the guide field is 50% of the reconnecting component. Spatial distributions of those energization effects reveal that reconnection exhausts, contracting islands, and island-merging regions are the three most important regions for compression and shear acceleration. This study connects particle energization by particle guiding-center drift motions with that due to background fluid motions, as in the energetic particle transport theory. It provides foundations for building particle transport models for large-scale reconnection acceleration such as those in solar flares.

Enrichment of Zinc in Galactic Chemodynamical Evolution Models

Abstract The heaviest iron-peak element Zinc (Zn) has been used as an important tracer of cosmic chemical evolution. Spectroscopic observations of the metal-poor stars in Local Group galaxies show an increasing trend of [Zn/Fe] ratios toward lower metallicity. However, the enrichment of Zn in galaxies is not well understood due to poor knowledge of astrophysical sites of Zn, as well as metal mixing in galaxies. Here we show possible explanations for the observed trend by taking into account electron-capture supernovae (ECSNe) as one of the sources of Zn in our chemodynamical simulations of dwarf galaxies. We find that the ejecta from ECSNe contribute to stars with [Zn/Fe] ≳ 0.5. We also find that scatters of [Zn/Fe] in higher metallicities originate from the ejecta of type Ia supernovae. On the other hand, it appears difficult to explain the observed trends if we do not consider ECSNe as a source of Zn. These results come from an inhomogeneous spatial metallicity distribution due to the inefficiency of the metal mixing. We find that the optimal value of the scaling factor for the metal diffusion coefficient is ∼0.01 in the shear-based metal mixing model in smoothed particle hydrodynamics simulations. These results suggest that ECSNe could be one of the contributors of the enrichment of Zn in galaxies.

Trojan Horse cross section measurements and their impact on primordial nucleosynthesis

Big Bang Nucleosynthesis (BBN) nucleosynthesis requires several nuclear physics inputs and, among them, an important role is played by nuclear reaction rates. They are among the most important input for a quantitative description of the early Universe. An up-to-date compilation of direct cross sections of d(d,p)t, d(d,n)3He and 3He(d,p)4He reactions is given, being these ones among the most uncertain bare-nucleus cross sections.An intense experimental effort has been carried on in the last decade to apply the Trojan Horse Method (THM) to study reactions of relevance for the BBN and measure their astrophysical S(E)-factor. The result of these recent measurements is reviewed and compared with the available direct data. The reaction rates and the relative error for the four reactions of interest are then numerically calculated in the temperature ranges of relevance for BBN (0.01<T9 <10) and compared with up-to-date reaction rate compilations. Their value were therefore used as input physics for primordial nucleosynthesis calculations in order to evaluate their impact on the calculated primordial abundances of D, 3,4He and 7Li. These ones were then compared with the observational primordial abundance estimates in different astrophysical sites. A comparison was also performed with calculations using other reaction rates compilations available in literature.

Direct measurements and detection techniques with low-energy RIBs

Astrophysical reactions involving radioactive isotopes (RI) are of importance for the stellar energy generation and nucleosynthesis especially in high-temperature astrophysical sites, such as X-ray bursts, core-collapse supernovae, and supermassive metalpoor stars. In spite of the essential diffculties in the experimental evaluation of those reaction rates, there are several successful approaches to study them, owing to the recent technical developments in the beam production, measurement method, and detectors. Among them, the measurements of α resonant scattering and (α, p) reactions using the thick-target method in inverse kinematics are discussed. The experiments at the low-energy RI beam separator CRIB, operated by Center for Nuclear Study (CNS), the University of Tokyo, are introduced as examples for such studies.

Metallicity fluctuation statistics in the interstellar medium and young stars – I. Variance and correlation

The distributions of a galaxy's gas and stars in chemical space encodes a tremendous amount of information about that galaxy's physical properties and assembly history. However, present methods for extracting information from chemical distributions are based either on coarse averages measured over galactic scales (e.g., metallicity gradients) or on searching for clusters in chemical space that can be identified with individual star clusters or gas clouds on $\sim 1$ pc scales. These approaches discard most of the information, because in galaxies gas and young stars are observed to be distributed fractally, with correlations on all scales, and the same is likely to be true of metals. In this paper we introduce a first theoretical model, based on stochastically-forced diffusion, capable of predicting the multi-scale statistics of metal fields. We derive the variance, correlation function, and power spectrum of the metal distribution from first principles, and determine how these quantities depend on elements' astrophysical origin sites and on the large-scale properties of galaxies. Among other results, we explain for the first time why the typical abundance scatter observed in the interstellar media of nearby galaxies is $\approx 0.1$ dex, and we predict that this scatter will be correlated on spatial scales of $\sim 0.5-1$ kpc, and over time scales of $\sim 100-300$ Myr. We discuss the implications of our results for future chemical tagging studies.

R -process nucleosynthesis from matter ejected in binary neutron star mergers

When binary systems of neutron stars merge, a very small fraction of their rest mass is ejected, either dynamically or secularly. This material is neutron-rich and its nucleosynthesis could provide the astrophysical site for the production of heavy elements in the universe, together with a kilonova signal confirming neutron-star mergers as the origin of short gamma-ray bursts. We perform full general-relativistic simulations of binary neutron-star mergers employing three different nuclear-physics EOSs, considering both equal- and unequal-mass configurations, and adopting a leakage scheme to account for neutrino radiative losses. Using a combination of techniques, we carry out an extensive and systematic study of the hydrodynamical, thermodynamical, and geometrical properties of the matter ejected dynamically, employing the WinNet nuclear-reaction network to recover the relative abundances of heavy elements produced by each configurations. Among the results obtained, three are particularly important. First, we find that both the properties of the dynamical ejecta and the nucleosynthesis yields are robust against variations of the EOS and masses, and match very well the observed chemical abundances. Second, using a conservative but robust criterion for unbound matter, we find that the amount of ejected mass is $\lesssim 10^{-3}\,M_{\odot}$, hence at least one order of magnitude smaller than what normally assumed in modelling kilonova signals. Finally, using a simplified and gray-opacity model we assess the observability of the infrared kilonova emission finding, that for all binaries the luminosity peaks around $\sim1/2$ day in the $H$-band, reaching a maximum magnitude of $-13$, and decreasing rapidly after one day. These rather low luminosities make the prospects for detecting kilonovae less promising than what assumed so far.

First measurement of the S34(p,γ)Cl35 reaction rate through indirect methods for presolar nova grains

Sulphur isotopic ratio measurements may help to establish the astrophysical sites in which certain presolar grains were formed. Nova model predictions of the S-34/S-32 ratio are, however, unreliable due to the lack of an experimental S-34(p,gamma)Cl-35 reaction rate. To this end, we have measured the S-34(He-3,d)Cl-35 reaction at 20 MeV using a high resolution quadrupole-dipole-dipole-dipolemagnetic spectrograph. Twenty-two levels over 6.2 MeV < E-x (Cl-35) < 7.4 MeV were identified, ten of which were previously unobserved. Proton-transfer spectroscopic factors have been measured for the first time over the energy range relevant for novae. With this new spectroscopic information a new S-34(p,gamma)Cl-35 reaction rate has been determined using a Monte Carlo method. Hydrodynamic nova model calculations have been performed using this new reaction rate. These models show that remaining uncertainties in the S-34(p,gamma) rate affect nucleosynthesis predictions by less than a factor of 1.4, and predict a S-34/S-32 isotopic ratio of 0.014-0.017. Since recent type II supernova models predict S-34/S-32 = 0.026-0.053, the S-34/S-32 isotopic ratio may be used, in conjunction with other isotopic signatures, to distinguish presolar grains from oxygen-neon nova and type II supernova origin. Our results address a key nuclear physics uncertainty on which recent considerations discounting the nova origin of several grains depend.