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

博文

Materials Studio官方教程:DMol3——计算化学反应的自由能【1】

已有 4835 次阅读 2021-10-21 17:33 |个人分类:科研干货|系统分类:科研笔记

做计算找华算.png

目的:介绍了DMol3预测分子体系热力学性质的方法。

所用模块:Materials Visualizer、DMol3

前提条件:利用DMol3中的离域内坐标对固体进行几何优化(Geometry optimization for solids using delocalized internal coordinates in DMol3)教程

背景

振动分析计算的结果可用于计算重要的热力学性质,如焓(H)、熵(S)、自由能(G)和恒压热容(Cp)与温度之间的函数关系。DMol3总能量在0 K时产生总电子能。各种平移、旋转和振动分量用于计算有限温度下的HSGCp,如下所述。

介绍

在本教程中,将使用DMol3计算简单反应,即1-丁烯异构化形成环丁烷的自由能:

00.jpg

测量了每种碳氢化合物的燃烧热,结果相减得出[1]:

000.png

在同一文件[1]中还给出了每种物质的熵,图片=-9.8 cal/mol/K。这些结果可用于计算T=298.15 K时的自由能,使用:

0000.jpg

因此:

00000.png

为了通过第一性原理计算确定反应焓或自由能,使用如下所示的热力学循环,如果已知T0处的反应焓以及两个温度之间的产物和反应物的热容,则可以计算温度T1处的反应焓:

11.png

得出下列方程:

111.png

使用DMol3进行几何优化得到的结果给出了体系的总电子和离子能量(Etot)。焓值之差由下式给出:

3.png

注意:在某些情况下,必须包括可能得体系体积变化。

体系的内能(U)由电子、振动、平移和旋转分量产生,由以下方程式给出:

12.png

吉布斯自由能差值由下式给出:

123.png

式中,ΔS由以下等式给出:

04.png

一旦确定基态结构,并使用统计力学获得3N - 6(5)简正模的振动频率,DMol3就可以计算上述所有项。

本教程包括如下部分:

  • 开始

  • 准备计算结构模型

  • 优化几何结构并运行振动分析计算

  • 根据DMol3输出文件中包含的数据计算有限温度下的自由能

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


1、开始

首先启动Materials Studio并创建一个新项目。

打开New Project对话框,输入butene作为项目名,单击OK按钮。

新项目将以butene为项目名显示于Project Explorer中。

2、准备计算结构模型

在本教程的这一部分中,将在单独的文档中构建反应物和产物结构模型。第一步是创建两个新的3D原子文档。

单击Standard工具条中的New按钮01.jpg的下拉箭头,从下拉列表中选择3D AtomisticDocument。重复该操作,创建两个新的3D原子文档。

确保第一个3D原子文档为当前文档。从Sketch工具条中选择Sketch Atom工具02.jpg。绘制由四个碳原子组成的链,然后按下ESC键停止绘制。单击其中一个位于末端的C-C键,将其变为双键。

为完成该结构,需要添加氢原子,并将结构中的原子排列为更加合理的初始几何构型。

单击Adjust Hydrogen按钮03.jpg,然后单击Sketch工具条上的Clean按钮04.jpg

右键单击3D Viewer,从弹出的快捷菜单中选择Display Style,打开Display Style对话框。在Atom选项卡中选择Ball and stick选项。

应对文档进行重命名,避免文档之间的混淆。

将新文档重命名为1-butene.xsd

现在构建反应中涉及的另一个分子的结构,即环丁烷。

单击3D Atomistic (2).xsd文件中的3D Viewer,使其成为当前文档。从Sketch工具条中单击Sketch Ring按钮05.jpg的下拉箭头,从下拉列表中选择4 Member。在文档中单击一次以插入环丁烷环。

单击Adjust Hydrogen按钮03.jpg,以及Clean按钮06.jpg

在Display Style对话框中,选择Ball and stick选项,关闭对话框。

应重命名文件。

将第二个新文档重命名为cyclobutane.xsd

单击3D Viewer工具条上的3D Viewer Selection Mode按钮07.jpg

建模过程结束后,应在单独的3D原子文档中看到下面所示的两个分子。

08.jpg

1-丁烯和环丁烷的结构

【系列教程】

Materials Studio官方教程:DMol3——搜索过渡态【1】

Materials Studio官方教程:DMol3——对固体进行几何优化【1】

Materials Studio官方教程:DMol3——能带结构和态密度的计算【1】

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

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

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

MS全媒体推广告.png




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

上一篇:Materials Studio官方教程:DMol3——搜索过渡态【2】
下一篇:Materials Studio官方教程:DMol3——计算化学反应的自由能【2】
收藏 IP: 113.104.248.*| 热度|

0

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

数据加载中...

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

GMT+8, 2024-5-11 18:31

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部