Difference between revisions of "エクササイス"

From salmon
Jump to: navigation, search
(→‎Exercise-3: Electron dynamics in C2H2 molecule under a pulsed electric field)
 
(45 intermediate revisions by the same user not shown)
Line 9: Line 9:
  
 
[[インストールと実行]]で述べたように、SALMONを実行するためには少なくとも
 
[[インストールと実行]]で述べたように、SALMONを実行するためには少なくとも
−
入力ファイルと擬ポテンシャルファイルが必要です。
+
インプットファイルと擬ポテンシャルファイルが必要です。
−
以下では、いくつかのサンプル計算に対する入力ファイルを示し、その入力ファイルに
+
以下では、いくつかのサンプル計算に対するインプットファイルを示し、そのインプットファイルに
 
現れるネームリスト変数の簡単な説明を行います。
 
現れるネームリスト変数の簡単な説明を行います。
−
その入力ファイルを修正して、あなた自身が行いたい計算を実行することができるでしょう。
+
そのインプットファイルを修正して、あなた自身が行いたい計算を実行することができるでしょう。
 
サンプルに含まれている元素の擬ポテンシャルも提供されています。
 
サンプルに含まれている元素の擬ポテンシャルも提供されています。
 
主要な出力ファイルの説明も行います。
 
主要な出力ファイルの説明も行います。
Line 19: Line 19:
  
 
最初の３つのエクササイス（エクササイス-1~3）は、孤立したアセチレン分子C2H2に対するものです。
 
最初の３つのエクササイス（エクササイス-1~3）は、孤立したアセチレン分子C2H2に対するものです。
−
もし孤立した系の電子ダイナミクス計算に興味がある場合は、これらのエクササイスをごらんください。
+
もし孤立した系の電子ダイナミクス計算に興味がある場合は、これらのエクササイスをご覧ください。
 
SALMONでは、通常最初に基底状態計算を行います。
 
SALMONでは、通常最初に基底状態計算を行います。
−
これは、[[#エクササイス-1:C2H2分子の基底状態計算|エクササイス-1]]で説明します。
+
これは、[[#エクササイス-1: C2H2分子の基底状態|エクササイス-1]]で説明します。
 
基底状態計算が終わったのち、２つの電子ダイナミクス計算が用意されています。
 
基底状態計算が終わったのち、２つの電子ダイナミクス計算が用意されています。
−
[[#エクササイス-2:C2H2分子の文曲率と光吸収|エクササイス-2]]では、分子の分極率と光吸収を得る
+
[[#エクササイス-2: C2H2分子の分極率と光吸収|エクササイス-2]]では、分子の分極率と光吸収を得る
 
実時間での線形光応答計算を示します。
 
実時間での線形光応答計算を示します。
−
[[#エクササイス-3:パルス電場中にあるC2H2分子の電子ダイナミクス|エクササイス-3]]では、
+
[[#エクササイス-3: パルス電場を印加したC2H2分子における電子ダイナミクス|エクササイス-3]]では、
−
パルス電場中にある分子の電子ダイナミクス計算を示します。
+
パルス電場が印加された分子における電子ダイナミクス計算を示します。
  
 
次の２つのエクササイス(エクササイス-4~5)は、結晶固体シリコンに対するものです。
 
次の２つのエクササイス(エクササイス-4~5)は、結晶固体シリコンに対するものです。
−
もし無限周期系における電子ダイナミクス計算を学ぶことに興味があれば、これらのエクササイスを見てください。
+
もし無限周期系における電子ダイナミクス計算を学ぶことに興味があれば、これらのエクササイスをご覧ください。
−
小さな単位セルを持つ系の基底状態計算は容易であり、時間発展計算の方がはるかに計算時間がかかるため、
+
周期系の計算では、小さな単位セルを持つ系の基底状態計算は容易であり、時間発展計算の方がはるかに計算時間がかかるため、
 
基底状態と時間発展の計算を一つのジョブとして実行することを推奨します。
 
基底状態と時間発展の計算を一つのジョブとして実行することを推奨します。
 
以下の２つのエクササイスも、そのようになっています。
 
以下の２つのエクササイスも、そのようになっています。
−
[[#エクササイス-4:結晶シリコンの誘電関数|エクササイス-4]]では、結晶シリコンの誘電関数を得るための
+
[[#エクササイス-4: シリコン結晶の誘電関数|エクササイス-4]]では、結晶シリコンの誘電関数を得るための
 
線形応答性質の計算を示します。
 
線形応答性質の計算を示します。
−
[[#エクササイス-5:パルス電場中にある結晶シリコンの電子ダイナミクス|エクササイス5]]
+
[[#エクササイス-5: パルス電場を印加したシリコン結晶における電子ダイナミクス|エクササイス-5]]
−
では、パルス電場により結晶シリコンに誘起される電子ダイナミクス計算を示します。
+
では、パルス電場を印加した結晶シリコンにおける電子ダイナミクス計算を示します。
  
−
最後のエクササイス(エクササイス-6)は、バルクシリコンにおけるパルス光の入射と伝播に対して、
+
最後のエクササイス(エクササイス-6)は、バルクシリコンへのパルス光の入射と伝播に対して、
 
パルス光の電磁場に対するマクスウェル方程式と単位セルの電子ダイナミクス計算を結合したものです。
 
パルス光の電磁場に対するマクスウェル方程式と単位セルの電子ダイナミクス計算を結合したものです。
−
ものです。
 
 
これは大きな計算時間を必要とするため、大規模並列スパコンの利用が推奨されます。
 
これは大きな計算時間を必要とするため、大規模並列スパコンの利用が推奨されます。
−
[[#エクササイス-6:シリコン薄膜中のパルス光の伝播|エクササイス-6]]で、
+
[[#エクササイス-6: パルス光のシリコン薄膜の伝播|エクササイス-6]]で、
 
バルクシリコンの表面に垂直に入射する直線偏光パルス光の入射計算を示します。
 
バルクシリコンの表面に垂直に入射する直線偏光パルス光の入射計算を示します。
  
Line 48: Line 47:
  
 
=== エクササイス-1: C2H2分子の基底状態 ===
 
=== エクササイス-1: C2H2分子の基底状態 ===
−
このエクササイスでは、静的コーン・シャム方程式を解くことにより、アセチレン分子(C2H2)の基底状態計算を学びます。
+
このエクササイスでは、静的コーン・シャム方程式を解くことによりアセチレン分子(C2H2)の基底状態を求める計算を学びます。
−
このエクササイスは、分子やナノ粒子のような孤立した系に対して、SALMONを用いた計算を行う準備をどのように
+
このエクササイスは、分子やナノ粒子のような孤立した系に対して、SALMONを用いた計算をどのようにすれば良いかを知るのに有用です。
−
すれば良いかを知るのに有用です。今のところ、SALMONで構造最適化を実行することはできないことに注意して
+
今のところ、SALMONで構造最適化を実行することはできないことに注意してください。
−
ください。分子の原子位置は入力ファイルて指定され、計算の過程では固定されます。
+
分子の原子位置はインプットファイルて指定され、計算の過程では固定されます。
  
 
==== 入力ファイル ====
 
==== 入力ファイル ====
−
: コードを実行するために、以下のファイルが必要です。
+
コードを実行するために、以下のファイルが必要です。
  
 
{| class="wikitable"
 
{| class="wikitable"
 
| file name      || description
 
| file name      || description
 
|-
 
|-
−
| ''C2H2_gs.inp'' || ネームリスト変数とそれらの値を含む入力ファイル          
+
| ''C2H2_gs.inp'' || ネームリスト変数とそれらの値を含むインプットファイル          
 
|-
 
|-
 
| ''C_rps.dat'' || 炭素原子の擬ポテンシャル
 
| ''C_rps.dat'' || 炭素原子の擬ポテンシャル
Line 66: Line 65:
 
|}
 
|}
  
−
* 上記の３つのファイル（圧縮ファイル）は、以下からダウンロードできます:
+
上記の３つのファイル（圧縮ファイル）は、以下からダウンロードできます:
  
−
  [[media:C2H2_gs_input.zip| 圧縮された入力ファイルと擬ポテンシャルファイルをダウンロード]]
+
  [[media:C2H2_gs_input.zip|圧縮されたインプットファイルと擬ポテンシャルファイルをダウンロード]]
  
−
* 入力ファイル''C2H2_gs.inp''では、ネームリスト変数が指定されています。それらの多くは、
+
インプットファイル''C2H2_gs.inp''では、ネームリスト変数が指定されています。それらの多くは、基底状態計算を実行するために必須です。それらの説明を以下に示します:
−
基底状態計算を実行するために必須です。それらの説明を以下に示します:
 
  
−
  [[入力ファイルの説明(C2H2分子の基底状態)]]
+
  [[インプットファイルの説明（C2H2分子の基底状態）]]
  
−
: あなたが計算をしたい他の系のための入力ファイルを準備するのに、この情報は有用です。
+
あなたが計算をしたい他の系のためのインプットファイルを準備するのに、この情報は有用です。
−
入力ファイルで用いることのできるネームリスト変数の完全なリストは、ダウンロードファイルの''SALMON/manual/input_variables.md''にあります。
+
インプットファイルで用いることのできるネームリスト変数の完全なリストは、ダウンロードファイルの''SALMON/manual/input_variables.md''にあります。
  
 
==== 出力ファイル ====
 
==== 出力ファイル ====
−
: 計算を実行した後、あなたがコードを実行したディレクトリに以下の出力ファイルが生成されています。
+
計算の終了後、あなたが計算を実行したディレクトリに以下の出力ファイルが生成されています。
  
 
{| class="wikitable"
 
{| class="wikitable"
Line 99: Line 97:
 
|}
 
|}
  
−
* 上記のファイル（バイナリファイル''C2H2_gs.bin''を除き、圧縮ファイル)は、以下からダウンロードすることができます:
+
上記のファイル（バイナリファイル''C2H2_gs.bin''を除いた圧縮ファイル）は、以下からダウンロードすることができます:
  
 
  [[media:C2H2_gs_output.zip|圧縮された出力ファイルのダウンロード]]
 
  [[media:C2H2_gs_output.zip|圧縮された出力ファイルのダウンロード]]
  
−
* 軌道エネルギーのような主要な計算結果は''C2H2_info.data''に含まれています。
+
軌道エネルギーのような主要な計算結果は''C2H2_info.data''に含まれています。''C2H2_info.data''と他の出力ファイルの内容を以下で説明します:
−
''C2H2_info.data''と他の出力ファイルは、以下で説明されています:
 
  
−
  [[出力ファイル(C2H2分子の基底状態)の説明]]
+
  [[出力ファイルの説明（C2H2分子の基底状態）]]
  
 
==== 画像イメージ ====
 
==== 画像イメージ ====
−
: 出力ファイルから作成されたいくつかの画像イメージを示します。
+
出力ファイルから作成されたいくつかの画像イメージを示します。
  
 
{| class="wikitable"
 
{| class="wikitable"
Line 121: Line 118:
 
|}
 
|}
  
−
===エクササイス-2: C2H2分子の分極率と光吸収 ===
+
=== エクササイス-2: C2H2分子の分極率と光吸収 ===
−
このエクササイスでは、時間依存コーン・シャム方程式を解いてアセチレン(C2H2)分子の線形応答計算を学びます。
+
このエクササイスでは、時間依存コーン・シャム方程式を解いてアセチレン(C2H2)分子の線形応答を調べる計算を学びます。
 
線形応答計算は、分子の分極率や振動子強度分布を与えます。
 
線形応答計算は、分子の分極率や振動子強度分布を与えます。
 
このエクササイスは、 [[#エクササイス-1: C2H2分子の基底状態|エクササイス-1]]で説明された基底状態計算を終えた後に実行する必要があります。
 
このエクササイスは、 [[#エクササイス-1: C2H2分子の基底状態|エクササイス-1]]で説明された基底状態計算を終えた後に実行する必要があります。
Line 131: Line 128:
 
振動数依存分極率の虚部は、振動子強度分布と光吸収断面積に比例します。
 
振動数依存分極率の虚部は、振動子強度分布と光吸収断面積に比例します。
  
−
==== 入力ファイル ====
+
==== インプットファイル ====
−
: 計算を実行するのに、線形応答計算のためのネームリスト変数とそれらの値を含む入力ファイル''C2H2_rt_response.inp''が必要です。
+
計算を実行するのに、線形応答計算のためのネームリスト変数とそれらの値を含むインプットファイル''C2H2_rt_response.inp''が必要です。
 
基底状態計算で生成したバイナリファイル''C2H2_gs.bin''と擬ポテンシャルファイルも必要です。
 
基底状態計算で生成したバイナリファイル''C2H2_gs.bin''と擬ポテンシャルファイルも必要です。
−
擬ポテンシャルファイルは、基底状態計算で用いたものと同じであることが必要です。
+
擬ポテンシャルファイルは、基底状態計算で用いたものと同じものを用いることが必要です。
  
 
{| class="wikitable"
 
{| class="wikitable"
 
| ファイル名      || 説明
 
| ファイル名      || 説明
 
|-
 
|-
−
| ''C2H2_rt_response.inp'' || ネームリスト変数とそれらの値を含む入力ファイル          
+
| ''C2H2_rt_response.inp'' || ネームリスト変数とそれらの値を含むインプットファイル          
 
|-
 
|-
 
| ''C_rps.dat'' || 炭素原子の擬ポテンシャル
 
| ''C_rps.dat'' || 炭素原子の擬ポテンシャル
Line 148: Line 145:
 
|}
 
|}
  
−
* 圧縮された''C2H2_rt_response.inp''ファイルを以下からダウンロードできます:
+
圧縮された''C2H2_rt_response.inp''ファイルを以下からダウンロードできます:
 +
 
 +
[[media:C2H2_rt_response_input.zip|圧縮されたインプットファイルのダウンロード]]
  
−
[[media:C2H2_rt_response_input.zip|圧縮された入力ファイルのダウンロード]]
+
インプットファイル''C2H2_rt_response.inp''では、ネームリスト変数が指定されています。それらの多くは、線形応答計算を実行するために必須です。それらの説明を以下に示します。
  
−
*入力ファイル''C2H2_rt_response.inp''において、ネームリスト変数が指定されます。それらのほとんどは線形応答計算を行うため必須です。それらの説明を以下で行います。
+
[[インプットファイルの説明（C2H2分子の分極率と光吸収）]]
  
−
[[入力ファイルの説明(C2H2分子の分極率と光吸収)]]
+
あなたが計算をしたい他の系のためのインプットファイルを準備するのに、この情報は有用です。 インプットファイルで用いることのできるネームリスト変数の完全なリストは、ダウンロードファイルのSALMON/manual/input_variables.mdにあります。
−
 
−
: これは、あなたが計算を行いたい他の系の入力ファイルを準備するのに役立つでしょう。入力ファイルで用いることができるネームリスト変数の完全なリストは、ダウンロードファイルの''SALMON/manual/input_variables.md''にあります。
 
  
 
==== 出力ファイル ====
 
==== 出力ファイル ====
−
: 計算の実行後、コードを実行したディレクトリに以下の出力ファイルが生成されます。
+
計算の終了後、あなたが計算を実行したディレクトリに以下の出力ファイルが生成されています。
  
 
{| class="wikitable"
 
{| class="wikitable"
Line 169: Line 166:
 
|}
 
|}
  
−
* 上記の圧縮ファイルを以下からダウンロードできます:
+
上記の圧縮ファイルを以下からダウンロードすることができます:
  
 
  [[media:C2H2_rt_response_output.zip|圧縮した出力ファイルのダウンロード]]
 
  [[media:C2H2_rt_response_output.zip|圧縮した出力ファイルのダウンロード]]
  
−
*出力ファイルの説明が以下に与えられています:
+
出力ファイルの内容を以下で説明します:
  
−
  [[出力ファイルの説明(C2H2分子の分極率と光吸収)]]
+
  [[出力ファイルの説明（C2H2分子の分極率と光吸収）]]
  
 
=== エクササイス-3: パルス電場を印加したC2H2分子における電子ダイナミクス ===
 
=== エクササイス-3: パルス電場を印加したC2H2分子における電子ダイナミクス ===
−
このエクササイスでは、時間依存コーン・シャム方程式を解くことにより、パルス電場を印加したアセチレン(C2H2)分子における電子ダイナミクスの計算を学ぶ。
+
 
−
計算の出力として、時間の関数としての全エネルギーや電気双極子モーメントなどの量が計算される。
+
このエクササイスでは、時間依存コーン・シャム方程式を解いてパルス電場を印加したアセチレン(C2H2)分子に起こる電子ダイナミクスの計算を学びます。
−
このエクササイスは、[[#エクササイス-1: C2H2分子の基底状態|エクササイス-1]]で説明した基底状態計算を終えた後に行う必要がある。
+
計算の出力として、時間の関数としての全エネルギーや電気双極子モーメントなどの量が得られます。
−
計算では、cos^2の包絡形を持つパルス電場が加えられる。パルス場を指定する強度、振動数、偏極の向き、キャリアエンベロープ位相などのパラメータは、入力ファイルにおいて指定される。
+
このエクササイスは、[[#エクササイス-1: C2H2分子の基底状態|エクササイス-1]]で説明された基底状態計算を終えた後に実行する必要があります。
 +
このエクササイスでは、cos^2の包絡形を持つパルス電場を用いています。パルス電場を特徴付ける強度、振動数、偏極の向き、キャリアエンベロープ位相などのパラメータは、インプットファイルで指定します。
  
 
==== 入力ファイル ====
 
==== 入力ファイル ====
−
: 計算を実行するために、以下のファイルを用いる。''C2H2.data''ファイルは、基底状態計算で生成される。
+
計算を実行するのに、以下のファイルを用います。''C2H2_gs.bin''は、基底状態計算で生成されたバイナリファイルです。
−
擬ポテンシャルファイルは、基底状態計算で既に用いたものである。従って、パルス電場計算のネームリスト変数とその値を指定する''C2H2_rt_pulse.inp''が、ユーザが用意する必要のある唯一のファイルである。
+
擬ポテンシャルファイルは、基底状態計算で用いたものと同じものを持ちることが必要です。
 +
従って、パルス電場計算のネームリスト変数とその値を指定する''C2H2_rt_pulse.inp''が、ユーザが用意する必要のある唯一のファイルです。
  
 
{| class="wikitable"
 
{| class="wikitable"
 
| ファイル名      || 説明
 
| ファイル名      || 説明
 
|-
 
|-
−
| ''C2H2_rt_pulse.inp'' || ネームリスト変数とそれらの値を含む入力ファイル          
+
| ''C2H2_rt_pulse.inp'' || ネームリスト変数とそれらの値を含むインプットファイル          
 
|-
 
|-
−
| ''C_rps.dat'' || pseodupotential file for Carbon
+
| ''C_rps.dat'' || 炭素原子に対する擬ポテンシャルファイル
 
|-
 
|-
−
| ''H_rps.dat'' || pseudopotential file for Hydrogen
+
| ''H_rps.dat'' || 水素原子に対する擬ポテンシャルファイル
 
|-
 
|-
−
| ''C2H2.data'' || binary file created in the ground state calculation
+
| ''C2H2_gs.bin'' || 基底状態計算で生成されたバイナリファイル
 
|}
 
|}
  
−
* You may download the ''C2H2_rt_pulse.inp'' file (zipped file) from:
+
圧縮された''C2H2_rt_pulse.inp''ファイルを以下からダウンロードできます:
  
−
  [[media:C2H2_rt_pulse_input.zip| Download zipped input file]]
+
  [[media:C2H2_rt_pulse_input.zip|圧縮されたインプットファイルのダウンロード]]
  
−
* In the input file ''C2H2_rt_pulse.inp'', namelists variables are specified. Most of them are mandatory to execute the calculation of electron dynamics induced by a pulsed electric field. We present explanations of the namelist variables that appear in the input file in:
+
インプットファイル''C2H2_rt_pulse.inp''では、ネームリスト変数が指定されています。
 +
それらの多くは、パルス電場により起こる電子ダイナミクス計算を実行するために必須です。
 +
それらの説明を以下に示します:
  
−
  [[Explanations of input files (C2H2 molecule under a pulsed electric field)]]
+
  [[インプットファイルの説明（パルス電場のもとでのC2H2分子）]]
−
 
−
: This will help you to prepare the input file for other systems and other pulsed electric fields that you want to calculate. A complete list of the namelist variables that can be used in the input file can be found in the downloaded file ''SALMON/manual/input_variables.md''.
 
  
−
==== Output files ====
+
あなたが計算をしたい他の系のためのインプットファイルを準備するのに、この情報は有用です。 インプットファイルで用いることのできるネームリスト変数の完全なリストは、ダウンロードファイルのSALMON/manual/input_variables.mdにあります。
−
: After the calculation, following output files are created in the directory that you run the code,
+
 
 +
==== 出力ファイル ====
 +
計算の終了後、あなたが計算を実行したディレクトリに以下の出力ファイルが生成されています。
  
 
{| class="wikitable"
 
{| class="wikitable"
−
| file name       || description
+
| ファイル名       || 説明
 
|-
 
|-
−
| ''C2H2_p.data'' || components of the electric dipole moment as functions of time              
+
| ''C2H2_p.data'' || 時間の関数としての電気双極子モーメントの成分              
 
|-
 
|-
−
| ''C2H2_ps.data'' || power spectrum that is obtained by a time-frequency Fourier transformation of the electric dipole moment
+
| ''C2H2_ps.data'' || 電気双極子モーメントを時間-振動数フーリェ変換することで得られるパワースペクトル
 
|}
 
|}
  
−
* You may download the above files (zipped file) from:
+
圧縮された上記のファイルを以下からダウンロードすることができます:
  
−
  [[media:C2H2_rt_pulse_output.zip|Download zipped output files]]
+
  [[media:C2H2_rt_pulse_output.zip|圧縮された出力ファイルをダウンロード]]
  
−
*Explanations of the files are described in:
+
出力ファイルの内容を以下で説明します:
  
−
  [[Explanations of output files (C2H2 molecule under a pulsed electric field)]]
+
  [[出力ファイルの説明（パルス電場のもとでのC2H2分子）]]
  
−
== Crystalline silicon (periodic solids) ==
+
== シリコン結晶(周期的な固体) ==
  
−
=== Exercise-4: Dielectric function of crystalline silicon===
+
=== エクササイス-4: シリコン結晶の誘電関数 ===
−
In this exercise, we learn the linear response calculation of the crystalline silicon of a diamond structure.
+
このエクササイスでは、ダイヤモンド構造を持つシリコン結晶の線形応答計算を学びます。
−
Calculation is done in a cubic unit cell that contains eight silicon atoms.
+
計算は、８個のシリコン原子を含む立方体の単位セルで行います。
−
Since the ground state calculation costs much less computational time than the time evolution calculation,
+
時間発展計算に比べて基底状態計算は計算時間がはるかに少ないため、両方の計算を連続して実行します。
−
both calculations are successively executed.
+
基底状態計算の終了後、単位セルの全ての電子に対して''z''方向に撃力の摂動を加えます。
−
After finishing the ground state calculation, an impulsive perturbation is applied to all electrons in the unit cell
+
ダイヤモンド構造では誘電関数は等方的なため、計算された誘電関数は摂動の方向には依存しません。
−
along ''z'' direction.
+
時間発展の間、単位セルの体積で平均を取った電子の流れの密度が計算されます。
−
Since the dielectric function is isotropic in the diamond structure, calculated dielectric function should not depend
+
電子の流れの密度に対して時間と振動数のフーリェ変換を行うと、振動数依存伝導度が得られます。
−
on the direction of the perturbation.
+
誘電関数は伝導度から通常の関係に従い得られます。
−
During the time evolution, electric current averaged over the unit cell volume is calculated.
 
−
A time-frequency Fourier transformation of the electric current gives us a frequency-dependent conductivity.
 
−
The dielectric function may be obtained from the conductivity using a standard relation.
 
  
−
==== Input files ====
+
==== 入力ファイル ====
−
: To run the code, following files are used:
+
計算を実行するのに、以下のファイルを用います:
  
 
{| class="wikitable"
 
{| class="wikitable"
−
| file name       || description
+
| ファイル名       || 説明
 
|-
 
|-
−
| ''Si_gs_rt_response.inp'' || input file that contain namelist variables and their values.         
+
| ''Si_gs_rt_response.inp'' || ネームリスト変数とそれらの値を含むインプットファイル       
 
|-
 
|-
−
| ''Si_rps.dat'' || pseodupotential file of silicon
+
| ''Si_rps.dat'' || シリコン原子の擬ポテンシャル
 
|}
 
|}
  
−
* You may download the above 2 files (zipped file) from:
+
圧縮された上記の２つのファイルを以下からダウンロードできます:
  
−
  [[media: Si_gs_rt_response_input.zip| Download zipped input and pseudopotential files]]
+
  [[media: Si_gs_rt_response_input.zip|圧縮されたインプットファイルと擬ポテンシャルファイルのダウンロード]]
  
−
* In the input file ''Si_gs_rt_response.inp'', namelists variables are specified. Most of them are mandatory to execute the calculation. We present explanations of the namelist variables that appear in the input file in:
+
インプットファイル''Si_gs_rt_response.inp''では、ネームリスト変数が指定されています。それらの多くは、計算を行うために必須です。
 +
それらの説明を以下に示します:
  
−
  [[Explanations of input files (dielectric function of crystalline silicon)]]
+
  [[インプットファイルの説明（シリコン結晶の誘電関数）]]
  
−
+
あなたが計算をしたい他の系のためのインプットファイルを準備するのに、この情報は有用です。 インプットファイルで用いることのできるネームリスト変数の完全なリストは、ダウンロードファイルのSALMON/manual/input_variables.mdにあります。
−
: This will help you to prepare the input file for other systems that you want to calculate. A complete list of the namelist variables that can be used in the input file can be found in the downloaded file ''SALMON/manual/input_variables.md''.
 
  
−
==== Output files ====
+
==== 出力ファイル ====
−
: After the calculation, following output files are created in the directory that you run the code,
+
計算の終了後、あなたが計算を実行したディレクトリに以下の出力ファイルが生成されています。
  
 
{| class="wikitable"
 
{| class="wikitable"
−
| file name       || description
+
| ファイル名       || 説明
 
|-
 
|-
−
| ''Si_gs_info.data'' || information of ground state calculation
+
| ''Si_gs_info.data'' || 基底状態計算の情報
 
|-
 
|-
−
| ''Si_eigen.data'' || energy eigenvalues of orbitals
+
| ''Si_eigen.data'' || 軌道のエネルギー固有値
 
|-
 
|-
−
| ''Si_k.data'' || information on k-points
+
| ''Si_k.data'' || k点に関する情報
 
|-
 
|-
−
| ''Si_rt.data'' || electric field, vector potential, and current as functions of time
+
| ''Si_rt.data'' || 時間の関数としての電場、ベクトルポテンシャル、カレント
 
|-
 
|-
−
| ''Si_force.data'' || force acting on atoms
+
| ''Si_force.data'' || 原子に働く力
 
|-
 
|-
−
| ''Si_lr.data'' || Fourier spectra of the dielectric functions
+
| ''Si_lr.data'' || 誘電関数のフーリェスペクトル
 
|-
 
|-
−
| ''Si_gs_rt_response.out'' || standard output file
+
| ''Si_gs_rt_response.out'' || 標準出力ファイル
 
|}
 
|}
  
−
* You may download the above files (zipped file) from:
+
圧縮された上記のファイルを以下からダウンロードすることができます:
  
−
  [[media:Si_gs_rt_response_output.zip|Download zipped output files]]
+
  [[media:Si_gs_rt_response_output.zip|圧縮した出力ファイルをダウンロード]]
  
−
* Explanations of the output files are described in:
+
出力ファイルの内容を以下で説明します:
  
−
  [[Explanation of output fiels (dielectric function of crystalline silicon)]]
+
  [[出力ファイルの説明（シリコン結晶の誘電関数）]]
  
−
=== Exercise-5: Electron dynamics in crystalline silicon under a pulsed electric field ===
+
=== エクササイス-5: パルス電場を印加したシリコン結晶における電子ダイナミクス ===
−
In this exercise, we learn the calculation of electron dynamics in a unit cell of crystalline silicon of a diamond structure.
+
このエクササイスでは、ダイヤモンド構造を持つシリコン結晶の単位セルにおける電子ダイナミクス計算を学びます。
−
Calculation is done in a cubic unit cell that contains eight silicon atoms.
+
計算は、８個のシリコン原子を含む立方体の単位セルで行います。
−
Since the ground state calculation costs much less computational time than the time evolution calculation,
+
時間発展計算に比べて基底状態計算は計算時間がはるかに少ないため、両方の計算を連続して実行します。
−
both calculations are successively executed.
+
基底状態計算を終えたのち、cos^2の包絡形を持つパルス電場を加えます。
−
After finishing the ground state calculation, a pulsed electric field that has cos^2 envelope shape is applied.
+
パルス場を特徴付ける強度、振動数、偏極の向き、キャリアエンベロープ位相などのパラメータは、インプットファイルで指定します。
−
The parameters that characterize the pulsed field such as magnitude, frequency, polarization, and carrier envelope phase
 
−
are specified in the input file.
 
  
−
==== Input files ====
+
==== 入力ファイル ====
−
: To run the code, following files are used:
+
計算を実行するのに、以下のファイルを用います:
  
 
{| class="wikitable"
 
{| class="wikitable"
−
| file name       || description
+
| ファイル名       || 説明
 
|-
 
|-
−
| ''Si_gs_rt_pulse.inp'' || input file that contain namelist variables and their values.         
+
| ''Si_gs_rt_pulse.inp'' || ネームリスト変数とそれらの値を含むインプットファイル       
 
|-
 
|-
−
| ''Si_rps.dat'' || pseodupotential file for Carbon
+
| ''Si_rps.dat'' || シリコン原子の擬ポテンシャルファイル
 
|}
 
|}
  
−
* You may download the above 2 files (zipped file) from:
+
圧縮された上記の２つのファイルを以下からダウンロードできます:
  
−
  [[media:Si_gs_rt_pulse_input.zip| Download zipped input and pseudopotential files]]
+
  [[media:Si_gs_rt_pulse_input.zip|圧縮されたインプットファイルと擬ポテンンシャルファイルのダウンロード]]
  
−
* In the input file ''Si_gs_rt_pulse.inp'', namelists variables are specified. Most of them are mandatory to execute the calculation. We present explanations of the namelist variables that appear in the input file in:
+
インプットファイル''Si_gs_rt_pulse.inp''では、ネームリスト変数が指定されています。
 +
それらの多くは、計算を行うために必須です。
 +
それらの説明を以下に示します:
  
−
  [[Explanation of input files (crystalline silicon under a pulsed electric field)]]
+
  [[インプットファイルの説明（パルス電場のもとでのシリコン結晶）]]
  
−
+
あなたが計算をしたい他の系のためのインプットファイルを準備するのに、この情報は有用です。 インプットファイルで用いることのできるネームリスト変数の完全なリストは、ダウンロードファイルのSALMON/manual/input_variables.mdにあります。
−
: This will help you to prepare the input file for other systems that you want to calculate. A complete list of the namelist variables that can be used in the input file can be found in the downloaded file ''SALMON/manual/input_variables.md''.
 
  
−
==== Output files ====
+
==== 出力ファイル ====
−
: After the calculation, following output files are created in the directory that you run the code,
+
計算の終了後、あなたが計算を実行したディレクトリに以下の出力ファイルが生成されています。
  
 
{| class="wikitable"
 
{| class="wikitable"
−
| file name       || description
+
| ファイル名       || 説明
 
|-
 
|-
−
| ''Si_gs_info.data'' || information of ground state calculation
+
| ''Si_gs_info.data'' || 基底状態計算の情報
 
|-
 
|-
−
| ''Si_eigen.data'' || energy eigenvalues of orbitals
+
| ''Si_eigen.data'' || 軌道のエネルギー固有値
 
|-
 
|-
−
| ''Si_k.data'' || information on k-points
+
| ''Si_k.data'' || k点に関する情報
 
|-
 
|-
−
| ''Si_rt.data'' || electric field, vector potential, and current as functions of time
+
| ''Si_rt.data'' || 時間の関数としての電場、ベクトルポテンシャル、カレント
 
|-
 
|-
−
| ''Si_force.data'' || force acting on atoms
+
| ''Si_force.data'' || 原子に働く力
 
|-
 
|-
−
| ''Si_lr.data'' || Fourier transformations of various quantities
+
| ''Si_lr.data'' || 様々な量のフーリェ変換
 
|-
 
|-
−
| ''Si_gs_rt_pulse.out'' || standard output file
+
| ''Si_gs_rt_pulse.out'' || 標準出力ファイル
 
|}
 
|}
  
−
* You may download the above files (zipped file) from:
+
圧縮された上記のファイルを以下からダウンロードすることができます:
  
−
  [[media:Si_gs_rt_pulse_output.zip|Download zipped output files]]
+
  [[media:Si_gs_rt_pulse_output.zip|圧縮した出力ファイルのダウンロード]]
  
−
* Explanations of the output files are described in:
+
出力ファイルの内容を以下で説明します:
  
−
  [[Explanation of output files (crystalline silicon under a pulsed electric field)]]
+
  [[出力ファイルの説明（パルス電場のもとでのシリコン結晶）]]
  
−
== Maxwell + TDDFT multiscale simulation ==
+
== マクスウェル+TDDFTマルチスケールシミュレーション ==
  
−
=== Exercise-6: Pulsed-light propagation through a silicon thin film ===
+
=== エクササイス-6: パルス光のシリコン薄膜の伝播 ===
−
In this exercise, we learn the calculation of the propagation of a pulsed light through a thin film of crystalline silicon.
+
このエクササイスでは、シリコン結晶の薄膜を伝播するパルス光の計算を学びます。
−
We consider a silicon thin film of ?? nm thickness, and an irradiation of a few-cycle, linearly polarized pulsed light normally on the thin film.
+
数サイクルの直線偏光したパルス光を薄膜に垂直入射します。
−
First, to set up initial orbitals, the ground state calculation is carried out.
+
始めに、軌道関数の初期状態を準備するため、基底状態計算を行います。
−
The pulsed light locates in the vacuum region in front of the thin film.
+
また、パルス光は薄膜の手前の真空領域に置きます。
−
The parameters that characterize the pulsed light such as magnitude and frequency are specified in the input file.
+
パルス光を特徴付ける、強さや振動数などのパラメータは、インプットファイルで指定します。
−
The calculation ends when the reflected and transmitted pulses reach the vacuum region.
+
反射波や透過波が真空領域に到達したら、計算を終了します。
  
−
==== Input files ====
+
==== 入力ファイル ====
−
: To run the code, following files are used:
+
計算を実行するのに、以下のファイルを用います:
  
 
{| class="wikitable"
 
{| class="wikitable"
−
| file name       || description
+
| ファイル名       || 説明
 
|-
 
|-
−
| ''Si_gs_rt_multiscale.inp'' || input file that contain namelist variables and their values.         
+
| ''Si_gs_rt_multiscale.inp'' || ネームリスト変数とそれらの値を含むインプットファイル   
 
|-
 
|-
−
| ''Si_rps.dat'' || pseodupotential file for silicon
+
| ''Si_rps.dat'' || シリコン原子の擬ポテンシャルファイル
 
|}
 
|}
  
−
* You may download the above two files (zipped file) from:
+
圧縮された上記の２つのファイルを以下からダウンロードできます:
  
−
  [[media: Si_gs_rt_multiscale_input.zip| Download zipped input and pseudopotential files]]
+
  [[media: Si_gs_rt_multiscale_input.zip|圧縮したインプットファイルと擬ポテンシャルファイルのダウンロード]]
  
 +
インプットファイル''Si_gs_rt_multiscale.inp''で、ネームリスト変数が指定されています。
 +
それらの多くは、計算を行うために必須です。 それらの説明を以下に示します:
  
−
* In the input file ''Si_gs_rt_multiscale.inp'', namelists variables are specified. Most of them are mandatory to execute the calculation. We present explanations of the namelist variables that appear in the input file in:
+
[[インプットファイルの説明（シリコン薄膜のパルス光伝播）]]
  
−
[[Explanation of input files (pulsed-light propagation through a silicon thin film)]]
+
あなたが計算をしたい他の系のためのインプットファイルを準備するのに、この情報は有用です。 インプットファイルで用いることのできるネームリスト変数の完全なリストは、ダウンロードファイルのSALMON/manual/input_variables.mdにあります。
  
−
+
==== 出力ファイル ====
−
: This will help you to prepare the input file for other systems that you want to calculate. A complete list of the namelist variables that can be used in the input file can be found in the downloaded file ''SALMON/manual/input_variables.md''.
+
計算の終了後、あなたが計算を実行したディレクトリに以下の出力ファイルが生成されています。
−
 
 
−
==== Output files ====
 
−
: After the calculation, following output files are created in the directory that you run the code,
 
  
 
{| class="wikitable"
 
{| class="wikitable"
−
| file name       || description
+
| ファイル名       || 説明
 
|-
 
|-
−
| ''Si_gs_info.data'' || results of the ground state as well as input parameters
+
| ''Si_gs_info.data'' || 入力パラメータと、基底状態の計算結果
 
|-
 
|-
−
| ''Si_eigen.data'' || orbital energies in the ground state calculation
+
| ''Si_eigen.data'' || 基底状態計算での軌道エネルギー
 
|-
 
|-
−
| ''Si_k.data'' || information on k-points
+
| ''Si_k.data'' || k点に関する情報
 
|-
 
|-
−
| ''Si_Ac_xxxxxx.data'' || various quantities at a time as functions of macroscopic position
+
| ''Si_Ac_xxxxxx.data'' || ある時刻における様々な量を、巨視的座標の関数として示す
 
|-
 
|-
−
| ''Si_Ac_M_xxxxxx.data'' || various quantities at a macroscopic point as functions of time
+
| ''Si_Ac_M_xxxxxx.data'' || ある巨視的座標における様々な量を、時間の関数として示す
 
|-
 
|-
−
| ''Si_Ac_vac.data'' || vector potential at vacuum position adjacent to the medium
+
| ''Si_Ac_vac.data'' || 物質に隣接した真空の位置におけるベウトルポテンシャル
 
|-
 
|-
−
| ''Si_gs_rt_multiscale.out'' || standard output file
+
| ''Si_gs_rt_multiscale.out'' || 標準出力ファイル
 
|}
 
|}
  
−
* You may download the above files (zipped file) from:
+
圧縮された上記のファイルを以下からダウンロードすることができます:
  
−
  [[media:Si_gs_rt_multiscale_output.zip|Download zipped output files]]
+
  [[media:Si_gs_rt_multiscale_output.zip|圧縮した出力ファイルのダウンロード]]
  
−
* Explanations of the output files are described in:
+
出力ファイルの内容を以下で説明します:
  
−
  [[Explanation of output files (pulsed-light propagation through a silicon thin film)]]
+
  [[出力ファイルの説明（シリコン薄膜のパルス光伝播）]]

Latest revision as of 18:21, 4 February 2018

はじめに

SALMONエクササイスへようこそ！

このエクササイスでは、SALMONの多岐にわたる応用をカバーする いくつかのサンプルを用いて、基本からSALMONの使い方を説明します。 以下では、UNIX/Linux OSの計算環境にあることを前提にします。 最初に、あなたの計算環境にSALMONをダウンロードしインストールすることが必要です。 まだ行っていない場合は、ダウンロードとインストールと実行に従い作業してください。

インストールと実行で述べたように、SALMONを実行するためには少なくとも インプットファイルと擬ポテンシャルファイルが必要です。 以下では、いくつかのサンプル計算に対するインプットファイルを示し、そのインプットファイルに 現れるネームリスト変数の簡単な説明を行います。 そのインプットファイルを修正して、あなた自身が行いたい計算を実行することができるでしょう。 サンプルに含まれている元素の擬ポテンシャルも提供されています。 主要な出力ファイルの説明も行います。

以下では６つのエクササイスを紹介します。

最初の３つのエクササイス（エクササイス-1~3）は、孤立したアセチレン分子C2H2に対するものです。 もし孤立した系の電子ダイナミクス計算に興味がある場合は、これらのエクササイスをご覧ください。 SALMONでは、通常最初に基底状態計算を行います。 これは、エクササイス-1で説明します。 基底状態計算が終わったのち、２つの電子ダイナミクス計算が用意されています。 エクササイス-2では、分子の分極率と光吸収を得る 実時間での線形光応答計算を示します。 エクササイス-3では、 パルス電場が印加された分子における電子ダイナミクス計算を示します。

次の２つのエクササイス(エクササイス-4~5)は、結晶固体シリコンに対するものです。 もし無限周期系における電子ダイナミクス計算を学ぶことに興味があれば、これらのエクササイスをご覧ください。 周期系の計算では、小さな単位セルを持つ系の基底状態計算は容易であり、時間発展計算の方がはるかに計算時間がかかるため、 基底状態と時間発展の計算を一つのジョブとして実行することを推奨します。 以下の２つのエクササイスも、そのようになっています。 エクササイス-4では、結晶シリコンの誘電関数を得るための 線形応答性質の計算を示します。 エクササイス-5 では、パルス電場を印加した結晶シリコンにおける電子ダイナミクス計算を示します。

最後のエクササイス(エクササイス-6)は、バルクシリコンへのパルス光の入射と伝播に対して、 パルス光の電磁場に対するマクスウェル方程式と単位セルの電子ダイナミクス計算を結合したものです。 これは大きな計算時間を必要とするため、大規模並列スパコンの利用が推奨されます。 エクササイス-6で、 バルクシリコンの表面に垂直に入射する直線偏光パルス光の入射計算を示します。

C2H2 (孤立分子)

エクササイス-1: C2H2分子の基底状態

このエクササイスでは、静的コーン・シャム方程式を解くことによりアセチレン分子(C2H2)の基底状態を求める計算を学びます。 このエクササイスは、分子やナノ粒子のような孤立した系に対して、SALMONを用いた計算をどのようにすれば良いかを知るのに有用です。 今のところ、SALMONで構造最適化を実行することはできないことに注意してください。 分子の原子位置はインプットファイルて指定され、計算の過程では固定されます。

入力ファイル

コードを実行するために、以下のファイルが必要です。

file name description
C2H2_gs.inp ネームリスト変数とそれらの値を含むインプットファイル
C_rps.dat 炭素原子の擬ポテンシャル
H_rps.dat 水素原子の擬ポテンシャル

上記の３つのファイル（圧縮ファイル）は、以下からダウンロードできます: 

圧縮されたインプットファイルと擬ポテンシャルファイルをダウンロード

インプットファイルC2H2_gs.inpでは、ネームリスト変数が指定されています。それらの多くは、基底状態計算を実行するために必須です。それらの説明を以下に示します: 

インプットファイルの説明（C2H2分子の基底状態）

あなたが計算をしたい他の系のためのインプットファイルを準備するのに、この情報は有用です。 インプットファイルで用いることのできるネームリスト変数の完全なリストは、ダウンロードファイルのSALMON/manual/input_variables.mdにあります。

出力ファイル

計算の終了後、あなたが計算を実行したディレクトリに以下の出力ファイルが生成されています。

file name description
C2H2_info.data 基底状態の解の情報
dns.cube 電子密度のキューブファイル
elf.cube 電子局在関数(ELF)
psi1.cube, psi2.cube, ... 電子軌道
dos.data 状態密度
pdos1.data, pdos2.data, ... 射影された状態密度
C2H2_gs.bin 実時間計算で用いるバイナリ出力ファイル

上記のファイル（バイナリファイルC2H2_gs.binを除いた圧縮ファイル）は、以下からダウンロードすることができます: 

圧縮された出力ファイルのダウンロード

軌道エネルギーのような主要な計算結果はC2H2_info.dataに含まれています。C2H2_info.dataと他の出力ファイルの内容を以下で説明します: 

出力ファイルの説明（C2H2分子の基底状態）

画像イメージ

出力ファイルから作成されたいくつかの画像イメージを示します。

image files used to create the image
最高占有軌道(HOMO) psi1.cube, psi2.cube, ...
電子密度 dns.cube
電子局在関数 elf.cube

エクササイス-2: C2H2分子の分極率と光吸収

このエクササイスでは、時間依存コーン・シャム方程式を解いてアセチレン(C2H2)分子の線形応答を調べる計算を学びます。 線形応答計算は、分子の分極率や振動子強度分布を与えます。 このエクササイスは、 エクササイス-1で説明された基底状態計算を終えた後に実行する必要があります。 計算では、分子軸をz軸として、C2H2分子の全ての電子に分子軸方向に撃力的な摂動を加えます。 その後に、外場を加えることなく時間発展計算を実行します。 計算の途中で、電気双極モーメントを保持します。 時間発展計算の後、電気双極モーメントに対して時間-振動数フーリェ変換を行い、振動数依存分極率を得ます。 振動数依存分極率の虚部は、振動子強度分布と光吸収断面積に比例します。

インプットファイル

計算を実行するのに、線形応答計算のためのネームリスト変数とそれらの値を含むインプットファイルC2H2_rt_response.inpが必要です。 基底状態計算で生成したバイナリファイルC2H2_gs.binと擬ポテンシャルファイルも必要です。 擬ポテンシャルファイルは、基底状態計算で用いたものと同じものを用いることが必要です。

ファイル名 説明
C2H2_rt_response.inp ネームリスト変数とそれらの値を含むインプットファイル
C_rps.dat 炭素原子の擬ポテンシャル
H_rps.dat 水素原子の擬ポテンシャル
C2H2_gs.bin 基底状態計算で生成されたバイナリファイル

圧縮されたC2H2_rt_response.inpファイルを以下からダウンロードできます: 

圧縮されたインプットファイルのダウンロード

インプットファイルC2H2_rt_response.inpでは、ネームリスト変数が指定されています。それらの多くは、線形応答計算を実行するために必須です。それらの説明を以下に示します。 

インプットファイルの説明（C2H2分子の分極率と光吸収）

あなたが計算をしたい他の系のためのインプットファイルを準備するのに、この情報は有用です。 インプットファイルで用いることのできるネームリスト変数の完全なリストは、ダウンロードファイルのSALMON/manual/input_variables.mdにあります。

出力ファイル

計算の終了後、あなたが計算を実行したディレクトリに以下の出力ファイルが生成されています。

ファイル名 説明
C2H2_lr.data エネルギーの関数としての分極率と振動子強度分布
C2H2_p.data 時間の関数としての双極モーメントの成分

上記の圧縮ファイルを以下からダウンロードすることができます: 

圧縮した出力ファイルのダウンロード

出力ファイルの内容を以下で説明します: 

出力ファイルの説明（C2H2分子の分極率と光吸収）

エクササイス-3: パルス電場を印加したC2H2分子における電子ダイナミクス

このエクササイスでは、時間依存コーン・シャム方程式を解いてパルス電場を印加したアセチレン(C2H2)分子に起こる電子ダイナミクスの計算を学びます。 計算の出力として、時間の関数としての全エネルギーや電気双極子モーメントなどの量が得られます。 このエクササイスは、エクササイス-1で説明された基底状態計算を終えた後に実行する必要があります。 このエクササイスでは、cos^2の包絡形を持つパルス電場を用いています。パルス電場を特徴付ける強度、振動数、偏極の向き、キャリアエンベロープ位相などのパラメータは、インプットファイルで指定します。

入力ファイル

計算を実行するのに、以下のファイルを用います。C2H2_gs.binは、基底状態計算で生成されたバイナリファイルです。 擬ポテンシャルファイルは、基底状態計算で用いたものと同じものを持ちることが必要です。 従って、パルス電場計算のネームリスト変数とその値を指定するC2H2_rt_pulse.inpが、ユーザが用意する必要のある唯一のファイルです。

ファイル名 説明
C2H2_rt_pulse.inp ネームリスト変数とそれらの値を含むインプットファイル
C_rps.dat 炭素原子に対する擬ポテンシャルファイル
H_rps.dat 水素原子に対する擬ポテンシャルファイル
C2H2_gs.bin 基底状態計算で生成されたバイナリファイル

圧縮されたC2H2_rt_pulse.inpファイルを以下からダウンロードできます: 

圧縮されたインプットファイルのダウンロード

インプットファイルC2H2_rt_pulse.inpでは、ネームリスト変数が指定されています。 それらの多くは、パルス電場により起こる電子ダイナミクス計算を実行するために必須です。 それらの説明を以下に示します: 

インプットファイルの説明（パルス電場のもとでのC2H2分子）

あなたが計算をしたい他の系のためのインプットファイルを準備するのに、この情報は有用です。 インプットファイルで用いることのできるネームリスト変数の完全なリストは、ダウンロードファイルのSALMON/manual/input_variables.mdにあります。

出力ファイル

計算の終了後、あなたが計算を実行したディレクトリに以下の出力ファイルが生成されています。

ファイル名 説明
C2H2_p.data 時間の関数としての電気双極子モーメントの成分
C2H2_ps.data 電気双極子モーメントを時間-振動数フーリェ変換することで得られるパワースペクトル

圧縮された上記のファイルを以下からダウンロードすることができます: 

圧縮された出力ファイルをダウンロード

出力ファイルの内容を以下で説明します: 

出力ファイルの説明（パルス電場のもとでのC2H2分子）

シリコン結晶(周期的な固体)

エクササイス-4: シリコン結晶の誘電関数

このエクササイスでは、ダイヤモンド構造を持つシリコン結晶の線形応答計算を学びます。 計算は、８個のシリコン原子を含む立方体の単位セルで行います。 時間発展計算に比べて基底状態計算は計算時間がはるかに少ないため、両方の計算を連続して実行します。 基底状態計算の終了後、単位セルの全ての電子に対してz方向に撃力の摂動を加えます。 ダイヤモンド構造では誘電関数は等方的なため、計算された誘電関数は摂動の方向には依存しません。 時間発展の間、単位セルの体積で平均を取った電子の流れの密度が計算されます。 電子の流れの密度に対して時間と振動数のフーリェ変換を行うと、振動数依存伝導度が得られます。 誘電関数は伝導度から通常の関係に従い得られます。

入力ファイル

計算を実行するのに、以下のファイルを用います:

ファイル名 説明
Si_gs_rt_response.inp ネームリスト変数とそれらの値を含むインプットファイル
Si_rps.dat シリコン原子の擬ポテンシャル

圧縮された上記の２つのファイルを以下からダウンロードできます: 

圧縮されたインプットファイルと擬ポテンシャルファイルのダウンロード

インプットファイルSi_gs_rt_response.inpでは、ネームリスト変数が指定されています。それらの多くは、計算を行うために必須です。 それらの説明を以下に示します: 

インプットファイルの説明（シリコン結晶の誘電関数）

あなたが計算をしたい他の系のためのインプットファイルを準備するのに、この情報は有用です。 インプットファイルで用いることのできるネームリスト変数の完全なリストは、ダウンロードファイルのSALMON/manual/input_variables.mdにあります。

出力ファイル

計算の終了後、あなたが計算を実行したディレクトリに以下の出力ファイルが生成されています。

ファイル名 説明
Si_gs_info.data 基底状態計算の情報
Si_eigen.data 軌道のエネルギー固有値
Si_k.data k点に関する情報
Si_rt.data 時間の関数としての電場、ベクトルポテンシャル、カレント
Si_force.data 原子に働く力
Si_lr.data 誘電関数のフーリェスペクトル
Si_gs_rt_response.out 標準出力ファイル

圧縮された上記のファイルを以下からダウンロードすることができます: 

圧縮した出力ファイルをダウンロード

出力ファイルの内容を以下で説明します: 

出力ファイルの説明（シリコン結晶の誘電関数）

エクササイス-5: パルス電場を印加したシリコン結晶における電子ダイナミクス

このエクササイスでは、ダイヤモンド構造を持つシリコン結晶の単位セルにおける電子ダイナミクス計算を学びます。 計算は、８個のシリコン原子を含む立方体の単位セルで行います。 時間発展計算に比べて基底状態計算は計算時間がはるかに少ないため、両方の計算を連続して実行します。 基底状態計算を終えたのち、cos^2の包絡形を持つパルス電場を加えます。 パルス場を特徴付ける強度、振動数、偏極の向き、キャリアエンベロープ位相などのパラメータは、インプットファイルで指定します。

入力ファイル

計算を実行するのに、以下のファイルを用います:

ファイル名 説明
Si_gs_rt_pulse.inp ネームリスト変数とそれらの値を含むインプットファイル
Si_rps.dat シリコン原子の擬ポテンシャルファイル

圧縮された上記の２つのファイルを以下からダウンロードできます: 

圧縮されたインプットファイルと擬ポテンンシャルファイルのダウンロード

インプットファイルSi_gs_rt_pulse.inpでは、ネームリスト変数が指定されています。 それらの多くは、計算を行うために必須です。 それらの説明を以下に示します: 

インプットファイルの説明（パルス電場のもとでのシリコン結晶）

あなたが計算をしたい他の系のためのインプットファイルを準備するのに、この情報は有用です。 インプットファイルで用いることのできるネームリスト変数の完全なリストは、ダウンロードファイルのSALMON/manual/input_variables.mdにあります。

出力ファイル

計算の終了後、あなたが計算を実行したディレクトリに以下の出力ファイルが生成されています。

ファイル名 説明
Si_gs_info.data 基底状態計算の情報
Si_eigen.data 軌道のエネルギー固有値
Si_k.data k点に関する情報
Si_rt.data 時間の関数としての電場、ベクトルポテンシャル、カレント
Si_force.data 原子に働く力
Si_lr.data 様々な量のフーリェ変換
Si_gs_rt_pulse.out 標準出力ファイル

圧縮された上記のファイルを以下からダウンロードすることができます: 

圧縮した出力ファイルのダウンロード

出力ファイルの内容を以下で説明します: 

出力ファイルの説明（パルス電場のもとでのシリコン結晶）

マクスウェル+TDDFTマルチスケールシミュレーション

エクササイス-6: パルス光のシリコン薄膜の伝播

このエクササイスでは、シリコン結晶の薄膜を伝播するパルス光の計算を学びます。 数サイクルの直線偏光したパルス光を薄膜に垂直入射します。 始めに、軌道関数の初期状態を準備するため、基底状態計算を行います。 また、パルス光は薄膜の手前の真空領域に置きます。 パルス光を特徴付ける、強さや振動数などのパラメータは、インプットファイルで指定します。 反射波や透過波が真空領域に到達したら、計算を終了します。

入力ファイル

計算を実行するのに、以下のファイルを用います:

ファイル名 説明
Si_gs_rt_multiscale.inp ネームリスト変数とそれらの値を含むインプットファイル
Si_rps.dat シリコン原子の擬ポテンシャルファイル

圧縮された上記の２つのファイルを以下からダウンロードできます: 

圧縮したインプットファイルと擬ポテンシャルファイルのダウンロード

インプットファイルSi_gs_rt_multiscale.inpで、ネームリスト変数が指定されています。 それらの多くは、計算を行うために必須です。 それらの説明を以下に示します: 

インプットファイルの説明（シリコン薄膜のパルス光伝播）

あなたが計算をしたい他の系のためのインプットファイルを準備するのに、この情報は有用です。 インプットファイルで用いることのできるネームリスト変数の完全なリストは、ダウンロードファイルのSALMON/manual/input_variables.mdにあります。

出力ファイル

計算の終了後、あなたが計算を実行したディレクトリに以下の出力ファイルが生成されています。

ファイル名 説明
Si_gs_info.data 入力パラメータと、基底状態の計算結果
Si_eigen.data 基底状態計算での軌道エネルギー
Si_k.data k点に関する情報
Si_Ac_xxxxxx.data ある時刻における様々な量を、巨視的座標の関数として示す
Si_Ac_M_xxxxxx.data ある巨視的座標における様々な量を、時間の関数として示す
Si_Ac_vac.data 物質に隣接した真空の位置におけるベウトルポテンシャル
Si_gs_rt_multiscale.out 標準出力ファイル

圧縮された上記のファイルを以下からダウンロードすることができます: 

圧縮した出力ファイルのダウンロード

出力ファイルの内容を以下で説明します: 

出力ファイルの説明（シリコン薄膜のパルス光伝播）
Retrieved from "https://salmon-tddft.jp/mediawiki/index.php?title=エクササイス&oldid=1586"