Deepast 发表于 2023-3-30 13:25 嗯嗯 十分感谢  |
隔壁老张sir 发表于 2023-3-30 11:53 这个是要算一系列单点能，然后构建全维势能面 结构可以通过编程生成 算了单点后可以根据一些算法拟合势能面，这个一般也要自己写 你想练手的话建议开始把计算精度调低点 我懂的也不多，只能说到这 |
Deepast 发表于 2023-3-30 11:44 唉 我是组里第一个做这个的人，导师也不懂这个方向（方向是大老板硬给我的） |
隔壁老张sir 发表于 2023-3-30 10:05 准经典轨线的研究是吧，这个研究不是靠这种扫描做的 做这个的比较少，建议多问问导师和师兄师姐 |
sobereva 发表于 2023-3-30 09:27 sob老师好。我是参考《Global Ab Initio Ground-State Potential Energy Surface of N4》这篇文献，尝试构建N4体系的势能面。我这里是选了文章中的一种N2+N2的构型，按文章中给的键长变化参数，尝试计算一下势能变化。 “casscf(12,12)/maug-cc-pvtz”是该文献中用的，用它的实际意义我也不太清楚；“int=ultrafine Guess=Permute SCF=maxcyc=500”是我为了收敛加的；过程中有共价键断裂，就又加了“guess(mix,always)”；为了防止扫描过程分子跳变，加的“nosymm”。 sob老师说的“尝试不同定义方式”是指什么呀？ 我一开始也想着不用GIC，但是没有找到其他办法使扫描过程中两氮分子的中心距离保持不变，sob老师有什么好的建议吗？ |
|
不知道你想算什么，这关键词完全莫名其妙，看不出这计算有任何实际意义 GIC的算法并不稳健，尝试不同定义方式，或者根据实际情况考虑不用GIC的可能 |
Deepast 发表于 2023-3-30 02:39 您说的“先续算看看能不能跑，不行就手动调下”可以麻烦您说的具体一点吗？ 我刚接触量化不久，是个小白，不太知道具体是要怎么做  |
|
最好把广义坐标也简化下，我不确定程序崩和这个有没有关系 比如冻结的话没必要写这么多 开方和单位转换那部分可以删 |
|
gic有时候就会这样，程序不知道该怎么优化 我不知道有什么好的解决方法，有时候甚至交换原子顺序就能跑 先续算看看能不能跑，不行就手动调下，再不行换别的扫描方法 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-21 21:09 , Processed in 0.169760 second(s), 25 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites