Ngc3603的MSX數(shù)據(jù)分析ppt課件

1、Ngc3603的MSX數(shù)據(jù)分析紅外數(shù)據(jù)處置和分析THANK YOU ALL!New year 曙光擬和勝利收到helpdesk初步回應(yīng),the first valuable response in 2 years!主要內(nèi)容數(shù)據(jù)處置: 一波三折擬和結(jié)果和新嘗試塵埃粒子數(shù)密度計(jì)算(進(jìn)展中)第一次修正自創(chuàng)X射線處置方法(BingJiang),把四個(gè)波段疊加后取背景根據(jù)背景的含義,運(yùn)用區(qū)域內(nèi)而不是區(qū)域外的.(discuss with BingJiang)挖點(diǎn)源后再取背景單溫度灰體擬和結(jié)果溫度約30K.x2 在(4-3)=1附近.非常好的結(jié)果.致命錯(cuò)誤!Thanks to BingJiang! 擬和代碼

2、中指數(shù)部分的單位同源自數(shù)據(jù)的scale的單位混淆,導(dǎo)致e指數(shù)上偏向10倍! exp(50u/t) - exp(500u/t) u: 1e13HZ,t: K擬和結(jié)果: T500K!雙溫灰體譜擬和確定大致溫度先(YangChen).Super Shell(YangSu):Jon M.Saken, J.Michael Shull,Catharine D. Garmany,APJ,397:537-541,1992IRAS: Jon M.Saken, Robert A.Fesen, AND J.Michael Shull,APJ Supplement Series,81:715-745,1992 -雙溫

3、灰體mathematica擬和公式:F(u)=V1*u*B(u,T1)+V2*u*B(u,T2)問題:1.怪異的0自在度擬和.2.什么是多維X2方法?3.實(shí)際推導(dǎo):雙溫模型的數(shù)學(xué)猜測(cè)(略)?擬和結(jié)果: (FindMinimum)高溫t2高達(dá)356K! !(0.6644092449814092, v1 - 6.391802841729306*9, v2 - 2.8711509228537344, t1 - 30.125079121945603, t2 - 356.11816589614733) 第二次修正系數(shù)修正:Pixels-立體角立體角計(jì)算錯(cuò)誤: fits頭中的CDELTA運(yùn)用的是度而不是弧

4、度. (/180*PI)2背景計(jì)算: 不計(jì)入輻射強(qiáng)度為0區(qū)的背景.(對(duì)結(jié)果只需千分之一量級(jí)的影響)發(fā)現(xiàn)經(jīng)典文章S.K.Ghosh and D.K.Ojha,A&A 388,326-334(2002) - methodKathleen E. Kraemer, Russell F. Shipman, Stephan D. Price,1 Donald R. Mizuno,Thomas Kuchar,3 and Sean J. Carey,AJ,126:1423-1450,2003 - summaryStephan D.Price,Michael P.Egan,Sean J.Carey,Donald

5、 R.Mizuno,and Thomas A.Kuchar, AJ,121:2819-2842,2001 - msxKathleen E.Kraemer et alSeveral sources of emission can be expected to contributeto the MSX bands from the star-forming regions observedwith SPIRIT III. Ghosh & Ojha (2002) have been able tomodel the band A (6.810.8 lm) and band C (11.113.2 l

6、m)emission from five compact H ii regions as a combination ofthermal (graybody) dust emission and emission from theunidentified infrared bands at 7.78.6 and 11.2 lm typicallyidentified with polycyclic aromatic hydrocarbons (PAHs).Bands D (13.515.9 lm) and E (18.225.1 lm) are dominatedby the rising c

7、ontinua of cool dusty objects. All fourbands contain forbidden lines from Ar and Ne, amongothers, although analyses by Ghosh & Ojha (2002) andVerstraete et al. (2001) suggest that their contribution to thein-band radiance will be negligible compared with the PAHsand the continuum emission.新擬和的系數(shù)頻率相關(guān)

8、的系數(shù)(灰體譜的系數(shù))選取 根據(jù)Z.-R.Wang,R.McCray,Y.Chen,and Q.-Y.Qu Astron.Astrophys.240,98-104(1990) 取頻率的1次方,不思索消光. Ghosh思索消光. 有人取1.5次方(Jon M.Saken et al).擬和公式(參考S.K.Ghosh et al.)選定區(qū)域(西北起逆時(shí)針adbc)選定區(qū)域(西北起逆時(shí)針adbc).上:AC,下:DEMathematica擬和(環(huán)形)(BingJiang,XinLianLuo,YingHeZhao)擬和結(jié)果:環(huán)形(129.55K)A(119.384K)B(129.837K)C(13

9、4.742K,X22)D(124.031K)流程圖數(shù)據(jù)處置中的遺留問題誤差計(jì)算: 1.helpdesk.-從相對(duì)誤差計(jì)算是對(duì)的.waiting expert 2.平均誤差含義猜測(cè).(BingJiang). Sqrt(Sum(I*e)2)/N.Then N=what? 3.樣本?個(gè)體?點(diǎn)源挖取的控制.helpdesk|wavelet|形狀學(xué)嘗試:除法(溫度的空間分布)實(shí)際根底猜測(cè): d,e波段是延續(xù)譜,由于溫度約100k,頻率1-3(1013Hz),所以exp(500u/T)1,流量近似表示為: f=Vu4exp(-500u/T).取定空間坐標(biāo)后作除法: f(d)/f(e) = (ud/ue)4

10、exp(-500(ud-ue)/T).只需T是局域變量. 可見對(duì)每個(gè)pixel做除法后，可以得到溫度的空間分布問題: 1.誤差棒?2.除法的真正實(shí)際根據(jù)?除法的數(shù)據(jù)處置CDELTA-6角秒/pixel,但空間分辨率18.3角秒(Price et al. 2001), - 做3x3的bin得到有效數(shù)據(jù). dmcopy(BingJiang)限制條件. 1.根據(jù)exp(500u/t)1的條件,限制-0f(d)/f(e) 1(T200k). 其他點(diǎn)都是無效點(diǎn)置0. 2.根據(jù)d,e單溫度擬和的結(jié)果(K),過濾溫度在120-150k之間的塵埃.除法的結(jié)果(ds9平滑,未去背景)左:120-150k,右:

11、logF =-500u/T+logV+4logu消去溫度: log(F(ud)ue-log(F(ue)ud其中重要的項(xiàng)是logV(ud-ue),可以得到logV的空間分布信息. Errorbar: dlogF=1/FdF.要求F獲得大一些errorbar才有意義. Bin?Ngc3603的尺度NGC3603 is located in the Carina spiral arm at a distance about 7kpc (D.E.A.Nurnberger A&A 404,255-265,2003)取直徑1度,可估算尺度100pc,符合supershells定義(Jon M.Saken et al)顆粒計(jì)算(Z.R.Wang et al)記f(nu,T)=nu*B(nu,T)F(nu) = V/(d2)*J(nu) J(nu) = n*a*nu*B(nu,T) = n*a* f(nu,T)故F(nu) = V/(d2)*n*a*f(nu,T)不思索單位,V/(d2)*n*對(duì)應(yīng)擬和結(jié)果中的系數(shù)V.a的值可以確定(Z.R.Wang et al),d取7kpc,只需求出塵埃區(qū)域體積V就可以得到塵埃顆粒的數(shù)密度n.體積的計(jì)算環(huán)形區(qū)域體積由于對(duì)稱,可以簡化為簡單幾何體的組合.矩形區(qū)域體