weiwei128

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1432.8
散金: 29
帖子: 53
在線: 17.4小時
蟲號: 468243
注冊: 2007-11-27
性別: MM
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[求助] 分子動力學(xué)計算求助

版友們好，我嘗試用MD實現(xiàn)模擬退火，計算體系結(jié)構(gòu)變化。遇到了些問題，提示如下：
-----------------------------------------------------------------------------
|                                                                            |
|          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 !!!          |
|                                                                            |
|    VASP found 315 degrees of freedom                                  |
|    the temperature will equal 2*E(kin)/ (degrees of freedom)             |
|    this differs from previous releases, where T was 2*E(kin)/(3 NIONS). |
|    The new definition is more consistent                               |
|                                                                            |
-----------------------------------------------------------------------------

-----------------------------------------------------------------------------
|                                                                            |
|          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 !!!          |
|                                                                            |
|    The distance between some ions is very small                         |
|    please check the nearest neigbor list in the OUTCAR file             |
|       I HOPE YOU KNOW, WHAT YOU ARE  DOING                            |
|                                                                            |
-----------------------------------------------------------------------------

POSCAR, INCAR and KPOINTS ok, starting setup
WARNING: wrap around errors must be expected
FFT: planning ...          1
reading WAVECAR
LAPACK: Routine ZPOTRF failed!          2
LAPACK: Routine ZPOTRF failed!          2
LAPACK: Routine ZPOTRF failed!          2
LAPACK: Routine ZPOTRF failed!          2
LAPACK: Routine ZPOTRF failed!          2
LAPACK: Routine ZPOTRF failed!          2
LAPACK: Routine ZPOTRF failed!          2
LAPACK: Routine ZPOTRF failed!          2
   請問這種情況是我的vasp沒編譯好，還是INCAR的問題？謝謝大家，我的INCAR如下：
   NWRITE =    1
PREC =    L
ISTART =    0
ICHARG =    2
# ISPIN  =    2
# MAGMOM =  102*1 1*3 2*1
#  LNONCOLLINEAR = F
#  LSORBIT =    F
#  INIWAV =    1
#  LASPH  =    F
#  METAGGA=    F

# Electronic Relaxation 1

ENCUT  =  400.0 eV
#  ENINI  =  400.0
NELM =    60
NELMIN=  4
NELMDL= -4
EDIFF  = 1E-5
LREAL  =    On
#  LCOMPAT=    F
#  LREAL_COMPAT= F
#  GGA_COMPAT  = T
#  VOSKOWN=    1
#  ROPT = 0.00000 0.00000 0.0
#  NBANDS =  128
#  IALGO  =    38
# GGA=91
ALGO = F
#  LDIAG  =    T

# PAW control tag, page 80
  LMAXPAW = 0
  LMAXMIX = 4
# ENAUG  =  605.4 eV
  ADDGRID = T

# Ionic relaxation
EDIFFG = -0.04
NSW = 400
  NBLOCK =    1
  KBLOCK = 100
  IBRION =    0
  ISIF =    0
  POTIM  = 2
  TEBEG  = 297;
  TEEND  = 500;
  SMASS  =  -1
#  PSTRESS= 0.0

#PC function
APACO=10
#NPACO=200




# POMASS =  52.00 65.39 16.00
# ZVAL = 6.00 12.00  6.00
# RWIGS  = 1.00  1.00  1.00
# NUPDOWN=    0 #for AFM?

DOS related values:

# EMIN = -10.00
# EMAX =  15
#  NEDOS = 800
ISMEAR = 0
SIGMA  = 0.2


#  Mixing tag  page 73 on manual
#  IMIX =    4
# AMIX    = 0.20
BMIX    = 2                #
# AMIX_MAG = 0.8
# BMIX_MAG = 0.0001
MAXMIX=    50
#  AMIN    = 0.10
#  WC =       1000
#  INIMIX=    1
#  MIXPRE=    1


# Write flags
LWAVE  =    F
LCHARG =    F
# LORBIT =    12
#  LVTOT  =    F
#  LELF =    F

# Monopole, Dipole and Quadruole  corr
# NELECT =    636.0000
#  IDIPOL =    4
#  LDIPOL =    F

# L(S)DA+U tag
# LDAU = T
# LDAUTYPE = 2
# LDAUL = 2 -1 -1
# LDAUU = 2.0 0.0 0
# LDAUJ = 0.5 0.0 0
# LDAUPRINT = 2
#LDAUL LDAUU LDAUJ must be set for eve
# Parallel parameter
LPLANE = .TRUE
  NPAR = 4
  ISYM = 0

回復(fù)此樓

» 收錄本帖的淘帖專輯推薦

學(xué)習(xí)

» 猜你喜歡

296求調(diào)劑已經(jīng)有5人回復(fù)
0703化學(xué)調(diào)劑已經(jīng)有13人回復(fù)
085600材料與化工求調(diào)劑已經(jīng)有15人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有4人回復(fù)
生物學(xué)調(diào)劑招人�。�！已經(jīng)有4人回復(fù)
0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章已經(jīng)有17人回復(fù)
一志愿武漢理工材料工程專碩調(diào)劑已經(jīng)有5人回復(fù)
0856調(diào)劑，是學(xué)校就去已經(jīng)有7人回復(fù)
申博26年已經(jīng)有3人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有10人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

【求助】請教vasp計算恒溫第一原理分子動力學(xué) 已經(jīng)有9人回復(fù)
【求助】關(guān)于用分子動力學(xué)來做光譜的問題已經(jīng)有14人回復(fù)
【求助】如何利用分子動力學(xué)加速DFT計算的收斂已經(jīng)有6人回復(fù)
【求助】請教分子動力學(xué)模擬已經(jīng)有8人回復(fù)
【求助】第一原理分子動力學(xué)都能計算什么性質(zhì)呢？已經(jīng)有14人回復(fù)
【求助】分子動力學(xué)的溶劑化分析已經(jīng)有3人回復(fù)

1樓 2011-05-24 10:34:00

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jsheng5059

銀蟲 (小有名氣)

1ST強帖: 1
應(yīng)助: 0 (幼兒園)
金幣: 225.4
散金: 110
紅花: 1
帖子: 184
在線: 77.1小時
蟲號: 779456
注冊: 2009-05-25
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

【答案】應(yīng)助回帖

★
zhang668(金幣+1): 多謝指導(dǎo) 2011-05-24 15:44:24
weiwei128(金幣+10): 謝謝已糾正 2011-05-24 17:14:39

你這里有兩個Warning，第一個可以忽略。但第二個提示原子位置太近了。你需要檢查下POSCAR。計算錯誤可能是原子位置太近造成的。分子動力學(xué)一般是IALGO=48.模擬退火時間步長可以取大一些。

贊一下(1人)

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

skplin

Wair

2樓2011-05-24 12:19:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

enola

捐助貴賓 (正式寫手)

1ST強帖: 2
應(yīng)助: 13 (小學(xué)生)
金幣: 1485.9
散金: 1549
紅花: 14
帖子: 847
在線: 102.5小時
蟲號: 947564
注冊: 2010-01-24
專業(yè): 半導(dǎo)體光電子器件

引用回帖:

Originally posted by weiwei128 at 2011-05-24 10:34:00:
版友們好，我嘗試用MD實現(xiàn)模擬退火，計算體系結(jié)構(gòu)變化。遇到了些問題，提示如下：
-----------------------------------------------------------------------------
| ...

最好把POSCAR用VESTA看一下，有沒有問題。。。有時很明顯的結(jié)構(gòu)不對的

贊一下

回復(fù)此樓

timeflies..

3樓2011-05-24 19:04:56

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

skplin

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 94.9
散金: 66
帖子: 112
在線: 66.6小時
蟲號: 4562736
注冊: 2016-04-01
性別: GG
專業(yè): 金屬材料的力學(xué)行為

送紅花一朵

引用回帖:

2樓: Originally posted by jsheng5059 at 2011-05-24 12:19:08
你這里有兩個Warning，第一個可以忽略。但第二個提示原子位置太近了。你需要檢查下POSCAR。計算錯誤可能是原子位置太近造成的。分子動力學(xué)一般是IALGO=48.模擬退火時間步長可以取大一些。

第一個為什么不用管�。慷夹枰嘶饐�？退火關(guān)鍵的參數(shù)是什么？

贊一下

回復(fù)此樓

4樓2016-04-22 16:09:59

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考博] 申博26年 +3	八6八68 2026-03-19	3/150	2026-03-19 19:43 by nxgogo
[考研] 317求調(diào)劑 +3	申子申申 2026-03-19	6/300	2026-03-19 14:16 by 申子申申
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +4	生物工程調(diào)劑 2026-03-16	12/600	2026-03-19 11:10 by 生物工程調(diào)劑
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +10	秋有木北 2026-03-14	10/500	2026-03-19 05:52 by anny19840123
[考研] 330求調(diào)劑 +3	小材化本科 2026-03-18	3/150	2026-03-18 21:55 by 無懈可擊111
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +8	Ncdx123456 2026-03-13	9/450	2026-03-18 14:40 by haxia
[考研] 0854可跨調(diào)劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287 +8	小李0854 2026-03-16	8/400	2026-03-18 14:35 by 搏擊518
[考研] 070300化學(xué)319求調(diào)劑 +6	錦鯉0909 2026-03-17	6/300	2026-03-18 13:22 by Iveryant
[考研] 302求調(diào)劑 +10	呼呼呼。。。。 2026-03-17	10/500	2026-03-18 12:45 by Linda Hu
[考研] 312求調(diào)劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3	Eember. 2026-03-17	9/450	2026-03-18 10:28 by Eember.
[考研] 0703化學(xué)336分求調(diào)劑 +6	zbzihdhd 2026-03-15	7/350	2026-03-18 09:53 by zhukairuo
[考研] 301求調(diào)劑 +9	yy要上岸呀 2026-03-17	9/450	2026-03-18 08:58 by 無際的草原
[考研] 326求調(diào)劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[考研] 材料與化工專碩調(diào)劑 +5	heming3743 2026-03-16	5/250	2026-03-17 14:03 by 勇敢太監(jiān)王公公
[考博] 26申博 +4	八6八68 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨
[考研] 304求調(diào)劑 +5	素年祭語 2026-03-15	5/250	2026-03-16 17:00 by 我的船我的海
[考研] 304求調(diào)劑 +3	曼殊2266 2026-03-14	3/150	2026-03-16 16:39 by houyaoxu
[考研] 一志愿哈工大材料324分求調(diào)劑 +5	閆旭東 2026-03-14	5/250	2026-03-14 14:53 by 木瓜膏
[考研] 0856材料與化工301求調(diào)劑 +5	奕束光 2026-03-13	5/250	2026-03-13 22:00 by 星空星月