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í)

» 猜你喜歡

一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑已經(jīng)有10人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有8人回復(fù)
一志愿中海洋材料工程專碩330分求調(diào)劑已經(jīng)有6人回復(fù)
一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑已經(jīng)有3人回復(fù)
332求調(diào)劑已經(jīng)有3人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有7人回復(fù)
281求調(diào)劑（0805）已經(jīng)有17人回復(fù)
267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑已經(jīng)有10人回復(fù)
0817調(diào)劑已經(jīng)有3人回復(fù)
焦慮已經(jīng)有13人回復(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ā)表

[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +4	守候夕陽CF 2026-03-18	4/200	2026-03-18 22:16 by li123456789.
[考研] 330求調(diào)劑 +3	小材化本科 2026-03-18	3/150	2026-03-18 21:55 by 無懈可擊111
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-17	7/350	2026-03-18 20:41 by Wangjingyue
[考研] 326求調(diào)劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4	張vvvv 2026-03-15	6/300	2026-03-17 17:15 by ruiyingmiao
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 26考研求調(diào)劑 +6	丶宏Sir 2026-03-13	6/300	2026-03-17 16:13 by 醉在風(fēng)里
[考研] 梁成偉老師課題組歡迎你的加入 +8	一鴨鴨喲 2026-03-14	10/500	2026-03-17 15:07 by 一鴨鴨喲
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3	木子momo青爭 2026-03-15	3/150	2026-03-17 07:52 by laoshidan
[考研] 一志愿211 0703方向310分求調(diào)劑 +3	努力奮斗112 2026-03-15	3/150	2026-03-16 16:44 by houyaoxu
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 285求調(diào)劑 +6	ytter 2026-03-12	6/300	2026-03-16 15:05 by njzyff
[考研] 22408總分284求調(diào)劑 +3	InAspic 2026-03-13	3/150	2026-03-15 11:10 by zhq0425
[考研] 材料080500調(diào)劑求收留 +3	一顆meteor 2026-03-13	3/150	2026-03-14 10:54 by peike
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機器學(xué)習(xí) 一志愿哈工大 +4	林yaxin 2026-03-12	4/200	2026-03-13 22:04 by 星空星月
[考研] ［0860］321分求調(diào)劑，ab區(qū)皆可 +4	寶貴熱 2026-03-13	4/200	2026-03-13 22:01 by 星空星月
[考研] 考研調(diào)劑 +4	芬達46 2026-03-12	4/200	2026-03-13 16:04 by ruiyingmiao
[考研] 工科材料085601 279求調(diào)劑 +8	困于星晨 2026-03-12	10/500	2026-03-13 15:42 by ms629
[考研] 289求調(diào)劑 +3	李政瑩 2026-03-12	3/150	2026-03-13 11:02 by 求調(diào)劑zz
[考研] 290求調(diào)劑 +3	ADT 2026-03-13	3/150	2026-03-13 10:19 by peike