0
  • 聊天消息
  • 系统消息
  • 评论与回复
登录后你可以
  • 下载海量资料
  • 学习在线课程
  • 观看技术视频
  • 写文章/发帖/加入社区
会员中心
创作中心

完善资料让更多小伙伴认识你,还能领取20积分哦,立即完善>

3天内不再提示

AuToFF应用案例——冰融化模拟及TIP4P/ICE水模型使用

鸿之微 来源:鸿之微 2023-06-03 10:51 次阅读

目前很多种水分子模型可以用于预测液态水的物理特性,或确定液态水的(未知)结构。AuToFF支持多种水模型,包括OPC, OPC3, SPC, SPC/E, SPC/Eb, TIP3P, TIPS3P, TIP3P-FB, TIP4P, TIP4P/Ew, TIP4P/2005, TIP4P/ICE, TIP4P/ε。其中,TIP4P/ICE是一种刚性平面的四点水模型,亦即除了1个O和2个H原子外还加入了1个虚原子M(在氧附近),虚原子仅具有负电荷,改善了水分子周围的静电分布。TIP4P/ICE是专门用来描述冰的一系列性质的水模型,对水的相图有较好的还原(1bar下冰-Ih熔点为272.2K),因此本案例将介绍利用TIP4P/ICE水模型进行冰融化的分子动力学模拟

a92584da-01a4-11ee-90ce-dac502259ad0.png

图1 TIP4P构型

力场辅助工具-AuToFF

确定起始构型

选择AuToFF-全原子力场模块,上传冰晶体cif结构文件,水模型力场选择TIP4P/ICE,即可自动转成四点水模型的结构。冰晶体的单胞结构文件如下:

a9344f06-01a4-11ee-90ce-dac502259ad0.png

图2 确定冰晶体结构

选择力场

TIP4P/ICE是专门用来描述冰的一系列性质的水模型,对水的相图有较好的还原(1bar下冰-Ih熔点为272.2K),因此模拟冰融化过程也采用TIP4P/ICE力场模型。

a945eaea-01a4-11ee-90ce-dac502259ad0.png

图3 根据力场选择原子类型

生成拓扑文件

根据力场的选择即可生成拓扑文件(top文件)的相关力场参数(itp文件,包括LJ、键、键角参数,原子电荷)。此外生成拓扑文件可支持多款计算软件,包括:GROMACS、LAMMPS、AMBER、Moltemplate、OpenMM、TINKER、CHARMM。本案例采用GROMACS软件进行分子动力学模拟,用户也可选择LAMMPS进行模拟,其中下载的文件夹中除了力场拓扑文件之外还包含力场参数的文献来源。针对晶体结构支持扩胞操作,此处设置的是5*5*5。

a9680f44-01a4-11ee-90ce-dac502259ad0.png

图4 生成冰晶体拓扑结构

MD模拟

在模拟过程中,模拟步长设为2fs,采用Verlet算法来计算运动方程。模拟体系的三个方向均考虑周期性,是体相的模拟。具体的冰融化过程模拟采用NPT系综,从0 K升温到350 K,使用V-rescale控温,参考温度298.15 K, Berendsen控压, 参考压力为 1.01325 bar 。完整的GROMACS的mdp文件输入如下:

define = integrator = md dt
= 0.002 nsteps
= 1000000 comm-grps = system energygrps = ; nstxout = 0 nstvout = 0 nstfout = 0 nstlog = 500 nstenergy = 500 nstxout-compressed = 1000 compressed-x-grps = system ; annealing = single annealing_npoints = 2 annealing_time = 0 2000 ;ps annealing_temp = 0 350 ; pbc = xyz cutoff-scheme = Verlet coulombtype = PME rcoulomb
= 0.9 vdwtype
= cut-off rvdw
= 0.9 DispCorr
= EnerPres ; Tcoupl = V-rescale tau_t
= 0.2 tc_grps = system ref_t = 298.15 ; Pcoupl
= Berendsen pcoupltype = isotropic tau_p = 0.5 ref_p = 1.01325 compressibility = 4.5e-5 ; gen_vel = no gen_temp = 298.15 gen_seed = -1 ; freezegrps = freezedim = constraints = hbonds

MD结果分析

冰融化过程的轨迹变化通过VMD作图如下,可以清晰的展现出冰融化过程的体系结构变化。

a97ad3c2-01a4-11ee-90ce-dac502259ad0.gif

动画5 冰融化轨迹变化过程

编辑:黄飞

声明:本文内容及配图由入驻作者撰写或者入驻合作网站授权转载。文章观点仅代表作者本人,不代表电子发烧友网立场。文章及其配图仅供工程师学习之用,如有内容侵权或者其他违规问题,请联系本站处理。 举报投诉
  • ICE
    ICE
    +关注

    关注

    0

    文章

    33

    浏览量

    18563
  • TOF
    TOF
    +关注

    关注

    9

    文章

    475

    浏览量

    36223
  • 晶体结构
    +关注

    关注

    0

    文章

    20

    浏览量

    220

原文标题:AuToFF应用案例——冰融化模拟及TIP4P/ICE水模型使用

文章出处:【微信号:hzwtech,微信公众号:鸿之微】欢迎添加关注!文章转载请注明出处。

收藏 人收藏

    评论

    相关推荐

    热管式蓄冷空调的运行特性及分析

    模型"该模型可为热管式蓄系统的设计提供理论依据!关键词&热管’蓄冷’空调’运行特性[hide][/hide]
    发表于 03-18 22:17

    飞机机翼覆融化也能用上石墨烯技术了!

    的科学家创建出一种全新的石墨烯纳米带环氧涂层,在被施加电压后,能通过产生的电热实现覆融化。  在James Tour教授的带领下,研究人员将环氧树脂涂层与石墨烯纳米带相结合。石墨烯纳米带是由单层碳原子
    发表于 01-29 11:16

    ICD4与REAL ICE那个好

    嗨,所以…我还没有看到ICD4和真实的比较。ICD4对PIC微控制器编程特别快,特别是DSPIC?如果你已经有了真正的,购买ICD4的好
    发表于 01-25 10:57

    阻止Real Ice为电路供电

    有没有办法阻止真正的给电路供电呢?我有一个开发板的外部电源,我想用我的电源开关在电路板上循环电源。现在,当我拔掉电路板时,它会被真正的供电,所以我也必须拔掉它。谢谢,John Vickers
    发表于 03-19 06:59

    REAL ICE无法获取设备ID

    ICE在连接到硬件上时找不到我的设备。为了得到设备ID,“真正的”进行了所有的测试,包括回路测试。买了一个全新的真正的,当然还有一个。我发现FPGA和OS有不同的版本。(见附加屏幕截图)。我怀疑
    发表于 04-01 07:53

    Real Ice调试器存在很大问题

    。是否有一个司机可以为PIC32 MZ工作的任何人使用真正的确认?输出是:连接到MPLAB实…当前加载的固件在真正的ICE固件套件版本…01.44.26固件类型……PIC32 MZTAVER电压检测
    发表于 04-18 15:45

    电源管理芯片ICE3A3565P资料推荐

    电源管理芯片ICE3A3565P资料下载内容包括:ICE3A3565P引脚功能ICE3A3565P内部方框图ICE3A3565P典型应用电路
    发表于 03-29 07:43

    电源管理芯片ICE3A5065P相关资料下载

    电源管理芯片ICE3A5065P资料下载内容主要介绍了:ICE3A5065P引脚功能ICE3A5065P内部方框图ICE3A5065P典型应用电路
    发表于 03-30 07:08

    tip122.pdf

    The TIP120, TIP121 and TIP122 are siliconEpitaxial-Base NPN power transistors inmonolithic
    发表于 03-29 11:38 70次下载

    tip122-tip127参数中文资料

    tip122-tip127参数 pdf资料 该TIP120,TIP121和TIP122疏外延基NPN达林顿功率晶体管采用TO-220塑料封装。与互补类型的
    发表于 09-19 16:44 168次下载
    <b class='flag-5'>tip122-tip</b>127参数中文资料

    改进Messinger覆模型导线防临界电流计算

    为了得到不同覆气象条件下导线防临界电流,基于电流防时导线表面膜流动,建立导线表面膜流动模型
    发表于 03-22 10:51 0次下载
    改进Messinger覆<b class='flag-5'>冰</b><b class='flag-5'>模型</b>导线防<b class='flag-5'>冰</b>临界电流计算

    输电线路模拟导线覆监测

    模拟导线覆监测系统采用模拟导线的线路覆监测,其思路是在线路上悬挂1根与线路导线同型号、同材质的模拟导线,以产生与线路相同气象条件下的覆
    发表于 10-28 14:11 740次阅读

    模拟导线覆监测装置-风河智能

    武汉风河智能科技自主研发生产的模拟导线覆监测装置,其思路是在线路杆塔上悬挂一根与线路导线同型号、同材质的模拟导线,以产生与线路相同气象条件下的覆情况,并通过测量
    的头像 发表于 12-29 10:33 1137次阅读
    <b class='flag-5'>模拟</b>导线覆<b class='flag-5'>冰</b>监测装置-风河智能

    新火种AI|小摊牌了!大模型已获国内备案,克隆人发布箭在弦上

    1月4日,小公司宣布,已于去年12月成功获得“小模型”的国内备案。
    的头像 发表于 01-09 10:23 374次阅读

    技嘉RTX 4070 Ti SUPER Eagle OC ICE猎鹰显卡评测

    3月正是开学季,技嘉作为专注板卡大厂,在这个时候带来了4款基于GeForce RTX 40 系列的EAGLE OC ICE猎鹰显卡,分别对应RTX 4070 Ti SUPER、RTX 4070 SUPER、RTX 4060 T
    的头像 发表于 03-06 10:00 2126次阅读
    技嘉RTX 4070 Ti SUPER Eagle OC <b class='flag-5'>ICE</b><b class='flag-5'>冰</b>猎鹰显卡评测