计算化学公社

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

[辅助/分析程序] 计算两个平面夹角的小程序

[复制链接 Copy URL]

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

跳转到指定楼层 Go to specific reply
楼主
2021-Mar-3补充:后来在Multiwfn里加入了对动力学轨迹计算指定的两个平面间夹角随帧号变化的功能,见《计算分子动力学轨迹中两个环平面间的距离和夹角》(http://sobereva.com/590)。


今天思想家公社QQ群有人问怎么计算两个平面的夹角。如果本身是二面角好办直接量就行了,但如果不属于这种情况就稍微麻烦,虽然有一些可视化程序如diamond可以量但终究为这点事安装一个也麻烦。这里提供个自写的不足挂齿的程序twoplane2angle来计算夹角。先输入定义第一个平面的三个点的坐标,再输入定义第二个平面的三个点的坐标,然后就会输出夹角。

twoplane2angle.rar (330.44 KB, 下载次数 Times of downloads: 359)

评分 Rate

参与人数
Participants 6
eV +22 收起 理由
Reason
哎呦呦 + 2 好物!
greatzdk + 4 好物!
kulaomega + 5 好物!
小范范1989 + 5 赞!
978142355 + 4 好物!
卡开发发 + 2

查看全部评分 View all ratings

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

326

帖子

0

威望

5312

eV
积分
5638

Level 6 (一方通行)

2#
发表于 Post on 2015-8-21 15:22:50 | 只看该作者 Only view this author
能否提供源码学习下

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

3#
 楼主 Author| 发表于 Post on 2015-8-21 15:25:55 | 只看该作者 Only view this author
mol 发表于 2015-8-21 15:22
能否提供源码学习下

源码非常简单,见这里
http://sobereva.com/302
北京科音自然科学研究中心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

1295

帖子

0

威望

6289

eV
积分
7584

Level 6 (一方通行)

4#
发表于 Post on 2015-8-21 17:56:17 | 只看该作者 Only view this author
sob老师好,我想咨询一下如下的问题:
比如说,我想看图中两个苯环的夹角,这个不能用高斯view的二面角来表示是吧?因为高斯view中二面角需要四个原子来确定另个面,显然这两个面不是我想要的。
所以借助老师刚刚写的这个程序,我输入了3个苯环1的坐标和另外3个苯环2的坐标,这样得到的结果就是两个面的夹角了吧?如下图。
不晓得这样对不对?谢谢老师。

{HN07GU)$~%AH3AVZGW]PWQ.png (70.95 KB, 下载次数 Times of downloads: 120)

1

1

27PF~]0Y[BT%W]@5OQS@73L.png (10.53 KB, 下载次数 Times of downloads: 104)

27PF~]0Y[BT%W]@5OQS@73L.png
https://www.x-mol.com/groups/fan_jianzhong

strive for greatness

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

5#
 楼主 Author| 发表于 Post on 2015-8-21 21:01:37 | 只看该作者 Only view this author
小范范1989 发表于 2015-8-21 17:56
sob老师好,我想咨询一下如下的问题:
比如说,我想看图中两个苯环的夹角,这个不能用高斯view的二面角来 ...

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

1295

帖子

0

威望

6289

eV
积分
7584

Level 6 (一方通行)

6#
发表于 Post on 2015-8-22 08:17:50 | 只看该作者 Only view this author
https://www.x-mol.com/groups/fan_jianzhong

strive for greatness

144

帖子

0

威望

2820

eV
积分
2964

Level 5 (御坂)

7#
发表于 Post on 2015-8-28 12:49:06 | 只看该作者 Only view this author
还可以考虑用VMD中的vector操作来进行。也是相当的容易实现。

10

帖子

0

威望

259

eV
积分
269

Level 3 能力者

8#
发表于 Post on 2020-11-5 22:03:23 | 只看该作者 Only view this author
yjmaxpayne 发表于 2015-8-28 12:49
还可以考虑用VMD中的vector操作来进行。也是相当的容易实现。

你好 请问VECTOR在哪个模块里呀 没有找到

144

帖子

0

威望

2820

eV
积分
2964

Level 5 (御坂)

9#
发表于 Post on 2020-11-6 13:35:11 | 只看该作者 Only view this author
紫薯糯米小丸子 发表于 2020-11-5 22:03
你好 请问VECTOR在哪个模块里呀 没有找到

TCL Console, 例子如下:https://www.ks.uiuc.edu/Research/vmd/vmd-1.7.1/ug/node163.html

本版积分规则 Credits rule

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

GMT+8, 2026-2-21 06:34 , Processed in 0.234250 second(s), 24 queries , Gzip On.

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