Pulsar's position, proper motion and parallax are important parameters in timing equations. It is a really challenging work to fit astrometric parameters accurately through pulsar timing, especially for pulsars that show irregular timing properties. As the fast development of related techniques, it is possible to measure astrometric parameters of more and more pulsars in a model$\textrm{-}$independent manner with the Very Long Baseline Interferometry (VLBI). In this work, we select 16 normal pulsars, whose parallax and proper motion have not been successfully fitted with timing observations or show obvious differences with corresponding latest VLBI solutions, and do further studies on their timing properties. After updating astrometric parameters in pulsar ephemerides with the latest VLBI measurements, we derive the latest rotation solutions of these pulsars with observation data at S and C$\textrm{-}$band obtained from the Shanghai Tian Ma Radio Telescope (TMRT). Compared with spin frequency $\nu$ inferred from previous rotation solutions, the newly$\textrm{-}$fitted $\nu$ show differences larger than 10$^{-9}$ Hz for most pulsars. The contribution of the Shklovsky effect to period derivative $\dot{P}$ can be properly removed taking advantages of accurate proper motion and distance of target pulsars measured by VLBI astrometry. This further leads to a precise estimate of intrinsic characteristic age $\tau_{\rm c}$. Differences between the newly$\textrm{-}$measured $\tau_{\rm c}$ and corresponding previous results are as large as 2% for some pulsars. VLBI astrometric parameter solutions also lead to better measurements of timing irregularities. For PSR B0154$+$61, the glitch epoch (MJD 58279.5) measured with previous ephemeris is about 13 d later than the result (MJD 58266.4) obtained with VLBI astrometric parameter solutions.

We present a first proof-of-principle study for using deep neural networks (DNNs) as a novel search method for continuous gravitational waves (CWs) from unknown spinning neutron stars. The sensitivity of current wide-parameter-space CW searches is limited by the available computing power, which makes neural networks an interesting alternative to investigate, as they are extremely fast once trained and have recently been shown to rival the sensitivity of matched filtering for black-hole merger signals. We train a convolutional neural network with residual (short-cut) connections and compare its detection power to that of a fully-coherent matched-filtering search using the WEAVE pipeline. As test benchmarks we consider two types of all-sky searches over the frequency range from $20\,\mathrm{Hz}$ to $1000\,\mathrm{Hz}$: an `easy' search using $T=10^5\,\mathrm{s}$ of data, and a `harder' search using $T=10^6\,\mathrm{s}$. Detection probability $p_\mathrm{det}$ is measured on a signal population for which matched filtering achieves $p_\mathrm{det}=90\%$ in Gaussian noise. In the easiest test case ($T=10^5\,\mathrm{s}$ at $20\,\mathrm{Hz}$) the DNN achieves $p_\mathrm{det}\sim88\%$, corresponding to a loss in sensitivity depth of $\sim5\%$ versus coherent matched filtering. However, at higher-frequencies and longer observation time the DNN detection power decreases, until $p_\mathrm{det}\sim13\%$ and a loss of $\sim 66\%$ in sensitivity depth in the hardest case ($T=10^6\,\mathrm{s}$ at $1000\,\mathrm{Hz}$). We study the DNN generalization ability by testing on signals of different frequencies, spindowns and signal strengths than they were trained on. We observe excellent generalization: only five networks, each trained at a different frequency, would be able to cover the whole frequency range of the search.

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.

We present integrated pulse profiles at 5~GHz for 71 pulsars, including eight millisecond pulsars (MSPs), obtained using the Shanghai Tian Ma Radio Telescope (TMRT). Mean flux densities and pulse widths are measured. For 19 normal pulsars and one MSP, these are the first detections at 5~GHz and for a further 19, including five MPSs, the profiles have a better signal-to-noise ratio than previous observations. Mean flux density spectra between 400~MHz and 9~GHz are presented for 27 pulsars and correlations of power-law spectral index are found with characteristic age, radio pseudo-luminosity and spin-down luminosity. Mode changing was detected in five pulsars. The separation between the main pulse and interpulse is shown to be frequency independent for six pulsars but a frequency dependence of the relative intensity of the main pulse and interpulse is found. The frequency dependence of component separations is investigated for 20 pulsars and three groups are found: in seven cases the separation between the outmost leading and trailing components decreases with frequency, roughly in agreement with radius-to-frequency mapping; in eleven cases the separation is nearly constant; in the remain two cases the separation between the outmost components increases with frequency. We obtain the correlations of pulse widths with pulsar period and estimate the core widths of 23 multi-component profiles and conal widths of 17 multi-component profiles at 5.0~GHz using Gaussian fitting and discuss the width-period relationship at 5~GHz compared with the results at at 1.0~GHz and 8.6~GHz.

Using photometric and spectroscopic data of Supernova (SN) 2012cu, a fairly reddened type Ia supernova, we derived its color excess curves and probed the dust extinction in its host galaxy, NGC 4772. In order to derive the extinction as a function of wavelength (i.e., $A_\lambda$), we model the color excess curves of SN 2012cu in terms of dust models consisting of silicate and carbonaceous (graphite or amorphous carbon) dust. The modeled extinction law towards SN 2012cu extends flatly to the far-ultraviolet (UV) bands, which is much flatter than those of the Milky Way and Magellanic Clouds, and the 2175 $\AA$ feature is very weak or absent. The flatness of the modeled extinction curve in the UV bands suggests a "grey" extinction law of the active galactic nucleus in the vicinity of the SN 2012cu-Earth line of sight. Our results indicate that the sizes of the dust in the ISM towards SN 2012cu in NGC 4772 are larger than those of the Milky Way and the Large Magellanic Cloud, and much larger than that of the Small Magellanic Cloud. The best fitting gives an observed visual extinction towards SN 2012cu of $A_V \approx 2.6$ mag, a reddening of $E(B-V) \approx 1.0$ mag, with a total-to-selective extinction ratio $R_V \approx 2.7$, consistent with previous results.

One large glitch was detected in PSR B1737$-$30 using data spanning from MJD 57999 to 58406 obtained with the newly built Shanghai Tian Ma Radio Telescope (TMRT). The glitch took place at the time around MJD 58232.4 when the pulsar underwent an increase in the rotation frequency of $\Delta \nu$ about 1.38$\times 10^{-6}$ Hz, corresponding to a fractional step change of $\Delta \nu / \nu$ $\thicksim$ 8.39$\times 10^{-7}$. Post$\textrm{-}$glitch $\nu$ gradually decreased to the pre$\textrm{-}$glitch value. The frequency derivative was observed to undergo a step change of about $-$9$\times 10^{-16}$ s$^{-2}$. Since July 1987, there are 36 glitches already reported in PSR B1737$-$30 including this one. According to our analysis, the glitch size distribution is well described by the power law with index of 1.13. The distribution of the interval between two adjacent glitches (waiting time $\Delta T$) follows a Poissonian probability density function. For PSR B1737$-$30, the interval is prone to be long after a large glitch. But no correlation is found between glitch size and the interval since previous glitch.

We have investigated the mode-changing properties of PSR B0329+54 using 31 epochs of simultaneous 13 cm/3 cm single-pulse observations obtained with Shanghai Tian Ma 65 m telescope. The pulsar was found in the abnormal emission mode 17 times, accounting for ~13% of the 41.6 hours total observation time. Single pulse analyses indicate that mode changes took place simultaneously at 13 cm/3 cm within a few rotational periods. We detected occasional bright and narrow pulses whose peak flux densities were 10 times higher than that of the integrated profile in both bands. At 3 cm, about 0.66% and 0.27% of single pulses were bright in the normal mode and abnormal mode respectively, but at 13 cm the occurrence rate was only about 0.007%. We divided the pulsar radiation window into three components (C1, C2 and C3) corresponding to the main peaks of the integrated profile. The bright pulses preferentially occurred at pulse phases corresponding to the peaks of C2 and C3. Fluctuation spectra showed that C2 had excess red noise in the normal mode, but broad quasi-periodic features with central frequencies around 0.12 cycles/period in the abnormal mode. At 3 cm, C3 had a stronger quasi-periodic modulation centered around 0.06 cycles/period in the abnormal mode. Although there were some asymmetries in the two-dimensional fluctuation spectra, we found no clear evidence for systematic subpulse drifting. Consistent with previous low-frequency observations, we found a very low nulling probability for B0329+54 with upper limits of 0.13% and 1.68% at 13 cm/3 cm respectively.

We have investigated the mode-changing properties of PSR B0329+54 using 31 epochs of simultaneous 13 cm/3 cm single-pulse observations obtained with Shanghai Tian Ma 65 m telescope. The pulsar was found in the abnormal emission mode 17 times, accounting for ~13% of the 41.6 hours total observation time. Single pulse analyses indicate that mode changes took place simultaneously at 13 cm/3 cm within a few rotational periods. We detected occasional bright and narrow pulses whose peak flux densities were 10 times higher than that of the integrated profile in both bands. At 3 cm, about 0.66% and 0.27% of single pulses were bright in the normal mode and abnormal mode respectively, but at 13 cm the occurrence rate was only about 0.007%. We divided the pulsar radiation window into three components (C1, C2 and C3) corresponding to the main peaks of the integrated profile. The bright pulses preferentially occurred at pulse phases corresponding to the peaks of C2 and C3. Fluctuation spectra showed that C2 had excess red noise in the normal mode, but broad quasi-periodic features with central frequencies around 0.12 cycles/period in the abnormal mode. At 3 cm, C3 had a stronger quasi-periodic modulation centered around 0.06 cycles/period in the abnormal mode. Although there were some asymmetries in the two-dimensional fluctuation spectra, we found no clear evidence for systematic subpulse drifting. Consistent with previous low-frequency observations, we found a very low nulling probability for B0329+54 with upper limits of 0.13% and 1.68% at 13 cm/3 cm respectively.

We report the detection of widespread CH$_2$OHCHO and HOCH$_2$CH$_2$OH emission in Galactic center giant molecular cloud Sagittarius B2 using the Shanghai Tianma 65m Radio Telescope. Our observations show for the first time that the spatial distribution of these two important prebiotic molecules extends over 15 arc-minutes, corresponding to a linear size of approximately 36 pc. These two molecules are not just distributed in or near the hot cores. The abundance of these two molecules seems to decrease from the cold outer region to the central region associated with star-formation activity. Results present here suggest that these two molecules are likely to form through a low temperature process. Recent theoretical and experimental studies demonstrated that prebiotic molecules can be efficiently formed in icy grain mantles through several pathways. However, these complex ice features cannot be directly observed, and most constraints on the ice compositions come from millimeter observations of desorbed ice chemistry products. These results, combined with laboratory studies, strongly support the existence of abundant prebiotic molecules in ices.

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.

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

We report Shanghai Tian Ma Radio Telescope detections of several long carbon-chain molecules at C and Ku band, including HC3N, HC5N, HC7N, HC9N, C3S, C6H and C8H toward the starless cloud Serpens South 1a. We detected some transitions (HC9N J=13-12 F=12-11 and F=14-13, H13CCCN J=2-1 F=1-0 and F=1-1, HC13CCN J=2-1 F=2-2, F=1-0 and F=1-1, HCC13CN J=2-1 F=1-0 and F=1-1) and resolved some hyperfine components (HC5N J=6-5 F=5-4, H13CCCN J=2-1 F=2-1) for the first time in the interstellar medium. The column densities of these carbon-chain molecules in a range of 10^{12}-10^{13} cm^{-2} are comparable to two carbon-chain molecule rich sources, TMC-1 and Lupus-1A. The abundance ratios are 1.00:(1.11\pm0.15):(1.47\pm0.18) for [H13CCCN]:[HC13CCN]:[HCC13CN]. This result implies that the 13C isotope is also concentrated in the carbon atom adjacent to the nitrogen atom in HC3N in Serpens south 1a, which is similar to TMC-1. The [HC3N]/[H13CCCN] ratio of 78\pm9, the [HC3N]/[HC13CCN] ratio of 70\pm8, and the [HC3N]/[HCC13CN] ratio of 53\pm4 are also comparable to those in TMC-1. In any case, Serpens South 1a proves a testing ground for understanding carbon-chain chemistry.

In this paper, we report radio observations of the Galactic Center magnetar PSR J1745-2900 at six epochs between June and October, 2014. These observations were carried out using the new Shanghai Tian Ma Radio Telescope at a frequency of 8.6 GHz. Both the flux density and integrated profile of PSR J1745-2900 show dramatic changes from epoch to epoch showing that the pulsar was in its "erratic" phase. On MJD 56836, the flux density of this magnetar was about 8.7 mJy, which was ten times large than that reported at the time of discovery, enabling a single-pulse analysis. The emission is dominated by narrow "spiky" pulses which follow a log-normal distribution in peak flux density. From 1913 pulses, we detected 53 pulses whose peak flux density is ten times greater than that of the integrated profile. They are concentrated in pulse phase at the peaks of the integrated profile. The pulse widths at the 50% level of these bright pulses was between 0.2 to 0.9 deg, much narrower than that of integrated profile (~12 deg). The observed pulse widths may be limited by interstellar scattering. No clear correlation was found between the widths and peak flux density of these pulses and no evidence was found for subpulse drifting. Relatively strong spiky pulses are also detected in the other five epochs of observation, showing the same properties as that detected in MJD 56836. These strong spiky pulses cannot be classified as "giant" pulses but are more closely related to normal pulse emission.

Winking (oscillating) filaments have been observed for many years. However, observations of successive winking filaments in one event have not been reported yet. In this paper, we present the observations of a chain of winking filaments and a subsequent jet that are observed right after the X2.1 are in AR11283. The event also produced an Extreme-ultraviolet (EUV) wave that has two components: upward dome-like wave (850 km/s) and lateral surface wave (554 km/s) which was very weak (or invisible) in imaging observations. By analyzing the temporal and spatial relationships between the oscillating filaments and the EUV waves, we propose that all the winking filaments and the jet were triggered by the weak (or invisible) lateral surface EUV wave. The oscillation of the filaments last for two or three cycles, and their periods, Doppler velocity amplitudes, and damping times are 11 - 22 minutes, 6 - 14 km/s, and 25 - 60 minutes, respectively. We further estimate the radial component magnetic field and the maximum kinetic energy of the filaments, and they are 5 - 10 Gauss and ~ 10^19 J, respectively. The estimated maximum kinetic energy is comparable to the minimum energy of ordinary EUV waves, suggesting that EUV waves can efficiently launch filament oscillations on their path. Based on our analysis results, we conclude that the EUV wave is a good agent for triggering and connecting successive but separated solar activities in the solar atmosphere, and it is also important for producing solar sympathetic eruptions.

We present an observational study of a quasi-periodic fast propagating (QFP) magnetosonic wave on 2012, April 23. The multiple wave trains were observed along an active region open loop system which has a divergence geometry. The wave trains were first observed in 171 A observations at a distance of 150 Mm from the footpoint of the guiding loop system and with a speed of 689 km/s, then they appeared in 193 A observations after their interaction with a perpendicular underlaying loop system on the path, in the meantime, the wave speed decelerated to 343 km/s quickly within a short timescale. The sudden deceleration of the wave trains and their appearance in 193 A observations caused by the interaction are interpreted through geometric effect and the density increase of the guiding loop system, respectively. On the other hand, with Wavelet and Fourier analysis methods we find that the wave trains has a common period of 80 s with the associated flare. In addition, a few low frequencies are also identified in the QFP wave. We propose that the generation of the period of 80 s was caused by the periodic releasing of energy busts through some nonlinear processes in magnetic reconnection or the so-called oscillatory reconnection mechanism, while the low frequencies detected in the QFP wave were possibly the manifestations of the leakage of pressure-driven oscillations from the photosphere or chromosphere, which could be an important source for driving QFP waves in the low corona. Our observational results also indicate that the properties of the guiding magnetic structure such as the distributions of magnetic field and density as well as geometry are crucial for modulating the propagation behaviors of QFP waves. Therefore, QFP waves could be used for remote diagnostics of the local physical properties of the solar corona.

We present observations of the diffraction, refraction, and reflection of a global extreme-ultraviolet (EUV) wave propagating in the solar corona. These intriguing phenomena are observed when the wave interacts with two remote active regions, and they together exhibit the wave property of this EUV wave. When the wave approached AR11465, it became weaker and finally disappeared in the active region, but a few minutes latter a new wavefront appeared behind the active region, and it was not concentric with the incoming wave. In addition, a reflected wave was also observed simultaneously on the wave incoming side. When the wave approached AR11459, it transmitted through the active region directly and without reflection. The formation of the new wavefront and the transmission could be explained with diffraction and refraction effects, respectively. We propose that the different behaviors observed during the interactions may caused by different speed gradients at the boundaries of the two active regions. For the origin of the EUV wave, we find that it formed ahead of a group of expanding loops a few minutes after the start of the loops' expansion, which represents the initiation of the associated coronal mass ejection (CME). Based on these results, we conclude that the EUV wave should be a nonlinear magnetosonic wave or shock driven by the associated CME, which propagated faster than the ambient fast-mode speed and gradually slowed down to an ordinary linear wave. Our observations support the hybrid model that includes both fast wave and slow non-wave components.