Search NASA⌕ Search

SEARCH · Search NASA

Results for “open-shell”

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 19 records

Open-shell coupled-cluster theory

An efficient formulation of a recently proposed open-shell singles and doubles coupled-cluster (OCCSD) method is presented. This formulation is in terms of spatial orbital one- and two-electron integrals. Our new OCCSD method is based on 'symmetric spin orbitals' and is thus symmetric (or antisymmetric) in the spin indices. It therefore contains about half the number of independent parameters in the coupled-cluster wave function compared to other open-shell CCSD methods. It is shown that the formulation presented here contains less than half the number of n exp 6 steps (where n is the number of molecular orbitals) of other recently proposed open-shell CCSD methods. A new approach by which amplitudes in our method may be compared with amplitudes in a previous OCCSD method is examined.

Jayatilaka, Dylan↗

The form of spin orbitals for open-shell restricted Hartree-Fock reference functions

A new set of spin orbitals is proposed for single-reference correlation procedures based on an open-shell restricted Hartree-Fock (RHF) reference function. The spatial parts of the new spin orbitals are from an open-shell RHF calculation. These spin orbitals form a complete, orthonormal set. The RHF single determinant based on these spin orbitals is an S sub 2 eigenfunction, but a linear combination of S sub z eigenfunctions. The advantage of these spin orbitals is that equations for perturbation theory or coupled-cluster theory are symmetric in the spin indices and hence the number of independent parameters to be determined is dramatically reduced.

Jayatilaka, Dylan↗

An open-shell restricted Hartree-Fock perturbation theory based on symmetric spin orbitals

A new open-shell perturbation theory is formulated in terms of symmetric spin orbitals. Only one set of spatial orbitals is required, thereby reducing the number of independent coefficients in the perturbed wavefunctions. For second order, the computational cost is shown to be similar to a closed-shell calculation. This formalism is therefore more efficient than the recently developed RMP, ROMP or RMP-MBPT theories. The perturbation theory described herein was designed to have a close correspondence with our recently proposed coupled-cluster theory based on symmetric spin orbitals. The first-order wavefunction contains contributions from only doubly excited determinants. Equilibrium structures and vibrational frequencies determined from second-order perturbation theory are presented for OH, NH, CH, 02, NH2 and CH2.

Lee, Timothy J.↗

Development of New Open-Shell Perturbation and Coupled-Cluster Theories Based on Symmetric Spin Orbitals

A new spin orbital basis is employed in the development of efficient open-shell coupled-cluster and perturbation theories that are based on a restricted Hartree-Fock (RHF) reference function. The spin orbital basis differs from the standard one in the spin functions that are associated with the singly occupied spatial orbital. The occupied orbital (in the spin orbital basis) is assigned the delta(+) = 1/square root of 2(alpha+Beta) spin function while the unoccupied orbital is assigned the delta(-) = 1/square root of 2(alpha-Beta) spin function. The doubly occupied and unoccupied orbitals (in the reference function) are assigned the standard alpha and Beta spin functions. The coupled-cluster and perturbation theory wave functions based on this set of "symmetric spin orbitals" exhibit much more symmetry than those based on the standard spin orbital basis. This, together with interacting space arguments, leads to a dramatic reduction in the computational cost for both coupled-cluster and perturbation theory. Additionally, perturbation theory based on "symmetric spin orbitals" obeys Brillouin's theorem provided that spin and spatial excitations are both considered. Other properties of the coupled-cluster and perturbation theory wave functions and models will be discussed.

Lee, Timothy J.↗

Recent Advances in Open-Shell Perturbation Theory and Coupled-Cluster Theory

Comparisons of various recently developed open-shell RHF perturbation theories will be presented. Among the aspects considered are spin-contamination, computational cost, and quality of numerical results. In addition, a new approach to avoid the disk storage and I/O bottlenecks in large scale coupled-cluster calculations will be discussed.

Lee, Timothy J.↗

Photoionization cross sections for atomic chlorine using an open-shell random phase approximation

The use of the Random Phase Approximation with Exchange (RPAE) for calculating partial and total photoionization cross sections and photoelectron angular distributions for open shell atoms is examined for atomic chlorine. Whereas the RPAE corrections in argon (Z=18) are large, it is found that those in chlorine (Z=17) are much smaller due to geometric factors. Hartree-Fock calculations with and without core relaxation are also presented. Sizable deviations from the close coupling results of Conneely are also found.

Starace, A. F.↗

Ionization processes in collisions of open-shell atoms. III - The autoionizing states of nitrogen

Results of a study of the energy spectra of electrons produced in collisions of N atoms with inert gases at low keV energies are reported. Ionization here is partly due to production of the (1Dnl) autoionizing states of nitrogen and partly due to another mechanism, which is presumably quasi molecular Auger ionization. A discussion of the assignments of the autoionizing states is presented.

Boumsellek, S.↗

On the uniqueness of the constrained space orbital variation (CSOV) technique

Several CSOV analyses are performed for the 1Sigma(+) state of NiCO, and it is shown that the importance of the CO sigma donation, Ni pi back donation, and interunit polarizations are virtually independent of the order of the CSOV steps, provided that the open-shell 3d sigma and 4s Ni orbitals are orthogonalized to the CO. This order of orthogonalization is consistent with the polarization of the Ni observed in the unconstrained SCF wavefunction. If instead the CO is orthogonalized to the open-shell Ni orbitals, the frozen orbital repulsion and entire CSOV analysis becomes unphysical. A comparison of the SCF and CAS SCF descriptions for the NiCO 1Sigma(+) state shows the importance of the s to d promotion and sd hybridization in reducing the repulsion and increasing the Ni to CO pi bonding. For LiF, CSOV analyses starting from both the neutral and ionic asymptotes show the bonding to be predominantly Li(+) - F(-). These examples show the uniqueness of the CSOV decomposition.

Bauschlicher, C. W., Jr.↗

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

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

Approximate spin projection of three-component UHF wavefunctions - The states of the pentachlorocyclopentadienyl cation and the croconate dianion, C5O5/2-/

The approximate spin projection method of Amos et al. is extended to handle UHF wave functions having three significant components of differing multiplicity. An expression is given for the energy after single annihilation which differs from that of Amos and Hall. The new expression reproduces the results obtained from a previous exact calculation for which the weights and energies of the components are known. The extended approximate projection method is applied to the pi-electron UHF wave functions for the ground states of the pentachlorocyclopentadienyl cation and the croconate dianion, C5O5(2-). The results indicate a triplet ground state for the former and a singlet ground state for the latter, in agreement with experimental ESR susceptibility measurements for these molecular ions. C5C15(-) cannont be treated by restricted Hartree-Fock theory, due to its open-shell ground state. Incorrect results are obtained for the croconate dianion, if restricted Hartree-Fock theory and singly excited configuration interactions are utilized.

Phillips, D. H.↗

Photoelectron angular distributions as a probe of anisotropic electron-ion interactions

Expressions are given for atomic photoelectron angular distributions in LS coupling in which the role of anisotropic final state electron-ion interactions emerges explicitly. Calculations of photoelectron angular distributions for atomic sulfur are presented in which these anisotropic interactions produce clear deviations from the predictions of the Cooper-Zare model. Such effects are expected to be a general feature of photoelectron angular distributions for most open-shell atoms.

Dill, D.↗

Theoretical studies of diatomic and triatomic systems containing the group IB atoms Cu, Ag, and Au

Selected portions of the ground state potential energy surfaces of the Cu3, Ag3, AgCu2, and AuCu2 trimers are studied at the single-reference singles plus doubles configuration interaction and couple pair functional levels correlating 33 electrons. The calculations use the effective core potentials of Hay and Wadt (1985) to replace the deep core levels. The Cu3 and Ag3 molecules are found to have 2B2 obtuse-angled ground states, with low-lying 2A1 acute-angled excited states. The AgCu2 and AuCu2 molecules have 2A1 acute-angled ground states. The Cu3 molecule has a smaller 3d population than Ag3, and Cu3 has a smaller s electron density in the open-shell orbital than does Ag3, in agreement with recent ESR experiments.

Walch, Stephen P.↗

Theoretical study of the homonuclear tetramers and pentamers of the group IB metals (Cu, Ag, and Au)

The optimal structures of the homonuclear tetramers and pentamers of Cu, Ag, and Au are found to be planar trapezoidal. However, on the basis of spin densities deduced from electron spin resonance experiments, it has been suggested that these pentamers have distorted trigonal bipyramidal structures. This apparent discrepancy is reconciled by the fact that Mulliken populational analyses indicate that the open-shell spin density distribution on the planar structures are also consistent with experiment. Since the vertical electron attachment energies for both the tetramers and pentamers are in good agreement with experiment, the negative ions probably also have planar trapezoidal structures. The trends in ionization potentials, atomization energies, and electron affinities with cluster size are discussed.

Bauschlicher, Charles W., Jr.↗

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