gmx sorient分析溶质周围的溶剂取向

gmx sorient [-f [<.xtc/.trr/...>]] [-s [<.tpr/.tpb/...>]] [-n [<.ndx>]]

[-o [<.xvg>]] [-no [<.xvg>]] [-ro [<.xvg>]] [-co [<.xvg>]]

[-rc [<.xvg>]] [-nice ] [-b ] [-e ] [-dt ]

[-[no]w] [-xvg ] [-[no]com] [-[no]v23] [-rmin ]

[-rmax ] [-cbin ] [-rbin ] [-[no]pbc]

gmx sorient用于分析溶质分子周围的溶剂分子的取向. 它可以计算从一个或多个参考位置到每个溶剂分子第一个原子的向量(向量 A⃗ A→)与另外两个向量之间的角度:

θ1θ1: 向量 A⃗ A→ 与从溶剂分子第一个原子到第二和第三个原子中点的向量之间的夹角

θ2θ2: 向量 A⃗ A→ 与由三个原子定义的溶剂分子平面的法线之间的夹角, 或者, 当使用-v23选项时, 向量 A⃗ A→ 与从原子2到原子3的向量之间的夹角.

参考位置可以是一组原子或是一组原子的质心. 溶剂原子组中的每个溶剂分子只能包含3个原子. 对每一帧, -o和-no选项只会考虑处于-rmin和-rmax之间的溶剂分子.

-o: rmin≤r≤rmaxrmin≤r≤rmax 范围内 cosθ1cos⁡θ1 的分布

-no: rmin≤r≤rmaxrmin≤r≤rmax 范围内 cosθ2cos⁡θ2 的分布

-ro: ⟨cosθ1⟩⟨cos⁡θ1⟩ 和 ⟨3cos2θ2−1⟩⟨3cos2⁡θ2−1⟩ 与距离 rr 的函数关系

-co: 对距离 rr 范围内所有溶剂分子的 cosθ1cos⁡θ1 和 3cos2θ2−13cos2⁡θ2−1 进行加和, 得到它们与 rr 的函数关系

-rc: 溶剂分子的分布与 rr 的函数关系

补充说明

此程序特别适用于计算溶质分子周围水分子的角度分布.

设溶质为单原子离子或分子质心Ref, 溶剂为水分子原子1为O, 原子2和3为H, 

则 θ1θ1 对应Ref至O的向量 A⃗ =R⃗ Ref−R⃗ OA→=R→Ref−R→O 与O至两个H连线中点的向量 R⃗ OH2+R⃗ OH32R→OH2+R→OH32 之间的夹角, 

后一向量的方向与水分子偶极矩的方向相同. 因此, θ1θ1 可视为溶质分子周围水分子偶极矩的取向. θ2θ2 对应 A⃗ A→ 与水分子平面法线的夹角. 

当使用-v23选项时, 则为 A⃗ A→ 与两个H连线 R⃗ H3−R⃗ H2R→H3−R→H2 之间的夹角.