The "ATOMS," standing for {\it ALMA Three-millimeter Observations of Massive Star-forming regions}, survey has observed 146 active star forming regions with ALMA Band 3, aiming to systematically investigate the spatial distribution of various dense gas tracers in a large sample of Galactic massive clumps, to study the roles of stellar feedback in star formation, and to characterize filamentary structures inside massive clumps. In this work, the observations, data analysis, and example science of the "ATOMS" survey are presented, using a case study for the G9.62+0.19 complex. Toward this source, some transitions, commonly assumed to trace dense gas, including CS $J = 2-1$, HCO$^+$ $J = 1-0$ and HCN $J = 1-0$, are found to show extended gas emission in low density regions within the clump; less than 25\% of their emission is from dense cores. SO, CH$_3$OH, H$^{13}$CN and HC$_3$N show similar morphologies in their spatial distributions and reveal well the dense cores. Widespread narrow SiO emission is present (over $\sim$1 pc), which may be caused by slow shocks from large--scale colliding flows or H{\sc ii} regions. Stellar feedback from an expanding H{\sc ii} region has greatly reshaped the natal clump, significantly changed the spatial distribution of gas, and may also account for the sequential high-mass star formation in the G9.62+0.19 complex. The ATOMS survey data can be jointly analyzed with other survey data, e.g., "MALT90", "Orion B", "EMPIRE", "ALMA\_IMF", and "ALMAGAL", to deepen our understandings of "dense gas" star formation scaling relations and massive proto-cluster formation.

We report studies of the relationships between the total bolometric luminosity ($L_{\rm bol}$ or $L_{\rm TIR}$) and the molecular line luminosities of $J=1-0$ transitions of H$^{13}$CN, H$^{13}$CO$^+$, HCN, and HCO$^+$ with data obtained from ACA observations in the "ATOMS" survey of 146 active Galactic star forming regions. The correlations between $L_{\rm bol}$ and molecular line luminosities $L'_{\rm mol}$ of the four transitions all appear to be approximately linear. Line emission of isotopologues shows as large scatters in $L_{\rm bol}$-$L'_{\rm mol}$ relations as their main line emission. The log($L_{\rm bol}$/$L'_{\rm mol}$) for different molecular line tracers have similar distributions. The $L_{\rm bol}$-to-$L'_{\rm mol}$ ratios do not change with galactocentric distances ($R_{\rm GC}$) and clump masses ($M_{\rm clump}$). The molecular line luminosity ratios (HCN-to-HCO$^+$, H$^{13}$CN-to-H$^{13}$CO$^+$, HCN-to-H$^{13}$CN and HCO$^+$-to-H$^{13}$CO$^+$) all appear constant against $L_{\rm bol}$, dust temperature ($T_{\rm d}$), $M_{\rm clump}$ and $R_{\rm GC}$. Our studies suggest that both the main lines and isotopologue lines are good tracers of the total masses of dense gas in Galactic molecular clumps. The large optical depths of main lines do not affect the interpretation of the slopes in star formation relations. We find that the mean star formation efficiency (SFE) of massive Galactic clumps in the "ATOMS" survey is reasonably consistent with other measures of the SFE for dense gas, even those using very different tracers or examining very different spatial scales.

Fast radio bursts (FRBs) are short pulses observed in radio band from cosmological distances, some of which emit repeating bursts. The physical origins of these mysterious events have been subject to wide speculations and heated debates. One class of models invoke soft gamma-ray repeaters (SGRs), or magnetars, as the sources of FRBs. Magnetars are rotating neutron stars with extremely strong magnetic field and can sporadically emit bursts from X-ray ($\sim$keV) to soft gamma-ray ($\sim$sub-MeV) with duration from $10^{-2}$~s to $10^2$~s. However, even though some bright radio bursts have been observed from some magnetars, no FRB-like events had been detected to be associated with any magnetar burst, including one giant flare, and no radio burst has been associated with any X-ray event from any magnetar. Therefore, there is still no observational evidence for magnetar-FRB association up to today. Recently, a pair of FRB-like bursts (FRB~200428 hereafter) separated by 30 milliseconds (ms) were detected from the general direction of the Galactic magnetar SGR~J1935+2154. Here we report the detection of a non-thermal X-ray burst in the 1-250 keV energy band with the Insight-HXMT satellite, which we identify as emitted from SGR~J1935+2154. The burst showed two hard peaks with a separation of $\sim30$ ms, consistent with the separation between the two bursts in FRB~200428. The delay time between the double radio and X-ray peaks is $\sim8.57$ s, fully consistent with the dispersion delay of FRB~200428. We thus identify the non-thermal X-ray burst is associated with FRB~200428 whose high energy counterpart is the two hard peaks in X-ray. Our results suggest that the non-thermal X-ray burst and FRB~200428 share the same physical origin in an explosive event from SGR~J1935+2154.

The CsI detectors of the High Energy X-ray Telescope of the Hard X-ray Modulation Telescope (HXMT/CsI) can be used for gamma-ray all sky monitoring and searching for the electromagnetic counterpart of gravitational wave source. The instrumental responses are mainly obtained by Monte Carlo simulation with the Geant4 tool and the mass model of both the satellite and all the payloads, which is updated and tested with the Crab pulse emission in various incident directions. Both the Energy-Channel relationship and the energy resolution are calibrated in two working modes (Normal-Gain mode & Low-Gain Mode) with the different detection energy ranges. The simulative spectral analyses show that HXMT/CsI can constrain the spectral parameters much better in the high energy band than that in the low energy band. The joint spectral analyses are performed to ten bright GRBs observed simultaneously with HXMT/CsI and other instruments (Fermi/GBM, Swift/BAT, Konus-Wind), and the results show that the GRB flux given by HXMT/CsI is systematically higher by $7.0\pm8.8\%$ than those given by the other instruments. The HXMT/CsI-Fermi/GBM joint fittings also show that the high energy spectral parameter can be constrained much better as the HXMT/CsI data are used in the joint fittings.

We use the In data collected during the 2019 outburst from X-ray pulsar 4U 1901+03 to complement the orbital parameters reported by Fermi/GBM. Using the Insight-HXMT, we examine the correlation between the derivative of the intrinsic spin frequency and bolometric flux based on accretion torque models. It was found that the pulse profiles significantly evolve during the outburst. The existence of two types of the profile's pattern discovered in the Insight-HXMT data indicates that this source experienced transition between a super-critical and a sub-critical accretion regime during its 2019 outburst. Based on the evolution of the pulse profiles and the torque model, we derive the distance to 4U 1901+03 as 12.4+-0.2 kpc.

We study the evolution of the temporal properties of MAXI 1820+070 during the 2018 outburst in its hard state from MJD 58190 to 58289 with Insight-HXMT in a broad energy band 1-150 keV. We find different behaviors of the hardness ratio, the fractional rms and time lag before and after MJD 58257, suggesting a transition occurred at around this point. The observed time lags between the soft photons in the 1-5 keV band and the hard photons in higher energy bands, up to 150 keV, are frequency-dependent: the time lags in the low-frequency range, 2-10 mHz, are both soft and hard lags with a timescale of dozens of seconds but without a clear trend along the outburst; the time lags in the high-frequency range, 1-10 Hz, are only hard lags with a timescale of tens of milliseconds; first increase until around MJD 58257 and decrease after this date. The high-frequency time lags are significantly correlated to the photon index derived from the fit to the quasi-simultaneous NICER spectrum in the 1-10 keV band. This result is qualitatively consistent with a model in which the high-frequency time lags are produced by Comptonization in a jet.

Glitches correspond to sudden jumps of rotation frequency ($\nu$) and its derivative ($\dot{\nu}$) of pulsars, the origin of which remains not well understood yet, partly because the jump processes of most glitches are not well time-resolved. There are three large glitches of the Crab pulsar, detected in 1989, 1996 and 2017, which were found to have delayed spin-up processes before the normal recovery processes. Here we report two additional glitches of the Crab pulsar occurred in 2004 and 2011 for which we discovered delayed spin up processes, and present refined parameters of the largest glitch occurred in 2017. The initial rising time of the glitch is determined as $<0.48$ hour. We also carried out a statistical study of these five glitches with observed spin-up processes. The two glitches occurred in 2004 and 2011 have delayed spin-up time scales ($\tau_{1}$) of $1.7\pm0.8$\,days and $1.6\pm0.4$\,days, respectively. We find that the $\Delta{\nu}$ vs. $|\Delta{\dot\nu}|$ relation of these five glitches is similar to those with no detected delayed spin-up process, indicating that they are similar to the others in nature except that they have larger amplitudes. For these five glitches, the amplitudes of the delayed spin-up process ($|\Delta{\nu}_{\rm d1}|$) and recovery process ($\Delta{\nu}_{\rm d2}$), their time scales ($\tau_{1}$, $\tau_{2}$), and permanent changes in spin frequency ($\Delta{\nu}_{\rm p}$) and total frequency step ($\Delta{\nu}_{\rm g}$) have positive correlations. From these correlations, we suggest that the delayed spin-up processes are common for all glitches, but are too short and thus difficult to be detected for most glitches.

No robust detection of prompt electromagnetic counterparts to fast radio bursts (FRBs) has yet been obtained, in spite of several multi-wavelength searches carried out so far. Specifically, X/gamma-ray counterparts are predicted by some models. We planned on searching for prompt gamma-ray counterparts in the Insight-Hard X-ray Modulation Telescope (Insight-HXMT) data, taking advantage of the unique combination of large effective area in the keV-MeV energy range and of sub-ms time resolution. We selected 39 FRBs that were promptly visible from the High-Energy (HE) instrument aboard Insight-HXMT. After calculating the expected arrival times at the location of the spacecraft, we searched for a significant excess in both individual and cumulative time profiles over a wide range of time resolutions, from several seconds down to sub-ms scales. Using the dispersion measures in excess of the Galactic terms, we estimated the upper limits on the redshifts. No convincing signal was found and for each FRB we constrained the gamma-ray isotropic-equivalent luminosity and the released energy as a function of emission timescale. For the nearest FRB source, the periodic repeater FRB180916.J0158+65, we find $L_{\gamma,iso}<5.5\times 10^{47}$ erg/s over 1 s, whereas $L_{\gamma,iso}<10^{49}-10^{51}$ erg/s for the bulk of FRBs. The same values scale up by a factor of ~100 for a ms-long emission. Even on a timescale comparable with that of the radio pulse itself no keV-MeV emission is observed. A systematic association with either long or short GRBs is ruled out with high confidence, except for subluminous events, as is the case for core-collapse of massive stars (long) or binary neutron star mergers (short) viewed off axis. Only giant flares from extra-galactic magnetars at least ten times more energetic than Galactic siblings are ruled out for the nearest FRB.

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

Multi-peaked features appear on the integrated pulse profiles of fast radio burst observed below 2.5 GHz and the instantaneous spectrum of many bursts observed between 4 and 8 GHz. The mechanism of pulse or spectrum shaping has attracted little attention. Here we show that these interference-like pulse profiles are mostly the instantaneous spectra near the source regions of fast radio bursts. The corresponding instantaneous spectra are coincident to the spectrum from a single electron passing through a tapered undulator. The multi-peaked spectrum observed between 4 and 8 GHz can also be explained consistently by this type of spectrum. The spectrum is invisible unless the particles in the radiation beam are bunched. The bunching effect is probably due to the acceleration of particles in the plasma wakefield.

A new black hole X-ray binary (BHXRB) MAXI J1535-571 was discovered by MAXI during its outburst in 2017. Using observations taken by the first Chinese X-ray satellite, the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT), we perform a joint spectral analysis (2-150 keV) in both energy and time domains. The energy spectra provide the essential input for probing the intrinsic Quasi-Periodic Oscillation (QPO) fractional rms spectra (FRS). Our results show that during the intermediate state, the energy spectra are in general consistent with those reported by Swift/XRT and NuSTAR. However, the QPO FRS become harder and the FRS residuals may suggest the presence of either an additional power-law component in the energy spectrum or a turn-over in the intrinsic QPO FRS at high energies.

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

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

We present a detailed timing study of the brightest persistent X-ray source Sco X-1 using the data collected by the Hard X-ray Modulation Telescope ($Insight$-HXMT) from July 2017 to August 2018. A complete $Z$-track hardness-intensity diagram (HID) is obtained. The normal branch oscillations (NBOs) at $\sim$ 6 Hz in the lower part of the normal branch (NB) and the flare branch oscillations (FBOs) at $\sim$ 16 Hz in the beginning part of the flaring branch (FB) are found in observations with the Low Energy X-ray Telescope (LE) and the Medium Energy X-ray Telescope (ME) of $Insight$-HXMT, while the horizontal branch oscillations (HBOs) at $\sim$ 40 Hz and the kilohertz quasi-periodic oscillations (kHz QPOs) at $\sim$ 800 Hz are found simultaneously up to 60 keV for the first time on the horizontal branch (HB) by the High Energy X-ray Telescope (HE) and ME. We find that for all types of the observed QPOs, the centroid frequencies are independent of energy, while the root mean square (rms) increases with energy; the centroid frequencies of both the HBOs and kHz QPOs increase along the $Z$-track from the top to the bottom of the HB; and the NBOs show soft phase lags increasing with energy. A continuous QPO transition from the FB to NB in $\sim$ 200 s are also detected. Our results indicate that the non-thermal emission is the origin of all types of QPOs, the innermost region of the accretion disk is non-thermal in nature, and the corona is nonhomogeneous geometrically.

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

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

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

Fast radio bursts are radio transients observed mainly around 1.5 GHz. Their peak frequency decreases at a rate of 100 ~ 500 MHz/s and some of them have a broader pulse with an exponentially decaying tail. Common assumptions for fast radio bursts include a dispersion effect resulting in the peak frequency drifting and a scattering effect resulting in pulse broadening. These assumptions attribute the abnormally large dispersion measure and scattering measure to the environmental medium of the host galaxy. Here we show that the radiation of fast radio bursts can be explained as an undulator radiation and the large dispersion measure can be due to a motion effect mainly from the rotation of the source which is probably variable stars. In our scenario, the pulse broadening is near-field effects and the pulse itself represents a Fresnel diffraction pattern sweeping the observer. Our work is the first analysis of properties of fast radio bursts in the context of a special mechanism of the radiation instead of a special propagation environment of the radiation.

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

A model for the evolution of low-luminosity radio galaxies is presented. In the model, the lobes inflated by low-power jets are assumed to expand in near pressure-balance against the external medium. Both cases of constant external pressure and decreasing external pressure are considered. Evolution of an individual source is described by the power-size track. The source appears as its lobe is inflated and radio luminosity increases to above the detection level; the source then moves along the track and eventually disappears as its luminosity drops below the detection limit. The power-size tracks are calculated including the combined energy losses due to synchrotron radiation, adiabatic expansion, and inverse Compton scattering. It is shown that in general, the constant-pressure model predicts an excess number of luminous, small-size sources while underpredicting large-size sources in the power-size diagram. The predicted spectra are steep for most sources, which is inconsistent with observations. By comparison, the pressure-limiting model fits observations better. In this model, low-luminosity sources undergo substantial expansion losses in the initial phase and as a result, it predicts fewer luminous, small-size sources. The resultant spectra are flat for most sources except for the oldest ones, which seems consistent with observations. The power-size tracks, in contrast to that of high-luminosity radio galaxies, are characterized by a slow increase in luminosity for most of the source's life, followed by a rapid decline when the synchrotron or inverse Compton scattering losses set in.

Cosmic-ray streaming instabilities at supernova shocks are discussed in the quasilinear diffusion formalism which takes into account the feedback effect of wave growth on the cosmic ray streaming motion. In particular, the nonresonant instability that leads to magnetic field amplification in the short wavelength regime is considered. The linear growth rate is calculated using kinetic theory for a streaming distribution. We show that the nonresonant instability is actually driven by a compensating current in the background plasma. The nonresonant instability can develop into a nonlinear regime generating turbulence. The saturation of the amplified magnetic fields due to particle diffusion in the turbulence is derived analytically. It is shown that the evolution of parallel and perpendicular cosmic-ray pressures is predominantly determined by nonresonant diffusion. However, the saturation is determined by resonant diffusion which tends to reduce the streaming motion through pitch angle scattering. The saturated level can exceed the mean background magnetic field.

We consider the motion of a charge in a large amplitude electrostatic wave with a triangular wave form relevant to an oscillating model of a pulsar magnetosphere. The (one-dimensional) orbit of a particle in such a wave is found exactly in terms of Weierstrass functions. The result is used to discuss linear acceleration emission (at both low and high frequencies) in an oscillating model for pulsars. An explicit expression for the emissivity is derived in an accompanying paper (Melrose & Luo 2009), and used here to derive an expression for the absorption coefficient at low frequencies. We show that absorption can be negative, corresponding to maser emission. For the large amplitude required to trigger pair creation in an oscillating model, the rate of the maser growth is too small to be effective. Effective growth requires smaller amplitudes, such that the maximum Lorentz factor gained by acceleration in the wave is $\lapprox10$.

The theory of linear acceleration emission is developed for a large amplitude electrostatic wave in which all particles become highly relativistic in much less than a wave period. An Airy integral approximation is shown to apply near the phases where the electric field passes through zero and the Lorentz factors of all particles have their maxima. The emissivity is derived for an individual particle and is integrated over frequency and solid angle to find the power radiated per particle. The result is different from that implied by the generalized Larmor formula which, we argue, is not valid in this case. We also discuss a mathematical inconsistency that arises when one evaluates the power spectrum by integrating the emissivity over solid angle. The correct power spectrum increases as the 4/3rd power of the frequency at low frequencies, and falls off exponentially above a characteristic frequency. We discuss application of linear acceleration emission to the emission of high frequency photons in an oscillating model for pulsars. We conclude that it cannot account for gamma-ray emission, but can play a role in secondary pair creation.

An analytical model for oscillating pair creation above the pulsar polar cap is presented in which the parallel electric field is treated as a large amplitude, superluminal, electrostatic wave. An exact formalism for such wave is derived in one-dimension and applied to both the low-density regime in which the pair plasma density is much lower than the corotating charge density and the high-density regime in which the pair plasma density is much higher than the corotating charge density. In the low-density regime, which is relevant during the phase leading to a pair cascade, a parallel electric field develops resulting in rapid acceleration of particles. The rapid acceleration leads to bursts of pair production and the system switches to the oscillatory phase, corresponding to the high density regime, in which pairs oscillate with net drift motion in the direction of wave propagation. Oscillating pairs lead to a current that oscillates with large amplitude about the Goldreich-Julian current. The drift motion can be highly relativistic if the phase speed of large amplitude waves is moderately higher than the speed of light. Thus, the model predicts a relativistic outflow of pairs, a feature that is required for avoiding overheating of the pulsar polar cap and is also needed for the pulsar wind.

Magnetic dissipation through decay of Alfven waves in the magnetar magnetosphere is discussed. Transport of magnetic fields in the star leads to dissipation of the magnetic energy through either direct internal heating or transferring of the energy in waves that decay in the magnetar magnetosphere. In the latter case, the Alfven waves are excited by crust dislocations or elastic waves underneath the star's surface. It is suggested that these Alfven waves can decay into ion sound waves which can be effectively damped leading to strong plasma heating. Hot plasmas expand producing transient X-rays.

It is proposed that radiation belts similar to the ones in the planetary magnetosphere can exist for a pulsar with a relatively long period and a strong magnetic field. In the belts located in the closed field line region near the light cylinder relativistic pairs are trapped and maintained at a density substantially higher than the local Goldreich-Julian corotation density. The trapped plasma can be supplied and replenished by either direct injection of relativistic pairs from acceleration of externally-supplied particles in a dormant outer gap or in situ ionization of the accreted neutral material in the trapping region. The radiation belts can be disrupted by waves that are excited in the region as the result of plasma instabilities or emitted from the surface due to starquakes or stellar oscillations. The disruption can cause an intermittent particle precipitation toward the star producing radio bursts. It is suggested that such bursts may be seen as rotating radio transients (RRATs).

The effect of photon-beam-induced turbulence on propagation of radio emission in a pulsar magnetosphere is discussed. Beamed radio emission with a high brightness temperature can generate low-frequency plasma waves in the pulsar magnetosphere and these waves scatter the radio beam. We consider this effect on propagation of radio emission both in the open field line region and in the closed field line region. The former is applicable to most cases of pulsar radio emission where the propagation is confined to the polar region; it is shown that the induced process is not effective for radio emission of moderately high brightness temperature but can have a severe effect on giant pulses. For giant pulses not to be affected by this process, they must be emitted very close to the light cylinder. We show that the induced process is efficient in the closed field line region, inhibiting propagation of the radio emission in this region.

Coherent synchrotron emission by particles moving along semi-infinite tracks is discussed, with a specific application to radio emission from air showers induced by high-energy cosmic rays. It is shown that in general, radiation from a particle moving along a semi-infinite orbit consists of usual synchrotron emission and modified impulsive bremsstrahlung. The latter component is due to the instantaneous onset of the curved trajectory of the emitting particle at its creation. Inclusion of the bremsstrahlung leads to broadening of the radiation pattern and a slower decay of the spectrum at the cut-off frequency than the conventional synchrotron emission. Possible implications of these features for air shower radio emission are discussed.

The evolution of Alfven turbulence due to three-wave interactions is discussed using kinetic theory for a collisionless, thermal plasma. There are three low-frequency modes, analogous to the three modes of compressible MHD. When only Alfven waves are considered, the known anisotropy of turbulence in incompressible MHD theory is reproduced. Inclusion of a fast mode wave leads to separation of turbulence into two regimes: small wave numbers where three-wave processes involving a fast mode is dominant, and large wave numbers where the three Alfven wave process is dominant. Possible application of the anisotropic Alfven turbulence to the interstellar medium and dissipation of magnetic energy in magnetars is discussed.

We present an empirical model for single pulses of radio emission from pulsars based on gaussian probability distributions for relevant variables. The radiation at a specific pulse phase is represented as the superposition of radiation in two (approximately) orthogonally polarized modes (OPMs) from one or more subsources in the emission region of the pulsar. For each subsource, the polarization states are drawn randomly from statistical distributions, with the mean and the variance on the Poincar\'e sphere as free parameters. The intensity of one OPM is chosen from a log-normal distribution, and the intensity of the other OPM is assumed to be partially correlated, with the degree of correlation also chosen from a gaussian distribution. The model is used to construct simulated data described in the same format as real data: distributions of the polarization of pulses on the Poincar\'e sphere and histograms of the intensity and other parameters. We concentrate on the interpretation of data for specific phases of PSR B0329+54 for which the OPMs are not orthogonal, with one well defined and the other spread out around an annulus on the Poincar\'e sphere at some phases. The results support the assumption that the radiation emerges in two OPMs with closely correlated intensities, and that in a statistical fraction of pulses one OPM is invisible.

Fast spinning magnetars are discussed as strong sources of high energy neutrinos. Pulsars may be born with a short rotation period of milliseconds with the magnetic field amplified through dynamo processes up to $\sim 10^{15}-10^{16} \rm G$. As such millisecond magnetars (MSMs) have an enormous spin-down power $\sim 10^{50} {\rm erg} {\rm s}^{-1}$, they can be potentially a strong, extragalactic high-energy neutrino source. Specifically, acceleration of ions and subsequent photomeson production within the MSM magnetosphere are considered. As in normal pulsars, particle acceleration leads to electron-positron pair cascades that constrains the acceleration efficiency. The limit on the neutrino power as a fraction of the spin-down power is calculated. It is shown that neutrinos produced in the inner magnetosphere have characteristic energy about a few $\times100$ GeV due to the constraint of cooling of charged pions through inverse Compton scattering. TeV neutrinos may be produced in the outer magnetosphere where ions can be accelerated to much higher energies and the pion cooling is less severe than in the inner magnetosphere. High energy neutrinos can also be produced from interactions between ultra-high energy protons accelerated in the magnetosphere and a diffuse thermal radiation from the ejecta or from the interaction region between the MSM wind and remnant shell. The detectability of neutrinos in the early spin-down phase by the current available and planned neutrino detectors is discussed.

A model for the main observational characteristics of the radio emission of pulsars with well organized drifting subpulses is presented. We propose that drifting subpulses result from the modulation of the radio emission mechanism due to long-wavelength drift waves in the magnetosphere. The drift waves are generated at shorter wavelengths, and their nonlinear evolution favours accumulation in a specific azimuthal eigenmode with an integral number, $m$, of nodes encircling the magnetic pole. The electric field of the drift waves is along the magnetic field lines, and this modulates the distribution for particles and hence the radio emission mechanism. The ratio of the frequency of the eigenmode to the rotation frequency of the star is insensitive to the magnetic field strength and the period of rotation, and is of order unity. The period, $P_3$, of the drifting subpulses is attributed to the mismatch between this frequency and the nearest harmonic of the rotation frequency of the star.

A time-dependent model for pair creation in a pulsar magnetosphere is developed. It is argued that the parallel electric field that develops in a charge-starved region (a gap) of a pulsar magnetosphere oscillates with large amplitude. Electrons and positrons are accelerated periodically and the amplitude of the oscillations is assumed large enough to cause creation of upgoing and downgoing pairs at different phases of the oscillation. With a charge-starved initial condition, we find that the oscillations result in bursts of pair creation in which the pair density rises exponentially with time. The pair density saturates at $N_\pm\simeq E_{0}^2/(8\pi m_ec^2\Gamma_{\rm thr})$, where $E_0$ is the parallel electric field in the charge-starved initial state, and $\Gamma_{\rm thr}$ is the Lorentz factor for effec tive pair creation. The frequency of oscillations following the pair creation burst is given roughly by $\omega_{\rm osc}=eE_0/(8m_ec\Gamma_{\rm thr})$. A positive feedback keeps the system stable, such that the average pair creation rate balances the loss rate due to pairs escaping the magnetosphere.

Single pulses of pulsar radio emission are modeled as superposition of radiation originating from many small subsources that are randomly distributed in the emission region. The individual subsources are given an intrinsic finite angular beaming. This model can produce fluctuations in intensity and pulse profiles that are similar to the microstructure observed in some pulsars. Statistics of the phase resolved flux density of a simulated single pulse can be approximated by a lognormal distribution, which is in good agreement with observations.

Refraction of wave propagation in a corotating pulsar magnetospheric plasma is considered as a possible interpretation for observed asymmetric pulse profiles with multiple components. The pulsar radio emission produced inside the magnetosphere propagates outward through the rotating magnetosphere, subject to refraction by the intervening plasma that is spatially inhomogeneous. Both effects of a relativistic distribution of the plasma and rotation on wave propagation are considered. It is shown that refraction coupled with rotation can produce asymmetric conal structures of the profile. The differential refraction due to the rotation can cause the conal structures to skew toward the rotation direction and lead to asymmetry in relative intensities between the leading and trailing components. Both of these features are potentially observable.

Radiative deceleration of relativistic jets in active galactic nuclei as the result of inverse Compton scattering of soft photons from accretion discs is discussed. The Klein-Nishina (KN) cross section is used in the calculation of the radiation force due to inverse Compton scattering. Our result shows that deceleration due to scattering in the KN regime is important only for jets starting with a bulk Lorentz factor larger than 1000. When the bulk Lorentz factor satisfies this condition, particles scattering in the Thomson regime contribute positively to the radiation force (acceleration), but those particles scattering in the KN regime are dominant and the overall effect is deceleration. In the KN limit, the drag due to Compton scattering, though less severe than in the Thomson limit, strongly constrains the bulk Lorentz factor. Most of the power from the deceleration goes into radiation and hence the ability of the jet to transport significant power (in particle kinetic energy) out of the subparsec region is severely limited. The deceleration efficiency decreases significantly if the jet contains protons and the proton to electron number density ratio satisfies the condition $n_p/n_{e0}>2\gamma_{\rm min}/\mu_p$ where $\gamma_{\rm min}$ is the minimum Lorentz factor of relativistic electrons (or positrons) in the jet frame and $\mu_p$ is the proton to electron mass ratio.

It is shown that for moderately hot polar caps (with effective temperature of $\sim 10^6$ K), the efficiency of polar gap acceleration is lower compared to the case in which the polar caps are relatively cool and inverse Compton scattering plays no role in controlling the gap. For young pulsars with superstrong magnetic fields ($\geq10^9 T$) and hot polar caps (with temperature of $\ge 5\times10^6 K$), because of the energy loss of electrons or positrons due to resonant inverse Compton scattering in the vicinity of polar caps, pair cascades occur at distances further away from the polar cap, and in this case, we have a relatively high acceleration efficiency with ions carrying most of the particle luminosity.

We consider the result of acceleration of heavy ions in the slot gap potential of a very young pulsar with a hot polar cap. Photodisintegration of the heavy ions in the radiation field of the polar cap and pulsar surface gives rise to a flux of energetic neutrons. Some fraction of these neutrons interact with target nuclei in the supernova shell to produce neutrino and gamma-ray signals which should be observable from very young supernova remnants in our galaxy for a range of pulsar parameters.