1 / 23

NGC3603 的红外辐射疑难

NGC3603 的红外辐射疑难. 尘埃密度温度计算中存在的问题. 内容. 体积计算中的错误 原模型中存在的问题 何去何从 ? 几个可能的方案 . 当前新模型及结果 下一步. 体积计算中的错误 1. 球环截球区域的体积 : 其中截面的内外边界的高度. 体积计算中的错误 2. 矩形区域计算 : 原来的计算中存在积分方向错误. 体积计算中的错误 2. 三种新计算方法 : 1.c 程序 2D 积分 . 步长为 pixel.z 方向可以通过球面方程用其他两维坐标表示 . 2.YangChen 近似 .z 方向直接用 ring 区域半径 r1,r2.

howard-avery
Download Presentation

NGC3603 的红外辐射疑难

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. NGC3603的红外辐射疑难 尘埃密度温度计算中存在的问题

  2. 内容 • 体积计算中的错误 • 原模型中存在的问题 • 何去何从?几个可能的方案. • 当前新模型及结果 • 下一步

  3. 体积计算中的错误1 • 球环截球区域的体积: 其中截面的内外边界的高度

  4. 体积计算中的错误2 矩形区域计算: 原来的计算中存在积分方向错误

  5. 体积计算中的错误2 三种新计算方法: 1.c程序2D积分.步长为pixel.z方向可以通过球面方程用其他两维坐标表示. 2.YangChen近似.z方向直接用ring区域半径r1,r2. 3.BingJiang近似.侧向补全法.截面*(z半短轴+z半长轴).

  6. 原模型细节(ZRWang et al)

  7. PAH的贡献

  8. 温度估算

  9. 温度疑难 拟和的温度约120K>>16K!需要恒星温度升高一个量级才能弥补这种差异. 如果Q~波长的^-2,适当升高恒星温度到3倍,可能得到100K的尘埃温度(YangChen). 于是决定更换模型,使用2次方的关系.

  10. Q的计算 根据Draine 1989 APJ 285,89,8-70微米的辐射主要贡献来自于Si颗粒的贡献,同前,只考虑球形颗粒,从si的Q图( 文中FIG. 5b)取q:

  11. GIMP取点

  12. 区域(西北起逆时针adbc)

  13. 拟和结果T~114K(Ring)

  14. AB,CD区T降低~10K

  15. 气尘比疑难 • 根据Drain et al,取尘埃半径为0.1微米, 计算得到 同Z.R.Wang et al比较,小了6个数量级！同样,气尘比的结果注定荒谬! • 考虑到Q~a(Draine et al),式(1)中,尘埃数密度反比于a^3,尘埃半径降低2个量级可以抵消这种差异. • 但是,气尘比的计算中,数密度同半径3次方刚好是乘积,调整半径不会影响气尘比的量级!

  16. 新方案.设想 • 1.数据. 背景重新计算? MSX新数据?(BingJiang) e波段使用laplace算子? Su~1/mu模型质疑，用积分模型代替? • 2.体积 如果非球状呢? NGC3603是超级扁的? 形状采用小块拼接或乘上一个空间弥散系数,降低有效体积. • 3.Q 根据draine的曲线分段取Q~波长关系再拟和? • 4.100K尘埃红外辐射的文章? 150K et al. • 5.拟和 两温度模型+nu^2模拟? cstat方法cxc.harvard.edu/ciao/ahelp/cstat.html D,E加入PAH拟和

  17. 当前新模型 • 考虑到核心问题是温度,目标:降温 • 拟和模型,D,E也加入PAH辐射. • Q? 波长^-2? • 问题: PAH辐射的系数怎样确定?综合多个PAH光谱(YangChen)

  18. Verstraete et al.

  19. M.Jura et al.

  20. G.Mulas 2006 A&A 446,537

  21. 当前的拟和结果0.8+30K • rr[u_, a_, v_, t_] = • Norm[{(u + f[v, t][rawdata[[1, 1]]] - rawdata[[1, 2]])/ • rawdata[[1, 3]], (a*u + f[v, t][rawdata[[2, 1]]] - rawdata[[2, 2]])/ • rawdata[[2, 3]], (0.8*a*u + f[v, t][rawdata[[3, 1]]] - • rawdata[[3, 2]])/ • rawdata[[3, 3]], (0.1*a*u + f[v, t][rawdata[[4, 1]]] - • rawdata[[4, 2]])/rawdata[[4, 3]]}] • f12 = FindMinimum[rr[u, 1.39216, v, t], {u, 1}, {v, 1}, {t, 50}, • MaxIterations -> 1000000] • \!\({0.841040117453336`, {u -> 3.2963920906093045`, • v -> 5.343542450965798`*^9, t -> 30.384445183392366`}}\)

  22. 下一步 • 用GIMP取点 • 转换为流量/频率 • 综合得到 • 比较cstat拟和的结果 • 估算气尘比

  23. Thank you!

More Related