The chemical Hamiltonian approach in density functional theory

The chemical Hamiltonian approach (CHA) for handling the basis set superposition error problem in intermolecular interactions has been implemented within density functional theory (DFT) using Gaussian atomic basis sets. As test examples, the potential curves of the water dimer were calculated using the Vosko-Wilk-Nusair, Becke-Perdew and Perdew exchange-correlation functionals. Comparisons with the counterpoise correction method show that CHA within DFT performs as well as previously for Hartree-Fock.     

First principles calculation of the mechanical compression of two organic molecular crystals

The mechanical compression curves for the organic molecular crystals 1,1-diamino-2,2-dinitroethylene and beta-octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (beta-HMX) are calculated using the Hartree-Fock approximation to the solutions of the many-body Schr?dinger equation for a periodic system as implemented in the computer program CRYSTAL. No correction was made for basis set superposition error. The equilibrium lattice parameters are reproduced to within 1% of reported experimental values. Pressure values on the isotherm also agree well with reported experimental values. To obtain accurate results, the relaxation of all the atomic coordinates as well as the lattice parameters under a fixed volume constraint was required.     

Density functional theory and molecular clusters

Density functional theory (DFT) methods, including nonlocal density gradient terms in the exchange and correlation energy functionals, were applied to various types of molecular clusters: H-bonded, ionic, electrostatic, and London. Reliable results on the structure and stabilization energy were obtained for the first two types of cluster as long as Becke3LYP and Becke3P86 functionals and basis sets of at least DZ + P quality were used. DFT methods with currently available functionals failed completely, however, for London-type clusters, for which no minimum was found on the potential energy surfaces. DFT interaction energy exhibits the same basis set extension dependence as the Hartree-Fock (HF) interaction energy. Therefore, the Boys-Bernardi function counterpoise procedure should be employed for elimination of the DFT basis set superposition error. © 1995 John Wiley & Sons, Inc.     

HX · MX3 molecules: heterodimer vapour ph ase complexes (M = B, Al; X = F, Cl) with halogen bridges An ab initio molecular orbital study

The equilibrium structures of HX · MX₃ vapour phase complexes are examined at the Hartree-Fock level using the DZP basis set. The complexes are found to be weak electron pair donor-acceptor adducts with long and weak intermolecular halogen bridges and with the hydrogens interacting with one of the halogen atoms of MX₃ in a cyclic arrangement. The energies of complex formation are calculated with and without correction for the basis set superposition error.     

Calculation of longitudinal polarizability and second hyperpolarizability of polyacetylene with the coupled perturbed Hartree-Fock/Kohn-Sham scheme: where it is shown how finite oligomer chains tend to the infinite periodic polymer

The longitudinal polarizability, α(xx), and second hyperpolarizability, γ(xxxx), of polyacetylene are evaluated by using the coupled perturbed Hartree-Fock/Kohn-Sham (HF/KS) scheme as implemented in the periodic CRYSTAL code and a split valence type basis set. Four different density functionals, namely local density approximation (LDA) (pure local), Perdew-Becke-Ernzerhof (PBE) (gradient corrected), PBE0, and B3LYP (hybrid), and the Hartree-Fock Hamiltonian are compared. It is shown that very tight computational conditions must be used to obtain well converged results, especially for γ(xxxx), that is, very sensitive to the number of k(->) points in reciprocal space when the band gap is small (as for LDA and PBE), and to the extension of summations of the exact exchange series (HF and hybrids). The band gap in LDA is only 0.01 eV: at least 300 k(->) points are required to obtain well converged total energy and equilibrium geometry, and 1200 for well converged optical properties. Also, the exchange series convergence is related to the band gap. The PBE0 band gap is as small as 1.4 eV and the exchange summation must extend to about 130 A? from the origin cell. Total energy, band gap, equilibrium geometry, polarizability, and second hyperpolarizability of oligomers -(C(2)H(2))(m)-, with m up to 50 (202 atoms), and of the polymer have been compared. It turns out that oligomers of that length provide an extremely poor representation of the infinite chain polarizability and hyperpolarizability when the gap is smaller than 0.2 eV (that is, for LDA and PBE). Huge differences are observed on α(xx) and γ(xxxx) of the polymer when different functionals are used, that is in connection to the well-known density functional theory (DFT) overshoot, reported in the literature about short oligomers: for the infinite model the ratio between LDA (or PBE) and HF becomes even more dramatic (about 500 for α(xx) and 10(10) for γ(xxxx)). On the basis of previous systematic comparisons of results obtained with various approaches including DFT, HF, Moller-Plesset (MP2) and coupled cluster for finite chains, we can argue that, for the infinite chain, the present HF results are the most reliable.     

On the counterpoise correction for the basis set superposition error in large systems

The appropriateness of the use of the counterpoise correction for the basis set superposition error in SCF calculations of the interaction energies for pairs of aliphatic amino acids is analyzed in this paper. Our results show that for this type of molecule where the magnitude of the basis set superposition error can become quite big, the use of the counterpoise method provides interaction energies in good agreement with near Hartree-Fock values. The inaccuracies associated with the counterpoise method are much less important compared with the basis set superposition error itself. It is shown that the use of a well-balanced minimal basis set together with the counterpoise method is a good compromise (quality versus computational cost) for calculating interaction energies in systems involving molecules of biological interest.     

Intermolecular potentials of the silane dimer calculated with Hartree-Fock theory, Møller-Plesset perturbation theory, and density functional theory

We have calculated the intermolecular interaction potentials of the silane dimer at the D3d conformation using the Hartree-Fock (HF) self-consistent theory, the correlation-corrected second-order M?ller-Plesset (MP2) perturbation theory, and the density functional theory (DFT) with 108 functionals chosen from the combinations of 9 exchange and 12 correlation functionals. Single-point coupled cluster [CCSD(T)] calculations have also been carried out to calibrate the correlation effect. The HF calculations yield unbound potentials largely because of the exchange-repulsion interaction. In the MP2 calculations, the basis set effects on the repulsion exponent, the equilibrium bond length, the binding energy, and the asymptotic behavior of the calculated intermolecular potentials have been thoroughly studied. We have employed basis sets from the Slater type orbitals fitted with Gaussian functions (STO-nG, n = 3 approximately 6), Pople's medium size basis sets [up to 6-311++G(3df,3pd)], to Dunning's correlation consistent basis sets (cc-pVXZ and aug-cc-pVXZ, X = D, T, Q). With increasing basis size, the repulsion exponent and the equilibrium bond length converge at the 6-31G** basis set and the 6-311++G(3d,3p) basis set, respectively, while a large basis set (aug-cc-pVTZ) is required to converge the binding energy at a chemical accuracy ( approximately 0.05 kcal/mol). Up to the largest basis set used, the asymptotic dispersion coefficient has not converged to the expected C6 value from molecular polarizability calculations. We attribute the slow convergence partly to the inefficacy of using the MP2 calculations with Gaussian type functions to model the asymptotic behavior. Both the basis set superposition error (BSSE) corrected and uncorrected results are presented to emphasize the importance of including such corrections. Only the BSSE corrected results systematically converge to the expected potential curve with increasing basis size. The DFT calculations generate a wide range of interaction patterns, from purely unbound to strongly bound, underestimating or overestimating the binding energy. The binding energies calculated using the OPTXHCTH147, PBEVP86, PBEP86, PW91TPSS, PW91PBE, and PW91PW91 functionals and the equilibrium bond lengths calculated using the MPWHCTH93, TPSSHCTH, PBEVP86, PBEP86, PW91TPSS, PW91PBE, and PW91PW91 functionals are close to the MP2 results using the 6-311++G(3df,3pd) basis set. A correlation between the calculated DFT potentials and the exchange and correlation enhancement factors at the low-density region has been elucidated. The asymptotic behaviors of the DFT potentials are also analyzed.     

Basis set dependence of solute-solvent interaction energy of benzene in water: a HF/DFT study

Solute-solvent interaction energies for the benzene molecule dissolved in water are computed using Hartree-Fock and B3LYP density functional theories. Explicit solvent molecules up to 14-A away from the dissolved benzene molecule are included in the calculation of interaction energies. Both basis set dependence and basis set superposition errors are carefully examined. It is found that the use of a larger basis set for the region near the solute together with a smaller basis set for the outer region gives results very close to what would have been obtained if the larger basis set had been used for the whole system. It is also shown that a correction for the basis set superposition error is a necessary component in this kind of calculations. With this correction, results obtained with different tested basis sets converge to within 1 kcal/mol.     

Intermolecular potentials of the methane dimer calculated with Møller-Plesset perturbation theory and density functional theory

We have calculated the intermolecular interaction potentials of the methane dimer at the minimum-energy D(3d) conformation using the Hartree-Fock (HF) self-consistent theory, the correlation-corrected second-order M?ller-Plesset (MP2) perturbation theory, and the density functional theory (DFT) with the Perdew-Wang (PW91) functional as the exchange or the correlation part. The HF calculations yield unbound potentials largely due to the exchange-repulsion interaction. In the MP2 calculations, the basis set effects on the repulsion exponent, the equilibrium bond length, the binding energy, and the asymptotic behavior of the calculated intermolecular potentials have been thoroughly studied. We have employed basis sets from the Slater-type orbitals fitted with Gaussian functions (STO-nG) (n=3-6) [Quantum Theory of Molecular and Solids: The Self-Consistent Field for Molecular and Solids (McGraw-Hill, New York, 1974), Vol. 4], Pople's medium size basis sets of Krishnan et al. [J. Chem. Phys. 72, 650 (1980)] [up to 6-311++G(3df,3pd)] to Dunning's correlation consistent basis sets [J. Chem. Phys. 90, 1007 (1989)] (cc-pVXZ and aug-cc-pVXZ) (X=D, T, and Q). With increasing basis size, the repulsion exponent and the equilibrium bond length converge at the 6-31G** basis set and the 6-311++G(2d,2p) basis set, respectively, while a large basis set (aug-cc-pVTZ) is required to converge the binding energy at a chemical accuracy (approximately 0.01 kcal/mol). Up to the largest basis set used, the asymptotic dispersion coefficient has not converged to the destined C6 value from molecular polarizability calculations. The slow convergence could indicate the inefficacy of using the MP2 calculations with Gaussian-type functions to model the asymptotic behavior. Both the basis set superposition error (BSSE) corrected and uncorrected results are presented to emphasize the importance of including such corrections. Only the BSSE corrected results systematically converge to the destined potential curve with increasing basis size. The DFT calculations generate a wide range of interaction patterns, from purely unbound to strongly bound, underestimating or overestimating the binding energy. The binding energy calculated using the PW91PW91 functional and the equilibrium bond length calculated using the PW91VP86 functional are close to the MP2 results at the basis set limit.     

去氢抗坏血酸分子振动光谱的理论研究

采用RHF, MP2, DFT(B3LYP)方法, 以6-311++G**为基组研究了去氢抗坏血酸分子(DHA)的平衡几何构型和振动光谱. 计算结果表明, 采用RHF, B3LYP以及MP2 方法优化得到的几何结构以及频率值是一致的. 采用B3LYP/6-311++G**计算了DHA分子平衡构型下的谐振动力场﹑振动频率和振动强度. 使用Wilson的GF矩阵方法对DHA分子进行了简正坐标分析, 依据所得的势能分布对DHA分子的振动基频进行了合理的理论归属.     

The HeHe Coulomb and exchange interaction energy

The Coulomb and exchange interaction energy between two helium atoms is calculated, increasing systematically the basis set, by a supermolecular method which excludes the superposition error by using non-orthogonal orbitals. At 5.6 bohr an energy limit of 9.69 K is found at the Hartree-Fock level; for correlated wavefunctions the corresponding value is 11.86 K using only s orbitals, 10.70 K including p orbitals and 10.58 K adding d orbitals.     

Hartree-Fock complete basis set limit properties for transition metal diatomics

Numerical Hartree-Fock (HF) energies accurate to at least 1 microhartree are reported for 27 diatomic transition-metal-containing species. The convergence of HF energies toward this numerical limit upon increasing the basis set size has been investigated, where standard nonrelativistic all-electron correlation consistent basis sets and augmented basis sets, developed by Balabanov and Peterson [J. Chem. Phys. 123, 064107 (2005)], were employed. Several schemes which enable the complete basis set (CBS) limit to be determined have been investigated, and the resulting energies have been compared to the numerical Hartree-Fock energies. When comparing basis set extrapolation schemes, those in the form of exponential functions perform well for our test set, with mean absolute deviations from numerical HF energies of 234 and 153 microE(h), when the CBS limit has been determined using a two-point fit as proposed by Halkier et al. [Chem. Phys. Lett. 302, 437 (1999)] on calculations of triple- and quadruple-zeta basis set qualities and calculations of quadruple- and quintuple-zeta basis set qualities, respectively. Overall, extrapolation schemes in the form of a power series are not recommended for the extrapolation of transition metal HF energies. The impact of basis set superposition error has also been examined.     

A theoretical study on the intermolecular interaction of energetic system-Nitromethane dimer

Three optimized geometries of nitromethane dimer have been obtained at the HF/6-31G level.Dimer binding energies have been corrected for the basis set superposition error (BSSE) and the zero point energy.Computed results indicate that the cyclic structure of (CH3NO2)2 is the most stable of three optimized geometries,whose corrected binding energyis 17.29 kJ mol-1 at the MP4SDTQ/6-31G//HF/6-31G level.In the optimized structures of nitromethane dimer,the inter-molecular hydrogen bond has not been found; and the charge-transfer interaction between CH3NO2 subsystems is weak; and the correlation interaction energy makes a little contribution to the intermolecular interaction energy of the dimer.     

Hydrogen bonding in the urea dimers and adenine–thymine DNA base pair: anharmonic effects in the intermolecular H-bond and intramolecular H-stretching vibrations

The equilibrium structures, binding energies, vibrational harmonic frequencies, and the anharmonic corrections for two different (cyclic and asymmetric) urea dimers and for the adenine–thymine DNA base pair system have been studied using the second-order Møller–Plesset perturbation theory (MP2) method and different density functional theory (DFT) exchange–correlation (XC) functionals (BLYP, B3LYP, PBE, HCTH407, KMLYP, and BH and HLYP) with the D95V, D95V**, and D95V++** basis sets. The widely used a posteriori Boys–Bernardi or counterpoise correction scheme for basis set superposition error (BSSE) has been included in the calculations to take into account the BSSE effects during geometry optimization (on structure), on binding energies and on the different levels of approximation used for calculating the vibrational frequencies. The results obtained with the ab initio MP2 method are compared with those calculated with different DFT XC functionals; and finally the suitability of these DFT XC functionals to describe intermolecular hydrogen bonds as well as harmonic frequencies and the anharmonic corrections is assessed and discussed.     

Ab initio study of the proton affinity of a number of ortho-substituted pyridines

The proton affinity, dipole moment, and first vertical ionization potential of a number of ortho-mono-substituted pyridines are computed at the restricted Hartree-Fock level, using a number of different basis sets. Some approximations concerning the geometry, the force constants, and the basis set superposition error are introduced. The effect of both the basis set and the approximations is investigated systematically. Two nonstandard basis sets, 4-31G(*) and 4-31+G(*), in which polarization functions are only applied to heteroatoms, are defined. A qualitative explanation for the effects of electron correlation on the proton affinity is given and underbuilt with the results of some MP2 calculations. While accuracy is maintained, computational cost is reduced very substantially.     

CH3NO2+X(X=H,OH, CH3, CH2[^3B1], O[^3P])→CH2NO2+XH吸氢反应过渡态 结 构和反应位垒的DFT计算研究

采用DFT(B3LYP)方法,分别在6-311g(d,p),6-311++g(d,p)和自洽相关基组cc-pVIZ水平上优化了基态硝基甲烷和自由基H,OH,CH3,CH2[^3B1]以及O[^3P]等发生吸氢反应时的过渡态结构,并计算了反应的位垒。研究表明,对同一反应,不同基组下优化得到的过渡态几何结构基本一致;反应位垒数值的大小也基本接近,经校正,硝基甲烷同自由基反应位垒的理论计算值同实验结果基本吻合。     

2-羟基吡啶与水氢键作用的理论研究

本文采用量子化学的Hatree-Fock方法和密度泛函理论(DFT)的B3LYP方法,在6-31G(d)水平上,研究了2-羟基吡啶分子(Hy)及其酮式互变异构体2(1H)-吡啶酮(Py)与水的相互作用。考察它们之间在形成Hy…H2O,Py…H2O,Hy…Hy,Py…Py和Hy…Py等复合物前后的能量变化和分子结构参数变化特点。计算结果表明,在这些复合物中都形成了较强的氢键作用,在水合物中,Py与水形成复合物时能量降低较多,与实验结果一致。经过零点振动能(ZPVE)和基组叠加误差(BSSE)校正后的复合物离解能分别为38.3,40.8,73.0,82.7和71.1 kJ/mol(B3LYP/6-31G(d)),水合物的离解能远小于二聚体复合物,而酮式结构的二聚体的离解能最大。     

Estimating the basis set superposition error in the CCSD(T)(F12) explicitly correlated method using the example of a water dimer

Changes in the basis set superposition errors upon transitioning from conventional CCSD(T) to the CCSD(T)(F12) explicitly correlated method is studied using the example of a water dimer. A comparison of the compensation errors for CCSD(T) and CCSD(T)(F12) reveals a substantial reduction in the superposition error upon use of the latter. Numerical experiments with water dimers show it is possible theoretically predict an equilibrium distance between oxygen atoms that is similar to the experimental data (2.946 Å), as is the predicted energy of dissociation of a dimer (5.4 ± 0.7 kcal/mol). It is found that the structural and energy parameters of hydrogen bonds in water dimers can be calculated precisely even with two-exponential correlation-consistent basis sets if we use the explicitly correlated approach and subsequently correct the basis set superposition error.     

Basis sets for molecular interactions. 2. Application to H3N?HF,H3N?HOH,H2O?HF, (NH3)2, and H3CH?OH2

Modifications of the standard 6-31G** basis set as recommended in the accompanying paper are found to markedly lower the basis set superposition error (BSSE) in the title complexes, in contrast to enlargement to a triple-ζ scheme or by addition of a diffuse sp shell or a second set of d-functions without prior optimization, all of which lead to BSSE increase. After appropriate correction for correlation and superposition effects, all basis sets (with the exception of the standard 6-31G** and 6-311G** with their very large BSSE) predict the cyclic geometry of NH₃ dimer to be more stable than the linear arrangement. Correlation and BSSE can shift the equilibrium intermolecular distance in H₃CH-OH₂ by up to 0.4 Å. Failure to correct for superposition error leads to a drastic exaggeration of both the SCF and MP2 components of the interaction energy in this complex. Much better estimates are furnished by our recommended basis sets with their smaller superposition errors.     

Basis set effects on the intermolecular interaction of the H2-H2 system obtained using ab initio molecular orbital calculations with the Møller-Plesset perturbation correction

The intermolecular interaction potential of the H₂-H₂ system was calculated by an ab initio molecular orbital method using several basis sets (up to 6-31 lG(3pd)) with inclusion of the electron correlation correction of the Møller-Plesset perturbation method and the basis set superposition error (BSSE) correction of the counterpoise method in order to evaluate the basis set effect. The calculated interaction energies depend strongly on the basis set used. Whereas the interaction energies of the repulsive and coulombic energy components calculated at the Hartree-Fock level are not affected by a change of basis set, the dispersion energy component depends strongly on the basis set used. Parameters of an exp-6-1 type non-bonding interaction potential were optimized on the basis of the MP4(SDTQ)/6-311G(3p) level intermolecular interaction energies of the H₂-H₂ system.