identation

金蟲 (著名寫手)

1ST強帖: 16
應(yīng)助: 46 (小學生)
貴賓: 2.516
金幣: 253.4
散金: 2603
紅花: 44
帖子: 1283
在線: 500.5小時
蟲號: 842155
注冊: 2009-09-06
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[求助] 請問siesta計算BornCharge時為什么輸出的張量不滿足對稱性要求呢？謝謝

我測試siesta，計算立方相BaTiO3的波恩有效電荷，結(jié)果輸出的張量中分量很多啊，立方晶每個原子的有效電荷應(yīng)該只有一個分量啊，另外，立方BaTiO3的三個O原子應(yīng)該是等價的啊，為什么計算出的有效電荷不一樣呢？
感激不盡！求指教！
輸入文件如下：

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
這個我覺得應(yīng)該是沒三行屬于相同的原子的，因為有效電荷是二階張量

回復(fù)此樓

» 猜你喜歡

288求調(diào)劑一志愿哈工大材料與化工已經(jīng)有16人回復(fù)
0805求調(diào)劑已經(jīng)有4人回復(fù)
求調(diào)劑，一志愿北林食品與營養(yǎng)095500，301分，已過六級，有科研經(jīng)歷已經(jīng)有3人回復(fù)
0856材料與化工調(diào)劑，339 已經(jīng)有5人回復(fù)
311（085601）求調(diào)劑已經(jīng)有10人回復(fù)
080200專業(yè)277分，求帶走！已經(jīng)有4人回復(fù)
359求調(diào)劑已經(jīng)有6人回復(fù)
342求調(diào)劑已經(jīng)有6人回復(fù)
070300化學279求調(diào)劑已經(jīng)有10人回復(fù)
086000生物與醫(yī)藥298調(diào)劑求助已經(jīng)有5人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

《張量分析》田宗若.西北工業(yè)大學出版社已經(jīng)有163人回復(fù)
什么是張量，基本思想是什么？已經(jīng)有14人回復(fù)
什么是張量？基本的思想是什么呢已經(jīng)有13人回復(fù)
期盼地球物理專業(yè)老師同學們進行矩張量反演方面的指點已經(jīng)有13人回復(fù)
通向凝聚態(tài)物理和量子化學的終極之路--密度矩陣重整化群，從矩陣乘積態(tài)到張量網(wǎng)態(tài) 已經(jīng)有19人回復(fù)
【求助】用vasp算壓電張量（系數(shù)），彈性張量需要注意的問題？已經(jīng)有9人回復(fù)
【求助】請問有人用WIEN2K算過光電導張量嗎？已經(jīng)有4人回復(fù)
【請教】關(guān)于電子順磁光譜EPR譜圖的g張量是如何求得的已經(jīng)有5人回復(fù)
【求助】張量的運算符號（塑性細觀力學一書p121）已經(jīng)有9人回復(fù)
【求助】張量運算請教已經(jīng)有7人回復(fù)
【求助】關(guān)于發(fā)表張量的一篇文章一般都找哪些雜志呀，想投EI類的已經(jīng)有3人回復(fù)

1樓 2011-08-02 18:59:00

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

identation

金蟲 (著名寫手)

1ST強帖: 16
應(yīng)助: 46 (小學生)
貴賓: 2.516
金幣: 253.4
散金: 2603
紅花: 44
帖子: 1283
在線: 500.5小時
蟲號: 842155
注冊: 2009-09-06
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

8樓: Originally posted by iamikaruk at 2011-08-03 00:47:57:
ok，好了。填應(yīng)助了。話說除了abinit外，難道就沒有一個第一性原理軟件可以把對稱性處理的很好嗎？

嘿嘿 siesta好像不能識別對稱性的 pwscf簡化計算量的時候能處理晶體結(jié)構(gòu)生成還是不行

贊一下

回復(fù)此樓

9樓2011-08-03 01:00:33

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 16 個回答

巴蒂LBJ

木蟲 (知名作家)

應(yīng)助: 5 (幼兒園)
金幣: 3028.3
散金: 1107
紅花: 12
沙發(fā): 38
帖子: 7098
在線: 321.5小時
蟲號: 1165239
注冊: 2010-12-08
性別: GG
專業(yè): 配位化學

幫忙頂貼

回復(fù)此樓

2樓2011-08-02 19:18:31

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

iamikaruk

木蟲 (著名寫手)

1ST強帖: 1
應(yīng)助: 11 (小學生)
金幣: 4950.3
散金: 118
紅花: 76
帖子: 1362
在線: 543.8小時
蟲號: 1315985
注冊: 2011-06-05
專業(yè): 無機非金屬類電介質(zhì)與電解

對稱性不太嚴格，但是刨除那些小的數(shù)值誤差，Born有效電荷張量基本上是對的

贊一下

回復(fù)此樓

[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]

3樓2011-08-02 21:00:38

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

identation

金蟲 (著名寫手)

1ST強帖: 16
應(yīng)助: 46 (小學生)
貴賓: 2.516
金幣: 253.4
散金: 2603
紅花: 44
帖子: 1283
在線: 500.5小時
蟲號: 842155
注冊: 2009-09-06
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

3樓: Originally posted by iamikaruk at 2011-08-02 21:00:38:
對稱性不太嚴格，但是刨除那些小的數(shù)值誤差，Born有效電荷張量基本上是對的

我這個體系如果保留小數(shù)點后三位，對于Ba和Ti是對的但是對于O原子就有兩個分量了這是怎么回事呢

贊一下

回復(fù)此樓

4樓2011-08-02 22:34:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 085602 307分求調(diào)劑 +10	不知道叫什么！ 2026-03-26	10/500	2026-03-31 19:53 by Dyhoer
[考研] 材料科學與工程求調(diào)劑 +13	深V宿舍吧 2026-03-29	13/650	2026-03-31 19:50 by Dyhoer
[考研] 317求調(diào)劑 +12	蛋黃咸肉粽 2026-03-26	13/650	2026-03-31 16:44 by asdfzly
[考研] 322求調(diào)劑 +3	熹僖XX 2026-03-31	3/150	2026-03-31 16:32 by 記事本2026
[考研] 266求調(diào)劑 +3	哇呼哼呼哼 2026-03-29	3/150	2026-03-31 10:06 by cal0306
[考研] 335求調(diào)劑 +3	321* 2026-03-31	3/150	2026-03-31 09:00 by 功夫瘋狂
[考研] 085601一志愿西北工業(yè)大學初試346 +4	085601初試346 2026-03-30	4/200	2026-03-31 07:47 by jp9609
[考研] 085701環(huán)境工程求調(diào)劑 +11	多久上課 2026-03-27	12/600	2026-03-30 21:21 by 研究僧導導
[考研] 0703一志愿9，初試成績：338，四六級已過，有科研經(jīng)歷，求調(diào)劑！ +7	Zuhui0306 2026-03-25	7/350	2026-03-30 19:01 by 源_2020
[考研] 295材料工程專碩求調(diào)劑 +10	1428151015 2026-03-27	10/500	2026-03-30 19:00 by 源_2020
[考研] 105500藥學求調(diào)劑，一志愿山東大學藥學，348分 +3	gr哈哈哈 2026-03-28	3/150	2026-03-30 18:56 by 源_2020
[考研] 334分一志愿武理材料求調(diào)劑 +16	李李不服輸 2026-03-26	16/800	2026-03-30 17:47 by wangjy2002
[考研] 317分一志愿南理工材料工程本科湖工大求調(diào)劑 +12	芋泥小鈴鐺 2026-03-28	12/600	2026-03-30 17:06 by wangjy2002
[考研] 394求調(diào)劑 +3	好事多磨靜候佳?/a> 2026-03-26	5/250	2026-03-28 14:24 by 唐沐兒
[考研] 339求調(diào)劑，想調(diào)回江蘇 +6	烤麥芽 2026-03-27	8/400	2026-03-28 10:40 by 烤麥芽
[考研] 0856調(diào)劑 +5	求求讓我有書讀?/a> 2026-03-26	6/300	2026-03-27 15:12 by caszguilin
[考研] 調(diào)劑 +4	柚柚yoyo 2026-03-26	4/200	2026-03-26 20:43 by fmesaito
[考研] 打過很多競賽，085406控制工程300分，求調(diào)劑 +3	askeladz 2026-03-26	3/150	2026-03-26 09:08 by 給你你注意休息
[考研] 07化學303求調(diào)劑 +5	睿08 2026-03-25	5/250	2026-03-25 22:46 by 418490947
[考研] 材料與化工304求B區(qū)調(diào)劑 +3	邱gl 2026-03-25	3/150	2026-03-25 19:03 by Ainin_