| 查看: 992 | 回复: 7 | ||
xj20120828木虫 (正式写手)
|
[求助]
频率验证 已有1人参与
|
|
大家好，我用vasp计算完过渡态想验证一下频率，设置 PREC=Normal ISTART=0 ICHARG=2 ISPIN=1 LREAL=AUTO LWAVE=.FALSE. LCHARG=.FALSE. ISYM=0 Electronic Relaxation ENCUT=500 NELM=60 EDIFF=0.1E-3 #NELMIN=4 GGA=PE VOSKOWN=1 Ionic relaxtion EDIFFG=-.1E-0 NSW=1 IBRION=5 ISIF=2 POTIM=0.2 DOS related values: ISMEAR=0 !Gaussian SIGMA=0.02 Magnetism #MAGMOM=8*0 18*0 18*0 3*2 LORBIT=10 #NUPDOWN=2 #For NEB calculation #IMAGES=4 #LCLIMB=.True. #SPRING=-5 #ICHAIN=0 #LSCALAPACK = .FALSE. 计算频率，为什么设置NSW=1，还是会进行许多步的离子位置优化？不是应该只有一步吗？ |
回复此楼» 猜你喜欢
【求调剂】085601材料工程专硕 | 总分272 |
已经有5人回复
-
275求调剂
已经有13人回复
-
0856，材料与化工321分求调剂
已经有8人回复
-
085701环境工程，267求调剂
已经有12人回复
-
320分，材料与化工专业，求调剂
已经有7人回复
-
085701环境工程求调剂
已经有8人回复
-
295材料工程专硕求调剂
已经有5人回复
-
265求调剂
已经有8人回复
-
322求调剂
已经有4人回复
-
085602 307分 求调剂
已经有6人回复
xj20120828
木虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 2320.6
- 散金: 476
- 红花: 1
- 帖子: 511
- 在线: 195.4小时
- 虫号: 2205360
- 注册: 2012-12-25
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
2楼2014-08-13 10:01:48
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
xj20120828: 金币+20, ★★★★★最佳答案 2014-08-14 15:40:50
xj20120828: 金币+20, ★★★★★最佳答案 2014-08-14 15:40:50
|
For IBRION=5, all atoms are displaced in all three Cartesian directions, resulting in a significant computational effort even for moderately sized high symmetry systems. 因為你開了 IBRION=5 VASP 會對每一顆"可允許移動"的原子 計算 x y z 三個方向移動 POTIM 距離後的受力 所以 你有 N 顆 "允許移動" 的原子 就會有 3N 個離子步計算 (非優化) |
3楼2014-08-13 15:48:56
xj20120828
木虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 2320.6
- 散金: 476
- 红花: 1
- 帖子: 511
- 在线: 195.4小时
- 虫号: 2205360
- 注册: 2012-12-25
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
4楼2014-08-14 15:40:13
5楼2014-08-14 15:53:08
xj20120828
木虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 2320.6
- 散金: 476
- 红花: 1
- 帖子: 511
- 在线: 195.4小时
- 虫号: 2205360
- 注册: 2012-12-25
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
|
EDIFFG要那么高的精度吗？还是EDIFF? 我优化的是EDIFFG=-0.01。我算出来四个虚频， 171 f/i= 0.389539 THz 2.447547 2PiTHz 12.993627 cm-1 1.611005 meV 172 f/i= 1.438074 THz 9.035684 2PiTHz 47.968978 cm-1 5.947397 meV 173 f/i= 2.256721 THz 14.179397 2PiTHz 75.276111 cm-1 9.333052 meV 174 f/i= 22.875331 THz 143.729946 2PiTHz 763.038898 cm-1 94.604800 meV 68个原子。请问大概多少频率就算比较大的虚频了？ |
6楼2014-08-14 16:22:11
|
??o?ε? PHONON ????????@?N??? EDIFFG EDIFFG=0.01 ????????????????????? ???????M??PHONON??????? ?????l????? ???е? 171, 172, 173 ????"????"???? 0 ??????????????l?????? 174 ?? (~94 meV) 68??????????????????????????????? > ??l??С?????????o?P ?f???????"?^??"??l...??~?????? ??????X 10 meV ????????????? 0 ???? ????l??С??l?g??????~??С????l"???"??????????l ????N??????????????l (????]??L??^??????????????????) ?z????l?? vibration mode (no. 174) ?? OUTCAR ??????????? frequency ???????? phonon mode ?? POSCAR????????l?? phonon mode ?????????????С??λ?? (???С??? ??? 0.02 ??????????????????????С.....) ??????????? ??r???? EDIFFG ? 1E-04 ????????????? phonon ??? ???ЙC????????l????? |
» 本帖已获得的红花（最新10朵）
xj201208287楼2014-08-14 17:05:37
xj20120828
木虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 2320.6
- 散金: 476
- 红花: 1
- 帖子: 511
- 在线: 195.4小时
- 虫号: 2205360
- 注册: 2012-12-25
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
8楼2014-08-14 19:45:38
20