LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Mailing List Archives
Re: [lammps-users] Fwd: Regarding Tersoff and Eam pot.
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [lammps-users] Fwd: Regarding Tersoff and Eam pot.


From: Axel Kohlmeyer <akohlmey@...24...>
Date: Mon, 21 Aug 2017 09:18:11 -0400

your input prints outlandishly bogus energies in the very first step.
since the input works with the single potential file, it cannot be the geometry, thus it must be the potential files or parameters.

the cutoff for morse is very strange. have you plotted the potential with these parameters? does this look reasonable to you?
have you tested the other parameter files individually?
as a quick hack, you can use:

pair_modify pair <substyle> compute no

to turn off the individual contributions.

at any rate, this looks like it is something that you should have easily caught yourself, had you done the proper diligent testing of each component of your input, rather than dumping it on the mailing list. i am particularly upset, that you didn't report this extremely unphysical energy (and pressure) at the very beginning.
that would have saved a *lot* of time and effort. it is also irritating, that even after been explicitly asked to follow the mailing list guide lines, you still don't adhere to guidelines 3., 8. and 9.

axel.

On Mon, Aug 21, 2017 at 4:37 AM, Sharma MD <sa.md.dtm@...24...> wrote:

---------- Forwarded message ----------
From: Sharma MD <sa.md.dtm@...24...>
Date: Sun, Aug 20, 2017 at 1:46 PM
Subject: Re: [lammps-users] Regarding Tersoff and Eam pot.
To: Axel Kohlmeyer <akohlmey@...24...>


Thank you for suggestion.
I used Be as Substrate material and Diamond as tool. Be is HCP and Diamond has cubic diamond structure. substrate has 3 layers of atoms: rigid layer at bottom, thermostat layer in the middle and newton layer at the top. Similarly for Diamond. rigid layer is for supporting the structure and thermo layer for diffusing the heat generated in the cutting zone and newton layer for actual dynamics of system. Thermo + Newton layer together forms mobile layer on which NVE dynamics is employed.
First the equilibration is carried out on this mobile layer at 293K for 10 ps before cutting starts. However, error occurs.

Files are attached. Kindly suggest.

On Sat, Aug 19, 2017 at 6:58 PM, Axel Kohlmeyer <akohlmey@...24...> wrote:


On Sat, Aug 19, 2017 at 9:12 AM, Sharma MD <sa.md.dtm@...24...> wrote:
It happened in the very first step.

System was fine.

there must be something bogus in your input that causes the overflow of the memory allocation during neighbor list setup.
but you are not providing sufficient information, so there is no way to track this down.

if you want to get constructive help, you have to follow the mailing list guidelines and provide sufficient information and detail
 

On Sat, Aug 19, 2017 at 5:35 PM, Axel Kohlmeyer <akohlmey@...24...> wrote:


--
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.



------------------------------------------------------------------------------
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot
_______________________________________________
lammps-users mailing list
lammps-users@...6297....sourceforge.net
https://lists.sourceforge.net/lists/listinfo/lammps-users




--
Dr. Axel Kohlmeyer  akohlmey@...12...24...  http://goo.gl/1wk0
College of Science & Technology, Temple University, Philadelphia PA, USA
International Centre for Theoretical Physics, Trieste. Italy.