One of the things still missing in the EFMO method is an interface to PCM. Here I attempt to sketch the method, based on the FMO-PCM and EFP-PCM interfaces.
The PCM electrostatic interaction free energy between solute and solvent is
For EFMO $\mathbf{V}$ is the electrostatic potential from the static multipoles and induced dipoles
The simplest approximation to $\mathbf{V}$ is
$$\begin{align*}
\mathbf{V}& =\sum_I^N \mathbf{V}_{I}\\
& = \sum_I^N (\mathbf{V}^{\text{mul}}_{I}+\mathbf{V}^{\mu}_{I})
\end{align*}$$
The potential at tesserae $j$ due to induced dipoles on fragment $I$ is given by:
$$\mathbf{V}_I^\mu (j)=\sum_{i \in I} (\mathbf{R}^T)_{ji}\boldsymbol{\mu_{i}}$$
The induced dipoles are obtained iteratively:
$$\boldsymbol{\mu}_i=\boldsymbol{\alpha}_i\left(\mathbf{F}^{\text{mul}}_i+\mathbf{F}^q_i-\sum_{i\neq j}\mathbf{D}_{ij}\boldsymbol{\mu}_j\right)$$
where $\boldsymbol{\alpha}_i$ is the dipole polarizability tensor at site $i$ and $\mathbf{F}^{\text{mul}}_i$ and $\mathbf{F}^q_i$ are the electrostatic fields due to all static multipoles and ASCs felt at point $i$.
Procedure:
1. Compute EFMO gas phase energy
2. Use gas phase static multipoles and $\boldsymbol{\mu}$ to construct $\mathbf{V}$
3. Solve $\mathbf{Cq=-V}$
This essentially means that the gas phase static multipoles and $\alpha$'s are corrected in step 2. E.g. for static monopoles ($q$)'s:
$$V_I^q(i)=\frac{q_i^{I}}{|r-r_i|}$$
and
$$\begin{align*}
V(i)&=V_I^q(i)+\sum^N_{J<I} (V_{IJ}^q(i)-V_I^q(i))\\
& = [q^{I}_i+\sum^N_{J<I}(q^{IJ}_i-q^{I}_i)]\frac{1}{|r-r_i|} \\
\end{align*}$$
All other steps are the same.
This work is licensed under a Creative Commons Attribution 4.0
The PCM electrostatic interaction free energy between solute and solvent is
$G_s=\frac{1}{2}\mathbf{V}^T \mathbf{q}$
$\mathbf{q}$ are the apparent surface charges (ASCs), which for large systems are obtained by solving this equation iteratively
$\mathbf{Cq=-V}$
For EFMO $\mathbf{V}$ is the electrostatic potential from the static multipoles and induced dipoles
$\mathbf{V=V^{\text{mul}}+V^{\mu}}$
The simplest approximation to $\mathbf{V}$ is
$$\begin{align*}
\mathbf{V}& =\sum_I^N \mathbf{V}_{I}\\
& = \sum_I^N (\mathbf{V}^{\text{mul}}_{I}+\mathbf{V}^{\mu}_{I})
\end{align*}$$
The potential at tesserae $j$ due to induced dipoles on fragment $I$ is given by:
$$\mathbf{V}_I^\mu (j)=\sum_{i \in I} (\mathbf{R}^T)_{ji}\boldsymbol{\mu_{i}}$$
The induced dipoles are obtained iteratively:
$$\boldsymbol{\mu}_i=\boldsymbol{\alpha}_i\left(\mathbf{F}^{\text{mul}}_i+\mathbf{F}^q_i-\sum_{i\neq j}\mathbf{D}_{ij}\boldsymbol{\mu}_j\right)$$
where $\boldsymbol{\alpha}_i$ is the dipole polarizability tensor at site $i$ and $\mathbf{F}^{\text{mul}}_i$ and $\mathbf{F}^q_i$ are the electrostatic fields due to all static multipoles and ASCs felt at point $i$.
Procedure:
1. Compute EFMO gas phase energy
2. Use gas phase static multipoles and $\boldsymbol{\mu}$ to construct $\mathbf{V}$
3. Solve $\mathbf{Cq=-V}$
4. Use $\mathbf{q}$ and $\boldsymbol{\mu}_i=\boldsymbol{\alpha}_i\left(\mathbf{F}^{\text{mul}}_i+\mathbf{F}^q_i-\sum_{i\neq j}\mathbf{D}_{ij}\boldsymbol{\mu}_j\right)$ to find new $\boldsymbol{\mu}$
5. Repeat steps 2-4 until self consistency
6. Compute $G_s$
For more accurate results one can approximate $\mathbf{V}$ as
$$\mathbf{V}=\sum_I^N \mathbf{V}_{I}+\sum^N_{I}\sum^N_{J<I} (\mathbf{V}_{IJ}-\mathbf{V}_{I}-\mathbf{V}_{J})$$This essentially means that the gas phase static multipoles and $\alpha$'s are corrected in step 2. E.g. for static monopoles ($q$)'s:
$$V_I^q(i)=\frac{q_i^{I}}{|r-r_i|}$$
and
$$\begin{align*}
V(i)&=V_I^q(i)+\sum^N_{J<I} (V_{IJ}^q(i)-V_I^q(i))\\
& = [q^{I}_i+\sum^N_{J<I}(q^{IJ}_i-q^{I}_i)]\frac{1}{|r-r_i|} \\
\end{align*}$$
All other steps are the same.
This work is licensed under a Creative Commons Attribution 4.0
No comments:
Post a Comment