差分

このページの2つのバージョン間の差分を表示します。

--- tabuchi:xrayutils [2026/02/26 09:27] – [1.4 回折スペクトルの計算] mtab
+++ tabuchi:xrayutils [2026/02/27 17:27] (現在) – [1.2 構造因子の計算] mtab
@@ 行 176: / 行 176: @@
 </code>
 ++
+  - 構造因子のエネルギー依存性の計算と表示 ++ ここをクリックすると表示 |
+<code python example-02-2.py>
+#
+# 構造因子の計算 (エネルギー依存)
+#
+import xrayutilities as xu
+import numpy
+import matplotlib.pyplot as plt
+GaAs = xu.materials.GaAs
+energy= numpy.linspace(500, 20000, 5000) # 500 - 20000 eV
+F = GaAs.StructureFactorForEnergy(GaAs.Q(1, 1, 1), energy)
+plt.figure(); plt.clf()
+plt.plot(energy, F.real, '-k', label='Re(F)')
+plt.plot(energy, F.imag, '-r', label='Imag(F)')
+plt.xlabel("Energy (eV)"); plt.ylabel("F"); plt.legend()
+plt.show()
+</code>
+{{:tabuchi:ex-02-2.png?600|}}
+++
 ==== - 回折パタンの計算と逆格子マップの表示 ====
@@ 行 275: / 行 302: @@
 {{:tabuchi:ex-04-1.png?600|}}\\
 ++
-  - AlGaN/GaN(Sub)の回折スペクトルの計算と表示 ++ ここをクリックすると表示 |
+  - AlGaN/GaN(Sub)の$\omega-2\theta$スペクトルの計算と表示 ++ ここをクリックすると表示 |
 <code python example-04-2.py>
 ##### プログラムはほとんど example-03-1.py と同じ #####
@@ 行 376: / 行 403: @@
 緩和していないと、(GaNの方が格子定数が長くて、AlNの方が短いので、AlNは横方向に引っ張られて)AlNのc軸方向の面間隔が短くなるので　$2d_{\rm c}\sin\theta = \lambda \Rightarrow \sin\theta = \lambda/2d_{\rm c}$ の関係から回折角 $\theta$ は大きくなるはずだが、実際そうなっていることがわかる。
 ++
-  - 仮に生成したAlGaN/GaN(Sub)の回折スペクトルに対するなんちゃってフィッティング ++ ここをクリックすると表示 |
+  - 仮に生成したAlGaN/GaN(Sub)の$\omega-2\theta$スペクトルに対するなんちゃってフィッティング ++ ここをクリックすると表示 |
 \\ この例では、真ん中あたりに出てくる
   # この結果をフィッティング対象と考えて AlN 組成x と、緩和率 relaxation をパラメータとして
@@ 行 555: / 行 582: @@
 ++
   - AlGaN/GaN(Sub)のロッキングカーブスペクトルの計算と表示 ++ ここをクリックすると表示 |
+後半に出てくる、
+  #### ここまでは先にやった omega-2theta のシミュレーションの時と同一 ####
+と書かれた行までは、2つ目の例と同じコードです。
 <code python example-04-4.py>
 import numpy as np
@@ 行 643: / 行 673: @@
 plt.show()
 </code>
-{{:tabuchi:ex-04-4.png?600|}}
+{{:tabuchi:ex-04-4.png?600|}}\\
+青線が基板(GaN)のロッキングカーブで、オレンジ線がAlGaN/GaNのロッキングカーブ。
+シータはそれぞれについてシミュレーションで計算された $\omega-2\theta$ スペクトルのピーク位置に合わせて計算されてます
+(単純に「一番値が大きいところ」なので、AlGaN/GaN も基板のピークを捕まえてるかも)。
+いずれにしても、AlGaN/GaN のロッキングカーブは、基板(GaN)のスペクトルを含んでしまってます。
+これは、$2\theta$を制限しない計算になっているためで、現実の測定で言うと、検出器前スリットを前回にして測定している状況に相当します。
+現実でも、この様な設定で測定を行うとロッキングカーブを見ているのか$\omega-2\theta$を見ているのかよくわからなくなりますが、
+この計算はそれを再現していることになります。
 ++
 ------