求问做六水合三价铁和氯离子的优化步骤

dbxhtw

请问在优化六水合三价铁和氯离子的时候，我是先用片段波函数#p PBE1PBE/gen  em=gd3bj guess(fragment=3) nosymm scf=(maxcyc=512) ，然后读取该文件做波函数稳定性优化#p PBE1PBE/gen  em=gd3bj guess=read  nosymm scf=(maxcyc=512)  geom=allcheck stable=opt  scrf(SMD,solvent=Water) ，然后在用得到的结果去做opt，请问这系列流程是正确的吗？还有就是第二步需要加隐式溶剂模型吗

sobereva

根本不需要片段组合波函数那么麻烦，这体系又不是什么难算的东西。没事甭把SCF迭代次数上限设那么大，看
解决SCF不收敛问题的方法
http://sobereva.com/61

stable=opt都没什么必要，恰当设定自旋多重度对这种简单体系直接就能收敛到稳定波函数

优化时用溶剂模型的必要性下文说了
谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327（http://bbs.keinsci.com/thread-3345-1-1.html）
不推荐opt过程用SMD

dbxhtw

谢谢老师的回答，我是需要分别求出三价铁和氯离子  以及  二价铁和氯自由基的  这两个体系的能量，在计算电荷和自旋多重度的时候，我发现这连个体系很难区分，都是2,2或者2,4，或者2,6.所以想试试片段波函数，

dbxhtw

另外还有个问题想请教，计算三价铁和氯离子的时候，应该是五水合三价铁再配位一个氯，还是六水合三价铁旁边放一个氯呢

sobereva

dbxhtw 发表于 2025-11-16 16:02
谢谢老师的回答，我是需要分别求出三价铁和氯离子  以及  二价铁和氯自由基的  这两个体系的能量，在计算电 ...

片段组合波函数设的仅仅是初猜，除非一个体系的SCF参数空间真的有不同极小点同时你又设得恰到好处，否则SCF过程该往哪收敛还往哪收敛。对于当前情况你的这个想法完全行不通
Fe的氧化态是多少，要算完了用Multiwfn按下文确认
使用Multiwfn通过LOBA方法计算氧化态
http://sobereva.com/362（http://bbs.keinsci.com/thread-5118-1-1.html）
而Cl是什么状态，应当考察自旋密度
谈谈自旋密度、自旋布居以及在Multiwfn中的绘制和计算
http://sobereva.com/353（http://bbs.keinsci.com/thread-4422-1-1.html）

sobereva

dbxhtw 发表于 2025-11-16 16:18
另外还有个问题想请教，计算三价铁和氯离子的时候，应该是五水合三价铁再配位一个氯，还是六水合三价铁旁边 ...

要算哪种状态是自己决定的不是别人告诉你的。如果你想算能量最低的配位情况，试图计算Cl-和水的配体交换能便知
顺带一提，中性水溶液中六水合铁(III)离子中的一个配位水记得用OH-形态，这是此物质最主要的存在状态

dbxhtw

谢谢sob老师！