 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3173 | 回復(fù): 16 | ||||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||||
[交流]
【求助】如何在siesta中設(shè)置ghost state原子的計算
|
||||||
|
用siesta計算體相金剛石中一個空位，形成能好像不準(zhǔn)，單胞越大，形成能也越大，沒有合理的解釋。 師兄說，siesta用原子軌道基組，可能存在鬼態(tài)原子，需要在空位處設(shè)置一個假想的原子，我查了一下手冊，沒有完全搞清楚如何設(shè)置？ 不知道哪位大俠有做過？ |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 siesta問題 | SIESTA | siesta基本概念 |
» 猜你喜歡
能源材料化學(xué)課題組招收碩士研究生8-10名
已經(jīng)有6人回復(fù)
-
296求調(diào)劑
已經(jīng)有5人回復(fù)
-
0703化學(xué)調(diào)劑 ，六級已過，有科研經(jīng)歷
已經(jīng)有8人回復(fù)
-
268求調(diào)劑
已經(jīng)有8人回復(fù)
-
材料，紡織，生物（0856、0710），化學(xué)招生啦
已經(jīng)有7人回復(fù)
-
301求調(diào)劑
已經(jīng)有5人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有7人回復(fù)
-
材料與化工求調(diào)劑
已經(jīng)有6人回復(fù)
-
299求調(diào)劑
已經(jīng)有4人回復(fù)
-
277調(diào)劑
已經(jīng)有6人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 為什么采用CASTEP計算CaTiO3的光學(xué)性質(zhì)其中的靜態(tài)介電常數(shù)與實(shí)驗(yàn)值相差很大
已經(jīng)有14人回復(fù)
- 關(guān)于Discover動力學(xué)計算的一個問題
已經(jīng)有26人回復(fù)
- 如何計算空位部分所帶的電子
已經(jīng)有5人回復(fù)
- 動力學(xué)計算過程：時間步長、步數(shù)一般設(shè)置為多少？
已經(jīng)有10人回復(fù)
- siesta 計算電子占據(jù)數(shù)的問題
已經(jīng)有3人回復(fù)
- 關(guān)于SIESTA中如何計算功函數(shù)的一種做法探討
已經(jīng)有8人回復(fù)
- siesta計算弱相互作用VDW
已經(jīng)有17人回復(fù)
- 【求助】siesta進(jìn)行ordern計算出錯求助
已經(jīng)有4人回復(fù)
- 【求助】siesta關(guān)于金屬計算的參數(shù)設(shè)置
已經(jīng)有5人回復(fù)
- 【求助】Ms計算磁性要不要設(shè)置每個原子的磁矩？？
已經(jīng)有10人回復(fù)
- 【求助】siesta 計算DOS應(yīng)該注意什么問題
已經(jīng)有3人回復(fù)
» 搶金幣啦！回帖就可以得到:
-
南京大學(xué)智能驅(qū)動與感知材料實(shí)驗(yàn)室誠招博士/碩士/科研助理/訪問學(xué)生
+1/280
-
華南師范大學(xué)（211）博士招生- 電子、自動化、機(jī)械、生物學(xué)、物理相關(guān)專業(yè)
+2/42
-
山東師范大學(xué)有機(jī)化學(xué)專業(yè)胡忠燕老師課題組招收2026屆碩士研究生以及調(diào)劑生
+1/40
-
中國科學(xué)院杭州醫(yī)學(xué)研究所覃江江課題組招聘博士后
+1/36
-
湘潭大學(xué)“過程強(qiáng)化與綠色化工”創(chuàng)新團(tuán)隊補(bǔ)招2026年秋入學(xué)博士生
+1/35
-
華南師范大學(xué)（211）申請審核制博士招生- 電子、自動化、機(jī)械、生物學(xué)、物理相關(guān)專業(yè)
+1/25
-
【預(yù)碩】【預(yù)博】青島大學(xué)生物醫(yī)藥創(chuàng)新轉(zhuǎn)化
+1/16
-
【2026博士招生/博后招聘】北京航空航天大學(xué)潘彪課題組——AI芯片設(shè)計方向
+1/11
-
福建師范大學(xué)海峽柔性電子學(xué)院招收2026調(diào)劑碩士
+1/9
-
國家級人才團(tuán)隊誠聘生物信息相關(guān)專業(yè)的科研助理
+1/8
-
哈爾濱理工大學(xué)材料與化學(xué)工程學(xué)院 納米功能材料與光電催化團(tuán)隊 招收碩士研究生
+1/7
-
中國民航飛行學(xué)院先進(jìn)連接團(tuán)隊招收航空宇航、航空機(jī)械專業(yè)碩士研究生（學(xué)碩專碩均可）
+1/7
-
26年博士招生
+1/6
-
【第三輪招生，31號截止】澳科大誠招2026秋全獎博士研究生（納米醫(yī)學(xué)/生物材料方向）
+1/6
-
上海第二工業(yè)大學(xué)-朱大海課題組招生（過線就能調(diào)劑）
+1/5
-
125603工業(yè)工程與管理調(diào)劑
+1/2
-
歡迎報考浙江工業(yè)大學(xué)�？巳芈�(lián)合學(xué)院碩士研究生
+1/1
-
香港科技大學(xué)（廣州） 招收溫室氣體監(jiān)測與減排方向的博士、博后、碩士、RA
+1/1
-
北京大學(xué)環(huán)境學(xué)院張子帥研究員招收科研助理
+1/1
-
【2026年春季補(bǔ)錄】廈門大學(xué)田中群院士/易駿教授納米光子學(xué)團(tuán)隊招收2026年入學(xué)博士生
+1/1
zhangguangping
木蟲 (著名寫手)
- 1ST強(qiáng)帖: 23
- 應(yīng)助: 71 (初中生)
- 貴賓: 0.031
- 金幣: 1973.1
- 帖子: 2681
- 在線: 2128.8小時
- 蟲號: 529624
★ ★ ★ ★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
zzy870720z(金幣+5): 謝謝zhang兄指教，說的太好了 2011-03-08 15:56:38
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
zzy870720z(金幣+5): 謝謝zhang兄指教，說的太好了 2011-03-08 15:56:38
|
其實(shí)說句實(shí)話，我是真不厲害，很菜的。不是謙虛，我是自我表白。我在論壇里面比較活躍不是顯擺我知道的多，而是恰是我知道少，但是身邊有沒有人用這個軟件。所以只能依靠emuch論壇的這個版塊。通過討論可以相互交流一下，哪怕我的是錯誤的，我也可以說一下，這樣高手看見了，可以來指正我，這也是一個自我的進(jìn)步。另外咱們這個版塊不是很活躍。貌似高手現(xiàn)在都不大來了，或者來了之后看到問題的幼稚就嗤之一笑，然后不懶的回復(fù)了。希望老蟲還是多來指點(diǎn)指點(diǎn)我們這些幼蟲們。 另外你的這個問題實(shí)際上可以通過多認(rèn)真的讀一讀手冊是可以得到解決的。使用siesta工作的是手邊要么有一本打印的手冊。要么在電腦上打開一份電子的手冊。有什么問題，直接搜索一下相關(guān)的關(guān)鍵字。比如你的問題：直接在手冊里面用ctrl+F搜索ghost,立馬會定位到ChemicalSpeciesLabel位置，再搜，會定位到AtomicMass位置，接下來在PS.KBprojectors會發(fā)現(xiàn)對ghost設(shè)置的一些條件，這個需要檢查一下產(chǎn)生的PS.KBprojectors對ghost atom是不是只有一個，如果是的話，那就是程序給你處理好了。如果不是的話，計算出來預(yù)料不到結(jié)果，需要考慮這個問題。PAO.Basis在這兒，解釋了程序?qū)host atom的基矢的處理辦法。 到此為止，應(yīng)該對siesta處理ghost atom有一個了解。如果不是量化高手的話，可能對手冊上的解釋還是知道如何下手，這個時候需要google了。在google里面輸入ghost atom in siesta,一搜，google大多數(shù)時候不會辜負(fù)你的。 http://www.google.com.hk/search? ... i=&aql=&oq= 這個頭兩條就是siesta前輩們討論關(guān)于怎么設(shè)置ghost。 很懊惱的是打不開，這個時候我來告訴你，網(wǎng)頁代理來幫助：http://www.7daili.com/，首先打開這個網(wǎng)頁代理，然后打開你要打開但是打不的那個網(wǎng)頁，從瀏覽器上把這個網(wǎng)址復(fù)制到你的網(wǎng)頁代理里面，然后讓它給你打開。不要嫌棄這個太慢，耐心點(diǎn)，能打開總比打不開好多了。燈打開之后，看到里面討論的時候，欣喜不由從心底產(chǎn)生。在之內(nèi)的鏈接可以直接打開了（其實(shí)還是網(wǎng)頁代理給你打開的）。不管你從這些討論中得沒得到你要的結(jié)果，但是起碼，你學(xué)習(xí)了，搜索了，看到之前也有人討論這個問題了。這就夠了。如果問題解決了，皆大歡喜，如果還是沒有解決，在走下一條通往羅馬的大道�？傊�，這個過程是一個需要的過程，我認(rèn)為也是發(fā)現(xiàn)問題，解決問題的作為學(xué)生的過程�？赡苣銜f，作為學(xué)生，時間短，論文壓力大。沒有時間花費(fèi)在這個上面，直接問問知道就行。前提是如果有知道在身邊什么都好說。如果沒有呢？難道咱們的問題不不解決了？向小木蟲丟出來一個求助，不一定有人立馬回復(fù)你。這個期間還得自己去走上面的過程，只不過小木蟲的討論給上面的過程多了一個希望而已。 這些只是我的感想而已。我也是初學(xué)者。來到小木蟲，首先，要搜索一下關(guān)于你的這個話題的帖子。最起碼，來到小木蟲定居之后，要把之前討論過的有關(guān)你用的這個軟件的所有帖子全部看一遍，知道別人都討論了些什么。結(jié)果怎么樣，哪怕自己用不到，說不定以后用的時候，就需要回來再次檢索。平時沒事，一定要天天登錄一下小木蟲看一下自己用的這個程序有什么新的問題，別人是怎么解決的。這樣對該軟件在小木蟲的歷史就有了了解。對歷史不懂得人，會走很多彎路的。比如最近yashumimi發(fā)的那個帖子transista的錯誤，其實(shí)，別人之前早就討論過了。只是她沒有關(guān)心，以為自己不會碰到這個問題，壓根就沒有看這個帖子。但是她就應(yīng)為自己對這個帖子的默然，也像mail-list求助了，也發(fā)帖子了。我也看到一段時間沒有人回復(fù)。估計期間她也費(fèi)了好大勁，做了不少計算，并且花了不少機(jī)時。這個都是對歷史不解或者上小木蟲之后沒有首先檢索惹的禍。 最后，小木蟲不是一個有問題就往上丟，沒問題的時候壓根不來的地方，而是一個同行之間相互了解相互交流的平臺，窗口。就像一些學(xué)術(shù)會議一樣。只不過小木蟲上天天再開“學(xué)術(shù)會議”，“作報告的人”都是一些資源的提供者，更或者對我們這些板塊，“作報告的人”是一些遇到問題的人，報告的內(nèi)容就是問題本身而已。 當(dāng)自己不是“作報告的人”時候要常來當(dāng)聽報告的。因?yàn)檫@個問題被別人討論之后，再重復(fù)的提出來，可能會使得一些高手失去回答的耐心。不但對自己沒有好處，然高手失去了對咱們新一代新手的信心，對別的新手來說，是一個莫大的損失。 我也看到有些蟲友都工作了，都當(dāng)教授了，還是常常忙里偷閑來答復(fù)一下新蟲的問題。我看了，確實(shí)有些感動。 最后，祝愿小木蟲SIESTA板塊，越來越紅火，能成就一些新手。希望老蟲多來哺養(yǎng)一下我們這些新蟲。 末了，言辭不妥之處，望新老蟲友見諒。  |
4樓2011-03-08 11:32:09
zhangguangping
木蟲 (著名寫手)
- 1ST強(qiáng)帖: 23
- 應(yīng)助: 71 (初中生)
- 貴賓: 0.031
- 金幣: 1973.1
- 帖子: 2681
- 在線: 2128.8小時
- 蟲號: 529624
★ ★
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ù)用一個負(fù)值就行，然后最好給他起一個別的元素名稱，贗勢也相應(yīng)的改變一下。這個時候他的原子質(zhì)量就會是1.0d30 a.u.如果你計算的是金剛石，那么gost原子的贗勢就用C的就行，和你就算中用的其他的C用一個贗勢就行，只不過最好起一個別的名字，然后對贗勢也重命名一下，將其元素序號用一個負(fù)值比如-6表示。 |
2樓2011-03-07 22:14:00
|
你好，這是我測試NV在金剛石2*2*2超單胞C62N中的測試輸入數(shù)據(jù)： 計算鬼態(tài)為Cg，在搭建結(jié)構(gòu)的時候，我在空位處加上一個C原子，寫成輸入文件以后，將其改成Cg，標(biāo)號為3。其中，C標(biāo)號為1，N標(biāo)號為2。贗勢方面，我將C.psf直接復(fù)制成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
zhangguangping
木蟲 (著名寫手)
- 1ST強(qiáng)帖: 23
- 應(yīng)助: 71 (初中生)
- 貴賓: 0.031
- 金幣: 1973.1
- 帖子: 2681
- 在線: 2128.8小時
- 蟲號: 529624
★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
zzy870720z(金幣+2): 繼續(xù)獎勵 2011-03-08 15:56:51
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
zzy870720z(金幣+2): 繼續(xù)獎勵 2011-03-08 15:56:51
|
其實(shí)，我之前從來沒有做過ghost state的計算，并且對ghost state的知識和用處不了解。僅僅是我之前看手冊的時候記得手冊中提到了ghost atom的處理。但是具體怎么處理，我是沒往心里記的。昨天晚上臨回宿舍前，上了一下小木蟲，發(fā)現(xiàn)了你的這個問題。怕你急于解答，所以我就臨時的翻看了一下手冊，并且到網(wǎng)上搜索了一下ghos state在siesta中的討論而已。具體的加入ghost atom之后計算結(jié)果應(yīng)該是什么樣子是預(yù)期，我不清楚，這個應(yīng)該你是高手。 需要指明的是，C_g.psf里面的內(nèi)容不需要改的。表頭部分的，程序是不讀的，只是給使用者一個參數(shù)介紹，讓使用者知道這個贗勢的類型和參數(shù)而已。 |
5樓2011-03-08 11:39:57
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 268求調(diào)劑 +7 | 好運(yùn)連綿不絕 2026-03-12 | 8/400 |
|
|---|---|---|---|---|
|
|
[考研] 085601材料工程專碩求調(diào)劑 +4 | 慕寒mio 2026-03-16 | 4/200 |
|
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
|
[基金申請] 被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了 +3 | beefly 2026-03-14 | 3/150 |
|
|
|
[考研] 085600材料與化工求調(diào)劑 +5 | 緒幸與子 2026-03-17 | 5/250 |
|
|
[考研] 289求調(diào)劑 +6 | 步川酷紫123 2026-03-11 | 6/300 |
|
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
|
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
|
[教師之家] 焦慮 +7 | 水冰月月野兔 2026-03-13 | 9/450 |
|
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請 +6 | zlingli 2026-03-13 | 6/300 |
|
|
|
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
|
[考研] 288求調(diào)劑 +4 | 奇點(diǎn)0314 2026-03-14 | 4/200 |
|
|
|
[考研] 328求調(diào)劑 +3 | 5201314Lsy！ 2026-03-13 | 6/300 |
|
|
|
[考研] 【0703化學(xué)調(diào)劑】-一志愿華中師范大學(xué)-六級475 +5 | Becho359 2026-03-11 | 5/250 |
|
|
|
[考研] 0703化學(xué)一志愿211 總分320求調(diào)劑 +5 | 瑪卡巴卡啊哈 2026-03-11 | 5/250 |
|
|
|
[考研] 301求調(diào)劑 +6 | Liyouyumairs 2026-03-11 | 6/300 |
|
|
|
[考研] 270求調(diào)劑 085600材料與化工專碩 +3 | YXCT 2026-03-11 | 3/150 |
|