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.

We report on the results of a multiwavelength monitoring campaign of the bright, nearby Seyfert galaxy, Ark 120 using a ~50-day observing programme with Swift and a ~4-month co-ordinated ground-based observing campaign, predominantly using the Skynet Robotic Telescope Network. We find Ark 120 to be variable at all optical, UV, and X-ray wavelengths, with the variability observed to be well-correlated between wavelength bands on short timescales. We perform cross-correlation analysis across all available wavelength bands, detecting time delays between emission in the X-ray band and the Swift V, B and UVW1 bands. In each case, we find that the longer-wavelength emission is delayed with respect to the shorter-wavelength emission. Within our measurement uncertainties, the time delays are consistent with the \tau ~ \lambda^{4/3} relation, as predicted by a disc reprocessing scenario. The measured lag centroids are \tau_{cent} = 11.90 +/- 7.33, 10.80 +/- 4.08, and 10.60 +/- 2.87 days between the X-ray and V, B, and UVW1 bands, respectively. These time delays are longer than those expected from standard accretion theory and, as such, Ark 120 may be another example of an active galaxy whose accretion disc appears to exist on a larger scale than predicted by the standard thin-disc model. Additionally, we detect further inter-band time delays: most notably between the ground-based I and B bands (\tau_{cent} = 3.46 +/- 0.86 days), and between both the Swift XRT and UVW1 bands and the I band (\tau_{cent} = 12.34 +/- 4.83 and 2.69 +/- 2.05 days, respectively), highlighting the importance of co-ordinated ground-based optical observations.

We report the lens mass and distance measurements of the nearby microlensing event TCP J05074264+2447555. We measure the microlens parallax vector ${\pi}_{\rm E}$ using Spitzer and ground-based light curves with constraints on the direction of lens-source relative proper motion derived from Very Large Telescope Interferometer (VLTI) GRAVITY observations. Combining this ${\pi}_{\rm E}$ determination with the angular Einstein radius $\theta_{\rm E}$ measured by VLTI GRAVITY observations, we find that the lens is a star with mass $M_{\rm L} = 0.495 \pm 0.063~M_{\odot}$ at a distance $D_{\rm L} = 429 \pm 21~{\rm pc}$. We find that the blended light basically all comes from the lens. The lens-source proper motion is $\mu_{\rm rel,hel} = 26.55 \pm 0.36~{\rm mas\,yr^{-1}}$, so with currently available adaptive-optics (AO) instruments, the lens and source can be resolved in 2021. This is the first microlensing event whose lens mass is unambiguously measured by interferometry + satellite parallax observations, which opens a new window for mass measurements of isolated objects such as stellar-mass black holes.

We present a model atom for C I - C II - C III - C IV using the most up-to-date atomic data and evaluated the non-local thermodynamic equilibrium (NLTE) line formation in classical 1D atmospheric models of O-B-type stars. Our models predict the emission lines of C II 9903~\AA\ and 18535~\AA\ to appear at effective temperature \Teff~$\geq$~17\,500~K, those of C II 6151~\AA\ and 6461~\AA\ to appear at \Teff~$>$~25\,000~K, and those of C III 5695, 6728--44, 9701--17~\AA\ to appear at \Teff~$\geq$~35\,000~K (log~$g$=4.0). Emission occurs in the lines of minority species, where the photoionization-recombination mechanism provides a depopulation of the lower levels to a greater extent than the upper levels. For C II 9903 and 18535~\AA, the upper levels are mainly populated from C III reservoir through the Rydberg states. For C III 5695 and 6728--44~\AA, the lower levels are depopulated due to photon losses in UV transitions at 885, 1308, and 1426--28~\AA\ which become optically thin in the photosphere. We analysed the lines of C I, C II, C III, and C IV for twenty-two O-B-type stars with temperature range between 15\,800 $\leq$~\Teff~$\leq$ 38\,000~K. Abundances from emission lines of C I, C II and C III are in agreement with those from absorption ones for most of the stars. We obtained log~$\epsilon_{\rm C}$=8.36$\pm$0.08 from twenty B-type stars, that is in line with the present-day Cosmic Abundance Standard. The obtained carbon abundances in 15~Mon and HD~42088 from emission and absorption lines are 8.27$\pm$0.11 and 8.31$\pm$0.11, respectively.

The typical blazar S5 0716$+$714 is very interesting due to its rapid and large amplitude variability and high duty cycle of micro-variability in optical band. We analyze the observations in I, R and V bands obtained with the $1.0m$ telescope at Weihai observatory of Shandong University from 2011 to 2018. The model of synchrotron radiation from turbulent cells in a jet has been proposed as a mechanism for explaining micro-variability seen in blazar light curves. Parameters such as the sizes of turbulent cells, the enhanced particle densities, and the location of the turbulent cells in the jet can be studied using this model. The model predicts a time lag between variations as observed in different frequency bands. Automatic model fitting method for micro-variability is developed, and the fitting results of our multi-frequency micro-variability observations support the model. The results show that both the amplitude and duration of flares decomposed from the micro-variability light curves confirm to the log-normal distribution. The turbulent cell size is within the range of about 5 to 55 AU, and the time lags of the micro-variability flares between the I-R and R-V bands should be several minutes. The time lags obtained from the turbulence model are consistent with the fitting statistical results, and the time lags of flares are correlated with the time lags of the whole light curve.

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.

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.

We collect long term $\gamma$-ray, optical and radio $15$ GHz light curves of quasar object PMN J2345-1555. The correlation analyses between them are performed via the local cross-correlation function (LCCF). We found that all the optical $V$, $R$ band and the infrared $J$ band are correlated with the radio 15 GHz at beyond $3\sigma$ significance level, and the lag times are $-221.81^{+6.26}_{-6.72}$, $-201.38^{+6.42}_{-6.02}$ and $-192.27^{+8.26}_{-7.37}$ days, respectively. The $\gamma$-ray is strongly correlated with optical, but weakly correlated with the radio. We present that time lags between different frequencies can be used as an alternative parameter to derive the core-shift measurement. For this target, the magnetic field and particle density at 1 parsec in jet are derived to be $0.61$ Gauss and $1533/\gamma_{\rm min}$ cm$^{-3}$, respectively. The black hole mass and the 15 GHz core position in jet are estimated to be $10^{8.44} {\rm M}_{\odot}$ and $30$ parsec, respectively. The lag times enable us to derive that the optical and the $\gamma$-ray emitting regions coincide, which are located at $4.26^{+0.83}_{-0.79}$ pc away from 15 GHz core position in jet and beyond the broad line region (BLR). We found that a $3\sigma$ correlation between the color index and the radio light curve, which indicates that opacity may play an important role in the variation. The $\delta V-\delta R$ behaviors are complex, while the $R-J$ shows a bluer when brighter trend. As hinted from radio images, we proposed a positional dependent spectral index model to explain the color index behaviors, which is complementary for the shock in jet model. The curvature effects and contribution from accretion disk may also affect variables of blazars in many aspects.

We report on a recent multi-band optical photometric and polarimetric observational campaign of the blazar OJ 287 which was carried out during September 2016 -- December 2017. We employed nine telescopes in Bulgaria, China, Georgia, Japan, Serbia, Spain and the United States. We collected over 1800 photometric image frames in BVRI bands and over 100 polarimetric measurements over ~175 nights. In 11 nights with many quasi-simultaneous multi-band (V, R, I) observations, we did not detect any genuine intraday variability in flux or color. On longer timescales, multiple flaring events were seen. Large changes in color with respect to time and in a color--magnitude diagram were seen, and while only a weak systematic variability trend was noticed in color with respect to time, the color--magnitude diagram shows a bluer-when-brighter trend. Large changes in the degree of polarization, and substantial swings in the polarization angle were detected. The fractional Stokes parameters of the polarization showed a systematic trend with time in the beginning of these observations, followed by chaotic changes and then an apparently systematic variation at the end. These polarization changes coincide with the detection and duration of the source at very high energies as seen by VERITAS. The spectral index shows a systematic variation with time and V-band magnitude. We briefly discuss possible physical mechanisms that could explain the observed flux, color, polarization, and spectral variability.

The enormous sizes and variability of emission of radio-loud AGNs arise from the relativistic flows of plasma along two oppositely directed jets. We use the Athena hydrodynamics code to simulate an extensive suite of 54 propagating three-dimensional relativistic jets with wide ranges of input jet velocities and jet-to-ambient matter density ratios. We determine which parameter sets yield unstable jets that produce jet dominated FR I type radio galaxy morphologies and which tend to produce stable jets with hot-spots and FR II morphologies. Nearly all our simulations involve jets with internal pressures matched to those of the ambient medium but we also consider over-pressured jets and discuss differences from the standard ones. We also show that the results are not strongly dependent on the adiabatic index of the fluid. We focus on simulations that remain stable for extended distances (60-240) times the initial jet radius. Scaled to the much smaller sizes probed by VLBI observations, the fluctuations in such simulated flows yield variability in observed emissivity on timescales from months. Adopting results for the densities, pressures and velocities from these simulations we estimate normalized rest frame synchrotron emissivities from individual cells in the jets. The observed emission from each cell is strongly dependent upon its variable Doppler boosting factor. We sum the fluxes from thousands of zones around the primary reconfinement shock. The light curves an power-spectra, with red-noise slopes between -2.1 and -2.5, so produced are similar to those observed from blazars.

Non-local thermodynamical equilibrium (NLTE) line formation for Mg I and Mg II lines is considered in classical 1D-LTE model atmospheres of the Sun and 17 stars with reliable atmospheric parameters and in a broad range of spectral types: 3900 $\le$ \Teff\ $\le$ 17500~K, 1.1 $\le$ log$g$ $\le$ 4.7, and $-$2.6 $\le$ [Fe/H] $\le$ +0.4. We find that, for each star, NLTE leads to smaller line-to-line scatter. For 10 stars, NLTE leads to consistent abundances from Mg I and Mg II, while the difference in the LTE abundance varies between $-0.21$ and +0.23~dex. We obtain an abundance discrepancy betweeen Mg I and Mg II in the two very metal-poor stars, HD 140283 and HD 84937. An origin of these abundance differences remains unclear. Our standard NLTE modelling predicts Mg I emission lines at 7.736, 11.789, 12.224, and 12.321 $\mu$m in the atmospheres with \Teff\ $\le$ 7000 K. We reproduce well the Mg I 12.2 and 12.3 $\mu$m emission lines in Procyon. However, for the Sun and 3 K-giants, the predicted Mg I emission lines are too weak compared with the observations. For stars with 7000 K $\le$ \Teff\ $\leq$ 17500 K, we recommend the Mg II 3848, 3850, 4384, 4390, 4427, and 4433 \AA\ lines for Mg abundance determinations even at the LTE assumption due to their small NLTE effects. The Mg I 4167, 4571, 4702, 5528, 5167, 5172, and 5183 \AA\ lines can be safely used in the LTE analysis of stars with 7000 K $<$ \Teff\ $\leq$ 8000 K. For the hotter stars, with \Teff\ from 8000 to 9500 K, the NLTE effects are minor only for Mg\ione\ 4167, 4702, and 4528 \AA.

We report the photometry of six transits of the hot Jupiter HAT-P-29b obtained from 2013 October to 2015 January. We analyze the new light curves, in combination with the published photometric, and Doppler velocimetric, and spectroscopic measurements, finding an updated orbital ephemeris for the HAT-P-29 system, $T_{\rm C}[0]= 2456170.5494(15)\,[\rm{BJD_{TDB}}]$ and $P=5.723390(13) \,{\rm days}$. It is $17.63\,{\rm s}$ ($4.0\,\sigma$) longer than the previously published value, amounting to errors exceeding $2.5\,\rm{hrs}$ at the time of writing (on UTC 2018 June 1). The measured transit mid-times for HAT-P-29b show no compelling evidence of timing anomalies from a linear model, which rules out the presence of a perturbers with masses greater than $0.6$, $0.7$, $0.5$, and $0.4\,{\rm M_\oplus}$ near the $1:2$, $2:3$, $3:2$, and $2:1$ resonances with HAT-P-29b, respectively.

We present a multi-wavelength spectral and temporal investigation of OJ 287 emission during its strong optical-to-X-ray activity between July 2016 - July 2017. The daily $\gamma$-ray fluxes from \emph{Fermi}-LAT are consistent with no variability. The strong optical-to-X-ray variability is accompanied by a change in power-law spectral index of the X-ray spectrum from $< 2$ to $>2$, with variations often associated with changes in optical polarization properties. Cross-correlations between optical-to-X-ray emission during four continuous segments show simultaneous optical-ultraviolet (UV) variations while the X-ray and UV/optical are simultaneous only during the middle two segments. In the first segment, the results suggest X-rays lag the optical/UV, while in the last segment X-rays lead by $\sim$ 5-6 days. The last segment also shows a systematic trend with variations appearing first at higher energies followed by lower energy ones. The LAT spectrum before the VHE activity is similar to preceding quiescent state spectrum while it hardens during VHE activity period and is consistent with the extrapolated VHE spectrum during the latter. Overall, the broadband spectral energy distributions (SEDs) during high activity periods are a combination of a typical OJ 287 SED and an HBL SED, and can be explained in a two-zone leptonic model, with the second zone located at parsec scales, beyond the broad line region, being responsible for the HBL-like spectrum. The change of polarization properties from systematic to chaotic and back to systematic, before, during and after the VHE activity, suggest dynamic roles for magnetic fields and turbulence.

Within the standard hierarchical structure formation scenario, Milky Way-mass dark matter haloes have hundreds of dark matter subhaloes with mass $\gtrsim 10^8 \, {\rm M_{\odot}}$. Over the lifetime of a galactic disc a fraction of these may pass close to the central region and interact with the disc. We extract the properties of subhaloes, such as their mass and trajectories, from a realistic cosmological simulation to study their potential effect on stellar discs. We find that massive subhalo impacts can generate disc heating, rings, bars, warps, lopsidedness as wells as spiral structures in the disc. Specifically, strong counter-rotating single-armed spiral structures form each time a massive subhalo passes through the disc. Such single-armed spirals wind up relatively quickly (over $1-2$ Gyrs) and are generally followed by co-rotating two-armed spiral structures that both develop and wind up more slowly. In our simulations self-gravity in the disc is not very strong and these spiral structures are found to be kinematic density waves. We demonstrate that there is a clear link between each spiral mode in the disc and a given subhalo that caused it, and by changing the mass of the subhalo we can modulate the strength of the spirals. Furthermore, we find that the majority of subhaloes interact with the disc impulsively, such that the strength of spirals generated by subhaloes is proportional to the total torque they exert. We conclude that only a handful of encounters with massive subhaloes is sufficient for re-generating and sustaining spiral structures in discs over their entire lifetime.

We present eight new light curves of the transiting extra-solar planet HAT-P-25b obtained from 2013 to 2016 with three telescopes at two observatories. We use the new light curves, along with recent literature material, to estimate the physical and orbital parameters of the transiting planet. Specifically, we determine the mid-transit times (T$_{C}$) and update the linear ephemeris, T$_{C[0]}$=2456418.80996$\pm$0.00025 [$\mathrm{BJD}_\mathrm{TDB}$] and P=3.65281572$\pm$0.00000095 days. We carry out a search for transit timing variations (TTVs), and find no significant TTV signal at the $\Delta T=$80 s-level, placing a limit on the possible strength of planet-planet interactions ($\mathrm{TTV_{G}}$). In the course of our analysis, we calculate the upper mass-limits of the potential nearby perturbers. Near the 1:2, 2:1, and 3:1 resonances with HAT-P-25b, perturbers with masses greater than 0.5, 0.3, and 0.5 $\mathrm{M_{\oplus}}$ respectively, can be excluded. Furthermore, based on the analysis of TTVs caused by light travel time effect (LTTE) we also eliminate the possibility that a long-period perturber exists with $M_{\rm p}> 3000 \,\mathrm{M_{J}}$ within $a=11.2\,{\rm AU}$ of the parent star.

We analyzed Kepler light curves of the blazar W2R 1926+42 that provided nearly continuous coverage from quarter 11 through quarter 17 (589 days between 2011 and 2013) and examined some of their flux variability properties. We investigate the possibility that the light curve is dominated by a large number of individual flares and adopt exponential rise and decay models to investigate the symmetry properties of flares. We found that those variations of W2R 1926+42 are predominantly asymmetric with weak tendencies toward positive asymmetry (rapid rise and slow decay). The durations (D) and the amplitudes (F0) of flares can be fit with log-normal distributions. The energy (E) of each flare is also estimated for the first time. There are positive correlations between logD and logE with a slope of 1.36, and between logF0 and logE with a slope of 1.12. Lomb-Scargle periodograms are used to estimate the power spectral density (PSD) shape. It is well described by a power law with an index ranging between -1.1 and -1.5. The sizes of the emission regions, R, are estimated to be in the range of 1.1*10^15 cm - 6.6*10^16 cm. The flare asymmetry is difficult to explain by a light travel time effect but may be caused by differences between the timescales for acceleration and dissipation of high-energy particles in the relativistic jet. A jet-in-jet model also could produce the observed log-normal distributions.

The optical images of near-Earth asteroid 4179 Toutatis acquired by Chang'e-2 spacecraft show that Toutatis has an elongated contact binary configuration, with the contact point located along the long axis. We speculate that such configuration may have resulted from a low-speed impact between two components. In this work, we performed a series of numerical simulations and compared the results with the optical images, to examine the mechanism and better understand the formation of Toutatis. Herein we propose an scenario that an assumed separated binary precursor could undergo a close encounter with Earth, leading to an impact between the primary and secondary, and the elongation is caused by Earth's tide. The precursor is assumed to be a doubly synchronous binary with a semi-major axis of 4 $R_p$ (radius of primary) and the two components are represented as spherical cohesionless self-gravitating granular aggregates. The mutual orbits are simulated in a Monte Carlo routine to provide appropriate parameters for our \emph{N}-body simulations of impact and tidal distortion. We employ the \emph{pkdgrav} package with a soft-sphere discrete element method (SSDEM) to explore the entire scenarios. The results show that contact binary configurations are natural outcomes under this scenario, whereas the shape of the primary is almost not affected by the impact of the secondary. However, our simulations further provide an elongated contact binary configuration best-matching to the shape of Toutatis at an approaching distance $r_p$ = 1.4 $\sim$ 1.5 $R_e$ (Earth radius), indicative of a likely formation scenario for configurations of Toutatis-like elongated contact binaries.

TianQin is a space-based laser interferometric gravitational wave detector aimed at detecting gravitational waves at low frequencies (0.1 mHz -- 1 Hz). It is formed by three identical drag-free spacecrafts in an equilateral triangular constellation orbiting around the Earth. The distance between each pair of spacecrafts is approximately $1.7 \times 10^{5} ~\rm{km}$. The spacecrafts are interconnected by infrared laser beams forming up to three Michelson-type interferometers. The detailed mission design and the study of science objectives for the TianQin project depend crucially on the orbit and the response of the detector. In this paper, we provide the analytic expressions for the coordinates of the orbit for each spacecraft in the heliocentric-ecliptic coordinate system to the leading orders. This enables a sufficiently accurate study of science objectives and data analysis, and serves as a first step to further orbit design and optimization. We calculate the response of a single Michelson detector to plane gravitational waves in arbitrary waveform which is valid in the full range of the sensitive frequencies. It is then used to generate the more realistic sensitivity curve of TianQin. We apply this model on a reference white-dwarf binary as a proof of principle.

Hydrogen-poor superluminous supernovae (SLSNe-I) have been predominantly found in low-metallicity, star-forming dwarf galaxies. Here we identify Gaia17biu/SN 2017egm as an SLSN-I occurring in a "normal" spiral galaxy (NGC 3191) in terms of stellar mass (several times 10^10 M_sun) and metallicity (roughly Solar). At redshift z=0.031, Gaia17biu is also the lowest redshift SLSN-I to date, and the absence of a larger population of SLSNe-I in dwarf galaxies of similar redshift suggests that metallicity is likely less important to the production of SLSNe-I than previously believed. With the smallest distance and highest apparent brightness for an SLSN-I, we are able to study Gaia17biu in unprecedented detail. Its pre-peak near-ultraviolet to optical color is similar to that of Gaia16apd and among the bluest observed for an SLSN-I while its peak luminosity (M_g = -21 mag) is substantially lower than Gaia16apd. Thanks to the high signal-to-noise ratios of our spectra, we identify several new spectroscopic features that may help to probe the properties of these enigmatic explosions. We detect polarization at the ~0.5% level that is not strongly dependent on wavelength, suggesting a modest, global departure from spherical symmetry. In addition, we put the tightest upper limit yet on the radio luminosity of an SLSN-I with <5.4x10^26 erg/s/Hz (at 10 GHz), which is almost a factor of 40 better than previous upper limits and one of the few measured at an early stage in the evolution of an SLSN-I. This limit largely rules out an association of this SLSNe-I with known populations of gamma-ray burst (GRB) like central engines.

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.

In asteroid rendezvous missions, the dynamical environment near the asteroid's surface should be made clear prior to the mission launch. However, most of the asteroids have irregular shapes, which lower the efficiency of calculating their gravitational field by adopting the traditional polyhedral method. In this work, we propose a method to partition the space near the asteroid adaptively along three spherical coordinates and use Chebyshev polynomials interpolation to represent the gravitational acceleration in each cell. Moreover, we compare four different interpolation schemes to obtain the best precision with the identical initial parameters. An error-adaptive octree division is combined to improve the interpolation precision near the surface. As an example, we take the typical irregular-shaped near-Earth asteroid 4179 Toutatis to show the advantage of this method, as a result, we show that the efficiency can be increased by hundreds to thousands times with our method. In a word, this method can be applicable to other irregular-shaped asteroids and can greatly improve the evaluation efficiency.

The Doppler tracking data of the Chang'e 3 lunar mission is used to constrain the stochastic background of gravitational wave in cosmology within the 1 mHz to 0.05 Hz frequency band. Our result improves on the upper bound on the energy density of the stochastic background of gravitational wave in the 0.02 Hz to 0.05 Hz band obtained by the Apollo missions, with the improvement reaching almost one order of magnitude at around 0.05 Hz. Detailed noise analysis of the Doppler tracking data is also presented, with the prospect that these noise sources will be mitigated in future Chinese deep space missions. A feasibility study is also undertaken to understand the scientific capability of the Chang'e 4 mission, due to be launched in 2018, in relation to the stochastic gravitational wave background around 0.01 Hz. The study indicates that the upper bound on the energy density may be further improved by another order of magnitude from the Chang'e 3 mission, which will fill the gap in the frequency band from 0.02 Hz to 0.1 Hz in the foreseeable future.

We carried out a pilot campaign of radio and optical band intra-day variability (IDV) observations of five blazars (3C66A, S5 0716+714, OJ287, B0925+504, and BL Lacertae) on December 18--21, 2015 by using the radio telescope in Effelsberg (Germany) and several optical telescopes in Asia, Europe, and America. After calibration, the light curves from both 5 GHz radio band and the optical R band were obtained, although the data were not smoothly sampled over the sampling period of about four days. We tentatively analyse the amplitudes and time scales of the variabilities, and any possible periodicity. The blazars vary significantly in the radio (except 3C66A and BL Lacertae with only marginal variations) and optical bands on intra- and inter-day time scales, and the source B0925+504 exhibits a strong quasi-periodic radio variability. No significant correlation between the radio- and optical-band variability appears in the five sources, which we attribute to the radio IDV being dominated by interstellar scintillation whereas the optical variability comes from the source itself. However, the radio- and optical-band variations appear to be weakly correlated in some sources and should be investigated based on well-sampled data from future observations.

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.

Here we report on the results of the WEBT photo-polarimetric campaign targeting the blazar S5~0716+71, organized in March 2014 to monitor the source simultaneously in BVRI and near IR filters. The campaign resulted in an unprecedented dataset spanning $\sim 110$\,h of nearly continuous, multi-band observations, including two sets of densely sampled polarimetric data mainly in R filter. During the campaign, the source displayed pronounced variability with peak-to-peak variations of about $30\%$ and "bluer-when-brighter" spectral evolution, consisting of a day-timescale modulation with superimposed hourlong microflares characterized by $\sim 0.1$\,mag flux changes. We performed an in-depth search for quasi-periodicities in the source light curve; hints for the presence of oscillations on timescales of $\sim 3$\,h and $\sim 5$\,h do not represent highly significant departures from a pure red-noise power spectrum. We observed that, at a certain configuration of the optical polarization angle relative to the positional angle of the innermost radio jet in the source, changes in the polarization degree led the total flux variability by about 2\,h; meanwhile, when the relative configuration of the polarization and jet angles altered, no such lag could be noted. The microflaring events, when analyzed as separate pulse emission components, were found to be characterized by a very high polarization degree ($> 30\%$) and polarization angles which differed substantially from the polarization angle of the underlying background component, or from the radio jet positional angle. We discuss the results in the general context of blazar emission and energy dissipation models.

We employ very high resolution simulations of isolated Milky Way-like galaxies to study the effect of triaxial dark matter haloes on exponential stellar discs. Non-adiabatic halo shape changes can trigger two-armed grand-design spiral structures which extend all the way to the edge of the disc. Their pattern speed coincides with the inner Lindblad resonance indicating that they are kinematic density waves which can persist up to several Gyr. In dynamically cold discs, grand-design spirals are swing amplified and after a few Gyr can lead to the formation of (multi-armed) transient recurrent spirals. Stellar discs misaligned to the principal planes of the host triaxial halo develop characteristic integral shaped warps, but otherwise exhibit very similar spiral structures as aligned discs. For the grand-design spirals in our simulations, their strength dependence with radius is determined by the torque on the disc, suggesting that by studying grand-design spirals without bars it may be possible to set constraints on the tidal field and host dark matter halo shape.

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.

TianQin is a proposal for a space-borne detector of gravitational waves in the millihertz frequencies. The experiment relies on a constellation of three drag-free spacecraft orbiting the Earth. Inter-spacecraft laser interferometry is used to monitor the distances between the test masses. The experiment is designed to be capable of detecting a signal with high confidence from a single source of gravitational waves within a few months of observing time. We describe the preliminary mission concept for TianQin, including the candidate source and experimental designs. We present estimates for the major constituents of the experiment's error budget and discuss the project's overall feasibility. Given the current level of technology readiness, we expect TianQin to be flown in the second half of the next decade.

This paper presents analysis of the rotational parameters of Toutatis based on the observational results from Chang'e-2's close flyby. The 3-D shape model derived from ground-based radar observation is used to calculate the 3-1-3 Euler angles at the flyby epoch, which are evaluated to be $-20.1^\circ\pm1^\circ$, $27.6^\circ\pm1^\circ$ and $42.2^\circ\pm1^\circ$. The large amplitude of Toutatis' tumbling attitude is demonstrated to be the result of the large deviation of the angular momentum axis and the rotational axis. Two rotational periods are evaluated to be $5.38\pm0.03$ days for rotation about the long axis and $7.40\pm0.03$ days for precession of the long axis about the angular momentum vector based on Fourier analysis. These results provide a further understanding of rotational state of Toutatis.

The occurrence of low-amplitude flux variations in blazars on hourly timescales, commonly known as microvariability, is still a widely debated subject in high-energy astrophysics. Several competing scenarios have been proposed to explain such occurrences, including various jet plasma instabilities leading to the formation of shocks, magnetic reconnection sites, and turbulence. In this letter we present the results of our detailed investigation of a prominent, five-hour-long optical microflare detected during recent WEBT campaign in 2014, March 2-6 targeting the blazar 0716+714. After separating the flaring component from the underlying base emission continuum of the blazar, we find that the microflare is highly polarized, with the polarization degree $\sim (40-60)\%$$\pm (2-10)\%$, and the electric vector position angle $\sim (10 - 20)$deg$\pm (1-8)$deg slightly misaligned with respect to the position angle of the radio jet. The microflare evolution in the $(Q,\,U)$ Stokes parameter space exhibits a looping behavior with a counter-clockwise rotation, meaning polarization degree decreasing with the flux (but higher in the flux decaying phase), and approximately stable polarization angle. The overall very high polarization degree of the flare, its symmetric flux rise and decay profiles, and also its structured evolution in the $Q-U$ plane, all imply that the observed flux variation corresponds to a single emission region characterized by a highly ordered magnetic field. As discussed in the paper, a small-scale but strong shock propagating within the outflow, and compressing a disordered magnetic field component, provides a natural, though not unique, interpretation of our findings.

Planets are known to orbit giant stars, yet there is a shortage of planets orbiting within ~0.5 AU (P<100 days). First-ascent giants have not expanded enough to engulf such planets, but tidal forces can bring planets to the surface of the star far beyond the stellar radius. So the question remains: are tidal forces strong enough in these stars to engulf all the missing planets? We describe a high-cadence observational program to obtain precise radial velocities of bright giants from Weihai Observatory of Shandong University. We present data on the planet host Beta Gem (HD 62509), confirming our ability to derive accurate and precise velocities; our data achieve an rms of 7.3 m/s about the Keplerian orbit fit. This planet-search programme currently receives ~100 nights per year, allowing us to aggressively pursue short-period planets to determine whether they are truly absent.

In this work, we investigate the rotational dynamics of the ginger-shaped near-Earth asteroid 4179 Toutatis, which was closely observed by Chang'e-2 at a distance of $770\pm120~$ meters from the asteroid's surface during the outbound flyby \citep{Huang2013} on 13 December 2012. A sequence of high-resolution images was acquired during the flyby mission. In combination with ground-based radar observations collected over the last two decades, we analyze these flyby images and determine the orientation of the asteroid at the flyby epoch. The 3-1-3 Euler angles of the conversion matrix from the J2000 ecliptic coordinate system to the body-fixed frame are evaluated to be $-20.1^\circ\pm1^\circ$, $27.6^\circ\pm1^\circ$ and $42.2^\circ\pm1^\circ$, respectively. The least-squares method is utilized to determine the rotational parameters and spin state of Toutatis. The characteristics of the spin-state parameters and angular momentum variations are extensively studied using numerical simulations, which confirm those reported by \citet{Takahashi2013}. The large amplitude of Toutatis' precession is assumed to be responsible for its tumbling attitude as observed from Earth. Toutatis' angular momentum orientation is determined to be described by $\lambda_{H}=180.2^{+0.2^\circ}_{-0.3^\circ}$ and $\beta_{H}=-54.75^{+0.15^\circ}_{-0.10^\circ}$, implying that it has remained nearly unchanged for two decades. Furthermore, using Fourier analysis to explore the change in the orientation of Toutatis' axes, we reveal that the two rotational periods are 5.38 and 7.40 days, respectively, consistent with the results of the former investigation. Hence, our investigation provides a clear understanding of the state of the rotational dynamics of Toutatis.

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.

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 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.

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

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.

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.

We report our investigation of the first transiting planet candidate from the YETI project in the young (~4 Myr old) open cluster Trumpler 37. The transit-like signal detected in the lightcurve of the F8V star 2M21385603+5711345 repeats every 1.364894+/-0.000015 days, and has a depth of 54.5+/-0.8 mmag in R. Membership to the cluster is supported by its mean radial velocity and location in the color-magnitude diagram, while the Li diagnostic and proper motion are inconclusive in this regard. Follow-up photometric monitoring and adaptive optics imaging allow us to rule out many possible blend scenarios, but our radial-velocity measurements show it to be an eclipsing single-lined spectroscopic binary with a late-type (mid-M) stellar companion, rather than one of planetary nature. The estimated mass of the companion is 0.15-0.44 solar masses. The search for planets around very young stars such as those targeted by the YETI survey remains of critical importance to understand the early stages of planet formation and evolution.

On 13 December 2012, Chang'e-2 conducted a successful flyby of the near-Earth asteroid 4179 Toutatis at a closest distance of 770 $\pm$ 120 meters from the asteroid's surface. The highest-resolution image, with a resolution of better than 3 meters, reveals new discoveries on the asteroid, e.g., a giant basin at the big end, a sharply perpendicular silhouette near the neck region, and direct evidence of boulders and regolith, which suggests that Toutatis may bear a rubble-pile structure. Toutatis' maximum physical length and width are (4.75 $\times$ 1.95 km) $\pm$10$\%$, respectively, and the direction of the +$z$ axis is estimated to be (250$\pm$5$^\circ$, 63$\pm$5$^\circ$) with respect to the J2000 ecliptic coordinate system. The bifurcated configuration is indicative of a contact binary origin for Toutatis, which is composed of two lobes (head and body). Chang'e-2 observations have significantly improved our understanding of the characteristics, formation, and evolution of asteroids in general.

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.

Context. The international whole earth blazar telescope (WEBT) consortium planned and carried out three days of intensive micro-variability observations of S5 0716+714 from February 22, 2009 to February 25, 2009. This object was chosen due to its bright apparent magnitude range, its high declination, and its very large duty cycle for micro-variations. Aims. We report here on the long continuous optical micro-variability light curve of 0716+714 obtained during the multi-site observing campaign during which the Blazar showed almost constant variability over a 0.5 magnitude range. The resulting light curve is presented here for the first time. Observations from participating observatories were corrected for instrumental differences and combined to construct the overall smoothed light curve. Methods. Thirty-six observatories in sixteen countries participated in this continuous monitoring program and twenty of them submitted data for compilation into a continuous light curve. The light curve was analyzed using several techniques including Fourier transform, Wavelet and noise analysis techniques. Those results led us to model the light curve by attributing the variations to a series of synchrotron pulses. Results. We have interpreted the observed microvariations in this extended light curve in terms of a new model consisting of individual stochastic pulses due to cells in a turbulent jet which are energized by a passing shock and cool by means of synchrotron emission. We obtained an excellent fit to the 72-hour light curve with the synchrotron pulse model.

3C 66A was monitored by the BATC(Beijing-Arizona-Taipei-Connecticut) telescope from 2005 to 2008,1994 observations were obtained on 89 nights. Detailed research and analysis was performed on these observations in this paper. A long term burst occurred in the whole light curve. No intra-day variability was claimed in our campaign by intra-night light curve analysis. Time lag of shorter wavelenth preceding longer wavelength was shown by correlation analysis. The results showed that the optical spectral shape turned flatter when the source brightened, and the spectral variability indicator was bigger on shorter time-scale as determined by the color indices variation analysis.

Long-duration gamma-ray bursts (GRBs) at z < 1 are in most cases found to be accompanied by bright, broad-lined Type Ic supernovae (SNe Ic-BL). The highest-energy GRBs are mostly located at higher redshifts, where the associated SNe are hard to detect observationally. Here we present early and late observations of the optical counterpart of the very energetic GRB 130427A. Despite its moderate redshift z = 0.3399+/-0.0002, GRB 130427A is at the high end of the GRB energy distribution, with an isotropic-equivalent energy release of Eiso ~ 9.6x10^53 erg, more than an order of magnitude more energetic than other GRBs with spectroscopically confirmed SNe. In our dense photometric monitoring, we detect excess flux in the host-subtracted r-band light curve, consistent with what expected from an emerging SN, ~0.2 mag fainter than the prototypical SN 1998bw. A spectrum obtained around the time of the SN peak (16.7 days after the GRB) reveals broad undulations typical of SNe Ic-BL, confirming the presence of a SN, designated SN 2013cq. The spectral shape and early peak time are similar to those of the high expansion velocity SN 2010bh associated with GRB 100316D. Our findings demonstrate that high-energy long-duration GRBs, commonly detected at high redshift, can also be associated with SNe Ic-BL, pointing to a common progenitor mechanism.

With an apparent cluster diameter of 1.5{\deg} and an age of ~4 Myr, Trumpler 37 is an ideal target for photometric monitoring of young stars as well as for the search of planetary transits, eclipsing binaries and other sources of variability. The YETI consortium has monitored Trumpler 37 throughout 2010 and 2011 to obtain a comprehensive view of variable phenomena in this region. In this first paper we present the cluster properties and membership determination as derived from an extensive investigation of the literature. We also compared the coordinate list to some YETI images. For 1872 stars we found literature data. Among them 774 have high probability of being member and 125 a medium probability. Based on infrared data we re-calculate a cluster extinction of 0.9-1.2 mag. We can confirm the age and distance to be 3-5 Myr and ~870 pc. Stellar masses are determined from theoretical models and the mass function is fitted with a power-law index of alpha=1.90 (0.1-0.4 M_sun) and alpha=1.12 (1-10 M_sun).

GM Cep in the young (~4 Myr) open cluster Trumpler 37 has been known to be an abrupt variable and to have a circumstellar disk with very active accretion. Our monitoring observations in 2009-2011 revealed the star to show sporadic flare events, each with brightening of < 0.5 mag lasting for days. These brightening events, associated with a color change toward the blue, should originate from an increased accretion activity. Moreover, the star also underwent a brightness drop of ~1 mag lasting for about a month, during which the star became bluer when fainter. Such brightness drops seem to have a recurrence time scale of a year, as evidenced in our data and the photometric behavior of GM Cep over a century. Between consecutive drops, the star brightened gradually by about 1 mag and became blue at peak luminosity. We propose that the drop is caused by obscuration of the central star by an orbiting dust concentration. The UX Orionis type of activity in GM Cep therefore exemplifies the disk inhomogeneity process in transition between grain coagulation and planetesimal formation in a young circumstellar disk.

We present our photometric observations of GM Cep, a solar type variable in the young (~4 Myr) open cluster Trumpler 37. The star is known as a classical T Tauri star with a circumstellar disk and active accretion. GM\,Cep was suspected to undergo an outburst, thus a candidate for an EXor-type variable. In our monitoring campaign observations in 2010-2011, GM Cep experienced a ~0.82 mag brightness decrease in the R band lasting for 39 days, and frequent, transient flare-like episodes with amplitude < 1 mag, each lasting for about 10 days. The brightening was accompanied with a bluer color, presumably arising from increased accretion activity. Interestingly, the star also turned bluer in the fading phase. Combining the AAVSO and literature data, we found a quasi-cyclic peroid of ~311 days for the fading event. A possible mechanism for the fading could be obscuration by a clump of dust around the star. We proposed that GM Cep therefore should be a UXor-type variable in the transition phase between grain coagulation and planetesimal formation process in the circumstellar disk.

We present the Young Exoplanet Transit Initiative (YETI), in which we use several 0.2 to 2.6m telescopes around the world to monitor continuously young (< 100 Myr), nearby (< 1 kpc) stellar clusters mainly to detect young transiting planets (and to study other variability phenomena on time-scales from minutes to years). The telescope network enables us to observe the targets continuously for several days in order not to miss any transit. The runs are typically one to two weeks long, about three runs per year per cluster in two or three subsequent years for about ten clusters. There are thousands of stars detectable in each field with several hundred known cluster members, e.g. in the first cluster observed, Tr-37, a typical cluster for the YETI survey, there are at least 469 known young stars detected in YETI data down to R=16.5 mag with sufficient precision of 50 milli-mag rms (5 mmag rms down to R=14.5 mag) to detect transits, so that we can expect at least about one young transiting object in this cluster. If we observe 10 similar clusters, we can expect to detect approximately 10 young transiting planets with radius determinations. The precision given above is for a typical telescope of the YETI network, namely the 60/90-cm Jena telescope (similar brightness limit, namely within +/-1 mag, for the others) so that planetary transits can be detected. For planets with mass and radius determinations, we can calculate the mean density and probe the internal structure. We aim to constrain planet formation models and their time-scales by discovering planets younger than 100 Myr and determining not only their orbital parameters, but also measuring their true masses and radii, which is possible so far only by the transit method. Here, we present an overview and first results. (Abstract shortened)

The aim of our observations is to investigate the intranight variability properties and the spectral variability of BL Lacertae. 799 optical multi-band observations were intensively made with the Beijing-Arizona-Taiwan-Connecticut (BATC) 60/90cm Schmidt telescope during the outburst composed of two subsequent flares in 2005. The second flare, whose rising phase lasted at least 44 days, was observed with amplitudes of more than 1.1 mag in three BATC optical bands. In general, on intranight timescale the amplitude of variability and the variation rate are larger at the shorter wavelength, and the variation rate is comparable in the rising and decaying phases within each band. A possible time lag between the light curves in bands e and m, around 11.6 minutes, was obtained. Based on the analysis of the colour index variation with the source brightness, the variability of BL Lacertae can be interpreted as having two components: a "strongly-chromatic" intranight component and a "mildly-chromatic" internight component, which may be the results of both intrinsic physical mechanism and geometric effects.

Many quasi-simultaneous optical observations of 17 blazars are obtained from previous papers published over the last 19 years in order to investigate the spectral slope variability and understand the radiation mechanism of blazars. The long-period dereddened optical spectral slopes are calculated. We analyse the average spectral slope distribution, which suggests that the spectra of flat spectrum radio quasars (FSRQs) and high energy peaked BL Lac objects (HBLs) are probably deformed by other emission components. The average spectral slopes of low energy peaked BL Lac objects(LBLs), which scatter around 1.5, show a good accordance with the synchrotron self-Compton (SSC) loss-dominated model. We present and discuss the variability between the spectral slope and optical luminosity. The spectra of all HBLs and LBLs get flatter when they turn brighter, while for FSRQs this trend does not exist or may even be reversed. This phenomenon may imply that there is a thermal contribution to the optical spectrum for FSRQs. For the FSRQ 1156+295, there is a hint that the slope gets flatter at both the brightest and faintest states. Our result shows that three subclasses locate in different regions in the pattern of slope variability indicator versus average spectral slope. The relativistic jet mechanism is supported by the significant correlation between the optical Doppler factor and the average spectral slope.