计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 量子化学 (Quantum Chemistry) › 请教C60的旋转如何进行势能面扫描？

请教C60的旋转如何进行势能面扫描？

查看数: 959 | 评论数: 13 | 收藏 Add to favorites 2

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

lr8616

发布时间: 2024-1-14 09:07

正文摘要:

请问想模仿卢老师这篇文章中C18旋转的势能面扫描，试了半天没搞清楚该怎么操作，难道是需要手动摆结构然后计算吗？

回复 Reply

ionexchangeC 发表于 Post on 2024-1-19 18:01:18

想要旋转至少得有个旋转轴吧，如果以之前这个位置的18碳环所在平面为“地面”，敢问你是要让C60绕“俯仰”、“滚转”还是“偏航”轴旋转？

lr8616 发表于 Post on 2024-1-17 20:01:33

chrinide 发表于 2024-1-17 15:50
我觉得你这样做势能面扫描的意义有点费解。
首先，扫描势能面只能完全不考虑两个复合分子的形变（扫描过程 ...

谢谢~我会认真考虑下的

lr8616 发表于 Post on 2024-1-17 20:01:19

丁越发表于 2024-1-17 10:10
用VMD中的平移+旋转变换就可以实现

首先用Multiwfn把体系旋转使得环处于XY平面，在300 -> 7 -> 11选项可 ...

好详细！太感谢啦！！

chrinide 发表于 Post on 2024-1-17 15:50:10

本帖最后由 chrinide 于 2024-1-17 16:51 编辑

我觉得你这样做势能面扫描的意义有点费解。
首先，扫描势能面只能完全不考虑两个复合分子的形变（扫描过程中）。
可行方案一：分别优化两个复合分子，按楼上的的方法，重合几何中心（这一步具有任意性），固定C60（假定为Ih对称），按固定角度固定平面扫描旋转环状分子。这个扫面得到的“势能面”意义有点费解，因为C60（假定Ih对称性）不是球面，只有C60是完美球面，这个“势能面”扫描才意义明确。
对偶思考一下就清楚明白：固定环状分子，重合几何中心，让C60来旋转，C60可以绕几何中心任意旋转，让C60在球坐标下旋转至一个非Ih对性的位置，再固定C60, 按照前述方法旋转环状分子来扫描，是不是又是一个不同的“势能面”？
方案二，优化符合物结构得到一个可能的全局最优结构，再来扫描，这个时候可能有个问题，如果几何优化是完全不约束的优化，可能最终的复合物结构几何中心不重叠，无法良好的定义扫描。
方案三，固定一个分子不变进行优化，同样会出现方案一方案二的问题。

所以，综上，这样做的“势能面”意义非常费解，可能会被审稿人质疑。

丁越 发表于 Post on 2024-1-17 10:10:08

本帖最后由丁越于 2024-1-17 12:41 编辑

用VMD中的平移+旋转变换就可以实现

首先用Multiwfn把体系旋转使得环处于XY平面，在300 -> 7 -> 11选项可实现这个目的。
然后把得到的结构保存为xyz格式，载入VMD进行如下操作：
1. 先预定义几个要频繁使用变量，第一个是选择体系所有原子，其次定义C60的质心

set sel [atomselect top all]
set com [measure center [atomselect top {same fragment as index 10}] weight mass] #假设C60上的一个C原子序号为10，根据fragment选择语法就可以选中整个C60

复制代码

2. 把整个体系平移，使得C60的质心挪到坐标原点