LAMOST恒星参数测量pipeline的改进课件

LAMOST的恒星参数LAMOST数据部2014用户培训大会2014/7/2LAMOST的恒星参数LAMOST数据部报告内容低分辨率光谱恒星参数的测量方案恒星参数的测量过程SNR对测量精度的影响LAMOST光谱的一致性

误差估计报告内容低分辨率光谱恒星参数的测量方案低分辨率光谱恒星参数的测量方案光谱特征（谱线）测量＋经验关系？连续谱的确定？低分辨率光谱恒星参数的测量方案光谱特征（谱线）测量＋经验关系低分辨率光谱恒星参数的测量方案特征测量在不同分辨率下系统差必然存在模版匹配对个别缺失谱线敏感程度低交叉相关匹配的视向速度测量不仅利用线心而是整个谱线上点的特征模版选择是理论的还是经验的？（不确定性？参数范围？）选择在模版网格中进行怎样的内插方式？交叉相关局部最有化问题和随机误差的确定低分辨率光谱恒星参数的测量方案特征测量在不同分辨率下系统差必LAMOST光谱的参数测量方案一致性？（空间和时间）及变形缺失特征

交叉相关寻找最优解是必然选择只利用谱线的交叉相关而去掉伪连续谱Ulyss是以Elodie为模版训练出来的交叉相关及内插的非线性内插程序选择初始值避免收敛到局部最有解LAMOST光谱的参数测量方案一致性？（空间和时间）及变形缺AFGK恒星参数的测量过程概述测量参数用的光谱特征：4000A～6800A范围内的谱线测量过程（pipeline）：1、CFI1(与NGS类似)2、Ulyss13、去除Continuum4、CFI25、Ulyss26、决策：比较Ulyss2和Ulyss1的结果去掉3sigma外的点，三个参数都一致的保留AFGK恒星参数的测量过程概述测量参数用的光谱特征：4000参数的初始假设

CFI1类似于SSPP中模板匹配的模块NGS，用其9阶多项式归一化的模板（KuruczNEWODF)。对观测光谱也采用了9阶多项式迭代计算连续谱。计算与模板的相关函数系数，按照Teff->[Fe/H]->logg的顺序计算最相关的模板作为初始假设CFI1参数的初始假设CFI1类似于SSPP中模板匹配的模块NG初始假设CFI1的可靠性初始假设CFI1的可靠性初始假设CFI1的补充检验鉴于SDSS模版（36条光谱）存在的问题，我们用LAMOST高银纬恒星的数据构建了LAMOST（183条光谱）的分类模版库，改正了1Dpipeline的分类错误（2014,AJ147,101）纠正了v2.6.4版本将G型星错分成F9，M型星错分成G或K型恒星等问题;可以作为参数初值CFI的检验初始假设CFI1的补充检验鉴于SDSS模版（36条光谱）存在利用ulyss测量参数构建非线性化模型函数TGM，可生成任意给定Teff,logg,[Fe/H]的恒星谱，最小化待测谱和模型谱间的拟合χ²值。构建内插模型，选用ELODIE光谱库2000多条光谱，按温度分为热OBA,温FGK,冷M三个区间，光谱的每个点由以Teff,logg,

[Fe/H]为变量的19～26项组合多项式近似表征。多项式各项的组合形式和温度划分边界及权重分配系数，都是通过在对实测谱和内插生成的谱之间的残差做细致调整而定。各项间正交，避免相互抵消，依贡献是否显著逐渐填加新项，残差是否降低。为了保证收敛到全局最优解，采用CFI进行限制，同时也减少ULySS的运算量和时间。而无需在参数空间分布中设定多个起始迭代点。以CFLIBstarHD4307为例：收敛域检测。利用ulyss测量参数利用Ulyss1结果去连续谱伪连续谱应由天文学家确定，才能保证保留所需参数而海量光谱很难完成，也无法做到统一一致折衷方案：1、计算出了一套网格更细的Kurucz，采用宋轶晗开发的类似于iraf的基于python的continuum软件人工画出所有连续谱2、用Ulyss1的结果找一个参数最接近的理论光谱，将此光谱叠加一个5阶多项式与理论光谱进行拟和，该理论光谱对应的人工划出的连续谱乘以5阶多项式就是观测光谱的连续谱利用Ulyss1结果去连续谱伪连续谱应由天文学家确定，才能保去除连续谱的例子去除连续谱的例子迭代：CFI2／Ulyss2对归一化后的光谱再次计算初始假设cfi利用初始假设，再次运行Ulyss迭代：CFI2／Ulyss2对归一化后的光谱再次计算初始假设迭代前后决策：迭代前后，如果参数结果差别很小，认为参数以及通过参数所找到的模板是正确的。参数是可信的。如果迭代前后，差别很大，可能会跨跃几个子类型，认为参数不可靠，人眼检查。对LAMOST的DR1光谱，迭代处理后，三个参数的3sigma外的点取并集，由人检查其参数确定的准确性。迭代前后决策：迭代前后，如果参数结果差别很小，认为参数以及通.迭代前后三参数的对比.迭代前后三参数的对比信噪比对参数的影响16信噪比对参数的影响16LAMOST光谱的一致性问题选源为同源的谱线不同LAMOST光谱的一致性问题选源为同源的谱线不同其它独立工作对LAMOST参数的对比WangW（ApJSubmitted）

比较LASP和高分辨光谱确定的参数(Brunttetal.2012;Thygesenetal.2012;Molenda-Żakowiczetal.2013),Teff：36±108K[Fe/H]：0.02±0.11dex比较logg和星震学测量的logg,0.02±0.13dex,只有3%的outliers.其它独立工作对LAMOST参数的对比WangW（ApJS其它独立工作对LAMOST参数的对比Dong,S（ApJL）

比较LASP和Buchhave高分辨光谱用Torresetal.(2012)方法确定的金属丰度[Fe/H]：0.02±0.015dex其它独立工作对LAMOST参数的对比Dong,S（ApJL其它独立工作对LAMOST参数的对比Gao,H.,Zhang,H.W.（RAASubmitted）

比较LASP和PASTELCatalog中的高分辨光谱确定的参数Teff:11±110Klogg:0.03±0.19dex[Fe/H]:0.01±0.11dexrv:-3.78±4.91km/s其它独立工作对LAMOST参数的对比Gao,H.,Zha其它独立工作对LAMOST参数的对比XiangM.S.（MNRASwillbesubmitted）Teff:9±135Klogg:0.02±0.23dex[Fe/H]:0.06±0.13dexrv：1.2±3.1km/s其它独立工作对LAMOST参数的对比XiangM.S.（M下一步工作扩充模板的参数空间范围扩充波长测量范围——考虑采用经过高分辨进行参数定标后的LAMOST观测的实测模版下一步工作扩充模板的参数空间范围

谢谢！lal@wuyue@dubing@yhsong@

愿孔啸早日康复！谢谢！LAMOST的恒星参数LAMOST数据部2014用户培训大会2014/7/2LAMOST的恒星参数LAMOST数据部报告内容低分辨率光谱恒星参数的测量方案恒星参数的测量过程SNR对测量精度的影响LAMOST光谱的一致性

误差估计报告内容低分辨率光谱恒星参数的测量方案低分辨率光谱恒星参数的测量方案光谱特征（谱线）测量＋经验关系？连续谱的确定？低分辨率光谱恒星参数的测量方案光谱特征（谱线）测量＋经验关系低分辨率光谱恒星参数的测量方案特征测量在不同分辨率下系统差必然存在模版匹配对个别缺失谱线敏感程度低交叉相关匹配的视向速度测量不仅利用线心而是整个谱线上点的特征模版选择是理论的还是经验的？（不确定性？参数范围？）选择在模版网格中进行怎样的内插方式？交叉相关局部最有化问题和随机误差的确定低分辨率光谱恒星参数的测量方案特征测量在不同分辨率下系统差必LAMOST光谱的参数测量方案一致性？（空间和时间）及变形缺失特征

交叉相关寻找最优解是必然选择只利用谱线的交叉相关而去掉伪连续谱Ulyss是以Elodie为模版训练出来的交叉相关及内插的非线性内插程序选择初始值避免收敛到局部最有解LAMOST光谱的参数测量方案一致性？（空间和时间）及变形缺AFGK恒星参数的测量过程概述测量参数用的光谱特征：4000A～6800A范围内的谱线测量过程（pipeline）：1、CFI1(与NGS类似)2、Ulyss13、去除Continuum4、CFI25、Ulyss26、决策：比较Ulyss2和Ulyss1的结果去掉3sigma外的点，三个参数都一致的保留AFGK恒星参数的测量过程概述测量参数用的光谱特征：4000参数的初始假设

CFI1类似于SSPP中模板匹配的模块NGS，用其9阶多项式归一化的模板（KuruczNEWODF)。对观测光谱也采用了9阶多项式迭代计算连续谱。计算与模板的相关函数系数，按照Teff->[Fe/H]->logg的顺序计算最相关的模板作为初始假设CFI1参数的初始假设CFI1类似于SSPP中模板匹配的模块NG初始假设CFI1的可靠性初始假设CFI1的可靠性初始假设CFI1的补充检验鉴于SDSS模版（36条光谱）存在的问题，我们用LAMOST高银纬恒星的数据构建了LAMOST（183条光谱）的分类模版库，改正了1Dpipeline的分类错误（2014,AJ147,101）纠正了v2.6.4版本将G型星错分成F9，M型星错分成G或K型恒星等问题;可以作为参数初值CFI的检验初始假设CFI1的补充检验鉴于SDSS模版（36条光谱）存在利用ulyss测量参数构建非线性化模型函数TGM，可生成任意给定Teff,logg,[Fe/H]的恒星谱，最小化待测谱和模型谱间的拟合χ²值。构建内插模型，选用ELODIE光谱库2000多条光谱，按温度分为热OBA,温FGK,冷M三个区间，光谱的每个点由以Teff,logg,

[Fe/H]为变量的19～26项组合多项式近似表征。多项式各项的组合形式和温度划分边界及权重分配系数，都是通过在对实测谱和内插生成的谱之间的残差做细致调整而定。各项间正交，避免相互抵消，依贡献是否显著逐渐填加新项，残差是否降低。为了保证收敛到全局最优解，采用CFI进行限制，同时也减少ULySS的运算量和时间。而无需在参数空间分布中设定多个起始迭代点。以CFLIBstarHD4307为例：收敛域检测。利用ulyss测量参数利用Ulyss1结果去连续谱伪连续谱应由天文学家确定，才能保证保留所需参数而海量光谱很难完成，也无法做到统一一致折衷方案：1、计算出了一套网格更细的Kurucz，采用宋轶晗开发的类似于iraf的基于python的continuum软件人工画出所有连续谱2、用Ulyss1的结果找一个参数最接近的理论光谱，将此光谱叠加一个5阶多项式与理论光谱进行拟和，该理论光谱对应的人工划出的连续谱乘以5阶多项式就是观测光谱的连续谱利用Ulyss1结果去连续谱伪连续谱应由天文学家确定，才能保去除连续谱的例子去除连续谱的例子迭代：CFI2／Ulyss2对归一化后的光谱再次计算初始假设cfi利用初始假设，再次运行Ulyss迭代：CFI2／Ulyss2对归一化后的光谱再次计算初始假设迭代前后决策：迭代前后，如果参数结果差别很小，认为参数以及通过参数所找到的模板是正确的。参数是可信的。如果迭代前后，差别很大，可能会跨跃几个子类型，认为参数不可靠，人眼检查。对LAMOST的DR1光谱，迭代处理后，三个参数的3sigma外的点取并集，由人检查其参数确定的准确性。迭代前后决策：迭代前后，如果参数结果差别很小，认为参数以及通.迭代前后三参数的对比.迭代前后三参数的对比信噪比对参数的影响39信噪比对参数的影响16LAMOST光谱的一致性问题选源为同源的谱线不同LAMOST光谱的一致性问题选源为同源的谱线不同其它独立工作对LAMOST参数的对比WangW（ApJSubmitted）

比较LASP和高分辨光谱确定的参数(Brunttetal.2012;Thygesenetal.2012;Molenda-Żakowiczetal.2013),Teff：36±108K[Fe/H]：0.02±0.11dex比较logg和星震学测量的logg,0.02±0.13dex,只有3%的outliers.其它独立工作对LAMOST参数的对比WangW（ApJS其它独立工作对LAMOST参数的对比Dong,S（ApJL）

比较LASP和Buchhave高分辨光谱用Torresetal.