求助，如何观看vasp结构优化的力的变化趋势

XIEH

高斯的结构优化可以中途下载log文件用GaussView的optimization功能，看能量和力的变化趋势。
vasp的OSZICAR可以看到每一个离子步的E0，那么力是看哪个呢？

乐平

Alston453 发表于 2024-11-7 15:50
File "check_VASP_convergency.py", line 203
    raise FileNotFoundError(f"The file {file_path.nam ...

你的 POSCAR, OUTCAR, OSZICAR 是不是改名了？

如果改名了，需要按我帖子最后写的，输入参数 和 新的文件名

http://bbs.keinsci.com/forum.php ... 65&fromuid=1532

1. python check_VASP_convergency.py -p 你的新POSCAR -o 你的新OUTCAR -s 你的新OSZICAR

复制代码

Alston453

乐平 发表于 2024-6-28 17:58
之前的有错误，我更新了一下脚本。

File "check_VASP_convergency.py", line 203
    raise FileNotFoundError(f"The file {file_path.name} does not exist.")
请问输入您说的运行命令后出现这个问题怎么解决呀大佬

乐平

阿宁 发表于 2024-10-26 17:15
您好，感谢分享。我在使用你的脚本时出现如下报错，百度找不到解决方法，想请教一下。
[qxn@Master ~/00 ...

报错里写得很清楚啊

你的 numpy 版本太高了（2.0版），里面没有 np.asfarray 这个函数了，需要把它换成 np.asarray

阿宁

乐平 发表于 2024-5-12 22:32
可以我写的 python 脚本

基本用法，输入

您好，感谢分享。我在使用你的脚本时出现如下报错，百度找不到解决方法，想请教一下。
[qxn@Master ~/001/yzcu4]$ python check_VASP_convergency.py
`np.asfarray` was removed in the NumPy 2.0 release. Use `np.asarray` with a proper dtype instead.
Traceback (most recent call last):
  File "/home/qxn/001/yzcu4/check_VASP_convergency.py", line 216, in <module>
    main()
  File "/home/qxn/001/yzcu4/check_VASP_convergency.py", line 209, in main
    F = grab_Force(poscar, outcar)
  File "/home/qxn/001/yzcu4/check_VASP_convergency.py", line 106, in grab_Force
    data = np.asfarray(data).reshape((count,atomNum,3))
  File "/home/qxn/.local/lib/python3.9/site-packages/numpy/__init__.py", line 397, in __getattr__
    raise AttributeError(
AttributeError: `np.asfarray` was removed in the NumPy 2.0 release. Use `np.asarray` with a proper dtype instead.

XIEH

Yun_ZL 发表于 2024-9-10 09:29
我发现编译VTSTcode之后，可以使用快捷命令grep FORCES：OUTCAR进行查看力的变化，没有编译vtstcode的情况 ...

谢谢你的分享，但是我的估计没有编译VTSTcode，输入命令没有反应

Yun_ZL

我发现编译VTSTcode之后，可以使用快捷命令grep FORCES：OUTCAR进行查看力的变化，没有编译vtstcode的情况下vasp的outcar并没有输出力的信息

lxt999

你好 我想请问一下 我在算AIMD分子动力学 查看势能 也可以用这个脚本吗

XIEH

KangBao 发表于 2024-7-27 19:24
我的方法是直接在终端输入"grep FORCES: OUTCAR"获取力的变化，如果有做图需要可以获取力的列进行做图。脚 ...

感谢回复，已解决

KangBao

我的方法是直接在终端输入"grep FORCES: OUTCAR"获取力的变化，如果有做图需要可以获取力的列进行做图。脚本更直观，但这样也可以很方便的浏览趋势

乐平

Minshen 发表于 2024-6-24 11:03
您好，请教一下，我用您的脚本，运行出来，力的图上只有一个点，改变横轴纵轴范围也是这样,是咋回事呢

之前的有错误，我更新了一下脚本。

Minshen

您好，请教一下，我用您的脚本，运行出来，力的图上只有一个点，改变横轴纵轴范围也是这样,是咋回事呢

XIEH

Weldingspock 发表于 2024-5-13 17:24
我也推荐一个，虽然不能画图，直接显示在屏幕上也能清楚看到
https://github.com/Ionizing/usefultools-fo ...

谢谢

Weldingspock

我也推荐一个，虽然不能画图，直接显示在屏幕上也能清楚看到
https://github.com/Ionizing/usefultools-for-vasp

XIEH

卡开发发 发表于 2024-5-13 05:50
可以使用python库ase。ase gui XDATCAR，会显示结构轨迹和能量变化曲线(i,e-E[-1]，此处i为步数，e为能量， ...

感谢回复