In recent work, two of us proposed a nonsingular Ekpyrotic cosmology making use of an S-brane which forms at the end of the phase of Ekpyrotic contraction. This S-Brane mediates a transition between contraction and expansion. Graviitational waves passing through the S-Brane acquire a roughly scale-invariant spectrum, and if the S-Brane has zero shear, then a roughly scale-invariant spectrum of cosmological perturbatiions results. Here, we study the production of gauge field fluctuations driven by the decay of the S-Brane, and we show that the reheating process via gauge field production will be efficient, leading to a radiation-dominated expanding phase.

Providing a possible connection between neutrino emission and gravitational-wave (GW) bursts is important to our understanding of the physical processes that occur when black holes or neutron stars merge. In the Daya Bay experiment, using data collected from December 2011 to August 2017, a search has been performed for electron-antineutrino signals coinciding with detected GW events, including GW150914, GW151012, GW151226, GW170104, GW170608, GW170814, and GW170817. We used three time windows of $\mathrm{\pm 10~s}$, $\mathrm{\pm 500~s}$, and $\mathrm{\pm 1000~s}$ relative to the occurrence of the GW events, and a neutrino energy range of 1.8 to 100 MeV to search for correlated neutrino candidates. The detected electron-antineutrino candidates are consistent with the expected background rates for all the three time windows. Assuming monochromatic spectra, we found upper limits (90% confidence level) on electron-antineutrino fluence of $(1.13~-~2.44) \times 10^{11}~\rm{cm^{-2}}$ at 5 MeV to $8.0 \times 10^{7}~\rm{cm^{-2}}$ at 100 MeV for the three time windows. Under the assumption of a Fermi-Dirac spectrum, the upper limits were found to be $(5.4~-~7.0)\times 10^{9}~\rm{cm^{-2}}$ for the three time windows.

Providing a possible connection between neutrino emission and gravitational-wave (GW) bursts is important to our understanding of the physical processes that occur when black holes or neutron stars merge. In the Daya Bay experiment, using data collected from December 2011 to August 2017, a search has been performed for electron-antineutrino signals coinciding with detected GW events, including GW150914, GW151012, GW151226, GW170104, GW170608, GW170814, and GW170817. We used three time windows of $\mathrm{\pm 10~s}$, $\mathrm{\pm 500~s}$, and $\mathrm{\pm 1000~s}$ relative to the occurrence of the GW events, and a neutrino energy range of 1.8 to 100 MeV to search for correlated neutrino candidates. The detected electron-antineutrino candidates are consistent with the expected background rates for all the three time windows. Assuming monochromatic spectra, we found upper limits (90% confidence level) on electron-antineutrino fluence of $(1.13~-~2.44) \times 10^{11}~\rm{cm^{-2}}$ at 5 MeV to $8.0 \times 10^{7}~\rm{cm^{-2}}$ at 100 MeV for the three time windows. Under the assumption of a Fermi-Dirac spectrum, the upper limits were found to be $(5.4~-~7.0)\times 10^{9}~\rm{cm^{-2}}$ for the three time windows.

Providing a possible connection between neutrino emission and gravitational-wave (GW) bursts is important to our understanding of the physical processes that occur when black holes or neutron stars merge. In the Daya Bay experiment, using data collected from December 2011 to August 2017, a search has been performed for electron-antineutrino signals coinciding with detected GW events, including GW150914, GW151012, GW151226, GW170104, GW170608, GW170814, and GW170817. We used three time windows of $\mathrm{\pm 10~s}$, $\mathrm{\pm 500~s}$, and $\mathrm{\pm 1000~s}$ relative to the occurrence of the GW events, and a neutrino energy range of 1.8 to 100 MeV to search for correlated neutrino candidates. The detected electron-antineutrino candidates are consistent with the expected background rates for all the three time windows. Assuming monochromatic spectra, we found upper limits (90% confidence level) on electron-antineutrino fluence of $(1.13~-~2.44) \times 10^{11}~\rm{cm^{-2}}$ at 5 MeV to $8.0 \times 10^{7}~\rm{cm^{-2}}$ at 100 MeV for the three time windows. Under the assumption of a Fermi-Dirac spectrum, the upper limits were found to be $(5.4~-~7.0)\times 10^{9}~\rm{cm^{-2}}$ for the three time windows.

Gravitational waves from coalescing neutron stars encode information about nuclear matter at extreme densities, inaccessible by laboratory experiments. The late inspiral is influenced by the presence of tides, which depend on the neutron star equation of state. Neutron star mergers are expected to often produce rapidly-rotating remnant neutron stars that emit gravitational waves. These will provide clues to the extremely hot post-merger environment. This signature of nuclear matter in gravitational waves contains most information in the 2-4 kHz frequency band, which is outside of the most sensitive band of current detectors. We present the design concept and science case for a neutron star extreme matter observatory (NEMO): a gravitational-wave interferometer optimized to study nuclear physics with merging neutron stars. The concept uses high circulating laser power, quantum squeezing and a detector topology specifically designed to achieve the high-frequency sensitivity necessary to probe nuclear matter using gravitational waves. Above one kHz, the proposed strain sensitivity is comparable to full third-generation detectors at a fraction of the cost. Such sensitivity changes expected event rates for detection of post-merger remnants from approximately one per few decades with two A+ detectors to a few per year, and potentially allows for the first gravitational-wave observations of supernovae, isolated neutron stars, and other exotica.

The microquasar V404 Cygni (also known as GS 2023+338) was previously reported to have possible GeV gamma-ray emission in two days during its 2015 outburst. In order to provide more detailed information at the high energy range for this black hole binary system, we conduct detailed analysis to the data obtained with the Large Area Telescope (LAT) onboard the Fermi Gamma Ray Space Telescope (Fermi). Both LAT database and source catalog used are the latest. From the analysis, we can not confirm the previous detection. Instead, we find one possible detection (~4 sigma) of the source at the end of the outburst in the time period of 2015 Aug. 17-19, and one convincing detection (~7 sigma) in 2016 Aug. 23-25. The latter shows that the gamma-ray emission of the source is soft with photon index $\Gamma\sim 2.9$, mostly detected below ~1.3 GeV with Fermi LAT. As gamma-ray emission from microquasars is likely associated with their jet activity, we discuss the results by comparing with those well studied cases, namely Cyg X-3 and Cyg X-1. The detection establishes V404 Cygni as one of four microquasars with detectable gamma-ray emission, and adds interesting features to the small group, or in a more general context to X-ray binaries with jets.

We examine morphology-separated color-mass diagrams to study the quenching of star formation in $\sim 100,000$ ($z\sim0$) Sloan Digital Sky Survey (SDSS) and $\sim 20,000$ ($z\sim1$) Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS) galaxies. To classify galaxies morphologically, we developed Galaxy Morphology Network (GaMorNet), a convolutional neural network that classifies galaxies according to their bulge-to-total light ratio. GaMorNet does not need a large training set of real data and can be applied to data sets with a range of signal-to-noise ratios and spatial resolutions. GaMorNet's source code as well as the trained models are made public as part of this work ( http://www.astro.yale.edu/aghosh/gamornet.html ). We first trained GaMorNet on simulations of galaxies with a bulge and a disk component and then transfer learned using $\sim25\%$ of each data set to achieve misclassification rates of $\lesssim5\%$. The misclassified sample of galaxies is dominated by small galaxies with low signal-to-noise ratios. Using the GaMorNet classifications, we find that bulge- and disk-dominated galaxies have distinct color-mass diagrams, in agreement with previous studies. For both SDSS and CANDELS galaxies, disk-dominated galaxies peak in the blue cloud, across a broad range of masses, consistent with the slow exhaustion of star-forming gas with no rapid quenching. A small population of red disks is found at high mass ($\sim14\%$ of disks at $z\sim0$ and $2\%$ of disks at $z \sim 1$). In contrast, bulge-dominated galaxies are mostly red, with much smaller numbers down toward the blue cloud, suggesting rapid quenching and fast evolution across the green valley. This inferred difference in quenching mechanism is in agreement with previous studies that used other morphology classification techniques on much smaller samples at $z\sim0$ and $z\sim1$.

Axion is a popular candidate for dark matter particles. Axionic dark matter may form Bose-Einstein condensate and may be gravitationally bound to form axion clumps. Under the presence of electromagnetic waves with frequency $\omega=m_{a}/2$, where $m_{a}$ is the axion mass, a resonant enhancement may occur, causing an instability of the axion clumps. With analytical and numerical approaches, we study the resonant instability of axionic dark matter clumps with infinite homogeneous mass distribution, as well as distribution with a finite boundary. After taking realistic astrophysical environments into consideration, including gravitational redshift and plasma effects, we obtain an instability region in the axion density-clump size parameter space with given mass and coupling of axions. In particular, we show that, for axion clumps formed by the QCD axions in equilibrium, no resonant instability will occur.

We present the Atacama Large Millimeter/submillimeter Array (ALMA) observation of the Sub-millimeter galaxy (SMG) ID 141 at z=4.24 in the [N II] 205 $\mu$m line (hereafter [N II]) and the underlying continuum at (rest-frame) 197.6 $\mu$m. Benefiting from lensing magnification by a galaxy pair at z=0.595, ID 141 is one of the brightest z$>4$ SMGs. At the angular resolutions of $\sim1.2''$ to $1.5''$ ($1'' \sim6.9$ kpc), our observation clearly separates, and moderately resolves the two lensed images in both continuum and line emission at $\rm S/N>5$ . Our continuum-based lensing model implies an averaged amplification factor of $\sim5.8$ and reveals that the de-lensed continuum image has the S\'ersic index $\simeq 0.95$ and the S\'ersic radius of $\sim0.18'' (\sim 1.24$ kpc). Furthermore, the reconstructed [N II] velocity field in the source plane is dominated by a rotation component with a maximum velocity of $\sim 300$ km/s at large radii, indicating a dark matter halo mass of $\sim 10^{12}M_{\odot}$. This, together with the reconstructed velocity dispersion field being smooth and modest in value ($<100$ km/s) over much of the outer parts of the galaxy, favours the interpretation of ID 141 being a disk galaxy dynamically supported by rotation. The observed [N II]/CO (7-6) and [N II]/[C II] 158 $\mu$m line luminosity ratios, which are consistent with the corresponding line ratio vs. far-infrared color correlation from local luminous infrared galaxies, imply a de-lensed star formation rate of ($1.8\pm 0.6)\times10^3M_\odot$/yr and provide an independent estimate on the size of the star-forming region $0.7^{+0.3}_{-0.3}$ kpc in radius.

Fast radio bursts (FRBs) are mysterious millisecond-duration radio transients observed from extragalactic distances. Magnetars have been long speculated as the possible engine to power repeating bursts from FRB sources, but no convincing evidence has been collected so far. Recently, a Galactic magnetar dubbed Soft Gamma-ray Repeater (SGR) J1935+2154 entered an active phase by emitting intense soft gamma-ray bursts. One fast radio burst with two peaks (hereafter FRB 200428) and a luminosity comparable to those of extragalactic FRBs was detected from the source in association with a soft $\gamma$-ray / hard X-ray flare. Here we report an eight-hour targeted observational campaign with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) in four observing sessions, assisted by multi-wavelength (optical and hard X-rays) observations. During the third session, twenty-nine SGR bursts were detected by the {\it Fermi} Gamma-ray Burst Monitor (GBM). Throughout the observing period, especially during the scrutinized epochs when SGR bursts reached the Earth, no single dispersed pulsed emission was detected by FAST. This places a fluence upper limit eight orders of magnitude deeper than the fluence of FRB 200428. Our results suggest that FRB -- SGR burst associations are rather rare. FRBs may be highly relativistic and geometrically beamed or may carry narrow spectra with the characteristic frequencies of the most bursts far outside the FAST band. More intriguingly, it is also possible that the physical condition to achieve coherent radiation in SGR bursts is contrived. Only under extreme conditions could an FRB be made in association with an SGR burst.

In this work, we collect quasi-simultaneous infrared, optical, X-ray and $\gamma$-ray data of 60 $Fermi$-4LAC flat spectrum radio quasars (FSRQs). In the framework of the conventional one-zone leptonic model, we investigate the physical properties of $Fermi$-4LAC FSRQs' jets by modeling their quasi-simultaneous spectral energy distributions (SEDs). Our main results are summarized as follows. (1) There is a linear correlation between synchrotron peak frequency and curvature of the electron energy distribution. As suggested by previous works, the slope of the best linear fitting equation of this correlation is consistent with statistic acceleration which needs a fluctuation of fractional acceleration gain. (2) The gamma-ray dissipation regions are located at the range from 0.1 to 10 pc away from the super-massive black hole, and located outside the broad-line region (BLR) and within the dusty torus (DT). (3) A size relation $P_{\rm e}$ (the kinetic power carried in relativistic electrons) $\sim$ $P_{\rm B}$ (Poynting flux) $\leq$ $P_{\rm r}$ (the radiative power ) $<$ $P_{\rm p}$ (the kinetic power in cold protons) is found in our modeling. Among them, $P_{\rm e}\sim P_{\rm B}$ suggests that SEDs of almost all FSRQs with parameters are close to equipartition between the magnetic field and the relativistic electrons. The $P_{\rm e} < P_{\rm r}$ suggest that the most energy of the relativistic electrons are dissipated by EC radiation for FSRQs. (4) There is an anti-correlation between the peak energy of SEDs ($\gamma_{\rm peak}$) and the jet power ($P_{\rm jet}$), which is consistent with the blazar sequence.

We discuss the possibility of producing a light dark photon dark matter through a coupling between the dark photon field and the inflaton. The dark photon with a large wavelength is efficiently produced due to the inflaton motion during inflation and becomes non-relativistic before the time of matter-radiation equality. We compute the amount of production analytically. The correct relic abundance is realized with a dark photon mass extending down to $10^{-21} \, \rm eV$.

The 21~cm intensity mapping experiments promise to measure the baryon acoustic oscillations (BAO) through the mapping of H{\small I} gas at the post-reionization epoch. Its success depends on the assumption that H{\small I} gas traces the matter density distribution {\it linearly} on large scales. However, both nonlinear halo clustering and nonlinear effects modulating H{\small I} gas in halos may spoil this. We employ three approaches to generate the mock H{\small I} density from a large-scale N-body simulation at low redshifts, and demonstrate that the assumption of H{\small I} linearity is valid at the scale corresponding to the first BAO peak, but breaks down at $k \gtrsim 0.1\,h\, {\rm Mpc}^{-1}$. The nonlinear effects of halo clustering and H{\small I} content modulation counteract each other at small scales, and their competition results in a model-dependent "sweet-spot" redshift near $z$=1 where the H{\small I} bias is scale-independent down to small scales. We also find that the linear H{\small I} bias scales approximately linearly with redshift for $z\le 3$.

In a recent paper we proposed a mechanism for a continuous transition between a contracting Ekpyrotic phase and the Standard Big Bang phase of expansion: the bounce is generated by an S-brane which represents the effects of higher mass string states in the low energy effective field theory. We showed that gravitational waves on cosmological scales obtain a nearly scale-invariant spectrum. Here, we study the cosmological fluctuations in this setup, assuming that the S-brane has zero shear. We find a nearly scale-invariant spectrum of cosmological perturbations with a slight red tilt. The scenario yields two consistency relations for cosmological observations, the first one relating the tensor to scalar ratio with the scalar spectral tilt, the second relating the tensor tilt to the scalar tilt. The predicted tensor to scalar ratio is within the reach of upcoming CMB observations. The tensor tilt is blue.

We propose a mechanism borrowed from string theory which yields a non-singular transition from a phase of Ekpyrotic contraction to the expanding phase of Standard Big Bang cosmology. The same mechanism converts the initial vacuum spectrum of cosmological fluctuations before the bounce into a scale-invariant one, and also changes the spectrum of gravitational waves into an almost scale-invariant one. The scalar and tensor tilts are predicted to be the same, in contrast to the predictions from the "String Gas Cosmology" scenario. The amplitude of the gravitational wave spectrum depends on the ratio of the string scale to the Planck scale and may be in reach of upcoming experiments.

We show that it is possible to use gravitational wave detectors to observe the occurrence of a first order phase transition in Pati-Salam extensions of the Standard Model. We show that the peak frequency of the expected gravitational wave signals ranges within $0.1-10$ Hz. We find amusing that the next generation of gravity waves detectors are able to explore time honoured extensions of the Standard Model occurring at energy scales inaccessible by present and future particle physics accelerators.

We present a cosmic ray muon tomographic experiment for measuring the muon flux and imaging the tunnel overburden structures in Changshu, China. The device used in this study is a tracking detector based on the plastic scintillator with SiPM technology, which can be conveniently operated in field works. The compact system with sensitive area of $6400 cm^2$ can measure the angular distribution of cosmic muons. It's able to image the overburden density length from the surface of overburden to the detector along the muon tracks. The open sky muon flux measurement outside the tunnel has a good agreement with the modified Gassier Formula model. The distributions of muon flux at three positions inside the tunnel are very similar to that of open sky. Assuming the average density of overburden compact sandstone is $2.65 g/cm^3$, the overburden thickness can be obtained from the density length derived from the difference of muon flux inside and outside the tunnel. Moreover, for known penetrated lengths (i.e., topography of overburden), the density anomalies of the overburden can also been obtained. This study suggests a potential application for imaging and detecting subsurface structures in civil engineering, tunnels or caverns with the cosmic ray muon telescope.

Violent multi-wavelength variabilities are observed in gamma-ray-selected blazars. We present an analysis of long-term light curves for eight bright blazars to explore the co-variation pattern in the gamma-ray and radio bands. We extract their gamma-ray light curves and spectra with data observed by the Fermi/LAT since 2008. We find diverse co-variation patterns between the gamma-ray and radio (at 43 GHz) fluxes in these sources. The gamma-ray and radio fluxes of 3C 454.3 and PKS 1633+382 are correlated without any time-lag, suggesting that they are from the same radiation region. Similar correlation is also observed in 3C 273 and PKS 1222+216, but the radio flux is lag behind the gamma-ray flux approximately ~160 days and ~290 days, respectively. This likely suggests that their gamma-ray emission regions are located at the upstream of their radio cores at 43 GHz. The gamma-ray and radio fluxes of the other four blazars are not correlated, implying that the gamma-ray and radio emission may be from different regions in their jets. The gamma-ray light curves of the eight blazars can be decomposed into some long timescale variability components and fast spike flares. We propose that they may be attributed to the central engine activity and the magnetic reconnection process or turbulence in the local emission region, respectively.

Very high energy photons from cosmological gamma-ray bursts (GRBs) are expected to interact with extragalactic background light (EBL) and produce electron-positron pairs when they propagate through intergalactic medium (IGM). These relativistic pairs will then up-scatter cosmic microwave background (CMB) photons and emit secondary GeV emission. Meanwhile, the motion of these pairs are deflected by intergalactic magnetic field (IGMF), so the secondary GeV photons arrive later than the primary emission. It has been suggested that the properties of the secondary GeV emission can be used to constrain IGMF. Recently, TeV gamma-ray emission has been detected, for the first time, from a GRB (GRB 190114C) by the MAGIC telescope and its steep ${\rm \gamma-ray}$ spectrum shows a clear evidence of absorption by EBL. We then constrain the IGMF with the GeV flux limit obtained from the $Fermi$-LAT observations. We find a limit of $>10^{-19.5}$ G for the coherence length of $\lambda \leq 1$ Mpc. Although this limit is weaker than that obtained by using blazars, it represents the first limit from ${\rm \gamma-ray}$ observations of GRBs, which provides an independent constraint on IGMF. We also find that, for transient ${\rm \gamma-ray}$ sources, one can choose a favorable time window to search for the echo emission at a particular energy.

4FGL J0935.3+0901 is a gamma-ray source detected by the Large Area Telescope (LAT) onboard the Fermi Gamma-Ray Space Telescope. We have conducted detailed analysis of the LAT data for this source and multi-wavelength studies of the source field. Its gamma-ray emission can be described with a power law ($\Gamma=2.0\pm0.2$) with an exponential cutoff ($E_c=2.9\pm1.6$ GeV), while the flux shows significant long-term variations. From analysis of archival Neil Gehrels Swift X-ray Telescope data, we find only one X-ray source in the LAT's 2$\sigma$ error region. Within a 3.7" radius error circle of the X-ray source, there is only one optical object down to $r'\sim$23 mag. Time-resolved photometry of the optical object indicates a likely 2.5~hr periodic modulation, while its spectrum shows double-peaked hydrogen and helium emission lines (similar to those seen in accretion discs in low-mass X-ray binaries). Combining these results, we conclude that we have discovered a compact X-ray emitting binary in likely association with 4FGL J0935.3+0901, i.e., a millisecond pulsar (MSP) binary. We discuss the implication of the optical spectral features: this binary could be a transitional MSP system at a sub-luminous disc state, although the other possibility, the binary in a rotation-powered state showing the optical emission lines due to intrabinary interaction processes, can not be excluded. Further observational studies will help determine detailed properties of this candidate MSP binary and thus clarify its current state.

Similar to the Schwarzschild coordinates for spherical black holes, the Baldwin, Jeffery and Rosen (BJR) coordinates for plane gravitational waves are often singular, and extensions beyond such singularities are necessary, before studying asymptotic properties of such spacetimes at the null infinity of the plane, on which the gravitational waves propagate. The latter is closely related to the studies of memory effects and soft graviton theorems. In this paper, we point out that in the BJR coordinates all the spacetimes are singular physically at the focused point $u = u_s$, except for the two cases: (1) $\alpha =1/2, \; \forall \; \chi_n$; and (2) $\alpha =1, \; \chi_i = 0\; (i = 1, 2, 3)$, where $\chi_n$ are the coefficients in the expansion $\chi \equiv \left[{\mbox{det}}\left(g_{ab}\right) \right]^{1/4} = \left(u - u_s\right)^{\alpha}\sum_{n = 0}^{\infty}\chi_n \left(u - u_s\right)^n$ with $\chi_0 \not= 0$, the constant $\alpha \in (0, 1]$ characterizes the strength of the singularities, and $g_{ab}$ denotes the reduced metric on the two-dimensional plane orthogonal to the propagation direction of the wave. Therefore, the hypersurfaces $u= u_s$ already represent the boundaries of such spacetimes, and the null infinity does not belong to them. As a result, they cannot be used to study properties of plane gravitational waves at null infinities, including memory effects and soft graviton theorems.

A comprehensively theoretical analysis on the broadband spectral energy distributions (SEDs) of large-scale jet knots in 3C 273 is presented for revealing their X-ray radiation mechanism. We show that these SEDs cannot be explained with a single electron population model when the Doppler boosting effect is either considered or not. By adding a more energetic electron (the leptonic model) or proton (the hadronic model) population, the SEDs of all knots are well represented. In the leptonic model, the electron population that contributes the X-ray emission is more energetic than the one responsible for the radio-optical emission by almost two orders of magnitude; the derived equipartition magnetic field strengths (B_eq) are ~0.1 mG. In the hadronic model, the protons with energy of ~20 PeV are required to interpret the observed X-rays; the B_eq values are several mG, larger than that in the leptonic model. Based on the fact that no resolved substructures are observed in these knots and the fast cooling-time of the high-energy electrons is difficult to explain the observed X-ray morphologies, we argue that two distinct electron populations accelerated in these knots are unreasonable and their X-ray emission would be attributed to the proton synchrotron radiation accelerated in these knots. In case of these knots have relativistic motion towards the observer, the super-Eddington issue of the hadronic model could be avoided. Multiwavelength polarimetry and the gamma-ray observations with high resolution may be helpful to discriminate these models.

In this paper we present our imaging observations on the CO(7-6) line and its underlying continuum emission of the young submillimeter galaxy LESS 073 at redshift 4.755, using the Atacama Large Millimeter/submillimeter Array (ALMA). At the achieved resolution of $\sim$$1^{\prime\prime}.2\times0^{\prime\prime}.9$ ($8\times6$~kpc$^2$), the CO(7-6) emission is largely unresolved (with a deconvolved size of $1^{\prime\prime}.1(\pm0^{\prime\prime}.5) \times 0^{\prime\prime}.9(\pm0^{\prime\prime}.8)$), and the continuum emission is totally unresolved. The CO(7-6) line emission has an integrated flux of $0.86\pm0.08$~Jy km/s, and a line width of $343\pm40$ km/s. The continuum emission has a flux density of 0.51 mJy. By fitting the observed far-infrared (FIR) spectral energy distribution of LESS 073 with a single-temperature modified blackbody function, we obtained a dust temperature $T_{\rm dust}=57.6\pm3.5$ K, 60-to-100 $\mu$m flux density ratio $f_{60}/f_{100}=0.86\pm0.08$, and total infrared luminosity $L_{\rm IR}=(5.8\pm0.9) \times 10^{12}~L_\odot$. The SED-fit-based $f_{60}/f_{100}$ is consistent with those estimated from various line ratios as advocated by our earlier work, indicating that those proposed line-ratio-based method can be used to practically derive $f_{60}/f_{100}$ for high-$z$ sources. The total molecular gas mass of LESS 073 is $(3.3\pm1.7) \times10^{10}~M_\odot$, and the inferred gas depletion time is about 43 Myr.

The Large Optical/infrared Telescope (LOT) is a ground-based 12m diameter optical/infrared telescope which is proposed to be built in the western part of China in the next decade. Based on satellite remote sensing data, along with geographical, logistical and political considerations, three candidate sites were chosen for ground-based astronomical performance monitoring. These sites include: Ali in Tibet, Daocheng in Sichuan, and Muztagh Ata in Xinjiang. Up until now, all three sites have continuously collected data for two years. In this paper, we will introduce this site testing campaign, and present its monitoring results obtained during the period between March 2017 and March 2019.

ZTF J153932.16+502738.8 (ZTFJ1539) is an eclipsing double-white-dwarf system with an orbital period of 6.91 minutes, and is a significant source of LISA detecting gravitational wave. However, the massive white dwarf (WD) with mass of about 0.61 M$_\odot$ has a high effective temperature (48900 K), and the lower mass WD with mass of about 0.21 M$_{\odot}$ has a low effective temperature($<$10000 K). It is challenging the popular theory of binary evolution. We investigate the formation of ZTFJ1539 via nova and Algol scenarios. Assuming that the massive WD in ZTFJ1539 just experiences a thermalnuclear runaway, nova scenario can explain the effective temperatures of two WDs in ZTFJ1539. However, in order to enlarging a semi-detached orbit of about 4---5 minutes to a detached orbit of about 7 minutes, nova scenario needs a much high kick velocity of about 200 km s$^{-1}$ during nova eruption. The high kick velocity can result in high eccentricity of about 0.2---0.6. Algol scenario can also produce ZTFJ1539 if we take a high efficient parameter for ejecting common envelope and enhance the mass-loss rate via stellar wind trigger by tidal effect.

Although habitability, defined as the general possibility of hosting life, is expected to occur under a broad range of conditions, the standard scenario to allow for habitable environments is often described through habitable zones (HZs). Previous work indicates that stellar binary systems typically possess S-type or P-type HZs, with the S-type HZs forming ring-type structures around the individual stars and P-type HZs forming similar structures around both stars, if considered a pair. However, depending on the stellar and orbital parameters of the system, typically, there are also regions within the systems outside of the HZs, referred to as dead zones (DZs). In this study, we will convey quantitative information on the width and location of DZs for various systems. The results will also depend on the definition of the stellar HZs as those are informed by the planetary climate models.

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 results from a search for a radio transient associated with the LIGO/Virgo source S190814bv, a likely neutron star-black hole (NSBH) merger, with the Australian Square Kilometre Array Pathfinder. We imaged a $30\,{\rm deg}^2$ field at $\Delta T$=2, 9 and 33 days post-merger at a frequency of 944\,MHz, comparing them to reference images from the Rapid ASKAP Continuum Survey observed 110 days prior to the event. Each epoch of our observations covers $89\%$ of the LIGO/Virgo localisation region. We conducted an untargeted search for radio transients in this field, resulting in 21 candidates. For one of these, \object[AT2019osy]{AT2019osy}, we performed multi-wavelength follow-up and ultimately ruled out the association with S190814bv. All other candidates are likely unrelated variables, but we cannot conclusively rule them out. We discuss our results in the context of model predictions for radio emission from neutron star-black hole mergers and place constrains on the circum-merger density and inclination angle of the merger. This survey is simultaneously the first large-scale radio follow-up of an NSBH merger, and the most sensitive widefield radio transients search to-date.

The formation of bimodal dust species (namely the silicate and amorphous carbon dust grains coexistent) in a nova eruption is an open problem. According to the nova model simulated by Modules for Experiments in Stellar Astrophysics code, we calculate the formation and growth of carbon (C) and forsterite (Mg2SiO4) dust grains purely in nova ejecta for the fee-expansion model and the radiative shock model by assuming spherical geometry of the nova ejecta. In these models, the chemical properties of pre-existing circumstellar medium are not taken into account. In the free-expansion model, the nova ejecta is not an idea environment for dust nucleation. However, it can efficiently produce dust in the radiative shock model. We estimate that every nova can produce C grains with an average mass of about $10^{-9}$ and $10^{-8}$ ${\rm M_\odot}$, and Mg$_2$SiO$_4$ grains with an average mass of about $10^{-8}$ and $10^{-7}$ ${\rm M_\odot}$. Based on the mass of ejected gas, the ratio of dust to gas is about 1\%. The C grains form first after several or tens of days of nova eruption. After that, the Mg$_2$SiO$_4$ grains begin to grow in tens of days, which is consistent with observations.

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

Subdwarf A-type stars (sdAs) are objects that have hydrogen-rich spectra with surface gravity similar to that of hot subdwarf stars but effective temperature below the zero-age horizontal branch (ZAHB). They are considered to be metal-poor main sequence (MS) stars or extremely low-mass white dwarfs (ELM WDs). In this work, using the stellar evolution code, Modules for Experiments in Stellar Astrophysics (MESA), we investigate the sdAs formed both by the evolution of (pre-)ELM WDs in double degenerate systems (DDs) and metal-poor main sequence with the single evolution models. We find that both of the evolutionary tracks of ELM WDs and metal-poor MSs can explain the observation properties of sdAs. However, the proportion between these two populations are uncertain. In this work, we adopt the method of binary population synthesis of both ELM WDs in the disk and metal-poor MS in the halo to obtain the populations of ELM WDs and metal-poor MSs at different stellar population ages and calculate their proportions. We find that the proportion of metal-poor MSs to sdA population for a stellar population of 10 Gyr is $\sim$98.5 percent, which is consistent with the conclusion that most sdAs ($>$95 percent) are metal-poor MSs. And the proportion of ELM WDs (metal-poor MSs) to sdA population increases (decreases) from 0.1\% (99.9\%) to 20\% (80\%) with the stellar population ages from 5 to 13.7 Gyr.

In our search for gamma-ray emission from Be X-ray binaries from analysis of the data obtained with the Large Area Telescope (LAT) on board the Fermi Gamma-Ray Space Telescope, we find likely detection of GRO J1008-57. The binary has an orbital period of 249.48 days, and it is only significantly detected in its orbital phase 0.8--0.9 (>4 sigma). Further light curve analysis indicates that the detection is probably largely due to an emitting event in one orbital cycle around year 2012--2013, following a giant X-ray outburst of the source. This property of having occasional gamma-ray emitting events is similar to that seen in another high-mass X-ray binary 4U 1036-56. However, models considering possible gamma-ray emission from an accreting neutron star have difficulty in explaining the observed 10^{34} erg s^{-1} luminosity of the source, unless the distance was largely over-estimated. Further observational studies are required, in order to more clearly establish the high-energy emission properties of GRO J1008-57 or similar high-mass X-ray binaries and find clues for understanding how gamma-ray emission is produced from them.

A high-energy muon neutrino event, IceCube-170922A, was recently discovered in both spatial and temporal coincidence with a gamma-ray flare of the blazar TXS 0506+056. It has been shown, with standard one-zone models, that neutrinos can be produced in the blazar jet via hadronic interactions, but with a flux which is mostly limited by the X-ray data. In this work, we explore the neutrino production from TXS 0506+056 by invoking two physically distinct emission zones in the jet, separated by the broad line region (BLR). Using the Doppler-boosted radiation of the BLR as the target photon field, the inner zone accounts for the neutrino and gamma-ray emission via $p\gamma$ interactions and inverse Compton scattering respectively, while the outer zone produces the optical and X-ray emission via synchrotron and synchrotron self-Compton processes. The different conditions of the two zones allow us to suppress the X-ray emission from the electromagnetic cascade, and set a much higher upper limit on the muon neutrino flux (i.e., $\sim 10^{-11}\rm erg~cm^{-2}s^{-1}$) than in one-zone models. We compare, in detail, our scenario with one-zone models discussed in the literature, and argue that differentiating between such scenarios will become possible with next generation neutrino telescopes, such as IceCube-Gen2.

Three- and five-minute oscillations of sunspots have different spatial distributions in the solar atmospheric layers. The spatial distributions are crucial to reveal the physical origin of sunspot oscillations and to investigate their propagation. In this study, six sunspots observed by Solar Dynamics Observatory/Atmospheric Imaging Assembly were used to obtain the spatial distributions of three- and five-minute oscillations. The fast Fourier transform method is applied to represent the power spectra of oscillation modes. We find that, from the temperature minimum to the lower corona, the powers of the five-minute oscillation exhibit a circle-shape distribution around its umbra, and the shapes gradually expand with temperature increase. However, the circle-shape is disappeared and the powers of the oscillations appear to be very disordered in the higher corona. This indicates that the five-minute oscillation can be suppressed in the high-temperature region. For the three-minute oscillations, from the temperature minimum to the high corona, their powers mostly distribute within an umbra, and part of them locate at the coronal fan loop structures. Moreover, those relative higher powers are mostly concentrated in the position of coronal loop footpoints.

It has been suggested that low energy effective field theories should satisfy given conditions in order to be successfully embedded into string theory. In the case of a single canonically normalized scalar field this translates into conditions on its potential and the derivatives thereof. In this Letter we revisit stochastic models of small field inflation and study the compatibility of the swampland constraints with entropy considerations. We show that stochastic inflation either violates entropy bounds or the swampland criterium on the slope of the scalar field potential. Furthermore, we illustrate that such models are faced with a graceful exit problem: any patch of space which exits the region of eternal inflation is either not large enough to explain the isotropy of the cosmic microwave background, or has a spectrum of fluctuations with an unacceptably large red tilt.

The Gamma-Ray Integrated Detectors (GRID) is a space mission concept dedicated to monitoring the transient gamma-ray sky in the energy range from 10 keV to 2 MeV using scintillation detectors onboard CubeSats in low Earth orbits. The primary targets of GRID are the gamma-ray bursts (GRBs) in the local universe. The scientific goal of GRID is, in synergy with ground-based gravitational wave (GW) detectors such as LIGO and VIRGO, to accumulate a sample of GRBs associated with the merger of two compact stars and study jets and related physics of those objects. It also involves observing and studying other gamma-ray transients such as long GRBs, soft gamma-ray repeaters, terrestrial gamma-ray flashes, and solar flares. With multiple CubeSats in various orbits, GRID is unaffected by the Earth occultation and serves as a full-time and all-sky monitor. Assuming a horizon of 200 Mpc for ground-based GW detectors, we expect to see a few associated GW-GRB events per year. With about 10 CubeSats in operation, GRID is capable of localizing a faint GRB like 170817A with a 90% error radius of about 10 degrees, through triangulation and flux modulation. GRID is proposed and developed by students, with considerable contribution from undergraduate students, and will remain operated as a student project in the future. The current GRID collaboration involves more than 20 institutes and keeps growing. On August 29th, the first GRID detector onboard a CubeSat was launched into a Sun-synchronous orbit and is currently under test.

The Gamma-Ray Integrated Detectors (GRID) is a space mission concept dedicated to monitoring the transient gamma-ray sky in the energy range from 10 keV to 2 MeV using scintillation detectors onboard CubeSats in low Earth orbits. The primary targets of GRID are the gamma-ray bursts (GRBs) in the local universe. The scientific goal of GRID is, in synergy with ground-based gravitational wave (GW) detectors such as LIGO and VIRGO, to accumulate a sample of GRBs associated with the merger of two compact stars and study jets and related physics of those objects. It also involves observing and studying other gamma-ray transients such as long GRBs, soft gamma-ray repeaters, terrestrial gamma-ray flashes, and solar flares. With multiple CubeSats in various orbits, GRID is unaffected by the Earth occultation and serves as a full-time and all-sky monitor. Assuming a horizon of 200 Mpc for ground-based GW detectors, we expect to see a few associated GW-GRB events per year. With about 10 CubeSats in operation, GRID is capable of localizing a faint GRB like 170817A with a 90% error radius of about 10 degrees, through triangulation and flux modulation. GRID is proposed and developed by students, with considerable contribution from undergraduate students, and will remain operated as a student project in the future. The current GRID collaboration involves more than 20 institutes and keeps growing. On August 29th, the first GRID detector onboard a CubeSat was launched into a Sun-synchronous orbit and is currently under test.

The DArk Matter Particle Explorer (DAMPE), a satellite-based cosmic ray and gamma-ray detector, was launched on December 17, 2015, and began its on-orbit operation on December 24, 2015. In this work we document the on-orbit calibration procedures used by DAMPE and report the calibration results of the Plastic Scintillator strip Detector (PSD), the Silicon-Tungsten tracKer-converter (STK), the BGO imaging calorimeter (BGO), and the Neutron Detector (NUD). The results are obtained using Galactic cosmic rays, bright known GeV gamma-ray sources, and charge injection into the front-end electronics of each sub-detector. The determination of the boundary of the South Atlantic Anomaly (SAA), the measurement of the live time, and the alignments of the detectors are also introduced. The calibration results demonstrate the stability of the detectors in almost two years of the on-orbit operation.

Isolated millisecond pulsars (IMSPs) are a topic of academic contention. There are various models to explain their formation. We explore the formation of IMSP via quark novae (QN). During this formation process, low-mass X-ray binaries (LMXBs) are disrupted when the mass of the neutron star (NS) reaches 1.8 $\rm M_\odot$. Using population synthesis, this work estimates that the Galactic birthrate of QN-produced IMSPs lies between $\sim 9.5\times10^{-6}$ and $\sim1.7\times10^{-4}$ $\rm yr^{-1}$. The uncertainties shown in our experiment model is due to the QN's kick velocity. Furthermore, our findings not only show that QN-produced IMSPs are statistically more significant than those produced by mergers, but also that millisecond pulsar binaries with a high eccentricity may originate from LMXBs that have been involved in, yet not disrupted by, a QN.

We report a ground X-ray calibration of two X-ray telescope prototypes at the PANTER X-ray Test Facility, of the Max-Planck-Institute for Extraterrestrial Physics, in Neuried, Germany. The X-ray telescope prototypes were developed by the Institute of Precision Optical Engineering (IPOE) of Tongji University, in a conical Wolter-I configuration, using thermal glass slumping technology. Prototype #1 with 3 layers and Prototype #2 with 21 layers were tested to assess the prototypes' on-axis imaging performance. The measurement of Prototype #1 indicates a Half Power Diameter (HPD) of 82" at 1.49 keV. As for Prototype #2, we performed more comprehensive measurements of on-axis angular resolution and effective area at several energies ranging from 0.5-10 keV. The HPD and effective area are 111" and 39 cm^2 at 1.49 keV, respectively, at which energy the on-axis performance of the prototypes is our greatest concern.

The formation of bimodal dust species (that is, both of silicate and amorphous carbon dust grains are observed in a nova eruption) in nova ejecta is still debated. Using the Modules for Experiments in Stellar Astrophysics code and considering the effects of WD's mass, mass-accretion rate and the chemical profiles of WD which are described by new parameter --- mixing depth on the chemical abundances of nova ejecta, we investigate the possibility that bimodal dust species are produced in a nova eruption. We find that $C/O$ (the ratio of carbon number density to oxygen number density) of nova ejecta is affected by the mixing depth. For the model with a small mixing depth, the $C/O$ of nova ejecta can evolute from lager than 1.0 to smaller than 1.0 in a whole eruption, which provides the chemical condition for the formation of bimodal dust species.

We develop a novel method to explore the galaxy-halo connection using the galaxy imaging surveys by modeling the projected two-point correlation function measured from the galaxies with reasonable photometric redshift measurements. By assuming a Gaussian form of the photometric redshift errors, we are able to simultaneously constrain the halo occupation distribution (HOD) models and the effective photometric redshift uncertainties. Tests with mock galaxy catalogs demonstrate that this method can successfully recover (within $\sim 1\sigma$) the intrinsic large-scale galaxy bias, as well as the HOD models and the effective photometric redshift uncertainty. This method also works well even for galaxy samples with 10 per cent catastrophic photometric redshift errors.

The Large High Altitude Air Shower Observatory (LHAASO) project is a new generation multi-component instrument, to be built at 4410 meters of altitude in the Sichuan province of China, with the aim to study with unprecedented sensitivity the spec trum, the composition and the anisotropy of cosmic rays in the energy range between 10$^{12}$ and 10$^{18}$ eV, as well as to act simultaneously as a wide aperture (one stereoradiant), continuously-operated gamma ray telescope in the energy range between 10$^{11}$ and $10^{15}$ eV. The experiment will be able of continuously surveying the TeV sky for steady and transient sources from 100 GeV to 1 PeV, t hus opening for the first time the 100-1000 TeV range to the direct observations of the high energy cosmic ray sources. In addition, the different observables (electronic, muonic and Cherenkov/fluorescence components) that will be measured in LHAASO will allow to investigate origin, acceleration and propagation of the radiation through a measurement of energy spec trum, elemental composition and anisotropy with unprecedented resolution. The remarkable sensitivity of LHAASO in cosmic rays physics and gamma astronomy would play a key-role in the comprehensive general program to explore the High Energy Universe. LHAASO will allow important studies of fundamental physics (such as indirect dark matter search, Lorentz invariance violation, quantum gravity) and solar and heliospheric physics. In this document we introduce the concept of LHAASO and the main science goals, providing an overview of the project.

We report the first experimental results on spin-dependent elastic weakly interacting massive particle (WIMP) nucleon scattering from the XENON1T dark matter search experiment. The analysis uses the full ton year exposure of XENON1T to constrain the spin-dependent proton-only and neutron-only cases. No significant signal excess is observed, and a profile likelihood ratio analysis is used to set exclusion limits on the WIMP-nucleon interactions. This includes the most stringent constraint to date on the WIMP-neutron cross section, with a minimum of $6.3\times10^{-42}$ cm$^2$ at 30 GeV/c${}^2$ and 90% confidence level. The results are compared with those from collider searches and used to exclude new parameter space in an isoscalar theory with an axial-vector mediator.

The data from the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope have recently been updated. We thus re-analyze the LAT data for the supernova remnant (SNR) SN 1006. Two parts of gamma-ray emission from the region is clearly resolved, which correspond to the north-east (NE) and south-west (SW) limbs of the SNR. The former was detected in the previous LAT data (Xing et al. 2016), but the latter is newly detected in this work. The detection of the two limbs are at a 4 sigma significance level, and the spectral results for the NE limb is consistent with those obtained in previous detection analyses. We construct the broadband spectral energy distribution (SED) for the SW limb. Different scenarios are considered for the SED in gamma-ray energies. We conclude very similar to that of the NE limb, the high-energy and very high-energy emission from the SW limb is likely dominated by the leptonic process, in which high-energy electrons accelerated from the shell region of the SNR inverse-Compton scatter background photons to gamma-rays.

Gamma-ray burst (GRB) central engines and jet production mechanisms are still open questions. Assuming that the shallow decay segments of canonical X-ray afterglow lightcurves of {\em Swift} GRBs are attributed to the magnetic dipole (MD) radiations of newly-born magnetars, we derive the parameters of the magnetars and explore their possible relations to jet and MD wind emission. We show that the magnetar initial spin period ($P_0$) are tightly correlated with the jet energy ($E_{\rm jet}$), which is almost proportional to the wind energy ($E_{\rm wind}$). Our least square fits yield $P_0\propto E^{-0.36\pm 0.03}_{\rm jet}$ and $E_{\rm wind}\propto E^{0.91\pm 0.07}_{\rm jet}$. These relations may imply that a magnetar with faster rotating speed can power a more energetic GRB, and energy partition between the jet and wind may be quasi-universal. Although the $P_0-E_{\rm jet}$ relation is driven by a few sub-energetic GRBs in our sample, our Monte Carlo simulation analysis shows that sample selection biases from instrumental flux limits and contaminations of the bright jet afterglows cannot make this correlation. Within this jet-wind paradigm, we propose that GRB 101225A-like ultra-long GRBs, whose prompt gamma-ray/X-ray lightcurves are featured as a long-lasting plateau with a sharp drop, may be the orphan MD wind emission being due to misalignment of their jet axis to the light of sight. Brief discussion on the orphan MD wind emission and its association with the gravitational wave radiation of newly-born magnetars is presented.

The Near-infrared Sky Brightness Monitor (NIRBM) aims to measure the middle infrared sky background in Antarctica. The NIRBM mainly consists of an InGaAs detector, a chopper, a reflector, a cooler and a black body. The reflector can rotate to scan the sky with a field of view ranging from 0{\deg} to 180{\deg}. Electromechanical control and weak signal readout functions are accomplished by the same circuit, whose core chip is a STM32F407VG microcontroller. Considering the environment is harsh for humans in Antarctica, a multi-level remote control software system is designed and implemented. A set of EPICS IOCs are developed to control each hardware module independently via serial port communication with the STM32 microcontroller. The tornado web framework and PyEpics are introduced as a combination where PyEpics is used to monitor or change the EPICS Process Variables, functioning as a client for the EPICS framework. Tornado is responsible for the specific operation process of inter-device collaboration, and expose a set of interfaces to users to make calls. Considering the high delay and low bandwidth of the network environment, the tornado back-end is designed as a master-and-agent architecture to improve domestic user experience. The master node is deployed in Antarctic while multiple agent nodes can be deployed domestic. The master and agent nodes communicate with each other through the WebSocket protocol to exchange latest information so that bandwidth is saved. The GUI is implemented in the form of single-page application based on the Vue framework which communicates with tornado through WebSocket and AJAX requests. The web page integrates device control, data curve drawing, alarm display, auto observation and other functions together.

We are starting a project to find $\gamma$-ray millisecond pulsars (MSPs) among the unidentified sources detected by the Large Area Telescope (LAT) onboard the {Fermi Gamma-Ray Space Telescope (Fermi)}, by radio observations. The selection of good candidates from analysis of the LAT data is an important part of the project. Given that there is more than 10 years worth of LAT data and the advent of the newly released LAT 8-year point source list (FL8Y), we have conducted a selection analysis, on the basis of our previous analysis, and report the results here. Setting the requirements for the unidentified sources in FL8Y of Galactic latitudes $|b|>$ 5$^\circ$ and curvature significances $>$3$\sigma$, there are 202 sources with detection signficances $>$6$\sigma$. We select 57 relatively bright ones (detection significances $>$15$\sigma$) and analyze their 10.2 years of LAT data. Their variability is checked to exclude variable sources (likely~blazars), test statistic maps are constructed to avoid contaminated sources, and curvature significances are re-obtained and compared to their $\gamma$-ray spectra to exclude non-significant sources. In the end, 48 candidates are found. Based on the available information, mostly from multi-wavelength studies, we discuss the possible nature of several of the candidates. Most of these candidates are currently being observed with the 65-meter Shanghai Tian Ma Radio Telescope.

We report our analysis of the data obtained with the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope (Fermi) for the SS 433/W50 region. The total data length is ten years. Different from the previous results reported for this region (for which an old version of database was used), we show that excess emission is detected around the w1 region in the western lobe of the jets from SS 433. The region is bright at X-rays due to the interaction between the jet and the ambient medium. This detection also matches well the recent results of the very-high-energy detection of SS 433/W50 with the High Altitude Water Cherenkov (HAWC). However, the eastern regions that are slightly brighter in HAWC's observation are not detected in the Fermi data. Constructing the broad-band spectral energy distribution (SED) for the western region, we compare with the HAWC results for the eastern regions and discuss the possible origin of the emission. A leptonic scenario can provide a fit to the Fermi GeV spectrum and HAWC TeV detection, where the former and latter are due to the synchrotron radiation and inverse-Compton process respectively. However, the model can not explain the X-ray and radio emission from the region simultaneously, which thus requires further observational and theoretical studies of the region in order to clarify the reasons.

We present first results on the scalar WIMP-pion coupling from 1 t$\times$yr of exposure with the XENON1T experiment. This interaction is generated when the WIMP couples to a virtual pion exchanged between the nucleons in a nucleus. In contrast to most non-relativistic operators, these pion-exchange currents can be coherently enhanced by the total number of nucleons, and therefore may dominate in scenarios where spin-independent WIMP-nucleon interactions are suppressed. Moreover, for natural values of the couplings, they dominate over the spin-dependent channel due to their coherence in the nucleus. Using the signal model of this new WIMP-pion channel, no significant excess is found, leading to an upper limit cross section of $6.4\times10^{-46}$ cm$^2$ (90 % confidence level) at 30 GeV/c$^2$ WIMP mass.

In Paper I and II, a comprehensive approach was utilized for the calculation of S-type and P-type habitable regions in stellar binary systems for both circular and elliptical orbits of the binary components. It considered a joint constraint including orbital stability and a habitable region for a possible system planet through the stellar radiative energy fluxes ("radiative habitable zone"; RHZ). Specifically, the stellar S-type and P-type RHZs are calculated based on the solution of a fourth order polynomial. However, in concurrent developments, mostly during 2013 and 2014, important improvements have been made in the computation of stellar habitable zones for single stars based on updated climate models given by R. K. Kopparapu and collaborators. These models entail considerable changes for the inner and outer limits of the stellar habitable zones. Moreover, regarding the habitability limit given by the runaway greenhouse effect, notable disparities were identified between Earth, Mars, and super-Earth planets due to differences in their atmospheric models, thus affecting their potential for habitability. It is the aim of this study to compute S-type and P-type habitable regions of binaries in response to the updated planetary models. Moreover, our study will also consider improved relationships between effective temperatures, radii, and masses for low-luminosity stars.

Initially unpolarized solar radiation becomes polarized by scattering in the Earth's atmosphere. In particular molecular scattering polarizes electromagnetic radiation, but also scattering of radiation at aerosols, cloud droplets and ice crystals polarizes. Each atmospheric constituent produces a characteristic polarization signal, thus spectro-polarimetric measurements are frequently employed for remote sensing of aerosol and cloud properties. Retrieval algorithms require efficient radiative transfer models. Usually, these apply the plane-parallel approximation, assuming that the atmosphere consists of horizontally homogeneous layers. For remote sensing applications, the radiance is considered constant over the instantaneous field-of-view of the instrument and each sensor element is treated independently in plane-parallel approximation, neglecting horizontal radiation transport between adjacent pixels. In order to estimate the errors due to the IPA approximation, three-dimensional (3D) vector radiative transfer models are required. So far, only a few such models exist. Therefore, the International Polarized Radiative Transfer (IPRT) working group of the International Radiation Commission (IRC) has initiated a model intercomparison project in order to provide benchmark results for polarized radiative transfer. The group has already performed an intercomparison for one-dimensional (1D) multi-layer test cases (Emde et al., 2015). This paper presents the continuation of the intercomparison project for 2D and 3D test cases: a step cloud, a cubic cloud, and a more realistic scenario including a 3D cloud field generated by a Large Eddy Simulation model and typical background aerosols. The commonly established benchmark results for 3D polarized radiative transfer are available at the IPRT website (http://www.meteo.physik.uni-muenchen.de/~iprt).

This paper uses population synthesis to investigate the possible origin of isolated millisecond pulsars (IMSPs) as born from the coalescence of a neutron star (NS) and a white dwarf (WD). Results show that the galactic birth - rate of IMSPs is likely to lie between $5.8\times10^{-5}$ $\rm yr^{-1}$ and $2.0\times10^{-4}$ $\rm yr^{-1}$, depending on critical variables, such as the stability of mass transfer via the Roche lobe and the value of kick velocity. In addition to this, this paper estimates that the solar mass of IMSPs can range from $1.5$ and 2.0 $M_{\odot}$, making them more massive than other 'normal' pulsars. Finally, the majority of IMSPs in our simulations have spin periods ranging from several to 20 milliseconds, which is consistent with previous observations.

The third-generation gravitational wave detectors are under development by operating the detector in cryogenic temperature to reduce the thermal noise. Silicon and sapphire are promising candidate materials for the test masses and suspension elements due to their remarkable mechanical and thermal properties at cryogenic temperature. Here we present the performances of the cryogenic thermal cycling and strength testing on hydroxide catalysis bonding between sapphire and silicon. Our results suggest that although these two materials have very different coefficients of thermal expansion, but if the flatness and the thermally grown $\mathrm{SiO_2}$ oxidation layer on the silicon surface are controlled well, the bonded samples can still survive thermal cycling from room temperature to 5.5 K. A breaking strength of 3.6$\pm 0.6$ MPa is measured for the bonds between sapphire and silicon with a 190 nm silicon oxidation thickness after cooling cycle. We construct a hybrid sapphire-silicon Fabry-Perot cavity with the developing bonding technique in our lab. The measurement results reveal that the cavity can survive repeated thermal cycling while maintaining a good finesse.

In this paper we present the enhanced X-ray Timing and Polarimetry mission - eXTP. eXTP is a space science mission designed to study fundamental physics under extreme conditions of density, gravity and magnetism. The mission aims at determining the equation of state of matter at supra-nuclear density, measuring effects of QED, and understanding the dynamics of matter in strong-field gravity. In addition to investigating fundamental physics, eXTP will be a very powerful observatory for astrophysics that will provide observations of unprecedented quality on a variety of galactic and extragalactic objects. In particular, its wide field monitoring capabilities will be highly instrumental to detect the electro-magnetic counterparts of gravitational wave sources. The paper provides a detailed description of: (1) the technological and technical aspects, and the expected performance of the instruments of the scientific payload; (2) the elements and functions of the mission, from the spacecraft to the ground segment.

Space weathering is an important surface process occurring on the Moon and other airless bodies, especially those that have no magnetic field. The optical effects of the Moon's space weathering have been largely investigated in the laboratory for lunar samples and lunar analogues. However, duplication of the pristine regolith on Earth is not possible. Here we report the space weathering from the unique perspective of the Chang'E-3's (CE-3) "Yutu" rover, building on our previous work (Wang et al. 2017; Wu and Hapke 2018). Measurement of the visually undisturbed uppermost regolith as well as locations that have been affected by rocket exhaust from the spacecraft by the Visible-Near Infrared Spectrometer (VNIS) revealed that the returned samples bring a biased information about the pristine lunar regolith. The uppermost surficial regolith is much more weathered than the regolith immediately below, and the finest fraction is rich in space weathered products. These materials are very dark and attenuated throughout the visible and near-infrared (VNIR) wavelengths, hence reduce the reflectance and mask the absorption features. The effects on the spectral slope caused by space weathering are wavelength-dependent: the visible and near-infrared continuum slope (VNCS) increases while the visible slope (VS) decreases. In the visible wavelength, the optical effects of space weathering and TiO2 are identical: both reduce albedo and blue the spectra. This suggests that developing new TiO2 abundance algorithm is needed. Optical maturity indices are composition related and hence only locally meaningful. Since optical remote sensing can only sense the uppermost few microns of regolith and since this surface tends to be very weathered, the interpretation of surface composition using optical remote sensing data needs to be carefully evaluated. Sampling the uppermost surface is suggested.

Highly precise pulsar timing is very important for understanding the nature of a neutron star, and it can even be used to detect gravitational waves. Unfortunately, the accuracy of the pulsar timing is seriously affected by the spin-down irregularities of pulsars.

Highly precise pulsar timing is very important for understanding the nature of neutron stars and the implied physics, even being used to detect directly gravitational waves. Unfortunately, the accuracy of the pulsar timing is seriously affected by the spin-down irregularities, such as spin fluctuations with a manifestation of low frequency structures (the so-called red noise processes), and various activities of magnetospheres. Except from the random timing noise, significant periodic or quasi-periodic components could also be found in some timing residuals of pulsars, indicating the presence of unmodelled deterministic effects. The physical origins of these effects still remain unexplained. In this paper, we suggest a new mechanism involving the de Haas-van Alphen magnetic oscillation, which could trigger the observed low frequency structures. We find that the proposed magnetic oscillation period is about 1-$10^{2}$ yr, which is about $10^{-4}$ times as long as the classical characteristic time scale of interior magnetic field evolution for a normal neutron star. Due to the magnetic oscillation, we estimate both the braking index range between $10^{-5}$ and $10^{5}$ and the residuals range between 4 to 906 ms for some specific samples of quasi-periodic pulsars. Those are consistent with the pulsar timing observations. Nonlinear phenomena of de Haas-van Alphen effect and the Condon domain structure are considered. Avalanche-like magnetized process of domain wall motion (rapidly magnetic energy release) could associated with some special emitting manifestations of neuron star radiation.

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.

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

Using atomistic simulations, we characterize the adsorption process of organic molecules on carbon nanoparticles, both of which have been reported to be abundant in the interstellar medium (ISM). It is found that the aromatic organics are adsorbed more readily than the aliphatic ones. This selectivity would favor the formation of polycyclic aromatic hydrocarbons (PAHs) or fullerene-like structures in the ISM due to structural similarity. It is also observed in our simulations that the molecules form a monolayer over the nanoparticle surface before stacking up in aggregates. This suggests a possible layer-by-layer formation process of onion-like nanostructures in the ISM. These findings reveal the possible role of carbon nanoparticles as selective catalysts that could provide reaction substrates for the formation of interstellar PAHs, high-fullerenes and soots from gas-phase molecules.

DArk Matter Particle Explorer (DAMPE) is a space-borne apparatus for detecting the high-energy cosmic-rays like electrons, $\gamma$-rays, protons and heavy-ions. Plastic Scintillator Detector (PSD) is the top-most sub-detector of the DAMPE. The PSD is designed to measure the charge of incident high-energy particles and it also serves as a veto detector for discriminating $\gamma$-rays from charged particles. In this paper, PSD on-orbit calibration procedure is described, which includes five steps of pedestal, dynode correlation, response to minimum-ionizing particles (MIPs), light attenuation function and energy reconstruction. A method for reconstructing the charge of incident high energy cosmic-ray particles is introduced. The detection efficiency of each PSD strip is verified to be above 99.5%, the total efficiency of the PSD for charged particles is above 99.99%.

The mineralogy of mare basalts reflects the chemical composition of the magma source, as well as the physical and chemical environment of the rock's formation. It is significant for understanding the thermal evolution of the Moon. In this study, the spatial distribution of the mineralogy of the lunar northern nearside basalts was mapped using the Moon Mineralogy Mapper (M3) data. The study area, which is an elongated mare, Mare Frigoris and northern Mare Imbrium, was mapped and characterized into 27 units based on multi-source data, including spectrum, terrain, and element abundance. We extracted 177 M3 spectra from fresh craters. Spectral parameters such as absorption center and band area ratio (BAR) were obtained through data processing. The variation of mafic mineral of this region was acquired by analyzing these parameters. The basaltic units in eastern Mare Frigoris, which are older, have been found to be dominated by clinopyroxene with lower CaO compared to the returned lunar samples; this is similar to older basaltic units in Mare Imbrium. The basaltic units of western Mare Frigoris and Sinus Roris which are younger have been found to be rich in olivine. The late-stage basalts in Oceanus Procellarum and Mare Imbrium show the same feature. These wide-spread olivine-rich basalts suggest the uniqueness in the evolution of the Moon. Geographically speaking, Mare Frigoris is an individual mare, but the magma source region have connections with surrounding maria in consideration of mineral differences between western and eastern Frigoris, as well as mineral similarities with maria at the same location.

We report on a search for Weakly Interacting Massive Particles (WIMPs) using 278.8 days of data collected with the XENON1T experiment at LNGS. XENON1T utilizes a liquid xenon time projection chamber with a fiducial mass of $(1.30 \pm 0.01)$ t, resulting in a 1.0 t$\times$yr exposure. The energy region of interest, [1.4, 10.6] $\mathrm{keV_{ee}}$ ([4.9, 40.9] $\mathrm{keV_{nr}}$), exhibits an ultra-low electron recoil background rate of $(82\substack{+5 \\ -3}\textrm{ (sys)}\pm3\textrm{ (stat)})$ events/$(\mathrm{t}\times\mathrm{yr}\times\mathrm{keV_{ee}})$. No significant excess over background is found and a profile likelihood analysis parameterized in spatial and energy dimensions excludes new parameter space for the WIMP-nucleon spin-independent elastic scatter cross-section for WIMP masses above 6 GeV/c${}^2$, with a minimum of $4.1\times10^{-47}$ cm$^2$ at 30 GeV/c${}^2$ and 90% confidence level.

We investigated 64 pairs of interacting-CME events identified by the simultaneous observations of SOHO and STEREO spacecraft from 2010 January to 2014 August, to examine the relationship between the large SEP events in the energy of 25-60MeV and the properties of the interacting CMEs. We found that during CME interactions the large SEP events in this study were all generated by CMEs with the presence of enhanced type II radio bursts, which also have wider longitudinal distributions comparing to events with the absence of type II radio burst or its enhancement (almost associated with small SEP events). It seems that the signature of type II radio bursts enhancement is a good discriminator between large SEP and small or none SEP event producers during CME interactions. The type II radio burst enhancement is more likely to be generated by CME interactions, with the main CME having larger speed (v), angular width (WD), mass (m) and kinetic energy (E_k), that taking over the preceding CMEs which also have higher v, WD, m and E_k, than those preceding CMEs in CME pairs missing the type II radio bursts or enhancements. Our analysis also revealed that the intensities of associated SEP events correlate negatively with the intersection height of the two CMEs. Most of type-II-enhanced events and SEP events are coincidentally and almost always made by the fast and wide main CMEs that sweeping fully over a relatively slower and narrower preceding CMEs. We suggest that a fast CME with enough energy completely overtaking a relatively narrower preceding CME, especially in low height, can drive more energetic shock signified by the enhanced type II radio bursts. The shock may accelerate ambient particles and lead to large SEP event more easily.

We present the results from our Atacama Large Millimeter/submillimeter Array (ALMA) imaging observations of the CO (7-6), [CI] 370 um (hereafter [CI]) and [NII] 205 um (hereafter [NI]I) lines and their underlying continuum emission of BRI 1335-0417, an infrared bright quasar at z = 4.407. At the achieved resolutions of 1.1" to 1.2" (or 7.5 to 8.2 kpc), the continuum at 205 and 372 um (rest-frame), the CO (7-6), and the [CI] emissions are at best barely resolved whereas the [NII] emission is well resolved with an ALMA beam de-convolved major axis of 1.3" (+/- 0.3") or 9 (+/-2) kpc. As a warm dense gas tracer, the CO (7-6) emission shows a more compact spatial distribution and a significantly higher peak velocity dispersion than the other two lines that probe lower density gas, a picture favoring a merger-triggered star formation (SF) scenario over an orderly rotating SF disk. The CO (7-6) data also indicate a possible QSO-driven gas outflow that reaches a maximum line-of-sight velocity of 500 to 600 km/s. The far-infrared (FIR) dust temperature (T_dust) of 41.5 K from a gray-body fit to the continuum agrees well with the average T_dust inferred from various line luminosity ratios. The resulting L_CO(7-6)/L_FIR luminosity ratio is consistent with that of local luminous infrared galaxies powered predominantly by SF. The CO(7-6) luminosity-inferred SF rate is 5.1 (+/-1.5) x 10^3 M_solar/yr . The system has an effective star-forming region of 1.7 (+1.7/-0.8) kpc in diameter and a molecular gas reservoir of ~5 x 10^{11} M_solar.

DAMPE (DArk Matter Particle Explorer) is a spaceborne high-energy cosmic ray and gamma-ray detector, successfully launched in December 2015. It is designed to probe astroparticle physics in the broad energy range from few GeV to 100 TeV. The scientific goals of DAMPE include the identification of possible signatures of Dark Matter annihilation or decay, the study of the origin and propagation mechanisms of cosmic-ray particles, and gamma-ray astronomy. DAMPE consists of four sub-detectors: a plastic scintillator strip detector, a Silicon-Tungsten tracKer-converter (STK), a BGO calorimeter and a neutron detector. The STK is composed of six double layers of single-sided silicon micro-strip detectors interleaved with three layers of tungsten for photon conversions into electron-positron pairs. The STK is a crucial component of DAMPE, allowing to determine the direction of incoming photons, to reconstruct tracks of cosmic rays and to estimate their absolute charge (Z). We present the in-flight performance of the STK based on two years of in-flight DAMPE data, which includes the noise behavior, signal response, thermal and mechanical stability, alignment and position resolution.

By a method of population synthesis we construct a model of dark matter (DM) accretion onto binary black holes (BHs), and investigate the merger rate of the binary BHs. We find that the merger rate can weakly increase (less than 10\%). However, the DM accretion can efficiently enhance the masses of binary BHs. In our model, the result for $Z = 0.01$ without the DM accretion can not explain the GW170104, GW170814, and GW150914, while with the DM accretion, it can cover all observations well. For the higher metallicity ($Z = 0.02$), our model can not explain the mergers of high-mass binary BHs like GW170104, GW170814, and GW150914. We estimate that the merger rate of binary BHs lies between $\rm55Gpc^{-3}yr^{-1}$ to $\rm197Gpc^{-3}yr^{-1}$.

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

Tibet is known as the third pole of the earth, as high as the South Pole and North Pole. The Ngari (Ali) observatory in Tibet has the advantage of plenty of photometric night, low precipitable water vapor, high transmittance, good seeing. It is a good site, and promising to be one of the best place for infrared and submillimeter observations in the world. However, there is no data available for sky background brightness in such place, which restrict the astronomical development of the sites. In the near infrared band of J, H, Ks, a NIR sky brightness monitor (NISBM) is designed based on InGaAs photoelectric diode. By using the method of chopper modulation and digital lock-in amplifier processing, the SNR (Signal Noise Ratio), detectivity and the data acquisition speed of the device is greatly improved. For each band of J, H, Ks, an independent instrument is designed and calibrated in laboratory. The NISBM has been installed in Ngari observatory in July of 2017 and obtained the first data of NIR sky brightness at Ngari observatory.

MARST (Multi-channel Antarctic Solar Telescope) aims to observe the Sun in multiple wavebands in Antarctica and will be China's first solar telescope in Antarctica. The telescope has two tubes, corresponding to Photosphere observation which uses 11 filters and Chromosphere observation in H{\alpha} waveband. The two tubes need to be coordinated to observe at the same time. The telescope will observe the Sun for a long time, so a self-guiding module is needed to improve sun tracking. Besides, performing solar specific flat-field exposure is necessary for analyzing. EPICS is introduced to control each hardware and an autonomous observation system based on RTS2 is designed under such demands. EPICS application modules are implemented for each device: telescope mount & focuser, filter wheel, H{\alpha} filter, dome with webcams, Andor CCD and PI CCD. We also integrate EPICS modules into RTS2 framework with an XML format configuration. To control these applications autonomously we have developed a RTS2 executor module where two plan classes are instantiated to control two sets of filters and CCDs, and to ensure only one could control the mount at the same time. Different types of observation plans are designed to describe different series of processes with different priorities. To improve sun tracking, we calculate the centroid of each image to get the offset, then apply the correction to the telescope during observation process. For frontend users, a GUI based on PyQt5 and QML is implemented and connects to rts2-httpd and rts2-proxy modules so that users can control devices, check images and get logs.

A near-infrared (NIR) measurement based on digital orthogonal vector lock-in amplifier (LIA) is present in this paper. NIR sky background radiation is very weak. To detect the signals obscured by noise, the best way achieved is to use a chopper to modulate the detected signal and using a LIA to demodulate. The effect of 1/f noise of detector, dark current and other noises can be reduced to get sufficient signal-to-noise ratio (SNR). The orthogonal vector LIA can avoid the phase shift on the accuracy of measurement by two orthogonal components. In order to simplify the system, a digital algorithm is adopted to realize the LIA which is operated in a microcontroller with ARM cortex-M4. Data is obtained through ADC and the signal of detector is amplified and filtered. Then the phase sensitive detection (PSD), low-pass filter (LPF) and amplitude phase calculation are performed. The digital method can greatly simplify the circuit, and conveniently adjust the time constant of the LPF to realize the different equivalent noise bandwidth (ENB). The algorithm has the specification of high precision, flexible usage, simple implementation and low computation resource. By using this method, the weak infrared signal submerged by the noise can be obtained, which extremely improves the detection capability of the system.

Using the stellar evolution code---Modules for Experiments in Stellar Astrophysics (MESA), we investigate the evolution of massive stars with different rotational velocities and metallicities towards Wolf-Rayet stars. In our simulations, the initial rotating velocities are taken as 0, 250, 500 and 650 km s$^{-1}$, and the metallicities equal to 0.02, 0.014, 0.008, 0.006, 0.004 and 0.002. We show our rapid rotation models in the HR diagram compared with the observations. We find that the rotational mixing is less efficient at high metallicity, and these stars become Wolf-Rayet (WR) stars when the helium in their center is ignited. However, rapid rotating massive stars at low metallicity can easily evolve into WR stars due to the rotation resulted in chemically homogeneous evolution. This can explain the origin of single WR stars in galaxy at low metallicity. In our models, the observed SMC WR stars are consistent with the single-star evolution models. However at higher metallicities our single-star evolution models can only explain the luminous, hydrogen-rich WN stars and O stars (which are classified as WR stars previously).

Very faint X-ray transients (VFXTs) are a group of X-ray binaries with low luminosities, the peak X-ray luminosities during their outbursts being only 10$^{34}$--10$^{36}$ erg\,s$^{-1}$. Using the $\gamma$-ray data obtained with the Large Area Telescope (LAT) onboard {\it Fermi Gamma-Ray Space Telescope}, we investigate their possible nature of containing rotation-powered pulsars, or more specifically being transitional millisecond pulsars (MSPs). Among more than 40 known VFXTs, we select 12 neutron star systems. We analyze the LAT data for the fields of the 12 VFXTs in 0.2--300 GeV energy range, but do not find any counterparts likely detected by {\it Fermi}. We obtain the luminosity upper limits for the 12 sources. While the distances to the sources are largely uncertain, the upper limits are comparable to the luminosities of two transitional systems PSR J1023$-$0038 and XSS J12270$-$4859. From our study, we conclude that no evidence is found at $\gamma$-rays for the suggestion that some of VFXTs could contain rotation-powered MSPs (or be transitional MSP systems).

We report the results from our analysis of 8 years of the data for the gamma-ray binary LS I +61 303, obtained with the Large Area Telescope onboard the Fermi Gamma-Ray Space Telescope. We find a significant dip around the binary's periastron in the superorbital light curves, and by fitting the light curves with a sinusoidal function, clear phase shifts are obtained. The superorbital modulation seen in the binary has been long known and different scenarios have been proposed. Based on our results, we suggest that the circumstellar disk around the Be companion of this binary has an elliptical shape and the major axis of the disk rotates at the superorbital period of 1667 days. As a result, the density of the ambient material around the compact star of the binary changes along the binary orbit over the superorbital period, causing the phase shifts in the modulation, and around periastron, the compact star probably enters the disk, causing the appearance of the dip. Numerical simulations may be conducted in order to study the detailed physical processes and verify our suggestion.

Pulsar-like compact stars provide us a unique laboratory to explore properties of dense matter at supra-nuclear densities. One of the models for pulsar-like stars is that they are totally composed of "strangeons", and in this paper we studied the pulsar glitches in a strangeon star model. Strangeon stars would be solidified during cooling, and the solid stars would be natural to have glitches as the result of starquakes. Based on the starquake model established before, we proposed that when the starquake occurs, the inner motion of the star which changes the moment of inertia and has impact on the glitch sizes, is divided into plastic flow and elastic motion. The plastic flow which is induced in the fractured part of the outer layer, would move tangentially to redistribute the matter of the star and would be hard to recover. The elastic motion, on the other hand, changes its shape and would recover significantly. Under this scenario, we could understand the behaviors of glitches without significant energy releasing, including the Crab and the Vela pulsars, in an uniform model. We derive the recovery coefficient as a function of glitch size, as well as the time interval between two successive glitches as the function of the released stress. Our results show consistency with observational data under reasonable ranges of parameters. The implications on the oblateness of the Crab and the Vela pulsars are discussed.

In this paper, we describe the XENON100 data analyses used to assess the target-intrinsic background sources radon ($^{222}$Rn), thoron ($^{220}$Rn) and krypton ($^{85}$Kr). We detail the event selections of high-energy alpha particles and decay-specific delayed coincidences. We derive distributions of the individual radionuclides inside the detector and quantify their abundances during the main three science runs of the experiment over a period of $\sim$ 4 years, from January 2010 to January 2014. We compare our results to external measurements of radon emanation and krypton concentrations where we find good agreement. We report an observed reduction in concentrations of radon daughters that we attribute to the plating-out of charged ions on the negatively biased cathode.

The Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) is the largest existing spectroscopic survey telescope, having 32 scientific charge-coupled-device (CCD) cameras for receiving spectra. Stability and automation of the camera-control software are essential, but cannot be provided by the existing system. The Remote Telescope System 2nd Version (RTS2) is an open-source and automatic observatory-control system. However, all previous RTS2 applications have concerned small telescopes. This paper focuses on implementation of an RTS2-based camera-control system for the 32 CCDs of LAMOST. A virtual camera module inherited from the RTS2 camera module is built as a device component working on the RTS2 framework. To improve the controllability and robustness, a virtualized layer is designed using the master-slave software paradigm, and the virtual camera module is mapped to the 32 real cameras of LAMOST. The new system is deployed in the actual environment and experimentally tested. Finally, multiple observations are conducted using this new RTS2-framework-based control system. The new camera-control system is found to satisfy the requirements for LAMOST automatic camera control. This is the first time that RTS2 has been applied to a large telescope, and provides a referential solution for full RTS2 introduction to the LAMOST observatory control system.

We continue to investigate the binary system Kepler-16, consisting of a K-type main-sequence star, a red dwarf, and a circumbinary Saturnian planet. As part of our study, we describe the system's habitable zone based on different climate models. We also report on stability investigations for possible Earth-mass Trojans while expanding a previous study by B. L. Quarles and collaborators given in 2012. For the climate models we carefully consider the relevance of the system's parameters. Furthermore, we pursue new stability simulations for the Earth-mass objects starting along the orbit of Kepler-16b. The eccentricity distribution as obtained prefers values close to circular, whereas the inclination distribution remains flat. The stable solutions are distributed near the co-orbital Lagrangian points, thus enhancing the plausibility that Earth-mass Trojans might be able to exist in the Kepler-16(AB) system.

We report on the response of liquid xenon to low energy electronic recoils below 15 keV from beta decays of tritium at drift fields of 92 V/cm, 154 V/cm and 366 V/cm using the XENON100 detector. A data-to-simulation fitting method based on Markov Chain Monte Carlo is used to extract the photon yields and recombination fluctuations from the experimental data. The photon yields measured at the two lower fields are in agreement with those from literature; additional measurements at a higher field of 366 V/cm are presented. The electronic and nuclear recoil discrimination as well as its dependence on the drift field and photon detection efficiency are investigated at these low energies. The results provide new measurements in the energy region of interest for dark matter searches using liquid xenon.

The unprecedented volume and rate of transient events that will be discovered by the Large Synoptic Survey Telescope (LSST) demands that the astronomical community update its followup paradigm. Alert-brokers -- automated software system to sift through, characterize, annotate and prioritize events for followup -- will be critical tools for managing alert streams in the LSST era. The Arizona-NOAO Temporal Analysis and Response to Events System (ANTARES) is one such broker. In this work, we develop a machine learning pipeline to characterize and classify variable and transient sources only using the available multiband optical photometry. We describe three illustrative stages of the pipeline, serving the three goals of early, intermediate and retrospective classification of alerts. The first takes the form of variable vs transient categorization, the second, a multi-class typing of the combined variable and transient dataset, and the third, a purity-driven subtyping of a transient class. While several similar algorithms have proven themselves in simulations, we validate their performance on real observations for the first time. We quantitatively evaluate our pipeline on sparse, unevenly sampled, heteroskedastic data from various existing observational campaigns, and demonstrate very competitive classification performance. We describe our progress towards adapting the pipeline developed in this work into a real-time broker working on live alert streams from time-domain surveys.

Non-attractor inflation is known as the only single field inflationary scenario that can violate non-Gaussianity consistency relation with the Bunch-Davies vacuum state and generate large local non-Gaussianity. However, it is also known that the non-attractor inflation by itself is incomplete and should be followed by a phase of slow-roll attractor. Moreover, there is a transition process between these two phases. In the past literature, this transition was approximated as instant and the evolution of non-Gaussianity in this phase was not fully studied. In this paper, we follow the detailed evolution of the non-Gaussianity through the transition phase into the slow-roll attractor phase, considering different types of transition. We find that the transition process has important effect on the size of the local non-Gaussianity. We first compute the net contribution of the non-Gaussianities at the end of inflation in canonical non-attractor models. If the curvature perturbations keep evolving during the transition - such as in the case of smooth transition or some sharp transition scenarios - the $\mathcal{O}(1)$ local non-Gaussianity generated in the non-attractor phase can be completely erased by the subsequent evolution, although the consistency relation remains violated. In extremal cases of sharp transition where the super-horizon modes freeze immediately right after the end of the non-attractor phase, the original non-attractor result can be recovered. We also study models with non-canonical kinetic terms, and find that the transition can typically contribute a suppression factor in the squeezed bispectrum, but the final local non-Gaussianity can still be made parametrically large.

The hybrid Tibet AS array was successfully constructed in 2014. It has 4500 m$^{2}$ underground water Cherenkov pools used as the muon detector (MD) and 789 scintillator detectors covering 36900 m$^{2}$ as the surface array. At 100 TeV, cosmic-ray background events can be rejected by approximately 99.99\%, according to the full Monte Carlo (MC) simulation for $\gamma$-ray observations. In order to use the muon detector efficiently, we propose to extend the surface array area to 72900 m$^{2}$ by adding 120 scintillator detectors around the current array to increase the effective detection area. A new prototype scintillator detector is developed via optimizing the detector geometry and its optical surface, by selecting the reflective material and adopting dynode readout. This detector can meet our physics requirements with a positional non-uniformity of the output charge within 10\% (with reference to the center of the scintillator), time resolution FWHM of $\sim$2.2 ns, and dynamic range from 1 to 500 minimum ionization particles.

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

We report the first dark matter search results from XENON1T, a $\sim$2000-kg-target-mass dual-phase (liquid-gas) xenon time projection chamber in operation at the Laboratori Nazionali del Gran Sasso in Italy and the first ton-scale detector of this kind. The blinded search used 34.2 live days of data acquired between November 2016 and January 2017. Inside the (1042$\pm$12) kg fiducial mass and in the [5, 40] $\mathrm{keV}_{\mathrm{nr}}$ energy range of interest for WIMP dark matter searches, the electronic recoil background was $(1.93 \pm 0.25) \times 10^{-4}$ events/(kg $\times$ day $\times \mathrm{keV}_{\mathrm{ee}}$), the lowest ever achieved in a dark matter detector. A profile likelihood analysis shows that the data is consistent with the background-only hypothesis. We derive the most stringent exclusion limits on the spin-independent WIMP-nucleon interaction cross section for WIMP masses above 10 GeV/c${}^2$, with a minimum of 7.7 $\times 10^{-47}$ cm${}^2$ for 35-GeV/c${}^2$ WIMPs at 90% confidence level.

We study the power spectrum of quasi-single field inflation where strong coupling is considered. The contribution from the massive propagator can be divided into local and non-local contributions. The local one is the leading contribution and is power-law suppressed as a function of mass, while the non-local contribution is exponentially suppressed in the large mass limit. For the local contribution, it is possible to use the effective field theory approach to study the power spectrum in the strongly coupled region of the parameter space. For the non-local contribution, we developed a partial effective field theory method to simplify the calculation: When there are multiple massive propagators, one can fully compute it after integrating out all but one massive propagator by effective field theory. The result retains the "standard clock" signal, which is interesting for probing the expansion history of the primordial universe and the physics of a "cosmological collider". The error involved compared to the full calculation is power law suppressed by the effective mass of the heavy field.

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.

We derive fitting formulae for the quick determination of the existence of S-type and P-type habitable zones in binary systems. Based on previous work, we consider the limits of the climatological habitable zone in binary systems (which sensitively depend on the system parameters) based on a joint constraint encompassing planetary orbital stability and a habitable region for a possible system planet. Additionally, we employ updated results on planetary climate models obtained by Kopparapu and collaborators. Our results are applied to four P-type systems (Kepler-34, Kepler-35, Kepler-413, and Kepler-1647) and two S-type systems (TrES-2 and KOI-1257). Our method allows to gauge the existence of climatological habitable zones for these systems in a straightforward manner with detailed consideration of the observational uncertainties. Further applications may include studies of other existing systems as well as systems to be identified through future observational campaigns.

We report new searches for the solar axions and galactic axion-like dark matter particles, using the first low-background data from PandaX-II experiment at China Jinping Underground Laboratory, corresponding to a total exposure of about $2.7\times 10^4$ kg$\cdot$day. No solar axion or galactic axion-like dark matter particle candidate has been identified. The upper limit on the axion-electron coupling ($g_{Ae}$) from the solar flux is found to be about $4.35 \times 10^{-12}$ in mass range from $10^{-5}$ to 1 keV/$c^2$ with 90\% confidence level, similar to the recent LUX result. We also report a new best limit from the $^{57}$Fe de-excitation. On the other hand, the upper limit from the galactic axions is on the order of $10^{-13}$ in the mass range from 1 keV/$c^2$ to 10 keV/$c^2$ with 90\% confidence level, slightly improved compared with the LUX.

Fiber bursts are a type of fine structure that is often superposed on type IV radio continuum emission during solar flares. Although studied for many decades, its physical exciter, emission mechanism, and association with the flare energy release remain unclear, partly due to the lack of simultaneous imaging observations. We report the first dynamic spectroscopic imaging observations of decimetric fiber bursts, which occurred during the rise phase of a long-duration eruptive flare on 2012 March 3, as obtained by the Karl G. Jansky Very Large Array in 1--2 GHz. Our results show that the fiber sources are located near and above one footpoint of the flare loops. The fiber source and the background continuum source are found to be cospatial and share the same morphology. It is likely that they are associated with nonthermal electrons trapped in the converging magnetic fields near the footpoint, as supported by a persistent coronal hard X-ray source present during the flare rise phase. We analyze three groups of fiber bursts in detail with dynamic imaging spectroscopy and obtain their mean frequency-dependent centroid trajectories in projection. By using a barometric density model and magnetic field based on a potential-field extrapolation, we further reconstruct the 3-D source trajectories of fiber bursts, for comparison with expectations from the whistler wave model and two MHD-based models. We conclude that the observed fiber burst properties are consistent with an exciter moving at the propagation velocity expected for whistler waves, or models that posit similar exciter velocities.

We investigate the long-term spectral variability in the ultra-luminous X-ray source Holmberg IX X--1. By analyzing the data from eight {\it Suzaku} and 13 {\it XMM-Newton} observations conducted between 2001 and 2015, we perform a detailed spectral modeling for all spectra with simple models and complex physical models. We find that the spectra can be well explained by a disc plus thermal Comptonization model. Applying this model, we unveil correlations between the X-ray luminosity ($L_{\rm X}$) and the spectral parameters. Among the correlations, a particular one is the statistically significant positive correlation between $L_{\rm X}$ and the photon index ($\Gamma$), while at the high luminosities of $> 2\times10^{40}\,{\rm~erg\ s}^{-1}$, the source becomes marginally hard and that results a change in the slope of the $\Gamma - L_{\rm X}$ correlation. Similar variability behavior is observed in the optical depth of the source around $L_{\rm X} \sim 2\times10^{40}\,{\rm~erg\ s}^{-1}$ as the source becomes more optically thick. We consider the scenario that a corona covers the inner part of the disc, and the correlations can be explained as to be driven by the variability of seed photons from the disc input into the corona. On the basis of the disc-corona model, we discuss the physical processes that are possibly indicated by the variability of the spectral parameters. Our analysis reveals the complex variability behavior of Holmberg IX X--1 and the variability mechanism is likely related to the geometry of the X-ray emitting regions.

The liquid drop model of 2-flavored ($u$ and $d$) nucleus is well known and successful, analogically, a similar drop model for 3-flavored ($u$, $d$ and $s$) nucleus is developed. A 3-flavored nucleus conjectured could be stable only if its baryon number is lager than a critical one, $A_{\rm c}$, in which strangeons are the constituent as an analogy of nucleons for nucleus. We try to model strangeon matter in a sense of phenomenological liquid drop, with two free parameters: the mass per bayron of a strangeon in vacuum, $M$, and potential deep between strangeons, $\epsilon$. It is found that, for $M\sim$ GeV and $\epsilon\sim 100$ MeV, strangeon matter could be stable and its critical number could be as low as $A_{\rm c}=300$.

The DArk Matter Particle Explorer (DAMPE), one of the four scientific space science missions within the framework of the Strategic Pioneer Program on Space Science of the Chinese Academy of Sciences, is a general purpose high energy cosmic-ray and gamma-ray observatory, which was successfully launched on December 17th, 2015 from the Jiuquan Satellite Launch Center. The DAMPE scientific objectives include the study of galactic cosmic rays up to $\sim 10$ TeV and hundreds of TeV for electrons/gammas and nuclei respectively, and the search for dark matter signatures in their spectra. In this paper we illustrate the layout of the DAMPE instrument, and discuss the results of beam tests and calibrations performed on ground. Finally we present the expected performance in space and give an overview of the mission key scientific goals.

1SWASP\,J162117.36$+$441254.2 was originally classified as an EW-type binary with a period of 0.20785\,days. However, it was detected to have undergone a stellar outburst on June 3, 2016. Although the system was latter classified as a cataclysmic variable (CV) and the event was attributed as a dwarf-nova outburst, the physical reason is still unknown. This binary has been monitored photometrically since April 19, 2016 and many light curves were obtained before, during and after the outburst. Those light and color curves observed before the outburst indicate that the system is a special CV. The white dwarf is not accreting material from the secondary and there are no accretion disks surrounding the white dwarf. By comparing the light curves obtained from April 19 to September 14, 2016, it was found that magnetic activity of the secondary is associated with the outburst. We show strong evidence that the $L_1$ region on the secondary was heavily spotted before and after the outburst and thus quench the mass transfer, while the outburst is produced by a sudden mass accretion of the white dwarf. These results suggest that J162117 is a good astrophysical laboratory to study stellar magnetic activity and its influences on CV mass transfer and mass accretion.

We present results of searches for vector and pseudo-scalar bosonic super-WIMPs, which are dark matter candidates with masses at the keV-scale, with the XENON100 experiment. XENON100 is a dual-phase xenon time projection chamber operated at the Laboratori Nazionali del Gran Sasso. A profile likelihood analysis of data with an exposure of 224.6 live days $\times$ 34\,kg showed no evidence for a signal above the expected background. We thus obtain new and stringent upper limits in the $(8-125)$\,keV/c$^2$ mass range, excluding couplings to electrons with coupling constants of $g_{ae} > 3\times10^{-13}$ for pseudo-scalar and $\alpha'/\alpha > 2\times10^{-28}$ for vector super-WIMPs, respectively. These limits are derived under the assumption that super-WIMPs constitute all of the dark matter in our galaxy.

Using the Modules for Experiments in Stellar Astrophysics code, we investigate the influences of irradiation on ultra-compact X-ray binary (UCXB) evolution. Although the persistent UCXBs have short orbital periods which result in high irradiation flux, the irradiation hardly affects the evolution of persistent sources because the WDs in these binaries have large masses which lead to very low irradiation depth. The irradiation has a significant effect on the transient sources during outburst phase. At the beginning of the outburst, high X-ray luminosity produces high radiation flux, which results in the significant expansion of WD. Then, the irradiation triggers high mass-transfer rates, which can last several days for the transient sources with WDs whose masses are larger than $\sim0.015 M_\odot$ or several hundred years for these sources with WDs whose masses are less than $\sim0.012 M_\odot$. The observed three persistent UCXBs, XTE J0929-314, 4U 1916-05 and SWIFT J1756.9-2508, may belong to the latter.

We collect a sample of 78 Fermi detected blazars with broad line region (BLR) data and the quasi-simultaneous multi-wave band data. By analyzing the sample, we find: (1) For whole blazar sample, there exist significant correlations between radio, $\gamma$-ray luminosity and BLR luminosity. The slope of the best-fit linear regression equation is close to the theoretical value 1. Our results provided another support of jet-disk symbiosis. (2) For FSRQ sample, we find a significant correlation between radio luminosity and optical luminosity, and weaker correlations between X-ray luminosity and radio, optical luminosity. It indicates that the X-ray band radiation mechanism is slightly different from radio and optical bands. The X-ray emission is diluted by inverse Compton scattering (IC). There is a significant correlation between BLR luminosity and $\gamma$-ray luminosity, but no correlations between X-ray luminosity and BLR, $\gamma$-ray luminosity, which suggests that the seed soft photons are different for IC process. From these results we suggest that the X-ray band emission is diluted by synchrotron self-Compton (SSC) process, $\gamma$-ray band emission is produced by external Compton (EC) process.

The XENON1T experiment at the Laboratori Nazionali del Gran Sasso (LNGS) is the first WIMP dark matter detector operating with a liquid xenon target mass above the ton-scale. Out of its 3.2t liquid xenon inventory, 2.0t constitute the active target of the dual-phase time projection chamber. The scintillation and ionization signals from particle interactions are detected with low-background photomultipliers. This article describes the XENON1T instrument and its subsystems as well as strategies to achieve an unprecedented low background level. First results on the detector response and the performance of the subsystems are also presented.

Bolometric luminosity is a basic physical parameter that is widely used in the study of blazars. Due to the lack of simultaneous full wavelength data, several estimators of the bolometric luminosity are being used in practice. In this paper, we study and evaluate the reliability and significance of six estimators, the $5GHz$ luminosity, the 1keV luminosity, the $\gamma$-ray luminosity, the $5GHz$ luminosity+the 1keV luminosity, the $5GHz$ luminosity+the $\gamma$-ray luminosity and the 1keV luminosity+the $\gamma$-ray luminosity, by analyzing the linear correlations between the integrated bolometric luminosity and them. Our main results are as follows. (i) All the six estimators are reliable in the sense that they are all significant correlated with the bolometric luminosity. (ii) Ranking from the higher significance of the reliability to lower one the six estimators are the the $5GHz$ luminosity+the $\gamma$-ray luminosity, the 1keV luminosity+the $\gamma$-ray luminosity, $\gamma$-ray luminosity, the $5GHz$ luminosity+the 1keV luminosity, the $5GHz$ luminosity and the 1keV luminosity. (iii) We suggest that the bolometric luminosity can be well estimated by the $\gamma$-ray luminosity using the best linear equation that given in this paper for Fermi FSRQs. (iv) According to the linear regressions obtained in the analysis, we provide calibration for each estimator.

We report on WIMP search results in the XENON100 detector using a non-relativistic effective field theory approach. The data from science run II (34 kg $\times$ 224.6 live days) was re-analyzed, with an increased recoil energy interval compared to previous analyses, ranging from $(6.6 - 240)~\mathrm{keV_\mathrm{nr}}$. The data is found to be compatible with the background-only hypothesis. We present 90% confidence level exclusion limits on the coupling constants of WIMP-nucleon effective operators using a binned profile likelihood method. We also consider the case of inelastic WIMP scattering, where incident WIMPs may up-scatter to a higher mass state, and set exclusion limits on this model as well.

We present the first constraints on the spin-dependent, inelastic scattering cross section of Weakly Interacting Massive Particles (WIMPs) on nucleons from XENON100 data with an exposure of 7.64$\times$10$^3$\,kg\,day. XENON100 is a dual-phase xenon time projection chamber with 62\,kg of active mass, operated at the Laboratori Nazionali del Gran Sasso (LNGS) in Italy and designed to search for nuclear recoils from WIMP-nucleus interactions. Here we explore inelastic scattering, where a transition to a low-lying excited nuclear state of $^{129}$Xe is induced. The experimental signature is a nuclear recoil observed together with the prompt de-excitation photon. We see no evidence for such inelastic WIMP-$^{129}$Xe interactions. A profile likelihood analysis allows us to set a 90\% C.L. upper limit on the inelastic, spin-dependent WIMP-nucleon cross section of $3.3 \times 10^{-38}$\,cm$^{2}$ at 100\,GeV/c$^2$. This is the most constraining result to date, and sets the pathway for an analysis of this interaction channel in upcoming, larger dual-phase xenon detectors.

Following our work presented in Dai et al. (2016), we report our detailed data analysis for another 38 Fermi gamma-ray un-associated sources. These sources are selected from the Fermi Large Area Telescope (LAT) third source catalog on the basis of the properties of known gamma-ray millisecond pulsars (MSPs) and for the purpose of finding likely candidate MSPs. From our analysis of the LAT data, we identify that among the 38 sources, 28 of them are single point-like sources with clean background and their spectra show significant curvature. We also conduct analysis of archival X-ray data available for 24 of the 28 sources. In the fields of 10 sources, there are at least one X-ray object, and in those of the other 14 sources, no X-ray object is detected but probably due to the X-ray observations being short. We discuss the possible MSP nature for these sources. Six of them(J0514.6-4406, J1035.7-6720, J1624.2-4041, J1744.1-7619, J1946.4-5403, and J2039.6-5618) are most likely associated with pulsars because of multi-wavelength identifications including direct radio or gamma-ray detection of pulsations. To firmly establish the associations or verify the MSP nature for other sources, deep X-ray and/or optical observations are needed.

The multi-TeV $\gamma$-rays from the Galactic Center (GC) have a cutoff at tens of TeV, whereas the diffuse emission has no such cutoff, which is regarded as an indication of PeV proton acceleration by the HESS experiment. It is important to understand the inconsistency and study the possibility that PeV cosmic-ray acceleration could account for the apparently contradictory point and diffuse $\gamma$-ray spectra. In this work, we propose that the cosmic rays are accelerated up to $>$PeV in GC. The interaction between cosmic rays and molecular clouds is responsible for the multi-TeV $\gamma$-ray emissions from both the point source and diffuse sources today. Enhanced by the small volume filling factor (VFF) of the clumpy structure, the absorption of the $\gamma$-rays leads to a sharp cutoff spectrum at tens of TeV produced in the GC. Away from galactic center, the VFF grows and the absorption enhancement becomes negligible. As a result, the spectra of $\gamma$-ray emissions for both point source and diffuse sources can be successfully reproduced under such self-consistent picture. In addition, a "surviving-tail" at $\sim$100 TeV is expected from the point source, which can be observed by future projects CTA and LHAASO. Neutrinos are simultaneously produced during proton-proton (PP) collision. With 5-10 years observations, the KM3Net experiment will be able to detect the PeV source according to our calculation.

Aqua MODIS, unlike its predecessor on board the Terra spacecraft, had always been thought to have been spared from significant deleterious impacts of electronic crosstalk on its imagery. However, recent efforts brought to our attention the presence of striping artifacts in Aqua MODIS images from band 24 (4.47$\mu$m), which upon further inspection proved to have a noticeable impact on the quality of the L1B product and to have been present since the beginning of the mission, in 2002. Using images of the Moon from scheduled lunar observations, we linked the artifacts with electronic crosstalk contamination of the response of detector 1 of band 24 by signal sent from the detector 10 of band 26 (1.375$\mu$m), a neighboring band in the same focal plane assembly. In this paper, we report on these findings, the artifact mitigation strategy adopted by us, and on our success in restoring band 24 detector 1 behavior and image quality.

he DArk Matter Particle Explorer (DAMPE) is a general purposed satellite-borne high energy $\gamma-$ray and cosmic ray detector, and among the scientific objectives of DAMPE are the searches for the origin of cosmic rays and an understanding of Dark Matter particles. As one of the four detectors in DAMPE, the Plastic Scintillator Detector (PSD) plays an important role in the particle charge measurement and the photons/electrons separation. The PSD has 82 modules, each consists of a long organic plastic scintillator bar and two PMTs at both ends for readout, in two layers and covers an overall active area larger than 82 cm $\times$ 82 cm. It can identify the charge states for relativistic ions from H to Fe, and the detector efficiency for Z=1 particles can reach 0.9999. The PSD has been successfully launched with DAMPE on Dec. 17, 2015. In this paper, the design, the assembly, the qualification tests of the PSD and some of the performance measured on the ground have been described in detail.

New constraints are presented on the spin-dependent WIMP-nucleon interaction from the PandaX-II experiment, using a data set corresponding to a total exposure of 3.3$\times10^4$ kg-days. Assuming a standard axial-vector spin-dependent WIMP interaction with $^{129}$Xe and $^{131}$Xe nuclei, the most stringent upper limits on WIMP-neutron cross sections for WIMPs with masses above 10 GeV/c$^{2}$ are set in all dark matter direct detection experiments. The minimum upper limit of $4.1\times 10^{-41}$ cm$^2$ at 90\% confidence level is obtained for a WIMP mass of 40 GeV/c$^{2}$. This represents more than a factor of two improvement on the best available limits at this and higher masses. These improved cross-section limits provide more stringent constraints on the effective WIMP-proton and WIMP-neutron couplings.

The Giant Radio Array for Neutrino Detection (GRAND) is a planned array of ~200 000 radio antennas deployed over ~200 000 km2 in a mountainous site. It aims primarly at detecting high-energy neutrinos via the observation of extensive air showers induced by the decay in the atmosphere of taus produced by the interaction of cosmic neutrinos under the Earth surface. GRAND aims at reaching a neutrino sensitivity of 5.10$^{11}$ E$^{-2}$ GeV$^{-1}$cm$^{-2}$s$^{-1}$sr$^{-1}$ above 3.10$^{16}$ eV. This ensures the detection of cosmogenic neutrinos in the most pessimistic source models, and ~50 events per year are expected for the standard models. The instrument will also detect UHECRs and possibly FRBs. Here we show how our preliminary design should enable us to reach our sensitivity goals, and discuss the steps to be taken to achieve GRAND.

We investigate the effects of the core-collapse supernova ejecta on a rapidly rotating and massive companion star. We show that the stripped mass raises by twice when compare with a massive but non-rotating companion star. In close binaries with orbital periods of about 1 day, the stripped masses reach up to $\sim 1 M_\odot$. By simulating the evolutions of the rotational velocities of the massive companion stars based on different stripped masses, we find that the rotational velocity decreases greatly for stripped mass that is higher than about $1 M_\odot$. Of all the known high mass X-ray binaries (HMXBs), Cygnus X-3 and 1WGA J0648.024418 have the shortest orbital periods of 0.2 and 1.55 days, respectively. The optical counterpart of the former is a Wolf-Rayet star, whereas it is a hot subdwarf for the latter. Applying our model to the two HMXBs, we suggest that the hydrogen-rich envelopes of their optical counterparts may have been stripped by CCSN ejecta.

We report the identification of seven transient X-ray sources in the nearby Magellanic-type galaxy NGC 4449 using the archival multi-epoch X-ray observations conducted with {\it Chandra}, {\it XMM-Newton} and {\it Swift} telescopes over year 2001--2013. Among them, two sources are classified as supersoft X-ray sources (SSSs) because of their soft X-ray color and rest of the sources are X-ray binaries (XRBs). Transient SSSs spectra can be fitted with a blackbody of effective temperature $\sim 80-105$ eV and luminosities were $\simeq 10^{37} - 10^{38} {\rm~erg\ s}^{-1}$ in 0.3--8 keV. These properties are consistent with the widely accepted model for SSSs, an accreting white dwarf with the steady nuclear burning on its surface, while the SSS emission has also been observed in many post-nova systems. Detailed analysis of one sufficiently bright SSS revealed the strong short-term variability, possibly showing a 2.3 hour periodic modulation, and long-term variability, detectable over 23 years with different X-ray telescopes before year 2003. The X-ray properties of four other transients are consistent with neutron star or black hole binaries in their hard state, while the remaining source is most likely an XRB with a quasi-soft X-ray spectrum. Analysis of archival {\it Hubble Space Telescope} image data was also conducted, and multiple massive stars were found as possible counterparts. We conclude that the X-ray transient properties in NGC 4449 are similar to those in other Magellanic-type galaxies.

We explore the astrobiological significance of F-type stars of spectral type between F5 V and F9.5 V, which possess Jupiter-type planets within or close to their climatological habitable zones. These planets, or at least a subset of those, may also possess rocky exomoons, which potentially offer habitable environments. Our work considers eight selected systems. The Jupiter-type planets in these systems are in notably differing orbits with eccentricities between 0.08 (about Mars) and 0.72. We consider the stellar UV environments provided by the photospheric stellar radiation, which allows us to compute the UV habitable zones for the systems. Following previous studies, DNA is taken as a proxy for carbon-based macromolecules following the paradigm that extraterrestrial biology might be based on hydrocarbons. We found that the damage inflicted on DNA is notably different for the range of systems studied, and also varies according to the orbit of the Jupiter-type planet, especially for systems of high ellipticity, as expected. For some systems excessive values of damage are attained compared to today's Earth or during the Archean eon. Considering that the detection of exomoons around different types of stars will remain challenging in the foreseeable future, we view our work also as an example and template for investigating the combined requirements of climatological and UV-based habitability for exosolar objects.

The Arizona-NOAO Temporal Analysis and Response to Events System (ANTARES) is a joint effort of NOAO and the Department of Computer Science at the University of Arizona to build prototype software to process alerts from time-domain surveys, especially LSST, to identify those alerts that must be followed up immediately. Value is added by annotating incoming alerts with existing information from previous surveys and compilations across the electromagnetic spectrum and from the history of past alerts. Comparison against a knowledge repository of properties and features of known or predicted kinds of variable phenomena is used for categorization. The architecture and algorithms being employed are described.

We work with a UV conformal U(1)' extension of the Standard Model, motivated by the hierarchy problem and recent collider anomalies. This model admits fermionic vector portal WIMP dark matter charged under the U(1)' gauge group. The asymptotically safe boundary conditions can be used to fix the coupling parameters, which allows the observed thermal relic abundance to constrain the mass of the dark matter particle. This highly restricts the parameter space, allowing strong predictions to be made. The parameter space of several UV conformal U(1)' scenarios will be explored, and both bounds and possible signals from direct and indirect detection observation methods will be discussed.

The LAMP (Lightweight Asymmetry and Magnetism Probe) X-ray telescope is a mission concept to measure the polarization of X-ray astronomical sources at 250 eV via imaging mirrors that reflect at incidence angles near the polarization angle, i.e., 45 deg. Hence, it will require the adoption of multilayer coatings with a few nanometers d-spacing in order to enhance the reflectivity. The nickel electroforming technology has already been successfully used to fabricate the high angular resolution imaging mirrors of the X-ray telescopes SAX, XMM-Newton, and Swift/XRT. We are investigating this consolidated technology as a possible technique to manufacture focusing mirrors for LAMP. Although the very good reflectivity performances of this kind of mirrors were already demonstrated in grazing incidence, the reflectivity and the scattering properties have not been tested directly at the unusually large angle of 45 deg. Other possible substrates are represented by thin glass foils or silicon wafers. In this paper we present the results of the X-ray reflectivity campaign performed at the BEAR beamline of Elettra - Sincrotrone Trieste on multilayer coatings of various composition (Cr/C, Co/C), deposited with different sputtering parameters on nickel, silicon, and glass substrates, using polarized X-rays in the spectral range 240 - 290 eV.

We present initial work about attaining fitting formulas for the quick determination of the existence of S-type and P-type habitable zones in binary systems. Following previous work, we calculate the limits of the climatological habitable zone in binary systems (which sensitively depend on the system parameters) based on a joint constraint encompassing planetary orbital stability and a habitable region for a possible system planet. We also consider updated results on planetary climate models previously obtained by Kopparapu and collaborators. Fitting equations based on our work are presented for selected cases.

We fit the spectral energy distributions (SEDs) of members of a large sample of Fermi 2LAC blazars to synchrotron and inverse Compton (IC) models. Our main results are as follows. (i) As suggested by previous works, the correlation between peak frequency and curvature can be explained by statistical or stochastic particle acceleration mechanisms. For BL Lacs, we find a linear correlation between synchrotron peak frequency and its curvature. The slope of the correlation is consistent with the stochastic acceleration mechanisms and confirm previous studies. For FSRQs, we also find a linear correlation, but its slope cannot be explained by previous theoretical models. (ii) We find a significant correlation between IC luminosity and synchrotron luminosity. The slope of the correlation of FSRQs is consistent with the EC process. And the slope of the correlation of BL Lac is consistent with the SSC process. (iii) We find several significant correlations between IC curvature and several basic parameters of blazars (black hole mass, broad line luminosity, the Lorentz factor of jet). We also find significant correlations between bolometric luminosity and these basic parameters of blazars which suggest that the origin of jet is a mixture of the mechanisms proposed by Blandford $\&$ Znajek and by Blandford $\&$ Payne.

Nested multilayer mirrors are commonly used in X-ray telescope structures to increase the collecting area. To balance the difficulty and cost of producing these mirrors, the classical Wolter-I structure has previously been replaced with a conical Wolter-I structure, but it can lead to significantly poorer angular resolution. In this paper, we consider changing one of the mirror shapes (paraboloid or hyperboloid) of the Wolter-I structure to a conical mirror shape, while the other mirror shape remains a quadric surface-type structure, which could thus ensure the imaging quality. Based on a theoretical analysis and the results of simulations, the cone-hyperboloid and paraboloid-cone structures are found to provide the best angular resolutions, and the cone-hyperboloid structure is shorter than the paraboloid-cone structure at the focal length. The cone-hyperboloid structure is then nested to obtain the best on-axis angular resolution and off-axis images.

We study the effects of the non-attractor initial conditions for the canonical single-field inflation. The non-attractor stage can last only several $e$-folding numbers, and should be followed by hilltop inflation. This two-stage evolution leads to large scale suppression in the primordial power spectrum, which is favored by recent observations. Moreover we give a detailed calculation of primordial non-Guassianity due to the "from non-attractor to slow-roll" transition, and find step features in the local and equilateral shapes. We conclude that a plateau-like inflaton potential with an initial non-attractor phase yields interesting features in both power spectrum and bispectrum.

eXTP is a science mission designed to study the state of matter under extreme conditions of density, gravity and magnetism. Primary targets include isolated and binary neutron stars, strong magnetic field systems like magnetars, and stellar-mass and supermassive black holes. The mission carries a unique and unprecedented suite of state-of-the-art scientific instruments enabling for the first time ever the simultaneous spectral-timing-polarimetry studies of cosmic sources in the energy range from 0.5-30 keV (and beyond). Key elements of the payload are: the Spectroscopic Focusing Array (SFA) - a set of 11 X-ray optics for a total effective area of about 0.9 m^2 and 0.6 m^2 at 2 keV and 6 keV respectively, equipped with Silicon Drift Detectors offering <180 eV spectral resolution; the Large Area Detector (LAD) - a deployable set of 640 Silicon Drift Detectors, for a total effective area of about 3.4 m^2, between 6 and 10 keV, and spectral resolution <250 eV; the Polarimetry Focusing Array (PFA) - a set of 2 X-ray telescope, for a total effective area of 250 cm^2 at 2 keV, equipped with imaging gas pixel photoelectric polarimeters; the Wide Field Monitor (WFM) - a set of 3 coded mask wide field units, equipped with position-sensitive Silicon Drift Detectors, each covering a 90 degrees x 90 degrees FoV. The eXTP international consortium includes mostly major institutions of the Chinese Academy of Sciences and Universities in China, as well as major institutions in several European countries and the United States. The predecessor of eXTP, the XTP mission concept, has been selected and funded as one of the so-called background missions in the Strategic Priority Space Science Program of the Chinese Academy of Sciences since 2011. The strong European participation has significantly enhanced the scientific capabilities of eXTP. The planned launch date of the mission is earlier than 2025.

Neutron stars are ideal astrophysical laboratories for testing theories of the de Haas-van Alphen (dHvA) effect and diamagnetic phase transition which is associated with magnetic domain formation. The "magnetic interaction" between delocalized magnetic moments of electrons (the Shoenberg effect), can result in an effect of the diamagnetic phase transition into domains of alternating magnetization (Condon's domains). Associated with the domain formation are prominent magnetic field oscillation and anisotropic magnetic stress which may be large enough to fracture the crust of magnetar with a super-strong field. Even if the fracture is impossible as in "low-field" magnetar, the depinning phase transition of domain wall motion driven by low field rate (mainly due to the Hall effect) in the randomly perturbed crust can result in a catastrophically variation of magnetic field. This intermittent motion, similar to the avalanche process, makes the Hall effect be dissipative. These qualitative consequences about magnetized electron gas are consistent with observations of magnetar emission, and especially the threshold critical dynamics of driven domain wall can partially overcome the difficulties of "low-field" magnetar bursts and the heating mechanism of transient, or "outbursting" magnetar.

According to the nova model from \citet{Yaron2005} and \citet{Jose1998} and using Monte Carlo simulation method, we investigate the contribution of chemical abundances in nova ejecta to the interstellar medium (ISM) of the Galaxy. We find that the ejected mass by classical novae (CNe) is about $2.7\times10^{-3}$ $ \rm M_\odot\ {\rm yr^{-1}}$. In the nova ejecta, the isotopic ratios of C, N and O, that is, $^{13}$C/$^{12}$C, $^{15}$N/$^{14}$N and $^{17}$O/$^{16}$O, are higher about one order of magnitude than those in red giants. We estimate that about 10$\%$, 5$\%$ and 20$\%$ of $^{13}$C, $^{15}$N and $^{17}$O in the ISM of the Galaxy come from nova ejecta, respectively. However, the chemical abundances of C, N and O calculated by our model can not cover all of observational values. This means that there is still a long way to go for understanding novae.

We report our analysis of the Fermi Large Area Telescope data for 39 millisecond pulsars (MSPs) listed in the second $\gamma$-ray pulsar catalog. Spectra of the pulsars are obtained. We fit the spectra with a function of a power law with exponential cutoff, and find the best-fit parameters of photon index $\Gamma = 1.54^{+0.10}_{-0.11}$ and cutoff energy $E_{c} = 3.70^{+0.95}_{-0.70}$ GeV. This spectral shape, which includes the intrinsic differences in the spectra of the MSPs, can be used for finding candidate MSPs and unidentified types of sources detected by Fermi at high Galactic latitudes. In one of the MSPs PSR J0614-3329, we find significant pulsed emission upto 200 GeV. The result has thus added this MSP to the group of the Crab and Vela pulsars that have been detected with >50 GeV pulsed emission. Comparing the $\gamma$-ray spectrum of PSR J0614-3329 with those of the Crab and Vela pulsars, we discuss possible emission mechanisms for the very high-energy component.

On 2012 March 7, the STEREO Ahead and Behind spacecraft, along with the near-earth spacecraft (e.g. SOHO, Wind) situated between the two STEREO spacecraft, observed an extremely large global solar energetic particle (SEP) event in Solar Cycle 24. Two successive coronal mass ejections (CMEs) have been detected close in time. From the multi-point in-situ observations, it can be found that this SEP event was caused by the first CME, and the second one was not involved. Using the velocity dispersion analysis (VDA), we find that for well magnetically connected point, the energetic protons and electrons are released nearly at the same time. The path lengths to STEREO-B(STB) of protons and electrons have distinct difference and deviate remarkably from the nominal Parker spiral path length, which is likely due to the presence of interplanetary magnetic structures situated between the source and the STB. Also the VDA method seems only to obtain reasonable results at well-connected locations and the inferred energetic particles release times in different energy channels are similar. We suggest that good-connection is crucial for obtaining both accurate release time and path length simultaneously, agreeing with the modeling result of Wang & Qin (2015).

We report the results from our detailed analysis of the Fermi Large Area Telescope (LAT) data for the pulsar binary PSR B1259-63/LS 2883. During the GeV flares that occurred when the pulsar was in the periastron passages, we have detected a 5--300 GeV component at $\simeq 5\sigma$ in emission from the binary. The detection verifies the presence of the component that has been marginally found in the previous studies of the binary. Furthermore, we have discovered that this component was present even in the quiescent state of the binary, specifically the mean anomaly phase 0.7--0.9. The component can be described by a power law with photon index $\Gamma\sim 1.4$, and the flux in the flares is approximately one order of magnitude higher than that in quiescence. We discuss the origin of this component. It likely arises from the inverse-Compton process: high-energy particles from the shock, due to the interaction between the winds from the pulsar and massive companion, scatter the seed photons from the companion to GeV/TeV energies. Based on this scenario, model fits to the broad-band X-ray--to--TeV spectra of the binary in the flaring and quiescent states are provided.

It was recently proposed that a giant flare of the blazar PKS B1424-418 at redshift $z=1.522$ is in association with a PeV-energy neutrino event detected by IceCube. Based on this association we here suggest that the flight time difference between the PeV neutrino and gamma-ray photons from blazar flares can be used to constrain the violations of equivalence principle (EP) and the Lorentz invariance for neutrinos. From the calculated Shapiro delay due to clusters or superclusters in the nearby universe, we find that violation of the equivalence principle for neutrinos and photons is constrained to an accuracy of at least $10^{-5}$, which is two orders of magnitude tighter than the constraint placed by MeV neutrinos from supernova 1987A. Lorentz invariance violation (LIV) arises in various quantum-gravity theories, which predicts an energy-dependent velocity of propagation in vacuum for photons and neutrinos. We find that the association of the PeV neutrino with the gamma-ray outburst set limits on the energy scale of possible LIV to $> 0.01 E_{pl}$ for linear LIV models and $>6\times 10^{-8} E_{pl}$ for quadratic order LIV models, where $E_{pl}$ is the Planck energy scale. These are the most stringent constraints on neutrino LIV for subluminal neutrinos.

We report the likely detection of gamma-ray emission from the northeast shell region of the historical supernova remnant (SNR) SN 1006. Having analyzed 7 years of Fermi LAT Pass 8 data for the region of SN 1006, we found a GeV gamma-ray source detected with 4 sigma significance. Both the position and spectrum of the source match those of HESS J1504-418 respectively, which is TeV emission from SN 1006. Considering the source as the GeV gamma-ray counterpart to SN~1006, the broadband spectral energy distribution is found to be approximately consistent with the leptonic scenario that has been proposed for the TeV emission from the SNR. Our result has likely confirmed the previous study of the SNRs with TeV shell-like morphology: SN 1006 is one of them sharing very similar peak luminosity and spectral shape.

We present the spectral and temporal properties of 15 candidate transient X-ray sources detected in archival {\it XMM-Newton} and {\it Chandra} observations of the nearby Magellanic-type, SB(s)m galaxy NGC 55. On the basis of an X-ray color classification scheme, the majority of the sources may be identified as X-ray binaries (XRBs), and six sources are soft, including a likely supernova remnant. We perform a detailed spectral and variability analysis of the data for two bright candidate XRBs. Both sources displayed strong short-term X-ray variability, and their X-ray spectra and hardness ratios are consistent with those of XRBs. These results, combined with their high X-ray luminosities ($\sim 10^{38}~\rm erg~s^{-1}$), strongly suggest that they are black hole (BH) binaries. Seven less luminous sources have spectral properties consistent with those of neutron star or BH XRBs in both normal and high-rate accretion modes, but one of them is the likely counterpart to a background galaxy (because of positional coincidence). From the spectral analysis, we find that the six soft sources are candidate super soft sources (SSSs), with dominant emission in the soft (0.3--2 keV) X-ray band. Archival {\it Hubble Space Telescope} optical images for 7 sources are available, and the data suggest that most of them are likely to be high-mass XRBs. Our analysis has revealed the heterogeneous nature of the transient population in NGC 55 (6 high-mass XRBs, 1 low-mass XRBs, 6 SSSs, 1 active galactic nucleus), helping establish the similarity of the X-ray properties of this galaxy to those of other Magellanic-type galaxies.

In recent decades, inorganic crystals have been widely used in dark matter direct search experiments. To contribute to the understanding of the capabilities of $CsI(Na)$ and $CaF_{2}(Eu)$ crystals, a mono-energetic neutron beam is utilized to study the properties of nuclear recoils, which are expected to be similar to signals of dark matter direct detection. The quenching factor of nuclear recoils in $CsI(Na)$ and $CaF_{2}(Eu)$, as well as an improved discrimination factor between nuclear recoils and ${\gamma}$ backgrounds in $CsI(Na)$, are reported.

We report our detailed data analysis for 39 $\gamma$-ray sources selected from the 992 unassociated sources in the \textit{Fermi} Large Area Telescope (LAT) third source catalog. The selection criteria, which were set for finding candidate millisecond pulsars (MSPs), are non-variables with curved spectra and $>$5$^{\circ}$ Galactic latitudes. From our analysis, 24 sources were found to be point-like sources not contaminated by background or nearby unknown sources. Three of them, J1544.6$-$1125, J1625.1$-$0021, and J1653.6$-$0158, have been previously studied, indicating that they are likely MSPs. The spectra of J0318.1+0252 and J2053.9+2922 do not have properties similar to that of known $\gamma$-ray MSPs, and we thus suggest that they are not MSPs. Analysis of archival X-ray data for most of the 24 sources were also conducted. Four sources were found with X-ray objects in their error circles, and 16 with no detection. The ratios between the $\gamma$-ray fluxes and X-ray fluxes or flux upper limits are generally lower than those of the known $\gamma$-ray MSPs, suggesting that if the $\gamma$-ray sources are MSPs, none of the X-ray objects are the counterparts. Deep X-ray or radio observations of these sources are needed in order to identify their MSP nature.

Gravitational waves from coalescences of neutron stars or stellar-mass black holes into intermediate-mass black holes (IMBHs) of $\gtrsim 100$ solar masses represent one of the exciting possible sources for advanced gravitational-wave detectors. These sources can provide definitive evidence for the existence of IMBHs, probe globular-cluster dynamics, and potentially serve as tests of general relativity. We analyse the accuracy with which we can measure the masses and spins of the IMBH and its companion in intermediate-mass ratio coalescences. We find that we can identify an IMBH with a mass above $100 ~ M_\odot$ with $95\%$ confidence provided the massive body exceeds $130 ~ M_\odot$. For source masses above $\sim200 ~ M_\odot$, the best measured parameter is the frequency of the quasi-normal ringdown. Consequently, the total mass is measured better than the chirp mass for massive binaries, but the total mass is still partly degenerate with spin, which cannot be accurately measured. Low-frequency detector sensitivity is particularly important for massive sources, since sensitivity to the inspiral phase is critical for measuring the mass of the stellar-mass companion. We show that we can accurately infer source parameters for cosmologically redshifted signals by applying appropriate corrections. We investigate the impact of uncertainty in the model gravitational waveforms and conclude that our main results are likely robust to systematics.

The magnetar 4U~0142+61 has been well studied at optical and infrared wavelengths and is known to have a complicated broad-band spectrum over the wavelength range. Here we report the result from our linear imaging polarimetry of the magnetar at optical $I$-band. From the polarimetric observation carried out with the 8.2-m Subaru telescope, we determine the degree of linear polarization $P=1.0\pm$3.4\%, or $P\leq$5.6\% (90\% confidence level). Considering models suggested for optical emission from magnetars, we discuss the implications of our result. The upper limit measurement indicates that different from radio pulsars, magnetars probably would not have strongly polarized optical emission if the emission arises from their magnetosphere as suggested.

The Lightweight Asymmetry and Magnetism Probe (LAMP) is a micro-satellite mission concept dedicated for astronomical X-ray polarimetry and is currently under early phase study. It consists of segmented paraboloidal multilayer mirrors with a collecting area of about 1300 cm^2 to reflect and focus 250 eV X-rays, which will be detected by position sensitive detectors at the focal plane. The primary targets of LAMP include the thermal emission from the surface of pulsars and synchrotron emission produced by relativistic jets in blazars. With the expected sensitivity, it will allow us to detect polarization or place a tight upper limit for about 10 pulsars and 20 blazars. In addition to measuring magnetic structures in these objects, LAMP will also enable us to discover bare quark stars if they exist, whose thermal emission is expected to be zero polarized, while the thermal emission from neutron stars is believed to be highly polarized due to plasma polarization and the quantum electrodynamics (QED) effect. Here we present an overview of the mission concept, its science objectives and simulated observational results.

As the evolutionary link between the radio millisecond pulsars (MSPs) and the low mass X-ray binaries or intermediate mass X-ray binaries, the millisecond X-ray pulsars (MSXPs) are important objects in testing theories of pulsar formation and evolution. In general, neutron stars in MSXPs can form via core collapse supernova (CC channel) of massive stars or accretion induced collapse (AIC channel) of an accreting ONeMg WD whose mass reaches the Chandrasekhar limit. Here, in addition to CC and AIC channels we also consider another channel, i.e., evolution induced collapse (EIC channel) of a helium star with mass between $1.4$ and $2.5 M_{\odot}$. Using a population synthesis code, we have studied MSXPs arising from three different evolutionary channels. We find that the Galactic birthrates of transient MSXPs and persistent MSXPs are about 0.7---$1.4\times 10^{-4}$ yr$^{-1}$. Our population synthesis calculations have shown that about 50\%---90\% of the MSXPs have undergone CC channel, about 10\%---40\% of them have undergone EIC channel, and the MSXPs via AIC channel are the least.

In gamma ray astronomy, the energy range from sub-100GeV to TeV is crucial due to where there is a gap between space experiments and ground-based ones. In addition, observations in this energy range are expected to provide more details about the high energy emission from GRBs,and thus to understand EBL. Based on the observation results and the related knowledge, scientists may be able to unveil the mysteries of galaxy formation and the evolution of early universe. One of the principal issues for next generation Imaging Atmospheric Cherenkov Telescopes (IACT) is to achieve larger field of view (FoV). In this work, we report a refractive water convex lens as light collector to test the feasibility of a new generation of IACT, and some preliminary test results on the optical properties (the focal length, spot size, transmittance, etc.) of a 0.9 m diameter water lens, the photodetectors and DAQ system of a prototype are presented and discussed.

The X-ray binary 4U 1543$-$624 has been provisionally identified as an ultracompact system with an orbital period of $\simeq$18~min. We have carried out time-resolved optical imaging of the binary to verify the ultra-short orbital period. Using 140\,min of high-cadence $r'$-band photometry we recover the previously-seen sinusoidal modulation and determine a period $P=18.20\pm0.09$\,min. In addition, we also see a 7.0$\times 10^{-4}$\,mag\,min$^{-1}$ linear decay, likely related to variations in the source's accretion activity. Assuming that the sinusoidal modulation arises from X-ray heating of the inner face of the companion star, we estimate a distance of 6.0--6.7\,kpc and an inclination angle of 34$^{\circ}$--61$^{\circ}$ (90\% confidence) for the binary. Given the stability of the modulation we can confirm that the modulation is orbital in origin and 4U 1543$-$624 is an ultracompact X-ray binary.

We report the results of a weakly-interacting massive particle (WIMP) dark matter search using the full 80.1\;live-day exposure of the first stage of the PandaX experiment (PandaX-I) located in the China Jin-Ping Underground Laboratory. The PandaX-I detector has been optimized for detecting low-mass WIMPs, achieving a photon detection efficiency of 9.6\%. With a fiducial liquid xenon target mass of 54.0\,kg, no significant excess event were found above the expected background. A profile likelihood analysis confirms our earlier finding that the PandaX-I data disfavor all positive low-mass WIMP signals reported in the literature under standard assumptions. A stringent bound on the low mass WIMP is set at WIMP mass below 10\,GeV/c$^2$, demonstrating that liquid xenon detectors can be competitive for low-mass WIMP searches.

High-energy neutrino astronomy will probe the working of the most violent phenomena in the Universe. The Giant Radio Array for Neutrino Detection (GRAND) project consists of an array of $\sim10^5$ radio antennas deployed over $\sim$200000km$^2$ in a mountainous site. It aims at detecting high-energy neutrinos via the measurement of air showers induced by the decay in the atmosphere of $\tau$ leptons produced by the interaction of the cosmic neutrinos under the Earth surface. Our objective with GRAND is to reach a neutrino sensitivity of $3\times10^{-11}E^{-2}$GeV$^{-1}$cm$^{-2}$s$^{-1}$sr$^{-1}$ above $3 \times10^{16}$eV. This sensitivity ensures the detection of cosmogenic neutrinos in the most pessimistic source models, and about 100 events per year are expected for the standard models. GRAND would also probe the neutrino signals produced at the potential sources of UHECRs. We show how our preliminary design should enable us to reach our sensitivity goals, and present the experimental characteristics. We assess the possibility to adapt GRAND to other astrophysical radio measurements. We discuss in this token the technological options for the detector and the steps to be taken to achieve the GRAND project.

The Spitzer-Cosmic Assembly Deep Near-Infrared Extragalactic Legacy Survey (S-CANDELS; PI G. Fazio) is a Cycle 8 Exploration Program designed to detect galaxies at very high redshifts (z > 5). To mitigate the effects of cosmic variance and also to take advantage of deep coextensive coverage in multiple bands by the Hubble Space Telescope Multi-Cycle Treasury Program CANDELS, S-CANDELS was carried out within five widely separated extragalactic fields: the UKIDSS Ultra-Deep Survey, the Extended Chandra Deep Field South, COSMOS, the HST Deep Field North, and the Extended Groth Strip. S-CANDELS builds upon the existing coverage of these fields from the Spitzer Extended Deep Survey (SEDS) by increasing the integration time from 12 hours to a total of 50 hours but within a smaller area, 0.16 square degrees. The additional depth significantly increases the survey completeness at faint magnitudes. This paper describes the S-CANDELS survey design, processing, and publicly-available data products. We present IRAC dual-band 3.6+4.5 micron catalogs reaching to a depth of 26.5 AB mag. Deep IRAC counts for the roughly 135,000 galaxies detected by S-CANDELS are consistent with models based on known galaxy populations. The increase in depth beyond earlier Spitzer/IRAC surveys does not reveal a significant additional contribution from discrete sources to the diffuse Cosmic Infrared Background (CIB). Thus it remains true that only roughly half of the estimated CIB flux from COBE/DIRBE is resolved.

The origin of dust in a galaxy is poorly understood. Recently, the surveys of the Large Magellanic Cloud (LMC) provide astrophysical laboratories for the dust studies. By a method of population synthesis, we investigate the contributions of dust produced by asymptotic giant branch (AGB) stars, common envelope (CE) ejecta and type II supernovae (SNe II) to the total dust budget in the LMC. Based on our models, the dust production rates (DPRs) of AGB stars in the LMC are between about $2.5\times10^{-5}$ and $4.0\times10^{-6}M_\odot{\rm yr^{-1}}$. The uncertainty mainly results from different models for the dust yields of AGB stars. The DPRs of CE ejecta are about $6.3\times10^{-6}$(The initial binary fraction is 50\%). These results are within the large scatter of several observational estimates. AGB stars mainly produce carbon grains, which is consistent with the observations. Most of dust grains manufactured by CE ejecta are silicate and iron grains. The contributions of SNe II are very uncertain. Compared with SNe II without reverse shock, the DPRs of AGB stars and CE ejecta are negligible. However, if only 2 \% of dust grains produced by SNe II can survive after reverse shock, the contributions of SNe II are very small. The total dust masses produced by AGB stars in the LMC are between $2.8\times10^4$ and $3.2\times10^5M_\odot$, and those produced by CE ejecta are about $6.3\times10^4$. They are much lower than the values estimated by observations. Therefore, there should be other dust sources in the LMC.

Providing an early warning of galactic supernova explosions from neutrino signals is important in studying supernova dynamics and neutrino physics. A dedicated supernova trigger system has been designed and installed in the data acquisition system at Daya Bay and integrated into the worldwide Supernova Early Warning System (SNEWS). Daya Bay's unique feature of eight identically-designed detectors deployed in three separate experimental halls makes the trigger system naturally robust against cosmogenic backgrounds, enabling a prompt analysis of online triggers and a tight control of the false-alert rate. The trigger system is estimated to be fully sensitive to 1987A-type supernova bursts throughout most of the Milky Way. The significant gain in sensitivity of the eight-detector configuration over a mass-equivalent single detector is also estimated. The experience of this online trigger system is applicable to future projects with spatially distributed detectors.

We report our discovery of orbitally modulated $\gamma$-ray emission from the black widow system PSR J1311-3430. We analyze the \textit{Fermi} Large Area Telescope data during the offpulse phase interval of the pulsar, and find the orbital modulation signal at a $\sim$3$\sigma$ confidence level. Further spectral analysis shows no significant differences for the spectra obtained during the bright and faint orbital phase ranges. A simple sinusoid-like function can describe the modulation. Given these properties, we suggest that the intrabinary $\gamma$-ray emission arises from the region close to the companion and the modulation is caused by the occultation of the emitting region by the companion, similar to that is seen in the transitional millisecond pulsar binary (MSP) PSR J1023+0038. Considering the X-ray detection of intrabinary shock emission from eclipsing MSP binaries recently reported, this discovery further suggests the general existence of intrabinary $\gamma$-ray emission from them.

We report on the likely detection of gamma-ray emission from the supernova remnant (SNR) Kesteven 27 (Kes 27). We analyze 5.7 yr Fermi Large Area Telescope data of the SNR region and find an unresolved source at a position consistent with the radio brightness peak and the X-ray knot of Kes 27, which is located in the eastern region of the SNR and caused by the interaction with a nearby HI cloud. The source's emission is best fit with a power-law spectrum with a photon index of 2.5 +/- 0.1 and a >0.2 GeV luminosity of 5.8*10^34 erg s^-1 assuming a distance of 4.3 kpc, as derived from radio observations of the nearby HI cloud. Comparing the properties of the source with that of other SNRs that are known to be interacting with nearby high-density clouds, we discuss the origin of the source's emission. The spectral energy distribution of the source can be described by a hadronic model that considers the interaction of energetic protons, escaping from the shock front of Kes 27, with a high-density cloud.

Gamma-ray bursts (GRBs) are proposed as candidate sources of ultra-high energy cosmic rays (UHECRs). We study the possibility that the PeV neutrinos recently observed by IceCube are produced by GRB cosmic rays interacting with the interstellar gas in the host galaxies. By studying the relation between the X-ray absorption column density N_H and the surface star-formation rate of GRB host galaxies, we find that N_H is a good indicator of the surface gas density of the host galaxies. Then we are able to calculate the neutrino production efficiency of CRs for GRBs with known N_H. We collect a sample of GRBs that have both measurements of N_H and accurate gamma-ray fluence, and attempt to calculate the accumulated neutrino flux based on the current knowledge about GRBs and their host galaxies. When the CR intensity produced by GRBs is normalized with the observed UHECR flux above $10^{19}{\rm eV}$, the accumulated neutrino flux at PeV energies is estimated to be about $(0.3\pm0.2)\times10^{-8} \rm{GeV\ cm^{-2}\ s^{-1}\ sr^{-1}} $ (per flavor) under the assumption that GRB energy production rate follows the cosmic star-formation rate and the favorable assumption about the CR diffusion coefficient. This flux is insufficient to account for the IceCube observations, but the estimate suffers from some assumptions in the calculation and thus we can not rule out this scenario at present.

We have simulated the Expanded Owens Valley Solar Array (EOVSA) radio images generated at multiple frequencies from a model solar active region, embedded in a realistic solar disk model, and explored the resulting datacube for different spectral analysis schemes to evaluate the potential for realizing one of EOVSA's most important scientific goals--coronal magnetography. In this paper, we focus on modeling the gyroresonance and free-free emission from an on-disk solar active region model with realistic complexities in electron density, temperature and magnetic field distribution. We compare the magnetic field parameters extrapolated from the image datacube along each line of sight after folding through the EOVSA instrumental profile with the original (unfolded) parameters used in the model. We find that even the most easily automated, image-based analysis approach (Level 0) provides reasonable quantitative results, although they are affected by systematic effects due to finite sampling in the Fourier (uv) plane. Finally, we note the potential for errors due to misidentified harmonics of the gyrofrequency, and discuss the prospects for applying a more sophisticated spectrally-based analysis scheme (Level 1) to resolve the issue in cases where improved uv coverage and spatial resolution are available.

We report the Fermi Large Area Telescope (LAT) detection of a $\gamma$-ray source at the position of SAX J1808.4$-$3658. This transient low-mass X-ray binary contains an accreting millisecond puslar, which is only seen during its month-long outbursts and likely switches to be rotation powered during its quiescent state. Emission from the $\gamma$-ray source can be described by a power law with an exponential cutoff, the characteristic form for pulsar emission. Folding the source's 2.0--300 GeV photons at the binary orbital period, a weak modulation is seen (with an H-test value of $\sim$17). In addition, three sets of archival XMM-Newton data for the source field are analyzed, and we find only one X-ray source with 3--4$\sigma$ flux variations in the 2$\sigma$ error circle of the $\gamma$-ray source. However based on the X-ray properties, this X-ray source is not likely a background AGN, the major class of Fermi sources detected by LAT. These results support the possible association between the $\gamma$-ray source and SAX J1808.4$-$3658 and thus the scenario that the millisecond pulsar is rotation powered in the quiescent state. Considering a source distance of 3.5 kpc for SAX J1808.4$-$3658, the 0.1--300 GeV luminosity is 5.7$\times 10^{33}$ erg s$^{-1}$, implying a $\gamma$-ray conversion efficiency of 63\% for the pulsar in this binary.

We report the outbursts from a newly discovered X-ray transient in the Magellanic-type, SB(s)m galaxy NGC 55. The transient source, XMMU J001446.81-391123.48, was undetectable in the 2001 \xmm{} and 2004 \chandra{} observations, but detected in a 2010 \xmm{} observation at a significance level of $9\sigma$ in the 0.3--8 keV energy band. The \xmm{} spectrum is consistent with a power law with photon index $\Gamma = 3.17^{+0.22}_{-0.20}$, but is better fit with a $kT_{in} = 0.70\pm0.06$\,keV disk blackbody. The luminosity was $\sim 10^{38}$\,erg\,s$^{-1}$, and the source displayed strong short-term X-ray variability. These results, combined with the hardness ratios of its emission, strongly suggest an X-ray binary nature for the source. The follow-up studies with \swift{} XRT observations revealed that the source exhibited recurrent outbursts with period about a month. The XRT spectra can be described by a power law ($\Gamma\sim 2.5$--2.9) or a disk blackbody ($kT_{in}\sim 0.8$--1.0\,keV), and the luminosity was in a range of 10$^{38}$--10$^{39}$\,erg\,s$^{-1}$, with no evidence showing any significant changes of the spectral parameters in the observations. Based on the X-ray spectral and temporal properties, we conclude that XMMU J001446.81-391123.48 is a new transient X-ray binary in NGC 55, which possibly contains a black hole primary.

We investigate the local stability of strongly magnetized relativistic tori orbiting Kerr black holes, for the case of a purely toroidal magnetic field topology. Our approach encompasses both solving the full dispersion relation numerically and solving a simplified analytic treatment. Both approaches indicate that such tori are subject to an unstable non-axisymmetric magnetorotational mode, regardless of such background properties as the gas-to-magnetic pressure ratio and specific angular momentum distribution. We demonstrate that our modal analysis matches well with what is seen in global, three-dimensional, general relativistic magnetohydrodynamic simulations.

We present CO $J$=2-1 observations towards 32 nearby gas-rich star-forming galaxies selected from the ALFALFA and WISE catalogs, using the Sub-millimeter Telescope. Our sample is selected to be dominated by intermediate-$M_{\rm *}$ galaxies. The scaling-relations between molecular gas, atomic gas and galactic properties (stellar mass, NUV$- r$ and WISE color W3$-$W2) are examined and discussed. Our results show that (1). In the galaxies with stellar mass $M_{\rm *}$ $\leqslant 10^{10}$ $M_{\odot}$, HI fraction ($f_{\rm HI}$ $\equiv$ $M_{\rm HI}$/$M_{\rm *}$) is significantly higher than that of more massive galaxies, while H$_2$ gas fraction ($f_{\rm H_2}$ $\equiv$ $M_{\rm H_2}$/$M_{\rm *}$) remain nearly unchanged. (2). Comparing with $f_{\rm H_2}$, $f_{\rm HI}$ correlates better with both $M_{\rm *}$ and NUV$- r$. (3). A new parameter, WISE color W3$-$W2 (12\,$\mu$m$-$4.6\,$\mu$m) is introduced, which is similar to NUV$- r$ in tracing star formation activity, and we find that W3$-$W2 has a tighter anti-correlation with log $f_{\rm H_2}$ than the anti-correlation of (NUV$- r$) - $f_{\rm HI}$, (NUV$- r$) - $f_{\rm H_2}$ and (W3$-$W2) - $f_{\rm HI}$. This indicates that W3$-$W2 can trace the H$_2$ fraction in galaxies. For gas ratio $M_{\rm H_2}$/$M_{\rm HI}$, only in the intermediate-$M_{\rm *}$ galaxies it appears to depend on $M_{\rm *}$ and NUV$- r$. We find a tight correlation between the molecular gas mass $M_{\rm H_2}$ and 12\,$\mu$m (W3) luminosities ($L_{\rm 12\,\mu m}$), and the slope is close to unity (1.03 $\pm$ 0.06) for the SMT sample. This correlation may reflect that the cold gas and dust are well mixed on global galactic scale. Using the all-sky 12\,$\mu$m (W3) data available in WISE, this correlation can be used to estimate CO flux for molecular gas observations and can even predict H$_2$ mass for star-forming galaxies.

Because of the disappearance of its accretion disk since the time period around 2012 November--December, XSS J12270-4859 has recently been identified as, in addition to PSR J1023+0038, another transitional millisecond pulsar binary. We have carried out detailed analysis of the Fermi Large Area Telescope data for the binary. While both spectra before and after the disk-disappearance transition are well described by an exponentially cut-off power law, typical for pulsars' emission in the Fermi's 0.2-300 GeV, a factor of 2 flux decrease related to the transition is detected. A weak orbital modulation is seen, but only detectable in the after-transition data, same to that found at X-rays. In the long-term light curve of the source before the transition, a factor of 3 flux variations are seen. Comparing to the properties of J1023+0038, we disucss the implications from these results. We suggest that since the modulation is aligned with that at X-rays in orbital phase, it possibly arises due to the occultation of the gamma-ray emitting region by the companion. The origin of the variations in the long-term light curve is not clear, because the source field also contains unidentified radio or X-ray sources and their contamination can not be excluded. Multi-wavelength observations of the source field will help identify the origin of the variations by detecting any related flux changes from the in-field sources.

The Fermi source 2FGL~J0523.3$-$2530 has recently been identified as a candidate millisecond pulsar binary with an orbital period of 16.5 hrs. We have carried out detailed studies of the source's emission properties by analyzing data taken with the Fermi Large Area Telescope in the 0.2--300 GeV energy range. Long-term, yearly variability from the source has been found, with a factor of 4 flux variations in 1--300 GeV. From spectral analysis, we find an extra spectral component at 2--3 GeV that causes the source brightening. While no orbital modulations have been found from the Fermi data over the whole period of 2008--2014, orbital modulation in the source's $>$2 GeV emission is detected during the last 1.5 yrs of the Fermi observation. Our results support the millisecond pulsar binary nature of 2FGL~J0523.3$-$2530. Multi-wavelength observations of the source are warranted in order to find any correlated flux variations and thus help determine the origin of the long-term variability, which currently is not understood.

We report on the first dark-matter (DM) search results from PandaX-I, a low threshold dual-phase xenon experiment operating at the China Jinping Underground Laboratory. In the 37-kg liquid xenon target with 17.4 live-days of exposure, no DM particle candidate event was found. This result sets a stringent limit for low-mass DM particles and disfavors the interpretation of previously-reported positive experimental results. The minimum upper limit, $3.7\times10^{-44}$\,cm$^2$, for the spin-independent isoscalar DM-particle-nucleon scattering cross section is obtained at a DM-particle mass of 49\,GeV/c$^2$ at 90\% confidence level.

We report on our searches for debris disks around seven relatively nearby radio pulsars, which are isolated sources and were carefully selected as the targets on the basis of our deep $K_s$-band imaging survey. The $K_s$ images obtained with the 6.5\,m Baade Magellan Telescope at Las Campanas Observatory are analyzed together with the \textit{Spitzer}/IRAC images at 4.5 and 8.0~$\mu$m and the \textit{WISE} images at 3.4, 4.6, 12 and 22~$\mu$m. No infrared (IR) counterparts of these pulsars are found, with flux upper limits of $\sim \mu$Jy at near-infrared ($\lambda<10\ \mu$m) and $\sim$10--1000\,$\mu$Jy at mid-infrared wavelengths ($\lambda>10\ \mu$m). The results of this search are discussed in terms of the efficiency of converting the pulsar spin-down energy to thermal energy and X-ray heating of debris disks, with comparison made to the two magnetars 4U~0142+61 and 1E~2259+586 which are suggested to harbor a debris disk.

The High Energy cosmic-Radiation Detection (HERD) facility onboard China's Space Station is planned for operation starting around 2020 for about 10 years. It is designed as a next generation space facility focused on indirect dark matter search, precise cosmic ray spectrum and composition measurements up to the knee energy, and high energy gamma-ray monitoring and survey. The calorimeter plays an essential role in the main scientific objectives of HERD. A 3-D cubic calorimeter filled with high granularity crystals as active material is a very promising choice for the calorimeter. HERD is mainly composed of a 3-D calorimeter (CALO) surrounded by silicon trackers (TK) from all five sides except the bottom. CALO is made of 9261 cubes of LYSO crystals, corresponding to about 55 radiation lengths and 3 nuclear interaction lengths, respectively. Here the simulation results of the performance of CALO with GEANT4 and FLUKA are presented: 1) the total absorption CALO and its absorption depth for precise energy measurements (energy resolution: 1\% for electrons and gamma-rays beyond 100 GeV, 20\% for protons from 100 GeV to 1 PeV); 2) its granularity for particle identification (electron/proton separation power better than $10^{-5}$); 3) the homogenous geometry for detecting particles arriving from every unblocked direction for large effective geometrical factor ($>$3 ${\rm m}^{2}{\rm sr}$ for electron and diffuse gamma-rays, $>$2 $ {\rm m}^{2}{\rm sr}$ for cosmic ray nuclei); 4) expected observational results such as gamma-ray line spectrum from dark matter annihilation and spectrum measurement of various cosmic ray chemical components.

We report the results from our analysis of the \textit{Fermi} Large Area Telescope data of the \textit{Fermi} unassociated source 2FGL J1906.5$+$0720, which is a high-ranked candidate pulsar. In order to better study our target, we first update the ephemeris for PSR J1907$+$0602, which is used to help remove any possible contamination due to strong emission from this nearby pulsar. From our analysis, 2FGL J1906.5$+$0720 is confirmed to have a significant low energy cutoff at $\sim$ 1~GeV in its emission (14$\sigma$--18$\sigma$ significance), consistent with those seen in young pulsars. We search for pulsations but no spin frequency signals are found in a frequency range of 0.5--32 Hz. No single model can fully describe the source's overall \fermi\ $\gamma$-ray spectrum, and the reason for this is the excess emission detected at energies of $\geq$4 GeV. The high-energy component possibly indicates emission from a pulsar wind nebula, when considering 2FGL J1906.5$+$0720 as a young pulsar. We conclude that 2FGL J1906.5$+$0720 is likely a pulsar based on the emission properties we have obtained, and observations at other energies are needed in order to confirm its pulsar nature.

We report on the results from our search for the Wide-field Infrared Survey Explorer detection of the Galactic low-mass X-ray binaries. Among 187 binaries catalogued in Liu et al. (2007), we find 13 counterparts and two candidate counterparts. For the 13 counterparts, two (4U~0614+091 and GX~339$-$4) have already been confirmed by previous studies to have a jet and one (GRS~1915+105) to have a candidate circumbinary disk, from which the detected infrared emission arose. Having collected the broad-band optical and near-infrared data in literature and constructed flux density spectra for the other 10 binaries, we identify that three (A0620$-$00, XTE J1118+480, and GX 1+4) are candidate circumbinary disk systems, four (Cen X-4, 4U 1700+24, 3A 1954+319, and Cyg X-2) had thermal emission from their companion stars, and three (Sco X-1, Her X-1, and Swift J1753.5$-$0127) are peculiar systems with the origin of their infrared emission rather uncertain. We discuss the results and WISE counterparts' brightness distribution among the known LMXBs, and suggest that more than half of the LMXBs would have a jet, a circumbinary disk, or the both.

Different pieces of observational evidence suggest the existence of disks around isolated neutron stars. Such disks could be formed from supernova fallback when neutron stars are born in core-collapse supernova explosions. Efforts have been made to search for disks around different classes of pulsars, which include millisecond pulsars, young neutron star classes (magnetars, central compact objects, and X-ray dim isolated neutron stars), and regular radio pulsars. We review the main results from observations at wavelengths of from optical to sub-millimeter/millimeter.

Recent observations strongly suggest that the millisecond pulsar binary PSR J1023+0038 has developed an accretion disk since 2013 June. We present a multi-wavelength analysis of PSR J1023+0038, which reveals that 1) its gamma-rays suddenly brightened within a few days in June/July 2013 and has remained at a high gamma-ray state for several months; 2) both UV and X-ray fluxes have increased by roughly an order of magnitude, and 3) the spectral energy distribution has changed significantly after the gamma-ray sudden flux change. Time variabilities associated with UV and X-rays are on the order of 100-500 seconds and 50-100 seconds, respectively. Our model suggests that a newly formed accretion disk due to the sudden increase of the stellar wind could explain the changes of all these observed features. The increase of UV is emitted from the disk, and a new component in gamma-rays is produced by inverse Compton scattering between the new UV component and pulsar wind. The increase of X-rays results from the enhancement of injection pulsar wind energy into the intra-binary shock due to the increase of the stellar wind. We also predict that the radio pulses may be blocked by the evaporated winds from the disk and the pulsar is still powered by rotation.

We have carried out $JHK_s$ imaging of 12 white dwarf debris disk candidates from the WIRED SDSS DR7 catalog, aiming to confirm or rule out disks among these sources. On the basis of positional identification and the flux density spectra, we find that seven white dwarfs have excess infrared emission, but mostly at WISE W1 and W2 bands, four are due to nearby red objects consistent with background galaxies or very low mass dwarfs, and one exhibits excess emission at $JHK_s$ consistent with an unresolved L0 companion at the correct distance. While our photometry is not inconsistent with all seven excesses arising from disks, the stellar properties are distinct from the known population of debris disk white dwarfs, making the possibility questionable. In order to further investigate the nature of these infrared sources, warm Spitzer imaging is needed, which may help resolve galaxies from the white dwarfs and provide more accurate flux measurements.

We report on the results from our $\gamma$-ray analysis of the supernova remnant (SNR) RCW 103 region. The data were taken with the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. An extended source is found at a position consistent with that of RCW 103, and its emission was only detected above 1 GeV (10$\sigma$ significance), having a power-law spectrum with a photon index of 2.0$\pm$0.1. We obtain its 1--300 GeV spectrum, and the total flux gives a luminosity of 8.3$\times 10^{33}$ erg s$^{-1}$ at a source distance of 3.3 kpc. Given the positional coincidence and property similarities of this source with other SNRs, we identify it as the likely Fermi $\gamma$-ray counterpart to RCW 103. Including radio measurements of RCW 103, the spectral energy distribution (SED) is modeled by considering emission mechanisms based on both hadronic and leptonic scenarios. We find that models in the two scenarios can reproduce the observed SED, while in the hadronic scenario the existence of SNR--molecular-cloud interaction is suggested as a high density of the target protons is required.

We present the results of our radio observational campaign on GRB 100418a, for which we used the Australia Telescope Compact Array (ATCA), Very Large Array (VLA) and the Very Long Baseline Array (VLBA). GRB 100418a was a peculiar GRB with unusual X-ray and optical afterglow profiles featuring a plateau phase with a very shallow rise. This observed plateau phase was believed to be due to a continued energy injection mechanism which powered the forward shock, giving rise to an unusual and long-lasting afterglow. The radio afterglow of GRB 100418a was detectable several weeks after the prompt emission. We conducted long-term monitoring observations of the afterglow and attempted to test the energy injection model advocating that the continuous energy injection is due to shells of material moving at a wide range of Lorentz factors. We obtained an upper limit of $\gamma$ $<$ 7 for the expansion rate of the GRB 100418a radio afterglow, indicating that the range-of-Lorentz factor model could only be applicable for relatively slow moving ejecta. A preferred explanation could be that continued activity of the central engine may have powered the long-lasting afterglow.

The origin of dust grains in the interstellar medium is still open problem. \cite{Nicholls2013} found the presence of a significant amount of dust around V1309 Sco which maybe originate from the merger of a contact binary. We investigate the origin of dust around V1309 Sco, and suggest that these dust grains are efficiently produced in the binary-merger ejecta. By means of \emph{AGBDUST} code, we estimate that $\sim 5.2\times10^{-4} M_\odot$ of dust grains are produced, and their radii are $\sim 10^{-5}$ cm. These dust grains mainly are composed of silicate and iron grains. Because the mass of the binary-merger ejecta is very small, the contribution of dust produced by binary-merger ejecta to the overall dust production in the interstellar medium is negligible. However, it is the most important that the discovery of a significant amount of dust around V1309 Sco offers a direct support for the idea---common-envelope ejecta provides an ideal environment for dust formation and growth. Therefore, we confirm that common-envelope ejecta can be important source of cosmic dust.

The electron gas inside a neutron star is highly degenerate and relativistic. Due to the electron-electron magnetic interaction, the differential susceptibility can equal or exceed 1, which causes the magnetic system of the neutron star to become metastable or unstable. The Fermi liquid of nucleons under the crust can be in a metastable state, while the crust is unstable to the formation of layers of alternating magnetization. The change of the magnetic stress acting on adjacent domains can result in a series of shifts or fractures in the crust. The releasing of magnetic free energy and elastic energy in the crust can cause the bursts observed in magnetars. Simultaneously, a series of shifts or fractures in the deep crust which is closed to the Fermi liquid of nucleons can trigger the phase transition of the Fermi liquid of nucleons from a metastable state to a stable state. The released magnetic free energy in the Fermi liquid of nucleons corresponds to the giant flares observed in some magnetars.

In order to constrain the broad-band spectral energy distribution of the afterglow of GRB 100621A, dedicated observations were performed in the optical/near-infrared with the 7-channel "Gamma-Ray Burst Optical and Near-infrared Detector" (GROND) at the 2.2m MPG/ESO telescope, in the sub-millimeter band with the large bolometer array LABOCA at APEX, and at radio frequencies with ATCA. Utilizing also Swift X-ray observations, we attempt an interpretation of the observational data within the fireball scenario. The afterglow of GRB 100621A shows a very complex temporal as well as spectral evolution. We identify three different emission components, the most spectacular one causing a sudden intensity jump about one hour after the prompt emission. The spectrum of this component is much steeper than the canonical afterglow. We interpret this component using the prescription of Vlasis et al. (2011) for a two-shell collision after the first shell has been decelerated by the circumburst medium. We use the fireball scenario to derive constraints on the microphysical parameters of the first shell. Long-term energy injection into a narrow jet seems to provide an adequate description. Another noteworthy result is the large ($A_V$ = 3.6 mag) line-of-sight host extinction of the afterglow in an otherwise extremely blue host galaxy.

We report on the discovery of an infrared cometary nebula around PSR J1549$-$4848 in our Spitzer survey of a few middle-aged radio pulsars. Following the discovery, multi-wavelength imaging and spectroscopic observations of the nebula were carried out. We detected the nebula in Spitzer IRAC 8.0, MIPS 24 and 70 $\mu$m imaging and in Spitzer IRS 7.5--14.4 $\mu$m spectroscopic observations, and also in the WISE all-sky survey at 12 and 22 $\mu$m.These data were analyzed in detail, and we find that the nebula can be described with a standard bow-shock shape, and that its spectrum contains polycyclic aromatic hydrocarbon and H$_2$ emission features. However, it is not certain which object drives the nebula. We analyze the field stars and conclude that none of them can be the associated object because stars with a strong wind or mass ejection that usually produce bow shocks are much brighter than the field stars. The pulsar is approximately 15\arcsec\ away from the region in which the associated object is expected to be located. In order to resolve the discrepancy, we suggest that a highly collimated wind could be emitted from the pulsar and produce the bow shock. X-ray imaging to detect the interaction of the wind with the ambient medium and high-spatial resolution radio imaging to determine the proper motion of the pulsar should be carried out, which will help verify the association of the pulsar with the bow shock nebula.

We report our search for \gamma-ray emission in the energy range from 100 MeV to 300 GeV from four Accreting Millisecond Pulsars (AMPs), SAX J1808.4-3658, IGR J00291+5934, XTE J1814-338, and XTE J0929-314. The data are from four-year observations carried out by Large Area Telescope (LAT) onboard the Fermi \gamma-ray Space Telescope. The AMPs were not detected, and their \gamma-ray luminosity upper limits we obtain are 5.1*10^33 ergs/s for SAX J1808.4-3658, 2.1*10^33 ergs/s for IGR J00291+5934, 1.2*10^34 ergs/s for XTE J1814-338, and 2.2*10^33 ergs/s for XTE J0929-314. We compare our results with \gamma-ray irradiation luminosities required for producing optical modulations seen from the companions in the AMPs, which has been suggested by Takata et al. (2012), and our upper limits have excluded \gamma-ray emission as the heating source in these systems except XTE J0929-314, the upper limit of which is not deep enough. Our results also do not support the model proposed by Takata et al. (2012) that relatively strong \gamma-ray emission could arise from the outer gap of a high-mass neutron star controlled by the photon-photon pair-creation for the AMPs. Two AMPs, SAX J1808.4-3658 and IGR J00291+5934, have the measurements of their spin-down rates, and we derive the upper limits of their \gamma-ray conversion efficiencies, which are 57% and 3%, respectively. We discuss the implications to the AMP systems by comparing the efficiency upper limit values with that of 20 \gamma-ray millisecond pulsars (MSP) detected by Fermi and the newly discovered transitional MSP binary J1023+0038.

The dynamics of collisionless galaxy can be described by the Vlasov-Poisson system. By the Jean's theorem, all the spherically symmetric steady galaxy models are given by a distribution of {\Phi}(E,L), where E is the particle energy and L the angular momentum. In a celebrated Doremus-Feix-Baumann Theorem, the galaxy model {\Phi}(E,L) is stable if the distribution {\Phi} is monotonically decreasing with respect to the particle energy E. On the other hand, the stability of {\Phi}(E,L) remains largely open otherwise. Based on a recent abstract instability criterion of Guo-Lin, we constuct examples of unstable galaxy models of f(E,L) and f(E) in which f fails to be monotone in E.

Properties of X-ray luminosities in low-mass X-ray binaries (LMXBs) mainly depend on donors. We have carried out a detailed study of donors in persistent neutron-star LMXBs (PLMXBs) by means of a population synthesis code. PLMXBs with different donors have different formation channels. Our numerical simulations show that more than 90% of PLMXBs have main sequence (MS) donors, and PLMXBs with red giant (RG) donors via stellar wind (Wind) are negligible. In our model, most of neutron stars (NSs) in PLMXBs with hydrogen-rich donors form via core-collapse supernovae, while more than 90% of NSs in PLMXBs with naked helium star (He) donors or white dwarf (WD) donors form via an evolution-induced collapse via helium star ($1.4 \leq M_{\rm He}/M_\odot \leq 2.5$) or an accretion-induced collapses for an accreting ONeMg WD.

Strange quark stars (SSs) may originate from accreting neutron stars (NSs) in low-mass X-ray binaries (LMXBs). Assuming that conversion of NS matter to SSs occurs when the core density of accreting NS reaches to the density of quark deconfinement, $\sim 5 \rho_0$, where $\rho_0\sim 2.7\times 10^{14}$g cm$^{-3}$ is nuclear saturation density, we investigate LMXBs with SSs (qLMXBs). In our simulations, about 1\permil --- 10% of LMXBs can produce SSs, which greatly depends on the masses of nascent NSs and the fraction of transferred matter accreted by the NSs. If the conversion does not affect binaries systems, LMXBs evolve into qLMXBs. We find that some observational properties (spin periods, X-ray luminosities and orbital periods) of qLMXBs are similar with those of LMXBs, and it is difficult to differ them. If the conversion disturbs the binaries systems, LMXBs can produce isolated SSs. These isolated SSs could be submillisecond pulsars, and their birthrate in the Galaxy is $\sim$5--70 per Myr.

In the WISE all-sky source catalogue there are 76 million mid-infrared (MIR) point sources that were detected at the first three WISE bands and have association with only one 2MASS near-IR source within 3 arcsec. We search for their identifications in the SIMBAD database and find 3.2 million identified sources. Based on these known sources, we establish three criteria for selecting candidate AGB stars in the Galaxy, which are three defined occupation zones in a color-color diagram, Galactic latitude |gb|< 20 deg, and "corrected" WISE third-band W3c < 11. Applying these criteria to the WISE+2MASS sources, 1.37 million of them are selected. We analyze the WISE third-band W3 distribution of the selected sources, and further establish that W3 < 8 is required in order to exclude a large fraction of normal stars in them. We therefore find 0.47 million candidate AGB stars in our Galaxy from the WISE source catalogue. Using W3c, we estimate their distances and derive their Galactic distributions. The candidates are generally located around the Galactic center uniformly, with 68% (1-sigma) of them within approximately 8 kpc. We discuss that optical spectroscopy can be used to verify the C-rich AGB stars in our candidates, and they will be good targets for the LAMOST survey that is planned to start from fall of 2012.

We report our multi-band infrared (IR) imaging of the transitional millisecond pulsar system J1023+0038, a rare pulsar binary known to have an accretion disk in 2000--2001. The observations were carried out with ground-based and space telescopes from near-IR to far-IR wavelengths. We detected the source in near-IR JH bands and Spitzer 3.6 and 4.5 $\mu$m mid-IR channels. Combined with the previously-reported optical spectrum of the source, the IR emission is found to arise from the companion star, with no excess emission detected in the wavelength range. Because our near-IR fluxes are nearly equal to those obtained by the 2MASS all-sky survey in 2000 Feb., the result indicates that the binary did not contain the accretion disk at the time, whose existence would have raised the near-IR fluxes to 2-times larger values. Our observations have thus established the short-term nature of the interacting phase seen in 2000--2001: the accretion disk at most existed for 2.5 yrs. The binary was not detected by the WISE all-sky survey carried out in 2010 at its 12 and 22 $\mu$m bands and our Herschel far-IR imaging at 70 and 160 $\mu$m. Depending on the assumed properties of the dust, the resulting flux upper limits provide a constraint of <3x10^{22}--3x10^{25} g on the mass of the dust grains that possibly exist as the remnant of the previously-seen accretion disk.

We report on time-resolved optical imaging of the X-ray binary SAX J1808.4-3658 during its quiescent state and 2008 outburst. The binary, containing an accretion-powered millisecond pulsar, has a large sinusoidal-like modulation in its quiescent optical emission. We employ a Markov chain Monte Carlo technique to fit our multi-band light curve data in quiescence with an irradiated star model, and derive a tight constraint of $50^{+6}_{-5}$ deg on the inclination angle $i$ of the binary system. The pulsar and its companion are constrained to have masses of $0.97^{+0.31}_{-0.22} M_{\sun}$ and $0.04^{+0.02}_{-0.01} M_{\sun}$ (both 1$\sigma$ ranges), respectively. The dependence of these results on the measurements of the companion's projected radial velocity is discussed. We also find that the accretion disk had nearly constant optical fluxes over a $\sim$500 day period in the quiescent state our data covered, but started brightening 1.5 months before the 2008 outburst. Variations in modulation during the outburst were detected in our four observations made 7-12 days after the start of the outburst, and a sinusoidal-like modulation with 0.2 mag amplitude changed to have a smaller amplitude of 0.1 mag. The modulation variations are discussed. We estimate the albedo of the companion during its quiescence and the outburst, which was approximately 0 and 0.8 (for isotropic emission), respectively. This large difference probably provides additional evidence that the neutron star in the binary turns on as a radio pulsar in quiescence.

The quasi-periodic oscillations (QPOs) in black hole (BH) systems of different scales are interpreted based on the magnetic reconnection of the large-scale magnetic fields generated by the toroidal electric currents flowing in the inner region of accretion disk, where the current density is assumed to be proportional to the mass density of the accreting plasma. The magnetic connection (MC) is taken into account in resolving the dynamic equations of the accretion disk, in which the MC between the inner and outer disk regions, the MC between the plunging region and the disk, and the MC between the BH horizon and the disk are involved. It turns out that the single QPO frequency of several BH systems of different scales can be fitted by invoking the magnetic reconnection due to the MC between the inner and outer regions of the disk, where the BH binaries XTE J1859+226, XTE J1650-500 and GRS 1915+105 and the massive BHs in NGC 5408 X-1 and RE J1034+396 are included. In addition, the X-ray spectra corresponding to the QPOs are fitted for these sources based on the typical disk-corona model.

Solving Newtonian steady-state wind equations with accurate weak interaction rates and magnetic fields (MFs) of young neutron stars considered, we study the dynamics and nucleosynthesis of neutrino-driven winds (NDWs) from proto neutron stars (PNSs). For a typical 1.4 M$_{\odot}$ PNS model, we find the nucleosynthesis products are closely related to the luminosity of neutrinos and anti-neutrinos. The lower the luminosity is, the larger effect to the NDWs caused by weak interactions and MFs is. At a high anti-neutrino luminosity of typically $8\times 10^{51}$ erg s$^{-1}$, neutrinos and anti-neutrinos dominate the processes in a NDW and the MFs hardly change the wind's properties. While at a low anti-neutrino luminosity of $10^{51}$ erg s$^{-1}$ at the late stage of a NDW, the mass of product and nucleosynthesis are changed significantly in the strong MFs, they are less important than those in the early stage when the anti-neutrino luminosity is high. Therefore for the most models considered for the NDWs from PNSs, based on our calculations the influences of MFs and the net weak interactions on the nucleosynthesis is not significant.

We present the results of a search for gravitational waves associated with 154 gamma-ray bursts (GRBs) that were detected by satellite-based gamma-ray experiments in 2009-2010, during the sixth LIGO science run and the second and third Virgo science runs. We perform two distinct searches: a modeled search for coalescences of either two neutron stars or a neutron star and black hole; and a search for generic, unmodeled gravitational-wave bursts. We find no evidence for gravitational-wave counterparts, either with any individual GRB in this sample or with the population as a whole. For all GRBs we place lower bounds on the distance to the progenitor, under the optimistic assumption of a gravitational-wave emission energy of 10^-2 M c^2 at 150 Hz, with a median limit of 17 Mpc. For short hard GRBs we place exclusion distances on binary neutron star and neutron star-black hole progenitors, using astrophysically motivated priors on the source parameters, with median values of 16 Mpc and 28 Mpc respectively. These distance limits, while significantly larger than for a search that is not aided by GRB satellite observations, are not large enough to expect a coincidence with a GRB. However, projecting these exclusions to the sensitivities of Advanced LIGO and Virgo, which should begin operation in 2015, we find that the detection of gravitational waves associated with GRBs will become quite possible.

A model of low-frequency quasi-periodic oscillations (LFQPOs) of black hole X-ray binaries (BHXBs) is proposed based on the perturbed magnetohydrodynamic (MHD) equations of accretion disk. It turns out that the LFQPOs frequencies of some BHXBs can be fitted by the frequencies of the toroidal Alfv\'en wave oscillation corresponding to the maximal radiation flux. In addition, the positive correlation of the LFQPO frequencies with the radiation flux from accretion disk is well interpreted.

We present the first imaging circular polarimetry of the anomalous X-ray pulsar (AXP) 4U 0142+61 at optical wavelengths. The AXP is the only magnetar that has been well studied at optical and infrared wavelengths and is known to have a complicated broad-band spectrum over the wavelength range. The optical polarimetric observation was carried out with the 8.2-m Subaru telescope at I-band. From the observation, the degree of circular polarization V was measured to be V=1.1+/-2.0%, or |V|<4.3% (90% confidence). The relatively large uncertainty was due to the faintness of the source (I=23.4-24.0). Considering the currently suggested models for optical emission from magnetars, our result is not sufficiently conclusive to discriminate the models. We suggest that because linear polarization is expected to be strong in the models, linear polarimetry of this magnetar should be conducted.