| 查看: 2649 | 回复: 15 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
请问siesta计算BornCharge时为什么输出的张量不满足对称性要求呢？谢谢
|
||
|
我测试siesta，计算立方相BaTiO3的波恩有效电荷，结果输出的张量中分量很多啊，立方晶每个原子的有效电荷应该只有一个分量啊，另外，立方BaTiO3的三个O原子应该是等价的啊，为什么计算出的有效电荷不一样呢？ 感激不尽！求指教！ 输入文件如下： SystemName Barium titanate (BaTiO3) cubic structure SystemLabel batio3 NumberOfSpecies 3 NumberOfAtoms 5 %block ChemicalSpeciesLabel 1 56 Ba 2 22 Ti 3 8 O %endblock ChemicalSpeciesLabel %block PS.lmax Ba 3 Ti 3 O 3 %endblock PS.lmax PAO.BasisType split %block PAO.Basis Ba 5 0.70 n=5 0 1 E 98.95 5.90 6.39899920099572 1.00000000000000 n=6 0 2 E 97.95 6.50 6.99956358549850 5.99949655050583 1.00000000000000 1.00000000000000 n=5 1 1 E 96.95 6.09 6.59922580635556 1.00000000000000 n=5 2 1 E 95.95 6.60 7.09995685034315 1.00000000000000 n=6 1 1 E 94.96 6.60 7.09996103032325 1.00000000000000 Ti 5 1.91 n=3 0 1 E 93.95 5.20 5.69946662616249 1.00000000000000 n=3 1 1 E 95.47 5.20 5.69941339465994 1.00000000000000 n=4 0 2 E 96.47 5.60 6.09996398975307 5.09944363262274 1.00000000000000 1.00000000000000 n=3 2 2 E 46.05 4.95 5.94327035784617 4.70009988294302 1.00000000000000 1.00000000000000 n=4 1 1 E 0.50 1.77 3.05365979938936 1.00000000000000 O 3 -0.28 n=2 0 2 E 40.58 3.95 4.95272270428712 3.60331408800389 1.00000000000000 1.00000000000000 n=2 1 2 E 36.78 4.35 4.99990228025066 3.89745395068600 1.00000000000000 1.00000000000000 n=3 2 1 E 21.69 0.93 2.73276990670788 1.00000000000000 %endblock PAO.Basis LatticeConstant 7.57 bohr %block LatticeVectors 1.00 0.00 0.00 0.00 1.00 0.00 0.00 0.00 1.00 %endblock LatticeVectors AtomicCoordinatesFormat Fractional %block AtomicCoordinatesAndAtomicSpecies 0.00 0.00 0.00 1 0.50 0.50 0.50 2 0.00 0.50 0.50 3 0.50 0.00 0.50 3 0.50 0.50 0.00 3 %endblock AtomicCoordinatesAndAtomicSpecies %block kgrid_Monkhorst_Pack 6 0 0 0.5 0 6 0 0.5 0 0 6 0.5 %endblock kgrid_Monkhorst_Pack MeshCutoff 200 Ry DM.NumberPulay 3 DM.UseSaveDM F DM.MixingWeight 0.0200 # New DM amount for next SCF cycle DM.Tolerance 1.d-4 # Tolerance in maximum difference # between input and output DM DM.NumberKick 15 DM.KickMixingWeight 1.000 ElectronicTemperature 100 K MaxSCFIterations 400 MD.TypeOfRun FC MD.FCDispl 0.04 Bohr BornCharge T %block PolarizationGrids 7 3 3 3 7 3 3 3 7 %endblock PolarizationGrids WriteCoorStep .true. WriteForces .true. 计算的有效电荷如下： BC matrix 2.2565879 0.0000003 0.0000000 0.0002761 2.2586398 -0.0000013 -0.0000002 0.0003220 2.2587290 7.0653476 -0.0000006 0.0000283 -0.0000712 7.0646728 -0.0000014 -0.0000002 0.0000143 7.0650252 -5.4552768 -0.0000003 -0.0001544 0.0000997 -1.9375166 0.0000012 -0.0000020 0.0000339 -1.9375039 -1.9375306 -0.0000008 0.0000191 -0.0000036 -5.4553489 -0.0000027 0.0000038 0.0000020 -1.9376244 -1.9376920 0.0000007 -0.0000649 -0.0000915 -1.9377835 -0.0000031 -0.0000028 -0.0000270 -5.4563491 这个我觉得应该是没三行属于相同的原子的，因为有效电荷是二阶张量 |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有94人回复
-
投稿咨询
已经有4人回复
-
时间戳变了，能看出什么问题？
已经有5人回复
-
93BebMhtakh前后11位开头都是大写
已经有5人回复
-
感觉是下周放榜了
已经有11人回复
-
filecode=后面第一个是大写字母
已经有10人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有15人回复
-
今天系统多次维护，明天很可能放榜！
已经有9人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 《张量分析》田宗若.西北工业大学出版社
已经有163人回复
- 什么是张量，基本思想是什么？
已经有14人回复
- 什么是张量？基本的思想是什么呢
已经有13人回复
- 期盼地球物理专业老师同学们进行矩张量反演方面的指点
已经有13人回复
- 通向凝聚态物理和量子化学的终极之路--密度矩阵重整化群，从矩阵乘积态到张量网态
已经有19人回复
- 【求助】用vasp算压电张量（系数），弹性张量需要注意的问题？
已经有9人回复
- 【求助】请问有人用WIEN2K算过光电导张量吗？
已经有4人回复
- 【请教】关于电子顺磁光谱EPR谱图的g张量是如何求得的
已经有5人回复
- 【求助】张量的运算符号（塑性细观力学一书p121）
已经有9人回复
- 【求助】张量运算请教
已经有7人回复
- 【求助】关于发表张量的一篇文章一般都找哪些杂志呀，想投EI类的
已经有3人回复
4楼2011-08-02 22:34:07
巴蒂LBJ
木虫 (知名作家)
- 应助: 5 (幼儿园)
- 金币: 3028.3
- 散金: 1107
- 红花: 12
- 沙发: 38
- 帖子: 7098
- 在线: 321.5小时
- 虫号: 1165239
- 注册: 2010-12-08
- 性别: GG
- 专业: 配位化学
帮忙顶贴 2楼2011-08-02 19:18:31
iamikaruk
木虫 (著名写手)
- 1ST强帖: 1
- 应助: 11 (小学生)
- 金币: 4950.3
- 散金: 118
- 红花: 76
- 帖子: 1362
- 在线: 543.8小时
- 虫号: 1315985
- 注册: 2011-06-05
- 专业: 无机非金属类电介质与电解
3楼2011-08-02 21:00:38
iamikaruk
木虫 (著名写手)
- 1ST强帖: 1
- 应助: 11 (小学生)
- 金币: 4950.3
- 散金: 118
- 红花: 76
- 帖子: 1362
- 在线: 543.8小时
- 虫号: 1315985
- 注册: 2011-06-05
- 专业: 无机非金属类电介质与电解
5楼2011-08-02 23:43:00
