Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Update 4.3ADMPQeqForce.md #150

Merged
merged 1 commit into from
Nov 9, 2023
Merged
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
4 changes: 2 additions & 2 deletions docs/user_guide/4.3ADMPQeqForce.md
Original file line number Diff line number Diff line change
Expand Up @@ -21,7 +21,7 @@ The total coulombic energy is composed of the followling four parts.
Here, the first term $`E_{pc}`$ is electrostatic energy of point-charge systems, which is typically done in Ewald/PME method. For isolated molecule/clusters, it can also be computed using cutoff scheme:

```math
E_{pc}\left(q_i\right) = \sum_{i<j}\frac{q_i q_j}{r_{ij}}
E_{pc}\left(q_i\right) = \sum_{i< j}\frac{q_i q_j}{r_{ij}}
```
Note that we currently only consider Ewad3D (so called Ewald3DC method), instead of rigorous Ewald2D.

Expand All @@ -32,7 +32,7 @@ The second term is the short-range damping term, which can be different for diff
* For example, for EQeq model (`JPCL 3, 2506`), the damping kernel is ($`K`$ is the dielectric constant):

```math
E_{sr}\left(q_i\right) = \sum_{i<j} f(r_{ij})\frac{q_i q_j}{r_{ij}} \\
E_{sr}\left(q_i\right) = \sum_{i< j} f(r_{ij})\frac{q_i q_j}{r_{ij}} \\
f(r) = - \exp\left[-\left(\frac{J_{ij}r}{K}\right)^2\right]\left(1-\frac{J_{ij}r}{K}+\left(\frac{J_{ij}r}{K}\right)^2\right) \\
J_{ij} = \sqrt{J_i J_j}
```
Expand Down
Loading