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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 请问siesta计算BornCharge时为什么输出的张量不满足对称性要求呢？谢谢
161/2返回列表12下一页
查看: 2196  |  回复: 15
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

identation

金虫 (著名写手)

[求助] 请问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
这个我觉得应该是没三行属于相同的原子的，因为有效电荷是二阶张量
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-08-02 18:59:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

巴蒂LBJ

木虫 (知名作家)
帮忙顶贴
回复此楼
2楼2011-08-02 19:18:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
对称性不太严格，但是刨除那些小的数值误差，Born有效电荷张量基本上是对的
赞一下 回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
3楼2011-08-02 21:00:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

identation

金虫 (著名写手)
引用回帖:
3楼: Originally posted by iamikaruk at 2011-08-02 21:00:38:
对称性不太严格，但是刨除那些小的数值误差，Born有效电荷张量基本上是对的

我这个体系如果保留小数点后三位，对于Ba和Ti是对的 但是对于O原子 就有两个分量了 这是怎么回事呢
赞一下 回复此楼
4楼2011-08-02 22:34:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
引用回帖:
4楼: Originally posted by identation at 2011-08-02 22:34:07:
我这个体系如果保留小数点后三位，对于Ba和Ti是对的 但是对于O原子 就有两个分量了 这是怎么回事呢

很正常的，O和Ti有两种关系，一种是平行于Ti-O键，一种是垂直于Ti-O键，所以有两个分量
赞一下 回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
5楼2011-08-02 23:43:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

identation

金虫 (著名写手)
引用回帖:
5楼: Originally posted by iamikaruk at 2011-08-02 23:43:00:
很正常的，O和Ti有两种关系，一种是平行于Ti-O键，一种是垂直于Ti-O键，所以有两个分量

你怎么不应助呢 没法给你金币啊 嘿嘿
赞一下 回复此楼
6楼2011-08-03 00:06:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

identation

金虫 (著名写手)
引用回帖:
6楼: Originally posted by identation at 2011-08-03 00:06:09:
你怎么不应助呢 没法给你金币啊 嘿嘿

当前回帖不是应助回帖，您不能对此进行金币奖励哦
赞一下 回复此楼
7楼2011-08-03 00:06:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)

【答案】应助回帖

identation(金币+15): 2011-08-03 00:59:48
ok，好了。填应助了。话说除了abinit外，难道就没有一个第一性原理软件可以把对称性处理的很好吗？
赞一下 回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
8楼2011-08-03 00:47:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

identation

金虫 (著名写手)
引用回帖:
8楼: Originally posted by iamikaruk at 2011-08-03 00:47:57:
ok，好了。填应助了。话说除了abinit外，难道就没有一个第一性原理软件可以把对称性处理的很好吗？

嘿嘿 siesta好像不能识别对称性的 pwscf简化计算量的时候能处理 晶体结构生成还是不行
赞一下 回复此楼
9楼2011-08-03 01:00:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
vasp的Born effective charge对称性也不好，本来三个张量分量根据对称性应该一样的，它硬是弄成不一样了。
赞一下 回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
10楼2011-08-03 09:40:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 identation 的主题更新
161/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +7 步川酷紫123 2026-03-13 7/350 2026-03-18 17:12 by 尽舜尧1
[考研] 085601专硕，总分342求调剂，地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 能源材料化学课题组招收硕士研究生8-10名 +4 脱颖而出 2026-03-16 9/450 2026-03-18 14:13 by 脱颖而出
[考研] 0703化学调剂 +4 pupcoco 2026-03-17 7/350 2026-03-18 12:14 by djl2006
[考研] 265求调剂 +3 梁梁校校 2026-03-17 3/150 2026-03-18 09:12 by zhukairuo
[考研] 268求调剂 +6 简单点0 2026-03-17 6/300 2026-03-18 09:04 by 无际的草原
[考研] 【0856】化学工程（085602）313 分，本科学科评估A类院校化学工程与工艺，诚求调剂 +7 小刘快快上岸 2026-03-11 8/400 2026-03-17 16:57 by ruiyingmiao
[考研] 有没有道铁/土木的想调剂南林，给自己招师弟中～ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[考研] 290求调剂 +6 孔志浩 2026-03-12 11/550 2026-03-17 14:41 by 周舟舟77
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
[考研] 一志愿，福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[基金申请] 国自科面上基金字体 +6 iwuli 2026-03-12 7/350 2026-03-16 21:18 by sculhf
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[基金申请] 现在如何回避去年的某一个专家，不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 255求调剂 +3 李嘉慧， 2026-03-12 4/200 2026-03-14 16:58 by 有只狸奴
[考研] 304求调剂 +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by 星空星月
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] ［0860］321分求调剂，ab区皆可 +4 宝贵热 2026-03-13 4/200 2026-03-13 22:01 by 星空星月
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生，2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
信息提示
关闭
请填处理意见
关闭