layout: post title: GROMACS分析教程:水分子序参数的计算 categories:
水分子的序参数有很多种, 不同序参数的计算方法也不尽相同. GROMACS自带的分析程序gmx h2order
可以计算其中的两种序参数, 对于这种可以直接使用已有命令进行计算的序参数, 就无需赘言了.
有些序参数的计算比较复杂, 没有直接可以使用的命令, 这时就需要我们自己进行计算. 如果你会写代码计算就容易些, 否则的话, 或者到网上寻找已有的代码/脚本, 或者咨询/求助下有经验的人.
如果你确定要自己写代码进行计算, 那我的建议还是和以前一样: 充分利用GROMACS已有的功能和自带的工具, 使用简单的脚本组合其输入输出, 整理所得的数据, 获得自己需要的结果. 这种做法比你从头写一个分析程序简单不少, 因为你不需要自己读取轨迹, 处理PBC, 查找原子, 也不需要计算距离, 角度等, 只需要想清楚计算流程, 将需要执行的步骤写在脚本中就可以了. 当然, 有时还是需要进行一些简单的中间处理, 但大多是准备脚本所需要的输入和输出.
本教程就以一个比较复杂的序参数来示例下具体的操作过程, 做法与以前氢键数目的计算类似, 只不过稍微复杂一些.
在水合物的模拟中, 常会使用角度序参数AOP(angular order parameter), 三体序参数F3(也称四面体序参数, tetrahedral order parameter)和四体序参数 F4(four-body order parameter)来表征水分子所处的状态. 其中AOP和F3是相同的序参数, 只是叫法不同而已, 所以不需要进行区分. 一般来说, 在文章中如果只讨论F3, 那大多会称其为AOP或F, 如果同时还讨论F4, 那多半会称其为F3, 而不称AOP或F.
AOP/F3是一个由三体(三个水分子)构型组成的参数, 表征了中心氧原子与周围3.5Å范围内的其他氧原子形成的四面体与正四面体的偏离程度. 其计算公式如下:
$$\alg \text{AOP}/F_3 &= {1\over n_i(ni-1)/2}\sum{j=1}^{ni-1} \sum{k=j+1}^{ni} \left( \left| \cos\theta{jik} \right|\cos\theta_{jik} + \cos^2(109.47^\circ) \right)^2 \ &= {1\over n_i(ni-1)/2}\sum{j=1}^{ni-1} \sum{k=j+1}^{ni} \left( \left| \cos\theta{jik} \right|\cos\theta_{jik} + {1\over9} \right)^2 \ealg$$
所涉及的角度为中心水分子的氧原子i与其周围3.5Å范围内任意两个其他水分子的氧原子j和k之间的夹角, 其中109.47度为正四面体中心-顶点连线之间夹角, 也是3维空间中4个向量间最小夹角的最大值 $$\arccos \left({-1 \over 3}\right)=2\arctan(\sqrt 2)≈109.4712°$$
有文献中认为109.47°是水分子的键角, 也有文献中对此角度使用水分子的键角, 这些都是不正确的.
本质上说, AOP/F3就是判断中心氧原子周围的氧原子是否处于四面体的角度位置. 但涉及具体的公式, 在文献上有几种不同的版本. 主要区别在于:
不同版本所得的AOP/F3对一些特殊体系的值是不同的. 比如, 文献上水合物AOP/F3常见的值就有两个, 0.1, 或0.01. 对本文最后所给的几篇参考文献进行统计, 所得结果如下
<table id='tab-0'><caption></caption> <tr> <th rowspan="1" colspan="1" style="text-align:center;">方法</th> <th rowspan="1" colspan="1" style="text-align:center;">冰</th> <th rowspan="1" colspan="1" style="text-align:center;">固态水</th> <th rowspan="1" colspan="1" style="text-align:center;">水合物</th> <th rowspan="1" colspan="1" style="text-align:center;">液态水</th> </tr> <tr> <td rowspan="1" colspan="1" style="text-align:center;">平方/平均/键角</td> <td rowspan="1" colspan="1" style="text-align:center;">0</td> <td rowspan="1" colspan="1" style="text-align:center;"></td> <td rowspan="1" colspan="1" style="text-align:center;"></td> <td rowspan="1" colspan="1" style="text-align:center;">0.1</td> </tr> <tr> <td rowspan="1" colspan="1" style="text-align:center;">平方/平均/键角</td> <td rowspan="1" colspan="1" style="text-align:center;"></td> <td rowspan="1" colspan="1" style="text-align:center;">0</td> <td rowspan="1" colspan="1" style="text-align:center;"></td> <td rowspan="1" colspan="1" style="text-align:center;">0.1</td> </tr> <tr> <td rowspan="1" colspan="1" style="text-align:center;">平方/平均</td> <td rowspan="1" colspan="1" style="text-align:center;"></td> <td rowspan="1" colspan="1" style="text-align:center;"></td> <td rowspan="1" colspan="1" style="text-align:center;"></td> <td rowspan="1" colspan="1" style="text-align:center;"></td> </tr> <tr> <td rowspan="1" colspan="1" style="text-align:center;">平方/平均</td> <td rowspan="1" colspan="1" style="text-align:center;"></td> <td rowspan="1" colspan="1" style="text-align:center;">0.01</td> <td rowspan="1" colspan="1" style="text-align:center;"></td> <td rowspan="1" colspan="1" style="text-align:center;">0.1</td> </tr> <tr> <td rowspan="1" colspan="1" style="text-align:center;">平方/平均</td> <td rowspan="1" colspan="1" style="text-align:center;">0</td> <td rowspan="1" colspan="1" style="text-align:center;"></td> <td rowspan="1" colspan="1" style="text-align:center;">0.1</td> <td rowspan="1" colspan="1" style="text-align:center;">0.8</td> </tr> <tr> <td rowspan="1" colspan="1" style="text-align:center;">平方/未平均</td> <td rowspan="1" colspan="1" style="text-align:center;">0</td> <td rowspan="1" colspan="1" style="text-align:center;"></td> <td rowspan="1" colspan="1" style="text-align:center;">0.1</td> <td rowspan="1" colspan="1" style="text-align:center;">0.8, 界面 0.4</td> </tr> <tr> <td rowspan="1" colspan="1" style="text-align:center;">平方/未平均</td> <td rowspan="1" colspan="1" style="text-align:center;">0</td> <td rowspan="1" colspan="1" style="text-align:center;">0.02</td> <td rowspan="1" colspan="1" style="text-align:center;"></td> <td rowspan="1" colspan="1" style="text-align:center;"></td> </tr> <tr> <td rowspan="1" colspan="1" style="text-align:center;">平方/未知</td> <td rowspan="1" colspan="1" style="text-align:center;"></td> <td rowspan="1" colspan="1" style="text-align:center;">0.01</td> <td rowspan="1" colspan="1" style="text-align:center;"></td> <td rowspan="1" colspan="1" style="text-align:center;">0.1</td> </tr> <tr> <td rowspan="1" colspan="1" style="text-align:center;">未平方/平均/10倍</td> <td rowspan="1" colspan="1" style="text-align:center;"></td> <td rowspan="1" colspan="1" style="text-align:center;"></td> <td rowspan="1" colspan="1" style="text-align:center;">0.1</td> <td rowspan="1" colspan="1" style="text-align:center;">0.9</td> </tr> <tr> <td rowspan="1" colspan="1" style="text-align:center;">未知</td> <td rowspan="1" colspan="1" style="text-align:center;"></td> <td rowspan="1" colspan="1" style="text-align:center;">0.01</td> <td rowspan="1" colspan="1" style="text-align:center;"></td> <td rowspan="1" colspan="1" style="text-align:center;">0.1</td> </tr> </table>可见, 大多数关于水合物模拟的文献采用的计算公式都与上文的定义相同. 在这种定义下, 固态水和水合物的AOP/F3为0.01, 而液态水的AOP/F3为0.1. 液态水AOP/F3为0.8的两篇文献来自同一个研究组, 但给出的公式却不一致, 可以推测他们没有注意到这个问题.
我们来比较下平方和未平方版本公式对角度的依赖性.
平方版本所得的AOP/F3在100度附近AOP/F3对角度变化不敏感, 这样力场, 温度等对其影响较小.
如果我们假定每个角度都是独立的, 那么根据空间中两随机向量间夹角的概率密度分布为 $p(\q)={1 \over 2}\sin\q$, 由此可以得到F3的平均值为
$$\alg \left< F_3 \right> &=\int_0^\p F_3(\q) p(\q) d\q=\int_0^\p {1\over 2} \left( |\cos\q| \cos\q+{1\over9} \right)^n \sin\q d\q \ &= \begin{cases} {1 \over 9} \approx 0.11111 & n=1 \ {86 \over 405} \approx 0.21235 & n=2 \end{cases} \ealg$$
由于所得结果与液态水的AOP/F3值存在很大差距, 我们可以认为在液态水中, 对应的夹角并不是随机分布的.
F4统计的是水体系中两个相邻水分子最外侧氢原子以及水分子中氧原子组成二面角:
$$F4 = {1\over n}\Sum{i=1}^n \cos 3\phi_i$$
一些体系的F4值: 冰 -0.4(Ice-Ih -0.5; Ice-II -0.3; Ice-IV -0.4); 液态水 -0.04; 水合物 0.7(sI型 0.89; sII型 0.96).
对正四面体的冰, 其F4精确值为
$$\alg 4F_4 &= \cos(3\times 120°) +\cos(3\times 180°) \ &+\cos(3\times -60°)+\cos(3\times 60°)=-2 \ealg$$
当然, 实际的冰有多种结构, 所以其F4值并不精确地等于-0.5. 这就是为什么文献上经常会使用-0.4的原因.
AOP/F3, F4的值对于固态水(水合物与冰)和液态水有较大差别, 都可用于判定某区域内水分子的相态情况. 但相对而言, F4的表征能力更好一些, 因为AOP/F3不包括氢原子的信息, 不容易区分各种固态结构的水. 因此很多文章直接使用F4, 而没有给出AOP/F3.
理解了定义, 那我们来看一下, 如何计算AOP/F3
最简单最直接, 也是最容易想到的方法就是, 对每一帧轨迹中每个需要考虑的水分子, 查找距离其氧原子3.5 Å之内的其他氧原子, 对于找到的这些氧原子, 计算其中任意两个与中心氧原子所成角度, 然后根据定义计算相应的AOP/F3值.
具体来说, 如果中心氧原子的编号为 $I_0$, 距离其3.5Å之内的其他氧原子的编号为 $I_1$, $I_2$, $I_3$, $I_4$, ..., $I_n$, 那么可能的角度总数为C(n,2)=n(n-1)/2, 这些角度的编号三联对为 $I_1 I_0 I_2$, $I_1 I_0 I_3$, $I_1 I_0 I_4$, ..., $I_1 I_0 I_n$, $I_2 I_0 I_3$, $I_2 I_0 I_4$, ..., $I_2 I_0 In$, ..., $I{n-1} I_0 I_n$.
所以, 我们需要做的工作包括下面几点
gmx select
可以做到gmx gangle
可以做到所以, 任务分解后, GROMACS自带的工具可以帮我们完成计算最复杂的部分, 我们要做的只是简单的处理输入输出.
假定我们的要计算的中心氧原子编号为$cntAtom
, gro结构文件名称为$conf
, 水中氧原子的名称为OW
,
对于任务1, 我们需要的执行的命令如下
<div class="highlight"><pre style="line-height:125%"><span style="color:#A2F">gmx</span> select <span style="color:#666">-f</span> $conf.gro <span style="color:#666">-on</span> $conf.ndx <span style="color:#666">-s</span> <span style="color:#666">-quiet</span> << EOS <span style="color:#A2F">r</span> = distance from atomnr $cntAtom <span style="color:#A2F">atomname</span> OW and 0. < r and r <.35 <span style="color:#A2F">EOS</span></pre></div>执行成功后, 得到的氧原子编号在输出的$conf.ndx
中.
对于任务2, 我们需要读取$conf.ndx
中的原子编号, 然后生成角度三联对. 这不是困难的任务, 只要一个双重循环就好了, 可以使用你熟悉的任何语言来实现. 我们将角度编号三联对放到另一个索引文件$pid.ndx
中.
对于任务3, 指定前面得到的索引文件中的角度三联对, 直接调用命令即可
<div class="highlight"><pre style="line-height:125%"><span style="color:#A2F">gmx</span> gangle <span style="color:#666">-f</span> $conf.gro <span style="color:#666">-n</span> $pid.ndx <span style="color:#666">-group1</span> 0 <span style="color:#666">-s</span> <span style="color:#666">-oall</span> <span style="color:#666">-quiet</span></pre></div>这样在默认的输出文件angles.xvg
中就得到了所有的角度值.
对于任务4, 我们需要自己处理下得到的角度. 就简单的方法就是将这个文件的内容追加到一个另外的数据文件中, 然后用自己熟悉的工具来计算AOP/F3. 当然, 也可以一步到位, 直接将AOP/F3输出到数据文件中.
F4的计算也是类似, 但可能更麻烦. 因为在定义中角度 $\f$ 对应距离最远的两个氢原子, 这意味着我们不仅要知道需要的氧原子, 还要确定出对应的4个氢原子对中哪一对的距离最远, 从而计算相应的二面角. 所以我们需要增加一步计算并比较距离的任务, 从而才能得到所需要的二面角编号四联对. 距离的计算最好还是借助gmx distance
, 避免自己写代码.
具体来说, 主要包括以下几点
第一步: 输出与中心氧原子, 选中氧原子对应的氢原子的距离编号对.
第二步: 根据上一步得到的氢原子编号对, 通过命令计算出4对氢原子的距离.
<div class="highlight"><pre style="line-height:125%"><span style="color:#A2F">gmx</span> distance <span style="color:#666">-f</span> $conf.gro <span style="color:#666">-n</span> $pid.ndx <span style="color:#666">-select</span> 0 <span style="color:#666">-s</span> <span style="color:#666">-oall</span> <span style="color:#666">-quiet</span></pre></div>第三步: 根据距离的最大值确定要计算的二面角, 并输出要计算二面角的编号四联对.
第四步: 直接调用命令计算二面角.
<div class="highlight"><pre style="line-height:125%"><span style="color:#A2F">gmx</span> gangle <span style="color:#666">-f</span> $conf.gro <span style="color:#666">-n</span> $pid.ndx <span style="color:#666">-g1</span> dihedral <span style="color:#666">-group1</span> 0 <span style="color:#666">-s</span> <span style="color:#666">-oall</span> <span style="color:#666">-quiet</span></pre></div>最后, 用自己熟悉的工具处理下得到的二面角, 将F4输出到数据文件中.
我们已经知道如何处理一帧构型, 那就将其变为一个脚本AOP.bsh
, 其命令行参数只有一个, 那就是每帧的编号, 这样就可以借助gmx trjconv -exec
功能处理所有的轨迹帧了.
总结一下流程就是:
gmx trjconv -exec
执行此脚本, 处理所有的帧. 当然, 你也可以先输出所有的帧, 然后脚本处理, 虽然麻烦一些.使用的时候只需要修改中心氧原子的编号cntAtom
即可. 具体来说
其中的3
表示使用的水模型为三位点水模型, 如SPC, SPC/E, TIP3P等. 如果你使用四位点水模型, 如TIP4P, 应相应地改为4
. 此外, 这种简单的循环方式只适用于连续的水分子, 且其中的原子排列顺序为OW HW1 HW2
. 如果你的所有水分子在整个体系的编号系统中被其他分子分成了好几个部分, 那你或者不使用这中简单的循环方式, 或者分别计算每部分水分子.
说明:
gmx gangle
是支持选区, 但我们没法直接提供要计算的角度编号三联对, 因为gmx select
动态选区语法不够强大. 如果gmx select
能够支持组合所选原子的编号, 那计算F3需要的角度只要简单的使用gmx gangle
就可以了.gmx gangle
计算时对角度编号三联对的数目有限制, 具体多大没有仔细考察. 但如果使用的截断距离扩大10倍, 达到35Å, 那在8GB内存的Windows下是无法处理的.gmx trjconv
. 实际上你也可以自己创建一个脚本, 自己执行.IceIh_Hayward
冰IceIh_Hayward
冰结构来自以前的博文常见冰的晶体结构及其cif文件. 下图为其晶体结构:
我们任意选取其中的两个中心氧原子进行计算, 例如: 10号氧原子(坐标3.38 2.60 5.06
), 70号氧原子(坐标5.63 6.50 5.06
). 使用Materials Studio(也可采用VMD, Chemcraft, GaussView等)软件测量周围原子与中心原子所成的角度和二面角.
修改AOP.bsh
脚本中的中心原子编号得到不同中心氧原子的计算结果, 将其与软件测量的结果进行对比:
$\q$, $\f$, F3和F4的分布如下
采用与示例1相同的方法和步骤, 我们计算了I型甲烷水合物结构中水分子的序参数. I型甲烷水合物中的水分子为四位点的TIP4P模型. 下图为I型甲烷水合物的晶体结构.
<figure><script>var Mol2=new ChemDoodle.TransformCanvas3D('Mol-2',502,376.5);Mol2.specs.shapes_color = '#fff';Mol2.specs.backgroundColor = 'black';Mol2.specs.projectionPerspective_3D = false;Mol2.specs.set3DRepresentation('Ball and Stick');Mol2.specs.compass_display = true;//Mol2.specs.atoms_resolution_3D = 15; //Mol2.specs.bonds_resolution_3D = 15; Mol2.handle = null;Mol2.timeout = 15;Mol2.specs.crystals_unitCellLineWidth = 1.5;Mol2.startAnimation = ChemDoodle._AnimatorCanvas.prototype.startAnimation;Mol2.stopAnimation = ChemDoodle._AnimatorCanvas.prototype.stopAnimation;Mol2.isRunning = ChemDoodle._AnimatorCanvas.prototype.isRunning;Mol2.dblclick = ChemDoodle.RotatorCanvas.prototype.dblclick;Mol2.nextFrame = function(delta){var matrix = [];ChemDoodle.lib.mat4.identity(matrix);var change = delta*Math.PI/15000;ChemDoodle.lib.mat4.rotate(matrix,change,[1,0,0]);ChemDoodle.lib.mat4.rotate(matrix,change,[0,1,0]);ChemDoodle.lib.mat4.rotate(matrix,change,[0,0,1]);ChemDoodle.lib.mat4.multiply(this.rotationMatrix, matrix)};var Fmol='data_Gallium_Rubidium_Oxygen_Manganese_Argon_Carbon_Silicon\n1472\n_symmetry_space_group_name_H-M \'P 1\'\n_cell_length_a 23.47\n_cell_length_b 23.47\n_cell_length_c 23.47\n_cell_angle_alpha 90\n_cell_angle_beta 90\n_cell_angle_gamma 90\nloop_\n_atom_site_label\n_atom_site_type_symbol\n_atom_site_fract_x\n_atom_site_fract_y\n_atom_site_fract_z\n1SOL O 0.545 0.700 0.604\n1SOL H 0.577 0.678 0.616\n1SOL H 0.513 0.678 0.616\n2SOL O 0.545 0.898 0.605\n2SOL H 0.545 0.861 0.622\n2SOL H 0.545 0.890 0.565\n3SOL O 0.545 0.700 0.991\n3SOL H 0.577 0.678 0.980\n3SOL H 0.513 0.678 0.980\n4SOL O 0.545 0.895 0.992\n4SOL H 0.577 0.917 0.981\n4SOL H 0.545 0.895 1.033\n5SOL Cl 0.604 0.545 0.701\n5SOL H 0.616 0.577 0.680\n5SOL H 0.628 0.545 0.735\n6SOL O 0.605 0.545 0.898\n6SOL H 0.622 0.545 0.861\n6SOL H 0.565 0.545 0.890\n7SOL O 0.991 0.545 0.701\n7SOL H 0.968 0.544 0.735\n7SOL H 0.979 0.512 0.680\n8SOL O 0.991 0.545 0.895\n8SOL H 0.980 0.577 0.917\n8SOL H 0.980 0.513 0.917\n9SOL O 0.698 0.605 0.545\n9SOL H 0.735 0.622 0.545\n9SOL H 0.705 0.565 0.545\n10SOL O 0.898 0.605 0.545\n10SOL H 0.861 0.622 0.545\n10SOL H 0.890 0.565 0.545\n11SOL O 0.700 0.992 0.545\n11SOL H 0.700 1.033 0.544\n11SOL H 0.679 0.981 0.512\n12SOL O 0.895 0.991 0.545\n12SOL H 0.916 0.979 0.577\n12SOL H 0.861 0.968 0.545\n13SOL O 0.954 0.798 0.738\n13SOL H 0.933 0.830 0.726\n13SOL H 0.987 0.798 0.715\n14SOL O 0.642 0.798 0.739\n14SOL H 0.642 0.798 0.780\n14SOL H 0.664 0.830 0.728\n15SOL O 0.951 0.798 0.859\n15SOL H 0.958 0.798 0.818\n15SOL H 0.988 0.798 0.876\n16SOL O 0.642 0.798 0.857\n16SOL H 0.663 0.830 0.869\n16SOL H 0.608 0.798 0.881\n17SOL O 0.798 0.856 0.642\n17SOL H 0.798 0.816 0.642\n17SOL H 0.830 0.867 0.664\n18SOL O 0.798 0.856 0.954\n18SOL H 0.798 0.816 0.954\n18SOL H 0.830 0.867 0.932\n19SOL O 0.798 0.738 0.643\n19SOL H 0.766 0.726 0.664\n19SOL H 0.830 0.726 0.664\n20SOL O 0.798 0.738 0.953\n20SOL H 0.830 0.726 0.931\n20SOL H 0.766 0.726 0.931\n21SOL O 0.856 0.954 0.798\n21SOL H 0.816 0.954 0.798\n21SOL H 0.867 0.932 0.830\n22SOL O 0.857 0.642 0.797\n22SOL H 0.869 0.663 0.765\n22SOL H 0.881 0.608 0.797\n23SOL O 0.738 0.954 0.797\n23SOL H 0.726 0.932 0.765\n23SOL H 0.715 0.987 0.797\n24SOL O 0.739 0.642 0.798\n24SOL H 0.780 0.642 0.798\n24SOL H 0.728 0.664 0.830\n25SOL O 0.637 0.637 0.637\n25SOL H 0.660 0.627 0.605\n25SOL H 0.662 0.661 0.660\n26SOL O 0.957 0.958 0.637\n26SOL H 0.991 0.937 0.627\n26SOL H 0.971 0.989 0.659\n27SOL O 0.958 0.637 0.958\n27SOL H 0.934 0.661 0.935\n27SOL H 0.935 0.627 0.990\n28SOL O 0.637 0.958 0.958\n28SOL H 0.662 0.935 0.935\n28SOL H 0.660 0.969 0.990\n29SOL O 0.890 0.890 0.705\n29SOL H 0.879 0.914 0.737\n29SOL H 0.914 0.915 0.682\n30SOL O 0.705 0.705 0.704\n30SOL H 0.683 0.736 0.718\n30SOL H 0.715 0.684 0.738\n31SOL O 0.890 0.705 0.891\n31SOL H 0.912 0.736 0.877\n31SOL H 0.880 0.684 0.857\n32SOL O 0.704 0.890 0.890\n32SOL H 0.738 0.880 0.911\n32SOL H 0.718 0.912 0.859\n33SOL O 0.958 0.958 0.957\n33SOL H 0.989 0.936 0.971\n33SOL H 0.937 0.968 0.991\n34SOL O 0.637 0.638 0.958\n34SOL H 0.627 0.604 0.937\n34SOL H 0.659 0.624 0.989\n35SOL O 0.638 0.958 0.638\n35SOL H 0.606 0.936 0.624\n35SOL H 0.658 0.968 0.604\n36SOL O 0.958 0.638 0.638\n36SOL H 0.968 0.604 0.658\n36SOL H 0.936 0.624 0.606\n37SOL O 0.705 0.705 0.890\n37SOL H 0.682 0.737 0.879\n37SOL H 0.680 0.682 0.914\n38SOL O 0.890 0.890 0.890\n38SOL H 0.914 0.859 0.879\n38SOL H 0.915 0.913 0.914\n39SOL O 0.705 0.890 0.705\n39SOL H 0.737 0.879 0.682\n39SOL H 0.682 0.914 0.680\n40SOL O 0.890 0.705 0.705\n40SOL H 0.914 0.737 0.716\n40SOL H 0.915 0.682 0.682\n41SOL O 0.545 0.798 0.671\n41SOL H 0.546 0.764 0.648\n41SOL H 0.579 0.796 0.694\n42SOL O 0.545 0.798 0.924\n42SOL H 0.545 0.766 0.949\n42SOL H 0.545 0.830 0.949\n43SOL O 0.671 0.545 0.798\n43SOL H 0.696 0.577 0.798\n43SOL H 0.696 0.513 0.798\n44SOL O 0.924 0.545 0.798\n44SOL H 0.947 0.543 0.832\n44SOL H 0.902 0.511 0.799\n45SOL O 0.798 0.671 0.545\n45SOL H 0.798 0.696 0.577\n45SOL H 0.798 0.696 0.513\n46SOL O 0.798 0.924 0.545\n46SOL H 0.764 0.947 0.546\n46SOL H 0.796 0.902 0.579\n47SOL O 0.038 0.701 0.604\n47SOL H 0.006 0.680 0.616\n47SOL H 0.038 0.735 0.628\n48SOL O 0.039 0.895 0.603\n48SOL H 0.072 0.917 0.614\n48SOL H 0.040 0.895 0.563\n49SOL O 0.039 0.701 0.991\n49SOL H 0.006 0.680 0.979\n49SOL H 0.038 0.735 0.968\n50SOL O 0.039 0.895 0.992\n50SOL H 0.072 0.917 0.981\n50SOL H 0.040 0.895 1.033\n51SOL O 0.098 0.545 0.700\n51SOL H 0.057 0.545 0.700\n51SOL H 0.109 0.577 0.679\n52SOL O 0.098 0.545 0.895\n52SOL H 0.057 0.545 0.895\n52SOL H 0.109 0.577 0.917\n53SOL O 0.486 0.545 0.701\n53SOL H 0.527 0.544 0.700\n53SOL H 0.475 0.512 0.679\n54SOL O 0.485 0.545 0.895\n54SOL H 0.473 0.577 0.916\n54SOL H 0.462 0.545 0.861\n55SOL O 0.195 0.604 0.544\n55SOL H 0.228 0.628 0.544\n55SOL H 0.173 0.616 0.512\n56SOL O 0.389 0.604 0.545\n56SOL H 0.410 0.616 0.577\n56SOL H 0.355 0.628 0.545\n57SOL O 0.195 0.991 0.545\n57SOL H 0.228 0.968 0.544\n57SOL H 0.174 0.979 0.512\n58SOL O 0.389 0.992 0.545\n58SOL H 0.411 0.981 0.577\n58SOL H 0.389 1.033 0.545\n59SOL O 0.445 0.798 0.737\n59SOL H 0.482 0.798 0.720\n59SOL H 0.452 0.798 0.778\n60SOL O 0.137 0.798 0.738\n60SOL H 0.159 0.766 0.726\n60SOL H 0.159 0.830 0.726\n61SOL O 0.448 0.798 0.857\n61SOL H 0.481 0.798 0.881\n61SOL H 0.427 0.830 0.869\n62SOL O 0.136 0.798 0.856\n62SOL H 0.136 0.797 0.816\n62SOL H 0.158 0.765 0.867\n63SOL O 0.292 0.856 0.642\n63SOL H 0.293 0.816 0.642\n63SOL H 0.325 0.867 0.664\n64SOL O 0.292 0.857 0.953\n64SOL H 0.259 0.869 0.931\n64SOL H 0.324 0.869 0.931\n65SOL O 0.292 0.738 0.643\n65SOL H 0.259 0.726 0.664\n65SOL H 0.324 0.726 0.664\n66SOL O 0.292 0.737 0.951\n66SOL H 0.292 0.778 0.958\n66SOL H 0.292 0.720 0.988\n67SOL O 0.352 0.954 0.797\n67SOL H 0.363 0.932 0.765\n67SOL H 0.375 0.987 0.797\n68SOL O 0.352 0.642 0.798\n68SOL H 0.363 0.663 0.830\n68SOL H 0.375 0.608 0.798\n69SOL O 0.233 0.954 0.798\n69SOL H 0.274 0.954 0.798\n69SOL H 0.222 0.932 0.830\n70SOL O 0.231 0.645 0.798\n70SOL H 0.271 0.637 0.798\n70SOL H 0.214 0.608 0.798\n71SOL O 0.132 0.637 0.638\n71SOL H 0.100 0.659 0.624\n71SOL H 0.153 0.627 0.604\n72SOL O 0.452 0.958 0.637\n72SOL H 0.485 0.937 0.627\n72SOL H 0.465 0.989 0.659\n73SOL O 0.452 0.637 0.958\n73SOL H 0.427 0.661 0.935\n73SOL H 0.429 0.627 0.990\n74SOL O 0.132 0.957 0.958\n74SOL H 0.121 0.991 0.937\n74SOL H 0.153 0.971 0.989\n75SOL O 0.385 0.890 0.705\n75SOL H 0.408 0.859 0.716\n75SOL H 0.409 0.914 0.682\n76SOL O 0.199 0.705 0.705\n76SOL H 0.210 0.682 0.737\n76SOL H 0.176 0.680 0.682\n77SOL O 0.385 0.704 0.890\n77SOL H 0.405 0.738 0.880\n77SOL H 0.353 0.718 0.912\n78SOL O 0.199 0.890 0.890\n78SOL H 0.176 0.859 0.879\n78SOL H 0.175 0.913 0.914\n79SOL O 0.452 0.958 0.958\n79SOL H 0.485 0.937 0.968\n79SOL H 0.466 0.989 0.936\n80SOL O 0.132 0.637 0.958\n80SOL H 0.100 0.660 0.969\n80SOL H 0.154 0.662 0.934\n81SOL O 0.132 0.958 0.637\n81SOL H 0.155 0.969 0.605\n81SOL H 0.156 0.935 0.660\n82SOL O 0.452 0.637 0.637\n82SOL H 0.463 0.606 0.660\n82SOL H 0.429 0.660 0.662\n83SOL O 0.199 0.705 0.890\n83SOL H 0.232 0.716 0.911\n83SOL H 0.212 0.683 0.859\n84SOL O 0.385 0.890 0.890\n84SOL H 0.374 0.914 0.859\n84SOL H 0.408 0.915 0.914\n85SOL O 0.199 0.890 0.705\n85SOL H 0.232 0.880 0.684\n85SOL H 0.213 0.912 0.736\n86SOL O 0.385 0.705 0.704\n86SOL H 0.406 0.736 0.718\n86SOL H 0.374 0.684 0.738\n87SOL O 0.039 0.798 0.671\n87SOL H 0.040 0.832 0.648\n87SOL H 0.073 0.799 0.694\n88SOL O 0.039 0.798 0.924\n88SOL H 0.040 0.831 0.947\n88SOL H 0.073 0.799 0.902\n89SOL O 0.166 0.545 0.798\n89SOL H 0.141 0.545 0.766\n89SOL H 0.141 0.545 0.830\n90SOL O 0.418 0.545 0.798\n90SOL H 0.441 0.543 0.764\n90SOL H 0.395 0.511 0.796\n91SOL O 0.292 0.671 0.545\n91SOL H 0.292 0.696 0.577\n91SOL H 0.292 0.696 0.513\n92SOL O 0.292 0.924 0.545\n92SOL H 0.326 0.947 0.546\n92SOL H 0.293 0.902 0.579\n93SOL O 0.545 0.192 0.605\n93SOL H 0.545 0.229 0.622\n93SOL H 0.545 0.199 0.565\n94SOL O 0.545 0.392 0.605\n94SOL H 0.545 0.355 0.622\n94SOL H 0.545 0.384 0.565\n95SOL O 0.545 0.195 0.991\n95SOL H 0.544 0.228 0.968\n95SOL H 0.512 0.174 0.980\n96SOL O 0.545 0.389 0.991\n96SOL H 0.577 0.410 0.980\n96SOL H 0.545 0.355 0.968\n97SOL O 0.603 0.039 0.701\n97SOL H 0.614 0.006 0.679\n97SOL H 0.563 0.038 0.700\n98SOL O 0.604 0.039 0.896\n98SOL H 0.616 0.006 0.917\n98SOL H 0.616 0.071 0.917\n99SOL O 0.990 0.039 0.698\n99SOL H 0.973 0.039 0.735\n99SOL H 1.030 0.039 0.705\n100SOL O 0.991 0.039 0.894\n100SOL H 0.979 0.006 0.916\n100SOL H 0.968 0.038 0.861\n101SOL O 0.700 0.098 0.545\n101SOL H 0.678 0.110 0.577\n101SOL H 0.678 0.110 0.513\n102SOL O 0.895 0.098 0.545\n102SOL H 0.895 0.057 0.545\n102SOL H 0.917 0.109 0.577\n103SOL O 0.700 0.485 0.545\n103SOL H 0.678 0.474 0.577\n103SOL H 0.678 0.473 0.513\n104SOL O 0.895 0.485 0.545\n104SOL H 0.861 0.462 0.544\n104SOL H 0.916 0.473 0.512\n105SOL O 0.953 0.292 0.738\n105SOL H 0.931 0.324 0.726\n105SOL H 0.931 0.259 0.726\n106SOL O 0.643 0.292 0.738\n106SOL H 0.664 0.259 0.726\n106SOL H 0.664 0.324 0.726\n107SOL O 0.951 0.292 0.859\n107SOL H 0.958 0.292 0.818\n107SOL H 0.988 0.292 0.876\n108SOL O 0.642 0.292 0.856\n108SOL H 0.642 0.291 0.816\n108SOL H 0.664 0.259 0.867\n109SOL O 0.798 0.351 0.642\n109SOL H 0.798 0.310 0.642\n109SOL H 0.830 0.362 0.664\n110SOL O 0.798 0.351 0.954\n110SOL H 0.797 0.310 0.954\n110SOL H 0.765 0.362 0.932\n111SOL O 0.798 0.232 0.642\n111SOL H 0.765 0.220 0.663\n111SOL H 0.797 0.209 0.608\n112SOL O 0.798 0.232 0.954\n112SOL H 0.830 0.220 0.933\n112SOL H 0.798 0.209 0.987\n113SOL O 0.856 0.448 0.798\n113SOL H 0.816 0.447 0.798\n113SOL H 0.867 0.426 0.830\n114SOL O 0.856 0.136 0.798\n114SOL H 0.816 0.136 0.797\n114SOL H 0.867 0.158 0.765\n115SOL O 0.738 0.447 0.798\n115SOL H 0.726 0.426 0.830\n115SOL H 0.726 0.425 0.766\n116SOL O 0.738 0.136 0.798\n116SOL H 0.726 0.157 0.830\n116SOL H 0.715 0.102 0.798\n117SOL O 0.638 0.132 0.637\n117SOL H 0.604 0.153 0.627\n117SOL H 0.624 0.100 0.659\n118SOL O 0.958 0.452 0.638\n118SOL H 0.937 0.463 0.604\n118SOL H 0.989 0.430 0.624\n119SOL O 0.958 0.132 0.958\n119SOL H 0.969 0.100 0.935\n119SOL H 0.934 0.154 0.934\n120SOL O 0.637 0.452 0.958\n120SOL H 0.627 0.485 0.937\n120SOL H 0.659 0.465 0.989\n121SOL O 0.890 0.385 0.705\n121SOL H 0.879 0.408 0.737\n121SOL H 0.914 0.409 0.682\n122SOL O 0.705 0.199 0.705\n122SOL H 0.716 0.176 0.737\n122SOL H 0.681 0.175 0.682\n123SOL O 0.890 0.199 0.891\n123SOL H 0.912 0.231 0.877\n123SOL H 0.879 0.178 0.857\n124SOL O 0.705 0.385 0.890\n124SOL H 0.682 0.353 0.879\n124SOL H 0.680 0.407 0.914\n125SOL O 0.958 0.452 0.958\n125SOL H 0.934 0.429 0.935\n125SOL H 0.935 0.463 0.990\n126SOL O 0.637 0.132 0.957\n126SOL H 0.606 0.153 0.971\n126SOL H 0.658 0.121 0.991\n127SOL O 0.638 0.452 0.637\n127SOL H 0.624 0.484 0.659\n127SOL H 0.604 0.431 0.627\n128SOL O 0.957 0.132 0.637\n128SOL H 0.971 0.100 0.659\n128SOL H 0.991 0.153 0.627\n129SOL O 0.705 0.199 0.890\n129SOL H 0.737 0.210 0.914\n129SOL H 0.682 0.176 0.915\n130SOL O 0.890 0.385 0.890\n130SOL H 0.911 0.352 0.879\n130SOL H 0.859 0.371 0.912\n131SOL O 0.705 0.385 0.705\n131SOL H 0.737 0.374 0.682\n131SOL H 0.682 0.408 0.680\n132SOL O 0.890 0.199 0.705\n132SOL H 0.859 0.210 0.682\n132SOL H 0.914 0.176 0.680\n133SOL O 0.545 0.292 0.671\n133SOL H 0.577 0.292 0.696\n133SOL H 0.513 0.292 0.696\n134SOL O 0.545 0.292 0.925\n134SOL H 0.577 0.292 0.900\n134SOL H 0.513 0.292 0.899\n135SOL O 0.671 0.039 0.798\n135SOL H 0.646 0.039 0.766\n135SOL H 0.647 0.039 0.830\n136SOL O 0.925 0.039 0.798\n136SOL H 0.899 0.007 0.798\n136SOL H 0.900 0.071 0.798\n137SOL O 0.798 0.166 0.545\n137SOL H 0.766 0.141 0.545\n137SOL H 0.830 0.141 0.545\n138SOL O 0.798 0.418 0.545\n138SOL H 0.764 0.441 0.546\n138SOL H 0.796 0.395 0.579\n139SOL O 0.039 0.195 0.604\n139SOL H 0.072 0.174 0.616\n139SOL H 0.040 0.228 0.628\n140SOL O 0.039 0.389 0.604\n140SOL H 0.072 0.410 0.616\n140SOL H 0.040 0.355 0.628\n141SOL O 0.039 0.195 0.992\n141SOL H 0.006 0.173 0.981\n141SOL H 0.038 0.195 1.033\n142SOL O 0.039 0.389 0.992\n142SOL H 0.006 0.411 0.981\n142SOL H 0.038 0.389 1.033\n143SOL O 0.098 0.038 0.701\n143SOL H 0.110 0.006 0.680\n143SOL H 0.122 0.038 0.735\n144SOL O 0.099 0.039 0.898\n144SOL H 0.059 0.039 0.890\n144SOL H 0.116 0.039 0.861\n145SOL O 0.485 0.039 0.701\n145SOL H 0.473 0.072 0.680\n145SOL H 0.462 0.040 0.735\n146SOL O 0.486 0.039 0.895\n146SOL H 0.527 0.040 0.895\n146SOL H 0.475 0.072 0.917\n147SOL O 0.195 0.098 0.545\n147SOL H 0.195 0.057 0.544\n147SOL H 0.173 0.109 0.512\n148SOL O 0.389 0.098 0.544\n148SOL H 0.355 0.122 0.544\n148SOL H 0.410 0.110 0.512\n149SOL O 0.195 0.486 0.545\n149SOL H 0.195 0.527 0.545\n149SOL H 0.173 0.475 0.577\n150SOL O 0.389 0.486 0.545\n150SOL H 0.389 0.527 0.544\n150SOL H 0.411 0.475 0.512\n151SOL O 0.447 0.292 0.738\n151SOL H 0.426 0.324 0.726\n151SOL H 0.426 0.259 0.726\n152SOL O 0.136 0.292 0.739\n152SOL H 0.136 0.291 0.780\n152SOL H 0.158 0.259 0.728\n153SOL O 0.448 0.292 0.856\n153SOL H 0.447 0.293 0.816\n153SOL H 0.426 0.325 0.867\n154SOL O 0.135 0.292 0.857\n154SOL H 0.157 0.325 0.869\n154SOL H 0.102 0.293 0.881\n155SOL O 0.292 0.352 0.642\n155SOL H 0.259 0.363 0.663\n155SOL H 0.291 0.375 0.608\n156SOL O 0.291 0.352 0.954\n156SOL H 0.259 0.363 0.932\n156SOL H 0.291 0.375 0.987\n157SOL O 0.292 0.233 0.642\n157SOL H 0.291 0.274 0.642\n157SOL H 0.259 0.222 0.664\n158SOL O 0.292 0.233 0.954\n158SOL H 0.291 0.274 0.954\n158SOL H 0.259 0.222 0.932\n159SOL O 0.352 0.447 0.798\n159SOL H 0.363 0.425 0.766\n159SOL H 0.363 0.426 0.830\n160SOL O 0.352 0.137 0.798\n160SOL H 0.363 0.159 0.830\n160SOL H 0.363 0.159 0.766\n161SOL O 0.231 0.445 0.798\n161SOL H 0.214 0.482 0.798\n161SOL H 0.271 0.452 0.798\n162SOL O 0.233 0.136 0.798\n162SOL H 0.274 0.136 0.798\n162SOL H 0.222 0.158 0.830\n163SOL O 0.132 0.132 0.637\n163SOL H 0.121 0.098 0.658\n163SOL H 0.153 0.118 0.606\n164SOL O 0.452 0.452 0.638\n164SOL H 0.484 0.430 0.624\n164SOL H 0.432 0.463 0.604\n165SOL O 0.452 0.132 0.958\n165SOL H 0.427 0.155 0.935\n165SOL H 0.429 0.121 0.990\n166SOL O 0.132 0.452 0.958\n166SOL H 0.118 0.484 0.936\n166SOL H 0.098 0.431 0.969\n167SOL O 0.385 0.385 0.705\n167SOL H 0.353 0.374 0.682\n167SOL H 0.407 0.408 0.680\n168SOL O 0.199 0.199 0.705\n168SOL H 0.210 0.176 0.737\n168SOL H 0.176 0.175 0.682\n169SOL O 0.385 0.199 0.890\n169SOL H 0.406 0.232 0.879\n169SOL H 0.353 0.213 0.912\n170SOL O 0.199 0.385 0.890\n170SOL H 0.210 0.408 0.859\n170SOL H 0.176 0.409 0.914\n171SOL O 0.452 0.452 0.958\n171SOL H 0.466 0.484 0.936\n171SOL H 0.485 0.431 0.969\n172SOL O 0.132 0.132 0.958\n172SOL H 0.098 0.153 0.969\n172SOL H 0.118 0.100 0.936\n173SOL O 0.132 0.452 0.637\n173SOL H 0.121 0.484 0.660\n173SOL H 0.155 0.429 0.662\n174SOL O 0.452 0.132 0.638\n174SOL H 0.484 0.153 0.624\n174SOL H 0.432 0.121 0.604\n175SOL O 0.199 0.199 0.890\n175SOL H 0.176 0.231 0.879\n175SOL H 0.175 0.176 0.914\n176SOL O 0.385 0.385 0.890\n176SOL H 0.353 0.374 0.914\n176SOL H 0.407 0.408 0.915\n177SOL O 0.199 0.385 0.704\n177SOL H 0.177 0.353 0.718\n177SOL H 0.210 0.405 0.738\n178SOL O 0.385 0.199 0.705\n178SOL H 0.353 0.210 0.682\n178SOL H 0.407 0.176 0.680\n179SOL O 0.039 0.292 0.671\n179SOL H 0.007 0.292 0.696\n179SOL H 0.071 0.292 0.696\n180SOL O 0.039 0.292 0.924\n180SOL H 0.039 0.259 0.949\n180SOL H 0.039 0.324 0.949\n181SOL O 0.165 0.039 0.798\n181SOL H 0.190 0.007 0.798\n181SOL H 0.190 0.071 0.798\n182SOL O 0.418 0.039 0.798\n182SOL H 0.441 0.040 0.832\n182SOL H 0.396 0.073 0.799\n183SOL O 0.292 0.165 0.545\n183SOL H 0.258 0.143 0.546\n183SOL H 0.291 0.188 0.579\n184SOL O 0.292 0.418 0.545\n184SOL H 0.259 0.443 0.545\n184SOL H 0.324 0.443 0.545\n185SOL O 0.545 0.701 0.098\n185SOL H 0.544 0.700 0.057\n185SOL H 0.512 0.679 0.109\n186SOL O 0.545 0.895 0.098\n186SOL H 0.577 0.916 0.110\n186SOL H 0.545 0.861 0.122\n187SOL O 0.545 0.700 0.486\n187SOL H 0.545 0.700 0.527\n187SOL H 0.577 0.679 0.475\n188SOL O 0.545 0.895 0.485\n188SOL H 0.544 0.861 0.462\n188SOL H 0.512 0.916 0.473\n189SOL O 0.604 0.545 0.194\n189SOL H 0.616 0.577 0.173\n189SOL H 0.616 0.513 0.173\n190SOL O 0.605 0.545 0.392\n190SOL H 0.622 0.545 0.355\n190SOL H 0.565 0.545 0.384\n191SOL O 0.990 0.545 0.192\n191SOL H 0.973 0.545 0.229\n191SOL H 1.030 0.545 0.199\n192SOL O 0.991 0.545 0.389\n192SOL H 0.980 0.577 0.411\n192SOL H 0.980 0.513 0.411\n193SOL O 0.701 0.604 0.039\n193SOL H 0.680 0.616 0.072\n193SOL H 0.735 0.628 0.040\n194SOL O 0.894 0.604 0.039\n194SOL H 0.916 0.616 0.072\n194SOL H 0.861 0.628 0.040\n195SOL O 0.701 0.991 0.039\n195SOL H 0.680 0.979 0.072\n195SOL H 0.735 0.968 0.040\n196SOL O 0.895 0.992 0.039\n196SOL H 0.917 0.981 0.072\n196SOL H 0.895 1.033 0.040\n197SOL O 0.951 0.798 0.231\n197SOL H 0.958 0.798 0.271\n197SOL H 0.988 0.798 0.214\n198SOL O 0.643 0.798 0.232\n198SOL H 0.664 0.766 0.220\n198SOL H 0.664 0.830 0.220\n199SOL O 0.953 0.798 0.352\n199SOL H 0.931 0.766 0.363\n199SOL H 0.931 0.830 0.363\n200SOL O 0.642 0.798 0.351\n200SOL H 0.642 0.797 0.310\n200SOL H 0.664 0.765 0.362\n201SOL O 0.798 0.857 0.136\n201SOL H 0.830 0.869 0.157\n201SOL H 0.798 0.881 0.102\n202SOL O 0.798 0.857 0.448\n202SOL H 0.798 0.881 0.481\n202SOL H 0.830 0.869 0.427\n203SOL O 0.798 0.739 0.136\n203SOL H 0.798 0.780 0.136\n203SOL H 0.830 0.728 0.158\n204SOL O 0.798 0.739 0.448\n204SOL H 0.798 0.780 0.447\n204SOL H 0.830 0.728 0.426\n205SOL O 0.859 0.951 0.292\n205SOL H 0.818 0.958 0.292\n205SOL H 0.876 0.988 0.292\n206SOL O 0.857 0.641 0.291\n206SOL H 0.869 0.663 0.259\n206SOL H 0.881 0.608 0.291\n207SOL O 0.738 0.954 0.292\n207SOL H 0.726 0.932 0.325\n207SOL H 0.715 0.987 0.293\n208SOL O 0.739 0.642 0.292\n208SOL H 0.780 0.642 0.293\n208SOL H 0.728 0.664 0.325\n209SOL O 0.637 0.637 0.132\n209SOL H 0.605 0.660 0.121\n209SOL H 0.660 0.662 0.155\n210SOL O 0.958 0.958 0.132\n210SOL H 0.934 0.935 0.156\n210SOL H 0.969 0.990 0.155\n211SOL O 0.958 0.637 0.452\n211SOL H 0.937 0.627 0.485\n211SOL H 0.989 0.659 0.465\n212SOL O 0.637 0.958 0.452\n212SOL H 0.606 0.935 0.463\n212SOL H 0.660 0.934 0.429\n213SOL O 0.890 0.890 0.199\n213SOL H 0.912 0.859 0.212\n213SOL H 0.879 0.911 0.232\n214SOL O 0.704 0.705 0.199\n214SOL H 0.738 0.715 0.179\n214SOL H 0.718 0.683 0.231\n215SOL O 0.890 0.705 0.385\n215SOL H 0.879 0.682 0.353\n215SOL H 0.914 0.680 0.407\n216SOL O 0.705 0.891 0.385\n216SOL H 0.684 0.857 0.374\n216SOL H 0.736 0.877 0.406\n217SOL O 0.958 0.957 0.452\n217SOL H 0.936 0.971 0.484\n217SOL H 0.968 0.991 0.432\n218SOL O 0.637 0.638 0.452\n218SOL H 0.659 0.624 0.484\n218SOL H 0.627 0.604 0.431\n219SOL O 0.637 0.958 0.132\n219SOL H 0.661 0.935 0.156\n219SOL H 0.627 0.990 0.155\n220SOL O 0.957 0.638 0.132\n220SOL H 0.971 0.606 0.153\n220SOL H 0.991 0.658 0.121\n221SOL O 0.705 0.705 0.385\n221SOL H 0.737 0.716 0.408\n221SOL H 0.682 0.682 0.409\n222SOL O 0.890 0.890 0.385\n222SOL H 0.879 0.914 0.353\n222SOL H 0.914 0.915 0.407\n223SOL O 0.704 0.890 0.199\n223SOL H 0.738 0.880 0.179\n223SOL H 0.718 0.912 0.231\n224SOL O 0.890 0.705 0.199\n224SOL H 0.914 0.737 0.210\n224SOL H 0.915 0.682 0.176\n225SOL O 0.545 0.798 0.166\n225SOL H 0.546 0.764 0.143\n225SOL H 0.579 0.796 0.188\n226SOL O 0.545 0.798 0.418\n226SOL H 0.546 0.764 0.441\n226SOL H 0.579 0.796 0.395\n227SOL O 0.671 0.545 0.292\n227SOL H 0.648 0.546 0.258\n227SOL H 0.694 0.578 0.291\n228SOL O 0.925 0.545 0.292\n228SOL H 0.947 0.543 0.326\n228SOL H 0.902 0.511 0.293\n229SOL O 0.798 0.671 0.039\n229SOL H 0.798 0.696 0.007\n229SOL H 0.798 0.696 0.071\n230SOL O 0.798 0.924 0.039\n230SOL H 0.799 0.902 0.005\n230SOL H 0.832 0.947 0.037\n231SOL O 0.039 0.701 0.098\n231SOL H 0.040 0.700 0.057\n231SOL H 0.072 0.679 0.109\n232SOL O 0.039 0.895 0.098\n232SOL H 0.006 0.916 0.110\n232SOL H 0.038 0.861 0.122\n233SOL O 0.039 0.698 0.484\n233SOL H 0.039 0.705 0.524\n233SOL H 0.039 0.735 0.467\n234SOL O 0.039 0.895 0.485\n234SOL H 0.006 0.916 0.473\n234SOL H 0.038 0.861 0.462\n235SOL O 0.098 0.545 0.195\n235SOL H 0.110 0.577 0.174\n235SOL H 0.122 0.545 0.228\n236SOL O 0.098 0.545 0.389\n236SOL H 0.057 0.544 0.389\n236SOL H 0.109 0.512 0.411\n237SOL O 0.486 0.545 0.195\n237SOL H 0.527 0.544 0.195\n237SOL H 0.475 0.512 0.173\n238SOL O 0.485 0.545 0.389\n238SOL H 0.473 0.577 0.410\n238SOL H 0.462 0.545 0.355\n239SOL O 0.195 0.603 0.039\n239SOL H 0.173 0.614 0.006\n239SOL H 0.195 0.563 0.038\n240SOL O 0.389 0.603 0.039\n240SOL H 0.411 0.614 0.072\n240SOL H 0.389 0.563 0.040\n241SOL O 0.195 0.992 0.039\n241SOL H 0.173 0.981 0.072\n241SOL H 0.195 1.033 0.040\n242SOL O 0.389 0.991 0.039\n242SOL H 0.410 0.979 0.006\n242SOL H 0.355 0.968 0.038\n243SOL O 0.445 0.798 0.231\n243SOL H 0.482 0.798 0.214\n243SOL H 0.452 0.798 0.271\n244SOL O 0.137 0.798 0.232\n244SOL H 0.159 0.766 0.220\n244SOL H 0.159 0.830 0.220\n245SOL O 0.448 0.798 0.352\n245SOL H 0.481 0.798 0.375\n245SOL H 0.427 0.830 0.363\n246SOL O 0.136 0.798 0.351\n246SOL H 0.136 0.797 0.310\n246SOL H 0.158 0.765 0.362\n247SOL O 0.292 0.857 0.135\n247SOL H 0.324 0.869 0.157\n247SOL H 0.293 0.881 0.102\n248SOL O 0.292 0.857 0.448\n248SOL H 0.293 0.881 0.481\n248SOL H 0.325 0.869 0.427\n249SOL O 0.292 0.737 0.139\n249SOL H 0.292 0.778 0.132\n249SOL H 0.292 0.720 0.102\n250SOL O 0.292 0.739 0.448\n250SOL H 0.293 0.780 0.447\n250SOL H 0.325 0.728 0.426\n251SOL O 0.351 0.954 0.292\n251SOL H 0.310 0.954 0.291\n251SOL H 0.362 0.932 0.259\n252SOL O 0.352 0.643 0.292\n252SOL H 0.363 0.664 0.324\n252SOL H 0.363 0.664 0.259\n253SOL O 0.232 0.953 0.292\n253SOL H 0.220 0.931 0.324\n253SOL H 0.220 0.931 0.259\n254SOL O 0.233 0.642 0.292\n254SOL H 0.274 0.642 0.293\n254SOL H 0.222 0.664 0.324\n255SOL O 0.132 0.637 0.132\n255SOL H 0.155 0.627 0.100\n255SOL H 0.156 0.661 0.154\n256SOL O 0.452 0.958 0.132\n256SOL H 0.484 0.936 0.118\n256SOL H 0.431 0.969 0.098\n257SOL O 0.452 0.637 0.452\n257SOL H 0.431 0.627 0.485\n257SOL H 0.484 0.659 0.465\n258SOL O 0.132 0.958 0.452\n258SOL H 0.100 0.935 0.463\n258SOL H 0.154 0.934 0.429\n259SOL O 0.385 0.890 0.199\n259SOL H 0.408 0.859 0.210\n259SOL H 0.409 0.914 0.176\n260SOL O 0.199 0.705 0.199\n260SOL H 0.232 0.716 0.178\n260SOL H 0.213 0.683 0.231\n261SOL O 0.385 0.705 0.385\n261SOL H 0.408 0.737 0.374\n261SOL H 0.409 0.682 0.408\n262SOL O 0.199 0.891 0.385\n262SOL H 0.178 0.857 0.374\n262SOL H 0.231 0.877 0.406\n263SOL O 0.452 0.957 0.452\n263SOL H 0.430 0.971 0.484\n263SOL H 0.463 0.991 0.432\n264SOL O 0.132 0.637 0.452\n264SOL H 0.098 0.658 0.463\n264SOL H 0.118 0.606 0.430\n265SOL O 0.132 0.958 0.132\n265SOL H 0.098 0.937 0.121\n265SOL H 0.118 0.989 0.153\n266SOL O 0.452 0.637 0.132\n266SOL H 0.463 0.606 0.155\n266SOL H 0.429 0.660 0.156\n267SOL O 0.199 0.705 0.385\n267SOL H 0.231 0.716 0.408\n267SOL H 0.176 0.681 0.409\n268SOL O 0.385 0.890 0.385\n268SOL H 0.374 0.914 0.353\n268SOL H 0.408 0.915 0.407\n269SOL O 0.199 0.890 0.199\n269SOL H 0.231 0.879 0.176\n269SOL H 0.176 0.914 0.175\n270SOL O 0.385 0.704 0.199\n270SOL H 0.405 0.738 0.210\n270SOL H 0.353 0.718 0.177\n271SOL O 0.039 0.798 0.166\n271SOL H 0.040 0.764 0.143\n271SOL H 0.073 0.796 0.188\n272SOL O 0.039 0.798 0.418\n272SOL H 0.007 0.798 0.394\n272SOL H 0.071 0.798 0.394\n273SOL O 0.165 0.545 0.292\n273SOL H 0.143 0.546 0.326\n273SOL H 0.188 0.579 0.293\n274SOL O 0.418 0.545 0.292\n274SOL H 0.441 0.546 0.258\n274SOL H 0.396 0.579 0.291\n275SOL O 0.292 0.671 0.039\n275SOL H 0.259 0.647 0.039\n275SOL H 0.324 0.647 0.039\n276SOL O 0.292 0.924 0.039\n276SOL H 0.291 0.902 0.005\n276SOL H 0.258 0.947 0.037\n277SOL O 0.545 0.192 0.099\n277SOL H 0.545 0.199 0.059\n277SOL H 0.545 0.229 0.116\n278SOL O 0.545 0.392 0.099\n278SOL H 0.545 0.384 0.059\n278SOL H 0.545 0.354 0.116\n279SOL O 0.544 0.195 0.485\n279SOL H 0.544 0.228 0.462\n279SOL H 0.512 0.173 0.473\n280SOL O 0.545 0.389 0.485\n280SOL H 0.577 0.411 0.474\n280SOL H 0.513 0.411 0.473\n281SOL O 0.603 0.039 0.195\n281SOL H 0.614 0.072 0.173\n281SOL H 0.563 0.040 0.195\n282SOL O 0.604 0.038 0.389\n282SOL H 0.616 0.006 0.410\n282SOL H 0.628 0.038 0.355\n283SOL O 0.991 0.039 0.195\n283SOL H 0.980 0.072 0.174\n283SOL H 0.968 0.040 0.228\n284SOL O 0.992 0.039 0.389\n284SOL H 0.981 0.072 0.411\n284SOL H 1.033 0.040 0.389\n285SOL O 0.701 0.098 0.039\n285SOL H 0.700 0.057 0.040\n285SOL H 0.679 0.109 0.072\n286SOL O 0.896 0.098 0.039\n286SOL H 0.917 0.110 0.006\n286SOL H 0.917 0.110 0.071\n287SOL O 0.701 0.486 0.039\n287SOL H 0.700 0.527 0.040\n287SOL H 0.679 0.475 0.071\n288SOL O 0.898 0.484 0.039\n288SOL H 0.890 0.524 0.039\n288SOL H 0.861 0.467 0.039\n289SOL O 0.954 0.292 0.233\n289SOL H 0.954 0.293 0.274\n289SOL H 0.932 0.325 0.222\n290SOL O 0.642 0.292 0.233\n290SOL H 0.642 0.291 0.274\n290SOL H 0.664 0.259 0.222\n291SOL O 0.954 0.292 0.352\n291SOL H 0.932 0.325 0.363\n291SOL H 0.987 0.293 0.375\n292SOL O 0.642 0.292 0.352\n292SOL H 0.663 0.259 0.363\n292SOL H 0.608 0.291 0.375\n293SOL O 0.797 0.352 0.136\n293SOL H 0.765 0.363 0.157\n293SOL H 0.797 0.375 0.102\n294SOL O 0.798 0.352 0.445\n294SOL H 0.798 0.369 0.482\n294SOL H 0.798 0.312 0.452\n295SOL O 0.798 0.231 0.139\n295SOL H 0.798 0.271 0.132\n295SOL H 0.798 0.214 0.102\n296SOL O 0.798 0.232 0.448\n296SOL H 0.798 0.209 0.481\n296SOL H 0.830 0.220 0.427\n297SOL O 0.856 0.448 0.292\n297SOL H 0.816 0.447 0.291\n297SOL H 0.867 0.426 0.259\n298SOL O 0.857 0.137 0.292\n298SOL H 0.869 0.159 0.324\n298SOL H 0.869 0.159 0.259\n299SOL O 0.738 0.448 0.292\n299SOL H 0.715 0.481 0.293\n299SOL H 0.726 0.427 0.325\n300SOL O 0.739 0.136 0.292\n300SOL H 0.780 0.136 0.293\n300SOL H 0.728 0.158 0.324\n301SOL O 0.637 0.132 0.132\n301SOL H 0.605 0.155 0.121\n301SOL H 0.660 0.156 0.155\n302SOL O 0.958 0.452 0.132\n302SOL H 0.969 0.485 0.153\n302SOL H 0.936 0.465 0.100\n303SOL O 0.958 0.132 0.452\n303SOL H 0.937 0.121 0.485\n303SOL H 0.989 0.153 0.466\n304SOL O 0.637 0.452 0.452\n304SOL H 0.627 0.484 0.429\n304SOL H 0.661 0.429 0.427\n305SOL O 0.890 0.385 0.199\n305SOL H 0.859 0.374 0.176\n305SOL H 0.914 0.408 0.175\n306SOL O 0.704 0.199 0.199\n306SOL H 0.738 0.210 0.179\n306SOL H 0.718 0.177 0.231\n307SOL O 0.890 0.199 0.385\n307SOL H 0.914 0.231 0.374\n307SOL H 0.915 0.176 0.408\n308SOL O 0.705 0.385 0.385\n308SOL H 0.684 0.352 0.374\n308SOL H 0.736 0.372 0.406\n309SOL O 0.958 0.452 0.452\n309SOL H 0.935 0.427 0.429\n309SOL H 0.990 0.429 0.463\n310SOL O 0.638 0.132 0.452\n310SOL H 0.604 0.152 0.463\n310SOL H 0.624 0.100 0.430\n311SOL O 0.637 0.452 0.132\n311SOL H 0.605 0.429 0.121\n311SOL H 0.660 0.427 0.155\n312SOL O 0.958 0.132 0.132\n312SOL H 0.935 0.156 0.155\n312SOL H 0.990 0.155 0.121\n313SOL O 0.705 0.199 0.385\n313SOL H 0.737 0.210 0.408\n313SOL H 0.682 0.176 0.409\n314SOL O 0.891 0.385 0.385\n314SOL H 0.857 0.374 0.406\n314SOL H 0.877 0.406 0.353\n315SOL O 0.705 0.385 0.199\n315SOL H 0.683 0.353 0.213\n315SOL H 0.716 0.406 0.232\n316SOL O 0.890 0.199 0.199\n316SOL H 0.911 0.232 0.210\n316SOL H 0.859 0.212 0.177\n317SOL O 0.545 0.292 0.165\n317SOL H 0.577 0.292 0.190\n317SOL H 0.513 0.292 0.190\n318SOL O 0.545 0.292 0.418\n318SOL H 0.543 0.326 0.441\n318SOL H 0.511 0.293 0.395\n319SOL O 0.671 0.039 0.292\n319SOL H 0.648 0.040 0.258\n319SOL H 0.694 0.073 0.291\n320SOL O 0.924 0.039 0.292\n320SOL H 0.947 0.040 0.326\n320SOL H 0.902 0.073 0.293\n321SOL O 0.798 0.166 0.039\n321SOL H 0.766 0.141 0.039\n321SOL H 0.830 0.141 0.039\n322SOL O 0.798 0.418 0.039\n322SOL H 0.796 0.395 0.005\n322SOL H 0.764 0.441 0.037\n323SOL O 0.039 0.195 0.098\n323SOL H 0.072 0.173 0.110\n323SOL H 0.040 0.228 0.122\n324SOL O 0.039 0.389 0.098\n324SOL H 0.006 0.411 0.110\n324SOL H 0.071 0.411 0.110\n325SOL O 0.039 0.192 0.484\n325SOL H 0.039 0.199 0.524\n325SOL H 0.039 0.229 0.467\n326SOL O 0.039 0.389 0.486\n326SOL H 0.040 0.389 0.527\n326SOL H 0.072 0.411 0.475\n327SOL O 0.099 0.039 0.192\n327SOL H 0.059 0.039 0.199\n327SOL H 0.116 0.039 0.229\n328SOL O 0.098 0.038 0.389\n328SOL H 0.110 0.006 0.410\n328SOL H 0.122 0.038 0.355\n329SOL O 0.485 0.039 0.194\n329SOL H 0.473 0.006 0.173\n329SOL H 0.473 0.071 0.173\n330SOL O 0.484 0.039 0.392\n330SOL H 0.524 0.039 0.384\n330SOL H 0.467 0.039 0.354\n331SOL O 0.195 0.098 0.038\n331SOL H 0.173 0.110 0.006\n331SOL H 0.228 0.122 0.038\n332SOL O 0.389 0.098 0.039\n332SOL H 0.389 0.057 0.040\n332SOL H 0.411 0.109 0.072\n333SOL O 0.194 0.485 0.039\n333SOL H 0.173 0.473 0.006\n333SOL H 0.173 0.473 0.071\n334SOL O 0.389 0.485 0.038\n334SOL H 0.410 0.473 0.006\n334SOL H 0.355 0.462 0.038\n335SOL O 0.448 0.292 0.233\n335SOL H 0.447 0.291 0.274\n335SOL H 0.426 0.259 0.222\n336SOL O 0.136 0.292 0.232\n336SOL H 0.157 0.259 0.220\n336SOL H 0.102 0.291 0.209\n337SOL O 0.447 0.292 0.352\n337SOL H 0.426 0.259 0.363\n337SOL H 0.425 0.324 0.363\n338SOL O 0.136 0.292 0.351\n338SOL H 0.136 0.291 0.310\n338SOL H 0.158 0.259 0.362\n339SOL O 0.292 0.352 0.136\n339SOL H 0.324 0.363 0.159\n339SOL H 0.259 0.363 0.159\n340SOL O 0.292 0.352 0.448\n340SOL H 0.291 0.375 0.481\n340SOL H 0.259 0.363 0.427\n341SOL O 0.292 0.231 0.139\n341SOL H 0.292 0.271 0.132\n341SOL H 0.292 0.214 0.102\n342SOL O 0.292 0.231 0.445\n342SOL H 0.292 0.214 0.482\n342SOL H 0.292 0.271 0.452\n343SOL O 0.352 0.448 0.292\n343SOL H 0.375 0.481 0.291\n343SOL H 0.363 0.427 0.259\n344SOL O 0.352 0.139 0.292\n344SOL H 0.312 0.132 0.292\n344SOL H 0.369 0.102 0.292\n345SOL O 0.231 0.445 0.292\n345SOL H 0.214 0.482 0.292\n345SOL H 0.271 0.452 0.292\n346SOL O 0.232 0.136 0.292\n346SOL H 0.220 0.159 0.259\n346SOL H 0.220 0.159 0.324\n347SOL O 0.132 0.132 0.132\n347SOL H 0.121 0.098 0.153\n347SOL H 0.153 0.118 0.100\n348SOL O 0.452 0.452 0.132\n348SOL H 0.484 0.430 0.118\n348SOL H 0.431 0.463 0.098\n349SOL O 0.452 0.132 0.452\n349SOL H 0.463 0.100 0.429\n349SOL H 0.429 0.154 0.427\n350SOL O 0.132 0.452 0.452\n350SOL H 0.155 0.463 0.484\n350SOL H 0.156 0.429 0.429\n351SOL O 0.385 0.385 0.199\n351SOL H 0.408 0.353 0.210\n351SOL H 0.409 0.407 0.176\n352SOL O 0.199 0.199 0.199\n352SOL H 0.231 0.210 0.176\n352SOL H 0.176 0.176 0.175\n353SOL O 0.385 0.199 0.385\n353SOL H 0.352 0.210 0.406\n353SOL H 0.372 0.177 0.353\n354SOL O 0.199 0.385 0.385\n354SOL H 0.177 0.353 0.371\n354SOL H 0.210 0.406 0.352\n355SOL O 0.452 0.452 0.452\n355SOL H 0.463 0.484 0.429\n355SOL H 0.429 0.429 0.427\n356SOL O 0.132 0.132 0.452\n356SOL H 0.098 0.153 0.463\n356SOL H 0.118 0.100 0.430\n357SOL O 0.132 0.452 0.132\n357SOL H 0.121 0.484 0.155\n357SOL H 0.155 0.429 0.156\n358SOL O 0.452 0.132 0.132\n358SOL H 0.484 0.155 0.121\n358SOL H 0.429 0.156 0.155\n359SOL O 0.199 0.199 0.385\n359SOL H 0.231 0.210 0.408\n359SOL H 0.176 0.176 0.409\n360SOL O 0.385 0.385 0.385\n360SOL H 0.352 0.374 0.406\n360SOL H 0.371 0.406 0.353\n361SOL O 0.199 0.385 0.199\n361SOL H 0.177 0.353 0.213\n361SOL H 0.210 0.406 0.232\n362SOL O 0.385 0.199 0.199\n362SOL H 0.352 0.210 0.178\n362SOL H 0.372 0.177 0.231\n363SOL O 0.039 0.292 0.165\n363SOL H 0.005 0.293 0.188\n363SOL H 0.037 0.326 0.143\n364SOL O 0.039 0.292 0.418\n364SOL H 0.040 0.326 0.441\n364SOL H 0.072 0.293 0.396\n365SOL O 0.165 0.039 0.292\n365SOL H 0.190 0.007 0.292\n365SOL H 0.190 0.071 0.292\n366SOL O 0.418 0.039 0.292\n366SOL H 0.395 0.005 0.291\n366SOL H 0.441 0.037 0.258\n367SOL O 0.292 0.165 0.039\n367SOL H 0.293 0.188 0.005\n367SOL H 0.326 0.143 0.037\n368SOL O 0.292 0.418 0.039\n368SOL H 0.258 0.441 0.040\n368SOL H 0.291 0.396 0.073\n';var cell=ChemDoodle.readCIF(Fmol, 1,1,1);Mol2.loadContent([cell.molecule], [cell.unitCell]);Mol2.startAnimation();Mol2.stopAnimation();var d2r=Math.PI/180, mat2 = [];ChemDoodle.lib.mat4.identity(mat2);ChemDoodle.lib.mat4.rotate(mat2,0*d2r,[1,0,0]);ChemDoodle.lib.mat4.rotate(mat2,0*d2r,[0,1,0]);ChemDoodle.lib.mat4.rotate(mat2,0*d2r,[0,0,1]);ChemDoodle.lib.mat4.multiply(Mol2.rotationMatrix, mat2);Mol2.repaint();var $=function(id){return document.getElementById(id)};function setProj2(yesPers){Mol2.specs.projectionPerspective_3D = yesPers;Mol2.setupScene();Mol2.repaint()}function setModel2(model){Mol2.specs.set3DRepresentation(model);Mol2.setupScene();Mol2.repaint()}function setSupercell2(){var cell=ChemDoodle.readCIF(Fmol, $("Mol2x").value, $("Mol2y").value, $("Mol2z").value);Mol2.loadContent([cell.molecule], [cell.unitCell]);Mol2.repaint()}</script><br><span class="meta">视图: <input type="radio" name="group2" onclick="setProj2(true)">投影 <input type="radio" name="group2" onclick="setProj2(false)" checked="">正交<br>模型: <input type="radio" name="model" onclick="setModel2('Ball and Stick')" checked="">球棍 <input type="radio" name="model" onclick="setModel2('van der Waals Spheres')">范德华球 <input type="radio" name="model" onclick="setModel2('Stick')">棍状 <input type="radio" name="model" onclick="setModel2('Wireframe')">线框 <input type="radio" name="model" onclick="setModel2('Line')">线型 <input type="checkbox" onclick="Mol2.specs.atoms_displayLabels_3D=this.checked;Mol2.repaint()">名称<br>超晶胞: X <input type="text" style="width:20px;" id="Mol2x" value="1"> Y <input type="text" style="width:20px;" id="Mol2y" value="1"> Z <input type="text" style="width:20px;" id="Mol2z" value="1"> <input type="button" value="创建" onclick="setSupercell2()"><br>左键: 转动 滚轮: 缩放 双击: 自动旋转开关 Alt+左键: 移动</span><br><figurecaption>Fig.2</figurecaption></figure>我们选取结构中的17号氧原子进行计算(坐标14.180 12.790 16.450
). 使用Materials Studio软件测量得到周围原子与17号氧原子所成角度和二面角大小.
与示例1相同, 修改AOP.bsh
脚本中的中心原子编号得到计算结果, 并将其与软件测量结果进行对比.
$\q$, $\f$, F3和F4的分布如下
gro
文件中坐标的精度有限, 对某些原子的坐标有时候会与原始文件中的存在差别. 如果gro
文件中输出的坐标位数多一些就不会存在这些问题(新版本GROMACS不再支持变精度的gro文件). 因此, 我们在验证过程中可以使用gmx trjconv
输出的构型文件来代替原始的构型文件.AOP.bsh
脚本计算整个模拟时间段内序参数的变化, 可能需要很长时间. 简单的解决方法是采用并行, 将轨迹复制多份, 每份只计算一段时间的序参数, 然后再合并结果.使用GROMACS自带的spc216.gro
作为初始构型, 采用TIP3P水模型运行200 ps模拟.
计算得到水合物生成过程中F3/F4的变化如下图.
可见, 我们得到的F3在0.1附近, 这与前面所讲的一致.
$\q$, $\f$, F3和F4的分布如下
脚本写起来简单, 因此上面这种组合GROMACS工具的方法有其方便之处, 但如果要计算的数据太过复杂, 脚本处理起来就比较麻烦, 而且计算速度也很慢, 因此, 组合脚本的方法不适用于计算太复杂的数据, 处理长时间的轨迹. 对于我们的情况, F3使用脚本计算还算简单, 但F4使用脚本计算就已经很麻烦了. 是时候考虑下使用编译型语言来提高计算速度了.
简单的提速方法是不再调用GROMACS工具来计算角度, 距离等, 而是直接使用脚本所需要的数据. 但这种方法受限于脚本的运行速度, 效果可能仍然不够理想. 因此, 我们需要考虑使用编译型语言来处理. 推荐使用C/C++, 或Fortran.
因为我熟悉Fortran, 所以就使用Fortran编写了序参数的计算代码.
先编译好AOP.f90
代码, 并在gmx trjconv
命令中进行调用
注意, index.ndx
中氧原子组的名称必须和程序参数指定的一致. AOP.exe
此项默认氧原子的名称为OW
.
此外, 你也可以先将要计算的水分子的坐标抽取出来, 保存为traj.gro
, 然后直接运行程序. 但在这种情况下需要使用traj.gro
重新创建index.ndx
, 因为traj.gro
中的氧原子编号可能会与原始轨迹中的不同.
我们通过从甲烷水合物轨迹中选取的几帧构型来验证bash脚本和fortran代码计算结果的一致性. 结果如下图所示.
<table id='tab-5'><caption></caption> <tr> <th rowspan="1" colspan="1" style="text-align:center;">序参数</th> <th rowspan="1" colspan="1" style="text-align:center;">方法</th> <th rowspan="1" colspan="1" style="text-align:center;">0 ps</th> <th rowspan="1" colspan="1" style="text-align:center;">5 ps</th> <th rowspan="1" colspan="1" style="text-align:center;">10 ps</th> <th rowspan="1" colspan="1" style="text-align:center;">15 ps</th> <th rowspan="1" colspan="1" style="text-align:center;">20 ps</th> <th rowspan="1" colspan="1" style="text-align:center;">25 ps</th> <th rowspan="1" colspan="1" style="text-align:center;">30 ps</th> </tr> <tr> <td rowspan="2" colspan="1" style="text-align:center;">F3</td> <td rowspan="1" colspan="1" style="text-align:center;">Bash</td> <td rowspan="1" colspan="1" style="text-align:center;">0.0566</td> <td rowspan="1" colspan="1" style="text-align:center;">0.0565</td> <td rowspan="1" colspan="1" style="text-align:center;">0.0567</td> <td rowspan="1" colspan="1" style="text-align:center;">0.0570</td> <td rowspan="1" colspan="1" style="text-align:center;">0.0569</td> <td rowspan="1" colspan="1" style="text-align:center;">0.0578</td> <td rowspan="1" colspan="1" style="text-align:center;">0.0591</td> </tr> <tr> <td rowspan="1" colspan="1" style="text-align:center;">Fortran</td> <td rowspan="1" colspan="1" style="text-align:center;">0.0566</td> <td rowspan="1" colspan="1" style="text-align:center;">0.0565</td> <td rowspan="1" colspan="1" style="text-align:center;">0.0567</td> <td rowspan="1" colspan="1" style="text-align:center;">0.0570</td> <td rowspan="1" colspan="1" style="text-align:center;">0.0569</td> <td rowspan="1" colspan="1" style="text-align:center;">0.0578</td> <td rowspan="1" colspan="1" style="text-align:center;">0.0591</td> </tr> <tr> <td rowspan="2" colspan="1" style="text-align:center;">F4</td> <td rowspan="1" colspan="1" style="text-align:center;">Bash</td> <td rowspan="1" colspan="1" style="text-align:center;">-0.0478</td> <td rowspan="1" colspan="1" style="text-align:center;">-0.0283</td> <td rowspan="1" colspan="1" style="text-align:center;">-0.0338</td> <td rowspan="1" colspan="1" style="text-align:center;">-0.0265</td> <td rowspan="1" colspan="1" style="text-align:center;">-0.0149</td> <td rowspan="1" colspan="1" style="text-align:center;">-0.0281</td> <td rowspan="1" colspan="1" style="text-align:center;">-0.0189</td> </tr> <tr> <td rowspan="1" colspan="1" style="text-align:center;">Fortran</td> <td rowspan="1" colspan="1" style="text-align:center;">-0.0479</td> <td rowspan="1" colspan="1" style="text-align:center;">-0.0283</td> <td rowspan="1" colspan="1" style="text-align:center;">-0.0338</td> <td rowspan="1" colspan="1" style="text-align:center;">-0.0264</td> <td rowspan="1" colspan="1" style="text-align:center;">-0.0149</td> <td rowspan="1" colspan="1" style="text-align:center;">-0.0281</td> <td rowspan="1" colspan="1" style="text-align:center;">-0.0190</td> </tr> </table>可见二者所得结果一致, 这说明代码中涉及PBC处理, 距离和角度计算的部分是正确的.
$\q$, $\f$, F3和F4的分布如下
在此示例中, 每帧构型中氧原子的数目为1610, 水分子的总原子数目为6440. 在测试机器上(操作系统: Windows10; 内存: 8GB; 处理器:Inter(R) Core(TM)i5-7400 CPU @3.0GHZ), 脚本处理一帧构型所需的时间为5分钟, 而Fortran代码仅需1秒左右. 可以看出Fortran的计算速度远远高于脚本.
在水合物计算相关的文章中, 通常根据F3/F4随时间的变化规律的来分析水合物的生成过程. 在此, 我们建立如下图所示的初始结构, 盒子左边为2×2×2的甲烷水合物晶胞, 右边为甲烷-水混合溶液.
<figure><script>var Mol3=new ChemDoodle.TransformCanvas3D('Mol-3',502,376.5);Mol3.specs.shapes_color = '#fff';Mol3.specs.backgroundColor = 'black';Mol3.specs.projectionPerspective_3D = false;Mol3.specs.set3DRepresentation('Ball and Stick');Mol3.specs.compass_display = true;//Mol3.specs.atoms_resolution_3D = 15; //Mol3.specs.bonds_resolution_3D = 15; Mol3.handle = null;Mol3.timeout = 15;Mol3.specs.crystals_unitCellLineWidth = 1.5;Mol3.startAnimation = ChemDoodle._AnimatorCanvas.prototype.startAnimation;Mol3.stopAnimation = ChemDoodle._AnimatorCanvas.prototype.stopAnimation;Mol3.isRunning = ChemDoodle._AnimatorCanvas.prototype.isRunning;Mol3.dblclick = ChemDoodle.RotatorCanvas.prototype.dblclick;Mol3.nextFrame = function(delta){var matrix = [];ChemDoodle.lib.mat4.identity(matrix);var change = delta*Math.PI/15000;ChemDoodle.lib.mat4.rotate(matrix,change,[1,0,0]);ChemDoodle.lib.mat4.rotate(matrix,change,[0,1,0]);ChemDoodle.lib.mat4.rotate(matrix,change,[0,0,1]);ChemDoodle.lib.mat4.multiply(this.rotationMatrix, matrix)};var Fmol='data_Built_with_Packmol\n4480\n_symmetry_space_group_name_H-M \'P 1\'\n_cell_length_a 23.75\n_cell_length_b 23.75\n_cell_length_c 60.65\n_cell_angle_alpha 90\n_cell_angle_beta 90\n_cell_angle_gamma 90\nloop_\n_atom_site_label\n_atom_site_type_symbol\n_atom_site_fract_x\n_atom_site_fract_y\n_atom_site_fract_z\n1SOL O 0.538 0.692 0.234\n1SOL H 0.570 0.670 0.238\n1SOL H 0.507 0.670 0.238\n2SOL O 0.538 0.887 0.234\n2SOL H 0.538 0.851 0.241\n2SOL H 0.538 0.880 0.219\n3SOL O 0.538 0.692 0.384\n3SOL H 0.570 0.670 0.379\n3SOL H 0.507 0.670 0.379\n4SOL O 0.539 0.885 0.384\n4SOL H 0.571 0.906 0.380\n4SOL H 0.539 0.885 0.400\n5SOL O 0.597 0.539 0.271\n5SOL H 0.609 0.571 0.263\n5SOL H 0.620 0.539 0.284\n6SOL O 0.598 0.538 0.347\n6SOL H 0.615 0.538 0.333\n6SOL H 0.559 0.538 0.345\n7SOL O 0.979 0.538 0.271\n7SOL H 0.956 0.537 0.284\n7SOL H 0.968 0.506 0.263\n8SOL O 0.979 0.538 0.346\n8SOL H 0.968 0.570 0.355\n8SOL H 0.968 0.507 0.355\n9SOL O 0.690 0.598 0.211\n9SOL H 0.726 0.615 0.211\n9SOL H 0.697 0.558 0.211\n10SOL O 0.887 0.598 0.211\n10SOL H 0.851 0.615 0.211\n10SOL H 0.880 0.558 0.211\n11SOL O 0.692 0.980 0.211\n11SOL H 0.692 1.021 0.210\n11SOL H 0.671 0.969 0.198\n12SOL O 0.884 0.979 0.211\n12SOL H 0.905 0.968 0.223\n12SOL H 0.851 0.956 0.211\n13SOL O 0.943 0.789 0.286\n13SOL H 0.922 0.821 0.281\n13SOL H 0.976 0.789 0.277\n14SOL O 0.635 0.789 0.286\n14SOL H 0.635 0.789 0.302\n14SOL H 0.656 0.821 0.282\n15SOL O 0.940 0.788 0.332\n15SOL H 0.947 0.788 0.317\n15SOL H 0.976 0.788 0.339\n16SOL O 0.634 0.789 0.332\n16SOL H 0.655 0.821 0.336\n16SOL H 0.601 0.789 0.341\n17SOL O 0.789 0.846 0.248\n17SOL H 0.789 0.806 0.248\n17SOL H 0.821 0.857 0.257\n18SOL O 0.789 0.846 0.369\n18SOL H 0.789 0.806 0.369\n18SOL H 0.821 0.857 0.361\n19SOL O 0.788 0.730 0.249\n19SOL H 0.757 0.718 0.257\n19SOL H 0.820 0.718 0.257\n20SOL O 0.788 0.730 0.369\n20SOL H 0.820 0.718 0.360\n20SOL H 0.757 0.718 0.360\n21SOL O 0.846 0.942 0.309\n21SOL H 0.806 0.942 0.309\n21SOL H 0.857 0.921 0.321\n22SOL O 0.847 0.634 0.308\n22SOL H 0.859 0.655 0.296\n22SOL H 0.870 0.601 0.308\n23SOL O 0.730 0.943 0.308\n23SOL H 0.718 0.921 0.296\n23SOL H 0.706 0.976 0.308\n24SOL O 0.731 0.635 0.309\n24SOL H 0.771 0.635 0.309\n24SOL H 0.720 0.656 0.321\n25SOL O 0.630 0.630 0.247\n25SOL H 0.653 0.619 0.234\n25SOL H 0.654 0.653 0.256\n26SOL O 0.946 0.947 0.247\n26SOL H 0.979 0.926 0.243\n26SOL H 0.960 0.978 0.255\n27SOL O 0.947 0.630 0.371\n27SOL H 0.923 0.653 0.362\n27SOL H 0.924 0.619 0.383\n28SOL O 0.630 0.947 0.371\n28SOL H 0.654 0.924 0.362\n28SOL H 0.653 0.957 0.383\n29SOL O 0.880 0.880 0.273\n29SOL H 0.869 0.903 0.285\n29SOL H 0.903 0.904 0.264\n30SOL O 0.696 0.696 0.273\n30SOL H 0.675 0.728 0.278\n30SOL H 0.707 0.676 0.286\n31SOL O 0.880 0.696 0.345\n31SOL H 0.901 0.728 0.339\n31SOL H 0.869 0.676 0.332\n32SOL O 0.696 0.880 0.345\n32SOL H 0.729 0.869 0.353\n32SOL H 0.710 0.901 0.332\n33SOL O 0.947 0.947 0.370\n33SOL H 0.978 0.925 0.376\n33SOL H 0.926 0.957 0.384\n34SOL O 0.630 0.630 0.371\n34SOL H 0.620 0.597 0.363\n34SOL H 0.651 0.617 0.383\n35SOL O 0.630 0.947 0.247\n35SOL H 0.599 0.925 0.242\n35SOL H 0.650 0.957 0.234\n36SOL O 0.947 0.630 0.247\n36SOL H 0.957 0.597 0.255\n36SOL H 0.925 0.617 0.234\n37SOL O 0.696 0.696 0.345\n37SOL H 0.674 0.728 0.340\n37SOL H 0.672 0.674 0.354\n38SOL O 0.880 0.880 0.345\n38SOL H 0.903 0.849 0.340\n38SOL H 0.904 0.902 0.354\n39SOL O 0.696 0.880 0.273\n39SOL H 0.728 0.869 0.264\n39SOL H 0.674 0.903 0.263\n40SOL O 0.880 0.696 0.273\n40SOL H 0.903 0.728 0.277\n40SOL H 0.904 0.674 0.264\n41SOL O 0.538 0.788 0.260\n41SOL H 0.539 0.755 0.251\n41SOL H 0.572 0.787 0.269\n42SOL O 0.538 0.788 0.358\n42SOL H 0.538 0.757 0.367\n42SOL H 0.538 0.820 0.367\n43SOL O 0.663 0.538 0.309\n43SOL H 0.688 0.570 0.309\n43SOL H 0.688 0.507 0.309\n44SOL O 0.913 0.538 0.309\n44SOL H 0.936 0.537 0.322\n44SOL H 0.891 0.505 0.309\n45SOL O 0.788 0.663 0.211\n45SOL H 0.788 0.688 0.223\n45SOL H 0.788 0.688 0.198\n46SOL O 0.788 0.913 0.211\n46SOL H 0.755 0.936 0.211\n46SOL H 0.787 0.891 0.224\n47SOL O 0.038 0.693 0.234\n47SOL H 0.006 0.672 0.238\n47SOL H 0.037 0.726 0.243\n48SOL O 0.038 0.884 0.233\n48SOL H 0.071 0.906 0.238\n48SOL H 0.039 0.885 0.218\n49SOL O 0.038 0.693 0.384\n49SOL H 0.006 0.672 0.379\n49SOL H 0.037 0.726 0.374\n50SOL O 0.038 0.885 0.384\n50SOL H 0.071 0.906 0.380\n50SOL H 0.039 0.885 0.400\n51SOL O 0.096 0.539 0.271\n51SOL H 0.056 0.539 0.271\n51SOL H 0.107 0.571 0.263\n52SOL O 0.096 0.539 0.346\n52SOL H 0.056 0.539 0.346\n52SOL H 0.107 0.571 0.355\n53SOL O 0.480 0.538 0.271\n53SOL H 0.520 0.537 0.271\n53SOL H 0.469 0.506 0.263\n54SOL O 0.480 0.539 0.346\n54SOL H 0.468 0.571 0.354\n54SOL H 0.456 0.539 0.333\n55SOL O 0.193 0.597 0.211\n55SOL H 0.226 0.621 0.210\n55SOL H 0.171 0.609 0.198\n56SOL O 0.384 0.597 0.211\n56SOL H 0.405 0.609 0.223\n56SOL H 0.351 0.621 0.211\n57SOL O 0.193 0.979 0.211\n57SOL H 0.226 0.956 0.210\n57SOL H 0.172 0.968 0.198\n58SOL O 0.384 0.980 0.211\n58SOL H 0.406 0.970 0.223\n58SOL H 0.384 1.021 0.211\n59SOL O 0.440 0.788 0.285\n59SOL H 0.476 0.788 0.279\n59SOL H 0.447 0.788 0.301\n60SOL O 0.135 0.788 0.286\n60SOL H 0.157 0.757 0.281\n60SOL H 0.157 0.820 0.281\n61SOL O 0.443 0.789 0.332\n61SOL H 0.476 0.789 0.341\n61SOL H 0.422 0.821 0.336\n62SOL O 0.134 0.788 0.331\n62SOL H 0.135 0.787 0.316\n62SOL H 0.156 0.756 0.336\n63SOL O 0.288 0.846 0.248\n63SOL H 0.289 0.806 0.248\n63SOL H 0.321 0.857 0.257\n64SOL O 0.288 0.847 0.369\n64SOL H 0.256 0.859 0.360\n64SOL H 0.320 0.859 0.360\n65SOL O 0.288 0.730 0.249\n65SOL H 0.256 0.718 0.257\n65SOL H 0.320 0.718 0.257\n66SOL O 0.288 0.728 0.368\n66SOL H 0.288 0.768 0.371\n66SOL H 0.288 0.712 0.382\n67SOL O 0.347 0.943 0.308\n67SOL H 0.359 0.921 0.296\n67SOL H 0.371 0.976 0.308\n68SOL O 0.347 0.634 0.309\n68SOL H 0.359 0.655 0.321\n68SOL H 0.371 0.601 0.309\n69SOL O 0.230 0.942 0.309\n69SOL H 0.270 0.942 0.309\n69SOL H 0.219 0.921 0.321\n70SOL O 0.229 0.637 0.309\n70SOL H 0.268 0.630 0.309\n70SOL H 0.212 0.600 0.309\n71SOL O 0.130 0.630 0.247\n71SOL H 0.099 0.651 0.242\n71SOL H 0.151 0.619 0.234\n72SOL O 0.446 0.947 0.247\n72SOL H 0.480 0.926 0.243\n72SOL H 0.460 0.978 0.255\n73SOL O 0.447 0.630 0.371\n73SOL H 0.422 0.653 0.362\n73SOL H 0.424 0.619 0.383\n74SOL O 0.130 0.946 0.371\n74SOL H 0.120 0.979 0.363\n74SOL H 0.152 0.960 0.383\n75SOL O 0.380 0.880 0.273\n75SOL H 0.403 0.849 0.277\n75SOL H 0.404 0.903 0.264\n76SOL O 0.197 0.696 0.273\n76SOL H 0.208 0.674 0.285\n76SOL H 0.173 0.672 0.264\n77SOL O 0.381 0.696 0.345\n77SOL H 0.400 0.729 0.340\n77SOL H 0.349 0.710 0.353\n78SOL O 0.197 0.880 0.345\n78SOL H 0.174 0.849 0.340\n78SOL H 0.173 0.902 0.354\n79SOL O 0.446 0.947 0.371\n79SOL H 0.480 0.926 0.375\n79SOL H 0.460 0.978 0.362\n80SOL O 0.130 0.630 0.371\n80SOL H 0.099 0.653 0.375\n80SOL H 0.152 0.654 0.362\n81SOL O 0.130 0.947 0.247\n81SOL H 0.153 0.957 0.234\n81SOL H 0.154 0.924 0.256\n82SOL O 0.447 0.630 0.247\n82SOL H 0.458 0.599 0.256\n82SOL H 0.424 0.652 0.256\n83SOL O 0.196 0.696 0.345\n83SOL H 0.229 0.707 0.353\n83SOL H 0.210 0.675 0.332\n84SOL O 0.380 0.880 0.345\n84SOL H 0.369 0.903 0.332\n84SOL H 0.403 0.904 0.354\n85SOL O 0.196 0.880 0.273\n85SOL H 0.229 0.869 0.265\n85SOL H 0.210 0.901 0.285\n86SOL O 0.380 0.696 0.273\n86SOL H 0.401 0.728 0.278\n86SOL H 0.370 0.676 0.286\n87SOL O 0.038 0.788 0.260\n87SOL H 0.040 0.822 0.251\n87SOL H 0.072 0.789 0.268\n88SOL O 0.038 0.788 0.358\n88SOL H 0.040 0.821 0.366\n88SOL H 0.072 0.789 0.349\n89SOL O 0.164 0.538 0.309\n89SOL H 0.139 0.538 0.296\n89SOL H 0.139 0.538 0.321\n90SOL O 0.413 0.538 0.309\n90SOL H 0.436 0.537 0.296\n90SOL H 0.391 0.505 0.308\n91SOL O 0.288 0.663 0.211\n91SOL H 0.288 0.688 0.223\n91SOL H 0.288 0.688 0.198\n92SOL O 0.288 0.913 0.211\n92SOL H 0.322 0.936 0.211\n92SOL H 0.290 0.891 0.224\n93SOL O 0.538 0.190 0.234\n93SOL H 0.538 0.226 0.241\n93SOL H 0.538 0.197 0.219\n94SOL O 0.538 0.387 0.234\n94SOL H 0.538 0.351 0.241\n94SOL H 0.538 0.380 0.219\n95SOL O 0.538 0.193 0.384\n95SOL H 0.537 0.226 0.374\n95SOL H 0.506 0.172 0.379\n96SOL O 0.539 0.384 0.384\n96SOL H 0.571 0.405 0.379\n96SOL H 0.539 0.351 0.374\n97SOL O 0.596 0.038 0.271\n97SOL H 0.607 0.006 0.263\n97SOL H 0.556 0.037 0.271\n98SOL O 0.597 0.038 0.347\n98SOL H 0.609 0.006 0.355\n98SOL H 0.609 0.070 0.355\n99SOL O 0.979 0.038 0.270\n99SOL H 0.962 0.038 0.284\n99SOL H 1.018 0.038 0.273\n100SOL O 0.979 0.038 0.346\n100SOL H 0.968 0.006 0.354\n100SOL H 0.956 0.037 0.333\n101SOL O 0.692 0.097 0.211\n101SOL H 0.670 0.109 0.223\n101SOL H 0.670 0.109 0.198\n102SOL O 0.885 0.096 0.211\n102SOL H 0.885 0.056 0.211\n102SOL H 0.906 0.107 0.223\n103SOL O 0.692 0.480 0.211\n103SOL H 0.670 0.468 0.223\n103SOL H 0.670 0.468 0.198\n104SOL O 0.884 0.480 0.211\n104SOL H 0.851 0.456 0.210\n104SOL H 0.905 0.468 0.198\n105SOL O 0.942 0.288 0.286\n105SOL H 0.920 0.320 0.281\n105SOL H 0.920 0.256 0.281\n106SOL O 0.635 0.288 0.286\n106SOL H 0.656 0.256 0.281\n106SOL H 0.656 0.320 0.281\n107SOL O 0.940 0.288 0.332\n107SOL H 0.947 0.288 0.317\n107SOL H 0.976 0.288 0.339\n108SOL O 0.634 0.288 0.331\n108SOL H 0.635 0.288 0.316\n108SOL H 0.656 0.256 0.336\n109SOL O 0.789 0.347 0.248\n109SOL H 0.789 0.307 0.248\n109SOL H 0.821 0.357 0.257\n110SOL O 0.788 0.347 0.369\n110SOL H 0.787 0.307 0.369\n110SOL H 0.756 0.357 0.361\n111SOL O 0.788 0.229 0.248\n111SOL H 0.756 0.218 0.257\n111SOL H 0.787 0.206 0.235\n112SOL O 0.789 0.229 0.369\n112SOL H 0.821 0.218 0.361\n112SOL H 0.789 0.206 0.382\n113SOL O 0.846 0.443 0.309\n113SOL H 0.806 0.442 0.309\n113SOL H 0.857 0.421 0.321\n114SOL O 0.846 0.134 0.309\n114SOL H 0.806 0.135 0.308\n114SOL H 0.857 0.156 0.296\n115SOL O 0.730 0.442 0.309\n115SOL H 0.718 0.421 0.321\n115SOL H 0.718 0.420 0.296\n116SOL O 0.730 0.134 0.309\n116SOL H 0.718 0.155 0.321\n116SOL H 0.707 0.101 0.309\n117SOL O 0.630 0.130 0.247\n117SOL H 0.597 0.151 0.243\n117SOL H 0.617 0.099 0.255\n118SOL O 0.947 0.447 0.247\n118SOL H 0.926 0.457 0.234\n118SOL H 0.978 0.425 0.242\n119SOL O 0.947 0.130 0.371\n119SOL H 0.957 0.099 0.362\n119SOL H 0.923 0.152 0.361\n120SOL O 0.630 0.446 0.371\n120SOL H 0.619 0.480 0.363\n120SOL H 0.651 0.460 0.383\n121SOL O 0.880 0.380 0.273\n121SOL H 0.869 0.403 0.285\n121SOL H 0.903 0.404 0.264\n122SOL O 0.696 0.197 0.273\n122SOL H 0.707 0.174 0.285\n122SOL H 0.673 0.173 0.264\n123SOL O 0.880 0.197 0.345\n123SOL H 0.901 0.228 0.339\n123SOL H 0.869 0.176 0.332\n124SOL O 0.696 0.380 0.345\n124SOL H 0.674 0.349 0.340\n124SOL H 0.672 0.403 0.354\n125SOL O 0.947 0.447 0.371\n125SOL H 0.923 0.424 0.362\n125SOL H 0.924 0.458 0.383\n126SOL O 0.630 0.130 0.370\n126SOL H 0.599 0.152 0.376\n126SOL H 0.650 0.120 0.384\n127SOL O 0.630 0.447 0.247\n127SOL H 0.617 0.478 0.255\n127SOL H 0.597 0.426 0.243\n128SOL O 0.946 0.130 0.247\n128SOL H 0.960 0.099 0.255\n128SOL H 0.979 0.151 0.243\n129SOL O 0.696 0.197 0.345\n129SOL H 0.728 0.208 0.354\n129SOL H 0.674 0.174 0.354\n130SOL O 0.880 0.381 0.345\n130SOL H 0.901 0.347 0.340\n130SOL H 0.849 0.367 0.353\n131SOL O 0.697 0.380 0.273\n131SOL H 0.728 0.369 0.264\n131SOL H 0.674 0.403 0.263\n132SOL O 0.880 0.197 0.273\n132SOL H 0.848 0.208 0.264\n132SOL H 0.903 0.174 0.263\n133SOL O 0.538 0.288 0.260\n133SOL H 0.570 0.288 0.269\n133SOL H 0.507 0.288 0.269\n134SOL O 0.538 0.288 0.358\n134SOL H 0.570 0.288 0.348\n134SOL H 0.507 0.288 0.348\n135SOL O 0.664 0.038 0.309\n135SOL H 0.639 0.038 0.296\n135SOL H 0.639 0.038 0.321\n136SOL O 0.914 0.038 0.309\n136SOL H 0.889 0.007 0.309\n136SOL H 0.889 0.070 0.309\n137SOL O 0.788 0.164 0.211\n137SOL H 0.757 0.139 0.211\n137SOL H 0.820 0.139 0.211\n138SOL O 0.788 0.413 0.211\n138SOL H 0.755 0.436 0.211\n138SOL H 0.787 0.391 0.224\n139SOL O 0.039 0.193 0.234\n139SOL H 0.071 0.172 0.238\n139SOL H 0.039 0.226 0.243\n140SOL O 0.039 0.384 0.234\n140SOL H 0.071 0.405 0.238\n140SOL H 0.039 0.351 0.243\n141SOL O 0.038 0.192 0.384\n141SOL H 0.006 0.171 0.380\n141SOL H 0.037 0.192 0.400\n142SOL O 0.038 0.384 0.384\n142SOL H 0.006 0.406 0.380\n142SOL H 0.037 0.384 0.400\n143SOL O 0.097 0.038 0.271\n143SOL H 0.109 0.006 0.263\n143SOL H 0.120 0.037 0.284\n144SOL O 0.098 0.038 0.347\n144SOL H 0.059 0.038 0.345\n144SOL H 0.115 0.038 0.333\n145SOL O 0.480 0.038 0.271\n145SOL H 0.468 0.071 0.263\n145SOL H 0.456 0.039 0.284\n146SOL O 0.480 0.038 0.346\n146SOL H 0.520 0.039 0.346\n146SOL H 0.469 0.071 0.355\n147SOL O 0.192 0.096 0.211\n147SOL H 0.192 0.056 0.210\n147SOL H 0.171 0.107 0.198\n148SOL O 0.384 0.097 0.211\n148SOL H 0.351 0.120 0.210\n148SOL H 0.405 0.109 0.198\n149SOL O 0.192 0.480 0.211\n149SOL H 0.192 0.520 0.211\n149SOL H 0.171 0.469 0.223\n150SOL O 0.384 0.480 0.211\n150SOL H 0.384 0.520 0.210\n150SOL H 0.406 0.469 0.198\n151SOL O 0.442 0.288 0.286\n151SOL H 0.421 0.320 0.281\n151SOL H 0.421 0.256 0.281\n152SOL O 0.134 0.288 0.286\n152SOL H 0.135 0.288 0.302\n152SOL H 0.156 0.256 0.282\n153SOL O 0.443 0.288 0.331\n153SOL H 0.442 0.289 0.316\n153SOL H 0.421 0.321 0.336\n154SOL O 0.134 0.289 0.332\n154SOL H 0.155 0.321 0.336\n154SOL H 0.101 0.289 0.341\n155SOL O 0.288 0.347 0.248\n155SOL H 0.256 0.359 0.257\n155SOL H 0.288 0.371 0.235\n156SOL O 0.288 0.347 0.369\n156SOL H 0.256 0.359 0.361\n156SOL H 0.288 0.371 0.382\n157SOL O 0.288 0.230 0.248\n157SOL H 0.288 0.270 0.248\n157SOL H 0.256 0.219 0.257\n158SOL O 0.288 0.230 0.369\n158SOL H 0.288 0.270 0.369\n158SOL H 0.256 0.219 0.361\n159SOL O 0.347 0.442 0.309\n159SOL H 0.359 0.420 0.296\n159SOL H 0.359 0.421 0.321\n160SOL O 0.347 0.135 0.309\n160SOL H 0.359 0.157 0.321\n160SOL H 0.359 0.157 0.296\n161SOL O 0.229 0.440 0.309\n161SOL H 0.212 0.476 0.309\n161SOL H 0.268 0.447 0.309\n162SOL O 0.230 0.135 0.309\n162SOL H 0.270 0.135 0.309\n162SOL H 0.219 0.156 0.321\n163SOL O 0.130 0.131 0.247\n163SOL H 0.119 0.097 0.255\n163SOL H 0.152 0.117 0.234\n164SOL O 0.446 0.447 0.247\n164SOL H 0.478 0.425 0.242\n164SOL H 0.427 0.457 0.234\n165SOL O 0.447 0.130 0.371\n165SOL H 0.422 0.153 0.362\n165SOL H 0.424 0.119 0.383\n166SOL O 0.131 0.447 0.371\n166SOL H 0.117 0.478 0.362\n166SOL H 0.097 0.426 0.375\n167SOL O 0.380 0.380 0.273\n167SOL H 0.349 0.369 0.264\n167SOL H 0.403 0.403 0.263\n168SOL O 0.197 0.197 0.273\n168SOL H 0.208 0.174 0.285\n168SOL H 0.174 0.173 0.264\n169SOL O 0.380 0.196 0.345\n169SOL H 0.401 0.229 0.340\n169SOL H 0.349 0.210 0.353\n170SOL O 0.197 0.380 0.345\n170SOL H 0.208 0.403 0.332\n170SOL H 0.174 0.404 0.354\n171SOL O 0.446 0.447 0.371\n171SOL H 0.460 0.478 0.362\n171SOL H 0.480 0.426 0.375\n172SOL O 0.131 0.130 0.371\n172SOL H 0.097 0.151 0.375\n172SOL H 0.117 0.099 0.362\n173SOL O 0.130 0.447 0.247\n173SOL H 0.119 0.478 0.256\n173SOL H 0.153 0.424 0.256\n174SOL O 0.446 0.130 0.247\n174SOL H 0.478 0.152 0.242\n174SOL H 0.427 0.120 0.234\n175SOL O 0.197 0.197 0.345\n175SOL H 0.174 0.228 0.340\n175SOL H 0.173 0.174 0.354\n176SOL O 0.380 0.380 0.345\n176SOL H 0.349 0.369 0.354\n176SOL H 0.403 0.403 0.354\n177SOL O 0.197 0.380 0.273\n177SOL H 0.175 0.349 0.278\n177SOL H 0.207 0.400 0.286\n178SOL O 0.380 0.197 0.273\n178SOL H 0.349 0.208 0.264\n178SOL H 0.403 0.173 0.263\n179SOL O 0.038 0.288 0.260\n179SOL H 0.007 0.288 0.269\n179SOL H 0.070 0.288 0.269\n180SOL O 0.038 0.288 0.358\n180SOL H 0.038 0.256 0.367\n180SOL H 0.038 0.320 0.367\n181SOL O 0.163 0.038 0.309\n181SOL H 0.188 0.007 0.309\n181SOL H 0.188 0.070 0.309\n182SOL O 0.413 0.038 0.309\n182SOL H 0.436 0.040 0.322\n182SOL H 0.391 0.072 0.309\n183SOL O 0.288 0.163 0.211\n183SOL H 0.255 0.141 0.211\n183SOL H 0.287 0.186 0.224\n184SOL O 0.288 0.413 0.211\n184SOL H 0.256 0.438 0.211\n184SOL H 0.320 0.438 0.211\n185SOL O 0.538 0.693 0.038\n185SOL H 0.537 0.692 0.022\n185SOL H 0.506 0.671 0.042\n186SOL O 0.539 0.884 0.038\n186SOL H 0.571 0.905 0.043\n186SOL H 0.539 0.851 0.047\n187SOL O 0.539 0.692 0.188\n187SOL H 0.539 0.692 0.204\n187SOL H 0.571 0.671 0.184\n188SOL O 0.538 0.884 0.188\n188SOL H 0.537 0.851 0.179\n188SOL H 0.506 0.905 0.183\n189SOL O 0.597 0.538 0.075\n189SOL H 0.609 0.570 0.067\n189SOL H 0.609 0.507 0.067\n190SOL O 0.598 0.538 0.152\n190SOL H 0.615 0.538 0.137\n190SOL H 0.558 0.538 0.149\n191SOL O 0.979 0.538 0.074\n191SOL H 0.962 0.538 0.089\n191SOL H 1.018 0.538 0.077\n192SOL O 0.979 0.538 0.151\n192SOL H 0.968 0.570 0.159\n192SOL H 0.968 0.507 0.159\n193SOL O 0.693 0.597 0.015\n193SOL H 0.672 0.609 0.028\n193SOL H 0.726 0.620 0.015\n194SOL O 0.884 0.597 0.015\n194SOL H 0.905 0.609 0.028\n194SOL H 0.851 0.620 0.015\n195SOL O 0.693 0.979 0.015\n195SOL H 0.672 0.968 0.028\n195SOL H 0.726 0.956 0.015\n196SOL O 0.885 0.980 0.015\n196SOL H 0.906 0.969 0.028\n196SOL H 0.885 1.021 0.015\n197SOL O 0.940 0.788 0.090\n197SOL H 0.947 0.788 0.105\n197SOL H 0.976 0.788 0.083\n198SOL O 0.635 0.788 0.090\n198SOL H 0.656 0.757 0.085\n198SOL H 0.656 0.820 0.085\n199SOL O 0.942 0.788 0.136\n199SOL H 0.920 0.757 0.141\n199SOL H 0.920 0.820 0.140\n200SOL O 0.634 0.788 0.136\n200SOL H 0.635 0.787 0.120\n200SOL H 0.656 0.756 0.140\n201SOL O 0.789 0.847 0.053\n201SOL H 0.821 0.859 0.061\n201SOL H 0.789 0.870 0.040\n202SOL O 0.789 0.847 0.173\n202SOL H 0.789 0.870 0.186\n202SOL H 0.821 0.859 0.165\n203SOL O 0.789 0.731 0.053\n203SOL H 0.789 0.771 0.053\n203SOL H 0.821 0.720 0.061\n204SOL O 0.788 0.731 0.173\n204SOL H 0.789 0.771 0.173\n204SOL H 0.821 0.720 0.165\n205SOL O 0.848 0.940 0.113\n205SOL H 0.808 0.947 0.113\n205SOL H 0.865 0.976 0.113\n206SOL O 0.847 0.634 0.113\n206SOL H 0.859 0.655 0.100\n206SOL H 0.870 0.601 0.113\n207SOL O 0.730 0.943 0.113\n207SOL H 0.718 0.921 0.126\n207SOL H 0.706 0.976 0.113\n208SOL O 0.731 0.634 0.113\n208SOL H 0.771 0.635 0.113\n208SOL H 0.720 0.656 0.126\n209SOL O 0.630 0.630 0.051\n209SOL H 0.598 0.653 0.047\n209SOL H 0.653 0.654 0.060\n210SOL O 0.947 0.947 0.051\n210SOL H 0.923 0.924 0.060\n210SOL H 0.957 0.978 0.060\n211SOL O 0.947 0.630 0.175\n211SOL H 0.926 0.619 0.188\n211SOL H 0.978 0.651 0.180\n212SOL O 0.630 0.947 0.175\n212SOL H 0.599 0.924 0.179\n212SOL H 0.652 0.923 0.166\n213SOL O 0.880 0.880 0.077\n213SOL H 0.901 0.849 0.082\n213SOL H 0.869 0.901 0.090\n214SOL O 0.696 0.696 0.077\n214SOL H 0.729 0.707 0.069\n214SOL H 0.710 0.675 0.089\n215SOL O 0.880 0.696 0.149\n215SOL H 0.869 0.674 0.137\n215SOL H 0.903 0.672 0.158\n216SOL O 0.696 0.880 0.149\n216SOL H 0.676 0.847 0.145\n216SOL H 0.728 0.867 0.157\n217SOL O 0.947 0.946 0.175\n217SOL H 0.925 0.960 0.187\n217SOL H 0.957 0.979 0.167\n218SOL O 0.630 0.630 0.175\n218SOL H 0.651 0.617 0.187\n218SOL H 0.619 0.597 0.167\n219SOL O 0.630 0.947 0.051\n219SOL H 0.653 0.924 0.060\n219SOL H 0.619 0.978 0.060\n220SOL O 0.946 0.630 0.051\n220SOL H 0.960 0.599 0.059\n220SOL H 0.979 0.650 0.047\n221SOL O 0.697 0.696 0.149\n221SOL H 0.728 0.707 0.158\n221SOL H 0.674 0.674 0.158\n222SOL O 0.880 0.880 0.149\n222SOL H 0.869 0.903 0.137\n222SOL H 0.903 0.904 0.158\n223SOL O 0.696 0.880 0.077\n223SOL H 0.729 0.869 0.069\n223SOL H 0.710 0.901 0.089\n224SOL O 0.880 0.696 0.077\n224SOL H 0.903 0.728 0.081\n224SOL H 0.904 0.674 0.068\n225SOL O 0.538 0.788 0.064\n225SOL H 0.539 0.755 0.055\n225SOL H 0.572 0.787 0.073\n226SOL O 0.538 0.788 0.162\n226SOL H 0.539 0.755 0.171\n226SOL H 0.572 0.787 0.153\n227SOL O 0.663 0.538 0.113\n227SOL H 0.641 0.539 0.100\n227SOL H 0.685 0.571 0.112\n228SOL O 0.914 0.538 0.113\n228SOL H 0.936 0.537 0.126\n228SOL H 0.891 0.505 0.113\n229SOL O 0.788 0.663 0.015\n229SOL H 0.788 0.688 0.003\n229SOL H 0.788 0.688 0.028\n230SOL O 0.788 0.913 0.015\n230SOL H 0.789 0.891 0.002\n230SOL H 0.822 0.936 0.015\n231SOL O 0.038 0.693 0.038\n231SOL H 0.039 0.692 0.022\n231SOL H 0.071 0.671 0.042\n232SOL O 0.038 0.884 0.038\n232SOL H 0.006 0.905 0.043\n232SOL H 0.037 0.851 0.047\n233SOL O 0.038 0.690 0.187\n233SOL H 0.038 0.697 0.203\n233SOL H 0.038 0.726 0.181\n234SOL O 0.038 0.884 0.188\n234SOL H 0.006 0.905 0.183\n234SOL H 0.037 0.851 0.179\n235SOL O 0.097 0.539 0.076\n235SOL H 0.109 0.571 0.067\n235SOL H 0.120 0.539 0.088\n236SOL O 0.096 0.538 0.150\n236SOL H 0.056 0.537 0.151\n236SOL H 0.107 0.506 0.159\n237SOL O 0.480 0.538 0.076\n237SOL H 0.520 0.537 0.075\n237SOL H 0.469 0.506 0.067\n238SOL O 0.480 0.539 0.150\n238SOL H 0.468 0.571 0.159\n238SOL H 0.456 0.539 0.138\n239SOL O 0.193 0.596 0.015\n239SOL H 0.171 0.607 0.002\n239SOL H 0.192 0.556 0.015\n240SOL O 0.384 0.596 0.015\n240SOL H 0.406 0.607 0.028\n240SOL H 0.384 0.556 0.015\n241SOL O 0.192 0.980 0.015\n241SOL H 0.171 0.969 0.028\n241SOL H 0.192 1.021 0.015\n242SOL O 0.384 0.979 0.015\n242SOL H 0.405 0.968 0.002\n242SOL H 0.351 0.956 0.015\n243SOL O 0.440 0.788 0.090\n243SOL H 0.476 0.788 0.083\n243SOL H 0.447 0.788 0.105\n244SOL O 0.135 0.788 0.090\n244SOL H 0.157 0.757 0.085\n244SOL H 0.157 0.820 0.085\n245SOL O 0.443 0.789 0.136\n245SOL H 0.476 0.789 0.145\n245SOL H 0.422 0.821 0.141\n246SOL O 0.134 0.788 0.136\n246SOL H 0.135 0.787 0.120\n246SOL H 0.156 0.756 0.140\n247SOL O 0.289 0.847 0.052\n247SOL H 0.320 0.859 0.061\n247SOL H 0.289 0.870 0.040\n248SOL O 0.289 0.847 0.173\n248SOL H 0.289 0.870 0.186\n248SOL H 0.321 0.859 0.165\n249SOL O 0.288 0.728 0.054\n249SOL H 0.288 0.768 0.051\n249SOL H 0.288 0.712 0.039\n250SOL O 0.288 0.731 0.173\n250SOL H 0.289 0.771 0.173\n250SOL H 0.321 0.720 0.165\n251SOL O 0.347 0.943 0.113\n251SOL H 0.307 0.942 0.113\n251SOL H 0.357 0.921 0.100\n252SOL O 0.347 0.635 0.113\n252SOL H 0.359 0.656 0.125\n252SOL H 0.359 0.656 0.100\n253SOL O 0.229 0.942 0.113\n253SOL H 0.218 0.920 0.125\n253SOL H 0.218 0.920 0.100\n254SOL O 0.230 0.634 0.113\n254SOL H 0.270 0.635 0.113\n254SOL H 0.219 0.656 0.125\n255SOL O 0.130 0.630 0.051\n255SOL H 0.153 0.619 0.039\n255SOL H 0.154 0.653 0.060\n256SOL O 0.447 0.947 0.051\n256SOL H 0.478 0.925 0.046\n256SOL H 0.426 0.957 0.038\n257SOL O 0.447 0.630 0.175\n257SOL H 0.426 0.619 0.188\n257SOL H 0.478 0.651 0.180\n258SOL O 0.130 0.947 0.175\n258SOL H 0.099 0.924 0.179\n258SOL H 0.152 0.923 0.166\n259SOL O 0.380 0.880 0.077\n259SOL H 0.403 0.849 0.081\n259SOL H 0.404 0.903 0.068\n260SOL O 0.196 0.696 0.077\n260SOL H 0.229 0.707 0.069\n260SOL H 0.210 0.675 0.089\n261SOL O 0.380 0.697 0.149\n261SOL H 0.403 0.728 0.145\n261SOL H 0.404 0.674 0.158\n262SOL O 0.197 0.880 0.149\n262SOL H 0.176 0.847 0.145\n262SOL H 0.228 0.867 0.157\n263SOL O 0.447 0.946 0.175\n263SOL H 0.425 0.960 0.187\n263SOL H 0.457 0.979 0.167\n264SOL O 0.131 0.630 0.175\n264SOL H 0.097 0.651 0.179\n264SOL H 0.117 0.599 0.167\n265SOL O 0.131 0.947 0.051\n265SOL H 0.097 0.926 0.047\n265SOL H 0.117 0.978 0.059\n266SOL O 0.447 0.630 0.051\n266SOL H 0.458 0.599 0.060\n266SOL H 0.424 0.652 0.060\n267SOL O 0.197 0.696 0.149\n267SOL H 0.228 0.707 0.158\n267SOL H 0.174 0.673 0.158\n268SOL O 0.380 0.880 0.149\n268SOL H 0.369 0.903 0.137\n268SOL H 0.403 0.904 0.158\n269SOL O 0.197 0.880 0.077\n269SOL H 0.228 0.869 0.068\n269SOL H 0.174 0.903 0.068\n270SOL O 0.380 0.696 0.077\n270SOL H 0.400 0.729 0.081\n270SOL H 0.349 0.709 0.069\n271SOL O 0.038 0.788 0.064\n271SOL H 0.040 0.755 0.055\n271SOL H 0.072 0.787 0.073\n272SOL O 0.038 0.788 0.162\n272SOL H 0.007 0.788 0.152\n272SOL H 0.070 0.788 0.152\n273SOL O 0.163 0.538 0.113\n273SOL H 0.141 0.539 0.126\n273SOL H 0.186 0.572 0.113\n274SOL O 0.413 0.538 0.113\n274SOL H 0.436 0.539 0.100\n274SOL H 0.391 0.572 0.112\n275SOL O 0.288 0.664 0.015\n275SOL H 0.256 0.639 0.015\n275SOL H 0.320 0.639 0.015\n276SOL O 0.288 0.913 0.015\n276SOL H 0.287 0.891 0.002\n276SOL H 0.255 0.936 0.015\n277SOL O 0.538 0.190 0.038\n277SOL H 0.538 0.197 0.023\n277SOL H 0.538 0.227 0.045\n278SOL O 0.538 0.387 0.038\n278SOL H 0.538 0.380 0.023\n278SOL H 0.538 0.350 0.045\n279SOL O 0.538 0.193 0.188\n279SOL H 0.537 0.226 0.179\n279SOL H 0.506 0.171 0.183\n280SOL O 0.538 0.385 0.188\n280SOL H 0.570 0.406 0.183\n280SOL H 0.507 0.406 0.183\n281SOL O 0.596 0.038 0.076\n281SOL H 0.607 0.071 0.067\n281SOL H 0.556 0.039 0.075\n282SOL O 0.597 0.038 0.150\n282SOL H 0.609 0.006 0.159\n282SOL H 0.621 0.037 0.138\n283SOL O 0.979 0.038 0.076\n283SOL H 0.968 0.071 0.067\n283SOL H 0.956 0.039 0.088\n284SOL O 0.980 0.038 0.151\n284SOL H 0.969 0.071 0.159\n284SOL H 1.021 0.039 0.151\n285SOL O 0.693 0.096 0.015\n285SOL H 0.692 0.056 0.015\n285SOL H 0.671 0.107 0.028\n286SOL O 0.885 0.097 0.015\n286SOL H 0.907 0.109 0.002\n286SOL H 0.907 0.109 0.028\n287SOL O 0.693 0.480 0.015\n287SOL H 0.692 0.520 0.015\n287SOL H 0.671 0.469 0.028\n288SOL O 0.887 0.479 0.015\n288SOL H 0.880 0.518 0.015\n288SOL H 0.851 0.462 0.015\n289SOL O 0.943 0.288 0.090\n289SOL H 0.942 0.289 0.106\n289SOL H 0.921 0.321 0.086\n290SOL O 0.634 0.288 0.090\n290SOL H 0.635 0.288 0.106\n290SOL H 0.656 0.256 0.086\n291SOL O 0.943 0.289 0.136\n291SOL H 0.921 0.321 0.141\n291SOL H 0.976 0.289 0.145\n292SOL O 0.634 0.288 0.136\n292SOL H 0.655 0.256 0.141\n292SOL H 0.601 0.288 0.145\n293SOL O 0.788 0.347 0.053\n293SOL H 0.756 0.359 0.061\n293SOL H 0.787 0.371 0.040\n294SOL O 0.788 0.348 0.172\n294SOL H 0.788 0.365 0.186\n294SOL H 0.788 0.309 0.175\n295SOL O 0.788 0.229 0.054\n295SOL H 0.788 0.268 0.051\n295SOL H 0.788 0.212 0.039\n296SOL O 0.789 0.229 0.173\n296SOL H 0.789 0.206 0.186\n296SOL H 0.821 0.218 0.165\n297SOL O 0.846 0.443 0.113\n297SOL H 0.806 0.442 0.113\n297SOL H 0.857 0.421 0.100\n298SOL O 0.847 0.135 0.113\n298SOL H 0.859 0.157 0.125\n298SOL H 0.859 0.157 0.100\n299SOL O 0.730 0.443 0.113\n299SOL H 0.707 0.476 0.113\n299SOL H 0.718 0.421 0.126\n300SOL O 0.731 0.134 0.113\n300SOL H 0.771 0.135 0.113\n300SOL H 0.720 0.156 0.125\n301SOL O 0.630 0.130 0.051\n301SOL H 0.598 0.153 0.047\n301SOL H 0.653 0.154 0.060\n302SOL O 0.947 0.446 0.051\n302SOL H 0.957 0.480 0.059\n302SOL H 0.925 0.460 0.039\n303SOL O 0.947 0.130 0.175\n303SOL H 0.926 0.120 0.188\n303SOL H 0.978 0.151 0.180\n304SOL O 0.630 0.447 0.175\n304SOL H 0.619 0.478 0.166\n304SOL H 0.653 0.424 0.165\n305SOL O 0.880 0.380 0.077\n305SOL H 0.848 0.369 0.068\n305SOL H 0.903 0.403 0.068\n306SOL O 0.696 0.197 0.077\n306SOL H 0.729 0.207 0.069\n306SOL H 0.710 0.175 0.089\n307SOL O 0.880 0.197 0.149\n307SOL H 0.903 0.228 0.145\n307SOL H 0.904 0.174 0.158\n308SOL O 0.696 0.381 0.149\n308SOL H 0.676 0.347 0.145\n308SOL H 0.728 0.367 0.157\n309SOL O 0.947 0.447 0.175\n309SOL H 0.924 0.422 0.166\n309SOL H 0.978 0.424 0.179\n310SOL O 0.630 0.130 0.175\n310SOL H 0.597 0.150 0.179\n310SOL H 0.617 0.099 0.167\n311SOL O 0.630 0.447 0.051\n311SOL H 0.598 0.424 0.047\n311SOL H 0.653 0.422 0.060\n312SOL O 0.947 0.130 0.051\n312SOL H 0.924 0.154 0.060\n312SOL H 0.978 0.153 0.047\n313SOL O 0.697 0.197 0.149\n313SOL H 0.728 0.208 0.158\n313SOL H 0.674 0.173 0.158\n314SOL O 0.880 0.380 0.149\n314SOL H 0.847 0.369 0.157\n314SOL H 0.867 0.402 0.137\n315SOL O 0.696 0.380 0.077\n315SOL H 0.675 0.349 0.082\n315SOL H 0.707 0.401 0.090\n316SOL O 0.880 0.196 0.077\n316SOL H 0.901 0.229 0.081\n316SOL H 0.849 0.210 0.069\n317SOL O 0.538 0.288 0.064\n317SOL H 0.570 0.288 0.074\n317SOL H 0.507 0.288 0.074\n318SOL O 0.538 0.288 0.162\n318SOL H 0.537 0.322 0.171\n318SOL H 0.505 0.290 0.153\n319SOL O 0.663 0.038 0.113\n319SOL H 0.641 0.040 0.100\n319SOL H 0.685 0.072 0.112\n320SOL O 0.913 0.038 0.113\n320SOL H 0.936 0.040 0.126\n320SOL H 0.891 0.072 0.113\n321SOL O 0.788 0.164 0.015\n321SOL H 0.757 0.139 0.015\n321SOL H 0.820 0.139 0.015\n322SOL O 0.788 0.413 0.015\n322SOL H 0.787 0.391 0.002\n322SOL H 0.755 0.436 0.015\n323SOL O 0.039 0.193 0.038\n323SOL H 0.071 0.171 0.043\n323SOL H 0.039 0.226 0.047\n324SOL O 0.038 0.385 0.038\n324SOL H 0.006 0.406 0.043\n324SOL H 0.070 0.406 0.043\n325SOL O 0.038 0.190 0.187\n325SOL H 0.038 0.197 0.203\n325SOL H 0.038 0.227 0.181\n326SOL O 0.038 0.384 0.188\n326SOL H 0.039 0.384 0.204\n326SOL H 0.071 0.406 0.184\n327SOL O 0.098 0.038 0.074\n327SOL H 0.059 0.038 0.077\n327SOL H 0.115 0.038 0.089\n328SOL O 0.097 0.038 0.150\n328SOL H 0.109 0.006 0.159\n328SOL H 0.120 0.037 0.138\n329SOL O 0.480 0.038 0.075\n329SOL H 0.468 0.006 0.067\n329SOL H 0.468 0.070 0.067\n330SOL O 0.479 0.038 0.152\n330SOL H 0.518 0.038 0.149\n330SOL H 0.462 0.038 0.137\n331SOL O 0.193 0.097 0.015\n331SOL H 0.171 0.109 0.002\n331SOL H 0.226 0.120 0.015\n332SOL O 0.384 0.096 0.015\n332SOL H 0.384 0.056 0.015\n332SOL H 0.406 0.107 0.028\n333SOL O 0.192 0.480 0.015\n333SOL H 0.171 0.468 0.002\n333SOL H 0.171 0.468 0.028\n334SOL O 0.384 0.480 0.015\n334SOL H 0.405 0.468 0.002\n334SOL H 0.351 0.456 0.015\n335SOL O 0.443 0.288 0.090\n335SOL H 0.442 0.288 0.106\n335SOL H 0.421 0.256 0.086\n336SOL O 0.134 0.288 0.090\n336SOL H 0.155 0.256 0.085\n336SOL H 0.101 0.288 0.081\n337SOL O 0.442 0.288 0.136\n337SOL H 0.421 0.256 0.141\n337SOL H 0.420 0.320 0.141\n338SOL O 0.134 0.288 0.136\n338SOL H 0.135 0.288 0.120\n338SOL H 0.156 0.256 0.140\n339SOL O 0.288 0.347 0.053\n339SOL H 0.320 0.359 0.061\n339SOL H 0.256 0.359 0.061\n340SOL O 0.288 0.347 0.173\n340SOL H 0.288 0.371 0.186\n340SOL H 0.256 0.359 0.165\n341SOL O 0.288 0.229 0.054\n341SOL H 0.288 0.268 0.051\n341SOL H 0.288 0.212 0.039\n342SOL O 0.288 0.229 0.172\n342SOL H 0.288 0.212 0.186\n342SOL H 0.288 0.268 0.175\n343SOL O 0.347 0.443 0.113\n343SOL H 0.371 0.476 0.113\n343SOL H 0.359 0.421 0.100\n344SOL O 0.348 0.137 0.113\n344SOL H 0.309 0.130 0.113\n344SOL H 0.365 0.101 0.113\n345SOL O 0.229 0.440 0.113\n345SOL H 0.212 0.476 0.113\n345SOL H 0.268 0.447 0.113\n346SOL O 0.229 0.135 0.113\n346SOL H 0.218 0.157 0.100\n346SOL H 0.218 0.157 0.125\n347SOL O 0.130 0.131 0.051\n347SOL H 0.119 0.097 0.059\n347SOL H 0.152 0.117 0.039\n348SOL O 0.447 0.447 0.051\n348SOL H 0.478 0.425 0.046\n348SOL H 0.426 0.458 0.038\n349SOL O 0.447 0.130 0.175\n349SOL H 0.458 0.099 0.166\n349SOL H 0.424 0.152 0.165\n350SOL O 0.130 0.447 0.175\n350SOL H 0.153 0.458 0.187\n350SOL H 0.154 0.424 0.166\n351SOL O 0.380 0.380 0.077\n351SOL H 0.403 0.349 0.081\n351SOL H 0.404 0.403 0.068\n352SOL O 0.197 0.197 0.077\n352SOL H 0.228 0.208 0.068\n352SOL H 0.174 0.173 0.068\n353SOL O 0.381 0.197 0.149\n353SOL H 0.347 0.208 0.157\n353SOL H 0.367 0.175 0.137\n354SOL O 0.197 0.380 0.149\n354SOL H 0.175 0.349 0.144\n354SOL H 0.208 0.401 0.136\n355SOL O 0.447 0.447 0.175\n355SOL H 0.458 0.478 0.166\n355SOL H 0.424 0.424 0.165\n356SOL O 0.131 0.130 0.175\n356SOL H 0.097 0.151 0.179\n356SOL H 0.117 0.099 0.167\n357SOL O 0.130 0.447 0.051\n357SOL H 0.119 0.478 0.060\n357SOL H 0.153 0.424 0.060\n358SOL O 0.447 0.130 0.051\n358SOL H 0.478 0.153 0.047\n358SOL H 0.424 0.154 0.060\n359SOL O 0.197 0.197 0.149\n359SOL H 0.228 0.208 0.158\n359SOL H 0.174 0.173 0.158\n360SOL O 0.381 0.380 0.149\n360SOL H 0.347 0.369 0.157\n360SOL H 0.367 0.402 0.137\n361SOL O 0.197 0.380 0.077\n361SOL H 0.175 0.349 0.082\n361SOL H 0.208 0.401 0.090\n362SOL O 0.381 0.197 0.077\n362SOL H 0.347 0.208 0.069\n362SOL H 0.367 0.175 0.089\n363SOL O 0.038 0.288 0.064\n363SOL H 0.005 0.290 0.073\n363SOL H 0.037 0.322 0.055\n364SOL O 0.038 0.288 0.162\n364SOL H 0.040 0.322 0.171\n364SOL H 0.072 0.290 0.153\n365SOL O 0.163 0.038 0.113\n365SOL H 0.188 0.007 0.113\n365SOL H 0.188 0.070 0.113\n366SOL O 0.413 0.038 0.113\n366SOL H 0.391 0.005 0.112\n366SOL H 0.436 0.037 0.100\n367SOL O 0.288 0.163 0.015\n367SOL H 0.290 0.186 0.002\n367SOL H 0.322 0.141 0.015\n368SOL O 0.288 0.413 0.015\n368SOL H 0.255 0.436 0.015\n368SOL H 0.287 0.391 0.028\n369MOL C 0.538 0.538 0.211\n369MOL H 0.523 0.496 0.215\n369MOL H 0.523 0.550 0.194\n369MOL H 0.523 0.570 0.223\n369MOL H 0.585 0.538 0.211\n370MOL C 0.788 0.788 0.309\n370MOL H 0.773 0.745 0.313\n370MOL H 0.773 0.800 0.292\n370MOL H 0.773 0.819 0.321\n370MOL H 0.835 0.788 0.309\n371MOL C 0.663 0.788 0.211\n371MOL H 0.648 0.745 0.215\n371MOL H 0.648 0.800 0.194\n371MOL H 0.648 0.819 0.223\n371MOL H 0.710 0.788 0.211\n372MOL C 0.913 0.788 0.211\n372MOL H 0.898 0.745 0.215\n372MOL H 0.898 0.800 0.194\n372MOL H 0.898 0.819 0.223\n372MOL H 0.960 0.788 0.211\n373MOL C 0.538 0.663 0.309\n373MOL H 0.523 0.621 0.313\n373MOL H 0.523 0.675 0.292\n373MOL H 0.523 0.695 0.321\n373MOL H 0.585 0.663 0.309\n374MOL C 0.538 0.913 0.309\n374MOL H 0.523 0.870 0.313\n374MOL H 0.523 0.925 0.292\n374MOL H 0.523 0.945 0.321\n374MOL H 0.585 0.913 0.309\n375MOL C 0.788 0.538 0.260\n375MOL H 0.773 0.496 0.264\n375MOL H 0.773 0.549 0.243\n375MOL H 0.773 0.570 0.272\n375MOL H 0.835 0.538 0.260\n376MOL C 0.788 0.538 0.358\n376MOL H 0.773 0.496 0.362\n376MOL H 0.773 0.549 0.341\n376MOL H 0.773 0.570 0.370\n376MOL H 0.835 0.538 0.358\n377MOL C 0.038 0.538 0.211\n377MOL H 0.023 0.496 0.215\n377MOL H 0.023 0.550 0.194\n377MOL H 0.023 0.570 0.223\n377MOL H 0.085 0.538 0.211\n378MOL C 0.288 0.788 0.309\n378MOL H 0.273 0.745 0.313\n378MOL H 0.273 0.800 0.292\n378MOL H 0.273 0.819 0.321\n378MOL H 0.336 0.788 0.309\n379MOL C 0.163 0.788 0.211\n379MOL H 0.148 0.745 0.215\n379MOL H 0.148 0.800 0.194\n379MOL H 0.148 0.819 0.223\n379MOL H 0.211 0.788 0.211\n380MOL C 0.413 0.788 0.211\n380MOL H 0.398 0.745 0.215\n380MOL H 0.398 0.800 0.194\n380MOL H 0.398 0.819 0.223\n380MOL H 0.461 0.788 0.211\n381MOL C 0.038 0.663 0.309\n381MOL H 0.023 0.621 0.313\n381MOL H 0.023 0.675 0.292\n381MOL H 0.023 0.695 0.321\n381MOL H 0.085 0.663 0.309\n382MOL C 0.038 0.913 0.309\n382MOL H 0.023 0.870 0.313\n382MOL H 0.023 0.925 0.292\n382MOL H 0.023 0.945 0.321\n382MOL H 0.085 0.913 0.309\n383MOL C 0.288 0.538 0.260\n383MOL H 0.273 0.496 0.264\n383MOL H 0.273 0.549 0.243\n383MOL H 0.273 0.570 0.272\n383MOL H 0.336 0.538 0.260\n384MOL C 0.288 0.538 0.358\n384MOL H 0.273 0.496 0.362\n384MOL H 0.273 0.549 0.341\n384MOL H 0.273 0.570 0.370\n384MOL H 0.336 0.538 0.358\n385MOL C 0.538 0.038 0.211\n385MOL H 0.523 -0.005 0.215\n385MOL H 0.523 0.050 0.194\n385MOL H 0.523 0.069 0.223\n385MOL H 0.585 0.038 0.211\n386MOL C 0.788 0.288 0.309\n386MOL H 0.773 0.245 0.313\n386MOL H 0.773 0.300 0.292\n386MOL H 0.773 0.320 0.321\n386MOL H 0.835 0.288 0.309\n387MOL C 0.663 0.288 0.211\n387MOL H 0.648 0.245 0.215\n387MOL H 0.648 0.300 0.194\n387MOL H 0.648 0.320 0.223\n387MOL H 0.710 0.288 0.211\n388MOL C 0.913 0.288 0.211\n388MOL H 0.898 0.245 0.215\n388MOL H 0.898 0.300 0.194\n388MOL H 0.898 0.320 0.223\n388MOL H 0.960 0.288 0.211\n389MOL C 0.538 0.163 0.309\n389MOL H 0.523 0.120 0.313\n389MOL H 0.523 0.175 0.292\n389MOL H 0.523 0.195 0.321\n389MOL H 0.585 0.163 0.309\n390MOL C 0.538 0.413 0.309\n390MOL H 0.523 0.371 0.313\n390MOL H 0.523 0.425 0.292\n390MOL H 0.523 0.445 0.321\n390MOL H 0.585 0.413 0.309\n391MOL C 0.788 0.038 0.260\n391MOL H 0.773 -0.005 0.264\n391MOL H 0.773 0.050 0.243\n391MOL H 0.773 0.070 0.272\n391MOL H 0.835 0.038 0.260\n392MOL C 0.788 0.038 0.358\n392MOL H 0.773 -0.005 0.362\n392MOL H 0.773 0.050 0.341\n392MOL H 0.773 0.070 0.370\n392MOL H 0.835 0.038 0.358\n393MOL C 0.038 0.038 0.211\n393MOL H 0.023 -0.005 0.215\n393MOL H 0.023 0.050 0.194\n393MOL H 0.023 0.069 0.223\n393MOL H 0.085 0.038 0.211\n394MOL C 0.288 0.288 0.309\n394MOL H 0.273 0.245 0.313\n394MOL H 0.273 0.300 0.292\n394MOL H 0.273 0.320 0.321\n394MOL H 0.336 0.288 0.309\n395MOL C 0.163 0.288 0.211\n395MOL H 0.148 0.245 0.215\n395MOL H 0.148 0.300 0.194\n395MOL H 0.148 0.320 0.223\n395MOL H 0.211 0.288 0.211\n396MOL C 0.413 0.288 0.211\n396MOL H 0.398 0.245 0.215\n396MOL H 0.398 0.300 0.194\n396MOL H 0.398 0.320 0.223\n396MOL H 0.461 0.288 0.211\n397MOL C 0.038 0.163 0.309\n397MOL H 0.023 0.120 0.313\n397MOL H 0.023 0.175 0.292\n397MOL H 0.023 0.195 0.321\n397MOL H 0.085 0.163 0.309\n398MOL C 0.038 0.413 0.309\n398MOL H 0.023 0.371 0.313\n398MOL H 0.023 0.425 0.292\n398MOL H 0.023 0.445 0.321\n398MOL H 0.085 0.413 0.309\n399MOL C 0.288 0.038 0.260\n399MOL H 0.273 -0.005 0.264\n399MOL H 0.273 0.050 0.243\n399MOL H 0.273 0.070 0.272\n399MOL H 0.336 0.038 0.260\n400MOL C 0.288 0.038 0.358\n400MOL H 0.273 -0.005 0.362\n400MOL H 0.273 0.050 0.341\n400MOL H 0.273 0.070 0.370\n400MOL H 0.336 0.038 0.358\n401MOL C 0.538 0.538 0.015\n401MOL H 0.523 0.496 0.019\n401MOL H 0.523 0.549 -0.002\n401MOL H 0.523 0.570 0.027\n401MOL H 0.585 0.538 0.015\n402MOL C 0.788 0.788 0.113\n402MOL H 0.773 0.745 0.117\n402MOL H 0.773 0.800 0.096\n402MOL H 0.773 0.819 0.125\n402MOL H 0.835 0.788 0.113\n403MOL C 0.663 0.788 0.015\n403MOL H 0.648 0.745 0.019\n403MOL H 0.648 0.800 -0.002\n403MOL H 0.648 0.820 0.027\n403MOL H 0.710 0.788 0.015\n404MOL C 0.913 0.788 0.015\n404MOL H 0.898 0.745 0.019\n404MOL H 0.898 0.800 -0.002\n404MOL H 0.898 0.820 0.027\n404MOL H 0.960 0.788 0.015\n405MOL C 0.538 0.663 0.113\n405MOL H 0.523 0.621 0.117\n405MOL H 0.523 0.675 0.096\n405MOL H 0.523 0.695 0.125\n405MOL H 0.585 0.663 0.113\n406MOL C 0.538 0.913 0.113\n406MOL H 0.523 0.870 0.117\n406MOL H 0.523 0.925 0.096\n406MOL H 0.523 0.945 0.125\n406MOL H 0.585 0.913 0.113\n407MOL C 0.788 0.538 0.064\n407MOL H 0.773 0.496 0.069\n407MOL H 0.773 0.549 0.047\n407MOL H 0.773 0.570 0.076\n407MOL H 0.835 0.538 0.064\n408MOL C 0.788 0.538 0.162\n408MOL H 0.773 0.496 0.167\n408MOL H 0.773 0.549 0.145\n408MOL H 0.773 0.570 0.174\n408MOL H 0.835 0.538 0.162\n409MOL C 0.038 0.538 0.015\n409MOL H 0.023 0.496 0.019\n409MOL H 0.023 0.549 -0.002\n409MOL H 0.023 0.570 0.027\n409MOL H 0.085 0.538 0.015\n410MOL C 0.288 0.788 0.113\n410MOL H 0.273 0.745 0.117\n410MOL H 0.273 0.800 0.096\n410MOL H 0.273 0.819 0.125\n410MOL H 0.336 0.788 0.113\n411MOL C 0.163 0.788 0.015\n411MOL H 0.148 0.745 0.019\n411MOL H 0.148 0.800 -0.002\n411MOL H 0.148 0.820 0.027\n411MOL H 0.211 0.788 0.015\n412MOL C 0.413 0.788 0.015\n412MOL H 0.398 0.745 0.019\n412MOL H 0.398 0.800 -0.002\n412MOL H 0.398 0.820 0.027\n412MOL H 0.461 0.788 0.015\n413MOL C 0.038 0.663 0.113\n413MOL H 0.023 0.621 0.117\n413MOL H 0.023 0.675 0.096\n413MOL H 0.023 0.695 0.125\n413MOL H 0.085 0.663 0.113\n414MOL C 0.038 0.913 0.113\n414MOL H 0.023 0.870 0.117\n414MOL H 0.023 0.925 0.096\n414MOL H 0.023 0.945 0.125\n414MOL H 0.085 0.913 0.113\n415MOL C 0.288 0.538 0.064\n415MOL H 0.273 0.496 0.069\n415MOL H 0.273 0.549 0.047\n415MOL H 0.273 0.570 0.076\n415MOL H 0.336 0.538 0.064\n416MOL C 0.288 0.538 0.162\n416MOL H 0.273 0.496 0.167\n416MOL H 0.273 0.549 0.145\n416MOL H 0.273 0.570 0.174\n416MOL H 0.336 0.538 0.162\n417MOL C 0.538 0.038 0.015\n417MOL H 0.523 -0.005 0.019\n417MOL H 0.523 0.050 -0.002\n417MOL H 0.523 0.069 0.027\n417MOL H 0.585 0.038 0.015\n418MOL C 0.788 0.288 0.113\n418MOL H 0.773 0.245 0.117\n418MOL H 0.773 0.300 0.096\n418MOL H 0.773 0.320 0.125\n418MOL H 0.835 0.288 0.113\n419MOL C 0.663 0.288 0.015\n419MOL H 0.648 0.245 0.019\n419MOL H 0.648 0.300 -0.002\n419MOL H 0.648 0.320 0.027\n419MOL H 0.710 0.288 0.015\n420MOL C 0.913 0.288 0.015\n420MOL H 0.898 0.245 0.019\n420MOL H 0.898 0.300 -0.002\n420MOL H 0.898 0.320 0.027\n420MOL H 0.960 0.288 0.015\n421MOL C 0.538 0.163 0.113\n421MOL H 0.523 0.120 0.117\n421MOL H 0.523 0.175 0.096\n421MOL H 0.523 0.195 0.125\n421MOL H 0.585 0.163 0.113\n422MOL C 0.538 0.413 0.113\n422MOL H 0.523 0.371 0.117\n422MOL H 0.523 0.425 0.096\n422MOL H 0.523 0.445 0.125\n422MOL H 0.585 0.413 0.113\n423MOL C 0.788 0.038 0.064\n423MOL H 0.773 -0.005 0.068\n423MOL H 0.773 0.050 0.047\n423MOL H 0.773 0.070 0.076\n423MOL H 0.835 0.038 0.064\n424MOL C 0.788 0.038 0.162\n424MOL H 0.773 -0.005 0.166\n424MOL H 0.773 0.050 0.145\n424MOL H 0.773 0.070 0.174\n424MOL H 0.835 0.038 0.162\n425MOL C 0.038 0.038 0.015\n425MOL H 0.023 -0.005 0.019\n425MOL H 0.023 0.050 -0.002\n425MOL H 0.023 0.069 0.027\n425MOL H 0.085 0.038 0.015\n426MOL C 0.288 0.288 0.113\n426MOL H 0.273 0.245 0.117\n426MOL H 0.273 0.300 0.096\n426MOL H 0.273 0.320 0.125\n426MOL H 0.336 0.288 0.113\n427MOL C 0.163 0.288 0.015\n427MOL H 0.148 0.245 0.019\n427MOL H 0.148 0.300 -0.002\n427MOL H 0.148 0.320 0.027\n427MOL H 0.211 0.288 0.015\n428MOL C 0.413 0.288 0.015\n428MOL H 0.398 0.245 0.019\n428MOL H 0.398 0.300 -0.002\n428MOL H 0.398 0.320 0.027\n428MOL H 0.461 0.288 0.015\n429MOL C 0.038 0.163 0.113\n429MOL H 0.023 0.120 0.117\n429MOL H 0.023 0.175 0.096\n429MOL H 0.023 0.195 0.125\n429MOL H 0.085 0.163 0.113\n430MOL C 0.038 0.413 0.113\n430MOL H 0.023 0.371 0.117\n430MOL H 0.023 0.425 0.096\n430MOL H 0.023 0.445 0.125\n430MOL H 0.085 0.413 0.113\n431MOL C 0.288 0.038 0.064\n431MOL H 0.273 -0.005 0.068\n431MOL H 0.273 0.050 0.047\n431MOL H 0.273 0.070 0.076\n431MOL H 0.336 0.038 0.064\n432MOL C 0.288 0.038 0.162\n432MOL H 0.273 -0.005 0.166\n432MOL H 0.273 0.050 0.145\n432MOL H 0.273 0.070 0.174\n432MOL H 0.336 0.038 0.162\n1SOL O 0.456 0.261 0.860\n1SOL H 0.490 0.263 0.851\n1SOL H 0.432 0.289 0.853\n2SOL O -0.000 0.396 0.959\n2SOL H 0.005 0.369 0.947\n2SOL H 0.010 0.377 0.972\n3SOL O 0.462 0.027 0.846\n3SOL H 0.432 0.054 0.845\n3SOL H 0.453 -0.001 0.835\n4SOL O 0.647 0.567 0.864\n4SOL H 0.664 0.534 0.858\n4SOL H 0.621 0.552 0.875\n5SOL O 0.293 0.338 0.874\n5SOL H 0.289 0.305 0.883\n5SOL H 0.333 0.339 0.871\n6SOL O 0.376 0.995 0.638\n6SOL H 0.362 0.990 0.653\n6SOL H 0.369 0.960 0.631\n7SOL O 0.567 0.206 0.600\n7SOL H 0.541 0.217 0.611\n7SOL H 0.545 0.210 0.587\n8SOL O 0.638 0.979 0.435\n8SOL H 0.633 1.000 0.422\n8SOL H 0.605 0.956 0.436\n9SOL O 0.813 0.615 0.730\n9SOL H 0.810 0.576 0.733\n9SOL H 0.784 0.632 0.738\n10SOL O 0.211 0.517 0.752\n10SOL H 0.216 0.534 0.766\n10SOL H 0.246 0.504 0.748\n11SOL O 0.829 0.328 0.456\n11SOL H 0.840 0.315 0.442\n11SOL H 0.811 0.296 0.463\n12SOL O 0.373 0.825 0.467\n12SOL H 0.341 0.820 0.458\n12SOL H 0.397 0.794 0.464\n13SOL O 0.673 0.795 0.577\n13SOL H 0.678 0.818 0.564\n13SOL H 0.638 0.807 0.583\n14SOL O 0.728 0.714 0.967\n14SOL H 0.720 0.746 0.976\n14SOL H 0.759 0.725 0.958\n15SOL O 0.943 0.637 0.985\n15SOL H 0.978 0.619 0.988\n15SOL H 0.944 0.649 0.971\n16SOL O 0.762 0.031 0.957\n16SOL H 0.760 0.000 0.967\n16SOL H 0.800 0.032 0.952\n17SOL O 0.307 0.675 0.554\n17SOL H 0.344 0.683 0.560\n17SOL H 0.300 0.637 0.558\n18SOL O 0.949 0.187 0.913\n18SOL H 0.945 0.149 0.908\n18SOL H 0.920 0.193 0.923\n19SOL O 0.079 0.636 0.953\n19SOL H 0.054 0.624 0.941\n19SOL H 0.091 0.672 0.948\n20SOL O 0.322 0.880 0.568\n20SOL H 0.317 0.848 0.577\n20SOL H 0.307 0.911 0.576\n21SOL O 0.367 0.886 0.609\n21SOL H 0.398 0.867 0.615\n21SOL H 0.382 0.907 0.597\n22SOL O 0.344 0.301 0.834\n22SOL H 0.338 0.299 0.819\n22SOL H 0.365 0.336 0.836\n23SOL O 0.208 0.152 0.962\n23SOL H 0.224 0.181 0.953\n23SOL H 0.227 0.119 0.957\n24SOL O 0.340 0.019 0.903\n24SOL H 0.312 -0.000 0.895\n24SOL H 0.345 0.054 0.895\n25SOL O 0.648 0.346 0.511\n25SOL H 0.639 0.373 0.523\n25SOL H 0.621 0.317 0.513\n26SOL O 0.587 0.319 0.458\n26SOL H 0.613 0.339 0.449\n26SOL H 0.578 0.285 0.450\n27SOL O 0.917 0.880 0.690\n27SOL H 0.951 0.867 0.683\n27SOL H 0.930 0.906 0.701\n28SOL O 0.820 0.744 0.816\n28SOL H 0.807 0.706 0.814\n28SOL H 0.804 0.756 0.829\n29SOL O 0.102 0.115 0.935\n29SOL H 0.089 0.124 0.921\n29SOL H 0.107 0.150 0.943\n30SOL O 0.474 0.024 0.578\n30SOL H 0.453 0.044 0.589\n30SOL H 0.509 0.044 0.577\n31SOL O 0.697 0.198 0.902\n31SOL H 0.717 0.229 0.908\n31SOL H 0.667 0.190 0.912\n32SOL O 0.243 0.190 0.695\n32SOL H 0.214 0.211 0.703\n32SOL H 0.239 0.153 0.700\n33SOL O 0.723 0.952 0.446\n33SOL H 0.737 0.945 0.432\n33SOL H 0.725 0.917 0.454\n34SOL O 0.917 0.370 0.779\n34SOL H 0.906 0.350 0.792\n34SOL H 0.913 0.343 0.768\n35SOL O 0.785 0.270 0.749\n35SOL H 0.766 0.296 0.739\n35SOL H 0.824 0.273 0.745\n36SOL O 0.870 0.764 0.761\n36SOL H 0.859 0.792 0.771\n36SOL H 0.907 0.774 0.756\n37SOL O 0.590 0.598 0.425\n37SOL H 0.608 0.584 0.412\n37SOL H 0.577 0.635 0.421\n38SOL O 0.434 0.235 0.783\n38SOL H 0.430 0.272 0.777\n38SOL H 0.399 0.227 0.790\n39SOL O 0.099 0.349 0.658\n39SOL H 0.077 0.357 0.645\n39SOL H 0.131 0.371 0.657\n40SOL O 0.156 0.932 0.584\n40SOL H 0.192 0.949 0.581\n40SOL H 0.144 0.952 0.597\n41SOL O 0.878 0.836 0.923\n41SOL H 0.879 0.868 0.933\n41SOL H 0.910 0.840 0.914\n42SOL O 0.632 0.095 0.765\n42SOL H 0.598 0.073 0.763\n42SOL H 0.657 0.070 0.772\n43SOL O 0.466 0.257 0.654\n43SOL H 0.470 0.297 0.653\n43SOL H 0.435 0.248 0.645\n44SOL O 0.521 0.133 0.681\n44SOL H 0.520 0.165 0.691\n44SOL H 0.495 0.107 0.688\n45SOL O 0.186 0.977 0.937\n45SOL H 0.204 0.955 0.948\n45SOL H 0.190 0.957 0.923\n46SOL O 0.415 0.609 0.832\n46SOL H 0.437 0.614 0.845\n46SOL H 0.425 0.573 0.826\n47SOL O 0.119 0.068 0.662\n47SOL H 0.124 0.036 0.653\n47SOL H 0.090 0.090 0.656\n48SOL O 0.710 0.535 0.675\n48SOL H 0.743 0.541 0.666\n48SOL H 0.699 0.498 0.671\n49SOL O 0.308 0.258 0.461\n49SOL H 0.321 0.278 0.474\n49SOL H 0.269 0.256 0.463\n50SOL O 0.943 0.575 0.922\n50SOL H 0.919 0.607 0.922\n50SOL H 0.918 0.544 0.919\n51SOL O 0.764 0.656 0.605\n51SOL H 0.759 0.653 0.590\n51SOL H 0.803 0.648 0.608\n52SOL O 0.385 0.094 0.924\n52SOL H 0.389 0.128 0.916\n52SOL H 0.422 0.080 0.925\n53SOL O 0.684 0.979 0.502\n53SOL H 0.650 1.000 0.504\n53SOL H 0.712 1.000 0.509\n54SOL O 0.468 0.900 0.452\n54SOL H 0.462 0.928 0.463\n54SOL H 0.501 0.912 0.444\n55SOL O 0.673 0.182 0.431\n55SOL H 0.647 0.155 0.436\n55SOL H 0.705 0.178 0.440\n56SOL O 0.491 0.301 0.602\n56SOL H 0.461 0.315 0.594\n56SOL H 0.474 0.281 0.614\n57SOL O 0.499 0.910 0.815\n57SOL H 0.463 0.922 0.820\n57SOL H 0.519 0.944 0.812\n58SOL O 0.147 0.960 0.512\n58SOL H 0.121 0.985 0.505\n58SOL H 0.180 0.981 0.514\n59SOL O 0.770 0.888 0.618\n59SOL H 0.786 0.864 0.607\n59SOL H 0.802 0.900 0.627\n60SOL O 0.524 0.040 0.888\n60SOL H 0.521 0.046 0.873\n60SOL H 0.523 0.000 0.890\n61SOL O 0.747 0.734 0.922\n61SOL H 0.725 0.749 0.910\n61SOL H 0.783 0.752 0.921\n62SOL O 0.435 0.134 0.650\n62SOL H 0.445 0.162 0.661\n62SOL H 0.399 0.121 0.655\n63SOL O 0.032 0.259 0.727\n63SOL H 0.001 0.241 0.734\n63SOL H 0.060 0.229 0.726\n64SOL O 0.629 0.278 0.748\n64SOL H 0.664 0.262 0.744\n64SOL H 0.633 0.291 0.762\n65SOL O 0.679 0.213 0.507\n65SOL H 0.686 0.245 0.517\n65SOL H 0.659 0.229 0.495\n66SOL O 0.448 0.900 0.650\n66SOL H 0.465 0.871 0.659\n66SOL H 0.464 0.934 0.655\n67SOL O 0.395 0.013 0.795\n67SOL H 0.431 0.000 0.800\n67SOL H 0.389 0.047 0.803\n68SOL O 0.887 0.459 0.791\n68SOL H 0.901 0.441 0.805\n68SOL H 0.920 0.462 0.782\n69SOL O 0.133 0.693 0.489\n69SOL H 0.173 0.684 0.488\n69SOL H 0.127 0.720 0.478\n70SOL O 0.699 0.517 0.591\n70SOL H 0.728 0.494 0.596\n70SOL H 0.667 0.507 0.600\n71SOL O 0.327 0.260 0.971\n71SOL H 0.320 0.267 0.987\n71SOL H 0.317 0.222 0.970\n72SOL O 0.488 0.937 0.746\n72SOL H 0.472 0.912 0.736\n72SOL H 0.486 0.973 0.739\n73SOL O 0.919 0.700 0.613\n73SOL H 0.929 0.669 0.622\n73SOL H 0.879 0.701 0.614\n74SOL O 0.723 0.877 0.992\n74SOL H 0.748 0.853 1.000\n74SOL H 0.688 0.858 0.993\n75SOL O 0.333 0.279 0.709\n75SOL H 0.326 0.240 0.708\n75SOL H 0.299 0.296 0.704\n76SOL O 0.990 0.802 0.935\n76SOL H 0.971 0.770 0.940\n76SOL H 0.975 0.832 0.944\n77SOL O 0.308 0.302 0.771\n77SOL H 0.306 0.277 0.783\n77SOL H 0.282 0.287 0.760\n78SOL O 0.149 0.380 0.902\n78SOL H 0.154 0.341 0.899\n78SOL H 0.136 0.397 0.889\n79SOL O 0.437 0.850 0.964\n79SOL H 0.450 0.856 0.949\n79SOL H 0.397 0.849 0.963\n80SOL O 0.976 0.099 0.940\n80SOL H 0.938 0.091 0.936\n80SOL H 0.994 0.063 0.941\n81SOL O 0.302 0.853 0.646\n81SOL H 0.313 0.839 0.632\n81SOL H 0.329 0.882 0.649\n82SOL O 0.151 0.268 0.540\n82SOL H 0.155 0.290 0.527\n82SOL H 0.147 0.296 0.551\n83SOL O 0.048 0.956 0.944\n83SOL H 0.054 0.962 0.959\n83SOL H 0.019 0.928 0.943\n84SOL O 0.327 0.115 0.962\n84SOL H 0.346 0.139 0.952\n84SOL H 0.339 0.078 0.959\n85SOL O 0.625 0.720 0.700\n85SOL H 0.623 0.683 0.694\n85SOL H 0.593 0.739 0.695\n86SOL O 0.750 0.978 0.740\n86SOL H 0.729 1.000 0.729\n86SOL H 0.788 0.989 0.738\n87SOL O 0.512 0.015 0.755\n87SOL H 0.479 0.000 0.761\n87SOL H 0.542 -0.000 0.763\n88SOL O 0.803 0.464 0.642\n88SOL H 0.788 0.438 0.652\n88SOL H 0.843 0.464 0.645\n89SOL O 0.394 0.429 0.862\n89SOL H 0.355 0.423 0.858\n89SOL H 0.403 0.401 0.872\n90SOL O 0.464 0.330 0.739\n90SOL H 0.490 0.343 0.750\n90SOL H 0.428 0.342 0.744\n91SOL O 0.598 0.607 0.803\n91SOL H 0.586 0.616 0.789\n91SOL H 0.572 0.579 0.808\n92SOL O 0.479 0.830 0.412\n92SOL H 0.444 0.835 0.420\n92SOL H 0.486 0.791 0.413\n93SOL O 0.147 0.394 0.853\n93SOL H 0.149 0.355 0.857\n93SOL H 0.149 0.395 0.837\n94SOL O 0.648 0.924 0.938\n94SOL H 0.681 0.915 0.930\n94SOL H 0.645 0.893 0.948\n95SOL O 0.232 0.270 0.858\n95SOL H 0.195 0.269 0.852\n95SOL H 0.253 0.291 0.847\n96SOL O 0.591 0.672 0.737\n96SOL H 0.629 0.669 0.732\n96SOL H 0.592 0.656 0.752\n97SOL O 0.835 0.496 0.422\n97SOL H 0.856 0.519 0.412\n97SOL H 0.852 0.459 0.420\n98SOL O 0.299 0.620 0.693\n98SOL H 0.305 0.627 0.678\n98SOL H 0.272 0.591 0.693\n99SOL O 0.069 0.069 0.996\n99SOL H 0.076 0.108 1.000\n99SOL H 0.102 0.058 0.989\n100SOL O 0.458 0.236 0.469\n100SOL H 0.456 0.218 0.455\n100SOL H 0.420 0.249 0.471\n101SOL O 0.395 0.341 0.443\n101SOL H 0.416 0.332 0.456\n101SOL H 0.360 0.323 0.445\n102SOL O 0.665 0.629 0.507\n102SOL H 0.635 0.651 0.513\n102SOL H 0.667 0.596 0.515\n103SOL O 0.940 0.028 0.454\n103SOL H 0.964 0.009 0.464\n103SOL H 0.929 0.000 0.443\n104SOL O 0.005 0.508 0.559\n104SOL H 0.023 0.491 0.572\n104SOL H 0.030 0.539 0.555\n105SOL O 0.507 0.028 0.667\n105SOL H 0.510 0.000 0.678\n105SOL H 0.468 0.035 0.665\n106SOL O 0.763 0.635 0.905\n106SOL H 0.800 0.640 0.900\n106SOL H 0.741 0.661 0.897\n107SOL O 1.000 0.151 0.449\n107SOL H 0.985 0.140 0.463\n107SOL H 1.000 0.117 0.441\n108SOL O 0.206 0.588 0.800\n108SOL H 0.196 0.617 0.811\n108SOL H 0.174 0.563 0.800\n109SOL O 0.251 0.835 0.860\n109SOL H 0.247 0.864 0.849\n109SOL H 0.226 0.846 0.872\n110SOL O 0.977 0.323 0.647\n110SOL H 1.000 0.317 0.660\n110SOL H 0.949 0.349 0.652\n111SOL O 0.254 0.429 0.896\n111SOL H 0.242 0.414 0.882\n111SOL H 0.239 0.403 0.907\n112SOL O 0.373 0.491 0.931\n112SOL H 0.371 0.483 0.916\n112SOL H 0.376 0.531 0.932\n113SOL O 0.359 0.776 0.766\n113SOL H 0.347 0.745 0.775\n113SOL H 0.340 0.808 0.772\n114SOL O 0.733 0.377 0.765\n114SOL H 0.710 0.349 0.772\n114SOL H 0.712 0.411 0.767\n115SOL O 0.869 0.976 0.954\n115SOL H 0.849 0.971 0.967\n115SOL H 0.900 1.000 0.957\n116SOL O 0.680 0.064 0.672\n116SOL H 0.674 0.067 0.657\n116SOL H 0.702 0.096 0.676\n117SOL O 0.818 0.760 0.466\n117SOL H 0.809 0.721 0.469\n117SOL H 0.842 0.759 0.454\n118SOL O 0.391 0.973 0.874\n118SOL H 0.413 1.000 0.866\n118SOL H 0.370 0.996 0.884\n119SOL O 0.342 0.035 0.844\n119SOL H 0.313 0.063 0.842\n119SOL H 0.348 0.035 0.859\n120SOL O 0.365 0.796 0.551\n120SOL H 0.386 0.830 0.550\n120SOL H 0.392 0.769 0.557\n121SOL O 0.910 0.227 0.451\n121SOL H 0.909 0.192 0.459\n121SOL H 0.875 0.229 0.444\n122SOL O 0.951 0.527 0.490\n122SOL H 0.979 0.531 0.501\n122SOL H 0.928 0.560 0.491\n123SOL O 0.632 0.827 0.813\n123SOL H 0.608 0.844 0.823\n123SOL H 0.639 0.790 0.818\n124SOL O 0.995 0.091 0.701\n124SOL H 0.995 0.058 0.710\n124SOL H 1.000 0.121 0.711\n125SOL O 0.539 0.208 0.759\n125SOL H 0.565 0.200 0.770\n125SOL H 0.528 0.246 0.761\n126SOL O 0.334 0.160 0.676\n126SOL H 0.334 0.144 0.691\n126SOL H 0.363 0.187 0.676\n127SOL O 0.448 0.316 0.965\n127SOL H 0.425 0.311 0.978\n127SOL H 0.472 0.347 0.969\n128SOL O 0.585 0.583 0.586\n128SOL H 0.603 0.553 0.577\n128SOL H 0.550 0.588 0.579\n129SOL O 0.229 0.053 0.766\n129SOL H 0.211 0.022 0.773\n129SOL H 0.267 0.051 0.770\n130SOL O 0.512 0.219 0.720\n130SOL H 0.551 0.209 0.722\n130SOL H 0.509 0.256 0.726\n131SOL O 0.387 0.909 0.924\n131SOL H 0.398 0.884 0.913\n131SOL H 0.351 0.922 0.920\n132SOL O 0.859 0.534 0.875\n132SOL H 0.823 0.520 0.879\n132SOL H 0.857 0.573 0.878\n133SOL O 0.403 0.533 0.714\n133SOL H 0.437 0.549 0.708\n133SOL H 0.381 0.565 0.718\n134SOL O 0.565 0.896 0.733\n134SOL H 0.605 0.893 0.732\n134SOL H 0.551 0.862 0.727\n135SOL O 0.761 0.155 0.656\n135SOL H 0.800 0.148 0.658\n135SOL H 0.757 0.178 0.644\n136SOL O 0.317 0.299 0.939\n136SOL H 0.293 0.266 0.937\n136SOL H 0.296 0.328 0.932\n137SOL O 0.466 0.653 0.562\n137SOL H 0.451 0.624 0.571\n137SOL H 0.467 0.685 0.571\n138SOL O 0.024 0.037 0.920\n138SOL H 0.053 0.036 0.931\n138SOL H 0.000 0.068 0.925\n139SOL O 0.085 0.289 0.945\n139SOL H 0.053 0.292 0.955\n139SOL H 0.114 0.273 0.954\n140SOL O 0.715 0.899 0.578\n140SOL H 0.686 0.912 0.568\n140SOL H 0.696 0.881 0.590\n141SOL O 0.770 0.299 0.508\n141SOL H 0.740 0.301 0.498\n141SOL H 0.795 0.330 0.504\n142SOL O 0.000 0.924 0.577\n142SOL H 0.011 0.886 0.575\n142SOL H 0.034 0.943 0.581\n143SOL O 0.827 0.103 0.461\n143SOL H 0.858 0.102 0.451\n143SOL H 0.809 0.067 0.459\n144SOL O 0.956 0.837 0.987\n144SOL H 0.936 0.816 0.976\n144SOL H 0.955 0.815 1.000\n145SOL O 0.784 0.512 0.940\n145SOL H 0.772 0.503 0.955\n145SOL H 0.819 0.492 0.939\n146SOL O 0.577 0.692 0.860\n146SOL H 0.602 0.674 0.870\n146SOL H 0.559 0.720 0.868\n147SOL O 0.236 0.680 0.514\n147SOL H 0.260 0.690 0.502\n147SOL H 0.246 0.707 0.526\n148SOL O 0.761 0.964 0.598\n148SOL H 0.797 0.972 0.592\n148SOL H 0.743 1.000 0.599\n149SOL O 0.933 0.414 0.590\n149SOL H 0.961 0.386 0.593\n149SOL H 0.944 0.446 0.599\n150SOL O 0.051 0.614 0.670\n150SOL H 0.038 0.613 0.685\n150SOL H 0.075 0.581 0.669\n151SOL O 0.383 0.935 0.760\n151SOL H 0.362 0.966 0.754\n151SOL H 0.356 0.914 0.768\n152SOL O 0.509 0.468 0.693\n152SOL H 0.480 0.447 0.686\n152SOL H 0.507 0.504 0.687\n153SOL O 0.669 0.360 0.420\n153SOL H 0.643 0.383 0.412\n153SOL H 0.704 0.365 0.412\n154SOL O 0.057 0.638 0.729\n154SOL H 0.075 0.616 0.719\n154SOL H 0.039 0.612 0.739\n155SOL O 0.029 0.462 0.642\n155SOL H 0.034 0.498 0.635\n155SOL H 0.067 0.452 0.646\n156SOL O 0.683 0.477 0.822\n156SOL H 0.703 0.443 0.820\n156SOL H 0.680 0.493 0.808\n157SOL O 0.172 0.854 0.915\n157SOL H 0.205 0.846 0.907\n157SOL H 0.146 0.826 0.911\n158SOL O 0.698 0.448 0.567\n158SOL H 0.707 0.449 0.552\n158SOL H 0.677 0.414 0.569\n159SOL O 0.765 0.571 0.622\n159SOL H 0.796 0.547 0.625\n159SOL H 0.732 0.549 0.625\n160SOL O 0.339 0.976 0.522\n160SOL H 0.373 0.997 0.519\n160SOL H 0.317 1.000 0.531\n161SOL O 0.329 0.147 0.472\n161SOL H 0.362 0.158 0.464\n161SOL H 0.337 0.160 0.487\n162SOL O 0.598 0.520 0.947\n162SOL H 0.630 0.513 0.938\n162SOL H 0.588 0.484 0.953\n163SOL O 0.096 0.809 0.693\n163SOL H 0.082 0.772 0.691\n163SOL H 0.065 0.833 0.690\n164SOL O 0.501 0.366 0.850\n164SOL H 0.524 0.365 0.863\n164SOL H 0.527 0.366 0.838\n165SOL O 0.986 0.469 0.528\n165SOL H 1.000 0.436 0.535\n165SOL H 1.000 0.499 0.537\n166SOL O 0.520 0.522 0.830\n166SOL H 0.499 0.489 0.832\n166SOL H 0.557 0.512 0.835\n167SOL O 0.661 0.172 0.570\n167SOL H 0.623 0.169 0.575\n167SOL H 0.657 0.175 0.554\n168SOL O 0.141 0.544 0.714\n168SOL H 0.172 0.520 0.718\n168SOL H 0.125 0.526 0.701\n169SOL O 0.697 0.749 0.761\n169SOL H 0.691 0.748 0.745\n169SOL H 0.672 0.778 0.766\n170SOL O 0.530 0.226 0.424\n170SOL H 0.530 0.253 0.412\n170SOL H 0.555 0.197 0.419\n171SOL O 0.865 0.877 0.883\n171SOL H 0.872 0.913 0.890\n171SOL H 0.837 0.861 0.892\n172SOL O 0.790 0.384 0.696\n172SOL H 0.813 0.368 0.707\n172SOL H 0.776 0.417 0.703\n173SOL O 0.574 0.664 0.556\n173SOL H 0.559 0.702 0.556\n173SOL H 0.577 0.655 0.541\n174SOL O 0.292 0.194 0.554\n174SOL H 0.321 0.202 0.564\n174SOL H 0.263 0.221 0.557\n175SOL O 0.547 0.604 0.859\n175SOL H 0.543 0.606 0.843\n175SOL H 0.532 0.567 0.862\n176SOL O 0.668 0.299 0.559\n176SOL H 0.689 0.321 0.569\n176SOL H 0.640 0.324 0.554\n177SOL O 0.615 0.172 0.701\n177SOL H 0.622 0.133 0.705\n177SOL H 0.648 0.191 0.706\n178SOL O 0.349 0.570 0.593\n178SOL H 0.344 0.540 0.583\n178SOL H 0.371 0.554 0.605\n179SOL O 0.474 0.707 0.532\n179SOL H 0.471 0.741 0.541\n179SOL H 0.510 0.709 0.526\n180SOL O 0.639 0.718 0.906\n180SOL H 0.614 0.686 0.906\n180SOL H 0.624 0.741 0.894\n181SOL O 0.625 0.740 0.502\n181SOL H 0.627 0.755 0.517\n181SOL H 0.587 0.745 0.498\n182SOL O 0.074 0.448 0.833\n182SOL H 0.043 0.468 0.826\n182SOL H 0.086 0.474 0.844\n183SOL O 0.520 0.578 0.543\n183SOL H 0.501 0.605 0.534\n183SOL H 0.552 0.568 0.535\n184SOL O 0.503 0.408 0.599\n184SOL H 0.480 0.425 0.588\n184SOL H 0.528 0.383 0.591\n185SOL O 0.309 0.653 0.997\n185SOL H 0.270 0.649 1.000\n185SOL H 0.312 0.685 0.988\n186SOL O 0.819 0.628 0.667\n186SOL H 0.853 0.615 0.673\n186SOL H 0.815 0.606 0.654\n187SOL O 0.147 0.226 0.921\n187SOL H 0.152 0.190 0.914\n187SOL H 0.150 0.253 0.910\n188SOL O 0.939 0.851 0.520\n188SOL H 0.917 0.883 0.522\n188SOL H 0.949 0.852 0.504\n189SOL O 0.072 0.031 0.747\n189SOL H 0.043 0.029 0.737\n189SOL H 0.091 0.065 0.744\n190SOL O 0.838 0.840 0.838\n190SOL H 0.843 0.829 0.853\n190SOL H 0.869 0.864 0.835\n191SOL O 0.867 0.774 0.530\n191SOL H 0.857 0.813 0.529\n191SOL H 0.873 0.767 0.545\n192SOL O 0.173 0.665 0.985\n192SOL H 0.184 0.692 0.996\n192SOL H 0.133 0.670 0.984\n193SOL O 0.371 0.649 0.611\n193SOL H 0.394 0.643 0.624\n193SOL H 0.398 0.648 0.599\n194SOL O 0.503 0.966 0.979\n194SOL H 0.469 0.974 0.971\n194SOL H 0.525 0.999 0.977\n195SOL O 0.577 0.866 0.622\n195SOL H 0.583 0.835 0.632\n195SOL H 0.602 0.894 0.628\n196SOL O 0.258 0.185 0.734\n196SOL H 0.288 0.173 0.743\n196SOL H 0.257 0.225 0.735\n197SOL O 0.000 0.637 0.454\n197SOL H 0.009 0.665 0.443\n197SOL H 0.024 0.606 0.451\n198SOL O 0.707 0.371 0.661\n198SOL H 0.708 0.333 0.666\n198SOL H 0.718 0.392 0.673\n199SOL O 0.794 0.128 0.913\n199SOL H 0.787 0.097 0.923\n199SOL H 0.768 0.156 0.918\n200SOL O 0.523 0.763 0.915\n200SOL H 0.536 0.766 0.900\n200SOL H 0.541 0.730 0.920\n201SOL O 0.767 0.027 0.553\n201SOL H 0.747 0.018 0.566\n201SOL H 0.754 0.000 0.542\n202SOL O 0.451 0.816 0.588\n202SOL H 0.487 0.813 0.594\n202SOL H 0.458 0.825 0.573\n203SOL O 0.967 0.567 0.956\n203SOL H 0.937 0.543 0.961\n203SOL H 1.000 0.547 0.960\n204SOL O 0.754 0.413 0.913\n204SOL H 0.790 0.430 0.912\n204SOL H 0.732 0.429 0.902\n205SOL O 0.040 0.304 0.992\n205SOL H 0.061 0.333 0.999\n205SOL H 0.003 0.307 0.998\n206SOL O 0.112 0.341 0.727\n206SOL H 0.109 0.330 0.712\n206SOL H 0.121 0.307 0.735\n207SOL O 0.498 0.424 0.986\n207SOL H 0.533 0.405 0.985\n207SOL H 0.502 0.444 1.000\n208SOL O 0.964 0.499 0.709\n208SOL H 0.995 0.520 0.704\n208SOL H 0.980 0.467 0.716\n209SOL O 0.979 0.427 0.929\n209SOL H 0.988 0.394 0.937\n209SOL H 1.000 0.456 0.936\n210SOL O 0.317 0.918 0.878\n210SOL H 0.292 0.903 0.889\n210SOL H 0.325 0.888 0.868\n211SOL O 0.209 0.698 0.635\n211SOL H 0.185 0.666 0.638\n211SOL H 0.194 0.714 0.622\n212SOL O 0.765 0.808 0.552\n212SOL H 0.767 0.807 0.536\n212SOL H 0.774 0.770 0.556\n213SOL O 0.640 1.000 0.819\n213SOL H 0.640 0.968 0.828\n213SOL H 0.606 0.998 0.811\n214SOL O 0.615 0.952 0.763\n214SOL H 0.602 0.976 0.751\n214SOL H 0.632 0.976 0.773\n215SOL O 0.713 0.144 0.822\n215SOL H 0.741 0.145 0.811\n215SOL H 0.703 0.106 0.823\n216SOL O 0.359 0.521 0.679\n216SOL H 0.370 0.493 0.668\n216SOL H 0.319 0.518 0.680\n217SOL O 0.301 0.158 0.436\n217SOL H 0.275 0.133 0.443\n217SOL H 0.297 0.151 0.421\n218SOL O 0.093 0.897 0.541\n218SOL H 0.094 0.931 0.549\n218SOL H 0.056 0.883 0.544\n219SOL O 0.243 0.087 0.884\n219SOL H 0.248 0.051 0.877\n219SOL H 0.240 0.112 0.872\n220SOL O 0.098 0.138 0.707\n220SOL H 0.110 0.141 0.722\n220SOL H 0.131 0.125 0.699\n221SOL O 0.885 0.732 0.951\n221SOL H 0.909 0.733 0.964\n221SOL H 0.852 0.752 0.956\n222SOL O 0.340 0.728 0.513\n222SOL H 0.362 0.721 0.526\n222SOL H 0.318 0.761 0.517\n223SOL O 0.529 0.754 0.452\n223SOL H 0.559 0.759 0.462\n223SOL H 0.498 0.774 0.458\n224SOL O 0.899 0.876 0.781\n224SOL H 0.926 0.890 0.791\n224SOL H 0.904 0.898 0.768\n225SOL O 0.641 0.811 0.723\n225SOL H 0.673 0.832 0.728\n225SOL H 0.645 0.811 0.707\n226SOL O 0.476 0.968 0.538\n226SOL H 0.505 0.978 0.548\n226SOL H 0.453 1.000 0.537\n227SOL O 0.803 0.933 1.000\n227SOL H 0.795 0.973 0.999\n227SOL H 0.843 0.931 1.000\n228SOL O 0.335 0.956 0.463\n228SOL H 0.300 0.937 0.461\n228SOL H 0.352 0.934 0.474\n229SOL O 0.079 0.734 0.737\n229SOL H 0.066 0.760 0.726\n229SOL H 0.116 0.747 0.740\n230SOL O 0.444 0.685 0.806\n230SOL H 0.463 0.657 0.798\n230SOL H 0.406 0.684 0.801\n231SOL O 0.840 0.941 0.660\n231SOL H 0.802 0.954 0.661\n231SOL H 0.854 0.957 0.647\n232SOL O 0.130 0.657 0.680\n232SOL H 0.132 0.681 0.667\n232SOL H 0.167 0.658 0.686\n233SOL O 0.303 0.636 0.778\n233SOL H 0.315 0.663 0.789\n233SOL H 0.335 0.614 0.775\n234SOL O 0.747 0.357 0.618\n234SOL H 0.769 0.357 0.631\n234SOL H 0.714 0.378 0.622\n235SOL O 0.407 1.000 0.428\n235SOL H 0.408 0.977 0.441\n235SOL H 0.444 0.998 0.422\n236SOL O 0.163 0.514 0.593\n236SOL H 0.157 0.481 0.602\n236SOL H 0.126 0.524 0.588\n237SOL O 0.338 0.424 0.512\n237SOL H 0.374 0.441 0.514\n237SOL H 0.313 0.455 0.511\n238SOL O 0.571 0.094 0.540\n238SOL H 0.557 0.059 0.534\n238SOL H 0.611 0.092 0.538\n239SOL O 0.845 0.441 0.596\n239SOL H 0.837 0.401 0.595\n239SOL H 0.814 0.459 0.589\n240SOL O 0.849 0.622 0.500\n240SOL H 0.843 0.638 0.486\n240SOL H 0.873 0.648 0.508\n241SOL O 0.446 0.863 0.770\n241SOL H 0.451 0.831 0.780\n241SOL H 0.451 0.895 0.779\n242SOL O 0.301 0.467 0.634\n242SOL H 0.317 0.442 0.644\n242SOL H 0.312 0.504 0.639\n243SOL O 0.147 0.000 0.879\n243SOL H 0.109 0.006 0.884\n243SOL H 0.152 0.021 0.866\n244SOL O 0.643 0.333 0.947\n244SOL H 0.616 0.328 0.959\n244SOL H 0.667 0.301 0.948\n245SOL O 0.958 0.236 0.691\n245SOL H 0.945 0.216 0.703\n245SOL H 0.998 0.236 0.692\n246SOL O 0.397 0.223 0.749\n246SOL H 0.392 0.249 0.737\n246SOL H 0.361 0.224 0.757\n247SOL O 0.503 0.619 0.699\n247SOL H 0.520 0.648 0.708\n247SOL H 0.528 0.588 0.701\n248SOL O 0.813 0.847 0.957\n248SOL H 0.825 0.832 0.971\n248SOL H 0.792 0.880 0.960\n249SOL O 0.720 0.468 0.420\n249SOL H 0.696 0.449 0.430\n249SOL H 0.703 0.504 0.418\n250SOL O 0.973 0.507 0.835\n250SOL H 0.971 0.548 0.834\n250SOL H 0.965 0.499 0.850\n251SOL O 0.047 0.286 0.559\n251SOL H 0.020 0.288 0.570\n251SOL H 0.069 0.253 0.562\n252SOL O 0.935 0.954 0.914\n252SOL H 0.914 0.988 0.911\n252SOL H 0.956 0.948 0.901\n253SOL O 0.244 0.388 0.834\n253SOL H 0.241 0.413 0.846\n253SOL H 0.283 0.377 0.833\n254SOL O 0.982 0.959 0.824\n254SOL H 0.989 0.977 0.838\n254SOL H 0.947 0.972 0.819\n255SOL O 0.395 0.640 0.449\n255SOL H 0.363 0.635 0.458\n255SOL H 0.421 0.613 0.454\n256SOL O 0.525 0.793 0.823\n256SOL H 0.546 0.767 0.832\n256SOL H 0.517 0.772 0.810\n257SOL O 0.820 0.543 0.819\n257SOL H 0.814 0.506 0.825\n257SOL H 0.848 0.539 0.808\n258SOL O 0.618 0.793 0.413\n258SOL H 0.647 0.765 0.412\n258SOL H 0.638 0.828 0.412\n259SOL O 0.534 0.258 0.499\n259SOL H 0.547 0.251 0.484\n259SOL H 0.499 0.277 0.497\n260SOL O 0.076 0.658 0.640\n260SOL H 0.096 0.689 0.633\n260SOL H 0.038 0.671 0.642\n261SOL O 0.253 0.187 0.613\n261SOL H 0.225 0.177 0.602\n261SOL H 0.261 0.226 0.610\n262SOL O 0.709 0.180 0.751\n262SOL H 0.688 0.169 0.764\n262SOL H 0.704 0.149 0.741\n263SOL O 0.429 0.024 0.883\n263SOL H 0.435 0.019 0.898\n263SOL H 0.430 0.064 0.881\n264SOL O 0.765 0.135 0.545\n264SOL H 0.792 0.131 0.533\n264SOL H 0.736 0.109 0.541\n265SOL O 0.979 0.211 0.487\n265SOL H 0.960 0.247 0.487\n265SOL H 1.000 0.210 0.500\n266SOL O 0.155 0.619 0.932\n266SOL H 0.140 0.584 0.928\n266SOL H 0.171 0.613 0.946\n267SOL O 0.488 0.905 0.878\n267SOL H 0.511 0.930 0.887\n267SOL H 0.451 0.908 0.885\n268SOL O 0.690 0.293 0.604\n268SOL H 0.687 0.256 0.599\n268SOL H 0.669 0.293 0.617\n269SOL O 0.606 0.471 0.421\n269SOL H 0.567 0.465 0.419\n269SOL H 0.609 0.480 0.437\n270SOL O 0.987 0.089 0.638\n270SOL H 0.947 0.095 0.637\n270SOL H 1.000 0.120 0.646\n271SOL O 0.592 0.828 0.539\n271SOL H 0.576 0.844 0.526\n271SOL H 0.603 0.859 0.548\n272SOL O 0.486 0.155 0.537\n272SOL H 0.485 0.125 0.548\n272SOL H 0.460 0.182 0.543\n273SOL O 0.518 0.969 0.617\n273SOL H 0.531 1.001 0.625\n273SOL H 0.520 0.938 0.627\n274SOL O 0.336 0.013 0.743\n274SOL H 0.343 0.049 0.749\n274SOL H 0.362 0.009 0.731\n275SOL O 0.355 0.546 0.989\n275SOL H 0.328 0.555 1.000\n275SOL H 0.363 0.507 0.991\n276SOL O 0.749 0.093 0.417\n276SOL H 0.753 0.101 0.432\n276SOL H 0.711 0.100 0.414\n277SOL O 0.575 0.788 0.760\n277SOL H 0.595 0.820 0.754\n277SOL H 0.589 0.757 0.752\n278SOL O 0.021 0.399 0.879\n278SOL H 0.020 0.437 0.885\n278SOL H 0.014 0.377 0.892\n279SOL O 0.879 0.740 0.662\n279SOL H 0.879 0.715 0.649\n279SOL H 0.856 0.771 0.657\n280SOL O 0.824 0.593 0.981\n280SOL H 0.851 0.613 0.990\n280SOL H 0.804 0.621 0.973\n281SOL O 0.703 0.924 0.797\n281SOL H 0.704 0.964 0.796\n281SOL H 0.665 0.915 0.800\n282SOL O 0.596 0.435 0.691\n282SOL H 0.620 0.466 0.687\n282SOL H 0.609 0.404 0.683\n283SOL O 1.001 0.855 0.755\n283SOL H 0.973 0.837 0.763\n283SOL H 0.987 0.893 0.753\n284SOL O 0.224 0.691 0.844\n284SOL H 0.251 0.688 0.832\n284SOL H 0.189 0.681 0.837\n285SOL O 1.000 0.760 0.900\n285SOL H 0.962 0.768 0.903\n285SOL H 1.000 0.723 0.894\n286SOL O 0.048 0.193 0.984\n286SOL H 0.060 0.215 0.997\n286SOL H 0.028 0.219 0.975\n287SOL O 0.944 0.783 0.701\n287SOL H 0.979 0.763 0.701\n287SOL H 0.918 0.758 0.694\n288SOL O 0.019 0.529 0.968\n288SOL H 0.022 0.546 0.982\n288SOL H 0.027 0.490 0.970\n289SOL O 0.394 0.386 0.484\n289SOL H 0.396 0.362 0.497\n289SOL H 0.432 0.398 0.482\n290SOL O 0.472 0.171 0.920\n290SOL H 0.487 0.138 0.912\n290SOL H 0.501 0.185 0.928\n291SOL O 0.621 0.507 0.538\n291SOL H 0.620 0.467 0.542\n291SOL H 0.651 0.510 0.528\n292SOL O 0.832 0.075 0.595\n292SOL H 0.805 0.074 0.607\n292SOL H 0.864 0.057 0.601\n293SOL O 0.706 0.801 0.942\n293SOL H 0.724 0.829 0.934\n293SOL H 0.696 0.819 0.956\n294SOL O 0.345 0.367 0.792\n294SOL H 0.373 0.375 0.803\n294SOL H 0.310 0.376 0.798\n295SOL O 0.177 0.623 0.748\n295SOL H 0.152 0.593 0.746\n295SOL H 0.165 0.652 0.738\n296SOL O 0.482 0.610 0.609\n296SOL H 0.475 0.600 0.624\n296SOL H 0.478 0.576 0.601\n297SOL O 0.065 0.197 0.626\n297SOL H 0.063 0.184 0.641\n297SOL H 0.052 0.235 0.627\n298SOL O 0.332 0.611 0.965\n298SOL H 0.360 0.601 0.955\n298SOL H 0.301 0.588 0.962\n299SOL O 0.994 0.717 0.655\n299SOL H 0.984 0.683 0.662\n299SOL H 0.971 0.745 0.662\n300SOL O 0.611 0.024 0.504\n300SOL H 0.626 0.003 0.492\n300SOL H 0.576 0.037 0.499\n301SOL O 0.208 0.521 0.857\n301SOL H 0.235 0.495 0.863\n301SOL H 0.210 0.516 0.841\n302SOL O 0.204 0.250 0.778\n302SOL H 0.170 0.236 0.785\n302SOL H 0.201 0.243 0.763\n303SOL O 0.265 0.599 0.462\n303SOL H 0.236 0.579 0.469\n303SOL H 0.293 0.572 0.459\n304SOL O 0.907 0.695 0.909\n304SOL H 0.909 0.672 0.896\n304SOL H 0.881 0.724 0.906\n305SOL O 0.113 0.502 0.425\n305SOL H 0.095 0.525 0.436\n305SOL H 0.110 0.523 0.412\n306SOL O 0.085 0.433 0.534\n306SOL H 0.053 0.423 0.525\n306SOL H 0.076 0.470 0.539\n307SOL O 0.184 0.135 0.774\n307SOL H 0.153 0.123 0.765\n307SOL H 0.213 0.147 0.764\n308SOL O 0.955 0.679 0.783\n308SOL H 0.962 0.670 0.768\n308SOL H 0.943 0.645 0.789\n309SOL O 0.390 0.576 0.544\n309SOL H 0.421 0.552 0.541\n309SOL H 0.399 0.612 0.538\n310SOL O 0.034 0.234 0.923\n310SOL H 0.019 0.213 0.910\n310SOL H 0.009 0.265 0.924\n311SOL O 0.526 0.197 0.804\n311SOL H 0.544 0.222 0.815\n311SOL H 0.488 0.197 0.809\n312SOL O 0.699 0.261 0.458\n312SOL H 0.727 0.267 0.447\n312SOL H 0.696 0.296 0.466\n313SOL O 0.103 0.413 0.693\n313SOL H 0.078 0.419 0.681\n313SOL H 0.134 0.437 0.690\n314SOL O 0.659 0.000 0.874\n314SOL H 0.641 0.026 0.864\n314SOL H 0.697 0.000 0.870\n315SOL O 0.258 0.838 0.540\n315SOL H 0.291 0.841 0.532\n315SOL H 0.259 0.801 0.546\n316SOL O 0.442 0.181 0.422\n316SOL H 0.438 0.213 0.412\n316SOL H 0.441 0.149 0.413\n317SOL O 0.199 0.869 0.654\n317SOL H 0.193 0.839 0.643\n317SOL H 0.178 0.900 0.648\n318SOL O 0.842 0.027 0.428\n318SOL H 0.860 0.008 0.416\n318SOL H 0.803 0.015 0.427\n319SOL O 0.515 0.294 0.545\n319SOL H 0.488 0.269 0.552\n319SOL H 0.547 0.290 0.554\n320SOL O 0.957 0.744 0.531\n320SOL H 0.987 0.742 0.520\n320SOL H 0.934 0.711 0.527\n321SOL O 0.982 0.265 0.530\n321SOL H 0.963 0.247 0.542\n321SOL H 0.983 0.304 0.534\n322SOL O 0.058 0.413 0.778\n322SOL H 0.074 0.378 0.782\n322SOL H 0.087 0.440 0.779\n323SOL O 0.684 0.702 0.575\n323SOL H 0.651 0.709 0.566\n323SOL H 0.714 0.703 0.564\n324SOL O 0.939 0.963 0.739\n324SOL H 0.965 0.957 0.728\n324SOL H 0.925 1.000 0.737\n325SOL O 0.242 0.870 0.690\n325SOL H 0.269 0.867 0.701\n325SOL H 0.207 0.875 0.697\n326SOL O 0.694 0.101 0.467\n326SOL H 0.667 0.125 0.474\n326SOL H 0.688 0.064 0.473\n327SOL O 0.596 0.475 0.997\n327SOL H 0.577 0.504 0.988\n327SOL H 0.614 0.452 0.986\n328SOL O 0.931 0.777 0.801\n328SOL H 0.909 0.803 0.809\n328SOL H 0.919 0.741 0.806\n329SOL O 0.891 0.664 0.455\n329SOL H 0.918 0.687 0.462\n329SOL H 0.907 0.656 0.441\n330SOL O 0.596 0.583 0.626\n330SOL H 0.622 0.595 0.637\n330SOL H 0.567 0.611 0.626\n331SOL O 0.175 0.637 0.880\n331SOL H 0.176 0.632 0.895\n331SOL H 0.213 0.643 0.876\n332SOL O 0.798 0.805 0.418\n332SOL H 0.837 0.798 0.420\n332SOL H 0.783 0.771 0.412\n333SOL O 0.786 0.164 0.873\n333SOL H 0.764 0.196 0.869\n333SOL H 0.760 0.137 0.878\n334SOL O 0.204 0.297 0.816\n334SOL H 0.219 0.328 0.808\n334SOL H 0.228 0.266 0.813\n335SOL O 0.647 0.337 0.901\n335SOL H 0.635 0.302 0.907\n335SOL H 0.683 0.343 0.907\n336SOL O 0.978 0.317 0.441\n336SOL H 0.946 0.300 0.434\n336SOL H 0.965 0.325 0.456\n337SOL O 0.307 0.799 0.984\n337SOL H 0.280 0.792 0.973\n337SOL H 0.339 0.778 0.979\n338SOL O 0.993 0.868 0.639\n338SOL H 1.001 0.884 0.653\n338SOL H 0.960 0.886 0.634\n339SOL O 0.158 0.184 0.560\n339SOL H 0.155 0.215 0.569\n339SOL H 0.125 0.162 0.562\n340SOL O 0.208 0.734 0.891\n340SOL H 0.233 0.711 0.899\n340SOL H 0.230 0.758 0.883\n341SOL O 0.576 0.951 0.516\n341SOL H 0.605 0.932 0.524\n341SOL H 0.573 0.987 0.523\n342SOL O 0.925 0.325 0.918\n342SOL H 0.901 0.341 0.929\n342SOL H 0.923 0.286 0.920\n343SOL O 0.439 0.513 0.787\n343SOL H 0.410 0.533 0.779\n343SOL H 0.461 0.542 0.793\n344SOL O 0.340 0.843 0.822\n344SOL H 0.349 0.816 0.834\n344SOL H 0.366 0.872 0.824\n345SOL O 0.810 0.392 0.742\n345SOL H 0.800 0.427 0.749\n345SOL H 0.821 0.367 0.753\n346SOL O 0.205 0.344 0.938\n346SOL H 0.218 0.328 0.952\n346SOL H 0.176 0.318 0.933\n347SOL O 0.061 0.897 1.000\n347SOL H 0.074 0.859 0.999\n347SOL H 0.081 0.917 0.989\n348SOL O 0.648 0.281 0.872\n348SOL H 0.649 0.302 0.858\n348SOL H 0.686 0.271 0.875\n349SOL O 0.866 0.311 0.627\n349SOL H 0.861 0.349 0.631\n349SOL H 0.885 0.293 0.638\n350SOL O 0.755 0.135 0.986\n350SOL H 0.765 0.144 0.971\n350SOL H 0.715 0.129 0.985\n351SOL O 0.555 0.377 0.542\n351SOL H 0.551 0.411 0.534\n351SOL H 0.534 0.384 0.555\n352SOL O 0.467 0.136 0.835\n352SOL H 0.500 0.116 0.840\n352SOL H 0.463 0.166 0.845\n353SOL O 0.174 0.763 0.673\n353SOL H 0.148 0.774 0.662\n353SOL H 0.200 0.793 0.675\n354SOL O 0.570 0.751 0.595\n354SOL H 0.560 0.768 0.608\n354SOL H 0.536 0.736 0.589\n355SOL O 0.000 0.335 0.748\n355SOL H 0.000 0.352 0.762\n355SOL H 0.022 0.359 0.739\n356SOL O 0.853 0.972 0.781\n356SOL H 0.829 0.953 0.770\n356SOL H 0.872 1.000 0.772\n357SOL O 0.586 0.310 0.594\n357SOL H 0.613 0.340 0.595\n357SOL H 0.603 0.283 0.585\n358SOL O 0.808 0.702 0.535\n358SOL H 0.815 0.664 0.537\n358SOL H 0.808 0.707 0.519\n359SOL O 0.679 0.612 0.783\n359SOL H 0.668 0.625 0.769\n359SOL H 0.713 0.592 0.780\n360SOL O 0.190 0.863 0.563\n360SOL H 0.182 0.888 0.552\n360SOL H 0.167 0.831 0.561\n361SOL O 0.017 0.485 0.690\n361SOL H 0.013 0.523 0.685\n361SOL H 0.041 0.486 0.702\n362SOL O 0.112 0.159 0.430\n362SOL H 0.127 0.185 0.440\n362SOL H 0.072 0.158 0.433\n363SOL O 0.822 0.602 0.939\n363SOL H 0.782 0.605 0.939\n363SOL H 0.831 0.582 0.925\n364SOL O 0.806 0.726 0.863\n364SOL H 0.779 0.699 0.858\n364SOL H 0.817 0.712 0.877\n365SOL O 0.295 0.608 0.628\n365SOL H 0.324 0.595 0.638\n365SOL H 0.284 0.644 0.634\n366SOL O 0.392 0.787 0.690\n366SOL H 0.412 0.807 0.679\n366SOL H 0.389 0.749 0.684\n367SOL O 0.829 0.962 0.440\n367SOL H 0.827 0.983 0.453\n367SOL H 0.827 0.924 0.444\n368SOL O 0.431 0.477 0.969\n368SOL H 0.405 0.451 0.963\n368SOL H 0.419 0.512 0.963\n369SOL O 0.200 0.378 0.988\n369SOL H 0.172 0.377 0.977\n369SOL H 0.195 0.343 0.996\n370SOL O 0.561 0.409 0.796\n370SOL H 0.541 0.420 0.809\n370SOL H 0.556 0.369 0.795\n371SOL O 0.309 0.704 0.861\n371SOL H 0.315 0.740 0.868\n371SOL H 0.346 0.688 0.859\n372SOL O 0.096 0.280 0.799\n372SOL H 0.079 0.308 0.808\n372SOL H 0.101 0.248 0.808\n373SOL O 0.125 0.552 0.510\n373SOL H 0.123 0.546 0.525\n373SOL H 0.162 0.566 0.507\n374SOL O 0.469 0.379 0.522\n374SOL H 0.491 0.378 0.509\n374SOL H 0.455 0.342 0.524\n375SOL O 0.181 0.072 0.910\n375SOL H 0.184 0.076 0.926\n375SOL H 0.195 0.035 0.908\n376SOL O 0.966 0.553 0.797\n376SOL H 0.992 0.553 0.785\n376SOL H 0.930 0.552 0.791\n377SOL O 0.262 0.097 0.492\n377SOL H 0.232 0.117 0.500\n377SOL H 0.257 0.059 0.497\n378SOL O 0.319 0.134 0.787\n378SOL H 0.285 0.156 0.786\n378SOL H 0.311 0.105 0.798\n379SOL O 0.146 0.784 0.952\n379SOL H 0.124 0.758 0.943\n379SOL H 0.184 0.779 0.947\n380SOL O 0.955 0.970 0.525\n380SOL H 0.954 1.000 0.515\n380SOL H 0.993 0.956 0.525\n381SOL O 0.970 0.939 0.956\n381SOL H 0.991 0.923 0.944\n381SOL H 0.953 0.908 0.963\n382SOL O 0.267 0.159 0.990\n382SOL H 0.240 0.145 1.000\n382SOL H 0.302 0.145 0.995\n383SOL O 0.392 0.291 0.548\n383SOL H 0.389 0.264 0.560\n383SOL H 0.407 0.323 0.555\n384SOL O 0.853 0.045 0.982\n384SOL H 0.820 0.064 0.987\n384SOL H 0.845 0.007 0.984\n385SOL O 0.467 0.295 0.924\n385SOL H 0.430 0.311 0.927\n385SOL H 0.459 0.261 0.916\n386SOL O 0.712 0.382 0.454\n386SOL H 0.750 0.381 0.460\n386SOL H 0.694 0.409 0.463\n387SOL O 0.531 0.988 0.855\n387SOL H 0.561 0.968 0.847\n387SOL H 0.498 0.978 0.847\n388SOL O 0.059 0.244 0.414\n388SOL H 0.098 0.239 0.412\n388SOL H 0.054 0.283 0.417\n389SOL O 1.000 0.334 0.973\n389SOL H 0.983 0.320 0.960\n389SOL H 0.969 0.344 0.982\n390SOL O 0.273 0.064 0.934\n390SOL H 0.293 0.085 0.923\n390SOL H 0.280 0.026 0.931\n391SOL O 0.848 0.965 0.695\n391SOL H 0.815 0.977 0.703\n391SOL H 0.864 1.000 0.690\n392SOL O 0.940 0.805 0.868\n392SOL H 0.945 0.845 0.870\n392SOL H 0.907 0.797 0.877\n393SOL O 0.424 0.782 0.817\n393SOL H 0.438 0.819 0.815\n393SOL H 0.401 0.776 0.804\n394SOL O 0.657 0.145 0.857\n394SOL H 0.673 0.120 0.867\n394SOL H 0.675 0.180 0.860\n395SOL O 0.222 0.345 0.432\n395SOL H 0.230 0.340 0.417\n395SOL H 0.190 0.322 0.435\n396SOL O 0.401 0.656 0.671\n396SOL H 0.394 0.677 0.658\n396SOL H 0.411 0.619 0.666\n397SOL O 0.400 0.482 0.565\n397SOL H 0.391 0.458 0.577\n397SOL H 0.440 0.481 0.564\n398SOL O 0.684 0.101 0.591\n398SOL H 0.666 0.069 0.586\n398SOL H 0.715 0.088 0.600\n399SOL O 0.568 0.061 0.643\n399SOL H 0.561 0.041 0.629\n399SOL H 0.541 0.092 0.642\n400SOL O 0.352 0.000 0.454\n400SOL H 0.322 0.001 0.444\n400SOL H 0.384 0.011 0.446\n401SOL O 0.360 0.955 0.954\n401SOL H 0.381 0.944 0.967\n401SOL H 0.367 0.994 0.952\n402SOL O 0.245 0.836 0.600\n402SOL H 0.226 0.805 0.606\n402SOL H 0.234 0.867 0.608\n403SOL O 0.595 0.509 0.774\n403SOL H 0.556 0.513 0.771\n403SOL H 0.613 0.541 0.768\n404SOL O 0.080 0.957 0.896\n404SOL H 0.060 0.979 0.886\n404SOL H 0.077 0.977 0.910\n405SOL O 0.791 0.254 0.625\n405SOL H 0.768 0.235 0.615\n405SOL H 0.822 0.229 0.627\n406SOL O 0.259 0.854 0.943\n406SOL H 0.267 0.879 0.931\n406SOL H 0.222 0.863 0.948\n407SOL O 0.473 0.136 0.961\n407SOL H 0.494 0.163 0.969\n407SOL H 0.441 0.128 0.970\n408SOL O 0.563 0.872 0.793\n408SOL H 0.592 0.877 0.783\n408SOL H 0.541 0.842 0.788\n409SOL O 0.619 0.132 0.965\n409SOL H 0.629 0.103 0.975\n409SOL H 0.581 0.124 0.961\n410SOL O 0.871 0.579 0.533\n410SOL H 0.893 0.554 0.524\n410SOL H 0.833 0.569 0.530\n411SOL O 0.626 0.125 0.804\n411SOL H 0.598 0.138 0.794\n411SOL H 0.621 0.148 0.817\n412SOL O 0.408 0.206 0.991\n412SOL H 0.388 0.179 1.000\n412SOL H 0.406 0.239 1.000\n413SOL O 0.397 0.356 0.656\n413SOL H 0.395 0.342 0.642\n413SOL H 0.405 0.395 0.655\n414SOL O 0.886 0.020 0.847\n414SOL H 0.915 0.021 0.836\n414SOL H 0.904 0.033 0.860\n415SOL O 0.579 0.862 0.442\n415SOL H 0.559 0.852 0.429\n415SOL H 0.600 0.830 0.445\n416SOL O 0.000 0.093 0.967\n416SOL H 0.035 0.102 0.960\n416SOL H -0.000 0.115 0.980\n417SOL O 0.802 0.487 0.548\n417SOL H 0.776 0.498 0.536\n417SOL H 0.792 0.510 0.560\n418SOL O 0.185 0.442 0.513\n418SOL H 0.195 0.403 0.510\n418SOL H 0.221 0.460 0.515\n419SOL O 0.020 0.217 0.796\n419SOL H 0.000 0.191 0.787\n419SOL H 0.000 0.251 0.794\n420SOL O 0.869 0.379 0.665\n420SOL H 0.833 0.361 0.665\n420SOL H 0.872 0.395 0.680\n421SOL O 0.189 0.564 0.992\n421SOL H 0.200 0.568 0.977\n421SOL H 0.223 0.570 1.000\n422SOL O 0.554 0.065 0.710\n422SOL H 0.532 0.045 0.721\n422SOL H 0.590 0.047 0.711\n423SOL O 0.746 0.821 0.801\n423SOL H 0.712 0.824 0.793\n423SOL H 0.735 0.808 0.816\n424SOL O 0.403 0.526 0.477\n424SOL H 0.440 0.541 0.475\n424SOL H 0.408 0.487 0.477\n425SOL O 0.072 0.172 0.765\n425SOL H 0.042 0.173 0.755\n425SOL H 0.091 0.207 0.764\n426SOL O 0.506 0.733 0.671\n426SOL H 0.491 0.757 0.660\n426SOL H 0.490 0.698 0.668\n427SOL O 0.580 0.033 0.837\n427SOL H 0.554 0.019 0.826\n427SOL H 0.582 0.072 0.834\n428SOL O 0.392 0.218 0.948\n428SOL H 0.390 0.217 0.933\n428SOL H 0.431 0.224 0.951\n429SOL O 0.900 0.198 0.869\n429SOL H 0.876 0.168 0.864\n429SOL H 0.879 0.231 0.866\n430SOL O 0.535 0.349 0.702\n430SOL H 0.557 0.330 0.713\n430SOL H 0.530 0.386 0.708\n431SOL O 0.228 0.489 0.672\n431SOL H 0.195 0.509 0.676\n431SOL H 0.232 0.495 0.656\n432SOL O 0.740 0.585 0.864\n432SOL H 0.771 0.610 0.867\n432SOL H 0.752 0.562 0.852\n433SOL O 0.227 0.160 0.648\n433SOL H 0.240 0.187 0.658\n433SOL H 0.189 0.152 0.651\n434SOL O 0.693 0.379 0.710\n434SOL H 0.680 0.416 0.707\n434SOL H 0.665 0.355 0.705\n435SOL O 0.708 0.710 0.513\n435SOL H 0.718 0.691 0.526\n435SOL H 0.724 0.747 0.515\n436SOL O 0.926 0.924 0.557\n436SOL H 0.941 0.938 0.571\n436SOL H 0.886 0.921 0.560\n437SOL O 1.000 0.408 0.487\n437SOL H 0.982 0.440 0.494\n437SOL H 1.000 0.414 0.472\n438SOL O 0.275 0.475 0.987\n438SOL H 0.251 0.451 0.978\n438SOL H 0.253 0.482 1.000\n439SOL O 0.323 0.451 0.471\n439SOL H 0.332 0.421 0.462\n439SOL H 0.310 0.480 0.462\n440SOL O 0.879 0.565 0.456\n440SOL H 0.857 0.553 0.468\n440SOL H 0.889 0.532 0.448\n441SOL O 0.628 0.353 0.821\n441SOL H 0.665 0.349 0.827\n441SOL H 0.613 0.386 0.827\n442SOL O 0.869 0.219 0.502\n442SOL H 0.868 0.234 0.487\n442SOL H 0.896 0.189 0.501\n443SOL O 0.910 0.509 0.555\n443SOL H 0.883 0.525 0.565\n443SOL H 0.933 0.486 0.564\n444SOL O 0.109 0.216 0.872\n444SOL H 0.119 0.254 0.872\n444SOL H 0.081 0.212 0.883\n445SOL O 0.291 1.000 0.920\n445SOL H 0.272 0.978 0.908\n445SOL H 0.285 0.979 0.933\n446SOL O 0.419 0.147 0.573\n446SOL H 0.424 0.171 0.586\n446SOL H 0.381 0.135 0.574\n447SOL O 0.383 0.403 0.606\n447SOL H 0.355 0.378 0.612\n447SOL H 0.418 0.392 0.612\n448SOL O 0.760 0.880 0.704\n448SOL H 0.752 0.912 0.712\n448SOL H 0.758 0.849 0.714\n449SOL O 0.878 0.648 0.576\n449SOL H 0.844 0.669 0.575\n449SOL H 0.901 0.664 0.564\n450SOL O 0.675 0.925 0.865\n450SOL H 0.649 0.950 0.872\n450SOL H 0.669 0.890 0.872\n451SOL O 0.050 0.560 0.830\n451SOL H 0.035 0.593 0.823\n451SOL H 0.028 0.556 0.843\n452SOL O 0.762 0.072 0.852\n452SOL H 0.795 0.087 0.844\n452SOL H 0.756 0.036 0.846\n453SOL O 0.611 0.611 0.943\n453SOL H 0.642 0.601 0.953\n453SOL H 0.595 0.643 0.950\n454SOL O 0.896 0.043 0.539\n454SOL H 0.936 0.046 0.539\n454SOL H 0.887 0.024 0.553\n455SOL O 0.435 0.684 0.497\n455SOL H 0.429 0.722 0.501\n455SOL H 0.419 0.663 0.509\n456SOL O 0.742 0.511 0.740\n456SOL H 0.723 0.547 0.741\n456SOL H 0.737 0.494 0.754\n457SOL O 0.467 0.513 0.740\n457SOL H 0.477 0.533 0.754\n457SOL H 0.495 0.485 0.739\n458SOL O 0.788 0.096 0.760\n458SOL H 0.784 0.127 0.771\n458SOL H 0.758 0.071 0.764\n459SOL O 0.559 0.472 0.580\n459SOL H 0.548 0.480 0.565\n459SOL H 0.589 0.446 0.579\n460SOL O 0.395 0.240 0.516\n460SOL H 0.358 0.225 0.514\n460SOL H 0.418 0.220 0.506\n461SOL O 0.977 0.105 0.498\n461SOL H 0.942 0.110 0.506\n461SOL H 1.000 0.084 0.508\n462SOL O 0.739 0.971 0.970\n462SOL H 0.751 1.000 0.960\n462SOL H 0.721 0.944 0.960\n463SOL O 0.113 0.882 0.830\n463SOL H 0.103 0.845 0.824\n463SOL H 0.083 0.905 0.825\n464SOL O 0.199 0.056 0.618\n464SOL H 0.211 0.030 0.629\n464SOL H 0.216 0.091 0.622\n465SOL O 0.402 0.307 0.892\n465SOL H 0.398 0.268 0.889\n465SOL H 0.439 0.317 0.888\n466SOL O 0.613 0.717 0.995\n466SOL H 0.592 0.749 1.000\n466SOL H 0.600 0.709 0.981\n467SOL O 0.543 0.839 0.978\n467SOL H 0.531 0.814 0.966\n467SOL H 0.583 0.835 0.978\n468SOL O 0.379 0.718 0.908\n468SOL H 0.400 0.752 0.908\n468SOL H 0.343 0.728 0.902\n469SOL O 0.544 0.544 0.454\n469SOL H 0.582 0.548 0.460\n469SOL H 0.535 0.506 0.457\n470SOL O 0.808 0.700 1.000\n470SOL H 0.768 0.693 1.000\n470SOL H 0.818 0.709 0.985\n471SOL O 0.203 0.269 0.632\n471SOL H 0.211 0.273 0.648\n471SOL H 0.187 0.233 0.631\n472SOL O 0.063 0.431 0.472\n472SOL H 0.036 0.402 0.475\n472SOL H 0.060 0.456 0.484\n473SOL O 0.233 0.788 0.818\n473SOL H 0.207 0.777 0.829\n473SOL H 0.269 0.777 0.824\n474SOL O 0.935 0.115 0.673\n474SOL H 0.950 0.152 0.674\n474SOL H 0.900 0.117 0.680\n475SOL O 0.891 0.407 0.550\n475SOL H 0.907 0.427 0.538\n475SOL H 0.909 0.371 0.550\n476SOL O 0.625 0.970 0.608\n476SOL H 0.603 1.000 0.603\n476SOL H 0.648 0.958 0.596\n477SOL O 0.228 0.117 0.563\n477SOL H 0.227 0.127 0.548\n477SOL H 0.204 0.085 0.564\n478SOL O 0.699 0.578 0.561\n478SOL H 0.735 0.595 0.563\n478SOL H 0.674 0.608 0.558\n479SOL O 0.305 0.383 0.964\n479SOL H 0.341 0.373 0.970\n479SOL H 0.312 0.412 0.954\n480SOL O 0.177 0.240 0.992\n480SOL H 0.212 0.233 0.984\n480SOL H 0.171 0.207 1.000\n481SOL O 0.593 0.800 0.939\n481SOL H 0.601 0.827 0.928\n481SOL H 0.616 0.768 0.936\n482SOL O 0.652 0.169 0.613\n482SOL H 0.633 0.147 0.624\n482SOL H 0.688 0.178 0.619\n483SOL O 0.974 0.343 0.708\n483SOL H 0.943 0.331 0.716\n483SOL H 0.959 0.367 0.697\n484SOL O 0.755 0.927 0.834\n484SOL H 0.736 0.891 0.831\n484SOL H 0.787 0.925 0.825\n485SOL O 0.807 0.667 0.427\n485SOL H 0.831 0.699 0.425\n485SOL H 0.829 0.636 0.422\n486SOL O 0.139 0.527 0.896\n486SOL H 0.123 0.500 0.886\n486SOL H 0.172 0.540 0.890\n487SOL O 0.445 0.542 0.869\n487SOL H 0.433 0.517 0.858\n487SOL H 0.411 0.557 0.875\n488SOL O 0.554 0.580 0.661\n488SOL H 0.515 0.582 0.664\n488SOL H 0.563 0.541 0.660\n489SOL O 0.165 0.038 0.516\n489SOL H 0.197 0.045 0.525\n489SOL H 0.175 0.005 0.508\n490SOL O 0.303 0.019 0.810\n490SOL H 0.275 0.000 0.819\n490SOL H 0.303 -0.001 0.797\n491SOL O 1.000 0.990 0.781\n491SOL H 0.991 0.955 0.788\n491SOL H 0.967 1.000 0.773\n492SOL O 0.238 0.426 0.698\n492SOL H 0.264 0.401 0.706\n492SOL H 0.215 0.441 0.709\n493SOL O 0.101 0.806 0.427\n493SOL H 0.088 0.800 0.412\n493SOL H 0.127 0.776 0.430\n494SOL O 0.195 0.783 0.710\n494SOL H 0.197 0.746 0.717\n494SOL H 0.203 0.807 0.723\n495SOL O 0.115 0.712 0.856\n495SOL H 0.094 0.733 0.866\n495SOL H 0.107 0.674 0.859\n496SOL O 0.669 0.648 0.604\n496SOL H 0.630 0.653 0.604\n496SOL H 0.676 0.610 0.608\n497SOL O 0.195 0.502 0.557\n497SOL H 0.206 0.464 0.560\n497SOL H 0.155 0.499 0.555\n498SOL O 0.930 0.189 0.751\n498SOL H 0.909 0.194 0.737\n498SOL H 0.947 0.153 0.749\n499SOL O 0.102 0.681 0.906\n499SOL H 0.132 0.698 0.915\n499SOL H 0.070 0.702 0.911\n500SOL O 0.258 0.702 0.422\n500SOL H 0.285 0.680 0.430\n500SOL H 0.275 0.739 0.421\n501SOL O 0.297 0.987 0.416\n501SOL H 0.301 0.994 0.432\n501SOL H 0.331 0.968 0.412\n502SOL O 0.276 0.587 0.733\n502SOL H 0.269 0.596 0.748\n502SOL H 0.239 0.580 0.727\n503SOL O 0.520 0.101 0.605\n503SOL H 0.560 0.111 0.606\n503SOL H 0.501 0.136 0.603\n504SOL O 0.397 0.398 0.702\n504SOL H 0.363 0.376 0.700\n504SOL H 0.389 0.432 0.694\n505SOL O 0.349 0.547 0.848\n505SOL H 0.313 0.528 0.848\n505SOL H 0.341 0.584 0.853\n506SOL O 0.424 0.648 0.736\n506SOL H 0.404 0.626 0.747\n506SOL H 0.450 0.623 0.730\n507SOL O 0.952 0.191 0.632\n507SOL H 0.961 0.221 0.643\n507SOL H 0.914 0.181 0.636\n508SOL O 0.463 0.469 0.514\n508SOL H 0.469 0.476 0.530\n508SOL H 0.457 0.506 0.508\n509SOL O 0.303 0.344 0.647\n509SOL H 0.280 0.328 0.636\n509SOL H 0.277 0.357 0.658\n510SOL O 0.892 0.021 0.633\n510SOL H 0.862 -0.000 0.639\n510SOL H 0.925 0.000 0.635\n511SOL O 0.760 0.591 0.438\n511SOL H 0.731 0.594 0.427\n511SOL H 0.756 0.553 0.444\n512SOL O 0.990 0.640 0.519\n512SOL H 1.000 0.620 0.506\n512SOL H 0.957 0.622 0.524\n513SOL O 0.873 0.266 0.781\n513SOL H 0.858 0.229 0.782\n513SOL H 0.841 0.290 0.781\n514SOL O 0.795 0.961 0.932\n514SOL H 0.801 0.981 0.919\n514SOL H 0.795 0.923 0.928\n515SOL O 0.721 0.267 0.705\n515SOL H 0.693 0.271 0.694\n515SOL H 0.742 0.301 0.704\n516SOL O 0.860 0.863 0.718\n516SOL H 0.881 0.843 0.729\n516SOL H 0.857 0.901 0.723\n517SOL O 0.804 0.283 0.894\n517SOL H 0.824 0.253 0.902\n517SOL H 0.780 0.299 0.905\n518SOL O 0.749 0.719 0.646\n518SOL H 0.754 0.683 0.652\n518SOL H 0.784 0.726 0.639\n519SOL O 0.879 0.771 0.493\n519SOL H 0.898 0.791 0.481\n519SOL H 0.897 0.735 0.493\n520SOL O 0.181 0.443 0.803\n520SOL H 0.165 0.472 0.794\n520SOL H 0.155 0.413 0.802\n521SOL O 0.239 0.981 0.863\n521SOL H 0.267 0.999 0.872\n521SOL H 0.259 0.956 0.854\n522SOL O 0.566 0.187 0.524\n522SOL H 0.561 0.211 0.536\n522SOL H 0.605 0.177 0.524\n523SOL O 0.132 0.952 0.852\n523SOL H 0.161 0.931 0.859\n523SOL H 0.139 0.990 0.857\n524SOL O 0.079 0.085 0.869\n524SOL H 0.090 0.121 0.875\n524SOL H 0.039 0.087 0.869\n525SOL O 0.565 0.530 0.502\n525SOL H 0.536 0.536 0.491\n525SOL H 0.588 0.563 0.501\n526SOL O 0.689 0.831 0.495\n526SOL H 0.696 0.866 0.487\n526SOL H 0.685 0.841 0.510\n527SOL O 0.893 0.286 0.962\n527SOL H 0.885 0.309 0.975\n527SOL H 0.868 0.256 0.964\n528SOL O 0.385 0.146 0.873\n528SOL H 0.406 0.168 0.883\n528SOL H 0.351 0.137 0.880\n529SOL O 0.027 0.226 0.515\n529SOL H 0.045 0.189 0.515\n529SOL H 0.048 0.248 0.525\n530SOL O 0.771 0.409 0.796\n530SOL H 0.779 0.376 0.804\n530SOL H 0.800 0.435 0.800\n531SOL O 0.856 0.149 0.414\n531SOL H 0.830 0.166 0.424\n531SOL H 0.844 0.112 0.412\n532SOL O 0.628 0.950 0.475\n532SOL H 0.603 0.926 0.483\n532SOL H 0.606 0.983 0.471\n533SOL O 0.206 0.615 0.582\n533SOL H 0.213 0.653 0.577\n533SOL H 0.243 0.599 0.584\n534SOL O 0.277 0.693 0.603\n534SOL H 0.296 0.672 0.592\n534SOL H 0.276 0.731 0.597\n535SOL O 0.646 0.032 0.807\n535SOL H 0.666 0.021 0.820\n535SOL H 0.624 0.000 0.803\n536SOL O 0.223 0.376 0.550\n536SOL H 0.247 0.397 0.540\n536SOL H 0.230 0.338 0.546\n537SOL O 0.840 0.228 0.536\n537SOL H 0.836 0.268 0.534\n537SOL H 0.808 0.217 0.544\n538SOL O 0.220 0.082 0.838\n538SOL H 0.238 0.108 0.828\n538SOL H 0.204 0.055 0.829\n539SOL O 0.491 0.457 0.874\n539SOL H 0.488 0.429 0.885\n539SOL H 0.528 0.473 0.876\n540SOL O 0.466 0.322 0.816\n540SOL H 0.432 0.304 0.812\n540SOL H 0.494 0.307 0.807\n541SOL O 0.957 0.962 0.674\n541SOL H 0.935 0.944 0.663\n541SOL H 0.950 1.001 0.672\n542SOL O 0.929 0.542 0.601\n542SOL H 0.895 0.562 0.598\n542SOL H 0.955 0.555 0.590\n543SOL O 0.221 0.880 0.472\n543SOL H 0.257 0.865 0.476\n543SOL H 0.199 0.877 0.485\n544SOL O 0.270 0.432 0.592\n544SOL H 0.271 0.413 0.607\n544SOL H 0.274 0.471 0.596\n545SOL O 0.514 0.455 0.638\n545SOL H 0.545 0.448 0.629\n545SOL H 0.485 0.429 0.634\n546SOL O 0.336 0.800 0.944\n546SOL H 0.344 0.829 0.933\n546SOL H 0.349 0.766 0.937\n547SOL O 0.906 0.571 0.752\n547SOL H 0.899 0.600 0.762\n547SOL H 0.877 0.544 0.755\n548SOL O 0.369 0.387 0.907\n548SOL H 0.392 0.393 0.920\n548SOL H 0.332 0.396 0.912\n549SOL O 0.294 0.605 0.817\n549SOL H 0.272 0.598 0.830\n549SOL H 0.325 0.627 0.822\n550SOL O 0.756 0.257 0.788\n550SOL H 0.731 0.244 0.776\n550SOL H 0.771 0.223 0.794\n551SOL O 0.463 0.728 0.940\n551SOL H 0.470 0.766 0.945\n551SOL H 0.427 0.719 0.946\n552SOL O 0.961 0.167 0.841\n552SOL H 0.933 0.179 0.831\n552SOL H 0.993 0.190 0.838\n1MOL C 0.763 0.357 0.544\n1MOL H 0.805 0.368 0.537\n1MOL H 0.727 0.373 0.534\n1MOL H 0.759 0.309 0.546\n1MOL H 0.760 0.377 0.561\n2MOL C 0.443 0.415 0.782\n2MOL H 0.477 0.435 0.771\n2MOL H 0.451 0.425 0.800\n2MOL H 0.400 0.433 0.777\n2MOL H 0.443 0.368 0.780\n3MOL C 0.459 0.113 0.781\n3MOL H 0.427 0.139 0.791\n3MOL H 0.436 0.076 0.773\n3MOL H 0.479 0.140 0.768\n3MOL H 0.492 0.096 0.793\n4MOL C 0.309 0.495 0.801\n4MOL H 0.316 0.464 0.787\n4MOL H 0.285 0.533 0.794\n4MOL H 0.283 0.474 0.815\n4MOL H 0.351 0.509 0.808\n5MOL C 0.851 0.290 0.577\n5MOL H 0.835 0.324 0.566\n5MOL H 0.814 0.267 0.586\n5MOL H 0.876 0.258 0.567\n5MOL H 0.879 0.309 0.591\n6MOL C 0.439 0.298 0.691\n6MOL H 0.454 0.341 0.685\n6MOL H 0.404 0.282 0.680\n6MOL H 0.475 0.267 0.691\n6MOL H 0.422 0.302 0.708\n7MOL C 0.266 0.206 0.900\n7MOL H 0.235 0.170 0.900\n7MOL H 0.286 0.211 0.883\n7MOL H 0.243 0.247 0.905\n7MOL H 0.299 0.197 0.913\n8MOL C 0.505 0.365 0.428\n8MOL H 0.482 0.404 0.434\n8MOL H 0.551 0.375 0.426\n8MOL H 0.501 0.330 0.441\n8MOL H 0.486 0.350 0.412\n9MOL C 0.709 0.505 0.487\n9MOL H 0.697 0.533 0.473\n9MOL H 0.725 0.532 0.501\n9MOL H 0.671 0.480 0.493\n9MOL H 0.744 0.475 0.481\n10MOL C 0.650 0.456 0.635\n10MOL H 0.625 0.498 0.634\n10MOL H 0.697 0.465 0.632\n10MOL H 0.643 0.435 0.651\n10MOL H 0.635 0.428 0.621\n11MOL C 0.465 0.551 0.422\n11MOL H 0.459 0.536 0.440\n11MOL H 0.431 0.531 0.411\n11MOL H 0.462 0.598 0.421\n11MOL H 0.509 0.537 0.416\n12MOL C 0.538 0.412 0.928\n12MOL H 0.551 0.454 0.921\n12MOL H 0.525 0.417 0.946\n12MOL H 0.500 0.395 0.919\n12MOL H 0.574 0.380 0.927\n13MOL C 0.235 0.254 0.507\n13MOL H 0.272 0.234 0.499\n13MOL H 0.245 0.263 0.525\n13MOL H 0.225 0.296 0.499\n13MOL H 0.197 0.225 0.506\n14MOL C 0.781 0.706 0.769\n14MOL H 0.824 0.686 0.774\n14MOL H 0.769 0.737 0.783\n14MOL H 0.786 0.730 0.753\n14MOL H 0.748 0.672 0.768\n15MOL C 0.701 0.757 0.460\n15MOL H 0.654 0.760 0.457\n15MOL H 0.710 0.750 0.478\n15MOL H 0.719 0.720 0.450\n15MOL H 0.723 0.797 0.454\n16MOL C 0.214 0.491 0.934\n16MOL H 0.220 0.539 0.934\n16MOL H 0.187 0.480 0.949\n16MOL H 0.192 0.477 0.919\n16MOL H 0.256 0.470 0.936\n17MOL C 0.681 0.647 0.828\n17MOL H 0.667 0.669 0.812\n17MOL H 0.643 0.627 0.836\n17MOL H 0.713 0.613 0.823\n17MOL H 0.701 0.679 0.839\n18MOL C 0.106 0.674 0.785\n18MOL H 0.145 0.683 0.795\n18MOL H 0.094 0.627 0.787\n18MOL H 0.069 0.701 0.791\n18MOL H 0.115 0.684 0.767\n19MOL C 0.296 0.571 0.523\n19MOL H 0.309 0.616 0.523\n19MOL H 0.280 0.559 0.540\n19MOL H 0.261 0.564 0.510\n19MOL H 0.334 0.543 0.519\n20MOL C 0.488 0.757 0.715\n20MOL H 0.459 0.725 0.706\n20MOL H 0.527 0.734 0.721\n20MOL H 0.501 0.791 0.703\n20MOL H 0.465 0.776 0.729\n21MOL C 0.563 0.161 0.873\n21MOL H 0.519 0.163 0.880\n21MOL H 0.563 0.173 0.854\n21MOL H 0.592 0.191 0.882\n21MOL H 0.581 0.116 0.874\n22MOL C 0.254 0.969 0.621\n22MOL H 0.287 1.001 0.626\n22MOL H 0.263 0.954 0.603\n22MOL H 0.256 0.931 0.633\n22MOL H 0.211 0.989 0.622\n23MOL C 0.643 0.223 0.659\n23MOL H 0.691 0.220 0.663\n23MOL H 0.625 0.259 0.669\n23MOL H 0.636 0.231 0.641\n23MOL H 0.622 0.181 0.664\n24MOL C 0.928 0.664 0.715\n24MOL H 0.887 0.646 0.708\n24MOL H 0.920 0.704 0.724\n24MOL H 0.948 0.632 0.727\n24MOL H 0.958 0.672 0.701\n25MOL C 0.061 0.588 0.607\n25MOL H 0.086 0.576 0.622\n25MOL H 0.042 0.549 0.599\n25MOL H 0.026 0.619 0.612\n25MOL H 0.090 0.609 0.595\n26MOL C 0.904 0.819 0.588\n26MOL H 0.870 0.852 0.592\n26MOL H 0.915 0.822 0.569\n26MOL H 0.943 0.827 0.598\n26MOL H 0.887 0.775 0.592\n27MOL C 0.079 0.845 0.632\n27MOL H 0.113 0.823 0.623\n27MOL H 0.073 0.823 0.649\n27MOL H 0.091 0.891 0.635\n27MOL H 0.037 0.844 0.623\n28MOL C 0.976 0.699 0.843\n28MOL H 1.000 0.720 0.829\n28MOL H 0.972 0.652 0.839\n28MOL H 1.000 0.705 0.859\n28MOL H 0.933 0.718 0.844\n29MOL C 0.620 0.389 0.751\n29MOL H 0.621 0.436 0.746\n29MOL H 0.623 0.386 0.770\n29MOL H 0.579 0.370 0.745\n29MOL H 0.657 0.367 0.743\n30MOL C 0.482 0.018 0.487\n30MOL H 0.483 0.066 0.488\n30MOL H 0.492 0.000 0.504\n30MOL H 0.437 0.004 0.482\n30MOL H 0.514 0.003 0.474\n31MOL C 0.189 0.599 0.429\n31MOL H 0.205 0.638 0.438\n31MOL H 0.177 0.566 0.442\n31MOL H 0.152 0.611 0.418\n31MOL H 0.225 0.582 0.418\n32MOL C 0.657 0.037 0.929\n32MOL H 0.681 0.078 0.927\n32MOL H 0.635 0.037 0.946\n32MOL H 0.625 0.032 0.916\n32MOL H 0.688 0.000 0.929\n33MOL C 0.463 0.666 0.890\n33MOL H 0.421 0.644 0.889\n33MOL H 0.466 0.687 0.907\n33MOL H 0.499 0.635 0.888\n33MOL H 0.467 0.700 0.877\n34MOL C 0.627 0.427 0.867\n34MOL H 0.642 0.461 0.878\n34MOL H 0.664 0.407 0.858\n34MOL H 0.605 0.392 0.877\n34MOL H 0.596 0.445 0.854\n35MOL C 0.942 0.074 0.790\n35MOL H 0.974 0.066 0.776\n35MOL H 0.952 0.115 0.798\n35MOL H 0.898 0.075 0.783\n35MOL H 0.945 0.038 0.802\n36MOL C 0.205 0.429 0.462\n36MOL H 0.200 0.388 0.472\n36MOL H 0.243 0.424 0.450\n36MOL H 0.165 0.437 0.452\n36MOL H 0.214 0.465 0.474\n37MOL C 0.772 0.165 0.491\n37MOL H 0.816 0.146 0.493\n37MOL H 0.739 0.134 0.497\n37MOL H 0.770 0.205 0.502\n37MOL H 0.764 0.177 0.473\n38MOL C 0.775 0.947 0.538\n38MOL H 0.773 0.900 0.536\n38MOL H 0.803 0.958 0.553\n38MOL H 0.792 0.967 0.523\n38MOL H 0.731 0.964 0.542\n39MOL C 0.261 0.694 0.945\n39MOL H 0.231 0.699 0.960\n39MOL H 0.256 0.731 0.933\n39MOL H 0.306 0.693 0.952\n39MOL H 0.251 0.653 0.936\n40MOL C 0.734 0.682 0.705\n40MOL H 0.726 0.727 0.710\n40MOL H 0.717 0.677 0.687\n40MOL H 0.781 0.673 0.705\n40MOL H 0.712 0.651 0.716\n41MOL C 0.341 0.426 0.737\n41MOL H 0.330 0.380 0.734\n41MOL H 0.331 0.451 0.721\n41MOL H 0.315 0.443 0.751\n41MOL H 0.388 0.430 0.741\n42MOL C 0.304 0.736 0.716\n42MOL H 0.285 0.714 0.731\n42MOL H 0.287 0.715 0.701\n42MOL H 0.293 0.782 0.716\n42MOL H 0.352 0.731 0.717\n43MOL C 0.876 0.436 0.984\n43MOL H 0.880 0.480 0.977\n43MOL H 0.912 0.428 0.996\n43MOL H 0.834 0.432 0.992\n43MOL H 0.879 0.405 0.969\n44MOL C 0.564 0.314 0.647\n44MOL H 0.538 0.346 0.637\n44MOL H 0.611 0.322 0.644\n44MOL H 0.555 0.318 0.665\n44MOL H 0.552 0.270 0.641\n45MOL C 0.645 0.638 0.462\n45MOL H 0.685 0.653 0.471\n45MOL H 0.628 0.672 0.451\n45MOL H 0.656 0.599 0.452\n45MOL H 0.613 0.626 0.475\n46MOL C 0.580 0.846 0.872\n46MOL H 0.607 0.821 0.860\n46MOL H 0.536 0.827 0.873\n46MOL H 0.600 0.844 0.889\n46MOL H 0.577 0.891 0.866\n47MOL C 0.738 0.820 0.881\n47MOL H 0.776 0.793 0.886\n47MOL H 0.731 0.852 0.895\n47MOL H 0.746 0.843 0.865\n47MOL H 0.699 0.792 0.880\n48MOL C 0.511 0.559 0.918\n48MOL H 0.466 0.567 0.912\n48MOL H 0.511 0.532 0.933\n48MOL H 0.534 0.535 0.904\n48MOL H 0.533 0.600 0.921\n49MOL C 0.298 0.065 0.658\n49MOL H 0.336 0.061 0.669\n49MOL H 0.310 0.089 0.642\n49MOL H 0.283 0.021 0.653\n49MOL H 0.263 0.088 0.667\n50MOL C 0.030 0.735 0.977\n50MOL H 0.015 0.708 0.963\n50MOL H 0.075 0.750 0.974\n50MOL H 0.002 0.774 0.979\n50MOL H 0.029 0.710 0.993\n51MOL C 0.431 0.797 0.863\n51MOL H 0.400 0.834 0.863\n51MOL H 0.466 0.806 0.850\n51MOL H 0.451 0.791 0.880\n51MOL H 0.408 0.757 0.858\n52MOL C 0.097 0.222 0.683\n52MOL H 0.113 0.186 0.672\n52MOL H 0.122 0.225 0.699\n52MOL H 0.101 0.263 0.673\n52MOL H 0.051 0.215 0.687\n53MOL C 0.150 0.879 0.785\n53MOL H 0.171 0.883 0.768\n53MOL H 0.169 0.911 0.797\n53MOL H 0.103 0.888 0.783\n53MOL H 0.156 0.835 0.792\n54MOL C 0.304 0.331 0.578\n54MOL H 0.307 0.304 0.562\n54MOL H 0.320 0.306 0.592\n54MOL H 0.330 0.371 0.576\n54MOL H 0.258 0.343 0.581\n55MOL C 0.315 0.573 0.895\n55MOL H 0.299 0.609 0.884\n55MOL H 0.282 0.563 0.908\n55MOL H 0.323 0.535 0.885\n55MOL H 0.356 0.587 0.904\n56MOL C 0.911 0.347 0.499\n56MOL H 0.885 0.320 0.511\n56MOL H 0.957 0.333 0.501\n56MOL H 0.896 0.340 0.482\n56MOL H 0.907 0.393 0.504\n57MOL C 0.029 0.795 0.462\n57MOL H 0.004 0.809 0.447\n57MOL H 0.000 0.771 0.473\n57MOL H 0.065 0.767 0.456\n57MOL H 0.046 0.834 0.471\n58MOL C 0.379 0.732 0.432\n58MOL H 0.354 0.723 0.448\n58MOL H 0.426 0.735 0.437\n58MOL H 0.373 0.697 0.420\n58MOL H 0.365 0.774 0.425\n59MOL C 0.401 0.120 0.721\n59MOL H 0.434 0.107 0.734\n59MOL H 0.399 0.087 0.707\n59MOL H 0.358 0.125 0.729\n59MOL H 0.414 0.163 0.713\n60MOL C 0.621 0.534 0.724\n60MOL H 0.655 0.501 0.724\n60MOL H 0.608 0.543 0.706\n60MOL H 0.582 0.517 0.733\n60MOL H 0.636 0.574 0.732\n61MOL C 0.885 0.410 0.857\n61MOL H 0.841 0.390 0.857\n61MOL H 0.917 0.376 0.852\n61MOL H 0.885 0.445 0.844\n61MOL H 0.895 0.429 0.874\n62MOL C 0.307 0.758 0.658\n62MOL H 0.309 0.713 0.665\n62MOL H 0.276 0.760 0.644\n62MOL H 0.293 0.787 0.672\n62MOL H 0.351 0.771 0.652\n63MOL C 0.481 0.927 0.580\n63MOL H 0.445 0.917 0.569\n63MOL H 0.523 0.912 0.573\n63MOL H 0.483 0.975 0.583\n63MOL H 0.475 0.905 0.597\n64MOL C 0.772 0.297 0.834\n64MOL H 0.783 0.271 0.849\n64MOL H 0.732 0.279 0.826\n64MOL H 0.765 0.343 0.839\n64MOL H 0.808 0.295 0.822\n65MOL C 0.363 0.201 0.619\n65MOL H 0.365 0.249 0.619\n65MOL H 0.343 0.184 0.604\n65MOL H 0.407 0.184 0.621\n65MOL H 0.337 0.187 0.634\n66MOL C 0.776 0.865 0.760\n66MOL H 0.764 0.819 0.760\n66MOL H 0.749 0.888 0.747\n66MOL H 0.769 0.884 0.777\n66MOL H 0.822 0.869 0.755\n67MOL C 0.573 0.422 0.486\n67MOL H 0.531 0.442 0.483\n67MOL H 0.603 0.432 0.472\n67MOL H 0.592 0.439 0.502\n67MOL H 0.568 0.374 0.488\n68MOL C 0.639 0.264 0.803\n68MOL H 0.597 0.285 0.798\n68MOL H 0.640 0.256 0.821\n68MOL H 0.643 0.223 0.793\n68MOL H 0.675 0.294 0.798\n69MOL C 0.357 0.957 0.705\n69MOL H 0.365 1.000 0.713\n69MOL H 0.344 0.964 0.687\n69MOL H 0.398 0.931 0.705\n69MOL H 0.323 0.934 0.714\n70MOL C 0.337 0.194 0.838\n70MOL H 0.321 0.155 0.846\n70MOL H 0.365 0.219 0.849\n70MOL H 0.363 0.181 0.822\n70MOL H 0.299 0.221 0.832\n71MOL C 0.165 0.381 0.620\n71MOL H 0.146 0.336 0.620\n71MOL H 0.182 0.390 0.603\n71MOL H 0.131 0.413 0.625\n71MOL H 0.200 0.384 0.633\n72MOL C 0.211 0.072 0.724\n72MOL H 0.202 0.067 0.706\n72MOL H 0.188 0.038 0.734\n72MOL H 0.258 0.069 0.727\n72MOL H 0.196 0.116 0.730\n73MOL C 0.706 0.813 0.631\n73MOL H 0.749 0.811 0.639\n73MOL H 0.711 0.811 0.612\n73MOL H 0.679 0.776 0.637\n73MOL H 0.684 0.854 0.636\n74MOL C 0.837 0.228 0.692\n74MOL H 0.851 0.188 0.702\n74MOL H 0.792 0.221 0.685\n74MOL H 0.837 0.265 0.704\n74MOL H 0.867 0.236 0.678\n75MOL C 0.098 0.315 0.484\n75MOL H 0.114 0.343 0.471\n75MOL H 0.051 0.307 0.482\n75MOL H 0.120 0.272 0.483\n75MOL H 0.107 0.334 0.501\n76MOL C 0.509 0.656 0.469\n76MOL H 0.474 0.685 0.463\n76MOL H 0.490 0.619 0.478\n76MOL H 0.533 0.640 0.454\n76MOL H 0.539 0.680 0.480\n77MOL C 0.971 0.048 0.579\n77MOL H 0.993 0.052 0.595\n77MOL H 0.955 0.003 0.576\n77MOL H 0.933 0.078 0.578\n77MOL H 1.001 0.060 0.565\n78MOL C 0.335 0.046 0.601\n78MOL H 0.362 0.083 0.608\n78MOL H 0.354 0.005 0.607\n78MOL H 0.291 0.051 0.608\n78MOL H 0.334 0.048 0.582\n79MOL C 0.977 0.266 0.584\n79MOL H 1.004 0.288 0.597\n79MOL H 1.004 0.235 0.574\n79MOL H 0.941 0.242 0.592\n79MOL H 0.958 0.298 0.572\n80MOL C 0.375 0.862 0.507\n80MOL H 0.368 0.888 0.523\n80MOL H 0.334 0.857 0.497\n80MOL H 0.407 0.884 0.496\n80MOL H 0.392 0.819 0.512\n81MOL C 0.524 0.973 0.926\n81MOL H 0.480 0.956 0.929\n81MOL H 0.524 1.000 0.910\n81MOL H 0.537 1.000 0.941\n81MOL H 0.555 0.936 0.924\n82MOL C 0.053 0.387 0.588\n82MOL H 0.080 0.348 0.585\n82MOL H 0.032 0.384 0.605\n82MOL H 0.081 0.427 0.587\n82MOL H 0.020 0.389 0.575\n83MOL C 0.864 0.631 0.840\n83MOL H 0.817 0.625 0.836\n83MOL H 0.886 0.652 0.826\n83MOL H 0.868 0.659 0.855\n83MOL H 0.884 0.588 0.843\n84MOL C 0.088 0.108 0.809\n84MOL H 0.084 0.151 0.800\n84MOL H 0.059 0.106 0.823\n84MOL H 0.075 0.074 0.796\n84MOL H 0.133 0.101 0.814\n85MOL C 0.387 0.057 0.521\n85MOL H 0.352 0.060 0.533\n85MOL H 0.429 0.064 0.529\n85MOL H 0.386 0.014 0.512\n85MOL H 0.380 0.091 0.507\n86MOL C 0.818 0.026 0.502\n86MOL H 0.791 0.000 0.491\n86MOL H 0.829 -0.000 0.517\n86MOL H 0.858 0.040 0.494\n86MOL H 0.793 0.065 0.508\n87MOL C 0.432 0.741 0.620\n87MOL H 0.453 0.711 0.632\n87MOL H 0.447 0.731 0.603\n87MOL H 0.384 0.736 0.621\n87MOL H 0.444 0.786 0.625\n88MOL C 0.072 0.445 0.933\n88MOL H 0.029 0.462 0.936\n88MOL H 0.085 0.417 0.948\n88MOL H 0.073 0.419 0.917\n88MOL H 0.104 0.481 0.931\n89MOL C 0.662 0.926 0.668\n89MOL H 0.688 0.951 0.680\n89MOL H 0.674 0.938 0.651\n89MOL H 0.669 0.879 0.671\n89MOL H 0.616 0.935 0.671\n90MOL C 0.613 0.696 0.652\n90MOL H 0.600 0.699 0.634\n90MOL H 0.613 0.741 0.659\n90MOL H 0.581 0.669 0.661\n90MOL H 0.656 0.677 0.654\n91MOL C 0.188 0.408 0.753\n91MOL H 0.210 0.426 0.768\n91MOL H 0.222 0.395 0.741\n91MOL H 0.162 0.370 0.758\n91MOL H 0.160 0.442 0.745\n92MOL C 0.505 0.829 0.500\n92MOL H 0.471 0.820 0.488\n92MOL H 0.503 0.796 0.514\n92MOL H 0.548 0.828 0.492\n92MOL H 0.498 0.873 0.508\n93MOL C 0.795 0.975 0.874\n93MOL H 0.785 0.932 0.881\n93MOL H 0.827 0.997 0.884\n93MOL H 0.813 0.969 0.856\n93MOL H 0.755 1.001 0.873\n94MOL C 0.141 0.798 0.521\n94MOL H 0.173 0.766 0.528\n94MOL H 0.124 0.783 0.505\n94MOL H 0.106 0.803 0.534\n94MOL H 0.163 0.840 0.519\n95MOL C 0.139 0.655 0.546\n95MOL H 0.108 0.662 0.560\n95MOL H 0.116 0.659 0.530\n95MOL H 0.173 0.688 0.547\n95MOL H 0.158 0.611 0.548\n96MOL C 0.821 0.511 0.702\n96MOL H 0.838 0.484 0.716\n96MOL H 0.776 0.523 0.705\n96MOL H 0.848 0.551 0.700\n96MOL H 0.824 0.485 0.686\n';var cell=ChemDoodle.readCIF(Fmol, 1,1,1);Mol3.loadContent([cell.molecule], [cell.unitCell]);Mol3.startAnimation();Mol3.stopAnimation();var d2r=Math.PI/180, mat3 = [];ChemDoodle.lib.mat4.identity(mat3);ChemDoodle.lib.mat4.rotate(mat3,0*d2r,[1,0,0]);ChemDoodle.lib.mat4.rotate(mat3,90*d2r,[0,1,0]);ChemDoodle.lib.mat4.rotate(mat3,0*d2r,[0,0,1]);ChemDoodle.lib.mat4.multiply(Mol3.rotationMatrix, mat3);Mol3.repaint();var $=function(id){return document.getElementById(id)};function setProj3(yesPers){Mol3.specs.projectionPerspective_3D = yesPers;Mol3.setupScene();Mol3.repaint()}function setModel3(model){Mol3.specs.set3DRepresentation(model);Mol3.setupScene();Mol3.repaint()}function setSupercell3(){var cell=ChemDoodle.readCIF(Fmol, $("Mol3x").value, $("Mol3y").value, $("Mol3z").value);Mol3.loadContent([cell.molecule], [cell.unitCell]);Mol3.repaint()}</script><br><span class="meta">视图: <input type="radio" name="group2" onclick="setProj3(true)">投影 <input type="radio" name="group2" onclick="setProj3(false)" checked="">正交<br>模型: <input type="radio" name="model" onclick="setModel3('Ball and Stick')" checked="">球棍 <input type="radio" name="model" onclick="setModel3('van der Waals Spheres')">范德华球 <input type="radio" name="model" onclick="setModel3('Stick')">棍状 <input type="radio" name="model" onclick="setModel3('Wireframe')">线框 <input type="radio" name="model" onclick="setModel3('Line')">线型 <input type="checkbox" onclick="Mol3.specs.atoms_displayLabels_3D=this.checked;Mol3.repaint()">名称<br>超晶胞: X <input type="text" style="width:20px;" id="Mol3x" value="1"> Y <input type="text" style="width:20px;" id="Mol3y" value="1"> Z <input type="text" style="width:20px;" id="Mol3z" value="1"> <input type="button" value="创建" onclick="setSupercell3()"><br>左键: 转动 滚轮: 缩放 双击: 自动旋转开关 Alt+左键: 移动</span><br><figurecaption>Fig.3</figurecaption></figure>依次进行EM, NVT, NPT和80 ns的MD模拟, 水合物的最终结构如下图.
计算得到水合物生成过程中F3/F4的变化如下图.
可以看出, 代码的计算结果与文献一致, 当水合物生成以后, 对应的F3值接近0.01, F4值接近0.7.
$\q$, $\f$, F3和F4的分布如下
我们使用Fortran处理gro文件时, 需要先进行一步转换步骤, 当轨迹很大时, 得到的gro文件会变得更大. 要想计算速度更快一些, 那么我们可以考虑直接处理xtc
轨迹文件. 为此, 我们需要使用GROMACS提供的XTC库.
可惜这里留给我空白太少了, 只得作罢.
如果你使用C/C++完成了上面的处理程序, 那么你可以更进一步, 将其直接集成到GROMACS的源代码中, 称为GROMACS的一个内置工具, 像GROMACS自带的其他工具一样. 这样你重新编译后就可以直接gmx aop
得到计算结果了.
如果你确定自己的代码没有问题, 那么可以将其发送给GROMACS开发人员, 让他们将其合并到新版本的发布中, 这样会有更多人受益的.
具体流程可以参考GROMACS团簇分析代码.
在前面的定义中, 中心水分子氧原子周围其他水分子氧原子的数目和编号都不是确定的, 因而要计算的角度的数目也是不确定的. 这样我们就必须一帧一帧的处理, 需要使用gmx select
动态地选取所需的水分子, 而不能使用固定的编号列表. 如果我们改变一下定义, 计算时只考虑距离中心水分子氧原子最近的4个其他水分子氧原子所组成的角度(实际上, 有些论文中的AOP/F3序参数就是这样定义的, 前面定义中所用的距离截断距离3.5Å, 来自SPC水模型的RDF峰值. 其对应的配位数大致是4), 那计算流程就可以简化一些. 因为我们可以先使用gmx trjorder
, 按照距离对中心氧原子周围的其他水分子进行排序, 这样我们所需要的4个其他氧原子编号就是固定的2 3 4 5
, 从而要计算的角度的编号三联对也使固定的. 这样我们就不需要再使用gmx select
, 也不需要自己写代码生成角度编号三联对了.
如果我们需要计算处于不同盒子位置的氧原子的序参数, 用以表征界面的位置, 那我们需要输出每个中心氧原子的位置.
我们可以将F3, F4输出到相应PDB文件的占据数列和温度因子列, 这样就可以直接使用VMD等可视化程序对每个水分子按照其相应的F3或F4进行着色了. 对PDB进行简单的处理也可以得到每个时刻不同位置水分子的序参数, 因为可以很容易地计算界面位置.
如果你的轨迹是由其他程序生成的, 那你可以先将其转换为.gro
格式的轨迹. 你可以使用VMD进行转换, 也可以自己写脚本或程序进行转换.
平方/未知: 固态水 0.01, 液态水 0.1
平方/平均/错误角度: 固态水 0.0, 液态水 0.1
平方/未平均: 冰 0, 固态水 0.02
平方/平均: 固态水 0.01, 液态水 0.1
平方/未平均: 冰 0, 水合物 0.1, 液态水 0.8, 界面 0.4
平方/平均:
未平方/平均/10倍: 水合物 0.1, 液态水 0.9
平方/平均: 冰 0, 水合物 0.1, 液态水 0.8
未知: 固态水 0.01, 液态水 0.1
平方/平均/错误角度: 冰 0, 液态水 0.1