版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門(mén)版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

MJFaber

銅蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 291
散金: 100
帖子: 46
在線: 22.9小時(shí)
蟲(chóng)號(hào): 3769269
注冊(cè): 2015-03-27
專業(yè): 力學(xué)中的基本問(wèn)題和方法

[求助] 彈性常數(shù)計(jì)算

最近在學(xué)習(xí)VASP計(jì)算彈性常數(shù)，所以試算了一下Pd的彈性常數(shù)與文獻(xiàn)里面計(jì)算的做對(duì)比。
但是結(jié)果差距好大，不知道是什么原因。請(qǐng)大神指導(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ù)

1樓 2017-10-25 15:36:56

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

bcwah717

禁蟲(chóng) (小有名氣)

本帖內(nèi)容被屏蔽

3樓2018-10-17 12:21:55

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

查看全部 3 個(gè)回答

孫瑩潔

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 17.5
帖子: 12
在線: 29.3小時(shí)
蟲(chóng)號(hào): 2312007
注冊(cè): 2013-03-02
專業(yè): 金屬材料的力學(xué)行為

如果是新版本的vasp可以直接一部計(jì)算的，還有你查看一下你的vasp.log文件看一下里面的warning,應(yīng)該會(huì)有幫助

發(fā)自小木蟲(chóng)Android客戶端

贊一下(2人)

回復(fù)此樓

2樓2017-11-20 23:04:38

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

查看全部 3 個(gè)回答

普通表情龍兔虎貓

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

[考研] 302求調(diào)劑一志愿華中師范大學(xué) +4	小江小江江江 2026-04-02	4/200	2026-04-02 14:49 by 5896
[考研] 初試301，代碼085701環(huán)境工程，本碩一致，四六級(jí)已過(guò)，有二區(qū)一作，共發(fā)表5篇論文 +6	axibli 2026-04-01	6/300	2026-04-02 13:42 by Ecowxq666！
[考研] 311求調(diào)劑 +14	藍(lán)月亮亮 2026-03-30	14/700	2026-04-02 12:18 by 1753564080
[考研] 一志愿北交大材料工程總分358 +8	cs0106 2026-04-01	9/450	2026-04-02 10:36 by 不吃魚(yú)的貓
[考研] 化學(xué)工程專碩324分，一志愿中國(guó)礦業(yè)大學(xué)求調(diào)劑 +7	耿耿1314 2026-04-01	7/350	2026-04-02 07:40 by 尚水閣主
[考研] 330分求調(diào)劑 +11	qzenlc 2026-03-29	11/550	2026-04-01 14:32 by chenqifeng666
[考研] 考研材料工程351分調(diào)劑 +5	整個(gè)好的 2026-03-31	5/250	2026-04-01 09:36 by topgun2009
[考研] 070300化學(xué)專業(yè)279調(diào)劑 +10	哈哈哈^_^ 2026-03-31	10/500	2026-03-31 23:13 by liu823948201
[考研] 080500-315分復(fù)試調(diào)劑 +9	上岸3821 2026-03-31	9/450	2026-03-31 17:29 by 唐沐兒
[考研] 340求調(diào)劑 +4	希望如此i 2026-03-31	4/200	2026-03-31 16:40 by 690616278
[考研] 物理學(xué)調(diào)劑 +4	小羊36 2026-03-30	4/200	2026-03-31 16:16 by lishahe
[考研] 一志愿浙江大學(xué)工科動(dòng)力工程370,數(shù)一121,專業(yè)課135，現(xiàn)在能去哪里 +3	080700調(diào)劑 2026-03-30	4/200	2026-03-31 12:00 by KLMY666
[考研] 085602化工求調(diào)劑（331分） +8	111@127 2026-03-30	8/400	2026-03-30 21:23 by 研究僧導(dǎo)導(dǎo)
[考研] 085701環(huán)境工程求調(diào)劑 +11	多久上課 2026-03-27	12/600	2026-03-30 21:21 by 研究僧導(dǎo)導(dǎo)
[考研] 297 地理學(xué)070500 復(fù)試求調(diào)劑 +3	小圓圈圈ooo 2026-03-30	3/150	2026-03-30 21:05 by 余震yz
[考研] 085701求調(diào)劑初試286分 +5	secret0328 2026-03-28	5/250	2026-03-30 12:54 by fangnagu
[考研] 348求調(diào)劑 +6	小懶蟲(chóng)不懶了 2026-03-28	6/300	2026-03-30 10:29 by Evan_Liu
[考研] 295求調(diào)劑 +5	wei-5 2026-03-26	5/250	2026-03-30 08:34 by 探123
[考研] 340求調(diào)劑 +6	Amber00 2026-03-26	6/300	2026-03-29 12:06 by 無(wú)際的草原
[考研] 332求92調(diào)劑 +8	蕉蕉123 2026-03-28	8/400	2026-03-29 10:46 by 周梓丹

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時(shí)熱門(mén)版塊排行榜

MJFaber

[求助] 彈性常數(shù)計(jì)算

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

» 猜你喜歡

bcwah717

孫瑩潔

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频