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

Quantum chemically calculated Abraham parameters for quantifying and predicting polymer hydrophobicity

The leakage and accumulation of plastic in the environment is a significant and growing problem with numerous detrimental impacts and has led to a push toward the design and development of more environmentally benign materials. To this end, we have developed a quantum chemistry-based model for predicting the mobility of polymer materials from molecular structure. Hydrophobicity is used as a surrogate for mobility given that hydrophobic interactions drive much of the partitioning of contaminants in and out of various environmentally relevant compartments. To model polymer hydrophobicity, we adjusted a previously developed Quantum Chemically Calculated Abraham Parameter model to calculate Abraham parameters of small molecules from molecular structure information. The resulting model predicted the octanol-water partition coefficient (K OW ) of polymer repeating units with a root mean square error (RMSE) of 0.48 (log scale). Additionally, the hydrophobicity of high molecular weight polymer materials was captured through solubility parameters and Nile red staining experiments from the literature and predicted with RMSEs of 1.21 (J/cc) 0.5 and 3.42 nm, respectively. Finally, to test the environmental applicability of the model, the relative adsorption capacity of three polymers was predicted and used to unify sorption isotherms across multiple sorbates and polymer sorbents.

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

Exploring isospin symmetry breaking in exotic nuclei: High-precision mass measurement of 23 Si and shell-model calculations of 𝑇 = 5/2 nuclei

Here, we present a high-precision mass measurement of the proton-rich nucleus 23 Si, performed with the LEBIT Penning trap at the Facility for Rare Isotope Beams (FRIB) utilizing the time-of-flight ion cyclotron resonance (TOF-ICR) technique. We determined a mass excess of 23362.9(5.8) keV, which agrees with a recent storage-ring measurement from the experimental Cooler-Storage Ring (CSRe) in Lanzhou but has a factor of 20 improved precision 23 Si is hence the nucleus with the most precisely known mass among all nuclei with an isospin projection of 𝑇 𝑧 = −5/2. We performed shell-model calculations with the USDC and USDCm Hamiltonians to study binding energy differences and Thomas-Ehrmann shifts in mirror systems with an isospin up to 𝑇 = 5/2. Our experimental result and other recently reported masses of neutron-deficient sd-shell nuclei agree well with the theoretical predictions, demonstrating that isospin symmetry breaking in sd-shell nuclei—even at high isospin values—is well described by modern shell-model calculations.

20 ≤ A ≤ 38↗
DOE OSTI2025

Numerical calculation of Coulomb corrections in forward elastic 𝑝↑⁢𝑝 and 𝑝↑⁢𝐴 scattering

The analysis of RHIC hydrogen gas jet target polarimeter measurements of transverse analyzing powers 𝐴 N ⁡(𝑡) in proton-nucleus scattering requires accurate Coulomb corrections to both spin-flip and nonflip amplitudes. These corrections must cover a wide range of nuclear charges 𝑍 and form factor slopes, with flexibility to vary form factors during data fitting. To avoid technically challenging calculations involving a small but finite fictitious photon mass, the Coulomb correction to the nonflip electromagnetic amplitude with an exponential form factor was related to the corresponding correction for the spin-flip amplitude. This approach allows soft photon contributions to all amplitudes, including those with nonexponential form factors, to be calculated in the massless photon limit using only analytical expressions and numerically stable integrals with nonsingular integrands and finite integration limits. In addition, an absorptive correction to the spin-flip electromagnetic amplitude, which plays a critical role in spin effects in forward polarized proton-nucleus scattering, was accurately evaluated.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Error detection on quantum computers improving the accuracy of chemical calculations

The ultimate goal of quantum error correction is to achieve the fault-tolerance threshold beyond which quantum computers can be made arbitrarily accurate. This requires extraordinary resources and engineering efforts. We show that even without achieving full fault-tolerance, quantum error detection is already useful on the current generation of quantum hardware. We demonstrate this by executing an end-to-end chemical calculation for the hydrogen molecule encoded in the [[4, 2, 2]] quantum error-detecting code. The encoded calculation with logical qubits significantly improves the accuracy of the molecular ground-state energy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Dual-cone variational calculation of the two-electron reduced density matrix

The computation of strongly correlated quantum systems is challenging because of its potentially exponential scaling in the number of electron configurations. Variational calculation of the two-electron reduced density matrix (2-RDM) without the many-electron wave function exploits the pairwise nature of the electronic Coulomb interaction to compute a lower bound on the ground-state energy with polynomial computational scaling. Recently, a dual-cone formulation of the variational 2-RDM calculation was shown to generate the ground-state energy, albeit not the 2-RDM, at a substantially reduced computational cost, especially for higher $\textit{N}$-representability conditions such as the T2 constraint. Here we generalize the dual-cone variational 2-RDM method to compute not only the ground-state energy but also the 2-RDM. The central result is that we can compute the 2-RDM from a generalization of the Hellmann-Feynman theorem. Specifically, we prove that in the Lagrangian formulation of the dual-cone optimization the 2-RDM is the Lagrange multiplier. Further, we apply the method to computing the energies and properties of strongly correlated electrons—including atomic charges, electron densities, dipole moments, and orbital occupations—in an illustrative hydrogen chain and the nitrogen-fixation catalyst FeMoco. The dual variational computation of the 2-RDM with T2 or higher N -representability conditions provides a polynomially scaling approach to strongly correlated molecules and materials with significant applications in atomic and molecular and condensed-matter chemistry and physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Convergent close-coupling calculations of electrons scattering on electronically excited molecular hydrogen

We use the adiabatic-nuclei molecular convergent close-coupling method to perform calculations of 0.01–1000 eV electrons scattering on the c 3 Π u , a 3 $Σ$$^{+}_{g}$, B 1 $Σ$$^{+}_{u}$, C 1 Π u , and EF 1 $Σ$$^{+}_{g}$ states of H 2 in the v=0 vibrational level. Elastic, superelastic, ionization, and grand-total cross sections are presented, as well as cross sections for excitation of the n=2–3 singlet and triplet states of H 2 (where n is the atomic-limit principle quantum number). Comparison with available theoretical results is made. Good agreement is found with the recent R-matrix results [J. Phys. B 53, 245203 (2020)] for most of the exchange and dipole-forbidden transitions, but not for the dipole-allowed transitions. Here, the sources of disagreement were found to be an unconverged partial-wave expansion and the utilization of the fixed-nuclei approximation (as opposed to adiabatic-nuclei) in the R-matrix calculations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Searching for optimal THz generation through calculations of the asymmetry of photoelectron momentum distributions by an improved strong-field-approximation method

From the solution of the time-dependent Schrödinger equation (TDSE), it has been shown that in a two-color laser pulse the relative phase that optimizes terahertz (THz) radiation also optimizes the asymmetry of the angular distribution of photoelectrons. Here, in this work, we show that a second-order strong-field approximation can accurately calculate the asymmetry of photoelectrons and thus can locate the phase delay that optimizes THz generation, doing so four orders of magnitude faster than TDSE calculations. We further trace that this is possible because THz emission originates from a free-free radiative transition between the rescattered electron and the target ion, similar to other electron-ion collisions in the laser field, e.g., high-harmonic generation and nonsequential double ionization. Our results pave the way to locate the optimal phase delay for any typical laser pulses for maximal THz generation in the laboratory.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Ab initio calculations of the rate of carrier trapping and release at dopant sites in NaI: Tl beyond the harmonic approximation

We present ab initio calculations of electron capture and release coefficients as a function of temperature for Tl dopants in the widely used scintillator NaI, which is a soft ionic crystal. The modelled capture and release events occur by transitions mediated by multiple phonon absorption and emission between states with significantly different local geometries around the trapping site. We demonstrate that such transitions are not well-described by the normal harmonic approximation to the nuclear dynamics. We avoid the harmonic approximation by numerically solving the vibrational Schrödinger equation for the motion along/of a single effective phonon coordinate. The localized trapped state is not correctly described by semilocal density functionals, so we employ a hybrid functional tuned to reproduce the band gap and lattice constant of the host material. Our calculations, which combine plane-wave and cluster electronic structure methods in a novel way, are in reasonable agreement with available experimental data for detrapping and predict an unusual temperature dependence for trapping rates.

Prange, Micah P.↗
DOE OSTI2020

Calculations of electronic excitation by protons and alpha particles in silicon

This paper presents an ab initio method of calculating the valence electron contribution to the stopping and interaction cross sections of swift ions in matter using linear response calculations of the dielectric function in crystals. We estimate the contribution to the stopping and cross sections by repeated electron capture and stripping of the projectile ion by the medium. We also predict the electronic stopping of protons and alpha particles in silicon and compare to measured values. The predictions generally agree with the data at ion speeds greater than 1 atomic unit. We find that electron capture and stripping is at most a small contribution to the stopping for both ions studied.

electron excitation, silicon, nuclear interactions↗
DOE OSTI2021

Electron correlations in the cubic paramagnetic perovskite Sr ( V , Mn ) O 3 : Results from fully self-consistent self-energy embedding calculations

In this study, we use the thermodynamically consistent and conserving self-energy embedding theory (SEET) to study the spectra of the prototypical undistorted cubic perovskites SrVO 3 and SrMnO 3 . In the strongly correlated metallic SrVO 3 we find that the usual attribution of the satellite peaks at –1.8 eV to Hund or Hubbard physics in the t 2 g orbitals is inconsistent with our calculations. In the strongly correlated insulator SrMnO 3 we recover insulating behavior due to a feedback effect between the strongly correlated orbitals and the weakly correlated environment. Our calculation shows a systematic convergence of spectral features as the space of strongly correlated orbitals is enlarged, paving the way to a systematic parameter-free study of correlated perovskites.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Effects of substitutional Ga on the phonons of δ -phase Pu based on density functional theory calculations

The delta phase of Pu is stabilized by Ga doping, but the mechanism of this stabilization remains an open question. Density functional theory calculations focused on how Ga doping affects the phonons sheds some light on the phonons' contribution to the stabilization. Here, the calculated phonon modes of Ga-doped delta phase Pu fall into two distinct types: localized, high frequency Ga-dominated phonon modes, and Pu-dominated modes at lower frequencies. Increasing the Ga concentration has an effect on the Pu-dominated phonon modes opposite to that of compression: higher-frequency modes soften, and lower-frequency modes stiffen. The latter provides an indication that the stabilization mechanism is not due to a thermodynamic contribution from the phonons. Furthermore, the stiffened phonon modes include candidate modes that describe possible pathways into low-temperature phases, suggesting that doping with Ga could impede such pathways.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Importance of nonuniform Brillouin zone sampling for ab initio Bethe-Salpeter equation calculations of exciton binding energies in crystalline solids

Excitons are prevalent in semiconductors and insulators, and their binding energies are critical for optoelectronic applications. The state-of-the-art method for first-principles calculations of excitons in extended systems is the ab initio GW-Bethe-Salpeter equation (BSE) approach, which can require a fine sampling of reciprocal space to accurately resolve solid-state exciton properties. Here, in this study, we show, for a range of semiconductors and insulators, that the commonly employed approach of uniformly sampling the Brillouin zone can lead to underconverged exciton binding energies, as impractical grid sizes are required to achieve adequate convergence. We further show that nonuniform sampling of the Brillouin zone, focused on the region of reciprocal space where the exciton wave function resides, enables efficient rapid numerical convergence of exciton binding energies at a given level of theory. We propose a well-defined convergence procedure, which can be carried out at relatively low computational cost and which in some cases leads to a correction of previous best theoretical estimates by almost a factor of 2, qualitatively changing the predicted exciton physics. These results call for the adoption of nonuniform sampling methods for ab initio GW-BSE calculations and for revisiting previously computed values for exciton binding energies of many systems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Predicting phase transitions in PbTi⁢O 3 using zentropy through quasiharmonic phonon calculations

According to x-ray diffraction (XRD) measurements, PbTi⁢O 3 undergoes a phase transition from a tetragonal ferroelectric (FE) phase to a cubic paraelectric phase at 763 K. However, x-ray absorption fine-structure (XAFS) measurements indicate that PbTi⁢O 3 is locally tetragonal even after the phase transition. The difference in these results is because XAFS measurements can probe local features of a structure, while XRD averages over such local features. For both measurements to be consistent, PbTi⁢O 3 is macroscopically cubic but locally tetragonal after the phase transition. Despite this, most models, such as the Landau-Ginsburg-Devonshire theory and effective Hamiltonians, are still unable to explain this phenomenon. Moreover, these methods involve model parameters fitted to experimental or theoretical data and do not consider other tetragonal configurations, such as domain walls, to predict the phase transition. In our previous study, we used our zentropy approach to predict the phase transition by considering the tetragonal FE ground-state configuration and the tetragonal 90° and 180° domain wall configurations with their total energies at 0 K. Here, in this paper, the Helmholtz energies of the three configurations are obtained from density functional theory calculations through energy-volume curves and phonon calculations. The predicted phase transition temperature using the meta-GGA 𝑟 2⁢ SCAN and revised multiplicities of configurations is 716 K, showing good agreement with the experimental value of 763 K.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Calculating fluctuations and self-correlations numerically for causal charge diffusion in relativistic heavy-ion collisions

Here, we study the propagation and diffusion of electric charge fluctuations in the Bjorken hydrodynamic model with both white and Catteneo noise using purely numerical methods. We show that a lattice of noise fluctuations is required to fully calculate the two-point correlators of charge. We introduce a numerical procedure to solve the stochastic differential equations that arise from the charge conservation equation on the lattice event by event. We explicitly identify the self-correlation term in the case of Catteneo noise and provide a physical interpretation. We provide a numerical recipe to remove this contribution from the full two-point correlators. Finally, we calculate the balance functions for charged hadrons. By limiting the speed of signal propagation, we observe the expected narrowing of the balance functions after removing the self-correlations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Dirac oscillator: An alternative basis for nuclear structure calculations

The isotropic harmonic oscillator supplemented by a strong spin-orbit interaction has been the cornerstone of nuclear structure since its inception more than seven decades ago. In this paper we introduce—or rather re-introduce—the “Dirac oscillator,” a fully relativistic basis that has all the desired attributes of the ordinary harmonic oscillator while naturally incorporating a strong spin-orbit coupling. To assess the power and flexibility of the Dirac oscillator basis in the solution of nuclear structure problems within the framework of covariant density-functional theory. Here, self-consistent calculations of binding energies and ground-state densities for a selected set of doubly magic nuclei are performed using the Dirac oscillator basis and are then compared against results obtained with the often-used Runge-Kutta method. Results obtained using the Dirac oscillator basis reproduce with high accuracy those derived using the Runge-Kutta method and suggest a clear path for a generalization to systems with axial symmetry. Although the harmonic oscillator with spin-orbit corrections has been the staple of the nuclear shell model since the beginning, the Dirac oscillator is practically unknown among the nuclear physics community. In this paper we illustrate the power and flexibility of the Dirac oscillator and suggest extensions to the study of systems without spherical symmetry, as required in constrained calculations of nuclear excitations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Skyrme pseudopotentials at next-to-next-to-leading order: Construction of local densities and first symmetry-breaking calculations

There is an ongoing quest to improve on the spectroscopic quality of nuclear energy density functionals (EDFs) of the Skyrme type through extensions of its traditional form. One direction for such activities is the inclusion of terms of higher order in gradients in the EDF. We report on exploratory symmetry-breaking calculations performed for an extension of the Skyrme EDF that includes central terms with four gradients at next-to-next-to-leading order (N2LO) and for which the high-quality parametrization SN2LO1 has been constructed recently. Up to now, the investigation of such functionals with higher-order terms was limited to infinite matter and spherically symmetric configurations of singly and doubly magic nuclei. We address here nuclei and phenomena that require us to consider axial and nonaxial deformation, both for reflection-symmetric and also reflection-asymmetric shapes, as well as the breaking of time-reversal invariance. Achieving these calculations demanded a number of formal developments. These all resulted from the formulation of the N2LO EDF requiring the introduction of new local densities with additional gradients that are not present in the EDF at NLO. Their choice is not unique, but can differ in the way the gradients are coupled. While designing a numerical implementation of N2LO EDFs in Cartesian three-dimensional coordinate-space representation, we have developed a novel definition and a new unifying notation for normal and pair densities that contain gradients at arbitrary order. Besides having mnemonic advantages, the new notation allows for the easy identification of redundancies and reducibilities in a given set of local densities, and the new definition makes it straightforward to construct densities that automatically adopt the symmetries of the many-body state they are constructed from. The resulting scheme resolves several issues with some of the choices that have been made for local densities in the past, in particular when breaking time-reversal symmetry. Guided by general practical considerations, we propose an alternative form of the N2LO contribution to the Skyrme EDF that is built from a different set of densities. It has exactly the same physics content, but is much more efficient to handle in formal discussions and, compared to the original formulation, leads to a substantial reduction of computational cost and memory requirements in deformed codes. As representative examples for the performance of SN2LO1, we have chosen the ground states of even-even Kr and Nd isotopes, the fission barrier of 240 Pu as well as the superdeformed rotational band of 194 Hg. Overall, for the nuclei and phenomena studied here, the SN2LO1 parametrization does not yet present a systematic improvement over standard NLO parametrizations. This finding calls for improved fit protocols that better discriminate between NLO and N2LO terms and better exploit the unique features of the additional degrees of freedom offered by the latter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Converged ab initio calculations of heavy nuclei

We propose a novel storage scheme for three-nucleon (3N) interaction matrix elements relevant for the normal-ordered two-body approximation used extensively in ab initio calculations of atomic nuclei. This scheme reduces the required memory by approximately two orders of magnitude, which allows the generation of 3N interaction matrix elements with the standard truncation of E 3max =28, well beyond the previous limit of 18. We demonstrate that this is sufficient to obtain the ground-state energy of 132 Sn converged to within a few MeV with respect to the E 3max truncation. In addition, we study the asymptotic convergence behavior and perform extrapolations to the un-truncated limit. Finally, we investigate the impact of truncations made when evolving free-space 3N interactions with the similarity renormalization group. We find that the contribution of blocks with angular momentum J rel > 9/2 to the ground-state energy is dominated by a basis-truncation artifact, which vanishes in the large-space limit, so these computationally expensive components can be neglected. For the two sets of nuclear interactions employed in this work, the resulting binding energy of 132 Sn agrees with the experimental value within theoretical uncertainties. This work enables converged ab initio calculations of heavy nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Ab initio coupled-cluster calculations of ground and dipole excited states in He 8

We perform coupled-cluster calculations of ground- and dipole excited-state properties of the 8 He halo nucleus with nucleon-nucleon and three-nucleon interactions from chiral effective field theory, both with and without explicit delta degrees of freedom. By increasing the precision in our coupled-cluster calculations via the inclusion of leading-order three-particle three-hole excitations in the cluster operator, we obtain a ground-state energy and a charge radius that are consistent with experiment, albeit with a slight underbinding. We also investigate the excited states induced by the electric dipole operator and present a discussion on the Thomas-Reiche-Kuhn and cluster sum rules. Lastly, we compute the electric dipole polarizability, providing a theoretical benchmark for future experimental determinations that will study this exotic nucleus.

6 ≤ A ≤ 19↗