24小时热门版块排行榜    

查看: 4296  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +9 majunge000 2026-08-11 11/550 2026-08-15 23:20 by botar
[基金申请] 时间戳又变了8-15 +12 archvillain 2026-08-15 23/1150 2026-08-15 20:47 by foyo404
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] filecode +8 cratir 2026-08-14 12/600 2026-08-15 18:08 by zyfgau
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +7 Ldrop2023 2026-08-13 7/350 2026-08-15 15:57 by mzhh2000
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 22/1100 2026-08-14 23:45 by Noways
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
请填处理意见