|
谢谢老师 |
zhuyoucai 发表于 2023-6-25 15:38 对 |
sobereva 发表于 2023-6-25 04:12 老师您好,我查了一些资料,我是否可以这样理解那句话。我首先通过下面的命令得到了稳定的波函数O3-T-ub3lyp-stableopt.chk,如下所示。 %chk=O3-T-ub3lyp-stableopt.chk %mem=24GB %nprocshared=24 # ub3lyp def2svp stable=opt 随后我将其log文件导出存成新的gjf,将输入改为下面,其中oldchk是上述chk文件的路径位置,下面是新生成的chk文件。那么我读取的应该是oldchk的波函数对吧,谢谢老师的解答。 %oldchk=/home/zhuyoucai/g09/BS-calculation/O3-T-ub3lyp-stableopt/O3-T-ub3lyp-stableopt.chk %chk=O3-T-ub3lyp-opt.chk %mem=60GB %nprocshared=20 # opt freq ub3lyp def2svp guess=read |
sobereva 发表于 2023-6-25 04:12 老师您好,下面这句话我有点没看懂。重新编辑一个新输入文件是什么意思呢,是加一个link吗老师 “现在开始计算自旋极化单重态·CH2-CH2-CH2-CH2·波函数。可以重新编辑一个新输入文件,也可以直接把上面的输入文件的route section改成 # ub3lyp/6-31g(d) guess=read nosymm 并进行运算(其它的都无需改动)。guess=read说明计算时从chk文件中读取片段组合波函数作为初猜,nosymm让程序在迭代过程中允许波函数对称性的改变。迭代结束后,会看到Mulliken自旋布居在两边的碳上分别接近1和-1,说明这正是所希望的态,从Multiwfn绘制的自旋密度的等值面图上能更直观地验证(等值面0.025,绿色和蓝色分别代表正值和负值):” |
|
看我的博文就足够了 谈谈片段组合波函数与自旋极化单重态 http://sobereva.com/82 没必要看官网上说的 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-24 09:26 , Processed in 0.192686 second(s), 25 queries , Gzip On.