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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】如何根据cell文件写poscar
121/2返回列表12下一页
查看: 1802  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

xiaowu787

木虫 (正式写手)

[交流] 【求助】如何根据cell文件写poscar 已有3人参与

刚开始学ASP,看论坛上高手说过可以通过MS写OSCAR。从MS 导出cell文件,还是不知道该如何写出基矢？请高手指点一下。原子的个数就是BLOCK POSITIONS_FRAC中的原子个数吗？谢谢！！

Al2O3
PBC    4.7590    4.7590   12.9910   90.0000   90.0000  120.0000 (P1)

cell文件如下：

%BLOCK LATTICE_CART
       2.747609931073429       0.000000000000000       4.330333333333333
      -1.373804965536714       2.379500000000001       4.330333333333333
      -1.373804965536714      -2.379500000000000       4.330333333333333
%ENDBLOCK LATTICE_CART

%BLOCK POSITIONS_FRAC
  O   0.5559999999999998  -0.0559999999999999   0.2500000000000000
  O   0.2500000000000000   0.5559999999999998  -0.0559999999999999
  O  -0.0559999999999999   0.2500000000000000   0.5559999999999998
  O  -0.5559999999999998   0.0559999999999999  -0.2500000000000000
  O  -0.2500000000000000  -0.5559999999999998   0.0559999999999999
  O   0.0559999999999999  -0.2500000000000000  -0.5559999999999998
Al   0.3519999999999999   0.3519999999999999   0.3519999999999999
Al   0.1480000000000001   0.1480000000000001   0.1480000000000001
Al  -0.3519999999999999  -0.3519999999999999  -0.3519999999999999
Al   0.8520000000000000   0.8520000000000000   0.8520000000000000
%ENDBLOCK POSITIONS_FRAC
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-11-02 22:53:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

bingmou

金虫 (著名写手)
★ ★
xiaowu787(金币+5):谢谢答复！ 2010-11-03 11:35:12
cenwanglai(金币+2):谢谢回复！ 2010-11-03 13:22:21
示例：
Al2O3
1.0
      2.747609931073429       0.000000000000000       4.330333333333333
      -1.373804965536714       2.379500000000001       4.330333333333333
      -1.373804965536714      -2.379500000000000       4.330333333333333
6 4 !前面为O的个数，后面为Al的个数，这个顺序要和POTCAR对应
Direct !MS的都是Direct
   0.5559999999999998  -0.0559999999999999   0.2500000000000000
   0.2500000000000000   0.5559999999999998  -0.0559999999999999
  -0.0559999999999999   0.2500000000000000   0.5559999999999998
  -0.5559999999999998   0.0559999999999999  -0.2500000000000000
  -0.2500000000000000  -0.5559999999999998   0.0559999999999999
   0.0559999999999999  -0.2500000000000000  -0.5559999999999998
   0.3519999999999999   0.3519999999999999   0.3519999999999999
   0.1480000000000001   0.1480000000000001   0.1480000000000001
  -0.3519999999999999  -0.3519999999999999  -0.3519999999999999
   0.8520000000000000   0.8520000000000000   0.8520000000000000
赞一下(2人) 回复此楼
3楼2010-11-03 11:17:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

xiaowu787

木虫 (正式写手)
这是从其他地方查到的信息,请高手指点一下，是不是问题太简单了？可是我没做过这方面的。

•Prototype: αAl2O3
•Pearson Symbol: hR10
•Strukturbericht Designation: D51
•Space Group: R3c (Cartesian and lattice coordinate listings available)
•Number: 167
•Reference: Villars and Calvert, Pearson's Handbook, Vol. I, pp. 970-1
•Primitive Vectors:
A1  =  b X + a Y + a Z
A2  =  a X + b Y + a Z
A3  =  a X + a Y + b Z
赞一下 回复此楼
2楼2010-11-03 11:07:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowu787

木虫 (正式写手)
引用回帖:
Originally posted by bingmou at 2010-11-03 11:17:57:
示例：
Al2O3
1.0
      2.747609931073429       0.000000000000000       4.330333333333333
      -1.373804965536714       2.379500000000001       4.330333333333333
      -1.373804965536714       ...

由于基失包含了晶格常数，只需把晶格缩放系数设为1是吧。再问一下，晶格缩放系数设为1和通过转换变成较大的数有区别吗？谢谢！
赞一下 回复此楼
4楼2010-11-03 11:34:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主
★ ★ ★
bingmou(金币+3):感谢分享 2010-11-03 14:06:48
xiaowu787(金币+5):谢谢分享！！ 2010-11-03 15:45:25
#Script for Linux. Font Code: UTF-8
#把该文件保存到/usr/bin/目录下，必要时运行dos2unix和chmod
#删除原有poscarfromcell,删除前询问
#*.cell为MS running后立即停止中生成的隐藏文件，显示隐藏文件夹后就可以看到。把该文件放到linux下工作目录内。
#运行本脚本前，需要删除.cell文件%ENDBLOCK POSITIONS_FRAC后面的所有内容,如有必要，运行dos2unix *.cell
#运行时确保当前目录下仅有一个*.cell文件。
#在命令提示符下输入cell2poscar回车，生成文件为poscarfromcell.
#有好的修改建议或脚本请发邮件给：cenwanglai@163.com

rm -i poscarfromcell

echo "SYSTEM COMMENTS">>poscarfromcell
echo 1 >>poscarfromcell

#查找并输出格矢坐标
cat *.cell | grep -e '^[ ]\+[0-9]\+'>>poscarfromcell

#数氧原子个数
L=$(grep -cw "O" *.cell)

#数铈原子个数
M=$(grep -cw "Ce" *.cell)

#数X原子个数，如果需要，删掉下一行开始的#，并替换“X”
#N=$(grep -cw "X" *.cell)

#输入各类原子个数,如果需要，下行中增加$N
echo $L $M >> poscarfromcell

echo "Selective Dynamic" >>poscarfromcell
echo "Direct" >>poscarfromcell

#找氧原子O坐标。根据实际情况替换"O"
cat *.cell | awk '$1=="O" {print $2,$3,$4,"T T T"}'>>poscarfromcell

#找Ce原子坐标。根据实际情况替换"Ce"
cat *.cell | awk '$1=="Ce" {print $2,$3,$4,"T T T"}'>>poscarfromcell

#找X原子坐标，如果需要，删掉下一行开头的#,根据实际情况替换"X"
#cat *.cell | awk '$1=="X" {print $2,$3,$4,"T T T"}'>>poscarfromcell
赞一下(1人) 回复此楼
5楼2010-11-03 13:26:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主
读一下这个脚本，根据你的情况修改一下就可以。
赞一下 回复此楼
6楼2010-11-03 13:28:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)
★
zzy870720z(金币+1):谢谢提示 2010-11-03 14:14:21
引用回帖:
Originally posted by xiaowu787 at 2010-11-03 11:34:49:



由于基失包含了晶格常数，只需把晶格缩放系数设为1是吧。再问一下，晶格缩放系数设为1和通过转换变成较大的数有区别吗？谢谢！

在计算能带的时候取K点的时候小心一点，理解好K点的坐标就行了
赞一下(1人) 回复此楼
7楼2010-11-03 14:07:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowu787

木虫 (正式写手)
引用回帖:
Originally posted by cenwanglai at 2010-11-03 13:26:17:
#Script for Linux. Font Code: UTF-8
#把该文件保存到/usr/bin/目录下，必要时运行dos2unix和chmod
#删除原有poscarfromcell,删除前询问
#*.cell为MS running后立即停止中生成的隐藏文件，显示隐藏文件夹后就 ...

这个脚本没有声明shell,脚本运行不必拷到/usr/bin/目录下，一般只需要更改一下运行权限就可以了。正准备测试一下,谢谢
赞一下 回复此楼
8楼2010-11-03 16:14:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by xiaowu787 at 2010-11-03 16:14:15:



这个脚本没有声明shell,脚本运行不必拷到/usr/bin/目录下，一般只需要更改一下运行权限就可以了。正准备测试一下,谢谢

呵呵，刚学习linux的时候编写的一个脚本，很久没有检查了。幸苦！你测试后帮修改再发一个吧。
赞一下(1人) 回复此楼
9楼2010-11-03 16:17:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowu787

木虫 (正式写手)
★
bingmou(金币+1):感谢交流 2010-11-03 19:54:13
引用回帖:
Originally posted by cenwanglai at 2010-11-03 16:17:16:



呵呵，刚学习linux的时候编写的一个脚本，很久没有检查了。幸苦！你测试后帮修改再发一个吧。

谢谢了，简单修改了一下测试通过，保存后改一下权限就会用了。
CODE:
#!/bin/sh
#Script for Linux. Font Code: UTF-8
#把该文件保存到/usr/bin/目录下，必要时运行dos2unix和chmod
#删除原有poscarfromcell,删除前询问
#*.cell为MS running后立即停止中生成的隐藏文件，显示隐藏文件夹后就可以看到。把该文件放到linux下工作目录内。
#运行本脚本前，需要删除.cell文件%ENDBLOCK POSITIONS_FRAC后面的所有内容,如有必要，运行dos2unix *.cell
#运行时确保当前目录下仅有一个*.cell文件。
#在命令提示符下输入cell2poscar回车，生成文件为poscarfromcell.
#有好的修改建议或脚本请发邮件给：[email]cenwanglai@163.com[/email]

rm -i poscarfromcell

echo "SYSTEM COMMENTS">>poscarfromcell
echo 1 >>poscarfromcell

#查找并输出格矢坐标
cat *.cell | grep -e '^[ ]\+[0-9]\+'>>poscarfromcell

#数氧原子个数
L=$(grep -cw "O" *.cell)

#数铈原子个数
M=$(grep -cw "Ce" *.cell)

#数X原子个数，如果需要，删掉下一行开始的#，并替换“X”
#N=$(grep -cw "X" *.cell)

#输入各类原子个数,如果需要，下行中增加$N
echo $L $M >> poscarfromcell

echo "Selective Dynamic" >>poscarfromcell
echo "Direct" >>poscarfromcell

#找氧原子O坐标。根据实际情况替换"O"
cat *.cell |awk '$1=="O" {printf "%20.16f %20.16f %20.16f %s",$2,$3,$4,"T T T\n"}'>>poscarfromcell
#找Ce原子坐标。根据实际情况替换"Ce"
cat *.cell |awk '$1=="Ce" {printf "%20.16f %20.16f %20.16f %s",$2,$3,$4,"T T T\n"}'>>poscarfromcell

#找X原子坐标，如果需要，删掉下一行开头的#,根据实际情况替换"X"
#cat *.cell | awk '$1=="X" {print $2,$3,$4,"T T T"}'>>poscarfromcell

暂时还有读基失时，负号的读不出来

[ Last edited by xiaowu787 on 2010-11-3 at 19:28 ]
赞一下(1人) 回复此楼
10楼2010-11-03 18:27:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaowu787 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料080500调剂求收留 +6 一颗meteor 2026-03-13 6/300 2026-03-20 10:41 by EBSD
[考研] 271材料工程求调剂 +7 .6lL 2026-03-18 7/350 2026-03-20 09:10 by xingguangj
[考研] 材料学硕318求调剂 +5 February_Feb 2026-03-19 5/250 2026-03-19 23:51 by 23Postgrad
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[考研] 085600材料与化工求调剂 +6 绪幸与子 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
[教师之家] 焦虑 +9 水冰月月野兔 2026-03-13 13/650 2026-03-19 09:50 by otani
[考研] 328求调剂，英语六级551，有科研经历 +3 生物工程调剂 2026-03-17 7/350 2026-03-18 20:41 by Wangjingyue
[考研] 311求调剂 +11 冬十三 2026-03-15 12/600 2026-03-18 14:36 by 星空星月
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 0854，计算机类招收调剂 +3 胡辣汤放糖 2026-03-15 6/300 2026-03-18 12:09 by 上岸上岸……..
[考研] 299求调剂 +5 △小透明* 2026-03-17 5/250 2026-03-18 11:49 by 尽舜尧1
[考研] 生物学071000 329分求调剂 +3 我爱生物生物爱� 2026-03-17 3/150 2026-03-18 10:12 by macy2011
[考研] 278求调剂 +5 烟火先于春 2026-03-17 5/250 2026-03-18 08:43 by 星空星月
[考研] 268求调剂 +8 一定有学上- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[考研] 301求调剂 +4 A_JiXing 2026-03-16 4/200 2026-03-17 17:32 by ruiyingmiao
[考研] 一志愿南京大学，080500材料科学与工程，调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
信息提示
关闭
请填处理意见
关闭