计算化学公社

 找回密码 Forget password
 注册 Register
Views: 9181|回复 Reply: 13
打印 Print 上一主题 Last thread 下一主题 Next thread

[ORCA] 如何用ORCA模拟固态发光的单三态之间的SOC矩阵元?

[复制链接 Copy URL]

131

帖子

0

威望

927

eV
积分
1058

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
各位老师好!
      最近在学习用Gaussian中的oniom模拟固态发光,高层的QM区域就一个分子,低层的MM层有很多同种分子。
现在想在考虑MM层和QM之间相互作用的前提下计算Sn和Tn之间的旋轨耦合矩阵元,请问可以用ORCA算吗?r如果可以的话,请问怎么处理MM区域的分子呢?

490

帖子

2

威望

5306

eV
积分
5836

Level 6 (一方通行)

2#
发表于 Post on 2019-11-26 19:44:26 | 只看该作者 Only view this author
高斯的话可以把MM部分当做背景电荷的方式来考虑,具体请参考charge关键词。
我记得ORCA也有类似于Gaussian的Charge这样的关键词,搜一搜手册便知。

131

帖子

0

威望

927

eV
积分
1058

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2019-11-26 21:21:09 | 只看该作者 Only view this author
让你变成回忆 发表于 2019-11-26 19:44
高斯的话可以把MM部分当做背景电荷的方式来考虑,具体请参考charge关键词。
我记得ORCA也有类似于Gaussian ...

嗯嗯,可以做就行,我翻翻手册,谢啦!

131

帖子

0

威望

927

eV
积分
1058

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2019-11-28 21:44:06 | 只看该作者 Only view this author
本帖最后由 暖空 于 2019-11-28 21:46 编辑
让你变成回忆 发表于 2019-11-26 19:44
高斯的话可以把MM部分当做背景电荷的方式来考虑,具体请参考charge关键词。
我记得ORCA也有类似于Gaussian ...

嗨,您好!
在ORCA中把MM层设为点电荷计算QM区域SOC矩阵元的计算我跑出来了,但是不清楚自己的做法是否正确,所以又来麻烦您啦。

我查找了ORCA 4.1.2 的手册在里面找到了关于点电荷设置的描述,并且做了相关的SOC计算,可以正常结束。

我翻看了ORCA的输出,运算过程中也计算了MM区,感觉应该没啥问题。不过在手册中写到设置点电荷的方式有两种(计算的原理不太一样):

第一种,如果设置的点电荷比较少的话,直接在坐标后面加点电荷和相应的坐标;
第二种,适合上千的点电荷,需要设置一个额外的电荷文件。但是我的理解是这种额外的电荷文件适合的是QM软件结合MM软件计算的QM/MM。

所以,虽然我的MM层也有两千多个原子,但是我还是用的第一种设置点电荷的方式(手册中说的适合设置较少点电荷的那种)。

手册中的后两种设置方法我也都尝试了,都报错了,我还单独计算了QM区域的SOC矩阵元(没设MM层的电荷),和设置了MM层的电荷的结果几乎没差别...

ORCA 4.1.2中关于设置点电荷的描述和我把MM层电荷直接加到QM坐标后面的输出已上传,麻烦大神有时间帮我瞄一眼我的做法和理解对不对,谢谢啦!
ORCA_point charges.pdf (267.88 KB, 下载次数 Times of downloads: 59)
orca-soc.rar (145.91 KB, 下载次数 Times of downloads: 69)

490

帖子

2

威望

5306

eV
积分
5836

Level 6 (一方通行)

5#
发表于 Post on 2019-11-29 08:27:24 | 只看该作者 Only view this author
暖空 发表于 2019-11-28 21:44
嗨,您好!
在ORCA中把MM层设为点电荷计算QM区域SOC矩阵元的计算我跑出来了,但是不清楚自己的做法是否 ...

不好意思,在ORCA中使用charge我也没做过,当初也是简单浏览过手册,没有做实际的测试。

131

帖子

0

威望

927

eV
积分
1058

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2019-11-29 08:46:30 | 只看该作者 Only view this author
让你变成回忆 发表于 2019-11-29 08:27
不好意思,在ORCA中使用charge我也没做过,当初也是简单浏览过手册,没有做实际的测试。

嗯嗯,那我再找找资料,感谢回复

7

帖子

0

威望

409

eV
积分
416

Level 3 能力者

7#
发表于 Post on 2021-8-29 22:13:18 | 只看该作者 Only view this author
你好,我最近在做和你这篇工作一样的计算,你的背景电荷是怎么得到的?

131

帖子

0

威望

927

eV
积分
1058

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2021-8-29 22:52:41 | 只看该作者 Only view this author
wugaiyan 发表于 2021-8-29 22:13
你好,我最近在做和你这篇工作一样的计算,你的背景电荷是怎么得到的?

我的MM层和QM层都是一样的分子,我是用RESP拟合了一个分子的原子电荷,然后用sob老师写的setcharge小程序再把电荷赋给每个分子,你搜一下RESP,sob老师写过一个帖子

7

帖子

0

威望

409

eV
积分
416

Level 3 能力者

9#
发表于 Post on 2021-8-30 08:55:46 | 只看该作者 Only view this author
暖空 发表于 2021-8-29 22:52
我的MM层和QM层都是一样的分子,我是用RESP拟合了一个分子的原子电荷,然后用sob老师写的setcharge小程序 ...

好的  谢谢

68

帖子

0

威望

2313

eV
积分
2381

Level 5 (御坂)

10#
发表于 Post on 2021-11-2 17:32:50 | 只看该作者 Only view this author
你好,问下我用类似方法加入了背景电荷,在做结构优化的时候,报错,Unknown element symbol: 'Q'
这个问题你遇到过吗?

1万

帖子

0

威望

9886

eV
积分
22130

Level 6 (一方通行)

11#
发表于 Post on 2021-11-2 18:31:00 | 只看该作者 Only view this author
annaqz 发表于 2021-11-2 10:32
你好,问下我用类似方法加入了背景电荷,在做结构优化的时候,报错,Unknown element symbol: 'Q'
这个问 ...

如果要带着背景电荷做结构优化,光用点电荷是不行的,必须用QM/MM,因为必须包括QM区和MM区之间的非静电相互作用。可以参见orca说明书的QM/MM章节
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

68

帖子

0

威望

2313

eV
积分
2381

Level 5 (御坂)

12#
发表于 Post on 2021-11-2 18:59:47 | 只看该作者 Only view this author
wzkchem5 发表于 2021-11-2 18:31
如果要带着背景电荷做结构优化,光用点电荷是不行的,必须用QM/MM,因为必须包括QM区和MM区之间的非静电 ...

好的,非常感谢。其实我是之前也翻阅了手册里QM/MM的一些资料,但是感觉看完有点迷糊。我现在想算一个团簇,中心分子2000多个原子,环境分子800多个原子。按照http://bbs.keinsci.com/thread-20960-1-1.html计算得到了CHELPG电荷。然后我想对这个体系进行结构优化。能否采用xtb方法呢?这个关键词要怎么写呢?!QMMM Opt xtb
%qmmm ?我找了一下手册中xtb部分,没看到相关的说明。

1万

帖子

0

威望

9886

eV
积分
22130

Level 6 (一方通行)

13#
发表于 Post on 2021-11-2 19:40:02 | 只看该作者 Only view this author
annaqz 发表于 2021-11-2 11:59
好的,非常感谢。其实我是之前也翻阅了手册里QM/MM的一些资料,但是感觉看完有点迷糊。我现在想算一个团 ...

这个中心区太大了,与其选这么大的中心区域然后用xtb,不如选200个原子左右的区域用DFT,外面用xtb
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

68

帖子

0

威望

2313

eV
积分
2381

Level 5 (御坂)

14#
发表于 Post on 2021-11-2 20:44:14 | 只看该作者 Only view this author
wzkchem5 发表于 2021-11-2 19:40
这个中心区太大了,与其选这么大的中心区域然后用xtb,不如选200个原子左右的区域用DFT,外面用xtb

好的,那我试一下。非常感谢!!

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-23 23:20 , Processed in 0.173128 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list