LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Mailing List Archives
Re: [lammps-users] Inquiry on Epsilon/Sigma Units in Quartic Bonds
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [lammps-users] Inquiry on Epsilon/Sigma Units in Quartic Bonds

From: Steve Plimpton <sjplimp@...24...>
Date: Thu, 12 Oct 2017 07:08:10 -0600

I think that's all that it's been used for, i.e. FENE chains.

However, I don't see a reason we couldn't allow sigma/epsilon to be specified
(optionally) as part of the bond coeff command.  In which case you would
have to carefully choose all the other params, as it describes in the doc

When monomers with sigma = epsilon = 1.0 are used, the following choice of parameters gives a quartic potential that looks nearly like the FENE potential: K = 1200, B1 = -0.55, B2 = 0.25, Rc = 1.3, and U0 = 34.6878. Different parameters can be specified using the bond_coeff command, but you will need to choose them carefully so they form a suitable bond potential.

to end with something that worked.


On Tue, Oct 3, 2017 at 11:25 AM, MASATO KOIZUMI <mkoizumi@...3709...> wrote:
Dear LAMMPS Users,

Hello. This is Masato Koizumi, and I am using LAMMPS for my molecular dynamics simulations. 

I would like to ask a question in regards to the quartic bond potential. If the sigma and epsilon parameters in the WCA terms of the quartic potential are set intrinsically to 1.0, I was curious if the bond_style quartic command can only be used for non-dimensional LJ units. Thank you.


Masato Koizumi

Check out the vibrant tech community on one of the world's most
engaging tech sites,!
lammps-users mailing list