----- Original Message ----- From: Sebastian Waltz <[email protected]> Date: Tuesday, June 1, 2010 21:33 Subject: [gmx-users] Forces on atom which is shifted out of the minimum To: [email protected]
> Hi all together, > > I need to calculate the first derivative of the force or > the second derivative of the potential. Therefore I need > the force on an atom which is shifted out of the minimum > along the constrained bond direction for each time step. > There in the Gromacs code is the force actually calculated? > Does someone has any idea where to start? Since GROMACS doesn't use the second derivatives (except perhaps in some implementations of EM), it's not routinely calculated. If you want to calculate them, then you will need to look at modifying a whole pile of code (but first see if anything useful is done with the EM integrators). Bonded and non-bonded interactions are treated differently. Making the non-bonded code calculate and output these second derivatives with any speed will be serious work. My usual observation is that you'll be hopelessly lost doing such far-reaching changes without having gotten a debugger and stepped through a GROMACS run on a toy example. Mark -- gmx-users mailing list [email protected] http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [email protected]. Can't post? Read http://www.gromacs.org/mailing_lists/users.php
