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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 怎么分析分子动力学模块cp.x的运算结果？
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 3603  |  回复: 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

铁杆木虫 (正式写手)
我是按照这上面的步骤操作的，怎么也会出错：
http://www.fisica.uniud.it/~giannozz/QE-Tutorial/tutorial_cp.html
（运行例子也是下载这个页面上的。）
赞一下 回复此楼
6楼2012-11-29 11:33:14
已阅   回复此楼   关注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 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600 材料与化工 329分求调剂 +13 Mr. Z 2026-03-25 13/650 2026-03-29 15:24 by 唐沐儿
[考研] 求调剂一志愿武汉理工大学材料工程（085601） +7 WW.' 2026-03-23 9/450 2026-03-29 12:42 by fmesaito
[考研] 318一志愿吉林大学生物与医药 求调剂 +4 笃行致远. 2026-03-28 4/200 2026-03-29 11:55 by longlotian
[考研] 085600，专业课化工原理，321分求调剂 +5 大馋小子 2026-03-28 5/250 2026-03-29 08:56 by qingfeng258
[考研] 329求调剂 +7 星野? 2026-03-26 7/350 2026-03-29 06:43 by 544594351
[考研] 081200-11408-276学硕求调剂 +6 崔wj 2026-03-26 6/300 2026-03-29 01:11 by hanserlol
[考研] 394求调剂 +3 好事多磨静候佳� 2026-03-26 5/250 2026-03-28 14:24 by 唐沐儿
[考研] 材料与化工（0856）304求B区调剂 +8 邱gl 2026-03-27 8/400 2026-03-28 12:42 by 唐沐儿
[考研] 材料277求调剂 +7 min3 2026-03-24 7/350 2026-03-28 11:39 by xuxiang
[考研] 药学105500求调剂 +3 Ssun。。 2026-03-28 3/150 2026-03-28 11:24 by lxf170613
[考研] 311求调剂 +3 希望上岸阿小杨 2026-03-23 3/150 2026-03-28 07:57 by 热情沙漠
[考研] 266求调剂 +11 阳阳哇塞 2026-03-27 12/600 2026-03-27 17:56 by yu221
[考研] 一志愿吉大071010，316分求调剂 +3 xgbiknn 2026-03-27 3/150 2026-03-27 10:36 by guoweigw
[考研] 一志愿陕师大生物学071000，298分，求调剂 +5 SYA！ 2026-03-23 5/250 2026-03-27 09:29 by 不吃魚的貓
[考研] 调剂求收留 +7 果然有我 2026-03-26 7/350 2026-03-27 00:26 by wxiongid
[考研] 0703化学求调剂 +3 丹青奶盖 2026-03-26 5/250 2026-03-26 20:11 by macy2011
[考研] 07化学303求调剂 +5 睿08 2026-03-25 5/250 2026-03-25 22:46 by 418490947
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研经历 +6 cccchenso 2026-03-23 6/300 2026-03-25 22:25 by 544594351
[考博] 26申博自荐 +3 whh869393 2026-03-24 3/150 2026-03-24 09:55 by 21018060
[考研] 333求调剂 +3 ALULU4408 2026-03-23 3/150 2026-03-23 19:04 by macy2011
信息提示
关闭
请填处理意见
关闭