24小时热门版块排行榜    

Znn3bq.jpeg
查看: 3393  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 268分085602化学工程调剂 +12 月照花林。 2026-04-09 12/600 2026-04-09 14:15 by sxx1988
[考研] 化学工程调剂289 +36 yang婷 2026-04-07 37/1850 2026-04-09 14:10 by zxchoshi
[硕博家园] 有没有学校材料专业收跨调(一志愿085410) +4 momo(上岸版) 2026-04-06 7/350 2026-04-09 13:50 by 叫我丢丢
[考研] 288求调剂,一志愿华南理工大学071005 +14 ioodiiij 2026-04-08 14/700 2026-04-09 10:50 by 三七七想上岸
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +20 努力奋斗112 2026-04-07 21/1050 2026-04-08 14:54 by screening
[考研] 336求调剂,一志愿中科大 +9 墨彧 yuyu 2026-04-06 9/450 2026-04-08 11:24 by 想读书的菌菌
[考研] 材料专硕 调剂 +17 CXN123456 2026-04-03 17/850 2026-04-08 11:10 by lijunpoly
[考研] 301求调剂 +10 细胞相关蛋白 2026-04-03 10/500 2026-04-08 10:36 by tjzhao
[考研] 求调剂 +28 111623 2026-04-04 33/1650 2026-04-08 09:24 by 泽润东方
[考研] 336材料与化工085600求调剂 +19 水星记infp 2026-04-05 22/1100 2026-04-07 21:11 by yongzhesheng
[考研] 生物调剂 +5 橙子橙子橙子啊 2026-04-05 9/450 2026-04-07 15:31 by 上岸快快
[考研] 285求调剂 +8 AZMK 2026-04-04 11/550 2026-04-06 13:56 by BruceLiu320
[考研] 324求调剂 +3 k可乐 2026-04-05 4/200 2026-04-06 09:54 by 蓝云思雨
[考研] 一志愿 江南大学 085602 化工专硕 338分求调剂 +15 路痴小琪 2026-04-05 15/750 2026-04-06 09:27 by cql1109
[考研] 302分 085601求调剂推荐 +11 zyx上岸! 2026-04-05 11/550 2026-04-05 22:13 by dongzh2009
[考研] 材料调剂 +7 dxy调剂 2026-04-04 7/350 2026-04-05 09:15 by 陌秋26
[考研] 277求调剂 +4 12A3 2026-04-02 5/250 2026-04-04 20:28 by 蓝云思雨
[考研] 342求调剂 +3 Liang7111 2026-04-04 5/250 2026-04-04 19:47 by dongzh2009
[考研] 一志愿中国石油大学化学工程323分求调剂 +4 化工专硕323分 2026-04-03 6/300 2026-04-03 22:12 by dongzh2009
[考研] 312求调剂 +4 赊月色 2026-04-02 5/250 2026-04-03 08:21 by fangshan711
信息提示
请填处理意见