版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

identation

金蟲 (著名寫手)

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

[求助] 請(qǐng)問siesta計(jì)算BornCharge時(shí)為什么輸出的張量不滿足對(duì)稱性要求呢？謝謝

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

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.

計(jì)算的有效電荷如下：

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
這個(gè)我覺得應(yīng)該是沒三行屬于相同的原子的，因?yàn)橛行щ姾墒嵌A張量

回復(fù)此樓

» 猜你喜歡

342求調(diào)劑已經(jīng)有3人回復(fù)
070300化學(xué)專業(yè)279調(diào)劑已經(jīng)有9人回復(fù)
求調(diào)劑：085600材料與化工，考材科基，總分319 已經(jīng)有20人回復(fù)
一志愿華東理工大學(xué)，080500學(xué)碩，317分，求調(diào)劑已經(jīng)有9人回復(fù)
一志愿a區(qū)211，085601-307分求調(diào)劑已經(jīng)有23人回復(fù)
合肥區(qū)域性重點(diǎn)一本招收調(diào)劑已經(jīng)有8人回復(fù)
材料調(diào)劑已經(jīng)有10人回復(fù)
288求調(diào)劑一志愿哈工大材料與化工已經(jīng)有13人回復(fù)
材料科學(xué)與工程339求調(diào)劑已經(jīng)有3人回復(fù)
化學(xué)308分調(diào)劑已經(jīng)有11人回復(fù)

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

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

巴蒂LBJ

木蟲 (知名作家)

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

幫忙頂貼

回復(fù)此樓

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

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

iamikaruk

木蟲 (著名寫手)

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

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

贊一下

回復(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強(qiáng)帖: 16
應(yīng)助: 46 (小學(xué)生)
貴賓: 2.516
金幣: 253.4
散金: 2603
紅花: 44
帖子: 1283
在線: 500.5小時(shí)
蟲號(hào): 842155
注冊(cè): 2009-09-06
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

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

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

贊一下

回復(fù)此樓

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

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

iamikaruk

木蟲 (著名寫手)

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

引用回帖:

4樓: Originally posted by identation at 2011-08-02 22:34:07:
我這個(gè)體系如果保留小數(shù)點(diǎn)后三位，對(duì)于Ba和Ti是對(duì)的但是對(duì)于O原子就有兩個(gè)分量了這是怎么回事呢

很正常的，O和Ti有兩種關(guān)系，一種是平行于Ti-O鍵，一種是垂直于Ti-O鍵，所以有兩個(gè)分量

贊一下

回復(fù)此樓

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

5樓2011-08-02 23:43:00

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

identation

金蟲 (著名寫手)

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

引用回帖:

5樓: Originally posted by iamikaruk at 2011-08-02 23:43:00:
很正常的，O和Ti有兩種關(guān)系，一種是平行于Ti-O鍵，一種是垂直于Ti-O鍵，所以有兩個(gè)分量

你怎么不應(yīng)助呢沒法給你金幣啊嘿嘿

贊一下

回復(fù)此樓

6樓2011-08-03 00:06:09

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

identation

金蟲 (著名寫手)

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

引用回帖:

6樓: Originally posted by identation at 2011-08-03 00:06:09:
你怎么不應(yīng)助呢沒法給你金幣啊嘿嘿

當(dāng)前回帖不是應(yīng)助回帖，您不能對(duì)此進(jìn)行金幣獎(jiǎng)勵(lì)哦

贊一下

回復(fù)此樓

7樓2011-08-03 00:06:24

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

iamikaruk

木蟲 (著名寫手)

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

【答案】應(yīng)助回帖

identation(金幣+15): 2011-08-03 00:59:48

ok，好了。填應(yīng)助了。話說除了abinit外，難道就沒有一個(gè)第一性原理軟件可以把對(duì)稱性處理的很好嗎？

贊一下

回復(fù)此樓

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

8樓2011-08-03 00:47:57

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

identation

金蟲 (著名寫手)

1ST強(qiáng)帖: 16
應(yīng)助: 46 (小學(xué)生)
貴賓: 2.516
金幣: 253.4
散金: 2603
紅花: 44
帖子: 1283
在線: 500.5小時(shí)
蟲號(hào): 842155
注冊(cè): 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外，難道就沒有一個(gè)第一性原理軟件可以把對(duì)稱性處理的很好嗎？

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

贊一下

回復(fù)此樓

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

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

iamikaruk

木蟲 (著名寫手)

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

vasp的Born effective charge對(duì)稱性也不好，本來三個(gè)張量分量根據(jù)對(duì)稱性應(yīng)該一樣的，它硬是弄成不一樣了。

贊一下

回復(fù)此樓

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

10樓2011-08-03 09:40:45

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 identation 的主題更新

返回列表

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

[考研] 合肥區(qū)域性重點(diǎn)一本招收調(diào)劑 +4	6266jl 2026-03-30	8/400	2026-03-31 18:43 by 6266jl
[考研] 285求調(diào)劑 +3	FZAC123 2026-03-30	3/150	2026-03-31 17:49 by 熱情沙漠
[考研] 271求調(diào)劑 +5	勒布朗@ 2026-03-31	7/350	2026-03-31 16:26 by hypershenger
[考研] 343求調(diào)劑 +8	愛羈絆 2026-03-28	8/400	2026-03-31 16:12 by 不吃魚的貓
[考研] 334求調(diào)劑 +7	Trying] 2026-03-31	7/350	2026-03-31 12:33 by 無際的草原
[考研] 福建理工大學(xué)材料學(xué)院先進(jìn)合金團(tuán)隊(duì)招收考研調(diào)劑學(xué)生 +3	大華金商都 2026-03-30	4/200	2026-03-31 01:04 by 方英俊602
[考研] 085602化工求調(diào)劑（331分） +8	111@127 2026-03-30	8/400	2026-03-30 21:23 by 研究僧導(dǎo)導(dǎo)
[考研] 抱歉 +3	田洪有 2026-03-30	3/150	2026-03-30 19:11 by 迷糊CCPs
[考研] 334分一志愿武理材料求調(diào)劑 +16	李李不服輸 2026-03-26	16/800	2026-03-30 17:47 by wangjy2002
[考研] 343求調(diào)劑085601 +3	要努力學(xué)習(xí)x 2026-03-29	3/150	2026-03-29 18:35 by wxiongid
[考研] 調(diào)劑考研 +3	王杰一 2026-03-29	3/150	2026-03-29 08:09 by fmesaito
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +8	邱gl 2026-03-27	8/400	2026-03-28 12:42 by 唐沐兒
[考研] 286求調(diào)劑 +12	PolarBear11 2026-03-26	12/600	2026-03-28 12:14 by zllcz
[考研] 292求調(diào)劑 +14	鵝鵝鵝額額額額?/a> 2026-03-25	15/750	2026-03-28 08:45 by WYUMater
[考研] 0856調(diào)劑 +5	求求讓我有書讀?/a> 2026-03-26	6/300	2026-03-27 15:12 by caszguilin
[考研] 298調(diào)劑 +3	jiyingjie123 2026-03-27	3/150	2026-03-27 11:57 by wxiongid
[論文投稿] Journal of Mechanical Science and Technology +3	Russ_ss 2026-03-25	5/250	2026-03-27 10:49 by 陸小果畫大餅
[考研] 一志愿吉大071010，316分求調(diào)劑 +3	xgbiknn 2026-03-27	3/150	2026-03-27 10:36 by guoweigw
[考研] 324求調(diào)劑 +5	hanamiko 2026-03-26	5/250	2026-03-27 10:33 by wangjy2002
[考研] 一志愿南京郵電大學(xué) 288分材料考研求調(diào)劑 +3	jl0720 2026-03-26	3/150	2026-03-26 13:39 by zzll406

24小時(shí)熱門版塊排行榜

[求助] 請(qǐng)問siesta計(jì)算BornCharge時(shí)為什么輸出的張量不滿足對(duì)稱性要求呢？謝謝

» 猜你喜歡

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

【答案】應(yīng)助回帖

[求助] 請(qǐng)問siesta計(jì)算BornCharge時(shí)為什么輸出的張量不滿足對(duì)稱性要求呢？謝謝

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