wzkchem5 发表于 2024-4-18 01:58 对老师，我确实是想尝试拟合经典分子动力学BMH势函数的参数，您有什么教程或者参考文献推荐么？ 我查过一些文献，但是感觉对计算这部分的描述很模糊。 如果说考虑实际体系，那是否需要考虑各成分间的配比？还是说包含这几种原子就可以。我是需要使用第一性原理在周期性体系研究还是仍然使用量子化学在孤立体系研究呢？ |
liu256136 发表于 2024-4-17 16:20 一样的问题，体系多参考性质太强，DFT不够可靠。如果体系比较大，多参考态方法算不动（或者虽然单点能算得动，但是结构优化/势能面扫描算不动），那可以用DFT优化结构，这个属于不得已而为之；但你这里只有两个原子，没道理不用多参考态方法。 如果是为了拟合势能函数而做的这种计算的话，个人建议是把重点放在计算实际体系上，例如MgO-Al2O3-MnS体系里面拉伸一根Mn-O键。相比之下，只有Mn、O两个原子的体系的势能曲线，不是那么重要。因为凡是涉及共价键/强配位键，尤其是可能涉及多重键的情况，多体效应都很显著，所以研究清楚任意两个原子之间的势能函数，并不足以阐释实际的多原子体系。 |
wzkchem5 发表于 2024-4-17 22:36 老师，我还想再请教您一下，我最初是用Gaussian做的柔性扫描，泛函分别使用了B3LYP和B2PLYP，使用DFT-D3矫正，基组是def2-TZVP，没有碰到不收敛的问题，那这个结果您认为可靠么？ 我想看一下MgO-Al2O3-MnS体系内各种原子间的势能随原子间距的变化情况，同时可能需要关注一下原子间的弱相互作用。所以其实最初就是想换一个高一点的理论方法，然后看了很多sob老师的帖子，所以切换到了ORCA，然后用的PWPB95+D3+def2-QZVPP这个组合。目前是除了Mn-O和Mn-S外，其他的都正常收敛了，如果我Mn-O和Mn-S更换了理论方法后，其他的是不是都需要重新计算？ |
wzkchem5 发表于 2024-4-17 22:36 谢谢老师的解答，我先去了解一下您说的这些，恶补一下专业知识。  |
liu256136 发表于 2024-4-17 15:20 至于你这个特定的体系，本来就是多参考特性极其强的体系，恐怕即使用PWPB95收敛了，结果也不可靠。建议改用多参考态方法，最少最少也得用MRCISD+Q，活性空间取整个价层空间加double shell。 |
本帖最后由 wzkchem5 于 2024-4-17 15:37 编辑 liu256136 发表于 2024-4-17 15:20 “线性相关”和“基组线性相关”是两码事。 SCF里面可能遇到的线性相关主要分为三种，一种是DIIS（或KDIIS等类似方法）矢量的线性相关，一种是轨道基线性相关，一种是辅助基线性相关。后两种统称为基组线性相关。 DIIS矢量线性相关的成因及解决方法见我那个链接，这里就不赘述了。 轨道基线性相关是因为存在过多的弥散函数，或原子间离得太近，具体成因、诊断方法及解决方法见我2楼的链接（注意这个链接是BDF手册的一部分，因此其中的关键词都是BDF的关键词，但原理和ORCA是相通的）。 辅助基线性相关和轨道基线性相关类似，只不过问题出现在辅助基组（例如/J、/JK或/C基组），具体成因、诊断方法及解决方法见orca手册的pivoted Cholesky decomposition章节。 |
|
为什么你认为这是基组线性相关问题？ 按https://bdf-manual.readthedocs.i ... tml#scfconvproblems判断你这里SCF不收敛的实际原因是什么，以及如何解决 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2024-11-24 23:22 , Processed in 0.162472 second(s), 25 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites