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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】如何根据cell文件写poscar
121/2返回列表12下一页
查看: 1830  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 273求调剂 +42 麦小叮当 2026-04-06 49/2450 2026-04-10 08:08 by kangsm
[考研] 已调剂 +18 柴郡猫_ 2026-04-09 19/950 2026-04-09 22:10 by 柴郡猫_
[考研] 一志愿中国科学院上海有机所，有机化学356分找调剂 +11 Nadiums 2026-04-09 11/550 2026-04-09 18:04 by lijunpoly
[考研] 312求调剂 +3 李鸿飞飞 2026-04-06 3/150 2026-04-09 17:32 by wp06
[考研] 332，085601求调剂 +12 ydfyh 2026-04-09 14/700 2026-04-09 17:28 by wp06
[考研] 求调剂材料科学与工程一志愿985初试365分 +5 材化李可 2026-04-08 5/250 2026-04-09 17:00 by Lilly_Li
[考研] 085801 总分275 本科新能源 求调剂 +8 bradoner 2026-04-08 9/450 2026-04-09 13:43 by only周
[考研] 0860004 求调剂 309分 +6 Yin DY 2026-04-09 6/300 2026-04-09 10:19 by 啊李999
[考研] 285求调剂 +20 哦呦呼o 2026-04-04 20/1000 2026-04-08 22:23 by yutian743
[考研] 一志愿211，化学学硕，310分，本科重点双非，求调剂 +10 努力奋斗112 2026-04-07 10/500 2026-04-08 15:01 by screening
[考研] 081200-11408-367学硕求调剂 +4 1_2_3111 2026-04-06 4/200 2026-04-07 08:13 by jp9609
[考研] 285求调剂 +15 哦呦呼o 2026-04-04 17/850 2026-04-06 23:02 by chenzhimin
[考研] 一志愿北交大材料工程总分358求调剂 +10 cs0106 2026-04-05 12/600 2026-04-06 19:41 by 无际的草原
[考研] 085405软件工程301分求调剂，专硕可跨专业，四六级已过 +3 静静想想 2026-04-05 3/150 2026-04-06 15:23 by nepu_uu
[考研] 材料工程310专硕调剂 +14 捞捞我…. 2026-04-04 15/750 2026-04-06 14:18 by lqwchd
[考研] 0857大类环境工程B区求调剂 +3 龚禹铭 2026-04-05 3/150 2026-04-06 10:22 by 蓝云思雨
[考研] 0860 求调剂 一志愿国科大 348 分 +3 WiiiP 2026-04-03 3/150 2026-04-05 17:43 by Ecowxq666！
[考研] 一志愿上海海洋大学083200食品学硕，求调剂，接受其他专业083200 +4 what张 2026-04-04 5/250 2026-04-05 14:07 by chw1980_0
[考研] 328分调剂 +6 门men 2026-04-04 6/300 2026-04-05 13:40 by imissbao
[考研] 281求调剂 +10 aaawhy 2026-04-03 10/500 2026-04-03 21:42 by lbsjt
信息提示
关闭
请填处理意见
关闭