We present here the analysis of giant micropulses from the Vela pulsar. A total of 4187 giant micropulses with peak flux density $>$2.5 Jy were detected during almost 4 hours of observations carried out with the Yunnan 40-m radio telescope at 6800 MHz. Nine of the giant micropulses arrived approximately 3 to 4 ms earlier than the peak of average pulse profile, longer than that at lower frequencies. The remaining giant micropulses were clustered into three distributions which correspond to three main emission regions, including four occurring on the trailing edge of averaged profile.

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.

By analysing the data acquired from the Parkes 64-m radio telescope at 1369 MHz, we report on the phase-stationary non-drift amplitude modulation observed in PSR J1048-5832. The high-sensitivity observations revealed that the central and trailing components of the pulse profile of this pulsar switch between a strong mode and a weak mode periodically. However, the leading component remains unchanged. Polarization properties of the strong and weak modes are investigated. Considering the similarity to mode changing, we argue that the periodic amplitude modulation in PSR J1048$-$5832 is periodic mode changing. The fluctuation spectral analysis showed that the modulation period is very short (~2.1 s or 17 P1), where P1 is the rotation period of the pulsar. We find that this periodic amplitude modulation is hard to explain by existing models that account for the periodic phenomena in pulsars like subpulse drifting.

We derive novel constraints on millicharged dark matter and ultralight axion-like particles using pulsar timing and fast radio burst observations. Millicharged dark matter affects the dispersion measure of the time of arrival of radio pulses in a way analogous to free electrons. Light pseudo-scalar dark matter, on the other hand, causes the polarization angle of radio signals to oscillate. We show that current and future data can set strong constraints in both cases. For dark matter particles of charge $\epsilon e$, these constraints are ${\epsilon}/{m_{\rm milli}} \lesssim 10^{-8}{\rm eV}^{-1}$, for masses $m_{\rm milli}\gtrsim 10^{-6}\,$eV. For axion-like particles, the analysis of signals from pulsars yields constraints in the axial coupling of the order of $g/m_a\lesssim 10^{-13} {\rm GeV}^{-1}/(10^{-22}{\rm eV})$. Both bounds scale as $(\rho/\rho_{\rm dm})^{1/2}$ if the energy density $\rho$ of the components is a fraction of the total dark matter energy density $\rho_{\rm dm}$. We do a detailed study of both effects using data from two samples of pulsars in the galaxy and in globular clusters, as well as data from FRB 121102 and PSR J0437$-$4715. We show that in both cases actual pulsar data constrain a new region of the parameter space for these models, and will improve with future pulsar-timing observations.

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

PSR J1825$-$0935 (PSR B1822$-$09) switches between radio-quiet (Q-mode) and radio-bright (B-mode) modes. The Q-mode is known to have a periodic fluctuation that modulates both the interpulse and the main pulse with the same period. Earlier investigators argued that the periodic Q-mode modulation is associated with drifting subpulses. We report on single-pulse observations of PSR J1825$-$0935 that were made using the Parkes 64-m radio telescope with a central frequency of 1369 MHz. The high-sensitivity observations revealed that the periodic Q-mode modulation is in fact a periodic longitude-stationary intensity modulation occurring in the interpulse and the main pulse. The fluctuation spectral analysis showed that the modulation period is about $43 P_1$, where $P_1$ is the rotation period of the pulsar. Furthermore, we confirm that the modulation patterns in the interpulse and the main pulse are phase-locked. Specifically, the intensities of the interpulse and the immediately following main pulse are more highly correlated than for the main pulse and interpulse at any other lag. Polarization properties of the strong and weak Q-mode states are different, even for the trailing part of the main pulse which does not show the periodic intensity modulation.

A primary aim of the Nuclear Spectroscopic Telescope Array (NuSTAR) mission is to find and characterize heavily obscured Active Galactic Nuclei (AGNs). Based on mid-infrared photometry from the Wide-Field Infrared Survey Explorer (WISE) and optical photometry from the Sloan Digital Sky Surveys, we have selected a large population of luminous obscured AGN (i.e., obscured quasars). Here we report NuSTAR observations of four WISE-selected heavily obscured quasars for which we have optical spectroscopy from the Southern African Large Telescope and W. M. Keck Observatory. Optical diagnostics confirm that all four targets are AGNs. With NuSTAR hard X-ray observations, three of the four objects are undetected, while the fourth has a marginal detection. We confirm that these objects have observed hard X-ray (10-40 keV) luminosities at or below ~ 10^43 erg s^-1. We compare X-ray and IR luminosities to obtain estimates of the hydrogen column densities (N_H) based on the suppression of the hard X-ray emission. We estimate N_H of these quasars to be at or larger than 10^25 cm^-2, confirming that WISE and optical selection can identify very heavily obscured quasars that may be missed in X-ray surveys, and do not contribute significantly to the cosmic X-ray background. From the optical Balmer decrements, we found that our three extreme obscured targets lie in highly reddened host environments. This galactic extinction is not adequate to explain the more obscured AGN, but it may imply a different scale of obscuration in the galaxy.

We presented timing results of PSR B2035+36 using ~9-yr observations with the Nanshan 25-m radio telescope. PSR B2035+36 was reported to exhibit significant changes in pulse profile correlated with spin-down state variations. We found that the pulsar underwent a glitch with a jump in the frequency of $\Delta{\nu}\sim12.4(5)\, \rm nHz$ around MJD 52950. Unusually, the spin-down rate increased persistently over $800$ days after the glitch, and the average spin-down rate of post-glitch was about $9.6\%$ larger than pre-glitch. Accompanied with the glitch activity, pulse profile became narrower. The pulsar began to switch between two emission modes after the glitch, with pulse width ($W_{\ 50 \rm mean}$) of $8.5(7)^{\circ}$ and $3.7(3)^{\circ}$, respectively. Besides that, the relatively narrow pulse profile gradually became dominant. All of the observations indicate that there should be connection between magnetospheric behavior and glitch activity. We discuss one possibility of magnetosphere fluctuation triggered by glitch event.

We report on single-pulse observations of the Galactic Center magnetar PSR J1745$-$2900 that were made using the Parkes 64-m radio telescope with a central frequency of 3.1 GHz at five observing epochs between 2013 July and August. The shape of the integrated pulse profiles was relatively stable across the five observations, indicating that the pulsar was in a stable state between MJDs 56475 and 56514. This extends the known stable state of this pulsar to 6.8 months. Short term pulse shape variations were also detected. It is shown that this pulsar switches between two emission modes frequently and that the typical duration of each mode is about ten minutes. No giant pulses or subpulse drifting were observed. Apparent nulls in the pulse emission were detected on MJD 56500. Although there are many differences between the radio emission of magnetars and normal radio pulsars, they also share some properties. The detection of mode changing and pulse nulling in PSR J1745$-$2900 suggests that the basic radio emission process for magnetars and normal pulsars is the same.

We present the results of an optical spectroscopic survey of 46 heavily obscured quasar candidates. Objects are selected using their mid-infrared (mid-IR) colours and magnitudes from the Wide-Field Infrared Survey Explorer (WISE) and their optical magnitudes from the Sloan Digital Sky Survey (SDSS). Candidate Active Galactic Nuclei (AGNs) are selected to have mid-IR colours indicative of quasar activity and lie in a region of mid-IR colour space outside previously published X-ray based selection regions. We obtain optical spectra for our sample using the Robert Stobie Spectrograph on the Southern African Large Telescope. Thirty objects (65%) have identifiable emission lines, allowing for the determination of spectroscopic redshifts. Other than one object at $z\sim2.6$, candidates have moderate redshifts ranging from $z=0.1$ to $0.8$ with a median of 0.3. Twenty-one (70%) of our objects with identified redshift (46% of the whole sample) are identified as AGNs through common optical diagnostics. We model the spectral energy distributions of our sample and found that all require a strong AGN component, with an average intrinsic AGN fraction at 8$\,\mu$m of 0.91. Additionally, the fits require large extinction coefficients with an average $E(B-V)_\textrm{AGN} = 17.8$ (average $A(V)_\textrm{AGN} = 53.4$). By focusing on the area outside traditional mid-IR photometric cuts, we are able to capture and characterise a population of deeply buried quasars that were previously unattainable through X-ray surveys alone.

Pulse arrival time measurements over fourteen years with the Nanshan 25-m and Parkes 64-m radio telescopes have been used to determine the average profile and timing properties for PSR J1757$-$2421 (B1754$-$24). Analysis of the radio profile data shows a large variation of spectral index across the profile and an unusual increase in component separation with increasing frequency. The timing observations show that PSR B1754$-$24 underwent a large glitch with a fractional increase in spin frequency of $\Delta\nu_{\rm g}/\nu\sim 7.8\times10^{-6}$ in May 2011. At the time of the glitch there was a large permanent jump in the rotational frequency accompanied by two exponential recovery terms with timescales of 15(6) and 97(15) days, respectively.

This paper reports on polarimetric radiation properties based on the switching modes of normal PSR B2020+28 by analysing the data acquired from the Nanshan 25-m radio telescope at 1556 MHz. With nearly 8 hours quasi-continuous observation, the data presented some striking and updated phenomena. The change of relative intensity between the leading and trailing components is the predominant feature of mode switching. The intensity ratio between the leading and trailing components are measured for the individual profiles averaged over 30 seconds. It is found that there is an excess of high ratios over the normal distribution, which indicates that two modes exist in the pulsar. The distribution of abnormal mode has a narrower width indicating that the abnormal mode is more stable than the normal mode. A total of 76 mode switching events are detected in our data. It spends 89\% in the normal mode and 11\% in the abnormal mode. The intrinsic distributions of mode timescales are constrained with power-law distributions. The significant difference in the index of the duration distribution between normal and abnormal modes possibly indicates that the timescale for the abnormal mode to get stable is shorter than that for the normal mode. The frequent switching between both modes may indicate that the oscillations between different magnetospheric states are rapid.

Aims. To make a detailed study of the nulling and subpulse drifting in PSR J1727$-$2739 for investigation of its radiation properties. Methods. The observations were carried out on 20 March, 2004 using the Parkes 64-m radio telescope, with a central frequency of 1518 MHz. A total of 5568 single pulses were analysed. Results. This pulsar shows well defined nulls with lengths lasting from 6 to 281 pulses and separated by burst phases ranging from 2 to 133 pulses. We estimate a nulling fraction of around 68\%. No emission in the average pulse profile integrated over all null pulses is detected with significance above $3\sigma$. Most transitions from nulls to bursts are within a few pulses, whereas the transitions from bursts to nulls exhibit two patterns of decay: decrease gradually or rapidly. In the burst phase, we find that there are two distinct subpulse drift modes with vertical spacing between the drift bands of $9.7 \pm 1.6$ and $5.2 \pm 0.9$ pulse periods, while sometimes there is a third mode with no subpulse drifting. Some mode transitions occur within a single burst, while others are separated by nulls. Different modes have different average pulse profiles. Possible physical mechanisms are discussed.

Polish doughnuts (PDs) are geometrically thick disks that rotate with super-Keplerian velocities in their innermost parts, and whose long and narrow funnels along rotation axes collimate the emerging radiation into beams. In this paper we construct an extremal family of PDs that maximize both geometrical thickness and radiative efficiency. We then derive upper limits for these quantities and subsequently for the related ability to collimate radiation. PDs with such extreme properties may explain the observed properties of the ultraluminous X-ray sources without the need for the black hole masses to exceed ~ 10 solar masses. However, we show that strong advective cooling, which is expected to be one of the dominant cooling mechanisms in accretion flows with super-Eddington accretion rates, tends to reduce the geometrical thickness and luminosity of PDs substantially. We also show that the beamed radiation emerging from the PD funnels corresponds to "isotropic" luminosities that linearly scale with the mass accretion rate, and do not obey the familiar and well-known logarithmic relation.

We investigate the structural properties of the underlying hosts of 34 blue compact dwarf (BCD) galaxies with deep near-infrared (NIR) photometry. The BCD sample is selected from the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey in the Great observatories origins Deep Survey North and South fields. We extract the surface brightness profile (SBP) in the optical F 435W and NIR F 160W bands. The SBPs of BCDs in the H band reach 26 mag arcsec^-2 at the 3\sigma level, which is so far the deepest NIR imaging of BCDs. Then we fit the SBPs with one- and two- component Sersic models. About half of the BCDs favour the two-component model which significantly improves the fit quality. The effective radii of the underlying hosts of BCDs in the B band are smaller than those of early-type dwarfs (dEs) and dwarf irregulars at a fixed luminosity. This discrepancy is similar to findings in many previous works. However, the difference in structural parameters between BCDs and other dwarf galaxies seems to be less significant in the H band. Furthermore, we find a remarkable agreement between the underlying hosts of BCDs and dEs. All dwarf galaxies seem to follow a similar luminosity-radius relationship which suggests a unified structural evolution for dwarf galaxies. We conclude that a possible evolution track from BCDs to dEs cannot be ruled out, with no significant change of structure needed in the evolutionary scenario.

We present high signal-to-noise ratio, multi-frequency polarization pulse profiles for 24 millisecond pulsars that are being observed as part of the Parkes Pulsar Timing Array (PPTA) project. The pulsars are observed in three bands, centred close to 730, 1400 and 3100 MHz, using a dual-band 10 cm/50 cm receiver and the central beam of the 20 cm multibeam receiver. Observations spanning approximately six years have been carefully calibrated and summed to produce high S/N profiles. This allows us to study the individual profile components and in particular how they evolve with frequency. We also identify previously undetected profile features. For many pulsars we show that pulsed emission extends across almost the entire pulse profile. The pulse component widths and component separations follow a complex evolution with frequency; in some cases these parameters increase and in other cases they decrease with increasing frequency. The evolution with frequency of the polarization properties of the profile is also non-trivial. We provide evidence that the pre- and post-cursors generally have higher fractional linear polarization than the main pulse. We have obtained the spectral index and rotation measure for each pulsar by fitting across all three observing bands. For the majority of pulsars, the spectra follow a single power-law and the position angles follow a $\lambda^2$ relation, as expected. However, clear deviations are seen for some pulsars. We also present phase-resolved measurements of the spectral index, fractional linear polarization and rotation measure. All these properties are shown to vary systematically over the pulse profile.

Intermittent pulsars are part-time radio pulsars. They have higher slow down rate in the on state (radio-loud) than in the off state (radio-quiet). This gives the evidence that particle wind may play an important role in pulsar spindown. The effect of particle acceleration is included in modeling the rotational energy loss rate of the neutron star. Applying the pulsar wind model to the three intermittent pulsars (PSR B1931+24, PSR J1841-0500, and PSR J1832+0029), their magnetic field and inclination angle are calculated simultaneously. The theoretical braking indices of intermittent pulsars are also given. In the pulsar wind model, the density of the particle wind can always be the Goldreich-Julian density. This may ensure that different on states of intermittent pulsars are stable. The duty cycle of particle wind can be determined from timing observations. It is consistent with the duty cycle of the on state. Inclination angle and braking index observations of intermittent pulsars may help to test different models of particle acceleration. At present, the inverse Compton scattering induced space charge limited flow with field saturation model can be ruled out.

We investigate the mass-metallicity relation and its dependence on galaxy physical properties with a sample of 703 Lyman-break analogues (LBAs) in local Universe, which have similar properties to high redshift star-forming galaxies. The sample is selected according to $\ha$ luminosity, $L(\ha)>10^{41.8}\,{\rm erg\,s^{-1}}$, and surface brightness, $I(\ha)>10^{40.5}\,{\rm erg\,s^{-1}\,kpc^{-2}}$, criteria. The mass-metallicity relation of LBAs harmoniously agrees with that of star-forming galaxies at $z \sim$ 1.4-1.7 in stellar mass range of $10^{8.5}M_{\odot}<M_{*}<10^{11}M_{\odot}$. The relation between stellar mass, metallicity and star formation rate of our sample is roughly consistent with the local fundamental metallicity relation. We find that the mass-metallicity relation shows a strong correlation with the 4000\AA\, break; galaxies with higher 4000\AA\, break typically have higher metallicity at a fixed mass, by 0.06 dex in average. This trend is independent of the methodology of metallicity. We also use the metallicity estimated by $T_{\rm e}$-method to confirm it. The scatter in mass-metallicity relation can be reduced from 0.091 to 0.077 dex by a three-dimensional relation between stellar mass, metallicity and 4000\AA\, break. The reduction of scatter in mass-metallicity relation suggests that the galaxy stellar age plays an important role as the second parameter in the mass-metallicity relation of LBAs.

Context. The supermassive black hole, Sagittarius (Sgr) A*, in the centre of our Galaxy has the largest angular size in the sky among all astrophysical black holes. Its shadow, assuming no rotation, spans ~ 50 microarcsec. Resolving such dimensions has long been out of reach for astronomical instruments until a new generation of interferometers being operational during this decade. Of particular interest is the Event Horizon Telescope (EHT) with resolution ~ 20 microarcsec in the millimeter-wavelength range 0.87 mm - 1.3 mm. Aims. We investigate the ability of the fully general relativistic Komissarov (2006) analytical magnetized torus model to account for observable constraints at Sgr A* in the centimeter and millimeter domains. The impact of the magnetic field geometry on the observables is also studied. Methods. We calculate ray-traced centimeter- and millimeter-wavelength synchrotron spectra and images of a magnetized accretion torus surrounding the central black hole in Sgr A*. We assume stationarity, axial symmetry, constant specific angular momentum and polytropic equation of state. A hybrid population of thermal and non-thermal electrons is considered. Results. We show that the torus model is capable of reproducing spectral constraints in the millimeter domain, and in particular in the observable domain of the EHT. However, the torus model is not yet able to fit the centimeter spectrum. 1.3 mm images at high inclinations are in agreement with observable constraints. Conclusions. The ability of the torus model to account for observations of Sgr A* in the millimeter domain is interesting in the perspective of the future EHT. Such an analytical model allows very fast computations. It will thus be a suitable test bed for investigating large domains of physical parameters, as well as non-black-hole compact object candidates and alternative theories of gravity.

A "pulsar timing array" (PTA), in which observations of a large sample of pulsars spread across the celestial sphere are combined, allows investigation of "global" phenomena such as a background of gravitational waves or instabilities in atomic timescales that produce correlated timing residuals in the pulsars of the array. The Parkes Pulsar Timing Array (PPTA) is an implementation of the PTA concept based on observations with the Parkes 64-m radio telescope. A sample of 20 millisecond pulsars is being observed at three radio-frequency bands, 50cm (~700 MHz), 20cm (~1400 MHz) and 10cm (~3100 MHz), with observations at intervals of 2 - 3 weeks. Regular observations commenced in early 2005. This paper describes the systems used for the PPTA observations and data processing, including calibration and timing analysis. The strategy behind the choice of pulsars, observing parameters and analysis methods is discussed. Results are presented for PPTA data in the three bands taken between 2005 March and 2011 March. For ten of the 20 pulsars, rms timing residuals are less than 1 microsec for the best band after fitting for pulse frequency and its first time derivative. Significant "red" timing noise is detected in about half of the sample. We discuss the implications of these results on future projects including the International Pulsar Timing Array (IPTA) and a PTA based on the Square Kilometre Array. We also present an "extended PPTA" data set that combines PPTA data with earlier Parkes timing data for these pulsars.

We report on variations in the mean position angle of the 20 millisecond pulsars being observed as part of the Parkes Pulsar Timing Array (PPTA) project. It is found that the observed variations are dominated by changes in the Faraday rotation occurring in the Earth's ionosphere. Two ionospheric models are used to correct for the ionospheric contribution and it is found that one based on the International Reference Ionosphere gave the best results. Little or no significant long-term variation in interstellar RM was found with limits typically about 0.1 rad m$^{-2}$ yr$^{-1}$ in absolute value. In a few cases, apparently significant RM variations over timescales of a few 100 days or more were seen. These are unlikely to be due to localised magnetised regions crossing the line of sight since the implied magnetic fields are too high. Most probably they are statistical fluctuations due to random spatial and temporal variations in the interstellar electron density and magnetic field along the line of sight.

Polarization profiles are presented for 20 millisecond pulsars that are being observed as part of the Parkes Pulsar Timing Array project. The observations used the Parkes multibeam receiver with a central frequency of 1369 MHz and the Parkes digital filterbank pulsar signal-processing system PDFB2. Because of the large total observing time, the summed polarization profiles have very high signal/noise ratios and show many previously undetected profile features. Thirteen of the 20 pulsars show emission over more than half of the pulse period. Polarization variations across the profiles are complex and the observed position angle variations are generally not in accord with the rotating-vector model for pulsar polarization. Never-the-less, the polarization properties are broadly similar to those of normal (non-millisecond) pulsars, suggesting that the basic radio emission mechanism is the same in both classes of pulsar. The results support the idea that radio emission from millisecond pulsars originates high in the pulsar magnetosphere, probably close to the emission regions for high-energy X-ray and gamma-ray emission. Rotation measures were obtained for all 20 pulsars, eight of which had no previously published measurements.