Theoretical investigation of gas-surface interactions

A Dirac-Hartree-Fock code was developed for polyatomic molecules. The program uses integrals over symmetry-adapted real spherical harmonic Gaussian basis functions generated by a modification of the MOLECULE integrals program. A single Gaussian function is used for the nuclear charge distribution, to ensure proper boundary conditions at the nuclei. The Gaussian primitive functions are chosen to satisfy the kinetic balance condition. However, contracted functions which do not necessarily satisfy this condition may be used. The Fock matrix is constructed in the scalar basis and transformed to a jj-coupled 2-spinor basis before diagonalization. The program was tested against numerical results for atoms with a Gaussian nucleus and diatomic molecules with point nuclei. The energies converge on the numerical values as the basis set size is increased. Full use of molecular symmetry (restricted to D sub 2h and subgroups) is yet to be implemented

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

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

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

The 2S(+) - 2P separation in KO

The 2S(+) - 2P 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 2p sub 3/2 state about 200 cm(exp -1) above the ground 2 sigma(+) state in agreement with our previous estimate

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

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

Relativistic double-zeta, triple-zeta, and quadruple-zeta basis sets for the lanthanides La–Lu

Relativistic basis sets of double-zeta, triple-zeta, and quadruple-zeta quality have been optimized for the lanthanide elements La-Lu. The basis sets include SCF exponents for the occupied spinors and for the 6p shell, exponents of correlating functions for the valence shells (4f, 5d and 6s) and the outer core shells (4d, 5s and 5p), and diffuse functions, including functions for dipole polarization of the 4f shell. A finite nuclear size was used in all optimizations. The basis sets are illustrated by calculations on YbF. Prescriptions are given for constructing contracted basis sets. The basis sets are available as an internet archive and from the Dirac program web site, http://dirac. chem. sdu. dk. © 2010 The Author(s)

The DIRAC code for relativistic molecular calculations

DIRAC is a freely distributed general-purpose program system for one-, two-, and four-component relativistic molecular calculations at the level of Hartree?Fock, Kohn?Sham (including range-separated theory), multiconfigurational self-consistent-field, multireference configuration interaction, electron propagator, and various flavors of coupled cluster theory. At the self-consistent-field level, a highly original scheme, based on quaternion algebra, is implemented for the treatment of both spatial and time reversal symmetry. DIRAC features a very general module for the calculation of molecular properties that to a large extent may be defined by the user and further analyzed through a powerful visualization module. It allows for the inclusion of environmental effects through three different classes of increasingly sophisticated embedding approaches: the implicit solvation polarizable continuum model, the explicit polarizable embedding model, and the frozen density embedding model.Fil: Saue, Trond. Université Paul Sabatier; Francia. Centre National de la Recherche Scientifique; FranciaFil: Bast, Radovan. Uit The Arctic University Of Norway; NoruegaFil: Gomes, André Severo Pereira. University Of Lille.; Francia. Centre National de la Recherche Scientifique; FranciaFil: Jensen, Hans Jorgen Aa.. University of Southern Denmark; DinamarcaFil: Visscher, Lucas. Vrije Universiteit Amsterdam; Países BajosFil: Aucar, Ignacio Agustín. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Nordeste. Instituto de Modelado e Innovación Tecnológica. Universidad Nacional del Nordeste. Facultad de Ciencias Exactas Naturales y Agrimensura. Instituto de Modelado e Innovación Tecnológica; Argentina. Universidad Nacional del Nordeste. Facultad de Ciencias Exactas y Naturales y Agrimensura. Departamento de Física; ArgentinaFil: Di Remigio, Roberto. Uit The Arctic University of Norway; NoruegaFil: Dyall, Kenneth G.. Dirac Solutions; Estados UnidosFil: Eliav, Ephraim. Universitat Tel Aviv.; IsraelFil: Fasshauer, Elke. Aarhus University. Department of Bioscience; DinamarcaFil: Fleig, Timo. Université Paul Sabatier; Francia. Centre National de la Recherche Scientifique; FranciaFil: Halbert, Loïc. Centre National de la Recherche Scientifique; Francia. University Of Lille.; FranciaFil: Hedegård, Erik Donovan. Lund University; SueciaFil: Helmich-Paris, Benjamin. Max-planck-institut Für Kohlenforschung; AlemaniaFil: Ilias, Miroslav. Matej Bel University; EslovaquiaFil: Jacob, Christoph R.. Technische Universität Braunschweig; AlemaniaFil: Knecht, Stefan. Eth Zürich, Laboratorium Für Physikalische Chemie; SuizaFil: Laerdahl, Jon K.. Oslo University Hospital; NoruegaFil: Vidal, Marta L.. Department Of Chemistry; DinamarcaFil: Nayak, Malaya K.. Bhabha Atomic Research Centre; IndiaFil: Olejniczak, Malgorzata. University Of Warsaw; PoloniaFil: Olsen, Jógvan Magnus Haugaard. Uit The Arctic University Of Norway; NoruegaFil: Pernpointner, Markus. Kybeidos Gmbh; AlemaniaFil: Senjean, Bruno. Universiteit Leiden; Países BajosFil: Shee, Avijit. Department Of Chemistry; Estados UnidosFil: Sunaga, Ayaki. Tokyo Metropolitan University; JapónFil: van Stralen, Joost N. P.. Vrije Universiteit Amsterdam; Países Bajo