This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision Next revision Both sides next revision | ||
quantumespresso:フォノンの状態密度と分散 [2019/01/07 19:49] koudai [電子状態計算] |
quantumespresso:フォノンの状態密度と分散 [2019/06/26 21:09] koudai [状態密度] |
||
---|---|---|---|
Line 9: | Line 9: | ||
<file - Si.scf.in> | <file - Si.scf.in> | ||
& | & | ||
- | calculation=' | + | |
- | | + | prefix |
- | | + | pseudo_dir |
- | | + | outdir |
/ | / | ||
&system | &system | ||
- | ibrav = 2, | + | |
- | | + | celldm(1) = 10.2656, |
- | | + | nat = 2, |
- | | + | ntyp = 1, |
- | | + | ecutwfc |
/ | / | ||
& | & | ||
+ | conv_thr = 1.0d-12 | ||
/ | / | ||
ATOMIC_SPECIES | ATOMIC_SPECIES | ||
Line 39: | Line 40: | ||
<file - Si.ph.in> | <file - Si.ph.in> | ||
phonons of Si | phonons of Si | ||
- | & | + | & |
- | tr2_ph=1.0d-14, | + | tr2_ph = 1.0d-14, |
- | prefix=' | + | prefix = ' |
- | amass(1)=28.08, | + | outdir = ' |
- | outdir=' | + | fildyn = ' |
- | fildyn=' | + | ldisp = .true., |
- | ldisp=.true., | + | nq1 = 4, |
- | nq1=4, | + | nq2 = 4, |
- | nq2=4, | + | nq3 = 4, |
- | nq3=4, | + | / |
- | / | + | |
</ | </ | ||
Line 58: | Line 58: | ||
実行は次のようにする。 | 実行は次のようにする。 | ||
- | $ ph.x < si.ph.in > si.ph.out | + | $ ph.x < Si.ph.in > Si.ph.out |
フォノンの計算には時間がかかるので注意。 | フォノンの計算には時間がかかるので注意。 | ||
Line 90: | Line 90: | ||
<file - matdyn.phdos.in> | <file - matdyn.phdos.in> | ||
&input | &input | ||
- | flfrc = 'si.fc' | + | flfrc = 'Si.fc' |
- | fldos = 'si.phdos' | + | fldos = 'Si.phdos' |
asr = ' | asr = ' | ||
dos = .true. | dos = .true. | ||
Line 112: | Line 112: | ||
$ matdyn.x < matdyn.phdos.in > matdyn.phdos.out | $ matdyn.x < matdyn.phdos.in > matdyn.phdos.out | ||
+ | |||
+ | 下図は Si.phdos に出力されたデータをプロットした。 | ||
{{: | {{: | ||
Line 118: | Line 120: | ||
===== 分散 ===== | ===== 分散 ===== | ||
- | 分散の場合は計算する波数点の経路を指定する。 | + | 分散の場合は計算する波数点の経路(デカルト座標で、単位は格子定数aを使って2π/ |
<file - matdyn.freq.in> | <file - matdyn.freq.in> | ||
&input | &input | ||
asr = ' | asr = ' | ||
- | flfrc = 'si.fc' | + | flfrc = 'Si.fc' |
- | flfrq = 'si.freq' | + | flfrq = 'Si.freq' |
q_in_band_form=.true. | q_in_band_form=.true. | ||
/ | / | ||
Line 145: | Line 147: | ||
$ plotband.x | $ plotband.x | ||
- | Input file > si.freq | + | Input file > Si.freq |
Reading 6 bands at 41 k-points | Reading 6 bands at 41 k-points | ||
Range: 0.0000 510.2482eV Emin, Emax > 0, 550 | Range: 0.0000 510.2482eV Emin, Emax > 0, 550 |