求助：CAM-B3LYP计算振动分辨的紫外吸收光谱出错

liaorongbao

求助：振动分辨的紫外吸收光谱CAM-B3LYP计算过程出错
依据振动分辨的紫外吸收光谱计算方法(http://sobereva.com/223)，但做了如下调整

第一步：
把计算激发态的# cis(root=1)/6-31G*          opt freq=saveNM
改成了密度泛函# CAM-B3LYP/6-31G* TD(ROOT=1) opt freq=saveNM

改写后完整的输入文件如下：
%rwf=anisole_exc.rwf
%nosave
%mem=30GB
%nprocshared=30
%chk=anisole_exc.chk
# maxdisk=500GB
# CAM-B3LYP/6-31G* TD(ROOT=1) opt freq=saveNM

anisole S1

0 1
C                  2.28445000    0.32691600    0.00003800
C                  1.34343300    1.38142800   -0.00010400
C                 -0.03741900    1.09055200   -0.00003400
C                 -0.44320000   -0.26595700   -0.00005800
C                  0.49751700   -1.32326400   -0.00002000
C                  1.87678500   -1.02037200    0.00017600
H                  3.33003700    0.56133100   -0.00002300
H                  1.67739200    2.39758400   -0.00018600
H                 -0.75550500    1.88022800    0.00018000
H                  0.12880400   -2.32514600   -0.00027400
H                  2.60234200   -1.80568700    0.00035800
O                 -1.73605800   -0.64937600   -0.00027400
C                 -2.82140400    0.30027500    0.00020300
H                 -3.71931200   -0.29382400    0.00075700
H                 -2.78766600    0.92165900    0.88463800
H                 -2.78859800    0.92139000   -0.88445500



在我的工作目录下使用pwd命令显示：home/jinfeng。(在这里我之前写错了，应该是/home/jinfeng。这就是出错说找不到文件的原因)
计算结束后在home/jinfeng目录下能看到文件anisole_exc.chk，且可以用gview正常打开anisole_exc.log文件


第二步：
把计算基态的# hf/6-31G*        opt freq=FC nosymm
改成密度泛函# CAM-B3LYP/6-31G* opt freq=FC nosymm

改写后完整的输入文件如下：
%rwf=anisole.rwf
%nosave
%mem=30GB
%nprocshared=30
%chk=anisole.chk
# maxdisk=500GB
# CAM-B3LYP/6-31G* opt freq=FC nosymm

anisole S0

0 1
C                 -2.27936800    0.32609200    0.00008100
C                 -1.34379900    1.33903700   -0.00003400
C                  0.01378100    1.05103000   -0.00016100
C                  0.43663600   -0.26434700   -0.00020100
C                 -0.50250800   -1.28814500   -0.00003400
C                 -1.84702200   -0.99409200    0.00014300
H                 -3.32624700    0.55310300    0.00026000
H                 -1.66092300    2.36324200    0.00002000
H                  0.71962800    1.85453200   -0.00023700
H                 -0.14729200   -2.29705700   -0.00003400
H                 -2.56259500   -1.79229300    0.00028200
O                  1.75151800   -0.65529700   -0.00024100
C                  2.80692300    0.31779600    0.00026800
H                  3.72320000   -0.24943200    0.00040700
H                  2.76733600    0.94326200   -0.88309400
H                  2.76687800    0.94280000    0.88395800

home/jinfeng/anisole_exc.chk(出错原因在这里，应该写成/home/jinfeng/anisole_exc.chk或者./anisole_exc.chk或者anisole_exc.chk。这就是出错说找不到文件的原因)


最后的出错症状如下：
             Generation of the Franck-Condon spectrum

**********************************************************************


     ==================================================
               Information on the Simulation
     ==================================================

Type of Spectroscopy: ONE-PHOTON ABSORPTION         
Model applied to the transition: ADIABATIC HESSIAN
Approx. of the electronic transition dipole moment: FC
Temperature effect are not taken into account.

     ==================================================
                  Treatment of Input Data
     ==================================================
Data for initial state taken from current calculation.
Data for   final state taken from checkpoint file "home/jinfeng/anisole_exc.chk"
Error termination in NtrErr:
ntran open failure returned to fopen.

abdoman

%chk=anisole.chk 错了吧。是不是%chk=anisole_exc.chk

liaorongbao

abdoman 发表于 2019-7-29 10:13
%chk=anisole.chk 错了吧。是不是%chk=anisole_exc.chk

正在使用bomsaude网友的方法重算，查找出错原因。谢谢您的关注。补充：刚刚发现原因了，home/jinfeng/anisole_exc.chk改成anisole_exc.chk就好了

sobereva

liaorongbao 发表于 2019-7-29 10:26
正在使用bomsaude网友的方法重算，查找出错原因。谢谢您的关注。补充：刚刚发现原因了，home/jinfeng/ani ...

不能写home，需要写/home。前者代表是当前目录下的home目录，显然和实际不符

liaorongbao

sobereva 发表于 2019-7-29 22:53
不能写home，需要写/home。前者代表是当前目录下的home目录，显然和实际不符

的确如Sob老师所说！