|
|
本帖最后由 likjvau 于 2022-3-3 04:24 编辑
各位老师好，学生在对掺稀土元素的配合物进行几何优化和激发态计算时，出现了激发能是负数的情况，且其余激发能也与实验不匹配
Excited State 1: 4.001-A -0.9545 eV -1298.90 nm f=-0.0000 <S**2>=3.753
首先使用大核赝势MWB57对稀土Er3+进行几何优化，由于赝化了全部4f电子，所以体系的自旋多重度为奇数，几何优化使用电荷和自旋为：0 1
由于小核赝势进行基态优化太困难，经过咨询，了解到，可以使用大核赝势优化基态结构，然后用小核赝势进行激发态计算
但是我通过使用小核赝势MWB28进行单点能计算（自旋只能是偶数），发现自旋为4时单点能最低，所以就认为基态自旋为4，所以直接使用了MWB28和自旋4在原本使用MWB57和自旋1优化出的基态结构上进行了激发态计算
起初在td计算中加上了stable=opt，然后激发态计算从爆出波函数不稳定后，整个计算就报错停了，调研后发现stable=opt不能与td计算一起使用。就把stable=opt去掉了，然后由于激发态计算过程中scf产生的不是基态真正波函数，计算结果中存在负激发能。
学生有疑问：
我是否需要在单点能计算时判断出基态自旋为4后，用mwb28和自旋4重新优化结构？（但是这个基本优化不出来，小核赝势优化非常难受。。）
我是否可以在基态优化使用自旋1，而在激发态计算直接换用小核赝势，并改用自旋4
出现负激发态这个问题是为什么？应该怎么解决？
我的几何优化和激发态计算预测吸收光谱的过程，都出现了哪些错误？
救救孩子吧，辛苦各位老师解答，万分感谢
|
|
发表于 Post on 2022-3-3 04:24:39
|
只看该作者 Only view this author
|只看大图 View big picture
|倒序浏览 Reverse view
|阅读模式 Reading model
收藏 Add to favorites
楼主 Author