Mie saçılması

**Mie saçılması** veya **Mie teorisi**, düzlem bir elektromanyetik dalganın (ışık) homojen bir küre tarafından saçılmasını ifade eder. Maxwell denklemlerinin **Lorenz–Mie–Debye çözümü** olarak da bilinmektedir. Denklemlerin çözümü sonsuz bir vektör küresel harmonik serisi şeklinde yazılır. Saçılma ismini fizikçi Gustav Mie'den almaktadır; analitik çözümü ilk kez 1908 yılında yayınlanmıştır.

"Mie saçılması" terimi aynı zamanda ışığın dalga boyunun kürenin yarıçapı ile yaklaşık eşit olduğu saçılma durumları için de kullanılmaktadır; bu saçılma özellikle atmosferde görülür. Genel Mie saçılmasının teknik olarak bir boyut limiti yoktur: saçılma sonuçları kürenin dalga boyundan küçük olduğunu durumlarda Rayleigh saçılmasına, büyük olduğu özel durumlarda ise geometrik optiğe yakınsar. Mie saçılması ve tesir kesitinin hesaplanması için birçok kod bulunmaktadır; Mie teorisinin konsentrik küreler, sonsuz silindir ve küre kümelerine uygulandığı kodlar da mevcuttur.

## Teori

Mie saçılmasında çarpan ve saçılan elektrik alanlar vektör küresel harmonikleri olarak yazılır. Helmholtz denkleminin küresel koordinat sisteminde değişkenlerin ayrımı ile çözülmesi ile elde edilir. Kürenin yüzeyinde elektromanyetik sınır koşulları sağlanarak saçılan dalga farklı katsayılarda çok kutuplu harmonikler cinsinden gösterilir. Işık saçılma oranı genellikle optik tesir kesiti ile ifade edilir ve bu katsayı Mie saçılması ile elde edilebilir. Teori kullanarak Mie rezonansları ve saçılma katsayıları da bulunabilmektedir.

Işığın dalga boyunun kürenin yarıçapına oranına ve malzeme özelliklerine göre farklı yaklaşım ve tanımlamalarda bulunulabilir. Rayleigh saçılması dalga boyundan çok daha küçük parçacıklardaki ışık saçılmasını belirtir. Bu dalga boyları için yarı-statik yaklaşım kullanılarak Laplace denklemi de çözülebilir. Rayleigh–Gans yaklaşımı ve ayrıksı kırınım teorisi (van de Hulst yaklaşımı) parçacığın kırılma indisinin bulunduğu ortamınkine çok yakın olduğu durumlarda kullanılabilir. Yaklaşımlar her ne kadar birçok uygulamada iş görse de atmosferdeki su damlacıklarındaki, hücrelerdeki ve emülsiyonlardaki saçılmaların hesaplanması için tam teorinin kullanılması esastır.

## Matematiksel temeller
### Mie saçılmasının türetilmesi

Mie saçılmasının standart analitik çözümü genelde *z*-ekseninde hareket eden ve *x*-ekseninde polarize olmuş bir düzlem dalga için yapılır. Parçacığın yalıtkanlık sabiti ve manyetik geçirgenliği \varepsilon_1 ve \mu_1 ile gösterilirken, parçacığın bulunduğu ortamdaki değerler \varepsilon ve \mu ile gösterilir.

Saçılma probleminin çözümü için öncelikle vektör Helmholtz denklemi küresel koordinatlarda yazılır. Helmholtz denklemi, elektrik ve manyetik alanlar için şu şekilde gösterilir:

:
\nabla^{2} \mathbf{E}+{k}^{2} \mathbf{E}=0, \ \ \ \ \nabla^{2} \mathbf{H}+{k}^{2} \mathbf{H}=0


Helmholtz denklemi dışında, alanların  \nabla \cdot \mathbf{E}=\nabla \cdot \mathbf{H}=0,  \nabla \times \mathbf{E}=i \omega\mu \mathbf{H} ve  \nabla \times \mathbf{H}=-i \omega\varepsilon \mathbf{E} koşullarını sağlaması gerekir. Daha sonra \mathbf{M} ve \mathbf{N} küresel harmoniklerinin *üreten fonksiyonu* olarak kabul edebileceğimiz skalar \psi fonksiyonu denkleme eklenir; vektör küresel harmonikler gerekli koşulları sağlamaktadır.

Küresel koordinatlarda açılan dalga denklemi şu şekilde ifade edilebilir:

:\frac{1}{r^2}\frac{\partial}{\partial r} \left(r^2 \frac{\partial \psi}{\partial r} \right) +
\frac{1}{r^2 sin \theta}\frac{\partial}{\partial \theta} \left(sin \theta \frac{\partial \psi}{\partial \theta} \right) +
\frac{1}{r^2 sin \theta} \left(\frac{\partial^2 \psi}{\partial \phi^2} \right) + k^2 \psi = 0


Bu denklem değişkenlerin ayrımı yöntemi ile çözülür. Denklemin çözümüne göre üreten fonksiyonunu sağlayan vektör küresel harmonikler şu şekilde ifade edilir:

: 
\mathbf{M}_{^e_o m n}=\nabla \times\left(\mathbf{r} \psi_{^e_o m n}\right)
 &nbsp;— manyetik harmonikler (TE)
: 
\mathbf{N}_{^e_o m n}=\frac{\nabla \times \mathbf{M}_{^e_o m n}}{\mathbf{k}}
 &nbsp;— elektrik harmonikler (TM)

ve
: 
{\psi_{e m n}=\cos m \varphi P_{n}^{m}(\cos \vartheta) z_{n}({k} r)}  
: 
{\psi_{o m n}=\sin m \varphi P_{n}^{m}(\cos \vartheta) z_{n}({k} r)}


P_{n}^{m}(\cos \theta)&nbsp;Asosiye Legendre polinomlarını ve  z_{n}({k} r)  ise küresel Bessel fonksiyonlarını gösterir.

Daha sonrasında küreye çarpan düzlem dalga vektör küresel harmonikler cinsinden açılır:
: 
\mathbf{E}_{inc}=E_0e^{ikr\cos\theta}\mathbf{e}_x=E_{0}\sum_{n=1}^{\infty} i^n\frac{2n+1}{n(n+1)}\left( \mathbf{M}^{(1)}_{o1n}(k, \mathbf{r})-i \mathbf{N}_{e1n}^{(1)}(k, \mathbf{r})\right)

: 
\mathbf{H}_{inc}=\frac{-k}{\omega\mu}E_{0}\sum_{n=1}^{\infty} i^n\frac{2n+1}{n(n+1)}\left( \mathbf{M}^{(1)}_{e1n}(k, \mathbf{r})+i \mathbf{N}_{o1n}^{(1)}(k, \mathbf{r})\right)

Buradaki (1) üst-imi \psi_{^e_omn} fonksiyonun radyal kısmının küresel Bessel fonksiyonu olduğunu belirtir. Açılım katsayıları ise şu integraller ile elde edilir:

:
\frac{\int_0^{2\pi}\int_0^{\pi}\mathbf{E}_{inc}\cdot\mathbf{M}^{(1)}_{^e_omn}\sin\theta d\theta d\varphi }{\int_0^{2\pi}\int_0^{\pi}|\mathbf{M}^{(1)}_{^e_omn}|^2\sin\theta d\theta d\varphi }

m\neq 1'ı sağlayan tüm katsayılar 0'dır.

Sonrasında şu koşullar çözümde göz önünde bulundurulur:

# Kürenin yüzeyinde elektromanyetik sınır koşulları sağlanmalıdır.

# Saçılma probleminin çözümü orijinde sınırlandırılmalıdır: bu nedenle kürenin içindeki alanların üreten fonksiyonu olarak küresel Bessel fonksiyonları seçilir.

# Saçılan alan, orijinden sonsuza doğru hareket etmelidir; bunun için birinci tip Hankel fonksiyonları seçilir.

Saçılan alanlar, daha sonra vektör harmonik açılımla şu şekilde yazılır:
: 
\mathbf{E}_{s}=\sum_{n=1}^{\infty} E_{n}\left(i a_{n} \mathbf{N}_{e 1n}^{(3)}(k, \mathbf{r})-b_{n} \mathbf{M}_{o 1 n}^{(3)}(k, \mathbf{r})\right)

: 
\mathbf{H}_{s}=\frac{k}{\omega\mu}\sum_{n=1}^{\infty} E_{n}\left(a_{n} \mathbf{M}_{e 1n}^{(3)}(k, \mathbf{r})+ib_{n} \mathbf{N}_{o 1 n}^{(3)}(k, \mathbf{r})\right)


Bu denklemlerde (3) üst-imi \psi_{^e_omn} fonksiyonun radyal kısmının küresel Hankel fonksiyonu olduğunu belirtir.

Kürenin içindeki alanlar ise şu şekilde açılabilir:

:
\mathbf{E}_{1}=\sum_{n=1}^{\infty} E_{n}\left(-i d_{n} \mathbf{N}_{e 1n}^{(1)}(k_1, \mathbf{r})+c_{n} \mathbf{M}_{o 1 n}^{(1)}(k_1, \mathbf{r})\right)

: 
\mathbf{H}_{1}=\frac{-k_1}{\omega\mu_1}\sum_{n=1}^{\infty} E_{n}\left(d_{n} \mathbf{M}_{e 1n}^{(1)}(k_1, \mathbf{r})+ic_{n} \mathbf{N}_{o 1 n}^{(1)}(k_1, \mathbf{r})\right)


k= \frac{\omega}{c}n ve k_1= \frac {\omega }{c}{n_1}&nbsp;küre dışındaki ve içindeki dalga vektörünü ifade eder. Kürenin dışı ve içindeki kırılma indisleri n ve n_1 ile gösterilir.

Sınır koşulları denklemlere uygulandıktan sonra Mie katsayıları elde edilir:

: 
c_n(\omega) = \frac {\mu_1\left[ \rho h_n(\rho)\right]'j_n(\rho) - \mu_1\left[ \rho j_n(\rho)\right]'h_n(\rho)}{\mu_1\left[ \rho h_n(\rho)\right]'j_n(\rho_1) - \mu\left[ \rho_1 j_n(\rho_1)\right]'h_n(\rho)} 

: 
d_n(\omega) = \frac {\mu_1n_1n\left[ \rho h_n(\rho)\right]'j_n(\rho) - \mu_1n_1n\left[ \rho j_n(\rho)\right]'h_n(\rho)}{\mu n_1^2\left[ \rho h_n(\rho)\right]'j_n(\rho_1) - \mu_1 n^2\left[ \rho_1 j_n(\rho_1)\right]'h_n(\rho)},

: 
b_n(\omega) = \frac {\mu_1\left[ \rho j_n(\rho)\right]'j_n(\rho_1) - \mu\left[ \rho_1 j_n(\rho_1)\right]'j_n(\rho)}{\mu_1\left[ \rho h_n(\rho)\right]'j_n(\rho_1) - \mu\left[ \rho_1 j_n(\rho_1)\right]'h_n(\rho)} 

: 
a_n(\omega) = \frac {\mu n_1^2\left[ \rho j_n(\rho)\right]'j_n(\rho_1) - \mu_1 n^2\left[ \rho_1 j_n(\rho_1)\right]'j_n(\rho)}{\mu n_1^2\left[ \rho h_n(\rho)\right]'j_n(\rho_1) - \mu_1 n^2\left[ \rho_1 j_n(\rho_1)\right]'h_n(\rho)},


:\rho=ka
:\rho_1=k_1a

Bu katsayılarda a kürenin yarıçapına, j_n ile h_n&nbsp;fonksiyonları da sırasıyla birince tip küresel Bessel ve Hankel fonksiyonlarına tekabül eder.

Mie saçılması nümerik olarak hesaplanırken sonsuz seri toplamının bir yerden sonra kesilmesi gerekir. Bununla ilgili en yaygın kriterlerden biri Wiscombe'un kriteridir ve L_{max}=x+4x^{1/3}+2 terimlerin doğruya yakın sonuçlar için yeterli olduğunu belirtir. Bu denklem de *x*, \rho ile eştir.

### Saçılma ve yok olma tesir kesitleri

Mie teorisinde ve saçılmada etkinlik katsayıları sıklıkla kullanılan parametrelerdendir. Bu katsayılar yok olma (Q_e), saçılma (Q_s) ve soğurma (Q_a) için tanımlanabilir. Bu katsayılar ilgili tesir kesitlerinin ışığın saçıldığı kürenin kesit alanına oranıdır. Örnek olarak, bir Mie saçılmasında saçılma etkinlik katsayısı  Q_s = \frac{\sigma_s}{\pi a^2}  formülü ile hesaplanabilir; bu denklemde \sigma_s saçılma tesir kesiti ve *a* da parçacık yarıçapıdır.

Yok olma kesiti, Mie teorisinde şu şekilde tanımlanır:
: \sigma_e = \sigma_s + \sigma_a and  Q_e = Q_s + Q_a 

Gene Mie teorisine göre, saçılma ve yok olma katsayıları sonsuz bir seri ile gösterilebilir:
: Q_s = \frac{2\pi}{k^2}\sum_{n=1}^\infty (2n+1)(|a_{n}|^2+|b_{n}|^2) 
: Q_e = \frac{2\pi}{k^2}\sum_{n=1}^\infty (2n+1)\Re(a_{n}+b_{n}) 

## Uygulamalar

Mie teorisi, atmosfer biliminde bulut ve tozlardaki ışık saçılmasının incelenmesi için sıkça kullanılmaktadır. Teorinin diğer kullanım alanları arasında biyomedikal sistemler ile radarlar bulunmaktadır. Mie saçılımı aynı zamanda bazı metamateryal ve plazmonik yapıların teorisinde yer edinmiştir.

## Kaynakça

## Kitap kaynakları
*
*
*
*

## Dış bağlantılar
* SCATTERLIB: Işık saçılması kodları koleksiyonu 
* Online Mie saçılması hesaplayıcısı 
* PyMieScatt , Mie saçılması Python kodu