コード入力
ノンシーケンシャルモードの体積にミー散乱の分布を適用するには、体積のObject PropertiesダイアログボックスのBulk ScatterタブにおいてMie.DLLファイルを選択します。
添付ファイル:
Fig1_MieInputsCropped.jpg [ 33.2 KiB | 表示回数: 645 回 ]
このDLLは下記五つの入力が必要です。
１．Particle Index
これは（実際の）散乱粒子の屈折率です。バックグランド媒質の屈折率、すなわち粒子が置かれる媒質の屈折率は、体積のMaterialコラムで指定されたものにより定義されます。粒子の屈折率がバックグラウンド媒質の屈折率と一致する場合（屈折率の差が10^-6以下の場合）、粒子がバックグラウンドと区別できないため散乱は起きません。粒子の屈折率が、<= 0.0に設定される場合、散乱は起きません（DLLは無視されます）。
２．Size (microns)
これは球形粒子の半径（R）です。この値を<= 0.0に設定した場合、散乱は起きません（DLLは無視されます）。許容可能な最大値は、正規化された球体サイズ(x = 2πr/λ)の許容可能な最大値に基づいており、現在は1934と設定されています。正規化された球体サイズが1934を超える場合、散乱は起きません。そのような大きい球体では、ZEMAXで球体自身をモデリングするほうがより適切です。もっと正確に言えば、球体に起因する実効散乱分布をモデリングすることが適切です。
３．Density (cm^-3)
これは体積の散乱粒子の密度（ns）です。この値を< 0.0に設定した場合、散乱は起きません（DLLは無視されます）。この値を> 0.0に設定した場合、体積における散乱の平均自由行程λmfp = 1/nsσsを求めるために、コードにより計算された散乱の横断面（σs）とともにこの値が利用されます。この密度を＝0.0に設定した場合、体積における散乱の平均自由行程がMean Pathパラメータ入力において指定あれた値と等しくなります（密度が> 0.0の場合、Mean Pathパラメータ入力において指定された値は無視されます）。
４．Min. thresh
ZEMAXにおける散乱計算は確率分布関数の積分を伴います（詳細は“How to Create a User-Defined Scattering Function”http://www.zemax.com/kb/articles/219/1/How-to-Create-a-User-Defined-Scattering-Function/Page1.html（英語）に記述されています）。積分は極性方向において0～π（方位角方向において0～2π）から通常は得られます。しかし、ピーク分布では、確立がゼロである角度空間の領域に対する積分を実行するには計算的に非効率です。最小閾値はユーザーにより極積分のより小さい上限を定義することを可能にします。上限は、確率分布がピーク分布のほんのわずかの割合と等しい（割合が最小閾値と等しい）際の極角空間における位置により算出されます。例えば、最小閾値が＝0.01と設定された場合、極積分の上限は、ピーク確率の0.01倍に等しい確率分布の位置になります。入力値が<= 0.0 または >= 0.5と設定された場合、デフォルト値の0.001が使われます。
この値を指定する際には注意が必要です。一般的に大きい粒子の確率分布はよりピーク形状であるため、最小閾値をより大きな値に設定しがちですが、散乱計算における角度分解能は固定（0.1度）であるため確率分布がよりピーク形状になれば1つの角度から次の角度までの確率が大幅に降下してしまいます。従って、積分において用いる十分なデータポイントを確保するためには、一般的に最小閾値を小さくする必要があります。
このパラメータの適切な値を判断する最良の方法は、このパラメータが異なった場合の結果を比較することです。0.01の値から始めて、10の倍数で減少させます。結果が大きく異なる場合は、結果が収束するまで値を減少させていきます。それ以上は、値を減少させてても計算時間を増やすだけで計算精度を上げることにはなりません。
５. Transmission
これは散乱現象を受ける全ての光線の密度を減らす減衰係数です。この係数はユーザーにより、散乱球体におけるエネルギー吸収を考慮することを可能にします。バックグラウンド媒質のエネルギー吸収は、体積のMaterialコラムにおいて定義された媒質のTransmission（透過）プロパティにより決定されます（詳細は“Fitting Index Data in ZEMAX”http://www.zemax.com/kb/articles/62/4/Fitting-Index-Data-in-ZEMAX/Page4.html（英語）に記述されています）。
散乱粒子に対する屈折率の架空コンポーネントが吸収に影響を与えますが（DLLにおけるTransmission入力により特徴付けられます）、このコンポーネントは散乱の確率分布にも影響を与えます。しかし、架空コンポーネントのマグニチュードが~ 10^-2またはそれ以上の（極めて吸収率の高い粒子に相当する）場合にだけ確率分布に影響を与えます。よって、散乱分布における架空コンポーネントの影響はここで紹介されるミーモデルでは無視されます。
最後に、AngleパラメータはこのDLLでは使用されませんがAngleパラメータの値が<= 0.0 または > 180.0として設定された場合、散乱はおきません（DLLは無視されます）。
ミー散乱分布は、散乱粒子の屈折率とバックグランド媒質の屈折率との比率(m = nparticle/nbackground)および正規化された球体サイズ(x = 2πr/λ)の二つの無次元パラメータにより特徴づけられます。ミー分布の個々の計算は、各（m, x）ペアに対し実施されなければなりませんが、これらの計算が完了した際には結果がメモリに保存され、任意の光線数を追跡するのに利用できます。最大10,000 ペアの(m,x) 値のデータをDLLで保存することができます。この限界を超えた場合には、それ以上のミーデータを保存することができず、それ以降の（m,x）ペアでは散乱が起きません。ほとんどのアプリケーションでは、10,000ペアで十分です。
