QHDXZJW的个人博客分享 http://blog.sciencenet.cn/u/QHDXZJW

博文

Materials Studio官方教程:Forcite——编辑力场【1】

已有 3051 次阅读 2022-1-4 10:42 |个人分类:科研干货|系统分类:科研笔记

微信图片_20211101094735.jpg

目的:说明了在Dreiding中编辑力场参数的工作流程,以改进计算和实验晶体结构数据之间的拟合情况。

所用模块:Materials Visualizer、Conformers、Forcite Plus, DMol3可选

前提条件:使用Conformers研究结构-能量之间的关系(Using Conformers to probe geometry-energy relationships)教程

背景

力场是经典模拟计算的核心,因为它体现了结构中每种类型的原子与其周围的原子相互作用的方式。对于体系中的每个原子,将指定一个力场类型,该类型描述原子的局部环境。力场包含各种描述性质的信息,如平衡键长和力场类型对之间的静电相互作用。通常,力场可能是普适的,以描述许多体系,或者是针对特定问题进行参数化。一个普适力场,如Dreiding,对于不同环境中的原子有许多参数。但是,对于特定情况,它的性能永远不会像针对该情况单独参数化的力场一样好。因此,能够编辑力场以便与尚未参数化的体系进行使用是非常有用的。

这在许多情况下都可能发生:可能想要参数化力场来模拟分子在分子筛框架中的吸附等温线,可能需要使用实验晶体数据来提高多晶型预测计算的准确性,或者可能需要更改力场以更好地描述聚合物的密度。在所有这些情况下,都需要编辑基础力场并观察参数更改的效果。

编辑力场是一项复杂的任务,因为许多相互作用是显式或隐式耦合的。例如,旋转扭转角时的能量分布不仅取决于扭矩项,还取决于非键能。因此,获得正确的参数可能需要相当长的时间,并涉及许多不同的计算。最好的方法是进行小的系统性更改,并查看对体系的影响。与所有形式的参数化一样,需要在获得可靠结果和过度拟合消耗资源之间权衡。

介绍

在本教程中,将修改小分子的单个扭转角,以改进优化后模型与晶体结构之间的匹配。最初,将使用Dreiding力场优化爆炸材料三氨基三硝基苯(TATB)的晶体结构。然后,将修改Dreiding力场,以改进晶体结构和使用力场优化的模型结构之间的匹配。将使用Conformers模块来探测扭矩-能量关系,使用DMol3提供从头算数据,并使用Forcite优化晶体结构。

本教程包括如下部分:

  • 开始

  • 优化晶体结构

  • 监测扭转角

  • 分析非键项影响

  • 修改范德华项

注意:为了和本教程中的参数保持一致,可以使用Settings Organizer对话框将工程中所有参数都设置为BIOVIA的默认值。

1、开始

首先启动Materials Studio并创建一个新工程。

打开New Project对话框,输入ff_edit作为工程名,单击OK按钮。

新工程将以ff_edit为工程名显示于Project Explorer中。下一步是导入要优化的爆炸材料TATNBZ的晶体结构。

打开Import Document对话框。导航至Examples\Documents\3D Model\文件夹,双击TATNBZ.xsd文件。将结构重命名为TATNBZ_crystal.xsd

2、优化晶体结构

在将修改后的力场应用于晶体结构以验证这些更改之前,将对孤立的分子执行所有力场编辑工作。应该复制一个孤立的分子。

Project Explorer中,右键单A击工程根目录,从弹出的快捷菜单中选择New | 3D Atomistic

TATNBZ_crystal.xsd,单击其中一个分子中的一个原子,右键并选择Select Fragment。按下CTRL + C键。在新建的原子文档中按下CTRL + V键。在Project Explorer中,将新的原子文档重命名为TATNBZ_molecule.xsd

如果检查晶体结构,会发现硝基在苯环所在的平面上。现在,将使用Forcite模块利用标准Dreiding力场优化结构。

确保TATNBZ_crystal.xsd文件为当前活动文档。单击Modules工具条上的Forcite按钮微信图片_20220104104525.jpg,在下拉列表中选择Calculation,打开Forcite Calculation对话框。

Task修改为Geometry Optimization。在Energy选项卡中,将Forcefield设置为Dreiding,并将Charges更改为Charge using Gasteiger

可以使用更复杂的电荷计算方法,例如使用DMol3应用ESP电荷,但Gasteiger电荷将适用于自定义力场的开发。

单击Run按钮。计算任务完成后,打开优化后的TATNBZ_crystal.xsd文档。

应该看到硝基不再和苯环在处于同一平面上。优化后的结构中硝基的离面性质会影响这些分子在晶体中聚集的方式,因此获得平面硝基非常重要。

硝基和苯环之间有一个单一的扭矩,这将控制基团的平面度。修改控制此扭矩的参数将是本教程的重点内容。

Materials Studio是久负盛名计算模拟软件,问世20余年来,经过不断地迭代优化,使其功能异常强大,极易上手,初学者只需通过简单的参数设置和点击鼠标就能完成DFT计算。其计算可靠性久经考验,备受Nature、Science等顶级期刊认可。

华算科技和Materials Studio官方代理深圳浦华系统联合推出Materials Studio建模、计算、分析课程。课程专为零基础学员设计,沿着理论讲解、模型搭建、性质计算、结果分析层层递进讲解,带你快速入门DFT计算。课程极度注重实操,全程线上直播,提供无限回放,课程群在线答疑。(详情点击下方图片跳转)

识别下方二维码报名,或者联系手机13005427160。

MS全媒体推广告.png






https://blog.sciencenet.cn/blog-2531381-1319518.html

上一篇:Materials Studio官方教程:Forcite——寻找分子在表面上的低能构型【2】
下一篇:Materials Studio官方教程:Forcite——编辑力场【2】
收藏 IP: 113.104.215.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-6-4 03:51

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部