|
|
本帖最后由 nusiew 于 2021-8-26 08:59 编辑
请问大家有没有听说过 hull 和 inhedron 两个程序？
各位老师、同学、前辈，大家晚上好！最近看文献，发现一种利用轨迹输出文件（由从头算分子模拟得到），结合计算几何，计算元素在两相间分配比的方法。对于其利用计算几何计算分配比的部分非常感兴趣，在此向大家求教。论文作者关注的领域是地核与岩浆中元素的分配，采用类似的方法一共发表了3篇文章，最近的一篇是2020年，在此举其2012年发表在PNAS上的一篇对其研究方法进行说明。
作者利用vasp进行第一性原理分子动力学模拟，想要计算碳元素在铁水和硅酸盐熔体两相间的分配比。模拟和计算的过程是：先将260个原子（包含125个O、51个Mg、44个Fe 、36个Si和4个C原子）随机均匀分布在模拟盒子中，跑若干步后盒子中的原子会合并成铁水和硅酸熔体两相（如下图所示），之后便是利用计算几何的方法，结合 hull 和 inhedron 两个程序对Fe相成分和硅酸盐熔体组成进行了计算。不知道大家有没有听说过 hull 和 inhedron 这两个程序？由于本人对计算几何不太了解，对论文中计算几何部分的操作步骤也不熟悉，不知道有没有相关方面的教程或书籍？希望大家多多指导。论文中具体的成分分析过程大致如下。
在平衡之后的每个 FPMD 步骤中，首先为 Fe 簇构建多面体，然后利用计算几何学中的概念和方法判断原子是否被包围在多面体中，从而获得包围在 Fe 相中的原子数。基本上，需要三个步骤。首先，计算 Fe 簇的质心并将其移动到模拟单元的中心。通过使用模拟单元的周期性边界条件，所有原子都相应地移动。这是使用计算几何的准备步骤。其次，通过使用程序“hull”，它采用随机增量算法。第三，通过使用随机射线交叉算法的程序“inhedron”确定原子是否在第二步中构造的多面体内部。在每个 FPMD 步骤中，Fe 相中的原子数在模拟的整个积累周期内相加并求平均值，以计算 Fe 相的组成。由于整个系统的整体组成是已知的，因此可以直接推导出硅酸盐熔体的组成。
模拟得到的两相，Fe相，初始成分和原子数，以及平衡两相成分和原子数如下：
具体的计算方法大家可以参考这三篇文献
2012：Carbon and other light element contents in the Earth’s core based on first-principles molecular dynamics
https://www.pnas.org/content/109/48/19579
2009：Partitioning of Si and O between liquid iron and silicate melt: A two-phase ab-initio molecular dynamics study（分析过程最详细）
https://agupubs.onlinelibrary.wi ... 0.1029/2009GL039751
2020：Magnesium partitioning between silicate melt and liquid iron using first-principles molecular dynamics: Implications for the early thermal history of the Earth's core
https://www.sciencedirect.com/sc ... 821X19306260#br0620
|
|
发表于 Post on 2021-8-24 19:58:23
|
只看该作者 Only view this author
|只看大图 View big picture
|倒序浏览 Reverse view
|阅读模式 Reading model
收藏 Add to favorites