1 / 65

A Density-Functional Theory for Covalent and Noncovalent Chemistry

A Density-Functional Theory for Covalent and Noncovalent Chemistry. Non-empirical and fast . • Review of the XDM (exchange-hole dipole moment) dispersion model of Becke and Johnson • Combining the model with non-empirical exchange-correlation GGAs (Kannemann and Becke)

aria
Download Presentation

A Density-Functional Theory for Covalent and Noncovalent Chemistry

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. A Density-Functional Theory for Covalent and Noncovalent Chemistry Non-empirical and fast  • • Review of the XDM (exchange-hole dipole moment) dispersion model of Becke and Johnson • • Combining the model with non-empirical exchange-correlation GGAs (Kannemann and Becke) • • Tests on standard bio-organic benchmark sets

  2. Dispersion Interactions from the Exchange-Hole Dipole Moment: The “XDM” model Axel D. Becke and Erin R. Johnson* Email: axel.becke@dal.ca Department of Chemistry, Dalhousie University, Halifax, Nova Scotia, Canada (*now at University of California, Merced)

  3. What is the source of the “instantaneous” multipole moments that generate the dispersion interaction? Suggestion: Becke and Johnson, J. Chem. Phys. 122, 154104 (2005) > Becke and Johnson, J. Chem. Phys. 127, 154108 (2007) < The dipole momentof the exchange hole! (Position, rather than time, dependent)

  4. In Hartree-Fock theory, the total energy of many-electron system is given by: The “exchange” energy

  5. The Exchange Hole The ( -spin) exchange energy can be rewritten as follows: where is called the “exchange hole”. Physical interpretation:each electron interacts with a “hole” whose shape (in terms of r2) depends on the electron’s position r1. When an electron is at r1, the hole measures the depletion of probability, with respect to the total electron density, of finding another electron of the same spin at r2. This arises from exchange antisymmetry.

  6. It is simple to prove that… • • The hole is always negative. • • The probability of finding another same-spin electron at r2=r1 (“on top” of the reference electron) is completely extinguished: Pauli or exchange “repulsion”! • The hole always contains exactly (minus) one electron:

  7. Dispersion Model: the Basic Idea An electron plus its exchange hole always has zero total charge but in general a non-zero dipole moment! This r1-dependent dipole moment is easily obtained by integrating over r2: Note that only occupied orbitals are involved. (This reduces, for the H atom, to the exact dipole moment of the H atom when the electron is at r1.)

  8. e- (r, W) dXs h+ (r - dXs , W) Nucleus In a spherical atom, consider the following simplified “2-point” picture: Notice that this picture generates higher multipole moments as well(with respect to the nucleus as origin)given by and that all these moments depend only on the magnitudeof the exchange-holedipolemoment.

  9. This is significant because the magnitude dXs(r) of the exchange-hole dipole moment can be approximated using local densities and the Becke-Roussel exchange-hole model [Phys. Rev. A 39, 3761 (1989)], a 2nd-order GGA: where b is the displacement from the reference point of the mean position of the BR model hole [Becke and Johnson, J. Chem. Phys. 123, 154101 (2005)].

  10. Therefore, the entire van der Waals theory that follows has two variants: Orbital (XX)based, or Density-functional (BR)based XX performs slightly better in rare-gas systems. BR performs better in intermolecular complexes. All our current work employs the BR (DFT) variant.

  11. The Dispersion Interaction: Spherical Atoms rA Vint(rA,rB) rB Vint(rA,rB) = multipole moments of electron+hole at rAinteracting with multipole moments of electron+hole at rB

  12. 2nd-Order Ground-State Perturbation Theory in the Closure (Ünsold) Approximation If the first-order, ground-state energy correction arising from a perturbation Vpert is zero: Then the second-order correction is approximately given by where the expectation values are in the ground state and is the average excitation energy.

  13. To evaluate the expectation value <Vint2> square the multipole-multipole interaction Vint(rA,rB): Vint2(rA,rB) = (dipole-dipole + dipole-quadrupole + dipole-octopole + quadrupole-quadrupole +…)2 Then integrate the squared interaction over all rAand rB. This is a “semiclassical” calculation of <Vint2>.

  14. The result is where with atomic moment integrals given by

  15. What about ? Assume that and that, for each atom, where α is its dipole polarizability. This easily follows from the same “semiclassical” 2nd-order perturbation theory applied to the polarizability of each atom. Thus our C6, C8, C10’s depend on atomic polarizabilities and moment integrations from Hartree-Fock (or KS)calculations! No fitted parameters or explicitly correlated wavefunctions!

  16. How well does it work? On the 21 pairs of the atoms H, He, Ne, Ar, Kr, Xe, the mean absolute percent errors are: C6 3.4 % C8 21.5 % C10 21.5 %

  17. From Free Atoms to Atoms in Molecules • Partition a molecular system into “atoms” using Hirshfeld weight functions: where is a spherical free atomic density placed at the appropriate nucleus and the n summation is over all nuclei. has value close to 1 at points near nucleus i and close to 0 elsewhere. Also, • Assume that an intermolecularCm(m=6,8,10) can be written as a sum of interatomicCm,ij : i in A j in B

  18. In the previous expressions for Cm replace A and B with i and j : with generalized to

  19. and, where is the effective polarizability of atom i in A. We propose that These are effective volume integrations. This is motivated by the well known qualitative (if not quantitative) general relationship between polarizability and volume. See Kannemann and Becke, JCP 136, 034109 (2012). All radii r in the above integrals and in the integrals for are with respect to the position of nucleus i.

  20. Atom-Molecule and Molecule-Molecule Dispersion Coefficients Test set: H2 and N2 with He, Ne, Ar, Kr, and Xe. Cl2 with He, Ne, Ar, Kr, and Xe (except C10). H2-H2 H2-N2 N2-N2 (only H2-H2 for C10). Fully-numerical Hartree-Fock calculations on the monomers using the NUMOL program (Becke and Dickson, 1989). MAPEs with respect to dispersion coefficients from frequency dependent MBPT polarizabilities: 12.7% for C6 16.5% for C8 11.9% for C10 (On a much more extensive test set of 178 intermolecular C6’s, the model has a MAPE of 9.1%)

  21. Everything, so far, has been about the asymptotic dispersion series between atom pairs, The asymptotic series needs to be damped in order to avoid divergences when R is small. i.e. need information about characteristic R values inside of which the asymptotic series is no longer valid. The usual approach is to use empirical vdW radii. However…

  22. “Critical” Interatomic Separation Since we can compute C6, C8, and C10 non-empirically, we can obtain non-empirical range information. There is a “critical” Rc,ij where the three dispersion terms are approximately equal: Take Rc,ij as the average of , , and The asymptotic dispersion series is obviously meaningless insideRc,ij. Therefore …

  23. General Dispersion Energy Formula … we use Rc,ij to damp the dispersion energy at small internuclear separations as follows: where Rvdw,ij is an effective van der Waals separation with only two universal fit parameters.Best-fita1and a2 values depend on the exchange-correlation theories with which the above is combined. … which brings us to the second part of the talk …

  24. Key XDM Publications A. D. Becke and E. R. Johnson, J. Chem. Phys. 127, 154108 (2007) Dispersion coefficients E. R. Johnson and A. D. Becke, J. Chem. Phys. 124, 174104 (2006) Damping functions A. D. Becke and E. R. Johnson, J. Chem. Phys. 123, 154101 (2005) Transform to a DFT

  25. Exchange-Correlation GGAs XDM Dispersion + Axel D. Becke and Felix O. Kannemann Email: axel.becke@dal.ca Department of Chemistry, Dalhousie University, Halifax, Nova Scotia, Canada

  26. It would be nice if we could use “standard” XC-GGAs (instead of, eg., Hartree-Fock) plus the XDM dispersion model to treat vdW interactions. Let’s look at the exchangepart first.

  27. Exchange GGA Functionals B86, B86b, B88 where the local (spin) density part is and the “reduced” (spin) density gradient is

  28. Exchange GGA Functionals PW86, PW91, PBE(96) where The “reduced” (spin) density gradient is and, for spin polarized systems,

  29. revPBE (Yang group)

  30. Now plot the exchange-only interaction energy in Ne2 Note that Hartree-Fock (exact exchange) is repulsive! Plot is from Kannemann and Becke, JCTC 5, 719 (2009)

  31. Ne2

  32. Can get anything! from artifactual binding to massive over-repulsion, depending on the choice of functional! Which exchange GGA best reproduces exact Hartree-Fock? Lacks and Gordon, PRA 47, 4681 (1993) Kannemann and Becke, JCTC 5, 719 (2009) Murray, Lee, and Langreth, JCTC 5, 2754 (2009) PW86, followed by B86b

  33. Exchange Enhancement Factor [Zhang, Pan, Yang, JCP 107, 7921 (1997)]

  34. PW86 is a completelynon-empirical exchange functional! Its 4 parameters are fit to a theoretical exchange-hole model. Perdew and Wang, PRB 33, 8800 (1986) That PW86 accurately reproduces Hartree-Fock repulsion energies is remarkable. The underlying theoretical model (truncated “GEA” hole) knows nothing about closed-shell atomic or molecular interactions! Could be a fortuitous accident?! Nevertheless, no parameters need to be fit to data 

  35. Now, what about (dynamical) correlation? Use the non-empirical “PBE” correlation functional: Perdew, Burke, and Ernzerhof, PRL 77, 3865 (1996) Therefore we have, with a1 and a2 to be determined.

  36. How to determine a1 and a2? Fit to the binding energies of the prototypicaldispersion-bound rare-gas systems He2 HeNe HeAr Ne2 NeAr Ar2 (reference data from Tang and Toennies, JCP 118, 4976 (2003)) At the CBS limit, we find the best-fit values a1=0.65 a2=1.68 (Kannemann and Becke, to be published)

  37. Gaussian 09, aug-cc-pV5Z, counterpoise, ultrafine grid, BEs in microHartree a1 = 0.65 a2 = 1.68Å RMS%E = 4.2%

  38. Note that in all subsequent benchmarking there is no (re)fitting of parameters! Our functional is, from here on, essentially nonempirical in all its parts.

  39. GGA + XDM Publications F.O. Kannemann and A.D. Becke, JCTC 5, 719 (2009) Rare-gas diatomics (numerical post-LDA) F.O. Kannemann and A. D. Becke, JCTC 6, 1081 (2010) Intermolecular complexes (numerical post-LDA) A.D. Becke, A.A. Arabi and F.O.Kannemann, Can. J. Chem. 88, 1057 (2010) Dunning aDZ and aTZ basis-set calculations (post-G09)

  40. Tests on Standard BioOrganic Benchmark Sets Axel D. Becke and Felix O. Kannemann Email: axel.becke@dal.ca Department of Chemistry, Dalhousie University, Halifax, Nova Scotia, Canada

  41. The “S22” and “S66” vdW Benchmark Sets of Hobza et al S22 S66 hydrogen bonding dispersion other noncovalent interactions

  42. References S22: Jurecka, Sponer, Cerny, Hobza, PCCP 8, 1985 (2006) S66(x8): Rezac, Riley, Hobza JCTC 7, 2427 (2011) The next slides contain Mean Absolute Deviations (MADs) for the S22 and the S66 benchmark sets in comparison with other popular DFT methods. Data for all methods other than ours are from Goerigk, Kruse & Grimme, CPC 12, 3421 (2011) All computations employ the def2-QZVP basis set

  43. PW86PBE-XDM no CP (0.30) CP (0.28)

  44. S66 PW86PBE-XDM no CP (0.27) CP (0.23)

  45. Aside Can we do better by combining XDM with hybrid functionals? (rather than the pure GGA, PW86+PBE)

  46. Burns, Vazquez-Mayagoitia, Sumpter, Sherrill, JCP 134, 084107 (2011) S22 MAD (kcal/mol) for B3LYP-XDM and other DFT methods

  47. However • The pure GGA, PW86+PBE, is completely nonempirical (B3LYP, with 3 fitted parameters, is not) • Density-fit basis sets can speed up pure GGA calculations, with no loss of accuracy, by an order of magnitude!

  48. Tests on the S66x8 Set (Nonequilibrium Geometries) Axel D. Becke and Felix O. Kannemann Email: axel.becke@dal.ca Department of Chemistry, Dalhousie University, Halifax, Nova Scotia, Canada

More Related