24小时热门版块排行榜    

Znn3bq.jpeg
查看: 3397  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

gynan2008

金虫 (正式写手)

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

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

» 猜你喜欢

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

苔花如米小,也学牡丹开!
已阅   回复此楼   关注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的回帖
查看全部 10 个回答

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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085400 328分 求调剂 +8 喂你一个大橙子 2026-04-09 12/600 2026-04-09 19:20 by zl8213662
[考研] 一志愿0703化学招61最终排名62化学求调剂 +14 招61排名62 2026-04-07 17/850 2026-04-09 18:25 by hy861222
[考研] 材料调剂 +10 18815505510 2026-04-09 11/550 2026-04-09 17:07 by 544594351
[考研] 327求调剂 +10 Xxjc1107. 2026-04-06 11/550 2026-04-09 01:21 by lature00
[考研] 270求调剂 +3 031127 2026-04-06 4/200 2026-04-08 21:00 by 逆水乘风
[考研] 266调剂 +8 daya sun 2026-04-07 9/450 2026-04-08 20:27 by yutian743
[考研] 265求调剂 +19 小木虫085600 2026-04-06 21/1050 2026-04-08 10:38 by 逆水乘风
[考研] 287求调剂 +6 Fnhc 2026-04-07 6/300 2026-04-08 10:05 by xingguangj
[考研] 计算机408|在校多次国家级竞赛获奖|申请调剂 +4 东山大白鹅 2026-04-05 4/200 2026-04-08 00:18 by chongya
[考研] 22408 318分求调剂 +4 勤奋的小笼包 2026-04-06 6/300 2026-04-07 15:05 by 纸鹤555
[考研] 333求调剂 +6 合乘杨习夕 2026-04-06 6/300 2026-04-07 09:44 by 猪会飞
[考研] 338求调剂 +4 我想上岸ii 2026-04-05 4/200 2026-04-06 21:04 by 木子君1218
[考研] 287求调剂 +3 通信学硕081000 2026-04-03 4/200 2026-04-06 21:03 by going home
[考研] 085405软件工程301分求调剂,专硕可跨专业,四六级已过 +3 静静想想 2026-04-05 3/150 2026-04-06 15:23 by nepu_uu
[考研] 283求调剂 +5 baiiyu 2026-04-05 6/300 2026-04-05 20:35 by 啵啵啵0119
[考研] 一志愿北交大材料工程总分358求调剂 +6 cs0106 2026-04-05 6/300 2026-04-05 16:34 by imissbao
[考研] 341求调剂 +3 洛多罗 2026-04-02 4/200 2026-04-04 21:36 by 智能智慧
[考研] 一志愿南农090401,268,求调剂 +5 一木鸟然 2026-04-04 5/250 2026-04-04 17:07 by babysonlkd
[考研] 338求调剂 +4 zzz,,r 2026-04-03 4/200 2026-04-03 16:39 by lijunpoly
[考研] 专硕 351 086100 也是考的材科基 本科也是材料 +8 202451007219 2026-04-02 8/400 2026-04-03 09:50 by 蓝云思雨
信息提示
请填处理意见