Search NASA⌕ Search

Engineering topics

Faegri, Knut, Jr.

Publications and source records attributed to Faegri, Knut, Jr..

Finite nucleus effects on relativistic energy corrections

The effect of using a finite nucleus model in quantum-chemical calculations is examined. Relativistic corrections from the first order Foldy-Wouthuysen terms are affected indirectly by the change in wavefunction, but also directly as a result of revised expressions for the Darwin and spin-orbit terms due to the change in nuclear potential. A calculation for the Rn atom indicates that the mass-velocity and Darwin corrections are much more sensitive to the finite nucleus than the non-relativistic total energy, but that the total contribution for these two terms is quite stable provided the revised form of the Darwin term is used. The spin-orbit interaction is not greatly affected by the choice of nuclear model.

Dyall, Kenneth G.↗

High quality Gaussian basis sets for fourth-row atoms

Energy optimized Gaussian basis sets of triple-zeta quality for the atoms Rb-Xe have been derived. Two series of basis sets are developed: (24s 16p 10d) and (26s 16p 10d) sets which were expanded to 13d and 19p functions as the 4d and 5p shells become occupied. For the atoms lighter than Cd, the (24s 16p 10d) sets with triple-zeta valence distributions are higher in energy than the corresponding double-zeta distribution. To ensure a triple-zeta distribution and a global energy minimum, the (26s 16p 10d) sets were derived. Total atomic energies from the largest basis sets are between 198 and 284 (mu)E(sub H) above the numerical Hartree-Fock energies.

Partridge, Harry↗

All-electron molecular Dirac-Hartree-Fock calculations - The group IV tetrahydrides CH4, SiH4, GeH4, SnH4, and PbH4

A basis-set-expansion Dirac-Hartree-Fock program for molecules is described. Bond lengths and harmonic frequencies are presented for the ground states of the group 4 tetrahydrides, CH4, SiH4, GeH4, SnH4, and PbH4. The results are compared with relativistic effective core potential (RECP) calculations, first-order perturbation theory (PT) calculations and with experimental data. The bond lengths are well predicted by first-order perturbation theory for all molecules, but none of the RECP's considered provides a consistent prediction. Perturbation theory overestimates the relativistic correction to the harmonic frequencies; the RECP calculations underestimate the correction.

Dyall, Kenneth G.↗

Heteronuclear transition metal diatomics - The bonding and electronic structure of ScNi, YNi, ScPd, and YPd

High quality ab initio calculations show that ScNi, YNi, ScPd, and YPd all have 2Sigma(+) ground states in agreement with electron spin resonance experiments. For ScNi and YNi, this is expected based on the lowest atomic asymptote. For ScPd and YPd, the lowest atomic asymptote would give the order of stability 2Delta greater than 2Pi equal to about 2Sigma(+), but the calculations show that mixing in of the excited asymptotes preferentially lowers the 2Sigma(+) state. The calculations show that the quartet states are about 20-30 kcal/mol above the ground state, and therefore probably do not contribute significantly to the unexpected g(vertical) values found in experiment. Calculations of excited states for YPd reveal some strong transitions that should be amenable to spectroscopic studies.

Faegri, Knut, Jr.↗

Kinetic balance and variational bounds failure in the solution of the Dirac equation in a finite Gaussian basis set

The paper investigates bounds failure in calculations using Gaussian basis sets for the solution of the one-electron Dirac equation for the 2p1/2 state of Hg(79+). It is shown that bounds failure indicates inadequacies in the basis set, both in terms of the exponent range and the number of functions. It is also shown that overrepresentation of the small component space may lead to unphysical results. It is concluded that it is important to use matched large and small component basis sets with an adequate size and exponent range.

Dyall, Kenneth G.↗

All-electron molecular Dirac-Hartree-Fock calculations: The group 4 tetrahydrides CH4, SiH4, GeH4, SnH4 and PbH4

A basis-set-expansion Dirac-Hartree-Fock program for molecules is described. Bond lengths and harmonic frequencies are presented for the ground states of the group 4 tetrahydrides, CH4, SiH4, GeH4, SnH4, and PbH4. The results are compared with relativistic effective core potential (RECP) calculations, first-order perturbation theory (PT) calculations and with experimental data. The bond lengths are well predicted by first-order perturbation theory for all molecules, but non of the RECP's considered provides a consistent prediction. Perturbation theory overestimates the relativistic correction to the harmonic frequencies; the RECP calculations underestimate the correction.

Dyall, Kenneth G.↗

Polyatomic molecular Dirac-Hartree-Fock calculations with Gaussian basis sets

Numerical methods have been used successfully in atomic Dirac-Hartree-Fock (DHF) calculations for many years. Some DHF calculations using numerical methods have been done on diatomic molecules, but while these serve a useful purpose for calibration, the computational effort in extending this approach to polyatomic molecules is prohibitive. An alternative more in line with traditional quantum chemistry is to use an analytical basis set expansion of the wave function. This approach fell into disrepute in the early 1980's due to problems with variational collapse and intruder states, but has recently been put on firm theoretical foundations. In particular, the problems of variational collapse are well understood, and prescriptions for avoiding the most serious failures have been developed. Consequently, it is now possible to develop reliable molecular programs using basis set methods. This paper describes such a program and reports results of test calculations to demonstrate the convergence and stability of the method.

Dyall, Kenneth G.↗