| 查看: 982 | 回復(fù): 10 | |||
[求助]
壓力下的電子結(jié)構(gòu)及彈性計(jì)算
|
| 小弟有個(gè)問題想請(qǐng)教各位蟲友,請(qǐng)問計(jì)算壓力下的電子結(jié)構(gòu)和彈性時(shí),在晶體結(jié)構(gòu)優(yōu)化完后,進(jìn)行態(tài)密度、能帶和彈性性質(zhì)計(jì)算時(shí),要不要加壓力?希望獲得各位蟲友的回答,謝謝! |
木蟲 (正式寫手)
木蟲 (著名寫手)
木蟲 (正式寫手)
榮譽(yù)版主 (知名作家)
老和山猥瑣派九段
|
可是我仔細(xì)對(duì)比我的彈性計(jì)算結(jié)果,在沒有應(yīng)變的情況下計(jì)算出的體系自由能與在計(jì)算彈性之前的結(jié)構(gòu)優(yōu)化后的體系自由能的數(shù)據(jù)是相同的,這又是怎么回事?下面是我的計(jì)算文件設(shè)置,麻煩大家?guī)臀铱匆幌拢?br /> (1)在壓力15GPa下晶格優(yōu)化,得到: free energy TOTEN = -27.116931 eV energy without entropy= -27.116931 energy(sigma->0) = -27.116931 (2)彈性常數(shù)C11、C12計(jì)算 script.sh文件參數(shù): ENCUT=520 ISTART=0 ICHARG=2 ISMEAR=0;SIGMA=0.2 NSW=80;IBRION=2 POTIM=0.2 IALGO=48 EDIFF=5E-06 ISIF=2 PREC=Normal LWAVE=.FALSE. LCHARG=.FALSE. ! echo "delta=$i"; mpirun -np 2 vasp5.2 cp CONTCAR pos.$i cp CONTCAR POSCAR cat > INCAR < SYSTEM= ENCUT=520 ISTART=0 ICHARG=2 ISMEAR=0 IALGO=48 EDIFF=5E-6 PREC=Normal LWAVE=.FALSE. LCHARG=.FALSE. 計(jì)算結(jié)果: -0.018 1 F= -.26572346E+02 E0= -.26572346E+02 -0.015 1 F= -.26674069E+02 E0= -.26674069E+02 -0.012 1 F= -.26770758E+02 E0= -.26770758E+02 -0.009 1 F= -.26863577E+02 E0= -.26863577E+02 -0.006 1 F= -.26952164E+02 E0= -.26952164E+02 -0.003 1 F= -.27036589E+02 E0= -.27036589E+02 0.00 1 F= -.27116931E+02 E0= -.27116931E+02 0.003 1 F= -.27193287E+02 E0= -.27193287E+02 0.006 1 F= -.27265715E+02 E0= -.27265715E+02 當(dāng)應(yīng)變超過0時(shí),體系總能仍在下降,這是為什么?明明在應(yīng)變?yōu)?的時(shí)候,體系的總能和優(yōu)化的總能相一致,可是隨著應(yīng)變的變化,體系的總能卻仍然降低。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考! +4 | lemonzzn 2026-03-17 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 一志愿211,0703化學(xué)310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-20 | 7/350 |
|
|
[考研] 求調(diào)劑 +4 | 要好好無聊 2026-03-21 | 4/200 |
|
|
[考研] 工科0856求調(diào)劑 +3 | 沐析汀汀 2026-03-21 | 3/150 |
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
[考研] 22 350 本科985求調(diào)劑,求老登收留 +3 | 李軼男003 2026-03-20 | 3/150 |
|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +5 | 守候夕陽CF 2026-03-18 | 5/250 |
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +6 | 想上岸的鯉魚 2026-03-18 | 7/350 |
|
|
[考研] 一志愿西南交大,求調(diào)劑 +5 | 材化逐夢(mèng)人 2026-03-18 | 5/250 |
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 3/150 |
|
|
[考研] 319求調(diào)劑 +3 | 小力氣珂珂 2026-03-20 | 3/150 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 一志愿福大288有機(jī)化學(xué),求調(diào)劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|