Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Deriving Neutron Star Equation of State from AdS/QCD

Wei Li wei_li_phys@163.com Department of Physics and Institute for Quantum Science and Technology, Shanghai University, 99 Shangda Road, Shanghai 200444, China    Jing-Yi Wu wujingyi222@mails.ucas.ac.cn School of Astronomy and Space Science, University of Chinese Academy of Sciences (UCAS), Beijing 100049, China    Kilar Zhang kilar@shu.edu.cn, Corresponding Author Department of Physics and Institute for Quantum Science and Technology, Shanghai University, 99 Shangda Road, Shanghai 200444, China Shanghai Key Lab for Astrophysics, Shanghai 200234, China Shanghai Key Laboratory of High Temperature Superconductors, Shanghai 200444, China
Abstract

Neutron stars are among the main targets for gravitational wave observatories, however, their equation of state is still not well established. Mainly phenomenological models with many parameters are widely used by far, while theoretical models are not so practical. In arXiv:1902.08477, a theoretical equation of state with only one parameter is derived from Sakai-Sugimoto model, as an application of AdS/QCD, where pointlike instanton case is taken into consideration. When the tidal deformability constraint from gravitational wave event is satisfied, the maximum mass is about 1.7 solar masses. Now we upgrade this model to instanton gas, with one more variable, the instanton width. This is not naively a free parameter, but a function of the chemical potential. Thus we end up with a more complicated and accurate model, but still with only one adjustable parameter. In this case, we find the maximum mass becomes 1.85 solar masses. This is an encouraging and exiting result, as a theoretically derived model.

I I. Introduction

The observations of neutron stars (NS) based on electromagnetic (EM) waves Margalit:2019dpi have a long history, yet the equation of state (EoS) of NS has remained as a puzzle. The precise derivation of EoS requires much higher density than neutron saturation level, the region where no perturbative methods work and almost impossible to settle. Nonetheless, there are still many phenomenological or approximate models. With EoS in hand, one can find the Mass-Radius relation and also the tidal Love number (TLN), by solving the Tolman-Oppenheimer-Volkoff (TOV) equations together Hinderer_2009 . That means, on the contrary, if we have the information of Mass-Radius or TLN, we can constrain the EoS, thus offering an alternative method to probe the high density region in nuclear physics Lattimer:2012nd .

Then, what is the situation on the observation side? We know it is easy to read off the masses of NS to high accuracy, especially for (binary) pulsars 2009Physics . While for the radius, realities are much less optimistic. Only through X-ray bursts Galloway:2007dn or thermal dynamics, could we find some limited information on the radius, with an error bar of several kilometers for a standard 11km NS 2020Sci…370.1450D . This could barely serve as a constraint, so there are more than one hundred common EoS models: no one violates the requirements.

Since the advent of Gravitational Waves observatories, especially after the first NS GW event, we enter the era of multi-messenger astronomy 2017Multi . The advance of GW data is not only its high efficiency in finding black holes (BH) and NS, but also an extra parameter that beyond the normal EM observation ability: the TLN. This TLN reflects the deformability of a compact star under an external quadratic force. In a sense, it shows the compactness of a star. Its theoretical value can actually be deduced from a second order perturbation of Einstein equation, along with the EoS. As a result, if we know the value of TLN from observation, inversely we could constrain the EoS. In fact, the TLN from GW offers very strong constraints, and even merely by the First Binary NS event GW170817 LIGOScientific:2017vwq ; LIGOScientific:2018hze . Half (including one marginal) of the 7 major EoS candidates are ruled out. So far the story is so good, yet no further strong TLN constraints are obtained, since TLN is 5 Post-Newtonian and difficult to obtain for distant events. For another NS events GW190425 LIGOScientific:2020aai , one have TLN much higher with only one upper limit, while for others, we don’t have any data. The window for EoS is still big enough for half of the original candidates to survive, say some 50 kinds.

Furthermore, people are not satisfied with pure phenomenological models, and the explore on theoretical derivations continues. Starting from the MIT bag model Rezaei:2024ydg , many efforts have been made. In Zhang:2019tqd , a new perspective is tried: EoS can be extracted from Sakai-Sugimoto model (SS model) Sakai:2005yt ; Sakai:2004cn , an application of the AdS/QCD dualitywitten1998anti ; Polchinski:2000uf , originally from superstring theory. The seeming surprising relation is in a sense very natural, since the initial goal of string theory is to explain strong interactions. The duality offers a way to deal with the strong coupled region in nuclear physics by the help of brane constructions. The main idea is to use D4/D8 branes, and the strings connecting different branes. In Zhang:2019tqd a simplified case is taken into consideration: the pointlike instanton case Li_2015 ; Bergman:2007wp . The EoS derived has only one adjusting parameter: the maximum separation distance between the branes. This is a great advantage compared to other models with numerous coefficients, that you do not induce any parameter by hand. By generating the M-R relation and TLN-M relations, it turns out that if the TLN constraint is satisfied, the maximal mass supported by this EoS is 1.7M1.7subscript𝑀direct-product1.7M_{\odot}1.7 italic_M start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT (solar masses), unfortunately not enough compared with the known observation of NS with over 2M2subscript𝑀direct-product2M_{\odot}2 italic_M start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT.

In this paper, we further upgrade the derivation to the more complicated instanton gas case Li_2015 ; Ghoroku:2012am , with one more parameter, the instanton width, which was fixed to zero in the pointlike instanton case. This introduces much more difficulties in calculation, since the once single variation now becomes a triple-variation, which needs more effort to solve.

A naive guess would be that this new parameter helps to enlarge the parameter space, thus the maximum mass could be lifted by combining the two parameters properly. However, it turns out that the instanton width is not a free parameter, but a function of the chemical potential. Thus, for instanton gas case, though the model is more accurate and complicated, we still end up with only one adjusting parameter, i.e. the brane separation distance. Then, by similar process, we generate the MR and TLN-M relation by applying this EoS, and find that when satisfying the TLN constraints, the maximum mass is 1.85M1.85subscript𝑀direct-product1.85M_{\odot}1.85 italic_M start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT. This is higher than the pointlike case discussed before, yet still below the criteria of 2M2subscript𝑀direct-product2M_{\odot}2 italic_M start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT from NS observation.

Nonetheless, this fact does cannot deny the power of this EoS. Considering that a real neutron star contains more than pure neutron, a multi layer or mixed model Zhang:2020dfi ; Zhang:2020pfh with different compositions should be more realistic, which could circumference the maximum mass problem.

View differently, although very inspiring, SS model still cannot grab the whole picture of NS EoS, as there are some minor flaws about this nice theory from the beginning, that need to think over more seriously. First,this duality is supposed to be exact only when Ncsubscript𝑁𝑐N_{c}italic_N start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT goes to infinity, yet in practice we say Nc=3subscript𝑁𝑐3N_{c}=3italic_N start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = 3 is large enough; in this instanton gas model we have second order phase transition, rather than the first order transition 111the homogeneous anstaz case does reproduce first order transition, but has no chiral restoration. If we want a more realistic model, we need improve the SS model itself, or construct a new holographic model from the beginning. We leave this to the future work.

The organization of the paper is as the following: After this introduction in section 1, we briefly review the main idea of SS model and introduce our constructions for EoS extraction in section 2. In section 3 the M-R and TLN-M relations are shown, and we apply the GW constraints. We conclude in Section 4. Some calculation details are collected in the appendix.

II II. Extracting Neutron Star Equation of State from Sakai-Sugimoto Model

II.1 II.1 Sakai-Sugimoto model

Our construction is based on the "top-down" model of holographic quantum chromodynamics originated from string theory: SS model. Let’s briefly collect the necessary basics, following the procedure in Li_2015 .

We are considering Nfsubscript𝑁𝑓N_{f}italic_N start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT D8- and Nfsubscript𝑁𝑓N_{f}italic_N start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT D8¯¯D8\overline{\rm D8}over¯ start_ARG D8 end_ARG-branes around Ncsubscript𝑁𝑐N_{c}italic_N start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT D4-branes background, connecting at a tip u=uc𝑢subscript𝑢𝑐u=u_{c}italic_u = italic_u start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT and maximally separated at u=𝑢u=\inftyitalic_u = ∞, with u𝑢uitalic_u the (dimensionless) holographic coordinate. Introducing the Kaluza-Klein Mass Mkksubscript𝑀𝑘𝑘M_{kk}italic_M start_POSTSUBSCRIPT italic_k italic_k end_POSTSUBSCRIPT Duff:1986hr , The asymptotic separation is \ellroman_ℓ, which is πMkk𝜋subscript𝑀𝑘𝑘\pi\over M_{kk}divide start_ARG italic_π end_ARG start_ARG italic_M start_POSTSUBSCRIPT italic_k italic_k end_POSTSUBSCRIPT end_ARG for confined geometry, and is a free parameter for deconfined case. We deal with the deconfined case in this letter, and the confined case can be obtained as a special case.

x4subscript𝑥4x_{4}italic_x start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT represents the coordinate of the compactified dimension, where the radius is 1Mkk1subscript𝑀𝑘𝑘1\over M_{kk}divide start_ARG 1 end_ARG start_ARG italic_M start_POSTSUBSCRIPT italic_k italic_k end_POSTSUBSCRIPT end_ARG, thus its periodicity is given as x4=x4+2πMkksubscript𝑥4subscript𝑥42𝜋subscript𝑀𝑘𝑘x_{4}=x_{4}+{2\pi\over M_{kk}}italic_x start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT = italic_x start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT + divide start_ARG 2 italic_π end_ARG start_ARG italic_M start_POSTSUBSCRIPT italic_k italic_k end_POSTSUBSCRIPT end_ARG. Notice that we are actually using the dimensionless counter parts as lower case terms Li_2015 , and all the equations should be understood in their dimensional form, like X4=X4+2πMkksubscript𝑋4subscript𝑋42𝜋subscript𝑀𝑘𝑘X_{4}=X_{4}+{2\pi\over M_{kk}}italic_X start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT = italic_X start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT + divide start_ARG 2 italic_π end_ARG start_ARG italic_M start_POSTSUBSCRIPT italic_k italic_k end_POSTSUBSCRIPT end_ARG. This separation boundary condition can be mathematically described as

=2uc𝑑ux42superscriptsubscriptsubscript𝑢𝑐differential-d𝑢superscriptsubscript𝑥4\displaystyle\ell=2\int_{u_{c}}^{\infty}dux_{4}^{\prime}roman_ℓ = 2 ∫ start_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT italic_d italic_u italic_x start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT (1)

where the derivative mark is taken with respect to u𝑢uitalic_u.

The Abelian part of the gauge fields are denoted by a^μsubscript^𝑎𝜇\hat{a}_{\mu}over^ start_ARG italic_a end_ARG start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT, and its 00 component a^0subscript^𝑎0\hat{a}_{0}over^ start_ARG italic_a end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT gives the chemical potential μ𝜇\muitalic_μ by its boundary value, μ=a^0()𝜇subscript^𝑎0\mu=\hat{a}_{0}({\infty})italic_μ = over^ start_ARG italic_a end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( ∞ ).

By Considering the gauge field action consisting a Dirac-Born-Infeld (DBI) and a Chern-Simons (CS) contribution, S=SDBI+SCS𝑆subscript𝑆𝐷𝐵𝐼subscript𝑆𝐶𝑆S=S_{DBI}+S_{CS}italic_S = italic_S start_POSTSUBSCRIPT italic_D italic_B italic_I end_POSTSUBSCRIPT + italic_S start_POSTSUBSCRIPT italic_C italic_S end_POSTSUBSCRIPT, the free energy density ΩΩ\Omegaroman_Ω can be written as

Ω𝒩=uc𝑑uu5/2ζ[1+g2+(nIQ)2u5]μnI,Ω𝒩superscriptsubscriptsubscript𝑢𝑐differential-d𝑢superscript𝑢52𝜁delimited-[]1subscript𝑔2superscriptsubscript𝑛𝐼𝑄2superscript𝑢5𝜇subscript𝑛𝐼\displaystyle\frac{\Omega}{\cal N}=\int_{u_{c}}^{\infty}du\,u^{5/2}\zeta\left[% 1+g_{2}+\frac{(n_{I}Q)^{2}}{u^{5}}\right]-\mu n_{I}\,,divide start_ARG roman_Ω end_ARG start_ARG caligraphic_N end_ARG = ∫ start_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT italic_d italic_u italic_u start_POSTSUPERSCRIPT 5 / 2 end_POSTSUPERSCRIPT italic_ζ [ 1 + italic_g start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + divide start_ARG ( italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT italic_Q ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_u start_POSTSUPERSCRIPT 5 end_POSTSUPERSCRIPT end_ARG ] - italic_μ italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , (2)

with the abbreviation

ζ1+g1+u3x42a^021+g2=u11/2x4k𝜁1subscript𝑔1superscript𝑢3superscriptsubscript𝑥42superscriptsubscript^𝑎021subscript𝑔2superscript𝑢112superscriptsubscript𝑥4𝑘\displaystyle\zeta\equiv\sqrt{\frac{1+g_{1}+u^{3}x_{4}^{\prime 2}-\hat{a}_{0}^% {\prime 2}}{1+g_{2}}}=\frac{u^{11/2}x_{4}^{\prime}}{k}italic_ζ ≡ square-root start_ARG divide start_ARG 1 + italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_u start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_x start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ 2 end_POSTSUPERSCRIPT - over^ start_ARG italic_a end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ 2 end_POSTSUPERSCRIPT end_ARG start_ARG 1 + italic_g start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG end_ARG = divide start_ARG italic_u start_POSTSUPERSCRIPT 11 / 2 end_POSTSUPERSCRIPT italic_x start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG start_ARG italic_k end_ARG (3)
=u5/2a^0nIQ=1+g11+g2k2u8+(nIQ)2u5.absentsuperscript𝑢52superscriptsubscript^𝑎0subscript𝑛𝐼𝑄1subscript𝑔11subscript𝑔2superscript𝑘2superscript𝑢8superscriptsubscript𝑛𝐼𝑄2superscript𝑢5\displaystyle=\frac{u^{5/2}\hat{a}_{0}^{\prime}}{n_{I}Q}=\frac{\sqrt{1+g_{1}}}% {\sqrt{1+g_{2}-\frac{k^{2}}{u^{8}}+\frac{(n_{I}Q)^{2}}{u^{5}}}}\,.= divide start_ARG italic_u start_POSTSUPERSCRIPT 5 / 2 end_POSTSUPERSCRIPT over^ start_ARG italic_a end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG start_ARG italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT italic_Q end_ARG = divide start_ARG square-root start_ARG 1 + italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG end_ARG start_ARG square-root start_ARG 1 + italic_g start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - divide start_ARG italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_u start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT end_ARG + divide start_ARG ( italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT italic_Q ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_u start_POSTSUPERSCRIPT 5 end_POSTSUPERSCRIPT end_ARG end_ARG end_ARG .

Here nIsubscript𝑛𝐼n_{I}italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT is the baryon number density, k is an auxiliary integration constant introduced during taking the equations of motion for a^0subscript^𝑎0\hat{a}_{0}over^ start_ARG italic_a end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. The detailed definition of other symbols g1subscript𝑔1g_{1}italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT, g2subscript𝑔2g_{2}italic_g start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT and Q𝑄Qitalic_Q are listed in appendix A.

In comparison with the pointlike instanton case considered in Zhang:2019tqd , we now upgrade the method to instanton gas case, with an extra parameter, the instanton width ρ𝜌\rhoitalic_ρ. This makes the model more realistic, though also more complicated to solve.

The temperature is set to zero to make life easier, and it is also confirmed in Zhang:2019tqd that a finite temperature only has minor effect.

By taking variations of the free energy, we obtain the following equations Li_2015 :

μ𝜇\displaystyle\muitalic_μ =\displaystyle== uc𝑑u[u5/22(g1nIζ1+g2nIζ)+a^0Q],superscriptsubscriptsubscript𝑢𝑐differential-d𝑢delimited-[]superscript𝑢522subscript𝑔1subscript𝑛𝐼superscript𝜁1subscript𝑔2subscript𝑛𝐼𝜁superscriptsubscript^𝑎0𝑄\displaystyle\int_{u_{c}}^{\infty}du\,\left[\frac{u^{5/2}}{2}\left(\frac{% \partial g_{1}}{\partial n_{I}}\zeta^{-1}+\frac{\partial g_{2}}{\partial n_{I}% }\zeta\right)+\hat{a}_{0}^{\prime}Q\right]\,,∫ start_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT italic_d italic_u [ divide start_ARG italic_u start_POSTSUPERSCRIPT 5 / 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 end_ARG ( divide start_ARG ∂ italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG ∂ italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT end_ARG italic_ζ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT + divide start_ARG ∂ italic_g start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG ∂ italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT end_ARG italic_ζ ) + over^ start_ARG italic_a end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_Q ] , (4a)
00\displaystyle 0 =\displaystyle== uc𝑑u[u5/22(g1ρζ1+g2ρζ)+nIa^0Qρ],superscriptsubscriptsubscript𝑢𝑐differential-d𝑢delimited-[]superscript𝑢522subscript𝑔1𝜌superscript𝜁1subscript𝑔2𝜌𝜁subscript𝑛𝐼superscriptsubscript^𝑎0𝑄𝜌\displaystyle\int_{u_{c}}^{\infty}du\,\left[\frac{u^{5/2}}{2}\left(\frac{% \partial g_{1}}{\partial\rho}\zeta^{-1}+\frac{\partial g_{2}}{\partial\rho}% \zeta\right)+n_{I}\hat{a}_{0}^{\prime}\frac{\partial Q}{\partial\rho}\right]\,,∫ start_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT italic_d italic_u [ divide start_ARG italic_u start_POSTSUPERSCRIPT 5 / 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 end_ARG ( divide start_ARG ∂ italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG ∂ italic_ρ end_ARG italic_ζ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT + divide start_ARG ∂ italic_g start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG ∂ italic_ρ end_ARG italic_ζ ) + italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT over^ start_ARG italic_a end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT divide start_ARG ∂ italic_Q end_ARG start_ARG ∂ italic_ρ end_ARG ] , (4b)
00\displaystyle 0 =\displaystyle== ucdu[u5/2ζg2(p+2uc)+ζ1g1(p2uc)2\displaystyle\int_{u_{c}}^{\infty}du\bigg{[}u^{5/2}\frac{\zeta g_{2}(p+\frac{2% }{u_{c}})+\zeta^{-1}g_{1}(p-\frac{2}{u_{c}})}{2}∫ start_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT italic_d italic_u [ italic_u start_POSTSUPERSCRIPT 5 / 2 end_POSTSUPERSCRIPT divide start_ARG italic_ζ italic_g start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_p + divide start_ARG 2 end_ARG start_ARG italic_u start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_ARG ) + italic_ζ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_p - divide start_ARG 2 end_ARG start_ARG italic_u start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_ARG ) end_ARG start_ARG 2 end_ARG
\displaystyle-- nIqa^0u3+2uc33u2ucαk6uc2c1(uuc)3/2+3uc2u1/2fcζ1g12].\displaystyle n_{I}q\hat{a}_{0}^{\prime}\frac{u^{3}+2u_{c}^{3}}{3u^{2}u_{c}}-% \frac{\alpha k}{6u_{c}^{2}c_{1}(u-u_{c})^{3/2}}+\frac{3u_{c}^{2}}{u^{1/2}f_{c}% }\frac{\zeta^{-1}g_{1}}{2}\bigg{]}\,.italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT italic_q over^ start_ARG italic_a end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT divide start_ARG italic_u start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT + 2 italic_u start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG start_ARG 3 italic_u start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_u start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_ARG - divide start_ARG italic_α italic_k end_ARG start_ARG 6 italic_u start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_u - italic_u start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 3 / 2 end_POSTSUPERSCRIPT end_ARG + divide start_ARG 3 italic_u start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_u start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT italic_f start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_ARG divide start_ARG italic_ζ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG ] .

Then by tedious numerical calculations, the variables as functions of μ𝜇\muitalic_μ can be obtained. Especially, we overcome the numerical difficulties accoutered in Li_2015 when μ𝜇\muitalic_μ is low, which is crucial in determining the low pressure behavior of the EoS. The overall frame for the numerical calculation is shown in appendix A.

After the free energy Ω(μ)Ω𝜇\Omega(\mu)roman_Ω ( italic_μ ) is obtained, the pressure pQCDsubscript𝑝𝑄𝐶𝐷p_{QCD}italic_p start_POSTSUBSCRIPT italic_Q italic_C italic_D end_POSTSUBSCRIPT and energy density ϵQCDsubscriptitalic-ϵ𝑄𝐶𝐷\epsilon_{QCD}italic_ϵ start_POSTSUBSCRIPT italic_Q italic_C italic_D end_POSTSUBSCRIPT can be determined through the standard thermodynamic relations (at zero temperature):

pQCD=Ω,ϵQCD=Ω+nIμ.formulae-sequencesubscript𝑝𝑄𝐶𝐷Ωsubscriptitalic-ϵ𝑄𝐶𝐷Ωsubscript𝑛𝐼𝜇\displaystyle p_{QCD}=-\Omega,\quad\epsilon_{QCD}=\Omega+n_{I}\mu\,.italic_p start_POSTSUBSCRIPT italic_Q italic_C italic_D end_POSTSUBSCRIPT = - roman_Ω , italic_ϵ start_POSTSUBSCRIPT italic_Q italic_C italic_D end_POSTSUBSCRIPT = roman_Ω + italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT italic_μ . (5)

II.2 II.2 Equation of state

With pQCD(μ)subscript𝑝𝑄𝐶𝐷𝜇p_{QCD}(\mu)italic_p start_POSTSUBSCRIPT italic_Q italic_C italic_D end_POSTSUBSCRIPT ( italic_μ ) and ϵQCD(μ)subscriptitalic-ϵ𝑄𝐶𝐷𝜇\epsilon_{QCD}(\mu)italic_ϵ start_POSTSUBSCRIPT italic_Q italic_C italic_D end_POSTSUBSCRIPT ( italic_μ ) in hand, we can then extract the EoS. Though the result is numerical, it is accurate enough to fit it into an analytic form. We choose the range pQCD[0,0.05]subscript𝑝𝑄𝐶𝐷00.05p_{QCD}\in[0,0.05]italic_p start_POSTSUBSCRIPT italic_Q italic_C italic_D end_POSTSUBSCRIPT ∈ [ 0 , 0.05 ], since both the core pressure of NS and the baryon density are typical after recovering the dimensions, as explained in Zhang:2019tqd . In this range, the EoS can be described by the following doubly polytropic function:

ϵQCD1=0.140pQCD10.429+3.896pQCD11.335,pQCD1[0,0.05].formulae-sequencesubscriptitalic-ϵ𝑄𝐶𝐷10.140superscriptsubscript𝑝𝑄𝐶𝐷10.4293.896superscriptsubscript𝑝𝑄𝐶𝐷11.335subscript𝑝𝑄𝐶𝐷100.05\displaystyle\epsilon_{QCD1}=0.140p_{QCD1}^{0.429}+3.896p_{QCD1}^{1.335},\quad p% _{QCD1}\in[0,0.05].italic_ϵ start_POSTSUBSCRIPT italic_Q italic_C italic_D 1 end_POSTSUBSCRIPT = 0.140 italic_p start_POSTSUBSCRIPT italic_Q italic_C italic_D 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 0.429 end_POSTSUPERSCRIPT + 3.896 italic_p start_POSTSUBSCRIPT italic_Q italic_C italic_D 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1.335 end_POSTSUPERSCRIPT , italic_p start_POSTSUBSCRIPT italic_Q italic_C italic_D 1 end_POSTSUBSCRIPT ∈ [ 0 , 0.05 ] . (6)

We then need to recover the dimension, by noticing that

p=c2𝒩7pQCD,ϵ=𝒩7ϵQCD,formulae-sequence𝑝superscript𝑐2𝒩superscript7subscript𝑝𝑄𝐶𝐷italic-ϵ𝒩superscript7subscriptitalic-ϵ𝑄𝐶𝐷\displaystyle p=c^{2}\mathcal{N}\ell^{-7}p_{QCD},\quad\epsilon=\mathcal{N}\ell% ^{-7}\epsilon_{QCD}\,,italic_p = italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT caligraphic_N roman_ℓ start_POSTSUPERSCRIPT - 7 end_POSTSUPERSCRIPT italic_p start_POSTSUBSCRIPT italic_Q italic_C italic_D end_POSTSUBSCRIPT , italic_ϵ = caligraphic_N roman_ℓ start_POSTSUPERSCRIPT - 7 end_POSTSUPERSCRIPT italic_ϵ start_POSTSUBSCRIPT italic_Q italic_C italic_D end_POSTSUBSCRIPT , (7)

with c𝑐citalic_c the speed of light. By selecting the typical values λYMNc24.9similar-to-or-equalssubscript𝜆𝑌𝑀subscript𝑁𝑐24.9\lambda_{YM}N_{c}\simeq 24.9italic_λ start_POSTSUBSCRIPT italic_Y italic_M end_POSTSUBSCRIPT italic_N start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ≃ 24.9, MKK949MeVsimilar-to-or-equalssubscript𝑀𝐾𝐾949MeVM_{KK}\simeq 949\mathrm{MeV}italic_M start_POSTSUBSCRIPT italic_K italic_K end_POSTSUBSCRIPT ≃ 949 roman_M roman_e roman_V, we take 𝒩=1.2×1010MeV4.𝒩1.2superscript1010superscriptMeV4\mathcal{N}=1.2\times 10^{10}\mathrm{MeV}^{4}.caligraphic_N = 1.2 × 10 start_POSTSUPERSCRIPT 10 end_POSTSUPERSCRIPT roman_MeV start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT . Thus, the only adjustable parameter left is \ellroman_ℓ.

In order to apply this EoS more conveniently to NS, we can rewrite it in terms of the astrophysical units:

r=GNM/c2,ϵ=M/r3,p=c2ϵ.formulae-sequencesubscript𝑟direct-productsubscript𝐺𝑁subscript𝑀direct-productsuperscript𝑐2formulae-sequencesubscriptitalic-ϵdirect-productsubscript𝑀direct-productsuperscriptsubscript𝑟direct-product3subscript𝑝direct-productsuperscript𝑐2subscriptitalic-ϵdirect-product\displaystyle r_{\odot}=G_{N}M_{\odot}/c^{2},\quad\epsilon_{\odot}=M_{\odot}/r% _{\odot}^{3},\quad p_{\odot}=c^{2}\epsilon_{\odot}\,.italic_r start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT = italic_G start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT / italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , italic_ϵ start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT = italic_M start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT / italic_r start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT , italic_p start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT = italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ϵ start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT . (8)

Then, (7) becomes:

p/p=𝒜pQCD,ϵ/ϵ=𝒜ϵQCD,formulae-sequence𝑝subscript𝑝direct-product𝒜subscript𝑝𝑄𝐶𝐷italic-ϵsubscriptitalic-ϵdirect-product𝒜subscriptitalic-ϵ𝑄𝐶𝐷\displaystyle p/p_{\odot}=\mathcal{A}\;p_{QCD},\quad\epsilon/\epsilon_{\odot}=% \mathcal{A}\;\epsilon_{QCD},italic_p / italic_p start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT = caligraphic_A italic_p start_POSTSUBSCRIPT italic_Q italic_C italic_D end_POSTSUBSCRIPT , italic_ϵ / italic_ϵ start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT = caligraphic_A italic_ϵ start_POSTSUBSCRIPT italic_Q italic_C italic_D end_POSTSUBSCRIPT , (9)

where 𝒜=1.8×105×7𝒜1.8superscript105superscript7\mathcal{A}=1.8\times 10^{-5}\times\ell^{-7}caligraphic_A = 1.8 × 10 start_POSTSUPERSCRIPT - 5 end_POSTSUPERSCRIPT × roman_ℓ start_POSTSUPERSCRIPT - 7 end_POSTSUPERSCRIPT.

Consequently, the dimensionless EoS (10) can be converted in the astrophysical units:

ϵ1ϵ=0.140𝒜0.571(p1p)0.429+3.896𝒜0.335(p1p)1.335.subscriptitalic-ϵ1subscriptitalic-ϵdirect-product0.140superscript𝒜0.571superscriptsubscript𝑝1subscript𝑝direct-product0.4293.896superscript𝒜0.335superscriptsubscript𝑝1subscript𝑝direct-product1.335\displaystyle{\epsilon_{1}\over\epsilon_{\odot}}=0.140\mathcal{A}^{0.571}\left% ({p_{1}\over p_{\odot}}\right)^{0.429}+3.896\mathcal{A}^{-0.335}\left({p_{1}% \over p_{\odot}}\right)^{1.335}.divide start_ARG italic_ϵ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_ϵ start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT end_ARG = 0.140 caligraphic_A start_POSTSUPERSCRIPT 0.571 end_POSTSUPERSCRIPT ( divide start_ARG italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_p start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT end_ARG ) start_POSTSUPERSCRIPT 0.429 end_POSTSUPERSCRIPT + 3.896 caligraphic_A start_POSTSUPERSCRIPT - 0.335 end_POSTSUPERSCRIPT ( divide start_ARG italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_p start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT end_ARG ) start_POSTSUPERSCRIPT 1.335 end_POSTSUPERSCRIPT . (10)

We also listed the pointlike result in Zhang:2019tqd , for comparison later.

ϵ2ϵ=0.131𝒜0.544(p2p)0.456+2.629𝒜0.192(p2p)1.192.subscriptitalic-ϵ2subscriptitalic-ϵdirect-product0.131superscript𝒜0.544superscriptsubscript𝑝2subscript𝑝direct-product0.4562.629superscript𝒜0.192superscriptsubscript𝑝2subscript𝑝direct-product1.192\displaystyle{\epsilon_{2}\over\epsilon_{\odot}}=0.131\mathcal{A}^{0.544}\left% ({p_{2}\over p_{\odot}}\right)^{0.456}+2.629\mathcal{A}^{-0.192}\left({p_{2}% \over p_{\odot}}\right)^{1.192}.divide start_ARG italic_ϵ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG italic_ϵ start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT end_ARG = 0.131 caligraphic_A start_POSTSUPERSCRIPT 0.544 end_POSTSUPERSCRIPT ( divide start_ARG italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG italic_p start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT end_ARG ) start_POSTSUPERSCRIPT 0.456 end_POSTSUPERSCRIPT + 2.629 caligraphic_A start_POSTSUPERSCRIPT - 0.192 end_POSTSUPERSCRIPT ( divide start_ARG italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG italic_p start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT end_ARG ) start_POSTSUPERSCRIPT 1.192 end_POSTSUPERSCRIPT . (11)

III III. Gravitation Wave Constraints

The inner structure of NS can be solved with the help of the TOV equations given below, which can be obtained from perturbating the Einstein equation:

dmdr=4πr2ϵ,dpdr=GN(ϵ+p/c2)dϕdr,dϕdr=m+4πr3p/c2r(r2GNm/c2),missing-subexpression𝑑𝑚𝑑𝑟4𝜋superscript𝑟2italic-ϵmissing-subexpression𝑑𝑝𝑑𝑟subscript𝐺𝑁italic-ϵ𝑝superscript𝑐2𝑑italic-ϕ𝑑𝑟missing-subexpression𝑑italic-ϕ𝑑𝑟𝑚4𝜋superscript𝑟3𝑝superscript𝑐2𝑟𝑟2subscript𝐺𝑁𝑚superscript𝑐2\displaystyle\begin{aligned} &\frac{dm}{dr}=4\pi r^{2}\epsilon\,,\\ &\frac{dp}{dr}=-G_{N}\left(\epsilon+p/c^{2}\right)\frac{d\phi}{dr}\,,\\ &\frac{d\phi}{dr}=\frac{m+4\pi r^{3}p/c^{2}}{r\left(r-2G_{N}m/c^{2}\right)}\,,% \end{aligned}start_ROW start_CELL end_CELL start_CELL divide start_ARG italic_d italic_m end_ARG start_ARG italic_d italic_r end_ARG = 4 italic_π italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ϵ , end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL divide start_ARG italic_d italic_p end_ARG start_ARG italic_d italic_r end_ARG = - italic_G start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT ( italic_ϵ + italic_p / italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) divide start_ARG italic_d italic_ϕ end_ARG start_ARG italic_d italic_r end_ARG , end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL divide start_ARG italic_d italic_ϕ end_ARG start_ARG italic_d italic_r end_ARG = divide start_ARG italic_m + 4 italic_π italic_r start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_p / italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_r ( italic_r - 2 italic_G start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT italic_m / italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG , end_CELL end_ROW (12)

where ϕitalic-ϕ\phiitalic_ϕ is the metric potential and can be decoupled from the above. Altogether we have 4 variables as functions or r𝑟ritalic_r, but with only three equations. It is the EoS that fit the missing piece, which shows the relation between the energy density and the pressure.

What is more, there is another important property of the star, the TLN, or tidal deformability, which is defined as the dimensionless coefficient ΛΛ\Lambdaroman_Λ introduced as

Qij=ΛM5ij,subscript𝑄𝑖𝑗Λsuperscript𝑀5subscript𝑖𝑗\displaystyle Q_{ij}=-\Lambda\,M^{5}\,\mathcal{E}_{ij},italic_Q start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT = - roman_Λ italic_M start_POSTSUPERSCRIPT 5 end_POSTSUPERSCRIPT caligraphic_E start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT , (13)

with M𝑀Mitalic_M the star mass, Qijsubscript𝑄𝑖𝑗Q_{ij}italic_Q start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT the induced quadrupole moment, and ijsubscript𝑖𝑗\mathcal{E}_{ij}caligraphic_E start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT the external tidal field strength. This show the deformability of a star under an external gravitational field. By higher perturbations of Einstein equation, the TLN can also be calculated, under given EoS. More details are given in appendix B.

With EoS in Eqs.(10) and Eqs.(11) and TOV equations (12), we plot Figure,1 by setting different values of adjustable parameter \ellroman_ℓ. Figure.1 shows the relationship between mass and radius. We can see five curves, labeled as 1,2,3 for instanton gas case while A, B for pointlike case. For curves 1,2,3, we choose values as 7=9900,10300superscript7990010300\ell^{-7}=9900,10300roman_ℓ start_POSTSUPERSCRIPT - 7 end_POSTSUPERSCRIPT = 9900 , 10300, and 10700107001070010700, respectively. For A and B, the values are 7=10300superscript710300\ell^{-7}=10300roman_ℓ start_POSTSUPERSCRIPT - 7 end_POSTSUPERSCRIPT = 10300 and 17000170001700017000. We see that the maximal mass for instanton gas model with the three values are ranging from 1.8M1.8subscript𝑀direct-product1.8M_{\odot}1.8 italic_M start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT to 1.9M1.9subscript𝑀direct-product1.9M_{\odot}1.9 italic_M start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT. In contrast, solving EoS for pointlike baryons with 7=17000superscript717000\ell^{-7}=17000roman_ℓ start_POSTSUPERSCRIPT - 7 end_POSTSUPERSCRIPT = 17000, the maximal mass of NS is lower which is about 1.7M1.7subscript𝑀direct-product1.7M_{\odot}1.7 italic_M start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT. It is seen from the tendency of maximal mass and \ellroman_ℓ that we can adjust the mass to exceed 2M2subscript𝑀direct-product2M_{\odot}2 italic_M start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT by lowering \ellroman_ℓ, which is closer to the upper limit of NS mass. Actually, there exists another curve for pointlike case whose maximal mass is up to 2.2M2.2subscript𝑀direct-product2.2M_{\odot}2.2 italic_M start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT. In principle any mass can be achieved by adjusting \ellroman_ℓ, since there exists an scaling symmetry in M-R and TLN-M relations, as shown in Maselli:2017vfi ; Zhang:2023hxd .

The reason we did not utilize EoS of more massive NS is to limit the radius of 1.4M1.4subscript𝑀direct-product1.4M_{\odot}1.4 italic_M start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT NS to below 13km, as confirmed by multi-messenger observations 2017Multi . On the other hand, we choose such special values of \ellroman_ℓ so that the tidal deformability can meet the constraint from the analysis of GW170817. There are also other GW NS event like GW190425, but the TLN window are much bigger, so we only need to apply the one from GW170817.

The boundary of TLN from GW170817 and the properties of the five EoS will be stated more detailed now. Figure.2 shows the TLN-M relation for these different NS. The curves are plotted by combining with EoS and Eqs.(13). Detailed derivation of ΛΛ\Lambdaroman_Λ is referred to Appendix B. As introduced above, the analysis of GW170817 impose a constraint of TLN. For the NS binary system in this event, we can get the combined deformability with the two mass M1,2subscript𝑀12M_{1,2}italic_M start_POSTSUBSCRIPT 1 , 2 end_POSTSUBSCRIPT (M1>M2subscript𝑀1subscript𝑀2M_{1}>M_{2}italic_M start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT > italic_M start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT) and their respective TLN Λ1,2subscriptΛ12\Lambda_{1,2}roman_Λ start_POSTSUBSCRIPT 1 , 2 end_POSTSUBSCRIPT

Λ¯=1613(M1+12M2)M14Λ1+(M2+12M1)M24Λ2(M1+M2)5.¯Λ1613subscript𝑀112subscript𝑀2superscriptsubscript𝑀14subscriptΛ1subscript𝑀212subscript𝑀1superscriptsubscript𝑀24subscriptΛ2superscriptsubscript𝑀1subscript𝑀25\displaystyle{\bar{\Lambda}}=\frac{16}{13}\frac{\left(M_{1}+12M_{2}\right)M_{1% }^{4}\Lambda_{1}+\left(M_{2}+12M_{1}\right)M_{2}^{4}\Lambda_{2}}{\left(M_{1}+M% _{2}\right)^{5}}.over¯ start_ARG roman_Λ end_ARG = divide start_ARG 16 end_ARG start_ARG 13 end_ARG divide start_ARG ( italic_M start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + 12 italic_M start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) italic_M start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT roman_Λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + ( italic_M start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + 12 italic_M start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) italic_M start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT roman_Λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG ( italic_M start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_M start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 5 end_POSTSUPERSCRIPT end_ARG . (14)

In LIGOScientific:2018hze , a window Λ¯=300230+420¯Λsubscriptsuperscript300420230\bar{\Lambda}=300^{+420}_{-230}over¯ start_ARG roman_Λ end_ARG = 300 start_POSTSUPERSCRIPT + 420 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 230 end_POSTSUBSCRIPT is given. So we choose three values for instanton gas model around Λ¯=700¯Λ700\bar{\Lambda}=700over¯ start_ARG roman_Λ end_ARG = 700 as shown in Figure.3. From the analysis of event GW170817, we know that the lowest mass of two NS is 1.1M1.1subscript𝑀direct-product1.1M_{\odot}1.1 italic_M start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT while the total mass is 2.7M2.7subscript𝑀direct-product2.7M_{\odot}2.7 italic_M start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT. So the mass ranges from 1.35M1.35subscript𝑀direct-product1.35M_{\odot}1.35 italic_M start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT to 1.6M1.6subscript𝑀direct-product1.6M_{\odot}1.6 italic_M start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT in Figure.3. Compared among all the three instanton gas EoS, we found that the green line is the allowed realistic model with maximum mass, because the curve satisfies the upper bound of Λ¯¯Λ\bar{\Lambda}over¯ start_ARG roman_Λ end_ARG. Compared with curve 3, the maximal mass of curve 2 is larger, which is closer to the predicted ceiling of NS mass. While for curve 1, though with higher maximum mass, its TLN exceeds the constraint. For further comparison, we set the same \ellroman_ℓ for curve A. It can be seen that the EoS violates the windows much from all the three figures. We also adjust Λ¯¯Λ\bar{\Lambda}over¯ start_ARG roman_Λ end_ARG of curve B to the similar with curve 2, with the value =1700017000\ell=17000roman_ℓ = 17000, as shown in the pointlike case Zhang:2019tqd . All those figures imply that, the instanton gas model indeed shows better results,reproducing a maximum 1.85M1.85subscript𝑀direct-product1.85M_{\odot}1.85 italic_M start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT, while satisfying the GW constraints on tidal deformability.

Refer to caption
Figure 1: Mass-Radius relations for different EoS
Refer to caption
Figure 2: TLN (ΛΛ\Lambdaroman_Λ) vs Mass (M) for the holographic stars of EoS with the same sets of values and labels for \ellroman_ℓ as in Fig. 1.
Refer to caption
Figure 3: Tidal deformability of the binary holographic stars of EoS vs one of the masses (Λ¯¯Λ\bar{\Lambda}over¯ start_ARG roman_Λ end_ARG vs M1subscriptM1\rm{M}_{1}roman_M start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT) for partial sets of the values for \ellroman_ℓ used in Fig. 1.

IV IV. Conclusions

The SS model derived from AdS/QCD offers an holographic way to deal with high density nuclear matter, which is difficult in traditional approach. In this letter, we apply the instanton gas structure in SS model to extract the relation between the neutron energy density and the pressure, i.e., the EoS. This is an update for the point-like case in Zhang:2019tqd , with the instanton width ρ𝜌\rhoitalic_ρ as an extra variable. One thing to emphasize is that, ρ𝜌\rhoitalic_ρ is found to be also a function of the chemical potential μ𝜇\muitalic_μ, rather than a free parameter. Thus, the EoS obtained still has only one adjustable parameter, the D8- and D8¯¯D8\overline{\rm D8}over¯ start_ARG D8 end_ARG-branes separation distance \ellroman_ℓ. By applying the constraint from GW observation data, especially the binary NS event GW170817 (other NS events offer less sharp windows for TLN), we find that the maximum support mass is about 1.85M1.85subscript𝑀direct-product1.85M_{\odot}1.85 italic_M start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT. This value is below the known criterion of known observation of 2M2subscript𝑀direct-product2M_{\odot}2 italic_M start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT NS, but still offer a brand new direction to finally solve the theoretical derivation of EoS in future.

Acknowledgements. The authors thank Prof. Feng-Li Lin and Zhoujian Cao for helpful discussions. KZ (Hong Zhang) is supported by a classified fund from Shanghai city.

Appendix A A. Numerical Calculation Setup

In practice, the process of numerical calculation is constructed as follows. To obtain the pressure p𝑝pitalic_p and energy density ϵitalic-ϵ\epsilonitalic_ϵ in the EoS, we need to solve the following three equations.

0=02zI1(z,nI,ρ,k)3(1+z2)23𝑑z,0superscriptsubscript02𝑧subscript𝐼1𝑧subscript𝑛𝐼𝜌𝑘3superscript1superscript𝑧223differential-d𝑧0=\int_{0}^{\infty}\frac{2zI_{1}(z,n_{I},\rho,k)}{3(1+z^{2})^{\frac{2}{3}}}\,dz,0 = ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT divide start_ARG 2 italic_z italic_I start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_z , italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ , italic_k ) end_ARG start_ARG 3 ( 1 + italic_z start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT divide start_ARG 2 end_ARG start_ARG 3 end_ARG end_POSTSUPERSCRIPT end_ARG italic_d italic_z , (15)
0=02zI2(z,nI,ρ,k)3(1+z2)23𝑑z,0superscriptsubscript02𝑧subscript𝐼2𝑧subscript𝑛𝐼𝜌𝑘3superscript1superscript𝑧223differential-d𝑧0=\int_{0}^{\infty}\frac{2zI_{2}(z,n_{I},\rho,k)}{3(1+z^{2})^{\frac{2}{3}}}\,dz,0 = ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT divide start_ARG 2 italic_z italic_I start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_z , italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ , italic_k ) end_ARG start_ARG 3 ( 1 + italic_z start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT divide start_ARG 2 end_ARG start_ARG 3 end_ARG end_POSTSUPERSCRIPT end_ARG italic_d italic_z , (16)
μ(nI,ρ,k)=02zI4(z,nI,ρ,k)3(1+z2)23𝑑z.𝜇subscript𝑛𝐼𝜌𝑘superscriptsubscript02𝑧subscript𝐼4𝑧subscript𝑛𝐼𝜌𝑘3superscript1superscript𝑧223differential-d𝑧\mu(n_{I},\rho,k)=\int_{0}^{\infty}\frac{2zI_{4}(z,n_{I},\rho,k)}{3(1+z^{2})^{% \frac{2}{3}}}\,dz.italic_μ ( italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ , italic_k ) = ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT divide start_ARG 2 italic_z italic_I start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT ( italic_z , italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ , italic_k ) end_ARG start_ARG 3 ( 1 + italic_z start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT divide start_ARG 2 end_ARG start_ARG 3 end_ARG end_POSTSUPERSCRIPT end_ARG italic_d italic_z . (17)

Then we can read off

p(nI,ρ,k)=2702zI5(z,nI,ρ,k)3(1+z2)23𝑑z+μ(nI,ρ,k)nI,𝑝subscript𝑛𝐼𝜌𝑘27superscriptsubscript02𝑧subscript𝐼5𝑧subscript𝑛𝐼𝜌𝑘3superscript1superscript𝑧223differential-d𝑧𝜇subscript𝑛𝐼𝜌𝑘subscript𝑛𝐼p(n_{I},\rho,k)=\frac{2}{7}-\int_{0}^{\infty}\frac{2zI_{5}(z,n_{I},\rho,k)}{3(% 1+z^{2})^{\frac{2}{3}}}\,dz+\mu(n_{I},\rho,k)\cdot n_{I},italic_p ( italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ , italic_k ) = divide start_ARG 2 end_ARG start_ARG 7 end_ARG - ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT divide start_ARG 2 italic_z italic_I start_POSTSUBSCRIPT 5 end_POSTSUBSCRIPT ( italic_z , italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ , italic_k ) end_ARG start_ARG 3 ( 1 + italic_z start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT divide start_ARG 2 end_ARG start_ARG 3 end_ARG end_POSTSUPERSCRIPT end_ARG italic_d italic_z + italic_μ ( italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ , italic_k ) ⋅ italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , (18)

and

ϵ(nI,ρ,k)=27+02zI5(z,nI,ρ,k)3(1+z2)23𝑑z.italic-ϵsubscript𝑛𝐼𝜌𝑘27superscriptsubscript02𝑧subscript𝐼5𝑧subscript𝑛𝐼𝜌𝑘3superscript1superscript𝑧223differential-d𝑧\epsilon(n_{I},\rho,k)=-\frac{2}{7}+\int_{0}^{\infty}\frac{2zI_{5}(z,n_{I},% \rho,k)}{3(1+z^{2})^{\frac{2}{3}}}\,dz.italic_ϵ ( italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ , italic_k ) = - divide start_ARG 2 end_ARG start_ARG 7 end_ARG + ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT divide start_ARG 2 italic_z italic_I start_POSTSUBSCRIPT 5 end_POSTSUBSCRIPT ( italic_z , italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ , italic_k ) end_ARG start_ARG 3 ( 1 + italic_z start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT divide start_ARG 2 end_ARG start_ARG 3 end_ARG end_POSTSUPERSCRIPT end_ARG italic_d italic_z . (19)

In which, for the convenience of numerical integration in the equations, we can change the variable u𝑢uitalic_u to z𝑧zitalic_z with z𝑧zitalic_z defined as u(z)=(1+z2)13𝑢𝑧superscript1superscript𝑧213u(z)=(1+z^{2})^{\frac{1}{3}}italic_u ( italic_z ) = ( 1 + italic_z start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG 3 end_ARG end_POSTSUPERSCRIPT before proceeding with the numerical integration.

In Eq. (15)-(19), the functions Iisubscript𝐼𝑖I_{i}italic_I start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT are respectively

I1(z,nI,ρ,k)=(3nIρ3(4+4u(z)3ρ2)8u(z)12(1+u(z)3+ρ2)72)×(u(z)4fc(z)ζ(z,nI,ρ,k)+ζ(z,nI,ρ,k))3nI2fc(z)ρ3u(z)3(2+2u(z)3+3ρ2)4u(z)52(1+u(z)3+ρ2)4×ζ(z,nI,ρ,k)4u(z)52(1+u(z)3+ρ2)4,subscript𝐼1𝑧subscript𝑛𝐼𝜌𝑘3subscript𝑛𝐼superscript𝜌344𝑢superscript𝑧3superscript𝜌28𝑢superscript𝑧12superscript1𝑢superscript𝑧3superscript𝜌272𝑢superscript𝑧4subscript𝑓𝑐𝑧𝜁𝑧subscript𝑛𝐼𝜌𝑘𝜁𝑧subscript𝑛𝐼𝜌𝑘3superscriptsubscript𝑛𝐼2subscript𝑓𝑐𝑧superscript𝜌3𝑢superscript𝑧322𝑢superscript𝑧33superscript𝜌24𝑢superscript𝑧52superscript1𝑢superscript𝑧3superscript𝜌24𝜁𝑧subscript𝑛𝐼𝜌𝑘4𝑢superscript𝑧52superscript1𝑢superscript𝑧3superscript𝜌24\begin{split}I_{1}(z,n_{I},\rho,k)=&\left(\frac{3n_{I}\rho^{3}(-4+4u(z)^{3}-% \rho^{2})}{8u(z)^{\frac{1}{2}}(-1+u(z)^{3}+\rho^{2})^{\frac{7}{2}}}\right)\\ &\times\left(\frac{u(z)^{4}}{f_{c}(z)\zeta(z,n_{I},\rho,k)}+\zeta(z,n_{I},\rho% ,k)\right)\\ &-\frac{3n_{I}^{2}f_{c}(z)\rho^{3}u(z)^{3}(-2+2u(z)^{3}+3\rho^{2})}{4u(z)^{% \frac{5}{2}}(-1+u(z)^{3}+\rho^{2})^{4}}\\ &\times\frac{\zeta(z,n_{I},\rho,k)}{4u(z)^{\frac{5}{2}}(-1+u(z)^{3}+\rho^{2})^% {4}},\end{split}start_ROW start_CELL italic_I start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_z , italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ , italic_k ) = end_CELL start_CELL ( divide start_ARG 3 italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT italic_ρ start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT ( - 4 + 4 italic_u ( italic_z ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT - italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG 8 italic_u ( italic_z ) start_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT ( - 1 + italic_u ( italic_z ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT + italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT divide start_ARG 7 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT end_ARG ) end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL × ( divide start_ARG italic_u ( italic_z ) start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG start_ARG italic_f start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ( italic_z ) italic_ζ ( italic_z , italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ , italic_k ) end_ARG + italic_ζ ( italic_z , italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ , italic_k ) ) end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL - divide start_ARG 3 italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_f start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ( italic_z ) italic_ρ start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_u ( italic_z ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT ( - 2 + 2 italic_u ( italic_z ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT + 3 italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG 4 italic_u ( italic_z ) start_POSTSUPERSCRIPT divide start_ARG 5 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT ( - 1 + italic_u ( italic_z ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT + italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL × divide start_ARG italic_ζ ( italic_z , italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ , italic_k ) end_ARG start_ARG 4 italic_u ( italic_z ) start_POSTSUPERSCRIPT divide start_ARG 5 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT ( - 1 + italic_u ( italic_z ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT + italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG , end_CELL end_ROW (20)
I2(z,nI,ρ,k)=3nIρ416u(z)12(1+u(z)3+ρ2)72×((7+8u(z)3+3ρ2)ζ(z,nI,ρ,k)u(z)7(2111ρ221u(z)3+5u(z)3ρ2)(1+u(z)3)2ζ(z,nI,ρ,k))3nI2ρ4(u(z)3+2)(2+2u(z)3+3ρ2)8u(z)52(1+u(z)3+ρ2)4×ζ(z,nI,ρ,k)8u(z)52(1+u(z)3+ρ2)4nI12(1+3nI4ρk2)124ρ12(u(z)1)32,subscript𝐼2𝑧subscript𝑛𝐼𝜌𝑘3subscript𝑛𝐼superscript𝜌416𝑢superscript𝑧12superscript1𝑢superscript𝑧3superscript𝜌27278𝑢superscript𝑧33superscript𝜌2𝜁𝑧subscript𝑛𝐼𝜌𝑘𝑢superscript𝑧72111superscript𝜌221𝑢superscript𝑧35𝑢superscript𝑧3superscript𝜌2superscript1𝑢superscript𝑧32𝜁𝑧subscript𝑛𝐼𝜌𝑘3superscriptsubscript𝑛𝐼2superscript𝜌4𝑢superscript𝑧3222𝑢superscript𝑧33superscript𝜌28𝑢superscript𝑧52superscript1𝑢superscript𝑧3superscript𝜌24𝜁𝑧subscript𝑛𝐼𝜌𝑘8𝑢superscript𝑧52superscript1𝑢superscript𝑧3superscript𝜌24superscriptsubscript𝑛𝐼12superscript13subscript𝑛𝐼4𝜌superscript𝑘2124superscript𝜌12superscript𝑢𝑧132\begin{split}I_{2}(z,n_{I},\rho,k)=&\frac{3n_{I}\rho^{4}}{16u(z)^{\frac{1}{2}}% (-1+u(z)^{3}+\rho^{2})^{\frac{7}{2}}}\\ &\times((7+8u(z)^{3}+3\rho^{2})\zeta(z,n_{I},\rho,k)\\ &-\frac{u(z)^{7}(21-11\rho^{2}-21u(z)^{3}+5u(z)^{3}\rho^{2})}{(-1+u(z)^{3})^{2% }\zeta(z,n_{I},\rho,k)})\\ &-\frac{3n_{I}^{2}\rho^{4}(u(z)^{3}+2)(-2+2u(z)^{3}+3\rho^{2})}{8u(z)^{\frac{5% }{2}}(-1+u(z)^{3}+\rho^{2})^{4}}\\ &\times\frac{\zeta(z,n_{I},\rho,k)}{8u(z)^{\frac{5}{2}}(-1+u(z)^{3}+\rho^{2})^% {4}}\\ &-\frac{n_{I}^{\frac{1}{2}}\left(1+\frac{3n_{I}}{4\rho}-k^{2}\right)^{\frac{1}% {2}}}{4\rho^{\frac{1}{2}}(u(z)-1)^{\frac{3}{2}}},\end{split}start_ROW start_CELL italic_I start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_z , italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ , italic_k ) = end_CELL start_CELL divide start_ARG 3 italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT italic_ρ start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG start_ARG 16 italic_u ( italic_z ) start_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT ( - 1 + italic_u ( italic_z ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT + italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT divide start_ARG 7 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT end_ARG end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL × ( ( 7 + 8 italic_u ( italic_z ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT + 3 italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) italic_ζ ( italic_z , italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ , italic_k ) end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL - divide start_ARG italic_u ( italic_z ) start_POSTSUPERSCRIPT 7 end_POSTSUPERSCRIPT ( 21 - 11 italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 21 italic_u ( italic_z ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT + 5 italic_u ( italic_z ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG ( - 1 + italic_u ( italic_z ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ζ ( italic_z , italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ , italic_k ) end_ARG ) end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL - divide start_ARG 3 italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ρ start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT ( italic_u ( italic_z ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT + 2 ) ( - 2 + 2 italic_u ( italic_z ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT + 3 italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG 8 italic_u ( italic_z ) start_POSTSUPERSCRIPT divide start_ARG 5 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT ( - 1 + italic_u ( italic_z ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT + italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL × divide start_ARG italic_ζ ( italic_z , italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ , italic_k ) end_ARG start_ARG 8 italic_u ( italic_z ) start_POSTSUPERSCRIPT divide start_ARG 5 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT ( - 1 + italic_u ( italic_z ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT + italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL - divide start_ARG italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT start_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT ( 1 + divide start_ARG 3 italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT end_ARG start_ARG 4 italic_ρ end_ARG - italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT end_ARG start_ARG 4 italic_ρ start_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT ( italic_u ( italic_z ) - 1 ) start_POSTSUPERSCRIPT divide start_ARG 3 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT end_ARG , end_CELL end_ROW (21)
I3(z,nI,ρ,k)=kζ(z,nI,ρ,k)u(z)112,subscript𝐼3𝑧subscript𝑛𝐼𝜌𝑘𝑘𝜁𝑧subscript𝑛𝐼𝜌𝑘𝑢superscript𝑧112I_{3}(z,n_{I},\rho,k)=\frac{k\zeta(z,n_{I},\rho,k)}{u(z)^{\frac{11}{2}}},italic_I start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT ( italic_z , italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ , italic_k ) = divide start_ARG italic_k italic_ζ ( italic_z , italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ , italic_k ) end_ARG start_ARG italic_u ( italic_z ) start_POSTSUPERSCRIPT divide start_ARG 11 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT end_ARG , (22)
I4(z,nI,ρ,k)=(3ρ48u(z)12(1+u(z)3+ρ2)52)×(u(z)7(u(z)31)ζ(z,nI,ρ,k)+ζ(z,nI,ρ,k))+nI(u(z)31)(2+2u(z)3+3ρ2)24u(z)52(1+u(z)3+ρ2)3×ζ(z,nI,ρ,k)4u(z)52(1+u(z)3+ρ2)3,subscript𝐼4𝑧subscript𝑛𝐼𝜌𝑘3superscript𝜌48𝑢superscript𝑧12superscript1𝑢superscript𝑧3superscript𝜌252𝑢superscript𝑧7𝑢superscript𝑧31𝜁𝑧subscript𝑛𝐼𝜌𝑘𝜁𝑧subscript𝑛𝐼𝜌𝑘subscript𝑛𝐼𝑢superscript𝑧31superscript22𝑢superscript𝑧33superscript𝜌224𝑢superscript𝑧52superscript1𝑢superscript𝑧3superscript𝜌23𝜁𝑧subscript𝑛𝐼𝜌𝑘4𝑢superscript𝑧52superscript1𝑢superscript𝑧3superscript𝜌23\begin{split}I_{4}(z,n_{I},\rho,k)=&\left(\frac{3\rho^{4}}{8u(z)^{\frac{1}{2}}% (-1+u(z)^{3}+\rho^{2})^{\frac{5}{2}}}\right)\\ &\times\left(\frac{u(z)^{7}}{(u(z)^{3}-1)\zeta(z,n_{I},\rho,k)}+\zeta(z,n_{I},% \rho,k)\right)\\ &+\frac{n_{I}(u(z)^{3}-1)(-2+2u(z)^{3}+3\rho^{2})^{2}}{4u(z)^{\frac{5}{2}}(-1+% u(z)^{3}+\rho^{2})^{3}}\\ &\times\frac{\zeta(z,n_{I},\rho,k)}{4u(z)^{\frac{5}{2}}(-1+u(z)^{3}+\rho^{2})^% {3}},\end{split}start_ROW start_CELL italic_I start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT ( italic_z , italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ , italic_k ) = end_CELL start_CELL ( divide start_ARG 3 italic_ρ start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG start_ARG 8 italic_u ( italic_z ) start_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT ( - 1 + italic_u ( italic_z ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT + italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT divide start_ARG 5 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT end_ARG ) end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL × ( divide start_ARG italic_u ( italic_z ) start_POSTSUPERSCRIPT 7 end_POSTSUPERSCRIPT end_ARG start_ARG ( italic_u ( italic_z ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT - 1 ) italic_ζ ( italic_z , italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ , italic_k ) end_ARG + italic_ζ ( italic_z , italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ , italic_k ) ) end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL + divide start_ARG italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT ( italic_u ( italic_z ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT - 1 ) ( - 2 + 2 italic_u ( italic_z ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT + 3 italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 4 italic_u ( italic_z ) start_POSTSUPERSCRIPT divide start_ARG 5 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT ( - 1 + italic_u ( italic_z ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT + italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL × divide start_ARG italic_ζ ( italic_z , italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ , italic_k ) end_ARG start_ARG 4 italic_u ( italic_z ) start_POSTSUPERSCRIPT divide start_ARG 5 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT ( - 1 + italic_u ( italic_z ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT + italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG , end_CELL end_ROW (23)

and

I5(z,nI,ρ,k)=u(z)52+u(z)52ζ(z,nI,ρ,k)×(1+g2(z,nI,ρ)+(nIQ(z,ρ))2u(z)5).subscript𝐼5𝑧subscript𝑛𝐼𝜌𝑘𝑢superscript𝑧52𝑢superscript𝑧52𝜁𝑧subscript𝑛𝐼𝜌𝑘1subscript𝑔2𝑧subscript𝑛𝐼𝜌superscriptsubscript𝑛𝐼𝑄𝑧𝜌2𝑢superscript𝑧5\begin{split}I_{5}(z,n_{I},\rho,k)=&-u(z)^{\frac{5}{2}}+u(z)^{\frac{5}{2}}% \cdot\zeta(z,n_{I},\rho,k)\\ &\times\left(1+g_{2}(z,n_{I},\rho)+\frac{(n_{I}\cdot Q(z,\rho))^{2}}{u(z)^{5}}% \right).\end{split}start_ROW start_CELL italic_I start_POSTSUBSCRIPT 5 end_POSTSUBSCRIPT ( italic_z , italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ , italic_k ) = end_CELL start_CELL - italic_u ( italic_z ) start_POSTSUPERSCRIPT divide start_ARG 5 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT + italic_u ( italic_z ) start_POSTSUPERSCRIPT divide start_ARG 5 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT ⋅ italic_ζ ( italic_z , italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ , italic_k ) end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL × ( 1 + italic_g start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_z , italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ ) + divide start_ARG ( italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT ⋅ italic_Q ( italic_z , italic_ρ ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_u ( italic_z ) start_POSTSUPERSCRIPT 5 end_POSTSUPERSCRIPT end_ARG ) . end_CELL end_ROW (24)

The u(z)52𝑢superscript𝑧52-u(z)^{\frac{5}{2}}- italic_u ( italic_z ) start_POSTSUPERSCRIPT divide start_ARG 5 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT term in I5subscript𝐼5I_{5}italic_I start_POSTSUBSCRIPT 5 end_POSTSUBSCRIPT, together with the 2727\frac{2}{7}divide start_ARG 2 end_ARG start_ARG 7 end_ARG terms in (18) (19), are important, which are introduced to cancel the divergences at infinity. Those do not affect the physical meaning, since the divergent parts never change.

The terms fc(z)subscript𝑓𝑐𝑧f_{c}(z)italic_f start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ( italic_z ), Q(z,ρ)𝑄𝑧𝜌Q(z,\rho)italic_Q ( italic_z , italic_ρ ), q(z,ρ)𝑞𝑧𝜌q(z,\rho)italic_q ( italic_z , italic_ρ ), g1(z,nI,ρ)subscript𝑔1𝑧subscript𝑛𝐼𝜌g_{1}(z,n_{I},\rho)italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_z , italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ ), g2(z,nI,ρ)subscript𝑔2𝑧subscript𝑛𝐼𝜌g_{2}(z,n_{I},\rho)italic_g start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_z , italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ ) and ζ(z,nI,ρ,k)𝜁𝑧subscript𝑛𝐼𝜌𝑘\zeta(z,n_{I},\rho,k)italic_ζ ( italic_z , italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ , italic_k ) in Eq. (20)-(24) are defined as

fc(z)=1u(z)3,subscript𝑓𝑐𝑧1𝑢superscript𝑧3f_{c}(z)=1-u(z)^{-3},italic_f start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ( italic_z ) = 1 - italic_u ( italic_z ) start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT , (25)
Q(z,ρ)=(u(z)32)(fc(z)12)(3ρ2+2u(z)32)2((1+u(z)3+ρ2)32),𝑄𝑧𝜌𝑢superscript𝑧32subscript𝑓𝑐superscript𝑧123superscript𝜌22𝑢superscript𝑧322superscript1𝑢superscript𝑧3superscript𝜌232Q(z,\rho)=\frac{(u(z)^{\frac{3}{2}})(f_{c}(z)^{\frac{1}{2}})(3\rho^{2}+2u(z)^{% 3}-2)}{2((-1+u(z)^{3}+\rho^{2})^{\frac{3}{2}})},italic_Q ( italic_z , italic_ρ ) = divide start_ARG ( italic_u ( italic_z ) start_POSTSUPERSCRIPT divide start_ARG 3 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT ) ( italic_f start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ( italic_z ) start_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT ) ( 3 italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 2 italic_u ( italic_z ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT - 2 ) end_ARG start_ARG 2 ( ( - 1 + italic_u ( italic_z ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT + italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT divide start_ARG 3 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT ) end_ARG , (26)
q(z,ρ)=9u(z)12ρ44fc(z)12(1+u(z)3+ρ2)52,𝑞𝑧𝜌9𝑢superscript𝑧12superscript𝜌44subscript𝑓𝑐superscript𝑧12superscript1𝑢superscript𝑧3superscript𝜌252q(z,\rho)=\frac{9\cdot u(z)^{\frac{1}{2}}\cdot\rho^{4}}{4\cdot f_{c}(z)^{\frac% {1}{2}}\cdot(-1+u(z)^{3}+\rho^{2})^{\frac{5}{2}}},italic_q ( italic_z , italic_ρ ) = divide start_ARG 9 ⋅ italic_u ( italic_z ) start_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT ⋅ italic_ρ start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG start_ARG 4 ⋅ italic_f start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ( italic_z ) start_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT ⋅ ( - 1 + italic_u ( italic_z ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT + italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT divide start_ARG 5 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT end_ARG , (27)
g1(z,nI,ρ)=u(z)12nIq(z,ρ)3fc(z)12,subscript𝑔1𝑧subscript𝑛𝐼𝜌𝑢superscript𝑧12subscript𝑛𝐼𝑞𝑧𝜌3subscript𝑓𝑐superscript𝑧12g_{1}(z,n_{I},\rho)=\frac{u(z)^{\frac{1}{2}}\cdot n_{I}\cdot q(z,\rho)}{3\cdot f% _{c}(z)^{\frac{1}{2}}},italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_z , italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ ) = divide start_ARG italic_u ( italic_z ) start_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT ⋅ italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT ⋅ italic_q ( italic_z , italic_ρ ) end_ARG start_ARG 3 ⋅ italic_f start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ( italic_z ) start_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT end_ARG , (28)
g2(z,nI,ρ)=fc(z)12nIq(z,ρ)3u(z)72,subscript𝑔2𝑧subscript𝑛𝐼𝜌subscript𝑓𝑐superscript𝑧12subscript𝑛𝐼𝑞𝑧𝜌3𝑢superscript𝑧72g_{2}(z,n_{I},\rho)=\frac{f_{c}(z)^{\frac{1}{2}}\cdot n_{I}\cdot q(z,\rho)}{3% \cdot u(z)^{\frac{7}{2}}},italic_g start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_z , italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ ) = divide start_ARG italic_f start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ( italic_z ) start_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT ⋅ italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT ⋅ italic_q ( italic_z , italic_ρ ) end_ARG start_ARG 3 ⋅ italic_u ( italic_z ) start_POSTSUPERSCRIPT divide start_ARG 7 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT end_ARG , (29)

and

ζ(z,nI,ρ,k)=(1+g1(z,nI,ρ))12(1+g2(z,nI,ρ)k2u(z)8+(nIQ(z,ρ))2u(z)5)12.𝜁𝑧subscript𝑛𝐼𝜌𝑘superscript1subscript𝑔1𝑧subscript𝑛𝐼𝜌12superscript1subscript𝑔2𝑧subscript𝑛𝐼𝜌superscript𝑘2𝑢superscript𝑧8superscriptsubscript𝑛𝐼𝑄𝑧𝜌2𝑢superscript𝑧512\zeta(z,n_{I},\rho,k)=\frac{\left(1+g_{1}(z,n_{I},\rho)\right)^{\frac{1}{2}}}{% \left(1+g_{2}(z,n_{I},\rho)-\frac{k^{2}}{u(z)^{8}}+\frac{(n_{I}\cdot Q(z,\rho)% )^{2}}{u(z)^{5}}\right)^{\frac{1}{2}}}.italic_ζ ( italic_z , italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ , italic_k ) = divide start_ARG ( 1 + italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_z , italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ ) ) start_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 + italic_g start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_z , italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT , italic_ρ ) - divide start_ARG italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_u ( italic_z ) start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT end_ARG + divide start_ARG ( italic_n start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT ⋅ italic_Q ( italic_z , italic_ρ ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_u ( italic_z ) start_POSTSUPERSCRIPT 5 end_POSTSUPERSCRIPT end_ARG ) start_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT end_ARG . (30)

Appendix B B. Tidal Love Number

Under Boyer-Lindquist coordinate (t,r,θ.ϕ)formulae-sequence𝑡𝑟𝜃italic-ϕ(t,r,\theta.\phi)( italic_t , italic_r , italic_θ . italic_ϕ ) 1967JMP…..8..265B , we describe the metric of the sphericlal static NS with

ds2=gαβ(0)dxαdxβ=eν(r)dt2+eλ(r)dr2+r2(dθ2+sin2θdϕ2).𝑑superscript𝑠2superscriptsubscript𝑔𝛼𝛽0𝑑superscript𝑥𝛼𝑑superscript𝑥𝛽superscript𝑒𝜈𝑟𝑑superscript𝑡2superscript𝑒𝜆𝑟𝑑superscript𝑟2superscript𝑟2𝑑superscript𝜃2superscript2𝜃𝑑superscriptitalic-ϕ2ds^{2}=g_{\alpha\beta}^{(0)}dx^{\alpha}dx^{\beta}=-e^{\nu(r)}dt^{2}+e^{\lambda% (r)}dr^{2}+r^{2}\left(d\theta^{2}+\sin^{2}\theta d\phi^{2}\right).italic_d italic_s start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = italic_g start_POSTSUBSCRIPT italic_α italic_β end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT italic_d italic_x start_POSTSUPERSCRIPT italic_α end_POSTSUPERSCRIPT italic_d italic_x start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT = - italic_e start_POSTSUPERSCRIPT italic_ν ( italic_r ) end_POSTSUPERSCRIPT italic_d italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_e start_POSTSUPERSCRIPT italic_λ ( italic_r ) end_POSTSUPERSCRIPT italic_d italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_d italic_θ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ italic_d italic_ϕ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) . (31)

Setting GN=c=1subscript𝐺𝑁𝑐1G_{N}=c=1italic_G start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT = italic_c = 1, ν(r)𝜈𝑟\nu(r)italic_ν ( italic_r ) and λ(r)𝜆𝑟\lambda(r)italic_λ ( italic_r ) in the above equation are obtained by ν(r)2ϕ(r)=λ(r)=ln(12m(r)r)𝜈𝑟2italic-ϕ𝑟𝜆𝑟12𝑚𝑟𝑟\nu(r)\equiv 2\phi(r)=-\lambda(r)=\ln(1-\frac{2m(r)}{r})italic_ν ( italic_r ) ≡ 2 italic_ϕ ( italic_r ) = - italic_λ ( italic_r ) = roman_ln ( start_ARG 1 - divide start_ARG 2 italic_m ( italic_r ) end_ARG start_ARG italic_r end_ARG end_ARG ).

The complete metric is combined with the two parts

gαβ=gαβ(0)+hαβ,subscript𝑔𝛼𝛽superscriptsubscript𝑔𝛼𝛽0subscript𝛼𝛽\displaystyle g_{\alpha\beta}=g_{\alpha\beta}^{(0)}+h_{\alpha\beta},italic_g start_POSTSUBSCRIPT italic_α italic_β end_POSTSUBSCRIPT = italic_g start_POSTSUBSCRIPT italic_α italic_β end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT + italic_h start_POSTSUBSCRIPT italic_α italic_β end_POSTSUBSCRIPT , (32)

where hαβsubscript𝛼𝛽h_{\alpha\beta}italic_h start_POSTSUBSCRIPT italic_α italic_β end_POSTSUBSCRIPT is a linearized metric perturbation.

Under Reege-Wheeler gauge, the perturbation has the expression Hinderer_2009

hαβ=diag[eν(r)H0(r),eλ(r)H2(r),\displaystyle h_{\alpha\beta}=\operatorname{diag}[e^{-\nu(r)}H_{0}(r),e^{% \lambda(r)}H_{2}(r),italic_h start_POSTSUBSCRIPT italic_α italic_β end_POSTSUBSCRIPT = roman_diag [ italic_e start_POSTSUPERSCRIPT - italic_ν ( italic_r ) end_POSTSUPERSCRIPT italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_r ) , italic_e start_POSTSUPERSCRIPT italic_λ ( italic_r ) end_POSTSUPERSCRIPT italic_H start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_r ) ,
r2K(r),r2sin2θK(r)]Y2m(θ,φ).\displaystyle r^{2}K(r),r^{2}\sin^{2}\theta K(r)]Y_{2m}(\theta,\varphi).italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_K ( italic_r ) , italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ italic_K ( italic_r ) ] italic_Y start_POSTSUBSCRIPT 2 italic_m end_POSTSUBSCRIPT ( italic_θ , italic_φ ) . (33)

Here H0(r)=H2(r)=H(r)subscript𝐻0𝑟subscript𝐻2𝑟𝐻𝑟H_{0}(r)=H_{2}(r)=H(r)italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_r ) = italic_H start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_r ) = italic_H ( italic_r ) and K(r)𝐾𝑟K(r)italic_K ( italic_r ) is also related to H(r)𝐻𝑟H(r)italic_H ( italic_r ) by solving perturbative Einstein equation. On the other hand, H(r)𝐻𝑟H(r)italic_H ( italic_r ) can be derived by solving the following equation

H′′+H[2r+eλ(2m(r)r2+4πr(pρ))]superscript𝐻′′superscript𝐻delimited-[]2𝑟superscript𝑒𝜆2𝑚𝑟superscript𝑟24𝜋𝑟𝑝𝜌\displaystyle H^{\prime\prime}+H^{\prime}\left[\frac{2}{r}+e^{\lambda}\left(% \frac{2m(r)}{r^{2}}+4\pi r(p-\rho)\right)\right]italic_H start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT + italic_H start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT [ divide start_ARG 2 end_ARG start_ARG italic_r end_ARG + italic_e start_POSTSUPERSCRIPT italic_λ end_POSTSUPERSCRIPT ( divide start_ARG 2 italic_m ( italic_r ) end_ARG start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + 4 italic_π italic_r ( italic_p - italic_ρ ) ) ]
+H[6eλr2+4πeλ(5ρ+9p+ρ+p(dp/dρ))ν2]=0,𝐻delimited-[]6superscript𝑒𝜆superscript𝑟24𝜋superscript𝑒𝜆5𝜌9𝑝𝜌𝑝𝑑𝑝𝑑𝜌superscript𝜈20\displaystyle+H\left[-\frac{6e^{\lambda}}{r^{2}}+4\pi e^{\lambda}\left(5\rho+9% p+\frac{\rho+p}{(dp/d\rho)}\right)-\nu^{\prime 2}\right]=0,+ italic_H [ - divide start_ARG 6 italic_e start_POSTSUPERSCRIPT italic_λ end_POSTSUPERSCRIPT end_ARG start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + 4 italic_π italic_e start_POSTSUPERSCRIPT italic_λ end_POSTSUPERSCRIPT ( 5 italic_ρ + 9 italic_p + divide start_ARG italic_ρ + italic_p end_ARG start_ARG ( italic_d italic_p / italic_d italic_ρ ) end_ARG ) - italic_ν start_POSTSUPERSCRIPT ′ 2 end_POSTSUPERSCRIPT ] = 0 , (34)

At the core of NS where r=0𝑟0r=0italic_r = 0, H(r)𝐻𝑟H(r)italic_H ( italic_r ) has a solution

H(r)=r2[12πr27(5ρ(0)+9p(0)+ρ(0)+p(0)(dp/dρ)(0))\displaystyle H(r)=\mathcal{B}r^{2}[1-\frac{2\pi r^{2}}{7}(5\rho(0)+9p(0)+% \frac{\rho(0)+p(0)}{(dp/d\rho)(0)})italic_H ( italic_r ) = caligraphic_B italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [ 1 - divide start_ARG 2 italic_π italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 7 end_ARG ( 5 italic_ρ ( 0 ) + 9 italic_p ( 0 ) + divide start_ARG italic_ρ ( 0 ) + italic_p ( 0 ) end_ARG start_ARG ( italic_d italic_p / italic_d italic_ρ ) ( 0 ) end_ARG )
+O(r3)],\displaystyle+O(r^{3})],+ italic_O ( italic_r start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT ) ] , (35)

with \mathcal{B}caligraphic_B a constant and dp/dρ𝑑𝑝𝑑𝜌dp/d\rhoitalic_d italic_p / italic_d italic_ρ depends on the EoS.

At the surface of NS where r𝑟ritalic_r is equal to the whole radius R𝑅Ritalic_R, we get the Tidal Love Number

k28C55(12C)2[2+2C(y1)y]×\displaystyle k_{2}\equiv\frac{8C^{5}}{5}(1-2C)^{2}[2+2C(y-1)-y]\timesitalic_k start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ≡ divide start_ARG 8 italic_C start_POSTSUPERSCRIPT 5 end_POSTSUPERSCRIPT end_ARG start_ARG 5 end_ARG ( 1 - 2 italic_C ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [ 2 + 2 italic_C ( italic_y - 1 ) - italic_y ] ×
{2C(63y+3C(5y8))+4C3[1311y\displaystyle\{2C(6-3y+3C(5y-8))+4C^{3}[13-11y{ 2 italic_C ( 6 - 3 italic_y + 3 italic_C ( 5 italic_y - 8 ) ) + 4 italic_C start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT [ 13 - 11 italic_y
+C(3y2)+2C2(1+y)]\displaystyle+C(3y-2)+2C^{2}(1+y)]+ italic_C ( 3 italic_y - 2 ) + 2 italic_C start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( 1 + italic_y ) ]
+3(12C)2[2y+2C(y1)]log(12C)}1.\displaystyle+3(1-2C)^{2}[2-y+2C(y-1)]\log(1-2C)\}^{-1}.+ 3 ( 1 - 2 italic_C ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [ 2 - italic_y + 2 italic_C ( italic_y - 1 ) ] roman_log ( start_ARG 1 - 2 italic_C end_ARG ) } start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT . (36)

Here the compactness of star CM/R𝐶𝑀𝑅C\equiv M/Ritalic_C ≡ italic_M / italic_R with the total mass M𝑀Mitalic_M. The quantity yRH(R)/H(R)𝑦𝑅superscript𝐻𝑅𝐻𝑅y\equiv RH^{\prime}(R)/H(R)italic_y ≡ italic_R italic_H start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ( italic_R ) / italic_H ( italic_R ) is set for simplicity in calculation. Actually, we can alternatively start from the evaluation of Y(r)rH(r)/H(r)𝑌𝑟𝑟superscript𝐻𝑟𝐻𝑟Y(r)\equiv rH^{\prime}(r)/H(r)italic_Y ( italic_r ) ≡ italic_r italic_H start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ( italic_r ) / italic_H ( italic_r ) directly, then the second order differential equation becomes first order, and easier to solve, as shown in Postnikov:2010yn .

The dimensionless tidal deformability ΛΛ\Lambdaroman_Λ is is related to the l=2𝑙2l=2italic_l = 2 tidal Love number k2subscript𝑘2k_{2}italic_k start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT by

Λ=23k2(RM)5Λ23subscript𝑘2superscript𝑅𝑀5\Lambda=\frac{2}{3}k_{2}\left(\frac{R}{M}\right)^{5}\text{. }roman_Λ = divide start_ARG 2 end_ARG start_ARG 3 end_ARG italic_k start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( divide start_ARG italic_R end_ARG start_ARG italic_M end_ARG ) start_POSTSUPERSCRIPT 5 end_POSTSUPERSCRIPT . (37)

References

  • [1] Ben Margalit and Brian D. Metzger. The Multi-Messenger Matrix: the Future of Neutron Star Merger Constraints on the Nuclear Equation of State. Astrophys. J. Lett., 880(1):L15, 2019.
  • [2] Tanja Hinderer. Erratum: “tidal love numbers of neutron stars” (2008, apj, 677, 1216). The Astrophysical Journal, 697(1):964, may 2009.
  • [3] James M. Lattimer. The nuclear equation of state and neutron star masses. Ann. Rev. Nucl. Part. Sci., 62:485–515, 2012.
  • [4] B. S Sathyaprakash and Bernard F Schutz. Physics, astrophysics and cosmology with gravitational waves. Living Reviews in Relativity, 12(1):2, 2009.
  • [5] Duncan Galloway, Feryal Ozel, and Dimitrios Psaltis. Biases for neutron-star mass, radius and distance measurements from Eddington-limited X-ray bursts. Mon. Not. Roy. Astron. Soc., 387:268, 2008.
  • [6] Tim Dietrich, Michael W. Coughlin, Peter T. H. Pang, Mattia Bulla, Jack Heinzel, Lina Issa, Ingo Tews, and Sarah Antier. Multimessenger constraints on the neutron-star equation of state and the Hubble constant. Science, 370(6523):1450–1453, December 2020.
  • [7] B. P. Abbott et al. Multi-messenger Observations of a Binary Neutron Star Merger. Astrophys. J. Lett., 848(2):L12, 2017.
  • [8] B. P. Abbott et al. GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral. Phys. Rev. Lett., 119(16):161101, 2017.
  • [9] B. P. Abbott et al. Properties of the binary neutron star merger GW170817. Phys. Rev. X, 9(1):011001, 2019.
  • [10] B. P. Abbott et al. GW190425: Observation of a Compact Binary Coalescence with Total Mass 3.4Msimilar-toabsent3.4subscript𝑀direct-product\sim 3.4M_{\odot}∼ 3.4 italic_M start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT. Astrophys. J. Lett., 892(1):L3, 2020.
  • [11] Amirhossein Rezaei and Mohammad Parsa Akrami. Exact and Efficient Numerical approaches to MIT Bag Model. 3 2024.
  • [12] Kilar Zhang, Takayuki Hirayama, Ling-Wei Luo, and Feng-Li Lin. Compact Star of Holographic Nuclear Matter and GW170817. Phys. Lett. B, 801:135176, 2020.
  • [13] Tadakatsu Sakai and Shigeki Sugimoto. More on a holographic dual of QCD. Prog. Theor. Phys., 114:1083–1118, 2005.
  • [14] Tadakatsu Sakai and Shigeki Sugimoto. Low energy hadron physics in holographic QCD. Prog. Theor. Phys., 113:843–882, 2005.
  • [15] Edward Witten. Anti-de sitter space, thermal phase transition, and confinement in gauge theories. arXiv preprint hep-th/9803131, 1998.
  • [16] Joseph Polchinski and Matthew J. Strassler. The String dual of a confining four-dimensional gauge theory. 3 2000.
  • [17] Si-wen Li, Andreas Schmitt, and Qun Wang. From holography towards real-world nuclear matter. Physical Review D, 92(2), July 2015.
  • [18] Oren Bergman, Gilad Lifschytz, and Matthew Lippert. Holographic Nuclear Physics. JHEP, 11:056, 2007.
  • [19] Kazuo Ghoroku, Kouki Kubo, Motoi Tachibana, Tomoki Taminato, and Fumihiko Toyoda. Holographic cold nuclear matter as dilute instanton gas. Phys. Rev. D, 87(6):066006, 2013.
  • [20] Kilar Zhang and Feng-Li Lin. Constraint on hybrid stars with gravitational wave events. Universe, 6(12):231, 2020.
  • [21] Kilar Zhang, Guo-Zhang Huang, Jie-Shiun Tsao, and Feng-Li Lin. GW170817 and GW190425 as hybrid stars of dark and nuclear matter. Eur. Phys. J. C, 82(4):366, 2022.
  • [22] the homogeneous anstaz case does reproduce first order transition, but has no chiral restoration.
  • [23] M. J. Duff, B. E. W. Nilsson, and C. N. Pope. Kaluza-Klein Supergravity. Phys. Rept., 130:1–142, 1986.
  • [24] Andrea Maselli, Pantelis Pnigouras, Niklas Gronlund Nielsen, Chris Kouvaris, and Kostas D. Kokkotas. Dark stars: gravitational and electromagnetic observables. Phys. Rev. D, 96(2):023005, 2017.
  • [25] Kilar Zhang, Ling-Wei Luo, Jie-Shiun Tsao, Chian-Shu Chen, and Feng-Li Lin. Dark stars and gravitational waves: Topical review. Results Phys., 53:106967, 2023.
  • [26] Robert H. Boyer and Richard W. Lindquist. Maximal Analytic Extension of the Kerr Metric. Journal of Mathematical Physics, 8(2):265–281, February 1967.
  • [27] Sergey Postnikov, Madappa Prakash, and James M. Lattimer. Tidal Love Numbers of Neutron and Self-Bound Quark Stars. Phys. Rev. D, 82:024016, 2010.