LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Mailing List Archives
Re: [lammps-users] how to introduce a repulsive term to Buckingham potential?
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [lammps-users] how to introduce a repulsive term to Buckingham potential?

From: Axel Kohlmeyer <akohlmey@...24...>
Date: Mon, 16 Jul 2018 12:27:45 -0400

On Mon, Jul 16, 2018 at 12:03 PM 河南师范大学王冬博 <417395440@...1204...> wrote:
Dear all,
         Recently, i am studying the silica system. Many articles refer to the BKS potential: U = Aexp(r/rho)-C/r**6, and some introduce a repulsive term to avoid the diverging at small interatomic distance, i think the potential like this: U = Aexp(r/rho)-(C/r**6)*(D/r**12). I want to konw which pair_style can describe this modified BKS potential?

​you have three options:
- write a custom pair style in C++ based on the existing buckingham potential
- use pair style table (please search the mailing list archives, suitable tables have been posted to the list multiple times)
- use pair style python to implement the interaction and then use pair_write to generate suitable tables (unless you don't care about performance).


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

Dr. Axel Kohlmeyer  akohlmey@...24...
College of Science & Technology, Temple University, Philadelphia PA, USA
International Centre for Theoretical Physics, Trieste. Italy.