 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 1343 | 回復(fù): 2 | |||
MJFaber銅蟲 (初入文壇)
|
[求助]
彈性常數(shù)計(jì)算
|
|
最近在學(xué)習(xí)VASP計(jì)算彈性常數(shù)，所以試算了一下Pd的彈性常數(shù)與文獻(xiàn)里面計(jì)算的做對比。 但是結(jié)果差距好大，不知道是什么原因。請大神指導(dǎo)啊  下面是結(jié)果和我的計(jì)算腳本。 C11 C12 C44 文獻(xiàn) 274.3 146.3 71.4 自己 204.7 150.8 57.8 #!/bin/bash -l #SBATCH -t 20:00:00 #SBATCH -p node #SBATCH -n 2 #SBATCH -J test-1 { echo "Cluster: " $NSC_RESOURCE_NAME echo "Running on: "$SLURM_NODELIST echo "Job ID: "$SLURM_JOBID } > vasp.log module load vasp module load intel/17.1 module load openmpi/2.0.2 cat /home/POTCARs/PBE/2015-09-21/Pd/POTCAR > POTCAR cwd=`pwd` #####結(jié)構(gòu)優(yōu)化 cat >INCAR <<! PREC = Accurate ENCUT = 350 ISTART = 0 ICHARG = 2 EDIFF = 1E-8 EDIFFG = -0.001 GGA=PE #NPAR=4 #### parallel part LREAL=.FALSE. LWAVE=.FALSE. LCHARG=.FALSE. ### Relaxation ISMEAR = 1 SIGMA = 0.1 POTIM = 0.1 NSW = 100 IBRION = 2 ISIF = 3 ! cat >POSCAR <<! Pd-fcc 1.0 3.8900001049 0.0000000000 0.0000000000 0.0000000000 3.8900001049 0.0000000000 0.0000000000 0.0000000000 3.8900001049 Pd 4 Direct 0.000000000 0.000000000 0.000000000 0.500000000 0.500000000 0.000000000 0.000000000 0.500000000 0.500000000 0.500000000 0.000000000 0.500000000 ! cat >KPOINTS <<! Automatic 0 Monkhorst Pack 20 20 20 0 0 0 ! mpirun /sw/apps/vasp/5.4.1_rackham_openmpi_intelmkl/vasp >> vasp1.log 2>&1 cp CONTCAR POSCAR #####靜態(tài)計(jì)算 cat >INCAR <<! PREC = Accurate ENCUT = 350 ISTART = 0 ICHARG = 2 EDIFF = 1E-8 EDIFFG = -0.001 GGA=PE #NPAR=4 #### parallel part LREAL=.FALSE. LWAVE=.TURE. LCHARG=.TURE. ### Relaxation ISMEAR = 1 NSW = 0 IBRION = -1 ISIF = 2 ! mpirun /sw/apps/vasp/5.4.1_rackham_openmpi_intelmkl/vasp >> vasp2.log 2>&1 cp CONTCAR POSCAR ####計(jì)算彈性常數(shù) cat >INCAR <<! System =Pd-fcc ISTART = 0 ICHARG = 2 ENCUT = 500 PREC = Accurate GGA=PE ISMEAR = 1 SIGMA = 0.1 POTIM = 0.015 NSW = 1 IBRION = 6 ISIF = 3 NFREE=4 EDIFF = 1E-8 EDIFFG = -0.001 SYMPREC=1E-10 LREAL=.FALSE. LWAVE=.FALSE. LCHARG=.FALSE. ! cat >KPOINTS <<! Automatic 0 Monkhorst Pack 40 40 40 0 0 0 ! mpirun /sw/apps/vasp/5.4.1_rackham_openmpi_intelmkl/vasp >> vasp3.log 2>&1 cp OUTCAR OUTCAR_3 |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 第一性原理計(jì)算 |
» 猜你喜歡
302求調(diào)劑一志愿華中師范大學(xué)
已經(jīng)有4人回復(fù)
-
288求調(diào)劑 一志愿哈工大 材料與化工
已經(jīng)有32人回復(fù)
-
一志愿鄭大材料工程290求調(diào)劑
已經(jīng)有20人回復(fù)
-
環(huán)境工程297分求調(diào)劑一志愿杭高院
已經(jīng)有12人回復(fù)
-
366求調(diào)劑一志愿東北大學(xué)
已經(jīng)有7人回復(fù)
-
材料調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿北京科技大學(xué)085601材料工程英一數(shù)二初試總分335求調(diào)劑
已經(jīng)有9人回復(fù)
-
材料調(diào)劑
已經(jīng)有5人回復(fù)
-
材料與化工306分找調(diào)劑
已經(jīng)有5人回復(fù)
-
301求調(diào)劑
已經(jīng)有5人回復(fù)
孫瑩潔 
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 17.5
- 帖子: 12
- 在線: 29.3小時(shí)
- 蟲號: 2312007
- 注冊: 2013-03-02
- 專業(yè): 金屬材料的力學(xué)行為
|
如果是新版本的vasp可以直接一部計(jì)算的，還有你查看一下你的vasp.log文件看一下里面的warning,應(yīng)該會(huì)有幫助 發(fā)自小木蟲Android客戶端 |
2樓2017-11-20 23:04:38
|
本帖內(nèi)容被屏蔽 |
3樓2018-10-17 12:21:55
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 MJFaber 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 302求調(diào)劑一志愿華中師范大學(xué) +4 | 小江小江江江 2026-04-02 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 085600 295分求調(diào)劑 +16 | W55j 2026-03-30 | 19/950 |
|
|
|
[考研] 一志愿廈門大學(xué)材料工程專碩354找調(diào)劑�。。� +7 | 貝唄鋇鋇 2026-03-30 | 7/350 |
|
|
|
[考研] 找調(diào)劑 +6 | 楚喬喬 2026-04-01 | 6/300 |
|
|
|
[考研] 考研調(diào)劑0702 +3 | 梅子味晚霞 2026-03-28 | 3/150 |
|
|
|
[考研] 考研調(diào)劑 +12 | Amber00 2026-03-31 | 12/600 |
|
|
|
[考研] 285求調(diào)劑 +11 | AZMK 2026-04-01 | 11/550 |
|
|
|
[考研] 292求調(diào)劑 +17 | 木蟲er12138 2026-04-01 | 17/850 |
|
|
|
[考研] 070300化學(xué)279求調(diào)劑 +15 | 哈哈哈^_^ 2026-03-31 | 17/850 |
|
|
|
[考研] 求調(diào)劑 +4 | DADA怪 2026-03-31 | 4/200 |
|
|
|
[考研] 070300求調(diào)劑306分 +5 | 26要上岸 2026-03-27 | 5/250 |
|
|
|
[考研] 333求調(diào)劑 +4 | 阿科逸 2026-03-31 | 4/200 |
|
|
|
[考研] 318求調(diào)劑 +3 | 篤行致遠(yuǎn). 2026-03-31 | 3/150 |
|
|
|
[考研] 0856求調(diào)劑 +9 | 楒桉 2026-03-28 | 9/450 |
|
|
|
[考研] 求化學(xué)調(diào)劑 +12 | wulanna 2026-03-28 | 12/600 |
|
|
|
[考研] 材料與化工272求調(diào)劑 +25 | 阿斯蒂芬2004 2026-03-28 | 25/1250 |
|
|
|
[考研] 286求調(diào)劑 +5 | 丟掉懶惰 2026-03-27 | 8/400 |
|
|
|
[考研] 262求調(diào)劑 +7 | ZZ..000 2026-03-30 | 8/400 |
|
|
|
[考研] 316求調(diào)劑 +7 | 江辭666 2026-03-26 | 7/350 |
|
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|