24小时热门版块排行榜    

查看: 2433  |  回复: 18

郭芳8888

金虫 (正式写手)

[求助] VSAP计算中出现这个提示是什么问题啊,求高手 已有2人参与

running on    4 nodes
distr:  one band on    4 nodes,    1 groups
vasp.4.6.36 17Feb09 complex
POSCAR found :  2 types and   16 ions
LDA part: xc-table for Ceperly-Alder, Vosko type interpolation para-ferro

-----------------------------------------------------------------------------
|                                                                             |
|           W    W    AA    RRRRR   N    N  II  N    N   GGGG   !!!           |
|           W    W   A  A   R    R  NN   N  II  NN   N  G    G  !!!           |
|           W    W  A    A  R    R  N N  N  II  N N  N  G       !!!           |
|           W WW W  AAAAAA  RRRRR   N  N N  II  N  N N  G  GGG   !            |
|           WW  WW  A    A  R   R   N   NN  II  N   NN  G    G                |
|           W    W  A    A  R    R  N    N  II  N    N   GGGG   !!!           |
|                                                                             |
|      Full k-point grid generated                                            |
|      Inversion symmetry is not applied                                      |
|                                                                             |
-----------------------------------------------------------------------------

WARNING: stress and forces are not correct
POSCAR, INCAR and KPOINTS ok, starting setup
FFT: planning ...           1
reading WAVECAR
charge-density read from file: FePt                                    
magnetization density read from file 1
entering main loop
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1     0.120370466258E+04    0.12037E+04   -0.18002E+04******   0.223E+02
DAV:   2    -0.265987527721E+01   -0.12064E+04   -0.11437E+04******   0.205E+02
DAV:   3    -0.106976558636E+03   -0.10432E+03   -0.86964E+02******   0.587E+01
DAV:   4    -0.113317149720E+03   -0.63406E+01   -0.60954E+01******   0.152E+01
DAV:   5    -0.113500973539E+03   -0.18382E+00   -0.18336E+00******   0.224E+00
DAV:   6    -0.113508621144E+03   -0.76476E-02   -0.76448E-02******   0.438E-01
DAV:   7    -0.113508858557E+03   -0.23741E-03   -0.23736E-03******   0.765E-02
DAV:   8    -0.113508869363E+03   -0.10806E-04   -0.10803E-04789194   0.166E-02
   1 F= -.11350887E+03 E0= -.11350887E+03  d E =0.000000E+00  mag=     0.0000     0.0000    17.4810
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mywai520

铁杆木虫 (著名写手)


【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
郭芳8888: 金币+5, ★★★★★最佳答案, 谢谢大神 2014-04-22 08:30:33
你非自洽计算的参数跟自洽计算不同:具体看以下英文说明:
if the k-mesh is changed, the old charge density is of course not a self-consistent solution with respect to the new mesh, therefore the stress and forces are not necessarily correct if the run is not converged to self-consistency itself.
also, the read k-points from WAVECAR (and the respective wavefunction coefficients) may differ from the newly generated k-mesh (and hence be inconsistent)
2楼2014-04-17 19:36:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

李天帮

木虫 (正式写手)

小小木虫

【答案】应助回帖

感谢参与,应助指数 +1
楼主所用非自洽计算的的参数跟自洽计算所用的参数不同导致了这样问题啊,建议楼主看看手册啊.
事虽难,做则成;路虽远,行则至。
3楼2014-04-18 08:00:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

郭芳8888

金虫 (正式写手)

引用回帖:
2楼: Originally posted by mywai520 at 2014-04-17 19:36:08
你非自洽计算的参数跟自洽计算不同:具体看以下英文说明:
if the k-mesh is changed, the old charge density is of course not a self-consistent solution with respect to the new mesh, therefore the stress ...

自洽就是共线计算,非自洽就是非共线计算么?我的导师只是说共线与非共线。共线时ICHARG=2 非共线时ICHARG=11  。是这样的么?还有就是哪个参数不一样,不好意思,我是初学者。
4楼2014-04-18 16:48:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mywai520

铁杆木虫 (著名写手)


引用回帖:
4楼: Originally posted by 郭芳8888 at 2014-04-18 16:48:49
自洽就是共线计算,非自洽就是非共线计算么?我的导师只是说共线与非共线。共线时ICHARG=2 非共线时ICHARG=11  。是这样的么?还有就是哪个参数不一样,不好意思,我是初学者。...

你把所有的参数贴出来,才知道。
5楼2014-04-18 17:06:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

郭芳8888

金虫 (正式写手)

引用回帖:
5楼: Originally posted by mywai520 at 2014-04-18 17:06:23
你把所有的参数贴出来,才知道。...

System =  FePt

NPAR = 1
VOSKOWN = 1

ISPIN = 2 spin: 1 paramegnetic; 2 spin polarized

ISTART = 0 read WAVECAR: 0-no;  1-yes constant energy cutoff;  2-yes constant basis
             [default: 1-if WAVECAR exists, 0-else]

ICHARG = 2 initial charge density: 0-calc from initial wave functions; 1-from CHGCAR;
            2-from atomic charges; +10 non-self-consistent
            #[default: 2-if ISTART=0, 0-else]

LWAVE  = .FALSE.     print out WAVECAR (huge!)                          [default: .TRUE.]
#LCHARG = .FALSE.      print out CHGCAR  (large, but required to restart) [default: .TRUE.]

# NBANDS = 60 number of bands included (p105)
              [default: # NELECT/2+NIONS/2]
         
# ENCUT= 800
         
PREC   =  Accurate precision of calc: Normal, Accurate

MAGMOM = 3 3 3 3 0 0 0 0 0 0 0 0 0 0 0 0
#          initial atom magnetic moment [default: NIONS]

RWIGS  =  1.111 1.226  # FCC: 0.390796318028403*a (A)
           atomic radii for each species; to calc local partial DOS
           (e.g. 4/3*Pi*R_3=a_3/4, R=a*0.5*(1.5/Pi)_1/3 for FCC,
            look for small overlap or covalent radius, look for Sum of
            V_shells=V_unit_cell)
   
Electronic diagonalization:
#  IALGO = 48     algorithm: 5-8 conjugate gradient (8->38), 44-48 residual minimization
                            [default: 38] 48-large system
# LDIAG = .TRUE.  sub-space diagonalisation [default: .TRUE.]
#  LREAL =  Auto  real-space projection: .FALSE.-in reciprocal; On-in real;

                 Auto for big (>20 atoms) systems [default: .FALSE.]


  ISMEAR = -5  smearing:
                1-OK,2 Methfessel-Paxton for metals+relax;
                0 Gaussian (semicond+insulat with big unit cell);
               -1 Fermi; -2 read from WAVCAR; -4 tetrahedron;
               -5-OK tetrahedron with corrections
                  (Bill Butler; acc. for E,DOS in metals+NO-relax; semicond+insulat)
                [default: 1]

# NELM   = 100 max number of electr self-consist steps [default: 60]
# NELMIN =   2 min number of electr self-consist steps [default: 2]
# NELMDL = -12  number of electr non-self-consist steps at the begining
            [default: -5 if # ISTART=0, INIWAV=1, IALGO=8
                     -12 if # ISTART=0, INIWAV=1, IALGO=48
                       0 else]


Parameters related to DOS:
NEDOS  =     901    number of gridpoints for DOS
EMIN   =    -5.0    min energy [eV] for DOS calc
EMAX   =    10.0    max energy [eV] for DOS calc

---------------------------------------------------------------------------------------
             Ionic Relaxation (pp. 13,57,93,111; result is in CONTCAR)
---------------------------------------------------------------------------------------

NSW    = 0  number of ionic steps (0-default; 10-OK)

IBRION =  -1   ionic relax: (-1)-no ionic steps; 0-MD; 1-quasi-Newton;
                              2-conjugate gradient-OK; 3-damped MD; 5-Hessian?
                              [default: (-1)-if # NSW=0,1; 0-else]

ISIF   =  0  stress tensor: 0-no   +relax ions
                             1-trace+relax ions
                             2-yes  +relax ions             -OK
                             3-yes  +relax ions shape volume-OK
                             4-yes  +relax ions shape
                             5-yes  +relax      shape
                             6-yes  +relax      shape volume
                             7-yes  +relax            volume
                             [default: 0-if # IBRION=0 (MD); 2-else]

# ENMAX = 550  plane wave cut off    [default: see POTCAR; 550 for acc. relax]
EDIFF = 1e-7 error in total energy [default: 1e-4;      1e-7 for acc. relax]

# POTIM = 0.1 reduce first trial step in the direction of force at relax [default: 1.0]
EDIFFG = -0.1e-2 condition to stop relax: # >0-if dE<EDIFFG; <0-if dFirce<Abs(EDIFFG)-OK;
                                            #  0-stop after NSW steps
                                            #  [default: EDIFF*10]
# IWAVPR = 1 Extrapol. wave f. at relax: 0-no; 1-simple; 2-of 2nd order; 3-mixter not recommended
                                        #[default-OK: 2-if IBRION=0; 1-if IBRION=1,2; 0-else]




# INIWAV = 1 initial wave functions: 0-jelium; 1-random-OK
            #[default: 1]

ISYM   =  0    symmetry: 0-no; 1-yes; 2-yes memory conserving

LPETIM =  F    write-flag & timer:


# SIGMA  = 0.2 broadening width [eV]
                0.2 for metals [default: 0.2] Amrit: 0.1
6楼2014-04-18 18:01:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

郭芳8888

金虫 (正式写手)

引用回帖:
5楼: Originally posted by mywai520 at 2014-04-18 17:06:23
你把所有的参数贴出来,才知道。...

System =  FePt

NPAR = 1
VOSKOWN = 1

ISPIN = 2 spin: 1 paramegnetic; 2 spin polarized

ISTART = 0 read WAVECAR: 0-no;  1-yes constant energy cutoff;  2-yes constant basis
             [default: 1-if WAVECAR exists, 0-else]

ICHARG = 11 initial charge density: 0-calc from initial wave functions; 1-from CHGCAR;
            2-from atomic charges; +10 non-self-consistent
            #[default: 2-if ISTART=0, 0-else]

LWAVE  = .FALSE.     print out WAVECAR (huge!)                          [default: .TRUE.]
LCHARG = .FALSE.      print out CHGCAR  (large, but required to restart) [default: .TRUE.]
LSORBIT = .TRUE.
SAXIX = 1 0 0
LMAXMIX = 4
NBANDS = 194
GGA_COMPAT = .FALSE.
# NBANDS = 60 number of bands included (p105)
              [default: # NELECT/2+NIONS/2]

# ENCUT= 800

PREC   =  Accurate precision of calc: Normal, Accurate
   
MAGMOM = 0 0 3 0 0 3 0 0 3 0 0 3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
#          initial atom magnetic moment [default: NIONS]

RWIGS  =  1.111 1.226  # FCC: 0.390796318028403*a (A)
           atomic radii for each species; to calc local partial DOS
           (e.g. 4/3*Pi*R_3=a_3/4, R=a*0.5*(1.5/Pi)_1/3 for FCC,
            look for small overlap or covalent radius, look for Sum of
            V_shells=V_unit_cell)

Electronic diagonalization:
#  IALGO = 48     algorithm: 5-8 conjugate gradient (8->38), 44-48 residual minimization
                            [default: 38] 48-large system
# LDIAG = .TRUE.  sub-space diagonalisation [default: .TRUE.]
#  LREAL =  Auto  real-space projection: .FALSE.-in reciprocal; On-in real;
                 Auto for big (>20 atoms) systems [default: .FALSE.]


  ISMEAR = -5  smearing:
                1-OK,2 Methfessel-Paxton for metals+relax;
                0 Gaussian (semicond+insulat with big unit cell);
               -1 Fermi; -2 read from WAVCAR; -4 tetrahedron;
               -5-OK tetrahedron with corrections
                  (Bill Butler; acc. for E,DOS in metals+NO-relax; semicond+insulat)
                [default: 1]

# NELM   = 100 max number of electr self-consist steps [default: 60]
# NELMIN =   2 min number of electr self-consist steps [default: 2]
# NELMDL = -12  number of electr non-self-consist steps at the begining
            [default: -5 if # ISTART=0, INIWAV=1, IALGO=8
                     -12 if # ISTART=0, INIWAV=1, IALGO=48
                       0 else]


Parameters related to DOS:
NEDOS  =     901    number of gridpoints for DOS
EMIN   =    -5.0    min energy [eV] for DOS calc
EMAX   =    10.0    max energy [eV] for DOS calc

---------------------------------------------------------------------------------------
             Ionic Relaxation (pp. 13,57,93,111; result is in CONTCAR)
---------------------------------------------------------------------------------------

NSW    = 0  number of ionic steps (0-default; 10-OK)

IBRION =  -1   ionic relax: (-1)-no ionic steps; 0-MD; 1-quasi-Newton;
                              2-conjugate gradient-OK; 3-damped MD; 5-Hessian?
                              [default: (-1)-if # NSW=0,1; 0-else]

ISIF   =  0  stress tensor: 0-no   +relax ions
                             1-trace+relax ions
                             2-yes  +relax ions             -OK
                             3-yes  +relax ions shape volume-OK
                             4-yes  +relax ions shape
                             5-yes  +relax      shape
                             6-yes  +relax      shape volume
                             7-yes  +relax            volume
                             [default: 0-if # IBRION=0 (MD); 2-else]

# ENMAX = 550  plane wave cut off    [default: see POTCAR; 550 for acc. relax]
# EDIFF = 1e-7 error in total energy [default: 1e-4;      1e-7 for acc. relax]

# POTIM = 0.1 reduce first trial step in the direction of force at relax [default: 1.0]
# EDIFFG = -0.1e-2 condition to stop relax: # >0-if dE<EDIFFG; <0-if dFirce<Abs(EDIFFG)-OK;
                                            #  0-stop after NSW steps
                                            #  [default: EDIFF*10]
# IWAVPR = 1 Extrapol. wave f. at relax: 0-no; 1-simple; 2-of 2nd order; 3-mixter not recommended
                                        #[default-OK: 2-if IBRION=0; 1-if IBRION=1,2; 0-else]




# INIWAV = 1 initial wave functions: 0-jelium; 1-random-OK
            #[default: 1]

ISYM   =  0    symmetry: 0-no; 1-yes; 2-yes memory conserving

LPETIM =  F    write-flag & timer:


# SIGMA  = 0.2 broadening width [eV]
                0.2 for metals [default: 0.2] Amrit: 0.1
7楼2014-04-18 18:02:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mywai520

铁杆木虫 (著名写手)


【答案】应助回帖

你的输入卡真复杂啊。初看你的磁矩设置就不一样。
8楼2014-04-18 18:22:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

郭芳8888

金虫 (正式写手)

引用回帖:
8楼: Originally posted by mywai520 at 2014-04-18 18:22:51
你的输入卡真复杂啊。初看你的磁矩设置就不一样。

能不能说的再详细一点啊??谢谢
9楼2014-04-21 09:27:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mywai520

铁杆木虫 (著名写手)


【答案】应助回帖

自洽MAGMOM = 3 3 3 3 0 0 0 0 0 0 0 0 0 0 0 0
非自洽:MAGMOM = 0 0 3 0 0 3 0 0 3 0 0 3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

完全搞不懂,你的自洽就这么12个原子,怎么在非自洽计算中就一下子这么多原子需要设置了?
10楼2014-04-21 09:39:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 郭芳8888 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料工程一志愿中科大总分312求调剂 +4 吃宵夜1 2026-02-28 6/300 2026-02-28 17:38 by 吃宵夜1
[考研] 311求调剂 +6 南迦720 2026-02-28 6/300 2026-02-28 17:17 by 刘茉莉
[考研] 高分子化学与物理调剂 +3 好好好1233 2026-02-28 4/200 2026-02-28 17:13 by sandychj
[考研] 化工专硕348,一志愿985求调剂 +3 弗格个 2026-02-28 5/250 2026-02-28 17:04 by sandychj
[考博] 博士自荐 +3 kkluvs 2026-02-28 3/150 2026-02-28 16:59 by StarAura
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 11/550 2026-02-28 16:57 by ichall
[考博] 博士推荐 +4 花儿笑? 2026-02-21 5/250 2026-02-28 16:55 by mumin1990
[考研] 0856材料求调剂 +8 hyf hyf hyf 2026-02-28 9/450 2026-02-28 16:27 by etapple
[考研] 285求调剂 +3 满头大汗的学生 2026-02-28 3/150 2026-02-28 16:22 by 无际的草原
[考研] 265分求调剂不调专业和学校有行学上就 +4 礼堂丁真258 2026-02-28 6/300 2026-02-28 16:18 by 求调剂zz
[考博] 26申博 +3 想申博! 2026-02-26 3/150 2026-02-28 16:07 by nxgogo
[考研] 材料类求调剂 +3 wana_kiko 2026-02-28 3/150 2026-02-28 15:03 by lature00
[考研] 295求调剂 +4 19171856320 2026-02-28 4/200 2026-02-28 13:39 by ms629
[考研] 304求调剂 +5 曼殊2266 2026-02-28 6/300 2026-02-28 12:44 by 迷糊CCPs
[硕博家园] 博士自荐 +6 科研狗111 2026-02-26 9/450 2026-02-28 12:32 by seaskyy
[考研] 272求调剂 +3 田智友 2026-02-28 3/150 2026-02-28 12:31 by 王加浩to
[考研] 298求调剂 +4 axyz3 2026-02-28 4/200 2026-02-28 11:21 by wang_dand
[基金申请] 什么是人一生最重要的? +10 瞬息宇宙 2026-02-21 10/500 2026-02-27 08:46 by tfang
[教师之家] 版面费该交吗 +14 苹果在哪里 2026-02-22 17/850 2026-02-26 11:55 by AGI智创机器人
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
请填处理意见