In an effort to probe the origin of surface brightness profile (SBP) breaks widely observed in nearby disk galaxies, we carry out a comparative study of stellar population profiles of 635 disk galaxies selected from the MaNGA spectroscopic survey. We classify our galaxies into single exponential (TI), down-bending (TII) and up-bending (TIII) SBP types, and derive their spin parameters and radial profiles of age/metallicity-sensitive spectral features. Most TII (TIII) galaxies have down-bending (up-bending) star formation rate (SFR) radial profiles, implying that abrupt radial changes of SFR intensities contribute to the formation of both TII and TIII breaks. Nevertheless, a comparison between our galaxies and simulations suggests that stellar migration plays a significant role in weakening down-bending $\Sigma_{\star}$ profile breaks. While there is a correlation between the break strengths of SBPs and age/metallicity-sensitive spectral features for TII galaxies, no such correlation is found for TIII galaxies, indicating that stellar migration may not play a major role in shaping TIII breaks, as is evidenced by a good correspondence between break strengths of $\Sigma_{\star}$ and surface brightness profiles of TIII galaxies. We do not find evidence for galaxy spin being a relevant parameter for forming different SBP types, nor do we find significant differences between the asymmetries of galaxies with different SBP types, suggesting that environmental disturbances or satellite accretion in the recent past do not significantly influence the break formation. By dividing our sample into early and late morphological types, we find that galaxies with different SBP types follow nearly the same tight stellar mass-$R_{25}$ relation, which makes the hypothesis that stellar migration alone can transform SBP types from TII to TI and then to TIII highly unlikely.

The "ATOMS," standing for {\it ALMA Three-millimeter Observations of Massive Star-forming regions}, survey has observed 146 active star forming regions with ALMA Band 3, aiming to systematically investigate the spatial distribution of various dense gas tracers in a large sample of Galactic massive clumps, to study the roles of stellar feedback in star formation, and to characterize filamentary structures inside massive clumps. In this work, the observations, data analysis, and example science of the "ATOMS" survey are presented, using a case study for the G9.62+0.19 complex. Toward this source, some transitions, commonly assumed to trace dense gas, including CS $J = 2-1$, HCO$^+$ $J = 1-0$ and HCN $J = 1-0$, are found to show extended gas emission in low density regions within the clump; less than 25\% of their emission is from dense cores. SO, CH$_3$OH, H$^{13}$CN and HC$_3$N show similar morphologies in their spatial distributions and reveal well the dense cores. Widespread narrow SiO emission is present (over $\sim$1 pc), which may be caused by slow shocks from large--scale colliding flows or H{\sc ii} regions. Stellar feedback from an expanding H{\sc ii} region has greatly reshaped the natal clump, significantly changed the spatial distribution of gas, and may also account for the sequential high-mass star formation in the G9.62+0.19 complex. The ATOMS survey data can be jointly analyzed with other survey data, e.g., "MALT90", "Orion B", "EMPIRE", "ALMA\_IMF", and "ALMAGAL", to deepen our understandings of "dense gas" star formation scaling relations and massive proto-cluster formation.

We report studies of the relationships between the total bolometric luminosity ($L_{\rm bol}$ or $L_{\rm TIR}$) and the molecular line luminosities of $J=1-0$ transitions of H$^{13}$CN, H$^{13}$CO$^+$, HCN, and HCO$^+$ with data obtained from ACA observations in the "ATOMS" survey of 146 active Galactic star forming regions. The correlations between $L_{\rm bol}$ and molecular line luminosities $L'_{\rm mol}$ of the four transitions all appear to be approximately linear. Line emission of isotopologues shows as large scatters in $L_{\rm bol}$-$L'_{\rm mol}$ relations as their main line emission. The log($L_{\rm bol}$/$L'_{\rm mol}$) for different molecular line tracers have similar distributions. The $L_{\rm bol}$-to-$L'_{\rm mol}$ ratios do not change with galactocentric distances ($R_{\rm GC}$) and clump masses ($M_{\rm clump}$). The molecular line luminosity ratios (HCN-to-HCO$^+$, H$^{13}$CN-to-H$^{13}$CO$^+$, HCN-to-H$^{13}$CN and HCO$^+$-to-H$^{13}$CO$^+$) all appear constant against $L_{\rm bol}$, dust temperature ($T_{\rm d}$), $M_{\rm clump}$ and $R_{\rm GC}$. Our studies suggest that both the main lines and isotopologue lines are good tracers of the total masses of dense gas in Galactic molecular clumps. The large optical depths of main lines do not affect the interpretation of the slopes in star formation relations. We find that the mean star formation efficiency (SFE) of massive Galactic clumps in the "ATOMS" survey is reasonably consistent with other measures of the SFE for dense gas, even those using very different tracers or examining very different spatial scales.

The iron line and the continuum-fitting methods are currently the two leading techniques for measuring black hole spins with electromagnetic radiation. They can be naturally extended for probing the spacetime geometry around black holes and testing general relativity in the strong field regime. In the past couple of years, there has been significant work to use the iron line method to test the nature of black holes. In this Letter, we use the continuum-fitting method and we show its capability of constraining the spacetime geometry around black holes by analyzing 17 RXTE data of the X-ray binary LMC X-1.

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

Light dark sectors in thermal contact with the Standard Model naturally produce the observed relic dark matter abundance and are the targets of a broad experimental search program. A key light dark sector model is the pseudo-Dirac fermion with a dark photon mediator. The dynamics of the fermionic excited states are often neglected. We consider scenarios in which a nontrivial abundance of excited states is produced and their subsequent de-excitation yields interesting electromagnetic signals in direct detection experiments. We study three mechanisms of populating the excited state: a primordial excited fraction, a component up-scattered in the sun, and a component up-scattered in the Earth. We find that the fractional abundance of primordial excited states is generically depleted to exponentially small fractions in the early universe. Nonetheless, this abundance can produce observable signals in current dark matter searches. MeV-scale dark matter with thermal cross sections and higher can be probed by down-scattering following excitation in the sun. Up-scatters of GeV-scale dark matter in the Earth can give rise to signals in current and upcoming terrestrial experiments and X-ray observations. We comment on the possible relevance of these scenarios to the recent excess in XENON1T.

The youngest, closest and most compact embedded massive star cluster known excites the supernebula in the nearby dwarf galaxy NGC 5253. It is a crucial target and test case for studying the birth and evolution of the most massive star clusters. We present observations of the ionized gas in this source with high spatial and spectral resolution. The data includes continuum images of free-free emission with ~0.15'' resolution made with the JVLA at 15, 22 and 33 GHz, and a full data cube of the [SIV]10.5 micron fine-structure emission line with ~4.5 km/s velocity resolution and 0.3'' beam, obtained with TEXES on Gemini North. We find that 1) the ionized gas extends out from the cluster in arms or jets, and 2) the ionized gas comprises two components offset both spatially and in velocity. We discuss mechanisms that may have created the observed velocity field; possibilities include large-scale jets or a subcluster falling onto the main source.

After almost 20 years of hunting, only about a dozen hot corinos, hot regions enriched in interstellar complex organic molecules (iCOMs), are known. Of them, many are binary systems with the two components showing drastically different molecular spectra. Two obvious questions arise. Why are hot corinos so difficult to find and why do their binary components seem chemically different? The answer to both questions could be a high dust opacity that would hide the molecular lines. To test this hypothesis, we observed methanol lines at centimeter wavelengths, where dust opacity is negligible, using the Very Large Array interferometer. We targeted the NGC 1333 IRAS 4A binary system, for which one of the two components, 4A1, has a spectrum deprived of iCOMs lines when observed at millimeter wavelengths, while the other component, 4A2, is very rich in iCOMs. We found that centimeter methanol lines are similarly bright toward 4A1 and 4A2. Their non-LTE analysis indicates gas density and temperature ($\geq2\times10^6$ cm$^{-3}$ and 100--190 K), methanol column density ($\sim10^{19}$ cm$^{-2}$) and extent ($\sim$35 au in radius) similar in 4A1 and 4A2, proving that both are hot corinos. Furthermore, the comparison with previous methanol line millimeter observations allows us to estimate the optical depth of the dust in front of 4A1 and 4A2, respectively. The obtained values explain the absence of iCOMs line emission toward 4A1 at millimeter wavelengths and indicate that the abundances toward 4A2 are underestimated by $\sim$30\%. Therefore, centimeter observations are crucial for the correct study of hot corinos, their census, and their molecular abundances.

We present 1.3 mm ALMA dust polarization observations at a resolution of $\sim$0.02 pc of three massive molecular clumps, MM1, MM4, and MM9, in the infrared dark cloud G28.34+0.06. With the sensitive and high-resolution continuum data, MM1 is resolved into a cluster of condensations. The magnetic field structure in each clump is revealed by the polarized emission. We found a trend of decreasing polarized emission fraction with increasing Stokes $I$ intensities in MM1 and MM4. Using the angular dispersion function method (a modified Davis-Chandrasekhar-Fermi method), the plane-of-sky magnetic field strength in two massive dense cores, MM1-Core1 and MM4-Core4, are estimated to be $\sim$1.6 mG and $\sim$0.32 mG, respectively. \textbf{The ordered magnetic energy is found to be smaller than the turbulent energy in the two cores, while the total magnetic energy is found to be comparable to the turbulent energy.} The total virial parameters in MM1-Core1 and MM4-Core4 are calculated to be $\sim$0.76 and $\sim$0.37, respectively, suggesting that massive star formation does not start in equilibrium. Using the polarization-intensity gradient-local gravity method, we found that the local gravity is closely aligned with intensity gradient in the three clumps, and the magnetic field tends to be aligned with the local gravity in MM1 and MM4 except for regions near the emission peak, which suggests that the gravity plays a dominant role in regulating the gas collapse. Half of the outflows in MM4 and MM9 are found to be aligned within 10$^{\circ}$ of the condensation-scale ($<$0.05 pc) magnetic field, indicating that the magnetic field could play an important role from condensation to disk scale in the early stage of massive star formation. We also found that the fragmentation in MM1-Core1 cannot be solely explained by thermal Jeans fragmentation or turbulent Jeans fragmentation.

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

We present observations of SO and $\rm SO_2$ lines toward the shocked regions along the L1157 chemically rich outflow, taken in the context of the Seeds Of Life In Space IRAM-NOrthern Extended Millimeter Array Large Program, and supported by data from Submillimeter Array and IRAM-30 m telescope at 1.1--3.6 mm wavelengths. We simultaneously analyze, for the first time, all of the brightest shocks in the blueshifted lobe, namely, B0, B1, and B2. We found the following. (1) SO and $\rm SO_2$ may trace different gas, given that the large(-scale) velocity gradient analysis indicates for $\rm SO_2$ a volume density ($\rm 10^5\text{--}10^6\,cm^{-3}$) denser than that of the gas emitting in SO by a factor up to an order of magnitude. (2) Investigating the 0.1 pc scale field of view, we note a tentative gradient along the path of the precessing jet. More specifically, $\rm \chi({SO/SO_2})$ decreases from the B0-B1 shocks to the older B2. (3) At a linear resolution of 500--1400 au, a tentative spatial displacement between the two emitting molecules is detected, with the SO peak closer (with respect to $\rm SO_2$) to the position where the recent jet is impinging on the B1 cavity wall. Our astrochemical modeling shows that the SO and $\rm SO_2$ abundances evolve on timescales less than about 1000 years. Furthermore, the modeling requires high abundances ($2\times10^{-6}$) of both $\rm H_2S/H$ and S/H injected in the gas phase due to the shock occurrence, so pre-frozen OCS only is not enough to reproduce our new observations.

We propose a new machine learning approach to Stokes inversion based on a convolutional neural network (CNN) and the Milne-Eddington (ME) method. The Stokes measurements used in this study were taken by the Near InfraRed Imaging Spectropolarimeter (NIRIS) on the 1.6 m Goode Solar Telescope (GST) at the Big Bear Solar Observatory. By learning the latent patterns in the training data prepared by the physics-based ME tool, the proposed CNN method is able to infer vector magnetic fields from the Stokes profiles of GST/NIRIS. Experimental results show that our CNN method produces smoother and cleaner magnetic maps than the widely used ME method. Furthermore, the CNN method is 4~6 times faster than the ME method, and is able to produce vector magnetic fields in near real-time, which is essential to space weather forecasting. Specifically, it takes ~50 seconds for the CNN method to process an image of 720 x 720 pixels comprising Stokes profiles of GST/NIRIS. Finally, the CNN-inferred results are highly correlated to the ME-calculated results and are closer to the ME's results with the Pearson product-moment correlation coefficient (PPMCC) being closer to 1 on average than those from other machine learning algorithms such as multiple support vector regression and multilayer perceptrons (MLP). In particular, the CNN method outperforms the current best machine learning method (MLP) by 2.6% on average in PPMCC according to our experimental study. Thus, the proposed physics-assisted deep learning-based CNN tool can be considered as an alternative, efficient method for Stokes inversion for high resolution polarimetric observations obtained by GST/NIRIS.

We examine the statistical properties of polarization maps from Planck 2018 within the patch of sky observed by the BICEP2/Keck experiment using the one point distribution function (1PDF), skewness, and kurtosis statistics. Our analysis is performed for the Q and U Stokes parameters and for the corresponding E- and B-modes of the CMB signals. We extend our analysis by studying the correlations between CMB polarization maps and residual maps (the difference between the full signal and the CMB map) for the frequency range of 100--217 GHz with both the Q/U and E/B approaches. Although all the CMB maps reveal almost Gaussian statistical properties for Q/U and E/B domains, we have detected very significant anomalies for cross-correlations with residuals at 100 GHz at the level of $3.7\sigma$ for the Commander map and $5.2\sigma$ for NILC, for both the Q and U parameters. Using the NILC--Commander difference, which does not contain a cosmological signal, we find a sub-dominant non-Gaussian component in Q skewness and kurtosis at the level of $4.3\sigma$ and $10\sigma$, respectively. For the B-mode we have found a very high level of cross-correlation (0.63--0.69) between the NILC/Commander maps and the 143 GHz total signal, which cannot be associated with the cosmological component. These strong deviations suggest that remnants of foregrounds, systematic effects, and component separation exist in the 2018 Planck CMB polarization maps in the BICEP2 sky area, which is far away from the Galactic plane. Our analysis also demonstrates the preferability of the Q/U domain over E/B for determination of the statistical properties of the derived CMB signals, due to non-locality of the transition Q/U $\rightarrow$ E/B.

GW Ori is a hierarchical triple system which has a rare circumtriple disk. We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of 1.3 mm dust continuum and 12CO J=2-1 molecular gas emission of the disk. For the first time, we identify three dust rings in the disk at ~46, 188, and 338 AU, with estimated dust mass of ~70-250 Earth masses, respectively. To our knowledge, the outer ring in GW Ori is the largest dust ring ever found in protoplanetary disks. We use visibility modelling of dust continuum to show that the disk has misaligned parts and the innermost dust ring is eccentric. The disk misalignment is also suggested by the CO kinematics modelling. We interpret these substructures as evidence of ongoing dynamical interactions between the triple stars and the circumtriple disk.

We report the discovery of a highly dispersed fast radio burst, FRB~181123, from an analysis of $\sim$1500~hr of drift-scan survey data taken using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The pulse has three distinct emission components, which vary with frequency across our 1.0--1.5~GHz observing band. We measure the peak flux density to be $>0.065$~Jy and the corresponding fluence $>0.2$~Jy~ms. Based on the observed dispersion measure of 1812~cm$^{-3}$~pc, we infer a redshift of $\sim 1.9$. From this, we estimate the peak luminosity and isotropic energy to be $\lesssim 2\times10^{43}$~erg~s$^{-1}$ and $\lesssim 2\times10^{40}$~erg, respectively. With only one FRB from the survey detected so far, our constraints on the event rate are limited. We derive a 95\% confidence lower limit for the event rate of 900 FRBs per day for FRBs with fluences $>0.025$~Jy~ms. We performed follow-up observations of the source with FAST for four hours and have not found a repeated burst. We discuss the implications of this discovery for our understanding of the physical mechanisms of FRBs.

We analyze multi-frequency, high-resolution continuum data obtained by ALMA and JVLA to study detailed structure of the dust distribution in the infant disk of a Class~0/I source, L1527 IRS. We find three clumps aligning in the north-south direction in the $7 {\rm \,mm}$ radio continuum image. The three clumps remain even after subtracting free-free contamination, which is estimated from the $1.3{\rm \,cm}$ continuum observations. The northern and southern clumps are located at a distance of $\sim 15{\rm \,au}$ from the central clump and are likely optically thick at $7{\rm \,mm}$ wavelength. The clumps have similar integrated intensities. The symmetric physical properties could be realized when a dust ring or spiral arms around the central protostar is projected to the plane of the sky. We demonstrates for the first time that such substructure may form even in the disk-forming stage, where the surrounding materials actively accrete toward a disk-protostar system.

Low-surface-brightness galaxies(LSBGs) are defined as galaxies that are fainter than dark night sky and are important for studying our universe. Particularly, edge-on galaxies are useful for the study of rotational velocity and dynamical properties of galaxies. Hence here we focus on searching for edge-on LSBGs. In order to find these edge-on dim galaxies, a series of effects caused by inclination, including the surface brightness profile, internal extinction, and scale length, have been corrected. In this work, we present a catalog of 281 edge-on LSBG candidates, which are selected from the cross-match between SDSS DR7 and the 40% ALFALFA catalog. We also present the properties of these edge-on LSBG candidates including absolute magnitude, central surface brightness, B-V color, scale length, and relative thickness. Our result suggests that the correction of inclination effects is very important for obtaining a complete sample of LSBGs.

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

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

A dark photon may kinetically mix with the Standard Model photon, leading to observable cosmological signatures. The mixing is resonantly enhanced when the dark photon mass matches the primordial plasma frequency, which depends sensitively on the underlying spatial distribution of electrons. Crucially, inhomogeneities in this distribution can have a significant impact on the nature of resonant conversions. We develop and describe, for the first time, a general analytic formalism to treat resonant oscillations in the presence of inhomogeneities. Our formalism follows from the theory of level crossings of random fields and only requires knowledge of the one-point probability distribution function (PDF) of the underlying electron number density fluctuations. We validate our formalism using simulations and illustrate the photon-to-dark photon conversion probability for several different choices of PDFs that are used to characterize the low-redshift Universe.

X-ray reflection spectroscopy is potentially a powerful tool to probe the spacetime geometry around astrophysical black holes and test general relativity in the strong field regime. However, precision tests of general relativity are only possible if we employ the correct astrophysical model and we can limit the systematic uncertainties. It is thus crucial to select the sources and the observations most suitable for these tests. In this work, we analyze simultaneous observations of XMM-Newton and NuSTAR of the supermassive black hole in Fairall 9. This source has a number of properties that make it a promising candidate for tests of general relativity using X-ray reflection spectroscopy. Nevertheless, we find that with the available data there is not a unique interpretation of the spectrum of Fairall 9, which prevents, for the moment, to use this source for robust tests of general relativity. This issue may be solved by future X-ray missions with a higher energy resolution near the iron line.

We compile a sample of about 157,000 spaxels from the Mapping Nearby Galaxies at the Apache Point Observatory survey to derive the average dust attenuation curve for subgalactic star-forming regions of local star-forming galaxies (SFGs) in the optical wavelength, following the method of \cite{Calzetti1994}. We obtain a $D_n(4000)$-independent average attenuation curve for spaxels with $1.1\leq D_n(4000)<1.3$, which is similar to the one derived from either local starbursts or normal SFGs. We examine whether and how the shape of the average attenuation curve changes with several local and global physical properties. For spaxels with $1.2\leq D_n(4000)<1.3$, we find no dependence on either local or global physical properties for the shape of the average attenuation curve. However, for spaxels with younger stellar population ($1.1\leq D_n(4000)<1.2$), shallower average attenuation curves are found for star-forming regions with smaller stellar mass surface density, smaller star formation rate surface density, or those residing in the outer region of galaxies. These results emphasize the risk of using one single attenuation curve to correct the dust reddening for all types of star-forming regions, especially for those with fairly young stellar population.

We report the results of a pilot program for a Green Bank Telescope (GBT) MUSTANG Galactic Plane survey at 3 mm (90 GHz), MGPS90. The survey achieves a typical $1\sigma$ depth of $1-2$ mJy beam$^{-1}$ with a 9" beam. We describe the survey parameters, quality assessment process, cataloging, and comparison with other data sets. We have identified 709 sources over seven observed fields selecting some of the most prominent millimeter-bright regions between $0\deg < \ell < 50\deg$ (total area $\approx 7.5 \deg^2$). The majority of these sources have counterparts at other wavelengths. By applying flux selection criteria to these sources, we successfully recovered several known hypercompact HII (HCHII) regions, but did not confirm any new ones. We identify 126 sources that have mm-wavelength counterparts but do not have cm-wavelength counterparts and are therefore candidate HCHII regions; of these, 10 are morphologically compact and are strong candidates for new HCHII regions. Given the limited number of candidates in the extended area in this survey compared to the relatively large numbers seen in protoclusters W51 and W49, it appears that most HCHII regions exist within dense protoclusters. Comparing the counts of HCHII to ultracompact HII (UCHII) regions, we infer the HCHII region lifetime is 16-46% that of the UCHII region lifetime. We additionally separated the 3 mm emission into dust and free-free emission by comparing with archival 870 $\mu$m and 20 cm data. In the selected pilot fields, most ($\gtrsim80$%) of the 3 mm emission comes from plasma, either through free-free or synchrotron emission.

A dark photon may kinetically mix with the ordinary photon, inducing oscillations with observable imprints on cosmology. Oscillations are resonantly enhanced if the dark photon mass equals the ordinary photon plasma mass, which tracks the free electron number density. Previous studies have assumed a homogeneous Universe; in this Letter, we introduce for the first time an analytic formalism for treating resonant oscillations in the presence of inhomogeneities of the photon plasma mass. We apply our formalism to determine constraints from Cosmic Microwave Background photons oscillating into dark photons, and from heating of the primordial plasma due to dark photon dark matter converting into low-energy photons. Including the effect of inhomogeneities demonstrates that prior homogeneous constraints are not conservative, and simultaneously extends current experimental limits into a vast new parameter space.

We report the discovery of a Tidal Disruption Event (TDE) candidate occurring in the Active Galactic Nucleus SDSS J022700.77-042020.6. A sudden increase in flux of J0227-0420 during the second half of 2009 is clearly shown in the long-term optical, UV, and NIR light curves. A plateau phase, following an initial decline, is seen in the NUV and optical u, g, r, i light curves. Moreover, we find possible evidence that the plateau phase in the NUV band may lag behind the optical ones by approximately 70-80 days with also a much shorter duration, i.e. $\sim$7-15 days against $\sim$40-50 days. The long-term NUV/optical (after the plateau phase), NIR and MIR light curves can be well characterized with a form of $L(t)\propto t^{-\beta}$, consistent with the expectation of a TDE. The plateaus can be explained if the stellar streams collide with the pre-existing AGN disk at different radii. Though the overall fallback rate decreases, the material in the outer disk gradually drifts inward and increases the local accretion rate at the inner region, producing the optical and UV plateaus. The possible lag between the optical and NUV plateaus can then be attributed to viscosity delay. The index $\beta$ of the NIR $J, H, K_s$ bands ($\sim1.4-3.3$) is steeper than that of the UV/optical ($\sim0.7-1.3$) and MIR bands ($\sim0.9-1.8$), which may suggest that a certain fraction of the dust in the inner region of the dusty torus may be sublimated during the TDE phase. Our results indicate that, due to collisions between stellar debris and pre-existing disk, the light curves of TDEs occurring in AGN may show distinctive features, which may shed new light on the accretion process.

Transition Radiation (TR) plays an important role in particle identification in high-energy physics and its characteristics provide a feasible method of energy calibration in the energy range up to 10 TeV, which is of interest for dark matter searches in cosmic rays. In a Transition Radiation Detector (TRD), the TR signal is superimposed onto the ionization energy loss signal induced by incident charged particles. In order to make the TR signal stand out from the background of ionization energy loss in a significant way, we optimized both the radiators and the detector. We have designed a new prototype of regular radiator optimized for a maximal TR photon yield, combined with the Side-On TRD which is supposed to improve the detection efficiency of TR. We started a test beam experiment with the Side-On TRD at Conseil Europ\'{e}en pour la Recherche Nucl\'{e}aire (CERN), and found that the experimental data is consistent with the simulation results.

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

We present Atacama Large Millimeter/submillimeter Array (ALMA) and Karl G. Jansky Very Large Array (JVLA) observations of the massive infrared dark cloud NGC 6334S (also known as IRDC G350.56+0.44), located at the southwestern end of the NGC 6334 molecular cloud complex. The H$^{13}$CO$^{+}$ and the NH$_{2}$D lines covered by the ALMA observations at a $\sim$3$^{\prime\prime}$ angular resolution ($\sim$0.02 pc) reveal that the spatially unresolved non-thermal motions are predominantly subsonic and transonic, a condition analogous to that found in low-mass star-forming molecular clouds. The observed supersonic non-thermal velocity dispersions in massive star forming regions, often reported in the literature, might be significantly biased by poor spatial resolutions that broaden the observed line widths due to unresolved motions within the telescope beam. Our 3~mm continuum image resolves 49 dense cores, whose masses range from 0.17 to 14 $M_{\odot}$. The majority of them are resolved with multiple velocity components. Our analyses of these gas velocity components find an anti-correlation between the gas mass and the virial parameter. This implies that the more massive structures tend to be more gravitationally unstable. Finally, we find that the external pressure in the NGC 6334S cloud is important in confining these dense structures, and may play a role in the formation of dense cores, and subsequently, the embedded young stars.

The role of ionization feedback on high-mass (>8 Msun) star formation (HMSF) is still highly debated. Questions remain concerning the presence of nearby HII regions changes the properties of early HMSF and whether HII regions promote or inhibit the formation of high-mass stars. To characterize the role of HII regions on the HMSF, we study the properties of a sample of candidates high-mass starless clumps (HMSCs), of which about 90% have masses larger than 100 Msun. These high-mass objects probably represent the earliest stages of HMSF; we search if (and how) their properties are modified by the presence of an HII region. We took advantage of the recently published catalog of HMSC candidates. By cross matching the HMSCs and HII regions, we classified HMSCs into three categories: 1) The HMSCs associated with HII regions both in the position in the projected plane of the sky and in velocity; 2) HMSCs associated in the plane of the sky, but not in velocity; and 3) HMSCs far away from any HII regions in the projected sky plane. We carried out comparisons between associated and nonassociated HMSCs based on statistical analyses of multiwavelength data from infrared to radio. Statistical analyses suggest that HMSCs associated with HII regions are warmer, more luminous, more centrally-peaked and turbulent. We also clearly show, for the first time, that the ratio of bolometric luminosity to envelope mass of HMSCs (L/M) could not be a reliable evolutionary probe for early HMSF due to the external heating effects of the HII regions. More centrally peaked and turbulent properties of HMSCs associated with HII regions may promote the formation of high-mass stars by limiting fragmentation. High resolution interferometric surveys toward HMSCs are crucial to reveal how HII regions impact the star formation process inside HMSCs.

Early shallow-decaying X-ray afterglows of gamma-ray bursts (GRBs) may be attributed to the dipole radiations of newly-born magnetars. Assuming that the GRB jets powered by magnetars are quasi-universal, we find that the jet structure can be parameterized as a uniform jet with a luminosity of $\log L_{\rm j}/{\rm erg\ s^{-1}}=52.68^{+0.76}_{-0.33}$ (1$\sigma$) and an opening angle $\theta_{\rm j}=2.10_{-1.28}^{+1.90}$ (50\% confidence level), surrounding by a power-law decay component with an index of ${-4.00^{+0.27}_{-0.37}}$ (1$\sigma$). The inferred local GRB rate is $\rho=9.6$ Gpc$^{-3}$ yr$^{-1}$ by including both the typical GRBs and LL-GRBs as the same population. The typical viewing angle is $3.3^{o}$, and may be $20^{o}\sim30^{o}$ for LL-GRBs. The X-ray luminosity function of the dipole radiation wind can be empirically described by a broken power-law function with indices $\beta_1=0.78^{+0.16}_{-0.15}$ and $\beta_2>1.6$ broken at $\log L_{b, w}/{\rm erg\ s^{-1}}=48.51^{+0.53}_{-0.65}$. In case of that the wind outflow is collimated and co-axial with the GRB jet, we find that the wind structure is similar to the GRB jet, i.e., $\log L_{\rm c, w}/{\rm erg\ s^{-1}}=48.38^{+0.30}_{-0.48}$, $\theta_{\rm c, w}={2.65^{o}}_{-1.19^{o}}^{+0.1.73^{o}}$, and $k_{\rm w}=4.57^{+1.21}_{-0.75}$. The observed correlation between the prompt gamma-ray luminosity and X-ray luminosity of the wind may be resulted from the viewing angle effect in such a jet-wind system. Discussion on survey with the X-ray instruments on board the {\em Einstein\ Probe} mission in the soft X-ray band for the jet and wind emission is also presented.

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

We present multiwavelength studies of a transient X-ray source, XMMSL1 J131952.3+225958, associated with the galaxy NGC 5092 at $z=0.023$ detected in the XMM-Newton SLew survey (XMMSL). The source brightened in the 0.2--2 keV band by a factor of $>20$ in 2005 as compared with previous flux limits and then faded by a factor of $>200$ as observed with it XMM-Newton in 2013 and with it Swift in 2018. At the flaring state, the X-ray spectrum can be modeled with a blackbody at a temperature of $\sim$ 60 eV and an overall luminosity of $\sim$ $1.5 \times 10^{43}$ erg s$^{-1}$. A UV flare and optical flare were also detected with the Galaxy Evolution Explorer and the Sloan Digital Sky Survey, respectively, within several months of the X-ray flare, whose nonstellar UV--optical spectrum can be described with a blackbody at a temperature of $\sim$ $(1-2) \times 10^4$ K and a luminosity of $\sim$ $(2-6) \times 10^{43}$ erg s$^{-1}$. Interestingly, mid-infrared monitoring observations of NGC 5092 with the Wide-field Infrared Survey Explorer 5--13 yr later show a continuous flux decline. These dramatic variability properties, from the X-ray through UV and optical to infrared, appear to be orderly, suggestive of a stellar tidal disruption event (TDE) by a massive black hole, confirming the postulation by Kanner et al.(2013). This TDE candidate belongs to a rare sample with contemporaneous bright emission detected in the X-ray, UV, and optical, which are later echoed by dust-reprocessed light in the mid-infrared. The black hole has a mass of $\sim$ $5 \times 10^{7} \rm M_{\odot}$, residing in a galaxy that is dominated by a middle-aged stellar population of 2.5 Gyr.

We study the properties of strange quark stars by employing a 3-flavor Nambu-Jona-Lasinio model with both scalar-isovector and vector-isovector interactions. Using the constraint on the vector-isoscalar interaction strength obtained from the elliptic flow splitting between particles and their antiparticles in relativistic heavy-ion collisions, we investigate the dependence of the properties of strange quark stars on the vector-isovector and the scalar-isovector interactions, and compare the results with the state-of-art astrophysical constraints on the compact star radius and mass as well as its tidal deformability from the GW170817 event. Results from our study reinforce the prospect of using both heavy-ion collisions and astrophysical observations to provide constraints on the isovector coupling strength in quark matter and thus the quark matter equation of state as well as the QCD phase structure at finite isospin chemical potentials.

We present two recurrent neural networks (RNNs), one based on gated recurrent units and the other based on long short-term memory, for predicting whether an active region (AR) that produces an M- or X-class flare will also produce a coronal mass ejection (CME). We model data samples in an AR as time series and use the RNNs to capture temporal information of the data samples. Each data sample has 18 physical parameters, or features, derived from photospheric vector magnetic field data taken by the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO). We survey M- and X-class flares that occurred from 2010 May to 2019 May using the Geostationary Operational Environmental Satellite's X-ray flare catalogs provided by the National Centers for Environmental Information (NCEI), and select those flares with identified ARs in the NCEI catalogs. In addition, we extract the associations of flares and CMEs from the Space Weather Database Of Notifications, Knowledge, Information (DONKI). We use the information gathered above to build the labels (positive versus negative) of the data samples at hand. Experimental results demonstrate the superiority of our RNNs over closely related machine learning methods in predicting the labels of the data samples. We also discuss an extension of our approach to predict a probabilistic estimate of how likely an M- or X-class flare will initiate a CME, with good performance results. To our knowledge this is the first time that RNNs have been used for CME prediction.

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

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

We study the origin of tail-like structures recently detected around the disk of SU Aurigae and several FU~Orionis-type stars. Dynamic protostellar disks featuring ejections of gaseous clumps and quiescent protoplanetary disks experiencing a close encounter with an intruder star were modeled using the numerical hydrodynamics code FEOSAD. Both the gas and dust dynamics were taken into account, including dust growth and mutual friction between the gas and dust components. Only plane-of-the-disk encounters were considered. Ejected clumps produce a unique type of tail that is characterized by a bow-shock shape. Such tails originate from the supersonic motion of ejected clumps through the dense envelope that often surrounds young gravitationally unstable protostellar disks. The ejected clumps either sit at the head of the tail-like structure or disperse if their mass is insufficient to withstand the head wind of the envelope. On the other hand, close encounters with quiescent protoplanetary disks produce three types of the tail-like structure; we define these as pre-collisional, post-collisional, and spiral tails. These tails can in principle be distinguished from one another by particular features of the gas and dust flow in and around them. We find that the brown-dwarf-mass intruders do not capture circumintruder disks during the encounter, while the subsolar-mass intruders can acquire appreciable circumintruder disks with elevated dust-to-gas ratios, which can ease their observational detection. However, this is true only for prograde collisions; the retrograde intruders fail to collect appreciable amounts of gas or dust from the disk of the target. The predicted mass of dust in the model tail-like structures is higher than what was inferred for similar structures in SU~Aur, FU~Ori, and Z~CMa, making their observational detection feasible. Abridged.

For filter-based magnetographs, the linear calibration method under the weak-field assumption is usually adopted; this leads to magnetic saturation effect in the regions with strong magnetic field. This article explores a new method to overcome the above disadvantage using a multilayer perceptron network, which we call MagMLP, based on a back-propagation algorithm with one input layer, five hidden layers, and one output layer. We use the data from the \textit{Spectropolarimeter} (SP) on board \textit{Hinode} to simulate single-wavelength observations for the model training, and take into account the influence of the Doppler velocity field and the filling factor. The training results show that the linear fitting coefficient (LFC) of the transverse field reaches above 0.91, and that of the longitudinal field is above 0.98. The generalization of the models is good because the corresponding LFCs are above 0.9 for the test subsets. Compared with the linear calibration method, the MagMLP is much more effective on dealing with the magnetic saturation effect. Analyzing an active region, the results of the linear calibration present an evident magnetic saturation effect in the umbra regions; the corresponding systematic error reaches values greater than 1000 G in most areas, or even exceeds 2000 G at some pixels. However, the results of MagMLP at these locations are very close to the inversion results, and the systematic errors are basically within 300 G. In addition, we find that there are many "bright spots" and "dark spots" on the inclination angle images from the inversion results of \textit{Hinode}/SP with values of 180 and 0 degrees, respectively, where the inversion is not reliable and does not produce a good result; the MagMLP handles these points well.

The Five-hundred-meter Aperture Spherical radio Telescope (FAST) passed national acceptance and is taking pilot cycle of 'Shared-Risk' observations. The 19-beam receiver covering 1.05-1.45 GHz was used for most of these observations. The electronics gain fluctuation of the system is better than 1\% over 3.5 hours, enabling enough stability for observations. Pointing accuracy, aperture efficiency and system temperature are three key parameters of FAST. The measured standard deviation of pointing accuracy is 7.9$''$, which satisfies the initial design of FAST. When zenith angle is less than 26.4$^\circ$, the aperture efficiency and system temperature around 1.4 GHz are $\sim$ 0.63 and less than 24 K for central beam, respectively. The measured value of these two parameters are better than designed value of 0.6 and 25 K, respectively. The sensitivity and stability of the 19-beam backend are confirmed to satisfy expectation by spectral HI observations toward N672 and polarization observations toward 3C286. The performance allows FAST to take sensitive observations in various scientific goals, from studies of pulsar to galaxy evolution.

We report the discovery of an extreme X-ray flux rise (by a factor of > 20) of the weak-line quasar SDSS J153913.47+395423.4 (hereafter SDSS J1539+3954) at z = 1.935. SDSS J1539+3954 is the most-luminous object among radio-quiet type 1 AGNs where such dramatic X-ray variability has been observed. Before the X-ray flux rise, SDSS J1539+3954 appeared X-ray weak compared with the expectation from its UV flux; after the rise, the ratio of its X-ray flux and UV flux is consistent with the majority of the AGN population. We also present a contemporaneous HET spectrum of SDSS J1539+3954, which demonstrates that its UV continuum level remains generally unchanged despite the dramatic increase in the X-ray flux, and its C iv emission line remains weak. The dramatic change only observed in the X-ray flux is consistent with a shielding model, where a thick inner accretion disk can block our line of sight to the central X-ray source. This thick inner accretion disk can also block the nuclear ionizing photons from reaching the high-ionization broad emission-line region, so that weak high-ionization emission lines are observed. Under this scenario, the extreme X-ray variability event may be caused by slight variations in the thickness of the disk. This event might also be explained by gravitational light-bending effects in a reflection model.

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

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

We study the distribution functions of the CMB polarization angle $\psi$, focusing on the Planck 2018 CMB maps. We extend the model of Preece & Battye (2014) of Gaussian correlated $Q$ and $U$ Stokes parameters to allow nonzero means. When the variances of $Q$ and $U$ are equal and their covariance and means are zero, the polarization angle is uniformly distributed. Otherwise the uniform distribution will be modulated by harmonics with $2\psi$ and $4\psi$ phases. These modulations are visible in the Planck 2018 CMB maps. Furthermore, the mean value of $U$ is peculiar compared to the power spectrum.

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.

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

The discovery of changing-look active galactic nuclei (CL AGNs), with appearance and disappearance of broad emission lines and/or with strong variation of line-of-sight column density within a few years, challenges the AGN unification model. We explore the physical mechanisms based on the X-ray spectral evolution for a sample of 15 CL AGNs. We find that the X-ray photon index, $\Gamma$, and Eddington-scaled X-ray luminosity, $L_{\rm 2-10 keV}/L_{\rm Edd}$, follow negative and positive correlations when $L_{\rm 2-10 keV}/L_{\rm Edd}$ is lower and higher than a critical value of $\sim 10^{-3}$. This different X-ray spectral evolution is roughly consistent with the prediction of the accretion-mode transition (e.g., clumpy cold gas or cold disk to advection dominated accretion flow, or vice visa). With quasi-simultaneous X-ray and optical spectrum observations within one year, we find that the CL AGNs observed with and without broad emission lines stay in the positive and negative part of the $\Gamma-L_{\rm 2-10 keV}/L_{\rm Edd}$ correlation respectively. Our result suggest that the change of the accretion mode may be the physical reason for the CL AGNs.

The anti- and positive correlations between X-ray photon index and Eddington-scaled X-ray luminosity were found in decay phase of X-ray binary outbursts and a sample of active galactic nuclei in former works. We further systematically investigate the evolution of X-ray spectral index, along the X-ray flux and Eddington ratio in eight outbursts of four black-hole X-ray binaries, where all selected outbursts have observational data from Rossi X-ray Timing Explorer in both rise and decay phases. In the initial rise phase, the X-ray spectral index is anti-correlated with the flux, and the X-ray spectrum quickly softens when the X-ray flux is approaching the peak value. In the decay phase, the X-ray photon index and the flux follow two different positive correlations, and they become anti-correlated again when the X-ray flux is below a critical value, where the anti-correlation part follow the same trend as that found in the initial rise phase. Compared with other X-ray binaries, GRO J1655-40 has an evident lower critical Eddington ratio for the anti- and positive transition, which suggests that its black-hole mass and distance are not well constrained or its intrinsic physic is different.

We report the discovery of large amplitude X-ray variability in the low luminosity AGN (LLAGN) MGC 7589, and present possible observational evidence for accretion mode transition in this source. Long-term X-ray flux variations by a factor of more than 50 are found using X-ray data obtained by Swift/XRT and XMM-Newton over 17 years. Results of long-term monitoring data in the UV, optical and infrared bands over ~20 years are also presented. The Eddington ratio increased from $10^{-3}$ to $\sim0.13$, suggesting a transition of the accretion flow from an ADAF to a standard thin accretion disc. Further evidence supporting the thin disc in the high luminosity state is found by the detection of a significant soft X-ray component in the X-ray spectrum. The temperature of this component ($\sim19^{+15}_{-7}$eV, fitted with a blackbody model) is in agreement with the predicted temperature of the inner region for a thin disc around a black hole (BH) with mass of $\sim10^{7}\,M_{\mathrm{Sun}}$. These results may indicate that NGC 7589 had experienced accretion mode transition over a timescale of a few years, suggesting the idea that similar accretion processes are at work for massive black hole and black hole X-ray binaries.

We present a pilot HI survey of 17 Planck Galactic Cold Clumps (PGCCs) with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). HI Narrow Self-Absorption (HINSA) is an effective method to detect cold HI being mixed with molecular hydrogen H$_2$ and improves our understanding of the atomic to molecular transition in the interstellar medium. HINSA was found in 58\% PGCCs that we observed. The column density of HINSA was found to have an intermediate correlation with that of $^{13}$CO, following $\rm log( N(HINSA)) = (0.52\pm 0.26) log(N_{^{13}CO}) + (10 \pm 4.1) $. HI abundance relative to total hydrogen [HI]/[H] has an average value of $4.4\times 10^{-3}$, which is about 2.8 times of the average value of previous HINSA surveys toward molecular clouds. For clouds with total column density N$\rm_H >5 \times 10^{20}$ cm$^{-2}$, an inverse correlation between HINSA abundance and total hydrogen column density is found, confirming the depletion of cold HI gas during molecular gas formation in more massive clouds. Nonthermal line width of $^{13}$CO is about 0-0.5 km s$^{-1}$ larger than that of HINSA. One possible explanation of narrower nonthermal width of HINSA is that HINSA region is smaller than that of $^{13}$CO. Based on an analytic model of H$_2$ formation and H$_2$ dissociation by cosmic ray, we found the cloud ages to be within 10$^{6.7}$-10$^{7.0}$ yr for five sources.

Precise measurements of the energy spectra of cosmic rays (CRs) show various kinds of features deviating from single power-laws, which give very interesting and important implications on their origin and propagation. Previous measurements from a few balloon and space experiments indicate the existence of spectral softenings around 10 TV for protons (and probably also for Helium nuclei). Very recently, the DArk Matter Particle Explorer (DAMPE) measurement about the proton spectrum clearly reveals such a softening with a high significance. Here we study the implications of these new measurements, as well as the groundbased indirect measurements, on the origin of CRs. We find that a single component of CRs fails to fit the spectral softening and the air shower experiment data simultaneously. In the framework of multiple components, we discuss two possible scenarios, the multiple source population scenario and the background plus nearby source scenario. Both scenarios give reasonable fits to the wide-band data from TeV to 100 PeV energies. Considering the anisotropy observations, the nearby source model is favored.

We observe the dust continuum at 225 GHz and CO isotopologue (12CO, 13CO, and C18O) J=2-1 emission lines toward the CR Cha protoplanetary disk using the Atacama Large Millimeter/Submillimeter Array (ALMA). The dust continuum image shows a dust gap-ring structure in the outer region of the dust disk. A faint dust ring is also detected around 120 au beyond the dust gap. The CO isotopologue lines indicate that the gas disk is more extended than the dust disk. The peak brightness temperature of the 13CO line shows a small bump around 130 au while 12CO and C18O lines do not show. We investigate two possible mechanisms for reproducing the observed dust gap-ring structure and a gas temperature bump. First, the observed gap structure can be opened by a Jupiter mass planet using the relation between the planet mass and the gap depth and width. Meanwhile, the radiative transfer calculations based on the observed dust surface density profile show that the observed dust ring could be formed by dust accumulation at the gas temperature bump, that is, the gas pressure bump produced beyond the outer edge of the dust disk.

The sunspot record is the only observational tracer of solar activity that provides a fundamental, multi-century reference. Its homogeneity has been largely maintained with a succession of long-duration visual observers. In this paper, we examine observations of one of the primary reference sunspot observers, Hisako Koyama. By consulting original archives of the National Museum of Nature and Science of Japan (hereafter, NMNS), we retrace the main steps of her solar-observing career, from 1945 to 1996. We also present the reconstruction of a full digital database of her sunspot observations at the NMNS, with her original drawings and logbooks. Here, we extend the availability of her observational data from 1947-1984 to 1945-1996. Comparisons with the international sunspot number (version 2) and with the group sunspot number series show a good global stability of Koyama's observations, with only temporary fluctuations over the main interval 1947-1982. Identifying drawings made by alternate observers throughout the series, we find that a single downward baseline shift in the record coincides with the partial contribution of replacement observers mostly after 1983. We determine the correction factor to bring the second part (1983-1996) to the same scale with Koyama's main interval (1947-1982). We find a downward jump by 9% after 1983, which then remains stable until 1996. Overall, the high quality of Koyama's observations with her life-long dedication leaves a lasting legacy of this exceptional personal achievement. With this comprehensive recovery, we now make the totality of this legacy directly accessible and exploitable for future research.

We have studied the high-mass star-forming complex NGC 6334 with ALMA in the continuum emission at a frequency of 87.6 GHz, achieving a spatial resolution of 1300 au. Detecting 142 compact sources distributed over the whole observed area, we then used machine learning algorithms to group the compact cores in different clusters. A total of four main clusters were identified: NGC 6334-E, NGC 6334-I, NGC 6334-I(N) and NGC 6334-I(NW). The typical separations between cluster members (4000-12000 au) together with the core masses (0.2-100 $M_{\odot}$) are in agreement with turbulent fragmentation at scales of 0.1 pc. We find that the CMFs (core mass functions) show an excess of massive cores compared to the IMF. Typical uncertainties in temperature or unresolved multiplicity may mimic the observed excess of massive cores. Evidence of mass segregation in NGC 6334-I and NGC 6334-I(N), with the most massive cores located closer to the center, was found. Correlations between the physical properties of the clusters and their evolutionary stage were searched for, and found to show a larger separation between cores in the more evolved clusters, favouring the role of gas expulsion and stellar ejection over time. From these results, we suggest that NGC 6334-I(N) was primordially segregated, while NGC 6334-I may have become mass-segregated over time due to dynamical effects. Finally, the lack of massive cores in the most evolved cluster suggests that the gas reservoir may be already exhausted, while the less evolved clusters still have access to a larger gas mass reservoir. In general, the fragmentation in NGC 6334 at large scales (about 1 pc) is governed by turbulent pressure, while at smaller scales (a few hundred au), thermal pressure regulates the fragmentation process.

FU Orionis objects are low-mass pre-main sequence stars characterized by dramatic outbursts of several magnitudes in brightness. These outbursts are linked to episodic accretion events in which stars gain a significant portion of their mass. The physical processes behind these accretion events are not yet well understood. The archetypical FU Ori system, FU Orionis, is composed of two young stars with detected gas and dust emission. The continuum emitting regions have not been resolved until now. Here, we present 1.3 mm observations of the FU Ori binary system with ALMA. The disks are resolved at 40 mas resolution. Radiative transfer modeling shows that the emission from FU Ori north (primary) is consistent with a dust disk with a characteristic radius of $\sim$11 au. The ratio between major and minor axes shows that the inclination of the disk is $\sim$37 deg. FU Ori south is consistent with a dust disk of similar inclination and size. Assuming the binary orbit shares the same inclination angle as the disks, the deprojected distance between north and south components is 0.6'', i.e. $\sim$250 au. Maps of $^{12}$CO emission show a complex kinematic environment with signatures disk rotation at the location of the northern component, and also (to a lesser extent) for FU Ori south. The revised disk geometry allows us to update FU Ori accretion models (Zhu et al.), yielding a stellar mass and mass accretion rate of FU Ori north of 0.6 M$_{\odot}$ and 3.8$\times10^{-5}$ M$_{\odot}$ yr$^{-1}$, respectively.

In order to study short timescale optical variability of $\gamma$-ray blazar S5 0716+714, quasi-simultaneous spectroscopic and multi-band photometric observations were performed from 2018 November to 2019 March with the 2.4 m optical telescope located at Lijiang Observatory of Yunnan Observatories. The observed spectra are well fitted with a power-law $F_{\lambda}=A\lambda ^{-\alpha}$ (spectral index $\alpha >0$). Correlations found between $\dot{\alpha}$, $\dot{A}$, $\dot{A}/A$, $\dot{F_{\rm{\lambda}}}$, and $\dot{F_{\rm{\lambda}}}/F_{\rm{\lambda}}$ are consistent with the trend of bluer-when-brighter (BWB). \textbf{The same case is for colors, magnitudes, color variation rates, and magnitude variation rates of photometric observations.} The variations of $\alpha$ lead those of $F_{\rm{\lambda}}$. Also, the color variations lead the magnitude variations. The observational data are mostly distributed in the I(+,+) and III(-,-) quadrants of coordinate system. Both of spectroscopic and photometric observations show BWB behaviors in S5 0716+714. The observed BWB may be explained by the shock-jet model, and its appearance may depend on the relative position of the observational frequency ranges with respect to the synchrotron peak frequencies, e.g., at the left of the peak frequencies. \textbf{Fractional variability amplitudes are $F_{\rm{var}}\sim 40\%$ for both of spectroscopic and photometric observations. Variations of $\alpha$ indicate variations of relativistic electron distribution producing the optical spectra. }

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.

As China's first X-ray astronomical satellite, the Hard X-ray Modulation Telescope (HXMT), which was dubbed as Insight-HXMT after the launch on June 15, 2017, is a wide-band (1-250 keV) slat-collimator-based X-ray astronomy satellite with the capability of all-sky monitoring in 0.2-3 MeV. It was designed to perform pointing, scanning and gamma-ray burst (GRB) observations and, based on the Direct Demodulation Method (DDM), the image of the scanned sky region can be reconstructed. Here we give an overview of the mission and its progresses, including payload, core sciences, ground calibration/facility, ground segment, data archive, software, in-orbit performance, calibration, background model, observations and some preliminary results.

The low energy (LE) X-ray telescope is one of the three main instruments of the Insight-Hard X-ray Modulation Telescope (Insight-HXMT). It is equipped with Swept Charge Device (SCD) sensor arrays with a total geometrical area of 384 cm2 and an energy band from 0.7 keV to 13 keV. In order to evaluate the particle induced X-ray background and the cosmic X-ray background simultaneously, LE adopts collimators to define four types of Field Of Views (FOVs). LE is constituted of three detector boxes (LEDs) and an electric control box (LEB) and achieves a good energy resolution of 140 eV at 5.9 keV, an excellent time resolution of 0.98 ms, as well as an extremely low pileup (<1% at 18000 cts/s). Detailed performance tests and calibration on the ground have been performed, including energy-channel relation, energy response, detection efficiency and time response.

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

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

Corona cooling was detected previously from stacking a series of short type-I bursts occurred during the low/had state of atoll outburst. Type-I bursts are hence regarded as sharp probe to our better understanding on the basic property of the corona. The launch of the first Chinese X-ray satellite Insight-HXMT has large detection area at hard X-rays which provide almost unique chance to move further in this research field. We report the first detection of the corona cooling by Insight-HXMT from single short type-I burst showing up during {\bf flare} of 4U 1636-536. This type-I X-ray burst has a duration of $\sim$13 seconds and hard X-ray shortage is detected with significance 6.2~$\sigma$ in 40-70 keV. A cross-correlation analysis between the lightcurves of soft and hard X-ray band, shows that the corona shortage lag the burst emission by 1.6 $\pm$1.2~s. These results are consistent with those derived previously from stacking a large amount of bursts detected by RXTE/PCA within a series of {\bf flares} of 4U 1636-536. Moreover, the broad bandwidth of Insight-HXMT allows as well for the first time to infer the burst influence upon the continuum spectrum via performing the spectral fitting of the burst, which ends up with the finding that hard X-ray shortage appears at around 40 keV in the continuum spectrum. These results suggest that the evolution of the corona along with the outburst{\bf /flare} of NS XRB may be traced via looking into a series of embedded type-I bursts by using Insight-HXMT.

The Medium Energy X-ray telescope (ME) is one of the three main telescopes on board the Insight Hard X-ray Modulation Telescope (Insight-HXMT) astronomy satellite. ME contains 1728 pixels of Si-PIN detectors sensitive in 5-30 keV with a total geometrical area of 952 cm2. Application Specific Integrated Circuit (ASIC) chips, VA32TA6, is used to achieve low power consumption and low readout noise. The collimators define three kinds of field of views (FOVs) for the telescope, 1{\deg}{\times}4{\deg}, 4{\deg}{\times}4{\deg}, and blocked ones. Combination of such FOVs can be used to estimate the in-orbit X-ray and particle background components. The energy resolution of ME is ~3 keV at 17.8 keV (FWHM) and the time resolution is 255 {\mu}s. In this paper, we introduce the design and performance of ME.

The Insight-Hard X-ray Modulation Telescope (Insight-HXMT) is a broad band X-ray and gamma-ray (1-3000 keV) astronomy satellite. The High Energy X-ray telescope (HE) is one of its three main telescopes. The main detector plane of HE is composed of 18 NaI(Tl)/CsI(Na) phoswich detectors, where NaI(Tl) serves as primary detector to measure ~ 20-250 keV photons incident from the field of view (FOV) defined by the collimators, and CsI(Na) is used as an active shield detector to NaI(Tl) by pulse shape discrimination. CsI(Na) is also used as an omnidirectional gamma-ray monitor. The HE collimators have a diverse FOV: 1.1{\deg}x 5.7{\deg} (15 units), 5.7{\deg}x 5.7{\deg} (2 units) and blocked (1 unit), thus the combined FOV of HE is about 5.7{\deg}x 5.7{\deg}. Each HE detector has a diameter of 190 mm, resulting in the total geometrical area of about 5100 cm_2. The energy resolution is ~15% at 60 keV. The timing accuracy is better than 10 {\mu}s and dead-time for each detector is less than 10 {\mu}s. HE is devoted to observe the spectra and temporal variability of X-ray sources in the 20-250 keV band either by pointing observations for known sources or scanning observations to unveil new sources, and to monitor the gamma-ray sky in 0.2-3 MeV. This paper presents the design and performance of the HE instruments. Results of the on-ground calibration experiments are also reported.

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

We present the observational results from a detailed timing analysis of the black hole candidate Swift J1658.2-4242 during its 2018 outburst with the observations of Hard X-ray Modulation Telescope (Insight-HXMT), Neutron Star Interior Composition Explorer (NICER) and AstroSat in 0.1-250keV. The evolution of intensity, hardness and integrated fractional root mean square (rms) observed by Insight-HXMT and NICER are presented in this paper. Type-C quasi-periodic oscillations (QPOs) observed by NICER (0.8-3.5Hz) and Insight-HXMT (1-1.6Hz) are also reported in this work. The features of the QPOs are analysed with an energy range of 0.5-50keV. The relations between QPO frequency and other characteristics such as intensity, hardness and QPO rms are carefully studied. The timing and spectral properties indicate that Swift J1658.2-4242 is a black hole binary system. Besides, the rms spectra of the source calculated from the simultaneous observation of Insight-HXMT, NICER and AstroSat support the Lense-Thirring origin of the QPOs. The relation between QPO phase lag and the centroid frequency of Swift J1658.2-4242 reveals a near zero constant when < 4Hz and a soft phase lag at 6.68 Hz. This independence follows the same trend as the high inclination galactic black hole binaries such as MAXI J1659-152.

We present a study of timing properties of the accreting pulsar 2S 1417-624 observed during its 2018 outburst, based on Swift/BAT, Fermi/GBM, Insight-HXMT and NICER observations. We report a dramatic change of the pulse profiles with luminosity. The morphology of the profile in the range 0.2-10.0keV switches from double to triple peaks at $\sim2.5$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$ and from triple to quadruple peaks at $\sim7$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$. The profile at high energies (25-100keV) shows significant evolutions as well. We explain this phenomenon according to existing theoretical models. We argue that the first change is related to the transition from the sub to the super-critical accretion regime, while the second to the transition of the accretion disc from the gas-dominated to the radiation pressure-dominated state. Considering the spin-up as well due to the accretion torque, this interpretation allows to estimate the magnetic field self-consistently at $\sim7\times 10^{12}$G.

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

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.

We report on analysis of observations of the bright transient X-ray pulsar \src obtained during its 2017-2018 giant outburst with Insight-HXMT, \emph{NuSTAR}, and \textit{Swift} observatories. We focus on the discovery of a sharp state transition of the timing and spectral properties of the source at super-Eddington accretion rates, which we associate with the transition of the accretion disk to a radiation pressure dominated (RPD) state, the first ever directly observed for magnetized neutron star. This transition occurs at slightly higher luminosity compared to already reported transition of the source from sub- to super-critical accretion regime associate with onset of an accretion column. We argue that this scenario can only be realized for comparatively weakly magnetized neutron star, not dissimilar to other ultra-luminous X-ray pulsars (ULPs), which accrete at similar rates. Further evidence for this conclusion is provided by the non-detection of the transition to the propeller state in quiescence which strongly implies compact magnetosphere and thus rules out magnetar-like fields.

For an orthogonal integral transform with complete dataset, any two components are linearly independent; however, when some data points are missing, there is going to be leakage from one component to another, which is referred to as the "leakage in integral transforms" in this work. A special case of this kind of leakage is the EB-leakage in detection of the cosmological gravitational wave background (CGWB). I first give the general solutions for all integral transforms, prove that they are the best solutions, and then apply them to the case of EB-leakage and detection of the CGWB. In the upcoming decade, most likely, new cosmic microwave background (CMB) data are from ground/balloon experiments, so they provide only partial sky coverage. Even in a fullsky mission, due to the Galactic foreground, part of the sky is still unusable. Within this context, the EB-leakage becomes inevitable. I show how to use the general solutions to achieve the minimal error bars of the EB-leakage, and use it to find out the maximum ability to detect the CGWB through CMB. The results show that, when focusing on the tensor-to-scalar ratio $r$ (at a pivot scale of 0.05 Mpc$^{-1}$), $1\%$ sky coverage ($f_{sky}=1\%$) is enough for a $5\sigma$-detection of $r\ge 10^{-2}$, but is barely enough for $r=10^{-3}$. If the target is to detect $r\sim 10^{-4}$ or $10^{-5}$, then $f_{sky}\ge 10\%$ is strongly recommended to enable a $5\sigma$-detection and to reserve some room for other errors.

We report high angular resolution observations of linearly polarized dust emission towards the Class 0 young stellar object (YSO) NGC1333 IRAS4A (hereafter, IRAS4A) using the Karl G. Jansky Very Large Array (JVLA) at K (11.5-16.7 mm), Ka (8.1-10.3 mm), and Q bands (6.3-7.9 mm), and using the Atacama Large Millimeter Array (ALMA) at Band 6 (1.2 mm) and Band 7 (0.85-0.89 mm). On 100-1000 AU scales, all of these observations consistently trace the hourglass shaped magnetic field topology as shown in the previous studies. In the innermost 100 AU region of IRAS4A1, the polarization position angles (E-field) detected at 6.3-16.7 mm are consistent, however, are nearly 90 degrees offset from those detected at 1.2 mm and 0.85-0.89 mm. Such a 90 degree offset may be explained by that the inner $\sim$100 AU area is optically thick at wavelengths shorter than $\sim$1.5 mm, whereby the observations probe the absorption of aligned dust against the weakly or unpolarized warm dust emission from the innermost region. This can also consistently explain why the highest angular resolution ALMA images at Band 7 show that the polarization percentage increases with dust brightness temperature in the inner $\sim$100 AU region of IRAS4A1. Following this interpretation and assuming that the dust grains are aligned with the magnetic fields, the inferred magnetic field position angle based on the 90$^{\circ}$ rotated at 6.3-7.9 mm in the central peak of IRAS4A1 is $\sim-22^{\circ}$, which is approximately consistent with the outflow direction $\sim-9^{\circ}$.

The dust properties in high-redshift galaxies provide clues to the origin of dust in the Universe. Although dust has been detected in galaxies at redshift $z>7$, it is difficult to constrain the dominant dust sources only from the total dust amount. Thus, we calculate the evolution of grain size distribution, expecting that different dust sources predict different grain size distributions. Using the star formation time-scale and the total baryonic mass constrained by the data in the literature, we calculate the evolution of grain size distribution. To explain the total dust masses in ALMA-detected $z>7$ galaxies, the following two solutions are possible: (i) high dust condensation efficiency in stellar ejecta, and (ii) efficient accretion (dust growth by accreting the gas-phase metals in the interstellar medium). We find that these two scenarios predict significantly different grain size distributions: in (i), the dust is dominated by large grains ($a\gtrsim 0.1~\mu$m, where $a$ is the grain radius), while in (ii), the small-grain ($a\lesssim 0.01~\mu$m) abundance is significantly enhanced by accretion. Accordingly, extinction curves are expected to be much steeper in (ii) than in (i). Thus, we conclude that extinction curves provide a viable way to distinguish the dominant dust sources in the early phase of galaxy evolution.

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.

We present ALMA observations of 12CO, 13CO, and C18O J=2--1 lines and the 230 GHz continuum for the FU Ori-type object (FUor) V900 Mon (d~1.5 kpc), for which the accretion burst was triggered between 1953 and 2009. We identified CO emission associated with a molecular bipolar outflow extending up to a ~10^4 au scale and a rotating molecular envelope extending over >10^4 au. The interaction with the hot energetic FUor wind, which was observed using optical spectroscopy, appears limited to a region within ~400 au of the star. The envelope mass and the collimation of the extended CO outflow suggest that the progenitor of this FUor is a low-mass Class I young stellar object (YSO). These parameters for V900 Mon, another FUor, and a few FUor-like stars are consistent with the idea that FUor outbursts are associated with normal YSOs. The continuum emission is marginally resolved in our observations with a 0."2x0."15 (~300x225 au) beam, and a Gaussian model provides a deconvolved FWHM of ~90 au. The emission is presumably associated with a dusty circumstellar disk, plus a possible contribution from a wind or a wind cavity close to the star. The warm compact nature of the disk continuum emission could be explained with viscous heating of the disk, while gravitational fragmentation in the outer disk and/or a combination of grain growth and their inward drift may also contribute to its compact nature.

Observations of young stellar objects (YSOs) in centimeter bands can probe the continuum emission from growing dust grains, ionized winds, and magnetospheric activity, which are intimately connected to the evolution of protoplanetary disks and the formation of planets. We have carried out sensitive continuum observations toward the Ophiuchus A star-forming region using the Karl G. Jansky Very Large Array (VLA) at 10 GHz over a field-of-view of 6$'$ with a spatial resolution of $\theta_{maj}$ $\times$ $\theta_{min}$ $\sim$ 0.4$''$ $\times$ 0.2$''$. We achieved a 5 $\mu$Jy beam$^{-1}$ root-mean-square noise level at the center of our mosaic field of view. Among the eighteen sources we detected, sixteen are YSOs (three Class 0, five Class I, six Class II, and two Class III) and two are extragalactic candidates. We find that thermal dust emission generally contributes less that 30% of the emission at 10 GHz. The radio emission is dominated by other types of emission such as gyro-synchrotron radiation from active magnetospheres, free-free emission from thermal jets, free-free emission from the outflowing photo-evaporated disk material, and/or synchrotron emission from accelerated cosmic-rays in jet or protostellar surface shocks. These different types of emission could not be clearly disentangled. Our non-detections towards Class II/III disks suggest that extreme UV-driven photoevaporation is insufficient to explain the disk dispersal, assuming that the contribution of UV photoevaporating stellar winds to radio flux does not evolve with time. The sensitivity of our data cannot exclude photoevaporation due to X-ray photons as an efficient mechanism for disk dispersal. Deeper surveys with the Square Kilometre Array will be able to provide strong constraints on disk photoevaporation.

We present initial results from the K-band focal plane array Examinations of Young STellar Object Natal Environments (KEYSTONE) survey, a large project on the 100-m Green Bank Telescope mapping ammonia emission across eleven giant molecular clouds at distances of $0.9-3.0$ kpc (Cygnus X North, Cygnus X South, M16, M17, MonR1, MonR2, NGC2264, NGC7538, Rosette, W3, and W48). This data release includes the NH$_3$ (1,1) and (2,2) maps for each cloud, which are modeled to produce maps of kinetic temperature, centroid velocity, velocity dispersion, and ammonia column density. Median cloud kinetic temperatures range from $11.4\pm2.2$ K in the coldest cloud (MonR1) to $23.0\pm6.5$ K in the warmest cloud (M17). Using dendrograms on the NH$_3$ (1,1) integrated intensity maps, we identify 856 dense gas clumps across the eleven clouds. Depending on the cloud observed, $40-100\%$ of the clumps are aligned spatially with filaments identified in H$_2$ column density maps derived from SED-fitting of dust continuum emission. A virial analysis reveals that 523 of the 835 clumps ($\sim63\%$) with mass estimates are bound by gravity alone. We find no significant difference between the virial parameter distributions for clumps aligned with the dust-continuum filaments and those unaligned with filaments. In some clouds, however, hubs or ridges of dense gas with unusually high mass and low virial parameters are located within a single filament or at the intersection of multiple filaments. These hubs and ridges tend to host water maser emission, multiple 70$\mu$m-detected protostars, and have masses and radii above an empirical threshold for forming massive stars.

The existence of millicharged dark matter (mDM) can leave a measurable imprint on 21-cm cosmology through mDM-baryon scattering. However, the minimal scenario is severely constrained by existing cosmological bounds on both the fraction of dark matter that can be millicharged and the mass of mDM particles. We point out that introducing a long-range force between a millicharged subcomponent of dark matter and the dominant cold dark matter (CDM) component leads to efficient cooling of baryons in the early universe, while also significantly extending the range of viable mDM masses. Such a scenario can explain the anomalous absorption signal in the sky-averaged 21-cm spectrum observed by EDGES, and leads to a number of testable predictions for the properties of the dark sector. The mDM mass can then lie between 10 MeV and a few hundreds of GeVs, and its scattering cross section with baryons lies within an unconstrained window of parameter space above direct detection limits and below current bounds from colliders. In this allowed region, mDM can make up as little as $10^{-8}$ of the total dark matter energy density. The CDM mass ranges from 10 MeV to a few GeVs, and has an interaction cross section with the Standard Model that is induced by a loop of mDM particles. This cross section is generically within reach of near-future low-threshold direct detection experiments.

The Solar System harbors deep unresolved mysteries despite centuries-long study. A highly intriguing case concerns anomalous non-dipolar and non-axisymmetric magnetic fields of Uranus and Neptune that have long eluded explanation by the prevailing theory. A thin-shell dynamo conjecture captures observed phenomena but leaves unexplained fundamental material basis and underlying mechanism. Here, we report the discovery of trihydrogen oxide (H3O) in metallic liquid state stabilized at extreme pressure and temperature conditions inside these icy planets. Calculated stability pressure field compared to known pressure-radius relation for Uranus and Neptune places metallic liquid H3O in a thin-shell zone near planetary cores. These findings from accurate quantum mechanical calculations rationalize the empirically conjectured thin-shell dynamo model and establish key physical benchmarks that are essential to elucidating the enigmatic magnetic-field anomaly of Uranus and Neptune, resolving a major mystery in planetary science.

L1551 IRS 5 is a FUor-like object located in the Taurus star forming region. We present ALMA 1.3 mm continuum observations using a wide range of baselines. The observations recovered the two circumstellar disks composing the system and, for the first time, resolved the circumbinary ring. We determined the geometry and estimated lower mass limits for the circumstellar disks using simple models. We calculated lower limits for the total mass of both circumstellar disks. After subtracting the two circumstellar disk models from the image, the residuals show a clearly resolved circumbinary ring. Using a radiative transfer model, we show that geometrical effects can explain some of the brightness asymmetries found in the ring. The remaining features are interpreted as enhancements in the dust density.

We report new Atacama Large Millimeter/submillimeter Array Band 3 (86-100 GHz; $\sim$80 mas angular resolution) and Band 4 (146-160 GHz; $\sim$50 mas angular resolution) observations of the dust continuum emission towards the archetypal and ongoing accretion burst young stellar object FU Ori, which simultaneously covered its companion, FU Ori S. In addition, we present near-infrared (2-2.45 $\mu$m) observations of FU Ori taken with the General Relativity Analysis via VLT InTerferometrY (GRAVITY; $\sim$1 mas angular resolution) instrument on the Very Large Telescope Interferometer (VLTI). We find that the emission in both FU Ori and FU Ori S at (sub)millimeter and near infrared bands is dominated by structures inward of $\sim$10 au radii. We detected closure phases close to zero from FU Ori with VLTI/GRAVITY, which indicate the source is approximately centrally symmetric and therefore is likely viewed nearly face-on. Our simple model to fit the GRAVITY data shows that the inner 0.4 au radii of the FU Ori disk has a triangular spectral shape at 2-2.45 $\mu$m, which is consistent with the H$_{2}$O and CO absorption features in a $\dot{M}\sim$10$^{-4}$ $M_{\odot}\,yr^{-1}$, viscously heated accretion disk. At larger ($\sim$0.4-10 au) radii, our analysis shows that viscous heating may also explain the observed (sub)millimeter and centimeter spectral energy distribution when we assume a constant, $\sim$10$^{-4}$ $M_{\odot}\,yr^{-1}$ mass inflow rate in this region. This explains how the inner 0.4 au disk is replenished with mass at a modest rate, such that it neither depletes nor accumulates significant masses over its short dynamic timescale. Finally, we tentatively detect evidence of vertical dust settling in the inner 10 au of the FU Ori disk, but confirmation requires more complete spectral sampling in the centimeter bands.

We present the photometric observations of blazars S5 0716+714 and 3C 273 with high temporal resolution (30--60s) in the $I$ or $R$ bands. The observations were performed with a 1.02 m optical telescope from 2007 March 07 to 2012 May 16. The $F$-test, one-way analysis of variance (ANOVA) test, and $z$-transformed discrete correlation function (ZDCF) cross-correlation zero lag test are used to search for intra-day variability (IDV). Four and five reliable IDVs survive three tests for S5 0716+714 and 3C 273, respectively. IDVs are found for S5 0716+714 and 3C 273. A flare on 2008 May 08 has $\Delta I \approx$ 0.06$\pm$0.01 mag in a duration of 0.54 hr for S5 0716+714. A flare on 2011 May 10 shows $\Delta R \approx$ 0.05$\pm$0.01 mag in a duration of 0.40 hr for 3C 273. Sharp dips appear on 2011 May 9 for 3C 273, and show $\Delta R \approx$ 0.05$\pm$0.01 mag. Under the assumptions that the IDV is tightly connected to black hole mass, $M_{\bullet}$, and that the flare durations are representative of the minimum characteristic timescales, we can estimate upper bounds to $M_{\bullet}$. In the case of the Kerr black holes, $M_{\bullet} \la 10^{8.91} M_{\odot}$ and $M_{\bullet} \la 10^{9.02} M_{\odot}$ are given for S5 0716+714 and 3C 273, respectively. These mass measurements are consistent with those measurements reported in the literatures. Also, we discuss the origins of optical variations found in this work.

In this work, we report the in-orbit demonstration of X-ray pulsar navigation with Insight-Hard X-ray Modulation Telescope (Insight-HXMT), which was launched on Jun. 15th, 2017. The new pulsar navigation method 'Significance Enhancement of Pulse-profile with Orbit-dynamics' (SEPO) is adopted to determine the orbit with observations of only one pulsar. In this test, the Crab pulsar is chosen and observed by Insight-HXMT from Aug. 31th to Sept. 5th in 2017. Using the 5-day-long observation data, the orbit of Insight-HXMT is determined successfully with the three telescopes onboard - High Energy X-ray Telescope (HE), Medium Energy X-ray Telescope (ME) and Low Energy X-ray Telescope (LE) - respectively. Combining all the data, the position and velocity of the Insight-HXMT are pinpointed to within 10 km (3 sigma) and 10 m/s (3 sigma), respectively.

Axion-like particles are a broad class of dark matter candidates which are expected to behave as a coherent, classical field with a weak coupling to photons. Research into the detectability of these particles with laser interferometers has recently revealed a number of promising experimental designs. Inspired by these ideas, we propose the Axion Detection with Birefringent Cavities (ADBC) experiment, a new axion interferometry concept using a cavity that exhibits birefringence between its two, linearly polarized laser eigenmodes. This experimental concept overcomes several limitations of the designs currently in the literature, and can be practically realized in the form of a simple bowtie cavity with tunable mirror angles. Our design thereby increases the sensitivity to the axion-photon coupling over a wide range of axion masses.

In observation of the cosmic microwave background (CMB) polarization, "$EB$~leakage" refers to the artificial $B$-mode signal coming from the leakage of $E$-mode signal when part of the sky is unavailable or excluded. Correction of such leakage is one of the preconditions for detecting primordial gravitational waves via the CMB $B$-mode signal. In this work, we design two independent methods for correcting the $EB$~leakage directly in the pixel domain using standard definitions of the $E$- and $B$-modes. The two methods give consistent results, and both are fast and easy to implement. Tests on a CMB simulation containing zero initial $B$-mode show an efficient suppression of the $EB$ leakage. When combined with the MASTER method to reconstruct the full-sky $B$-mode spectrum in simulations with a relatively simple mask, the error from EB-leakage is suppressed further by more than one order of magnitude at the recombination bump, and up to three orders of magnitude at higher multipoles, compared to a "pure" MASTER scheme under the same conditions. Meanwhile, although the final power spectrum estimation benefits from apodization, the pixel domain correction itself is done without apodization, and thus the methods offer more freedom in choosing an apodization based on specific requirements.

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.

A new, complete sample of 14,584 broad-line AGNs at $z<0.35$ is presented, which are uncovered homogeneously from the complete database of galaxies and quasars observed spectroscopically in the Sloan Digital Sky Survey Seventh Data Release. The stellar continuum is properly removed for each spectrum with significant host absorption line features, and careful analyses of the emission-line spectra, particularly in the H$\rm \alpha$ and H$\rm \beta$ wavebands, are carried out. The broad Balmer emission line, particularly, H$\rm \alpha$, is used to indicate the presence of an AGN. The broad H$\rm \alpha$ lines have luminosities in a range of $10^{38.5}$-$10^{44.3}$ erg s$^{-1}$, and line widths (FWHMs) of 500-34,000 km s$^{-1}$. The virial black hole masses, estimated from the broad line measurements, span a range of $10^{5.1}$-$10^{10.3}$ $M_\odot$, and the Eddington ratios vary from $-3.3$ to $1.3$ in logarithmic scale. Other quantities such as multi-wavelength photometric properties and flags denoting peculiar line profiles are also included in this catalog. We describe the construction of this catalog and briefly discuss its properties. The catalog is publicly available online. This homogeneously selected AGN catalog, along with the accurately measured spectral parameters, provide the most updated, largest AGN sample data, which will enable further comprehensive investigations of the properties of the AGN population in the low-redshift universe.

The recently discovered neutron star transient Swift J0243.6+6124 has been monitored by {\it the Hard X-ray Modulation Telescope} ({\it Insight-\rm HXMT). Based on the obtained data, we investigate the broadband spectrum of the source throughout the outburst. We estimate the broadband flux of the source and search for possible cyclotron line in the broadband spectrum. No evidence of line-like features is, however, found up to $\rm 150~keV$. In the absence of any cyclotron line in its energy spectrum, we estimate the magnetic field of the source based on the observed spin evolution of the neutron star by applying two accretion torque models. In both cases, we get consistent results with $B\rm \sim 10^{13}~G$, $D\rm \sim 6~kpc$ and peak luminosity of $\rm >10^{39}~erg~s^{-1}$ which makes the source the first Galactic ultraluminous X-ray source hosting a neutron star.

The discovery of nongravitational interactions between dark matter and the Standard Model would be an important step in unraveling the nature of dark matter. If such an interaction exists, it would have profound implications on how dark matter is produced in both the early universe and in collider experiments. In addition, it would also allow dark matter to deposit energy into Standard Model particles in unexpected ways. This thesis details some recent progress made in understanding these implications, including (i) a new freezeout mechanism for thermal dark matter dominated by a 3-to-2 process within a vector portal dark sector model; (ii) a study of how the existence of dark sector bound states can influence collider, direct and indirect searches for dark matter; (iii) a new axion dark matter interferometric search using a cavity that is sensitive to the axion-induced rotation of linearly polarized light; (iv) a definitive assessment of the potential contribution of dark matter annihilation and decay to cosmic reionization; (v) new constraints on dark matter annihilation rates and decay lifetimes from 21-cm cosmology, and (vi) a new numerical code, DarkHistory, which significantly improves the computation of the ionization and thermal histories of the universe in the presence of exotic sources of energy injection. These novel ideas span length scales ranging from table-top experiments to the entire cosmos, and represent just a few of the myriad of ways in which dark matter may yet surprise us.

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.

Time delay in Schwarzschild spacetime for null and timelike signals with arbitrary velocity $v$ is studied. The total travel time $t_\mathrm{if}$ is evaluated both exactly and approximately in the weak field limit, with the result given as functions of signal velocity, source-lens and lens-observer distances, angular position of the source and lens mass. Two time delays, $\Delta t_v$ between signals with different velocities but coming from same side of the lens and $\Delta t_\mathrm{p}$ between signals from different sides of the lens, as well as the difference $\Delta t_{\mathrm{p}v}$ between two $\Delta t_\mathrm{p}$'s are calculated. These time delays are applied to the gravitational-lensed supernova neutrinos and gravitational waves (GW). It is shown that the $\Delta t_v$ between different mass eigenstates of supernova neutrinos can be related to the mass square difference of these eigenstates and therefore could potentially be used to discriminate neutrino mass orderings, while the difference $\Delta t_{\mathrm{p}v}$ between neutrino and optical signals can be correlated with the absolute mass of neutrinos. The formula for time delay in a general lens mass profile is derived and the result is applied to the singular isothermal sphere case. For GWs, it is found that the difference $\Delta t_{\mathrm{p}v}$ between GW and GRB can only reach $1.45\times 10^{-5}$ second for very large source distance ($2\times 10^4$ [Mpc]) and source angle (10 [as]) if $v_{GM}=(1-3\times 10^{-15})c$. This time difference is at least three order smaller than uncertainties in time measurement of the recently observed GW/GRB signals and thus calls for improvement if $\Delta t_{\mathrm{p}v}$ is to be used to further constrain the GW velocity.

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 present our observational results of the 1.1 mm continuum and the HCO+ (3-2) line in HL Tau at angular resolutions of 0.1" obtained with ALMA and our data analysis of the 2.9 mm and 1.1 mm continuum and the HCO+ (3-2) and (1-0) lines of the HL Tau disk. The Keplerian rotation of the HL Tau disk is well resolved in the HCO+ (3-2) emission, and the stellar mass is estimated to be 2.1+/-0.2 Msun with a disk inclination angle of 47 deg. The radial profiles of the HCO+ column density and excitation temperature are measured with the LTE analysis of the two transitions of the HCO+ emission. An HCO+ gas gap at a radius of 30 au, where the column density drops by a factor of 4-8, is found in the HCO+ column density profile, coincident with the dust gap traced by the continuum emission. No other clear HCO+ gas gap is seen. This HCO+ gas gap can be opened by a planet with mass of 0.5-0.8 Mj, which is comparable to the planet mass adopted in numerical simulations to form the dust gap at the same radius in the HL Tau disk. In addition to the disk component, a one-arm spiral with a length of ~3" (520 au) stretching out from the inner disk is observed in the HCO+ (3-2) emission. The observed velocity structures along the spiral suggest an infalling and rotating gas stream toward the inner disk.

Foreground removal techniques for CMB analyses make specific assumptions about the properties of foregrounds in temperature and in polarization. By investigating the statistics of foreground components more understanding about the degree to which these assumptions are valid can be obtained. In this work we investigate $E$- and $B$-mode maps of the two strongest polarized foregrounds, synchrotron and thermal dust emission, with regards to their similarity with Gaussian processes, their spectral variations and cross-correlations. We perform tests in patches of $\sim3.7^\circ$ size collectively covering the full sky and find most of them to conform with their Gaussian expectation according to the statistics in use. Correlations exhibit distinct differences in $E$- and $B$-mode signals which point towards necessities in foreground removal methods. We discuss potential consequences and possible further directions.

We continue exploring the constraining capabilities of X-ray reflection spectroscopy to test the Kerr-nature of astrophysical black holes and we present the results of our analysis of two NuSTAR observations of Cygnus X-1 in the soft state. We find that the final measurement can strongly depend on the assumption of the intensity profile. We conclude that Cygnus X-1 is not suitable for tests of general relativity using X-ray reflection spectroscopy and we discuss the desired properties of a source to be a good candidate for our studies.

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.

The broadband spectral energy distribution (SED) of Active Galactic Nuclei (AGN) is investigated for a well-selected sample composed of $23$ Seyfert 1 galaxies observed simultaneously in the optical/UV and X-ray bands with the Neil Gehrels {\it Swift} Observatory. The optical to UV continuum spectra are modeled, for the first time, with emission from an accretion disk with a generalized radial temperature profile, in order to account for the intrinsic spectra which are found to be generally redder than the model prediction of the standard Shakura-Sunyaev disk (SSD) ($F_\nu\propto\nu^{+1/3}$). The power-law indices of the radial temperature profile ($T_{\rm eff}(R)\propto R^{-p}$, $R$ is the radius of the accretion disk) are inferred to be $p=0.5$ -- $0.75$ (a median of $0.63$), deviating from the canonical $p=0.75$ for the SSD model as widely adopted in previous studies. A marginal correlation of a flatter radial temperature profile (a smaller $p$ value) with increasing the Eddington ratio is suggested. Such a model produces generally a lower peak of accretion disk emission and thus a smaller bolometric luminosity in some of the AGN, particularly those with high Eddington ratios, than that based on the SSD model by a factor of several. The broadband SED, the bolometric correction factors and their dependence on some of the AGN parameters are re-visited. We suggest that such non-standard SSD disks may operate in AGN and are at least partly responsible for the reddened optical/UV spectra as observed. One possible explanation for these flattened temperature profiles is the mass loss process in form of disk winds/outflows.

The existence of possible massive white dwarfs more than the Chandrasekhar limit ($1.45M_{\odot }$, in which $M_{\odot }$ is mass of the sun) is a challenging topic. In this regard and motivated by the important effect of massive graviton on the structure of white dwarfs we study the white dwarfs in Vegh's massive gravity which is known as one of theories of de Rham, Gabadadze and Tolley (dRGT) like massive gravity. First, we consider the modified Tolman-Oppenheimer-Volkoff equation in this theory of massive gravity and solve it numerically by using the Chandrasekhar's equation of state. Our results show that the maximum mass of white dwarfs in massive gravity can be more than the Chandrasekhar limit ($M>1.45M_{\odot }$ ), and this result imposes some constraints on parameters of massive gravity. Then, we investigate the effects of various parameters on other properties of the white dwarfs such as mass-radius relation, mass-central density relation, Schwarzschild radius, average density and Kretschmann scalar. Next, we study dynamical stability condition for super-Chandrasekhar white dwarfs and show that these massive compact objects enjoy dynamical stability. Finally, in order to have a better insight, we compare the super-Chandrasekhar white dwarfs with the obtained massive neutron stars in dRGT like massive theory of gravity.

We present new ALMA observations of the molecular gas and far-infrared continuum around the brightest cluster galaxy (BCG) in the cool-core cluster MACS 1931.8-2635. Our observations reveal $1.9 \pm 0.3 \times 10^{10}$ M$_{\odot}$ of molecular gas, on par with the largest known reservoirs of cold gas in a cluster core. We detect CO(1-0), CO(3-2), and CO(4-3) emission from both diffuse and compact molecular gas components that extend from the BCG center out to $\sim30$ kpc to the northwest, tracing the UV knots and H$\alpha$ filaments observed by HST. Due to the lack of morphological symmetry, we hypothesize that the $\sim300$ km s$^{-1}$ velocity of the CO in the tail is not due to concurrent uplift by AGN jets, rather we may be observing the aftermath of a recent AGN outburst. The CO spectral line energy distribution suggests that molecular gas excitation is influenced by processes related to both star formation and recent AGN feedback. Continuum emission in Bands 6 and 7 arises from dust and is spatially coincident with young stars and nebular emission observed in the UV and optical. We constrain the temperature of several dust clumps to be $\lesssim 10$ K, which is too cold to be directly interacting with the surrounding $\sim 4.8$ keV intracluster medium (ICM). The cold dust population extends beyond the observed CO emission and must either be protected from interacting with the ICM or be surrounded by local volumes of ICM that are several keV colder than observed by Chandra.

Using the new equipment of the Shanghai Tian Ma Radio Telescope, we have searched for carbon-chain molecules (CCMs) towards five outflow sources and six Lupus I starless dust cores, including one region known to be characterized by warm carbon-chain chemistry (WCCC), Lupus I-1 (IRAS 15398-3359), and one TMC-1 like cloud, Lupus I-6 (Lupus-1A). Lines of HC3N J=2-1, HC5N J=6-5, HC7N J=14-13, 15-14, 16-15 and C3S J=3-2 were detected in all the targets except in the outflow source L1660 and the starless dust core Lupus I-3/4. The column densities of nitrogen-bearing species range from 10$^{12}$ to 10$^{14}$ cm$^{-2}$ and those of C$_3$S are about 10$^{12}$ cm$^{-2}$. Two outflow sources, I20582+7724 and L1221, could be identified as new carbon-chain--producing regions. Four of the Lupus I dust cores are newly identified as early quiescent and dark carbon-chain--producing regions similar to Lup I-6, which together with the WCCC source, Lup I-1, indicate that carbon-chain-producing regions are popular in Lupus I which can be regard as a Taurus like molecular cloud complex in our Galaxy. The column densities of C3S are larger than those of HC7N in the three outflow sources I20582, L1221 and L1251A. Shocked carbon-chain chemistry (SCCC) is proposed to explain the abnormal high abundances of C3S compared with those of nitrogen-bearing CCMs. Gas-grain chemical models support the idea that shocks can fuel the environment of those sources with enough $S^+$ thus driving the generation of S-bearing CCMs.

We detected carbon-chain molecules (CCMs) HC2n+1N (n=1-3) and C3S in Ku band as well as {high-energy excitation lines including C4H N=9-8, J=17/2-15/2, 19/2-17/2, and CH3CCH J=5-4, K=2 in the 3 mm band toward a starless core called the eastern molecular core (EMC) of L1489 IRS. Maps of all the observed lines were also obtained. Comparisons with a number of early starless cores and WCCC source L1527 show that the column densities of C4H and CH3CCH are close to those of L1527, and the CH3CCH column densities of the EMC and L1527 are slightly higher than those of TMC-1. The EMC and L1527 have similar C3S column densities, but they are much lower than those of all the starless cores, with only 6.5% and 10% of the TMC-1 value, respectively. The emissions of the N-bearing species of the EMC and L1527 are at the medium level of the starless cores. These comparisons show that the CCM emissions in the EMC are similar to those of L1527, though L1527 contains a protostar. Although dark and quiescent, the EMC is warmer and at a later evolutionary stage than classical carbon-chain--producing regions in the cold, dark, quiescent early phase. The PACS, SPIRE, and SCUBA maps evidently show that the L1489 IRS seems to be the heating source of the EMC. Although it is located at the margins of the EMC, its bolometric luminosity and bolometric temperature are relatively high. Above all, the EMC is a rather particular carbon-chain-producing region and is quite significant for CCM science.

We use APEX mapping observations of 13CO, and C18O (2-1) to investigate the internal gas kinematics of the filamentary cloud G350.54+0.69, composed of the two distinct filaments G350.5-N and G350.5-S. G350.54+0.69 as a whole is supersonic and gravitationally bound. We find a large-scale periodic velocity oscillation along the entire G350.5-N filament with a wavelength of ~1.3 pc and an amplitude of ~0.12 km/s. Comparing with gravitational-instability induced core formation models, we conjecture that this periodic velocity oscillation could be driven by a combination of longitudinal gravitational instability and a large-scale periodic physical oscillation along the filament. The latter may be an example of an MHD transverse wave. This hypothesis can be tested with Zeeman and dust polarization measurements.

We present a long short-term memory (LSTM) network for predicting whether an active region (AR) would produce a gamma-class flare within the next 24 hours. We consider three gamma classes, namely >=M5.0 class, >=M class, and >=C class, and build three LSTM models separately, each corresponding to a gamma class. Each LSTM model is used to make predictions of its corresponding gamma-class flares. The essence of our approach is to model data samples in an AR as time series and use LSTMs to capture temporal information of the data samples. Each data sample has 40 features including 25 magnetic parameters obtained from the Space-weather HMI Active Region Patches (SHARP) and related data products as well as 15 flare history parameters. We survey the flare events that occurred from 2010 May to 2018 May, using the GOES X-ray flare catalogs provided by the National Centers for Environmental Information (NCEI), and select flares with identified ARs in the NCEI flare catalogs. These flare events are used to build the labels (positive vs. negative) of the data samples. Experimental results show that (i) using only 14-22 most important features including both flare history and magnetic parameters can achieve better performance than using all the 40 features together; (ii) our LSTM network outperforms related machine learning methods in predicting the labels of the data samples. To our knowledge, this is the first time that LSTMs have been used for solar flare prediction.

We present NKBB, an XSPEC model for the thermal spectrum of thin accretion disks in parametric black hole spacetimes. We employ the Novikov-Thorne model for the description of the accretion disk and the formalism of the transfer function proposed by Cunningham for storing all the relativistic effects of the spacetime metric. The current version of the model assumes the Johannsen metric, but it can be easily extended to any stationary, axisymmetric, and asymptotically-flat spacetime without pathological properties. The model can be used within the XSPEC package to test the Kerr nature of astrophysical black holes with the continuum-fitting method.

We reexamine the results presented in a recent work by Nielsen et al. [1], in which the properties of the noise residuals in the 40\,ms chirp domain of GW150914 were investigated. This paper confirmed the presence of strong (i.e., about 0.80) correlations between residual noise in the Hanford and Livingston detectors in the chirp domain as previously seen by us [2] when using a numerical relativity template given in [3]. It was also shown in [1] that a so-called maximum likelihood template can reduce these statistically significant cross-correlations. Here, we demonstrate that the reduction of correlation and statistical significance is due to (i) the use of a peculiar template which is qualitatively different from the properties of GW150914 originally published by LIGO, (ii) a suspicious MCMC chain, (iii) uncertainties in the matching of the maximum likelihood (ML) template to the data in the Fourier domain, and (iv) a biased estimation of the significance that gives counter-intuitive results. We show that rematching the maximum likelihood template to the data in the 0.2\,s domain containing the GW150914 signal restores these correlations at the level of $60\%$ of those found in [1]. With necessary corrections, the probability given in [1] will decrease by more than one order of magnitude. Since the ML template is itself problematic, results associated with this template are illustrative rather than final.

The structured jet is involved to explain the afterglows and even the prompt emission of GRB 170817A. In this paper, we stress that for a precessing jet, the jet structure in the prompt emission phase and that in the afterglow phase may be different. The jet structure in the afterglow phase can be non-uniform even if a narrow-uniform jet is presented in the prompt emission phase. We estimate the jet structure in the afterglow phase under the situation that a narrow-uniform-precessing jet is launched from the central engine of gamma-ray burst. With different precession angles, it is found that the structured jet can be roughly described as follows: a narrow uniform core with power-law wings and sharp cut-off edges, a Gaussian profile, a ring shape, or other complex profile in energy per solid angle. Correspondingly, the afterglows for our obtained structured jets are also estimated. We find that the estimates of the intrinsic kinetic energy, the electron index, and the jet opening angle based on the afterglows formed in a precessing system may be incorrect. Our obtained structured jet is likely to be revealed by future observations for a fraction of gravitational wave detected merging compact binary systems (e.g., black hole-neutron star mergers).

Previous (sub)millimeter observations have found that the spectral indices of dust emission from some young stellar objects are lower than that of the black body emission in the Rayleigh-Jeans limit (i.e., 2.0). In particular, the recent Atacama Large Millimeter Array observations have spatially resolved that the innermost regions of the protoplanetary disks TW\,Hya and HD\,163296 present anomalously low (i.e., $<$2.0) millimeter spectral indices. In some previous works, such anomalously low millimeter spectral indices were considered unphysical and were attributed to measurement errors. The present work clarifies that if the albedo is high and is increasing with frequency, it is possible to reproduce such anomalously low spectral indices when the emission source is optically thick. In addition, to yield lower than 2.0 spectral index at (sub)millimeter bands, the required dust maximum grain size \amax is on the order of 10-100 \micron, which is very well consistent with the previously derived \amax values based on multi-wavelength dust polarimetric observations. In light of this, measuring Stokes I spectral index may also serve as an auxiliary approach for assessing whether the observed dust polarization is mainly due to dust scattering or is due to the aligned dust grains.

We report the discovery of extreme X-ray variability in a type 1 quasar: SDSS J$075101.42+291419.1$. It has a black hole mass of $1.6\times 10^7~\rm M_\odot$ measured from reverberation mapping (RM), and the black hole is accreting with a super-Eddington accretion rate. Its XMM-Newton observation in 2015 May reveals a flux drop by a factor of $\sim 22$ with respect to the Swift observation in 2013 May when it showed a typical level of X-ray emission relative to its UV/optical emission. The lack of correlated UV variability results in a steep X-ray-to-optical power-law slope ($\alpha_{\rm OX}$) of -1.97 in the low X-ray flux state, corresponding to an X-ray weakness factor of 36.2 at rest-frame 2 keV relative to its UV/optical luminosity. The mild UV/optical continuum and emission-line variability also suggest that the accretion rate did not change significantly. A single power-law model modified by Galactic absorption describes well the $0.3-10$ keV spectra of the X-ray observations in general. The spectral fitting reveals steep spectral shapes with $\Gamma\approx3$. We search for active galactic nuclei (AGNs) with such extreme X-ray variability in the literature and find that most of them are narrow-line Seyfert 1 galaxies and quasars with high accretion rates. The fraction of extremely X-ray variable objects among super-Eddington accreting AGNs is estimated to be $\approx 15-24\%$. We discuss two possible scenarios, disk reflection and partial covering absorption, to explain the extreme X-ray variability of SDSS J$075101.42+291419.1$. We propose a possible origin for the partial covering absorber, which is the thick inner accretion disk and its associated outflow in AGNs with high accretion rates.

We study the problem of EB-leakage that is associated with incomplete polarized CMB sky. In the blind case that assumes no additional information about the statistical properties and amplitudes of the signal from the missing sky region, we prove that the recycling method (Liu et al.~2018) gives the unique best estimate of the EB-leakage. Compared to the previous method, this method reduces the uncertainties in the BB power spectrum due to EB-leakage by more than one order of magnitude in the most interesting domain of multipoles, where $\ell$ is between $80$ and $200$. This work also provides a useful guideline for observational design of future CMB experiments.

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

We present a new public Python package, DarkHistory, for computing the effects of dark matter annihilation and decay on the temperature and ionization history of the early universe. DarkHistory simultaneously solves for the evolution of the free electron fraction and gas temperature, and for the cooling of annihilation/decay products and the secondary particles produced in the process. Consequently, we can self-consistently include the effects of both astrophysical and exotic sources of heating and ionization, and automatically take into account backreaction, where modifications to the ionization/temperature history in turn modify the energy-loss processes for injected particles. We present a number of worked examples, demonstrating how to use the code in a range of different configurations, in particular for arbitrary dark matter masses and annihilation/decay final states. Possible applications of DarkHistory include mapping out the effects of dark matter annihilation/decay on the global 21cm signal and the epoch of reionization, as well as the effects of exotic energy injections other than dark matter annihilation/decay. The code is available at https://github.com/hongwanliu/DarkHistory with documentation at https://darkhistory.readthedocs.io . Data files required to run the code can be downloaded at https://doi.org/10.7910/DVN/DUOUWA .

We describe PSR J1926-0652, a pulsar recently discovered with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Using sensitive single-pulse detections from FAST and long-term timing observations from the Parkes 64-m radio telescope, we probed phenomena on both long and short time scales. The FAST observations covered a wide frequency range from 270 to 800 MHz, enabling individual pulses to be studied in detail. The pulsar exhibits at least four profile components, short-term nulling lasting from 4 to 450 pulses, complex subpulse drifting behaviours and intermittency on scales of tens of minutes. While the average band spacing P3 is relatively constant across different bursts and components, significant variations in the separation of adjacent bands are seen, especially near the beginning and end of a burst. Band shapes and slopes are quite variable, especially for the trailing components and for the shorter bursts. We show that for each burst the last detectable pulse prior to emission ceasing has different properties compared to other pulses. These complexities pose challenges for the classic carousel-type models.

In CMB science, the simplest idea to remove a contaminated sky region is to multiply the sky map with a mask that is 0 for the contaminated region and 1 elsewhere, which is also called a top-hat masking. Although it is easy to use, such top-hat masking is known to suffer from various leakage problems. Therefore, we want to extend the top-hat masking to a series of semi-analytic functions called the P-filters. Most importantly, the P-filters can seamlessly realize the core idea of masking in CMB science, and, meanwhile, guarantee continuity up to the first derivative everywhere. The P-filters can significantly reduce many leakage problems without additional cost, including the leakages due to low-, high-, and band-pass filtering, and the E-to-E, B-to-B, B-to-E, and E-to-B leakages. The workings of the P-filter are illustrated by using the WMAP and Planck polarization sky maps. By comparison to the corresponding WMAP/Planck masks, we show that the P-filter performs much better than top-hat masking, and meanwhile, has the potential to supersede the principal idea of masking in CMB science. Compared to mask apodization, the P-filter is ``outward'', that tends to make proper use of the region that was marked as 0; whereas apodization is ``inward'', that always kills more signal in the region marked as 1.

In recent years, the disk populations in a number of young star-forming regions have been surveyed with ALMA. Understanding the disk properties and their correlation with those of the central star is critical to understand planet formation. In particular, a decrease of the average measured disk dust mass with the age of the region has been observed. We conducted high-sensitivity continuum ALMA observations of 43 Class II young stellar objects in CrA at 1.3 mm (230 GHz). The typical spatial resolution is 0.3". The continuum fluxes are used to estimate the dust masses of the disks, and a survival analysis is performed to estimate the average dust mass. We also obtained new VLT/X-Shooter spectra for 12 of the objects in our sample. 24 disks are detected, and stringent limits have been put on the average dust mass of the non-detections. Accounting for the upper limits, the average disk mass in CrA is $6\pm3\,\rm M_\oplus$, significantly lower than that of disks in other young (1-3 Myr) star forming regions (e.g. Lupus) and appears consistent with the 5-10 Myr old Upper Sco. The position of the stars in our sample on the HR diagram, however, seems to confirm that that CrA has age similar to Lupus. Neither external photoevaporation nor a lower than usual stellar mass distribution can explain the low disk masses. On the other hand, a low-mass disk population could be explained if the disks are small, which could happen if the parent cloud has a low temperature or intrinsic angular momentum, or if the the angular momentum of the cloud is removed by some physical mechanism such as magnetic braking. In order to fully explain and understand the dust mass distribution of protoplanetary disks and their evolution, it may also be necessary to take into consideration the initial conditions of star and disk formation process, which may vary from region to region, and affect planet formation.

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

The Five-hundred-meter Aperture Spherical radio Telescope(FAST) is expected to complete its commissioning in 2019. FAST will soon begin the Commensal Radio Astronomy FasT Survey(CRAFTS), a novel and unprecedented commensal drift scan survey of the entire sky visible from FAST. The goal of CRAFTS is to cover more than 20000 $deg^{2}$ and reach redshift up to about 0.35. We provide empirical measurements of the beam size and sensitivity of FAST across the 1.05 to 1.45 GHz frequency range of the FAST L-band Array of 19-beams(FLAN). Using a simulated HI-galaxy catalogue based on the HI Mass Function(HIMF), we estimate the number of galaxies that CRAFTS may detect. At redshifts below 0.35, over $6\, \times \, 10^{5}$ HI galaxies may be detected. Below the redshift of 0.07, the CRAFTS HIMF will be complete above a mass threshold of $10^{9.5}\,M_{\odot}$. FAST will be able to investigate the environmental and redshift dependence of the HIMF to an unprecedented depth, shedding light onto the missing baryon and missing satellite problems.

The Five-hundred-meter Aperture Spherical radio Telescope (FAST) was completed with its main structure installed on September 25, 2016, after which it entered the commissioning phase. This paper aims to introduce the commissioning progress of the FAST over the past two years. To improve its operational reliability and ensure effective observation time, FAST has been equipped with a real-time information system for the active reflector system and hierarchical commissioning scheme for the feed support system, which ultimately achieves safe operation of the two systems. For meeting the high-performance indices, a high-precision measurement system was set up based on the effective control methods that were implemented for the active reflector system and feed support system. Since the commissioning of the FAST, a low-frequency ultra-wideband receiver and 19-beam 1.05-1.45 GHz receiver have been mainly used. Telescope efficiency, pointing accuracy, and system noise temperature were completely tested and ultimately achieved the acceptance indices of the telescope. The FAST has been in the process of national acceptance preparations and has begun to search for pulsars. In the future, it will still strive to improve its capabilities and expand its application prospects.

To assist with the commissioning (Jiang et al. 2019) of the Five-hundred-meter Aperture Spherical radio Telescope (FAST), we performed a pulsar search, with the primary goal of developing and testing the pulsar data acquisition and processing pipelines. We tested and used three pipelines, two (P1 and P2 hereafter) searched for the periodic signature of pulsars whereas the other one was used to search for bright single pulses (P3 hereafter). A pulsar candidate was discovered in the observation on the 22nd August, 2017, and later confirmed by the Parkes radio telescope on the 10th September, 2017. This pulsar, named PSR J1900-0134, was the first pulsar discovered by FAST. The pulsar has a pulse period of 1.8 s and a dispersion measure (DM) of 188\,pc\,cm$^{-3}$.

Filaments are ubiquitous structures in molecular clouds and play an important role in the mass assembly of stars. We present results of dynamical stability analyses for filaments in the infrared dark cloud G14.225$-$0.506, where a delayed onset of massive star formation was reported in the two hubs at the convergence of multiple filaments of parsec length. Full-synthesis imaging is performed with the Atacama Large Millimeter/submillimeter Array (ALMA) to map the $\mathrm{N_2H^+} \; (1-0)$ emission in two hub-filament systems with a spatial resolution of $\sim 0.034 \; \mathrm{pc}$. Kinematics are derived from sophisticated spectral fitting algorithm that accounts for line blending, large optical depth, and multiple velocity components. We identify five velocity coherent filaments and derive their velocity gradients with principal component analysis. The mass accretion rates along the filaments are up to $10^{-4} \; \mathrm{M_\odot \, \mathrm{yr^{-1}}}$ and are significant enough to affect the hub dynamics within one free-fall time ($\sim 10^5 \; \mathrm{yr}$). The $\mathrm{N_2H^+}$ filaments are in equilibrium with virial parameter $\alpha_\mathrm{vir} \sim 1.2$. We compare $\alpha_\mathrm{vir}$ measured in the $\mathrm{N_2H^+}$ filaments, $\mathrm{NH_3}$ filaments, $870 \; \mu\mathrm{m}$ dense clumps, and $3 \; \mathrm{mm}$ dense cores. The decreasing trend in $\alpha_\mathrm{vir}$ with decreasing spatial scales persists, suggesting an increasingly important role of gravity at small scales. Meanwhile, $\alpha_\mathrm{vir}$ also decreases with decreasing non-thermal motions. In combination with the absence of high-mass protostars and massive cores, our results are consistent with the global hierarchical collapse scenario.

We present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other.

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of the CO ($J$=2--1) line emission from the protoplanetary disk around T-Tauri star SU Aurigae (hereafter SU Aur). Previous observations in optical and near infrared wavelengths find a unique structure in SU Aur. One of the highlights of the observational results is that an extended tail-like structure is associated with the disk, indicating mass transfer from or into the disk. Here we report the discovery of the counterpart of the tail-like structure in CO gas extending more than 1000 au long. Based on geometric and kinematic perspectives, both of the disk and the tail-like structure components physically connect to each other. Several theoretical studies predicted the observed tail-like structure via the following possible scenarios, 1) a gaseous stream from the molecular cloud remnant, 2) collision with a (sub)stellar intruder or a gaseous blob from the ambient cloud, and 3) ejection of a planetary or brown dwarf mass object due to gravitational instability via multi-body gravitational interaction. Since the tail-like structures associated with the SU Aur disk is a new example following RW Aurigae, some disks may experience the internal or external interaction and drastically lose mass during disk evolution.

We investigate star formation at very early evolutionary phases in five massive clouds in the inner 500 pc of the Galaxy, the Central Molecular Zone. Using interferometer observations of H$_2$O masers and ultra-compact H II regions, we find evidence of ongoing star formation embedded in cores of 0.2 pc scales and $\gtrsim$10$^5$ cm$^{-3}$ densities. Among the five clouds, Sgr C possesses a high (9%) fraction of gas mass in gravitationally bound and/or protostellar cores, and follows the dense ($\gtrsim$10$^4$ cm$^{-3}$) gas star formation relation that is extrapolated from nearby clouds. The other four clouds have less than 1% of their cloud masses in gravitationally bound and/or protostellar cores, and star formation rates 10 times lower than predicted by the dense gas star formation relation. At the spatial scale of these cores, the star formation efficiency is comparable to that in Galactic disk sources. We suggest that the overall inactive star formation in these Central Molecular Zone clouds could be because there is much less gas confined in gravitationally bound cores, which may be a result of the strong turbulence in this region and/or the very early evolutionary stage of the clouds when collapse has only recently started.

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

We present a spectral analysis of cool and cooling gas in 45 cool-core clusters and groups of galaxies obtained from Reflection Grating Spectrometer (RGS) XMM-$Newton$ observations. The high-resolution spectra show FeXVII emission in many clusters, which implies the existence of cooling flows. The cooling rates are measured between the bulk Intracluster Medium (ICM) temperature and 0.01 keV and are typically weak, operating at less than a few tens of $\rm M_{\odot}\rm yr^{-1}$ in clusters, and less than 1 $\rm M_{\odot}\rm yr^{-1}$ in groups of galaxies. They are 10-30% of the classical cooling rates in the absence of heating, which suggests that AGN feedback has a high level of efficiency. If cooling flows terminate at 0.7 keV in clusters, the associated cooling rates are higher, and have a typical value of a few to a few tens of $\rm M_{\odot}\rm yr^{-1}$. Since the soft X-ray emitting region, where the temperature $kT<1$ keV, is spatially associated with H$\alpha$ nebulosity, we examine the relation between the cooling rates above 0.7 keV and the H$\alpha$ nebulae. We find that the cooling rates have enough energy to power the total UV-optical luminosities, and are 5 to 50 times higher than the observed star formation rates for low luminosity objects. In 4 high luminosity clusters, the cooling rates above 0.7 keV are not sufficient and an inflow at a higher temperature is required. Further residual cooling below 0.7 keV indicates very low complete cooling rates in most clusters.

HL Tau is a Class I-II protostar embedded in an infalling and rotating envelope and possibly associated with a planet forming disk, and it is co-located in a 0.1 pc molecular cloud with two nearby young stellar objects. Our ALMA observations revealed two arc-like structures on a 1000 au scale connected to the disk, and their kinematics could not be explained with any conventional model of infalling and rotational motions. In this work, we investigate the nature of these arc-like structures connected to the HL Tau disk. We conducted new observations in the 13CO and C18O (3-2; 2-1) lines with JCMT and IRAM 30m, and obtained the ACA data with the 7-m array. With the single-dish, ACA, and ALMA data, we analyzed the gas motions on both 0.1 pc and 1000 au scales in the HL Tau region. We constructed new kinematical models of an infalling and rotating envelope with the consideration of relative motion between HL Tau and the envelope. By including the relative motion between HL Tau and its protostellar envelope, our kinematical model can explain the observed velocity features in the arc-like structures. The morphologies of the arc-like structures can also be explained with an asymmetric initial density distribution in our model envelope. In addition, our single-dish results support that HL Tau is located at the edge of a large-scale (0.1 pc) expanding shell driven by the wind or outflow from XZ Tau, as suggested in the literature. The estimated expanding velocity of the shell is comparable to the relative velocity between HL Tau and its envelope in our kinematical model. These results hints that the large-scale expanding motion likely impacts the protostellar envelope around HL Tau and affects its gas kinematics. We found that the mass infalling rate from the envelope onto the HL Tau disk can be decreased by a factor of two due to this impact by the large-scale expanding shell.

Fuzzy Dark Matter (FDM), motivated by string theory, has recently become a hot candidate for dark matter. The rest mass of FDM is believed to be $\sim 10^{-22}$eV and the corresponding de-Broglie wave length is $\sim 1$kpc. Therefore, the quantum effect of FDM plays an important role in structure formation. In order to study the cosmological structure formation in FDM model, several simulation techniques have been introduced. We review the current status and challenges in the cosmological simulation for the FDM model in this paper.

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

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

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.

We present combined observations of the NH3 (J,K) = (1,1) and (2,2) inversion transitions towards OMC1 in Orion A obtained by the Karl G. Jansky Very Large Array (VLA) and the 100 m Robert C. Byrd Green Bank Telescope (GBT). With an angular resolution of 6" (0.01 pc), these observations reveal with unprecedented detail the complex filamentary structure extending north of the active Orion BN/KL region in a field covering 6' x 7'. We find a 0.012 pc wide filament within OMC1, with an aspect ratio of ~37:1, that was missed in previous studies. Its orientation is directly compared to the relative orientation of the magnetic field from the James Clerk Maxwell Telescope BISTRO survey in Orion A. We find a small deviation of ~11 deg between the mean orientation of the filament and the magnetic field, suggesting that they are almost parallel to one another. The filament's column density is estimated to be 2-3 orders of magnitude larger than the filaments studied with Herschel and is possibly self-gravitating given the low values of turbulence found. We further produce maps of the gas kinematics by forward modeling the hyperfine structure of the NH3 (J,K) = (1,1) and (2,2) lines. The resulting distribution of velocity dispersions peaks at ~0.5 km/s, close to the subsonic regime of the gas. This value is about 0.2 km/s smaller than previously measured in single-dish observations of the same region, suggesting that higher angular and spectral resolution observations will identify even lower velocity dispersions that might reach the subsonic turbulence regime in dense gas filaments.

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.

Tidal disruption events (TDEs) in active galactic nuclei (AGNs) have been overlooked for a long time but tentatively been investigated recently. We report the discovery of a long-lasting luminous mid-infrared (mid-IR) flare in PS1-10adi, which is a newly-identified highly energetic transient event occurred in AGN. The IR luminosity of PS1-10adi, as well as other analogous events, are at least one order of magnitude higher than all known supernova, but can be well interpreted as the dust echoes of TDEs, whose ultra-high IR energy is reprocessed from the dusty torus around the black hole. The torus dust is sublimating during the early stage of the outburst and probably lead to the observed rapid emergence of Fe II lines. Moreover, the UV-optical rebrightening and contemporaneous X-ray onset after $\sim1500$ rest-frame days since the optical peak is also an intriguing feature of PS1-10adi, which could be attributed to the interaction between the high-velocity outflow and torus. We suggest that the luminous IR echo is a very typical phenomenon of TDEs in AGNs and may provide us an ideal opportunity to explore the torus properties.

We report new, $\sim$1000 AU spatial resolution observations of 225 GHz dust continuum emission towards the OB cluster-forming molecular clump G33.92+0.11. On parsec scales, this molecular clump presents a morphology with several arm-like dense gas structures surrounding the two central massive ($\gtrsim$100 $M_{\odot}$) cores. From the new, higher resolution observations, we identified 28 localized, spatially compact dust continuum emission sources, which may be candidates of young stellar objects. Only one of them is not embedded within known arm-like (or elongated) dense gas structures. The spatial separations of these compact sources can be very well explained by Jeans lengths. We found that G33.92+0.11 may be consistently described by a marginally centrifugally supported, Toomre unstable accretion flow which is approximately in a face-on projection. The arm-like overdensities are natural consequence of the Toomre instability, which can fragment to form young stellar objects in shorter time scales than the timescale of the global clump contraction. On our resolved spatial scales, there is not yet evidence that the fragmentation is halted by turbulence, magnetic field, or stellar feedback.

The filamentary infrared dark cloud L1287 is actively forming a dense cluster of low-mass YSOs at its inner $\sim$0.1 pc region. To help understand the origin of this low-mass YSO cluster, the present work aims at resolving the gas structures and kinematics. We have performed $\sim$1$"$ angular resolution ($\sim$930 AU) SMA observations at $\sim$1.3 mm. From a $\sim$2$"$ resolution 1.3 mm continuum image we identified six dense cores, namely SMA1-6 with masses in the range of $\sim0.4-4$ M$_\odot$. From a $\sim$1$"$ resolution 1.3 mm continuum image, we find a high fragmentation level, with 14 compact millimeter sources within 0.1 pc (two of them associated with the known accretion outburst YSOs RNO 1C and RNO 1B). The dense gas tracer DCN (3--2) traces well the dust continuum emission and shows a clear velocity gradient along the NW-SE direction centered at SMA3. There is another velocity gradient with opposite direction around the most luminous YSO IRAS 00338+6312. The fragmentation within 0.1 pc in L1287 is very high compared to other regions at the same spatial scales. The incoherent motions of dense gas flows are sometimes interpreted by being influenced by (proto)stellar feedback (e.g., outflows), which is not yet ruled out in this particular target source. The directions of the velocity gradients traced by DCN are approximately perpendicular to those of the dominant CO outflow(s). Therefore, we alternatively hypothesize that the velocity gradients revealed by DCN trace the convergence from the $\gtrsim$0.1 pc scales infalling motion towards the rotational motions around the more compact ($\sim0.02$ pc) sources. This global molecular gas converging flow may feed the formation of the dense low-mass YSO cluster. IRAS 00338+6312 is the most likely powering source of the dominant CO outflow. A compact blue-shifted outflow from RNO 1C is also identified.

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

Better understanding of Galactic foregrounds is one of the main obstacles to detection of primordial gravitational waves through measurement of the B mode in the polarized microwave sky. We generalize the method proposed in \cite{2011PhRvD..83h3003B} and decompose the polarization signals into the E and B families directly in the domain of the Stokes $Q$, $U$ parameters as $(Q,U)\equiv(Q_E, U_E)+(Q_B,U_B)$. This also enables an investigation of the morphology and the frequency dependence of these two families, which has been done in the WMAP K, Ka (tracing synchrotron emission) and Planck 2015 HFI maps (tracing thermal dust). The results reveal significant differences in spectra between the E and B families. The spectral index of the E family fluctuates less across the sky than that of the B family, and the same tendency occurs for the polarization angles of the dust and synchrotron channels. The new insight from WMAP and Planck data on the North Polar Spur and BICEP2 zones through our method clearly indicates that these zones are characterized by very low polarization intensity of the B family compared to the E family. We have detected global structure of the B family polarization angles at high Galactic latitudes which cannot be attributed to the cosmic microwave background or instrumental noise. However, we cannot exclude instrumental systematics as a partial contributor to these anomalies.

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of the dust continuum emission at 1.3 mm and 12CO J=2-1 line emission of the transitional disk around DM Tau. DM Tau's disk is thought to possess a dust-free inner cavity inside a few au, from the absence of near-infrared excess on its spectral energy distribution (SED). Previous submillimeter observations were, however, unable to detect the cavity; instead, a dust ring ~20 au in radius was seen. The excellent angular resolution achieved in the new ALMA observations, 43 x 31 mas, allows discovery of a 4 au radius inner dust ring, confirming previous SED modeling results. This inner ring is symmetric in continuum emission, but asymmetric in 12CO emission. The known (outer) dust ring at ~20 au is recovered and shows azimuthal asymmetry with a strong-weak side contrast of ~1.3. The gap between these two rings is depleted by a factor of ~40 in dust emission relative to the outer ring. An extended outer dust disk is revealed, separated from the outer ring by another gap. The location of the inner ring is comparable to that of the main asteroid belt in the solar system. As a disk with a "proto-asteroid belt," the DM Tau system offers valuable clues to disk evolution and planet formation in the terrestrial planet forming region.

The metallicity and its relationship with other galactic properties is a fundamental probe of the evolution of galaxies. In this work, we select about 750,000 star-forming spatial pixels from 1122 blue galaxies in the MaNGA survey to investigate the global stellar mass - local stellar mass surface density - gas-phase metallicity ($M_*$ - $\Sigma_*$ - $Z$ ) relation. At a fixed $M_*$, the metallicity increases steeply with increasing $\Sigma_*$. Similarly, at a fixed $\Sigma_*$, the metallicity increases strongly with increasing $M_*$ at low mass end, while this trend becomes less obvious at high mass end. We find the metallicity to be more strongly correlated to $\Sigma_*$ than to $M_*$. Furthermore, we construct a tight (0.07 dex scatter) $M_*$ - $\Sigma_*$ - $Z$ relation, which reduces the scatter in the $\Sigma_*$ - $Z$ relation by about 30$\%$ for galaxies with $7.8 < {\rm log}(M_*/M_\odot) < 11.0$, while the reduction of scatter is much weaker for high-mass galaxies. This result suggests that, especially for low-mass galaxies, the $M_*$ - $\Sigma_*$ - $Z$ relation is largely more fundamental than the $M_*$ - $Z$ and $\Sigma_*$ - $Z$ relations, meaning that both $M_*$ and $\Sigma_*$ play important roles in shaping the local metallicity. We also find that the local metallicity is probably independent on the local star formation rate surface density at a fixed $M_*$ and $\Sigma_*$. Our results are consistent with the scenario that the local metallicities in galaxies are shaped by the combination of the local stars formed in the history and the metal loss caused by galactic winds.

More than one hundred galaxy-scale strong gravitational lens systems have been found by searching for the emission lines coming from galaxies with redshifts higher than the lens galaxies. Based on this spectroscopic-selection method, we introduce the deep Residual Networks (ResNet, a kind of deep Convolutional Neural Networks) to search for the galaxy-Ly$\alpha$ emitter (LAE) lens candidates by recognizing the Ly$\alpha$ emission lines coming from high redshift galaxies ($2 < z < 3$) in the spectra of early-type galaxies (ETGs) at middle redshift ($z\sim 0.5$). The spectra of the ETGs come from the Data Release 12 (DR12) of the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey \uppercase\expandafter{\romannumeral3} (SDSS-\uppercase\expandafter{\romannumeral3}). In this paper, we first build a 28 layers ResNet model, and then artificially synthesize 150,000 training spectra, including 140,000 spectra without Ly$\alpha$ lines and 10,000 ones with Ly$\alpha$ lines, to train the networks. After 20 training epochs, we obtain a near-perfect test accuracy at 0.9954. The corresponding loss is 0.0028 and the completeness is 93.6\%. We finally apply our ResNet model to our predictive data with 174 known lens candidates. We obtain 1232 hits including 161 of the 174 known candidates (92.5\% discovery rate). Apart from the hits found in other works, our ResNet model also find 536 new hits. We then perform several subsequent selections on these 536 hits and present 5 most believable lens candidates.

In this paper we explore the effects of self-obscuration in protostellar disks with a radially decreasing temperature gradient and a colder midplane. We are motivated by recent reports of resolved dark lanes (`hamburgers') and (sub)mm spectral indices systematically below the ISM value for optically thin dust $\alpha_{\rm ISM} =3.7$. We explore several model grids, scaling disk mass and varying inclination angle $i$ and observing frequency $\nu$ from the VLA Ka band ($\sim 37$ GHz) to ALMA Band 8 ($\sim 405$ GHz). We also consider the effects of decreasing the index of the (sub-)mm dust opacity power law $\beta$ from 1.7 to 1. We find that a distribution of disk masses in the range $M_{\rm disk} = 0.01-2~M_\odot$ is needed to reproduce the observed distribution of spectral indices, and that assuming a fixed $\beta =1.7$ gives better results than $\beta=1$. A wide distribution of disk masses is also needed to produce some cases with $\alpha <2$, as reported for some sources in the literature. Such extremely low spectral indices arise naturally when the selected observing frequencies sample the appropriate change in the temperature structure of the optically thick model disk. Our results show that protostellar disk masses could often be underestimated by $> \times10$, and are consistent with recent hydrodynamical simulations. Although we do not rule out the possibility of some grain growth occurring within the short protostellar timescales, we conclude that self-obscuration needs to be taken into account.

Context: Currently, detection of the primordial gravitational waves by the B-mode of Cosmic Microwave Background (CMB) is primarily limited by our knowledge of the polarized microwave foreground emissions. Thus improvements of the foreground analysis are necessary. As revealed in~\cite{2018arXiv180410382L}, the E-mode and B-mode of the polarized foreground have noticeable different properties, both in morphology and frequency spectrum, suggesting that they arise from different physical processes, and need to be studied separately. Aims: I will study the polarized emission from Galactic loops, especially Loop I, and mainly focus on the following issues: Does it contribute predominantly to the E-mode or B-mode? In which frequency bands and in which sky regions can it be identified? Methods: Based on a well known result about the magnetic field alignment in supernova explosions, a theoretical expectation is established that the loop polarizations should be predominantly E-mode. In particular, the expected polarization angles of Loop I are compared with those from the real microwave band data of WMAP and Planck. Results and conclusions: The comparison between model and data shows remarkable consistency between data and expectation at all bands and for a large area of the sky. This result suggests that the polarized emission of Galactic Loop I is a major polarized component in all microwave bands from 23 to 353 GHz, and a considerable part of the polarized foreground is likely originated from a local bubble associated with Loop I, instead of the far more distant Galactic emission. The result also provides a possible way to explain the reported E-to-B excess~\citep{2016A&A...586A.133P} by contribution of the loops. Finally, this work may also provide the first geometrical evidence that the Earth was hit by a supernova explosion.

It is commonly thought that the final product of gravitational collapse is a black hole. Nevertheless, theoretical studies have not yet provided a final answer to the question whether black holes are the only possible outcome or whether naked singularities are also allowed. Observational tests may thus represent a complementary approach. In the present paper, we consider the Janis-Newman-Winicour metric, which describes a rotating source with a surface-like naked singularity. We calculate iron line shapes in the reflection spectrum of a putative disk around a Janis-Newman-Winicour singularity and we compare our results with the iron line shapes expected in the spectrum of a Kerr black hole. While it is difficult to distinguish the two spacetimes from the iron line shape in general, it seems that Janis-Newman-Winicour singularities cannot mimic fast-rotating Kerr black holes observed at a low or moderate inclination angle. Our analysis thus suggests that available observations of specific sources can already constrain the possible existence of Janis-Newman-Winicour singularities in the Universe.

We present near-IR imaging polarimetry of five classical FU Ori-type objects (FU Ori, V1057 Cyg, V1515 Cyg, V1735 Cyg, Z CMa) with a $\sim$0\farcs1 resolution observed using HiCIAO+AO188 at Subaru Telescope. We observed scattered light associated with circumstellar dust around four of them (i.e., all but V1515 Cyg). Their polarized intensity distribution shows a variety of morphologies with arms, tails or streams, spikes and fragmented distributions, many of which were reported in our previous paper. The morphologies of these reflection nebulae significantly differ from many other normal young stellar objects (Class I-II objects). These structures are attributed to gravitationally unstable disks, trails of clump ejections, dust blown by a wind or a jet, and a stellar companion. We can consistently explain our results with the scenario that their accretion outbursts (FUor outbursts) are triggered by gravitationally fragmenting disks, and with the hypothesis that many low-mass young stellar objects experience such outbursts.

During the evolution of proto-planetary disc, photo-evaporations of both central and external stars play important roles. Considering the complicated radiation surroundings in the clusters, where the star formed, the proto-planetary discs survive in different lifetimes due to flyby events. In this paper, we mainly focus on the disc around a T Tauri star, which encounters with another main-sequence star with different temperatures in hyperbolic orbits with different peri-center distances, eccentricities and inclinations. We find the criterion for gap-opening due to photo-evaporation of central star after the flyby event. A gap is opened in the late stage of gas disc, and induce that the gap only influence the planet formation and migration limitedly. If the flyby orbit has a moderate value of peri-center distance, which weakly depends on the eccentricity and inclination, the external photo-evaporation lead to a maximum mass loss during the flyby event. Flyby stars in orbits with smaller eccentricities or larger inclinations induce larger mass loss. Adopting a simple multiple flyby models, we conclude that in open clusters, gas discs usually survive in typical lifetimes between 1 and 10 Myr, except there are many massive stars in dense open clusters. In globular clusters, discs disperse very quickly and hardly produce the gas giant planets. The fast-depleted discs are probably responsible for the null detection of giant planets in globular clusters.

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

With recent Lyman-alpha forest data from BOSS and XQ-100, some studies suggested that the lower mass limit on the fuzzy dark matter (FDM) particles is lifted up to $10^{-21}\,\mathrm{eV}$. However, such a limit was obtained by $\Lambda$CDM simulations with the FDM initial condition and the quantum pressure (QP) was not taken into account which could have generated non-trivial effects in large scales structures. We investigate the QP effects in cosmological simulations systematically, and find that the QP leads to further suppression of the matter power spectrum at small scales, as well as the halo mass function in the low mass end. We estimate the suppressing effect of QP in the 1D flux power spectrum of Lyman-alpha forest and compare it with data from BOSS and XQ-100. The rough uncertainties of thermal gas properties in the flux power spectrum model calculation were discussed. We conclude that more systematic studies, especially with QP taken into account, are necessary to constrain FDM particle mass using Lyman-alpha forest.

We report Submillimeter Array molecular line observations in the 345 GHz band of five massive dense cores, Cyg-N38, Cyg-N43, Cyg-N48, Cyg-N51, and Cyg-N53 in the DR21 filament. The molecular line data reveal several dynamical features of the cores: (1) prominent outflows in all cores seen in the CO and SiO lines, (2) significant velocity gradients in Cyg-N43 and Cyg-N48 seen in the H13CN and H13CO+ lines suggesting 0.1-pc-scale rotational motions, and (3) possible infalls in Cyg-N48 found in the SiO and SO lines. Comparing the molecular line data and our dust polarization data in Ching et al. (2017), we find that the gradients of line-of-sight velocities appear to be randomly oriented relative to the plane-of-sky magnetic fields. Our simulations suggest that this random alignment implies parallel or random alignment between the velocity gradients and magnetic fields in the three dimensional space. The linewidths of H13CN emission are consistently wider than those of H13CO+ emission in the 3''-10'' detectable scales, which can be explained by the existence of ambipolar diffusion with maximum plane-of-the-sky magnetic field strengths of 1.9 mG and 5.1 mG in Cyg-N38 and Cyg-N48, respectively. Our results suggest that the gas dynamics may distort the magnetic fields of the cores of into complex structures and ambipolar diffusion could be important in dissipating the magnetic energies of the cores.

Active galactic nuclei (AGN) release a huge amount of energy into the intracluster medium (ICM) with the consequence of offsetting cooling and star formation (AGN feedback) in the centers of cool core clusters. The Phoenix cluster is among the most massive clusters of galaxies known in the Universe. It hosts a powerful starburst of several hundreds of Solar masses per year and a large amount of molecular gas in the center. In this work we use the high-resolution Reflection Grating Spectrometer (RGS) on board XMM-Newton to study the X-ray emitting cool gas in the Phoenix cluster and heating-cooling balance. We detect for the first time evidence of O VIII and Fe XXI-XXII emission lines, the latter demonstrating the presence of gas below 2 keV. We find a cooling rate of 350 (-200,+250) Msun/year below 2 keV (at the 90% confidence level), which is consistent with the star formation rate in this object. This cooling rate is high enough to produce the molecular gas found in the filaments via instabilities during the buoyant rising time. The line broadening indicates that the turbulence (~ 300 km/s or less) is below the level required to produce and propagate the heat throughout the cool core. This provides a natural explanation to the coexistence of large amounts of cool gas, star formation and a powerful AGN in the core. The AGN activity may be either at a young stage or in a different feedback mode, due to a high accretion rate.

The spectral evolution and spectral lag behavior of 92 bright pulses from 84 gamma-ray bursts (GRBs) observed by the Fermi GBM telescope are studied. These pulses can be classified into hard-to-soft pulses (H2S, 64/92), H2S-dominated-tracking pulses (21/92), and other tracking pulses (7/92). We focus on the relationship between spectral evolution and spectral lags of H2S and H2S-dominated-tracking pulses. %in hard-to-soft pulses (H2S, 64/92) and H2S-dominating-tracking (21/92) pulses. The main trend of spectral evolution (lag behavior) is estimated with $\log E_p\propto k_E\log(t+t_0)$ (${\hat{\tau}} \propto k_{\hat{\tau}}\log E$), where $E_p$ is the peak photon energy in the radiation spectrum, $t+t_0$ is the observer time relative to the beginning of pulse $-t_0$, and ${\hat{\tau}}$ is the spectral lag of photons with energy $E$ with respect to the energy band $8$-$25$ keV. For H2S and H2S-dominated-tracking pulses, a weak correlation between $k_{{\hat{\tau}}}/W$ and $k_E$ is found, where $W$ is the pulse width. We also study the spectral lag behavior with peak time $t_{\rm p_E}$ of pulses for 30 well-shaped pulses and estimate the main trend of the spectral lag behavior with $\log t_{\rm p_E}\propto k_{t_p}\log E$. It is found that $k_{t_p}$ is correlated with $k_E$. We perform simulations under a phenomenological model of spectral evolution, and find that these correlations are reproduced. We then conclude that spectral lags are closely related to spectral evolution within the pulse. The most natural explanation of these observations is that the emission is from the electrons in the same fluid unit at an emission site moving away from the central engine, as expected in the models invoking magnetic dissipation in a moderately-high-$\sigma$ outflow.

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

In this contribution we explore the new science that a Next Generation Very Large Array (ngVLA) would be able to perform on the topics of the youngest HII regions and (proto)stellar jets. Free-free continuum and radio recombination line (RRL) emission are often the only way of peering into the dense envelopes surrounding (proto)stars of all masses, and trace their initial feedback in the form of `radio jets', `hypercompact HII regions', or photoevaporating, partially-ionized flows. Properly disentangling free-free from dust emission is also mandatory in studies of protoplanetary and accretion disks. Current VLA research has reached an impasse in which a population of faint ionized radio sources, probably corresponding to the above mentioned objects, is detected, but their nature is mostly unknown. The ngVLA would allow us to resolve the density structure and kinematics of such sources, revolutionizing our knowledge of star formation across the entire stellar-mass spectrum.

Using the Wald formalism, we investigate the thermodynamics of charged black holes in D-dimensional stationary spacetimes with or without rotations in Einstein-\ae ther-Maxwell theory. In particular, assuming the existence of a scaling symmetry of the action, we obtain the Smarr integral formula, which can be applied to both Killing and universal horizons. When restricted to 4-dimensional spherically symmetric spacetimes, previous results obtained by a different method are re-derived.

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

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

Einstein's theory of general relativity was proposed over 100 years ago and has successfully passed a large number of observational tests in the weak field regime. However, the strong field regime is largely unexplored, and there are many modified and alternative theories that have the same predictions as Einstein's gravity for weak fields and present deviations when gravity becomes strong. RELXILL_NK is the first relativistic reflection model for probing the spacetime metric in the vicinity of astrophysical black holes and testing Einstein's gravity in the strong field regime. Here we present our current constraints on possible deviations from Einstein's gravity obtained from the black holes in 1H0707-495, Ark 564, GX 339-4, and GS 1354-645.

We report on the identification of a new {\gamma}-ray-emitting narrow-line Seyfert 1 (NLS1) galaxy, SDSS J211852.96-073227.5 (hereinafter J2118-0732). The galaxy, at a redshift of 0.26, is associated with a radio source of flat/inverted spectrum at high radio frequencies. The analysis of its optical spectrum obtained in the Sloan Digital Sky Survey (SDSS) revealed a small linewidth of the broad component of the H{\beta} line (full width at half-maximum = 1585 km s^{-1}), making it a radio-loud NLS1 galaxy -- an intriguing class of active galactic nuclei with exceptional multiwavelength properties. A new {\gamma}-ray source centred at J2118-0732 was sporadically detected during 2009--2013 in form of flares by the Fermi-LAT. Our XMM-Newton observations revealed a flat X-ray spectrum described by a simple power law, and a flux variation by a factor of ~2.5 in five months. The source also shows intraday variability in the infrared band. Its broad-band spectral energy distribution can be modelled by emission from a simple one-zone leptonic jet model, and the flux drop from infrared to X-rays in five months can be explained by changes of the jet parameters, though the exact values may be subject to relatively large uncertainties. With the NLS1-blazar composite nucleus, the clear detection of the host galaxy, and the synchronous variations in the multiwavelength fluxes, J2118-0732 provides a new perspective on the formation and evolution of relativistic jets under the regime of relatively small black hole masses and high accretion rates.

Spatially resolved structures in protoplanetary disks hint at unseen planets. Previous imaging observations of the transitional disk around MWC 758 revealed an inner cavity, a ring-like outer disk, emission clumps, and spiral arms, all possibly generated by companions. We present ALMA dust continuum observations of MWC 758 at 0.87 millimeter (mm) wavelength with 43$\times$39 mas angular resolution (6.9$\times$6.2 AU) and 20 $\mu$Jy beam$^{-1}$ rms. The central sub-mm emission cavity is revealed to be eccentric; once deprojected, its outer edge can be well-fitted by an ellipse with an eccentricity of 0.1 and one focus on the star. The broad ring-like outer disk is resolved into three narrow rings with two gaps in between. The outer two rings tentatively show the same eccentricity and orientation as the innermost ring bounding the inner cavity. The two previously known dust emission clumps are resolved in both the radial and azimuthal directions, with radial widths equal to $\sim$4$\times$ the local scale height. Only one of the two spiral arms previously imaged in near-infrared (NIR) scattered light is revealed in ALMA dust emission, at a slightly larger stellocentric distance owing to projection effects. We also submit evidence of disk truncation at $\sim$100 AU based on comparing NIR imaging observations with models. The spirals, the north clump, and the truncated disk edge are all broadly consistent with the presence of one companion exterior to the spirals at roughly 100 AU.

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

Based on the 850 $\mu$m dust continuum data from SCUBA-2 at James Clerk Maxwell Telescope (JCMT), we compare overall properties of Planck Galactic Cold Clumps (PGCCs) in the $\lambda$ Orionis cloud to those of PGCCs in the Orion A and B clouds. The Orion A and B clouds are well known active star-forming regions, while the $\lambda$ Orionis cloud has a different environment as a consequence of the interaction with a prominent OB association and a giant Hii region. PGCCs in the $\lambda$ Orionis cloud have higher dust temperatures ($Td=16.13\pm0.15$ K) and lower values of dust emissivity spectral index ($ \beta=1.65\pm0.02$) than PGCCs in the Orion A (Td=13.79$\pm 0.21$K, $\beta=2.07\pm0.03$) and Orion B ($Td=13.82\pm0.19$K, $\beta=1.96\pm0.02$) clouds. We find 119 sub-structures within the 40 detected PGCCs and identify them as cores. Of total 119 cores, 15 cores are discovered in the $\lambda$ Orionis cloud, while 74 and 30 cores are found in the Orion A and B clouds, respectively. The cores in the $\lambda$ Orionis cloud show much lower mean values of size R=0.08 pc, column density N(H2)=$(9.5\pm1.2) \times 10^{22}$ cm$^{-2}$, number density n(H2)=$(2.9 \pm 0.4)\times10^{5}$ cm$^{-3}$, and mass $M_{core}$=$1.0\pm0.3$ M$_{\odot}$ compared to the cores in the Orion A (R=0.11pc, $N(H2)=(2.3\pm0.3) \times 10^{23}$ cm$^{-2}$, n(H2)=$(3.8\pm0.5) \times 10^{5}$cm$^{-3}$, and $M_{core}$=$2.4 \pm 0.3$ M$_{\odot}$) and Orion B (R=0.16pc, N(H2)=$(3.8 \pm 0.4) \times 10^{23}$cm$^{-2}$, n(H2)=$(15.6\pm1.8)\times10^{5}$ cm$^{-3}$, and $M_{core}$= $2.7\pm0.3$ M$_{\odot}$) clouds. These core properties in the $\lambda$ Orionis cloud can be attributed to the photodissociation and external heating by the nearby Hii region, which may prevent the PGCCs from forming gravitationally bound structures and eventually disperse them. These results support the idea of negative stellar feedback on core formation.

Based on the 850 $\mu$m dust continuum data from SCUBA-2 at James Clerk Maxwell Telescope (JCMT), we compare overall properties of Planck Galactic Cold Clumps (PGCCs) in the $\lambda$ Orionis cloud to those of PGCCs in the Orion A and B clouds. The Orion A and B clouds are well known active star-forming regions, while the $\lambda$ Orionis cloud has a different environment as a consequence of the interaction with a prominent OB association and a giant Hii region. PGCCs in the $\lambda$ Orionis cloud have higher dust temperatures ($Td=16.13\pm0.15$ K) and lower values of dust emissivity spectral index ($ \beta=1.65\pm0.02$) than PGCCs in the Orion A (Td=13.79$\pm 0.21$K, $\beta=2.07\pm0.03$) and Orion B ($Td=13.82\pm0.19$K, $\beta=1.96\pm0.02$) clouds. We find 119 sub-structures within the 40 detected PGCCs and identify them as cores. Of total 119 cores, 15 cores are discovered in the $\lambda$ Orionis cloud, while 74 and 30 cores are found in the Orion A and B clouds, respectively. The cores in the $\lambda$ Orionis cloud show much lower mean values of size R=0.08 pc, column density N(H2)=$(9.5\pm1.2) \times 10^{22}$ cm$^{-2}$, number density n(H2)=$(2.9 \pm 0.4)\times10^{5}$ cm$^{-3}$, and mass $M_{core}$=$1.0\pm0.3$ M$_{\odot}$ compared to the cores in the Orion A (R=0.11pc, $N(H2)=(2.3\pm0.3) \times 10^{23}$ cm$^{-2}$, n(H2)=$(3.8\pm0.5) \times 10^{5}$cm$^{-3}$, and $M_{core}$=$2.4 \pm 0.3$ M$_{\odot}$) and Orion B (R=0.16pc, N(H2)=$(3.8 \pm 0.4) \times 10^{23}$cm$^{-2}$, n(H2)=$(15.6\pm1.8)\times10^{5}$ cm$^{-3}$, and $M_{core}$= $2.7\pm0.3$ M$_{\odot}$) clouds. These core properties in the $\lambda$ Orionis cloud can be attributed to the photodissociation and external heating by the nearby Hii region, which may prevent the PGCCs from forming gravitationally bound structures and eventually disperse them. These results support the idea of negative stellar feedback on core formation.

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

We investigate the global properties of the straight and isolated filamentary cloud G350.54+0.69 using Herschel continuum and APEX molecular line data. The overall straight morphology is similar to two other well studied nearby filaments (Musca and Taurus-B211/3) while the isolated nature of G350.54+0.69 appears similar to Musca. G350.54+0.69 is composed of two distinct filaments with a length ~5.9 pc for G350.5-N (~2.3 pc for G350.5-S), a total mass of ~810 $M_{\odot}$ (~ 110 $M_{\odot}$), and a mean temperature of ~ 18.2 K (~17.7 K). We identify 9 dense and gravitationally bound cores in the whole cloud G350.54+0.69. The separations between cores and the line mass of the whole cloud appear to follow the predictions of the "sausage" instability theory, which suggests that G350.54+0.69 could have undergone radial collapse and fragmentation. The presence of young protostars is consistent with this hypothesis. The line masses of the two filaments (~120 $M_{\odot}$ pc$^{-1}$ for G350.5-N, and ~45 $M_{\odot}$ pc$^{-1}$ for G350.5-S), mass-size distributions of the dense cores, and low-mass protostars collectively suggest that G350.54+0.69 is a site of ongoing low-mass star formation. Based on the above evidence, we place G350.54+0.69 in an intermediate evolutionary state between Musca and Taurus-B211/3. We suggest that investigations into straight (and isolated) versus those distributed inside molecular clouds may provide important clues into filament formation and evolution.

We performed the new JVLA full polarization observations at 40-48 GHz (6.3-7.5 mm) towards the nearby ($d$ $=$147$\pm$3.4 pc) Class 0 YSO IRAS 16293-2422, and compare with the previous SMA observations reported by Rao et al. (2009; 2014). We observed the quasar J1407+2827 which is weakly polarized and can be used as a leakage term calibrator for $<$9 GHz observations, to gauge the potential residual polarization leakage after calibration. We did not detect Stokes Q, U, and V intensities from the observations of J1407+2827, and constrain (3-$\sigma$) the residual polarization leakage after calibration to be $\lesssim$0.3\%. We detect linear polarization from one of the two binary components of our target source, IRAS\,16293-2422\,B. The derived polarization position angles from our observations are in excellent agreement with those detected from the previous observations of the SMA, implying that on the spatial scale we are probing ($\sim$50-1000 au), the physical mechanisms for polarizing the continuum emission do not vary significantly over the wavelength range of $\sim$0.88-7.5 mm. We hypothesize that the observed polarization position angles trace the magnetic field which converges from large scale to an approximately face-on rotating accretion flow. In this scenario, magnetic field is predominantly poloidal on $>$100 au scales, and becomes toroidal on smaller scales. However, this interpretation remains uncertain due to the high dust optical depths at the central region of IRAS\,16293-2422\,B and the uncertain temperature profile. We suggest that dust polarization at wavelengths comparable or longer than 7\,mm may still trace interstellar magnetic field. Future sensitive observations of dust polarization in the fully optically thin regime will have paramount importance for unambiguously resolving the magnetic field configuration.

A new sample of 204 low-mass black holes (LMBHs) in active galactic nuclei (AGNs) is presented with black hole masses in the range of (1-20) * 10^5 M_sun. The AGNs are selected from a systematic search among galaxies in the Seventh Data Release (DR 7) of the Sloan Digital Sky Survey (SDSS), and careful analyses of their optical spectra and precise measurement of spectral parameters. Combining them with our previous sample selected from the SDSS DR 4 makes it the largest LMBH sample so far, totaling over 500 objects. Some of the statistical properties of the combined LMBH AGN sample are briefly discussed, in the context of exploring the low-mass end of the AGN population. Their X-ray luminosities follow the extension of the previously known correlation with the [O III] luminosity. The effective optical-to-X-ray spectral indices \alpha_OX, albeit with a large scatter, are broadly consistent with the extension of the relation with the near-UV luminosity L_2500\AA. Interestingly, a correlation of \alpha_OX with black hole mass is also found in the sense that \alpha_OX is statistically flatter (stronger X-ray relative to optical) for lower black hole mass. Only 26 objects, mostly radio loud, were detected in radio at 20 cm in the FIRST survey, giving a radio loud fraction of 4%. The host galaxies of LMBHs have stellar masses in the range of 10^8.8-10^12.4 M_sun and optical colors typical of Sbc spirals. They are dominated by young stellar populations that seem to have undergone a continuous star formation history.

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

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

We present our optical photometric and spectroscopical observations of GRB 140629A. A redshift of $z=2.275\pm0.043$ is measured through the metal absorption lines in our spectroscopic data. Using our photometric data and multiple observational data observed with other telescopes, we show that its optical light curve is well interpreted with the standard forward shock models in the thin shell case. Its optical-X-ray afterglow spectrum is jointly fitted with a single power-law function, yielding a photon index of $-1.90\pm 0.05$. The optical extinction and neutral hydrogen absorption of the GRB host galaxy are negligible. The fit to the light curve with the standard models shows that the ambient density is $60\pm 9$ cm$^{-3}$ and the GRB radiating efficiency is as low as $\sim 0.24\%$, likely indicating a baryonic-dominated ejecta of this GRB. This burst agrees well with the $L_{\rm p, iso}-E_p^{'}-\Gamma_0$ relation, but confidently violates those empirical relations involving geometric corrections (or jet break time). This gives rise to an issue of possible selection effect on these relations since the jet opening angle of this GRB is extremely narrow (0.04 rad).

We present combined $\approx 14-37~\rm ks$ Chandra observations of seven $z = 1.6-2.7$ broad absorption line (BAL) quasars selected from the Large Bright Quasar Survey (LBQS). These seven objects are high-ionization BAL (HiBAL) quasars, and they were undetected in the Chandra hard band ($2-8$ keV) in previous observations. The stacking analyses of previous Chandra observations suggested that these seven objects likely contain some candidates for intrinsically X-ray weak BAL quasars. With the new Chandra observations, six targets are detected. We calculate their effective power-law photon indices and hard-band flux weakness, and find that two objects, LBQS $1203+1530$ and LBQS $1442-0011$, show soft/steep spectral shapes ($\Gamma_{\rm eff}= 2.2^{+0.9}_{-0.9}$ and $1.9_{-0.8}^{+0.9}$) and significant X-ray weakness in the hard band (by factors of $\approx$ 15 and 12). We conclude that the two HiBAL quasars are good candidates for intrinsically X-ray weak BAL quasars. The mid-infrared-to-UV spectral energy distributions (SEDs) of the two candidates are consistent with those of typical quasars. We constrain the fraction of intrinsically X-ray weak AGNs among HiBAL quasars to be $\approx 7-10\%$ ($2/29-3/29$), and we estimate it is $\approx 6- 23\%$ ($2/35-8/35$) among the general BAL quasar population. Such a fraction is considerably larger than the fraction among non-BAL quasars, and we suggest that intrinsically X-ray weak quasars are preferentially observed as BAL quasars. Intrinsically X-ray weak AGNs likely comprise a small minority of the luminous type 1 AGN population, and they should not affect significantly the completeness of these AGNs found in deep X-ray surveys.

Measurements of the temperature of the baryons at the end of the cosmic dark ages can potentially set very precise constraints on energy injection from exotic sources, such as annihilation or decay of the dark matter. However, additional effects that lower the gas temperature can substantially weaken the expected constraints on exotic energy injection, whereas additional radiation backgrounds can conceal the effect of an increased gas temperature in measurements of the 21-cm hyperfine transition of neutral hydrogen. Motivated in part by recent claims of a detection of 21-cm absorption from a redshift of 17 by the EDGES experiment, we derive the constraints on dark matter annihilation and decay that can be placed in the presence of extra radiation backgrounds or effects that modify the gas temperature, such as dark matter-baryon scattering and early baryon-photon decoupling. We find that if the EDGES observation is confirmed, then constraints on light dark matter decaying or annihilating to electrons will in most scenarios be stronger than existing state-of-the-art limits from the cosmic microwave background, potentially by several orders of magnitude. More generally, our results allow mapping any future measurement of the global 21-cm signal into constraints on dark matter annihilation and decay, within the broad range of scenarios we consider.

We present ALMA 0.87 mm continuum, HCO+ J=4--3 emission line, and CO J=3--2 emission line data of the disk of material around the young, Sun-like star PDS 70. These data reveal the existence of a possible two component transitional disk system with a radial dust gap of 0."2 +/- 0."05, an azimuthal gap in the HCO+ J=4--3 moment zero map, as well as two bridge-like features in the gas data. Interestingly these features in the gas disk have no analogue in the dust disk making them of particular interest. We modeled the dust disk using the Monte Carlo radiative transfer code HOCHUNK3D (Whitney et al. 2013) using a two disk components. We find that there is a radial gap that extends from 15-60 au in all grain sizes which differs from previous work.

We investigate the spatially-resolved star formation main sequence in star-forming galaxies (SFGs) using Integral Field Spectroscopic (IFS) observations from the Mapping Nearby Galaxies at the Apache Point Observatory (MaNGA) survey. We demonstrate that the correlation between the stellar mass surface density ($\Sigma_*$) and star formation rate surface density ($\Sigma_{\mathrm{SFR}}$) holds down to sub-galactic scale, leading to the Sub-Galactic Main Sequence (SGMS). By dividing galaxies into two populations based on their recent mass assembly modes, we find the resolved main sequence in galaxies with 'outside-in' mode is steeper than that in galaxies with 'inside-out' mode. This is also confirmed on a galaxy-by-galaxy level, where we find the distributions of SGMS slopes for individual galaxies are clearly separated for the two populations. When normalizing and stacking the SGMS of individual galaxies on one panel for the two populations, we find the inner regions of galaxies with 'inside-out' mode statistically exhibit a suppression in star formation, with a less significant trend in the outer regions of galaxies with 'outside-in' mode. In contrast, the inner regions of galaxies with 'outside-in' mode and the outer regions of galaxies with 'inside-out' mode follow a slightly sub-linear scaling relation with a slope $\sim$0.9, which is in good agreement with previous findings, suggesting that they are experiencing a universal regulation without influences of additional physical processes.

We study the degeneracy of theoretical gravitational waveforms for binary black hole mergers using an aligned-spin effective-one-body model. After appropriate truncation, bandpassing, and matching, we identify regions in the mass--spin parameter space containing waveforms similar to the template proposed for GW150914, with masses $m_1 = 36^{+5}_{-4} M_\odot$ and $m_2 = 29^{+4}_{-4} M_\odot$, using the cross-correlation coefficient as a measure of the similarity between waveforms. Remarkably high cross-correlations are found across broad regions of parameter space. The associated uncertanties exceed these from LIGO's Bayesian analysis considerably. We have shown that waveforms with greatly increased masses, such as $m_1 = 70 M_\odot$ and $m_2 = 35 M_\odot$, and strong anti-aligned spins ($\chi_1=0.95$ and $\chi_2=-0.95$) yield almost the same signal-to-noise ratio in the strain data for GW150914.

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

We propose a blind, template-free method for the extraction of a common signal between the Hanford and Livingston detectors and apply it especially to the GW150914 event. We construct a log-likelihood method that maximizes the cross-correlation between each detector and the common signal and minimizes the cross-correlation between the residuals. The reliability of this method is tested using simulations with an injected common signal. Finally, our method is used to assess the quality of theoretical gravitational wave templates for GW150914.

We explore the possibility that bound states involving dark matter particles could be detected by resonance searches at the LHC, and the generic implications of such scenarios for indirect and direct detection. We demonstrate that resonance searches are complementary to mono-jet searches and can probe dark matter masses above 1 TeV with current LHC data. We argue that this parameter regime, where the bound-state resonance channel is the most sensitive probe of the dark sector, arises most naturally in the context of non-trivial dark sectors with large couplings, nearly-degenerate dark-matter-like states, and multiple force carriers. The presence of bound states detectable by the LHC implies a minimal Sommerfeld enhancement that is appreciable, and potentially also radiative bound state formation in the Galactic halo, leading to large signals in indirect searches. We calculate these complementary constraints, which favor either models where the bound-state-forming dark matter constitutes a small fraction of the total density, or models where the late-time annihilation is suppressed at low velocities or late times. We present concrete examples of models that satisfy all these constraints and where the LHC resonance search is the most sensitive probe of the dark sector.

The long-term evolution of a circumstellar disk starting from its formation and ending in the T Tauri phase was simulated numerically with the purpose of studying the evolution of dust in the disk with distinct values of viscous \alpha-parameter and dust fragmentation velocity v_frag. We solved numerical hydrodynamics equations in the thin-disk limit, which are modified to include a dust component consisting of two parts: sub-micron-sized dust and grown dust with a maximum radius a_r. The former is strictly coupled to the gas, while the latter interacts with the gas via friction. The conversion of small to grown dust, dust growth, and dust self-gravity are also considered. We found that the process of dust growth known for the older protoplanetary phase also holds for the embedded phase of disk evolution. The dust growth efficiency depends on the radial distance from the star - a_r is largest in the inner disk and gradually declines with radial distance. In the inner disk, a_r is limited by the dust fragmentation barrier. The process of small-to-grown dust conversion is very fast once the disk is formed. The total mass of grown dust in the disk (beyond 1 AU) reaches tens or even hundreds of Earth masses already in the embedded phase of star formation and even a greater amount of grown dust drifts in the inner, unresolved 1 AU of the disk. Dust does not usually grow to radii greater than a few cm. A notable exception are models with \alpha <= 10^{-3}, in which case a zone with reduced mass transport develops in the inner disk and dust can grow to meter-sized boulders in the inner 10 AU. Grown dust drifts inward and accumulates in the inner disk regions. This effect is most pronounced in the \alpha <= 10^{-3} models where several hundreds of Earth masses can be accumulated in a narrow region of several AU from the star by the end of embedded phase. (abridged).

Filamentary structures are ubiquitous in high-mass star-forming molecular clouds. Their relation with high-mass star formation is still to be understood. Here we report interferometric observations toward 8 filamentary high-mass star-forming clouds. A total of 50 dense cores are identified in these clouds, most of which present signatures of high-mass star formation. Five of them are not associated with any star formation indicators, hence are prestellar core candidates. Evolutionary phases of these cores and their linewidths, temperatures, NH$_3$ abundances, and virial parameters are found to be correlated. In a sub-sample of 4 morphologically well-defined filaments, we find that their fragmentation can not be solely explained by thermal or turbulence pressure support. We also investigate distributions of gas temperatures and non-thermal motions along the filaments, and find a spatial correlation between non-thermal linewidths and star formation activities. We find evidence of gas flows along these filaments, and derive an accretion rate along filaments of $\sim$10$^{-4}$ M$_\odot$ yr$^{-1}$. These results suggest a strong relationship between massive filaments and high-mass star formation, through i) filamentary fragmentation in very early evolutionary phases to form dense cores, ii) accretion flows along filaments that are important for the growth of dense cores and protostars, and iii) enhancement of non-thermal motion in the filaments by the feedback or accretion during star formation.

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

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

We simulated the $R$-band contribution of the host galaxy of TeV $\gamma$-ray BL Lac object Mrk 501 in different aperture sizes and seeing conditions . The intensive observations were run with the 1.02 m optical telescope at Yunnan Observatories from 2010 May 15 to 18. Based on the host subtraction data presented in Nilsson et al. (2007), the subtraction of host galaxy contamination results in significant seeing-brightness correlations. These correlations would lead to illusive large amplitude variations at short timescales, which will mask the intrinsic micro variability, thus gives rise to difficulty in detecting the intrinsic micro variability. Both aperture size and seeing condition influence the flux measurements, but aperture size impact the result more significantly. Based on the parameters of elliptical galaxy provided in Nilsson et al. (1999), we simulated the host contributions of Mrk 501 in the different aperture sizes and seeing conditions. Our simulation data of the host galaxy obviously weaken these significant seeing-brightness correlations for the host-subtracted brightness of Mrk 501, and can help us discover the intrinsic short timescale micro variability. The pure nuclear flux is $\sim$ 8.0 $\rm{mJy}$ in \emph{R} band, i.e., AGN has a magnitude of $R\sim13^{m}_{\cdot}96$.