第一性原理计算在新材料研究中的应用

第一性原理 报告 软件系统测试报告下载sgs报告如何下载关于路面塌陷情况报告535n,sgs报告怎么下载竣工报告下载 USPEXMSVASPGaussianPd基Heusler合金Pd2CrGa晶体结构预测贵金属纳米粒子团簇结构预测FeZn13弹性常数的第一原理计算USPEX 一种专门用于结构预测的进化算法，结合外部电子结构计算程序（VASP，SIESTA和GULP等）来寻找从头算全局自由能最小值。采用遗传操作、原子置换、晶格突变等算子来得到尝试子结构，并对子结构进行筛选，只保留其中若干比例“优质”的结构参与下一代的遗传进化过程，从而快速地收敛到稳定相。USPEX 优点 成功地实现了对于任意给定温度、压强条件下，仅从材料化学成分组成进行晶体结构预测 无需实验数据，仅从材料的化学成分出发预测晶体结构，特别适用于高温、高压等极限条件下的结构预测。 支持各种晶胞结构的搜索。可以由实验得到的晶胞结构开始搜索，如晶胞参数、晶胞形状、晶胞体积等；也可以由已知和假设结构开始搜索。 通过几片结构的空间粘连，部分保留并考虑了原子的局域排布信息。反映了晶体中强的短程相互作用和当前一代的信息。对于处理较大的体系具有明显优势。 置换算法提供了用户自定义哪种原子相互交换的功能，特别适用于具有长程化学相似的不同种原子构成的体系。MaterialsStudio MaterialsStudio是专门为材料科学领域研究者开发的一款可运行在PC上的模拟软件. MS包含多种模块，其核心模块是visualizer，可运行于客户端PC，支持windows和Linux等系统。 MS包含有许多重要的功能模块，例如进行第一性原理计算的CASTEP模块，进行分子动力学计算的DMol3模块，进行几何机构优化与预测的GULP模块等等。MaterialsStudio文件窗口图形窗口属性窗口计算状态窗口菜单窗口Castep模块MaterialsStudioMaterialsStudioMaterialsStudio能带结构态密度VASP VASP是维也纳大学Hafner小组开发的进行电子结构计算和量子力学-分子动力学模拟软件包。 VASP通过近似求解Schrodinger方程得到体系的电子态和能量，既可以在密度泛函理论（DFT）框架内求解Kohn-Sham方程（已实现了混合泛函计算），也可以在Hartree-Fock（HF）的近似下求解Roothaan方程。 VASP采用周期性边界条件(或超原胞模型)处理原子、分子、团簇、纳米线（或管）、薄膜、晶体、准晶和无定性材料，以及表面体系和固体的问题。只能计算。建模和结果分析要依靠其他程序。纯计算软件，只有Linux版本终端窗口VASPVASPINCAR：控制计算的参数POSCAR：材料的描述(晶胞的尺寸，原子的位置）KPOINTS：布里渊区积分设置POTCAR：每个元素的赝势DOSCAR:包含了态密度的信息EIGENVAL:包含了能带的信息VASPVASP 优点 缺点 Castep 图形界面，操作方便建模和结果分析容易 在Windows下执行速度较慢。可选的计算 方法 快递客服问题件处理详细方法山木方法pdf计算方法pdf华与华方法下载八字理论方法下载 较少，结果不够精确 Vasp 代码稳定，执行速度快可以使用更高级的计算方法，使结果更精确 命令行操作，且Linux系统入门较难。建模和结果分析依赖其他软件。Gaussian Gaussian是一个功能强大的量子化学综合软件包。主要功能有对分子结构和能量、过渡态的能量和结构、化学键和反应能量、分子轨道、原子电荷和电势、红外和拉曼光谱等的计算。 研究周期性体系：Gaussian扩展了化学体系的研究范围，它可以用周期性边界条件的方法（PBC）模拟周期性体系 预测光谱：Gaussian可以计算各种光谱和光谱特性。包括：IR和Raman；预共振Raman；紫外-可见等等.第一性原理应用实例一 团簇是指由有限数目（几个至上千个）的原子、分子或离子通过一定的键合方式（包括物理和化学结合方式）构成的相对稳定的微观和亚微观聚集体，其空间尺度在10-10-10-8m之间。由于团簇特殊的空间尺度，其物理和化学性质随其包含原子数目的不同有很大的差异，产生许多奇异的特性。因此，团簇可视为是介于原子、分子与宏观固体之间物质结构的新层次或新凝聚态——即介观层次，代表了凝聚态物质的初始形态，常被称为“物质第五态”.贵金属纳米粒子团簇结构预测贵金属纳米粒子团簇结构预测碳原子碳60团簇金刚石晶体贵金属纳米粒子团簇结构预测由于团簇特殊的空间尺度，其物理和化学性质随其包含原子数目的不同有很大的差异，产生许多奇异的特性。 类型 典型代表 键合方式 平均结合能 范德瓦尔斯团簇 Arn、(N2)n 弱静电作用 0.3eV 分子团簇 In、(C6H6)n 弱静电作用 0.3-1eV 氢键团簇 (HF)n、(H2O)n 电荷转移特征 0.3-0.5eV 离子键团簇 (NaCl)n、(MgO)n 静电作用 2-4eV 共价键团簇 Cn、Sin、Gen 共价键结合 1-4eV 金属键团簇 Agn、Nan、Cun 电子结合 0.5-3eV贵金属纳米粒子团簇结构预测 贵金属团簇 由金属Au、Ag、Cu、Pt形成的金属键团簇称为贵金属团簇 近年来，贵金属团簇和纳米颗粒由于其独特的光学、电子学和催化性能而在生物学、医学、光学、催化和纳米电子学等领域引起了广泛的兴趣。贵金属纳米粒子团簇结构预测 贵金属团簇结构的理论研究 在各种团簇体系中，金团簇体系是目前研究得非常热门的领域之一。 金原子由于其相对原子质量较大，内层轨道电子的速度可以与光速相比拟，其4f、5d轨道电子又接近全充满，相对论效应显著，使其具备很多特殊性质。 大体系金团簇结构以高对称性的富勒烯结构或管状中空笼结构最为稳定。贵金属纳米粒子团簇结构预测 贵金属团簇的光学性质研究 由于贵金属团簇的特殊光学性能，其常被制成纳米复合结构材料进行研究 如Au/SiO2纳米复合薄膜，Ag/BaO复合薄膜 实际应用 表面增强拉曼信号SERS 生物传感器进行肿瘤定位及原位检测 纳米颗粒沉积提高太阳能电池效率贵金属纳米粒子团簇结构预测 存在问题 新的团簇结构主要靠人为排列组合进行探索，效率低下，且优化后得到的稳定结构常为局域能量最低点，需寻找高效的全局搜索方法得到最稳定基态结构。 国内外研究集中于单一金属团簇的电子性质，随着实验水平的提高，贵金属小团簇的吸收光谱得以研究，急需对其光学性质进行理论研究进行支持。 混合团簇研究较少，需研究掺杂原子的影响。贵金属纳米粒子团簇结构预测贵金属纳米粒子团簇结构预测具体操作 流程 快递问题件怎么处理流程河南自建厂房流程下载关于规范招聘需求审批流程制作流程表下载邮件下载流程设计 GaussView贵金属纳米粒子团簇结构预测计算结果Aun（2-20）团簇的结构贵金属纳米粒子团簇结构预测Aun（21-40）团簇的结构贵金属纳米粒子团簇结构预测Cun（2-20）团簇的结构贵金属纳米粒子团簇结构预测Cun（21-40）团簇的结构贵金属纳米粒子团簇结构预测AunCu1团簇的结构(2-20)贵金属纳米粒子团簇结构预测CunAu1团簇的结构(2-20)第一性原理应用二Pd基Heusler合金Pd2CrGa晶体结构预测* Heusler合金 Heusler合金因其独特的磁性形状记忆效应，而成为备受关注的一种新型功能材料。Heusler合金兼有大恢复应变、大输出应力、高响应频率等优良性能,是一种理想的驱动与传感材料,有着广泛的应用前景。 Heusler合金为高度有序排列的三元金属间化合物,它的晶格结构可以看作由四个互相贯通的fcc次晶格组成。A、B、C、D四个位置分别被不同的原子占据可以形成两种不同的结构,分别对应:Cu2MnAl型（空间群Fm-3m）和Hg2CuTi型(空间群F-43m)Heusler合金。Pd基Heusler合金Pd2CrGa晶体结构预测Heusler合金结构Pd基Heusler合金Pd2CrGa晶体结构预测 采用基于基因遗传算法软件包USPEX对Pd2CrGa进行晶体结构预测。 步骤：首先建立Pd2CrGa的L21结构模型，并采用自选极化方式对Pd2CrGa晶胞进行结构优化，体积优化，通过对一系列不同体积和能量的拟合，确定基态能量，进而获得结构的晶格常数。 通过基于密度泛函理论（DFT）的VASP软件包，采用经相对论校正的投影缀加波PAW(ProjectorAugmentedWave)方法，交换关联能采用广义梯度近似(GGA),平面波截断能选取500eV,K点网格采用12×12×12,计算过程均采用自旋极化的处理方式对Pd2CrGa晶体结构的磁性、态密度、弹性常数进行了计算。Pd基Heusler合金Pd2CrGa晶体结构预测 Pd2MGa（M=Cr,Fe）的结构（a）L21结构；（b）四方结构(-Pd原子-M原子-Ga原子)Pd基Heusler合金Pd2CrGa晶体结构预测* 磁性Pd2CrGa呈现铁磁性，Cr原子是Pd2CrGa总磁矩的主要来源。Pd基Heusler合金Pd2CrGa晶体结构预测从图中可以看出，Pd2CrGa在两种状态下均表现处很强的自旋极化，而且在费米面附近尤为明显Cr.原子的态密度差异是Pd2CrGa总态密度差异的主要来源，而Pd原子、Ga原子的自旋向上和自旋向下态密度的对称性较高，对总磁矩的贡献有限.所以，Cr原子是Pd2CrGa磁性的主要贡献者态密度Pd基Heusler合金Pd2CrGa晶体结构预测弹性常数 对于Pd2CrGa四方结构包含六个独立的弹性常数C11，C12，C13，C33，C44和C66. 四方结构的体积模量可由公式：B=1/9(2C11+2C12+4C13+C33)得到。结果四方结构的弹性常数满足四方相稳定判据：Pd2MnGa的四方结构满足稳定性判据.能稳定存在。 C11 C12 C13 C33 C44 C66 B Pd2CrGa 218.07 70.89 124.36 184.99 23.74 75.53 140.04第一性原理计算实例三FeZn13弹性常数的第一原理计算 自20世纪90年代起，热镀锌在汽车上、大型建筑上的应用迅速发展，新应用对镀锌板质量提出了更高的要求。热镀锌的几个相层中，含铁量为5%—6%的ζ相[1]（FeZn13相）是热镀锌过程中一个非常重要的相层。提高镀锌钢材的力学性能，改善锌镀层粘附性能，都需要深入细致的研究FeZn13相。 应用第一原理方法计算材料的电荷分布和弹性常数，可以获得材料的特征参数，从而能够表征、预测甚至设计材料的结构与性能，为以后对FeZn13的验研究和理论研究提供理论依据。FeZn13弹性常数的第一原理计算 材料弹性是人们选择和使用材料的依据之一；作为减振元件或结构则要求材料应变能要高。近代航空、航天、无线电及其精密仪器仪表工业对材料的弹性有更高的要求，不仅要有高的弹性模量，而且还要恒定，同时，准确测定材料的弹性常数，对于研究材料间原子的相互作用和相变都具有工程和理论意义。FeZn13弹性常数的第一原理计算FeZn13晶胞模型FeZn13晶胞的拓展图FeZn13晶格常数a=13.394Å,b=7.598Å,c=5.066Å,α=γ=90o,β=127.23º单斜立方结构，（空间群12：C2/m），晶胞中一共有28个原子，其中Z方向的四条棱的中间和单斜立方的中心各有一个Fe原子，晶胞中和面心位置一共有26个Zn原子。FeZn13弹性常数的第一原理计算Fe原子和Zn原子连接的4种键长晶胞中个原子的位置 x y z Fe1 0.0000 0.0000 ,0.5000 Zn2 0.0000 0.0000 0.0000 Zn3 0.2227 0.0000 0.0693 Zn4 0.1144 0.5000 0.2929 Zn5 0.1764 0.1792 0.5517 Zn6 0.0753 0.2898 0.8360 a=1.32313nm,b=0.76326nm,c=0.50567nm β=126.966o,Vol=0.2041nm3FeZn13弹性常数的第一原理计算 弹性常数是表征材料弹性的量。对大多数晶体来说，弹性常数有21个非零独立分量。而高对称性的晶体弹性常数还可以进一步减少。对于单斜晶体来说，一共有13个非零独立分量，要想得到这些弹性常数，只需对未变形的晶胞结构施加微小的变形，得到变形后晶胞的受力情况。 弹性常数的第一性原理计算可以分为两个阶段：(1)对晶体施加弹性极限范围内的应力。在计算中以矩阵的形式体现，并且针对每个独立的弹性常数分量矩阵的形式也会有所不同。 (2)计算应力作用下晶体发生变形前后的应力及能量的变化，根据胡克定律即可得到相应的弹性常数。FeZn13弹性常数的第一原理计算 利用基于密度泛函理论DFT[13][14]（Densityfunctionaltheory）的Abinit软件，采用经相对论校正的投影缀加波[15]（projectoraugmentedwave,PAW）方法。PAW方法精度高于一般的赝势计算而且计算时间少于全电子方法。交换关联相互作用采用GGA中的PBE形式，计算中对FeZn13的截断能取600eV，K点选取4×4×6。FeZn13弹性常数的第一原理计算 力-应变法是一种计算材料弹性常数的经典理论，根据胡克定理，弹性常数与应力张量应变张量之间的关系可表示为： 写成矩阵 FeZn13弹性常数的第一原理计算对于单斜晶体应力-应变关系为：对单斜晶体施加如下变形矩阵：FeZn13弹性常数的第一原理计算 将p1对应的应变带入到单斜晶体的应力-应变关系中，得到，即可求出C22带入其它的变形矩阵则可以得到其它的对应的弹性常数。 计算得到的FeZn13的13个独立的弹性常数如表所示： C11 C12 C13 C16 161.81 54.07 53.13 7.17 C22 C23 C26 C33 158.79 53.65 5.67 112.24 C36 C44 C45 C55 6.70 41.24 1.55 50.70 C66 34.62THANKYOU!