TY - JOUR
T1 - Incorporation of charge transfer into the explicit polarization fragment method by grand canonical density functional theory
AU - Isegawa, Miho
AU - Gao, Jiali
AU - Truhlar, Donald G.
N1 - Funding Information:
The authors are grateful to Hannah Leverentz for assistance with the calculations. This work is supported by the NIH (Gant No. NIGMS/1RC1GM091445) and NSF (Grant No. CHE09-56776).
PY - 2011/8/28
Y1 - 2011/8/28
N2 - Molecular fragmentation algorithms provide a powerful approach to extending electronic structure methods to very large systems. Here we present a method for including charge transfer between molecular fragments in the explicit polarization (X-Pol) fragment method for calculating potential energy surfaces. In the conventional X-Pol method, the total charge of each fragment is preserved, and charge transfer between fragments is not allowed. The description of charge transfer is made possible by treating each fragment as an open system with respect to the number of electrons. To achieve this, we applied Mermins finite temperature method to the X-Pol wave function. In the application of this method to X-Pol, the fragments are open systems that partially equilibrate their number of electrons through a quasithermodynamics electron reservoir. The number of electrons in a given fragment can take a fractional value, and the electrons of each fragment obey the Fermi-Dirac distribution. The equilibrium state for the electrons is determined by electronegativity equalization with conservation of the total number of electrons. The amount of charge transfer is controlled by re-interpreting the temperature parameter in the Fermi-Dirac distribution function as a coupling strength parameter. We determined this coupling parameter so as to reproduce the charge transfer energy obtained by block localized energy decomposition analysis. We apply the new method to ten systems, and we show that it can yield reasonable approximations to potential energy profiles, to charge transfer stabilization energies, and to the direction and amount of charge transferred.
AB - Molecular fragmentation algorithms provide a powerful approach to extending electronic structure methods to very large systems. Here we present a method for including charge transfer between molecular fragments in the explicit polarization (X-Pol) fragment method for calculating potential energy surfaces. In the conventional X-Pol method, the total charge of each fragment is preserved, and charge transfer between fragments is not allowed. The description of charge transfer is made possible by treating each fragment as an open system with respect to the number of electrons. To achieve this, we applied Mermins finite temperature method to the X-Pol wave function. In the application of this method to X-Pol, the fragments are open systems that partially equilibrate their number of electrons through a quasithermodynamics electron reservoir. The number of electrons in a given fragment can take a fractional value, and the electrons of each fragment obey the Fermi-Dirac distribution. The equilibrium state for the electrons is determined by electronegativity equalization with conservation of the total number of electrons. The amount of charge transfer is controlled by re-interpreting the temperature parameter in the Fermi-Dirac distribution function as a coupling strength parameter. We determined this coupling parameter so as to reproduce the charge transfer energy obtained by block localized energy decomposition analysis. We apply the new method to ten systems, and we show that it can yield reasonable approximations to potential energy profiles, to charge transfer stabilization energies, and to the direction and amount of charge transferred.
UR - http://www.scopus.com/inward/record.url?scp=80052402208&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=80052402208&partnerID=8YFLogxK
U2 - 10.1063/1.3624890
DO - 10.1063/1.3624890
M3 - Article
C2 - 21895159
AN - SCOPUS:80052402208
VL - 135
JO - Journal of Chemical Physics
JF - Journal of Chemical Physics
SN - 0021-9606
IS - 8
M1 - 084107
ER -