24小时热门版块排行榜    

查看: 3880  |  回复: 20

xiejiancqu

金虫 (初入文坛)

[求助] 用CASTEP计算五氧化二钒的态密度图与文献用VASP计算得到的态密度不一样??? 已有2人参与

请教 :我用CASTEP计算的V2O5的态密度,与文献相差很多-就是峰没有那么尖锐,没那么多峰,是我自己在设置参数出错了吗?还是其它原因?
文献用的是VASP ,GGA-PBE ,K为4X6X6,ECUT=500eV设置参数直接用文献的优化结果晶格常数误差很大,我GGA-PBE进行收敛测试得到K点为4x4x4,和ECUT为1270eV,优化后晶格常数c值误差有将近10%,后面用改变泛函为GGA-PBESOL.优化后晶格常数与实验值比较接近。进行计算的结果图所示:

用CASTEP计算五氧化二钒的态密度图与文献用VASP计算得到的态密度不一样???
V2O5-CASTEP.png


用CASTEP计算五氧化二钒的态密度图与文献用VASP计算得到的态密度不一样???-1
V2O5-VASP.png
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

1、DFT_D好像是在在6.1版本以后才能够在CASTEP使用,如果手头没买这个版本以上的只好不考虑;
2、峰的高低取决于计算过程被计入的价电子数目,如果原子数目一样,那么最有可能的就是赝势的差异。赝势方法并不计算所有电子,而VASP的赝势和CASTEP应该是有差异的,PAW方法更是另外一回事。
3、只要价带峰的导带底与文献接近即可。至于导带的形状、大小是否一样,完全可以不用理会。
4、U对于某些特定体系可以是带隙劈裂更大,不过与你现在的问题无关,如果文献加上了U,原则上作比较的话也应该考虑。这次的结果差不多其实已经接近文献了,接下来还可以做些其他的改善。
不一定挂在论坛,计算问题问题欢迎留言。
12楼2015-01-17 17:59:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
13楼: Originally posted by xiejiancqu at 2015-01-17 19:55:51
非常谢谢您的耐心指点!我用的是MS.6.0,文献用的VASP,赝势用PAW方法。我用CASTEP采用的是ultrasoft。呵呵......那我接下来应该如何加以改善呢?比较困惑,不好意思(是只新虫),还需要您再指点指点!万分感谢 ...

KalaShayminS前辈说的有些道理,我当时也没注意,右侧少个峰有可能是“能级不够多”,可以计算能带结构的时候加些空带数目,默认12肯定不够,到了6eV是电子态直接“掉下去”了。

峰值高低除了赝势或PAW的价电子组态不同,也有可能是你采用PBEsol和PBE不同造成的,目前尚不确定,一定要找出原因可以对照两者价电子组态看看有多大差异。展宽很大的时候电子态弥散开来,所以大的展宽的DOS不怎么好看。

我忘记了6.0能不能做DFTD,如果不能也就只有这样了,能+D的话对于VDW结合的方向的晶格常数会有些改善。可以加上个U看看会不会改善能带的结果,U值怎么取合理就得再想办法了。
不一定挂在论坛,计算问题问题欢迎留言。
15楼2015-01-18 08:15:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

tsglss

铁杆木虫 (著名写手)

其实已经很好了

[ 发自手机版 http://muchong.com/3g ]
做好自己
2楼2015-01-15 22:57:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiejiancqu

金虫 (初入文坛)

引用回帖:
2楼: Originally posted by tsglss at 2015-01-15 22:57:41
其实已经很好了

谢谢您的回复!我就是在疑惑是不是我的计算有误导致数据不对,所以峰没那么多,也没那么尖锐,还是因为所以软件不同造成的...... 望具体指点,谢谢!
3楼2015-01-16 09:27:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

KalaShayminS

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
ljw4010: 金币+2, 谢谢交流! 2015-01-18 21:05:04
xiejiancqu: 金币+4, ★★★很有帮助 2015-04-26 09:57:55
从图上对比来看应该是绘图时的设置问题,DOS绘图需要展宽,展宽越大,每个峰就越不尖锐,叠加后峰越少。
4楼2015-01-16 10:30:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiejiancqu

金虫 (初入文坛)

引用回帖:
4楼: Originally posted by KalaShayminS at 2015-01-16 10:30:28
从图上对比来看应该是绘图时的设置问题,DOS绘图需要展宽,展宽越大,每个峰就越不尖锐,叠加后峰越少。

谢谢指点!可文献的也是展宽的,我用Origin画图的横坐标设置与文献一样,如果是绘图设置的问题那我又应该怎样设置呢?(呵呵 Origin不是很熟悉使用)再问一个比较二的问题,我这个图要是写进文章里能用么?谢谢!
5楼2015-01-16 12:01:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
5楼: Originally posted by xiejiancqu at 2015-01-16 12:01:50
谢谢指点!可文献的也是展宽的,我用Origin画图的横坐标设置与文献一样,如果是绘图设置的问题那我又应该怎样设置呢?(呵呵 Origin不是很熟悉使用)再问一个比较二的问题,我这个图要是写进文章里能用么?谢谢 ...

不是绘图软件(数据处理软件Origin照成的)看一下计算使用的smearing和性质中的smearing是不是都太大了,这个smearing就是展宽。计算的时候smearing过大可能也会造成结果不可靠,该取多少最好由测试决定。

晶格常数不同有可能是赝势导致的,或是其他一些细节引起的,这个问题暂时不得而知。你做的是层状结构吧,也许考虑GrimmeD有可能有好的结果。

» 本帖已获得的红花(最新10朵)

不一定挂在论坛,计算问题问题欢迎留言。
6楼2015-01-16 23:51:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiejiancqu

金虫 (初入文坛)

送红花一朵
引用回帖:
6楼: Originally posted by 卡开发发 at 2015-01-16 23:51:09
不是绘图软件(数据处理软件Origin照成的)看一下计算使用的smearing和性质中的smearing是不是都太大了,这个smearing就是展宽。计算的时候smearing过大可能也会造成结果不可靠,该取多少最好由测试决定。

晶格 ...

谢谢指导!我调一下smearing试试,我做的是层状结构,您说的GrimmeD具体是什么呢(未了解)?麻烦您再具体指点一下,万分感谢!
7楼2015-01-17 09:50:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
7楼: Originally posted by xiejiancqu at 2015-01-17 09:50:41
谢谢指导!我调一下smearing试试,我做的是层状结构,您说的GrimmeD具体是什么呢(未了解)?麻烦您再具体指点一下,万分感谢!...

DFT-D那一项,弱作用DFT(LDA/GGA)的描述比较差。
不一定挂在论坛,计算问题问题欢迎留言。
8楼2015-01-17 12:16:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hhbingf

木虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
ljw4010: 金币+2, 谢谢交流! 2015-01-18 21:05:21
1,尖峰应该是因为你们用的晶胞不同,文献中所用的晶胞原子应该比你的所用的多;2,你所算出来的dos没有带隙,可能是因为你没有加U。
9楼2015-01-17 12:41:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiejiancqu

金虫 (初入文坛)

引用回帖:
8楼: Originally posted by 卡开发发 at 2015-01-17 12:16:36
DFT-D那一项,弱作用DFT(LDA/GGA)的描述比较差。...

谢谢您的指点!DFT-D那一项,用CASTEP具体如何考虑呢?不是很明白,麻烦您再指点一下,万分感谢!
呵呵,我将smearing调小,峰确实要多些尖锐些,可是与文献价带的最高峰比导带的最高峰还要高,而且在费米能级的右侧还少一个峰(见下图),这是又什么原因所致呢?麻烦您再指点,谢谢!
用CASTEP计算五氧化二钒的态密度图与文献用VASP计算得到的态密度不一样???-2
V2O5-CASTEP.png


用CASTEP计算五氧化二钒的态密度图与文献用VASP计算得到的态密度不一样???-3
V2O5-VASP.png

10楼2015-01-17 16:37:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiejiancqu 的主题更新
信息提示
请填处理意见