We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~ 1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multi-scale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~ 1 pc and remains continuous from the scales of filaments (~ 0.1 pc) to that of protostellar envelopes (~ 0.005 pc or ~ 1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data towards the filament network.

In LTE, the two inner satellite lines (ISLs) and the two outer satellite lines (OSLs) of the NH$_{3}$ (1,1) transition are each predicted to have equal intensities. However, hyperfine intensity anomalies (HIAs) are observed to be omnipresent in star formation regions, which is still not fully understood. In addressing this issue, we find that the computation method of the HIA by the ratio of the peak intensities may have defects, especially when being used to process the spectra with low velocity dispersions. Therefore we define the integrated HIAs of the ISLs (HIA$_{\rm IS}$) and OSLs (HIA$_{\rm OS}$) by the ratio of their redshifted to blueshifted integrated intensities and develop a procedure to calculate them. Based on this procedure, we present a systematic study of the integrated HIAs in the northern part of the Orion A MC. We find that integrated HIA$_{\rm IS}$ and HIA$_{\rm OS}$ are commonly present in the Orion A MC and no clear distinction is found at different locations of the MC. The medians of the integrated HIA$_{\rm IS}$ and HIA$_{\rm OS}$ are 0.921$\pm$0.003 and 1.422$\pm$0.009, respectively, which is consistent with the HIA core model and inconsistent with the CE model. Selecting those 170 positions where both integrated HIAs deviate by more than 3-$\sigma$ from unity, most (166) are characterized by HIA$_{\rm IS}$<1 and HIA$_{\rm OS}$>1, which suggests that the HIA core model plays a more significant role than the CE model. The remaining four positions are consistent with the CE model. We compare the integrated HIAs with the para-NH$_{3}$ column density ($N$(para-NH$_{3}$)), kinetic temperature ($T_{\rm K}$), total velocity dispersion ($\sigma_{\rm v}$), non-thermal velocity dispersion ($\sigma_{\rm NT}$), and the total opacity of the NH$_{3}$ (1,1) line ($\tau_{0}$). Their correlations can not be fully explained by neither the HIA core nor the CE model.

Fast radio bursts (FRBs) are short pulses observed in radio band from cosmological distances, some of which emit repeating bursts. The physical origins of these mysterious events have been subject to wide speculations and heated debates. One class of models invoke soft gamma-ray repeaters (SGRs), or magnetars, as the sources of FRBs. Magnetars are rotating neutron stars with extremely strong magnetic field and can sporadically emit bursts from X-ray ($\sim$keV) to soft gamma-ray ($\sim$sub-MeV) with duration from $10^{-2}$~s to $10^2$~s. However, even though some bright radio bursts have been observed from some magnetars, no FRB-like events had been detected to be associated with any magnetar burst, including one giant flare, and no radio burst has been associated with any X-ray event from any magnetar. Therefore, there is still no observational evidence for magnetar-FRB association up to today. Recently, a pair of FRB-like bursts (FRB~200428 hereafter) separated by 30 milliseconds (ms) were detected from the general direction of the Galactic magnetar SGR~J1935+2154. Here we report the detection of a non-thermal X-ray burst in the 1-250 keV energy band with the Insight-HXMT satellite, which we identify as emitted from SGR~J1935+2154. The burst showed two hard peaks with a separation of $\sim30$ ms, consistent with the separation between the two bursts in FRB~200428. The delay time between the double radio and X-ray peaks is $\sim8.57$ s, fully consistent with the dispersion delay of FRB~200428. We thus identify the non-thermal X-ray burst is associated with FRB~200428 whose high energy counterpart is the two hard peaks in X-ray. Our results suggest that the non-thermal X-ray burst and FRB~200428 share the same physical origin in an explosive event from SGR~J1935+2154.

We use the In data collected during the 2019 outburst from X-ray pulsar 4U 1901+03 to complement the orbital parameters reported by Fermi/GBM. Using the Insight-HXMT, we examine the correlation between the derivative of the intrinsic spin frequency and bolometric flux based on accretion torque models. It was found that the pulse profiles significantly evolve during the outburst. The existence of two types of the profile's pattern discovered in the Insight-HXMT data indicates that this source experienced transition between a super-critical and a sub-critical accretion regime during its 2019 outburst. Based on the evolution of the pulse profiles and the torque model, we derive the distance to 4U 1901+03 as 12.4+-0.2 kpc.

We study the evolution of the temporal properties of MAXI 1820+070 during the 2018 outburst in its hard state from MJD 58190 to 58289 with Insight-HXMT in a broad energy band 1-150 keV. We find different behaviors of the hardness ratio, the fractional rms and time lag before and after MJD 58257, suggesting a transition occurred at around this point. The observed time lags between the soft photons in the 1-5 keV band and the hard photons in higher energy bands, up to 150 keV, are frequency-dependent: the time lags in the low-frequency range, 2-10 mHz, are both soft and hard lags with a timescale of dozens of seconds but without a clear trend along the outburst; the time lags in the high-frequency range, 1-10 Hz, are only hard lags with a timescale of tens of milliseconds; first increase until around MJD 58257 and decrease after this date. The high-frequency time lags are significantly correlated to the photon index derived from the fit to the quasi-simultaneous NICER spectrum in the 1-10 keV band. This result is qualitatively consistent with a model in which the high-frequency time lags are produced by Comptonization in a jet.

Glitches correspond to sudden jumps of rotation frequency ($\nu$) and its derivative ($\dot{\nu}$) of pulsars, the origin of which remains not well understood yet, partly because the jump processes of most glitches are not well time-resolved. There are three large glitches of the Crab pulsar, detected in 1989, 1996 and 2017, which were found to have delayed spin-up processes before the normal recovery processes. Here we report two additional glitches of the Crab pulsar occurred in 2004 and 2011 for which we discovered delayed spin up processes, and present refined parameters of the largest glitch occurred in 2017. The initial rising time of the glitch is determined as $<0.48$ hour. We also carried out a statistical study of these five glitches with observed spin-up processes. The two glitches occurred in 2004 and 2011 have delayed spin-up time scales ($\tau_{1}$) of $1.7\pm0.8$\,days and $1.6\pm0.4$\,days, respectively. We find that the $\Delta{\nu}$ vs. $|\Delta{\dot\nu}|$ relation of these five glitches is similar to those with no detected delayed spin-up process, indicating that they are similar to the others in nature except that they have larger amplitudes. For these five glitches, the amplitudes of the delayed spin-up process ($|\Delta{\nu}_{\rm d1}|$) and recovery process ($\Delta{\nu}_{\rm d2}$), their time scales ($\tau_{1}$, $\tau_{2}$), and permanent changes in spin frequency ($\Delta{\nu}_{\rm p}$) and total frequency step ($\Delta{\nu}_{\rm g}$) have positive correlations. From these correlations, we suggest that the delayed spin-up processes are common for all glitches, but are too short and thus difficult to be detected for most glitches.

The Kepler mission has discovered thousands of exoplanets around various stars with different spectral types (M, K, G, and F) and thus different masses and effective temperatures. Previous studies have shown that the planet occurrence rate, in terms of average number of planets per star, drops with increasing stellar effective temperature (Teff). In this paper, with the final Kepler Data Release (DR25) catalog, we revisit the relation between stellar effective temperature (as well as mass) and planet occurrence, but in terms of the fraction of stars with planets and the number of planets per planetary system (i.e. planet multiplicity). We find that both the fraction of stars with planets and planet multiplicity decrease with increasing stellar temperature and mass. Specifically, about 75% late-type stars (Teff<5000 K) have Kepler-like planets with an average planet multiplicity of ~2.8; while for early-type stars (Teff>6500 K), this fraction and the average multiplicity fall down to ~35% and ~1.8, respectively. The decreasing trend in the fraction of stars with planets is very significant with ${\Delta}$AIC> 30, though the trend in planet multiplicity is somewhat tentative with ${\Delta}$AIC~5. Our results also allow us to derive the dispersion of planetary orbital inclinations in relationship with stellar effective temperature. Interestingly, it is found to be similar to the well-known trend between obliquity and stellar temperature, indicating that the two trends might have a common origin.

No robust detection of prompt electromagnetic counterparts to fast radio bursts (FRBs) has yet been obtained, in spite of several multi-wavelength searches carried out so far. Specifically, X/gamma-ray counterparts are predicted by some models. We planned on searching for prompt gamma-ray counterparts in the Insight-Hard X-ray Modulation Telescope (Insight-HXMT) data, taking advantage of the unique combination of large effective area in the keV-MeV energy range and of sub-ms time resolution. We selected 39 FRBs that were promptly visible from the High-Energy (HE) instrument aboard Insight-HXMT. After calculating the expected arrival times at the location of the spacecraft, we searched for a significant excess in both individual and cumulative time profiles over a wide range of time resolutions, from several seconds down to sub-ms scales. Using the dispersion measures in excess of the Galactic terms, we estimated the upper limits on the redshifts. No convincing signal was found and for each FRB we constrained the gamma-ray isotropic-equivalent luminosity and the released energy as a function of emission timescale. For the nearest FRB source, the periodic repeater FRB180916.J0158+65, we find $L_{\gamma,iso}<5.5\times 10^{47}$ erg/s over 1 s, whereas $L_{\gamma,iso}<10^{49}-10^{51}$ erg/s for the bulk of FRBs. The same values scale up by a factor of ~100 for a ms-long emission. Even on a timescale comparable with that of the radio pulse itself no keV-MeV emission is observed. A systematic association with either long or short GRBs is ruled out with high confidence, except for subluminous events, as is the case for core-collapse of massive stars (long) or binary neutron star mergers (short) viewed off axis. Only giant flares from extra-galactic magnetars at least ten times more energetic than Galactic siblings are ruled out for the nearest FRB.

We report on our analysis of the 2019 outburst of the X-ray accreting pulsar 4U 1901+03 observed with Insight-HXMT and NICER. Both spectra and pulse profiles evolve significantly in the decaying phase of the outburst. Dozens of flares are observed throughout the outburst. They are more frequent and brighter at the outburst peak. We find that the flares, which have a duration from tens to hundreds of seconds, are generally brighter than the persistent emission by a factor of $\sim$ 1.5. The pulse profile shape during the flares can be significantly different than that of the persistent emission. In particular, a phase shift is clearly observed in many cases. We interpret these findings as direct evidence of changes of the pulsed beam pattern, due to transitions between the sub- and super-critical accretion regimes on a short time scale. We also observe that at comparable luminosities the flares' pulse profiles are rather similar to those of the persistent emission. This indicates that the accretion on the polar cap of the neutron star is mainly determined by the luminosity, i.e., the mass accretion rate.

We perform a systematic search of quasi-periodic variabilities in $\gamma$-ray active galactic nuclei (AGNs) in the third \emph{Fermi} Large Area Telescope source catalog (3FGL). We employe two techniques, Lomb-Scargle Periodogram (LSP) and Weighted Wavelet Z-transform (WWZ), to obtain power spectra of $\gamma$-ray light curves covering from 2008 August to 2016 December. The results show that besides several objects that have been reported in previous works, an additional source, the FSRQ PKS 0601-70 has a possible quasi-periodic variability of 450 days in its $\gamma$-ray light curves with the significance of $>3\sigma$. The physical implications of our findings are discussed.

A new black hole X-ray binary (BHXRB) MAXI J1535-571 was discovered by MAXI during its outburst in 2017. Using observations taken by the first Chinese X-ray satellite, the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT), we perform a joint spectral analysis (2-150 keV) in both energy and time domains. The energy spectra provide the essential input for probing the intrinsic Quasi-Periodic Oscillation (QPO) fractional rms spectra (FRS). Our results show that during the intermediate state, the energy spectra are in general consistent with those reported by Swift/XRT and NuSTAR. However, the QPO FRS become harder and the FRS residuals may suggest the presence of either an additional power-law component in the energy spectrum or a turn-over in the intrinsic QPO FRS at high energies.

We study the 2018 outburst of Aql X-1 via the monitor of all sky X-ray image (MAXI) data. We show that the outburst starting in February 2018 is a member of short-low class in the frame of outburst duration and the peak count rate although the outburst morphology is slightly different from the other fast-rise-exponential-decay (FRED) type outbursts with a milder rising stage. We study the partial accretion in the weak propeller stage of Aql X-1 via the MAXI data of the 2018 outburst. We report on the spectral analysis of 3 observations of Aquila X-1 obtained by Insight - hard X-ray modulation telescope (Insight-HXMT) during the late decay stage of the 2018 outburst. We discuss that the data taken by Insight-HXMT is just after the transition to the weak propeller stage. Our analysis shows the necessity of a comptonization component to take into account the existence of an electron cloud resulting photons partly up-scattered.

We search through $\gamma$-ray data obtained with the Large Area Telescope (LAT) onboard the Fermi Gamma-Ray Space Telescope and find 24 blazars (or candidates) that have a single clear flare event in their 9.5 year long-term light curves. We define these events as loner flares since each flare stands out significantly above the relatively stable, low-flux light curve. We analyze the LAT data in detail for these 24 sources. The flares in ten of them are primarily due to a single sharp peak, for which we study by fitting with two different analytic functions. The time durations thus determined for the sharp peaks are in a range of 4-25 days. The $\gamma$-ray spectra of the 24 blazar sources can be described with a power-law or a log-parabola function. We obtain their spectral properties in the flaring and quiescent states, and find that in the flares 16 of the sources have harder emission and three have softer emission while the other five keep the same emission. We discuss a possible correlation between the differences in photon index in the quiescent and flaring states and photon indices in quiescence. In addition, the sharp peak flares seem to have a tendency of having long time durations and hard emission, possibly related to their physical origin in a blazar jet. Studies of more similar flares will help establish these possible features.

We present a detailed timing study of the brightest persistent X-ray source Sco X-1 using the data collected by the Hard X-ray Modulation Telescope ($Insight$-HXMT) from July 2017 to August 2018. A complete $Z$-track hardness-intensity diagram (HID) is obtained. The normal branch oscillations (NBOs) at $\sim$ 6 Hz in the lower part of the normal branch (NB) and the flare branch oscillations (FBOs) at $\sim$ 16 Hz in the beginning part of the flaring branch (FB) are found in observations with the Low Energy X-ray Telescope (LE) and the Medium Energy X-ray Telescope (ME) of $Insight$-HXMT, while the horizontal branch oscillations (HBOs) at $\sim$ 40 Hz and the kilohertz quasi-periodic oscillations (kHz QPOs) at $\sim$ 800 Hz are found simultaneously up to 60 keV for the first time on the horizontal branch (HB) by the High Energy X-ray Telescope (HE) and ME. We find that for all types of the observed QPOs, the centroid frequencies are independent of energy, while the root mean square (rms) increases with energy; the centroid frequencies of both the HBOs and kHz QPOs increase along the $Z$-track from the top to the bottom of the HB; and the NBOs show soft phase lags increasing with energy. A continuous QPO transition from the FB to NB in $\sim$ 200 s are also detected. Our results indicate that the non-thermal emission is the origin of all types of QPOs, the innermost region of the accretion disk is non-thermal in nature, and the corona is nonhomogeneous geometrically.

The evidences for the influence of thermonuclear (type-I) X-ray bursts upon the surrounding environments in neutron star low-mass X-ray binaries (LMXB) were detected previously via spectral and timing analyses. Benefitting from a broad energy coverage of Insight-HXMT, we analyze one photospheric radius expansion (PRE) burst, and find an emission excess at soft X-rays. Our spectral analysis shows that, such an excess is not likely relevant to the disk reflection induced by the burst emission and can be attributed to an enhanced pre-burst/persistent emission. We find that the burst and enhanced persistent emissions sum up to exceed Eddington luminosity by $\sim$ 40 percentages. We speculate that the enhanced emission is from a region beyond the PRE radius, or through the Comptonization of the corona.

The Hard X-ray Modulation Telescope (HXMT) named Insight is China's first X-ray astronomical satellite. The Low Energy X-ray Telescope (LE) is one of its main payloads onboard. The detectors of LE adopt swept charge device CCD236 with L-shaped transfer electrodes. Charges in detection area are read out continuously along specific paths, which leads to a time response distribution of photons readout time. We designed a long exposure readout mode to measure the time response distribution. In this mode, CCD236 firstly performs exposure without readout, then all charges generated in preceding exposure phase are read out completely. Through analysis of the photons readout time in this mode, we obtained the probability distribution of photons readout time.

The long-term evolution of the centroid energy of the CRSF in Her X-1 is still a mystery. We report a new measurement from a campaign between {\sl Insight}-HXMT and {\sl NuSTAR} performed in February 2018. Generally, the two satellites show well consistent results of timing and spectral properties. The joint spectral analysis confirms that the previously observed long decay phase has ended, and that the line energy instead keeps constant around 37.5 keV after flux correction.

Considering that the existence of relativistic particles in the protostellar jet has been confirmed by the detection of linearly polarized radio emission from the HH 80-81 jet, we search for $\gamma$-rays (which could be produced by the relativistic particles) from the HH 80-81 system using ten-year {\it Fermi}-LAT observations. A significant $\gamma$-ray excess is found in the direction of the HH 80-81 system with Test-Statistic (TS) value$>$100 , which is likely produced in the HH 80-81 jet. The $\gamma$-ray spectrum extends only to 1 GeV with a photon index of 3.5. No significant variability is found in the $\gamma$-ray emission. Combining the $\gamma$-ray and previous radio observations, we may exclude the relativistic Bremsstrahlung process for the production of the $\gamma$-rays under the assumption of a magnetic field strength of 0.1 mG. $\gamma$-rays can be produced through inverse-Compton (IC) scattering of relativistic electrons or inelastic proton-proton ($pp$) interaction. The maximum energies of electrons and protons are constrained to be $<0.1\ $TeV and $\sim\ $20 GeV, respectively.

Cosmic microwave background (CMB) measurements are fundamentally limited by photon statistics. Therefore, ground-based CMB observatories have been increasing the number of detectors that are simultaneously observing the sky. Thanks to the advent of monolithically fabricated transition edge sensor (TES) arrays, the number of on-sky detectors has been increasing exponentially for over a decade. The next-generation experiment CMB-S4 will increase this detector count by more than an order of magnitude from the current state-of-the-art to ~500,000. The readout of such a huge number of exquisitely precise sub-Kelvin sensors is feasible using an existing technology: frequency-domain multiplexing (fMux). To further optimize this system and reduce complexity and cost, we have recently made significant advances including the elimination of 4 K electronics, a massive decrease of parasitic in-series impedances, and a significant increase in multiplexing factor.

We selected 90 massive star-forming clumps with strong N2H+, HCO+, HCN, and HNC emission from the Millimetre Astronomy Legacy Team 90 GHz survey. We obtained Herschel data for all 90 sources and NRAO VLA Sky Survey data for 51 of them. We convolved and regridded all images to the same resolution and pixel size and derived the temperature, H2 column density, molecules' abundances and abundance, and ratios of each pixel. Our analysis yields three main conclusions. First, the abundances of N2H+, HCO+, HCN, and HNC increase when the column density decreases and the temperature increases, with spatial variations in their abundances dominated by changes in the H2 column density. Second, the abundance ratios between N2H+, HCO+, HCN, and HNC also display systemic variations as a function of the column density due to the chemical properties of these molecules. Third, the sources associated with the 20 cm continuum emission can be classified into four types based on the behavior of the abundances of the four molecules considered here as a function of this emission. The variations of the first three types could also be attributed to the variation of the H2 column density.

A total of 188 high-mass outflows have been identified from a sample of 694 clumps from the Millimetre Astronomy Legacy Team 90 GHz survey, representing a detection rate of approximately 27%. The detection rate of outflows increases from the protostellar stage to the H II stage, but decreases again at the photodissociation (PDR) stage suggesting that outflows are being switched off during the PDR stage. An intimate relationship is found between outflow action and the presence of masers, and water masers appear together with 6.7 GHz methanol masers. Comparing the infall detection rate of clumps with and without outflows, we find that outflow candidates have a lower infall detection rate. Finally, we find that outflow action has some influence on the local environment and the clump itself, and this influence decreases with increasing evolutionary time as the outflow action ceases.

We derive molecular-gas-phase $^{12}$C/$^{13}$C isotope ratios for the central few 100 pc of the three nearby starburst galaxies NGC 253, NGC 1068, and NGC 4945 making use of the $\lambda$ $\sim$ 3 mm $^{12}$CN and $^{13}$CN $N$ = 1--0 lines in the ALMA Band 3. The $^{12}$C/$^{13}$C isotopic ratios derived from the ratios of these lines range from 30 to 67 with an average of 41.6 $\pm$ 0.2 in NGC 253, from 24 to 62 with an average of 38.3 $\pm$ 0.4 in NGC 1068, and from 6 to 44 with an average of 16.9 $\pm$ 0.3 in NGC 4945. The highest $^{12}$C/$^{13}$C isotopic ratios are determined in some of the outskirts of the nuclear regions of the three starburst galaxies. The lowest ratios are associated with the northeastern and southwestern molecular peaks of NGC 253, the northeastern and southwestern edge of the mapped region in NGC 1068, and the very center of NGC 4945. In case of NGC 1068, the measured ratios suggest inflow from the outer part of NGC 1068 into the circum-nuclear disk through both the halo and the bar. Low $^{12}$C/$^{13}$C isotopic ratios in the central regions of these starburst galaxies indicate the presence of highly processed material.

A multi-wavelength analysis of the large Galactic infrared bubble N 24 is has been presented in this paper in order to investigate the molecular and star formation environment around expanding HII regions. Using archival data from Herschel and ATLASGAL, the distribution and physical properties of the dust over the entire bubble are studied. Twenty three dense clumps are identified using the Clumpfind2d algorithm with sizes and masses in the range 0.65-1.73 pc and 600-16300 Msun, respectively. To analyse the molecular environment in N 24, observations of NH3 (1,1) and (2,2) were carried out using the Nanshan 26m radio telescope. Analysis of the kinetic temperature and gravitational stability of these clumps suggests gravitational collapse in several of them. The mass-size distributions of the clumps and the presence of massive young protostars indicate that the shell of N 24 is a region of ongoing massive star formation. The compatibility of the dynamical and fragmentation time-scales and the overabundance of YSOs and clumps on the rim suggest that the collect and collapse mechanism is in play at the boundary of the bubble, but the existence of the IRDC at the edge of bubble indicates that radiation-driven implosion mechanism may also have played a role there.

We present the POL-2 850 $\mu$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We then use the polarization data at 850 $\mu$m to infer the plane-of-sky orientation of the large-scale magnetic field in Barnard 1. This magnetic field runs North-South across most of the cloud, with the exception of B1-c where it turns more East-West. From the dispersion of polarization angles, we calculate a turbulence correlation length of $5.0 \pm 2.5$ arcsec ($1500$ au), and a turbulent-to-total magnetic energy ratio of $0.5 \pm 0.3$ inside the cloud. We combine this turbulent-to-total magnetic energy ratio with observations of NH$_3$ molecular lines from the Green Bank Ammonia Survey (GAS) to estimate the strength of the plane-of-sky component of the magnetic field through the Davis-Chandrasekhar-Fermi method. With a plane-of-sky amplitude of $120 \pm 60$ $\mu$G and a criticality criterion $\lambda_c = 3.0 \pm 1.5$, we find that Barnard 1 is a supercritical molecular cloud with a magnetic field nearly dominated by its turbulent component.

The formaldehyde $\rm H_2CO(1_{10} - 1_{11})$ absorption line and H$110\alpha$ radio recombination line (RRL) have been observed toward the Aquila Molecular Cloud using the Nanshan 25 m telescope operated by the Xinjiang Astronomical Observatory CAS. These first observations of the $\rm H_2CO$ $(1_{10} - 1_{11})$ absorption line determine the extent of the molecular regions that are affected by the ongoing star formation in the Aquila molecular complex and show some of the dynamic properties. The distribution of the excitation temperature $T_{ex}$ for $\rm H_2CO$ identifies the two known star formation regions W40 and Serpens South as well as a smaller new region Serpens 3. The intensity and velocity distributions of $\rm H_2CO$ and $\rm ^{13}CO(1-0)$ do not agree well with each other, which confirms that the $\rm H_2CO$ absorption structure is mostly determined by the excitation of the molecules resulting from the star formation rather than by the availability of molecular material as represented by the distribution. Some velocity-coherent linear $\rm ^{13}CO(1-0)$ structures have been identified in velocity channel maps of $\rm H_2CO$ and it is found that the three star formation regions lie on the intersect points of filaments. The $\rm H110\alpha$ emission is found only at the location of the W40 H II region and spectral profile indicates a redshifted spherical outflow structure in the outskirts of the H II region. Sensitive mapping of $\rm H_2CO$ absorption of the Aquila Complex has correctly identified the locations of star-formation activity in complex molecular clouds and the spectral profiles reveal the dominant velocity components and may identify the presence of outflows.

We present the 850 $\mu$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scale ($\lesssim 1.0$ pc) hub-filament structure (HFS) located at the end of a parsec-scale filament. To investigate whether or not the observed polarization traces the magnetic field in the HFS, we analyze the dependence between the observed polarization fraction and total intensity using a Bayesian approach with the polarization fraction described by the Rice likelihood function, which can correctly describe the probability density function (PDF) of the observed polarization fraction for low signal-to-noise ratio (SNR) data. We find a power-law dependence between the polarization fraction and total intensity with an index of 0.56 in $A_V\sim$ 20--300 mag regions, suggesting that the dust grains in these dense regions can still be aligned with magnetic fields in the IC5146 regions. Our polarization maps reveal a curved magnetic field, possibly dragged by the contraction along the parsec-scale filament. We further obtain a magnetic field strength of 0.5$\pm$0.2 mG toward the central hub using the Davis-Chandrasekhar-Fermi method, corresponding to a mass-to-flux criticality of $\sim$ $1.3\pm0.4$ and an Alfv\'{e}nic Mach number of $<$0.6. These results suggest that gravity and magnetic field is currently of comparable importance in the HFS, and turbulence is less important.

This paper reports on the discovery that an eclipsing binary system, EPIC 202843107 , has a {\delta} Scuti variable component. The phased light curve from Kepler space telescope presents a detached configuration. The binary modelling indicates that the two component stars have almost the same radius and may have experienced orbital circularization. Frequency analyses are performed for the residual light curve after subtracting the binary variations. The frequency spectrum reveals that one component star is a {\delta} Scuti variable. A large frequency separation is cross-identified with the histogram graph, the Fourier transform, and the echelle diagram method. The mean density of the {\delta} Scuti component is estimated to be 0.09 g/cm3 based on the large separation and density relation. Systems like EPIC 202843107 are helpful to study the stellar evolution and physical state for binary stars.

We report 850~$\mu$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $\rho$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic field across the $\sim$0.1~pc core shows a predominant northeast-southwest orientation centering between $\sim$40$^\circ$ to $\sim$100$^\circ$, indicating that the field in the core is well aligned with the magnetic field in lower-density regions of the cloud probed by near-infrared observations and also the cloud-scale magnetic field traced by Planck observations. The polarization percentage ($P$) decreases with an increasing total intensity ($I$) with a power-law index of $-$1.03 $\pm$ 0.05. We estimate the plane-of-sky field strength ($B_{\mathrm{pos}}$) using modified Davis-Chandrasekhar-Fermi (DCF) methods based on structure function (SF), auto-correlation (ACF), and unsharp masking (UM) analyses. We find that the estimates from the SF, ACF, and UM methods yield strengths of 103 $\pm$ 46 $\mu$G, 136 $\pm$ 69 $\mu$G, and 213 $\pm$ 115 $\mu$G, respectively. Our calculations suggest that the Ophiuchus C core is near magnetically critical or slightly magnetically supercritical (i.e. unstable to collapse). The total magnetic energy calculated from the SF method is comparable to the turbulent energy in Ophiuchus C, while the ACF method and the UM method only set upper limits for the total magnetic energy because of large uncertainties.

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.

We report first results from the CHinese Exoplanet Searching Program from Antarctica (CHESPA)---a wide-field high-resolution photometric survey for transiting exoplanets carried out using telescopes of the AST3 (Antarctic Survey Telescopes times 3) project. There are now three telescopes (AST3-I, AST3-II, and CSTAR-II) operating at Dome A---the highest point on the Antarctic Plateau---in a fully automatic and remote mode to exploit the superb observing conditions of the site, and its long and uninterrupted polar nights. The search for transiting exoplanets is one of the key projects for AST3. During the Austral winters of 2016 and 2017 we used the AST3-II telescope to survey a set of target fields near the southern ecliptic pole, falling within the continuous viewing zone of the TESS mission \citep{Ricker10}. The first data release of the 2016 data, including images, catalogs and lightcurves of 26578 bright stars ($7.5\le i \le15$) was presented in \citet{Zhang18}. The best precision, as measured by the RMS of the lightcurves at the optimum magnitude of the survey ($i=10$), is around 2\,mmag. We detect 222 objects with plausible transit signals from these data, 116 of which are plausible transiting exoplanet candidates according to their stellar properties as given by the TESS Input Catalog \citep{Stassun17}, Gaia DR2 \citep{Gaia18} and TESS-HERMES spectroscopy \citep{Sharma18}. With the first data release from TESS expected in late 2018, this candidate list will be a timely for improving the rejection of potential false positives.

Located at Dome A, the highest point of the Antarctic plateau, the Chinese Kunlun station is considered to be one of the best ground-based photometric sites because of its extremely cold, dry, and stable atmosphere(Saunders et al. 2009). A target can be monitored from there for over 40 days without diurnal interruption during a polar winter. This makes Kunlun station a perfect site to search for short-period transiting exoplanets. Since 2008, an observatory has been built at Kunlun station and three telescopes are working there. Using these telescopes, the AST3 project has been carried out over the last six years with a search for transiting exoplanets as one of its key programs (CHESPA). In the austral winters of 2016 and 2017, a set of target fields in the Southern CVZ of TESS (Ricker et al. 2009) were monitored by the AST3-II telescope. In this paper, we introduce the CHESPA and present the first data release containing photometry of 26,578 bright stars (m_i < 15). The best photometric precision at the optimum magnitude for the survey is around 2 mmag. To demonstrate the data quality, we also present a catalog of 221 variables with a brightness variation greater than 5 mmag from the 2016 data. Among these variables, 179 are newly identified periodic variables not listed in the AAVSO databasea), and 67 are listed in the Candidate Target List(Stassun et al. 2017). These variables will require careful attention to avoid false-positive signals when searching for transiting exoplanets. Dozens of new transiting exoplanet candidates will be also released in a subsequent paper(Zhang et al. 2018b).

Planetary systems can survive the stellar evolution, as evidenced by the atmospheric metal pollution and dusty disks of single white dwarfs. Recent observations show that 1 to 4 percent of single white dwarfs are accompanied by dusty disks, while the occurrence rate of metal pollution is about 25 to 50 percent . The dusty disks and metal pollution have been associated with accretion of remnant planetary systems around white dwarfs, yet the relation between these two phenomena is still unclear. Here we suggest an evolutionary scenario to link the two observational phenomena. By analyzing a sample of metal polluted white dwarfs, we find that the mass accretion rate onto the white dwarf generally follows a broken power law decay, which matches well with the theoretical prediction, if assuming dust accretion is primarily driven by Poynting-Robertson drag and the dust source is primarily delivered via dynamically falling asteroids perturbed by a Jovian planet. The presence of disks is mainly at the early stage (about 0.1 to 0.7 billion years) of the whole process of metal pollution, which is detectable until about 8 billion years, naturally explaining the fraction (about 2-16 percent) of metal-polluted white dwarfs having dusty disks. The success of this scenario also implies that the configuration of an asteroid belt with an outer gas giant might be common around stars of several solar masses.

A relativistic electron-positron pair beam can be produced in the interaction of TeV photons from a blazar with the extragalactic background light (EBL). The relativistic $e^{\pm}$ pairs would loss energy through inverse-Compton scattering (ICS) photons of cosmic microwave background (CMB) or plasma instabilities. The dominant energy-loss process is under debate. Based on the assumption that the dominant energy-loss process is ICS, the resulted cascade GeV radiation is usually used to constrain the intergalactic magnetic field (IGMF). Here, we include the energy-loss due to plasma oblique instability in the calculation of cascade gamma-ray flux, and investigate the impact of the plasma instability on the constraint of IGMF. The up-to-date GeV data and archival TeV data of the blazar 1ES 0229+200 are used. The results indicate that even if the oblique instability cooling is dominating over ICS cooling, the cascade flux could be still used to constrain the IGMF. It is found that with the ratio between the cooling rates of the oblique instability and the ICS varying from 0.1, 1 to 10, the lower limit of the IGMF putted by the cascade flux and the gamma-ray data changes from $8\times10^{-18}\ $G, $5\times10^{-18}\ $G to $10^{-18}\ $G. If the ratio between the two cooling rates is 30, the estimate of IGMF based on the cascade flux is invalid.

Since 2016 October, the active galaxy PKS 2247-131 has undergone a gamma-ray outburst, which we studied using data obtained with the Fermi Gamma-ray Space Telescope. The emission arises from a relativistic jet in PKS 2247-131, as an optical spectrum only shows a few weak absorption lines, typical of the BL Lacertae sub-class of the blazar class of active galactic nuclei. Here we report a ~34.5 day quasi-periodic oscillation (QPO) in the emission after the initial flux peak of the outburst. Compared to one-year time-scale QPOs, previously identified in blazars in Fermi energies, PKS 2247-131 exhibits the first clear case of a relatively short, month-like oscillation. We show that this QPO can be explained in terms of a helical structure in the jet, where the viewing angle to the dominant emission region in the jet undergoes periodic changes. The time scale of the QPO suggests the presence of binary supermassive black holes in PKS 2247-131.

The Globular Cluster 47 Tucanae was firstly detected in gamma-rays by the Large Area Telescope (LAT) onboard the \emph{Fermi} Gamma-ray Space Telescope, and the gamma-ray emission has been widely attributed to the millisecond pulsars. In this work, we analyze the Fermi-LAT pass 8 data ranging from 2008 August to 2017 May and report the detection of a modulation with a period of $18.416\pm0.008$ hours at a significance level of $\sim 4.8\sigma$. This is the first time to detect a significant modulation with a period much longer than that of millisecond pulsars in gamma-rays from Globular Clusters. The periodic modulation signal displays in the {\it Swift}-BAT data as well. The phase-folded Chandra X-ray light curve of a point source may have provided an additional clue. These facts indicate the detection of the first young pulsar binary in gamma-rays in Globular Clusters.

The mechanism for producing gamma-ray quasi-periodic oscillation (QPO) in blazar is unknown. One possibility is the geometric model, in which without the need for intrinsic quasi-periodic variation, the relativistic Doppler factor changes periodically, resulting in observed gamma-ray QPO. We propose a method to test this geometric model. We analyze the {\it Fermi}-LAT data of PG 1553+113 spanning from 2008 August until 2018 February. According to 29 four-month average spectral energy distributions (SEDs) in the energy range of 0.1-300 GeV, we split the {\it Fermi}-LAT energy range into three bands: 0.1-1 GeV, 1-10 GeV, and 10-300 GeV. The spectrum in each energy range can be successfully fitted by a power-law. The light curves and photon indices in the three energy ranges are obtained. Then, light curves in three narrow energy ranges, i.e., 0.2-0.5 GeV, 2-5 GeV and 20- 40 GeV, are constructed, and the relative variability amplitudes in the three narrow energy ranges are calculated. A discrete-correlation analysis is performed for the light curves. Our results indicate that (i) the light curves in the different energy ranges follow the same pattern showed in the light curve above 0.1 GeV; (ii) the three groups of photon indices in the energy ranges of 0.1-1 GeV, 1-10 GeV, and 10-300 GeV keep nearly constant; (iii) the ratio between relative variability amplitudes in different narrow energy ranges are equal (within their errors) to the prediction by the Doppler effect. Our results support the scenario of the relativistic boost producing the gamma-ray QPO for PG 1553+113.

An unbiased search of molecular outflows within the region of the COHRS survey has identified 157 high-mass outflows from a sample of 770 ATLASGAL clumps with a detection rate of 20\%. The detection rate of outflows increases for clumps with higher M$_{clump}$, L$_{bol}$, L$_{bol}$/M$_{clump}$, N$_{H_{2}}$, and T$_{dust}$ compared to the clumps with no outflow. The detection rates of the outflow increases from protostellar (8\%) to YSO clump (17\%) and to MSF clump (29\%). The detection rate 26\% for quiescent clump is preliminary, because the sample of quiescent clumps is small. A statistical relation between the outflow and clump masses for our sample is $\log(M_{out}/M_{\bigodot}) = (-1.1\pm0.21) + (0.9\pm0.07)\log(M_{clump}/M_{\bigodot})$. The detection rate of outflows and the outflow mass-loss rate show an increase with increasing M$_{clump}$, L$_{bol}$, N$_{H_{2}}$, and T$_{dust}$, which indicates that clumps with outflow with higher parameter values are at a more advanced evolutionary stage. The outflow mechanical force increases with increasing bolometric luminosities. No clear evidence has yet been found that higher mass outflows have different launching conditions than low-mass outflows.

We present the X-ray timing results of the new black hole candidate (BHC) MAXI J1535-571 during its 2017 outburst from Hard X-ray Modulation Telescope (\emph{Insight}-HXMT) observations taken from 2017 September 6 to 23. Following the definitions given by \citet{Belloni2010}, we find that the source exhibits state transitions from Low/Hard state (LHS) to Hard Intermediate state (HIMS) and eventually to Soft Intermediate state (SIMS). Quasi-periodic oscillations (QPOs) are found in the intermediate states, which suggest different types of QPOs. With the large effective area of \emph{Insight}-HXMT at high energies, we are able to present the energy dependence of the QPO amplitude and centroid frequency up to 100 keV which is rarely explored by previous satellites. We also find that the phase lag at the type-C QPOs centroid frequency is negative (soft lags) and strongly correlated with the centroid frequency. By assuming a geometrical origin of type-C QPOs, the source is consistent with being a high inclination system.

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.

Although the formation and the properties of hot Jupiters (with orbital periods $P<10\,$d) have attracted a great deal of attention, the origins of warm Jupiters ($10<P<100\,$d) are less well-studied. Using a transit timing analysis, we present the orbital parameters of five planetary systems containing warm Jupiters, Kepler 30, Kepler 117, Kepler 302, Kepler 487 and Kepler 418. Three of them, Kepler-30 c($M_p=549.4{\pm{5.6}}M_{\oplus}$), Kepler-117 c($M_p=702{\pm{63}}M_{\oplus}$) and Kepler 302 c($M_p=933\pm 527M_{\oplus}$), are confirmed to be real warm Jupiters based on their mass. Insights drawn from the radius-temperature relationship lead to the inference that hot Jupiters and warm Jupiters can be roughly separated by ${T_{\rm eff,c}} =1123.7\pm3.3$ K. Also, ${T_{\rm eff,c}}$ provides a good separation for Jupiters with companion fraction consistent with zero($T_{\rm eff}>T_{\rm eff,c}$) and those with companion fraction significantly different from zero ($T_{\rm eff}<T_{\rm eff,c}$).

We present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope (JCMT) in Hawaii. This work is part of the B-fields In Star-forming Region Observations (BISTRO) survey initiated to understand the role of magnetic fields in star formation for nearby star-forming molecular clouds. We present a first look at the geometry and strength of magnetic fields in Oph-B. The field geometry is traced over $\sim$0.2 pc, with clear detection of both of the sub-clumps of Oph-B. The field pattern appears significantly disordered in sub-clump Oph-B1. The field geometry in Oph-B2 is more ordered, with a tendency to be along the major axis of the clump, parallel to the filamentary structure within which it lies. The degree of polarization decreases systematically towards the dense core material in the two sub-clumps. The field lines in the lower density material along the periphery are smoothly joined to the large scale magnetic fields probed by NIR polarization observations. We estimated a magnetic field strength of 630$\pm$410 $\mu$G in the Oph-B2 sub-clump using a Davis-Chandeasekhar-Fermi analysis. With this magnetic field strength, we find a mass-to-flux ratio $\lambda$= 1.6$\pm$1.1, which suggests that the Oph-B2 clump is slightly magnetically supercritical.

Stars are commonly formed in binary systems, which provide a natural laboratory for studying planet formation in extreme conditions. In our first paper (Paper I) of a series Xie et al. (2011), we have shown that the intermediate stage - from planetesimals to planetary embryos/cores - of planet formation can proceed even in highly inclined binaries. Following Paper I, here we numerically study the late stage of terrestrial planet formation, i.e., from embryos to full planets, in binary systems of various orbital configurations. We identify an orbital alignment or anti-alignment effect; namely, although an inclined binary generally misaligns the planetary orbits with respect to the spin axis of the primary host star (i.e., causing large obliquity), it could align or anti-align the planetary orbits with respect to the binary orbit. Such an orbital (anti-)alignment effect is caused by the combination of orbital differential precession and self-damping, and it is mostly significant in cases of intermediate binary separations, i.e., a_B~40-200 AU for terrestrial planet formation around 1 AU from the primary stars. In such intermediate separation binaries, somewhat contrary to intuition, the binary companion can aid planet growth by having increased the rate of collisions, forming significantly more massive but fewer planets. In the other two ends, the companion is ether too close thus plays a violently disruptive role or too wide to have significant effect on planet formation. Future observations, which can discover more planet-bearing binary star systems and constrain their masses and 3-D orbital motions will test our numerical findings.

The three Antarctic Survey Telescopes (AST3) aim to carry out time domain imaging survey at Dome A, Antarctica. The first of the three telescopes (AST3-1) was successfully deployed on January 2012. AST3-1 is a 500\,mm aperture modified Schmidt telescope with a 680\,mm diameter primary mirror. AST3-1 is equipped with a SDSS $i$ filter and a 10k $\times$ 10k frame transfer CCD camera, reduced to 5k $\times$ 10k by electronic shuttering, resulting in a 4.3 deg$^2$ field-of-view. To verify the capability of AST3-1 for a variety of science goals, extensive commissioning was carried out between March and May 2012. The commissioning included a survey covering 2000 deg$^2$ as well as the entire Large and Small Magellanic Clouds. Frequent repeated images were made of the center of the Large Magellanic Cloud, a selected exoplanet transit field, and fields including some Wolf-Rayet stars. Here we present the data reduction and photometric measurements of the point sources observed by AST3-1. We have achieved a survey depth of 19.3\,mag in 60 s exposures with 5\,mmag precision in the light curves of bright stars. The facility achieves sub-mmag photometric precision under stable survey conditions, approaching its photon noise limit. These results demonstrate that AST3-1 at Dome A is extraordinarily competitive in time-domain astronomy, including both quick searches for faint transients and the detection of tiny transit signals.

In order to understand the initial conditions and early evolution of star formation in a wide range of Galactic environments, we carried out an investigation of 64 \textit{Planck} Galactic Cold Clumps (PGCCs) in the second quadrant of the Milky Way. Using the $^{13}$CO and C$^{18}$O $J = 1 - 0$ lines, and 850\,$\mu$m continuum observations, we investigated cloud fragmentation and evolution associated with star formation. We extracted 468 clumps and 117 cores from the $^{13}$CO line and 850\,$\mu$m continuum maps, respectively. We make use of the Bayesian Distance Calculator and derived the distances of all 64 PGCCs. We found that in general, the mass-size plane follows a relation of $m\sim r^{1.67}$. At a given scale, the masses of our objects are around 1/10 of that of typical Galactic massive star-forming regions. Analysis of the clump and core masses, virial parameters, densities, and mass-size relation suggests that the PGCCs in our sample have a low core formation efficiency ($\sim$3.0\%), and most PGCCs are likely low-mass star-forming candidates. Statistical study indicates that the 850\,$\mu$m cores are more turbulent, more optically thick, and denser than the $^{13}$CO clumps for star formation candidates, suggesting that the 850\,$\mu$m cores are likely more appropriate future star-formation candidates than the $^{13}$CO clumps.

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 present 850 $\mu$m imaging polarimetry data of the $\rho$ Oph-A core taken with the Submillimeter Common-User Bolometer Array-2 (SCUBA-2) and its polarimeter (POL-2), as part of our ongoing survey project, BISTRO (B-fields In STar forming RegiOns). The polarization vectors are used to identify the orientation of the magnetic field projected on the plane of the sky at a resolution of 0.01 pc. We identify 10 subregions with distinct polarization fractions and angles in the 0.2 pc $\rho$ Oph A core; some of them can be part of a coherent magnetic field structure in the $\rho$ Oph region. The results are consistent with previous observations of the brightest regions of $\rho$ Oph-A, where the degrees of polarization are at a level of a few percents, but our data reveal for the first time the magnetic field structures in the fainter regions surrounding the core where the degree of polarization is much higher ($> 5 \%$). A comparison with previous near-infrared polarimetric data shows that there are several magnetic field components which are consistent at near-infrared and submillimeter wavelengths. Using the Davis-Chandrasekhar-Fermi method, we also derive magnetic field strengths in several sub-core regions, which range from approximately 0.2 to 5 mG. We also find a correlation between the magnetic field orientations projected on the sky with the core centroid velocity components.

Exoplanets discovered over the last decades have provided a new sample of giant exoplanets, hot Jupiters. For lack of enough materials in current locations of hot Jupiters, they are perceived to form outside snowline. Then, migrate to the locations observed through interactions with gas disks or high-eccentricity mechanisms. We examined the efficiencies of different high-eccentricity mechanisms to form hot Jupiters in near coplaner multi-planet systems. These mechanisms include planet-planet scattering, Kozai-Lidov mechanism, coplanar high-eccentricity migration, secular chaos, as well as other two new mechanisms we find in this work, which can produce hot Jupiters with high inclinations even retrograde. We find Kozai-Lidov mechanism plays the most important role in producing hot Jupiters among these mechanisms. Secular chaos is not the usual channel for the formation of hot Jupiters due to the lack of angular momentum deficit within 10^7 Tin (periods of the inner orbit). According to comparisons between the observations and simulations, we speculate that there are at least two populations of hot Jupiters. One population migrates into the boundary of tidal effects due to interactions with gas disk, such as ups And b, WASP-47 b and HIP 14810 b. These systems usually have at least two planets with lower eccentricities, and keep dynamical stable in compact orbital configurations. Another population forms through high-eccentricity mechanisms after the excitation of eccentricity due to dynamical instability. This kind of hot Jupiters usually has Jupiterlike companions in distant orbits with moderate or high eccentricities.

We discover a population of short-period, Neptune-size planets sharing key similarities with hot Jupiters: both populations are preferentially hosted by metal-rich stars, and both are preferentially found in Kepler systems with single transiting planets. We use accurate LAMOST DR4 stellar parameters for main-sequence stars to study the distributions of short-period 1d < P < 10d Kepler planets as a function of host star metallicity. The radius distribution of planets around metal-rich stars is more "puffed up" as compared to that around metal-poor hosts. In two period-radius regimes, planets preferentially reside around metal-rich stars, while there are hardly any planets around metal-poor stars. One is the well-known hot Jupiters, and the other is a population of Neptune-size planets (2 R_Earth <~ R_p <~ 6 R_Earth), dubbed as "Hoptunes". Also like hot Jupiters, Hoptunes occur more frequently in systems with single transiting planets though the fraction of Hoptunes occurring in multiples is larger than that of hot Jupiters. About 1% of solar-type stars host "Hoptunes", and the frequencies of Hoptunes and hot Jupiters increase with consistent trends as a function of [Fe/H]. In the planet radius distribution, hot Jupiters and Hoptunes are separated by a "valley" at approximately Saturn size (in the range of 6 R_Earth <~ R_p <~ 10 R_Earth), and this "hot-Saturn valley" represents approximately an order-of-magnitude decrease in planet frequency compared to hot Jupiters and Hoptunes. The empirical "kinship" between Hoptunes and hot Jupiters suggests likely common processes (migration and/or formation) responsible for their existence.

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.

Proxima Centauri is known as the closest star from the Sun. Recently, radial velocity observations revealed the existence of an Earth-mass planet around it. With an orbital period of ~11 days, the surface of Proxima Centauri b is temperate and might be habitable. We took a photometric monitoring campaign to search for its transit, using the Bright Star Survey Telescope at the Zhongshan Station in Antarctica. A transit-like signal appearing on 2016 September 8th, is identified tentatively. Its midtime, $T_{C}=2,457,640.1990\pm0.0017$ HJD, is consistent with the predicted ephemeris based on RV orbit in a 1$\sigma$ confidence interval. Time-correlated noise is pronounced in the light curve of Proxima Centauri, affecting detection of transits. We develop a technique, in a Gaussian process framework, to gauge the statistical significance of potential transit detection. The tentative transit signal reported here, has a confidence level of $2.5\sigma$. Further detection of its periodic signals is necessary to confirm the planetary transit of Proxima Centauri b. We plan to monitor Proxima Centauri in next Polar night at Dome A in Antarctica, taking the advantage of continuous darkness. \citet{Kipping17} reported two tentative transit-like signals of Proxima Centauri b, observed by the Microvariability and Oscillation of Stars space Telescope in 2014 and 2015, respectively. The midtransit time of our detection is 138 minutes later than that predicted by their transit ephemeris. If all the signals are real transits, the misalignment of the epochs plausibly suggests transit timing variations of Proxima Centauri b induced by an outer planet in this system.

We aim to directly determine the kinetic temperature and spatial density with formaldehyde for the $\sim$100 brightest ATLASGAL-selected clumps at 870 $\mu$m representing various evolutionary stages of high-mass star formation. Ten transitions ($J$ = 3-2 and 4-3) of ortho- and para-H$_2$CO near 211, 218, 225, and 291 GHz were observed with the APEX 12 m telescope. Using non-LTE models with RADEX, we derive the gas kinetic temperature and spatial density using the measured p-H$_2$CO 3$_{21}$-2$_{20}$/3$_{03}$-2$_{02}$, 4$_{22}$-3$_{21}$/4$_{04}$-3$_{03}$, and 4$_{04}$-3$_{03}$/3$_{03}$-2$_{02}$ ratios. The gas kinetic temperatures derived from the p-H$_2$CO 3$_{21}$-2$_{20}$/3$_{03}$-2$_{02}$ and 4$_{22}$-3$_{21}$/4$_{04}$-3$_{03}$ line ratios are high, ranging from 43 to $>$300 K with an unweighted average of 91 $\pm$ 4 K. Deduced $T_{\rm kin}$ values from the $J$ = 3-2 and 4-3 transitions are similar. Spatial densities of the gas derived from the p-H$_2$CO 4$_{04}$-3$_{03}$/3$_{03}$-2$_{02}$ line ratios yield 0.6-8.3 $\times$ 10$^6$ cm$^{-3}$ with an unweighted average of 1.5 ($\pm$0.1) $\times$ 10$^6$ cm$^{-3}$. A comparison of kinetic temperatures derived from p-H$_2$CO, NH$_3$, and the dust emission indicates that p-H$_2$CO traces a distinctly higher temperature than the NH$_3$ (2,2)/(1,1) transitions and the dust, tracing heated gas more directly associated with the star formation process. The H$_2$CO linewidths are found to be correlated with bolometric luminosity and increase with the evolutionary stage of the clumps, which suggests that higher luminosities tend to be associated with a more turbulent molecular medium. It seems that the spatial densities measured with H$_2$CO do not vary significantly with the evolutionary stage of the clumps. However, averaged gas kinetic temperatures derived from H$_2$CO increase with time through the evolution of the clumps.

We investigate the $\gamma$-ray and X-ray properties of the flat spectrum radio quasar (FSRQ) \4c50 at redshift $z= 1.517$. The {\it Fermi} Large Area Telescope (LAT) data indicate that this source was in an active state since 2013 July. During the active period, the $\gamma$-ray flux increased by more than a factor of three and two distinct flares were detected with the variability timescale as short as several hours. The $\gamma$-ray spectra can be well fitted by a log-parabola. From the fitting, we find a correlation between the peak energy and spectral curvature for the $\gamma$-ray spectra, which is the first time seen in $\gamma$-ray emission from a blazar. The {\it Swift} X-ray Telescope (XRT) data show that the source was variable at X-ray energies, but no evidence shows flux or spectral changes related to the $\gamma$-ray activity. The broad-band X-ray spectrum obtained with {\it Swift} XRT and {\it NuSTAR} is well described by a broken PL model, with an extremely hard spectrum ($\Gamma_{1} \sim 0.1$) below the break energy, $E_{\rm break} \sim 2.1 {\rm keV}$, and $\Gamma_{2} \sim 1.5$ above the break energy. The spectral steepening below $\sim 3$ keV is likely due to the low energy cut-off in the energy distribution of the photon-emitting electron population. Both the $\gamma$-ray and X-ray emission appear harder when brighter. The broad-band spectral energy distribution (SED) is constructed for the source, and we provide a model fit to the SED. Our modeling suggests that the emission region should be outside the broad line region, and the properties of the region indicate a jet with 42\% of the Eddington power during the active state.

The LIGO detection of gravitational waves (GW) from merging black holes in 2015 marked the beginning of a new era in observational astronomy. The detection of an electromagnetic signal from a GW source is the critical next step to explore in detail the physics involved. The Antarctic Survey Telescopes (AST3), located at Dome A, Antarctica, is uniquely situated for rapid response time-domain astronomy with its continuous night-time coverage during the austral winter. We report optical observations of the GW source (GW~170817) in the nearby galaxy NGC 4993 using AST3. The data show a rapidly fading transient at around 1 day after the GW trigger, with the $i$-band magnitude declining from $17.23\pm0.13$ magnitude to $17.72\pm0.09$ magnitude in $\sim 1.8$ hour. The brightness and time evolution of the optical transient associated with GW~170817 are broadly consistent with the predictions of models involving merging binary neutron stars. We infer from our data that the merging process ejected about $\sim 10^{-2}$ solar mass of radioactive material at a speed of up to $30\%$ the speed of light.

The discovery of the first electromagnetic counterpart to a gravitational wave signal has generated follow-up observations by over 50 facilities world-wide, ushering in the new era of multi-messenger astronomy. In this paper, we present follow-up observations of the gravitational wave event GW170817 and its electromagnetic counterpart SSS17a/DLT17ck (IAU label AT2017gfo) by 14 Australian telescopes and partner observatories as part of Australian-based and Australian-led research programs. We report early- to late-time multi-wavelength observations, including optical imaging and spectroscopy, mid-infrared imaging, radio imaging, and searches for fast radio bursts. Our optical spectra reveal that the transient source afterglow cooled from approximately 6400K to 2100K over a 7-day period and produced no significant optical emission lines. The spectral profiles, cooling rate, and photometric light curves are consistent with the expected outburst and subsequent processes of a binary neutron star merger. Star formation in the host galaxy probably ceased at least a Gyr ago, although there is evidence for a galaxy merger. Binary pulsars with short (100 Myr) decay times are therefore unlikely progenitors, but pulsars like PSR B1534+12 with its 2.7 Gyr coalescence time could produce such a merger. The displacement (about 2.2 kpc) of the binary star system from the centre of the main galaxy is not unusual for stars in the host galaxy or stars originating in the merging galaxy, and therefore any constraints on the kick velocity imparted to the progenitor are poor.

We report here the results of searching for inelastic scattering of dark matter (initial and final state dark matter particles differ by a small mass splitting) with nucleon with the first 79.6-day of PandaX-II data (Run 9). We set the upper limits for the spin independent WIMP-nucleon scattering cross section up to a mass splitting of 300 keV/c$^2$ at two benchmark dark matter masses of 1 and 10 TeV/c$^2$.

We mapped the kinetic temperature structure of the Orion molecular cloud 1 with para-H2CO(303-202, 322-221, and 321-220) using the APEX 12m telescope. This is compared with the temperatures derived from the ratio of the NH3(2,2)/(1,1) inversion lines and the dust emission. Using the RADEX non-LTE model, we derive the gas kinetic temperature modeling the measured averaged line ratios of para-H2CO 322-221/303-202 and 321-220/303-202. The gas kinetic temperatures derived from the para-H2CO line ratios are warm, ranging from 30 to >200 K with an average of 62 K at a spatial density of 10$^5$ cm$^{-3}$. These temperatures are higher than those obtained from NH3(2,2)/(1,1) and CH3CCH(6-5) in the OMC-1 region. The gas kinetic temperatures derived from para-H2CO agree with those obtained from warm dust components measured in the mid infrared (MIR), which indicates that the para-H2CO(3-2) ratios trace dense and warm gas. The cold dust components measured in the far infrared (FIR) are consistent with those measured with NH3(2,2)/(1,1) and the CH3CCH(6-5) line series. With dust at MIR wavelengths and para-H2CO(3-2) on one side and dust at FIR wavelengths, NH3(2,2)/(1,1), and CH3CCH(6-5) on the other, dust and gas temperatures appear to be equivalent in the dense gas of the OMC-1 region, but provide a bimodal distribution, one more directly related to star formation than the other. The non-thermal velocity dispersions of para-H2CO are positively correlated with the gas kinetic temperatures in regions of strong non-thermal motion (Mach number >2.5) of the OMC-1, implying that the higher temperature traced by para-H2CO is related to turbulence on a 0.06 pc scale. Combining the temperature measurements with para-H2CO and NH3(2,2)/(1,1) line ratios, we find direct evidence for the dense gas along the northern part of the OMC-1 10 km s$^{-1}$ filament heated by radiation from the central Orion nebula.

We report a new search of weakly interacting massive particles (WIMPs) using the combined low background data sets in 2016 and 2017 from the PandaX-II experiment in China. The latest data set contains a new exposure of 77.1 live day, with the background reduced to a level of 0.8$\times10^{-3}$ evt/kg/day, improved by a factor of 2.5 in comparison to the previous run in 2016. No excess events were found above the expected background. With a total exposure of 5.4$\times10^4$ kg day, the most stringent upper limit on spin-independent WIMP-nucleon cross section was set for a WIMP with mass larger than 100 GeV/c$^2$, with the lowest exclusion at 8.6$\times10^{-47}$ cm$^2$ at 40 GeV/c$^2$.

Most protostars have luminosities that are fainter than expected from steady accretion over the protostellar lifetime. The solution to this problem may lie in episodic mass accretion -- prolonged periods of very low accretion punctuated by short bursts of rapid accretion. However, the timescale and amplitude for variability at the protostellar phase is almost entirely unconstrained. In "A JCMT/SCUBA-2 Transient Survey of Protostars in Nearby Star Forming Regions", we are monitoring monthly with SCUBA-2 the sub-mm emission in eight fields within nearby (<500 pc) star forming regions to measure the accretion variability of protostars. The total survey area of ~1.6 sq.deg. includes ~105 peaks with peaks brighter than 0.5 Jy/beam (43 associated with embedded protostars or disks) and 237 peaks of 0.125-0.5 Jy/beam (50 with embedded protostars or disks). Each field has enough bright peaks for flux calibration relative to other peaks in the same field, which improves upon the nominal flux calibration uncertainties of sub-mm observations to reach a precision of ~2-3% rms, and also provides quantified confidence in any measured variability. The timescales and amplitudes of any sub-mm variation will then be converted into variations in accretion rate and subsequently used to infer the physical causes of the variability. This survey is the first dedicated survey for sub-mm variability and complements other transient surveys at optical and near-IR wavelengths, which are not sensitive to accretion variability of deeply embedded protostars.

We report a nominally high-confidence $\gamma$-ray quasi-periodic modulation in the blazar PKS 0301$-$243. For this target, we analyze its \emph{Fermi}-LAT Pass 8 data covering from 2008 August to 2017 May. Two techniques, i.e., the maximum likelihood optimization and the exposure-weighted aperture photometry, are used to build the $\gamma$-ray light curves. Then both the Lomb-Scargle Periodogram and the Weighted Wavelet Z-transform are applied to the light curves to search for period signals. A quasi-periodicity with a period of $2.1\pm0.3$ yr appears at the significance level of $\sim5\sigma$, although it should be noted that this putative quasi-period variability is seen in a data set barely four times longer. We speculate that this $\gamma$-ray quasi-periodic modulation might be evidence of a binary supermassive black hole.

We present ten $R$-band photometric observations of eight different transits of the hot Jupiter HAT-P-33b, which has been targeted by our Transiting Exoplanet Monitoring Project (TEMP). The data were obtained by two telescopes at the Xinglong Station of National Astronomical Observatories of China (NAOC) from 2013 December through 2016 January, and exhibit photometric scatter of $1.6-3.0\,\rm{mmag}$. After jointly analyzing the previously published photometric data, radial-velocity (RV) measurements, and our new light curves, we revisit the system parameters and orbital ephemeris for the HAT-P-33b system. Our results are consistent with the published values except for the planet-to-star radius ratio ($R_{P}/R_{*}$), the ingress/egress duration ($\tau$) and the total duration ($T_{14}$), which together indicate a slightly shallower and shorter transit shape. Our results are based on more complete light curves, whereas the previously published work had only one complete transit light curve. No significant anomalies in Transit Timing Variations (TTVs) are found, and we place upper mass limits on potential perturbers, largely supplanting the loose constraints provided by the extant RV data. The TTV limits are stronger near mean-motion resonances, especially for the low-order commensurabilities. We can exclude the existence of a perturber with mass larger than 0.6, 0.3, 0.5, 0.5, and $0.3\,{\rm M_\oplus}$ near the 1:3, 1:2, 2:3, 3:2, and 2:1 resonances, respectively.

We present the first results from the B-fields In STar-forming Region Observations (BISTRO) survey, using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions which the survey will aim to answer. The most important of these is the role of magnetic fields in the star formation process on the scale of individual filaments and cores in dense regions. We describe the data acquisition and reduction processes for POL-2, demonstrating both repeatability and consistency with previous data. We present a first-look analysis of the first results from the BISTRO survey in the OMC 1 region. We see that the magnetic field lies approximately perpendicular to the famous 'integral filament' in the densest regions of that filament. Furthermore, we see an 'hour-glass' magnetic field morphology extending beyond the densest region of the integral filament into the less-dense surrounding material, and discuss possible causes for this. We also discuss the more complex morphology seen along the Orion Bar region. We examine the morphology of the field along the lower-density north-eastern filament. We find consistency with previous theoretical models that predict magnetic fields lying parallel to low-density, non-self-gravitating filaments, and perpendicular to higher-density, self-gravitating filaments.

We present direct detection constraints on the absorption of hidden-photon dark matter with particle masses in the range 1.2-30 eV$c^{-2}$ with the DAMIC experiment at SNOLAB. Under the assumption that the local dark matter is entirely constituted of hidden photons, the sensitivity to the kinetic mixing parameter $\kappa$ is competitive with constraints from solar emission, reaching a minimum value of 2.2$\times$$10^{-14}$ at 17 eV$c^{-2}$. These results are the most stringent direct detection constraints on hidden-photon dark matter in the galactic halo with masses 3-12 eV$c^{-2}$ and the first demonstration of direct experimental sensitivity to ionization signals $<$12 eV from dark matter interactions.

The Large Magellanic Cloud (LMC), the closest star forming galaxy with low metallicity, provides an ideal laboratory to study star formation in such an environment. The classical dense molecular gas thermometer NH3 is rarely available in a low metallicity environment because of photoionization and a lack of nitrogen atoms. Our goal is to directly measure the gas kinetic temperature with formaldehyde toward six star-forming regions in the LMC. Three rotational transitions of para-H2CO near 218 GHz were observed with the APEX 12m telescope toward six star forming regions in the LMC. Those data are complemented by C18O 2-1 spectra. Using non-LTE modeling with RADEX, we derive the gas kinetic temperature and spatial density, using as constraints the measured para-H2CO 321-220/303-202 and para-H2CO 303-202/C18O 2-1 ratios. Excluding the quiescent cloud N159S, where only one para-H2CO line could be detected, the gas kinetic temperatures derived from the preferred para-H2CO 321-220/303-202 line ratios range from 35 to 63 K with an average of 47 K. Spatial densities of the gas derived from the paraH2CO 303-202/C18O 2-1 line ratios yield 0.4-2.9x10^5 cm^-3 with an average of 1.5x10^5 cm^-3. Temperatures derived from the para-H2CO line ratio are similar to those obtained with the same method from Galactic star forming regions and agree with results derived from CO in the dense regions of the LMC. A comparison of kinetic temperatures derived from para-H2CO with those from the dust also shows good agreement. This suggests that the dust and para-H2CO are well mixed in the studied star forming regions. A correlation between the gas kinetic temperatures derived from para-H2CO and infrared luminosity, represented by the 250um flux, suggests that the kinetic temperatures traced by para-H2CO are correlated with the ongoing massive star formation in the LMC.

AST3-1 is the second-generation wide-field optical photometric telescope dedicated to time domain astronomy at Dome A, Antarctica. Here we present the results of $i$ band images survey from AST3-1 towards one Galactic disk field. Based on time-series photometry of 92,583 stars, 560 variable stars were detected with $i$ magnitude $\leq$ 16.5 mag during eight days of observations; 339 of these are previously unknown variables. We tentatively classify the 560 variables as 285 eclipsing binaries (EW, EB, EA), 27 pulsating variable stars ($\delta$~Scuti, $\gamma$~Doradus, $\delta$~Cephei variable and RR Lyrae stars) and 248 other types of variables (unclassified periodic, multi-periodic and aperiodic variable stars). Among the eclipsing binaries, 34 show O'Connell effects. One of the aperiodic variables shows a plateau light curve and another one shows a secondary maximum after peak brightness. We also detected a complex binary system with RS CVn-like light curve morphology; this object is being followed-up spectroscopically using the Gemini South telescope.

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.

GAIA leads us to step into a new era with a high astrometry precision of 10 uas. Under such a precision, astrometry will play important roles in detecting and characterizing exoplanets. Specially, we can identify planet pairs in mean motion resonances(MMRs) via astrometry, which constrains the formation and evolution of planetary systems. In accordance with observations, we consider two Jupiters or two super-Earths systems in 1:2, 2:3 and 3:4 MMRs. Our simulations show the false alarm probabilities(FAPs) of a third planet are extremely small while the real two planets can be good fitted with signal-to-noise ratio(SNR)> 3. The probability of reconstructing a resonant system is related with the eccentricities and resonance intensity. Generally, when SNR >= 10, if eccentricities of both planets are larger than 0.01 and the resonance is quite strong, the probabilities to reconstruct the planet pair in MMRs >= 80%. Jupiter pairs in MMRs are reconstructed more easily than super-Earth pairs with similar SNR when we consider the dynamical stability. FAPs are also calculated when we detect planet pairs in or near MMRs. FAPs for 1:2 MMR are largest, i.e., FAPs > 15% when SNR <= 10. Extrapolating from the Kepler planet pairs near MMRs and assuming SNR to be 3, we will discover and reconstruct a few tens of Jupiter pairs and hundreds of super-Earth pairs in 2:3 and 1:2 MMRs within 30 pc. We also compare the differences between even and uneven data cadence and find that planets are better measured with more uniform phase coverage.

Context. Infrared dark clouds (IRDCs) are ubiquitous in the Milky Way, yet they play a crucial role in breeding newly-formed stars. Aims. With the aim of further understanding the dynamics, chemistry, and evolution of IRDCs, we carried out multi-wavelength observations on a small sample. Methods. We performed new observations with the IRAM 30 m and CSO 10.4 m telescopes, with tracers ${\rm HCO^+}$, HCN, ${\rm N_2H^+}$, ${\rm C^{18}O}$, DCO$^+$, SiO, and DCN toward six IRDCs G031.97+00.07, G033.69-00.01, G034.43+00.24, G035.39-00.33, G038.95-00.47, and G053.11+00.05. Results. We investigated 44 cores including 37 cores reported in previous work and seven newly-identified cores. Toward the dense cores, we detected 6 DCO$^+$, and 5 DCN lines. Using pixel-by-pixel spectral energy distribution (SED) fits of the $\textit{Herschel}$ 70 to 500 $\mu$m, we obtained dust temperature and column density distributions of the IRDCs. We found that ${\rm N_2H^+}$ emission has a strong correlation with the dust temperature and column density distributions, while ${\rm C^{18}O}$ showed the weakest correlation. It is suggested that ${\rm N_2H^+}$ is indeed a good tracer in very dense conditions, but ${\rm C^{18}O}$ is an unreliable one, as it has a relatively low critical density and is vulnerable to freezing-out onto the surface of cold dust grains. The dynamics within IRDCs are active, with infall, outflow, and collapse; the spectra are abundant especially in deuterium species. Conclusions. We observe many blueshifted and redshifted profiles, respectively, with ${\rm HCO^+}$ and ${\rm C^{18}O}$ toward the same core. This case can be well explained by model "envelope expansion with core collapse (EECC)".

We present results of a dark matter search performed with a 0.6 kg day exposure of the DAMIC experiment at the SNOLAB underground laboratory. We measure the energy spectrum of ionization events in the bulk silicon of charge-coupled devices down to a signal of 60 eV electron equivalent. The data are consistent with radiogenic backgrounds, and constraints on the spin-independent WIMP-nucleon elastic-scattering cross section are accordingly placed. A region of parameter space relevant to the potential signal from the CDMS-II Si experiment is excluded using the same target for the first time. This result obtained with a limited exposure demonstrates the potential to explore the low-mass WIMP region (<10 GeV/$c^{2}$) of the upcoming DAMIC100, a 100 g detector currently being installed in SNOLAB.

We report a measurement of the ionization efficiency of silicon nuclei recoiling with sub-keV kinetic energy in the bulk silicon of a charge-coupled device (CCD). Nuclear recoils are produced by low-energy neutrons ($<$24 keV) from a $^{124}$Sb-$^{9}$Be photoneutron source, and their ionization signal is measured down to 60 eV electron equivalent. This energy range, previously unexplored, is relevant for the detection of low-mass dark matter particles. The measured efficiency is found to deviate from the extrapolation to low energies of the Lindhard model. This measurement also demonstrates the sensitivity to nuclear recoils of CCDs employed by DAMIC, a dark matter direct detection experiment located in the SNOLAB underground laboratory.

For a general understanding of the physics involved in the star formation process, measurements of physical parameters such as temperature and density are indispensable. The chemical and physical properties of dense clumps of molecular clouds are strongly affected by the kinetic temperature. Therefore, this parameter is essential for a better understanding of the interstellar medium. Formaldehyde, a molecule which traces the entire dense molecular gas, appears to be the most reliable tracer to directly measure the gas kinetic temperature.We aim to determine the kinetic temperature with spectral lines from formaldehyde and to compare the results with those obtained from ammonia lines for a large number of massive clumps.Three 218 GHz transitions (JKAKC=303-202, 322-221, and 321-220) of para-H2CO were observed with the 15m James Clerk Maxwell Telescope (JCMT) toward 30 massive clumps of the Galactic disk at various stages of high-mass star formation. Using the RADEX non-LTE model, we derive the gas kinetic temperature modeling the measured para-H2CO 322-221/303-202and 321-220/303-202 ratios. The gas kinetic temperatures derived from the para-H2CO (321-220/303-202) line ratios range from 30 to 61 K with an average of 46 K. A comparison of kinetic temperature derived from para-H2CO, NH3, and the dust emission indicates that in many cases para-H2CO traces a similar kinetic temperature to the NH3 (2,2)/(1,1) transitions and the dust associated with the HII regions. Distinctly higher temperatures are probed by para-H2CO in the clumps associated with outflows/shocks. Kinetic temperatures obtained from para-H2CO trace turbulence to a higher degree than NH3 (2,2)/(1,1) in the massive clumps. The non-thermal velocity dispersions of para-H2CO lines are positively correlated with the gas kinetic temperature. The massive clumps are significantly influenced by supersonic non-thermal motions.

Many exoplanets have been found in orbits close to their host stars and thus they are subject to the effects of photo-evaporation. Previous studies have shown that a large portion of exoplanets detected by the Kepler mission have been significantly eroded by photo-evaporation. In this paper, we numerically study the effects of photo-evaporation on the orbital evolution of a hypothesized moon system around a planet. We find that photo-evaporation is crucial to the stability of the moon system. Photo-evaporation can erode the atmosphere of the planet thus leading to significant mass loss. As the planet loses mass, its Hill radius shrinks and its moons increase their orbital semi-major axes and eccentricities. When some moons approach their critical semi-major axes, global instability of the moon system would be triggered, which usually ends up with two, one or even zero surviving moons. Some lost moons could escape from the moon system to become a new planet orbiting the star or run away further to become a free-floating object in the Galaxy. Given the destructive role of photo-evaporation, we speculate that exo-moons are less common for close-in planets ($<0.1$ AU), especially those around M-type stars, because they are more X-ray luminous and thus enhancing photo-evaporation. The lessons we learn in this study may be helpful for the target selection of on-going/future exomoon searching programs.

The nearly circular (mean eccentricity <e>~0.06) and coplanar (mean mutual inclination <i>~3 deg) orbits of the Solar System planets motivated Kant and Laplace to put forth the hypothesis that planets are formed in disks, which has developed into the widely accepted theory of planet formation. Surprisingly, the first several hundred extrasolar planets (mostly Jovian) discovered using the Radial Velocity (RV) technique are commonly on eccentric orbits (<e> ~ 0.3). This raises a fundamental question: Are the Solar System and its formation special? The Kepler mission has found thousands of transiting planets dominated by sub-Neptunes, but most of their orbital eccentricities remain unknown. By using the precise spectroscopic host star parameters from the LAMOST observations, we measure the eccentricity distributions for a large (698) and homogeneous Kepler planet sample with transit duration statistics. Nearly half of the planets are in systems with single transiting planets (singles), while the other half are multiple-transiting planets (multiples). We find an eccentricity dichotomy: on average, Kepler singles are on eccentric orbits with <e>~0.3, while the multiples are on nearly circular (<e> = 0.04^{+0.03}_{-0.04}) and coplanar (<i> = 1.4^{+0.8}_{-1.1} deg) orbits similar to the Solar System planets. Our results are consistent with previous studies of smaller samples and individual systems. We also show that Kepler multiples and solar system objects follow a common relation <e>~(1-2)x<i> between mean eccentricities and mutual inclinations. The prevalence of circular orbits and the common relation may imply that the solar system is not so atypical in the galaxy after all.

The Chinese Small Telescope ARray (CSTAR) is the first Chinese astronomical instrument placed in Antarctica. It is a group of four identical, fully automatic $14.5\,\rm{cm}$ telescopes, with an field of view (FOV) of $20\,\rm{deg^2}$ centered on the South Celestial Pole. Placed at Antarctic Dome A, CSTAR is designed to provide high-cadence photometry for site monitoring and variable sources detection. During the 2008 observing season, CSTAR has taken high-precision photometric data for 18,145 stars around the South Celestial Pole. At $i\,=\,7.5$ and $12$, the photometric precision reaches $\sim 8$ mmag and $\sim 30$ mmag with a cadence of 20s or 30s, respectively. Using robust detection method, we have found 15 stellar flares on 13 sources, including two classified variables. We have also found a linear relation between the decay times and the total durations of the stellar flares. The details of all detected flares along with their stellar properties are presented in this work.

We determine the proportions of two mixed crater populations distinguishable by size distributions on the Moon. A "multiple power-law" model is built to formulate crater size distribution $N(D) \propto D^{-\alpha}$ whose slope $\alpha$ varies with crater diameter $D$. Fitted size distribution of lunar highland craters characterized by $\alpha = 1.17 \pm 0.04$, $1.88 \pm 0.07$, $3.17 \pm 0.10$ and $1.40 \pm 0.15$ for consecutive $D$ intervals divided by 49, 120 and 251 km and that of lunar Class 1 craters with a single slope $\alpha = 1.96 \pm 0.14$, are taken as Population 1 and 2 crater size distribution respectively, whose sum is then fitted to the size distribution of global lunar craters with $D$ between 10 and 100 km. Estimated crater densities of Population 1 and 2 are $44 \times 10^{-5}$ and $5 \times 10^{-5}$ km$^{-2}$ respectively, leading to the proportion of the latter $10 \%$. The results underlines the need for considering the Population 1 craters and the relevant impactors, the primordial main-belt asteroids which dominated the late heavy bombardment.

We formulate the lunar cratering distribution and verify the cratering asymmetries generated by the main-belt asteroids (MBAs) as well as the near-Earth objects (NEOs). Based on a planar model that excludes the terrestrial and lunar gravitations on the impactors and assuming the impactor encounter speed with Earth $v_{\rm{enc}}$ is higher than the lunar orbital speed $v_{\rm{M}}$, we rigorously integrated the lunar cratering distribution, and derived its approximation to the first order of $v_{\rm{M}}/v_{\rm{enc}}$. Numerical simulations of lunar bombardment by the MBAs during the late heavy bombardment were performed with an Earth-Moon distance $a_{\rm{M}}$ = 20--60 Earth radii in five cases. The analytical model directly proves the existence of a leading/trailing asymmetry and the absence of near/far asymmetry. The approximate form of the leading/trailing asymmetry is $(1 + A_1 \cos\beta)$, which decreases as the apex distance $\beta$ increases. The numerical simulations show evidence of a pole/equator asymmetry as well as the leading/trailing asymmetry, and the former is empirically described as $(1 + A_2 \cos2\varphi)$, which decreases as the latitude modulus $|\varphi|$ increases. The amplitudes $A_{1,2}$ are reliable measurements of asymmetries. Our analysis explicitly indicates the quantitative relations between cratering distribution and bombardment conditions (impactor properties and the lunar orbital status) like $A_1 \propto v_{\rm{M}}/v_{\rm{enc}}$, resulting in a method for reproducing the bombardment conditions through measuring the asymmetry. Mutual confirmation between analytical model and numerical simulations is found in terms of the cratering distribution and its variation with $a_{\rm{M}}$. Estimates of $A_1$ for crater density distributions generated by the MBAs and the NEOs are 0.101--0.159 and 0.117, respectively.

In this work, we aim to characterise high-mass clumps with infall motions. We selected 327 clumps from the Millimetre Astronomy Legacy Team 90-GHz (MALT90) survey, and identified 100 infall candidates. Combined with the results of He et al. (2015), we obtained a sample of 732 high-mass clumps, including 231 massive infall candidates and 501 clumps where infall is not detected. Objects in our sample were classified as pre-stellar, proto-stellar, HII or photo-dissociation region (PDR). The detection rates of the infall candidates in the pre-stellar, proto-stellar, HII and PDR stages are 41.2%, 36.6%, 30.6% and 12.7%, respectively. The infall candidates have a higher H$_{2}$ column density and volume density compared with the clumps where infall is not detected at every stage. For the infall candidates, the median values of the infall rates at the pre-stellar, proto-stellar, HII and PDR stages are 2.6$\times$10$^{-3}$, 7.0$\times$10$^{-3}$, 6.5$\times$10$^{-3}$ and 5.5$\times$10$^{-3}$ M$_\odot$ yr$^{-1}$, respectively. These values indicate that infall candidates at later evolutionary stages are still accumulating material efficiently. It is interesting to find that both infall candidates and clumps where infall is not detected show a clear trend of increasing mass from the pre-stellar to proto-stellar, and to the HII stages. The power indices of the clump mass function (ClMF) are 2.04$\pm$0.16 and 2.17$\pm$0.31 for the infall candidates and clumps where infall is not detected, respectively, which agree well with the power index of the stellar initial mass function (2.35) and the cold Planck cores (2.0).

In the third Fermi catalogue (3FGL) there are sixteen gamma-ray globular clusters. After analyzing the recent released Pass 8 data of Fermi Large Area Telescope (LAT), we report the discovery of significant gamma-ray emission from M 15 and NGC 6397, confirm that NGC 5904 is a gamma-ray emitter and find evidence of gamma-ray emission from NGC 6218 and NGC 6139. Inter- estingly, in the globular clusters M 15, NGC 6397 and NGC 5904, millisecond pulsars (MSPs) have been found in radio or X-rays, which are strongly in support of the MSP origin of the gamma-ray emission. However, due to the relative low luminosity of the gamma-ray emission we do not find any evidence for the gamma-ray pulsation or flux variability of these sources.

The DAMIC experiment uses fully depleted, high resistivity CCDs to search for dark matter particles. With an energy threshold $\sim$50 eV$_{ee}$, and excellent energy and spatial resolutions, the DAMIC CCDs are well-suited to identify and suppress radioactive backgrounds, having an unrivaled sensitivity to WIMPs with masses $<$6 GeV/$c^2$. Early results motivated the construction of a 100 g detector, DAMIC100, currently being installed at SNOLAB. This contribution discusses the installation progress, new calibration efforts near the threshold, a preliminary result with 2014 data, and the prospects for physics results after one year of data taking.

We present the discovery of Balmer line absorption from H$\alpha$ to H$\gamma$ in an iron low-ionization broad absorption line (FeLoBAL) quasar SDSS J152350.42+391405.2 (hereafter J1523), by the quasi-simultaneous optical and near-infrared spectroscopy. The Balmer line absorption is at $z_{absor}$ = 0.6039 +/-0.0021 and blueshifted by v=10,353 km/s with respect to the Balmer emission lines. All Balmer BALs have uniform absorption profile with the widths of $\Delta$ v ~12,000 km/s. We also found the absorption trough in He 1* $\lambda$10830 with the same velocity and width in the H-band TripleSpec spectrum of J1523. This object is only the tenth active galactic nucleus known to exhibit non-stellar Balmer absorption, and also the case with the highest velocity and broadest Balmer absorption lines which have ever been found. A CLOUDY analysis shows that the absorbers require an gas density of $log_{10} n_ e (cm^{-3})=9$ and an ionization parameter of $log_{10} U=-1.0$. They locate at a distance of ~0.2 pc from the central ionizing source which is slightly farther than that of BELRs. Furthermore, J1523 is one of the brightest Balmer-BAL quasar ever reported, with unique iron absorption variations, making it as the most promising candidate for follow up high-resolution spectroscopy, multi-band observations, and long-term monitoring.

Resolution, usually defined by the Rayleigh criterion or the Full Width at Half Maximum of a Point Spread Function, is a basic property of an image. Here, we present a new statistical definition of image resolution based on the cross-correlation properties of the pixels in an image. It is shown that the new definition of image resolution depends not only on the PSF of an imaging device, but also on the signal-to-noise ratio of the data and on the structures of an object. In an image, the resolution does not have to be uniform. Our new definition is also suitable for the interpretation of the result of a deconvolution. We illustrate this, in this paper, with a Wiener deconvolution. It is found that weak structures can be extracted from low signal-to-noise ratio data, but with low resolution; a high-resolution image was obtained from high signal-to-noise ratio data after a Wiener deconvolution. The new definition can also be used to compare various deconvolution algorithms on their processing effects, such as resolution, sensitivity and sidelobe level, etc.

The DAMIC (Dark Matter in CCDs) experiment uses high resistivity, scientific grade CCDs to search for dark matter. The CCD's low electronic noise allows an unprecedently low energy threshold of a few tens of eV that make it possible to detect silicon recoils resulting from interactions of low mass WIMPs. In addition the CCD's high spatial resolution and the excellent energy response results in very effective background identification techniques. The experiment has a unique sensitivity to dark matter particles with masses below 10 GeV/c$^2$. Previous results have demonstrated the potential of this technology, motivating the construction of DAMIC100, a 100 grams silicon target detector currently being installed at SNOLAB. In this contribution, the mode of operation and unique imaging capabilities of the CCDs, and how they may be exploited to characterize and suppress backgrounds will be discussed, as well as physics results after one year of data taking.

We present measurements of radioactive contamination in the high-resistivity silicon charge-coupled devices (CCDs) used by the DAMIC experiment to search for dark matter particles. Novel analysis methods, which exploit the unique spatial resolution of CCDs, were developed to identify $\alpha$ and $\beta$ particles. Uranium and thorium contamination in the CCD bulk was measured through $\alpha$ spectroscopy, with an upper limit on the $^{238}$U ($^{232}$Th) decay rate of 5 (15) kg$^{-1}$ d$^{-1}$ at 95% CL. We also searched for pairs of spatially correlated electron tracks separated in time by up to tens of days, as expected from $^{32}$Si-$^{32}$P or $^{210}$Pb-$^{210}$Bi sequences of $\beta$ decays. The decay rate of $^{32}$Si was found to be $80^{+110}_{-65}$ kg$^{-1}$ d$^{-1}$ (95% CI). An upper limit of $\sim$35 kg$^{-1}$ d$^{-1}$ (95% CL) on the $^{210}$Pb decay rate was obtained independently by $\alpha$ spectroscopy and the $\beta$ decay sequence search. These levels of radioactive contamination are sufficiently low for the successful operation of CCDs in the forthcoming 100 g DAMIC detector.

The Chinese Small Telescope ARray (CSTAR) is the first telescope facility built at Dome A, Antarctica. During the 2008 observing season, the installation provided long-baseline and high-cadence photometric observations in the i-band for 18,145 targets within 20 deg2 CSTAR field around the South Celestial Pole for the purpose of monitoring the astronomical observing quality of Dome A and detecting various types of photometric variability. Using sensitive and robust detection methods, we discover 274 potential variables from this data set, 83 of which are new discoveries. We characterize most of them, providing the periods, amplitudes and classes of variability. The catalog of all these variables is presented along with the discussion of their statistical properties.

The Chinese Small Telescope ARray (CSTAR) has observed an area around the Celestial South Pole at Dome A since 2008. About $20,000$ light curves in the i band were obtained lasting from March to July, 2008. The photometric precision achieves about 4 mmag at i = 7.5 and 20 mmag at i = 12 within a 30 s exposure time. These light curves are analyzed using Lomb--Scargle, Phase Dispersion Minimization, and Box Least Squares methods to search for periodic signals. False positives may appear as a variable signature caused by contaminating stars and the observation mode of CSTAR. Therefore the period and position of each variable candidate are checked to eliminate false positives. Eclipsing binaries are removed by visual inspection, frequency spectrum analysis and locally linear embedding technique. We identify 53 eclipsing binaries in the field of view of CSTAR, containing 24 detached binaries, 8 semi-detached binaries, 18 contact binaries, and 3 ellipsoidal variables. To derive the parameters of these binaries, we use the Eclipsing Binaries via Artificial Intelligence (EBAI) method. The primary and the secondary eclipse timing variations (ETVs) for semi-detached and contact systems are analyzed. Correlated primary and secondary ETVs confirmed by false alarm tests may indicate an unseen perturbing companion. Through ETV analysis, we identify two triple systems (CSTAR J084612.64-883342.9 and CSTAR J220502.55-895206.7). The orbital parameters of the third body in CSTAR J220502.55-895206.7 are derived using a simple dynamical model.

The Hard X-ray Modulation Telescope (HXMT) will perform an all-sky survey in hard X-ray band as well as deep imaging of a series of small sky regions. We expect various compact objects to be detected in these imaging observations. Point source detection performance of HXMT imaging observation depends not only on the instrument but also on its data analysis since images are reconstructed from HXMT observed data with numeric methods. Denoising technique plays an import part in HXMT imaging data analysis pipeline alongside with demodulation and source detection. In this paper we have implemented several methods for denoising HXMT data and evaluated the point source detection performances in terms of sensitivities and location accuracies. The results show that direct demodulation with 1-fold cross correlation should be the default reconstruction and regularization methods, although both sensitivity and location accuracy could be further imporved by selecting and tuning numerical methods in data analysis of HXMT imaging observations.

Massive star-forming regions with observed infall motions are good sites for studying the birth of massive stars. In this paper, 405 compact sources have been extracted from the APEX Telescope Large Area Survey of the Galaxy (ATLASGAL) compact sources that also have been observed in the Millimetre Astronomy Legacy Team 90 GHz (MALT90) survey during Years 1 and 2. These observations are complemented with Spitzer GLIMPSE/MIPSGAL mid-IR survey data to help classify the elected star-forming clumps into three evolutionary stages: pre-stellar, proto-stellar and UCHII regions. The results suggest that 0.05 g cm$^{-2}$ is a reliable empirical lower bound for the clump surface densities required for massive-star formation to occur. The optically thick HCO$^{+}$(1-0) and HNC(1-0) lines, as well as the optically thin N$_{2}$H$^{+}$(1-0) line were used to search for infall motions toward these sources. By analyzing the asymmetries of the optically thick HCO$^{+}$(1-0) and HNC(1-0) lines and the mapping observations of HCO$^{+}$(1-0), a total of 131 reliable infall candidates have been identified. The HCO$^{+}$(1-0) line shows the highest occurrence of obvious asymmetric features, suggesting that it may be a better infall motion tracer than other lines such as HNC(1-0). The detection rates of infall candidates toward pre-stellar, proto-stellar and UCHII clumps are 0.3452, 0.3861 and 0.2152, respectively. The relatively high detection rate of infall candidates toward UCHII clumps indicates that many UCHII regions are still accreting matter. The peak column densities and masses of the infall candidates, in general, display a increasing trend with progressing evolutionary stages. However, the rough estimates of the mass infall rate show no obvious variation with evolutionary stage.

Globular clusters (GCs) are the oldest stellar system in the Galaxy, in which the millisecond pulsars are widely believed to be the only steady {\gamma}-ray emitters. So far 9 {\gamma}-ray GCs have been identified and a few candidates such as 2MS-GC01 and IC 1257 have been suggested. In this work, after analyzing the publicly-available Fermi-LAT data we confirm the significant {\gamma}-ray emission from 2MS-GC01 and IC 1257 and report the discovery for {\gamma}-ray emission from NGC 5904 and NGC 6656 within their tidal radii. Also a strong evidence of significant {\gamma}-ray emission is found from FSR 1735. From the observed {\gamma}-ray luminosities, the numbers of MSPs that are expected to be present in these GCs are estimated.

Analysis of the transit timing variations (TTVs) of candidate pairs near mean-motion resonances (MMRs) is an effective method to confirm planets. Hitherto, 68 planets in 34 multi-planet systems have been confirmed via TTVs. We analyze the TTVs of all candidates from the most recent Kepler data with a time span as long as about 1350 days (Q0-Q15). The anti-correlations of TTV signals and the mass upper limits of candidate pairs in the same system are calculated, using an improved method suitable for long-period TTVs. If the false alarm probability (FAP) of a candidate pair is less than 10^{-3} and the mass upper limit for each candidate is less than 13 M_J, we confirm them as planets in the same system. Finally, 8 planets in 4 multi-planet systems are confirmed via analysis of their TTVs. All of the 4 planet pairs are near first-order MMRs, including KOI-2672 near 2:1 MMR, KOI-1236, KOI-1563 and KOI-2038 near 3:2 MMR. Four planets have relatively long orbital periods (> 35 day). KOI-2672.01 has a period of 88.51658 days and a fit mass of 17 M_\oplus. It is the longest-period planet confirmed, to date, near first-order MMR via TTVs.

We present the optical spectroscopic observations of X Per from 1999 to 2013 with the 2.16m telescope at Xinglong Station and the 2.4m telescope at Lijiang station, National Astronomical Observatories of China. Combining with the public optical photometric data, we find certain epochs of anti-correlations between the optical brightness and the intensity of the H{\alpha} and HeI 6678 lines, which may be attributed to the mass ejections from the Be star. Alternative explanations are however also possible. The variability of FeII 6317 line in the spectra of X Per might be also caused by the shocked waves formed after the mass ejections from the Be star. The X-ray activities of the system might also be connected with the mass ejection events from the Be star. When the ejected materials were transported from the surface of the Be star to the orbit of neutron star, an X-ray flare could be observed in its X-ray light curves. We use the neutron star as a probe to constrain the motion of the ejected material in the circumstellar disk. With the diffusion time of the ejected material from the surface of Be star to the orbit of neutron star, the viscosity parameter {\alpha} of the circumstellar disk is estimated to be 0.39 and 0.28 for the different time, indicating that the disk around Be star may be truncated by the neutron star at the 2:1 resonance radius and Type I X-ray outburst is unlikely to be observed in X Per.

Spherical harmonic moments are well-suited for capturing anisotropy at any scale in the flux of cosmic rays. An unambiguous measurement of the full set of spherical harmonic coefficients requires full-sky coverage. This can be achieved by combining data from observatories located in both the northern and southern hemispheres. To this end, a joint analysis using data recorded at the Telescope Array and the Pierre Auger Observatory above $10^{19}$ eV is presented in this work. The resulting multipolar expansion of the flux of cosmic rays allows us to perform a series of anisotropy searches, and in particular to report on the angular power spectrum of cosmic rays above $10^{19}$ eV. No significant deviation from isotropic expectations is found throughout the analyses performed. Upper limits on the amplitudes of the dipole and quadrupole moments are derived as a function of the direction in the sky, varying between 7% and 13% for the dipole and between 7% and 10% for a symmetric quadrupole.

We describe the method devised to reconstruct inclined cosmic-ray air showers with zenith angles greater than $60^\circ$ detected with the surface array of the Pierre Auger Observatory. The measured signals at the ground level are fitted to muon density distributions predicted with atmospheric cascade models to obtain the relative shower size as an overall normalization parameter. The method is evaluated using simulated showers to test its performance. The energy of the cosmic rays is calibrated using a sub-sample of events reconstructed with both the fluorescence and surface array techniques. The reconstruction method described here provides the basis of complementary analyses including an independent measurement of the energy spectrum of ultra-high energy cosmic rays using very inclined events collected by the Pierre Auger Observatory.

Transit times around single stars can be described well by a linear ephemeris. However, transit times in circumbinary systems are influenced both by the gravitational perturbations and the orbital phase variations of the central binary star. Adopting a coplanar analog of Kepler-16 as an example, we find that circumbinary planets can transit the same star more than once during a single planetary orbit, a phenomenon we call "tight transits". In certain geometric, the projected orbital velocity of the planet and the secondary star can approach zero and change sign, resulting in very long transits and/or 2-3 transits during a single binary orbit. Whether tight transits are possible for a particular system depends primarily on the binary mass ratio and the orbital architecture of both the binary and the planet. We derive a time-dependent criterion to judge when tight transits are possible for any circumbinary system. These results are verified with full dynamical integrations that also reveal other tight transit characteristics, i.e., the transit durations and the intervals between tight transits. For the seven currently known circumbinary systems, we estimate these critical parameters both analytically and numerically. Due to the mutual inclination between the planet and the binary, tight transits can only occur across the less massive star B in Kepler-16, -34, -35, and -47 (for both planets). The long-term average frequency of tight transits (compared to typical transits) for Kepler-16, -34, and -35 are estimated to be several percent. Using full numerical integrations, the next tight transit for each system is predicted and the soonest example appears to be Kepler-47b and -47c, which are likely to have tight transits before 2025. These unique and valuable events often deserve special observational scrutiny.

A transformed auto-correlation method is presented here, where a received signal is transformed based on a priori reflecting model, and then the transformed signal is cross-correlated to its original one. If the model is correct, after transformation, the reflected signal will be coherent to the transmitted signal, with zero delay. A map of transformed auto-correlation function with zero delay can be generated in a given parametric space. The significant peaks in the map may indicate the possible reflectors nearby the central transmitter. The true values of the parameters of reflectors can be estimated at the same time.

We introduce the fully-depleted charge-coupled device (CCD) as a particle detector. We demonstrate its low energy threshold operation, capable of detecting ionizing energy depositions in a single pixel down to 50 eVee. We present results of energy calibrations from 0.3 keVee to 60 keVee, showing that the CCD is a fully active detector with uniform energy response throughout the silicon target, good resolution (Fano ~0.16), and remarkable linear response to electron energy depositions. We show the capability of the CCD to localize the depth of particle interactions within the silicon target. We discuss the mode of operation and unique imaging capabilities of the CCD, and how they may be exploited to characterize and suppress backgrounds. We present the first results from the deployment of 250 um thick CCDs in SNOLAB, a prototype for the upcoming DAMIC100. DAMIC100 will have a target mass of 0.1 kg and should be able to directly test the CDMS-Si signal within a year of operation.

The Hard X-ray Modulation Telescope (HXMT) is a Chinese space telescope mission. It is scheduled for launch in 2015. The telescope will perform an all-sky survey in hard X-ray band (1 - 250 keV), a series of deep imaging observations of small sky regions as well as pointed observations. In this work we present a conceptual method to reconstruct light curves from HXMT imaging observation directly, in order to monitor time-varying objects such as GRB, AXP and SGR in hard X-ray band with HXMT imaging observations.

A flux of neutrons from an astrophysical source in the Galaxy can be detected in the Pierre Auger Observatory as an excess of cosmic-ray air showers arriving from the direction of the source. To avoid the statistical penalty for making many trials, classes of objects are tested in combinations as nine "target sets", in addition to the search for a neutron flux from the Galactic Center or from the Galactic Plane. Within a target set, each candidate source is weighted in proportion to its electromagnetic flux, its exposure to the Auger Observatory, and its flux attenuation factor due to neutron decay. These searches do not find evidence for a neutron flux from any class of candidate sources. Tabulated results give the combined p-value for each class, with and without the weights, and also the flux upper limit for the most significant candidate source within each class. These limits on fluxes of neutrons significantly constrain models of EeV proton emission from non-transient discrete sources in the Galaxy.

Measurements of air showers made using the hybrid technique developed with the fluorescence and surface detectors of the Pierre Auger Observatory allow a sensitive search for point sources of EeV photons anywhere in the exposed sky. A multivariate analysis reduces the background of hadronic cosmic rays. The search is sensitive to a declination band from -85{\deg} to +20{\deg}, in an energy range from 10^17.3 eV to 10^18.5 eV. No photon point source has been detected. An upper limit on the photon flux has been derived for every direction. The mean value of the energy flux limit that results from this, assuming a photon spectral index of -2, is 0.06 eV cm^-2 s^-1, and no celestial direction exceeds 0.25 eV cm^-2 s^-1. These upper limits constrain scenarios in which EeV cosmic ray protons are emitted by non-transient sources in the Galaxy.

Aims. The process of gravitational fragmentation in the L1482 molecular filament of the California molecular cloud is studied by combining several complementary observations and physical estimates. We investigate the kinematic and dynamical states of this molecular filament and physical properties of several dozens of dense molecular clumps embedded therein. Methods. We present and compare molecular line emission observations of the J=2--1 and J=3--2 transitions of 12CO in this molecular complex, using the KOSMA 3-meter telescope. These observations are complemented with archival data observations and analyses of the 13CO J=1--0 emission obtained at the Purple Mountain Observatory 13.7-meter radio telescope at Delingha Station in QingHai Province of west China, as well as infrared emission maps from the Herschel Space Telescope online archive, obtained with the SPIRE and PACS cameras. Comparison of these complementary datasets allow for a comprehensive multi-wavelength analysis of the L1482 molecular filament. Results. We have identified 23 clumps along the molecular filament L1482 in the California molecular cloud. All these molecular clumps show supersonic non-thermal gas motions. While surprisingly similar in mass and size to the much better known Orion molecular cloud, the formation rate of high-mass stars appears to be suppressed in the California molecular cloud relative to that in the Orion molecular cloud based on the mass-radius threshold derived from the static Bonnor Ebert sphere. Our analysis suggests that these molecular filaments are thermally supercritical and molecular clumps may form by gravitational fragmentation along the filament. Instead of being static, these molecular clumps are most likely in processes of dynamic evolution.

We observed the H2CO(110-111) absorption lines and H110{\alpha} radio recombination lines (RRL) toward 180 NH3 sources using the Nanshan 25-m radio telescope. In our observation, 138 sources were found to have H2CO lines and 36 have H110{\alpha} RRLs. Among the 138 detected H2CO sources, 38 sources were first detected. The detection rates of H2CO have a better correlation with extinction than with background continuum radiation. Line center velocities of H2CO and NH3 agree well. The line width ratios of H2CO and NH3 are generally larger than 1 and are similar to that of 13CO. The correlation between column densities of H2CO and extinction is better than that between NH3 and extinction. These line width relation and column density relation indicate H2CO is distributed on a larger scale than that of NH3, being similar to the regions of 13CO. The abundance ratios between NH3 and H2CO were found to be different in local clouds and other clouds.

We present large scale observations of C18O(1-0) towards four massive star forming regions for MON R2, S156, DR17/L906 and M17/M18. The transitions of H2CO(110-111), C18O(1-0) and 6 cm continuum were compared towards the four regions. Analysis of observation and Non--LTE model results shows that the brightness temperature of the formaldehyde absorption line is strongest in background continuum temperature range of about 3 - 8 K. The excitation of the H2CO absorption line is affected by strong background continuum emission. From the comparison of H2CO and C18O maps, we found that the extent of H2CO absorption is broader than that of C18O emission in the four regions. Except for the DR17 region, the H2CO absorption maximum is located at the same position with the C18O peak. The good correlation between intensities and widths of H2CO absorption and C18O emission lines indicate that the H2CO absorption line can trace dense and warm regions of the molecular cloud. Finding that N(H2CO) was well correlated with N(C18O) in the four regions and that the average column density ratio is <H2CO/N(C18O)> ~ 0.03.

We search for systems hosting eclipsing discs using a complete sample of eclipsing binaries (EBs); those previously identified in the third phase of the Optical Gravitational Lensing Experiment (OGLE-III). Within a subsample of 2,823 high-cadence, high-photometric precision and large eclipsing depth detached EBs previously identified in the Large Magellanic Cloud (LMC), we find that the skewness and kurtosis of the light curves magnitude distribution within the primary eclipse can distinguish EBs hosting a disc from those without. Two systems with previously identified eclipsing discs (OGLE-LMC-ECL-11893 and OGLE-LMC-ECL-17782) are identified with near zero skewness ($|S|<0.5$) and positive kurtosis. No additional eclipsing disc systems were found in the OGLE-III LMC, Small Magellanic Cloud (SMC) or Galactic Disc (GD) EB light curves. We estimate that the fraction of detached near main-sequence LMC EBs (which have a primary with an $I$-band magnitude brighter than $\simeq 19\,\rm{mag}$) that host a disc is approximately 1/1000. As circumstellar disc lifetimes are short, we expected to primarily find eclipsing discs around young stars. In addition, as there is more room for a disc in a widely separated binary and because a disk close to a luminous star would be above the dust sublimation temperature, we expected to primarily find eclipsing discs in long period binaries. However, OGLE-LMC-ECL-17782 is a 13.3 day period B star system with a transient and hot ($\sim 6000\,\rm{K}$, $\sim 0.1\,\rm{AU}$ radius) disc and Scott et al. (in prep.) estimate an age of $150\,\rm{Myr}$ for OGLE-LMC-ECL-11893. Both discs are unexpected in the EB sample and impel explanation.

The Chinese Small Telescope ARray (CSTAR) is a group of four identical, fully automated, static 14.5 cm telescopes. CSTAR is located at Dome A, Antarctica and covers 20 square degree of sky around the South Celestial Pole. The installation is designed to provide high-cadence photometry for the purpose of monitoring the quality of the astronomical observing conditions at Dome A and detecting transiting exoplanets. CSTAR has been operational since 2008, and has taken a rich and high-precision photometric data set of 10,690 stars. In the first observing season, we obtained 291,911 qualified science frames with 20-second integrations in the i-band. Photometric precision reaches about 4 mmag at 20-second cadence at i=7.5, and is about 20 mmag at i=12. Using robust detection methods, ten promising exoplanet candidates were found. Four of these were found to be giants using spectroscopic follow-up. All of these transit candidates are presented here along with the discussion of their detailed properties as well as the follow-up observations.

Intriguing work on observations of 4.83 GHz formaldehyde (H2CO) absorptions and 4.87 GHz H110a radio recombination lines (RRLs) towards 6.7 GHz methanol (CH3OH) maser sources is presented. Methanol masers provide ideal sites to probe the earliest stages of massive star formation, while 4.8 GHz formaldehyde absorptions are accurate probes of physical conditions in dense $(10^{3} - 10^{5} cm^{-3})$ and low temperature molecular clouds towards massive star forming regions. The work is aimed at studying feature similarities between the formaldehyde absorptions and the methanol masers so as to expand knowledge of events and physical conditions in massive star forming regions. A total of 176 methanol maser sources were observed for formaldehyde absorptions, and formaldehyde absorptions were detected 138 of them. 53 of the formaldehyde absorptions were newly detected. We noted a poor correlation between the methanol and formaldehyde intensities, an indication that the signals (though arise from about the same regions) are enhanced by different mechanisms. Our results show higher detection rates of the formaldehyde lines for sources with stronger methanol signals. The strongest formaldehyde absorptions were associated with IRAS sources and IRDCs that have developed HII regions, and that do not have EGOs.

We present the results of our recent study on the interactions between a giant planet and a self-gravitating gas disk. We investigate how the disk's self-gravity affects the gap formation process and the migration of the giant planet. Two series of 1-D and 2-D hydrodynamic simulations are performed. We select several surface densities and focus on the gravitationally stable region. To obtain more reliable gravity torques exerted on the planet, a refined treatment of disk's gravity is adopted in the vicinity of the planet. Our results indicate that the net effect of the disk's self-gravity on the gap formation process depends on the surface density of the disk. We notice that there are two critical values, \Sigma_I and \Sigma_II. When the surface density of the disk is lower than the first one, \Sigma_0 < \Sigma_I, the effect of self-gravity suppresses the formation of a gap. When \Sigma_0 > \Sigma_I, the self-gravity of the gas tends to benefit the gap formation process and enlarge the width/depth of the gap. According to our 1-D and 2-D simulations, we estimate the first critical surface density \Sigma_I \approx 0.8MMSN. This effect increases until the surface density reaches the second critical value \Sigma_II. When \Sigma_0 > \Sigma_II, the gravitational turbulence in the disk becomes dominant and the gap formation process is suppressed again. Our 2-D simulations show that this critical surface density is around 3.5MMSN. We also study the associated orbital evolution of a giant planet. Under the effect of the disk's self-gravity, the migration rate of the giant planet increases when the disk is dominated by gravitational turbulence. We show that the migration timescale associates with the effective viscosity and can be up to 10^4 yr.

DAMIC (Dark Matter in CCDs) is a novel dark matter experiment that has unique sensitivity to dark matter particles with masses below 10 GeV. Due to its low electronic readout noise (R.M.S. ~3 e-) this instrument is able to reach a detection threshold below 0.5 keV nuclear recoil energy, making the search for dark matter particles with low masses possible. We report on early results and experience gained from a detector that has been running at SNOLAB from Dec 2012. We also discuss the measured and expected backgrounds and present the plan for future detectors to be installed in 2014.

The Pierre Auger Observatory is the world's largest cosmic ray observatory. Our current exposure reaches nearly 40,000 km$^2$ str and provides us with an unprecedented quality data set. The performance and stability of the detectors and their enhancements are described. Data analyses have led to a number of major breakthroughs. Among these we discuss the energy spectrum and the searches for large-scale anisotropies. We present analyses of our X$_{max}$ data and show how it can be interpreted in terms of mass composition. We also describe some new analyses that extract mass sensitive parameters from the 100% duty cycle SD data. A coherent interpretation of all these recent results opens new directions. The consequences regarding the cosmic ray composition and the properties of UHECR sources are briefly discussed.