[Openmd-users] Update

Dan Gezelter gezelter at nd.edu
Tue Jul 27 09:57:01 EDT 2010


> 
>  You had me worried, so I just checked by printing out inside the LJ code that calculates the potential.  With cutoffRadius and switchingRadius set to the same value, and electrostaticSummationMethod set to "SHIFTED_FORCE", you do indeed get the shifted LJ potential.


Clarification:  You get the FORCE-shifted LJ potential.   V(r) = V_{LJ}(r) - V_{rcut} -  V'_{LJ}(rcut) * (r - rcut)

 --D

***********************************************
  J. Daniel Gezelter
  Associate Professor of Chemistry
  Director of Graduate Admissions
  Department of Chemistry and Biochemistry
  251 Nieuwland Science Hall
  University of Notre Dame
  Notre Dame, IN 46556-5670

  phone:  +1 (574) 631-7595
  fax:    +1 (574) 631-6652
  e-mail: gezelter at nd.edu
  web:    http://www.nd.edu/~gezelter
************************************************





More information about the Openmd-users mailing list