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.

The interstellar oxygen isotopic ratio of 18O/17O can reflect the relative amount of the secular enrichment by ejecta from high-mass versus intermediate-mass stars. Previous observations found a Galactic gradient of 18O/17O, i.e., low ratios in the Galactic center and large values in the Galactic disk, which supports the insideout formation scenario of our Galaxy. However, the observed objects are not many and, in particular, not so many at large galactocentric distances. Thus we started a systematic study on Galactic interstellar 18O/17O, through observations of C18O and C17O multi-transition lines toward a large sample of 286 sources (at least one order of magnitude larger than previous ones), from the Galactic center region to the far outer Galaxy (~22 kpc). Here we present our observations of J=1-0 lines of C18O and C17O, with the ARO12m and the IRAM 30m telescope. We detected successfully both C18O and C17O 1-0 lines for 34 sources among our IRAM30m sample of 50 targets and for 166 sources among our ARO12m sample of 260 targets. The C18O optical depth effect on our ratio results, evaluated by fitting results of C17O spectra with hyperfine components and our RADEX non-LTE model calculation for the strongest source, was found to be insignificant. Beam dilution does not seem to be a problem either, which was supported by the fact of no systematic variation between the isotopic ratio and the heliocentric distance, and consistent measured ratios from two telescopes for most of those detected sources. Our results, though there are still very few detections made for sources in the outer Galaxy, confirm the apparent 18O/17O gradient of 18O/17O = (0.10+-0.03)R_GC+(2.95+-0.30), with a Pearson's rank correlation coefficient R = 0.69. This is supported by the newest Galactic chemical evolution model including the impact of massive stellar rotators and novae.

The origin and evolution of spectral curvature in blazar spectral energy distribution (SED) is still unclear. Since the observed SED curvature is related to an intrinsic curvature in emitting electron energy distribution (EED), we study this question by employing a log-parabolic EED with a curvature parameter and peak energy to model the quasi-simultaneous broadband SEDs of selected blazars in Fermi-LAT Bright AGN Sample (LBAS) using synchrotron and inverse Compton (IC) processes. We find that log-parabolic IC model can successfully explain the emission in all blazars in our sample. On average, FSRQs have higher magnetic field, Doppler factor, and curvature than BL Lac objects. The BL Lac objects show an anticorrelation between the curvature parameter of the EED and its peak energy, which is a signature of stochastic acceleration. FSRQs do not manifest such correlation and rather show a mild positive relationship between these parameters. This suggests that the evolution of spectral curvature in the BL Lac objects is dominated by a strong stochastic acceleration component, whereas the curvature in FSRQs evolves in a cooling dominated regime due to an additional external Compton (EC) component. The strong cooling in FSRQs not only restricts the electron peak energy but also adds extra curvature to the high energy tail of emitting EED. Since the curvature decreases from FSRQs toward high peak BL Lac objects (HBLs), opposite to peak energy, the curvature parameter can be considered a third parameter of the blazar sequence in addition to peak frequency and luminosity.

We present optical-infrared photometric and spectroscopic observations of Gaia18dvy, located in the Cygnus OB3 association at a distance of 1.88 kpc. The object was noted by the Gaia alerts system when its lightcurve exhibited a $\gtrsim$4 mag rise in 2018-2019. The brightening was also observable at mid-infared wavelengths. The infrared colors of Gaia18dvy became bluer as the outburst progressed. Its optical and near-infrared spectroscopic characteristics in the outburst phase are consistent with those of bona fide FU Orionis-type young eruptive stars. The progenitor of the outburst is probably a low-mass K-type star with an optical extinction of $\sim$3 mag. A radiative transfer modeling of the circumstellar structure, based on the quiescent spectral energy distribution, indicates a disk with a mass of $4{\times}10^{-3}\,M_{\odot}$. Our simple accretion disk modeling implies that the accretion rate had been exponentially increasing for more than 3 years until mid-2019, when it reached a peak value of $6.9 \times 10^{-6}\,M_{\odot}$yr$^{-1}$. In many respects, Gaia18dvy is similar to the FU Ori-type object HBC 722

In Gaia DR2, the unprecedented high-precision level reached in sub-mas for astrometry and mmag for photometry. Using cluster members identified with these astrometry and photometry in Gaia DR2, we can obtain a reliable determination of cluster properties. However, because of the shortcoming of Gaia spectroscopic observation in dealing with densely crowded cluster region, the number of radial velocity and metallicity for cluster member stars from Gaia DR2 is still lacking. In this study, we aim to improve the cluster properties by combining the LAMOST spectra. In particular, we provide the list of cluster members with spectroscopic parameters as an add-value catalog in LAMOST DR5, which can be used to perform detailed study for a better understanding on the stellar properties, by using their spectra and fundamental properties from the host cluster. We cross-matched the spectroscopic catalog in LAMOST DR5 with the identified cluster members in Cantat-Gaudin et al.2018 and then used members with spectroscopic parameters to derive statistical properties of open clusters. We obtained a list of 8811 members with spectroscopic parameters and a catalog of 295 cluster properties. In addition, we study the radial and vertical metallicity gradient and age-metallicity relation with the compiled open clusters as tracers, finding slopes of -0.053$\pm$0.004 dex kpc$^{-1}$, -0.252$\pm$0.039 dex kpc$^{-1}$ and 0.022$\pm$0.008 dex Gyr$^{-1}$, respectively. Both slopes of metallicity distribution relation for young clusters (0.1 Gyr < Age < 2 Gyr) and the age-metallicity relation for clusters within 6 Gyr are consistent with literature results. In order to fully study the chemical evolution history in the disk, more spectroscopic observations for old and distant open clusters are needed for further investigation.

Dark energy (DE) plays an important role in the expansion history of our universe. In the earlier works, universe including DE with a constant or time-varying equation of state (EOS) has been well studied. Different from predecessors, we investigate the cosmic expansion with multicomponent DE. Today's total energy density of DE is separated into several parts equally and every part has a constant EOS $w_i$. We modify the Friedmann equation and the parameterized post-Friedmann description of dark energy, then put constraints on $w_i$ from Planck 2018 TT,TE,EE$+$lowE$+$lensing, BAO data, PANTHEON samples and SH0ES measurement. We find that the largest and the smallest values of $w_i$ have no overlap at $95\%$C.L.. Our results indicate that DE is composed of more than one candidate. Besides, the density evolution of multicomponent DE is tick-like and different from the other DE models'. That is to say, multicomponent DE makes a contribution to late and early universe.

Newton's constant has observational effects on both the CMB power spectra and the light curves of SNIa. We use Planck data, BAO data and the SNIa measurement to constrain the varying Newton's constant $G$ during the CMB epoch and the redshift ranges of PANTHEON samples, and find no evidence indicating that $G$ is varying with redshift. By extending the $\Lambda$CDM model with one free parameter $G$, we get $G =(6.65635_{-0.18560}^{+0.18766} ) \times 10^{-11} \rm m^3kg^{-1}s^{-2}$ and $H_0=67.62^{+1.24}_{-1.25} $ km s$^{-1}$ Mpc$^{-1}$ at 68$\%$ C.L. from Planck$+$BAO$+$uncalibrated PANTHEON. The results show the value of $G$ is consistent with CODATA 2018, but the $H_0$ tension can't be solved in this way.

In this work, we investigated the X-ray emission for a sample of young radio AGNs by combining their data from Chandra/XMM-Newton and at other wavebands. We find strong correlations between the X-ray luminosity $L_{\rm X}$ in 2$-$10 keV and the radio luminosities $L_{\rm R}$ at 5 GHz for VLBI radio-core, VLA radio-core and FIRST component, indicating that both pc- and kpc-scale radio emission strongly correlate with X-ray emission in these sources. We find approximately linear dependence of radio on X-ray luminosity in the sources with radiative efficient accretion flows (i.e., the Eddington ratio $R_{\rm edd} \gtrsim 10^{-3} $) with b $\sim$ 1 ($L_{\rm R}$ $\propto$ $L_{\rm X} ^{b}$) and $\xi_{\rm RX}$ $\sim$ 1 in fundamental plane using VLBI data, where the dependence is consistent with the re-analysed result on the previous study in \cite{2016ApJ...818..185F} at $R_{\rm edd} \gtrsim 10^{-3}$, however is significantly deviated from the theoretical prediction of accretion flow as the origin of X-ray emission. In contrast to radio-quiet quasars, there is no significant correlation between $\Gamma$ and Eddington ratio. Our results seem to indicate that the X-ray emission of high-accreting young radio AGNs may be from jet. We constructed the SEDs for 18 sources (most are in radiative efficient accretion) including 9 galaxies and 9 quasars with high-quality X-ray data, and find that the X-ray emission of most quasars is more luminous than that of normal radio-quiet quasars. This is clearly seen from the quasar composite SED, of which the X-ray emission is apparently higher than that of radio-quiet quasars, likely supporting the jet-related X-ray emission in young radio AGNs. The scenario that the X-ray emission is from self-synchrotron Compton (SSC) is discussed.

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.

FU Orionis-type stars are young stellar objects showing large outbursts due to highly enhanced accretion from the circumstellar disk onto the protostar. FUor-type outbursts happen in a wide variety of sources from the very embedded ones to those with almost no sign of extended emission beyond the disk. The subsequent eruptions might gradually clear up the obscuring envelope material and drive the protostar on its way to become a disk-only T Tauri star. We used VLT/VISIR to obtain the first spectra that cover the 8-13 $\mu$m mid-infrared wavelength range in low-resolution of five recently discovered FUors. Four objects from our sample show the 10 $\mu$m silicate feature in emission. We study the shape and strength of the silicate feature in these objects and find that they mostly contain large amorphous grains, suggesting that large grains are typically not settled to the midplane in FUor disks. This is a general characteristic of FUors, as opposed to regular T Tauri-type stars whose disks display anything from pristine small grains to significant grain growth. We classify our targets by determining whether the silicate feature is in emission or in absorption, and confront them with the evolutionary scenarios on the dispersal of the envelopes around young stars. In our sample, all Class II objects exhibit silicate emission, while for Class I objects, the appearance of the feature in emission or absorption depends on the viewing angle with respect to the outflow cavity. This highlights the importance of geometric effects when interpreting the silicate feature.

Since September 2018, LAMOST starts a new 5-year medium-resolution spectroscopic survey (MRS) using bright/gray nights. We present the scientific goals of LAMOST-MRS and propose a near optimistic strategy of the survey. A complete footprint is also provided. Not only the regular medium-resolution survey, but also a time-domain spectroscopic survey is being conducted since 2018 and will be end in 2023. According to the detailed survey plan, we expect that LAMOST-MRS can observe about 2 million stellar spectra with ~7500 and limiting magnitude of around G=15 mag. Moreover, it will also provide about 200 thousand stars with averagely 60-epoch observations and limiting magnitude of G~14 mag. These high quality spectra will give around 20 elemental abundances, rotational velocities, emission line profiles as well as precise radial velocity with uncertainty less than 1 km/s. With these data, we expect that LAMOST can effectively leverage sciences on stellar physics, e.g. exotic binary stars, detailed observation of many types of variable stars etc., planet host stars, emission nebulae, open clusters, young pre-main-sequence stars etc.

We present the first statistical study of magnetic structures and associated energy dissipation observed during a single period of turbulent magnetic reconnection, by using the \textit{in-situ} measurements of the Magnetospheric Multiscale mission in the Earth's magnetotail on July 26, 2017. The structures are selected by identifying a bipolar signature in the magnetic field and categorized as plasmoids or current sheets via an automated algorithm which examines current density and plasma flow. The size of the plasmoids forms a decaying exponential distribution ranging from sub-electron up to ion scales. The presence of substantial number of current sheets is consistent with a physical picture of dynamic production and merging of plasmoids during turbulent reconnection. The magnetic structures are locations of significant energy dissipation via electric field parallel to the local magnetic field, while dissipation via perpendicular electric field dominates outside of the structures. Significant energy also returns from particles to fields.

We model clusters of axions with spherically symmetric momentum but arbitrary spatial distributions and study the directional profile of photos produced in their evolution through spontaneous and stimulated decay of axions via the process $a \rightarrow \gamma + \gamma$. Several specific examples are presented.

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.

Identifying the sources of the highest energy cosmic rays requires understanding how they are deflected by the stochastic, spatially intermittent intergalactic magnetic field. Here we report measurements of energetic charged-particle propagation through a laser-produced magnetized plasma with these properties. We characterize the diffusive transport of the particles experimentally. The results show that the transport is diffusive and that, for the regime of interest for the highest-energy cosmic rays, the diffusion coefficient is unaffected by the spatial intermittency of the magnetic field.

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.

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

We present a systematic analysis of the phase lags associated with the type-C QPOs in GRS 1915+105 using RXTE data. Our sample comprises of 620 RXTE observations with type-C QPOs ranging from ~0.4 Hz to ~6.3 Hz. Based on our analysis, we confirm that the QPO phase lags decrease with QPO frequency, and change sign from positive to negative at a QPO frequency of ~2 Hz. In addition, we find that the slope of this relation is significantly different between QPOs below and above 2 Hz. The relation between the QPO lags and QPO rms can be well fitted with a broken line: as the QPO lags go from negative to positive, the QPO rms first increases, reaching its maximum at around zero lag, and then decreases. The phase-lag behaviour of the subharmonic of the QPO is similar to that of the QPO fundamental, where the subharmonic lags decrease with subharmonic frequency and change sign from positive to negative at a subharmonic frequency of ~1 Hz; on the contrary, the second harmonic of the QPO shows a quite different phase-lag behaviour, where all the second harmonics show hard lags that remain more or less constant. For both the QPO and its (sub)harmonics, the slope of the lag-energy spectra shows a similar evolution with frequency as the average phase lags. This suggests that the lag-energy spectra drives the average phase lags. We discuss the possibility for the change in lag sign, and the physical origin of the QPO lags.

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

We present our results from a mid-infrared interferometric survey targeted at the planet-forming region in the circumstellar disks around low- and intermediate-mass young stars. Our sample consist of 82 objects, including T Tauri stars, Herbig Ae stars, and young eruptive stars. Our main results are: 1) Disks around T Tauri stars are similar to those around Herbig Ae stars, but are relatively more extended once we account for stellar luminosity. 2) From the distribution of the sizes of the mid-infrared emitting region we find that inner dusty disk holes may be present in roughly half of the sample. 3) Our analysis of the silicate spectral feature reveals that the dust in the inner $\sim$1 au region of disks is generally more processed than that in the outer regions. 4) The dust in the disks of T Tauri stars typically show weaker silicate emission in the N band spectrum, compared to Herbig Ae stars, which may indicate a general difference in the disk structure. Our data products are available at VizieR, and at the following web page: http://konkoly.hu/MIDI_atlas

Analysis of model fit results of 15,210 electron velocity distribution functions (VDFs), observed within $\pm$2 hours of 52 interplanetary (IP) shocks by the Wind spacecraft near 1 AU, is presented as the third and final part on electron VDFs near IP shocks. The core electrons and protons dominate in the magnitude and change in the partial-to-total thermal pressure ratio, with the core electrons often gaining as much or more than the protons. Only a moderate positive correlation is observed between the electron temperature and the kinetic energy change across the shock, while weaker, if any, correlations were found with any other macroscopic shock parameter. No VDF parameter correlated with the shock normal angle. The electron VDF evolves from a narrowly peaked core with flaring suprathermal tails in the upstream to either a slightly hotter core with steeper tails or much hotter flattop core with even steeper tails downstream of the weaker and strongest shocks, respectively. Both quasi-static and fluctuating fields are examined as possible mechanisms modifying the VDF but neither is sufficient alone. For instance, flattop VDFs can be generated by nonlinear ion acoustic wave stochastic acceleration (i.e., inelastic collisions) while other work suggested they result from the combination of quasi-static and fluctuating fields. This three-part study shows that not only are these systems not thermodynamic in nature, even kinetic models may require modification to include things like inelastic collision operators to properly model electron VDF evolution across shocks or in the solar wind.

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.

In recent years, several optimal dynamos have been discovered. They minimize the magnetic energy dissipation or, equivalently, maximize the growth rate at a fixed magnetic Reynolds number. In the optimal dynamo of Willis (2012, Phys. Rev. Lett. 109, 251101), we find mean-field dynamo action for planar averages. One component of the magnetic field grows exponentially while the other decays in an oscillatory fashion near onset. This behavior is different from that of an alpha^2 dynamo, where the two non-vanishing components of the planar averages are coupled and have the same growth rate. For the Willis dynamo, we find that the mean field is excited by a negative turbulent magnetic diffusivity, which has a non-uniform spatial profile near onset. The temporal oscillations in the decaying component are caused by the corresponding component of the diffusivity tensor being complex when the mean field is decaying and, in this way, time-dependent. The growing mean field can be modeled by a negative magnetic diffusivity combined with a positive magnetic hyperdiffusivity. In two other classes of optimal dynamos of Chen et al. (2015, J. Fluid Mech. 783, 23), we find, to some extent, similar mean-field dynamo actions. When the magnetic boundary conditions are mixed, the two components of the planar averaged field grow at different rates when the dynamo is 15% supercritical. When the mean magnetic field satisfies homogeneous boundary conditions (where the magnetic field is tangential to the boundary), mean-field dynamo action is found for one-dimensional averages, but not for planar averages. Despite having different spatial profiles, both dynamos show negative turbulent magnetic diffusivities. Our finding suggests negative turbulent magnetic diffusivities may support a broader class of dynamos than previously thought, including these three optimal dynamos.

The Herbig Ae star HD 169142 is known to have a gaseous disk with a large inner hole, and also a photometrically variable inner dust component in the sub-au region. Following up our previous analysis, we further studied the temporal evolution of inner dust around HD 169142, which may provide information on the evolution from late-stage protoplanetary disks to debris disks. We used near-infrared interferometric observations obtained with VLTI/PIONIER to constrain the dust distribution at three epochs spanning six years. We also studied the photometric variability of HD 169142 using our optical-infrared observations and archival data. Our results indicate that a dust ring at ~0.3 au formed at some time between 2013 and 2018, and then faded (but did not completely disappear) by 2019. The short-term variability resembles that observed in extreme debris disks, and is likely related to short-lived dust of secondary origin, though variable shadowing from the inner ring could be an alternative interpretation. If confirmed, this is the first direct detection of secondary dust production inside a protoplanetary disk.

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

An essential step to understanding protoplanetary evolution is the study of disks that contain gaps or inner holes. The pretransitional disk around the Herbig star HD 169142 exhibits multi-gap disk structure, differentiated gas and dust distribution, planet candidates, and near-infrared fading in the past decades, which make it a valuable target for a case study of disk evolution. Using near-infrared interferometric observations with VLTI/PIONIER, we aim to study the dust properties in the inner sub-au region of the disk in the years 2011-2013, when the object is already in its near-infrared faint state. We first performed simple geometric modeling to characterize the size and shape of the NIR-emitting region. We then performed Monte-Carlo radiative transfer simulations on grids of models and compared the model predictions with the interferometric and photometric observations. We find that the observations are consistent with optically thin gray dust lying at Rin ~ 0.07 au, passively heated to T ~ 1500 K. Models with sub-micron optically thin dust are excluded because such dust will be heated to much higher temperatures at similar distance. The observations can also be reproduced with a model consisting of optically thick dust at Rin ~ 0.06 au, but this model is plausible only if refractory dust species enduring ~2400 K exist in the inner disk.

Axion remains the most convincing solution to the strong-CP problem and a well-motivated dark matter candidate, causing the search for axions and axion-like particles(ALPs) to attract attention continually. The exchange of such particles may cause anomalous spin-dependent forces, inspiring many laboratory ALP searching experiments based on the detection of macroscopic monopole-dipole interactions between polarized electrons/nucleons and unpolarized nucleons. Since there is no exact proof of the existence of these interactions, to detect them is still of great significance. In the present paper, we study the electron-neucleon monopole-dipole interaction with a new method, in which a hybrid spin-nanocantilever optomechanical system consisting of a nitrogen-vacancy(NV) center and a nanocantilever resonator is used. With a static magnetic field and a pump microwave beam and a probe microwave beam applied, a probe absorption spectrum could be obtained. Through specific peaks appearing in the spectrum, we can identify this monopole-dipole interaction. And we also provide a prospective constraint to constrain the interaction. Furthermore, because our method can also be applied to the detection of some other spin-dependent interactions, this work provides new ideas for the experimental searches of the anomalous spin-dependent interactions.

The radial velocity (RV) is a basic physical quantity which can be determined through Doppler shift of the spectrum of a star. The precision of RV measurement depends on the resolution of the spectrum we used and the accuracy of wavelength calibration. In this work, radial velocities of LAMOST-II medium resolution (R ~ 7500) spectra are measured for 1,594,956 spectra (each spectrum has two wavebands) through matching with templates. A set of RV standard stars are used to recalibrate the zero point of the measurement, and some reference sets with RVs derived from medium/high-resolution observations are used to evaluate the accuracy of the measurement. Comparing with reference sets, the accuracy of our measurement can get 0.0227 km s/1 with respect to radial velocities standard stars. The intrinsic precision is estimated with the multiple observations of single stars, which can achieve to 1.36 km s/1,1.08 km s/1, 0.91 km s/1 for the spectra at signal-to-noise levels of 10, 20, 50, respectively.

We search through $\gamma$-ray data obtained with the Large Area Telescope (LAT) onboard the Fermi Gamma-Ray Space Telescope and find 24 blazars (or candidates) that have a single clear flare event in their 9.5 year long-term light curves. We define these events as loner flares since each flare stands out significantly above the relatively stable, low-flux light curve. We analyze the LAT data in detail for these 24 sources. The flares in ten of them are primarily due to a single sharp peak, for which we study by fitting with two different analytic functions. The time durations thus determined for the sharp peaks are in a range of 4-25 days. The $\gamma$-ray spectra of the 24 blazar sources can be described with a power-law or a log-parabola function. We obtain their spectral properties in the flaring and quiescent states, and find that in the flares 16 of the sources have harder emission and three have softer emission while the other five keep the same emission. We discuss a possible correlation between the differences in photon index in the quiescent and flaring states and photon indices in quiescence. In addition, the sharp peak flares seem to have a tendency of having long time durations and hard emission, possibly related to their physical origin in a blazar jet. Studies of more similar flares will help establish these possible features.

EX Lup is the prototype of a class of pre-main sequence eruptive stars defined by their repetitive outbursts lasting several months. In 2008 January-September EX Lup underwent its historically largest outburst, brightening by about 4 magnitudes in visual light. In previous studies we discovered on-going silicate crystal formation in the inner disk during the outburst, but also noticed that the measured crystallinity fraction started decreasing after the source returned to the quiescent phase. Here we present new observations of the 10 $\mu$m silicate feature, obtained with the MIDI and VISIR instruments at Paranal Observatory. The observations demonstrate that within five years practically all crystalline forsterite disappeared from the surface of the inner disk. We reconstruct this process by presenting a series of parametric axisymmetric radiative transfer models of an expanding dust cloud that transports the crystals from the terrestrial zone to outer disk regions where comets are supposed to form. Possibly the early Sun also experienced similar flare-ups, and the forming planetesimals might have incorporated crystalline silicate material produced by such outbursts. Finally, we discuss how far the location of the dust cloud could be constrained by future JWST observations.

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

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 Hard X-ray Modulation Telescope ($Insight$-HXMT) was successfully launched on June 15th, 2017. It performs broad band X-ray scan survey of the Galactic Plane to detect new black holes and other objects in active states. It also observes X-ray binaries to study their X-ray variabilities. Here we will introduce the Science Operations of $Insight$-HXMT, which is responsible for collecting and evaluating observation proposals, scheduling observations, and monitoring the working status of the payloads.

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

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

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

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

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

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

Using the very recently reported mass $2.14^{+0.10}_{-0.09}M_\odot$ of PSR J0740+6620 together with the data of finite nuclei and the constraints on the equation of state of symmetric nuclear matter at suprasaturation densities from flow data in heavy-ion collisions, we show that the symmetry energy $E_{\rm sym}(n)$ cannot be supersoft so that it becomes negative at suprasaturation densities in neutron stars (NSs) and thus may make the NS have a pure neutron matter core. This is in contrast to the fact that using the mass $2.01 \pm 0.04 M_\odot$ of PSR J0348+0432 as the NS maximum mass cannot rule out the supersoft high-density $E_{\rm sym}(n)$. Furthermore, we find the stiffer high-density $E_{\rm sym}(n)$ based on the existence of $2.14M_\odot$ NSs leads to a strong constraint of $\Lambda_{1.4} \ge 348^{+88}_{-51}$ for the dimensionless tidal deformability of the canonical $1.4M_\odot$ NS.

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.

A statistical analysis of 15,210 electron velocity distribution function (VDF) fits, observed within $\pm$2 hours of 52 interplanetary (IP) shocks by the $Wind$ spacecraft near 1 AU, is presented. This is the second in a three-part series on electron VDFs near IP shocks. The electron velocity moment statistics for the dense, low energy core, tenuous, hot halo, and field-aligned beam/strahl are a statistically significant list of values illustrated with both histograms and tabular lists for reference and baselines in future work. The beam/strahl fit results in the upstream are currently the closest thing to a proper parameterization of the beam/strahl electron velocity moments in the ambient solar wind. This work will also serve as a 1 AU baseline and reference for missions like $Parker \ Solar \ Probe$ and $Solar \ Orbiter$. The median density, temperature, beta, and temperature anisotropy values for the core(halo)[beam/strahl] components, with subscripts $ec$($eh$)[$eb$], of all fit results respectively are $n{\scriptstyle_{ec(h)[b]}}$ $\sim$ 11.3(0.36)[0.17] $cm^{-3}$, $T{\scriptstyle_{ec(h)[b], tot}}$ $\sim$ 14.6(48.4)[40.2] $eV$, $\beta{\scriptstyle_{ec(h)[b], tot}}$ $\sim$ 0.93(0.11)[0.05], and $\mathcal{A}{\scriptstyle_{ec(h)[b]}}$ $\sim$ 0.98(1.03)[0.93]. The nuanced details of the fitting method and data product description were published in Paper I and the detailed analysis of the results will be shown in Paper III.

Both the neutron skin thickness $\Delta R_{np}$ of atomic nuclei and the low-energy neutrino-nucleon ($\nu N$) interactions are of fundamental importance in nuclear and particle physics, astrophysics as well as new physics beyond the standard model (SM) but largely uncertain currently, and the coherent elastic neutrino-nucleus scattering (CE$\nu$NS) provides a clean way to extract their information. New physics beyond the SM may cause effectively a shift of the SM weak mixing angle $\theta_W$ in low-energy $\nu N$ interactions, leading to an effective weak mixing angle $\theta^*_W$. By analyzing the CE$\nu$NS data of the COHERENT experiment, we find that while a one-parameter fit to the COHERENT data by varying $\Delta R_{np}$ produces $\Delta R^{\rm{CsI}}_{np} \simeq 0.68^{+0.91}_{-1.13}$ fm for CsI with an unrealistically large central value by fixing $\sin^2 \theta^*_W$ at the low-energy SM value of $\sin^2\theta_W^{\rm{SM}} = 0.23857$, a two-dimensional fit by varying $\Delta R_{np}$ and $\sin^2 \theta^*_W$ leads to a strong positive correlation between $\Delta R_{np}$ and $\sin^2 \theta^*_W$ with significantly smaller central values of $\Delta R^{\rm {CsI}}_{np} \simeq 0.24_{-2.03}^{+2.30}$ fm and $\sin^2 \theta^*_W = 0.21_{-0.10}^{+0.13}$. Although the uncertainty is too large to claim a determination of $\Delta R^{\rm{CsI}}_{np}$ and $\sin^2 \theta^*_W$, the present study suggests that the multi-dimensional fit is important in future analyses of high-precision CE$\nu$NS data. The implication of the possible deviation of $\sin^2 \theta^*_W$ from $\sin^2\theta_W^{\rm{SM}}$ on new physics beyond the SM is also discussed.

The CFHT Large Area U-band Deep Survey (CLAUDS) uses data taken with the MegaCam mosaic imager on CFHT to produce images of 18.60 deg2 with median seeing of FWHM=0.92 arcsec and to a median depth of U = 27.1 AB (5 sigma in 2 arcsec apertures), with selected areas that total 1.36 deg2 reaching a median depth of U=27.7 AB. These are the deepest U-band images assembled to date over this large an area. These data are located in four fields also imaged to comparably faint levels in grizy and several narrowband filters as part of the Hyper Suprime-Cam (HSC) Subaru Strategic Program (HSC-SSP). These CFHT and Subaru datasets will remain unmatched in their combination of area and depth until the advent of the Large Synoptic Survey Telescope (LSST). This paper provides an overview of the scientific motivation for CLAUDS and gives details of the observing strategy, observations, data reduction, and data merging with the HSC-SSP. Three early applications of these deep data are used to illustrate the potential of the dataset: deep U-band galaxy number counts, z~3 Lyman break galaxy (LBG) selection, and photometric redshifts improved by adding CLAUDS U to the Subaru HSC grizy photometry.

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

Chorus-like whistler-mode waves that are known to play a fundamental role in driving radiation-belt dynamics are excited on the Large Plasma Device by the injection of a helical electron beam into a cold plasma. The mode structure of the excited whistler wave is identified using a phase-correlation technique showing that the waves are excited through a combination of Landau resonance, cyclotron resonance and anomalous cyclotron resonance. The dominant wave mode excited through cyclotron resonance is quasi-parallel propagating, whereas wave modes excited through Landau resonance and anomalous cyclotron resonance propagate at oblique angles that are close to the resonance cone. An analysis of the linear wave growth rates captures the major observations in the experiment. The results have important implications for the generation process of whistler waves in the Earth's inner magnetosphere.

A multi-institutional, multi-national science team will soon submit a NASA proposal to build a constellation of spacecraft to fly into the near-Earth solar wind in a swarm spanning a multitude of scales in order to obtain critically needed measurements that will reveal the underlying dynamics of magnetized turbulence. This white paper, submitted to the Plasma 2020 Decadal Survey Committee, provides a brief overview of turbulent systems that constitute an area of compelling plasma physics research, including why this mission is needed, and how this mission will achieve the goal of revealing how energy is transferred across scales and boundaries in plasmas throughout the universe.

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 present new catalogs of M giant and M dwarf stars from the LAMOST DR5. In total, 39,796 M giants and 501,152 M dwarfs are identified from the classification pipeline. The template-fitting results contain M giants with 7 temperature subtypes from M0 to M6, M dwarfs with 18 temperature subtypes from K7.0 to M8.5 and 12 metallicity subclasses from dMr to usdMp. We cross-matched our M-type catalog with the 2MASS and WISE catalog to obtain infrared magnitude and colors. Adopting the distances derived from the parallaxes in \gaia{} DR2, the M_G vs. (G_bp-G_rp)_0 diagram shows that there are also early-type stars and white dwarf-M dwarf binaries included in our M type stars sample, with a contamination rate of about 4.6% for M giants and 0.48% for M dwarfs. We found that CaH spectral indices are an efficient selection criteria for carbon stars. A total of 289 carbon stars were identified from the M giants sample, and further confirmed by LAMOST spectra.

Due to ionosphere absorption and the interference by natural and artificial radio emissions, ground observation of the sky at the decameter or longer is very difficult. This unexplored part of electromagnetic spectrum has the potential of great discoveries, notably in the study of cosmic dark ages and dawn, but also in heliophysics and space weather, planets, cosmic ray and neutrinos, pulsar and interstellar medium, extragalactic radio sources, and even SETI. At a forum organized by the International Space Science Institute-Beijing (ISSI-BJ), we discussed the prospect of opening up this window for astronomical observations by using a constellation of small or micro-satellites. We discussed the past experiments and the current ones such as the low frequency payload on Chang'e-4 mission lander, relay satellite and the Longjiang satellite, and also the future DSL mission, which is a linear array on lunar orbit which can make synthesized map of the whole sky as well as measure the global spectrum. We also discuss the synergy with other experiments, including ground global experiments such as EDGES, SARAS, SCI-HI and High-z, PRIZM/Albatros, ground imaging facillities such as LOFAR and MWA, and space experiments such as SUNRISE, DARE/DAPPER and PRATUSH. We also discussed some technical aspects of the DSL concept.

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 Ultra-Long Wavelength (ULW) regime of longer than 10 m (corresponding to frequencies below 30 MHz), remains as the last virtually unexplored window in radio astronomy, and is presently attracting considerable attention as an area of potentially rewarding studies. However, the opaqueness of the Earth's ionosphere makes the ULW celestial radio emission very difficult to detect with ground-based instrumentation. The impact of the ionosphere on ULW radio emission depends on the Solar cycle activities and varies with time. In addition, the ULW spectrum region is densely populated by intensive artificial radio frequency interference (RFI). An obvious solution of these problems is to place an ULW radio telescope in space. However, this solution is expensive and poses non-negligible technological challenges. An alternative approach is triggered by recent studies showing that the period of post 2020 will be most suitable for exploratory ground-based ULW radio observations due to the expected 'calm' state of the ionosphere; the ionospheric cutoff-frequency could be well below 10 MHz, even in the day time. In anticipation of this upcoming opportunity, we propose and present in this paper a concept of an experimental ULW radio array, with the intention of setting it up in Inner Mongolia, China. This ULW facility will use the infrastructure of the currently operational Miangtu Spectra Radio Heliograph (MUSER). The proposed ULW array covers the frequency range from 1 to 72 MHz. This experimental array will be used for exploratory studies of celestial radio emission in the ULW range of the spectrum.

We studied the long-term evolution of the spectral-temporal correlated properties of the black hole candidate Swift J1753.5-0127 from the onset of its outburst until 2011 with the Rossi X-ray Timing Explorer (RXTE). The source stayed most of its lifetime during hard state, with occasionally transitioned to the hard intermediate state. Similar to typical black hole transients, Swift J1753.5-0127 traces a clear hard line in absolute rms - intensity diagram during the low hard state, with expected highest absolute rms, while shows a clear turn during the hard intermediate state, accompanied by lower absolute rms. Different from Cyg X-1, we found that frequency-dependent time lag increased significantly in the 0.02 -- 3.2 Hz band during state transition in this source. The X-ray time lags in 0.02 -- 3.2 Hz can therefore be used as indicators of state transition in this source. Type-C QPO frequency is positively related with its fractional rms and X-ray photon index, suggesting a moving inwards disk/corona scenario. We discussed the physical interpretation of our results in our paper.

While the equation of state (EOS) of symmetric nuclear matter (SNM) at suprasaturation densities has been relatively well constrained from heavy-ion collisions, the EOS of high-density neutron-rich matter is still largely uncertain due to the poorly known high-density behavior of the symmetry energy. Using the constraints on the EOS of SNM at suprasaturation densities from heavy-ion collisions together with the data of finite nuclei and the existence of $2M_\odot$ neutron stars from electromagnetic (EM) observations, we show that the high-density symmetry energy cannot be too soft, which leads to lower bounds on dimensionless tidal deformability of $\Lambda_{1.4} \ge 193$ and radius of $R_{1.4} \ge 11.1$ km for $1.4M_\odot$ neutron star. Furthermore, we find that the recent constraint of $\Lambda_{1.4} \le 580$ from the gravitational wave signal GW170817 detected from the binary neutron star merger by the LIGO and Virgo Collaborations rules out too stiff high-density symmetry energy, leading to an upper limit of $R_{1.4} \le 13.3$ km. All these terrestrial nuclear experiments and astrophysical observations based on strong, EM and gravitational measurements together put stringent constraints on the high-density symmetry energy and the EOS of SNM, pure neutron matter and neutron star matter.

The lattice Hamiltonian method is developed for solving the Vlasov equation with nuclear mean-field based on the Skyrme pseudopotential up to next-to-next-to-next-to leading order. The ground states of nuclei are obtained through varying the total energy with respect to the density distribution of nucleons. Owing to the self-consistent treatment of initial nuclear ground state and the exact energy conservation in the lattice Hamiltonian method, the present framework of solving the Vlasov equation exhibits very stable nuclear ground state evolution. As a first application of the new lattice Hamiltonian Vlasov method, we explore the iso-scalar giant monopole and iso-vector giant dipole modes of finite nuclei. The obtained results are shown to be comparable to that from random-phase approximation and consistent with the experimental data, indicating the capability of the present method in dealing with the long-time near-equilibrium nuclear dynamics.

Analysis of 15314 electron velocity distribution functions (VDFs) within $\pm$2 hours of 52 interplanetary (IP) shocks observed by the \emph{Wind} spacecraft near 1 AU are introduced. The electron VDFs are fit to the sum of three model functions for the cold dense core, hot tenuous halo, and field-aligned beam/strahl component. The best results were found by modeling the core as either a bi-kappa or a symmetric (or asymmetric) bi-self-similar velocity distribution function, while both the halo and beam/strahl components were best fit to bi-kappa velocity distribution function. This is the first statistical study to show that the core electron distribution is better fit to a self-similar velocity distribution function than a bi-Maxwellian under all conditions. The self-similar distribution deviation from a Maxwellian is a measure of inelasticity in particle scattering from waves and/or turbulence. The range of values defined by the lower and upper quartiles for the kappa exponents are $\kappa{\scriptstyle_{ec}}$ $\sim$ 5.40--10.2 for the core, $\kappa{\scriptstyle_{eh}}$ $\sim$ 3.58--5.34 for the halo, and $\kappa{\scriptstyle_{eb}}$ $\sim$ 3.40--5.16 for the beam/strahl. The lower-to-upper quartile range of symmetric bi-self-similar core exponents are $s{\scriptstyle_{ec}}$ $\sim$ 2.00--2.04, and asymmetric bi-self-similar core exponents are $p{\scriptstyle_{ec}}$ $\sim$ 2.20--4.00 for the parallel exponent, and $q{\scriptstyle_{ec}}$ $\sim$ 2.00--2.46 for the perpendicular exponent. The nuanced details of the fit procedure and description of resulting data product are also presented. The statistics and detailed analysis of the results are presented in Paper II and Paper III of this three-part study.

We studied the disc of the unclassified B[e] star HD 50138, in order to explore its structure, and to find indications for the evolutionary status of this system, whether it is a young Herbig Be or a post-main-sequence star. Using high spatial resolution interferometric measurements from MIDI instrument (N-band) on the Very Large Telescope Interferometer, we analysed the disc size, the time-variability of the disc's thermal emission, and the spectral shape of the $10~\mu$m silicate feature. By fitting simple disc models, we determined the inclination and the mid-infrared size of the disc, confirming earlier results based on a lower number of observations. We searched for mid-infrared temporal variability of different regions of the disc, and concluded that its morphology is not experiencing significant changes over the observed epochs. We characterized the mid-infrared silicate feature by determining the feature amplitude and the $11.3/9.8~\mu$m flux ratio. The latter parameter is a good indicator of the grain size. The shape of the feature suggests the presence of crystalline silicate grains in the disc. The interferometric data revealed a strong radial trend in the mineralogy: while the disc's innermost region seems to be dominated by forsterite grains, at intermediate radii both forsterite and enstatite may be present. The outer disc may predominantly contain amorphous silicate particles. A comparison of the observed spectral shape with that of a sample of intermediate-mass stars (supergiants, Herbig Ae/Be stars, unclassified B[e] stars) implied that the evolutionary state of HD 50138 cannot be unambiguously decided from mid-IR spectroscopy.

In order to study the stellar population and possible substructures in the outskirts of Double Cluster $h$ and $\chi$ Persei, we investigate using the GAIA DR2 data a sky area of about 7.5 degrees in radius around the Double Cluster cores. We identify member stars using various criteria, including their kinematics (viz, proper motion), individual parallaxes, as well as photometric properties. A total of 2186 member stars in the parameter space were identified as members. Based on the spatial distribution of the member stars, we find an extended halo structure of $h$ and $\chi$ Persei, about 6 - 8 times larger than their core radii. We report the discovery of filamentary substructures extending to about 200 pc away from the Double Cluster. The tangential velocities of these distant substructures suggest that they are more likely to be the remnants of primordial structures, instead of a tidally disrupted stream from the cluster cores. Moreover, the internal kinematic analysis indicates that halo stars seems to be experiencing a dynamic stretching in the RA direction, while the impact of the core components is relatively negligible. This work also suggests that the physical scale and internal motions of young massive star clusters may be more complex than previously thought.

In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to study one common aspect of these objects: their often transient nature. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

The strongly isospin-dependent tensor force leads to short-range correlations (SRC) between neutron-proton (deuteron-like) pairs much stronger than those between proton-proton and neutron-neutron pairs. As a result of the short-range correlations, the single-nucleon momentum distribution develops a high-momentum tail above the Fermi surface. Because of the strongly isospin-dependent short-range correlations, in neutron-rich matter a higher fraction of protons will be depleted from its Fermi sea and populate above the Fermi surface compared to neutrons. This isospin-dependent nucleon momentum distribution may have effects on: (1) nucleon spectroscopic factors of rare isotopes, (2) the equation of state especially the density dependence of nuclear symmetry energy, (3) the coexistence of a proton-skin in momentum space and a neutron-skin in coordinate space (i.e., protons move much faster than neutrons near the surface of heavy nuclei). In this talk, we discuss these features and their possible experimental manifestations. As an example, SRC effects on the nuclear symmetry energy are discussed in detail using a modified Gogny-Hartree-Fock (GHF) energy density functional (EDF) encapsulating the SRC-induced high momentum tail (HMT) in the single-nucleon momentum distribution.

Based on an extended Skyrme interaction that includes the terms in relative momenta up to sixth order, corresponding to the so-called Skyrme pseudopotential up to next-to-next-to-next-to leading order (N3LO), we derive the expressions of Hamiltonian density and single nucleon potential under general non-equilibrium conditions which can be applied in transport model simulations of heavy-ion collisions induced by neutron-rich nuclei. While the conventional Skyrme interactions, which include the terms in relative momenta up to second order, predict an incorrect behavior as a function of energy for nucleon optical potential in nuclear matter, the present extended N3LO Skyrme interaction can give a nice description for the empirical nucleon optical potential. We also construct three interaction sets with different high-density behaviors of the symmetry energy, by fitting both the empirical nucleon optical potential up to energy of $1$ GeV and the empirical properties of isospin asymmetric nuclear matter. These extended N3LO Skyrme interactions will be useful in transport model simulations of heavy-ion collisions induced by neutron-rich nuclei at intermediate and high energies, and they can also be useful in nuclear structure studies within the mean-field model.

Since 2016 October, the active galaxy PKS 2247-131 has undergone a gamma-ray outburst, which we studied using data obtained with the Fermi Gamma-ray Space Telescope. The emission arises from a relativistic jet in PKS 2247-131, as an optical spectrum only shows a few weak absorption lines, typical of the BL Lacertae sub-class of the blazar class of active galactic nuclei. Here we report a ~34.5 day quasi-periodic oscillation (QPO) in the emission after the initial flux peak of the outburst. Compared to one-year time-scale QPOs, previously identified in blazars in Fermi energies, PKS 2247-131 exhibits the first clear case of a relatively short, month-like oscillation. We show that this QPO can be explained in terms of a helical structure in the jet, where the viewing angle to the dominant emission region in the jet undergoes periodic changes. The time scale of the QPO suggests the presence of binary supermassive black holes in PKS 2247-131.

The Space Technology 7 Disturbance Reduction System (ST7-DRS) is a NASA technology demonstration payload that operated from January 2016 through July of 2017 on the European Space Agency's LISA Pathfinder spacecraft. The joint goal of the NASA and ESA missions was to validate key technologies for a future space-based gravitational wave observatory targeting the source-rich milliHertz band. The two primary components of ST7-DRS are a micropropulsion system based on colloidal micro-Newton thrusters (CMNTs) and a control system that simultaneously controls the attitude and position of the spacecraft and the two free-flying test masses (TMs). This paper presents our main experimental results and summarizes the overall the performance of the CMNTs and control laws. We find that the CMNT performance to be consistent with pre-flight predictions, with a measured system thrust noise on the order of $100\,\textrm{nN}/\sqrt{\textrm{Hz}}$ in the $1\,\textrm{mHz}\leq f \leq 30\,\textrm{mHz}$ band. The control system maintained the TM-spacecraft separation with an RMS error of less than 2$\,$nm and a noise spectral density of less than $3\,\textrm{nm}/\sqrt{\textrm{Hz}}$ in the same band. Thruster calibration measurements yield thrust values consistent with the performance model and ground-based thrust-stand measurements, to within a few percent. We also report a differential acceleration noise between the two test masses with a spectral density of roughly $3\,\textrm{fm}/\textrm{s}^2/\sqrt{\textrm{Hz}}$ in the $1\,\textrm{mHz}\leq f \leq 30\,\textrm{mHz}$ band, slightly less than twice as large as the best performance reported with the baseline LISA Pathfinder configuration and below the current requirements for the Laser Interferometer Space Antenna (LISA) mission.

We present the distance priors from the finally released $Planck~ \textrm{TT,TE,EE}+\textrm{lowE}$ data in 2018. The uncertainties are around $40\%$ smaller than those from $Planck$ 2015 TT$+$lowP. In order to check the validity of these new distance priors, we adopt the distance priors to constrain the cosmological parameters in different dark energy models, including the $\Lambda$CDM model, the $w$CDM model and the CPL model, and conclude that the distance priors provide consistent constraints on the relevant cosmological parameters compared to those from the full $Planck$ 2018 data release.

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.

Radio detection of Ultra High Energetic Cosmic Rays and Neutrinos (UHECRv) which hit the Moon has been investigated in recent years. In preparation for near-future lunar science missions, we discuss technical requirements for radio experiments onboard lunar orbiters or on a lunar lander. We also develop an analysis of UHECRv aperture by including UHECv events occurring in the sub-layers of lunar regolith. It is verified that even using a single antenna onboard lunar orbiters or a few meters above the Moon's surface, dozens of lunar UHECRv events are detectable for one-year of observation at energy levels of $10^{18}$ eV to $10^{23}$ eV. Furthermore, it is shown that an antenna 3 meters above the Moon's surface could detect lower energy lunar UHECR events at the level of $10^{15}$ eV to $10^{18}$ eV which might not be detectable from lunar orbiters or ground-based observations.

We present a study of decimetric radio activity, using the first high time cadence (0.5 s) images from the ${\textit{Giant Meterwave Radio Telescope}}$ (${\textit{GMRT}}$) at 610 MHz, associated with ${\textit{GOES}}$ C1.4 and M1.0 class solar flares, and a coronal mass ejection (CME) onset that occurred on 20 June 2015. The high spatial resolution images from ${\textit{GMRT}}$ show a strong radio source during the C1.4 flare, located ~500$''$ away from the flaring site with no corresponding bright footpoints or coronal features nearby. In contrast, however, strong radio sources are found near the flaring site during the M1.0 flare and around the CME onset time. Weak radio sources, located near the flaring site, are also found during the maximum of the C1.4 flare activity, which show a temporal association with metric type III bursts identified by the Solar Broadband Radio Spectrometer at Yunnan Astronomical Observatory. Based on a multi-wavelength analysis and magnetic potential field source surface extrapolations, we suggest that the source electrons of ${\textit{GMRT}}$ radio sources and metric type III bursts were originated from a common electron acceleration site. We also show that the strong ${\textit{GMRT}}$ radio source is generated by a coherent emission process and its apparent location far from the flaring site is possibly due to the wave-ducting effect.

The temperatures of dust grains play important roles in the chemical evolution of molecular clouds. Unlike large grains, the temperature fluctuations of small grains induced by photons may be significant. Therefore, if the grain size distribution is included in astrochemical models, the temperatures of small dust grains may not be assumed to be constant. We simulate a full gas-grain reaction network with a set of dust grain radii using the classical MRN grain size distribution and include the temperature fluctuations of small dust grains. Monte Carlo method is used to simulate the real-time dust grain's temperature fluctuations which is caused by the external low energy photons and the internal cosmic ray induced secondary photons. The increase of dust grains radii as ice mantles accumulate on grain surfaces is also included in our models. We found that surface CO$_2$ abundances in models with grain size distribution and temperature fluctuations are more than one order of magnitude larger than those with single grain size. Small amounts of terrestrial complex organic molecules (COMs) can also be formed on small grains due to the temperature spikes induced by external low energy photons. However, cosmic ray induced secondary photons overheat small grains so that surface CO sublime and less radicals are formed on grains surfaces, thus the production of surface CO$_2$ and COMs decreases by about one order of magnitude. The overheating of small grains can be offset by grain growth so that the formation of surface CO$_2$ and COMs becomes more efficient.

Context. Protoplanetary disks show large diversity regarding their morphology and dust composition. With mid-infrared interferometry the thermal emission of disks can be spatially resolved, and the distribution and properties of the dust within can be studied. Aims. Our aim is to perform a statistical analysis on a large sample of 82 disks around low- and intermediate-mass young stars, based on mid-infrared interferometric observations. We intend to study the distribution of disk sizes, variability, and the silicate dust mineralogy. Methods. Archival mid-infrared interferometric data from the MIDI instrument on the VLTI are homogeneously reduced and calibrated. Geometric disk models are used to fit the observations to get spatial information about the disks. An automatic spectral decomposition pipeline is applied to analyze the shape of the silicate feature. Results. We present the resulting data products in the form of an atlas, containing N band correlated and total spectra, visibilities, and differential phases. The majority of our data can be well fitted with a continuous disk model, except for a few objects, where a gapped model gives a better match. From the mid-infrared size--luminosity relation we find that disks around T Tauri stars are generally colder and more extended with respect to the stellar luminosity than disks around Herbig Ae stars. We find that in the innermost part of the disks ($r \lesssim 1$~au) the silicate feature is generally weaker than in the outer parts, suggesting that in the inner parts the dust is substantially more processed. We analyze stellar multiplicity and find that in two systems (AB Aur and HD 72106) data suggest a new companion or asymmetric inner disk structure. We make predictions for the observability of our objects with the upcoming MATISSE instrument, supporting the practical preparations of future MATISSE observations of T Tauri stars.

Based on broadband spectral energy distribution (SEDs), we estimate the jet physical parameters of 1392 $\gamma$-ray-loud active galactic nuclei (AGNs), the largest sample so far. The (SED) jet power and magnetization parameter are derived for these AGNs. Out of these sources, the accretion disk luminosity of 232 sources and (extended) kinetic jet powers of 159 sources are compiled from archived papers. We find the following, (1) Flat-spectrum radio quasars (FSRQs) and BL Lacs are well separated by $\Gamma=-0.127\log L_{\gamma}+8.18$, in $\gamma$-ray luminosity versus photon index plane with a success rate of 88.6\%. (2) Most FSRQs present a (SED) jet power larger than the accretion power, which suggests that the relativistic jet-launching mechanism is dominated by the Blandford-Znajek process. This result confirms previous findings. (3) There is asignificant anticorrelation between jet magnetization and a ratio of the (SED) jet power to the (extended) kinetic jet power, which, for the first time, provides supporting evidence for the jet energy transportation theory: a high-magnetization jet may be more easily to transport energy to a large scale than a low-magnetization jet.

We compile the radio-optical-X-ray spectral energy distributions (SEDs) of 65 knots and 29 hotspots in 41 active galactic nucleus jets to examine their high energy radiation mechanisms. Their SEDs can be fitted with the single-zone leptonic models, except for the hotspot of Pictor A and six knots of 3C 273. The X-ray emission of one hotspot and 22 knots is well explained as synchrotron radiations under the equipartition condition; they usually have lower X-ray and radio luminosities than the others, which may be due to a lower beaming factor. An inverse Compton (IC) process is involved for explaining the X-ray emission of the other SEDs. Without considering the equipartition condition, their X-ray emission can be attributed to the synchrotron-self-Compton (SSC) process, but the derived jet power (P_jet) are not correlated with L_k and most of them are larger than L_k with more than three orders of magnitude, where L_k is the jet kinetic power estimated with their radio emission. Under the equipartition condition, the X-ray emission is well interpreted with the IC process to the cosmic microwave background photons (IC/CMB). In this scenario, the derived P_jet of knots and hotspots are correlated with and comparable to L_k. These results suggest that the IC/CMB model may be the promising interpretation of their X-ray emission. In addition, a tentative knot-hotspot sequence in the synchrotron peak-energy--peak-luminosity plane is observed, similar to the blazar sequence, which may be attributed to their different cooling mechanisms of electrons.

We report simultaneous observations at different energy bands in radio, optical, UV, X-rays and $\gamma$ rays of the flat-spectrum radio-quasar 3C 345. We built the light curve of the source at different frequencies from 2008, the beginning of the \textit{Fermi} all-sky survey, to 2016, using new data and public archives. In particular we obtained several optical spectra, to study the behavior of emission lines and the continuum in different activity states and to derive the black hole mass. 3C 345 showed two flaring episodes in 2009, which occurred simultaneously in $\gamma$ ray, optical/UV and X-rays, and were later followed in radio. The source shows an inverse Compton dominated spectral energy distribution, which moved from higher to lower frequencies from the high to the low state. The reverberation of emission lines during one outburst event allowed us to constrain the location of production of $\gamma$ rays very close to the broad-line region, and possibly in the jet-base. We report the observation of an increased accretion after the outburst, possibly induced by the decrease of magnetic field intensity with respect to the low state.

We present a new analysis of the LAMOST DR1 survey spectral database performed with the code SP_Ace, which provides the derived stellar parameters T$_{\rm eff}$, log (g), [Fe/H], and [$\alpha$/Fe] for 1,097,231 stellar objects. We tested the reliability of our results by comparing them to reference results from high spectral resolution surveys. The expected errors can be summarized as $\sim$120 K in T$_{\rm eff}$, $\sim$0.2 in log (g), $\sim$0.15 dex in [Fe/H], and $\sim$0.1 dex in [$\alpha$/Fe] for spectra with S/N$>$40, with some differences between dwarf and giant stars. SP_Ace provides error estimations consistent with the discrepancies observed between derived and reference parameters. Some systematic errors are identified and discussed. The resulting catalog is publicly available at the LAMOST and CDS websites.

In radio astronomy, the Ultra-Long Wavelengths (ULW) regime of longer than 10 m (frequencies below 30 MHz), remains the last virtually unexplored window of the celestial electromagnetic spectrum. The strength of the science case for extending radio astronomy into the ULW window is growing. However, the opaqueness of the Earth's ionosphere makes ULW observations by ground-based facilities practically impossible. Furthermore, the ULW spectrum is full of anthropogenic radio frequency interference (RFI). The only radical solution for both problems is in placing an ULW astronomy facility in space. We present a concept of a key element of a space-borne ULW array facility, an antenna that addresses radio astronomical specifications. A tripole-type antenna and amplifier are analysed as a solution for ULW implementation. A receiver system with a low power dissipation is discussed as well. The active antenna is optimized to operate at the noise level defined by the celestial emission in the frequency band 1 - 30 MHz. Field experiments with a prototype tripole antenna enabled estimates of the system noise temperature. They indicated that the proposed concept meets the requirements of a space-borne ULW array facility.

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

We analyze an asteroseismic binary system: KIC 7107778, a non-eclipsing, unresolved target, with solar-like oscillations in both components. We used \emph{Kepler} short cadence time series spanning nearly two years to obtain the power spectrum. Oscillation mode parameters were determined using Bayesian inference and a nested sampling Monte Carlo algorithm with the \texttt{DIAMONDS} package. The power profiles of the two components fully overlap, indicating their close similarity. We modelled the two stars with \texttt{MESA} and calculated oscillation frequencies with \texttt{GYRE}. Stellar fundamental parameters (mass, radius and age) were estimated by grid modelling with atmospheric parameters and the oscillation frequencies of $l=0,2$ modes as constraints. Most $l=1$ mixed modes were identified with models searched using a bisection method. Stellar parameters for the two sub-giant stars are: $M_A=1.42\pm0.06$ $M_{\astrosun}$, $M_B=1.39\pm0.03$ $M_{\astrosun}$, $R_A=2.93\pm0.05$ $R_{\astrosun}$, $R_B=2.76\pm0.04$ $R_{\astrosun}$, $t_A=3.32\pm0.54$ Gyr and $t_B=3.51\pm0.33$ Gyr. The mass difference of the system is $\sim$1\%. The results confirm their simultaneous birth and evolution, as is expected from binary formation. KIC 7107778 comprises almost identical twins, and is the first asteroseismic sub-giant binary to be detected.

We model simultaneous or quasi-simultaneous multi-band spectral energy distributions (SEDs) for a sample of 25 blazars that have radio core-shift measurements, where a one-zone leptonic model and Markov chain Monte Carlo technique are adopted. In the SED fitting for 23 low-synchrotron-peaked (LSP) blazars, the seed photons from broad line region (BLR) and molecular torus are considered respectively in the external Compton process. We find that the SED fitting with the seed photons from the torus are better than those utilizing BLR photons, which suggests that the $\gamma$-ray emitting region may be located outside the BLR. Assuming the magnetic field strength in the $\gamma$-ray emitting region as constrained from the SED fitting follows the magnetic field distribution as derived from the radio core-shift measurements (i.e., $B(R)\simeq B_{\rm 1 pc}(R/\rm 1pc)^{-1}$, $R$ is the distance from the central engine and $B_{\rm 1 pc}$ is the magnetic field strength at 1 pc), we further calculate the location of the $\gamma$-ray emitting region, $R_{\gamma}$, for these blazars. We find that $R_{\gamma}\sim2\times10^{4}R_{\rm S}\simeq10 R_{\rm BLR}$ ($R_{\rm S}$ is the Schwarzschild radius and $R_{\rm BLR}$ is the BLR size), where $R_{\rm BLR}$ is estimated from the broad line luminosities using the empirical correlations obtained using the reverberation mapping methods.

Various kinds of isovector nucleon effective masses are used in the literature to characterize the momentum/energy dependence of the nucleon symmetry potential or self-energy due to the space/time non-locality of the underlying isovector strong interaction in neutron-rich nucleonic matter. The multifaceted studies on nucleon isovector effective masses are multi-disciplinary in nature. Besides structures, masses and low-lying excited states of nuclei as well as nuclear reactions, studies of the isospin dependence of short-range correlations in nuclei from scatterings of high-energy electrons and protons on heavy nuclei also help understand nucleon effective masses especially the so-called E-mass in neutron-rich matter. A thorough understanding of all kinds of nucleon effective masses has multiple impacts on many interesting issues in both nuclear physics and astrophysics. We review some of the significant progresses made in recent years by the nuclear physics community in resolving some of the hotly debated and longstanding issues regarding nucleon effective masses especially in dense neutron-rich matter.

The Xuyi Schmidt Telescope Photometric Survey of the Galactic Anti-center (XSTPS-GAC) is a photometric sky survey that covers nearly 6 000 deg^2 towards Galactic anti-center in g r i bands. Half of its survey field locates on the Galactic Anti-center disk, which makes XSTPS-GAC highly suitable for searching new open clusters in the GAC region. In this paper, we report new open cluster candidates discovered in this survey, as well as properties of these open cluster candidates, such as age, distance and reddening, derived by isochrone fitting in the color-magnitude diagram (CMD). These open cluster candidates are stellar density peaks detected in the star density maps by applying the method from Koposov et al. (2008). Each candidate is inspected on its true color image composed from XSTPS-GAC three band images. Then its CMD is checked, in order to identify whether the central region stars have a clear isochrone-like trend differing from the background stars. The parameters derived from isochrone fitting for these candidates are mainly based on three band photometry of XSTPS-GAC. Meanwhile, when these new candidates are able to be seen clearly on 2MASS, their parameters are also derived based on 2MASS (J-H, J) CMD. Finally, there are 320 known open clusters rediscovered and 24 new open cluster candidates discovered in this work. Further more, the parameters of these new candidates, as well as another 11 known recovered open clusters, are properly determined for the first time.

If axion-like particles (ALPs) exist, photons can convert to ALPs on passage through regions containing magnetic fields. The magnetised intracluster medium of large galaxy clusters provides a region that is highly efficient at ALP-photon conversion. X-ray observations of Active Galactic Nuclei (AGNs) located within galaxy clusters can be used to search for and constrain ALPs, as photon-ALP conversion would lead to energy-dependent quasi-sinusoidal modulations in the X-ray spectrum of an AGN. We use Chandra observations of the central AGN of the Perseus Cluster, NGC1275, to place bounds on massive ALPs up to m_a ~ 10^-11 eV, extending previous work that used this dataset to constrain massless ALPs.

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.

Data from the NASA Magnetospheric Multiscale (MMS) mission are used to investigate asymmetric magnetic reconnection at the dayside boundary between the Earth's magnetosphere and the solar wind (the magnetopause). High-resolution measurements of plasmas, electric and magnetic fields, and waves are used to identify highly localized (~15 electron Debye lengths) standing wave structures with large electric-field amplitudes (up to 100 mV/m). These wave structures are associated with spatially oscillatory dissipation, which appears as alternatingly positive and negative values of J dot E (dissipation). For small guide magnetic fields the wave structures occur in the electron stagnation region at the magnetosphere edge of the EDR. For larger guide fields the structures also occur near the reconnection x-line. This difference is explained in terms of channels for the out-of-plane current (agyrotropic electrons at the stagnation point and guide-field-aligned electrons at the x-line).

A prediction of the steady-state reconnection electric field in asymmetric reconnection is obtained by maximizing the reconnection rate as a function of the opening angle made by the upstream magnetic field on the weak magnetic field (magnetosheath) side. The prediction is within a factor of two of the widely examined asymmetric reconnection model [Cassak and Shay, Phys. Plasmas 14, 102114, 2007] in the collisionless limit, and they scale the same over a wide parameter regime. The previous model had the effective aspect ratio of the diffusion region as a free parameter, which simulations and observations suggest is on the order of 0.1, but the present model has no free parameters. In conjunction with the symmetric case [Liu et al., Phys. Rev. Lett. 118, 085101, 2017], this work further suggests that this nearly universal number 0.1, essentially the normalized fast reconnection rate, is a geometrical factor arising from maximizing the reconnection rate within magnetohydrodynamic (MHD)-scale constraints.

Based on systematic analyses of several popular non-relativistic energy density functionals with mean-field approximation, we estimate the value of the fourth-order symmetry energy $E_{\text{sym,4}}(\rho)$ at nuclear normal density $\rho_0$ and its density dependence, and explore the correlation between $E_{\text{sym,4}}(\rho_0)$ and other macroscopic quantities of nuclear matter properties. We use the empirical values of some nuclear macroscopic quantities to construct model parameter sets by Monte Carlo method for the conventional Skyrme-Hartree-Fock (SHF) model, the extended Skyrme-Hartree-Fock (eSHF) model, the Gogny-Hartree-Fock (GHF) model, and the momentum-dependent interaction (MDI) model. The value of $E_{\text{sym,4}}(\rho_0)$ is estimated to be $1.02\pm0.49$ MeV for the SHF model, $1.02\pm0.50$ MeV for the eSHF model, $0.70\pm0.60$ MeV for the GHF model, and $0.74\pm0.63$ MeV for the MDI model. Moreover, our results indicate that the density dependence of $E_{\text{sym,4}}(\rho)$ is model dependent, especially at higher densities. Furthermore, we find that the $E_{\text{sym},4}(\rho_0)$ has strong positive (negative) correlation with isoscalar (isovector) nucleon effective mass $m_{s,0}^*$ ($m_{v,0}^*$) at $\rho_0$. In particular, for the SHF and eSHF models, the $E_{\text{sym,4}}(\rho)$ is completely determined by the isoscalar and isovector nucleon effective masses $m_{s}^*(\rho)$ and $m_{v}^*(\rho)$, and the analytical expression is given. In the mean-field models, the magnitude of $ E_{\text{sym,4}}(\rho_0)$ is generally less than $2$ MeV, and its density dependence depends on models, especially at higher densities. $ E_{\text{sym,4}}(\rho_0)$ is strongly correlated with $m_{s,0}^*$ and $m_{v,0}^*$.

There are two major theoretical issues for the star formation law (the relation between the surface densities of molecular gas and star formation rate on a galaxy scale): (i) At low metallicity, it is not obvious that star-forming regions are rich in H$_2$ because the H$_2$ formation rate depends on the dust abundance; and (ii) whether or not CO really traces H$_2$ is uncertain, especially at low metallicity. To clarify these issues, we use a hydrodynamic simulation of an isolated disc galaxy with a spatial resolution of a few tens parsecs. The evolution of dust abundance and grain size distribution is treated consistently with the metal enrichment and the physical state of the interstellar medium. We compute the H$_2$ and CO abundances using a subgrid post-processing model based on the dust abundance and the dissociating radiation field calculated in the simulation. We find that when the metallicity is $\lesssim 0.4$ Z$_\odot$ ($t<1$ Gyr), H$_2$ is not a good tracer of star formation rate because H$_2$-rich regions are limited to dense compact regions. At $Z\gtrsim 0.8$ Z$_\odot$, a tight star formation law is established for both H$_2$ and CO. At old ($t \sim 10$ Gyr) ages, we also find that adopting the so-called MRN grain size distribution with an appropriate dust-to-metal ratio over the entire disc gives reasonable estimates for the H$_2$ and CO abundances. For CO, improving the spatial resolution of the simulation is important while the H$_2$ abundance is not sensitive to sub-resolution structures at $Z\gtrsim 0.4$ Z$_\odot$.

The synchrotron peak of the X-ray bright High Energy Peaked Blazar (HBL) PKS 2155$-$304 occurs in the UV-EUV region and hence its X-ray emission (0.6--10 keV) lies mostly in the falling part of the synchrotron hump. We aim to study the X-ray emission of PKS 2155$-$304 during different intensity states in 2009$-$2014 using XMM$-$Newton satellite. We studied the spectral curvature of all of the observations to provide crucial information on the energy distribution of the non-thermal particles. Most of the observations show curvature or deviation from a single power-law and can be well modeled by a log parabola model. In some of the observations, we find spectral flattening after 6 keV. In order to find the possible origin of the X-ray excess, we built the Multi-band Spectral Energy distribution (SED). We find that the X-ray excess in PKS 2155--304 is difficult to fit in the one zone model but, could be easily reconciled in the spine/layer jet structure. The hard X-ray excess can be explained by the inverse Comptonization of the synchrotron photons (from the layer) by the spine electrons.

Within the nonlinear relativistic mean field (NL-RMF) model, we show that both the pressure of symmetric nuclear matter at supra-saturation densities and the maximum mass of neutron stars are sensitive to the skewness coefficient $J_0$ of symmetric nuclear matter. Using experimental constraints on the pressure of symmetric nuclear matter at supra-saturation densities from flow data in heavy ion collisions and the astrophysical observation of a large mass neutron star PSR J0348+0432, with the former favoring a smaller $J_0$ while the latter a larger $J_0$, we extract a constraint of $-494 \mathrm{MeV}\leq J_0\leq -10 \mathrm{MeV}$ based on the NL-RMF model. This constraint is compared with the results obtained in other analyses.

Finding the electromagnetic (EM) counterpart of binary compact star merger, especially the binary neutron star (BNS) merger, is critically important for gravitational wave (GW) astronomy, cosmology and fundamental physics. On Aug. 17, 2017, Advanced LIGO and \textit{Fermi}/GBM independently triggered the first BNS merger, GW170817, and its high energy EM counterpart, GRB 170817A, respectively, resulting in a global observation campaign covering gamma-ray, X-ray, UV, optical, IR, radio as well as neutrinos. The High Energy X-ray telescope (HE) onboard \textit{Insight}-HXMT (Hard X-ray Modulation Telescope) is the unique high-energy gamma-ray telescope that monitored the entire GW localization area and especially the optical counterpart (SSS17a/AT2017gfo) with very large collection area ($\sim$1000 cm$^2$) and microsecond time resolution in 0.2-5 MeV. In addition, \textit{Insight}-HXMT quickly implemented a Target of Opportunity (ToO) observation to scan the GW localization area for potential X-ray emission from the GW source. Although it did not detect any significant high energy (0.2-5 MeV) radiation from GW170817, its observation helped to confirm the unexpected weak and soft nature of GRB 170817A. Meanwhile, \textit{Insight}-HXMT/HE provides one of the most stringent constraints (~10$^{-7}$ to 10$^{-6}$ erg/cm$^2$/s) for both GRB170817A and any other possible precursor or extended emissions in 0.2-5 MeV, which help us to better understand the properties of EM radiation from this BNS merger. Therefore the observation of \textit{Insight}-HXMT constitutes an important chapter in the full context of multi-wavelength and multi-messenger observation of this historical GW event.

We update the constraints on the cosmological parameters by adopting the Planck data released in 2015 and Baryon Acoustic Oscillation (BAO) measurements including the new DR14 quasar sample measurement at redshift $z=1.52$, and we conclude that the based six-parameter $\Lambda$CDM model is preferred. Exploring some extensions to the $\Lambda$CDM models, we find that the equation of state of dark energy reads $w=-1.036\pm 0.056$ in the $w$CDM model, the effective relativistic degrees of freedom in the Universe is $N_\text{eff}=3.09_{-0.20}^{+0.18}$ in the $N_\text{eff}+\Lambda$CDM model and the spatial curvature parameter is $\Omega_k=(1.8\pm1.9)\times 10^{-3}$ in the $\Omega_k+\Lambda$CDM model at $68\%$ confidence level (C.L.), and the $95\%$ C.L. upper bounds on the sum of three active neutrinos masses are $\sum m_\nu<0.16$ eV for the normal hierarchy (NH) and $\sum m_\nu<0.19$ eV for the inverted hierarchy (IH) with $\Delta\chi^2\equiv \chi^2_\text{NH}-\chi^2_\text{IH}=-1.25$.

The symmetry energy characterizes the isospin dependent part of the equation of state of isospin asymmetric strong interaction matter and it plays a critical role in many issues of nuclear physics and astrophysics. In this talk, we briefly review the current status on the determination of the symmetry energy in nucleon (nuclear) and quark matter. For nuclear matter, while the subsaturation density behaviors of the symmetry energy are relatively well-determined and significant progress has been made on the symmetry energy around saturation density, the determination of the suprasaturation density behaviors of the symmetry energy remains a big challenge. For quark matter, which is expected to appear in dense matter at high baryon densities, we briefly review the recent work about the effects of quark matter symmetry energy on the properties of quark stars and the implication of possible existence of heavy quark stars on quark matter symmetry energy. The results indicate that the $u$ and $d$ quarks could feel very different interactions in isospin asymmetric quark matter, which may have important implications on the isospin effects of partonic dynamics in relativistic heavy-ion collisions.

We establish a relation between the equation of state (EOS) of nuclear matter and the fourth-order symmetry energy $a_{\rm{sym,4}}(A)$ of finite nuclei in a semi-empirical nuclear mass formula by self-consistently considering the bulk, surface and Coulomb contributions to the nuclear mass. Such a relation allows us to extract information on nuclear matter fourth-order symmetry energy $E_{\rm{sym,4}}(\rho_0)$ at normal nuclear density $\rho_0$ from analyzing nuclear mass data. Based on the recent precise extraction of $a_{\rm{sym,4}}(A)$ via the double difference of the "experimental" symmetry energy extracted from nuclear masses, for the first time, we estimate a value of $E_{\rm{sym,4}}(\rho_0) = 20.0\pm4.6$ MeV. Such a value of $E_{\rm{sym,4}}(\rho_0)$ is significantly larger than the predictions from mean-field models and thus suggests the importance of considering the effects of beyond the mean-field approximation in nuclear matter calculations.

For the first time, Kataoka \& Stawarz reported a clear detection of a hard X-ray excess, above $\gtrsim$20 keV, in the high-synchrotron-peaked BL Lac object Mrk 421. We find that this feature may not be produced by the low-energy part of the same electron population that produced the Fermi/LAT $\gamma$-ray. Because of that it is required that the power-law electron energy go down to $\gamma_{\rm min}\approx19$, which predicts a very strong radio emission (radio flux larger than the observed) even considering the synchrotron self-absorption effect. We investigate the possibility of this excess being produced from the spine/layer jet structure, which has been clearly detected in Mrk 421. We find that (1) similar to one-zone modeling, the spine emissions provide good modeling of the broadband spectral energy distribution, except for the hard X-ray excess; and (2) the hard X-ray excess can be well represented by the synchrotron photons (from the layer) being inverse Compton scattered by the spine electrons.

We present the evolution of the phase lags associated with the type-C QPO in GX 339--4 during the rising phase of the 2006/2007 outburst. We find that the phase lags at the QPO frequency are always positive (hard), and show very different behavior between QPOs with frequencies below and above $\sim1.7$ Hz: when the QPO frequency is below $\sim1.7$ Hz, the phase lags increase both with QPO frequency and energy, while when the QPO frequency is above $\sim1.7$ Hz, the phase lags remain more or less constant. When the QPO frequency is higher than $\sim1.7$ Hz, a broad feature is always present in the lag-energy spectra at around 6.5 keV, suggesting that the reflection component may have a significant contribution to the phase lags. Below $\sim1.7$ Hz, the QPO rms first decreases with energy and then turns to almost flat, while above $\sim1.7$ Hz, the QPO rms increases with energy. During the transition from the low-hard state to the hard-intermediate state, the second harmonic and subharmonic of this QPO appear in the power density spectra. The second-harmonic and subharmonic phase lags show very similar evolution with their centroid frequencies. However, the energy dependence of the second-harmonic and subharmonic phase lags are quite different. Our results suggest that, at different phases of the outburst, different mechanisms may be responsible for the phase lags of the QPO. We briefly discuss the possible scenarios for producing the lags.

The measurement of neutron star mass and radius is one of the most direct way to distinguish between various dense matter equations of state. The mass and radius of accreting neutron stars hosted in low mass X-ray binaries can be constrained by several methods, including photospheric radius expansion from type-I X-ray bursts and from quiescent spectra. In this paper, we apply for the first time these two methods simultaneously to constrain the mass and radius of Aql X-1, as a reliable distance estimation, high signal-to-noise ratio quiescent spectra from Chandra and XMM-Newton, and photospheric radius expansion bursts from RXTE are available. This is also used to verify the consistency between the two methods, and to narrow down the uncertainties of the neutron star mass and radius. It is found that the distance to Aql X-1 should be in the range of $4.0-5.75$ kpc, based on the overlapping confidence regions between photospheric radius expansion burst and quiescent spectra methods. In addition, we show that the mass and radius determined for the compact star in Aql X-1 are compatible with strange star equations of state and conventional neutron star models.

The observation and research of the solar radio emission have unique scientific values in solar and space physics and related space weather forecasting applications, since the observed spectral structures may carry important information about energetic electrons and underlying physical mechanisms. In this study, we present the design of a novel dynamic spectrograph that is installed at the Chashan solar radio station operated by Laboratory for Radio Technologies, Institute of Space Sciences at Shandong University. The spectrograph is characterized by the real-time storage of digitized radio intensity data in the time domain and its capability to perform off-line spectral analysis of the radio spectra. The analog signals received via antennas and amplified with a low-noise amplifier are converted into digital data at a speed reaching up to 32 k data points per millisecond. The digital data are then saved into a high-speed electronic disk for further off-line spectral analysis. Using different word length (1 k - 32 k) and time cadence (5 ms - 10 s) for the off-line fast Fourier transform analysis, we can obtain the dynamic spectrum of a radio burst with different (user-defined) temporal (5 ms - 10 s) and spectral (3 kHz ~ 320 kHz) resolution. This brings a great flexibility and convenience to data analysis of solar radio bursts, especially when some specific fine spectral structures are under study.

The Doppler tracking data of the Chang'e 3 lunar mission is used to constrain the stochastic background of gravitational wave in cosmology within the 1 mHz to 0.05 Hz frequency band. Our result improves on the upper bound on the energy density of the stochastic background of gravitational wave in the 0.02 Hz to 0.05 Hz band obtained by the Apollo missions, with the improvement reaching almost one order of magnitude at around 0.05 Hz. Detailed noise analysis of the Doppler tracking data is also presented, with the prospect that these noise sources will be mitigated in future Chinese deep space missions. A feasibility study is also undertaken to understand the scientific capability of the Chang'e 4 mission, due to be launched in 2018, in relation to the stochastic gravitational wave background around 0.01 Hz. The study indicates that the upper bound on the energy density may be further improved by another order of magnitude from the Chang'e 3 mission, which will fill the gap in the frequency band from 0.02 Hz to 0.1 Hz in the foreseeable future.

We systematically investigate the thermodynamic properties of homogeneous nuclear matter with light clusters at low densities and finite temperatures using a generalized nonlinear relativistic mean-field (gNL-RMF) model, in which light clusters up to $\alpha $ ($1 \le A \le 4$) are included as explicit degrees of freedom and treated as point-like particles with their interactions described by meson exchanges and the medium effects on the cluster binding energies are described by density- and temperature-dependent energy shifts with the parameters obtained by fitting the experimental cluster Mott densities. We find that the composition of low density nuclear matter with light clusters is essentially determined by the density- and temperature-dependence of the cluster binding energy shifts. Compared with the values of the conventional (second-order) symmetry energy, symmetry free energy and symmetry entropy, their fourth-order values are found to be significant at low densities ($\rho \sim 10^{-3}$ fm$^{-3}$) and low temperatures ($T \lesssim 3$ MeV), indicating the invalidity of the empirical parabolic law for the isospin asymmetry dependence of these nuclear matter properties. Our results indicate that in the density region of $\rho \gtrsim 0.02$ fm$^{-3}$, the clustering effects become insignificant and the nuclear matter is dominated by nucleon degree of freedom. In addition, we compare the gNL-RMF model predictions with the corresponding experimental data on the symmetry energy and symmetry free energy at low densities and finite temperatures extracted from heavy-ion collisions, and reasonable agreement is found.

Context. DG Tau is a low-mass pre-main sequence star, whose strongly accreting protoplanetary disk exhibits a so-far enigmatic behavior: its mid-infrared thermal emission is strongly time-variable, even turning the 10 $\mu$m silicate feature from emission to absorption temporarily. Aims. We look for the reason for the spectral variability at high spatial resolution and at multiple epochs. Methods. We study the temporal variability of the mid-infrared interferometric signal, observed with the VLTI/MIDI instrument at six epochs between 2011 and 2014. We fit a geometric disk model to the observed interferometric signal to obtain spatial information about the disk. We also model the mid-infrared spectra by template fitting to characterize the profile and time dependence of the silicate emission. We use physically motivated radiative transfer modeling to interpret the mid-infrared interferometric spectra. Results. The inner disk (r<1-3 au) spectra exhibit a 10 $\mu$m absorption feature related to amorphous silicate grains. The outer disk (r>1-3 au) spectra show a crystalline silicate feature in emission, similar to the spectra of comet Hale-Bopp. The striking difference between the inner and outer disk spectral feature is highly unusual among T Tauri stars. The mid-infrared variability is dominated by the outer disk. The strength of the silicate feature changed by more than a factor of two. Between 2011 and 2014 the half-light radius of the mid-infrared-emitting region decreased from 1.15 to 0.7 au. Conclusions. For the origin of the absorption we discuss four possible explanations: a cold obscuring envelope, an accretion heated inner disk, a temperature inversion on the disk surface and a misaligned inner geometry. The silicate emission in the outer disk can be explained by dusty material high above the disk plane, whose mass can change with time, possibly due to turbulence in the disk.

We present a new tool for colour-magnitude diagram (CMD) studies, $Powerful~CMD$. This tool is built on the basis of the advanced stellar population synthesis (ASPS) model, in which single stars, binary stars, rotating stars, and star formation history have been taken into account. Via $Powerful~CMD$, the distance modulus, colour excess, metallicity, age, binary fraction, rotating star fraction, and star formation history of star clusters can be determined simultaneously from observed CMDs. The new tool is tested via both simulated and real star clusters. Five parameters of clusters NGC6362, NGC6652, NGC6838 and M67 are determined and compared to other works. It is shown that this tool is useful for CMD studies, in particular for those with the data of the Hubble Space Telescope (HST). Moreover, we find that the inclusion of binaries in theoretical stellar population models may lead to smaller colour excess compared to the case of single star population models.

Based on a covariant coalescence model with a blast-wave-like parametrization for the phase-space configuration of constituent particles at freeze-out, we derive an approximate analytical formula for the yields of clusters produced in relativistic heavy-ion collisions. Compared to previous existing formulae, the present work additionally considers the contributions from the longitudinal dimension in momentum space, the relativistic corrections and the finite size effects of the produced clusters relative to the spatial distribution of constituent particles at freeze-out. The new analytical coalescence formula provides a useful tool to evaluate the yield of produced clusters, such as light nuclei from nucleon coalescence and hadrons from quark coalescence, in heavy-ion collisions. As a first application of the new analytical formula, we explore the strangeness population factor $S_3 = ^3_{\Lambda}$H/($^3$He$\times \Lambda/$p) based on nucleon/$\Lambda$ coalescence as well as the production of exotic hadrons based on quark coalescence, in central Pb+Pb collisions at $\sqrt{s_{NN}}=2.76$ TeV. The results are compared with the predictions from other models.

We performed a global study on the timing and spectral properties of type-B quasi-periodic oscillations (QPOs) in the outbursts of black hole X-ray binaries. The sample is built based on the observations of {\it Rossi X-ray Timing Explorer}, via searching in the literature in RXTE era for all the identified type-B QPOs. To enlarge the sample, we also investigated some type-B QPOs that are reported but not yet fully identified. Regarding to the time lag and hard/soft flux ratio, we found that the sources with type-B QPOs behave in two subgroups. In one subgroup, type-B QPO shows a hard time lag that firstly decreases and then reverse into a soft time lag along with softening of the energy spectrum. In the other subgroup, type-B QPOs distribute only in a small region with hard time lag and relatively soft hardness. These findings may be understood with a diversity of the homogeneity showing up for the hot inner flow of different sources. We confirm the universality of a positive relation between the type-B QPO frequency and the hard component luminosity in different sources. We explain the results by considering that the type-B QPO photons are produced in the inner accretion flow around the central black hole, under a local Eddington limit. Using this relationship, we derived an mass estimation of $9.3-27.1$ $\rm{M_{\odot}}$ for the black hole in H 1743-322.

We analyze the relation between the emission radii of twin kilohertz quasi-periodic oscillations (kHz QPOs) and the co-rotation radii of the 12 neutron star low mass X-ray binaries (NS-LMXBs) which are simultaneously detected with the twin kHz QPOs and NS spins. We find that the average co-rotation radius of these sources is r_co about 32 km, and all the emission positions of twin kHz QPOs lie inside the corotation radii, indicating that the twin kHz QPOs are formed in the spin-up process. It is noticed that the upper frequency of twin kHz QPOs is higher than NS spin frequency by > 10%, which may account for a critical velocity difference between the Keplerian motion of accretion matter and NS spin that is corresponding to the production of twin kHz QPOs. In addition, we also find that about 83% of twin kHz QPOs cluster around the radius range of 15-20 km, which may be affected by the hard surface or the local strong magnetic field of NS. As a special case, SAX J1808.4-3658 shows the larger emission radii of twin kHz QPOs of r about 21-24 km, which may be due to its low accretion rate or small measured NS mass (< 1.4 solar mass).

We use Planck data released in 2015 to constrain the abundance of primordial black holes (PBHs) in dark matter in two different reionization models (one is the instantaneous reionization and the other is the asymmetric reionization), and significantly improve the upper limits on the abundance of PBHs from WMAP 3-year data by around two orders of magnitude. Furthermore, these new limits imply that the event rates of mergers of PBH binaries (Gpc$^{-3}$ yr$^{-1}$) are less than $0.002$ for $M_\text{pbh}=30M_\odot$, $5$ for $M_\text{pbh}=10M_\odot$ and $2000$ for $M_\text{pbh}=2M_\odot$ at $95\%$ confidence level (C.L.), and thus GW150914 seems very unlikely produced by the merger of a PBH binary.

Recent progress in theory, experiment and observation challenges the mean field models using the conventional Skyrme interaction, suggesting that the extension of the conventional Skyrme interaction is necessary. In this work, by fitting the experimental data of a number of finite nuclei together with a few additional constraints on nuclear matter using the simulated annealing method, we construct three Skyrme interaction parameter sets, namely, eMSL07, eMSL08 and eMSL09, based on an extended Skyrme interaction which includes additional momentum and density dependent two-body forces to effectively simulate the momentum dependence of the three-body force. The three new interactions can reasonably describe the ground-state properties and the isoscalar giant monopole resonance energies of various spherical nuclei used in the fit as well as the ground-state properties of many other spherical nuclei, nicely conform to the current knowledge on the equation of state of asymmetric nuclear matter, eliminate the notorious unphysical instabilities of symmetric nuclear matter and pure neutron matter up to a very high density of $1.2$ fm$^{-3}$, and simultaneously support heavier neutron stars with mass larger than two times solar mass. One important difference of the three new interactions is about the prediction of the symmetry energy at supra-saturation densities, and these new interactions are thus potentially useful for the determination of the largely uncertain high-density symmetry energy in future. In addition, a comparison is made for the predictions of nuclear matter, finite nuclei and neutron stars with the three new interactions versus those with three typical interactions BSk22, BSk24 and BSk26 from Brussels group.

Neutron-skins in coordinate space and proton-skins in momentum space are predicted to coexist in heavy nuclei and their correlation is governed by Liouville's theorem and Heisenberg's uncertainty principle. An analysis of their correlation within a further extended Thomas-Fermi (ETF$^+$) approximation incorporating effects of nucleon short-range correlations reveals generally that protons move faster than neutrons in neutron-skins of heavy nuclei.

Nuclear fusion reactions are the most important processes in nature to power stars and produce new elements, and lie at the center of the understanding of nucleosynthesis in the universe. It is critically important to study the reactions in full plasma environments that are close to true astrophysical conditions. By using laser-driven counter-streaming collisionless plasmas, we studied the fusion D$+$D$\rightarrow n +^3$He in a Gamow-like window around 27 keV. The results show that astrophysical nuclear reaction yield can be modulated significantly by the self-generated electromagnetic fields and the collective motion of the plasma. This plasma-version mini-collider may provide a novel tool for studies of astrophysics-interested nuclear reactions in plasma with tunable energies in earth-based laboratories.

In this talk, we first briefly review the isospin dependence of the total nucleon effective mass $M^{\ast}_{J}$ inferred from analyzing nucleon-nucleus scattering data within an isospin dependent non-relativistic optical potential model, and the isospin dependence of the nucleon E-mass $M^{\ast,\rm{E}}_{J}$ obtained from applying the Migdal-Luttinger theorem to a phenomenological single-nucleon momentum distribution in nuclei constrained by recent electron-nucleus scattering experiments. Combining information about the isospin dependence of both the nucleon total effective mass and E-mass, we then infer the isospin dependence of nucleon k-mass using the well-known relation $M^{\ast}_{J}=M^{\ast,\rm{E}}_{J}\cdot M^{\ast,\rm{k}}_{J}$. Implications of the results on the nucleon mean free path (MFP) in neutron-rich matter are discussed.

III Zw 2 is the prototype of radio-intermediate quasars. Although there is the evidence of possessing strong jet, significant $\gamma$-ray emission has not been reported before. In this work, we carry out a detailed analysis of the latest {\it Fermi}-LAT {\it Pass} 8 data. No significant $\gamma$-ray signal has been detected in the time-averaged 7-year {\it Fermi}-LAT data of III Zw 2, we however have identified two distinct $\gamma$-ray flares with isotropic luminosities of $\sim 10^{45}$ erg $\rm s^{-1}$. Multiwavelength data analysis (including also the optical photometric observations from Yunnan Observatories) are presented and the main finding is the simultaneous optical and $\gamma$-ray flares of III Zw 2 appearing in Nov. 2009. Violent $\gamma$-ray variability with doubling timescale of 2.5 hours was detected in another $\gamma$-ray flare in May 2010, for which the 3-hour $\gamma$-ray peak flux is $\sim 250$ times of the average flux in 7 years. Rather similar behaviors are observed in blazars and the blazar model can reasonably reproduce the spectral energy distribution of III Zw 2 in a wide energy range, strongly suggesting that its central engine resembles that of blazars. In view of its core which shares radio similarities with young radio sources together with the weak extended radio lobe emission, we suggest that III Zw 2 harbors a recurrent activity core and serves as a valuable target for investigating the fuelling and triggering of the activity in radio loud AGNs.

We demonstrate that the observation of old strange quark stars (SQSs) can set important limits on the scattering cross sections $\sigma_q$ between the light quarks and the non-interacting scalar dark matter (DM). By analyzing a set of 1403 of solitary pulsarlike compact stars in the Milky Way, we find the old solitary pulsar PSR J1801-0857D can set the most stringent upper limits on $\sigma_q$ or the DM-proton scattering cross sections $\sigma_p$. By converting $\sigma_q$ into $\sigma_p$ based on effective operator analyses, we show the resulting $\sigma_p$ limit by assuming PSR J1801-0857D to be a SQS could be comparable with that of the current direct detection experiments but much weaker (by several orders of magnitude) than that obtained by assuming PSR J1801-0857D to be a neutron star (NS), which requires an extremely small $\sigma_p$ far beyond the limits of direct detection experiments. Our findings imply that the old pulsars are favored to be SQSs rather than NSs if the scalar DM were observed by future terrestrial experiments.

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

We show that the measured yield ratio $^3_{\Lambda}$H/$^3$He ($^3_{\overline{\Lambda}}\overline{\text{H}}$/$^3\overline{\text{He}}$) in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV and in Pb+Pb collisions at $\sqrt{s_{NN}}=2.76$ TeV can be understood within a covariant coalescence model if (anti-)$\Lambda$ particles freeze out earlier than (anti-)nucleons but their relative freezeout time is closer at $\sqrt{s_{NN}}=2.76$ TeV than at $\sqrt{s_{NN}}=200$ GeV. The earlier (anti-)$\Lambda$ freezeout can significantly enhance the yield of (anti)hypernucleus $^4_{\Lambda}$H ($^4_{\overline{\Lambda}}\overline{\text{H}}$), leading to that $^4_{\overline{\Lambda}}\overline{\text{H}}$ has a comparable abundance with $^4\overline{\text{He}}$ and thus provides an easily measured antimatter candidate heavier than $^4\overline{\text{He}}$. The future measurement on $^4_{\Lambda}$H ($^4_{\overline{\Lambda}}\overline{\text{H}}$) would be very useful to understand the (anti-)$\Lambda$ freezeout dynamics and the production mechanism of (anti)hypernuclei in relativistic heavy-ion collisions.

PKS 1424+240 is a distant very high energy gamma-ray BL Lac object with redshift $z=0.601$. It was found that pure synchrotron self-Compton (SSC) process normally need extreme input parameters (e.g., very low magnetic field intensity and extraordinarily large Doppler factor) to explain its multi-wavelength spectral energy distributions (SEDs). To avoid the extreme model parameters, different models have been proposed (e.g., two-zone SSC model or lepto-hadronic model). In this work, we employ the traditional one-zone leptonic model after including a weak external Compton component to re-explore the simultaneous multi-wavelength SEDs of PKS 1424+240 in both high (2009) and low (2013) states. We find that the input parameters of magnetic field and Doppler factor are roughly consistent with those of other BL Lacs if a weak external photon field from either broad line region (BLR) or the dust torus. However, the required energy density of seed photons from BLR or torus is about 3 orders of magnitude less than that constrained in luminous quasars (e.g., flat-spectrum radio quasars, FSRQs). This result suggests that the BLR/torus in BL Lacs is much weaker than that of luminous FSRQs (but not fully disappear), and the inverse-Compton of external photons from BLR/torus may still play a role even in high synchrotron peaked blazars.

Stellar rotation, age spread and binary stars are thought to be three most possible causes of the peculiar color-magnitude diagrams (CMDs) of some star clusters, which exhibit extended main-sequence turn-offs (eMSTOs). It is far from getting a clear answer. This paper studies the effects of three above causes on the CMDs of star clusters systematically. A rapid stellar evolutionary code and a recently published database of rotational effects of single stars have been used, via an advanced stellar population synthesis technique. As a result, we find a consistent result for rotation to recent works, which suggests that rotation is able to explain, at least partially, the eMSTOs of clusters, if clusters are not too old ($<$ 2.0\,Gyr). In addition, an age spread of 200 to 500\,Myr reproduces extended turn-offs for all clusters younger than 2.5\,Gyr, in particular, for those younger than 2.2\,Gyr. Age spread also results in extended red clumps (eRCs) for clusters younger than 0.5\,Gyr. The younger the clusters, the clearer the eRC structures. Moreover, it is shown that binaries (including interactive binaries) affect the spread of MSTO slightly for old clusters, but they can contribute to the eMSTOs of clusters younger than 0.5\,Gyr. Our result suggests a possible way to disentangle the roles of stellar rotation and age spread, i.e., checking the existence of CMDs with both eMSTO and eRC in clusters younger than 0.5\,Gyr.

Using a spectroscopically confirmed sample of M-giants, M-dwarfs and quasars from the LAMOST survey, we assess how well WISE $\&$ 2MASS color-cuts can be used to select M-giant stars. The WISE bands are very efficient at separating M-giants from M-dwarfs and we present a simple classification that can produce a clean and relatively complete sample of M-giants. We derive a new photometric relation to estimate the metallicity for M-giants, calibrated using data from the APOGEE survey. We find a strong correlation between the $(W1-W2)$ color and $\rm [M/H]$, where almost all of the scatter is due to photometric uncertainties. We show that previous photometric distance relations, which are mostly based on stellar models, may be biased and devise a new empirical distance relation, investigating trends with metallicity and star formation history. Given these relations, we investigate the properties of M-giants in the Sagittarius stream. The offset in the orbital plane between the leading and trailing tails is reproduced and, by identifying distant M-giants in the direction of the Galactic anti-center, we confirm that the previously detected debris in the outer halo is the apocenter of the trailing tail. We also find tentative evidence supporting an existing overdensity near the leading tail in the Northern Galactic hemisphere, possibly an extension to the trailing tail (so-called Branch C). We have measured the metallicity distribution along the stream, finding a clear metallicity offset between the leading and trailing tails, in agreement with models for the stream formation. We include an online table of M-giants to facilitate further studies.

Fast electron beams (FEBs) are common products of solar active phenomena. Solar radio bursts are an important diagnostic tool in the understanding of FEBs as well as the solar plasma environment in which they are propagating along solar magnetic fields. In particular, the evolutions of the energy spectrum and velocity distribution of FEBs due to the interaction with the ambient plasma and field when propagating can significantly influence the efficiency and property of their emissions. In this paper, we discuss some possible evolutions of the energy spectrum and velocity distribution of FEBs due to the energy loss processes and the pitch-angle effect caused by the magnetic field inhomogeneity, and analyze the effects of these evolutions on electron cyclotron maser (ECM) emission, which is one of the most important mechanisms of producing solar radio bursts by FEBs. The results show that the growth rates all decrease with the energy loss factor $Q$, but increase with the magnetic mirror ratio $\sigma$ as well as with the steepness index $\delta$. Moreover, the evolution of FEBs also can significantly influence the fastest growing mode and the fastest growing phase angle. This leads to the change of the polarization sense of ECM emission. In particular, our results also reveal that FEB that undergoes different evolution processes will generate different types of ECM emission. We believe the present results to be very helpful on more comprehensive understanding of dynamic spectra of solar radio bursts.

Based on mean field calculations with Skyrme interactions, we extract a constraint on the isovector effective mass in nuclear matter at saturation density $\rho_0$, i.e., $m_{v}^{\ast}(\rho_0)=(0.77\pm0.03) m$ by combining the experimental data of the centroid energy of the isovector giant dipole resonance (IVGDR) and the electric dipole polarizability $\alpha_{\mathrm{D}}$ in $^{208}$Pb. Meanwhile, the isoscalar effective mass at $\rho_0$ is determined to be $m_{s}^{\ast}(\rho_0)=(0.91\pm0.05) m$ by analyzing the measured excitation energy of the isoscalar giant quadrupole resonance (ISGQR) in $^{208}$Pb. From the constrained $m_{s}^{\ast}(\rho_0)$ and $m_{v}^{\ast}(\rho_0)$, we obtain the isospin splitting of nucleon effective mass in asymmetric nuclear matter of isospin asymmetry $\delta$ at $\rho_0$ as $[m_n^{\ast}(\rho_0,\delta)-m_p^{\ast}(\rho_0,\delta)]/m = \Delta m^*_1(\rho_0) \delta + O(\delta^3)$ with the linear isospin splitting coefficient $\Delta m^*_1(\rho_0) = 0.33\pm0.16$. We notice that using the recently corrected data on the $\alpha_{\mathrm{D}}$ in $^{208}$Pb with the contribution of the quasideuteron effect subtracted slightly enhances the isovector effective mass to $m_{v}^{\ast}(\rho_0)=(0.80\pm0.03) m$ and reduces the linear isospin splitting coefficient to $\Delta m^*_1(\rho_0) = 0.27\pm0.15$. Furthermore, the constraints on $m_{v}^{\ast}(\rho)$, $m_{s}^{\ast}(\rho)$ and $\Delta m^*_1(\rho)$ at other densities are obtained from the similar analyses and we find that the $\Delta m^*_1(\rho)$ increases with the density.

Energetic electrons with power-law spectrum are most commonly observed in astrophysics. This paper investigates electron cyclotron maser emission (ECME) from the power-law electrons, in which strong pitch-angle anisotropy is emphasized. The electron distribution function proposed in this paper can describe various types of pitch-angle anisotropy. Results show that the emission properties of ECME, including radiation growth, propagation, and frequency properties, depend considerably on the types of electron pitch-angle anisotropy, and different wave modes show different dependences on the pitch angle of electrons. In particular, the maximum growth rate of X2 mode rapidly decreases with respect to the electron pitch-angle cosine $\mu_0$ at which the electron distribution peaks, while the growth rates for other modes (X1, O1, O2) initially increase before decreasing as $\mu_0$ increases. Moreover, the O mode as well as the X mode can be the fastest growth mode, in terms of not only the plasma parameter but also the type of electron pitch-angle distribution. This result presents a significant extension of the recent researches on ECME driven by the lower-energy cutoff of power-law electrons, in which the X mode is generally the fastest growth mode.

A quark model with running coupling and running strange quark mass, which is thermodynamically self-consistent at both high and lower densities, is presented and applied to study properties of strange quark matter and structure of compact stars. An additional term to the thermodynamic potential density is determined by meeting the fundamental differential equation of thermodynamics. It plays an important role in comparatively lower density and ignorable at extremely high density, acting as a chemical-potential dependent bag constant. In this thermodynamically enhanced perturbative QCD model, strange quark matter still has the possibility of being absolutely stable, while the pure quark star has a sharp surface with a maximum mass as large as about 2 times the solar mass and a maximum radius of about 11 kilometers.

A small fraction($<10\%$) of SDSS main sample galaxies(MGs) have not been targeted with spectroscopy due to the the fiber collision effect. These galaxies have been compiled into the input catalog of the LAMOST extra-galactic survey and named as the complementary galaxy sample. In this paper, we introduce the project and the status of the spectroscopies of the complementary galaxies in the first two years of the LAMOST spectral survey(till Sep. of 2014). Moreover, we present a sample of 1,102 galaxy pairs identified from the LAMOST complementary galaxies and SDSS MGs, which are defined as that the two members have a projected distance smaller than 100 kpc and the recessional velocity difference smaller than 500 $\rm kms^{-1}$. Compared with the SDSS only selected galaxy pairs, the LAMOST-SDSS pairs take the advantages of not being biased toward large separations and therefor play as a useful supplement to the statistical studies of galaxy interaction and galaxy merging.

Studies of pre-transitional disks, with a gap region between the inner infrared-emitting region and the outer disk, are important to improving our understanding of disk evolution and planet formation. Previous infrared interferometric observations have shown hints of a gap region in the protoplanetary disk around the Herbig Ae star HD~144432. We study the dust distribution around this star with two-dimensional radiative transfer modeling. We compare the model predictions obtained via the Monte-Carlo radiative transfer code RADMC-3D with infrared interferometric observations and the {\SED} of HD~144432. The best-fit model that we found consists of an inner optically thin component at $0.21\enDash0.32~\AU$ and an optically thick outer disk at $1.4\enDash10~\AU$. We also found an alternative model in which the inner sub-AU region consists of an optically thin and an optically thick component. Our modeling suggests an optically thin component exists in the inner sub-AU region, although an optically thick component may coexist in the same region. Our modeling also suggests a gap-like discontinuity in the disk of HD~144432.

TianQin is a proposal for a space-borne detector of gravitational waves in the millihertz frequencies. The experiment relies on a constellation of three drag-free spacecraft orbiting the Earth. Inter-spacecraft laser interferometry is used to monitor the distances between the test masses. The experiment is designed to be capable of detecting a signal with high confidence from a single source of gravitational waves within a few months of observing time. We describe the preliminary mission concept for TianQin, including the candidate source and experimental designs. We present estimates for the major constituents of the experiment's error budget and discuss the project's overall feasibility. Given the current level of technology readiness, we expect TianQin to be flown in the second half of the next decade.

We use near-infrared observations obtained as part of the {\sl Visible and Infrared Survey Telescope for Astronomy} (VISTA) Survey of the Magellanic Clouds (VMC), as well as two complementary {\sl Hubble Space Telescope} ({\sl HST}) data sets, to study the luminosity and mass functions as a function of clustercentric radius of the main-sequence stars in the Galactic globular cluster 47 Tucanae. The {\sl HST} observations indicate a relative deficit in the numbers of faint stars in the central region of the cluster compared with its periphery, for $18.75\leq m_{\rm F606W}\leq 20.9$ mag (corresponding to a stellar mass range of $0.55<m_\ast/{M_\odot}<0.73$). The stellar number counts at $6.7'$ from the cluster core show a deficit for $17.62\leq m_{\rm F606W}\leq 19.7$ mag (i.e., $0.65<m_\ast/{M_\odot}<0.82$), which is consistent with expectations from mass segregation. The VMC-based stellar mass functions exhibit power-law shapes for masses in the range $0.55<m_\ast/{M_\odot}< 0.82$. These power laws are characterized by an almost constant slope, $\alpha$. The radial distribution of the power-law slopes $\alpha$ thus shows evidence of the importance of both mass segregation and tidal stripping, for both the first- and second-generation stars in 47 Tuc.

One century after its formulation, Einstein's general relativity has made remarkable predictions and turned out to be compatible with all experimental tests. Most of these tests probe the theory in the weak-field regime, and there are theoretical and experimental reasons to believe that general relativity should be modified when gravitational fields are strong and spacetime curvature is large. The best astrophysical laboratories to probe strong-field gravity are black holes and neutron stars, whether isolated or in binary systems. We review the motivations to consider extensions of general relativity. We present a (necessarily incomplete) catalog of modified theories of gravity for which strong-field predictions have been computed and contrasted to Einstein's theory, and we summarize our current understanding of the structure and dynamics of compact objects in these theories. We discuss current bounds on modified gravity from binary pulsar and cosmological observations, and we highlight the potential of future gravitational wave measurements to inform us on the behavior of gravity in the strong-field regime.

Combining the covariant coalescence model and a blast-wave-like analytical parametrization for (anti-)nucleon phase-space freezeout configuration, we explore light (anti-)nucleus production in central Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV. Using the nucleon freezeout configuration (denoted by FO1) determined from the measured spectra of protons (p), deutrons (d) and $^{3}$He, we find the predicted yield of $^{4}$He is significantly smaller than the experimental data. We show this disagreement can be removed by using a nucleon freezeout configuration (denoted by FO2) in which the nucleons are assumed to freeze out earlier than those in FO1 to effectively consider the effect of large binding energy value of $^{4}$He. Assuming the binding energy effect also exists for the production of $^5\text{Li}$, $^5\overline{\text{Li}}$, $^6\text{Li}$ and $^6\overline{\text{Li}}$ due to their similar binding energy values as $^{4}$He, we find the yields of these heavier (anti-)nuclei can be enhanced by a factor of about one order, implying that although the stable (anti-)$^6$Li nucleus is unlikely to be observed, the unstable (anti-)$^5$Li nucleus could be produced in observable abundance in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV where it may be identified through the p-$^4\text{He}$ ($\overline{\text{p}}$-$^4\overline{\text{He}}$) invariant mass spectrum. The future experimental measurement on (anti-)$^5\text{Li}$ would be very useful to understand the production mechanism of heavier antimatter.

Exploring nucleon drip lines and astrophysical rapid neutron capture process (r-process) paths in the nuclear landscape is extremely challenging in nuclear physics and astrophysics. While various models predict similar proton drip line, their predictions for neutron drip line and the r-process paths involving heavy neutron-rich nuclei exhibit a significant variation which hampers our accurate understanding of the r-process nucleosynthesis mechanism. Using microscopic density functional theory with a representative set of non-relativistic and relativistic interactions, we demonstrate for the first time that this variation is mainly due to the uncertainty of nuclear matter symmetry energy $E_{\rm{sym}}(\rho_{\rm{sc}})$ at the subsaturation cross density $\rho_{\rm{sc}}=0.11/0.16\times\rho_0$ ($\rho_0$ is saturation density), which reflects the symmetry energy of heavy nuclei. Using the recent accurate constraint on $E_{\rm{sym}}(\rho_{\rm{sc}})$ from the binding energy difference of heavy isotope pairs, we obtain quite precise predictions for the location of the neutron drip line, the r-process paths and the number of bound nuclei in the nuclear landscape. Our results have important implications on extrapolating the properties of unknown neutron-rich rare isotopes from the data on known nuclei.

We infer the emission positions of twin kilohertz quasi-periodic oscillations (kHz QPOs) in neutron star low mass X-ray binaries (NS-LMXBs) based on the Alfven wave oscillation model (AWOM). For most sources, the emission radii of kHz QPOs cluster around a region of 16-19 km with the assumed NS radii of 15 km. Cir X-1 has the larger emission radii of 23-38 km than those of the other sources, which may be ascribed to its large magnetosphere-disk radius or strong NS surface magnetic field. SAX J1808.4-3658 is also a particular source with the relative large emission radii of kHz QPOs of 20 - 23 km, which may be due to its large inferred NS radius of 18 - 19 km. The emission radii of kHz QPOs for all the sources are larger than the NS radii, and the possible explanations of which are presented. The similarity of the emission radii of kHz QPOs (16-19 km) for both the low/high luminosity Atoll/Z sources is found, which indicates that both sources share the similar magnetosphere- disk radii.

It is currently a big challenge to accurately determine the symmetry energy $E_{\text{sym}}(\rho)$ and the pure neutron matter equation of state $E_{\text{PNM}}(\rho)$, even their values around saturation density $\rho_0 $. We find that the electric dipole polarizability $\alpha _ {\text{D}}$ in $^{208}$Pb can be determined uniquely by the magnitude of the $E_{\text{sym}}(\rho)$ or almost equivalently the $E_{\text{PNM}}(\rho)$ at subsaturation densities around $\rho_0/3 $, shedding a light upon the genuine correlation between the $\alpha _ {\text{D}}$ and the $E_{\text{sym}}(\rho)$. By analyzing the experimental data of the $\alpha _ {\text{D}}$ in $^{208}$Pb from RCNP using a number of non-relativistic and relativistic mean-field models, we obtain very stringent constraints on $E_{\text{sym}}(\rho)$ and $E_{\text{PNM}}(\rho)$ around $\rho_0/3 $. The obtained constraints are found to be in good agreement with the results extracted in other analyses. In particular, our results provide for the first time the experimental constraints on $E_{\text{PNM}}(\rho)$ around $\rho_0/3 $, which are in harmony with the recent determination of $E_{\text{PNM}}(\rho)$ from microscopic theoretical studies and potentially useful in constraining the largely uncertain many-nucleon interactions in microscopic calculations of neutron matter.

A $^{238}$U projectile beam was used to create cadmium isotopes via abrasion-fission at 410 MeV/u in a beryllium target at the entrance of the in-flight separator FRS at GSI. The fission fragments were separated with the FRS and injected into the isochronous storage ring ESR for mass measurements. The Isochronous Mass Spectrometry (IMS) was performed under two different experimental conditions, with and without B$\rho$-tagging at the dispersive central focal plane of the FRS. In the experiment with B$\rho$-tagging the magnetic rigidity of the injected fragments was determined by an accuracy of $2\times 10^{-4}$. A new method of data analysis, using a correlation matrix for the combined data set from both experiments, has provided mass values for 25 different isotopes for the first time. The high selectivity and sensitivity of the experiment and analysis has given access even to rare isotopes detected with a few atoms per week. In this letter we present for the $^{129,130,131}$Cd isotopes mass values directly measured for the first time. The Cd results clearly show a very pronounced shell effect at $N=82$ which is in agreement with the conclusion from $\gamma$-ray spectroscopy of $^{130}$Cd and confirms the assumptions of modern shell-model calculations.

We explore the systematics of the density dependence of nuclear matter symmetry energy in the ambit of microscopic calculations with various energy density functionals, and find that the symmetry energy from subsaturation density to supra-saturation density can be well determined by three characteristic parameters of the symmetry energy at saturation density $\rho_0 $, i.e., the magnitude $E_{\text{sym}}({\rho_0 })$, the density slope $L$ and the density curvature $K_{\text{sym}}$. This finding opens a new window to constrain the supra-saturation density behavior of the symmetry energy from its (sub-)saturation density behavior. In particular, we obtain $L=46.7 \pm 12.8$ MeV and $K_{\text{sym}}=-166.9 \pm 168.3$ MeV as well as $E_{\text{sym}}({2\rho _{0}}) \approx 40.2 \pm 12.8$ MeV and $L({2\rho _{0}}) \approx 8.9 \pm 108.7$ MeV based on the present knowledge of $E_{\text{sym}}({\rho_{0}}) = 32.5 \pm 0.5$ MeV, $E_{\text{sym}}({\rho_c}) = 26.65 \pm 0.2$ MeV and $L({\rho_c}) = 46.0 \pm 4.5$ MeV at $\rho_{\rm{c}}= 0.11$ fm$^{-3}$ extracted from nuclear mass and the neutron skin thickness of Sn isotopes. Our results indicate that the symmetry energy cannot be stiffer than a linear density dependence.In addition, we also discuss the quark matter symmetry energy since the deconfined quarks could be the right degree of freedom in dense matter at high baryon densities.

We report a timing analysis of the black hole binary GRS 1915+105 with the NuSTAR observatory. A strong type-C QPO below 2 Hz appears in the power density spectrum during the whole observation, whose frequency is correlated with the 3-25 keV count rate. The QPO shows a sudden increase in frequency along with an increase in flux and a softening of the spectrum. We discuss the possible origin of the QPO and the reasons that lead to the QPO frequency variation. It is suggested that the reflection component has little influence on QPO frequency and the increase in QPO frequency could be associated with the inward motion of the outer part of the disk.

In this work, we provide 2189 photometric- and kinematic-selected member candidates of 24 star clusters from the LAMOST DR2 catalog. We perform two-step membership identification: selection along the stellar track in the color-magnitude diagram, i.e., photometric identification, and the selection from the distribution of radial velocities, i.e. the kinematic identification. We find that the radial velocity from the LAMOST data are very helpful in the membership identification. The mean probability of membership is 40\% for the radial velocity selected sample. With these 24 star clusters, we investigate the performance of the radial velocity and metallicity estimated in the LAMOST pipeline. We find that the systematic offset in radial velocity and metallicity are $0.85\pm1.26$\,\kms\ and $-0.08\pm0.04$\,dex, with dispersions of $5.47_{-0.71}^{+1.16}$\,\kms\ and $0.13_{-0.02}^{+0.04}$\,dex, respectively. Finally, we propose that the photometric member candidates of the clusters covered by the LAMOST footprints should be assigned higher priority so that more member stars can be observed.

We systematically analyze the 6-year {\it Fermi}/LAT data of the lobe-dominated quasars (LDQs) in the complete LDQ sample from 3CRR survey and report the discovery of high-energy $\gamma$-ray emission from 3C 275.1. The $\gamma$-ray emission of 3C 207 is confirmed and significant variability of the lightcurve is identified. We do not find statistically significant $\gamma$-ray emission from other LDQs. 3C 275.1 is the known $\gamma$-ray quasar with the lowest core dominance parameter (i.e., $R=0.11$). We also show that both the northern radio hotspot and parsec jet models can reasonably reproduce the $\gamma$-ray data. The parsec jet model, however, is favored by the potential $\gamma$-ray variability at the timescale of months. We suggest that some dimmer $\gamma$-ray LDQs will be detected in the future and LDQs could contribute non-negligibly to the extragalactic $\gamma$-ray background.

We perform a discrimination procedure with the spectral index diagram of TiO5 and CaH2+CaH3 to separate M giants from M dwarfs. Using the M giant spectra identified from the LAMOST DR1 with high signal-to-noise ratio (SNR), we have successfully assembled a set of M giant templates, which show more reliable spectral features. Combining with the M dwarf/subdwarf templates in Zhong et al. (2015), we present an extended M-type templates library which includes not only M dwarfs with well-defined temperature and metallicity grid but also M giants with subtype from M0 to M6. Then, the template-fit algorithm were used to automatically identify and classify M giant stars from the LAMOST DR1. The result of M giant stars catalog is cross-matched with 2MASS JHKs and WISE W1/W2 infrared photometry. In addition, we calculated the heliocentric radial velocity of all M giant stars by using the cross-correlation method with the template spectrum in a zero-velocity restframe. Using the relationship between the absolute infrared magnitude MJ and our classified spectroscopic subtype, we derived the spectroscopic distance of M giants with uncertainties of about 40%. A catalog of 8639 M giants is provided. As an additional search result, we also present 101690 M dwarfs/subdwarfs catalog which were classified by our classification pipeline.

We report on searching for Classical B-type emission-line (CBe) stars from the first data release (DR1) of the Large Sky Area Multi-Object fiber Spectroscopic Telescope (LAMOST; also named the Guoshoujing Telescope). A total of 192 (12 known CBes) objects were identified as CBe candidates with prominent He~I~$\lambda4387$, He~I~$\lambda4471$, and Mg~II~$\lambda4481$ absorption lines, as well as H$\beta$~$\lambda4861$ and H$\alpha$~$\lambda6563$ emission lines. These candidates significantly increases current CBe sample of about 8\%. Most of the CBe candidates are distributed at the Galactic Anti-Center due to the LAMOST observing strategy. Only two of CBes are in the star clusters with ages of 15.8 and 398~Myr, respectively.

We present a method to estimate distances to stars with spectroscopically derived stellar parameters. The technique is a Bayesian approach with likelihood estimated via comparison of measured parameters to a grid of stellar isochrones, and returns a posterior probability density function for each star's absolute magnitude. This technique is tailored specifically to data from the Large Sky Area Multi-object Fiber Spectroscopic Telescope (LAMOST) survey. Because LAMOST obtains roughly 3000 stellar spectra simultaneously within each ~5-degree diameter "plate" that is observed, we can use the stellar parameters of the observed stars to account for the stellar luminosity function and target selection effects. This removes biasing assumptions about the underlying populations, both due to predictions of the luminosity function from stellar evolution modeling, and from Galactic models of stellar populations along each line of sight. Using calibration data of stars with known distances and stellar parameters, we show that our method recovers distances for most stars within ~20%, but with some systematic overestimation of distances to halo giants. We apply our code to the LAMOST database, and show that the current precision of LAMOST stellar parameters permits measurements of distances with ~40% error bars. This precision should improve as the LAMOST data pipelines continue to be refined.

The Large sky Area Multi-Object Spectroscopic Telescope (LAMOST) General Survey is a spectroscopic survey that will eventually cover approximately half of the celestial sphere and collect 10 million spectra of stars, galaxies and QSOs. Objects both in the pilot survey and the first year general survey are included in the LAMOST First Data Release (DR1). The pilot survey started in October 2011 and ended in June 2012, and the data have been released to the public as the LAMOST Pilot Data Release in August 2012. The general survey started in September 2012, and completed its first year of operation in June 2013. The LAMOST DR1 includes a total of 1202 plates containing 2,955,336 spectra, of which 1,790,879 spectra have observed signal-to-noise S/N >10. All data with S/N>2 are formally released as LAMOST DR1 under the LAMOST data policy. This data release contains a total of 2,204,696 spectra, of which 1,944,329 are stellar spectra, 12,082 are galaxy spectra and 5,017 are quasars. The DR1 includes not only spectra, but also three stellar catalogues with measured parameters: AFGK-type stars with high quality spectra (1,061,918 entries), A-type stars (100,073 entries), and M stars (121,522 entries). This paper introduces the survey design, the observational and instrumental limitations, data reduction and analysis, and some caveats. Description of the FITS structure of spectral files and parameter catalogues is also provided.

The neutron-proton effective mass splitting in neutron-rich nucleonic matter reflects the space-time nonlocality of the isovector nuclear interaction. It affects the neutron/proton ratio during the earlier evolution of the Universe, cooling of protoneutron stars, structure of rare isotopes and dynamics of heavy-ion collisions. While there is still no consensus on whether the neutron-proton effective mass splitting is negative, zero or positive and how it depends on the density as well as the isospin-asymmetry of the medium, significant progress has been made in recent yeas in addressing these issues. We first recall the connections among the neutron-proton effective mass splitting, the momentum dependence of the isovector potential and the density dependence of the symmetry energy. We then make a few observations about the progress in calculating the neutron-proton effective mass splitting using various nuclear many-body theories and its effects on the isospin-dependence of in-medium nucleon-nucleon cross sections. Perhaps, our most reliable knowledge so far about the neutron-proton effective mass splitting at saturation density of nuclear matter comes from optical model analyses of huge sets of nucleon-nucleus scattering data accumulated over the last five decades. The momentum dependence of the symmetry potential from these analyses provide a useful boundary condition at saturation density for calibrating nuclear many-body calculations. Several observables in heavy-ion collisions have been identified as sensitive probes of the neutron-proton effective mass splitting in dense neutron-rich matter based on transport model simulations. We review these observables and comment on the latest experimental findings.

The color-magnitude diagram (CMD) of globular cluster NGC1651 has special structures including a broad main sequence, an extended main sequence turn-off and an extended red giant clump. The reason for such special CMDs remains unclear. In order to test how different the results from various stellar population assumptions are, we study a high-quality CMD of NGC1651 from the Hubble Space Telescope archive via eight kinds of models. Distance modulus, extinction, age ranges, star formation mode, fraction of binaries, and fraction of rotational stars are determined and then compared. The results show that stellar populations both with and without age spread can reproduce the special structure of the observed CMD. A composite population with extended star formation from 1.8\,Gyrs ago to 1.4\,Gyrs ago, which contains 50 per cent binaries and 70 per cent rotational stars, fits the observed CMD best. Meanwhile, a 1.5\,Gyr-old simple population that consists of rotational stars can also fit the observed CMD well. The results of CMD fitting are shown to depend strongly on stellar population type (simple or composite), and fraction of rotators. If the member stars of NGC1651 formed in a single star burst, the effect of stellar rotation should be very important for the explanation of observed CMDs. Otherwise, the effect may be small. It is also possible that the special observed CMD is a result of the combined effects of stellar binarity, rotation and age spread. Therefore, further work on stellar population type and fraction of rotational stars of intermediate-age clusters are necessary to understand their observed CMDs.

We compiled the radio, optical, and X-ray data of blazars from the Sloan Digital Sky Survey (SDSS) database, and presented the distribution of luminosities and broad band spectral indices. The distribution of luminosities shows that the averaged luminosity of flat-spectral radio quasars (FSRQs) is larger than that of BL Lacs objects. On the other hand, the broad band spectral energy distribution reveals that FSRQs and low energy peaked BL Lac objects (LBLs) objects have similar spectral properties, but high energy peaked BL Lac objects (HBLs) have a distinct spectral property. This may be due to that different subclasses of blazars have different intrinsic environments and are at different cooling levels. Even so, a unified scheme also is revealed from the color-color diagram, which hints that there are similar physical processes operating in all objects under a range of intrinsic physical conditions or beaming parameter.

Isospin-violating dark matter (IVDM), which couples differently with protons and neutrons, provides a promising mechanism to ameliorate the tension among recent direct detection experiments. Assuming DM is non-interacting bosonic asymmetric IVDM, we investigate how the existence of old neutron stars limits the DM-proton scattering cross-section $\sigma_{\rm p}$, especially the effects of the isospin violating DM-nucleon interactions and the symmetry energy in the equation of state of isospin asymmetric nuclear matter. Our calculations are completely based on general relativity and especially the structure of neutron stars is obtained by solving the Tolman-Oppenheimer-Volkoff equations with nuclear matter equation of state constrained by terrestrial experiments. We find that, by considering the more realistic neutron star model rather than a simple uniform neutron sphere as usual, the $\sigma_{\rm p}$ bounds from old neutron stars can be varied by more than an order of magnitude depending on the specific values of the DM neutron-to-proton coupling ratio $f_{\rm n}/f_{\rm p}$, and they can be further varied by more than a factor of two depending on the density dependence of the symmetry energy. In particular, we demonstrate that the observed nearby isolated old neutron star PSR B1257+12 can set a very strong limit on $\sigma_{\rm p}$ for low-mass DM particles ($ \le 20 \, {\rm GeV}$) that reaches a sensitivity beyond the current best limits from direct detection experiments and disfavors the DM interpretation of previously-reported positive experimental results, including the IVDM.

We investigate properties of strange quark matter in the framework of SU(3) Nambu--Jona-Lasinio(NJL) model with vector interaction under strong magnetic fields. The effects of vector-isoscalar and vector-isovector interaction on the equation of state of strange quark matter are investigated, and it is found that the equation of state is not sensitive to the vector-isovector interaction, however, a repulsive interaction in the vector-isoscalar channel gives a stiffer equation of state for cold dense quark matter. In the presence of magnetic field, gluons will be magnetized via quark loops, and the contribution from magnetized gluons to the equation of state is also estimated. The sound velocity square is a quantity to measure the hardness or softness of dense quark matter, and in the NJL model without vector interaction at zero magnetic field the sound velocity square is always less than 1/3. It is found that a repulsive vector-isoscalar interaction and a positive pressure contribution from magnetized gluons can enhance the sound velocity square, which can even reach 1. To construct quark magnetars under strong magnetic fields, we consider anisotropic pressures and use a density-dependent magnetic field profile to mimic the magnetic field distribution in a quark star. We also analyze the parameter region for the magnitude of vector-isoscalar interaction and the contribution from magnetized gluons in order to produce 2 solar mass quark magnetars.

Nuclear structure observables usually most effectively probe the properties of nuclear matter at subsaturation densities rather than at saturation density. We demonstrate that the electric dipole polarizibility $\alpha _ {\text{D}}$ in $^{208}$Pb is sensitive to both the magnitude $E_{\text{sym}}(\rho_{\text{c}})$ and density slope $L(\rho_{\text{c}})$ of the symmetry energy at a subsaturation cross density $\rho_{\text{c}} = 0.11$ fm$^{-3}$. Using the experimental data of $\alpha _ {\text{D}}$ in $^{208}$Pb from RCNP and the recent accurate constraint of $E_{\text{sym}}(\rho_{\text{c}})$ from the binding energy difference of heavy isotope pairs, we extract a value of $L(\rho_{\text{c}}) = 47.3 \pm 7.8$ MeV. The implication of the present constraint of $L(\rho_{\text{c}})$ to the symmetry energy at saturation density, the neutron skin thickness of $^{208}$Pb and the core-crust transition density in neutron stars is discussed.

We studied the horizontal branch oscillations (HBO) and the band-limited components observed in the power spectra of the transient neutron star low-mass X-ray binary XTE J1701-462 and the persistent "Sco-like" Z source GX 17+2. These two components were studied based on the state-resolved spectra. We found that the frequencies of XTE J1701-462 lie on the known correlations (WK and PBK), showing consistency with other types of X-ray binaries (black holes, atoll sources and millisecond X-ray pulsars). However, GX 17+2 is shifted from the WK correlation like other typical Z sources. We suggest that the WK/PBK main track forms a boundary which separates persistent sources from transient sources. The characteristic frequencies of break and HBO are independent of accretion rate in both sources, although it depends on spectral models. We also report the energy dependence of the HBO and break frequencies in XTE J1701-462 and how the temporal properties change with spectral state in XTE J1701-462 and GX 17+2. We studied the correlation between rms at the break and the HBO frequency. We suggest that HBO and break components for both sources probably arise from a similar physical mechanism: Comptonization emission from the corona. These two components could be caused by same kind of oscillation in a corona who with uneven density, and they could be generated from different areas of corona. We further suggest that different proportions of the Comptonization component in the total flux cause the different distribution between GX 17+2 and XTE J1701-462 in the $rms_{\rm{break}}$-$rms_{\rm{HBO}}$ diagram.

Photospheric radius expansion (PRE) bursts have already been used to constrain the masses and radii of neutron stars. RXTE observed three PRE bursts in 4U 1746-37, all with low touchdown fluxes. We discuss here the possibility of low mass neutron star in 4U 1746-37 because the Eddington luminosity depends on stellar mass. With typical values of hydrogen mass fraction and color correction factor, a Monte-Carlo simulation was applied to constrain the mass and radius of neutron star in 4U 1746-37. 4U 1746-37 has a high inclination angle. Two geometric effects, the reflection of the far side accretion disc and the obscuration of the near side accretion disc have also been included in the mass and radius constraints of 4U 1746-37. If the reflection of the far side accretion disc is accounted, a low mass compact object (mass of $0.41\pm0.14~M_{\odot}$ and radius of $8.73\pm1.54~\rm km$ at 68% confidence) exists in 4U 1746-37. If another effect operated, 4U 1746-37 may contain an ultra low mass and small radius object ($M=0.21\pm0.06~M_{\odot},~R=6.26\pm0.99~\rm km$ at 68% confidence). Combined all possibilities, the mass of 4U 1746-37 is $0.41^{+0.70}_{-0.30}~M_\odot$ at 99.7% confidence. For such low mass NS, it could be reproduced by a self-bound compact star, i.e., quark star or quark-cluster star.

We fit the (quasi-)simultaneous multi-waveband spectral energy distributions (SEDs) for a sample of low-synchrotron-peaked (LSP) blazars with a one-zone leptonic model. The seed photons that predominantly come from broad line region (BLR) and infrared (IR) molecular torus are considered respectively in external Compton process. We find that the modeling with IR seed photons is systematically better than that with BLR photons based on a $\chi^2$ test, which suggest that $\gamma$-ray emitting region most possibly stay outside the BLR. The minimum electron Lorentz factor, $\gamma_{\rm min}$, is constrained from the modeling for these LSP blazars with good soft X-ray data (ranges from 5 to 160 with a median value of 55), which plays a key role in jet power estimation. Assuming one-to-one ratio of proton and electron, we find that the jet power for LSP blazars is systematically higher than that of FR II radio galaxies at given 151 MHz radio luminosity, $L_{\rm 151MHz}$, even though FR IIs are regarded as same as LSP blazars in unification scheme except the jet viewing angle. The possible reason is that there are some $e^{\pm}$ pairs in the jet of these blazars. If this is the case, we find the number density of $e^{\pm}$ pairs should be several times higher than that of $e^{-}-p$ pairs by assuming the jet power is the same for LSP blazars and FR IIs at given $L_{\rm 151MHz}$.

We quantify and correct systematic errors in PPMXL proper motions using extragalactic sources from the first two LAMOST data releases and the Veron-Cetty & Veron Catalog of Quasars. Although the majority of the sources are from the Veron catalog, LAMOST makes important contributions in regions that are not well-sampled by previous catalogs, particularly at low Galactic latitudes and in the south Galactic cap. We show that quasars in PPMXL have measureable and significant proper motions, which reflect the systematic zero-point offsets present in the catalog. We confirm the global proper motion shifts seen by Wu, Ma, & Zhou (2011), and additionally find smaller-scale fluctuations of the QSO-derived corrections to an absolute frame. We average the proper motions of 158,106 extragalactic objects in bins of 3x3 degrees and present a table of proper motion corrections.

Within the confined isospin- and density-dependent mass model, we study the properties of strange quark matter (SQM) and quark stars (QSs) under strong magnetic fields. The equation of state of SQM under a constant magnetic field is obtained self-consistently and the pressure perpendicular to the magnetic field is shown to be larger than that parallel to the magnetic field, implying that the properties of magnetized QSs generally depend on both the strength and the orientation of the magnetic fields distributed inside the stars. Using a density-dependent magnetic field profile which is introduced to mimic the magnetic field strength distribution in a star, we study the properties of static spherical QSs by assuming two extreme cases for the magnetic field orientation in the stars, i.e., the radial orientation in which the local magnetic fields are along the radial direction and the transverse orientation in which the local magnetic fields are randomly oriented but perpendicular to the radial direction. Our results indicate that including the magnetic fields with radial (transverse) orientation can significantly decrease (increase) the maximum mass of QSs, demonstrating the importance of the magnetic field orientation inside the magnetized compact stars.

Isospin-violating dark matter (IVDM) provides a possible mechanism to ameliorate the tension among recent direct detection experiments. For IVDM, we demonstrate that the results of direct detection experiments based on neutron-rich target nuclei may depend strongly on the density dependence of the symmetry energy which is presently largely unknown and controls the neutron skin thickness that reflects the relative difference of neutron and proton form factors in the neutron-rich nuclei. In particular, using the neutron and proton form factors obtained from Skyrme-Hartree-Fock calculations by varying the symmetry energy within the uncertainty region set by the latest model-independent measurement of the neutron skin thickness of $^{208}$Pb from PREX experiment at JLab, we find that, for IVDM with neutron-to-proton coupling ratio fixed to $f_n/f_p=-0.7$, the form factor effect may enhance the sensitivity of Xe-based detectors (e.g., XENON100 and LUX) to the DM-proton cross section by a factor of $3$ in the DM mass region constrained by CMDS-II(Si) and even by more than an order of magnitude for heavy DM with mass larger than $80$ GeV, compared with the results using the empirical Helm form factor. Our results further indicate that the form factor effect can significantly modify the recoil spectrum of Xe-based detectors for heavy IVDM with $f_n/f_p=-0.7$.

BL Lacertae (BL Lac) objects at high redshifts ($z\geq 2$) are rarely detected. Through careful analysis of the SDSS spectrum, BZB J1450+5201 is confirmed to be a high-$z$ BL Lac object with $z\geq$ 2.471 by identifying the Ly$\alpha$ 1216 and CIV 1548/1550 absorption lines. This indicates that BZB J1450+5201 is the most distant BL Lac object discovered to date. Careful analysis of the five-year \fermi data of 2FGL J1451.0+5159 shows that its $\gamma$-ray emission is robust with confidence level of 6.2$\sigma$ at 1-3 GeV and 6.7$\sigma$ at 3-10 GeV, and that the confusion of bright neighbor is negligible, which can not be fixed in the analysis of the two-year data. Meanwhile, 2FGL J1451.0+5159 is confirmed to be associated with BZB J1450+5201 using the five-year data. The analysis of multiwavelength data, from radio to $\gamma$-ray energies, indicates BZB J1450+5201 is an intermediate synchrotron peaked (ISP) source and consistent with distributions of other ISP sources at lower redshifts in the second LAT AGN catalog. The pure SSC model seems to be disfavoured, while scattering of weak external emission plus SSC process can provide a satisfactory description of the broadband emission.

Langmuir waves (LWs), which are believed to play a crucial role in the plasma emission of solar radio bursts, can be excited by streaming instability of energetic electron beams. However, solar hard X-ray observations imply that the energetic flare electrons usually have a power-law energy distribution with a lower energy cutoff. In this paper, we investigate LWs driven by the power-law electrons. The results show that power-law electrons with the steepness cutoff behavior can excite LWs effectively because of the population inversion distribution below the cutoff energy ($E_c$). The growth rate of LWs increases with the steepness index ($\delta$) and decreases with the power-law index ($\alpha$). The wave number of the fastest growing LWs ($k\lambda_D$), decreases with the characteristic velocity of the power-law electrons ($v_{c}=\sqrt{2E_{c}/m_{e}}$) and increases with the thermal velocity of ambient electrons ($v_T$). This can be helpful for us to understand better the physics of LWs and the dynamics of energetic electron beams in space and astrophysical plasmas.

We report the discovery of 28 candidate high-velocity stars (HVSs) at heliocentric distances of less than 3 kpc, based on the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) Data Release 1. Our sample of HVS candidates covers a much broader color range than the equivalent ranges discussed in previous studies and comprises the first and largest sample of HVSs in the solar neighborhood. The sufficiently accurate observed and derived parameters for all candidates allow us to ascertain their nature as genuine HVSs, while a subset of 12 objects represents the most promising candidates. Our results also highlight the great potential of discovering statistically large numbers of HVSs of different spectral types in LAMOST survey data. This will ultimately enable us to achieve a better understanding of the nature of Galactic HVSs and their ejection mechanisms, and to constrain the structure of the Galaxy.

The SED of blazars show significant curvature. In this paper, we study the curvature properties for a large sample of Fermi/LAT bright blazars based on quasi-simultaneous SED. Both SEDs of synchrotron and inverse Compton (IC) components are fitted by a log-parabolic law in log_v-log_vfv diagram. The second-degree term of log-parabola measures the curvature of SED. We find a statistically significant correlation between synchrotron peak frequency and its curvature. This result is in agreement with the theoretical prediction, and confirms previous studies, which dealt with single source with various epoch observations or a small sample. If a broken power-law is employed to fit the SED, the difference between the two spectral indexes (alpha2-alpha1) can be considered as a "surrogate" of the SED curvature. We collect spectral parameters of a sample blazars from literature, and find a correlation between the synchrotron peak frequency and the spectral difference. We do not find a significant correlation between the IC peak frequency and its curvature, which may be caused by complicated seed photon field. It is also found that the synchrotron curvatures are on average larger than that of IC curvatures, and there is no correlation between these two parameters. As suggested by previous works in literature, both the log-parabolic law of SED and above correlation can be explained by statistical and/or stochastic particle accelerations. Stochastic particle acceleration predicts a different slope of the correlation from that of statistical one, and our result seems favor stochastic acceleration mechanisms and emission processes. Some of other evidences also seem to support that the electron energy distribution (and/or synchrotron SED) may be log-parabolic, which include SED modeling, particle acceleration simulation, and comparisons between some predictions and empirical relations/correlations.

The New Vacuum Solar Telescope (NVST) is a 1 meter vacuum solar telescope that aims to observe the fine structures on the Sun. The main tasks of NVST are high resolution imaging and spectral observations, including the measurements of solar magnetic field. NVST is the primary ground-based facility of Chinese solar community in this solar cycle. It is located by the Fuxian Lake of southwest China, where the seeing is good enough to perform high resolution observations. In this paper, we first introduce the general conditions of Fuxian Solar Observatory and the primary science cases of NVST. Then, the basic structures of this telescope and instruments are described in detail. Finally, some typical high resolution data of solar photosphere and chromosphere are also shown.

Solar coronal shocks are very common phenomena in the solar atmosphere and are believed to be the drivers of solar type II radio bursts. However, the microphysical nature of these emissions is still an open problem. This paper proposes that electron cyclotron maser (ECM) emission is responsible for the generation of radiations from the coronal shocks. In the present model, an energetic ion beam accelerated by the shock excites first Alfv\'en wave (AW) and then the excited AW leads to the formation of a density-depleted duct along the foreshock boundary of the shock. In this density-depleted duct, the energetic electron beam produced via the shock acceleration can effectively excite radio emission by the ECM instability. Our results show that this model may have potential application to solar type II radio bursts.

A thorough understanding of many astrophysical phenomena associated with compact objects requires reliable knowledge about both the equation of state (EOS)of super-dense nuclear matter and the theory of strong-field gravity simultaneously because of the EOS-gravity degeneracy. Currently, variations of the neutron star (NS) mass-radius correlation from using alternative gravity theories are much larger than those from changing the NS matter EOS within known constraints. At least two independent observables are required to break the EOS-gravity degeneracy. Using model EOSs for hybrid stars and a Yukawa-type non-Newtonian gravity, we investigate both the mass-radius correlation and pulsating frequencies of NSs. While the maximum mass of NSs increases, the frequencies of the $f$, $p_1$, $p_2$, and $w_I$ pulsating modes are found to decrease with the increasing strength of the Yukawa-type non-Newtonian gravity, providing a useful reference for future determination simultaneously of both the strong-field gravity and the supranuclear EOS by combining data of x-ray and gravitational wave emissions of neutron stars.

We study the temporal and energy spectral properties of the unique neutron star low-mass X-ray binary XTE J1701-462. After assuming the HB/NB vertex as a reference position of accretion rate, the horizontal branch oscillation (HBO) of the HB/NB vertex is roughly 50 Hz. It indicates that the HBO is independent with the accretion rate or the source intensity. The spectral analysis shows $R_{\rm{in}}\propto\dot{M}_{\rm{Disk}}^{2.9\pm0.09}$ in the HB/NB vertex and $R_{\rm{in}}\propto\dot{M}_{\rm{Disk}}^{1.7\pm0.06}$ in the NB/FB vertex, which implies that different accretion rate may be produced in the HB/NB vertex and the NB/FB vertex. The Comptonization component could be fitted by constrained broken power law (CBPL) or nthComp. Different with GX 17+2, the frequencies of HBO positively correlate with the inner disk radius, which contradict with the prediction of Lense-Thirring precession model. XTE J1701-462, both in the Cyg-like phase and in the Sco-like phase, follows a positive correlation between the break frequency of broad band noise and the characteristic frequency of HBO, which is called the W-K relation. An anticorrelation between the frequency of HBO and photon energy is observed. Moreover, the rms of HBO increases with photon energy till ~10 keV. We discuss the possible origin of HBO from corona in XTE J1701-462.

Using all the observations from Rossi X-ray Timing Explorer for accreting neutron star XTE J1701-462, we carry out a systematic study on the cross-correlation between its soft and hard light curves. The anti-correlations appear at the hard vertex and on the upper normal branch (NB) in the Cyg-like interval, occur on the horizontal branch (HB) and/or upper NB in the first and second Sco-like intervals, and display at the upper flaring branch (FB) in the third Sco-like interval. We suggest that the anti-correlation might evolve with luminosity, because with decreasing of luminosity, the source evolved from a Cyg-like Z source, via a Sco-like Z source, to an atoll source. In the Cyg-like interval, the positive correlations are mostly distributed on the HB, which is not consistent with that of the Cyg-like Z source GX 5-1 and Cyg X-2 whose HBs host ambiguous correlations and anti-correlations. From our spectral analyses, it is found that the ratio of the hard emission to the soft emission basically keeps unvaried from the HB to the hard vertex, which might result in the positive correlation. Additionally, we compare the spectra in the third Sco-like interval with those of atoll source 4U 1735-44 and find that their spectral evolution along the tracks in the color-color diagrams are similar, indicating that in this interval, XTE J1701-462 approaches an atoll source. Both hard time lags and soft time lags are detected and the possible origins of them are discussed.

The single-particle potentials for nucleons and hyperons in neutron-rich matter generally depends on the density and isospin asymmetry of the medium as well as the momentum and isospin of the particle. It further depends on the temperature of the matter if the latter is in thermal equilibrium. We review here the extension of a Gogny-type isospin- and momentum-dependent interaction in several aspects made in recent years and their applications in studying intermediate-energy heavy ion collisions, thermal properties of asymmetric nuclear matter and properties of neutron stars. The importance of the isospin- and momentum-dependence of the single-particle potential, especially the momentum dependence of the isovector potential, is clearly revealed throughout these studies.

S5 0716+714 is a typical BL Lacertae object. In this paper we present the analysis and results of long term simultaneous observations in the radio, near-infrared, optical, X-ray and $\gamma$-ray bands, together with our own photometric observations for this source. The light curves show that the variability amplitudes in $\gamma$-ray and optical bands are larger than those in the hard X-ray and radio bands and that the spectral energy distribution (SED) peaks move to shorter wavelengths when the source becomes brighter, which are similar to other blazars, i.e., more variable at wavelengths shorter than the SED peak frequencies. Analysis shows that the characteristic variability timescales in the 14.5 GHz, the optical, the X-ray, and the $\gamma$-ray bands are comparable to each other. The variations of the hard X-ray and 14.5 GHz emissions are correlated with zero-lag, so are the V band and $\gamma$-ray variations, which are consistent with the leptonic models. Coincidences of $\gamma$-ray and optical flares with a dramatic change of the optical polarization are detected. Hadronic models do not have the same nature explanation for these observations as the leptonic models. A strong optical flare correlating a $\gamma$-ray flare whose peak flux is lower than the average flux is detected. Leptonic model can explain this variability phenomenon through simultaneous SED modeling. Different leptonic models are distinguished by average SED modeling. The synchrotron plus synchrotron self-Compton (SSC) model is ruled out due to the extreme input parameters. Scattering of external seed photons, such as the hot dust or broad line region emission, and the SSC process are probably both needed to explain the $\gamma$-ray emission of S5 0716+714.

By cross-correlating an archival sample of 170 BL Lacs with 2 year \fermilat AGN sample, we have compiled a sample of 100 BL Lacs with \fermi detection (FBLs), and a sample of 70 non-\fermi BL Lacs (NFBLs). We compared various parameters of FBLs with those of NFBLs, including the redshift, the low frequency radio luminosity at 408 MHz ($L_{\rm 408MHz}$), the absolute magnitude of host galaxies ($M_{\rm host}$), the polarization fraction from NVSS survey ($P_{\rm NVSS}$), the observed arcsecond scale radio core flux at 5 GHz ($F_{\rm core}$) and jet Doppler factor; all the parameters are directly \textbf{measured} or derived from available data from literatures. We found that the Doppler factor is on average larger in FBLs than in NFBLs, and the $Fermi~ \gamma$-ray detection rate is higher in sources with higher Doppler factor. In contrast, there are no significant differences in terms of the intrinsic parameters of redshift, $ L_{\rm 408MHz}$, $ M_{\rm host}$ and $ P_{\rm NVSS}$. FBLs seem to have a higher probability of exhibiting measurable proper motion. These results strongly indicate a higher beaming effect in FBLs compared to NFBLs. The radio core flux is found to be strongly correlated with $\gamma$-ray flux, which remains after excluding the common dependence of the Doppler factor. At the fixed Doppler factor, FBLs have systematically larger radio core flux than NFBLs, implying lower $\gamma$-ray emission in NFBLs since the radio and $\gamma$-ray flux are significantly correlated. Our results indicate that the Doppler factor is an important parameter of $\gamma$-ray detection, the non-detection of $\gamma$-ray emission in NFBLs is likely due to low beaming effect, and/or low intrinsic $\gamma$-ray flux, and the gamma-rays are likely produced co-spatially with the arcsecond-scale radio core radiation and mainly through the SSC process.

We extend the confined-density-dependent-mass (CDDM) model to include isospin dependence of the equivalent quark mass. Within the confined-isospin-density-dependent-mass (CIDDM) model, we study the quark matter symmetry energy, the stability of strange quark matter, and the properties of quark stars. We find that including isospin dependence of the equivalent quark mass can significantly influence the quark matter symmetry energy as well as the properties of strange quark matter and quark stars. While the recently discovered large mass pulsars PSR J1614-2230 and PSR J0348+0432 with masses around $2M_{\odot}$ cannot be quark stars within the CDDM model, they can be well described by quark stars in the CIDDM model. In particular, our results indicate that the two-flavor $u$-$d$ quark matter symmetry energy should be at least about twice that of a free quark gas or normal quark matter within conventional Nambu-Jona-Lasinio model in order to describe the PSR J1614-2230 and PSR J0348+0432 as quark stars.

We collect the data of twin kilohertz quasi-periodic oscillations (kHz QPOs) published before 2012 from 26 neutron star (NS) low-mass X-ray binary (LMXB) sources, then we analyze the centroid frequency (\nu) distribution of twin kHz QPOs (lower frequency \nu_1 and upper frequency \nu_2) both for Atoll and Z sources. For the data without shift-and-add, we find that Atoll and Z sources show the different distributions of \nu_1, \nu_2 and \nu_2/\nu_1, but the same distribution of \Delta\nu (difference of twin kHz QPOs), which indicates that twin kHz QPOs may share the common properties of LXMBs and have the same physical origins. The distribution of \Delta\nu is quite different from constant value, so is \nu_2/\nu_1 from constant ratio. The weighted mean values and maxima of \nu_1 and \nu_2 in Atoll sources are slightly higher than those in Z sources. We also find that shift-and-add technique can reconstruct the distribution of \nu_1 and \Delta\nu. The K-S test results of \nu_1 and \Delta\nu between Atoll and Z sources from data with shift-and-add are quite different from those without it, and we think that this may be caused by the selection biases of the sample. We also study the properties of the quality factor (Q) and the root-mean-squared (rms) amplitude of 4U 0614+09 with the data from the two observational methods, but the errors are too big to make a robust conclusion. The NS spin frequency (\nu_s) distribution of 28 NS-LMXBs show a bigger mean value (about 408Hz) than that (about 281 Hz) of the radio binary millisecond pulsars (MSPs), which may be due to the lack of the spin detections from Z sources (systematically lower than 281 Hz). Furthermore, on the relations between the kHz QPOs and NS spin frequency \nu_s, we find the approximate correlations of the mean values of \Delta\nu with NS spin and its half, respectively.

We test the relativistic precession model (RPM) and the MHD Alfven wave oscillation model (AWOM) for the kHz QPOs by the sources with measured NS masses and twin kHz QPO frequencies. For RPM, the derived NS mass of Cyg X-2 (SAX J1808.4-3658 and 4U 1820-30) is 1.96 +/- 0.10 solar masses (2.83 +/- 0.04 solar masses and 1.85 +/- 0.02 solar masses), which is 30% (100% and 40%) higher than the measured result 1.5 +/- 0.3 solar masses (< 1.4 solar masses and 1.29 + 0.19 / - 0.07 solar masses). For AWOM, where the free parameter of model is the density of star, we infer the NS radii to be around 10 - 20 km for the above three sources, based on which we can infer the matter compositions inside NSs with the help of the equations of state (EOSs). In particular, for SAX J1808.4-3658, AWOM shows a lower mass density of its NS than those of the other known kHz QPO sources, with the radius range of 17 - 20 km, which excludes the strange quark matter inside its star.

We show that the neutron skin thickness $\Delta r_{np}$ of heavy nuclei is uniquely fixed by the symmetry energy density slope $L({\rho})$ at a subsaturation cross density $\rho_c \approx 0.11$ fm$^{-3}$ rather than at saturation density $\rho_0$, while the binding energy difference $\Delta E$ between a heavy isotope pair is essentially determined by the magnitude of the symmetry energy $E_{\text{sym}}({\rho})$ at the same $\rho_c$. Furthermore, we find a value of $L({\rho_c})$ leads to a negative $E_{\text{sym}}({\rho_{0}})$-$L({\rho_{0}})$ correlation while a value of $E_{\text{sym}}({\rho_{c})}$ leads to a positive one. Using data on $\Delta r_{np}$ of Sn isotopes and $\Delta E$ of a number of heavy isotope pairs, we obtain simultaneously $E_{\text{sym}}({\rho_{c})}=26.65\pm0.20$ MeV and $L({\rho_c})= 46.0\pm4.5$ MeV at 95% confidence level, whose extrapolation gives $E_{\text{sym}}({\rho_{0}})=32.3\pm1.0$ MeV and $L({\rho_{0}})=45.2\pm10.0$ MeV. The implication of these new constraints on the $\Delta r_{np}$ of $^{208}$Pb and the core-crust transition density in neutron stars is discussed.

According to the Hugenholtz-Van Hove theorem, the symmetry energy and its density slope can be decomposed uniquely in terms of the single-nucleon potential in asymmetric nuclear matter which, at normal density, can be constrained by the nucleon optical model potential extracted from analyzing the nucleon-nucleus scattering data. To more accurately extract information on the symmetry energy and its density slope at normal density from neutron-nucleus scattering data, going beyond the well-known Lane potential, we include consistently the second order terms in isospin asymmetry in both the optical model potential and the symmetry energy decomposition. We find that the strength of the second-order symmetry potential $U_{sym,2}$ in asymmetric nuclear matter is significant compared to the first-order one U_{sym,1}$, especially at high nucleon momentum. While the $U_{sym,1}$ at normal density decreases with nucleon momentum, the $U_{sym,2}$ is found to have the opposite behavior. Moreover, we discuss effects of the $U_{sym,1}$ and $U_{sym,2}$ on determining the density dependence of the symmetry energy, and we find that the available neutron-nucleus scattering data favor a softer density dependence of the symmetry energy at normal density with a negligible contribution from the $U_{sym,2}$.

We investigate the cross spectrum of XTE J1701-462 in various types of Neutron Star-Low Mass X-ray Binary subclasses, during its 2006-2007 outburst. We analyze the relation between the time lags and temporal variabilities. We find that the hard time lags accompany with the horizontal branch oscillations (HBOs) and the soft time lags dominate the noise in the low frequency range 0.1-10 Hz on HB. In the Cyg-like phase, the time lags decrease on the middle normal branch (NB) from HB/NB vertex to NB/FB vertex, whereas the time lags are roughly invariant in the Sco-like source. We discuss that the Compton upscattering by the corona introduces the soft lag in low frequency noise. We interpret that the variation of the Comptonization component from the disk emission lead to the HBOs' time lags evolution along the Z tracks. We also report the rms amplitude spectrum and phase lag spectrum for the NBO. The $\sim$ $160^{\circ}$ phase lag is found. We present that the rms amplitude of both the Cyg-like and the Sco-like NBO linearly increase with the photon energy in low energy bands, and it will drop in the highest energy band.

Effects of the non-Newtonian gravity on properties of finite nuclei are studied by consistently incorporating both the direct and exchange contribution of the Yukawa potential in the Hartree-Fock approach using a well-tested Skyrme force for the strong interaction. It is shown for the first time that the strength of the Yukawa term in the non-Newtonian gravity is limited to $\log(|\alpha|)<1.75/[\lambda(\rm fm)]^{0.54} + 33.6$ within the length scale of $\lambda=1-10$ fm in order for the calculated properties of finite nuclei not to be in conflict with accurate experimental data available.

About seventy percent of intermediate-age star clusters in the Large Magellanic Clouds have been confirmed to have broad main sequence, multiple or extended turn-offs and dual red giant clumps. The observed result seems against the classical idea that such clusters are simple stellar populations. Although many models have been used for explaining the results via factors such as prolonged star formation history, metallicity spread, differential redenning, selection effect, observational uncertainty, stellar rotation, and binary interaction, the reason for the special color-magnitude diagrams is still uncertain. We revisit this question via the combination of stellar rotation and binary effects. As a result, it shows "golf club" color-magnitude diagrams with broad or multiple turn-offs, dual red clump, blue stragglers, red stragglers, and extended main sequences. Because both binaries and massive rotators are common, our result suggests that most color-magnitude diagrams including extended turn-off or multiple turn-offs can be explained using simple stellar populations including both binary and stellar rotation effects, or composite populations with two components.

Binaries are very common in galaxies, and more than half of Galactic hot subdwarf stars, which are thought as a possible origin of UV-upturn of old stellar populations, are found in binaries. Previous works showed that binary evolution can make the spectra of binary star populations significantly different from those of single star populations. However, the effect of binary evolution has not been taken into account in most works of spectral fitting of galaxies. This paper studies the role of binary evolution in spectral fitting of early-type galaxies, via a stellar population synthesis model including both single and binary star populations. Spectra from ultraviolet to optical band are fitted to determine a few parameters of galaxies. The results show that the inclusion of binaries in stellar population models may lead to obvious change in the determination of some galaxy parameters and therefore it is potentially important for spectral studies. In particular, the ages of young components of composite stellar populations become much older when using binary star population models instead of single star population models. This implies that binary star population models will measure significantly different star formation histories (SFHs) for galaxies compared to single star population models. In addition, stellar population models with binary interactions measure larger dust extinctions than single star population models on average. It suggests that when using binary star population models instead of single star population models, negative extinctions are possibly unnecessary in spectral fitting of early-type galaxies.

Significant progress has been made in recent years in constraining nuclear symmetry energy at and below the saturation density of nuclear matter using data from both terrestrial nuclear experiments and astrophysical observations. However, many interesting questions remain to be studied especially at supra-saturation densities. In this lecture note, after a brief summary of the currently available constraints on nuclear symmetry energy near the saturation density we first discuss the relationship between the symmetry energy and the isopin and momentum dependence of the single-nucleon potential in isospin-asymmetric nuclear medium. We then discuss several open issues regarding effects of the tensor force induced neutron-proton short-range correlation (SRC) on nuclear symmetry energy. Finally, as an example of the impacts of nuclear symmetry energy on properties of neutron stars and gravitational waves, we illustrate effects of the high-density symmetry energy on the tidal polarizability of neutron stars in coalescing binaries.

We summarize the current status on constraining the density dependence of the symmetry energy from terrestrial laboratory measurements and astrophysical observations. While the value $E_{sym}({\rho_{0}})$ and density slope $L$ of the symmetry energy at saturation density $\rho_{0}$ can vary largely depending on the data or methods, all the existing constraints are essentially consistent with $E_{sym}({\rho_{0}}) = 31 \pm 2$ MeV and $L = 50 \pm 20$ MeV. The determination of the supra-saturation density behavior of the symmetry energy remains a big challenge.

The pulsating DA white dwarfs (ZZ Ceti stars) are $g$-mode non-radial pulsators. Asteroseismology provides strong constraints on their global parameters and internal structure. Since all the DA white dwarfs falling in the ZZ Ceti instability strip do pulsate, the internal structure derived from asteroseismology brings knowledge for the DA white dwarfs as a whole group. HS 0507+0434B is one of the ZZ Ceti stars which lies approximately in the middle of the instability strip for which we have undertaken a detailed asteroseismological study. We carried out multisite observation campaigns in 2007 and from December 2009 to January 2010. In total, 206 hours of photometric time-series have been collected. They have been analysed by means of Fourier analysis and simultaneous multi-frequency sine-wave fitting. In total, 39 frequency values are resolved including 6 triplets and a number of linear combinations. We identify the triplets as $\ell$=1 $g$-modes split by rotation. We derived the period spacing, the rotational splitting and the rotation rate. From the comparison of the observed periods with the theoretical periods of a series of models we estimate the fundamental parameters of the star: its total mass M$_{*}$/M$_{\odot}$ = 0.675, its luminosity L/L$_{\odot}$=3.5$\times 10^{-3}$, and its hydrogen mass fraction M$_{H}$/M$_{*}$= 10$^{-8.5}$.

We take the recently published data of twin kHz quasi-period oscillations (QPOs) in neutron star (NS) lowmass X-ray binaries (LMXBs) as the samples, and investigate the morphology of the samples, which focuses on the quality factor, peak frequency of kHz QPOs, and try to infer their physical mechanism. We notice that: (1) The quality factors of upper kHz QPOs are low (2 ~ 20 in general) and increase with the kHz QPO peak frequencies for both Z and Atoll sources. (2) The distribution of quality factor versus frequency for the lower kHz QPOs are quite different between Z and Atoll sources. For most Z source samples, the quality factors of lower kHz QPOs are low (usually lower than 15) and rise steadily with the peak frequencies except for Sco X-1, which drop abruptly at the frequency of about 750 Hz. While for most Atoll sources, the quality factors of lower kHz QPOs are very high (from 2 to 200) and usually have a rising part, a maximum and an abrupt drop. (3) There are three Atoll sources (4U 1728-34, 4U 1636-53 and 4U 1608-52) of displaying very high quality factors for lower kHz QPOs. These three sources have been detected with the spin frequencies and sidebands, in which the source with higher spin frequency presents higher quality factor of lower kHz QPOs and lower difference between sideband frequency and lower kHz QPO frequency.

A general gyrokinetic dispersion relation is gotten and is applied to analysis linear kinetic coupling of anisotropic firehose (or, kinetic Alfven wave) and mirror mode. Nyquist stability analysis is also given.

In this talk, we report our recent study of constraining the non-Newtonian gravity at femtometer scale. We incorporate the Yukawa-type non-Newtonian gravitational potential consistently to the Skyrme functional form using the exact treatment for the direct contribution and density-matrix expansion method for the exchange contribution. The effects from the non-Newtonian potential on finite nuclei properties are then studied together with a well-tested Skyrme force. Assuming that the framework without non-Newtonian gravity can explain the binding energies and charge radii of medium to heavy nuclei within 2% error, we set an upper limit for the strength of the non-Newtonian gravitational potential at femtometer scale.

Within a density-dependent relativistic mean-field model using in-medium meson-hadron coupling constants and meson masses, we explore effects of in-medium hyperon interactions on properties of neutron stars. It is found that the hyperonic constituents in large-mass neutron stars can not be simply ruled out, while the recently measured mass of the millisecond pulsar J1614-2230 can constrain significantly the in-medium hyperon interactions. Moreover, effects of nuclear symmetry energy on hyperonization in neutron stars are also discussed.

We describe the current plans for a spectroscopic survey of millions of stars in the Milky Way galaxy using the Guo Shou Jing Telescope (GSJT, formerly the Large Area Multi-Object Spectroscopic Telescope - LAMOST). The survey will obtain spectra for 2.5 million stars brighter than $r<19$ during dark/grey time, and 5 million stars brighter than $r<17$ or $J<16$ on nights that are moonlit or have low transparency. The survey will begin in fall of 2012, and will run for at least four years. The telescope design constrains the optimal declination range for observations to $10^\circ<\delta<50^\circ$, and site conditions lead to an emphasis on stars in the direction of the Galactic anticenter. The survey is divided into three parts with different target selection strategies: disk, anticenter, and spheroid. The resulting dataset will be used to study the merger history of the Milky Way, the substructure and evolution of the disks, the nature of the first generation of stars through identification of the lowest metallicity stars, and star formation through study of open clusters and the OB associations. Detailed design of the LEGUE survey will be completed after a review of the results of the pilot survey in summer 2012.

The weather at Xinglong Observing Station, where the Guo Shou Jing Telescope (GSJT) is located, is strongly affected by the monsoon climate in north-east China. The LAMOST survey strategy is constrained by these weather patterns. In this paper, we present a statistics on observing hours from 2004 to 2007, and the sky brightness, seeing, and sky transparency from 1995 to 2011 at the site. We investigate effects of the site conditions on the survey plan. Operable hours each month shows strong correlation with season: on average there are 8 operable hours per night available in December, but only 1-2 hours in July and August. The seeing and the sky transparency also vary with seasons. Although the seeing is worse in windy winters, and the atmospheric extinction is worse in the spring and summer, the site is adequate for the proposed scientific program of LAMOST survey. With a Monte Carlo simulation using historical data on the site condition, we find that the available observation hours constrain the survey footprint from 22h to 16h in right ascension; the sky brightness allows LAMOST to obtain the limit magnitude of V = 19.5mag with S/N = 10.

We describe the target selection algorithm for the low latitude disk portion of the LAMOST Pilot Survey, which aims to test systems in preparation for the LAMOST spectroscopic survey. We use the PPMXL (Roeser et al. 2010) astrometric catalog, which provides positions, proper motions, B/R/I magnitudes (mostly) from USNO-B (Monet et al. 2003) and J/H/Ks from The Two Micron All Sky Survey (2MASS, see Skrutskie et al. 2006) as well. We chose 8 plates along the Galactic plane, in the region $0^\circ<\alpha<67^\circ$ and $42^\circ<\delta<59^\circ$, that cover 22 known open clusters with a range of ages. Adjacent plates may have small overlapping. Each plate covers an area $2.5^\circ$ in radius,with central star (for Shack-Hartmann guider) brighter than $\sim8^{\rm th}$ magnitude. For each plate, we create an input catalog in the magnitude range $11.3<Imag<16.3$ and $Bmag$ available from PPMXL. The stars are selected to satisfy the requirements of the fiber positioning system and have a uniform distribution in the $I$ vs. $B-I$ color-magnitude diagram. Our final input catalog consists of 12,000 objects on each of 8 plates that are observable during the winter observing season in Xinglong Station of the National Astronomical Observatory of China.

We use samples of local main-sequence stars to show that the radial gradient of [Fe/H] in the thin disk of the Milky Way decreases with mean effective stellar temperature. Many of these stars are visiting the solar neighborhood from the inner and outer Galaxy. We use the angular momentum of each star about the Galactic center to determine the guiding center radius and to eliminate the effects of epicyclic motion, which would otherwise blur the estimated gradients. We interpret the effective temperature as a proxy for mean age, and conclude that the decreasing gradient is consistent with the predictions of radial mixing due to transient spiral patterns. We find some evidence that the trend of decreasing gradient with increasing mean age breaks to a constant gradient for samples of stars whose main-sequence life-times exceed the likely age of the thin disk.

We study the sub-AU-scale circumstellar environment of the Herbig Ae star HD144432 with near-infrared (NIR) VLTI/AMBER observations to investigate the structure of its inner dust disk. The interferometric observations were carried out with the AMBER instrument in the H and K band. We interpret the measured H- and K-band visibilities, the near- and mid-infrared visibilities from the literature, and the SED of HD144432 by using geometric ring models and ring-shaped temperature-gradient disk models with power-law temperature distributions. We derived a K-band ring-fit radius of 0.17 \pm 0.01 AU and an H-band radius of 0.18 \pm 0.01 AU (for a distance of 145 pc). This measured K-band radius of \sim0.17 AU lies in the range between the dust sublimation radius of \sim0.13 AU (predicted for a dust sublimation temperature of 1500 K and gray dust) and the prediction of models including backwarming (\sim0.27 AU). We found that an additional extended halo component is required in both the geometric and temperature-gradient modeling. In the best temperature- gradient model, the disk consists of two components. The inner part of the disk is a thin ring with an inner radius of \sim0.21 AU, a temperature of \sim1600 K, and a ring thickness \sim0.02 AU. The outer part extends from \sim1 AU to \sim10 AU with an inner temperature of \sim400 K. We find that the disk is nearly face-on with an inclination angle of < 28 degree. Our temperature-gradient modeling suggests that the NIR excess is dominated by emission from a narrow, bright rim located at the dust sublimation radius, while an extended halo component contributes \sim6% to the total flux at 2 {\mu}m. The MIR model emission has a two-component structure with \sim20% flux from the inner ring and the rest from the outer part. This two-component structure suggests a disk gap, which is possibly caused by the shadow of a puffed-up inner rim.

We explore effects of the light vector $U$-boson, which is weakly coupled to nucleons, on the transition density $\rho_{t}$ and pressure $P_{t}$ at the inner edge separating the liquid core from the solid crust of neutron stars. Three methods, i.e., the thermodynamical approach, the curvature matrix approach and the Vlasov equation approach are used to determine the transition density $\rho_{t}$ with the Skyrme effective nucleon-nucleon interactions. We find that the $\rho_{t}$ and $P_{t}$ depend on not only the ratio of coupling strength to mass squared of the $U$-boson $g^{2}/\mu ^{2}$ but also its mass $\mu $ due to the finite range interaction from the $U$-boson exchange. In particular, our results indicate that the $\rho_{t}$ and $P_{t}$ are sensitive to both $g^{2}/\mu ^{2}$ and $\mu $ if the $U$-boson mass $\mu $ is larger than about 2 MeV. Furthermore, we show that both $g^{2}/\mu ^{2}$ and $\mu $ can have significant influence on the mass-radius relation and the crustal fraction of total moment of inertia of neutron stars. In addition, we study the exchange term contribution of the $U$-boson based on the density matrix expansion method, and demonstrate that the exchange term effects on the nuclear matter equation of state as well as the $\rho_{t}$ and $P_{t}$ are generally negligible.

The nuclear symmetry energy $E_{sym}(\rho)$ and its density slope $L(\rho)$ can be decomposed analytically in terms of the single-nucleon potential in isospin asymmetric nuclear matter. Using three popular nuclear effective interaction models which have been extensively used in nuclear structure and reaction studies, namely, the isospin and momentum dependent MDI interaction model, the Skyrme Hartree-Fock approach and the Gogny Hartree-Fock approach, we analyze the contribution of different terms in the single-nucleon potential to the $E_{sym}(\rho)$ and $L(\rho)$. Our results show that the observed different density behaviors of $E_{sym}(\rho)$ for different interactions are essentially due to the variation of the symmetry potential $U_{sym,1}(\rho,k)$. Furthermore, we find that the contribution of the second-order symmetry potential $U_{sym,2}(\rho,k)$ to the $L(\rho)$ generally cannot be neglected. Moreover, our results demonstrate that the magnitude of the $U_{sym,2}(\rho,k)$ is generally comparable with that of $U_{sym,1}(\rho,k)$, indicating the second-order symmetry potential $U_{sym,2}(\rho,k)$ may have significant corrections to the widely used Lane approximation to the single-nucleon potential in extremely neutron(proton)-rich nuclear matter.

Within the nonlinear relativistic mean field model, we derive the analytical expression of the nuclear matter fourth-order symmetry energy $E_{sym,4}(\rho)$. Based on two accurately calibrated interactions FSUGold and IU-FSU, our results show that the value of $E_{sym,4}(\rho)$ at normal nuclear matter density $\rho_{0}$ is generally less than 1 MeV, confirming the empirical parabolic approximation to the equation of state for asymmetric nuclear matter at $\rho_{0}$. On the other hand, we find that the $E_{sym,4}(\rho)$ may become nonnegligible at high densities. Furthermore, the analytical form of the $E_{sym,4}(\rho)$ provides the possibility to study the higher-order effects on the isobaric incompressibility of asymmetric nuclear matter, i.e., $K_{sat}(\delta)=K_{0}+K_{sat,2}\delta ^{2}+K_{sat,4}\delta ^{4}+\mathcal{O}(\delta ^{6})$ where $\delta =(\rho_{n}-\rho_{p})/\rho $ is the isospin asymmetry, and we find that the value of $K_{sat,4}$ is generally small compared with that of the $K_{sat,2}$. In addition, we study the effects of the $E_{sym,4}(\rho)$ on the proton fraction $x_{p}$ and the core-crust transition density $\rho_{t}$ and pressure $P_{t}$ in neutron stars. Interestingly, we find that, compared with the results from the empirical parabolic approximation, including the $E_{sym,4}(\rho)$ contribution can significantly enhance the $x_{p}$ at high densities and strongly reduce the $\rho_{t}$ and $P_{t}$ in neutron stars, demonstrating that the widely used empirical parabolic approximation may cause large errors in determining the $x_{p}$ at high densities as well as the $\rho_{t}$ and $P_{t}$ in neutron stars within the nonlinear relativistic mean field model, consistent with previous nonrelativistic calculations.

Using the Hugenholtz-Van Hove theorem, we derive analytical expressions for the nuclear symmetry energy $E_{sym}(\rho)$ and its density slope $L(\rho)$ in terms of the Lorentz covariant nucleon self-energies in isospin asymmetric nuclear matter. These general expressions are useful for determining the density dependence of the symmetry energy and understanding the Lorentz structure and the microscopic origin of the nuclear symmetry energy in relativistic covariant formulism. As an example, we analyze the Lorentz covariant nucleon self-energy decomposition of $E_{sym}(\rho)$ and $L(\rho)$ and derive the corresponding analytical expressions within the nonlinear $\sigma$-$\omega$-$\rho$-$\delta$ relativistic mean field model.

Using the recently published data of twin kHz quasi-period oscillations (QPOs) in neutron star low-mass X-ray binaries (LMXBs), we study the different profiles between bright Z sources and less luminous Atoll sources in a statistical way. We find the quality factors of upper kHz QPOs show a narrow distribution both for Z sources and Atoll sources, which concentrate at 7.98 and 9.75 respectively, the quality factors of lower kHz QPOs show a narrow distribution for Z sources and a broader distribution for Atoll sources, which concentrate at 5.25 and 86.22 respectively. In order to investigate the relation between the quality factor and the peakfrequency of kHz QPOs, we fit the data with power-law, linear and exponential functions, respectively. There is an obvious trend that the quality factors increase with the peak frequencies both for upper and lower QPOs. The implications of our results are discussed.

We present a new technique to calculate the spectral lags of gamma-ray bursts (GRBs). Unlike previous processing methods, we first smooth the light curves of gamma-ray bursts in high and low energy bands using the "Loess" filter, then, we directly define the spectral lags as such to maximize the cross-correlation function (CCF) between two smoothed light curves. This method is suitable for various shapes of CCF; it effectively avoids the errors caused by manual selections for the fitting function and fitting interval. Using the method, we have carefully measured the spectral lags of individual pulses contained in BAT/Swift gamma-ray bursts with known redshifts, and confirmed the anti-correlation between the spectral lag and the isotropy luminosity. The distribution of spectral lags can be well fitted by four Gaussian components, with the centroids at 0.03 s, 0.09 s, 0.15 s, and 0.21 s, respectively. We find that some spectral lags of the multi-peak GRBs seem to evolve with time.

We use the accretion disk/corona+jet model to fit the multi-band spectral energy distributions (SEDs) of two unusual radio-intermediate/quiet quasars. It is found that the optical/UV emission of III Zw 2 is probably dominated by the emission from the accretion disk. The X-ray emission should be dominated by the radiation from the jet, while the contribution of the disk corona is negligible. The optical/UV component in the SED of PG 1407+265 can be well modeled as the emission from the accretion disk, while the IR component is attributed to the thermal radiation from the dust torus with an opening angle ~ 50\circ. If the X-ray continuum emission is dominated by the synchrotron emission of the jet, the source should be a "high peak frequency blazar", which obviously deviates the normal blazar sequence. The observed SED can also be fitted quite well by the accretion disk/corona model with the viscosity parameter ? = 0:5. The spectrum of the accretion disk/corona in PG 1407+265 satisfies the weak line quasar criterion suggested in Laor & Davis.

Based on microscopic Hartree-Fock + random phase approximation calculations with Skyrme interactions, we study the correlations between the nuclear breathing mode energy $E_{ISGMR}$ and properties of asymmetric nuclear matter with a recently developed analysis method. Our results indicate that the $E_{ISGMR}$ of $^{208}$Pb exhibits moderate correlations with the density slope $L$ of the symmetry energy and the isoscalar nucleon effective mass $m_{s,0}^{\ast}$ besides a strong dependence on the incompressibility $K_{0}$ of symmetric nuclear matter. Using the present empirical values of $L=60\pm 30$ MeV and $m_{s,0}^{\ast}=(0.8\pm 0.1)m$, we obtain a theoretical uncertainty of about $\pm 16$ MeV for the extraction of $K_{0}$ from the $E_{ISGMR}$ of $^{208}$Pb. Furthermore, we find the $E_{ISGMR}$ difference between $^{100}$Sn and $^{132}$Sn strongly correlates with $L$ and thus provides a potentially useful probe of the symmetry energy.

A wide field of view Cherenkov/fluorescence telescope array is one of the main components of the Large High Altitude Air Shower Observatory project. To serve as Cherenkov and fluorescence detectors, a flexible and mobile design is adopted for easy reconfiguring of the telescope array. Two prototype telescopes have been constructed and successfully run at the site of the ARGO-YBJ experiment in Tibet. The features and performance of the telescopes are presented.

In this paper, we construct a new phenomenological isospin dependent global neutron-nucleus optical model potential. Based on the existing experimental data of elastic scattering angular distributions for neutron as projectile, we obtain a set of the isospin dependent global neutron-nucleus optical model potential parameters, which can basically reproduce the experimental data for target nuclei from $^{24}$Mg to $^{242}$Pu with the energy region up to 200 MeV.

We compile 23 Gamma-ray Bursts (GRBs) and 21 blazars with estimated Doppler factors, and the Doppler factors of GRBs are estimated from their Lorentz factors by assuming their jet viewing angles \theta->0 degree. Using the conventional assumption that the prompt emission of GRBs is dominated by the synchrotron radiation, we calculate the synchrotron luminosity of GRBs from their total isotropic energy and burst duration. Intriguingly, we discover a uniform correlation between the synchrotron luminosity and Doppler factor, L_syn \propto D^3.1, for GRBs and blazars, which suggests that they may share some similar jet physics. One possible reason is that GRBs and blazars have, more or less, similar intrinsic synchrotron luminosities and both of them are strongly enhanced by the beaming effect. After Doppler and redshift-correction, we find that the intrinsic peak energy of the GRBs ranges from 0.1 to 3 keV with a typical value of 1 keV. We further correct the beaming effect for the observed luminosity of GRBs and find that there exists a positive correlation between the intrinsic synchrotron luminosity and peak energy for GRBs, which is similar to that of blazars. Our results suggest that both the intrinsic positive correlation and the beaming effect may be responsible for the observed tight correlation between the isotropic energy and the peak energy in GRBs (so called "Amati" relation).

In this paper, we use the quasi-simultaneous spectra of Fermi bright blazars and Fermi detected narrow line Seyfert 1 (NLS1) to study the blazar sequence and inverse Compton (IC) models. I. The synchrotron peak luminosities (L_{s}) significantly inverse correlate with the synchrotron peak frequencies (\nu_{s}), L_{s}\propto\nu_{s}^{-0.44}, which is consistent with the blazar sequence. In addition to the correlation, there are some blazars showing low \nu_{s} and low L_{s}. To study the relation between these low \nu_{s} low L_{s} blazars and the blazar sequence, we present correlations of the parameter L_{s}\nu_{s}^{1/4} with the ratio of Compton to synchrotron peak frequencies (r_{Cs}\equiv\nu_{C}/\nu_{s}) and with the ratio of Compton to synchrotron luminosities (CD\equiv L_{C}/L_{s}). The results indicate that both correlations are significant with a Pearson's probability for null correlation of p=0.0218 and p=0.0286 respectively. This does not support the idea that the low \nu_{s} low L_{s} blazars are sources with less beaming. Another possibility, as suggested by Ghisellini & Tavecchio, is that these blazars have relative lower black hole masses. To test this, we collect the black hole masses of 30 blazars from archives, and find that the hole mass correlates with the parameter L_{s}\nu_{s}^{0.44} (p=0.0344). Therefore, the black hole masses of low \nu_{s} low L_{s} blazars are statistically small. The NLS1s are thought to have lower black hole masses. We find that the four NLS1s detected by Fermi have low \nu_{s} and low L_{s}. This supports the above result. II. The ratio r_{Cs} correlates with CD significantly (p=0.00375). The external Compton (EC) model can naturally explain this correlation, while synchrotron self Compton (SSC) model can not. This agrees with the findings of many authors that the EC process dominates the gamma-ray emission of Flat Spectrum Radio Quasars.

The bulk parameters characterizing the energy of symmetric nuclear matter and the symmetry energy defined at normal nuclear density $\rho_0 $ provide important information on the equation of state (EOS) of isospin asymmetric nuclear matter. While significant progress has been made in determining some lower order bulk characteristic parameters, such as the energy $E_0(\rho_0)$ and incompressibility $K_0$ of symmetric nuclear matter as well as the symmetry energy $E_{sym}(\rho_0)$ and its slope parameter $L$, yet the higher order bulk characteristic parameters are still poorly known. Here, we analyze the correlations between the lower and higher order bulk characteristic parameters within the framework of Skyrme Hartree-Fock energy density functional and then estimate the values of some higher order bulk characteristic parameters. In particular, we obtain $J_0=-355 \pm 95$ MeV and $I_0=1473 \pm 680$ MeV for the third-order and fourth-order derivative parameters of symmetric nuclear matter at $\rho_0 $ and $K_{sym} = -100 \pm 165$ MeV, $J_{sym} = 224 \pm 385$ MeV, $I_{sym} = -1309 \pm 2025$ MeV for the curvature parameter, third-order and fourth-order derivative parameters of the symmetry energy at $\rho_0 $, using the empirical constraints on $E_0(\rho_0)$, $K_0$, $E_{sym}(\rho_0)$, $L$, and the isoscalar and isovector nucleon effective masses. Furthermore, our results indicate that the three parameters $E_0(\rho_0)$, $K_0$, and $J_0$ can reasonably characterize the EOS of symmetric nuclear matter up to $2\rho_0 $ while the symmetry energy up to $2\rho_0 $ can be well described by $E_{sym}(\rho_0)$, $L$, and $K_{sym}$.

We show in the Skyrme-Hartree-Fock approach that unambiguous correlations exist between observables of finite nuclei and nuclear matter properties. Using this correlation analysis to existing data on the neutron skin thickness of Sn isotopes, we find important constraints on the value E_{sym}(rho_0) and density slope L of the nuclear symmetry energy at saturation density. Combining these constraints with those from recent analyses of isospin diffusion and double neutron/proton ratio in heavy ion collisions leads to a value of L=58\pm 18 MeV approximately independent of E_{sym}(\rho_0).

Significant progress has been made in recent years in constraining the density dependence of nuclear symmetry energy using terrestrial nuclear laboratory data. Around and below the nuclear matter saturation density, the experimental constraints start to merge in a relatively narrow region. At supra-saturation densities, there are, however, still large uncertainties. After summarizing the latest experimental constraints on the density dependence of nuclear symmetry energy, we highlight a few recent studies examining imprints of nuclear symmetry energy on the binding energy, energy release during hadron-quark phase transitions as well as the $w$-mode frequency and damping time of gravitational wave emission of neutron stars.

Within the Skyrme-Hartree-Fock (SHF) approach, we show that for a fixed mass number A, both the symmetry energy coefficient a_{sym}(A) in the semi-empirical mass formula and the nuclear matter symmetry energy E_{sym}(\rho_A) at a subsaturation reference density rho_A can be determined essentially by the symmetry energy E_{sym}(rho_0) and its density slope L at saturation density rho_0. Meanwhile, we find the dependence of a_{sym}(A) on E_{sym}(rho_0) or L is approximately linear and is very similar to the corresponding linear dependence displayed by E_{sym}(\rho_A), providing an explanation for the relation E_{sym}(\rho_A) \approx a_{sym}(A). Our results indicate that a value of E_{sym}(\rho_A) leads to a linear correlation between E_{sym}(rho_0) and L and thus can put important constraints on E_{sym}(rho_0) and L. Particularly, the values of E_{sym}(rho_0)= 30.5 +- 3 MeV and L= 52.5 +- 20 MeV are simultaneously obtained by combining the constraints from recently extracted E_{sym}(\rho_A=0.1 fm^{-3}) with those from recent analyses of neutron skin thickness of Sn isotopes in the same SHF approach.

We investigate the evolution of young star clusters using N-body simulations. We confirm that subvirial and fractal-structured clusters will dynamically mass segregate on a short timescale (within 0.5 Myr). We adopt a modified minimum-spanning-tree (MST) method to measure the degree of mass segregation, demonstrating that the stars escaping from a cluster's potential are important for the temporal dependence of mass segregation in the cluster. The form of the initial velocity distribution will also affect the degree of mass segregation. If it depends on radius, the outer parts of the cluster would expand without undergoing collapse. In velocity space, we find 'inverse mass segregation,' which indicates that massive stars have higher velocity dispersions than their lower-mass counterparts.

Probably No! As an example, using soft EOSs consistent with existing terrestrial nuclear laboratory experiments for hybrid neutron stars containing a quark core described with MIT bag model using reasonable parameters, we show that the recently discovered new holder of neutron star maximum mass PSR J1614-2230 of $1.97\pm0.04M_{\odot}$ can be well described by incorporating a Yukawa gravitational correction that is consistent with existing constraints from neutron-proton and neutron-lead scatterings as well as the spectroscopy of antiproton atoms.

We present an {\sl XMM-Newton} observation of the eclipsing binary Algol which contains an X-ray dark B8V primary and an X-ray bright K2IV secondary. The observation covered the optical secondary eclipse and captured an X-ray flare that was eclipsed by the B star. The EPIC and RGS spectra of Algol in its quiescent state are described by a two-temperature plasma model. The cool component has a temperature around 6.4$\times 10^{6}$ K while that of the hot component ranges from 2 to 4.0$\times 10^{7}$ K. Coronal abundances of C, N, O, Ne, Mg, Si and Fe were obtained for each component for both the quiescent and the flare phases, with generally upper limits for S and Ar, and C, N, and O for the hot component. F-tests show that the abundances need not to be different between the cool and the hot component and between the quiescent and the flare phase with the exception of Fe. Whereas the Fe abundance of the cool component remains constant at $\sim$0.14, the hot component shows an Fe abundance of $\sim$0.28, which increases to $\sim$0.44 during the flare. This increase is expected from the chromospheric evaporation model. The absorbing column density $N_H$ of the quiescent emission is 2.5$\times10^{20}$ cm$^{-2}$, while that of the flare-only emission is significantly lower and consistent with the column density of the interstellar medium. This observation substantiates earlier suggestions of the presence of X-ray absorbing material in the Algol system.

AGNs with hard $\gamma$-ray emission identified so far are radio-loud. III Zw 2 is a radio intermediate AGN with a relativistic jet. We study its spectral energy distribution (SED) and find that the broad band emissions are dominated by the non-thermal emissions from the jet. We model its SED through a synchrotron + inverse Compton (IC) model. The results show that the IC component of III Zw 2 peaks at a few MeV, and the flux density drops rapidly at higher energy with photon index $\Gamma\approx3.3$ above 0.1 GeV. The predicted flux is slightly over the sensitivity of $Fermi$/LAT, but it is not included in the first $Fermi$/LAT AGN catalog. The reason for this may be: 1) that the IC peak is low and the spectrum is very steep above 0.1 GeV, 2) that III Zw 2 is in the low state during the period of the $Fermi$/LAT operation. We also find that III Zw 2 follows similar jet processes as those in $\gamma$-ray AGNs, e.g., the relation between jet power and radiation power, called the blazar sequence. We suggest that III Zw 2 may be a young source at an earlier stage of jet activity.

Expressing explicitly the parameters of the standard Skyrme interaction in terms of the macroscopic properties of asymmetric nuclear matter, we show in the Skyrme-Hartree-Fock approach that unambiguous correlations exist between observables of finite nuclei and nuclear matter properties. We find that existing data on neutron skin thickness $\Delta r_{np}$ of Sn isotopes give an important constraint on the symmetry energy $E_{sym}({\rho _{0}})$ and its density slope $L$ at saturation density ${\rho _{0}}$. Combining these constraints with those from recent analyses of isospin diffusion and double neutron/proton ratio in heavy-ion collisions at intermediate energies leads to a more stringent limit on $L$ approximately independent of $E_{sym}({\rho _{0}})$. The implication of these new constraints on the $\Delta r_{np}$ of $^{208}$Pb as well as the core-crust transition density and pressure in neutron stars is discussed.

The Large High Altitude Air Shower Observatory project is proposed to study high energy gamma ray astronomy ( 40 GeV-1 PeV ) and cosmic ray physics ( 20 TeV-1 EeV ). The wide field of view Cherenkov telescope array, as a component of the LHAASO project, will be used to study energy spectrum and compositions of cosmic ray by measuring the total Cherenkov light generated by air showers and shower maximum depth. Two prototype telescopes have been in operation since 2008. The pointing accuracy of each telescope is crucial to the direction reconstruction of the primary particles. On the other hand the primary energy reconstruction relies on the shape of the Cherenkov image on the camera and the unrecorded photons due to the imperfect connections between photomultiplier tubes. UV bright stars are used as point-like objects to calibrate the pointing and to study the optical properties of the camera, the spot size and the fractions of unrecorded photons in the insensitive areas of the camera.

Based on the Hugenholtz-Van Hove theorem, it is shown that both the symmetry energy E$_{sym}(\rho)$ and its density slope $L(\rho)$ at normal density $\rho_0$ are completely determined by the global nucleon optical potentials that can be extracted directly from nucleon-nucleus scatterings, (p,n) charge exchange reactions and single-particle energy levels of bound states. Adopting a value of $m^*/m=0.7$ for the nucleon effective k-mass in symmetric nuclear matter at $\rho_0$ and averaging all phenomenological isovector nucleon potentials constrained by world data available in the literature since 1969, the best estimates of $E_{sym}(\rho_0)=31.3$ MeV and $L(\rho_0)=52.7$ MeV are simultaneously obtained. Uncertainties involved in the estimates are discussed.

PG 1407+265 is a radio quiet quasar but has a relativistic jet. In this report, we show some peculiar properties of its optical and X-ray emissions, which indicate their possible non-thermal origins produced from jet. We use a simple synchrotron + synchrotron self Compton (SSC) model to fit the emissions with different ratios of energy densities between magnetic field and electrons (\eta=UB/Ue), which predict different {\gamma}-ray luminosity. The First LAT AGN Catalog (1LAC) did not include PG 1407+265, which indicates an upper limit of {\gamma}-ray luminosity. This upper limit constrains the ratio unreasonable large (\eta>10^4-10^5). This inversely indicates that the optical and X-ray emissions may be not produced from the beaming jet. We discuss the physical implications of these results.

Improved wavelength calibrators for high-resolution astrophysical spectrographs will be essential for precision radial velocity (RV) detection of Earth-like exoplanets and direct observation of cosmological deceleration. The astro-comb is a combination of an octave-spanning femtosecond laser frequency comb and a Fabry-P\'erot cavity used to achieve calibrator line spacings that can be resolved by an astrophysical spectrograph. Systematic spectral shifts associated with the cavity can be 0.1-1 MHz, corresponding to RV errors of 10-100 cm/s, due to the dispersive properties of the cavity mirrors over broad spectral widths. Although these systematic shifts are very stable, their correction is crucial to high accuracy astrophysical spectroscopy. Here, we demonstrate an \emph{in-situ} technique to determine the systematic shifts of astro-comb lines due to finite Fabry-P\'erot cavity dispersion. The technique is practical for implementation at a telescope-based spectrograph to enable wavelength calibration accuracy better than 10 cm/s.

The isospin- and momentum-dependent MDI interaction, which has been extensively used in intermediate-energy heavy-ion reactions to study the properties of asymmetric nuclear matter, is extended to include the nucleon-hyperon and hyperon-hyperon interactions by assuming same density, momentum and isospin dependence as for the nucleon-nucleon interaction. The parameters in these interactions are determined from the empirical hyperon single-particle potentials in symmetric nuclear matter at saturation density. The extended MDI interaction is then used to study in the mean-field approximation the equation of state of hypernuclear matter and also the properties of hybrid stars by including the phase transition from the hypernuclear matter to the quark matter at high densities. In particular, the effects of attractive and repulsive $\Sigma$N interactions and different values of symmetry energies on the hybrid star properties are investigated.

Using the momentum-dependent MDI effective interaction for nucleons, we have studied the transition density and pressure at the boundary between the inner crust and liquid core of hot neutron stars. We find that their values are larger in neutrino-trapped neutron stars than in neutrino-free neutron stars. Furthermore, both are found to decrease with increasing temperature of a neutron star as well as increasing slope parameter of the nuclear symmetry energy, except that the transition pressure in neutrino-trapped neutron stars for the case of small symmetry energy slope parameter first increases and then decreases with increasing temperature. We have also studied the effect of the nuclear symmetry energy on the critical temperature above which the inner crust in a hot neutron star disappears and found that with increasing value of the symmetry energy slope parameter, the critical temperature decreases slightly in neutrino-trapped neutron stars but first decreases and then increases in neutrino-free neutron stars.

Using the Hugenholtz-Van Hove theorem, we derive general expressions for the quadratic and quartic symmetry energies in terms of single-nucleon potentials in isospin asymmetric nuclear matter. These analytical relations are useful for gaining deeper insights into the microscopic origins of the uncertainties in our knowledge on nuclear symmetry energies especially at supra-saturation densities. As examples, the formalism is applied to two model single-nucleon potentials that are widely used in transport model simulations of heavy-ion reactions.

The broad Balmer emission-line profiles resulting from the reprocessing of UV/X-ray radiation from a warped accretion disc induced by the Bardeen-Petterson effect are studied. We adopt a thin warped disc geometry and a central ring-like illuminating source in our model. We compute the steady-state shape of the warped disc numerically, and then use it in the calculation of the line profile. We find that, from the outer radius to the inner radius of the disc, the warp is twisted by an angle of $\sim\pi$ before being flattened efficiently into the equatorial plane. The profiles obtained depend weakly on the illuminating source radius in the range from $3r_{g}$ to $10r_g$, but depend strongly on this radius when it approaches the marginally stable orbit of an extreme Kerr black hole. Double- or triplet-peaked line profiles are present in most cases when the illuminating source radius is low. The triplet-peaked line profiles observed from the Sloan Digital Sky Survey may be a {"}signature" of a warped disc.

The observed radio-optical-X-ray spectral energy distributions (SEDs) of 22 hot spots and 45 knots in the jets of 35 active galactic nuclei are complied from literature and modeled with single-zone lepton models. It is found that the observed luminosities at 5 GHz (L_5GHz) and at 1 keV (L_1keV) are tightly correlated, and the two kinds of sources can be roughly separated with a division of L_1keV=L_5GHz. Our SED fits show that the mechanisms of the X-rays are diverse. Considering the sources at rest, the synchrotron-self-Compton (SSC) scattering would dominate the IC process. This model can interpret the X-rays of some hot spots with a magnetic field strength (B_ssc^delta=1) being consistent with the equipartition magnetic field (B_eq^delta=1) in one order of magnitude, but an unreasonably low B_ssc^delta=1 is required to model the X-rays for all knots. The ratio R_B=B_eq^delta=1/B_ssc^delta=1 is greater than 1 and it is tightly anti-correlated with ratio R_L= L_1keV/L_5GHz for both the knots and the hot spots. We propose that the deviation may be due to the neglect of the relativistic bulk motion for these sources. Considering this effect, we show that the IC/CMB model well explains the X-ray emission for most sources. Both B_eq' and delta are tentatively correlated with R_L. Corrected by the beaming effect, the L'_5GHz-L'_1keV relations for the two kinds of sources are even tighter than the observed ones. These facts suggest that, under the equipartition condition, the observational differences of the X-rays from the knots and hot spots may be mainly due the differences on the Doppler boosting effect and the co-moving magnetic field of the two kinds of sources. Our IC scattering models predict a prominent GeV-TeV component in the SEDs for some sources, which are detectable with H.E.S.S. and Fermi/LAT.

Based on the joint-observations of the radio broadband spectral emissions of solar eclipse on August 1, 2008 at Jiuquan (total eclipse) and Huairou (partial eclipse) at the frequencies of 2.00 -- 5.60 GHz (Jiuquan), 2.60 -- 3.80 GHZ (Chinese solar broadband radiospectrometer, SBRS/Huairou), and 5.20 -- 7.60 GHz (SBRS/Huairou), the authors assemble a successive series of broadband spectrum with a frequency of 2.60 -- 7.60 GHz to observe the solar eclipse synchronously. This is the first attempt to analyze the solar eclipse radio emission under the two telescopes located at different places with broadband frequencies in the periods of total and partial eclipse. With these analyses, the authors made a new semiempirical model of the coronal plasma density of the quiet Sun and made a comparison with the classic models.

Using the nuclear symmetry energy that has been recently constrained by the isospin diffusion data in intermediate-energy heavy ion collisions, we have studied the transition density and pressure at the inner edge of neutron star crusts, and they are found to be 0.040 fm$^{-3}$ $\leq \rho_{t}\leq 0.065$ fm$^{-3}$ and 0.01 MeV/fm$^{3}$ $\leq P_{t}\leq 0.26$ MeV/fm$^{3}$, respectively, in both the dynamical and thermodynamical approaches. We have also found that the widely used parabolic approximation to the equation of state of asymmetric nuclear matter gives significantly higher values of core-crust transition density and pressure, especially for stiff symmetry energies. With these newly determined transition density and pressure, we have obtained an improved relation between the mass and radius of neutron stars.

Considering the non-Newtonian gravity proposed in the grand unification theories, we show that the stability and observed global properties of neutron stars can not rule out the super-soft nuclear symmetry energies at supra-saturation densities. The degree of possible violation of the Inverse-Square-Law of gravity in neutron stars is estimated using an Equation of State (EOS) of neutron-rich nuclear matter consistent with the available terrestrial laboratory data.

We show that the widely used parabolic approximation to the Equation of State (EOS) of asymmetric nuclear matter leads systematically to significantly higher core-crust transition densities and pressures. Using an EOS for neutron-rich nuclear matter constrained by the isospin diffusion data from heavy-ion reactions in the same sub-saturation density range as the neutron star crust, the density and pressure at the inner edge separating the liquid core from the solid crust of neutron stars are determined to be 0.040 fm$^{-3}$ $\leq \rho_{t}\leq 0.065$ fm$^{-3}$ and 0.01 MeV/fm$^{3}$ $\leq P_{t}\leq 0.26$ MeV/fm$^{3}$, respectively. Implications of these constraints on the Vela pulsar are discussed.

Analytical expressions for the saturation density as well as the binding energy and incompressibility at the saturation density of asymmetric nuclear matter are given exactly up to 4th-order in the isospin asymmetry delta =(rho_n - rho_p)/rho using 11 characteristic parameters defined at the normal nuclear density rho_0. Using an isospin- and momentum-dependent modified Gogny (MDI) interaction and the SHF approach with 63 popular Skyrme interactions, we have systematically studied the isospin dependence of the saturation properties of asymmetric nuclear matter, particularly the incompressibility $K_{sat}(\delta )=K_{0}+K_{sat,2}\delta ^{2}+K_{sat,4}\delta ^{4}+O(\delta ^{6})$ at the saturation density. Our results show that the magnitude of the high-order $K_{sat,4}$ parameter is generally small compared to that of the $K_{\sat,2}$ parameter. The latter essentially characterizes the isospin dependence of the incompressibility at the saturation density and can be expressed as $K_{sat,2}=K_{sym}-6L-\frac{J_{0}}{K_{0}}L$, Furthermore, we have constructed a phenomenological modified Skyrme-like (MSL) model which can reasonably describe the general properties of symmetric nuclear matter and the symmetry energy predicted by both the MDI model and the SHF approach. The results indicate that the high-order $J_{0}$ contribution to $K_{sat,2}$ generally cannot be neglected. In addition, it is found that there exists a nicely linear correlation between $K_{sym}$ and $L$ as well as between $J_{0}/K_{0}$ and $K_{0}$. These correlations together with the empirical constraints on $K_{0}$, $L$, $E_{sym}(\rho_{0})$ and the nucleon effective mass lead to an estimate of $K_{sat,2}=-370\pm 120$ MeV.

Pictor A is a nearby Fanaroff-Riley class II (FR II) radio galaxy with a bright hot spot, the western hot spot. Observation of high polarization in the optical emission of the hot spot indicates that the optical emission could be synchrotron radiation of relativistic electrons in the hot spot. These electrons may be able to produce high energy gamma-ray photons through inverse Compton (IC) scattering. We use single-zone and multi-zone synchrotron + synchrotron-self-Compton (SSC) models to fit the observed spectral energy distribution (SED) from the radio to the X-ray band of the hot spot. Our esults show that in the case of a much weaker magnetic field strength than the equipartition magnetic field, both the single-zone and multi-zone models can fit the SED, but the multi-zone model significantly improves the fit. The two models predict the hot spot as a GeV-TeV source, which might be marginally detectable with Fermi/LAT and HESS. The inverse Compton scattering of cosmic microwave background (IC/CMB) is also considered, but its contribution to GeV-TeV emission is negligible. Note that under the equipartition condition, the SED can also be fit with the multi-zone model, but the predicted flux at 10^22 Hz is too weak to be detectable. The detection of TeV gamma-rays from this FR II radio galaxy, if confirmed, would establish a new subclass of extragalactic source in this energy regime since most of the AGNs detected to date at TeV energies are high-energy-peaked BL Lac objects.

We consider a thin accretion disc warped due to the Bardeen-Petterson effect, presenting both analytical and numerical solutions for the situation that the two viscosity coefficients vary with radius as power law, with the two power law indices not necessarily equal. The analytical solutions are compared with numerical ones, showing that our new analytical solution is more accurate than previous one, which overestimates the inclination changing in the outer disc. Our new analytical solution is appropriate for moderately warped discs, while for extremely misaligned disc, only numerical solution is appropriate.

Within both dynamical and thermodynamical approaches using the equation of state for neutron-rich nuclear matter constrained by the recent isospin diffusion data from heavy-ion reactions in the same sub-saturation density range as the neutron star crust, the density and pressure at the inner edge separating the liquid core from the solid crust of neutron stars are determined to be 0.040 fm$^{-3}$ $\leq \rho_{t}\leq 0.065$ fm$^{-3}$ and 0.01 MeV/fm$^{3}$ $\leq P_{t}\leq 0.26$ MeV/fm$^{3}$, respectively. These together with the observed minimum crustal fraction of the total moment of inertia allow us to set a new limit for the radius of the Vela pulsar significantly different from the previous estimate. It is further shown that the widely used parabolic approximation to the equation of state of asymmetric nuclear matter leads systematically to significantly higher core-crust transition densities and pressures, especially with stiffer symmetry energy functionals.

The transition density $\rho_{t}$ and pressure $P_{t}$ at the inner edge separating the liquid core from the solid crust of neutron stars are systematically studied using a modified Gogny (MDI) and 47 popular Skyrme interactions within well established dynamical and thermodynamical methods. It is shown that the widely used parabolic approximation to the full Equation of State (EOS) of isospin asymmetric nuclear matter may lead to huge errors in estimating the \rho_{t} and P_{t}, especially for stiffer symmetry energy functionals $E_{sym}(\rho)$. The \rho_{t} and P_{t} decrease roughly linearly with the increasing slope parameter $L$ of the $E_{sym}(\rho)$ using the full EOS within both methods. It is also shown that the thickness, fractional mass and moment of inertia of neutron star crust are all very sensitive to the parameter $L$ through the $\rho_{t}$. Moreover, it is shown that the $E_{sym}(\rho)$ constrained in the same sub-saturation density range as the neutron star crust by the isospin diffusion data in heavy-ion collisions at intermediate energies limits the transition density and pressure to 0.040$ fm^-3}< \rho_{t} < 0.065$ fm^-3 and 0.01 MeV/fm^3 < P_{t} < 0.26$ MeV/fm^3, respectively. These constrained values for the transition density and pressure are significantly lower than their fiducial values currently used in the literature. Furthermore, the mass-radius relation and several other properties closely related to the neutron star crust are studied by using the MDI interaction. It is found that the newly constrained $\rho_t$ and $P_t$ together with the earlier estimate of $\Delta I/I>0.014$ for the crustal fraction of the moment of inertia of the Vela pulsar impose a stringent constraint of R>= 4.7+4.0M/M_sun km for the radius $R$ and mass $M$ of neutron stars.

Heavy-ion reactions especially those induced by radioactive beams provide useful information about the density dependence of the nuclear symmetry energy, thus the Equation of State of neutron-rich nuclear matter, relevant for many astrophysical studies. The latest developments in constraining the symmetry energy at both sub- and supra-saturation densities from analyses of the isopsin diffusion and the $\pi^-/\pi^+$ ratio in heavy-ion collisions using the IBUU04 transport model are discussed. Astrophysical ramifications of the partially constrained symmetry energy on properties of neutron star crusts, gravitational waves emitted by deformed pulsars and the w-mode oscillations of neutron stars are presented briefly.

Within an isospin- and momentum-dependent hadronic transport model it is shown that the recent FOPI data on the $\pi^-/\pi^+$ ratio in central heavy-ion collisions at SIS/GSI energies (Willy Reisdorf {\it et al.}, NPA {\bf 781}, 459 (2007)) provide circumstantial evidence suggesting a rather soft nuclear symmetry energy \esym at $\rho\geq 2\rho_0$ compared to the Akmal-Pandharipande-Ravenhall prediction. Some astrophysical implications and the need for further experimental confirmations are discussed.

We performed a Cr-K emission line survey in young supernova remnants (SNRs) with the Chandra archival data. Our sample includes W49B, Cas A, Tycho and Kepler. We confirmed the existence of the Cr line in W49B and discovered this emission line in the other three SNRs. The line center energies, equivalent widths (EWs) and fluxes of the Cr lines are given. The Cr in Cas A is in a high ionization state while that in Tycho and Kepler is in a much lower one. We find a good positive correlation between Cr and Fe line center energies, suggesting a common origin of Cr and Fe in the nucleosynthesis, which is consistent with the theoretical predictions. We propose that the EW ratio between Cr and Fe can be used as a supplementary constraint on the progenitors' properties and the explosion mechanism.

By adopting the chemical evolution model of the Milky Way disk, we have studied the star formation and chemical evolution history for M31 galaxy disk. We mainly concentrated on the global properties of the M31 disk. The model has been scaled to the related disk parameters, mainly the disk scale length and total disk baryon mass of M31, which we have adopted to be $r_d$=5.5kpc and $M_{tot} = 7\times 10^{10}$ \ms. It is found that, when the classical Kennicutt star formation law was applied, the obtained radial profiles of gas surface density and star formation rate (SFR) have great difference from the observed results in M31 disk. Then we have adopted modified SFR as we did for the Milky Way galaxy, that is the SFR is radial dependent. Detailed calculations show that by adjusting the star formation efficiency, it is possible to get reasonable gas and abundance profiles, but the total disk SFR is a factor of 2-3 higher than that estimated from observations. And also the predicted SFR radial profile is also much higher than what GALEX observed in the outer part. Possible reasons could be that the M31 disk has been interacted by other factors which seriously altered the star formation history, or the observed SFR is underestimated due to inappropriate extinction correction.

Recent analyses of several isospin effects in heavy-ion reactions have allowed us to constrain the density dependence of nuclear symmetry energy at sub-saturation densities within a narrow range. Combined with constraints on the Equation of State (EOS) of symmetric nuclear matter obtained previously from analyzing the elliptic flow in relativistic heavy-ion collisions, the EOS of neutron-rich nuclear matter is thus partially constrained. Here we report effects of the partially constrained EOS of neutron-rich nuclear matter on the mass-radius correlation, moment of inertia, elliptical deformation and gravitational radiation of (rapidly) rotating neutron stars.

About 10TeV gamma-ray emission within 10 pc region from the Galactic Center had been reported by 4 independent groups. Considering that this TeV gamma-ray emission is produced via a hadronic model, and the relativistic protons came from the tidal disruption of stars by massive black holes, we investigate the spectral nature of the injected relativistic protons required by the hadronic model. The calculation was carried on the tidal disruption of the different types of stars and the different propagation mechanisms of protons in the interstellar medium. Compared with the observation data from HESS, we find for the best fitting that the power-law index of the spectrum of the injected protons is about -1.9, when a red giant star is tidally disrupted, and the effective confinement of protons diffusion mechanism is adopted.

The ultimate goal of studying isospin physics via heavy-ion reactions with neutron-rich, stable and/or radioactive nuclei is to explore the isospin dependence of in-medium nuclear effective interactions and the equation of state of neutron-rich nuclear matter, particularly the isospin-dependent term in the equation of state, i.e., the density dependence of the symmetry energy. Because of its great importance for understanding many phenomena in both nuclear physics and astrophysics, the study of the density dependence of the nuclear symmetry energy has been the main focus of the intermediate-energy heavy-ion physics community during the last decade, and significant progress has been achieved both experimentally and theoretically. In particular, a number of phenomena or observables have been identified as sensitive probes to the density dependence of the nuclear symmetry energy. Experimental studies have confirmed some of these interesting isospin-dependent effects and allowed us to constrain relatively stringently the symmetry energy at sub-saturation densities. The impacts of this constrained density dependence of the symmetry energy on the properties of neutron stars have also been studied, and they were found to be very useful for the astrophysical community. With new opportunities provided by the various radioactive beam facilities being constructed around the world, the study of isospin physics is expected to remain one of the forefront research areas in nuclear physics. In this report, we review the major progress achieved during the last decade in isospin physics with heavy ion reactions and discuss future challenges to the most important issues in this field.

We present high spatial resolution X-ray spectroscopy of supernova remnant Cassiopeia A with the {\sl Chandra} observations. The X-ray emitting region of this remnant was divided into 38 $\times$ 34 pixels with a scale of 10$\arcsec$ $\times$ 10$\arcsec$ each. Spectra of 960 pixels were created and fitted with an absorbed two component non-equilibrium ionization model. With the spectral analysis results we obtained maps of absorbing column density, temperatures, ionization ages, and the abundances for Ne, Mg, Si, S, Ca and Fe. The Si, S and possibly Ca abundance maps show obviously jet structures, while Fe doesn't follow the jet but seems to be distributed perpendicular to it. In the range of about two orders of magnitude, the abundances of Si, S and Ca show tight correlations between each other, suggesting them to be ejecta from explosive O-burning and incomplete Si-burning. Meanwhile, Ne abundance is well correlated with that of Mg, indicating them to be the ashes of explosive C/Ne burning. The Fe abundance is positively correlated with that of Si when Si abundance is lower than 3 solar abundances, but a negative correlation appears when the Si abundance is higher. We suggest that such a two phase correlation is the results of different ways in which Fe is synthesized.

We review our work on Galactic open clusters in recent years, and introduce our proposed large program for the LOCS (LAMOST Open Cluster Survey). First, based on the most complete open clusters sample with metallicity, age and distance data as well as kinematic information, some preliminary statistical analysis regarding the spatial and metallicity distributions is presented. In particular, a radial abundance gradient of - 0.058$\pm$ 0.006 dex kpc$^{-1}$ was derived, and by dividing clusters into age groups we show that the disk abundance gradient was steeper in the past. Secondly, proper motions, membership probabilities, and velocity dispersions of stars in the regions of two very young open clusters are derived. Both clusters show clear evidence of mass segregation, which provides support for the ``primordial'' mass segregation scenarios. Based on the great advantages of the forthcoming LAMOST facility, we have proposed a detailed open cluster survey with LAMOST (the LOCS). The aim, feasibility, and the present development of the LOCS are briefly summarized.

Using various relativistic mean-field models, including the nonlinear ones with meson field self-interactions, those with density-dependent meson-nucleon couplings, and the point-coupling models without meson fields, we have studied the isospin-dependent bulk and single-particle properties of asymmetric nuclear matter. In particular, we have determined the density dependence of nuclear symmetry energy from these different relativistic mean-field models and compare the results with the constraints recently extracted from analyses of experimental data on isospin diffusion and isotopic scaling in intermediate-energy heavy ion collisions as well as from measured isotopic dependence of the giant monopole resonances in even-A Sn isotopes. Among the 23 parameter sets in the relativistic mean-filed model that are commonly used for nuclear structure studies, only a few are found to give symmetry energies that are consistent with the empirical constraints. We have also studied the nuclear symmetry potential and the isospin-splitting of the nucleon effective mass in isospin asymmetric nuclear matter. We find that both the momentum dependence of the nuclear symmetry potential at fixed baryon density and the isospin-splitting of the nucleon effective mass in neutron-rich nuclear matter depend not only on the nuclear interactions but also on the definition of the nucleon optical potential.

We derive the proper motions, membership probabilities, and velocity dispersions of stars in the regions of the young (about 2-4 Myr-old) open clusters NGC 2244 (the central cluster in the Monoceros R2 association) and NGC 6530 (the dominant cluster in the Sgr OB1 association) from photographic plate material obtained at Shanghai Astronomical Observatory, with time baselines of 34 and 87 years, respectively. Both clusters show clear evidence of mass segregation, but they do not exhibit any significant velocity-mass (or, equivalently, a velocity-luminosity) dependence. This provides strong support for the suggestion that the observed mass segregation is -- at least partially -- due to the way in which star formation has proceeded in these complex star-forming regions (``primordial'' mass segregation). Based on arguments related to the clusters' published initial mass functions, in conjunction with our new measurements of their internal velocity dispersions (35 and 8 km/s for NGC 2244 and NGC 6530, respectively), we provide strong arguments in favor of the dissolution of NGC 2244 on very short time-scales, while we speculate that NGC 6530 may be more stable against the effects of internal two-body relaxation. However, this latter object may well be destroyed by the strong tidal field prevalent at its location in the Galactic plane in the direction of the Galactic Center.

The liquid-gas phase transition in hot neutron-rich nuclear matter is investigated within a self-consistent thermal model using an isospin and momentum dependent interaction (MDI) constrained by the isospin diffusion data in heavy-ion collisions, a momentum-independent interaction (MID), and an isoscalar momentum-dependent interaction (eMDYI). The boundary of the phase-coexistence region is shown to be sensitive to the density dependence of the nuclear symmetry energy with a softer symmetry energy giving a higher critical pressure and a larger area of phase-coexistence region. Compared with the momentum-independent MID interaction, the isospin and momentum-dependent MDI interaction is found to increase the critical pressure and enlarge the area of phase-coexistence region. For the isoscalar momentum-dependent eMDYI interaction, a limiting pressure above which the liquid-gas phase transition cannot take place has been found and it is shown to be sensitive to the stiffness of the symmetry energy.

Heavy-ion reactions provide a unique means to investigate the equation of state (EOS) of neutron-rich nuclear matter, especially the density dependence of the nuclear symmetry energy $E_{sym}(\rho)$. The latter plays an important role in understanding many key issues in both nuclear physics and astrophysics. Recent analyses of heavy-ion reactions have already put a stringent constraint on the $E_{sym}(\rho)$ around the saturation density. This subsequently allowed us to constrain significantly the radii and cooling mechanisms of neutron stars as well as the possible changing rate of the gravitational constant G.

Heavy-ion reactions induced by neutron-rich nuclei provide a unique means to investigate the equation of state of isospin-asymmetric nuclear matter, especially the density dependence of the nuclear symmetry energy. In particular, recent analyses of the isospin diffusion data in heavy-ion reactions have already put a stringent constraint on the nuclear symmetry energy around the nuclear matter saturation density. We review this exciting result and discuss its implications on nuclear effective interactions and the neutron skin thickness of heavy nuclei. In addition, we also review the theoretical progress on probing the high density behaviors of the nuclear symmetry energy in heavy-ion reactions induced by high energy radioactive beams.

Within a self-consistent thermal model using an isospin and momentum dependent interaction (MDI) constrained by the isospin diffusion data in heavy-ion collisions, we investigate the temperature dependence of the symmetry energy $E_{sym}(\rho, T)$ and symmetry free energy $F_{sym}(\rho, T) $ for hot, isospin asymmetric nuclear matter. It is shown that the symmetry energy $E_{sym}(\rho, T)$ generally decreases with increasing temperature while the symmetry free energy $F_{sym}(\rho, T)$ exhibits opposite temperature dependence. The decrement of the symmetry energy with temperature is essentially due to the decrement of the potential energy part of the symmetry energy with temperature. The difference between the symmetry energy and symmetry free energy is found to be quite small around the saturation density of nuclear matter. While at very low densities, they differ significantly from each other. In comparison with the experimental data of temperature dependent symmetry energy extracted from the isotopic scaling analysis of intermediate mass fragments (IMF's) in heavy-ion collisions, the resulting density and temperature dependent symmetry energy $E_{sym}(\rho, T) $ is then used to estimate the average freeze-out density of the IMF's.used to estimate the average freeze-out density of the IMF's.

Using the relativistic impulse approximation with the Love-Franey \textsl{NN} scattering amplitude developed by Murdock and Horowitz, we investigate the low-energy (100 MeV$\leq E_{\mathrm{kin}}\leq 400$ MeV) behavior of the nucleon Dirac optical potential, the Schr\"{o}dinger-equivalent potential, and the nuclear symmetry potential in isospin asymmetric nuclear matter. We find that the nuclear symmetry potential at fixed baryon density decreases with increasing nucleon energy. In particular, the nuclear symmetry potential at saturation density changes from positive to negative values at nucleon kinetic energy of about 200 MeV. Furthermore,the obtained energy and density dependence of the nuclear symmetry potential is consistent with those of the isospin- and momentum-dependent MDI interaction with $x=0$, which has been found to describe reasonably both the isospin diffusion data from heavy-ion collisions and the empirical neutron-skin thickness of $^{208} $Pb.

Recent analysis of the isospin diffusion data from heavy-ion collisions based on an isospin- and momentum-dependent transport model with in-medium nucleon-nucleon cross sections has led to the extraction of a value of $L=88\pm 25$ MeV for the slope of the nuclear symmetry energy at saturation density. This imposes stringent constraints on both the parameters in the Skyrme effective interactions and the neutron skin thickness of heavy nuclei. Among the 21 sets of Skyrme interactions commonly used in nuclear structure studies, the 4 sets SIV, SV, G$_\sigma$, and R$_\sigma$ are found to give $L$ values that are consistent with the extracted one. Further study on the correlations between the thickness of the neutron skin in finite nuclei and the nuclear matter symmetry energy in the Skyrme Hartree-Fock approach leads to predicted thickness of the neutron skin of $0.22\pm 0.04$ fm for $^{208}$Pb, $0.29\pm 0.04$ fm for $^{132}$Sn, and $0.22\pm 0.04$ fm for $^{124}$Sn.

Nuclear reactions induced by stable and/or radioactive neutron-rich nuclei provide the opportunity to pin down the equation of state of neutron-rich matter, especially the density ($\rho$) dependence of its isospin-dependent part, i.e., the nuclear symmetry energy $E_{\rm sym}$. A conservative constraint, $32(\rho /\rho_{0})^{0.7} < E_{\rm sym}(\rho ) < 32(\rho /\rho _{0})^{1.1}$, around the nuclear matter saturation density $\rho_0$ has recently been obtained from the isospin diffusion data in intermediate energy heavy-ion collisions. We review this exciting result and discuss its consequences and implications on nuclear effective interactions, radii and cooling mechanisms of neutron stars.

Nucleon-nucleon (NN) cross sections are evaluated in neutron-rich matter using a scaling model according to nucleon effective masses. It is found that the in-medium NN cross sections are not only reduced but also have a different isospin dependence compared with the free-space ones. Because of the neutron-proton effective mass splitting the difference between nn and pp scattering cross sections increases with the increasing isospin asymmetry of the medium. Within the transport model IBUU04, the in-medium NN cross sections are found to influence significantly the isospin transport in heavy-ion reactions. With the in-medium NN cross sections, a symmetry energy of $E_{sym}(\rho)\approx 31.6(\rho /\rho_{0})^{0.69}$ was found most acceptable compared with both the MSU isospin diffusion data and the presently acceptable neutron-skin thickness in $^{208}$Pb. The isospin dependent part $K_{asy}(\rho _{0})$ of isobaric nuclear incompressibility was further narrowed down to $-500\pm 50$ MeV. The possibility of determining simultaneously the in-medium NN cross sections and the symmetry energy was also studied. The proton transverse flow, or even better the combined transverse flow of neutrons and protons, can be used as a probe of the in-medium NN cross sections without much hindrance from the uncertainties of the symmetry energy.

The latest development in determining the density dependence of the nuclear symmetry energy using heavy-ion collisions is reviewed. Within the IBUU04 version of an isospin- and momentum-dependent transport model using a modified Gogny effective interaction, recent experimental data from NSCL/MSU on isospin diffusion are found to be consistent with a nuclear symmetry energy of $E_{sym}(\rho)\approx 31.6(\rho /\rho_{0})^{1.05}$ at subnormal densities. Predictions on several observables sensitive to the density dependence of the symmetry energy at supranormal densities accessible at GSI and the planned Rare Isotope Accelerator (RIA) are also made.

Large-voltage and high-temperature plasma columns produced by pinch discharge can generate gamma-ray flashes with energy spectra and spectral evolution consistent with that observed in gamma-ray bursts (GRBs). The inverse Compton scattering (ICS) during the discharge process can produce high linear polarization. The calculation indicates that the observed polarization depends on the angle between the observer's line-of-sight to the GRB and the direction of the pinch discharge, but only weakly depends on observed gamma-ray energy.

With an isospin- and momentum-dependent transport model, we find that the degree of isospin diffusion in heavy ion collisions at intermediate energies is affected by both the stiffness of the nuclear symmetry energy and the momentum dependence of the nucleon potential. Using a momentum dependence derived from the Gogny effective interaction, recent experimental data from NSCL/MSU on isospin diffusion are shown to be consistent with a nuclear symmetry energy given by $E_{\text{sym}}(\rho)\approx 31.6(\rho /\rho_{0})^{1.05} $ at subnormal densities. This leads to a significantly constrained value of about -550 MeV for the isospin-dependent part of the isobaric incompressibility of isospin asymmetric nuclear matter.

Using the data acquired in the Time To Spill (TTS) mode for long gamma-ray bursts (GRBs) collected by the Burst and Transient Source Experiment on board the Compton Gamma Ray Observatory (BATSE/CGRO), we have carefully measured spectral lags in time between the low (25-55 keV) and high (110-320 keV) energy bands of individual pulses contained in 64 multi-peak GRBs. We find that the temporal lead by higher-energy gamma-ray photons (i.e., positive lags) is the norm in this selected sample set of long GRBs. While relatively few in number, some pulses of several long GRBs do show negative lags. This distribution of spectral lags in long GRBs is in contrast to that in short GRBs. This apparent difference poses challenges and constraints on the physical mechanism(s) of producing long and short GRBs. The relation between the pulse peak count rates and the spectral lags is also examined. Observationally, there seems to be no clear evidence for systematic spectral lag-luminosity connection for pulses within a given long GRB.

A technique for timescale analysis of spectral lags performed directly in the time domain is developed. Simulation studies are made to compare the time domain technique with the Fourier frequency analysis for spectral time lags. The time domain technique is applied to studying rapid variabilities of X-ray binaries and $\gamma$-ray bursts. The results indicate that in comparison with the Fourier analysis the timescale analysis technique is more powerful for the study of spectral lags in rapid variabilities on short time scales and short duration flaring phenomena.

Slim disks have been received much attention because of the increasing evidence for the super-critical accreting objects. In this paper, we make an attempt to construct a unified model, in which the viscosity and the dimensionless accretion rate can span rather wide ranges. We replace blackbody radiation under diffusion approximation with a bridged formula, which accounts for both blackbody radiation and thermal bremsstrahlung in optically-thick and -thin cases, respectively. Thus this allows us to investigate the structures of and the emergent spectra from slim disks in a wider parameter space, covering transition regions from optically thick to optically thin. We show that there is a maximum transition radius, roughly $R_{\rm tr, max}/R_g\sim 50$ when $\dot{M}/\dot{M}_{\rm C}\sim 15$. The emergent spectra from the unified model of the accretion disk have been calculated. A simple model of hot corona above the slim disk is taken into account for the hard X-ray spectrum in this paper based on Wang & Netzer (2003). We have applied the present model to the microquasar GRS 1915+105, narrow line Seyfert 1 galaxies RE J1034+396 and Akn 564. Our model can explain well the broadband X-ray spectra of narrow line Seyfert 1 galaxies, microquasars and possible ultra-luminous compact X-ray sources. The present model can be widely applied to the candidates of super-critical accreting objects.

We present the results from a detailed deprojection analysis of Abell 1835 as observed by {\it XMM-Newton}. If we fit the spectra with an isothermal plasma model, the deprojected temperature profile is flat in the outer region around 7.6 keV and decreases to $\sim$ 5.6 keV in the center, which may be connected with the gas cooling. In the central part, a two-component thermal plasma model can fit the spectrum significantly better. Moreover, the cool component (T $\sim$ 1.8 keV) has a much lower metal abundance than the hot component (T $\sim$ 8 keV), which may be due to the longer cooling time for the cool gas with lower abundance. In addition, it was found that without a main isothermal component, the standard cooling flow model cannot fit the spectrum satisfactorily. From the isothermal model fitting results we also derived the electron density $n_e$, and fitted its radial distribution with a double-$\beta$ model. The $n_e$ profile inferred with the double-$\beta$ model and the deprojected X-ray gas temperature profile were then combined to derive the total mass and the total projected mass of the cluster. The projected mass is lower than that derived from the weak lensing method. However, assuming that the cluster extends to a larger radius $\sim15'$ as found by Clowe & Schneider (2002), the two results are consistent within the error bars. Furthermore, we calculated the projected mass within the radius of $\sim$ 153 kpc implied by the presence of a gravitational lensing arc, which is about half of the mass determined from the optical lensing.

We investigated the time lags and the evolution of the cross spectra of Z source GX~5-1, observed by the Rossi X-ray Timing Explorer (RXTE), when it is in the horizontal branch oscillations. We showed that the time lags of 3 horizontal branch oscillations are related to the position on the hardness intensity diagram. All of the three QPOs were shown to have hard time lags. However on the cross spectra, one is in a `dip', one in a `bump', the other has no so obvious characteristic. The time lags of two of the QPOs decrease with QPO's frequency, while the other has a trend increasing with its frequency. Moreover, in the normal branch, we found no significant time lags in the present observational data.

The coalescence model based on nucleon distribution functions from an isospin-dependent transport model is used to study the production of light clusters such as deuteron, triton, and $^{3}$He from heavy-ion collisions induced by neutron-rich nuclei at intermediate energies. It is found that the emission time of light clusters depends on their masses. For clusters with the same momentum per nucleon, heavier ones are emitted earlier. Both the yield and energy spectrum of light clusters are sensitive to the density dependence of nuclear symmetry energy, with more light clusters produced in the case of a stiff symmetry energy. On the other hand, effects due to the stiffness of the isoscalar part of nuclear equation of state and the medium dependence of nucleon-nucleon cross sections on light cluster production are unimportant. We have also studied the correlation functions of clusters, and they are affected by the density dependence of nuclear symmetry energy as well, with the stiff symmetry energy giving a stronger anti-correlation of light clusters, particularly for those with large kinetic energies. Dependence of light cluster production on the centrality and incident energy of heavy ion collisions as well as the mass of the reaction system is also investigated.

We have compiled two new Galactic open cluster catalogues. The first one has 119 objects with age, distance and metallicity data, while the second one has 144 objects with both absolute proper motion and radial velocity data, of which 45 clusters also with metallicity data available. An iron radial gradient of about $-0.063\pm$0.008 \dexkpc is obtained from our sample. By dividing clusters into age groups, we show that iron gradient was steeper in the past. A disk age-metallicity relation could very probably exist based on the current sample.

When thermal relativistic electrons with isotropic distribution of velocities move in a gas region, or impinge upon the surface of a cloud that consists of a dense gas or doped dusts, the Cerenkov effect produces peculiar atomic or ionic emission lines -- the Cerenkov line-like radiation. This newly recognized emission mechanism may find wide applications in high-energy astrophysics. In this paper, we tentatively adopt this new line emission mechanism to discuss the origin of iron K$_{\alpha}$ feature of AGNs. Motivation of this research is to attempt a solution to a problem encountered by the ``disk-fluorescence line'' model, i.e. the lack of temporal response of the observed iron K$_{\alpha}$ line flux to the changes of the X-ray continuum flux. If the Cerenkov line emission is indeed responsible significantly for the iron K$_{\alpha}$ feature, the conventional scenario around the central supermassive black holes of AGNs would need to be modified to accommodate more energetic, more violent and much denser environments than previously thought.

Using an isospin-dependent transport model for heavy-ion collisions induced by neutron-rich nuclei at intermediate energies, we study the production of light clusters such as deuteron, triton, and $^{3}$He via coalescence of nucleons. We find that both the yields and energy spectra of these light clusters are affected significantly by the density dependence of nuclear symmetry energy, with a stiffer symmetry energy giving a larger yield.

Using an isospin-dependent transport model, we study isospin effects on two-nucleon correlation functions in heavy-ion collisions induced by neutron-rich nuclei at intermediate energies. We find that these correlation functions are sensitive to the density dependence of nuclear symmetry energy, but not to the incompressibility of symmetric nuclear matter and the medium dependence of nucleon-nucleon cross sections. This sensitivity is mainly due to effects of nuclear symmetry energy on the emission times of neutrons and protons as well as their relative emission sequence. We also study the variations of the symmetry energy effects on nucleon-nucleon correlations with respect to the impact parameter, incident energy, and mass number of heavy ion collisions.

Using an isospin-dependent transport model, we study the effects of nuclear symmetry energy on two-nucleon correlation functions in heavy ion collisions induced by neutron-rich nuclei. We find that the density dependence of the nuclear symmetry energy affects significantly the nucleon emission times in these collisions, leading to larger values of two-nucleon correlation functions for a symmetry energy that has a stronger density dependence. Two-nucleon correlation functions are thus useful tools for extracting information about the nuclear symmetry energy from heavy ion collisions.