NGC3603的红外辐射疑难.ppt

NGC3603的红外辐射疑难 尘埃密度温度计算中存在的问题 内容 体积计算中的错误 原模型中存在的问题 何去何从?几个可能的方案. 当前新模型及结果 下一步 体积计算中的错误1 球环截球区域的体积: 其中截面的内外边界的高度 体积计算中的错误2 体积计算中的错误2 三种新计算方法: 1.c程序2D积分.步长为pixel.z方向可以通过球面方程用其他两维坐标表示. 2.YangChen近似.z方向直接用ring区域半径r1,r2. 3.BingJiang近似.侧向补全法.截面*(z半短轴+z半长轴). 原模型细节(ZRWang et al) PAH的贡献 温度估算 温度疑难 拟和的温度约120K16K!需要恒星温度升高一个量级才能弥补这种差异. 如果Q~波长的^-2,适当升高恒星温度到3倍,可能得到100K的尘埃温度(YangChen). 于是决定更换模型,使用2次方的关系. Q的计算 GIMP取点 区域(西北起逆时针adbc) 拟和结果T~114K(Ring) AB,CD区T降低~10K 气尘比疑难 根据Drain et al,取尘埃半径为0.1微米, 计算得到 同Z.R.Wang et al比较,小了6个数量级！同样,气尘比的结果注定荒谬! 考虑到Q~a(Draine et al),式(1)中,尘埃数密度反比于a^3,尘埃半径降低2个量级可以抵消这种差异. 但是,气尘比的计算中,数密度同半径3次方刚好是乘积,调整半径不会影响气尘比的量级! 新方案.设想 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拟和 当前新模型 考虑到核心问题是温度,目标:降温 拟和模型,D,E也加入PAH辐射. Q? 波长^-2? 问题: PAH辐射的系数怎样确定?综合多个PAH光谱(YangChen) Verstraete et al. M.Jura et al. G.Mulas 2006 AA 446,537 当前的拟和结果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`}}\) 下一步 用GIMP取点 转换为流量/频率 综合得到 比较cstat拟和的结果 估算气尘比 Thank you! * 矩形区域计算: 原来的计算中存在积分方向错误 根据Draine 1989 APJ 285,89,8-70微米的辐射主要贡献来自于Si颗粒的贡献,同前,只考虑球形颗粒,从si的Q图( 文中FIG. 5b)取q: *