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.     

Preprocessing photospheric vector magnetograms for nonlinear force-free field extrapolation of the global corona

Nonlinear force-free magnetic field(NLFFF) extrapolation based on the observed photospheric magnetic field is the most important method to obtain the coronal magnetic field nowadays.However, raw photospheric magnetograms contain magnetic forces and small-scale noises, and fail to be consistent with the force-free assumption of NLFFF models. The procedure for removing the forces and noises in observed data is called preprocessing. In this paper, we extend the preprocessing code of Jiang Feng to spherical coordinates for a full sphere. We first smooth the observed data with Gaussian smoothing, and then split the smoothed magnetic field into a potential field and a non-potential field.The potential part is computed by a numerical potential field model, and the non-potential part is preprocessed using an optimization method to minimize the magnetic forces and magnetic torques. Applying the code to synoptic charts of the vector magnetic field from SDO/HMI, we find it can effectively reduce the noises and forces, and improve the quality of data for a better input which will be used for NLFFF extrapolations applied to the global corona.     

On the internal structure and magnetic fields of

First, a sequence of four-zone models for the interior of Venus is constructed under the assumption of hydrostatic equilibrium. While the equation of state for each zone is taken to be the Bullen's relation with its coefficients consistent with the PREM Earth model (Dziewonski and Anderson, 1981), the position of core-mantle boundary is determined by matching solutions of the Emden's equation in different regions. The results of hydrostatic models indicate the presence of a reasonably large molten iron core in Venus, broadly similar to the Earth. It is also found that the position of the core-mantle interface is nearly model-independent. Second, we focus on the question why Venus does not possess a significant global magnetic field and on what we can learn from this fact. Solutions of magnetohydrodynamic equations appropriate for the molten core of Venus are discussed. It is argued that, because the Elsasser number measuring the relative importance of Coriolis and Lorentz forces satisfies 1, equations for the problem of thermal convection in the Venusian fluid core must be nearly uncoupled with the dynamo equation. The existence of a global magnetic field, though small, then suggests that the size of the magnetic Reynolds numberR _m must beR _m =O(10), sustaining a dynamo action near its marginal state but not an active dynamo in the Venusian molten core. On the basis of asymptotic relations for finite amplitude convection, a useful constraint on important physical parameters for the liquid core of Venus is derived and discussed.     

Numerical analysis of the symmetric methods

Aimed at the initial value problem of the particular second-order ordinary differential equations,y =f(x, y), the symmetric methods (Quinlan and Tremaine, 1990) and our methods (Xu and Zhang, 1994) have been compared in detail by integrating the artificial earth satellite orbits in this paper. In the end, we point out clearly that the integral accuracy of numerical integration of the satellite orbits by applying our methods is obviously higher than that by applying the same order formula of the symmetric methods when the integration time-interval is not greater than 12000 periods.