3DT-CM: A Low-complexity Cross-matching Algorithm for Large Astronomical Catalogues Using 3d-tree Approach

Location-based cross-matching is a preprocessing step in astronomy that aims to identify records belonging to the same celestial body based on the angular distance formula. The traditional approach involves comparing each record in one catalog with every record in the other catalog, resulting in a one-to-one comparison with high computational complexity. To reduce the computational time, index partitioning methods are used to divide the sky into regions and perform local cross-matching. In addition, cross-matching algorithms have been adopted on highperformance architectures to improve their efficiency. But the index partitioning methods and computation architectures only increase the degree of parallelism, and cannot decrease the complexity of pairwise-based crossmatching algorithm itself. A better algorithm is needed to further improve the performance of cross-matching algorithm. In this paper, we propose a 3d-tree-based cross-matching algorithm that converts the angular distance formula into an equivalent 3dEuclidean distance and uses 3d-tree method to reduce the overall computational complexity and to avoid boundary issues. Furthermore, we demonstrate the superiority of the 3d-tree approach over the 2d-tree method and implement it using a multi-threading technique during both the construction and querying phases. We have experimentally evaluated the proposed 3d-tree cross-matching algorithm using publicly available catalog data. The results show that our algorithm applied on two 32-core CPUs achieves equivalent performance than previous experiments conducted on a six-node CPU-GPU cluster.

NUV Star Catalog from the Lunar-based Ultraviolet Telescope Survey: First Release

We present a star catalog extracted from the Lunar-based Ultraviolet Telescope （LUT） survey program. LUT＇s observable sky area is a circular belt around the Moon＇s north pole, and the survey program covers a preferred area of about 2400 deg2 which includes a region of the Galactic plane. The data are processed with an automatic pipeline which copes with stray light contamination, artificial sources, cosmic rays, flat field calibration, photometry and so on. In the first release version, the catalog provides high confidence sources which have been cross-identified with the Tycho-2 catalog. All the sources have signalto-noise ratio larger than 5, and the corresponding magnitude limit is typically 14.4 mag, but can be as deep as -16 mag if stray light contamination is at the lowest level. A total of 86 467 stars are recorded in the catalog. The full catalog in electronic form is available online.

A statistical classification of the unassociated gamma-ray sources in the second Fermi Large Area Telescope Catalog

With the assistance of the identified/associated sources in the second Fermi Large Area Telescope （LAT） catalog, we analyze and resolve the spatial distribution and the distributions of the gamma-ray spectral and variability indices of the remaining 575 unassociated Fermi LAT sources. Consequently, it is suggested that the unassociated sources could statistically consist of Galactic supernova remnants/pulsar wind nebulae, BL Lacertae objects, fiat spectrum radio quasars and other types of active galaxies with fractions of 25%, 29%, 41% and 5%, respectively.

Correcting the Contamination of Second-order Spectra:Improving HαMeasurements in Reverberation Mapping Campaigns

Long-term spectroscopic monitoring campaigns on active galactic nuclei(AGNs)provide a wealth of information about its interior structure and kinematics.However,a number of the observations suffer from the contamination of second-order spectra(SOS)which will introduce some undesirable uncertainties at the red side of the spectra.In this paper,we test the effect of SOS and propose a method to correct it in the time domain spectroscopic data using the simultaneously observed comparison stars.Based on the reverberation mapping(RM)data of NGC 5548 in2019,one of the most intensively monitored AGNs by the Lijiang 2.4 m telescope,we find that the scientific object,comparison star,and spectrophotometric standard star can jointly introduce up to~30%SOS for Grism 14.This irregular but smooth SOS significantly affects the flux density and profile of the emission line,while having little effect on the light curve.After applying our method to each spectrum,we find that the SOS can be corrected effectively.The deviation between corrected and intrinsic spectra is~2%,and the impact of SOS on time lag is very minor.This method makes it possible to obtain the HαRM measurements from archival data provided that the spectral shape of the AGN under investigation does not have a large change.

On the Relation Between Symmetry of Radio Galaxies and Their Physical Parameters

Gravity as a fundamental force plays a dominant role in the formation and evolution of cosmic objects and leaves its effect in the emergence of symmetric and asymmetric structures.Thus,analyzing the symmetry criteria allows us to uncover mechanisms behind the gravity interaction and understand the underlying physical processes that contribute to the formation of large-scale structures such as galaxies.We use a segmentation process using intensity thresholding and the k-means clustering algorithm to analyze radio galaxy images.We employ a symmetry criterion and explore the relation between morphological symmetry in radio maps and host galaxy properties.Optical properties(stellar mass,black hole mass,optical size(R_(50)),concentration,stellar mass surface density(μ_(50)),and stellar age)and radio properties(radio flux density,radio luminosity,and radio size)are considered.We found that there is a correlation between symmetry and radio size,indicating larger radio sources have smaller symmetry indices.Therefore,size of radio sources should be considered in any investigation of symmetry.Weak correlations are also observed with other properties,such as R_(50)for FRI galaxies and stellar age.We compare the symmetry differences between FRI and FRII radio galaxies.FRII galaxies show higher symmetry in 1.4 GHz and 150 MHz maps.Investigating the influence of radio source sizes,we discovered that this result is independent of the sizes of radio sources.These findings contribute to our understanding of the morphological properties and analyses of radio galaxies.

The Fifth Edition of the General Catalogue of Variable Stars:experiences in the constellation Centaurus

We have recently announced that the General Catalogue of Variable Stars enters the stage of its fifth, purely electronic edition （GCVS 5.1）. We have included 1408 variable stars in the constellation Centaurus in this new version, GCVS 5.1. Working on this revision, we applied current possibilities from data mining, suggested new variability types for many variable stars and found new light elements for a large number of periodic variables. This paper describes the work completed during the preparation of GCVS 5.1 for Centaurus and discusses in detail a number of the most astrophysically significant cases.

The RFI Fast Mitigation Algorithm Based on Block LMS Filter

The radio telescope possesses high sensitivity and strong signal collection capabilities.While receiving celestial radiation signals,it also captures Radio Frequency Interferences(RFIs)introduced by human activities.RFI,as signals originating from sources other than the astronomical targets,significantly impacts the quality of astronomical data.This paper presents an RFI fast mitigation algorithm based on block Least Mean Square(LMS)algorithm.It enhances the traditional adaptive LMS filter by grouping L adjacent time-sampled points into one block and applying the same filter coefficients for filtering within each block.This transformation reduces multiplication calculations and enhances algorithm efficiency by leveraging the time-domain convolution theorem.The algorithm is tested using baseband data from the Parkes 64 m radio telescope's pulsar observations and simulated data.The results confirm the algorithm's effectiveness,as the pulsar profile after RFI mitigation closely matches the original pulsar profile.

The Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST)

The Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST, also called the Guo Shou Jing Telescope) is a special reflecting Schmidt telescope. LAMOST’s special design allows both a large aperture (effective aperture of 3.6 m–4.9 m) and a wide field of view (FOV) (5°). It has an innovative active reflecting Schmidt configuration which continuously changes the mirror’s surface that adjusts during the observation process and combines thin deformable mirror active optics with segmented active optics. Its primary mirror (6.67m×6.05 m) and active Schmidt mirror (5.74m×4.40 m) are both segmented, and composed of 37 and 24 hexagonal sub-mirrors respectively. By using a parallel controllable fiber positioning technique, the focal surface of 1.75 m in diameter can accommodate 4000 optical fibers. Also, LAMOST has 16 spectrographs with 32 CCD cameras. LAMOST will be the telescope with the highest rate of spectral acquisition. As a national large scientific project, the LAMOST project was formally proposed in 1996, and approved by the Chinese government in 1997. The construction started in 2001, was completed in 2008 and passed the official acceptance in June 2009. The LAMOST pilot survey was started in October 2011 and the spectroscopic survey will launch in September 2012. Up to now, LAMOST has released more than 480 000 spectra of objects. LAMOST will make an important contribution to the study of the large-scale structure of the Universe, structure and evolution of the Galaxy, and cross-identification of multiwaveband properties in celestial objects.

Optical fiber positioning based on four-quadrant detector with Gaussian fitting method

With the development of large-scale spectral surveys, fiber positioning technology has been developing rapidly. Because of the performance advantages of a four-quadrant（4Q） detector, a fiber positioning and real-time monitoring system based on the 4Q detector is proposed. The detection accuracy of this system is directly determined by the precision of the center of the spot. A Gaussian fitting algorithm based on the 4Q detector is studied and applied in the fiber positioning process to improve the calculated accuracy of the spot center. The relationship between the center position of the incident spot and the detector output signal is deduced. An experimental platform is built to complete the simulated experiment. Then we use the Gaussian fitting method to process experimental data, compare the fitting value with the theoretical one and calculate the corresponding error.

Brief Introduction to the "Catalogue of Recent Chinese Earthquakes (1912～1990, M_S≥4.7)

This paper briefly introduces the principles and methods for compiling the "Catalog of Recent Chinese Earthquakes (1912 ～ 1990, MS≥4.7)." This new catalog is compiled by revising and supplementing the 1983 version of the Catalog of Chinese Earthquakes.

An Approach of Tropospheric Correction for VLBI Phase-Referencing using GPS Data

The dominant source of error in VLBI phase-referencing is the troposphere at observing frequencies above 5 GHz. We compare the tropospheric zenith delays derived from VLBI and GPS data at VLBA stations collocated with GPS antennas. The systematic biases and standard deviations both are at the level of sub-centimeter. Based on this agreement, we suggest a new method of tropospheric correction in phase-referencing using combined VLBI and GPS data.

LAMOST CCD camera-control system based on RTS2

The Large Sky Area Multi-Object Fiber Spectroscopic Telescope（LAMOST） is the largest existing spectroscopic survey telescope, having 32 scientific charge-coupled-device（CCD） cameras for acquiring spectra. Stability and automation of the camera-control software are essential, but cannot be provided by the existing system. The Remote Telescope System 2 nd Version（RTS2） is an open-source and automatic observatory-control system. However, all previous RTS2 applications were developed for small telescopes. This paper focuses on implementation of an RTS2-based camera-control system for the 32 CCDs of LAMOST. A virtual camera module inherited from the RTS2 camera module is built as a device component working on the RTS2 framework. To improve the controllability and robustness, a virtualized layer is designed using the master-slave software paradigm, and the virtual camera module is mapped to the 32 real cameras of LAMOST. The new system is deployed in the actual environment and experimentally tested. Finally, multiple observations are conducted using this new RTS2-frameworkbased control system. The new camera-control system is found to satisfy the requirements for automatic camera control in LAMOST. This is the first time that RTS2 has been applied to a large telescope, and provides a referential solution for full RTS2 introduction to the LAMOST observatory control system.

Reconstructing the landing trajectory of the CE-3 lunar probe by using images from the landing camera

An accurate determination of the landing trajectory of Chang＇e-3 （CE-3） is significant for verifying orbital control strategy, optimizing orbital planning, accu- rately determining the landing site of CE-3 and analyzing the geological background of the landing site. Due to complexities involved in the landing process, there are some differences between the planned trajectory and the actual trajectory of CE-3. The land- ing camera on CE-3 recorded a sequence of the landing process with a frequency of 10 frames per second. These images recorded by the landing camera and high-resolution images of the lunar surface are utilized to calculate the position of the probe, so as to reconstruct its precise trajectory. This paper proposes using the method of trajectory reconstruction by Single Image Space Resection to make a detailed study of the hov- ering stage at a height of 100 m above the lunar surface. Analysis of the data shows that the closer CE-3 came to the lunar surface, the higher the spatial resolution of im- ages that were acquired became, and the more accurately the horizontal and vertical position of CE-3 could be determined. The horizontal and vertical accuracies were 7.09 m and 4.27 m respectively during the hovering stage at a height of 100.02 m. The reconstructed trajectory can reflect the change in CE-3＇s position during the powered descent process. A slight movement in CE-3 during the hovering stage is also clearly demonstrated. These results will provide a basis for analysis of orbit control strategy, and it will be conducive to adjustment and optimization of orbit control strategy in follow-up missions.

Automatic Normalization and Equivalent-Width Measurement of High-Resolution Stellar Spectra

It is well known that normalization, radial velocity correction and equivalent-width measurement of high-resolution stellar spectra are time-consuming work. In order to improve the efficiency we present an automatic method for these routines. The continuum is determined by fitting the ＇high points＇ in the spectrum. After continuum normalization, the program automatically searches for the position of the Ha line and obtains a rough radial velocity, then computes an accurate radial velocity by cross-correlation between the given spectrum and the solar spectrum. In this method, the equivalent-width is automatically measured using Gaussian fitting. A comparison between our results and those from traditional analysis shows that the typical error for equivalent width is around 3.8% in our method. Developing such automatic routines does not mean to replace the interactive reduction method： it is just for a quick extraction of information from the spectra, especially those obtained in large sky surveys.

A planetary companion orbiting the intermediate-mass G Giant HD 173416

The discovery of a planetary companion to the intermediate-mass late-type giant star HD173416 from precise Doppler surveys of G and K giants at Xinglong station and Okayama Astrophysical Observatory （OAO） is presented in this letter. The planet has a minimum mass of 2.7 MJ, an eccentricity of 0.21, a semimajor axis of 1.16 AU and an orbital period of 324 days.

The HI gas and star formation in star-forming galaxies selected from ALFALFA

This paper presents the results of Hαimaging of 169 galaxies randomly selected from the α.40-SDSS catalog.The sample has excluded all low surface brightness galaxies (LSBGs) whose central surface brightness in B band (μ0(B)) fainter than 22.5 mag arcsec-2.It can be used as the counterparts sample to LSBGs.We observed their Hα and R band images by using the 2.16 m telescope at the Xinglong Observatory of the National Astronomical Observatories,Chinese Academy of Sciences (NAOC).The main goal of this work is to present the properties of those galaxies,together with Hαflux and star formation-,gas-,stellar mass-surface density.In addition,we confirm the correlations among HI content,stellar mass and star formation in ALFALFA galaxies.The HI mass increases with stellar mass,and the slope slows down at the higher stellar mass.The overall trend was that the specific star formation rate (s SFR) decreases with stellar mass,and the sSFR dropped sharply when the stellar mass is close to 1010.3~1010.5M⊙.The weak correlation between SFR/MHIand MHIimplies the HI contribute little to star formation.Our sample,which are mostly star-forming galaxies,follows the revisited Kennicutt-Schmidt law as well as the Kennicutt-Schmidt law.

Discovery of 2716 hot emission-line stars from LAMOST DR5

We present a catalog of 3339 hot emission-line stars(ELSs)identified from 451695 O,B and A type spectra,provided by LAMOST Data Release 5(DR5).We developed an automated Python routine that identified 5437 spectra having a peak between 6561 and 6568.False detections and bad spectra were removed,leaving 4138 good emission-line spectra of 3339 unique ELSs.We re-estimated the spectral types of 3307 spectra as the LAMOST Stellar Parameter Pipeline(LASP)did not provide accurate spectral types for these emission-line spectra.As Herbig Ae/Be stars exhibit higher excess in near-infrared and mid-infrared wavelengths than classical Ae/Be stars,we relied on 2 MASS and WISE photometry to distinguish them.Finally,we report 1089 classical Be,233 classical Ae and 56 Herbig Ae/Be stars identified from LAMOST DR5.In addition,928 B[em]/A[em]stars and 240 CAe/CBe potential candidates are identified.From our sample of 3339 hot ELSs,2716 ELSs identified in this work do not have any record in the SIMBAD database and they can be considered as new detections.Identification of such a large homogeneous set of emission-line spectra will help the community study the emission phenomenon in detail without worrying about the inherent biases when compiling from various sources.

Manifestations of bright points observed in G-band and CaⅡH by Hinode/SOT

An algorithm was developed for identifying and tracking a magnetic bright point, or bright point （BP） for short, observed in both the photosphere （G-band） and chromosphere （Ca II H）, as well as for pairing a photospheric BP （PBP） with its conjugate chromospheric BP （CBP）. Two sets of data observed by Hinode/SOT in the quiet Sun near the disk center were analyzed. About 278 PBP-CBP pairs were identified and tracked. Lifetimes of both the PBPs and CBPs follow an exponential distribution with average lifetimes of 174 s and 163 s, respectively. We found that the differences in appearance time, in disappearance time and in lifetime of the two kinds of BPs all follow Gaussian distributions,which may indicate that the mechanisms of PBP and CBP formation/disintegration are different. However, the lifetimes of PBPs and CBPs are positively correlated with one another, with a correlation coefficient of 0.8. Furthermore, we calculated the horizontal displacement between the PBP and its conjugate CBP, which follows a Gaussian function with an average and standard deviation of （67.7 ± 38.5）km. We also calculated the amplitude of the flux tube shape change which might be caused by MHD waves propagating along the flux tube, and found that it follows an exponential distribution very well.

Analysis of Numerical Results in High Temperature Congealment and Chemistry Non-equilibrium Flow Field

Using the air plasma ignition technique, physicochemical process of burning can be accelerated; concentration limit ofretrofires both can be extended; reliability of retrofires and stability of burning can be improved. In this paper, using internalequivalent heat area in place of electric are that created Ohm heat, the flow fields of thermodynamic equilibrium chemistry con-gealment and chemistry non-equilibrium in the plasma generator were simulated. The influences of the inlet prerotation angleof air, the inlet total pressure of air and the airflow compression angle of spray nozzle on the temperature on the surface of elec-

Estimating Photometric Redshifts with Artificial Neural Networks and Multi-Parameters

We calculate photometric redshifts from the Sloan Digital Sky Survey Data Release 2 （SDSS DR2） Galaxy Sample using artificial neural networks （ANNs）. Different input sets based on various parameters （e.g. magnitude, color index, flux information） are explored. Mainly, parameters from broadband photometry are utilized and their performances in redshift prediction are compared. While any parameter may be easily incorporated in the input, our results indicate that using the dereddened magnitudes often produces more accurate photometric redshifts than using the Petrosian magnitudes or model magnitudes as input, but the model magnitudes are superior to the Petrosian magnitudes. Also, better performance resuits when more effective parameters are used in the training set. The method is tested on a sample of 79 346 galaxies from the SDSS DR2. When using 19 parameters based on the dereddened magnitudes, the rms error in redshift estimation is σz = 0.020184. The ANN is highly competitive tool compared to the traditional template-fitting methods when a large and representative training set is available.