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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】如何根据cell文件写poscar
121/2返回列表12下一页
查看: 1803  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 工科材料085601 279求调剂 +7 困于星晨 2026-03-17 9/450 2026-03-20 17:38 by 无懈可击111
[考研] 一志愿中南化学337求调剂 +4 niko- 2026-03-19 5/250 2026-03-20 13:49 by 促天成
[考研] 能源材料化学课题组招收硕士研究生8-10名 +5 脱颖而出 2026-03-16 14/700 2026-03-20 09:30 by kkcoco25
[考研] 求调剂 +3 暗涌afhb 2026-03-16 3/150 2026-03-20 00:28 by 河南大学校友
[考研] 一志愿中国海洋大学，生物学，301分，求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考研] 307求调剂 +9 冷笙123 2026-03-17 9/450 2026-03-19 22:44 by 学员8dgXkO
[考研] 266求调剂 +5 阳阳哇塞 2026-03-14 10/500 2026-03-19 15:08 by 阳阳哇塞
[考研] 求调剂，一志愿:南京航空航天大学大学 ，080500材料科学与工程学硕，总分289分 +3 @taotao 2026-03-19 3/150 2026-03-19 14:07 by peike
[考研] 化学求调剂 +3 临泽境llllll 2026-03-17 4/200 2026-03-19 13:59 by houyaoxu
[考研] 0703化学调剂 +5 pupcoco 2026-03-17 8/400 2026-03-19 13:58 by houyaoxu
[考研] 一志愿福大288有机化学，求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
[考研] 一志愿中海洋材料工程专硕330分求调剂 +7 小材化本科 2026-03-18 7/350 2026-03-19 10:46 by Linda Hu
[考研] 330求调剂 +3 小材化本科 2026-03-18 3/150 2026-03-18 21:55 by 无懈可击111
[考研] 一志愿武理材料305分求调剂 +5 想上岸的鲤鱼 2026-03-18 6/300 2026-03-18 17:53 by 无际的草原
[考研] 0854可跨调剂，一作一项核心论文五项专利，省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[考研] 321求调剂 +5 大米饭！ 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 0703化学调剂 290分有科研经历，论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
信息提示
关闭
请填处理意见
关闭