Basis set superposition error

Last updated

In quantum chemistry, calculations using finite basis sets are susceptible to basis set superposition error (BSSE). As the atoms of interacting molecules (or of different parts of the same molecule - intramolecular BSSE) [1] [2] approach one another, their basis functions overlap. Each monomer "borrows" functions from other nearby components, effectively increasing its basis set and improving the calculation of derived properties such as energy. [3] If the total energy is minimised as a function of the system geometry, the short-range energies from the mixed basis sets must be compared with the long-range energies from the unmixed sets, and this mismatch introduces an error.

Other than using infinite basis sets, two methods exist to eliminate the BSSE. In the chemical Hamiltonian approach (CHA), [4] [5] basis set mixing is prevented a priori, by replacing the conventional Hamiltonian with one in which all the projector-containing terms that would allow mixing have been removed. In the counterpoise method (CP), [6] [7] [8] the BSSE is calculated by re-performing all the calculations using the mixed basis sets, and the error is then subtracted a posteriori from the uncorrected energy. (The mixed basis sets are realised by introducing "ghost orbitals", basis set functions which have no electrons or protons. [2] It however has been shown that there is an inherent danger in using counterpoise corrected energy surfaces, due to the inconsistent effect of the correction in different areas of the energy surface. [3] ) Though conceptually very different, the two methods tend to give similar results. [9] It also has been shown that the error is often larger when using the CP method since the central atoms in the system have much greater freedom to mix with all of the available functions compared to the outer atoms. Whereas in the CHA model, those orbitals have no greater intrinsic freedom and therefore the correction treats all fragments equally. [10] The errors inherent in either BSSE correction disappear more rapidly than the total value of BSSE in larger basis sets. [11]

Related Research Articles

Computational chemistry is a branch of chemistry that uses computer simulation to assist in solving chemical problems. It uses methods of theoretical chemistry, incorporated into efficient computer programs, to calculate the structures and properties of molecules and solids. It is necessary because, apart from relatively recent results concerning the hydrogen molecular ion, the quantum many-body problem cannot be solved analytically, much less in closed form. While computational results normally complement the information obtained by chemical experiments, it can in some cases predict hitherto unobserved chemical phenomena. It is widely used in the design of new drugs and materials.

Density-functional theory (DFT) is a computational quantum mechanical modelling method used in physics, chemistry and materials science to investigate the electronic structure of many-body systems, in particular atoms, molecules, and the condensed phases. Using this theory, the properties of a many-electron system can be determined by using functionals, i.e. functions of another function. In the case of DFT, these are functionals of the spatially dependent electron density. DFT is among the most popular and versatile methods available in condensed-matter physics, computational physics, and computational chemistry.

Chemistry at Harvard Macromolecular Mechanics (CHARMM) is the name of a widely used set of force fields for molecular dynamics, and the name for the molecular dynamics simulation and analysis computer software package associated with them. The CHARMM Development Project involves a worldwide network of developers working with Martin Karplus and his group at Harvard to develop and maintain the CHARMM program. Licenses for this software are available, for a fee, to people and groups working in academia.

Coupled cluster (CC) is a numerical technique used for describing many-body systems. Its most common use is as one of several post-Hartree–Fock ab initio quantum chemistry methods in the field of computational chemistry, but it is also used in nuclear physics. Coupled cluster essentially takes the basic Hartree–Fock molecular orbital method and constructs multi-electron wavefunctions using the exponential cluster operator to account for electron correlation. Some of the most accurate calculations for small to medium-sized molecules use this method.

MOLPRO

MOLPRO is a software package used for accurate ab initio quantum chemistry calculations. It is developed by Peter Knowles at Cardiff University and Hans-Joachim Werner at Universität Stuttgart in collaboration with other authors.

The extended Hückel method is a semiempirical quantum chemistry method, developed by Roald Hoffmann since 1963. It is based on the Hückel method but, while the original Hückel method only considers pi orbitals, the extended method also includes the sigma orbitals.

Møller–Plesset perturbation theory (MP) is one of several quantum chemistry post–Hartree–Fock ab initio methods in the field of computational chemistry. It improves on the Hartree–Fock method by adding electron correlation effects by means of Rayleigh–Schrödinger perturbation theory (RS-PT), usually to second (MP2), third (MP3) or fourth (MP4) order. Its main idea was published as early as 1934 by Christian Møller and Milton S. Plesset.

Per-Olov Löwdin was a Swedish physicist, professor at the University of Uppsala from 1960 to 1983, and in parallel at the University of Florida until 1993.

In computational chemistry, post–Hartree–Fock methods are the set of methods developed to improve on the Hartree–Fock (HF), or self-consistent field (SCF) method. They add electron correlation which is a more accurate way of including the repulsions between electrons than in the Hartree–Fock method where repulsions are only averaged.

A basis set in theoretical and computational chemistry is a set of functions that is used to represent the electronic wave function in the Hartree–Fock method or density-functional theory in order to turn the partial differential equations of the model into algebraic equations suitable for efficient implementation on a computer.

Spartan (chemistry software)

Spartan is a molecular modelling and computational chemistry application from Wavefunction. It contains code for molecular mechanics, semi-empirical methods, ab initio models, density functional models, post-Hartree–Fock models, and thermochemical recipes including G3(MP2) and T1. Quantum chemistry calculations in Spartan are powered by Q-Chem.

Ab initio quantum chemistry methods are computational chemistry methods based on quantum chemistry. The term ab initio was first used in quantum chemistry by Robert Parr and coworkers, including David Craig in a semiempirical study on the excited states of benzene. The background is described by Parr. Ab initio means "from first principles" or "from the beginning", implying that the only inputs into an ab initio calculation are physical constants. Ab initio quantum chemistry methods attempt to solve the electronic Schrödinger equation given the positions of the nuclei and the number of electrons in order to yield useful information such as electron densities, energies and other properties of the system. The ability to run these calculations has enabled theoretical chemists to solve a range of problems and their importance is highlighted by the awarding of the Nobel prize to John Pople and Walter Kohn.

Quantum chemistry composite methods are computational chemistry methods that aim for high accuracy by combining the results of several calculations. They combine methods with a high level of theory and a small basis set with methods that employ lower levels of theory with larger basis sets. They are commonly used to calculate thermodynamic quantities such as enthalpies of formation, atomization energies, ionization energies and electron affinities. They aim for chemical accuracy which is usually defined as within 1 kcal/mol of the experimental value. The first systematic model chemistry of this type with broad applicability was called Gaussian-1 (G1) introduced by John Pople. This was quickly replaced by the Gaussian-2 (G2) which has been used extensively. The Gaussian-3 (G3) was introduced later.

Minnesota Functionals (Myz) are a group of highly parameterized approximate exchange-correlation energy functionals in density functional theory (DFT). They are developed by the group of Prof. Donald Truhlar at the University of Minnesota. These functionals are based on the meta-GGA approximation, i.e. they include terms that depend on the kinetic energy density, and are all based on complicated functional forms parametrized on high-quality benchmark databases. These functionals can be used for traditional quantum chemistry and solid-state physics calculations. The Myz functionals are widely used and tested in the quantum chemistry community.

BigDFT is a free software package for physicists and chemists, distributed under the GNU General Public License, whose main program allows the total energy, charge density, and electronic structure of systems made of electrons and nuclei to be calculated within density functional theory (DFT), using pseudopotentials, and a wavelet basis.

DMol3 is a commercial software package which uses density functional theory with a numerical radial function basis set to calculate the electronic properties of molecules, clusters, surfaces and crystalline solid materials from first principles. DMol3 can either use gas phase boundary conditions or 3d periodic boundary conditions for solids or simulations of lower-dimensional periodicity. It has also pioneered the use of the conductor-like screening model COSMO Solvation Model for quantum simulations of solvated molecules and recently of wetted surfaces. DMol3 permits geometry optimisation and saddle point search with and without geometry constraints, as well as calculation of a variety of derived properties of the electronic configuration. DMol3 development started in the early eighties with B. Delley then associated with A.J.Freeman and D.E.Ellis at Northwestern University. In 1989 DMol3 appeared as DMol, the first commercial density functional package for industrial use by Biosym Technologies now Accelrys. Delley's 1990 publication was cited more than 3000 times.

Poul Jørgensen is professor of chemistry at the Department of Chemistry, Aarhus University (AU), Denmark and director of the qLEAP Center for Theoretical Chemistry at AU, which was established in April 2012. Jørgensen has made seminal contributions to the field of electronic structure theory. He is also one of the main authors of the DALTON program and a member of the International Academy of Quantum Molecular Science.

Pople diagram

A Pople diagram or Pople's Diagram is a diagram which describes the relationship between various calculation methods in computational chemistry. It was initially introduced in January 1965 by Sir John Pople,, during the Symposium of Atomic and Molecular Quantum Theory in Florida. The Pople Diagram can be either 2-dimensional or 3-dimensional, with the axes representing ab inito methods, basis sets and treatment of relativity. The diagram attempts to balance calculations by giving all aspects of a computation equal weight.

Symmetry-adapted perturbation theory or SAPT is a methodology in electronic structure theory developed to describe non-covalent interactions between atoms and/or molecules. SAPT is a member of the family of methods known as energy decomposition analysis (EDA). Most EDA methods decompose a total interaction energy that is computed via a supermolecular approach, such that:

FHI-aims Molecular dynamics modelling software

FHI-aims is a shared-source software package for computational molecular and materials science written in Fortran. It uses density functional theory and many-body perturbation theory to simulate chemical and physical properties of atoms, molecules, nanostructures, solids, and surfaces. Originally developed at the Fritz Haber Institute in Berlin the ongoing development of the FHI-aims source code is now driven by a world-wide community of collaborating research institutions.

References

  1. Balabin, Roman M. (2008). "Enthalpy difference between conformations of normal alkanes: Intramolecular basis set superposition error (BSSE) in the case of n-butane and n-hexane". J. Chem. Phys. 129 (16): 164101. Bibcode:2008JChPh.129p4101B. doi:10.1063/1.2997349. PMID   19045241.CS1 maint: discouraged parameter (link)
  2. 1 2 Hobza, Pavel; Müller-Dethlefs, Klaus (2010). Non-covalent Interactions: Theory and Experiment (PDF). Cambridge, England: Royal Society of Chemistry. p. 13. ISBN   978-1-84755-853-4. Archived from the original (PDF) on 2011-06-06. Retrieved 2010-05-21.
  3. 1 2 Liedl, Klaus R. (1998). "Dangers of counterpoise corrected hypersurfaces. Advantages of basis set superposition improvement". The Journal of Chemical Physics. 108 (8): 3199–3204. Bibcode:1998JChPh.108.3199L. doi:10.1063/1.475715.
  4. Mayer, I.; Valiron, P. (1998). "Second order Møller–Plesset perturbation theory without basis set superposition error". J. Chem. Phys. 109 (9): 3360–3373. Bibcode:1998JChPh.109.3360M. doi:10.1063/1.476931.
  5. Bende, Attila. "THE CHEMICAL HAMILTONIAN APPROACH (CHA)" . Retrieved 14 May 2010.CS1 maint: discouraged parameter (link)
  6. Van Duijneveldt, Frans B.; van Duijneveldt-van de Rijdt, Jeanne G. C. M.; van Lenthe, Joop H. (1994). "State of the art in counterpoise theory". Chem. Rev. 94 (7): 1873–1885. doi:10.1021/cr00031a007.
  7. Rösch, N. (2003). "Counterpoise Correction". Technical University of Munich, Quantum Chemistry Laboratory. Archived from the original on 18 April 2015. Retrieved 14 May 2010.CS1 maint: discouraged parameter (link)
  8. Sedano, Pedro Salvador (2000). "Counterpoise Corrected Potential Energy Surfaces". University of Girona. Retrieved 14 May 2010.CS1 maint: discouraged parameter (link)
  9. Paizs, Béla; Suhai, Sándor (1998). "Comparative study of BSSE correction methods at DFT and MP2 levels of theory". J. Comput. Chem. 19 (6): 575–584. doi:10.1002/(SICI)1096-987X(19980430)19:6<575::AID-JCC1>3.0.CO;2-O.
  10. Mentel, Lukasz; Baerends, Evert Jan (2013). "Can the Counterpoise Correction for Basis Set Superposition Effect Be Justified?". J. Comput. Chem. 10 (1): 252–267. doi:10.1021/ct400990u. PMID   26579908.
  11. Mayer, I. (2004). "Interrelations between the a priori and a posteriori BSSE correction schemes". Int. J. Quantum Chem. 100 (4): 559–566. doi:10.1002/qua.10827.

See also