版塊導(dǎo)航: 正在加載中...

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

24小時(shí)熱門版塊排行榜

返回列表

donggaomin

金蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 673.8
散金: 819
紅花: 9
帖子: 565
在線: 244.3小時(shí)
蟲號(hào): 853446
注冊(cè): 2009-09-21
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[求助] 加壓后的體系怎么計(jì)算彈性常數(shù)？（或其它性質(zhì)）

請(qǐng)問：侯博士教程上的計(jì)算彈性常數(shù)的應(yīng)力-應(yīng)變法和能量—應(yīng)變法，是否適用于加了靜水壓的情況？是否適用于加了單軸壓力的情況？
加了壓力之后的計(jì)算與沒有加壓的計(jì)算不同在哪里？
謝謝！

我現(xiàn)在的思路：
我用PSTRESS是給體系加壓（PSTRESS=1500;ISIF=3）得到加了壓強(qiáng)之后的晶格常數(shù)及原子坐標(biāo)；然后就把它當(dāng)成poscar, 計(jì)算后面的應(yīng)力-應(yīng)變關(guān)系。

步驟：
一、用PSTRESS=1500;ISIF=3...得加壓下的CONTCAR;（后面的計(jì)算中都不再設(shè)置PSTRESS）
二、cp CONTCAR OLDPOS; 用defvector.x 施加應(yīng)變得到應(yīng)變后的POSCAR；
三、固定晶格常數(shù)馳豫原子位置（ISIF=2，IBRION=2）；
四、靜態(tài)計(jì)算得精確的應(yīng)力或能量，用侯博士所說的應(yīng)力-應(yīng)變法擬合出cij.(注：侯博士介紹的應(yīng)力-應(yīng)變法沒有加壓)

請(qǐng)問各位前輩，我的計(jì)算思路或步驟有問題嗎？謝謝！

步驟一中所用INCAR：
SYSTE=AL
ENCUT = 800
ISTART = 0
ICHARG = 2
ISMEAR = 0; SIGMA = 0.2
NSW = 60; IBRION = 2
EDIFF = 1E-5
EDIFFG = -1E-2
ISIF = 3
POTIM = 0.2
PSTRESS=1500
PREC = Accurate
LWAVE = .FALSE.
LCHARG =.FALSE.

步驟二中的defvector.f(.x)太長(zhǎng)就不列出來了

步驟三和四中的INCAR：
for i in -0.003 -0.0006 -0.0003 0.00 \
0.0003 0.0006 0.003 0.005
do
echo $i | ./defvector.x
cp fort.3 POSCAR
####
cat > INCAR < SYSTEM = Al-fcc
ENCUT = 800
ISTART = 0
ICHARG = 2
ISMEAR = 0; SIGMA = 0.2
NSW = 60; IBRION = 2
EDIFF = 1E-5
EDIFFG = -1E-2
ISIF = 2
POTIM = 0.2
PREC = Accurate
LWAVE = .FALSE.
LCHARG =.FALSE.
!
echo "delta = $i ";
$BINDIR/mpirun -np $NSLOTS -machinefile $TMPDIR/machines -nolocal \
$BINDIR/vasp5.2mpi_pgi
cp CONTCAR pos.$i
cp CONTCAR POSCAR
cat > INCAR < SYSTEM = Al-fcc
ENCUT = 800
ISTART = 0
ICHARG = 2
ISMEAR = -5
EDIFF = 1E-5
PREC = Accurate
LWAVE = .FALSE.
LCHARG =.FALSE.
!
若我還有哪些說得不清楚，請(qǐng)大家告訴我，我加上。

[ Last edited by donggaomin on 2011-4-20 at 09:53 ]

回復(fù)此樓

好好學(xué)習(xí)

1樓 2011-04-19 23:05:22

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

xh512

金蟲 (正式寫手)

應(yīng)助: 2 (幼兒園)
金幣: 149.1
散金: 1339
紅花: 11
帖子: 309
在線: 923.6小時(shí)
蟲號(hào): 461655
注冊(cè): 2007-11-18
專業(yè): 金屬材料的合金相、相變及

【答案】應(yīng)助回帖

★
donggaomin(金幣+10): 謝謝您！您說的是不是步驟三？步驟四（靜態(tài)計(jì)算）呢？ 2011-04-20 10:00:37
donggaomin(金幣+20): 謝謝您提供的參考文獻(xiàn)PRB 79, 134102，上面應(yīng)該有詳細(xì)的EOS方法介紹，我去學(xué)學(xué) 2011-04-20 11:02:03
zsl_321123(金幣+1): 感謝建議 2011-04-29 14:17:05

施加應(yīng)變的時(shí)候也得加上PSTRESS才行。另外，提醒一下，PSTRESS方式加壓可能
不大靠譜。侯博士提過，PSTRESS搜索空間有限。

贊一下(1人)

回復(fù)此樓

2樓2011-04-20 09:54:03

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

donggaomin

金蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 673.8
散金: 819
紅花: 9
帖子: 565
在線: 244.3小時(shí)
蟲號(hào): 853446
注冊(cè): 2009-09-21
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

大家都來說說��！

回復(fù)此樓

好好學(xué)習(xí)

3樓2011-04-20 09:55:38

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

donggaomin

金蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 673.8
散金: 819
紅花: 9
帖子: 565
在線: 244.3小時(shí)
蟲號(hào): 853446
注冊(cè): 2009-09-21
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

Originally posted by xh512 at 2011-04-20 09:54:03:
施加應(yīng)變的時(shí)候也得加上PSTRESS才行。另外，提醒一下，PSTRESS方式加壓可能
不大靠譜。侯博士提過，PSTRESS搜索空間有限。

謝謝！
請(qǐng)問：您說的是不是步驟三也要加上PSTRESS?
另外，不用PSTRESS加壓的話，您指的應(yīng)該是通過改變晶格常數(shù)來加壓，但據(jù)說改變晶格常數(shù)只能給立方結(jié)構(gòu)加壓，請(qǐng)您再指點(diǎn)一二！

贊一下

回復(fù)此樓

好好學(xué)習(xí)

4樓2011-04-20 10:05:13

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

xh512

金蟲 (正式寫手)

應(yīng)助: 2 (幼兒園)
金幣: 149.1
散金: 1339
紅花: 11
帖子: 309
在線: 923.6小時(shí)
蟲號(hào): 461655
注冊(cè): 2007-11-18
專業(yè): 金屬材料的合金相、相變及

★
zsl_321123(金幣+1): 鼓勵(lì)交流 2011-04-29 14:17:26

引用回帖:

Originally posted by donggaomin at 2011-04-20 10:05:13:
謝謝！
請(qǐng)問：您說的是不是步驟三也要加上PSTRESS?
另外，不用PSTRESS加壓的話，您指的應(yīng)該是通過改變晶格常數(shù)來加壓，但據(jù)說改變晶格常數(shù)只能給立方結(jié)構(gòu)加壓，請(qǐng)您再指點(diǎn)一二！

非立方結(jié)構(gòu)也能通過EOS方法加壓，文獻(xiàn) PRB 79, 134102中對(duì)比了這兩種方法。PSTRESS也不一定不行，看你的體系了，只是目前EOS用的較多而已。

贊一下(1人)

回復(fù)此樓

5樓2011-04-20 10:42:57

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

donggaomin

金蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 673.8
散金: 819
紅花: 9
帖子: 565
在線: 244.3小時(shí)
蟲號(hào): 853446
注冊(cè): 2009-09-21
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

Originally posted by xh512 at 2011-04-20 10:42:57:
非立方結(jié)構(gòu)也能通過EOS方法加壓，文獻(xiàn) PRB 79, 134102中對(duì)比了這兩種方法。PSTRESS也不一定不行，看你的體系了，只是目前EOS用的較多而已。

您好！請(qǐng)教一下EOS方法加壓：
我下載了兩個(gè)EOS擬合程序，有可視界面操作（見附件），分別是P-V和P-V-T，界面如下圖所示：

我的理解是：選好擬合函數(shù)（如Vinet），填入V0，輸入所加壓強(qiáng)P，則得壓強(qiáng)P下的體積V，從而得出此時(shí)的晶格常數(shù)。

疑問一：我上面的理解是否就是EOS方法加壓？如果不是，那應(yīng)該是怎樣的？
疑問二：輸入?yún)?shù)中有一個(gè)是體彈性模量K0，它也會(huì)影響所得的V，請(qǐng)問K0如何定？
疑問三：這些加壓擬合是否適用于所有的元素及元素的組合？（如PbTiO3等）

贊一下

回復(fù)此樓

好好學(xué)習(xí)

6樓2011-04-26 22:43:31

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

xh512

金蟲 (正式寫手)

應(yīng)助: 2 (幼兒園)
金幣: 149.1
散金: 1339
紅花: 11
帖子: 309
在線: 923.6小時(shí)
蟲號(hào): 461655
注冊(cè): 2007-11-18
專業(yè): 金屬材料的合金相、相變及

★
zsl_321123(金幣+1): 感謝建議 2011-04-29 14:18:13

引用回帖:

Originally posted by donggaomin at 2011-04-26 22:43:31:
您好！請(qǐng)教一下EOS方法加壓：
我下載了兩個(gè)EOS擬合程序，有可視界面操作（見附件），分別是P-V和P-V-T，界面如下圖所示：

我的理解是：選好擬合函數(shù)（如Vi ...

不好意思，才注意到回帖。
沒了解過這個(gè)程序。
我感覺K0可能是平衡條件下的體彈性模量，相應(yīng)的K0'就是偏導(dǎo)體模量了。
最好是看說明書吧。或者輸入不同的值驗(yàn)證下。
我覺得這個(gè)程序應(yīng)該對(duì)任何體系都適用。
既然是P-V關(guān)系，輸入壓強(qiáng)，當(dāng)然得到對(duì)應(yīng)壓強(qiáng)下的體積了。
不確定的話，就想程序作者寫信吧，
祝你好運(yùn)。

贊一下(1人)

回復(fù)此樓

7樓2011-04-27 18:57:42

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

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

返回列表

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

[考研] 279 分求調(diào)劑 +3	睡個(gè)好覺_16 2026-03-24	3/150	2026-03-26 15:43 by zzll406
[考研] 263求調(diào)劑 +6	yqdszhdap－ 2026-03-22	10/500	2026-03-26 13:11 by 公瑾逍遙
[考研] 292求調(diào)劑 +8	鵝鵝鵝額額額額?/a> 2026-03-25	9/450	2026-03-26 12:21 by 張凱十八號(hào)
[考研] 一志愿北京化工大學(xué)材料與化工 264分各科過A區(qū)國(guó)家線 +6	哈哈157349 2026-03-21	6/300	2026-03-26 10:39 by 醉在風(fēng)里
[考研] 347求調(diào)劑 +4	L when 2026-03-25	4/200	2026-03-25 13:37 by cocolv
[考研] 0854人工智能方向招收調(diào)劑 +4	章小魚567 2026-03-24	4/200	2026-03-25 13:29 by 2177681040
[考研] 求調(diào)劑 +3	李李不服輸 2026-03-25	3/150	2026-03-25 13:03 by cmz0325
[考研] 292求調(diào)劑 +4	鵝鵝鵝額額額額?/a> 2026-03-24	4/200	2026-03-24 16:41 by peike
[考研] 279分求調(diào)劑一志愿211 +18	chaojifeixia 2026-03-19	20/1000	2026-03-24 10:34 by dolphin_ycj
[考研] 一志愿國(guó)科過程所081700，274求調(diào)劑 +3	三水研0水立方 2026-03-23	3/150	2026-03-23 23:11 by MajorWen
[考研] 工科0856求調(diào)劑 +5	沐析汀汀 2026-03-21	5/250	2026-03-23 17:56 by 海瑟薇-
[考研] 一志愿070300浙大化學(xué)358分，求調(diào)劑！ +4	酥酥魚.. 2026-03-21	4/200	2026-03-23 08:12 by Iveryant
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +5	困困困困坤坤 2026-03-20	6/300	2026-03-22 17:41 by hxsm
[考研] 324求調(diào)劑 +6	lucky呀呀呀鴨 2026-03-20	6/300	2026-03-22 16:01 by ColorlessPI
[考研] 一志愿華中科技大學(xué)071000，求調(diào)劑 +4	沿岸有貝殼6 2026-03-21	4/200	2026-03-22 07:21 by ilovexiaobin
[考研] 0703化學(xué)調(diào)劑 +4	妮妮ninicgb 2026-03-21	4/200	2026-03-21 18:39 by 學(xué)員8dgXkO
[考研] 求助 +5	夢(mèng)里的無言 2026-03-21	6/300	2026-03-21 17:51 by 學(xué)員8dgXkO
[考研] 336求調(diào)劑 +5	rmc8866 2026-03-21	5/250	2026-03-21 17:24 by 學(xué)員8dgXkO
[考研] 一志愿西南交通專碩材料355 本科雙非求調(diào)劑 +5	西南交通專材355 2026-03-19	5/250	2026-03-20 21:10 by JourneyLucky
[考研] 261求B區(qū)調(diào)劑，科研經(jīng)歷豐富 +3	牛奶很忙 2026-03-20	4/200	2026-03-20 19:34 by JourneyLucky

24小時(shí)熱門版塊排行榜

[求助] 加壓后的體系怎么計(jì)算彈性常數(shù)？（或其它性質(zhì)）

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

» 猜你喜歡

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

【答案】應(yīng)助回帖