【分子动力学】lammps模拟碳纳米管的弹性模量

大佬们，想请教各位一个问题，我用lammps模拟碳纳米管的弹性模量时所得结果比别人论文中结果大10倍左右。以下是我的模拟步骤和in文件。恳请大家百忙之中帮助确定我的错误所在，万分感激！
模拟步骤：
1、 在MS中进行（8，8）单壁碳纳米管的构建，力场选用cvff；
2、 导出为data文件；
3、 在lammps中按照data文件中的力场设置势函数，然后进行压缩，输出应力应变；
4、 在python中作出应力应变图，计算弹性模量。
#------------模型基本设置------------------------
units metal #单位制
dimension 3
atom_style full #原子类型
boundary p p p

#-------------设置常数-------------------------
variable T equal 300 #初始温度
variable dt equal 0.0001 #模拟步长

#-----------设置邻居参数-----------------------
neighbor 2.0 bin
neigh_modify every 1 delay 0 check yes

#--------------势函数设置---------------------
pair_style lj/cut/coul/long 10
bond_style harmonic
angle_style harmonic
dihedral_style harmonic
improper_style cvff
kspace_style pppm/cg 1.0e-5
read_data cnt8.cvff.data
kspace_modify gewald 0.4

#-------------能量最小化--------------------------
min_style sd
minimize 1.0e-15 1.0e-15 10000 10000

#--------------系综设置----------------------------
fix myber all temp/berendsen ${T} ${T} 0.01 #nve系综,berendsen方式控温
fix mynve all nve
dump mydmp all atom 500 nve.xyz #输出轨迹文件
timestep ${dt} #模拟步长
thermo 500 #热力学输出
run 10000 #弛豫10000步
#去掉全局控温
unfix myber
undump mydmp

#--------------------------应力应变----------------------------------
variable tmp equal “lz”
variable L0 equal ${tmp}
variable strain equal “(lz - v_L0)/v_L0”
variable stressz equal “-pzz/10000”

#--------------重新设置轨迹文件输出-----------------
dump mydmp all atom 1000 run.xyz

#-------------------------应力应变保存到文件---------------------------
fix def3 all print 100 “${strain} ${stressz}” file stress-strain-cvff.dat screen no

#-------------设置压缩，压缩速率为-0.01-------------------------
fix 1 all nve
fix 2 all deform 100 z erate -0.01 units box remap x
run 30000