Nonlinear analytic solution for perturbed relative motion using differential equinoctial elements

This paper develops a nonlinear analytic solution for satellite relative motion in J₂-perturbed elliptic orbits by using the geometric method that can avoid directly solving the complex differential equations. The differential equinoctial elements (DEEs) are used to remove any singularities for zero-eccentricity or zero-inclination orbits. Based on the relationship between the relative states and the DEEs, state transition tensors (STTs) for transforming the osculating DEEs and propagating the mean DEEs have been derived. The formulation of these STTs has been split into a set of vector and matrix operations, which avoids directly expanding the complex second-order terms, and thus, the obtained STTs could be easy-to-understand and easy-to-code. Numerical results show that the proposed nonlinear solution is valid for zero-eccentricity and zero-inclination reference orbit and is more accurate than the previous linear or nonlinear methods for the long-term prediction of satellite relative motion.     

Observation and Research of the Transits of Extrasolar Planets

In order to make more accurate estimates of the physical parameters of extrasolar planets, the observation and research on the phenomena of transits of the 7 known stars with the planetary systems TrES-1, TrES-3, XO-2, WASP-1, WASP-2, WASP-3 and HAT-P-7 are carried out with the 1-m reflecting telescope at Shandong University Weihai Astronomical Observatory/Weihai Municipal Astronomical Observatory. The fundamental conditions of the observation and data processing are introduced, and the transit light curves as well as the parameters of some planets derived from them are given. At the same time, when the results are summarized and analyzed, the prospects for the more in-depth and detailed researches which will be further carried out are also described.     

Simulation of the Unusual Solar Minimum with 3D SIP-CESE MHD Model by Comparison with Multi-Satellite Observations

The observations both near the Sun and in the heliosphere during the activity minimum between solar cycles 23 and 24 exhibit different phenomena from those typical of the previous solar minima. In this paper, we have chosen Carrington rotation 2070 in 2008 to investigate the properties of the background solar wind by using the three-dimensional (3D) Solar?CInterPlanetary Conservation Element/Solution Element Magnetohydrodynamic (MHD) model. We also study the effects of polar magnetic fields on the characteristics of the solar corona and the solar wind by conducting simulations with an axisymmetric polar flux added to the observed magnetic field. The numerical results are compared with the observations from multiple satellites, such as the Solar and Heliospheric Observatory (SOHO), Ulysses, Solar Terrestrial Relations Observatory (STEREO), Wind and the Advanced Composition Explorer (ACE). The comparison demonstrates that the first simulation with the observed magnetic fields reproduces some observed peculiarities near the Sun, such as relatively small polar coronal holes, the presence of mid- and low-latitude holes, a tilted and warped current sheet, and the broad multiple streamers. The numerical results also capture the inconsistency between the locus of the minimum wind speed and the location of the heliospheric current sheet, and predict slightly slower and cooler polar streams with a relatively smaller latitudinal width, broad low-latitude intermediate-speed streams, and globally weak magnetic field and low density in the heliosphere. The second simulation with strengthened polar fields indicates that the weak polar fields in the current minimum play a crucial role in determining the states of the corona and the solar wind.     

Multi-Wavelength Variability in PKS 2155-304

We study multi-wavelength variability in BL Lacertae object PKS 2155-304 in the frame of the time dependent one-zone synchrotron self-Compton (SSC) model, where stochastic particle acceleration is taken into account. In this model, a homogeneously and isotropically spherical structure is assumed, the Fokker–Planck type equation which describes the evolution of the particles energy is numerically solved, and the synchrotron and self-Compton components from the spherical blob are calculated. Our results can reproduce observed spectra energy distribution (SED) and give definite predictions for the flux and spectral variability of PKS 2155-304. We find that particle injection rate, magnetic field and Doppler factor in the acceleration zone are important parameters for explaining its flaring behaviour.     

SONG-China Project: A Global Automated Observation Networktwo

Driven by the technological advancements and scientific objectives, the data acquisition in observational astronomy has been changed greatly in recent years. Fully automated or even autonomous ground-based network of telescopes has now become a tendency for time-domain observational projects. The Stellar Observations Network Group (SONG) is an international collaboration with the participation and contribution of the Chinese astronomy community. The scientific goal of SONG is time-domain astrophysics, such as the astroseismology and the research of variable stars in open clusters. The SONG project aims to build a global network of one-meter telescopes equipped with high-precision and high-resolution spectrographs, and two-channel lucky-imaging cameras. It is the Chinese initiative to install a 50 cm binocular photometry telescope at each SONG node to share the network platform and infrastructure. This work is focused on the design and implementation of SONG/50BiN in technology and methodology, for the ground-based network composed of multiple sites and a variety of instruments.     

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.     

Binary Star Fractions from the LAMOST DR4

Stellar systems composed of single, double, triple or higher-order systems are rightfully regarded as the fundamental building blocks of the Milky Way. Binary stars play an important role in formation and evolution of the Galaxy. Through comparing the radial velocity variations from multiepoch observations, we analyze the binary fraction of dwarf stars observed with LAMOST. Effects of different model assumptions, such as orbital period distributions on the estimate of binary fractions,are investigated. The results based on log-normal distribution of orbital periods reproduce the previous complete analyses better than the power-law distribution. We find that the binary fraction increases with Teff and decreases with [Fe/H]. We first investigate the relation between α-elements and binary fraction in such a large sample as provided by LAMOST. The old stars with high [α/Fe] dominate with a higher binary fraction than young stars with low [α/Fe]. At the same mass, earlier forming stars possess a higher binary fraction than newly forming ones, which may be related with evolution of the Galaxy.     

On the intrinsic shape of the gamma-ray spectrum for Fermi blazars

The curvature of the γ-ray spectrum in blazars may reflect the intrinsic distribution of emitting electrons, which will further give some information on the possible acceleration and cooling processes in the emitting region. The γ-ray spectra of Fermi blazars are normally fitted either by a single power-law(PL) or a log-normal(call Logarithmic Parabola, LP) form. The possible reason for this difference is not clear. We statistically explore this issue based on the different observational properties of 1419 Fermi blazars in the 3 LAC Clean Sample. We find that the γ-ray flux(100 Me V–100 Ge V) and variability index follow bimodal distributions for PL and LP blazars, where the γ-ray flux and variability index show a positive correlation. However, the distributions of γ-ray luminosity and redshift follow a unimodal distribution. Our results suggest that the bimodal distribution of γ-ray fluxes for LP and PL blazars may not be intrinsic and all blazars may have an intrinsically curved γ-ray spectrum, and the PL spectrum is just caused by the fitting effect due to less photons.     

Fermi/LAT observations of lobe-dominant radio galaxy 3C 207 and possible radiation region of γ-rays

3 C 207 is a lobe-dominant radio galaxy with a one sided jet and bright knots, spanning a kpcMpc scale, which have been resolved in the radio, optical and X-ray bands. This target was confirmed as a γ-ray emitter with Fermi/LAT, but it is uncertain whether the γ-ray emission region is the core or knots due to the low spatial resolution of Fermi/LAT. We present an analysis of its Fermi/LAT data acquired during the past 9 years. Different from the radio and optical emission from the core, it is found that theγ-ray emission is steady without detection of flux variation at over a 2σ confidence level. This likely implies that the γ-ray emission is from its knots. We collect the radio, optical and X-ray data of knot-A,the closest knot from the core at 1.4′′, and compile its spectral energy distribution(SED). Although the single-zone synchrotron+SSC+IC/CMB model that assumes knot-A is at rest can reproduce the SED in the radio-optical-X-ray band, the predicted γ-ray flux is lower than the LAT observations and the derived magnetic field strength deviates from the equipartition condition by 3 orders of magnitude. Assuming that knot-A is moving relativistically, its SED from radio to γ-ray bands would be represented well with the single-zone synchrotron+SSC+IC/CMB model under the equipartition condition. These results likely suggest that the γ-ray emission may be from knot-A via the IC/CMB process and the knot should have relativistical motion. The jet power derived from our model parameters is also roughly consistent with the kinetic power estimated with radio data.