The fundamental strong interaction determines the nature of pulsar-like compact stars which are essentially in the form of bulk strong matter. From an observational point of view, it is proposed that bulk strong matter could be composed of strangeons, i.e. quark-clusters with there-light-flavor symmetry of quarks, and therefore pulsar-like compact objects could actually be strangeon stars. The equation of state (EOS) of strangeon stars is described in a Lennard-Jones model for the purpose of constraining the EOS by both the tidal deformability $\Lambda$ of GW170817 and $M_{\rm TOV}$. It is found that the allowed parameter space is quite large as most of the Lennard-Jones EOS models satisfy the tidal deformability constraint by GW170817. The future GW detections for smaller values of $\Lambda$ and mass measurement for larger values of $M_{\rm TOV}$ will help a better constraint on the strangeon star model.

In this White Paper we present the potential of the Enhanced X-ray Timing and Polarimetry (eXTP) mission for determining the nature of dense matter; neutron star cores host an extreme density regime which cannot be replicated in a terrestrial laboratory. The tightest statistical constraints on the dense matter equation of state will come from pulse profile modelling of accretion-powered pulsars, burst oscillation sources, and rotation-powered pulsars. Additional constraints will derive from spin measurements, burst spectra, and properties of the accretion flows in the vicinity of the neutron star. Under development by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

X-ray-dim isolated neutron stars (XDINSs), also known as the Magnificent Seven, exhibits a Planck-like soft X-ray spectrum. In the optical/ultraviolet(UV) band, there is an excess of radiation compared to the extrapolation from the X-ray spectrum. A model of bremsstrahlung emission from a nonuniform plasma atmosphere is proposed in the regime of a strangeon star to explain the optical/UV excess and its spectral deviation. The nonuniform atmosphere is formed by the accretion of interstellar matter which moves along the magnetic field lines to polar caps, and these particles may spread out of the pole regions. The modeled electron temperatures are $\sim100-200$ eV with radiation radii $R_{\rm opt}^{\infty}\sim5-14$ km. The spectra of five sources (RX J0720.4--3125, RX J0806.4--4123, RX J1308.6+2127, RX J1605.3+3249, RX J1856.5--3754) from optical/UV to X-ray bands could be well fitted by the radiative model, and exhibit gaussian absorption lines at $\sim 100-500$ eV as would be expected. Furthermore, the surroundings (i.e., fallback disks or dusty belts) of XDINSs could be tested by future observations with the Large Optical Telescope (LOT).

Pulsar-like compact stars provide us a unique laboratory to explore properties of dense matter at supra-nuclear densities. One of the models for pulsar-like stars is that they are totally composed of "strangeons", and in this paper we studied the pulsar glitches in a strangeon star model. Strangeon stars would be solidified during cooling, and the solid stars would be natural to have glitches as the result of starquakes. Based on the starquake model established before, we proposed that when the starquake occurs, the inner motion of the star which changes the moment of inertia and has impact on the glitch sizes, is divided into plastic flow and elastic motion. The plastic flow which is induced in the fractured part of the outer layer, would move tangentially to redistribute the matter of the star and would be hard to recover. The elastic motion, on the other hand, changes its shape and would recover significantly. Under this scenario, we could understand the behaviors of glitches without significant energy releasing, including the Crab and the Vela pulsars, in an uniform model. We derive the recovery coefficient as a function of glitch size, as well as the time interval between two successive glitches as the function of the released stress. Our results show consistency with observational data under reasonable ranges of parameters. The implications on the oblateness of the Crab and the Vela pulsars are discussed.

The state of supranuclear matter in compact star remains puzzling, and it is argued that pulsars could be strangeon stars. What 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 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 surface by 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 binary 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 (GRB), is calculated. It is found that, after a promote burst, an X-ray plateau could follow in a timescale of $10^{2-3}$ s. Certainly, the results could be tested also by further observational synergies between gravitational wave detectors (e.g., aLIGO) and X-ray telescopes (e.g., Chinese HXMT and eXTP), and especially if the detected gravitational wave form is checked by peculiar equation of state provided by the numerical relativistical simulation.

The state equation of strangeon matter is very stiff due to the non-relativistic nature of and the repulsing interaction between the particles, and pulsar masses as high as $\sim 3M_\odot$ would be expected. However, an adiabatic sound speed, $c_s=\sqrt{\partial P/\partial \rho}$, is usually superluminal for strangeon matter, and dynamic response of strangeon star (e.g., binary merger) could not be tractable in calculation. We examine signal propagation in strangeon matter, and calculate the propagation speed, $c_{\rm signal}$, in reality. It is found that as the causality condition is satisfied, i.e., $c_{\rm signal}<c$, and the signal speed as a function of stellar radius is presented.

The neutrino burst detected during supernova SN1987A is explained in a strangeon star model, in which it is proposed that a pulsar-like compact object is composed of strangeons (strangeon: an abbreviation of "strange nucleon"). A nascent strangeon star's initial internal energy is calculated, with the inclusion of pion excitation (energy around 10^53 erg, comparable to the gravitational binding energy of a collapsed core). A liquid-solid phase transition at temperature ~ 1-2 MeV may occur only a few ten-seconds after core-collapse, and the thermal evolution of strangeon star is then modeled. It is found that the neutrino burst observed from SN 1987A could be re-produced in such a cooling model.

The nature of pulsar-like compact stars is essentially a central question of the fundamental strong interaction (explained in quantum chromo-dynamics) at low energy scale, the solution of which still remains a challenge though tremendous efforts have been tried. This kind of compact objects could actually be strange quark stars if strange quark matter in bulk may constitute the true ground state of the strong-interaction matter rather than Fe^56 (the so-called Witten's conjecture). From astrophysical points of view, however, it is proposed that strange cluster matter could be absolutely stable and thus those compact stars could be strange cluster stars in fact. This proposal could be regarded as a general Witten's conjecture: strange matter in bulk could be absolutely stable, in which quarks are either free (for strange quark matter) or localized (for strange cluster matter). Strange cluster with three-light-flavor symmetry is renamed strangeon, being coined by combining "strange nucleon" for the sake of simplicity. A strangeon star can then be thought as a 3-flavored gigantic nucleus, and strangeons are its constituent as an analogy of nucleons which are the constituent of a normal (micro) nucleus. The observational consequences of strangeon stars show that different manifestations of pulsar-like compact stars could be understood in the regime of strangeon stars, and we are expecting more evidence for strangeon star by advanced facilities (e.g., FAST, SKA, and eXTP).

Aims. We aim to probe the internal structure and investigate more detailed information of the $\delta$ Scuti star HD 50844 with asteroseismology. Methods. We analyse the observed frequencies of the $\delta$ Scuti star HD 50844 obtained by Balona (2014), and search for possible multiplets based on the rotational splitting law of g-mode. We tried to disentangle the frequency spectra of HD 50844 by means of the rotational splitting only. We then compare them with theoretical pulsation modes, which correspond to stellar evolutionary models with various sets of initial metallicity and stellar mass, to find the best-fitting model. Results. There are three multiplets including two complete triplets and one incomplete quintuplet, in which mode identifications for spherical harmonic degree $l$ and azimuthal number $m$ are unique. The corresponding rotational period of HD 50844 is found to be 2.44$^{+0.13}_{-0.08}$ days. The physical parameters of HD 50844 are well limited in a small region by three modes identified as nonradial ones ($f_{11}$, $f_{22}$, and $f_{29}$) and by the fundamental radial mode ($f_{4}$). Our results show that the three nonradial modes ($f_{11}$, $f_{22}$, and $f_{29}$) are all mixed modes, which mainly represent the property of the helium core. The fundamental radial mode ($f_{4}$) mainly represents the property of the stellar envelope. In order to fit these four pulsation modes, both the helium core and the stellar envelope must be matched to the actual structure of HD 50844. Finally, the mass of the helium core of HD 50844 is estimated to be 0.173 $\pm$ 0.004 $M_{\odot}$ for the first time. The physical parameters of HD 50844 are determined to be $M=$ 1.81 $\pm$ 0.01 $M_{\odot}$, $Z=$ 0.008 $\pm$ 0.001. $T_{\rm eff}=$ 7508 $\pm$ 125 K, log$g=$ 3.658 $\pm$ 0.004, $R=$ 3.300 $\pm$ 0.023 $R_{\odot}$, $L=$ 30.98 $\pm$ 2.39 $L_{\odot}$.

Pulsar-like compact stars usually have strong magnetic fields, with the strength from $\sim 10^8$ to $\sim 10^{12}$ Gauss on surface. How such strong magnetic fields can be generated and maintained is still an unsolved problem, which is, in principle, related to the interior structure of compact stars, i.e., the equation of state of cold matter at supra-nuclear density. In this paper we are trying to solve the problem in the regime of solid quark-cluster stars. Inside quark-cluster stars, the extremely low ratio of number density of electrons to that of baryons $n_e/n_b$ and the screening effect from quark-clusters could reduce the long-range Coulomb interaction between electrons to short-range interaction. In this case, the Stoner's model could apply, and we find that the condition for ferromagnetism is consistent with that for validity of Stoner's model. Under the screened Coulomb repulsion, the electrons inside the stars could spontaneously magnetized and become ferromagnetic, and hence would contribute non-zero net magnetic momentum to the whole star. We conclude that, for most cases in solid quark-cluster stars, the amount of net magnetic momentum, which is proportional to the amount of unbalanced spins $\xi=(n_+-n_-)/n_e$ and depends on the number density of electrons $n_e=n_++n_-$, could be significant with non-zero $\xi$. The net magnetic moments of electron system in solid quark-cluster stars could be large enough to induce the observed magnetic fields for pulsars with $B\sim 10^{11}$ to $\sim 10^{13}$ Gauss.

Pulsar glitches, i.e. the sudden spin-ups of pulsars, have been detected for most pulsars that we known. The mechanism giving rise to this kind of phenomenon is uncertain, although a large data set has been built. In the framework of star-quake model, based on~\cite{Baym1971}, the glitch-sizes (the relative increases of spin-frequencies during glitches) $\Delta \Omega/\Omega$ depend on the released energies during glitches, with less released energies corresponding to smaller glitch sizes. On the other hand, as one of dark matter candidates, our Galaxy might be filled with the so called strange nuggets (SNs) which are the relics from the early Universe. In this case the collisions between pulsars and SNs is inevitable, and these collisions would lead to glitches when enough elastic energy has been accumulated during the spin-down process. The SNs-triggered glitches could release less energy, because the accumulated elastic energy would be less than that in the scenario of glitches without SNs. Therefore, if a pulsar is hit frequently by SNs, it would tend to have more small size glitches, whose values of $\Delta \Omega/\Omega$ are smaller than that in standard star-quake model (with larger amounts of released energy). Based on the assumption that in our Galaxy the distribution of SNs is similar to that of dark matter, as well as on the glitch data in ATNF Pulsar Catalogue and Jodrell Bank glitch table, we find that in our Galaxy the incidences of small size glitches exhibit tendencies consistent with the collision rates of pulsars and strange nuggets. Further test of this scenario is expected by detecting more small glitches (e.g., by the Square Kilometre Array).

In this paper we are trying to solve the problem of the origin of strong magnetic fields in the framework of solid quark-cluster stars. We propose that, under the Coulomb repulsion, the electrons inside the stars could spontaneously magnetized and become ferromagnetic, and hence would contribute non-zero net magnetic momentum to the whole star. The leading order approximation in our simple model shows that the magnetism of solid quark-cluster matter could be possible. For most cases in solid quark-cluster stars, the amount of net magnetic momentum, could be significant, and the net magnetic moments of electron system in solid quark-cluster stars could be large enough to induce the observed magnetic fields for pulsars with $B\sim 10^{11}$ to $\sim 10^{13}$ Gauss. This mechanism of generating magnetic field is not valid for so-called magnetars.

The state of super-dense matter is essential for us to understand the nature of pulsars, but the non- perturbative quantum chromodynamics (QCD) makes it very difficult for direct calculations of the state of cold matter at realistic baryon number densities inside compact stars. Nevertheless, from an observational point of view, it is conjectured that pulsars could be made up of quark clusters since the strong coupling between quarks might render quarks grouped in clusters. We are trying an effort to find an equation of state of condensed quark-cluster matter in a phenomenological way. Supposing that the quark-clusters could be analogized to inert gases, we apply here the corresponding-state approach to derive the equation of state of quark-cluster matter, as was similarly demonstrated for nuclear and neutron-star matter in 1970s. According to the calculations presented, the quark-cluster stars, which are composed of quark-cluster matter, could then have high maximum mass that is consistent with observations and, in turn, further observations of pulsar mass would also put constraints to the properties of quark-cluster matter. Moreover, the melting heat during solid-liquid phase conversion and the related astrophysical consequences are also briefly discussed.

The study of dense matter at ultra-high density has a very long history, which is meaningful for us to understand not only cosmic events in extreme circumstances but also fundamental laws of physics. It is well known that the state of cold matter at supra-nuclear density depends on the non-perturbative nature of quantum chromo-dynamics (QCD) and is essential for modeling pulsars. A so-called H-cluster matter is proposed in this paper as the nature of dense matter in reality. In compact stars at only a few nuclear densities but low temperature, quarks could be interacting strongly with each other there. That might render quarks grouped in clusters, although the hypothetical quark-clusters in cold dense matter has not been confirmed due to the lack of both theoretical and experimental evidence. Motivated by recent lattice QCD simulations of the H-dibaryons (with structure uuddss), we are therefore considering here a possible kind of quark-clusters, H-clusters, that could emerge inside compact stars during their initial cooling, as the dominant components inside (the degree of freedom could then be H-clusters there).Taking into account the in-medium stiffening effect, we find that at baryon densities of compact stars $H$-cluster matter could be more stable than nuclear matter. We also find that for the H-cluster matter with lattice structure, the equation of state could be so stiff that it would seem to be "superluminal" in most dense region. However, the real sound speed for H-cluster matter is in fact hard to calculate, so at this stage we do not put constraints on our model from the usual requirment of causality.

Based on the turbulent convection model (TCM), we investigate chemical mixing in the bottom overshooting region of the convective envelope of intermediate-mass stars, focusing on its influence on the formation and extension of blue loops in the Hertzsprung-Russell (HR) diagram. A diffusive mixing model is adopted during the Red Giant Branch (RGB) phase. The properties of the blue loop are changed by modification of the element profiles above the H-burning shell, which results from the incomplete mixing in the bottom overshooting region when stellar model evolves up along the RGB. Such modification of the element profiles will lead to an increase of the opacity in the region just above the H-burning shell and a decrease of the opacity in the outer homogeneous convection zone, which will result in a quick decrease of the H-shell nuclear luminosity $L_{H}$ when the stellar model evolves from the RGB tip to its bottom and, finally, a much weaker and smaller convection zone will be obtained in the stellar envelope. This helps to form a longer blue loop. The extension of the blue loop is very sensitive to the parameters $(C_{X}$ and $\alpha_{TCM})$ of the diffusive mixing model and of the TCM. The results mainly show that: 1) comparing to the results of the classical model with the mixing-length theory, the lengths of the obtained blue loops with different combinations of the values of $C_{X}$ and $\alpha_{TCM}$ are all increased and the length of the blue loop increases with the values of parameters $C_{X}$ and $\alpha_{TCM}$; 2) the diffusive mixing model can significantly extend the time of stellar models lingering on the blue side of the HR diagram, even though the length of the blue loop for the 7$M_{\odot}$ star has a less prominent difference between the classical and diffusive mixing model;

The nature of pulsar-like compact stars is still in controversy although the first pulsar was found more than 40 years ago. Generally speaking, conventional neutron stars and non-mainstream quark stars are two types of models to describe the inner structure of pulsars, with the former composed mainly of hadrons and the latter of a peculiar kind of matter whose state equation should be understood in the level of quarks rather than hadrons. To construct a more realistic model from both theoretical and observational points of view, we conjecture that pulsars could be "quark-cluster stars" which are composed of quark-clusters with almost equal numbers of up, down and strange quarks. Clustering quark matter could be the result of strong coupling between quarks inside realistic compact stars. The lightest quark clusters could be of H-dibaryons, while quark clusters could also be heavier with more quarks. Being essentially related to the non-perturbative quantum-chromo dynamics (QCD), the state of supra-nuclear condensed matter is really difficult to obtain strictly by only theoretical QCD-calculations, and we expect, nevertheless, that astrophysical observations could help us to have a final solution.

Based on the turbulent convection model (TCM) of Li & Yang (2007), we have studied the characteristics of turbulent convection in the envelopes of 2 and 5M stars at the RGB and AGB phases. The TCM has been applied successfully in the whole convective envelopes including the convective unstable zone and the overshooting regions. We find that the convective motions become stronger and stronger when the stellar models are located upper and upper along the Hayashi line. In the convective unstable zone, we find that the turbulent correlations are proportional to functions of a common factor $(\nabla-\nabla_{ad})\bar{T}$, which explains similar distributions of those correlations. For the TCM we find that if the obtained stellar structure of temperature is close to that of the MLT, the convective motion will be at a much larger velocity and thus more violent. However, if the turbulent velocity is adjusted to close to that of the MLT, the superadiabatic convection zone is much more extended inward and a lower effective temperature of the stellar model will be obtained. For the overshooting distance we find that the e-folding lengths of the turbulent kinetic energy k in both the top and bottom overshooting regions decrease as the stellar model is located up along the Hayashi line, but both the extents of the decrease are not obvious. And the overshooting distances of the turbulent correlation $\bar{u_{r}^{'}T^{'}}$ are almost the same for the different stellar models with a same set of TCM's parameters. ...

It is conventionally thought that the state equation of dense matter softens and thus cannot result in high maximum mass if pulsars are quark stars, and that a recently discovered $2M_\odot$ pulsar (PSR J1614-2230) may make quark stars to be unlikely. However, this standard point of view would be revisited and updated if quark clustering could occur in cold quark matter because of the strong coupling between quarks at realistic baryon densities of compact stars, and it is addressed that the state equation of clustering quark matter stiffs to support compact stars with maximum mass $M_{\rm max}>2M_\odot$. In this brief note, it is demonstrated that large parameter spaces are allowed for $M_{\rm max}>2M_\odot$ in a Lennard-Jones model of clustered quark matter, and the newly measured highest mass of PSR J1614-2230 would be meaningful to constrain the number of quarks inside single quark-cluster, to be $N_q<\sim 10^3$.

Strange quark nuggets (SQNs) could be the relics of the cosmological QCD phase transition, and they could very likely be the candidate of cold quark matter if survived the cooling of the later Universe, although the formation and evolution of these SQNs depend on the physical state of the hot QGP (quark-gluon plasma) phase and the state of cold quark matter. We reconsider the possibility of SQNs as cold dark matter, and find that the formation of black holes in primordial halos could be significantly different from the standard scenario. In a primordial halo, the collision between gas and SQNs could be frequent enough, and thus the viscosity acting on each SQN would decrease its angular momentum and make it to sink into the center of the halo, as well as heat the gas. The SQNs with baryon numbers less than $10^{35}$ could assemble in the center of the halo before the formation of primordial stars. A black hole could form by merger of these SQNs, and then its mass could quickly become about $10^3\ M_\odot$ or higher, by accreting the surrounding SQNs or gas. The black holes formed in this way could be the seeds for the supermassive black holes at redshift as high as $z\sim 6$.

It depends on the state of matter at supra-nuclear density to model pulsar's structure, which is unfortunately not certain due to the difficulties in physics. In cold quark matter at realistic baryon densities of compact stars (with an average value of $\sim 2-3\rho_0$), the interaction between quarks is so strong that they would condensate in position space to form quark-clusters. We argue that quarks in quark stars are grouped in clusters, then we apply two phenomenological models for quark stars, the polytropic model and Lennard-Jones model. Both of the two models have stiffer EoS, and larger maximum mass for quark stars (larger than 2 $M_\odot$). The gravitational energy releases during the AIQ process could explain the observed energy of three supergiant flares from soft gamma-ray repeaters ($\sim 10^{47}$ ergs).

Quark clustering could occur in cold quark matter because of the strong coupling between quarks at realistic baryon densities of compact stars. Although one may still not be able to calculate this conjectured matter from first principles, the inter-cluster interaction might be analogized to the interaction between inert molecules. Cold quark matter would then crystallize in a solid state if the inter-cluster potential is deep enough to trap the clusters in the wells. We apply the Lennard-Jones potential to describe the inter-cluster potential, and derive the equations of state, which are stiffer than that derived in conventional models (e.g., MIT bag model). If quark stars are composed of Lennard-Jones matter, they could have high maximum masses ($>2M_{\odot}$) as well as very low masses ($<10^{-3}M_{\odot}$). These features could be tested by observations.

A polytropic quark star model is suggested in order to establish a general framework in which theoretical quark star models could be tested by observations. The key difference between polytropic quark stars and the polytropic model studied previously for normal (i.e., non-quarkian) stars is related to two issues: (i) a constant term representing the contribution of vacuum energy may be added in the energy density and the pressure for a quark star, but not for a normal star; (ii) the quark star models with non-vanishing density at the stellar surface are not avoidable due to the strong interaction between quarks. The first one implies that the vacuum inside a quark star is different from that outside, while the second one is relevant to the effect of color confinement. The polytropic equations of state are stiffer than that derived in conventional realistic models (e.g., the bag model) for quark matter, and pulsar-like stars calculated with a polytropic equation of state could then have high maximum masses (> 2 M_sun). Quark stars can also be very low massive, and be still gravitationally stable even if the polytropic index, n, is greater than 3. All these would result in different mass-radius relations, which could be tested by observations. In addition, substantial strain energy would develop in a solid quark star during its accretion/spindown phase, and could be high enough to take a star-quake. The energy released during star-quakes could be as high as ~ 10^{47} ergs if the tangential pressure is ~ 10^{-6} higher than the radial one.