 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3168 | 回復(fù): 16 | ||||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||||
hw45888792金蟲 (正式寫手)
|
[交流]
【求助】如何在siesta中設(shè)置ghost state原子的計(jì)算
|
|||||
|
用siesta計(jì)算體相金剛石中一個(gè)空位，形成能好像不準(zhǔn)，單胞越大，形成能也越大，沒有合理的解釋。 師兄說，siesta用原子軌道基組，可能存在鬼態(tài)原子，需要在空位處設(shè)置一個(gè)假想的原子，我查了一下手冊(cè)，沒有完全搞清楚如何設(shè)置？ 不知道哪位大俠有做過？ |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 siesta問題 | SIESTA | siesta基本概念 |
» 猜你喜歡
303求調(diào)劑
已經(jīng)有3人回復(fù)
-
有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～
已經(jīng)有7人回復(fù)
-
梁成偉老師課題組歡迎你的加入
已經(jīng)有10人回復(fù)
-
一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng)
已經(jīng)有8人回復(fù)
-
070300化學(xué)319求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有4人回復(fù)
-
268求調(diào)劑
已經(jīng)有7人回復(fù)
-
290求調(diào)劑
已經(jīng)有11人回復(fù)
-
材料與化工求調(diào)劑
已經(jīng)有4人回復(fù)
-
293求調(diào)劑
已經(jīng)有11人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 為什么采用CASTEP計(jì)算CaTiO3的光學(xué)性質(zhì)其中的靜態(tài)介電常數(shù)與實(shí)驗(yàn)值相差很大
已經(jīng)有14人回復(fù)
- 關(guān)于Discover動(dòng)力學(xué)計(jì)算的一個(gè)問題
已經(jīng)有26人回復(fù)
- 如何計(jì)算空位部分所帶的電子
已經(jīng)有5人回復(fù)
- 動(dòng)力學(xué)計(jì)算過程：時(shí)間步長(zhǎng)、步數(shù)一般設(shè)置為多少？
已經(jīng)有10人回復(fù)
- siesta 計(jì)算電子占據(jù)數(shù)的問題
已經(jīng)有3人回復(fù)
- 關(guān)于SIESTA中如何計(jì)算功函數(shù)的一種做法探討
已經(jīng)有8人回復(fù)
- siesta計(jì)算弱相互作用VDW
已經(jīng)有17人回復(fù)
- 【求助】siesta進(jìn)行ordern計(jì)算出錯(cuò)求助
已經(jīng)有4人回復(fù)
- 【求助】siesta關(guān)于金屬計(jì)算的參數(shù)設(shè)置
已經(jīng)有5人回復(fù)
- 【求助】Ms計(jì)算磁性要不要設(shè)置每個(gè)原子的磁矩？？
已經(jīng)有10人回復(fù)
- 【求助】siesta 計(jì)算DOS應(yīng)該注意什么問題
已經(jīng)有3人回復(fù)
» 搶金幣啦！回帖就可以得到:
-
南京大學(xué)智能驅(qū)動(dòng)與感知材料實(shí)驗(yàn)室誠(chéng)招博士/碩士/科研助理/訪問學(xué)生
+1/280
-
深圳技術(shù)大學(xué)集成電路與光電芯片學(xué)院郝俊杰課題組誠(chéng)聘博士后、研究助理、訪問學(xué)生
+2/128
-
南京醫(yī)科大學(xué)生殖醫(yī)學(xué)與子代健康國(guó)重實(shí)驗(yàn)室-董飛宏課題組-學(xué)術(shù)碩士招生
+1/97
-
物理學(xué) 調(diào)劑
+1/88
-
愿遇見你，遇見不止擦肩而過的緣分
+1/79
-
華南師范大學(xué)（211）博士招生- 電子、自動(dòng)化、機(jī)械、生物學(xué)、物理相關(guān)專業(yè)
+2/42
-
催化方向推薦一個(gè)極好的博導(dǎo)
+1/39
-
華南師范大學(xué)（211）博士招生- 電子、自動(dòng)化、機(jī)械、生物學(xué)、物理相關(guān)專業(yè)
+2/38
-
生殖醫(yī)學(xué)與子代健康全國(guó)重點(diǎn)實(shí)驗(yàn)室遺傳學(xué)課題組招收研究生(長(zhǎng)期有效)
+1/34
-
教研論文SCI期刊投稿選刊
+1/34
-
【武漢高�！俊拘履茉床牧吓c器件課題組】 招收優(yōu)秀碩士調(diào)劑生！�。�
+5/20
-
學(xué)術(shù)/專業(yè)碩士調(diào)劑——化學(xué)、資源與環(huán)境
+1/15
-
南京林業(yè)大學(xué)化學(xué)工程學(xué)院柏惺峰課題組招收2026年申請(qǐng)考核制博士2名
+1/12
-
重慶大學(xué)誠(chéng)招2026年生物材料方向博士生
+1/8
-
中國(guó)民航飛行學(xué)院先進(jìn)連接團(tuán)隊(duì)招收航空宇航、航空機(jī)械專業(yè)碩士研究生（學(xué)碩專碩均可）
+1/7
-
課題組招收環(huán)境及相關(guān)專業(yè)調(diào)劑碩士研究生（歡迎優(yōu)秀學(xué)生加入）
+1/6
-
上海第二工業(yè)大學(xué)-朱大海課題組招生（過線就能調(diào)劑）
+1/5
-
北京工業(yè)大學(xué)材料工程調(diào)劑03
+1/3
-
上海大學(xué)微電子學(xué)院李夢(mèng)姣教授團(tuán)隊(duì)-招收集成電路/電子/材料/微納制造博士生1名
+1/2
-
北京大學(xué)環(huán)境學(xué)院張子帥研究員招收科研助理
+1/1
hw45888792
金蟲 (正式寫手)
- 1ST強(qiáng)帖: 1
- 應(yīng)助: 5 (幼兒園)
- 金幣: 848.5
- 帖子: 339
- 在線: 282.9小時(shí)
- 蟲號(hào): 711113
7樓2011-03-08 13:21:23
zhangguangping
木蟲 (著名寫手)
- 1ST強(qiáng)帖: 23
- 應(yīng)助: 71 (初中生)
- 貴賓: 0.031
- 金幣: 1973.1
- 帖子: 2681
- 在線: 2128.8小時(shí)
- 蟲號(hào): 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
|
參照手冊(cè)上的ChemicalSpeciesLabel的說明，對(duì)于gost原子的原子序數(shù)用一個(gè)負(fù)值就行，然后最好給他起一個(gè)別的元素名稱，贗勢(shì)也相應(yīng)的改變一下。這個(gè)時(shí)候他的原子質(zhì)量就會(huì)是1.0d30 a.u.如果你計(jì)算的是金剛石，那么gost原子的贗勢(shì)就用C的就行，和你就算中用的其他的C用一個(gè)贗勢(shì)就行，只不過最好起一個(gè)別的名字，然后對(duì)贗勢(shì)也重命名一下，將其元素序號(hào)用一個(gè)負(fù)值比如-6表示。 |
2樓2011-03-07 22:14:00
hw45888792
金蟲 (正式寫手)
- 1ST強(qiáng)帖: 1
- 應(yīng)助: 5 (幼兒園)
- 金幣: 848.5
- 帖子: 339
- 在線: 282.9小時(shí)
- 蟲號(hào): 711113
|
你好，這是我測(cè)試NV在金剛石2*2*2超單胞C62N中的測(cè)試輸入數(shù)據(jù)： 計(jì)算鬼態(tài)為Cg，在搭建結(jié)構(gòu)的時(shí)候，我在空位處加上一個(gè)C原子，寫成輸入文件以后，將其改成Cg，標(biāo)號(hào)為3。其中，C標(biāo)號(hào)為1，N標(biāo)號(hào)為2。贗勢(shì)方面，我將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ī)計(jì)算，就不用鬼態(tài)，輸入文件如下： SystemName C62N SystemLabel C62N NumberOfSpecies 2 %block ChemicalSpeciesLabel 1 6 C 2 7 N %endblock ChemicalSpeciesLabel 不知道是否正確？ 計(jì)算過程沒有任何問題，對(duì)比結(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 計(jì)算結(jié)果顯示，加入鬼態(tài)后，總能量略低0.1eV左右，而且鬼態(tài)原子附近有大約0.2單自旋電荷分布，其它C原子是2個(gè)電荷分布。 |
3樓2011-03-08 11:11:00
zhangguangping
木蟲 (著名寫手)
- 1ST強(qiáng)帖: 23
- 應(yīng)助: 71 (初中生)
- 貴賓: 0.031
- 金幣: 1973.1
- 帖子: 2681
- 在線: 2128.8小時(shí)
- 蟲號(hào): 529624
★ ★ ★ ★ ★ ★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
zzy870720z(金幣+5): 謝謝zhang兄指教，說的太好了 2011-03-08 15:56:38
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
zzy870720z(金幣+5): 謝謝zhang兄指教，說的太好了 2011-03-08 15:56:38
|
其實(shí)說句實(shí)話，我是真不厲害，很菜的。不是謙虛，我是自我表白。我在論壇里面比較活躍不是顯擺我知道的多，而是恰是我知道少，但是身邊有沒有人用這個(gè)軟件。所以只能依靠emuch論壇的這個(gè)版塊。通過討論可以相互交流一下，哪怕我的是錯(cuò)誤的，我也可以說一下，這樣高手看見了，可以來指正我，這也是一個(gè)自我的進(jìn)步。另外咱們這個(gè)版塊不是很活躍。貌似高手現(xiàn)在都不大來了，或者來了之后看到問題的幼稚就嗤之一笑，然后不懶的回復(fù)了。希望老蟲還是多來指點(diǎn)指點(diǎn)我們這些幼蟲們。 另外你的這個(gè)問題實(shí)際上可以通過多認(rèn)真的讀一讀手冊(cè)是可以得到解決的。使用siesta工作的是手邊要么有一本打印的手冊(cè)。要么在電腦上打開一份電子的手冊(cè)。有什么問題，直接搜索一下相關(guān)的關(guān)鍵字。比如你的問題：直接在手冊(cè)里面用ctrl+F搜索ghost,立馬會(huì)定位到ChemicalSpeciesLabel位置，再搜，會(huì)定位到AtomicMass位置，接下來在PS.KBprojectors會(huì)發(fā)現(xiàn)對(duì)ghost設(shè)置的一些條件，這個(gè)需要檢查一下產(chǎn)生的PS.KBprojectors對(duì)ghost atom是不是只有一個(gè)，如果是的話，那就是程序給你處理好了。如果不是的話，計(jì)算出來預(yù)料不到結(jié)果，需要考慮這個(gè)問題。PAO.Basis在這兒，解釋了程序?qū)host atom的基矢的處理辦法。 到此為止，應(yīng)該對(duì)siesta處理ghost atom有一個(gè)了解。如果不是量化高手的話，可能對(duì)手冊(cè)上的解釋還是知道如何下手，這個(gè)時(shí)候需要google了。在google里面輸入ghost atom in siesta,一搜，google大多數(shù)時(shí)候不會(huì)辜負(fù)你的。 http://www.google.com.hk/search? ... i=&aql=&oq= 這個(gè)頭兩條就是siesta前輩們討論關(guān)于怎么設(shè)置ghost。 很懊惱的是打不開，這個(gè)時(shí)候我來告訴你，網(wǎng)頁代理來幫助：http://www.7daili.com/，首先打開這個(gè)網(wǎng)頁代理，然后打開你要打開但是打不的那個(gè)網(wǎng)頁，從瀏覽器上把這個(gè)網(wǎng)址復(fù)制到你的網(wǎng)頁代理里面，然后讓它給你打開。不要嫌棄這個(gè)太慢，耐心點(diǎn)，能打開總比打不開好多了。燈打開之后，看到里面討論的時(shí)候，欣喜不由從心底產(chǎn)生。在之內(nèi)的鏈接可以直接打開了（其實(shí)還是網(wǎng)頁代理給你打開的）。不管你從這些討論中得沒得到你要的結(jié)果，但是起碼，你學(xué)習(xí)了，搜索了，看到之前也有人討論這個(gè)問題了。這就夠了。如果問題解決了，皆大歡喜，如果還是沒有解決，在走下一條通往羅馬的大道�？傊�，這個(gè)過程是一個(gè)需要的過程，我認(rèn)為也是發(fā)現(xiàn)問題，解決問題的作為學(xué)生的過程。可能你會(huì)說，作為學(xué)生，時(shí)間短，論文壓力大。沒有時(shí)間花費(fèi)在這個(gè)上面，直接問問知道就行。前提是如果有知道在身邊什么都好說。如果沒有呢？難道咱們的問題不不解決了？向小木蟲丟出來一個(gè)求助，不一定有人立馬回復(fù)你。這個(gè)期間還得自己去走上面的過程，只不過小木蟲的討論給上面的過程多了一個(gè)希望而已。 這些只是我的感想而已。我也是初學(xué)者。來到小木蟲，首先，要搜索一下關(guān)于你的這個(gè)話題的帖子。最起碼，來到小木蟲定居之后，要把之前討論過的有關(guān)你用的這個(gè)軟件的所有帖子全部看一遍，知道別人都討論了些什么。結(jié)果怎么樣，哪怕自己用不到，說不定以后用的時(shí)候，就需要回來再次檢索。平時(shí)沒事，一定要天天登錄一下小木蟲看一下自己用的這個(gè)程序有什么新的問題，別人是怎么解決的。這樣對(duì)該軟件在小木蟲的歷史就有了了解。對(duì)歷史不懂得人，會(huì)走很多彎路的。比如最近yashumimi發(fā)的那個(gè)帖子transista的錯(cuò)誤，其實(shí)，別人之前早就討論過了。只是她沒有關(guān)心，以為自己不會(huì)碰到這個(gè)問題，壓根就沒有看這個(gè)帖子。但是她就應(yīng)為自己對(duì)這個(gè)帖子的默然，也像mail-list求助了，也發(fā)帖子了。我也看到一段時(shí)間沒有人回復(fù)。估計(jì)期間她也費(fèi)了好大勁，做了不少計(jì)算，并且花了不少機(jī)時(shí)。這個(gè)都是對(duì)歷史不解或者上小木蟲之后沒有首先檢索惹的禍。 最后，小木蟲不是一個(gè)有問題就往上丟，沒問題的時(shí)候壓根不來的地方，而是一個(gè)同行之間相互了解相互交流的平臺(tái)，窗口。就像一些學(xué)術(shù)會(huì)議一樣。只不過小木蟲上天天再開“學(xué)術(shù)會(huì)議”，“作報(bào)告的人”都是一些資源的提供者，更或者對(duì)我們這些板塊，“作報(bào)告的人”是一些遇到問題的人，報(bào)告的內(nèi)容就是問題本身而已。 當(dāng)自己不是“作報(bào)告的人”時(shí)候要常來當(dāng)聽報(bào)告的。因?yàn)檫@個(gè)問題被別人討論之后，再重復(fù)的提出來，可能會(huì)使得一些高手失去回答的耐心。不但對(duì)自己沒有好處，然高手失去了對(duì)咱們新一代新手的信心，對(duì)別的新手來說，是一個(gè)莫大的損失。 我也看到有些蟲友都工作了，都當(dāng)教授了，還是常常忙里偷閑來答復(fù)一下新蟲的問題。我看了，確實(shí)有些感動(dòng)。 最后，祝愿小木蟲SIESTA板塊，越來越紅火，能成就一些新手。希望老蟲多來哺養(yǎng)一下我們這些新蟲。 末了，言辭不妥之處，望新老蟲友見諒。  |
4樓2011-03-08 11:32:09
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[碩博家園] 深圳大學(xué)碩士招生（2026秋，傳感器方向，僅錄取第一志愿） +4 | xujiaoszu 2026-03-11 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +6 | SUICHILD 2026-03-12 | 6/300 |
|
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
|
[考研] 085601求調(diào)劑 +3 | Du.11 2026-03-16 | 3/150 |
|
|
|
[考研] 311求調(diào)劑 +5 | 26研0 2026-03-15 | 5/250 |
|
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +5 | d如愿上岸 2026-03-12 | 8/400 |
|
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-12 | 6/300 |
|
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
|
[考研] 080500，材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識(shí)可樂 2026-03-14 | 5/250 |
|
|
|
[基金申請(qǐng)] 現(xiàn)在如何回避去年的某一個(gè)專家，不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
|
[考研] 0805，333求調(diào)劑 +3 | 112253525 2026-03-10 | 3/150 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-12 | 7/350 |
|
|
|
[考研] 0703，333分求調(diào)劑 一志愿鄭州大學(xué)-物理化學(xué) +3 | 李魔女斗篷 2026-03-11 | 3/150 |
|
|
|
[考研] ［0860］321分求調(diào)劑，ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|
|
|
[考研] 0856材料與化工301求調(diào)劑 +5 | 奕束光 2026-03-13 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +3 | 程雨杭 2026-03-12 | 3/150 |
|
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|
|
|
[考博] 26讀博 +4 | Rui135246 2026-03-12 | 10/500 |
|