Accurate calculations of intermolecular interaction energies using explicitly correlated coupled cluster wave functions and a dispersion-weighted MP2 method.

Explicitly correlated coupled-cluster calculations of intermolecular interaction energies for the S22 benchmark set of Jurecka, Sponer, Cerny, and Hobza (Chem. Phys. Phys. Chem. 2006, 8, 1985) are presented. Results obtained with the recently proposed CCSD(T)-F12a method and augmented double-zeta basis sets are found to be in very close agreement with basis set extrapolated conventional CCSD(T) results. Furthermore, we propose a dispersion-weighted MP2 (DW-MP2) approximation that combines the good accuracy of MP2 for complexes with predominately electrostatic bonding and SCS-MP2 for dispersion-dominated ones. The MP2-F12 and SCS-MP2-F12 correlation energies are weighted by a switching function that depends on the relative HF and correlation contributions to the interaction energy. For the S22 set, this yields a mean absolute deviation of 0.2 kcal/mol from the CCSD(T)-F12a results. The method, which allows obtaining accurate results at low cost, is also tested for a number of dimers that are not in the training set.

Measures of ionicity of alkaline-earth oxides from the analysis of ab initio cluster wave functions

We present an analysis of the M-O chemical bonding in the binary oxides MgO, CaO, SrO, BaO, and Al2O3 based on ab initio wave functions. The model used to represent the local environment of a metal cation in the bulk oxide is an MO6 cluster which also includes the effect of the lattice Madelung potential. The analysis of the wave functions for these clusters leads to the conclusion that all the alkaline-earth oxides must be regarded as highly ionic oxides; however, the ionic character of the oxides decreases as one goes from MgO, almost perfectly ionic, to BaO. In Al2O3 the ionic character is further reduced; however, even in this case, the departure from the ideal, fully ionic, model of Al3+ is not exceptionally large. These conclusions are based on three measures, a decomposition of the Mq+-Oq- interaction energy, the number of electrons associated to the oxygen ions as obtained from a projection operator technique, and the analysis of the cation core-level binding energies. The increasing covalent character along the series MgO, CaO, SrO, and BaO is discussed in view of the existing theoretical models and experimental data.

Optimization of trial wave functions for use in quantum Monte Carlo with application to LiH

Methods for both partial and full optimization of wavefunction parameters are explored, and these are applied to the LiH molecule. A partial optimization can be easily performed with little difficulty. But to perform a full optimization we must avoid a wrong minimum, and deal with linear-dependency, time step-dependency and ensemble-dependency problems. Five basis sets are examined. The optimized wavefunction with a 3-function set gives a variational energy of -7.998 + 0.005 a.u., which is comparable to that (-7.990 + 0.003) 1 of Reynold's unoptimized \fin ( a double-~ set of eight functions). The optimized wavefunction with a double~ plus 3dz2 set gives ari energy of -8.052 + 0.003 a.u., which is comparable with the fixed-node energy (-8.059 + 0.004)1 of the \fin. The optimized double-~ function itself gives an energy of -8.049 + 0.002 a.u. Each number above was obtained on a Bourrghs 7900 mainframe computer with 14 -15 hrs CPU time.

Wave functions for the methane molecule

If the potential field due to the nuclei in the methane molecule is expanded in terms of a set of spherical harmonics about the carbon nucleus, only the terms involving s, f, and higher harmonic functions differ from zero in the equilibrium configuration. Wave functions have been calculated for the equilibrium configuration, first including only the spherically symmetric s term in the potential, and secondly including both the s and the f terms. In the first calculation the complete Hartree-Fock S.C.F. wave functions were determined; in the second calculation a variation method was used to determine the best form of the wave function involving f harmonics. The resulting wave functions and electron density functions are presented and discussed

HELICITY WAVE-FUNCTIONS FOR MASSLESS AND MASSIVE SPIN-2 PARTICLES

Starting from general properties of a spin-2 field, we construct helicity wave functions in the framework of the Weyl-van der Waerden spinor formalism. We discuss here the cases of massless and massive spin-2 particles.

Comment on Wave functions for a Duffin-Kemmer-Petiau particle in a time-dependent potential [J. Math. Phys. 48, 073515 (2007)]

It is shown that the paper Wave functions for a Duffin-Kemmer-Petiau particle in a time-dependent potential by Merad and Bensaid [J. Math. Phys. 48, 073515 (2007)] is not correct in using inadvertently a non-Hermitian Hamiltonian in a formalism that does require Hermitian Hamiltonians.

Classes of exact wave functions for general time-dependent Dirac Hamiltonians in 1+1 dimensions

The construction of two classes of exact solutions for the most general time-dependent Dirac Hamiltonian in 1+1 dimensions was discussed. The extension of solutions by introduction of a time-dependent mass was elaborated. The possibility of existence of a generalized Lewis-Riesenfeld invariant connected with such solutions was also analyzed.

Compact Singlet S Helium Wave functions(corrected)

This is a corrected version of the Phys. Rev. A 74,14501 (2006) article. The result is improved slightly from that in the original paper.

Cyclobutadiene automerization and rotation of ethylene: Energetics of the barriers by using Spin-polarized wave functions

Spin-projected spin polarized Møller–Plesset and spin polarized coupled clusters calculations have been made to estimate the cyclobutadiene automerization, the ethylene torsion barriers in their ground state, and the gap between the singlet and triplet states of ethylene. The results have been obtained optimizing the geometries at MP4 and/or CCSD levels, by an extensive Gaussian basis set. A comparative analysis with more complex calculations, up to MP5 and CCSDTQP, together with others from the literature, have also been made, showing the efficacy of using spin-polarized wave functions as a reference wave function for Møller–Plesset and coupled clusters calculations, in such problems.

Prolate spheroidal wave functions for electromagnetic theory /

"Contract AF33(616)-6079 Project No. 9-(13-6278) Task 40572. Sponsored by: Wright Air Development Center"

Chebyshev real wave packet propagation: H+O-2 (J=0) state-to-state reactive scattering calculations

In this paper we explore the relative performance of two recently developed wave packet methodologies for reactive scattering, namely the real wave packet Chebyshev domain propagation of Gray and Balint-Kurti [J. Chem. Phys. 108, 950 (1998)] and the Lanczos subspace wave packet approach of Smith [J. Chem. Phys. 116, 2354 (2002); Chem. Phys. Lett. 336, 149 (2001)]. In the former method, a modified Schrodinger equation is employed to propagate the real part of the wave packet via the well-known Chebyshev iteration. While the time-dependent wave packet from the modified Schrodinger equation is different from that obtained using the standard Schrodinger equation, time-to-energy Fourier transformation yields wave functions which differ only trivially by normalization. In the Lanczos subspace approach the linear system of equations defining the action of the Green operator may be solved via either time-dependent or time-independent methods, both of which are extremely efficient due to the simple tridiagonal structure of the Hamiltonian in the Lanczos representation. The two different wave packet methods are applied to three dimensional reactive scattering of H+O-2 (total J=0). State-to-state reaction probabilities, product state distributions, as well as initial-state-resolved cumulative reaction probabilities are examined. (C) 2002 American Institute of Physics.

Bargmann-Wigner method and (6S + 1)-component wave equations

A modified Bargmann-Wigner method is used to derive (6s + 1)-component wave equations. The relation between different forms of these equations is shown.

New predictions for inclusive heavy-quarkonium P-wave decays

We show that some nonrelativistic quantum chromodynamics color-octet matrix elements can be written in terms of (derivatives of) wave functions at the origin and of nonperturbative universal constants once the factorization between the soft and ultrasoft scales is achieved by using an effective field theory where only ultrasoft degrees of freedom are kept as dynamical entities. This allows us to derive a new set of relations between inclusive heavy-quarkonium P-wave decays into light hadrons with different principal quantum numbers and with different heavy flavors. In particular, we can estimate the ratios of the decay widths of bottomonium P-wave states from charmonium data.

New predictions for inclusive heavy-quarkonium p-wave decays

We show that some nonrelativistic quantum chromodynamics color-octet matrix elements can be written in terms of (derivatives of) wave functions at the origin and of nonperturbative universal constants once the factorization between the soft and ultrasoft scales is achieved by using an effective field theory where only ultrasoft degrees of freedom are kept as dynamical entities. This allows us to derive a new set of relations between inclusive heavy-quarkonium P-wave decays into light hadrons with different principal quantum numbers and with different heavy flavors. In particular, we can estimate the ratios of the decay widths of bottomonium P-wave states from charmonium data.