Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CALCULATION METHODS”

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 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Persistent Protein Motions in a Rugged Energy Landscape Revealed by Normal Mode Ensemble Analysis

Proteins are allosteric machines that couple motions at distinct, often distant, sites to control biological function. Low-frequency structural vibrations are a mechanism of this long-distance connection and are often used computationally to predict correlations, but experimentally identifying the vibrations associated with specific motions has proved challenging. Spectroscopy is an ideal tool to explore these excitations, but measurements have been largely unable to identify important frequency bands. The result is at odds with some previous calculations and raises the question what methods could successfully characterize protein structural vibrations. Here we show the lack of spectral structure arises in part from the variations in protein structure as the protein samples the energy landscape. However, by averaging over the energy landscape as sampled using an aggregate 18.5 μs of all-atom molecular dynamics simulation of hen egg white lysozyme and normal-mode analyses, we find vibrations with large overlap with functional displacements are surprisingly concentrated in narrow frequency bands. These bands are not apparent in either the ensemble averaged vibrational density of states or isotropic absorption. However, in the case of the ensemble averaged anisotropic absorption, there is persistent spectral structure and overlap between this structure and the functional displacement frequency bands. We systematically lay out heuristics for calculating the spectra robustly, including the need for statistical sampling of the protein and inclusion of adequate water in the spectral calculation. The results show the congested spectrum of these complex molecules obscures important frequency bands associated with function and reveal a method to overcome this congestion by combining structurally sensitive spectroscopy with robust normal mode ensemble analysis.

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

Non-Hermitian quantum mechanics approach for extracting and emulating continuum physics based on bound-state-like calculations: Detailed description

Here, this work applies a reduced basis method to study the continuum physics of a finite quantum system—either few or many-body. Specifically, I develop reduced-order models, or emulators, for the underlying inhomogeneous Schrödinger equation and train the emulators against the equation's bound-state-like solutions at complex energies. The emulators rapidly and accurately interpolate and extrapolate the matrix elements of the Hamiltonian resolvent operator (Green's function) across a parameter space that includes both complex energy and other real-valued physical inputs in the Schrödinger equation. The spectra, discretized and compressed as the result of emulation, and the associated resolvent matrix elements (or amplitudes), have the defining characteristics of non-Hermitian quantum mechanics calculations, featuring complex eigenenergies with negative imaginary parts and branch cuts moved below the real axis in the complex energy plane. Therefore, one now has a method that extracts continuum physics from bound-state-like calculations and emulates those extractions in the input parameter space. Building on a prior Letter [Zhang, Phys. Rev. Lett. 135, 242501 (2025)], this article provides the full theoretical details, a comprehensive analysis of the method's performance, and a brief discussion of how it can be coupled with existing continuum approaches to perform emulations in their input parameter spaces.

ab initio calculations↗
DOE OSTI2021

Elaborated thermochemical treatment of HF, CO, N 2 , and H 2 O: Insight into HEAT and its extensions

Here, empirical, highly accurate non-relativistic electronic total atomization energies (eTAEs) are established by combining experimental or computationally converged treatments of the nuclear motion and relativistic contributions with the total atomization energies of HF, CO, N 2 , and H 2 O obtained from the Active Thermochemical Tables. These eTAEs, which have estimated (2σ) uncertainties of less than 10 cm -1 (0.12 kJ mol -1 ), form the basis for an analysis of high-level ab initio quantum chemical calculations that aim at reproducing these eTAEs for the title molecules. The results are then employed to analyze the performance of the high-accuracy extrapolated ab initio thermochemistry, or High-Accuracy Extrapolated Ab Initio Thermochemistry (HEAT), family of theoretical methods. The method known as HEAT-345(Q), in particular, is found to benefit from fortuitous error cancellation between its treatment of the zero-point energy, extrapolation errors in the Hartree-Fock and coupled cluster contributions, neglect of post-(T) core-correlation, and the basis-set error involved in higher-level correlation corrections. In addition to shedding light on a longstanding curiosity of the HEAT protocol—where the cheapest HEAT-345(Q) performs comparably to the theoretically more complete HEAT-456QP procedure—this study lays the foundation for extended HEAT variants that offer substantial improvements in accuracy relative to the established approaches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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

Convergence acceleration of Monte Carlo many-body perturbation methods by direct sampling

In the Monte Carlo many-body perturbation (MC-MP) method, the conventional correlation-correction formula, which is a long sum of products of low-dimensional integrals, is first recast into a short sum of high-dimensional integrals over electron-pair and imaginary-time coordinates. These high-dimensional integrals are then evaluated by the Monte Carlo method with random coordinates generated by the Metropolis–Hasting algorithm according to a suitable distribution. The latter algorithm, while advantageous in its ability to sample nearly any distribution, introduces autocorrelation in sampled coordinates, which in turn increases the statistical uncertainty of the integrals and thus the computational cost. It also involves wasteful rejected moves and an initial “burn-in” step as well as displays hysteresis. Here, an algorithm is proposed that directly produces a random sequence of electron-pair coordinates for the same distribution used in the MC-MP method, which is free from autocorrelation, rejected moves, a burn-in step, or hysteresis. Furthermore, this direct-sampling algorithm is shown to accelerate second- (MC-MP2) and third-order Monte Carlo many-body perturbation (MC-MP3) calculations by up to 222% and 38%, respectively.

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

Nucleon-Pion Spectroscopy from Sparsened Correlators

Neutrino oscillation experiments require accurate reconstructions of neutrino energies, which depend in part on a theoretical understanding of the axial $N\rightarrow\Delta$ transition form factors. A lattice QCD study of this transition will require construction of all hadronic states with energies up to $m_\Delta$, which at the physical point includes $N\pi$ and $N\pi\pi$. Building interpolating operators from sparse grids at the source and sink is a versatile method that allows construction of a wide range of diagram topologies and has successfully been used in other multi-hadron calculations. We will discuss application of this method to nucleon-pion systems and present preliminary results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Assessment of DFT functionals for a minimal nitrogenase [Fe(SH) 4 H]- model employing state-of-the-art ab initio methods

We have designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S. We show that the energy difference between these two isomers (ΔE) is hard to estimate with quantum-mechanical (QM) methods. For example, different density functional theory (DFT) methods give ΔE estimates that vary by almost 140 kJ/mol, mainly depending on the amount of exact Hartree–Fock included (0%–54%). The model is so small that it can be treated by many high-level QM methods, including coupled-cluster (CC) and multiconfigurational perturbation theory approaches. With extrapolated CC series (up to fully connected coupled-cluster calculations with singles, doubles, and triples) and semistochastic heat-bath configuration interaction methods, we obtain results that seem to be converged to full configuration interaction results within 5 kJ/mol. Our best result for ΔE is 101 kJ/mol. With this reference, we show that M06 and B3LYP-D3 give the best results among 35 DFT methods tested for this system. Brueckner doubles coupled cluster with perturbaitve triples seems to be the most accurate coupled-cluster approach with approximate triples. CCSD(T) with Kohn–Sham orbitals gives results within 4–11 kJ/mol of the extrapolated CC results, depending on the DFT method. Single-reference CC calculations seem to be reasonably accurate (giving an error of ~5 kJ/mol compared to multireference methods), even if the D1 diagnostic is quite high (0.25) for one of the two isomers.

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

New Hybrid Model for Evaluating the Frequency-Dependent Leakage Inductance of a Variable Inductance Transformer (VIT)

Skin and proximity effects can cause a significant drop in the effective leakage inductance of a transformer when the operating frequency is increased. Although the magnetic image method-based double-2-D model can calculate the lowfrequency leakage inductance with sufficient accuracy, it is inherently a frequency-independent model. While Dowell’s 1-D model uses frequency-dependent relations to account for both skin and proximity effects, its accuracy is severely affected by the assumed winding geometry. In this paper, a hybrid model is proposed that uses superposition to combine a modified Dowell’s model with the double-2-D model. The proposed model is investigated on a variable inductance transformer (VIT)—a partially-filled transformer whose leakage inductance can be varied by moving one of the windings mechanically. The frequency-dependent leakage inductances of the VIT evaluated using the hybrid model are in excellent agreement with the corresponding finite element method (FEM) simulated and experimentally measured values, thereby validating the proposed hybrid model.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Reduced-basis method for few-body bound-state emulation

Recent advances in both theoretical and computational methods have enabled large-scale, precision calculations of the properties of atomic nuclei. With the growing complexity of modern nuclear theory, however, also comes the need for novel methods to perform systematic studies and quantify the uncertainties of models when confronted with experimental data. Here, this study presents an application of such an approach, the reduced basis method, to substantially lower computational costs by constructing a significantly smaller Hamiltonian subspace informed by previous solutions. Our method shows comparable efficiency and accuracy to other dimensionality reduction techniques on an artificial three-body bound system while providing a richer representation of physical information in its projection and training subspace. This methodological advancement can be applied in other contexts and has the potential to greatly improve our ability to systematically explore theoretical models and thus enhance our understanding of the fundamental properties of nuclear systems.

cluster models↗
DOE OSTI2020

Feasibility of using thermal response methods for nonintrusive compressed air flow measurement

Robust, low-cost nonintrusive flow meters are of interest in many industries. Particularly, a reliable nonintrusive flow measurement for the diagnosis of air leaks in compressed air systems is desirable. Measurement of the air flow due to leaks in the system ensures an accurate estimation of potential cost and energy savings. This study evaluates a novel method of using thermal responses to nonintrusively measure leakage rates in compressed air lines. The method uses heat and the resulting thermal response to calculate the flow rate inside the compressed air line. Compared to the current methods for flow measurement, this method can simplify flow measurement while decreasing the sensitivity to errors when measuring flow rates. In this study, the methodology of the proposed method is explained along with the potential advantages to the design. Two approaches are evaluated: a dynamic step response and sinusoidal frequency response. Simulated tests evaluate the feasibility of the proposed methods, followed by experiments that validate the simulation results. A clear correlation between the thermal step response and the flow rate indicate viability of the proposed method in simulation. Experimental results yielded similar results, confirming the validity of the proposed method. The results of a field test in an industrial environment demonstrate the capability of the approach to other flow rate measurement techniques.

42 ENGINEERING↗
DOE OSTI2025

Generalized approach for rapid entropy calculation of liquids and solids

We build a comprehensive methodology for the fast computation of entropy across both solid and liquid phases. The proposed method utilizes a single trajectory of molecular dynamics (MD) to facilitate the calculation of entropy, which is composed of three components. The electronic entropy is determined through the temporal average acquired from density functional theory MD simulations. The vibrational entropy, typically the predominant contributor to the total entropy, even within the liquid state, is evaluated by computing the phonon density of states via the velocity autocorrelation function. The most arduous component to quantify, the configurational entropy, is assessed by probability analysis of the local structural arrangement and atomic distribution. We illustrate, through a variety of examples, that this method is both a versatile and valid technique for characterizing the thermodynamic states of both solids and liquids. Furthermore, this method is employed to expedite the calculation of melting temperatures, demonstrating its practical utility in computational thermodynamics.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Water Catalysis of the Reaction of Methanol with OH Radical in the Atmosphere is Negligible

Abstract Faced with the contradictory results of two recent experimental studies [Jara‐Toro et al., Angew. Chem. Int. Ed . 2017 , 56 , 2166 and Chao et al., Angew. Chem. Int. Ed . 2019 , 58 , 5013] of the possible catalytic effect of water vapor on CH 3 OH + OH reaction, we report calculations that corroborate the conclusion made by Chao et al. and extend the rate constant evaluation down to 200 K. The rate constants of the CH 3 OH + OH reaction catalyzed by a water molecule are computed as functions of temperature and relative humidity using high‐level electronic structure and kinetics calculations. The Wuhan–Minnesota Scaling (WMS) method is used to provide accurate energetics to benchmark a density functional for direct dynamics. Both high‐frequency and low‐frequency anharmonicities are included. Variational and tunneling effects are treated by canonical variational transition state theory with multidimensional small‐curvature tunneling. And, most significantly, we include multistructural effects in the rate constant calculations. Our calculations show that the catalytic effect of water vapor is not observable at 200–400 K.

Wu, Junjun↗
DOE OSTI2021

Using projection operators with maximum overlap methods to simplify challenging self-consistent field optimization

Maximum overlap methods are effective tools for optimizing challenging ground- and excited-state wave functions using self-consistent field models such as Hartree-Fock and Kohn-Sham density functional theory. Nevertheless, such models have shown significant sensitivity to the user-defined initial guess of the target wave function. Here in this work, a projection operator framework is defined and used to provide a metric for non-aufbau orbital selection in maximum-overlap-methods. The resulting algorithms, termed the Projection-based Maximum Overlap Method (PMOM) and Projection-based Initial Maximum Overlap Method (PIMOM), are shown to perform exceptionally well when using simple user-defined target solutions based on occupied/virtual molecular orbital permutations. This work also presents a new metric that provides a simple and conceptually convenient measure of agreement between the desired target and the current or final SCF results during a calculation employing a maximum-overlap method.

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

Comparison of interband related optical transitions and excitons in ZnGeN 2 and GaN

The optical dielectric function of ZnGeN 2 is calculated from the interband transitions using the energy bands calculated in the quasiparticle self-consistent (QS)GW^ method using two different levels of approximation: the independent particle approximation (IPA) and the Bethe-Salpeter Equation (BSE) approach. The first allows us to relate peaks in ε 2 to specific bands and k-points but does not include electron-hole interaction effects. The second includes electron-hole interaction or excitonic effects. The corresponding changes in the shape of ε 2 (ω) are found to be similar to those in GaN. The screened Coulomb interaction W^ is here calculated including electron-hole interactions in the polarization function and gives a band structure already going beyond the random phase approximation. The static dielectric constants including only electronic screening, commonly referred to as ε∞ , were calculated separately by extrapolating the wave vector dependent macroscopic dielectric function in the long-wavelength limit. Below the quasiparticle gap, we find three bound excitons optically active for different polarization. The convergence of these bound excitons with respect to the density of the k-mesh used in the BSE is studied and found to require a fine mesh. It is also found that these bound excitons originate from only the lowest conduction band and the top three valence bands. To incorporate the lattice screening, we include a scaling factor (ε∞/ε 0 ) 2 , which allows us to obtain exciton binding energies of the correct order of magnitude similar to those in GaN. Here, the excitons are related to each of the three fold split valence bands and the splittings of the latter are also studied as function of strain. Finally, a relation between the anisotropic effective masses and the valence band splitting is pointed out and explained.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Acceleration of Solvation Free Energy Calculation via Thermodynamic Integration Coupled with Gaussian Process Regression and Improved Gelman–Rubin Convergence Diagnostics

The determination of the solvation free energy of ions and molecules holds profound importance across a spectrum of applications spanning chemistry, biology, energy storage, and the environment. Molecular dynamics simulations are powerful tools for computing this critical parameter. Nevertheless, the accurate and efficient calculation of the solvation free energy becomes a formidable endeavor when dealing with complex systems characterized by potent Coulombic interactions and sluggish ion dynamics and, consequently, slow transition across various metastable states. Here, in the present study, we expose limitations stemming from the conventional calculation of the statistical inefficiency g in the thermodynamic integration method, a factor that can hinder the determination of convergence of the solvation free energy and its associated uncertainty. Instead, we propose a robust scheme based on Gelman–Rubin convergence diagnostics. We leverage this improved estimation of uncertainties to introduce an innovative accelerated thermodynamic integration method based on the Gaussian Process regression. This methodology is applied to the calculation of the solvation free energy of trivalent rare-earth elements immersed in ionic liquids, a scenario in which the aforementioned challenges render standard approaches ineffective. The proposed method proves to be effective in computing solvation free energy in situations where traditional thermodynamic integration methods fall short.

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

First-principles residual resistivity using a locally self-consistent multiple scattering method

The locally self-consistent multiple scattering (LSMS) method can perform efficient first-principles calculations of systems with a large number of atoms. Here, in this paper, we combine the Kubo-Greenwood equation with LSMS, enabling us to calculate the first-principles residual resistivity of large systems. This has been implemented in the open-source code LSMS . We apply this method to selected pure elements and binary random alloys. The results compare well with experiment, and with values obtained from a first-principles effective medium technique (the Korringa-Kohn-Rostoker coherent potential approximation). We discuss future applications of this method to complex systems where other methods are not applicable.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

CI-MBPT and Intensity-Based Lifetime Calculations for Th II

Lifetime calculations of Th II J = 1.5 and 2.5 odd states are performed with configuration–interaction many-body perturbation theory (CI-MBPT). For many J = 2.5 states, lifetimes are quite accurate, but two pairs of J = 2.5 odd states and many groups of J = 1.5 states are strongly mixed, making theoretical predictions unreliable. To solve this problem, a method based on intensities is used. To relate experimental intensities to lifetimes, two parameters, one an overall coefficient of proportionality for transition rates and one temperature of the Boltzmann distribution of populations, are introduced and fitted to minimize the deviation between theoretical and intensity-derived lifetimes. For strongly mixed groups of states, the averaged lifetimes obtained from averaged transition rates were used instead of individual lifetimes in the fit. Close agreement is obtained. Then intensity branching ratios are used to extract individual lifetimes for the strongly mixed states. The resulting lifetimes are compared to available directly measured lifetimes and reasonable agreement is found, considering limited accuracy of intensity measurements. The method of intensity-based lifetime calculations with fit to theoretical lifetimes is quite general and can be applied to many complex atoms where strong mixing between multiple states exists.

74 ATOMIC AND MOLECULAR PHYSICS↗