24小时热门版块排行榜    

CyRhmU.jpeg
查看: 3202  |  回复: 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 的主题更新
信息提示
请填处理意见