计算化学公社

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

[新手求助] 求助计算MECP的方法

[复制链接 Copy URL]

32

帖子

0

威望

319

eV
积分
351

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
大家好,我最近看文献,看到计算最小能量交叉点用的是Harvey's algorithm的ORCA程序,我知道sob老师编译了一个计算MECP的程序,我是计算磷光的,我想知道两种程序在计算S0和T态势能交叉点的能量有没有差别?希望大家能为我解答疑惑,谢谢大家。

1

帖子

0

威望

195

eV
积分
196

Level 3 能力者

21#
发表于 Post on 2024-1-18 11:44:02 | 只看该作者 Only view this author
Flower 发表于 2016-9-6 15:17
谢谢sob老师,辛苦您了,我有个问题想请教您。老师,我使用sobMECP得到的MECP的能量跟文献上使用ORCA程序 ...

您好,这篇文献可以分享一下吗?

21

帖子

0

威望

781

eV
积分
802

Level 4 (黑子)

20#
发表于 Post on 2023-6-29 09:56:33 | 只看该作者 Only view this author
本帖最后由 liushu 于 2023-6-29 14:00 编辑

已解决

32

帖子

0

威望

319

eV
积分
351

Level 3 能力者

19#
 楼主 Author| 发表于 Post on 2016-9-12 08:07:11 | 只看该作者 Only view this author

谢谢sob老师。

6万

帖子

99

威望

5万

eV
积分
120424

管理员

公社社长

18#
发表于 Post on 2016-9-11 10:02:55 | 只看该作者 Only view this author
Flower 发表于 2016-9-10 08:39
sob老师,祝您教师节快乐。谢谢sob老师,老师,TDE=5.d-5,5.d-5代表多少值呢。

5*10^-5,即5E-5
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

32

帖子

0

威望

319

eV
积分
351

Level 3 能力者

17#
 楼主 Author| 发表于 Post on 2016-9-10 08:39:56 | 只看该作者 Only view this author
sobereva 发表于 2016-9-10 01:56
parameter的括号里这些就是默认值。你要改默认值就直接改这里的

sob老师,祝您教师节快乐。谢谢sob老师,老师,TDE=5.d-5,5.d-5代表多少值呢。

6万

帖子

99

威望

5万

eV
积分
120424

管理员

公社社长

16#
发表于 Post on 2016-9-10 01:56:44 | 只看该作者 Only view this author
Flower 发表于 2016-9-9 21:22
谢谢sob老师,我想了一下午还是不会,老师,我也看到您给别人回复的改变源代码使能量的收敛标准总为 ...

parameter的括号里这些就是默认值。你要改默认值就直接改这里的
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

32

帖子

0

威望

319

eV
积分
351

Level 3 能力者

15#
 楼主 Author| 发表于 Post on 2016-9-9 21:22:43 | 只看该作者 Only view this author
sobereva 发表于 2016-9-9 12:15
就是改.lt.后面的阈值就行了
.lt.是less than的意思

谢谢sob老师,我想了一下午还是不会,老师,我也看到您给别人回复的改变源代码使能量的收敛标准总为真,可我看源代码有五个项,不知道是都要改还是只给其中一个;老师,lt是less than, 可是后面的TGMAX(例如)代表受力最大值,那么sobmecp设置默认的最大值是多少呢?(老师,五个值我其实都想知道默认是多少);最后一个问题,PARAMETER (TDE=5.d-5,TDXMax=4.d-3,TDXRMS=2.5d-3,TGMax=7.d-4,TGRMS=5.d-4) 括号里面代表是什么意思,希望老师能够解答我的疑惑。谢谢老师。真是辛苦您了。

6万

帖子

99

威望

5万

eV
积分
120424

管理员

公社社长

14#
发表于 Post on 2016-9-9 12:15:50 | 只看该作者 Only view this author
Flower 发表于 2016-9-9 11:27
谢谢sob老师,源代码MECP.f应该是下面这些吧,老师,您的帖子也看了,可是不知道该怎么改才能使得收敛变 ...

就是改.lt.后面的阈值就行了
.lt.是less than的意思
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

32

帖子

0

威望

319

eV
积分
351

Level 3 能力者

13#
 楼主 Author| 发表于 Post on 2016-9-9 11:27:13 | 只看该作者 Only view this author
sobereva 发表于 2016-9-9 01:13
源代码很短,看看就知道。帖子里我记得也提过

谢谢sob老师,源代码MECP.f应该是下面这些吧,老师,您的帖子也看了,可是不知道该怎么改才能使得收敛变得严格一点,希望能得到老师的指点,谢谢老师。
      IF (GMax .lt. TGMax) THEN
          PConv(1) = .true.
          flags(1) = "YES"
      END IF
      IF (GRMS .lt. TGRMS) THEN
          PConv(2) = .true.
          flags(2) = "YES"
      END IF
      IF (DXMax .lt. TDXMax) THEN
          PConv(3) = .true.
          flags(3) = "YES"
      END IF
      IF (DXRMS .lt. TDXRMS) THEN
          PConv(4) = .true.
          flags(4) = "YES"
      END IF
      IF (DE .lt. TDE) THEN
          PConv(5) = .true.
          flags(5) = "YES"
      END IF

6万

帖子

99

威望

5万

eV
积分
120424

管理员

公社社长

12#
发表于 Post on 2016-9-9 01:13:32 | 只看该作者 Only view this author
Flower 发表于 2016-9-8 16:25
谢谢sobl老师,能量差确实比我找的大很多,初猜结构我已经换过几次了。老师,我比较菜,能否告诉我源代码 ...

源代码很短,看看就知道。帖子里我记得也提过
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

32

帖子

0

威望

319

eV
积分
351

Level 3 能力者

11#
 楼主 Author| 发表于 Post on 2016-9-8 16:25:26 | 只看该作者 Only view this author
sobereva 发表于 2016-9-8 11:34
如果文献里MECP和MC结构差得比这明显要大,可能是SobMECP用的收敛标准对当前体系或当前初猜结构不合适导 ...

谢谢sobl老师,能量差确实比我找的大很多,初猜结构我已经换过几次了。老师,我比较菜,能否告诉我源代码的收敛限该怎么改。非常感谢老师。

6万

帖子

99

威望

5万

eV
积分
120424

管理员

公社社长

10#
发表于 Post on 2016-9-8 11:34:32 | 只看该作者 Only view this author
Flower 发表于 2016-9-8 09:16
抱歉,让老师困惑了。我做的是铱配合物磷光材料,MC态是mental-centered,我理解的是自旋密度大部分分布 ...

如果文献里MECP和MC结构差得比这明显要大,可能是SobMECP用的收敛标准对当前体系或当前初猜结构不合适导致提前判断为找到MECP了,可以尝试改代码调节收敛限,或者试试别的初猜结构。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

32

帖子

0

威望

319

eV
积分
351

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2016-9-8 09:16:32 | 只看该作者 Only view this author
sobereva 发表于 2016-9-7 23:56
不清楚你的意思,MC态指什么,MECP要求哪两个态之间的?什么叫能量差太小?要同时看结构。

ORCA手册里 ...

抱歉,让老师困惑了。我做的是铱配合物磷光材料,MC态是mental-centered,我理解的是自旋密度大部分分布在金属中心,MECP求的是三态势能面和基态的交叉点,失活过程,T1 -> ts -> mc -> mecp -> So。能量差太小指的是MECP与MC之间的能量差,我得到的差只有0.1kcal/mol,离我想要的1.7kcal/mol差的远,构型的话键长最大拉伸了0.03埃。老师,我的想法是如果能通过sobMECP计算得到我想要的结果就不用去用ORCA计算了,所以才想知道是不是我哪里可以改进然后再计算。谢谢老师了。

6万

帖子

99

威望

5万

eV
积分
120424

管理员

公社社长

8#
发表于 Post on 2016-9-7 23:56:52 | 只看该作者 Only view this author
Flower 发表于 2016-9-7 20:54
谢谢sob老师。老师,我还有几个问题想请教您。(1)我用sobmecp计算得到的MECP的能量与MC态之间的能量差 ...

不清楚你的意思,MC态指什么,MECP要求哪两个态之间的?什么叫能量差太小?要同时看结构。

ORCA手册里就有其MECP功能的描述和例子。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

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

GMT+8, 2025-8-27 16:54 , Processed in 0.289687 second(s), 21 queries , Gzip On.

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