Designing Gaussian basis sets to the theoretical study of the piezoelectric effect of perovskite (BaTiO3)

Gaussian basis sets (24s14p, 30s19p14d, and 33s21p14d for O (P-3), Ti (S-5), and Ba (S-1) atoms, respectively), are designed with the strategy of the Generator Coordinate Hartree-Fock method. The basis sets are then contracted to [6s4p], [10s5p4d], and [16s9p5d] to O, Ti, and Ba atoms, respectively, and used in calculations of total and orbital energies of (TiO+2)-Ti-1 and (BaO)-Ba-1 fragments for quality evaluation in molecular studies. For O atom, the [6s4p] basis set is enriched with d polarization function and used along with the [10s5p4d] and [16s9p5d] basis sets for the theoretical study of the piezoelectric effect of perovskite (BaTiO3). The results of this work evidence that the piezoelectric properties in BaTiO3 can be caused by electrostatic interactions. (C) 2003 Elsevier B.V. All rights reserved.

Development of basis sets to calculations of the electronic structure of YMnO3

The generator coordinate Hartree-Fock method was used to develop 20s17p, 30s20p14d, and 30s21p16d Gaussian basis sets for the O ((3)p), Mn (S-6), and Y (D-2) atoms, respectively. The Gaussian basis sets were contracted to 20s17p/9s7p, 30s20p14d/11s7p7d, and 30s21p16d/14s7p7d and utilized in calculations of total energy and orbital energies of the (MnO1+)-Mn-5 and (YO1+)-Y-3 fragments to evaluate its quality in molecular studies. Finally, the contracted basis set for O atom was supplemented with one polarization function of d symmetry and used along with the other contracted basis sets (for Mn and Y) to calculate dipole moments, total energy, and total atomic charges in YMnO3 in space group D-6h. The analysis of those properties showed that is reasonable to believe that YMnO3 present behavior of piezoelectric material. (C) 2003 Elsevier B.V. All rights reserved.

Gaussian basis sets by generator coordinate Hartree-Fock method to ab initio calculations of electron affinities of enolates

The Generator Coordinate Hartree-Fock (GCHF) method is employed to generate uncontracted 15s and 18s11p gaussian basis sets for the H, C and O atoms, respectively. These basis sets are then contracted to 3s and 4s H atom and 6s5p, for C and O atoms by a standard procedure. For quality evaluation of contracted basis sets in molecular calculations, we have accomplished calculations of total and orbital energies in the Hartree-Fock-Roothaaii (HFR) approach for CH, C(2) and CO molecules. The results obtained with the uncontracted basis sets are compared with values obtained with the standard D95, 6-311G basis sets and with values reported in the literature. The 4s and 6s5p basis sets are enriched with polarization and diffuse functions for atoms of the parent neutral systems and of the enolates anions (cycloheptanone enolate, 2,5-dimethyleyelopentanone enolate, 4-heptanone enolate, and di-isopropyl ketone enolate) from the literature, in order to assess their performance in ab initio molecular calculations, and applied for calculations of electron affinities of the enolates. The calculations were performed at the DFT (BLYP and B3LYP) and HF levels and compared with the corresponding experimental values and with those obtained by using other 6-3 1 + +G((*)) and 6-311 + +G((*)) basis sets from literature. For the enolates studied, the differences between the electron affinities obtained with GCHF basis sets, at the B3LYP level, and the experimental values are -0.001, -0,014, -0.001, and -0.001 eV. (C) 2002 Elsevier B.V. B.V. All rights reserved.

Contracted GTF basis sets applied to the theoretical interpretation of the Raman spectrum of hexaaquachromium(III) ion

Contracted GTF basis sets designed with aid of the Generator Coordinate Hartree-Fock (GCHF) method for H(2S), O2-(1S), and Cr3+(4F) atomic species are applied to perform theoretical interpretation of the Raman spectrum of hexaaquachromium(III) ion. The 16s, 16s 10p, and 24s17p13d GTF basis sets were contracted to [4s] for H atom, [6s4p], and [9s6p3d] for O2- and Cr3+, respectively, by Dunning's scheme. For Cr3+, the [9s6p3d] basis set was enriched with f polarization function and used in combination com [4s] and [6s4p] in the study of our interest. The results obtained in this report show that the contracted GTF basis sets used are a useful alternative for the theoretical interpretation of Raman spectrum of hexaaquachromium(III) ion and that GCHF method is an effective alternative to selection of GTF basis sets for theoretical study of vibrational properties of poliatomic species. © 2003 Elsevier Science B.V. All rights reserved.

Investigating the Nature of Noble Gas-Copper Bonds by the Quantum Theory of Atoms in Molecules

We investigated noble gas copper bonds in linear complexes represented by the NgCuX general formula in which Ng and X stand for a noble gas (neon, argon, krypton, or xenon) and a halogen (fluorine, chlorine or bromine), respectively, by coupled cluster methods and modified cc-pVQZ basis sets. The quantum theory of atoms in molecules (QTAIM) shows a linear relation between the dissociation energy or noble gas-copper bonds and the amount of electronic charge transferred mainly from the noble gas to copper during complexation. Large changes in the QTAIM quadrupole moments of copper and noble gases resulting from this bonding and a comparison between NgCuX and NgNaCl systems indicate that these noble gas-copper bonds should be better interpreted as predominantly covalent. Finally, QTAIM atomic dipoles of noble gases in NgNaCl systems agree satisfactorily with atomic dipoles given by a simple model for these NgNa van der Waals bonds.

Intramolecular hydrogen bond: Can it be part of the basis set of valence internal coordinates in normal mode analysis?

It has been shown earlier1] that the relaxed force constants (RFCs) could be used as a measure of bond strength only when the bonds form a part of the complete valence internal coordinates (VIC) basis. However, if the bond is not a part of the complete VIC basis, its RFC is not necessarily a measure of bond strength. Sometimes, it is possible to have a complete VIC basis that does not contain the intramolecular hydrogen bond (IMHB) as part of the basis. This means the RFC of IMHB is not necessarily a measure of bond strength. However, we know that IMHB is a weak bond and hence its RFC has to be a measure of bond strength. We resolve this problem of IMHB not being part of the complete basis by postulating `equivalent' basis sets where IMHB is part of the basis at least in one of the equivalent sets of VIC. As long as a given IMHB appears in one of the equivalent complete VIC basis sets, its RFC could be used as a measure of bond strength parameter.

Systematic study of the static electrical properties of the CO molecule: influence of the basis set size and correlation energy

The influence of the basis set size and the correlation energy in the static electrical properties of the CO molecule is assessed. In particular, we have studied both the nuclear relaxation and the vibrational contributions to the static molecular electrical properties, the vibrational Stark effect (VSE) and the vibrational intensity effect (VIE). From a mathematical point of view, when a static and uniform electric field is applied to a molecule, the energy of this system can be expressed in terms of a double power series with respect to the bond length and to the field strength. From the power series expansion of the potential energy, field-dependent expressions for the equilibrium geometry, for the potential energy and for the force constant are obtained. The nuclear relaxation and vibrational contributions to the molecular electrical properties are analyzed in terms of the derivatives of the electronic molecular properties. In general, the results presented show that accurate inclusion of the correlation energy and large basis sets are needed to calculate the molecular electrical properties and their derivatives with respect to either nuclear displacements or/and field strength. With respect to experimental data, the calculated power series coefficients are overestimated by the SCF, CISD, and QCISD methods. On the contrary, perturbation methods (MP2 and MP4) tend to underestimate them. In average and using the 6-311 + G(3df) basis set and for the CO molecule, the nuclear relaxation and the vibrational contributions to the molecular electrical properties amount to 11.7%, 3.3%, and 69.7% of the purely electronic μ, α, and β values, respectively

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

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

Ab initio benchmark study for the oxidative addition of CH4 to Pd: importance of basis-set flexibility and polarization

To obtain a state-of-the-art benchmark potential energy surface (PES) for the archetypal oxidative addition of the methane C-H bond to the palladium atom, we have explored this PES using a hierarchical series of ab initio methods (Hartree-Fock, second-order Møller-Plesset perturbation theory, fourth-order Møller-Plesset perturbation theory with single, double and quadruple excitations, coupled cluster theory with single and double excitations (CCSD), and with triple excitations treated perturbatively [CCSD(T)]) and hybrid density functional theory using the B3LYP functional, in combination with a hierarchical series of ten Gaussian-type basis sets, up to g polarization. Relativistic effects are taken into account either through a relativistic effective core potential for palladium or through a full four-component all-electron approach. Counterpoise corrected relative energies of stationary points are converged to within 0.1-0.2 kcal/mol as a function of the basis-set size. Our best estimate of kinetic and thermodynamic parameters is -8.1 (-8.3) kcal/mol for the formation of the reactant complex, 5.8 (3.1) kcal/mol for the activation energy relative to the separate reactants, and 0.8 (-1.2) kcal/mol for the reaction energy (zero-point vibrational energy-corrected values in parentheses). This agrees well with available experimental data. Our work highlights the importance of sufficient higher angular momentum polarization functions, f and g, for correctly describing metal-d-electron correlation and, thus, for obtaining reliable relative energies. We show that standard basis sets, such as LANL2DZ+ 1f for palladium, are not sufficiently polarized for this purpose and lead to erroneous CCSD(T) results. B3LYP is associated with smaller basis set superposition errors and shows faster convergence with basis-set size but yields relative energies (in particular, a reaction barrier) that are ca. 3.5 kcal/mol higher than the corresponding CCSD(T) values

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.

Potential interstellar molecules : Formation of neutral C6CO from C6CO- in the gas phase

Computations at the RCCSD(T)/aug-cc-pVDZ//B3LYP/6-31G* level of theory indicate that neutral C6CO is a stable species. The ground state of this neutral is the singlet cumulene oxide :C=C=C=C=C=C=C=O. The adiabatic electron affinity and dipole moment of singlet C6CO are 2.47 eV and 4.13 D, respectively, at this level of theory. The anion (C6CO)(-.) should be a possible precursor to this neutral. It has been formed by an unequivocal synthesis in the ion source of a mass spectrometer by the S(N)2(Si) reaction between (CH3)(3)Si-C=C-C=C-C=C-CO-CMe3 and F- to form C-=C-C=C-C=C-CO-CMe3 which loses Me3C in the source to form C6CO-.. Charge stripping of this anion by vertical Franck-Condon oxidation forms C6CO, characterised by the neutralisation-reionisation spectrum (-NR+) of C6CO-., which is stable during the timeframe of this experiment (10(-6) s), Copyright (C) 2000 John Wiley & Sons, Ltd.

Generation of two isomers of C5H from the corresponding anions. A theoretically motivated mass spectrometric study

Molecular orbital calculations have predicted the stability of a range of connectivities for the radical C5H potential surface. The most energetically favorable of these include the linear C4CH geometry and two ring-chain structures HC2C3 and C2C3H The corresponding anions are also shown to be theoretically stable, and furthermore, a fourth isomer, C2CHC2, is predicted to be the most stable anion connectivity. These results have motivated experimental efforts. Methodologies for the generation of the non-ring-containing isomeric anions C4CH and C2CHC2 have been developed utilizing negative ion mass spectrometry. The absolute connectivities of the anions have been established using deuterium labeling, charge reversal, and neutralization reionization techniques. The success of the latter experiment confirms theoretical predictions of stability of the corresponding neutral species. This is the first reported observation of the neutral C2CHC2 species that calculations predict to be substantially less stable than the C4CH connectivity but still bound relative to isomerization processes.