哪位同学帮忙查一下铑Rh的晶格常数~~谢谢啦:)
贵金属铑和铱的电子结构和物理性质Electronic Structures and Physical Properties of Rhodium and Iridium关键词：Rh,Ir,电子结构,晶格常数,结合能作者：彭坤,杨昕昕,谢佑卿概述：应用纯金属单原子(OA)理论确定了面心立方(fcc)结构贵金属Rh和Ir的电子结构分别为[Kr](4dn)1.9000(4dc)6.7400(5sc)0.0000(5sf)0.3600和[Xe](5dn)0.7100(5dc)7.2900(6sc)0.2700(6sf)0.7300,并计算了它们的势能曲线、晶格常数、结合能、弹性和膨胀系数随温度的变化.另外对其体心立方(bcc),密排六方(hcp)结构初态特征晶体和初态液体的电子结构进行了研究.从而提供了不同晶体结构Rh和Ir的结构参数和性能参数,为Rh和Ir金属元素及其相关材料的理论设计提供了完整的数据资料.真懒.
为您推荐：
其他类似问题
扫描下载二维码小木虫 --- 500万硕博科研人员喜爱的学术科研平台
&&查看话题
【求助】关于 final bulk modulus
大家好，求助一下有关final bulk modulus 是什么，我在一篇文章中看到的，说The clulated final bulk modulus is 280 GPa。在 vasp计算结果中怎么找到这个压强的数值？请高人指点
是在平衡状态对应的晶格常数左右附近取几个值，通过birch-Murnaghan equation of state 算出相应的体积与能量进行拟合吗？
研究生必备与500万研究生在线互动！
扫描下载送金币
浏览器进程
打开微信扫一扫
随时随地聊科研 上传我的文档
 下载
 收藏
该文档贡献者很忙，什么也没留下。
 下载此文档
正在努力加载中...
vasp弹性常数计算过程
下载积分：1500
内容提示：vasp弹性常数计算过程，vasp 静态计算，vasp计算能带，弹性常数，化学平衡常数计算，rc时间常数计算，耦合常数计算，平衡常数的计算，化学平衡常数的计算，晶格常数计算公式
文档格式：DOC|
浏览次数：53|
上传日期： 19:39:40|
文档星级：
该用户还上传了这些文档
vasp弹性常数计算过程
官方公共微信小木虫 --- 500万硕博科研人员喜爱的学术科研平台
&&查看话题
求助CASTEP中构建α-Ti (N)固溶体的晶体结构模型
现已知：N、O等轻原子在低温下在六方结构的过渡金属Ti 或Nb 或Zr或Hf中以八面体或四面体间隙位置存在，保持hexagonal -Ti晶格 （P 63/m m c，194）
& &（1）例：Ti0.83 N0.17&&ICSD #108614 钛氮固溶体
Atom OX&&SITE& &x& && && && && && &y& && && && && && &&&z& && & SOF& &
Ti 1& &+0&&2c& & 0.3333& && && &0.6667& && && && &0.25& &0.83
N 1&&+0&&2c& &&&0.3333& && && & 0.6667& && && && &0.25& & 0.17
请高手看MS中采用两种方式构建有无问题：
A: 在加原子时occupancy分别设为 SOF对应值 0.83 0.17，但提示同一位置不能加原子，如何解决？还有SOF与occupancy是一个概念吗？
B：采用VCA混合原子法，先构造Ti的晶格，composition设置为Ti 83% N17%,但计算出来的弹性常数值很小，好像不对，请高手分析？VCA方法适用范围？
****如何构建金属晶格中的间隙原子模型******** 超晶胞？？求八面体或四面体间隙在MS中坐标计算方法？
& &（2）例：TiO0.325&&钛氧固溶体
Atom #& & OX& && & SITE& && && &x& && && && & y& && && & z& && && & SOF
Ti 1& && && &&&+0& & 2c& && && & 0.3333& &&&0.& && & 1. 0
O 1& && && &&&+0& &2a& && && &&&0& && && && & 0& && && &&&0& && && && &0.48
Ti1位置应该就是原始Ti晶格的Ti位置，O的坐标是否为hexagonal -Ti在MS中的八面体间隙位置呢？
*************超晶胞构建的方法，如何安放N/O等间隙原子？***********
QQ图片32.jpg
非常感谢，
（1）那删除原子有什么规则吗？上图为2×2×2的超胞，共16个原子，比如构建Ti16O5时，这5个氧原子如何分布？
（2）几何优化时，是否加对称性，加完后对称性变了，不是六方了，这样就不满足固溶体晶格不变了。
（3）可否改变氧个数，比较形成能说明溶解度；
谢谢，补充一下，第一种情况中是指构造超胞后，在部分取代Ti，还是将N插入其T间隙？
另，第一个式子中n/ti为0.2，这与写成TiN0.2有什么区别，因为在findit中没有类似TiN0.2的写法，还有pdf卡片中TiN0.26，TiN0.32等又如何构造呢？
1 这是个很困难的问题，一般需要尝试多种可能性。
2 固溶体是个长程无序结构，计算只能用周期性体系，理论上需要非常大的超晶胞才能保持删除原子后有完全一致的对称性，但实际上很多性质不太依赖周期性，这一点取决于结构构造，也需要尝试。
3 可以，但如上所述正确构造结构十分困难，请阅读有关文献看看怎么做。
1 取代，因为结构明确表示Ti和N的位置是一样的。
2 只是习惯性写法不同。
所有小数的都需要用超晶胞删除/插入/替换来构造。
1。取代完后，几何优化提示加对称性，能提高速度，但这会改变原始的晶系，如六方变三斜等。是否加对称性。
2。几何优化选项中是否要选cell opt，优化完单胞后对称性也变了。
优化完后的晶系不是原来固溶体晶格结构了，这样正确吗？
1。取代完后，几何优化提示加对称性，能提高速度，但这会改变原始的晶系，如六方变三斜等。是否加对称性。
2。几何优化选项中是否要选cell opt，优化完单胞后对称性也变了。
优化完后的晶系不是原来固溶体晶格结构了，这样正确吗？
研究生必备与500万研究生在线互动！
扫描下载送金币
浏览器进程
打开微信扫一扫
随时随地聊科研小木虫 --- 500万硕博科研人员喜爱的学术科研平台
&&查看话题
关于用VASP计算弹性常数的两个问题
最近在学习用VASP 5.3.5计算弹性常数，设置IBRION = 6, ISIF = 4计算。发现有两个问题想请教大家：
1. 当我用16核计算时，设置NCORE = 4，走完第一个离子步后报错，停止计算：
这是否意味着要直接取消掉NCORE一类的相关参数？取消之后是warning，提示我设置NCORE。
2. 由于计算时晶胞有一定的变形，导致K点网格和晶格不匹配，所以往往产生以下的提示，但是会继续计算：
VERY BAD NEWS! internal error in subroutine IBZKPT:
Reciprocal lattice and k-lattice belong to different class of lattices. Often results are still useful...& && &96
不知道是否对结果有什么影响？
谢谢大家了。
非常感谢你的回答，您说的手动算确实是一个选择，我现在尝试的是面心立方的铝的弹性常数的计算，像这样高对称的结构如果用VASP自带的功能计算而出现我说的提示不知道是否要紧？因为它会继续计算。
另外，不太明白您说的“parallel的开关只能选择默认值”，好像NCORE不存在默认值吧？是NCORE =1的意思？
最后，非常感谢你附的参考文献。:hand:
有关parallelisation，主要是关系到multi-core multi-node的运算的效率
NCORE和NPAR是差不多的，INCAR内声明一个即可。NCORE默认值为1。parallelisation只会影响运算效率，一般不影响运算的结果。
但IBRION=5，6貌似都只能用默认值算(也就是你不在INCAR声明NPAR或者NCORE，或者把这些注释掉)，然后你用了默认值，输出文件里面肯定会提示你NPAR或者NCORE使用默认值计算的效率不高，此时忽略警告即可（因为你只能用默认值）。
还有，IBRION=6默认应该只取两个点（NFREE=2），加上strain=0的情况总共3个点，个人觉得不是很准确，因为拟合直线的话一般多取一些比较好（例如NFREE=4），然后调整好合适的strain的大小（POTIM默认为0.015Angst，根据你的体系调整）
parallelisation&&http://cms.mpi.univie.ac.at/vasp/vasp/Parallelisation_NPAR_NCORE_LPLANE_KPAR_tag.html
NCORE http://cms.mpi.univie.ac.at/wiki/index.php/NCORE
IBRION=6、NFREE、POTIM： http://cms.mpi.univie.ac.at/vasp/vasp/IBRION_5_IBRION_6.html
第一个问题我觉得不能设置并行参数，只能把NCORE注释掉；第二个暂时还没有什么好方法，而且最近算的弹性常数老是错的，全是负数:(
你好，我最近使用VASP计算D3h晶体的弹性常数。然后输出的结果不是很懂。
不知道如何从VASP的OUTCAR中获得 Cij, 请多指教。
TOTAL ELASTIC MODULI (kBar)
Direction& & XX& && && & YY& && && & ZZ& && && & XY& && && & YZ& && && & ZX
--------------------------------------------------------------------------------
XX& && && &614.9821& & 162.6485& &&&-1.1175& &&&-0.5210& && &0.0000& && &0.0000
YY& && && &162.6485& & 615.4451& &&&-1.0843& && &0.9220& && &0.0000& && &0.0000
ZZ& && && & -1.1175& &&&-1.0843& &&&-1.7237& && &0.0287& && &0.0000& && &0.0000
XY& && && & -0.5210& && &0.9220& && &0.0287& & 227.1156& && &0.0000& && &0.0000
YZ& && && &&&0.0000& && &0.0000& && &0.0000& && &0.0000& &&&-0.4616& && &0.0064
ZX& && && &&&0.0000& && &0.0000& && &0.0000& && &0.0000& && &0.0064& &&&-0.4690
另外如何根据晶体对称性来确定Cij中独立变量，
如何从VASP的OUTCAR文件中获得Cij数值？
还请不吝指教。
另外你有即时沟通工具吗？ 有的话，方便提供一下吗？
你好，我最近使用VASP计算D3h晶体的弹性常数。然后输出的结果不是很懂。
不知道如何从VASP的OUTCAR中获得 Cij, 请多指教。
TOTAL ELASTIC MODULI (kBar)
Direction& & XX& && && & YY& && && & ZZ& && && & XY& && && & YZ& && && & ZX
--------------------------------------------------------------------------------
XX& && && &614.9821& & 162.6485& &&&-1.1175& &&&-0.5210& && &0.0000& && &0.0000
YY& && && &162.6485& & 615.4451& &&&-1.0843& && &0.9220& && &0.0000& && &0.0000
ZZ& && && & -1.1175& &&&-1.0843& &&&-1.7237& && &0.0287& && &0.0000& && &0.0000
XY& && && & -0.5210& && &0.9220& && &0.0287& & 227.1156& && &0.0000& && &0.0000
YZ& && && &&&0.0000& && &0.0000& && &0.0000& && &0.0000& &&&-0.4616& && &0.0064
ZX& && && &&&0.0000& && &0.0000& && &0.0000& && &0.0000& && &0.0064& &&&-0.4690
另外如何根据晶体对称性来确定Cij中独立变量，
如何从VASP的OUTCAR文件中获得Cij数值？
还请不吝指教。
另外你有即时沟通工具吗？ 有的话，方便提供一下吗？
你好，我最近使用VASP计算D3h晶体的弹性常数。然后输出的结果不是很懂。
不知道如何从VASP的OUTCAR中获得 Cij, 请多指教。
TOTAL ELASTIC MODULI (kBar)
Direction& & XX& && && & YY& && && & ZZ& && && & XY& && && & YZ& && && & ZX
--------------------------------------------------------------------------------
XX& && && &614.9821& & 162.6485& &&&-1.1175& &&&-0.5210& && &0.0000& && &0.0000
YY& && && &162.6485& & 615.4451& &&&-1.0843& && &0.9220& && &0.0000& && &0.0000
ZZ& && && & -1.1175& &&&-1.0843& &&&-1.7237& && &0.0287& && &0.0000& && &0.0000
XY& && && & -0.5210& && &0.9220& && &0.0287& & 227.1156& && &0.0000& && &0.0000
YZ& && && &&&0.0000& && &0.0000& && &0.0000& && &0.0000& &&&-0.4616& && &0.0064
ZX& && && &&&0.0000& && &0.0000& && &0.0000& && &0.0000& && &0.0064& &&&-0.4690
另外如何根据晶体对称性来确定Cij中独立变量，
如何从VASP的OUTCAR文件中获得Cij数值？
还请不吝指教。
另外你有即时沟通工具吗？ 有的话，方便提供一下吗？
确定Cij的独立变量看文献就行，或者弹塑性力学上有，这个还是比较容易查找到的。
OUTCAR中给出的是一个二维六阶矩阵，XX-XX就是C11，XX-YY是C22，以此类推，不同的是XY-XY对应的是C66，我看你在另一个帖子中回复了，相信你也已经知道啦。
QQ什么的还是算啦，因为我自己就是一个菜鸟。
你贴出来的这个TOTAL ELATIC MODULI就是弹性常数了啊，注意XY=6 XZ=5 YZ=4 就行了
一般对称性越低 独立变量越多，具体看 J. F. Nye Physical Properties of Crystals&&P140-141 的图表有详细的说明
研究生必备与500万研究生在线互动！
扫描下载送金币
浏览器进程
打开微信扫一扫
随时随地聊科研