Search NASA⌕ Search

SEARCH · Search NASA

Results for “nuclear energy density functional”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Theoretical Prediction and Interpretation of 237 Np Mössbauer Isomer Shifts

Different theoretical approaches for the calculation of 237 Np Mössbauer isomer shifts are investigated. The traditional contact density route is compared to a previously proposed alternative approach that uses energy derivatives with respect to the nuclear radius. Both approaches yield similar results as long as suitable basis sets augmented with large exponents and relativistic methods are used. Density functional theory (DFT) calculations do not show a strong dependency of the 237 Np isomer shift on the chosen functional. Wavefunction calculations show that dynamic electron correlation can be important when covalent bonding influences the isomer shift. Effects from spin–orbit coupling are small. The isomer shifts of ionic solids and Np(III) organometallic complexes are largely governed by the oxidation state of Np. Isomer shifts of organometallic Np(IV) complexes are strongly affected by donation bonding. Furthermore, detailed analysis of the wavefunction results with different active spaces demonstrates that correlation among the outer core Np and occupied ligand frontier orbitals contributes significantly to isomer shifts of Np(IV) compounds.

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

237 Np Mössbauer Isomer Shifts: A Lesson About the Balance of Static and Dynamic Electron Correlation in Heavy Element Complexes

A large set of neptunium compounds with different oxidation states (III to VII) was assembled to study the Mössbauer isomer shift by wave function calculations and better understand covalency in f-elements complexes. The contact density approach was used to calculate the isomer shift using complete active space self-consistent field (CASSCF) multiconfiguration wave functions, as well as matrix product states [from Density Matrix Renormalization Group (DMRG) algorithms] for large active spaces. Dynamic correlation effects for the isomer shifts were treated via CASPT2 energy derivatives with respect to the nuclear radius. The CASSCF calculations appear to produce different orbital overlocalization errors for low and high Np oxidation states. For compounds with low Np oxidation numbers, the errors can be attributed to the overlocalization of the 5f orbitals. For the compounds with high Np oxidation numbers, the main errors arise from an overlocalization of ligand orbitals and concomitant to weak donation bonding. Attempts to mitigate the overlocalization errors with large active spaces using DMRG were only partially successful, showing that explicit treatment of dynamic correlation is necessary for accurate predictions of Mössbauer isomer shifts. The CASPT2 calculations perform very satisfactorily. For a subset of Np compounds, both static and dynamic correlation effects were substantial. Furthermore, a rational active space selection based on orbital entanglement diagrams proved beneficial for determining the optimal reference wave function.

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

Graph-theoretical KMC Framework for Calculating Effective Diffusivity in TPBAR Components: Effective Diffusivity of Tritium in α-Zr Grain Boundaries

We report the development of a computer simulation tool based on the Graph-theoretical kinetic Monte Carlo (GT-KMC) approach that can be used to simulate diffusion and derive effective diffusivity of species in complex structures of TPBAR components. The GT-KMC framework was implemented in AKSOME, an on-lattice self-learning kinetic Monte Carlo tool developed at PNNL and was previously used to study solute diffusion in metal alloys. This newly developed tool will be helpful in simulating atom diffusion and extracting their effective diffusivities in grain boundaries and interfaces in various TPBAR components. Once benchmarked, the tool will be used to simulate the diffusion of tritium (hydrogen) along grain boundaries in α-Zr using the activation energy barrier data previously obtained from Density Functional Theory calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Charge Radius of Neutron-Deficient Ni 54 and Symmetry Energy Constraints Using the Difference in Mirror Pair Charge Radii

he nuclear root-mean-square charge radius of Ni 54 was determined with collinear laser spectroscopy to be R ( Ni 54 ) = 3.737 ( 3 ) fm . In conjunction with the known radius of the mirror nucleus Fe 54 , the difference of the charge radii was extracted as Δ R ch = 0.049 ( 4 ) fm . Based on the correlation between Δ R ch and the slope of the symmetry energy at nuclear saturation density ( L ), we deduced 21 ≤ L ≤ 88 MeV . The present result is consistent with the L from the binary neutron star merger GW170817, favoring a soft neutron matter EOS, and barely consistent with the PREX-2 result within 1 σ error bands. Our result indicates the neutron-skin thickness of Ca 48 as 0.15–0.21 fm.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Implications of PREX-2 on the electric dipole polarizability of neutron-rich nuclei

Background: The recent announcement by the PREX collaboration of an unanticipated thick neutron skin in 208 Pb ($R^{208}_{skin}$) has challenged the understanding of neutron-rich matter in the vicinity of nuclear saturation density. Whereas earlier constraints indicate that the symmetry energy is relatively soft, the PREX-2 result seems to suggest the opposite. Purpose: To confront constraints on the symmetry energy obtained from measurements of the electric dipole polarizability against those informed by the PREX-2 measurement of $R^{208}_{skin}$ and by the correlations that it entails. Methods: Here, covariant energy density functionals informed by the properties of finite nuclei are used to compute the electric dipole response of 48 Ca, 68 Ni, 132 Sn, and 208 Pb. The set of functionals used in this work are consistent with experimental data, yet are flexible enough in that they span a wide range of values of $R^{208}_{skin}$ Results: It is found that theoretical predictions of the electric dipole polarizability that are consistent with the PREX-2 measurement systematically overestimate the corresponding values extracted from the direct measurements of the distribution of electric dipole strength. Conclusions: The neutron skin thickness of 208 Pb extracted from parity violating electron scattering and the electric dipole polarizability measured in photoabsorption experiments are two of the cleanest experimental tools used to constrain the symmetry energy around nuclear saturation density. However, the recent value of $R^{208}_{skin}$ that suggests a fairly stiff symmetry energy stands in stark contrast to the conclusions derived from the electric dipole polarizability. At present, I offer no solution to this dilemma.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Constraining the nuclear spin distribution using improved 197 Au neutron resonance parameters

New neutron transmission data at resonance energies using a 197 Au sample were measured using an early version of the Device for Indirect Capture Experiments on Radionuclides (DICER), which is under development at the Los Alamos Neutron Science Center (LANSCE). These data were combined with previous neutron transmission and capture data in a simultaneous R-matrix analysis to extract improved neutron resonance parameters for this nuclide. As a result, total radiation widths, Γ γ , were obtained for 33 J=1 and 44 J=2 197 Au+n resonances. Γ γ distributions for these two spins states were compared to distributions calculated according to the nuclear statistical model using published nuclear level density (NLD) and photon strength functions (PSF) measured using the Oslo technique. The calculated distributions were found to be narrower and the average values for the two spins states closer together than the data. The calculation can be brought into agreement with the data by substantial modifications to the spin distribution in 198 Au as a function of excitation energy. As far as we know, the spin distribution currently is otherwise poorly constrained. The modified spin distribution changes the shapes of the NLD and PSF extracted using the Oslo technique and so could have broad implications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Probing the Gluonic Structure of the Deuteron with J / ψ Photoproduction in d + Au Ultraperipheral Collisions

Understanding gluon density distributions and how they are modified in nuclei are among the most important goals in nuclear physics. In recent years, diffractive vector meson production measured in ultraperipheral collisions (UPCs) at heavy-ion colliders has provided a new tool for probing the gluon density. In this Letter, we report the first measurement of J/psi photoproduction off the deuteron in UPCs at the center-of-mass energy root √S-NN = 200 GeV in d + Au collisions. The differential cross section as a function of momentum transfer -t is measured. In addition, data with a neutron tagged in the deuteron-going zero-degree calorimeter is investigated for the first time, which is found to be consistent with the expectation of incoherent diffractive scattering at low momentum transfer. Theoretical predictions based on the color glass condensate saturation model and the leading twist approximation nuclear shadowing model are compared with the data quantitatively. A better agreement with the saturation model has been observed. With the current measurement, the results are found to be directly sensitive to the gluon density distribution of the deuteron and the deuteron breakup process, which provides insights into the nuclear gluonic structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Mitigating Green's function Monte Carlo signal-to-noise problems using contour deformations

The Green's function Monte Carlo (GFMC) method provides accurate solutions to the nuclear many-body problem and predicts properties of light nuclei starting from realistic two- and three-body interactions. Controlling the GFMC fermion-sign problem is crucial, as the signal-to-noise ratio decreases exponentially with Euclidean time, requiring significant computing resources. Inspired by similar scenarios in lattice quantum field theory and spin systems, in this work, we employ integration contour deformations to improve the GFMC signal-to-noise ratio. Machine learning techniques are used to select optimal contours with minimal variance from parameterized families of deformations. As a proof of principle, we consider the deuteron binding energies and Euclidean density response functions. We only observe mild signal-to-noise improvement for the binding energy case. On the other hand, we achieve an order of magnitude reduction of the variance for Euclidean density response functions, paving the way for computing electron- and neutrino-nucleus cross-sections of larger nuclei.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Observation of $t\bar{t}$ Production in Pb+Pb Collisions at $\sqrt{s_{NN}}$ = 5.02 TeV with the ATLAS Detector

Top-quark pair production is observed in lead–lead (Pb+Pb) collisions at $\sqrt{s_{NN}}$ = 5.02 TeV at the Large Hadron Collider with the ATLAS detector. The data sample was recorded in 2015 and 2018, amounting to an integrated luminosity of 1.9 nb −1 . Events with exactly one electron and one muon and at least two jets are selected. Top-quark pair production is measured with an observed (expected) significance of 5.0 (4.1) standard deviations. The measured top-quark pair production cross section is 𝜎 𝑡⁢$\tilde{i}$ = 3.6⁢$^{+1.0}_{−0.9}$⁢(stat)⁢ $^{+0.8}_{−0.5}$⁢(syst) μ⁢b, with a total relative uncertainty of 31%, and is consistent with theoretical predictions using a range of different nuclear parton distribution functions. The observation of this process consolidates the evidence of the existence of all quark flavors in the preequilibrium stage of the quark-gluon plasma at very high energy densities, similar to the conditions present in the early Universe.

Quark-gluon plasma↗
DOE OSTI2026

Enhancements in Laser-Direct-Drive Nuclear Performance with Target Radius

Inertial confinement fusion is likely to require significant improvements in technology and design for large targets to implode at high driver energies. Here, to assess the potential for benefits, we report on nuclear performance as a function of target radius 𝑅 t in direct-drive cryogenic implosions as performed on the OMEGA laser. The neutron yield and areal density are found to increase as 𝑅 5.0±0.2 t and 𝑅 1.8±0.2 t , respectively, and exhibit a much stronger dependence on target radius than previous studies have assumed. As this paper demonstrates, these types of sensitivities should be expected when implosions are unstable, and have been degraded by a range of multidimensional effects in 3D. If the imperfections in a target and laser system are fixed in magnitude, it follows that larger implosions have higher relative quality, and approach criteria to ignite much more quickly than commonly appreciated.

inertial confinement fusion↗
DOE OSTI2022

Hydrogen Localization and Cluster Formation in α-Zr from First-Principles Investigations

We report alpha-zirconium (α-Zr) is a hydrogen atom host material used in nuclear energy and semiconducting applications due to Zr's excellent properties such as corrosion resistance, low thermal neutrons cross-section, and thermodynamically favorable interstitial sequestering of H atoms. Herein we use density functional theory simulations to systematically analyze the initial stages of H insertion in the Zr matrix, and determine optimal hydrogen localization, distribution, and clustering tendencies from a thermodynamics and kinetics perspective. Our results suggest that H occupation of tetrahedral sites is strongly favored, consistent with previous experimental and theoretical studies. Further, a cooperative effect in H localization is shown where H atoms favor clustering. The H pair is found to be the most stable by occupying tetrahedral sites with pairs at second nearest-neighbor positions (TTn2). Among the four H atom clusters, the most favorable configuration has two TTn2 pairs forming "box " or "zig-zag" shaped-cluster with the pairs separated by 3.03 or 3.25 Å, respectively. Nudged elastic band simulations suggest that the formation of TTn2 is facile and is limited by the small diffusion barrier of the monomer of less than 0.1 eV. Further, the formation routes for the 4H atom clusters are also presented. The current study can serve as a foundation to systematically model H atom clustering in metals and alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Stabilization of Charge-Transfer Excited States in Biological Systems: A Computational Focus on the Special Pair in Photosystem II Reaction Centers

Charge-transfer (CT) excited states play an important role in many biological processes. However, many computational approaches often inadequately address the equilibration effects of nuclear and environmental degrees of freedom on these states. One prominent example of systems in which CT states are of utmost importance is reaction centers (RC) in photosystems. Here we use a multiscale approach combined with time-dependent density functional theory to explore the lowest CT excited state of the special pair P D1 –P D2 in the Photosystem II-RC of a cyanobacterium. We find that the nonequilibrium CT excited state resides near the Soret band, making an exciton the lowest-energy excited state. However, accounting for nuclear and state-specific dielectric equilibration along the CT potential energy surface (PES), the CT state P D1 – –P D2 + stabilizes energetically below the excitonic state. Importantly, this underscores the crucial role of state-specific solvation in mapping the PES of CT states, as demonstrated in a simplified dimer model.

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

Bayesian inference of fine features of the nuclear equation of state from future neutron star radius measurements to 0.1 km accuracy

To more precisely constrain the equation of state (EOS) of supradense neutron-rich nuclear matter, future high-precision x-ray and gravitational wave observatories are proposed to measure the radii of neutron stars (NSs) with an accuracy better than about 0.1 km. However, it remains unclear what particular aspects (other than the stiffness generally spoken of in the literature) of the EOS and to what precision they will be better constrained. In this work, within a Bayesian framework using a metamodel EOS for NSs, we infer the posterior probability distribution functions (PDFs) of incompressibility K 0 and skewness J 0 of symmetric nuclear matter (SNM) as well as the slope L, curvature K sym , and skewness J sym characterizing the density dependence of nuclear symmetry energy E sym ⁡(ρ), respectively, from mean values of NS radii consistent with existing observations and an expected accuracy Δ⁢R ranging from about 1.0 to 0.1 km. Here, we found that (1) the Δ⁢R has little effect on inferring the stiffness of SNM at suprasaturation densities, (2) smaller Δ⁢R reveals more accurately not only the PDFs but also pairwise correlations among parameters characterizing high-density E sym ⁡(ρ), (3) a double-peak feature of the PDF(K sym ) corresponding to the strong K sym – J sym and K sym – L anticorrelations is revealed when Δ⁢R is less than about 0.2 km, and the locations of the two peaks are sensitive to the maximum value of J sym reflecting the stiffness of E sym ⁡(ρ) above about 3 times the saturation density ρ 0 of SNM, and (4) the high-precision radius measurement for canonical NSs is more useful than that for massive ones for constraining the EOS of nucleonic matter around (2–3)⁢ρ 0 .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Nuclear Computational Low Energy Initiative (NUCLEI)

The NUCLEI project, as defined by the scope of work, developed, implemented and run codes for large-scale computations of many topics in low-energy nuclear physics. Physics studied include the properties of nuclei and nuclear decays, nuclear structure and reactions, and the properties of nuclear matter. The computational techniques used include Quantum Monte Carlo, Configuration Interaction, Coupled Cluster, and Density Functional methods. The research program emphasized areas of high interest to current and possible future DOE nuclear physics facilities, including ATLAS and FRIB (nuclear structure and reactions, and nuclear astrophysics), TJNAF (neutron distributions in nuclei, few body systems, and electroweak processes), NIF (thermonuclear reactions), MAJORANA and FNPB (neutrinoless double-beta decay and physics beyond the Standard Model), and LANSCE (fission studies).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Weighted active space protocol for multireference machine-learned potentials

Multireference methods such as multiconfiguration pair-density functional theory accurately capture electronic correlation in systems with strong multiconfigurational character, but their cost precludes direct use in molecular dynamics. Combining these methods with machine-learned interatomic potentials (MLPs) can extend their reach. However, the sensitivity of multireference calculations to the choice of the active space complicates the consistent evaluation of energies and gradients across structurally diverse nuclear configurations. To overcome this limitation, we introduce the weighted active space protocol (WASP), a systematic approach to assign a consistent active space for a given system across uncorrelated configurations. By integrating WASP with MLPs and enhanced sampling techniques, we propose a data-efficient active learning cycle that enables the training of an MLP on multireference data. We demonstrated the approach on the TiC + -catalyzed C–H activation of methane, a reaction that poses challenges for Kohn–Sham density functional theory due to its significant multireference character. This framework enables accurate and efficient modeling of catalytic dynamics, establishing a paradigm for simulating complex reactive processes beyond the limits of conventional electronic-structure methods.

enhanced sampling↗
DOE OSTI2020

Efficient computational search for lanthanide-binding additive dopants for advanced U-Zr based fuels

A rational design strategy combining theory and experiment is highly desirable for the development of U-Zr based fuels that are resistant to lanthanide-induced fuel-cladding chemical interaction (FCCI). Here, we performed rapid computational screening of FCCI-mitigating (via strong lanthanide-binding) elements across the periodic table by utilizing three simple criteria and density functional theory (DFT) calculated total energies of relevant elements and compounds retrieved from high-throughput DFT databases. In addition to successfully identifying previously demonstrated additives, our search leads to the discovery of Bi as a highly effective element among all 68 elements screened, and its effectiveness was experimentally confirmed in this study.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Symmetry Breaking in the Lowest-Lying Excited-State of CCl4: Valence Shell Spectroscopy in the 5.0–10.8 eV Photon Energy Range

We report absolute high-resolution vacuum ultraviolet (VUV) photoabsorption cross-sections of carbon tetrachloride (CCl4) in the photon energy range 5.0–10.8 eV (248–115 nm). The molecular spectrum and electronic structure have been comprehensively investigated together with quantum chemical calculations, providing geometries, bond lengths, vertical excitation energies and oscillator strengths. The major electronic excitations have been assigned to valence and Rydberg transitions which are also accompanied by vibrational excitation assigned to degenerate stretching, v3′t2 and degenerate deformation v4′t2 modes. The rather complex nuclear dynamics along the degenerate deformation mode, v4′t2, have been thoroughly investigated by Time-Dependent Density Functional Theory (TD-DFT) method. The relevant Jahn–Teller distortion operative within the lowest-lying electronic excited-state is shown here for the first time in order to yield a weak absorption feature at 6.156 eV. Further calculations on the potential energy curves for the singlet excited-states along the C–Cl stretching coordinate show the relevance of efficient C–Cl bond excision.

Biochemistry & Molecular Biology↗
DOE OSTI2024

Water–Hydrocarbon Interactions in Anionic Pyrene Monohydrate

Interactions between water and polycyclic aromatic hydrocarbons are essential in many aspects of chemistry, from interstellar and atmospheric processes to interfacial hydrophobicity and wetting phenomena. Despite their growing importance, the intermolecular potentials of the water-hydrocarbon interactions are underdeveloped compared to water-water potentials, and there are similarly few experimental probes that are sensitive to the details of the water-hydrocarbon potential. We present a combined experimental and computational study of anionic pyrene monohydrate, one of the simplest water/hydrocarbon clusters. The action spectrum in the OH region of the mass-selected cluster ion provides a rigorous benchmark for intermolecular potentials and computational methodologies. We identify missing intermolecular interactions and shortcomings in conventional dynamics calculations by comparing experimental data to density functional theory and classical molecular dynamics calculations. Kinetic trapping is prevalent, even for one water molecule and one pyrene molecule, leading to slow equilibration in conventional molecular dynamics calculations, even on nanosecond timescales and at low temperatures (50 K). At constant energy, temperature fluctuations for the pair of molecules are substantial. Immersing the system in a bath of soft spheres and employing parallel tempering alleviates kinetic trapping and dampens temperature fluctuations, bringing the system closer to the thermodynamic limit. With such augmented sampling, a simple, flexible water model reproduces the linewidth and the asymmetric broadening of the symmetric OH stretching mode, which we assign to spectral diffusion. In the OH stretching region, dynamics calculations predict a more intense antisymmetric peak than experiments observe but do not predict the bimodal split symmetric peak that the experiments show. Furthermore, our work suggests that electronic polarization, missing in the empirical force field, is responsible for the first discrepancy and that quantum nuclear effects, captured neither in density functional theory nor in classical dynamics, may be responsible for the second.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗