计算化学公社

 找回密码 Forget password
 注册 Register

使用Multiwfn计算分子的动力学直径

查看数: 12554 | 评论数: 11 | 收藏 Add to favorites 10
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2019-8-19 20:31

正文摘要:

使用Multiwfn计算分子的动力学直径Using Multiwfn to calculate dynamic diameter of molecules 文/Sobereva@北京科音First release: 2019-Aug-19  Last update: 2024-May-4 1 相关知识&计算原理 ...

回复 Reply

sobereva 发表于 Post on 2024-7-23 03:32:55
chen0201 发表于 2024-7-18 15:19
请问老师,我按照您的方法计算了分子动力学直径,但是对于苯乙炔这种分子来说,动力学直径是用躺着的苯环 ...

躺着的
文中CO2的例子体现出各个方向跨度不同时怎么考虑
chen0201 发表于 Post on 2024-7-18 15:19:27
本帖最后由 chen0201 于 2024-7-18 15:22 编辑
sobereva 发表于 2024-5-6 12:51
今天看到一篇Nature文章用Multiwfn通过本文介绍的方法计算了动力学直径
Xiangjie Zhang, Tao Yan, Huaming ...

请问老师,我按照您的方法计算了分子动力学直径,但是对于苯乙炔这种分子来说,动力学直径是用躺着的苯环(分子动力学直径小~4)还是用立着的苯环(~6.5)来描述比较合适呢?或者有更好定量描述类似苯乙炔分子空间位阻效应的方法吗?
sobereva 发表于 Post on 2024-5-6 12:51:18
今天看到一篇Nature文章用Multiwfn通过本文介绍的方法计算了动力学直径
Xiangjie Zhang, Tao Yan, Huaming Hou, et al., Regioselective hydroformylation of propene catalysed by rhodium-zeolite, Nature (2024) https://www.nature.com/articles/s41586-024-07342-y

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
njust-lbc + 2 牛!

查看全部评分 View all ratings

sobereva 发表于 Post on 2024-5-5 23:41:43
zhangweiwei 发表于 2024-5-5 05:14
卢老师您好!@sobereva 这个方法我已经基本掌握,并验算了CO2。但是我按照您这个帖子,尝试验证了一下CH4的 ...

Multiwfn里旋转后,导出gjf文件,Gaussian重新算个单点(加上nosymm避免Gaussian自动修改朝向)得到新的波函数文件,重新再按照博文里的步骤计算直径就完了
zhangweiwei 发表于 Post on 2024-5-5 05:14:43
本帖最后由 zhangweiwei 于 2024-5-5 22:25 编辑

卢老师您好!@sobereva 这个方法我已经基本掌握,并验算了CO2。但是我按照您这个帖子,尝试验证了一下CH4的计算,发现有较大偏差(算出来4.537,原文4.046)。请问是我步骤不对,还是没有领悟原文/帖子的精髓?在第二步一定需要fch文件吗?因为如果我旋转了,就只能输出pdb/xyz文件,然后用pdb/xyz文件计算0.0015的等势面。
我的步骤如下:
gausssian view新建一个CH4.gjf,如帖子中一样操作,计算出chk转成fch。由于此时C-H键并不是平行于xyz轴,于是按照【Multiwfn中非常实用的几何操作和坐标变换功能介绍http://bbs.keinsci.com/forum.php ... 24674&fromuid=62697】,将C-H键旋转与z轴平行,然后生成pdb/xyz文件。然后将这个pdb/xyz文件按照您帖子里的步骤计算生成vtx.pdb得到如下信息:
Minimum X is   -2.20800000 Angstrom, at atom      40(C )
Minimum Y is   -2.20800000 Angstrom, at atom    1710(C )
Minimum Z is   -1.99600000 Angstrom, at atom    2491(C )
Maximum X is    2.44900000 Angstrom, at atom    4995(C )
Maximum Y is    2.44900000 Angstrom, at atom    2302(C )
Maximum Z is    2.54100000 Angstrom, at atom    2493(C )
可以看到最短方向Z轴,长度为4.537,与原文4.046相差很大。
文件附在后面。如果是您,您会怎么做?恳请解答,谢谢!

vtx.pdb

391.4 KB, 下载次数 Times of downloads: 0

sobereva 发表于 Post on 2024-5-4 17:01:42
zhangweiwei 发表于 2024-5-4 00:22
卢老师好!@sobereva 您说:“如JPCA文中表2所示,用0.0015 a.u.作为等值面数值计算后,原理上最好再乘上拟 ...

除以
zhangweiwei 发表于 Post on 2024-5-4 00:22:44
卢老师好!@sobereva 您说:“如JPCA文中表2所示,用0.0015 a.u.作为等值面数值计算后,原理上最好再乘上拟合出来的系数1.025。”,文中得到的拟合曲线斜率是1.025,如图。y是QM结果,x是breck实际结果。应该是除以1.025才对吧?

6{@91}3GW_FIAQ7MPX{T)2X.png (30.51 KB, 下载次数 Times of downloads: 60)

6{@91}3GW_FIAQ7MPX{T)2X.png
sobereva 发表于 Post on 2023-7-7 10:14:12
ydba 发表于 2023-7-7 05:07
sob老师,可以通过http://sobereva.com/426这篇帖子中提到的“将分子摆正”后的new.pdb文件,再读取表面顶 ...

你可以用Multiwfn摆正后的结构计算波函数文件,按本文做法计算表面顶点位置
ydba 发表于 Post on 2023-7-7 05:07:29
sob老师,可以通过http://sobereva.com/426这篇帖子中提到的“将分子摆正”后的new.pdb文件,再读取表面顶点的最大、最小坐标吗得到动力学直径吗?
Shine剪水 发表于 Post on 2019-8-19 20:44:26
http://sobereva.com/190好像和这篇博文的方法3,4很类似
youyno 发表于 Post on 2019-8-19 20:42:25
大神棒棒哒!

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

GMT+8, 2024-11-23 10:12 , Processed in 0.205843 second(s), 33 queries , Gzip On.

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