ATAT-mcsqs生成准随机结构(SQS)更新

    通常使用第一性原理计算某些多元素占据原胞中同一位置的结构会优先考虑使用准随机结构(special quasirandom structure,SQS)来进行模拟建模。此篇教程意在整理一个较为简便的操作流程,以供参考。

合金理论自动化工具包(ATAT)1是一个通用名称,指的是合金理论的集合 Axel van de Walle2与多个研究小组合作开发的工具。而由于高熵合金内部多元、混合等特性,ATAT工具包内的mcsqs模块所涉及的特殊准随机结构Special quasi-random structure(SQS)于此有了非常大的应用需求。

部分有关mcsqs模块的使用教程如下

1、建无序格点参杂&合金模形工具包ATAT的安装与其各模块的详细使用例子

2、ATAT-mcsqs快速建立大体系高熵合金快速操作过程

3、高熵合金的mcsqs建模(https://www.bilibili.com/video/BV1QK4y1h752?spm_id_from=333.337.search-card.all.click&vd_source=4a437e68a7f05cd2b6b8232814f71b42)

4、ATAT(合金理论自动化程序)教程(上)

5、ATAT(合金理论自动化程序)教程(下)

首先,准备原始单一元素(共格的多元素合并为一种)晶胞,并通过vaspkit获得其最原始晶胞。

例如这里准备了一个面心立方的结构,组分单一,面心和顶角元素相同,这里设置为Si。

图片

Si                                      1.00000000000000          4    0.0000000000000000    0.0000000000000000     0.0000000000000000   4    0.0000000000000000     0.0000000000000000    0.0000000000000000    4   Si4Direct  0.0000000000000000  0.0000000000000000  0.0000000000000000  0.0000000000000000  0.5000000000000000  0.5000000000000000  0.5000000000000000  0.0000000000000000  0.5000000000000000  0.5000000000000000  0.5000000000000000  0.0000000000000000

通过vaspkit 602功能获得其原胞

vaspkit -task 602  +---------------------------------------------------------------+ |         VASPKIT Standard Edition 1.3.5 (03 Jul. 2022)         | |            Running VASPKIT Under Command-Line Mode            | +---------------------------------------------------------------+ -->> (01) Reading Structural Parameters from POSCAR File... +-------------------------- Summary ----------------------------+                           Prototype: A           Total Atoms in Input Cell:   4     Lattice Constants in Input Cell:   4.000   4.000   4.000        Lattice Angles in Input Cell:  90.000  90.000  90.000       Total Atoms in Primitive Cell:   1 Lattice Constants in Primitive Cell:   2.828   2.828   2.828    Lattice Angles in Primitive Cell:  60.000  60.000  60.000                      Crystal System: Cubic                       Crystal Class: m-3m                     Bravais Lattice: cF                         Space Group: 225                         Point Group: 32 [ Oh ]                       International: Fm-3m                 Symmetry Operations: 192 +---------------------------------------------------------------+ -->> (02) Written PRIMCELL.vasp file. +---------------------------------------------------------------+

原胞信息如下,只存在一个wyckoff点位。

Primitive Cell  1.000000    0.00000000000000    2.00000000000000    2.00000000000000    2.00000000000000    0.00000000000000    2.00000000000000    2.00000000000000    2.00000000000000    0.00000000000000  Si   1DIRECT    0.0000000000000000    0.0000000000000000    0.0000000000000000    Si1

图片

将vaspkit生成的原胞文件PRIMCELL.vasp复制为POSCAR,然后通过vaspkit 的414功能(VASPKIT版本需低于1.5)将POSCAR转化为ATAT的输入文件lat.in

图片

vaspkit -task 414  +---------------------------------------------------------------+ |         VASPKIT Standard Edition 1.3.5 (03 Jul. 2022)         | |            Running VASPKIT Under Command-Line Mode            | +---------------------------------------------------------------+ -->> (01) Reading Structural Parameters from POSCAR File... -->> (02) Written lat.in File! +---------------------------------------------------------------+

lat.in文件内容如下,后续计算需要重命名为rndstr.in

图片

根据实际需要的某一wyckoff点位的元素比例修改rndstr.in文件。这里修改为Si和C各占50%。修改后的rndstr.in文件如下

     2.82843     2.82843     2.82843   60.00   60.00   60.00    1.00000000    0.00000000    0.00000000    0.00000000    1.00000000    0.00000000    0.00000000    0.00000000    1.00000000     0.00000000    0.00000000    0.00000000  Si=0.5,C=0.5

注意:各点位各元素之间的比例与扩包比例相关联,最后确保扩包后各元素原子数量总数为整数。

因VASPKIT版本问题,新版本不再支持转换lat.in文件,建议使用ATOMKIT进行格式转换。

注:笔者关注到目前众多使用者的需求是获得晶格矢量正交的惯用胞超胞,在此特意更新注明,有两个途径来达成此需求。

一、跳过刚才的一步获得原胞的过程,直接基于惯用晶胞来生成rndstr.in文件。

  4 4 4 90 90 90      1.00000000    0.00000000    0.00000000    0.00000000    1.00000000    0.00000000    0.00000000    0.00000000    1.00000000   0.0000000000000000  0.0000000000000000  0.0000000000000000  Si=0.5,C=0.5  0.0000000000000000  0.5000000000000000  0.5000000000000000 Si=0.5,C=0.5  0.5000000000000000  0.0000000000000000  0.5000000000000000 Si=0.5,C=0.5  0.5000000000000000  0.5000000000000000  0.0000000000000000 Si=0.5,C=0.5

在后续计算过程中,会提示非最基本原胞的警告,但不影响计算和生成结构。

二、更改sqscell.out文件(请先阅读后续流程,熟悉晶格矢量变化后自行换算)

这里选择扩包方式为2×2×2,最后超胞总原子数量为8。

先使用ATAT获得in文件的对称性信息,-2=x需要根据结构实际情况修改,过小会导致ATAT出现段错误。

corrdump -nop -noe -2=4 -ro -l=rndstr.in  -clus ; getclus

再根据超胞的总原子数获得结构信息,-n=后面的数值与超胞数量相同

mcsqs -n=8

ATAT产生的超胞结构信息文件为sqscell.out文件

图片

根据既定扩包方式修改为:

1
2 0 00 2 00 0 2

生成超胞准随机结构:

mcsqs -rc

产生的准随机结构文件为:bestsqs.out

图片

导出为cif文件

str2cif < bestsqs.out > pos.cif

图片

图片

推荐:使用VESTA进行cif向POSCAR的转换.

通过vaspkit 将cif文件转化为POSCAR

vaspkit -task 105  +---------------------------------------------------------------+ |         VASPKIT Standard Edition 1.3.5 (03 Jul. 2022)         | |            Running VASPKIT Under Command-Line Mode            | +---------------------------------------------------------------+ Input the filename in cif format (e.g., POSCAR.cif): ------------>>pos.cif
 -->> (01) Written POSCAR file!

更为简便的方法可使用qvasp软件一键转化

qvasp -c2p

当然也可以通过VESTA导出。

POSCAR文件如下

图片

<think>好的,我现在需要回答用户关于Alloy Theoretic Automated Toolkit (ATAT)的使用和文档的问题。首先,我需要确认用户的需求是什么。用户可能是一位刚开始接触ATAT的研究人员,想要了解如何安装、使用这个工具包,以及相关的文档资源在哪里可以找到。 首先,我应该回忆一下ATAT的基本信息。根据用户提供的引用[^1],ATAT是由Axel van de Walle开发的,用于合金理论研究,包括相稳定性、热力学性质的计算和模拟。接下来,我需要整理ATAT的主要用途,比如处理合金的相图、原子间势能拟合、蒙特卡洛模拟等,这样用户能明白它能解决哪些问题。 然后是安装步骤。用户提到环境备,可能需要说明系统要求,比如Unix/Linux系统,需要安装哪些依赖,比如C++编译器、MPI(如果有并行计算需求)、Perl、Python等。安装步骤可能包括下载源码、解压、编译,或者使用预编译的版本。需要确保步骤清晰,避免用户遇到常见问题。 接下来是使用入门。用户可能需要知道如何运行ATAT,比如基本的命令行操作,可能需要举例说明几个常用命令,比如生成结构、运行能量计算、使用mcsqs进行团簇展开等。这部分可能需要结合文档中的示例,或者参考官方教程。 文档资源方面,用户需要知道在哪里可以找到详细的指南。ATAT的官方网站可能有手册、教程和示例,用户需要访问这些资源。同时,用户社区和邮件列表也是重要的支持渠道,可以推荐用户加入以获取帮助。 然后,我需要生成相关问题。这些问题应该围绕用户可能进一步关心的内容,比如具体模块的使用方法、如何备输入文件、如何分析输出结果、常见错误处理等。此外,可能还有关于ATAT在不同操作系统上的兼容性,或者如何扩展其功能的问题。 需要注意的是,用户可能对具体的命令行参数或脚本编写有疑问,因此在回答中应建议查阅官方文档,并可能提供一些典型用例的步骤。同时,确保所有技术术语和步骤用中文清晰表达,避免歧义。 最后,检查是否符合之前的系统级指令:使用正确的中文,数学公式用$...$,独立公式用$$...$$,但在这个问题中可能数学公式较少,但如果有提到相关公式,比如原子间势能的表达式,需要正确应用格式。同时,引用标识要正确添加,确保用户能追溯到引用来源。 总结回答的结构应该是:概述ATAT的用途,安装步骤,使用入门,文档资源,相关问题。确保各部分信息确,参考引用[^1]中的内容,并引导用户查阅官方文档获取最新信息。</think>### Alloy Theoretic Automated Toolkit (ATAT) 使用与文档指南 #### 1. **ATAT 概述** ATAT 是一个专注于合金理论研究的工具集合,主要用于计算合金的相稳定性、热力学性质及原子尺度模拟。其核心功能包括: - **团簇展开**:用于构建合金能量模型。 - **蒙特卡洛模拟**:预测合金的平衡态或非平衡态行为。 - **相图计算**:分析多组分系统的稳定相。 - **结构生成**:枚举合金的可能原子构型。 --- #### 2. **安装与环境备** **系统要求**: - 类 Unix 系统(如 Linux/macOS)。 - 需安装 C++ 编译器(如 `gcc`)、Perl 和 Python 解释器。 - 部分工具需要 MPI 支持并行计算(如 `map` 模块)。 **安装步骤**: 1. 从 [ATAT 官网](https://www.brown.edu/Departments/Engineering/Labs/avdw/atat/) 下载源码包。 2. 解压并编译: ```bash tar -zxvf atat_X.X.X.tar.gz cd atat_X.X.X make ``` 3. 将 `bin` 目录添加到系统路径: ```bash export PATH=$PATH:/path/to/atat/bin ``` --- #### 3. **快速入门** **示例:生成随机合金结构** 1. 使用 `gensquare` 生成初始结构: ```bash gensquare Al Ni > str.out ``` 2. 运行能量计算(需与第一性原理软件如 VASP 结合): ```bash runstruct vasp ``` 3. 使用 `mcsqs` 生成特殊随机结构SQS): ```bash mcsqs -n=16 ``` **关键命令**: - `map`:并行任务调度器。 - `phb`:计算相图。 - `clusterexpansion`:拟合团簇展开参数。 --- #### 4. **文档与资源** - **官方手册**:包含所有模块的详细说明,位于安装目录的 `doc` 文件夹。 - **教程与示例**:官网提供常见工作流程的案例(如团簇展开、相图计算)。 - **社区支持**:通过 [ATAT 邮件列表](https://www.brown.edu/Departments/Engineering/Labs/avwd/atat/mailing_list.html) 提问或查阅历史讨论。 ---
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值