24小时热门版块排行榜    

查看: 804  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 梦想 +5 myrtle 2026-08-26 5/250 2026-08-31 14:40 by hahaboy
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 国社科又开始会评了,不知道这次命运如何 +6 雨打竹帘 2026-08-30 10/500 2026-08-31 09:29 by huixian257
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 能否申诉? +6 echo8914667 2026-08-30 7/350 2026-08-31 08:41 by kingkocxr
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见