计算化学公社

标题: orca6.1.0计算配合物S1和T1的MECP关键词请教 [打印本页]

作者
Author: yzh 时间: 2026-4-29 19:31
标题: orca6.1.0计算配合物S1和T1的MECP关键词请教
本帖最后由 yzh 于 2026-4-30 12:21 编辑

请教一下各位老师：
orca6.1.0计算配合物S1和T1的MECP关键词如何输入，尝试多次都提示关键词存在错误。特别是在%mecp部分调试了好久都是在这里报错。
关键词具体如下：
PBE0 D3BJ def2-SVP RIJCOSX AutoAux TightSCF SurfCrossOpt
%pal nprocs 16 end
%maxcore 2000
%basis
newGTO Cu "def2-TZVP" end
end
%mecp
State1 1 1
State2 3 1
end
%tddft
nroots 5
tda false
end。
* XYZFILE 0 1 cu5-S1.xyz

其中%mecp部分尝试如下关键词也都是Unknown identifier in MECP block line
%surfcross multA 1 rootA 1 # S1 multB 3 rootB 1 # T1 end
%mecp MultA 1 RootA 1 # S1 MultB 3 RootB 0 # T1end
%mecp StatesA {1, 1, 1} (S1) StatesB {3, 1, 1} (T1)end

%surfcross StatesA {1, 1, 1.0} StatesB {3, 1, 1.0}end

谢谢！

作者
Author: wzkchem5 时间: 2026-4-30 11:09
但凡说提示存在错误，必须上传完整的报错信息。你都拿到报错信息了，但是故意不给大家看，这不是凭空增加大家的工作量吗

作者
Author: yzh 时间: 2026-4-30 12:14

wzkchem5 发表于 2026-4-30 11:09
但凡说提示存在错误，必须上传完整的报错信息。你都拿到报错信息了，但是故意不给大家看，这不是凭空增加大 ...

王老师，不好意思，谢谢！

作者
Author: wzkchem5 时间: 2026-4-30 13:34

yzh 发表于 2026-4-30 12:14
王老师，不好意思，谢谢！

报错信息说得很清楚了
Unknown identifier in MECP block line 8 :
Last token : STATE1
所以报错和STATE1有关，应当搜索手册MECP章节有关STATE1关键词的讲解，并仔细对照你的写法和手册的写法的不同，精确到标点符号等

作者
Author: yzh 时间: 2026-4-30 13:42
本帖最后由 yzh 于 2026-4-30 17:05 编辑

wzkchem5 发表于 2026-4-30 13:34
报错信息说得很清楚了
Unknown identifier in MECP block line 8 :
Last token : STATE1

谢谢王老师指导，后来我看了MECP内容，重新修改如下，能跑了。具体如下：

! PBE0 D3BJ def2-SVP RIJCOSX AutoAux TightSCF SurfCrossOpt
%pal nprocs 16 end
%maxcore 2000
%basis
newGTO Cu "def2-TZVP" end
end
%mecp
Mult 3
end
* XYZFILE 0 1 cu5-MECP-S1.xyz
请教一下，以S1-opt结构作为输入坐标，Mult 3，是否就是S1和T1的MECP？如果我需要S0和T1的MECP就是输入坐标为S0-opt的结构么？
谢谢！

作者
Author: wzkchem5 时间: 2026-5-1 16:03

yzh 发表于 2026-4-30 13:42
谢谢王老师指导，后来我看了MECP内容，重新修改如下，能跑了。具体如下：

! PBE0 D3BJ def2-SVP RIJCO ...

你这样写是S0和T1的MECP。注意，算的是哪个态，和你的初始结构是哪个态的平衡结构，没有丝毫关系。程序不知道你的结构是不是S1优化的结构，优化的到底是哪个态和T1的MECP，完全取决于输入文件其他部分怎么写。你要优化S1和T1的MECP，势必有和计算S1相关的关键词，也就是有%tddft block。再仔细看手册

作者
Author: yzh 时间: 2026-5-1 16:11

wzkchem5 发表于 2026-5-1 16:03
你这样写是S0和T1的MECP。注意，算的是哪个态，和你的初始结构是哪个态的平衡结构，没有丝毫关系。程序不 ...

非常感谢，我再去仔细看一下手册。

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