The Weakly Interacting Massive Particle (WIMP) has been one of the most attractive candidates for Dark Matter (DM), and the lightest neutralino ($\widetilde{\chi}^0_1$) in the Next-to-Minimal Supersymmetric Standard Model (NMSSM) is an interesting realization of WIMP. The Galactic Center Excess (GCE) can be explained by WIMP DM annihilations in the sky. In this work we consider the $Z_3$-NMSSM where the singlet $S$ and Singlino $\widetilde{S}^0$ components play important roles in the Higgs and DM sector. Guided by our analytical arguments, we perform a numerical scan over the NMSSM parameter space for the GCE explanation by considering various observables such as the Standard Model (SM) Higgs data measured by the ATLAS and CMS experiments, and the $B$-physics observables $BR(B_s\rightarrow X_s\gamma)$ and $BR(B_s\rightarrow \mu^+\mu^-)$. We find that the correlation between the coupling $C_{A_1 b\bar{b}}$ in $\langle \sigma_{b\bar{b}} v \rangle _{0}$ and the coupling $C_{Z \widetilde{\chi}^0_1 \widetilde{\chi}^0_1}$ in DM-neutron Spin Dependent (SD) scattering rate $\sigma^{SD}_{\widetilde{\chi}^0_1-N}$ makes all samples we obtain for GCE explanation get excluded by the PandaX-II results. Although the DM resonant annihilation scenarios may be beyond the reach of our analytical approximations and scan strategy, the aforementioned correlation can be a reasonable motivation for future experiments such as PandaX-nT to further test the NMSSM interpretation of GCE.

Left-Right symmetric model (LRSM) has been an attractive extension of the Standard Model (SM) which can address the origin of parity violation in the SM electroweak (EW) interactions, generate tiny neutrino masses, accommodate dark matter (DM) candidates and provide a natural framework for baryogenesis through leptogenesis. In this work we utilize the minimal LRSM to study the recently reported DAMPE results of cosmic $e^+e^-$ spectrum which exhibits a tentative peak around 1.4 TeV, while satisfying the current neutrino data. We propose to explain the DAMPE peak with a complex scalar DM $\chi$ in two scenarios: 1) $\chi\chi^* \to H_1^{++}H_1^{--} \to \ell_i^+\ell_i^+\ell_j^-\ell_j^-$; 2) $\chi\chi^* \to H_{k}^{++}H_{k}^{--} \to \ell_i^+\ell_i^+\ell_j^-\ell_j^-$ accompanied by $\chi\chi^* \to H_1^+ H_1^- \to \ell_i^+ \nu_{\ell_i} \ell_j^- \nu_{\ell_j}$ with $\ell_{i,j}=e,\mu,\tau$ and $k=1,2$. We fit the theoretical prediction on $e^+e^-$ spectrum to relevant experimental data to determine the scalar mass spectrum favored by the DAMPE excess. We also consider various constraints from theoretical principles, collider experiments as well as DM relic density and direct search experiments. We find that there are ample parameter space which can interpret the DAMPE data while passing the constraints. Our explanations, on the other hand, usually imply the existence of other new physics at the energy scale ranging from $10^7 {\rm GeV}$ to $10^{11} {\rm GeV}$. Collider tests of our explanations are also discussed.

Recently, DAMPE experiment released the new measurement of the total cosmic $e^+e^-$ flux between 25 GeV and 4.6 TeV which indicates a spectral softening at around 0.9 TeV and a tentative peak at around 1.4 TeV. We utilize the scalar dark matter (DM) annihilation scenario to explain the DAMPE peak by extending $G_{SM}\equiv SU(3)_C \times SU(2)_L \times U(1)_Y$ with additional $U(1)$ gauge symmetries while keeping anomaly free to generate $\chi \chi \to Z^\prime Z^\prime \to \ell\bar{\ell}\ell^\prime\overline{\ell^\prime}$, where $\chi, Z^\prime, \ell^{(^\prime)}$ denote the scalar DM, the new gauge boson and $\ell^{(^\prime)}=e,\mu,\tau$, respectively, with $m_\chi \sim m_{Z^\prime} \sim 2 \times 1.5$ (TeV). We first illustrate that the minimal framework $G_{SM} \times U(1)_{Y^\prime}$ with the above mass choices can explain the DAMPE excess but has been excluded by LHC constraints from the $Z^\prime$ searches. Then we study a non-minimal framework $G_{SM} \times U(1)_{Y^\prime} \times U(1)_{Y^{\prime \prime}}$ in which $U(1)_{Y^{\prime \prime}}$ mixes with $U(1)_{Y^\prime}$. We show that such a framework can interpret the DAMPE data while passing other constraints including the DM relic abundance, DM direct detection and collider bounds. We also investigate the predicted $e^+e^-$ spectrum in this framework and find that the mass splitting $\Delta m = m_\chi - m_{Z'}$ should be less than about 17 GeV to produce the peak-like structure.

Integrated pulse profiles at 8.6~GHz obtained with the Shanghai Tian Ma Radio Telescope (TMRT) are presented for a sample of 26 pulsars. Mean flux densities and pulse width parameters of these pulsars are estimated. For eleven pulsars these are the first high-frequency observations and for a further four, our observations have a better signal-to-noise ratio than previous observations. For one (PSR J0742-2822) the 8.6~GHz profiles differs from previously observed profiles. A comparison of 19 profiles with those at other frequencies shows that in nine cases the separation between the outmost leading and trailing components decreases with frequency, roughly in agreement with radius-to-frequency mapping, whereas in the other ten the separation is nearly constant. Different spectral indices of profile components lead to the variation of integrated pulse profile shapes with frequency. In seven pulsars with multi-component profiles, the spectral indices of the central components are steeper than those of the outer components. For the 12 pulsars with multi-component profiles in the high-frequency sample, we estimate the core width using gaussian fitting and discuss the width-period relationship.

The emission geometries, e.g. the emission region height, the beam shape, and radius-to-frequency mapping, are important predictions of pulsar radiation model. The multi-band radio observations carry such valuable information. In this paper, we study two bright pulsars, (PSRs B0329+54 and B1642-03) and observe them in high frequency (2.5 GHz, 5 GHz, and 8 GHz). The newly acquired data together with historical archive provide an atlas of multi-frequency profiles spanning from 100 MHz to 10 GHz. We study the frequency evolution of pulse profiles and the radiation regions with the these data. We firstly fit the pulse profiles with Gaussian functions to determine the phase of each component, and then calculate the radiation altitudes of different emission components and the radiation regions. We find that the inverse Compton scattering (ICS) model can reproduce the radiation geometry of these two pulsars. But for PSR B0329+54 the radiation can be generated in either annular gap (AG) or core gap (CG), while the radiation of PSR B1642-03 can only be generated in the CG. This difference is caused by the inclination angle and the impact angle of these two pulsars. The relation of beaming angle (the angle between the radiation direction and the magnetic axis) and the radiation altitudes versus frequency is also presented by modelling the beam-frequency evolution in the ICS model. The multi-band pulse profiles of these two pulsars can be described well by the ICS model combined with the CG and AG.

The pulse profile of PSR B1133+16 is usually regarded as a conal-double structure. However, its multifrequency profiles cannot simply be fitted with two Gaussian functions, and a third component is always needed to fit the bridge region (between two peaks). This would introduce additional, redundant parameters. In this paper, through a comparison of five fitting functions (Gaussian, von Mises, hyperbolic secant, square hyperbolic secant, and Lorentz), it is found that the square hyperbolic secant function can best reproduce the profile, yielding an improved fit. Moreover, a symmetric 2D radiation beam function, instead of a simple 1D Gaussian function, is used to fit the profile. Each profile with either well-resolved or not-so-well-resolved peaks could be fitted adequately using this beam function, and the bridge emission between the two peaks does not need to be a new component. Adopting inclination and impact angles based on polarization measurements, the opening angle ({\theta}_{\mu}0) of the radiation beam in a certain frequency band is derived from beam-function fitting. The corresponding radiation altitudes are then calculated. Based on multi-frequency profiles, we also computed the Lorentz factors of the particles and their dispersion at those locations in both the curvature-radiation (CR) and inverse-Compton-scattering (ICS) models. We found that the Lorentz factors of the particles decrease rapidly as the radiation altitude increases. Besides, the radiation prefers to be generated in annular region rather than core region, and this needs further validation.

We present the results of six-epoch Very Long Baseline Array (VLBA) observations of 3C~66A. The high-resolution Very Long Baseline Interferometer (VLBI) maps obtained at multi-frequency (2.3, 8.4, and 22.2 GHz) simultaneously enabled us to identify the brightest compact component with the core. We find that the spectrum of the core can be reasonably fitted by the synchrotron self-absorption model. Our VLBA maps show that the jet of 3C~66A has two bendings at about 1.2 and 4 mas from the core. We also give possible identifications of our jet components with the components in previous VLBA observations by analysing their proper motions. We find consistent differences of the position from the core in one component between different frequencies at six epochs.

We investigate a phase-reference VLBI observation that was conducted at 15.4 GHz by fast switching VLBA antennas between the compact steep-spectrum radio source 3C 138 and the quasar PKS 0528+134 which are about 4$^\circ$ away on the sky. By comparing the phase-reference mapping with the conventional hybrid mapping, we demonstrate the feasibility of high precision astrometric measurements for sources separated by 4$^\circ$. VLBI phase-reference mapping preserves the relative phase information, and thus provides an accurate relative position between 3C 138 and PKS 0528+134 of $\Delta\alpha=-9^m46^s.531000\pm0^s.000003$ and $\Delta\delta=3^\circ6^\prime26^{\prime\prime}.90311\pm0^{\prime\prime}.00007$ (J2000.0) in right ascension and declination, respectively. This gives an improved position of the nucleus (component A) of 3C 138 in J2000.0 to be RA=$05^h 21^m 9^s.885748$ and Dec=$16^\circ 38' 22''.05261$ under the assumption that the position of calibrator PKS 0528+134 is correct. We further made a hybrid map by performing several iterations of CLEAN and self-calibration on the phase-referenced data with the phase-reference map as an input model for the first phase self-calibration. Compared with the hybrid map from the limited visibility data directly obtained from fringe fitting 3C 138 data, this map has a similar dynamic range, but a higher angular resolution. Therefore, phase-reference technique is not only a means of phase connection, but also a means of increasing phase coherence time allowing self-calibration technique to be applied to much weaker sources.

We present results of Very Long Baseline Array (VLBA) observations of PKS 0528+134 at five frequencies (2.3, 5.0, 8.4, 15.4, and 22.2 GHz). These quasi-simultaneous data enable us to study the spectral distribution of Very Long Baseline Interferometer (VLBI) components for the first time in this highly variable source, from which the central compact core is identified. Our observations indicate that there are two bendings for the jet motion at parsec scale. We provide an approximate spatial fit to the curved jet trajectory using the Steffen et al. (\cite{Steffen95}) helical model. We further investigate the proper motions of three jet components, which all show superluminal motion. At high frequencies (15.4 and 22.2 GHz) we detected a new component, which is estimated to be related to a radio burst peaking at about 2000.

We present the results from the first quasi-simultaneous multi-frequency (2.3, 5.0, 8.4 and 15 GHz) Very Long Baseline Interferometry (VLBI) observations of a compact steep spectrum (CSS) superluminal source 3C138. For the first time, the spectral distribution of the components within its central 10 milli-arcsecond (mas) region was obtained. This enables us to identify the component at the western end as the location of the nuclear activity on the assumption that the central engine is associated with one of the detected components. The possibility that none of these visible components is the true core is also discussed. The new measurements further clarify the superluminal motions of its inner jet components. The multi-frequency data reveal a convex spectrum in one jet component, inferring the existence of free-free absorption by the ambient dense plasma.