雑多な記録

User Tools

Site Tools

Trace: • フォントの設定 • itensor • インストール • インストール
alamode:フォノンの分散と状態密度

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
alamode:フォノンの分散と状態密度 [2021/06/30 02:44]
koudai [概要] 		alamode:フォノンの分散と状態密度 [2021/08/17 02:49] (current)
koudai [フォノンの状態密度]
Line 3: Line 3:
   * Siを例にフォノンの分散と状態密度を計算します。   * Siを例にフォノンの分散と状態密度を計算します。
   * 計算には超格子の作成が必要です。   * 計算には超格子の作成が必要です。
-  * 化合物で見られる極性半導体ではTO-LO分裂というのがあるので、それは別のページで解説予定です。+  * 化合物で見られる極性半導体ではLO-TO分裂というものを考えないといけませんが、それは別のページで解説します。
   * 事前に cif2cell をインストールしてください。   * 事前に cif2cell をインストールしてください。
     * [[QuantumESPRESSO:CIF2Cell]]      * [[QuantumESPRESSO:CIF2Cell]] 
Line 23: Line 23:
 <code> <code>
 $ cif2cell -p pwscf -f 9008565.cif -o Si.pw.in $ cif2cell -p pwscf -f 9008565.cif -o Si.pw.in
-$ cif2cell --supercell=[3,3,3] -p pwscf -f 9008565.cif -o Si333.pw.in+$ cif2cell --supercell=[3,3,3] -p pwscf -f 9008565.cif -o Si_333.pw.in
 </code> </code>
  
Line 32: Line 32:
 また、原子に加わる力の大きさを出力するために tprnfor=.true. を忘れないようにします。 また、原子に加わる力の大きさを出力するために tprnfor=.true. を忘れないようにします。
  
-<file - Si333.pw.in>+<file - Si_333.pw.in>
 &control &control
    calculation='scf',    calculation='scf',
-   prefix='Si',+   prefix='Si_333',
    pseudo_dir = './pseudo/',    pseudo_dir = './pseudo/',
    outdir='./tmp/'    outdir='./tmp/'
Line 78: Line 78:
 超格子の原子位置を少しだけずらしたQEの入力ファイルを作成します。 超格子の原子位置を少しだけずらしたQEの入力ファイルを作成します。
 <nowiki>&cell</nowiki>フィールドと<nowiki>&position</nowiki>フィールドには、さきほどの超格子のインプットファイルで得られた基本並進ベクトルと原子位置をはりつけてください。 <nowiki>&cell</nowiki>フィールドと<nowiki>&position</nowiki>フィールドには、さきほどの超格子のインプットファイルで得られた基本並進ベクトルと原子位置をはりつけてください。
-<file - Si.alm0.in>+<file - Si_333.harmonic.in>
 &general &general
-  PREFIX = Si333+  PREFIX = Si_333
   MODE = suggest   MODE = suggest
   NAT = 54; NKD = 1   NAT = 54; NKD = 1
Line 111: Line 111:
 できたらALMを実行します。 できたらALMを実行します。
 <code> <code>
-$ alm Si.alm0.in > Si.alm0.log+$ alm Si_333.harmonic.in > Si_333.harmonic.log
 </code> </code>
  
-すると Si333.pattern_HARMONIC というファイルができるので、ALAMODEに付属のPythonスクリプトを使ってQEの入力ファイルを作成します。+すると Si_333.pattern_HARMONIC というファイルができるので、ALAMODEに付属のPythonスクリプトを使ってQEの入力ファイルを作成します。
 <code> <code>
-$ python -m displace --QE Si333.pw.in --mag 0.01 --prefix disp -pf Si333.pattern_HARMONIC+$ python -m displace --QE Si_333.pw.in --mag 0.01 --prefix disp -pf Si_333.pattern_HARMONIC
 </code> </code>
   * mag ... 原子をどの程度ずらすかで、0.01から0.04程度の値（単位はオングストローム）が推奨されていますので、いろいろ試してください   * mag ... 原子をどの程度ずらすかで、0.01から0.04程度の値（単位はオングストローム）が推奨されていますので、いろいろ試してください
Line 130: Line 130:
 disp1.pw.outに各原子にはたらく力が出力されていますので、これをALAMODEに付属のPythonスクリプトで整理します。 disp1.pw.outに各原子にはたらく力が出力されていますので、これをALAMODEに付属のPythonスクリプトで整理します。
 <code> <code>
-$ python -m extract --QE Si333.pw.in *.pw.out > DFSET_harmonic+$ python -m extract --QE Si_333.pw.in *.pw.out > DFSET_harmonic
 </code> </code>
  
 最後にALMを使って原子間力定数を計算します。 最後にALMを使って原子間力定数を計算します。
-Si.alm0.inをコピーしてSi.alm1.inを生成し、次の点を変更します +Si_333.harminic.inをコピーしてSi.harmonic_opt.inを生成し、次の点を変更します 
-  * <nowiki>&general</nowiki>フィールドでMODE=optimizeに変更+  * <nowiki>&general</nowiki>フィールドで MODE=optimize に変更
   * <nowiki>&optimize</nowiki>フィールドを作成して、中に DFSET = DFSET_harmonic と記述   * <nowiki>&optimize</nowiki>フィールドを作成して、中に DFSET = DFSET_harmonic と記述
-<file - Si.alm1.in>+<file - Si_333.harmonic_opt.in>
 &general &general
-  PREFIX = Si333+  PREFIX = Si_333
   MODE = optimize   MODE = optimize
   NAT = 54; NKD = 1   NAT = 54; NKD = 1
Line 174: Line 174:
 できたらALMを実行します できたらALMを実行します
 <code> <code>
-$ alm Si.alm1.in > Si.alm1.log+$ alm Si_333.harmonic_opt.in > Si.harmonic_opt.log
 </code> </code>
  
-Si333.fcsとSi333.xmlが出力されたら成功です。+Si_333.fcs と Si_333.xml が出力されたら成功です。
  
  
 ===== フォノンの分散 ===== ===== フォノンの分散 =====
  
-<file - Si_phband.in>+<file - Si_333.phband.in>
 &general &general
-  PREFIX = Si333+  PREFIX = Si_333
   MODE   = phonons   MODE   = phonons
-  FCSXML = Si333.xml+  FCSXML = Si_333.xml
  
   NKD = 1; KD = Si   NKD = 1; KD = Si
Line 211: Line 211:
  
 <code> <code>
-$ anphon Si_phband.in > Si_phband.out+$ anphon Si_333.phband.in > Si_333.phband.log
 </code> </code>
-フォノンの分散は Si333.bands に出力されます+フォノンの分散は Si_333.bands に出力されます
  
 {{ :alamode:bands.png?400 |}} {{ :alamode:bands.png?400 |}}
Line 221: Line 221:
 ===== フォノンの状態密度 ===== ===== フォノンの状態密度 =====
  
-Si_phband.in をコピーして si_phdos.in というファイルを作成します+Si_333.phband.in をコピーして Si_333.phdos.in というファイルを作成します
   * <nowiki>&general</nowiki>フィールドでDELTA_E（エネルギーの刻み幅。単位はcm^-1。デフォルトは10でかなり粗い）を指定します。   * <nowiki>&general</nowiki>フィールドでDELTA_E（エネルギーの刻み幅。単位はcm^-1。デフォルトは10でかなり粗い）を指定します。
   * <nowiki>&kpoint</nowiki>フィールドを状態密度の計算用に書き換えます。   * <nowiki>&kpoint</nowiki>フィールドを状態密度の計算用に書き換えます。
  
-<file - Si_phdos.in>+<file - Si_333.phdos.in>
 &general &general
-  PREFIX = Si333+  PREFIX = Si_333
   MODE   = phonons   MODE   = phonons
-  FCSXML = Si333.xml+  FCSXML = Si_333.xml
  
   NKD = 1; KD = Si   NKD = 1; KD = Si
Line 253: Line 253:
  
 <code> <code>
-$ anphon Si_phdos.in > Si_phdos.out+$ anphon Si_333.phdos.in > Si_333.phdos.log
 </code> </code>
  
-状態密度は Si333.dos に出力されます+状態密度は Si_333.dos に出力されます
  
 {{ :alamode:dos.png?400 |}} {{ :alamode:dos.png?400 |}}
  
-ついでに Si333.thermo に格子の熱力学量（比熱、エントロピー、内部エネルギー、自由エネルギー）の温度依存性も出力されます。+フォノンを自由ボソンとしたときの熱力学量（比熱、エントロピー、内部エネルギー、自由エネルギー）の温度依存性も Si_333.thermo に出力されます。
 温度の刻みはデフォルトで DT = 10 (K)ですので、より細かく取りたい場合は<nowiki>&general</nowiki>フィールドで DT の値を指定してください。 温度の刻みはデフォルトで DT = 10 (K)ですので、より細かく取りたい場合は<nowiki>&general</nowiki>フィールドで DT の値を指定してください。
alamode/フォノンの分散と状態密度.1624988676.txt.gz · Last modified: 2021/06/30 02:44 by koudai

Page Tools

Except where otherwise noted, content on this wiki is licensed under the following license: CC Attribution-Share Alike 4.0 International
CC Attribution-Share Alike 4.0 International Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki