计算化学公社

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

[ORCA] 关于ORCA计算体系分子动力学的问题

[复制链接 Copy URL]

82

帖子

0

威望

512

eV
积分
594

Level 4 (黑子)

ORCA的手册中提到:simple Berendsen thermostat is implemented in the ORCA MD module,是不是表示在模拟过程中能够控制体系的温度?然后手册中又提到:The cells can be defi ned as elastic, such that their size adapts to the system. This enables to run simulations under constant pressure,表示在模拟过程中能够给定体系的压力。那么是不是表示ORCA能够做NPT系综分子动力学。我想做的是:将一定温度压力下的氖气进行分子动力学模拟,然后再进行TDDFT模拟,计算荧光光谱。老师们能给些建议么?比如到底应该用哪个软件做我这种模型的分子动力学比较好,谢谢!

6万

帖子

99

威望

5万

eV
积分
120159

管理员

公社社长

2#
发表于 Post on 2019-3-4 16:22:51 | 只看该作者 Only view this author
NPT可以,但是模拟你目前这个问题应该没戏。氖气的密度太低了,TDDFT要算也得是凝聚相或者团簇
北京科音自然科学研究中心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

82

帖子

0

威望

512

eV
积分
594

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2019-3-4 17:13:49 | 只看该作者 Only view this author
sobereva 发表于 2019-3-4 16:22
NPT可以,但是模拟你目前这个问题应该没戏。氖气的密度太低了,TDDFT要算也得是凝聚相或者团簇

我的体系是一个温稠密体系,密度大概有1g/cm3左右。

6万

帖子

99

威望

5万

eV
积分
120159

管理员

公社社长

4#
发表于 Post on 2019-3-4 17:22:33 | 只看该作者 Only view this author
huangyoyo 发表于 2019-3-4 17:13
我的体系是一个温稠密体系,密度大概有1g/cm3左右。

我没见过这样体系的TDDFT计算,是否靠谱也很难说。如果要算的话,先看看有没有类似的计算文章
北京科音自然科学研究中心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

82

帖子

0

威望

512

eV
积分
594

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2019-3-5 08:23:41 | 只看该作者 Only view this author
sobereva 发表于 2019-3-4 17:22
我没见过这样体系的TDDFT计算,是否靠谱也很难说。如果要算的话,先看看有没有类似的计算文章

恩恩!不过我们这边有实验数据。另外我想请问算分子动力学是不是跟几何优化一样,对基组的要求不高,只要在相同水平下算频率,保证系统处于基态,最后根据你的“Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法”文章中提到的方法计算荧光光谱。

6万

帖子

99

威望

5万

eV
积分
120159

管理员

公社社长

6#
发表于 Post on 2019-3-5 19:34:46 | 只看该作者 Only view this author
huangyoyo 发表于 2019-3-5 08:23
恩恩!不过我们这边有实验数据。另外我想请问算分子动力学是不是跟几何优化一样,对基组的要求不高,只要 ...

如果你是要做AIMD,本身对计算级别要求不高,但也得至少能定性描述动力学行为(比如你的体系,不能用B3LYP,用的话得加色散校正)
但这种体系用经典力场的MD其实就够了。
北京科音自然科学研究中心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

82

帖子

0

威望

512

eV
积分
594

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2019-3-6 08:48:53 | 只看该作者 Only view this author
sobereva 发表于 2019-3-5 19:34
如果你是要做AIMD,本身对计算级别要求不高,但也得至少能定性描述动力学行为(比如你的体系,不能用B3LY ...

我的研究是一个强耦合的等离子体体系,也需要加色散校正?另外我打算用6-31+G*做分子动力学,看了你的文章“乱谈激发态的计算方法”后,我选择用ma-def2-TZVP(涉及里德堡激发)进行激发态计算,这样的基组选择可行么?

6万

帖子

99

威望

5万

eV
积分
120159

管理员

公社社长

8#
发表于 Post on 2019-3-7 02:59:31 | 只看该作者 Only view this author
huangyoyo 发表于 2019-3-6 08:48
我的研究是一个强耦合的等离子体体系,也需要加色散校正?另外我打算用6-31+G*做分子动力学,看了你的文 ...

DFT-D色散校正加了没坏处,也不增加耗时
动力学用那个级别可以,不过用6-311G*代替往往更省时间,结果还可能更好。对于里德堡激发计算,ma-TZVP只能算是刚踏入门槛,要定量准确的话至少aug-cc-pVTZ。还有余力的话提升到daug-cc-pVTZ。
北京科音自然科学研究中心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

82

帖子

0

威望

512

eV
积分
594

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2019-3-7 09:13:40 | 只看该作者 Only view this author
sobereva 发表于 2019-3-7 02:59
DFT-D色散校正加了没坏处,也不增加耗时
动力学用那个级别可以,不过用6-311G*代替往往更省时间,结果还 ...

ORCA里没有aug-cc-PVTZ合适的辅助基组,是不是不能结合RICOSX来算呢?(体系有大概两百多个原子)

6万

帖子

99

威望

5万

eV
积分
120159

管理员

公社社长

10#
发表于 Post on 2019-3-8 05:48:50 | 只看该作者 Only view this author
huangyoyo 发表于 2019-3-7 09:13
ORCA里没有aug-cc-PVTZ合适的辅助基组,是不是不能结合RICOSX来算呢?(体系有大概两百多个原子)

有自带的辅助基组

北京科音自然科学研究中心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

82

帖子

0

威望

512

eV
积分
594

Level 4 (黑子)

11#
 楼主 Author| 发表于 Post on 2019-3-8 09:24:52 | 只看该作者 Only view this author
本帖最后由 huangyoyo 于 2019-3-8 09:33 编辑
sobereva 发表于 2019-3-8 05:48
有自带的辅助基组

惰性气体轨道之间的耦合是根据JK耦合的,跃迁规则为deta_K=deta_J=-1,0,1。在做完分子动力学之后我还想算Ne的18.72638130对应的那个能级(如果按照LS耦合来看这对应的是一个单重态P轨道)到16.61906995对应的那个能级(三重态J=2的P轨道)的发射能,那么用ORCA能算么?如果假设用100个原子的话,18.72638130对应的那个态大概都是TDDFT结果中10000之后的某个态了,光想想都觉得不可行呢。

Ne.png (72.69 KB, 下载次数 Times of downloads: 42)

Ne.png

82

帖子

0

威望

512

eV
积分
594

Level 4 (黑子)

12#
 楼主 Author| 发表于 Post on 2019-3-15 15:06:33 | 只看该作者 Only view this author
sobereva 发表于 2019-3-4 16:22
NPT可以,但是模拟你目前这个问题应该没戏。氖气的密度太低了,TDDFT要算也得是凝聚相或者团簇

ORCA中的DFT-MD考虑电子效应的影响么(电子电离)?还是稳定后结构中的电子还是处于基态?请问能推荐一些关于将DFT-MD的书或者文献么?

6万

帖子

99

威望

5万

eV
积分
120159

管理员

公社社长

13#
发表于 Post on 2019-3-15 16:48:07 | 只看该作者 Only view this author
huangyoyo 发表于 2019-3-15 15:06
ORCA中的DFT-MD考虑电子效应的影响么(电子电离)?还是稳定后结构中的电子还是处于基态?请问能推荐一些 ...

除非你指定的是用计算激发态的方法,或者自旋多重度设成激发态的自旋多重度,否则MD过程中一直都在基态上
北京科音自然科学研究中心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

82

帖子

0

威望

512

eV
积分
594

Level 4 (黑子)

14#
 楼主 Author| 发表于 Post on 2019-3-18 10:29:51 | 只看该作者 Only view this author
sobereva 发表于 2019-3-15 16:48
除非你指定的是用计算激发态的方法,或者自旋多重度设成激发态的自旋多重度,否则MD过程中一直都在基态上

您好,我在一篇计算温稠密物质光学性质的文章中看到如下几句话:The electron system is described using the ab initio plane wave code ABINIT.
To incorporate effects of the ion system at elevated temperatures, molecular dynamics (MD) simulation was also carried out. All simulation results presented in the present work were obtained using 108 atoms in a cubic simulation box with periodic boundary conditions. For two-temperature MD simulations the ion and electron temperatures were controlled independently. While the electron temperature was fixed for every timestep in the simulation, the ion temperature was controlled by a Nos′e-Hoover thermostat.
看这字面意思,好像这个分子动力学能够考虑电离以及电子的行为?

6万

帖子

99

威望

5万

eV
积分
120159

管理员

公社社长

15#
发表于 Post on 2019-3-18 23:56:14 | 只看该作者 Only view this author
huangyoyo 发表于 2019-3-18 10:29
您好,我在一篇计算温稠密物质光学性质的文章中看到如下几句话:The electron system is described using ...

貌似是。我对此程序了解不多
北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2025-8-16 13:06 , Processed in 0.295651 second(s), 23 queries , Gzip On.

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