We study the chiral condensates in neutron star matter from nuclear to quark matter domain. We describe nuclear matter with a parity doublet model (PDM), quark matter with the Nambu--Jona-Lasino (NJL) model, and a matter at the intermediate density b
y interpolating nuclear and quark matter equations of state. The model parameters are constrained by nuclear physics and neutron star observations. Various condensates in the interpolated domain are estimated from the chemical potential dependence of the condensates at the boundaries of the interpolation. The use of the PDM with substantial chiral invariant mass ($m_0 gtrsim 500$ MeV, which is favored by the neutron star observations) predicts the mild chiral restoration, and the significant chiral condensate remains to baryon density $n_B sim 2-3n_0$ ($n_0simeq 0.16,{rm fm}^{-3}$: nuclear saturation density), smoothly approaching the NJL predictions for the color-flavor-locked phase at $n_B gtrsim 5n_0$. The same method is applied to estimate diquark condensates, number densities of up-, down- and strange-quarks, and the lepton fraction. In our descriptions the chiral restoration in the interpolated domain proceeds with two conceptually distinct chiral restoration effects; the first is associated with the positive scalar density in a nucleon, relevant in dilute regime, and the other primarily arises from the modification of the quark Dirac sea, which is triggered by the growth of the quark Fermi sea. We discuss several qualitative conjectures to interpolate the microphysics in nuclear and quark matter.
We construct an equation of state (EOS) for neutron stars by interpolating hadronic EOS at low density and quark EOS at high density. A hadronic model based on the parity doublet structure is used for hadronic matter and a quark model of Nambu--Jona-
Lasinio type is for quark matter. We assume crossover between hadronic matter and quark matter in the the color-flavor locked phase. The nucleon mass of the parity doublet model has a mass associated with the chiral symmetry breaking, and a chiral invariant mass $m_0$ which is insensitive to the chiral condensate. The value of $m_0$ affects the nuclear EOSs at low density, and has strong correlations with the radii of neutron stars. Using the constraint to the radius obtained by LIGO-Virgo and NICER, we find that $m_0$ is restricted as $600,mathrm{MeV}lesssim m_0 lesssim 900,mathrm{MeV}$.
