X

Search Results

Searching....

11.5 Effective Fragment Potential Method

11.5.2 Theoretical Background

(July 14, 2022)

The total energy of the system consists of the interaction energy of the effective fragments (Eef-ef) and the energy of the ab initio (i.e., QM) region in the field of the fragments. The former includes electrostatics, polarization, dispersion and exchange-repulsion contributions (the charge transfer term, which might be important for description of the ionic and highly polar species, is omitted in the current implementation):

Eef-ef=Eelec+Epol+Edisp+Eex-rep. (11.62)

The QM-EF interactions are computed as follows. The electrostatics (Coulomb) and polarization parts of the EFP potential contribute to the quantum Hamiltonian via one-electron terms,

Hpq=Hpq+p|V^Coul+V^pol|q (11.63)

whereas dispersion and exchange-repulsion QM-EF interactions are treated as additive corrections to the total energy.

The electrostatic component of the EFP energy accounts for Coulomb interactions. In molecular systems with hydrogen bonds or polar molecules, this is the leading contribution to the total intermolecular interaction energy. 144 Buckingham A. D.
Q. Rev. Chem. Soc.
(1959), 13, pp. 183.
Link
An accurate representation of the electrostatic potential is achieved by using multipole expansion (obtained from the Stone’s distributed multipole analysis) around atomic centers and bond midpoints (i.e., the points with high electronic density) and truncating this expansion at octupoles. 1153 Stone A. J.
Chem. Phys. Lett.
(1981), 83, pp. 233.
Link
, 1152 Stone A. J., Alderton M.
Mol. Phys.
(1985), 56, pp. 1047.
Link
, 267 Day P. N. et al.
J. Chem. Phys.
(1996), 105, pp. 1968.
Link
, 411 Gordon M. S. et al.
J. Phys. Chem. A
(2001), 105, pp. 293.
Link
The fragment-fragment electrostatic interactions consist of charge–charge, charge–dipole, charge–quadrupole, charge–octupole, dipole–dipole, dipole–quadrupole, and quadrupole–quadrupole terms, as well as terms describing interactions of electronic multipoles with the nuclei and nuclear repulsion energy.

The multipole representation of the electrostatic density of a fragment breaks down when the fragments are too close. The multipole interactions become too repulsive due to significant overlap of the electronic densities and the charge-penetration effect. The magnitude of the charge-penetration effect is usually around 15% of the total electrostatic energy in polar systems, however, it can be as large as 200% in systems with weak electrostatic interactions. 1109 Slipchenko L. V., Gordon M. S.
J. Comput. Chem.
(2007), 28, pp. 276.
Link
To account for the charge-penetration effect, the simple exponential damping of the charge-charge term is used. 349 Freitag M. A. et al.
J. Chem. Phys.
(2000), 112, pp. 7300.
Link
, 1109 Slipchenko L. V., Gordon M. S.
J. Comput. Chem.
(2007), 28, pp. 276.
Link
The charge-charge screened energy between the expansion points k and l is given by the following expression, where αk and αl are the damping parameters associated with the corresponding expansion points:

Eklch-ch={[1-(1+αkRkl/2)e-αkRkl]qkqlRklif αk=αl[1-(αl2αl2-αk2)e-αkRkl-(αk2αk2-αl2)e-αlRkl]qkqlRklif αkαl. (11.64)

Alternatively, one can obtain the short-range charge-penetration energy using the spherical Gaussian overlap (SGO) approximation: 1110 Slipchenko L. V., Gordon M. S.
Mol. Phys.
(2009), 107, pp. 999.
Link

Eklpen=-2(1-2ln|Skl|)1/2Skl2Rkl (11.65)

where Skl is the overlap integral between localized MOs k and l, calculated for the exchange-repulsion term, Eq. (11.77). This charge-penetration energy is calculated and printed separately from the rest of the electrostatic energy. Using overlap-based damping generally results in a more balanced description of intermolecular interactions and is recommended, except possibly where short-range repulsive interactions are of interest. 168 Carter-Fenk K., Herbert J. M.
Mol. Phys.
(2022), 121, pp. e2055504.
Link

Electrostatic interaction between an effective fragment and the QM part is described by perturbation V^Coul of the ab initio Hamiltonian [see Eq.  (11.63)]. The perturbation enters the one-electron part of the Hamiltonian as a sum of contributions from the expansion points of the effective fragments. Contribution from each expansion point consists of four terms originating from the electrostatic potential of the corresponding multipole (charge, dipole, quadrupole, and octupole):

VkCoul(x)=qkT(rkx)-ax,y,zμakTa(rkx)+13abx,y,zΘabkTab(rkx)-115abcx,y,zΩabckTabc(rkx), (11.66)

where q, μ, Θ and Ω are the net charge, dipole, quadrupole and octupole, respectively, located at multipole expansion points k (all atoms and bond midpoints). T are the electrostatic tensors of ranks zero to three. Interaction of the QM electronic density with multipole charges can be also augmented by gaussian-type damping function, 430 Gurunathan P. K. et al.
J. Phys. Chem. B
(2016), 120, pp. 6562.
Link
, 1034 Rojas C. I. V., Slipchenko L. V.
J. Chem. Theory Comput.
(2020), 16, pp. 6408.
Link
such that Eq. (11.66) becomes

VkCoul(x)=[qknuc+qkele(1-exp(-βkrkx2)]T(rkx)-ax,y,zμakTa(rkx)+13abx,y,zΘabkTab(rkx)-115abcx,y,zΩabckTabc(rkx), (11.67)

where qknuc is the nuclear charge and qkele is the electronic charge on multipole point k, respectively. Thus, only the electronic charges are damped (smeared) by gaussians. Damping parameters α and β in Eqs. (11.64) and (11.67) are determined by minimizing the difference between the electrostatic potentials from the damped multipole expansion and the electronic wave function in the parameter-generating calculation for each fragment. (See Section 11.5.7.)

Polarization accounts for the intramolecular charge redistribution in response to external electric field. It is the major component of many-body interactions responsible for cooperative molecular behavior. EFP employs distributed polarizabilities placed at the centers of valence LMOs. Unlike the isotropic total molecular polarizability tensor, the distributed polarizability tensors are anisotropic.

The polarization energy of a system consisting of an ab initio and effective fragment regions is 267 Day P. N. et al.
J. Chem. Phys.
(1996), 105, pp. 1968.
Link

Epol=-12kax,y,zμak(Famult,k+Faai,nuc,k)+12kax,y,zμ¯akFaai,elec,k (11.68)

where μk and μ¯k are the induced dipole and the conjugated induced dipole at the distributed point k; Fmult,k is the external field due to static multipoles and nuclei of other fragments, and Fai,elec,k and Fai,nuc,k are the fields due to the electronic density and nuclei of the ab initio part, respectively.

The induced dipoles at each polarizability point k are computed as

μk=αkFtotal,k (11.69)

where αk is the distributed polarizability tensor at k. The total field Ftotal,k comprises from the static field and the field due to other induced dipoles, Fkind, as well as the field due to nuclei and electronic density of the ab initio region:

Fai,total,k=Fmult,k+Find,k+Fai,elec,k+Fai,nuc,k. (11.70)

As follows from the above equation, the induced dipoles on a particular fragment depend on the values of the induced dipoles of all other fragments. Moreover, the induced dipoles on the effective fragments depend on the ab initio electronic density, which, in turn, is affected by the field created by these induced dipoles through a one electron contribution to the Hamiltonian:

V^kpol(x)=12ax,y,z(μak+μ¯ak)Ta(rkx) (11.71)

where k are the polarizability expansion points. The total polarization contribution is computed self-consistently using a two level iterative procedure. The objectives of the higher and lower levels are to converge the wave function and induced dipoles for a given fixed wave function, respectively. In the absence of the ab initio region, the induced dipoles of the EF system are iterated until self-consistency with each other.

Self-consistent treatment of polarization accounts for many-body interaction effects. Polarization energy between EFP fragments is augmented by gaussian-like damping functions with default parameter α=β=0.6, applied to electric field F: 1110 Slipchenko L. V., Gordon M. S.
Mol. Phys.
(2009), 107, pp. 999.
Link

F=F0fdamp (11.72)
fdamp=1.0-(1+αβr2)e-(αβ)1/2r2 (11.73)

Dispersion provides a leading contribution to van der Waals and π-stacking interactions, 167 Carter-Fenk K., Herbert J. M.
Phys. Chem. Chem. Phys.
(2020), 22, pp. 24870.
Link
and its leading-order contribution is

Edisp=nC6R6 (11.74)

where coefficients C6 are derived from the frequency-dependent distributed polarizabilities with expansion points located at the LMO centroids, i.e., at the same centers as the polarization expansion points. The higher-order dispersion terms (induced dipole-induced quadrupole, induced quadrupole/induced quadrupole, etc.) are approximated as 1/3 of the C6 term. 21 Adamovic I., Gordon M. S.
Mol. Phys.
(2005), 103, pp. 379.
Link

For small distances between effective fragments, dispersion interactions are corrected for charge penetration and electronic density overlap effect either with the Tang-Toennies damping formula 1176 Tang K. T., Toennies J. P.
J. Chem. Phys.
(1984), 80, pp. 3726.
Link
with parameter b=1.5,

C6kl(1-e-bRk=06(bR)kk!)C6kl, (11.75)

or else using interfragment overlap (so-called overlap-based damping): 1110 Slipchenko L. V., Gordon M. S.
Mol. Phys.
(2009), 107, pp. 999.
Link

C6kl(1-Skl2(1-2log|Skl|+2log2|Skl|))C6kl. (11.76)

Note:  QM/EFP dispersion interactions are currently disabled.

Exchange-repulsion originates from the Pauli exclusion principle, which states that the wave function of two identical fermions must be anti-symmetric. In traditional classical force fields, exchange-repulsion is introduced as a positive (repulsive) term, e.g., R-12 in the Lennard-Jones potential. In contrast, EFP uses a wave function-based formalism to account for this inherently quantum-mechanical effect.

The exchange-repulsion interaction is derived as an expansion in the intermolecular overlap, truncated at the quadratic term, 554 Jensen J. H., Gordon M. S.
Mol. Phys.
(1996), 89, pp. 1313.
Link
, 555 Jensen J. H., Gordon M. S.
J. Chem. Phys.
(1998), 108, pp. 4772.
Link
which requires that each effective fragment carries a basis set that is used to calculate overlap and kinetic one-electron integrals for each interacting pair of fragments. The exchange-repulsion contribution from each pair of localized orbitals i and j belonging to fragments A and B, respectively, is:

Eijexch=-4-2ln|Sij|πSij2Rij-2Sij(kAFikASkj+lBFjlBSil-2Tij)+2Sij2(JB-ZJRiJ-1+2lBRil-1+IA-ZIRIj-1+2kARkj-1-Rij-1) (11.77)

where i, j, k and l are the LMOs, I and J are the nuclei, S and T are the intermolecular overlap and kinetic energy integrals, and F is the Fock matrix element.

The expression for the Eijexch involves overlap and kinetic energy integrals between pairs of localized orbitals. In addition, since Eq. (11.77) is derived within an infinite basis set approximation, it requires a reasonably large basis set to be accurate. [The 6-311++G(3df,2p) basis set is recommended and 6-31+G is considered to be the smallest acceptable basis set.] These factors make exchange-repulsion the most computationally expensive part of an EFP energy calculation for systems of moderate size.

Large systems require additional considerations. Since total exchange-repulsion energy is given by a sum of terms in Eq. (11.77) over all the fragment pairs, its computational cost formally scales as 𝒪(N2) with the number of effective fragments N. However, exchange-repulsion is a short-range interaction; the overlap and kinetic energy integrals decay exponentially with the inter-fragment distance. Therefore, by employing a distance-based screening, the number of overlap and kinetic energy integrals scales as 𝒪(N). Consequently, for large systems exchange-repulsion may become less computationally expensive than the long-range components of EFP (such as Coulomb interactions).

Note:  The QM/EFP exchange-repulsion energy is currently disabled.