Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “self consistent field 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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Bootstrap embedding for interacting electrons in phonon coherent-state mean field

Here, we develop a Fermi–Bose bootstrap embedding framework for the ground state of interacting electrons coupled to a phonon mean field. The method combines bootstrap embedding for correlated electrons with a self-consistent coherent-state mean-field treatment for phonons. This method models the interacting electron–phonon problem as a system of correlated electrons traveling in a self-consistently specified potential landscape, allowing for efficient treatment of large lattice systems. Convergence of the methods for fragment size and total system size is demonstrated for the one-dimensional Hubbard–Holstein model for up to 350 sites. Finite-size scaling is performed to extrapolate to the infinite system size. Benchmarking against the density matrix renormalization group for a small 8-site system at half- and quarter-filling shows an orders-of-magnitude runtime advantage. The comparison further reveals that the method performs best in regimes dominated by localization, such as the Mott insulating phase and the strong-coupling tiny polaron regime, where the local embedding ansatz is still valid. However, due to the mean-field treatment for phonons, we find limitations of our methods in the weakly coupled delocalized region and at the Peierls transition, where quantum phonon fluctuations and long-range kinetic correlations become substantial.

Islam, Shariful [North Carolina State University, ↗
DOE OSTI2024

Exact-Two-Component Complete Active Space Method with Variational Treatment of Magnetic Field and Spin–Orbit Coupling: Application to X-ray Magnetic Circular Dichroism Spectroscopy

We introduce an exact-two-component complete active space self-consistent-field (X2C-CASSCF) method formulated under the restricted-magnetic-balance condition. This framework allows for the nonperturbative treatment of static magnetic fields using gauge-including atomic orbitals (GIAOs). The GIAO-X2C-CASSCF methodology effectively captures all microstates within the same 2J + 1-degenerate manifold and their splitting in a static magnetic field, which are not accessible through single-reference-based methods. We also present mathematical recursive expressions for evaluating one-electron relativistic integrals by using GIAOs in the presence of a finite magnetic field. Benchmark studies include oxygen and nitrogen K-edge X-ray magnetic circular dichroism spectroscopy (XMCD) for closed-shell organic compounds, as well as L-edge XMCD spectroscopy for the high-spin open-shell transition metal ion Mn 2+ and the tetrahedral Mn(II)O 4 6– complex.

Chemical calculations↗
NASA NTRS1994

The Fraternal Twins of Quartet O4(+)

Eleven stationary geometries of quartet O4(+) have been studied by ab initio methods. The geometries were optimized at the Complete Active Space Self-Consistent Field (CASSCF) level of theory and the energies were calculated by the multiconfigurational second order pertubation method (CASPT2), using Double-Zeta Plus polarization (DZP), Triple-Zeta Plus Double Polarization (TZ2P), average Atomic Natural Orbital (ANO) (5s4p2d) and average ANO (6s5p3d2f) basis sets. The rectangular and trans-planar structures are found to be the most stable, with an energy barrier to conversion between the two at the threshold of dissociation. Both have a delocalized hole and are stable relative to separated 02 and 02(+)by 11.0 and 11.5 kcal/mol for the rectangular and the 2 trans-planar structure, respectively, compared with the experimentally deduced energy in the range of 9.2 to 10.8 kcal/mol. The adiabatic ionization potentials of 04 and 02 are computed to be 11.67 and 12.21 eV, while experimental values are 11.66 and 12.07 eV, respectively. The vibrational frequencies have been computed for all degrees of freedom at the CASSCF level of theory. Symmetry breaking is found to be a particular problem in the computation of the antisymmetric stretch frequency for the delocalized structures at the CASSCF level of theory. Attempts to rectify these problems using the restricted active space self-consistent field (RASSCF) method leads to additional difficulties, but further analysis yields insight into the symmetry breaking and problems with earlier calculations. Finally, a nonorthogonal configuration interaction (CI) calculation based on the interaction of localized CASSCF wave functions using the Complete Active Space State Interation (CASSI) method leads to a balanced treatment of the antisymmetric stretch which is free from symmetry breaking. The study explains the four most prominent absorption frequencies observed in the partially unassigned IR spectrum of O4(+) isolated in solid neon as the antisymmetric OO stretch, and the combination band of the symmetric and antisymmetric 00 stretch of both the rectangular and trans-planar structures.

Lindh, Roland↗
NASA NTRS1971

The closure approximation in the hierarchy equations.

The expectation of the solution process in a stochastic operator equation can be obtained from averaged equations only under very special circumstances. Conditions for validity are given and the significance and validity of the approximation in widely used hierarchy methods and the ?self-consistent field' approximation in nonequilibrium statistical mechanics are clarified. The error at any level of the hierarchy can be given and can be avoided by the use of the iterative method.

Adomian, G.↗
DOE OSTI2021

Site-Selective Real-Time Observation of Bimolecular Electron Transfer in a Photocatalytic System Using L-Edge X-Ray Absorption Spectroscopy

Here, time-resolved X-ray absorption spectroscopy has been utilized to monitor the bimolecular electron transfer in a photocatalytic water splitting system. This has been possible by uniting the local probe and element specific character of X-ray transitions with insights from high-level ab initio calculations. The specific target has been a heteroleptic [Ir III (ppy) 2 (bpy)] + photosensitizer, in combination with triethylamine as a sacrificial reductant and Fe 3 (CO) 12 as a water reduction catalyst. The relevant molecular transitions have been characterized via high-resolution Ir L-edge X-ray absorption spectroscopy on the picosecond time scale and restricted active space self-consistent field calculations. The presented methods and results will enhance our understanding of functionally relevant bimolecular electron transfer reactions and thus will pave the road to rational optimization of photocatalytic performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2001

Ab Initio and Improved Empirical Potentials for the Calculation of the Anharmonic Vibrational States and Intramolecular Mode Coupling of N-Methylacetamide

The second-order Moller-Plesset ab initio electronic structure method is used to compute points for the anharmonic mode-coupled potential energy surface of N-methylacetamide (NMA) in the trans(sub ct) configuration, including all degrees of freedom. The vibrational states and the spectroscopy are directly computed from this potential surface using the Correlation Corrected Vibrational Self-Consistent Field (CC-VSCF) method. The results are compared with CC-VSCF calculations using both the standard and improved empirical Amber-like force fields and available low temperature experimental matrix data. Analysis of our calculated spectroscopic results show that: (1) The excellent agreement between the ab initio CC-VSCF calculated frequencies and the experimental data suggest that the computed anharmonic potentials for N-methylacetamide are of a very high quality; (2) For most transitions, the vibrational frequencies obtained from the ab initio CC-VSCF method are superior to those obtained using the empirical CC-VSCF methods, when compared with experimental data. However, the improved empirical force field yields better agreement with the experimental frequencies as compared with a standard AMBER-type force field; (3) The empirical force field in particular overestimates anharmonic couplings for the amide-2 mode, the methyl asymmetric bending modes, the out-of-plane methyl bending modes, and the methyl distortions; (4) Disagreement between the ab initio and empirical anharmonic couplings is greater than the disagreement between the frequencies, and thus the anharmonic part of the empirical potential seems to be less accurate than the harmonic contribution;and (5) Both the empirical and ab initio CC-VSCF calculations predict a negligible anharmonic coupling between the amide-1 and other internal modes. The implication of this is that the intramolecular energy flow between the amide-1 and the other internal modes may be smaller than anticipated. These results may have important implications for the anharmonic force fields of peptides, for which N-methylacetamide is a model.

Gregurick, Susan K.↗
NASA NTRS2002

Gas-Phase Reaction Pathways and Rate Coefficients for the Dichlorosilane-Hydrogen and Trichlorosilane-Hydrogen Systems

As part of NASA Ames Research Center's Integrated Process Team on Device/Process Modeling and Nanotechnology our goal is to create/contribute to a gas-phase chemical database for use in modeling microelectronics devices. In particular, we use ab initio methods to determine chemical reaction pathways and to evaluate reaction rate coefficients. Our initial studies concern reactions involved in the dichlorosilane-hydrogen (SiCl2H2--H2) and trichlorosilane-hydrogen (SiCl2H-H2) systems. Reactant, saddle point (transition state), and product geometries and their vibrational harmonic frequencies are determined using the complete-active-space self-consistent-field (CASSCF) electronic structure method with the correlation consistent polarized valence double-zeta basis set (cc-pVDZ). Reaction pathways are constructed by following the imaginary frequency mode of the saddle point to both the reactant and product. Accurate energetics are determined using the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations (CCSD(T)) extrapolated to the complete basis set limit. Using the data from the electronic structure calculations, reaction rate coefficients are obtained using conventional and variational transition state and RRKM theories.

Dateo, Christopher E.↗
DOE OSTI2020

Fast and Universal Kohn-Sham Density Functional Theory Algorithm for Warm Dense Matter to Hot Dense Plasma

Understanding many processes, e.g., fusion experiments, planetary interiors, and dwarf stars, depends strongly on microscopic physics modeling of warm dense matter and hot dense plasma. This complex state of matter consists of a transient mixture of degenerate and nearly free electrons, molecules, and ions. This regime challenges both experiment and analytical modeling, necessitating predictive ab initio atomistic computation, typically based on quantum mechanical Kohn-Sham density functional theory (KS-DFT). However, cubic computational scaling with temperature and system size prohibits the use of DFT through much of the warm dense matter regime. A recently developed stochastic approach to KS-DFT can be used at high temperatures, with the exact same accuracy as the deterministic approach, but the stochastic error can converge slowly and it remains expensive for intermediate temperatures (< 50 eV). Here we have developed a universal mixed stochastic-deterministic algorithm for DFT at any temperature. This approach leverages the physics of KS-DFT to seamlessly integrate the best aspects of these different approaches. We demonstrate that this method significantly accelerated self-consistent field calculations for temperatures from 3 to 50 eV, while producing stable molecular dynamics and accurate diffusion coefficients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

Accuracy of the Generalized Self-Consistent Method in Modelling the Elastic Behaviour of Periodic Composites

Local stress and strain fields in the unit cell of an infinite, two-dimensional, periodic fibrous lattice have been determined by an integral equation approach. The effect of the fibres is assimilated to an infinite two-dimensional array of fictitious body forces in the matrix constituent phase of the unit cell. By subtracting a volume averaged strain polarization term from the integral equation we effectively embed a finite number of unit cells in a homogenized medium in which the overall stress and strain correspond to the volume averaged stress and strain of the constrained unit cell. This paper demonstrates that the zeroth term in the governing integral equation expansion, which embeds one unit cell in the homogenized medium, corresponds to the generalized self-consistent approximation. By comparing the zeroth term approximation with higher order approximations to the integral equation summation, both the accuracy of the generalized self-consistent composite model and the rate of convergence of the integral summation can be assessed. Two example composites are studied. For a tungsten/copper elastic fibrous composite the generalized self-consistent model is shown to provide accurate, effective, elastic moduli and local field representations. The local elastic transverse stress field within the representative volume element of the generalized self-consistent method is shown to be in error by much larger amounts for a composite with periodically distributed voids, but homogenization leads to a cancelling of errors, and the effective transverse Young's modulus of the voided composite is shown to be in error by only 23% at a void volume fraction of 75%.

Walker, Kevin P.↗
DOE OSTI2020

A Shift Selection Strategy for Parallel Shift-invert Spectrum Slicing in Symmetric Self-consistent Eigenvalue Computation

The central importance of large-scale eigenvalue problems in scientific computation necessitates the development of massively parallel algorithms for their solution. Recent advances in dense numerical linear algebra have enabled the routine treatment of eigenvalue problems with dimensions on the order of hundreds of thousands on the world’s largest supercomputers. In cases where dense treatments are not feasible, Krylov subspace methods offer an attractive alternative due to the fact that they do not require storage of the problem matrices. However, demonstration of scalability of either of these classes of eigenvalue algorithms on computing architectures capable of expressing massive parallelism is non-trivial due to communication requirements and serial bottlenecks, respectively. In this work, we introduce the SISLICE method: a parallel shift-invert algorithm for the solution of the symmetric self-consistent field (SCF) eigenvalue problem. The SISLICE method drastically reduces the communication requirement of current parallel shift-invert eigenvalue algorithms through various shift selection and migration techniques based on density of states estimation and k-means clustering, respectively. This work demonstrates the robustness and parallel performance of the SISLICE method on a representative set of SCF eigenvalue problems and outlines research directions that will be explored in future work.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Diffusion behavior of lanthanide fission products in bcc Fe cladding: A first-principles study

Fuel-cladding chemical interaction poses significant challenges in nuclear reactors, where fission products generated from nuclear fuel interact with Fe-based cladding materials, potentially compromising their structural integrity. This study investigates the diffusion behavior of lanthanide fission products, Lanthanum (La), Cerium (Ce), Praseodymium (Pr), and Neodymium (Nd), within body-centered cubic (bcc) Fe cladding using the density functional theory, nudged elastic band method, and self-consistent mean field theory. Our results reveal significant vacancy binding energies, particularly with the 1st and 2nd nearest neighbors, which diminish beyond the 5th nearest neighbor, with La exhibiting the strongest binding affinity, followed by Nd, Ce, and Pr. The nudged elastic band calculations indicate significant high barriers for the dissociation of 1st nearest neighbor vacancy-solute pairs for all fission products. The tracer diffusion coefficients of these fission products were derived in an Arrhenius form, with a magnetic correction that accounts for the high-temperature paramagnetic state. The significant trapping effect of vacancies caused by a very dilute concentration of fission products reduces vacancy mobility, potentially leading to modifications in point defect supersaturation, void nucleation, and swelling under irradiation. These represent critical challenges for irradiated cladding materials. The tracer diffusion coefficients indicate that Nd diffuses the fastest, followed by La, Ce, and Pr. Furthermore, this study provides essential insights for understanding fission product transport in cladding materials and informs future design strategies to mitigate fuel-cladding chemical interaction, ultimately enhancing nuclear reactor safety and performance.

Diffusion↗
DOE OSTI2022

Dynamical electron correlation and the chemical bond. II. Recoupled pair bonds in the $a$ 4 Σ - states of CH and CF

We extended our studies of the effect of dynamical electron correlation on the covalent bonds in the AH and AF series (A = B–F) to the recoupled pair bonds in the excited $a$ 4 Σ - states of the CH and CF molecules. Dynamical correlation is energetically less important in the $a$ 4 Σ - states than in the corresponding X 2 Π states for both molecules, which is reflected in smaller changes in bond energies ($D$ e ). Changes in the equilibrium bond distance ($R$ e ) and vibrational frequency ($ω$ e ), on the other hand, are influenced by the changes in the slope and curvature of the dynamical electron correlation energy as a function of the internuclear distance $R$, $E$ DEC ($R$). In the CH($a$ 4 Σ - ) state, these changes are much smaller than in the CH(X 2 Π) state, but in the CF($a$ 4 Σ - ) state, they are larger, reflecting a significant difference in the shapes of $E$ DEC ($R$) curves. At large $R$, the shape of $E$ DEC ($R$) curves for covalent and recoupled pair bonds is similar although different in magnitude. For the CH($a$ 4 Σ - ) state, $E$ DEC ($R$) has a minimum at $R$ = $R$ e + 0.72 Å as the orbitals associated with the formation of the recoupled pair bond switch places. $E$ DEC (R) for the CF($a$ 4 Σ - ) state decreases continuously throughout the bound region of the potential energy curve because the dynamical correlation energy associated with the electrons in the lone pair orbitals is increasing. In conclusion, these results support our earlier conclusion that we still have much to learn about the nature of dynamical electron correlation in molecules.

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

Diffusion Behavior of Oversized Fission Products in bcc Fe Cladding: A First-Principles Study

Fuel-Cladding Chemical Interaction (FCCI) poses significant challenges in nuclear reactors, where fission products from nuclear fuel interact with Fe-based cladding materials, potentially compromising their structural integrity. This study investigates the diffusion behavior of oversized fission products, Pr, Nd, Ce, and La, within bcc Fe cladding using density functional theory (DFT), nudged elastic band (NEB) method, and self-consistent mean field (SCMF) theory. Our results reveal significant long-range vacancy binding energies, particularly up to the 6th nearest neighbor, with La exhibiting the strongest binding affinity, followed by Nd, Ce, and Pr. The NEB calculations indicate significant high barriers for the dissociation of 1nn vacancy-solute pairs for all fission products. The tracer diffusion coefficients of these fission products was derived in Arrhenius form. The significant trapping effect of vacancies by a very dilute amount of fission products reduces vacancy mobility, leading to an oversaturation of point defects, void nucleation, and swelling. These are critical issues for irradiated cladding materials. The tracer diffusion coefficients indicate that Nd diffuses the fastest, followed by La, Ce, and Pr. This study provides essential insights for developing advanced cladding materials and design strategies to mitigate FCCI, ultimately enhancing nuclear reactor safety and performance.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1991

An overview of self-consistent methods for fiber-reinforced composites

The Walker et al. (1989) self-consistent method to predict both the elastic and the inelastic effective material properties of composites is examined and compared with the results of other self-consistent and elastically based solutions. The elastic part of their method is shown to be identical to other self-consistent methods for non-dilute reinforced composite materials; they are the Hill (1965), Budiansky (1965), and Nemat-Nasser et al. (1982) derivations. A simplified form of the non-dilute self-consistent method is also derived. The predicted, elastic, effective material properties for fiber reinforced material using the Walker method was found to deviate from the elasticity solution for the v sub 31, K sub 12, and mu sub 31 material properties (fiber is in the 3 direction) especially at the larger volume fractions. Also, the prediction for the transverse shear modulus, mu sub 12, exceeds one of the accepted Hashin bounds. Only the longitudinal elastic modulus E sub 33 agrees with the elasticity solution. The differences between the Walker and the elasticity solutions are primarily due to the assumption used in the derivation of the self-consistent method, i.e., the strain fields in the inclusions and the matrix are assumed to remain constant, which is not a correct assumption for a high concentration of inclusions.

Gramoll, Kurt C.↗
NASA NTRS2001

Ab Initio Calculations of Anharmonic Vibrational Spectroscopy for Hydrogen Fluoride (HF)n (n=3,4) and Mixed Hydrogen Fluoride/Water (HF)n(H20)n (n=1,2,4) Clusters

Anharmonic vibrational frequencies and intensities are computed for hydrogen fluoride clusters (HF)n with n=3,4 and mixed clusters of hydrogen fluoride with water (HF)n(H2O)n where n=1,2. For the (HF)4(H2O)4 complex, the vibrational spectra are calculated at the harmonic level, and anharmonic effects are estimated. Potential energy surfaces for these systems are obtained at the MP2/TZP level of electronic structure theory. Vibrational states are calculated from the potential surface points using the correlation-corrected vibrational self-consistent field (CC-VSCF) method. The method accounts for the anharmonicities and couplings between all vibrational modes and provides fairly accurate anharmonic vibrational spectra that can be directly compared with experimental results without a need for empirical scaling. For (HF)n, good agreement is found with experimental data. This agreement shows that the MP2 potential surfaces for these systems are reasonably reliable. The accuracy is best for the stiff intramolecular modes, which indicates the validity of MP2 in describing coupling between intramolecular and intermolecular degrees of freedom. For (HF)n(H2O)n experimental results are unavailable. The computed intramolecular frequencies show a strong dependence on cluster size. Intensity features are predicted for future experiments.

Chaban, Galina M.↗
NASA NTRS2003

Anharmonic Vibrational Spectroscopy of the F-(H20)n, complexes, n=1,2

We report anharmonic vibrational spectra (fundamentals, first overtones) for the F-(H(sub 2)O) and F-(H(sub 2)O)2 clusters computed at the MP2 and CCSD(T) levels of theory with basis sets of triple zeta quality. Anharmonic corrections were estimated via the correlation-corrected vibrational self-consistent field (CC-VSCF) method. The CC-VSCF anharmonic spectra obtained on the potential energy surfaces evaluated at the CCSD(T) level of theory are the first ones reported at a correlated level beyond MP2. We have found that the average basis set effect (TZP vs. aug-cc-pVTZ) is on the order of 30-40 cm(exp -1), whereas the effects of different levels of electron correlation [MP2 vs. CCSD(T)] are smaller, 20-30 cm(exp -1). However, the basis set effect is much larger in the case of the H-bonded O-H stretch of the F-(H(sub 2)O) cluster amounting to 100 cm(exp -1) for the fundamentals and 200 cm (exp -1) for the first overtones. Our calculations are in agreement with the limited available set of experimental data for the F-(H(sub 2)O) and F-(H(sub 2)O)2 systems and provide additional information that can guide further experimental studies.

Chaban, Galina M.↗
NASA NTRS2004

Vibrational Spectroscopy of Ions and Radicals Present in the Interstellar Medium and in Planetary Atmospheres: A Theoretical Study

Anharmonic vibrational frequencies and intensities are calculated for OH(H2O)n and H(H2O)n radicals (that form on icy particles of the interstellar medium), HCO radical (the main intermediate in the synthesis of organic molecules in space), NH2(-) and C2H(-) anions, H5(+) cation, and other systems relevant to interstellar chemistry. In addition to pure ions and radicals, their complexes with water are studied to assess the effects of water environment on infrared spectra. The calculations are performed using the correlation-corrected vibrational self-consistent field (CC-VSCF) method with ab initio potential surfaces at the MP2 and CCSD(T) levels. Fundamental, overtone, and combination excitations are computed. The results are in good agreement with available experimental data and provide reliable predictions for vibrational excitations not yet measured in laboratory experiments. The data should be useful for interpretation of astronomically observed spectra and identification of ions and radicals present in the interstellar medium and in planetary atmospheres.

Chaban, Galina M.↗