當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

identation

金蟲 (著名寫手)

[交流] 自旋體系的聲子計(jì)算應(yīng)該用什么軟件呢？

我接觸過的幾款軟件的體會如下：
PWSCF
pwscf對待自旋體系有兩種方法：1，采用starting_magnetization指定每個(gè)原子的初始自旋，通過scf計(jì)算得到體系總磁矩，這種方法把體系當(dāng)成是metal，必須指定smearing以及degauss；然而，指定了degauss的體系計(jì)算Gamma點(diǎn)聲子的時(shí)候就不能計(jì)算Born effective charge了dielectric tensor了，很多時(shí)候我想關(guān)注的就是有效電荷；當(dāng)然，對于磁性體系，按照PWSCF的說明，是可以通過Berry phase方法先計(jì)算介電極化，通過介電極化對原子位移求導(dǎo)來得到有效電荷的，而且pwscf中只要degauss<0.01時(shí)，是可以計(jì)算介電極化的，但是PWSCF有一額外要求，就是nbnd=nelec/2，也就是帶數(shù)目是電子數(shù)目的一半，對于指定自旋的體系，因?yàn)閜wscf視之為金屬，所以還是不能計(jì)算介電極化；2，pwscf采用fixed spin moment的方法處理自旋體系，此時(shí)需要指定tot_magnetization，即指定總自旋，此時(shí)pwscf是按絕緣體處理待計(jì)算體系的，但是計(jì)算聲子的時(shí)候又難以收斂，只要加上Born有效電荷計(jì)算選項(xiàng)，計(jì)算就不能收斂了；而且Gamma點(diǎn)的計(jì)算結(jié)果采用指定tot_magnetization的結(jié)果和采用tarting_magnetization的結(jié)果又不一致，我之前試過的一個(gè)體系是一個(gè)有虛頻一個(gè)沒有，因此很困惑。
ABINIT
ABINIT計(jì)算聲子譜確實(shí)沒有PWSCF的這些問題，但是abinit的磁性計(jì)算很郁悶，我不知道是我輸入文件設(shè)定錯(cuò)了還是什么原因（我在版上問了很多次），當(dāng)采用絕緣體的形式，即occopt=0時(shí)，計(jì)算得到體系的磁矩總是幾乎為零；和實(shí)際根本不符，而當(dāng)我使用occopt=3時(shí)，計(jì)算得到的結(jié)果又是沒有帶隙，我一直很困惑，究竟是怎么回事呢？
我卻有看到很多文章使用ABINIT和PWSCF計(jì)算體系介電極化和聲子的，而且是自旋體系，究竟是什么原因我計(jì)算的結(jié)果不理想呢？
有沒有開源的活免費(fèi)的聲子計(jì)算軟件推薦呢？求大家賜教！

回復(fù)此樓

» 猜你喜歡

265求調(diào)劑已經(jīng)有8人回復(fù)
311求調(diào)劑已經(jīng)有4人回復(fù)
一志愿廈門大學(xué)化學(xué)學(xué)碩307求調(diào)劑已經(jīng)有13人回復(fù)
一志愿211，335分，0856，求調(diào)劑院校和導(dǎo)師已經(jīng)有15人回復(fù)
275學(xué)碩081000服從調(diào)劑到其他專業(yè)，保不住本專業(yè)了已經(jīng)有4人回復(fù)
274求調(diào)劑已經(jīng)有8人回復(fù)
367求調(diào)劑已經(jīng)有8人回復(fù)
一志愿9初試366 本雙非求調(diào)劑已經(jīng)有3人回復(fù)
314求調(diào)劑已經(jīng)有9人回復(fù)
298求調(diào)劑已經(jīng)有3人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:

VASP接FROPHO計(jì)算H2時(shí)聲子譜有虛頻怎么辦已經(jīng)有20人回復(fù)
今天用CASTEP計(jì)算聲子譜，出現(xiàn)下列提示，有能看懂的大俠幫忙解釋下已經(jīng)有11人回復(fù)
聲子的自旋已經(jīng)有41人回復(fù)
為何用晶格常數(shù)的實(shí)驗(yàn)值計(jì)算聲子？已經(jīng)有18人回復(fù)
關(guān)于用超胞方法計(jì)算聲子的一些經(jīng)驗(yàn) 已經(jīng)有21人回復(fù)
MS計(jì)算晶體聲子譜的時(shí)候，橫坐標(biāo)有負(fù)能量值，為什么？已經(jīng)有11人回復(fù)
castep計(jì)算Fe2O3聲子譜失敗，原因求教~~ 已經(jīng)有10人回復(fù)
【求助】有人用VASP5.2計(jì)算過聲子譜嗎？已經(jīng)有8人回復(fù)
計(jì)算Sb2Te3的聲子色散曲線，請教用什么軟件合適,VASP+phonopy怎么樣已經(jīng)有7人回復(fù)
【求助】關(guān)于PWSCF介電常數(shù)和聲子譜的計(jì)算已經(jīng)有7人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2011-07-03 21:29:14

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

leimiao_hit

木蟲之王 (文學(xué)泰斗)

★ ★
identation(金幣+3):謝謝參與
franch(金幣+1): 謝謝回帖交流 2011-07-04 11:03:06
identation(金幣+5): 2011-07-04 13:31:00

用elk或者exciting

贊一下(2人)

回復(fù)此樓

3樓2011-07-04 00:30:34

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 6 個(gè)回答

iamikaruk

木蟲 (著名寫手)

★ ★ ★
identation(金幣+3):謝謝參與
franch(金幣+2): 謝謝回帖交流 2011-07-04 11:02:59
identation(金幣+5): 2011-07-04 13:30:55

沒有用abinit算過自旋體系，但是絕緣體的確得用occopt<3的選項(xiàng)，occopt>=3是用于金屬的。
你用的贗勢是什么？如果你計(jì)算資源夠的話可以試試用elk或者exciting

贊一下(2人)

回復(fù)此樓

2樓2011-07-04 00:12:40

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

identation

金蟲 (著名寫手)

我發(fā)現(xiàn)，abinit和pwscf在處理磁性體系的時(shí)候總是和金屬很“曖昧”，所以磁性體系的處理感覺都很麻煩，elk用起來確實(shí)不錯(cuò)，磁性并不與金屬和絕緣體有什么聯(lián)系，可是elk計(jì)算太慢了

贊一下

回復(fù)此樓

5樓2011-07-04 13:33:23

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 296材料專碩求調(diào)劑 +3	202451007219 2026-04-02	4/200	2026-04-02 09:11 by syjjj0321
[考研] 301求調(diào)劑 +13	A_JiXing 2026-04-01	13/650	2026-04-02 09:01 by sanrepian
[考研] 266求調(diào)劑 +3	學(xué)員97LZgn 2026-04-02	3/150	2026-04-02 08:54 by Jaylen.
[考研] 085601一志愿中山大學(xué)深圳材料工程330求調(diào)劑 +7	pipiver 2026-03-30	7/350	2026-04-02 07:58 by 尚水閣主
[考研] 生物學(xué)308分求調(diào)劑（一志愿華東師大） +4	相信必會光芒萬?/a> 2026-03-31	4/200	2026-04-02 06:52 by ilovexiaobin
[考研] 085600 一志愿9 總分351 求調(diào)劑學(xué)校 +7	czhcz 2026-03-31	9/450	2026-04-01 19:24 by 無際的草原
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +15	柳溪邊 2026-03-26	16/800	2026-04-01 17:07 by 七度不信任
[考研] 350求調(diào)劑 +7	阿佳～ 2026-03-31	7/350	2026-04-01 16:12 by yanflower7133
[考研] 調(diào)劑 +3	好好讀書。 2026-04-01	5/250	2026-04-01 15:39 by 好好讀書。
[考研] 調(diào)劑 +5	好好讀書。 2026-03-28	7/350	2026-04-01 15:32 by 王亮_大連醫(yī)科大
[考研] 0856初試324分求調(diào)劑 +5	想上學(xué)求調(diào) 2026-04-01	5/250	2026-04-01 15:28 by 考研上岸呀一定?/a>
[考研] 一志愿鄭大085600，310分求調(diào)劑 +6	李瀟可 2026-03-26	6/300	2026-04-01 14:44 by chenqifeng666
[考研] 求調(diào)劑，一志愿南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +10	@taotao 2026-03-29	10/500	2026-04-01 09:30 by oooqiao
[考研] 土木304求調(diào)劑 +5	頂級擦擦 2026-03-31	5/250	2026-04-01 08:15 by fdcxdystjk￥
[考研] 309分085801求調(diào)劑 +7	學(xué)員Gtwj7W 2026-03-31	7/350	2026-04-01 02:36 by BruceLiu320
[考研] 合肥區(qū)域性重點(diǎn)一本招收調(diào)劑 +4	6266jl 2026-03-30	8/400	2026-03-31 18:43 by 6266jl
[考研] 一志愿大連理工大學(xué)，機(jī)械工程學(xué)碩，341 +3	西瓜田的守望者 2026-03-30	3/150	2026-03-31 11:08 by asdfzly
[考研] 083000環(huán)境科學(xué)與工程調(diào)劑，總分281 +4	橙子（勝意） 2026-03-30	4/200	2026-03-31 00:44 by Linzejun
[考研] 305求調(diào)劑 +5	哇盧卡庫 2026-03-26	5/250	2026-03-27 14:01 by laoshidan
[考研] 292求調(diào)劑 +4	求求了收下我吧?/a> 2026-03-26	4/200	2026-03-27 10:37 by zhshch