| 查看: 690 | 回复: 5 | |||
| 当前主题已经存档。 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 405963512 的 5 个金币 | |||
[交流]
【求助】用脚本优化晶格常数，随着常数变大，能量一直变低？！
|
|||
|
请问个问题，我在用脚本优化晶格常数时，随着晶格常数的变大，能量一直降低，文献中优化得到的常数为8.64，我将自己的脚本文件粘在下面，请大家指教。 这是脚本文件 #!bin/sh rm WAVECAR for a in 8.0 8.2 8.4 8.6 8.8 9.0 9.2 9.4 9.6 9.8 10 do cat >POSCAR < SYSTEM=MMX 1 20.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 20.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 $a 12 8 8 1 2 direct 0.2789786990926944 0.3897263371950115 0.4699673057663697 0.2777454151996646 0.3897263416549008 0.2838224217128855 0.2515833168178175 0.4657577422761720 0.3776045659763455 -0.2789786990926944 -0.3897263371950115 0.4699673057663697 -0.2777454151996646 -0.3897263416549008 0.2838224217128855 -0.2515833168178175 -0.4657577422761721 0.3776045659763455 -0.3897263371950115 0.2789786990926944 0.4699673057663697 -0.3897263416549008 0.2777454151996646 0.2838224217128855 -0.4657577422761720 0.2515833168178175 0.3776045659763455 0.3897263371950115 -0.2789786990926944 0.4699673057663697 0.3897263416549008 -0.2777454151996646 0.2838224217128855 0.4657577422761721 -0.2515833168178175 0.3776045659763455 0.3590474978089331 0.4466374993324282 0.3747565341853553 0.2873549982905391 0.4213825013488531 0.3766565395259330 -0.3590474978089331 -0.4466374993324282 0.3747565341853553 -0.2873549982905391 -0.4213825013488531 0.3766565395259330 -0.4466374993324282 0.3590474978089331 0.3747565341853553 -0.4213825013488531 0.2873549982905391 0.3766565395259330 0.4466374993324281 -0.3590474978089331 0.3747565341853553 0.4213825013488531 -0.2873549982905391 0.3766565395259330 0.3850824974477290 0.4807125024497511 0.5192265344523849 0.4008712537586691 0.4405299997422842 0.2300715340041570 -0.3850824974477290 -0.4807125024497511 0.5192265344523849 -0.4008712537586691 -0.4405299997422842 0.2300715340041570 -0.4807125024497511 0.3850824974477290 0.5192265344523849 -0.4405299997422841 0.4008712537586691 0.2300715340041570 0.4807125024497511 -0.3850824974477291 0.5192265344523849 0.4405299997422841 -0.4008712537586691 0.2300715340041570 0.5000000000000000 0.5000000000000000 0.8072715116404946 0.5000000000000000 0.5000000000000000 0.5091715408229286 0.5000000000000000 0.5000000000000000 0.2414615345858988 ! echo "a=$a" vasp E=`grep "1 F" OSZICAR | tail -1 | awk '{printf "%12.6f \n", $5 }'` echo $a $E >>comment done 我只想优化c轴。 下面是输出文件 8.0 -143.617860 8.2 -148.206810 8.4 -152.032360 8.6 -155.171910 8.8 -157.694190 9.0 -159.657590 9.2 -161.128170 9.4 -162.164350 9.6 -162.817840 9.8 -163.140250 10 -163.165850 我在原来计算简单的晶格测试时，会得到不错的结果，但是都是该变量在第二行，这次是某一个边的长度，不知是否这样 请大家帮忙看下，感激！ |
回复此楼» 猜你喜欢
一志愿211，化学学硕，310分，本科重点双非，求调剂
已经有9人回复
-
一志愿211，化学学硕，310分，本科重点双非，求调剂
已经有17人回复
-
070300化学学硕311分求调剂
已经有15人回复
-
材料调剂
已经有13人回复
-
材料与化工专硕306分找合适调剂
已经有19人回复
-
化学工程调剂289
已经有19人回复
-
化工学硕 285求调剂
已经有22人回复
-
298求调剂
已经有6人回复
-
368化学求调剂
已经有3人回复
-
326求调剂
已经有6人回复
gavinliu7390
木虫 (著名写手)
叶落鹰飞
- 1ST强帖: 3
- 应助: 51 (初中生)
- 贵宾: 0.481
- 金币: 1564.7
- 散金: 2975
- 红花: 34
- 帖子: 1278
- 在线: 241.4小时
- 虫号: 675997
- 注册: 2008-12-17
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
2楼2009-10-25 12:26:43
3楼2009-10-25 15:15:57
4楼2009-10-25 15:17:59
wangzg9292
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3933
- 散金: 56
- 帖子: 348
- 在线: 89.6小时
- 虫号: 723844
- 注册: 2009-03-16
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
5楼2009-10-25 22:12:37
stractor
金虫 (著名写手)
- 1ST强帖: 6
- 应助: 91 (初中生)
- 金币: 1487.6
- 散金: 666
- 红花: 30
- 沙发: 1
- 帖子: 2224
- 在线: 426.8小时
- 虫号: 394087
- 注册: 2007-06-06
- 专业: 半导体材料
6楼2009-10-30 21:28:34
