Rotating Massive Strangeon Stars and X-Ray Plateau of Short GRBs

Strangeon stars,which are proposed to describe the nature of pulsar-like compact stars,have passed various observational tests.The maximum mass of a non-rotating strangeon star could be high,which implies that the remnants of binary strangeon star mergers could even be long-lived massive strangeon stars.We study rigidly rotating strangeon stars in the slowly rotating approximation,using the Lennard-Jones model for the equation of state.Rotation can significantly increase the maximum mass of strangeon stars with unchanged baryon numbers,enlarging the mass-range of long-lived strangeon stars.During spin-down after merger,the decrease of radius of the remnant will lead to the release of gravitational energy.Taking into account the efficiency of converting the gravitational energy luminosity to the observed X-ray luminosity,we find that the gravitational energy could provide an alternative energy source for the plateau emission of X-ray afterglow.The fitting results of X-ray plateau emission of some short gamma-ray bursts suggest that the magnetic dipole field strength of the remnants can be much smaller than that of expected when the plateau emission is powered only by spin-down luminosity of magnetars.

Free Energy of Anisotropic Strangeon Stars

Can pulsar-like compact objects release further huge free energy besides the kinematic energy of rotation?This is actually relevant to the equation of state of cold supra-nuclear matter,which is still under hot debate.Enormous energy is surely needed to understand various observations,such asγ-ray bursts,fast radio bursts and softγ-ray repeaters.In this paper,the elastic/gravitational free energy of solid strangeon stars is revisited for strangeon stars,with two anisotropic models to calculate in general relativity.It is found that huge free energy(>10^(46)erg)could be released via starquakes,given an extremely small anisotropy((p_(t)-p_(r))/p_(r)~10^(-4),with pt/pr the tangential/radial pressure),implying that pulsar-like stars could have great potential of free energy release without extremely strong magnetic fields in the solid strangeon star model.

Probing the Internal Physics of Neutron Stars through the Observed Braking Indices and Magnetic Tilt Angles of Several Young Pulsars

The braking indices of pulsars may contain important information about the internal physics of neutron stars(NSs),such as neutron superfluidity and internal magnetic fields.As a subsequent paper of Cheng et al.,we perform the same analysis as that done in the previous paper to other young pulsars with a steady braking index,n.Combining the timing data of these pulsars with the theory of magnetic field decay,and using their measured magnetic tilt angles,we can set constraints on the number of precession cycles,ξ,which represents the interactions between superfluid neutrons and other particles in the NS interior.For the pulsars considered in this paper,the results show thatξis within the range of a few×10~3 to a few×10~6.Interestingly,for the Crab and Vela pulsars,the constraints onξobtained with our method are generally consistent with that derived from modeling of the glitch rise behaviors of the two pulsars.Furthermore,we find that the internal magnetic fields of pulsar with n<3 may be dominated by the toroidal components.Our results may not only help to understand the interactions between the superfluid neutrons and other particles in the interior of NSs but also be important for the study of continuous gravitational waves from pulsars.

The EOSs and the Blatant Discrepancy in Modelling Massive Neutron Stars: Origin and a Possible Solution Method

Exploring the state of ultra-cold supranuclear dense matter that makes up the cores of massive neutron stars is one of the greatest unresolved problems in modern physics. In this letter, we show that when the interiors of pulsars are made of compressible and dissipative normal matter, the commonly used solution procedures combined with the known EOSs yield widely scattered solutions and poorly determined radii. A remarkable agreement emerges, however, if pulsars harbour cores that are made of incompressible entropy-free superfluids (SuSu-matter) embedded in flat spacetimes. Such supranuclear dense matter should condensate to form false vacua as predicated by non-perterbative QCD vacuum. The solutions here are found to be physically consistent and mathematically elegant, irrespective of the object’s mass. Based thereon, we conclude that the true masses of massive NSs may differ significantly from those revealed by direct observation.

A Fermi-LAT Study of Globular Cluster Dynamical Evolution in the Milky Way:Millisecond Pulsars as the Probe

Using archival Fermi-LAT data with a time span of~12 yr,we study the population of Millisecond Pulsars(MSPs)in Globular Clusters(GlCs)and investigate their dependence on cluster dynamical evolution in the Milky Way.We show that theγ-ray luminosity(L_(γ))and emissivity(i.e.,ε_(γ)=L_(γ)/M,with M the cluster mass)are good indicators of the population and abundance of MSPs in GlCs,and they are highly dependent on the dynamical evolution history of the host clusters.Specifically speaking,the dynamically older GlCs with more compact structures are more likely to have larger L_(γ)andε_(γ),and these trends can be summarized as strong correlations with cluster stellar encounter rateΓand the specific encounter rate(Λ=Γ/M),with L_(γ)∝Γ^(0.7±0.11)andε_(γ)∝Λ^(0.73±0.13)for dynamically normal GlCs.However,as GlCs evolve into deep core collapse,these trends are found to be reversed,implying that strong encounters may have lead to the disruption of Low-Mass X-ray Binaries and ejection of MSPs from core-collapsed systems.Besides,the GlCs are found to exhibit largerε_(γ)with increasing stellar mass function slope(ε_(γ)∝10^((0.52±0.1)α)),decreasing tidal radius(ε_(γ)∝R_(t)^(-10±0.22))and distances from the Galactic Center(GC,ε_(γ)∝R_(gc)^(-1.13±0.21)).These correlations indicate that,as GlCs losing kinetic energy and spiral in toward the GC,tidal stripping and mass segregation have a preference in leading to the loss of normal stars from GlCs,while MSPs are more likely to concentrate to cluster center and be deposited into the GC.Moreover,we gaugeε_(γ)of GlCs is~10-1000 times larger than the Galactic bulge,the latter is thought to reside thousands of unresolved MSPs and may be responsible for the GC 7-ray excess,which supports that GlCs are generous contributors to the population of MSPs in the GC.

The Origin of Glitches in Pulsars—Phase Oscillation between Anisotropic Superfluid and Normal State of Neutrons in Neutron Stars

Considering neutron star heating by magnetic dipole radiation from ^3PF2 superfluid neutron vortices inside the star, we propose a neutron phase oscillation model between the normal neutron Fermi fluid and the ^3PF2 superfluid neutron vortices at the transition temperature of Ttrans = （2 - 3）× 10^8 K. With this model we can qualitatively explain most of the observations on pulsar glitches up to date.

Spin Evolution of Neutron Stars in OB/X-ray Binaries

We have investigated the relation between the orbital period Porb and the spin period Ps of neutron stars in OB/X-ray binaries. By simulating the time-development of the mass loss rate and radius expansion of a 20M☉ donor star, we have calculated the detailed spin evolution of the neutron star before steady wind accretion occurs (that is, when the break spin period is reached), or when the OB star begins evolving off the main sequence or has filled its Roche lobe. Our results are compatible with the observations of OB/X-ray binaries. We find that in relatively narrow systems with orbital periods less than tens of days, neutron stars with initial magnetic field B0 stronger than about 3×1012 G can reach the break spin period to allow steady wind accretion in the main sequence time, whereas neutron stars with B0 < 3×1012 G and/or in wide systems would still be in one of the pulsar, rapid rotator or propeller phases when the companion evolves off the main sequence or fills its Roche lobe. Our results may help understand the various characteristics of the observed OB/neutron star binaries along with their distributions in the Ps-Porb diagram.

Resonant cyclotron scattering in pulsar magnetospheres and its application to isolated neutron stars

Resonant cyclotron scattering （RCS） in pulsar magnetospheres is considered. The photon diffusion equation （Kompaneets equation） for RCS is derived. The photon system is modeled three dimensionally. Numerical calculations show that there exist not only up scattering but also down scattering of RCS, depending on the parameter space. RCS＇s possible applications to spectral energy distributions of magnetar candidates and radio quiet isolated neutron stars （INSs） are pointed out. The optical/UV excess of INSs may be caused by the down scattering of RCS. The calculations for RX J1856.5-3754 and RX J0720.4-3125 are presented and compared with their observational data. In our model, the INSs are proposed to be normal neutron stars, although the quark star hypothesis is still possible. The low pulsation amplitude of INSs is a natural consequence in the RCS model.

Merging strangeon stars

The state of supranuclear matter in compact stars remains puzzling, and it is argued that pul- sars could be strangeon stars. What would happen if binary strangeon stars merge？ This kind of merger could result in the formation of a hyper-massive strangeon star, accompanied by bursts of gravitational waves and electromagnetic radiation （and even a strangeon kilonova explained in the paper）. The tidal polarizability of binary strangeon stars is different from that of binary neutron stars, because a strangeon star is self-bound on the surface by the fundamental strong force while a neutron star by the gravity, and their equations of state are different. Our calculation shows that the tidal polarizability of merging bi- nary strangeon stars is favored by GW170817. Three kinds of kilonovae （i.e., of neutron, quark and strangeon） are discussed, and the light curve of the kilonova AT 2017gfo following GW170817 could be explained by considering the decaying strangeon nuggets and remnant star spin-down. Additionally, the energy ejected to the fireball around the nascent remnant strangeon star, being manifested as a gamma-ray burst, is calculated. It is found that, after a prompt burst, an X-ray plateau could follow in a timescale of 102 - 103 s. Certainly, the results could be tested also by further observational synergies between gravitational wave detectors （e.g., Advanced LIGO） and X-ray telescopes （e.g., the Chinese HXMT satellite and eXTP mission）, and especially if the detected gravitational wave form is checked by peculiar equations of state provided by the numerical relativistical simulation.

Systematic comparison of initial velocities for neutron stars in different models

I have studied the initial velocity（Maxwellian and exponential distributions） and the scale height of isolated old（aged≥10^9yr） neutron stars（NSs） at different Galactocentric distances R in three population models. The smooth time-independent 3-D axisymmetric gravitational potentials（MiyamotoNagai and Paczy n′ski models） were also used. The correlation between these quantities significantly affects the shapes of the profiles and distributions of the simulated sample, because the differences in the initial kick can arise from differences in the formation and evolution of NSs with other physical parameters. The scale height of the density distribution increases systematically with R. I have also shown that the distribution of old NSs in these population models agrees with the observed structure of the Galaxy in terms of initial velocities（1-D and 3-D）, as well as the scale height distributions. These distributions tend to have an asymptotic behavior at the point R = 2.75 kpc. This means that the quality of the models can be described in terms of a mean of the fitted Gaussian, and this could also give an overall perspective of the phase space properties of nearby old NSs on a given gravitational potential.

Singlet pairing gaps of neutrons and protons in hyperonic neutron stars

The ^1 S0 nucleonic superfluids are investigated within the relativistic meanfield model and Bardeen-Cooper-Schrieffer theory in hyperonic neutron stars. The ^1 S0 pairing gaps of neutrons and protons are calculated based on the Reid soft-core interaction as the nucleon-nucleon interaction. In particular, we have studied the influence of degrees of freedom for hyperons on the ^1 S0 nucleonic pairing gap in neutron star matter. It is found that the appearance of hyperons has little impact on the baryonic density range and the size of the ^1S0 neutronic pairing gap; the ^1S0 protonic pairing gap also decreases slightly in this region where ρB = 0.0-0.393 fm^-3. However, if baryonic density becomes greater than 0.393 fm^-3, the ^1S0 protonic pairing gap obviously increases. In addition, the possible range for a protonic superfluid is obviously enlarged due to the presence of hyperons. In our results, the hyperons change the 1 So protonic pairing gap, which must change the cooling properties of neutron stars.

Microscopic Magnetic Dipole Radiation in Neutron Stars

There is a ^3P2 neutron superfluid region in NS （neutron star） interior. For a rotating NS the ^3P2 superfluid region is like a system of rotating magnetic dipoles. It will give out electromagnetic radiation, which may provide a new heating mechanism of NSs. This mechanism plus some cooling agent may give a sound explanation to NS glitches.

Transient Phenomena from Gas Erupting Neutron Stars

Starquakes probably occur in rapidly spinning or ultra high field neutron stars. In this article, we argue that highly compressed gas containing electron–positron pairs could evaporate and erupt from inside the neutron star when a crack forms and then heals during a starquake. Under the influence of the existing oscillation modes of the star, the crack may have sufficiently large size and long lifetime. An appropriate amount of gas can erupt into the magnetosphere with relativistic and nearly uniform velocity, producing various transient and bursting phenomena.

The optical/ultraviolet excess of isolated neutron stars in the resonant cyclotron scattering model

X-ray dim isolated neutron stars are peculiar pulsar-like objects, characterized by their Planck-like spectrum. In studying their spectral energy distributions, optical/ultraviolet （UV） excess is a long standing problem. Recently Kaplan et al. measured the optical/UV excess for all seven sources, which is understandable in the resonant cyclotron scattering （RCS） model previously addressed. The RCS model calculations show that the RCS process can account for the observed optical/UV excess for most sources. The flat spectrum of RX J2143.0＋0654 may be due to contributions from the bremsstrahlung emission of the electron system in addition to the RCS process.

First-order phase transitions in rotating hybrid stars and pulsar glitches

The first order deconfinement phase transitions in rotating hybrid stars are studied and it is found that if the surface tension is sufficiently large, the transition from metastable hadron matter to stable mixed hadron-quark matter during the spindown history of a hybrid star can cause a glitch.

Single-pulse Emission Variation of Two Pulsars Discovered by FAST

We investigate the single-pulse emission variations of two pulsars,PSRs J0211+4235 and J0553+4111,observed with the Five-hundred-meter Aperture Spherical radio Telescope at the 1.25 GHz central frequency.The observation sessions span from 2020 December to 2021 July,with 21 and 22 observations for them respectively.The integrated pulse profile of PSR J0211+4235 shows that there is a weak pulse component following the main component,and PSR J0553+4111 displays a bimodal profile with a bridge component in the middle.PSR J0211+4235 presents significant nulling phenomenon with nulling duration lasting from 2 to 115 pulses and burst duration lasting from 2 to 113 pulses.The NF of each observation is determined to be 45%-55%.No emission greater than threeσis found in the mean integrated profile of all nulling pulses.In most cases,the pulse energy changes abruptly during the transition from null to burst,while in the transition from burst to null there are two trends:abrupt and gradual.We find that the nulling phenomenon of PSR J0211+4235 is periodic by the Fourier transform of the null and burst state.In addition,the single-pulse modulation characteristics of these two pulsars are investigated,and the distributions of modulation index,LRFS and 2DFS are analyzed with PSRSALSA.The left peak of PSR J0553+4111 has intensity modulation.Finally,the polarization properties of these two pulsars are obtained through polarization calibration,and their characteristics are analyzed.The possible physical mechanisms of these phenomena are discussed.

Searching for Pulsars with Phase Characteristics

We present a method by using the phase characteristics of radio observation data for pulsar search and candidate identification.The phase characteristics are relations between the pulsar signal and the phase correction in the frequency-domain,and we regard it as a new search diagnostic characteristic.Based on the phase characteristics,a search method is presented:calculating dispersion measure(DM)—frequency data to select candidate frequencies,and then confirming of candidates by using the broadband characteristics of pulsar signals.Based on this method,we performed a search test on short observation data of M15 and M71,which were observed by Five-hundredmeter Aperture Spherical radio Telescope,and some of the Galactic Plane Pulsar Snapshot survey data.Results show that it can get similar search results to PRESTO(Pulsa R Exploration and Search TOolkit)while having a faster processing speed.

The Analyses of Globular Cluster Pulsars and Their Detection Efficiency

Up to 2022 November,267 pulsars had been discovered in 36 globular clusters(GCs).In this paper,we present our studies on the distribution of GC pulsar parameters and the detection efficiency.The power law relation between average dispersion measure(■)and dispersion measure difference(ΔDM)of known pulsars in GCs is lgΔDM∝1.52lg■.The sensitivity could be the key to finding more pulsars.As a result,several years after the construction of a large radio telescope facility,the number of known GC pulsars will likely be increased accordingly.We suggest that currently GCs in the southern hemisphere could have higher possibilities for finding new pulsars.

FAST Pulsar Database.I.Polarization Profiles of 682 Pulsars

Pulsar polarization profiles form a very basic database for understanding the emission processes in a pulsar magnetosphere.After careful polarization calibration of the 19-beam L-band receiver and verification of beamoffset observation results,we obtain polarization profiles of 682 pulsars from observations by the Five-hundredmeter Aperture Spherical radio Telescope(FAST)duringthe Galactic Plane Pulsar Snapshot survey and other normal FAST projects.Among them,polarization profiles of about 460 pulsars are observed for the first time.The profiles exhibit diverse features.Some pulsars have a polarization position angle curve with a good S-shaped swing,some with orthogonal modes;some have components with highly linearly polarized components or strong circularly polarized components;some have a very wide profile,coming from an aligned rotator,and some have an interpulse from a perpendicular rotator;some wide profiles are caused by interstellar scattering.We derive geometric parameters for 190 pulsars from the S-shaped position angle curves or with orthogonal modes.We find that the linear and circular polarization or the widths of pulse profiles have various frequency dependencies.Pulsars with a large fraction of linear polarization are more likely to have a large Edot.

Constraining the Orbital Inclination and Companion Properties of Three Black Widow Pulsars Detected by FAST

Black widows(BWs)are millisecond pulsars ablating their companion stars.The out-flowing material from the companion can block the radio emission of the pulsar,resulting in eclipses.In this paper,we construct a model for the radio eclipse by calculating the geometry of the bow shock between the winds of the pulsar and companion,where the shock shapes the eclipsing medium but had not been described in detail in previous works.The model is further used to explain the variations of the flux density and dispersion measure of three BW pulsars(i.e.,PSR B1957+20,J2055+3829,and J2051-0827)detected by the Five-hundred-meter Aperture Spherical radio Telescope.Consequently,we constrained the parameters of the three BW systems such as the inclination angles and true anomalies of the observer as well as the mass-loss rates and wind velocity of the companion stars.With the help of these constraints,it is expected that magnetic fields of companion stars and even masses of pulsars could further be determined as some extra observation can be achieved in the future.