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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » 其他 » 【其他】在线答疑:VASP的问题
594/6返回列表上一页123456下一页
查看: 4965  |  回复: 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下一页
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料与化工300求调剂 +34 肖开文 2026-04-09 37/1850 2026-04-10 18:58 by HPUCZ
[考研] 一志愿北理工298英一数二求调剂 +13 Reframe 2026-04-10 15/750 2026-04-10 17:25 by laoshidan
[考研] 药学专硕调剂 +5 ? 一路生?花? 2026-04-10 6/300 2026-04-10 11:54 by pengliang8036
[考研] 269求调剂 +10 啊啊我我 2026-04-07 10/500 2026-04-10 10:38 by 高维春
[考研] 一志愿鲁东大学071000生物学学硕初试分数276求调剂 +20 慕绝cc 2026-04-09 24/1200 2026-04-10 10:22 by xujun0624
[考研] 青岛科技大学材料学院，环境学院调剂补录4月10日以前都可以 +3 1青科大。 2026-04-09 5/250 2026-04-10 09:58 by 翩翩一书生
[考研] 08600生物与医药-327 +10 18755400796 2026-04-05 10/500 2026-04-10 08:14 by kangsm
[考研] 考研调剂-材料类-284 +28 想换手机不想解� 2026-04-08 28/1400 2026-04-09 20:08 by 倒数321?
[考研] 085600材料与化工，求调剂 +11 won_qii 2026-04-07 11/550 2026-04-09 17:03 by luoyongfeng
[考研] 085400电子信息类（川大控制工程）求调剂可跨专业 求老师联系 +3 626776879 2026-04-08 3/150 2026-04-09 16:05 by 猪会飞
[考研] 308求调剂 +17 墨墨漠 2026-04-06 17/850 2026-04-09 09:25 by 壹往無前
[考研] 材料调剂 +14 一样YWY 2026-04-06 14/700 2026-04-08 23:00 by 猪会飞
[考研] 一志愿南昌大学，085600，344分求调剂 +11 调剂上岸玘 2026-04-05 12/600 2026-04-08 16:17 by luoyongfeng
[考研] 化工学硕 285求调剂 +26 Wisjxn 2026-04-07 26/1300 2026-04-08 14:42 by screening
[考研] 机械工程264学硕求调剂 +3 qiushangxian 2026-04-06 3/150 2026-04-08 01:53 by Linzejun
[考研] 085404 293求调剂 +8 勇远库爱314 2026-04-06 9/450 2026-04-07 13:05 by flydream1314
[考研] 324求调剂 +3 k可乐 2026-04-05 4/200 2026-04-06 09:54 by 蓝云思雨
[考研] 材料专硕322分 +10 哈哈哈吼吼吼哈 2026-04-04 10/500 2026-04-05 21:22 by 学员8dgXkO
[考研] 313求调剂 +5 海日海日 2026-04-04 7/350 2026-04-05 13:58 by imissbao
[考研] 328分调剂 +6 门men 2026-04-04 6/300 2026-04-05 13:40 by imissbao
信息提示
关闭
请填处理意见
关闭