Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “calculator”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

EPR-based uncertainty validation of the calculated external doses for population exposed in the urals region

Tooth enamel Electron Paramagnetic Resonance (EPR) spectroscopy was used as a method for external dosimetry in the territories contaminated in the 1950s by PA ‘Mayak’ (Urals region) to validate the mean dose estimates predicted by the Techa River Dosimetry System (TRDS). The purpose of this study is to validate the uncertainties of TRDS doses. Ninety percent confidence intervals (90% confidence interval, CI) of dose estimated with both methods were compared for 220 people. All data were grouped according to the width of 90%CI, viz.: (1) 90%CI of EPR-based dose ≤ 90%CI of TRDS prediction (38 cases); (2) 90%CI of EPR-based dose > 90%CI of TRDS prediction (182 cases). About 91% of 90%CIs overlap. In group 1, 100% cases overlap. In group 2, 80% of the cases were non-contradictive (the calculated 90%CI is completely within the measured one). Interval comparison of doses predicted retrospectively and estimated based on individual measurements are non-contradictory and demonstrate a good agreement.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2026

Ab initio many-fermion structure calculations on a quantum computer

To overcome the limitations of existing algorithms for solving self-bound quantum many-body problems—such as those encountered in nuclear and particle physics—that access only a restricted subset of energy levels and provide limited structural information, we introduce and demonstrate a novel quantum-classical approach capable of resolving the complete bound-state spectrum. This method also provides the total angular momentum 𝐽 associated with each eigenstate. Here, our approach is based on expressing the Hamiltonian in second-quantized form within a novel input model combined with a scan scheme, enabling broad applicability to configuration-interaction calculations across diverse fields. We apply this hybrid method to compute, for the first time, the bound-state spectrum together with corresponding 𝐽 values of 20 O using a realistic strong-interaction Hamiltonian. Our approach applies to hadron spectra and 𝐽 values solved in the relativistic basis light-front quantization approach.

Du, Weijie [Chinese Academy of Sciences (CAS), Lan↗
DOE OSTI2025

Tunable energy landscape of screw dislocation cores by compositional fluctuations in bcc high-entropy alloys from first-principles calculations

The energy landscape of screw dislocation cores plays a central role in dislocation-mediated deformation mechanisms in body-centered cubic (bcc) metals. In bcc high-entropy alloys (HEAs), this energy landscape is modulated by local compositional fluctuations, which has important implications for deformation processes in these materials. Through first-principles calculations, this study investigates high-symmetry screw dislocation core structures in NbTaMoW and NbTaTiHf bcc HEAs. The results show that alloying group IV transition metals lead to large local lattice distortions at dislocation cores, which is demonstrated to be an important factor governing fluctuations in core configurations along a dislocation line. Importantly, group IV elements near the core induce features in the energy landscape that are exclusive for HEAs, specifically lowering the energy of core configurations that are unstable in elemental bcc metals. A combined influence of these chemical effects with crystallographic details enables the activation of glide planes, a feature that has been linked to ductility improvements in bcc HEAs. These findings provide new insights into the atomic-scale mechanisms underlying dislocation mobility in bcc HEAs, offering a pathway for designing materials with tailored mechanical properties.

Borges, Pedro P P O↗
DOE OSTI2026

Quantum Monte Carlo calculation of {delta}C in the superallowed beta decay of 10C

We perform an ab initio quantum Monte Carlo calculation of the isospin-symmetry-breaking correction δC to the superallowed β decay of ¹⁰C. Using both phenomenological and chiral nuclear interactions, we evaluate the Fermi matrix element and quantify its deviation from the canonical √2 value. The resulting δC values lie in the range ≈ 0.15–0.25% and are consistent, within sizable uncertainties (approximately 34%–65% relative), across Hamiltonians, indicating no statistically significant dependence on the choice of nuclear interaction. The extracted values of Vud are also found to be compatible with current determinations within these uncertainties.

Piarulli, M↗
DOE OSTI2025

Self-consistent microscopic calculations for electron captures on nuclei in core-collapse supernovae

Calculations for electron capture rates on nuclei with atomic numbers between 𝑍 = 20 and 𝑍 = 52 are performed in a self-consistent finite-temperature covariant energy density functional theory within the relativistic quasiparticle random-phase approximation. Electron captures on these nuclei contribute most to reducing the electron fraction during the collapse phase of core-collapse supernovae. The rates include contributions from allowed (Gamow-Teller) and first-forbidden (FF) transitions, and it is shown that the latter become dominant at high stellar densities and temperatures. Temperature-dependent effects such as Pauli unblocking and transitions from thermally excited states are also included. The new rates are implemented in a spherically symmetric one-dimensional simulation of the core-collapse phase. The results indicate that the increase in electron capture rates, due to inclusion of FF transitions, leads to reductions of the electron fraction at nuclear saturation density, the peak neutrino luminosity, and enclosed mass at core bounce. The new rates reaffirm that the most relevant nuclei for the deleptonization situate around the 𝑁 = 50 and 82 shell closures, but, compared to previous simulations, nuclei are less proton rich. Here, the new rates developed in this work are available, and will be of benefit to improve the accuracy of multidimensional supernova simulations.

Electron & muon capture↗
DOE OSTI2026

Machine-learned quantum molecular dynamics calculations of warm dense equation of state and ionic transport coefficients of deuterated water

White dwarf models require accurate equations of state and ionic transport coefficients in the warm dense matter regime, where kinetic theory models and tabulated equations of state are often inaccurate. In this work, spectral-partitioned density functional theory and machine-learned interatomic potentials are combined to perform large-scale, first-principles quantum molecular dynamics simulations of deuterated water (D 2 O) near the principal Hugoniot. This approach retains Kohn-Sham accuracy while achieving orders-of-magnitude speedup, yielding converged equation of state and transport properties over a broad pressure and temperature range. The results reveal the thermodynamic conditions under which ionic transport models for interdiffusivity and shear viscosity converge and identify those in closest agreement with density functional theory benchmarks at temperatures in the warm dense matter regime. The present framework extends first-principles transport calculations to higher temperatures than previously achieved, and provides an efficient, scalable, and general approach for studying transport properties in complex multicomponent mixtures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Fast and accurate calculation of EXAFS Debye-Waller factors in U⁢O2 using the dynamical matrix method

Theoretical modeling of bonding dynamics in metal oxides is required for predicting their thermal conductivity, catalytic activity, and mechanical properties. A primary challenge is the scarcity of experimental methods for validating theoretical predictions of these atomic-scale dynamics. This work presents a workflow that uses experimental extended x-ray absorption fine structure (EXAFS) data collected at high temperatures to validate an interatomic force field for uranium dioxide (UO2), an important model material. The validated force field is then used to drive computationally intensive molecular dynamics (MD) simulations and as input for the much faster dynamical matrix Debye-Waller (DMDW) method. The predicted values of the Debye-Waller factors from the DMDW calculations are in good agreement with those obtained from the MD simulations, with residual pair-specific differences attributable to quantum zero-point motion at low temperatures and lattice anharmonicity at high temperatures. We further show that theoretical EXAFS spectra constructed directly from DMDW-derived Debye-Waller factors reproduce the experimental data (at relatively low temperatures) with accuracy comparable to full MD-EXAFS, providing an additional validation of the choice of the potential. This study establishes a validated, rapid computational pathway for modeling bond dynamics, naturally incorporating quantum nuclear\\\\r\\\\nstatistics absent in classical simulations, which are essential for the mechanistic understanding of complex oxide materials.

58 GEOSCIENCES↗
DOE OSTI2026

Quantum Monte Carlo calculation of 𝛿 NS in 10 C using an effective field theory approach

Here, we compute radiative corrections to the superallowed 𝛽 decay of 10 C in an effective field theory approach using nuclear matrix elements obtained from quantum Monte Carlo calculations. These corrections are an important ingredient in the extraction of the Cabibbo-Kobayashi-Masakawa quark mixing matrix element 𝑉 𝑢⁢𝑑 , and the role of this work is to illuminate the uncertainties arising from nuclear structure. Our results provide good agreement with both the traditional extraction of 𝑉 𝑢⁢𝑑 , as well as with a more recent evaluation performed using the no-core shell model and a dispersion formalism. The dominant uncertainty in this approach is the presence of two unknown low-energy constants that enter into the relevant nuclear matrix elements. Future determinations of these low-energy constants—either from QCD or modeling them with two nucleon amplitudes—would improve the precision of the extraction in this formalism.

beta decay↗
DOE OSTI2026

From closed shells to open shells: Coupled-cluster calculations of atomic nuclei

Coupled-cluster theory is a powerful tool for first-principles calculations of atomic nuclei, enabling accurate predictions of nuclear observables across the Segrè chart. While coupled-cluster computations are especially efficient at shell closures, extensions have been developed to tackle open-shell nuclei, by exploiting the equation-of-motion method or by expanding the coupled-cluster wave function on top of a symmetry-breaking (either deformed or superfluid) reference state. In this study, we provide a comprehensive comparison of these different formulations applied to the calcium and nickel isotopes using nuclear two-and three-body interactions from chiral effective field theory. Here, based on ground-state energies, two-neutron separation energies, and two-neutron shell gaps, different coupled-cluster computations—based on symmetry-broken reference states and equationof-motion techniques— offer consistent descriptions of bulk properties across medium-mass isotopic chains.

Marino, Francesco [Johannes Gutenberg-Universität ↗
DOE OSTI2023

Elastic scattering and rotational excitation of H 2 by electron impact: Convergent close-coupling calculations

We apply the adiabatic-nuclei molecular convergent close-coupling (MCCC) method to the study of elastic scattering and rotational excitation of H 2 by 0.01–20-eV electrons. Integral cross sections are presented for all rotational transitions with |Δ⁢N| ≤ 2 and N i = 0–31 within the v = 0 vibrational level and differential cross sections for a selection of transitions. Agreement with the available measurements and previous calculations is mostly excellent, depending on the transition and incident energy. Finally, we suggest possible reasons, and argue for the accuracy of the MCCC data.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Image charge interaction correction in charged-defect calculations

Charged-defect calculation using a periodic supercell is a significant class of problems in solid state physics. However, the finite supercell size induces an undesirable long-range image charge Coulomb interaction. Although a variety of methods have been proposed to eliminate such image Coulomb interaction, most of the previous schemes are based on a rough approximation of the defect charge screening. In this work, we present a rigorous derivation of the image charge interaction with a defect screening model where the use of a bulk macroscopic dielectric constant can be avoided. Here, we have verified this approach in comparison with a widely used approach for 12 different defects. Our correction scheme offers a much faster convergence concerning the supercell size for cases with considerable image charge interactions. In those cases, we also found that the nonlinear dielectric screening might play an important role. Our proposed defect screening model will also shed new light on understanding the defect screening properties and can be applied to other defect systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Fast and robust all-electron density functional theory calculations in solids using orthogonalized enriched finite elements

Here, we present a computationally efficient approach to perform systematically convergent real-space all-electron Kohn-Sham density functional theory calculations for solids using an enriched finite element (FE) basis. The enriched FE basis is constructed by augmenting the classical FE basis with atom-centered numerical basis functions, comprising of atomic solutions to the Kohn-Sham problem. Notably, to improve the conditioning, we orthogonalize the enrichment functions with respect to the classical FE basis, without sacrificing the locality of the resultant basis. In addition to improved conditioning, this orthogonalization procedure also renders the overlap matrix block diagonal, greatly simplifying its inversion. Subsequently, we use a Chebyshev polynomial based filtering technique to efficiently compute the occupied eigenspace in each self-consistent field iteration. We demonstrate the accuracy and efficiency of the proposed approach on periodic unit cells and supercells. The benchmark studies show a staggering 130× speedup of the orthogonalized enriched FE basis over the classical FE basis. We also present a comparison of the orthogonalized enriched FE basis with the linearized augmented plane-wave + local orbitals basis, both in terms of accuracy and efficiency. Notably, we demonstrate that the orthogonalized enriched FE basis can handle large system sizes of ~10 000 electrons. Finally, we observe good parallel scalability of our implementation with 92% efficiency at 22× speedup for a system with 620 electrons

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

First-principles calculation of ballistic current from electron-hole interaction

The bulk photovoltaic effect (BPVE) has attracted increasing interest due to its potential to overcome the efficiency limit of traditional photovoltaics, and much effort has been devoted to understanding its underlying physics. However, previous work has shown that theoretical models of the shift current and the phonon-assisted ballistic current in real materials do not fully account for the experimental BPVE photocurrent, so other mechanisms should be investigated in order to obtain a complete picture of BPVE. As such, in this paper, we demonstrate two approaches that enable the ab initio calculation of the ballistic current originating from the electron-hole interaction in semiconductors. Using BaTiO 3 and monolayer MoS 2 as two examples, we show clearly that for them the asymmetric scattering from an electron-hole interaction is less appreciable than that from an electron-phonon interaction, indicating more scattering processes need to be included to further improve the BPVE theory. Moreover, we provide a first-principles approach for material prediction and design in order to search for materials with a larger ballistic current due to electron-hole interactions.

14 SOLAR ENERGY↗
DOE OSTI2021

Type-II Dirac cones and electron-phonon interaction in monolayer biphenylene from first-principles calculations

In this work, we report a first-principles investigation of electronic structure, topological bands, and electron-phonon interactions in metallic biphenylene sheets. Biphenylene is a recently synthesized sp 2 -bonded carbon allotrope. We find coupling of electrons at the Fermi surface to very high frequency carbon-derived phonons, analogous to superconducting MgB 2 . This leads to low-temperature weak coupling superconductivity due to an unusual combination of exceptionally large logarithmically averaged phonon frequency ω log =1369 K and moderate electron-phonon coupling. The electronic structure shows a two-band Fermi surface dominated by C p z orbitals and a pair of type-II tilted Dirac cones along the Γ-Y line at the Brillouin zone boundary. Berry curvature and edge-state calculations show that monolayer biphenylene is a two-dimensional $\mathbb{Z}$ 2 topological material. Thus, monolayer biphenylene is predicted to be a topological superconductor based on C p orbitals and high-frequency phonons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Electronic and vibrational properties of bulk Cr 2 Ge 2 Te 6 from first-principles calculations

Here we use density functional theory based ab initio calculations to investigate the structural, vibrational, magnetic, and electronic properties of the layered ferromagnet Cr 2 Ge 2 Te 6 (CGT) that has attracted attention for potential spintronic applications. We optimize the structure for a fixed c/a ratio. Our results are in excellent agreement with experimental data on structure, phonons, and electronic properties. The use of van der Waals interactions and relativistic spin-orbit coupling yields accurate lattice constants and interlayer distances. Simultaneously, the adopted theoretical methods lead to a rigorous description of the vibrational normal modes, as well as the valence bands, that are in excellent agreement with Raman spectroscopy and angle-resolved photoemission spectroscopy spectra, respectively. The magnetic moment is slightly overestimated, and the magnetic anisotropy has the correct sign but is greater in magnitude than that in experiment. We find that inclusion of the on-site Coulomb repulsion on Cr−3d orbitals worsens agreement with experiment, in contrast to previous studies. Our results provide a path toward ab initio analysis of magnetic heterostructures based on CGT.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Thermodynamic and electron transport properties of Ca 3 Ru 2 O 7 from first-principles phonon calculations and Boltzmann transport theory

Here, this work demonstrates a first-principles-based approach to obtaining finite temperature thermal and electronic transport properties which can be employed to model and understand mesoscale structural evolution during electronic, magnetic, and structural phase transitions. A computationally tractable model was introduced to estimate electron relaxation time and its temperature dependence. The model is applied to Ca 3 Ru 2 O 7 with a focus on understanding its electrical resistivity across the electronic phase transition at 48 K. A quasiharmonic phonon approach to the lattice vibrations was employed to account for thermal expansion while the Boltzmann transport theory including spin-orbit coupling was used to calculate the electron-transport properties, including the temperature dependence of electrical conductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Calculating n-point charge correlations in evolving systems

Background: In dynamic systems such as heavy-ion collisions, charge susceptibilities and local charge correlations change with time. These changes are accompanied by nonlocal correlations, which spread diffusively with time and are constrained by local charge conservation. Such correlations have been measured at the Relativistic Heavy-Ion Collider (RHIC) and have been linked to the chemical evolution of the matter and possible phase separation or critical phenomena. The evolution of such correlations has been modeled for two-point correlators superimposed onto hydrodynamics, but not for three- or four-point correlators. These higher-order correlators represent an essential basis for calculating the kurtosis and skewness of charge distributions measured at RHIC. For short-range correlations, such as those between charges on the same particle, and the associated charge-balance correlations a theoretical formalism was lacking for handling higher-order correlations. In particular, one needed to understand how higher-order charge fluctuations was split onto lower numbers of particles. For example, correlations of order three have contributions with charge all on one particle, split onto two particle, or split onto three. By assuming local chemical equilibrium of short-range correlations, the sources of the various correlators were found to be uniquely determined. The relevant theoretical foundation is presented here, including a diagrammatic technique to correctly account for all the possible terms to higher-order correlations. Here, we found a consistent theoretical treatment and its applicability, viability, and tractability are discussed. The formalism derived here enables realistic and quantitative modeling of the three- and four-point charge correlations necessary for understanding measurements of the kurtosis and skewness of charge distributions at RHIC. These observables have been promoted as signals of phase transitions or measures of chemical evolution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Microscopic calculations of nuclear level densities with the Lanczos method

A new method for computing the density of states in nuclei making use of an extrapolated form of the tridiagonal matrix obtained from the Lanczos method is presented. It will be shown that the global, average properties of the entire Lanczos matrix can be predicted from just four Lanczos iterations. The extrapolated Lanczos matrix (ELM) approach provides for an accurate computation of the density of states described within the configuration space, which, in some cases, is sufficient to accurately calculate the density of states at, or near, the neutron separation energy. Comparisons between theory and experiment are shown for $^{57}Fe$, $^{74}Ge$, and $^6Ge$. In addition, we show results for the J -dependence of moments and the level density for these three nuclei.

Shell model↗