计算化学公社

标题: 求助：CAM-B3LYP计算振动分辨的紫外吸收光谱出错 [打印本页]

作者
Author: liaorongbao 时间: 2019-7-28 20:19
标题: 求助：CAM-B3LYP计算振动分辨的紫外吸收光谱出错
本帖最后由 liaorongbao 于 2019-7-29 17:06 编辑

求助：振动分辨的紫外吸收光谱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.

作者
Author: abdoman 时间: 2019-7-29 10:13
%chk=anisole.chk 错了吧。是不是%chk=anisole_exc.chk

作者
Author: liaorongbao 时间: 2019-7-29 10:26
本帖最后由 liaorongbao 于 2019-7-29 11:30 编辑

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

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

作者
Author: sobereva 时间: 2019-7-29 22:53

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

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

作者
Author: liaorongbao 时间: 2019-8-10 07:23

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

的确如Sob老师所说！

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