Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nuclear Reactions”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

Systematic exploration of heavy element nucleosynthesis in protomagnetar outflows

ABSTRACT We study the nucleosynthesis products in neutrino-driven winds from rapidly rotating, highly magnetized and misaligned protomagnetars using the nuclear reaction network SkyNet. We adopt a semi-analytic parametrized model for the protomagnetar and systematically study the capabilities of its neutrino-driven wind for synthesizing nuclei and eventually producing ultra-high energy cosmic rays (UHECRs). We find that for neutron-rich outflows (Ye < 0.5), synthesis of heavy elements ($\overline{A}\sim 20-65$) is possible during the first $\sim 10\, {\rm s}$ of the outflow, but these nuclei are subjected to composition-altering photodisintegration during the epoch of particle acceleration at the dissipation radii. However, after the first $\sim 10\, {\rm s}$ of the outflow, nucleosynthesis reaches lighter elements ($\overline{A}\sim 10-50$) that are not subjected to subsequent photodisintegration. For proton-rich (Ye ≥ 0.5) outflows, synthesis is more limited ($\overline{A}\sim 4-15$). These suggest that while protomagnetars typically do not synthesize nuclei heavier than second r-process peak elements, they are intriguing sources of intermediate/heavy mass UHECRs. For all configurations, the most rapidly rotating protomagnetars are more conducive for nucleosynthesis with a weaker dependence on the magnetic field strength.

79 ASTRONOMY AND ASTROPHYSICS↗

Radioisotope production using lasers: From basic science to applications

The discovery of chirped pulse amplification has led to great improvements in laser technology, enabling energetic laser beams to be compressed to pulse durations of tens of femtoseconds and focused to a few micrometers. Protons with energies of tens of MeV can be accelerated using, for instance, target normal sheath acceleration and focused on secondary targets. Under such conditions, nuclear reactions can occur, with the production of radioisotopes suitable for medical application. The use of high-repetition lasers to produce such isotopes is competitive with conventional methods mostly based on accelerators. In this paper, we study the production of 67Cu, 63Zn, 18F, and 11C, which are currently used in positron emission tomography and other applications. At the same time, we study the reactions 10B(p,α)7Be and 70Zn(p,4n)67Ga to put further constraints on the proton distributions at different angles, as well as the reaction 11B(p,α)8Be relevant for energy production. The experiment was performed at the 1 PW laser facility at Vega III in Salamanca, Spain. Angular distributions of radioisotopes in the forward (with respect to the laser direction) and backward directions were measured using a high purity germanium detector. Our results are in reasonable agreement with numerical estimates obtained following the approach of Kimura and Bonasera [Nucl. Instrum. Methods Phys. Res., Sect. A 637, 164–170 (2011)].

11C↗

An exploratory study on helium mobility in amorphous and crystallized bulk metallic glasses

In this study, Zr 52.5 Cu 17.9 Ni 14.6 Al 10 Ti 5 and Cu 60 Zr 20 Hf 10 Ti 10 amorphous and crystallized bulk metallic glasses (BMGs) were implanted at room temperature by 150 keV 3 He + ions to fluences of 2 × 10 15 cm -2 and 5 × 10 15 cm -2 , producing peak implanted He concentrations of ~1,100 and 3,500 appm, respectively. Nuclear reaction analysis (NRA) performed on post-implantation annealed samples at various temperatures between 250 and 600 °C (~0.30 to 0.75 of the absolute melting temperature T m of the BMGs) revealed no appreciable He migration. Complementary thermal desorption spectroscopy (TDS) on as-implanted samples, which involved in situ annealing, showed appreciable He release only at temperatures higher than ~700 °C (~0.76-0.90 T m ) in all the samples. Half or more of the implanted He did not desorb up to the maximum investigated TDS annealing temperature of 770 °C (~0.85-0.97 T m ). The lack of pronounced He diffusion up to 600 °C as observed from the NRA and TDS data was attributed to He trapping in free volume sites present in the BMGs, which seems to play a similar role to vacancies in crystalline materials. It was also observed that in both the amorphous and crystallized forms of the Zr BMG, the sample implanted to the lower fluence released more He as compared to the sample implanted with a higher fluence. The crystalline forms of both the Cu and Zr BMGs released more He as compared to their amorphous counterparts.

36 MATERIALS SCIENCE↗

Gamma-ray signatures of 𝑟-process radioactivity from the collapse of magnetized white dwarfs

We predict the gamma-ray line emission from 𝑟-process nuclei synthesized in the ejecta of the accretion-induced collapse (AIC) of a magnetized, rapidly rotating white dwarf. Using ejecta from a two-dimensional general-relativistic neutrino-magnetohydrodynamic simulation, further evolved with a radiation-hydrodynamics code coupled to an in situ nuclear reaction network, we construct angle-dependent gamma-ray spectra in the 0.01–10 MeV band via composition-dependent ray tracing through the ejecta. The emission between ∼1 and 10 d is dominated by 132 I (𝑡 1/2 = 2.3 h), continuously replenished by the decay of its parent 132 Te (𝑡 1/2 = 3.2 d), with additional contributions from 131 I, 133 Xe, and 132 Te . At 𝑡 ≳ 20 d, 56 Co (from 56 Ni decay) becomes the primary emitter. The simultaneous presence of 𝑟 process and iron-peak gamma-ray lines is distinctive of AIC ejecta and absent in binary neutron star mergers, where iron-peak nuclei are generally not synthesized. Comparing with the 3⁢𝜎 continuum sensitivities of planned MeV gamma-ray telescopes (COSI, AMEGO-X, e-ASTROGAM, GRAMS, GammaTPC), we find the brightest 𝑟-process lines detectable to ∼10 Mpc by GammaTPC and GRAMS, with the signal approaching their sensitivity threshold at 30 Mpc. As a result, the 𝑟-process spectral features survive time integration over ∼30 d exposures, demonstrating robustness against the long observation times required by gamma-ray detectors.

Nuclear reactions↗

Deep learning for NLTE spectral opacities

Computer simulations of high energy density science experiments are computationally challenging, consisting of multiple physics calculations including radiation transport, hydrodynamics, atomic physics, nuclear reactions, laser–plasma interactions, and more. To simulate inertial confinement fusion (ICF) experiments at high fidelity, each of these physics calculations should be as detailed as possible. However, this quickly becomes too computationally expensive even for modern supercomputers, and thus many simplifying assumptions are made to reduce the required computational time. Much of the research has focused on acceleration techniques for the various packages in multiphysics codes. In this work, we explore a novel method for accelerating physics packages via machine learning. The non-local thermodynamic equilibrium (NLTE) package is one of the most expensive calculations in the simulations of indirect drive inertial confinement fusion, taking several tens of percent of the total wall clock time. We explore the use of machine learning to accelerate this package, by essentially replacing the physics calculation with a deep neural network that has been trained to emulate the physics code. Overall, we demonstrate the feasibility of this approach on a simple problem and perform a side-by-side comparison of the physics calculation and the neural network inline in an ICF Hohlraum simulation. We show that the neural network achieves a 10× speed up in NLTE computational time while achieving good agreement with the physics code for several quantities of interest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Classification of events from α -induced reactions in the MUSIC detector via statistical and ML methods

The Multi-Sampling Ionization Chamber (MUSIC) detector is typically used to measure nuclear reaction cross sections relevant for nuclear astrophysics, fusion studies, and other applications. From the MUSIC data produced in one experiment scientists carefully extract an order of 10 3 events of interest from about 10 9 total events, where each event can be represented by an 18-dimensional vector. However, the standard data classification process is based on expert driven, manually intensive data analysis techniques that require several months to identify patterns and classify the relevant events from the collected data. Here, to address this issue, we present a method for the classification of events originating from specific α-induced reactions by combining statistical and machine learning methods that require significantly less input from the domain scientist, relative to the standard technique. Here, we applied the new method to two experimental data sets and compared our results with those obtained using traditional methods. With few exceptions, the number of events classified by our method agrees within ±20% with the results obtained using traditional methods. With the present method, which is the first of its kind for the MUSIC data, we have established the foundation for the automated extraction of physical events of interest from experiments using the MUSIC detector.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Green's function knockout formalism

Knockout nuclear reactions, in which a nucleon is removed from a nucleus as a result of the collision with another nucleus, have been widely used as an experimental tool, both to populate isotopes further removed from stability and to obtain information about the single-particle nature of the nuclear spectrum. In order to fully exploit the experimental information, theory is needed for the description of both the structure of the nuclei involved, and the dynamics associated with the nucleon removal mechanisms. Here, the standard approach, using theoretical shell-model spectroscopic factors for the structure description coupled with an eikonal model of reaction, has been successful when used in the context of the removal of valence nucleons in nuclei close to stability. However, it has been argued that the reaction theory might need to be revisited in the case of exotic nuclei, more specifically for highly asymmetric nuclei in which the deficient species (neutrons or protons) is being removed. We present here a new formalism for the nucleon-removal and -addition reaction through knockout and transfer reactions, that treats consistently structure and reaction properties using dispersive optical potentials. In particular, our formalism includes the dynamical effects associated with the removal of a nucleon from the projectile, which might explain the long standing puzzle of the quenching of spectroscopic factors in nuclei with extreme neutrons-to-protons ratios.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Bayesian quantification of observability and equation of state of twin stars

The possibility of discovering twin stars, two neutron stars (NSs) with the same mass but different radii, is usually studied in forward modelings by using a restricted number of NS matter equations of state (EOSs) encapsulating a first-order phase transition from hadronic to quark matter (QM). Informing our likelihood function with the NS radius data from GW170817 and using a metamodel with nine parameters capable of mimicking most NS EOSs available in the literature, we conduct a Bayesian quantification of the observability and underlying EOSs of twin stars. Of the accepted EOSs, between 12 and 18% yield twin stars, depending on the restrictions we place on the second branch. The possibility of twin stars remains robust even under recent observational constraints. Here, we show that many of these twin star scenarios are observable with currently available levels of accuracy in measuring NS radii. We also present the marginalized posterior probability density functions (PDFs) of every EOS parameter for each of four mass-radius correlation topologies. We find that the inferred EOS depends sensitively on not only whether twin stars are present, but also the category of twin stars, indicating that the observation of twin stars would provide a strong constraint on the underlying EOS. In particular, for two coexisting hybrid stars having QM cores at different densities, the PDF for QM speed of sound squared 𝑐$^2_{qm}$ has two peaks, one below and another above the conformal limit 𝑐$^2_{qm}$ = 1/3 predicted by perturbative QCD.

QCD in nuclear reactions↗

Extension of the Hauser-Feshbach fission fragment decay model to multichance fission

The Hauser-Feshbach fission fragment decay model, HF 3 D, which calculates the statistical decay of fission fragments, has been expanded to include multichance fission, up to neutron incident energies of 20 MeV. The deterministic decay takes as input prescission quantities—fission probabilities and the average energy causing fission—and postscission quantities—yields in mass, charge, total kinetic energy, spin, and parity. From these fission fragment initial conditions, the full decay is followed through both prompt and delayed particle emissions, allowing for the calculation of prompt neutron and γ properties, such as multiplicity and energy distributions, both independent and cumulative fission yields, and delayed neutron observables. Here, we describe the implementation of multichance fission in the HF 3 D model, and show an example of prompt and delayed quantities beyond first-chance fission, using the example of neutron-induced fission on 235 U. This expansion represents significant progress in consistently modeling the emission of prompt and delayed particles from fissile systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Measurement of the Ar ( e , e ' p ) and Ti ( e , e ' p ) cross sections in Jefferson Lab Hall A

The E12-14-012 experiment, performed in Jefferson Lab Hall A, has collected exclusive electron-scattering data ( e , e ' p ) in parallel kinematics using natural argon and natural titanium targets. In this study we report the first results of the analysis of the data set corresponding to beam energy 2222 GeV, electron scattering angle 21 . 5 ° , and proton emission angle – 50 ° . The differential cross sections, measured with ≈ 4 % uncertainty, have been studied as a function of missing energy and missing momentum, and compared to the results of Monte Carlo simulations, obtained from a model based on the distorted-wave impulse approximation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Angular correlation between the fission fragment intrinsic spins

The generation of fission fragment (FF) spins and of their relative orbital angular momentum has been debated for more than six decades, and no consensus has been achieved so far. The interpretation of recent experimental results of Wilson et al., has been challenged by several recent theoretical studies, which are not in agreement with one another. According to the interpretation by Wilson et al., the FF spins emerge very long after scission occurrs. Randrup and Vogt, while agreeing that the FF spins are uncorrelated, conclude on the basis of a phenomenological model that these spins are uncorrelated already before scission. Bulgac et al., in a fully microscopic study demonstrate that the primordial FF spin final values are defined before the emission of prompt neutrons and statistical gammas and are strongly correlated with a relative angle between spins close to 2π/3, a result in full agreement with the present independent analysis. The prompt neutrons and statistical gammas carry a significant amount of angular momentum according to the study of Stetcu et al. which can lead to a decorrelation of the FF bandhead spins of the yrast lines measured by Wilson et al., and which also provides arguments why these measured spins are so different from the primordial FF spins evaluated by Bulgac et al.. Here, I show that the unexpected character of the angular correlation between the primordial FF intrinsic spins, recently evaluated by Bulgac et al., which favors FF intrinsic spins pointing predominantly in opposite directions, can be understood by using simple general phase space arguments. The observation by Wilson et al., that the FF spins are uncorrelated follows from both the results of the microscopic calculation of Bulgac et al., and the present analysis of the full correlated probability distribution of the FF spins together with the relative orbital angular momentum of the primordial FFs. These arguments may apply also to heavy-ion collisions and, since there is no use of specifics of the particle interactions, the present results might apply to atomic and molecular systems as well.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Modeling fission dynamics at the barrier in a discrete-basis formalism

A configuration-interaction model is presented for the barrier region of induced fission. The configuration space is composed of seniority-zero configurations constructed from self-consistent mean-field wave functions. The Hamiltonian matrix elements between configurations include diabatic and pairing interactions between particles. Other aspects of the Hamiltonian are treated statistically, guided by phenomenological input of compound-nucleus transmission coefficients. In this exploratory study the configuration space is restricted to neutron excitations only. A key observable calculated in the model is the fission-to-capture branching ratio. We find that both pairing and diabatic interactions are important for achieving large branching to the fission channels. In accordance with the transition-state theory of fission, the calculated branching ratio is found to be quite insensitive to the fission decay widths of the pre-scission configurations. Furthermore, the barrier-top dynamics appear to be quite different from transition-state theory in that the transport is distributed over many excited configurations at the barrier top.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Entanglement entropy, single-particle occupation probabilities, and short-range correlations

For quantum many-body systems with short-range correlations (SRCs), the intimate relationship between their magnitude, the behavior of the single-particle occupation probabilities at momenta larger than the Fermi momentum, and the entanglement entropy is a new qualitative aspect not studied and exploited yet. A large body of recent condensed matter studies indicates that the time evolution of the entanglement entropy describes the nonequilibrium dynamics of isolated and strongly interacting many-body systems, in a manner similar to the Boltzmann entropy, which is strictly defined for dilute and weakly interacting many-body systems. Both theoretical and experimental studies in nuclei and cold atomic gases have shown that the fermion momentum distribution has a generic behavior n(k)=C/k 4 at momenta larger than the Fermi momentum, due to the presence of SRCs, with approximately 20% of the particles having momenta larger than the Fermi momentum. Further, the presence of the long momentum tails in the presence of SRCs changes the textbook relation between the single-particle kinetic energy and occupation probabilities, n mf ⁡(k) = 1/{1+ exp ⁡β[ϵ⁡(k)-μ]} for momenta very different form the Fermi momentum, particularly for dynamics processes. SRCs induced high-momentum tails of the single-particle occupation probabilities increase the entanglement entropy of fermionic systems, which in its turn affects the dynamics of many nuclear reactions, such as heavy-ion collisions and nuclear fission.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Dynamical model of J/ψ photoproduction on the nucleon

Here, a dynamical model based on a phenomenological charm quark-nucleon (c – N) potential v cN and the Pomeronexchange mechanism is constructed to investigate the J/ψ photoproduction on the nucleon from threshold to invariant mass W = 300 GeV. The J/ψ – N potential, V J/ψN (r), is constructed by folding v cN into the wave function φ J/ψ ($c\bar{c}$) of J/ψ within a constituent quark model (CQM) of Segovia et al. [Int. J. Mod. Phys. E 22, 1330026 (2013)]. A photoproduction amplitude is also generated by v cN by a $c\bar{c}$–loop integration over the γ → $c\bar{c}$ vertex function and φ J/ψ ($c\bar{c}$). No commonly used vector meson dominance assumption is used to define this photoproduction amplitude which is needed to describe the data near the threshold. The c – N potential v cN (r) is parameterized in a form such that the predicted V J/ψN (r) at large distances has the same Yukawa potential form extracted from a lattice QCD (LQCD) calculation of Kawanai and Sasaki, [Phys. Rev. D 82, 091501(R) (2010)]. The parameters of vcN are determined by fitting the total cross-section data of Jefferson Laboratory (JLab) by performing calculations that include J/ψ – N final-state interactions (FSI). The resulting differential cross sections dσ /dt are found in good agreements with the data. It is shown that the FSI effects dominate the cross section in the very near-threshold region, allowing for sensitive testing of the predicted J/ψ – N scattering amplitudes. By imposing the constraints of J/ψ – N potential extracted from the LQCD calculation of Kawanai and Sasaki, [Phys. Rev. D 82, 091501(R) (2010)], we have obtained three J/ψ – N potentials which fit the JLab data equally well. The resulting J/ψ – N scattering lengths are in the range of a = [-0.05, -0.25] fm. With the determined v cN (r) and the wave functions generated from the same CQM, the constructed model is used to predict the cross sections of photoproduction of η c (1S) and ψ(2S) mesons for future experimental tests.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Single-channel and single-energy partial-wave analysis with continuity improved through minimal phase constraints

Single-energy partial-wave analysis has often been applied as a way to fit data with minimal model dependence. However, remaining unconstrained, partial waves at neighboring energies will vary discontinuously because the overall amplitude phase cannot be determined through single-channel measurements. This problem can be mitigated through the use of a constraining penalty function based on an associated energy-dependent fit. However, the weight given to this constraint results in a biased fit to the data. In this paper, for the first time, we explore a constraining function which does not influence the fit to data. The constraint comes from the overall phase found in multichannel fits which, in the present study, are the Bonn-Gatchina and Jülich-Bonn multichannel analyses. The data are well reproduced and weighting of the penalty function does not influence the result. The method is applied to K⁢Λ photoproduction data and all observables can be maximally well reproduced. While the employed multichannel analyses display very different multipole amplitudes, we show that the major difference between two sets of multipoles can be related to the different overall phases.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Two-neutrino double-𝛽 decay in pionless effective field theory from a Euclidean finite-volume correlation function

Two-neutrino double-β decay of certain nuclear isotopes is one of the rarest Standard Model processes observed in nature. Its neutrinoless counterpart is an exotic lepton-number nonconserving process that is widely searched to determine if the neutrinos are Majorana fermions. In order to connect the rate of these processes to the Standard Model and beyond the Standard Model interactions, it is essential that the corresponding nuclear matrix elements are constrained reliably from theory. Lattice quantum chromodynamics (LQCD) and low-energy effective field theories (EFTs) are expected to play an essential role in constraining the matrix element of the two-nucleon subprocess, which could in turn provide the input into ab initio nuclear-structure calculations in larger isotopes. Focusing on the two-neutrino process $nn\rightarrow pp(ee\bar{v}_{e}\bar{v}_e)$, the amplitude is constructed in this work in pionless EFT at next-to-leading order, demonstrating the emergence of a renormalization-scale independent amplitude and the absence of any new low-energy constant at this order beyond those present in the single-weak process. Most importantly, it is shown how a LQCD four-point correlation function in Euclidean and finite-volume spacetime can be used to constrain the Minkowski infinite-volume amplitude in the EFT. The same formalism is provided for the related single-weak process, which is an input to the double-β decay formalism. The LQCD-EFT matching procedure outlined for the double-weak amplitude paves the road toward constraining the two-nucleon matrix element entering the neutrinoless double-β decay amplitude with a light Majorana neutrino.

79 ASTRONOMY AND ASTROPHYSICS↗