Basis set superposition error-counterpoise corrected potential energy surfaces : application to hydrogen peroxide ... X (X=F-, Cl-, Br-, Li+, Na+) complexes

Møller-Plesset (MP2) and Becke-3-Lee-Yang-Parr (B3LYP) calculations have been used to compare the geometrical parameters, hydrogen-bonding properties, vibrational frequencies and relative energies for several X- and X+ hydrogen peroxide complexes. The geometries and interaction energies were corrected for the basis set superposition error (BSSE) in all the complexes (1-5), using the full counterpoise method, yielding small BSSE values for the 6-311 + G(3df,2p) basis set used. The interaction energies calculated ranged from medium to strong hydrogen-bonding systems (1-3) and strong electrostatic interactions (4 and 5). The molecular interactions have been characterized using the atoms in molecules theory (AIM), and by the analysis of the vibrational frequencies. The minima on the BSSE-counterpoise corrected potential-energy surface (PES) have been determined as described by S. Simón, M. Duran, and J. J. Dannenberg, and the results were compared with the uncorrected PES

A chemical Hamiltonian approach study of the basis set superposition error changes on electron densities and one- and two-center energy components

A comparision of the local effects of the basis set superposition error (BSSE) on the electron densities and energy components of three representative H-bonded complexes was carried out. The electron densities were obtained with Hartee-Fock and density functional theory versions of the chemical Hamiltonian approach (CHA) methodology. It was shown that the effects of the BSSE were common for all complexes studied. The electron density difference maps and the chemical energy component analysis (CECA) analysis confirmed that the local effects of the BSSE were different when diffuse functions were present in the calculations

Effect of basis set superposition error on the electron density of molecular complexes

The effect of basis set superposition error (BSSE) on molecular complexes is analyzed. The BSSE causes artificial delocalizations which modify the first order electron density. The mechanism of this effect is assessed for the hydrogen fluoride dimer with several basis sets. The BSSE-corrected first-order electron density is obtained using the chemical Hamiltonian approach versions of the Roothaan and Kohn-Sham equations. The corrected densities are compared to uncorrected densities based on the charge density critical points. Contour difference maps between BSSE-corrected and uncorrected densities on the molecular plane are also plotted to gain insight into the effects of BSSE correction on the electron density

Second-order Møller-Plesset perturbation theory without basis set superposition error : II : open-shell systems

The basis set superposition error-free second-order MØller-Plesset perturbation theory of intermolecular interactions was studied. The difficulties of the counterpoise (CP) correction in open-shell systems were also discussed. The calculations were performed by a program which was used for testing the new variants of the theory. It was shown that the CP correction for the diabatic surfaces should be preferred to the adiabatic ones

C-H···O H-bonded complexes: how does basis set superposition error change their potential-energy surfaces?

Geometries, vibrational frequencies, and interaction energies of the CNH⋯O3 and HCCH⋯O3 complexes are calculated in a counterpoise-corrected (CP-corrected) potential-energy surface (PES) that corrects for the basis set superposition error (BSSE). Ab initio calculations are performed at the Hartree-Fock (HF) and second-order Møller-Plesset (MP2) levels, using the 6-31G(d,p) and D95++(d,p) basis sets. Interaction energies are presented including corrections for zero-point vibrational energy (ZPVE) and thermal correction to enthalpy at 298 K. The CP-corrected and conventional PES are compared; the unconnected PES obtained using the larger basis set including diffuse functions exhibits a double well shape, whereas use of the 6-31G(d,p) basis set leads to a flat single-well profile. The CP-corrected PES has always a multiple-well shape. In particular, it is shown that the CP-corrected PES using the smaller basis set is qualitatively analogous to that obtained with the larger basis sets, so the CP method becomes useful to correctly describe large systems, where the use of small basis sets may be necessary

How does basis set superposition error change the potential surfaces for hydrogen-bonded dimers?

We describe a simple method to automate the geometric optimization of molecular orbital calculations of supermolecules on potential surfaces that are corrected for basis set superposition error using the counterpoise (CP) method. This method is applied to the H-bonding complexes HF/HCN, HF/H2O, and HCCH/H2O using the 6-31G(d,p) and D95 + + (d,p) basis sets at both the Hartree-Fock and second-order Møller-Plesset levels. We report the interaction energies, geometries, and vibrational frequencies of these complexes on the CP-optimized surfaces; and compare them with similar values calculated using traditional methods, including the (more traditional) single point CP correction. Upon optimization on the CP-corrected surface, the interaction energies become more negative (before vibrational corrections) and the H-bonding stretching vibrations decrease in all cases. The extent of the effects vary from extremely small to quite large depending on the complex and the calculational method. The relative magnitudes of the vibrational corrections cannot be predicted from the H-bond stretching frequencies alone

Basis set superposition error effects, excited-state potential energy surface and photodynamics of thymine

En aquesta tesi he estudiat l'efecte de l'error de superposició de base (BSSE) en la planaritat d'algunes molècules. He observat que l'ús d'alguns mètodes de càlcul amb determinades funcions de base descriuen mínims d'energia no planars per les bases nitrogenades de l'ADN. He demostrat que aquests problemes es poden arreglar utilitzant el mètode Counterpoise per corregir el BSSE en els càlculs. En aquesta tesi també he estudiat la fotofísica de la timina i els resultats mostren que existeixen dos camins de relaxació des de l'estat excitat que permeten la regeneració de l'estructura inicial de forma ultraràpida.

Implementation and application of basis set superposition error-correction schemes to the theoretical modeling of weak intermolecular interactions

This thesis deals with the so-called Basis Set Superposition Error (BSSE) from both a methodological and a practical point of view. The purpose of the present thesis is twofold: (a) to contribute step ahead in the correct characterization of weakly bound complexes and, (b) to shed light the understanding of the actual implications of the basis set extension effects in the ab intio calculations and contribute to the BSSE debate. The existing BSSE-correction procedures are deeply analyzed, compared, validated and, if necessary, improved. A new interpretation of the counterpoise (CP) method is used in order to define counterpoise-corrected descriptions of the molecular complexes. This novel point of view allows for a study of the BSSE-effects not only in the interaction energy but also on the potential energy surface and, in general, in any property derived from the molecular energy and its derivatives A program has been developed for the calculation of CP-corrected geometry optimizations and vibrational frequencies, also using several counterpoise schemes for the case of molecular clusters. The method has also been implemented in Gaussian98 revA10 package. The Chemical Hamiltonian Approach (CHA) methodology has been also implemented at the RHF and UHF levels of theory for an arbitrary number interacting systems using an algorithm based on block-diagonal matrices. Along with the methodological development, the effects of the BSSE on the properties of molecular complexes have been discussed in detail. The CP and CHA methodologies are used for the determination of BSSE-corrected molecular complexes properties related to the Potential Energy Surfaces and molecular wavefunction, respectively. First, the behaviour of both BSSE-correction schemes are systematically compared at different levels of theory and basis sets for a number of hydrogen-bonded complexes. The Complete Basis Set (CBS) limit of both uncorrected and CP-corrected molecular properties like stabilization energies and intermolecular distances has also been determined, showing the capital importance of the BSSE correction. Several controversial topics of the BSSE correction are addressed as well. The application of the counterpoise method is applied to internal rotational barriers. The importance of the nuclear relaxation term is also pointed out. The viability of the CP method for dealing with charged complexes and the BSSE effects on the double-well PES blue-shifted hydrogen bonds is also studied in detail. In the case of the molecular clusters the effect of high-order BSSE effects introduced with the hierarchical counterpoise scheme is also determined. The effect of the BSSE on the electron density-related properties is also addressed. The first-order electron density obtained with the CHA/F and CHA/DFT methodologies was used to assess, both graphically and numerically, the redistribution of the charge density upon BSSE-correction. Several tools like the Atoms in Molecules topologycal analysis, density difference maps, Quantum Molecular Similarity, and Chemical Energy Component Analysis were used to deeply analyze, for the first time, the BSSE effects on the electron density of several hydrogen bonded complexes of increasing size. The indirect effect of the BSSE on intermolecular perturbation theory results is also pointed out It is shown that for a BSSE-free SAPT study of hydrogen fluoride clusters, the use of a counterpoise-corrected PES is essential in order to determine the proper molecular geometry to perform the SAPT analysis.

Interaccions moleculars: correcció de l'error de superposició de base i propietat de les funcions flotants

A primary interest of this thesis is to obtain a powerful tool for determining structural properties, electrical and reactivity of molecules. A second interest is the study of fundamental error based on complex overlay of bridge hydrogen. One way to correct this error, using Counterpoise correction proposed by Boys and Bernardi. Usually the Counterpoise correction is applied promptly on the geometries previously optimized. Our goal was to find areas of potential which had all the points fixed with CP. These surfaces have a minimum corresponding to a surface other than corrected, ie, the geometric parameters will be different. The curvature of this minimum will also be different, therefore the vibrational frequency will also change when they are corrected with BSSE. Once constructed these surfaces have been studied various complex. It has also been investigated as the method for calculating the error influenced on the basis superposition.

Ab initio studies of the structure and vibrational spectra of methyl isocyanide-borane and acetonitrile-borane

The SCF/DZP and MP2/DZP methods of ab initio quantum chemistry have been utilized to study the structure, vibrational spectra, binding energy, and barrier to internal rotation of methyl isocyanide-borane and acetonitrile-borane adducts. The eclipsed conformation of the complexes was predicted to be a minimum, and the staggered form is a transition state with a barrier height of about 10 cal/mol. The vibrational analyses of CH3NC-BH3 and CH3CN-BH3 and several of their isotopomers have been carried out by the GF matrix method. Computations have also been carried out for free CH3NC and CH3CN in order to investigate the changes in CH3NC and CH3CN as a result of their complex formation with BH3. To obtain an acceptable set of force constants, a recently proposed procedure ''RECOVES'' has been utilized. The increase in the N=C/C=N stretching force constant of CH3NC/CH3CN on adduct formation is interpreted with the help of Parr and Borkman's model. The binding energies for the two adducts have been determined taking basis set superposition error (BSSE) into consideration. The effect of the BSSE on structure, dipole moment, and vibrational frequencies of CH3CN and CH3NC is also evaluated. The predicted infrared band intensities for the two complexes are in good agreement with the experimentally observed features, and they have been utilized in the assignment of vibrational frequencies.

Quantum chemical study on influence of intermolecular hydrogen bonding on the geometry, the atomic charges and the vibrational dynamics of 2,6-dichlorobenzonitrile

FT-IR (4000-400 cm(-1)) and FT-Raman (4000-200 cm(-1)) spectral measurements on solid 2,6-dichlorobenzonitrile (2,6-DCBN) have been done. The molecular geometry, harmonic vibrational frequencies and bonding features in the ground state have been calculated by density functional theory at the B3LYP/6-311++G (d,p) level. A comparison between the calculated and the experimental results covering the molecular structure has been made. The assignments of the fundamental vibrational modes have been done on the basis of the potential energy distribution (PED). To investigate the influence of intermolecular hydrogen bonding on the geometry, the charge distribution and the vibrational spectrum of 2,6-DCBN; calculations have been done for the monomer as well as the tetramer. The intermolecular interaction energies corrected for basis set superposition error (BSSE) have been calculated using counterpoise method. Based on these results, the correlations between the vibrational modes and the structure of the tetramer have been discussed. Molecular electrostatic potential (MEP) contour map has been plotted in order to predict how different geometries could interact. The Natural Bond Orbital (NBO) analysis has been done for the chemical interpretation of hyperconjugative interactions and electron density transfer between occupied (bonding or lone pair) orbitals to unoccupied (antibonding or Rydberg) orbitals. UV spectrum was measured in methanol solution. The energies and oscillator strengths were calculated by Time Dependent Density Functional Theory (TD-DFT) and matched to the experimental findings. TD-DFT method has also been used for theoretically studying the hydrogen bonding dynamics by monitoring the spectral shifts of some characteristic vibrational modes involved in the formation of hydrogen bonds in the ground and the first excited state. The C-13 nuclear magnetic resonance (NMR) chemical shifts of the molecule were calculated by the Gauge independent atomic orbital (GIAO) method and compared with experimental results. Standard thermodynamic functions have been obtained and changes in thermodynamic properties on going from monomer to tetramer have been presented. (C) 2013 Elsevier B.V. All rights reserved.

Theoretical investigation on the adsorption of Ag+ and hydrated Ag+ cations on clean Si(111) surface

In this paper, the adsorption of Ag+ and hydrated Ag+ cations on clean Si(111) surface were investigated by using cluster (Gaussian 03) and periodic (DMol(3)) ab initio calculations. Si(111) surface was described with cluster models (Si14H17 and Si22H21) and a four-silicon layer slab with periodic boundary conditions. The effect of basis set superposition error (BSSE) was taken into account by applying the counterpoise correction. The calculated results indicated that the binding energies between hydrated Ag+ cations and clean Si(111) surface are large, suggesting a strong interaction between hydrated Ag+ cations and the semiconductor surface. With the increase of number, water molecules form hydrogen bond network with one another and only one water molecule binds directly to the Ag+ cation. The Ag+ cation in aqueous solution will safely attach to the clean Si(111) surface.

Theoretical study on adsorption of Na+ and Na+(H2O)(n) (n=1-6) on a clean Si(111) surface

To model the adsorption of Na+ in aqueous solution on the semiconductor surface, the interactions of Na+ and Na+(H2O)(n) (n = 1-6) with a clean Si(111) surface were investigated by using hybrid density functional theory (B3LYP) and Moller-Plesset second-order perturbation (MP2) methods. The Si(111) surface was described with Si8H12, Si16H20, and Si22H21 Cluster models. The effect of the basis set superposition error (BSSE) was taken into account by applying the counterpoise (CP) correction. The calculated results indicated that the interactions between the Na+ cation and the dangling bonds of the Si(111) surface are primarily electrostatic with partial orbital interactions. The magnitude of the binding energies depends weakly on the adsorption sites and the size of the clusters. When water molecules are present, the interaction between the Nal and Si(I 11) surfaces weakens and the binding energy has the tendency to saturate. On a Si22H21 cluster described surface, the optimized Na+-surface distance for Na+(H2O)(5) adsorbed at on-top site is 4.16 angstrom and the CP-corrected binding energy (MP2) is -35.4 kJ/mol, implying a weakly adsorption of hydrated Na+ cation on clean Si(111) surface.

12C2H2-Ar van der Waals complex

New theoretical and experimental results on the acetylene-Ar van der Waals complex are presented and the literature is reviewed. New ab initio calculations at the MP2 level were performed using large basis sets with diffuse functions and taking into account the basis set superposition error. It was found that the structure of acetylene is not significantly altered by the complexation and that its vibrational frequencies are only slightly lowered. Finally, it was observed that the calculated properties of the complex (structure, vibrational spectrum, bond dissociation energy) are not sensitive to the structure imposed on acetylene. Experimentally, acetylene-Ar was produced in a supersonic expansion under experimental conditions corresponding to 9 K rotational temperature. Thanks to the performances of CW-CRDS detection, the Ka = 0 ← 1, 1 ← 0, and 2 ← 1 sub- bands of the v1 + v3 band could be recorded and resolved and most of their lines assigned. Upper-state rotational constants were fitted, however not including the upper Ka = 2 state, which shows K-doubling the opposite of the expected. The Lorentzian width of most line profiles sets the mean lifetime to some 7.5 ns. Local perturbations affecting line positions and/or line widths are demonstrated. Additional series of lines tentatively attributed to acetylene-Ar are discussed.© 2009 American Chemical Society.

Accurate ground-state potential energy surfaces of the C2H2–Kr and C2H2–Xe van der Waals complexes

Accurate ab initio intermolecular potential energy surfaces (IPES) have been obtained for the first time for the ground electronic state of the C 2H2-Kr and C2H2-Xe van der Waals complexes. Extensive tests, including complete basis set and all-electron scalar relativistic results, support their calculation at the CCSD(T) level of theory, using small-core relativistic pseudopotentials for the rare-gas atoms and aug-cc-pVQZ basis sets extended with a set of 3s3p2d1f1g mid-bond functions. All results are corrected for the basis set superposition error. The importance of the scalar relativistic and rare-gas outer-core (n.1)d correlation effects is investigated. The calculated IPES, adjusted to analytical functions, are characterized by global minima corresponding to skew T-shaped geometries, in which the Jacobi vector positioning the rare-gas atom with respect to the center of mass of the C2H2 moiety corresponds to distances of 4.064 and 4.229Å, and angles of 65.22° and 68.67° for C 2H2-Kr and C2H2-Xe, respectively. The interaction energy of both complexes is estimated to be -151.88 (1.817 kJ mol-1) and -182.76 cm-1 (2.186 kJ mol-1), respectively. The evolution of the topology of the IPES as a function of the rare-gas atom, from He to Xe, is also discussed. © 2012 Taylor and Francis.