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

Density functional energy decomposition into one- and two-atom contributions

The present work provides a generalization of Mayer's energy decomposition for the density-functional theory (DFT) case. It is shown that one- and two-atom Hartree-Fock energy components in Mayer's approach can be represented as an action of a one-atom potential VA on a one-atom density ρ A or ρ B. To treat the exchange-correlation term in the DFT energy expression in a similar way, the exchange-correlation energy density per electron is expanded into a linear combination of basis functions. Calculations carried out for a number of density functionals demonstrate that the DFT and Hartree-Fock two-atom energies agree to a reasonable extent with each other. The two-atom energies for strong covalent bonds are within the range of typical bond dissociation energies and are therefore a convenient computational tool for assessment of individual bond strength in polyatomic molecules. For nonspecific nonbonding interactions, the two-atom energies are low. They can be either repulsive or slightly attractive, but the DFT results more frequently yield small attractive values compared to the Hartree-Fock case. The hydrogen bond in the water dimer is calculated to be between the strong covalent and nonbonding interactions on the energy scale

Spin polarization, orbital occupation and band gap opening in vanadium dioxide: the effect of screened Hartree–Fock exchange

The metal–insulator transition of VO2 so far has evaded an accurate description by density functional theory. The screened hybrid functional of Heyd, Scuseria and Ernzerhof leads to reasonable solutions for both the low-temperature monoclinic and high-temperature rutile phases only if spin polarization is excluded from the calculations. We explore whether a satisfactory agreement with experiment can be achieved by tuning the fraction of Hartree Fock exchange (a) in the density functional. It is found that two branches of locally stable solutions exist for the rutile phase for 12:5% 6 a 6 20%. One is metallic and has the correct stability as compared to the monoclinic phase, the other is insulating with lower energy than the metallic branch. We discuss these observations based on the V 3d orbital occupations and conclude that a ¼ 10% is the best possible choice for spin-polarized VO2 calculations.

Probing the (6,7)Li nucleon densities through a new break-up process approach

In this work we present a double folding optical model analysis of new near-barrier quasi-elastic experimental data for the (6,7)Li + (120)Sn systems. From the analysis, it was possible to confirm the ground-state nucleon densities assumed for the weakly bound (6,7)Li isotopes. The apparent discrepancies between the experimental densities and those based on Dirac-Hartree-Fock Bogoliubov (DHB) calculations were removed. A new approach that simulates the projectile break-up and a positive polarization from couplings of (6,7)Li bound states with the continuum was considered in the reaction mechanism. (C) 2010 Elsevier B.V. All rights reserved.

Experimental and theoretical investigation of molecular structure and conformation of the 4-isopropylthioxanthone

In this study, the molecular structure and conformational analyses of the 4-isopropylthioxanthone (4-ITX) are reported according to experimental and theoretical results. The compound crystallizes in the centrosymmetric P (1) over bar space group with only one molecule in the asymmetric unit, presenting the most stable conformation, in which the three fused-rings adopt a planar geometry, and the isopropyl group assumes a torsional angle with less sterical hindrance. The structural and conformational analyses were performed using theoretical calculations such as Hartree-Fock (HF), DFT method in combination with 6-311G(d,p) and 6-31++G(d,p) and the results were compared with infrared spectroscopy (FT-IR) and X-ray diffraction (XRD). The supramolecular assembly of 4-ITX is kept by non-classical C-H center dot center dot center dot O hydrogen bonds and weak interactions such as pi-pi stacking. 4-ITX was also studied by (1)H and (13)C NMR spectroscopy. UV-Vis absorption spectroscopic properties of the 4-ITX showed the long-wavelength maximum shifts towards high energy when the solvent polarity increases. (C) 2011 Elsevier B.V. All rights reserved.

Exploring new molecular species on the (1)[H, Se, F] singlet potential energy surface: Energetics, structures, IR spectra, and heats of formation

Structural, energetic, and vibrational properties of new molecular species, HSeF and HFSe, the associated transition state, and dissociation fragments are investigated using a state-of-the-art theoretical approach, CCSD(T)/CBS. HSeF is a normal covalently bonded molecule 38.98 kcal mol (1) more stable than the complex HF-Se, which shows an unusual structure with a central fluorine atom and a bond angle of 101.8 degrees.A barrier (Delta G(#)) of 49.01 kcal mol (1) separates the two species. Vibrational frequencies are also quite distinct. Heats of formation are evaluated for the diatomic fragments and HSeF. Final Delta(f)H values depend on the experimental accuracy of those of Se(g) and H(2)Se. (c) 2009 Elsevier B.V. All rights reserved.

A quantum chemical study on a set of non-imidazole H(3) antihistamine molecules

Molecular orbital calculations were carried out on a set of 28 non-imidazole H(3) antihistamine compounds using the Hartree-Fock method in order to investigate the possible relationships between electronic structural properties and binding affinity for H3 receptors (pK(i)). It was observed that the frontier effective-for-reaction molecular orbital (FERMO) energies were better correlated with pK(i) values than highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) energy values. Exploratory data analysis through hierarchical cluster (HCA) and principal component analysis (PCA) showed a separation of the compounds in two sets, one grouping the molecules with high pK(i) values, the other gathering low pK(i) value compounds. This separation was obtained with the use of the following descriptors: FERMO energies (epsilon(FERMO)), charges derived from the electrostatic potential on the nitrogen atom (N(1)), electronic density indexes for FERMO on the N(1) atom (Sigma((FERMO))c(i)(2)). and electrophilicity (omega`). These electronic descriptors were used to construct a quantitative structure-activity relationship (QSAR) model through the partial least-squares (PLS) method with three principal components. This model generated Q(2) = 0.88 and R(2) = 0.927 values obtained from a training set and external validation of 23 and 5 molecules, respectively. After the analysis of the PLS regression equation and the values for the selected electronic descriptors, it is suggested that high values of FERMO energies and of Sigma((FERMO))c(i)(2), together with low values of electrophilicity and pronounced negative charges on N(1) appear as desirable properties for the conception of new molecules which might have high binding affinity. 2010 Elsevier Inc. All rights reserved.

GCHF basis sets and their application in the electronic structure study of PrMnO3

The approach called generator coordinate Hartree-Fock (GCHF) method is used in the selection of Gaussian basis set [25s18p for O ((3)p), 31s21p14d for Mn (S-6), and 33s22p16d9f for Pr ((4)J)] for atoms. The role of the weight functions in the assessment of the numerical integration range of the GCHF equations is shown. These basis sets are contracted to (25s18p/9s5p), (31s21p14d/9s6p4d), and (33s22pl6d9f118sl2p5d3f) by segmented contraction scheme of Dunning and they are utilized in calculations of Restricted-Open-HF (ROHF) Total and Orbital energies of the (MnO+1)-Mn-3 and (PrO+1)-Pr-1 fragments, to evaluate their quality in molecular studies. The addition of one d polarization function in the contracted (9s5p) basis set for O(P-3) atom and their application with the contracted (9s6p4d), (18s21p5d3f) basis sets for Mn (S-6) and Pr-Pr ((4)j) atoms lead to the electronic structure study of PrMnO3. The dipole moment, the total energy, and total atomic charges properties were calculated and were carried out at ROHF level with the [PrMnO3](2) fragment. The calculated values show that PrMnO3 does not present piezoelectric properties. (C) 2004 Elsevier B.V. All rights reserved.

Gaussian basis sets to the theoretical study of the electronic structure of perovskite (LaMnO3)

The Generator Coordinate Hartree-Fock (GCHF) method is applied to generate extended (20s14p), (30s19p13d), and (31s23p18d) Gaussian basis sets for the 0, Mn, and La atoms, respectively. The role of the weight functions (WFs) in the assessment of the numerical integration range of the GCHF equations is shown. These basis sets are then contracted to [5s3p] and [11s6p6d] for 0 and Mn atoms, respectively, and [17s11p7d] for La atom 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-Roothaan (HFR) method for (MnO1+)-Mn-5 and (LaO1+)-La-1 fragments. The results obtained with the contracted basis sets are compared with values obtained with the extended basis sets. The addition of one d polarization function in the contracted basis set for 0 atom and its utilization with the contracted basis sets for Mn and La atoms leads to the calculations of dipole moment and total atomic charges of perovskite (LaMnO3). The calculations were performed at the HFR level with the crystal [LaMnO3](2) fragment in space group C-2v the values of dipole moment, total energy, and total atomic charges showed that it is reasonable to believe that LaMnO3 presents behaviour of piezoelectric material. (C) 2003 Elsevier B.V. All rights reserved.

Investigation of piezoelectricity in perovskite (LaFeO3): A theoretical study

In a previous report we studied theoretically the piezoelectric effect in barium titanate (BaTiO3) [O. Treu Filho, J.C. Pinheiro, R.T. Kondo, J. Mol. Struct. (THEOCHEM), 671 (2004) 71]. In this article we applied the Hartree-Fock (HF) theory in the investigation of piezoelectricity in LaFeO3. Initially, the generator coordinate HF (GCHF) method was used to build 22s14p, 30s19p13d, and 32s24p17d Gaussian basis sets for O(3p), Fe(D-5), and La(D-2) atoms. Then those basis sets were contracted to [7s6p], [13s8p6d], and [18s13p7d], respectively. The quality of the contracted basis sets in polyatomic calculations was evaluated through calculations of total and orbital energies (HOMO and HOMO-1) of (FeO1+)-Fe-2 and (LaO1+)-La-1. Finally, the contracted basis sets were supplemented with polarization and diffuse functions and used to investigate the piezoelectricity in LaFeO3. The calculated properties were dipole moment, total energy, and atomic charges and the analysis of those properties showed that covalent bonds constitute the electronic structure of [LaFCO3](2) fragment. Therefore, it is reasonable to believe that LaFeO3 does not present piezoelectric properties. (C) 2006 Elsevier B.V. All rights reserved.

Theoretical analysis on TiO(2)(110)/V surface

Theoretical analysis based on the Hartree-Fock method were performed in order to study the stoichiometric TiO(2) (110) surface and the vanadium substituted system. The Pople with polarization 3-21G* basis set level was used. The TiO(2) (110) surface was modeled using a (TiO(2))(15) cluster model. In order to take into account the finite size of the cluster, we have studied two different models: the point charge and the hydrogen saturated methodologies. The charge values used in the point charge calculations were optimized. The density of states, orbital self-consistend field (SCF) energies, and Mulliken charge values were analyzed. The method and model's dependence on the analyzed results are discussed. The theoretical results are compared with available experimental data. (C) 2001 John Wiley & Sons, Inc.

Design of Gaussian basis sets to the theoretical interpretation of IR-spectrum of hexaaquachromium (III) ion, tetraoxochromium (IV) ion, and tetraoxochromium (VI) ion

The Generator Coordinate Hartree-Fock (GCHF) method is employed to design 16s, 16s10p, 24s17p13d, 25s17p13d, and 26s17p Gaussian basis sets for the H ((2)S), O ((3)P), O(2-) ((1)S), Cr(3+) ((4)F), Cr(4+) ((3)F), and Cr(6+) ((1)S) atomic species. These basis sets are then contracted to (4s) for H ((2)S), (6s4p) for O ((3)P), and O(2-) ((1)S), (986p3d) for Cr(3+) ((4)F), (10s8p3d) for Cr(4+) ((3)F), and (13s7p) for Cr(6+) (1S) by a standard procedure. For evaluation of the quality of those basis sets in molecular calculations, we have accomplished studies of total and orbital (HOMO and HOMO-1) energies at the HF-Roothaan level for the molecular species of our interest. The results obtained with the contracted basis sets are compared to the values obtained with our extended basis sets and to the standard 6-311G basis set from literature. Finally, the contracted basis sets are enriched with polarization function and then utilized in the theoretical interpretation of IR-spectrum of hexaaquachromium (III) ion, [Cr(H(2)O)(6)](3+), tetraoxochromium (IV) ion, [CrO(4)](4-), and tetraoxochromium (VI) ion, [CrO(4)](2-). The respective theoretical harmonic frequencies and IR-intensities were computed at the density functional theory (DFT) level. In the DFT calculations we employed the Becke's 1988 functional using the LYP correlation functional. The comparison between the results obtained and the corresponding experimental values indicates a very good description of the IR-spectra of the molecular ions studied, and that the GCHF method is still a legitimate alternative for selection of Gaussian basis sets. (C) 2003 Elsevier B.V. All rights reserved.

Basis sets applied to the theoretical study of the vibrational structure of hexaaquaaluminum(III) ion

Gaussian basis sets were developed with the Generator Coordinate Hartree-Fock (GCHF) method for the atoms from H (14s), O (23s16p), and Al (29sl9p) in the ground state. These basis sets were then contracted to 3s (12,1,1), 5s3p (18,2,1,1,1/14,1,1), and 7s5p (20,3,2,1,1,1,1/14,2,1,1,1) for H, O and Al atoms, respectively, by a standard procedure. The quality of contracted basis sets in molecular calculations was evaluated through studies of the total and orbital (epsilon(HOMO) and epsilon(HOMO-1)) energies at the HF level for the hexaaquaaluminum(III) ion, [AI(H(2)O)(6)](3+). For the O atom, the 5s3p was supplemented with d polarization function and it was used in combination with 3s, and 7s5p for H and Al atoms was used to the theoretical interpretation of the Infrared (IR) spectrum of hexaaquaaluminum(III) ion. The calculations of the IR-spectrum were also performed at the HF level and it showed that the basis sets obtained with the aid of GCHF method lead to the selection of useful contracted Gaussian basis sets for the theoretical study of vibrational property of ionic specie of our interest. (C) 2004 Elsevier B.V. All rights reserved.

All electron fully relativistic Dirac-Fock calculation for darmstadtium carbide using prolapse free basis set

For the first time, ab inito all electron fully relativistic and correlated Dirac-Fock calculations with prolapse free basis set are reported for a Super Heavy Element. We investigated the relativistic effects on bonding and on some spectroscopic constants for the darmstadtium carbide and our results at DF/CCSD(T) with a prolapse free basis set suggest for R-e, omega(e) and D-e the values of 174 pm, 1114 cm(-1) and 7.29 eV, respectively. These values are very similar to the values for PtC found on literature. It was also found that prolapse free basis set may be important to estimate the dissociation energy using Relativistic 4-components correlated methods. (c) 2007 ELsevier B.V. All rights reserved.

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.