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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

NASA NTRS1995

The Ozonide Anion: A Theoretical Study

Complete active space self-consistent field (CASSCF) and CASSCF second order perturbation theory (CASPT2) methods have been used for the geometry optimization and calculation of harmonic and fundamental frequencies of the ozonide ion O3(-) and the ozonide lithium complex Li(+)O3(-). For O3(-) harmonic frequencies have also been obtained using the coupled-cluster method, CCSD(T). Infrared intensities are computed from dipole moment derivatives at the CASSCF level. The predicted equilibrium geometry for O3(-) is: Roo = 1.361 Angstroms and angle ooo = 115.4 degrees, and the fundamental frequencies are: nu(sub 1) = 989 per centimeter, nu(sub 2) 556 per centimeter, nu(sub 3) = 870 per centimeter (experimental values are: Roo =1.36 plus or minus 0.02 Angstroms, angle (ooo) = 111.8 plus or minus 2.0 degrees, nu(sub 1) = 975(50) per centimeter, nu(sub 2) = 550(50) per centimeter, nu(sub 3) = 880(50) per centimeter). Corresponding data for the lithium ozonide complex have also been obtained. The presented data contradict the previous interpretation of the IR and Raman spectrum obtained after deposition of ozone in N2, argon, or neon matrices with atomic beams of alkali metals. The presence of the lithium cation raises the asymmetric stretch frequency to about 940 per centimeter, which is contradictory to assumptions made in the assignments of the matrix spectra. Calculations made in a dielectric medium strongly suggest that the effect of the matrix on the IR spectrum is small for O3(-) itself. The dissociation and atomization energies of O3(-) are found to be in agreement with experiment.

Borowski, Poitr↗

NASA NTRS1991

An ab initio study of the C3(+) cation using multireference methods

The energy difference between the linear 2 sigma(sup +, sub u) and cyclic 2B(sub 2) structures of C3(+) has been investigated using large (5s3p2d1f) basis sets and multireference electron correlation treatments, including complete active space self consistent fields (CASSCF), multireference configuration interaction (MRCI), and averaged coupled-pair functional (ACPF) methods, as well as the single-reference quadratic configuration interaction (QCISD(T)) method. Our best estimate, including a correction for basis set incompleteness, is that the linear form lies above the cyclic from by 5.2(+1.5 to -1.0) kcal/mol. The 2 sigma(sup +, sub u) state is probably not a transition state, but a local minimum. Reliable computation of the cyclic/linear energy difference in C3(+) is extremely demanding of the electron correlation treatment used: of the single-reference methods previously considered, CCSD(T) and QCISD(T) perform best. The MRCI + Q(0.01)/(4s2p1d) energy separation of 1.68 kcal/mol should provide a comparison standard for other electron correlation methods applied to this system.

Taylor, Peter R.↗

NASA NTRS1991

The computation of C-C and N-N bond dissociation energies for singly, doubly, and triply bonded systems

The bond dissociation energies (De) of C2H2, C2H4, C2H6, N2, N2H2, and N2H4 are studied at various levels of correlation treatment. The convergence of De with respect to the one-particle basis is studied at the single-reference modified coupled-pair-functional (MCPF)level. At all levels of correlation treatment, the errors in the bond dissociation energies increase with the degree of multiple bond character. The multireference configuration-interaction (MRCI) De values, corrected for an estimate of higher excitations, are in excellent agreement with those determined using the size-extensive averaged-coupled-pair-functional (ACPF) method. The full-valence complete-active-space self-consistent-field (CASSCF)/MRCI calculations are reproduced very well by MRCI calculations based on a CASSCF calculation that includes in the active space only those electrons involved in the C-C or N-N bonds.

Langhoff, Stephen R.↗

NASA NTRS1985

Fast and optimal solution to the Rankine-Hugoniot problem

A new, definitive, reliable and fast iterative method is described for determining the geometrical properties of a shock (i.e., theta sub Bn, yields N, V sub s and M sub A), the conservation constants and the self-consistent asymptotic magnetofluid variables, that uses the three dimensional magnetic field and plasma observations. The method is well conditioned and reliable at all theta sub Bn angles regardless of the shock strength or geometry. Explicit proof of uniqueness of the shock geometry solution by either analytical or graphical methods is given. The method is applied to synthetic and real shocks, including a bow shock event and the results are then compared with those determined by preaveraging methods and other iterative schemes. A complete analysis of the confidence region and error bounds of the solution is also presented.

Vinas, A. F.↗

NASA NTRS1986

Fast and optimal solution to the 'Rankine-Hugoniot problem'

Vinas, A. F.↗

NASA NTRS1989

The computation of C-C and N-N bond dissociation energies for singly, doubly, and triply bonded systems

The bond dissociation energies (D sub e) of C2H2, C2H4, C2H6, N2, N2H2, and N2H4 are studied at various levels of correlation treatment. The convergence of D sub e with respect to the one particle basis is studied at the single reference modified coupled-pair functional (MCPF) level. At all levels of correlation treatment, the errors in the bond dissociation energies increase with the degree of multiple bond character. The multireference configuration interaction (MRCI) D sub e values, corrected for an estimate of higher excitations, are in excellent agreement with those determined using the size extensive averaged coupled pair functional (ACPF) method. It was found that the full valence complete active space self consistent field (CASSCF)/MRCI calculations are reproduced very well by MRCI calculations based on a CASSCF calculation that includes in the active space only those electrons involved in the C-C or N-N bonds. To achieve chemical accuracy (1 kcal/mole) for the D sub e values of the doubly bonded species C2H4 and N2H2 requires one particle basis sets including up through h angular momentum functions (l = 5) and a multireference treatment of electron correlation: still higher levels of calculation are required to achieve chemical accuracy for the triply bonded species C2H2 and N2.

Langhoff, Stephen R.↗

NASA NTRS1989

Three-dimensional equilibria for the extended magnetotail and the generation of field-aligned current sheets

Using the magnetotail equilibrium theory and a solution method outlined by Birn (1987), self-consistent three-dimensional models were constructed for the quiet average magnetotail beyond about 20 RE distance but earthward of a potential distant neutral line, which take into account the decrease of the tail flaring with distance. It is found that this effect is coupled with the presence of magnetic shear and thus with field-aligned electric currents. These currents have the signature of region 1 currents, toward the earth on the dawnside and away on the duskside, and contribute about 500,000 A to the total Birkeland current.

Birn, J.↗

NASA NTRS2018

Overview of Methods for Deriving the Radiative Transfer Theory from the Maxwell Equations. I: Approach Based on the Far-Field Foldy Equations

In this paper, we revisit, with further enhancements and clarifications, the self-consistent first-principles approach developed previously for deriving the vector radiative transfer theory for a discrete random medium with a sparse concentration of particles. We specifically consider the case of a plane-parallel particulate layer embedded in an otherwise homogeneous unbounded medium. The solution method is based on the far-field Foldy equations, an order-of-scattering expansion for the total field derived under the Twersky approximation, the computation of the coherent field by assuming that the positions of the particles are uncorrelated, and the ladder approximation for the coherency dyadic. The latter yields an integral equation for the diffuse specific coherency dyadic, defined through an angular spectrum repre- sentation for the coherency dyadic, which in turn, gives the vector radiative transfer equation for the diffuse specific intensity column vector. We analyze specifically the computation of the coherent field for inhomogeneous particulate media and multiple species of particles, the continuous extension of the far-field representation to the near field, the Foldy approximation, and the Foldy integral equation for the coherent field. Finally, we discuss the transition from the vector to the scalar radiative transfer equation.

Doicu, Adrian↗

NASA NTRS1994

Relativistic and Correlation Effects in CuH, AgH and AuH: Comparison of Various Relativistic Methods

The effects of relativity on the bond lengths, dissociation energies, and harmonic vibrational frequencies of the 1Epsilon(+) electronic ground states of the group IB hydrides CuH, AgH and AuH have been evaluated with a variety of ab initio methods. These properties were investigated with moderately-sized basis sets at the self-consistent field Hartree Fock (SCF HF) level and with second-order Moller-Plesset (MP2) perturbation theory for electron correlation. Comparisons were made between all-electron results using the nonrelativistic Hamiltonian, perturbation theory (PT) at first-order with only the one-electron non-fine structure terms of the Breit-Pauli Hamiltonian, the spin-free Douglas-Kroll (DK) transformed Dirac Hamiltonian and the untransformed Dirac Hamiltonian, and results using two sets of relativistic effective core potentials (RECPs). The expected trends of bond length decrease, dissociation energy increase and harmonic frequency increase with both relativity and correlation are found. Both sets of RECPs are shown to give good results, if accompanied by a reasonable basis set. The DK method is demonstrated to be an inexpensive, reliable approximation to the DHF method.

Collins, Charlene L.↗

NASA NTRS1992

Ab initio study of the molecular structure and vibrational spectrum of nitric acid and its protonated forms

The equilibrium structures, harmonic vibrational frequencies, IR intensities, and relative energetics of HNO3 and its protonated form H2NO3+ were investigated using double-zeta plus polarization and triple-zeta plus polarization basis sets in conjunction with high-level ab initio methods. The latter include second-order Moller-Plesset perturbation theory, the single and double excitation coupled cluster (CCSD) methods, a perturbational estimate of the effects of connected triple excitations (CCSD(T)), and the self-consistent field. To determine accurate energy differences CCSD(T) energies were computed using large atomic natural orbital basis sets. Four different isomers of H2NO3+ were considered. The lowest energy form of protonated nitric acid was found to correspond to a complex between H2O and NO2+, which is consistent with earlier theoretical and experimental studies.

Lee, Timothy J.↗

NASA NTRS2012

Plasma Diffusion in Self-Consistent Fluctuations

The problem of particle diffusion in position space, as a consequence ofeleclromagnetic fluctuations is addressed. Numerical results obtained with a self-consistent hybrid code are presented, and a method to calculate diffusion coefficient in the direction perpendicular to the mean magnetic field is proposed. The diffusion is estimated for two different types of fluctuations. The first type (resuiting from an agyrotropic in itiai setting)is stationary, wide band white noise, and associated to Gaussian probability distribution function for the magnetic fluctuations. The second type (result ing from a Kelvin-Helmholtz instability) is non-stationary, with a power-law spectrum, and a non-Gaussian probabi lity distribution function. The results of the study allow revisiting the question of loading particles of solar wind origin in the Earth magnetosphere.

Smets, R.↗

NASA NTRS1997

Self-Consistent Superthermal Electron Effects on Plasmaspheric Refilling

The effects of self-consistently including superthermal electrons in the definition of the ambipolar electric field are investigated for the case of plasmaspheric refilling after a geomagnetic storm. By using the total electron population in the hydrodynamic equations, a method for incorporating superthermal electron parameters in the electric field and electron temperature calculation is developed. Also, the ambipolar electric field is included in the kinetic equation for the superthermal electrons through a change of variables using the total energy and the first adiabatic invariant. Calculations based on these changes are performed by coupling time-dependent models of the thermal plasma and superthermal electrons. Results from this treatment of the electric field and the self-consistent development of the solution are discussed in detail. Specifically, there is a decreased thermal electron density in the plasmasphere during the first few minutes of refilling, a slightly accelerated proton shock front, and a decreased superthermal electron flux due to the deceleration by the electric field. The timescales of plasmaspheric refilling are discussed and determined to be somewhat shorter than previously calculated for the thermal plasma and superthermal electron population due to the effects of the field-aligned potential.

Liemohn, M. W.↗

NASA NTRS1997

The Vibrational Frequencies of CaO2, ScO2, and TiO2: A Comparison of Theoretical Methods

The vibrational frequencies of several states of CaO2, ScO2, and TiO2 are computed at using density functional theory (DFT), the Hatree-Fock approach, second order Moller-Plesset perturbation theory (MP2), and the complete-active-space self-consistent-field theory. Three different functionals are used in the DFT calculations, including two hybrid functionals. The coupled cluster singles and doubles approach including the effect of unlinked triples, determined using perturbation theory, is applied to selected states. The Becke-Perdew 86 functional appears to be the cost effective method of choice, although even this functional does not perform well for one state of CaO2. The MP2 approach is significantly inferior to the DFT approaches.

Rosi, Marzio↗

NASA NTRS1986

Theoretical research program to study chemical reactions in AOTV bow shock tubes

Progress in the development of computational methods for the characterization of chemical reactions in aerobraking orbit transfer vehicle (AOTV) propulsive flows is reported. Two main areas of code development were undertaken: (1) the implementation of CASSCF (complete active space self-consistent field) and SCF (self-consistent field) analytical first derivatives on the CRAY X-MP; and (2) the installation of the complete set of electronic structure codes on the CRAY 2. In the area of application calculations the main effort was devoted to performing full configuration-interaction calculations and using these results to benchmark other methods. Preprints describing some of the systems studied are included.

Taylor, P.↗

NASA NTRS1987

The role and behavior of spin in gravitational physics

A self-consistent method of introducing spin into any Lagrangian based theory of gravitation was developed. The metric variation of the Lagrangian in the theory leads to an improved energy-momentum tensor which represents the source term in the gravitational field equations. The goal of the research is the construction of a theory general enough to be used to investigate spin effects in astrophysical objects and cosmology, and also to serve as a basis for discussion of the theoretical ideas tested by the NASA Gyroscope Experiment (aboard Gravity Probe B). Specific accomplishments in the following areas are summarized: the inclusion of electromagnetism into the variational principle for spinning matter, formulation of a self-consistent theory for the case of a fluid in which particle production processes occur, and the derivation of the Raychaudhuri equation in the case of spinning matter.

Ray, John R.↗

NASA NTRS1986

Benchmark full configuration-interaction calculations on H2O, F- and F

Full configuration-interaction calculations are reported, and compared to other methods, for H2O at its equilibrium geometry and at two geometries with the H-O bonds stretched. Since the percentage of the self-consistent field (SCF) reference in the full configuration-interaction (FCI) wave function decreases greatly with the bond elongation, the accuracy of techniques based on a single reference do not compare well with the FCI results. However, the results from a complete active space SCF/multireference configuration-interaction (CASSCF/MRCI) treatment are in good agreement with the FCI. Correlation effects in F compared to Ne are far more similar than for F- compared to Ne, despite F- and Ne being isoelectronic. Since the importance of higher than double excitations is more important for F- than F, a very high percentage of the correlation must be obtained to accurately compute the electron affinity. In a CASSCF/MRCI treatment the higher than quadruple excitations contribute 0.02 eV to the electron affinity (EA), even for modest basis sets.

Bauschlicher, C. W., Jr.↗

NASA NTRS1986

Benchmark full configuration-interaction calculations on H2O, F, and

Full configuration-interaction calculations are reported, and compared to other methods, for H2O at its equilibrium geometry and at two geometries with the H-O bonds stretched. Since the percentage of the self-consistent field (SCF) reference in the full configuration-interaction (FCI) wave function decreases greatly with the bond elongation, the accuracy of techniques based on a single reference do not compare well with the FCI results. However, the results from a complete active space SCF/multireference configuration-interaction (CASSCF/MRCI) treatment are in good agreement with the FCI. Correlation effects in F compared to Ne are far more similar than for Fcompared to Ne, despite F- and Ne being isoelectronic. Since the importance of higher than double excitations is more important for F- than F, a very high percentage of the correlation must be obtained to accurately compute the electron affinity. In a CASSCF/MRCI treatment the higher than quadruple excitations contribute 0.02 eV to the electron affinity (EA), even for modest basis sets.

Bauschlicher, C. W., Jr.↗

NASA NTRS1993

A new look at rainfall fluctuations and scaling properties of spatial rainfall using orthogonal wavelets

It has been observed that the finite-dimensional distribution functions of rainfall cannot obey simple scaling laws due to rainfall intermittency (mixed distribution with an atom at zero) and the probability of rainfall being an increasing function of area. Although rainfall fluctuations do not suffer these limitations, it is interesting to note that very few attempts have been made to study them in terms of their self-similarity characteristics. This is due to the lack of unambiguous definition of fluctuations in multidimensions. This paper shows that wavelet transforms offer a convenient and consistent method for the decomposition of inhomogeneous and anisotropic rainfall fields in two dimensions and that the components of this decomposition can be looked at as fluctuations of the rainfall field. It is also shown that under some mild assumptions, the component fields can be treated as homogeneous and thus are amenable to second-order analysis, which can provide useful insight into the nature of the process. The fact that wavelet transforms are a space-scale method also provides a convenient tool to study scaling characteristics of the process. Orthogonal wavelets are used, and these properties are investigated for a squall-line storm to study the presence of self-similarity.

Kumar, Praveen↗