In this paper, we present the first light curve synthesis and orbital period change analysis of nine contact binaries around the short period limit. It is found that all these systems are W-subtype contact binaries. One of them is a medium contact system while the others are shallow contact ones. Four of them manifest obvious O'Connell effect explained by a dark spot or hot spot on one of the component stars. Third light was detected in three systems. By investigating orbital period variations, we found that four of the targets display a secular period decrease while the others exhibit a long-term period increase. The secular period decrease is more likely caused by angular momentum loss while the long-term period increase is due to mass transfer from the less massive component to the more massive one. Based on the statistic of 19 ultrashort period contact binaries with known orbital period changes, we found that seven of them display long-term decrease (three of them also exhibit cyclic variations), ten of them manifest long-term increase while two of them only show cyclic variation and that most of them are shallow contact binaries supporting the long timescale angular momentum loss theory suggested by Stepien. For the three deep contact systems, we found that they are probably triple systems. The tertiary companion plays an essential role during their formation and evolution.

Context: The identification and characterisation of populations of young massive stars in (giant) HII regions provides important constraints on i) the formation process of massive stars and their early feedback on the environment, and ii) the initial conditions for population synthesis models predicting the evolution of ensembles of stars. Aims: We identify and characterise the stellar populations of the following young giant HII regions: M8, G333.6-0.2, and NGC6357. Methods: We acquired H- and K-band spectra of around 200 stars using The K-band KMOS on the ESO Very Large Telescope. The targets for M8 and NGC6357 were selected from the MYStIX project, which combines X-ray observations with near-infrared and mid-infrared data. For G333.6-0.2, the sample selection is based on the near-infrared colours combined with X-ray data. We introduce an automatic spectral classification method in order to obtain temperatures and luminosities for the observed stars. We analyse the stellar populations using their photometric, astrometric, and spectroscopic properties and compared the position of the stars in the Hertzprung-Russell diagram with stellar evolution models to constrain their ages and mass ranges. Results: We confirm the presence of candidate ionising sources in the three regions and report new ones, including the first spectroscopically identified O stars in G333.6-0.2. In M8 and NGC6357, two populations are identified: i) OB main-sequence stars ($M > 5~\rm{M_{\odot}}$) and ii) pre-main sequence stars ($M\approx0.5-5~\rm{M_{\odot}}$). The ages of the clusters are $\sim$1-3~Myr, $< 3$~Myr, and $\sim$0.5-3~Myr for M8, G333.6-0.2, and NGC6357, respectively. We show that MYStIX selected targets have $>$ 90\% probability of being members of the HII region, whereas a selection based on near infrared (NIR) colours leads to a membership probability of only $\sim$70\%.

The precise measurement of the spectrum of protons, the most abundant component of the cosmic radiation, is necessary to understand the source and acceleration of cosmic rays in the Milky Way. This work reports the measurement of the cosmic ray proton fluxes with kinetic energies from 40 GeV to 100 TeV, with two and a half years of data recorded by the DArk Matter Particle Explorer (DAMPE). This is the first time an experiment directly measures the cosmic ray protons up to ~100 TeV with a high statistics. The measured spectrum confirms the spectral hardening found by previous experiments and reveals a softening at ~13.6 TeV, with the spectral index changing from ~2.60 to ~2.85. Our result suggests the existence of a new spectral feature of cosmic rays at energies lower than the so-called knee, and sheds new light on the origin of Galactic cosmic rays.

By analyzing two sets of complete $BVR_cI_c$ light curves of V342 UMa and three sets of complete $BVR_cI_c$ light curves of V509 Cam, we determined that the two systems are both W-subtype contact binaries and that V342 UMa shows a shallow contact configuration, while V509 Cam exhibits a medium contact configuration. Since both of them are totally eclipsing binaries, the physical parameters derived only by the photometric light curves are reliable.Meanwhile, the period changes of the two targets were analyzed based on all available eclipsing times. We discovered that V342 UMa shows long-term period decrease with a rate of $-1.02(\pm0.54)\times10^{-7}$ days/year and that V509 Cam exhibits long-term period increase with a rate of $3.96(\pm0.90)\times10^{-8}$ days/year. Both the conservative mass transfer and AML via magnetic stellar winds can interpret the long-term period decrease of V342 UMa. The long-term period increase of V509 Cam can be explained by mass transfer from the less massive star to the more massive one. The absolute parameters of the two binaries were estimated according to the Gaia distances and our derived photometric solution results. This method can be extended to other contact binaries without radial velocities but with reliable photometric solutions. The evolutionary states of them were discussed, we found that they reveal identical properties of other W-subtype contact systems.

GeV gamma ray is an important observation target of DArk Matter Particle Explorer (DAMPE) for indirect dark matter searching and high energy astrophysics. We present in this work a set of accurate instrument response functions of DAMPE (DmpIRFs) including the effective area, point-spread function and energy dispersion that are crucial for the gamma-ray data analysis based on the high statistics simulation data. A dedicated software named DmpST is developed to facilitate the scientific analyses of DAMPE gamma-ray data. Considering the limited number of photons and the angular resolution of DAMPE, the maximum likelihood method is adopted in the DmpST to better disentangle different source components. The basic mathematics and the framework regarding this software are also introduced in this paper.

The period distribution of contact binaries exhibits a very sharp short period cut-off at 0.22 days. In order to provide valuable information on this short period limit, we observed ten totally eclipsing contact binaries with orbital periods near this cut-off. By detailed analysis using W-D code, we determined that two of these systems are A-subtype contact binaries while the others are W-subtype contact binaries and all the targets show shallow contact configurations. Half of the targets exhibit stellar spot activity and four of them have third light. A statistical work on well studied USPCBs was carried out, and physical parameters of fifty-five USPCBs were obtained. Some common properties for these systems were derived. Due to the study of the period-color diagram of USPCBs, we found that the period-color relation of USPCBs is different from other W UMa type contact binaries and USPCBs are metal poor old stars. In addition, the evolutionary states of these systems were discussed by constructing the color-density diagram. We derived that the evolutionary states of the two components of USPCBs show identical characteristics of other contact binaries despite slower evolutionary status caused by smaller mass. We suggested that both the fully convective limit claimed by Rucinski (1992) and the AML theory recommended by Stepien (2006) can produce the short period cut-off and a tertiary companion is very important during the formation of the short period contact binaries by removing angular momentum from the host eclipsing pairs under certain circumstances.

One of the main purposes of the DArk Matter Particle Explorer (DAMPE) is to measure the cosmic ray nuclei up to several tens of TeV or beyond, whose origin and propagation remains a hot topic in astrophysics. The Plastic Scintillator Detector (PSD) on top of DAMPE is designed to measure the charges of cosmic ray nuclei from H to Fe and serves as a veto detector for discriminating gamma-rays from charged particles. We propose in this paper a charge reconstruction procedure to optimize the PSD performance in charge measurement. Essentials of our approach, including track finding, alignment of PSD, light attenuation correction, quenching and equalization correction are described detailedly in this paper after a brief description of the structure and operational principle of the PSD. Our results show that the PSD works very well and almost all the elements in cosmic rays from H to Fe are clearly identified in the charge spectrum.

DArk Matter Particle Explorer (DAMPE) is a general purpose high energy cosmic ray and gamma ray observatory, aiming to detect high energy electrons and gammas in the energy range 5 GeV to 10 TeV and hundreds of TeV for nuclei. This paper provides a method using machine learning to identify electrons and separate them from gammas,protons,helium and heavy nuclei with the DAMPE data from 2016 January 1 to 2017 June 30, in energy range from 10 to 100 GeV.

We have studied the variability of \object{S5 0716+714} at radio 15 GHz and $\gamma$-ray band using three different methods. A possible periodicity of $P_{15 GHz}=266.0\pm11.5$ and $P_{\gamma}=344.0\pm16.4$ days are obtained for radio 15 GHz and $\gamma$-ray light curves, respectively. The variability may be related to the intrinsically emission mechanism. The difference between the variability timescales of radio 15 GHz and $\gamma$-ray may be due to that the emission of radio 15 GHz is produced via the synchrotron process, while the $\gamma$-ray is produced by both the SSC and EC processes.

The DArk Matter Particle Explorer (DAMPE) is one of the four satellites within the Strategic Pioneer Research Program in Space Science of the Chinese Academy of Science (CAS). The Silicon-Tungsten Tracker (STK), which is composed of 768 singled-sided silicon microstrip detectors, is one of the four subdetectors in DAMPE, providing track reconstruction and charge identification for relativistic charged particles. The charge response of DAMPE silicon microstrip detectors is complicated, depending on the incident angle and impact position. A new charge reconstruction algorithm for the DAMPE silicon microstrip detector is introduced in this paper. This algorithm can correct the complicated charge response, and was proved applicable by the ion test beam.

The Kepler-9 system harbors three known transiting planets. The system holds significant interest for several reasons. First, the outer two planets exhibit a period ratio that is close to a 2:1 orbital commensurability, with attendant dynamical consequences. Second, both planets lie in the planetary mass "desert" that is generally associated with the rapid gas agglomeration phase of the core accretion process. Third, there exist attractive prospects for accurately measuring both the sky-projected stellar spin-orbit angles as well as the mutual orbital inclination between the planets in the system. Following the original \textit{Kepler} detection announcement in 2010, the initially reported orbital ephemerides for Kepler-9~b and c have degraded significantly, due to the limited time base-line of observations on which the discovery of the system rested. Here, we report new ground-based photometric observations and extensive dynamical modeling of the system. These efforts allow us to photometrically recover the transit of Kepler-9~b, and thereby greatly improve the predictions for upcoming transit mid-times. Accurate ephemerides of this system are important in order to confidently schedule follow-up observations of this system, for both in-transit Doppler measurements as well as for atmospheric transmission spectra taken during transit.

The DArk Matter Particle Explorer (DAMPE), also known as Wukong in China, launched on December 17, 2015, is a new high energy cosmic ray and {\gamma}-ray satellite-borne observatory in space. One of the main scientific goals of DAMPE is to observe GeV-TeV high energy {\gamma}-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 {\gamma}-rays, it is challenging to identify {\gamma}-rays with sufficiently high efficiency minimizing the amount of charged cosmic ray contamination. In this work we present a method to identify {\gamma}-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 {\gamma}-ray events at $\sim10$ GeV amounts to less than 1% of the selected sample. Finally, we use flight data to verify the effectiveness of the method by highlighting known {\gamma}-ray sources in the sky and by reconstructing preliminary light curves of the Geminga pulsar.

High energy cosmic ray electrons plus positrons (CREs), which lose energy quickly during their propagation, provide an ideal probe of Galactic high-energy processes and may enable the observation of phenomena such as dark-matter particle annihilation or decay. The CRE spectrum has been directly measured up to $\sim 2$ TeV in previous balloon- or space-borne experiments, and indirectly up to $\sim 5$ TeV by ground-based Cherenkov $\gamma$-ray telescope arrays. Evidence for a spectral break in the TeV energy range has been provided by indirect measurements of H.E.S.S., although the results were qualified by sizeable systematic uncertainties. Here we report a direct measurement of CREs in the energy range $25~{\rm GeV}-4.6~{\rm TeV}$ by the DArk Matter Particle Explorer (DAMPE) with unprecedentedly high energy resolution and low background. The majority of the spectrum can be properly fitted by a smoothly broken power-law model rather than a single power-law model. The direct detection of a spectral break at $E \sim0.9$ TeV confirms the evidence found by H.E.S.S., clarifies the behavior of the CRE spectrum at energies above 1 TeV and sheds light on the physical origin of the sub-TeV CREs.

The DArk Matter Particle Explorer (DAMPE) is one of the four satellites within Strategic Pioneer Research Program in Space Science of the Chinese Academy of Science (CAS). DAMPE can detect electrons, photons in a wide energy range (5 GeV to 10 TeV) and ions up to iron (100GeV to 100 TeV). Silicon-Tungsten Tracker (STK) is one of the four subdetectors in DAMPE, providing photon-electron conversion, track reconstruction and charge identification for ions. Ion beam test was carried out in CERN with 60GeV/u Lead primary beams. Charge reconstruction and charge resolution of STK detectors were investigated.

The DArk Matter Particle Explorer (DAMPE) is one of the four satellites within Strategic Pioneer Research Program in Space Science of the Chinese Academy of Science (CAS). DAMPE can detect electrons, photons and ions in a wide energy range (5 GeV to 10 TeV) and ions up to iron (100GeV to 100 TeV). Plastic Scintillator Detector (PSD) is one of the four payloads in DAMPE, providing e/{\gamma} separation and charge identification up to Iron. An ion beam test was carried out for the Qualification Model of PSD in CERN with 40GeV/u Argon primary beams. The Birk's saturation and charge resolution of PSD were investigated.

LX Ser is a deeply eclipsing cataclysmic variable with an orbital period of $0.^d 1584325$. Sixty two new eclipse times were determined by our observations and the AAVSO International Data base. Combining all available eclipse times, we analyzed the O-C behavior of LX Ser. We found that the O-C diagram of LX Ser shows a sinusoidal oscillation with a period of 22.8 yr and an amplitude of 0.00035 days. Two mechanisms (i.e., the Applegate mechanism and the light travel time effect) are applied to explain the cyclic modulation. We found that the Applegate mechanism is difficult to explain the cyclic oscillation in the orbital period. Therefore, the cyclic period change is most likely to be caused by the light travel time effect due to the presence of a third body. The mass of the tertiary component was determined to be $M_3\sim7.5 M_{Jup}$. We supposed that the tertiary companion is plausible a giant planet. The stability of the giant planet was checked, and we found that the multiple system is stable.

The Weihai Echelle Spectrograph (WES) is the first fiber-fed echelle spectrograph for astronomical observation in China. It is primarily used for chemical abundance and asteroseismology studies of nearby bright stars, as well as radial velocity detections for exoplanets. The optical design of WES is based on the widely demonstrated and well-established white-pupil concept. We describe the WES in detail and present some examples of its performance. A single exposure echelle image covers the spectral region 371-1,100 nm in 107 spectral orders over the rectangular CCD. The spectral resolution $R=\lambda/\Delta\lambda$ changes from 40,600 to 57,000 through adjusting the entrance slit width from full to 2.2 pixels sampling at the fiber-exit. The limiting magnitude scales to $V=8$ with a signal-to-noise ratio (SNR) of more than 100 in $V$ for an hour exposure, at the spectral resolution R$\approx$40,000 in the median seeing of 1.7$^{\prime\prime}$ at Weihai Observatory (WHO) for the 1-meter telescope. The radial velocity (RV) measurement accuracy of WES is estimated to be $<$10 m/s in 10 months (302 days) and better than 15 m/s in 4.4 years (1,617 days) in the recent data processing.

We have assembled the historical variability data of Mrk 421 at radio 15 GHz, X-ray and gamma-ray bands, spanning about 6.3, 10.3 and 7.5 yr, respective- ly. We analyzed the variability by using three methods. The results indicated that there is a period of 287.6 +- 4.4 days for 15 GHz, 309.5 +- 5.8 days for X- ray and 283.4+- 4.7 days for gamma-ray, respectively. This period can be reasonably explained by the nonballistic helical motion of the emitting material. The cor- relation analysis suggested that the variabilities of radio 15 GHz, X-ray and 0-ray are remarkable correlated, and the emission of radio 15 GHz lags behind that of X-ray, and the X-ray flux lags behind the gamma-ray. This suggests that the gamma-ray derives from inverse Compton (IC) scattering of the synchrotron photons, supporting the synchrotron self-Compton (SSC) model. Moreover, the time delay between different wavebands could be explained by the shock-in-jet models, in which a moving emission region produces the radio to gamma-ray activity, implying that the emission region of gamma-ray is closer to the center than ones of X-ray and radio emission.

In this paper, new determined BVR_cI_c light curves and radial velocities of V781 Tau are presented. By analyzing the light curves and radial velocities simultaneously, we found that V781 Tau is a W-subtype medium contact binary star with a mass ratio of q=2.207+-0.005 and a contact degree of f=21.6(+-1.0)%. The difference between the two light maxima was explained by a dark spot on the less massive primary component. The orbital period change of V781 Tau was also investigated. A secular decrease at a rate of $-6.01(+-2.28)*10^{-8} d/yr and a cyclic modulation with a period of 44.8+-5.7 yr and an amplitude of 0.0064+-0.0011 day were discovered. The continuous period decrease may be caused by angular momentum loss due to magnetic stellar wind. Applegate mechanism failed to explain the cyclic modulation. It is highly possible that the cyclic oscillation is the result of the light travel time effect by a third companion.

We present $GalevNB$ (Galev for $N$-body simulations), an utility that converts fundamental stellar properties of $N$-body simulations into observational properties using the $GALEV$ (GAlaxy EVolutionary synthesis models) package, and thus allowing direct comparisons between observations and $N$-body simulations. It works by converting fundamental stellar properties, such as stellar mass, temperature, luminosity and metallicity into observational magnitudes for a variety of filters of mainstream instruments/telescopes, such as HST, ESO, SDSS, 2MASS, etc.), and into spectra that spans from far-UV (90 $\rm \AA$) to near-IR (160 $\rm \mu$m). As an application, we use $GalevNB$ to investigate the secular evolution of spectral energy distribution (SED) and color-magnitude diagram (CMD) of a simulated star cluster over a few hundred million years. With the results given by $GalevNB$ we discover an UV-excess in the SED of the cluster over the whole simulation time. We also identify four candidates that contribute to the FUV peak, core helium burning stars, thermal pulsing asymptotic giant branch (TPAGB) stars, white dwarfs and naked helium stars.

The Active Particle-induced X-ray Spectrometer (APXS) is one of the payloads on board the Yutu rover of Chang'E-3 mission. In order to assess the instrumental performance of APXS, a ground verification test was done for two unknown samples (basaltic rock, mixed powder sample). In this paper, the details of the experiment configurations and data analysis method are presented. The results show that the elemental abundance of major elements can be well determined by the APXS with relative deviations < 15 wt. % (detection distance = 30 mm, acquisition time = 30 min). The derived detection limit of each major element is inversely proportional to acquisition time and directly proportional to detection distance, suggesting that the appropriate distance should be < 50mm.

X-ray spectrometer is one of the satellite payloads on Chang'E-2 satellite. The soft X-ray detector is one of the device on X-ray spectrometer which is designed to detect the major rock-forming elements within 0.5-10keV range on lunar surface. In this paper, energy linearity and energy resolution calibration is done using a weak Fe55 source, while temperature and time effect is considered not take big error. The total uncertainty is estimated to be within 5% after correction.

Under the steady state condition, the spectrum of electrons is investigated by solving the continuity equation under the complex radiation of both the synchrotron and Compton processes. The resulted GRB spectrum is a broken power law in both the fast and slow cooling phases. On the basis of this electron spectrum, the spectral indices of the Band function in four different phases are presented. In the complex radiation frame, the detail investigation on physical parameters reveals that both the reverse shock photosphere model and the forward shock with strong coupling model can answer the $\alpha \sim -1$ problem. A possible marginal to fast cooling phase transition in GRB 080916C is discussed. The time resolved spectra in different pulses of GRB 100724B, GRB 100826A and GRB 130606B are investigated. We found that the flux is proportional to the peak energy in almost all pulses. The phases for different pulses are determined according to the spectral index revolution. We found the strong correlations between spectral indices and the peak energy in GRB 100826A, which can be explained by the Compton effect in the fast cooling phase. However, the complex scenario predicts a steeper index for the injected electrons, which challenges the acceleration mechanism in GRBs.

We investigate the correlation between infrared (JHKL) and optical (B) fluxes of the variable nucleus of the Seyfert galaxy NGC 4151 using partially published data for the last 6 years (2008-2013.). Here we are using the same data as in Oknyansky et al. (2014), but include also optical (B) data from Guo et al. We find that the lag of flux in all the infrared bands is the same, 40 +- 6 days, to within the measurement accuracy. Variability in the J and K bands is not quite simultaneous, perhaps due to the differing contributions of the accretion disk in these bands. The lag found for the K band compared with the B band is not significantly different from earlier values obtained for the period 2000-2007. However, finding approximately the same lags in all IR bands for 2008-2013 differs from previous results at earlier epochs when the lag increased with increasing wavelength. Examples of almost the same lag in different IR bands are known for some other active nuclei. In the case of NGC 4151 it appears that the relative lags between the IR bands may be different in different years. The available data, unfortunately, do not allow us to investigate a possible change in the lags during the test interval. We discuss our results in the framework of the standard model where the variable infrared radiation is mainly due to thermal re-emission from the part of the dusty torus closest to the central source. There is also a contribution of some IR emission from the accretion disk, and this contribution increases with decreasing wavelength. Some cosmological applications of obtained results are discussed.

When observing the blazars, we identified a new eclipsing binary GSC 03517-00663. The light curves of GSC 03517-00663 are typical EW-type light curves. Based on the observation using the 1m telescope at Weihai Observatory of Shandong University, complete VRI light curves were determined. Then, we analyzed the multiple light curves using the W-D program. It is found that GSC 03517-00663 has a mass ratio of q=0.164 and a contact degree of f=69.2%. GSC 03517-00663 is a deep, low mass ratio overcontact binary. The light curves of GSC 03517-00663 show strong O'Connell effect, it was explained by employing a dark spot on the secondary component.

The first photometric analysis of IR Cas was carried out based on the new observed BVRI light curves. The symmetric light curves and nearly flat secondary minimum indicate that very precise photometric results can be determined. We found that IR Cas is a near contact binary with the primary component filling its Roche lobe. An analysis of the O-C diagram based on all available times of light minimum reveals evidence for a periodic change with a semiamplitude of 0.0153 days and a period of 39.7 years superimposed on a secular decrease at a rate of dp/dt=-1.28(\pm0.09)\times10^{-7} d yr^{-1}. The most reasonable explanation for the periodic change is the light time-travel effect due to a third body. The period decrease may be caused by mass transfer from the primary component to the secondary. With the decreasing period, IR Cas would eventually evolve into a contact system.

When observing the transiting extrasolar planets, we discovered a new W UMa type binary GSC 03553-00845. Following observation was carried out using the 1m telescope at Weihai Observatory of Shandong University. Complete BVR light curves were determined. Using the W-D program, we analyzed the light curves. Photometric solution reveals that GSC 03553-00845 is a W-subtype W UMa type binary with a mass ratio of q=2.904, it is an overcontact binary system by a contact degree of f=29.5% with a small temperature difference between the components (\Delta T=206 K indicating a good thermal contact between the components. More observation of GSC 03553-00845 is needed in order to analyze the light curve variation and orbital period change.

In this paper, a total of about 28000 images in $V$ and $R$ band obtained on 161 nights using the one-meter optical telescope at Weihai Observatory (WHO) of Shandong University since 2008 to 2012 have been processed to measure the sky brightness. It provides us with an unprecedented database, which can be used to study the variation of the sky brightness with the sky position, the moonlight contribution, and the twilight sky brightness. The darkest sky brightness is about 19.0 and 18.6 $mag$ $arcsec^{-2}$ in $V$ and $R$ band, respectively. An obvious darkening trend is found at the first half of the night at WHO, and the variation rate is much larger in summer than that in other seasons. The sky brightness variation depends more on the azimuth than on the altitude of the telescope pointing for WHO. Our results indicate that the sky brightness at WHO is seriously influenced by the urban light.

We reported the monitoring results of OT 546 in V, R and I bands, observed on 22 nights from February 16 to July 1 in 2009 at Weihai Observatory of Shandong University. During our monitoring, its variability amplitude was small and a possible microvariability was detected on one night using both C and F tests

A new generation detector for the high energy cosmic ray - the DAMPE(DArk Matter Particle Explorer) is a satellite based project. Its main object is the measurement of energy spectrum of cosmic ray nuclei from 100GeV to 100TeV, the high energy electrons and gamma ray from 5GeV to 10TeV. A silicon matrix detector described in this paper, is employed for the sea level cosmic ray energy and position detection while the prototype testing of the DAMPE. This matrix is composed by the 180 silicon PIN detectors, which covers an area of 32*20 cm2. The primary testing results are shown including MIPs energy spectrum and the position sensitive map.

We report B, V, and R band CCD photometry of the Seyfert galaxy NGC 4151 obtained with the 1.0-m telescope at Weihai Observatory of Shandong University and the 1.56-m telescope at Shanghai Astronomical Observatory from 2005 December to 2013 February. Combining all available data from literature, we have constructed a historical light curve from 1910 to 2013 to study the periodicity of the source using three different methods (the Jurkevich method, the Lomb-Scargle periodogram method and the Discrete Correlation Function method). We find possible periods of P_1=4\pm0.1, P_2=7.5\pm0.3 and P_3=15.9\pm0.3 yr.

Well-sampled optical and radio light curves of BL Lacertae in B, V, R, I bands and 4.8, 8.0, 14.5 GHz from 1968 to 2014 were presented in this paper. A possible $1.26 \pm 0.05$ yr period in optical bands and a $7.50 \pm 0.15$ yr period in radio bands were detected based on discrete correlation function, structure function as well as Jurkevich method. Correlations among different bands were also analyzed and no reliable time delay was found between optical bands. Very weak correlations were detected between V band and radio bands. However, in radio bands the variation at low frequency lagged that at high frequency obviously. The spectrum of BL Lacertae turned mildly bluer when the object turned brighter, and stronger bluer-when-brighter trends were found for short flares. A scenario including a precessing helical jet and periodic shocks was put forward to interpret the variation characteristics of BL Lacertae.

A radiation gene box (RGB) onboard the SJ-10 satellite is a device carrying mice and drosophila cells to determine the biological effects of space radiation environment. The shielded fluxes of different radioactive sources were calculated and the linear energy transfers of gamma-rays, electrons, protons and alpha-particles in tissue were acquired using A-150 tissue-equivalent plastic. Then, a conceptual model of a space radiation instrument employing three semiconductor sub-detectors for deriving the charged and uncharged radiation environment of the RGB was designed. The energy depositions in the three sub-detectors were classified into fifteen channels (bins) in an algorithm derived from the Monte Carlo method. The physical feasibility of the conceptual instrument was also verified by Monte Carlo simulations.

The one-meter telescope at Weihai Observatory of Shandong University is an f/8 Cassegrain telescope. Three sets of filters, including Johnson--Cousins UBVRI, Sloan Digital Sky Survey u$'$g$'$r$'$i$'$z$'$ and Str\"{o}mgren uvby, are installed in a dual layer filterwheel. The photometric system and the CCD camera are introduced in this paper, followed by detailed analysis of their performances, and determination of the relevant parameters, including gain, readout noise, dark current and linearity of the CCD camera. In addition, the characteristics of the site astro-climate condition, including typical seeing, clear nights statistics, and average sky brightness were studied systematically based on data gathered from Sep. 2007 to Aug. 2013, and were reported in this work. Photometric calibrations were done using 8 nights Landolt standard star observations, which yielded transformation coefficients, photometry precision and system throughput. The limiting magnitudes were simulated using the derived calibration parameters and classic observation conditions at WHO.

The BL Lac object Mrk 421 was observed in optical bands from 2009 April to 2012 May with the 1.0 m telescope at Weihai Observatory of Shandong University. Microvariability was analysed by C and F tests, but no significant microvariability was detected during our observations.

OI 090.4 was monitored on 21 nights from 2006 to 2012 for studying the variability. Strong variations occurred during the past 6 years. The long-term variability amplitude is consistent with previous results. Microvariability was analyzed for 43 intra-night light curves. 30 out of 43 light curves showed microvariability by C and F tests analysis.