Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nuclei”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Investigating the crust of neutron stars with neural-network quantum states

An accurate description of low-density nuclear matter is crucial for explaining the physics of neutron star crusts. In the density range between approximately 0.01 fm −3 and 0.1 fm −3 , matter transitions from neutron-rich nuclei to various higher-density pasta shapes, before ultimately reaching a uniform liquid. In this work, we introduce a variational Monte Carlo method based on a neural Pfaffian-Jastrow quantum state, which allows us to model the transition from the liquid phase to neutron-rich nuclei microscopically. At low densities, nuclear clusters dynamically emerge from the microscopic interactions among protons and neutrons, which we model based on pionless effective field theory. Our variational Monte Carlo approach represents a significant improvement over the state-of-the-art auxiliary-field diffusion Monte Carlo method, which is severely hindered by the fermion-sign problem in this low-density regime and cannot capture the onset of clusters. In addition to computing the energy per particle of symmetric nuclear matter and pure neutron matter, we analyze an intermediate isospin-asymmetry configuration to elucidate the formation of nuclear clusters. We also provide evidence that the presence of such nuclear clusters influences the amount of protons in the crust compared to protons in beta-equilibrated, neutrino-transparent matter.

Nuclear astrophysics↗
DOE OSTI2025

Discovering nuclear models from symbolic machine learning

Numerous phenomenological nuclear models have been proposed to describe specific observables within different regions of the nuclear chart. However, developing a unified model that describes the complex behavior of all nuclei remains an open challenge. Here, we explore whether symbolic Machine Learning (ML) can rediscover traditional nuclear physics models or identify alternatives with improved simplicity, fidelity, and predictive power. To address this challenge, we developed a Multi-objective Iterated Symbolic Regression approach that handles symbolic regressions over multiple target observables, accounts for experimental uncertainties and is robust against high-dimensional problems. As a proof of principle, we applied this method to describe the nuclear binding energies and charge radii of light and medium mass nuclei. Our approach identified simple analytical relationships based on the number of protons and neutrons, providing interpretable models with precision comparable to state-of-the-art nuclear models. Additionally, we integrated this ML-discovered model with an existing complementary model to estimate the limits of nuclear stability. These results highlight the potential of symbolic ML to develop accurate nuclear models and guide our description of complex many-body problems.

Nuclear structure↗
DOE OSTI2023

The quest for superheavy elements and the limit of the periodic table

The borders of the periodic table of the elements and of the chart of nuclides are not set in stone. The desire to explore the properties of atoms and their nuclei in a regime of very large numbers of electrons, protons and neutrons has motivated new experimental facilities to create new elements and nuclides at the limits of atomic number and mass. But the small production rates and short lifetimes of superheavy nuclei and their atoms mean that ‘atom-at-a-time’ studies are the only experimental way to probe them. The physical and chemical data obtained so far, augmented by theoretical calculations, indicate significant deviations from extrapolations from lighter elements and isotopes. This situation raises the following question: how much further can one push the limits of the periodic table? Here, in this Review, we describe the major challenges in the field of the superheavy elements and speculate about future directions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

O Corona, where art thou? eROSITA’s view of UV-optical-IR variability-selected massive black holes in low-mass galaxies

Finding massive black holes (MBHs,M BH ≈ 10 4 –10 7 M ⊙ ) in the nuclei of low-mass galaxies $\left( {{M_*}\mathop {\mathop < \limits_ }\limits_ {{10}^{10}}{M_ \odot }} \right)$ is crucial to constrain seeding and growth of black holes over cosmic time, but it is particularly challenging due to their low accretion luminosities. Variability selection via long-term photometric ultraviolet, optical, or infrared (UVOIR) light curves has proved effective and identifies lower-Eddington ratios compared to broad and narrow optical spectral lines searches. In the inefficient accretion regime, X-ray and radio searches are effective, but they have been limited to small samples. Therefore, differences between selection techniques have remained uncertain. Here, we present the first large systematic investigation of the X-ray properties of a sample of known MBH candidates in dwarf galaxies. We extracted X-ray photometry and spectra of a sample of ~200 UVOIR variability-selected MBHs and significantly detected 17 of them in the deepest available SRG/eROSITA image, of which four are newly discovered X-ray sources and two are new secure MBHs. This implies that tens to hundreds of LSST MBHs will have SRG/eROSITA counterparts, depending on the seeding model adopted. Surprisingly, the stacked X-ray images of the many non-detected MBHs are incompatible with standard disk-corona relations, typical of active galactic nuclei, inferred from both the optical and radio fluxes. They are instead compatible with the X-ray emission predicted for normal galaxies. After careful consideration of potential biases, we identified that this X-ray weakness needs a physical origin. A possibility is that a canonical X-ray corona might be lacking in the majority of this population of UVOIR-variability selected low-mass galaxies or that unusual accretion modes and spectral energy distributions are in place for MBHs in dwarf galaxies. This result reveals the potential for severe biases in occupation fractions derived from data from only one waveband combined with SEDs and scaling relations of more massive black holes and galaxies.

Astronomy & Astrophysics↗
DOE OSTI2020

Correlations in the Initial Conditions of Heavy-Ion Collisions

Ultracentral collisions of heavy nuclei, in which the impact parameter is nearly zero, are especially sensitive to the details of the initial state model and the microscopic mechanism for collective flow. In a hydrodynamic “flow” picture, the final state momentum correlations are a direct response to the fluctuating initial geometry, although models of the initial geometry differ widely. Alternatively, dynamical mechanisms based in the color glass condensate (CGC) formalism can naturally lead to many-body correlations with very different systematics. Here we present a calculation of event-by-event elliptic flow in both the hydrodynamic and CGC paradigms and show that they can be qualitatively distinguished in ultracentral collisions of deformed nuclei. Specifically, the multiplicity dependence in such collisions is qualitatively opposite, with the CGC correlations increasing with multiplicity while the hydrodynamic correlations decrease. The consistency of the latter with experimental data on UU collisions appears to rule out a CGC-mediated explanation. We find that these qualitative features also persist in small deformed systems and can therefore be a valuable test of the microscopic physics in that regime.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Following Fission Products in Explosive Astrophysical Environments

The astrophysical process by which the heaviest elements are formed in the universe is known as the rapid neutron capture process, or r process, of nucleosynthesis. The r process is characterized by the neutron capture and β - decay of short-lived, neutron-rich atomic nuclei; in suitably extreme environments, nuclear fission can also play a major role in determining the ensuing nucleosynthesis. In this work, we present the application of our recently developed nucleosynthesis tracing framework to precisely quantify the impact that neutron-induced and β - -delayed fission processes have in r-process environments that produce fissioning nuclei.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Building bridges between nuclear structure and reactions

The past couple of decades have seen tremendous advances in nuclear structure and reaction theory. Innovative theory frameworks for describing the nuclear many-body system, increasingly powerful computers, and opportunities for confronting theory predictions with data on unstable nuclei, have been driving the field. An important goal is to move from phenomenological ingredients in reaction calculations to predictive theories based on microscopic frameworks. We discuss ongoing efforts aimed at integrating microscopic descriptions of nuclear structure into reaction predictions for medium-mass and heavy nuclei. This contribution highlights areas where Eric Bauge, a champion for building bridges, has made important contributions by encouraging and enabling collaborations between communities with complementary expertise.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Microscopic-macroscopic level densities for low excitation energies

Level density ρ(E,Q) is derived within the micro-macroscopic approximation (MMA) for a system of strongly interacting Fermi particles with the energy E and additional integrals of motion Q , in line with several topics of the universal and fruitful activity of A. S. Davydov. Within the extended Thomas Fermi and semiclassical periodic orbit theory beyond the Fermi-gas saddle-point method, we obtain ρ ∝ I ν (S)/S ν , where I ν (S) is the modified Bessel function of the entropy S . For small shell-structure contribution, one finds ν = κ/2 + 1, where κ is the number of additional integrals of motion. This integer number is a dimension of Q, Q = { N, Z , …} for the case of two-component atomic nuclei, where N and Z are the numbers of neutrons and protons, respectively. For much larger shell structure contributions, one obtains ν = κ /2 + 2. The MMA level density ρ reaches the well-known Fermi gas asymptote for large excitation energies and the finite micro-canonical combinatoric limit for low excitation energies. Further, the additional integrals of motion can also be the projection of the angular momentum of a nuclear system for nuclear rotations of deformed nuclei, number of excitons for collective dynamics, and so on. Fitting the MMA total level density ρ( E , Q) for a set of the integrals of motion Q = { N, Z }, to experimental data on a long nuclear isotope chain for low excitation energies, one obtains the results for the inverse level-density parameter K , which differs significantly from those of neutron resonances due to shell, isotopic asymmetry, and pairing effects.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Ion irradiation and modification: The role of coupled electronic and nuclear energy dissipation and subsequent nonequilibrium processes in materials

Understanding material responses to energy deposition from energetic charged particles is important for defect engineering, ion-beam processing, ion-beam analysis and modification, geologic aging, space exploration, and nuclear applications. As an incident ion penetrates a solid, its energy is transferred to electrons and to atomic nuclei of the solid. Much of this electronic energy deposition is subsequently transferred to the atomic structure via electron–phonon (e–ph) coupling, leading to local inelastic thermal spikes in which energy dissipation is influenced by the local environment. In addition, intense ionization can lead to high densities of localized electronic excitations in wide-bandgap materials and ceramics that can affect defect dynamics and atomic mobility. Specifically, energy exchange between electrons and atomic nuclei, along with localized electronic excitations, can lead to substantial competitive (ionization-induced annealing), additive (both electronic and nuclear energy depositions contributing to damage production), and synergistic (more damage than the sums of separate damage processes) effects. Although nonmonotonic effects of the e–ph coupling strength and athermal processes are demonstrated for pre-existing defects and residual damage during ion–solid interactions, there is limited understanding of when such electronic effects must be considered in atomic-scale models of damage production and evolution in a broad variety of materials. Complex ceramics and chemically disordered solid solution alloys with different constituent elements allow a systematic evaluation of defect dynamics and irradiation performance with increasing complexity. Current knowledge regarding tuning of bonding characteristics and chemical disorder to control atomic-level dynamics is reviewed. Although a lack of fundamental understanding obstructs the advancement of reliable predictions for ion beam material modification, it highlights challenges and opens research opportunities. Insights into the complex electronic and atomic correlations with extreme energy deposition will strengthen our ability to design materials and predict ion-irradiation-induced damage in a radiation environment, and they may pave the way to better control fundamental processes and design new material functionalities for advanced technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Stochastic density functional theory: Real- and energy-space fragmentation for noise reduction

Stochastic density functional theory (sDFT) is becoming a valuable tool for studying ground-state properties of extended materials. The computational complexity of describing the Kohn-Sham orbitals is replaced by introducing a set of random (stochastic) orbitals leading to linear and often sub-linear scaling of certain ground-state observables at the account of introducing a statistical error. Schemes to reduce the noise are essential, for example, for determining the structure using the forces obtained from sDFT. Recently, we have introduced two embedding schemes to mitigate the statistical fluctuations in the electron density and resultant forces on the nuclei. Both techniques were based on fragmenting the system either in real space or slicing the occupied space into energy windows, allowing for a significant reduction in the statistical fluctuations. For chemical accuracy, further reduction of the noise is required, which could be achieved by increasing the number of stochastic orbitals. However, the convergence is relatively slow as the statistical error scales as 1 / N &#x3C7; according to the central limit theorem, where N χ is the number of random orbitals. In this paper, we combined the embedding schemes mentioned above and introduced a new approach that builds on overlapped fragments and energy windows. The new approach significantly lowers the noise for ground-state properties, such as the electron density, total energy, and forces on the nuclei, as demonstrated for a G-center in bulk silicon.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Beam particle identification and tagging of incompletely stripped heavy beams with HEIST

A challenge preventing successful inverse kinematics measurements with heavy nuclei that are not fully stripped is identifying and tagging the beam particles. For this purpose, the HEavy ISotope Tagger (HEIST) has been developed. HEIST utilizes two micro-channel plate timing detectors to measure the time-of-flight, a multi-sampling ion chamber to measure energy loss, and a high-purity germanium detector to identify isomer decays and calibrate the isotope identification system. HEIST has successfully identified 198 Pb and other nearby nuclei at energies of about 75 MeV/A. In the experiment discussed, a typical cut containing 89% of all 198 Pb 80+ in the beam had a purity of 86%. We examine the issues of charge state contamination. Here, the observed charge state populations of these ions are presented and, using an adjusted beam energy, are well described by the charge state model GLOBAL.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

An in-silico NMR laboratory for nuclear magnetic shieldings computed via finite fields: Exploring nucleus-specific renormalizations of MP2 and MP3

Here we developed and implemented a method-independent, fully numerical, finite difference approach to calculating nuclear magnetic resonance shieldings, using gauge-including atomic orbitals. The resulting capability can be used to explore non-standard methods, given only the energy as a function of finite-applied magnetic fields and nuclear spins. For example, standard second-order Møller-Plesset theory (MP2) has well-known efficacy for 1H and 13C shieldings and known limitations for other nuclei such as 15N and 17O. It is, therefore, interesting to seek methods that offer good accuracy for 15N and 17O shieldings without greatly increased compute costs, as well as exploring whether such methods can further improve 1H and 13C shieldings. Using a small molecule test set of 28 species, we assessed two alternatives: κ regularized MP2 (κ-MP2), which provides energy-dependent damping of large amplitudes, and MP2.X, which includes a variable fraction, X, of third-order correlation (MP3). The aug-cc-pVTZ basis was used, and coupled cluster with singles and doubles and perturbative triples [CCSD(T)] results were taken as reference values. Our κ-MP2 results reveal significant improvements over MP2 for 13C and 15N, with the optimal κ value being element-specific. κ-MP2 with κ = 2 offers a 30% rms error reduction over MP2. For 15N, κ-MP2 with κ = 1.1 provides a 90% error reduction vs MP2 and a 60% error reduction vs CCSD. On the other hand, MP2.X with a scaling factor of 0.6 outperformed CCSD for all heavy nuclei. These results can be understood as providing renormalization of doubles amplitudes to partially account for neglected triple and higher substitutions and offer promising opportunities for future applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Relativistic Douglas–Kroll–Hess calculations of hyperfine interactions within first-principles multireference methods

A relativistic magnetic hyperfine interaction Hamiltonian based on the Douglas–Kroll–Hess (DKH) theory up to the second order is implemented within the ab initio multireference methods, including spin–orbit coupling in the Molcas/OpenMolcas package. This implementation is applied to calculate relativistic hyperfine coupling (HFC) parameters for atomic systems and diatomic radicals with valence s or d orbitals by systematically varying active space size in the restricted active space self-consistent field formalism with restricted active space state interaction for spin–orbit coupling. The DKH relativistic treatment of the hyperfine interaction reduces the Fermi contact contribution to the HFC due to the presence of kinetic factors that regularize the singularity of the Dirac delta function in the nonrelativistic Fermi contact operator. This effect is more prominent for heavier nuclei. As the active space size increases, the relativistic correction of the Fermi contact contribution converges well to the experimental data for light and moderately heavy nuclei. The relativistic correction, however, does not significantly affect the spin-dipole contribution to the hyperfine interaction. In addition to the atomic and molecular systems, the implementation is applied to calculate the relativistic HFC parameters for large trivalent and divalent Tb-based single-molecule magnets (SMMs), such as Tb(III)Pc2 and Tb(II)(CpiPr5)2 without ligand truncation using well-converged basis sets. In particular, for the divalent SMM, which has an unpaired valence 6s/5d hybrid orbital, the relativistic treatment of HFC is crucial for a proper description of the Fermi contact contribution. Even with the relativistic hyperfine Hamiltonian, the divalent SMM is shown to exhibit strong tunability of HFC via an external electric field (i.e., strong hyperfine Stark effect).

Chemistry↗
DOE OSTI2025

Multistage nucleation pathway in LiF molten salt mirrors the crystal–melt interface structure

Despite over a century of studies, fundamental questions remain about the processes governing crystal nucleation from melts or solutions. Research over the past three decades has presented mounting evidence for kinetic pathways of crystal nucleation that are more complex than envisioned by the simplest forms of classical theory. Such observations have been presented for colloidal and elemental systems with covalent and metallic bonding. Despite the technological and geochemical importance of molten salts, similar studies for these ionically bonded systems are currently lacking. Here we develop a machine learning interatomic potential for a model ionic system: LiF. The potential features quantum-level accuracy for both liquid and multiple solid polymorphs over wide temperature and pressure ranges and accurately reproduces experimentally measured properties. Thanks to the efficiency of the potential, which enables microsecond-scale molecular dynamics simulations, induction times for nucleation of LiF solids from their melts are computed over a range of undercoolings. With the aid of a set of robust local order parameters established here, the simulations reveal that homogeneous crystal nucleation in undercooled melts preferentially initiates from liquid regions showing slow dynamics and high bond orientational order simultaneously, and the second-shell order of both precritical nuclei and the surface of postcritical nuclei is dominated by hexagonal close packing and body-centered cubic local structure, even though the nucleus core is dominated by face-centered cubic structure corresponding to the stable rocksalt crystal structure. Finally, we establish a connection between the crystallization pathway and the equilibrium crystal-melt interface structure.

Applied Physical Sciences↗
DOE OSTI2023

Photonuclear Production of 67 Cu From Gallium

Copper-67 is a radioisotope of interest for medical imaging and therapy as well as for under- standing stellar and interstellar evolution pertaining to the formation of proton-rich nuclei. Since 67 Cu decays 100% to 67 Zinc, understanding this reaction can shed light on the abundance of this and other p-nuclei elements in the universe. Here, the photonuclear production of 67 Cu from 71 Ga and natural gallium is examined as an alternative to its photoproduction from zinc. Two research and development production runs were performed at Thomas Jefferson National Accelerator Facility using an electron linac. During the first run, an 805-W, 30.9-MeV beam was used to irradiate a 1-mm tungsten radiator to create a bremsstrahlung flux. The resulting gamma photons irradiated 50.9 g of natural gallium encased in a graphite crucible for 24.2 h; 7.02 Bq/W∙s∙kg of 67 Cu activity was produced. During the second run, a 4380-W, 31.5-MeV beam was used for 12.0 h on the same target containing 60 g of natural gallium; 6.41 Bq/W∙s∙kg of 67 Cu activity was produced. Because of the difficulties in spectroscopically differentiating 67 Cu from 67 Ga, prior to each run, an isotopically pure 71 Ga disk was irradiated using a 100-W beam for 1 h, at the same respective energies. Finally, these baseline irradiations allowed for separation of 67 Cu from 67 Ga in the spectroscopic measurements of the natural gallium targets.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2020

Benchmark calculations of electron impact electronic excitation of the hydrogen molecule

We present benchmark integrated and differential cross-sections for electron collisions with H 2 using two different theoretical approaches, namely, the R-matrix and molecular convergent close-coupling. This is similar to comparative studies conducted on electron–atom collisions for H, He and Mg. Electron impact excitation to the ${\mathrm{b}}^{3}{{\Sigma}}_{\mathrm{u}}^{+}$, ${\mathrm{a}}^{3}{{\Sigma}}_{\mathrm{g}}^{+}$, ${\mathrm{B}}^{1}{{\Sigma}}_{\mathrm{u}}^{+}$, c3Πu, ${\mathrm{E}\mathrm{F}}^{1}{{\Sigma}}_{\mathrm{g}}^{+}$, C1Πu, ${\mathrm{e}}^{3}{{\Sigma}}_{\mathrm{u}}^{+}$, ${\mathrm{h}}^{3}{{\Sigma}}_{\mathrm{g}}^{+}$, ${\mathrm{B}}^{\prime 1}{{\Sigma}}_{\mathrm{u}}^{+}$ and d3Πu excited electronic states are considered. Calculations are presented in both the fixed nuclei and adiabatic nuclei approximations, where the latter is shown only for the ${\mathrm{b}}^{3}{{\Sigma}}_{\mathrm{u}}^{+}$ state. Good agreement is found for all transitions presented. Where available, we compare with existing experimental and recommended data.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Exact sum rules with approximate ground states

Electromagnetic and weak transitions tell us a great deal about the structure of atomic nuclei. Yet modeling transitions can be difficult: it is often easier to compute the ground state, if only as an approximation, than excited states. One alternative is through transition sum rules, in particular the non-energy-weighted and energy-weighted sum rules, which can be computed as expectation values of operators. We investigate by computing sum rules for a variety of nuclei, comparing the numerically exact full configuration-interaction shell model, as a reference, to Hartree-Fock, projected Hartree-Fock, and the nucleon pair approximation. Furthermore, these approximations yield reasonable agreement, which we explain by prior work on the systematics of transition moments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Rotational bands beyond the Elliott model

Rotational bands are commonplace in the spectra of atomic nuclei. Inspired by early descriptions of these bands by quadrupole deformations of a liquid drop, Elliott constructed discrete nucleon representations of SU(3) from fermionic creation and annihilation operators. Ever since, Elliott's model has been foundational to descriptions of rotation in nuclei. Later work, however, suggested the symplectic extension Sp(3, R) provides a more unified picture. We decompose no-core shell-model nuclear wave functions into symmetry-defined subspaces for several beryllium isotopes, as well as 20 Ne, using the quadratic Casimirs of both Elliott's SU(3) and Sp(3, R). The band structure, delineated by strong B(E2) values, has a more consistent description in Sp(3, R) rather than SU(3). In particular, we confirm previous work finding in some nuclides strongly connected upper and lower bands with the same underlying symplectic structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗