Arxiv:1905.05296V1

Home , Information geometry

arXiv:1905.05296v1 [nucl-th] 13 May 2019 4 3 eutn itnei nivrerlto with relation equilibrium inverse between probability in fluctuation is the distance resulting in,ltrrfie nrf 6 ydtriigthe determining by “distance” [6] tensor, ref. the metric in of refined later measure statisti- tions, configura- for a equilibrium thermal metric different between a i.e. of systems, definition an cal 5] with [4, started initial thermodynamics, to methods equilibrium geometry information of application course, of [1–3]. physics and, optimiza- mathematics as processing, neuroscience, such tion, also signal sciences, but information learning, of inference machine areas statistical wider in in only not roles prob- of distributions. manifold ability a dually in with connections together affine coupled metric de- Riemannian It a geometrical inference. fines has statistical invariant [4] in of involved Rao structure studies by paper from seminal emerged the with started study to tool powerful systems. a statistical considered now is etry, INTRODUCTION I. Ca Mh,05.07 Pa,11.38 24.10 numbers: PACS densit baryon deconfinement/confinem small and and restoration/breaking temperature symmetry geometri high the between at comparison Chromodynamics The th masses. describes quark reliably finite for geometry Thermodynamic model. (NJL) nrp density. entropy nttt fPril n ula hsc,Fclyo Mathe of Faculty Physics, Nuclear and Particle of Institute iatmnod iia nvriàd oa“aSpez” P Sapienza”, “La Roma di Università Fisica, di Dipartimento h omls fReana emtyi ple osuythe study to applied is geometry Riemannian of formalism The 1 hroyai emty(G,aspecific a (TG), geometry Thermodynamic important play structures geometric These which [1–3], geometry information Indeed, geom- differential particular in and Geometry, hsdfiiinof definition This iatmnod iia nvriàd aai,VaSnaSo Santa Via Catania, di Università Fisica, di Dipartimento hroyai emtyo ab-Jn aii model Lasinio Jona Nambu- of Geometry Thermodynamic g µν hog h esa fthe of Hessian the through , 2 g oeoikc ,100Pau ,CehRepublic Czech 8, Prague 18000 2, Holešovičkách V µν NN ein iCtna -52 aai,Italy. Catania, I-95123 Catania, di Sezione INFN, .Castorina P. scuilsnethe since crucial is 2 , 4 .Lanteri D. , Abstract 1 esv aibe ihuiso oue eval- and volume, metric) a the of by units uated with variable, tensive e ic h J oe ie la indications clear gives per model interesting NJL only not the is since study se, car- This is model out. (NJL) Lasinio ried Jona - Nambu the to Res- models. Hadron the Gas with onance results the compare to and ped- rtcldcnnmn temperature to deconfinement the evaluate critical particular, to (pseudo-) density, baryon in low and and, perature tem- large theories at (QCD) field Quantum-Chromodynamics to plied [13]. holes transitions AdS black Hawking-Page Gauss-Bonnet-AdS charged the in in dyonic [12]; in hole black first transition like phase liquid-gas the order in liquid-liquid [11]; [9, and transitions fluids phase systems Lennard-Jones ferromagnetic the for for 10], [8], Argon Neon Hydrogen, and Helium, for coexistence phase in the to justification of theoretical meaning a physical offers the- and fluctuations ory classical the [7], generalizes developed which be can thermo- theory consistent fluctuation and dynamic covariant a Indeed, tem. orlto length, correlation bouevleo h clrcurvature scalar the of the value between absolute correspondence “interaction the the i.e to hypothesis”, leads it moreover, and, states nti ae ytmtcapiaino TD of application systematic a paper this In oercnl 1,1] Dhsbe ap- been has TD 15], [14, recently More systems: different many in tested been has TG 1 , n regimes. ent 2 a td fNLmdlado 21 Quantum (2+1) of and model NJL of study cal hs iga,bt ntecia ii and limit chiral the in both diagram, phase e .Mancani S. , hw la oncinbtenchiral between connection clear a shows y hs rniin nNmu-oaLasinio -Jona Nambu in transitions phase aisadPyis hre University Charles Physics, and matics azl loMr ,015Rm,Italy Roma, 00185 2, Moro Aldo iazzale a6,I913Ctna Italy. Catania, I-95123 64, fia ξ ftetemdnmcsys- thermodynamic the of , 3 R . ξ 3 h ueo the of cube the , R a in- (an T c on some dynamical mechanism, as chiral sym- where ∆aµ = aµ aµ and − U metry, for low energy QCD but also because a QCD fundamental property, quark confinement, ∂2s gµν = µ ν (3) is missing in NJL model with some interesting −∂a ∂a a=aU consequences on the geometrical description.

The TD approach is recalled in Sec. II and in defines a positive-definite Riemannian metric on Sec. III the phase diagram of the Nambu-Jona the space of thermodynamic states as the Hes- Lasinio model is discussed. Sec. IV is devoted sian of the entropy density, s, with respect its to the thermodynamic geometry description of natural variables aµ. chiral symmetry restoration in NJL model in the One can show [7] that previous formulas give chiral limit and for finite fermion masses. The ge- the probability of the spontaneous fluctuations ometrical difference in describing QCD and NJL between equilibrium states. Indeed, by expand- phase transitions is considered in Sec. V and ing eq. (1) up to second order for a aU , the ≃ Sec. VI contains our comments and conclusions. maximum entropy state, one finds the classical gaussian normalized fluctuation probability density: II. THERMODYNAMIC GEOMETRY n 2 n V P (a, aU ) d a = √gU In this section the procedure to define the 2 π × thermodynamic metric is briefly recalled (the de- V µ ν n exp gµν ∆a ∆a d a , tails are in ref. [7, 16]) and the description of × − 2 phase transitions by the scalar curvature, R, is (4) discussed, making also use of the application to n real fluids. being g the determinant of gµν and √g d x the usual invariant volume on a Riemannian manifold. A. Thermodynamic metric In the analysis of the phase transitions in NJL model by thermodynamic geometry we shall con- Let AU be a large thermodynamic system sider a two dimensional manifold, where the in- (universe) and let us consider an open subsystem tensive coordinates are β = 1/T and γ = µ/T , 0 − A with thermodynamic coordinates a , the inter- with µ chemical potential. Moreover the met- i nal energy density, and a , the number densities ric (3) turns out to be related with the deriva- of particles of different species. The probability tives of the potential φ = p/T , where p is the density to find A in the “point” a = (a0, a1, ) · · · pressure [16]: is given by φ φ g = ,ββ ,βγ , (5) P (a, a ) dna = C eSU (a,aU ) dna , (1) µν φ φ U ,βγ ,γγ being C a normalization constant, with the usual comma notation for derivatives. 0 i aU = (a , a , ) denotes the state of the The scalar curvature R simply becomes U U · · · universe and SU its total entropy, formally φ φ φ,γγ regarded as an exact function of the parameters 1 ,ββ ,βγ R = φ φ φ . (6) of A and A . 2 g2 ,βββ ,ββγ ,βγγ U φ φ φ ,ββγ ,βγγ ,γγγ On the basis of the maximum entropy prin- ciple and in the framework of Consistent and Covariant Fluctuation Theory (CCFT) [7], the B. Phase transition in thermodynamic thermodynamic properties of A can be studied geometry through the introduction of a quadratic form, The main results of the thermodynamic ge- 2 µ ν (∆ℓ) = gµν ∆a ∆a , (2) ometry within Ruppeiner’s formulation [7] are:

2 Figure 1: Schematic picture of the meaning of ξ: the intricate line represents the surface of ρ(r)= ρ0, i.e. that separating two sides with local mean densities ρ>ρ0 and ρ<ρ0. By tracing any straight line, the intersection points are separated by an average distance equal to ξ. Figure from [18]. Figure 2: Schematic pictures of different possible particle arrangements: (a) cluster of particles with 1) the (inverse) relation between the line element volume R pulled together by the attractive part of | | and the fluctuation probability between equilib- the interparticle interaction (R< 0); (b) a repulsive rium states; 2) the, so called, Interaction hypoth- solid-like cluster held up by hard-core particle repul- esis: the absolute value of the scalar curvature sion (R > 0); (c-d) a fluid in two phases near the critical point: the bottom half is a liquid phase con- R is proportional to a power of the correlation taining vapor droplets with volume Rl . The top half length, i.e. R ξd, where d is the effective spa- | | | |∼ is a coexisting vapor phase containing liquid droplets tial dimension of the underling thermodynamic with volume Rv . In (c) Rv = Rl and the droplets system. are commensurate,| | in (d)| liquid| | and| vapor phases The meaning of the correlation length and have incommensurate droplets; (e) liquid phase; (f) of the scalar curvature can be represented as in solid phase with R> 0. Figure from [18]. Fig.1 (a schematic picture due to Widom [17]): the intricate line represents what the surface of applications see Tab. I of Ref. [21]. density ρ(r)= ρ0 might look at any instant. This surface separates two sides with local mean den- The relation between R and ξd is easy to | | sities ρ > ρ0 and ρ < ρ0. By tracing any straight verify for second-order phase transitions, since R line, the intersection points with the surface ρ0 diverges, but the criterium to define a new phase are separated by an average distance equal to ξ. in term of the curvature R for a first order phase Because such points are separated by the same transition or a crossover is less clear. mean distance ξ, whatever the direction of the The approach called R-Crossing Method line, it is convenient to think that regions as vol- (RCM) [8] is often applied to define first order d ume elements (“droplets”) of dimension R ξ . phase transitions. It is based on the continuity ∼ Figure 2 shows a schematic summary of different of the scalar curvature: knowing the thermody- configurations. namic quantities in the two phases, i.e. R, one The interaction hypothesis has been con- can build up the transition curve by imposing firmed by the study of the classical ideal gas the continuity of R. The RCM, coherent with (R = 0 [6]) and of the van der Waals gas [7], Widom’s microscopic description of the liquid- for which, near the liquid-vapor critical point, gas coexistence region (i.e. with the idea that −2 Tc, the curvature is R (T Tc) /Tc . the correlation lengths of the two phases must ∼ | − | Other confirmations come from the study be the same at the transition) has been tested in of the Takahashi Gas [7], the Curie-Weiss systems with different features: vapor-liquid co- model [19], the ferromagnetic monodimensional existence line for the Lennard-Jones fluids [9, 10], Ising model [20]. For a more complete list of first and second order phase transitions of mean-

3 field Curie-Weiss model (ferromagnetic systems), ρ = N/V is the particle density) and the corre- liquid-liquid phase transitions [11], phase transi- sponding thermodynamic line element is given tions of cosmological interest as the liquid-gas- by [18] like first order phase transition in dyonic charged 1 ∂2f 1 ∂2f AdS black hole [12]. Another criterion, applied ∆ℓ2 = ∆T 2 + ∆ρ2 −T ∂T 2 T ∂ρ2 in the study of first order phase transitions in ρ T real fluids [18] and Lennard-Jones systems [10] (7) is a first kind discontinuity in R. The scalar curvature turns out to be Finally, two different phases can be linked 1 ∂ 1 ∂gρρ ∂ 1 ∂gT T by a crossover, as for the QCD deconfinement R = + , √g ∂T √g ∂T ∂ρ √g ∂ρ transition. Also in this case there is no defini- (8) tive conclusion on the behavior of R, although, with it has been recently shown [14] that the condition R = 0 predicts a temperature for the transi- 1 ∂2f 1 ∂2f gT T = , gρρ = . tion from QCD to the Hadron resonance Gas at −T ∂T 2 T ∂ρ2 ρ T low baryon density in agreement with freeze out (9) curve [22–24] and (within 10%) with lattice data and g = gT T gρρ. [25, 26]. In Ref. [18, 29, 33] the real fluid free energy is Another interesting aspect of the geometrical modeled on the NIST Chemistry WebBook and approach to phase transitions is that the sign of R is evaluated in the liquid and vapor phases and the scalar curvature brings information on the along the liquid-vapor coexistence curve ending microscopic interactions, since R turns out to be at the critical point Tc. positive for fermi statistical interactions and neg- At the critical point R with a power → −∞ ative in the bosonic case [27, 28]. Therefore a law behavior and in the asymptotic critical re- change in sign of R is an indication of the bal- gion, i.e very close to the critical temperature, ance between effective interactions, even when the values of the scalar curvature evaluated in the no transition occurs, and theoretical curves with two phases coincide. However in other regions of R = 0 in pure fluids identify some anomalous the thermodynamic parameter space the values behaviors observed in the experimental data of of R in the liquid and the vapor phases [18] are several substances (in particular, water) [18, 29]. quite different and mesoscopic fluctuating struc- A transition from R > 0 to R < 0 has been tures of different sizes occur in the two phases also shown for the Lennard-Jones system [9, 10] (see fig. 2.d). and Anyon gas [30, 31]. For black holes [32], the In the phase diagram of fluids, R is generally the change in sign of the curvature occurs at the found to be negative since the average molecu- Hawking-Page transition temperature, therefore lar distances are such that the attractive part of associated with the condition R = 0. the intermolecular potential dominates. However In the next sections we shall apply the ther- different anomalous regions, i.e. with R> 0, ex- modynamic geometry approach to NJL phase ist (see fig.4 in ref. [29]). They are localized: (a) diagram both in the chiral limit and for finite in the supercritical liquid region, near the melt- fermion mass. The behavior of the scalar curva- ing line; (b) in the liquid phase near the triple ture in the quantitative description of the critical point (for water); (c) in the vapor phase, in some line in the T µ plane will be pointed out. regions called “repulsive clusters” [29]. − The thermodynamic states for cases (a) and (b), named solid-like-liquid states, emerge when C. An example: real fluids the liquid organizes into solid-like structures at large densities, with a small intermolecular av- The geometrical study of fluids is based on erage separation. The states in “repulsive clus- the Helmholtz free energy per volume, f, in ter” areas (case c), are characterized by values terms of (T, ρ) coordinates (T is the temperature, of R much larger than the volume of a single

4 molecule and by low density and have been ob- with served in 97 different fluids (except those con- Ψ(T,µ) = 1 n (µ) n−(µ) (15) sisting of the simplest molecules) along the sat- − + − urated vapor phase curve. and 1 n±(µ)= . (16) E±µ III. NAMBU - JONA LASINIO MODEL 1+exp T n o For three flavors with mu = md = m, and ms = In Nambu–Jona Lasinio (NJL) model with 6 m, one has Mu = Md = Ms, and the SU(3) 6 two flavors (f = u, d), the SU(2) lagrangian [34– lagrangian is [36] 36] is given by = ψ (i ∂/ m) ψ + + , (17) LSU(3) − L4 L6 =ψ (i ∂/ m) ψf + LSU(2) f − where 2 2 b + G ψ ψf + ψ iγ5 τ ψf , 2 2 f f −→ = G ψλa ψ + ψiγ λaψ (18) L4 5 h i(10) a X h i and the ’t Hooft interaction, 6, is given by being G a dimensionful coupling, m the current L quark mass (m = 0 is the chiral limit) and −→τ the 6 = K det ψ (1 + γ5) ψ + det ψ (1 γ5) ψ , Pauli matrices. In mean-field approximation the L − " − # thermodynamic potential, Ω, at finite tempera- (19) T ture and chemical potential turns out to be [36] with ψ = (u, d, s) , m = diag (m,m,ms), λ0 = 2/3 1 × , being 1 × the 3 3 identity matrix, 2 3 3 3 3 × (Mf m) and where λ (a = 1, ..., 8) are the Gell-Mann Ω(Mf )= − + Nf Ωf , (11) p a b 4 G matrices and K and G dimensionful couplings. with The gap equations,

3 Mi = mi 4 G ψ ψi + 2 K ψ ψj ψ ψk d p − i j k Ωf = 2 Nc Ef − (2 π)3 − (j, k = i) , Z 6 3 E +µ d p − f f (20) 2 Nc T ln 1+ e T − (2 π)3 − are coupled with the quark condensates Z − 3 Ef µf 3 d p − d p Mi 2 Nc T 1+ e T , ψ ψ = 2 N Ψ , (21) − (2 π)3 i i c 3 i Z − (2 π) Ei (12) Z where where Mf is the dynamically generated mass, 1 1 Ψ = 1 − , (22) i Ei+µi Ei µi 2 2 − T − T Ef = p + Mf , Nc and Nf are the number 1+ e 1+ e of colorsq and flavors respectively, µf is the quark and the mean-field thermodynamic potential Ω f chemical potential and the integrals are regu- turns out to be [36] lated by a cutoff Λ. For mu = md, µ = µu = µd, 2 Ω= Ωf + 2 G ψf ψf the generated quark mass is M = Mu = Md . − f=u,d,s f=u,d,s (23) To evaluate the minimum of Ω by eq. (11), X X 4 K uu dd ss , one has to solve the self-consistent gap equation − h i h i with Ωf in eq. (12). M = m 2 G ψψ , (13) − Finally, the potential we need for the thermodynamic geometry calculations is where ψψ is the quark-antiquark condensate: P φ(β, γ)= = Ω(β, γ) β , (24) d3p M T − ψψ = 2 Nc Nf Ψ(T,µ) , (14) − (2 π)3 E where P = Ω is the pressure. Z −

5 IV. THERMODYNAMIC GEOMETRY OF (a) CHIRAL SYMMETRY RESTORATION IN 300 NJL MODEL 250

) 200 1. Two flavors in the chiral limit MeV ( 150 M Let us first discuss the chiral limit (m = 0) 100 for two flavors, starting from the breaking of chi- 50 ral symmetry at T = µ = 0, with the value of 0 the dynamical mass M0(0, 0) = 300 MeV , cor- 0 50 100 150 200 responding to Λ = 650 MeV and G = 5.01 Temperature T (MeV) − − × (b) 10 6 MeV 2 [34, 35]. 300 The well known solution M(T,µ) of the gap equation (13), for different values of the tempera- 250 ) ture and of the quark chemical potential, is plot- 200 MeV ted in Fig. 3.a and 3.b. The restoration of the ( 150 chiral symmetry is a first order phase transition M 100 at large chemical potential and a second order one at low µ. 50 The study of the critical line of the symmetry 0 0 100 200 300 400 restoration, T (µ), by thermodynamic geometry Chemical Potential μ (MeV) requires the, straightforward but laborious, calculation of the scalar curvature R, reported in Figure 3: a) The dynamically generated mass, M, appendix A. in the NJL model with two favors in the chiral limit It turns out that R diverges at the critical (mu = md = 0 MeV) and again the temperature. | | Black line is for µ =0 MeV; the others are for growing temperature, i.e. there is a second order phase ⋆ µ, up to µ = 300 MeV and with step of ∆µ = 20 MeV. transition, for µ <µ 290 MeV, as shown b) M as a function of the chemical potential µ. Black ≃ ⋆ in fig. 4 for µ = 0. For µ>µ there is, in- line is for T = 10 MeV; the others are for growing T , stead, a first order phase transition. The dynam- up to T = 170 MeV and with step of ∆T = 20 MeV. ically generated mass, M, now takes the charac- teristic behavior plotted in figure 5, where the attractive, suggesting that the chiral symmetry black curves (both the continuous and the dot- restoration is due to thermal fluctuations. ted) are for T = 30 MeV and the two light-gray On the other hand, at large chemical po- lines define the spinodal points. Between the two tential R turns out to be positive, indicating a spinodal (light-gray) lines one can evaluate three screening of the potential and an increase of the different scalar curvatures: the first one for the repulsive interaction at large density. higher-mass branch (black curve in figure 5); the The complete critical line obtained by ther- second one for M = 0 MeV and the last one is re- modynamic geometry is depicted in Figure 6 lated to the M-branch that interpolates between where the continuous line shows the II order M = 0 and the upper M-curve (dotted curve in phase transition and the dashed lines the spin- figure 5). At fixed temperature and between the odal curves of the first order one. The green band spinodal lines (see fig. 5), there is a discontinuity is the region of negative R. in R which identifies the two dashed curves in | | fig. 6. The crossing temperature from the I order 2. Two flavors with chiral masses phase transition to the II order turns out to be about 58 MeV. With finite chiral quark masses, at high tem- For small µ and near the transition the cur- perature and low chemical potential, there is vature is negative, i.e. the interaction is mostly a smooth crossover rather than a second-order

6 μ=0 MeV 1.5 )

fm 150 ( 1.0 3  1 R |

0.5 ) MeV

( 100 0.0 T ■ R<0 100 150 200 250 Temperature T (MeV)

Figure 4: R from µ = 0 MeV: second order phase transition. 50

300 T=30 MeV 250

) 200 0 50 100 150 200 250 300 350 MeV ( 150 μ (MeV) M 100 Figure 6: The transition temperature: continuous 50 line is for II order phase transition and the dashed 0 ones for the first order one. The transition point is 270 280 290 300 310 320 330 ⋆ ⋆ at µχ = 290 MeV and Tχ = 58 MeV. The green band Chemical Potential μ (MeV) is the region of R< 0 Figure 5: The dynamically generated mass M in the 2 flavors NJL chiral model and temperature T = 30 MeV. agreement with that one obtained by chiral susceptibility (see eq. (A23) in appendix A). phase transition. Moreover, the first-order phase boundary ends in a second-order endpoint [36]. For low temperature and large chemical po- The solution of the gap equation (13) (with tential, the scalar curvature R has the same be- Λ = 650 MeV and G = 5.01 10−6 MeV−2 and havior previously discussed in the chiral limit, × m0 = 5.5 MeV) as a function of T and µ is shown i.e. a first order phase transition. in Figs. 7.a and 7.b. The critical point, (T ⋆,µ⋆) between the To clarify the effect of the chiral mass in the crossover and the first order phase transition de- calculation of the scalar curvature, Fig. 8 shows pends on m0 and for (the generally accepted that R diverges in the chiral limit but for m0 = 0, value) m = 5.5 MeV one has µ⋆ 329 MeV 6 0 ≃ near the transition temperature, it has a mini- and T ⋆ 32 MeV. mum, corresponding to a maximum of R , i.e. ∼ | | to a finite correlation length. Therefore, m = 0 Figure 9 shows the critical line for m = 0 6 0 changes the behavior of R near the critical tem- 5.5 MeV: the continuous line is obtained by the perature: the divergence of the II order phase maximum of R and the dashed ones are the | | transition turns into a minimum in the negative spinodal curves. The black circle is at µ⋆ = R region and the transition temperature eval- 329 MeV and T ⋆ = 32 MeV. The green band uated by the maximum of R is completely in is the region of R< 0. | |

7 (a) 0.4 300

250 0.2 ) 200 0.0 MeV (

150 ) 3 M

fm -0.2 100 ( R μ=0 MeV 50 -0.4 m0=0 MeV 0 0 50 100 150 200 250 300 350 -0.6 m0=2.5 MeV Temperature T (MeV) m0=5.5 MeV (b) -0.8 300 120 140 160 180 200 220 240 250 Temperature T (MeV)

) 200

MeV Figure 8: R from µ = 0 MeV and diﬀerent val- ( 150

M ues of the bare mass m0: continuous line is from 100 m0 = 0 MeV (the chiral limit) and R shows a negative divergence. Dashed line is from m0 = 2.5 MeV 50 and the dotted from m0 = 5.5 MeV; both show a 0 finite region with negative R around the transition 0 100 200 300 400 Chemical Potential μ (MeV) temperature, which corresponds to the local maximum of R . | | Figure 7: a) The dynamical generated mass, M, in the NJL model with two favors of identical mass behavior of the scalar curvature is essentially the (mu = md = 5.5 MeV) and again the temperature. Black line is for µ =0 MeV; the others are for grow- same of the previous case with two flavors and ing µ, up to µ = 340 MeV and with step of ∆µ = 20 physical masses. MeV. b) M as a function of the chemical potential µ. In figure 10 the ratios χs/χsmax (dashed line), Black line is for T = 10 MeV; the others are for grow- χu/χumax (dotted line) and R / R max (continu- | | | | ing T , up to T = 400 MeV and with step of ∆T = 20 ous line) are depicted to visualize that the max- MeV. imum in R corresponds to the peak of chiral | | susceptibilities. 3. Three flavors Figure 11 shows the transition temperature by the evaluation of R: the continuous line is Three flavor NJL model is studied with the again obtained by the maximum of R and the | | parameter values [37] dashed ones are the spinodal curves. The black circle is at µ⋆ 335 MeV and T ⋆ 35 MeV. The ∼ ∼ Λ = 631.4 MeV G Λ2 = 1.835 green band is the region of negative R. K Λ5 = 9.29 (25)

m = 5.5 MeV ms = 135.7 MeV A. Thermal geometric definition of the phase transitions in NJL model: summary and only one chemical potential (µ = µd = µu, µs = 0). The dynamically generated masses It is useful to conclude this section by sum- Mu = Md and Ms are now solutions of the sys- marizing the geometrical definition of the phase tem of eq. (20) and eq. (21). Their behavior is transitions: similar to that one depicted in fig. 7, but with different values for light and strange quarks. Also a II order phase transition occurs for two • in this case there is a crossover at low chemi- flavors in the chiral limit (m = 0) at low cal potential and large T and a first order phase chemical potential. This transition is char- transition at low temperature and large µ. The acterized by a divergent scalar curvature;

8 250 250

■ R<0

■ R<0 200 200

) 150 ) 150 MeV MeV ( ( T T

100 100

50 50

0 0 0 50 100 150 200 250 300 350 0 50 100 150 200 250 300 350 μ (MeV) μ (MeV)

Figure 9: The transition temperature by the R con- Figure 11: The transition temperature by the R ditions and from m0 = 5.5 MeV: continuous line is conditions: continuous line is obtained by the lo- obtained by the local maximum of R , the dashed cal maximum of R , the dashed ones indicate the ones indicate the spinodal lines. The| | circle is at spinodal lines. The| | circle is at µ⋆ 335 MeV and µ⋆ = 329 MeV and T ⋆ = 32 MeV. The green band is T ⋆ 35 MeV.The green band is the∼ region of R< 0 the region of R< 0. ∼

for chiral masses, there is a crossover, both • for two and three ﬂavors, at low chemical potential and large T . The transition tem- 1.0 perature is deﬁned as the maximum of R | | in the negative-R region and it is in agree- 0.8 ment with the chiral susceptibility analysis χ [38] (eqs. (A23), (B36) and (B37) in 0.6 appendix);

0.4 there exists a I order phase transition at χs/χs max • low temperature and large µ, both with χu/χu max 0.2 two and three ﬂavors and both in the chiral |R|/|Rmax| limit or with chiral masses. This transition 0.0 180 190 200 210 220 is related with a discontinuity in R. T (MeV)

Figure 10: The ratio χs/χsmax (dashed line), V. NJL MODEL AND QCD CROSSOVER χu/χumax (dotted line) and R / R max (continuous | | | | line) at µ =0 MeV. As seen in Sec. IV, the NJL crossover can be identiﬁed by a local maximum of R . However, | |

9 200 VI. COMMENTS AND CONCLUSIONS

180 Thermodynamic geometry reliably describes ms/ml= χ for 20 the phase diagram of NJL model, both in χ for ms/ml=27 160 the chiral limit and for finite mass, and indicates a geometrical interplay between chiral symmetry restoration/breaking and deconfine- 140 ment/confinement regimes. Moreover in a very recent paper [43] the chi- χ 0 120 ral phase transition temperature Tc , corresponding to a “true” chiral transition in the limit ms/ml >> 1, turns out to be about 25 MeV less 100 than the pseudo-critical temperature. Fig 12 suggests that a small variation from 80 ms/ml = 20 to ms/ml = 27 changes the maximum of chiral susceptibility from R = 0 to a 60 finite value of R , as in NJL model. It could | | be possible that considering the effective chiral -1.0 -0.8 -0.6 -0.4 -0.2 0.0 limit, i.e. ms/ml >> 1 one recovers by thermo- 3 R (fm ) dynamic geometry a “true” chiral phase transi- Figure 12: The chiral susceptibility χ at µ =0 MeV tion at lower temperature, with typical scaling and as a function of the scalar curvature R for laws. The role of color confinement in QCD in physical value of the strange quark mass, ms, and terms of thermodynamic geometry will be dis- ms/mℓ = 20 (dotted line) or ms/mℓ = 27 (continu- cussed in different models in a forthcoming pa- ous line). per.

Acknowledgements The authors thank H.Satz for useful comments. NJL model misses color confinement and therefore there is no a priori reason to apply the same geometric criterium for non perturbative QCD Appendix A: NJL model with two quarks dynamics. To evaluate the scalar curvature R one needs Indeed, in the thermodynamic geometry de- the derivatives of the potential φ, up to third scription of QCD deconfinement transition in order, which can be written in terms of the dy- ref. [14] the criterium R = 0 has been used. It namical generated mass M. Therefore, the solu- indicates the transition from a mostly fermionic tion of the GAP equation uniquely determines all system (as the quark-gluon plasma) to an essen- those functions. Indeed, after a straightforward tially bosonic one (as the hadron resonance gas) calculation, one gets (a comma indicates partial but, as shown in fig. 12, it exactly corresponds to derivative) the maximum of chiral susceptibility, confirming the well known [39–42] interplay between con- b M finement and chiral symmetry breaking. M = 1 , (A1) ,β 1 f f M 2 − 1 − 2 The transition temperature evaluated by R = 0 is in agreement with the freeze-out hadronization curve and with the pseudo-critical g M M = 1 , (A2) temperature by lattice data within 10% [14]. ,γ 1 f f M 2 − 1 − 2 10 Λ n− (1 2n−)+ n (1 2n ) M =d b M + (b + f ) M + 2 ,β +,β + ,ββ 3 2 1,β ,β b4 =κM dpp − 2 − , " 0 T E Z (A13) + (b + 2 Tf ) M 2 M + f M M 2 , 4 2 ,β 3 ,β Λ # 2 n− (1 n−) n+ (1 n+) g1 = κM dp p − − − , (A3) 0 E Z (A14)

Λ Ψγ 2 g = κ dp p4 , (A15) M,γγ =d g3 M + (g2 + f1,γ) M,γ + g4 M M,γ+ 2 M E3 " Z0

2 + f3 M M,γ , Λ 2 n−,γ (1 2n−) n+,γ (1 2n+) # g3 =κM dpp − − − , (A4) 0 E Z (A16)

Λ 2 n−,γ (1 2n−)+ n+,γ (1 2n+) 2 g4 =κM dpp − 2 − , M,βγ =d g5 M + b2 M,γ + f1,γM,β + g4 M M,β+ 0 T E " Z (A17)

+ f T M 2 M + f M M M , Λ 2 ,γ 3 ,β ,γ 2 # g5 =κM dpp [n−,γ (1 2n−)+ n+,γ (1 2n+)] , 0 − − (A5) Z (A18) with Nc Nf κM = 2 G (A19) −1 π2 d = 1 f f M 2 , (A6) − 1 − 2 and n in Eq. (16). ± Λ 4 By deriving Eq. (24) and Eq. (11) and defin- p Ψ ing f1 = κM dp 3 , (A7) 0 E Z κ κ = M , (A20) Ω 2 G Λ 2 n− (1 n−)+ n+ (1 n+) f2 = κM dp p − 2 − , one gets 0 T E Z (A8) Λ (M m)2 φ = κ dp p2 E Ψ , (A21) Λ ,β Ω − n− (1 n−)+ n (1 n ) 0 − 4 G 4 + + Z f3 = κM dp p − 4 − , 0 T E Z (A9) Λ 2 Λ φ,γ = κΩ dp p (n+ n−) . (A22) 2 0 − b1 = κM dp p [n− (1 n−)+ n+ (1 n+)] , Z 0 − − Z (A10) The calculation of second and third order derivatives is straightforward. Λ Ψ Finally, the two flavors chiral susceptibility, b = κ dp p4 β , (A11) 2 M E3 χ, is defined as [38] Z0 2f ∂M 1 Λ χ = = 2 = 2 ∂m 1 f1 f2 M b3 =κM dpp [n−,β (1 2n−)+ n+,β (1 2n+)] , − − (A23) − − M,β M,γ Z0 = = . (A12) b1 M g1 M

11 Appendix B: Three ﬂavors Mu,γ (δ bu Mu ǫ)= cu ǫ Ms,γ ζ − −  , In a three ﬂavors systems the derivatives of   (cs θ cu λ) (δ bu Mu ǫ) cu ǫ bu λ Mu the dynamically generated mass Mu = Md and Ms,γ = − − − ( η bs Ms θ) (δ bu Mu ǫ) bu λ Mu ζ M are − − − s  (B2)  Mu,β (δ bu Mu ǫ)= au ǫ Ms,β ζ − −  ,  (as θ au λ) (δ bu Mu ǫ) au ǫ bu λ Mu Ms,β = − − − ( η bs Ms θ) (δ bu Mu ǫ) bu λ Mu ζ − − −  (B1) 

Mu,ββ (δ bu Mu ǫ)= du ǫ + Au,β ǫ,β (Mu,βδ,β + Ms,βζ,β + ǫ Du,β Mu Mu,β) Ms,ββ ζ − − −   (ds θ du λ + As,β θ,β Au,β λ,β + λ Du,β Mu Mu,β θ Ds,β Ms Ms,β) (δ bu Mu ǫ)  Ms,ββ = − − − − ,  (η bs Ms θ) (δ bu Mu ǫ) ζ λ bu Mu −  − − − [du ǫ + Au,β ǫ,β (Mu,βδ,β + Ms,βζ,β + ǫ Du,β Mu Mu,β)] λ bu Mu −  − (η bs Ms θ) (δ bu Mu ǫ) ζ λ bu Mu  − − −  (B3) 

Mu,γγ (δ bu Mu ǫ)= eu ǫ + Au,γ ǫ,γ (Mu,γδ,γ + Ms,γζ,γ + ǫ Du,γ Mu Mu,γ) Ms,γγ ζ − − −   (es θ eu λ + As,γ θ,γ Au,γ λ,γ + λ Du,γ Mu Mu,γ θ Ds,γ Ms Ms,γ) (δ bu Mu ǫ)  Ms,γγ = − − − −  (η bs Ms θ) (δ bu Mu ǫ) ζ λ bu Mu −  − − − [eu ǫ + Au,γ ǫ,γ (Mu,γδ,γ + Ms,γζ,γ + ǫ Du,γ Mu Mu,γ)] λ b M −  u u  − (η bs Ms θ) (δ bu Mu ǫ) ζ λ bu Mu  − − − (B4)  and

Mu,βγ (δ bu Mu ǫ)= fu ǫ + Au,β ǫ,γ (Mu,βδ,γ + Ms,βζ,γ + ǫ Du,β Mu Mu,γ) Ms,βγ ζ − − −   (fs θ fu λ + As,β θ,γ Au,β λ,γ + λ Du,β Mu Mu,γ θ Ds,β Ms Ms,γ) (δ bu Mu ǫ)  Ms,βγ = − − − − ,  (η bs Ms θ) (δ bu Mu ǫ) ζ λ bu Mu −  − − − [fu ǫ + Au,β ǫ,γ (Mu,βδ,γ + Ms,βζ,γ + ǫ Du,β Mu Mu,γ)] λ bu Mu −  − (η bs Ms θ) (δ bu Mu ǫ) ζ λ bu Mu  − − −  (B5) 

Λ where Nc 2 n−f (1 n−f ) n+f (1 n+f ) cf = 2 dp p − − − , Λ π 0 Ef Nc 2 Z − − (B8) af = 2 dp p [n f (1 n f )+ n+f (1 n+f )] , π 0 − − Z (B6)

Λ Nc 2 n−f (1 n−f )+ n+f (1 n+f ) bf = dp p − − , π2 T E2 Z0 f (B7)

12 Λ 2 Nc 2 2 M M,β 2 M,β − − df = 2 dp p 2 + p 4 [n f (1 n f )+ n+f (1 n+f )]+ π 0 ( E T E ! − − Z (B9) M M,β + 1+ [(1 2 n−f ) n−f,β + (1 2 n f ) n f,β] , T E2 − − + + )

Λ 2 Nc 2 2 M,γ ef = dp p p [n−f (1 n−f )+ n f (1 n f )] + π2 T E4 − + − + Z0 ( M M,γ + [(1 2 n−f ) n−f,γ + (1 2 n f ) n f,γ]+ (B10) T E2 − − + +

(1 2 n−f ) n−f,γ (1 2 n f ) n f,γ + − − − + + , E )

Λ Nc 2 M M,γ 2 M,βM,γ − − ff = 2 dp p 2 + p 4 [n f (1 n f )+ n+f (1 n+f )] + π 0 ( E T E − − Z (B11) M M,β + 1+ [(1 2 n−f ) n−f,γ + (1 2 n f ) n f,γ] , T E2 − − + + )

Af,β = af + bf Mf Mf,β , (B12) θ(Mu, Ms) = 4 G (1 F uB) Ms , (B22) − 1 λ(M , M ) = 4 K u M , (B23) Af,γ = cf + bf Mf Mf,γ , (B13) u s u

Λ N Ψf C = d + b M M , (B14) F = c dp p4 , (B24) f,ββ f f f f,ββ 1f π2 E3 Z0 f

Cf,γγ = ef + bf Mf Mf,γγ , (B15) 1 nf± = (B25) p2 M 2±µ q + f f C = f + b M M , (B16) 1+exp T f,βγ f f f f,βγ ( )

2 and u uu , s ss . K 2 ≡ h i ≡ h i B(Mu, Ms) = 4 G 2 u 2 K s , (B17) About the thermodynamic potential φ = − G − Ω β one has − Λ δ(Mu, Ms)=(1 F1u B) , (B18) Nc 2 2 − φ,β = dp p Ef Ψf + 2 Gs + π2 f u,d,s 0 =X Z ζ(M , M ) = 2 K u , (B19) u s + u (Mu mu)+ s (Ms ms)= − − N Λ c dp p2 E Ψ + K u2 s+ ǫ(Mu, Ms)= B Mu , (B20) 2 f f π 0 f=Xu,d,s Z s (Ms ms) η(Mu, Ms) = (1 4GF s) (1 F uB) + u (Mu mu)+ − , − 1 − 1 − − 2 2 2 (B21) 8K u F u , (B26) − 1

13 Λ Λ Nc 2 Nc 2 − − φ,γ = 2 dp p (n+f n f ) (B27) φ,βγγ = 2 dp p (n+f,βγ n f,βγ) π 0 − π 0 − f=Xu,d,s Z f=Xu,d,s Z (B32)

Λ Nc Λ 2 Nc 2 φ,ββ= dp p Ef (n+f,β + n−f,β) φ = dp p (n n− ) − π2 ,γγγ 2 +f,γγ f,γγ f=u,d,s Z0 π 0 − X f=Xu,d,s Z (B28) (B33)

Λ Nc 2 Λ − Nc 2 φ,βγ = 2 dp p (n+f,β n f,β) − π − φ,βββ= 2 dp p Ef (n+f,ββ + n f,ββ)+ f=u,d,s 0 − π X Z f=u,d,s Z0 (B29) X + (af + bf Mf Mf,β) Mf Mf,β Λ Nc 2 f=u,d,s φ = dp p (n n− ) X ,γγ 2 +f,γ f,γ (B34) π 0 − f=Xu,d,s Z (B30) Finally, by deﬁning Λ Nc 2 2 − H = F + b M , (B35) φ,ββγ = 2 dp p (n+f,ββ n f,ββ) f 1f f f π 0 − f=Xu,d,s Z (B31) the chiral susceptibilities are

∂Mu 1 4 GHs χu = χd = = − 2 2 2 (B36) ∂mu 1 4 G(Hu + Hs) + 4 Hu Hs(4 G 2 KGs K u ) + 2 KsHu − − − and

∂Ms 1 (4 G 2 Ks) Hu χs = == − −2 2 2 (B37) ∂ms 1 4 G(Hu + Hs) + 4 Hu Hs(4 G 2 KGs 2 K u ) + 2 KsHu − − −

[1] N. Ay, J. Jost, H. Van Le, and L. Schwachhofer. mannian geometric model. Phys. Rev. A, Information Geometry, A Series of Modern Sur- 20:1608–1613, Oct 1979. veys in Mathematics Volume 64, Springer. [7] George Ruppeiner. Riemannian geometry in [2] Shun ichi Amari. Information Geometry and thermodynamic ﬂuctuation theory. Rev. Mod. Its Applications, Applied Mathematical Sciences Phys., 67:605–659, 1995. [Erratum: Rev. Mod. Volume 194, Springer. Phys.68,313(1996)]. [3] M. Suzuki. Information geometry and statistical [8] George Ruppeiner, Anurag Sahay, Tapobrata manifold, arXiv:1410.3369. Sarkar, and Gautam Sengupta. Thermodynamic [4] R. C. Rao. Information and accuracy attain- Geometry, Phase Transitions, and the Widom able in the estimation of statistical parameters. Line. Phys. Rev., E86:052103, 2012. Bulletin of the Calcutta Mathematical Society, [9] Helge-Otmar May and Peter Mausbach. Rie- 37:81–91, 1945. mannian geometry study of vapor-liquid phase [5] F. Weinhold. Metric geometry of equilibrium equilibria and supercritical behavior of the thermodynamics. J. Chem. Phys., 63:2479, lennard-jones ﬂuid. Phys. Rev. E, 85:031201, 2484, 2488, 2496., 1975. Mar 2012. [6] George Ruppeiner. Thermodynamics: A rie- [10] Helge-Otmar May, Peter Mausbach, and George

14 Ruppeiner. Thermodynamic curvature for at- [26] A. Bazavov et al. The QCD Equation of State 6 tractive and repulsive intermolecular forces. to (µB ) from Lattice QCD. Phys. Rev., Phys. Rev. E, 88:032123, Sep 2013. D95(5):054504,O 2017. [11] Anshuman Dey, Pratim Roy, and Tapobrata [27] H. Janyszek. Riemannian geometry and sta- Sarkar. Information geometry, phase transi- bility of thermodynamical equilibrium systems. tions, and the Widom line: Magnetic and liquid Journal of Physics A: Mathematical and Gen- systems. Physica, A392:6341–6352, 2013. eral, 23(4):477–490, feb 1990. [12] Pankaj Chaturvedi, Anirban Das, and Gau- [28] M. R. Ubriaco. The role of curvature in quan- tam Sengupta. Thermodynamic Geometry and tum statistical mechanics. J. Phys.: Conf. Ser., Phase Transitions of Dyonic Charged AdS Black 766:012007, 2016. Holes. Eur. Phys. J., C77(2):110, 2017. [29] R. Ruppeiner, N. Dyjack, A. McAloon, and J. [13] Anurag Sahay and Rishabh Jha. Geometry Stoops. Solid-like features in dense vapors near of criticality, supercriticality and Hawking-Page the fluid critical point. The Journal of Chemical transitions in Gauss-Bonnet-AdS black holes. Physics, 146:224501, 2017. Phys. Rev., D96(12):126017, 2017. [30] Behrouz Mirza and Hosein Mohammadzadeh. [14] P. Castorina, M. Imbrosciano, and D. Lanteri. Ruppeiner Geometry of Anyon Gas. Phys. Rev., Thermodynamic Geometry of Strongly Interact- E78:021127, 2008. ing Matter. Phys. Rev., D98(9):096006, 2018. [31] M. R. Ubriaco. Stability and anyonic behavior of [15] Paolo Castorina, Mauro Imbrosciano, and systems with m-statistics. Physica A: Statistical Daniele Lanteri. Thermodynamic Geometry Mechanics and its Applications, 392(20):4868– and Deconfinement Temperature. Eur. Phys. J. 4873, 2013. Plus, 134(4):164, 2019. [32] Anurag Sahay, Tapobrata Sarkar, and Gautam [16] George Ruppeiner. Riemannian geometric ap- Sengupta. On the Thermodynamic Geometry proach to critical points: General theory. Phys. and Critical Phenomena of AdS Black Holes. Rev. E, 57:5135–5145, May 1998. JHEP, 07:082, 2010. [17] B. Widom. The critical point and scaling theory. [33] G. Ruppeiner, P. Mausbach, and H.-0. May. Physica, 73(1):107 – 118, 1974. Thermodynamic r-diagrams reveal solid-like [18] George Ruppeiner. Thermodynamic curvature fluid states. Physics Letters A, 379(7):646–649, from the critical point to the triple point. Phys. 2015. Rev. E, 86:021130, Aug 2012. [34] P. Zhuang, J. Hufner, and S. P. Klevansky. [19] H. Janyszek and R. Mrugala. Riemannian ge- Thermodynamics of a quark - meson plasma in ometry and the thermodynamics of model mag- the Nambu-Jona-Lasinio model. Nucl. Phys., netic systems. Phys. Rev. A, 39:6515–6523, Jun A576:525–552, 1994. 1989. [35] T. M. Schwarz, S. P. Klevansky, and G. Papp. [20] H. Janyszek and R. Mrugala. Riemannian and The Phase diagram and bulk thermodynamical finslerian geometry and fluctuations of ther- quantities in the NJL model at finite tempera- modynamic systems. Advances in Thermody- ture and density. Phys. Rev., C60:055205, 1999. namics, Vol. 3, Nonequilibrium Theory and Ex- [36] Michael Buballa. NJL model analysis of quark tremum Principles, pages 159–174, 1990. matter at large density. Phys. Rept., 407:205– [21] G. Ruppeiner. Thermodynamic curvature mea- 376, 2005. sures interactions. Am. J. Phys., 78:1170, 2010. [37] A. Barducci, R. Casalbuoni, G. Pettini, and L. [22] Michele Floris. Hadron yields and the phase Ravagli. A NJL-based study of the QCD critical diagram of strongly interacting matter. Nucl. line. Phys. Rev., D72:056002, 2005. Phys., A931:103–112, 2014. [38] Yue Zhao, Lei Chang, Wei Yuan, and Yu-xin [23] Sabita Das. Identified particle production and Liu. Chiral susceptibility and chiral phase tran- freeze-out properties in heavy-ion collisions at sition in Nambu-Jona-Lasinio model. Eur. Phys. RHIC Beam Energy Scan program. 2014. [EPJ J., C56:483–492, 2008. Web Conf.90,08007(2015)]. [39] Aharon Casher. Chiral Symmetry Breaking [24] L. Adamczyk et al. Bulk Properties of the in Quark Confining Theories. Phys. Lett., Medium Produced in Relativistic Heavy-Ion 83B:395–398, 1979. Collisions from the Beam Energy Scan Program. [40] Tom Banks and A. Casher. Chiral Symmetry Phys. Rev., C96(4):044904, 2017. Breaking in Confining Theories. Nucl. Phys., [25] Patrick Steinbrecher. The QCD crossover at B169:103–125, 1980. zero and non-zero baryon densities from Lattice [41] P. Castorina and P. Cea. he role of chirality in QCD. Nucl. Phys., A982:847–850, 2019. demarcating confinement from liberation. Phys.

15 Lett., 112B:151–153, 1982. [43] H. T. Ding, P. Hegde, F. Karsch, A. Lahiri, S. T. [42] S. Digal, E. Laermann, and H. Satz. Decon- Li, S. Mukherjee, and P. Petreczky. Chiral phase finement through chiral symmetry restoration in transition of (2+1)-flavor QCD. Nucl. Phys., two flavor QCD. Eur. Phys. J., C18:583–586, A982:211–214, 2019. 2001.