计算化学公社

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

[xtb] 我模拟几个有机分子形成晶体的结构,xtb怎么算周期性?

[复制链接 Copy URL]

72

帖子

0

威望

387

eV
积分
459

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
我有几个有机分子形成晶体,预先预测能不能形成,用xtb怎么算周期性的结构?

331

帖子

3

威望

1764

eV
积分
2155

Level 5 (御坂)

2#
发表于 Post on 2025-3-18 15:28:16 | 只看该作者 Only view this author
可以试试直接用cp2k算。cp2k不仅支持GFN0-xTB和GFN1-xTB方法,还能用PBE-D3做精度更好的结构优化。
可以在windows用wsl很方便的安装cp2k:cp2k-2025.1在Debian-wsl2中的安装(附所需所有第三方库)
cp2k的输入文件可以借助multiwfn生成:使用Multiwfn非常便利地创建CP2K程序的输入文件
cp2k的系统学习见:第4届北京科音CP2K第一性原理计算培训班第一轮通知

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
FRRMM + 1 赞!

查看全部评分 View all ratings

515

帖子

1

威望

4699

eV
积分
5234

Level 6 (一方通行)

3#
发表于 Post on 2025-3-19 09:07:21 | 只看该作者 Only view this author
这是晶体结构搜索问题,在论坛里搜索,这个问题不是简单算算单点能或者结合能的问题

32

帖子

1

威望

1604

eV
积分
1656

Level 5 (御坂)

4#
发表于 Post on 2025-3-19 09:50:52 | 只看该作者 Only view this author
可以用dftbplus, 内有GFN1/2-xTB方法,可算周期性结构

6万

帖子

99

威望

5万

eV
积分
120169

管理员

公社社长

5#
发表于 Post on 2025-3-19 13:22:48 | 只看该作者 Only view this author
xtb自己计算周期性的功能只支持GFN-FF力场,北京科音高级量子化学培训班(http://www.keinsci.com/KAQC)里详细xtb的时候我还专门给了个例子。其它的参考前面的回答。
GFN-FF对于有机分子间相互作用能的描述不明显逊于GFN1或2-xTB。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

72

帖子

0

威望

387

eV
积分
459

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2025-3-19 15:49:24 | 只看该作者 Only view this author
tjuptz 发表于 2025-3-19 09:07
这是晶体结构搜索问题,在论坛里搜索,这个问题不是简单算算单点能或者结合能的问题

是的,就是做构型搜索。我先用molclus产生500个初始构型,然后用crest计算结合后的结构。我用已知文献的一篇文献来算,如果不用周期性算,前面50个结构都没有文献的晶体结构。不知道用什么方法来预测了?

72

帖子

0

威望

387

eV
积分
459

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2025-3-19 16:11:53 | 只看该作者 Only view this author
本帖最后由 wilely 于 2025-3-19 16:29 编辑
sobereva 发表于 2025-3-19 13:22
xtb自己计算周期性的功能只支持GFN-FF力场,北京科音高级量子化学培训班(http://www.keinsci.com/KAQC)里 ...

我现在用molclus产生500个初始构型,然后这500结构要怎么处理成周期性结构,再用GFN-FF一个个的优化吗,可以用什么软件自动全优化吗?

72

帖子

0

威望

387

eV
积分
459

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2025-3-19 17:32:51 | 只看该作者 Only view this author
sobereva 发表于 2025-3-19 13:22
xtb自己计算周期性的功能只支持GFN-FF力场,北京科音高级量子化学培训班(http://www.keinsci.com/KAQC)里 ...

刚查了一下,现在ASE可以调用xtb算周期性,还有tblite可以直接用。

515

帖子

1

威望

4699

eV
积分
5234

Level 6 (一方通行)

9#
发表于 Post on 2025-3-20 17:20:09 | 只看该作者 Only view this author
wilely 发表于 2025-3-19 15:49
是的,就是做构型搜索。我先用molclus产生500个初始构型,然后用crest计算结合后的结构。我用已知文献的 ...

团簇没有晶胞相关参数,算不上晶体,可以了解下uspex和calypso

72

帖子

0

威望

387

eV
积分
459

Level 3 能力者

10#
 楼主 Author| 发表于 Post on 2025-3-20 17:48:04 | 只看该作者 Only view this author
tjuptz 发表于 2025-3-20 17:20
团簇没有晶胞相关参数,算不上晶体,可以了解下uspex和calypso

不是团族,是几个有机分子形成的共晶。我用vesta取出晶胞做优化,vasp计算总是结束后输出NaN,没结果,不知道哪里错误

515

帖子

1

威望

4699

eV
积分
5234

Level 6 (一方通行)

11#
发表于 Post on 2025-3-20 18:17:04 | 只看该作者 Only view this author
wilely 发表于 2025-3-20 17:48
不是团族,是几个有机分子形成的共晶。我用vesta取出晶胞做优化,vasp计算总是结束后输出NaN,没结果,不 ...

你都有晶胞了还预测个啥劲儿

72

帖子

0

威望

387

eV
积分
459

Level 3 能力者

12#
 楼主 Author| 发表于 Post on 2025-3-20 20:10:14 | 只看该作者 Only view this author
tjuptz 发表于 2025-3-20 18:17
你都有晶胞了还预测个啥劲儿

我重复文献的一个晶体,然后变换一个分子,看能不能形成晶体?

515

帖子

1

威望

4699

eV
积分
5234

Level 6 (一方通行)

13#
发表于 Post on 2025-3-21 09:56:10 | 只看该作者 Only view this author
wilely 发表于 2025-3-20 20:10
我重复文献的一个晶体,然后变换一个分子,看能不能形成晶体?

那你用文献里的晶胞替换分子也不太合理,这种是用uspex或者calypso做晶体结构搜索的研究范畴

72

帖子

0

威望

387

eV
积分
459

Level 3 能力者

14#
 楼主 Author| 发表于 Post on 2025-3-21 10:19:44 | 只看该作者 Only view this author
tjuptz 发表于 2025-3-21 09:56
那你用文献里的晶胞替换分子也不太合理,这种是用uspex或者calypso做晶体结构搜索的研究范畴

好的。这两个软件用过吗?复杂吗?

515

帖子

1

威望

4699

eV
积分
5234

Level 6 (一方通行)

15#
发表于 Post on 2025-3-21 10:40:15 | 只看该作者 Only view this author
wilely 发表于 2025-3-21 10:19
好的。这两个软件用过吗?复杂吗?

我没用过,这种工作量比较大。你可以看手册了解下选用哪个

本版积分规则 Credits rule

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

GMT+8, 2025-8-16 20:38 , Processed in 0.180585 second(s), 21 queries , Gzip On.

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