| 查看: 1334 | 回复: 11 | ||
| 当前主题已经存档。 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
jugengfans金虫 (正式写手)
|
[交流]
【求助】为什么静态计算和原先的磁矩差异如此之大！
|
|
|
RT，在结构优化后我把NSW等注释掉，进行静态计算， 但是，为什么得出的磁矩根刚才计算时的出入非常大！ 结构优化时我的磁矩是4.1，静态计算后变成了6.7,这是怎么回事？ 恳请各位帮忙看一下，谢谢了！ 我作静态计算的的INCAR文件： SYSTEM = layer ISTART = 0 ICHARG = 2 ENCUT = 520 ISPIN = 2 #MAGMOM = 8*0.0 8*1.0 24*0.0 EDIFF = 1E-4 EDFIIG = -1E-2 ISMEAR = 2 NPAR = 4 #NELM=200 #NELMDL=8 #WEIMIN=0 AMIX = 0.2 BMIX = 0.00001 AMIX_MAG = 0.8 BMIX_MAG = 0.00001 LREAL = .TRUE. Dynamic: #IBRION = 2 #NSW = 500 #POTIM = 0.2 #LORBIT = 11 #LDAU = .TRUE. #LDAUTYPE = 2 #LDAUL = 2 -1 -1 #LDAUU = 3.00 0.00 0.00 #LDAUJ = 0.87 0.00 0.00 #LDAUPRINT = 2 LMAXMIX = 4 ! Important: mix paw occupancies up to L=4 |
回复此楼
六198
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1542.4
- 散金: 20
- 红花: 2
- 帖子: 697
- 在线: 92小时
- 虫号: 470500
- 注册: 2007-11-29
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2009-08-07 08:51:52
xirainbow
木虫 (正式写手)
- 1ST强帖: 2
- 应助: 0 (幼儿园)
- 贵宾: 0.03
- 金币: 2132.6
- 红花: 2
- 帖子: 775
- 在线: 6.1小时
- 虫号: 431211
- 注册: 2007-08-11
- 专业: 光谱学
4楼2009-08-07 08:54:20
nkleof
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2093.2
- 散金: 649
- 红花: 5
- 帖子: 813
- 在线: 165.9小时
- 虫号: 152581
- 注册: 2006-01-02
- 专业: 凝聚态物性 II ：电子结构
★ ★ ★
jugengfans(金币+1):谢谢参与
yjr(金币+2,VIP+0):谢谢您的热心！ 8-7 12:34
jugengfans(金币+1):谢谢参与
yjr(金币+2,VIP+0):谢谢您的热心！ 8-7 12:34
|
你在做静态计算的时候没有输入MAGMOM值，所以vasp默认按照铁磁序优化，有可能和你原来的磁结构不相同，而恰好体系在铁磁结构也有个稳定态，就得到了不同的结果。vasp说明里面提到可以在INCAR里设置ICHARG=1让程序从前面的CHGCAR文件读取电荷密度，同时也读取磁矩的信息，但是几乎很难以实用，因为改变了k点或者NBANDS等等参数后，可能造成读取不到合适的电荷密度，所以我一直都在静态计算的时候也设置MAGMOM值，这样和优化时得到的磁矩基本上相同。 |
5楼2009-08-07 09:49:48
jugengfans
金虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 773.9
- 散金: 151
- 红花: 3
- 帖子: 537
- 在线: 205.6小时
- 虫号: 291609
- 注册: 2006-10-29
- 专业: 凝聚态物性 II ：电子结构
6楼2009-08-07 12:35:01
