Eliminating spin contamination in auxiliary-field quantum Monte Carlo: realistic potential energy curve of F2

The use of an approximate reference state wave function |Phi_r> in electronic many-body methods can break the spin symmetry of Born-Oppenheimer spin-independent Hamiltonians. This can result in significant errors, especially when bonds are stretched or broken. A simple spin-projection method is introduced for auxiliary-field quantum Monte Carlo (AFQMC) calculations, which yields spin-contamination-free results, even with a spin-contaminated |Phi_r>. The method is applied to the difficult F2 molecule, which is unbound within unrestricted Hartree-Fock (UHF). With a UHF |Phi_r>, spin contamination causes large systematic errors and long equilibration times in AFQMC in the intermediate, bond-breaking region. The spin-projection method eliminates these problems, and delivers an accurate potential energy curve from equilibrium to the dissociation limit using the UHF |Phi_r>. Realistic potential energy curves are obtained with a cc-pVQZ basis. The calculated spectroscopic constants are in excellent agreement with experiment.Comment: 8 pages, 6 figures, submitted to J. Chem. Phy

Three-body non-additive forces between spin-polarized alkali atoms

Three-body non-additive forces in systems of three spin-polarized alkali atoms (Li, Na, K, Rb and Cs) are investigated using high-level ab initio calculations. The non-additive forces are found to be large, especially near the equilateral equilibrium geometries. For Li, they increase the three-atom potential well depth by a factor of 4 and reduce the equilibrium interatomic distance by 0.9 A. The non-additive forces originate principally from chemical bonding arising from sp mixing effects.Comment: 4 pages, 3 figures (in 5 files

Comment on “The diatomic dication CuZn2+ in the gas phase” [J. Chem. Phys. 135, 034306 (2011)]

In this Comment, the density functional theory (DFT) calculations carried out by Diez et al. [J. Chem. Phys. 135, 034306 (2011)10.1063/1.3613624] are revised within the framework of the coupled-cluster single double triple method. These more sophisticated calculations allow us to show that the 2Σ+ electronic ground state of CuZn2+, characterized as the metastable ground state by DFT calculations, is a repulsive state instead. The 2Δ and 2Π metastable states of CuZn2+, on the other hand, should be responsible for the formation mechanism of the dication through the near-resonant electron transfer CuZn+ + Ar+ → CuZn2+ + Ar reaction.Fil: Pis Diez, Reinaldo. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - La Plata. Centro de Química Inorgánica "Dr. Pedro J. Aymonino". Universidad Nacional de La Plata. Facultad de Ciencias Exactas. Centro de Química Inorgánica "Dr. Pedro J. Aymonino"; ArgentinaFil: Franzreb, Klaus. Arizona State University; Estados UnidosFil: Alonso, Julio A.. Universidad de Valladolid; Españ

Core repulsion effects in alkali trimers

The present paper is related to a talk presented during the Symposium on Coherent Control and Ultracold Chemistry held during the Sixth Congress of the International Society for Theoretical Chemical Physics (ISTCP-VI, July 2008). The talk was entitled "Electronic structure properties of alkali dimers and trimers. Prospects for alignment of ultracold molecules". Here we report on the electrostatic repulsion forces of the ionic cores at short separation, involved when the potential energy surfaces of alkali trimers are calculated with a quantum chemistry approach based on effective large-core potentials for ionic core description. We demonstrate that such forces in the triatomic molecule can be obtained as the sum of three pairwise terms. We illustrate our results on the lowest electronic states of Cs$_3$, which are computed for the first time within a full configuration interaction based on a large Gaussian basis set. As a preliminary section, we also propose a brief introduction about the importance of alkali trimer systems in the context of cold and ultracold molecules

Mechanisms for the Oxonolysis of Ethene and Propene: Reliability of Quantum Chemical Predictions

Reactions of ozone with ethene and propene leading to primary ozonide (concerted and stepwise ozonolysis) or epoxide and singlet molecular oxygen (partial ozonolysis) are studied theoretically. The mechanism of concerted ozonolysis proceeds via a single transition structure which is a partial diradical. The transition structures and intermediates in the stepwise ozonolysis and partial ozonolysis mechanisms are singlet diradicals. Spin-restricted and unrestricted density functional methods are employed to calculate the structures of the closed-shell and diradical species. Although the partial diradicals exhibit moderate to pronounced instability in their RDFT and RHF solutions, RDFT is required to locate the transition structure for concerted ozonolysis. Spin projected fourth-order Møller–Plesset theory (PMP4) was used to correct the DFT energies. The calculated pre-exponential factors and activation energies for the concerted ozonolysis of ethene and propene are in good agreement with experimental values. However, the PMP4//DFT procedure incorrectly predicts the stepwise mechanism as the favored channel. UCCSD(T) predicts the concerted mechanism as the favored channel but significantly overestimates the activation energies. RCCSD(T) is found to be more accurate than UCCSD(T) for the calculation of the concerted mechanism but is not applicable to the diradical intermediates. The major difficulty in accurate prediction of the rate constant data for these reactions is the wide range of spin contamination for the reference UHF wave functions and UDFT solutions across the potential energy surface. The possibility of the partial ozonolysis mechanism being the source of epoxide observed in some experiments is discussed

Relativistic coupled cluster calculations of spectroscopic and chemical properties for element 120

The coupled cluster calculations with accounting for relativistic effects to study spectroscopic and chemical properties of element 120 (E120) are performed. Similar calculations for Ba are also done and they are in a good agreement with the experimental data. Dissociation energies of diatomic X-H and X-Au molecules, where X=E120, Ba, are calculated; for E120 they are found to be $1.5\div2$ times smaller than those for Ba

Bond breaking with auxiliary-field quantum Monte Carlo

Bond stretching mimics different levels of electron correlation and provides a challenging testbed for approximate many-body computational methods. Using the recently developed phaseless auxiliary-field quantum Monte Carlo (AF QMC) method, we examine bond stretching in the well-studied molecules BH and N$_2$, and in the H$_{50}$ chain. To control the sign/phase problem, the phaseless AF QMC method constrains the paths in the auxiliary-field path integrals with an approximate phase condition that depends on a trial wave function. With single Slater determinants from unrestricted Hartree-Fock (UHF) as trial wave function, the phaseless AF QMC method generally gives better overall accuracy and a more uniform behavior than the coupled cluster CCSD(T) method in mapping the potential-energy curve. In both BH and N$_2$, we also study the use of multiple-determinant trial wave functions from multi-configuration self-consistent-field (MCSCF) calculations. The increase in computational cost versus the gain in statistical and systematic accuracy are examined. With such trial wave functions, excellent results are obtained across the entire region between equilibrium and the dissociation limit.Comment: 8 pages, 3 figures and 3 tables. Submitted to JC