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.

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.

A search for evidence of small-scale inhomogeneities in dense cores from line profile analysis

In order to search for intensity fluctuations on the HCN（1-0） and HCO＋（1-0） line pro- files, which could arise due to possible small-scale inhomogeneous structure, long-term observations of high-mass star-forming cores S140 and S199 were carried out. The data were processed by the Fourier filtering method. Line temperature fluctuations that exceed the noise level were detected. Assuming the cores consist of a large number of randomly moving small thermal fragments, the total number of frag- ments is - 4 × 106 for the region with linear size - 0.1 pc in S140 and - 106 for the region with linear size - 0.3 pc in S 199. Physical parameters of fragments in S 140 were obtained from detailed modeling of the HCN emission in the framework of the clumpy cloud model.

Example-based super-resolution for single-image analysis from the Chang'e-1 Mission

Due to the low spatial resolution of images taken from the Chang＇e-1 （CE-I） orbiter, the details of the lunar surface are blurred and lost. Considering the limited spatial resolution of image data obtained by a CCD camera on CE-1, an example-based super-resolution （SR） algorithm is employed to obtain high- resolution （HR） images. SR reconstruction is important for the application of image data to increase the resolution of images. In this article, a novel example-based algorithm is proposed to implement SR reconstruction by single-image analysis, and the computational cost is reduced compared to other example-based SR methods. The results show that this method can enhance the resolution of images using SR and recover detailed information about the lunar surface. Thus it can be used for surveying HR terrain and geological features. Moreover, the algorithm is significant for the HR processing of remotely sensed images obtained by other imaging systems.

Time-resolved spectral analysis of prompt emission from long gamma-ray bursts with GeV emission

We performed detailed time-resolved spectroscopy of bright tong gamma- ray bursts （GRBs） which show significant GeV emissions （GRB 080916C, GRB 090902B and GRB 090926A）. In addition to the standard Band model, we also use a model consisting of a black body and a power law to fit the spectra. We find that for the latter model there are indications of an additional soft component in the spectra. While previous studies have shown that such models are required for GRB 090902B, here we find that a composite spectral model consisting of two blackbodies and a power law adequately fits the data of all the three bright GRBs. We investigate the evolution of the spectral parameters and find several interesting features that appear in all three GRBs, like （a） temperatures of the blackbodies are strongly correlated with each other, （b） fluxes in the black body components are strongly correlated with each other, （c） the temperatures of the black body trace the profile of the individual pulses of the GRBs, and （d） the characteristics of power law components like the spectral index and the delayed onset bear a close similarity to the emission characteristics in the GeV regions. We discuss the implications of these results and the possibility of identifying the radiation mechanisms during the prompt emission of GRBs.

X-Ray Properties of PSR J1811-1925 by Nu STAR

We analyzed the spectral properties and pulse profile of PSR J1811-1925,a pulsar located in the center of composite supernova remnant(SNR)G11.2-0.3,by using high timing resolution archival data from the Nuclear Spectroscopic Telescope Array Mission(NuSTAR).Analysis of archival Chandra data over different regions rules out the SNR shell as the site of the hard X-ray emission while spectral analysis indicates that the NuSTAR photons originate in the pulsar and its nebula.The pulse profile exhibits a broad single peak up to 35 keV.The jointed spectrum by combining NuSTAR and Chandra can be well fitted by a power-law model with a photon index ofΓ=1.58±0.04.The integrated flux of jointed spectrum over 1-10 keV is 3.36×10^(-12)erg cm^(-2)s^(-1).The spectrum of pulsar having photon indexΓ=1.33±0.06 and a 1-10 keV flux of 0.91×10^(-12)erg cm^(-2)s^(-1).We also performed the phase-resolved spectral analysis by splitting the whole pulse-on phase into five phase bins.The photon indices of the bins are all around 1.4,indicating that the photon index does not evolve with the phase.

An algorithm to resolve γ-rays from charged cosmic rays with DAMPE

The DArk Matter Particle Explorer（DAMPE）,also known as Wukong in China,which was launched on 2015 December 17,is a new high energy cosmic ray and γ-ray satellite-borne observatory.One of the main scientific goals of DAMPE is to observe Ge V-Te V high energy γ-rays with accurate energy,angular and time resolution,to indirectly search for dark matter particles and for the study of high energy astrophysics. Due to the comparatively higher fluxes of charged cosmic rays with respect to γ-rays,it is challenging to identify γ-rays with sufficiently high efficiency,minimizing the amount of charged cosmic ray contamination. In this work we present a method to identify γ-rays in DAMPE data based on Monte Carlo simulations,using the powerful electromagnetic/hadronic shower discrimination provided by the calorimeter and the veto detection of charged particles provided by the plastic scintillation detector. Monte Carlo simulations show that after this selection the number of electrons and protons that contaminate the selected γ-ray events at～10 Ge V amounts to less than 1% of the selected sample.Finally,we use flight data to verify the effectiveness of the method by highlighting known γ-ray sources in the sky and by reconstructing preliminary light curves of the Geminga pulsar.

Wavelet Analysis of Several Important Periodic Properties in the Relative Sunspot Numbers

We investigate the wavelet transform of yearly mean relative sunspot number series from 1700 to 2002. The curve of the global wavelet power spectrum peaks at 11-yr, 53-yr and 101-yr periods. The evolution of the amplitudes of the three periods is studied. The results show that around 1750 and 1800, the amplitude of the 53-yr period was much higher than that of the the 11-yr period, that the ca. 53-yr period was apparent only for the interval from 1725 to 1850, and was very low after 1850, that around 1750, 1800 and 1900, the amplitude of the 101-yr period was higher than that of the 11-yr period and that, from 1940 to 2000, the 11-yr period greatly dominates over the other two periods.

The W UMa binaries USNO-A2.0 1350-17365531,V471 Cas,V479 Lac and V560 Lac:light curve solutions and global parameters based on Gaia distances

We present photometric observations in Sloan filters g′, i′of the eclipsing W UMa stars USNOA2.0 1350-17365531, V471 Cas, V479 Lac and V560 Lac. The sinusoidal-like O-C diagram of V471 Cas indicates the presence of a third body with mass 0.12 M_⊙(a red dwarf) at distance 897 R_⊙. The O-C diagram of V479 Lac reveals a period decrease of d P/dt =-1.69 × 10-6d yr-1. The results of the light curve solutions are:(i) the targets are overcontact binaries with small fill-out factors;(ii) their components are F–K stars, comparable in size, whose temperature differences are below 80 K;(iii) all targets undergo partial eclipses and to limit the possible mass ratios we carried out two-step q-search analysis. The target global parameters(luminosities, radii, masses) were obtained on the basis of their Gaia distances and the results of our light curve solutions. The obtained total mass of V560 Lac turns out to be smaller than the lower mass limit for presently known W UMa binaries of 1.0-1.2 M_⊙, i.e. this target is a peculiar overcontact system.

Relative Spectral Lag: a New Redshift Indicator of Gamma-ray Bursts

Using 64 ms count data of long gamma-ray bursts （T90 〉 2.6 s）, we analyze the quantity named relative spectral lag （RSL）, T31/FWHM（1）. We investigated in detail all the correlations between the RSL and other parameters for a sample of nine long bursts, using the general cross-correlation technique that includes the lag between two different energy bands. We conclude that the distribution of RSLs is normal and has a mean value of 0.1; that the RSLs are weakly correlated with the FWHM, the asymmetry, peak flux （Fp）, peak energy （Ep） and spectral indexes （α and β）, while they are uncorrelated with τ31, the hardness- ratio （HR31） and the peak time （tm）. Our important discovery is that redshift （z） and peak luminosity （Lp） are strongly correlated with the RSL, which can be measured easily and directly, making the RSL a good redshift and peak luminosity indicator.

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.

An unusual timing and spectral state of a black hole microquasar XTE J1550-564

An unusual timing and spectral state of a black hole microquasar XTE J1550- 564 observed with RXTE is analyzed. Millisecond variabilities are found, which are significantly shorter than the minimum possible time scale in the light curves of black hole binaries, as suggested by Sunyaev ＆ Revnivtsev （2000）. The X-ray spectral fitting result indicates that there is an unusual soft component in the spectrum, which may be responsible for the millisecond variabilities. The millisecond variabilities as well as the unusual soft spectral component should be produced from some small, but independent active regions in the accretion disk.

Timing in the Time Domain: Cygnus X-1

Quantities characterizing temporal property, e.g., power density, coherence, and time lag, can be defined and calculated directly in the time domain without using the Fourier transformation. Spectral hardness, variability duration, and correlation between different characteristic quamtities on different time scale can be studied in the time domain as well. The temporal analysis technique in the time domain is a powerful tool, particularly in studying rapid variability on short time scales (or in high frequencies). Results of studying variabilities of X-rays from Cyg X-1 with the analysis technique in the time domain and RXTE data reveal valuable clues to understanding production and propagation processes of X-rays and structure of accretion disk in the black hole system.

Automated Separation of Stars and Normal Galaxies Based on Statistical Mixture Modeling with RBF Neural Networks

For LAMOST, the largest sky survey program in China, the solution of the problem of automatic discrimination of stars from galaxies by spectra has shown that the results of the PSF test can be significantly refined. However, the problem is made worse when the redshifts of galaxies are not available. We present a new automatic method of star/(normal) galaxy separation, which is based on Statistical Mixture Modeling with Radial Basis Function Neural Networks (SMM-RBFNN). This work is a continuation of our previous one, where active and non-active celestial objects were successfully segregated. By combining the method in this paper and the previous one, stars can now be effectively separated from galaxies and AGNs by their spectra-a major goal of LAMOST, and an indispensable step in any automatic spectrum classification system. In our work, the training set includes standard stellar spectra from Jacoby's spectrum library and simulated galaxy spectra of EO, SO, Sa, Sb types with redshift ranging from 0 to 1.2, and the test set of stellar spectra from Pickles' atlas and SDSS spectra of normal galaxies with SNR of 13. Experiments show that our SMM-RBFNN is more efficient in both the training and testing stages than the BPNN (back propagation neural networks), and more importantly, it can achieve a good classification accuracy of 99.22% and 96.52%, respectively for stars and normal galaxies.

New open cluster candidates discovered in the XSTPS-GAC survey

The Xuyi Schmidt Telescope Photometric Survey of the Galactic Anti-center（XSTPS-GAC）is a photometric sky survey that covers nearly 6000 deg^2 towards the Galactic Anti-center（GAC） in the g,r,i bands. Half of its survey field is located on the Galactic Anti-center disk,which makes XSTPSGAC highly suitable to search for new open clusters in the GAC region. In this paper,we report new open cluster candidates discovered in this survey,as well as properties of these open cluster candidates,such as age,distance and reddening,derived by isochrone fitting in the color-magnitude diagram（CMD）.These open cluster candidates are stellar density peaks detected in the star density maps by applying the method from Koposov et al. Each candidate is inspected in terms of its true color image composed from three XSTPS-GAC band images. Then its CMD is checked,in order to identify whether the central region stars have a clear isochrone-like trend differing from background stars. The parameters derived from isochrone fitting for these candidates are mainly based on three band photometry of XSTPS-GAC.Moreover,when these new candidates are able to be seen clearly in 2 MASS data,their parameters are also derived based on the 2 MASS（J-H,J） CMD. There are a total of 320 known open clusters rediscovered and 24 new open cluster candidates discovered in this work. Furthermore,the parameters of these new candidates,as well as another 11 previously known open clusters,are properly determined for the first time.

Weekly inter-technique combination of SLR,VLBI,GPS and DORIS at the solution level

Constructing and maintaining a stable terrestrial reference frame （TRF） is one of the key objectives of fundamental astronomy and geodesy. The datum realization for all the global TRF versions, such as ITRF2014 and its predecessor ITRF2008, assumes linear time evolution for transformation parameters and then imposes some conditions on these Helmert transformation parameters. In this paper, we investigate a new approach, which is based on weekly estimation of station positions and Helmert transformation parameters from a combination of the solutions of four space-geodetic techniques, i.e., Satellite Laser Ranging （SLR）, Very Long Baseline Interferometry （VLBI）, Global Positioning System （GPS） and Doppler Orbitography and Radiopositioning Integrated by Satellite （DORIS）. For this study, an interval of one week is chosen because the arc length of the SLR solutions is seven days. The major advantage of this weekly estimated reference frame is that both the non-linear station motions and the non-linear origin motion are implicitly taken into account. In order to study the non-linear behavior of station motions and physical parameters, ITRF2008 is used as a reference. As for datum definition of weekly reference frame, on one hand SLR is the unique technique to realize the origin and determine the scale together with VLBI, and on the other hand the orientation is realized via no net rotation with respect to ITRF2005 on a subset of core stations. Given the fact that without enough collocations an inter-technique combined TRF could not exist, the selection and relative weight of local ties surveyed at co-location sites are critical issues. To get stable results, we first assume that, if there were no events such as equipment changes between the measurement epoch of the local tie and that of the space- geodetic solution, the relative position between the two co-located stations should be invariant and this local tie could be used for computing the inter-technique combined reference flame in those weeks during the stable period of this tie. The resulting time series of both station positions and transformation parameters are studied in detail and are compared with ITRF2008. The residual station positions in the weekly combined reference frame are usually in the range of two millimeters without any periodic characteristic, but the residual station positions, when subtracting the regularized station position in ITRF2008, may reach a magnitude of a few centimeters and seem to have a significant annual signal. The physical parameter series between the weekly reference frame and ITRF2008 also show the obvious existence of an annual signal and reach a magnitude of one centimeter for origin motion and two parts per billion （ppb） for scale.

Evaluation of Quantum Chemical Methods and Basis Sets Applied in the Molecular Modeling of Artemisinin

In this paper, we evaluate semiempirical methods (AM1, PM3, and ZINDO), HF and DFT (B3LYP) in different basis sets to determine which method best describes the sign and magnitude of the geometrical parameters of artemisinin in the region of the endoperoxide ring compared to crystallographic data. We also classify these methods using statistical analysis. The results of PCA were based on three main components, explaining 98.0539% of the total variance, for the geometrical parameters C3O13, O1O2C3, O13C12C12a, and O2C3O13C12. The DFT method (B3LYP) corresponded well with the experimental data in the hierarchical cluster analysis (HCA). The experimental and theoretical angles were analyzed by simple linear regression, and statistical parameters (correlation coefficients, significance, and predictability) were evaluated to determine the accuracy of the calculations. The statistical analysis exhibited a good correlation and high predictive power for the DFT (B3LYP) method in the 6-31G** basis set.

Autoregressive Spectral Estimation for Quasi-Periodic Oscillations

Modern methods of spectral estimation based on parametric time-series models are useful tools in power spectral analysis. We apply the autoregressive （AR） model to study quasi-periodic oscillations （QPOs）. An empirical formula to estimate the expectation and standard deviation of the noise AR power densities is derived, which can be used to estimate the statistical significance of an apparent QPO peak in an AR spectrum. An iterative adding-noise algorithm in AR spectral analysis is proposed and applied to studying QPOs in the X-ray binary Cir X-1.

R-band host galaxy contamination of TeV γ-ray blazar Mrk 501： effects of aperture size and seeing

We simulated the R-band contribution of the host galaxy of TeV γ-ray BL Lac object Mrk 501 in different aperture sizes and seeing conditions. An intensive set of observations was acquired with the 1.02m optical telescope, managed by Yunnan Observatories, from 2010 May 15 to 18. Based on the host subtraction data usually used in the literature, the subtraction of host galaxy contamination results in significant seeing-brightness correlations. These correlations would lead to illusive large amplitude variations at short timescales, which will mask the intrinsic microvariability, thus giving rise to difficulty in detecting the intrinsic microvariability. Both aperture size and seeing condition influence the flux measurements, but the aperture size impacts the result more significantly. Based on the parameters of an elliptical galaxy provided in the literature, we simulated the host contributions of Mrk 501 in different aperture sizes and seeing conditions. Our simulation data of the host galaxy obviously weaken these significant seeing-brightness correlations for the host-subtracted brightness of Mrk 501, and can help us discover the intrinsic short timescale microvariability. The pure nuclear flux is -8.0 mJy in the R band, i.e., the AGN has a magnitude of R - 13.96 mag.

Time-dependent Behaviour of the Low Amplitudeδ Scuti Star HD 52788

Wavelet transform is applied to reanalyze the low amplitude δ Scuti star HD 52788, which exhibits complex light variations with uncertain frequency solutions. We gain an insight into the strong instability of pulsation present in the star. Based on an estimate of the star's physical parameters, its evolutionary status is determined. An attempt of asteroseismic modelling failed to predict the observed dense frequencies. Because of its varying pulsation spectrum, HD 52788 is a distinctive and very interesting object among δ Sct stars for testing current models of stellar evolution and pulsation.