skip to main content

Title: Nuclear reaction network unveils novel reaction patterns based on stellar energies
Abstract Despite the advances in discovering new nuclei, modeling microscopic nuclear structure, nuclear reactors, and stellar nucleosynthesis, we still lack a systemic tool, such as a network approach, to understand the structure and dynamics of over 70 thousands reactions compiled in JINA REACLIB. To this end, we develop an analysis framework, under which it is simple to know which reactions generally are possible and which are not, by counting neutrons and protons incoming to and outgoing from any target nucleus. Specifically, we assemble here a nuclear reaction network in which a node represents a nuclide, and a link represents a direct reaction between nuclides. Interestingly, the degree distribution of nuclear network exhibits a bimodal distribution that significantly deviates from the common power-law distribution of scale-free networks and Poisson distribution of random networks. Based on the dynamics from the cross section parameterizations in REACLIB, we surprisingly find that the distribution is universal for reactions with a rate below the threshold, λ < e − T γ , where T is the temperature and γ ≈ 1.05. Moreover, we discover three rules that govern the structure pattern of nuclear reaction network: (i) reaction-type is determined by linking choices, (ii) network distances between the reacting nuclides on 2D grid of Z vs N of nuclides are short, and (iii) each node in- and out-degrees are close to each other. By incorporating these three rules, our model universally unveils the underlying nuclear reaction patterns hidden in a large and dense nuclear reaction network regardless of nuclide chart expansions. It enables us to predict missing links that represent possible new nuclear reactions not yet discovered.  more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
New Journal of Physics
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract During the X-ray bursts of GS 1826−24, a “clocked burster”, the nuclear reaction flow that surges through the rapid-proton capture process path has to pass through the NiCu cycles before reaching the ZnGa cycles that moderate further hydrogen burning in the region above the germanium and selenium isotopes. The 57 Cu(p, γ ) 58 Zn reaction that occurs in the NiCu cycles plays an important role in influencing the burst light curves found by Cyburt et al. We deduce the 57 Cu(p, γ ) 58 Zn reaction rate based on the experimentally determined important nuclear structure information, isobaric-multiplet-mass equation, and large-scale shell-model calculations. Based on the isobaric-multiplet-mass equation, we propose a possible order of 1 1 + - and 2 3 + -dominant resonance states and constrain the resonance energy of the 1 2 + state. The latter reduces the contribution of the 1 2 + -dominant resonance state. The new reaction rate is up to a factor of 4 lower than the Forstner et al. rate recommended by JINA REACLIB v2.2 at the temperature regime sensitive to clocked bursts of GS 1826−24. Using the simulation from the one-dimensional implicit hydrodynamic code K epler to model the thermonuclear X-ray bursts of the GS 1826−24 clocked burster, we find that the new 57 Cu(p, γ ) 58 Zn reaction rate, coupled with the latest 56 Ni(p, γ ) 57 Cu and 55 Ni(p, γ ) 56 Cu reaction rates, redistributes the reaction flow in the NiCu cycles and strongly influences the burst ash composition, whereas the 59 Cu(p, α ) 56 Ni and 59 Cu(p, γ ) 60 Zn reactions suppress the influence of the 57 Cu(p, γ ) 58 Zn reaction and diminish the impact of nuclear reaction flow that bypasses the important 56 Ni waiting point induced by the 55 Ni(p, γ ) 56 Cu reaction on the burst light curve. 
    more » « less
  2. An accurate description of the evolution of dislocation networks is an essential part of discrete and continuum dislocation dynamics models. These networks evolve by motion of the dislocation lines and by forming junctions between these lines via cross slip, annihilation and junction reactions. In this work, we introduce these dislocation reactions into continuum dislocation models using the theory of de Rham currents. We introduce dislocations on each slip system as potentially open lines whose boundaries are associated with junction points and, therefore, still create a network of collectively closed lines that satisfy the classical relations and for the dislocation density tensor and the plastic distortion . To ensure this, we leverage Frank’s second rule at the junction nodes and the concept of virtual dislocation segments. We introduce the junction point density as a new state variable that represents the distribution of junction points within the crystal containing the dislocation network. Adding this information requires knowledge of the global structure of the dislocation network, which we obtain from its representation as a graph. We derive transport relations for the dislocation line density on each slip system in the crystal, which now includes a term that corresponds to the motion of junction points. We also derive the transport relations for junction points, which include source terms that reflect the topology changes of the dislocation network due to junction formation. 
    more » « less
  3. Liu, W. ; Wang, Y. ; Guo, B. ; Tang, X. ; Zeng, S. (Ed.)
    In Type-I X-ray bursts (XRBs), the rapid-proton capture (rp-) process passes through the NiCu and ZnGa cycles before reaching the region above Ge and Se isotopes that hydrogen burning actively powers the XRBs. The sensitivity study performed by Cyburt et al . [1] shows that the 57 Cu(p, γ ) 58 Zn reaction in the NiCu cycles is the fifth most important rp-reaction influencing the burst light curves. Langer et al . [2] precisely measured some low-lying energy levels of 58 Zn to deduce the 57 Cu(p, γ ) 58 Zn reaction rate. Nevertheless, the order of the 1 + 1 and 2 + 3 resonance states that dominate at 0:2 ≲ T (GK) ≲ 0:8 is not confirmed. The 1 + 2 resonance state, which dominates at the XRB sensitive temperature regime 0:8 ≲ T (GK) ≲ 2 was not detected. Using isobaric-multipletmass equation (IMME), we estimate the order of the 1 + 1 and 2 + 3 resonance states and estimate the lower limit of the 1 + 2 resonance energy. We then determine the 57 Cu(p, γ ) 58 Zn reaction rate using the full pf -model space shell model calculations. The new rate is up to a factor of four lower than the Forstner et al . [3] rate recommended by JINA REACLIBv2.2. Using the present 57 Cu(p, γ ) 58 Zn, the latest 56 Ni(p, γ ) 57 Cu and 55 Ni(p, γ ) 56 Cu reaction rates, and 1D implicit hydrodynamic K epler code, we model the thermonuclear XRBs of the clocked burster GS 1826–24. We find that the new rates regulate the reaction flow in the NiCu cycles and strongly influence the burst-ash composition. The 59 Cu(p, γ ) 56 Ni and 59 Cu(p, α ) 60 Zn reactions suppress the influence of the 57 Cu(p, γ ) 58 Zn reaction. They strongly diminish the impact of the nuclear reaction flow that bypasses the 56 Ni waiting point induced by the 55 Ni(p, γ ) 56 Cu reaction on burst light curve. 
    more » « less
  4. Context. The γ -process nucleosynthesis in core-collapse supernovae is generally accepted as a feasible process for the synthesis of neutron-deficient isotopes beyond iron. However, crucial discrepancies between theory and observations still exist: the average yields of γ -process nucleosynthesis from massive stars are still insufficient to reproduce the solar distribution in galactic chemical evolution calculations, and the yields of the Mo and Ru isotopes are a factor of ten lower than the yields of the other γ -process nuclei. Aims. We investigate the γ -process in five sets of core-collapse supernova models published in the literature with initial masses of 15, 20, and 25 M ⊙ at solar metallicity. Methods. We compared the γ -process overproduction factors from the different models. To highlight the possible effect of nuclear physics input, we also considered 23 ratios of two isotopes close to each other in mass relative to their solar values. Further, we investigated the contribution of C–O shell mergers in the supernova progenitors as an additional site of the γ -process. Results. Our analysis shows that a large scatter among the different models exists for both the γ -process integrated yields and the isotopic ratios. We find only ten ratios that agree with their solar values, all the others differ by at least a factor of three from the solar values in all the considered sets of models. The γ -process within C–O shell mergers mostly influences the isotopic ratios that involve intermediate and heavy proton-rich isotopes with A  > 100. Conclusions. We conclude that there are large discrepancies both among the different data sets and between the model predictions and the solar abundance distribution. More calculations are needed; particularly updating the nuclear network, because the majority of the models considered in this work do not use the latest reaction rates for the γ -process nucleosynthesis. Moreover, the role of C–O shell mergers requires further investigation. 
    more » « less
  5. Signaling and metabolic pathways, which consist of a series of reactions producing target molecules from source compounds, are cornerstones of cellular biology. The cellular reaction networks containing such pathways can be precisely modeled by directed hypergraphs, where each reaction corresponds to a hyperedge, directed from its set of reactants to its set of products. Given such a network represented by a directed hypergraph, inferring the most likely set of reactions that produce a given target from a given set of sources corresponds to finding a shortest hyperpath, which is NP-complete. The best methods currently available for shortest hyperpaths either offer no guarantee of optimality, or exclude hyperpaths containing cycles even though cycles are abundant in real biological pathways. We derive a novel graph-theoretic characterization of hyperpaths, leveraged in a new formulation of the general shortest hyperpath problem as an integer linear program that for the first time handles hyperpaths containing cycles, and present a novel cutting-plane algorithm that can solve this integer program to optimality in practice. This represents a major advance over the best prior exact algorithm, which was limited to acyclic hyperpaths (and hence fails to find a solution for the many biological instances where all hyperpaths are in fact cyclic). In comprehensive experiments over thousands of instances from the standard NCI-PID and Reactome databases, we demonstrate that our cutting-plane algorithm quickly finds an optimal hyperpath, with a median running-time of under ten seconds and a maximum time of around thirty minutes, even on large instances with many thousands of reactions. Source code implementing our cutting-plane algorithm for shortest hyperpaths in a new tool called Mmunin is available free for research use at 
    more » « less