计算化学公社

标题: 我用mtp势函数来跑lammps压缩Ag2S，过程不报错但是也不输出 [打印本页]

作者
Author: alv 时间: 2025-10-10 15:59
标题: 我用mtp势函数来跑lammps压缩Ag2S，过程不报错但是也不输出
两个附件一个是输入文件，一个是mtp势函数文件，为方便上传都改为txt文件。

mpirun -np 20 lmp_mpi -in Ag2S-fuxian.imp
LAMMPS (22 Jul 2025 - Update 1)
Reading data file ...
  triclinic box = (0 0 0) to (106.35932 89.659726 94.530658) with tilt (0 -32.89871 0)
  5 by 2 by 2 MPI processor grid
  reading atoms ...
  50700 atoms
  reading velocities ...
  50700 velocities
  read_data CPU = 0.303 seconds
WARNING: No fixes with time integration, atoms won't move
For more information see https://docs.lammps.org/err0028 (../verlet.cpp:60)
Wrapper initialization
abinitio = false
calculate_efs = true
check_errors = false
extrapolation_control = true
fit = false
lotf = false
mlip = true
select = false
write_cfgs = false
write_cfgs:skip = 0
mlip:load_from = pot.mtp08
MLIP type is non-linear MTP
Selection: loading state from "pot.mtp08"
No selection data provided with MLIP
这个是输入完出现的，然后就什么都没了。是因为我的输入文件有问题吗？求指导

(, 下载次数 Times of downloads: 6) (, 下载次数 Times of downloads: 5)

作者
Author: 海棠依旧cdut 时间: 2025-10-16 10:49
虽然我和你做的领域不一样，代码也看不出来什么问题。但是你说既不报错也不出结果的问题我也遇到过，我只能从个人经验上分享一下我的原因：minimize消耗了一定的步数，导致后续分析没有跑够输出一次的步骤就结束了。我的处理方式是在能量最小化后添加一个 reset_timestep 0。你可以试试这样设置后有没有改变
# 能量最小化
min_style cg
minimize 1e-10 1e-10 5000 10000
reset_timestep 0

欢迎光临计算化学公社 (http://bbs.keinsci.com/)