计算化学公社

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

[CP2K] CP2K优化亚铁螯合物框架结构的参数设置问题

[复制链接 Copy URL]

9

帖子

0

威望

527

eV
积分
536

Level 4 (黑子)

本帖最后由 默默 于 2025-5-8 22:21 编辑

各位老师好,
我最近用了类似图片中的亚铁螯合物单体合成了两例螯合物的共价有机框架材料,目前已经解析出了晶体结构,现在想用CP2K做点计算,主要是弱相互作用分析那些。
在结构优化的时候我遇到两个问题。

(1)高斯算单体的自旋多重度为5的时候能量最低,但是材料实验上测出来是顺磁的,是不是自旋多重度设为1就行?

(2)Fe(II)的体系是否需要DFT+U(因为是新材料我没找到类似体系参考文献的U值,MOF也没有相近的)?   
不+U进行 晶胞+结构优化(自旋多重度设置为1)前后F和N配位键键长基本一致(实验值1.93,优化后1.92)
如果用杂化泛函的话,结构优化是否也需要用杂化泛函?(一个体系约500个原子/另一个体系约750个原子)







6万

帖子

99

威望

5万

eV
积分
120101

管理员

公社社长

2#
发表于 Post on 2025-5-9 02:33:11 | 只看该作者 Only view this author
1 既然都有了晶体结构,显然直接用CP2K算就完了
材料是顺磁的绝对不代表各个位点都是闭壳层;若每个Fe位点都有自旋磁矩,但方向随意,整体也是顺磁。当成闭壳层计算明显不合理

2 我没看见体系里哪有F
CP2K用杂化泛函算这东西晶胞的单点轻轻松松,优化太耗时。纯泛函优化+杂化泛函单点是常用组合
北京科音自然科学研究中心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

9

帖子

0

威望

527

eV
积分
536

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2025-5-9 18:03:22 | 只看该作者 Only view this author
sobereva 发表于 2025-5-9 02:33
1 既然都有了晶体结构,显然直接用CP2K算就完了
材料是顺磁的绝对不代表各个位点都是闭壳层;若每个Fe位点 ...

非常感谢老师的回复,我还想请教您一下我应该怎么设置参数呢,还是用  MULTIPLICITY    1  然后打开UKS吗?
1个体系晶胞内有4个Fe(II),一个体系有8个Fe(II),各个原子初猜的磁矩应该怎么设置呢?

6万

帖子

99

威望

5万

eV
积分
120101

管理员

公社社长

4#
发表于 Post on 2025-5-9 18:54:56 | 只看该作者 Only view this author
默默 发表于 2025-5-9 18:03
非常感谢老师的回复,我还想请教您一下我应该怎么设置参数呢,还是用  MULTIPLICITY    1  然后打开UKS吗 ...

光写MULTIPLICITY    1和UKS毫无意义,结果和RKS完全一样,初猜对称性没有打破

CP2K没法对各个位点有磁矩但是整体是顺磁性(自旋磁矩无序)的状态进行计算。实际计算上,只能是把各种能算的自旋磁矩排列的情况都算算(顺磁性、铁磁性...通过恰当设置MAGNETIZATION实现),取能量最低的态
北京科音自然科学研究中心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

9

帖子

0

威望

527

eV
积分
536

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2025-5-14 21:37:46 | 只看该作者 Only view this author
sobereva 发表于 2025-5-9 18:54
光写MULTIPLICITY    1和UKS毫无意义,结果和RKS完全一样,初猜对称性没有打破

CP2K没法对各个位点有 ...

sob老师打扰了,我这几天思考了下,按我的理解不同自旋态的能量也是要在杂化泛函下算单点,但是我的体系有一个轴只有9埃,没法用ADMM做杂化泛函计算,耗时可能比较大。

由于我做计算的初衷只想可视化弱相互作用力,对于这种体系是不是直接用准分子近似的NCI或者IGM方法显示下材料堆积的弱相互作用力就可以了(分子力学优化H,弱相互作用力主要是pi-pi堆积和C-H…pi)

6万

帖子

99

威望

5万

eV
积分
120101

管理员

公社社长

6#
发表于 Post on 2025-5-14 23:25:03 | 只看该作者 Only view this author
默默 发表于 2025-5-14 21:37
sob老师打扰了,我这几天思考了下,按我的理解不同自旋态的能量也是要在杂化泛函下算单点,但是我的体系 ...

“只有9埃,没法用ADMM做杂化泛函计算” 前后根本没有直接的因果关系

如果你不想基于波函数做弱相互作用可视化分析,用今年我新提出的mIGM,比promolecular的NCI和IGM都好得多,介绍看
Tian Lu, Graphically revealing weak interactions in dynamic environments using amIGM method, Structure and Bonding. Springer, Berlin, Heidelberg. DOI: 10.1007/430_2025_95
Tian Lu, Visualization Analysis of Covalent and Noncovalent Interactions in Real Space, Angew. Chem. Int. Ed. (2025) DOI: 10.1002/anie.202504895
例子在最新版Multiwfn手册里搜mIGM。用法和下文的IGMH没有任何区别,只不过选IGMH的地方改成选mIGM
使用Multiwfn做IGMH分析非常清晰直观地展现化学体系中的相互作用
http://sobereva.com/621http://bbs.keinsci.com/thread-28147-1-1.html

北京科音自然科学研究中心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

9

帖子

0

威望

527

eV
积分
536

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2025-5-14 23:56:08 | 只看该作者 Only view this author
sobereva 发表于 2025-5-14 23:25
“只有9埃,没法用ADMM做杂化泛函计算” 前后根本没有直接的因果关系

如果你不想基于波函数做弱相互作 ...

收到,谢谢sob老师,我学习学习

236

帖子

0

威望

1227

eV
积分
1463

Level 4 (黑子)

实验组内的DFT计算、第一性原理、MD模拟爱好者

8#
发表于 Post on 2025-7-2 14:42:32 | 只看该作者 Only view this author
本帖最后由 JCenter 于 2025-7-2 14:43 编辑
sobereva 发表于 2025-5-9 02:33
1 既然都有了晶体结构,显然直接用CP2K算就完了
材料是顺磁的绝对不代表各个位点都是闭壳层;若每个Fe位点 ...

卢老师,我有个问题想请教下您:在gaussian中,非周期性有机分子往往同时是结构优化和振动分析同时做的。振动分析的一个目的是看有无有无虚频,是否为极小值点。在cp2k中周期性物质的结构优化和振动分析没法同时做,但是计算级别也往往要保持一致。在cp2k中,是不是每次做完opt或者cell_opt后都要进行振动分析,查看有无虚频呢。
目前专攻:
基于DFT、MD模拟的自由能计算
基于过渡态理论的自由能垒和反应速率常数计算
基于MD模拟的交联聚合物计算

6万

帖子

99

威望

5万

eV
积分
120101

管理员

公社社长

9#
发表于 Post on 2025-7-3 03:56:45 | 只看该作者 Only view this author
JCenter 发表于 2025-7-2 14:42
卢老师,我有个问题想请教下您:在gaussian中,非周期性有机分子往往同时是结构优化和振动分析同时做的。 ...

拿不准是不是极小点的情况应该这样。
北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2025-8-13 22:02 , Processed in 0.189703 second(s), 23 queries , Gzip On.

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