| MOLPRO extractor (ME) |
功能
| 1. 从一组MOLPRO自旋-轨道耦合（SOC）计算的输出文件中，提取Ω态的能量，打印到omiga.dat输出文件中。
2. 对奇数个电子的线型分子SOC计算，判断MOLPRO产生的Ω态的对称性。由于MOLPRO只能用C2v和D2h群处理线型分子，所以SOC计算输出的Ω对称性为E1/2（C2v）或E1/2_g和E1/2_u（D2h）。ME程序可以判断C2v双点群下的对称性（D2h目前暂不支持）。
3. 从一组MOLPRO的MULTI计算输出文件中，提取MCSCF能量，打印到multi.dat输出文件中。
4. 从一组MOLPRO的MRCI或CASPTn计算输出文件中，提取MCSCF零级参考态能量，打印到ciref.dat输出文件中。
5. 从一组MOLPRO的MRCI/CASPT2/CASPT3计算输出文件中，提取MRCI/MRCI+Q/CASPT2/CASPT3的能量，打印到ci.dat文件中。
6. 从一组MOLPRO自旋-轨道耦合（SOC）计算的输出文件中，提取Omega态之间的跃迁偶极矩（TDM）的六个分量（xyz分量，包括实部和虚部），把TDM六个分量以及TDM的大小打印到tdmomega.dat输出文件中。输出的TDM未考虑电子态的简并度。因此对于简并态间的跃迁，需按以下公式计算总TDM：
d1,d2
TDM_tot ** 2 = Σ (TDM_i,j ** 2) / d1
i,j
d1,d2：上、下态简并度；i=1,...,d1; j=1,...,d2；程序打印的TDM对应公式的TDM_i,j。
7. 从一组MOLCAS的CASPT2计算输出文件中，提取CASPT2或MS-CASPT2能量，打印到caspt.dat输出文件中。
8. 所有的.dat数据输出文件均采用逗号为分隔符，可以用MS EXCEL、ORIGIN打开后画图。 |
| 输入说明 | 参见使用手册 |
注意事项
| 目前，Ω态对称性的判断存在以下限制：
1. 必须是从C2v映射到C∞v，且必须是奇数个电子的线型分子体系。D2h-->D∞h目前暂不支持；
2. 仅支持二重和四重ΛS态。
3. 轨道角动量支持Σ-，Σ+，Π，Δ，Ф，Г。因此目前可能得到的最大Ω=11/2（来自4Г）。MOLPRO的SOC计算中不能含有其它轨道量子数，否则出错。
4. 最多处理50个Ω态（若不考虑二度简并，为100个），若超过则ME报错。这是因为MOLPRO只能打印2位整数的态序号。要处理更多的态，需要先修改MOLPRO源代码，或者先对MOLPRO SOC的输出文件进行二次处理（人工或另写程序），然后修改ME的源代码。 |
| 下载 | 源代码 |
| | 可执行程序（Windows）和示例 |
| 运行方法（Win） | |
| | 1. 编辑批处理文件“run.bat”，根据需要，可修改其中的输入文件名。
2. 编辑输入文件。
3. 用鼠标双击run.bat，运行程序。
4. 计算结果位于“*.dat”文件中。例如，test1.inp的输出为：
Title: TlBr, CASPT2 + SO-CASSCF (8e,8MO), BS: ECP
E --- cm^-1
R, A1 , A1 , A1 , ...
1.800000000, 0.00, 20665.45, 24749.45, ...
1.900000000, 0.00, 21605.05, 25599.51, ... |
test3.inp的输出为： ****** L-S states ******
~~~~~~~~~~~~~~~~~~~~~~~~
2A1:
2De 2De 2S+ 2S+ 2De 2Ga 2S+ 2De 2S+
2B1:
2Pi 2Fi 2Pi 2Pi 2Pi 2Fi 2Pi
2A2:
2De 2S- 2De 2De 2S- 2Ga 2De
4A1:
4De
4B1:
4Pi 4Fi 4Pi
4A2:
4De 4S- 4S-
***** Omega states *****
~~~~~~~~~~~~~~~~~~~~~~~~
1/2 (*18):
3 9 11 17 25 31 39 45 57 59 61 69 73 77 79 83 85 99
3/2 (*15):
1 7 13 21 27 37 41 49 55 63 67 71 75 81 89
5/2 (*10):
5 15 23 35 47 51 53 65 87 93
7/2 (* 5):
19 33 43 91 97
9/2 (* 2):
29 95 |
|
升级：
| 2005.10.21
1. 修改了拼写错误（omiga --> omega） ；
2. 删除了$omegasymm选项组的nom选项，由程序自动计算该值。2006.1.5
1. 提取MOLPRO的CI/CI+Q计算结果；
2. 添加对MOLCAS计算结果的支持，为此把$job选项组的itype选项由数值改为数组，进行控制；
3. 提取MOLCAS的CASPT2/MS-CASPT2计算结果。 2006.3.8
1. 提取MOLPRO计算的CASPT2（来自rs2c/rs2/rs3模块）/CASPT3（来自rs3模块）能量。 2006.5.14
1. 提取MOLPRO计算的Ω态之间的跃迁偶极矩。 2007.2.13
1. 提取MOLPRO计算的AQCC，ACPF能量。 2007.7.1
1. 提取MOLPRO计算的MCSCF轨道因子。
1. 提取DIRAC 04计算的轨道因子。
3. 由于Intel Fortran与GNU g95编译器不兼容，导致对输入文件的格式要求不同。为避免此问题，$omegasymm选项组的逻辑型变量lsprint改为整型lprint。
4. 为了输入格式统一，在Omega态对称性的计算中，输入文件改由$files选项组控制，并且每次只分析一个。
5. 修正了个别语法错误。 2007.7.14
1. 提取GAMESS-US或PC-GAMESS计算的MCQDPT2能量。 2007.11.15
1. 对于奇数个电子的线型分子（z为分子轴），打印垂直方向和平行方向的跃迁偶极矩。 |
|