[gpumd] Re: How can I set the cutoff and maximum number of neighbors in xyz.in for graphite systems?

  • From: Bruce Fan <brucenju@xxxxxxxxx>
  • To: gpumd@xxxxxxxxxxxxx
  • Date: Thu, 29 Apr 2021 16:47:15 +0800

Setting the maximum number of neighbors to *600 *is a good choice (the
maximum number allowed is *1024 *and you can even use this), but setting
the cutoff to 15 A is not good (it will slow down your simulation
significantly) if your LJ potential has a cutoff of 9 A. I have already
suggested a value: 9 A + 1 A = 10 A.

Zheyong

On Thu, Apr 29, 2021 at 3:52 PM "소순성" <soonsung2001@xxxxxxxxxx> wrote:

Hi,

Sorry for answering now the questions I asked long before .


1. maximum number of neighbors  & 2. cutoff

I set up the *maximum number of neighbors and cutoff *one by one (for
rule of thumb).

It is time-consuming process but for safer simulation, it seems to be
necessary.


With my xyz.in which have the coordinates of

74.03172420 0.00000000 0.00000000)

-37.01586210 64.11335384 0.00000000

0.00000000 0.00000000 6.70800000


I set the *maximum number of neighbors=600 and **cutoff=15. *and it
worked.



--------- Original Mail ---------

*Sender* : Bruce Fan <brucenju@xxxxxxxxx>
*Recipient* : <gpumd@xxxxxxxxxxxxx>
*Received Date* : 2021/Feb/26(Fri) 22:38:03
*Subject* : [gpumd] Re: How can I set the cutoff and maximum number of
neighbors in xyz.in for graphite systems?
1. maximum number of neighbors
In the manual, it means literally *maximum number of neighbors* any atom
can ever has.
As I calculated for graphene monolayer, I set it as 3 because just an atom
meets three adjacent atoms.
However, when I must consider the VdW interaction between each layer, how
can I estimate it? Just Do I test it?

*If you don't know a safe value for this parameter, just set it to the
maximal number that is allowed in GPUMD: 1024. This is a quite safe number.
When you simulate a very large system, say, a system with one million atoms
or more, this number is crucial for determining the memory consumption of
the code. In short: if this number is too large, it might waste memory; if
this number is too small, it might lead to unexpected errors. *

2. cutoff
In the manual, it means *initial cutoff distance used for building the
neighbor list *
and said "if you really don't know how to choose the value of "cutoff",
you can simply set it to cutoff distance of the maximal force cutoff in
your system plus 1.0 A."
Is it mean that if I use lj potential for graphene-graphene system
of which cutoff distance is about 9.0 A, I can set the "cutoff" in xyz.in
to 10.0 A?

*Yes, this is a very good choice.*

*Reminder: you may try to enable neighbor list updating to make your
simulation "safer", but one of our recent studies revealed an artifact for
heat transport calculations in these layered systems (with weak LJ
inter-layer interactions) that can be alleviated by disabling the neighbor
list updating! Sounds counterintuitive? You can check this paper: *

*Thermal conductivity reduction in carbon nanotube by fullerene
encapsulation: A molecular dynamics study*
*https://www.sciencedirect.com/science/article/abs/pii/S0008622320301378?via%3Dihub
<https://www.sciencedirect.com/science/article/abs/pii/S0008622320301378?via%3Dihub>*

How can I get the same results with triclinic system?

*There are no special things about triclinic systems compared to
rectangular ones regarding thermal transport calculations. The only caveat
in GPUMD is that one CANNOT use pressure control (i.e., NPT ensemble) when
using a triclinic box. This is a restriction for the time being and we may
try to improve it in the future. *


Zheyong


On Fri, Feb 26, 2021 at 10:19 AM "소순성" <soonsung2001@xxxxxxxxxx> wrote:

Hi, Fan.
Thank you for GPUMD and today, I have questions.

I'd like to HNEMD simulation for graphene bilayer or graphite multilayer.

So I prepared the xyz.in file with triclinic system and
here I have a question when it comes to what I set in the xyz.in file:

1. maximum number of neighbors
In the manual, it means literally *maximum number of neighbors* any atom
can ever has.
As I calculated for graphene monolayer, I set it as 3 because just an
atom meets three adjacent atoms.
However, when I must consider the VdW interaction between each layer, how
can I estimate it? Just Do I test it?


2. cutoff
In the manual, it means *initial cutoff distance used for building the
neighbor list *
and said "if you really don't know how to choose the value of "cutoff",
you can simply set it to cutoff distance of the maximal force cutoff in
your system plus 1.0 A."
Is it mean that if I use lj potential for graphene-graphene system
of which cutoff distance is about 9.0 A, I can set the "cutoff" in xyz.in
to 10.0 A?




How can I get the same results with triclinic system?

P.S I attached my input files for triclinic system.
















Other related posts: