materialsstudio案例

《materialsstudio案例》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、MaterialsStudio案例1：Au（111）表面自组装单层膜结构优化目的：用MaterialsStudio（MS）软件对金表面自组装膜的结构进行优化。模块：MinimizerMSDiscover结构优化原理分子的势能一般为键合（键长、键角、二面角、扭转角等）和非键合相互作用（静电作用、范德华作用等）能量项的加和，总势能是各类势能之和，如下式：除了一些简单的分子以外，大多数的势能是分子中一些复杂形势的势能的组合。势能为分子中原子坐标的函数，由原子不同的坐标所得到的势能构成势能面（PotentialEnergySurface，

2、PES）。势能越低，构象越稳定，在系统中出现的机率越大；反之，势能越高，构象越不稳定，在系统中出现的机率越小。通常势能面可得到许多极小值的位置，其中对应于最低能量的点称为全局最小值（GlobalEnergyMinimum），相当于分子最稳定的构象。由势能面求最低极小值的过程称为能量最小化（EnergyMinimum），其所对应的结构为最优化结构（OptimizedStructure），能量最小化过程，亦是结构优化的过程。通过最小化算法进行结构优化时，应避免陷入局部最小值（localminimum），也就是避免仅得到某一构象附近的相

3、对稳定的构象，而力求得到全局最小值，即实现全局优化。分子力学的最小化算法能较快进行能量优化，但它的局限性在于易陷入局部势阱，求得的往往是局部最小值，而要寻求全局最小值只能采用系统搜寻法或分子动力学法。在MaterialsStudio的Discover模块中，能量最小化算法有以下四种：1）最陡下降法（SteepestDescent），为一经典的方法，通过迭代求导，对多变量的非线性目标函数极小化，按能量梯度相反的方向对坐标添加一位移，即能12量函数的负梯度方向是目标函数最陡下降的方向，所以称为最陡下降法。此法计算简单，速度快，但在极小

4、值附近收敛性不够好，造成移动方向正交。最陡下降法适用于优化的最初阶段。2）共轭梯度法（ConjugateGradient），在求导时，目标函数下降方向不是仅选取最陡下降法所采用的能量函数的负梯度方向，而是选取两个共轭梯度方向，即前次迭代时的能量函数负梯度方向与当前迭代时的能量函数负梯度方向的线性组合。此法收敛性较好，但对分子起始结构要求较高，因此常与最陡下降法联合使用，先用最陡下降法优化，再用共轭梯度法优化至收敛。3）牛顿方法（Newton），以二阶导数方法求得极小值。此法的收敛很迅速，也常与最陡下降法联合使用。4）综合法（Sma

5、rtMinimizer），该方法可以混合最陡下降法，共轭梯度法和牛顿法进行结构优化，在MS中是可选择的。点开各种方法后面的More，可设定收敛精度（Convergence），算法（Algorithm）和一维搜索（Linesearch，指每一次迭代中的精度）等。12SmartMinimizer中，只有共轭梯度法和牛顿法才可以选择不同的算法：SmartMinimizer中，牛顿法可以设定最大的原子数，如果体系的原子数大于所设定的值，则计算是会自动地转为前面设定的收敛法（共轭梯度法或最陡下降法），收敛精度会改为共轭梯度法的默认收敛精度值

6、。Au（111）表面磺酸甜菜碱封端封端癸烷基硫醇自组装单层膜MS结构优化过程盒子尺寸：34.609×29.970×91.91874Å水层：20Å，水密度为1g/cm3真空层：40Å（以避免镜像重叠，更好地模拟表面）12分子力学模拟设置分子力场采用COMPASS，截断半径为12.5Ǻ，计算范德华作用和库仑力作用的加和方法采用AtomBased，迭代方法采用SmartMinimizer。计算结果12计算结果：计算步数：27967总势能：-13756.136539kcal/mol非键作用能：-12624.680243kcal/mol结果

7、文件包括（1）Status文件（2）SimulationEnergies文件（3）参数文件（4）能量文件结果文件（1）Status文件（txt文件格式）Operation:minimizeMinimizationMethod/Stage:CG(Fletcher-Reeves).Step27967Maximumderivative:0.252CPUTime4:08:38-13756.14Potential(kcal/mol)---------14204.88Electrostatic11215.29vdWRepulsive-9635

8、.09vdWDispersive0.00HydrogenBond-12624.68Nonbond(Total)12--------872.94Angle387.32Bond-1872.15Torsion0.00Out-of-Plane-54.60R