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
quantumespresso:フォノンの状態密度と分散 [2020/11/19 02:49]
koudai [固有振動モードの可視化]
quantumespresso:フォノンの状態密度と分散 [2021/06/27 22:04] (current)
Line 56: Line 56:
 |nq1,nq2,nq3|x,y,z方向のメッシュの数| |nq1,nq2,nq3|x,y,z方向のメッシュの数|
  
-実行は次のようにす+実行は次のようにします。
  
   $ ph.x < Si.ph.in > Si.ph.out   $ ph.x < Si.ph.in > Si.ph.out
  
-フォノンの計算には時間がかかるので、メッシュのとりすぎに注意。+フォノンの計算には時間がかかるので、メッシュのとりすぎに注意します 
 +ただし、メッシュの数による収束は必ず確認すること。 
  
 ===== 実空間へのフーリエ変換 ===== ===== 実空間へのフーリエ変換 =====
Line 169: Line 171:
 ==== 固有振動モードの可視化 ==== ==== 固有振動モードの可視化 ====
  
-  * [[https://www.materialscloud.org/work/tools/interactivephonon|Interactive phonon visualizer]]を使えば、フォノンの各エネルギーに対応する固有振動モードを可視化することができます。 +  * [[http://henriquemiranda.github.io/phononwebsite/|Phonon website]]を使えば、フォノンの各エネルギーに対応する固有振動モードを可視化することができます。
-  * Si.scf.in と matdyn.modes をアップロードしてください。matdyn.modesはmatdyn.xを実行すると自動的に生成されます。 +
-===== 注意 =====+
  
-構造が不安定な系では対応するモードのフォノンの波数が負なる。 +ホームディレクトリに戻りGithubから必要なソフトをダウンロードします 
-もし有限温度で構造相転移があるのであれば、pw.x計算する際に&systemsmearing='fd'して、degauss温度依存性調べとができる。 +<code> 
-1[Ry]=158000[K]なの、例えばdegaus=0.001とれば158K構造が安定か調べることができる。+$ cd 
 +$ git clone https://github.com/henriquemiranda/phononwebsite 
 +</code> 
 + 
 +インストールします 
 +<code> 
 +$ cd phononwebsite 
 +$ python setup.py install --user 
 +</code> 
 + 
 +フォノン計算作業ディレクトリに入り、SCF計算の入力ファイルをprefix.scf、フォノン分散計算で得られたmatdyn.modesをprefix.modesという名前コピーします。 
 +prefix部分はなんも良いです。 
 +   * 結晶の基本並進ベクトル CELL_PARAMETERS を指定してい場合は、ボーア半径指定してください。 
 +   * 原子位置 (ATOMIC_POSITIONS) は分率座標 (crystal) 指定してください。 
 + 
 +そして次のコマンド実行します。 
 +<code> 
 +$ read_qe_phonon.py prefix 
 +</code> 
 +ると prefix.json というファイル得られます。 
 + 
 +さきほどダウンロードしたディレクトリの中にある phononwebsite/phonon.html か、あるいは https://henriquemiranda.github.io/phononwebsite/phonon.html を開き、Custom file の [Browse] で prefix.json を開きま。 
 + 
 +  * エネルギー分散をクリックすること、そのエネルギーに対応する固有振動モードのアニメーションを見ることができます 
 +  * 振動のアニメーションを保存したい場合は Export movie の [gif] をクリックすと録画が開始します 
 +===== 注意 =====
  
-(ただしpw.x計算k点の数はなりたくさん必要なる+構造が不安定な系では対応するモードフォノンの振動数が複素数なります(出力ファイルは負振動で表現されます)。 
 +これ簡単にはフォノンのバネ定数が負にってしまう(つま変形しほうがエネルギー的安定ことに対応します。
  
quantumespresso/フォノンの状態密度と分散.1605721751.txt.gz · Last modified: 2021/06/27 22:00 (external edit)