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

登录注册

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

24小时热门版块排行榜

北京石油化工学院2026年研究生招生接收调剂公告

返回列表

Gina88

木虫 (正式写手)

1ST强帖: 1
应助: 9 (幼儿园)
金币: 3328.8
散金: 124
红花: 8
帖子: 395
在线: 113小时
虫号: 1047684
注册: 2010-06-26
专业: 理论和计算化学

[求助] 总能---长度测试不是抛物线形状怎么办呢？

下面贴出了计算Si纳米线的INCAR和产生POSCAR的文件。
计算得到了一个能量-长度图。我有如下问题：
1.
按照网上所说，这个图应该是抛物线的，这个二次曲线系数则和杨氏模量相关。
可是我计算出来不是抛物线的，这是怎么回事啊？
2.
还有，测试的时候，我的长度方向的伸长和压缩量应该选多少啊。
我选的是15%。是不是这个百分比太大，所以导致不是抛物线形的啊？
3.
我这里将ISTART设置成了1，不知道这个设置有没有问题？
看手册上说，要是开始计算时没有WAVECAR，计算时，程序会自动将ISTART当成0啊，后面的计算有WAVECAR，则ISTART=1.

INCAR内容：
SYSTEM=Si111_relax
ISTART=1;ICHARG=2

PREC=Accurate
ENCUT=300
GGA=91;VOSKOWN=1
ISMEAR=0;SIGMA=0.08
EDIFF=1E-4
EDIFFG=-0.08
NELMIN=6

LREAL=A
ALGO=Fast
ISYM=2

ISIF=2
NSW=200
IBRION=2
POTIM=0.5

!LWAVE=.FALSE.
LCHARG=.FALSE.

产生POSCAR的文件：
rm WAVECAR >& /dev/null
for i in 0.85 0.88 0.91 0.94 0.97 1.00 1.03 1.06 1.09 1.12 1.15
do
c=`echo "9.405986241642619*$i" |bc`
cat >POSCAR< Si1112_1cell
1
22.680000000000000  -0.000000000000000  -0.000000000000000
-11.340000000000000  19.641456157831101  0.000000000000000
-0.000000000000019  0.000000000000010  $c
38
Direct
0.387122727907410  0.274245455814824  0.546291382498555
0.330684091861117  0.330684091861119  0.212948761497599
0.387122727907409  0.274245455814822  0.796298348249276
0.330684091861118  0.330684091861118  0.462955727248317
0.499999999999997  0.330684091861118  0.212948761497597
0.499999999999998  0.330684091861117  0.462955727248317
0.274245455814827  0.387122727907411  0.879634003499515
0.387122727907413  0.443561363953705  0.546291382498562
0.330684091861121  0.499999999999999  0.212948761497605
0.274245455814826  0.387122727907412  0.129613106247365
0.387122727907413  0.443561363953704  0.796298348249276
0.330684091861122  0.499999999999999  0.462955727248324
0.443561363953704  0.387122727907410  0.879634003499515
0.556438636046290  0.443561363953704  0.546291382498556
0.500000000000000  0.500000000000000  0.212948761497604
0.443561363953704  0.387122727907410  0.129613106247357
0.556438636046290  0.443561363953703  0.796298348249276
0.500000000000001  0.500000000000000  0.462955727248324
0.612877272092583  0.387122727907411  0.879634003499515
0.669315908138876  0.499999999999998  0.212948761497597
0.612877272092582  0.387122727907411  0.129613106247357
0.669315908138877  0.499999999999998  0.462955727248316
0.387122727907413  0.612877272092585  0.546291382498564
0.387122727907413  0.612877272092583  0.796298348249276
0.443561363953708  0.556438636046291  0.879634003499515
0.556438636046294  0.612877272092585  0.546291382498562
0.500000000000001  0.669315908138879  0.212948761497604
0.443561363953708  0.556438636046292  0.129613106247364
0.556438636046293  0.612877272092584  0.796298348249276
0.500000000000001  0.669315908138880  0.462955727248323
0.612877272092587  0.556438636046290  0.879634003499515
0.725754544185174  0.612877272092586  0.546291382498563
0.669315908138879  0.669315908138879  0.212948761497604
0.612877272092586  0.556438636046292  0.129613106247364
0.725754544185174  0.612877272092585  0.796298348249276
0.669315908138880  0.669315908138879  0.462955727248323
0.612877272092587  0.725754544185170  0.879634003499515
0.612877272092586  0.725754544185172  0.129613106247364
!
echo "c=9.405986241642619*$i angstrom"
time mpirun -np $NP -machinefile $CURDIR/nodelist $HOME/vasp.5.2.11-intel9.1/vasp.5.2.11
en=`grep "entropy=" OUTCAR|tail -1|awk '{printf "%12.6f \n",$7 }'`
E=`grep "TOTEN" OUTCAR|tail -1|awk '{printf "%12.6f \n",$5 }'`
V=`grep "volume" OUTCAR|tail -1|awk '{printf "%12.4f \n",$5 }'`
echo "$c  $en" >>commentVE
cp CONTCAR CONTCAR$i
done

[ Last edited by Gina88 on 2011-11-10 at 15:48 ]

回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

1楼 2011-11-09 15:10:23

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

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子，希望有用哦~

关于弹性常数分析问题，C44特别大可能是什么原因造成的呢已经有3人回复
在药企里面有生产型的企业，外包性质的企业，药品报批的企业，请教哪一个比较有前途呢已经有4人回复
求助这种情况我该怎么办呢已经有5人回复
硕士毕业了，去事业单位好呢还是企业单位好呢已经有25人回复
【求助】做聚氨酯时，用来滴加异氰酸酯的恒压漏斗活塞拔不出来了，该怎么办啊？？已经有12人回复
【交流】做洗涤产品为什么一般是阴离子和非离子复配，阳离子主要是干什么的呢已经有10人回复
急问：点“检查保护”后让输密码要输吗？输了别人怎么打开看呢已经有4人回复
【求助/交流】PBS缓冲溶液高温灭菌后怎么变浑浊了呢？已经有19人回复
【原创】发现教授论文有问题怎么办？已经有11人回复
搞催化的在有机反应路线上创新投什么杂志好呢?Communiations之类的已经有5人回复
请问文章并列第一作者是不是很不好？已经有55人回复
【请教】培养晶体得到的是无机盐，该怎么办呢？已经有9人回复
《寻求帮助》在985高校工作不到半个月的我，要不要开始异地恋呢？~~ 已经有32人回复
急！！土壤中的硝酸盐的来源都有那些个呢？？？尿素施入后怎么变成硝酸盐的呢？？已经有4人回复

点击这里搜索更多相关资源

科研从小木虫开始，人人为我，我为人人

相关版块跳转我要订阅楼主 Gina88 的主题更新