Engineering topics

Dyall, Kenneth G.

Publications and source records attributed to Dyall, Kenneth G..

At least 19 records

DOE OSTI2024

A new computational framework for spinor-based relativistic exact two-component calculations using contracted basis functions

Here, a new computational framework for spinor-based relativistic exact two-component (X2C) calculations is developed using contracted basis sets with a spin–orbit contraction scheme. Generally contracted, j-adapted basis sets of p-block elements using primitive functions in the correlation-consistent basis sets are constructed for the X2C Hamiltonian with atomic mean-field spin–orbit integrals (the X2CAMF scheme). The contraction coefficients are taken from atomic X2CAMF Hartree–Fock spinors, thereby following the simple concept of a linear combination of atomic orbitals. Benchmark calculations of spin–orbit splittings, equilibrium bond lengths, and harmonic vibrational frequencies demonstrate the accuracy and efficacy of the j-adapted spin–orbit contraction scheme.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2001

Is the Lamb shift chemically significant?

The contribution of the Lamb shift to the atomization energies of some prototype molecules, BF3, AlF3, and GaF3, is estimated by a perturbation procedure. It is found to be in the range of 3-5% of the one-electron scalar relativistic contribution to the atomization energy. The maximum absolute value is 0.2 kcal/mol for GaF3. These sample calculations indicate that the Lamb shift is probably small enough to be neglected for energetics of molecules containing light atoms if the target accuracy is 1 kcal/mol, but for higher accuracy calculations and for molecules containing heavy elements it must be considered.

Dyall, Kenneth G.↗

NASA NTRS1999

Bond Dissociation Energies of the Tungsten Fluorides and Their Singly-Charged Ions: A Density Functional Survey

The dissociation of WF6 and the related singly-charged cations and anions into the lower fluorides and fluorine atoms has been investigated theoretically using density functional theory (B3LYP) and relativistic effective core potentials, with estimates of spin-orbit effects included using a simple model. The inclusion of spin-orbit is essential for a correct description of the thermochemistry. The total atomization energy of the neutral and anionic WF6 is reproduced to within 25 kcal/mol, but comparison of individual bond dissociation energies with available experimental data shows discrepancies of up to 10 kcal/mol. The results are nevertheless useful to help resolve discrepancies in experimental data and provide estimates of missing data.

Dyall, Kenneth G.↗

NASA NTRS1995

MP2 Studies of Relativistic Effects on the Linear Stationary Points of the H+Cl2 Yields HCl+Cl and Cl+HCl Yields ClH+Cl Reactions

The influence of relativistic effects on the linear stationary points of the potential energy surfaces of the ClHCl and HCl2 systems is studied. Scalar relativistic effects have little influence on the geometries and the energies of these points. Spin-orbit effects have no influence on the barrier of the forward reaction of H with Cl2, but increase the reaction energy. The total relativistic effect is an increase of the reaction energy by 2.3 kJ /mol. The barrier for the iso-energetic exchange reaction of HCl with Cl is decreased by 0.5 kJ /mol due to the scalar relativistic effects. This small decrease is cancelled by the larger increase of 2.5 kJ /mol due to spin-orbit effects, giving a total relativistic effect of 2.0 kJ /mol.

Visscher, L.↗

NASA NTRS1995

Kramers' Restricted Closed Shell CCSD Theory

A Kramers' restricted version of the closed shell coupled cluster singles doubles theory is presented. The theory may be used in conjunction with 2 or 4-component relativistic reference wavefunctions. The intrinsic treatment of the spin-orbit coupling doubles the number of independent quantities (amplitudes and integrals) relative to a spin-independent formalism. The number of operations required to evaluate the equations is four times larger than in the optimal spin-independent closed shell formalism.

Visscher, Lucas↗

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 NTRS1994

Theoretical investigation of gas-surface interactions

The goal of this project was to develop computational tools for the calculation of the electronic structure of molecules containing heavy atoms, and to use these tools in the study of catalytic processes, with the overall objective of gaining an understanding of the catalytic process which could be used to design more efficient catalysts. The main catalytic system of interest was the combustion of hydrogen on platinum surfaces. Under this project, a flexible Dirac-Hartree-Fock (DHF) program has been developed, a code to calculate DHF correlation energies at the second-order Moeller-Plesset perturbation (PP2) level is almost complete, and code to include correlation at the MCSCF and MCSCF/MP2 level is planned. The tools so far developed have been validated and used to calibrate some more approximate methods, and applied to investigate the importance of relativistic effects in the bonding of hydrogen to platinum.

Dyall, Kenneth G.↗

NASA NTRS1994

An Exact Separation of the Spin-Free and Spin-Dependent Terms of the Dirac-Coulomb-Breit Hamiltonian

The Dirac Hamiltonian is transformed by extracting the operator (sigma x p)/2mc from the small component of the wave function and applying it to the operators of the original Hamiltonian. The resultant operators contain products of Paull matrices that can be rearranged to give spin-free and spin-dependent operators. These operators are the ones encountered in the Breit-Pauli Hamiltonian, as well as some of higher order in alpha(sup 2). However, since the transformation of the original Dirac Hamiltonian is exact, the new Hamiltonian can be used in variational calculations, with or without the spin-dependent terms. The new small component functions have the same symmetry properties as the large component. Use of only the spin-free terms of the new Hamiltonian permits the same factorization over spin variables as in nonrelativistic theory, and therefore all the post-Self-Consistent Field (SCF) machinery of nonrelativistic calculations can be applied. However, the single-particle functions are two-component orbitals having a large and small component, and the SCF methods must be modified accordingly. Numerical examples are presented, and comparisons are made with the spin-free second-order Douglas-Kroll transformed Hamiltonian of Hess.

Dyall, Kenneth G.↗

NASA NTRS1994

Theoretical Studies of PNC in Molecules: current Capabilites, Future Prospects, and Issues

The current capabilities of relativistic electronic structure codes for calculations on molecules such as TlF and YbF, which are of interest in the study of parity non-conservation, will be presented. The issues of basis set size and correlation methods will be discussed, and planned developments of the codes which would extend the range of their capabilities will be outlined. Issues of the sensitivity of the desired properties to the form of the nuclear potential and the flexibility of the wave function near the nucleus will be raised to stimulate discussion.

Dyall, Kenneth G.↗

NASA NTRS1994

On The Choice of a Zeroth-Order Hamiltonian for Second-Order Perturbation Theory with A CASSCF Reference Function

A new approach to perturbation theory based on a CASSCF reference function has been developed. The key to the approach is the definition of the zeroth order Hamiltonian, H(sub 0), which includes the full CI Hamiltonian for the active space. In the the inactive and secondary spaces, operators may be chosen which reduce to the usual Moller-Plesset or Epstein-Nesbet forms in the limit of a null active space. These operators are diagonal in the orbital indices and permit the block-diagonalization of H(sub 0). The reference is an eigenfunction of Ho without N-particle projection. H(sub 0) automatically incorporates denominator shifts in the style of those appearing in recent open-shell perturbation theories. Comparative results are presented for a few test cases.

Dyall, Kenneth G.↗

NASA NTRS1994

Second-Order Moller-Plesset Perturbation Theory for Molecular Dirac-Hartree-Fock Wave Functions

Moller-Plesset perturbation theory is developed to second order for a selection of Kramers restricted Dirac-Hartree-Fock closed and open-shell reference wave functions. The open-shell wave functions considered are limited to those with no more than two electrons in open shells, but include the case of a two-configuration SCF reference. Denominator shifts are included in the style of Davidson's OPT2 method. An implementation which uses unordered integrals with labels is presented, and results are given for a few test cases.

Dyall, Kenneth G.↗

NASA NTRS1994

Some Aspects of the Implementation of Double Group Symmetry and Electron Correlation in Molecular 4-Component Calculations

The efficient implementation of method for electron correlation in molecular 4-component calculations demands that symmetry be exploited where possible. Algorithms for the construction of matrices and the transformation of integrals over symmetry-adapted basis functions, where the point group is restricted to D(sub 2h) and subgroups, will be presented. The merits of keeping the primitive integrals in the scalar basis will be compared with those of transforming them to the 2-spinor basis.

Dyall, Kenneth G.↗

NASA NTRS1993

Relativistic effects on the bonding and properties of the hydrides of platinum

The ground state of PtH2 and several low-lying states of PtH(+) and PtH have been studied at the all-electron self-consistent-field level of theory to examine the importance of relativistic effects. The results of calculations based on Dirac-Hartree-Fock theory, nonrelativistic theory, and the spin-free no-pair relativistic approximation of Hess are compared to separate the effects of the spin-free terms and the spin-orbit terms of the Hamiltonian on the relativistic corrections to the molecular properties. Comparison is also made between first-order perturbation theory including the one-electron spin-free terms and the method of Hess to determine the size of effects beyond first order. It is found that the spin-orbit interaction significantly affects the properties and energetics of these molecules because of the participation of the Pt 5d orbitals in the bonding, and that effects beyond first order in perturbation theory are large. Any treatment of Pt compounds will have to include both the spin-free and spin-orbit interactions for an accurate description.

Dyall, Kenneth G.↗

NASA NTRS1993

Relativistic Corrections to the Properties of the Alkali Fluorides

Relativistic corrections to the bond lengths, dissociation energies and harmonic frequencies of KF, RbF and CsF have been obtained at the self-consistent field level by dissociating to ions. The relativistic corrections to the bond lengths, harmonic frequencies and dissociation energies to the ions are very small, due to the ionic nature of these molecules and the similarity of the relativistic and nonrelativistic ionic radii.

Dyall, Kenneth G.↗

NASA NTRS1993

All-electron molecular Dirac-Hartree-Fock calculations - Properties of the group IV monoxides GeO, SnO, and PbO

Dirac-Hartree-Fock calculations have been carried out on the ground states of the group IV monoxides GeO, SnO and PbO. Geometries, dipole moments and infrared data are presented. For comparison, nonrelativistic, first-order perturbation and relativistic effective core potential calculations have also been carried out. Where appropriate the results are compared with the experimental data and previous calculations. Spin-orbit effects are of great importance for PbO, where first-order perturbation theory including only the mass-velocity and Darwin terms is inadequate to predict the relativistic corrections to the properties. The relativistic effective core potential results show a larger deviation from the all-electron values than for the hydrides, and confirm the conclusions drawn on the basis of the hydride calculations.

Dyall, Kenneth G.↗

NASA NTRS1993

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.↗

NASA NTRS1992

The 2Sigma(+)-2Pi separation in KO

The 2Sigma(+)-2Pi separation in KO is investigated using large basis sets and high levels of correlation treatment. Relativistic effects are included at the Dirac-Fock level and reduce the separation only slightly. The basis set superposition error is considered in detail. On the basis of these calculations, our best estimate places the 2Pi3/2 state about 200/cm above the ground 2Sigma(+) state, in agreement with our previous estimate.

Bauschlicher, Charles W., Jr.↗

NASA NTRS1992

Theoretical investigation of gas-surface interactions

The investigation into the appearance of intruder states from the negative continuum when some of the two-electron integrals were omitted was completed. The work shows that, provided all integrals involving core contracted functions in an atomic general contraction are included, or that the core functions are radially localized, meaningful results are obtained and intruder states do not appear. In the area of program development, the Dirac-Hartree-Fock (DHF) program for closed-shell polyatomic molecules was extended to permit Kramers-restricted open-shell DHF calculations with one electron in an open shell or one hole in a closed shell, or state-averaged DHF calculations over several particle or hole doublet states. One application of the open-shell code was to the KO molecule. Another major area of program development is the transformation of integrals from the scalar basis in which they are generated to the 2-spinor basis employed in parts of the DHF program, and hence to supermatrix form. Particularly concerning the omission of small component integrals, and with increase in availability of disk space, it is now possible to consider transforming the integrals. The use of ordered integrals, either in the scalar basis or in the 2-spinor basis, would considerably speed up the construction of the Fock matrix, and even more so if supermatrices were constructed. A considerable amount of effort was spent on analyzing the integral ordering and tranformation for the DHF program. The work of assessing the reliability of the relativistic effective core potentials (RECPs) was continued with calculation of the group IV monoxides. The perturbation of the metal atom provided by oxygen is expected to be larger than that provided by hydrogen and thus provide a better test of the qualification of the RECPs. Calculations on some platinum hydrides were carried out at nonrelativistic (NR), perturbation theory (PT) and DHF levels. Reprints of four papers describing this work are included.

Dyall, Kenneth G.↗