Thank you for kind answer.
It will be really helpful comments for me!
I have another question.
When it comes to HNEMD and NEMD examples in GPUMD example folder,
I found that the number of atoms in xyz.in file and the compute_shc command of
HNEMD and NEMD are different each other.
NEMD -> compute_shc 2 250 1 1000 400.0 group 0 4
-> xyz.in 40400
HNEMD -> compute_shc 2 250 1 1000 400.0 group 0 0
-> xyz.in 24000
I wonder there are special reasons for that.
If the same logic can be applied to the large system I considered,
can I adjust 'compute_shc' command to 'group 0 0' in HNEMD to save
--------- 원본 메일 ---------보낸사람 : Bruce Fan <brucenju@xxxxxxxxx>
받는사람 : <gpumd@xxxxxxxxxxxxx>
받은날짜 : 2020/Oct/15(Thu) 18:02:25
제목 : [gpumd] Re: Large size NEMD simulation.
1. Your results look reasonable. Within 10 nanosecond, your huge system is far
from reaching a steady state. It may require 100 nanosecond.
2. In the beginning, the rates of energy change in source and sink regions can
be different because the system is still in a transient state. The sink part
has larger rate of energy chage, because the local thermal conductivity (or
diffusivity) is larger there due to lower temperature.
3. You are calculating the SHC for group 0, which is fixed (freezed). This is
not a wise choice. You can make your groups smaller and number of groups
larger, and calculate SHC for any group other than group 0. You can minimize
the size of group 0, to minimize the waste of computing time.
4. Your other 8 groups are of the same size, with 4 micron for each. This means
you have 8 microns for source and sink regions and only 24 microns for the
"sample". You can make the source and sinl3 regions smaller to make your
5. You set the cutoff to 3 A and maximum neighbor to 10. This just wates memory
and time. Change them to 2.1 A and 3, repectively!
6. You are brave to consider this huge size. Good luck!
"소순성" <soonsung2001@xxxxxxxxxx> 于 2020年10月14日周三 07:13写道：
Thank you Zheyong Fan for using this good GPUMD program.
I have questions of this week.
I tested the NEMD calculation of very large size of graphene sheet,
(20nm x 30micro, about 22 million C atoms, x,z free boundary and y periodic
boundary, and biased in y direction.)
and found that the degree of energy of the thermostat coupling to the heat sink
region which is increasing is different with that to the sink region.
The absolute values of the slopes of the lines are different although total
energy seems stable.
It is my run.in file and I attached the part of information (It is continuing
because of very large size calculation!)
potential /GPUMD-2.5.1/potentials/tersoff/Graphene_Lindsay_2010_modified.txt 0
ensemble nvt_ber 300 300 0.005
ensemble heat_lan 300 500 100 1 8
compute 0 100 10 temperature
compute_shc 2 250 1 1000 500.0 group 0 0
Number of atoms is 22816728.
Maximum number of neighbors is 12.
Initial cutoff for neighbor list is 3 A.
Use orthogonal box.
Do not specify initial velocities here.
Have 1 grouping method(s).
Box lengths are
Lx = 1.9926000000e+02 A
Ly = 2.9999772000e+05 A
Lz = 3.3540000000e+00 A
Use free boundary conditions along x.
Use periodic boundary conditions along y.
Use free boundary conditions along z.
There are 9 groups of atoms in grouping method 0.
136728 atoms in group 0.
2835000 atoms in group 1.
2835000 atoms in group 2.
2835000 atoms in group 3.
2835000 atoms in group 4.
2835000 atoms in group 5.
2835000 atoms in group 6.
2835000 atoms in group 7.
2835000 atoms in group 8.
There is only one atom type.
22816728 atoms of type 0.
How can I interpret it? (Is it right simulation?)
Did the number of atoms in group 0 affect to the result?