24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

yjhqdm

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 752
散金: 15
红花: 1
帖子: 196
在线: 64.6小时
虫号: 289465
注册: 2006-10-22
专业: 物理学II

[求助] lammps新手超简单程序疑问已有3人参与

写了一个简单的Lammps程序，但是问题多多，实在无法解答，求助。
我上传了附件，求大侠下载试试！
文件在附件中，执行操作如下：
E:\lammps>lmp_win_no-mpi -in example.input
LAMMPS (12 Apr 2013)
Reading data file ...
  orthogonal box = (-2 -2 -2) to (2 2 2)
  1 by 1 by 1 MPI processor grid
ERROR: Unknown identifier in data file: 1 12.01 (../read_data.cpp:343)
提示在Masses的第一行时错了，好奇怪，求大侠解疑。。。新手。
我在data文件中写了Header了，还是错误。。。。哪位下载下来直接运行看看我哪里错了

Lammps程序为：
variable vz index 10 20
log log_${vz}.txt
units metal
dimension 3
boundary f p s
atom_style atomic
read_data example.data
pair_style tersoff
pair_coeff * * SiCGe.tersoff C C C
group a id == 1 3 5
group 1 id == 2 4 6
group a1 union a 1
neighbor 2.0 bin
neigh_modify delay 0
minimize 1.0e-6 1.0e-6 1000 1000
fix 1 all nve
thermo 100
compute 1 all temp
compute 2 all ke/atom
compute 3 all pe/atom
variable et atom c_2[]+c_3[]
fix 2 all ave/time 1 1 10 c_1 file temp_${vz}.dat
dump 1 all xyz 100 config_${vz}.dat
dump 2 all custom 100 etotal_${vz}.dat v_et
velocity all create 400.0 273237 units box
run 1000
clear
next vz
jump example.input

example.dat文件为：
LAMMPS data file

6 atoms

3 atoms types

-2 2 xlo xhi
-2 2 ylo yhi
-2 2 zlo zhi

Masses

1 12.01
2 12.01
3 12.01

Atoms

1 1  1.000000 0.000000 0.000000
2 2  0.500000 0.866025 0.000000
3 3 -0.500000 0.866025 0.000000
4 3 -1.000000 0.000000 0.000000
5 2 -0.500000 -0.866025 0.000000
6 1  0.500000 -0.866025 0.000000

以上data文件我写了Header了。

回复此楼