Eğer incelemek istediğimiz fiziksel sistemin potansiyeli küresel simetrik bir dağılıma sahipse, Schrödinger Denklemi‘ni küresel koordinatlarda çalışmak oldukça pratik olacaktır. (Bkz. Küresel koordinatlar) Çözümler sırasında bütün diferansiyel denklemlerin detaylarına girmeyeceğim, fakat bunların çözümlerinin getirdiği sonuçlar kesinlikle dikkat çekici olacaktır. Bu yüzden sonucun neden böyle çıktığını daha iyi anlamak adına bu tipteki diferansiyellerin çözümüne olabildiğince değinmeye çalıştım. Fakat elbette tümünü bir başlık altında çözmemiz mümkün değil, bu sebeple bir kısmını sizin ilginize bırakıyorum. Yazının sonunda kuantum sayıları olan l ve m nin küresel koordinatlarda schrödinger denklemi çözümü sırasında nasıl elde edildiğini göreceğiz.
Schrödinger Denklemi;


Klasik mekanikte kinetik enerji ile potansiyel enerjinin toplamını ifade eden Hamiltonyen operatörü olarak gösterdiğimiz H’ı aşağıdaki şekilde açabiliriz.


Burada


Eğer buradaki p ifadesinin karesini alacak olursak,


Şimdi elimizdekileri (1) no’lu Schrödinger denkleminde yerleştirmeye başlayabiliriz. Öncelikle (2) no’lu denklem ile verilen Hamiltonyen operatörünü (1) no’lu Schrödinger denkleminde yerine yazalım.


Böylelikle potansiyel enerji (V) ve dalga fonksiyonu


yazılır. Böylelikle


sağlanır. Zamana bağımlı Schrödinger denkleminin çözümü ise c_n sabitler olmak üzere


şeklinde ifade edilir.
Potansiyel, yalnızca orijinden olan uzaklığın bir fonksiyonu olduğundan küresel koordinatları seçmek işleri kolaylaştıracaktır. Çünkü küresel simetrik bir yapıda merkezden orijine olan koordinatlardaki değişimi kartezyen koordinat sisteminde (x,y,z)’nin tümüyle, küresel koordinatlarda ise (r,Θ,Φ)’den yalnızca r’ye bağlı olarak ele almak mümkündür. Bunu yapmak için öncelikle 7 no’lu denklemde ∇² ifadesi yerine küresel koordinatlardaki Laplace operatörünü yazmalıyız. Küresel koordinatlarda Laplace operatörü


olarak verildiğinden (9) no’lu denklemi (7) no’lu denklemde yerine yazarsak


Denklem bir hayli karışık görünüyor. Fakat elimizde mantıksal olarak kategorilendirilebilecek gibi duran değişkenler var. Bu yüzden değişkenlere ayırma metodunu kullanarak


Şimdi (11) no’lu denklemi, (10) no’lu denklemde yerine yazabiliriz. Burada yaptığımız tek şey \psi fonksiyonunu iki ayrı parçaya ayırmak oldu. (10) no’lu denklemde ilgili türevlerin olduğu yere ilgili parçaların geldiğine dikkat edin. Örneğin ilk parçada dD(r)/dr geliyor, bu aslında d(D(r)Y(Θ,Φ))/dr’dir. Fakat Y(Θ,Φ) parçası r’ye bağlı olmadığından bu kısımdan bir etki gelmez. Uygun şekilde düzenleyip yazarsak


Karmaşık görünen bu denklemi bir iki matematik hareketiyle ihtişamlı görüntüsüne kavuşturabiliriz. Öncelikle


Şimdi de son bir hamle olarak D(r)Y(Θ,Φ) ile bölelim


Böylelikle denklemimiz tamamen değişkenlerine ayrılmış oldu. Denklemin bazı parçaları yalnızca r’ye bağlı iken, geri kalanlar da yalnızca Θve Φ’ye bağlı. Daha anlaşılır şekilde düzenlersek


Bu son durumda üst tarafın yalnızca r’nin, alt tarafın ise yalnızca Θve Φ‘nin bir fonksiyonu olduğuna dikkat edin. İki farklı parametreye ait parçaların toplamı sıfıra eşit olduğuna göre, bu parçalar birbirinin zıt işaretlisi olan iki sabite eşittir. Bu sabite x gibi herhangi bir sayı değeri atamak yerine l(l+1) diyeceğiz. Bunun sebebi çözümü getiren Legendre polinomlarıdır. Elbette biz sonucun ne olduğunu bildiğimiz için böyle bir yol izliyoruz, çözümü yeniden keşfetmediğimizi unutmayın. Direkt doğru yöntemi uygulayarak sonuca ulaşıyoruz. O yüzden kafanız “Nasıl da bunun geleceğini bildiler?” diyerek karışmasın. Bu şekilde yeniden düzenlersek


Böylelikle denklemimizi istediğimiz şekilde düzenlemiş olduk. Bundan sonra yapmamız gereken bu iki ayrı parçayı matematiksel metotlar izleyerek çözmek. Öncelikle üstteki denklem parçacığı ile başlayalım. Bu kısım yalnızca r’ye bağlı olduğundan dikine (radyal) parçadır. Çözümü yapmak için öncelikle f(r)=rD(r) şeklinde bir fonksiyon tanımlamamız gerekiyor. Böylesi bir fonksiyonu, denklemdeki parçalarını elde ederek yerine koyduğumuzda istediğimiz tipte bir denkleme ulaşacağız.
Öncelikle (16) no’lu denklemde üstte yer alan dikine (radyal) ifadeyi, D(r) ile çarpıyoruz.


Denklemin sağındaki sabit terim yerine D(r) ifadesini ekleyerek işleri biraz karmaşıklaştırmışız gibi görünebilir, fakat amacımız denklemin başındaki türev ifadesini yalın hale getirmekti. Şimdi tanımladığımız f(r) fonksiyonu üzerinden türev alıp, bu ifadede yerine yazabiliriz.


Bu durumda denklemde yer alan aşağıdaki ifade f(r) cinsinden


şekline gelir. Şimdi bunu (17) no’lu denklemde yerine yazarak


Burada 1/r ifadesi biraz rahatsız edici duruyor. Bu yüzden denklemi r ile çarpıyoruz.


Şimdi denklemimiz çok daha sade bir görüntüye kavuştu. Son bir hamleyle tek boyuttaki Schrödinger denklemine benzetmemiz mümkün görünüyor. Tek boyutta Schrödinger denklemi aşağıdaki şekilde veriliyordu.


(21) no’lu denklemi


Şimdi yapmamız gereken denklemin sağ tarafındaki ifadeyi sola geçirip, E ifadesini de sağa geçirmek


Böylelikle elde ettiğimiz denklem şekil olarak (22) no’lu tek boyuttaki Schrödinger denklemiyle aynıdır. (22) no’lu denklemle kıyasladığımız zaman büyük parantez içerisinde kalan ifadenin aslında tamamının potansiyel olduğunu görürüz. Burada potansiyel V’ye ek olarak
Sırada (16) no’lu denklemimizde alt kısımda verilen Θve Φ‘ye bağlı açısal kısım var. Bu kısmın sonucunu matematiksel olarak daha etkileyici bulduğumdan sona bırakmak istedim. Kolayca görebilmek adına (16) no’lu denklemi tekrar yazalım.


Artık alttaki denklem ile ilgileniyoruz. Bu denklemi Y(Θ,Φ)\sin²(Θ) ile çarparsak


Şimdi yine, daha önce yaptığımız gibi değişkenlere ayırma metodunu izleyeceğiz. Y(Θ,Φ) fonksiyonu Θ ve Φ‘nin fonksiyonları olarak ayrı ayrı yazılabilir.


Burada kolayca anlaşılsın diye Θ (teta) için T harfini, Φ (fi) için de F harfini seçtim. Karmaşık sembollerle donatmak yerine, böyle okumanın daha kolay olacağını düşünüyorum. Elbette siz istediğiniz bir harfi veya sembolü atayabilirsiniz. (27) no’lu ifadeyi (26)’da yerine yazalım


İçeride yer alan türevleri aşağıdaki şekilde çözebiliriz.


Sağdaki ifade yalnızca Φ’ye bağlı bir fonksiyonun Θ’ya göre değişimi 0 olacağından gider. Böylece denklem parçası


halini alır. Bunu yerine koymadan önce ilgilenenmeniz gereken diğer denklem parçası var.


Eşitliğin solda kalan parçası yalnızca Θ’ya bağlı fonksiyonun Φ’ye göre değişimi 0 olduğundan gider.


Bir kez daha bu ifadenin Φ’ye göre değişimini alırsak


Yine aynı sebepten denklem parçası sadeleşerek aşağıdaki şekilde düzenlenebilir.


(30) ve (34) no’lu denklemler (28) no’lu denklemde yerine yazılırsa


Eşitliği T(Θ)F(Φ) ile bölüp, sağdaki ifadeyi de sola alırsak denklem


şeklini alır. Böylelikle denklemin ilk terimi yalnızca Θ, ikinci terimi ise yalnızca Φ cinsinden yazılarak denklem değişkenlerine ayrıştırılmış oldu. Yine daha önceki sebeplerimizden ötürü bu sefer bir tarafı m²’ye diğer tarafı -m²’ye eşitleyeceğiz.


Alttaki diferansiyel denklemin çözümü oldukça kolaydır.


Bu tipteki bir denklemde


Buradaki türev ifadesini aşağıdaki şekilde açıp yerine yazalım.


Burada


olarak karşımıza çıkar. Bu durumda yalnızca Φ’ye bağımlı fonksiyonumuz


şeklinde yazılabilir. Burada negatif olan çözümün nereye gittiğini sorabilirsiniz. Burayı pozitif alarak m’nin negatif değerleri için bu koşulu da sağladığımıza dikkat edin. Böylelikle m hem negatif hem pozitif olabilir. Aslında buradaki çözüm tam olarak
Küresel koordinatlarda Φ’yi 2π kadar artırdığımızda aynı noktaya geri döndüğümüze göre


yazabiliriz. Bunu alıp (42) no’lu denklemde yerine yazarsak


Ve şimşeklerin çaktığı an! Böylesi bir eşitliğin sağlanabilmesi için
Şimdi Θ denklemine dönebiliriz.
Bu tipteki diferansiyel denklemin çözümü diğeri kadar kolay değildir. Çözüm
olarak verilir. Birleşik Legendre fonksiyonu ise aşağıdaki şekilde tanımlanır.
Burada verilen Pₗ(x) l-inci Legendre polinomudur. n-inci derecen bir Legendre polinomu Rodrigues formülü ile aşağıdaki şekilde tanımlanır.
İlk birkaç Legendre polinomu aşağıdaki şekildedir.


Legendre polinomlarını bildiğimize göre birleşik (assosiye) legendre fonksiyonlarını da bulabiliriz. Burada m=0 için
m=0 ve l=0,1,2,3 değerleri için
m=1 ve l=1,2,3 değerleri için
m=2 ve l=2,3,4 değerleri için
Burada dikkatinizi bir şeyin çekmiş olması gerekiyor. Daima |m|≤l koşulundaki sonuçlardan bahsettik. Bunun sebebi |m|>l durumunda (d/dx)’in mertebesinin, önüne geldiği l-inci Legendre polinomunu sıfır yapmasıdır. Böylesine bir durum için m’nin alabileceği (2l+1) tane değer vardır. Bu değerler
m=-l, -l+1,…,-1,0,1,…,l-1,l aralığındadır. (l=0,1,2,3,…) olmak üzere.


Şimdi normalizasyon işlemini uygulayabiliriz. Hacim elemanı üzerinden integral alacak olursak
Bunları ayrı ayrı da normlayabileceğimiz açıktır.
Normlanmış açısal dalga fonksiyonlarının matematiksel çok özel bir karşılığı vardır ve fizikte sıklıkla karşımıza çıkar. Bunlara küresel harmonikler diyoruz. Aşağıdaki şekilde tanımlanır.
Buradaki ε aşağıdaki aşağıdaki ilk birkaç küresel harmonikte de deneyimleyeceğiniz üzere m≥ 0 için
İlk birkaç küresel harmonik
Bazı ‘lerin grafik gösterimi
İşte burada gördüğümüz l açısal momentum kuantum sayısını (ya da azimutal kuantum sayısını), m de manyetik kuantum sayısını ifade eder. Böylelikle lisede görüp de nereden geldiğine pek de anlam veremediğimiz bu iki kuantum sayısı bu şekilde elde edilir.
Ögetay Kayalı
Kaynaklar
1. Kuantum Mekaniğine Giriş – II Ders Notları – Cemal Parlak – Ege Üniversitesi
2. Kuantum Mekaniğine Giriş – David J. Griffiths – Çeviri: Haluk Özbek, Sondan Durukanoğlu Feyiz
3. http://www.rpi.edu/dept/phys/Courses/phys410/lct4.pdf
4. http://hyperphysics.phy-astr.gsu.edu/hbase/quantum/sch3d.html
5. http://mathworld.wolfram.com/LegendrePolynomial.html
6. http://hyperphysics.phy-astr.gsu.edu/hbase/math/legend.html
7. http://mathworld.wolfram.com/AssociatedLegendrePolynomial.html
8. http://mathworld.wolfram.com/SphericalHarmonic.html