 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3178 | 回復: 16 | ||||||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||||||
[交流]
【求助】如何在siesta中設置ghost state原子的計算
|
||||||
|
用siesta計算體相金剛石中一個空位，形成能好像不準，單胞越大，形成能也越大，沒有合理的解釋。 師兄說，siesta用原子軌道基組，可能存在鬼態(tài)原子，需要在空位處設置一個假想的原子，我查了一下手冊，沒有完全搞清楚如何設置？ 不知道哪位大俠有做過？ |
回復此樓» 收錄本帖的淘帖專輯推薦
 siesta問題 | SIESTA | siesta基本概念 |
» 猜你喜歡
環(huán)境領域全國重點實驗室招收博士1-2名
已經(jīng)有3人回復
-
材料專碩306英一數(shù)二
已經(jīng)有10人回復
-
301求調(diào)劑
已經(jīng)有6人回復
-
一志愿天津大學化學工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有7人回復
-
302求調(diào)劑
已經(jīng)有6人回復
-
26博士申請
已經(jīng)有3人回復
-
268求調(diào)劑
已經(jīng)有3人回復
-
311求調(diào)劑
已經(jīng)有10人回復
-
被我言中：新模板不強調(diào)格式了，假專家開始管格式了
已經(jīng)有4人回復
-
303求調(diào)劑
已經(jīng)有5人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 為什么采用CASTEP計算CaTiO3的光學性質(zhì)其中的靜態(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人回復
» 搶金幣啦！回帖就可以得到:
-
西湖大學工學院智能微型機器人實驗室博士后招聘啟事
+4/228
-
溫州大學呂晶晶課題組2026年招生
+1/84
-
太原理工大學博士招生 電池負極材料 固態(tài)電解質(zhì)方向 有工作經(jīng)驗者優(yōu)先
+1/79
-
博士后出站經(jīng)驗請教
+1/78
-
光學工程學碩調(diào)劑信息
+2/52
-
Call for paper: [Fermentation – IF 3.3]
+2/38
-
催化方向推薦一個極好的博導
+1/36
-
上海交通大學電催化方向博后招聘-有鋰電、高分子、燃料電池背景者優(yōu)先
+1/32
-
北京某研究院生物化學與分子生物學相關專業(yè)申請考核制博士招生（春季）
+1/15
-
福建師范大學化學與材料學院杜克釗團隊博士/碩士招生
+1/14
-
南京醫(yī)科大學-膜蛋白結(jié)構(gòu)與功能課題組-招收2026級博士研究生
+1/11
-
【全獎招生】北師港浸大ESLAS實驗室招收密碼工程/網(wǎng)絡安全/計算機視覺博士/博士后
+1/10
-
英國�？巳卮髮W & 法國巴黎薩克雷大學聯(lián)合培養(yǎng)博士
+1/7
-
國家杰青低維材料與器件力學團隊2026年招收博士研究生
+1/6
-
南方科技大學生命科學學院基礎免疫與微生物學系招收博士后
+1/4
-
江蘇南京林業(yè)大學-高分子功能材料和器件-招2名碩士研究生
+1/3
-
求推薦性價比好用的球磨罐
+1/2
-
2026年浙江清華長三角研究院智慧能源中心招聘博士后-電氣方向
+1/2
-
福建師范大學環(huán)境微生物技術課題組接受2026級碩士調(diào)劑
+1/2
-
招收理論凝聚態(tài)物理/納米光學/量子計算方向博士、碩士研究生/博士后
+1/1
|
你好，這是我測試NV在金剛石2*2*2超單胞C62N中的測試輸入數(shù)據(jù)： 計算鬼態(tài)為Cg，在搭建結(jié)構(gòu)的時候，我在空位處加上一個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 不知道是否正確？ 計算過程沒有任何問題，對比結(jié)果我也大概寫在下面： 加入鬼態(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 計算結(jié)果顯示，加入鬼態(tài)后，總能量略低0.1eV左右，而且鬼態(tài)原子附近有大約0.2單自旋電荷分布，其它C原子是2個電荷分布。 |
3樓2011-03-08 11:11:00
★ ★
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ù)用一個負值就行，然后最好給他起一個別的元素名稱，贗勢也相應的改變一下。這個時候他的原子質(zhì)量就會是1.0d30 a.u.如果你計算的是金剛石，那么gost原子的贗勢就用C的就行，和你就算中用的其他的C用一個贗勢就行，只不過最好起一個別的名字，然后對贗勢也重命名一下，將其元素序號用一個負值比如-6表示。 |
2樓2011-03-07 22:14: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是不是只有一個，如果是的話，那就是程序給你處理好了。如果不是的話，計算出來預料不到結(jié)果，需要考慮這個問題。PAO.Basis在這兒，解釋了程序?qū)host 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)頁代理給你打開的）。不管你從這些討論中得沒得到你要的結(jié)果，但是起碼，你學習了，搜索了，看到之前也有人討論這個問題了。這就夠了。如果問題解決了，皆大歡喜，如果還是沒有解決，在走下一條通往羅馬的大道。總之，這個過程是一個需要的過程，我認為也是發(fā)現(xiàn)問題，解決問題的作為學生的過程�？赡苣銜f，作為學生，時間短，論文壓力大。沒有時間花費在這個上面，直接問問知道就行。前提是如果有知道在身邊什么都好說。如果沒有呢？難道咱們的問題不不解決了？向小木蟲丟出來一個求助，不一定有人立馬回復你。這個期間還得自己去走上面的過程，只不過小木蟲的討論給上面的過程多了一個希望而已。 這些只是我的感想而已。我也是初學者。來到小木蟲，首先，要搜索一下關于你的這個話題的帖子。最起碼，來到小木蟲定居之后，要把之前討論過的有關你用的這個軟件的所有帖子全部看一遍，知道別人都討論了些什么。結(jié)果怎么樣，哪怕自己用不到，說不定以后用的時候，就需要回來再次檢索。平時沒事，一定要天天登錄一下小木蟲看一下自己用的這個程序有什么新的問題，別人是怎么解決的。這樣對該軟件在小木蟲的歷史就有了了解。對歷史不懂得人，會走很多彎路的。比如最近yashumimi發(fā)的那個帖子transista的錯誤，其實，別人之前早就討論過了。只是她沒有關心，以為自己不會碰到這個問題，壓根就沒有看這個帖子。但是她就應為自己對這個帖子的默然，也像mail-list求助了，也發(fā)帖子了。我也看到一段時間沒有人回復。估計期間她也費了好大勁，做了不少計算，并且花了不少機時。這個都是對歷史不解或者上小木蟲之后沒有首先檢索惹的禍。 最后，小木蟲不是一個有問題就往上丟，沒問題的時候壓根不來的地方，而是一個同行之間相互了解相互交流的平臺，窗口。就像一些學術會議一樣。只不過小木蟲上天天再開“學術會議”，“作報告的人”都是一些資源的提供者，更或者對我們這些板塊，“作報告的人”是一些遇到問題的人，報告的內(nèi)容就是問題本身而已。 當自己不是“作報告的人”時候要常來當聽報告的。因為這個問題被別人討論之后，再重復的提出來，可能會使得一些高手失去回答的耐心。不但對自己沒有好處，然高手失去了對咱們新一代新手的信心，對別的新手來說，是一個莫大的損失。 我也看到有些蟲友都工作了，都當教授了，還是常常忙里偷閑來答復一下新蟲的問題。我看了，確實有些感動。 最后，祝愿小木蟲SIESTA板塊，越來越紅火，能成就一些新手。希望老蟲多來哺養(yǎng)一下我們這些新蟲。 末了，言辭不妥之處，望新老蟲友見諒。  |
4樓2011-03-08 11:32:09
★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
zzy870720z(金幣+2): 繼續(xù)獎勵 2011-03-08 15:56:51
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
zzy870720z(金幣+2): 繼續(xù)獎勵 2011-03-08 15:56:51
|
其實，我之前從來沒有做過ghost state的計算，并且對ghost state的知識和用處不了解。僅僅是我之前看手冊的時候記得手冊中提到了ghost atom的處理。但是具體怎么處理，我是沒往心里記的。昨天晚上臨回宿舍前，上了一下小木蟲，發(fā)現(xiàn)了你的這個問題。怕你急于解答，所以我就臨時的翻看了一下手冊，并且到網(wǎng)上搜索了一下ghos state在siesta中的討論而已。具體的加入ghost atom之后計算結(jié)果應該是什么樣子是預期，我不清楚，這個應該你是高手。 需要指明的是，C_g.psf里面的內(nèi)容不需要改的。表頭部分的，程序是不讀的，只是給使用者一個參數(shù)介紹，讓使用者知道這個贗勢的類型和參數(shù)而已。 |
5樓2011-03-08 11:39:57
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +8 | 一鴨鴨喲 2026-03-14 | 10/500 |
|
|
[考研] 278求調(diào)劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
|
[考研] 286求調(diào)劑 +3 | lemonzzn 2026-03-16 | 5/250 |
|
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
|
[考研] 一志愿華中師范071000，325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
|
[考博] 東華理工大學化材專業(yè)26屆碩士博士申請 +6 | zlingli 2026-03-13 | 6/300 |
|
|
|
[考研] 中科大材料與化工319求調(diào)劑 +3 | 孟鑫材料 2026-03-14 | 3/150 |
|
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
|
[考研] 255求調(diào)劑 +3 | 李嘉慧， 2026-03-12 | 4/200 |
|
|
|
[考研] 【0703化學調(diào)劑】-一志愿華中師范大學-六級475 +5 | Becho359 2026-03-11 | 5/250 |
|
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
|
[考研] 材料371求調(diào)劑 +9 | 鱷魚? 2026-03-11 | 11/550 |
|
|
|
[考研] 一志愿西南交大，材料專碩317求調(diào)劑 +5 | lx8568 2026-03-11 | 5/250 |
|
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|
|
|
[考研] 一志愿211化學學碩310分求調(diào)劑 +8 | 努力奮斗112 2026-03-12 | 9/450 |
|
|
|
[論文投稿]
投稿問題
5+4
|
星光燦爛xt 2026-03-12 | 6/300 |
|
|
|
[考研] 290求調(diào)劑 +3 | ADT 2026-03-13 | 3/150 |
|
|
|
[考研] 321求調(diào)劑（食品/專碩） +3 | xc321 2026-03-12 | 6/300 |
|
|
|
[考博] 2026年博士申請 +3 | QwQwQW10 2026-03-11 | 3/150 |
|