求助：MS转data文件警告以及lammps运行报错

cdcdcd

各位老师好，我想请教一下关于将MS导出模型用msi2lmp小程序转data文件时出现的警告以及运行报错的问题
1.将ms中的模型在cvff力场导出为car文件和mdf文件，用msi2lmp程序转data文件时出现这样的警告
  WARNING inconsistent # of connects on atom %d type c
  请问这种这种警告会影响后面的计算吗？如果会产生影响,用什么方法进行修改？
2.运行时报错为ERROR on proc 36: Bond atoms 281 303 missing on proc 36 at step 1822728 (../ntopo_bond_all.cpp:63)
   我查看手册中的错误信息是这样的
   Bond atoms %d %d missing on proc %d at step %ld
   The second atom needed to compute a particular bond is missing on this processor. Typically this is because the pairwise cutoff is set too short or the bond has blown apart and an atom is too far away.
   修改了很多次截断半径，要么依然报这样的错，要么报错为Neighbor list overflow,boost neigh_modify one.
   请问出现这样的错误应该怎样修改呢？与上面转data时的警告有关吗？

请各位老师指教！

Lacrimosa

可能是在生成data文件的时候缺少了相关的参数，可以查看一下相应的原子成键参数在in文件里是否存在。例如报错里提到的atom 281 303

cdcdcd

Lacrimosa 发表于 2020-10-27 18:30
可能是在生成data文件的时候缺少了相关的参数，可以查看一下相应的原子成键参数在in文件里是否存在。例如报 ...

我的in文件里面只有读取data文件的命令，没有键角的详细信息
但是在data文件中可以找到atom 281 303 的键的信息

Lacrimosa

cdcdcd 发表于 2020-10-27 19:53
我的in文件里面只有读取data文件的命令，没有键角的详细信息
但是在data文件中可以找到atom 281 303 的 ...

bond_style也设置了么？

cdcdcd

Lacrimosa 发表于 2020-10-27 20:24
bond_style也设置了么？

设置了
pair_style        lj/cut 12.5
atom_style  full
bond_style  harmonic
angle_style  harmonic  
dihedral_style  harmonic
improper_style  cvff
boundary        p p p
read_data       CO2-130.data
pair_coeff        * * 100.0 3.0
timestep        0.005

Lacrimosa

cdcdcd 发表于 2020-10-27 20:34
设置了
pair_style        lj/cut 12.5
atom_style  full

确认一下data文件中bond_coeff是否齐全，timestep单位如果是fs可能有些太大了，设为1-2fs比较合适。再查看一下报错的部分原子是不是距离太远了

cdcdcd

Lacrimosa 发表于 2020-10-27 20:38
确认一下data文件中bond_coeff是否齐全，timestep单位如果是fs可能有些太大了，设为1-2fs比较合适。再查 ...

老师，您好，data文件中包含了bond_coeff信息，我修改了timestep的数值，改为0.001了，前面弛豫的部分可以正常算出来了，但后面计算的部分报错了
报错为
dump            21 2 custom 100000 dump.atom id type x y z vx vy
variable        msd11 atom c_18[1]
variable        msd12 atom c_18[2]
variable        msd13 atom c_18[3]
variable        msd14 atom c_18[4]
dump            22 2 custom 100000 dump.msd.atom id v_msd11 v_msd12 v_msd13 v_msd14
dump            23 2 xyz 100000 dump.XYZ.xyz
run                          10000000000
ERROR: Invalid run command N value (../run.cpp:112)
Last command: run                          10000000000
是运行的时间步长设置的不合理，我去lammps官网查的错误信息是这样的
Invalid run command N value
The number of timesteps must fit in a 32-bit integer. If you want to run for more steps than this, perform multiple shorter runs.
如果想将run命令的步长拆分计算，应该要怎样设置命令呢？谢谢老师指教

Lacrimosa

cdcdcd 发表于 2020-10-28 11:32
老师，您好，data文件中包含了bond_coeff信息，我修改了timestep的数值，改为0.001了，前面弛豫的部分可 ...

你这个时间都10微秒了，太长了吧？先看看是否有必要跑这么久吧。run的N value应该是不能超过2147483647，所以可以分别设为多个run 每个的值小于上限就好啦。你可以先跑200ps让体系驰豫一下，然后把timestep改成0.002试试能不能成功，如果可以，steps就可以减少一倍了。然后设5个run，每个1000000000即可。
例如:
timestep 0.002
run 1000000000

timestep 0.002
run 1000000000
...

cdcdcd

Lacrimosa 发表于 2020-10-28 14:13
你这个时间都10微秒了，太长了吧？先看看是否有必要跑这么久吧。run的N value应该是不能超过2147483647， ...

老师，我的in文件是这样的
units                lj
pair_style        lj/cut 15
atom_style  full
bond_style  harmonic
angle_style  harmonic  
dihedral_style  harmonic
improper_style  cvff
boundary        p p p
read_data       h2-1500.data
pair_coeff        * * 100.0 3.0
timestep        0.001
group       1 type 1
group       2 type 2
neighbor        1 bin
neigh_modify        every 1 delay 0 check yes
minimize        1.0e-12 1.0e-12 50000 50000
velocity        all create 298 872877 dist gaussian
fix                11 all npt temp 298 298 0.1 iso 15 15 1 drag 0.2
thermo          1000
thermo_style    custom step temp pe ke etotal vol press
run             1000000
unfix           11

fix             12 all nvt temp 298 298 0.5 drag 0.2
compute         13 all pe/atom
compute         14 all ke/atom
compute         15 all coord/atom cutoff 15
compute         16 all temp
compute         17 all pressure 16
compute         18 2 msd
compute         19 2 rdf 100 cutoff 15
fix             20 2 ave/time 100000 1 100000 c_19[1] c_19[2] c_19[3] file dump.rdf.h2.rdf mode vector

dump            21 2 custom 100000 dump.atom id type x y z vx vy
variable        msd11 atom c_18[1]
variable        msd12 atom c_18[2]
variable        msd13 atom c_18[3]
variable        msd14 atom c_18[4]
dump            22 2 custom 100000 dump.msd.atom id v_msd11 v_msd12 v_msd13 v_msd14
dump            23 2 xyz 100000 dump.XYZ.xyz
run                          10000000
一共需要模拟10ns，您说的单位我不是很理解，我这里的步长0.001，它指的是0.001fs,还是1fs呀？

Lacrimosa

cdcdcd 发表于 2020-10-30 09:15
老师，我的in文件是这样的
units                lj
pair_style        lj/cut 15

你用的单位是lj，应该不适合你的力场，你最好查看一下力场参数的单位，然后再看看lammps手册里units那一部分，确定好单位制后就知道timestep的单位了，
举个例子，如果
units real
timestep 1 ；这里timestep的单位是fs，所以1就是1fs
run  10000000 ；这里可知总时长=10000000*1fs=10ns

如果untis是metal
units metal
timestep 0.01 ；这里单位是ps，0.01ps即为1fs
run 10000000 ;总时长=10000000*0.01ps=10ns

cdcdcd

老师，该怎么去查看力场参数的单位呢？
我这个使用的力场是lj/cut。手册里面的介绍是这样的
ϵ  (energy units)
σ  (distance units)
LJ cutoff (distance units)
然后查看了units部分，不知道选哪一种是合适的，lj单位是这样的
mass = mass or m
distance =σ , where x* = x / σ
time = tau, where t* = t (ϵ  / m / σ ^2)^1/2
energy = ϵ, where E* = E /ϵ
velocity = σ /tau, where v* = v tau /σ
temperature = reduced LJ temperature, where T* = T Kb / ϵ
pressure = reduced LJ pressure, where P* = P σ^3 / ϵ
想请问这里的单位和力场的单位是合适的吗？
*号在这里表示什么意思？
下降的LJ温度，下降的LJ压力如何理解？
谢谢老师指点。

悟小空

请问你解决了吗，我是用pcff算环氧树脂，跟你出现类似的问题

爱模拟的bao

我也遇到了这个问题，想问下你这解决了吗？

Lacrimosa

cdcdcd 发表于 2020-11-10 15:18
老师，该怎么去查看力场参数的单位呢？
我这个使用的力场是lj/cut。手册里面的介绍是这样的
ϵ  (en ...

lj是约化单位, 你如果是用msi2lmp生成的力场参数应该看看msi2lmp中有没有提供相关的信息

ohlalaoh

大哥  我也遇到第一个问题  你咋解决的 有影响吗