 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3190 | 回復: 16 | ||||||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||||||
[交流]
【求助】如何在siesta中設置ghost state原子的計算
|
||||||
|
用siesta計算體相金剛石中一個空位，形成能好像不準，單胞越大，形成能也越大，沒有合理的解釋。 師兄說，siesta用原子軌道基組，可能存在鬼態(tài)原子，需要在空位處設置一個假想的原子，我查了一下手冊，沒有完全搞清楚如何設置？ 不知道哪位大俠有做過？ |
回復此樓» 收錄本帖的淘帖專輯推薦
 siesta問題 | SIESTA | siesta基本概念 |
» 猜你喜歡
本人考085602 化學工程 專碩
已經(jīng)有23人回復
-
求調(diào)劑院校信息
已經(jīng)有4人回復
-
085600材料與化工306
已經(jīng)有4人回復
-
286求調(diào)劑
已經(jīng)有10人回復
-
328求調(diào)劑，英語六級551，有科研經(jīng)歷
已經(jīng)有9人回復
-
一志愿北京化工大學070300 學碩336求調(diào)劑
已經(jīng)有4人回復
-
286分人工智能專業(yè)請求調(diào)劑愿意跨考！
已經(jīng)有8人回復
-
資源與環(huán)境 調(diào)劑申請(333分)
已經(jīng)有5人回復
-
280求調(diào)劑
已經(jīng)有12人回復
-
269專碩求調(diào)劑
已經(jīng)有5人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 為什么采用CASTEP計算CaTiO3的光學性質其中的靜態(tài)介電常數(shù)與實驗值相差很大
已經(jīng)有14人回復
- 關于Discover動力學計算的一個問題
已經(jīng)有26人回復
- 如何計算空位部分所帶的電子
已經(jīng)有5人回復
- 動力學計算過程：時間步長、步數(shù)一般設置為多少？
已經(jīng)有10人回復
- siesta 計算電子占據(jù)數(shù)的問題
已經(jīng)有3人回復
- 關于SIESTA中如何計算功函數(shù)的一種做法探討
已經(jīng)有8人回復
- siesta計算弱相互作用VDW
已經(jīng)有17人回復
- 【求助】siesta進行ordern計算出錯求助
已經(jīng)有4人回復
- 【求助】siesta關于金屬計算的參數(shù)設置
已經(jīng)有5人回復
- 【求助】Ms計算磁性要不要設置每個原子的磁矩？？
已經(jīng)有10人回復
- 【求助】siesta 計算DOS應該注意什么問題
已經(jīng)有3人回復
» 搶金幣啦！回帖就可以得到:
-
青島高校2026招生碩士研究生
+5/770
-
福建師范大學葉應祥課題組招收2026級博士研究生
+2/150
-
福建農(nóng)林大學材料工程學院高分子材料課題組招生
+1/79
-
愿遇見你，遇見不止擦肩而過的緣分
+1/74
-
飄在美國，有點想家
+1/51
-
PBAT溶解做成乳液
+1/36
-
0854電子信息調(diào)劑，閩南師大光電芯片研發(fā)實驗室，柯少穎教授課題組，集成電路
+1/35
-
西安某大學能源化工新材料研發(fā)項目組接收2026級調(diào)劑研究生-計算機類
+1/32
-
武漢工程大學資源生物技術課題組生物/環(huán)境/化工/礦物專業(yè)碩士、博士研究生招生及調(diào)劑
+2/30
-
春眠不覺曉
+1/29
-
同濟大學環(huán)境學院 肖倩特聘研究員課題組 招聘碩士/博士（長期有效）
+1/22
-
招聘科研助理——儲能方向（儲能、流體傳熱、CFD模擬）
+1/16
-
【預碩】【預博】青島大學生物醫(yī)藥創(chuàng)新轉化
+1/15
-
湖南大學微生物結構與功能實驗室2026年計劃招收博士研究生
+1/15
-
南昌一本市政工程接受調(diào)劑
+1/12
-
國家級人才團隊誠聘生物信息相關專業(yè)的科研助理
+1/4
-
招收調(diào)劑生 國家重點實驗室（廣西某師范大學、B區(qū)）
+1/3
-
27年博士申請
+1/3
-
211/雙一流石河子大學化學化工學院--電催化，CO2捕集和轉化
+1/2
-
211/雙一流石河子大學化學化工學院-電池/電催化研究方向
+1/1
8樓2011-03-08 13:28:01
★ ★
zzy870720z(金幣+2): 謝謝指教 2011-03-07 22:20:35
hw45888792(金幣+4): 好厲害呀，好像回答過我好多問題啊~ 2011-03-08 10:18:18
zzy870720z(金幣+2): 謝謝指教 2011-03-07 22:20:35
hw45888792(金幣+4): 好厲害呀，好像回答過我好多問題啊~ 2011-03-08 10:18:18
|
參照手冊上的ChemicalSpeciesLabel的說明，對于gost原子的原子序數(shù)用一個負值就行，然后最好給他起一個別的元素名稱，贗勢也相應的改變一下。這個時候他的原子質量就會是1.0d30 a.u.如果你計算的是金剛石，那么gost原子的贗勢就用C的就行，和你就算中用的其他的C用一個贗勢就行，只不過最好起一個別的名字，然后對贗勢也重命名一下，將其元素序號用一個負值比如-6表示。 |
2樓2011-03-07 22:14:00
|
你好，這是我測試NV在金剛石2*2*2超單胞C62N中的測試輸入數(shù)據(jù)： 計算鬼態(tài)為Cg，在搭建結構的時候，我在空位處加上一個C原子，寫成輸入文件以后，將其改成Cg，標號為3。其中，C標號為1，N標號為2。贗勢方面，我將C.psf直接復制成Cg.psf，然后將Cg.psf里面的元素名稱改成Cg，不知是否可行？ SystemName C62N_ghost SystemLabel C62N_ghost NumberOfSpecies 3 %block ChemicalSpeciesLabel 1 6 C 2 7 N 3 -6 Cg %endblock ChemicalSpeciesLabel 如果是常規(guī)計算，就不用鬼態(tài)，輸入文件如下： SystemName C62N SystemLabel C62N NumberOfSpecies 2 %block ChemicalSpeciesLabel 1 6 C 2 7 N %endblock ChemicalSpeciesLabel 不知道是否正確？ 計算過程沒有任何問題，對比結果我也大概寫在下面： 加入鬼態(tài)： siesta: Program's energy decomposition (eV): siesta: Eions = 16680.576136 siesta: Ena = 3653.202863 siesta: Ekin = 7163.588024 siesta: Enl = -820.400928 siesta: DEna = -182.458272 siesta: DUscf = 15.575306 siesta: DUext = 0.000000 siesta: Exc = -3054.613637 siesta: eta*DQ = 0.000000 siesta: Emadel = 0.000000 siesta: Ekinion = 0.000000 siesta: Eharris = -9905.679557 siesta: Etot = -9905.682781 siesta: FreeEng = -9905.682783 siesta: Final energy (eV): siesta: Kinetic = 7163.588024 siesta: Hartree = 894.608032 siesta: Ext. field = 0.000000 siesta: Exch.-corr. = -3054.613637 siesta: Ion-electron = -8393.758885 siesta: Ion-ion = -6515.506314 siesta: Ekinion = 0.000000 siesta: Total = -9905.682781 不加鬼態(tài) siesta: Program's energy decomposition (eV): siesta: Eions = 16680.576136 siesta: Ena = 3653.202863 siesta: Ekin = 7164.385562 siesta: Enl = -820.554764 siesta: DEna = -182.943454 siesta: DUscf = 15.589130 siesta: DUext = 0.000000 siesta: Exc = -3054.694636 siesta: eta*DQ = 0.000000 siesta: Emadel = 0.000000 siesta: Ekinion = 0.000000 siesta: Eharris = -9905.587400 siesta: Etot = -9905.591435 siesta: FreeEng = -9905.591437 siesta: Final energy (eV): siesta: Kinetic = 7164.385562 siesta: Hartree = 894.861106 siesta: Ext. field = 0.000000 siesta: Exch.-corr. = -3054.694636 siesta: Ion-electron = -8394.637153 siesta: Ion-ion = -6515.506314 siesta: Ekinion = 0.000000 siesta: Total = -9905.591435 計算結果顯示，加入鬼態(tài)后，總能量略低0.1eV左右，而且鬼態(tài)原子附近有大約0.2單自旋電荷分布，其它C原子是2個電荷分布。 |
3樓2011-03-08 11:11:00
★ ★ ★ ★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
zzy870720z(金幣+5): 謝謝zhang兄指教，說的太好了 2011-03-08 15:56:38
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
zzy870720z(金幣+5): 謝謝zhang兄指教，說的太好了 2011-03-08 15:56:38
|
其實說句實話，我是真不厲害，很菜的。不是謙虛，我是自我表白。我在論壇里面比較活躍不是顯擺我知道的多，而是恰是我知道少，但是身邊有沒有人用這個軟件。所以只能依靠emuch論壇的這個版塊。通過討論可以相互交流一下，哪怕我的是錯誤的，我也可以說一下，這樣高手看見了，可以來指正我，這也是一個自我的進步。另外咱們這個版塊不是很活躍。貌似高手現(xiàn)在都不大來了，或者來了之后看到問題的幼稚就嗤之一笑，然后不懶的回復了。希望老蟲還是多來指點指點我們這些幼蟲們。 另外你的這個問題實際上可以通過多認真的讀一讀手冊是可以得到解決的。使用siesta工作的是手邊要么有一本打印的手冊。要么在電腦上打開一份電子的手冊。有什么問題，直接搜索一下相關的關鍵字。比如你的問題：直接在手冊里面用ctrl+F搜索ghost,立馬會定位到ChemicalSpeciesLabel位置，再搜，會定位到AtomicMass位置，接下來在PS.KBprojectors會發(fā)現(xiàn)對ghost設置的一些條件，這個需要檢查一下產(chǎn)生的PS.KBprojectors對ghost atom是不是只有一個，如果是的話，那就是程序給你處理好了。如果不是的話，計算出來預料不到結果，需要考慮這個問題。PAO.Basis在這兒，解釋了程序對ghost atom的基矢的處理辦法。 到此為止，應該對siesta處理ghost atom有一個了解。如果不是量化高手的話，可能對手冊上的解釋還是知道如何下手，這個時候需要google了。在google里面輸入ghost atom in siesta,一搜，google大多數(shù)時候不會辜負你的。 http://www.google.com.hk/search? ... i=&aql=&oq= 這個頭兩條就是siesta前輩們討論關于怎么設置ghost。 很懊惱的是打不開，這個時候我來告訴你，網(wǎng)頁代理來幫助：http://www.7daili.com/，首先打開這個網(wǎng)頁代理，然后打開你要打開但是打不的那個網(wǎng)頁，從瀏覽器上把這個網(wǎng)址復制到你的網(wǎng)頁代理里面，然后讓它給你打開。不要嫌棄這個太慢，耐心點，能打開總比打不開好多了。燈打開之后，看到里面討論的時候，欣喜不由從心底產(chǎn)生。在之內(nèi)的鏈接可以直接打開了（其實還是網(wǎng)頁代理給你打開的）。不管你從這些討論中得沒得到你要的結果，但是起碼，你學習了，搜索了，看到之前也有人討論這個問題了。這就夠了。如果問題解決了，皆大歡喜，如果還是沒有解決，在走下一條通往羅馬的大道�？傊@個過程是一個需要的過程，我認為也是發(fā)現(xiàn)問題，解決問題的作為學生的過程�？赡苣銜f，作為學生，時間短，論文壓力大。沒有時間花費在這個上面，直接問問知道就行。前提是如果有知道在身邊什么都好說。如果沒有呢？難道咱們的問題不不解決了？向小木蟲丟出來一個求助，不一定有人立馬回復你。這個期間還得自己去走上面的過程，只不過小木蟲的討論給上面的過程多了一個希望而已。 這些只是我的感想而已。我也是初學者。來到小木蟲，首先，要搜索一下關于你的這個話題的帖子。最起碼，來到小木蟲定居之后，要把之前討論過的有關你用的這個軟件的所有帖子全部看一遍，知道別人都討論了些什么。結果怎么樣，哪怕自己用不到，說不定以后用的時候，就需要回來再次檢索。平時沒事，一定要天天登錄一下小木蟲看一下自己用的這個程序有什么新的問題，別人是怎么解決的。這樣對該軟件在小木蟲的歷史就有了了解。對歷史不懂得人，會走很多彎路的。比如最近yashumimi發(fā)的那個帖子transista的錯誤，其實，別人之前早就討論過了。只是她沒有關心，以為自己不會碰到這個問題，壓根就沒有看這個帖子。但是她就應為自己對這個帖子的默然，也像mail-list求助了，也發(fā)帖子了。我也看到一段時間沒有人回復。估計期間她也費了好大勁，做了不少計算，并且花了不少機時。這個都是對歷史不解或者上小木蟲之后沒有首先檢索惹的禍。 最后，小木蟲不是一個有問題就往上丟，沒問題的時候壓根不來的地方，而是一個同行之間相互了解相互交流的平臺，窗口。就像一些學術會議一樣。只不過小木蟲上天天再開“學術會議”，“作報告的人”都是一些資源的提供者，更或者對我們這些板塊，“作報告的人”是一些遇到問題的人，報告的內(nèi)容就是問題本身而已。 當自己不是“作報告的人”時候要常來當聽報告的。因為這個問題被別人討論之后，再重復的提出來，可能會使得一些高手失去回答的耐心。不但對自己沒有好處，然高手失去了對咱們新一代新手的信心，對別的新手來說，是一個莫大的損失。 我也看到有些蟲友都工作了，都當教授了，還是常常忙里偷閑來答復一下新蟲的問題。我看了，確實有些感動。 最后，祝愿小木蟲SIESTA板塊，越來越紅火，能成就一些新手。希望老蟲多來哺養(yǎng)一下我們這些新蟲。 末了，言辭不妥之處，望新老蟲友見諒。  |
4樓2011-03-08 11:32:09
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600材料與化工306 +4 | z1z2z3879 2026-03-21 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 廣西大學材料導師推薦 +3 | 夏夏夏小正 2026-03-17 | 5/250 |
|
|
|
[考研] 化學工程321分求調(diào)劑 +18 | 大米飯！ 2026-03-15 | 22/1100 |
|
|
|
[考研] 326求調(diào)劑 +5 | 諾貝爾化學獎覬?/a> 2026-03-15 | 8/400 |
|
|
|
[考研]
 一志愿東華大學控制學碩320求調(diào)劑
+3
|
Grand777 2026-03-21 | 3/150 |
|
|
|
[考研] 326求調(diào)劑 +4 | mlpqaz03 2026-03-15 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +4 | 要好好無聊 2026-03-21 | 4/200 |
|
|
|
[考研] 297求調(diào)劑 +11 | 戲精丹丹丹 2026-03-17 | 12/600 |
|
|
|
[考研] 生物學一志愿985，分數(shù)349求調(diào)劑 +3 | zxts12 2026-03-21 | 3/150 |
|
|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
|
[考研] 華東師范大學-071000生物學-293分-求調(diào)劑 +3 | 研究生何瑤明 2026-03-18 | 3/150 |
|
|
|
[考研] 324分 085600材料化工求調(diào)劑 +4 | llllkkkhh 2026-03-18 | 4/200 |
|
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學-求調(diào)劑 +9 | 手機用戶 2026-03-17 | 9/450 |
|
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
|
[考研] 【同濟軟件】軟件（085405）考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
|
[考研]
0854，計算機類招收調(diào)劑
+3
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|