Because of the development of many-body theories of nuclear matter, the long-standing, open problem of the equation of state (EOS) of dense matter may be understood in the near future through the confrontation of theoretical calculations with laboratory measurements of nuclear properties \& reactions and increasingly accurate observations in astronomy. In this review, we focus on the following six aspects: 1) providing a survey of the quark mean-field (QMF) model, which consistently describes a nucleon and many-body nucleonic system from a quark potential; 2) applying QMF to both nuclear matter and neutron stars; 3) extending QMF formalism to the description of hypernuclei and hyperon matter, as well as hyperon stars; 4) exploring the hadron-quark phase transition and hybrid stars by combining the QMF model with the quark matter model characterized by the sound speed; 5) constraining interquark interactions through both the gravitational wave signals and electromagnetic signals of binary merger event GW170817; and 6) discussing further opportunities to study dense matter EOS from compact objects, such as neutron star cooling and pulsar glitches.

The unknown equation of state (EoS) of neutron stars (NSs) is puzzling because of rich non-perturbative effects of strong interaction there. A method to constrain the EoS by using the detected X-ray plateaus of gamma-ray bursts (GRBs) is proposed in this paper. Observations show some GRB X-ray plateaus may be powered by strongly magnetized millisecond NSs. The properties of these NSs should then satisfy: (i) the spin-down luminosity of these NSs should be brighter than the observed luminosity of the X-ray plateaus; (ii) the total rotational energy of these NSs should be larger than the total energy of the X-ray plateaus. Through the case study of GRB 170714A, the moment of inertia of NSs is constrained as $I>1.0\times 10^{45}\left ( \frac{P_{\rm cri}}{1\;\rm ms} \right )^{2} \;\rm g\cdot cm^{2}$, where $P_{\rm cri}$ is the critical rotational period that an NS can achieve. The constraint of the radii of NSs according to GRB 080607 is shown in Table 1.

We revisit the constraint on the maximum mass of cold spherical neutron stars coming from the observational results of GW170817. We develop a new framework for the analysis by employing both energy and angular momentum conservation laws as well as solid results of latest numerical-relativity simulations and of neutron stars in equilibrium. The new analysis shows that the maximum mass of cold spherical neutron stars can be only weakly constrained as $M_{\rm max} \alt 2.3M_\odot$. Our present result illustrates that the merger remnant neutron star at the onset of collapse to a black hole is not necessarily rapidly rotating and shows that we have to take into account the angular momentum conservation law to impose the constraint on the maximum mass of neutron stars.

Rapidly and differentially rotating compact stars are believed to be formed in binary neutron star merger events, according to both numerical simulations and the multi-messenger observation of GW170817. The lifetime and evolution of such a differentially rotating star, is tightly related to the observations in the post-merger phase. Various studies on the maximum mass of differentially rotating neutron stars have been done in the past, most of which assume the so-called $j$-const law as the rotation profile inside the star and consider only neutron star equations of state. In this paper, we extend the studies to strange star models, as well as to a new rotation profile model. Significant differences are found between differentially rotating strange stars and neutron stars, with both differential rotation laws. A moderate differential rotation rate for neutron stars is found to be too large for strange stars, resulting in a rapid drop in the maximum mass as the differential rotation degree is increased further from $\hat{A}\sim2.0$, where $\hat{A}$ is a parameter characterizing the differential rotation rate for $j$-const law. As a result the maximum mass of a differentially rotating self-bound star drops below the uniformly rotating mass shedding limit for a reasonable degree of differential rotation. The continuous transition to the toroidal sequence is also found to happen at a much smaller differential rotation rate and angular momentum than for neutron stars. In spite of those differences, $\hat{A}$-insensitive relation between the maximum mass for a given angular momentum is still found to hold, even for the new differential rotation law. Astrophysical consequences of these differences and how to distinguish between strange star and neutron star models with future observations are also discussed.

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.

Matter state inside neutron stars is an exciting problem in astrophysics, nuclear physics and particle physics. The equation of state (EOS) of neutron stars plays a crucial role in the present multimessenger astronomy, especially after the event of GW170817. We propose a new neutron star EOS "QMF18" from the quark level, which describes well robust observational constraints from free-space nucleon, nuclear matter saturation, heavy pulsar measurements and the tidal deformability of the very recent GW170817 observation. For this purpose, we employ the quark-mean-field (QMF) model, allowing one to tune the density dependence of the symmetry energy and study effectively its correlations with the Love number and the tidal deformability. We provide tabulated data for the new EOS and compare it with other recent EOSs from various many-body frameworks.

The LIGO/VIRGO detection of the gravitational waves from a binary merger system, GW170817, has put a clean and strong constraint on the tidal deformability of the merging objects. From this constraint, deep insights can be obtained in compact star equation of states, which has been one of the most puzzling problems for nuclear physicists and astrophysicists. Employing one of the most widely-used quark star EOS model, we characterize the star properties by the strange quark mass ($m_s$), an effective bag constant ($B_{\rm eff}$), the perturbative QCD correction ($a_4$), as well as the gap parameter ($\Delta$) when considering quark pairing, and investigate the dependences of the tidal deformablity on them. We find that the tidal deformability is dominated by $B_{\rm eff}$, and insensitive to $m_s$, $a_4$. We discuss the correlation between the tidal deformability and the maximum mass ($M_\mathrm{TOV}$) of a static quark star, which allows the model possibility to rule out the existence of quark stars with future gravitational wave observations and mass measurements. The current tidal deformability measurement implies $M_\mathrm{TOV} \le2.18\,M_\odot$ ($2.32\,M_\odot$ when pairing is considered) for quark stars. Combining with two-solar-mass pulsar observations, we also make constraints on the poorly known gap parameter $\Delta$ for color-flavor-locked quark matter.

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.

Quasi-equilibrium models of uniformly rotating axisymmetric and triaxial quark stars are computed in general relativistic gravity scenario. The Isenberg-Wilson-Mathews (IWM) formulation is employed and the Compact Object CALculator (COCAL) code is extended to treat rotating stars with finite surface density and new equations of state (EOSs). Besides the MIT bag model for quark matter which is composed of de-confined quarks, we examine a new EOS proposed by Lai and Xu that is based on quark clustering and results in a stiff EOS that can support masses up to $3.3M_\odot$ in the case we considered. We perform convergence tests for our new code to evaluate the effect of finite surface density in the accuracy of our solutions and construct sequences of solutions for both small and high compactness. The onset of secular instability due to viscous dissipation is identified and possible implications are discussed. An estimate of the gravitational wave amplitude and luminosity based on quadrupole formulas is presented and comparison with neutron stars is discussed.

TThe glitch of anomalous X-ray pulsars \& soft gamma repeaters (AXP/SGRs) usually accompanied with detectable energy releases manifesting as X-ray bursts or outbursts, while the glitch of some pulsars like Vela release negligible energy. We find that these two types of glitches can naturally correspond to two types of starquake of solid strange stars. By applying the EoS of quark cluster star and some realistic pulsar parameters, we can reproduce consistent results compared with previous constraints and observations.

Glitch (sudden spinup) is a common phenomenon in pulsar observations. However, the physical mechanism of glitch is still a matter of debate because it depends on the puzzle of pulsar's inner structure, i.e., the equation of state of dense matter. Some pulsars (e.g., Vela-like) show large glitches ({\Delta}{\nu}/{\nu}~10^-6) but release negligible energy, whereas the large glitches of AXPs/SGRs (anomalous X-ray pulsars/soft gamma repeaters) are usually (but not always) accompanied with detectable energy releases manifesting as X-ray bursts or outbursts. We try to understand this aspect of glitches in a starquake model of solid quark stars. There are actually two kinds of glitches in this scenario: bulk-invariable (Type I) and bulk-variable (Type II) ones. The total stellar volume changes (and then energy releases) significantly for the latter but not for the former. Therefore, glitches accompanied with X-ray bursts (e.g., that of AXP/SGRs) could originate from Type II starquakes induced probably by accretion, while the others without evident energy release (e.g., that of Vela pulsar) would be the result of Type I starquakes due to, simply, a change of stellar ellipticity.

Glitch (sudden spinup) is a common phenomenon in pulsar observations. However, the physical mechanism of glitch is still a matter of debate because it depends on the puzzle of pulsar's inner structure, i.e., the equation of state of dense matter. Some pulsars (e.g., Vela-like) show large glitches ({\Delta}{\nu}/{\nu}~10^-6) but release negligible energy, whereas the large glitches of AXPs/SGRs (anomalous X-ray pulsars/soft gamma repeaters) are usually (but not always) accompanied with detectable energy releases manifesting as X-ray bursts or outbursts. We try to understand this aspect of glitches in a starquake model of solid quark stars. There are actually two kinds of glitches in this scenario: bulk-invariable (Type I) and bulk-variable (Type II) ones. The total stellar volume changes (and then energy releases) significantly for the latter but not for the former. Therefore, glitches accompanied with X-ray bursts (e.g., that of AXP/SGRs) could originate from Type II starquakes induced probably by accretion, while the others without evident energy release (e.g., that of Vela pulsar) would be the result of Type I starquakes due to, simply, a change of stellar ellipticity.

In this paper we construct an expanding phase with phantom matter, in which the scale factor expands very slowly but the Hubble parameter increases gradually, and assume that this expanding phase could be matched to our late observational cosmology by the proper mechanism. We obtain the nearly scale-invariant spectrum of adiabatic fluctuations in this scenario, different from the simplest inflation and usual ekpyrotic/cyclic scenario, the tilt of nearly scale-invariant spectrum in this scenario is blue. Although there exists an uncertainty surrounding the way in which the perturbations propagate through the transition in our scenario, which is dependent on the detail of possible "bounce" physics, compared with inflation and ekpyrotic/cyclic scenario, our work may provide another feasible cosmological scenario generating the nearly scale-invariant perturbation spectrum.