Izingo modelis

Paskelbta gruodžio 9, 2010.

Izingo modelis – tai apibendrintas feromagnetizmo matematinis modelis statistinėje fizikoje. Šiame modelyje magnetinį sukinį turinčios dalelės yra patalpinamos į apibendrinto grafo viršūnes. Bendru atveju grafo struktūra gali būti įvairi, bet įprastu atveju apsiribojama įvairių dimensijų gardelėmis. Modelinės sistemos elgesys yra stebimas esant įvairioms temperatūroms ir yra ieškoma fazinio virsmo taško. Taigi pagrindinė šio modelio problema ir tikslas yra supaprastinto fazinio virsmo modeliuojamoje struktūroje atkūrimas. Norint suprasti ir pažinti Izingo modelio interpretacijų įvairovę vertėtų paskaityti 1, 2 darbus, nes šiame aprašyme mes apsiribosime tik viena galima, skaitmenine, modelio interpretacija – sąlyčio su termostatu interpretacija.

Svarbu pastebėti, kad šis modelis puikiai tinka ne tik įprastiems taikymams statistinėje fizikoje, bet ir turi tarpdisciplininių taikymų. Juk į grafus galima patalpinti realių žmonių elgesį imituojančius agentus, kurių nuomones vienu ar kitu klausimu galima būtų apibendrintai vaizduoti kaip sukinius, ir taip bandyti imituoti žmogišką bandos jausmą (pvz. 3, 4, 5). Vieną pavyzdžiais pateiktų modelių 5 jau buvome pristatę ankstenėje Rizikos fizikos svetainėje (žr. http://mokslasplius.itpa.lt/rizikos-fizika/apibendrintas-sukiniu-modelis-rinkai).

Apie Izingo modelį

Kiekviena modeliuojamos sistemos dalelė gali sąveikauti tiek su kitomis dalelėmis, tiek su išoriniu magnetiniu lauku, taigi šios sistemos Hamiltonianas gali būti užrašytas kaip

$a390289cc85ca69ad8260c388e30cb1f T 000000 0 ordinary Izingo modelis interactive models general models statistinė fizika seni modeliai magnetizmas Ising$

čia $409cd168358cde98a1a14fef5971cc98 T 000000 0 inline Izingo modelis interactive models general models statistinė fizika seni modeliai magnetizmas Ising$ sąveikos tarp dviejų dalelių, i-tosios ir j-tosios, stipris, $90d787b653612059c0066b25313ddae2 T 000000 0 inline Izingo modelis interactive models general models statistinė fizika seni modeliai magnetizmas Ising$ i-tosios dalelės sukinys, $98b337a12c772c46050314c613de0031 T 000000 0 inline Izingo modelis interactive models general models statistinė fizika seni modeliai magnetizmas Ising$ magnetinio lauko stipris.

Šis Hamiltonianas yra ganėtinai bendras, o originaliame Izingo modelyje jis yra šiek tiek supaprastinamas. Visu pirma tariama, kad dalelių sukiniai gali išsidėstyti tik lygiagrečiai pasirinktai ašiai, tad $90d787b653612059c0066b25313ddae2 T 000000 0 inline Izingo modelis interactive models general models statistinė fizika seni modeliai magnetizmas Ising$ projekcija į tą ašį, $e406ac4d7c470823a8619c13dd7101be T 000000 0 inline Izingo modelis interactive models general models statistinė fizika seni modeliai magnetizmas Ising$ , gali įgyti tik $967ffa3ca82c4b8aad1075067fb3fec5 T 000000 0 inline Izingo modelis interactive models general models statistinė fizika seni modeliai magnetizmas Ising$ reikšmes. Padarius šią prielaidą apie sukinius galima atsižvelgti tik į vieną magnetinio lauko $98b337a12c772c46050314c613de0031 T 000000 0 inline Izingo modelis interactive models general models statistinė fizika seni modeliai magnetizmas Ising$ projekciją, projekciją į pasirinktą ašį, $2510c39011c5be704182423e3a695e91 T 000000 0 inline Izingo modelis interactive models general models statistinė fizika seni modeliai magnetizmas Ising$ . Taip pat atkreipkime dėmesį į tai, kad sąveikos stipris tarp dviejų dalelių turėtų greitai mažėti didėjant atstumui tarp dalelių. Taigi paprasčiausiu atveju galima padaryti prielaidą, kad dalelė sąveikauja tik su savo artimiausioms kaimynėmis. Dvimatės gardelės, mūsų pasirinktos struktūros, atveju šių kaimynių yra keturios. Taigi Hamiltonianas yra supaprastinamas iki

$2f065c9d423f9c5ad5fadcf5027a68c7 T 000000 0 ordinary Izingo modelis interactive models general models statistinė fizika seni modeliai magnetizmas Ising$

čia $45b32ea04698cee08d688a80f4a9c2ad T 000000 0 inline Izingo modelis interactive models general models statistinė fizika seni modeliai magnetizmas Ising$ ir $5fdb0734a2a8679264029c65df7a492b T 000000 0 inline Izingo modelis interactive models general models statistinė fizika seni modeliai magnetizmas Ising$ yra dvimatės gardelės koordinatės (sumavimus atliekamas per jas).

Originalaus Izingo modelio atveju sąveikos stipris tarp kaimynių dalelių, $ff44570aca8241914870afbc310cdb85 T 000000 0 inline Izingo modelis interactive models general models statistinė fizika seni modeliai magnetizmas Ising$ , yra laikomas lygiu 1 – modeliuojama medžiaga žemiau kritinės temperatūros turi feromagnetinių savybių, o aukščiau – paramagnetinių. Taigi iš Hamiltoniano išraiškos galima būtų išimti $ff44570aca8241914870afbc310cdb85 T 000000 0 inline Izingo modelis interactive models general models statistinė fizika seni modeliai magnetizmas Ising$ , bet kadangi programoje realizavome ir priešingą, antiferomagnetinį ( $7043bc98118b720f14c9d064a7322d15 T 000000 0 inline Izingo modelis interactive models general models statistinė fizika seni modeliai magnetizmas Ising$ ), atvejį, tai $ff44570aca8241914870afbc310cdb85 T 000000 0 inline Izingo modelis interactive models general models statistinė fizika seni modeliai magnetizmas Ising$ Hamiltoniano išraiškoje palikome. Visgi toliau tekste apsiribosime tik feromagnetinio atvejo analize, nes antiferomagnetiniu atveju fazinis virsmas nėra stebimas.

Temperatūra į modelį yra įvedama naudojant termostatą. Kas kartą atsitiktinai pasirinkus dalelę, ji yra sujungiama su pasirinktos temperatūros termostatu. Sujungtai sistemai leidžiama nusistovėti – t.y. nauja dalelės būsena yra parenkama atsitiktinai pagal energijos stebėjimo tikimybę. Kaip žinome ši tikimybė yra nustatoma iš Bolcmano pasiskirstymo – tikimybė stebėti būseną su energija $3a3ea00cfc35332cedf6e5e9a32e94da T 000000 0 inline Izingo modelis interactive models general models statistinė fizika seni modeliai magnetizmas Ising$ yra proporcinga $f601ba3e73539ef1b1a591c4682e3625 T 000000 0 inline Izingo modelis interactive models general models statistinė fizika seni modeliai magnetizmas Ising$ .

Modelio fazės ir fazinis virsmas modelyje

Iš tokio temperatūros įvedimo seka, kad modelio elgsena bus skirtinga esant skirtingoms temperatūroms – egzistuos dvi skirtingos fazės. Aukštos temperatūros atveju dominuos netvarka – dalelės didins savo energiją išorinės šilumos sąskaita. O jei temperatūra bus žema, tai dalelės negalės pasiimti energijos iš aplinkos, tad dominuos jų polinkis minimizuoti savo energiją – dominuos tvarka. Tvarkos, feromagnetinės, fazės atveju iš pradžių susiformuos ir plėsis vienodų sukinių domenai, kol galų gale vienas domenas apims visą gardelę. Įdomus yra tarpinis atvejis – kai už tvarką ir netvarką atsakingi procesai atsveria vienas kitą. Tokiu atveju, krizinės temperatūros aplinkoje, stebėsime besiformuojančius domenus ir tai kaip jie nyksta dėl atsitiktinių sukinio „apsivertimų“. Žemiau, 1 pav., demonstruojame tipinius gardelių vaizdus šiais trim atvejais.

phases Izingo modelis interactive models general models statistinė fizika seni modeliai magnetizmas Ising

1 pav. Trys skirtingos Izingo modelio fazės: paramagnetinė, krizinė, feromagnetinė.

Norint skaitmeniškai rasti krizinę temperatūrą, reikia perrinkti galimas temperatūras stebint kaip kinta įmagnetėjimas ir sistemos energija. Žemiau, 2 pav., mes demonstruojame kaip kinta šie sistemos stebimieji kintamieji keičiant modelio parametrą $65919ef78270ec07b1d1fec959d7fe41 T 000000 0 inline Izingo modelis interactive models general models statistinė fizika seni modeliai magnetizmas Ising$ ( $a412e23d2ab4696295db4aec6e2851c7 T 000000 0 inline Izingo modelis interactive models general models statistinė fizika seni modeliai magnetizmas Ising$ ), kuris turi atvirkštinės temperatūros išmatuotos sukinių energijos vienetais prasmę. Kaip matome iš paveikslų fazinis virsmas vyksta intervale $963fb2441f258b0d2d1a96e7dd1b4f4f T 000000 0 inline Izingo modelis interactive models general models statistinė fizika seni modeliai magnetizmas Ising$ .

temperature ising metrapolis 2D Izingo modelis interactive models general models statistinė fizika seni modeliai magnetizmas Ising

2 pav. Sistemos absoliutus įmagnetėjimas (a) ir energija (b) stebima po 10⁷ sukinių apvertimų esant skirtingoms modelį valdančioms parametro ΔE/kT vertėms.

Lyg šiol laikėme, kad dalelių neveikia išorinis magnetinis laukas. Kas pasikeistų jei įjungtume išorinį magnetinį lauką? Išorinis magnetinis laukas iš esmės nieko nepakeičia, o tik pažeidžia sukinių simetriją – sukiniai linksta pasisukti magnetinio lauko kryptimi. Žemiau, 3 pav., pateikiame h-1/T fazinę diagramą.

magnetic field temperature ising phases Izingo modelis interactive models general models statistinė fizika seni modeliai magnetizmas Ising

3 pav. h-1/T fazinė diagrama: raudonoje ir mėlynoje srityse įsivyrauja vienos krypties sukiniai (s=1 raudonoje, s=-1 mėlynoje), violetinė sritis atitinka paramagnetinę fazę, dryžuota sritis apie h=0 atitinka feromagnetinę fazę, o žalia sritis krizinę būseną. h=0 sritis, dėl vizualumo, yra padidinta likusių fazinės erdvės sričių atžvilgiu.

Java programėlė

Šioje programėlėje realizuotas aukščiau aptartas Izingo modelis 50×50 dvimatėje sukinių gardelėje. Programėlė išveda keturis paveikslus – sukinių paveikslas, energijos, entropijos ir įmagnetėjimo kitimo grafikai per paskutines 100 iteracijų. Viena iteracija atitinka tūkstantį bandymų apversti dalelės sukinį. Pabaigus kiekvieną iteraciją atnaujinami visi keturi paveikslai. Skaičiavimai pradedami nuspaudus mygtuką „Pradėti“ ir yra tęsiami kol vartotojas nenuspaudžia mygtuko „Stabdyti“. Sustabdžius skaičiavimą galima jį vėliau pratęsti arba pradėti skaičiavimą iš naujo.

Aukščiau šio užrašo Jūs turėtumėte matyti Java naršyklės programėlę. Jei jos nematote, patikrinkite ar turite įsidiegę JRE ir ar Jūsų naršyklėje Java yra įjungta. Taip pat įsitikinkite ar turite naujausią prieinamą JRE versiją. Naujausią JRE versiją galima parsisiųsti iš http://java.com/getjava.

Literatūra

R. J. Baxter. Exactly solved models in statistical mechanics. Academic Press Inc., Harcourt Brace Jovanovich Publishers, London, 1982. http://tpsrv.anu.edu.au/Members/baxter/book.
J. P. Sethna. Statistical Mechanics: Entropy, Order Parameters and Complexity. Clarendon Press, Oxford, 2009. http://pages.physics.cornell.edu/sethna/StatMech/.
S. Bornholdt. Expectation bubbles in a spin model of markets: Intermittency from frustration across scales. International Journal of Modern Physics C 12 (5), 2001, psl. 667-674.
T. Kaizoji, S. Bornholdt, Y. Fujiwara. Dynamics of price and trading volume in a spin model of stock markets with heterogeneous agents. Physica A 316, 2002, psl. 441-452.
S. H. Yook, H. J. Kim, Y. Kim. Agent-based generalized spin model for financial markets on two-dimensional lattices. Journal of the Korean Physical Society 52, 2008, psl. S150-S153.

Pagrindinės temos Bendri modeliai, Interaktyvūs modeliai. Raktažodžiai Ising, magnetizmas, seni modeliai, statistinė fizika.

Palikti atsiliepimą (RSS)

Rizikos fizika

Izingo modelis

Apie Izingo modelį

Modelio fazės ir fazinis virsmas modelyje

Java programėlė

Literatūra

Palikti atsiliepimą Atšaukti atsiliepimą

Pagrindinės temos

Paskutiniai 5 darbai DB

Rizikos fizika

Izingo modelis

Apie Izingo modelį

Modelio fazės ir fazinis virsmas modelyje

Java programėlė

Literatūra

Palikti atsiliepimą Atšaukti atsiliepimą

Pasidalinti

Pagrindinės temos

Raktažodžiai

Paskutiniai 5 darbai DB