费米维基

用户工具

站点工具

您的足迹:
atk:生成无定形结构

差别

这里会显示出您选择的修订版和当前版本之间的差别。

到此差别页面的链接

两侧同时换到之前的修订记录前一修订版
后一修订版		前一修订版
atk:生成无定形结构 [2016/10/10 16:29] – [修改脚本] dong.dongatk:生成无定形结构 [2018/03/20 22:23] (当前版本) – liu.jun
行 1: 行 1:
-======生成无定形结构======+======创建无定形材料结构模型======
  
 =====引言===== =====引言=====
  
  
-无定形固体材料在广泛的技术领域里的重要性越来越高。晶体结构通常都能够很好的定义和描述，而要对特定无定形材料进行原子尺度的模拟，建立结构模型是相当困难的一件事。本例示范如何使用 VNL 创建无定形结构，并使用 ATK 在不同精度下进行分子动力学模拟。+无定形固体材料在广泛的技术领域里的重要性越来越高。晶体结构通常都能够很好的定义和描述，而要对特定无定形材料进行原子尺度的模拟，建立结构模型是相当困难的一件事。本例示范如何使用 VNL 创建无定形结构，并使用 QuantumATK 在不同精度下进行分子动力学模拟。
  
 无定形结构，顾名思义就是组分原子、分子没有任何周期性顺序的结构。无定形结构常见于各种材料中，例如玻璃、聚合物。其力学性质、电子性质，通常与晶体区别很大。一个典型的例子是无定形二氧化硅的热胀系数比对应的石英晶体要小一个数量级。 无定形结构，顾名思义就是组分原子、分子没有任何周期性顺序的结构。无定形结构常见于各种材料中，例如玻璃、聚合物。其力学性质、电子性质，通常与晶体区别很大。一个典型的例子是无定形二氧化硅的热胀系数比对应的石英晶体要小一个数量级。
行 16: 行 16:
 本例模拟无定形二氧化硅，也就是熔融的二氧化硅（石英）。这是应用最广泛的无定形材料，比如用于常规或特种玻璃，或者半导体工业。本例还会介绍如何创建晶体区域和无定形区域之间的界面。 本例模拟无定形二氧化硅，也就是熔融的二氧化硅（石英）。这是应用最广泛的无定形材料，比如用于常规或特种玻璃，或者半导体工业。本例还会介绍如何创建晶体区域和无定形区域之间的界面。
  
-如果你不熟悉分子动力学模拟，那么可以从 [[atk:分子动力学基础|【教程：分子动力学基础】]] 开始，这样可以帮助你熟悉 VNL-ATK 的分子动力学功能。+如果你不熟悉分子动力学模拟，那么可以从 [[atk:分子动力学基础|【教程：分子动力学基础】]] 开始，这样可以帮助你熟悉 QuantumATK 的分子动力学功能。
  
 <WRAP center box 100%> <WRAP center box 100%>
行 29: 行 29:
 您应该从 SiO<sub>2</sub> 晶体出发，本例使用方石英晶体结构。之所以选用方石英，是因为方石英的晶体结构很容易转换成正交盒子，这样进行模拟会更方便，石英就没有这么方便。另外，方石英的密度与实验测量的无定形二氧化硅密度很接近，为 2.2g/cm<sup>3</sup>。 您应该从 SiO<sub>2</sub> 晶体出发，本例使用方石英晶体结构。之所以选用方石英，是因为方石英的晶体结构很容易转换成正交盒子，这样进行模拟会更方便，石英就没有这么方便。另外，方石英的密度与实验测量的无定形二氧化硅密度很接近，为 2.2g/cm<sup>3</sup>。
  
-打开 Builder，点击 Add ‣ From Database，搜索 SiO2。选择 cristobalite，点击右下角“➕”将其添加到 Builder 窗口中：+打开 Builder，点击 Add ‣ From Database，搜索 SiO2。选择 cristobalite，点击右下角“＋”将其添加到 Builder 窗口中：
  
 {{ :atk:amorphous01.png?600 |}} {{ :atk:amorphous01.png?600 |}}
行 40: 行 40:
 将结构拖到 **ScriptGenerator**，添加一个 NewCalculator，以及 Optimize ‣ MolecularDynamics 到 Script 面板。双击 NewCalculator，选择 **ATK-Classical** 计算器，使用 Pedone_Fe2_2006 力场((A. Pedone, et al: A New Self-Consistent Empirical Interatomic Potential Model for Oxides, Silicates, and Silica-Based Glasses. J. Phys. Chem. B 110, 11780, 2006))。去掉 Print 和 Save 的 ☑️，点击 OK 结束 NewCalculator 设置。 将结构拖到 **ScriptGenerator**，添加一个 NewCalculator，以及 Optimize ‣ MolecularDynamics 到 Script 面板。双击 NewCalculator，选择 **ATK-Classical** 计算器，使用 Pedone_Fe2_2006 力场((A. Pedone, et al: A New Self-Consistent Empirical Interatomic Potential Model for Oxides, Silicates, and Silica-Based Glasses. J. Phys. Chem. B 110, 11780, 2006))。去掉 Print 和 Save 的 ☑️，点击 OK 结束 NewCalculator 设置。
  
-双击 **MolecularDynamics**，为了熔化晶体，需要将温度设置的非常高，比如 5000K。因为这个初始的模拟主要的目的是为了将原子的排布随机化，因此这一步不需要固定压强。因此 Type 选择 NVT 系综（关于 ATK 分子动力学模拟的系综，更详细的文档参考【[[atk:分子动力学基础]]】）。NVT 系综中，Langevin 是一个很好的选择，系统与热浴紧密耦合，因此不会出现某个原子过热而导致系统不稳定的情况。+双击 **MolecularDynamics**，为了熔化晶体，需要将温度设置的非常高，比如 5000K。因为这个初始的模拟主要的目的是为了将原子的排布随机化，因此这一步不需要固定压强。因此 Type 选择 NVT 系综（关于 QuantumATK 分子动力学模拟的系综，更详细的文档参考【[[atk:分子动力学基础]]】）。NVT 系综中，Langevin 是一个很好的选择，系统与热浴紧密耦合，因此不会出现某个原子过热而导致系统不稳定的情况。
  
 设置 Steps 为 20000，Log interval 为 5000，保存轨迹的文件名为 SiO2_5000K_low_density_traj.nc。在 Initial Velocity 栏，选择Maxwell-Boltzmann，并设置温度为 5000K。 设置 Steps 为 20000，Log interval 为 5000，保存轨迹的文件名为 SiO2_5000K_low_density_traj.nc。在 Initial Velocity 栏，选择Maxwell-Boltzmann，并设置温度为 5000K。
行 257: 行 257:
 最后，这一节在技术上，提供其他产生无定形结构的例子，包括三氧化二铝、二氧化钛、二氧化铪。 最后，这一节在技术上，提供其他产生无定形结构的例子，包括三氧化二铝、二氧化钛、二氧化铪。
 ====Al2O3==== ====Al2O3====
-创建三氧化二铝无定形结构，基本流程可以参考文献<sup>3</sup> 和((G. Gutierrez and B. Johansson: Molecular dynamics study of structural properties of amorphous Al2O3. Phys. Rev. B 65, 104202, 2002))。使用相同的经典力场，可以选择 Matsui_CaMgAlSiO_1994，或者可以使用比较早的不包含该力场的 VNL-ATK 版本，下载该操作向导附录的python文件。如果将其拷贝并粘贴到脚本的最前面，模拟需要的力场就有了。+创建三氧化二铝无定形结构，基本流程可以参考文献<sup>3</sup> 和((G. Gutierrez and B. Johansson: Molecular dynamics study of structural properties of amorphous Al2O3. Phys. Rev. B 65, 104202, 2002))。使用相同的经典力场，可以选择 Matsui_CaMgAlSiO_1994，或者可以使用比较早的不包含该力场的 QuantumATK 版本，下载该操作向导附录的python文件。如果将其拷贝并粘贴到脚本的最前面，模拟需要的力场就有了。
  
 为了设置这个体系，你需要导入 α-Al<sub>2</sub>O<sub>3</sub>，因此需要下载文件 {{ :atk:al2o3_corundum_ams_data.zip |}}。在 **Builder** 里面，点击Add ‣ From Files，选中该文件。为了更方便的控制 MD 模拟，三角晶系盒子需要改为立方盒子。有几种方式。可能的一种：Builders ‣ Surface(Cleave) tool。选择（10-10）米勒指数点击 Next 并选择 3x1 表面晶格。点击 Next 选择 Periodic (bulk-like) ，厚度为 4 层。点击 Finish 的到正交的超胞结构，包含 360 原子。 为了设置这个体系，你需要导入 α-Al<sub>2</sub>O<sub>3</sub>，因此需要下载文件 {{ :atk:al2o3_corundum_ams_data.zip |}}。在 **Builder** 里面，点击Add ‣ From Files，选中该文件。为了更方便的控制 MD 模拟，三角晶系盒子需要改为立方盒子。有几种方式。可能的一种：Builders ‣ Surface(Cleave) tool。选择（10-10）米勒指数点击 Next 并选择 3x1 表面晶格。点击 Next 选择 Periodic (bulk-like) ，厚度为 4 层。点击 Finish 的到正交的超胞结构，包含 360 原子。
行 280: 行 280:
 为了人为地降低密度，需要再次将三个晶格常数改为原先的 1.1 倍。 为了人为地降低密度，需要再次将三个晶格常数改为原先的 1.1 倍。
  
-第一个 MolecularDynamics 产生 3000K 的液态结构，密度被降低。基于这个体系，密度需要被调整到目标密度 3.175 g/cm<sup>3</sup>。通常收工计算微观单胞的宏观密度，以及找到晶格矢量的精确比例因子比较麻烦。通过插入如下ATK单元引擎代码到脚本中，第二个 MolecularDynamics 之前，可以很方便的实现（也就是bulk_configuration = md_trajectory.lastImage()这一行之后的内容）：+第一个 MolecularDynamics 产生 3000K 的液态结构，密度被降低。基于这个体系，密度需要被调整到目标密度 3.175 g/cm<sup>3</sup>。通常收工计算微观单胞的宏观密度，以及找到晶格矢量的精确比例因子比较麻烦。通过插入如下QuantumATK单元引擎代码到脚本中，第二个 MolecularDynamics 之前，可以很方便的实现（也就是bulk_configuration = md_trajectory.lastImage()这一行之后的内容）：
  
 <code pyhton> <code pyhton>
行 338: 行 338:
  
 ====HfO2==== ====HfO2====
-对于二氧化铪，可以使用 Wang //et al// 力场((Y. Wang, F. Zahid, J. Wang, H. Guo: Structure and dielectric properties of amorphous high-k oxides: HfO2 , ZrO2 , and their alloys. Phys. Rev. B 85, 224110, 2012))，在 VNL－ATK 中名 为Wang_HfOZr_2012。+对于二氧化铪，可以使用 Wang //et al// 力场((Y. Wang, F. Zahid, J. Wang, H. Guo: Structure and dielectric properties of amorphous high-k oxides: HfO2 , ZrO2 , and their alloys. Phys. Rev. B 85, 224110, 2012))，在 QuantumATK 中名 为Wang_HfOZr_2012。
  
 同样的，可以从 database 中添加 HfO2 结构。通过 Bulk tools ‣ Supercell plugin，点击 Conventional 和 Transform，将 FCC 晶格转为常规单胞。对体系进行 3x3x3 重复，得到合适的超胞尺寸。 同样的，可以从 database 中添加 HfO2 结构。通过 Bulk tools ‣ Supercell plugin，点击 Conventional 和 Transform，将 FCC 晶格转为常规单胞。对体系进行 3x3x3 重复，得到合适的超胞尺寸。
atk/生成无定形结构.1476088165.txt.gz · 最后更改: 2016/10/10 16:29 由 dong.dong

页面工具

© 2014-2022 费米科技（京ICP备14023855号）