24小时热门版块排行榜    

查看: 3857  |  回复: 8

xiaoxiao6726

银虫 (初入文坛)

[求助] 【求助】VASP做分子动力学模拟NVT系综下为啥内能竟然是升高的?

很奇怪,在用VASP做分子动力学模拟时,同样的INCAR设置下有些体系内能是降低的,有些体系竟然是升高的,如图所示,请教各位高手这可能是什么原因造成的?非常感谢!
回复此楼
爱生活爱宝宝
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

叶落鹰飞

【答案】应助回帖

★ ★
uuv2010(金币+2): 欢迎参加交流 2011-08-01 12:09:34
xiaoxiao6726(金币+10): 多谢啦,还想请教一个问题,在什么样的误差范围内可以认为内能的变化是可以接受的呢? 2011-08-01 15:02:27
正常,主要是由于计算误差导致的。   
改善的办法是减小MD的步长。
真理是一点点接近的!
2楼2011-08-01 10:13:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Gina88

木虫 (正式写手)

引用回帖:
2楼: Originally posted by gavinliu7390 at 2011-08-01 10:13:03:
正常,主要是由于计算误差导致的。   
改善的办法是减小MD的步长。

请问大侠,分子动力学的时候出现以下警告怎么办啊?

Successfully completed.

Resource usage summary:

    CPU time   :16096231.00 sec.
    Max Memory :     22146 MB
    Max Swap   :     25954 MB

    Max Processes  :        46
    Max Threads    :        79

The output (if any) follows:

running on   32 nodes
distr:  one band on    1 nodes,   32 groups
vasp.5.2.2 15Apr09 complex
POSCAR found :  1 types and   74 ions
LDA part: xc-table for Ceperly-Alder, Vosko type interpolation para-ferro

-----------------------------------------------------------------------------
|                                                                             |
|           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    219 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                                  |
|                                                                             |
-----------------------------------------------------------------------------

POSCAR, INCAR and KPOINTS ok, starting setup
FFT: planning ...(            1 )
reading WAVECAR
prediction of wavefunctions initialized - no I/O
entering main loop
       N       E                     dE             d eps       ncg     rms          rms(c)
RMM:   1     0.258916673329E+04    0.25892E+04   -0.47678E+04   768   0.587E+02
RMM:   2     0.830026842535E+03   -0.17591E+04   -0.20517E+04   768   0.134E+02
RMM:   3     0.239121326438E+03   -0.59091E+03   -0.94025E+03   768   0.957E+01
RMM:   4    -0.155929959519E+03   -0.39505E+03   -0.39046E+03   768   0.775E+01
RMM:   5    -0.280110078607E+03   -0.12418E+03   -0.11962E+03   768   0.478E+01
RMM:   6    -0.329937330070E+03   -0.49827E+02   -0.41162E+02   768   0.299E+01
RMM:   7    -0.348843301236E+03   -0.18906E+02   -0.16211E+02   768   0.185E+01
RMM:   8    -0.357263286004E+03   -0.84200E+01   -0.71033E+01   768   0.125E+01
RMM:   9    -0.364251388673E+03   -0.69881E+01   -0.66987E+01  1946   0.831E+00
3楼2011-08-01 10:44:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

叶落鹰飞

【答案】应助回帖

★ ★
franch(金币+2): 谢谢回帖交流 2011-08-01 21:30:28
引用回帖:
3楼: Originally posted by Gina88 at 2011-08-01 10:44:00:
请问大侠,分子动力学的时候出现以下警告怎么办啊?

Successfully completed.

Resource usage summary:

    CPU time   :16096231.00 sec.
    Max Memory :     22146 MB
    Max Swap   :     2595 ...

可以忽略,不用管它。
真理是一点点接近的!
4楼2011-08-01 14:35:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

enola

捐助贵宾 (正式写手)

引用回帖:
1楼: Originally posted by xiaoxiao6726 at 2011-08-01 09:08:50:
很奇怪,在用VASP做分子动力学模拟时,同样的INCAR设置下有些体系内能是降低的,有些体系竟然是升高的,如图所示,请教各位高手这可能是什么原因造成的?非常感谢!

xiaoxiao6726,你好,请问你在做NVT模拟时,怎么样来设置热浴呢?使用smear>0吗?
谢谢
timeflies..
5楼2011-08-03 16:20:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoxiao6726

银虫 (初入文坛)

引用回帖:
5楼: Originally posted by enola at 2011-08-03 16:20:34:
xiaoxiao6726,你好,请问你在做NVT模拟时,怎么样来设置热浴呢?使用smear>0吗?
谢谢

SMASS=0就可以吧,你可以看OUTCAR文件,它会自动计算并设置成一个合适的值。

» 本帖已获得的红花(最新10朵)

爱生活爱宝宝
6楼2011-08-04 11:23:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

enola

捐助贵宾 (正式写手)

引用回帖:
6楼: Originally posted by xiaoxiao6726 at 2011-08-04 11:23:44:
SMASS=0就可以吧,你可以看OUTCAR文件,它会自动计算并设置成一个合适的值。

好的,谢谢你阿。。。我去试试看
timeflies..
7楼2011-08-05 08:47:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

凌云志123

金虫 (正式写手)

你好,请问你的内能一直升高的问题解决了吗?我也遇到同样的问题
8楼2015-11-27 15:45:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

10133202

铁虫 (小有名气)

送红花一朵
引用回帖:
6楼: Originally posted by xiaoxiao6726 at 2011-08-04 11:23:44
SMASS=0就可以吧,你可以看OUTCAR文件,它会自动计算并设置成一个合适的值。...

xiaoxiao您好,请问在outcar文件里怎么看smass优化的结果呢?您说的意思是这样吗:假定模拟NVT系综,SMASS设定大于等于零的任意一个值,然后OUTCAR中会将SMASS的合理的值给出来吗?那怎么看呢?能否具体讲解下?谢谢亲。
9楼2018-12-08 13:15:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaoxiao6726 的主题更新
信息提示
请填处理意见