| 查看: 3632 | 回复: 15 | ||||||
| 本帖产生 1 个 1ST强帖 ，点击这里进行查看 | ||||||
langren24新虫 (小有名气)
|
[求助]
spin-orbit计算求助！
|
|||||
|
最近用spin-orbit计算体系的磁性总是出错，现将计算过程中的一些问题和疑惑贴出来，求大虾解答： 1. 用spin-orbit与spin-polarizedINCAR文件的不同有哪些？除了需要加 LSORBIT=.TRUE. SAXIS =X Y Z NBANDS变为共线的两倍，MAGMON也变为xyz方向之外 需不需要加 GGA_COMPAT=.FALSE. ISYM =0 呢？（这是在VASP手册上看到的） 2.自旋量子轴SAXIS =X Y Z 应该怎么设置啊？应该与体系晶格结构有关，轴的设置应该根据实空间晶格设置呢，还是根据倒空间晶格设置呢，求举例！ 3.根据VASP手册上提供的方法，计算spin-orbit之前应该先计算共线情况下的WAVECAR和CHGCAR,然后再开始spin-orbit的计算。 但是我在计算的时候导入之前计算的WAVECAR和CHGCAR，计算总会出错，不知道是什么原因？是不是应该计算非共线情况下的WAVECAR和CHGCAR，然后再计算spin-orbit呢？ 4.计算spin-orbit需不需要使用超胞呢？（掺杂或者不掺杂的情况下） 使用超胞的计算结果该怎么分析呢？ 5.计算spin-orbit对赝势（譬如GGA-PBE等）有没有要求，对算法（譬如ALGO=Normal、ALGO=Damped等），对计算精度等有没有要求？ [ Last edited by langren24 on 2011-12-13 at 15:28 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 vasp计算资料 | first principle | 自旋轨道耦合 | First principle |
» 猜你喜欢
一志愿天大材料与化工（085600）总分338
已经有4人回复
-
085700资源与环境308求调剂
已经有3人回复
-
求材料调剂
已经有8人回复
-
294求调剂材料与化工专硕
已经有5人回复
-
一志愿华中科技大学，080502，354分求调剂
已经有4人回复
-
一志愿吉林大学材料学硕321求调剂
已经有6人回复
-
085410人工智能专硕317求调剂（0854都可以）
已经有3人回复
-
330求调剂
已经有3人回复
-
一志愿中海洋材料工程专硕330分求调剂
已经有5人回复
-
304求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- molpro计算自旋轨道耦合输入求助
已经有3人回复
- 求助高斯能量计算问题
已经有6人回复
- spin-orbit怎么用啊？
已经有6人回复
- 【求助】用Orca计算一个片簇的单点能，不收敛。
已经有7人回复
- 【求助】vasp静态计算出现的WARNING
已经有18人回复
- 【讨论】ADF能否z在考虑spin orbit相对论效应下计算吸收光谱
已经有8人回复
- 【求助】wien2k中K-point个数
已经有9人回复
- 【求助】算磁各向异性能
已经有4人回复
- 【求助】spin-orbit coupling
已经有16人回复
- 【求助成功】请教：castep计算结果中2*Integrated Spin Density的意义？
已经有10人回复
- 【求助】PBE0 计算Si、Ge能带结构
已经有9人回复
- 【求助】按照计算Dmol分子轨道的例子计算的结果不同？
已经有16人回复
- 【求助】G09能否用dft计算旋轨耦合
已经有5人回复
- 求助：翻译PRB中理论的一部分
已经有4人回复
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
youzhizhe(金币+5, 1ST强帖+1): 谢谢交流。授予EPI一枚。计入专家考核存档。 2011-12-13 20:30:15
langren24(金币+2): 2011-12-14 14:13:18
感谢参与，应助指数 +1
youzhizhe(金币+5, 1ST强帖+1): 谢谢交流。授予EPI一枚。计入专家考核存档。 2011-12-13 20:30:15
langren24(金币+2): 2011-12-14 14:13:18
|
1. ISYM =0在打开SO时要加，用GGA交换泛函时需要用GGA_COMPAT=.FALSE.改善精度。 2. SAXIS和MAGMOM有两种常用设置，手册上有介绍，默认是0 0 1方向，是相对实空间笛卡尔坐标系的。一种是SAXIX为0 0 1方向，此时MAGMOM为任意方向，另外一种是MAGMOM均为 0 0 z或0 0 -z方向，SAXIX为任意方向，第二种允许读取先前存在的共线或非共线电荷密度和波函数，且精度更高。 3. 导入之前计算的WAVECAR和CHGCAR出错原因是KPOINTS在SO打开时对称性发生变化。你可以用SO生成的IBZKPT,然后拷贝为KPOINTS,重新计算WAVECAR和CHGCAR，然后再读取就可以了。 4. 不需要用超胞！ 5. SO只支持PAW赝势，对精度要求高，能量精度EDIFF一般为10^-7eV. |
3楼2011-12-13 20:17:54
fzx2008
荣誉版主 (著名写手)
- 1ST强帖: 9
- 应助: 272 (大学生)
- 贵宾: 1.045
- 金币: 21814.1
- 散金: 710
- 红花: 33
- 沙发: 1
- 帖子: 1808
- 在线: 1601小时
- 虫号: 621341
- 注册: 2008-10-09
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 第一性原理
2楼2011-12-13 15:32:44
雪狼1210
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 919.7
- 散金: 43
- 帖子: 128
- 在线: 155.9小时
- 虫号: 939826
- 注册: 2010-01-09
- 专业: 凝聚态物性 II ：电子结构
4楼2011-12-13 21:02:05
spatial
荣誉版主 (职业作家)
- 应助: 1 (幼儿园)
- 贵宾: 17.95
- 金币: 2915.7
- 帖子: 3871
- 在线: 44.7小时
- 虫号: 46929
- 注册: 2004-05-27
- 专业: 凝聚态物性 II ：电子结构
5楼2011-12-14 12:41:05
langren24
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 127.5
- 散金: 11
- 帖子: 193
- 在线: 96.4小时
- 虫号: 1136177
- 注册: 2010-11-01
- 专业: 新型信息器件
6楼2011-12-14 14:24:52
langren24
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 127.5
- 散金: 11
- 帖子: 193
- 在线: 96.4小时
- 虫号: 1136177
- 注册: 2010-11-01
- 专业: 新型信息器件
7楼2011-12-14 14:52:02
8楼2011-12-14 16:08:28
雪狼1210
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 919.7
- 散金: 43
- 帖子: 128
- 在线: 155.9小时
- 虫号: 939826
- 注册: 2010-01-09
- 专业: 凝聚态物性 II ：电子结构
9楼2011-12-14 16:22:04
10楼2011-12-14 16:49:34
