Differences

This shows you the differences between two versions of the page.

--- quantumespresso:フォノンの状態密度と分散 [2018/11/04 19:34]
koudai 作成
+++ quantumespresso:フォノンの状態密度と分散 [2019/01/31 17:25]
koudai [分散]
@@ Line 9: / Line 9: @@
 <file - Si.scf.in>
 &control
-   calculation='scf',
+  calculation = 'scf',
-   prefix='Si'
+  prefix      = 'Si'
-   pseudo_dir = './',
+  pseudo_dir  = './pseudo',
-   outdir='./tmp/'
+  outdir      = './tmp/'
 /
 &system
-   ibrav=2, celldm(1)=10.2, nat=2, ntyp=1,
+  ibrav     = 2,
-   ecutwfc=18.0
+  celldm(1) = 10.2656,
+  nat       = 2,
+  ntyp      = 1,
+  ecutwfc   = 50.0
 /
 &electrons
+  conv_thr = 1.0d-12
 /
 ATOMIC_SPECIES
- Si  28.086  Si.pz-vbc.UPF
+ Si  28.0855  Si.pbesol-n-kjpaw_psl.1.0.0.UPF
 ATOMIC_POSITIONS (alat)
  Si 0.00 0.00 0.00
@@ Line 36: / Line 40: @@
 <file - Si.ph.in>
 phonons of Si
- &inputph
+&inputph
-  tr2_ph=1.0d-14,
+  tr2_ph = 1.0d-14,
-  prefix='Si',
+  prefix = 'Si',
-  amass(1)=28.08,
+  outdir = './tmp/',
-  outdir='./tmp/',
+  fildyn = 'Si.dyn',
-  fildyn='Si.dyn',
+  ldisp  = .true.,
-  ldisp=.true.,
+  nq1    = 4,
-  nq1=4,
+  nq2    = 4,
-  nq2=4,
+  nq3    = 4,
-  nq3=4,
+/
- /
 </file>
@@ Line 55: / Line 58: @@
 実行は次のようにする。
-  $ ph.x < si.ph.in > si.ph.out
+  $ ph.x < Si.ph.in > Si.ph.out
 フォノンの計算には時間がかかるので注意。
@@ Line 74: / Line 77: @@
 出力先はSi.fcになる。 実行は次のようにする。
-  $ q2r.x < q2r.in > q2r.out
+  $ mpirun q2r.x < q2r.in > q2r.out
 以上で準備は終了である。
@@ Line 87: / Line 90: @@
 <file - matdyn.phdos.in>
 &input
-  flfrc = 'si.fc'
+  flfrc = 'Si.fc'
-  fldos = 'si.phdos'
+  fldos = 'Si.phdos'
   asr = 'crystal'
   dos = .true.
@@ Line 108: / Line 111: @@
 実行は次のようにする。
-  $ matdyn.x < matdyn.phdos.in > matdyn.phdos.out
+  $ mpirun matdyn.x < matdyn.phdos.in > matdyn.phdos.out
 {{:quantumespresso:phdos.png?400|}}
@@ Line 115: / Line 118: @@
 ===== 分散 =====
-分散の場合は計算する波数点の経路を指定する。
+分散の場合は計算する波数点の経路（デカルト座標で、単位は格子定数aを使って2π/a）を指定する。
 <file - matdyn.freq.in>
 &input
   asr = 'crystal'
-  flfrc = 'si.fc'
+  flfrc = 'Si.fc'
-  flfrq = 'si.freq'
+  flfrq = 'Si.freq'
   q_in_band_form=.true.
 /
@@ Line 136: / Line 139: @@
 実行はさきほどと同様。
-  $ matdyn.x < matdyn.freq.in > matdyn.freq.out
+  $ mpirun matdyn.x < matdyn.freq.in > matdyn.freq.out
 結果はSi.freqに出力される。
@@ Line 142: / Line 145: @@
   $ plotband.x
-  Input file > si.freq
+  Input file > Si.freq
   Reading 6 bands at 41 k-points
   Range: 0.0000 510.2482eV Emin, Emax > 0, 550

雑多な記録

User Tools

Site Tools

Differences

Page Tools