Quasi-bound states of coupled Morse oscillators

Quasi-bound (resonance) states are present in the continuous spectrum of the Hamiltonian of two coupled Morse oscillators. Two different methods for approximating these as localized states are compared. The algebraic approach is shown to be in very good accord with the other method which is formulated in coordinate space and hence is differential in character. For these highly excited states an intermultiplet mixing term is included in the algebraic Hamiltonian.     

Perturbational approach to the electron correlation cusp applied to helium‐like atoms

A recently proposed perturbational approach to the electron correlation cusp problem 1 is tested in the context of three spherically symmetrical two‐electron systems: helium atom, hydride anion, and a solvable model system. The interelectronic interaction is partitioned into long‐ and short‐range components. The long‐range interaction, lacking the singularities responsible for the electron correlation cusp, is included in the reference Hamiltonian. Accelerated convergence of orbital‐based methods for this smooth reference Hamiltonian is shown by a detailed partial wave analysis. Contracted orbital basis sets constructed from atomic natural orbitals are shown to be significantly better for the new Hamiltonian than standard basis sets of the same size. The short‐range component becomes the perturbation. The low‐order perturbation equations are solved variationally using basis sets of correlated Gaussian geminals. Variational energies and low‐order perturbation wave functions for the model system are shown to be in excellent agreement with highly accurate numerical solutions for that system. Approximations of the reference wave functions, described by fewer basis functions, are tested for use in the perturbation equations and shown to provide significant computational advantages with tolerable loss of accuracy. Lower bounds for the radius of convergence of the resulting perturbation expansions are estimated. The proposed method is capable of achieving sub‐μHartree accuracy for all systems considered here. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2003     

Combination of perturbative and variational methods for calculating molecular spectra: calculation of the upsilon=3-5 CH stretch overtone spectrum of CHF3

Highly excited states of the CHF3 molecule belonging to the third, fourth, and fifth Fermi polyad are calculated using a combination of the Van Vleck perturbation theory and a variational treatment. The perturbation theory preconditions the Hamiltonian matrix by transforming away all couplings except those between nearly degenerate states. This transformation is implemented so that eigenvalues can be found with significantly smaller matrices than that which would be needed in the original normal mode representation. Even with preconditioning, at the energies as high as 3-5 quanta in the CH stretch, it is not possible to directly diagonalize the Hamiltonian matrix due to the large basis sets required. Iterative methods, particularly the block-Davidson method, are explored for finding the eigenvalues. The methods are compared and the advantages discussed.     

Convergence of the perturbation theory expansion for spin-spin coupling constants

The convergence of the full perturbation theory expansion for the Fermi contact term in the spin-spin coupling constant of hydrogen fluoride has been studied. By examining the contribution of higher virtual states as the basis set is expanded to 18 atomic orbitals, it is shown that at this level of approximation the expansion does not converge. The need to also establish convergence before accepting values calculated from configuration interaction wavefunctions is discussed.The Puerto Rico Nuclear Center is operated by the University of Puerto Rico for the United States Atomic Energy Commission under Contract No. AT(40-1)-1833.     

Consistent generalization of the Møller-Plesset partitioning to open-shell and multiconfigurational SCF reference states in many-body perturbation theory

The Møller-Plesset partitioning of the many-electron Hamiltonian is generalized for arbitrary open-shell and multiconfigurational SCF reference states. The method has been initially implemented at the second and third orders for two-configuration generalized valence bond reference wavefunctions in which only two electrons are correlated. Potential curves for the hydrogen molecule, lithium hydride and the helium dication agree excellently with full CI results.     

Calculation of excited singlet and triplet state polarizabilities using the CNDO/S CI method. An application to naphthalene

A procedure is outlined for the calculation of molecular static electric polarizabilities in excited singlet and triplet states using the ”finite perturbation theory“ in conjunction with the CNDO/S CI method. Numerical results for the ground and the three lowest excited singlet and triplet states of naphthalene are presented. It turns out that the generalized Hellmann-Feynman theorem is approximately valid for the CNDO/S CI wavefunctions and that triplet and singlet state polarizabilities in states of the same symmetry may strongly differ.     

Lithium atom spin density from the Hiller-Sucher-Feinberg identity

Summary The electronic spin density, which determines the observed Fermi contact hyperfine splitting, is usually represented by a delta function operator at the nucleus. Approximate wavefunctions determined by overall energetic considerations may show large errors for such a highly localized property. Hiller, Sucher, and Feinberg (HSF) have shown that the delta function operator can be replaced by a global operator. The possibility that this may lead to an improved method for calculation of the spin density is examined for the ground and first excited states of the lithium atom. Particular attention is given to simple spin polarization wavefunctions that provide the leading contributions to the spin density. It is found that the delta function and HSF formulations give very nearly the same results when the wavefunctions are determined by essentially exact numerical methods. However, the HSF approach shows clear advantages in calculations carried out with finite Slater or contracted Gaussian type basis sets.     

Energy shift in (dtμ)e due to the finite size of the muonic molecular ion (dtμ)+

The energy shift due to the presence of the extended (dtμ)₁₁ pseudonucleus (in its first excited state with one unit of angular momentum) in the quasihydrogenlike system (dtμ)₁₁ e _1s is estimated using perturbation theory up to second order with two choices of zeroth order electron wavefunctions. The energy shift is found to be 0.50 meV using pure Coulomb electron wavefunctions and 0.58 meV using electron wavefunctions calculated with a potential modified to take partial account of the finite size of (dtμ)₁₁. In both cases, the perturbation Hamiltonian is expanded in multipoles, retaining terms up to and including octupole terms.     

Application of miller's classical path approximation to the calculation of the shape of one-dimensional anharmonic probability distributions

Miller's classical-limit expression for the thermally averaged probability distribution of an ensemble of one-dimensional oscillators is applied to a model anharmonic system, the Morse oscillator. It is shown that this method, which requires no knowledge of the stationary state wavefunctions or eigen values, accurately characterizes the shape of the distribution function at very low as well as higher temperatures.     

Characterization of vibrational transition modes by use of normal forms

Summary In this paper we use the Birkhoff-Gustavson perturbation theory to analyze the vibrational modes of two linearly coupled Morse oscillators in the transition region from normal modes to local modes. Our study is based on: truncation of the Hamiltonian written in normal mode coordinates at the 4th order, transformation to normal form and analytical study; construction and use of the approximate integrals of motion of the exact Hamiltonian according to Birkhoff and Gustavson theory. By comparison with a previous analytical study, we demonstrate that perturbation theory, based either on local or normal modes can be used to accurately describe transition modes.     

Internally contracted multiconfiguration-reference configuration interaction calculations for excited states

Summary The calculation of electronically excited states with the internally contracted multiconfiguration-reference configuration interaction (CMRCI) method is discussed. A straightforward method, in which contracted functions for all states are included in the basis, is shown to be very accurate and stable even in cases of narrow avoided crossings. However, the expense strongly increases with the number of states. A new method is proposed, which employs different contracted basis sets for each state, and in which eigensolutions of the Hamiltonian are found using an approximate projection operator technique. The computational effort for this method scales only linearly with the number of states. The two methods are compared for various applications.Dedicated in honor of Prof. Klaus Ruedenberg     

Improved perturbative treatment of electronic energies from a minimal-norm approach to many-body perturbation theory

We propose a zeroth-order Hamiltonian for many-body perturbation theory based on the unitary decomposition of the two-particle reduced Hamiltonian. For the zeroth-order Hamiltonian constrained to be diagonal in the Hartree-Fock basis set, the two-particle reduced perturbation matrix is chosen to have a minimal Frobenius norm. When compared with the M?ller-Plesset partitioning, the method yields more accurate second-order energies.     

The connected‐moments polynomial approach for Hamiltonian eigenvalues calculation and its application to the one‐particle systems

The new connected‐moments polynomial approach (CMP) is developed for evaluation of Hamiltonian eigenvalues. It is based on properties of specially designed polynomial and does not use any basis set and variational procedure. Like all the methods based on hamiltonain moments knowledge, the CMP is conceptually simple but is less tedious and is usually convergent even for very “crude” trial functions. This method is applicable not only to the ground state energy calculation but also to the excited states. The formalism is presented in two modifications: non‐local (integral) and local (integral‐free) ones. An accuracy of both versions is illustrated by numerical examples of Hamiltonian eigenvalues calculations for harmonic and anharmonic oscillators. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

The relationship between matrix elements of position and momentum in vibronic coupling

A general expression is derived that relates the matrix elements of position and momentum allowing for the possibility of a mixed basis. It is shown that when Born-Oppenheimer wavefunctions for two different electronic states are used as basis functions, the use of the usual expression relating matrix elements of position and momentum can lead to results that may not be of the correct order of magnitude.     

New relativistic atomic natural orbital basis sets for lanthanide atoms with applications to the Ce diatom and LuF3

New basis sets of the atomic natural orbital (ANO) type have been developed for the lanthanide atoms La-Lu. The ANOs have been obtained from the average density matrix of the ground and lowest excited states of the atom, the positive ions, and the atom in an electric field. Scalar relativistic effects are included through the use of a Douglas-Kroll-Hess Hamiltonian. Multiconfigurational wave functions have been used with dynamic correlation included using second-order perturbation theory (CASSCF/CASPT2). The basis sets are applied in calculations of ionization energies and some excitation energies. Computed ionization energies have an accuracy better than 0.1 eV in most cases. Two molecular applications are included as illustration: the cerium diatom and the LuF3 molecule. In both cases it is shown that 4f orbitals are not involved in the chemical bond in contrast to an earlier claim for the latter molecule.     

Second order explicitly correlated R12 theory revisited: a second quantization framework for treatment of the operators' partitionings

Second order R12 theory is presented and derived alternatively using the second quantized hole-particle formalism. We have shown that in order to ensure the strong orthogonality between the R12 and the conventional part of the wave function, the explicit use of projection operators can be easily avoided by an appropriate partitioning of the involved operators to parts which are fully describable within the computational orbital basis and complementary parts that involve imaginary orbitals from the complete orbital basis. Various Hamiltonian splittings are discussed and computationally investigated for a set of nine molecules and their atomization energies. If no generalized Brillouin condition is assumed, with all relevant partitionings the one-particle contribution arising in the explicitly correlated part of the first order wave function has to be considered and has a significant role when smaller atomic orbital basis sets are used. The most appropriate Hamiltonian splitting results if one follows the conventional perturbation theory for a general non-Hartree-Fock reference. Then, no couplings between the R12 part and the conventional part arise within the first order wave function. The computationally most favorable splitting when the whole complementary part of the Hamiltonian is treated as a perturbation fails badly. These conclusions also apply to MP2-F12 approaches with different correlation factors.     

Ivar Waller 90 years on June 11, 1988

Linkage properties of the diagrammatic representation of the energies obtained in the multireference many-body perturbation calculations with respect to the incompleteness or completeness of the model space are discussed. The case of not completely degenerate model space is considered for which a comparison with the standard single-reference many-body perturbation expansion is possible. The Hose–Kaldor type of graphical representation of the perturbation expansion for the effective Hamiltonian is used in this comparison. It is shown that for an incomplete model space the perturbation expansion is not size-extensive. In this case, for a truncated expansion of the effective Hamiltonian, the energies obtained by diagonalization of the effective Hamiltonian matrix are represented by both linked and unlinked irreducible contributions. The unlinked ones do not appear when the complete model space is used.     

On the quantum field theory of photoionisation and electron scattering reactions on atoms

Brillouin-Wigner perturbation theory, formulated in the Schrödinger picture of quantum field theory, is employed to derive a perturbative scheme for the scattering matrix for photoionisation and electron scattering reactions on atoms. It is important to note that the intermediate states appearing in this series are physical, i.e. fully correlated, eigenstates of the total Hamiltonian. The scheme is amenable to numerical analysis: The key point is the use of an “energy-optimised”g-Hartree basis which yields an efficient treatment of atomic correlations and makes Brillouin-Wigner and Rayleigh-Schrödinger perturbation theories coincide.