| 5 | 1/1 | 返回列表 |
| 查看: 2253 | 回復(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)太長就不列出來了 步驟三和四中的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ā)表 | |
|---|---|---|---|---|
|
[考研] 329求調(diào)劑 +8 | 星野? 2026-03-26 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 322求調(diào)劑 +8 | 宋明欣 2026-03-27 | 8/400 |
|
|
[考研] 277跪求調(diào)劑 +7 | 1915668 2026-03-27 | 11/550 |
|
|
[考研] 085601材料工程找調(diào)劑 +16 | oatmealR 2026-03-29 | 17/850 |
|
|
[考研] 337求調(diào)劑 +6 | 《樹》 2026-03-29 | 6/300 |
|
|
[考研] 求調(diào)劑 +10 | 張zz111 2026-03-27 | 11/550 |
|
|
[考研] 318一志愿吉林大學(xué)生物與醫(yī)藥 求調(diào)劑 +5 | 篤行致遠(yuǎn). 2026-03-28 | 5/250 |
|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +14 | Mr. Z 2026-03-25 | 14/700 |
|
|
[考研] 數(shù)一英一271專碩(085401)求調(diào)劑,可跨 +7 | 前行必有光 2026-03-28 | 8/400 |
|
|
[考研] 11408軟件工程求調(diào)劑 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
[考研] 347求調(diào)劑 +3 | 山頂見α 2026-03-25 | 3/150 |
|
|
[考研] 材料277求調(diào)劑 +7 | min3 2026-03-24 | 7/350 |
|
|
[考研] 一志愿上海理工能源動(dòng)力(085800)310分求調(diào)劑 +3 | zhangmingc 2026-03-27 | 4/200 |
|
|
[考研] 085600,材料與化工321分,求調(diào)劑 +9 | 大饞小子 2026-03-27 | 9/450 |
|
|
[考研] 調(diào)劑推薦 +5 | 清酒714 2026-03-26 | 6/300 |
|
|
[考研] 一志愿哈工大,085400,320,求調(diào)劑 +4 | gdlf9999 2026-03-24 | 4/200 |
|
|
[考研] 303求調(diào)劑 +6 | 藍(lán)山月 2026-03-25 | 6/300 |
|
|
[考研] 材料專碩 335 分求調(diào)劑 +4 | 拒絕冷暴力 2026-03-25 | 4/200 |
|
|
[考研] 考研一志愿蘇州大學(xué)初始315(英一)求調(diào)劑 +3 | sbdksD 2026-03-24 | 4/200 |
|
|
[考研]
|
13659058978 2026-03-24 | 4/200 |
|