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

fzuvivi

金蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 1168
帖子: 224
在線: 330.1小時(shí)
蟲號(hào): 881678

[交流] 【求助】關(guān)于擬合得到體彈性模量&自旋極化&LDA+U等問(wèn)題，望大家指點(diǎn)指點(diǎn)

最近跑vasp，遇到了一些問(wèn)題，

（1）關(guān)于Si 體彈性模量的計(jì)算：

   根據(jù)體模量的定義，我通過(guò)單體積下充分弛豫后得到一系列的不同晶格常數(shù)和對(duì)應(yīng)能量值的關(guān)系（E-A），然后通過(guò)一個(gè)腳本murn.x利用Murnaghan方程來(lái)擬合得到最低能量對(duì)應(yīng)的晶格常數(shù)大�。╝lat）以及體模量值（B0）；可是發(fā)現(xiàn)，雖然晶格常數(shù)可以與實(shí)驗(yàn)值吻合的很好（誤差在1%左右）；但體模量卻出現(xiàn)較的偏差（誤差達(dá)10%以上）。后來(lái)也通過(guò)計(jì)算彈性常數(shù)來(lái)得到體模量，可是準(zhǔn)確性沒(méi)什么改觀。。。請(qǐng)問(wèn)，是不是有什么原因會(huì)很大程度的影響體模量的準(zhǔn)確性。還是我擬合體模量的方法有問(wèn)題，是否有其他的方法能得到較精確的模量值？

附上Si的原子優(yōu)化腳本：

#!/bin/sh
rm WAVECAR
for i in 5.15 5.20 5.25 5.30 5.35 5.40 5.43 5.455.50 5.53 5.55 5.60 5.65 5.70
do
cat >POSCAR< Si-fcc
$i
0.0 0.5 0.5
0.5 0.0 0.5
0.5 0.5 0.0
2
Direct
0.0 0.0 0.0
0.25 0.25 0.25
!
echo "a = $i angstrom"; time vasp

E=$(grep "TOTEN" OUTCAR |tail -1| awk '{printf "%12.6f \n", $5}')
V=$(grep "volume" OUTCAR |tail -1| awk '{printf "%12.4f \n", $5}')
echo $i $V $E >>EtVo.dat

done

INCAR：

SYSTEM = Si-fcc
ENCUT = 450
ISTART = 0; ICHARG = 2
ISMEAR = -5
IBRION = 2； NSW = 60
ISIF = 2
EDIFF = 1E-4； EDIFFG = -1E-3
PREC = Accurate

（2）關(guān)于體材料的自旋極化設(shè)置問(wèn)題：

   計(jì)算Co基磁性材料時(shí)，在INCAR中設(shè)置了參數(shù)ISPIN=2；可是靜態(tài)自洽scf計(jì)算、能帶結(jié)構(gòu)band計(jì)算等得到的OUTCAR中并未記錄magnetization信息；只有在態(tài)密度dos計(jì)算的OUTCAR文件里有s、p、d 以及Total的magnetization？而且在同種Co基磁性體系的不同結(jié)構(gòu)中，有的結(jié)構(gòu)似乎沒(méi)有自旋極化（ Total magnetization=0.0）；有的結(jié)構(gòu)卻有明顯的自旋極化現(xiàn)象（Total magnetization=2.0）；

   附上其中一個(gè)結(jié)構(gòu)的算DOS的INCAR：

SYSTEM = Co-work
ENCUT = 450
ISTART = 1; ICHARG = 11
ISMEAR = 0; SIGMA = 0.1
GGA = 91; VOSKOWN =1
ISPIN = 2
RWIGS = 2.163 1.512
PREC = Accurate

（3）關(guān)于LDA+U的設(shè)置：

   在計(jì)算Co基材料時(shí)，由于有3d電子的存在，所以想對(duì)GGA勢(shì)進(jìn)行“+U”計(jì)算。可是不知道其中的U值和J值該怎么設(shè)置？（LDAU、LDAUTYPE、LDAUL、LDAUU、LDAUJ）。。vasp手冊(cè)說(shuō)要參考實(shí)驗(yàn)文獻(xiàn)；具體要參考什么實(shí)驗(yàn)值？論壇上有人說(shuō)帶隙值，那如果是金屬，怎么辦？還有，如果沒(méi)有實(shí)驗(yàn)文獻(xiàn)值，那以什么來(lái)作為判據(jù)，選擇正確的U和J值。

附：上次一位高手的解答，但是還是很疑惑

引用回帖:

Originally posted by gavinliu7390 at 2010-12-23 20:31:58:
1，可以擬合物態(tài)方程來(lái)求體彈模量。
2，我也不是很懂。
3， U和J確實(shí)是用實(shí)驗(yàn)的值來(lái)擬合的。一般說(shuō)來(lái)是改變U和J來(lái)對(duì)比實(shí)驗(yàn)上的物態(tài)方程和平衡體積。

[ Last edited by fzuvivi on 2010-12-29 at 08:54 ]

回復(fù)此樓

» 收錄本帖的淘帖專輯推薦

科研相關(guān)資料

» 猜你喜歡

考研求調(diào)劑已經(jīng)有4人回復(fù)
085601專碩，總分342求調(diào)劑，地區(qū)不限已經(jīng)有4人回復(fù)
085601材料工程專碩求調(diào)劑已經(jīng)有5人回復(fù)
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有6人回復(fù)
能源材料化學(xué)課題組招收碩士研究生8-10名已經(jīng)有6人回復(fù)
296求調(diào)劑已經(jīng)有5人回復(fù)
0703化學(xué)調(diào)劑，六級(jí)已過(guò)，有科研經(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ù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

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

還是決定了13年考研，希望得到點(diǎn)大家的經(jīng)驗(yàn)，謝謝！已經(jīng)有229人回復(fù)
希望文章能順利接受已經(jīng)有29人回復(fù)
求助：origin高斯擬合已經(jīng)有10人回復(fù)
擬合過(guò)過(guò)彈性模量的沒(méi)？已經(jīng)有43人回復(fù)
提供一個(gè) Birch-Murnaghan 擬合的算法（更新1）已經(jīng)有95人回復(fù)
丙烯酸酯彈性體已經(jīng)有23人回復(fù)
金屬電沉積的極化曲線已經(jīng)有33人回復(fù)
為什么說(shuō)，單個(gè)高分子鏈已具備了高彈性？已經(jīng)有5人回復(fù)
對(duì)數(shù)據(jù)進(jìn)行曲線擬合后如何得到擬合誤差已經(jīng)有30人回復(fù)
鋼為什么會(huì)有彈性？已經(jīng)有5人回復(fù)
彈性常數(shù)計(jì)算出現(xiàn)詭異。。。。來(lái)看看嘍已經(jīng)有5人回復(fù)
如何做出擬合曲線求值已經(jīng)有24人回復(fù)
求助橢偏擬合已經(jīng)有3人回復(fù)
石墨烯中的自旋軌道耦合已經(jīng)有7人回復(fù)
超精細(xì)相互作用下的自旋翻轉(zhuǎn) 已經(jīng)有7人回復(fù)

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

查看全部散金貼

1樓 2010-12-27 22:21:28

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

gougou531

禁蟲 (正式寫手)

本帖內(nèi)容被屏蔽

4樓2011-01-05 09:14:18

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

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

fzuvivi

金蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 1168
帖子: 224
在線: 330.1小時(shí)
蟲號(hào): 881678

給自己頂一個(gè)吧。

回復(fù)此樓

2樓2010-12-29 08:56:47

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

銳利的碎片

木蟲 (正式寫手)

1ST強(qiáng)帖: 1
應(yīng)助: 136 (高中生)
金幣: 3637.1
帖子: 988
在線: 1224.9小時(shí)
蟲號(hào): 961933

★ ★
ellsaking(金幣+2):謝謝回帖交流~ 2010-12-29 09:55:09
fzuvivi(金幣+2): 2010-12-29 17:20:37

彈性模量可以用pbesol或者AM05，會(huì)有一定改善，也可以用雜化泛函。
U值可以用constrained lda算，一般都是取幾個(gè)來(lái)比較實(shí)驗(yàn)值。金屬可以看d帶的寬度，這個(gè)可以通過(guò)比較光電子譜分析。
磁性的話，有的結(jié)構(gòu)基態(tài)是反鐵磁的，所以要分別考慮各種情況。

贊一下(1人)

回復(fù)此樓

3樓2010-12-29 09:44:18

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

zzy870720z

榮譽(yù)版主 (文壇精英)

引用回帖:

Originally posted by gougou531 at 2011-01-05 09:14:18:

在CASTEP中怎樣設(shè)置LDA+U，在優(yōu)化時(shí)不能用LDA+U來(lái)優(yōu)化晶胞，那么實(shí)在計(jì)算中用嗎？

好像ms5.5可以支持+U優(yōu)化

贊一下

回復(fù)此樓

5樓2011-01-05 14:15:49

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

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

[考研] 材料專碩306英一數(shù)二 +7	z1z2z3879 2026-03-16	9/450	2026-03-17 17:31 by ccjequ
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3	1喜春8 2026-03-17	5/250	2026-03-17 17:18 by ber川cool子
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4	張vvvv 2026-03-15	6/300	2026-03-17 17:15 by ruiyingmiao
[考研] 303求調(diào)劑 +3	睿08 2026-03-17	3/150	2026-03-17 15:24 by 哦哦123
[考研] 302求調(diào)劑 +4	小賈同學(xué)123 2026-03-15	8/400	2026-03-17 10:33 by 小賈同學(xué)123
[考研] 藥學(xué)383 求調(diào)劑 +3	藥學(xué)chy 2026-03-15	4/200	2026-03-16 20:51 by 元子^0^
[考研] 復(fù)試調(diào)劑 +3	呼呼？~+123456 2026-03-14	3/150	2026-03-14 16:53 by WTUChen
[考研] 266求調(diào)劑 +4	學(xué)員97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo
[考研] 招收0805（材料）調(diào)劑 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、
[考研] 材料專碩288分求調(diào)劑一志愿211 +4	在家想你 2026-03-11	4/200	2026-03-13 22:49 by JourneyLucky
[考研] 0703，333分求調(diào)劑一志愿鄭州大學(xué)-物理化學(xué) +3	李魔女斗篷 2026-03-11	3/150	2026-03-13 22:24 by JourneyLucky
[考研] 求材料調(diào)劑 +5	隔壁陳先生 2026-03-12	5/250	2026-03-13 22:03 by 星空星月
[考研] 315求調(diào)劑 +9	小羊小羊_ 2026-03-11	10/500	2026-03-13 21:13 by SXNU李老師
[考研] 281求調(diào)劑 +9	Koxui 2026-03-12	11/550	2026-03-13 20:50 by Koxui
[考研] 材料與化工085600調(diào)劑求老師收留 +9	jiaanl 2026-03-11	9/450	2026-03-13 20:22 by JourneyLucky
[考研] 材料專碩350 求調(diào)劑 +4	王金科 2026-03-12	4/200	2026-03-13 16:02 by ruiyingmiao
[考研] 工科278分求調(diào)劑 +5	周慢熱啊 2026-03-12	7/350	2026-03-13 15:49 by JourneyLucky
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8	努力奮斗112 2026-03-12	9/450	2026-03-13 15:41 by JourneyLucky
[考研] 308求調(diào)劑 +3	是Lupa啊 2026-03-12	3/150	2026-03-13 14:30 by 求調(diào)劑zz
[考研] 大連大學(xué)化學(xué)專業(yè)研究生調(diào)劑 +3	琪久. 2026-03-10	8/400	2026-03-11 10:02 by 琪久.