24小时热门版块排行榜    

Znn3bq.jpeg
查看: 4066  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026山东省优青 +3 luckinging 2026-05-14 4/200 2026-05-19 16:15 by 52wls
[论文投稿] 求助大佬sci投稿哪个好中 +4 江沅188 2026-05-12 5/250 2026-05-19 15:50 by CFD_Helper
[文学芳草园] 献血感触 +4 呀呀好傻 2026-05-19 4/200 2026-05-19 15:09 by seaskyy
[论文投稿] Sci. Bull. 悲剧经验 +4 jyang1999 2026-05-16 4/200 2026-05-19 10:56 by maolC
[考博] 找博士生导师 +6 小代想上岸 2026-05-15 7/350 2026-05-19 10:22 by free_fisher
[基金申请] 同样的基金本子,换个专家直接从C变A! +3 国自然国社科中 2026-05-19 3/150 2026-05-19 08:50 by Equinoxhua
[考博] 2026博士还有哪些学校有名额 +7 小王求读研 2026-05-15 8/400 2026-05-19 08:27 by zhyzzh
[考博] 云南大学材料与能源学院解琳课题组钙钛矿博士招生 +3 光伏爱好者 2026-05-17 4/200 2026-05-19 01:19 by 风起于浮萍
[考博] 博士申请 +5 星…… 2026-05-18 6/300 2026-05-18 23:49 by 糊糊涂涂好
[基金申请] 面上本子正文33页,违规吗?会被低分嘛? +8 1234567wang 2026-05-17 10/500 2026-05-18 18:52 by zzahkj
[基金申请] 今年审到国自然15份,谈谈感受 +16 国自然国社科中 2026-05-17 16/800 2026-05-18 14:58 by gy116024
[硕博家园] 考博自荐 +5 科研狗111 2026-05-13 6/300 2026-05-18 11:22 by 糊糊涂涂好
[基金申请] 青C资助名额大幅增加! +12 西葫芦炒鸡蛋 2026-05-13 16/800 2026-05-18 10:02 by Equinoxhua
[硕博家园] 我在等一个没有答案的答案 +3 Love_MH 2026-05-17 3/150 2026-05-18 02:22 by 竹林孤影
[考博] 光量子物理方向 博士招生 1人(2026.09) +3 sandyworld 2026-05-15 4/200 2026-05-17 14:38 by sandyworld
[有机交流] 求助2,4-二氯-5-嘧啶甲醛的合成方法 20+3 光吃不拉 2026-05-14 6/300 2026-05-16 19:46 by Equinoxhua
[有机交流] 如何实现卤原子转化 +3 BT20230424 2026-05-15 5/250 2026-05-16 16:20 by czyzsu
[文学芳草园] 风把牡丹吹跑了 +5 myrtle 2026-05-12 9/450 2026-05-15 15:27 by myrtle
[考博] 26应届毕业生考博求助 +3 wo一定上岸 2026-05-13 3/150 2026-05-14 21:47 by 明海天涯
[考博] 材料类只有一篇综述能申博么 +4 乐逍遥谷 2026-05-13 4/200 2026-05-14 12:05 by zhyzzh
信息提示
请填处理意见