24小时热门版块排行榜    

查看: 794  |  回复: 4

ddwy

木虫 (小有名气)

[求助] HISTORY

大家好:
我用DL_POLY做模拟过程中遇见一个问题:CONTROL文件中明明写的是从第2000001步开始写HISTORY,每200步记录一次,写到第5000000步结束,正常情况下,HISTORY文件中一共要记录3000000步,也就是15000个snapshots,但是在HISTORY中却只记录了10548个,是怎么回事?下面是相应的CONTROL、HISTORY文件的主要部分:
temperature 298.15
# pressure 0.001  
ensemble nvt hoover  0.5

integrat velocity
steps 5000000
equilibration 2000000
mult  1
scale  5
stats 20
rdf            10
print rdf

traj 2000001 200 0

densvar 10

timestep 0.001
cutoff        17.5

ewald precision 1.0E-6

minimise energy 50 5.0E-03

shake 1.0E-4
mxshak  1000
quaternion tolerance  1.0E-05

job time       99999999999999.0
close time     200000.0

finish



HISTORY文件:

This is the title line. Change this at dl_field.control. Reads only 80 c
         0         1      5075                10548            107104394
timestep   2000001      5075 0 1            0.001000         2000.001000
       36.5411224570        0.0000000000        0.0000000000
        0.0000000000       36.5411224570        0.0000000000
        0.0000000000        0.0000000000       36.5411224570
n3               1   14.006700   -0.300000    0.001321
    -14.39204560        -7.129439892         6.655407644
hn               2    1.007970    0.330000    0.025434
    -13.71943858        -6.342743032         6.833563250
hn               3    1.007970    0.330000    0.012115
    -13.90251563        -8.029155915         6.825459467
hn               4    1.007970    0.330000    0.013617
    -15.16152579        -7.073129145         7.361173963
c1               5   12.011000    0.210000    0.006457
    -14.92887974        -6.947868975         5.279145164
h3               6    1.007970    0.100000    0.015502
    -15.51765202        -7.813452274         5.001153200
c3               7   12.011000   -0.270000    0.007139
    -15.81760489        -5.681052690         5.372765110
h1               8    1.007970    0.090000    0.031123
    -16.52371735        -5.607495283         4.492964758
h1               9    1.007970    0.090000    0.045888
    -16.58945209        -5.725914663         6.226788007
h1              10    1.007970    0.090000    0.055990
    -15.22617482        -4.803296129         5.430600240
c6              11   12.011000    0.340000    0.008470
    -13.85244286        -6.772854713         4.181272071
oc3             12   15.999400   -0.670000    0.003495
    -12.64690293        -6.612563759         4.513710555
oc3             13   15.999400   -0.670000    0.004525
    -14.19533041        -6.789528726         2.962201713
n3              14   14.006700   -0.300000    0.007100
     9.083714187         13.12434088         3.230931917
hn              15    1.007970    0.330000    0.015408
     9.799614544         13.85093581         3.004330322
hn              16    1.007970    0.330000    0.013673
     9.209553990         12.78405456         4.214747430
hn              17    1.007970    0.330000    0.011592
     8.123385047         13.54147225         3.109655461
c1              18   12.011000    0.210000    0.004365
     9.218579513         11.99469629         2.284462493
h3              19    1.007970    0.100000    0.013782
     10.00071147         11.32116123         2.636253965
c3              20   12.011000   -0.270000    0.004674
     7.854011719         11.28476128         2.251172830
h1              21    1.007970    0.090000    0.060129
     7.842466044         10.41942858         1.444509710
h1              22    1.007970    0.090000    0.068361
     7.608257112         10.87184244         3.169401447
h1              23    1.007970    0.090000    0.000650
     7.062671444         12.00666112         1.941293063
c6              24   12.011000    0.340000    0.003390
     9.547663706         12.45036844        0.8421650024
oc3             25   15.999400   -0.670000    0.003469




谢谢各位
回复此楼
我喜欢收集微笑没有原因……
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nansh13

金虫 (正式写手)

猜想  你可以把jobtime放大    扩大100倍

另外    你是怎么分析输出文件的  比如  history   rdf   vaf  zdf的
苦心人天不负卧薪尝胆三千越甲可吞吴
2楼2014-09-22 16:27:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lanlinger13

金虫 (小有名气)

引用回帖:
2楼: Originally posted by nansh13 at 2014-09-22 16:27:55
猜想  你可以把jobtime放大    扩大100倍

另外    你是怎么分析输出文件的  比如  history   rdf   vaf  zdf的

dlpoly的GUI貌似可以分析这些的
3楼2014-11-27 10:29:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ddwy

木虫 (小有名气)

引用回帖:
2楼: Originally posted by nansh13 at 2014-09-22 16:27:55
猜想  你可以把jobtime放大    扩大100倍

另外    你是怎么分析输出文件的  比如  history   rdf   vaf  zdf的

HISTORY文件我们是提取溶液中溶质的结构,就是去掉水分子,还可以提取不同原子数目的团簇,以及存在时间最长的团簇,这些自己可以写一些小程序
我喜欢收集微笑没有原因……
4楼2014-11-28 09:20:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ddwy

木虫 (小有名气)

引用回帖:
3楼: Originally posted by lanlinger13 at 2014-11-27 10:29:05
dlpoly的GUI貌似可以分析这些的...

嗯,我试着分析一下,谢谢
我喜欢收集微笑没有原因……
5楼2014-11-28 09:21:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ddwy 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] fileCode有新解读? +7 Tide man 2026-08-08 13/650 2026-08-08 09:47 by zhengxianzb
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 08:59 by 6vVgjDL4CnGu
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 4/200 2026-08-08 08:27 by 6vVgjDL4CnGu
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 KXLV3nuBVBY7 2026-08-07 4/200 2026-08-08 08:21 by 6vVgjDL4CnGu
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 3/150 2026-08-08 03:47 by 6vVgjDL4CnGu
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
信息提示
请填处理意见