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

» 猜你喜歡

288求調(diào)劑一志愿哈工大材料與化工已經(jīng)有9人回復(fù)
化學(xué)工程085602 305分求調(diào)劑已經(jīng)有31人回復(fù)
070300化學(xué)專業(yè)279調(diào)劑已經(jīng)有3人回復(fù)
304求調(diào)劑已經(jīng)有6人回復(fù)
085404 22408 315分已經(jīng)有5人回復(fù)
一志愿a區(qū)211，085601-307分求調(diào)劑已經(jīng)有10人回復(fù)
求調(diào)劑已經(jīng)有7人回復(fù)
290分調(diào)劑求助已經(jīng)有7人回復(fù)
一志愿華東理工大學(xué)，080500學(xué)碩，317分，求調(diào)劑已經(jīng)有3人回復(fù)
考研調(diào)劑求助已經(jīng)有4人回復(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的回帖

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的回帖

查看全部 4 個回答

jsheng5059

銀蟲 (小有名氣)

1ST強(qiáng)帖: 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強(qiáng)帖: 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的回帖

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

[考研] 311求調(diào)劑一志愿合肥工業(yè)大學(xué) +9	秋二十二 2026-03-30	9/450	2026-03-31 13:15 by Dyhoer
[考研] 本科211生物醫(yī)學(xué)工程085409求調(diào)劑339分 +5	里子木yy 2026-03-29	5/250	2026-03-31 10:44 by WWM123
[考研] 一志愿南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +5	@taotao 2026-03-30	5/250	2026-03-31 09:41 by zhshch
[考研] 一志愿西交288化工專碩求調(diào)劑 +4	好運好運接接 2026-03-24	4/200	2026-03-31 09:18 by 西京學(xué)院土木工?/a>
[考研] 本科211總分289，08工學(xué)真心求調(diào)劑 +3	utopiaE 2026-03-30	3/150	2026-03-30 23:42 by ms629
[考研] 085701環(huán)境工程求調(diào)劑 +11	多久上課 2026-03-27	12/600	2026-03-30 21:21 by 研究僧導(dǎo)導(dǎo)
[考研] 317求調(diào)劑 +8	十閑wx 2026-03-24	8/400	2026-03-30 21:00 by peike
[考研] 調(diào)劑 +4	GK72 2026-03-30	4/200	2026-03-30 20:32 by dick_runner
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +9	hanamiko 2026-03-27	9/450	2026-03-30 20:10 by 無際的草原
[考研] 各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩 +10	Koxui 2026-03-28	10/500	2026-03-30 19:33 by 源_2020
[考研] 材料與化工304求B區(qū)調(diào)劑 +4	邱gl 2026-03-26	7/350	2026-03-30 08:39 by 探123
[考研] 一志愿北京工業(yè)大學(xué)，324分求調(diào)劑 +6	零八# 2026-03-28	6/300	2026-03-29 21:20 by nanaliuyun
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +5	hanamiko 2026-03-29	5/250	2026-03-29 16:46 by 學(xué)員8dgXkO
[考研] 壓國家一區(qū)線，求導(dǎo)師收留，有恩必謝！ +7	迷人的哈哈 2026-03-28	7/350	2026-03-28 16:47 by 催化大白
[考研] 本科新能源科學(xué)與工程，一志愿華理能動285求調(diào)劑 +3	AZMK 2026-03-27	5/250	2026-03-28 16:19 by xxxsssccc
[考研] 304求調(diào)劑 +6	曼殊2266 2026-03-27	6/300	2026-03-28 14:10 by 唐沐兒
[考研] 299求調(diào)劑 +7	嗯嗯嗯嗯2 2026-03-27	7/350	2026-03-28 13:09 by 唐沐兒
[考研] 0703化學(xué)求調(diào)劑，各位老師看看我�。�！ +5	祁祺祺 2026-03-25	5/250	2026-03-27 21:44 by 東方豬豬
[考研] 266求調(diào)劑 +11	陽陽哇塞 2026-03-27	12/600	2026-03-27 17:56 by yu221
[考研] 08開頭275求調(diào)劑 +4	拉誰不重要 2026-03-26	4/200	2026-03-27 14:12 by Delta2012