24小时热门版块排行榜    

查看: 689  |  回复: 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的回帖

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的回帖
查看全部 5 个回答

nansh13

金虫 (正式写手)

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

另外    你是怎么分析输出文件的  比如  history   rdf   vaf  zdf的
苦心人天不负卧薪尝胆三千越甲可吞吴
2楼2014-09-22 16:27:55
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 北科281学硕材料求调剂 +7 tcxiaoxx 2026-03-20 7/350 2026-03-23 09:28 by haoshis
[考研] 276求调剂 +3 YNRYG 2026-03-21 4/200 2026-03-23 08:31 by 醉在风里
[考研] 石河子大学(211、双一流)硕博研究生长期招生公告 +3 李子目 2026-03-22 3/150 2026-03-22 21:01 by 怎么释怀
[考研] 324求调剂 +6 lucky呀呀呀鸭 2026-03-20 6/300 2026-03-22 16:01 by ColorlessPI
[考研] 0856材料专硕353求调剂 +4 NIFFFfff 2026-03-20 4/200 2026-03-22 09:49 by 2026paper
[考研] 材料学硕301分求调剂 +7 Liyouyumairs 2026-03-21 7/350 2026-03-21 22:31 by peike
[考研] 297求调剂 +3 喜欢还是不甘心 2026-03-20 3/150 2026-03-21 18:33 by 学员8dgXkO
[考研] 材料学硕333求调剂 +3 北道巷 2026-03-18 3/150 2026-03-21 18:17 by 学员8dgXkO
[考研] 302求调剂 +12 呼呼呼。。。。 2026-03-17 12/600 2026-03-21 17:29 by ColorlessPI
[考研] 306求0703调剂一志愿华中师范 +5 纸鱼ly 2026-03-21 5/250 2026-03-21 17:11 by 学员8dgXkO
[考研] 求调剂 +6 Mqqqqqq 2026-03-19 6/300 2026-03-21 08:04 by JourneyLucky
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂! +3 不想理你 2026-03-16 3/150 2026-03-21 03:59 by JourneyLucky
[考研] 332求调剂 +4 ydfyh 2026-03-17 4/200 2026-03-21 02:20 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +3 晨昏线与星海 2026-03-18 3/150 2026-03-21 00:46 by JourneyLucky
[考研] 294求调剂材料与化工专硕 +15 陌の森林 2026-03-18 15/750 2026-03-20 23:28 by JourneyLucky
[考研] A区线材料学调剂 +5 周周无极 2026-03-20 5/250 2026-03-20 21:33 by laoshidan
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
[考研] 0703化学调剂 +5 pupcoco 2026-03-17 8/400 2026-03-19 13:58 by houyaoxu
[考研] 0854可跨调剂,一作一项核心论文五项专利,省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 277调剂 +5 自由煎饼果子 2026-03-16 6/300 2026-03-17 19:26 by 李leezz
信息提示
请填处理意见