这件事说起来容易做起来难。
要得到光度的分布，只有测光是不够的，还需要获得距离信息。而距离信息是最难的目标。去年我的工作中用BHB恒星定标距离会产生１０％的误差，再加上自行的不确定度，难度可以想象。现在，得到光度函数有２个途径。１是经验模型，比如利用star counts，这种方法用于银晕只能通过球状星团获得距离，用于银盘只能通过太阳临域有视差的恒星来获得距离。2是理论模型，通过恒星演化理论计算理论演化轨迹。
但是这里面有很多非常不靠谱的事情。比如恒星演化理论，针对大质量星没有很好的结果；恒星计数方法针对不同光谱型的恒星会得到不同的星族成分的scale height，差别很大；太阳临域有视差观测的样本很有限，其color-M关系外推的误差相当大；球状星团的光度函数似乎也存在差异，而核心附近的成员星又无法观测到。。。
根本上说，光度函数是一个难点和终点，围绕着这个话题天文学家争论了一个世纪了，今天还在争论。相对好的做法是把1、2两个途径结合起来，也就是所谓的semi-empirical model，两家取长补短，用恒星演化理论模拟光度，用星团和恒星计数结果模拟极端恒星样本。目前读到的两篇文献，都是这样的处理思路。一篇是台湾中央大学天文所2000年的博士论文，另一篇是TRILEGAL这个基于web的模拟程序的描述，PI是意大利的。
我利用后一篇文章的程序简单演示了一下前一篇文章中的结果，思路基本上整理顺了。这里要感谢所里的CTO，Peter先生，给他发邮件只有用一句话，他就帮我搞定了一切安装MATLAB的过程，甚至调试好了JAVA平台，直接给我弄到桌面上，让我可以轻松使用正版MATLAB 2007a系统。这件事要是自己捣鼓，一个星期不一定搞得定，而且还不合法。
下一篇blog打算详细介绍一下ARI研究所的情况，顺便感谢几位几天之内就给我极大帮助的先生。
待续。。。
Recent Comments