| 5 | 1/1 | 返回列表 |
| 查看: 2254 | 回復(fù): 6 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
donggaomin金蟲 (正式寫手)
|
[求助]
加壓后的體系怎么計(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 ] |
學(xué)習(xí) |

金蟲 (正式寫手)

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

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿 南京航空航天大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +4 | @taotao 2026-03-30 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 085600 295分求調(diào)劑 +7 | W55j 2026-03-30 | 9/450 |
|
|
[考研] 286分調(diào)劑 +10 | Faune 2026-03-30 | 11/550 |
|
|
[考研] 一志愿:西北大學(xué),英一數(shù)一408-284分求調(diào)劑 +5 | 12.27 2026-03-27 | 5/250 |
|
|
[考研] 一志愿鄭州大學(xué),080500學(xué)碩,總分317分求調(diào)劑 +10 | 舉個(gè)栗子oi 2026-03-24 | 11/550 |
|
|
[考研] 0703化學(xué) +14 | 妮妮ninicgb 2026-03-27 | 14/700 |
|
|
[考研] 291求調(diào)劑 +5 | Y-cap 2026-03-29 | 6/300 |
|
|
[考研] 305求調(diào)劑 +8 | RuiFairyrui 2026-03-28 | 8/400 |
|
|
[考研] 本科雙非材料,跨考一志愿華電085801電氣,283求調(diào)劑,任何專業(yè)都可以 +6 | 芝士雪baoo 2026-03-28 | 8/400 |
|
|
[考研] 321求調(diào)劑 +7 | 璞玉~~ 2026-03-25 | 8/400 |
|
|
[考研] 312,生物學(xué)求調(diào)劑 +3 | 小譯同學(xué)abc 2026-03-28 | 3/150 |
|
|
[考研] 086000調(diào)劑 +3 | 7901117076 2026-03-26 | 3/150 |
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
[考研] 341求調(diào)劑 +7 | 青檸檬1 2026-03-26 | 7/350 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
[考研] 機(jī)械學(xué)碩310分,數(shù)一英一,一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|
|
[考研]
|
WWW西西弗斯 2026-03-24 | 8/400 |
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|
|
[考研]
|
13659058978 2026-03-24 | 4/200 |
|