应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Abinit » 钛酸钡晶格参数哪个正确？
51/1返回列表
查看: 3201  |  回复: 9
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

gynan2008

金虫 (正式写手)

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

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
苔花如米小，也学牡丹开！
1楼2012-10-31 19:18:58
已阅   回复此楼   关注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的回帖
查看全部 10 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭