计算化学公社

 找回密码 Forget password
 注册 Register
Views: 573|回复 Reply: 6
打印 Print 上一主题 Last thread 下一主题 Next thread

[辅助/分析程序] 修改cube格点文件中的分子结构的脚本

[复制链接 Copy URL]

908

帖子

37

威望

5435

eV
积分
7083

Level 6 (一方通行)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 ggdh 于 2024-11-30 20:44 编辑

有些时候,出于各种奇葩原因,需要修改cube中的分子结构。
http://sobereva.com/125 喂给gpt后,让其帮忙写了个程序,用法如下:
1. 用gaussview打开cube文件(这里设为abc.cub)
2. 进行你需要的修改
3. 保存成gjf文件(这里设为xyz.gjf),和cub文件放到一起
4. 把脚本modcube.py和cub.gjf放到同一位置
5. 打开命令行并运行(需要安装python,并且把python加入环境变量中)
  1. python modcube.py abc.cub xyz.gjf
复制代码
产生的output.cube即为修改后的cube文件

modcube.py

5.62 KB, 下载次数 Times of downloads: 13

评分 Rate

参与人数
Participants 2
eV +7 收起 理由
Reason
xiaobo-haust + 2 精品内容
卡开发发 + 5 233333

查看全部评分 View all ratings

421

帖子

5

威望

1678

eV
积分
2199

Level 5 (御坂)

鸩羽

2#
发表于 Post on 2024-11-30 20:39:41 | 只看该作者 Only view this author
奇葩原因是指
某不知名实验组从苞米地里长出来的计算选手

449

帖子

0

威望

1561

eV
积分
2010

Level 5 (御坂)

3#
发表于 Post on 2024-11-30 21:18:31 | 只看该作者 Only view this author
gpt用于写程序太强大了
现在喂进去一本书都可以写一个可执行文件

2

帖子

0

威望

25

eV
积分
27

Level 2 能力者

4#
发表于 Post on 2025-1-1 19:02:48 | 只看该作者 Only view this author
奇葩原因之一,我计算的电荷密度差文件,其分子因为周期性边界条件导致边界的一些原子不显示,因此想修改一下cube文件中的原子的分布。

101

帖子

1

威望

645

eV
积分
766

Level 4 (黑子)

5#
发表于 Post on 2025-1-1 19:11:35 | 只看该作者 Only view this author
xiaobo-haust 发表于 2025-1-1 19:02
奇葩原因之一,我计算的电荷密度差文件,其分子因为周期性边界条件导致边界的一些原子不显示,因此想修改一 ...

那其实ase package很适合你,读进去直接替换atoms之后write就可以了😂,转换成chgcar也可以的

2

帖子

0

威望

25

eV
积分
27

Level 2 能力者

6#
发表于 Post on 2025-1-5 00:19:05 | 只看该作者 Only view this author
五十八 发表于 2025-1-1 19:11
那其实ase package很适合你,读进去直接替换atoms之后write就可以了😂,转换成chgcar也可以的

谢谢,我去尝试一下

101

帖子

1

威望

645

eV
积分
766

Level 4 (黑子)

7#
发表于 Post on 2025-2-17 23:23:31 | 只看该作者 Only view this author
xiaobo-haust 发表于 2025-1-5 00:19
谢谢,我去尝试一下

注意哦,ase的可能因为版本会有两处bug,注释掉就行了

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-18 06:36 , Processed in 0.268383 second(s), 32 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list