| 查看: 3665 | 回复: 9 | ||
gynan2008金虫 (正式写手)
|
[求助]
钛酸钡晶格参数哪个正确？
|
|
文献里面有看到立方相晶格参数是3.95*****，前几天用abinit计算时，别人给的常数是7.901******，我用不同的acell计算总能，也是七点几左右。 两个问题，1哪个正确？ 2、我用不同acel去计算总能，是不是就是在优化晶格参数？ 3.nband的设置我还是不明白，不是说价电子数/2么，可是设置的要比物质的价电子数多很多啊。 |
回复此楼» 猜你喜欢
让我中一个面上吧！
已经有10人回复
-
科研孤儿太难了
已经有14人回复
-
言之凿凿，明天官网会通知放榜？
已经有17人回复
-
时间戳又变了
已经有17人回复
-
今天系统多次维护，明天很可能放榜！
已经有11人回复
-
时间戳变了，能看出什么问题？
已经有23人回复
-
只有每年这种时候来逛逛小木虫
已经有22人回复
-
放榜前的不淡定
已经有12人回复
-
重要消息，中午系统在维护
已经有12人回复
-
今天基金会出结果吗？20260819
已经有15人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 对XRD进行分析，计算钛酸钡的晶格常数
已经有9人回复
- 晶胞参数影响钛酸钡基陶瓷居里温度的内在原因？
已经有3人回复
- 已知晶格参数，怎么构建一个立方体出来啊
已经有7人回复
- 求助vasp_晶格参数扫描
已经有5人回复
- 做掺杂后超晶胞晶格参数和体积的变化
已经有5人回复
- 求问用X-ray diffraction 能定量的测定出晶体的晶格参数吗？
已经有5人回复
- CaIn2O4的Pnam空间群结构晶格参数和CaBi2O4的C2/c空间群晶格参数
已经有5人回复
- 【求助】在晶体表面外延生长的ZnO薄膜晶格参数
已经有6人回复
- 【求助】VASP怎么做晶格参数的结构优化
已经有17人回复
- 【求助】vasp优化晶格参数a
已经有13人回复
- 【讨论】如何得出晶格参数a,b,c的值
已经有13人回复
- 【求助】如何从Murnaghan方程的结果中得到优化的晶格参数？
已经有11人回复
- 【求助】XRD算晶格参数
已经有9人回复
- 【求助】lammps如何测量模拟以后材料的晶格参数？
已经有7人回复
- 【求助】怎么用jade计算粒径和晶格参数
已经有6人回复
gynan2008
金虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 1621.4
- 散金: 100
- 红花: 4
- 帖子: 367
- 在线: 94.4小时
- 虫号: 1421581
- 注册: 2011-09-28
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
2楼2012-10-31 19:20:16
qqxiwuzhiren
铁杆木虫 (著名写手)
- 应助: 53 (初中生)
- 金币: 8526.2
- 散金: 334
- 红花: 12
- 帖子: 1076
- 在线: 261.6小时
- 虫号: 1466233
- 注册: 2011-10-29
- 性别: GG
- 专业: 功能陶瓷
3楼2012-10-31 20:13:00
gynan2008
金虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 1621.4
- 散金: 100
- 红花: 4
- 帖子: 367
- 在线: 94.4小时
- 虫号: 1421581
- 注册: 2011-09-28
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
|
#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
Toapollo
铁杆木虫 (正式写手)
西部牛仔
- 应助: 18 (小学生)
- 金币: 7428.6
- 散金: 80
- 红花: 7
- 帖子: 603
- 在线: 263.7小时
- 虫号: 499729
- 注册: 2008-02-13
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
gynan2008: 金币+8, ★★★很有帮助 2012-11-01 10:43:54
sunyang1988: 金币+2, 谢谢交流 2012-11-01 19:17:01
感谢参与，应助指数 +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
tangosnow
木虫 (小有名气)
- 应助: 37 (小学生)
- 金币: 4531.9
- 散金: 5
- 红花: 4
- 帖子: 295
- 在线: 284.6小时
- 虫号: 325625
- 注册: 2007-03-17
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
6楼2012-11-01 05:29:59
gynan2008
金虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 1621.4
- 散金: 100
- 红花: 4
- 帖子: 367
- 在线: 94.4小时
- 虫号: 1421581
- 注册: 2011-09-28
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
7楼2012-11-01 10:43:28
gynan2008
金虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 1621.4
- 散金: 100
- 红花: 4
- 帖子: 367
- 在线: 94.4小时
- 虫号: 1421581
- 注册: 2011-09-28
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
是不是还得学学量化？ 8楼2012-11-01 10:47:46
Toapollo
铁杆木虫 (正式写手)
西部牛仔
- 应助: 18 (小学生)
- 金币: 7428.6
- 散金: 80
- 红花: 7
- 帖子: 603
- 在线: 263.7小时
- 虫号: 499729
- 注册: 2008-02-13
- 专业: 凝聚态物性 II ：电子结构
9楼2012-11-02 01:30:46
10楼2018-03-16 15:56:04
