Deriving Neutron Star Equation of State from AdS/QCD
Wei Li
wei_li_phys@163.comDepartment of Physics and Institute for Quantum Science and Technology, Shanghai University, 99 Shangda Road, Shanghai 200444, China
Jing-Yi Wu
wujingyi222@mails.ucas.ac.cnSchool of Astronomy and Space Science, University of
Chinese Academy of Sciences (UCAS), Beijing 100049, China
Kilar Zhang
kilar@shu.edu.cn, Corresponding AuthorDepartment 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 (solar masses), unfortunately not enough compared with the known observation of NS with over .
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 . This is higher than the pointlike case discussed before, yet still below the criteria of 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 goes to infinity, yet in practice we say 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 D8- and -branes around D4-branes background, connecting at a tip and maximally separated at , with the (dimensionless) holographic coordinate. Introducing the Kaluza-Klein Mass Duff:1986hr , The asymptotic separation is , which is 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.
represents the coordinate of the compactified dimension, where the radius is , thus its periodicity is given as . 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 .
This separation boundary condition can be mathematically described as
(1)
where the derivative mark ′ is taken with respect to .
The Abelian part of the gauge fields are denoted by , and its component gives the chemical potential by its boundary value, .
By Considering the gauge field action consisting a Dirac-Born-Infeld (DBI) and a Chern-Simons (CS) contribution,
,
the free energy density can be written as
(2)
with the abbreviation
(3)
Here is the baryon number density, k is an auxiliary integration constant introduced during taking the equations of motion for . The detailed definition of other symbols , and 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 . 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 :
(4a)
(4b)
Then by tedious numerical calculations, the variables as functions of can be obtained. Especially, we overcome the numerical difficulties accoutered in Li_2015 when 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 is obtained,
the pressure and energy density can be determined through the standard thermodynamic relations (at zero temperature):
(5)
II.2 II.2 Equation of state
With and 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 , 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:
(6)
We then need to recover the dimension, by noticing that
(7)
with the speed of light.
By selecting the typical values , , we take
Thus, the only adjustable parameter left is .
In order to apply this EoS more conveniently to NS, we can rewrite it in terms of the astrophysical units:
Consequently, the dimensionless EoS (10) can be converted in the astrophysical units:
(10)
We also listed the pointlike result in Zhang:2019tqd , for comparison later.
(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:
(12)
where is the metric potential and can be decoupled from the above. Altogether we have 4 variables as functions or , 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 introduced as
(13)
with the star mass, the induced quadrupole moment, and 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 .
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 , and , respectively. For A and B, the values are and . We see that the maximal mass for instanton gas model with the three values are ranging from to . In contrast, solving EoS for pointlike baryons with , the maximal mass of NS is lower which is about . It is seen from the tendency of maximal mass and that we can adjust the mass to exceed by lowering , which is closer to the upper limit of NS mass. Actually, there exists another curve for pointlike case whose maximal mass is up to . In principle any mass can be achieved by adjusting , 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 NS to below 13km, as confirmed by multi-messenger observations 2017Multi . On the other hand, we choose such special values of 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 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 () and their respective TLN
(14)
In LIGOScientific:2018hze , a window is given. So we choose three values for instanton gas model around as shown in Figure.3. From the analysis of event GW170817, we know that the lowest mass of two NS is while the total mass is . So the mass ranges from to 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 . 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 for curve A. It can be seen that the EoS violates the windows much from all the three figures. We also adjust of curve B to the similar with curve 2, with the value , as shown in the pointlike case Zhang:2019tqd . All those figures imply that, the instanton gas model indeed shows better results,reproducing a maximum , while satisfying the GW constraints on tidal deformability.
Figure 1: Mass-Radius relations for different EoSFigure 2: TLN () vs Mass (M) for the holographic stars of EoS
with the same sets of values and labels for as in Fig. 1.Figure 3: Tidal deformability of the binary holographic stars of EoS vs one of the
masses ( vs ) for partial sets of the values for 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 as an extra variable. One thing to emphasize is that, is found to be also a function of the chemical potential , rather than a free parameter. Thus, the EoS obtained still has only one adjustable parameter, the D8- and -branes separation distance .
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 . This value is below the known criterion of known observation of 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 and energy density in the EoS, we need to solve the following three equations.
(15)
(16)
(17)
Then we can read off
(18)
and
(19)
In which, for the convenience of numerical integration in the equations, we can change the variable to with defined as before proceeding with the numerical integration.
The term in , together with the 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 , , , , and in Eq. (20)-(24) are defined as
(25)
(26)
(27)
(28)
(29)
and
(30)
Appendix B B. Tidal Love Number
Under Boyer-Lindquist coordinate 1967JMP…..8..265B , we describe the metric of the sphericlal static NS with
(31)
Setting , and in the above equation are obtained by .
The complete metric is combined with the two parts
(32)
where is a linearized metric perturbation.
Under Reege-Wheeler gauge, the perturbation has the expression Hinderer_2009
(33)
Here and is also related to by solving
perturbative Einstein equation. On the other hand, can be derived by solving the following equation
(34)
At the core of NS where , has a solution
(35)
with a constant and depends on the EoS.
At the surface of NS where is equal to the whole radius , we get the Tidal Love Number
(36)
Here the compactness of star with the total mass . The quantity is set for simplicity in calculation. Actually, we can alternatively start from the evaluation of directly, then the second order differential equation becomes first order, and easier to solve, as shown in Postnikov:2010yn .
The dimensionless tidal deformability is is related to the tidal Love number by
(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 .
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.
[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.