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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » 其他 » 【其他】在线答疑:VASP的问题
594/6返回列表上一页123456下一页
查看: 4931  |  回复: 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下一页
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 297求调剂 +7 戏精丹丹丹 2026-03-17 7/350 2026-03-18 14:20 by 星空星月
[考研] 一志愿985，本科211，0817化学工程与技术319求调剂 +6 Liwangman 2026-03-15 6/300 2026-03-18 13:21 by 尽舜尧1
[考研] 302求调剂 +10 呼呼呼。。。。 2026-03-17 10/500 2026-03-18 12:45 by Linda Hu
[考研] 299求调剂 +5 △小透明* 2026-03-17 5/250 2026-03-18 11:49 by 尽舜尧1
[考研] 工科材料085601 279求调剂 +6 困于星晨 2026-03-17 6/300 2026-03-18 10:21 by kkcoco25
[考研] 268求调剂 +8 一定有学上- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[考研] 302求调剂 +9 负心者当诛 2026-03-11 9/450 2026-03-17 17:13 by ruiyingmiao
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 26考研求调剂 +6 丶宏Sir 2026-03-13 6/300 2026-03-17 16:13 by 醉在风里
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
[考研] 一志愿，福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 321求调剂 +5 大米饭！ 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
[考研] 招收0805（材料）调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
[硕博家园] 085600 260分求调剂 +3 天空还下雨么 2026-03-13 5/250 2026-03-13 18:46 by 天空还下雨么
[考研] 307求调剂 +5 超级伊昂大王 2026-03-12 5/250 2026-03-13 15:56 by 棒棒球手
信息提示
关闭
请填处理意见
关闭