动力天文中常用多步方法的比较（Ⅱ）

本文从几个方面给出了动力天文中三种常用的多步积分方法的比较结果,为积分器的选择提供了依据。     

射电天文毫米波接收机强度校准

射电天文强度校准的目的是将接收设备对天文观测源的响应转换为天文意义上的流量.在常用的射电天文强度校准方法中,厘米波校准主要使用噪声注入模式,就是将1个标准噪声信号在馈源和极化器之后注入到接收机内部进行校准.由于毫米波微波器件的小型化导致噪声注入模式不易实现,加之注入模式可能引入噪声,因此BTL (Bell Telephone Laboratory)最早提出使用斩波轮技术进行毫米波校准,就是在馈源口面交替放置或者移除常温黑体进行校准;之后BIMA (Berkeley Illinois Maryland Association Array)又提出使用常温、热负载进行校准; ALMA (Atacama Large Millimeter Array)对单、双负载校准方式的精度进行计算后,认为双温度负载校准方式有潜力实现1%的校准精度,并最终设计出机械臂式双温度负载校准机构;此后, GBT (Green Bank Telescope) 4 mm波段制冷接收机设计出旋转盘式双温度校准机构; OSO (Onsala Space Observatory)最新研制的3 mm波段制冷接收机设计出波束切换式双温度校准机构.中国科学院新疆天文台QTT(Qi Tai Telescope)项目的启动推动了毫米波接收机研制进程,为提高毫米波强度校准精度,相关的技术预研已经开始.     

卫星双向定时精度分析

介绍了Compass卫星双向定时的基本原理,在讨论双向定时计算模型的基础上,给出了考虑Sagnac效应项的详细计算模型,并计算了Sagnac效应项在几个代表性地区的影响,结果表明:对于距离中心站较远地区的用户,这种影响可达几十个纳秒。最后,对Compass卫星双向定时理论精度进行了分析,并采用几个地区的试验结果进行了验证。     

天文知识讲座(4)——天球坐标与行星的视运动

一、天球坐标系 天体在天球上的视位置,可以用天体在天球的球面坐标来表示,天文常用的球面坐标系有地平坐标系、赤道坐标系和黄道坐标系。 1.地平坐标系 地平坐杯系有两个参量:地平高度和方位角。观测者的头顶方向与天球相交点叫天顶点(Z)。从观测者脚底方向与天球的交点叫天底(Z′)。与天顶和天底连线     

深度学习在天文学中的应用与改进

近年来,深度学习和人工智能技术迅猛发展,在多个学科领域得到了广泛关注和应用。天文学研究也不甘落后,涌现出一大批应用深度学习进行数据分析的工作。总结了深度学习在天文中的应用情况和趋势、天文数据类型和机器学习任务、天文中常用的深度学习网络模型和方法,以及深度学习在天文研究中的代表性应用和进展,并探讨和提出了其未来在天文学领域中的应用和改进建议。     

基于SLR精密轨道的天文定位精度分析

利用全球卫星激光测距服务系统(ILRS,International Laser Ranging Service)标准点资料对Ajisai卫星进行精密定轨,残差均方根(RMS)优于3 cm,得到该星的精密轨道.进而对长春站40 cm空间碎片光电望远镜获得的Ajisai卫星的天文定位资料进行精度分析,外符合精度约3″左右.单独利用天文定位数据进行轨道改进,内符合精度优于3″.改进轨道的x、y、z坐标3分量在观测数据覆盖范围内的精度在100 m之内.同样地对Jason-1卫星作数据分析,结果和Ajisai卫星精度相当.分析各个弧段的精度变化,发现定标星个数减少,会导致天文定位精度下降.据此提出可以把最少定标星比例作为评定数据质量的参考指标之一.     

星表系统差及其计算方法

本文对星表系统差的性质进行了讨论，并分析和评估了目前常用的星表系统差计算方法的优、缺点，从而提出了一个新的星表系统差计算方法．同时也提出了星表系统差计算精度的检验标准．最后对用不同方法计算的系统差进行了分析比较．结果表明，新方法的系统差计算精度好于目前常用的解析方法．     

大气折射的计算

把普尔科沃大气折射表的第4、第5版进行完全计算机化，用多项式拟合大气折射表，直至天顶距83°为止，公式和表列值完全符合，同时用实际计算进行了精度比较。中丹水平子午环大气折射计算采用了计算机化的算法。大气折射是精密天体定位的重要研究课题之一，它除了影响光学天文定位外，还严重地影响射电天文、激光精密定位及GPS比对精度，随着天文定位精度的提高，大气折射的研究更显出其重要性。     

适用于古天文研究的计算机软件

天文软件是进行古天文研究非常有用的辅助工具。对一些常见的商业软件与专业软件的计算结果进行了分析比较,发现它们在两三百年的短期内都是足够可靠的,但是若干商业软件并不适用于涉及几千年跨度的古天文研究。分析表明:无论从数据精确度和长期稳定性,还是从功能的强大程度来说,SkyMap都是值得推广的天文计算软件。     

马克苏托夫望远镜的一种设计方法

马克苏托夫望远镜由球面反射镜和弯月形透镜组成,弯月形透镜主要的作用是抵消球面反射镜的球差,采用一定厚度的弯月形透镜可以不引进色差,适当选择弯月形透镜到球面反射镜的距离还可以消去彗差,并且这种系统的象散也是非常小的,但场曲是存在的,焦面大致等于以焦距为半径的一个球面,只要把底片弯成和焦面一样,那末所拍的照片在5°甚至8°之内星象都是很优良的.这种望远镜是有畸变的,不过并不太大,对天文工     

Presence of LaO, ScO and VO Molecular Lines in Sunspot Umbral Spectra

High-resolution Fourier Transform Spectrometer sunspot umbral spectra of the National Solar Observatory/National Optical Astronomy Observatory at Kitt Peak were used to detect rotational lines from 19 electronic transition bands of the molecules LaO, ScO and VO, in the wavenumber range of 11 775 to 20 600 cm⁻¹. The presence of lines from the following transitions is confirmed: A ² Π _r1/2 – X ² Σ ⁺(0, 0; 0, 1), A ² Π _r3/2 – X ² Σ ⁺(1, 0), B ² Σ ⁺ – X ² Σ ⁺(0, 0; 0, 1; 1, 0) and C ² Π _r1/2 – A′²Δ _r3/2(0, 0; 1, 1) of LaO; A ² Π _r3/2 – X ² Σ ⁺(0, 0), A ² Π _r1/2 – X ² Σ ⁺(0, 0) and B ² Σ ⁺ – X ² Σ ⁺(0, 0) of ScO; and C ⁴ Σ ⁻ – X ⁴ Σ ⁻(0, 1; 1, 0; 0, 2) and (2, 0) of VO. However, the presence of A ² Π _r3/2 – X ² Σ ⁺(0, 0) and C ² Π _r3/2 – A′²Δ _r5/2(0, 0; 1, 1) of LaO and C ⁴ Σ ⁻ – X ⁴ Σ ⁻(0, 0) of VO are found to be doubtful because the lines are very weak, and detections are difficult owing to heavy blending by strong rotational lines of other molecules. Equivalent widths are measured for well-resolved lines and, thereby, the effective rotational temperatures are estimated for the systems for which the presence is confirmed.     

Sources and distribution of trace elements in Estonian peat

This paper presents the results of the distribution of trace elements in Estonian mires. Sixty four mires, representative of the different landscape units, were analyzed for the content of 16 trace elements (Cr, Mn, Ni, Cu, Zn, and Pb using AAS; Cd by GF-AAS; Hg by the cold vapour method; and V, Co, As, Sr, Mo, Th, and U by XRF) as well as other peat characteristics (peat type, degree of humification, pH and ash content). The results of the research show that concentrations of trace elements in peat are generally low: V 3.8 ± 0.6, Cr 3.1 ± 0.2, Mn 35.1 ± 2.7, Co 0.50 ± 0.05, Ni 3.7 ± 0.2, Cu 4.4 ± 0.3, Zn 10.0 ± 0.7, As 2.4 ± 0.3, Sr 21.9 ± 0.9, Mo 1.2 ± 0.2, Cd 0.12 ± 0.01, Hg 0.05 ± 0.01, Pb 3.3 ± 0.2, Th 0.47 ± 0.05, U 1.3 ± 0.2 μg g^− 1 and S 0.25 ± 0.02%. Statistical analyses on these large database showed that Co has the highest positive correlations with many elements and ash content. As, Ni, Mo, ash content and pH are also significantly correlated. The lowest abundance of most trace elements was recorded in mires fed only by precipitation (ombrotrophic), and the highest in mires fed by groundwater and springs (minerotrophic), which are situated in the flood plains of river valleys. Concentrations usually differ between the superficial, middle and bottom peat layers, but the significance decreases depending on the type of mire in the following order: transitional mires – raised bogs – fens. Differences among mire types are highest for the superficial but not significant for the basal peat layers.The use of peat with high concentrations of trace elements in agriculture, horticulture, as fuel, for water purification etc., may pose a risk for humans: via the food chain, through inhalation, drinking water etc.     

Climatic implications of correlated Upper Pleistocene glacial and fluvial deposits on the Cinca and Gállego Rivers (NE Spain) based on OSL dating and soil stratigraphy

We correlate Upper Pleistocene glacial and fluvial deposits of the Cinca and Gállego River valleys (south central Pyrenees and Ebro basin, Spain) using geomorphic position, luminescence dates, and time-related trends in soil development. The ages obtained from glacial deposits indicate glacial periods at 85 ± 5 ka, 64 ± 11 ka, and 36 ± 3 ka (from glacial till) and 20 ± 3 ka (from loess). The fluvial drainage system, fed by glaciers in the headwaters, developed extensive terrace systems in the Cinca River valley at 178 ± 21 ka, 97 ± 16 ka, 61 ± 4 ka, 47 ± 4 ka, and 11 ± 1 ka, and in the Gállego River valley at 151 ± 11 ka, 68 ± 7 ka, and 45 ± 3 ka. The times of maximum geomorphic activity related to cold phases coincide with Late Pleistocene marine isotope stages and Heinrich events. The maximum extent of glaciers during the last glacial occurred at 64 ± 11 ka, and the terraces correlated with this glacial phase are the most extensive in both the Cinca (61 ± 4 ka) and Gállego (68 ± 7 ka) valleys, indicating a strong increase in fluvial discharge and availability of sediments related to the transition to deglaciation. The global Last Glacial Maximum is scarcely represented in the south central Pyrenees owing to dominantly dry conditions at that time. Precipitation must be controlled by the position of the Iberian Peninsula with respect to the North Atlantic atmospheric circulation system. The glacial systems and the associated fluvial dynamic seem sensitive to 1) global climate changes controlled by insolation, 2) North Atlantic thermohaline circulation influenced by freshwater pulses into the North Atlantic, and 3) anomalies in atmospheric circulation in the North Atlantic controlling precipitation on the Iberian Peninsula. Our scenario of glacial and fluvial evolution during the Late Pleistocene in northern Spain could be extrapolated to other glaciated mountainous areas in southern Europe.     

Stability of motion in the Sitnikov 3-body problem

We study the stability of motion in the 3-body Sitnikov problem, with the two equal mass primaries (m ₁ = m ₂ = 0.5) rotating in the x, y plane and vary the mass of the third particle, 0 ≤ m ₃ < 10⁻³, placed initially on the z-axis. We begin by finding for the restricted problem (with m ₃ = 0) an apparently infinite sequence of stability intervals on the z-axis, whose width grows and tends to a fixed non-zero value, as we move away from z = 0. We then estimate the extent of “islands” of bounded motion in x, y, z space about these intervals and show that it also increases as |z| grows. Turning to the so-called extended Sitnikov problem, where the third particle moves only along the z-axis, we find that, as m ₃ increases, the domain of allowed motion grows significantly and chaotic regions in phase space appear through a series of saddle-node bifurcations. Finally, we concentrate on the general 3-body problem and demonstrate that, for very small masses, m ₃ ≈ 10⁻⁶, the “islands” of bounded motion about the z-axis stability intervals are larger than the ones for m ₃ = 0. Furthermore, as m ₃ increases, it is the regions of bounded motion closest to z = 0 that disappear first, while the ones further away “disperse” at larger m ₃ values, thus providing further evidence of an increasing stability of the motion away from the plane of the two primaries, as observed in the m ₃ = 0 case.     

A Study on the CO Isotopic Lines of the Star Forming Region AFGL 5157

By the mapping observations simultaneously at the ¹²CO (J=1-0), ¹³CO (J=1-0), and C¹⁸O (J=1-0) lines on the area of 24’×24’ (12 pc×12 pc) of the star forming region AFGL 5157, we have obtained the distribution and averaged physical parameters for the respective ¹³CO and C¹⁸O cores of this molecu- lar cloud. At the edge of the molecular cloud, the isotopic abundance ratio is X [(¹³CO)/(C¹⁸O)] ≈ 10, close to the ratio of a giant molecular cloud. The viral masses of the ¹³CO and C¹⁸O cores are less than the masses of the molecu-lar cloud cores, so the molecular cloud cores are gravitationally unstable, and the C¹⁸O molecular cloud core is more easy to collapse. The column density distributions of the C¹⁸O molecular cloud core in the northeast and southwest directions are, respectively, 1.1 × 10²³× z^−0.43 and 4.6 × 10²⁵× z^−0.58, where z is the distance from the center of the molecular cloud core. The high velocity molecular out?ow has been con?rmed from our ¹²CO spectra, the mass loss rate of the out?ow has been estimated, and the mass-velocity relation of the out?ow is ?tted by a power-law function of m ∝ v^−1.8. The star formation rate of the ¹³CO molecular cloud core is as high as 23%, probably, under the in?uence of     

X-ray observation of the Jovian impacts of comet Shoemaker - Levy 9

An ASCA observation of the Jovian impact of the comet Shoemaker-Levy 9 is reported. Four impacts of H, L, Q1 and R were observed and four impacts of B, C, G, and Q2 were observed within 60 minutes after their impacts. No significant flaring of X-ray emission was observed. Upper limit X-ray fluxes of 90 % confidence level, averaged 5 minutes just after the impacts, were 2.4 × 10^–13 erg sec^–1 cm^–2, 3.5 × 10^–13 erg sec^–1 cm^–2, 1.6 × 10^–13 erg sec^–1 cm^–2 and 2.9 × 10^–13 erg sec^–1 cm^–2 for the impacts of H, L, Q1 and R, respectively, in the 0.5(0.7 for H and Q1)–10 keV energy range. However, a hint of X-ray enhancement around Jupiter from July 17 to July 19 was detected with about 2 6 × 10^–14 erg sec^–1 cm^–2 in the 0.5–10 keV energy range.     

Thermonuclear reaction rates from (p,n) reactions on targets withA = 92–122

From the experimentally measured (p,n) cross-sections for^92,94Zr,⁹³Nb,^95,96,98Mo,¹⁰³Rh,^107,109Ag,¹¹⁰Pd,¹¹⁵In,^117,122Sn nuclides, for proton energies below 7 MeV, thermonuclear reaction rates in the temperature range 1≤T ₉≤5(T₉=10⁹ K) have been extracted. These reaction rates have been fitted to a three parameter empirical expression proposed by Fowler.     

Production and detection of carbon dioxide on Iapetus

Cassini VIMS detected carbon dioxide on the surface of Iapetus during its insertion orbit. We evaluated the CO₂ distribution on Iapetus and determined that it is concentrated almost exclusively on Iapetus’ dark material. VIMS spectra show a 4.27-μm feature with an absorption depth of 24%, which, if it were in the form of free ice, requires a layer 31 nm thick. Extrapolating for all dark material on Iapetus, the total observable CO₂ would be 2.3 × 10⁸ kg.Previous studies note that free CO₂ is unstable at 10 AU over geologic timescales. Carbon dioxide could, however, be stable if trapped or complexed, such as in inclusions or clathrates. While complexed CO₂ has a lower thermal volatility, loss due to photodissociation by UV radiation and gravitational escape would occur at a rate of 2.6 × 10⁷ kg year⁻¹. Thus, Iapetus’ entire inventory of surface CO₂ could be lost within a few decades.The high loss/destruction rate of CO₂ requires an active source. We conducted experiments that generated CO₂ by UV radiation of simulated icy regolith under Iapetus-like conditions. The simulated regolith was created by flash-freezing degassed water, crushing it into sub-millimeter sized particles, and then mixing it with isotopically labeled amorphous carbon (¹³C) dust. These samples were placed in a vacuum chamber and cooled to temperatures between 50 K and 160 K. The samples were irradiated with UV light, and the products were measured using a mass spectrometer, from which we measured ¹³CO₂ production at a rate of 2.0 × 10¹² mol s⁻¹. Extrapolating to Iapetus and adjusting for the solar UV intensity and Iapetus’ surface area, we calculated that CO₂ production for the entire surface would be 1.1 × 10⁷ kg year⁻¹, which is only a factor of two less than the loss rate. As such, UV photochemical generation of CO₂ is a plausible source of the detected CO₂.     

Removal of Titan's noble gases by their trapping in its haze

Titan's haze, formed by photolysis of C₂H₂, C₂H₄ and HCN, was found experimentally to trap Ar, Kr and Xe with efficiencies of 3.5 × 10⁻⁴, 1.9 × 10⁻³ and 6.5 × 10⁻² [noble gas atom]/[carbon atom] in the polymer, respectively. The rate of aerosol formation and settling down of 3 × 10⁻¹³ kg m⁻² s⁻¹, as inferred from our experiments on CH₄ photolysis in the far UV [Podolak, M., Bar-Nun, A., 1979. Icarus 39, 272-276], is sufficient to reduce the mixing ratios of ³⁶Ar and ⁴⁰Ar to their low values of (2.8 ± 0.3) × 10⁻⁷ and (4.3 ± 0.1) × 10⁻³, respectively, and those of Kr and Xe to below the detection limit of 10⁻⁸.     

HCN fluorescence on Titan

The HCN emission features near 3 μm recently detected by Geballe et al. (2003, Astrophys. J. 583, L39) are analyzed with a model for fluorescence of sunlight in the ν₃ band of HCN. The emission spectrum is consistent with current knowledge of the atmospheric temperature profile and the HCN distribution inferred from millimeter-wave observations. The spectrum is insensitive to the abundance of HCN in the thermosphere and the thousand-fold enhancement relative to photochemical models suggested by Geballe et al. (2003, Astrophys. J. 583, L39) is not required to explain the observations. We find that the spectrum can be matched with temperatures from 130 to 200 K, with slightly better fits at high temperature, contrary to the temperature determination of 130±10 K of Geballe et al. (2003, Astrophys. J. 583, L39). The HCN emission spectrum is sensitive to the collisional de-excitation probability, P₁₀, for the ν₃ state and we determine a value of 10⁻⁵ with an accuracy of about a factor of two. Analysis of absorption lines in the C₂H₂ν₃ band near 3 μm, detected in the same spectrum, indicate a C₂H₂ mole fraction near 0.01 μbar of 10⁻⁵ for P₁₀=10⁻⁴. The derived mole fraction, however, is dependent upon the value adopted for P₁₀ and lower values are required if P₁₀ at Titan temperatures is less than its room temperature value.