计算化学公社

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

[Multiwfn使用咨询] 求助计算每一帧的自由体积以及原子分布的描述

[复制链接 Copy URL]

7

帖子

0

威望

41

eV
积分
48

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
老师您好,如果我想算轨迹文件每一帧的自由体积,Multiwfn我该怎么设置,或者有其他办法?另外,自由体积周围原子分布的描述(比如周围原子种类数量分布,配位数等)有没有什么好的表征方法或者思路,应该怎么操作,感谢帮助

6万

帖子

99

威望

5万

eV
积分
120061

管理员

公社社长

2#
发表于 Post on 2022-12-1 07:00:41 | 只看该作者 Only view this author
先把轨迹文件拆分成单帧结构文件(比如对于GROMACS,可以用gmx trjconv结合-sep)。怎么写脚本批量计算,看此文
详谈Multiwfn的命令行方式运行和批量运行的方法
http://sobereva.com/612http://bbs.keinsci.com/thread-24929-1-1.html

“自由体积周围原子分布”指代不明。至少拿具体体系、实际问题说事
北京科音自然科学研究中心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

7

帖子

0

威望

41

eV
积分
48

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2022-12-26 21:57:55 | 只看该作者 Only view this author
sobereva 发表于 2022-12-1 07:00
先把轨迹文件拆分成单帧结构文件(比如对于GROMACS,可以用gmx trjconv结合-sep)。怎么写脚本批量计算,看 ...

老师,自由体积的分布,比如模型空洞的大小及不同大小孔洞的数量分布,这个该怎么分析?"自由体积周围的原子分布"这个我想算的是玻璃模型,比如钠钙硅玻璃,想知道空隙周围的原子分布状态,类似于径向分布函数分析?我该怎么做?谢谢老师

6万

帖子

99

威望

5万

eV
积分
120061

管理员

公社社长

4#
发表于 Post on 2022-12-26 23:11:13 | 只看该作者 Only view this author
xinyoulingxi 发表于 2022-12-26 21:57
老师,自由体积的分布,比如模型空洞的大小及不同大小孔洞的数量分布,这个该怎么分析?"自由体积周围的原 ...

下文提到的功能可以给出孔洞数目、各个孔洞的尺寸,允许用gro、带晶胞信息的pdb文件等作为输入文件
使用Multiwfn可视化分子孔洞并计算孔洞体积
http://sobereva.com/408http://bbs.keinsci.com/thread-9331-1-1.html

“自由体积周围原子分布”问题得自己写分析脚本,比如VMD里的tcl脚本。会写脚本就没有分析不了的
北京科音自然科学研究中心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

7

帖子

0

威望

41

eV
积分
48

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2022-12-27 21:32:54 | 只看该作者 Only view this author
sobereva 发表于 2022-12-26 23:11
下文提到的功能可以给出孔洞数目、各个孔洞的尺寸,允许用gro、带晶胞信息的pdb文件等作为输入文件
使用 ...

好的,谢谢老师,我试试

7

帖子

0

威望

41

eV
积分
48

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2023-1-2 17:56:15 | 只看该作者 Only view this author
sobereva 发表于 2022-12-26 23:11
下文提到的功能可以给出孔洞数目、各个孔洞的尺寸,允许用gro、带晶胞信息的pdb文件等作为输入文件
使用 ...

老师还有个问题,同一个模型用Multiwfn中300的功能算出来的自由体积和使用200功能算出来的域孔洞体积的和是不是相等的?因为我是先用300的功能做了自由体积可视化的图,然后再想分析不同尺寸自由体积的累积体积百分数,这两个结果不一样的话结果可能会不可靠。。。

6万

帖子

99

威望

5万

eV
积分
120061

管理员

公社社长

7#
发表于 Post on 2023-1-3 03:52:59 | 只看该作者 Only view this author
xinyoulingxi 发表于 2023-1-2 17:56
老师还有个问题,同一个模型用Multiwfn中300的功能算出来的自由体积和使用200功能算出来的域孔洞体积的和 ...

“300的功能” 、“200功能” 指代不明。说清楚什么主功能的什么子功能
北京科音自然科学研究中心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

7

帖子

0

威望

41

eV
积分
48

Level 2 能力者

8#
 楼主 Author| 发表于 Post on 2023-1-3 18:07:55 | 只看该作者 Only view this author
sobereva 发表于 2023-1-3 03:52
“300的功能” 、“200功能” 指代不明。说清楚什么主功能的什么子功能

《使用Multiwfn图形化展示分子动力学模拟体系中的孔洞、自由区域》(http://sobereva.com/539)和《使用Multiwfn计算晶体结构中自由区域的体积、图形化展现自由区域》(http://sobereva.com/617)就是老师这两篇文章里用到的功能

6万

帖子

99

威望

5万

eV
积分
120061

管理员

公社社长

9#
发表于 Post on 2023-1-6 08:59:31 | 只看该作者 Only view this author
xinyoulingxi 发表于 2023-1-3 18:07
《使用Multiwfn图形化展示分子动力学模拟体系中的孔洞、自由区域》(http://sobereva.com/539)和《使用M ...

这两篇文章都只用了主功能300的子功能1
北京科音自然科学研究中心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

7

帖子

0

威望

41

eV
积分
48

Level 2 能力者

10#
 楼主 Author| 发表于 Post on 2023-1-10 23:30:09 | 只看该作者 Only view this author
sobereva 发表于 2023-1-6 08:59
这两篇文章都只用了主功能300的子功能1

说错了老师,我想问的是这个文章《使用Multiwfn可视化分子孔洞并计算孔洞体积》(http://sobereva.com/408)中使用200中功能计算的孔洞体积的和和这个文章《使用Multiwfn图形化展示分子动力学模拟体系中的孔洞、自由区域》(http://sobereva.com/539)中使用300功能计算的孔洞自由体积是不是相等?这两个有什么联系,之前问的时候写错了,不好意思

6万

帖子

99

威望

5万

eV
积分
120061

管理员

公社社长

11#
发表于 Post on 2023-1-12 06:56:38 | 只看该作者 Only view this author
xinyoulingxi 发表于 2023-1-10 23:30
说错了老师,我想问的是这个文章《使用Multiwfn可视化分子孔洞并计算孔洞体积》(http://sobereva.com/40 ...

不相等。
主功能300那个计算自由体积的时候是根据原子范德华半径定义的体系的范德华表面,在其以外就是自由区域。
而《使用Multiwfn可视化分子孔洞并计算孔洞体积》里是基于原子孤立状态密度叠加而成的promolecular密度等值面定义的范德华表面。
因此两种方法的范德华表面定义不同,自由体积数值肯定有差异。
北京科音自然科学研究中心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

7

帖子

0

威望

41

eV
积分
48

Level 2 能力者

12#
 楼主 Author| 发表于 Post on 2023-1-12 11:36:36 | 只看该作者 Only view this author
sobereva 发表于 2023-1-12 06:56
不相等。
主功能300那个计算自由体积的时候是根据原子范德华半径定义的体系的范德华表面,在其以外就是 ...

那老师如何用范德华半径定义的体系的范德华表面按200的功能分析空洞尺寸大小和不同空洞的数量分布,就是相当于对300计算的自由体积结果做一个对于空洞尺寸和不同尺寸空洞数量的精细分析

6万

帖子

99

威望

5万

eV
积分
120061

管理员

公社社长

13#
发表于 Post on 2023-1-13 06:21:33 | 只看该作者 Only view this author
xinyoulingxi 发表于 2023-1-12 11:36
那老师如何用范德华半径定义的体系的范德华表面按200的功能分析空洞尺寸大小和不同空洞的数量分布,就是 ...

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

威望

543

eV
积分
625

Level 4 (黑子)

14#
发表于 Post on 2023-3-27 15:44:50 | 只看该作者 Only view this author
兄弟 解决了吗?

本版积分规则 Credits rule

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

GMT+8, 2025-8-12 21:13 , Processed in 2.068352 second(s), 20 queries , Gzip On.

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