|
本帖最后由 李海恩 于 2024-8-1 15:28 编辑
前情提要:本人太菜,算个过渡金属氢氧化物切完表面怎么都不收敛;大佬说要不你试试CHGNet优化完再丢进vasp可以算快一点
分享一下过程。。万一有人用呢.jpg
因为不涉及MD所以直接在WSL里使用了,请各位大佬批评指正
1.CHGnet简介:在LAMMPS中使用CHGNethttp://bbs.keinsci.com/forum.php?mod=viewthread&tid=40402&fromuid=17159
官方:https://chgnet.lbl.gov
2.Windows下使用WSL:https://zhuanlan.zhihu.com/p/698816488,我装的是Ubuntu-22.04
3.安装chgnet:
(1)安装pip:https://www.php.cn/faq/660949.html
(2)按照官网教程输入pip install chgnet(第一次装会下载一堆需要的软件
(3)装完以后可以去官网Usage部分按需复制需要的脚本使用
ps.苯人只用Structure Optimization部分做预优化,复制官网脚本运行以后发现会报numpy不兼容,所以按照提示重新装了屏幕所示版本
希望大家都不报错 阿门
4.使用
解决兼容性问题以后发现官网脚本直接运行报错(别骂了别骂了,不会敲代码寸步难行QAQ)
于是去找了万能的gpt4.0,改来改去终于改出一版可以自动输出final structure为笛卡尔坐标并保存为poscar文件的脚本
(如果是分数坐标就不用改这么多,我这个版本的numpy好像默认输出分数坐标)
随便从电脑里找一个Co(OH)2运行脚本测试(本结构先用VASP优化过,所以两步就跑完了)
发现屏幕直接输出有坐标和磁矩
大师 我(好像)练成了!
附件是gpt写的脚本,调用当前目录的POSCAR优化,并自动输出opt.poscar
ps.chgnet算得蛮快,扩胞以后90个原子一秒竟然可以跑5步!
优化完再丢进vasp比直接用vasp算快一半时间(仅针对本体系)
以上就是计算机小白使用chgnet过程,特别鸣谢chatgpt敲代码,祝大家计算都不报错!
|
-
-
gpt3.py
1.5 KB, 下载次数 Times of downloads: 1
gpt修改输出笛卡尔坐标
-
-
opt.py
609 Bytes, 下载次数 Times of downloads: 1
默认可能输出分数坐标
评分 Rate
-
查看全部评分 View all ratings
|