计算化学公社

标题: 请问如何计算环糊精金属有机框架（CD-MOF-1）体系的原子电荷 [打印本页]

作者
Author: clh 时间: 2023-4-4 10:54
标题: 请问如何计算环糊精金属有机框架（CD-MOF-1）体系的原子电荷
尝试模拟环糊精金属有机框架（CD-MOF-1）与小分子的相互作用，根据文献Nanoporous Carbohydrate MetalOrganic Frameworks (doi: 10.1021/ja208224f)中所给出的CCDC号(773709)下载得到cif文件。目前仅会一点gromacs，ORCA和CP2K都没整明白。
一开始是抠出单个环糊精（不知道将与钾离子配位的氧补成水分子的做法对不对）用ORCA的RESP2懒人脚本算的，发现小破笔记本算不动。然后是再把环糊精里的葡萄糖抠出来打算计算完再补全环糊精的原子电荷（不知道是不是瞎整）。算了半天也不知道算不算的出来，下含当前的inp.out文件。
然后看公社以往的帖子，发现有说用Multiwfn里内含的EEM算的，试了一下发现没有钾的参数

。
就是想整出原子电荷后续分析一下MOF和小分子的相互作用，也不需要非常精确的结果。限于本人的水平，求教有没有简单的方法得到说得过去的原子电荷数据。

作者
Author: sobereva 时间: 2023-4-4 12:52
最耗时的过程是几何优化，还可能碰上难收敛
鉴于你是基于晶体结构，质量还说得过去，对于要求不高的情况，直接拿第二个图里的结构在B3LYP/6-31G*下算个单点得到波函数文件，然后给Multiwfn算RESP电荷就完。小破笔记本也算得动

作者
Author: clh 时间: 2023-4-4 21:33
本帖最后由 clh 于 2023-4-5 08:52 编辑

感谢老师，原子电荷算出来了。

作者
Author: bobosiji 时间: 2023-4-5 09:59

sobereva 发表于 2023-4-4 12:52
最耗时的过程是几何优化，还可能碰上难收敛
鉴于你是基于晶体结构，质量还说得过去，对于要求不高的情况， ...

楼主把葡萄糖和K+一起算，应该是要考虑它们之间的电荷转移了。请教3个问题：
1、看楼主的截图，算电荷时候应该没考虑环糊精金属有机框架（CD-MOF-1）的mof和环糊精之间的电荷转移；这个近似是否合理？
2、楼主最后一个截图，相当于从周期性的环糊精钾离子复合物里截取一个葡萄糖钾离子单元来算电荷，这个葡萄糖用H还是-CH3来封端呢？
3、进一步近似，直接算葡萄糖（不包含K+；即不考虑葡萄糖与K+之间的电荷转移）的电荷，然后用该电荷模拟环糊精钾离子相互作用，能否得到合理结果？
多谢提点~

欢迎光临计算化学公社 (http://bbs.keinsci.com/)