24小时热门版块排行榜    

查看: 2548  |  回复: 9

黑暗游侠

铜虫 (正式写手)

[求助] 求晶格常数优化小程序(a,b固定,只优化c轴) 已有2人参与

哪位大神有晶格常数优化小程序,给我发一份吧。急需~~~谢谢各位了~~~
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dengc09

禁虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
黑暗游侠: 金币+10, ★★★很有帮助 2014-05-28 15:56:08
本帖内容被屏蔽

2楼2014-05-28 15:40:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

May.Be

新虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
黑暗游侠: 金币+10, ★★★★★最佳答案 2014-05-28 15:55:57
给你个脚本,我就是用这个优化的

» 本帖附件资源列表

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

黑暗游侠

铜虫 (正式写手)

引用回帖:
2楼: Originally posted by dengc09 at 2014-05-28 15:40:06
可以这样吗?INCAR中ISIF=2,循环脚本用这个:
for i in x1 x2 x3......
do
cat >POSCAR<<!
XXXX
1.0
a  0  0
0  b  0
0  0  $i
X
direct
XXX

楼主试一下

请问 那几个X都表示什么?非常感谢
4楼2014-05-28 15:50:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

黑暗游侠

铜虫 (正式写手)

引用回帖:
3楼: Originally posted by May.Be at 2014-05-28 15:47:18
给你个脚本,我就是用这个优化的

请问 里面的原子坐标是对应c轴多大的时候?是第一个c轴取值吗?
5楼2014-05-28 15:54:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

黑暗游侠

铜虫 (正式写手)

引用回帖:
3楼: Originally posted by May.Be at 2014-05-28 15:47:18
给你个脚本,我就是用这个优化的

我懂了,是分数坐标~~~多谢
6楼2014-05-28 16:01:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dengc09

禁虫 (正式写手)

本帖内容被屏蔽

7楼2014-05-28 16:05:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

May.Be

新虫 (初入文坛)

引用回帖:
6楼: Originally posted by 黑暗游侠 at 2014-05-28 16:01:03
我懂了,是分数坐标~~~多谢...

你就换一下这个一行的数值就行for i in  6.76 6.78 6.80 6.82 6.84
这就是你估计的c,算出来后找最低能量对应的数值,其实算出来以后有一个拟合的,但我一般就用直接对应的这个了,不用拟合差不了多少。拟合需要个小程序比较麻烦
8楼2014-05-28 16:09:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

黑暗游侠

铜虫 (正式写手)

引用回帖:
8楼: Originally posted by May.Be at 2014-05-28 16:09:22
你就换一下这个一行的数值就行for i in  6.76 6.78 6.80 6.82 6.84
这就是你估计的c,算出来后找最低能量对应的数值,其实算出来以后有一个拟合的,但我一般就用直接对应的这个了,不用拟合差不了多少。拟合需要个 ...

这个只是对每个c在原有的POSCAR上进行静态计算吗?用不用对每个c都进行一次晶格弛豫,然后再静态计算?
9楼2014-05-28 16:41:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

黑暗游侠

铜虫 (正式写手)

引用回帖:
7楼: Originally posted by dengc09 at 2014-05-28 16:05:40
我写的是我对你的问题的想法,X是你自己的POSCAR的其他部分啊,楼主难道没用过VASP?楼上的run我看了,应该是对的,你再设置一下ISIF=2,把路径改一下,其余文件不变,就可以符合你的要求...

非常感谢
10楼2014-05-28 16:41:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 黑暗游侠 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 明天应该可查了!? +5 chengyan1220 2026-08-23 5/250 2026-08-24 23:28 by 我4大白菜
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +7 scalable 2026-08-24 7/350 2026-08-24 23:01 by anata1209
[基金申请] 没有任何消息-是不是就凉了 +7 图啦图啦 2026-08-24 8/400 2026-08-24 22:00 by maomao_da
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 0XLacIJUOj8D 2026-08-24 3/150 2026-08-24 21:59 by BZKMTicpDhFj
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 0XLacIJUOj8D 2026-08-24 5/250 2026-08-24 21:51 by BZKMTicpDhFj
[基金申请] 人气不行了 +10 fansofjerry 2026-08-21 10/500 2026-08-24 21:03 by zhanghaozhu
[基金申请] 我面上完蛋了 +9 且听虎啸 2026-08-20 10/500 2026-08-24 20:08 by 苦寒来香
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 19/950 2026-08-24 16:52 by iaeyuan
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 什么时候开奖? +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见