24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 3981  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料考研求调剂总分280 +11 mkjlz1 2026-04-06 11/550 2026-04-06 21:17 by Equinoxhua
[考研] 085404 293求调剂 +3 勇远库爱314 2026-04-06 3/150 2026-04-06 17:34 by 土木硕士招生
[考研] 复试调剂 +5 asdasdassda 2026-04-05 5/250 2026-04-06 09:32 by dongzh2009
[考研] 一志愿武汉理工大学-085601材料工程(专硕)-总分353求调剂 +3 2626262626li 2026-04-02 3/150 2026-04-06 09:08 by 无际的草原
[考研] 085600,320分求调剂 +16 大馋小子 2026-04-04 17/850 2026-04-06 07:58 by MOF_Catal
[考研] 生物与医药调剂 +4 十七sa 2026-04-05 4/200 2026-04-05 20:05 by lys0704
[考研] 调剂 +3 李广火 2026-04-05 3/150 2026-04-05 18:57 by 蓝云思雨
[考研] 生物工程求调剂 +6 喜欢还是不甘心 2026-04-05 6/300 2026-04-05 10:28 by 唐沐儿
[考研] 285求调剂 +11 哦呦呼o 2026-04-04 11/550 2026-04-05 08:15 by 544594351
[考研] 294求调剂 +6 Grey_Ey 2026-04-02 9/450 2026-04-04 22:07 by hemengdong
[考研] 调剂 +4 是可乐不是可乐 2026-04-04 4/200 2026-04-04 19:41 by 唐沐儿
[考研] 一志愿东北大学085901土木专硕345求调剂 +3 zxt11111 2026-04-04 3/150 2026-04-04 14:21 by 土木硕士招生
[考研] 本9一志愿2 0854低分专硕286求调剂 +9 芒种111 2026-04-04 9/450 2026-04-04 11:01 by tangruihua
[考研] 化工求调剂 +11 荔香芝士椰奶 2026-04-03 11/550 2026-04-03 22:06 by 啵啵啵0119
[考研] 求调剂 +4 压力??大 2026-04-03 4/200 2026-04-03 21:36 by 啵啵啵0119
[考研] 338求调剂 +4 zzz,,r 2026-04-03 4/200 2026-04-03 16:39 by lijunpoly
[考研] 英一数一408,总分284,二战真诚求调剂 +13 12.27 2026-03-30 15/750 2026-04-03 14:41 by 氮气气气
[考研] 366求调剂 +7 sbdnd 2026-04-03 7/350 2026-04-03 12:40 by cymywx
[考研] 325分化学调剂 +5 15771691647 2026-04-02 5/250 2026-04-03 09:58 by ChemPharm
[考研] 318求调剂,计算材料方向 +10 吸喵有害笙命 2026-04-01 11/550 2026-04-02 16:29 by oooqiao
信息提示
请填处理意见