level 1
尴尬力学男
楼主
请教一个关于分子动力学的问题。
我想模拟在单壁碳纳米管中氦原子作用,C原子为1类型原子,He原子为2类型原子。
其中希望C C采用airebo势, C He 和 He He采用LJ/cut 势。
但是我运行的时候错误如下:
Incorrect args for pair coefficients(../pair_airebo.cpp:192).
本人分子动力学入门中,特来请教。
感激任何回复。
相关的input文件和data文件如下:
_______________________________________________________________________________
Input files:
variable simname index cswnt_he
log log.${simname}.txt
unitsmetal
boundary s s s
lattice fcc 1.0
timestep 0.0002
atom_modify map array
atom_stylefull
read_data${simname}.data
#atom type 1 C
#atom type 2 He
pair_style hybrid airebo 3.0 lj/cut 10
pair_coeff * * airebo CH.airebo C Null
pair_coeff 1 2 lj/cut 0.00165375 3.191
pair_coeff 2 2 lj/cut 0.00093800 2.633
neighbor 2.0 bin
neigh_modifydelay 10
_______________________________________________________________________________
Data files:
LAMMPS data file. msi2lmp v3.9.8 / 06 Oct 2016 / CGCMM for SWNT_He
25200 atoms
2 atom types
0.003312122 37.250312122 xlo xhi
0.293685400 37.540685400 ylo yhi
-0.022014429 245.929185571 zlo zhi
Masses
1 12.011150 # cp
2 4.002600 # he
Atoms # full
1 1 1 0 19.19045576 35.53785354 0.616519345 0 0 0 # cp
2 1 1 0 21.99976872 35.23551897 0.614444202 0 0 0 # cp
3 1 1 0 26.01861866 33.87070195 0.614559762 0 0 0 # cp
4 1 1 0 23.37425783 34.89458376 0.61399558 0 0 0 # cp
5 1 1 0 27.26880791 33.19805089 0.614558771 0 0 0 # cp
6 1 1 0 29.59596562 31.5748403 0.613086376 0 0 0 # cp
7 1 1 0 32.44801397 28.41996606 0.614755423 0 0 0 # cp
8 1 1 0 30.6289291 30.5988577 0.613605121 0 0 0 # cp
9 1 1 0 33.22275244 27.22843876 0.614977618 0 0 0 # cp
10 1 1 0 34.47710228 24.68289273 0.615899395 0 0 0 # cp
11 1 1 0 35.45917005 20.54740305 0.614980059 0 0 0 # cp
12 1 1 0 34.92700818 23.33331907 0.615855631 0 0 0 # cp
13 1 1 0 35.55347936 19.132148 0.614872746 0 0 0 # cp
14 1 1 0 35.41197355 16.30051037 0.614947948 0 0 0 # cp
15 1 1 0 34.30592923 12.19942918 0.614949394 0 0 0 # cp
16 1 1 0 35.16116906 14.90368659 0.61498736 0 0 0 # cp
17 1 1 0 33.70986141 10.91155188 0.614906372 0 0 0 # cp
18 1 1 0 32.20140024 8.511419879 0.614667969 0 0 0 # cp
2017年07月31日 09点07分
1
我想模拟在单壁碳纳米管中氦原子作用,C原子为1类型原子,He原子为2类型原子。
其中希望C C采用airebo势, C He 和 He He采用LJ/cut 势。
但是我运行的时候错误如下:
Incorrect args for pair coefficients(../pair_airebo.cpp:192).
本人分子动力学入门中,特来请教。
感激任何回复。
相关的input文件和data文件如下:
_______________________________________________________________________________
Input files:
variable simname index cswnt_he
log log.${simname}.txt
unitsmetal
boundary s s s
lattice fcc 1.0
timestep 0.0002
atom_modify map array
atom_stylefull
read_data${simname}.data
#atom type 1 C
#atom type 2 He
pair_style hybrid airebo 3.0 lj/cut 10
pair_coeff * * airebo CH.airebo C Null
pair_coeff 1 2 lj/cut 0.00165375 3.191
pair_coeff 2 2 lj/cut 0.00093800 2.633
neighbor 2.0 bin
neigh_modifydelay 10
_______________________________________________________________________________
Data files:
LAMMPS data file. msi2lmp v3.9.8 / 06 Oct 2016 / CGCMM for SWNT_He
25200 atoms
2 atom types
0.003312122 37.250312122 xlo xhi
0.293685400 37.540685400 ylo yhi
-0.022014429 245.929185571 zlo zhi
Masses
1 12.011150 # cp
2 4.002600 # he
Atoms # full
1 1 1 0 19.19045576 35.53785354 0.616519345 0 0 0 # cp
2 1 1 0 21.99976872 35.23551897 0.614444202 0 0 0 # cp
3 1 1 0 26.01861866 33.87070195 0.614559762 0 0 0 # cp
4 1 1 0 23.37425783 34.89458376 0.61399558 0 0 0 # cp
5 1 1 0 27.26880791 33.19805089 0.614558771 0 0 0 # cp
6 1 1 0 29.59596562 31.5748403 0.613086376 0 0 0 # cp
7 1 1 0 32.44801397 28.41996606 0.614755423 0 0 0 # cp
8 1 1 0 30.6289291 30.5988577 0.613605121 0 0 0 # cp
9 1 1 0 33.22275244 27.22843876 0.614977618 0 0 0 # cp
10 1 1 0 34.47710228 24.68289273 0.615899395 0 0 0 # cp
11 1 1 0 35.45917005 20.54740305 0.614980059 0 0 0 # cp
12 1 1 0 34.92700818 23.33331907 0.615855631 0 0 0 # cp
13 1 1 0 35.55347936 19.132148 0.614872746 0 0 0 # cp
14 1 1 0 35.41197355 16.30051037 0.614947948 0 0 0 # cp
15 1 1 0 34.30592923 12.19942918 0.614949394 0 0 0 # cp
16 1 1 0 35.16116906 14.90368659 0.61498736 0 0 0 # cp
17 1 1 0 33.70986141 10.91155188 0.614906372 0 0 0 # cp
18 1 1 0 32.20140024 8.511419879 0.614667969 0 0 0 # cp