ShangChien 发表于 2021-10-10 19:46
sob老师，我好像才明白您的意思。也就是说gaussian会自动识别“D * * * * F”中的通配符*喽。学到了！！

是的
优化任务一开始会输出一个冗余内坐标列表，从里面也能直接看到当前哪些被冻结了

sobereva 发表于 2021-10-9 16:09
根本就不需要任何脚本，直接opt=modredundant，末尾空一行写D * * * * F就能冻结所有二面角

sob老师，我好像才明白您的意思。也就是说gaussian会自动识别“D * * * * F”中的通配符*喽。学到了！！

sobereva 发表于 2021-10-9 16:09
根本就不需要任何脚本，直接opt=modredundant，末尾空一行写D * * * * F就能冻结所有二面角

sob老师授课辛苦了！我这里的需求就是，批量找出确定每个二面角"D * * * * F"表达式的四个原子序号。我的课题需要高通量计算大量分子，高斯手动固定太慢了。

根本就不需要任何脚本，直接opt=modredundant，末尾空一行写D * * * * F就能冻结所有二面角