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

24小时热门版块排行榜    

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

mika

金虫 (正式写手)

纳米宝贝

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

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

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

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

» 猜你喜欢

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

mika

金虫 (正式写手)

纳米宝贝
引用回帖:
2楼: Originally posted by gfei2008 at 2012-05-02 23:37:23:
有个简单的指令你试试：grep "POSITION" -n OUTCAR
这个n表示你的原子数目，你把得到的数据存到一个文件下面，然后用ultraedit 去掉一些不用的提示就可以得到你想要的数据。
例如：grep "POSITI ...

谢谢，我去试试。
回复此楼
8楼2012-05-03 16:30:16
已阅   回复此楼   关注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 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 342求调剂 +12 Mary Keen 2026-03-28 13/650 2026-04-01 21:02 by 流情牧豪
[考研] 307分求调剂 +14 (o~o) 2026-03-31 15/750 2026-04-01 20:43 by longlotian
[考研] 070300一志愿211，312分求调剂院校 +14 小黄鸭宝 2026-03-30 14/700 2026-04-01 20:19 by 赖春艳
[考研] 一志愿北交大材料工程总分358 +7 cs0106 2026-04-01 8/400 2026-04-01 18:34 by 记事本2026
[考研] 安徽大学专硕生物与医药专业(086000)324分，英语已过四六级，六级521，求调剂 +10 美味可乐鸡翅 2026-03-26 12/600 2026-04-01 11:17 by syh9288
[考研] 267求调剂 +13 uiybh 2026-03-31 13/650 2026-04-01 10:25 by 探123
[考研] 0710生物学求调剂！ +6 叙述文 2026-03-31 6/300 2026-04-01 09:39 by JourneyLucky
[考研] 085701环境工程，267求调剂 +17 minht 2026-03-26 17/850 2026-04-01 09:11 by xiayizhi
[考研] 080200学硕，机械工程专业277分，求带走！ +4 瓶子PZ 2026-03-31 4/200 2026-03-31 20:16 by vgtyfty
[考研] 一志愿北京科技大学材料学硕328分求调剂 +4 1段时间 2026-03-31 5/250 2026-03-31 19:56 by wxiongid
[考研] 考研调剂 +9 小蜡新笔 2026-03-29 10/500 2026-03-31 19:52 by Dyhoer
[考研] 已决定调剂院校 +8 JKSOIID 2026-03-26 8/400 2026-03-31 19:51 by mg1014
[考研] 一志愿南昌大学324求调剂 +6 hanamiko 2026-03-29 6/300 2026-03-31 16:35 by hypershenger
[考研] 343求调剂 +8 爱羁绊 2026-03-28 8/400 2026-03-31 16:12 by 不吃魚的貓
[考研] 266分，求材料冶金能源化工等调剂 +8 哇呼哼呼哼 2026-03-27 10/500 2026-03-31 13:35 by Huaxue_Wang
[考研] 105500药学求调剂，一志愿山东大学药学，348分 +3 gr哈哈哈 2026-03-28 3/150 2026-03-30 18:56 by 源_2020
[考研] 305求调剂 +8 RuiFairyrui 2026-03-28 8/400 2026-03-29 08:22 by fmesaito
[考研] 394求调剂 +3 好事多磨静候佳� 2026-03-26 5/250 2026-03-28 14:24 by 唐沐儿
[考研] 340求调剂 +5 jhx777 2026-03-27 5/250 2026-03-28 04:18 by fmesaito
[考研] 一志愿上海理工能源动力（085800）310分求调剂 +3 zhangmingc 2026-03-27 4/200 2026-03-27 19:01 by 给你你注意休息
信息提示
关闭
请填处理意见
关闭