计算化学公社

标题: 关于使用UDFT方法计算对称破缺态能量的问题 [打印本页]

作者
Author: 林守龙 时间: 2024-5-21 21:44
标题: 关于使用UDFT方法计算对称破缺态能量的问题
各位老师好，对于高斯UDFT计算自旋极化单重态的能量问题，采用了两种方法，为什么计算出的能量大小不同，请问问题出在了哪里呢？
方法（1）：# stable=opt ub3lyp guess=mix def2tzvp 获得结构的波函数文件.chk，然后%oldchk=scan.chk，# nosymm guess=read def2tzvp um062x计算能量，结果是-499.493383 hartree.
方法（2）：# um062x guess=mix stable=opt def2tzvp直接计算能量，结果是-499.500072 hartree. 明显低于方法（1）算出的能量。
请问两个方法的结果相差这么大(4.2kcal/mol) 的原因是什么？问题出在了哪里？

作者
Author: zjxitcc 时间: 2024-5-21 21:59
在（1）中，若前一步写了nosymm，这一步读取轨道时也要写nosymm；在UM062X计算中要再次检验波函数稳定性。

作者
Author: 林守龙 时间: 2024-5-21 22:29

zjxitcc 发表于 2024-5-21 21:59
在（1）中，若前一步写了nosymm，这一步读取轨道时也要写nosymm；在UM062X计算中要再次检验波函数稳定性。

请问老师nosymm关键词的作用是什么？另外没看懂您说的意思，（1）我在获得波函数时没有写nosymm，在计算能量时我用了nosymm，这样还要检验波函数的稳定性吗？

作者
Author: sobereva 时间: 2024-5-21 23:42

林守龙发表于 2024-5-21 22:29
请问老师nosymm关键词的作用是什么？另外没看懂您说的意思，（1）我在获得波函数时没有写nosymm，在计算 ...

你之前是B3LYP波函数，M06-2X读取它收敛的波函数当初猜不保证这次能收敛到M06-2X的稳定波函数，拿不准就重新做稳定性测试

作者
Author: 林守龙 时间: 2024-5-22 01:12

sobereva 发表于 2024-5-21 23:42
你之前是B3LYP波函数，M06-2X读取它收敛的波函数当初猜不保证这次能收敛到M06-2X的稳定波函数，拿不准就 ...

好的谢谢老师

欢迎光临计算化学公社 (http://bbs.keinsci.com/)