计算化学公社

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

[新手求助] 如何正确选取原子测定分子结构中两个环的二面角

[复制链接 Copy URL]

588

帖子

0

威望

3310

eV
积分
3898

Level 5 (御坂)


如图分子,我想求红圈中两个环的二面角,应该在GaussView中选取哪四个原子或者还有什么其它方法?谢谢指教!





222

帖子

5

威望

2601

eV
积分
2923

Level 5 (御坂)

2#
发表于 Post on 2021-1-19 10:17:41 | 只看该作者 Only view this author
本帖最后由 liuyuje714 于 2021-1-19 10:36 编辑

我们一般说的二面角是指四个原子构成的那种,你这种就是个两平面夹角。我不知道gaussview能不能量,但自己可以写个小程序。
  • 如果这两个环你当做纯刚性处理(所有原子共面),则只需要任意选着三点,代表平面然后求的平面法向量,最后两个环的法向量即可求出夹角。
  • 如果你认为两个环并不是纯刚性,那就需要先对每一个环的所有原子做最小二乘拟合得到平面方程(得法向量),然后同上。
一般简化处理就用1,此法的参考可见sob博文http://sobereva.com/302
复杂的拟合用2,此法参考https://jerkwin.github.io/2018/08/29/Pi-Pi%E5%A0%86%E7%A7%AF%E8%B7%9D%E7%A6%BB%E5%92%8C%E5%A0%86%E7%A7%AF%E8%A7%92%E5%BA%A6%E7%9A%84%E8%AE%A1%E7%AE%97/

1万

帖子

0

威望

9900

eV
积分
22154

Level 6 (一方通行)

3#
发表于 Post on 2021-1-19 11:06:03 | 只看该作者 Only view this author
这个问题不能用gaussview直接测量的方法来解决,因为gaussview要求选定的4个原子中,中间的那两个原子必须在两个面相交的棱上。建议查一下,有两个面,已知每个面上面的各3个点的坐标,如何求二面角。或者自己用立体几何推一下也行
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
本团队长期招收研究生,有意者可私信联系

1245

帖子

1

威望

3860

eV
积分
5125

Level 6 (一方通行)

4#
发表于 Post on 2021-1-19 11:09:45 | 只看该作者 Only view this author
得先找到两个平面相交的线,提取出来原子坐标用算一下应该就有了

135

帖子

1

威望

4135

eV
积分
4290

Level 6 (一方通行)

5#
发表于 Post on 2021-1-19 11:28:06 | 只看该作者 Only view this author
既然无法用 GaussView 直接量角,那直接计算也不复杂。

计算方法是,1. 在这两个苯环上各选三个点,计算出两个法向量,2. 求这两个法向量的夹角。(因为法向量的方向跟选取的三个点有关,所以二面角也可能是法向量夹角的补角。)
十八介姑娘一蕾花呀,白白介牙齿、红红介嘴唇,得人惜

334

帖子

1

威望

6163

eV
积分
6517

Level 6 (一方通行)

6#
发表于 Post on 2021-1-19 11:51:55 | 只看该作者 Only view this author
GaussView没法测量不具有2个共同原子的两个环的夹角。通用的方法就如楼上所说用到空间解析几何的方法
1. 选取平面上三点得到平面方程 2. 根据平面方程可以求出两个平面二面角  


自己如果觉得算着麻烦,这里补充一下两个在线的计算器。可以借助下面两个在线的计算器得到夹角。

https://keisan.casio.com/exec/system/1223596129  输入三点坐标得到平面方程的 在线计算器  


https://onlinemschool.com/math/a ... rdinate/plane_angl/   根据两个平面方程求二面角的在线工具

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
ABetaCarw + 5 谢谢分享

查看全部评分 View all ratings

588

帖子

0

威望

3310

eV
积分
3898

Level 5 (御坂)

7#
 楼主 Author| 发表于 Post on 2021-1-19 12:07:01 | 只看该作者 Only view this author
王二葛 发表于 2021-1-19 11:28
既然无法用 GaussView 直接量角,那直接计算也不复杂。

计算方法是,1. 在这两个苯环上各选三个点,计算 ...

感谢回复

588

帖子

0

威望

3310

eV
积分
3898

Level 5 (御坂)

8#
 楼主 Author| 发表于 Post on 2021-1-19 12:08:43 | 只看该作者 Only view this author
paramecium86 发表于 2021-1-19 11:51
GaussView没法测量不具有2个共同原子的两个环的夹角。通用的方法就如楼上所说用到空间解析几何的方法
1.  ...

非常感谢

1043

帖子

0

威望

4188

eV
积分
5231

Level 6 (一方通行)

9#
发表于 Post on 2021-1-19 12:22:20 来自手机 | 只看该作者 Only view this author
合理做法是:分别确定两个环的最小二乘平面的法向矢量(用Multiwfn好像可以做),然后求这两个法矢的点积的反余弦就得交角。

330

帖子

0

威望

3018

eV
积分
3348

Level 5 (御坂)

计算化学路人甲

10#
发表于 Post on 2021-1-19 14:41:57 | 只看该作者 Only view this author
你可以导出成xyz或者sdf格式,然后在Mercury中测量一下,Mercury的有些常用功能是免费的。

541

帖子

3

威望

6624

eV
积分
7225

Level 6 (一方通行)

11#
发表于 Post on 2021-1-19 14:49:09 | 只看该作者 Only view this author
Novice 发表于 2021-1-19 14:41
你可以导出成xyz或者sdf格式,然后在Mercury中测量一下,Mercury的有些常用功能是免费的。

可以具体说说 什么模块能够做这种测量吗?

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
Novice + 1 见下帖

查看全部评分 View all ratings

恍惚月余,深谙人与人之间的差距。以后还应努力学习,才能与强者比肩。

330

帖子

0

威望

3018

eV
积分
3348

Level 5 (御坂)

计算化学路人甲

12#
发表于 Post on 2021-1-19 15:04:13 | 只看该作者 Only view this author
本帖最后由 Novice 于 2021-1-19 15:25 编辑

选中位于一个平面上的至少三个原子,利用calculate菜单下面的plane功能产生一个平面,同样的步骤产生另一个平面。然后将工具栏中的“picking mode"改为”measure Angles",点击产生的两个平面,即在图上显示出了两个面的夹角。下图是我随便产生的二苯甲烷的二面角:



评分 Rate

参与人数
Participants 3
eV +10 收起 理由
Reason
共和国面壁者 + 3 好东西
biogon + 5
王二葛 + 2 赞!

查看全部评分 View all ratings

1万

帖子

0

威望

9900

eV
积分
22154

Level 6 (一方通行)

13#
发表于 Post on 2021-1-19 19:23:02 | 只看该作者 Only view this author
Novice 发表于 2021-1-19 15:04
选中位于一个平面上的至少三个原子,利用calculate菜单下面的plane功能产生一个平面,同样的步骤产生另一个 ...

原来真的可以这样做,我之前只在mercury里量过一个原子到一个平面的距离
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
本团队长期招收研究生,有意者可私信联系

330

帖子

0

威望

3018

eV
积分
3348

Level 5 (御坂)

计算化学路人甲

14#
发表于 Post on 2021-1-19 20:26:40 | 只看该作者 Only view this author
嗯嗯,图形操作很方便,越做功能越丰富,好像新版必须要有序列号?我也没有下载最新版,不太清楚。

6万

帖子

99

威望

6万

eV
积分
125208

管理员

公社社长

15#
发表于 Post on 2021-1-20 04:58:38 | 只看该作者 Only view this author
我写过现成的
计算两个平面夹角的小程序
http://sobereva.com/302http://bbs.keinsci.com/thread-1623-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

本版积分规则 Credits rule

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

GMT+8, 2026-2-27 18:42 , Processed in 0.858184 second(s), 24 queries , Gzip On.

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