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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 请问siesta计算BornCharge时为什么输出的张量不满足对称性要求呢？谢谢
161/2返回列表12下一页
查看: 2194  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料与化工 +5 安全上岸！ 2026-03-16 5/250 2026-03-18 15:33 by cmz0325
[考研] 286求调剂 +6 lemonzzn 2026-03-16 9/450 2026-03-18 15:31 by vgtyfty
[考研] 085601材料工程专硕求调剂 +6 慕寒mio 2026-03-16 6/300 2026-03-18 14:26 by 007_lilei
[考研] 一志愿西南交大，求调剂 +4 材化逐梦人 2026-03-18 4/200 2026-03-18 14:22 by 007_lilei
[考研] 材料专硕306英一数二 +10 z1z2z3879 2026-03-16 13/650 2026-03-18 14:20 by 007_lilei
[考研] 302求调剂 +10 呼呼呼。。。。 2026-03-17 10/500 2026-03-18 12:45 by Linda Hu
[考研] 一志愿中国海洋大学，生物学，301分，求调剂 +3 1孙悟空 2026-03-17 3/150 2026-03-18 10:28 by macy2011
[考研] 生物学071000 329分求调剂 +3 我爱生物生物爱� 2026-03-17 3/150 2026-03-18 10:12 by macy2011
[考研] 334求调剂 +3 志存高远意在机� 2026-03-16 3/150 2026-03-18 08:34 by lm4875102
[考研] 277调剂 +5 自由煎饼果子 2026-03-16 6/300 2026-03-17 19:26 by 李leezz
[考研] 268求调剂 +8 一定有学上- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 一志愿南京大学，080500材料科学与工程，调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 求调剂 +7 18880831720 2026-03-11 7/350 2026-03-13 16:10 by JourneyLucky
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
[考研] 290求调剂 +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生，2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
信息提示
关闭
请填处理意见
关闭