451
|
Cisneros GA, Piquemal JP, Darden TA. Intermolecular electrostatic energies using density fitting. J Chem Phys 2007; 123:044109. [PMID: 16095348 PMCID: PMC2693352 DOI: 10.1063/1.1947192] [Citation(s) in RCA: 53] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A method is presented to calculate the electron-electron and nuclear-electron intermolecular Coulomb interaction energy between two molecules by separately fitting the unperturbed molecular electron density of each monomer. This method is based on the variational Coulomb fitting method which relies on the expansion of the ab initio molecular electron density in site-centered auxiliary basis sets. By expanding the electron density of each monomer in this way the integral expressions for the intermolecular electrostatic calculations are simplified, lowering the operation count as well as the memory usage. Furthermore, this method allows the calculation of intermolecular Coulomb interactions with any level of theory from which a one-electron density matrix can be obtained. Our implementation is initially tested by calculating molecular properties with the density fitting method using three different auxiliary basis sets and comparing them to results obtained from ab initio calculations. These properties include dipoles for a series of molecules, as well as the molecular electrostatic potential and electric field for water. Subsequently, the intermolecular electrostatic energy is tested by calculating ten stationary points on the water dimer potential-energy surface. Results are presented for electron densities obtained at four different levels of theory using two different basis sets, fitted with three auxiliary basis sets. Additionally, a one-dimensional electrostatic energy surface scan is performed for four different systems (H2O dimer, Mg2+-H2O, Cu+-H2O, and n-methyl-formamide dimer). Our results show a very good agreement with ab initio calculations for all properties as well as interaction energies.
Collapse
Affiliation(s)
- G Andrés Cisneros
- Laboratory of Structural Biology, National Institute of Environmental Health Sciences, Research Triangle Park, North Carolina 27707, USA.
| | | | | |
Collapse
|
452
|
Gerber IC, Angyán JG. London dispersion forces by range-separated hybrid density functional with second order perturbational corrections: The case of rare gas complexes. J Chem Phys 2007; 126:044103. [PMID: 17286458 DOI: 10.1063/1.2431644] [Citation(s) in RCA: 70] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
A satisfactory account of the van der Waals (vdW) (London dispersion) forces is, in general not possible by the Kohn-Sham method using standard local, semilocal generalized gradient approximation (GGA), or meta-GGA density functionals. The recently proposed range-separated hybrid (RSH) approach, supplemented by second order perturbational corrections (MP2) to include long-range dynamic correlation effects, offers a physically consistent, seamless description of dispersion forces. It is based on a rigorous generalization of the Kohn-Sham method, where long-range exchange and correlation effects are treated by wave function methods, while short-range electron exchange and correlation are handled by local or semilocal functionals. The method is tested on a series of rare gas dimers in comparison with standard wave function theory and density functional theory approaches. In contrast to the most successful exchange correlation functionals, which describe at best the vdW minimum, the RSH+MP2 approach is valid also in the asymptotic region and the potential curve displays the correct 1/R(6) behavior at large internuclear separations. In contrast to usual MP2 calculations, the basis set superposition error is considerably reduced, making RSH+MP2 an ideal tool for exploring the potential energy surface of weakly bound molecular complexes.
Collapse
Affiliation(s)
- I C Gerber
- Laboratoire de Cristallographie et de Modélisation des Matériaux Minéraux et Biologiques, UMR 7036, CNRS, Nancy Université, B.P. 239, F-54506 Vandoeuvre-lès-Nancy, France
| | | |
Collapse
|
453
|
|
454
|
Pluhácková K, Jurecka P, Hobza P. Stabilisation energy of C6H6⋯C6X6(X = F, Cl, Br, I, CN) complexes: complete basis set limit calculations at MP2 and CCSD(T) levels. Phys Chem Chem Phys 2007; 9:755-60. [PMID: 17268688 DOI: 10.1039/b615318f] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Stabilisation energies of stacked structures of C(6)H(6)...C(6)X(6) (X = F, Cl, Br, CN) complexes were determined at the CCSD(T) complete basis set (CBS) limit level. These energies were constructed from MP2/CBS stabilisation energies and a CCSD(T) correction term determined with a medium basis set (6-31G**). The former energies were extrapolated using the two-point formula of Helgaker et al. from aug-cc-pVDZ and aug-cc-pVTZ Hartree-Fock energies and MP2 correlation energies. The CCSD(T) correction term is systematically repulsive. The final CCSD(T)/CBS stabilisation energies are large, considerably larger than previously calculated and increase in the series as follows: hexafluorobenzene (6.3 kcal mol(-1)), hexachlorobenzene (8.8 kcal mol(-1)), hexabromobenzene (8.1 kcal mol(-1)) and hexacyanobenzene (11.0 kcal mol(-1)). MP2/SDD** relativistic calculations performed for all complexes mentioned and also for benzene[dot dot dot]hexaiodobenzene have clearly shown that due to relativistic effects the stabilisation energy of the hexaiodobenzene complex is lower than that of hexabromobenzene complex. The decomposition of the total interaction energy to physically defined energy components was made by using the symmetry adapted perturbation treatment (SAPT). The main stabilisation contribution for all complexes investigated is due to London dispersion energy, with the induction term being smaller. Electrostatic and induction terms which are attractive are compensated by their exchange counterparts. The stacked motif in the complexes studied is very stable and might thus be valuable as a supramolecular synthon.
Collapse
Affiliation(s)
- Kristýna Pluhácková
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic and Centre for Biomolecules and Complex Molecular Systems, Flemingovo nám 2, 166 10, Prague 6, Czech Republic
| | | | | |
Collapse
|
455
|
Podeszwa R, Bukowski R, Rice BM, Szalewicz K. Potential energy surface for cyclotrimethylene trinitramine dimer from symmetry-adapted perturbation theory. Phys Chem Chem Phys 2007; 9:5561-9. [DOI: 10.1039/b709192c] [Citation(s) in RCA: 67] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
456
|
Tekin A, Jansen G. How accurate is the density functional theory combined with symmetry-adapted perturbation theory approach for CH–π and π–π interactions? A comparison to supermolecular calculations for the acetylene–benzene dimer. Phys Chem Chem Phys 2007; 9:1680-7. [PMID: 17396179 DOI: 10.1039/b618997k] [Citation(s) in RCA: 87] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Five different orientations of the acetylene-benzene dimer including the T-shaped global minimum structure are used to assess the accuracy of the density functional theory combined with symmetry adapted perturbation theory (DFT-SAPT) approach in its density-fitting implementation (DF-DFT-SAPT) for the study of CH-pi and pi-pi interactions. The results are compared with the outcome of counterpoise corrected supermolecular calculations employing second-order Møller-Plesset (MP2), spin-component scaled MP2 (SCS-MP2) and single and double excitation coupled cluster theory including perturbative triple excitations (CCSD(T)). For all considered orientations MP2 predicts much deeper potential energy curves with considerably shifted minima compared to CCSD(T) and DFT-SAPT. In spite of being an improvement over the results of MP2, SCS-MP2 tends to underestimate the well depth while DFT-SAPT, employing an asymptotically corrected hybrid exchange-correlation potential in conjunction with the adiabatic local density approximation for the exchange-correlation kernel, is found to be in excellent agreement with CCSD(T). Furthermore, DFT-SAPT provides a detailed understanding of the importance of the electrostatic, induction and dispersion contributions to the total interaction energy and their repulsive exchange corrections.
Collapse
Affiliation(s)
- Adem Tekin
- Fachbereich Chemie, Universität Duisburg-Essen, Universitätsstrasse 5, 45117, Essen, Germany
| | | |
Collapse
|
457
|
Grimme S, Antony J, Schwabe T, Mück-Lichtenfeld C. Density functional theory with dispersion corrections for supramolecular structures, aggregates, and complexes of (bio)organic molecules. Org Biomol Chem 2007; 5:741-58. [PMID: 17315059 DOI: 10.1039/b615319b] [Citation(s) in RCA: 519] [Impact Index Per Article: 30.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Kohn-Sham density functional theory (KS-DFT) is nowadays the most widely used quantum chemical method for electronic structure calculations in chemistry and physics. Its further application in e.g. supramolecular chemistry or biochemistry has mainly been hampered by the inability of almost all current density functionals to describe the ubiquitous attractive long-range van der Waals (dispersion) interactions. We review here methods to overcome this defect, and describe in detail a very successful correction that is based on damped -C(6).R(-6) potentials (DFT-D). As examples we consider the non-covalent inter- and intra-molecular interactions in unsaturated organic molecules (so-called pi-pi stacking in benzenes and dyes), in biologically relevant systems (nucleic acid bases/pairs, proteins, and 'folding' models), between fluorinated molecules, between curved aromatics (corannulene and carbon nanotubes) and small molecules, and for the encapsulation of methane in water clusters. In selected cases we partition the interaction energies into the most relevant contributions from exchange-repulsion, electrostatics, and dispersion in order to provide qualitative insight into the binding character.
Collapse
Affiliation(s)
- Stefan Grimme
- Theoretische Organische Chemie, Organisch-Chemisches Institut der Universität Münster, Corrensstrasse 40, D-48149 Münster, Germany.
| | | | | | | |
Collapse
|
458
|
Hill JG, Platts JA. Spin-Component Scaling Methods for Weak and Stacking Interactions. J Chem Theory Comput 2006; 3:80-5. [DOI: 10.1021/ct6002737] [Citation(s) in RCA: 132] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- J. Grant Hill
- School of Chemistry, Cardiff University, Park Place, Cardiff CF10 3AT, U.K
| | - James A. Platts
- School of Chemistry, Cardiff University, Park Place, Cardiff CF10 3AT, U.K
| |
Collapse
|
459
|
Korona T, Jeziorski B. One-electron properties and electrostatic interaction energies from the expectation value expression and wave function of singles and doubles coupled cluster theory. J Chem Phys 2006; 125:184109. [PMID: 17115740 DOI: 10.1063/1.2364489] [Citation(s) in RCA: 57] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
One-electron density matrices resulting from the explicitly connected commutator expansion of the expectation value were implemented at the singles and doubles coupled cluster (CCSD) level. In the proposed approach the one-electron density matrix is obtained at a little extra cost in comparison to the calculation of the CCSD correlation energy. Therefore, in terms of the computational time the new method is significantly less demanding than the conventional linear-response CCSD theory which requires additionally an expensive calculation of the left-hand solution of the CCSD equations. The quality of the new density matrices was investigated by computing a set of one-electron properties for a series of molecules of varying sizes and comparing the results with data obtained using the full configuration interaction method or higher level coupled cluster theory. It has been found that the results obtained using the new approach are of the same quality as those predicted by the linear-response CCSD method. The novel one-electron density matrices have also been applied to study the energy of the electrostatic interaction for a number of van der Waals complexes, including the benzene and azulene dimers.
Collapse
Affiliation(s)
- Tatiana Korona
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland.
| | | |
Collapse
|
460
|
Patkowski K, Szalewicz K, Jeziorski B. Third-order interactions in symmetry-adapted perturbation theory. J Chem Phys 2006; 125:154107. [PMID: 17059239 DOI: 10.1063/1.2358353] [Citation(s) in RCA: 89] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
We present an extension of many-body symmetry-adapted perturbation theory (SAPT) by including all third-order polarization and exchange contributions obtained with the neglect of intramonomer correlation effects. The third-order polarization energy, which naturally decomposes into the induction, dispersion, and mixed, induction-dispersion components, is significantly quenched at short range by electron exchange effects. We propose a decomposition of the total third-order exchange energy into the exchange-induction, exchange-dispersion, and exchange-induction-dispersion contributions which provide the quenching for the corresponding individual polarization contributions. All components of the third-order energy have been expressed in terms of molecular integrals and orbital energies. The obtained formulas, valid for both dimer- and monomer-centered basis sets, have been implemented within the general closed-shell many-electron SAPT program. Test calculations for several small dimers have been performed and their results are presented. For dispersion-bound dimers, the inclusion of the third-order effects eliminates the need for a hybrid SAPT approach, involving supermolecular Hartree-Fock calculations. For dimers consisting of strongly polar monomers, the hybrid approach remains more accurate. It is shown that, due to the extent of the quenching, the third-order polarization effects should be included only together with their exchange counterparts. Furthermore, the latter have to be calculated exactly, rather than estimated by scaling the second-order values.
Collapse
Affiliation(s)
- Konrad Patkowski
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA.
| | | | | |
Collapse
|
461
|
Ringer AL, Figgs MS, Sinnokrot MO, Sherrill CD. Aliphatic C−H/π Interactions: Methane−Benzene, Methane−Phenol, and Methane−Indole Complexes. J Phys Chem A 2006; 110:10822-8. [PMID: 16970377 DOI: 10.1021/jp062740l] [Citation(s) in RCA: 165] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Noncovalent C-H/pi interactions are prevalent in biochemistry and are important in molecular recognition. In this work, we present potential energy curves for methane-benzene, methane-phenol, and methane-indole complexes as prototypes for interactions between C-H bonds and the aromatic components of phenylalanine, tyrosine, and tryptophan. Second-order perturbation theory (MP2) is used in conjunction with the aug-cc-pVDZ and aug-cc-pVTZ basis sets to determine the counterpoise-corrected interaction energy for selected complex configurations. Using corrections for higher-order electron correlation determined with coupled-cluster theory through perturbative triples [CCSD(T)] in the aug-cc-pVDZ basis set, we estimate, through an additive approximation, results at the very accurate CCSD(T)/aug-cc-pVTZ level of theory. Symmetry-adapted perturbation theory (SAPT) is employed to determine the physically significant components of the total interaction energy for each complex.
Collapse
Affiliation(s)
- Ashley L Ringer
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | | | | | | |
Collapse
|
462
|
Podeszwa R, Bukowski R, Szalewicz K. Potential Energy Surface for the Benzene Dimer and Perturbational Analysis of π−π Interactions. J Phys Chem A 2006; 110:10345-54. [PMID: 16928128 DOI: 10.1021/jp064095o] [Citation(s) in RCA: 297] [Impact Index Per Article: 16.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a complete 6-dimensional potential energy surface for the benzene dimer obtained using symmetry-adapted perturbation theory (SAPT) of intermolecular interactions based on Kohn-Sham's description of monomers. Ab initio calculations were performed for 491 dimer geometries in a triple-zeta-quality basis set supplemented by bond functions. An accurate analytic fit to the ab initio results has been developed and low-energy stationary points on the potential energy surface have been found. We have determined that there are three minima on the surface. Two of them, the tilted T-shape and the parallel-displaced, are nearly isoenergetic with interaction energies of -2.77 and -2.74 kcal/mol, respectively. The third minimum, a twisted edge-to-edge conformation, is significantly less attractive, with the interaction energy equal to -1.82 kcal/mol. Both the T-shape and sandwich geometries, sometimes assumed to be minima, are shown to be only saddle points. The potential energy surface is extremely flat between the two lowest minima, the barrier being only 0.10 kcal/mol above the global minimum. The second-virial coefficient obtained with the new potential agrees well with experimental results over a wide range of temperatures. The SAPT approach rigorously decomposes the interaction energy into physical components. The relative importance of these components has been analyzed.
Collapse
Affiliation(s)
- Rafał Podeszwa
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| | | | | |
Collapse
|
463
|
Sinnokrot MO, Sherrill CD. High-Accuracy Quantum Mechanical Studies of π−π Interactions in Benzene Dimers. J Phys Chem A 2006; 110:10656-68. [PMID: 16970354 DOI: 10.1021/jp0610416] [Citation(s) in RCA: 588] [Impact Index Per Article: 32.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Although supramolecular chemistry and noncovalent interactions are playing an increasingly important role in modern chemical research, a detailed understanding of prototype noncovalent interactions remains lacking. In particular, pi-pi interactions, which are ubiquitous in biological systems, are not fully understood in terms of their strength, geometrical dependence, substituent effects, or fundamental physical nature. However, state-of-the-art quantum chemical methods are beginning to provide answers to these questions. Coupled-cluster theory through perturbative triple excitations in conjunction with large basis sets and extrapolations to the complete basis set limit have provided definitive results for the binding energy of several configurations of the benzene dimer, and benchmark-quality ab initio potential curves are being used to calibrate new density functional and force-field models for pi-pi interactions. Studies of substituted benzene dimers indicate flaws in the conventional wisdom about substituent effects in pi-pi interactions. Three-body and four-body interactions in benzene clusters have also been examined.
Collapse
Affiliation(s)
- Mutasem Omar Sinnokrot
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | | |
Collapse
|
464
|
Multicentered integrated QM:QM methods for weakly bound clusters: An efficient and accurate 2-body:many-body treatment of hydrogen bonding and van der Waals interactions. Chem Phys Lett 2006. [DOI: 10.1016/j.cplett.2006.06.021] [Citation(s) in RCA: 58] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
465
|
Bukowski R, Szalewicz K, Groenenboom G, van der Avoird A. Interaction potential for water dimer from symmetry-adapted perturbation theory based on density functional description of monomers. J Chem Phys 2006; 125:44301. [PMID: 16942136 DOI: 10.1063/1.2220040] [Citation(s) in RCA: 61] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A new six-dimensional interaction potential for the water dimer has been obtained by fitting interaction energies computed at 2510 geometries using a variant of symmetry-adapted perturbation theory (SAPT) based on density functional theory (DFT) description of monomers, referred to as SAPT(DFT). The stationary points, second virial coefficient, vibration-rotation-tunneling spectrum, and structure of liquid water computed with the new potential are in very good agreement with experiment and advanced ab initio calculations, confirming the high level of accuracy provided by SAPT(DFT).
Collapse
Affiliation(s)
- Robert Bukowski
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| | | | | | | |
Collapse
|
466
|
Torheyden M, Jansen G. A new potential energy surface for the water dimer obtained from separate fits ofab initioelectrostatic, induction, dispersion and exchange energy contributions. Mol Phys 2006. [DOI: 10.1080/00268970600679188] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
467
|
Söderhjelm P, Karlström G, Ryde U. Comparison of overlap-based models for approximating the exchange-repulsion energy. J Chem Phys 2006; 124:244101. [PMID: 16821967 DOI: 10.1063/1.2206182] [Citation(s) in RCA: 33] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Different ways of approximating the exchange-repulsion energy with a classical potential function have been investigated by fitting various expressions to the exact exchange-repulsion energy for a large set of molecular dimers. The expressions involve either the orbital overlap or the electron-density overlap. For comparison, the parameter-free exchange-repulsion model of the effective fragment potential (EFP) is also evaluated. The results show that exchange-repulsion energy is nearly proportional to both the orbital overlap and the density overlap. For accurate results, a distance-dependent correction is needed in both cases. If few parameters are desired, orbital overlap is superior to density overlap, but the fit to density overlap can be significantly improved by introducing more parameters. The EFP performs well, except for delocalized pi systems. However, an overlap expression with a few parameters seems to be slightly more accurate and considerably easier to approximate.
Collapse
Affiliation(s)
- Pär Söderhjelm
- Department of Theoretical Chemistry, Chemical Center, Lund University, P.O. Box 124, S-22100 Lund, Sweden.
| | | | | |
Collapse
|
468
|
Beran GJO, Head-Gordon M, Gwaltney SR. Second-order correction to perfect pairing: An inexpensive electronic structure method for the treatment of strong electron-electron correlations. J Chem Phys 2006; 124:114107. [PMID: 16555874 DOI: 10.1063/1.2176603] [Citation(s) in RCA: 32] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
We have formulated a second-order perturbative correction for perfect-pairing wave functions [PP2] based on similarity-transformed perturbation techniques in coupled cluster theory. The perfect-pairing approximation is used to obtain a simple reference wave function which can qualitatively describe bond breaking, diradicals, and other highly correlated systems, and the perturbative correction accounts for the dynamical correlation. An efficient implementation of this correction using the resolution of the identity approximation enables PP2 to be computed at a cost only a few times larger than that of canonical MP2 for systems with hundreds of active electrons and tens of heavy atoms. PP2 significantly improves on MP2 predictions in various systems with a challenging electronic structure.
Collapse
Affiliation(s)
- Gregory J O Beran
- Department of Chemistry, University of California, Berkeley, California 94720-1460, USA
| | | | | |
Collapse
|
469
|
Misquitta AJ, Stone AJ. Distributed polarizabilities obtained using a constrained density-fitting algorithm. J Chem Phys 2006; 124:024111. [PMID: 16422575 DOI: 10.1063/1.2150828] [Citation(s) in RCA: 73] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
A computationally efficient method for obtaining distributed polarizabilities of arbitrary rank using a constrained density-fitting algorithm is demonstrated on the hydrogen, carbon dioxide, formamide, and N-methylpropanamide molecules. A description of the molecular polarization in terms of local polarizabilities without charge-flow terms is obtained when the nonlocal components of the polarizability tensor are transformed away using the localization method of Le Sueur and Stone [Mol. Phys. 83, 293 (1994)]. The resulting local polarizabilities are shown to be stable with respect to basis set used, exhibiting none of the artifacts of earlier basis-space partitioning methods. We also investigate the transferability of the resulting local polarizability models.
Collapse
Affiliation(s)
- Alston J Misquitta
- The University Chemical Laboratory, Lensfield Road, Cambridge CB2 1EW, United Kingdom.
| | | |
Collapse
|
470
|
Hill JG, Platts JA, Werner HJ. Calculation of intermolecular interactions in the benzene dimer using coupled-cluster and local electron correlation methods. Phys Chem Chem Phys 2006; 8:4072-8. [PMID: 17028695 DOI: 10.1039/b608623c] [Citation(s) in RCA: 181] [Impact Index Per Article: 10.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Potential energy curves for the parallel-displaced, T-shaped and sandwich structures of the benzene dimer are computed with density fitted local second-order Møller-Plesset perturbation theory (DF-LMP2) as well as with the spin-component scaled (SCS) variant of DF-LMP2. While DF-LMP2 strongly overestimates the dispersion interaction, in common with canonical MP2, the DF-SCS-LMP2 interaction energies are in excellent agreement with the best available literature values along the entire potential energy curves. The DF-SCS-LMP2 dissociation energies for the three structures are also compared with new complete basis set estimates of the interaction energies obtained from accurate coupled cluster (CCSD(T)) and DF-SCS-MP2 calculations. Since LMP2 is essentially free of basis set superposition errors, counterpoise corrections are not required. As a result, DF-SCS-LMP2 is computationally inexpensive and represents an attractive method for the study of larger pi-stacked systems such as truncated sections of DNA.
Collapse
Affiliation(s)
- J Grant Hill
- School of Chemistry, Cardiff University, Park Place, Cardiff, UK CF10 3AT
| | | | | |
Collapse
|
471
|
Antony J, Grimme S. Density functional theory including dispersion corrections for intermolecular interactions in a large benchmark set of biologically relevant molecules. Phys Chem Chem Phys 2006; 8:5287-93. [DOI: 10.1039/b612585a] [Citation(s) in RCA: 403] [Impact Index Per Article: 22.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
472
|
|
473
|
Parac M, Etinski M, Peric M, Grimme S. A Theoretical Investigation of the Geometries and Binding Energies of Molecular Tweezer and Clip Host−Guest Systems. J Chem Theory Comput 2005; 1:1110-8. [DOI: 10.1021/ct050122n] [Citation(s) in RCA: 60] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
474
|
Bukowski R, Podeszwa R, Szalewicz K. Efficient calculation of coupled Kohn–Sham dynamic susceptibility functions and dispersion energies with density fitting. Chem Phys Lett 2005. [DOI: 10.1016/j.cplett.2005.08.048] [Citation(s) in RCA: 98] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
475
|
Piacenza M, Grimme S. Van der Waals Complexes of Polar Aromatic Molecules: Unexpected Structures for Dimers of Azulene. J Am Chem Soc 2005; 127:14841-8. [PMID: 16231938 DOI: 10.1021/ja053613q] [Citation(s) in RCA: 87] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Full geometry optimizations at the dispersion corrected DFT-BLYP/TZV2P level of theory have been performed for dimers of azulene that may serve as a model system for the van der Waals complexes of polar pi systems. The structures and binding energies for 11 dimers are investigated in detail. The DFT-D interaction energies have been successfully checked against results from the accurate SCS-MP2/aug-cc-pVTZ approach. Out of the nine investigated stacked complexes, eight have binding energies larger than 7.4 kcal/mol (SCS-MP2) that exceed the value of 7.1 kcal/mol for the best naphthalene dimer. T-shaped arrangements (CH...pi) are significantly less stable. Two out of the three best structures have an antiparallel alignment of the monomer dipole moments in the complex, although the best ones with a parallel orientation are only about 0.5 kcal/mol less strongly bound which points to a minor importance of dipole-dipole interactions to binding. Quite surprisingly, the energetically lowest structure (DeltaE = -9.2 kcal/mol) corresponds to a situation where the two seven-membered rings are located almost on top of each other (7-7) and the long molecular axes are rotated against each other by 130 degrees. The 7-7 structural motif is found also in other energetically low-lying structures, and the expected 5-7 (two-side) arrangement is less strongly bound by about 2 kcal/mol. This can be explained by the electrostatic potential of azulene that only partially reflects the charge separation according to the common 4n + 2 pi electron rule. General rules for predicting stable van der Waals complexes of polar pi systems are discussed.
Collapse
Affiliation(s)
- Manuel Piacenza
- Theoretische Organische Chemie, Organisch-Chemisches Institut der Universität Münster, Corrensstrasse 40, D-48149 Münster, Germany
| | | |
Collapse
|
476
|
Podeszwa R, Szalewicz K. Accurate interaction energies for argon, krypton, and benzene dimers from perturbation theory based on the Kohn–Sham model. Chem Phys Lett 2005. [DOI: 10.1016/j.cplett.2005.07.029] [Citation(s) in RCA: 89] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
477
|
Hesselmann A. Local correlation potentials from Brueckner coupled-cluster theory. J Chem Phys 2005; 122:244108. [PMID: 16035747 DOI: 10.1063/1.1947167] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Local correlation potentials have been obtained from the nonlocal Brueckner coupled-cluster correlation potentials for the rare-gas atoms He, Ne, and Ar and the CO molecule. It is shown that the local correlation potential can mainly be expressed as a sum of two components: a "pure" correlation part and a relaxation contribution. While the total correlation potentials show an oscillating behavior near the nuclei, indicating the atomic shell structure, their components decrease rather monotonously, with a step structure in case of Ne and Ar. By looking at the determinantal overlap and one-electron properties it has been found that the orbitals obtained from these local potentials form a determinant which very well corresponds with the Brueckner determinant. Thus the previously found closeness between the Hartree-Fock determinant and the exchange-only Kohn-Sham determinant [Della Sala and Gorling, J. Chem. Phys. 115, 5718 (2001)] is confirmed also for the correlated case.
Collapse
Affiliation(s)
- A Hesselmann
- Center of Theoretical Chemistry, Arhus University, Denmark.
| |
Collapse
|