Imposing an external magnetic field in short-pulse intense laser-plasma interaction is of broad scientific interest in related plasma research areas. We propose a simple method using a virtual current layer by introducing an extra current density term to simulate the external magnetic field, and demonstrate it with three-dimensional particle-in-cell simulations. The field distribution and its evolution in sub-picosecond time scale are obtained. The magnetization process takes a much longer time than that of laser-plasma interaction due to plasma diamagnetism arising from collective response. The long-time evolution of magnetic diffusion and diamagnetic current can be predicted based on a simplified analytic model in combination with simulations.

The bright radio pulsar B1727$-$47 with a characteristic age of 80 kyr is among the first pulsars discovered 50 yr ago. Using its regular timing observations and interferometric positions at three epochs, we measured, for the first time, the pulsar proper motion of $151 \pm 19$ mas yr$^{-1}$. At the dispersion measure distance of $\gtrsim 2.7$ kpc, this would suggest a record transverse velocity of the pulsar $\gtrsim 1900$ km s$^{-1}$. However, a backward extrapolation of the pulsar track to its birth epoch points remarkably close to the center of the evolved nearby supernova remnant RCW 114, which suggests genuine association of the two objects. In this case, the pulsar is substantially closer ($\sim 0.6$ kpc) and younger ($\sim 50$ kyr), and its velocity ($\sim400$ km s$^{-1}$) is compatible with the observed pulsar velocity distribution. We also identified two new glitches of the pulsar. We discuss implications of our results on the pulsar and remnant properties.

PSR B0919+06 is known for its abnormal emission phenomenon, where the pulse emission window occasionally shifts progressively in longitude and returns afterwards. The physical mechanism behind this phenomenon is still under investigation. In this paper, we present our ultra-wideband observation of this pulsar using the Five-hundred-meter Aperture Spherical radio Telescope (FAST), with simultaneous measurements in the frequency ranges 280-780 MHz and 1250-1550 MHz. We have identified three abnormal events, each of which becomes less apparent as the frequency decreases. At 1400 MHz, the averaged profile slightly shifted after the first and third abnormal events, implying a relationship between abnormal event and profile variation. We also found a linear trend in the left-edge position of the averaged profiles between the first and third event as well as after the third event, suggesting the existence of a slow-drifting mode between the two major events. The second event has a comparatively small shift in phase and is thus categorized as a `small flare state'. During the third event, a sequence of approximately nine pulses was seen to significantly weaken in all frequency bands, likely associated with the pseudo-nulling observed at 150 MHz. A three-component de-composition analysis of the normal averaged profiles shows that the trailing component is dominant at our observing frequencies, while the centre component has a comparatively steeper spectrum. We found the overall flux density in an abnormal event to slightly differ from that in an ordinary state, and the difference shows a frequency dependence. A comparison of the normal, abnormal and dimmed averaged profile indicates that the leading component is likely to be stable in all states.

With the largest dish Five-hundred-meter Aperture Spherical radio Telescope (FAST), both the mean and single pulses of PSR B2016$+$28, especially including the single-pulse structure, are investigated in detail in this study. The mean pulse profiles at different frequencies can be well fitted in a conal model, and the peak separation of intensity-dependent pulse profiles increases with intensity. The integrated pulses are obviously frequency dependent (pulse width decreases by $\sim\,20\%$ as frequency increases from 300 MHz to 750 MHz), but the structure of single pulses changes slightly (the corresponding correlation scale decreases by only $\sim\,1\%$). This disparity between mean and single pulses provides independent evidence for the existence of the RS-type vacuum inner gap, indicating a strong bond between particles on the pulsar surface. Diffused drifting sub-pulses are analyzed. The results show that the modulation period along pulse series ($P_3$) is positively correlated to the separation between two adjacent sub-pulses ($P_2$). This correlation may hint a rough surface on the pulsar, eventually resulting in the irregular drift of sparks. All the observational results may have significant implications in the dynamics of pulsar magnetosphere and are discussed extensively in this paper.

Rotating radio transients (RRATs) are peculiar astronomical objects whose emission mechanism remains under investigation. In this paper, we present observations of three RRATs, J1538+2345, J1854+0306 and J1913+1330, observed with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Specifically, we analyze the mean pulse profiles and temporal flux density evolutions of the RRATs. Owing to the high sensitivity of FAST, the derived burst rates of the three RRATs are higher than those in previous reports. RRAT J1854+0306 exhibited a time-dynamic mean pulse profile, whereas RRAT J1913+1330 showed distinct radiation and nulling segments on its pulse intensity trains. The mean pulse profile variation with frequency is also studied for RRAT J1538+2345 and RRAT J1913+1330, and the profiles at different frequencies could be well fitted with a cone-core model and a conal-beam model, respectively.

On the strength of a double Nobel prize winning experiment (Super)Kamiokande and an extremely successful long baseline neutrino programme, the third generation Water Cherenkov detector, Hyper-Kamiokande, is being developed by an international collaboration as a leading worldwide experiment based in Japan. The Hyper-Kamiokande detector will be hosted in the Tochibora mine, about 295 km away from the J-PARC proton accelerator research complex in Tokai, Japan. The currently existing accelerator will be steadily upgraded to reach a MW beam by the start of the experiment. A suite of near detectors will be vital to constrain the beam for neutrino oscillation measurements. A new cavern will be excavated at the Tochibora mine to host the detector. The experiment will be the largest underground water Cherenkov detector in the world and will be instrumented with new technology photosensors, faster and with higher quantum efficiency than the ones in Super-Kamiokande. The science that will be developed will be able to shape the future theoretical framework and generations of experiments. Hyper-Kamiokande will be able to measure with the highest precision the leptonic CP violation that could explain the baryon asymmetry in the Universe. The experiment also has a demonstrated excellent capability to search for proton decay, providing a significant improvement in discovery sensitivity over current searches for the proton lifetime. The atmospheric neutrinos will allow to determine the neutrino mass ordering and, together with the beam, able to precisely test the three-flavour neutrino oscillation paradigm and search for new phenomena. A strong astrophysical programme will be carried out at the experiment that will detect supernova neutrinos and will measure precisely solar neutrino oscillation.

Practical application of the harmonic summing technique in the power-spectrum analysis for searching pulsars has exhibited the technique's effectiveness. In this paper, theoretical verification of harmonic summing considering power's noise-signal probability distribution is given. With the top-hat and the modified von Mises pulse profile models, contours along which spectra total power is expected to exceed the 3\,$\sigma$ detection threshold with 0.999 confidence corresponding to $m=1, 2, 4, 8, 16$, or 32 harmonics summed are given with respect to the mean pulse amplitude and the pulse duty cycle. Optimized numbers of harmonics summed relative to the duty cycles are given. The routine presented builds a theoretical estimate of the minimum detectable mean flux density, i.e. sensitivity, under the power-spectrum searching method.

Determination of detection sensitivity in a number of previous pulsar search programmes was done via the straightfoward use of the radiometer equation. In the same surveys, the Fourier domain method was used to search for pulsars. As detection sensitivity is partially a function of the searching method, the straightfoward use of the radiometer equation for detection sensitivity determination is not consistent with the Fourier searching method. In this proceeding, I clarify the problem and note the way for sensitivity determination which is consistent with the Fourier searching method. More details can be found in Yu (2018).

Based on the interstellar CO/H2 ratio, carbon monoxide-based censuses of protoplanetary disks in Lupus, sigma Orionis, and Chamaeleon I found no disks more massive than the minimum-mass solar nebula, which is inconsistent with the existence of exoplanets more massive than Jupiter. Observations and models are converging on the idea that ionization-driven chemistry depletes carbon monoxide in T-Tauri disks. Yet the extent of CO depletion depends on the incident flux of ionizing radiation, and some T-Tauri stars may have winds strong enough to shield their disks from cosmic rays. There is also a range of X-ray luminosities possible for a given stellar mass. Here we use a suite of chemical models, each with a different incident X-ray or cosmic-ray flux, to assess whether CO depletion is a typical outcome for T-Tauri disks. We find that CO dissociation in the outer disk is a robust result for realistic ionization rates, with abundance reductions between 70% and 99.99% over 2 Myr of evolution. Furthermore, after the initial dissociation epoch, the inner disk shows some recovery of the CO abundance from CO2 dissociation. In highly ionized disks, CO recovery in the inner disk combined with depletion in the outer disk creates a centrally peaked CO abundance distribution. The emitting area in rare CO isotopologues may be an indirect ionization indicator: in a cluster of disks with similar ages, those with the most compact CO isotopologue emission see the highest ionization rates.

Neutron star glitches are observed via artificially scheduled pulsar pulse arrival-time observations. Detection probability density of glitch events for a given data set is essentially required knowledge for realizing glitch detectability with specified observing system and schedule. In Yu & Liu, the detection probability density was derived for the Yu et al. data set. In this proceeding, further discussions are presented.

We report preliminary results of the analysis of the proper motion of the bright radio pulsar B1727$-$47. Using archival Parkes timing data, as well as original and archival ATCA interferometry observations, we, for the first time, constrain the pulsar proper motion at the level of 148$\pm$11 mas yr$^{-1}$. The backward extrapolation of the proper motion vector to the pulsar birth epoch points at the center of the Galactic supernova remnant RCW 114 suggesting the genuine association between the two objects. We discuss the implications of the association and argue that the distance to the system is less than $1$ kpc. This value is at least two times lower than the dispersion measure distance estimates. This suggests that the existing Galaxy electron density models are incomplete in the direction to the pulsar.

Turbulence is the leading candidate for angular momentum transport in protoplanetary disks and therefore influences disk lifetimes and planet formation timescales. However, the turbulent properties of protoplanetary disks are poorly constrained observationally. Recent studies have found turbulent speeds smaller than what fully-developed MRI would produce (Flaherty et al. 2015, 2017). However, existing studies assumed a constant CO/H2 ratio of 0.0001 in locations where CO is not frozen-out or photo-dissociated. Our previous studies of evolving disk chemistry indicate that CO is depleted by incorporation into complex organic molecules well inside the freeze-out radius of CO. We consider the effects of this chemical depletion on measurements of turbulence. Simon et al. (2015) suggested that the ratio of the peak line flux to the flux at line center of the CO J=3-2 transition is a reasonable diagnostic of turbulence, so we focus on that metric, while adding some analysis of the more complex effects on spatial distribution. We simulate the emission lines of CO based on chemical evolution models presented in Yu et al. (2016), and find that the peak-to-trough ratio changes as a function of time as CO is destroyed. Specifically, a CO-depleted disk with high turbulent velocity mimics the peak-to-trough ratios of a non-CO-depleted disk with lower turbulent velocity. We suggest that disk observers and modelers take into account the possibility of CO depletion when using line peak-to-trough ratios to constrain the degree of turbulence in disks. Assuming that CO/H2 = 0.0001 at all disk radii can lead to underestimates of turbulent speeds in the disk by at least 0.2 km/s.

Neutron stars are observed to undergo small, abrupt rotational speed-up. This phenomenon is known as glitch. In pulsar timing observations, detection of a neutron star glitch is constrained by the time of occurrence of the event relative to entire observing span and observing cadences, time of occurrence of preceding/subsequent glitches relative to observing cadences and the strength of timing noise. Using the Yu et al. (2013) data sets, in this paper, we analyse the observational selection in terms of detection probability. We define partial probabilities for the constraints and use Monte Carlo method with assuming glitches distribute uniformly to solve the complete probability formula for both group case involving 157 pulsars and individual cases for each of the seven pulsars with glitch numbers $\geqslant 5$. In the simulations, numerical Bayesian inference is used for glitch identification. With the derived detection probability density and observed results, we uncover glitch size probability distribution embedded in the data for both the group and individual cases. We find the most prominent correction occurred for PSR~J1341$-$6220, in which exponent of the power-law model varies from the observed $+0.7^{+1.4}_{-0.7}$ to $-0.4^{+1.0}_{-0.4}$. We suggest observers determine the detection probability for glitch theories, e.g. the self-organised criticality.

Gas in protostellar disks provides the raw material for giant planet formation and controls the dynamics of the planetesimal-building dust grains. Accurate gas mass measurements help map the observed properties of planet-forming disks onto the formation environments of known exoplanets. Rare isotopologues of carbon monoxide (CO) have been used as gas mass tracers for disks in the Lupus star-forming region, with an assumed interstellar CO/H$_2$ abundance ratio. Unfortunately, observations of T-Tauri disks show that CO abundance is not interstellar---a finding reproduced by models that show CO abundance decreasing both with distance from the star and as a function of time. Here we present radiative transfer simulations that assess the accuracy of CO-based disk mass measurements. We find that the combination of CO chemical depletion in the outer disk and optically thick emission from the inner disk leads observers to underestimate gas mass by more than an order of magnitude if they use the standard assumptions of interstellar CO/H$_2$ ratio and optically thin emission. Furthermore, CO abundance changes on million-year timescales, introducing an age/mass degeneracy into observations. To reach factor of a few accuracy for CO-based disk mass measurements, we suggest that observers and modelers adopt the following strategies: (1) select the low-$J$ transitions; (2) observe multiple CO isotopologues and use either intensity ratios or normalized line profiles to diagnose CO chemical depletion; and (3) use spatially resolved observations to measure the CO abundance distribution.

The gas near the midplanes of planet-forming protostellar disks remains largely unprobed by observations due to the high optical depth of commonly observed molecules such as CO and H$_2$O. However, rotational emission lines from rare molecules may have optical depths near unity in the vertical direction, so that the lines are strong enough to be detected, yet remain transparent enough to trace the disk midplane. Here we present a chemical model of an evolving T-Tauri disk and predict the optical depths of rotational transitions of $^{12}$C$^{16}$O, $^{13}$C$^{16}$O, $^{12}$C$^{17}$O and $^{12}$C$^{18}$O. The MRI-active disk is primarily heated by the central star due to the formation of the dead zone. CO does not freeze out in our modeled region within $70$AU around a sunlike star. However, the abundance of CO decreases because of the formation of complex organic molecules (COM), producing an effect that can be misinterpreted as the "snow line". These results are robust to variations in our assumptions about the evolution of the gas to dust ratio. The optical depths of low-order rotational lines of C$^{17}$O are around unity, making it possible to see into the disk midplane using C$^{17}$O. Combining observations with modeled C$^{17}$O$/$H$_2$ ratios, like those we provide, can yield estimates of protoplanetary disks' gas masses.

Alfv\'{e}nic turbulent cascade perpendicular and parallel to the background magnetic field is studied accounting for anisotropic dispersive effects and turbulent intermittency. The perpendicular dispersion and intermittency make the perpendicular-wavenumber magnetic spectra steeper and speed up production of high ion-cyclotron frequencies by the turbulent cascade. On the contrary, the parallel dispersion makes the spectra flatter and decelerate the frequency cascade above the ion-cyclotron frequency. Competition of the above factors results in spectral indices distributed in the interval [-2,-3], where -2 is the index of high-frequency space-filling turbulence, and -3 is the index of low-frequency intermittent turbulence formed by tube-like fluctuations. Spectra of fully intermittent turbulence fill a narrower range of spectral indices [-7/3,-3], which almost coincides with the range of indexes measured in the solar wind. This suggests that the kinetic-scale turbulent spectra are shaped mainly by dispersion and intermittency. A small mismatch with measured indexes of about 0.1 can be associated with damping effects not studied here.

Linear properties of kinetic Alfv\'{e}n waves (KAWs) and kinetic slow waves (KSWs) are studied in the framework of two-fluid magnetohydrodynamics. We obtain the wave dispersion relations that are valid in a wide range of the wave frequency $\omega $ and plasma-to-magnetic pressure ratio $\beta $. The KAW frequency can reach and exceed the ion cyclotron frequency at ion kinetic scales, whereas the KSW frequency remains sub-cyclotron. At $\beta \sim 1$, the plasma and magnetic pressure perturbations of both modes are in anti-phase, so that there is nearly no total pressure perturbations. However, these modes exhibit also several opposite properties. At high $% \beta $, the electric polarization ratios of KAWs and KSWs are opposite at the ion gyroradius scale, where KAWs are polarized in sense of electron gyration (right-hand polarized) and KSWs are left-hand polarized. The magnetic helicity $\sigma \sim 1$ for KAWs and $\sigma \sim -1$ for KSWs, and the ion Alfv\'{e}n ratio $R_{Ai}\ll 1$ for KAWs and $R_{Ai}\gg 1$ for KSWs. We also found transition wavenumbers where KAWs change their polarization from left- to right-hand. These new properties can be used to discriminate KAWs and KSWs when interpreting kinetic-scale electromagnetic fluctuations observed in various solar-terrestrial plasmas. This concerns, in particular, identification of modes responsible for kinetic-scale pressure-balanced fluctuations and turbulence in the solar wind.

Parkes pulse arrival-time data for 165 radio pulsars spanning from 1990 to 2011 have been searched for period glitches. Forty-six events out of the detected 107 glitches were found to be new contributions to the entire glitch population which currently contains approximately 400 events.

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.

Reconnection of the self-generated magnetic fields in laser-plasma interaction was first investigated experimentally by Nilson {\it et al.} [Phys. Rev. Lett. 97, 255001 (2006)] by shining two laser pulses a distance apart on a solid target layer. An elongated current sheet (CS) was observed in the plasma between the two laser spots. In order to more closely model magnetotail reconnection, here two side-by-side thin target layers, instead of a single one, are used. It is found that at one end of the elongated CS a fan-like electron outflow region including three well-collimated electron jets appears. The ($>1$ MeV) tail of the jet energy distribution exhibits a power-law scaling. The enhanced electron acceleration is attributed to the intense inductive electric field in the narrow electron dominated reconnection region, as well as additional acceleration as they are trapped inside the rapidly moving plasmoid formed in and ejected from the CS. The ejection also induces a secondary CS.

We present a theoretical model for the thermal X-ray emission and cooling of isolated pulsars, assuming that pulsars are solid quark stars. We calculate the heat capacity for such a quark star, and the results show that the residual thermal energy cannot sustain the observed thermal X-ray luminosities seen in typical isolated X-ray pulsars. We conclude that other heating mechanisms must be in operation if the pulsars are in fact solid quark stars. Two possible heating mechanisms are explored. Firstly, for pulsars with little magnetospheric activities, accretion from the interstellar medium or from the material in the associated supernova remnants may power the observed thermal emission. In the propeller regime, a disk-accretion rate ${\dot M}\sim$1% of the Eddington rate with an accretion onto the stellar surface at a rate of $\sim 0.1% {\dot M}$ could explain the observed emission luminosities of the dim isolated neutron stars and the central compact objects. Secondly, for pulsars with significant magnetospheric activities, the pulsar spindown luminosities may have been as the sources of the thermal energy via reversing plasma current flows. A phenomenological study between pulsar bolometric X-ray luminosities and the spin energy loss rates presents the probable existence of a 1/2-law or a linear law, i.e. $L_{\rm bol}^{\infty}\propto\dot{E}^{1/2}$ or $L_{\rm bol}^{\infty}\propto\dot{E}$. This result together with the thermal properties of solid quark stars allow us to calculate the thermal evolution of such stars. Thermal evolution curves, or cooling curves, are calculated and compared with the `temperature-age' data obtained from 17 active X-ray pulsars. It is shown that the bolometric X-ray observations of these sources are consistent with the solid quark star pulsar model.

We present here a phenomenological solid quark star pulsar model to interpret the observed thermal X-ray emission of isolated pulsars. The heat capacity for solid quark stars was found to be quite small, so that the residual internal stellar heat gained at the birth of the star could be dissipated in an extremely short timescale. However, the bombardment induced by backflowing plasma at the poles of solid quark stars would get the stars be reheated, so that long term soft X-ray emission can be sustained. Such a scenario could be used for those X-ray pulsars with significant magnetospheric activities, and their cooling processes would thus be established. Dim X-ray isolated neutron stars (XDINs) as well as compact central objects (CCOs) have been observed with dominant soft X-ray radiation combined with little magnetospheric manifestations. Such sources could be solid quark stars accreting in the propeller regime.