donggaomin

金蟲 (正式寫手)

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

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

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

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

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

請問各位前輩，我的計算思路或步驟有問題嗎？謝謝！

步驟一中所用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)太長就不列出來了

步驟三和四中的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.
!
若我還有哪些說得不清楚，請大家告訴我，我加上。

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

回復此樓

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

學習

» 猜你喜歡

289求調(diào)劑已經(jīng)有4人回復
292求調(diào)劑已經(jīng)有4人回復
080500求調(diào)劑已經(jīng)有3人回復
291求調(diào)劑已經(jīng)有3人回復
材料292調(diào)劑已經(jīng)有8人回復
一志愿211 初試270分求調(diào)劑已經(jīng)有6人回復
求調(diào)劑已經(jīng)有3人回復
274求調(diào)劑已經(jīng)有4人回復
300求調(diào)劑，材料科學英一數(shù)二已經(jīng)有5人回復
招08考數(shù)學已經(jīng)有15人回復

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

【求助】castep計算彈性常數(shù)問題已經(jīng)有14人回復
【求助】CASTEP計算彈性常數(shù) 已經(jīng)有15人回復
【求助】CASTEP計算彈性常數(shù)得到應力應變曲線已經(jīng)有17人回復
【求助】MS彈性常數(shù)計算分析問題已經(jīng)有5人回復
【求助】CASTEP服務器彈性常數(shù)的計算已經(jīng)有18人回復
【討論】用vasp計算過彈性常數(shù)的請進！！更新中��！已經(jīng)有52人回復

好好學習

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

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

xh512

金蟲 (正式寫手)

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

【答案】應助回帖

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

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

贊一下(1人)

回復此樓

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

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

donggaomin

金蟲 (正式寫手)

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

大家都來說說啊！

回復此樓

好好學習

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

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

donggaomin

金蟲 (正式寫手)

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

引用回帖:

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

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

贊一下

回復此樓

好好學習

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

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

xh512

金蟲 (正式寫手)

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

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

引用回帖:

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

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

贊一下(1人)

回復此樓

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

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

donggaomin

金蟲 (正式寫手)

應助: 1 (幼兒園)
金幣: 673.8
散金: 819
紅花: 9
帖子: 565
在線: 244.3小時
蟲號: 853446
注冊: 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方法加壓，文獻 PRB 79, 134102中對比了這兩種方法。PSTRESS也不一定不行，看你的體系了，只是目前EOS用的較多而已。

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

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

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

贊一下

回復此樓

好好學習

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

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

xh512

金蟲 (正式寫手)

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

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

引用回帖:

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

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

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

贊一下(1人)

回復此樓

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

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

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

[考研] 300求調(diào)劑，材料科學英一數(shù)二 +5	leaflight 2026-03-24	5/250	2026-03-24 16:25 by laoshidan
[考研] 招08考數(shù)學 +6	laoshidan 2026-03-20	15/750	2026-03-24 16:22 by laoshidan
[考研] 321求調(diào)劑 +4	Ymlll 2026-03-24	4/200	2026-03-24 14:44 by sprinining
[考研] 0856材料專碩353求調(diào)劑 +5	NIFFFfff 2026-03-20	5/250	2026-03-24 11:46 by 544594351
[考研] 一志愿山東大學藥學學碩求調(diào)劑 +3	開開心心沒煩惱 2026-03-23	4/200	2026-03-24 00:06 by 開開心心沒煩惱
[論文投稿] 急發(fā)核心期刊論文 +3	賢達問津 2026-03-23	5/250	2026-03-23 17:13 by 妹子不好惹
[考研] 接收2026碩士調(diào)劑(學碩+專碩) +4	allen-yin 2026-03-23	6/300	2026-03-23 15:04 by 汪��？！
[考研] 求調(diào)劑材料學碩080500，總分289分 5+3	@taotao 2026-03-19	21/1050	2026-03-23 10:17 by 冠c哥
[考研] 一志愿070300浙大化學358分，求調(diào)劑！ +4	酥酥魚.. 2026-03-21	4/200	2026-03-23 08:12 by Iveryant
[考研] 323求調(diào)劑 +6	洼小桶 2026-03-18	6/300	2026-03-23 00:29 by king123！
[考研] 260求調(diào)劑 +3	朱芷琳 2026-03-20	4/200	2026-03-22 15:12 by 朱芷琳
[考研] 269專碩求調(diào)劑 +6	金恩貝 2026-03-21	6/300	2026-03-22 14:31 by ColorlessPI
[考研] 0805 316求調(diào)劑 +3	大雪深藏 2026-03-18	3/150	2026-03-21 18:55 by 學員8dgXkO
[考研] 0703化學297求調(diào)劑 +3	Daisy☆ 2026-03-20	3/150	2026-03-21 17:45 by ColorlessPI
[考研] 302求調(diào)劑 +12	呼呼呼。。。。 2026-03-17	12/600	2026-03-21 17:29 by ColorlessPI
[考研] 0805材料320求調(diào)劑 +3	深海物語 2026-03-20	3/150	2026-03-21 15:46 by 無際的草原
[考研] 材料與化工（0856）304求 B區(qū) 調(diào)劑 +3	邱gl 2026-03-21	3/150	2026-03-21 13:47 by lature00
[考研] 一志愿西北大學，070300化學學碩，總分287，雙非一本，求調(diào)劑。 +4	晨昏線與星海 2026-03-19	4/200	2026-03-20 22:15 by JourneyLucky
[考研] 材料與化工 322求調(diào)劑 +4	然11 2026-03-19	4/200	2026-03-20 22:12 by luoyongfeng
[考研] 一志愿西安交通大學學碩 354求調(diào)劑211或者雙一流 +3	我想要讀研究生 2026-03-20	3/150	2026-03-20 20:13 by JourneyLucky