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 "sample" longer.
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!)
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?