UNIVERZITET U BEOGRADU MATEMATICˇKI FAKULTET Dusˇan Zˇ. Onic´ TERMALNO ZRACˇENJE OSTATAKA SUPERNOVIH U RADIO-PODRUCˇJU doktorska disertacija Beograd, 2013 UNIVERSITY OF BELGRADE FACULTY OF MATHEMATICS Dusˇan Zˇ. Onic´ THERMAL RADIATION OF SUPERNOVA REMNANTS IN RADIO DOMAIN Doctoral Dissertation Belgrade, 2013 Mentor: prof. dr Dejan Urosˇevic´, vanredni profesor, Univerzitet u Beogradu, Matematicˇki fakultet Clmnovi komisive: prof. dr Denis Lihi, Univerzitet u Kalgariju, Katedra za fiziku i astronomiju doc. dr Bojan Arbutina, docent, Univerzitet u Beogradu, Matematicˇki fakultet dr Branislav Vukotic´, naucˇni saradnik Astronomske opservatorije u Beogradu Zˇeleo bih da izrazim veliku zahvalnost mentoru, prof. dr Dejanu Urosˇevic´u, jer me je uveo u veoma interesantnu oblast istrazˇivanja ostataka supernovih, kao i za iskazano veliko strpljenje i poverenje tokom izrade ove doktorske disertacije. Tako -de, zˇeleo bih da istaknem veliku pomoc´ i podrsˇku, kao i uvek dobronamerne kritike doc. dr Bojana Arbutine i dr Branislava Vukotic´a. Istakao bih i veoma bitne savete prof. dr Denisa Lihija, bez kojih ova teza ne bi bila potpuna. Posebnu zahvalnost dugujem prof. dr Vladimiru Cˇadezˇu i prof. dr Trajku Angelovu na korisnim savetima i podrsˇci tokom izrade ove doktorske disertacije. Ovim putem bih se zahvalio i na podrsˇci svih kolega sa Katedre za astronomiju. Konacˇno, ova doktorska disertacija svakako ne bi nastala bez podrsˇke i razumevanja porodice i najblizˇih prijatelja. U Beogradu, maja 2013. D.O. TERMALNO ZRACˇENJE OSTATAKA SUPERNOVIH U RADIO-PODRUCˇJU Sazˇetak Evolucija ostataka supernovih vezana je za kretanje jednog, bezsudarnog udarnog talasa, nastalog prilikom prostiranja odbacˇenog zvezdanog materijala nakon eksplo- zije supernove, velikom brzinom kroz me -duzvezdanu sredinu. Teorijska razmatranja udarnih talasa su neophodna radi boljeg razumevanja emisionih maglina i, gene- ralno, procesa u me -duzvezdanoj sredini. U okviru ove disertacije, izme -du ostalog, predstavljeni su rezultati magnetohidrodinamicˇke (MHD) teorije udarnih talasa sa posebnim osvrtom na moguc´a resˇenja u svetlu modela idealne radijativne MHD za opticˇki gust slucˇaj. Specijalno, diskutovana su moguc´a svojstva udarnih talasa kada je dozvoljen skok adijabatskog indeksa i/ili parametra koji odre -duje odnos gasnog prema ukupnom pritisku. Osnovna hipoteza, razmatrana u okviru ove disertacije, odnosi se na moguc´nost znacˇajnog ucˇesˇc´a termalnog zakocˇnog zracˇenja u neprekidnom radio-spektru poje- dinih Galakticˇkih ostataka supernovih. Pozitivan odgovor na ovo pitanje bi mogao predstavljati prirodno objasˇnjenje blago konveksnih neprekidnih radio-spektara po- jedinih Galakticˇkih ostataka koji se prostiru kroz sredinu koja je nesˇto gusˇc´a od prosecˇne. U tom smislu je vazˇno utvrditi postojanje indikatora termalnog ansambla elektrona, na dovoljno niskim temperaturama i dovoljno velikim gustinama, kako bi se zakocˇno zracˇenje moglo detektovati u radio-podrucˇju (blizina, interakcija ili pro- stiranje ostatka kroz molekulski oblak, prisustvo ohla -denih elektrona koji su ranije zracˇili u X-podrucˇju, u fazama nakon Sedov-Tejlorove, detekcija niskofrekventnog krivljenja radio-spektra usled termalne apsorpcije vezane za sam ostatak, detekcija u H , identifikacija radio-rekombinacionih linija vezanih za sam ostatak, itd). Znacˇajno prisustvo termalne komponente u principu mozˇe objasniti i radio- spektralne indekse manje od 0.5 odre -dene kod nekolicine evolutivno starijih ostataka (uglavnom mesˇane morfologije), koji se sˇire kroz gusˇc´u sredinu. Zapravo, u tom slucˇaju bi manji spektralni indeksi, pod pretpostavkom prostog stepenog zakona, predstavljali prirodnu manifestaciju znacˇajnog udela termalnog zakocˇnog zracˇenja u neprekidnom radio-spektru ostataka supernovih. Ipak, trenutno dostupni neprekidni radio-spektri Galakticˇkih ostataka nisu do- voljno precizno odre -deni da bi se mogao dati konacˇan zakljucˇak u svetlu termalne hipoteze. Detaljnija rasprava je jedino moguc´a u slucˇaju tri Galakticˇka ostatka (3C396, IC443, 3C391) za koje je udeo termalne komponente, bez obzira na znacˇajnu neodre -denost, procenjen u okviru ove disertacije. Nova posmatranja c´e, u bliskoj buduc´nosti, stvoriti uslove za donosˇenje cˇvrsˇc´ih zakljucˇaka u smislu postojanja tzv. radio-termalno aktivnih ostataka supernovih. Konacˇno, u okviru ove disertacije je naglasˇen i znacˇaj posmatranja ostataka su- pernovih u X i -podrucˇju, kao i multifrekvencione analize. Pored toga, predlozˇena je moguc´a identifikacija -zracˇenja ostatka supernove 3C434.1 na osnovu preliminarne analize posmatranja sa svemirskog teleskopa Fermi. Kljucˇne recˇi: ostaci supernovih, me -duzvezdana materija, udarni talasi, magneto- hidrodinamika Naucˇna oblast: Astronomija Uzˇa naucˇna oblast: Astrofizika UDK broj: 524.35:52-6(043.3) THERMAL RADIATION OF SUPERNOVA REMNANTS IN RADIO DOMAIN Abstract The evolution of supernova remnants is linked to the propagation of a collisionless shock wave, formed during the initial expansion of high-velocity supernova ejecta through the interstellar environment. Theoretical studies of shock waves are very important for the analysis of supernova remnants, as well as processes in the in- terstellar medium in general. In this doctoral dissertation, some theoretical results based on the magnetohydrodynamical theory of shock waves are presented, with special emphasis on ideal radiative magnetohydrodynamics for the optically thick case. Particularly, solutions for the case when jump in adiabatic index and/or ratio of gas to total pressure is allowed, are discussed. The main hypothesis of this dissertation is that thermal bremsstrahlung radia- tion at radio continuum frequencies can provide a significant contribution in the case of several Galactic supernova remnants. This hypothesis can give a natural expla- nation for nearly concave up radio continuum spectra of several Galactic supernova remnants that are expanding in the environment with higher than average density. In this context, it is important to identify the existence of the possible indicators of ensemble of thermal electrons at sufficiently low temperatures and sufficiently high densities so that the thermal bremsstrahlung radiation linked to a particular remnant could be observed at radio continuum frequencies (vicinity, interaction or expansion through the molecular cloud, presence of the cooled thermal X-ray electrons during the post Sedov-Taylor phases, detection of low-frequency turnovers associated with thermal absorption linked to the remnant, detection in H , identification of radio recombination lines linked to the remnant, etc). The significant presence of thermal component could theoretically explain radio- spectral indices less than 0.5 measured for several evolutionary older supernova rem- nants, (mainly of mixed-morphology class) that expand in the high density region. Actually, these smaller radio-spectral indices, under the assumption of simple power law, would represent a natural manifestation of a significant fraction of thermal emission at radio continuum frequencies. However, present knowledge of the radio continuum spectra of Galactic super- nova remnants is still not determined precisely enough for any definite conclusions to be made about the inherent thermal radio-emission from supernova remnants. A thorough analysis is only possible in the case of three Galactic supernova remnants (3C396, IC443, 3C391) for which the thermal contribution is determined despite high associated uncertainties. New observations in the near future will lay the groundwork for making firmer conclusions about the existence of the so-called radio thermally active supernova remnants. This dissertation highlights the importance of observations of supernova rem- nants in X and -rays, and multiwavelength analysis is general. Besides, it suggests a possible detection of -rays from supernova remnant 3C434.1 based on the obser- vations made by Fzrmi. Keywords: supernova remnants, interstellar medium, shock waves, magneto- hydrodynamics Scientific area: Astronomy Scientific field: Astrophysics UDC number: 524.35:52-6(043.3) Sadrzˇaj 1 Uvod 1 1.1 Osnovne karakteristike me -duzvezdane sredine . . . . . . . . . . . . . 2 1.2 Od supernovih ka ostacima supernovih . . . . . . . . . . . . . . . . . 4 2 Magnetohidrodinamicˇka teorija udarnih talasa 7 2.1 Osnovni pojmovi fizike plazme . . . . . . . . . . . . . . . . . . . . . . 7 2.2 Model idealne magnetohidrodinamike . . . . . . . . . . . . . . . . . . 14 2.3 Prostiranje poremec´aja kroz fluid . . . . . . . . . . . . . . . . . . . . 16 2.4 Udarni talasi u idealnoj MHD . . . . . . . . . . . . . . . . . . . . . . 17 2.4.1 Granicˇne oblasti . . . . . . . . . . . . . . . . . . . . . . . . . . 19 2.4.2 Udarni talasi . . . . . . . . . . . . . . . . . . . . . . . . . . . 19 2.5 Udarni talasi u slabo-neidealnoj plazmi visoke temperature . . . . . . 22 2.5.1 Jednacˇine za skokove fizicˇkih parametara . . . . . . . . . . . . 24 2.5.2 Slucˇaj kada i/ili w trpe skok . . . . . . . . . . . . . . . . . . 30 2.6 Bezsudarni udarni talasi . . . . . . . . . . . . . . . . . . . . . . . . . 34 2.6.1 Ubrzavanje cˇestica na udarnim talasima . . . . . . . . . . . . 35 3 Dinamicˇka evolucija ostataka supernovih 37 3.1 Osnovne faze evolucije OSN . . . . . . . . . . . . . . . . . . . . . . . 37 3.2 Pulsarske magline . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45 4 Zracˇenje ostataka supernovih 46 4.1 Posmatracˇka klasifikacija OSN . . . . . . . . . . . . . . . . . . . . . . 48 4.1.1 OSN mesˇane morfologije . . . . . . . . . . . . . . . . . . . . . 51 4.2 Radio-zracˇenje ostataka supernovih . . . . . . . . . . . . . . . . . . . 53 4.2.1 Predlozˇeni modeli radio-zracˇenja OSN sa Q 0:5 . . . . . . . 61 4.2.2 Radio-rekombinacione linije . . . . . . . . . . . . . . . . . . . 69 4.3 Infracrveno zracˇenje OSN . . . . . . . . . . . . . . . . . . . . . . . . 70 4.4 OSN u vidljivoj svetlosti . . . . . . . . . . . . . . . . . . . . . . . . . 70 4.5 Ultraljubicˇasto zracˇenje OSN . . . . . . . . . . . . . . . . . . . . . . 72 4.6 Zracˇenje OSN u X-podrucˇju . . . . . . . . . . . . . . . . . . . . . . . 72 4.7 Gama zracˇenje ostataka supernovih . . . . . . . . . . . . . . . . . . . 74 4.8 Molekuli . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 78 5 O postojanju ,,radio-termalno aktivnih”ostataka supernovih 80 5.1 Teorijska razmatranja . . . . . . . . . . . . . . . . . . . . . . . . . . . 80 5.1.1 Termalni ansambl . . . . . . . . . . . . . . . . . . . . . . . . . 84 5.1.2 Niskofrekventno krivljenje radio-spektra . . . . . . . . . . . . 88 5.1.3 Termalna radio-luminoznost . . . . . . . . . . . . . . . . . . . 91 5.2 Termalni model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 94 5.3 Analiza i rezultati . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97 5.3.1 OSN 3C396 (G39.2-0.3) . . . . . . . . . . . . . . . . . . . . . 98 5.3.2 IC443 (G189.1+3.0) i 3C391 (G31.9+0.0) . . . . . . . . . . . . 104 5.3.3 Preostali Galakticˇki OSN . . . . . . . . . . . . . . . . . . . . . 110 5.4 Diskusija . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 114 6 Zakljucˇak i planovi za buduc´i rad 122 Literatura 125 Prilog: Tenzori u R3 139 Biografija autora 146 1 Uvod Dugi niz godina su posmatranja u vidljivom delu elektromagnetnog spektra pruzˇala jedine informacije o vasionskim objektima. Tek relativno skoro, razvojem prvo radio-astronomije, a zatim nakon tzv. kosmicˇke revolucije1, omoguc´eno je da se elektromagnetno zracˇenje nebeskih tela detaljnije proucˇi. Danas, multifrekvenciona analiza predstavlja standardnu proceduru prilikom izucˇavanja razlicˇitih kosmicˇkih fenomena2. Veliki broj fizicˇki razlicˇitih vasionskih objekata/procesa je otkriven do danas. Odgonetanjem znacˇajnog broja zagonetki otvoreno je mnosˇtvo novih pitanja. Izu- cˇavanje kosmosa predstavlja modernu oblast nauke. O tome svedocˇi i nezanemarljiv broj dodeljenih Nobelovih nagrada za fiziku upravo vezanih za doprinos na polju astronomije3. Tako -de, u zˇelji da se objasne razlicˇiti vasionski fenomeni, u poslednje vreme se razvio veliki broj novih disciplina koje na odgovarajuc´i nacˇin kombinuju znanja iz razlicˇitih oblasti nauke (astrohemija, astrobiologija, arheoastronomija, itd). Odlicˇan pregled istorije astronomije dat je npr. u knjizi Longaira (Longair 2006). U okviru ove teze se, pre svega, razmatra zracˇenje ostataka supernovih (dalje OSN, eng. supzrnovv rzmnvnt B hcg), i to konkretno u radio-podrucˇju. U uvod- nom poglavlju ove disertacije bic´e, ukratko, recˇi o najvazˇnijim karakteristikama me -duzvezdane sredine kroz koju se prostiru ovi znacˇajni vasionski objekti kao i o njihovim roditeljskim supernovama. )Pocˇetak kosmicˇke revolucije se obicˇno vezuje za 04. 10. 1957. godine, kada je lansiran prvi vesˇtacˇki satelit. 2Naravno, u zavisnosti od konkretnog problema, katkad dopunjena npr. metodama astro- seizmologije i neutrinske astronomije. 3kixtor Frvnxis Hzss, 1936., bvrtin gylz i Vntony Hzwish, 1974., Vrno Vllvn eznzivs i gowzrt looyrow lilson, 1978., huwrvhmvnyvn Chvnyrvszkhvr i lillivm Vlfrzy Fowlzr, 1983., gusszll VC Hulsz i Joszph Hooton ivylorA Jr, 1993., gvymony DvvisA Jr, bvsvtoshi Koshiwv i gixxvryo Givxxoni, 2002., John CC bvthzr i Gzorgz FC hmoot, 2006., hvul ezrlmuttzr, Wrivn eC hxhmiyt i Vyvm GC gizss, 2011. godine, kao i mnogi drugi naucˇnici koji su svojim otkric´ima doprineli razvoju astronomije. 1 1.1 Osnovne karakteristike me -duzvezdane sredine Dugo godina se smatralo da je prostor izme -du zvezda u Galaksiji prazan4. Danas je poznato da ga ispunjava tzv. me -duzvezdana materija (dalje MZM) koja cˇini oko 10% Galakticˇke mase (Vukic´evic´-Karabin & Atanackovic´-Vukmanovic´ 2004). Za prosecˇnog astronoma MZM predstavlja veliku smetnju koja onemoguc´ava da se po- smatraju slabi vasionski izvori elektromagnetnog zracˇenja. Intenzivno izucˇavanje MZM otpocˇelo je tek u drugoj polovini dvadesetog veka, kada su se i posmatracˇke i teorijske tehnike razvile do tog stepena da je mogao da se uvidi susˇtinski znacˇaj MZM za sve grane astronomije. Izucˇavanje fizike MZM danas predstavlja modernu i veoma aktivnu naucˇnu disciplinu. Me -duzvezdanu materiju cˇine razlicˇite komponente koje se uzajamno prozˇimaju. Hemijski sastav gasne komponente je slicˇan sastavu zvezdanih omotacˇa: oko 67% po masi je vodonik, oko 31% helijum i oko 2% tezˇi elementi. Cˇesticˇna komponenta (prasˇina) cˇini oko 1% ukupne mase MZM. Iako je uglavnom izmesˇana sa gasnom komponentom, prasˇina mozˇe formirati i manje ili vec´e lokalne koncentracije koje znacˇajno slabe svetlost zvezda. Hladna i gusta gasna komponenta se najcˇesˇc´e javlja u formi atomskih ili molekulskih oblaka (superoblaci, gigantski molekulski oblaci, oblaci atomskog vodonika, Bokove globule, itd). Iako zauzima svega oko 1% za- premine ona ima znacˇajan udeo u ukupnoj masi MZM. Standardni tipovi i karak- teristicˇne vrednosti nekih fizicˇkih parametara gasne komponente MZM predstavljeni su u tabeli 1. Specificˇne strukture MZM, magline, predstavljaju oblasti lokalnog zgusˇnjenja MZM i uglavnom su nehomogene. Uobicˇajena je posmatracˇka klasifikacija maglina na tamne i svetle (za detalje videti npr. Vukic´evic´-Karabin & Atanackovic´-Vukmano- vic´ 2004). OSN se tako svrstavaju u svetle emisione magline. MZM je prozˇeta magnetnim i gravitacionim poljem, uz prisustvo kosmicˇkih zraka5 (visokoenergetskih cˇestica) i elektromagnetnog zracˇenja. Interesantno je za- 4Kompletnosti radi, korisno je spomenuti i hipotezu o ,,etru”(supstanciji koja bi ispunjavala sav vasionski prostor), koja je napusˇtena nakon cˇuvenog negativnog rezultata Majkelson-Morlijevog eksperimenta iz 1887. godine (Vlwzrt Vwrvhvm bixhzlson, 1852-1931 – dobitnik Nobelove nagrade za fiziku 1907. godine; Zywvry lillivms borlzy, 1838-1923). 5Za detetalje videti npr. knjigu Sˇlikajzera (Schlickeiser 2002). 2 paziti i da su gustine energije elektromagnetnog zracˇenja, magnetnog polja, kosmi- cˇkih zraka i toplotnog kretanja MZM istog reda velicˇine, oko 1 eVRcm3 u okolini Sunca. Magnetno polje u Galaksiji je u principu homogeno na skali duzˇina reda velicˇine kpc. Srednja vrednost Galakticˇkog magnetnog polja (magnetne indukcije) iznosi oko 5 G. Neke od pomenutih komponenti MZM su u ravnotezˇi pritisaka (detaljnije o ra- zmatranju toplotnih faza MZM videti u dole navedenoj literaturi). Dok se npr. molekulski oblaci odrzˇavaju sopstvenom gravitacijom i spoljasˇnjim pritiskom, ma- gline poput H ii regiona i OSN se sˇire kroz me -duzvezdanu sredinu (dalje MZS). Sve do sada navedeno odnosi se na onu MZM koja je dostupna merenjima. Trenutno preovladava miˇsljenje da gravitacione anomalije na velikim skalama mogu da se objasne ,,tamnom materijom”. Kakva su svojstva i koliki je udeo te ne- luminozne materije josˇ uvek nije dovoljno dobro poznato (Vukic´evic´-Karabin & Atanackovic´-Vukmanovic´ 2004). Danas je literatura posvec´ena izucˇavanju MZM sveobuhvatna i detaljna6. Tabela 1: Standardni tipovi i karakteristicˇne vrednosti nekih fizicˇkih parametara gasne komponente MZM prema knjizi Dopite i Saterlenda (Dopita & Sutherland 2003). Tip i [K] n [cm−3] Molekulska MZM 20 S 103 Hladna neutralna MZM ∼ 100 20− 60 Topla neutralna MZM† ∼ 6000 ∼ 0:3 Topla jonizovana MZM‡ (6− 12)× 103 S 1 Vrela jonizovana MZMr S 106 Q 10−2 † Recˇ je o vrednostima koje vazˇe za oblasti dovoljno udaljene od mladih zvezda. U slucˇaju fotodisocijativnih regiona koji okruzˇuju H ii regione koncentracija mozˇe biti znacˇajno vec´a. ‡ Iako je ova sredina uobicˇajena za H ii regione, ona se vezuje i za pojedine regione daleko od oblasti gde nastaju zvezde (gde koncentracija mozˇe biti oko 0:3 cm−3 pri tem- peraturi reda 6000 K), tzv. slojevi Rejnoldsa. r Ova sredina je zagrejana udarnim talasima OSN ili jakim zvezdanim vetrovima. 6Videti npr. knjige Kaplana (Kaplan 1966), Dopite i Saterlenda (Dopita & Sutherland 2003), Lakea (Lequeux 2005), Kvoka (Kwok 2007) kao i pregledni cˇlanak Koksa (Cox 2005). 3 1.2 Od supernovih ka ostacima supernovih Posmatrajuc´i golim okom, supernova7 (dalje SN) predstavlja pojavu nove zvzzyz na nebu. Posmatrani fenomen, u roku od nekoliko dana, dostizˇe svoj maksimalni sjaj da bi ga zatim, polako, i po godinu dana smanjivao. Nasuprot posmatracˇkoj klasi- fikaciji SN, zasnovanoj na analizi opticˇkih spektara, supernove su fizicˇki prouzroko- vane ili kolapsom jezgra masivnih zvezda (mase zvezde roditelja vec´e od 8M ) ili termonuklearnom eksplozijom C/O belog patuljka u tesno dvojnim sistemima (za detalje videti Arbutina 2005 i tamo navedenu literaturu). Danas se smatra da su- pernove tipa Ia (koje u svom opticˇkom spektru nemaju linije vodonika) poticˇu od starije zvezdane populacije II (termonuklearne SN), dok supernove tipa II i Ib/c vode poreklo od mladih i masivnih zvezda (kolapsirajuc´e SN). Dok se pod kompaktnim ostacima podrazumevaju produkti kolapsa zvezdanog jezgra (npr. neutronske zvezde), ostatke supernovih oznacˇavaju udarni talasi koji raznose zvezdane omotacˇe kroz MZS. Samo je mali broj OSN doveden u direktnu vezu sa (posmatranim) supernovama od kojih poticˇu (npr. SN koju je posmatrao Tiho Brahe8 1574. godine – OSN i–xho i SN koju je posmatrao Johan Kepler9 1604. godine – OSN Kepler). Posmatrani Galakticˇki OSN uglavnom su previˇse stari da bi postojao zapis o njihovim supernovama. Naime, dok fenomen supernove traje relativno kratko (oko godinu dana), ostatak supernove se mozˇe posmatrati reda stotina hiljada godina. Iako su pojedini, stariji OSN, vidljivi kao svetle emisione magline na noc´nom nebu (delovi OSN X–gnus aoop cˇine maglinu Veo u sazvezˇ -du Labud, OSN IC443 je poznat kao maglina Meduza u sazvezˇ -du Blizanci, itd), tek je razvoj radio-astronomije pruzˇio moguc´nost njihove ozbiljnije analize. To nije cˇudno ukoliko se ima u vidu da je prividno najjacˇi radio-izvor na nebu na niskim radio-frekvencijama, Cas A, upravo ostatak supernove. Svakako, tek je razvoj astronomije na svim talasnim duzˇinama pruzˇio moguc´nost stvaranja potpune predstave o fizicˇkim procesima vezanim za ove vasionske objekte. 7Izraz supzrnovv, prvi put su primenili Bade i Cviki (lilhzlm Hzinrixh lvltzr Wvvyz, 1893- 1960; Fritz owixky, 1898-1974), 1934. godine (Baade & Zwicky 1934). 0iyxho Wrvhz, 1546-1601. 1Johvnnzs Kzplzr, 1571-1630. 4 Supernove predaju ogromnu kolicˇinu energije, reda 1051 erg (1 erg = 10−7 J), okolnoj MZS i zajedno sa svojim ostacima cˇine glavne regulatore gasnog sadrzˇaja galaksija. Naravno, ovde je recˇ o kineticˇkoj energiji eksplozije, pri cˇemu neutrini, koji veoma slabo interaguju sa okolnom materijom, mogu odnositi i oko 1053 erg. Tako -de, danas se smatra da je u slucˇaju tzv. hipernovih oslobo -dena energija vec´a za bar red velicˇine od prethodno pomenutih. Inacˇe, svega oko 1% pomenute energije se izracˇi u vidu elektromagnetnog zracˇenja odgovornog za vidljivi fenomen supernove. Materijal odbacˇen eksplozijom supernove (eng. zjzxtv) se efikasno ubacuje u MZS i tako je owogv(xujz tezˇim elementima (od kojih vec´ina upravo mozˇe nastati jedino u supernovama). Naravno i MZM uticˇe na evoluciju OSN, o cˇemu c´e biti viˇse recˇi u narednim poglavljima ove teze. Na slici 1 je predstavljen OSN S147 (maglina Sˇpageti) u vidljivom delu spektra. Recˇ je o veoma starom OSN uglovnog precˇnika od oko 3 stepena a stvarnog precˇnika od oko 42 pc (Green 2009). 5 Slika 1: Maglina Sˇpageti (OSN S147) koja se protezˇe kroz delove sazvezˇ -da Bika i Kocˇijasˇa (izvor: Dvviyz Dz bvrtin ; thz EhVDEhdDcVhV ehotoshop FIih aiwzB rvtor). 6 2 Magnetohidrodinamicˇka teorija udarnih talasa Udarni talasi predstavljaju jednu od najznacˇajnih pojava u me -duzvezdanoj sre- dini. Njihov uticaj na okolinu je viˇsestruk. Oni mogu dovesti do zagrevanja, te disocijacije i jonizacije, kao i kompresije gasa odnosno plazme kroz koju se pro- stiru. Analiza udarnih talasa obuhvata vazˇan deo velikog broja naucˇnih istrazˇivanja u astronomiji: od objasˇnjavanja mehanizma formiranja zvezda pa sve do pitanja vezanih za ubrzavanje cˇestica u MZS. Kako su ostaci supernovih vezani za udarni talas koji se prostire kroz MZS, na samom pocˇetku ove disertacije, bic´e koncizno predstavljena magnetohidrodinamicˇka (dalje MHD) teorija udarnih talasa. Zapravo, recˇ je o najjednostavnijem pristupu proucˇavanja udarnih talasa koji se prostiru kroz plazmu u prisustvu spoljasˇnjeg magnetnog polja (Landau & Lifˇsic 1984; Shu 1992; Vukic´evic´-Karabin 1994; Kivel- son & Russell 1995; Cravens 1997; Goedbloed et al. 2010). Na samom pocˇetku ovog poglavlja bic´e definisani najvazˇniji parametri plazme i predstavljene njihove vrednosti u slucˇaju razlicˇitih vasionskih plazmi (videti poglavlje 2.1). Zatim sledi sazˇeto izlaganje modela idealne magnetohidrodinamike, ostavljajuc´i detalje navede- noj literaturi (videti poglavlje 2.2). Konacˇno, bic´e ukratko predstavljene osobine i tipovi udarnih talasa u svetlu modela idealne MHD (videti poglavlje 2.3) pri cˇemu c´e se ukratko razmotriti i rezultati idealne radijativne MHD u slucˇaju opticˇki guste sredine (videti poglavlje 2.4). Na samom kraju, bic´e predstavljeni osnovni pravci istrazˇivanja udarnih talasa u smislu primene metoda koje izlaze iz okvira idealne magnetohidrodinamike, konkretno, metoda kineticˇke teorije (videti poglavlje 2.5). 2.1 Osnovni pojmovi fizike plazme Sistem elektricˇno10 naelektrisanih cˇestica predstavlja plazmu11 ako je njegovo fi- zicˇko ponasˇanje odre -deno kolektivnom interakcijom zasnovanom na Kulonovim sila- )(U okviru ove disertacije ne postoji moguc´nost zabune vezane za terminologiju, kako se ne razmatraju fenomeni u kojima do izrazˇaja dolazi tzv. boja naelektrisanje. U daljem tekstu bic´e izostavljena recˇ zlzktrixno. ))Langmuir (Irving avngmuir, 1881-1957), dobitnik Nobelove nagrade za hemiju 1932. godine, prvi je iskoristio izraz plvzmv pri opisu jonizovanog gasa 1927. godine (Goossens 2003). 7 ma12 (Milic´ 1977). Posredstvom elektromagnetnog polja, koje poticˇe od svih cˇestica plazme zajedno, svaka pojedinacˇna cˇestica istovremeno interaguje sa svim ostalim cˇesticama. Dakle, recˇ je o kolektivnoj interakciji prouzrokovanoj Kulonovim13 si- lama14. Kolektivna interakcija u plazmi, usled cˇinjenice da je Kulonova sila dugog dometa, osz(xv se vec´ pri relativno niskim stepenima jonizacije i manifestuje se kroz razlicˇite fenomene15. Najupadljivija posledica postojanja kolektivne interakcije jeste tendencija plazme ka elektricˇnoj neutralnosti, odnosno, stanju u kome je ukupna zapreminska gustina naelektrisanja (/el) jednaka nuli. Lokalno nastajanje viˇska pozitivnih ili negativnih naelektrisanja, do kojeg mozˇe doc´i usled toplotnog kretanja cˇestica u plazmi, prac´eno je uspostavljanjem elektricˇnog polja velikog intenziteta koje sputava svako kretanje koje bi dovodilo do daljeg povec´anja prostornog naelektirsanja. Neka se npr. razma- tra plazma koja sacˇinjava Suncˇevu koronu i neka se pretpostavi da se ista sastoji iz elektrona i samo jedne vrste pozitivnih jona (npr. protona). Koncentracije cˇestica (np ≈ ne) su reda velicˇine 1014 m−3 (Goossens 2003). Temperatura u Suncˇevoj ko- roni je reda velicˇine 106 K, sˇto odgovara srednjoj energiji termalnog kretanja od oko 86.2 eV po cˇestici (i [K] ≈ 11600 ki [eV], gde je k tzv. Bolcmanova16 konstanta). Ilustrativno je odrediti koliki bi se potencijal uspostavio na povrsˇini proizvoljne sfere radijusa 1 m, u ovoj plazmi, ako bi se desilo da, usled termalnog kretanja, iz pos- matrane oblasti iza -de 1% elektrona. Lako se mozˇe pokazati, na osnovu izraza za potencijal (<) na povrsˇini sfere radijusa g: < = 1 4."0 f g P (2:1) )2U okviru ove disertacije bic´e razmatrana samo tzv. gasna plazma. Inacˇe, pored gasne plazme, postoje i drugi fizicˇki sistemi kod kojih je dominantna kolektivna interakcija bazirana na Kulonovoj sili. )3ChvrlzsBVugustin yz Coulomw, 1736-1806. )4Pojam plazme se generaliˇse i na slucˇajeve kada je dominantna kolektivna interakcija neelektro- magnetne prirode, npr. slucˇaj kvark-gluon plazme. )5Svakom agregatnom stanju odgovara odre -deni stepen ure -denosti (simetrije) odnosno karakte- risticˇna srednja energija jedne cˇestice (0:01 eV za cˇvrsto stanje, 0:1 eV za tecˇno i 1 eV za gasovito). Ako bi se gasnom stanju dovodila energija do trenutka u kojem bi srednja energija po jednoj cˇestici bila reda 10 eV mogao bi se ocˇekivati fazni prelaz u novo agregatno stanje - gasnu plazmu. Kako je za vec´inu atoma, energija jonizacije oko 10 eV vazˇna osobina plazme je parcijalna (ili potpuna) jonizacija sistema. )6auywig Zyuvry Woltzmvnn, 1844-1906. 8 gde je "0 dielektricˇna propustljivost vakuuma, a f ukupno naelektrisanje u sferi f = /el · k = z(np − ne) · 4 3 .g3P (2:2) da je recˇ o velikim vrednostima potencijala koje bi morali savladati elektroni (∼ 6 kV). Simbol z oznacˇava tzv. elementarno naelektrisanje. Dakle, u sferi radijusa 1 m nije moguc´e da razlika izme -du broja jona i elektrona dostigne 1%. Posmatrani sistem je prakticˇno elektroneutralan. Za proizvoljnu sferu radijusa 10 cm, dobilo bi se da bi za 1% razlike izme -du broja jona i elektrona uspostavljeni potencijal bio reda 600 V, sˇto je opet vec´e od srednje energije termalnog kretanja, pa to opet ne bi bilo moguc´e da se realizuje. Jasno se primec´uje tendencija ka makroskopskoj elektroneutralnosti. Minimalni radijus, iznad kojeg se sistem mozˇe tretirati kao makroskopski elektro- neutralan, predstavlja jedan od najvazˇnijih parametara plazme i naziva se (elektron- ski) Debajev17 radijus (Goossens 2003; Bellan 2004; Bittencourt 2004; Goedbloed & Poedts 2004): rD = √ "0ki nz2 : (2:3) Van tzv. Debajeve sfere, kolektiv naelektrisanih cˇestica c´e biti elektroneutralan, tako da se mozˇe definisati prvi kriterijum plazmenog stanja: a rDP (2:4) gde je a karakteristicˇna dimenzija sistema. Za Suncˇevu koronu, Debajev radijus je reda 6:9 × 10−3 m (videti tabelu 2.1). U slucˇajevima znatnijeg odstupanja od ter- modinamicˇke ravnotezˇe (dalje TDR) mozˇe se uvesti tzv. efektivni Debajev radijus (za detalje videti Milic´ 1977). Inacˇe, poznato je da je pojam termodinamicˇke tem- perature vezan za postojanje Maksvelove18 raspodele po brzinama cˇestica sistema. Ukoliko cˇestice sistema odstupaju od Maksvelove raspodele, tada termodinamicˇka temperatura nije definisana. )7Debaj (eztzr Joszph lillivm Dzwyz, 1884-1966) je dobitnik Nobelove nagrade za hemiju 1936. godine. )0Jvmzs Clzrk bvxwzll, 1831-1879. 9 Svojevrstan odziv plazme na lokalno narusˇavanje makroskopske elektroneutralno- sti (na skalama vec´im od Debajevog radijusa) predstavljaju tzv. (elektronske) plaz- mene oscilacije, frekvencije: !p = √ nz2 "0me P (2:5) gde je same oznacˇena masa elektrona. U tom smislu se mozˇe uvesti i drugi kriterijum plazmenog stanja: e  epP (2:6) gde je sa e oznacˇena karakteristicˇna vremenska skala posmatranog procesa, dok je sa ep obelezˇen period plazmenih oscilacija. Naravno, moguc´e je definisati i tzv. jonsku plazmenu frekvenciju koja odgovara jonskim oscilacijama, koje se srazmerno tezˇe pobu -duju (usled razlike u masi elektrona i jona). Tako -de, recˇ je i o znatno nizˇim vrednostima frekvencija. Josˇ jedna manifestacija kolektivne interakcije jeste tzv. elektrostaticˇko ekrani- ranje (zaklanjanje) na rastojanjima reda Debajevog radijusa. Pojava elektrostaticˇkog ekraniranja se ogleda u tome sˇto potencijal elektrostaticˇkog polja u blizini odre -denog jona tipa , kada se on nalazi u plazmi, nije dat Kulonovim izrazom. Usled privlacˇe- nja naelektrisanja suprotnog znaka i odbijanja naelektrisanja istog znaka, oko svakog jona u plazmi se obrazuje neravnomerna raspodela naelektrisanja u kojoj dominiraju naelektrisanja suprotnog znaka - formira se tzv. owlvk nvzlzktrisvnjv. Unutar tog owlvkv je elektroneutralnost narusˇena. Najbitnija reperkusija formiranja ovog owlvkv jeste pojava zaklanjanja (ekraniranja) polja centralnog jona. Naelektrisane cˇestice se unutar razmatranog owlvkv krec´u toplotnim brzinama usled cˇega se njegova struktu- ra i oblik menjaju u vremenu, no ipak se mozˇe posmatrati, u srednjem, u smislu sfernosimetricˇne strukture. Za slucˇaj termodinamicˇki ravnotezˇnih plazmi izveden je tzv. Debaj-Hikelov19 potencijal (Chen 1974; Vukic´evic´-Karabin 1994; Bittencourt 2004):