Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Calculation”

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 OSTI2026

Real time heat load calculation software based on EPICS for Fermilab PIP-II CM tests

Fermilab has a project to improve the proton beam energy which is called PIP-II (the 2nd Proton Improvement Plan). There is a superconducting linear accelerator, LINAC, to improve the proton beam power and the LINAC consists of 5 types of cryomodules (CM), 1 HWR CM, 2 SSR1 CM, 4 SSR2 CM, LB650 CM, and HB650 CM. The prototypes of these cryomodules are being tested at Fermilab’s CryoModule Test Facility (CMTF). Heat load measurements are an important part of the prototype CM testing. The CMTF cryogenic control system was developed based on the ACNET (Accelerator Control NETwork) for CM testing for other projects, but the PIP-II cryogenic control system will be implemented using the Experimental Physics and Industrial Control System (EPICS). As part of the prototype CM testing campaign, an EPICS based control system has been implemented at CMTF. This EPICS cryogenic control system includes real-time heat load calculation software utilizing the Fortran implementation of Hepak. This paper details the real time heat load calculation software developed for the prototype CM testing including the first results from the HB650 CM.

Yoon, S. [Fermilab]↗
DOE OSTI2026

Untangling the heavy-flavor mess: status of the Fermilab-MILC calculation of the $B_{(s)}\to D^{(\ast)}_{(s)}\ellν$ form factors

We present the status of calculations of the form factors of the most relevant heavy-to-heavy and heavy-to-light decay channels. Using seven $N_f = 2+1+1$ HISQ ensembles, with lattice spacings ranging from 0.15 fm down to 0.06 fm, we calculate the form factors of the decays, including correlations among them. More than half of our ensembles feature physical pion masses, and the heavy quarks are simulated at their physical masses using the Wilson-clover action with the Fermilab interpretation. Even though we have recently seen huge qualitative and quantitative leaps in the characterization of heavy-to-heavy decays, these advances have failed to translate into improvements for the inclusive vs exclusive question, or the matter of the Lepton Flavor Universality ratios. In particular, in the $B\to D^{\ast}\ellν$ channel, the current situation of the lattice-QCD form factors is far from clear. Further, the latest lattice-QCD results on the heavy-to-light form factors display unexplained tensions that must urgently be resolved. The work presented here is an attempt to address these issues.

Vaquero, Alejandro [Zaragoza U.]↗
DOE OSTI2025

Log-N/Period Power Offset Sensitivity Calculation

This study focused on the calculation of sensitivity values for Channels 4 and 5 in the Log-N/Period system at the Advanced Test Reactor (ATR). The Log-N/Period system provides an indication of reactor power generation between criticality and full power. This indication is determined by the power offset and the current produced by the detector. Sensitivity values for Channels 4 and 5 were calculated using existing equations, documents, reactor startup data.

42 - ENGINEERING↗
DOE OSTI2025

Systematic many-fermion Hamiltonian input scheme and spectral calculations on quantum computers

We present a novel input scheme for general second-quantized Hamiltonians of relativistic or non-relativistic many-fermion systems. This input scheme incorporates the fermionic anticommutation relations, particle number variations, and respects the symmetries of the Hamiltonian. Based on our input scheme, we propose a hybrid quantum-classical framework for spectral calculations on future quantum hardwares. We provide explicit circuit designs and the associated gate cost. We demonstrate our hybrid framework by solving the low-lying spectra of 42 Ca and 46 Ca. Our input scheme provides new pathways to solving the spectra and dynamics of the relativistic and nonrelativistic many-fermion systems via first-principles approaches.

Hybrid spectral calculation framework↗
DOE OSTI2025

Diffusion of alloying and fission atoms in α-U investigated by first-principles calculation

Stability and diffusion of small defect clusters containing alloying or fission atoms (“F” is used to represent both atoms) in an α-uranium (U) crystal are investigated using first-principles calculations. Here, results indicate that different stable defect structures are observed for an interstitial defect and a cluster (F n V m ) containing different alloying or fission atoms and vacancies. An interstitial defect is found to be most stable either at the center of a pyramid pentahedron or in a 〈010〉 mixed dumbbell configuration and the latter one is reported for the first time. Formation energy of a small F n V m cluster depends on the distribution of vacancies. Based on these stable defects, new migration energies and paths have been explored in this work. Anisotropic migration of an interstitial defect by jumping between nearest centers of two pyramid pentahedrons or jumping between two mixed U-Zr dumbbells has been suggested. Furthermore, a new migration-rotation mechanism has been explored for the first time for a FV 2 cluster, resulting in a 3D diffusion. Finally, the effect of substitutional Pu atoms on the migration of a FV 2 cluster suggests that optimizing the Pu concentration may prohibit the diffusion of some alloying or fission atoms, increasing the safe performance of metallic fuels.

3D diffusion↗
DOE OSTI2025

Ab initio calculations of overlap integrals for μ → e conversion in nuclei

The rate for μ → e conversion in nuclei is set to provide the most stringent test of lepton-flavor symmetry and a window into physics beyond the Standard Model. However, to disentangle new lepton-flavor-violating interactions, in combination with information from μ → ey and μ → 3e, it is critical that uncertainties at each step of the analysis be controlled and fully quantified. In this regard, nuclear response functions related to the coupling to neutrons are notoriously problematic, since they are not directly constrained by experiment. We address these shortcomings by combining ab initio calculations with a recently improved determination of charge distributions from electron scattering by exploiting strong correlations among charge, point-proton, and point-neutron radii and densities. We present overlap integrals for 27 Al, 48 Ca, and 48 Ti including full covariance matrices, allowing, for the first time, for a comprehensive consideration of nuclear structure uncertainties in the interpretation of μ → e experiments.

ab initio calculations↗
DOE OSTI2025

Ab initio calculations of the carbon and oxygen isotopes: Energies, correlations, and superfluid pairing

We perform ab initio nuclear lattice calculations of the neutron-rich carbon and oxygen isotopes using high-fidelity chiral interactions. We find good agreement with the observed binding energies and compute correlations associated with each two-nucleon interaction channel. For the isospin T = 1 channels, we show that the dependence on T z provides a measure of the correlations among the extra neutrons in the neutron-rich nuclei. For the spin-singlet S-wave channel, we observe that any paired neutron interacts with the nuclear core as well as its neutron pair partner, while any unpaired neutron interacts primarily with only the nuclear core. For the other partial waves, the correlations among the extra neutrons grow more slowly and smoothly with the number of neutrons. These general patterns are observed in both the carbon and oxygen isotopes and may be universal features that appear in many neutron-rich nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Real-Space Pseudopotential Method for the Calculation of Third-Row Elements X-ray Photoelectron Spectroscopic Signatures

X-ray photoelectron spectroscopy (XPS) is a powerful characterization technique that unveils subtle chemical environment differences via core– electron binding energy (CEBE) analysis. We extend the development of realspace pseudopotential methods to calculating 1s, 2s, and 2p 3/2 CEBEs of third-row elements (S, P, and Si) within the framework of Kohn–Sham density-functional theory (KS-DFT). The new approach systematically prevents variational collapse and simplifies core-excited orbital selection within dense energy level distributions. However, careful error cancellation analysis is required to achieve accuracy comparable to all-electron methods and experiments. Combined with real-space KS-DFT implementation, this development enables large-scale simulations with both Dirichlet boundary conditions and periodic boundary conditions.

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

Efficient Mixed-Precision Matrix Factorization of the Inverse Overlap Matrix in Electronic Structure Calculations with AI-Hardware and GPUs

In recent years, a new kind of accelerated hardware has gained popularity in the artificial intelligence (AI) community which enables extremely high-performance tensor contractions in reduced precision for deep neural network calculations. In this article, we exploit Nvidia Tensor cores, a prototypical example of such AI-hardware, to develop a mixed precision approach for computing a dense matrix factorization of the inverse overlap matrix in electronic structure theory, S –1 . This factorization of S –1 , written as ZZT = S –1 , is used to transform the general matrix eigenvalue problem into a standard matrix eigenvalue problem. Here we present a mixed precision iterative refinement algorithm where Z is given recursively using matrix–matrix multiplications and can be computed with high performance on Tensor cores. To understand the performance and accuracy of Tensor cores, comparisons are made to GPU-only implementations in single and double precision. Additionally, we propose a nonparametric stopping criteria which is robust in the face of lower precision floating point operations. The algorithm is particularly useful when we have a good initial guess to Z, for example, from previous time steps in quantum-mechanical molecular dynamics simulations or from a previous iteration in a geometry optimization.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Bond Dissociation Energies and Electronic Calculations on the Actinide Halides ThX and UX (X = Cl, Br, I)

Resonant two-photon ionization spectroscopy has been used to locate predissociation thresholds in the spectra of the actinide halides ThX and UX, where X = Cl, Br, and I. These predissociation thresholds are identified as the bond dissociation energies (BDEs) of the molecules. The resulting values show very similar BDEs for the corresponding ThX and UX species, with the thorium molecules being slightly more strongly bound: D 0 (ThCl) = 5.077(6) eV, D 0 (ThBr) = 4.391(4) eV, D 0 (ThI) = 3.537(8) eV, D 0 (UCl) = 4.989(3) eV, D 0 (UBr) = 4.313(3) eV, and D 0 (UI) = 3.449(8) eV. Here, the estimated error limit is given in parentheses in units of the last reported digit. Spinor-based coupled cluster calculations have also been carried out on the halides of this work, including also ThF and UF. Here, the final D 0 values after including contributions due to basis set incompleteness, outer-core-correlation, picture-change, and QED effects are within 0.04 eV of the present experimental values in each case.

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

Emulation of the calculations of final r -process abundance patterns with a neural network

This work explores the construction of a fast emulator for the calculation of the final pattern of nucleosynthesis in the rapid neutron capture process (the r-process). An emulator is built using a feed-forward artificial neural network (ANN). We train the ANN with nuclear data and relative abundance patterns. We take as input the β-decay half-lives and the one-neutron separation energy of the nuclei in the rare-earth region. The output is the final isotopic abundance pattern. In this work, we focus on the nuclear data and abundance patterns in the rare-earth region to reduce the dimension of the input and output space. We show that the ANN can capture the effect of the changes in the nuclear physics inputs on the final r-process abundance pattern in the adopted astrophysical conditions. We employ the deep ensemble method to quantify the prediction uncertainty of the neural network emulator. The emulator achieves a speed-up by a factor of about 20 000 in obtaining a final abundance pattern in the rare-earth region. The emulator may be utilized in statistical analyses such as uncertainty quantification, inverse problems, and sensitivity analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Calculation of rf-induced temporal jitter in ultrafast electron diffraction

A significant contribution to the temporal resolution of an ultrafast electron diffraction (UED) instrument is arrival time jitter caused by amplitude and phase variation of radio-frequency (rf) cavities. In this paper, we present a semianalytical approach for calculating rf-induced temporal jitter from klystron and rf cavity parameters. Our approach allows fast estimation of temporal jitter for MeV-UED beamlines and can serve as a virtual timing tool when shot-to-shot measurements of rf amplitude and phase jitters are available. A simulation study for the SLAC MeV-UED instrument is presented, and the temporal resolution of several beamline configurations is compared.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Melting curves of atomic hydrogen and deuterium calculated using path-integral Monte Carlo

We calculate the melting line of atomic hydrogen and deuterium up to 900 GPa with path-integral Monte Carlo using a machine-learned interatomic potential. We improve upon previous simulations of melting by treating the electrons with reptation quantum Monte Carlo, and by performing solid and liquid simulations using isothermal-isobaric path-integral Monte Carlo. Here, the resulting melting line for atomic hydrogen is higher than previous estimates. There is a small but resolvable decrease in the melting temperature as pressure is increased, which can be attributed to quantum effects.

08 HYDROGEN↗
DOE OSTI2025

Updating Critical Temperature Calculations for Several Secondary High Explosives

In designing experiments where high explosives (HEs) are heated, it is important to have an understanding of where thermal runaway may occur. This determination is often done by using the Frank-Kamenetskii (FK) equation. For several highly-studied HEs, the parameters necessary for FK calculations are typically referenced from decades-old literature, and more recent experimental data have shown that these values require adjustment. For example, some of the legacy parameters dangerously overpredict the critical temperature by tens of degrees relative to more recent observed values. Herein, we first summarize historical results and highlight the importance of insulation effects when estimating critical temperatures for small samples, which is especially relevant to small-scale aging, compatibility, and characterization experiments. Here, we then present updated FK parameters for several common secondary CHNO HEs like RDX, HMX, PETN, and TATB, as well as the first reported values for 2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105). Our updated parameters produce critical-temperature curves that are consistent with both mm-scale differentical scanning calorimetry (DSC) measurements and cm-scale data from the Lawrence Livermore National Laboratory one-dimensional time-to-explosion (ODTX) experiment. In analyzing the critical temperature versus sample size curves derived from the transcendental FK equation, we discovered that these curves are remarkably well described by a simple power-law function with a universal power of 0.142 and a HE-specific scaling factor.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2024

Lattice QCD calculation of the pion distribution amplitude with domain wall fermions at physical pion mass

We present a direct lattice QCD calculation of the x-dependence of the pion distribution amplitude (DA), which is performed using the quasi-DA in large momentum effective theory on a domain-wall fermion ensemble at physical quark masses and spacing a ≈ 0.084 fm. The bare quais-DA matrix elements are renormalized in the hybrid scheme and matched to $\overline{MS}$ with a subtraction of the leading renormalon in the Wilson-line mass. For the first time, we include threshold resummation in the perturbative matching onto the light-cone DA, which resums the large logarithms in the soft gluon limit at next-to-next-to-leading log. The resummed results show controlled scale-variation uncertainty within the range of momentum fraction x ϵ [0.25, 0.75] at the largest pion momentum P z ≈ 1.85 GeV. In addition, we apply the same analysis to quasi-DAs from a highly-improved-staggered-quark ensemble at physical pion mass and a = 0.076 fm. By comparison we find with 2σ confidence level that the DA obtained from chiral fermions is flatter and lower near x = 0.5.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Chemical trends favoring interstitial cluster formation in bcc high-entropy alloys from first-principles calculations

Achieving high strength and ductility is a common goal in structural alloy design. Body-centered cubic high-entropy alloys (HEAs) commonly highlight the conflict between these properties, with stronger alloys being brittle and vice versa. Recent reports suggest interstitial solutes can be used to overcome this trade-off, in some cases providing both strength and ductility enhancements. This effect has been correlated with interstitial cluster formation, although the conditions favoring their formation remain incompletely understood. Using first-principles calculations of solution energies and diffusivities, we provide insights into thermodynamic and kinetic factors favoring interstitial solute clusters. Among C, N and O solutes, O interstitials display most desirable diffusion kinetics. Further, the results highlight the importance of local composition fluctuations in the HEAs to enable the formation of clusters of appreciable size. The results are explained in terms of bonding and distortion trends across solutes and HEA compositions to provide guidelines for alloy design.

Borges, Pedro P P O↗
DOE OSTI2026

Extending SLUSCHI for Automated Diffusion Calculations

We present an extension of the SLUSCHI package (Solid and Liquid in Ultra Small Coexistence with Hovering Interfaces) to enable automated diffusion calculations from first-principles molecular dynamics. While the original SLUSCHI workflow was designed for melting temperature estimation via solid-liquid coexistence, we adapt its input and output handling to isolate the volume search stage and generate one production trajectory suitable for diffusion analysis. Post-processing tools parse VASP outputs, compute mean-square displacements (MSD), and extract tracer diffusivities using the Einstein relation with robust error estimates through block averaging. Diagnostic plots, including MSD curves, running slopes, and velocity autocorrelations, are produced automatically to help identify diffusive regimes. The method has been validated through representative case studies: self-diffusion in Al-Cu liquid alloys, sublattice melting in Li7La3Zr2O12 and Er2O3, interstitial oxygen transport in bcc and fcc Fe, and oxygen diffusivity in Fe-O liquids with variable Si and Al contents. Viscosity and diffusivity are linked through the Stokes-Einstein relation, with composition dependence assessed via simple linear mixing. This capability broadens SLUSCHI from melting-point predictions to transport property evaluation, enabling high-throughput, fully first-principles datasets of diffusion coefficients and viscosities across metals and oxides.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Hydrogen and water interactions with CrMnFeCoNi alloy from density functional theory calculations

High entropy alloys (HEAs) are a promising class of materials with remarkable mechanical and catalytic properties. Among these, the quinary CrMnFeCoNi alloy (also called “Cantor alloy”) has attracted considerable attention given its thermodynamic stability and remarkable mechanical properties under different temperatures. Given that various degradation mechanisms involve multiple contaminants, such as hydrogen and water in hydrogen embrittlement and surface poisoning, respectively, understanding their interactions with the Cantor alloy is critical for its practical applications as structural, nuclear, or hydrogen storage material. In this work, we perform first-principles calculations based on Density Functional Theory (DFT) to investigate such interactions when considering various microstructures, including bulk materials and those containing certain defects, such as grain boundaries, stacking faults, and vacancies. We also employ Global Sensitivity Analysis to identify the importance of different factors in the stability of the impurities. We find that the accuracy of the H formation energy is significantly affected by spin polarization and chemical short-range order. The study also identifies a strong tendency for hydrogen interstitials to segregate to Σ5(210)/[001] symmetric tilt grain boundary, even when H concentrations are high, suggesting that a certain type of grain boundaries acts as H sinks within the alloy. Further, this result is reinforced by the low formation energy of vacancy-hydrogen complexes, which can contain multiple hydrogen atoms. Finally, the surface reactivity analysis reveals that the adsorption energy of oxygen and hydroxyl groups is highly sensitive to the specific metal atom involved in the binding, with a clear preference for chromium atoms, which could have implications for the alloy’s oxidation and corrosion behavior.

36 MATERIALS SCIENCE↗