应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » 其他 » 【其他】在线答疑:VASP的问题
594/6返回列表上一页123456下一页
查看: 4832  |  回复: 58
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

wuchenwf

荣誉版主 (职业作家)
★ ★
aylayl08(金币+2,VIP+0):谢谢提示，新年快乐 1-2 16:26
引用回帖:
Originally posted by wo8165827 at 2009-12-30 13:21:
请问下磁性计算和态度计算的过程是不是一样的？

也就是说先要结构优化，然后静态计算，最后是磁性计算和态密度计算。

对，不过结构优化中你就需要设置ISPIN,MAGMOM
赞一下(10人) 回复此楼
31楼2010-01-01 21:23:57
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

lishufei

金虫 (正式写手)

h原子poscar问题

★
小木虫(金币+0.5):给个红包，谢谢回帖交流
我在一个文件里看到H原子poscar
atom 1
     15.00000    .00000    .00000
       .00000  15.00000    .00000
       .00000    .00000  15.00000
   1
cart
0    0    0
为什么是这个呀，清高手帮助解释一下
赞一下(1人) 回复此楼
32楼2010-01-03 14:26:58
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

lishufei

金虫 (正式写手)

最近刚学VASP，有一点问题想请教各位，

★
小木虫(金币+0.5):给个红包，谢谢回帖交流
我对於VASP OUTCAR中输出的的能量项不太了解，想请各位帮我解释一下


举例来说我对氢原子做了计算，得到从OUTCAR 中得到下列数据

   Free energy of the ion-electron system (eV)
   ---------------------------------------------------
   alpha Z PSCENC = 0.00621465
   Ewald energy TEWEN = -1.38027565
   -1/2 Hartree DENC = -6.91107031
   -V(xc)+E(xc) XCENC = 2.32601856
   PAW double counting = 0.70286470 -0.71256934
   entropy T*S EENTRO = 0.00434843
   eigenvalues EBANDS = -7.71257941
   atomic energy EATOM = 12.52153358
   ---------------------------------------------------
   free energy TOTEN = -1.15551480 eV

   energy without entropy = -1.15986323 energy(sigma->0) = -1.15696428

在我的概念中所谓的总能应该是：核动能+电子动能+核-核斥能+电子-电子斥能+核-电子位能＝系统总能

可是VASP给的数据我看不太懂，这边的free energy TOTEN，代表什麽？直接翻译的话是自由能，但是是指何种自由能？

有人说free energy TOTEN 这项就是总能了，可是0K 时H原子基态能量不是-13.6ev吗？

还有atomic energy EATOM 、eigenvalues EBANDS、ntropy T*S EENTRO、  -V(xc)+E(xc) XCENC、1/2 Hartree DENC、Ewald energy TEWEN、alpha Z PSCENC

这几项代表什麽的能量？

我之前的使用的计算软体GAUSSIAN都是分别给出：核动能、电子动能、核-核斥能、电子-电子斥能、核-电子吸引能、总能、交换能、关系能等，很一目了然。

可是VASP的我真得看不太懂，

想请教各位
赞一下(1人) 回复此楼
33楼2010-01-03 15:12:53
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

yunml

铜虫 (小有名气)

请教

★
小木虫(金币+0.5):给个红包，谢谢回帖交流
如何在vasp输出结果中看费米能级值。 我想看通过DOSCAR的数据算d band center.

谢谢！
赞一下(1人) 回复此楼
一以贯之
34楼2010-01-03 17:48:47
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

yangfangwi

新虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
你好，我计算zn单原子的能量，魇势选用的是gga，可是计算的结果是正值.这个怎么解决？用lda魇势结果为负。但是我总觉的gga也应该是正确的啊！不知道是不是自己哪弄错啦。。。
poscar
ZN 1
20
1       0        0
0      1        0
0      0        1
1
Direct
0 0 0
incar
SYSTEM=zn
ENCUT=350
NELMDL=5
ISMEAR=0
SIGMA=0.1
k
1 1 1
谢谢啦！
赞一下(1人) 回复此楼
35楼2010-01-04 20:21:27
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

lishufei

金虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
k-point   1:       0.0000  0.0000    0.0000
band No. band energies   occupation
赞一下(1人) 回复此楼
36楼2010-01-05 09:37:21
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

lishufei

金虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
Steepest descent step on ions:
trial-energy change:    0.229038  1 .order    0.190139   -0.112971    0.493250
  (g-gl).g = 0.113E+00      g.g   = 0.113E+00  gl.gl    = 0.000E+00
g(Force)  = 0.113E+00   g(Stress)= 0.000E+00 ortho     = 0.000E+00
gamma     =   0.00000
trial     =   1.00000
opt step  =   0.14000  (harmonic =   0.18635) maximal distance =0.00730821
next E    =    -9.861023   (d E  =  -0.00780)
Steepest descent step on ions:
trial-energy change:    0.229038  1 .order    0.190139   -0.112971    0.493250
  (g-gl).g = 0.113E+00      g.g   = 0.113E+00  gl.gl    = 0.000E+00
g(Force)  = 0.113E+00   g(Stress)= 0.000E+00 ortho     = 0.000E+00
gamma     =   0.00000
trial     =   1.00000
opt step  =   0.14000  (harmonic =   0.18635) maximal distance =0.00730821
next E    =    -9.861023   (d E  =  -0.00780)


我在算o原子二聚体时是什么意思阿
赞一下(1人) 回复此楼
37楼2010-01-05 12:53:22
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

lishufei

金虫 (正式写手)

我看co频率计算有一段这样写是什么意思啦？请教前辈

★
小木虫(金币+0.5):给个红包，谢谢回帖交流
Finite differences progress:
  Degree of freedom:   2/  2
  Displacement:        2/  2
  Total:               4/  4

FORCES
------

INITIAL FORCE
  0.000000   0.000000  -0.066860
  0.000000   0.000000   0.066860
DOF:   1 ATOM:   1 AXIS:   3 DISPLACEMENT:   1
  0.000000   0.000000  -2.526029
  0.000000   0.000000   2.526029

DOF:   1 ATOM:   1 AXIS:   3 DISPLACEMENT:   2
  0.000000   0.000000   2.069092
  0.000000   0.000000  -2.069092

DOF:   2 ATOM:   2 AXIS:   3 DISPLACEMENT:   1
  0.000000   0.000000   2.062261
  0.000000   0.000000  -2.062261

DOF:   2 ATOM:   2 AXIS:   3 DISPLACEMENT:   2
  0.000000   0.000000  -2.511030
  0.000000   0.000000   2.511030

  
DYNMAT
------
    2    1    1    2
16.000 12.011
    1    1  0.0000  0.0000  0.0200
  0.000000   0.000000  -2.297561
  0.000000   0.000000   2.297561
    2    1  0.0000  0.0000  0.0200
  0.000000   0.000000   2.286646
  0.000000   0.000000  -2.286646
  

SECOND DERIVATIVES (NOT SYMMETRIZED)
------------------------------------
               1Z          2Z
  1Z  -114.878038  114.878038
  2Z   114.332280 -114.332280


SECOND DERIVATIVES (SYMMETRYZED)
--------------------------------
               1Z          2Z
  1Z  -114.878038  114.605159
  2Z   114.605159 -114.332280


MASS-WEIGHTED SECOND DERIVATIVES
--------------------------------
               1Z          2Z
  1Z    -7.179877    8.267127
  2Z     8.267127   -9.518964


Eigenvectors and eigenvalues of the dynamical matrix
----------------------------------------------------


   1 f  =   63.884308 THz   401.396943 2PiTHz 2130.951061 cm-1   264.204354 meV
             X         Y         Z           dx          dy          dz
      0.000000  0.000000  0.000000            0           0   -0.655715  
      0.000000  0.000000  1.143000            0           0    0.755008  

   2 f/i=    0.075306 THz     0.473159 2PiTHz    2.511925 cm-1     0.311439 meV
             X         Y         Z           dx          dy          dz
      0.000000  0.000000  0.000000            0           0   -0.755008  
      0.000000  0.000000  1.143000            0           0   -0.655715  

Eigenvectors after division by SQRT(mass)

Eigenvectors and eigenvalues of the dynamical matrix
----------------------------------------------------


   1 f  =   63.884308 THz   401.396943 2PiTHz 2130.951061 cm-1   264.204354 meV
             X         Y         Z           dx          dy          dz
      0.000000  0.000000  0.000000            0           0   -0.163929  
      0.000000  0.000000  1.143000            0           0    0.217852  

   2 f/i=    0.075306 THz     0.473159 2PiTHz    2.511925 cm-1     0.311439 meV
             X         Y     





    Z           dx          dy          dz
      0.000000  0.000000  0.000000            0           0   -0.188752  
      0.000000  0.000000  1.143000            0           0   -0.189202  

Finite differences POTIM=  2.000000000000000E-002
  LATTYP: Found a simple cubic cell.
ALAT       =     8.0000000000
  
  Lattice vectors:
  
A1 = (   8.0000000000,   0.0000000000,   0.0000000000)
A2 = (   0.0000000000,   8.0000000000,   0.0000000000)
A3 = (   0.0000000000,   0.0000000000,   8.0000000000)
Subroutine PRICEL returns:
Original cell was already a primitive cell.


Analysis of symmetry for initial positions (statically):

Routine SETGRP: Setting up the symmetry group for a
simple cubic supercell.
赞一下(1人) 回复此楼
38楼2010-01-05 22:33:24
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

wuchenwf

荣誉版主 (职业作家)
★
aylayl08(金币+1,VIP+0):谢谢解答 1-10 09:07
引用回帖:
Originally posted by lishufei at 2010-1-3 14:26:
我在一个文件里看到H原子poscar
atom 1
     15.00000    .00000    .00000
       .00000  15.00000    .00000
       .00000    .00000  15.00000
   1
cart
0    0    0
为什么是这个呀，清高手帮助解 ...

这个是计算单原子的能量，15A是个大盒子。
赞一下(10人) 回复此楼
39楼2010-01-09 22:53:06
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

wuchenwf

荣誉版主 (职业作家)
★
aylayl08(金币+1,VIP+0):谢谢解答 1-10 09:07
引用回帖:
Originally posted by yunml at 2010-1-3 17:48:
如何在vasp输出结果中看费米能级值。 我想看通过DOSCAR的数据算d band center.

谢谢！

在 outcar最后找，里面说的挺清楚的。
或者grep fermi “OUTCAR” 也可以的
赞一下(10人) 回复此楼
40楼2010-01-09 22:54:30
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wuchenwf 的主题更新
594/6返回列表上一页123456下一页
信息提示
关闭
请填处理意见
关闭