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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 从outcar中提取原子坐标的脚本
51/1返回列表
查看: 5252  |  回复: 14
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 2 个 1ST强帖 ，点击这里进行查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

mika

金虫 (正式写手)

纳米宝贝

[求助] 从outcar中提取原子坐标的脚本

想从outcar里面提取每一步自洽后指定原子的xyz坐标，请问有没有相关的脚本啊，谢谢
回复此楼

» 收录本帖的淘帖专辑推荐

VASP 第一性原理相关文档 第一性原理计算辅助工具 第一性原理和电化学
纳米及第一性原理方向 lywiailyw vasp从结构优化到光学性质 phystt

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-05-02 20:13:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fzx2008

荣誉版主 (著名写手)

优秀版主优秀版主
dxf_dxf_dxf 辛苦了！
赞一下 回复此楼
5楼2012-05-03 13:35:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答

gfei2008

木虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
WDD880227: 金币+2, 感谢您的详细解答~~很有帮助 2012-05-03 08:24:54
有个简单的指令你试试：grep "POSITION" -n OUTCAR
这个n表示你的原子数目，你把得到的数据存到一个文件下面，然后用ultraedit 去掉一些不用的提示就可以得到你想要的数据。
例如：grep "POSITION" -n OUTCAR>data
然后得到的坐标和受力在这个data文件中，然后用ultraedit处理一下就可以啦。
赞一下(1人) 回复此楼
2楼2012-05-02 23:37:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxf_dxf_dxf

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
fzx2008: 金币+5, 1ST强帖+1, 谢谢指教，详细 2012-05-03 13:04:02
你好！我用官方的例子来解释。
首先，vasp进行几何优化时（IBRION=2 or 3，ISIF=2，3，5等）时，才会移动原子。vasp会根据ISIF的值改变晶格常数或移动原子位置，求出体系能量最小或受力最小对应的几何结构。当设置了IBRION和ISIF参数之后，vasp实际会进行三层循环（见mannual09版7.1节Fig3）：
最内层循环，是用来矩阵对角化，求解薛定谔方程。
vasp会先给出一个初始解，带入薛定谔方程。
如果没有达到收敛条件，则修改后再次带入方程
中层循环，是用来求出自洽的电荷密度。
vasp会先猜测一个初始的电荷密度，由此构造薛定谔方程，求解得到波函数，再根据新的波函数算出（输出的）电荷密度，和（输入的）初始电荷密度对比。
如果差别没有达到收敛条件，则二者混合后作为新的输入电荷密度，重复上一步。
在mannual6.7节NWRITE-tag中，这一步对应于“electronic step”
在OUTCAR中，对应类似与“----------- Iteration    3(   2)  ----------------”中的“2”
外层循环，是用来求出能量最低、受力最小的几何结构。
vasp会根据初始POSCAR文件求出体系的总能和受力情况。根据算法改变晶格常数或离子位置，得到新的结构，再次求出总能和受力。
如果差别或绝对值未达到收敛条件，则继续改变晶格常数或离子位置。
在mannual6.7节NWRITE-tag中，这一步对应于“ion step”
在OUTCAR中，对应类似与“----------- Iteration    3(   2)  ----------------”中的“3”

如果在INCAR文件中没有设置NWRITE参数，则在OUTCAR文件中可以看到每一次ion step结束后vasp会给出优化情况：

vasp.5.2.2 15Apr09 complex
executed on             LinuxIFC date 2012.05.03  12:04:41
running on    8 nodes
distr:  one band on    1 nodes,    8 groups
--------------------------------------------------------------------------------------------------------
INCAR:
   ...
--------------------------------------------------------------------------------------------------------
...
...
----------------------------------------- Iteration    1(   1)  ---------------------------------------
    POTLOK:  cpu time    0.01: real time    0.01
    ...
----------------------------------------- Iteration    1(   2)  ---------------------------------------
     ...
--------------------------------------------------------------------------------------------------------
...
------------------------aborting loop because EDIFF is reached----------------------------------------
...
VOLUME and BASIS-vectors are now :
-----------------------------------------------------------------------------
...
direct lattice vectors reciprocal lattice vectors 0.000000000 2.750000000 2.750000000 -0.181818182 0.181818182 0.181818182 2.750000000 0.000000000 2.750000000 0.181818182 -0.181818182 0.181818182 2.750000000 2.750000000 0.000000000 0.181818182 0.181818182 -0.181818182
...
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      4.81250      4.81250      4.81250         0.000000      0.000000      0.000000
      0.68750      0.68750      0.68750         0.000000      0.000000      0.000000
-----------------------------------------------------------------------------------
...
----------------------------------------- Iteration    2(   1)  ---------------------------------------
  ...
------------------------ aborting loop because EDIFF is reached----------------------------------------
...
VOLUME and BASIS-vectors are now :
-----------------------------------------------------------------------------
...

direct lattice vectors                 reciprocal lattice vectors
     0.000000000  2.721754370  2.721754370    -0.183705042  0.183705042  0.183705042
     2.721754370  0.000000000  2.721754370     0.183705042 -0.183705042  0.183705042
     2.721754370  2.721754370  0.000000000     0.183705042  0.183705042 -0.183705042
...
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      4.76307      4.76307      4.76307         0.000000      0.000000      0.000000
      0.68044      0.68044      0.68044         0.000000      0.000000      0.000000
-----------------------------------------------------------------------------------

...
----------------------------------------- Iteration    3(   1)  ---------------------------------------
...
------------------------ aborting loop because EDIFF is reached ----------------------------------------
...
VOLUME and BASIS-vectors are now :
-----------------------------------------------------------------------------
...

  direct lattice vectors                 reciprocal lattice vectors
     0.000000000  2.732453619  2.732453619    -0.182985723  0.182985723  0.182985723
     2.732453619  0.000000000  2.732453619     0.182985723 -0.182985723  0.182985723
     2.732453619  2.732453619  0.000000000     0.182985723  0.182985723 -0.182985723
...
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      4.78179      4.78179      4.78179         0.000000      0.000000      0.000000
      0.68311      0.68311      0.68311         0.000000      0.000000      0.000000
-----------------------------------------------------------------------------------
...
--------------------------------------------------------------------------------------------------------
    POTLOK:  cpu time    0.01: real time    0.01
--------------------------------------------------------------------------------------------------------

--------------------------------------------------------------------------------------------------------
reached required accuracy - stopping structural energy minimisation
...

由上可知，整个优化过程，外层循环执行了3次，共产生了四个结构：
初始结构
引用回帖:
原胞基矢
0.000000000 2.750000000 2.750000000
2.750000000 0.000000000 2.750000000
2.750000000 2.750000000 0.000000000
原子位置（相对于原胞基矢）
0.87500000 0.87500000 0.87500000
0.12500000 0.12500000 0.12500000
原子位置（相对于笛卡尔系）
4.81250000 4.81250000 4.81250000
0.68750000 0.68750000 0.68750000

第一次外层循环结束后的结构
引用回帖:
原胞基矢
     0.000000000  2.750000000  2.750000000
     2.750000000  0.000000000  2.750000000
     2.750000000  2.750000000  0.000000000
原子位置（相对于笛卡尔系）
      4.81250      4.81250      4.81250
      0.68750      0.68750      0.68750

第二次外层循环结束后的结构
引用回帖:
原胞基矢
     0.000000000  2.721754370  2.721754370
     2.721754370  0.000000000  2.721754370
     2.721754370  2.721754370  0.000000000
原子位置（相对于笛卡尔系）
     4.76307      4.76307      4.76307
      0.68044      0.68044      0.68044

第三次外层循环结束后的结构
引用回帖:
原胞基矢
     0.000000000  2.732453619  2.732453619
     2.732453619  0.000000000  2.732453619
     2.732453619  2.732453619  0.000000000
原子位置（相对于笛卡尔系）
      4.78179      4.78179      4.78179
      0.68311      0.68311      0.68311
赞一下(4人) 回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : 2_4_diamondSi_vol_rex.tar.bz2
    • 2012-05-03 12:15:48, 1.39 M
3楼2012-05-03 12:59:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxf_dxf_dxf

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
fzx2008: 1ST强帖+1, 谢谢指教 2012-05-03 13:34:14
uuv2010: 金币+10, 多谢分享~ 2012-05-04 18:36:54
由此给出思路及相应脚本：
  首先，要知道体系有多少个原子。这个可以从POSCAR中读取。
  然后，要告诉脚本，你关心的是哪个原子。我们把这个原子在POSCAR中的序号写到一个叫“in”的文件里面，让脚本去读就好了！
  再然后，利用grep从OUTCAR中读取我们关心的信息，并加以处理。
  依据本例，体系共有2个原子，假设我们关心POSCAR中的第二个原子。
  首先，让grep去找关键字“POSITION”，我们关心的数字在这个关键字之后的1（“---"一行）+2（2个原子）行。所以用
CODE:
grep -A3 "POSITION" OUTCAR

命令，把冗长的OUTCAR简化成如下的模样：
引用回帖:
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      4.81250      4.81250      4.81250         0.000000      0.000000      0.000000
      0.68750      0.68750      0.68750         0.000000      0.000000      0.000000
--
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      4.76307      4.76307      4.76307         0.000000      0.000000      0.000000
      0.68044      0.68044      0.68044         0.000000      0.000000      0.000000
--
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      4.78179      4.78179      4.78179         0.000000      0.000000      0.000000
      0.68311      0.68311      0.68311         0.000000      0.000000      0.000000

这样还不够。我们要把”POSITON“行、”-------“行、分隔符”--“删掉：
这3种类型的行的共同特点是，如果以空格（默认的）为分隔符，则最多只有2项，第三项为空。于是用awk语句的$3>0来过滤掉它们，然后用print $1“ ”$2" "$3打印有用的前三个数字（我们不关心受力）
CODE:
grep -A3 "POSITION" OUTCAR | awk '$3>0 {print $1“  ”$2"  "$3}'

于是结果变成这个样子：
引用回帖:
4.81250  4.81250  4.81250
0.68750  0.68750  0.68750
4.76307  4.76307  4.76307
0.68044  0.68044  0.68044
4.78179  4.78179  4.78179
0.68311  0.68311  0.68311

这样还不够。我们需要的只是我们关心的那一个原子。怎么把它调出来呢？
由于我们保留了所有的原子的位置，因此到目前为止，输出的行数应该=体系原子数x外层循环次数。我们只需要找到行数对体系原子数取余数后结果为我们关心的原子对应的序数的行就行了：
NR%2==0(第2个原子，也就是余数为2，即余数为0)
CODE:
grep -A3 "POSITION" OUTCAR | awk '$3>0 {print $1"  "$2"  "$3}' | awk 'NR%2==0 {print $0}'

结果如下：
引用回帖:
0.68750  0.68750  0.68750
0.68044  0.68044  0.68044
0.68311  0.68311  0.68311

大功告成！
所以语句如下：
引用回帖:
grep -A3 "POSITION" OUTCAR | awk '$3>0 {print $1"  "$2"  "$3}' | awk 'NR%2==0 {print $0}'

当然，你需要明白“A3”中3的含义，还有"NR%2==0"中2和0的含义。否则生搬硬套会得到错误的结果~
赞一下(5人) 回复此楼
4楼2012-05-03 13:26:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 8rmuugja8q 2026-02-22 8/400 2026-02-23 12:22 by alian_214
[基金申请] 体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低 +10 瞬息宇宙 2026-02-20 13/650 2026-02-23 11:23 by holypower
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 khieu8v8m0 2026-02-22 8/400 2026-02-23 09:35 by w4l55oybr1
[教师之家] 为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？ +5 欢乐颂叶蓁 2026-02-21 8/400 2026-02-23 09:19 by 欢乐颂叶蓁
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 w89i99eaeh 2026-02-22 5/250 2026-02-23 08:04 by w4l55oybr1
[博后之家] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 khieu8v8m0 2026-02-22 6/300 2026-02-23 07:59 by w4l55oybr1
[博后之家] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +6 3dfhjxgsh7 2026-02-22 9/450 2026-02-23 07:49 by w4l55oybr1
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 khieu8v8m0 2026-02-22 4/200 2026-02-23 06:46 by jsjzfl
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 khieu8v8m0 2026-02-22 5/250 2026-02-23 06:29 by w4l55oybr1
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 khieu8v8m0 2026-02-22 8/400 2026-02-23 06:24 by w4l55oybr1
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +5 3dfhjxgsh7 2026-02-22 6/300 2026-02-23 02:04 by 5jlh3qtdvx
[教师之家] 版面费该交吗 +7 苹果在哪里 2026-02-22 8/400 2026-02-22 22:37 by otani
[基金申请] 面上可以超过30页吧？ +4 阿拉贡aragon 2026-02-22 4/200 2026-02-22 21:22 by 山西悬空寺空悬�
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 usprnugpzw 2026-02-21 6/300 2026-02-22 19:48 by w89i99eaeh
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 3dfhjxgsh7 2026-02-22 4/200 2026-02-22 16:52 by khieu8v8m0
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 usprnugpzw 2026-02-22 3/150 2026-02-22 16:37 by khieu8v8m0
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 usprnugpzw 2026-02-21 4/200 2026-02-22 16:27 by khieu8v8m0
[基金申请] “人文社科而论，许多学术研究还没有达到民国时期的水平” +4 苏东坡二世 2026-02-18 5/250 2026-02-22 16:07 by liangep1573
[基金申请] 什么是人一生最重要的？ +4 瞬息宇宙 2026-02-21 4/200 2026-02-22 11:44 by huagongfeihu
[基金申请] 今年春晚有几个节目很不错，点赞！ +11 瞬息宇宙 2026-02-16 12/600 2026-02-21 21:14 by lq493392203
信息提示
关闭
请填处理意见
关闭