The Milky Way (MW) mass $M_{\rm vir}$ is a fundamental quantity in astronomy. Although it has been measured extensively, it is still uncertain to more than a factor of two due to limited number or spatial coverage of kinematic tracers. Here we use a novel method to constrain $M_{\rm vir}$ based on the properties of the MW corona. We build a hydrostatic corona model with non-thermal pressure support and a physically-motivated density profile, and derive the temperature distribution, which depends on $M_{\rm vir}$. While the temperature profile decreases substantially with radius, the X-ray-emission-weighted average temperature $T_{\rm em}$ is quite uniform toward different sight lines, consistent with X-ray observations. Using available measurements of $T_{\rm em}$, we find that $M_{\rm vir}$=$1.60^{+1.35}_{-0.41}$ $\times$ $10^{12}$ $M_{\rm sun}$ assuming an Navarro-Frenk-White (NFW) total matter distribution. This estimate is independent of the uncertain total corona mass and the gas temperature at very large radii, and is on the high mass side of current measurements. Adopting a $\beta$ model for the corona density distribution or a total matter distribution contributed by an NFW dark matter distribution and a central cold baryonic matter distribution leads to similar estimates of $M_{\rm vir}$. Non-thermal pressure support, which likely exists in the corona, leads to higher values of $M_{\rm vir}$.

The Fermi bubbles are two giant bubbles in gamma rays lying above and below the Galactic center (GC). Despite numerous studies on the bubbles, their origin and emission mechanism remain elusive. Here we use a suite of hydrodynamic simulations to study the scenario where the cosmic rays (CRs) in the bubbles are mainly accelerated at the forward shocks driven by a pair of opposing jets from Sgr A*. We find that an active galactic nucleus (AGN) jet event happened $5-6$ Myr ago can naturally reproduce the bilobular morphology of the bubbles, and the postshock gas temperature in the bubbles is heated to $\sim0.4$ keV, consistent with recent X-ray observations. The forward shocks compress the hot halo gas, and at low latitudes, the compressed gas shows an X-shaped structure, naturally explaining the biconical X-ray structure in the ROSAT 1.5 keV map in both morphology and X-ray surface brightness. CR acceleration is most efficient in the head regions of the bubbles during the first 2 Myrs. The opposing jets release a total energy of $\sim 10^{55}$ erg with an Eddington ratio of $\sim 10^{-3}$, which falls well in the range of the hot accretion flow mode for black holes. Our simulations further show that the forward shocks driven by spherical winds at the GC typically produce bubbles with much wider bases than observed, and could not reproduce the biconical X-ray structure at low latitudes. This suggests that starburst or AGN winds are unlikely the origin of the bubbles in the shock scenario.

The second observational campaign of gravitational waves organized by the LIGO/Virgo Collaborations has led to several breakthroughs such as the detection of gravitational wave signals from merger systems involving black holes or neutrons stars. During O2,14 gravitational wave alerts were sent to the astronomical community with sky regions covering mostly over hundreds of square degrees. Among them, 6 have been finally confirmed as real astrophysical events. Since 2013, a new set of ground-based robotic telescopes called GWAC and its pathfinder mini-GWAC have been developed to contribute to the various challenges of themulti-messenger and time domain astronomy. The GWAC system is built up in the framework of the ground-segment system of the SVOM mission that will be devoted to the study of the multi-wavelength transient sky in the next decade. During O2, only the mini-GWAC telescopenetwork was fully operational. Due to the wide field of view and fast automatic follow-up capabilities of the mini-GWAC telescopes, they were well adapted to efficiently cover the sky localization areas of the gravitational wave event candidates. In this paper, we present the mini-GWAC pipeline we have set up to respond to the GW alerts and we report our optical follow-up observations of 8 GW alerts detected during the O2 run. Our observations provided the largest coverage of the GW localization areas in a short latency made by any optical facility. We found tens of optical transient candidates in our images, but none of those could be securely associated with any confirmed black hole-black hole merger event. Based on this first experience and the near future technical improvements of our network system, we will be more competitive to detect the optical counterparts from some gravitational wave events that will be detected during the upcoming O3 run, especially those emerging from binary neutron star mergers.

Lijiang 2.4-meter Telescope(LJT), the largest common-purpose optical telescope in China, has been applied to the world-wide astronomers since 2008. It is located at Gaomeigu site, Lijiang Observatory(LJO), the southwest of China. The site has very good observational conditions. Since 10-year operation, several instruments have been equipped on the LJT. Astronomers can perform both photometric and spectral observations. The main scientific goals of LJT include photometric and spectral evolution of supernova, reverberation mapping of active galactic nucleus, physical properties of binary star and near-earth object(comet and asteroid), identification of exoplanet, and all kinds of transients. Until now, the masses of 41 high accretion rate black holes have been measured, and more than 168 supernova have been identified by the LJT. More than 190 papers related to the LJT have been published. In this paper, the general observation condition of the Gaomeigu site is introduced at first. Then, the LJT structure is described in detail, including the optical, mechanical, motion and control system. The specification of all the instruments, and some detailed parameters of the YFOSC is also presented. Finally, some important scientific results and future expectations are summarized.

The Earth observation satellites (EOSs) scheduling is of great importance to achieve efficient observation missions. The agile EOSs (AEOS) with stronger attitude maneuvering capacity can greatly improve observation efficiency while increasing scheduling complexity. The multiple AEOSs, oversubscribed targets scheduling problem with multiple observations are addressed, and the potential observation missions are modeled as nodes in the complex networks. To solve the problem, an improved feedback structured heuristic is designed by defining the node and target importance factors. On the basis of a real world Chinese AEOS constellation, simulation experiments are conducted to validate the heuristic efficiency in comparison with a constructive algorithm and a structured genetic algorithm.

We examine the effects that dynamical instability has on shaping the orbital properties of exoplanetary systems. Using N-body simulations of non-EMS (Equal Mutual Separation), multi-planet systems we find that the lower limit of the instability timescale $t$ is determined by the minimal mutual separation $K_{\rm min}$ in units of the mutual Hill radius. Planetary systems showing instability generally include planet pairs with period ratio $<1.33$. Our final period ratio distribution of all adjacent planet pairs shows dip-peak structures near first-order mean motion resonances similar to those observed in the \kepler\ planetary data. Then we compare the probability density function (PDF) of the de-biased \kepler\ period ratios with those in our simulations and find a lack of planet pairs with period ratio $> 2.1$ in the observations---possibly caused either by inward migration before the dissipation of the disk or by planet pairs not forming with period ratios $> 2.1$ with the same frequency they do with smaller period ratios. By comparing the PDF of the period ratio between simulation and observation, we obtain an upper limit of 0.03 on the scale parameter of the Rayleigh distributed eccentricities when the gas disk dissipated. Finally, our results suggest that a viable definition for a `packed' or `compact' planetary system be one that has at least one planet pair with a period ratio less than 1.33. This criterion would imply that 4\% of the \kepler\ systems (or 6\% of the systems with more than two planets) are compact.

The LAMOST-Kepler (LK-) project was initiated to use the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) to make spectroscopic follow-up observations for the targets in the field of the Kepler mission. The Kepler field is divided into 14 subfields that are adapted to the LAMOST circular field with a diameter of 5 degrees. During the regular survey phase of LAMOST, the LK-project took data from 2012 June to 2017 June and covered all the 14 subfields at least twice. In particular, we describe in this paper the second Data Release of the LK-project, including all spectra acquired through 2015 May to 2017 June together with the first round observations of the LK-project from 2012 June to 2014 September. The LK-project now counts 227 870 spectra of 156 390 stars, among which we have derived atmospheric parameters (log g, T eff and [Fe/H]) and heliocentric radial velocity (RV) for 173 971 spectra of 126 172 stars. These parameters were obtained with the most recent version of the LAMOST Stellar Parameter Pipeline v 2.9.7. Nearly one half, namely 76 283 targets, are observed both by LAMOST and Kepler telescopes. These spectra, establishing a large spectroscopy library, will be useful for the entire astronomical community, particularly for planetary science and stellar variability on Kepler targets.

We present our photometric monitoring of a flat spectrum radio quasar (FSRQ) 3C 454.3 at Yunnan observatories from 2006 to 2011. We find that the optical color of 3C 454.3 shows obvious redder-when-brighter trend, which reaches a saturation stage when the source is brighter than 15.15 mag at V band. We perform a simulation with multiple values of disk luminosity and spectral index to reproduce the magnitude-color diagram. The results show that the contamination caused by the disk radiation alone is difficult to produce the observed color variability. The variability properties during the outburst in December 2009 are also compared with $\gamma$-ray data derived from Fermi $\gamma$-ray space telescope. The flux variation of these two bands follow a linear relation with $F_{\gamma} \propto F_R^{1.14\pm0.07}$, which provides an observational evidence for external Compton process in 3C 454.3. Meanwhile, this flux correlation indicates that electron injection is the main mechanism for variability origin. We also explore the variation of the flux ratio $F_{\gamma}/F_R$ and the detailed structures in the lightcurves, and discuss some possible origins for the detailed variability behaviors.

All of the 14 subfields of the Kepler field have been observed at least once with the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST, Xinglong Observatory, China) during the 2012-2014 observation seasons. There are 88,628 reduced spectra with SNR$_g$ (signal-to-noise ratio in g band) $\geq$ 6 after the first round (2012-2014) of observations for the LAMOST-Kepler project (LK-project). By adopting the upgraded version of the LAMOST Stellar Parameter pipeline (LASP), we have determined the atmospheric parameters ($T_{\rm eff}$ , $\log g$, and $\rm [Fe/H]$) and heliocentric radial velocity $v_{\rm rad}$ for 51,406 stars with 61,226 spectra. Compared with atmospheric parameters derived from both high-resolution spectroscopy and asteroseismology method for common stars in Huber et al. (2014), an external calibration of LASP atmospheric parameters was made, leading to the determination of external errors for the giants and dwarfs, respectively. Multiple spectroscopic observations for the same objects of the LK-project were used to estimate the internal uncertainties of the atmospheric parameters as a function of SNR$_g$ with the unbiased estimation method. The LASP atmospheric parameters were calibrated based on both the external and internal uncertainties for the giants and dwarfs, respectively. A general statistical analysis of the stellar parameters leads to discovery of 106 candidate metal-poor stars, 9 candidate very metal-poor stars, and 18 candidate high-velocity stars. Fitting formulae were obtained segmentally for both the calibrated atmospheric parameters of the LK-project and the KIC parameters with the common stars. The calibrated atmospheric parameters and radial velocities of the LK-project will be useful for studying stars in the Kepler field.

The Chinese Small Telescope ARray (hereinafter CSTAR) is the first Chinese astronomical instrument on the Antarctic ice cap. The low temperature and low pressure testing of the data acquisition system was carried out in a laboratory refrigerator and on the 4500m Pamirs high plateau, respectively. The results from the final four nights of test observations demonstrated that CSTAR was ready for operation at Dome A, Antarctica. In this paper we present a description of CSTAR and the performance derived from the test observations.

We reassess the evidence that WMAP temperature maps contain a statistically significant "cold spot" by repeating the analysis using simple circular top-hat (disk) weights, as well as Gaussian weights of varying width. Contrary to previous results that used Spherical Mexican Hat Wavelets, we find no significant signal at any scale when we compare the coldest spot from our sky to ones from simulated Gaussian random, isotropic maps. We trace this apparent discrepancy to the fact that WMAP cold spot's temperature profile just happens to favor the particular profile given by the wavelet. Since randomly generated maps typically do not exhibit this coincidence, we conclude that the original cold spot significance originated at least partly due to a fortuitous choice of using a particular basis of weight functions. We also examine significance of a more general measure that returns the most significant result among several choices of the weighting function, angular scale of the spot, and the statistics applied, and again find a null result.

Recent fits of cosmological parameters by the Wilkinson Microwave Anisotropy Probe (WMAP) measurement favor a primordial scalar spectrum with varying index. This result, if stands, could severely constrain inflation model buildings. Most extant slow-roll inflation models allow for only a tiny amount of scale variations in the spectrum. We propose in this paper an extra-dimensional inflation model which is natural theoretically and can generate the required variations of the spectral index as implied by the WMAP for suitable choices of parameters.

The Friedmann law on the brane generically depends quadratically on the brane energy density and involves a ``dark radiation'' term due to the bulk Weyl tensor. Despite its unfamiliar form, we show how it can be derived from a standard four-dimensional Brans-Dicke theory at low energy. In particular, the dark radiation term is found to depend linearly on the brane energy densities. For any equation of state on the branes, the radion evolves such as to generate radiation-dominated cosmology. The radiation-dominated era is conventional and consistent with nucleosynthesis.

We analyse the cosmological implications of brane-antibrane systems in string-theoretic orbifold and orientifold models. In a class of realistic models, consistency conditions require branes and antibranes to be stuck at different fixed points, and so their mutual attraction generates a potential for one of the radii of the underlying torus or the 4D string dilaton. Assuming that all other moduli have been fixed by string effects, we find that this potential leads naturally to a period of cosmic inflation with the radion or dilaton field as the inflaton. The slow-roll conditions are satisfied more generically than if the branes were free to move within the space. The appearance of tachyon fields at certain points in moduli space indicates the onset of phase transitions to different non-BPS brane systems, providing ways of ending inflation and reheating the corresponding observable brane universe. In each case we find relations between the inflationary parameters and the string scale to get the correct spectrum of density perturbations. In some examples the small numbers required as inputs are no smaller than 0.01, and are the same small quantities which are required to explain the gauge hierarchy.

We show how the motion through the extra dimensions of a gas of branes and antibranes can, under certain circumstances, produce an era of inflation as seen by observers trapped on a 3-brane, with the inflaton being the inter-brane separation. Although most of our discussion refers to arbitrary p-branes, when we need to be specific we assume that they are D-branes of Type II or Type I string theory. For realistic brane couplings, such as those arising in string theory, the inter-brane potentials are too steep to inflate the universe for acceptably long times. However, for special regions of the parameter space of brane-antibrane positions the brane motion is slow enough for there to be sufficient inflation. Inflation would be more generic in models where the inter-brane interactions are much weaker. The spectrum of primordial density fluctuations predicted has index n slightly less than 1, and an acceptable amplitude, provided that the extra dimensions have linear size 1/r ~ 10^{12} GeV. Reheating occurs as in hybrid inflation, with the tachyonic instability of the brane-antibrane system taking over for small separations. The tachyon field can induce a cascade mechanism within which higher-dimension branes annihilate into lower-dimension ones. We argue that such a cascade naturally stops with the production of 3-branes in 10-dimensional string theory.

In the Kaluza-Klein (KK) scenario with n large extra dimensions where gravity propagates in the 4+n dimensional bulk of spacetime while gauge and matter fields are confined to a four dimensional subspace, the light graviton KK modes can be produced in the Sun, red giants and supernovae. We study the energy-loss rates through photon-photon annihilation, electron-positron annihilation, gravi-Compton-Primakoff scattering, gravi-bremsstrahlung and nucleon-nucleon bremsstrahlung, and derive lower limits to the string scale M_S. The most stringent lower limit obtained from SN1987A leads to $M_S> 30 - 130$ TeV (2.1-9.2 TeV) for the case of two (three) large extra dimensions.