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

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

Two Mpc-size contact discontinuities have recently been identified in the XMM-Newton and Suzaku X-ray observations of the outskirts of the Perseus cluster (Walker et al. 2020). These structures have been tentatively interpreted as "sloshing cold fronts", which are customarily associated with differential motions of the cluster gas perturbed by a merger. In this study we consider an alternative scenario, namely, that the most prominent discontinuity near the cluster virial radius is the result of the collision between the accretion shock and the "runaway" merger shock. We also discuss the possible origin of the second discontinuity at ~1.2 Mpc.

We investigated the pulsar radio luminosity ($L$), emission efficiency (ratio of radio luminosity to its spin-down power $\dot{E}$), and death line in the diagram of magnetic field (B) versus spin period (P), and found that the dependence of pulsar radio luminosity on its spin-down power ($L-\dot{E}$) is very weak, shown as $L\sim\dot{E}^{0.06}$, which deduces an equivalent inverse correlation between emission efficiency and spin-down power as $\xi\sim \dot{E}^{-0.94}$. Furthermore, we examined the distributions of radio luminosity of millisecond and normal pulsars, and found that, for the similar spin-down powers, the radio luminosity of millisecond pulsars is about one order of magnitude lower than that of the normal pulsars. The analysis of pulsar radio flux suggests that this correlations are not due to a selective effect, but are intrinsic to the pulsar radio emission physics. Their radio radiations may be dominated by the different radiation mechanisms. The cut-off phenomenon of currently observed radio pulsars in B-P diagram is usually referred as the "pulsar death line", which corresponds to $\dot{E}\approx 10^{30}$ erg/s and is obtained by the cut-off voltage of electron acceleration gap in the polar cap model of pulsar proposed by Ruderman and Sutherland. Observationally, this death line can be inferred by the actual observed pulsar flux $S\approx $1mJy and 1kpc distance, together with the maximum radio emission efficiency of 1\%. At present, the actual observed pulsar flux can reach 0.01mJy by FAST telescope, which will arise the observational limit of spin-down power of pulsar as low as $\dot{E}\approx 10^28$ erg/s. This means that the new death line is downward shifted two orders of magnitude, which might be favorably referred as the "observational limit-line", and accordingly the pulsar theoretical model for the cut-off voltage of gap should be heavily modified.

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

In this brief report, we investigate the existence of 4-dimensional static spherically symmetric black holes (BHs) in the Einstein-complex-scalar-Gauss-Bonnet (EcsGB) gravity with an arbitrary potential $V(\phi)$ and a coupling $f(\phi)$ between the scalar field $\phi$ and the Gauss-Bonnet (GB) term. We find that static regular BH solutions with complex scalar hairs do not exist. This conclusion does not depend on the coupling between the GB term and the scalar field, nor on the scalar potential $V(\phi)$ and the presence of a cosmological constant $\Lambda$ (which can be either positive or negative), as longer as the scalar field remains complex and is regular across the horizon.

We studied the filament structures and dense cores in OMC-2,3 region in Orion A North molecular cloud using the high-resolution N2H+ (1-0) spectral cube observed with the Atacama Large Millimeter/Submillimeter Array (ALMA). The filament network over a total length of 2 pc is found to contain 170 intersections and 128 candidate dense cores. The dense cores are all displaced from the infrared point sources (possible young stars), and the major fraction of cores (103) are located around the intersections. Towards the intersections, there is also an increasing trend for the total column density Ntot as well as the the power-law index of the column-density Probability Distribution Function (N-PDF), suggesting that the intersections would in general have more significant gas assembly than the other part of the filament paths. The virial analysis shows that the dense cores mostly have virial mass ratio of alpha_vir=M_vir/M_gas<1.0, suggesting them to be bounded by the self gravity. In the mean time, only about 23 percent of the cores have critical mass ratio of alpha_crit=M_crit/M_gas<1.0, suggesting them to be unstable against core collapse. Combining these results, it shows that the major fraction of the cold starless and possible prestellar cores in OMC-2,3 are being assembled around the intersections, and currently in a gravitationally bound state. But more extensive core collapse and star formation may still require continuous core-mass growth or other perturbatio

The study of planet occurrence as a function of stellar mass is important for a better understanding of planet formation. Estimating stellar mass, especially in the red giant regime, is difficult. In particular, stellar masses of a sample of evolved planet-hosting stars based on spectroscopy and grid-based modelling have been put to question over the past decade with claims they were overestimated. Although efforts have been made in the past to reconcile this dispute using asteroseismology, results were inconclusive. In an attempt to resolve this controversy, we study four more evolved planet-hosting stars in this paper using asteroseismology, and we revisit previous results to make an informed study of the whole ensemble in a self-consistent way. For the four new stars, we measure their masses by locating their characteristic oscillation frequency, $\mathrm{\nu}_{\mathrm{max}}$, from their radial velocity time series observed by SONG. For two stars, we are also able to measure the large frequency separation, $\mathrm{\Delta\nu}$, helped by extended SONG single-site and dual-site observations and new TESS observations. We establish the robustness of the $\mathrm{\nu}_{\mathrm{max}}$-only-based results by determining the stellar mass from $\mathrm{\Delta\nu}$, and from both $\mathrm{\Delta\nu}$ and $\mathrm{\nu}_{\mathrm{max}}$. We then compare the seismic masses of the full ensemble of 16 stars with the spectroscopic masses from three different literature sources. We find an offset between the seismic and spectroscopic mass scales that is mass-dependent, suggesting that the previously claimed overestimation of spectroscopic masses only affects stars more massive than about 1.6 M$_\mathrm{\odot}$.

Quark matter with only $u$ and $d$ quarks ($ud$QM) might be the ground state of baryonic matter at large baryon number $A>A_{\rm min}$. With $A_{\rm min}\gtrsim 300$, this has no direct conflict with the stability of ordinary nuclei. An intriguing test of this scenario is to look for quantum nucleation of $ud$QM inside neutron stars due to their large baryon densities. In this paper, we study the transition rate of cold neutron stars to $ud$ quark stars ($ud$QSs) and the astrophysical implications, considering the relevant theoretical uncertainties and observational constraints. It turns out that a large portion of parameter space predicts an instantaneous transition, and so the observed neutron stars are mostly $ud$QSs. We find this possibility still viable under the recent gravitational wave and pulsar observations, although there are debates on its compatibility with some observations that involve complicated structure of quark matter. The tension could be partially relieved in the two-families scenario, where the high-mass stars ($M\gtrsim2 M_{\odot}$) are all $ud$QSs and the low-mass ones ($M\sim1.4\, M_{\odot}$) are mostly hadronic stars. In this case, the slow transition of the low-mass hadronic stars points to a very specific class of hadronic models with moderately stiff EOSs, and $ud$QM properties are also strongly constrained.

The discovery of gravitational waves, which are ripples of space-time itself, opened a new window to test general relativity, because it predicts that there are only plus and cross polarizations for gravitational waves. For alternative theories of gravity, there may be up to six polarizations. The measurement of the polarization is one of the major scientific goals for future gravitational wave detectors. To evaluate the capability of the detector, we need to use the frequency dependent response functions averaged over the source direction and polarization angle. We derive the full analytical formulas of the averaged response functions for all six possible polarizations and present their asymptotic behaviors based on these analytical formulas. Compared with the numerical simulation, the full analytical formulas are more efficient and valid for any equal-arm interferometric gravitational wave detector without optical cavities in the arms and for a time-delay-interferometry Michelson combination.

The discovery of gravitational waves, which are ripples of space-time itself, opened a new window to test general relativity, because it predicts that there are only plus and cross polarizations for gravitational waves. For alternative theories of gravity, there may be up to six polarizations. The measurement of the polarization is one of the major scientific goals for future gravitational wave detectors. To evaluate the capability of the detector, we need to use the frequency dependent response functions averaged over the source direction and polarization angle. We derive the full analytical formulas of the averaged response functions for all six possible polarizations and present their asymptotic behaviors based on these analytical formulas. Compared with the numerical simulation, the full analytical formulas are more efficient and valid for any equal-arm interferometric gravitational wave detector without optical cavities in the arms and for a time-delay-interferometry Michelson combination.

We carry out a systematic study of the density structure of gas in the Central Molecular Zone (CMZ) in the Galactic center by extracting clumps from the APEX Telescope Large Area Survey of the Galaxy survey at 870 $\mu$m. We find that the clumps follow a scaling of $m = \rho_0 r^3$ which corresponds to a characteristic density of $n_{\rm H_2} = 1.6 \times 10^3\,\rm cm^{-3}$ ($\rho_0 =112\;M_{\odot}\;\rm pc^{-3}$) with a variation of $\approx 0.5\,\rm dex$, where we assumed a gas-to-dust mass ratio of 100. This characteristic density can be interpreted as the result of thermal pressure equilibrium between the molecular gas and the warm ambient interstellar medium. Such an equilibrium can plausibly be established since shear has approximately the same strength as self-gravity. Our findings may explain the fact that star formation in the CMZ is highly inefficient compared to the rest of the Milky Way disk. We also identify a population of clumps whose densities are two orders of magnitudes higher in the vicinity of the Sgr B2 region, which we propose are produced by collisions between the clumps of lower densities. For these collisions to occur, processes such as compressive tides probably have created the appropriate condition by assembling the clumps together.

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

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.

To investigate the use of saliency-map analysis to aid in searches for transient signals, such as fast radio bursts and individual pulses from radio pulsars. We aim to demonstrate that saliency maps provide the means to understand predictions from machine learning algorithms and can be implemented in piplines used to search for transient events. We have implemented a new deep learning methodology to predict whether or not any segment of the data contains a transient event. The algorithm has been trained using real and simulated data sets. We demonstrate that the algorithm is able to identify such events. The output results are visually analysed via the use of saliency maps. We find that saliency maps can produce an enhanced image of any transient feature without the need for de-dispersion or removal of radio frequency interference. Such maps can be used to understand which features in the image were used in making the machine learning decision and to visualise the transient event. Even though the algorithm reported here was developed to demonstrate saliency-map analysis, we have detected, in archival data, a single burst event with dispersion measure of $41$\,cm$^{-3}$pc that is not associated with any currently known pulsar.

Distance measurements of supernova remnants (SNRs) are essential and important. Accurate estimates of physical size, dust masses, and some other properties of SNRs depend critically on accurate distance measurements. However, the determination of SNR distances is still a tough task. Red clump stars (RCs) have a long history been used as standard candles. In this work, we take RCs as tracers to determine the distances to a large group of SNRs in the inner disk. We first select RC stars based on the near-infrared (IR) color-magnitude diagram (CMD). Then, the distance to and extinction of RC stars are calculated. To extend the measurable range of distance, we combine near-IR photometric data from the 2MASS survey with the deeper UKIDSS and VVV surveys. With the help of the Gaia parallaxes, we also remove contaminants including dwarfs and giants. Because an SN explosion compresses the surrounding interstellar medium, the SNR region would become denser and exhibit higher extinction than the surroundings. The distance of a SNR is then recognized by the position where the extinction and its gradient is higher than that of the ambient medium. A total of 63 SNRs' distances in the Galactic inner disk are determined and divided into three Levels A, B, and C with decreasing reliability. The distances to 43 SNRs are well determined with reliability A or B. The diameters and dust masses of SNRs are estimated with the obtained distance and extinction.

The initial stage of star formation is a complex area study because of its high density and low temperature. Under such conditions, many molecules become depleted from the gas phase by freezing out onto dust grains. However, the deuterated species could remain gaseous and are thus ideal tracers. We investigate the gas dynamics and NH$_2$D chemistry in eight massive pre/protocluster clumps. We present NH$_2$D 1$_{11}$-1$_{01}$ (at 85.926 GHz), NH$_3$ (1, 1) and (2, 2) observations in the eight clumps using the PdBI and the VLA, respectively. We find that the distribution between deuterium fractionation and kinetic temperature shows a number density peak at around $T_{\rm kin}=16.1$ K, and the NH$_2$D cores are mainly located at a temperature range of 13.0 to 22.0 K. We detect seven instances of extremely high deuterium fractionation of $1.0 \leqslant D_{\rm frac} \leqslant 1.41$. We find that the NH$_2$D emission does not appear to coincide exactly with either dust continuum or NH$_3$ peak positions, but often surrounds the star-formation active regions. This suggests that the NH$_{2}$D has been destroyed by the central young stellar object (YSO) due to its heating. The detected NH$_2$D lines are very narrow with a median width of $\rm 0.98\pm0.02 km/s$. The extracted NH$_2$D cores are gravitationally bound ($\alpha_{\rm vir} < 1$), are likely prestellar or starless, and can potentially form intermediate-mass or high-mass stars. Using NH$_3$ (1, 1) as a dynamical tracer, we find very complicated dynamical movement, which can be explained by a combined process with outflow, rotation, convergent flow, collision, large velocity gradient, and rotating toroids. High deuterium fractionation strongly depends on the temperature condition. NH$_2$D is a poor evolutionary indicator of high-mass star formation in evolved stages, but a useful tracer in the starless and prestellar cores.

Several types/classes of shocks naturally arise during formation and evolution of galaxy clusters. One such class is represented by accretion shocks, associated with deceleration of infalling baryons. Such shocks, characterized by a very high Mach number, are present even in 1D models of cluster evolution. Another class is composed of "runaway merger shocks", which appear when a merger shock, driven by a sufficiently massive infalling subcluster, propagates away from the main-cluster center. We argue that, when the merger shock overtakes the accretion shock, a new long-living shock is formed that propagates to large distances from the main cluster (well beyond its virial radius) affecting the cold gas around the cluster. We refer to these structures as Merger-accelerated Accretion shocks (MA-shocks) in this paper. We show examples of such MA-shocks in 1D and 3D simulations and discuss their characteristic properties. In particular, (1) MA-shocks shape the boundary separating the hot intracluster medium (ICM) from the unshocked gas, giving this boundary a "flower-like" morphology. In 3D, MA-shocks occupy space between the dense accreting filaments. (2) Evolution of MA-shocks highly depends on the Mach number of the runaway merger shock and the mass accretion rate parameter of the cluster. (3) MA-shocks may lead to the misalignment of the ICM boundary and the splashback radius.

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

The possibility that in the mass range around $10^{-12}\ M_\odot$ most of dark matter constitutes of primordial black holes (PBHs) is a very interesting topic. To produce PBHs with this mass, the primordial scalar power spectrum needs to be enhanced to the order of 0.01 at the scale $k\sim 10^{12}\ \text{Mpc}^{-1}$. The enhanced power spectrum also produces large secondary gravitational waves at the mHz band. A phenomenological delta function power spectrum is usually used to discuss the production of PBHs and secondary gravitational waves. Based on G and k inflations, we propose a new mechanism to enhance the power spectrum at small scales by introducing a non-canonical kinetic term $[1-2G(\phi)]X$ with the function $G(\phi)$ having a peak. Away from the peak, $G(\phi)$ is negligible and we recover the usual slow-roll inflation which is constrained by the cosmic microwave background anisotrpy observations. Around the peak, the slow-roll inflation transiently turns to ultra slow-roll inflation. The enhancement of the power spectrum can be obtained with generic potentials, and there is no need to fine tune the parameters in $G(\phi)$. The energy spectrum $\Omega_{GW}(f)$ of secondary gravitational waves have the characteristic power law behaviour $\Omega_{GW}(f)\sim f^{n}$ and is testable by pulsar timing array and space based gravitational wave detectors.

In this paper, we systematically study spherically symmetric static spacetimes in the framework of Einstein-aether theory, and pay particular attention to the existence of black holes (BHs). In the theory, one scalar and one vector two gravitational modes additionally appear. To avoid the vacuum gravi-\v{C}erenkov radiation, they must all propagate with speeds greater than or at least equal to the speed of light. In the spherical case, only the scalar mode is relevant, so BH horizons are defined by this mode. In the present studies we first clarify several subtle issues. In particular, we find that, out of the five non-trivial field equations, only three are independent, so the problem is well-posed, as now there are only three unknown functions, two, $F(r)$ and $B(r)$, describe the spacetime and one, $A(r)$, describes the aether field. In addition, the two second-order differential equations for $A$ and $F$ are independent of $B$, and once they are found, $B$ is given simply by an algebraic expression of $F,\; A$ and their derivatives. To simplify the problem further, we explore the symmetry of field redefinitions, and work first with the redefined metric and aether field, and then obtain the physical ones by the inverse transformations. These clarifications significantly simplify the computational labor, which is important, as the problem is highly involved mathematically. In fact, it is exactly because of these, we find various numerical BH solutions with an accuracy that is at least two orders higher than previous ones. More important, these BH solutions are the only ones that satisfy the self-consistent conditions and meantime are consistent with all the observational constraints obtained so far. The locations of universal horizons are also identified, together with several other observationally interesting quantities.

A large number of observations from the Parkes 64\,m-diameter radio telescope, recorded with high time resolution, are publicly available. We have re-processed all of the observations obtained during the first four years (from 1997 to 2001) of the Parkes Multibeam receiver system in order to identify transient events and have built a database that records the 568,736,756 pulse candidates generated during this search. We have discovered a new fast radio burst (FRB), FRB~010305, with a dispersion measure of 350$\pm$5\,\,cm$^{-3}\,$pc and explored why so few FRBs have been discovered in data prior to 2001. After accounting for the dispersion smearing across the channel bandwidth and the sky regions surveyed, the number of FRBs is found to be consistent with model predictions. We also present five single pulse candidates from unknown sources, but with Galactic dispersion measures. We extract a diverse range of sources from the database, which can be used, for example, as a training set of data for new software being developed to search for FRBs in the presence of radio frequency interference.

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.

Considering the shortage of comparative studies on pulse temporal characteristics between single-peaked short gamma-ray bursts (sGRBs) and double-peaked sGRBs, we examine the pulse properties of a sample of 100 BATSE sGRBs using the BATSE Time-Tagged Event (TTE) data with a 5-ms resolution. 243 individual pulses of the single-peaked, double-peaked and triple-peaked sGRBs are fitted to get their statistical properties such as the pulse asymmetry, amplitude, peak time, and pulse width, etc. For the double-peaked sGRBs, according to the overlapping ratio between two adjacent peaks, we first define two kinds of double-peaked sGRBs as M-loose and M-tight types and find that most of the first pulses are similar to the single-peaked ones. We study the dependence of the Full Width at Half Maximum (FWHM) on the photon energy among different energy bands. Interestingly, it is found that a power-law relation with an index of -0.4 does exist between the FWHM and the photon energy for the single- and double-peaked sGRBs. More interestingly, we notice that the power-law relation with a positive index also exist for several special short bursts. In view of the three typical timescales of pulses, namely the angular spreading timescale, the dynamic timescale, and the cooling timescale, we propose that the diverse power-law indexes together with the profile evolution of GRB pulse can be used as an effective probe to diagnose the structure and evolution of the relativistically jetted outflows.

The Large Underground Xenon (LUX) dark matter search was a 250-kg active mass dual-phase time projection chamber that operated by detecting light and ionization signals from particles incident on a xenon target. In December 2015, LUX reported a minimum 90% upper C.L. of 6e-46 cm^2 on the spin-independent WIMP-nucleon elastic scattering cross section based on a 1.4e4 kg*day exposure in its first science run. Tension between experiments and the absence of a definitive positive detection suggest it would be prudent to search for WIMPs outside the standard spin-independent/spin-dependent paradigm. Recent theoretical work has identified a complete basis of 14 independent effective field theory (EFT) operators to describe WIMP-nucleon interactions. In addition to spin-independent and spin-dependent nuclear responses, these operators can produce novel responses such as angular-momentum-dependent and spin-orbit couplings. Here we report on a search for all 14 of these EFT couplings with data from LUX's first science run. Limits are placed on each coupling as a function of WIMP mass.

Testing gravity by binary pulsars nowadays becomes a key issue. Screened modified gravity is a kind of scalar-tensor theory with screening mechanism in order to satisfy the tight Solar System tests. In this paper, we investigate how the screening mechanism affects the orbital dynamics of binary pulsars, and calculate in detail the five post-Keplerian (PK) parameters in this theory. These parameters differ from those of general relativity (GR), and the differences are quantified by the scalar charges, which lead to the dipole radiation in this theory. We combine the observables of PK parameters for the ten binary pulsars, respectively, to place the constraints on the scalar charges and possible deviations from GR. The dipole radiation in the neutron star (NS) - white dwarf (WD) binaries leads to more stringent constraints on deviations from GR. The most constraining systems for the scalar charges of NS and WD are PSR~B1913$+$16 and PSR~J1738$+$0333, respectively. The results of all tests exclude significant strong-field deviations and show good agreement with GR.

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.

The BICEP/Keck (BK) experiment aims to detect the imprint of primordial gravitational waves in the Cosmic Microwave Background polarization, which would be direct evidence of the inflation theory. While the tensor-to-scalar ratio has been constrained to be r_0.05 < 0.06 at 95% c.l., further improvements on this upper limit are hindered by polarized Galactic foreground emissions and removal of gravitational lensing polarization. The 30/40 GHz receiver of the BICEP Array (BA) will deploy at the end of 2019 and will constrain the synchrotron foreground with unprecedented accuracy within the BK sky patch. We will show the design of the 30/40 GHz detectors and test results summarizing its performance. The low optical and atmospheric loading at these frequencies requires our TES detectors to have low saturation power in order to be photon-noise dominated. To realize the low thermal conductivity required from a 250 mK base temperature, we developed new bolometer leg designs. We will present the relevant measured detector parameters: G, Tc, Rn, Psat , and spectral bands, and noise spectra. We achieved a per bolometer NEP including all noise components of 2.07E-17 W/sqrt(Hz), including an anticipated photon noise level 1.54E-17 W/sqrt(Hz).

Branches of cosmic inflationary models, such as slow-roll inflation, predict a background of primordial gravitational waves that imprints a unique odd-parity B-mode pattern in the Cosmic Microwave Background (CMB) at amplitudes that are within experimental reach. The BICEP/Keck (BK) experiment targets this primordial signature, the amplitude of which is parameterized by the tensor-to-scalar ratio r, by observing the polarized microwave sky through the exceptionally clean and stable atmosphere at the South Pole. B-mode measurements require an instrument with exquisite sensitivity, tight control of systematics, and wide frequency coverage to disentangle the primordial signal from the Galactic foregrounds. BICEP Array represents the most recent stage of the BK program, and comprises four BICEP3-class receivers observing at 30/40, 95, 150 and 220/270 GHz. The 30/40 GHz receiver will be deployed at the South Pole during the 2019/2020 austral summer. After 3 full years of observations with 30,000+ detectors, BICEP Array will measure primordial gravitational waves to a precision $\sigma(r)$ between 0.002 and 0.004, depending on foreground complexity and the degree of lensing removal. In this paper we give an overview of the instrument, highlighting the design features in terms of cryogenics, magnetic shielding, detectors and readout architecture as well as reporting on the integration and tests that are ongoing with the first receiver at 30/40 GHz.

Families of cosmic inflation models predict a primordial gravitational-wave background that imprints B-mode polarization pattern in the Cosmic Microwave Background (CMB). High sensitivity instruments with wide frequency coverage and well-controlled systematic errors are needed to constrain the faint B-mode amplitude. We have developed antenna-coupled Transition Edge Sensor (TES) arrays for high-sensitivity polarized CMB observations over a wide range of millimeter-wave bands. BICEP Array, the latest phase of the BICEP/Keck experiment series, is a multi-receiver experiment designed to search for inflationary B-mode polarization to a precision $\sigma$(r) between 0.002 and 0.004 after 3 full years of observations, depending on foreground complexity and the degree of lensing removal. We describe the electromagnetic design and measured performance of BICEP Array low-frequency 40-GHz detector, their packaging in focal plane modules, and optical characterization including efficiency and beam matching between polarization pairs. We summarize the design and simulated optical performance, including an approach to improve the optical efficiency due to mismatch losses. We report the measured beam maps for a new broad-band corrugation design to minimize beam differential ellipticity between polarization pairs caused by interactions with the module housing frame, which helps minimize polarized beam mismatch that converts CMB temperature to polarization ($T \rightarrow P$) anisotropy in CMB maps.

The BICEP/Keck experiment (BK) is a series of small-aperture refracting telescopes observing degree-scale Cosmic Microwave Background (CMB) polarization from the South Pole in search of a primordial $B$-mode signature. This $B$-mode signal arises from primordial gravitational waves interacting with the CMB, and has amplitude parametrized by the tensor-to-scalar ratio $r$. Since 2016, BICEP3 and the Keck Array have been observing with 4800 total antenna-coupled transition-edge sensor detectors, with frequency bands spanning 95, 150, 220, and 270 GHz. Here we present the optical performance of these receivers from 2016 to 2019, including far-field beams measured in situ with an improved chopped thermal source and instrument spectral response measured with a field-deployable Fourier Transform Spectrometer. As a pair differencing experiment, an important systematic that must be controlled is the differential beam response between the co-located, orthogonally polarized detectors. We generate per-detector far-field beam maps and the corresponding differential beam mismatch that is used to estimate the temperature-to-polarization leakage in our CMB maps and to give feedback on detector and optics fabrication. The differential beam parameters presented here were estimated using improved low-level beam map analysis techniques, including efficient removal of non-Gaussian noise as well as improved spatial masking. These techniques help minimize systematic uncertainty in the beam analysis, with the goal of constraining the bias on $r$ induced by temperature-to-polarization leakage to be subdominant to the statistical uncertainty. This is essential as we progress to higher detector counts in the next generation of CMB experiments.

Recently, it was shown that quark matter with only $u$ and $d$ quarks ($ud$QM) can be the ground state of matter for baryon numbers $A>A_\textrm{min}$ with $A_{\rm min}\gtrsim 300$. In this paper, we explore $ud$ quark stars ($ud$QSs) that are composed of $ud$QM, in the context of the two-families scenario in which $ud$QSs and hadronic stars (HSs) can coexist. Distinct signatures are discussed compared to the conventional study regarding strange quark stars (SQSs). We show that the requirements of $A_{\rm min}\gtrsim 300$ and the most massive compact star observed being a $ud$QS together may put stringent constraints on the allowed parameter space of $ud$QSs. Then, we study the related gravitational-wave probe of the tidal deformability in binary star mergers, including the $ud$QS-$ud$QS and $ud$QS-HS cases. The obtained values of the tidal deformability at 1.4 solar masses and the average tidal deformability are all in good compatibility with the experimental constraints of GW170817. This study points to a new possible interpretation of the GW170817 binary merger event, where $ud$QS may be at least one component of the binary system detected.

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 describe an on-sky demonstration of a microwave-multiplexing readout system in one of the receivers of the Keck Array, a polarimetry experiment observing the cosmic microwave background at the South Pole. During the austral summer of 2018-2019, we replaced the time-division multiplexing readout system with microwave-multiplexing components including superconducting microwave resonators coupled to radio-frequency superconducting quantum interference devices at the sub-Kelvin focal plane, coaxial-cable plumbing and amplification between room temperature and the cold stages, and a SLAC Microresonator Radio Frequency system for the warm electronics. In the range 5-6 GHz, a single coaxial cable reads out 528 channels. The readout system is coupled to transition-edge sensors, which are in turn coupled to 150-GHz slot-dipole phased-array antennas. Observations began in April 2019, and we report here on an initial characterization of the system performance.

Asteroseismology offers the possibility of probing stellar interiors and testing evolutionary and seismic models. Precise photometry and spectroscopy obtained during multi-site campaigns on young open clusters allows discovering rich samples of pulsating stars and using them in a simultaneous seismic modelling called ensemble asteroseismology. The aim of this study is to obtain the age of the open cluster NGC 6910 by means of ensemble asteroseismology of the early-type pulsating members, to derive their stellar parameters, and to classify the excited modes. We used time-series analysis, performed photometric and spectroscopic mode identification, and calculated grids of evolutionary and seismic models to apply the procedure of ensemble asteroseismology for nine pulsating members of NGC 6910. With two iterations of the procedure of ensemble asteroseismology applied to nine pulsating stars we derived an age of 10.6$^{+0.9}_{-0.8}$ Myr for NGC 6910. Of the nine pulsating stars examined in the paper, eight are $\beta$ Cep stars, including three that are hybrid $\beta$ Cep and slowly pulsating B-type (SPB) pulsators, and one is an SPB star. Interestingly, the least massive $\beta$ Cep star, NGC 6910-38, has a mass of about 5.6 M$_\odot$. The present theory does not predict unstable $p$ modes in B-type stars with such a low mass. The $g$ modes with relatively high frequencies ($>3.5$ d$^{-1}$), observed in three members of the cluster, are also stable according to seismic modelling. Both findings pose a challenge for theoretical calculations and prompt a revision of the opacities. The procedure of ensemble asteroseismology was found to be successful for NGC 6910 and $\chi$ Per on the basis of pulsating B-type stars and can therefore be applied to other young open clusters that are rich in such stars.

We present a novel analysis technique for liquid xenon time projection chambers that allows for a lower threshold by relying on events with a prompt scintillation signal consisting of single detected photons. The energy threshold of the LUX dark matter experiment is primarily determined by the smallest scintillation response detectable, which previously required a 2-fold coincidence signal in its photomultiplier arrays, enforced in data analysis. The technique presented here exploits the double photoelectron emission effect observed in some photomultiplier models at vacuum ultraviolet wavelengths. We demonstrate this analysis using an electron recoil calibration dataset and place new constraints on the spin-independent scattering cross section of weakly interacting massive particles (WIMPs) down to 2.5 GeV/c$^2$ WIMP mass using the 2013 LUX dataset. This new technique is promising to enhance light WIMP and astrophysical neutrino searches in next-generation liquid xenon experiments.

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.

To deepen our understanding of the chemical properties of the Planck Galactic Cold Clump (PGCC) G168.72-15.48, we performed observations of nine molecular species, namely, \ce{c-C3H}, \ce{H2CO}, \ce{HC5N}, \ce{HC7N}, \ce{SO}, \ce{CCH}, \ce{N2H+}, \ce{CH3OH}, and \ce{CH3CCH}, toward two dense cores in PGCC G168.72-15.48 using the Tianma Radio Telescope and Purple Mountain Observatory Telescope. We detected \ce{c-C3H}, \ce{H2CO}, \ce{HC5N}, \ce{N2H+}, \ce{CCH}, and \ce{CH3OH} in both G168-H1 and G168-H2 cores, whereas \ce{HC7N} and \ce{CH3CCH} were detected only in G168-H1 and SO was detected only in G168-H2. Mapping observations reveal that the \ce{CCH}, \ce{N2H+}, \ce{CH3OH}, and \ce{CH3CCH} emissions are well coupled with the dust emission in G168-H1. Additionally, \ce{N2H+} exhibits an exceptionally weak emission in the denser and more evolved G168-H2 core, which may be attributed to the \ce{N2H+} depletion. We suggest that the \ce{N2H+} depletion in G168-H2 is dominated by \ce{N2} depletion, rather than the destruction by CO. The local thermodynamic equilibrium calculations indicate that the carbon-chain molecules of \ce{CCH}, \ce{HC5N}, \ce{HC7N}, and \ce{CH3CCH} are more abundant in the younger G168-H1 core. We found that starless core G168-H1 may have the properties of cold dark clouds based on its abundances of carbon-chain molecules. While, the prestellar core G168-H2 exhibits lower carbon-chain molecular abundances than the general cold dark clouds. With our gas-grain astrochemical model calculations, we attribute the observed chemical differences between G168-H1 and G168-H2 to their different gas densities and different evolutionary stages.

The evolution of neutron star (NS) magnetic field (B-field) has long been an important topic, which is still not yet settled down. Here, we analyze the NS B-fields inferred by the cyclotron resonance scattering features (CRSFs) for the high mass X-ray binaries (HMXBs) and by the magnetic dipole model for the spin-down pulsars. We find that the B-fields of both the 32 NSHMXBs and 28 young pulsars with the supernova remnants follow the log-normal distributions, with the average values of 3.4 * 10^12 G and 4.1 * 10^12 G respectively, which are further verified to come from the same continuous distribution by the statistical tests. These results declaim that the two methods of measuring NS B-fields are reliable for the above two groups of samples. In addition, since the NS-HMXBs have experienced the spin-down phase as the normal pulsars without accretion and then the spin-up phase by accretion, their ages should be about million years (Myrs). Our statistical facts imply that the B-fields of NS-HMXBs have little decayed in their non-accretion spin-down phases of ~ Myrs, as well as in their accretion phases of ~ 0.1Myrs.

We describe the design and assembly of the LUX-ZEPLIN experiment, a direct detection search for cosmic WIMP dark matter particles. The centerpiece of the experiment is a large liquid xenon time projection chamber sensitive to low energy nuclear recoils. Rejection of backgrounds is enhanced by a Xe skin veto detector and by a liquid scintillator Outer Detector loaded with gadolinium for efficient neutron capture and tagging. LZ is located in the Davis Cavern at the 4850' level of the Sanford Underground Research Facility in Lead, South Dakota, USA. We describe the major subsystems of the experiment and its key design features and requirements.

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.

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.

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.

To track the formation and evolution links of the millisecond pulsars (MSPs) powered by accretion and rotation in the Galactic field, we investigate the spin period (P) and spin-down power (Edot) distributions of the MSPs observed at the wavebands of radio, X-ray, and {\gamma}-ray. We find that all but one (119/120) of the {\gamma}-ray MSPs have been detected with the radio signals (radio+{\gamma} MSPs); on the contrary, nearly half of the radio MSPs (118/237) have not been detected with {\gamma}-rays (radio-only MSPs). In addition, the radio+{\gamma} MSPs are shown to be the relative faster and more energetic objects (<P> ~ 3.28 ms and <Edot> ~ 4.5 * 10^34 erg s^-1) compared with the radio-only MSPs (<P> ~ 4.70 ms and <Edot> ~ 1.0 * 10^34 erg s^-1), while the spin periods of these two MSP populations are compatible with the log-normal distributions by the statistical tests. Most rotation-powered MSPs (RMSPs) with the radio eclipsing (31/34) exhibit the radio+{\gamma} signals, which share the faster spin (<P> ~ 2.78 ms) and larger spin-down power (<Edot> ~ 4.1 * 10^34 erg s^-1) distributions than the non-eclipsing ones (<P> ~ 4.19 ms, <Edot> ~ 2.4 * 10^34 erg s^-1), implying the radio+{\gamma} MSPs to be younger than the radio-only MSPs. It is noticed that the spin distribution of the accretion-powered X-ray MSPs shows a clustering phenomenon around ~ 1.6-2.0 ms, which is not observed in RMSPs, hinting that the RMSPs may experience the multiple possible origins. Particularly, all the three super-fast spinning RMSPs with P ~ 1.4-1.6 ms exhibit the non-eclipsing, and we argue that they may be the distinctive sources formed by the accretion induced collapse (AIC) of white dwarfs.

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

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

The scattering of dark matter (DM) particles with sub-GeV masses off nuclei is difficult to detect using liquid xenon-based DM search instruments because the energy transfer during nuclear recoils is smaller than the typical detector threshold. However, the tree-level DM-nucleus scattering diagram can be accompanied by simultaneous emission of a Bremsstrahlung photon or a so-called "Migdal" electron. These provide an electron recoil component to the experimental signature at higher energies than the corresponding nuclear recoil. The presence of this signature allows liquid xenon detectors to use both the scintillation and the ionization signals in the analysis where the nuclear recoil signal would not be otherwise visible. We report constraints on spin-independent DM-nucleon scattering for DM particles with masses of 0.4-5 GeV/c$^2$ using 1.4$\times10^4$ kg$\cdot$day of search exposure from the 2013 data from the Large Underground Xenon (LUX) experiment for four different classes of mediators. This analysis extends the reach of liquid xenon-based DM search instruments to lower DM masses than has been achieved previously.

N131 is a typical infrared dust bubble showing an expanding ringlike shell. We study what kinds of CO line ratios can be used to trace the interaction in the expanding bubble. We carry out new $\rm CO\,(3-2)$ observations towards bubble N131 using the JCMT 15-m telescope, and derive line ratios by combining with our previous $\rm CO\,(2-1)$ and $\rm CO\,(1-0)$ data from the IRAM 30-m telescope observations. To trace the interaction between the molecular gas and the ionized gas in the HII region, we use RADEX to model the dependence of CO line ratios on kinetic temperature and H$_2$ volume density, and examine the abnormal line ratios based on other simulations. We present $\rm CO\,(3-2)$, $\rm CO\,(2-1)$, and $\rm CO\,(1-0)$ integrated intensity maps convolved to the same angular resolution (22.5$"$). The three different CO transition maps show apparently similar morphology. The line ratios of $W_{\rm CO\,(3-2)}$/$W_{\rm CO\,(2-1)}$ mostly range from 0.2 to 1.2 with a median of $0.54\pm0.12$, while the line ratios of $W_{\rm CO\,(2-1)}$/$W_{\rm CO\,(1-0)}$ range from 0.5 to 1.6 with a median of $0.84\pm0.15$. The high CO line ratios $W_{\rm CO\,(3-2)}$/$W_{\rm CO\,(2-1)}\gtrsim 0.8 $ and $W_{\rm CO\,(2-1)}$/$W_{\rm CO\,(1-0)}\gtrsim 1.2$ are beyond the threshold predicted by numerical simulations based on the assumed density-temperature structure for the inner rims of ringlike shell, where are the compressed areas in bubble N131. These high CO integrated intensity ratios, such as $W_{\rm CO\,(3-2)}$/$W_{\rm CO\,(2-1)}\gtrsim0.8$ and $W_{\rm CO\,(2-1)}$/$W_{\rm CO\,(1-0)}\gtrsim1.2$, can be used as a tracer of gas compressed regions with a relatively high temperature and density. This further suggests that the non-Gaussian part of the line-ratio distribution can be used to trace the interaction between the molecular gas and the hot gas in the bubble.

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.

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

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

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.

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.

In this paper, we study the gravitational waveforms, polarizations and radiation powers of the first relativistic triple systems PSR J0337 + 1715, observed in 2014, by using the post-Newtonian approximations to their lowest order. Although they cannot be observed either by current or next generation of the detectors, they do provide useful information to test different theories of gravity. In particular, we carry out the studies in three different theories, general relativity (GR), Brans-Dicke (BD) gravity, and Einstein-aether ($\ae$) theory. The tensor modes $h_{+}$ and $h_{\times}$ exist in all three theories and have almost equal amplitudes. Their frequencies are all peaked at two locations, $ f^{(+, \times)}_1 \simeq 0.0686656 \mu$Hz and $ f^{(+, \times)}_2 \simeq 14.2138 \mu$Hz, which are about twice of the outer and inner orbital frequencies of the triple system, as predicted in GR. In $\ae$-theory, all the six polarization modes are different from zero, but the breathing $h_b$ and longitudinal $h_L$ modes are not independent and also peaked at two frequencies, but at the frequencies, $ f^{(b, L, \ae)}_{1} \simeq 0.0457771 \mu$Hz and $f^{(b, L, \ae)}_{ 2} \simeq 7.09545 \mu$Hz. A similar phenomenon is also observed in BD gravity, in which only the three modes $h_{+},\; h_{\times}$ and $h_{b}$ exit. We also study the radiation powers, and find that the quadrupole emission in each of the three theories has almost the same amplitude, but the dipole emission can be as big as the quadrupole emission in $\ae$-theory. This can provide a very promising window to obtain severe constraints on $\ae$-theory by the multi-band gravitational wave astronomy.

The massive galaxy population above the characteristic Schechter mass M*~10^10.6 Msun contributes to about half of the total stellar mass in the local Universe. These massive galaxies usually reside in hot dark matter haloes above the critical shock-heating mass ~10^12 Msun, where the external cold gas supply to these galaxies is expected to be suppressed. When galaxies run out of their cold gas reservoir, they become dead and quiescent. Therefore, massive quiescent galaxies living in hot haloes are commonly believed to be gas-poor. Based on the data from SDSS, ALFALFA, GASS and COLD GASS surveys, here we show that the vast majority of the massive, quiescent, central disk galaxies in the nearby Universe have a remarkably large amount of cold atomic hydrogen gas, surprisingly similar to star-forming galaxies. Both star-forming and quiescent disk galaxies show identical symmetric double-horn HI spectra, indicating similar regularly rotating HI disks. Relative to their star-forming counterparts, massive quiescent central disk galaxies are quenched because of their significantly reduced molecular gas content, lower dust content, and lower star formation efficiency. Our findings reveal a new picture, which clearly demonstrates the detailed star-formation quenching process in massive galaxies and provides a stringent constraint on the physical mechanism of quenching.

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.

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.

It is well known that aligned, aspherical dust grains emit polarized radiation and that the degree of polarization depends on the angle $\psi$ between the interstellar magnetic field and the line of sight. However, anisotropy of the dust absorption cross sections also modulates the $total\ intensity$ of the radiation as the viewing geometry changes. We report a detection of this effect in the high Galactic latitude $Planck$ data, finding that the 353 GHz dust intensity per $N_{\rm HI}$ is smaller when the Galactic magnetic field is mostly in the plane of the sky and larger when the field is mostly along the line of sight. These variations are of opposite sign and roughly equal magnitude as the changes in polarized intensity per $N_{\rm HI}$ with $\psi$, as predicted. In principle, the variation in intensity can be used in conjunction with the dust polarization angle to constrain the full 3D orientation of the Galactic magnetic field.

Space-based gravitational wave detectors cannot keep rigid structures and precise arm length equality, so the precise equality of detector arms which is required in a ground-based interferometer to cancel the overwhelming laser noise is impossible. The time-delay interferometry method is applied to unequal arm lengths to cancel the laser frequency noise. We give analytical formulas of the averaged response functions for tensor, vector, breathing and longitudinal polarizations in different TDI combinations, and obtain their asymptotic behaviors. At low frequencies, $f\ll f_*$, the averaged response functions of all TDI combinations increase as $f^2$ for all six polarizations. The one exception is the averaged response functions of $\zeta$ for all six polarizations increase as $f^4$ in the equilateral-triangle case. At high frequencies, $f\gg f_*$, the averaged response functions of all TDI combinations for the tensor and breathing modes fall off as $1/f^2$, the averaged response functions of all TDI combinations for the vector mode fall off as $\ln(f)/f^2$ , and the averaged response functions of all TDI combinations for the longitudinal mode fall as $1/f$. We also give LISA and TianQin sensitivity curves in different TDI combinations for tensor, vector, breathing and longitudinal polarizations.

Space-based gravitational wave detectors cannot keep rigid structures and precise arm length equality, so the precise equality of detector arms which is required in a ground-based interferometer to cancel the overwhelming laser noise is impossible. The time-delay interferometry method is applied to unequal arm lengths to cancel the laser frequency noise. We give analytical formulas of the averaged response functions for tensor, vector, breathing and longitudinal polarizations in different TDI combinations, and obtain their asymptotic behaviors. At low frequencies, $f\ll f_*$, the averaged response functions of all TDI combinations increase as $f^2$ for all six polarizations. The one exception is the averaged response functions of $\zeta$ for all six polarizations increase as $f^4$ in the equilateral-triangle case. At high frequencies, $f\gg f_*$, the averaged response functions of all TDI combinations for the tensor and breathing modes fall off as $1/f^2$, the averaged response functions of all TDI combinations for the vector mode fall off as $\ln(f)/f^2$ , and the averaged response functions of all TDI combinations for the longitudinal mode fall as $1/f$. We also give LISA and TianQin sensitivity curves in different TDI combinations for tensor, vector, breathing and longitudinal polarizations.

We present a numerical framework for the variability of active galactic nuclei (AGN), which links the variability of AGN over a broad range of timescales and luminosities to the observed properties of the AGN population as a whole, and particularly the Eddington ratio distribution function (ERDF). We have implemented our framework on GPU architecture, relying on previously published time series generating algorithms. After extensive tests that characterise several intrinsic and numerical aspects of the simulations, we describe some applications used for current and future time domain surveys and for the study of extremely variable sources (e.g., "changing look" or flaring AGN). Specifically, we define a simulation setup which reproduces the AGN variability observed in the PTF/iPTF survey, and use it to forward model longer light curves of the kind that may be observed within the LSST main survey. Thanks to our effcient implementations, these simulations are able to cover for example over 1 Myr with a roughly weekly cadence. We envision that this framework will become highly valuable to prepare for, and best exploit, data from upcoming time domain surveys, such as for example LSST.

Precision measurements of cosmic microwave background (CMB) polarization require extreme control of instrumental systematics. In a companion paper we have presented cosmological constraints from observations with the BICEP2 and Keck Array experiments up to and including the 2015 observing season (BK15), resulting in the deepest CMB polarization maps to date and a statistical sensitivity to the tensor-to-scalar ratio of $\sigma(r) = 0.020$. In this work we characterize the beams and constrain potential systematic contamination from main beam shape mismatch at the three BK15 frequencies (95, 150, and 220 GHz). Far-field maps of 7,360 distinct beam patterns taken from 2010-2015 are used to measure differential beam parameters and predict the contribution of temperature-to-polarization leakage to the BK15 B-mode maps. In the multifrequency, multicomponent likelihood analysis that uses BK15, Planck, and WMAP maps to separate sky components, we find that adding this predicted leakage to simulations induces a bias of $\Delta r = 0.0027 \pm 0.0019$. Future results using higher-quality beam maps and improved techniques to detect such leakage in CMB data will substantially reduce this uncertainty, enabling the levels of systematics control needed for BICEP Array and other experiments that plan to definitively probe large-field inflation.

Moderately strong shocks arise naturally when two subclusters merge. For instance, when a smaller subcluster falls into the gravitational potential of a more massive cluster, a bow shock is formed and moves together with the subcluster. After pericenter passage, however, the subcluster is decelerated by the gravity of the main cluster, while the shock continues moving away from the cluster center. These shocks are considered as promising candidates for powering radio relics found in many clusters. The aim of this paper is to explore the fate of such shocks when they travel to the cluster outskirts, far from the place where the shocks were initiated. In a uniform medium, such a "runaway" shock should weaken with distance. However, as shocks move to large radii in galaxy clusters, the shock is moving down a steep density gradient that helps the shock to maintain its strength over a large distance. Observations and numerical simulations show that, beyond $R_{500}$, gas density profiles are as steep as, or steeper than, $\sim r^{-3}$, suggesting that there exists a "Habitable zone" for moderately strong shocks in cluster outskirts where the shock strength can be maintained or even amplified. A characteristic feature of runaway shocks is that the strong compression, relative to the initial state, is confined to a narrow region just behind the shock. Therefore, if such a shock runs over a region with a pre-existing population of relativistic particles, then the boost in radio emissivity, due to pure adiabatic compression, will also be confined to a narrow radial shell.

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.

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.

We explore the interrelationships between the galaxy group halo mass and various observable group properties. We propose a simple scenario that describes the evolution of the central galaxies and their host dark matter halos. Star formation quenching is one key process in this scenario, which leads to the different assembly histories of blue groups (group with a blue central) and red groups (group with a red central). For blue groups, both the central galaxy and the halo continue to grow their mass. For red groups, the central galaxy has been quenched and its stellar mass remains about constant, while its halo continues to grow by merging smaller halos. From this simple scenario, we speculate about the driving properties that should strongly correlate with the group halo mass. We then apply the machine learning algorithm the Random Forest (RF) regressor to blue groups and red groups separately in the semianalytical model L-GALAXIES to explore these nonlinear multicorrelations and to verify the scenario as proposed above. Remarkably, the results given by the RF regressor are fully consistent with the prediction from our simple scenario and hence provide strong support for it. As a consequence, the group halo mass can be more accurately determined from observable galaxy properties by the RF regressor with a 50% reduction in error. A halo mass more accurately determined in this way also enables more accurate investigations on the galaxy-halo connection and other important related issues, including galactic conformity and the effect of halo assembly bias on galaxy assembly.

We present our piggyback search for fast radio bursts using the Nanshan 26m Radio Telescope and the Kunming 40m Radio Telescope. The observations are performed in the L-band from 1380 MHz to 1700 MHz at Nanshan and S-band from 2170 MHz to 2310 MHz at Kunming. We built the \textsc{Roach2}-based FFT spectrometer and developed the real-time transient search software. We introduce a new radio interference mitigation technique named \emph{zero-DM matched filter} and give the formula of the signal-to-noise ratio loss in the transient search. Though we have no positive detection of bursts in about 1600 and 2400 hours data at Nanshan and Kunming respectively, an intriguing peryton was detected at Nanshan, from which hundreds of bursts were recorded. Perytons are terrestrial radio signals that mimic celestial fast radio bursts. They were first reported at Parkes and identified as microwave oven interferences later. The bursts detected at Nanshan show similar frequency swept emission and have double-peaked profiles. They appeared in different sky regions in about tens of minutes observations and the dispersion measure index is not exactly 2, which indicates the terrestrial origin. The peryton differs drastically from the known perytons detected at Parkes, because it appeared in a precise period of $p=1.71287\pm 0.00004$ s. Its origin remains unknown.

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

Environment is one of the key external drivers of the galaxies, while active galactic nucleus (AGN) is one of the key internal drivers. Both of them play fundamental roles in regulating the formation and evolution of galaxies. We explore the interrelationship between environment and AGN in SDSS. At a given stellar mass, the specific star formation rate distribution of the AGN host galaxies remains unchanged with over-density, with the peak of the distribution around the Green Valley. We show that, at a given stellar mass, the AGN fraction that has been commonly used in previous studies (defined as the number of AGNs relative to all galaxies including passive and star forming ones) does decrease with increasing over-density for satellites. This is largely due to the fact that the fraction of passive galaxies strongly depends on environment. In order to investigate the intrinsic correlation between AGN and environment, especially under the assumption that AGN feedback is responsible for star formation quenching, the AGN fraction should be defined as the number of AGNs relative to the star-forming galaxies only. With the new definition, we find little dependence of AGN fraction on over-density, central/satellite, and group halo mass. There is only marginal evidence that AGN may prefer denser regions, which is possibly due to more frequent interaction of galaxies or higher merger rate in groups. Our results support the scenario that internal secular evolution is the predominant mechanism of triggering AGN activity, while external environment related processes only play a minor role.

In the early stages of star formation, a protostar is deeply embedded in an optically thick envelope such that it is not directly observable. Variations in the protostellar accretion rate, however, will cause luminosity changes that are reprocessed by the surrounding envelope and are observable at submillimeter wavelengths. We searched for submillimeter flux variability toward 12 Planck Galactic Cold Clumps detected by the James Clerk Maxwell Telescope (JCMT)-SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE) survey. These observations were conducted at 850 um using the JCMT/SCUBA-2. Each field was observed three times over about 14 months between 2016 April and 2017 June. We applied a relative flux calibration and achieved a calibration uncertainty of ~ 3.6% on average. We identified 136 clumps across 12 fields and detected four sources with flux variations of ~ 30%. For three of these sources, the variations appear to be primarily due to large-scale contamination, leaving one plausible candidate. The flux change of the candidate may be associated with low- or intermediate-mass star formation assuming a distance of 1.5 kpc, although we cannot completely rule out the possibility that it is a random deviation. Further studies with dedicated monitoring would provide a better understanding of the detailed relationship between submillimeter flux and accretion rate variabilities while enhancing the search for variability in star-forming clumps farther away than the Gould Belt.

Turbulence evolution in a density-stratified medium differs from that of homogeneous isotropic turbulence described by the Kolmogorov picture. We evaluate the degree of this effect in the intracluster medium (ICM) with hydrodynamical simulations. We find that the buoyancy effect induced by ICM density stratification introduces qualitative changes to the turbulence energy evolution, morphology, and the density fluctuation - turbulence Mach number relation, and likely explains the radial dependence of the ICM turbulence amplitude as found previously in cosmological simulations. A new channel of energy flow between the kinetic and the potential energy is opened up by buoyancy. When the gravitational potential is kept constant with time, this energy flow leaves oscillations to the energy evolution, and leads to a balanced state of the two energies where both asymptote to power-law time evolution with slopes shallower than that for the turbulence kinetic energy of homogeneous isotropic turbulence. We discuss that the energy evolution can differ more significantly from that of homogeneous isotropic turbulence when there is a time variation of the gravitational potential. Morphologically, ICM turbulence can show a layered vertical structure and large horizontal vortical eddies in the central regions with the greatest density stratification. In addition, we find that the coefficient in the linear density fluctuation - turbulence Mach number relation caused by density stratification is in general a variable with position and time.

Radio surveys are widely used to study active galactic nuclei. Radio interferometric observations typically trade-off surface brightness sensitivity for angular resolution. Hence, observations using a wide range of baseline lengths are required to recover both bright small-scale structures and diffuse extended emission. We investigate if generative adversarial networks (GANs) can extract additional information from radio data and might ultimately recover extended flux from a survey with a high angular resolution and vice versa. We use a GAN for the image-to-image translation between two different data sets, namely the Faint Images of the Radio Sky at Twenty-Centimeters (FIRST) and the NRAO VLA Sky Survey (NVSS) radio surveys. The GAN is trained to generate the corresponding image cutout from the other survey for a given input. The results are analyzed with a variety of metrics, including structural similarity as well as flux and size comparison of the extracted sources. RadioGAN is able to recover extended flux density within a $20\%$ margin for almost half of the sources and learns more complex relations between sources in the two surveys than simply convolving them with a different synthesized beam. RadioGAN is also able to achieve subbeam resolution by recognizing complicated underlying structures from unresolved sources. RadioGAN generates over a third of the sources within a $20\%$ deviation from both original size and flux for the FIRST to NVSS translation, while for the NVSS to FIRST mapping it achieves almost $30\%$ within this range.

We report results from an extensive set of measurements of the \b{eta}-decay response in liquid xenon.These measurements are derived from high-statistics calibration data from injected sources of both $^{3}$H and $^{14}$C in the LUX detector. The mean light-to-charge ratio is reported for 13 electric field values ranging from 43 to 491 V/cm, and for energies ranging from 1.5 to 145 keV.

Massive clumps tend to fragment into clusters of cores and condensations, some of which form high-mass stars. In this work, we study the structure of massive clumps at different scales, analyze the fragmentation process, and investigate the possibility that star formation is triggered by nearby HII regions. We present a high angular resolution study of a sample of 8 massive proto-cluster clumps. Combining infrared data, we use few-arcsecond resolution radio- and millimeter interferometric data to study their fragmentation and evolution. Our sample is unique in the sense that all the clumps have neighboring HII regions. Taking advantage of that, we test triggered star formation using a novel method where we study the alignment of the centres of mass traced by dust emission at multiple scales. The eight massive clumps have masses ranging from 228 to 2279 $M_\odot$. The brightest compact structures within infrared bright clumps are typically associated with embedded compact radio continuum sources. The smaller scale structures of $R_{\rm eff}$ $\sim$ 0.02 pc observed within each clump are mostly gravitationally bound and massive enough to form at least a B3-B0 type star. Many condensations have masses larger than 8 $M_\odot$ at small scale of $R_{\rm eff}$ $\sim$ 0.02 pc. Although the clumps are mostly infrared quiet, the dynamical movements are active at clump scale ($\sim$ 1 pc). We studied the spatial distribution of the gas conditions detected at different scales. For some sources we find hints of external triggering, whereas for others we find no significant pattern that indicates triggering is dynamically unimportant. This probably indicates that the different clumps go through different evolutionary paths. In this respect, studies with larger samples are highly desired.

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 aim to investigate the impact of the ionized radiation from the M16 HII region on the surrounding molecular cloud and on its hosted star formation. To present comprehensive multi-wavelength observations towards the M16 HII region, we used new CO data and existing infrared, optical, and submillimeter data. The 12CO J=1-0, 13CO J=1-0, and C18O J=1-0 data were obtained with the Purple Mountain Observatory (PMO) 13.7m radio telescope. To trace massive clumps and extract young stellar objects (YSOs) associated with the M16 HII region, we used the ATLASGAL and GLIMPSE I catalogs, respectively. From CO data, we discern a large-scale filament with three velocity components. Because these three components overlap with each other in both velocity and space, the filament may be made of three layers. The M16 ionized gas interacts with the large-scale filament and has reshaped its structure. In the large-scale filament, we find 51 compact cores from the ATLASGAL catalog, 20 of them being quiescent. The mean excitation temperature of these cores is 22.5 K, while this is 22.2 K for the quiescent cores. This high temperature observed for the quiescent cores suggests that the cores may be heated by M16 and do not experience internal heating from sources in the cores. Through the relationship between the mass and radius of these cores, we obtain that 45% of all the cores are massive enough to potentially form massive stars. Compared with the thermal motion, the turbulence created by the nonthermal motion is responsible for the core formation. For the pillars observed towards M16, the H II region may give rise to the strong turbulence.

Gravitational waves are perturbations of the metric of space-time. Six polarizations are possible, although general relativity predicts that only two such polarizations, tensor plus and tensor cross are present for gravitational waves. We give the analytical formulas for the antenna response functions for the six polarizations which are valid for any equal-arm interferometric gravitational-wave detectors without optical cavities in the arms.The response function averaged over the source direction and polarization angle decreases at high frequencies which deteriorates the signal-to-noise ratio registered in the detector. At high frequencies, the averaged response functions for the tensor and breathing modes fall of as $1/f^2$, the averaged response function for the longitudinal mode falls off as $1/f$ and the averaged response function for the vector mode falls off as $\ln(f)/f^2$.

Gravitational waves are perturbations of the metric of space-time. Six polarizations are possible, although general relativity predicts that only two such polarizations, tensor plus and tensor cross are present for gravitational waves. We give the analytical formulas for the antenna response functions for the six polarizations which are valid for any equal-arm interferometric gravitational-wave detectors without optical cavities in the arms.The response function averaged over the source direction and polarization angle decreases at high frequencies which deteriorates the signal-to-noise ratio registered in the detector. At high frequencies, the averaged response functions for the tensor and breathing modes fall of as $1/f^2$, the averaged response function for the longitudinal mode falls off as $1/f$ and the averaged response function for the vector mode falls off as $\ln(f)/f^2$.

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

We have derived new physical quantities for several High-Mass X-ray Binaries (HMXBs) with supergiant (SG) companions through their cyclotron lines. The parameters are: the terminal velocity of the wind, the mass loss rate of the donor, the effective temperature and the magnetic fields. These parameters influence significantly the improvement of the model of accretion. In spite of the variety of their observational properties, the corresponding magnetic field is around B ~ 10^12 G. This result can be constrained by the effects on stellar evolution. In addition, we have performed a segmentation in the parameter space of donors intended for several SG-HMXB listed in our sample set. The parameter space can be categorized into five regimes depending on the possibility of disk formation associated with accretion from the stellar wind. This can give a quantitative clarification of the observed variability and the properties of these objects. We show that, when these systems come into the direct accretion region, systems with corresponding parameters can emit X-rays.

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.

PSR B0919+06 is known for its abnormal emission phenomenon, where the pulse emission window occasionally shifts progressively in longitude and returns afterwards. The physical mechanism behind this phenomenon is still under investigation. In this paper, we present our ultra-wideband observation of this pulsar using the Five-hundred-meter Aperture Spherical radio Telescope (FAST), with simultaneous measurements in the frequency ranges 280-780 MHz and 1250-1550 MHz. We have identified three abnormal events, each of which becomes less apparent as the frequency decreases. At 1400 MHz, the averaged profile slightly shifted after the first and third abnormal events, implying a relationship between abnormal event and profile variation. We also found a linear trend in the left-edge position of the averaged profiles between the first and third event as well as after the third event, suggesting the existence of a slow-drifting mode between the two major events. The second event has a comparatively small shift in phase and is thus categorized as a `small flare state'. During the third event, a sequence of approximately nine pulses was seen to significantly weaken in all frequency bands, likely associated with the pseudo-nulling observed at 150 MHz. A three-component de-composition analysis of the normal averaged profiles shows that the trailing component is dominant at our observing frequencies, while the centre component has a comparatively steeper spectrum. We found the overall flux density in an abnormal event to slightly differ from that in an ordinary state, and the difference shows a frequency dependence. A comparison of the normal, abnormal and dimmed averaged profile indicates that the leading component is likely to be stable in all states.

We aim to reveal the physical properties and chemical composition of the cores in the California molecular cloud (CMC), so as to better understand the initial conditions of star formation. We made a high-resolution column density map (18.2") with Herschel data, and extracted a complete sample of the cores in the CMC with the \textsl{fellwalker} algorithm. We performed new single-pointing observations of molecular lines near 90 GHz with the IRAM 30m telescope along the main filament of the CMC. In addition, we also performed a numerical modeling of chemical evolution for the cores under the physical conditions. We extracted 300 cores, of which 33 are protostellar and 267 are starless cores. About 51\% (137 of 267) of the starless cores are prestellar cores. Three cores have the potential to evolve into high-mass stars. The prestellar core mass function (CMF) can be well fit by a log-normal form. The high-mass end of the prestellar CMF shows a power-law form with an index $\alpha=-0.9\pm 0.1$ that is shallower than that of the Galactic field stellar mass function. Combining the mass transformation efficiency ($\varepsilon$) from the prestellar core to the star of $15\pm 1\%$ and the core formation efficiency (CFE) of 5.5\%, we suggest an overall star formation efficiency of about 1\% in the CMC. In the single-pointing observations with the IRAM 30m telescope, we find that 6 cores show blue-skewed profile, while 4 cores show red-skewed profile. [$\rm {HCO}^{+}$]/[HNC] and [$\rm {HCO}^{+}$]/$\rm [N_{2}H^{+}]$ in protostellar cores are higher than those in prestellar cores; this can be used as chemical clocks. The best-fit chemical age of the cores with line observations is $\sim 5\times 10^4$~years.

We study a merger of the NGC 4839 group with the Coma cluster using X-ray observations from the XMM-Newton and Chandra telescopes. X-ray data show two prominent features: (i) a long (~600 kpc in projection) and bent tail of cool gas trailing (towards south-west) the optical center of NGC 4839, and ii) a "sheath" region of enhanced X-ray surface brightness enveloping the group, which is due to hotter gas. While at first glance the X-ray images suggest that we are witnessing the first infall of NGC 4839 into the Coma cluster core, we argue that a post-merger scenario provides a better explanation of the observed features and illustrate this with a series of numerical simulations. In this scenario, the tail is formed when the group, initially moving to the south-west, reverses its radial velocity after crossing the apocenter, the ram pressure ceases and the ram-pressure-displaced gas falls back toward the center of the group and overshoots it. Shortly after the apocenter passage, the optical galaxy, dark matter and gaseous core move in a north-east direction, while the displaced gas continues moving to the south-west. The "sheath" is explained as being due to interaction of the re-infalling group with its own tail of stripped gas mixed with the Coma gas. In this scenario, the shock, driven by the group before reaching the apocenter, has already detached from the group and would be located close to the famous relic to the south-west of the Coma cluster.

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

Asteroseismic analysis of solar-like stars allows us to determine physical parameters such as stellar mass, with a higher precision compared to most other methods. Even in a well-studied cluster such as the Hyades, the masses of the red giant stars are not well known, and previous mass estimates are based on model calculations (isochrones). The four known red giants in the Hyades are assumed to be clump (core-helium-burning) stars based on their positions in colour-magnitude diagrams, however asteroseismology offers an opportunity to test this assumption. Using asteroseismic techniques combined with other methods, we aim to derive physical parameters and the evolutionary stage for the planet hosting star epsilon Tau, which is one of the four red giants located in the Hyades. We analysed time-series data from both ground and space to perform the asteroseismic analysis. By combining high signal-to-noise (S/N) radial-velocity data from the ground-based SONG network with continuous space-based data from the revised Kepler mission K2, we derive and characterize 27 individual oscillation modes for epsilon Tau, along with global oscillation parameters such as the large frequency separation and the ratio between the amplitude of the oscillations measured in radial velocity and intensity as a function of frequency. The latter has been measured previously for only two stars, the Sun and Procyon. Combining the seismic analysis with interferometric and spectroscopic measurements, we derive physical parameters for epsilon Tau, and discuss its evolutionary status.

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.

The Zwicky Transient Facility (ZTF), a public-private enterprise, is a new time domain survey employing a dedicated camera on the Palomar 48-inch Schmidt telescope with a 47 deg$^2$ field of view and 8 second readout time. It is well positioned in the development of time domain astronomy, offering operations at 10% of the scale and style of the Large Synoptic Survey Telescope (LSST) with a single 1-m class survey telescope. The public surveys will cover the observable northern sky every three nights in g and r filters and the visible Galactic plane every night in g and r. Alerts generated by these surveys are sent in real time to brokers. A consortium of universities which provided funding ("partnership") are undertaking several boutique surveys. The combination of these surveys producing one million alerts per night allows for exploration of transient and variable astrophysical phenomena brighter than r $\sim$ 20.5 on timescales of minutes to years. We describe the primary science objectives driving ZTF including the physics of supernovae and relativistic explosions, multi-messenger astrophysics, supernova cosmology, active galactic nuclei and tidal disruption events, stellar variability, and Solar System objects.

Gravitationally bound hierarchies containing three or more components are very common in our Universe. In this paper we study {\em periodic} gravitational wave (GW) form, their polarizations, response function, its Fourier transform, and energy loss rate of a triple system through three different channels of radiation, the scalar, vector and tensor modes, in Einstein-aether theory of gravity. In the weak-field approximations and with the recently obtained constraints of the theory, we first analyze the energy loss rate of a binary system, and find that the dipole contributions from the scalar and vector modes could be of the order of ${\cal{O}}\left(c_{14}\right){\cal{O}}\left(G_Nm/d\right)^2$, where $c_{14} \; (\equiv c_{1} + c_{4})$ is constrained to $c_{14} \lesssim {\cal{O}}\left(10^{-5}\right)$ by current observations, where $c_i$'s are the four coupling constants of the theory. On the other hand, the "strong-field" effects for a binary system of neutron stars are about six orders lower than that of GR. So, in this paper we ignore these "strong-field" effects and first develop the general formulas to the lowest post-Newtonian order, by taking the coupling of the aether field with matter into account. Within this approximation, we find that the scalar breather mode and the scalar longitudinal mode are all suppressed by a factor of ${\cal{O}}\left(c_{14}\right)$ with respect to the transverse-traceless modes ($h_{+}$ and $h_{\times}$), while the vectorial modes $(h_{X}$ and $h_{Y}$) are suppressed by a factor of $c_{13} \lesssim {\cal{O}}\left(10^{-15}\right)$. Applying the general formulas to a triple system with periodic orbits, we find that the corresponding GW form, response function, and its Fourier transform depend sensitively on both the configuration of the triple system and their orientations with respect to the detectors.

We present long-term, multi-color photometric monitoring and spectroscopic observations of the short-period M-type eclipsing binary BX Tri. Six flare events were recorded in 4 years from 2014 to 2017. Three of them were detected on one night within an orbital cycle. The strongest one was identified on December 23, 2014. With the amplitudes of $\Delta$B=0.48 mag, $\Delta$V=0.28 mag, $\Delta$R=0.10 mag, $\Delta$I=0.02 mag, the total energy due to the flare event was measured to be 4.08 ($\pm$0.24) $\times$$10^{34}$ erg, exceeding the superflare level ($10^{34}$). Based on the observations, the evolutionary status of the binary system as well as the long-term orbital period changes were analyzed. It reveals that BX Tri is probably a semi-detached system with the primary component filling its Roche lobe. The extremely high occurrence of flare events of the binary could be related to the rapid mass transfer between components.

A survey of C2H N=1-0 and N2H+ J=1-0 toward Planck Galactic cold clumps (PGCCs) was performed using the Purple Mountain Observatory's 13.7 m telescope. C2H and N2H+ were chosen to study the chemical evolutionary states of PGCCs. Among 121 observed molecular cores associated with PGCCs, 71 and 58 are detected with C2H N=1-0 and N2H+ J=1-0, respectively. The detected lines of most sources can be fitted with a single component with compatible Vlsr and line widths, which confirms that these PGCC cores are very cold (with gas temperatures 9-21 K) and quiescent while still dominanted by turbulence. The ratio between the column densities of C2H and N2H+ (N(C2H)/N(N2H+)) is found to be a good tracer for the evolutionary states of PGCC cores. Gas-grain chemical model can reproduce the decreasing trend of N(C2H)/N(N2H+) as a function of time. The cores with the lowest abundances of N2H+ (X[N2H+] < 10^{-10}) are the youngest, and have nearly constant abundances of C2H. In evolved cores with X[N2H+] ~ 1E-9, abundances of C2H drop quickly as the exhaustion of carbon atoms. Although these PGCC cores are in different evolutionary states, they are all quite young (<5E5 yr) with N(C2H) > N(N2H+). Mapping observations are carried out toward 20 PGCC cores. The PGCC cores in Cepheus have lower N(C2H)/N(N2H+) and larger line widths compared with those in Taurus. This implies that PGCC cores in Taurus are less chemically evolved than those in Cepheus.

We study how star formation (SF) is quenched in low-redshift disk galaxies with integral-field spectroscopy. We select 131 face-on spiral galaxies with stellar mass greater than $\rm 3\times10^{10}M_\odot$, and with spatially resolved spectrum from MaNGA DR13. We subdivide the sample into four groups based on the offset of their global specific star formation rate (SFR) from the star-forming main sequence and stack the radial profiles of stellar mass and SFR. By comparing the stacked profiles of quiescent and star-forming disk galaxies, we find that the decrease of the global SFR is caused by the suppression of SF at all radii, but with a more significant drop from the center to the outer regions following an inside-out pattern. As the global specific SFR decreases, the central stellar mass, the fraction of disk galaxies hosting stellar bars, and active galactic nuclei (AGNs; including both LINERs and Seyferts) all increase, indicating dynamical processes and AGN feedback are possible contributors to the inside-out quenching of SF in the local universe. However, if we include only Seyferts, or AGNs with ${\rm EW(H\alpha)>3\AA}$, the increasing trend of AGN fraction with decreasing global sSFR disappears. Therefore, if AGN feedback is contributing to quenching, we suspect that it operates in the low-luminosity AGN mode, as indicated by the increasing large bulge mass of the more passive disk galaxies.

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.

Context. Generative models open up the possibility to interrogate scientific data in a more data-driven way. Aims: We propose a method that uses generative models to explore hypotheses in astrophysics and other areas. We use a neural network to show how we can independently manipulate physical attributes by encoding objects in latent space. Methods: By learning a latent space representation of the data, we can use this network to forward model and explore hypotheses in a data-driven way. We train a neural network to generate artificial data to test hypotheses for the underlying physical processes. Results: We demonstrate this process using a well-studied process in astrophysics, the quenching of star formation in galaxies as they move from low-to high-density environments. This approach can help explore astrophysical and other phenomena in a way that is different from current methods based on simulations and observations.

The BICEP and Keck Array experiments are a suite of small-aperture refracting telescopes observing the microwave sky from the South Pole. They target the degree-scale B-mode polarization signal imprinted in the Cosmic Microwave Background (CMB) by primordial gravitational waves. Such a measurement would shed light on the physics of the very early universe. While BICEP2 observed for the first time a B-mode signal at 150 GHz, higher frequencies from the Planck satellite showed that it could be entirely due to the polarized emission from Galactic dust, though uncertainty remained high. Keck Array has been observing the same region of the sky for several years, with an increased detector count, producing the deepest polarized CMB maps to date. New detectors at 95 GHz were installed in 2014, and at 220 GHz in 2015. These observations enable a better constraint of galactic foreground emissions, as presented here. In 2015, BICEP2 was replaced by BICEP3, a 10 times higher throughput telescope observing at 95 GHz, while Keck Array is now focusing on higher frequencies. In the near future, BICEP Array will replace Keck Array, and will allow unprecedented sensitivity to the gravitational wave signal. High resolution observations from the South Pole Telescope (SPT) will also be used to remove the lensing contribution to B-modes.

X-ray observations of merging clusters provide many examples of bow shocks leading merging subclusters. While the Mach number of a shock can be estimated from the observed density jump using Rankine-Hugoniot condition, it reflects only the velocity of the shock itself and is generally not equal to the velocity of the infalling subcluster dark matter halo or to the velocity of the contact discontinuity separating gaseous atmospheres of the two subclusters. Here we systematically analyze additional information that can be obtained by measuring the standoff distance, i.e. the distance between the leading edge of the shock and the contact discontinuity that drives this shock. The standoff distance is influenced by a number of additional effects, e.g. (1) the gravitational pull of the main cluster (causing acceleration/deceleration of the infalling subcluster), (2) the density and pressure gradients of the atmosphere in the main cluster, (3) the non-spherical shape of the subcluster, and (4) projection effects. The first two effects tend to bias the standoff distance in the same direction, pushing the bow shock closer to (farther away from) the subcluster during the pre- (post-)merger stages. Particularly, in the post-merger stage, the shock could be much farther away from the subcluster than predicted by a model of a body moving at a constant speed in a uniform medium. This implies that a combination of the standoff distance with measurements of the Mach number from density/temperature jumps can provide important information on the merger, e.g. differentiating between the pre- and post-merger stages.

We present results from an analysis of all data taken by the BICEP2/Keck CMB polarization experiments up to and including the 2015 observing season. This includes the first Keck Array observations at 220 GHz and additional observations at 95 & 150 GHz. The $Q/U$ maps reach depths of 5.2, 2.9 and 26 $\mu$K$_{cmb}$ arcmin at 95, 150 and 220 GHz respectively over an effective area of $\approx 400$ square degrees. The 220 GHz maps achieve a signal-to-noise on polarized dust emission approximately equal to that of Planck at 353 GHz. We take auto- and cross-spectra between these maps and publicly available WMAP and Planck maps at frequencies from 23 to 353 GHz. We evaluate the joint likelihood of the spectra versus a multicomponent model of lensed-$\Lambda$CDM+$r$+dust+synchrotron+noise. The foreground model has seven parameters, and we impose priors on some of these using external information from Planck and WMAP derived from larger regions of sky. The model is shown to be an adequate description of the data at the current noise levels. The likelihood analysis yields the constraint $r_{0.05}<0.07$ at 95% confidence, which tightens to $r_{0.05}<0.06$ in conjunction with Planck temperature measurements and other data. The lensing signal is detected at $8.8 \sigma$ significance. Running maximum likelihood search on simulations we obtain unbiased results and find that $\sigma(r)=0.020$. These are the strongest constraints to date on primordial gravitational waves.

Various dark matter models predict annual and diurnal modulations of dark matter interaction rates in Earth-based experiments as a result of the Earth's motion in the halo. Observation of such features can provide generic evidence for detection of dark matter interactions. This paper reports a search for both annual and diurnal rate modulations in the LUX dark matter experiment using over 20 calendar months of data acquired between 2013 and 2016. This search focuses on electron recoil events at low energies, where leptophilic dark matter interactions are expected to occur and where the DAMA experiment has observed a strong rate modulation for over two decades. By using the innermost volume of the LUX detector and developing robust cuts and corrections, we obtained a stable event rate of 2.3$\pm$0.2~cpd/keV$_{\text{ee}}$/tonne, which is among the lowest in all dark matter experiments. No statistically significant annual modulation was observed in energy windows up to 26~keV$_{\text{ee}}$. Between 2 and 6~keV$_{\text{ee}}$, this analysis demonstrates the most sensitive annual modulation search up to date, with 9.2$\sigma$ tension with the DAMA/LIBRA result. We also report no observation of diurnal modulations above 0.2~cpd/keV$_{\text{ee}}$/tonne amplitude between 2 and 6~keV$_{\text{ee}}$.

The dust sub-millimetre polarisation of star-forming clouds carries information on dust and the role of magnetic fields in cloud evolution. With observations of a dense filamentary cloud G035.39-00.33, we aim to characterise the dust emission properties and the variations of the polarisation fraction. JCMT SCUBA-2/POL-2 data at 850um are combined with Planck 850um (353GHz) data to map polarisation fractions. With previous SCUBA-2 observations (450um and 850um) and Herschel data, the column densities are determined via modified blackbody fits and via radiative transfer modelling. Models are constructed to examine how the polarisation angles and fractions depend on potential magnetic field geometries and grain alignment. POL-2 data show clear changes in the magnetic field orientation. The filament has a peak column density of N(H2)~7 10^22 cm-2, a minimum dust temperature of T~12 K, and a mass of some 4300Msun for the area N(H2)> 5 10^21 cm-2. The estimated average value of the dust opacity spectral index is beta ~ 1.9. The ratio of sub-millimetre and J band optical depths is tau(250 um)/tau(J) ~ 2.5 10^-3, more than four times the typical values for diffuse medium. The polarisation fraction decreases as a function of column density to p ~ 1% in the central filament. Because of noise, the observed decrease of p(N) is significant only at N(H2)>2 10^22 cm-2. The observations suggest that the grain alignment is not constant. Although the data can be explained with a complete loss of alignment at densities above ~ 10^4 cm-3 or using the predictions of radiative torques alignment, the uncertainty of the field geometry and the spatial filtering of the SCUBA-2 data prevent strong conclusions. G035.39-00.33 shows strong signs of dust evolution and the low polarisation fraction is suggestive of a loss of polarised emission from its densest parts.

We have taken the kHz QPOs as a tool to probe the correlation between the tracks of X-ray color-color diagram (CCD) and magnetosphere-disk positions for the atoll source 4U 1728-34, based on the assumptions that the upper kHz QPO is ascribed to the Keplerian orbital motion and the neutron star (NS) magnetosphere is defined by the dipole magnetic field. We find that from the island to the banana state, the inner accretion disk gradually approaches the NS surface with the radius decreasing from r ~33.0km to ~15.9 km, corresponding to the magnetic field from B(r) ~4.8*10^6 G to ~4.3*10^7 G. In addition, we note the characteristics of some particular radii of magnetosphere-disk -r are: firstly, the whole atoll shape of the CCD links the disk radius range of ~15.9 - 33.0 km, which is just located inside the corotation radius of 4U 1728-34 -r_co ( ~34.4 km), implying that the CCD shape is involved in the NS spin-up state. Secondly, the island and banana states of CCD correspond to the two particular boundaries: (I)-near the corotation radius at r ~27.2 - 33.0 km, where the source lies in the island state; (II)-near the NS surface at r ~15.9 - 22.3 km, where the source lies in both the island and banana states. Thirdly, the vertex of the atoll shape in CCD, where the radiation transition from the hard to soft photons occurs, is found to be near the NS surface at r ~16.4 km. The above results suggest that both the magnetic field and accretion environment are related to the CCD structure of atoll track, where the corotation radius and NS hard surface play the significant roles in the radiation distribution of atoll source.

By using a non-local and time-dependent convection theory, we have calculated radial and low-degree non-radial oscillations for stellar evolutionary models with $M=1.4$--3.0\,$\mathrm{M}_\odot$. The results of our study predict theoretical instability strips for $\delta$ Scuti and $\gamma$ Doradus stars, which overlap with each other. The strip of $\gamma$ Doradus is slightly redder in colour than that of $\delta$ Scuti. We have paid great attention to the excitation and stabilization mechanisms for these two types of oscillations, and we conclude that radiative $\kappa$ mechanism plays a major role in the excitation of warm $\delta$ Scuti and $\gamma$ Doradus stars, while the coupling between convection and oscillations is responsible for excitation and stabilization in cool stars. Generally speaking, turbulent pressure is an excitation of oscillations, especially in cool $\delta$ Scuti and $\gamma$ Doradus stars and all cool Cepheid- and Mira-like stars. Turbulent thermal convection, on the other hand, is a damping mechanism against oscillations that actually plays the major role in giving rise to the red edge of the instability strip. Our study shows that oscillations of $\delta$ Scuti and $\gamma$ Doradus stars are both due to the combination of $\kappa$ mechanism and the coupling between convection and oscillations, and they belong to the same class of variables at the low-luminosity part of the Cepheid instability strip. Within the $\delta$ Scuti--$\gamma$ Doradus instability strip, most of the pulsating variables are very likely hybrids that are excited in both p and g modes.

We have computed linear non-adiabatic oscillations of luminous red giants using a non-local and anisotropic time-dependent theory of convection. The results show that low-order radial modes can be self-excited. Their excitation is the result of radiation and the coupling between convection and oscillations. Turbulent pressure has important effects on the excitation of oscillations in red variables.

Starting from hydrodynamic equations, we have established a set of hydrodynamic equations for average flow and a set of dynamic equations of auto- and cross-correlations of turbulent velocity and temperature fluctuations, following the classic Reynold's treatment of turbulence. The combination of the two sets of equations leads to a complete and self-consistent mathematical expressions ready for the calculations of stellar structure and oscillations. In this paper, non-locality and anisotropy of turbulent convection are concisely presented, together with defining and calibrating of the three convection parameters ($c_1$, $c_2$ and $c_3$) included in the algorithm. With the non-local theory of convection, the structure of the convective envelope and the major characteristics of non-adiabatic linear oscillations are demonstrated by numerical solutions. Great effort has been exercised to the choice of convection parameters and pulsation instabilities of the models, the results of which show that within large ranges of all three parameters ($c_1$, $c_2$ and $c_3$) the main properties of pulsation stability keep unchanged.

Within the framework of non-local time-dependent stellar convection theory, we study in detail the effect of turbulent anisotropy on stellar pulsation stability. The results show that anisotropy has no substantial influence on pulsation stability of g modes and low-order (radial order $n_\mathrm{r}<5$) p modes. The effect of turbulent anisotropy increases as the radial order increases. When turbulent anisotropy is neglected, most of high-order ($n_\mathrm{r}>5$) p modes of all low-temperature stars become unstable. Fortunately, within a wide range of the anisotropic parameter $c_3$, stellar pulsation stability is not sensitive to the specific value of $c_3$. Therefore it is safe to say that calibration errors of the convective parameter $c_3$ do not cause any uncertainty in the calculation of stellar pulsation stability.