这里会显示出您选择的修订版和当前版本之间的差别。
两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版 | ||
atk:添加_组合_修改经典势函数 [2019/04/29 10:02] – [无定形氧化物的势函数] xie.congwei | atk:添加_组合_修改经典势函数 [2019/04/29 10:26] (当前版本) – [设置系统] xie.congwei | ||
---|---|---|---|
行 185: | 行 185: | ||
===== Tersoff 和 Lennard-Jones 的组合势函数 ===== | ===== Tersoff 和 Lennard-Jones 的组合势函数 ===== | ||
+ | |||
+ | 在本节中,您将学习如何结合 Tersoff 势和 Lennard-Jones(LJ) 势模拟石墨晶胞中的单个锂原子。这种材料组合特别有趣,因为它提供了一个非常简单的锂离子电池中典型阳极的模型。Tersoff 势和 LJ 势的组合通常用于多组分材料,其中各个组分间相互作用较弱,主要是排斥力和范德华力。在本例的石墨中,组分中的强键合相互作用可以由 Tersoff 势或对于该材料的另一种合适的势描述。 | ||
+ | |||
+ | 在本例中,您将采用文献< | ||
==== 设置系统 ==== | ==== 设置系统 ==== | ||
+ | |||
+ | 为设置测试系统,打开 **Builder**,点击 Add {{: | ||
+ | |||
+ | 如下图所示选中一个包含 3 个原子(一个在第二层,两个在第三层)的三角形。 | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | 单击工具栏(如图)中的 {{: | ||
+ | |||
+ | 使用 {{: | ||
+ | |||
+ | 首先,找到元素列表。将列表中最后一个元素的名称从 '' | ||
+ | |||
+ | 如果您现在看一下计算器的部分,将会发现 Tersoff 势的详细定义。这是 //Show defaults// 选项的缘故。通过 '' | ||
+ | |||
+ | <code python> | ||
+ | potentialSet.addParticleType(ParticleType(symbol=' | ||
+ | mass=12.0107*atomic_mass_unit, | ||
+ | charge=None, | ||
+ | # Sigma and epsilon from Ref [3]. | ||
+ | sigma=3.3611*Ang, | ||
+ | epsilon=0.004207*eV)) | ||
+ | potentialSet.addParticleType(ParticleType(symbol=' | ||
+ | mass=6.941*atomic_mass_unit, | ||
+ | charge=None, | ||
+ | # Sigma and epsilon from Ref [3]. | ||
+ | sigma=0.826*Ang, | ||
+ | epsilon=0.271115*eV)) | ||
+ | </ | ||
+ | |||
+ | 对于碳,已添加 LJ 参数 $\epsilon$ 和 $sigma$ 并将其设置为文献< | ||
+ | |||
+ | 仅通过指定粒子参数,LJ 势不会变得活跃。您还必须添加相应的势函数。为实现此目的,添加以下两行: | ||
+ | |||
+ | <code python> | ||
+ | potentialSet.addPotential(LennardJonesPotential(particleType1=' | ||
+ | particleType2=' | ||
+ | r_cut=9.0*Ang)) | ||
+ | potentialSet.addPotential(LennardJonesPotential(particleType1=' | ||
+ | particleType2=' | ||
+ | r_cut=9.0*Ang)) | ||
+ | </ | ||
+ | |||
+ | 这将在碳和锂之间以及锂和锂之间添加 LJ 势。注意,碳原子之间没有 LJ 势,因为这些相互作用由 Tersoff 势充分定义(LJ 只是一个微小的校正)。 | ||
+ | |||
+ | 在 **ATK-ForceField** 里的 Lennard-Jones 势的定义为 | ||
+ | |||
+ | $$V_{ij}(r) = 4\epsilon_{ij} \left [ \left(\frac{\sigma_{ij}}{r} \right)^{12} - \left(\frac{\sigma_{ij}}{r} \right)^{6} \right]$$ | ||
+ | |||
+ | 利用以下组合规则获得不同元素间的 $\sigma_{ij}$ 和 $\epsilon_{ij}$ (参见 [[https:// | ||
+ | |||
+ | 您应该始终检查为所选 LJ 参数推荐的组合规则是否与此定义一致。如果没有,您应该使用更为通用的 LennardJonesMNPotential。 | ||
+ | |||
+ | 在 [[https:// | ||
+ | |||
+ | |||
+ | |||
+ | {{ : | ||
+ | |||
===== MoS2 中的层内和层间内聚力 ===== | ===== MoS2 中的层内和层间内聚力 ===== | ||
+ | |||
+ | <WRAP center important 100%> | ||
+ | === 注意 === | ||
+ | 如果您正在使用 **ATK2014**,请确保已升级到 **ATK2014.3** (或更高版本)以运行本教程的这一部分,因为早期的版本可能会产生一些问题。 | ||
+ | </ | ||
+ | |||
+ | 在最后一个示例中,您将学习如何使用标记仅在所选原子上定义势。我们考虑一种辉钼矿(MoS< | ||
+ | |||
+ | 为设置系统,请打开 // | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | 将构型发送到 **ScriptGenerator**。 | ||
+ | |||
+ | 在 **ScriptGenerator** 中添加一个 NewCalculator 和 Optimize {{: | ||
+ | |||
+ | 将最终脚本发送给到 **Editor**。 | ||
+ | |||
+ | 在 python 脚本中,找到定义 Stillinger-Weber 势的行。复制该行,并按以下方式调整这两行: | ||
+ | |||
+ | <code python> | ||
+ | sw_layer1 = StillingerWeber_MoS_2013(tags=' | ||
+ | sw_layer2 = StillingerWeber_MoS_2013(tags=' | ||
+ | </ | ||
+ | |||
+ | 这将仅对属于相应标记组的原子间定义每个势。此外,您必须定义只在不同组的硫原子间起作用的 Lennard-Jones 势。这可以通过添加以下代码行来完成: | ||
+ | |||
+ | <code python> | ||
+ | # Define a new potential for the interlayer interaction. | ||
+ | lj_interlayer_potential = TremoloXPotentialSet(name=" | ||
+ | # Add particle type definitions for both types. | ||
+ | lj_interlayer_potential.addParticleType(ParticleType.fromElement(Molybdenum)) | ||
+ | lj_interlayer_potential.addParticleType(ParticleType.fromElement(Sulfur, | ||
+ | # Add Lennard-Jones potentials between the sulfur atoms of different layers. | ||
+ | lj_interlayer_potential.addPotential(LennardJonesPotential(' | ||
+ | lj_interlayer_potential.setTags([' | ||
+ | </ | ||
+ | |||
+ | Lennard-Jones 参数取自参考文献< | ||
+ | |||
+ | 最后,所有三个势的集合必须在 TremoloX 计算器中组合。为完成此操作,请添加以下行: | ||
+ | |||
+ | <code python> | ||
+ | # Combine all 3 potential sets in a single calculator. | ||
+ | calculator = TremoloXCalculator(parameters=[sw_layer1, | ||
+ | </ | ||
+ | |||
+ | 完成定义势之后,将计算器添加到块体构型上并执行几何体优化,脚本照常继续。您可以在文件 [[https:// | ||
+ | |||
+ | 如果您运行脚本,将会获得优化的晶格常数,其与参考文献 < | ||
+ | |||
===== 参考 ===== | ===== 参考 ===== | ||
+ | |||
+ | * [1] (1, 2, 3, 4, 5) M Matsui and M. Akaogi: Molecular dynamics simulations of the structural and physical properties of the four polymorphs of TiO2. [[https:// | ||
+ | * [2] (1, 2) C.J. Fennell and J.D. Gezelter: Is the Ewald summation still necessary? Pairwise alternatives to the accepted standard for long-range electrostatics. [[https:// | ||
+ | * [3] (1, 2, 3, 4, 5) J.P. Trinastic, R. Hamdan, Y. Wu, L. Zhang, H.-P. Cheng: Unified interatomic potential and enery barrier distributions for amorphous oxides. [[https:// | ||
+ | * [4] (1, 2, 3) O. N. Kalugin, V. V. Chaban, O. V. Prezhdo: Microscopic Structure and Dynamics of Molecular Liquids and Electrolyte Solutions Confined by Carbon NanoTubes: Molecular Dynamics Simulations. [[https:// | ||
+ | * [5] (1, 2) J.-W. Jiang, H.S. Park, T. Rabczuk: Molecular dynamics simulations of single-layer molybdenum disulphide (MoS2): Stillinger-Weber parameterization, | ||
+ | * 英文原文:[[https:// |