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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
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.↗
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.↗
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 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 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.↗
NASA NTRS2005

Hydrogen Bonding and Vibrational Spectroscopy: A Theoretical Study

Effects of hydrogen bonding on vibrational spectra are studied for several hydrogen-bonded complexes, in which hydrogen bonding ranges from weak (<5 kcal/mol) to very strong (>25 kcal/mol). The systems studied include complexes of inorganic acids and salts with water and ammonia, as well as complexes of several organic molecules (nitriles and amino acids) with water. Since anharmonic effects are very strong in hydrogen-bonded systems, anharmonic vibrational frequencies and infrared intensities are computed using the correlation-corrected vibrational self-consistent field (CC-VSCF) method with ab initio potential surfaces at the MP2 and CCSD(T) levels. The most common spectral effects induced by hydrogen bonding are red shifts of stretching vibrational frequencies ranging from approx.200/cm to over 2000/cm and significant increases of infrared intensities for those bonds that participate in hydrogen bonding. However, some systems (e.g. nitrile-water complexes) exhibit shifts in the opposite direction (to the blue) upon formation of hydrogen bonds.

Chaban, Galina M.↗
NASA NTRS1991

Theoretical study of the NO beta system

A theoretical determination of the transition moment functions (TMFs) for the beta system of NO is presented. High levels of correlation treatment are required to show the changing degree of Rydberg character in the B2II with decreasing r values. The state-averaged complete-active-space self-consistent-field multireference configuration-interaction method is used for the determination. Previous lifetime measurements made with laser-induced fluorescence, varying between 2 and 0.85 microns, are discussed in terms of the calculated lifetimes for v-prime values 0-6, which vary from 2.12-1.17 microns. When larger r values are used for the transition moment function, the calculated lifetimes correlate with experimental lifetimes. The Einstein coefficients agree with experimental results, although limitations in the calibration of the spectrometer can account for systematic differences. The correlation with earlier experimental results suggests that radiative lifetimes are in the range of 1-2 microns.

Langhoff, Stephen R.↗
NASA NTRS1992

H-N2 interaction energies, transport cross sections, and collision integrals

The energies for the interaction of a hydrogen atom with a nitrogen molecule have been calculated for large separation distances using a complete-active-space self-consistent-field/externally contracted configuration interaction method. H-N2 transport cross sections and collision integrals have been calculated using sudden approximations and a semiclassical description of the scattering. The values of these quantities are found to be close to the corresponding values determined from the average (isotropic) potential energy. The collision integrals are applied to determine diffusion and viscosity coefficients; the theoretical diffusion agrees well with the measured data available from experiments at low temperatures.

Stallcop, James R.↗
NASA NTRS1991

Theoretical study of the spectroscopy of the alkali oxides LiO, NaO, and KO

A state-averaged complete-active-space self-consistent-field multireference configuration-interaction method is presented to characterize the bound-bound emission from the CPi-2 state into the two lowest ionic states of LiO, NaO, and KO. Ab initio calculations use the experimental results obtained by Woodward et al. (1989) of the emitting state as CPi-2, but indicate that the tentative experimental band assignments are incorrect.

Langhoff, Stephen R.↗
DOE OSTI2025

Bridging the gap between molecules and materials in quantum chemistry with localized active spaces

The number of materials that “bridge the gap” between single molecules and extended solids, such as metal-organic frameworks and organic semiconductors, has been increasing. Consequently, there is a growing need for modeling approaches that effectively integrate the real-space molecular perspective employed by computational chemists and the reciprocal-space dispersive perspective employed by computational physicists. Here, we propose the localized active space (LAS) approach as a promising method to successfully bridge this gap. The LAS approach extends the active space concept from multiconfigurational methods such as complete active space self-consistent field theory to multiple molecular fragments via a product-form wave function ansatz. Here, we apply this method to solid state phenomena by treating each unit cell as a fragment with different sets of local quantum numbers (e.g., charge and excitation number). State interaction between these LAS states (LASSI) thus provides a comprehensive basis for the study of charge and energy transfer, meeting and surpassing the capabilities of single-reference fragmentation approaches such as constrained density functional theory (cDFT). Most centrally, we show how combining this LASSI approach with multiconfigurational pair-density functional theory (MC-PDFT) provides an elegant and efficient method to compute band structures that capture multiconfigurational character. We apply the LASSI band structure approach to the computation of band gaps in stretched hydrogen chain, polyacetylene, and bulk nickel oxide (NiO), finding good or excellent quantitative agreement with reference values in all cases. Additionally, we use the LAS basis in one-dimensional model systems to demonstrate its ability to treat difficult solid-state phenomena such as exciton transfer and excitation at p-n junctions.

method development↗