User Tools

Site Tools


quantumespresso:フォノンの状態密度と分散

Differences

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

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
Next revision Both sides next revision
quantumespresso:フォノンの状態密度と分散 [2019/01/07 19:51]
koudai [電子状態計算]
quantumespresso:フォノンの状態密度と分散 [2019/01/07 21:02]
koudai [分散]
Line 9: Line 9:
 <file - Si.scf.in> <file - Si.scf.in>
 &control &control
-  calculation='scf', +  calculation = 'scf', 
-  prefix='Si' +  prefix      = 'Si' 
-  pseudo_dir = './pseudo', +  pseudo_dir  = './pseudo', 
-  outdir='./tmp/'+  outdir      = './tmp/'
 / /
 &system &system
Line 22: Line 22:
 / /
 &electrons &electrons
-  conv_thr         = 1.0d-12+  conv_thr = 1.0d-12
 / /
 ATOMIC_SPECIES ATOMIC_SPECIES
Line 40: Line 40:
 <file - Si.ph.in> <file - Si.ph.in>
 phonons of Si 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> </file>
  
Line 59: Line 58:
 実行は次のようにする。 実行は次のようにする。
  
-  $ ph.x < si.ph.in > si.ph.out+  $ ph.x < Si.ph.in > Si.ph.out
  
 フォノンの計算には時間がかかるので注意。 フォノンの計算には時間がかかるので注意。
Line 78: Line 77:
 出力先はSi.fcになる。 実行は次のようにする。 出力先はSi.fcになる。 実行は次のようにする。
  
-  $ q2r.x < q2r.in > q2r.out+  $ mpirun q2r.x < q2r.in > q2r.out
  
 以上で準備は終了である。 以上で準備は終了である。
Line 91: 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 = 'crystal'   asr = 'crystal'
   dos = .true.   dos = .true.
Line 112: Line 111:
 実行は次のようにする。 実行は次のようにする。
  
-  $ matdyn.x < matdyn.phdos.in > matdyn.phdos.out+  $ mpirun matdyn.x < matdyn.phdos.in > matdyn.phdos.out
  
 {{:quantumespresso:phdos.png?400|}} {{:quantumespresso:phdos.png?400|}}
Line 124: Line 123:
 &input &input
   asr = 'crystal'   asr = 'crystal'
-  flfrc = 'si.fc' +  flfrc = 'Si.fc' 
-  flfrq = 'si.freq'+  flfrq = 'Si.freq'
   q_in_band_form=.true.   q_in_band_form=.true.
 / /
Line 140: Line 139:
 実行はさきほどと同様。 実行はさきほどと同様。
  
-  $ matdyn.x < matdyn.freq.in > matdyn.freq.out+  $ mpirun matdyn.x < matdyn.freq.in > matdyn.freq.out
  
 結果はSi.freqに出力される。 結果はSi.freqに出力される。
Line 146: Line 145:
  
   $ 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
quantumespresso/フォノンの状態密度と分散.txt · Last modified: 2021/06/27 22:04 (external edit)