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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 请问siesta计算BornCharge时为什么输出的张量不满足对称性要求呢？谢谢
161/2返回列表12下一页
查看: 2193  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求材料调剂 +7 隔壁陈先生 2026-03-12 7/350 2026-03-18 16:14 by 枫桥ZL
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 0703化学336分求调剂 +6 zbzihdhd 2026-03-15 7/350 2026-03-18 09:53 by zhukairuo
[考研] 277调剂 +5 自由煎饼果子 2026-03-16 6/300 2026-03-17 19:26 by 李leezz
[考研] 考研化学学硕调剂，一志愿985 +4 张vvvv 2026-03-15 6/300 2026-03-17 17:15 by ruiyingmiao
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[论文投稿] 有没有大佬发小论文能带我个二作 +3 增锐漏人 2026-03-17 4/200 2026-03-17 09:26 by xs74101122
[考研] 一志愿，福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[考研] 283求调剂 +10 小楼。 2026-03-12 14/700 2026-03-16 16:08 by 13811244083
[考研] 285求调剂 +6 ytter 2026-03-12 6/300 2026-03-16 15:05 by njzyff
[考研] 0703化学调剂 290分有科研经历，论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 材料工程调剂 +4 咪咪空空 2026-03-11 4/200 2026-03-13 19:57 by JourneyLucky
[硕博家园] 085600 260分求调剂 +3 天空还下雨么 2026-03-13 5/250 2026-03-13 18:46 by 天空还下雨么
[考研] 求b区学校调剂 +3 周56 2026-03-11 3/150 2026-03-13 16:20 by JourneyLucky
[考研] 295求调剂 +3 小匕仔汁 2026-03-12 3/150 2026-03-13 15:17 by vgtyfty
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂！ +3 不想理你 2026-03-12 3/150 2026-03-13 14:18 by JourneyLucky
信息提示
关闭
请填处理意见
关闭