热固性聚合物的玻璃过渡温度

../_images/polymerglasstransitiontemp-0.png

在该Reaxff教程中,福彩3d字迷玻璃过渡温度,Tg从聚合物密度的温度依赖性解释,解释了聚合物密度的温度。

这里提出的系统和工作流程最初在出版物中描述 化学结构对环氧聚合物热力学性能的影响:加速Reaxff模拟和实验的比较, 聚合物159,354-368(2018).

教程包括以下步骤:

笔记

福彩3d字迷是 福彩3d字迷苛刻. 为了最佳性能,建议在福彩3d字迷群集上并行执行。 这可以最好地完成 吉伊的远程工作管理

进口聚合物结构

本教程将使用粘合剂升压方法产生的交联环氧聚合物(另见 债券提升教程 ).

首先下载聚合物并进入AMSInput:

点击 here 下载.xyz文件 detda-exy.xyz.
在AMSInput中导入坐标:
文件→导入坐标
../_images/polymerglasstransitiontemp-1.png

模拟退火

为了准确模拟密度的温度依赖性,我们必须确保聚合物结构尽可能均匀。聚合物结构内的大孔是一个非均匀结构的指示器,将导致非常嘈杂和难以分析结果。

模拟退火 是一个简单的优化策略,旨在在复杂和高维势能表面上产生全局最小能量结构。在退火期间,MD模拟期间的温度逐渐增加,希望该系统能够克服越来越多的能量障碍并有效地探索能量景观。加热相之后是冷却时间,其中系统被冷却回到初始温度。由于不能保证它是最佳的,所得到的结构在能量中显着降低了能量。

首先选择MD设置:

1. 在主面板中,选择 任务→分子动力学
2. 选择力字段 chon2017_weak.ff
../_images/polymerglasstransitiontemp-2.png
3. 点击 更多 旁边 任务:分子动力学 去MD细节
4. 设定 步数970000
5. 设定 样本频率5000
../_images/polymerglasstransitiontemp-3.png

要允许单位单元格缩小或扩展我们进一步添加Barostat:

1. 点击 更多 旁边 Barostat.
2. 选择 Berendsen. 来自 Barostat. dropdown menu.
3. 设定所需的 压力1.0 atm.
4. 设定 阻尼常数500 fs.
../_images/polymerglasstransitiontemp-4.png

模拟退火在恒温器面板中定义

1. 点击 更多 旁边 MD主要选择
2. 点击 更多 旁边 恒温器
3. 选择 恒温器→Berendsen.
4. 设定 阻尼常数100 fs
../_images/polymerglasstransitiontemp-5.png

接下来,我们定义所需的温度程序。 对于从福彩3d字迷中提取密度,我们将增加298.15至598.15的温度,这对于模拟退火仍然很低,但足以为我们的目的而充足。

我们想要实现的温度程序很简单:

  • 在30000步的过程中加热25 k
  • 在10000步的过程中对密度进行样本
  • 重复直到达到598.15 k,然后再次逆转到298.15 k
../_images/polymerglasstransitiontemp-5b.png

该程序的各个步骤可以输入输入字段 温度持久. 在我们的情况下,最容易将它们复制和粘贴到GUI中。

1. 从以下列表中复制各个温度
温度::
298.15 298.15 323.15 323.15 348.15 348.15 373.15 373.15 398.15 398.15 423.15 423.15 448.15 448.15 473.15 473.15 498.15 498.15 523.15 523.15 548.15 548.15 573.15 573.15 598.15 598.15 573.15 573.15 548.15 548.15 523.15 523.15 498.15 498.15 473.15 473.15 448.15 448.15 423.15 423.15 398.15 398.15 373.15 373.15 348.15 348.15 323.15 323.15 298.15 298.15
2. 粘贴他们进入 温度 field
3. 从以下列表中复制各个范围
持续时间::
10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000
4. 粘贴他们进入 持续时间 field
../_images/polymerglasstransitiontemp-6.png

我们现在准备开始模拟退火福彩3d字迷

1. 文件→另存为... 并给予它一个合适的名称(例如“SimulatedAnnealing”)
2. 文件→运行

提取密度与温度型材

从轨迹中提取密度和温度,以便在处理后我们利用Python脚本。

1. Download the script densities.py from here
2. 将其放在与simulated analing.ams InputFile中相同的文件夹中

接下来,您需要打开命令行以将脚本与AMS Python解释器执行。

小费

Windows和Mac用户应该打开一个 来自GUI的预先配置的命令行

在命令行中,脚本可以如下执行

$AMSBIN/plams densities.py -v resultsdir=SimulatedAnnealing.results

福彩3d字迷玻璃化转变温度

密度轮廓可以用任何曲线图绘图软件绘制,例如, gnuplot或libreoffice calc。 将文件导入选择的绘图软件,并绘制两列相互互动。 该图将属于两类之一:

../_images/polymerglasstransitiontemp-7.png

如果初始结构的最终密度和平衡结构显着差异,实验密度(如果可用)远远差,则密度是 不融合 并且应该进行另一个模拟退火运行。 如果最终和初始密度接近,则可以考虑仿真会聚,并且可以进行T的福彩3d字迷g.

福彩3d字迷T.g,通过数据点执行两个线性拟合,如下图所示。 这两条线之间的交叉点标志着玻璃化转变温度。

../_images/polymerglasstransitiontemp-8.png

应选择线性拟合点的亚群,使得最大r2获得价值。 If the R2值太低,即数据过于嘈杂,或者曲线不相交,我们建议尝试以下其中一个 疑难解答:

  • 运行另一个模拟退火模拟
  • 延长温度范围
  • 减少或增加密度中的平均参数.py脚本