计算化学公社

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

[新手求助] 在跑高温分子动力学时杂原子的手性是否能够保持?

[复制链接 Copy URL]

200

帖子

1

威望

1325

eV
积分
1545

Level 5 (御坂)

尊敬的各位老师好:
  我想用Molclus调用xTB去跑带有三级胺,三级膦,三级硫的分子动力学产生traj.xyz文件用于构象分析,三级胺的翻转能垒比较低,但是三级膦的翻转能垒比较高,室温下三级膦应该是可能出现手性的。而高温下,例如100摄氏度,手性膦将可以自由翻转,所以在跑MD取样的时候,我是不是不能把温度设的很高,不然三级膦就会翻转了,这个与实际室温下不符合?但是另一方面如果温度不设高,会不会采样不完全,导致不一定能找到最优的构象呢?
                                                                                                                                                                                                                 谢谢各位老师!


计算小白 水平很菜
随时可喷 万望指点

1665

帖子

5

威望

4788

eV
积分
6553

Level 6 (一方通行)

喵星人

2#
发表于 Post on 2022-3-1 03:08:23 | 只看该作者 Only view this author
如果你的AIMD不用特殊姿势(增强采样),在你能跑得动的时间尺度内基本上看不到翻转

你可以先用常规静态计算看看翻转能垒的范围,在想象在你模拟的时间尺度内发生的概率

200

帖子

1

威望

1325

eV
积分
1545

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2022-3-1 09:31:53 | 只看该作者 Only view this author
喵星大佬 发表于 2022-3-1 03:08
如果你的AIMD不用特殊姿势(增强采样),在你能跑得动的时间尺度内基本上看不到翻转

你可以先用常规静态计 ...

老师 请问这样的效果是三级胺可以正常翻转,三级膦不会翻转吗?还是想要三级胺翻转的。。。
计算小白 水平很菜
随时可喷 万望指点

1665

帖子

5

威望

4788

eV
积分
6553

Level 6 (一方通行)

喵星人

4#
发表于 Post on 2022-3-1 10:07:55 | 只看该作者 Only view this author
WaterMirror 发表于 2022-3-1 09:31
老师 请问这样的效果是三级胺可以正常翻转,三级膦不会翻转吗?还是想要三级胺翻转的。。。

如果你按照一般的温度正常的构型去用DFT跑MD,那你在你能跑得动的时间尺度内是看不到翻转的,无论是三级胺还是膦

你应该先用正常的量子化学找过渡态的方法先确认势垒的大致高度

200

帖子

1

威望

1325

eV
积分
1545

Level 5 (御坂)

5#
 楼主 Author| 发表于 Post on 2022-3-1 10:38:15 | 只看该作者 Only view this author
喵星大佬 发表于 2022-3-1 10:07
如果你按照一般的温度正常的构型去用DFT跑MD,那你在你能跑得动的时间尺度内是看不到翻转的,无论是三级 ...

哦哦,谢谢老师解答!
那就是应该要选是一个让胺可以翻转但是膦不能翻转的合适温度去跑MD是吗
计算小白 水平很菜
随时可喷 万望指点

1万

帖子

0

威望

9876

eV
积分
22118

Level 6 (一方通行)

6#
发表于 Post on 2022-3-1 16:55:05 | 只看该作者 Only view this author
WaterMirror 发表于 2022-3-1 03:38
哦哦,谢谢老师解答!
那就是应该要选是一个让胺可以翻转但是膦不能翻转的合适温度去跑MD是吗

理论上是的,但是需要注意一点,你跑动力学是ps量级的时间尺度,所以应该选择一个胺用1ps能翻转、膦用1ps不能翻转的温度。所以不能简单查实验文献来用,因为实验一般是测的一个温度下能不能用毫秒、秒甚至分钟量级的时间翻转。
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
本团队长期招收研究生,有意者可私信联系

200

帖子

1

威望

1325

eV
积分
1545

Level 5 (御坂)

7#
 楼主 Author| 发表于 Post on 2022-3-1 22:21:30 | 只看该作者 Only view this author
wzkchem5 发表于 2022-3-1 16:55
理论上是的,但是需要注意一点,你跑动力学是ps量级的时间尺度,所以应该选择一个胺用1ps能翻转、膦用1ps ...

谢谢老师那这要自己摸索了
计算小白 水平很菜
随时可喷 万望指点

6万

帖子

99

威望

6万

eV
积分
125150

管理员

公社社长

8#
发表于 Post on 2022-3-2 03:39:49 | 只看该作者 Only view this author
很多程序可以外加限制势,对锥形部分加上二面角限制势就可以避免出现动力学过程中的翻转
北京科音自然科学研究中心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

1

帖子

0

威望

69

eV
积分
70

Level 2 能力者

9#
发表于 Post on 2022-4-18 18:09:21 | 只看该作者 Only view this author
你好,请问如何计算手性分子的能垒差?

本版积分规则 Credits rule

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

GMT+8, 2026-2-23 06:11 , Processed in 0.180032 second(s), 26 queries , Gzip On.

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