版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

bluegirlll

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 105.9
帖子: 57
在线: 20.6小时
虫号: 1046082

[交流] 【求助】VASP怎么做晶格参数的结构优化

求助，怎么做晶格参数的优化，有没有例子脚本？非常感谢

回复此楼

» 收录本帖的淘帖专辑推荐

VASP	VASP软件从入门到高深

» 猜你喜欢

上海工程技术大学【激光智能制造】课题组招收硕士已经有4人回复
想换工作。大多数高校都是评职称时认可5年内在原单位取得的成果吗？已经有5人回复
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？已经有6人回复
带资进组求博导收留已经有10人回复
最近几年招的学生写论文不引自己组发的文章已经有11人回复
需要合成515-64-0，50g，能接单的留言已经有3人回复
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）已经有3人回复
临港实验室与上科大联培博士招生1名已经有8人回复

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

关于VASP一些琐事的求助（K点，吸附优化等）已经有12人回复
晶体结构优化和画能带图的方法是否保持一致？？？已经有12人回复
单轴施加应变结构优化的方法有哪些？已经有7人回复
求助vasp_晶格参数扫描已经有5人回复
vasp计算DOS,以及如何判断参数的设置是合理的，达到要求的精度已经有7人回复
vasp计算中，如何选择正确的K点已经有7人回复
【活动or资源】vasp晶体结构优化已经有142人回复
【求助】晶格常数大了以后 VASP的计算值准确吗已经有5人回复
【求助】vasp优化晶格参数a 已经有13人回复
【求助】简单问题----关于vasp晶格优化已经有7人回复
【求助】侯老师的vasp入门的晶格优化部分没看懂，问一下已经有25人回复
【求助】VASP结构优化已经有7人回复
【求助成功】VASP中，如何固定磁矩做结构优化已经有14人回复
【原创】晶格优化中ISIF取值的探究（六方MgO为例）已经有9人回复
【讨论】vasp结构优化之前的测试已经有11人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2010-12-07 09:44:41

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

bluegirlll

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 105.9
帖子: 57
在线: 20.6小时
虫号: 1046082

我已经按您说的改了，可是出来这样的结果，它去了OSIZCAR文件中的好多数据
SUMMARY.fcc内容如下：

3.5 DAV: 9 0.616371635786E+02 -0.21638E+01 -0.36614E+00 887 0.211E+01 0.379E+00 DAV: 10 0.598839270162E+02 -0.17532E+01 -0.22497E+00 946 0.164E+01 0.135E+00 DAV: 11 0.599095645060E+02 0.25637E-01 -0.14704E-01 931 0.481E+00 0.625E-01 DAV: 12 0.599206379143E+02 0.11073E-01 -0.57937E-02 824 0.307E+00 0.358E-01 DAV: 13 0.599130268782E+02 -0.76110E-02 -0.27145E-02 953 0.152E+00 0.107E-01 DAV: 14 0.599127542034E+02 -0.27267E-03 -0.25634E-03 1209 0.508E-01 0.870E-02 DAV: 15 0.599126009216E+02 -0.15328E-03 -0.17119E-04 615 0.166E-01 0.464E-02 DAV: 16 0.599123843746E+02 -0.21655E-03 -0.15291E-04 685 0.131E-01 0.242E-02 DAV: 17 0.599123887861E+02 0.44115E-05 -0.17822E-05 398 0.428E-02 1 F= 0.59912389E+02 E0= 0.59912389E+02 d E =0.000000E+00 mag= -0.0009
3.6 DAV: 8 0.486136355673E+02 -0.35626E+01 -0.34687E+00 890 0.236E+01 0.522E+00 DAV: 9 0.469129254896E+02 -0.17007E+01 -0.18613E+00 940 0.131E+01 0.343E+00 DAV: 10 0.461778722637E+02 -0.73505E+00 -0.13600E+00 876 0.131E+01 0.144E+00 DAV: 11 0.462037274218E+02 0.25855E-01 -0.10867E-01 882 0.402E+00 0.832E-01 DAV: 12 0.461914507126E+02 -0.12277E-01 -0.69380E-02 829 0.287E+00 0.323E-01 DAV: 13 0.461915837635E+02 0.13305E-03 -0.55215E-03 1206 0.675E-01 0.121E-01 DAV: 14 0.461917513859E+02 0.16762E-03 -0.15565E-03 1084 0.526E-01 0.796E-02 DAV: 15 0.461913353655E+02 -0.41602E-03 -0.36541E-04 788 0.163E-01 0.584E-02 DAV: 16 0.461912802755E+02 -0.55090E-04 -0.28427E-04 663 0.170E-01 1 F= 0.46191280E+02 E0= 0.46191280E+02 d E =0.000000E+00 mag= 0.0122
3.7 DAV: 21 0.345758323297E+02 -0.93737E-03 -0.94467E-04 929 0.306E-01 0.402E-01 DAV: 22 0.345795760961E+02 0.37438E-02 -0.48287E-03 922 0.658E-01 0.264E-01 DAV: 23 0.345802496016E+02 0.67351E-03 -0.12725E-03 968 0.325E-01 0.202E-01 DAV: 24 0.345826284269E+02 0.23788E-02 -0.54723E-04 796 0.218E-01 0.161E-01 DAV: 25 0.345871012547E+02 0.44728E-02 -0.30116E-03 938 0.545E-01 0.130E-01 DAV: 26 0.345872092251E+02 0.10797E-03 -0.30063E-04 728 0.187E-01 0.972E-02 DAV: 27 0.345858613196E+02 -0.13479E-02 -0.55383E-04 816 0.274E-01 0.609E-02 DAV: 28 0.345862324715E+02 0.37115E-03 -0.52166E-04 748 0.350E-01 0.271E-02 DAV: 29 0.345861649687E+02 -0.67503E-04 -0.13802E-05 395 0.461E-02 1 F= 0.34586165E+02 E0= 0.34586165E+02 d E =0.000000E+00 mag= 0.42683.8 DAV: 7 0.279114265739E+02 -0.16376E+01 -0.56337E+00 1022 0.204E+01 0.665E+00 DAV: 8 0.254025058027E+02 -0.25089E+01 -0.34546E+00 870 0.240E+01 0.374E+00 DAV: 9 0.248851128671E+02 -0.51739E+00 -0.47666E-01 992 0.551E+00 0.273E+00 DAV: 10 0.244305976128E+02 -0.45452E+00 -0.56551E-01 893 0.739E+00 0.128E+00 DAV: 11 0.243774407325E+02 -0.53157E-01 -0.99253E-02 1009 0.289E+00 0.572E-01 DAV: 12 0.243763728670E+02 -0.10679E-02 -0.99008E-03 878 0.101E+00 0.177E-01 DAV: 13 0.243764253820E+02 0.52515E-04 -0.16325E-03 1012 0.268E-01 0.420E-02 DAV: 14 0.243762900673E+02 -0.13531E-03 -0.12378E-04 502 0.108E-01 0.291E-02 DAV: 15 0.243762466643E+02 -0.43403E-04 -0.59665E-05 425 0.931E-02 1 F= 0.24376247E+02 E0= 0.24376247E+02 d E =0.000000E+00 mag= 0.00023.9 DAV: 6 0.186939152777E+02 -0.20100E+00 -0.56562E+00 1138 0.238E+01 0.105E+01 DAV: 7 0.175924468941E+02 -0.11015E+01 -0.48486E+00 1076 0.181E+01 0.604E+00 DAV: 8 0.159033008957E+02 -0.16891E+01 -0.34831E+00 842 0.257E+01 0.342E+00 DAV: 9 0.157252883319E+02 -0.17801E+00 -0.39195E-01 1270 0.648E+00 0.242E+00 DAV: 10 0.153946901067E+02 -0.33060E+00 -0.28033E-01 986 0.450E+00 0.127E+00 DAV: 11 0.153219436031E+02 -0.72747E-01 -0.12657E-01 974 0.338E+00 0.493E-01 DAV: 12 0.153217217700E+02 -0.22183E-03 -0.12968E-02 1019 0.155E+00 0.230E-01 DAV: 13 0.153228266672E+02 0.11049E-02 -0.25340E-03 1071 0.660E-01 0.796E-02 DAV: 14 0.153229062451E+02 0.79578E-04 -0.32392E-04 758 0.153E-01 1 F= 0.15322906E+02 E0= 0.15322906E+02 d E =0.000000E+00 mag= 0.00034.0 DAV: 6 0.931139208363E+01 0.56404E+00 -0.56782E+00 1180 0.241E+01 0.957E+00 DAV: 7 0.878331383251E+01 -0.52808E+00 -0.41165E+00 1150 0.170E+01 0.519E+00 DAV: 8 0.775427396386E+01 -0.10290E+01 -0.20414E+00 854 0.190E+01 0.279E+00 DAV: 9 0.755153875548E+01 -0.20274E+00 -0.26699E-01 1247 0.457E+00 0.199E+00 DAV: 10 0.734397362079E+01 -0.20757E+00 -0.13448E-01 918 0.347E+00 0.112E+00 DAV: 11 0.728475733198E+01 -0.59216E-01 -0.81750E-02 1000 0.288E+00 0.398E-01 DAV: 12 0.728444227142E+01 -0.31506E-03 -0.36268E-03 925 0.636E-01 0.201E-01 DAV: 13 0.728455095831E+01 0.10869E-03 -0.14001E-03 1122 0.340E-01 0.352E-02 DAV: 14 0.728458422560E+01 0.33267E-04 -0.81814E-05 452 0.917E-02 1 F= 0.72845842E+01 E0= 0.72845842E+01 d E =0.000000E+00 mag= -0.00014.1 DAV: 5 0.167530005868E+00 -0.14240E-01 -0.14230E-01 988 0.256E+00 0.177E+01 DAV: 6 0.125671077861E+01 0.10892E+01 -0.57949E+00 1187 0.243E+01 0.874E+00 DAV: 7 0.11036......................................................
请问这是怎么回事啊？非常感谢！

引用回帖:

Originally posted by zhang_jaj at 2010-12-08 10:36:06:

把引号改成··，就是tab键上的那个键，试试吧

赞一下

回复此楼

7楼2010-12-08 11:05:12

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

查看全部 18 个回答

first_yg

银虫 (小有名气)

应助: 1 (幼儿园)
金币: 463.5
帖子: 132
在线: 81.8小时
虫号: 374955

★
sunyang1988(金币+1):谢谢指点 2010-12-07 09:55:44
bluegirlll(金币+1): 2010-12-08 16:18:50

侯博士的中文说明里就有，你可以查查看

赞一下(1人)

回复此楼

2楼2010-12-07 09:51:38

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

ustbwgs

木虫 (正式写手)

应助: 14 (小学生)
金币: 1449.6
帖子: 720
在线: 165.5小时
虫号: 1141993

★ ★
zzy870720z(金币+2):谢谢提示 2010-12-07 14:18:47
bluegirlll(金币+1): 2010-12-08 09:27:38

先根据文献或实验参数写个差不多的，然后再在这个取值附近的范围内，写个脚本计算单点能，做能量收敛测试。

赞一下(1人)

回复此楼

3楼2010-12-07 14:11:34

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

bluegirlll

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 105.9
帖子: 57
在线: 20.6小时
虫号: 1046082

谢谢您啊，我用的是这个个脚本：

#! /bin/bash
BIN=/home/jhm/VASP/vasp
for i in 3.5 3.6 3.7 3.8 3.9 4.0 4.1 4.2 4.3 ; do
cat >POSCAR < fcc:
$i
0.0 0.5 0.5
0.5 0.0 0.5
0.5 0.5 0.0
1 2
cartesian
0 0 0
0.25 0.25 0.25
0.75 0.75 0.75
!
echo "a= $i" ; $BIN
E=·tail -l OSZICAR·; echo $i $E >>SUMMARY.fcc
done
计算得到的SUMMARY.fcc文件的内容为：

3.5 tail -l OSZICAR
3.6 tail -l OSZICAR
3.7 tail -l OSZICAR
3.8 tail -l OSZICAR
3.9 tail -l OSZICAR
4.0 tail -l OSZICAR
4.1 tail -l OSZICAR
4.2 tail -l OSZICAR
4.3 tail -l OSZICAR

是我脚本那个地方错了吗？请高人指点

引用回帖:

Originally posted by ustbwgs at 2010-12-07 14:11:34:
先根据文献或实验参数写个差不多的，然后再在这个取值附近的范围内，写个脚本计算单点能，做能量收敛测试。

赞一下

回复此楼

4楼2010-12-08 09:34:25

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

查看全部 18 个回答