计算化学公社

标题: 用dalton基于原子平均场近似法计算旋轨耦合常数输入文件怎么输入 [打印本页]

作者
Author: weixue 时间: 2021-12-10 08:13
标题: 用dalton基于原子平均场近似法计算旋轨耦合常数输入文件怎么输入
大家好！想问一下如果想用dalton基于原子平均场近似法计算S1和T1态之间的旋轨耦合常数，应该如何计算，输入文件怎么输入

作者
Author: sobereva 时间: 2021-12-10 09:25
此例给出S1~3与T1~3之间的SOC矩阵元

**DALTON INPUT
.RUN RESPONS
**INTEGRALS
.MNF-SO
**WAVE FUNCTIONS
.DFT
B3LYP
**RESPONS
*QUADRATIC
.DOUBLE RESIDUE
.ISPABC
1 0 1
.PROPRT
X1MNF-SO
.PROPRT
Y1MNF-SO
.PROPRT
Z1MNF-SO
.ROOTS
3
**END OF DALTON INPUT

作者
Author: weixue 时间: 2021-12-10 09:45

sobereva 发表于 2021-12-10 09:25
此例给出S1~3与T1~3之间的SOC矩阵元

**DALTON INPUT

您好！我刚才用这个试了一下，作业直接掉了，没有结果是什么问题呢

作者
Author: sobereva 时间: 2021-12-10 10:03

weixue 发表于 2021-12-10 09:45
您好！我刚才用这个试了一下，作业直接掉了，没有结果是什么问题呢

看输出信息判断

作者
Author: weixue 时间: 2021-12-10 10:09
是这样的报错信息

作者
Author: weixue 时间: 2021-12-10 10:10

sobereva 发表于 2021-12-10 10:03
看输出信息判断

报错信息回复在楼下了

作者
Author: weixue 时间: 2021-12-10 10:31

sobereva 发表于 2021-12-10 10:03
看输出信息判断

您好！我看了报错信息，是下面这样的，显示对H原子没有定义

作者
Author: abin 时间: 2021-12-10 18:32
错误在@Reason行，
需要修改源码，重新编译。
具体看dalton论坛。
我有处理好的。

不过时间有点久，不记得修改了哪里。

我记得手册有说明。

作者
Author: sobereva 时间: 2021-12-11 04:01

weixue 发表于 2021-12-10 10:31
您好！我看了报错信息，是下面这样的，显示对H原子没有定义

根本就不是这个问题

提示里写得很清楚，超过了maxorbs所致
自己改源代码把maxorbs有关的变量/数组改大，重新编译。搞不清楚就google搜

作者
Author: weixue 时间: 2022-4-28 08:33

sobereva 发表于 2021-12-11 04:01
根本就不是这个问题

提示里写得很清楚，超过了maxorbs所致

您好！请问我的分子有金属，那用您上面给的例子，是不是应该加上Douglas-Kroll这个关键词，但是作业迟迟没有结果文件，老师，能给一个具体模板吗

作者
Author: beefly 时间: 2022-4-28 16:27

weixue 发表于 2022-4-28 08:33
您好！请问我的分子有金属，那用您上面给的例子，是不是应该加上Douglas-Kroll这个关键词，但是作业迟迟 ...

你对相对论计算中的很多坑完全不了解。如果仅是要个数据，没必要在这个问题上耗下去，用赝势基组+SO经验参数最方便了，只要是Ru有28个芯电子的赝势，比如def2系列，默认的SO参数就可以用

作者
Author: weixue 时间: 2022-4-28 18:42

beefly 发表于 2022-4-28 16:27
你对相对论计算中的很多坑完全不了解。如果仅是要个数据，没必要在这个问题上耗下去，用赝势基组+SO经验 ...

您好！我现在只想要B3LYP/ANO-DK3/CC-PVDZ方法的计算结果，请问dal文件应该怎么输入呢

作者
Author: beefly 时间: 2022-4-30 22:56

weixue 发表于 2022-4-28 18:42
您好！我现在只想要B3LYP/ANO-DK3/CC-PVDZ方法的计算结果，请问dal文件应该怎么输入呢

根据出错提示改源代码，取消各种上限。ANO-DK3和STO-nG差不多，属于最小基组，这种计算毫无意义

作者
Author: 张晓婷 时间: 2023-4-26 19:21
本帖最后由张晓婷于 2023-4-26 19:22 编辑

sobereva 发表于 2021-12-10 09:25
此例给出S1~3与T1~3之间的SOC矩阵元

**DALTON INPUT

你好，sob老师，我想计算S1-5 到T1-5的旋轨耦合矩阵元，用您这个命令之后把.ROOTS改成5就可以了吗？提交的构型是优化过的S1对吗？加溶剂模型pcm是这么写吗？写的位置有要求吗？
*PCM
.SOLVNT
ACETONIT
*PCMCAV
.NEQRSP

老师，这个命令.NEQRSP有必要加吗？命令的意思是：通过关键词.NEQRSP，我们指定使用溶剂对响应计算的非平衡贡献。
谢谢sob老师

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