Computation of the atmosphere-less light intensity curve during a total solar eclipse by using Lunar Reconnaissance Orbiter topography data and the DE430 astronomical ephemeris

Observations of the sky irradiation intensity in the visible wavelengths during a solar eclipse permit to model the Sun diameter,a key number to constrain the internal structure of our star.In this paper,we present an algorithm that takes advantage of the precise Moon topography from Lunar Reconnaissance Orbiter to compute,with a high resolution in time,the geometrical part(i.e.top-of-atmosphere,and for a given wavelength)of the sky irradiation at any given location on the Earth during these events.The algorithm is also able to model the Baily’s beads.We give as an application the theoretical computation of the light curve corresponding to the solar eclipse observed at Lakeland(Queensland,North Australia)on 2012 November 13.The application to real data,with the introduction of atmospheric and instrumental passbands,will be considered in a forthcoming paper.     

Investigating variable stars in the open cluster NGC 1912 and its surrounding field

In this work,we studied the variable stars in the open cluster NGC 1912 based on the photometric observations and Gaia DR2 data.More than 3600 CCD frames in B,V,R filters were reduced,and we obtained the light curves that span about 63 hours.By analyzing these light curves,we detected 24 variable stars,including 16 periodic variable stars,seven eclipsing binaries and one star whose type is unclear.Among these 24 variable stars,11 are newly discovered,which are classified as sixγDoradus stars,oneδScuti star,three detached binaries and one contact binary.We also confirmed 13 previously known variable stars.Based on cluster members identified by Cantat-Gaudin et al.(2018),we inferred cluster memberships for these detected variable stars.Using Gaia DR2 data,we plotted a new color-magnitude diagram for NGC1912,and showed the nature of variable cluster members in kinematical properties and heliocentric distance.Among the 24 variable stars,seven variables are probable cluster members,which show homogeneity in kinematic characters and space position with the established cluster members.Four of the seven variable cluster members are the previously discovered stars,consisting of twoγDor stars and twoδSct stars.The remaining three variable cluster members,which are allγDor stars,are firstly detected in this work.The main physical parameters of these variable cluster members estimated from the color-magnitude diagram are log(age/yr)=8.75,[Fe/H]=-0.1,m-M=10.03 mag,and E(B-V)=0.307.     

Detection of a high-confidence quasi-periodic oscillation in radio light curve of the high redshift FSRQ PKS J0805–0111

In this work,we have searched for quasi-periodic oscillations(QPOs)in the 15 GHz light curve of the FSRQ PKS J0805-0111 monitored by the Owens Valley Radio Observatory(OVRO)40 m telescope during the period from 2008 January 9 to 2019 May 9,using the weighted wavelet Z-transform(WWZ)and the Lomb-Scargle Periodogram(LSP)techniques.This is the first time to search for a periodic radio signal in the FSRQ PKS J0805-0111 by these two methods.Both methods consistently reveal a repeating signal with a periodicity of 3.38±0.8 yr(>99.7%confidence level).In order to determine the significance of the periods,the false alarm probability method was applied,and a large number of Monte Carlo simulations were performed.As possible explanations,we discuss a number of scenarios including the thermal instability of thin disks scenario,the spiral jet scenario and the binary supermassive black hole scenario.We expect that the binary black hole scenario,where the QPO is caused by the precession of binary black holes,is the most likely explanation.FSRQ PKS J0805-0111 thus could be a good binary black hole candidate.In the binary black hole scenario,the distance between the primary black hole and the secondary black hole is about 1.71×10¹⁶ cm.     

A star-based method for precise wavelength calibration of the Chinese Space Station Telescope(CSST) slitless spectroscopic survey

The Chinese Space Station Telescope(CSST)spectroscopic survey aims to deliver high-quality low-resolution(R>200)slitless spectra for hundreds of millions of targets down to a limiting magnitude of about 21 mag,distributed within a large survey area(17500 deg2)and covering a wide wavelength range(255-1000 nm by three bands GU,GV,and GI).As slitless spectroscopy precludes the usage of wavelength calibration lamps,wavelength calibration is one of the most challenging issues in the reduction of slitless spectra,yet it plays a key role in measuring precise radial velocities of stars and redshifts of galaxies.In this work,we propose a star-based method that can monitor and correct for possible errors in the CSST wavelength calibration using normal scientific observations,taking advantage of the facts that(ⅰ)there are about ten million stars with reliable radial velocities now available thanks to spectroscopic surveys like LAMOST,(ⅱ)the large field of view of CSST enables efficient observations of such stars in a short period of time,and(ⅲ)radial velocities of such stars can be reliably measured using only a narrow segment of CSST spectra.We demonstrate that it is possible to achieve a wavelength calibration precision of a few km s^-1 for the GU band,and about 10 to 20 kms^-1 for the GV and GI bands,with only a few hundred velocity standard stars.Implementations of the method to other surveys are also discussed.     

硬X射线成像仪量能器测试系统的设计与实现

硬X射线成像仪(Hard X-ray Imager, HXI)是先进天基太阳天文台(Advanced Space-based Solar Observatory, ASO-S)的3大载荷之一, 其中量能器作为其重要组成部分, 承担着观测30--200keV能段的太阳硬X射线的任务. 在卫星发射之前, 需要开展大量的测试工作, 以确保HXI量能器的各项功能和性能满足设计需求. HXI量能器通道数众多, 内含99个溴化镧探测器, 分别由8块相同的前端电子学板控制. 除了对各个通道的性能进行测试外, 地检系统还需模拟量能器在轨面对不同太阳活动时的运行情况, 对量能器进行全面完备的测试. 此外, 地检系统还需足够稳定, 能满足量能器在单机测试、环境试验、热真空与振动等多个不同测试项目的长时间测试需求. 为此, 设计了地检板与上位机软件, 结合放射源、直流电源、高压模块等组成一套HXI量能器的地检系统, 对8块前端电子学板实现同步配置与管理, 能高效完成指令发送与数据接收, 满足量能器最大数据输出带宽400Mbps的需求. 利用该系统, 在地面完成了HXI量能器的功能、性能验证, 获得了量能器的线性、死时间、能量分辨率等各项性能指标, 为HXI量能器的在轨高性能运行提供了保障.     

明安图射电频谱日像仪图像的强度定标方法

明安图射电频谱日像仪(Mingantu Spectral Radioheliograph, MUSER)能够在0.4--15GHz超宽频带内实现高时间、高空间、高频率分辨率的太阳射电成像. 而射电亮温度是描述太阳物理过程的一个重要的参数, 在研究不同射电辐射机制、太阳磁场以及太阳爆发过程中非热粒子加速等问题上有着非常重要的作用, 因此必须对MUSER观测的图像进行亮温度定标. 将介绍一种适用于射电日像仪图像强度定标的方法. 太阳射电图像中包含着太阳圆盘的结构信息, 利用射电日像仪短基线的可视度函数拟合第一类贝塞尔函数, 可以得到图像中宁静太阳圆盘的射电半径和强度, 再利用瑞利-金斯定律和每天的太阳射电流量可以计算得到每天图像的定标因子$G_c$, 从而实现对MUSER图像强度的定标. 将该方法应用到MUSER的实际观测数据中, 包括宁静太阳和有太阳射电爆发等不同的情况, $G_c$的误差基本不超过10%, 得到的宁静太阳亮温度与其他宁静太阳的结果具有较高的相关性, 表明了此方法的可行性和有效性.     

慧眼-HXMT观测到的黑洞双星系统中的W-K 关系

通过几十年的观测研究, 黑洞X射线双星(X-Ray Binary, XRB)部分特征被揭示. 然而, 吸积盘结构尚不确定. 黑洞XRB功率密度谱的截断频率与准周期振荡(Quasi Periodic Oscillation, QPO)的相关性质(W-K关系)可以限制吸积盘结构. 利用慧眼-HXMT (Hard X-ray Modulation Telescope)观测到的5个黑洞XRB的数据, 对黑洞XRB的W-K关系进行了研究, 结果表明在慧眼-HXMT观测的3个探测器能段中W-K关系成立. 此外在MAXI J1535-571之中存在截断频率和吸积盘内半径的相关关系, 这和截断的吸积盘结构一致. 如果观测到的功率密度谱来自质量吸积率的扰动传播, 可以推测吸积盘内半径接近最内圆形稳定轨道, 此黑洞可能是高自旋系统.     

利用奇异谱分析提高极移中长期预报精度

由于空间大地观测数据传输耗时及处理过程复杂, 导致极移测量值的获取存在时延, 无法满足对高精度的极移预报值有重大需求的应用领域. 针对极移复杂的时变特性, 提出一种基于奇异谱分析(singular spectrum analysis, SSA)的预报方法. 首先用SSA分离提取极移时序中的高频组分与低频组分; 其次建立最小二乘(least square, LS)外推与自回归(AutoreGressive, AR)模型对极移高频和低频组分进行组合预报. 结果表明, SSA方法能够准确地分离和提取极移低频和高频组分, 对低频和高频组分组合预报可以显著改善极移的中长期(30--365d)预报精度, 与国际地球自转和参考系服务局(International Earth Rotation and Reference Systems Service, IERS)提供的A公报中的极移预报值相比, SSA方法对极移PMX分量(本初子午线方向)和PMY分量(西90$^\circ$子午线方向)的中长期预报精度改进最高分别可达45.97%和62.44%. 研究结果验证了SSA方法对极移中长期预报改进的有效性.     

大型射电望远镜主反射面调整方法研究

随着毫米波天文学和空间通信的重要性日益提高, 对天线性能提出了越来越高的要求, 而天线性能往往受到其反射器表面精度的限制. 微波全息技术是一种快速有效的检测反射面天线表面轮廓的测量技术. 通过微波全息测量得到天线口径场, 计算天马65m射电望远镜反射面与理想抛物面的偏差. 天马65m射电望远镜的主反射面板是放射状的, 有14圈. 面板的每个角都固定在面板下方促动器的螺栓上进行上下移动, 且相邻面板交点处的拐角共用一个促动器. 采用平面拟合的方法可以计算各块面板拐角处的调整值, 但是同一个促动器会得到4个不同的调整量. 通过平面拟合, 同时以天线照明函数为权重的平差计算方法得到相邻面板拐角的一个平差值, 即天马65m射电望远镜1104个促动器的最佳调整值. 通过多次调整和新算法的应用, 天马65m射电望远镜反射面的面形精度逐渐提高到了0.24mm.     

高层大气密度的午夜极大值特征分析与建模

利用CHAMP (CHAllenging Minisatellite Payload)、GRACE-A (Gravity Recovery and Climate Experiment-A)、SWARM-C (The Earth''s Magnetic Field and Environment Explorers-C)等3颗极轨卫星的资料, 研究360—480km高层大气密度在低纬度区域的午夜极大值(Midnight Density Maximum, MDM)现象. MDM一般出现在23:00- 02:00 LT (Local Time)之间,峰值位置在低纬度15°以内,谷值位置在中纬度35°-45°附近,整体略偏向南半球,振幅约为平均密度的26%.随着高度增大以及太阳辐射水平的增强,MDM振幅呈减小趋势;冬至和夏至日附近的季节效应会减弱MDM振幅,在春秋分日的振幅最大.用3个主流大气模型DTM2000 (Drag Temperature Model 2000), NRLMSISE00 (US Naval Research Laboratory, Mass Spectrom-eter and Incoherent Scatter radar)和JB2008 (Jacchia- Bowman 2008 model)对MDM进行模拟,JB2008没有刻画出MDM现象;另两个模型低估了MDM效应,在360km和480km两个高度DTM2000模型的振幅仅为观测的46%和53%, NRLMSISE00模型仅为观测的33%和26%;模型没有准确刻画出MDM与高度、辐射水平和季节的关系.联合3颗卫星的资料,研究了-种基于地理纬度的6阶勒让德多项式,同时融合地方时和高度因素的经验函数,在振幅和相位上可以较好地刻画MDM特征,相关系数达到0.923,可为大气密度模型的修正提供借鉴,服务于低轨道航天器高精度轨道预报.