 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2217 | 回復: 15 | ||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
[求助]
請問siesta計算BornCharge時為什么輸出的張量不滿足對稱性要求呢？謝謝
|
||
|
我測試siesta，計算立方相BaTiO3的波恩有效電荷，結(jié)果輸出的張量中分量很多啊，立方晶每個原子的有效電荷應該只有一個分量啊，另外，立方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 這個我覺得應該是沒三行屬于相同的原子的，因為有效電荷是二階張量 |
回復此樓» 猜你喜歡
288求調(diào)劑 一志愿哈工大 材料與化工
已經(jīng)有16人回復
-
0805求調(diào)劑
已經(jīng)有4人回復
-
求調(diào)劑，一志愿北林食品與營養(yǎng)095500，301分，已過六級，有科研經(jīng)歷
已經(jīng)有3人回復
-
0856材料與化工調(diào)劑，339
已經(jīng)有5人回復
-
311（085601）求調(diào)劑
已經(jīng)有10人回復
-
080200專業(yè)277分，求帶走！
已經(jīng)有4人回復
-
359求調(diào)劑
已經(jīng)有6人回復
-
342求調(diào)劑
已經(jīng)有6人回復
-
070300化學279求調(diào)劑
已經(jīng)有10人回復
-
086000生物與醫(yī)藥298調(diào)劑求助
已經(jīng)有5人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 《張量分析》田宗若.西北工業(yè)大學出版社
已經(jīng)有163人回復
- 什么是張量，基本思想是什么？
已經(jīng)有14人回復
- 什么是張量？基本的思想是什么呢
已經(jīng)有13人回復
- 期盼地球物理專業(yè)老師同學們進行矩張量反演方面的指點
已經(jīng)有13人回復
- 通向凝聚態(tài)物理和量子化學的終極之路--密度矩陣重整化群，從矩陣乘積態(tài)到張量網(wǎng)態(tài)
已經(jīng)有19人回復
- 【求助】用vasp算壓電張量（系數(shù)），彈性張量需要注意的問題？
已經(jīng)有9人回復
- 【求助】請問有人用WIEN2K算過光電導張量嗎？
已經(jīng)有4人回復
- 【請教】關于電子順磁光譜EPR譜圖的g張量是如何求得的
已經(jīng)有5人回復
- 【求助】張量的運算符號（塑性細觀力學一書p121）
已經(jīng)有9人回復
- 【求助】張量運算請教
已經(jīng)有7人回復
- 【求助】關于發(fā)表張量的一篇文章一般都找哪些雜志呀，想投EI類的
已經(jīng)有3人回復
7樓2011-08-03 00:06:24
巴蒂LBJ
木蟲 (知名作家)
- 應助: 5 (幼兒園)
- 金幣: 3028.3
- 散金: 1107
- 紅花: 12
- 沙發(fā): 38
- 帖子: 7098
- 在線: 321.5小時
- 蟲號: 1165239
- 注冊: 2010-12-08
- 性別: GG
- 專業(yè): 配位化學
幫忙頂貼 2樓2011-08-02 19:18:31
iamikaruk
木蟲 (著名寫手)
- 1ST強帖: 1
- 應助: 11 (小學生)
- 金幣: 4950.3
- 散金: 118
- 紅花: 76
- 帖子: 1362
- 在線: 543.8小時
- 蟲號: 1315985
- 注冊: 2011-06-05
- 專業(yè): 無機非金屬類電介質(zhì)與電解
3樓2011-08-02 21:00:38
4樓2011-08-02 22:34:07
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 08工科275分求調(diào)劑 +5 | AaAa7420 2026-03-31 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 375求調(diào)劑 +7 | 雨夏整夜 2026-03-29 | 7/350 |
|
|
|
[考研] 合肥區(qū)域性重點一本招收調(diào)劑 +4 | 6266jl 2026-03-30 | 8/400 |
|
|
|
[考研] 086000調(diào)劑 +5 | 7901117076 2026-03-26 | 5/250 |
|
|
|
[考研] 調(diào)劑 +4 | 好好讀書。 2026-03-28 | 5/250 |
|
|
|
[考研] 本科211安全工程，初試290分，求調(diào)劑 +3 | 2719846834 2026-03-28 | 3/150 |
|
|
|
[考研] 0703化學 +20 | 妮妮ninicgb 2026-03-27 | 20/1000 |
|
|
|
[考研] 272求調(diào)劑，接受跨專業(yè)調(diào)劑！ +3 | 閑魚盧 2026-03-31 | 3/150 |
|
|
|
[考研] 0817化工學碩調(diào)劑 +7 | 努力上岸中！ 2026-03-31 | 7/350 |
|
|
|
[考研] 085601一志愿西北工業(yè)大學初試346 +4 | 085601初試346 2026-03-30 | 4/200 |
|
|
|
[考研] 370求調(diào)劑 +3 | 080700調(diào)劑 2026-03-30 | 3/150 |
|
|
|
[有機交流]
 甲基亞磺磺酸鈉和甲基磺酸酯反應機理
 10+3
|
kaobao456 2026-03-29 | 4/200 |
|
|
|
[考研] 303求調(diào)劑 +7 | DLkz1314. 2026-03-30 | 7/350 |
|
|
|
[考研] 340求調(diào)劑 +6 | Amber00 2026-03-26 | 6/300 |
|
|
|
[考研] 356求調(diào)劑 +3 | gysy?s?a 2026-03-28 | 3/150 |
|
|
|
[考研] 316求調(diào)劑 +7 | 江辭666 2026-03-26 | 7/350 |
|
|
|
[考研] 315分求調(diào)劑 +7 | 26考研上岸版26 2026-03-26 | 7/350 |
|
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫 2026-03-26 | 5/250 |
|
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
|
[考研] 一志愿 南京郵電大學 288分 材料考研 求調(diào)劑 +3 | jl0720 2026-03-26 | 3/150 |
|