版塊導航: 正在加載中...

登錄注冊

應《網(wǎng)絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

fzuvivi

金蟲 (小有名氣)

應助: 4 (幼兒園)
金幣: 1168
帖子: 224
在線: 330.1小時
蟲號: 881678

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

最近跑vasp，遇到了一些問題，

（1）關于Si 體彈性模量的計算：

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

附上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）關于體材料的自旋極化設置問題：

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

   附上其中一個結(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）關于LDA+U的設置：

   在計算Co基材料時，由于有3d電子的存在，所以想對GGA勢進行“+U”計算�？墒遣恢榔渲械腢值和J值該怎么設置？（LDAU、LDAUTYPE、LDAUL、LDAUU、LDAUJ）。。vasp手冊說要參考實驗文獻；具體要參考什么實驗值？論壇上有人說帶隙值，那如果是金屬，怎么辦？還有，如果沒有實驗文獻值，那以什么來作為判據(jù)，選擇正確的U和J值。

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

引用回帖:

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

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

回復此樓

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

科研相關資料

» 猜你喜歡

材料，紡織，生物（0856、0710），化學招生啦已經(jīng)有7人回復
301求調(diào)劑已經(jīng)有5人回復
281求調(diào)劑（0805）已經(jīng)有7人回復
材料與化工求調(diào)劑已經(jīng)有6人回復
299求調(diào)劑已經(jīng)有4人回復
277調(diào)劑已經(jīng)有6人回復
341求調(diào)劑已經(jīng)有7人回復
328求調(diào)劑，英語六級551，有科研經(jīng)歷已經(jīng)有8人回復
293求調(diào)劑已經(jīng)有12人回復
070300化學319求調(diào)劑已經(jīng)有4人回復

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

» 本主題相關價值貼推薦，對您同樣有幫助:

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

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

查看全部散金貼

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

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

fzuvivi

金蟲 (小有名氣)

應助: 4 (幼兒園)
金幣: 1168
帖子: 224
在線: 330.1小時
蟲號: 881678

給自己頂一個吧。

回復此樓

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

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

銳利的碎片

木蟲 (正式寫手)

1ST強帖: 1
應助: 136 (高中生)
金幣: 3637.1
帖子: 988
在線: 1224.9小時
蟲號: 961933

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

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

贊一下(1人)

回復此樓

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

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

gougou531

禁蟲 (正式寫手)

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

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

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

zzy870720z

榮譽版主 (文壇精英)

引用回帖:

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

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

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

贊一下

回復此樓

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

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

相關版塊跳轉(zhuǎn) 我要訂閱樓主 fzuvivi 的主題更新

返回列表

普通表情龍兔虎貓高級回復 (可上傳附件)

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

[考研] 341求調(diào)劑 +5	搗蛋豬豬 2026-03-11	7/350	2026-03-17 19:09 by 搗蛋豬豬
[考研] 26考研求調(diào)劑 +6	丶宏Sir 2026-03-13	6/300	2026-03-17 16:13 by 醉在風里
[考研] 303求調(diào)劑 +3	睿08 2026-03-17	3/150	2026-03-17 15:24 by 哦哦123
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +3	share_joy 2026-03-16	3/150	2026-03-17 13:41 by houyaoxu
[基金申請] 國自科面上基金字體 +6	iwuli 2026-03-12	7/350	2026-03-16 21:18 by sculhf
[文學芳草園] 伙伴們，祝我生日快樂吧 +17	myrtle 2026-03-10	26/1300	2026-03-16 18:32 by 青橙Ln
[考研] 環(huán)境工程調(diào)劑 +6	大可digkids 2026-03-16	6/300	2026-03-16 17:16 by barlinike
[考研] 0703化學調(diào)劑，六級已過，有科研經(jīng)歷 +7	曦熙兮 2026-03-15	7/350	2026-03-16 16:34 by houyaoxu
[基金申請] NSFC申報書里申請人簡歷中代表性論著還需要在申報書最后的附件里面再上傳一遍嗎 20+5	NSFC2026我來了 2026-03-10	14/700	2026-03-15 23:53 by 不負韶華的虎
[考研] 22408總分284求調(diào)劑 +3	InAspic 2026-03-13	3/150	2026-03-15 11:10 by zhq0425
[考研] 085601材料工程315分求調(diào)劑 +3	yang_0104 2026-03-15	3/150	2026-03-15 10:58 by peike
[考研] 材料080500調(diào)劑求收留 +3	一顆meteor 2026-03-13	3/150	2026-03-14 10:54 by peike
[考研] 285 求調(diào)劑資源與環(huán)境一志愿北京化工大學 +3	未名考生 2026-03-10	3/150	2026-03-13 23:04 by JourneyLucky
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8	dream…… 2026-03-12	8/400	2026-03-13 22:17 by 星空星月
[考研] 307求調(diào)劑 +5	超級伊昂大王 2026-03-12	5/250	2026-03-13 15:56 by 棒棒球手
[考研] 295求調(diào)劑 +3	小匕仔汁 2026-03-12	3/150	2026-03-13 15:17 by vgtyfty
[考研] 求調(diào)劑 +3	程雨杭 2026-03-12	3/150	2026-03-13 15:06 by JourneyLucky
[考研] 328化工專碩求調(diào)劑 +4	。，。，。，。i 2026-03-12	4/200	2026-03-13 14:44 by JourneyLucky
[考研] 0856化學工程280分求調(diào)劑 +4	shenzxsn 2026-03-11	4/200	2026-03-13 11:55 by ymwdoctor
[考研] 081200-11408-276學碩求調(diào)劑 +3	崔wj 2026-03-12	4/200	2026-03-12 19:33 by 求調(diào)劑zz

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

fzuvivi

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

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

» 猜你喜歡

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

» 本主題相關價值貼推薦，對您同樣有幫助:

» 搶金幣啦！回帖就可以得到: 查看全部散金貼

fzuvivi

銳利的碎片

gougou531

zzy870720z

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

» 本主題相關價值貼推薦，對您同樣有幫助:

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

查看全部散金貼