用pair_style dpd/tstat

好的，这几天来我这看到了这个命令，但还是一直不确定，因为之前控温都是用的fix命令，而且复现的文献中说采用的LJ势，是不是可以用 pair_style hybrid  lj/cut  ……  dpd/tstat ……命令

引用回帖:

4楼: Originally posted by luyaobao at 2020-11-15 07:37:59
是
...

谢谢，那对应于pair_style dpd/tstat的pair_coeff 参数是什么呢，我看手册说分别是A (force units)、γ(force/velocity units)、cutoff (distance units)，其中的γ是阻尼系数吗，这个数应该怎么得到呢，我做的是氩原子的。

引用回帖:

6楼: Originally posted by luyaobao at 2020-11-16 12:37:35
温控没有系数A，γ不同液体粘度不同，看你的体系是什么了
...

units                  metal
dimension                     3
boundary                     f f p

atom_style                   atomic
neighbor                     2.0 bin
neigh_modify                 every 1  delay 5 check yes

region                             box block 0 169 0 168 0 400
create_box                   2 box

region                       wall_left block 0 4 0 168 0 400
region                       wall_right block 164 169 0 168 0 400
region                 wall_down block 4 164 0 4 0 400
region                 flow block 4 164 4 44 0 400

lattice                             fcc  5.738
create_atoms                 1 region flow

lattice                             fcc  3.483
create_atoms                 1 region wall_left
create_atoms                 1 region wall_right
create_atoms                 1 region wall_down

mass                               1 39.948
mass                               2 39.948

pair_style                   hybrid  lj/cut 8.5  dpd/tstat 86 86 8.5 34387
pair_coeff                   1 1 lj/cut 0.0103125 3.4
pair_coeff                   1 2 lj/cut 0.00515625 3.4
pair_coeff                   2 2 lj/cut 0.00515625 3.4
pair_coeff                   1 1 dpd/tstat 0.0672 8.5
comm_modify            vel yes

group                         1 region wall_left
group                         2 region wall_right
group                         3 region wall_down
group                         wall union 1 2 3
group                  flow region flow
group                         Ar subtract all wall
set                           group wall type 2

velocity               flow create 86 482748 mom yes rot yes
fix                           1 all nvt temp 86 86 0.02
compute                1 Ar chunk/atom bin/2d x lower 2 y lower 2 units box

velocity               wall set 0.0 0.0 0.0
fix                           2 wall setforce 0.0 0.0 0.0
fix                    3 flow ave/chunk 100 5 500 1 density/mass file c_0.5.profile

timestep                     0.005
thermo                             2000
thermo_style           custom step temp pe press density etotal

dump                   1 all custom 1000 c_0.5.dump id type x y z
run                    310000


想问您以下为什么我这个in文件输出的结果从始至终原子都不动，我文件中也设置了velocity命令了，您知道吗
，

引用回帖:

7楼: Originally posted by TW九九 at 2020-11-20 13:34:05
units                  metal
dimension                     3
boundary                     f f p

atom_style                   atomic
neighbor                     2.0 bin
neigh_modify                 every 1  delay 5 check yes
...

写错了，我的in文件没有fix nvt那一条命令，然后原子始终不动，后来我试着加上nvt，又动的太快了，几千步就原子丢失了