arXiv:0901.0111v1 [astro-ph.GA] 31 Dec 2008
Pro eedingsIAUSymposiumNo.259,2009
K.G.Strassmeier,A.G.Kosovi hev&J.E.Be kman,eds.
2009InternationalAstronomi alUnion DOI:00.0000/X000000000000000X
Cosmic-ray driven dynamo
in galactic disks
M. Hanasz
1, K. Otmianowska-Mazur
2, H. Lesch
3, G. Kowal
2;4,
M. Soida
2, D. W´
olta´
nski
1K. Kowalik
1, R.K. Paw laszek
1B. Kulesza- ˙
Zydzik
2 1Centrefor Astronomy,Ni holasCoperni usUniversity,PL-87148Piwni e/Torun,Poland, mhanaszastri.uni.torun.pl
2
Astronomi alObservatory,JagiellonianUniversity,ul.Orla171, 30-244Krakow,
3
Astronomi alObservatory,Muni hUniversity,S heinerstr.1,D-81679, Germany,
4
DepartmentofPhysi sandAstronomy,M MasterUniversity,1280MainSt.W., Hamilton,ONL8S4M1,Canada.
Abstract.
We present newdevelopmentsonthe Cosmi {Ray driven,gala ti dynamo, mod-eledbymeansofdire t,resistiveCR{MHD simulations,performed withZEUSandPIERNIK odes.Thedynamo a tion,leading tothe ampli ationoflarge{s ale gala ti magneti elds on gala ti rotation times ales,appearsas a result of gala ti dierential rotation, buoyan y ofthe osmi ray omponentandresistivedissipationofsmall{s aleturbulentmagneti elds. Ournewresultsin ludedemonstrationoftheglobal{gala ti dynamoa tiondrivenbyCosmi Rayssuppliedinsupernovaremnants.Anessentialout omeofthenewseriesofglobalgala ti dynamomodelsistheequipartitionofthegasturbulentenergywithmagneti eldenergyand osmi rayenergy,insaturatedstatesofthedynamoonlargegala ti s ales.Keywords.
Galaxies:ISM{magneti elds{ISM: osmi rays{magneti elds{MHD1. Introduction
Numerous pro esses have been proposed to explain initial magneti elds in early galaxies. However all these pro esses, like phase transitions in the early universe or Biermann battery in protogala ti obje ts provide only very week magneti elds of the order of B 10
20
G at the beginning of gala ti evolution. On the other hand, Rees(1987)proposedthatinitialmagneti eldsingalaxiesmighthavebeengenerated in rst stars and then s attered in the interstellar medium (ISM) by SN explosions, and subsequentlyamplied in plerioni (Crab{type) supernova remnants(SNRs). The meanmagneti eldonthegala ti s alehasbeenestimatedtobeoftheorderof10
9 G. Radioobservations(seeBe kBe k2009,thisvolume)indi atethattypi al ontemporary magneti eldsinspiralgalaxiesareoftheorderoffewuptofewtensofG,whi hmeans that magneti eldshaveundergoneampli ationbyatleastfourordersofmagnitudes withinthegala ti lifetime.Therefore,amodelforeÆ ientmagneti eldampli ation duringgala ti evolutionisne essary.Thestandardmodelofmagneti eldampli ation indiskgalaxiesisbasedonthetheoryofturbulentmeanelddynamo(seeWidrow2002 forare entreview).
Themean{elddynamotheoryhasbeensu essfulinexplainingmagneti eldsin var-iousastrophysi alobje ts,however,the lassi al,kinemati dynamo seemsto berather slowingalaxies.Magneti eldampli ation times alet
dynamo
(0:51)10 9
too longto explain 1G magneti elds ingalaxies whi h are only few 10 9
yrs old (z 13) (Wolfe,Lanzetta &Oren1992, Kronbergetal.2008). These ir umstan es seem to indi ate a need of an alternativeapproa h.We suggest that dire tnumeri al MHDmodelingofISMdynami s,through onsiderationsofdetailedphysi alpro esses, involvingmajor ISM omponents:gas, magneti eldsand osmi rays,is ne essaryto followmagneti historyofgalaxies.
2. Cosmic Rays in the interstellar medium
The dynami alrole of CRswas re ognized forthe rst time by Parker (1966), who noti edthataverti allystratiedISM, onsistingofthermalgasmagneti eldandCRs is unstable due to buoyan y of the weightless omponents: magneti elds and CRs. The CR omponentappears to be an important ingredient of this pro ess. A ording to thediusivesho ka elerationmodelsCRsare ontinuouslysuppliedto ISMbySN remnants.Therefore,strongbuoyan yee ts duetoCRsareunavoidable.
Theories ofdiusivesho ka elerationpredi t thatabout10% ofthe 10 51
ergof theSN II explosionenergyis onvertedto theCRenergy(see e.g.Joneset al1998and referen estherein).Observationaldataindi atethatgas,magneti eldsandCRsappear inapproximateenergeti equipartition,whi h meansthatallthe three omponentsare dynami ally oupled.Moreover,numeri alexperimentsbyGia alone&Jokipii (1999) indi ate that CRs diuse anisotropi ally in the ISM, along the mean magneti eld dire tion.
3. Cosmic Ray transport
Toin orporateCRpropagationinMHD onsiderationsweusethediusion-adve tion equation(e.g.S hli keiser&Ler he1985)
e r t +r(e r v)= p r rv+r( ^ Kre r ) (3.1) +CR sour es(SNremnants)
where the diusion term in written in the tensorial form to a ount for anisotropi diusivity of CRs. The sour es of CRs on the rhs. of Eqn. 3.1 orrespond to the CR produ tioninSupernovaremnants.
Thedynami sofISMin ludingtheenergeti sofCRs anbenowdes ribedbytheCR transportEqn. (3.1) andthe system of MHD equations.Wenote that inthe presen e of osmi rays additionalsour e term: rP
CR
shouldbe in ludedinthe gasequation ofmotion (seee.g.Berezinskiet al.1990),inordertoin orporatetheee tsof CRson gasdynami s.The osmi raydiusion{adve tionEqn.(3.1)hasbeensupplementedto theMHDalgorithmofZeus{3DMHD odebyHanasz&Les h(2003a)withtheaimof studyingISMphysi alpro esses,e.g.Parkerinstability,inwhi hCRsplayadynami ally importantrole.
4. CR-driven dynamo
CR{drivendynamohasbeenraisedbyParker(1992).Ourrstlo al(shearing{box)CR{ drivendynamonumeri alexperiments(Hanaszetal.2004,2006;Otmianowska-Mazuret al.2007)relyonthefollowingingredients:(1.)the osmi ray omponentdes ribedbythe diusion{adve tiontransportequation (3.1),in ludinglo alizedsour esof osmi rays | supernovaeremnants,explodingrandomlyinthediskvolume,whileSNsho ksand thermalee tsarenegle ted(seeGresseletal.2008a,2008b)forare entdynamomodel relyingonthethermalenergyoutputfromSNe);(2.)resistivityoftheISM(seeHanasz, Otmianowska{MazurandLes h2002),leadingtotopologi alevolutionofmagneti elds; (3) shearing boundary onditions (Hawley, Gammie and Balbus 1995) together with Coriolisand tidalfor es, aimed at modelingof dierentiallyrotatingdisksinthe lo al approximation;(4)realisti verti aldiskgravityandrotationfollowingthemodelofISM intheMilkyWaybyFerriere(1998).
In the most re ent paper (Hanaszet al. 2009) we presenta parameterstudy of the CR{drivendynamomodelbyexaminingthedependen eofmagneti eldampli ation on magneti diusivity, supernovae rate determining the CR inje tion rate, temporal modulation of SN a tivity, grid resolution, and CR diusion oeÆ ients. We nd the dominatingin uen eofmagneti diusivity(treatedasa freeparameter),amongother parameters,(Fig.1)ontheeÆ ien yof magneti eldampli ation.
Figure 1.
Timeevolutionofazimuthalmagneti uxandtotalmagneti energyfordierent valuesofmagneti diusivityinsimulationseriesA.The urvesrepresent respe tively asesof =0(A1),=1(A2),=10(A3),=100(A4)and=1000(A5)inunits p2 Myr
1 .
Thefastest magneti eldampli ation,observedintheexperiments, oin ideswith the magneti diusivity omparable to '
1 3 turb (obs), where turb (obs) ' 1 3 100p 10kms 1 ' 10 26 ms 1
. The question of whi h physi alpro ess an be responsible forthat largemagneti diusivity oeÆ ientspointsouttowardsthe urrentstudieson magneti re onne tion(Lazarianetal.2004,Cassaket al.2006)
Toinvestigateobservationalpropertiesof urrentdynamomodelsweperformaseries of numeri al simulations, pla ing the lo al shearing{boxes at dierent gala to entri radii. We repli ate the ontents of lo al ubes into rings, and then ombine rings to build upa disk(Otmianowska{Mazuret al.2009).A syntheti polarizationradio{map of syn hrotronemission onstru ted, on thebase of simulatedmagneti elds andCR
Figure 2.
Syntheti radio{mapofpolarizedradioemissionfor ourmodelatthein lination i=80Æ
andamapofrealgalaxyNGC5775
Apolarizedemissionmaptogetherwithsuperposedpolarizationve torsofsyn hrotron emission of an edge{on galaxy NGC 5775 is shown for omparison. As it is apparent the X{shaped stru turesof magneti eld, presentinthe realgalaxy,formalso inthe syntheti radioimage.ThepresentedresultsdemonstratethattheCR{MHDsimulations provideawayofveri ationofCR{dynamomodelsagainsttheobservationaldata.
Tosummarizetheoverallout omeofthelo al simulationsofCR{drivendynamo,we notethatthemodelprovideseÆ ientampli ationoflarge{s alemagneti elds,whi h attain,inthesaturatedstate,equipartitionwithturbulentgasmotions.However,wenote alsoadrawba kofthelo almodels,sin etheCRenergydensityremainsmorethanan order of magnitude larger than gas and magneti energy densities. We re ognize that improvementsofthe urrentmodelsarene essary,sin ethepossiblereasonoftheex ess of osmi rays may resultfrom theee t of trapping of osmi raysby predominantly horizontalmagneti eldinthe omputationalboxwhi hisperiodi inhorizontal dire -tions. Therefore,globalgala ti disksimulationsshouldbeperformed toenablees ape ofCRsalongthepredominantlyhorizontalmagneti eldinthegala ti plane,inorder to avoidtheCRex ess.
5. Global disk simulations
Re entlywestartedanewseriesofglobalsimulationsofCR{drivendynamowiththe aid of PIERNIKMHD ode (seeHanasz et al.2009a, 2009b, 2009 , 2009d),whi h is a grid{basedMPIparallelized,resistiveMHD ode basedontheRelaxingTVD (RTVD) s hemeby Jin& Xin(1995) andPen et al.(2003). Theoriginals heme isextended to dealwiththediusiveCR omponent(seeHanasz&Les h2003a).
In Fig. 3 we show results of one of the rst semi{globalsimulations of CR{driven dynamo.Thesimulationshavebeenperformedforaquarterofgala ti diskresembling MilkyWay,withparallelCRdiusion oeÆ ientK
k =310 28 m 2 s 1 ,ina omputa-tionaldomain20kp 20kp 20kp ,andresolutionof400400160grid ells,in 16MPIblo ks.
Figure 3.
Horizontal and verti al sli es through the omputational box for a simulation of CR{drivendynamoinasemi{globalgala ti disk.The olormapsre e t:1:gasdensity,2:CR energydensity,3:magneti eldstrength,4.theratioofmagneti toCRenergydensitiesand5. syn hrotronemissivityforasaturatedstateofthedynamo.Thegray{s alesinallpanels,ex ept 4.are hosenarbitrarily.Thedarkershades orrespond tolargervaluesofdisplayedquantities. Thegreys aleofpanel4.issu hthatbla kareasrepresentvaluesofemag=e r'1.eld: growth time of the mean magneti ux is approximately 210 Myr in the whole galaxy,whi his loseto thegala ti rotationperiodneartheorbitofSun.Asexpe ted, the globalsimulationsintrodu e a signi antimprovementwith respe t to the former shearing{boxsimulations.It is apparentinFig. 3that theratioof magneti to osmi ray energy is lose to one in the syn hrotronemitting volume( ompare panels 4.and 5.), whileover{equipartition of CRenergy with respe t to magneti energy still holds outsidethedisk.
ThelatestdevelopmentofourCR{dynamomodelisafullyglobalgala ti disk simula-tion(seethe omplementarypaperbyHanasz,etal.thisvolume),wherewedemonstrate thatdipolarmagneti eldssuppliedonsmallSN{remnants ales, anbeamplied expo-nentiallybytheCR{drivendynamotothepresentequipartitionvalues,andtransformed simultaneouslyto largegala ti {s ales.
6. Conclusions
WehaveshownthattheCR ontributiontothedynami sofISM,studiedbymeansof CR{MHDsimulations,inbothlo alandglobals ales,leadstoaveryeÆ ientmagneti eldampli ationingala ti disks.
The osmi ray driven dynamo amplies eÆ iently gala ti magneti eldson a ti-mes aleofgala ti rotation.Theresultinggrowthofthelarge{s alemagneti eldby4 ordersofmagnitudewithin 2Gyr(Hanaszatal.2004),isfastenoughtoexpe t1G magneti eldingalaxies atz13.
Wepointouttheadvantageofglobaldiskmodels,whi honthe ontrarytoshearing{ box models, provide CR energy equipartition with magneti eld in the syn hrotron emittingpartofthedisks.
A knowledgements
ThisworkwassupportedbyPolishMinistryofS ien eandHigherEdu ationthroughthe grants92/N{ASTROSIM/2008/0and PB0656/P03D/2004/26and byNi olaus
Coper-References
Be k R. 2008, in Cosmi Magneti Fields: From Planets, to Stars and Galaxies, K.G. Strassmeier,A.G. Kosovi hev &J.E. Be kman,eds., Pro .IAUSymp.259 (in press, arXiv:astro-ph/0812.4925)
Berezinskii,V.S., Bulanov, S.V.,Dogiel,V.A., Ginzburg, V.L.,Ptuskin,V.S. Astrophysi s of osmi rays,Amsterdam:North-Holland,1990.
Cassak,P. A.,Drake,J.F.,Shay, M.A.2006, ApJ644, L145 Gia alone, J., Jokipii,R.J.1999, ApJ520, 204
Ferriere,K. 1998,ApJ497, 759
Gressel,O.,Ziegler, U., Elstner,D., R udiger,G. 2008,AN,329, 61 Gressel,O.,Elstner, D., Ziegler,U., R udiger,G. 2008,A&A,486L, 35 Hanasz,M., Otmianowska{Mazur, K.,Les h,H.2002,A&A 386, 347 Hanasz,M., Les h, H.2003,A&A412, 331
Hanasz,M., Kowal, G.,Otmianowska{Mazur, K.,&Les h,H.2004, ApJ605, L33 Hanasz,M., Kowal, G.,Otmianowska{Mazur, K.,&Les h,H.2006, AN327, 469
Hanasz, M., Otmianowska{Mazur, K., Kowal, G., & Les h, H. 2009, A&A (in press, arXiv:astro-ph/0812.3906)
Hanasz,M., Kowalik, K., Woltanski, D., Paw laszek, R.K. 2009,in:K. Gozdziewskietal. (eds.) Extrasolar planets in multi{body systems:theory and observations,EAS Publi- ationsSeries(submitted,arXiv:astro-ph/0812.2161)
Hanasz,M., Kowalik, K., Woltanski, D., Paw laszek, R.K. 2009,in:K. Gozdziewskietal. (eds.) Extrasolar planets in multi{body systems:theory and observations,EAS Publi- ationsSeries(submitted,arXiv:astro-ph/0812.2799)
Hanasz, M., Kowalik, K., Woltanski, D., Paw laszek, R. K., 2009, in: M. de Avillez et al. (eds.), The Role of Disk{Halo Intera tion in Galaxy Evolution: Out ow vs Infall?, (submitted,arXiv:astro-ph/0812.4839)
Hanasz, M., Kowalik, K., Woltanski, D., Paw laszek,R. K., 2009, in: M. de Avillez et al. (eds.), The Role of Disk{Halo Intera tion in Galaxy Evolution: Out ow vs Infall?, (submitted,arXiv:astro-ph/0901.0104)
Hawley,J.F.,Gammie,C.F.,Balbus,S.A. 1995,ApJ440,442
Jones, T.W., Rudni k, L.,Jun,B.{I.; Borkowski,K. J.; Dubner, G.,Frail,D. A., Kang, H., Kassim,N.E., M Cray,R.1998,PASP 110,125
Jin, S.,Xin, Z.1995,Comm.PureAppl. Math.48,235
Kronberg,P.P., Bernet,M.L.,Miniati, F.,Lilly,S.J.,Short,M.B.,Higdon,D.M.2008, ApJ676, 70
Lazarian,A.,Vishnia , E.T.,Cho,J.2004,ApJ603, 180L
Otmianowska-Mazur, K.,Kowal,G.,Hanasz,M., 2007,ApJ, 668,1100 Otmianowska{Mazur,K.,Soida,M.,Kulesza{
_
Zydzik,B.,Hanasz,M.,Kowal,G.2009,ApJ (in press,arXiv:astro-ph/0812.2150)
Parker, E.N.1966,ApJ145,811 Parker, E.N.1992,ApJ401,137
Pen,U.{L., Arras,P., Wong,S. 2003,ApJS149, 447 Rees,M.J.1987, QJRAS28,197
Ryu,D., Kim, J.,Hong, S.S.,Jones,T.W. 2003,ApJ589, 338 S hli keiser, R.&Ler he,I.1985,A&A 151,151
Widrow,L.M. 2002,Rev.Mod.Phys. 74,775