实验I 铜材料体相性质的理论表征.docVIP

实验1 铜材料体相性质的理论表征 实验要求：了解固体材料晶相点阵结构和空间群；掌握材料系统建模软件MS(Material Studio)构建材料晶格素单胞和复单胞的方法，并会给晶胞参数赋值和获取原子坐标数据信息；会使用Linux操作系统和文本编辑器Vi的常用操作命令，完成相应指令和任务；掌握用PW91相关泛函计算自由原子和凝聚相的能量，并计算金属铜的内聚能；会利用态方程(EOS)拟合二次曲线计算面心立方的晶格参数；掌握三阶Birch-Murnaghan态方程计算材料的体积模量 实验目的：利用基于密度泛函理论的第一性原理计算铜体相的几何结构、内聚能和体积模量，并与仪器实验结果和以前理论研究结果比较，获得计算面心立方晶格材料体相性质的常用理论方法 实验原理：很多实际的应用要求精确地测定晶格参数，如测定相变、膨胀系数、弹性、热电效应、组分变化、应力等。理论计算晶格常数多采用固体的态方程(Equation of State,EOS),它是指晶体体系的总能量E或感知的压强p随体系体积V的变化，即具有E(V)或p(V)形式。E(V)曲线的行为如图1所示。 晶体的状态方程在基础和应用科学中具有重要的意义。比如平衡体积(V0)，体弹性模量(B0)及它的一阶导数(B0)。这些可测的物理性质直接与晶体的状态方程有关。在高压下的状态方程有好几种不同函数形式来描述，比如Murnaghan方程、Brich-Murnaghan(BM)方程和普适方程等Murnaghan状态方程由Murnaghan在1944年提出，总能量(E)和体积(V)的函数关系在Murnaghan态方程中描述为 其中E0，V0，B0和B0′分别为体系在零压下的能量，体积和弹性模量及其弹性模量对压强的一阶倒数。体弹性模量定义为 其中T为所研究体系温度，即材料在外力作用下产生单位弹性变形所需要的应力弹性模量可视为衡量材料产生弹性变形难易程度的指标，其值越大，使材料发生一定弹性变形的应力也越大，即材料刚度越大，亦即在一定应力作用下，发生弹性变形越小。Birch-Murnaghan(BM)态方程是由Brich把Gibb自由能按Eulerian应变进行展开得到的 其中， 由方程(1)可看出，到V=V0时，E=E0。也就是说，我们可以通过Murnaghan状态方程,在平衡晶格常数a附近，选取若干E-a数据点经曲线拟合得体系最低能量时对应的格子常数即为平衡晶格常数。根据计算得到的E-V数据点按上面的等式采用最小平方差拟合得到B0，B0和V。 内聚能是衡量聚集态物质间作用力的参数 其中，材料组成原子M(g)和M(s)分别处于气态(gas)和固态(solid)。故有 其中，H=E+pV表示相应体系的焓，其值可通过计算E和V而得。对于固体中的原子，体系能量仅由电子能量Eelec(s)和原子振动能Ev所贡献。Eelec(s)直接由第一性原理计算得到，而Ev可由统计热力学给出： 其中(i表示晶格原子的振动模式i的频率，h为普朗克常数，kB为玻尔兹曼常数。对气体原子来说，体系能量没有振动贡献，所以由电子能量Eelec(g)、转动能量Erot和能量平动Etrans组成，即 固体的体积V(s)可由单包体积算出，而气体的体积V(g)对焓贡献可由气体状态方程pV=kT而得。 实验步骤： 实验任务一 利用态方程EOS拟合计算铜的晶格常数a 1a)在Cu晶格格子常数实验值a(exptl)=3.615 ?附近取5个点，即a值分别为： 3.515 ?, 3.565 ?, 3.615 ?, 3.675 ?和3.725 ?。 1b)利用MS软件构造具有相应晶胞参数的铜素单包，如图2所示。并生成扩展名为cell的坐标文件。 1c)在Linux系统里，用vi创建坐标输入文件POSCAR,其内容如下。 bulk cu 1. -0.000000000000001 1.807500000000000 1.807500000000000 1.807500000000000 -0.000000000000001 1.807500000000000 1.807500000000000 1.807500000000000 0.000000000000001 cu 1 selective dynamics direct 0.0000000000000000 0.0000000000000000 0.0000000000000000 T T T 1d)创建倒空间输入文件KPOINTS,其内容如下。 k-Points sampling 0 Monkhorst-Pack 15 15 15 0 0 0 1e)创建密度泛函理论计算参数文件