Monte Carlo Simülasyonu ve Atomistik Spin Modelleri
Günümüzde fiziğin çeşitli dallarında yapılan araştırmalar teorik, deneysel ve
bilgisayara dayalı (computational) olmak üzere üç kısıma ayrılabilir. Monte
Carlo simülasyon tekniği 1953’de geliştirilmiş olup [1] temelde manyetik sistemleri
karakterize eden spin Hamiltonyenlerinin modellenmesi olmak üzere çok geniş bir
uygulama yelpazesine sahiptir. Bunların arasında sürekli spin sistemleri,
akışkanlar ve polimerlerin simülasyonlarını saymak mümkündür. Bu çalışmada manyetik
sistemlerin modellenmesi üzerinde duracağız.
Önemli Konfigürasyonların Örneklenmesi (Importance Sampling)
İstatistiksel mekanikte genellikle üzerinde çalıştığımız sistemde iç enerji ve
manyetizasyon gibi termodinamik ve manyetik özelliklerin ortalama değerlerini
hesaplamaya sıklıkla ihtiyaç duyarız.
(1)

yer aldığı bir örgüde bu ortalama ifadesi

katkısını barındırır. Tipik olarak

hesaplamak neredeyse imkansız bir hal alır. Bu zorluğu aşmanın bir yolu,
başlangıçta rastgele bir
spin şekillenimi seti'' oluşturarak aranan kesin çözüme
en yakın çözümü belli fiziksel kurallar yardımıyla elde etmektir. Ancak bu
şekillenim setinde bazı şekillenimler bölüşüm fonksiyonuna yeterince büyük katkı
verirken bazı şekillenimlerden ise neredeyse hiç katkı gelmez. Ayrıca söz konusu
şekillenim setinde herhangi bir şekillenimin sahip olduğu enerji gereği (yani
Boltzman ağırlığı
), bölüşüm fonksiyonuna katkı verecek olan
şekillenim sayısı da çok azdır. Bu nedenle bölüşüm fonksiyonu hesaplanırken tüm
şekillenimleri hesaba katmak yerine
ağırlığına göre yeterince
büyük katkı veren şekillenimleri hesaba katmak işi kolaylaştırır.
Importancesampling” adını verdiğimiz bu yöntemde ilk aşama bir \textbf{Markov zinciri}
oluşturmaktır. Sistemin


şekillenimden elde ediliyorsa, bu şekilde üretilen

oluşturduğu sete Markov zinciri adı verilir. Markov zincirlerinde ergodiklik
hipotezi gereğince, belli bir şekillenimden diğerine sonlu adımda ulaşılmalıdır.
Markov Zincirleri ve Master Denklemi
Bir sistem için tanımlanmış bir süreçte belli bir alt durumunun, sürecin tahmin
edilebilir eylemlerinden ve rastgele elemanlardan yola çıkılarak
kestirilebildiği, ancak davranışı deterministik olmayan süreçlere stokastik
süreçler denir. ,
,
kesikli zaman aralıklarında
,
,
ile verilen sonlu sayıdaki olası durumlar seti ile
tanımlanmış stokastik bir süreç ele alalım ve sistemin zamanında bulunduğu
durumu ile gösterelim.
koşullu olasılığına göre
(2)
yazabiliriz. Denklem (2)’den sistemin bir önceki zaman adımındaki


tür süreçlere Markov süreçleri diyoruz. Denklem (2) ile verilen koşullu
olasılık, sistemin tüm durumlarından bağımsız ise ve son durumun olasılığı
sadece bir önceki durumun olasılığına bağlı ise,
(3)
yazabiliriz. Bu durumda denklem (3)’e karşılık gelen

durumlarının dizisine Markov zinciri denir. Böylece, denklem (3) ile
verilen koşullu olasılık
(4)
ile verilen,


tanımlanabilir. Denklem (4) de tüm geçiş olasılıkları gibi
(5)
şartlarını sağlamalıdır. Böylece, sistemin


durumunda bulunma olasılığı

(6)
şeklinde tanımlayabiliriz.
için denklem (6)’nın zamana göre
değişimi master denklemi olarak adlandırılır.
(7)
Denklem (7)’deki olasılık ifadesi korunumludur. Ayrıca,



şekilde



nedeniyle denklem (7)’de tanımlanan master denklemi aynı zamanda bir
süreklilik denklemidir. Basitçe, bu denklem






Metropolis Algoritması
Klasik Metropolis yöntemine göre \cite{metropolis} bir zamanındaki şekillenimi
zamanındaki şekillenimden türetmek
için, ilk ve son durumlar arasındaki enerjilerin farkına bağlı olarak yazılan
bir geçiş olasılığı kullanılır. Üretilen durumlar zaman sırası ile dizilmiş,
ancak deterministik olmayan bir dizi oluşturur. Böyle bir sistemin zamana bağlı
davranışı
(8)
ile verilen master denklemi kullanılarak belirlenebilir. Denklem (8)’de






sağındaki iki terim birbirine eşit olur. Bu durumda elde edilen eşitlik
ayrıntılı denge koşulu (detailed balance condition) olarak adlandırılır ve
(9)
ile verilir. Klasik bir sistemde

(10)
şeklindedir. Paydadaki bölüşüm fonksiyonundan ötürü bu olasılığın değerini kesin
olarak belirlemek genellikle olanaksızdır. Bu zorluğu aşmak için bir Markov
zinciri oluşturulur. Başka bir deyişle, her bir yeni durum, doğrudan bir önceki
durumdan türetilir.


belirlemek için denklem (10)’a göre eski ve yeni durumlar arasındaki
enerji farkını bilmek yeterli olur:
(11)
Denklem (9) ile verilen detaylı denge koşulunu sağlayan her geçiş
olasılığı kabul edilebilir. Bu olasılığın ilk tanımı Metropolis tarafından
(12)
şeklinde yapılmıştır. Burada

(spin-flip) için geçen süredir ve genellikle 1 seçilir.
Metropolis algoritmasının uygulama adımları aşağıdaki gibidir [2,3]
- Başlangıç şekillenimi oluşturulur
-
örgü noktası seçilir
örgü noktasındaki değişkenin işaret değiştirmesi durumunda oluşan enerji farkı
hesaplanır
- rastgele bir
sayısı üretilir,
- Eğer
ise yeni şekillenim kabul edilir
- Bir sonraki örgü noktasına geçilir ve (3)’den itibaren prosedür tekrar edilir
- Markov zincirinin her bir adımında termodinamik değişkenlerin ortalaması hesaplanır.
- İlk birkaç iterasyonda sistem henüz dengeye ulaşmamış olabilir. Bu nedenle ortalama hesaplarken bu adımların ürettiği şekillenimler hesaba katılmaz.
- Sistemi dengeye getirmek için gerekli zaman şu şekilde bulunur: Simülasyon, aynı parametreler için birkaç farklı şekillenimden başlatılarak her bir başlangıç şekillenimi için sistemin manyetizasyonunun ortak saturasyon değerine kadar olan zaman adımları ihmal edilir.
- Genellikle
kadar adımda sistemi dengeye getirmek yeterli olur [3]
Denklem (9)’un tek çözümü Metropolis fonksiyonu değildir. Alternatif
olarak Glauber geçiş fonksiyonu [4]
(13)
şeklindedir. Yüksek sıcaklıklarda (paramanyetik fazda) Metropolis algoritması
her bir adımda seçilen spinin çevrilmesine izin verir ve sistem iki durum
arasında salınım yapar. Bu nedenle Metropolis algoritması yeterince yüksek
sıcaklıklarda ergodiklik özelliğini yitirirken Glauber dinamiği ise tüm sıcaklık
spektrumu için ergodik kalır.
Ising modeli: Manyetik ve Termodinamik Özelliklerin Hesaplanması
(14)
Manyetizasyon:
(15)
İç enerji:
(16)
Duygunluk:
(17)
Isı kapasitesi:
(18)
Binder kümülantı:
(19)




Şekil 1 (a): Ferromanyetik düzende farklı sıcaklıklar için şekillenimler.
,
,
ve
.




Şekil 1 (b): Antiferromanyetik düzende farklı sıcaklıklar için
şekillenimler. ,
,
ve
. Kırmızı
kareler , mavi kareler
spin yönelimlerini göstermektedir.
Sayısal Sonuçlar
Denklem (14) Hamiltonyeni ile tanımlanmış olan iki boyutlu Ising modelin
manyetik ve termodinamik özelliklerinin örgü parametresi ve sıcaklık ile
değişimi (15)-(19) eşitlikleri kullanılarak hesaplanabilir.
Kare örgü için sistemin çeşitli sıcaklıklardaki ferromanyetik
ve antiferromanyetik
şekillenimleri Şekil-(??)’de
görülmektedir. Ferromanyetizma ve antiferromanyetizma, Pauli ilkesi
ve Coulomb etkileşmelerine dayanan olgulardır. ve
gibi iki spin
ele alırsak, bu iki spin arasındaki değiş-tokuş etkileşmesi
olacaktır. Bu durumda için değiş-tokuş enerjisi
negatif
değer alır. Bu nedenle spinler birbirleri ile paralel olarak yönelmeyi tercih
ederler ve kritik sıcaklığın altında sistemin net manyetizasyonu sıfırdan farklı
olur. için ise
enerjisi pozitif olacağı için spinler,
minimum enerji koşulunu sağlamak adına birbirleri ile anti-paralel olarak
yönelirler ve durum antiferromanyetizma olarak isimlendirilir. Antiferromanyetik
sistemlerin taban durumda net manyetizasyonu sıfırdır.
Ferromanyetizma ve antiferromanyetizama için Şekil-1’de gösterilen örgü
şekillenimleri sabit sıcaklıktaki Markov zincirinin bir halkası
olarak düşünülebilir. Bu şekillenimlerden üretilen Markov zincirlerinin sabit
sıcaklıktaki evinimleri yani zaman serileri Şekil-2’de
görülmektedir [5].
Düşük sıcaklıklarda manyetizasyon saturasyon değeri civarındadır. Artan
sıcaklıkla birlikte sistemin manyetizasyonu saturasyon değerinden uzaklaşır ve
sistemin sıcaklığı kritik sıcaklık civarına ulaştığında ve
değerleri
arasında büyük dalgalanmalar ortaya çıkar. Şekil-1’den de görüldüğü gibi kritik
sıcaklık civarında sonlu ve değişken ebatlarda irili ufaklı kümelenmeler oluşur.
Sistem sıcaklığı kritik sıcaklığın üstüne çıktığında ise artık net manyetizasyon
sıfır değeri etrafında salınım yapar. Bu nedenle ortalama manyetizasyon değeri
sıfır olur ve sistem paramanyetik faza geçmiş olur.

Şekil 2: Ferromanyetik sistem için farklı sıcaklıklarda Markov zincirleri
Bir Markov zinciri boyunca sabit sıcaklıkta elde edilen termodinamik ve manyetik
özelliklerin ortalamaları ise o sıcaklıktaki ortalama değerleri
verir. Şekil-3 ‘de kare örgü için manyetizasyon, iç enerji, ısı sığası ve manyetik duygunluğun sıcaklıkla
değişiminin farklı örgü parametreleri için davranışı görülmektedir.
değeri büyüdükçe hesaba katılan spin sayısı artar ve sistemin davranışı
termodinamik limite yaklaşır. Benzer şekilde, arttıkça ısı sığası ve manyetik duygunluk kritik sıcaklıkta
ıraksama davranışı sergiler.


Şekil 3: Kare örgü için manyetizasyon, iç enerji, ısı sığası ve manyetik duygunluk eğrilerinin sıcaklıkla değişimi. Sol paneldeki kesikli çizgi kesin çözümü göstermektedir.
Kritik Üsteller
Manyetizasyon ve korelasyon uzunluğu kritik üstelleri:
(20)
(21)
Bu iki denklem birleştirilirse
(22)
elde edilir. Termodinamik limitte


Ancak üzerinde çalıştığımız sistem sonlu (

kritik sıcaklık noktasında korelasyon uzunluğu

Dolayısıyla lineer boyutları


(23)
yazılabilir. Bu iki eşitlik oranlanırsa
(24)
olduğu kolayca görülebilir.


(25)
Denklem (25)’in sol tarafı farklı


çizilirse


kestiği yer kritik sıcaklığa,

karşılık gelir. Kritik üstel hesaplamanın bir diğer yolu denklem (22)’de



Kritik sıcaklığı yüksek hassasiyetle hesaplamanın bir yolu ise denklem
(19) ile verilen Binder kümülantının

farklı

kümülantın sıcaklıkla değişimi görülmektedir. Farklı

eğrilerin kesim noktası sistemin kritik sıcaklığına karşılık gelir. Sistemin
kritik sıcaklığı belirlendikten sonra
(26)
ile verilen eşitliklerden yararlanılarak net manyetizasyon, manyetik duygunluk
ve ısı kapasitesi için



edilebilir. Kritik üstellerin analitik değerleri,

(27)
ve
(28)

Şekil 4: Sol panel: denklem (19)’den hesaplanan Binder kümülantının
farklı değerleri için sıcaklıkla değişimi. Sağ panel: denklem (25)
ile verilen fonksiyonunun sıcaklıkla değişimi. Eğrilerin birbirini
kestiği noktanın bileşeni sistemin kritik sıcaklığını,
bileşeni ise
üstelini verir.
eşitlikleri ile birbirine bağlıdır. Denklem (26)’dan manyetizasyon ve
manyetik duygunluk için elde edilen sonuçlar Şekil-5 ve
Şekil-6’da gösterilmiştir. Diğer yandan, Şekil-7’den de
görüldüğü gibi, ısı sığasının ile değişimi logaritmik olmadığı için
düzlemindeki davranışı doğrusal değil, eğrisel bir yapıya
sahiptir. Bu durum olmasından kaynaklanır [6].
Kritik üsteller için kabul edilerek elde edilen sonuçlar
Tablo-1’de verilmiştir.

Şekil 5: Net manyetizasyonun ‘de örgü parametresi
ile değişimi. Sol
panel doğrusal, sağ panel ise logaritmik ölçekte çizilmiştir.

Şekil 6: Manyetik duygunluğun ‘de örgü parametresi
ile değişimi.
Sol panel doğrusal, sağ panel ise logaritmik ölçekte çizilmiştir.

Şekil 7: Isı kapasitesinin ‘de örgü parametresi
ile değişimi. Sol
panel doğrusal, sağ panel ise logaritmik ölçekte çizilmiştir.

Tablo 1: 2 boyutlu Ising modelde termal ve manyetik özelliklerin kritik
üstelleri. MC-1: denklem (26)’dan MC-2 ise denklem (25)’den elde
edilen sonuçları göstermektedir.
Kaynaklar
[1] N. Metropolis, A. W. Rosenbluth, M. N. Rosenbluth, A. H. Teller, E. Teller, J. Chem. Phys. 21 (1953) 1087.
[2] D. P. Landau, K. Binder, A Guide to Monte Carlo Simulations in Statistical Physics, SE (2005) p.68
[3] M. E. J. Newman, G. T. Barkema, Monte Carlo Methods in Statistical Physics, Clarendon Press (Oxford) (2001)
[4]R. J. Glauber, J. Math. Phys. 4 (1963) 294.
[5] N. J. Giordano, Computational Physics, (1997) p.217
[6]H. Gould, J. Tobochnik, An Introduction to Computer Simulation Methods, SE (1996) p.572.