计算化学公社

标题: 求助[FeCl4]2- L-edge X射线吸收谱NDO分析 [打印本页]

作者
Author: conanatm 时间: 2019-7-26 21:55
标题: 求助[FeCl4]2- L-edge X射线吸收谱NDO分析
大家好。我正在按照ORCA的手册,使用ROCIS重复[FeCl4]2- L-edge X射线吸收谱的计算（9.23.2, "Natural Transition Orbitals / Natural Difference Orbitals")。手册中给出了NDO分析，激发态14对应2pz到3dyz的跃迁（符合L-edge)。我的计算可以重复NDO46[1]，但是NDO49[0]看起来也是d轨道（应是p轨道）。希望大家可以提些建议。
此外我注意到Multiwfn并不支持读取ORCA的ROCIS输出文件，读取molden文件提示找不到"excited state"关键词。
此处附上手册中NDO的截图，我得到的NDO，绘图参数，输入文件，以及输出文件的一部分：

手册给出的NDO：

(, 下载次数 Times of downloads: 30)

重复出的NDO46[1]：

(, 下载次数 Times of downloads: 26)

无法重复的NDO49[0]：

(, 下载次数 Times of downloads: 31)

orca_plot绘图参数：

(, 下载次数 Times of downloads: 24)

输入文件：

!B3LYP def2-TZVP def2-TZVP/C (手册中为def2/J) Conv TightSCF LargePrint PAL4

%ROCIS

NRoots 40

PrintLevel 3

MaxCore 4000

MaxDim 360

SOC true

DoRI true

DoNDO true

NDOThresh 1e-4

NDOStates= 1,2,3,4,5,6,7,8,9,10,13,14,15

DoLowerMult true

DoHigherMult true

DoDFTCIS true DFTCIS_c = 0.18, 0.20, 0.40

OrbWin = 6,8,0,2000

end

* xyz -2 5

Fe -17.84299991694815 -0.53096694321123 6.09104775508499

Cl -19.84288422845700 0.31089495619796 7.04101319789001

Cl -17.84298666758073 0.11868125024595 3.81067954087770

Cl -17.84301352218429 -2.87052442818457 6.45826391412877

Cl -15.84311566482982 0.31091516495189 7.04099559201853

部分输出文件：

-------------------------------------------------

NATURAL DIFFERENCE ORBITALS FOR STATE 14

-----------------------------------------------

done

Solving eigenvalue problem for the Occupied space ... done

Solving eigenvalue problem for the Acceptor space ... done

Natural Difference Orbitals were saved in ndo.14.ndo

Threshold for printing occupation numbers 1.0e-04

E= 25.568527 au 695.755 eV 5611643.0 cm**-1 (手册中为692.469 eV 5585137.0 cm**-1)

49[0] -> 46[1] : n= 0.66696849 (手册中为0.81173217)

48[0] -> 47[1] : n= 0.27076736 (手册中为0.17903699)

47[0] -> 48[1] : n= 0.06831675

46[0] -> 49[1] : n= 0.00101007

谢谢大家。

作者
Author: sobereva 时间: 2019-7-26 22:32
“读取molden文件提示找不到"excited state"关键词”这句话意义不明。molden文件里本来就没excited state信息，读这些信息是从输出文件里读的
Multiwfn目前不官方支持ROCIS任务的输出，但是基于此时的molden看轨道是应当没问题的

我没做过这个例子。不排除程序更新后手册没有相应地更新导致对不上

作者
Author: conanatm 时间: 2019-7-26 22:42

sobereva 发表于 2019-7-26 22:32
“读取molden文件提示找不到"excited state"关键词”这句话意义不明。molden文件里本来就没excited state信 ...

是的，可以从ORCA生成的molden文件中读取轨道，但是不能做NTO分析（提示找不到“excited state“关键词）。

作者
Author: conanatm 时间: 2019-7-28 21:19
问题已解决。Avogadro会自动读取在同一目录下的其他输出文件（如图所示），我之前错将MO当成NDO绘图。事实上orca_plot生成的每个cube文件只包含一个NDO轨道，用Avogadro绘制相应的等值面即可。

(, 下载次数 Times of downloads: 17)

我有一个追加问题：NDO是在引入SOC前计算的，如manual中的state 14是ROCIS的本征向量。而SOC会改变激发态数目与能量。那么如果我们对模拟所得谱图中的某个峰感兴趣，我们如何得知应选择哪些激发态进行NDO分析？

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