The mean stellar alpha-to-iron abundance ratio ([$\alpha$/Fe]) of a galaxy is an indicator of galactic star formation timescale. It is important for understanding the star formation history of early-type galaxies (ETGs) as their star formation processes have basically stopped. Using the model templates which are made by Vazdekis et al., we apply the pPXF based spectral fitting method to estimate the [$\alpha$/Fe] of 196 high signal-to-noise ratio ETGs from the MaNGA survey. The velocity dispersions within 1R$_e$ ($\sigma_{e}$) range from 27 to 270 km/s. We find a flat relation between the mean [$\alpha$/Fe] within the 1R$_e^{maj}$ ellipses and log($\sigma_{e}$), even if limiting to the massive sample with log($\sigma_{e}$/km s$^{-1}$)$>$1.9. However, the relation becomes positive after we exclude the Mg$_1$ feature in our fits, which agrees with the results from the previous work with other stellar population models, albeit with relatively large scatter. It indicates that the Vazdekis models have self-inconsistency and do not actually cover the observations in Mg$_1$. For reference, we also measure the stellar population radial gradients within 1R$_e^{maj}$ ellipses. Due to the low resolution of age estimations for old objects and the Mg$_1$ issue, the uncertainties of these gradients cannot be neglected.

A new look at the structure of the electron diffusion region in collisionless magnetic reconnection is presented. The research is based on a particle-in-cell simulation of asymmetric magnetic reconnection, which include a temperature gradient across the current layer in addition to density and magnetic field gradient. We find that none of X-point, flow stagnation point, and local current density peak coincide. Current and energy balance analyses around the flow stagnation point and current density peak show consistently that current dissipation is associated with the divergence of nongyrotropic electron pressure. Furthermore, the same pressure terms, when combined with shear-type gradients of the electron flow velocity, also serve to maintain local thermal energy against convective losses. These effects are similar to those found also in symmetric magnetic reconnection. In addition, we find here significant effects related to the convection of current, which we can relate to a generalized diamagnetic drift by the nongyrotropic pressure divergence. Therefore, only part of the pressure force serves to dissipate the current density. However, the prior conclusion that the role of the reconnection electric field is to maintain the current density, which was obtained for a symmetric system, applies here as well. Finally, we discuss related features of electron distribution function in the EDR.

Energetic electrons accelerated by solar flares often give rise to type III radio bursts at a broad waveband and even interplanetary type III bursts (IT3) if the wavelength extends to decameter-kilometer. In this Letter, we investigate the probability of the flares that produce IT3, based on the sample of 2272 flares above M-class observed from 1996 to 2016. It is found that only 49.6% of the flares are detected to be accompanied with IT3. The duration, peak flux, and fluence of the flares with and without IT3 both present power-law distributions in the frequency domain, but the corresponding spectral indices for the former (2.06$\pm$0.17, 2.04$\pm$0.18, and 1.55$\pm$0.09) are obviously smaller than that for the latter (2.82$\pm$0.22, 2.51$\pm$0.19, and 2.40$\pm$0.09), showing that the flares with IT3 have longer durations and higher peak fluxes. We further examine the relevance of coronal mass ejections (CMEs) to the two groups of flares. It is found that 58% (655 of 1127) of the flares with IT3 but only 19% (200 of 1078) of the flares without IT3 are associated with CMEs, and that the associated CMEs for the flares with IT3 are inclined to be wider and faster. This indicates that CMEs may also play a role in producing IT3, speculatively facilitating the escape of accelerated electrons from the low corona to the interplanetary space.

Active galactic nuclei (AGNs) have been attracting research attention due to their special observable properties. Specifically, a majority of AGNs are detected by Fermi-LAT missions, but not by Fermi-LAT, which raises the question of whether any differences exist between the two. To answer this issue, we compile a sample of 291 superluminal AGNs (189 FDSs and 102 non-FDSs) from available multi-wavelength radio, optical, and X-ray (or even $\gamma$-ray) data and Doppler factors and proper motion ($\mu$) (or apparent velocity ($\beta_{\rm{app}}$)); calculated the apparent velocity from their proper motion, Lorentz factor ($\Gamma$), viewing angle ($\phi$) and co-moving viewing angle ($\phi_{co}$) for the sources with available Doppler factor ($\delta$); and performed some statistical analyses for both types. Our study indicated that1. In terms of average values, FDSs have higher proper motions ($\mu$), apparent velocities ($\beta_{\rm app}$), Doppler factor ($\delta$), Lorentz factor ($\Gamma$), and smaller viewing angle ($\phi$). Nevertheless, there is no clear difference in co-moving viewing angles ($\phi_{\rm co}$).

Convolutional Neutral Networks have been successfully applied in searching for strong lensing systems, leading to discoveries of new candidates from large surveys. On the other hand, systematic investigations about their robustness are still lacking. In this paper, we first construct a neutral network, and apply it to $r$-band images of Luminous Red Galaxies (LRGs) of the Kilo Degree Survey (KiDS) Data Release 3 to search for strong lensing systems. We build two sets of training samples, one fully from simulations, and the other one using the LRG stamps from KiDS observations as the foreground lens images. With the former training sample, we find 48 high probability candidates after human-inspection, and among them, 27 are newly identified. Using the latter training set, about 67\% of the aforementioned 48 candidates are also found, and there are 11 more new strong lensing candidates identified. We then carry out tests on the robustness of the network performance with respect to the variation of PSF. With the testing samples constructed using PSF in the range of 0.4 to 2 times of the median PSF of the training sample, we find that our network performs rather stable, and the degradation is small. We also investigate how the volume of the training set can affect our network performance by varying it from 0.1 millions to 0.8 millions. The output results are rather stable showing that within the considered range, our network performance is not very sensitive to the volume size.

Magnetic reconnection in strongly magnetized astrophysical plasma environments is believed to be the primary process for fast energy release and particle energization. Currently there is strong interest in relativistic magnetic reconnection, in that it may provide a new explanation for high-energy particle acceleration and radiation in strongly magnetized astrophysical systems. We review recent advances in particle acceleration and reconnection physics in the magnetically-dominated regime. More discussion is focused on the physics of particle acceleration, power-law formation as well as the reconnection rate problem. In addition, we provide an outlook for studying reconnection acceleration mechanisms and kinetic physics in the next step.

The Coronal Global Evolutionary Model (CGEM) provides data-driven simulations of the magnetic field in the solar corona to better understand the build-up of magnetic energy that leads to eruptive events. The CGEM project has developed six capabilities. CGEM modules (1) prepare time series of full-disk vector magnetic field observations to (2) derive the changing electric field in the solar photosphere over active-region scales. This local electric field is (3) incorporated into a surface flux transport model that reconstructs a global electric field that evolves magnetic flux in a consistent way. These electric fields drive a (4) 3D spherical magneto-frictional (SMF) model, either at high-resolution over a restricted range of solid angle or at lower resolution over a global domain, to determine the magnetic field and current density in the low corona. An SMF-generated initial field above an active region and the evolving electric field at the photosphere are used to drive (5) detailed magneto-hydrodynamic (MHD) simulations of active regions in the low corona. SMF or MHD solutions are then used to compute emissivity proxies that can be compared with coronal observations. Finally, a lower-resolution SMF magnetic field is used to initialize (6) a global MHD model that is driven by an SMF electric-field time series to simulate the outer corona and heliosphere, ultimately connecting Sun to Earth. As a demonstration, this report features results of CGEM applied to observations of the evolution of NOAA Active Region 11158 in February 2011.

In this paper, we present a cosmological model-independent test for the cosmic opacity at high redshifts ($z\sim5$). We achieve this with the opacity-dependent luminosity distances derived from nonlinear relation between X-ray and UV emissions of quasars, combined with two types of opacity-independent luminosity distances derived from the Hubble parameter measurements and simulated gravitational wave (GW) events achievable with the Einstein Telescope (ET). In the framework of two phenomenological parameterizations adopted to describe cosmic opacity at high redshifts, our main results show that a transparent universe is supported by the current observational data at 2$\sigma$ confidence level. However, the derived value of the cosmic opacity is slightly sensitive to the parametrization of $\tau(z)$, which highlights the importance of choosing a reliable parametrization to describe the optical depth $\tau(z)$ in the early universe. Compared with the previous works, the combination of the quasar data and the $H(z)$/GW observations in similar redshift ranges provides a novel way to confirm a transparent universe ($\epsilon=0$ at higher redshifts $z\sim 5$), with an accuracy of $\Delta \epsilon\sim 10^{-2}$. More importantly, our findings indicate that a strong degeneracy between the cosmic opacity parameter and the parameters characterizing the $L_{UV}-L_X$ relation of quasars, which reinforces the necessity of proper calibration for such new type of high-redshift standard candle (in a cosmological model-independent way).

The authors of this white paper met on 16-17 January 2020 at the New Jersey Institute of Technology, Newark, NJ, for a 2-day workshop that brought together a group of heliophysicists, data providers, expert modelers, and computer/data scientists. Their objective was to discuss critical developments and prospects of the application of machine and/or deep learning techniques for data analysis, modeling and forecasting in Heliophysics, and to shape a strategy for further developments in the field. The workshop combined a set of plenary sessions featuring invited introductory talks interleaved with a set of open discussion sessions. The outcome of the discussion is encapsulated in this white paper that also features a top-level list of recommendations agreed by participants.

As an important source for large geomagnetic storms, an "ICME-in-sheath" is a completely shocked interplanetary coronal mass ejection (ICME) stuck in the sheath between a shock and host ejecta. Typical characteristics are identified from coordinated multi-sets of observations: (1) it is usually short in duration and lasts a few hours at 1 AU; (2) its solar wind parameters, in particular the magnetic field, seem to keep enhanced for a large range of distances; and (3) common ICME signatures are often lost. The host ejecta could be a single ICME or a complex ejecta, being fast enough to drive a shock. These results clarify previous misinterpretations of this phenomenon as a normal part of a sheath region. The "ICME-in-sheath" phenomenon, together with a preconditioning effect, produced an extreme set of the magnetic field, speed and density near 1 AU in the 2012 July 23 case, all around their upper limits at the same time. This is probably the most extreme solar wind driving at 1 AU and enables us to estimate the plausible upper limit for geomagnetic storm activity. With an appropriate modification in the southward field, we suggest that a geomagnetic storm with a minimum $D_{\rm st}$ of about $-2000$ nT could occur in principle. The magnetopause would be compressed to about 3.3 Earth radii from the Earth's center, well inside the geosynchronous orbit.

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.

We present new Atacama Large Millimeter/submillimeter Array (ALMA) observations for three protoplanetary disks in Taurus at 2.9\,mm and comparisons with previous 1.3\,mm data both at an angular resolution of $\sim0.''1$ (15\,au for the distance of Taurus). In the single-ring disk DS Tau, double-ring disk GO Tau, and multiple-ring disk DL Tau, the same rings are detected at both wavelengths, with radial locations spanning from 50 to 120\,au. To quantify the dust emission morphology, the observed visibilities are modeled with a parametric prescription for the radial intensity profile. The disk outer radii, taken as 95\% of the total flux encircled in the model intensity profiles, are consistent at both wavelengths for the three disks. Dust evolution models show that dust trapping in local pressure maxima in the outer disk could explain the observed patterns. Dust rings are mostly unresolved. The marginally resolved ring in DS Tau shows a tentatively narrower ring at the longer wavelength, an observational feature expected from efficient dust trapping. The spectral index ($\alpha_{\rm mm}$) increases outward and exhibits local minima that correspond to the peaks of dust rings, indicative of the changes in grain properties across the disks. The low optical depths ($\tau\sim$0.1--0.2 at 2.9\,mm and 0.2--0.4 at 1.3\,mm) in the dust rings suggest that grains in the rings may have grown to millimeter sizes. The ubiquitous dust rings in protoplanetary disks modify the overall dynamics and evolution of dust grains, likely paving the way towards the new generation of planet formation.

The recent detections of the gravitational wave (GW) event GW170817 and its electromagnetic counterpart GRB 170817A produced by a binary neutron star (NS) merger is a new milestone of multimessenger astronomy. The time interval between these two signals has attracted widespread attention from physicists. In the braneworld scenario, GWs could propagate through the bulk while electromagnetic waves (EMWs) are bounded on the brane, i.e., our Universe. Therefore, the trajectories of GWs and EMWs may follow different pathes. If GWs and EMWs are originated simultaneously from the same source on the brane, they are expected to arrive at the observer successively. Consequently, the time delay between GW170817 and GRB 170817A may carry the information of the extra dimension. In this paper, we try to investigate the phenomenon in the context of a five-dimensional dS ($\text{dS}_5$) spacetime. We first study two special Universe models, i.e., de Sitter and Einstein-de Sitter models, and calculate the gravitation horizon radius in each case. For the real Universe, we then consider the $\Lambda$CDM model. Our results show that for the de Sitter model of the Universe, the $\text{dS}_5$ radius could not contribute to the time delay. With the data of the observation, we constrain the $\text{dS}_5$ radius to $\ell\gtrsim7.5\times10^{2}\,\text{Tpc}$ for the Einstein-de Sitter model and $\ell\gtrsim2.4\times10^{3}\,\text{Tpc}$ for the $\Lambda$CDM model. After considering the uncertainty in the source redshift and the time-lags given by different astrophysical processes of the binary NS merger, we find that our constraints are not sensitive to the redshift in the range of (0.005, 0.01) and the time-lag in the range of (-100s, 1.734s).

We apply the independent component analysis (ICA) to the real data from a gravitational wave detector for the first time. Specifically we use the iKAGRA data taken in April 2016, and calculate the correlations between the gravitational wave strain channel and 35 physical environmental channels. Using a couple of seismic channels which are found to be strongly correlated with the strain, we perform ICA. Injecting a sinusoidal continuous signal in the strain channel, we find that ICA recovers correct parameters with enhanced signal-to-noise ratio, which demonstrates usefulness of this method. Among the two implementations of ICA used here, we find the correlation method yields the optimal result for the case environmental noises act on the strain channel linearly.

The first charge-coupled device photometric light curves in B, V, R${_c}$, and I${_c}$ bands of the short-period contact binary star OQ Cam are presented. Through analyzing the light curves with the Wilson-Devinney code, it is found that OQ Cam is a W-type shallow contact binary star with a contact degree of f =20.2%. Based on the O-C analyzing, the orbit period displays a long term increasing with a rate of $dP/dt=4.40\times10{^{-7}} day \cdot year{^{-1}}$. The increasing in orbit period can be interpreted by mass transfer from the less massive component to the more massive one. As the orbit period increasing, OQ Cam would evolve from the present shallow contact configuration to a none contact stage. So it may be a potential candidate to confirm the thermal relaxation oscillation theory.

The cosmic curvature, a fundamental parameter for cosmology could hold deep clues to inflation and cosmic origins. We propose an improved model-independent method to constrain the cosmic curvature by combining the constructed Hubble diagram of high-redshift quasars with galactic-scale strong lensing systems expected to be seen by the forthcoming LSST survey. More specifically, the most recent quasar data are used as a new type of standard candles in the range $0.036<z<5.100$, whose luminosity distances can be directly derived from the non-linear relation between X-ray and UV luminosities. Compared with other methods, the proposed one involving the quasar data achieves constraints with higher precision ($\Delta \Omega_k\sim 10^{-2}$) at high redshifts ($z\sim 5.0$). We also investigate the influence of lens mass distribution in the framework of three types of lens models extensively used in strong lensing studies (SIS model, power-law spherical model, and extended power-law lens model), finding the strong correlation between the cosmic curvature and the lens model parameters. When the power-law mass density profile is assumed, the most stringent constraint on the cosmic curvature $\Omega_k$ can be obtained. Therefore, the issue of mass density profile in the early-type galaxies is still a critical one that needs to be investigated further.

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.

Context: Instrumental profile (IP) is the basic property of a spectrograph. Accurate IP characterisation is the prerequisite of accurate wavelength solution. It also facilitates new spectral acquisition methods such as the forward modeling and deconvolution. Aims: We investigate an IP modeling method for the fibre-fed echelle spectrograph with the emission lines of the ThAr lamp, and explore the method to evaluate the accuracy of IP characterisation. Methods: The backbone-residual (BR) model is put forward and tested on the fibre-fed High Resolution Spectrograph (HRS) at the Chinese Xinglong 2.16-m Telescope, which is the sum of the backbone function and the residual function. The backbone function is a bell-shaped function to describe the main component and the spatial variation of IP. The residual function, which is expressed as the cubic spline function, accounts for the difference between the bell-shaped function and the actual IP. The method of evaluating the accuracy of IP characterisation is based on the spectral reconstruction and Monte Carlo simulation. Results: The IP of HRS is characterised with the BR model, and the accuracy of the characterised IP reaches 0.006 of the peak value of the backbone function. This result demonstrates that the accurate IP characterisation has been achieved on HRS with the BR model, and the BR model is an excellent choice for accurate IP characterisation of fibre-fed echelle spectrographs.

The primary task of the 1.26-m telescope jointly operated by the National Astronomical Observatory and Guangzhou University is photometric observations of the g, r, and i bands. A data processing pipeline system was set up with mature software packages, such as IRAF, SExtractor, and SCAMP, to process approximately 5 GB of observational data automatically every day. However, the success ratio was significantly reduced when processing blurred images owing to telescope tracking error; this, in turn, significantly constrained the output of the telescope. We propose a robust automated photometric pipeline (RAPP) software that can correctly process blurred images. Two key techniques are presented in detail: blurred star enhancement and robust image matching. A series of tests proved that RAPP not only achieves a photometric success ratio and precision comparable to those of IRAF but also significantly reduces the data processing load and improves the efficiency.

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

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

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

This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process.

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.

The column density probability density functions (PDFs) of the mass density is a tool to study turbulence and self-gravity in molecular clouds. Self-absorption significantly makes the intensity PDFs different from the column density ones. In this work, we study the effects of self-absorption on the intensity PDFs utilizing three synthetic emission lines of CO isotopologs 12CO(1-0), 13CO(1-0), and C18O(1-0). We find that for supersonic turbulence, the intensity PDFs change their shape from the log-normal to a distribution that has a power-law part depending on the optical depth. We compute the intensity PDFs for massive star-forming region Vela C and find the change of intensity PDFs in observation agrees with our numerical results. We identify the gravitational collapsing regions on Ophiuchus molecular cloud through the PDFs of H2 column density. We see that the change of CO isotopologs' intensity PDFs is independent of self-gravity, which makes the intensity PDFs less reliable in identifying gravitational collapsing regions using emission lines.

A comprehensive analysis of plasma and composition characteristics inside magnetic clouds (MCs) observed by the Advanced Composition Explorer (ACE) spacecraft from 1998 February to 2011 August is presented. The results show that MCs have specific interior structures, and MCs of different speeds show differences in composition and structure. Compared with the slow MCs, fast MCs have enhanced mean charge states of iron, oxygen, silicon, magnesium, $\mathrm{O^{7+}/O^{6+}}$, $\mathrm{C^{6+}/C^{5+}}$, $\mathrm{C^{6+}/C^{4+}}$ and $\mathrm{Fe^{\geq16+}/Fe_{total}}$ values. For ionic species in fast MCs, a higher atomic number represents a greater enhancement of mean charge state than slow MCs. We also find that both the fast and slow MCs display bimodal structure distribution in the mean iron charge state ($\mathrm{\langle Q\rangle Fe}$), which suggests that the existence of flux rope prior to the eruption is common. Furthermore, the $\mathrm{\langle Q\rangle Fe} $, $\mathrm{Fe^{\geq16+}/Fe_{total}}$, and $\mathrm{O^{7+}/O^{6+}}$ ratio distribution inside fast MCs have the feature that the posterior peak is higher than the anterior one. This result agrees with the "standard model" for CME/flares, by which magnetic reconnection occurs beneath the flux rope, thereby ionizing the ions of the posterior part of flux rope sufficiently by high-energy electron collisions or by direct heating in the reconnection region.

In this work, we investigate what role the redshift drift data of Square Kilometre Array (SKA) will play in the cosmological parameter estimation in the future. To test the constraint capability of the redshift drift data of SKA-only, the $\Lambda$CDM model is chosen as a reference model. We find that using the SKA1 mock data, the $\Lambda$CDM model can be loosely constrained, while the model can be well constrained when the SKA2 mock data are used. When the mock data of SKA are combined with the data of the European Extremely Large Telescope (E-ELT), the constraints can be significantly improved almost as good as the data combination of the type Ia supernovae observation (SN), the cosmic microwave background observation (CMB), and the baryon acoustic oscillations observation (BAO). Furthermore, we explore the impact of the redshift drift data of SKA on the basis of SN+CMB+BAO+E-ELT in the $\Lambda$CDM model, the $w$CDM model, the CPL model, and the HDE model. We find that the redshift drift measurement of SKA could help to significantly improve the constraints on dark energy and could break the degeneracy between the cosmological parameters. Therefore, we conclude that redshift-drift observation of SKA would provide a good improvement in the cosmological parameter estimation in the future and has the enormous potential to be one of the most competitive cosmological probes in constraining dark energy.

In this work, we propose a powerful probe of neutrino effects on the large-scale structure (LSS) of the Universe, i.e., Minkowski functionals (MFs). The morphology of LSS can be fully described by four MFs. This tool, with strong statistical power, is robust to various systematics and can comprehensively probe all orders of N-point statistics. By using a pair of high-resolution N-body simulations, for the first time, we comprehensively studied the subtle neutrino effects on the morphology of LSS. For an ideal LSS survey of volume $\sim1.73$ Gpc$^3$/$h^3$, neutrino signals are mainly detected from void regions with a significant level up to $\thicksim 10\sigma$ and $\thicksim 300\sigma$ for CDM and total matter density fields, respectively. This demonstrates its enormous potential for much improving the neutrino mass constraint in the data analysis of up-coming ambitious LSS surveys.

The Galactic plane scanning survey is one of the main scientific objectives of the Hard X-ray Modulation Telescope (known as Insight-HXMT). During the two-year operation of Insight-HXMT, more than 1000 scanning observations have been performed and the whole Galactic plane ($\rm 0^{\circ}<l<360^{\circ}$, $\rm -10^{\circ}<b<10^{\circ}$) has been covered completely. We summarize the Galactic plane scanning survey of Insight-HXMT for two years, including the characteristics of the scanning data, the data analysis process and the preliminary results of the Low-Energy telescope, the Medium-Energy telescope and the High-Energy telescope. With the light curve PSF fitting method, the fluxes of the known sources in the scanned area as well as the flux errors are obtained for each scanning observation. From the relationships of SNRs and fluxes, the $5\sigma$ sensitivities of three telescopes of Insight-HXMT are estimated as $\rm \sim7.6\times10^{-11}~erg cm^{-2}~s^{-1}$ ($\rm 3 mCrab,~1-6 keV$), $\rm \sim4.0\times10^{-10}~erg~cm^{-2}~s^{-1}$ ($\rm 20~mCrab,~7-40 keV$) and $\rm \sim2.6\times10^{-10}~erg cm^{-2}~s^{-1}$ ($\rm 18 mCrab,~25-100 keV$) for an individual scanning observation of $2-3$ hours, respectively. Up to September 2019, more than 800 X-ray sources with various types are monitored by the three telescopes and their long-term light curves with three energy bands are obtained to make further scientific analyses.

In the third catalog of active galactic nuclei detected by the $Fermi$ Large Area Telescope Clean (3LAC) sample, there are 402 blazars candidates of uncertain type (BCU). The proposed analysis will help to evaluate the potential optical classification flat spectrum radio quasars (FSRQs) versus BL Lacertae (BL Lacs) objects of BCUs, which can help to understand which is the most elusive class of blazar hidden in the Fermi sample. By studying the 3LAC sample, we found some critical values of $\gamma$-ray photon spectral index ($\Gamma_{\rm ph}$), variability index (VI) and radio flux (${\rm F_R}$) of the sources separate known FSRQs and BL Lac objects. We further utilize those values to defined an empirical "high-confidence" candidate zone that can be used to classify the BCUs. Within such a zone ($\Gamma_{\rm ph}<2.187$, log${\rm F_R}<2.258$ and ${ \rm logVI <1.702}$), we found that 120 BCUs can be classified BL Lac candidates with a higher degree of confidence (with a misjudged rate $<1\%$). Our results suggest that an empirical "high confidence" diagnosis is possible to distinguish the BL Lacs from the Fermi observations based on only on the direct observational data of $\Gamma_{\rm ph}$, VI and ${\rm F_R}$.

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

The far-infrared (FIR) regime is one of the few wavelength ranges where no astronomical data with sub-arcsecond spatial resolution exist. Neither of the medium-term satellite projects like SPICA, Millimetron nor O.S.T. will resolve this malady. For many research areas, however, information at high spatial and spectral resolution in the FIR, taken from atomic fine-structure lines, from highly excited carbon monoxide (CO), light hydrids, and especially from water lines would open the door for transformative science. A main theme will be to trace the role of water in proto-planetary disks, to observationally advance our understanding of the planet formation process and, intimately related to that, the pathways to habitable planets and the emergence of life. Furthermore, key observations will zoom into the physics and chemistry of the star-formation process in our own Galaxy, as well as in external galaxies. The FIR provides unique tools to investigate in particular the energetics of heating, cooling and shocks. The velocity-resolved data in these tracers will reveal the detailed dynamics engrained in these processes in a spatially resolved fashion, and will deliver the perfect synergy with ground-based molecular line data for the colder dense gas.

We show, through an explicit calculation of the relevant Green's functions, that the transverse-traceless (TT) portion of the gravitational perturbations of Minkowski spacetime and of spatially flat cosmologies with a constant equation-of-state $w$ receive contributions from their isolated matter source(s) outside the past null cone of the observer. This implies the TT gravitational wave (GW) cannot be a standalone observable -- despite widespread (apparent) claims in the gravitational wave literature to the contrary. About a Minkowski background, all 4 of the gauge-invariant variables -- the two scalars, one vector and tensor -- play crucial roles to ensure the spatial tidal forces encoded within the gauge-invariant linearized Riemann tensor are causal. These gravitational tidal forces do not depend solely on the TT graviton but rather on the causal portion of its acceleration. However, in the far zone radiative limit, the flat spacetime `TT' graviton Green's function does reduce to the causal `tt' ones, which are the ones commonly used to compute gravitational waveforms. Similar remarks apply to the spin-1 photon; for instance, the electric field does not depend solely on the photon, but is the causal part of its velocity. As is known within the quantum theory of photons and linearized gravitons, there are obstacles to the construction of simultaneously gauge-invariant and Lorentz-covariant descriptions of these massless spin-1 and spin-2 states. Our results transparently demonstrate that the quantum operators associated with the helicity-1 photon and helicity-2 linear graviton both violate micro-causality: namely, they do not commute outside the light cone in flat and cosmological spacetimes.

We make the in-orbit calibration to the point-spread functions (PSFs) of the collimators of the Hard X-ray Modulation Telescope with the scanning observation of the Crab. We construct the empirical adjustments to the theoretically calculated geometrical PSFs. The adjustments contain two parts: a rotating matrix to adjust the directional deviation of the collimators and a paraboloidal function to correct the inhomogeneity of the real PSFs. The parameters of the adjusting matrices and paraboloidal functions are determined by fitting the scanning data with lower scanning speed and smaller intervals during the calibration observations. After the PSF calibration, the systematic errors in source localization in the Galactic plane scanning survey are 0.010 deg, 0.015 deg, 0.113 deg for the Low-Energy Telescope (LE), the Medium-Energy telescope (ME) and the High-Energy telescope (HE), respectively; meanwhile, the systematic errors in source flux estimation are 1.8%, 1.6%, 2.7% for LE, ME and HE, respectively.

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

Multi-peaked features appear on the integrated pulse profiles of fast radio burst observed below 2.5 GHz and the instantaneous spectrum of many bursts observed between 4 and 8 GHz. The mechanism of pulse or spectrum shaping has attracted little attention. Here we show that these interference-like pulse profiles are mostly the instantaneous spectra near the source regions of fast radio bursts. The corresponding instantaneous spectra are coincident to the spectrum from a single electron passing through a tapered undulator. The multi-peaked spectrum observed between 4 and 8 GHz can also be explained consistently by this type of spectrum. The spectrum is invisible unless the particles in the radiation beam are bunched. The bunching effect is probably due to the acceleration of particles in the plasma wakefield.

Astrophysical observations provide a unique opportunity to test possible signatures of Lorentz Invariance Violation (LIV), due to the high energies and long distances involved. In quantum theory of gravity, one may expect the modification of the dispersion relation between energy and momentum for photons, which can be probed with the time-lag (the arrival time delay between light curves in different energy bands) of Gamma-ray bursts (GRBs). In this paper, by using the detailed time-delay measurements of GRB 160625B at different energy bands, as well as 23 time-delay GRBs covering the redshifts range of $z=0.168-2.5$ (which were measured at different energy channels from the light curves), we propose an improved model-independent method (based on the newly-compiled sample of $H(z)$ measurements) to probe the energy-dependent velocity due to the modified dispersion relation for photons. In the framework of a more complex and reasonable theoretical expression to describe the time delays, our results imply that the intrinsic time lags can be better described with more GRBs time delay data. More importantly, through direct fitting of the time-delay measurements of a sample of GRBs, our limit on the LIV energy scale is comparable to that with unknown constant for the intrinsic time lag, much lower than the Planck energy scale in both linear LIV and quadratic LIV cases.

A deep learning network, Long-Short Term Memory (LSTM) network, is used in this work to predict whether the maximum flare class an active region (AR) will produce in the next 24 hours is class $\Gamma$. We considered $\Gamma$ are $\ge M$, $\ge C$ and any flare class. The essence of using LSTM, which is a recurrent neural network, is its capability to capture temporal information of the data samples. The input features are time sequences of 20 magnetic parameters from SHARPs - Space-weather HMI Active Region Patches. We analyzed active regions from June 2010 to Dec 2018, using the Geostationary Operational Environmental Satellite (GOES) X-ray flare catalogs and label the data samples with identified ARs in the GOES X-ray flare catalogs. Our results (i) shows consistent skill scores with recently published results using LSTMs and better than the previous work using single time input (eg. DeFN) (ii) The skill scores from the model show essential differences when different years of data was chosen for training and testing.

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

We provide evidence that gravitational radiation in a 4D radiation-dominated universe, with equation-of-state $w=1/3$, consists of two components: helicity-2 gravitons and massless scalar acoustic waves. On physical grounds, we would expect the homogeneous solution of the Weyl tensor components to be a good approximation to its inhomogeneous counterparts, whenever the observer is located well in the far zone of an isolated astrophysical source of cosmological gravitational radiation. We show explicitly that these homogeneous and inhomogeneous solutions both receive contributions from the gauge-invariant tensor {\it and} the two Bardeen (acoustic) scalars. Comparison of these Weyl tensor computations thus allows us to not only identify, in the high frequency limit, the corresponding gravitational tensor and scalar radiation; but also their oscillatory polarization patterns.

It is known that a five-dimensional free vector field $A_{M}$ cannot be localized on Randall-Sundrum (RS)-like thick branes, namely, the thick branes embedded in asymptotic Anti-de Sitter (AdS) spacetime. In order to localize a vector field on the RS-like thick brane, an extra coupling term should be introduced. In this paper, we generalize the geometrical coupling mechanism by adding two mass terms ($\alpha Rg^{MN}A_{M}A_{N}+\beta R^{MN}A_{M}A_{N}$) into the action. We decompose the fundamental vector field $A_{M}$ into three parts: transverse vector part $\hat{A}_{\mu}$, scalar parts $\phi$ and $A_{5}$. Then, we find that the transverse vector part $\hat{A}_{\mu}$ decouples from the scalar parts. In order to eliminate the tachyonic modes of $\hat{A}_{\mu}$, the two coupling parameters $\alpha$ and $\beta$ should satisfy a relation. Combining the restricted condition, we can get a combination parameter as $\gamma=\frac{3}{2}\pm\sqrt{1+12\alpha}$. Only if $\gamma>1/2$, the zero mode of $\hat{A}_{\mu}$ can be localized on the RS-like thick brane. We also investigate the resonant character of the vector part $\hat{A}_{\mu}$ for the general RS-like thick brane with the warp factor $A(z)=-\ln(1+k^2z^2)/2$ by choosing the relative probability method. The result shows that, only for $\gamma>3$, the massive resonant Kaluza-Klein modes can exist. The number of resonant Kaluza-Klein states increases with the combination parameter $\gamma$, and the lifetime of the first resonant state can be long enough as the age of our universe. This indicates that the vector resonances might be considered as one of the candidates of dark matter.

We report a detailed observational study of two quasi-periodic fast-propagating (QFP) magnetosonic wave events occurred on 2011 March 09 and 10, respectively. Interestingly, both the two events have two wave trains (WTs): one main and strong (WT-1) whereas the second appears small and weak (WT-2). Peculiar and common characteristics of the two events are observed, namely: 1) the two QFP waves are accompanied with brightenings during the whole stage of the eruptions; 2) both the two main wave trains are nearly propagating along the same direction; 3) EUV waves are found to be associated with the two events. Investigating various aspects of the target events, we argue that: 1) the second event is accompanied with a flux rope eruption during the whole stage; 2) the second event eruption produces a new filament-like (FL) dark feature; 3) the ripples of the two WT-2 QFP waves seem to result from different triggering mechanisms. Based on the obtained observational results, we propose that the funnel-like coronal loop system is indeed playing an important role in the two WT-1 QFP waves. The development of the second WT-2 QFP wave can be explained as due to the dispersion of the main EUV front. The co-existence of the two events offer thereby a significant opportunity to reveal what driving mechanisms and structures are tightly related to the waves.

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.

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

Cosmological applications of HII galaxies (HIIGx) and giant extragalactic HII regions (GEHR) to construct the Hubble diagram at higher redshifts require knowledge of the "$L$--$\sigma$" relation of the standard candles used. In this paper, we study the properties of a large sample of 156 sources (25 high-$z$ HII galaxies, 107 local HII galaxies, and 24 giant extragalactic HII regions) compiled by Terlevich et al.(2015). Using the the cosmological distances reconstructed through two new cosmology-independent methods, we investigate the correlation between the H$\beta$ emission-line luminosity $L$ and ionized-gas velocity dispersion $\sigma$. The method is based on non-parametric reconstruction using the measurements of Hubble parameters from cosmic clocks, as well as the simulated data of gravitational waves from the third-generation gravitational wave detector (the Einstein Telescope, ET), which can be considered as standard sirens. Assuming the emission-line luminosity versus ionized gas velocity dispersion relation, $\log L ($H$\beta) = \alpha \log \sigma($H$\beta)+\kappa$, we find the full sample provides a tight constraint on the correlation parameters. However, similar analysis done on three different sub-samples seems to support the scheme of treating HII galaxies and giant extragalactic HII regions with distinct strategies. Using the corrected "$L$--$\sigma$" relation for the HII observational sample beyond the current reach of Type Ia supernovae, we obtain a value of the matter density parameter, $\Omega_{m}=0.314\pm0.054$ (calibrated with standard clocks) and $\Omega_{m}=0.311\pm0.049$ (calibrated with standard sirens), in the spatially flat $\Lambda$CDM cosmology.

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

In this paper, we present the application of a new method measuring Hubble parameter $H(z)$ by using the anisotropy of luminosity distance($d_{L}$) of the gravitational wave(GW) standard sirens of neutron star(NS) binary system. The method has never been put into practice so far due to the lack of the ability of detecting GW. However, LIGO's success in detecting GW of black hole(BH) binary system merger announced the potential possibility of this new method. We apply this method to several GW detecting projects, including Advanced LIGO(aLIGO), Einstein Telescope(ET) and DECIGO, and evaluate its constraint ability on cosmological parameters of $H(z)$. It turns out that the $H(z)$ by aLIGO and ET is of bad accuracy, while the $H(z)$ by DECIGO shows a good one. We simulate $H(z)$ data at every 0.1 redshift span using the error information of $H(z)$ by DECIGO, and put the mock data into the forecasting of cosmological parameters. Compared with the previous data and method, we get an obviously tighter constraint on cosmological parameters by mock data, and a concomitantly higher value of Figure of Merit(FoM, the reciprocal of the area enclosed by the $2\sigma$ confidence region). For a 3-year-observation by standard sirens of DECIGO, the FoM value is as high as 170.82. If a 10-year-observation is launched, the FoM could reach 569.42. For comparison, the FoM of 38 actual observed $H(z)$ data(OHD) is 9.3. We also investigate the undulant universe, which shows a comparable improvement on the constraint of cosmological parameters. These improvement indicates that the new method has great potential in further cosmological constraints.

In this paper, we present a scheme to investigate the opacity of the Universe in a cosmological-model-independent way, with the combination of current and future measurements of type Ia supernova sample and galactic-scale strong gravitational lensing systems with SNe Ia acting as background sources. The observational data include the current newly-compiled SNe Ia data (Pantheon sample) and simulated sample of SNe Ia observed by the forthcoming Large Synoptic Survey Telescope (LSST) survey, which are taken for luminosity distances ($D_L$) possibly affected by the cosmic opacity, as well as strongly lensed SNe Ia observed by the LSST, which are responsible for providing the observed time-delay distance ($D_{\Delta t}$) unaffected by the cosmic opacity. Two parameterizations, $\tau(z)=2\beta z$ and $\tau(z)=(1+z)^{2\beta}-1$ are adopted for the optical depth associated to the cosmic absorption. Focusing on only one specific type of standard cosmological probe, this provides an original method to measure cosmic opacity at high precision. Working on the simulated sample of strongly lensed SNe Ia observed by the LSST in 10 year $z$-band search, our results show that, with the combination of the current newly-compiled SNe Ia data (Pantheon sample), there is no significant deviation from the transparency of the Universe at the current observational data level. Moreover, strongly lensed SNe Ia in a 10 year LSST $z$-band search would produce more robust constraints on the validity of cosmic transparency (at the precision of $\Delta\beta=10^{-2}$), with a larger sample of unlensed SNe Ia detected in future LSST survey. We have also discussed the ways in which our methodology could be improved, with the combination of current and future available data in gravitational wave (GW) and electromagnetic (EM) domain.

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.

So far most studies on the structure of coronal mass ejections (CMEs) are conducted through white-light coronagraphs, which demonstrate about one third of CMEs exhibit the typical three-part structure in the high corona (e.g., beyond 2 Rs), i.e., the bright front, the dark cavity and the bright core. In this paper, we address the CME structure in the low corona (e.g., below 1.3 Rs) through extreme-ultraviolet (EUV) passbands and find that the three-part CMEs in the white-light images can possess a similar three-part appearance in the EUV images, i.e., a leading edge, a low-density zone, and a filament or hot channel. The analyses identify that the leading edge and the filament or hot channel in the EUV passbands evolve into the front and the core later within several solar radii in the white-light passbands, respectively. What's more, we find that the CMEs without obvious cavity in the white-light images can also exhibit the clear three-part appearance in the EUV images, which means that the low-density zone in the EUV images (observed as the cavity in white-light images) can be compressed and/or transformed gradually by the expansion of the bright core and/or the reconnection of magnetic field surrounding the core during the CME propagation outward. Our study suggests that more CMEs can possess the clear three-part structure in their early eruption stage. The nature of the low-density zone between the leading edge and the filament or hot channel is discussed.

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.

During magnetically dominated relativistic reconnection, inflowing plasma depletes the initial relativistic pressure at the x-line and collisionless plasma heating inside the diffusion region is insufficient to overcome this loss. The resulting pressure drop causes a collapse at the x-line, essentially a localization mechanism of the diffusion region necessary for fast reconnection. The extension of this low-pressure region further explains the bursty nature of anti-parallel reconnection because a once opened outflow exhaust can also collapse, which repeatedly triggers secondary tearing islands. However, a stable single x-line reconnection can be achieved when an external guide field exists, since the reconnecting magnetic field component rotates out of the reconnection plane at outflows, providing additional magnetic pressure to sustain the opened exhausts.

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

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

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

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

The Kelvin-Helmholtz (KH) instability is commonly found in many astrophysical, laboratory, and space plasmas. It could mix plasma components of different properties and convert dynamic fluid energy from large scale structure to smaller ones. In this study, we combined the ground-based New Vacuum Solar Telescope (NVST) and the Solar Dynamic Observatories (SDO) / Atmospheric Imaging Assembly (AIA) to observe the plasma dynamics associated with active region 12673 on 09 September 2017. In this multi-temperature view, we identified three adjacent layers of plasma flowing at different speeds, and detected KH instabilities at their interfaces. We could unambiguously track a typical KH vortex and measure its motion. We found that the speed of this vortex suddenly tripled at a certain stage. This acceleration was synchronized with the enhancements in emission measure and average intensity of the 193 \AA{} data. We interpret this as evidence that KH instability triggers plasma heating. The intriguing feature in this event is that the KH instability observed in the NVST channel was nearly complementary to that in the AIA 193 \AA{}. Such a multi-thermal energy exchange process is easily overlooked in previous studies, as the cold plasma component is usually not visible in the extreme ultraviolet channels that are only sensitive to high temperature plasma emissions. Our finding indicates that embedded cold layers could interact with hot plasma as invisible matters. We speculate that this process could occur at a variety of length scales and could contribute to plasma heating.

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

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

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

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

Strongly gravitational lensing systems (SGL) encodes cosmology information in source/lens distance ratios $\mathcal{D}_{\rm obs}=\mathcal{D}_{\rm ls}/\mathcal{D}_{\rm s}$, which can be used to precisely constrain cosmological parameters. In this paper, based on future measurements of 390 strong lensing systems from the forthcoming LSST survey, we have successfully reconstructed the distance ratio $\mathcal{D}_{\rm obs}$ (with the source redshift $z_s\sim 4.0$), directly from the data without assuming any parametric form. A recently developed method based on model-independent reconstruction approach, Gaussian Processes (GP) is used in our study of these strong lensing systems. Our results show that independent measurement of the matter density parameter ($\Omega_m$) could be expected from such strong lensing statistics. More specifically, one can expect $\Omega_m$ to be estimated at the precision of $\Delta\Omega_m\sim0.015$ in the concordance $\Lambda$CDM model, which provides comparable constraint on $\Omega_m$ with Planck 2015 results. In the framework of the modified gravity theory (DGP), 390 detectable galactic lenses from future LSST survey would lead to stringent fits of $\Delta\Omega_m\sim0.030$. \textbf{Finally, we have discussed three possible sources of systematic errors (sample incompleteness, the determination of length of lens redshift bin, and the choice of lens redshift shells), and quantified their effects on the final cosmological constraints. Our results strongly indicate that future strong lensing surveys, with the accumulation of a larger and more accurate sample of detectable galactic lenses, will considerably benefit from the methodology described in this analysis. }

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

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

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

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.

Understanding planet formation requires to discern how dust grows in protoplanetary disks. An important parameter to measure in disks is the maximum dust grain size present. This is usually estimated through measurements of the dust opacity at different millimeter wavelengths assuming optically thin emission and dust opacity dominated by absorption. However, ALMA observations have shown that these assumptions might not be correct in the case of protoplanetary disks, leading to overestimation of particle sizes and to underestimation of the disk's mass. Here, we present an analysis of high quality ALMA and VLA images of the HL Tau protoplanetary disk, covering a wide range of wavelengths, from 0.8 mm to 1 cm, and with a physical resolution of $\sim$7.35 au. We describe a procedure to analyze a set of millimeter images without any assumption about the optical depth of the emission, and including the effects of absorption and scattering in the dust opacity. This procedure allows us to obtain the dust temperature, the dust surface density and the maximum particle size at each radius. In the HL Tau disk, we found that particles have already grown up to a few millimeters in size. We detect differences in the dust properties between dark and bright rings, with dark rings containing low dust density and small dust particles. Different features in the HL Tau disk seem to have different origins. Planet-disk interactions can explain substructure at the external half of the disk, but the internal rings seem to be associated to the presence of snow lines of several molecules.

Kilonova explosions typically release $\sim 10^{50-51}$ erg in kinetic energy, which is sufficient to constitute a kilonova remnant (KNR); however, it has not yet been confirmed. In this work, we investigate the probable association between G4.8+6.2 and the guest star of AD 1163, which is recorded by the Korea ancient astronomers. Although the evidence available is insufficient to draw a definite conclusion, it is at least theoretically self-consistent that the guest star of AD 1163 was a historical kilonova associated with G4.8+6.2, considering the possible short visible timescale of AD 1163, the relatively high Galactic latitude of G4.8+6.2, and that G4.8+6.2 is spatially coincident with the guest star of AD 1163. Further observation of G4.8+6.2 is needed to test our hypothesis. If our interpretation is correct, our results indicate that young KNRs should have a large diameter and low surface brightness, unlike other young supernova remnants.

Turbulence is ubiquitously observed in nearly collisionless heliospheric plasmas, including the solar wind and corona and the Earth's magnetosphere. Understanding the collisionless mechanisms responsible for the energy transfer from the turbulent fluctuations to the particles is a frontier in kinetic turbulence research. Collisionless energy transfer from the turbulence to the particles can take place reversibly, resulting in non-thermal energy in the particle velocity distribution functions (VDFs) before eventual collisional thermalization is realized. Exploiting the information contained in the fluctuations in the VDFs is valuable. Here we apply a recently developed method based on VDFs, the field-particle correlation technique, to a $\beta=1$, solar-wind-like, low-frequency Alfv\'enic turbulence simulation with well resolved phase space to identify the field-particle energy transfer in velocity space. The field-particle correlations reveal that the energy transfer, mediated by the parallel electric field, results in significant structuring of the ion and electron VDFs in the direction parallel to the magnetic field. Fourier modes representing the length scales between the ion and electron gyroradii show that energy transfer is resonant in nature, localized in velocity space to the Landau resonances for each Fourier mode. The energy transfer closely follows the Landau resonant velocities with varying perpendicular wavenumber $k_\perp$ and plasma $\beta$. This resonant signature, consistent with Landau damping, is observed in all diagnosed Fourier modes that cover the dissipation range of the simulation.

An unexpected strong geomagnetic storm occurred on 2018 August 26, which was caused by a slow coronal mass ejection (CME) from a gradual eruption of a large quiet-region filament. We investigate the eruption and propagation characteristics of this CME in relation to the strong geomagnetic storm with remote sensing and in situ observations. Coronal magnetic fields around the filament are extrapolated and compared with EUV observations. We determine the propagation direction and tilt angle of the CME flux rope near the Sun using a graduated cylindrical shell (GCS) model and the Sun-to-Earth kinematics of the CME with wide-angle imaging observations from STEREO A. We reconstruct the flux-rope structure using a Grad-Shafranov technique based on the in situ measurements at the Earth and compare it with those from solar observations and the GCS results. Our conclusions are as follows: (1) the eruption of the filament was unusually slow and occurred in the regions with relatively low critical heights of the coronal field decay index; (2) the axis of the CME flux rope rotated in the corona as well as in interplanetary space, which tended to be aligned with the local heliospheric current sheet; (3) the CME was bracketed between slow and fast solar winds, which enhanced the magnetic field inside the CME at 1 AU; (4) the geomagnetic storm was caused by the enhanced magnetic field and a southward orientation of the flux rope at 1 AU from the rotation of the flux rope.

The ubiquitous solar jets or jet-like activities are generally regarded as an important source of energy and mass input to the upper solar atmosphere and the solar wind. However, questions about their triggering and driving mechanisms are not completely understood. By taking advantage of high temporal and high spatial resolution stereoscopic observations taken by the Solar Dynamic Observatory (SDO) and the Solar Terrestrial Relations Observatory (STEREO), we report an intriguing two-sided-loop jet occurred on 2013 June 02, which was dynamically associated with the eruption of a mini-filament below an overlying large filament, and two distinct reconnection processes are identified during the formation stage. The SDO observations reveals that the two-sided-loop jet showed a concave shape with a projection speed of about 80 - 136. From the other view angle, the STEREO ahead observations clearly showed that the trajectory of the two arms of the two-sided-loop were along the cavity magnetic field lines hosting the large filament. Contrary to the well-accepted theoretical model, the present observation sheds new light on our understanding of the formation mechanism of two-sided-loop jets. Moreover, the eruption of the two-sided-loop jet not only supplied mass to the overlying large filament, but also provided a rare opportunity to diagnose the magnetic structure of the overlying large filament via the method of three-dimensional reconstruction.

In order to have a comprehensive view of the propagation and evolution of coronal mass ejections (CMEs) from the Sun to deep interplanetary space beyond 1 au, we carry out a kinematic analysis of 7 CMEs in solar cycle 23. The events are required to have coordinated coronagraph observations, interplanetary type II radio bursts, and multi-point in-situ measurements at the Earth and Ulysses. A graduated cylindrical shell model, an analytical model without free parameters and a magnetohydrodynamic model are used to derive CME kinematics near the Sun, to quantify the CME/shock propagation in the Sun-Earth space, and to connect in-situ signatures at the Earth and Ulysses, respectively. We find that each of the 7 CME-driven shocks experienced a major deceleration before reaching 1 au and thereafter propagated with a gradual deceleration from the Earth to larger distances. The resulting CME/shock propagation profile for each case is roughly consistent with all the data, which verifies the usefulness of the simple analytical model for CME/shock propagation in the heliosphere. The statistical analysis of CME kinematics indicates a tendency that the faster the CME, the larger the deceleration, and the shorter the deceleration time period within 1 au. For several of these events, the associated geomagnetic storms were mainly caused by the southward magnetic fields in the sheath region. In particular, the interaction between a CME-driven shock and a preceding ejecta significantly enhanced the preexisting southward magnetic fields and gave rise to a severe complex geomagnetic storm.

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

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

The impact of stellar multiplicity on the evolution of planet-forming disks is still the subject of debate. Here we present and analyze disk structures around ten multiple stellar systems that were included in an unbiased, high spatial resolution survey performed with ALMA of 32 protoplanetary disks in the Taurus star-forming region. At the unprecedented spatial resolution of ~0.12" we detect and spatially resolve the disks around all primary stars, and those around eight secondary and one tertiary star. The dust radii of disks around multiple stellar systems are smaller than those around single stars in the same stellar mass range and in the same region. The disks in multiple stellar systems also show a steeper decay of the millimeter continuum emission at the outer radius than disks around single stars, suggestive of the impact of tidal truncation on the shape of the disks in multiple systems. However, the observed ratio between the dust disk radii and the observed separation of the stars in the multiple systems is consistent with analytic predictions of the effect of tidal truncation only if the eccentricities of the binaries are rather high (typically >0.5), or if the observed dust radii are a factor of two smaller than the gas radii, as is typical for isolated systems. Similar high-resolution studies targeting the gaseous emission from disks in multiple stellar systems are required to resolve this question.

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

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

The spreading of the X-line out of the reconnection plane under a strong guide field is investigated using large-scale three-dimensional (3D) particle-in-cell (PIC) simulations in asymmetric magnetic reconnection. A simulation with a thick, ion-scale equilibrium current sheet (CS) reveals that the X-line spreads at the ambient ion/electron drift speeds, significantly slower than the Alfv\'en speed based on the guide field $V_{Ag}$. Additional simulations with a thinner, sub-ion-scale CS show that the X-line spreads at $V_{Ag}$ (Alfv\'enic spreading), much higher than the ambient species drifts. An Alfv\'enic signal consistent with kinetic Alfv\'en waves develops and propagates, leading to CS thinning and extending, which then results in reconnection onset. The continuous onset of reconnection in the signal propagation direction manifests as Alfv\'enic X-line spreading. The strong dependence on the CS thickness of the spreading speeds, and the X-line orientation are consistent with the collisionless tearing instability. Our simulations indicate that when the collisionless tearing growth is sufficiently strong in a thinner CS such that $\gamma/\Omega_{ci}\gtrsim\mathcal{O}(1)$, Alfv\'enic X-line spreading can take place. Our results compare favorably with a number of numerical simulations and recent magnetopause observations. A key implications is that the magnetopause CS is typically too thick for Alfv\'enic X-line spreading to effectively take place.

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.

We present a high-resolution ($\sim0.''12$, $\sim16$ au, mean sensitivity of $50~\mu$Jy~beam$^{-1}$ at 225 GHz) snapshot survey of 32 protoplanetary disks around young stars with spectral type earlier than M3 in the Taurus star-forming region using Atacama Large Millimeter Array (ALMA). This sample includes most mid-infrared excess members that were not previously imaged at high spatial resolution, excluding close binaries and highly extincted objects, thereby providing a more representative look at disk properties at 1--2 Myr. Our 1.3 mm continuum maps reveal 12 disks with prominent dust gaps and rings, 2 of which are around primary stars in wide binaries, and 20 disks with no resolved features at the observed resolution (hereafter smooth disks), 8 of which are around the primary star in wide binaries. The smooth disks were classified based on their lack of resolved substructures, but their most prominent property is that they are all compact with small effective emission radii ($R_{\rm eff,95\%} \lesssim 50$ au). In contrast, all disks with $R_{\rm eff,95\%}$ of at least 55 au in our sample show detectable substructures. Nevertheless, their inner emission cores (inside the resolved gaps) have similar peak brightness, power law profiles, and transition radii to the compact smooth disks, so the primary difference between these two categories is the lack of outer substructures in the latter. These compact disks may lose their outer disk through fast radial drift without dust trapping, or they might be born with small sizes. The compact dust disks, as well as the inner disk cores of extended ring disks, that look smooth at the current resolution will likely show small-scale or low-contrast substructures at higher resolution. The correlation between disk size and disk luminosity correlation demonstrates that some of the compact disks are optically thick at millimeter wavelengths.

We investigate a global extreme-ultraviolet (EUV) wave associated with a coronal mass ejection (CME)-driven shock on 2017 September 10. The EUV wave is transmitted by north- and south-polar coronal holes (CHs), which is observed by the Solar Dynamics Observatory (SDO) and Solar Terrestrial Relations Observatory A (STEREO-A) from opposite sides of the Sun. We obtain key findings on how the EUV wave interacts with multiple coronal structures, and on its connection with the CME-driven shock: (1) the transmitted EUV wave is still connected with the shock that is incurvated to the Sun, after the shock has reached the opposite side of the eruption; (2) the south CH transmitted EUV wave is accelerated inside an on-disk, low-density region with closed magnetic fields, which implies that an EUV wave can be accelerated in both open and closed magnetic field regions; (3) part of the primary EUV wavefront turns around a bright point (BP) with a bipolar magnetic structure when it approaches a dim, low-density filament channel near the BP; (4) the primary EUV wave is diffused and apparently halted near the boundaries of remote active regions (ARs) that are far from the eruption, and no obvious AR related secondary waves are detected; (5) the EUV wave extends to an unprecedented scale of ~360{\deg} in latitudes, which is attributed to the polar CH transmission. These results provide insights into the effects of coronal density and magnetic field distributions on the evolution of an EUV wave, and into the connection between the EUV wave and the associated CME-driven shock.

While a growing body of research indicates that relativistic magnetic reconnection is a prodigious source of particle acceleration in high-energy astrophysical systems, the dominant acceleration mechanism remains controversial. Using a combination of fully kinetic simulations and theoretical analysis, we demonstrate that Fermi-type acceleration within the large-scale motional electric fields dominates over direct acceleration from non-ideal electric fields within small-scale diffusion regions. This result has profound implications for modeling particle acceleration in large-scale astrophysical problems, since it opens up the possiblity of modeling the energetic spectra without resolving microscopic diffusion regions.

By analyzing two sets of complete $BVR_cI_c$ light curves of V342 UMa and three sets of complete $BVR_cI_c$ light curves of V509 Cam, we determined that the two systems are both W-subtype contact binaries and that V342 UMa shows a shallow contact configuration, while V509 Cam exhibits a medium contact configuration. Since both of them are totally eclipsing binaries, the physical parameters derived only by the photometric light curves are reliable.Meanwhile, the period changes of the two targets were analyzed based on all available eclipsing times. We discovered that V342 UMa shows long-term period decrease with a rate of $-1.02(\pm0.54)\times10^{-7}$ days/year and that V509 Cam exhibits long-term period increase with a rate of $3.96(\pm0.90)\times10^{-8}$ days/year. Both the conservative mass transfer and AML via magnetic stellar winds can interpret the long-term period decrease of V342 UMa. The long-term period increase of V509 Cam can be explained by mass transfer from the less massive star to the more massive one. The absolute parameters of the two binaries were estimated according to the Gaia distances and our derived photometric solution results. This method can be extended to other contact binaries without radial velocities but with reliable photometric solutions. The evolutionary states of them were discussed, we found that they reveal identical properties of other W-subtype contact systems.

In this paper,we report a detailed analysis of recurrent jets originated from a location with emerging, canceling and converging negative magnetic field at the east edge of NOAA active region AR11166 from 2011 March 09 to 10. The event presented several interesting features. First, a satellite sunspot appeared and collided with a pre-existing opposite polarity magnetic field and caused a recurrent solar jet event. Second, the evolution of the jets showed blowout-like nature and standard characteristics. Third, the satellite sunspot exhibited a motion toward southeast of AR11166 and merged with the emerging flux near the opposite polarity sunspot penumbra, which afterward, due to flux convergence and cancellation episodes, caused recurrent jets. Fourth, three of the blowout jets associated with coronal mass ejections (CMEs), were observed from field of view of the Solar Terrestrial Relations Observatory. Fifth, almost all the blowout jet eruptions were accompanied with flares or with more intense brightening in the jet base region, while almost standard jets did not manifest such obvious feature during eruptions. The most important, the blowout jets were inclined to faster and larger scale than the standard jets. The standard jets instead were inclined to relative longer-lasting. The obvious shearing and twisting motions of the magnetic field may be interpreted as due to the shearing and twisting motions for a blowout jet eruption. From the statistical results, about 30% blowout jets directly developed into CMEs. It suggests that the blowout jets and CMEs should have a tight relationship.

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

We propose a novel approach to reconstruct the surface magnetic helicity density on the Sun or sun-like stars. The magnetic vector potential is determined via decomposition of vector magnetic field measurements into toroidal and poloidal components. The method is verified using data from a non-axisymmetric dynamo model. We apply the method to vector field synoptic maps from Helioseismic and Magnetic Imager (HMI) onboard of Solar Dynamics Observatory (SDO) to study evolution of the magnetic helicity density during solar cycle 24. It is found that the mean helicity density of the non-axisymmetric magnetic field of the Sun evolves in a way which is similar to that reported for the current helicity density of the solar active regions. It has predominantly the negative sign in the northern hemisphere, and it is positive in the southern hemisphere. Also, the hemispheric helicity rule for the non-axisymmetric magnetic field showed the sign inversion at the end of cycle 24. Evolution of magnetic helicity density of large-scale axisymmetric magnetic field is different from that expected in dynamo theory. On one hand, the mean large- and small-scale components of magnetic helicity density display the hemispheric helicity rule of opposite sign at the beginning of cycle 24. However, later in the cycle, the two helicities exhibit the same sign in contrast with the theoretical expectations.

We investigate the chemical potential and baryon number density of the hadron-quark phase transition in neutron star matter. The hadron matter is described with relativistic mean field theory, and the quark matter is described with the Dyson-Schwinger equation approach of QCD. In order to study the first-order phase transition, we develop the sound speed interpolation scheme to construct the equation of state in the middle density region where the hadron phase and quark phase coexist. The phase transition chemical potential is constrained with the maximum mass, the tidal deformability and the radius of neutrons stars. And the most probable value of the phase transition chemical potential is found.

We investigate whether the current Hubble parameter $H(z)$ measurements could help improve the constraints on dark energy on the basis of the mainstream cosmological probes including the type Ia supernovae (SN) observation, the cosmic microwave background anisotropies (CMB) observation, and the baryon acoustic oscillations (BAO) observation. For the current $H(z)$ data, we use 30 data points measured by using a differential age method. Furthermore, we also consider the future $H(z)$ measurements based on the Sandage-Loeb (SL) test by means of the E-ELT in construction, and thus we also use 30 simulated $H(z)$ data according to a 10-year SL test observation. In this work, we choose four typical dark energy models as examples, i.e., the $\Lambda$CDM model, the $w$CDM model, the $\alpha$DE model, and the GCG model, to complete the analysis. We find that, when only the current $H(z)$ data are added, the constraints on these models are not improved compared to the cases using the SN+CMB+BAO data; but when further adding the 10-year SL test data, the constraint results are tremendously improved for all the four models. Therefore, we conclude that, although the current $H(z)$ measurements could not provide an evident improvement on the basis of the current mainstream cosmological probes, the future $H(z)$ measurements from the SL test would have enormous potential to change the status of the Hubble parameter measurements in constraining dark energy.

PDS 70b is the most robustly detected young planet imaged in the gap of a transition disk so far, found at a separation of ~195 mas (~22 au) from its host star and at a position angle of about 155 deg. This system is therefore a unique laboratory to characterize the properties of young planetary systems at the stage of their formation. We aim to trace direct and indirect imprints of PDS 70b on the gas and dust emission of the circumstellar disk in order to study the properties of this ~5 Myr young planetary system. We obtained ALMA band 7 observations of PDS 70 in dust continuum and $^{12}CO$ (3-2) and combined them with archival data resulting in an unprecedented angular resolution of about 70 mas (~8 au). We derive an upper limit on circumplanetary material at the location of PDS 70b of ~0.01 $M_{\oplus}$ and find a highly structured circumstellar disk in both dust and gas. The outer dust ring peaks at ~0.65" (74 au) and reveals a possible second unresolved peak at ~0.53" (60 au). The integrated intensity of CO also shows evidence of a depletion of emission at ~0.2" (23 au) with a width of ~0.1" (11 au). The gas kinematics show evidence of a deviation from Keplerian rotation inside $\lesssim$0.8" (91 au). This implies a pressure gradient that can account for the location of the dust ring well beyond the location of PDS 70b. Further in, we detect an inner disk which appears to be connected to the outer disk by a possible bridge feature in the North-West region in both gas and dust. We compare the observations to hydrodynamical simulations including a planet with different masses covering the estimated mass range previously derived from NIR photometry (~5-9 M$_{Jup}$). We find that even a planet with a mass of 10 M$_{Jup}$ may not be sufficient to explain the extent of the wide gap and an additional low-mass companion may be needed to account for the observed disk morphology.

The globular cluster (GC) specific frequency ($S_N$), defined as the number of GCs per unit galactic luminosity, represents the efficiency of GC formation (and survival) compared to field stars. Despite the naive expectation that star cluster formation should scale directly with star formation, this efficiency varies widely across galaxies. To explore this variation we measure the z-band GC specific frequency ($S_{N,z}$) for 43 early-type galaxies (ETGs) from the Hubble Space Telescope (HST)/Advanced Camera for Surveys (ACS) Fornax Cluster Survey. Combined with the homogenous measurements of $S_{N,z}$ in 100 ETGs from the HST/ACS Virgo Cluster Survey from Peng et al. (2008), we investigate the dependence of $S_{N,z}$ on mass and environment over a range of galaxy properties. We find that $S_{N,z}$ behaves similarly in the two galaxy clusters, despite the clusters' order-of-magnitude difference in mass density. The $S_{N,z}$ is low in intermediate-mass ETGs ($-20<M_z<-23$), and increases with galaxy luminosity. It is elevated at low masses, on average, but with a large scatter driven by galaxies in dense environments. The densest environments with the strongest tidal forces appear to strip the GC systems of low-mass galaxies. However, in low-mass galaxies that are not in strong tidal fields, denser environments correlate with enhanced GC formation efficiencies. Normalizing by inferred halo masses, the GC mass fraction, $\eta=(3.36\pm0.2)\times10^{-5}$, is constant for ETGs with stellar masses $\mathcal{M}_\star \lesssim 3\times10^{10}M_\odot$, in agreement with previous studies. The lack of correlation between the fraction of GCs and the nuclear light implies only a weak link between the infall of GCs and the formation of nuclei.

We propose a scheme to determine the chemical potential and baryon number density of the hadron-quark phase transition in cold dense strong interaction matter (compact star matter). The hadron matter is described with the relativistic mean field theory, and the quark matter is described with the Dyson-Schwinger equation approach of QCD. To study the first-order phase transition, we take the sound speed as the interpolation objective to construct the equation of state in the middle density region. With the maximum mass, the tidal deformability and the radius of neutron stars being taken as calibration quantities, the phase transition chemical potential is constrained to a quite small range. And the most probable value of the phase transition chemical potential is found.

ALMA has observed a plethora of ring-like structures in planet forming discs at distances of 10-100 au from their host star. Although several mechanisms have been invoked to explain the origin of such rings, a common explanation is that they trace new-born planets. Under the planetary hypothesis, a natural question is how to reconcile the apparently high frequency of gap-carving planets at 10-100 au with the paucity of Jupiter mass planets observed around main sequence stars at those separations. Here, we provide an analysis of the new-born planet population emerging from observations of gaps in discs, under the assumption that the observed gaps are due to planets. We use a simple estimate of the planet mass based on the gap morphology, and apply it to a sample of gaps recently obtained by us in a survey of Taurus with ALMA. We also include additional data from recent published surveys, thus analysing the largest gap sample to date, for a total of 48 gaps. The properties of the purported planets occupy a distinctively different region of parameter space with respect to the known exo-planet population, currently not accessible through planet finding methods. Thus, no discrepancy in the mass and radius distribution of the two populations can be claimed at this stage. We show that the mass of the inferred planets conforms to the theoretically expected trend for the minimum planet mass needed to carve a dust gap. Finally, we estimate the separation and mass of the putative planets after accounting for migration and accretion, for a range of evolutionary times, finding a good match with the distribution of cold Jupiters.

In this paper we present a curated dataset from the NASA Solar Dynamics Observatory (SDO) mission in a format suitable for machine learning research. Beginning from level 1 scientific products we have processed various instrumental corrections, downsampled to manageable spatial and temporal resolutions, and synchronized observations spatially and temporally. We illustrate the use of this dataset with two example applications: forecasting future EVE irradiance from present EVE irradiance and translating HMI observations into AIA observations. For each application we provide metrics and baselines for future model comparison. We anticipate this curated dataset will facilitate machine learning research in heliophysics and the physical sciences generally, increasing the scientific return of the SDO mission. This work is a direct result of the 2018 NASA Frontier Development Laboratory Program. Please see the appendix for access to the dataset.

We collected the information of 304 gamma-ray bursts (GRBs) from the literature, and analyzed the correlations among the host galaxy offsets (the distance from the site of the GRB to the center of its host galaxy), $T_{\rm 90,i}$ (the duration $T_{\rm 90}$ in rest-frame), $T_{\rm R45,i}$ (the duration $T_{\rm R45}$ in rest-frame), $E_{\rm \gamma,iso}$ (the isotropic equivalent energy), $L_{\rm \gamma,iso}$ ($=E_{\rm \gamma,iso}/T_{\rm 90,i}$, the isotropic equivalent luminosity) and $L_{\rm pk}$ (peak luminosity). We found that $ T_{\rm 90,i}$, $ T_{\rm R45,i}$, $ E_{\rm \gamma,iso}$, $L_{\rm pk}$ have negative correlation with $ {\rm offset}$, which is consistent with origin of short GRBs (SGRBs) and long GRBs (LGRBs). On separate analysis, we found similar results for $\log E_{\rm \gamma,iso}$ - $\log {\rm (offset)}$ and $\log L_{\rm pk}$ - $\log {\rm (offset)}$ relations in case of SGRBs only, while no obvious relation for LGRBs. There is no correlations between offset and $L_{\rm \gamma,iso}$. We also put the special GRB 170817A {and GRB 060218A} on the plots. {The two GRBs both have low luminosity and small offset.} In the $ \log ({\rm offset})- \log T_{\rm 90,i}$ plot, we found GRB 170817A locates in between the two regions of SGRBs and LGRBs and it is the outlier in the $ {\rm offset}-E_{\rm \gamma, iso}$, $ {\rm offset}-L_{\rm \gamma, iso}$ and $ {\rm offset}-L_{\rm pk}$ plots. Together with GRB 060218A being an outlier in all plots, it indicates the speciality of GRBs 170817A and 060218A, and might imply more subgroups of the GRB samples.

In order to obtain an overview of the gamma-ray bursts (GRBs), we need a full sample. In this paper, we collected 6289 GRBs (from GRB 910421 to GRB 160509A) from the literature, including prompt emission, afterglow and host galaxy properties. We hope to use this large sample to reveal the intrinsic properties of GRB. We have listed all the data in machine readable tables, including the properties of the GRBs, correlation coefficients and linear regression results of two arbitrary parameters, and linear regression results of any three parameters. These machine readable tables could be used as a data reservoir for further studies on the classifications or correlations. One may find some intrinsic properties from these statistical results. With this comprehensive table, it is possible to find relations between different parameters, and to classify the GRBs into different kinds of sub-groups. With the completion, it may reveal the nature of GRBs and may be used as tools like pseudo-redshift indicators, standard candles, etc. All the machine readable data and statistical results are available on the website of the journal.

Infrared spectroscopy with medium to high spectral resolution is essential to characterize the gas content of circumstellar disks. Unfortunately, conducting continuum and line radiative transfer of thermochemical disk models is too time-consuming to carry out large parameter studies. Simpler approaches using a slab model to fit continuum-subtracted spectra require the identification of either the global or local continuum. Continuum subtraction, particularly when covering a broad wavelength range, is challenging but critical in rich molecular spectra as hot (several hundreds K) molecular emission lines can also produce a pseudo continuum. In this work, we present CLIcK, a flexible tool to simultaneously fit the continuum and line emission. The DDN01 continuum model (Dullemond et al. 2001) and a plane-parallel slab of gas in local thermodynamic equilibrium are adopted to simulate the continuum and line emission respectively, both of them are fast enough for homogeneous studies of large disk samples. We applied CLIcK to fit the observed water spectrum of the AA Tau disk and obtained water vapor properties that are consistent with literature results. We also demonstrate that CLIcK properly retrieves the input parameters used to simulate the water spectrum of a circumstellar disk. CLIcK will be a versatile tool for the interpretation of future James Webb Space Telescope spectra.

It is paramount from both scientific and societal perspectives to understand the generation of extreme space weather. We discuss the formation of solar superstorms based on a comparative study of the 2012 July 23 and 2017 July 23 eruptions. The first one is Carrington-class, and the second could rival the 1989 March event that caused the most intense geomagnetic storm of the space age. Observations of these events in the historically weak solar cycle 24 indicate that a solar superstorm can occur in any solar cycle and at any phase of the cycle. Recurrent patterns are identified in both cases, including the long-lived eruptive nature of the active region, a complex event composed of successive eruptions from the same active region, and in-transit interaction between the successive eruptions resulting in exceptionally strong ejecta magnetic fields at 1 AU. Each case also shows unique characteristics. Preconditioning of the upstream solar wind leading to unusually high solar wind speeds at 1 AU is observed in the first case whereas not in the latter. This may suggest that the concept of "preconditioning" appears to be necessary for making a Carrington-class storm. We find a considerable deflection by nearby coronal holes in the second case but not in the first. On the basis of these results, we propose a hypothesis for further investigation that superstorms are "perfect storms" in nature, i.e., a combination of circumstances that results in an event of unusual magnitude. Historical records of some extreme events seem to support our hypothesis.

Using 3D particle-in-cell (PIC) simulations, we study magnetic reconnection with the x-line being spatially confined in the current direction. We include thick current layers to prevent reconnection at two ends of a thin current sheet that has a thickness on an ion inertial (di) scale. The reconnection rate and outflow speed drop significantly when the extent of the thin current sheet in the current direction is < O(10 di). When the thin current sheet extent is long enough, we find it consists of two distinct regions; an inactive region (on the ion-drifting side) exists adjacent to the active region where reconnection proceeds normally as in a 2D case. The extent of this inactive region is ~ O(10 di), and it suppresses reconnection when the thin current sheet extent is comparable or shorter. The time-scale of current sheet thinning toward fast reconnection can be translated into the spatial-scale of this inactive region; because electron drifts inside the ion diffusion region transport the reconnected magnetic flux, that drives outflows and furthers the current sheet thinning, away from this region. This is a consequence of the Hall effect in 3D. While this inactive region may explain the shortest possible azimuthal extent of dipolarizing flux bundles at Earth, it may also explain the dawn-dusk asymmetry observed at the magnetotail of Mercury, that has a global dawn-dusk extent much shorter than that of Earth.

In this paper, we assemble a well-defined sample of early-type gravitational lenses extracted from a large collection of 158 systems, and use the redshift distribution of galactic-scale lenses to test the standard cosmological model ($\Lambda$CDM) and the modified gravity theory (DGP). Two additional sub-samples are also included to account for possible selection effect introduced by the detectability of lens galaxies. Our results show that independent measurement of the matter density parameter ($\Omega_m$) could be expected from such strong lensing statistics. Based on future measurements of strong lensing systems from the forthcoming LSST survey, one can expect $\Omega_m$ to be estimated at the precision of $\Delta\Omega_m\sim 0.006$, which provides a better constraint on $\Omega_m$ than \textit{Planck} 2015 results. Moreover, use the lens redshift test is also used to constrain the characteristic velocity dispersion of the lensing galaxies, which is well consistent with that derived from the optical spectroscopic observations. A parameter $f_E$ is adopted to quantify the relation between the lensing-based velocity dispersion and the corresponding stellar value. Finally, the accumulation of detectable galactic lenses from future LSST survey would lead to more stringent fits of $\Delta f_E\sim10^{-3}$, which encourages us to test the global properties of early-type galaxies at much higher accuracy.

We report the first unambiguous quasi-periodic large-scale extreme-ultraviolet (EUV) wave or shock that was detected by the Atmospheric Imaging Assembly on board the Solar Dynamics Observatory. During the whip-like unwinding eruption of a small filament on 2012 April 24, multiple consecutive large-scale wavefronts emanating from AR11467 were observed simultaneously along the solar surface and a closed transequatorial loop system. In the meantime, an upward propagating dome-shaped wavefront was also observed, whose initial speed and deceleration are about 1392 km/s and 1.78 km/s^2, respectively. Along the solar surface, the quasi-peridoic wavefronts had a period of about 163 +/- 21 seconds and propagated at a nearly constant speed of 747 +/- 26 km/s; they interacted with active region AR11469 and launched a sympathetic upward propagating secondary EUV wave. The wavefronts along the loop system propagated at a speed of 897 km/s, and they were reflected back at the southern end of the loop system at a similar speed. In addition to the propagating waves, a standing kink wave was also present in the loop system simultaneously. Periodicity analysis reveals that the period of the wavefronts was consistent with that of the unwinding helical structures of the erupting filament. Based on these observational facts, we propose that the observed quasi-periodic EUV wavefronts were most likely excited by the periodic unwinding motion of the filament helical structures. In addition, two different seismological methods are applied to derive the magnetic field strength of the loop system, and for the first time the reliability of these inversion techniques are tested with the same magnetic structure.

We report the first unambiguous quasi-periodic large-scale extreme-ultraviolet (EUV) wave or shock that was detected by the Atmospheric Imaging Assembly on board the Solar Dynamics Observatory. During the whip-like unwinding eruption of a small filament on 2012 April 24, multiple consecutive large-scale wavefronts emanating from AR11467 were observed simultaneously along the solar surface and a closed transequatorial loop system. In the meantime, an upward propagating dome-shaped wavefront was also observed, whose initial speed and deceleration are about 1392 km/s and 1.78 km/s^2, respectively. Along the solar surface, the quasi-peridoic wavefronts had a period of about 163 +/- 21 seconds and propagated at a nearly constant speed of 747 +/- 26 km/s; they interacted with active region AR11469 and launched a sympathetic upward propagating secondary EUV wave. The wavefronts along the loop system propagated at a speed of 897 km/s, and they were reflected back at the southern end of the loop system at a similar speed. In addition to the propagating waves, a standing kink wave was also present in the loop system simultaneously. Periodicity analysis reveals that the period of the wavefronts was consistent with that of the unwinding helical structures of the erupting filament. Based on these observational facts, we propose that the observed quasi-periodic EUV wavefronts were most likely excited by the periodic unwinding motion of the filament helical structures. In addition, two different seismological methods are applied to derive the magnetic field strength of the loop system, and for the first time the reliability of these inversion techniques are tested with the same magnetic structure.

KAGRA is a second-generation interferometric gravitational-wave detector with 3-km arms constructed at Kamioka, Gifu in Japan. It is now in its final installation phase, which we call bKAGRA (baseline KAGRA), with scientific observations expected to begin in late 2019. One of the advantages of KAGRA is its underground location of at least 200 m below the ground surface, which brings small seismic motion at low frequencies and high stability of the detector. Another advantage is that it cools down the sapphire test mass mirrors to cryogenic temperatures to reduce thermal noise. In April-May 2018, we have operated a 3-km Michelson interferometer with a cryogenic test mass for 10 days, which was the first time that km-scale interferometer was operated at cryogenic temperatures. In this article, we report the results of this "bKAGRA Phase 1" operation. We have demonstrated the feasibility of 3-km interferometer alignment and control with cryogenic mirrors.

The Earths core formation process has decisive effect in the chemical differentiation between the Earths core and its mantle. Here, we propose a new core formation model which is caused by a special giant impact. This model suggests that the impactors core can be kept intact by its own sticky mantle under appropriate impacting conditions and let it merge into the targets core without contact with the targets mantle. We call this special giant impact that caused the new core formation mode as glue ball impact model (GBI). By simulating hundreds of giant impacts with the sizes from planetesimals to planets, the conditions that can lead to GBI have been found out. If with small impact angle (i.e., less than 20 degree), small impact velocity and small impactors mass but larger than 0.07 Mearth, there is a good chance to produce a GBI at the final stage of the Earths accretion. We find that it will be much easier to have GBIs at the late stage of the Earths accretion rather than at the early stage of it. The GBI model will pose a great challenge to many problems between the equilibrium of Earths core and mantle. It provides an additional source for the excess of highly siderophile elements in the Earths mantle and also brings excessive lithophile elements to the Earths core. The GBI model may shed light on the study of Moon-formation and chemical differentiations of the pro-Earth.

Unlike the ordinary supernovae (SNe) some of which are hydrogen and helium deficient (called Type Ic SNe), broad-lined Type Ic SNe (SNe Ic-bl) are very energetic events, and all SNe coincident with bona fide long duration gamma-ray bursts (LGRBs) are of Type Ic-bl. Understanding the progenitors and the mechanism driving SN Ic-bl explosions vs those of their SNe Ic cousins is key to understanding the SN-GRB relationship and jet production in massive stars. Here we present the largest set of host-galaxy spectra of 28 SNe Ic and 14 SN Ic-bl, all discovered before 2013 by the same untargeted survey, namely the Palomar Transient Factory (PTF). We carefully measure their gas-phase metallicities, stellar masses (M*s) and star-formation rates (SFRs) by taking into account recent progress in the metallicity field and propagating uncertainties correctly. We further re-analyze the hosts of 10 literature SN-GRBs using the same methods and compare them to our PTF SN hosts with the goal of constraining their progenitors from their local environments by conducting a thorough statistical comparison, including upper limits. We find that the metallicities, SFRs and M*s of our PTF SN Ic-bl hosts are statistically comparable to those of SN-GRBs, but significantly lower than those of the PTF SNe Ic. The mass-metallicity relations as defined by the SNe Ic-bl and SN-GRBs are not significantly different from the same relations as defined by the SDSS galaxies, in contrast to claims by earlier works. Our findings point towards low metallicity as a crucial ingredient for SN Ic-bl and SN-GRB production since we are able to break the degeneracy between high SFR and low metallicity. We suggest that the PTF SNe Ic-bl may have produced jets that were choked inside the star or were able break out of the star as unseen low-luminosity or off-axis GRBs.

Radial velocity is one of key measurements in understanding the fundamental properties of stars, stellar clusters and the Galaxy. A plate of stars in the Kepler field were observed in May of 2018 with the medium-resolution spectrographs of LAMOST, aiming to test the performance of this new system which is the upgraded equipment of LAMOST after the first five-year regular survey.We present our analysis on the radial velocity measurements (RVs) derived from these data. The results show that slight and significant systematic errors exist among the RVs obtained from the spectra collected by different spectrographs and exposures, respectively. After correcting the systematic errors with different techniques, the precision of RVs reaches ~1.3, ~1.0, ~0.5 and ~0.3 km/s at S/Nr = 10, 20, 50, and 100, respectively. Comparing with the RVs of the standard stars of the APOGEE survey, our RVs are calibrated with a zero-point shift of ~7 km/s. The results indicate that the LAMOST medium-resolution spectroscopic system may provide RVs in a reasonable accuracy and precision for the selected targets.

Quasi-periodic fast-propagating (QFP) magnetosonic waves and extreme ultraviolet (EUV) waves were proposed to be driven by solar flares and coronal mass ejections (CMEs), respectively. In this Letter, we present a detailed analysis of an interesting event in which we find that both QFP magnetosonic waves and EUV waves are excited simultaneously in one solar eruption event. The co-existence of the two wave phenomena offers an excellent opportunity to explore their driving mechanisms. The QFP waves propagate in a funnel-like loop system with a speed of 682--837 \speed{} and a lifetime of 2 minutes. On the contrary, the EUV waves, which present a faster component and a slower component, propagate in a wide angular extent, experiencing reflection and refraction across a magnetic quasi-separatrix layer. The faster component of the EUV waves travels with a speed of 412--1287 \speed{}, whereas the slower component travels with a speed of 246--390 \speed{}. The lifetime of the EUV waves is $\sim$15 minutes. It is revealed that the faster component of the EUV waves is cospatial with the first wavefront of the QFP wave train. Besides, The QFP waves have a period of about 45$\pm$5 seconds, which is absent in the associated flares. All these results imply that QFP waves can also be excited by mass ejections, including CMEs or jets.

In this paper, we present a detailed analysis of a coronal blowout jet eruption which was associated with an obvious extreme-ultraviolet (EUV) wave and one complicated coronal mass ejection (CME) event based on the multi-wavelength and multi-view-angle observations from {\sl Solar Dynamics Observatory} and {\sl Solar Terrestrial Relations Observatory}. It is found that the triggering of the blowout jet was due to the emergence and cancellation of magnetic fluxes on the photosphere. During the rising stage of the jet, the EUV wave appeared just ahead of the jet top, lasting about 4 minutes and at a speed of 458 - \speed{762}. In addition, obvious dark material is observed along the EUV jet body, which confirms the observation of a mini-filament eruption at the jet base in the chromosphere. Interestingly, two distinct but overlapped CME structures can be observed in corona together with the eruption of the blowout jet. One is in narrow jet-shape, while the other one is in bubble-shape. The jet-shaped component was unambiguously related with the outwardly running jet itself, while the bubble-like one might either be produced due to the reconstruction of the high coronal fields or by the internal reconnection during the mini-filament ejection according to the double-CME blowout jet model firstly proposed by Shen et al. (2012b), suggesting more observational evidence should be supplied to clear the current ambiguity based on large samples of blowout jets in future studies.

In this paper we describe the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to accretion flows in the strong field gravity regime around both stellar-mass and supermassive black-holes. eXTP has the unique capability of using advanced 'spectral-timing-polarimetry' techniques to analyze the rapid variations with three orthogonal diagnostics of the flow and its geometry, yielding unprecedented insight into the inner accreting regions, the effects of strong field gravity on the material within them and the powerful outflows which are driven by the accretion process.

A powerful coronal mass ejection (CME) occurred on 2017 September 10 near the end of the declining phase of the historically weak solar cycle 24. We obtain new insights concerning the geometry and kinematics of CME-driven shocks in relation to their heliospheric impacts from the optimal, multi-spacecraft observations of the eruption. The shock, which together with the CME driver can be tracked from the early stage to the outer corona, shows a large oblate structure produced by the vast expansion of the ejecta. The expansion speeds of the shock along the radial and lateral directions are much larger than the translational speed of the shock center, all of which increase during the flare rise phase, peak slightly after the flare maximum and then decrease. The near simultaneous arrival of the CME-driven shock at the Earth and Mars, which are separated by 156.6$^{\circ}$ in longitude, is consistent with the dominance of expansion over translation observed near the Sun. The shock decayed and failed to reach STEREO A around the backward direction. Comparison between ENLIL MHD simulations and the multi-point in situ measurements indicates that the shock expansion near the Sun is crucial for determining the arrival or non-arrival and space weather impact at certain heliospheric locations. The large shock geometry and kinematics have to be taken into account and properly treated for accurate predictions of the arrival time and space weather impact of CMEs.

The recent detections of gravitational waves (GWs) reported by LIGO/Virgo collaborations have made significant impact on physics and astronomy. A global network of GW detectors will play a key role to solve the unknown nature of the sources in coordinated observations with astronomical telescopes and detectors. Here we introduce KAGRA (former name LCGT; Large-scale Cryogenic Gravitational wave Telescope), a new GW detector with two 3-km baseline arms arranged in the shape of an "L", located inside the Mt. Ikenoyama, Kamioka, Gifu, Japan. KAGRA's design is similar to those of the second generations such as Advanced LIGO/Virgo, but it will be operating at the cryogenic temperature with sapphire mirrors. This low temperature feature is advantageous for improving the sensitivity around 100 Hz and is considered as an important feature for the third generation GW detector concept (e.g. Einstein Telescope of Europe or Cosmic Explorer of USA). Hence, KAGRA is often called as a 2.5 generation GW detector based on laser interferometry. The installation and commissioning of KAGRA is underway and its cryogenic systems have been successfully tested in May, 2018. KAGRA's first observation run is scheduled in late 2019, aiming to join the third observation run (O3) of the advanced LIGO/Virgo network. In this work, we describe a brief history of KAGRA and highlights of main feature. We also discuss the prospects of GW observation with KAGRA in the era of O3. When operating along with the existing GW detectors, KAGRA will be helpful to locate a GW source more accurately and to determine the source parameters with higher precision, providing information for follow-up observations of a GW trigger candidate.

Gap-like structures in protoplanetary disks are likely related to planet formation processes. In this paper, we present and analyze high resolution (0.17*0.11 arcsec) 1.3 mm ALMA continuum observations of the protoplanetary disk around the Herbig Ae star MWC 480. Our observations for the first time show a gap centered at ~74au with a width of ~23au, surrounded by a bright ring centered at ~98au from the central star. Detailed radiative transfer modeling of both the ALMA image and the broadband spectral energy distribution is used to constrain the surface density profile and structural parameters of the disk. If the width of the gap corresponds to 4~8 times the Hill radius of a single forming planet, then the putative planet would have a mass of 0.4~3 M_Jup. We test this prediction by performing global three-dimensional smoothed particle hydrodynamic gas/dust simulations of disks hosting a migrating and accreting planet. We find that the dust emission across the disk is consistent with the presence of an embedded planet with a mass of ~2.3 M_Jup at an orbital radius of ~78au. Given the surface density of the best-fit radiative transfer model, the amount of depleted mass in the gap is higher than the mass of the putative planet, which satisfies the basic condition for the formation of such a planet.

We present the Data Release 4&5 quasar catalog from the quasar survey by Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST), which includes quasars observed between September 2015 and June 2017. There are a total of 19,253 quasars identified by visual inspections of the spectra. Among them, 11,458 are independently discovered by LAMOST, in which 3296 were reported by SDSS DR12 and DR14 quasar catalog after our survey began, while the rest 8162 are new discoveries of LAMOST. We provide the emission line measurements for the Halpha, Hbeta, MgII and/or CIV for 18100 quasars. Since LAMOST does not have absolute flux calibration information, we obtain the monochromatic continuum luminosities by fitting the SDSS photometric data using the quasar spectra, and then estimate the black hole masses. The catalog and spectra for these quasars are available online. This is the third installment in the series of LAMOST quasar survey which has released spectra for totally ~43,000 quasars hitherto. There are 24,772 independently discovered quasars, 17,128 of which are newly discovered. In addition to the great supplement to the new quasar discoveries, LAMOST has also provided a large database (overlapped with SDSS) for investigating the quasar spectral variability and discovering unusual quasars, including changing-look quasars, with ongoing and upcoming large surveys.

We study the magnetic field evolution in the active region (AR) 12673 that produced the largest solar flare in the past decade on 2017 September 6. Fast flux emergence is one of the most prominent features of this AR. We calculate the magnetic helicity from photospheric tangential flows that shear and braid field lines (shear-helicity), and from normal flows that advect twisted magnetic flux into the corona (emergence-helicity), respectively. Our results show that the emergence-helicity accumulated in the corona is $-1.6\times10^{43}~Mx^2$ before the major eruption, while the shear-helicity accumulated in the corona is $-6\times10^{43}~Mx^2$, which contributes about 79\% of the total helicity. The shear-helicity flux is dominant throughout the overall investigated emergence phase. Our results imply that the emerged fields initially contain relatively low helicity. Much more helicity is built up by shearing and converging flows acting on preexisted and emerging flux. Shearing motions are getting stronger with the flux emergence, and especially on both sides of the polarity inversion line of the core field region. The evolution of the vertical currents shows that most of the intense currents do not appear initially with the emergence of the flux, which implies that most of the emerging flux is probably not strongly current-carrying. The helical magnetic fields (flux rope) in the core field region are probably formed by long-term photospheric motions. The shearing and converging motions are continuously generated driven by the flux emergence. AR 12673 is a representative as photospheric motions contribute most of the nonpotentiality in the AR with vigorous flux emergence.

We study the magnetic field evolution in the active region (AR) 12673 that produced the largest solar flare in the past decade on 2017 September 6. Fast flux emergence is one of the most prominent features of this AR. We calculate the magnetic helicity from photospheric tangential flows that shear and braid field lines (shear-helicity), and from normal flows that advect twisted magnetic flux into the corona (emergence-helicity), respectively. Our results show that the emergence-helicity accumulated in the corona is $-1.6\times10^{43}~Mx^2$ before the major eruption, while the shear-helicity accumulated in the corona is $-6\times10^{43}~Mx^2$, which contributes about 79\% of the total helicity. The shear-helicity flux is dominant throughout the overall investigated emergence phase. Our results imply that the emerged fields initially contain relatively low helicity. Much more helicity is built up by shearing and converging flows acting on preexisted and emerging flux. Shearing motions are getting stronger with the flux emergence, and especially on both sides of the polarity inversion line of the core field region. The evolution of the vertical currents shows that most of the intense currents do not appear initially with the emergence of the flux, which implies that most of the emerging flux is probably not strongly current-carrying. The helical magnetic fields (flux rope) in the core field region are probably formed by long-term photospheric motions. The shearing and converging motions are continuously generated driven by the flux emergence. AR 12673 is a representative as photospheric motions contribute most of the nonpotentiality in the AR with vigorous flux emergence.

We examine a method to detect the light speed variation from gamma ray burst data observed by the Fermi Gamma-ray Space Telescope (FGST). We suggest new criteria to determine the characteristic time for low energy photons by the energy curve and the average energy curve, and obtain similar results compared with those from the light curve. We offer a new criterion with both the light curve and the average energy curve to determine the characteristic time for low energy photons. We then apply the new criteria to the GBM NaI data, the GBM BGO data, and the LAT LLE data, and obtain consistent results for three different sets of low energy photons from different FERMI detectors.

Timing analysis of PSR J1705$-$1906 using data from Nanshan 25-m and Parkes 64-m radio telescopes, which span over fourteen years, shows that the pulsar exhibits significant proper motion, and rotation instability. We updated the astrometry parameters and the spin parameters of the pulsar. In order to minimize the effect of timing irregularities on measuring its position, we employ the Cholesky method to analyse the timing noise. We obtain the proper motion of $-$77(3) \,mas\,yr$^{-1}$ in right ascension and $-$38(29) \,mas\,yr$^{-1}$ in declination. The power spectrum of timing noise is analyzed for the first time, which gives the spectral exponent $\alpha=-5.2$ for the power-law model indicating that the fluctuations in spin frequency and spin-down rate dominate the red noise. We detect two small glitches from this pulsar with fractional jump in spin frequency of $\Delta \nu/\nu\sim2.9\times10^{-10}$ around MJD~55199 and $\Delta \nu/\nu\sim2.7\times10^{-10}$ around MJD~55953. Investigations of pulse profile at different time segments suggest no significant changes in the pulse profiles around the two glitches.

A fundamental property of molecular clouds is that they are turbulent, but how this turbulence is generated and maintained is unknown. One possibility is that stars forming within the cloud regenerate turbulence via their outflows, winds and radiation ("feedback"). Disentangling motions created by feedback from the initial cloud turbulence is challenging, however. Here we confront the relationship between stellar feedback and turbulence by identifying and separating the local and global impact of stellar winds. We analyze magnetohydrodyanamic simulations in which we track wind material as it interacts with the ambient cloud. By distinguishing between launched material, gas entrained by the wind and pristine gas we show energy is transferred away from the sources via magnetic waves excited by the expanding wind shells. This action at a distance enhances the fraction of stirring motion compared to compressing motion and produces a flatter velocity power spectrum. We conclude stellar feedback accounts for significant energy transfer within molecular clouds, an impact enhanced by magnetic waves, which have previously been neglected by observations. Altogether, stellar feedback can partially offset global turbulence dissipation.

Rings are the most frequently revealed substructure in ALMA dust observations of protoplanetary disks, but their origin is still hotly debated. In this paper, we identify dust substructures in 12 disks and measure their properties to investigate how they form. This subsample of disks is selected from a high-resolution ($\sim0.12''$) ALMA 1.33 mm survey of 32 disks in the Taurus star-forming region, which was designed to cover a wide range of sub-mm brightness and to be unbiased to previously known substructures. While axisymmetric rings and gaps are common within our sample, spiral patterns and high contrast azimuthal asymmetries are not detected. Fits of disk models to the visibilities lead to estimates of the location and shape of gaps and rings, the flux in each disk component, and the size of the disk. The dust substructures occur across a wide range of stellar mass and disk brightness. Disks with multiple rings tend to be more massive and more extended. The correlation between gap locations and widths, the intensity contrast between rings and gaps, and the separations of rings and gaps could all be explained if most gaps are opened by low-mass planets (super-Earths and Neptunes) in the condition of low disk turbulence ($\alpha=10^{-4}$). The gap locations are not well correlated with the expected locations of CO and N$_2$ ice lines, so condensation fronts are unlikely to be a universal mechanism to create gaps and rings, though they may play a role in some cases.

We use multi-band imagery data from the Sloan Digital Sky Survey (SDSS) to measure projected distances of 302 supernova type Ia (SNIa) from the centre of their host galaxies, normalized to the galaxy's brightness scale length, with a Bayesian approach. We test the hypothesis that SNIas further away from the centre of their host galaxy are less subject to dust contamination (as the dust column density in their environment is smaller) and/or come from a more homogeneous environment. Using the Mann-Whitney U test, we find a statistically significant difference in the observed colour correction distribution between SNIas that are near and those that are far from the centre of their host. The local p-value is 3 x 10^{-3}, which is significant at the 5 per cent level after look-elsewhere effect correction. We estimate the residual scatter of the two subgroups to be 0.073 +/- 0.018 for the far SNIas, compared to 0.114 +/- 0.009 for the near SNIas -- an improvement of 30 per cent, albeit with a low statistical significance of 2sigma. This confirms the importance of host galaxy properties in correctly interpreting SNIa observations for cosmological inference.

Ongoing and future space missions aim to identify potentially habitable planets in our Solar System and beyond. Planetary habitability is determined not only by a planet's current stellar insolation and atmospheric properties, but also by the evolutionary history of its climate. It has been suggested that icy planets and moons become habitable after their initial ice shield melts as their host stars brighten. Here we show from global climate model simulations that a habitable state is not achieved in the climatic evolution of those icy planets and moons that possess an inactive carbonate-silicate cycle and low concentrations of greenhouse gases. Examples for such planetary bodies are the icy moons Europa and Enceladus, and certain icy exoplanets orbiting G and F stars. We find that the stellar fluxes that are required to overcome a planet's initial snowball state are so large that they lead to significant water loss and preclude a habitable planet. Specifically, they exceed the moist greenhouse limit, at which water vapour accumulates at high altitudes where it can readily escape, or the runaway greenhouse limit, at which the strength of the greenhouse increases until the oceans boil away. We suggest that some icy planetary bodies may transition directly to a moist or runaway greenhouse without passing through a habitable Earth-like state.

We study scalar condensations around asymptotically Anti-de Sitter(AdS) regular reflecting shells. We show that the charged scalar field can condense around charged reflecting shells with both Dirichlet and Neumann boundary conditions. In particular, the radii of the asymptotically AdS hairy shells are discrete, which is similar to cases in asymptotically flat spacetimes. We also provide upper bounds for the radii of the hairy Dirichlet reflecting shells and above the bound, the scalar field cannot condense around the shell.

We demonstrate a compact 29.3 GHz visible astro-comb covering the spectrum from 560nm to 700nm. A 837 MHz Yb:fiber laser frequency comb phase locked to a Rb clock served as the seed comb to ensure the frequency stability and high side mode suppression ratio. After the visible super-continuum generation, a cavity-length-fixed Fabry-Perot cavity made by ultra-low expansion glass was utilized to filter the comb teeth for eliminating the rapid active dithering. The mirrors were home-made complementary chirped mirrors pair with zero dispersion and high reflection to guarantee no mode skipping. These filtered comb teeth were clearly resolved in an astronomical spectrograph of 49,000 resolution, exhibiting sharp linetype, zero noise floor, and uniform exposure amplitude.

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

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

The solar eruption on 2012 January 27 resulted in a wide-spread solar energetic particle (SEP) event observed by STEREO A and the near-Earth spacecraft (separated by 108\degree). The event was accompanied by an X-class flare, extreme-ultraviolet (EUV) wave and fast coronal mass ejection (CME). We investigate the particle release by comparing the release times of particles at the spacecraft and the times when magnetic connectivity between the source and the spacecraft was established. The EUV wave propagating to the magnetic footpoint of the spacecraft in the lower corona and the shock expanding to the open field line connecting the spacecraft in the upper corona are thought to be responsible for the particle release. We track the evolution of the EUV wave and model the propagation of the shock using EUV and white-light observations. No obvious evidence indicates that the EUV wave reached the magnetic footpoint of either STEREO A or L1-observers. Our shock modeling shows that the release time of the particles observed at L1 was consistent with the time when the shock first established contact with the magnetic field line connecting L1-observers. The release of the particles observed by STEREO A was delayed relative to the time when the shock was initially connected to STEREO A via the magnetic field line. We suggest that the particle acceleration efficiency of the portion of the shock connected to the spacecraft determines the release of energetic particles at the spacecraft.

Using the 135-second cadence of the photospheric vector data provided by the Helioseismic and Magnetic Imager telescope on board the Solar Dynamic Observatory, we examined the time-evolution of magnetic helicity fluxes across the photosphere during 16 flares with the energy class lower than M5.0. During the flare in 4 out of 16 events, we found impulsive changes in the helicity fluxes. This indicates that even the flare with less energy could be associated with anomalistic transportation of the magnetic helicity across the photosphere. Accompanying the impulsive helicity fluxes, the poynting fluxes across the photosphere evolved from positive to negative. As such, the transportations of magnetic energy across the photosphere were toward solar interior during these flares. In each of the 4 events, the impulsive change in the helicity flux was always mainly contributed by abrupt change in horizontal velocity field on a sunspot located near the flaring polarity inversion line. The velocity field on each sunspot shows either an obvious vortex patten or an shearing patten relative to the another magnetic polarity, which tended to relax the magnetic twist or shear in the corona. During these flares, abrupt change in the Lorentz force acting on these sunspots were found. The rotational motions and shearing motions of these sunspots always had the same directions with the resultant Lorentz forces. These results support the view that the impulsive helicity transportation during the flare could be driven by the change in the Lorentz force applied on the photosphere.

Using the 135-second cadence of the photospheric vector data provided by the Helioseismic and Magnetic Imager telescope on board the Solar Dynamic Observatory, we examined the time-evolution of magnetic helicity fluxes across the photosphere during 16 flares with the energy class lower than M5.0. During the flare in 4 out of 16 events, we found impulsive changes in the helicity fluxes. This indicates that even the flare with less energy could be associated with anomalistic transportation of the magnetic helicity across the photosphere. Accompanying the impulsive helicity fluxes, the poynting fluxes across the photosphere evolved from positive to negative. As such, the transportations of magnetic energy across the photosphere were toward solar interior during these flares. In each of the 4 events, the impulsive change in the helicity flux was always mainly contributed by abrupt change in horizontal velocity field on a sunspot located near the flaring polarity inversion line. The velocity field on each sunspot shows either an obvious vortex patten or an shearing patten relative to the another magnetic polarity, which tended to relax the magnetic twist or shear in the corona. During these flares, abrupt change in the Lorentz force acting on these sunspots were found. The rotational motions and shearing motions of these sunspots always had the same directions with the resultant Lorentz forces. These results support the view that the impulsive helicity transportation during the flare could be driven by the change in the Lorentz force applied on the photosphere.

We present the observations of an extreme ultraviolet (EUV) wave, a quasi-periodic fast-propagating (QFP) magnetosonic wave, and a kink wave that were simultaneously associated with the impingement of a coronal jet upon a group of coronal loops. After the interaction, the coronal loop showed obvious kink oscillation that had a period of about 428 seconds. In the meantime, a large-scale EUV wave and a QFP wave are observed on the west of the interaction position. It is interesting that the QFP wave showed refraction effect during its passing through two strong magnetic regions. The angular extent, speed, and lifetime of the EUV (QFP) wave were about 140 (40) degree, 423 (322) km/s, and 6 (26) minutes, respectively. It is measured that the period of the QFP wave was about 390 +/- 100 second. Based on the observational analysis results, we propose that the kink wave was probably excited by the interaction of the jet; the EUV was probably launched by the sudden expansion of the loop system due to the impingement of the coronal jet; and the QFP wave was possibly formed through the dispersive evolution of the disturbance caused by the jet-loop interaction.

We report the observations of a moderate but relatively intense geo-effective solar eruption on 2015 November 4 from the peripheral diffusive polarities of active region 12443. We use space-borne Solar Dynamics Observatory and ACE observations. EUV images identified helical pattern along a filament channel and we regard this channel as flux-rope structure. Flow velocity derived from tracked magnetograms infers converging motion along the polarity inversion line beneath the filament channel. An associated magnetic cancellation process was detected in the converging region. Further, the pre-eruptive EUV brightening was observed in the converging region, the most intense part of which appeared in the magnetic cancellation region. These observations imply that the converging and cancelling flux probably contributed to the formation of the helical magnetic fields associated with the flux rope. A filament-height estimation method suggests that the middle part of the filament probably lies at a low altitude and was consistent with the initial place of the eruption. A thick current channel associated with the flux rope is also determined. For a expanding thick current channel, the critical height of the decay index for torus instability lies in the range of 37 - 47 Mm. Southward magnetic fields in the sheath and the ejecta induced a geomagnetic storm with a Dst global minimum of ~-90 nT.

Well-sampled optical light curves of 50 gamma-ray bursts (GRBs) with plateau features are compiled from the literature. By empirical fitting, we obtained the parameters of the optical plateaus, such as the decay slopes ($\alpha_{\rm 1}$ and $\alpha_{\rm 2}$), the break times ($T_{\rm b}$), and the corresponding optical fluxes ($F_{\rm b}$) at the break times. The break time of optical plateaus ranges from tens of seconds to $10^6$ seconds, with a typical value about $10^4$ seconds. We have calculated the break luminosity, and it mainly ranges from $10^{44}$ erg $s^{-1}$ to $10^{47}$ erg $s^{-1}$, which is generally two or three orders of magnitude less than the corresponding break luminosity of the X-ray afterglow plateaus. We reanalyzed the optical plateaus and also found that a significantly tighter correlation exists when we added the isotropic equivalent energy of GRBs $E_{\rm \gamma,iso}$ into the $L_{\rm b,z}-T_{\rm b,z}$ relation. The best fit correlation is obtained to be $L_{\rm b,z}\propto T_{\rm b,z}^{-0.9}E_{\rm \gamma,iso}^{0.4}$. We next explored the possible correlations among $L_{\rm b,z}$, $T_{\rm b,z}$ and $E_{\rm p,i}$, and found there is also a tight correlation between them, which takes the form of $L_{\rm b,z}\propto T_{\rm b,z}^{-0.9}E_{\rm p,i}^{0.5}$. We argue that these two tight $L_{\rm b,z}-T_{\rm b,z}-E_{\rm \gamma,iso}$ and $L_{\rm b,z}-T_{\rm b,z}-E_{\rm p,i}$ correlations are more physical, and it may be directly related to radiation physics of GRBs. The tight correlations are possible to be used as standard candles.

We present two events to study the driving mechanism of extreme-ultraviolet (EUV) waves that are not associated with coronal mass ejections (CMEs), by using high resolution observations taken by the Atmospheric Imaging Assembly (AIA) on board Solar Dynamics Observatory. Observational results indicate that the observed EUV waves were accompanied by ares and coronal jets, but without CMEs that were regarded as the driver of most EUV waves in previous studies. In the first case, it is observed that a coronal jet ejected along a transequatorial loop system at a plane-of-the-sky (POS) speed of $335 \pm 22$ km/s, in the meantime, an arc-shaped EUV wave appeared on the eastern side of the loop system. In addition, the EUV wave further interacted with another interconnecting loop system and launched a fast propagating (QFP) magnetosonic wave along the loop system, which had a period of 200 s and a speed of $388 \pm 65$ km/s, respectively. In the second case, we also observed a coronal jet ejected at a POS speed of $282 \pm 44$ km/s along a transequatorial loop system and the generation of bright EUV wave on the eastern side of the loop system. Based on the observational results, we propose that the observed EUV waves on the eastern side of the transequatorial loop systems are fast-mode magnetosonic waves, and they were driven by the sudden lateral expansion of the transequatorial loop systems due to the direct impingement of the associated coronal jets, while the QFP wave in the first case formed due to the dispersive evolution of the disturbance caused by the interaction between the EUV wave and the interconnecting coronal loops. It is noted that EUV waves driven by sudden loop expansions have shorter lifetimes than those driven by CMEs.

The rare case of changing-look (CL) AGNs, with the appearance or disappearance of broad Balmer emission lines within a few years, challenges our understanding of the AGN unified model. We present a sample of 21 new CL AGNs at $0.08<z<0.58$, which doubles the number of such objects known to date. These new CL AGNs were discovered by several ways, from (1) repeat spectra in the SDSS, (2) repeat spectra in the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) and SDSS, and (3) photometric variability and new spectroscopic observations. We use the photometric data from surveys, including the SDSS imaging survey, the Pan-STARRS1, the DESI Legacy imaging survey, the Wide-field Infrared Survey Explorer (WISE), the Catalina Real-time Transient Survey, and the Palomar Transient Factory. The estimated upper limits of transition timescale of the CL AGNs in this sample spans from 0.9 to 13 years in the rest frame. The continuum flux in the optical and mid-infrared becomes brighter when the CL AGNs turn on, or vice versa. Variations of more than 0.2 mag in $W1$ band were detected in 15 CL AGNs during the transition. The optical and mid-infrared variability is not consistent with the scenario of variable obscuration in 10 CL AGNs at more than $3\sigma$ confidence level. We confirm a bluer-when-brighter trend in the optical. However, the mid-infrared WISE colors $W1-W2$ become redder when the objects become brighter in the $W1$ band, possibly due to a stronger hot dust contribution in the $W2$ band when the AGN activity becomes stronger. The physical mechanism of type transition is important for understanding the evolution of AGNs.

We demonstrate the dragging of the magnetic field by the super-Alfvenic shear flows out of the reconnection plane can strongly localize the reconnection x-line in collisionless plasmas, reversing the current direction at the x-line. Reconnection with this new morphology, which is impossible in resistive-magnetohydrodynamic (MHD), is enabled by electron inertia. Surprisingly, the quasi-steady reconnection rate remains of order 0.1 even though the aspect ratio of the local x-line geometry is larger than unity. We explain this by examining the transport of the reconnected magnetic flux and the opening angle made by the upstream magnetic field, concluding that the reconnection rate is still limited by the constraint imposed at the inflow region. This study further suggests the nearly universal fast rate value of order 0.1 cannot be explained by the physics of tearing modes, nor can it be explained by a universal localization mechanism.

The detailed observational analysis of a homologous Extreme-ultraviolet (EUV) wave event is presented to study the driving mechanism and the physical property of the EUV waves, combining high resolution data taken by the Solar Dynamics Observatory and the Solar TErrestrial RElations Observatory. It is observed that four homologous EUV waves originated from the same active region AR11476 within about one hour, and the time separations between consecutive waves were of 8 - 20 minutes. The waves showed narrow arc-shaped wavefronts and propagated in the same direction along a large-scale transequatorial loop system at a speed of 648 - 712 km/s and a deceleration of 0.985 - 1.219 km/s2. The EUV waves were accompanied by weak flares, coronal jets, and radio type III bursts, in which the EUV waves were delayed with respect to the start times of the radio type III bursts and coronal jets about 2 - 13 and 4 - 9 minutes, respectively. Different to previous studies of homologous EUV waves, no coronal mass ejections were found in the present event. Based on the observational results and the close temporal the spatial relationship between the EUV waves and the coronal jets, for the first time, we propose that the observed homologous EUV waves were large-amplitude nonlinear fast-mode magnetosonic waves or shocks driven by the associated recurrent coronal jets, resemble the generation mechanism of a piston shock in a tube. In addition, it is found that the recurrent jets were tightly associated with the alternating flux cancellation and emergence in the eruption source region and radio type III bursts.

The detailed observational analysis of a homologous Extreme-ultraviolet (EUV) wave event is presented to study the driving mechanism and the physical property of the EUV waves, combining high resolution data taken by the Solar Dynamics Observatory and the Solar TErrestrial RElations Observatory. It is observed that four homologous EUV waves originated from the same active region AR11476 within about one hour, and the time separations between consecutive waves were of 8 - 20 minutes. The waves showed narrow arc-shaped wavefronts and propagated in the same direction along a large-scale transequatorial loop system at a speed of 648 - 712 km/s and a deceleration of 0.985 - 1.219 km/s2. The EUV waves were accompanied by weak flares, coronal jets, and radio type III bursts, in which the EUV waves were delayed with respect to the start times of the radio type III bursts and coronal jets about 2 - 13 and 4 - 9 minutes, respectively. Different to previous studies of homologous EUV waves, no coronal mass ejections were found in the present event. Based on the observational results and the close temporal the spatial relationship between the EUV waves and the coronal jets, for the first time, we propose that the observed homologous EUV waves were large-amplitude nonlinear fast-mode magnetosonic waves or shocks driven by the associated recurrent coronal jets, resemble the generation mechanism of a piston shock in a tube. In addition, it is found that the recurrent jets were tightly associated with the alternating flux cancellation and emergence in the eruption source region and radio type III bursts.

We present the observations of a double-decker filament to study its formation, triggering, and eruption physics. It is observed that the double-decker filament was formed by splitting of an original single filament. During the splitting process, intermittent bright point bursts are observed in the filament channel, which resulted in the generation of the upper filament branch. The eruption of the newly formed double-decker filament was possibly triggered by two recurrent two-sided loop jets in the filament channel and the continuous mass unloading from the upper filament body. The interaction between the first jet and the filament directly resulted in the unstable of the lower branch and the fast rising phase of the upper branch. The second jet occurred at the same site about three hours after the first one, which further disturbed and accelerated the rising of the lower filament branch. It is interesting that the rising lower branch overtook the upper one, and then the two branches probably merged into one filament. Finally, the whole filament erupted violently and caused a large-scale coronal mass ejection, leaving behind a pair of flare ribbons and two dimming regions on the both sides of the filament channel. We think that the intermittent bursts may directly result in the rearrangement of the filament magnetic field and therefore the formation of the double-decker filament, then the recurrent jets further caused the fully eruption of the entire filament system. The study provides convincing evidence for supporting the scenario that a double-decker filament can be formed by splitting a single filament into two branches.

We investigate the kinetic properties of a typical fast-mode shock inside an interplanetary coronal mass ejection (ICME) observed on 1998 August 6 at 1 AU, including particle distributions and wave analysis with the in situ measurements from Wind. Key results are obtained concerning the shock and the shock-ICME interaction at kinetic scales: (1) gyrating ions, which may provide energy dissipation at the shock in addition to wave-particle interactions, are observed around the shock ramp; (2) despite the enhanced proton temperature anisotropy of the shocked plasma, the low plasma beta inside the ICME constrains the shocked plasma under the thresholds of the ion cyclotron and mirror-mode instabilities; (3) whistler heat flux instabilities, which can pitch--angle scatter halo electrons through a cyclotron resonance, are observed around the shock, and can explain the disappearance of bidirectional electrons inside the ICME together with normal betatron acceleration; (4) whistler waves near the shock are likely associated with the whistler heat flux instabilities excited at the shock ramp, which is consistent with the result that the waves may originate from the shock ramp; (5) the whistlers share a similar characteristic with the shocklet whistlers observed by Wilson et al, providing possible evidence that the shock is decaying because of the strong magnetic field inside the ICME.

The orientation and stability of the reconnection x-line in asymmetric geometry is studied using three-dimensional (3D) particle-in-cell simulations. We initiate reconnection at the center of a large simulation domain to minimize the boundary effect. The resulting x-line has sufficient freedom to develop along an optimal orientation, and it remains laminar. Companion 2D simulations indicate that this x-line orientation maximizes the reconnection rate. The divergence of the non-gyrotropic pressure tensor breaks the frozen-in condition, consistent with its 2D counterpart. We then design 3D simulations with one dimension being short to fix the x-line orientation, but long enough to allow the growth of the fastest growing oblique tearing modes. This numerical experiment suggests that reconnection tends to radiate secondary oblique tearing modes if it is externally (globally) forced to proceed along an orientation not favored by the local physics. The development of oblique structure easily leads to turbulence inside small periodic systems.

We analyze the polarization content of gravitational waves in Horndeski theory. Besides the familiar plus and cross polarizations in Einstein's General Relativity, there is one more polarization state which is the mixture of the transverse breathing and longitudinal polarizations.The additional mode is excited by the massive scalar field. In the massless limit, the longitudinal polarization disappears, while the breathing one persists. The upper bound on the graviton mass severely constrains the amplitude of the longitudinal polarization, which makes its detection highly unlikely by the ground-based or space-borne interferometers in the near future. However, pulsar timing arrays might be able to detect the polarization excited by the massive scalar field. Since additional polarization states appear in alternative theories of gravity, the measurement of the polarizations of gravitational waves can be used to probe the nature of gravity. In addition to the plus and cross states, the detection of the breathing polarization means that gravitation is mediated by massless spin 2 and spin 0 fields, and the detection of both the breathing and longitudinal states means that gravitation is propagated by the massless spin 2 and massive spin 0 fields.

We investigate how a weak coronal mass ejection (CME) launched on 2016 October 8 without obvious signatures in the low corona produced a relatively intense geomagnetic storm. Remote sensing observations from SDO, STEREO and SOHO and in situ measurements from WIND are employed to track the CME from the Sun to the Earth. Using a graduated cylindrical shell (GCS) model, we estimate the propagation direction and the morphology of the CME near the Sun. CME kinematics are determined from the wide-angle imaging observations of STEREO A and are used to predict the CME arrival time and speed at the Earth. We compare ENLIL MHD simulation results with in situ measurements to illustrate the background solar wind where the CME was propagating. We also apply a Grad--Shafranov technique to reconstruct the flux rope structure from in situ measurements in order to understand the geo-effectiveness associated with the CME magnetic field structure. Key results are obtained concerning how a weak CME can generate a relatively intense geomagnetic storm: (1) there were coronal holes at low latitudes, which could produce high speed streams (HSSs) to interact with the CME in interplanetary space; (2) the CME was bracketed between a slow wind ahead and a HSS behind, which enhanced the southward magnetic field inside the CME and gave rise to the unexpected geomagnetic storm.

We present analysis of the magnetic field in seven solar flare regions accompanied by the pulsations of hard X-ray (HXR) emission. These flares were studied by Kuznetsov et al. (2016) (Paper~I), and chosen here because of the availability of the vector magnetograms for their parent active regions (ARs) obtained with the SDO/HMI data. In Paper~I, based on the observations only, it was suggested that a magnetic flux rope (MFR) might play an important role in the process of generation of the HXR pulsations. The goal of the present paper is to test this hypothesis by using the extrapolation of magnetic field with the non-linear force-free field (NLFFF) method. Having done this, we found that before each flare indeed there was an MFR elongated along and above a magnetic polarity inversion line (MPIL) on the photosphere. In two flare regions the sources of the HXR pulsations were located at the footpoints of different magnetic field lines wrapping around the central axis, and constituting an MFR by themselves. In five other flares the parent field lines of the HXR pulsations were not a part of an MFR, but surrounded it in the form of an arcade of magnetic loops. These results show that, at least in the analyzed cases, the "single flare loop" models do not satisfy the observations and magnetic field modeling, while are consistent with the concept that the HXR pulsations are a consequence of successive episodes of energy release and electron acceleration in different magnetic flux tubes (loops) of a complex AR. An MFR could generate HXR pulsations by triggering episodes of magnetic reconnection in different loops in the course of its non-uniform evolution along an MPIL. However, since three events studied here were confined flares, actual eruptions may not be required to trigger sequential particle acceleration episodes in the magnetic systems containing an MFR.

The observational analysis is performed to study the excitation mechanism and the propagation properties of a quasi-periodic fast-propagating (QFP) magnetosonic wave. The QFP wave was associated with the eruption of a nearby mini-filament and a small B4 GOES flare, which may indicate that the generation of a QFP wave do not need too much flare energy. The propagation of the QFP wave was along a bundle of funnel-shaped open loops with a speed of about 1100+/-78, and an acceleration of -2.2+/-1.1. Periodicity analysis indicates that the periods of the QFP wave are 43+/-6, 79+/-18 second. For the first time, we find that the periods of the QFP wave and the accompanying flare are inconsistent, which is different from the findings as reported in previous studies. We propose that the present QFP wave was possibly caused by the mechanism of dispersive evolution of an initially broadband disturbance resulted from the nearby mini-filament eruption.

As an essential part of China’s Gravitational Waves Program, the Ali CMB Polarization Telescope (AliCPT) is a ground-based experiment aiming at the Primordial Gravitational Waves (PGWs) by measuring B-mode polarization of Cosmic Microwave Background (CMB). First proposed in 2014 and currently in fast construction phase, AliCPT is China’s first CMB project that plans for commissioning in 2019. Led by the Institute of High Energy Physics (IHEP) under the Chinese Academy of Sciences (CAS), the project is a worldwide collaboration of more than fifteen universities and research institutes. Ali CMB Project is briefly introduced.

Coronal mass ejections (CMEs) are large-scale eruptions of plasma from the coronae of stars. Understanding the plasma processes involved in CME initiation has applications to space weather forecasting and laboratory plasma experiments. James et al. (Sol. Phys. 292, 71, 2017) used EUV observations to conclude that a magnetic flux rope formed in the solar corona above NOAA Active Region 11504 before it erupted on 14 June 2012 (SOL2012-06-14). In this work, we use data from the Solar Dynamics Observatory to model the coronal magnetic field of the active region one hour prior to eruption using a nonlinear force-free field extrapolation, and find a flux rope reaching a maximum height of 150 Mm above the photosphere. Estimations of the average twist of the strongly asymmetric extrapolated flux rope are between 1.35 and 1.88 turns, depending on the choice of axis, although the erupting structure was not observed to kink. The decay index near the apex of the axis of the extrapolated flux rope is comparable to typical critical values required for the onset of the torus instability, so we suggest that the torus instability drove the eruption.

In this paper, we will give a general introduction to the project of Ali CMB Polarization Telescope (AliCPT), which is a Sino-US joint project led by the Institute of High Energy Physics (IHEP) and has involved many different institutes in China. It is the first ground-based Cosmic Microwave Background (CMB) polarization experiment in China and an integral part of China's Gravitational Waves Program. The main scientific goal of AliCPT project is to probe the primordial gravitational waves (PGWs) originated from the very early Universe. The AliCPT project includes two stages. The first stage referred to as AliCPT-1, is to build a telescope in the Ali region of Tibet with an altitude of 5,250 meters. Once completed, it will be the worldwide highest ground-based CMB observatory and open a new window for probing PGWs in northern hemisphere. AliCPT-1 telescope is designed to have about 7,000 TES detectors at 90GHz and 150GHz. The second stage is to have a more sensitive telescope (AliCPT-2) with the number of detectors more than 20,000. Our simulations show that AliCPT will improve the current constraint on the tensor-to-scalar ratio $r$ by one order of magnitude with 3 years' observation. Besides the PGWs, the AliCPT will also enable a precise measurement on the CMB rotation angle and provide a precise test on the CPT symmetry. We show 3 years' observation will improve the current limit by two order of magnitude.

The sky brightness is a critical parameter for estimating the coronal observation conditions for solar observatory. As part of a site-survey project in Western China, we measured the sky brightness continuously at the Lijiang Observatory in Yunnan province in 2011. A sky brightness monitor (SBM) was adopted to measure the sky brightness in a region extending from 4.5 to 7.0 apparent solar radii based on the experience of the Daniel K. Inouye Solar Tele- scope (DKIST) site survey. Every month, the data were collected manually for at least one week. We collected statistics of the sky brightness at four bandpasses located at 450, 530, 890, and 940 nm. The results indicate that aerosol scattering is of great importance for the diurnal variation of the sky brightness. For most of the year, the sky brightness remains under 20 millionths per airmass before local Noon. On average, the sky brightness is less than 20 millionths accounts for 40.41 % of the total observing time in a clear day. The best observation time is from 9:00 to 13:00 (Beijing time). Lijiang Observatory is therefore suitable for coronagraphs investigating the structures and dynamics of the corona.

We present an investigation of the rotation and non-radial motion of a coronal mass ejection (CME) from AR 12468 on 2015 December 16 using observations from SDO, SOHO, STEREO A and Wind. The EUV and HMI observations of the source region show that the associated magnetic flux rope (MFR) axis pointed to the east before the eruption. We use a nonlinear fore-free field (NLFFF) extrapolation to determine the configuration of the coronal magnetic field and calculate the magnetic energy density distributions at different heights. The distribution of the magnetic energy density shows a strong gradient toward the northeast. The propagation direction of the CME from a Graduated Cylindrical Shell (GCS) modeling deviates from the radial direction of the source region by about 45 degr in longitude and about 30 degr in latitude, which is consistent with the gradient of the magnetic energy distribution around the AR.The MFR axis determined by the GCS modeling points southward, which has rotated counterclockwise by about 95 degr compared with the orientation of the MFR in the low corona.The MFR reconstructed by a Grad-Shafranov (GS) method at 1 AU has almost the same orientation as the MFR from the GCS modeling, which indicates that the MFR rotation occurred in the low corona. It is the rotation of the MFR that caused the intense geomagnetic storm with the minimum Dst of -155 nT. These results suggest that the coronal magnetic field surrounding the MFR plays a crucial role in the MFR rotation and propagation direction.

We present an investigation of the rotation and non-radial motion of a coronal mass ejection (CME) from AR 12468 on 2015 December 16 using observations from SDO, SOHO, STEREO A and Wind. The EUV and HMI observations of the source region show that the associated magnetic flux rope (MFR) axis pointed to the east before the eruption. We use a nonlinear fore-free field (NLFFF) extrapolation to determine the configuration of the coronal magnetic field and calculate the magnetic energy density distributions at different heights. The distribution of the magnetic energy density shows a strong gradient toward the northeast. The propagation direction of the CME from a Graduated Cylindrical Shell (GCS) modeling deviates from the radial direction of the source region by about 45 degr in longitude and about 30 degr in latitude, which is consistent with the gradient of the magnetic energy distribution around the AR.The MFR axis determined by the GCS modeling points southward, which has rotated counterclockwise by about 95 degr compared with the orientation of the MFR in the low corona.The MFR reconstructed by a Grad-Shafranov (GS) method at 1 AU has almost the same orientation as the MFR from the GCS modeling, which indicates that the MFR rotation occurred in the low corona. It is the rotation of the MFR that caused the intense geomagnetic storm with the minimum Dst of -155 nT. These results suggest that the coronal magnetic field surrounding the MFR plays a crucial role in the MFR rotation and propagation direction.

Using high temporal and high spatial resolution observations taken by the Atmospheric Imaging Assembly onboard the Solar Dynamics Observatory, we present the detailed observational analysis of a high quality quasi-periodic fast- propagating (QFP) magnetosonic wave that was associated with the eruption of a magnetic flux rope and a GOES C5.0 flare. For the first time, we find that the QFP wave lasted during the entire flare lifetime rather than only the rising phase of the accompanying flare as reported in previous studies. In addition, the propagation of the different parts of the wave train showed different kinematics and morphologies. For the southern (northern) part, the speed, duration, intensity variation are about 875 +/- 29 (1485 +/- 233) km/s, 45 (60) minutes, and 4% (2%), and the pronounced periods of them are 106 +/- 12 and 160 +/- 18 (75 +/- 10 and 120 +/- 16) seconds, respectively. It is interesting that the northern part of the wave train showed obvious refraction effect when they pass through a region of strong magnetic field. Periodicity analysis result indicates that all the periods of the QFP wave can be found in the period spectrum of the accompanying flare, suggesting their common physical origin. We propose that the quasi-periodic nonlinear magnetohydrodynamics process in the magnetic reconnection that produces the accompanying flare should be important for exciting of QFP wave, and the different magnetic distribution along different paths can account for the different speeds and morphology evolution of the wave fronts.

As a dedicated solar radio interferometer, the MingantU SpEctral RadioHeliograph (MUSER) generates massive observational data in the frequency range of 400 MHz -- 15 GHz. High-performance imaging forms a significantly important aspect of MUSER's massive data processing requirements. In this study, we implement a practical high-performance imaging pipeline for MUSER data processing. At first, the specifications of the MUSER are introduced and its imaging requirements are analyzed. Referring to the most commonly used radio astronomy software such as CASA and MIRIAD, we then implement a high-performance imaging pipeline based on the Graphics Processing Unit (GPU) technology with respect to the current operational status of the MUSER. A series of critical algorithms and their pseudo codes, i.e., detection of the solar disk and sky brightness, automatic centering of the solar disk and estimation of the number of iterations for clean algorithms, are proposed in detail. The preliminary experimental results indicate that the proposed imaging approach significantly increases the processing performance of MUSER and generates images with high-quality, which can meet the requirements of the MUSER data processing.

We investigate the equation of state(EoS) and the effect of the hadron-quark phase transition of strong interaction matter in compact stars. The hadron matter is described with the relativistic mean field theory,and the quark matter is described with the Dyson-Schwinger equation approach of QCD. The complete EoS of the hybrid star matter is constructed with not only the Gibbs construction but also the 3-window interpolation. The mass-radius relation of hybrid stars is also investigated. We find that, although the EoSs of both the hadron matter with hyperon and $\Delta$-baryon and the quark matter are generally softer than that of the nucleon matter, the 3-window interpolation construction may provide an EoS stiff enough for a hybrid star with mass exceeding 2$M_{\odot}^{}$ and, in turn, solve the so called "hyperon puzzle".

Major construction and initial-phase operation of a second-generation gravitational-wave detector KAGRA has been completed. The entire 3-km detector is installed underground in a mine in order to be isolated from background seismic vibrations on the surface. This allows us to achieve a good sensitivity at low frequencies and high stability of the detector. Bare-bones equipment for the interferometer operation has been installed and the first test run was accomplished in March and April of 2016 with a rather simple configuration. The initial configuration of KAGRA is named {\it iKAGRA}. In this paper, we summarize the construction of KAGRA, including the study of the advantages and challenges of building an underground detector and the operation of the iKAGRA interferometer together with the geophysics interferometer that has been constructed in the same tunnel.

The Baryon Acoustic Oscillations (BAO) refer to the ripples of material density in the Universe. As the most direct density tracers in the universe, galaxies have been commonly used in studies of BAO peak detection. The spatial number density of galaxies, to a certain extent, reflects the distribution of the material density of our Universe. Using galaxies as matter tracers, we can construct more overlapping empty spheres (DT voids) than the matter tracers, via Delaunay Triangulation technique. We show that their radii excellently reflect the galaxy number density round them, and they can serve as reliable different density region tracers. Using the data from an unprecedented large-scale $N$-body simulation "TianNu", we conduct some fundamental statistical studies and clustering analysis of the DT voids. We discuss in detail the representative features of two-point correlation functions of different DT void populations. We show that the peak, the position of which corresponds to the average radius of data samples, is the most representative feature of the two-point correlation function of the DT voids. In addition, we also construct another voids, the disjoint voids, and investigate their some statistical properties and clustering properties. And we find that the occupied space of all disjoint voids accounts for about $45\%$ of the volume of the simulation box, regardless of the number density of mock galaxies. We also investigate the BAO detections based on different tracers, i.e. mock galaxies, low-density region tracers, and high-density region tracers respectively. Our results show that BAO intensities detected by low/high-density region tracers are enhanced significantly compared to the BAO detection by mock galaxies, for the mock galaxy catalogue with the number density of $7.52\times10^{-5}$ $h^3$ Mpc$^{-3}$.

The low braking-index pulsar PSR J1734$-$3333 could be born with superhigh internal magnetic fields $B_{\rm in}\sim10^{15}-10^{16}$ G, and undergo a supercritical accretion soon after its formation in a supernova explosion. The buried multipole magnetic fields will merger into a dipole magnetic field. Since the magnetic flow transfers from the core to the crust of the pulsar, its surface dipole field grows quickly at a power-law form assumed until it saturates at the level of internal dipole field. The increase in surface dipole magnetic field results in the observed low braking index of $n=0.9(2)$. Keeping an average field growth index $\varepsilon=1.34(6)$, this pulsar will become a magnetar with surface dipole magnetic field at the equator $B_{\rm d}\sim 2.6(1)\times 10^{14}$\,G and $\sim 5.3(2)\times 10^{14}$\,G after next 50\,kyrs and 100\,kys, respectively.