Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “theory - 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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Theory of gradient drift instabilities in low-temperature, partially magnetised plasmas

A fluid dispersion theory in partially magnetised plasmas is analysed to examine the conditions under which large-wavelength modes develop in Penning-type configurations, that is, where an electric field is imposed perpendicular to a homogeneous magnetic field. The fluid dispersion relation assuming a slab geometry shows that two types of low-frequency, gradient drift instabilities occur in the direction of the E×B and diamagnetic drifts. One type of instability, observed when the equilibrium electric field and plasma density gradient are in the same direction, is similar to the classic modified Simon–Hoh instability. A second instability is found for conditions in which (i) the diamagnetic drift is in the direction opposite to the E×B drift and (ii) the magnitude of the diamagnetic drift is sufficiently larger than the electron thermal speed. The present fluid dispersion theory suggests that the rotating spokes driven by such fluid instabilities propagate in the same direction as the diamagnetic drift, which can be in the same direction as or opposite to the E×B drift, depending on the plasma conditions. This finding may account for the observation, in some plasma devices, of the rotation of large-scale structures in both the E×B and -E×B directions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Bounce-averaged theory in arbitrary multi-well plasmas: solution domains and the graph structure of their connections

Bounce-averaged theories provide a framework for simulating relatively slow processes, such as collisional transport and quasilinear diffusion, by averaging these processes over the fast periodic motions of a particle on a closed orbit. This procedure dramatically increases the characteristic time scale and reduces the dimensionality of the modelled system. The natural coordinates for such calculations are the constants of motion (COM) of the fast particle motion, which by definition do not change during an orbit. However, for sufficiently complicated fields – particularly in the presence of local maxima of the electric potential and magnetic field – the COM are not sufficient to specify the particle trajectory. In such cases, multiple domains in COM space must be used to solve the problem, with boundary conditions enforced between the domains to ensure continuity and particle conservation. Previously, these domains have been imposed by hand, or by recognising local maxima in the fields, limiting the flexibility of bounce-averaged simulations. Here, we present a general set of conditions for identifying consistent domains and the boundary condition connections between the domains, allowing the application of bounce-averaged theories in arbitrarily complicated and dynamically evolving electromagnetic field geometries. We also show how the connections between the domains can be represented by a directed graph, which can help to succinctly represent the trajectory bifurcation structure.

fusion plasma↗
DOE OSTI2021

Measurement and Theory of Gas-Phase Ion Mobility Shifts Resulting from Isotopomer Mass Distribution Changes

The unanticipated discovery of recent ultra-high-resolution ion mobility spectrometry (IMS) measurements revealing that isotopomers-compounds that differ only in the isotopic substitution sites-can be separated has raised questions as to the physical basis for their separation. A study comparing IMS separations for two isotopomer sets in conjunction with theory and simulations accounting for ion rotational effects provides the first-ever prediction of rotation-mediated shifts. The simulations produce observable mobility shifts due to differences in gas–ion collision frequency and translational-to-rotational energy transfer. These differences can be attributed to distinct changes in the moment of inertia and center of mass between isotopomers. The simulations are in broad agreement with the observed experiments and consistent with relative mobility differences between isotopomers. These results provide a basis for refining IMS theory and a new foundation to obtain additional structural insights through IMS.

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

Applicability of Allen–Heine–Cardona Theory on MO x Metal Oxides and ABO 3 Perovskites: Toward High-Temperature Optoelectronic Applications [plus supplemental information]

We report characterizing the temperature-dependent electronic structures of optoelectronic materials is vital for their application in high-temperature processes. Allen–Heine–Cardona (AHC) theory is an established theory to determine the impact of electron–phonon coupling on the temperature renormalization of the electronic structure. A question of fundamental interest is how applicable this approach can be in a broad range of materials. Herein, we establish the consensus on the temperature dependence of band gaps between AHC predictions and experimental measurements, in a versatile set of nine MO x metal oxides and six ABO 3 perovskites that are utilized as gas sensors. Both calculations and measurements point to the appreciable band gap downshifts within the temperature range relevant to typical high-temperature applications in all examined materials. The correlation between the band gaps and temperatures was quantified at zero- and finite-temperature limits using the analytical parameters of the O’Donnell model. We found fairly comparative agreements between AHC- and experiment-derived parameters, especially for the parameters (S and $\langle\hbar$ω$\rangle$) defining the temperature dependence behavior of band gaps at the finite-temperature regime. Our findings provide important implications for the first-principles modeling of the electronic structures and their temperature dependency, which should be the focus toward high-temperature optoelectronic applications.

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

Theory-Guided Inelastic Neutron Scattering of Crystalline Alkaline Aluminate Salts Bearing Principal Motifs of Solution-State Species

Aluminate salts precipitated from caustic alkaline solutions exhibit a correlation between the anionic speciation and the identity of the alkali cation in the precipitate, with the aluminate ions occurring either in monomeric (Al(OH) 4 – ) or dimeric (Al 2 O(OH) 6 2– ) forms. The origin of this correlation is poorly understood as are the roles that oligomeric aluminate species play in determining the solution structure, prenucleation clusters, and precipitation pathways. Characterization of aluminate solution speciation with vibrational spectroscopy results in spectra that are difficult to interpret because the ions access a diverse and dynamic configurational space. To investigate the Al(OH) 4 – and Al 2 O(OH) 6 2– anions within a well-defined crystal lattice, inelastic neutron scattering (INS) and Raman spectroscopic data were collected and simulated by density functional theory for K 2 [Al 2 O(OH) 6 ], Rb 2 [Al 2 O(OH) 6 ], and Cs[Al(OH) 4 ]·2H 2 O. These structures capture archetypal solution aluminate species: the first two salts contain dimeric Al 2 O(OH) 6 2– anions, while the third contains the monomeric Al(OH) 4 – anion. Here, comparisons were made to the INS and Raman spectra of sodium aluminate solutions frozen in a glassy state. In contrast to solution systems, the crystal lattice of the salts results in well-defined vibrations and associated resolved bands in the INS spectra. The use of a theory-guided analysis of the INS of this solid alkaline aluminate series revealed that differences were related to the nature of the hydrogen-bonding network and showed that INS is a sensitive probe of the degree of completeness and strength of the bond network in hydrogen-bonded materials. Results suggest that the ionic size may explain cation-specific differences in crystallization pathways in alkaline aluminate salts.

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

Predicting Mössbauer Parameters of Nonheme Diiron Complexes with Density Functional Theory

Mössbauer spectroscopy provides significant insights into the electronic structure and environment of the metal centers. In this work, we investigate the electronic structures of a set of nonheme diiron complexes by evaluating two key parameters pertaining to Mössbauer spectroscopy, namely, the isomer shift (δ) and quadrupole splitting (|ΔE Q |), using different levels of density functional theory (DFT). The diiron systems investigated here span diverse oxidation states, bridging motifs, and spin coupling patterns, which present a challenging case for theoretical predictions. We demonstrate that the combination of B97-D3/def2-TZVP is an efficient approach in modeling both the δ and |ΔEQ| values with high accuracy for the representative nonheme diiron complexes. We also show that δ is accurately predicted irrespective of the choice of approximate density functional while the |ΔE Q | is sensitive to the level of theory employed. Further investigation shows that the present methodology assessed using synthetic nonheme diiron complexes could be extended to nonheme diiron enzyme active sites, featuring both ferromagnetic and antiferromagnetic coupling between the iron centers.

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

Unmasking Static Correlation Error in Hybrid Kohn–Sham Density Functional Theory

To uncover the way Hartree–Fock exchange brings static correlation error into hybrid Kohn–Sham density functional theory, we compare the potential energy curves of four diatomic molecules, namely, H 2 , F 2 , HF, and NaF, using both restricted and unrestricted Kohn–Sham theory. Here, we find that increasing the percentage of Hartree–Fock exchange significantly increases the static correlation error because Hartree–Fock exchange replaces the ability of local exchange in density functionals to account for the localization effects associated with static correlation.

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

Predicting Excitation Energies of Twisted Intramolecular Charge-Transfer States with the Time-Dependent Density Functional Theory: Comparison with Experimental Measurements in the Gas Phase and Solvents Ranging from Hexanes to Acetonitrile

Electronically excited states characterized by intramolecular charge transfer play an essential role in many biological processes and optical devices. The ability to make quantitative ab initio predictions of the relative energetics involved is a challenging yet desirable goal, especially for large molecules in solution. In this work, we present a data set of 61 experimental measurements of absorption and emission processes, both in the gas phase and in solvents representing a broad range of polarities, which involve intramolecular charge transfer mediated by a nonzero, "twisted" dihedral angle between one or more donor and acceptor subunits. Among a variety of density functionals investigated within the framework of linear-response theory, the "optimally tuned" LRC-ωPBE functional, which utilizes a system-specific yet nonempirical procedure to specify the range-separation parameter, emerges as the preferred choice. For the entire set of excitation energies, involving changes in dipole moment ranging from 4 to >20 Debye, the mean signed and absolute errors are 0.02 and 0.18 eV, respectively (compared, e.g., to -0.30 and 0.30 for PBE0, 0.44 and 0.47 for LRC-ωPBEh, 0.83 and 0.83 for ωB97X-V). We analyze the performance of polarizable continuum solvation models available in Q-Chem that partition the solvent response into fast and slow time scales, and clear trends emerge when measurements corresponding to the four small 4-(dimethylamino)benzonitrile (DMABN)-like molecules and a charged species are excluded. We make the case that the large errors found only for small molecules in the gas phase and weak solvents cannot be expected to improve via the optimal tuning procedure, which enforces a condition that is exact only in the well-separated donor-acceptor limit, and present empirical evidence implicating the outsized importance for small donor-acceptor systems of relaxation effects that cannot be accounted for by the linear-response time-dependent density functional theory within the adiabatic approximation. Finally, we demonstrate the utility of the optimally tuned density functional approach by targeting the charge-transfer states of a large biomimetic model system for light-harvesting structures in Photosystem II.

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

Green’s Function Formulation of Quantum Defect Embedding Theory

Here, we present a Green’s function formulation of the quantum defect embedding theory (QDET) where a double counting scheme is rigorously derived within the $G_0W_0$ approximation. We then show the robustness of our methodology by applying the theory with the newly derived scheme to several defects in diamond. Additionally, we discuss a strategy to obtain converged results as a function of the size and composition of the active space. Our results show that QDET is a promising approach to investigate strongly correlated states of defects in solids.

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

Tensor Hypercontraction of Cluster Perturbation Theory: Quartic Scaling Perturbation Series for the Coupled Cluster Singles and Doubles Ground-State Energies

Even though Cluster Perturbation has been shown to be a robust non-iterative alternative to Coupled Cluster Theory, it is still plagued by high order polynomial computational scaling and the storage of higher order tensors. Here we present a proof-of-concept strategy for implementing Cluster Perturbation Theory ground state energy series for the coupled cluster singles and doubles energy with N 4 computational scaling using Tensor Hypercontraction (THC). The reduction in computational scaling by two orders is achieved by decomposing two electron repulsion integrals, doubles amplitudes and multipliers, as well as selected doubles intermediates to Tensor Hypercontraction format. Using the outlined strategy, we showcase that the Tensor Hypercontraction pilot implementations retain numerical accuracy to within 1 kcal/mol relative to corresponding conventional and density fitting implementations and we empirically verify the N 4 scaling.

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

Excited-State Densities from Time-Dependent Density Functional Response Theory

While the variational principle for excited-state energies leads to a route to obtaining excited-state densities from time-dependent density functional theory, relatively little attention has been paid to the quality of the resulting densities in real space obtained with different exchange-correlation functional approximations or how nonadiabatic approximations developed for energies of states of double-excitation character perform for their densities. Here we derive an expression directly in real space for the excited-state density, which includes the case of nonadiabatic kernels and consequently is able, for the first time, to yield densities of states of double-excitation character. Under some well-defined simplifications, we compare the performance of the local-density approximation and exact-exchange approximation, which are in a sense at the opposite extremes of the fundamental functional approximations, on local and charge-transfer excitations in one-dimensional model systems and show that the dressed Time-Dependent Density Functional Theory (TDDFT) approach gives good densities of double excitations.

approximation↗
DOE OSTI2020

Spectral Narrowing Accompanies Enhanced Spatial Resolution in Saturated Coherent Anti-Stokes Raman Scattering (CARS): Comparisons of Experiment and Theory

Here, we demonstrate theoretically and confirm experimentally the mechanism by which spectral narrowing accompanies enhanced spatial resolution in a saturated coherent anti-Stokes Raman scattering (CARS) signal that is demodulated at the third harmonic (3$f$) of the pump modulation frequency ($f$). Under these modulation conditions, theory predicts a narrowing of the full-width at half-maximum (FWHM) of the CARS spectrum by a factor of 2.0 with respect to that of the spectrum obtained by demodulation at the fundamental frequency. Theory also predicts an improvement of spatial resolution by a factor of 1.7. Experimentally, narrowing of the FWHM of the CARS spectrum of 1,4-bis((E)-2-methylstyryl) benzene (MSB) crystals by a factor of 2.5 is observed upon saturation. Further experimental confirmation is provided from investigating diamond particles, for which spectral narrowing was enhanced by a factor of 2.8 and spatial resolution was enhanced by a factor of 2. Details of the mechanism and execution of the saturated CARS experiment are elucidated and limits to its applicability are suggested, one of which is the conclusion that the saturation approach is not suitable for extraction of harmonics beyond 3$f$ . In this work we have developed a more comprehensive understanding of the correlation between the observed experimental results and experimental factors than has been previously reported.

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

Reactions of U + with H 2 , D 2 , and HD Studied by Guided Ion Beam Tandem Mass Spectrometry and Theory

The kinetic energy-dependent reactions of the atomic actinide uranium cation (U + ) with H 2 , D 2 , and HD were examined by guided ion beam tandem mass spectrometry. An average 0 K bond dissociation energy of D 0 (U + – H) = 2.48 ± 0.06 eV is obtained by analysis of the endothermic product ion cross sections. Quantum chemistry calculations were performed for comparison with experimental thermochemistry, including high-level CASSCF–CASPT2–RASSI calculations of the spin–orbit corrections. CCSD(T) and the CASSCF levels show excellent agreement with experiment, whereas B3LYP and PBE0 slightly overestimate and the M06 approach badly underestimates the bond energy for UH + . Theory was also used to investigate the electronic structures of the reaction intermediates and potential energy surfaces. The experimental product branching ratio for the reaction of U + with HD indicates that these reactions occur primarily via a direct reaction mechanism, despite the presence of a deep-well for UH 2 + formation according to theory. The reactivity and hydride bond energy for U + are compared with those for transition metal, lanthanide, and actinide cations, and periodic trends are discussed. Furthermore, these comparisons suggest that the 5f electrons on uranium are largely core and uninvolved in the reactive chemistry.

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

Resolving Discrepancies between State-of-the-Art Theory and Experiment for HO 2 + HO 2 via Multiscale Informatics

Recent high-level theoretical calculations predict a mild temperature dependence for HO 2 + HO 2 inconsistent with state-of-the-art experimental determinations that upheld the stronger temperature dependence observed in early experiments. Via MultiScale Informatics analysis of the theoretical and experimental data, we identified an alternative interpretation of the raw experimental data that uses HO 2 + HO 2 rate constants nearly identical to theoretical predictions---implying that the theoretical and experimental data are actually consistent, at least when considering the raw data from experimental studies. Here, similar analyses of typical signals from low-temperature experiments indicate that an HOOOOH intermediate---identified by recent theory but absent from earlier interpretations---yields modest effects that are smaller than, but may have contributed to, the scatter in data among different experiments. More generally, the findings demonstrate that modern chemical theories and experiments have progressed to a point where meaningful comparison requires joint consideration of their data simultaneously.

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

Hydration Free Energies of Polypeptides from Popular Implicit Solvent Models versus All-Atom Simulation Results Based on Molecular Quasichemical Theory

Calculating the hydration free energy of a macromolecule in all-atom simulations has long remained a challenge, necessitating the use of models wherein the effect of the solvent is captured without explicit account of solvent degrees of freedom. This situation has changed with developments in the molecular quasi-chemical theory (QCT)-an approach that enables calculation of the hydration free energy of macromolecules within all-atom simulations at the same resolution as is possible for small molecular solutes. The theory also provides a rigorous and physically transparent framework to conceptualize and model interactions in molecular solutions and thus provides a convenient framework to investigate the assumptions in implicit solvent models. In this study, we compare the results using molecular QCT versus predictions from EEF1, ABSINTH, and GB/SA implicit solvent models for polyglycine and polyalanine solutes covering a range of chain lengths and conformations. The hydration free energies or the differences in hydration free energies between conformers obtained from the implicit solvent models do not agree with explicit solvent results, with the deviations being largest for the group additive EEF1 and ABSINTH models. GB/SA does better in capturing the qualitative trends seen in explicit solvent results. Finally, analysis founded on QCT reveals the critical importance of the cooperativity of hydration that is inherent in the hydrophilic and hydrophobic contributions to hydration-physics that is not well captured in additive models but somewhat better accounted for by means of a dielectric in the GB/SA approach.

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

Hole Polaron Transport in Bismuth Vanadate BiVO 4 from Hybrid Density Functional Theory

In this paper, we address the chemical character and the mobility of hole polarons in BiVO 4 (BVO). Two distinct structures for hole polarons in BVO have been reported in the literature to date. “Standard” hybrid density functional theory (DFT) calculations have predicted hole polarons to be (BiO 8 ) dodecahedron-centered polarons, while Hubbard U-DFT (DFT + U) calculations have predicted single oxygen atom-centered polarons. Resolving these contrasting findings is critical to validate theory in materials research. We investigated the structure of hole polarons in BVO using hybrid DFT with varying fractions (α) of the exact exchange interaction that enters hybrid functionals. For values of the fraction α, from α = 0.25 to α = 0.45, we obtained both structures, the h + (BiO 8 ) hole structure and the h + (O) hole structure. For the smaller values of α, the h + (BiO 8 ) structure is lower in energy, while for the larger values of α, the h+(O) structure is more stable. The piece-wise linearity of DFT energies E(q) with respect to a point defect charge q is maintained satisfactorily in all cases, although the deviation is the smallest for hybrid DFT with a small exact exchange fraction. For α = 0.25, we determined the diffusivity of h + (BiO 8 ) holes from the Marcus/Holstein two-state model. The activation barrier for Bi-to-Bi hops (ΔG # ~76 meV) and mobility (μ ~0.029 cm2 V –1 s –1 ) are found to be in close agreement with the published THz experimental data (ΔG # ~90 meV and μ exp ~0.02 cm2 V –1 s –1 , respectively). In h + (BiO 8 ) structures, nearly 25% of the hole charge and unpaired spin reside on the Bi atom and the rest is divided among the eight O atoms. The hole is best described as a Bi-centered hole, rather than a Bi 6s hole. In contrast, in h + (O) structures, about 75% of the unpaired spin population resides on the O atom. These findings suggest that if one wishes to resolve this theoretical dichotomy, there remains a need for further experimental characterization of hole polarons in BVO with techniques that could identify the oxidation state of Bi and O atoms in operando conditions and help validate theoretical results of structure and transport.

14 SOLAR ENERGY↗
DOE OSTI2022

Characterization of Lanthanum Monazite Surface Chemistry and Crystal Morphology through Density Functional Theory and Experimental Approaches

Monazite is a rare earth element (REE)-containing mineral that consists of (REE)PO 4 formal units and is one of the most important sources of these critical materials. The concentration of REEs from mined monazite ore often involves froth flotation, which is a beneficiation process that enhances the efficiency of downstream processing. The effectiveness of froth flotation is largely governed by the ability of collector agents to selectively bind to monazite particles. Thus, a molecular-level understanding of monazite interfacial chemistry is integral to the design of effective collector agents. To address this need, we performed density functional theory (DFT) calculations and a variety of experimental techniques to characterize La-monazite and elucidate its crystal morphology. Interestingly, we find minimal differences in the predicted morphologies of La-monazite for hydrous and anhydrous environments, which are largely dominated by low-index facets (e.g., {110}, {100}, and {010}). Indexing of synthesized La-monazite crystals via X-ray diffraction also uncovers {110} and {100} as the predominant facets. The average surface energies of 0% and 100% water coverage La-monazite crystals were predicted to be 0.87 and 0.76 J/m 2 , respectively, while calorimetry suggests values of 1.30 and 1.15 J/m 2 , respectively. The apparent discrepancies between the theoretical and experimental values are expected and attributed to defects present in physical crystals, in contrast to the perfect mineral surfaces in simulations. The difference in surface energy between the 0% and 100% water coverage morphologies predicted by theory is consistent with the value measured via calorimetry. DFT reveals a wide range of adsorption energies for water across the studied facets, but in all cases, water is predicted to strongly bind to monazite surfaces with an average adsorption energy of -92.7 kJ/mol for a La-monazite single crystal. Finally, this study provides the groundwork necessary for the rational design of froth flotation collector agents by granting molecular-level insight into the predominant facets of monazite.

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

Excitons in Hematite Fe 2 O 3 : Short-Time Dynamics from TD-DFT and Non-Adiabatic Dynamics Theories

We present a first-principles study of the short-time dynamics of excitons in hematite Fe 2 O 3 . We used time-dependent density functional theory (TD-DFT) with an underlying DFT+U treatment of electron interactions to characterize the electronic structure of excitons and nonadiabatic molecular dynamics theory (NA-MD) to determine their recombination (electronic ground-state recovery) and relaxation dynamics. Decoherence-corrected trajectory surface hopping approaches in NA-MD simulations yielded recovery times of ~1.1 to 1.8 ns and “higher-lying” exciton relaxation times of ~60 to 70 fs, in accord with experimentally derived lifetimes. With hematite phonons in the range of ~100 to 700 cm –1 , higher-lying excitons relax within one or two oscillations of the phonons before getting trapped into an electron–hole pair Exc-3 structure on the first excited state potential energy surface. This structure resembles already a pair of polarons (electron polaron plus hole polaron) with associated lattice distortions three (3) basal planes away. On longer time scales, the electron–hole bipolaronic pair hops to structures Exc-5, then Exc-7, then Exc-9, ... with the electron polaron and hole polaron separated by 5, 7, 9, ... basal planes in a process of charge separation. The largest frequency phonon ~672 cm –1 for the Exc-3 exciton structure is associated with the electron polaron moiety of the exciton. This phonon is a good candidate for giving rise to the recently observed and reported postexcitation transient IR absorption peak.

36 MATERIALS SCIENCE↗