利用MS的CASTEP模块计算铁磁性铁声子谱
声子是固体物理中的一个重要概念，它提供了使用大量重要性能的方法，例如比热、热膨胀、热传导、电子－声子相互作用、电阻系数和超导电性等。密度泛函理论(DFT)方法能预测这些性能，CASTEP提供了必要的功能。在点阵动力学计算中有两种主要的方法：密度泛函微扰理论(DFPT)和有限位移方法。第一种方法通常比较快速并且更加精确，但是它的执行是有问题的，并且受到一系列的限制。目前，在CASTEP中DFPT只能用于具有fixed occupancies(绝缘体)、没有自旋极化并且只能对norm-conserving赝势使用。这样，对广泛的材料种类 (包括磁性材料和金属)，声子计算只能通过采用有限位移算法来执行。
注意：本教程中的计算对CPU时间和内存需求方面的要求很高。
介绍：
在本教程中，将学习如何使用CASTEP来执行有限位移运算以获得磁性金属的声子散射和态密度。
1. 优化Fe晶胞的结构
从输入Fe的结构开始，它包含在Materials Studio提供的结构库中。
从菜单栏选择File | Import...，定位到Structures/metals/pure metals并选择Fe.msi。
通过将晶胞转化为原胞通常能减少计算时间。
从菜单栏选择Select Build | Symmetry | Primitive Cell。
显示Fe原胞。
现在使用CASTEP 优化Fe的几何结构。
从工具栏选择CASTEP工具，然后选择Calculation，或者从菜单栏选择Modules | CASTEP | Calculation。
显示CASTEP Calculation对话框。
在Setup选项卡上，将Task从Energy改为Geometry Optimization，设置Quality为Medium，设置Functional为LDA。选中Spin polarized复选框，取消选择Use formal spin as initial。设置Initial spin值为2。
LDA/CA-PZ局域交换关联泛函被认为是可获得的最准确的描述之一，将initial spin值设为2是因为我们正在模拟铁磁性的Fe晶体。
几何优化的默认值不包含晶胞的优化。
点击More...按钮，在CASTEP Geometry Optimization对话框上，选中Optimize cell，关闭对话框。
选择Electronic选项卡并设置Pseudopotentials值为Ultrasoft。
选择Job Control选项卡，为job选择Gateway位置并设置Runtime optimization为Memory。
点击More...按钮，打开CASTEP Job Control Options对话框，在Live updates区域取消选择所有选项，关闭对话框。
点击Run按钮启动job。
2. 计算声子散射和声子态密度(DOS)
为了计算声子散射和声子态密度，必须执行单点能量计算，并为计算选择适当的性能。
确保Fe CASTEP GeomOpt目录中的Fe.xsd是活动文档。
在CASTEP Calculation对话框上选择Setup选项卡，并设置Task为Energy。
在Properties选项卡上，选择Phonons复选框，选择Both 选项以选中态密度和散射。取消选择Calculate LO-TO splitting，将Method选为Finite displacement。
纵向光学横向光学(LO)劈裂不能对金属进行计算，因为它们在Γ-point是相同的。有限位移方案被设计用于金属和自旋极化系统(以及为那些具有高效超软势能的系统)。对铁磁性Fe这是理想的计算声子势能的方法。
点击More...按钮，显示CASTEP Phonon Properties Setup对话框。确保Method为Finite displacement。设置Supercell defined by cutoff radius值为3.6 Å。将Dispersion和Density of states的Quality值都设置为Fine。关闭CASTEP Phonon Properties Setup对话框。
注意：Cutoff radius的选择对有限位移运算是至关重要的参数。当使用较大的cutoff radius值时精度较高，因为这时考虑了更多的近邻。然而，随着该值的增加，计算时间增加的非常迅速。出于实际原因，在本教程中，对该参数选择了较小的值。声子频率的收敛作为cutoff radius的函数在执行有意义的实验计算时应该被研究。
选择Job Control选项卡并为计算选择Gateway。
点击More...按钮，打开CASTEP Job Control Options 对话框并选中所有Live updates选项，关闭对话框。
点击Run按钮。
Job被提交并开始运行。一个叫做Fe CASTEP Energy的新目录被建立在Fe CASTEP GeomOpt目录中。当能量计算结束时，新的结果文件被放置在该目录中，包括Fe_PhonDisp.castep和Fe_PhonDOS.castep。
3. 显示声子散射和态密度
从Materials Studio菜单栏选择Modules | CASTEP | Analysis，从性能列表中选择Phonon dispersion，确保Results file选择器显示Fe_PhononDisp.castep。
从Units下拉列表选择cm-1，从Graph style下拉列表选择Line。
点击View按钮。
一个新的图表文件－Fe Phonon Dispersion.xcd被建立在结果目录中。它看起来如下图所示。
声子散射的实验值如下图所示。
总体来说，计算的精度是可以接受的。可以通过用更大的cutoff radius值运行计算来得到与实验结果之间更好的一致性。
现在生成声子态密度曲线。
从Materials Studio菜单栏选择Modules | CASTEP | Analysis ，从性能列表中 选择Phonon density of states，确保Results file选择器显示Fe_PhonDOS.castep。
设置DOS display为Full，点击More... 按钮，打开CASTEP Phonon DOS Analysis Options对话框，选择Interpolation作为Integration method，设置Accuracy level为Fine。点击OK按钮。
在CASTEP Analysis对话框上点击View按钮。
建立了一个新的图表文件－Fe_PhonDOS Phonon DOS.xcd。如下图所示。
最后，欢迎大家通过微信公众号联系我们。
微信公众号：320科技工作室。
利用MS的CASTEP模块计算铁磁性铁声子谱相关推荐
- 基于MS的Dmol3模块计算离子液体的电子转移和成键
离子液体是由氢键供体(HBDs)和氢键受体(HBAs)按照一定比例混和组成的低共熔混合物,常温下呈液态.氯化胆碱由于其价格低廉.来源广泛.可生物降解等特点成为了最常用的氢键受体之一,形成了一系列胆碱类 ...
- MS CASTEP模块的实战策略
1.模型选定 我们在进行材料物理模拟所需要做的第一步 (也是很重要的一步) 是模型的选定或建构.CASTEP 虽然内建了很多功能来预测晶胞参数 (边长,夹角) 与原子位置,但仍然仰赖使用者告诉它 &q ...
- Materials Studio中的Castep模块金属矿物计算与分析
金属矿物属于大型周期性体系,因此要运用密度泛函理论研究其矿物性质要运用Mterials Studio中的Castep模块.本文主要讲述氧化锌矿物晶胞模型的构建以及完全解理面的计算和对氧化锌矿物表面性质 ...
- pythontime模块计算时长_用python的time模块查看你出生多长时间了
time模块的应用查看出生了多长时间 查看银行卡解封时间 1.查看出生了多长时间 这里我们使用time模块计算我们出生了多长时间 import time # 导入时间模块 b = input('你的出 ...
- python分数类_利用标准库fractions模块让Python支持分数类型的方法详解
前言 你可能不需要经常处理分数,但当你需要时,Python的Fraction类会给你很大的帮助.本文将给大家详细介绍关于利用标准库fractions模块让Python支持分数类型的相关内容,分享出来供 ...
- python利用tushare下载数据并计算当日收益率
python利用tushare下载数据并计算当日收益率 计算股票收益率的程序主要有以下几部分构成: 1.获取股票接口数据函数:pro_daily_stock() 2.计算收益率函数:cal_stock ...
- 体重计c语言程序,利用MSP430F149和HX711模块制作体重秤？
已结贴√ 问题点数:20 回复次数:1 利用MSP430F149和HX711模块制作体重秤? 以下代码是我从另一个论坛大拿那里看来的,如何改动才能只保留数码管显示体重呢?不需要蓝牙模块等等.救救毕设做 ...
- python基础教程：Python中利用sqrt()方法进行平方根计算的教程
这篇文章主要介绍了Python中利用sqrt()方法进行平方根计算的教程,是Python学习的基础知识,需要的朋友可以参考下 sqrt()方法返回x的平方根(x>0). 语法 以下是sqrt() ...
- 基于MS的Forcite模块模拟氢在钨表面的物理吸附
目的:介绍Forcite在生成分子在表面的优化结构中的使用. 模块:Forcite 背景: 在Materials Studio中Forcite是一个经典分子力学工具,由Accelrys科学家和软件工程 ...
最新文章
- 如何使用Nexus搭建Maven私服
- 云炬随笔20211002
- leetcode-7-整数翻转
- VMware如何进入安全模式
- [python作业AI毕业设计博客]大数据Hadoop工具python教程１-HDFS Hadoop分布式文件系统...
- 用Asp.net写自己的服务框架
- iOS TableView多级列表
- Python实现桌面程序：PyQt5 + QtDesigner -- 界面设计与逻辑编写
- SOTIF很快将会取代ISO 26262？为您详细解读SOTIF标准ISO/PAS 21448
- 怎样把ue4官方文档下载下来_ue4官方文档下载及翻译
- 多媒体计算机特性,多媒体计算机的基本特性
- 离散数学 2.3 Functions
- [MATLAB]matlab鼠标操作画两圆,做出两圆的公切线
- 单元测试中不能加入参数
- Linux系统配置网卡ip地址
- 中国日报聚焦游戏陪玩 直播平台加速拓展百亿市场
- 通过不同的域名来访问根目录和根目录下的子目录
- 彩灯控制器课程设计vhdl_基于VHDL的彩灯控制器设计与实现毕业设计(DOC毕业设计论文)...
- 一个初学Hibernate很好的网址
- 中国最伟大的杀人犯（海龙哥）