24小时热门版块排行榜    

Znn3bq.jpeg
查看: 3424  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 337求调剂 +3 jyz04 2026-04-18 3/150 2026-04-20 12:24 by 研可安
[教师之家] 又一批高校组建人工智能学院 师资行吗 不是骗人吗 +3 yexuqing 2026-04-19 3/150 2026-04-20 11:29 by jurkat.1640
[考博] 湖南大学刘巧玲课题组2026年第二批次博士研究生招生信息 +3 南风观火 2026-04-18 5/250 2026-04-20 10:13 by 南风观火
[考研] 通信工程求调剂!!! +7 zlb770521 2026-04-14 7/350 2026-04-19 20:56 by Equinoxhua
[考研] 297,工科调剂? +11 河南农业大学-能 2026-04-14 11/550 2026-04-19 20:07 by Equinoxhua
[考研] 291求调剂 +11 关忆北. 2026-04-14 11/550 2026-04-19 17:16 by 中豫男
[考研] 304求调剂 +8 castLight 2026-04-16 8/400 2026-04-19 17:14 by 中豫男
[考研] 求调剂推荐 +9 小聂爱学习 2026-04-14 9/450 2026-04-19 17:03 by 中豫男
[考研] 求调剂 +10 小聂爱学习 2026-04-16 12/600 2026-04-19 16:51 by 中豫男
[考研] 求调剂 +6 苦命人。。。 2026-04-18 7/350 2026-04-19 16:27 by 中豫男
[考研] 300求调剂 +12 橙a777 2026-04-15 12/600 2026-04-18 23:51 by 路病情
[考研] 297,工科调剂? +5 河南农业大学-能 2026-04-14 5/250 2026-04-18 15:17 by Equinoxhua
[考研] 22408 312求调剂 +24 门路摸摸 2026-04-14 26/1300 2026-04-18 13:04 by wunaiy88
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
[考研] 急需调剂 +9 绝不放弃22 2026-04-15 10/500 2026-04-18 08:09 by chixmc
[考研] 一志愿华中农业071010,320求调剂 +17 困困困困坤坤 2026-04-14 19/950 2026-04-17 20:08 by 关一盏灯cd
[考研] 295分求调剂 +5 ?要上岸? 2026-04-17 5/250 2026-04-17 16:51 by fenglj492
[考研] 求调剂学校 +14 不会吃肉 2026-04-13 16/800 2026-04-15 21:59 by noqvsozv
[考研] 一志愿A区211,22408 321求调剂 +6 随心所欲☆ 2026-04-15 7/350 2026-04-15 21:45 by lbsjt
[考研] 245求调剂 +6 冰糖橘?汽水 2026-04-13 10/500 2026-04-14 10:49 by jyl0317
信息提示
请填处理意见