where FC is a conservative force, FD is a dissipative force, and FR is +
where Fc is a conservative force, Fd is a dissipative force, and Fr is a random force. Rij is a unit vector in the direction Ri - Rj, Vij is the vector difference in velocities of the two atoms = Vi - Vj, alpha is a Gaussian random number with zero mean and unit variance, dt is @@ -41,6 +41,9 @@ the timestep size, and w(r) is a weighting factor that varies between 0 and 1. Rc is the cutoff. Sigma is set equal to sqrt(2 T gamma), where T is a parameter in the pair_style command.
+The pairwise energy associated with this potential is only due to the +conservative force term Fc. +
The following coefficients must be defined for each pair of atoms types via the pair_coeff command as in the examples above, or in the data file or restart files read by the diff --git a/doc/pair_dpd.txt b/doc/pair_dpd.txt index 0509b988b8..1cb74e67e3 100644 --- a/doc/pair_dpd.txt +++ b/doc/pair_dpd.txt @@ -30,7 +30,7 @@ on atom I due to atom J is given as a sum of 3 terms :c,image(Eqs/pair_dpd.jpg) -where FC is a conservative force, FD is a dissipative force, and FR is +where Fc is a conservative force, Fd is a dissipative force, and Fr is a random force. Rij is a unit vector in the direction Ri - Rj, Vij is the vector difference in velocities of the two atoms = Vi - Vj, alpha is a Gaussian random number with zero mean and unit variance, dt is @@ -38,6 +38,9 @@ the timestep size, and w(r) is a weighting factor that varies between 0 and 1. Rc is the cutoff. Sigma is set equal to sqrt(2 T gamma), where T is a parameter in the pair_style command. +The pairwise energy associated with this potential is only due to the +conservative force term Fc. + The following coefficients must be defined for each pair of atoms types via the "pair_coeff"_pair_coeff.html command as in the examples above, or in the data file or restart files read by the