force-field.md

@@ -13,7 +13,7 @@ r_{ij} = \sqrt{(x_j-x_i)^2 + (y_j-y_i)^2 + (z_j-z_i)^2}
 Atomic species are distinguished between electrode atoms and electrolyte atoms depending on their charge distribution. Electrode species are characterized by a Gaussian charge distribution of amplitude $`Q_i`$ and width $`\eta_i^{-1}`$
 
 ```math
-\rho_i (\mathbf{r}) = Q_i \eta_i^3 \pi ^{-3/2} \exp(-\eta_i^2 |\mathbf{r} - \mathbf{r_i}|^2)
+\rho_i (\mathbf{r}) = Q_i \eta_i^3 \pi ^{-3/2} \exp(-\eta_i^2 |\mathbf{r} - \mathbf{r}_i|^2)
 ```
 
 The value of the charge $`Q_i`$ can be kept constant or can fluctuate under a constant potential constraint. If so, the charges on the electrode atoms are computed at each timestep (for more information see [here](constant-potential)).
@@ -23,7 +23,7 @@ The electrodes are always treated as rigid bodies and can move only by consideri
 On the other hand, electrolyte species have a point charge distribution 
 
 ```math
-\rho_i (\mathbf{r}) = q_i \delta (\mathbf{r} - \mathbf{r_i})
+\rho_i (\mathbf{r}) = q_i \delta (\mathbf{r} - \mathbf{r}_i)
 ```
 
 They are usually mobile and can be polarizable (in polarizable force fields) and deformable (in the aspherical ion model).
@@ -43,10 +43,21 @@ and of course we have
 \mathbf{F_{ij}} = - \mathbf{F_{ji}}
 ```
 
+* [Electrostatic interactions](#elec) 
+   - [Coulomb potential](#coulomb)
+
+<a name="elec">
+
 # Electrostatic interactions
 
+</a>
+
+<a name="coulomb">
+
 ## Coulomb potential
 
+</a>
+
 The Coulomb potential for a set of point charges and Gaussian charges $`\rho=\sum_i \rho_i`$ is given by
 ```math
 V[\rho]=\frac{1}{2}\int d^3 \mathbf{r} d^3 \mathbf{r}'\frac{\rho(\mathbf{r})\rho(\mathbf{r}')}{\mid \mathbf{r}-\mathbf{r}'\mid}