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.

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.

The origin of the 6.67 hr period X-ray source, 1E161348-5055, in the young supernova remnant RCW 103 is puzzling. We propose that it may be the descendant of a Thorne-Zytkow Object (TZO). A TZO may at its formation have a rapidly spinning neutron star as a core, and a slowly rotating envelope. We found that the core could be braked quickly to an extremely long spin period by the coupling between its magnetic field and the envelope, and that the envelope could be disrupted by some powerful bursts or exhausted via stellar wind. If the envelope is disrupted after the core has spun down, the core will become an extremely long-period compact object, with a slow proper motion speed, surrounded by a supernova-remnant-like shell. These features all agree with the observations of 1E161348-5055. TZOs are expected to have produced extraordinary high abundances of lithium and rapid proton process elements that would remain in the remnants and could be used to test this scenario.

PSR B1821$-$24 is a solitary millisecond pulsar (MSP) which radiates multi-wavelength pulsed photons. It has complex radio, X-ray and $\gamma$-ray pulse profiles with distinct peak phase-separations that challenge the traditional caustic emission models. Using the single-pole annular gap model with suitable magnetic inclination angle ($\alpha=40^\circ$) and viewing angle ($\zeta=75^\circ$), we managed to reproduce its pulse profiles of three wavebands. It is found that the middle radio peak is originated from the core gap region at high altitudes, and the other two radio peaks are originated from the annular gap region at relatively low altitudes. Two peaks of both X-ray and $\gamma$-ray wavebands are fundamentally originated from annular gap region, while the $\gamma$-ray emission generated from the core gap region contributes somewhat to the first $\gamma$-ray peak. Precisely reproducing the multi-wavelength pulse profiles of PSR B1821$-$24 enables us to understand emission regions of distinct wavebands and justify pulsar emission models.

The very small braking index of PSR J1734-3333, $n=0.9\pm0.2$, challenges the current theories of braking mechanisms in pulsars. We present a possible interpretation that this pulsar is surrounded by a fall-back disk and braked by it. A modified braking torque is proposed based on the competition between the magnetic energy density of a pulsar and the kinetic energy density of a fall-back disk. With this torque, a self-similar disk can fit all the observed parameters of PSR J1734-3333 with natural initial parameters. In this regime, the star will evolve to the region having anomalous X-ray pulsars and soft gamma repeaters in the $P-\dot{P}$ diagram in about 20000 years and stay there for a very long time. The mass of the disk around PSR J1734-3333 in our model is about $10M_{\oplus}$, similar to the observed mass of the disk around AXP 4U 0142+61.

(adapted)Considering recent observations challenging the traditional magnetar model, we explore the wind braking of magnetars. There is evidence for strong multipole magnetic fields in active magnetars, but the dipole field inferred from spin down measurements may be strongly biased by a particle wind. Recent challenging observations of magnetars may be explained naturally in the wind braking scenario: (1) The supernova energies of magnetars are of normal value; (2) The non-detection in Fermi observations of magnetars; (3) The problem posed by the low-magnetic field soft gamma-ray repeaters; (4) The relation between magnetars and high magnetic field pulsars; (5) A decreasing period derivative during magnetar outbursts. Transient magnetars may still be magnetic dipole braking. This may explain why low luminosity magnetars are more likely to have radio emissions. In the wind braking scenario, magnetars are neutron stars with strong multipole field. For some sources, a strong dipole field may be no longer needed. A magnetism-powered pulsar wind nebula and a braking index smaller than three are the two predictions of the wind braking model.

Pulsed gamma-ray emission from millisecond pulsars (MSPs) has been detected by the sensitive Fermi, which sheds light on studies of the emission region and mechanism. In particular, the specific patterns of radio and gamma-ray emission from PSR J0101-6422 challenge the popular pulsar models, e.g. outer gap and two-pole caustic (TPC) models. Using the three dimension (3D) annular gap model, we have jointly simulated radio and gamma-ray light curves for three representative MSPs (PSR J0034-0534, PSR J0101-6422 and PSR J0437-4715) with distinct radio phase lags and present the best simulated results for these MSPs, particularly for PSR J0101-6422 with complex radio and gamma-ray pulse profiles and for PSR J0437-4715 with a radio interpulse. It is found that both the gamma-ray and radio emission originate from the annular gap region located in only one magnetic pole, and the radio emission region is not primarily lower than the gamma-ray region in most cases. In addition, the annular gap model with a small magnetic inclination angle instead of "orthogonal rotator" can account for MSPs' radio interpulse with a large phase separation from the main pulse. The annular gap model is a self-consistent model not only for young pulsars but also MSPs, and multi-wavelength light curves can be fundamentally explained by this model.

The concept of a "magnetar" was proposed mainly because of two factors. First, the X-ray luminosity of Anomalous X-ray Pulsars (AXPs) and Soft Gamma-Ray Repeaters (SGRs) is larger than the rotational energy loss rate, and second, the magnetic field strength calculated from "normal method" is super strong. It is proposed that the radiation energy of magnetar comes from its magnetic fields. Here it is argued that the magnetic field strength calculated through the normal method is incorrect when X-ray luminosity is larger than rotational energy loss rate, because the wind braking is not taken into account. Besides, the "anti-magnetar" and some other X-ray and radio observations are difficult to understand with a magnetar model. Instead of the magnetar, we propose a "quarctar", which is a crusted quark star in an accretion disk, to explain the observations. In this model, the persistent X-ray emission, burst luminosity, spectrum of AXPs and SGRs can be understood naturally. The radio-emitting AXPs, which are challenging the magnetar, can also be explained by the quarctar model.

Timing observations from the Parkes 64-m radio telescope for 165 pulsars between 1990 and 2011 have been searched for period glitches. A total of 107 glitches were identified in 36 pulsars, where 61 have previously been reported and 46 are new discoveries. Glitch parameters were measured by fitting the timing residual data. Observed relative glitch sizes \Delta\nu_g/\nu range between 10^-10 and 10^-5, where \nu = 1/P is the pulse frequency. We confirm that the distribution of \Delta\nu_g/\nu is bimodal with peaks at approximately 10^-9 and 10^-6. Glitches are mostly observed in pulsars with characteristic ages between 10^3 and 10^5 years, with large glitches mostly occurring in the younger pulsars. Exponential post-glitch recoveries were observed for 27 large glitches in 18 pulsars. The fraction Q of the glitch that recovers exponentially also has a bimodal distribution. Large glitches generally have low Q, typically a few per cent, but large Q values are observed in both large and small glitches. Observed time constants for exponential recoveries ranged between 10 and 300 days with some tendency for longer timescales in older pulsars. Shorter timescale recoveries may exist but were not revealed by our data which typically have observation intervals of 2 - 4 weeks. For most of the 36 pulsars with observed glitches, there is a persistent linear increase in \dot\nu in the inter-glitch interval. Where an exponential recovery is also observed, the effects of this are superimposed on the linear increase in \dot\nu. In some cases, the slope of the linear recovery changes at the time of a glitch. The \ddot\nu values characterising the linear changes in \dot\nu are almost always positive and, after subtracting the magnetospheric component of the braking, are approximately proportional to the ratio of |\dot\nu| and the inter-glitch interval, as predicted by vortex-creep models.

The multi-wavelength pulsed emission from young pulsars and millisecond pulsars can be well modeled with the single-pole 3-dimension annular gap and core gap model. To distinguish our single magnetic pole model from two-pole models (e.g. outer gap model and two-pole caustic model), the convincing values of the magnetic inclination angle and the viewing angle will play a key role.

To investigate the missing compact star of Supernova 1987A, we analyzed both the cooling and the heating processes of a possible compact star based on the upper limit of observational X-ray luminosity. From the cooling process we found that a solid quark-cluster star, which has a stiffer equation of state than that of conventional liquid quark star, has a heat capacity much smaller than a neutron star. It can cool down quickly, which can naturally explain the non-detection of a point source (neutron star or quark star) in X-ray band. On the other hand, we consider the heating process from magnetospheric activity and possible accretion, and obtain some constraints to the parameters of a possible pulsar. We conclude that a solid quark-cluster star can be fine with the observational limit in a large and acceptable parameter space. A pulsar with a short period and a strong magnetic field (or with a long period and a weak field) would has luminosity higher than the luminosity limit if the optical depth is not large enough to hide the compact star. The constraints of the pulsar parameters can be tested if the central compact object in 1987A is discovered by advanced facilities in the future.

In the framework of the three-dimensional (3D) annular gap model with reasonable parameters (the magnetic inclination angle \alpha = 45 deg and the view angle \zeta = 63 deg), we first use the latest hight energy data to self-consistently calculate radio, X-ray, gamma-ray and TeV (MAGIC and VERITAS) light curves, phase-averaged spectrum and phase-resolved spectra for the Crab pulsar. It is found that the acceleration electric field and potential in the annular gap and core gap are huge enough in the several tens of neutron star radii. The pulsed emission of radio, X-ray, gamma-ray and TeV are mainly generated from the emission of primary particles or secondary particles with different emission mechanisms in the nearly similar region of the annular gap located in the only one magnetic pole, which leads to the nearly "phase-aligned" multi-wavelength light curves. The emission of peak 1 (P1) and peak 2 (P2) is originated from the annular gap region near the null charge surface, while the emission of bridge is mainly originated from the core gap region. The phase-averaged spectrum and phase-resolved spectra of the Crab pulsar from soft X-ray to TeV band are produced by four components: synchrotron radiation from CR-induced and ICS-induced pairs dominates the X-ray band to soft gamma-ray band (100 eV to 10 MeV); curvature radiation and synchrotron radiation from the primary particles mainly contribute to gamma-ray band (10 MeV to \sim 20 GeV); ICS from the pairs significantly contributes to the TeV gamma-ray band (\sim 20 GeV to 400 GeV). The multi-wavelength pulsed emission from the Crab pulsar has been well modeled with the annular gap and core gap model. To distinguish our single magnetic pole model from two-pole models, the convincing values of the magnetic inclination angle and the viewing angle will play a key role.

The relationship between pulsar-like compact stars and gravitational waves is briefly reviewed. Due to regular spins, pulsars could be useful tools for us to detect ~nano-Hz low-frequency gravitational waves by pulsar-timing array technique; besides, they would also be ~kilo-Hz high-frequency gravitational wave radiators because of their compactness. The wave strain of an isolate pulsar depends on the equation state of cold matter at supra-nuclear densities. Therefore, a real detection of gravitational wave should be very meaningful in gravity physics, micro-theory of elementary strong interaction, and astronomy.

The inverse Compton scattering (ICS) model can explain various pulse profile shapes and diversity of pulse profile evolution based on the mechanism that the radio emission is generated through inverse Compton scattering between secondary relativistic particles and radio waves from polar gap avalanches. In this paper, we study the parameter space of ICS model for 15 pulsars, which share the common pulse profile evolution phenomena that the pulse profiles are narrower at higher observing frequencies. Two key parameters, the initial Lorentz factor and the energy loss factor of secondary particles are constrained using the least square fitting method, where we fit the theoretical curve of the "beam-frequency mapping" of the ICS model to the observed pulse widths at multiple frequencies. The uncertainty of the inclination and viewing angles are taken into account in the fitting process. It is found that the initial Lorentz factor is larger than 4000, and the energy loss factor is between 20 and 560. The Lorentz factor is consistent with the prediction of the inner vacuum gap model. Such high energy loss factors suggest significant energy loss for secondary particles at altitudes of a few tens to hundreds of kilometers.

The Vela pulsar represents a distinct group of {\gamma}-ray pulsars. Fermi {\gamma}-ray observations reveal that it has two sharp peaks (P1 and P2) in the light curve with a phase separation of 0.42 and a third peak (P3) in the bridge. The location and intensity of P3 are energy-dependent. We use the 3D magnetospheric model for the annular gap and core gap to simulate the {\gamma}-ray light curves, phase-averaged and phase-resolved spectra. We found that the acceleration electric field along a field line in the annular gap region decreases with heights. The emission at high energy GeV band is originated from the curvature radiation of accelerated primary particles, while the synchrotron radiation from secondary particles have some contributions to low energy {\gamma}-ray band (0.1 - 0.3 GeV). The {\gamma}-ray light curve peaks P1 and P2 are generated in the annular gap region near the altitude of null charge surface, whereas P3 and the bridge emission is generated in the core gap region. The intensity and location of P3 at different energy bands depend on the emission altitudes. The radio emission from the Vela pulsar should be generated in a high-altitude narrow regions of the annular gap, which leads to a radio phase lag of ~ 0.13 prior to the first {\gamma}-ray peak.

Though pulsars spin regularly, the differences between the observed and predicted ToA (time of arrival), known as "timing noise", can still reach a few milliseconds or more. We try to understand the noise in this paper. As proposed by Xu & Qiao in 2001, both dipole radiation and particle emission would result in pulsar braking. Accordingly, possible fluctuation of particle current flow is suggested here to contribute significant ToA variation of pulsars. We find that the particle emission fluctuation could lead to timing noise which can't be eliminated in timing process, and that a longer period fluctuation would arouse a stronger noise. The simulated timing noise profile and amplitude are in accord with the observed timing behaviors on the timescale of years.

Determining radiation location observationally plays a very important role in testing the pulsar radiation models. One-photon pair production in the strong magnetic field, $\gamma-e^{+}e^{1}$, is one of the important physical processes in pulsar radiation mechanisms. Photons near pulsar surface with sufficient energy will be absorbed in the magnetosphere and the absorption optical depth for these GeV $\gamma$-ray photons is usually large. In this paper, we include the aberrational, rotational and general relativistic effects and calculate the $\gamma$-B optical depth for $\gamma$-ray photons. Then we use the derived optical depth to determine the radiation altitude lower bounds for photons with given energies. As a case study, we calculate the lower bounds of radiation altitudes of Crab pulsar for photons with energy from 5 GeV to 1 TeV.

Anomalous X-ray Pulsars and Soft Gamma-Ray Repeaters have been generally recognized as neutron stars with super strong magnetic fields, namely "magnetars". The "magnetars" manifest that the luminosity in X-ray band are larger than the rotational energy loss rate, i.e. $L_{X}>\dot {E}_{\rm rot}$, and then the radiation energy is coming from the energy of magnetic field. Here it is argued that magnetars may not really exist. Some X-ray and radio observational results are contradicted with the magnetar model. (1) The X-ray luminosity of PSR J1852+0040 is much larger than the rotational energy loss rate ($L_{X}/\dot {E}_{\rm rot}\simeq 18)$, but the magnetic field is just $3.1\times 10^{11}$ G. Does this X-ray radiation energy come from the magnetic field? (2) In contrast to the above, the magnetic fields of radio pulsars J1847-0130 and PSR J1718-3718 are higher than that of AXP 1E 2259+586, why is the radiation energy of those two radio pulsars still coming from rotational energy? Furthermore, the magnetic field of the newly discovered SGR 0418+5729 with the lowest magnetic field is 3.0 e13 G, lower than the critical magnetic field $B_{\rm C}=4.414$ e13 G) (Esposito et al. 2010). (3) Some "magnetars" also emit normal transient radio pulses, what is the essential difference between radio pulsars and the "magnetars"? The observational fact arguments will be presented at first, then we discuss in what situation the conventional method to obtain magnetic field could not be correct.

Pulsed high energy radiation from pulsars is not yet completely understood. In this paper, we use the 3D self-consistent annular gap model to study light curves for both young and millisecond pulsars observed by the Fermi Gamma-ray Space Telescope. The annular gap can generate high energy emission for short-period pulsars. The annular gap regions are so large that they have enough electric potential drop to accelerate charged particles to produce gamma-ray photons. For young pulsars, the emission region is from the neutron star surface to about half of the light cylinder radius, and the peak emissivity is in the vicinity of the null charge surface. The emission region for the millisecond pulsars is located much lower than that of the young pulsars. The higher energy gamma-ray emission comes from higher altitudes in the magnetosphere. We present the simulated light curves for three young pulsars (the Crab, the Vela, the Geminga) and three millisecond pulsars (PSR J0030+0451, PSR J0218+4232, PSR J0437-3715) using the annular gap model. Our simulations can reproduce the main properties of observed light curves.

Pulsars have been recognized as normal neutron stars, but sometimes argued as quark stars. {\it Sub-millisecond pulsars, if detected, would play an essential and important role in distinguishing quark stars from neutron stars.} We focus on the formation of such sub-millisecond pulsars in this paper. A new approach to form a sub-millisecond pulsar (quark star) via accretion induced collapse (AIC) of a white dwarf is investigated here. Under this AIC process, we found that: (1) almost all the newborn quark stars could have an initial spin period of $\sim 0.1$ ms; (2) the nascent quark stars (even with a low mass) have sufficiently high spin-down luminosity and satisfy the conditions for pair production and sparking process to be as sub-millisecond radio pulsars; (3) in most cases, the timescales of newborn quark stars in the phase of spin period $< 1$ (or $<0.5$) ms can be long enough to be detected. As a comparison, an accretion spin-up process (for both neutron and quark stars) is also investigated. It is found that, quark stars formed through AIC process can have shorter periods ($\leq$ 0.5 ms); while the periods of neutron stars formed in accretion spin-up process must be longer than 0.5ms. Thus if a pulsar with a period less than 0.5 ms can be identified in the future, it should be a quark star.

Pulsars have been recognized as normal neutron stars or quark stars. Sub-millisecond pulsars, if detected, would play an essential and important role in distinguishing quark stars from neutron stars. A key question is how sub-millisecond pulsars could form. Both sub-Keplerian (for neutron and quark stars) and super-Keplerian cases (only for quark stars, which are bound additionally by strong interaction) have been discussed in this paper in order to investigate possible ways of forming sub-millisecond pulsars. In the sub-Keplerian case, the equilibrium periods of both neutron and quark stars could be as low as ~0.4 ms when they are spun up through accretion in binary systems. In the super-Keplerian case, pulsars could very likely have an initial period of ~0.1 ms if quark stars with different mass could be formed from accretion-induced collapse (AIC) of white dwarfs. The timescale for a sub-millisecond pulsar to keep its period <1 ms is restricted due to gravitational wave radiation. We found that the timescales of neutron stars within sub-millisecond periods are approximately ~10 yr, but the timescales estimated (>10 Myr) for low-mass quark stars could be long enough for us to detect.

Many models for the pulsar radio and $\gamma$-ray emissions have been developed. The tests for these models using observational data are very important. Tests for the pulsar radio emission models using frequency-altitude relation are presented in this paper. In the radio band, the mean pulse profiles evolve with observing frequencies. There are various styles of pulsar profile - frequency evolutions (which we call as "beam evolution" figure), e.g. some pulsars show that mean pulse profiles are wider and core emission is higher at higher frequencies than that at lower frequencies, but some other pulsars show completely the contrary results. We show that all these "beam evolution" figures can be understood by the Inverse Compton Scattering(ICS) model (see Qiao at al.2001 also). An important observing test is that, for a certain observing frequency different emission components are radiated from the different heights. For the $\gamma$-ray pulsars, the geometrical method (Wang et al. 2006) can be used to diagnose the radiation location for the $\gamma$-ray radiation. As an example, Wang et al. (2006) constrain the $\gamma$-ray radiation location of PSR B1055-52 to be the place near the null charge surface. Here we show that Wang's result matches the proposed radiation locations by the annular gap model as well as the outer gap models.

The consistence of the distributions of pulsar's kick velocities from the model of GRB and from the pulsar observations is tested based on the supernova-GRB ($\gamma$-ray burst) association and under the assumption that the GRB asymmetric explosions produce pulsars. The deduced distribution of kick velocity from the model of GRB and the observed kick distribution of radio pulsars are checked by K-S test. These two distributions are found to come from a same parent population. This result may indicate that GRBs could be really related to supernova, and that the asymmetry of GRB associated with supernova would cause the kick of pulsars.

If the binding energy of the pulsar's surface is not so high (the case of a neutron star), both the negative and positive charges will flow out freely from the surface of the star. The annular free flow model for $\gamma$-ray emission of pulsars is suggested in this paper. It is emphasized that: (1). Two kinds of acceleration regions (annular and core) need to be taken into account. The annular acceleration region is defined by the magnetic field lines that cross the null charge surface within the light cylinder. (2). If the potential drop in the annular region of a pulsar is high enough (normally the cases of young pulsars), charges in both the annular and the core regions could be accelerated and produce primary gamma-rays. Secondary pairs are generated in both regions and stream outwards to power the broadband radiations. (3). The potential drop in the annular region grows more rapidly than that in the core region. The annular acceleration process is a key point to produce wide emission beams as observed. (4). The advantages of both the polar cap and outer gap models are retained in this model. The geometric properties of the $\gamma$-ray emission from the annular flow is analogous to that presented in a previous work by Qiao et al., which match the observations well. (5). Since charges with different signs leave the pulsar through the annular and the core regions, respectively, the current closure problem can be partially solved.

The main idea and procedure are presented for a new synthetical method on constraining 3-D structure of emission regions of pulsars. With this method the emission regions can be synthetically constrained by fitting multi-wavelength observations features, e.g. pulse widths, phase offsets between radio pulse profiles and high energy light curves, and radio polarization properties. The main technique is based on numerically calculating the emission directions along any open field line via including aberration and retardation effects under a certain magnetic field configuration, e.g. vacuum dipole field. It can be easily extended by involving some further effects or magnetic field line configuration. The role of observed linear polarization properties in the method is discussed. The future application to known X/gamma-ray and radio pulsars will be helpful to test and improve current theoretical models on pulsar emission.

The pulse profile of pulsar gives geometric information about pulsar's radiation model. After investigating the pulse profiles of PSR B1642-03 and PSR B0950+08, we calculate the ratios of beam width and the emission height between different frequencies. We find that the ratios are almost constants as inclination angle $\alpha$ changes from $0^{\circ}$ to $90^{\circ}$. The ratios can be used to test pulsar's radiation model. In particular it can well constrain the parameter space of the inverse Compton scattering (ICS) model. These constrained parameters indicate some physical implication for the ICS model.

The large-scale magnetic field of our Galaxy can be probed in three dimensions using Faraday rotation of pulsar signals. We report on the determination of 223 rotation measures from polarization observations of relatively distant southern pulsars made using the Parkes radio telescope. Combined with previously published observations these data give clear evidence for large-scale counterclockwise fields (viewed from the north Galactic pole) in the spiral arms interior to the Sun and weaker evidence for a counterclockwise field in the Perseus arm. However, in interarm regions, including the Solar neighbourhood, we present evidence that suggests that large-scale fields are clockwise. We propose that the large-scale Galactic magnetic field has a bisymmetric structure with reversals on the boundaries of the spiral arms. Streaming motions associated with spiral density waves can directly generate such a structure from an initial inwardly directed radial field. Large-scale fields increase toward the Galactic Center, with a mean value of about 2~$\mu$G in the Solar neighbourhood and 4~$\mu$G at a Galactocentric radius of 3 kpc.

The aligned pulsars whose rotation axes and magnetic dipole axes are parallel should be positively charged. The total charge of pulsars is calculated after considering the electromagnetic field in and out the star under a specific condition. The statistical relation between the pulsar's rotation energy loss rate (or the period derivative) and the period may hint that the millisecond radio pulsars with small periods could be low-mass bare strange stars.

The global structure of the current flows in pulsar magnetospheres is investigated, with rough calculations of the elements in the magnetospheric circuit. It is emphasized that the potential of critical field lines is the same as that of interstellar medium, and that the pulsars whose rotation axes and magnetic dipole axes are parallel should be positively charged, in order to close the pulsar's current flows. The statistical relation between the radio luminosity and pulsar's electric charge (or the spindown power) may hint that the millisecond pulsars could be low-mass bare strange stars.

A new drifting pulsar, PSR J0815+09, was discovered in the Arecibo drift-scan searches. An intriguing feature of this source is that within the four pulse components in the integrated pulse profile, the sub-pulse drifting direction in the two leading components is opposite from that in the two trailing components. In view that the leading theoretical model for pulsar sub-pulse drifting can only interpret one-direction sub-pulse drifting, the observed "bi-drifting" phenomenon from PSR J0815+09 poses a great challenge to the pulsar theory. The inner annular gap (IAG), a new type of inner particle accelerator, was recently proposed to explain both $\gamma$-ray and radio emission from pulsars. Here we show that the coexistence of the IAG and the conventional inner core gap (ICG) offers a natural interpretation to the bi-drifting phenomenon. In particular, the peculiar drifting behavior in PSR J0815+09 can be reproduced within the inverse Compton scattering (ICS) model for pulsar radio emission.

The inner annular gap (IAG), a new type of inner gap whose magnetic field lines intersect the null charge surface (NCS), is proposed to explain $\gamma$-ray and radio emission from pulsars. The IAG can be an important source for high-energy particles. The particles can radiate between the NCS and the IAG. Some observational characteristics in both $\gamma$-ray and radio bands, such as the $\gamma$-ray emission beams of Crab-like, Vela-like and Geminga-like, can be reproduced by numerical method. It is predicted that the view angle $\zeta$ should be larger than the inclination angle ($\zeta>\alpha$), otherwise the $\gamma$-ray radiation will have little possibility to be observed. Whether the inner annular gap (or cap) is sparking (or free flow) depends on the surface binding energy of the pulsar. In stead of neutron star models, the scenario of the IAG is favorable for bare strange star models, because bare strange stars can easily satisfy the requisite condition to form an IAG for both pulsars ($\vec{\Omega}\cdot \vec{B}<0$) and anti-pulsars ($\vec\Omega \cdot \vec{B}>0$).

Pulsars are among the most highly polarized sources in the universe. The NVSS has catalogued 2 million radio sources with linear polarization measurements, from which we have selected 253 sources, with polarization percentage greater than 25%, as targets for pulsar searches. We believe that such a sample is not biased by selection effects against ultra-short spin or orbit periods. Using the Parkes 64m telescope, we conducted searches with sample intervals of 0.05 ms and 0.08 ms, sensitive to submillisecond pulsars. Unfortunately we did not find any new pulsars.

The long-term, highly periodic and correlated variations in both the pulse shape and the rate of slow-down of two isolated pulsars (PSRs) PSR B1828-11 and PSR B1642-03 were discovered recently. This phenomenon may provide evidence for "free precession" as suggested in the literature. Some authors presented various kinds of models to explain this phenomenon within the framework of free precession. Here we present an accretion disk model for this precession phenomenon instead. Under reasonable parameters, the observed phenomenon can be explained by an isolated pulsar with a fossil disk. This may link radio pulsars and anomalous X-ray pulsars (AXPs) and present an indirect evidence for the existence of the fossil disk in nature.

The radiation of pulsars have been observed for many years. A few pulsars are discovered to have both radio and gamma-ray emission. Many models on pulsar radiation have been developed, but so far we are still lacking an elaborate model which can explain the emission from radio to gamma-rays in detail. In this paper we present a joint model for radio and gamma-ray emission, in which both the dominate emission mechanisms are inverse Compton scattering. The pulse profiles at radio and gamma-ray bands are reproduced for the Crab-like, Vela-like and Geminga-like pulsars, in good agreement with observations.

Although pulsars can radiate electromagnetic wave from radio to gamma ray bands, we still have no a united model to understand the multi-band emission. In this paper the effort for a joint model is presented. The inverse Compton scattering (ICS) and a second acceleration process near the null charge surface are involved to account for the radio and the gama-ray emission, respectively. Various kind of pulse profiles and other observational properties can be reproduced.

Pulsars provide unique probes for the {\it large-scale} interstellar magnetic field in the Galactic disk. The Parkes multibeam pulsar survey has discovered many distant pulsars which enables us for the first time to explore the magnetic field in most of the nearby half of the Galactic disk. The fields are found to be {\it coherent in direction} over a linear scale of $\sim 10$ kpc between the Carina-Sagittarius and Crux-Scutum arms from $l\sim45\degr$ to $l\sim305\degr$. The coherent spiral structures and field direction reversals, including the newly determined counterclockwise field near the Norma arm, are consistent with bi-symmetric spiral model for the disk field. However, the antisymmetric rotation measure sky from the Galactic halo and the dipole field in the Galactic center suggest that the A0 dynamo is operating there.

In a paper by Sanwal et al. (2002), it is supposed to be very difficult to interpret the absorption features in term of cyclotron lines. However, we would like to address here that the possibility of the absorption being cyclotron resonance can not be ruled out.We propose that the isolate neutron star, 1E 1207.4-5209 in the center of supernova remnant PKS 1209-51/52, has a debris disk and is in a propeller phase, with an accretion rate 6x10^{-11} M_solar/year.

The fastest millisecond pulsar PSR B1937+21 presents an interpulse separated from the main pulse by nealy 180^o at radio frequencies. Recently, the ASCA observations (Takahashi et al. 2001) detected pulsed X-ray emission from this pulsar. Only a single narrow X-ray pulse is observed, which is coincident with the radio interpulse in phase. We investigate the possible origin of the pulsed X-rays from the polar cap (PC) accelerators or the outer gap (OG) accelerators in the frame of PC model and OG model, respectively, by assuming a dipolar magnetic field structure and the same radio emission pattern from its poles for the pulsar. The OG model can naturally explain the main observational facts. For the PC model, the coincidence between the X-ray pulse and the radio interpulse can not be reproduced in the assumed case. However when considering possible deviation from our assumption, PC model may still be valid for this pulsar in some cases.

According to two estimated relations between the initial period and the dynamo-generated magnetic dipole field of pulsars, we calculate the statistical distributions of pulsar initial periods. It is found that proto-pulsars are very likely to have rotation periods between 20 and 30 ms, and that most of the pulsars rotate initially at a period < 60 ms.

Pulsars provide unique probes of the large-scale interstellar magnetic field in the Galactic disk. Up to now, the limited Galactic distribution of the known pulsar population has restricted these investigations to within a few kiloparsec of the Sun. The Parkes multibeam pulsar survey has discovered many more-distant pulsars which enables us for the first time to explore the magnetic field in most of the nearby half of the Galactic disk. Here we report the detection of counterclockwise magnetic fields in the Norma spiral arm using pulsar rotation measures. The fields are coherent in direction over a linear scale of $\sim 5$ kpc along the arm and have a strength of $-4.4\pm0.9$ $\mu$G. The magnetic field between the Carina-Sagittarius and Crux-Scutum arms is confirmed to be coherent from $l\sim45\degr$ to $l\sim305\degr$ over a length of $\sim 10$ kpc. These results strengthen arguments for a bisymmetric spiral model for the field configuration in the Galactic disk.

Pulsar braking torques due to magnetodipole radiation and to unipolar generator are considered, which results in braking index being less than 3 and could be employed to test the emission models. Improved equations to obtain pulsar braking index and magnetic field are presented if we deem that the rotation energy loss rate equals to the sum of the dipole radiation energy loss rate and that of relativistic particles powered by unipolar generator. The magnetic field calculated by conventional way could be good enough but only modified by a factor of ~0.6 at most. Both inner and outer gaps may coexist in the magnetosphere of the Vela pulsar.

The shapes of pulse profiles, especially their variations with respect to observing frequencies, are very important to understand emission mechanisms of pulsars, while no previous attempt has been made in interpreting the complicated phenomenology. In this paper, we present theoretical simulations for the integrated pulse profiles and their frequency evolution within the framework of the inverse Compton scattering (ICS) model proposed by Qiao (1988) and Qiao & Lin (1998). Using the phase positions of the pulse components predicted by the ``beam-frequency figure'' of the ICS model, we present Gaussian fits to the multi-frequency pulse profiles for some pulsars. It is shown that the model can reproduce various types of the frequency evolution behaviors of pulse profiles observed.

Extremely relativistic particles play an important role in various radiative mechanisms of pulsar radio as well as high-energy emission. It is thus very essential to estimate the Lorentz factor of these particles in the radiation regions both observationally and theoretically. In this paper, according to the frequency dependence of component separation and considering some radiation processes, we obtain the Lorentz factor and its variation in the radiation region observationally for seven multi-frequency observed pulsars. It is found that the Lorentz factor decays substantially in the emission region, and the decay can not be well understood by all known energy loss processes. This result hints that there are some unknown processes which are responsible to loosing energy effectively.

It is of great important to study pulsar beam shapes if we are concerned with emission theories and pulsar birth rate. Both observations and/or the inverse Compton scattering model show that different emission components are emitted from different heights. The relative longitude phase shifts due to different heights of emission components and to the toroidal velocity of electron are considered in this paper. Several possible observational features arising from the phase shift effects are presented. The emission beams may not have circular cross sections although the emission regions may be symmetric with respect to the magnetic axes.

During supernova explosions, strange stars with almost bare quark surfaces may be formed. Under certain conditions, these stars could be rapidly spun down by the torque exerted by the fossil disks formed from the fall-back materials. They may also receive large kicks and hence, have large proper motion velocities. When these strange stars pass through the spherical ``Oort'' comet cloud formed during the pre-supernova era, they will capture some small-scale comet clouds and collide with some comet-like objects occasionally. These impacts can account for the repeating bursts as observed from the soft gamma repeaters (SGRs). According to this picture, it is expected that SGR 1900+14 will become active again during 2004-2005.

The Orthogonal Polarization Modes (OPM) have been reported observationally and accepted widely by pulsar researchers. However, no acceptable theory can show the origin of the OPM, which becomes a mystery in pulsar research field. Here a possible way to solve this mystery is presented. We ask a question: Does there exist any real so-called OPM in pulsar radiation? It is proposed in this paper that the `observed OPM' in individual pulses could be the results of depolarization of pulsar radiation and the observational uncertainties originated from polarimeter in observation. A possible way to check this idea is suggested. If the idea is verified, the pulsar research would be influenced significantly in theory and in observation.

The possibility and the implications of the idea, that pulsars are born as strange stars, are explored. Strange stars are very likely to have atmospheres with typical mass of $\sim 5\times 10^{-15}M_\odot$ but bare polar caps almost throughout their lifetimes, if they are produced during supernova explosions. A direct consequence of the bare polar cap is that the binding energies of both positively and negatively charged particles at the bare quark surface are nearly infinity, so that the vacuum polar gap sparking scenario as proposed by Ruderman & Sutherland should operate above the cap, regardless of the sense of the magnetic pole with respect to the rotational pole. Heat can not accumulate on the polar cap region due to the large thermal conductivity on the bare quark surface. We test this ``bare polar cap strange star'' (BPCSS) idea with the present broad band emission data of pulsars, and propose several possible criteria to distinguish BPCSSs from neutron stars.

Some characteristics of the inverse Compton scattering (ICS) model are reviewed. At least the following properties of radio pulsars can be reproduced in the model: core or central emission beam, one or two hollow emission cones, different emission heights of these components, diverse pulse profiles at various frequencies, linear and circular polarization features of core and cones.

We propose that the violent breakdown of the Ruderman-Sutherland type vacuum gap above the pulsar polar cap region may produce copious primary particles and make the charge density above the gap in excess of the Goldreich-Julian density. Such an excess residue of charge density may arise some re-acceleration regions above the gap due to the space-charge-limited flow mechanism. Such re-acceleration regions, if exist, may accelerate particles to carry away more spin-down energy, and may have some implications for pulsar $\gamma$-ray emission.

We carried out the Gaussian fitting to the profile of PSR B1237+25 and found that six components rather than five are necessary to make a good fit. In the central part, we found that the core emission is not filled pencil beam but is a small hollow cone. This implies that the impact angle could be $\beta<0.5^\circ$. The ``hollow core'' is in agreement with Inverse Compton Scattering model of radio pulsars.

Polarization profiles of several hundred pulsars have been published recently (eg. Gould & Lyne 1998, GL98; Weisberg et al. 1999). In this report, we summarize the characteristics of circular and high linear polarizations of pulsar profiles, based on all previously published data.

Qiao and his collaborators recently proposed an inverse Compton Scattering (ICS) model to explain radio emission of pulsars. In this paper, we investigate the polarization properties of pulsar emission in the model. First of all, using the lower frequency approximation, we derived the analytical amplitude of inverse Compton scattered wave of a single electron in strong magnetic field. We found that the out-going radio emission of a single relativistic electron scattering off the ``low frequency waves'' produced by gap-sparking should be linearly polarized and have no circular polarization at all. However, considering the coherency of the emission from a bunch of electrons, we found that the out-going radiation from the inner part of emission beam, i.e., that from the lower emission altitudes, prefers to have circular polarization. Computer simulations show that the polarization properties, such as the sense reversal of circular polarization near the pulse center, S-shape of position angle swing of the linear polarization, strong linear polarization in conal components, can be reproduced in the ICS model.

Many theoretical efforts were made to understand the core and conal emission identified from observations by Rankin (1983) and Lyne and Manchester (1988). One of them, named as inverse Compton scattering (ICS) model (Qiao & Lin 1998), has been proposed. It is found in the model that: there are central or `core' emission beam, and one or two hollow conical emission beams; the different emission components are emitted at different heights; owing to different radiation components emitted from different height, the observed emission beams can be shifted from each other due to retardation and aberration effects; the sizes of emission components change with frequencies. Recent developments of the model include: simulations of pulse profiles at different frequencies; studying the basic polarization properties of inverse Compton scattering in strong magnetic fields; computing the polarizations and spectrum of core and cones. A new classification system was also proposed. The main results calculated from the model are consistent with the observations.

Similar to the terrestrial collision accelerators of electron-antielectron, another kind of accelerator is above a positively or negatively charged pulsar polar cap. In the case of pulsars with magnetic axis parallel (anti-parallel) to rotational axis, relativistic antielectron (electron) with Lorentz factor about 10^6 hit the electrons in the polar caps. These scenarios are investigated both for pulsars being bare strange stars and for pulsars being neutron stars. Such a study may be valuable to differentiate neutron stars and bare strange stars observationally.

It is suggested that the inner energetic engine of Gamma ray burst (GRB) may be the result of the transition of normal hadron to quark-gluon plasma (QGP) in rapidly-rotating and spin-down newborn neutron star. When such a nascent neutron star slows down through dipole electromagnetic and quadruple gravitational radiation, the increasing center density may reach the QCD transition density, i.e., 5-10 nuclear density. Such kind of energy release from the phase transition would be responsible for GRB and its possible beaming effect. The relative dense gaseous environment of GRB location and the iron line observed in the X-ray afterglow support this idea. Some predictions in this model are given. keywords: Gamma-ray Bursts, Supernova, Quark-gluon plasma

Using the Thomas-Fermi model, we investigated the electric characteristics of a static non-magnetized strange star without crust in this paper. The exact solutions of electron number density and electric field above the quark surface are obtained. These results are useful if we are concerned about physical processes near the quark matter surfaces of strange stars.

We have obtained 63 rotation measures (RMs) from polarization observations of southern pulsars, of which 54 are new measurements and 3 are varied from previous values. The new pulsar RM data at high Galactic latitudes are mostly consistent with the antisymmetric RM distribution found previously. For the Galactic disc, evidence for a field reversal near the Perseus arm, and possibly another beyond it, is presented. Inside the Solar Circle, in addition to the two known field reversals in or near the Carina-Sagittartus arm and the Crux-Scutum arm, a further reversal in the Norma arm is tentatively identified. These reversals, together with the pitch angle derived from pulsar RM and stellar polarization distributions, are consistent with bisymmetric spiral (BSS) models for the large-scale magnetic field structure in the disc of our Galaxy. However, discrimination between models is complicated by the presence of smaller-scale irregularities in the magnetic field, as well as uncertainties in the theoretical modelling.

Based on definitions, two joint-criteria, namely, the optical-thin constraint and the energy budget constraint, are proposed to judge whether the emission nature of radio pulsars is incoherent or obligatory to be coherent. We find that the widely accepted criterion, $kT_B \le \epsilon$, is not a rational criterion to describe the optical-thin condition, even for the simplest case. The energy budget constraint could be released by introducing a certain efficient radiation mechanism (e.g. the inverse Compton scattering, QL98) with emission power of a single particle as high as a critical value $P_{sing,c} \sim 10^{-4}-10^{-3} {erg s}^{-1}$. This in principle poses the possibilities to interpret high luminosities of pulsars in terms of incoherent emission mechanisms, if the optical-thin constraint could be released by certain mechanism as well. Coherence may not be an essential condition to account for pulsar radio emission.

It is proposed that the `bare' strange matter stars might not be bare, and radio pulsars might be in fact `bare' strange stars. As strange matter stars being intensely magnetized rotate, the induced unipolar electric fields would be large enough to construct magnetospheres. This situation is very similar to that discussed by many authors for rotating neutron stars. Also, the strange stars with accretion crusts in binaries could act as X-ray pulsars or X-ray bursters. There are some advantages if radio pulsars are `bare' strange stars.

A model is suggested for the bursting x-ray pulsar, GRO J1744-28, where the compact object of the pulsar is proposed to be a strange matter star with crusts. It is suggested that two envelope crusts shield the two polar caps of an accreting strange stars owing to strong polar-cap magnetic fields. GRO J1744-28 might not have a unified crust, but two polar crusts. Bursts are the result of phase transition of polar crusts if one accreting crust is heavier than that the Coulomb force can support. It is shown that there should be phase-lags when bursting crust is expanded about one hundred meters. Many calculated properties of this model, such as the phase-lags, the bursting luminosity and the characteristic photon energy, are in agreement with observations.

We present a systematic study of the circular polarization in pulsar integrated profiles, based on published polarization data. For core components, we find no significant correlation between the sense-change of circular polarization and the sense of linear position angle variation. Circular polarization is not restricted to core components and, in some cases, reversals of circular polarization sense are observed across the conal emission. In conal double profiles, the sense of circular polarization is found to be correlated with the sense of position-angle variation. Pulsars with a high degree of linear polarization often have one hand of circular polarization across the whole profile. For most pulsars, the sign of circular polarization is the same at 50-cm and 20-cm wavelength, and the degree of polarization is similar, albeit with a wide scatter. However, at least two cases of frequency-dependent sign reversals are known. This diverse behaviour may require more than one mechanism to generate circular polarization.

The importance of the photon-photon pair production process ($\gamma+ \gamma^{\prime}\to e^{+}+e^{-}$) to form pair production cascades in pulsar polar caps is investigated within the framework of the Ruderman-Sutherland vacuum gap model. It is found that this process is unimportant if the polar caps are not hot enough, but will play a non-negligible role in the pair formation cascades when the polar cap temperatures are in excess of the critical temperatures, $T_{cri}$, which are around $4\times 10^6K$ when $P=0.1$s and will slowly increase with increasing periods. Compared with the $\gamma-B$ process, it is found that the two-photon annihilation process may ignite a central spark near the magnetic pole, where $\gamma-B$ sparks can not be formed due to the local weak curvatures. This central spark is large if the gap is dominated by the ``resonant ICS mode''. The possible connection of these central sparks with the observed pulsar ``core'' emission components is discussed.

Although the unified formula for gamma-ray absorption process involving both the magnetic field and a perpendicular electric field derived by Daugherty & Lerche (1975) is correct, we argued in this paper that their conclusion that the induced electric fields are important in the pair formation process in the pulsar magnetospheres is wrong and misleading. The key point is that usually the direction of a gamma photon at the emission point observed in the laboratory frame should be (v/c, 0, [1-(v/c)^2]^{1/2}) rather than (0, 0, 1), where v is the co-rotating velocity. This emission direction is just the one which results in zero attenuation coefficient of the gamma photon. Calculation shows that after the photon has moved a distance, its direction lead to the result that the induced electric field is also of minor importance. Thus only gamma-B process is the important mechanism for the pair production in the pulsar magnetospheres. The implications of the modification by ejecting the induced electric field are also discussed.

It is suggested in this paper that the `bare' strange star might be not bare, and there could be a magnetosphere around it. As a strange star might be an intensely magnetized rotator, the induced unipolar electric field would be large enough to construct a magnetosphere around the strange matter core. This kind of magnetosphere is very similar to that of the rotating magnetized neutron stars discussed by many authors. A magnetosphere will be established very soon through pair production by gamma-B or two photon processes after a strange star was born in a supernova explosion. It is emphasized that the fact that the strange star surface can not supply charged particles does not stop the formation of a space charge separated magnetosphere around the bare strange star. An accretion crust is quite difficult to come into being around an isolated strange star. Therefore the observed radio signals of an rotation-powered pulsar may come from a bare strange stars rather than a neutron stars or a strange star with an accretion crust. The idea, that the radio pulsars are the strange stars without crusts, is supported by some observations. The fact, that the iron emission lines have been observed in many X-ray pulsars but never been reported in X-ray emission of radio pulsars, is difficult to be understood if the radio pulsars are the neutron stars or strange star with crust where the surface binding energy of iron ions is too low to avoid a ion free-flow from the surface.

Core and conal emission beams of radio pulsars have been identified observationally (Rankin 1983,1993; Lyne and Manchester 1988), In the inner gap model (Ruderman and Sutherland 1975, hereafter RS75), the gap continually breaks down (sparking) by forming electron-positron pairs on a time scale of a few microseconds (RS75). This makes a large amplitude low frequency wave to be produced and which would be scattered by relativistic particles moving out from the gap. Under this assumption, Qiao (1988a,1992) presented an Inverse Compton Scattering (ICS) model for both core and conal emissions. This paper presented a development of the model. Retardation and aberration effects due to different radiation components emitted from different heights are considered. The luminosity of pulsar radio emission via the ICS process is discussed. Coherent emission by bunches of particles is adopted and which is adequate to explain pulsar radiation. The theoretical results are in agreement with observations very well.