应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 怎么分析分子动力学模块cp.x的运算结果？
52/1返回列表
查看: 3972  |  回复: 20
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

huazhorg

铁杆木虫 (正式写手)

[求助] 怎么分析分子动力学模块cp.x的运算结果？

如题。espresso 里example 18 以及随后的几个例子就是运用cp.x进行动力学计算。比如18里，对sio2的计算输入文件sio2.cp.start.in：
&control
    calculation='cp',
    restart_mode='from_scratch',
    nstep=20, iprint=20, isave=20,
    dt=5.0,
    ndr=90, ndw=91,
    pseudo_dir='/home/huazhong/espresso/pseudo/',
    outdir='/home/huazhong/tmp/',
/
&system
    ibrav=8, celldm(1)=9.28990, celldm(2)=1.73206, celldm(3)=1.09955,
    nat=18, ntyp=2, nbnd=48, nspin=1,
    ecutwfc=20.0, ecutrho=150.0,
    nr1b=16, nr2b=16, nr3b=16,
    qcutz=150., q2sigma=2.0, ecfixed=16.0,
/
&electrons
    electron_dynamics='damp', electron_damping=0.2,
    startingwfc='random', ampre=0.01,
    emass=700., emass_cutoff=3.,
/
&ions
    ion_dynamics='none',
    ion_radius(1)=1.0, ion_radius(2)=1.0,
/
ATOMIC_SPECIES
   O  16.00 O.pz-rrkjus.UPF
  Si  28.00 Si.pz-vbc.UPF
ATOMIC_POSITIONS
O  3.18829368  14.83237039   1.22882961
O  7.83231469   6.78704039   1.22882961
O  2.07443467   5.99537992   4.73758250
O  6.72031366  14.04231898   4.73758250
O  3.96307134  11.26989826   7.87860582
O  8.60802134   3.22295920   7.87860582
O  3.96307134   4.81915267   9.14625133
O  8.60802134  12.86448267   9.14625133
O  3.18736469   1.25668055   5.58029607
O  7.83324368   9.30201055   5.58029607
O  2.07536366  10.09206195   2.07358613
O  6.71938467   2.04673195   2.07358613
Si  0.28891589   8.04533000   3.40456284
Si  4.93386589   0.00000000   3.40456284
Si  2.13389003  12.27717358  -0.04188031
Si  6.77884003   4.23184358  -0.04188031
Si  2.13389003   3.81348642   6.85202747
Si  6.77884003  11.85881642   6.8520274

这个输入文件可以用xcrysden查看结构。其对应的输出文件为sio2.cp.start.out，可是这个输出文件用xcrysden不能查看结构，运行完后几何结构变成什么样了，也不知道。那么怎么分析这个输出文件？请各位仁兄、大侠前来指点。(输出文件sio2.cp.start.out见附件)[ Last edited by huazhorg on 2012-11-27 at 22:06 ]
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : sio2.cp.start.out
    • 2012-11-27 22:06:16, 15.82 K

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼 2012-11-27 22:00:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huazhorg

铁杆木虫 (正式写手)
引用回帖:
12楼: Originally posted by valenhou001 at 2012-11-30 08:37:47
正如所看到的信息一样，你应该指出你所用的cp或quantum espresso版本。
最新版本的cpv，不需要去修改它。

在你的输入文件中试试增加 “saverho = .TRUE." 再看。

谢谢侯老师！saverho = .TRUE.也不起作用，后来把老版本（4.2.1）全部删除，重新安装最新版5.0.2后，才没有问题，cppp.x可以运行了。运行后会产生一个xsf文件，这个文件可以用xcrysden查看结构，但这个文件并不是一个动画，而是一个静态结构，我现在想知道，这个结构是优化后的结构呢，还是中间某一步的结构呢？要想查看动力学优化中的结构变化过程，该怎么算呢？
赞一下 回复此楼
13楼2012-11-30 16:26:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 21 个回答

valenhou001

至尊木虫 (职业作家)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
huazhorg: 金币+2, ★有帮助, 谢谢指点 2012-11-28 16:38:08
huazhorg: 回帖置顶 2012-11-29 17:13:58
建议你先看一些有关cpmd或cp.x的网上教程。
赞一下 回复此楼
2楼2012-11-28 16:07:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huazhorg

铁杆木虫 (正式写手)
引用回帖:
2楼: Originally posted by valenhou001 at 2012-11-28 16:07:21
建议你先看一些有关cpmd或cp.x的网上教程。

谢谢侯老师，很高兴得到您的指点。看了一些教程，似乎略有点明白。它的优化坐标不是放在out文件里，而是放在outdir里一个‘prefix’.pos文件里吧，比如我算了一个水分子的优化，它的坐标是不是这个文件：
    150  0.01456168
  0.500990E+01  0.500990E+01  0.500000E+01
  0.675162E+01  0.468700E+01  0.509543E+01
  0.481703E+01  0.681514E+01  0.494747E+01
    160  0.01577113
  0.500990E+01  0.500990E+01  0.500000E+01
  0.653486E+01  0.438670E+01  0.541771E+01
  0.495464E+01  0.677315E+01  0.477299E+01
    170  0.01698057
  0.500990E+01  0.500990E+01  0.500000E+01
  0.650987E+01  0.424898E+01  0.556919E+01
  0.502422E+01  0.677687E+01  0.469369E+01
    180  0.01819001
  0.500990E+01  0.500990E+01  0.500000E+01
  0.652258E+01  0.419161E+01  0.564499E+01
  0.507221E+01  0.679087E+01  0.465831E+01
    190  0.01939945
  0.500990E+01  0.500990E+01  0.500000E+01
  0.653114E+01  0.418351E+01  0.567759E+01
  0.511155E+01  0.680055E+01  0.464628E+01
    200  0.02060889
  0.500990E+01  0.500990E+01  0.500000E+01
  0.653185E+01  0.420141E+01  0.568724E+01
  0.514636E+01  0.680329E+01  0.464624E+01
    210  0.02181834
  0.500990E+01  0.500990E+01  0.500000E+01
  0.653089E+01  0.422882E+01  0.568727E+01
  0.517752E+01  0.680156E+01  0.465167E+01
    220  0.02302778
  0.500990E+01  0.500990E+01  0.500000E+01
  0.653328E+01  0.425626E+01  0.568513E+01
  0.520506E+01  0.679849E+01  0.465888E+01
    230  0.02423722
  0.500990E+01  0.500990E+01  0.500000E+01
  0.654017E+01  0.427980E+01  0.568387E+01
  0.522899E+01  0.679605E+01  0.466601E+01
    240  0.02544666
  0.500990E+01  0.500990E+01  0.500000E+01
  0.654995E+01  0.429893E+01  0.568389E+01
  0.524946E+01  0.679487E+01  0.467228E+01
    250  0.02665611
  0.500990E+01  0.500990E+01  0.500000E+01
  0.656023E+01  0.431446E+01  0.568456E+01
  0.526672E+01  0.679468E+01  0.467751E+01
    260  0.02786555
  0.500990E+01  0.500990E+01  0.500000E+01
  0.656927E+01  0.432735E+01  0.568528E+01
  0.528110E+01  0.679498E+01  0.468180E+01
    270  0.02907499
  0.500990E+01  0.500990E+01  0.500000E+01
  0.657642E+01  0.433824E+01  0.568574E+01
  0.529294E+01  0.679531E+01  0.468528E+01
    280  0.03028443
  0.500990E+01  0.500990E+01  0.500000E+01
  0.658176E+01  0.434747E+01  0.568594E+01
  0.530259E+01  0.679548E+01  0.468808E+01
    290  0.03149387
  0.500990E+01  0.500990E+01  0.500000E+01
  0.658574E+01  0.435520E+01  0.568601E+01
  0.531038E+01  0.679547E+01  0.469031E+01
    300  0.03270332
  0.500990E+01  0.500990E+01  0.500000E+01
  0.658877E+01  0.436158E+01  0.568607E+01
  0.531662E+01  0.679534E+01  0.469205E+01
    310  0.03391276
  0.500990E+01  0.500990E+01  0.500000E+01
  0.659119E+01  0.436678E+01  0.568619E+01
  0.532158E+01  0.679517E+01  0.469337E+01
    320  0.03512220
  0.500990E+01  0.500990E+01  0.500000E+01
  0.659315E+01  0.437096E+01  0.568638E+01
  0.532550E+01  0.679501E+01  0.469434E+01
    330  0.03633164
  0.500990E+01  0.500990E+01  0.500000E+01
  0.659475E+01  0.437434E+01  0.568664E+01
  0.532858E+01  0.679489E+01  0.469502E+01
    340  0.03754108
  0.500990E+01  0.500990E+01  0.500000E+01
  0.659604E+01  0.437706E+01  0.568694E+01
  0.533098E+01  0.679480E+01  0.469545E+01

这个坐标是不是它的优化过程？另外怎么可视化这个过程呢，是不是要写个脚本把它转化成xyz的格式，然后再用xmakemol等这样的小工具来看优化过程？感觉这样好麻烦呢。。
赞一下 回复此楼
3楼2012-11-28 17:04:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
huazhorg: 金币+3, ★有帮助, 谢谢！ 2012-11-29 11:27:58
liliangfang: 金币+2, 应助指数+1, 专家考核, 谢谢交流 2012-12-01 20:21:36
该程序提供了一个cppp.x的模块，可进行数据处理。有关cppp.x模块的使用，请看CPV/Doc里的手册。
赞一下 回复此楼
4楼2012-11-29 08:20:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 21 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 19:43 by lXgU0RvNTDC9
[公派出国] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 7/350 2026-08-23 19:10 by lXgU0RvNTDC9
[论文投稿] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 8/400 2026-08-23 18:58 by lXgU0RvNTDC9
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-21 5/250 2026-08-23 14:33 by KM5EcsNQRBPn
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 什么时候开奖？ +9 CrisMessi 2026-08-18 10/500 2026-08-23 12:03 by 丶昵称占用
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 9/450 2026-08-23 10:44 by IXZuIJ2Q7OVy
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 5/250 2026-08-23 04:07 by OEbVnUOu01ol
[论文投稿] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 04:04 by OEbVnUOu01ol
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] 范进中举一文的中心思想 +3 炎黄贵胄 2026-08-22 4/200 2026-08-22 20:02 by mycsru
[基金申请] 今天基金会出结果吗？20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 今天系统多次维护，明天很可能放榜！ +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
关闭
请填处理意见
关闭