 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 1354 | 回復(fù): 3 | ||
[求助]
siesta 計(jì)算電子占據(jù)數(shù)的問題
|
|
我用siesta計(jì)算石墨烯納米帶的H和C原子的占據(jù)數(shù)發(fā)現(xiàn)，用SZ基組算出來的和用DZP計(jì)算出來的結(jié)果有很大的區(qū)別，DZP計(jì)算的結(jié)果很不合理。接著我用水分子做了一下測試，結(jié)果如下： SZ:Species: H Atom Qatom Qorb 1s 1 0.840 0.840 2 0.840 0.840 3 0.840 0.840 4 0.840 0.840 Species: C Atom Qatom Qorb 2s 2py 2pz 2px 5 4.639 1.216 1.141 1.141 1.141 DZP： Atom Qatom Qorb 1s 1s 1Ppy 1Ppz 1Ppx 1 1.203 0.831 0.272 0.034 0.033 0.033 2 1.203 0.831 0.272 0.033 0.034 0.033 3 1.204 0.833 0.271 0.033 0.033 0.033 4 1.204 0.833 0.271 0.033 0.033 0.034 2s 2s 2py 2pz 2px 2py 2pz 2px 2Pdxy 2Pdyz 2Pdz2 2Pdxz 2Pdx2-y2 5 3.187 -0.047 0.750 0.194 0.194 0.192 0.622 0.622 0.622 0.004 0.008 0.009 0.004 0.011 很明顯DZP 的結(jié)果不對，是不是siesta程序本身的原因？ |
回復(fù)此樓» 猜你喜歡
材料調(diào)劑
已經(jīng)有3人回復(fù)
-
考研調(diào)劑
已經(jīng)有5人回復(fù)
-
26申博
已經(jīng)有4人回復(fù)
-
調(diào)劑申請
已經(jīng)有7人回復(fù)
-
環(huán)境工程297分求調(diào)劑一志愿杭高院
已經(jīng)有8人回復(fù)
-
070300化學(xué)調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿085600中科院寧波所276分求調(diào)劑
已經(jīng)有15人回復(fù)
-
288資源與環(huán)境專碩求調(diào)劑，不限專業(yè)，有學(xué)上就行
已經(jīng)有22人回復(fù)
-
273求調(diào)劑
已經(jīng)有8人回復(fù)
-
267求調(diào)劑
已經(jīng)有9人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 用siesta畫聲子譜的問題
已經(jīng)有12人回復(fù)
- siesta計(jì)算dipole moment（電偶極矩）
已經(jīng)有7人回復(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 計(jì)算DOS應(yīng)該注意什么問題
已經(jīng)有3人回復(fù)
- 【求助】用siesta計(jì)算出現(xiàn)的問題
已經(jīng)有3人回復(fù)
- 【求助】siesta中關(guān)于幾何結(jié)構(gòu)優(yōu)化的問題【完結(jié)】
已經(jīng)有17人回復(fù)
- 【求助】關(guān)于siesta中贗勢價(jià)電子軌道的問題
已經(jīng)有10人回復(fù)
|
我的甲烷的input文件為： SystemName ch4 SystemLabel Ch4 NumberOfAtoms 5 NumberOfSpecies 2 MeshCutoff 50 Ry %block ChemicalSpeciesLabel 1 1 H # Species index, atomic number, species label 2 6 C %endblock ChemicalSpeciesLabel PAO.BasisSize DZP XC.functional GGA XC.authors PBE WriteMullikenPop 1 AtomicCoordinatesFormat Ang %block AtomicCoordinatesAndAtomicSpecies -0.380000123 -1.038179245 0.278179290 1 -0.379999607 0.278179338 -1.038179422 1 -0.379999607 0.760000084 0.760000084 1 1.139999981 0.000000000 0.000000000 1 0.000000000 0.000000000 0.000000000 2 %endblock AtomicCoordinatesAndAtomicSpecies save-rho T save-delta-rho T save-total-potential T save-neutral-atom-potential T save-hs T |
2樓2011-12-08 10:25:25
3樓2011-12-09 14:31:06
zhangguangping
木蟲 (著名寫手)
- 1ST強(qiáng)帖: 23
- 應(yīng)助: 71 (初中生)
- 貴賓: 0.031
- 金幣: 1973.1
- 散金: 9059
- 紅花: 79
- 帖子: 2681
- 在線: 2128.8小時(shí)
- 蟲號(hào): 529624
- 注冊: 2008-03-20
- 性別: GG
- 專業(yè): 原子和分子物理
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
凍凍2007: 金幣+5, ★★★很有幫助 2012-05-04 15:52:07
凍凍2007: 金幣+5 2012-05-15 22:28:57
感謝參與，應(yīng)助指數(shù) +1
凍凍2007: 金幣+5, ★★★很有幫助 2012-05-04 15:52:07
凍凍2007: 金幣+5 2012-05-15 22:28:57
|
首先，需要聲明的是對于采用原子軌道基組計(jì)算，隨著基組增大，計(jì)算結(jié)果并不一定要更好，但是趨勢是好的，這一點(diǎn)和平面波基組是不同的。 第二，關(guān)于你的計(jì)算的問題。存在至少兩個(gè)疑問。一個(gè)是你的MeshCutoff 50 Ry我認(rèn)為是小的，建議使用150Ry甚至200Ry。另外，一個(gè)問題是你的基組的問題，建議加入這個(gè)選項(xiàng)：PAO.EnergyShift 100 meV，并且使用我的這個(gè)值，默認(rèn)是272 meV。 估計(jì)這兩個(gè)選項(xiàng)可以改善你的問題。 再有就是你的贗勢的問題�？梢宰詈罂紤]它的因素。你先用官方網(wǎng)站上的贗勢計(jì)算試試。 |
4樓2012-02-28 15:03:10
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 凍凍2007 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考博] 26申博 +4 | 加油沖啊！ 2026-03-26 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 調(diào)劑推薦 +7 | 清酒714 2026-03-26 | 8/400 |
|
|
[考研] 086000生物與醫(yī)藥298調(diào)劑求助 +3 | 元元青青 2026-03-31 | 5/250 |
|
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-29 | 10/500 |
|
|
|
[考研] 309求調(diào)劑 +17 | 誰不是少年 2026-03-29 | 17/850 |
|
|
|
[考研] 311求調(diào)劑一志愿合肥工業(yè)大學(xué) +11 | 秋二十二 2026-03-30 | 11/550 |
|
|
|
[考研] 材料專碩調(diào)劑 +13 | 椰椰。 2026-03-29 | 13/650 |
|
|
|
[考研]
 一志愿武漢理工，總分321，英一數(shù)二，求老師收留。
+12
|
nnnnnnn5 2026-03-25 | 12/600 |
|
|
|
[考研] 085601一志愿西北工業(yè)大學(xué)初試346 +4 | 085601初試346 2026-03-30 | 4/200 |
|
|
|
[考研] 370求調(diào)劑 +3 | 080700調(diào)劑 2026-03-30 | 3/150 |
|
|
|
[考研] 本科211總分289，08工學(xué)真心求調(diào)劑 +3 | utopiaE 2026-03-30 | 3/150 |
|
|
|
[考研] 總分293求調(diào)劑 +8 | 加一一九 2026-03-25 | 11/550 |
|
|
|
[考研] 343求調(diào)劑 +6 | 愛羈絆 2026-03-29 | 6/300 |
|
|
|
[考研] 394求調(diào)劑 +3 | 好事多磨靜候佳?/a> 2026-03-26 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +6 | 蘆lty 2026-03-25 | 7/350 |
|
|
|
[考研] 070300求調(diào)劑306分 +4 | 26要上岸 2026-03-27 | 4/200 |
|
|
|
[考研] 藥學(xué)105500求調(diào)劑 +3 | Ssun。。 2026-03-28 | 3/150 |
|
|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|
|
[考研] 復(fù)試調(diào)劑，一志愿南農(nóng)083200食品科學(xué)與工程 +5 | XQTJZ 2026-03-26 | 5/250 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|