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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » QE计算relax，怎样迅速得到计算后的原子位置？
51/1返回列表
查看: 2428  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

小白flame

新虫 (初入文坛)

[求助] QE计算relax，怎样迅速得到计算后的原子位置？ 已有2人参与

刚刚上手QE，计算relax以后得到的输出文件内容太多，目前是手动把最后的原子位置提取出来，实在太麻烦，而且以后想做批量处理，肯定不能挨个手动提取。请问有没有什么好方法快速得到优化后的原子位置？
回复此楼

» 猜你喜欢
1楼 2016-11-01 10:57:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

magicmonk

至尊木虫 (著名写手)
引用回帖:
3楼: Originally posted by souledge at 2016-11-01 13:50:01
对于calculation = 'relax'的，假设输出文件为file.out，可以用这样的bash脚本：

#!/bin/bash
  NAT=`grep 'number of atoms/cell' file.out | awk '{print $5}'`
  LINE=$NAT
  let "$((LINE++))" ...

请问要将输出的坐标替换原来输入文件的坐标，应该怎么定向？
赞一下 回复此楼
我们最深的恐惧不是我们能力不够，我们最深的恐惧是能力超越了极限。
8楼2017-06-12 20:52:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

liqizuiyang

木虫 (著名写手)
★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
漫天飘雪: 金币+2, 谢谢交流 2016-11-01 22:18:59
之前写过一个类似的程序，楼主参考下。

这个程序只适用于手写基矢（ibrav = 0）的情形，不过不难改写。
赞一下(1人) 回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : out2pd.py
    • 2016-11-01 11:23:40, 3.49 K
2楼2016-11-01 11:25:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)
★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
漫天飘雪: 金币+2, 谢谢交流 2016-11-01 22:19:06
对于calculation = 'relax'的，假设输出文件为file.out，可以用这样的bash脚本：
CODE:
#!/bin/bash
  NAT=`grep 'number of atoms/cell' file.out | awk '{print $5}'`
  LINE=$NAT
  let "$((LINE++))"
  grep -A$(NAT) "ATOMIC_POSITIONS' file.out | tail -$(LINE)

如果需要输出，最后一句加个定向就好了~
赞一下(1人) 回复此楼
思想重于技巧，内涵重于表象
3楼2016-11-01 13:50:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小白flame

新虫 (初入文坛)
引用回帖:
3楼: Originally posted by souledge at 2016-11-01 13:50:01
对于calculation = 'relax'的，假设输出文件为file.out，可以用这样的bash脚本：

#!/bin/bash
  NAT=`grep 'number of atoms/cell' file.out | awk '{print $5}'`
  LINE=$NAT
  let "$((LINE++))" ...

多谢，还请问QE可以将relax的输出转化为scf计算的输入不？就是把两个先后的计算相连，直接得到最终结果
赞一下 回复此楼
4楼2016-11-01 15:14:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 本人女孩 +8 吼吼， 2026-04-10 10/500 2026-04-17 11:28 by 4444si
[考研] 279求调剂 +13 张番茄不炒蛋 2026-04-11 13/650 2026-04-17 10:38 by cuisz
[考研] 300求调剂 +11 橙a777 2026-04-15 11/550 2026-04-16 22:43 by cfdbai
[考研] 22专硕求调剂 +10 haoyun上岸 2026-04-11 12/600 2026-04-16 22:21 by 猪会飞
[考研] 急需调剂 +7 绝不放弃22 2026-04-15 7/350 2026-04-16 22:09 by SUSE_CL
[考研] 327求调剂 +26 Xxjc1107. 2026-04-13 29/1450 2026-04-16 10:52 by Espannnnnol
[考研] 0854调剂 +13 长弓傲 2026-04-12 16/800 2026-04-15 13:45 by fenglj492
[考研] 310求调剂 +16 666真好 2026-04-11 18/900 2026-04-15 13:28 by 黑科技矿业
[考研] 调剂求收留 +34 果然有我 2026-04-10 35/1750 2026-04-15 13:05 by 西北望—风沙
[考研] 272分材料子求调剂 +41 Loy0361 2026-04-10 54/2700 2026-04-14 18:00 by lhj2009
[考研] 366求调剂 +11 不知名的小卅 2026-04-11 11/550 2026-04-14 15:50 by zs92450
[考研] 105500药学求调剂 +4 x_skys 2026-04-12 4/200 2026-04-14 13:37 by rndfc
[考研] 305求调剂 +8 玛卡巴卡boom 2026-04-11 8/400 2026-04-14 09:04 by pengliang8036
[考研] 293求调剂 +16 我爱高数高数爱� 2026-04-12 18/900 2026-04-13 21:47 by 学员JpLReM
[考研] 材料考研调剂 +29 云木达达 2026-04-11 31/1550 2026-04-13 13:32 by lyh鲁老师
[考研] 一志愿浙大生物325分求调剂 +9 zysheng 2026-04-12 9/450 2026-04-12 22:31 by yuyin1233
[考研] 344 材料专业 求调剂211 无地域要求 +8 hualkop 2026-04-11 8/400 2026-04-12 22:24 by fqwang
[考研] 调剂 +5 文道星台 2026-04-11 5/250 2026-04-11 15:01 by 凯凯要变帅
[考研] 0854调剂 +5 音像店听花鼓戏 2026-04-10 5/250 2026-04-11 10:49 by qingpingzhu
[考研] 263能源动力专硕求调剂 +3 加大号饭盒袋 2026-04-10 3/150 2026-04-10 22:23 by 286640313
信息提示
关闭
请填处理意见
关闭