24小时热门版块排行榜    

查看: 3798  |  回复: 9

gynan2008

金虫 (正式写手)

[求助] 钛酸钡晶格参数哪个正确?

文献里面有看到立方相晶格参数是3.95*****,前几天用abinit计算时,别人给的常数是7.901******,我用不同的acell计算总能,也是七点几左右。
两个问题,1哪个正确?
                  2、我用不同acel去计算总能,是不是就是在优化晶格参数?
                 3.nband的设置我还是不明白,不是说价电子数/2么,可是设置的要比物质的价电子数多很多啊。
回复此楼

» 猜你喜欢

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

苔花如米小,也学牡丹开!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gynan2008

金虫 (正式写手)

因为第三个问题已经问过了,没太明白,所以再请教一次
苔花如米小,也学牡丹开!
2楼2012-10-31 19:20:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qqxiwuzhiren

铁杆木虫 (著名写手)

【答案】应助回帖

★
感谢参与,应助指数 +1
sunyang1988: 金币+1, 谢谢交流 2012-11-01 19:16:51
3.9多,一般的铁电材料的晶格参数都是3.9左右
3楼2012-10-31 20:13:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gynan2008

金虫 (正式写手)

#Definition of the unit cell  
ndtset 51
acell: 3.5  3.5  3.5
acell+ 0.1 0.1 0.1
#Definition of the atom types
ntypat  3                #定义原胞中原子的类别的数目,本例子中只有1类原子
znucl  56  22  8               #定义原胞中原子的核电荷数
natom  5              #定义原胞中原子的总个数,本例子中有2个原子
typat  1  2  3  3  3            #定义每类原子的个数,本例子中第一类原子有2个


                   #下面定义了原胞中原子的坐标
             xred    0.0000000000E+00  0.0000000000E+00  0.0000000000E+00
                     5.0000000000E-01  5.0000000000E-01  5.0000000000E-01
                     5.0000000000E-01  5.0000000000E-01  0.0000000000E+00
                     5.0000000000E-01  0.0000000000E+00  5.0000000000E-01
                     0.0000000000E+00  5.0000000000E-01  5.0000000000E-01
#Gives the number of band, explicitely (do not take the default)
nband  35        #定义了要计算的能带的数目,最好按这样来设置:
                         # nband > 原胞中总的价电子数目/2 + 10

#Exchange-correlation functional
ixc  1               #定义交换关联函数,本例子中,采用的是Teter Pade 参数化的LDA形式

#Definition of the planewave basis set
ecut    10.0      #定义了平面波的切断动能      

#Definition of the k-point grid
ngkpt  8  8  8            #下面定义了k点网格取样的大小
nshiftk  1        
               
shiftk  0.5  0.5  0.5
      
#Definition of the SCF procedure
nstep   30                   #电子自洽迭代的最大步数
diemac 12.0               #介电常数设置
toldfe  1.0d-10         #电子自洽收敛的标准
enunit   2
rprim:  1   0   0      #下面三行定义了原胞的基矢,本例子中Si是fcc结构
         0   1    0   
         0   0    1

我算了一下,结果是下面的样子,参数设置有问题么?

没办法传图片了,是一条收敛的曲线
苔花如米小,也学牡丹开!
4楼2012-10-31 20:32:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Toapollo

铁杆木虫 (正式写手)

西部牛仔

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
gynan2008: 金币+8, ★★★很有帮助 2012-11-01 10:43:54
sunyang1988: 金币+2, 谢谢交流 2012-11-01 19:17:01
3.9*和7.9*都对的,只是单位不同而已。一个是埃,一个是波尔。
用不同的acel计算总能,是晶胞优化的一种方法,不过效率太低而且太复杂了。但是,你这个代码的功能是体积膨胀,不是结构优化!
nband的设置按章规则来就行了,不懂楼主“设置的要比物质的价电子数多很多”啥意思?对于绝缘体,完全没必要加10个导带能级,那是心理安慰而已。
建议楼主先好好看一下帮助文档~
偷得一日闲便是福,枉作千年计岂非愚!
5楼2012-11-01 01:50:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tangosnow

木虫 (小有名气)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
gynan2008: 金币+2, ★★★很有帮助 2012-11-01 10:44:07
sunyang1988: 金币+1, 谢谢交流 2012-11-01 19:17:11
楼主可以用optcell,ionmov两个关键字去进行优化。
一般情况下,包括你现在的情况,可以不考虑nband,但是在某些特定情况下,需要进行nband的收敛测试
6楼2012-11-01 05:29:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gynan2008

金虫 (正式写手)

引用回帖:
6楼: Originally posted by tangosnow at 2012-11-01 05:29:59
楼主可以用optcell,ionmov两个关键字去进行优化。
一般情况下,包括你现在的情况,可以不考虑nband,但是在某些特定情况下,需要进行nband的收敛测试

哦,谢谢,我再看看教程边学边算吧
苔花如米小,也学牡丹开!
7楼2012-11-01 10:43:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gynan2008

金虫 (正式写手)

引用回帖:
5楼: Originally posted by Toapollo at 2012-11-01 01:50:43
3.9*和7.9*都对的,只是单位不同而已。一个是埃,一个是波尔。
用不同的acel计算总能,是晶胞优化的一种方法,不过效率太低而且太复杂了。但是,你这个代码的功能是体积膨胀,不是结构优化!
nband的设置按章规则 ...

是不是还得学学量化?
苔花如米小,也学牡丹开!
8楼2012-11-01 10:47:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Toapollo

铁杆木虫 (正式写手)

西部牛仔

引用回帖:
8楼: Originally posted by gynan2008 at 2012-11-01 10:47:46
是不是还得学学量化?...

这要看楼主的目标了。
如果计算是为了应用,那么应该是需要啥学啥,这是大多数人的情况;
如果是想从事物性研究,可能就需要相当好的理论物理和化学的功底了,否则是很难弄懂那么多复杂公式所表达的物理化学意义的,也就很难从本质上预测新结构的新机理了~
偷得一日闲便是福,枉作千年计岂非愚!
9楼2012-11-02 01:30:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ttyang

银虫 (初入文坛)

引用回帖:
7楼: Originally posted by gynan2008 at 2012-11-01 10:43:28
哦,谢谢,我再看看教程边学边算吧...

楼主可以给我发一下教程链接吗
10楼2018-03-16 15:56:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gynan2008 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:52 by Z9YWQ5EAO3qp
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:45 by Z9YWQ5EAO3qp
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:33 by Z9YWQ5EAO3qp
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 6/300 2026-09-30 16:14 by Z9YWQ5EAO3qp
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 13:35 by C3vABopYbZJS
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意! +10 zju2000 2026-09-26 13/650 2026-09-30 10:26 by yexuqing
[有机交流] 同一个分子,一条来自文献,一条来自AI——不告诉你答案,你会选哪条? +4 tianxiaoxian 2026-09-28 7/350 2026-09-30 08:34 by 88817753
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-29 19:27 by cqwQDCxMcL3I
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:37 by etmYJ6d2rquH
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:46 by etmYJ6d2rquH
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 GDBe8tDZqE8z 2026-09-28 6/300 2026-09-29 14:41 by etmYJ6d2rquH
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:32 by ez6fGg9abYaj
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 18:51 by 9lS3ad5oOymn
信息提示
请填处理意见