转载&如何用VASP计算晶格常数
我们用Pd金属作为例子。
Pd金属的实验上的晶格常数为3.89A。在这里，我们用VASP计算它的晶格常数。
首先将Pd所对应的POTCAR文件拷贝到目录下。然后准备好INCAR和KPOINTS文件。POSCAR文件我们将通过一个tcsh的script来产生。
&& KPOINTS文件可以如下：
&& Monkhorst Pack
&& Monkhorst Pack
&& 11 11 11
&& INCAR文件可以如下：
&& SYSTEM = Pd bulk
calculation
&& Startparameter for this
&& PREC = Accurate
&& ISTART = 0 job : 0-new 1-cont
&& ICHARG = 2 charge: 1-file
2-atom 10-const
&& ISPIN = 1 spin polarized
calculation?
&& Electronic Relaxation 1
&& EDIFF = 0.1E-03
stopping-criterion for ELM
&& LREAL = .FALSE. real-space
projection
&& Ionic relaxation
&& EDIFFG = 0.1E-02
stopping-criterion for IOM
&& NSW = 0 number of steps for
&& IBRION = 2 ionic relax: 0-MD
1-quasi-New 2-CG
&& ISIF = 2 stress and
relaxation
&& POTIM = 0.10 time-step for
ionic-motion
&& TEIN = 0.0 initial
temperature
&& TEBEG = 0.0; TEEND = 0.0
temperature during run
&& DOS related values:
&& ISMEAR = 0 ; SIGMA = 0.05
gaussian smear
&& Electronic relaxation 2
&& Write flags
&& LWAVE = F write WAVECAR
&& LCHARG = F write CHGCAR
&& 产生POSCAR和计算晶格常数的工作可以用以下的PBS
script来完成。
&& #!/bin/tcsh
&& #PBS -S /bin/sh
&& #PBS -l
nodes=4:athlon:ppn=2
&& #PBS -l cput=384:00:00
&& #PBS -m ae
&& #PBS -o output
&& #PBS -e error.log
&& # set parameter
&& set EXEC = 'vasp'
&& set SRC =
'/usr/common/executable'
&& # change working
&& cd $PBS_O_WORKDIR
&& # copy fresh executable from
depository
&& cp -f $SRC/$EXEC .
&& # execute mpi program
&& foreach a (3.3 3.4 3.5 3.6
&& echo "a= $a"
&& cat &POSCAR
&& 0.5 0.5 0.0
&& 0.0 0.5 0.5
&& 0.5 0.0 0.5
&& 0.0 0.0 0.0
&& 0.25 0.25 0.25
&& mpiexec -nostdin ./$EXEC
&& set E=`tail -2 OSZICAR`
&& echo $a $E
&& # remove executable
&& rm -f $EXEC
&& 如果不用不需要用PBS
script，则更加简单，如下即可。将其命名为lattice。
&& #!/bin/tcsh
&& foreach a (3.5 3.6 3.7 3.8 3.9
4.0 4.1 4.2)
&& echo "a= $a"
&& cat &POSCAR
&& 0.5 0.5 0.0
&& 0.0 0.5 0.5
&& 0.5 0.0 0.5
&& cartesian
&& 0.0 0.0 0.0
&& set E=`tail -1 OSZICAR`
&& echo $a $E
&& 用chmod +x
lattice，将其改为可执行文件。然后在命令行里键入./lattice 即可。
以下是用USPP-LDA运行完后的SUMMARY文件。每个计算用时13秒。 （在USPP中Pd的截断能量是198.955）
&& 3.5 1 F= -. E0=
-. d E =-.
&& 3.6 1 F= -. E0=
-. d E =-.
&& 3.7 1 F= -. E0=
-. d E =-.
&& 3.8 1 F= -. E0=
-. d E =-.
&& 3.9 1 F= -. E0=
-. d E =-.
&& 4.0 1 F= -. E0=
-. d E =-.
&& 4.1 1 F= -. E0=
-. d E =-.
&& 4.2 1 F= -. E0=
-. d E =-.
用抛物线拟和得到的晶格常数为$3.888\AA$,固体中每个原子的能量是$E_{bulk}=-6.4257$。
以下是采用PAW-LDA势运行完以后的SUMMARY文件。每个计算用时20秒。所以相对来说PAW势所需要的时间多一些，这是因为PAW势的
energy cutoff相对比较高（在PAW中Pd的截断能量是250.832）。
&& 3.5 1 F= -. E0=
-. d E =-.
&& 3.6 1 F= -. E0=
-. d E =-.
&& 3.7 1 F= -. E0=
-. d E =-.
&& 3.8 1 F= -. E0=
-. d E =-.
&& 3.9 1 F= -. E0=
-. d E =-.
&& 4.0 1 F= -. E0=
-. d E =-.
&& 4.1 1 F= -. E0=
-. d E =-.
&& 4.2 1 F= -. E0=
-. d E =-.
用抛物线拟和得到的晶格常数为$3.875\AA$,固体中每个原子的能量E_bulk=-6.4185eV
可见，PAW-LDA和USPP-LDA给出的晶格常数都和实验吻合的非常好，两者之间的差别也很小。在以下所有的计算中，如果没有特殊声明，我们都默认采用PAW-LDA的势。
&& 结合能(cohesive
energy)的定义如下：
&& -E_coh =
[E_bulk-N*E_atom]/N
所以我们要将固体中每个原子的能量减去单个Pd原子的能量，才能和实验的结合能相比较。对于过渡金属原子，计算单个原子的能量要特别注意。VASP的网页
上给出了求结合能所需的单个原子能量的修正值（详见VASP手册Pseudopotentials supplied with the
package一章）。可以在上面查到，Pd每个原子LDA的修正值为1.46eV。所以我们得到LDA近似下Pd的结合能为4.998eV。此值和实验
值比严重偏大，这是因为LDA
通常成键过强的关系。如果我们改用GGA的赝势，可以得到和实验比较吻合的结果。&
已投稿到：
以上网友发言只代表其个人观点，不代表新浪网的观点或立场。晶格常数 具体指什么 ？或者怎么算
以氯化钠为例 它指的是两个正离子或负离子之间的距离 当然是最近的正离子或负离子
年来R、ESR、NMR及Mssbauer谱等近代仪器方法也被用于配合物稳定常数测定的研究...（诚心为您解答，希望给予【好评】，感谢~）
⑤时间常数τ：热敏电阻器是有热惯性的，时间常数，就是一个描述热敏电阻器热惯性的参数
一个数学常数是指一个数值不变的常量，与之相反的是变量。跟大多数物理常数不一样的地方是，数学常数的定义是独立于所有物理测量的。
数学常数通常是实数或复数域的...
刀锋山最北边的那块图，西边有个叫巴什么的码头，里面全是虚空人，杀他们掉一个能够给你传送到另一个时空的东西，使用之后有点像指环王里带上魔戒的感觉，周围的东西都还在...
就是常数的运算形成的非线性的值域。如果您采纳的话或者觉得正确，麻烦给我好评哦，谢谢。
大家还关注
确定举报此问题
举报原因(必选)：
广告或垃圾信息
激进时政或意识形态话题
不雅词句或人身攻击
侵犯他人隐私
其它违法和不良信息
报告，这不是个问题
报告原因(必选)：
这不是个问题
这个问题分类似乎错了
这个不是我熟悉的地区博客访问： 149044
博文数量： 23
博客积分： 1410
博客等级： 上尉
技术积分： 261
注册时间：
IT168企业级官微
微信号：IT168qiye
系统架构师大会
微信号：SACC2013
分类： LINUX
我们用Pd金属作为例子。&&Pd金属的实验上的晶格常数为3.89A。在这里，我们用VASP计算它的晶格常数。&&首先将Pd所对应的POTCAR文件拷贝到目录下。然后准备好INCAR和KPOINTS文件。POSCAR文件我们将通过一个tcsh的script来产生。&&KPOINTS文件可以如下：&&Monkhorst Pack&&0&&Monkhorst Pack&& 11 11 11&& 0 0 0&&INCAR文件可以如下：&& SYSTEM = Pd bulk calculation&& Startparameter for this run:&& PREC = Accurate&& ISTART = 0 job : 0-new 1-cont 2-samecut&& ICHARG = 2 charge: 1-file 2-atom 10-const&& ISPIN = 1 spin polarized calculation?&& Electronic Relaxation 1&& EDIFF = 0.1E-03 stopping-criterion for ELM&& LREAL = .FALSE. real-space projection&& Ionic relaxation&& EDIFFG = 0.1E-02 stopping-criterion for IOM&& NSW = 0 number of steps for IOM&& IBRION = 2 ionic relax: 0-MD 1-quasi-New 2-CG&& ISIF = 2 stress and relaxation&& POTIM = 0.10 time-step for ionic-motion&& TEIN = 0.0 initial temperature&& TEBEG = 0.0; TEEND = 0.0 temperature during run&& DOS related values:&& ISMEAR = 0 ; SIGMA = 0.05 gaussian smear&& Electronic relaxation 2 (details)&& Write flags&& LWAVE = F write WAVECAR&& LCHARG = F write CHGCAR&&&&产生POSCAR和计算晶格常数的工作可以用以下的PBS script来完成。&&#!/bin/tcsh&&#PBS -S /bin/sh&&#PBS -l nodes=4:athlon:ppn=2&&#PBS -l cput=384:00:00&&#PBS -m ae&&#PBS -o output&&#PBS -e error.log&&# set parameter&&set EXEC = 'vasp'&&set SRC = '/usr/common/executable'&&# change working directory&&cd $PBS_O_WORKDIR&&# copy fresh executable from depository&&cp -f $SRC/$EXEC .&&# execute mpi program&&foreach a (3.3 3.4 3.5 3.6 3.7)&&echo "a= $a"&&cat >POSCAR <
阅读(4388) | 评论(1) | 转发(0) |
相关热门文章
给主人留下些什么吧！~~
不错 很实用 呵呵
请登录后评论。造价与管理
&&&&&&&&&&&&金属的晶格常数
金属的晶格常数
金属的常数一般为:&&&&&&1×10-10&m～7×10-10&m。&&&
&&&&不同元素组成的金属因晶格形式及晶格常数的不同,表现出不同的物理、化学和力学性能。金属的晶体结构可用X射线结构分析技术进行测定。
注：评论内容不得超过140个字&&
知识点记录
浏览次数：
编辑次数：1次&
最近更新：
[金属材料]热门知识君，已阅读到文档的结尾了呢~~
晶格常数的精确测定
扫扫二维码，随身浏览文档
手机或平板扫扫即可继续访问
晶格常数的精确测定
举报该文档为侵权文档。
举报该文档含有违规或不良信息。
反馈该文档无法正常浏览。
举报该文档为重复文档。
推荐理由：
将文档分享至：
分享完整地址
文档地址：
粘贴到BBS或博客
flash地址：
支持嵌入FLASH地址的网站使用
html代码：
&embed src='/DocinViewer-4.swf' width='100%' height='600' type=application/x-shockwave-flash ALLOWFULLSCREEN='true' ALLOWSCRIPTACCESS='always'&&/embed&
450px*300px480px*400px650px*490px
支持嵌入HTML代码的网站使用
您的内容已经提交成功
您所提交的内容需要审核后才能发布，请您等待！
3秒自动关闭窗口}