Re: [lammps-users] Question about using r-respa integrator
Axel Kohlmeyer <akohlmey@...24...>
Thu, 12 Apr 2018 04:47:34 -0400
On Thu, Apr 12, 2018 at 12:13 AM, Pengyu Huang <hughh.py@...24...> wrote:
>> On Wed, Apr 11, 2018 at 1:40 PM, Pengyu Huang <hughh.py@...24...> wrote:
>> > Thank you very much, Axel, for your suggestions.
>> > I thought by defaults the bond, angle, dihedral and improper forces are
>> > calculated at the innermost timestep, and the multi-cut off setups are
>> > for
>> > the pair forces calculation, as the documents says " If not specified,
>> > the
>> > defaults are that bond forces are computed at level 1 (innermost loop),
>> > angle forces are computed where bond forces are, dihedral forces are
>> > computed where angle forces are, improper forces are computed where
>> > dihedral
>> > forces are, pair forces are computed at the outermost level, and kspace
>> > forces are computed where pair forces are. " Thanks for pointing out
>> > this,
>> > I will double check the setup.
yes, but even though there are multiple options available, you do not
have to use them all.
for things that are "tricky" i prefer to not rely on defaults, but
make things explicit.
>> > Also thanks for your suggestions for the co2. I also want to model co2
>> > with
>> > flexible bonds and angle, so in this case I do not need a rigid
>> > integrator and using r-RESPA might help?
there is a problem, though.
an angle potential for a linear molecule has a problem for its
minimum, since the direction of the forces is not well defined. most
implementations will likely error out or produce invalid numbers.
have you researched whether there is an angle style that could be used
for this and whether this is available in LAMMPS?
>> > Regards,
>> > Pengyu
Dr. Axel Kohlmeyer akohlmey@...24... http://goo.gl/1wk0
College of Science & Technology, Temple University, Philadelphia PA, USA
International Centre for Theoretical Physics, Trieste. Italy.