Magnetic fields in the intracluster medium (ICM) affect the structure and the evolution of galaxy clusters. However, their properties are largely unknown, and measuring magnetic fields in galaxy clusters is challenging, especially on large-scales outside of individual radio sources. In this work, we probe the plane-of-the-sky orientation of magnetic fields in clusters using the intensity gradients. The technique is a branch of the Gradient Technique (GT) that employs emission intensity maps from turbulent gas. We utilize the Chandra X-ray images of the Perseus, M 87, Coma, and A2597 galaxy clusters, and the VLA radio observations of the synchrotron emission from Perseus. We find that the fields predominantly follow the spiral arms in Perseus, which is in agreement with numerical simulations. The GT-predicted magnetic field shows signatures of magnetic draping around rising bubbles driven by supermassive black hole (SMBH) feedback in the centers of cool-core clusters, as well as draping around substructures merging with the Coma cluster. We calculate the mean-field orientation with respect to the radial direction in these clusters. The orientation is consistent with the prediction of heat-flux-driven buoyancy instability (HBI) in cluster cool cores, and possibly the magnetothermal instability (MTI) in cluster outskirts. There is a broad agreement between the magnetic field of Perseus predicted using the X-ray and radio data. Further numerical studies and better future observations with higher resolution and the larger effective area will help reduce the uncertainties of this method.

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.

The nature of the time correlated noise component (the 1/f noise) of single dish radio telescopes is critical to the detectability of the HI signal in intensity mapping experiments. In this paper, we present the 1/f noise properties of the MeerKAT receiver system using South Celestial Pole (SCP) tracking data. We estimate both the temporal power spectrum density and the 2D power spectrum density for each of the antennas and polarizations. We apply Singular Value Decomposition (SVD) to the dataset and show that, by removing the strongest components, the 1/f noise can be drastically reduced, indicating that it is highly correlated in frequency. Without SVD mode subtraction, the knee frequency over a $20\,$MHz integration is higher than $0.1\,\rm Hz$; with just $2$~mode subtraction, the knee frequency is reduced to $\sim 3\times 10^{-3}\,{\rm Hz}$, indicating that the system induced 1/f-type variations are well under the thermal noise fluctuations over a few hundred seconds time scales. The 2D power spectrum shows that the 1/f-type variations are restricted to a small region in the time-frequency space, either with long wavelength correlations in frequency or in time. This gives a wide range of cosmological scales where the 21cm signal can be measured without further need to calibrate the gain time fluctuations. Finally, we demonstrate that a simple power spectrum parameterization is sufficient to describe the data and provide fitting parameters for both the 1D and 2D power spectrum.

It arises a puzzle in \NGC\, how to secularly maintain the counter-rotating disc from $0.2$ to $7\,$pc unambiguously detected by recent ALMA observations of molecular gas. Upon further analysis of disc dynamics, we find that the Kelvin-Helmholtz (KH) instability (KHI) results in an unavoidable catastrophe of the disc developed at the interface between the reversely rotating parts, and demonstrate that a close binary of supermassive black holes provides tidal torques as the unique external sources to prevent the disc from the KH catastrophe. We are led to the inescapable conclusion that there must be a binary black hole at the center of NGC 1068, to prevent it from the KH catastrophe. The binary is composed of black holes with a separation of $0.1\,$pc from GRAVITY/VLTI observations, a total mass of $1.3\times 10^{7}\:M_{\odot}$ and a mass ratio of $\sim 0.3$ estimated from the angular momentum budge of the global system. The KHI gives rise to forming a gap without cold gas at the velocity interface which overlaps with the observed gap of hot and cold dust regions. Releases of kinematic energies from the KHI of the disc are in agreement with observed emissions in radio and $\gamma$-rays. Such a binary is shrinking with a timescale much longer than the local Hubble time via gravitational waves, however, the KHI leads to an efficient annihilation of the orbital angular momentum and speed up merge of the binary, providing a new paradigm of solving the long term issue of "final parsec problem". Future observations of GRAVITY+/VLTI are expected to be able to spatially resolve the CB-SMBHs suggested in this paper.

A portion of light scalar dark matter, especially axions, may organize into gravitationally bound clumps (stars) and be present in large number in the galaxy today. It is therefore of utmost interest to determine if there are novel observational signatures of this scenario. Work has shown that for moderately large axion-photon couplings, such clumps can undergo parametric resonance into photons, for clumps above a critical mass $M^{\star}_c$ determined precisely by some of us in Ref. [1]. In order to obtain a clump above the critical mass in the galaxy today would require mergers. In this work we perform full 3-dimensional simulations of pairs of axion clumps and determine the conditions under which mergers take place through the emission of scalar waves, including analyzing head-on and non-head-on collisions, phase dependence, and relative velocities. Consistent with other work in the literature, we find that the final mass from the merger $M^{\star}_{\text{final}}\approx 0.7(M^{\star}_1+M^{\star}_2)$ is larger than each of the original clump masses (for $M^{\star}_1\sim M^{\star}_2$). Hence, it is possible for sub-critical mass clumps to merge and become super-critical and therefore undergo parametric resonance into photons. We find that mergers are expected to be kinematically allowed in the galaxy today for high Peccei-Quinn scales, which is strongly suggested by unification ideas, although the collision rate is small. While mergers can happen for axions with lower Peccei-Quinn scales due to statistical fluctuations in relative velocities, as they have a high collision rate. We estimate the collision and merger rates within the Milky Way galaxy today. We find that a merger leads to a flux of energy on earth that can be appreciable and we mention observational search strategies.

A cool phase of the interstellar medium has been observed in many giant elliptical galaxies, but its origin remains unclear. We propose that uneven heating from Type Ia supernovae (SNe Ia), together with radiative cooling, can lead to the formation of the cool phase. The basic idea is that since SNe Ia explode randomly, gas parcels which are not directly heated by SN shocks will cool, forming multiphase gas. We run a series of idealized high-resolution numerical simulations, and find that cool gas develops even when the overall SNe heating rate $H$ exceeds the cooling rate $C$ by a factor as large as 1.4. We also find that the time for multiphase gas development depends on the gas temperature. When the medium has a temperature $T = 3\times 10^6$ K, the cool phase forms within one cooling time \tc; however, the cool phase formation is delayed to a few times \tc\ for higher temperatures. The main reason for the delay is turbulent mixing. Cool gas formed this way would naturally have a metallicity lower than that of the hot medium. For constant $H/C$, there is more turbulent mixing for higher temperature gas. We note that this mechanism of producing cool gas cannot be captured in cosmological simulations, which usually fail to resolve individual SN remnants.

We examine the properties of the circumgalactic medium (CGM) at low redshift in a range of simulated Milky Way mass halos. The sample is comprised of seven idealized simulations, an adaptive mesh refinement cosmological zoom-in simulation, and two groups of 50 halos with star forming or quiescent galaxies taken from the IllustrisTNG100 simulation. The simulations have very different setups, resolution, and feedback models, but are analyzed in a uniform manner. By comparing median radial profiles and mass distributions of CGM properties, we isolate key similarities and differences. In doing so, we advance the efforts of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project that aims to understand the inherently multiscale galaxy formation process. In the cosmological simulations, the CGM exhibits nearly flat temperature distributions, and broad pressure and radial velocity distributions. In the idealized simulations, similar distributions are found in the inner CGM ($\lesssim 0.5 \, r_{\rm 200c}$) when strong galactic feedback models are employed, but the outer CGM ($\gtrsim 0.5 \, r_{\rm 200c}$) has a much less prominent cold phase, and narrower pressure and velocity distributions even in models with strong feedback. This comparative analysis demonstrates the dominant role feedback plays in shaping the inner CGM and the increased importance of cosmological effects, such as nonspherical accretion and satellite galaxies, in the outer CGM. Furthermore, our findings highlight that while cosmological simulations are required to capture the multiphase structure of the CGM at large radii, idealized simulations provide a robust framework to study how galactic feedback interacts with the inner CGM and thereby provide a reliable avenue to constrain feedback prescriptions.

Observations of the time-dependent cosmic-ray Sun shadow have been proven as a valuable diagnostic for the assessment of solar magnetic field models. In this paper, seven years of IceCube data are compared to solar activity and solar magnetic field models. A quantitative comparison of solar magnetic field models with IceCube data on the event rate level is performed for the first time. Additionally, a first energy-dependent analysis is presented and compared to recent predictions. We use seven years of IceCube data for the Moon and the Sun and compare them to simulations on data rate level. The simulations are performed for the geometrical shadow hypothesis for the Moon and the Sun and for a cosmic-ray propagation model governed by the solar magnetic field for the case of the Sun. We find that a linearly decreasing relationship between Sun shadow strength and solar activity is preferred over a constant relationship at the 6.4sigma level. We test two commonly used models of the coronal magnetic field, both combined with a Parker spiral, by modeling cosmic-ray propagation in the solar magnetic field. Both models predict a weakening of the shadow in times of high solar activity as it is also visible in the data. We find tensions with the data on the order of $3\sigma$ for both models, assuming only statistical uncertainties. The magnetic field model CSSS fits the data slightly better than the PFSS model. This is generally consistent with what is found previously by the Tibet AS-gamma Experiment, a deviation of the data from the two models is, however, not significant at this point. Regarding the energy dependence of the Sun shadow, we find indications that the shadowing effect increases with energy during times of high solar activity, in agreement with theoretical predictions.

Type Ia supernovae (SNe Ia) provide unique and important feedback in quiescent galaxies, but their impact has been underappreciated. In this paper, we analyze a series of high-resolution simulations to examine the energetics and turbulence of the medium under SNe Ia. We find that when SN remnants are resolved, their effects differ distinctly from a volumetric heating term, as is commonly assumed in unresolved simulations. First, the net heating is significantly higher than expected, by 30$\pm$10\% per cooling time. This is because a large fraction of the medium is pushed into lower densities which cool inefficiently. Second, the medium is turbulent; the root-mean-squared (RMS) velocity of the gas to 20-50 km s$^{-1}$ on a driving scale of tens of parsec. The velocity field of the medium is dominated by compressional modes, which are larger than the solenoidal components by a factor of 3-7. Third, the hot gas has a very broad density distribution. The ratio between the density fluctuations and the RMS Mach number, parameterized as $b$, is 2-20. This is in contrast to previous simulations of turbulent media, which have found $b\lesssim$ 1. The reason for the difference is mainly caused by the \textit{localized} heating of SNe Ia, which creates a large density contrast. Last, the typical length scale of a density fluctuation grows with time, forming increasingly larger bubbles and filamentary ridges. These underlying density fluctuations need to be included when X-ray observations are interpreted.

We present an archival analysis of Chandra X-ray observations for twelve nearby early-type galaxies hosting radio sources with radio power $>10^{23} \, \rm{W}~\rm{Hz}^{-1}$ at 1.4 GHz, similar to the radio power of the radio source in NGC 4261. Previously, in a similar analysis of eight nearby X-ray and optically-bright elliptical galaxies, Werner et al. 2012, found that NGC 4261 exhibited unusually low central gas entropy compared to the full sample. In the central 0.3 kpc of NGC 4261, the ratio of cooling time to freefall time ($t_{\rm{cool}}/t_{\rm ff}$) is less than $10$, indicating that cold clouds may be precipitating out of the hot ambient medium and providing fuel for accretion in the central region. NGC 4261 also hosts the most powerful radio source in the original sample. Because NGC 4261 may represent an important phase during which powerful feedback from a central active galactic nucleus (AGN) is fueled by multiphase condensation in the central kpc, we searched the Chandra archive for analogs to NGC 4261. We present entropy profiles of those galaxies as well as profiles of $t_{\rm{cool}}/t_{\rm ff}$. We find that one of them, IC 4296, exhibits properties similar to NGC 4261, including the presence of only single phase gas outside of $r \sim 2$ kpc and a similar central velocity dispersion. We compare the properties of NGC 4261 and IC 4296 to hydrodynamic simulations of AGN feedback fueled by precipitation. Over the course of those simulations, the single phase galaxy has an entropy gradient that remains similar to the entropy profiles inferred from our observations.

The past decade has seen significant progress in understanding galaxy formation and evolution using large-scale cosmological simulations. While these simulations produce galaxies in overall good agreement with observations, they employ different sub-grid models for galaxies and supermassive black holes (BHs). We investigate the impact of the sub-grid models on the BH mass properties of the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA simulations, focusing on the M_BH-M_star relation and the BH mass function. All simulations predict tight M_BH-M_star relations, and struggle to produce the lowest (M_BH< 10^7.5 Msun) in galaxies of M_star~10^10.5-10^11.5 Msun. While the time evolution of the mean M_BH-M_star relation is mild (<1 dex in BH mass for 0<z<5) for all the simulations, its linearity (shape) and normalization varies from simulation to simulation. The strength of SN feedback has a large impact on the linearity and time evolution for M_star<10^10.5 Msun. We find that the low-mass end is a good discriminant of the simulation models, and highlights the need for new observational constraints. At the high-mass end, strong AGN feedback can suppress the time evolution of the relation normalization. Compared with the observations of the local universe, we find an excess of BHs with M_BH>10^9 Msun in most of the simulations. The BH mass function is dominated by efficiently accreting BHs (log10 f_Edd >-2$) at high redshifts, and transitions progressively from the high-mass to the low-mass end to be governed by inactive BHs. The transition time and the contribution of active BHs are different among the simulations, and can be used to evaluate models against observations.

Star formation in the universe's most massive galaxies proceeds furiously early in time but then nearly ceases. Plenty of hot gas remains available but does not cool and condense into star-forming clouds. Active galactic nuclei (AGN) release enough energy to inhibit cooling of the hot gas, but energetic arguments alone do not explain why quenching of star formation is most effective in high-mass galaxies. In fact, optical observations show that quenching is more closely related to a galaxy's central stellar velocity dispersion ($\sigma_v$) than to any other characteristic. Here, we show that high $\sigma_v$ is critical to quenching because a deep central potential well maximizes the efficacy of AGN feedback. In order to remain quenched, a galaxy must continually sweep out the gas ejected from its aging stars. Supernova heating can accomplish this task as long as the AGN sufficiently reduces the gas pressure of the surrounding circumgalactic medium (CGM). We find that CGM pressure acts as the control knob on a valve that regulates AGN feedback and suggest that feedback power self-adjusts so that it suffices to lift the CGM out of the galaxy's potential well. Supernova heating then drives a galactic outflow that remains homogeneous if $\sigma_v \gtrsim 240 \, {\rm km \, s^{-1}}$. AGN feedback can effectively quench galaxies with a comparable velocity dispersion, but feedback in galaxies with a much lower velocity dispersion tends to result in convective circulation and accumulation of multiphase gas within the galaxy.

The $^{12}$C+$^{12}$C fusion reaction plays a crucial role in stellar evolution and explosions. Its open reaction channels mainly include $\alpha$, $p$, $n$, and ${}^{8}$Be. Despite more than a half century of efforts, large discrepancies remain among the experimental data measured using various techniques. In this work, we analyze the existing data using the statistical model. Our calculation shows: 1) the relative systematic uncertainties of the predicted branching ratios get smaller as the predicted ratios increase; 2) the total modified astrophysical S-factors (S$^*$ factors) of the $p$ and $\alpha$ channels can each be obtained by summing the S$^*$ factors of their corresponding ground-state transitions and the characteristic $\gamma$ rays while taking into account the contributions of the missing channels to the latter. After applying corrections based on branching ratios predicted by the statistical model, an agreement is achieved among the different data sets at ${E}_{cm}>$4 MeV, while some discrepancies remain at lower energies suggesting the need for better measurements in the near future. We find that the recent S$^*$ factor obtained from an indirect measurement is inconsistent with the direct measurement at energies below 2.6 MeV. We recommend upper and lower limits for the ${}^{12}$C+${}^{12}$C S$^*$ factor based on the existing models. A new $^{12}$C+$^{12}$C reaction rate is also recommended.

The infrared Mg I lines near 12 microns are a pair of emission lines which are magnetically sensitive and have been used to measure solar magnetic fields. Here we calculate the response of the Mg I 12.32 $\mu$m line during a flare and find that in our modeling this line has a complicated behavior. At the beginning of the flare heating, this line shows an intensity dimming at the line center. The intensity then increases when heating continues, with increasing contributions from the heated layers in the chromosphere. The line formation height and the line width also increase as a result. As for the polarized line profiles, we find that flare heating tends to decrease the Zeeman splitting width and attenuates the Stokes $V$ lobe intensity. The wider features in the Stokes $V$ profiles are more pronounced during flare heating, which should be considered when performing magnetic field inversions.

The Tianlai Cylinder Pathfinder is a radio interferometer array designed to test techniques for 21 cm intensity mapping in the post-reionization Universe, with the ultimate aim of mapping the large scale structure and measuring cosmological parameters such as the dark energy equation of state. Each of its three parallel cylinder reflectors is oriented in the north-south direction, and the array has a large field of view. As the Earth rotates, the northern sky is observed by drift scanning. The array is located in Hongliuxia, a radio-quiet site in Xinjiang, and saw its first light in September 2016. In this first data analysis paper for the Tianlai cylinder array, we discuss the sub-system qualification tests, and present basic system performance obtained from preliminary analysis of the commissioning observations during 2016-2018. We show typical interferometric visibility data, from which we derive the actual beam profile in the east-west direction and the frequency band-pass response. We describe also the calibration process to determine the complex gains for the array elements, either using bright astronomical point sources, or an artificial on site calibrator source, and discuss the instrument response stability, crucial for transit interferometry. Based on this analysis, we find a system temperature of about 90 K, and we also estimate the sensitivity of the array.

Detailed modelling of stellar oscillations is able to give precise estimates for stellar ages, but the inferred results typically depend on the adopted model parameters used for the age inference. High-quality asteroseismic data with precise measurements of mixed modes are available for 36 Kepler subgiants. To obtain a handle on the robustness of the ages for these stars, we first study the dependencies of seismic ages on three model input parameters. We find that inferred ages do not change systematically with the helium fraction (Y) or the mixing-length parameter but depend strongly on the metallicity ([M/H]) of the model. The results indicate that age estimates of subgiants have less model dependence and hence are more reliable than those of main-sequence stars or red giants. We then model individual oscillation frequencies of the same 36 Kepler subgiants, using observed metallicities, and obtain their ages with an average precision of ~15%. The comparison with previous age estimates with different stellar codes or input physics shows good agreement (mostly within 2#). We hence suggest that seismology-determined ages of subgiants are not greatly model-dependent.

Distance measurements for extragalactic objects are a fundamental problem in astronomy and cosmology. In the era of precision cosmology, we urgently need better measurements of cosmological distances to observationally test the increasing $H_{0}$ tension of the Hubble constant measured from different tools. Using spectroastrometry, GRAVITY at The Very Large Telescope Interferometer successfully revealed the structure, kinematics and angular sizes of the broad-line region (BLR) of 3C 273 with an unprecedentedly high spatial resolution. Fortunately, reverberation mapping (RM) of active galactic nuclei (AGNs) reliably provides linear sizes of their BLRs. Here we report a joint analysis of spectroastrometry and RM observations to measure AGN distances. We apply this analysis to 3C 273 observed by both GRAVITY and an RM campaign,and find an angular distance of $551.5_{-78.7}^{+97.3}\, {\rm Mpc}$ and $H_{0}=71.5_{-10.6}^{+11.9}\,{\rm km\,s^{-1}\,Mpc^{-1}}$. Advantages of the analysis are 1) its pure geometrical measurements and 2) it simultaneously yields mass of the central black hole in the BLR. Moreover, we can conveniently repeat measurements of selected AGNs to efficiently reduce the statistical and systematic errors. Future observations of a reasonably sized sample ($\sim 30$ AGNs) will provide distances of the AGNs and hence a new way of measuring $H_{0}$ with a high precision $\left(\lesssim 3\%\right)$ to test the $H_{0}$ tension.

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.

Gamma Ray Burst (GRBs) are empirically classified as long-duration GRBs (LGRBs, $>$ 2s) and short-duration GRBs (SGRBs, $<$ 2s). Physically they can be grouped into two distinct progenitor categories: those originating from collapse of massive stars (also known as Type II) and those related to mergers of compact stars (also known as Type I). Even though most LGRBs are Type II and most SGRBs are Type I, the duration criterion is not always reliable to determine the physical category of a certain GRB. Based on our previous comprehensive study of the multi-wavelength properties of long and short GRBs, here we utilize the Naive Bayes method to physically classify GRBs as Type I and Type II GRBs based on multi-wavelength criteria. It results in 0.5\% training error rate and 1\% test error rate. Moreover, there is a gap [-1.2, -0.16] in the distribution of the posterior Odds, $\log O({\rm II:I})$, the Type II to Type I probability ratio. Therefore, we propose to use ${\cal O} = \log O({\rm II:I})+0.7$ as the parameter to classify GRBs into Type I ($<0$) or Type II ($>0$). The only confirmed Type I GRB, GRB 170817A, has log $O({\rm II:I})=-10$. According to this criterion, the supernova-less long GRBs 060614 and 060505 belong to Type I, and two controversial short GRBs 090426 and 060121 belong to Type II.

With telescope apertures becoming larger and larger, the deployment of large-field telescopes is becoming increasingly popular. However, optical path calibration is necessary to ensure the image quality of large-field and large-diameter telescopes. In particular, focal plane attitude calibration is an essential optical path calibration technique that has a direct impact on image quality. In this paper, a focal plane attitude detection method using eight acquisition cameras is proposed based on the calibration requirements of the wide-field telescope, LAMOST. Comparison of simulation and experimental results shows that the detection accuracy of the proposed method can reach 30 arcsec. With additional testing and verification, this method could be used to facilitate regular focal plane attitude calibration for LAMOST as well as other large-field telescopes.

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.

The omnipresence of super-Earths suggests that they are able to be retained in natal disks around low-mass stars, whereas exoplanets' mass distribution indicates that some cores have transformed into gas giants through runaway gas accretion at 1AU from solar-type stars. In this paper, we show that transition to runaway gas accretion by cores may be self-impeded by an increase of the grain opacity in their envelope after they have acquired sufficient mass (typically 10Mearth) to enter a pebble-isolation phase. The accumulation of mm-m size pebbles in their migration barriers enhances their local fragmentation rates. The freshly produced sub-mm grains pass through the barrier, elevate the effective dust opacity and reduce the radiative flux in the cores envelope. These effects alone are adequate to suppress the transition to runaway accretion and preserve super-Earths in the stellar proximity (0.1 AU), albeit entropy advection between the envelope and the disk can further reduce the accretion rate. At intermediate distance (1AU) from their host stars, the escalation in the dust opacity dominates over entropy advection in stalling the transition to runaway accretion for marginally pebble-isolated cores. Beyond a few AU, the transformation of more massive cores to gas giants is reachable before severe depletion of disk gas. This requirement can be satisfied either in extended disks with large scale height via orderly accretion of migrating pebbles or through the mergers of oligarchic protoplanetary embryos, and can account for the correlated occurrence of long-period gas giants and close-in super-Earths.

3C 273 is an intensively monitored flat-spectrum radio quasar with both a beamed jet and blue bump together with broad emission lines. The coexistence of the comparably prominent jet and accretion disk leads to complicated variability properties. Recent reverberation mapping monitoring for 3C 273 revealed that the optical continuum shows a distinct long-term trend that does not have a corresponding echo in the Hbeta fluxes. We compile multi-wavelength monitoring data from the Swift archive and other ground-based programs and clearly find two components of emissions at optical wavelength. One component stems from the accretion disk itself and the other component can be ascribed to the jet contribution, which also naturally accounts for the non-echoed trend in reverberation mapping data. We develop an approach to decouple the optical emissions from the jet and accretion disk in 3C 273 with the aid of multi-wavelength monitoring data. By assuming the disk emission has a negligible polarization in consideration of the low inclination of the jet, the results show that the jet contributes a fraction of ~10% at the minimum and up to ~40% at the maximum to the total optical emissions. This is the first time to provide a physical interpretation to the "detrending" manipulation conventionally adopted in reverberation mapping analysis. Our work also illustrates the importance of appropriately analyzing variability properties in cases of coexisting jets and accretion disks.

Asteroseismology is a powerful tool for probing the internal structures of stars by using their natural pulsation frequencies. It relies on identifying sequences of pulsation modes that can be compared with theoretical models, which has been done successfully for many classes of pulsators, including low-mass solar-type stars, red giants, high-mass stars and white dwarfs. However, a large group of pulsating stars of intermediate mass--the so-called delta Scuti stars--have rich pulsation spectra for which systematic mode identification has not hitherto been possible. This arises because only a seemingly random subset of possible modes are excited, and because rapid rotation tends to spoil the regular patterns. Here we report the detection of remarkably regular sequences of high-frequency pulsation modes in 60 intermediate-mass main-sequence stars, allowing definitive mode identification. Some of these stars have space motions that indicate they are members of known associations of young stars, and modelling of their pulsation spectra confirms that these stars are indeed young.

The presence of mixed modes makes subgiants excellent targets for asteroseismology, providing a probe for the internal structure of stars. Here we study 36 \emph{Kepler} subgiants with solar-like oscillations and report their oscillation mode parameters. We performed a so-called peakbagging exercise, i.e. estimating oscillation mode frequencies, linewidths, and amplitudes with a power spectrum model, fitted in the Bayesian framework and sampled with a Markov Chain Monte Carlo algorithm. The uncertainties of the mode frequencies have a median value of 0.180 $\mu$Hz. We obtained seismic parameters from the peakbagging, analysed their correlation with stellar parameters, and examined against scaling relations. The behaviour of seismic parameters (e.g. $\Delta\nu$, $\nu_{\rm max}$, $\epsilon_p$) is in general consistent with theoretical predictions. We presented the observational p--g diagrams: $\gamma_1$--$\Delta\nu$ for early subgiants and $\Delta\Pi_1$--$\Delta\nu$ for late subgiants, and demonstrate their capability to estimate stellar mass. We also found a $\log g$ dependence on the linewidths and a mass dependence on the oscillation amplitudes and the widths of oscillation excess. This sample will be valuable constraints for modelling stars and studying mode physics such as excitation and damping.

We performed an intensive accretion disk reverberation mapping campaign on the high accretion rate active galactic nucleus Mrk 142 in early 2019. Mrk 142 was monitored with the Neil Gehrels Swift Observatory for 4 months in X-rays and 6 UV/optical filters. Ground-based photometric monitoring was obtained from the Las Cumbres Observatory, Liverpool Telescope and Dan Zowada Memorial Observatory in ugriz filters and the Yunnan Astronomical Observatory in V. Mrk 142 was highly variable throughout, displaying correlated variability across all wavelengths. We measure significant time lags between the different wavelength light curves, finding that through the UV and optical the wavelength-dependent lags, $\tau(\lambda)$, generally follow the relation $\tau(\lambda) \propto \lambda^{4/3}$, as expected for the $T\propto R^{-3/4}$ profile of a steady-state optically-thick, geometrically-thin accretion disk, though can also be fit by $\tau(\lambda) \propto \lambda^{2}$, as expected for a slim disk. The exceptions are the u and U band, where an excess lag is observed, as has been observed in other AGN and attributed to continuum emission arising in the broad-line region. Furthermore, we perform a flux-flux analysis to separate the constant and variable components of the spectral energy distribution, finding that the flux-dependence of the variable component is consistent with the $f_\nu\propto\nu^{1/3}$ spectrum expected for a geometrically-thin accretion disk. Moreover, the X-ray to UV lag is significantly offset from an extrapolation of the UV/optical trend, with the X-rays showing a poorer correlation with the UV than the UV does with the optical. The magnitude of the UV/optical lags is consistent with a highly super-Eddington accretion rate.

Before massive stars heavier than $(8 \cdots 10)$ solar masses evolve to the phase of a gravitational core collapse, they will emit a huge number of MeV-energy neutrinos that are mainly produced in the thermal processes and nuclear weak interactions. The detection of such pre-supernova (pre-SN) neutrinos could provide an important and independent early warning for the optical observations of core-collapse SNe. In this paper, we investigate the prospects of future large liquid-scintillator detectors for the observation of pre-SN neutrinos in both $\bar{\nu}^{}_e + p \to e^+ + n$ and $\nu (\bar{\nu}) + e^- \to \nu (\bar{\nu}) + e^-$ reaction channels, where $\nu$ ($\bar{\nu}$) denotes neutrinos (antineutrinos) of all three flavors. We propose a quantitative assessment of the capability in terms of three working criteria, namely, how far the SN distance can be covered, how long the early warning before the core collapse can be sent out, and how well the direction pointing to the SN can be determined. The dependence of the final results on the different models of progenitor stars, neutrino flavor conversions and the relevant backgrounds is also discussed.

Obtaining measurements of chromospheric and photometric activity of stars with near-solar fundamental parameters and rotation periods is important for a better understanding of solar-stellar connection. We select a sample of 2603 stars with near-solar fundamental parameters from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST)-Kepler field and use LAMOST spectra to measure their chromospheric activity and Kepler light curves to measure their photospheric activity (i.e., the amplitude of the photometric variability). While the rotation periods of 1556 of these stars could not be measured due to the low amplitude of the photometric variability and highly irregular temporal profile of light curves, 254 stars were further identified as having near-solar rotation periods. We show that stars with near-solar rotation periods have chromospheric activities that are systematically higher than stars with undetected rotation periods. Furthermore, while the solar level of photospheric and chromospheric activity appears to be typical for stars with undetected rotation periods, the Sun appears to be less active than most stars with near-solar rotation periods (both in terms of photospheric and chromospheric activity).

We investigate whether the hard X-ray photon index (${\Gamma}$) versus accretion rate correlation for super-Eddington accreting quasars is different from that for sub-Eddington accreting quasars. We construct a sample of 113 bright quasars from the Sloan Digital Sky Survey Data Release 14 quasar catalog, including 38 quasars as the super-Eddington subsample and 75 quasars as the sub-Eddington subsample. We derive black-hole masses using a simple-epoch virial mass formula based on the ${\rm H\beta}$ lines, and we use the standard thin disk model to derive the dimensionless accretion rates ($\dot{\mathscr{M}}$) for our sample. The X-ray data for these quasars are collected from the Chandra and XMM-Newton archives. We fit the hard X-ray spectra using a single power-law model to obtain ${\Gamma}$ values. We find a statistically significant ($R_{\rm S}=0.43$, $p=7.75\times{10}^{-3}$) correlation between ${\Gamma}$ and $\dot{\mathscr{M}}$ for the super-Eddington subsample. The ${\Gamma}$-$\dot{\mathscr{M}}$ correlation for the sub-Eddington subsample is also significant, but weaker ($R_{\rm S}=0.30$, $p=9.98\times{10}^{-3}$). Linear regression analysis shows that ${\rm \Gamma}=(0.34\pm0.11){\rm log}{\dot{\mathscr{M}}}+(1.71\pm0.17)$ and ${\rm \Gamma}=(0.09\pm0.04){\rm log}{\dot{\mathscr{M}}}+(1.93\pm0.04)$ for the super- and sub-Eddington subsamples, respectively. The ${\Gamma}$-$\dot{\mathscr{M}}$ correlations of the two subsamples are different, suggesting different disk-corona connections in these two types of systems. We propose one qualitative explanation of the steeper ${\Gamma}$-$\dot{\mathscr{M}}$ correlation in the super-Eddington regime that involves larger seed photon fluxes received by the compact coronae from the thick disks in super-Eddington accreting quasars.

Recent localizations of Fast Radio Bursts and identifications of their host galaxies confirmed their extragalactic origin. While FRB 121102 resides in the bright region of a dwarf star forming galaxy, other FRBs reside in more massive galaxies and are related to older stellar populations. We compare the host galaxy properties of the five FRBs with those of several types of stellar transients: from young to old population, long duration gamma ray bursts (LGRBs), superluminous supernovae (SLSNe), Type Ib/Ic supernovae (SN Ibc), Type II supernovae (SN II), type Ia supernovae (SN Ia), and short duration gamma ray bursts (SGRBs). We find that as a whole sample, the stellar mass and star formation rate of the five FRB host galaxies prefer a medium to old population, and are against a young population similar to LGRBs and SLSNe by a null probability 0.05. Individually, the FRB 121102 host is consistent with that of young population objects, especially SLSNe; the FRB 180924 environment is similar to that of SGRBs; and the FRB 190523 environment is similar to those of SN Ia or SGRBs. These results are consistent with the magnetar engine model for FRBs, if both magnetars produced from extreme explosions (GRBs/SLSNe) and from regular channels (e.g. those producing Galactic magnetars) can produce FRBs.

We present temporal variations of the Si IV line profiles at the flare ribbons in three solar flares observed by the Interface Region Imaging Spectrograph (IRIS). In the M1.1 flare on 2014 September 6 and the X1.6 flare on 2014 September 10, the Si IV line profiles evolve from wholly redshifted to red-wing enhanced with the flare development. However, in the B1.8 flare on 2016 December 2, the Si IV line profiles are wholly redshifted throughout the flare evolution. We fit the wholly redshifted line profiles with a single Gaussian function but the red-asymmetric ones with a double Gaussian function to deduce the corresponding Doppler velocities. In addition, we find that hard X-ray emission above 25 keV shows up in the two large flares, implying a nonthermal electron beam heating. In the microflare, there only appears weak hard X-ray emission up to 12 keV, indicative of a thermal heating mostly. We interpret the redshifts or red asymmetries of the Si IV line at the ribbons in the three flares as spectral manifestations of chromospheric condensation. We propose that whether the line appears to be wholly redshifted or red-asymmetric depends on the heating mechanisms and also on the propagation of the condensation.

Galaxy clusters are the most massive collapsed structures in the universe whose potential wells are filled with hot, X-ray emitting intracluster medium. Observations however show that a significant number of clusters (the so-called cool-core clusters) also contain large amounts of cold gas in their centres, some of which is in the form of spatially extended filaments spanning scales of tens of kiloparsecs. These findings have raised questions about the origin of the cold gas, as well as its relationship with the central active galactic nucleus (AGN), whose feedback has been established as a ubiquitous feature in such galaxy clusters. Here we report a radiation hydrodynamic simulation of AGN feedback in a galaxy cluster, in which cold filaments form from the warm, AGN-driven outflows with temperatures between $10^4$ and $10^7$ K as they rise in the cluster core. Our analysis reveals a new mechanism, which, through the combination of radiative cooling and ram pressure, naturally promotes outflows whose cooling time is shorter than their rising time, giving birth to spatially extended cold gas filaments. Our results strongly suggest that the formation of cold gas and AGN feedback in galaxy clusters are inextricably linked and shed light on how AGN feedback couples to the intracluster medium.

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.

Recent observations provide evidence that some cool-core clusters (CCCs) host quasars in their brightest cluster galaxies (BCGs). Motivated by these findings we use 3D radiation-hydrodynamic simulations with the code Enzo to explore the joint role of the kinetic and radiative feedback from supermassive black holes (SMBHs) in BCGs. We implement kinetic feedback as sub-relativistic plasma outflows and model radiative feedback using the ray-tracing radiative transfer or thermal energy injection. In our simulations the central SMBH transitions between the radiatively efficient and radiatively inefficient states on timescales of a few Gyr, as a function of its accretion rate. The timescale for this transition depends primarily on the fraction of power allocated to each feedback mode, and to a lesser degree on the overall feedback luminosity of the active galactic nucleus (AGN). Specifically, we find that (a) kinetic feedback must be present at both low and high accretion rates in order to prevent the cooling catastrophe, and (b) its contribution likely accounts for > 10% of the total AGN feedback power, since below this threshold simulated BCGs tend to host radio-loud quasars most of the time, in apparent contrast with observations. We also find a positive correlation between the AGN feedback power and the mass of the cold gas filaments in the cluster core, indicating that observations of H$\alpha$ filaments can be used as a measure of AGN feedback.

The IRIS Mg II k line serves as a very good tool to diagnose the heating processes in solar flares. Recent studies have shown that apart from the usual red asymmetries which are interpreted as the result of condensation downflows, this line could also show a blue-wing enhancement. To investigate how such a blue asymmetry is formed, we perform a grid of radiative hydrodynamic simulations and calculate the corresponding line profiles. We find that such a spectral feature is likely to originate from the upward plasma motion in the upper chromosphere. However, the formation region that is responsible for the blue-wing enhancement could be located in an evaporation region, in an upward moving blob, and even an upward moving condensation region. We discuss how the electron beam parameters affect these different dynamics of the atmosphere.

This paper reports the discovery of an Algol system KIC 10736223 that just past the rapid mass transfer stage. From the light curve and radial-velocity modelling we find KIC 10736223 to be a detached Algol system with the less-massive secondary nearly filling its Roche lobe. Based on the short-cadence Kepler data, we analyzed intrinsic oscillations of the pulsator and identified six secured independent $\delta$ Scuti-type pulsation modes ($f_{1}$, $f_3$, $f_{9}$, $f_{19}$, $f_{42}$, and $f_{48}$). We compute two grids of theoretical models to reproduce the $\delta$ Scuti freqiencies, and find fitting results of mass-accreting models meet well with those of single-star evolutionary models. The fundamental parameters of the primary star yielded with asteroseismology are $M$ = $1.57^{+0.05}_{-0.09}$ $M_{\odot}$, $Z$ = 0.009 $\pm$ 0.001, $R$ = $1.484^{+0.016}_{-0.028}$ $R_{\odot}$, $\log g$ = $4.291^{+0.004}_{-0.009}$, $T_{\rm eff}$ = $7748^{+230}_{-378}$ K, $L$ = $7.136^{+1.014}_{-1.519}$ $L_{\odot}$. The asteroseismic parameters match well with the dynamical parameters derived from the binary model. Moreover, our asteroseismic results show that the pulsator is an almost unevolved star with an age between 9.46-11.65 Myr for single-star evolutionary models and 2.67-3.14 Myr for mass-accreting models. Thereofore, KIC 10736223 may be an Algol system that has just undergone the rapid mass-transfer process.

The separation and elongation motions of solar flare ribbons have extensively been investigated. The degradation and disappearance of ribbons have rarely been explored. In this paper, we report our multiwavelength observations of a C5.5 circular-ribbon flare associated with two jets (jet1 and jet2) on 2014 August 24, focusing on the fast degradation of the outer circular ribbon (CR). The flare, consisting of a short inner ribbon (IR) and outer CR, was triggered by the eruption of a minifilament. The brightness of IR and outer CR reached their maxima simultaneously at $\sim$04:58 UT in all AIA wavelengths. Subsequently, the short eastern part of CR faded out quickly in 1600 {\AA} but gradually in EUV wavelengths. The long western part of CR degraded in the counterclockwise direction and experienced a deceleration. The degradation was distinctly divided into two phases: phase I with faster apparent speeds (58$-$69 km s$^{-1}$) and phase II with slower apparent speeds (29$-$35 km s$^{-1}$). The second phase stopped at $\sim$05:10 UT when the western CR totally disappeared. Besides the outward propagation of jet1, the jet spire experienced untwisting motion in the counterclockwise direction during 04:55$-$05:00 UT. We conclude that the event can be explained by the breakout jet model. The coherent brightenings of the IR and CR at $\sim$04:58 UT may result from the impulsive interchange reconnection near the null point, whereas sub-Alfv\'{e}nic slipping motion of the western CR in the counterclockwise direction indicates the occurrence of slipping magnetic reconnection. Another possible explanation of the quick disappearance of the hot loops connecting to the western CR is that they are simply reconnected sequentially without the need for significant slippage after the null point reconnection.

We present a survey of molecular outflows across the dark cloud complex in the Cygnus region, based on 46.75 deg^2 field of CO isotopologues data from Milky Way Imaging Scroll Painting (MWISP) survey. A supervised machine learning algorithm, Support Vector Machine (SVM), is introduced to accelerate our visual assessment of outflow features in the data cube of 12CO and 13CO J = 1-0 emission. A total of 130 outflow candidates are identified, of which 77 show bipolar structures and 118 are new detections. Spatially, these outflows are located inside dense molecular clouds and some of them are found in clusters or in elongated linear structures tracing the underlying gas filament morphology. Along the line of sight, 97, 31, and 2 candidates reside in the Local, Perseus, and Outer arm, respectively. Young stellar objects as outflow drivers are found near most outflows, while 36 candidates show no associated source. The clusters of outflows that we detect are inhomogeneous in their properties; nevertheless, we show that the outflows cannot inject turbulent energy on cloud scales. Instead, at best, they are restricted to affecting the so called "clump" and "core" scales, and this only on short (~0.3 Myr) estimated timescales. Combined with outflow samples in the literature, our work shows a tight outflow mass-size correlation.

Using the IceCube Neutrino Observatory, we search for high-energy neutrino emission coincident with compact binary mergers observed by the LIGO and Virgo gravitational wave (GW) detectors during their first and second observing runs. We present results from two searches targeting emission coincident with the sky localization of each gravitational wave event within a 1000 second time window centered around the reported merger time. One search uses a model-independent unbinned maximum likelihood analysis, which uses neutrino data from IceCube to search for point-like neutrino sources consistent with the sky localization of GW events. The other uses the Low-Latency Algorithm for Multi-messenger Astrophysics, which incorporates astrophysical priors through a Bayesian framework and includes LIGO-Virgo detector characteristics to determine the association between the GW source and the neutrinos. No significant neutrino coincidence is seen by either search during the first two observing runs of the LIGO-Virgo detectors. We set upper limits on the time-integrated neutrino emission within the 1000 second window for each of the 11 GW events. These limits range from 0.02-0.7 $\mathrm{GeV~cm^{-2}}$. We also set limits on the total isotropic equivalent energy, $E_{\mathrm{iso}}$, emitted in high-energy neutrinos by each GW event. These limits range from 1.7 $\times$ 10$^{51}$ - 1.8 $\times$ 10$^{55}$ erg. We conclude with an outlook for LIGO-Virgo observing run O3, during which both analyses are running in real time.

The Cosmology Large Angular Scale Surveyor (CLASS) is a telescope array that observes the cosmic microwave background (CMB) over 75% of the sky from the Atacama Desert, Chile, at frequency bands centered near 40, 90, 150, and 220 GHz. CLASS measures the large angular scale ($1^\circ\lesssim\theta\leqslant 90^\circ$) CMB polarization to constrain the tensor-to-scalar ratio at the $r\sim0.01$ level and the optical depth to last scattering to the sample variance limit. This paper presents the optical characterization of the 40 GHz telescope during its first observation era, from 2016 September to 2018 February. High signal-to-noise observations of the Moon establish the pointing and beam calibration. The telescope boresight pointing variation is $<0.023^\circ$ ($<1.6$% of the beam's full width at half maximum (FWHM)). We estimate beam parameters per detector and in aggregate, as in the CMB survey maps. The aggregate beam has an FWHM of $1.579^\circ\pm.001^\circ$ and a solid angle of $838 \pm 6\ \mu{\rm sr}$, consistent with physical optics simulations. The corresponding beam window function has a sub-percent error per multipole at $\ell < 200$. An extended $90^\circ$ beam map reveals no significant far sidelobes. The observed Moon polarization shows that the instrument polarization angles are consistent with the optical model and that the temperature-to-polarization leakage fraction is $<10^{-4}$ (95% C.L.). We find that the Moon-based results are consistent with measurements of M42, RCW 38, and Tau A from CLASS's CMB survey data. In particular, Tau A measurements establish degree-level precision for instrument polarization angles.

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

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.

Pulsar wind nebulae (PWNe) are main gamma-ray emitters in the Galactic plane. They are diffuse nebulae that emit non-thermal radiation. Pulsar winds, relativistic magnetized outflows from the central star, shocked in the ambient medium produce a multiwavelength emission from the radio through gamma rays. Although the leptonic scenario is able to explain most PWNe emission, a hadronic contribution cannot be excluded. A possible hadronic contribution to the high-energy gamma-ray emission inevitably leads to the production of neutrinos. Using 9.5 years of all-sky IceCube data, we report results from a stacking analysis to search for neutrino emission from 35 PWNe that are high-energy gamma-ray emitters. In the absence of any significant correlation, we set upper limits on the total neutrino emission from those PWNe and constraints on hadronic spectral components.

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

We report a quasi-periodic pulsation (QPP) event simultaneously detected from the spatial displacements of coronal loop at both EUV images and microwave emission during the preflare phase of a C1.1 flare on 2016 March 23. Using the motion magnification technique, a low-amplitude transverse oscillation with the growing period is discovered in a diffuse coronal loop in Atmospheric Imaging Assembly (AIA) image sequences at wavelength of 171 A, and the initial oscillation period is estimated to be ~397 s with a slow growth rate of 0.045. At the same time, a QPP with growing periods from roughly 300 s to nearly 500 s is discovered in the microwave flux in the same active region. Based on the imaging observations measured at EUV wavelengths by the AIA and at microwave 17 GHz by Nobeyama Radioheliograph, the diffuse coronal loop and the microwave radiation source are found to be connected through a hot loop seen in AIA images at wavelength of 94 A. The growing period of the QPP should be related to the modulation of LRC-circuit oscillating process in a current-carrying plasma loop. The existence of electric currents may imply the non-potentialities in the source region during the preflare phase.

We report pulsar timing observations carried out in L-band with NTSC's 40-meter Haoping Radio Telescope (HRT), which was constructed in 2014. The observations were carried out using the pulsar machine we developed. Timing observations toward millisecond pulsar J0437-4715 obtains a timing residual (r.m.s) of 397ns in the time span of 284 days. And our observations successfully detected Crab pulsar's glitch that happened on July 23rd, 2019.

In this work, we study how the dust coagulation/fragmentation will influence the evolution and observational appearances of vortices induced by a massive planet embedded in a low viscosity disk by performing global 2D high-resolution hydrodynamical simulations. Within the vortex, due to its higher gas surface density and steeper pressure gradients, dust coagulation, fragmentation and drift (to the vortex center) are all quite efficient, producing dust particles ranging from micron to $\sim 1.0\ {\rm cm}$, as well as overall high dust-to-gas ratio (above unity). In addition, the dust size distribution is quite non-uniform inside the vortex, with the mass weighted average dust size at the vortex center ($\sim 4.0$ mm) being a factor of $\sim10$ larger than other vortex regions. Both large ($\sim$ mm) and small (tens of micron) particles contribute strongly to affect the gas motion within the vortex. As such, we find that the inclusion of dust coagulation has a significant impact on the vortex lifetime and the typical vortex lifetime is about 1000 orbits. After the initial gaseous vortex is destroyed, the dust spreads into a ring with a few remaining smaller gaseous vortices with a high dust concentration and a large maximum size ($\sim$ mm). At late time, the synthetic dust continuum images for the coagulation case show as a ring inlaid with several hot spots at 1.33 mm band, while only distinct hot spots remain at 7.0 mm.

ART$^2$ is a 3D multi-wavelength Monte Carlo radiative transfer (RT) code that couples continuum and emission lines to track the propagation of photons and their interactions with the interstellar medium (ISM). The previous version of ART$^2$, which included continuum and Ly$\alpha$ line, has been extensively applied to hydrodynamics simulations to study multi-band properties of galaxies and ISM. Here, we describe new implementations of non-local thermodynamic equilibrium RT of molecular and atomic fine structure emission lines, and the parallelization of the code using a number of novel methods. The new ART$^2$ can efficiently and self-consistently produce a full spectrum that includes both continuum and lines such as [CII], [NII], [OIII], Ly$\alpha$, and CO. These essential features, together with the multi-phase ISM model and the adaptive grid, make ART$^2$ a multi-purpose code to study multi-wavelength properties of a wide range of astrophysical systems from planetary disks to large-scale structures. To demonstrate the capability of the new ART$^2$, we applied it to two hydrodynamics simulations: the zoom-in Milky Way Simulation to obtain panchromatic properties of individual galaxies, and the large-scale IllustrisTNG100 Simulation to obtain global properties such as the line intensity mappings. These products are essential for a broad array of studies such as the correlations between physical and panchromatic properties and their evolution. By enabling direct comparison between numerical simulations and multi-band observations, ART$^2$ provides a crucial theoretical framework for the interpretations of existing observations, the plan for future surveys, and the synergy between multi-band galaxy surveys and line intensity mappings. Therefore, ART$^2$ is a powerful and versatile tool to bridge the gap between theories and observations of cosmic structures.

We run three-dimensional numerical simulations for the accretion flow around the white dwarf (WD) in the progenitor system of Tycho's supernova (SN). The mass of the WD, mass of the companion star, and the orbital period are set to be 1$M_{\odot}$, 1.6$M_{\odot}$, and 0.794 day, respectively, based on theoretical and observational researches of Tycho's SN remnant (SNR). We find that when the magnetic field in the accreted material is negligible, outflowing wind is concentrated near the equatorial plane. When the magnetic field has energy equipartition with internal energy, polar wind is comparable with the equatorial wind. A carefully chosen magnetic field between the above two cases ($B=5.44\times10^3 \rm{G}$) can roughly reproduce the latitude-dependent wind required to form the peculiar periphery of Tycho's SNR.

GW190425 is the second neutron star merger event detected by the Advanced LIGO/Virgo detectors. If interpreted as a double neutron star merger, the total gravitational mass is substantially larger than that of the binary systems identified in the Galaxy. In this work we analyze the gravitational-wave data within the neutron star$-$black hole merger scenario. For the black hole, we yield a mass of $2.40^{+0.36}_{-0.32}M_\odot$ and an aligned spin of $0.141^{+0.067}_{-0.064}$. As for the neutron star we find a mass of $1.15^{+0.15}_{-0.13}M_\odot$ and the dimensionless tidal deformability of $1.4^{+3.8}_{-1.2}\times 10^{3}$. These parameter ranges are for 90\% credibility. The inferred masses of the neutron star and the black hole are not in tension with current observations and we suggest that GW190425 is a viable candidate of a neutron star$-$black hole merger event. Benefitting from the continual enhancement of the sensitivities of the advanced gravitational detectors and the increase of the number of the observatories, similar events are anticipated to be much more precisely measured in the future and the presence of black holes below the so-called mass gap will be unambiguously clarified. If confirmed, the mergers of neutron stars with (quickly rotating) low-mass black holes are likely important production sites of the heaviest r-process elements.

Monitoring chromospheric and photospheric indexes of magnetic activity can provide valuable information, especially the interaction between different parts of the atmosphere and their response to magnetic fields. We extract chromospheric indexes, S and Rhk+, for 59,816 stars from LAMOST spectra in the LAMOST-Kepler program, and photospheric index, Reff, for 5575 stars from Kepler light curves. The log Reff shows positive correlation with log Rhk+. We estimate the power-law indexes between Reff and Rhk+ for F-, G-, and K-type stars, respectively. We also confirm the dependence of both chromospheric and photospheric activity on stellar rotation. Ca II H and K emissions and photospheric variations generally decrease with increasing rotation periods for stars with rotation periods exceeding a few days. The power-law indexes in exponential decay regimes show different characteristics in the two activity-rotation relations. The updated largest sample including the activity proxies and reported rotation periods provides more information to understand the magnetic activity for cool stars.

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

We develop an estimator to extract the mean radial velocity profile of cosmic voids via the kinematic Sunyaev-Zeldovich effect of pairs of galaxies surrounding them. The estimator is tested with simulated pure kSZ map and void catalogue data from the same simulation. The results show that the recovered signal could be attenuated by low angular resolution of the map or large aperture photometry filter radius size, but the mean radial velocity profile can be fully recovered with our estimator. By applying the estimator to the Planck 2D-ILC CMB map, with galaxy and void catalogues from BOSS DR12, we find that the estimated void velocity profile is $3.31\sigma$ apart from null detection for for voids with continuously rising density profiles asymptoting to the mean density; and $1.75\sigma$ for voids with positive density contrast shell surrounded. By fitting the reconstructed data to the theoretical profile, we find the reduced $\chi^{2}$ to be $1.19$ and $0.62$ for the two types of void, respectively, indicating a good fit of the model to the data. We then forecast the detectability of the radial velocity profile of cosmic voids with future CMB surveys, including SPT-3G, AdvACT, and Simons Observatory. We find that the contamination effect from CMB residuals is negligible with survey area over $2,000~{\rm deg}^2$, especially with aperture photometry size less than $1\,{\rm arcmin}$. But the effect from instrumental noise is non-negligible. For future SPT-3G, AdvACT and Simons Observatory, the detection is potentially achievable from $3\sigma$ to $10\sigma$ C.L., depending on specific instrumental parameters. This opens a new window of probing dynamics of the cosmic structures from the kinematic Sunyaev-Zeldovich effect.

The circumgalactic medium (CGM) is closely linked to galaxy formation and evolution, but difficult to characterize observationally and typically poorly resolved in cosmological simulations. We use spherically-symmetric, idealized, high-resolution simulations of the CGM in $10^{12}M_\odot$ and $10^{11}M_\odot$ dark matter halos to characterize the gas pressure, turbulent and radial velocities, and degree of thermal and effective dynamic pressure support in the overall CGM as well as in its high- and low-temperature phases. We find that the $10^{12}M_\odot$ halo contains a CGM mostly formed of a hot gas halo in hydrostatic equilibrium out of which cold gas condenses and falls onto the central galaxy, while the $10^{11}M_\odot$ halo's CGM is not in hydrostatic equilibrium, has a wider spread of properties at a given galactocentric radius, does not have a clear separation of hot and cold phases, and is dominated by bulk motions. We also find that the degree of pressure support in the $10^{11}M_\odot$ halo is strongly dependent on the parameters of the galactic winds of the central galaxy. These results promote the idea that there is no "average" CGM and care must be taken when setting the initial conditions for a small-box simulation of a patch of the CGM.

We report observations of a transient source \fermi\ from radio to \grs. \fermi\ was discovered by the {\it Fermi-LAT} in May 2017. Follow-up {\it Swift-XRT} observations revealed three flaring episodes through March 2018, and the peak X-ray flux is about $10^3$ higher than the {\it ROSAT all-sky survey (RASS)} flux upper limit. Optical spectral measurements taken by the {\it Magellan 6.5-m telescope} and the {\it Lick-Shane telescope} both show a largely featureless spectrum, strengthening the BL Lac interpretation first proposed by \citet{Bruni18}. The optical and mid-infrared (MIR) emission goes to a higher state in 2018, when the flux in high energies goes down to a lower level. Our {\it RATAN-600m} measurements at 4.8~GHz and 8.2~GHz do not indicate any significant radio flux variation over the monitoring seasons in 2017 and 2018, nor deviate from the archival {\it NVSS} flux level. During GeV flaring times, the spectrum is very hard ($\Gamma_\gamma\sim$1.7) in the GeV band and at times also very hard (($\Gamma_{\rm X}\lesssim2$) in the X-rays, similar to a high-synchrotron-peak (or even an extreme) BL Lac object, making \fermi\ a good target for ground-based {\it Cherenkov telescopes}.

Cosmic voids are biased tracers of the large-scale structure of the universe. Separate universe simulations (SUS) enable accurate measurements of this biasing relation by implementing the peak-background split (PBS). In this work, we apply the SUS technique to measure the void bias parameters. We confirm that the PBS argument works well for underdense tracers. The response of the void size distribution depends on the void radius. For voids larger (smaller) than the size at the peak of the distribution, the void abundance responds negatively (positively) to a long wavelength mode. The linear bias from the SUS is in good agreement with the cross power spectrum measurement on large scales. Using the SUS, we have detected the quadratic void bias for the first time in simulations. We find that $ b_2 $ is negative when the magnitude of $ b_1 $ is small, and that it becomes positive and increases rapidly when $ |b_1| $ increases. We compare the results from voids identified in the halo density field with those from the dark matter distribution, and find that the results are qualitatively similar, but the biases generally shift to the larger voids sizes.

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.

The formation histories of lithium-rich and carbon-rich red giants are not yet understood. It has been proposed that the merger of a helium-core white dwarf with a red giant branch star might provide a solution. We have computed an extended grid of post-merger evolution models and combined these with predictions of binary star population synthesis. The results strongly support the proposal that the merger of a helium white dwarf with a red giant branch star can provide the progenitors of both lithium-rich red clump stars and early-R carbon stars. The distribution of post-merger models in $T_{\rm eff}$, $\log g$, $\log L$, the surface abundances of lithium and carbon, and the predicted space densities agree well with the observed distributions of these parameters for Li-rich and early-R stars in the Galaxy.

We perform a search for gravitational waves (GWs) from several supermassive binary black hole (SMBBH) candidates (NGC 5548, Mrk 231, OJ 287, PG 1302-102, NGC 4151, Ark 120 and 3C 66B) in long-term timing observations of the pulsar PSR J1909$-$3744 obtained using the Parkes radio telescope. No statistically significant signals were found. We constrain the chirp masses of those SMBBH candidates and find the chirp mass of NGC 5548 and 3C 66B to be less than $2.4 \times 10^9\,\rm M_{\odot}$ and $2.5 \times 10^9\,\rm M_{\odot}$ (with 95% confidence), respectively. Our upper limits remain a factor of 3 to 370 above the likely chirp masses for these candidates as estimated from other approaches. The observations processed here provide upper limits on the GW strain amplitude that improve upon the results from the first Parkes Pulsar Timing Array data release by a factor of 2 to 7. We investigate how information about the orbital parameters can help improve the search sensitivity for individual SMBBH systems. Finally, we show that these limits are insensitive to uncertainties in the Solar System ephemeris model.

As one of the most interesting Seyfert 1 galaxies, PG 2130+099 has been the target of several reverberation mapping (RM) campaigns over the years. However, its measured broad H$\beta$ line responses have been inconsistent, with time lags of $\sim$200 days, $\sim$25 days, and $\sim$10 days being reported for different epochs while its optical luminosity changed no more than 40%. To investigate this issue, we conducted a new RM-campaign with homogenous and high cadence (about $\sim$3 days) for two years during 2017--2019 to measure the kinematics and structure of the ionized gas. We successfully detected time lags of broad H$\beta$, He II, He I, and Fe II lines with respect to the varying 5100\AA continuum, revealing a stratified structure that is likely virialized with Keplerian kinematics in the first year of observations, but an inflow kinematics of the broad-line region from the second year. With a central black hole mass of $0.97_{-0.18}^{+0.15}\times 10^7~M_{\odot}$, PG 2130+099 has an accretion rate of $10^{2.1\pm0.5}L_{\rm Edd}c^{-2}$, where $L_{\rm Edd}$ is the Eddington luminosity and $c$ is speed of light, implying that it is a super-Eddington accretor and likely possesses a slim, rather than thin, accretion disk. The fast changes of the ionization structures of the three broad lines remain puzzling.

We consider possible effects of neutrino electric charge (millicharge) and charge radius on the neutrino-atom interaction processes such as (i) atomic ionization by neutrino impact and (ii) coherent elastic neutrino-nucleus scattering. The bounds on the neutrino millicharge and charge radius that follow from, respectively, the GEMMA and COHERENT experiments are presented and discussed.

In this work we re-analyze the data of GW170817, the first binary neutron star (BNS) merger event, in two ways, including the parameterized equation of state (EoS) method and gravitational wave (GW) parameter estimation analysis. Three sets of data, including the combined LIGO/Virgo (HLV) data, the LIGO-Hanford (H1) data and the LIGO-Livingston (L1) data, have been analyzed, respectively. We confirm the bimodal probability distribution of the tidal deformability of GW170817 with the HLV data and resolve the different contributions of H1 and L1 detectors. Our simulation reveals a tendency of ``overestimating" the tidal parameter with the decreasing signal-to-noise ratio (SNR). Such a tendency is likely caused by the dependence of the result on the prior (which is not well understood and has been assumed to be widely-distributed) in the low SNR case. Though this effect is interesting, the slightly higher SNR of H1 than L1 at frequencies above $\sim 100$ Hz is most-likely not the main reason for the lower tidal parameter region favored by H1. In view of the large fluctuation of the expected tide parameter in the case of low SNR, the different probability distributions favored by the L1 and H1 data of GW170817 are reasonable. We have also explored the dependence of expected difference of SNRs of H1 and L1 detectors on the source locations in the O3 and full sensitivity runs. Such maps are helpful in evaluating the relative powers of individual detectors on measuring the tidal deformabilities for the new double neutron star merger events.

Tidal interactions between the embedded planets and their surrounding protoplanetary disks are often postulated to produce the observed complex dust substructures, including rings, gaps, and asymmetries. In this Letter, we explore the consequences of dust coagulation on the dust dynamics and ring morphology. Coagulation of dust grains leads to dust size growth which, under typical disk conditions, produces faster radial drifts, potentially threatening the dust ring formation. Utilizing 2D hydrodynamical simulations of protoplanetary disks which include a full treatment of dust coagulation, we find that if the planet does not open a gap quickly enough, the formation of an inner ring is impeded due to dust coagulation and subsequent radial drift. Furthermore, we find that a "buildup" of sub-mm sized grains often appears in the dust emission at the outer edge of the dust disk.

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.

An accurate prediction of the geo-neutrino signal from the crust serves as a necessary prerequisite in the determination of the geo-neutrino flux from the mantle. In this work we report the non-negligible effect associated to the exact three-flavor antineutrino survival probability in the calculation of the crustal geo-neutrino signal, which was usually approximated as a constant average in previous studies. A geo-neutrino signal underestimation of about 1-2 TNU is observed as a result of the oscillatory behaviour within the local crustal region extending for about 300 km from the experimental site. We also estimated that the Mikheyev-Smirnov-Wolfenstein matter oscillation is responsible for a $0.1\%$-$0.3\%$ increase of the local crustal signal, depending on the detector location. This work reminds that the exact oscillation possibility in matter should be considered for future prediction of the local crustal geo-neutrino signal.

Two bright X-ray transients were reported from the Chandra Deep Field South archival data, namely CDF-S XT1 and XT2. Whereas the nature of the former is not identified, the latter was suggested as an excellent candidate for a rapidly spinning magnetar born from a binary neutron star (BNS) merger. Here we propose a unified model to interpret both transients within the framework of the BNS merger magnetar model. According to our picture, CDF-S XT2 is observed from the "free zone" where the magnetar spindown powered X-ray emission escapes freely, whereas CDF-S XT1 originates from the "trapped zone" where the X-ray emission is initially blocked by the dynamical ejecta and becomes transparent after the ejecta is pushed to a distance where Thomson optical depth drops below unity. We fit the magnetar model to the light curves of both transients and derived consistent parameters for the two events, with magnetic field, initial spin period and X-ray emission efficiency being ($B_p=10^{16}\,G$, $P=1.2\,\rm ms$, $\eta = 0.001$) and ($B_p=10^{15.8}\,G$, $P=4.4\, \rm ms$, $\eta = 0.001$) for XT1 and XT2, respectively. The "isotropic equivalent" ejecta mass of XT1 is $M_{\rm ej} \sim 10^{-3}$ $M_{\odot}$, while it is not constrained for XT2. Our results suggest that more extreme magnetar parameters are required to have XT1 detected from the trapped zone. The model parameters for both events are generally consistent with those derived from SGRB X-ray plateau observations. The host galaxy properties of both transients are also consistent with those of SGRBs. The event rate densities of both XT1 and XT2 are consistent with that of BNS mergers.

A common problem in cosmology is to integrate the product of two or more spherical Bessel functions (sBFs) with different configuration-space arguments against the power spectrum or its square, weighted by powers of wavenumber. Naively computing them scales as $N_{\rm g}^{p+1}$ with $p$ the number of configuration space arguments and $N_{\rm g}$ the grid size, and they cannot be done with Fast Fourier Transforms (FFTs). Here we show that by rewriting the sBFs as sums of products of sine and cosine and then using the product to sum identities, these integrals can then be performed using 1-D FFTs with $N_{\rm g} \log N_{\rm g}$ scaling. This "rotation" method has the potential to accelerate significantly a number of calculations in cosmology, such as perturbation theory predictions of loop integrals, higher order correlation functions, and analytic templates for correlation function covariance matrices. We implement this approach numerically both in a free-standing, publicly-available \textsc{Python} code and within the larger, publicly-available package \texttt{mcfit}. The rotation method evaluated with direct integrations already offers a factor of 6-10$\times$ speed-up over the naive approach in our test cases. Using FFTs, which the rotation method enables, then further improves this to a speed-up of $\sim$$1000-3000\times$ over the naive approach. The rotation method should be useful in light of upcoming large datasets such as DESI or LSST. In analysing these datasets recomputation of these integrals a substantial number of times, for instance to update perturbation theory predictions or covariance matrices as the input linear power spectrum is changed, will be one piece in a Monte Carlo Markov Chain cosmological parameter search: thus the overall savings from our method should be significant.

Pre-Big-Bang models in string cosmology predict a relic background of gravitational wave radiation in the early universe. The spectrum of this background shows that the energy density rises rapidly with frequency, which is an interesting target for high-frequency (i.e., kilohertz) detectors. In this paper, we discussed the constraining power of multiple configurations of current and future gravitational wave detector (GWD) networks to the stochastic background predicted in string cosmology. The constraining power is jointly determined by the overlap reduction function and the sensitivity curves of multiple detectors. And we further elaborated on the possible contribution of a future Chinese detector and a kilohertz detector to the constraining power of detector network for stochastic background in string cosmology. Our results show that the detectability of the GWD network for the string cosmology gravitational wave background will improve considerably with the joining of a Chinese detector. This is because a Chinese detector (e.g. located at Wuhan ) together with KAGRA, has a better overlap reduction function than the laser interferometer gravitational wave observatory detector pair, and therefore lead to more stringent limits for stochastic background detection. And with ideal overlap reduction function, namely, colocated detectors, a kilohertz sensitivity curve has better performance than previous detectors for stochastic background detection. Finally, the results are compared with the limitations given by the observational constraint of the Big Bang nucleosynthesis bound.

In this work we parameterize the Equation of State of dense neutron star (NS) matter with four pressure parameters of $\{\hat{p}_1, \hat{p}_2, \hat{p}_3, \hat{p}_4\}$ and then set the combined constraints with the data of GW 170817 and the data of 6 Low Mass X-ray Binaries (LMXBs) with thermonuclear burst or alternatively the symmetry energy of the nuclear interaction. We find that the nuclear data effectively narrow down the possible range of $\hat{p}_1$, the gravitational wave data plays the leading role in bounding $\hat{p}_2$, and the LMXB data as well as the lower bound on maximal gravitational mass of non-rotating NSs govern the constraints on $\hat{p}_3$ and $\hat{p}_4$. Using posterior samples of pressure parameters and some universal relations, we further investigate how the current data sets can advance our understanding of tidal deformability ($\Lambda$), moment of inertia ($I$) and binding energy ($BE$) of NSs. For a canonical mass of $1.4M_\odot$, we have $I_{1.4} = {1.43}^{+0.30}_{-0.13} \times 10^{38}~{\rm kg \cdot m^2}$, $\Lambda_{1.4} = 390_{-210}^{+280}$ , $R_{1.4} = 11.8_{-0.7}^{+1.2}~{\rm km}$ and $BE_{1.4} = {0.16}^{+0.01}_{-0.02} M_{\odot}$ if the constraints from the nuclear data and the gravitational wave data have been jointly applied. For the joint analysis of gravitational wave data and the LMXB data, we have $I_{1.4} = {1.28}^{+0.15}_{-0.08} \times 10^{38}~{\rm kg \cdot m^2}$, $\Lambda_{1.4} = 220_{-90}^{+90}$, $R_{1.4} = 11.1_{-0.6}^{+0.7}~{\rm km}$ and $BE_{1.4} = {0.18}^{+0.01}_{-0.01} M_{\odot}$. These results suggest that the current constraints on $\Lambda$ and $R$ still suffer from significant systematic uncertainties, while $I_{1.4}$ and $BE_{1.4}$ are better constrained.

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

A sample of Cataclysmic Variables (CVs) is presented including spectroscopically identified 380 spectra of 245 objects, of which 58 CV candidates are new discoveries. The BaggingTopPush and the Random Forest algorithms are applied to the Fifth Data Release (DR5) of LAMOST to retrieve CVs with strong emission lines and with broad absorption lines respectively. Based on spectroscopic classification, 134 dwarf novae, 41 nova-like variables and 19 magnetic CVs are identified from the sample. In addition, 89 high--inclination systems and 33 CVs showing companion stars are recognized and discussed for their distinct spectral characteristics. Comparisons between CVs from LAMOST and from published catalogs are made in spatial and magnitude distribution, and the difference of their locus in Gaia color--absolute magnitude diagram (CaMD) are also investigated. More interestingly, for two dwarf novae observed through LAMOST and SDSS in different epoch, their spectra both in quiescence phase and during outburst are exhibited.

We propose a revision of the system developed by L'epine et al. (2007) for spectroscopic M subdwarf classification. Based on an analysis of subdwarf spectra and templates from Savcheva et al. (2014), we show thatthe CaH1 feature originally proposed by Gizis (1997) is important in selecting reliable cool subdwarf spectra. This index should be used in combination with the [TiO5, CaH2+CaH3] relation provided by L\'epine et al. (2007) to avoid misclassification results. In the new system, the dwarf-subdwarf separators are first derived from a sample of more than 80,000 M dwarfs and a "labeled" subdwarf subsample, these objects being all visually identified from their optical spectra. Based on these two samples, we re-fit the initial [TiO5, CaH1] relation, and propose a new [CaOH, CaH1] relation supplementing the [TiO5, CaH1] relation to reduce the impact of uncertainty in flux calibration on classification accuracy. In addition, we recalibrate the $\zeta_{TiO/CaH}$ parameter defined in L'epine et al. (2007) to enable its successful application to LAMOST spectra. Using this new system, we select candidates from LAMOST Data Release 4 and finally identify a set of 2791 new M subdwarf stars, covering the spectral sequence from type M0 to M7. This sample contains a large number of objects located at low Galactic latitudes, especially in the Galactic anti-center direction, expanding beyond previously published halo- and thick disk-dominated samples. Besides, we detect magnetic activity in 141 objects. We present a catalog for this M subdwarf sample, including radial velocities, spectral indices and errors, activity flags, with a compilation of external data (photometric and GAIA DR2 astrometric parameters). The catalog is provided on-line, and the spectra can be retrieved from the LAMOST Data Release web portal.

We report our photometric and spectroscopic follow-up observations of the superoutbursts of three dwarf novae (GWAC\,180415A, GWAC\,181017A and GWAC\,181211A) identified independently by the Ground Wide-angle Cameras system, one of the ground-based instruments of the China-France SVOM mission. Based on a combination of our photometry and that taken from the AAVSO, our period analysis of the superhumps enables us to determine the mass ratios to be 0.0967-0.1163, 0.1879-0.1883 and 0.0981-0.1173 for GWAC\,180415A, GWAC\,181017A and GWAC\,181211A, respectively. GWAC\,180415A can be firmly identified as a WZ sge-type dwarf novae due to its long duration ($\sim2$ weeks) multiple rebrightenings with amplitudes of 3-4 magnitudes, the early superhump associated with a double-wave modulation and the low mass ratio. The inferred low mass ratio and location in the $\varepsilon-P_{\mathrm{orb}}$ diagram suggest that GWAC\,181211A is a WZ sge-type dwarf novae candidate. The measured Balmer decrements suggest the Balmer line emission is produced from an optical thick region in GWAC\,180415A and GWAC\,181017A, and from an optical thin region in GWAC\,181211A.

We computed a grid of theoretical models to fit the $\delta$ Scuti frequencies of CoRoT 100866999 detected earlier from the CoRoT timeserials. The pulsating primary star is determined to be a main sequence star with a rotation period of $4.1^{+0.6}_{-0.5}$ days, rotating slower than the orbital motion. The fundamental parameters of the primary star are determined to be $M$ = $1.71^{+0.13}_{-0.04}$ $M_{\odot}$, $Z=0.012^{+0.004}_{-0.000}$, $f_{\rm ov}$ = $0.02^{+0.00}_{-0.02}$, $T_{\rm eff}$ = $8024^{+249}_{-297}$ K, $L$ = $11.898^{+2.156}_{-1.847}$ $L_{\odot}$, $\log g$ = $4.166^{+0.013}_{-0.002}$, $R$ = $1.787^{+0.040}_{-0.016}$ $R_{\odot}$, and $X_{\rm c}$ = 0.488$^{+0.011}_{-0.020}$, matching well those obtained from the eclipsing light curve analysis. Based on the model fittings, $p_1$ and $p_5$ are suggested to be two dipole modes, and $p_3$, $p_4$, $p_6$, and $p_7$ to be four quadrupole modes. In particular, $p_4$ and $p_7$ are identified as two components of one quintuplet. Based on the best-fitting model, we find that $p_1$ is a g mode and the other nonradial modes have pronounced mixed characters, which give strong constraints on the convective core. Finally, the relative size of the convective core of CoRoT 100866999 is determined to $R_{\rm conv}/R$ = $0.0931^{+0.0003}_{-0.0013}$.

Supermassive black holes (SMBHs) are thought to provide energy that prevents catastrophic cooling in the centers of massive galaxies and galaxy clusters. However, it remains unclear how this "feedback" process operates. We use high-resolution optical data to study the kinematics of multi-phase filamentary structures by measuring the velocity structure function (VSF) of the filaments over a wide range of scales in the centers of three nearby galaxy clusters: Perseus, Abell 2597 and Virgo. We find that the motions of the filaments are turbulent in all three clusters studied. There is a clear correlation between features of the VSFs and the sizes of bubbles inflated by SMBH driven jets. Our study demonstrates that SMBHs are the main driver of turbulent gas motions in the centers of galaxy clusters and suggests that this turbulence is an important channel for coupling feedback to the environment. Our measured amplitude of turbulence is in good agreement with Hitomi Doppler line broadening measurement and X-ray surface brightness fluctuation analysis, suggesting that the motion of the cold filaments is well-coupled to that of the hot gas. The smallest scales we probe are comparable to the mean free path in the intracluster medium (ICM). Our direct detection of turbulence on these scales provides the clearest evidence to date that isotropic viscosity is suppressed in the weakly-collisional, magnetized intracluster plasma.

We developed a spectroscopic monitoring project to investigate the kinematics of the broad-line region (BLR) in active galactic nuclei (AGN) with ultra-fast outflows (UFOs). Mrk~79 is a radio-quiet AGN with UFOs and warm absorbers, had been monitored by three reverberation mapping (RM) campaigns, but its BLR kinematics is not understood yet. In this paper, we report the results from a new RM-campaign of Mrk~79, which was undertaken by Lijiang 2.4-m telescope. Mrk~79 is seeming to come out the faint state, the mean flux approximates a magnitude fainter than historical record. We successfully measured the lags of the broad emission lines including H$\beta~\lambda4861$, H$\gamma~\lambda4340$, He II $\lambda4686$ and He I $\lambda5876$ with respect to the varying AGN continuum. Based on the broad H$\beta~\lambda4861$ line, we measured black hole (BH) mass of $M_{\bullet}=5.13^{+1.57}_{-1.55}\times10^{7}M_{\odot}$, estimated accretion rates of ${\dot{M}_{\bullet}}=(0.05\pm0.02)~L_{\rm Edd}~c^{-2}$, indicating that Mrk~79 is a sub-Eddington accretor. We found that Mrk~79 deviates from the canonical Radius$-$Luminosity relationship. The marginal blueshift of the broad He II $\lambda4686$ line detected from rms spectrum indicates outflow of high-ionization gas. The velocity-resolved lag profiles of the broad H$\gamma~\lambda4340$, H$\beta~\lambda4861$, and He I $\lambda5876$ lines show similar signatures that the largest lag occurs in the red wing of the lines then the lag decreases to both sides. These signatures should suggest that the BLR of Keplerian motion probably exists the outflow gas motion. All findings including UFOs, warm absorbers, and the kinematics of high- and low-ionization BLR, may provide an indirect evidence that the BLR of Mrk~79 probably originates from disk wind.

We report circular polarization measurements from the first two years of observation with the 40 GHz polarimeter of the Cosmology Large Angular Scale Surveyor (CLASS). CLASS is conducting a multi-frequency survey covering 75% of the sky from the Atacama Desert designed to measure the cosmic microwave background (CMB) linear E and B polarization on angular scales $1^\circ \lesssim \theta \leq 90^\circ$, corresponding to a multipole range of $2 \leq \ell \lesssim 200$. The modulation technology enabling measurements of linear polarization at the largest angular scales from the ground, the Variable-delay Polarization Modulator, is uniquely designed to provide explicit sensitivity to circular polarization (Stokes $V$). We present a first detection of circularly polarized atmospheric emission at 40 GHz that is well described by a dipole with an amplitude of $124\pm4\,\mathrm{\mu K}$ when observed at an elevation of $45^\circ$, and discuss its potential impact as a foreground to CMB experiments. Filtering the atmospheric component, CLASS places a 95% C.L. upper limit of $0.4\,\mathrm{\mu K}^2$ to $13.5\,\mathrm{\mu K}^2$ on $\ell(\ell+1)C_\ell^{VV}/(2\pi)$ between $1 \leq \ell \leq 120$, representing a two-orders-of-magnitude improvement over previous limits.

We present a comprehensive study of a series of recurrent jets that occurred at the periphery of the NOAA active region 12114 on 2014 July 7. These jets were found to share the same source region and exhibited rotational motions as they propagated outward. The multi-wavelength imaging observations made by the AIA and {\it IRIS} telescopes reveal that some of the jets contain cool plasma only, while some others contain not only cool but also hot plasma. The Doppler velocities calculated from the {\it IRIS} spectra show a continuous evolution from blue to red shifts as the jet motions change from upward to downward. Additionally, some jets exhibit opposite Doppler shifts on their both sides, indicative of rotating motions along their axes. The inclination angle and three-dimensional velocity of the largest jet were inferred from the imaging and spectroscopic observations, which show a high consistence with those derived from the stereoscopic analysis using dual-perspective observations by {\it SDO}/AIA and {\it STEREO}-B/EUVI. By relating the jets to the local UV/EUV and full-disk {\it GOES} X-ray emission enhancements, we found that the previous five small-scale jets were triggered by five bright points while the last/largest one was triggered by a C1.6 solar flare. Together with a number of type III radio bursts generated during the jet eruptions as well as a weak CME that was observed in association with the last jet, our observations provide evidences in support of multi-scale magnetic reconnection processes being responsible for the production of jet events.

Given observations of the standard candles and the cosmic chronometers, we apply Pad\'{e} parameterization to the comoving distance and the Hubble paramter to find how stringent the constraint is set to the curvature parameter by the data. A weak informative prior is introduced in the modeling process to keep the inference away from the singularities. Bayesian evidence for different order of Pad\'{e} parameterizations is evaluated during the inference to select the most suitable parameterization in light of the data. The data we used prefer a parameterization form of comoving distance as $D_{01}(z)=\frac{a_0 z}{1+b_1 z}$ as well as a competitive form $D_{02}(z)=\frac{a_0 z}{1+b_1 z + b_2 z^2}$. Similar constraints on the spatial curvature parameter are established by those models and given the Hubble constant as a byproduct: $\Omega_k = 0.25^{+0.14}_{-0.13}$ (68\% confidence level [C.L.]), $H_0 = 67.7 \pm 2.0$ km/s/Mpc (68\% C.L.) for $D_{01}$, and $\Omega_k = -0.01 \pm 0.13$ (68\% C.L.), $H_0 = 68.8 \pm 2.0$ km/s/Mpc (68\% C.L.) for $D_{02}$. The evidence of different models demonstrates the qualitative analysis of the Pad\'{e} parameterizations for the comoving distance.

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.

We study the long-term variability in the optical monitoring database of Ark~120, a nearby radio-quiet active galactic nucleus (AGN) at a distance of 143 Mpc (z=0.03271). We compiled the historical archival photometric and spectroscopic data since 1974 and conducted a new two-year monitoring campaign in 2015-2017, resulting in a total temporal baseline over four decades. The long-term variations in the optical continuum exhibit a wave-like pattern and the Hbeta integrated flux series varies with a similar behavior. The broad Hbeta profiles have asymmetric double peaks, which change strongly with time and tend to merge into a single peak during some epochs. The period in the optical continuum determined from various period-search methods is about 20 yr and the estimated false alarm probability with null hypothesis simulations is about 1*10^-3. The overall variations of the broad Hbeta profiles also follow the same period. However, the present database only covers two cycles of the suggested period, which strongly encourages continued monitoring to track more cycles and confirm the periodicity. Nevertheless, in light of the possible periodicity and the complicated Hbeta profile, Ark~120 is one candidate of the nearest radio-quiet AGNs with possible periodic variability, and it is thereby a potential candidate host for a sub-parsec supermassive black hole binary.

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

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

This work uses a combination of a variational auto-encoder and generative adversarial network to compare different dark energy models in light of observations, e.g., the distance modulus from type Ia supernovae. The network finds an analytical variational approximation to the true posterior of the latent parameters in the models, yielding consistent model comparison results with those derived by the standard Bayesian method, which suffers from a computationally expensive integral over the parameters in the product of the likelihood and the prior. The parallel computational nature of the network together with the stochastic gradient descent optimization technique leads to an efficient way to compare the physical models given a set of observations. The converged network also provides interpolation for a dataset, which is useful for data reconstruction.

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

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

We present 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 investigate the impact of a highly eccentric 10 $M_{\rm \oplus}$ (where $M_{\rm \oplus}$ is the Earth mass) planet embedded in a dusty protoplanetary disk on the dust dynamics and its observational implications. By carrying out high-resolution 2D gas and dust two-fluid hydrodynamical simulations, we find that the planet's orbit can be circularized at large radii. After the planet's orbit is circularized, partial gap opening and dust ring formation happen close to the planet's circularization radius, which can explain the observed gaps/rings at the outer region of disks. When the disk mass and viscosity become low, we find that an eccentric planet can even open gaps and produce dust rings close to the pericenter and apocenter radii before its circularization. This offers alternative scenarios for explaining the observed dust rings and gaps in protoplanetary disks. A lower disk viscosity is favored to produce brighter rings in observations. An eccentric planet can also potentially slow down the dust radial drift in the outer region of the disk when the disk viscosity is low ($\alpha \lesssim2\times10^{-4}$) and the circularization is faster than the dust radial drift.

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

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

We study black hole - host galaxy correlations, and the relation between the over-massiveness (the distance from the average $M_{BH}-\sigma$ relation) of super-massive black holes (SMBHs) and star formation histories of their host galaxies in the Illustris and TNG100 simulations. We find that both simulations are able to produce black hole scaling relations in general agreement with observations at $z=0$, but with noticeable discrepancies. Both simulations show an offset from the observations for the $M_{BH}-\sigma$ relation, and the relation between $M_{BH}$ and the Sersic index. The relation between $M_{BH}$ and stellar mass $M_*$ is tighter than the observations, especially for TNG100. For massive galaxies in both simulations, the hosts of over-massive SMBHs (those above the mean $M_{BH}-\sigma$ relation) tend to have larger Sersic indices and lower baryon conversion efficiency, suggesting a multidimensional link between SMBHs and properties of their hosts. In Illustris, the hosts of over-massive SMBHs have formed earlier and have lower present-day star formation rates, in qualitative agreement with the observations for massive galaxies with $\sigma>100 \rm km/s$. For low-mass galaxies, such a correlation still holds in Illustris but does not exist in the observed data. For TNG100, the correlation between SMBH over-massiveness and star formation history is much weaker. The hosts of over-massive SMBHs generally have consistently larger star formation rates throughout history. These galaxies have higher stellar mass as well, due to the strong $M_{BH}-M_*$ correlation. Our findings show that simulated SMBH scaling relations and correlations are sensitive to features in the modeling of SMBHs.

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.

The formation and evolution of local brightest cluster galaxies (BCGs) is investigated by determining the stellar populations and dynamics from the galaxy core, though the outskirts and into the intracluster light (ICL). Integral spectroscopy of 23 BCGs observed out to 4 r_e is collected and high signal-to-noise regions are identified. Stellar population synthesis codes are used to determine the age, metallicity, velocity, and velocity dispersion of stars within each region. The intracluster light (ICL) spectra are best modeled with populations that are younger and less metal-rich than those of the BCG cores. The average BCG core age of the sample is 13.3$\pm$ 2.8 Gyr and the average metallicity is [Fe/H] = 0.30 $\pm$ 0.09, whereas for the ICL the average age is 9.2$\pm$3.5 Gyr and the average metallicity is [Fe/H] = 0.18$\pm$0.16. The velocity dispersion profile is seen to be rising or flat in most of the sample (17/23), and those with rising values reach the value of the host cluster's velocity dispersion in several cases. The most extended BCGs are closest to the peak of the cluster's X-ray luminosity. The results are consistent with the idea that the BCG cores and inner regions formed quickly and long ago, with the outer regions and ICL forming more recently, and continuing to assemble through minor merging. Any recent star formation in the BCGs is a minor component, and is associated with the cluster cool core status.

The Quijote simulations are a set of 43100 full N-body simulations spanning more than 7000 cosmological models in the $\{\Omega_{\rm m}, \Omega_{\rm b}, h, n_s, \sigma_8, M_\nu, w \}$ hyperplane. At a single redshift the simulations contain more than 8.5 trillions of particles over a combined volume of 43100 $(h^{-1}{\rm Gpc})^3$. Billions of dark matter halos and cosmic voids have been identified in the simulations, whose runs required more than 35 million core hours. The Quijote simulations have been designed for two main purposes: 1) to quantify the information content on cosmological observables, and 2) to provide enough data to train machine learning algorithms. In this paper we describe the simulations and show a few of their applications. We also release the Petabyte of data generated, comprising hundreds of thousands of simulation snapshots at multiple redshifts, halo and void catalogs, together with millions of summary statistics such as power spectra, bispectra, correlation functions, marked power spectra, and estimated probability density functions.

The typical blazar S5 0716$+$714 is very interesting due to its rapid and large amplitude variability and high duty cycle of micro-variability in optical band. We analyze the observations in I, R and V bands obtained with the $1.0m$ telescope at Weihai observatory of Shandong University from 2011 to 2018. The model of synchrotron radiation from turbulent cells in a jet has been proposed as a mechanism for explaining micro-variability seen in blazar light curves. Parameters such as the sizes of turbulent cells, the enhanced particle densities, and the location of the turbulent cells in the jet can be studied using this model. The model predicts a time lag between variations as observed in different frequency bands. Automatic model fitting method for micro-variability is developed, and the fitting results of our multi-frequency micro-variability observations support the model. The results show that both the amplitude and duration of flares decomposed from the micro-variability light curves confirm to the log-normal distribution. The turbulent cell size is within the range of about 5 to 55 AU, and the time lags of the micro-variability flares between the I-R and R-V bands should be several minutes. The time lags obtained from the turbulence model are consistent with the fitting statistical results, and the time lags of flares are correlated with the time lags of the whole light curve.

The analogy of the host galaxy of the repeating fast radio burst (FRB) source FRB 121102 and those of long gamma-ray bursts (GRBs) and super-luminous supernovae (SLSNe) has led to the suggestion that young magnetars born in GRBs and SLSNe could be the central engine of repeating FRBs. We test such a hypothesis by performing dedicated observations of the remnants of six GRBs with evidence of having a magnetar central engine using the Arecibo telescope and the Robert C. Byrd Green Bank Telescope (GBT). A total of $\sim 20$ hrs of observations of these sources did not detect any FRB from these remnants. Under the assumptions that all these GRBs left behind a long-lived magnetar and that the bursting rate of FRB 121102 is typical for a magnetar FRB engine, we estimate a non-detection probability of $8.9\times10^{-6}$. Even though these non-detections cannot exclude the young magnetar model of FRBs, we place constraints on the burst rate and luminosity function of FRBs from these GRB targets.

We search for host galaxy candidates of nearby fast radio bursts (FRBs), FRB 180729.J1316+55, FRB 171020, FRB 171213, FRB 180810.J1159+83, and FRB 180814.J0422+73 (the second repeating FRB). We compare the absolute magnitudes and the expected host dispersion measure $\rm DM_{host}$ of these candidates with that of the first repeating FRB, FRB 121102, as well as those of long gamma ray bursts (LGRBs) and superluminous supernovae (SLSNe), the proposed progenitor systems of FRB 121102. We find that while the FRB 121102 host is consistent with those of LGRBs and SLSNe, the nearby FRB host candidates, at least for FRB 180729.J1316+55, FRB 171020, and FRB180814.J0422+73, either have a smaller $\rm DM_{host}$ or are fainter than FRB121102 host, as well as the hosts of LGRBs and SLSNe. In order to avoid the uncertainty in estimating $\rm DM_{host}$ due to the line-of-sight effect, we propose a galaxy-group-based method to estimate the electron density in the inter-galactic regions, and hence, $\rm DM_{IGM}$. The result strengthens our conclusion. We conclude that the host galaxy of FRB 121102 is atypical, and LGRBs and SLSNe are likely not the progenitor systems of at least most nearby FRB sources. {The recently reported two FRB hosts differ from the host of FRB 121102 and also the host candidates suggested in this paper. This is consistent with the conclusion of our paper and suggests that the FRB hosts are very diverse. }

We forecast the cosmological constraints of the neutral hydrogen (HI) intensity mapping (IM) technique with radio telescopes by assuming 1-year of observational time. The current and future radio telescopes we consider here are FAST (Five hundred meter Aperture Spherical Telescope), BINGO (Baryon acoustic oscillations In Neutral Gas Observations), and SKA-I (Square Kilometre Array phase I) single-dish experiment. We also forecast the combined constraints of the three radio telescopes with Planck. We find that, the $1 \sigma$ errors of $(w_{0}, w_{a})$ for BINGO, FAST and SKA-I with respect to the fiducial values are respectively, $(0.9293, 3.5792), (0.4083, 1.5878), (0.3158, 0.4622)$. This is equivalent to $(56.04\%, 55.64\%)$ and $(66.02\%, 87.09\%)$ improvements in constraining $(w_{0}, w_{a})$ for FAST and SKA-I relative to BINGO. Simulations further show that SKA-I will put more stringent constraints than both FAST and BINGO when each of the experiment is combined with Planck measurement. For the $9$ cosmological parameters in consideration, we find that, there is a trade-off between SKA-I and FAST in constraining cosmological parameters, with each experiment being more superior in constraining a particular set of parameters.

We study the coherence of the near-infrared and X-ray background fluctuations and the X-ray spectral properties of the sources producing it. We use data from multiple Spitzer and Chandra surveys, including the UDS/SXDF surveys, the Hubble Deep Field North, the EGS/AEGIS field, the Chandra Deep Field South and the COSMOS surveys, comprising $\sim$2275 Spitzer/IRAC hours and $\sim$~16 Ms of Chandra data collected over a total area of $\sim$~1~deg$^2$. We report an overall $\sim$5$\sigma$ detection of a cross-power signal on large angular scales $>$ 20$''$ between the 3.6 and 4.5\mum\ and the X-ray bands, with the IR vs [1-2] keV signal detected at 5.2$\sigma$. The [0.5-1] and [2-4] keV bands are correlated with the infrared wavelengths at a $\sim$1$-$3$\sigma$ significance level. The hardest X-ray band ([4-7] keV) alone is not significantly correlated with any infrared wavelengths due to poor photon and sampling statistics. We study the X-ray SED of the cross-power signal. We find that its shape is consistent with a variety of source populations of accreting compact objects, such as local unabsorbed AGNs or high-z absorbed sources. We cannot exclude that the excess fluctuations are produced by more than one population. Because of poor statistics, the current relatively broad photometric bands employed here do not allow distinguishing the exact nature of these compact objects or if a fraction of the fluctuations have instead a local origin.

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.

Recent observations have successfully detected UV or infrared flux from galaxies at the epoch of reionization. However, the origin of their radiative properties has not been fully understood yet. Combining cosmological hydrodynamic simulations and radiative transfer calculations, we present theoretical predictions of multi-wavelength radiative properties of the first galaxies at z=6-15. We find that most of the gas and dust are ejected from star-forming regions due to supernova (SN) feedback, which allows UV photons to escape. We show that the peak of SED rapidly shifts between UV and infrared wavelengths on a timescale of 100 Myr due to intermittent star formation and feedback. When dusty gas covers the star-forming regions, the galaxies become bright in the observed-frame sub-millimeter wavelengths. In addition, we find that the escape fraction of ionizing photons also changes between 1-40% at z>10. The mass fraction of HII region changes with the star formation history, resulting in the fluctuations of metal lines and Lyman-alpha line luminosities. In the starbursting phase of galaxies with the halo mass $\sim 10^{11}\,{\rm M_{\odot}}$ ($10^{12}\,{\rm M_{\odot}}$), the simulated galaxy has $L_{\rm [O\,III]} \sim 10^{42}\, (10^{43})\, {\rm erg\,s^{-1}}$, which is consistent with the observed star-forming galaxies at z>7. Our simulations suggest that deep [C II] observation with ALMA can trace the distribution of neutral gas extending over $\sim 20$ physical kpc. We also find that the luminosity ratio $L_{\rm [O\,III]}/L_{\rm [C\,II]}$ decreases with bolometric luminosity due to metal enrichment. Our simulations show that the combination of multi-wavelength observations by ALMA and JWST will be able to reveal the multi-phase ISM structure and the transition from starbursting to outflowing phases of high-z galaxies.

We report an experiment to test quantum interference, entanglement and nonlocality using two dissimilar photon sources, the Sun and a semiconductor quantum dot on the Earth, which are separated by 150 million kilometers. By making the otherwise vastly distinct photons indistinguishable all degrees of freedom, we observe time-resolved two-photon quantum interference with a raw visibility of 0.796(17), well above the 0.5 classical limit, providing the first evidence of quantum nature of thermal light. Further, using the photons with no common history, we demonstrate post-selected two-photon entanglement with a state fidelity of 0.826(24), and a violation of Bell's inequality by 2.20(6). The experiment can be further extended to a larger scale using photons from distant stars, and open a new route to quantum optics experiments at an astronomical scale.

Matter evolved under influence of gravity from minuscule density fluctuations. Non-perturbative structure formed hierarchically over all scales, and developed non-Gaussian features in the Universe, known as the Cosmic Web. To fully understand the structure formation of the Universe is one of the holy grails of modern astrophysics. Astrophysicists survey large volumes of the Universe and employ a large ensemble of computer simulations to compare with the observed data in order to extract the full information of our own Universe. However, to evolve trillions of galaxies over billions of years even with the simplest physics is a daunting task. We build a deep neural network, the Deep Density Displacement Model (hereafter D$^3$M), to predict the non-linear structure formation of the Universe from simple linear perturbation theory. Our extensive analysis, demonstrates that D$^3$M outperforms the second order perturbation theory (hereafter 2LPT), the commonly used fast approximate simulation method, in point-wise comparison, 2-point correlation, and 3-point correlation. We also show that D$^3$M is able to accurately extrapolate far beyond its training data, and predict structure formation for significantly different cosmological parameters. Our study proves, for the first time, that deep learning is a practical and accurate alternative to approximate simulations of the gravitational structure formation of the Universe.

We investigate the possible presence of neutral hydrogen (HI) in intergalactic filaments at very low redshift ($z\sim 0.08$), by stacking a set of 274,712 2dFGRS galaxy pairs over 21-cm maps obtained with dedicated observations conducted with the Parkes radio telescope, over a total sky area of approximately 1,300 square degrees covering two patches in the northern and in the southern Galactic hemispheres. The stacking is performed by combining local maps in which each pair is brought to a common reference frame; the resulting signal from the edge galaxies is then removed to extract the filament residual emission. We repeat the analysis on maps cleaned removing either 10 or 20 foreground modes in a principal component analysis. Our study does not reveal any clear HI excess in the considered filaments in either case; we determine upper limits on the total filament HI brightness temperature at $T_{\rm b} \lesssim 10.3 \,\mu\text{K}$ for the 10-mode and at $T_{\rm b} \lesssim 4.8 \,\mu\text{K}$ for the 20-mode removed maps at the 95% confidence level. These estimates translate into upper limits for the local filament HI density parameter, $\Omega_{\rm HI}^{\rm (f)} \lesssim 7.0\times10^{-5}$ and $\Omega_{\rm HI}^{\rm (f)} \lesssim 3.2\times10^{-5}$ respectively, and for the HI column density, $N_{\rm HI} \lesssim 4.6\times10^{15}\,\text{cm}^{-2}$ and $N_{\rm HI} \lesssim 2.1\times10^{15}\,\text{cm}^{-2}$ respectively. These column density constraints are consistent with previous detections of HI in the warm-hot intergalactic medium obtained observing broad Ly $\alpha$ absorption systems. The present work shows for the first time how such constraints can be achieved using the stacking of galaxy pairs on 21-cm maps.

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

The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021 and start a five year survey in 2022. SO has 287 collaborators from 12 countries and 53 institutions, including 85 students and 90 postdocs. The SO experiment in its currently funded form ('SO-Nominal') consists of three 0.4 m Small Aperture Telescopes (SATs) and one 6 m Large Aperture Telescope (LAT). Optimized for minimizing systematic errors in polarization measurements at large angular scales, the SATs will perform a deep, degree-scale survey of 10% of the sky to search for the signature of primordial gravitational waves. The LAT will survey 40% of the sky with arc-minute resolution. These observations will measure (or limit) the sum of neutrino masses, search for light relics, measure the early behavior of Dark Energy, and refine our understanding of the intergalactic medium, clusters and the role of feedback in galaxy formation. With up to ten times the sensitivity and five times the angular resolution of the Planck satellite, and roughly an order of magnitude increase in mapping speed over currently operating ("Stage 3") experiments, SO will measure the CMB temperature and polarization fluctuations to exquisite precision in six frequency bands from 27 to 280 GHz. SO will rapidly advance CMB science while informing the design of future observatories such as CMB-S4.

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.

SPB star HD 50230, in fact a hybrid B-type pulsator, has been observed by CoRoT least 137 days. A nearly equidistant period spacing pattern are found among 8 modes which are extracted from the oscillation spectrum with more than 500 frequencies. However, it is thought to be most likely accidental by Szewczuk et al. (2014). In the present work, we analyze the 8 modes in depth with the $\chi^2$-matching method. Based on the best fitting model (model MA), we find that they can be well explained as a sequences of consecutive dipolar $(l,~m)=(1,~0)$. The period discrepancies between observations and the best fitting model are within 100 s except for the outlier which is up to 300 s. Based on the calculated CMMs, we find that, for pure g-mode oscillations, the buoyancy radius $\Lambda_{0}$ can be precisely measured with the $\chi^2$-matching method between observations and calculations. It represents the "Propagation time" of the g-mode from stellar surface to center. It is of $\Lambda_{0}=245.78\pm0.59~\mu$Hz with a precision of 0.24\%. In addition, we also find that HD 50230 is a metal-rich ($Z_{\rm init}=0.034-0.043$) star with a mass of $M=6.15-6.27~{\rm M_{\odot}}$. It is still located on hydrogen-burning phase with central hydrogen $X_{\rm C}=0.298-0.316$ (or $X_{\rm C}=0.306^{+0.010}_{-0.008}$), therefore has a convective core with a radius of $R_{\rm cc}=0.525-0.536~{\rm R_{\odot}}$ (or $R_{\rm cc}=0.531^{+0.005}_{-0.006}~{\rm R_{\odot}}$). In order to well interpret the structure of observed period spacing pattern, the convective core overshooting ($f_{\rm ov}=0.0175-0.0200$) and the extra diffusion mixing ($\log D_{\rm mix}=3.7-3.9$) should be taken into account in theoretical models.

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

Helioseismic observations have revealed many properties of the Sun: the depth and the helium abundance of the convection zone, the sound-speed and the density profiles in the solar interior. Those constraints have been used to judge the stellar evolution theory. With the old solar composition (e.g., GS98), the solar standard model is in reasonable agreement with the helioseismic constraints. However, a solar model with revised composition (e.g., AGSS09) with low abundance $Z$ of heavy elements cannot be consistent with those constraints. This is the so-called "solar abundance problem", standing for more than ten years even with the recent upward revised Ne abundance. Many mechanisms have been proposed to mitigate the problem. However, there is still not a low-$Z$ solar model satisfying all helioseismic constraints. In this paper, we report a possible solution to the solar abundance problem. With some extra physical processes that are not included in the standard model, solar models can be significantly improved. Our new solar models with convective overshoot, the solar wind, and an early mass accretion show consistency with helioseismic constraints, the solar Li abundance, and observations of solar neutrino fluxes.

In the present paper, we explore opportunities of applying the GRAVITY at the Very Large Telescope Interferometry (VLTI) with unprecedented spatial resolution to identify close binaries of supermassive black holes (CB-SMBHs) in active galactic nuclei (AGNs). Each SMBH is assumed to be separately surrounded by their own broad-line regions (BLRs) composed of clouds with virialized motion. Composition of the binary orbital motion and the virial motion of clouds in each BLR determines the projected velocity fields and hence differential phase curves, which are obviously different from that of a single BLR. We calculate emission line profiles and differential phase curves of CB-SMBHs for the GRAVITY. For the simplest case where angular momentums of two BLRs and orbital motion are parallel, a phase plateau generally appears in the phase curves. For other combinations of the angular momentum, the plateau is replaced by new peaks and valleys variously depending on the situations. Given a combination, phase curves are also sensitive to changes of parameters of CB-SMBHs. All these features are easily distinguished from the well-known $S$-shaped phase curves of a single BLR so that the GRAVITY is expected to reveal signals of CB-SMBH from candidates of AGNs. With joint analysis of observations of reverberation mapping campaigns, we can reliably identify CB-SMBHs, and measure their orbital parameters in the meanwhile. This independent measurement of the orbital parameters also has implications to analysis of Pulsar Timing Array (PTA) observations for properties of low-frequency gravitational waves in future.

Supermassive black hole feedback is thought to be responsible for the lack of star formation, or quiescence, in a significant fraction of galaxies. We explore how observable correlations between the specific star formation rate (sSFR), stellar mass (M$_{\rm{star}}$), and black hole mass (M$_{\rm{BH}}$) are sensitive to the physics of black hole feedback in a galaxy formation model. We use the IllustrisTNG simulation suite, specifically the TNG100 simulation and ten model variations that alter the parameters of the black hole model. Focusing on central galaxies at $z = 0$ with M$_{\rm{star}} > 10^{10}$ M$_{\odot}$, we find that the sSFR of galaxies in IllustrisTNG decreases once the energy from black hole kinetic winds at low accretion rates becomes larger than the gravitational binding energy of gas within the galaxy stellar radius. This occurs at a particular M$_{\rm{BH}}$ threshold above which galaxies are found to sharply transition from being mostly star-forming to mostly quiescent. As a result of this behavior, the fraction of quiescent galaxies as a function of M$_{\rm{star}}$ is sensitive to both the normalization of the M$_{\rm{BH}}$-M$_{\rm{star}}$ relation and the M$_{\rm{BH}}$ threshold for quiescence in IllustrisTNG. Finally, we compare these model results to observations of 91 central galaxies with dynamical M$_{\rm{BH}}$ measurements with the caveat that this sample is not representative of the whole galaxy population. While IllustrisTNG reproduces the observed trend that quiescent galaxies host more massive black holes, the observations exhibit a broader scatter in M$_{\rm{BH}}$ at a given M$_{\rm{star}}$ and show a smoother decline in sSFR with M$_{\rm{BH}}$.

Axion Like Particles (ALPs) are one promising kind of dark matter candidate particles that are predicted to couple with photons in the presence of magnetic fields. The oscillations between photons and ALPs travelling in the magnetic fields have been used to constrain ALP properties. In this work, we obtain some new constraints on the ALP mass $m_{\rm a}$ and the photon-ALP coupling constant $g$ with two different magnetic field models through TeV photons from PKS 2155-304. One is the discrete-$\varphi$ model that the magnetic field has the orientation angle $\varphi$ changes discretely and randomly from one coherent domain to the next, another is the linearly-continuous-$\varphi$ model that the magnetic field orientation angle $\varphi$ varies continuously across neighboring coherent domains. For the discrete-$\varphi$ model, we can obtain the best constraints on the ALP mass $m_{1}=m_{\rm a}/(1\ {\rm neV})=0.1$ and on the photon-ALP coupling constant $g_{11}=g/(10^{-11}\ {\rm GeV^{-1}})=5$, the reasonable range of the ALP mass $m_{1}$ is 0.08 $\thicksim$ 0.2 when $g_{11}$=5, and the only reasonable value of the photon-ALP coupling constant is $g_{11}$=5 when $m_{1}$=0.1. For the linearly-continuous-$\varphi$ model, we can obtain the best constraints on the ALP mass $m_{1}=0.1$ and on the photon-ALP coupling constant $g_{11}=0.7$, the reasonable range of the ALP mass $m_{1}$ is 0.05 $\thicksim$ 0.4 when $g_{11}$=0.7, and the reasonable range of the photon-ALP coupling constant $g_{11}$ is 0.5 $\thicksim$ 1 when $m_{1}$=0.1. All the results are consistent with the upper bound ($g<6.6\times10^{-11}\ {\rm GeV^{-1}}$, i.e. $g_{11}<6.6$) set by the CAST experiment.

The solar Ly$\alpha$ line is the strongest line in the ultraviolet waveband, and is greatly enhanced during solar flares. Here we present radiative hydrodynamic simulations of solar flares under different heating models, and calculate the response of this line taking into account non-equilibrium ionization of hydrogen and partial frequency redistribution. We find that in non-thermal heating models, the Ly$\alpha$ line can show a red or blue asymmetry corresponding to the chromospheric evaporation or condensation, respectively. The asymmetry may change from red to blue if the electron beam flux is large enough to produce a significant chromospheric condensation region. In the Ly$\alpha$ intensity lightcurve, there appears a dip when the change of asymmetry occurs. In thermal models, the Ly$\alpha$ line intensity peaks quickly and then falls, and the profile has an overall red asymmetry, which is similar to the profiles from heating by a soft electron beam. The Ly$\alpha$ profile shows a single red peak at the end of thermal heating, and the whole line is formed in a very small height range.

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

A multi-messenger, black hole (BH) - neutron star (NS) merger event still remains to be detected. The tidal (dynamical) ejecta from such an event, thought to produce a kinonova, is concentrated in the equatorial plane and occupies only part of the whole azimuthal angle. In addition, recent simulations suggest that the outflow or wind from the post-merger remnant disk, presumably anisotropic, can be a major ejecta component responsible for a kilonova. For any ejecta whose photosphere shape deviates from the spherical symmetry, the electron scattering at the photosphere causes a net polarization in the kilonova light. Recent observational and theoretical polarization studies have been focused to the NS-NS merger kilonova AT2017gfo. We extend those work to the case of a BH-NS merger kilonova. We show that the degree of polarization at the first $\sim 1$ hr can be up to $\sim$ 3\% if a small amount ($10^{-4} M_{\odot}$) of free neutrons have survived in the fastest component of the dynamical ejecta, whose beta-decay causes a precursor in the kilonova light. The polarization degree can be $\sim$ 0.6\% if free neutrons survived in the fastest component of the disk wind. Future polarization detection of a kilonova will constrain the morphology and composition of the dominant ejecta component, therefore help to identify the nature of the merger.

During the in-spiral stage of a compact binary, a wind bubble could be blown into interstellar medium, if the electromagnetic radiation due to the binary orbital motion is strong enough. Therefore, short gamma-ray bursts (SGRBs) due to double neutron star mergers would in principle happen in a wind bubble environment, which can influence the propagation of the SGRB jet and the consequent afterglow emission. By calculating the dynamics and synchrotron radiation of the jet-driven external shock, we reveal that an abrupt jump could appear in the afterglow light curves of SGRBs and the observational time of the jump is dependent on the viewing angle. This light curve jump provides an observational signature to constrain the radius of the wind bubble and then the power of the binary electromagnetic radiation, by combining with gravitational wave detection.

We present different signatures of chromospheric evaporation in two solar flares observed by the Interface Region Imaging Spectrograph (IRIS). In the B1.6 flare on 2016 December 6 (SOL2016-12-06T10:40), the transition region Si IV line and the chromospheric C II and Mg II lines show blueshifts with low velocities up to 20 km/s at the flare loop footpoints in the rise phase, indicative of a gentle chromospheric evaporation. While in the C1.6 flare on 2015 December 19 (SOL2015-12-19T10:51), the Si IV, C II, and Mg II lines exhibit redshifts with velocities from several to tens of km/s at the footpoints, which might suggest an explosive chromospheric evaporation. Explosive evaporation has been observed in many flares that were captured by IRIS; however, gentle evaporation, especially manifested as blueshifts in the cool Si IV, C II, and Mg II lines, has scarcely been reported. Our results bring some new insights into chromospheric evaporation in the IRIS era.

In recent years, the discovery of massive quasars at z~7 has provided a striking challenge to our understanding of the origin and growth of supermassive black holes in the early Universe. Mounting observational and theoretical evidence indicates the viability of massive seeds, formed by the collapse of supermassive stars, as a progenitor model for such early, massive accreting black holes. Although considerable progress has been made in our theoretical understanding, many questions remain regarding how (and how often) such objects may form, how they live and die, and how next generation observatories may yield new insight into the origin of these primordial titans. This review focusses on our present understanding of this remarkable formation scenario, based on discussions held at the Monash Prato Centre from November 20--24, 2017, during the workshop "Titans of the Early Universe: The Origin of the First Supermassive Black Holes."

We measure the bivariate luminosity function (BLF) of galaxy pairs and use it to probe and characterize the galaxy-galaxy interaction between pair members. The galaxy pair sample is selected from the main galaxy sample of Sloan Digital Sky Survey and supplied with a significant number of redshifts from the LAMOST spectral and GAMA surveys. We find the BLFs depend on the projected distance $d_{\text{p}}$ between pair members. At large separation $d_{\text{p}} > 150 h^{-1}\ \text{kpc}$, the BLF degenerates into a luminosity function (LF) of single galaxies, indicating few interactions between pair members. At $100 h^{-1}\ \text{kpc} \leq d_{\text{p}} \leq 150 h^{-1}\ \text{kpc}$, the BLF starts to show the correlation between pair members, in the sense that the shape of the conditional luminosity function (CLF) of one member galaxy starts to depend on the luminosity of the other member galaxy. Specifically, the CLF with a brighter companion has a steeper faint-end slope, which becomes even more significant at $50 h^{-1}\ \text{kpc} \leq d_{\text{p}} \leq 100 h^{-1}\ \text{kpc}$. This behavior is consistent with the scenario, \textit{and also is the observational evidence}, that dynamic friction drives massive major merger pairs to merge more quickly. At close distance $d_{\text{p}} \leq 50 h^{-1}\ \text{kpc}$, besides the merging time-scale effect, the BLF also shows an overall brightening of $\Delta M_r \geq 0.04$ mag, which reveals the enhanced star formation of the close-pair phase. By combining another statistical conclusion that the star formation rate of late-type galaxies in close pairs is enhanced at a level of about 40\%, we further conclude that the average starburst time-scale of close pairs is as long as 0.4 Gyr.

We compiled a sample of 73 active galactic nuclei (AGNs) with reverberation mapping (RM) observations from RM campaigns including our ongoing campaign of monitoring super-Eddington accreting massive black holes (SEAMBHs). This sample covers a large range of black hole (BH) mass $(M_{\bullet}=10^{6-9}~M_\odot)$, dimensionless accretion rates $(\dot{\mathscr{M}}=10^{-2.7}-10^{2.7})$ and 5100~\AA~luminosity $(L_{5100}=10^{42-46}~\rm erg~s^{-1})$, allowing us to systematically study the AGN variability and their relations with BH mass, accretion rates, and optical luminosity. We employed the damped random walk (DRW) model to delineate the optical variability of continuum at 5100~\AA~and obtained damped variability timescale ($\tau_{\rm d}$) and amplitude ($\sigma_{\rm d}$) using a Markov Chain Monte Carlo (MCMC) method. We also estimated the traditional variability amplitudes ($F_{\rm var}$), which provide a model-independent measure and therefore are used to test the DRW results. We found that AGN variability characteristics are generally correlated with $(M_{\bullet},\dot{\mathscr{M}},L_{5100})$. These correlations are smooth from sub-Eddington to super-Eddington accretion AGNs, probably implying that the AGN variability may be caused by the same physical mechanism.

The tomographic Alcock-Paczynski (AP) method can result in tight cosmological constraints by using small and intermediate clustering scales of the large scale structure (LSS) of the galaxy distribution. By focusing on the redshift dependence, the AP distortion can be distinguished from the distortions produced by the redshift space distortions (RSD). In this work, we combine the tomographic AP method with other recent observational datasets of SNIa+BAO+CMB+$H_0$ to reconstruct the dark energy equation-of-state $w$ in a non-parametric form. The result favors a dynamical DE at $z\lesssim1$, and shows a mild deviation ($\lesssim2\sigma$) from $w=-1$ at $z=0.5-0.7$. We find the addition of the AP method improves the low redshift ($z\lesssim0.7$) constraint by $\sim50\%$.

Using the Purple Mountain Observatory Delingha 13.7 m telescope, we conducted a large-scale $^{12}$CO $\left(J=1 \rightarrow 0\right)$ outflow survey (over $\sim$ 110 deg$^2$) toward the W3/4/5 complex and its surroundings. In total, 459 outflow candidates were identified. Approximately 62\% (284) were located in the Perseus arm, including W3/4/5 complex and its surroundings, while $\sim$ 35\% (162) were located in the Local arm, $\sim$ 1\% (5) in the Outer arm, and $\sim$ 2\% (8) in two interarm regions. This result indicated that star formation was concentrated in the Galactic spiral arms. The detailed spatial distribution of the outflow candidates showed that the Perseus arm presented the most active star formation among the study regions. The W3/4/5 complex is a great region to research massive star formation in a triggered environment. A key region, which has been well-studied by other researches, is in the eastern high-density W3 complex that neighbors the W4 complex. Two shell-like structures in the Local arm contain candidates that can be used to study the impact on star formation imposed by massive or intermediate-mass stars in relatively isolated systems. The majority of outflow candidates in the two interarm regions and the Outer arm are located in filamentary structures.

The growth of solids from sub-micron to millimeter and centimeter sizes is the early step toward the formation of planets inside protoplanetary disks (PPDs). However, such processes and their potential impact on the later stages of solid growth are still poorly understood. In this work, we test the hypothesis that most disks contain at least one ringed structure with a relative small radius. We have carried out a large family of 1D two-fluid (gas+dust) hydrodynamical simulations by evolving the gas and dust motion self-consistently while allowing dust size to evolve via coagulation and fragmentation. We investigate the joint effects of ringed structures and dust size growth on the overall sub-millimeter and millimeter (mm) flux and spectral index of PPDs. Ringed structures slow down the dust radial drift and speed up the dust growth. In particular, we find that those unresolved disks with a high fragmentation velocity ($\sim10\ {\rm m~s^{-1}}$) and a high dust surface density ($\sim10\ {\rm g\ cm^{-2}}$ in the ring) can have mm spectral indices as low as $\sim2.0$, consistent with mm observations of faint disks in nearby star forming regions. Furthermore, disks with more than one ringed structure can potentially reproduce brighter disks with spectral indices lower than $\sim2.5$. Future multi-wavelength high-resolution observations of these low spectral index sources can be used to test the existence of the ringed structures in the unresolved disks and differentiate the effects of dust size growth from optical depth.

NGC 4477 is a low-mass lenticular galaxy in the Virgo Cluster, residing at 100\,kpc to the north of M87. Using a total of 116\,ks {\sl Chandra} observations, we study the interplay between its hot ($\sim$0.3\,keV) gas halo and the central supermassive black hole. A possible cool core is indicated by the short cooling time of the gas at the galaxy centre. We identify a pair of symmetric cavities lying 1.1\,kpc southeast and 0.9\,kpc northwest of the galaxy centre with diameters of 1.3\,kpc and 0.9\,kpc, respectively. We estimate that these cavities are newly formed with an age of $\sim$4\,Myr. No radio emission is detected at the positions of the cavities with the existing VLA data. The total energy required to produce the two cavities is $\sim$$10^{54}$\,erg, at least two orders of magnitude smaller than that of typical X-ray cavities. NGC 4477 is arguably far the smallest system and the only lenticular galaxy in which AGN X-ray cavities have been found. It falls on the scaling relation between the cavity power and the AGN radio luminosity, calibrated for groups and clusters. Our findings suggest that AGN feedback is universal among all cool core systems. Finally, we note the presence of molecular gas in NGC~4477 in the shape of a regular disk with ordered rotation, which may not be related to the feedback loop.

An increasing number of white dwarf stars show atmospheric chemical composition polluted by heavy elements accreted from debris disk material. The existence of such debris disks strongly suggests the presence of one or more planet(s) whose gravitational interaction with rocky planetesimals is responsible for their disruption by tidal effect. The ZZ Ceti pulsator and polluted DAZ white dwarf GD 133 is a good candidate for searching for such a potential planet. We started in 2011 a photometric follow-up of its pulsations. As a result of this work in progress, we used the data gathered from 2011 to 2015 to make an asteroseismological analysis of GD 133, providing the star parameters from a best fit model with $M$/$M_{\odot}$ = 0.630 $\pm$ 0.002, $T_{\rm eff}$ = 12400 K $\pm$ 70 K, log($M_{\rm He}/M$) = -2.00 $\pm$ 0.02, log($M_{\rm H}/M$) = -4.50 $\pm$ 0.02 and determining a rotation period of $\approx$ 7 days.

We report results of the first reverberation mapping campaign of I Zwicky 1 during $2014$-$2016$, which showed unambiguous reverberations of the broad H$\beta$ line emission to the varying optical continuum. From analysis using several methods, we obtain a reverberation lag of $\tau_{\rm H\beta}=37.2^{+4.5}_{-4.9}\,$ days. Taking a virial factor of $f_{_{\rm BLR}}=1$, we find a black hole mass of $M_{\bullet}=9.30_{-1.38}^{+1.26}\times 10^6 M_{\odot}$ from the mean spectra. The accretion rate is estimated to be $203.9_{-65.8}^{+61.0}\,L_{\rm Edd}c^{-2}$, suggesting a super-Eddington accretor, where $L_{\rm Edd}$ is the Eddington luminosity and $c$ is the speed of light. By decomposing {\it Hubble Space Telescope} images, we find that the stellar mass of the bulge of its host galaxy is $\log (M_{\rm bulge}/M_{\odot}) = \rm 10.92\pm 0.07$. This leads to a black hole to bulge mass ratio of $\sim 10^{-4}$, which is significantly smaller than that of classical bulges and elliptical galaxies. After subtracting the host contamination from the observed luminosity, we find that I Zw 1 follows the empirical $R_{\rm BLR}\propto L_{5100}^{1/2}$ relation.

Neutron star-neutron star mergers are known to be associated with short gamma-ray bursts. If the neutron star equation of state is sufficiently stiff, at least some of such mergers will leave behind a supramassive or even a stable neutron star that spins rapidly with a strong magnetic field (i.e., a magnetar). Such a magnetar signature may have been observed as the X-ray plateau following a good fraction (up to 50%) of short gamma-ray bursts, and it has been expected that one may observe short gamma-ray burst-less X-ray transients powered by double neutron star mergers. A fast X-ray transient (CDF-S XT1) was recently found to be associated with a faint host galaxy whose redshift is unknown. Its X-ray and host-galaxy properties allow several possibleexplanations including a short gamma-ray burst seen off axis, a low-luminosity gamma-ray burst at high redshift, or a tidal disruption event involving an intermediate mass black hole and a white dwarf. Here we report a second X-ray transient, CDF-S XT2, that is associated with a galaxy at redshift z = 0.738. The light curve is fully consistent with being powered by a millisecond magnetar. More intriguingly, CDF-S XT2 lies in the outskirts of its star-forming host galaxy with a moderate offset from the galaxy center, as short bursts often do. The estimated event rate density of similar X-ray transients, when corrected to the local value, is consistent with the double neutron star merger rate density inferred from the detection of GW170817.

We present observations of the C-band $1_{10}-1_{11}$ (4.8 GHz) and Ku-band $2_{11}-2_{12}$ (14.5 GHz) K-doublet lines of H$_2$CO and the C-band $1_{10}-1_{11}$ (4.6 GHz) line of H$_2$$^{13}$CO toward a large sample of Galactic molecular clouds, through the Shanghai Tianma 65-m radio telescope (TMRT). Our sample with 112 sources includes strong H$_2$CO sources from the TMRT molecular line survey at C-band and other known H$_2$CO sources. All three lines are detected toward 38 objects (43 radial velocity components) yielding a detection rate of 34\%. Complementary observations of their continuum emission at both C- and Ku-bands were performed. Combining spectral line parameters and continuum data, we calculate the column densities, the optical depths and the isotope ratio H$_2$$^{12}$CO/H$_2$$^{13}$CO for each source. To evaluate photon trapping caused by sometimes significant opacities in the main isotopologue's rotational mm-wave lines connecting our measured K-doublets, and to obtain $^{12}$C/$^{13}$C abundance ratios, we used the RADEX non-LTE model accounting for radiative transfer effects. This implied the use of the new collision rates from \citet{Wiesenfeld2013}. Also implementing distance values from trigonometric parallax measurements for our sources, we obtain a linear fit of $^{12}$C/$^{13}$C = (5.08$\pm$1.10)D$_{GC}$ + (11.86$\pm$6.60), with a correlation coefficient of 0.58. D$_{GC}$ refers to Galactocentric distances. Our $^{12}$C/$^{13}$C ratios agree very well with the ones deduced from CN and C$^{18}$O but are lower than those previously reported on the basis of H$_2$CO, tending to suggest that the bulk of the H$_2$CO in our sources was formed on dust grain mantles and not in the gas phase.

The tilt, rotation, or offset of each CCD with respect to the focal plane, as well as the distortion of the focal plane itself, cause shape distortions to the observed objects, an effect typically known as field distortion (FD). We point out that FD provides a unique way of quantifying the accuracy of cosmic shear measurement. The idea is to stack the shear estimators from galaxies that share similar FD-induced shape distortions. Given that the latter can be calculated with parameters from astrometric calibrations, the accuracy of the shear estimator can be directly tested on real images. It provides a way to calibrate the multiplicative and additive shear recovery biases within the scientific data itself, without requiring simulations or any external data sets. We use the CFHTLenS images to demonstrate the accuracy of the Fourier_Quad shear recovery method. We highlight some details in our image processing pipeline, including background removal, source identification and deblending, astrometric calibration, star selection for PSF reconstruction, noise reduction, etc.. We show that in the shear ranges of -0.005 < g_1 < 0.005 and -0.008 < g_2 < 0.008, the multiplicative biases are at the level of < 0.04. Slight additive biases on the order of 5E-4 (6 sigma) are identified for sources provided by the official CFHTLenS catalog (not using its shear catalog), but are minor (4 sigma) for source catalog generated by our Fourier_Quad pipeline.

LAMOST DR5 released more than 200,000 low resolution spectra of early-type stars with S/N>50. Searching for metallic-line (Am) stars in such a large database and study of their statistical properties are presented in this paper. Six machine learning algorithms were experimented with using known Am spectra, and both the empirical criteria method(Hou et al. 2015) and the MKCLASS package(Gray et al. 2016) were also investigated. Comparing their performance, the random forest (RF) algorithm won, not only because RF has high successful rate but also it can derives and ranks features. Then the RF was applied to the early type stars of DR5, and 15,269 Am candidates were picked out. Manual identification was conducted based on the spectral features derived from the RF algorithm and verified by experts. After manual identification, 9,372 Am stars and 1,131 Ap candidates were compiled into a catalog. Statistical studies were conducted including temperature distribution, space distribution, and infrared photometry. The spectral types of Am stars are mainly between F0 and A4 with a peak around A7, which is similar to previous works. With the Gaia distances, we calculated the vertical height Z from the Galactic plane for each Am star. The distribution of Z suggests that the incidence rate of Am stars shows a descending gradient with increasing jZj. On the other hand, Am stars do not show a noteworthy pattern in the infrared band. As wavelength gets longer, the infrared excess of Am stars decreases, until little or no excess in W1 and W2 bands.

During the last decade, numerous and varied observations, along with increasingly sophisticated numerical simulations, have awakened astronomers to the central role the circumgalactic medium (CGM) plays in regulating galaxy evolution. It contains the majority of the baryonic matter associated with a galaxy, along with most of the metals, and must continually replenish the star forming gas in galaxies that continue to sustain star formation. And while the CGM is complex, containing gas ranging over orders of magnitude in temperature and density, a simple emergent property may be governing its structure and role. Observations increasingly suggest that the ambient CGM pressure cannot exceed the limit at which cold clouds start to condense out and precipitate toward the center of the potential well. If feedback fueled by those clouds then heats the CGM and causes it to expand, the pressure will drop and the "rain" will diminish. Such a feedback loop tends to suspend the CGM at the threshold pressure for precipitation. The coming decade will offer many opportunities to test this potentially fundamental principle of galaxy evolution.

Understanding the processes that drive galaxy formation and shape the observed properties of galaxies is one of the most interesting and challenging frontier problems of modern astrophysics. We now know that the evolution of galaxies is critically shaped by the energy injection from accreting supermassive black holes (SMBHs). However, it is unclear how exactly the physics of this feedback process affects galaxy formation and evolution. In particular, a major challenge is unraveling how the energy released near the SMBHs is distributed over nine orders of magnitude in distance throughout galaxies and their immediate environments. The best place to study the impact of SMBH feedback is in the hot atmospheres of massive galaxies, groups, and galaxy clusters, which host the most massive black holes in the Universe, and where we can directly image the impact of black holes on their surroundings. We identify critical questions and potential measurements that will likely transform our understanding of the physics of SMBH feedback and how it shapes galaxies, through detailed measurements of (i) the thermodynamic and velocity fluctuations in the intracluster medium (ICM) as well as (ii) the composition of the bubbles inflated by SMBHs in the centers of galaxy clusters, and their influence on the cluster gas and galaxy growth, using the next generation of high spectral and spatial resolution X-ray and microwave telescopes.

Recent observations of giant ellipticals and brightest cluster galaxies (BCGs) provide tentative evidence for a correlation between the luminosity of the H$\alpha$ emitting gas filaments and the strength of feedback associated with the active galactic nucleus (AGN). Motivated by this, we use 3D radiation-hydrodynamic simulations with the code Enzo to examine and quantify the relationship between the observable properties of the H$\alpha$ filaments and the kinetic and radiative feedback from supermassive black holes in BCGs. We find that the spatial extent and total mass of the filaments show positive correlations with AGN feedback power and can therefore be used as probes of the AGN activity. We also examine the relationship between the AGN feedback power and velocity dispersion of the H$\alpha$ filaments and find that the kinetic luminosity shows statistically significant correlation with the component of the velocity dispersion along the jet axis but not the components perpendicular to it.

The model-independent reconstruction of the energy spectra of $\overline{\nu}^{}_e$, $\nu^{}_e$ and $\nu^{}_x$ (i.e., $\nu^{}_\mu$, $\nu^{}_\tau$ and their antiparticles) from the future observation of a galactic core-collapse supernova (SN) is of crucial importance to understand the microscopic physics of SN explosions. To this end, we propose a practically useful method to combine the multi-channel detection of SN neutrinos in a large liquid-scintillator detector (e.g., JUNO), namely, the inverse beta decay $\overline{\nu}^{}_e + p \to e^+ + n$, the elastic neutrino-proton scattering $\nu + p \to \nu + p$ and the elastic neutrino-electron scattering $\nu + e^- \to \nu + e^-$, and reconstruct the energy spectra of $\overline{\nu}^{}_e$, $\nu^{}_e$ and $\nu^{}_x$ by making the best use of the observational data in those three channels. In addition, the neutrino energy spectra from the numerical simulations of the delayed neutrino-driven SN explosions are implemented to demonstrate the robustness of our method. Taking the ordinary matter effects into account, we also show how to extract the initial neutrino energy spectra in the presence of neutrino flavor conversions.

Despite many decades of study, the kinematics of the broad-line region of 3C~273 are still poorly understood. We report a new, high signal-to-noise, reverberation mapping campaign carried out from November 2008 to March 2018 that allows the determination of time lags between emission lines and the variable continuum with high precision. The time lag of variations in H$\beta$ relative to those of the 5100 Angstrom continuum is $146.8_{-12.1}^{+8.3}$ days in the rest frame, which agrees very well with the Paschen-$\alpha$ region measured by the GRAVITY at The Very Large Telescope Interferometer. The time lag of the H$\gamma$ emission line is found to be nearly the same as for H$\beta$. The lag of the Fe II emission is $322.0_{-57.9}^{+55.5}$ days, longer by a factor of $\sim$2 than that of the Balmer lines. The velocity-resolved lag measurements of the H$\beta$ line show a complex structure which can be possibly explained by a rotation-dominated disk with some inflowing radial velocity in the H$\beta$-emitting region. Taking the virial factor of $f_{\rm BLR} = 1.3$, we derive a BH mass of $M_{\bullet} = 4.1_{-0.4}^{+0.3} \times 10^8 M_{\odot}$ and an accretion rate of $9.3\,L_{\rm Edd}\,c^{-2}$ from the H$\beta$ line. The decomposition of its $HST$ images yields a host stellar mass of $M_* = 10^{11.3 \pm 0.7} M_\odot$, and a ratio of $M_{\bullet}/M_*\approx 2.0\times 10^{-3}$ in agreement with the Magorrian relation. In the near future, it is expected to compare the geometrically-thick BLR discovered by the GRAVITY in 3C 273 with its spatially-resolved torus in order to understand the potential connection between the BLR and the torus.

We develop methods for investigating baryon acoustic oscillation (BAO) features in cosmological models with non-trivial (but slowly varying) averaged spatial curvature: models that are not necessarily flat, close to flat, nor with constant spatial curvature. The class of models to which our methods apply include Lemaitre-Tolman-Bondi models, modified gravity cosmologies, and inhomogeneous cosmologies with backreaction - in which we do not have a prediction of the shape of the spatial 2-point correlation function, but where we nevertheless expect to see a BAO feature in the present-day galaxy distribution, in form of an excess in the galaxy 2-point correlation function. We apply our methods to the Baryon Oscillation Spectroscopic Survey (BOSS) dataset, investigating both the Lambda Cold Dark Matter ($\Lambda$CDM) and timescape cosmological models as case studies. The correlation functions measured in the two fiducial models contain a similarly-pronounced BAO feature. We use the relative tangential and radial BAO scales to measure the anisotropic Alcock-Paczy\'nski distortion parameter, $\epsilon$, which is independent of the underlying BAO preferred scale. We find that $\epsilon$ is consistent with zero in both fiducial cosmologies, indicating that models with a different spatial curvature behaviour can account for the relative positions of the tangential and radial BAO scale. We validate our methods using $\Lambda$CDM mocks.

The Simons Observatory (SO) is a new cosmic microwave background experiment being built on Cerro Toco in Chile, due to begin observations in the early 2020s. We describe the scientific goals of the experiment, motivate the design, and forecast its performance. SO will measure the temperature and polarization anisotropy of the cosmic microwave background in six frequency bands: 27, 39, 93, 145, 225 and 280 GHz. The initial configuration of SO will have three small-aperture 0.5-m telescopes (SATs) and one large-aperture 6-m telescope (LAT), with a total of 60,000 cryogenic bolometers. Our key science goals are to characterize the primordial perturbations, measure the number of relativistic species and the mass of neutrinos, test for deviations from a cosmological constant, improve our understanding of galaxy evolution, and constrain the duration of reionization. The SATs will target the largest angular scales observable from Chile, mapping ~10% of the sky to a white noise level of 2 $\mu$K-arcmin in combined 93 and 145 GHz bands, to measure the primordial tensor-to-scalar ratio, $r$, at a target level of $\sigma(r)=0.003$. The LAT will map ~40% of the sky at arcminute angular resolution to an expected white noise level of 6 $\mu$K-arcmin in combined 93 and 145 GHz bands, overlapping with the majority of the LSST sky region and partially with DESI. With up to an order of magnitude lower polarization noise than maps from the Planck satellite, the high-resolution sky maps will constrain cosmological parameters derived from the damping tail, gravitational lensing of the microwave background, the primordial bispectrum, and the thermal and kinematic Sunyaev-Zel'dovich effects, and will aid in delensing the large-angle polarization signal to measure the tensor-to-scalar ratio. The survey will also provide a legacy catalog of 16,000 galaxy clusters and more than 20,000 extragalactic sources.

The detection of the high-energy ($\sim290$ TeV) neutrino coincident with the flaring blazar TXS 0506+056, the first and only $3\sigma$ neutrino-source association to date, provides new, multimessenger tests of the weak equivalence principle (WEP) and Lorentz invariance. Assuming that the flight time difference between the TeV neutrino and gamma-ray photons from the blazar flare is mainly caused by the gravitational potential of the Laniakea supercluster of galaxies, we show that the deviation from the WEP for neutrinos and photons is conservatively constrained to have an accuracy of $10^{-6}-10^{-7}$, which is 3--4 orders of magnitude better than previous results placed by MeV neutrinos from supernova 1987A. In addition, we demonstrate that the association of the TeV neutrino with the blazar flare sets limits on the energy scales of quantum gravity for both linear and quadratic violations of Lorentz invariance (LIV) to $E_{\rm QG, 1}>3.2\times10^{15}-3.7\times10^{16}$ GeV and $E_{\rm QG, 2}>4.0\times10^{10}-1.4\times10^{11}$ GeV. These improve previous limits on both linear and quadratic LIV energy scales in neutrino propagation by 5--7 orders of magnitude.

Asymptotic Giant Branch (AGB) winds from evolved stars not only provide a non-trivial amount of mass and energy return, but also produce dust grains in massive elliptical galaxies. Due to the fast stellar velocity and the high ambient temperature, the wind is thought to form a comet-like tail, similar to Mira in the Local Bubble. Many massive elliptical galaxies and cluster central galaxies host extended dusty cold filaments. The fate of the cold dusty stellar wind and its relation to cold filaments are not well understood. In this work, we carry out both analytical and numerical studies of the interaction between an AGB wind and the surrounding hot gas. We find that the cooling time of the tail is inversely proportional to the ambient pressure. In the absence of cooling, or in low pressure environments (e.g., the outskirts of elliptical galaxies), AGB winds are quickly mixed into the hot gas, and all the AGB winds have similar appearance and head-to-tail ratio. In high pressure environments, such as the Local Bubble and the central regions of massive elliptical galaxies, some of the gas in the mixing layer between the stellar wind and the surrounding hot gas can cool efficiently and cause the tail to become longer. Our simulated tail of Mira itself has similar length and velocity to that observed, and appears similar to the simulated AGB tail in the central regions of massive galaxies. We speculate that instead of thermal instability, the induced condensation at the mixing layer of AGB winds may be the origin of cold filaments in massive galaxies and galaxy clusters. This naturally explains the existence of dust and PAH in the filaments.

We proposed a model-independent method to constrain the cosmological parameters using the Distance Sum Rule of the FLRW metric by combining the time delay distances and the comoving distances through a multi-messenger approach. The time delay distances are measured from lensed gravitational wave~(GW) signals together with their corresponding electromagnetic wave~(EM) counterpart, while the comoving distances are obtained from a parametrized fitting approach with independent supernova observations. With a series of simulations based on Einstein Telescope, Large Synoptic Survey Telescope and The Dark Energy Survey, we find that only 10 lensed GW+EM systems can achieve the constraining power comparable to and even stronger than 300 lensed quasar systems due to more precise time delay from lensed GW signals. Specifically, the cosmological parameters can be constrained to ~$k=0.01_{-0.05}^{+0.05}$ and ~$H_0=69.7_{-0.35}^{+0.35}$ (1$\sigma$). Our results show that more precise time delay measurements could provide more stringent cosmological parameter values, and lensed GW+EM systems therefore can be applied as a powerful tool in the future precision cosmology.

We present a catalog of stellar age and mass estimates for a sample of 640\,986 red giant branch (RGB) stars of the Galactic disk from the LAMOST Galactic Spectroscopic Survey (DR4). The RGB stars are distinguished from the red clump stars utilizing period spacing derived from the spectra with a machine learning method based on kernel principal component analysis (KPCA). Cross-validation suggests our method is capable of distinguishing RC from RGB stars with only 2 per cent contamination rate for stars with signal-to-noise ratio (SNR) higher than 50. The age and mass of these RGB stars are determined from their LAMOST spectra with KPCA method by taking the LAMOST - $Kepler$ giant stars having asteroseismic parameters and the LAMOST-TGAS sub-giant stars based on isochrones as training sets. Examinations suggest that the age and mass estimates of our RGB sample stars with SNR $>$ 30 have a median error of 30 per cent and 10 per cent, respectively. Stellar ages are found to exhibit positive vertical and negative radial gradients across the disk, and the age structure of the disk is strongly flared across the whole disk of $6<R<13$\,kpc. The data set demonstrates good correlations among stellar age, [Fe/H] and [$\alpha$/Fe]. There are two separate sequences in the [Fe/H] -- [$\alpha$/Fe] plane: a high--$\alpha$ sequence with stars older than $\sim$\,8\,Gyr and a low--$\alpha$ sequence composed of stars with ages covering the whole range of possible ages of stars. We also examine relations between age and kinematic parameters derived from the Gaia DR2 parallax and proper motions. Both the median value and dispersion of the orbital eccentricity are found to increase with age. The vertical angular momentum is found to fairly smoothly decrease with age from 2 to 12\,Gyr, with a rate of about $-$50\,kpc\,km\,s$^{-1}$\,Gyr$^{-1}$. A full table of the catalog is public available online.

A grid of theoretical models are computed to fit the 9 oscillation modes of IP Per detected earlier from a multi-site ground-based campaign. Fitting results show that there are two sets of theoretical models that could reproduce the observed frequencies more or less euqally well. In view of other available spectroscopic and photometric measurements, our best fitting stellar parameters for IP Per are $\upsilon_{\rm e} = 91^{+5}_{-3}$ km s$^{-1}$, $Z=0.009^{+0.004}_{-0.001}$, $M = 1.64^{+0.10}_{-0.04}$ $M_{\odot}$, $T_{\rm eff} = 7766^{+348}_{-173}$ K, $\log L/L_{\odot} = 1.125^{+0.094}_{-0.046}$, $\log g = 4.041^{+0.008}_{-0.003}$ dex, $R = 2.022^{+0.042}_{-0.018}$ $R_{\odot}$, $\tau_{0} = 8711^{+68}_{-35}$ s, age = $7.39^{+0.76}_{-0.46}$ Myr. Meanwhile, IP Per is found to be a pre-main sequence star where CN cycle has not yet reached the equilibrium state. At present fourteen percent of C12 have been turned into N14. Based on the best-fitting model, $f_{6}$ is identified as a radial mode, $f_{1}$ and $f_{2}$ as two dipole modes, and $f_{3}$, $f_{4}$, $f_{5}$, $f_{7}$, $f_{8}$, and $f_{9}$ as six quadrupole modes.

Optimal analyses using the 2-point functions of large-scale structure probes require accurate covariance matrices. A covariance matrix of the 2-point function comprises the disconnected part and the connected part. While the connected covariance only becomes important on small scales, the disconnected covariance is dominant on large scales, where the survey window has a significant impact. In this work, we develop an analytical method to compute the disconnected covariance, accounting for the window effect. Derived under the flat-sky approximation, our formalism is applicable to wide surveys by swapping in the curved-sky window functions. Our method works for both the power spectrum and the correlation function, and applies to the covariances of various probes including the multipoles and the wedges of 3D clustering, the angular and the projected statistics of clustering and shear, as well as the cross covariances between different probes. We verify the analytic covariance against the sample covariance from the galaxy mock simulations in two test cases: (1) the power spectrum multipole covariance, and (2) the joint covariance of the projected correlation function and the correlation function multipoles. Our method achieve good agreement with the mocks, while at a negligible computational cost. Unlike mocks, our analytic covariance is free of sampling noise, which often leads to numerical problems and the need to inflate the errors. In addition, our method can use the best-fit power spectrum as input, in contrast to the standard procedure of using a fiducial model that may deviate significantly from the truth. We also show that a naive diagonal power spectrum covariance underestimates the signal-to-noise ratio compared to our analytic covariance. The code that accompanies this paper is available at https://github.com/eelregit/covdisc.

We present the results of a search for galaxy clusters and groups in the $\sim2$ square degree of the COSMOS field using all available X-ray observations from the XMM-Newton and Chandra observatories. We reach an X-ray flux limit of $3\times10^{-16}\;ergs\;cm^{-2}\;s^{-1}$ in 0.5--2 keV range, and identify 247 X-ray groups with $M_{200c}=8\times10^{12}-3\times10^{14}\;M_{\odot}$ at a redshift range of $0.08\leq z<1.53$, using the multiband photometric redshift and the master spectroscopic redshift catalogues of the COSMOS. The X-ray centres of groups are determined using high-resolution Chandra imaging. We investigate the relations between the offset of the brightest group galaxies (BGGs) from halo X-ray centre and group properties and compare with predictions from semi-analytic models and hydrodynamical simulations. We find that BGG offset decreases with both increasing halo mass and decreasing redshift with no strong dependence on the X-ray flux and SNR. We show that the BGG offset decreases as a function of increasing magnitude gap with no considerable redshift dependent trend. The stellar mass of BGGs in observations extends over a wider dynamic range compared to model predictions. At $z<0.5$, the central dominant BGGs become more massive than those with large offsets by up to 0.3dex, in agreement with model prediction. The observed and predicted lognormal scatter in the stellar mass of both low- and large-offset BGGs at fixed halo mass is $\sim0.3$dex.

Many analyses have concluded that the accretion disc sizes measured from the microlensing variability of quasars are larger than the expectations from the standard thin disc theory by a factor of $\sim4$. We propose a simply model by invoking a strong wind from the disc to flatten its radial temperature profile, which can then reconcile the size discrepancy problem. This wind model has been successfully applied to several microlensed quasars with a wind strength $s\lesssim1.3$ by only considering the inward decreasing of the mass accretion rate (where $s$ is defined through $\dot{M}(R)\propto({R}/{R_{0}})^{s}$ ). After further incorporating the angular momentum transferred by the wind, our model can resolve the disc size problem with an even lower wind parameter. The corrected disc sizes under the wind model are correlated with black hole masses with a slope in agreement with our modified thin disc model.

Most results of the reverberation monitoring of active galaxies showed a universal scaling of the time delay of the Hbeta emission region with the monochromatic flux at 5100 A, with very small dipersion. Such a scaling favored the dust-based formation mechanism of the Broad Line Region (BLR). Recent reverberation measurements showed that actually a significant fraction of objects exhibits horter lags than the previously found scaling. Here we demonstrate that these shorter lags can be explained by the old concept of scaling of the BLR size with the ionization parameter. Assuming a universal value of this parameter and universal value of the cloud density reproduces the distribution of observational points in the time delay vs. monochromatic flux plane, provided that a range of black hole spins is allowed. However, a confirmation of the new measurements for low/moderate Eddington ratio sources is strongly needed before the dust-based origin of the BLR can be excluded.

This is the eighth in a series of papers reporting on a large reverberation mapping campaign to measure black hole (BH) mass in high accretion rate active galactic nuclei (AGNs). We employ the recently developed dynamical modeling approach for broad-line regions (BLRs) based on the method of Pancoast et al. to analyze the reverberation mapping dataset of Mrk 142 observed in the first monitoring season. In this approach, continuum variations are reconstructed using a damped random walk process, and BLR structure is delineated using a flexible disk-like geometry, in which BLR clouds move around the central BH with Keplerian orbits or inflow/outflow motion. The approach also includes the possibilities of anisotropic emission of BLR clouds, non-linear response of the line emission to the continuum, and different long-term trends in the continuum and emission-line variations. We implement the approach in a Bayesian framework that is apt for parallel computation and use a Markov Chain Monte Carlo technique to recover the parameters and uncertainties for the modeling, including mass of the central BH. We apply three BLR models with different prescriptions of BLR clouds distributions and find that the best model for fitting the data of Mrk 142 is a two-zone BLR model, consistent with the theoretical BLR model surrounding slim accretion disks. The best model yields a BH mass of $\log (M_\bullet/M_\odot)=6.23_{-0.45}^{+0.26}$, resulting in a virial factor of $\log f=-0.36_{-0.54}^{+0.33}$ for the full width at half maximum of the H$\beta$ line measured from the mean spectrum. The virial factors for the other measures of the H$\beta$ line width are also presented.

Overshooting and semiconvection are among the most uncertainties in the evolution of massive stars. Complete mixing over a certain distance beyond the convective boundary (Stothers \& Chin 1985) and an exponentially decaying diffusion outside the convection zone (Herwig 2000) are commonly adopted for the overshoot mixing. Recently, Li (2012, 2017) developed the $k$-$\omega$ model, which can be used in both convection zones and overshooting regions. We incorporated it in calculations of $30M_{\odot}$ stellar models. We find that in the main sequence stage, models with the $k$-$\omega$ model are almost identical to models with complete mixing in the overshooting region beyond the convective core, and the overshooting in the $k$-$\omega$ model is equivalent to an overshooting distance of about $0.15H_P$. In the post main sequence stage, we find that the overshooting below the bottom of the intermediate convection zone beyond the hydrogen-burning shell can significantly restrict the size of the hydrogen-depleted core, and can penetrate effectively into the hydrogen-burning shell. These two effects are crucial for the evolution of the core helium burning stage. During the core helium burning stage, we find that the overshooting model based on the $k$-$\omega$ model results in a similar complete mixing region but a much wider partial mixing region than the overshooting model based on Herwig (2000). In particular, the overshooting distance in the core helium burning stage may be significantly smaller than that in the main sequence phase for massive stars.

Sgr A* represents a unique laboratory for the detailed study of accretion processes around a low-luminosity supermassive black hole (SMBH). Recent X-ray observations have allowed for spatially resolved modeling of the emission from the accretion flow around the SMBH, placing tight constraints on the flux and spectral shape of the accretion from the inner region with $r<10^3~r_g$, where $r_g$ is the gravitational radius of the black hole. We present here the first modeling of the multi-band spectral energy distribution (SED) of this inner region to better constrain the physical condition of the innermost accretion flow. Our modeling uses the Markov chain Monte Carlo (MCMC) method to fit the SED, accounting for the limitations on the accretion rate at the outer radius of $10^3~r_g$ from the earlier works and the domination of the accretion flow within $30~r_g$ to the sub-mm bump. It is found that the fitting results of the outflow index could be very different. If only the most luminous part of the SED, the sub-mm bump, is considered, the outflow index is about 0, while if low-frequency radio data and X-ray data are also included, the outflow index could be 0.37 or even higher. The great difference of the fitting results indicates that the outflow index should be variable along radius, with a strong outflow in the outer region and a weak outflow in the innermost region. Such weak outflow agrees with numerical simulations and is possible to explain the multi-band SED even better.

We present an empirical stellar spectra library created using spectra from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) DR5. This library represents a uniform data set ranging from 3750 through 8500K in effective temperature (T_eff), from -2.5 through +1.0 dex in metallicity ([Fe/H]), and from 0 to 5.0 dex in gravity (log g). The spectra in the library have resolution R~1800, with well-calibrated fluxes and rest-framed wavelengths. Using a large number of red stars observed by LAMOST, we generated denser K type templates to fill in data missing from current empirical spectral libraries, particularly the late K type. For K giants, we calibrated the spectroscopic surface gravities against the asteroseismic surface gravities. To verify the reliability of the parameters labeled for this library, we performed an internal cross-validation by using a chi^2 minimization method to interpolate parameters of each individual spectrum using the remaining spectra in the library. We obtained precisions of 41 K, 0.11 dex, and 0.05 dex for T_eff, log g, and [Fe/H], respectively, which means the templates are labeled with correct stellar parameters. Through external comparisons, we confirmed that measurements of the stellar parameters through this library can achieve accuracies of approximately 125K in T_eff, 0.1 dex in [Fe/H] and 0.20 dex in log g without systematic offset. This empirical library is useful for stellar parameter measurements because it has large parameter coverage and full wavelength coverage from 3800 to 8900A.

The model of holographic dark energy in which dark energy interacts with dark matter is investigated in this paper. In particular, we consider the interacting holographic dark energy model in the context of a perturbed universe, which was never investigated in the literature. To avoid the large-scale instability problem in the interacting dark energy cosmology, we employ the generalized version of the parameterized post-Friedmann approach to treat the dark energy perturbations in the model. We use the current observational data to constrain the model. Since the cosmological perturbations are considered in the model, we can then employ the redshift-space distortions (RSD) measurements to constrain the model, in addition to the use of the measurements of expansion history, which was either never done in the literature. We find that, for both the cases with $Q=\beta H\rho_{\rm c}$ and $Q=\beta H_0\rho_{\rm c}$, the interacting holographic dark energy model is more favored by the current data, compared to the holographic dark energy model without interaction. It is also found that, with the help of the RSD data, a positive coupling $\beta$ can be detected at the $2.95\sigma$ statistical significance for the case of $Q=\beta H_0\rho_{\rm c}$.

We have started a long-term reverberation mapping project using the Wyoming Infrared Observatory 2.3 meter telescope titled "Monitoring AGNs with H\beta\ Asymmetry" (MAHA). The motivations of the project are to explore the geometry and kinematics of the gas responsible for complex H\beta\ emission-line profiles, ideally leading to an understanding of the structures and origins of the broad-line region (BLR). Furthermore, such a project provides the opportunity to search for evidence of close binary supermassive black holes. We describe MAHA and report initial results from our first campaign, from December 2016 to May 2017, highlighting velocity-resolved time lags for four AGNs with asymmetric H\beta\ lines. We find that 3C 120, Ark 120, and Mrk 6 display complex features different from the simple signatures expected for pure outflow, inflow, or a Keplerian disk. While three of the objects have been previously reverberation mapped, including velocity-resolved time lags in the cases of 3C 120 and Mrk 6, we report a time lag and corresponding black hole mass measurement for SBS 1518+593 for the first time. Furthermore, SBS 1518+593, the least asymmetric of the four, does show velocity-resolved time lags characteristic of a Keplerian disk or virialized motion more generally. Also, the velocity-resolved time lags of 3C 120 have significantly changed since previously observed, indicating an evolution of its BLR structure. Future analyses of the data for these objects and others in MAHA will explore the full diversity of H\beta\ lines and the physics of AGN BLRs.

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

We outline a full-scale search for galaxies exhibiting double-peaked profiles of promi- nent narrow emission lines, motivated by the prospect of finding objects related to merging galaxies, and even dual active galactic nuclei candidates as by-product, from the Large Sky Area Multi-object Fiber Spectroscopic Telescope (LAMOST) Data Re- lease 4. We assemble a large sample of 325 candidates with double-peaked or strong asymmetric narrow emission lines, with 33 objects therein appearing optically resolved dual-cored structures, close companions or signs of recent interaction on the Sloan Dig- ital Sky Survey images. A candidate from LAMOST (J074810.95+281349.2) is also stressed here based on the kinematic and spatial decompositions of the double-peaked narrow emission line target, with analysis from the cross-referenced Mapping Nearby Galaxies at the Apache Point Observatory (MaNGA) survey datacube. MaNGA en- ables us to constrain the origin of double peaks for these sources, and with the IFU data we infer that the most promising origin of double-peaked profiles for LAMOST J074810.95+281349.2 is the `Rotation Dominated + Disturbance' structure.

Recent observations have successfully detected UV-bright and infrared-bright galaxies in the epoch of reionization. However, the origin of their radiative properties has not been understood yet. Combining cosmological hydrodynamic simulations and radiative transfer calculations, we present predictions of multi-wavelength radiative properties of the first galaxies at $z\sim 6-15$. Using zoom-in initial conditions, we investigate three massive galaxies and their satellites in different environment and halo masses: $M_{\rm h} = 2.4\times 10^{10}{\rm M_{\odot}}$ (Halo-10), $1.6\times 10^{11}{\rm M_{\odot}}$ (Halo-11) and $0.7\times 10^{12}{\rm M_{\odot}}$ (Halo-12) at $z=6$. We find that most of gas and dust are ejected from star-forming regions by supernova feedback, which allows UV photons to escape. We show that the peak of the spectral energy distribution (SED) rapidly changes between UV and infrared wavelengths on a time-scale of $\sim$\,100\,Myrs due to intermittent star formation and feedback, and the escape fraction of UV photons fluctuates in the range of $0.2-0.8$ at $z<10$ with a time-averaged value of 0.3. When dusty gas covers the star-forming regions, the galaxies become bright in the observed-frame sub-millimeter wavelengths. We predict the detectability of high-$z$ galaxies with the Atacama Large Millimeter Array (ALMA). For a sensitivity limit of $0.1\,{\rm mJy}$ at $850\,{\rm \mu m}$, the detection probability of galaxies in halos $M_{\rm h} \gtrsim 10^{11}\,{\rm M_{\odot}}$ at $z\lesssim 7$ exceeds fifty per cent. We argue that supernova feedback can produce the observed diversity of SEDs for high-$z$ galaxies.

We develop a novel approach in exploring the joint dependence of halo bias on multiple halo properties using Gaussian process regression. Using a $\Lambda$CDM $N$-body simulation, we carry out a comprehensive study of the joint bias dependence on halo structure, formation history and environment. We show that the bias is a multivariate function of halo properties that falls into three regimes. For massive haloes, halo mass explains the majority of bias variation. For early-forming haloes, bias depends sensitively on the recent mass accretion history. For low-mass and late-forming haloes, bias depends more on the structure of a halo such as its shape and spin. Our framework enables us to convincingly prove that $V_\mathrm{max}/V_\mathrm{vir}$ is a lossy proxy of formation time for bias modelling, whereas the mass, spin, shape and formation time variables are non-redundant with respect to each other. Combining mass and formation time largely accounts for the mass accretion history dependence of bias. Combining all the internal halo properties fully accounts for the density profile dependence inside haloes, and predicts the clustering variation of individual haloes to a $20\%$ level at $\sim 10\mathrm{Mpc}h^{-1}$. When an environmental density is measured outside $1\mathrm{Mpc}h^{-1}$ from the halo centre, it outperforms and largely accounts for the bias dependence on the internal halo structure, explaining the bias variation above a level of $30\%$.

The enormous sizes and variability of emission of radio-loud AGNs arise from the relativistic flows of plasma along two oppositely directed jets. We use the Athena hydrodynamics code to simulate an extensive suite of 54 propagating three-dimensional relativistic jets with wide ranges of input jet velocities and jet-to-ambient matter density ratios. We determine which parameter sets yield unstable jets that produce jet dominated FR I type radio galaxy morphologies and which tend to produce stable jets with hot-spots and FR II morphologies. Nearly all our simulations involve jets with internal pressures matched to those of the ambient medium but we also consider over-pressured jets and discuss differences from the standard ones. We also show that the results are not strongly dependent on the adiabatic index of the fluid. We focus on simulations that remain stable for extended distances (60-240) times the initial jet radius. Scaled to the much smaller sizes probed by VLBI observations, the fluctuations in such simulated flows yield variability in observed emissivity on timescales from months. Adopting results for the densities, pressures and velocities from these simulations we estimate normalized rest frame synchrotron emissivities from individual cells in the jets. The observed emission from each cell is strongly dependent upon its variable Doppler boosting factor. We sum the fluxes from thousands of zones around the primary reconfinement shock. The light curves an power-spectra, with red-noise slopes between -2.1 and -2.5, so produced are similar to those observed from blazars.

We present the X-UDS survey, a set of wide and deep Chandra observations of the Subaru-XMM Deep/UKIDSS Ultra Deep Survey (SXDS/UDS) field. The survey consists of 25 observations that cover a total area of 0.33 deg$^{2}$. The observations are combined to provide a nominal depth of ~600 ksec in the central 100 arcmin$^{2}$ region of the field that has been imaged with Hubble/WFC3 by the CANDELS survey and $\sim$200 ksec in the remainder of the field. In this paper, we outline the survey's scientific goals, describe our observing strategy, and detail our data reduction and point source detection algorithms. Our analysis has resulted in a total of 868 band-merged point sources detected with a false-positive Poisson probability of $<1\times10^{-4}$. In addition, we present the results of an X-ray spectral analysis and provide best-fitting neutral hydrogen column densities, $N_{\rm H}$, as well as a sample of 51 Compton-thick active galactic nucleus candidates. Using this sample, we find the intrinsic Compton-thick fraction to be 30-35% over a wide range in redshift ($z=0.1-3$), suggesting the obscured fraction does not evolve very strongly with epoch. However, if we assume that the Compton-thick fraction is dependent on luminosity, as is seen for Compton-thin sources, then our results are consistent with a rise in the obscured fraction out to $z\sim3$. Finally, an examination of the host morphologies of our Compton-thick candidates shows a high fraction of morphological disturbances, in agreement with our previous results. All data products described in this paper are made available via a public website.

We present the X-UDS survey, a set of wide and deep Chandra observations of the Subaru-XMM Deep/UKIDSS Ultra Deep Survey (SXDS/UDS) field. The survey consists of 25 observations that cover a total area of 0.33 deg$^{2}$. The observations are combined to provide a nominal depth of ~600 ksec in the central 100 arcmin$^{2}$ region of the field that has been imaged with Hubble/WFC3 by the CANDELS survey and $\sim$200 ksec in the remainder of the field. In this paper, we outline the survey's scientific goals, describe our observing strategy, and detail our data reduction and point source detection algorithms. Our analysis has resulted in a total of 868 band-merged point sources detected with a false-positive Poisson probability of $<1\times10^{-4}$. In addition, we present the results of an X-ray spectral analysis and provide best-fitting neutral hydrogen column densities, $N_{\rm H}$, as well as a sample of 51 Compton-thick active galactic nucleus candidates. Using this sample, we find the intrinsic Compton-thick fraction to be 30-35% over a wide range in redshift ($z=0.1-3$), suggesting the obscured fraction does not evolve very strongly with epoch. However, if we assume that the Compton-thick fraction is dependent on luminosity, as is seen for Compton-thin sources, then our results are consistent with a rise in the obscured fraction out to $z\sim3$. Finally, an examination of the host morphologies of our Compton-thick candidates shows a high fraction of morphological disturbances, in agreement with our previous results. All data products described in this paper are made available via a public website.

The Simons Observatory (SO) will measure the cosmic microwave background (CMB) in both temperature and polarization over a wide range of angular scales and frequencies from 27-270 GHz with unprecedented sensitivity. One technology for coupling light onto the $\sim$50 detector wafers that SO will field is spline-profiled feedhorns, which offer tunability between coupling efficiency and control of beam polarization leakage effects. We will present efforts to scale up feedhorn production for SO and their viability for future CMB experiments, including direct-machining metal feedhorn arrays and laser machining stacked Si arrays.

NGC 5548 has been well spectroscopically monitored for reverberation mapping of the central kinematics by 19 campaigns. Using the maximum entropy method in this Letter, we build up a high-quality velocity-delay map of the H$\beta$ emission line in the light curves of the continuum and the line variations observed between 2015-2016. The map shows the response strength and lags of the velocity fields of the H$\beta$ emitting regions. The velocity-delay structure of the map is generally symmetric, with strong red and blue wings at time lag $\tau \leq 15$ days, a narrower velocity distribution at $\tau \geq 15$ days, and a deficit of response in the core. This is suggestive of a disk geometry of the broad-line region (BLR). The relatively weaker H$\beta$ response at the longer lags in the red side indicates anisotropic emission from the outer part of the BLR. We also recover the velocity-delay maps of NGC 5548 from the historical data of 13 years to investigate the long-term variability of its BLR. In general, the BLR of NGC 5548 was switching between the inflow and virialized phases in the past years. The resultant maps of seven years reveal inflow signatures and show decreasing lags, indicating that the changes in the BLR size are related to the infalling BLR gas. The other four maps show potential disk signatures which are similar to our map.

Magnetic clouds (MCs) are transient magnetic structures giving the strongest southward magnetic field (Bz south) in the solar wind. The sheath regions of MCs may also carry southward magnetic field. Southward magnetic field is responsible for causing space-weather disturbances. We report a comprehensive analysis of MCs and Bz components in their sheath regions during 1995 to 2017. Eighty-five percent of 303 MCs contain a south Bz up to 50 nT. Sheath Bz during the 23 years may reach as high as 40 nT. The MCs of strongest magnetic magnitude and Bz south occur in the declining phase of the solar cycle. The bipolar MCs have solar-cycle dependence in their polarity, but not in the occurrence frequency. Unipolar MCs show solar-cycle dependence in their occurrence frequency but not in their polarity. MCs with the highest speeds, largest total B magnitudes and sheath Bz south are from source regions closer to the solar disk center. About 80% of large Dst storms are caused by MC events. The combinations of south Bz in the sheath and the south-first MCs in close succession have given the largest storms. The solar-cycle dependence of bipolar MCs is extended to 2017, spanning 42 years. We find that the bipolar MC Bz polarity solar-cycle dependence is given by MCs originated from quiescent filaments in decayed active regions and a group of weak MCs of unclear sources, while the polarity of bipolar MCs with active-region flares always has mixed Bz polarity without solar-cycle dependence and is therefore the least predictable for Bz forecasting.

The Gravitational wave high-energy Electromagnetic Counterpart All-sky Monitor (GECAM) project is the planned Chinese space telescope for detecting the X and gamma-ray counterpart. It consists of two micro-satellites in low earth orbit with the advantages of instantaneous full-sky coverage, low energy threshold down to 6 keV and can be achieved within a short period and small budget. Due to the limitation of size, weight and power consumption of micro-satellites, silicon photomultipliers (SiPMs) are used to replace the photomultiplier tubes (PMTs) to assemble a novel gamma-ray detector. A prototype of a SiPM array with LaBr3 crystal is built and tested, and it shows a high detection efficiency (70% at 5.9 keV) and an acceptable uniformity. The low-energy X-ray of 5.9 keV can be detected by a simply readout circuit, and the energy resolution is 6.5% (FWHM) at 662 keV. The design and performance of the detector are discussed in detail in this paper.

The Slowly Pulsating B-type (SPB) stars are the upper main-sequence stars on the HR diagram. Their oscillations are high-order, low-degree g-mode and can be used to probe the structure of the radiative zone where is near the outer boundary of the convective core and constrain the chemical mixing in stellar interiors. In SPB stars, the period spacing periodically varies with periods. It has regarded as a signature of the chemical composition gradient beyond the convective core. Based on theoretical calculations, we find that the variation frequency of the period spacings ($f_{\Delta P}$) is related to the width of the $\mu$-gradient region on the buoyancy radius ($\Lambda_{\mu}$) with the relation of $f_{\Delta P}\sim0.5\Lambda_{\mu}$. This indicates that the variation frequency $f_{\Delta P}$ is sensitive to the central hydrogen mass fraction $X_{\rm C}$ (i.e., the evolutionary status). Finally, we find that the variation frequency $f_{\Delta P}$ and the means of the period spacings $\langle\Delta P\rangle$ can be used to construct a new C-D-like diagram ($f_{\Delta P}$ vs. $\langle\Delta P\rangle$) which can be used to roughly decide the stellar evolutionary stages and to approximately determine stellar mass for SPB stars.