INDICE INTRODUCCIÓN PEQUEÑA BIOGRAFÍA SU INFANCIA SU EDUCACIÓN LOS 12 TRABAJOS.
SU X - conference-indico (Indico)LinkedClusterExpansionfor SU(3) Spin ModelwithaSignProblem Quang...
1
Linked Cluster Expansion for SU(3) Spin Model with a Sign Problem Quang Pham, Jangho Kim, Owe Philipsen, Jonas Scheunert Institute for Theoretical Physics, Goethe University Frankfurt, Frankfurt am Main Introduction T μ chirally restored quark matter vacuum nuclear matter early universe RHIC FAIR ? ? ? Recently, a lot of attention has been paid to whether a critical point exists on the phase diagram of strongly interacting matter (quark matter). The main tool for extracting quantitative results is based on Monte Carlo algorithms. Unfortunately, at nite baryon densities all Monte Carlo algorithms suer from the sign problem. We revisit a well-known method in statistical mechanics called Series Expansion Method, which provides an alternative way to overcome the sign problem. An SU(3) spin model is used as a rst step to deal with QCD at nite chemical potential. SU(3) spin model The SU(3) spin model is a three dimensional eecti- ve theory of QCD at strong coupling limit and large quark mass limit. The partition function of the model reads Z = Z D[U 0 ] exp τ X h ~ x,~ yi (L ~ x L * ~ y + L * ~ x L ~ y ) +κ X ~ x (e aμ L ~ x + e -aμ L * ~ x ) ! , where L ~ x ’s are Polyakov loops, κ ∼ m -2 q . The rst term presents the nearest neighbour interaction of Polyakov loops, while the second term couples the system to quarks at nite chemical potential μ. A particular observable that was used in our compu- tation is the susceptibility of equilibrium Polyakov loop L, which diverges near critical points: L = 1 V X ~ x (L ~ x + L * ~ x ), χ L = 1 V ∂ 2 log Z[J] ∂J 2 Conclusions • Analytic calculations for the SU(3) spin mo- del via LCE are done for the rst time. • Our results agree with those from simulati- on with ux representation algorithm within errors. • To include the parameter μ it is needed to compute higher orders. • Our next goal is to apply the series expansion method to improved eective theories. Linked Cluster Expansion (LCE) Similar to the generating functional for connected diagrams in QFT, we dene the grand canonical potential as W [J, ν]= log Z[J, ν] where Z[J, ν] is the partition function coupled to an external source J. A linked cluster expansion (LCE) of the grand canonical potential is then dened as the Taylor expansion with respect to the coupling v(x,y) around the free theory, W [J, ν]= exp ν( ~ x, ~ y) δ δ ˆ ν( ~ x, ~ y) W [J,ˆ ν] ˆ ν=0 There is a very simple expression of W through graph representation, In each term of the graph representation, the number of bonds corresponds to the order of the coupling ν. In additions, to complete the computation one also needs to calculate integrals at each vertices of graphs, e.g. I mn = Z dU 0 L m L *n Series Analysis The Polyakov loop susceptibility computed up to order K = 13 is a series of the interaction coupling τ while the coecients are functions of κ and μ, χ L = K X n=0 c n (κ, μ)τ n . These polynomials are obviously not divergent near critical points. Thus, the (generalised) Pade’ approxi- mant is needed for extracting information about cri- tical points and critical exponents. Example: for one coupling series (κ = 0), the Polyakov loop suscepti- bility near τ c behaves as χ L ∼ 1 (τ - τ c ) γ → d dτ log χ L ∼ -γ τ - τ c ∼ [X/Y ]= a 0 + a 1 τ + ··· + a X τ X 1 + b 1 τ + ··· + b Y τ Y References [1] F. Karsch, H.W. Wyld, Complex Langevin simula- tion of the SU(3) spin model with nonzero chemical potential, Phys. Rev. Lett.55(1985) 2242. [2] Y. D. Mercado, C. Gattringer, Monte Carlo simu- lation of the SU(3) spin model with chemical po- tential in a ux representation, Nucl.Phys. B862 (2012) 737-750 [3] G. Aarts, F. A. James, Complex Langevin dyna- mics in the SU(3) spin model at nonzero chemical potential revisited, JHEP 1201 (2012) 118 [4] M. Wortis, in Phase Transitions and Critical Phe- nomena, vol. 3, C. Domb and M. S. Greeneds. (Academic Press, London, 1974), p. 113. [5] G. A. Baker, Jr., Application of the Padé approxi- mant method to the investigation of some magne- tic properties of the Ising Model, Phys. Rev. 124 768-74 (1961) [6] J. Kim, Q. Pham, O. Philipsen, J. Scheunert, to appear Results Estimation of τ c and critical exponents for SU(2), SU(3). γ SU(3) c ≈ 0.712 which is unknown universality class ∼ rst order phase transition 0.135 0.14 0.145 0.15 0.155 0.16 0.165 0.17 0.175 0.18 2 4 6 8 10 12 τ c n Pade result Simulation result 0 0.5 1 1.5 2 SU2 SU3 γ c SU(2) SU(3) Ising Unbiased estimations of critical points (κ c , τ c ) at μ = 0 from our series τ κ (κ c , τ c )
Transcript of SU X - conference-indico (Indico)LinkedClusterExpansionfor SU(3) Spin ModelwithaSignProblem Quang...
-
LinkedClusterExpansionforSU(3)SpinModelwithaSignProblem
Quang Pham, Jangho Kim, Owe Philipsen, Jonas ScheunertInstitute for Theoretical Physics, Goethe University Frankfurt, Frankfurt am Main
Introduction
T
μ
chirally restored quark matter
vacuumnuclearmatter
earl
y un
iver
se
RH
IC
FAIR ? ?
?
Recently, a lot of attention has been paid to whethera critical point exists on the phase diagram ofstrongly interacting matter (quark matter). Themain tool for extracting quantitative results is basedon Monte Carlo algorithms. Unfortunately, at �nitebaryon densities all Monte Carlo algorithms su�erfrom the sign problem.
We revisit a well-known method in statisticalmechanics called Series Expansion Method, whichprovides an alternative way to overcome the signproblem.
An SU(3) spin model is used as a �rst step to dealwith QCD at �nite chemical potential.
SU(3) spin modelThe SU(3) spin model is a three dimensional e�ecti-ve theory of QCD at strong coupling limit and largequark mass limit. The partition function of the modelreads
Z =
∫D[U0] exp
τ ∑〈~x,~y〉
(L~xL∗~y + L
∗~xL~y)
+κ∑~x
(eaµL~x + e−aµL∗~x)
),
where L~x’s are Polyakov loops, κ ∼ m−2q . The �rstterm presents the nearest neighbour interaction ofPolyakov loops, while the second term couples thesystem to quarks at �nite chemical potential µ.
A particular observable that was used in our compu-tation is the susceptibility of equilibrium Polyakovloop L, which diverges near critical points:
L =1V
∑~x
(L~x + L∗~x),
χL =1V
∂2 logZ[J]∂J2
Conclusions• Analytic calculations for the SU(3) spin mo-
del via LCE are done for the �rst time.
• Our results agree with those from simulati-on with �ux representation algorithm withinerrors.
• To include the parameter µ it is needed tocompute higher orders.
• Our next goal is to apply the series expansionmethod to improved e�ective theories.
Linked Cluster Expansion (LCE)Similar to the generating functional for connected diagrams in QFT, we de�ne the grand canonical potentialas
W[J,ν] = logZ[J,ν]
where Z[J,ν] is the partition function coupled to an external source J. A linked cluster expansion (LCE) ofthe grand canonical potential is then de�ned as the Taylor expansion with respect to the coupling v(x,y)around the free theory,
W[J,ν] =(exp
(ν(~x,~y)
δ
δν̂(~x,~y)
)W[J, ν̂]
) ∣∣∣ν̂=0
There is a very simple expression ofW through graph representation,
In each term of the graph representation, the number of bonds corresponds to the order of the coupling ν. Inadditions, to complete the computation one also needs to calculate integrals at each vertices of graphs, e.g.
Imn =
∫dU0LmL∗n
Series AnalysisThe Polyakov loop susceptibility computed up toorder K = 13 is a series of the interaction couplingτ while the coe�cients are functions of κ and µ,
χL =
K∑n=0
cn(κ,µ)τn.
These polynomials are obviously not divergent nearcritical points. Thus, the (generalised) Pade’ approxi-mant is needed for extracting information about cri-tical points and critical exponents. Example: for onecoupling series (κ = 0), the Polyakov loop suscepti-bility near τc behaves as
χL ∼1
(τ− τc)γ→ d
dτlog χL ∼
−γ
τ− τc
∼ [X/Y] =a0 + a1τ+ · · ·+ aXτX
1+ b1τ+ · · ·+ bYτY
References[1] F. Karsch, H.W. Wyld, Complex Langevin simula-
tion of the SU(3) spinmodel with nonzero chemicalpotential, Phys. Rev. Lett.55(1985) 2242.
[2] Y. D. Mercado, C. Gattringer, Monte Carlo simu-lation of the SU(3) spin model with chemical po-tential in a �ux representation, Nucl.Phys. B862(2012) 737-750
[3] G. Aarts, F. A. James, Complex Langevin dyna-mics in the SU(3) spin model at nonzero chemicalpotential revisited, JHEP 1201 (2012) 118
[4] M. Wortis, in Phase Transitions and Critical Phe-nomena, vol. 3, C. Domb and M. S. Greeneds.(Academic Press, London, 1974), p. 113.
[5] G. A. Baker, Jr., Application of the Padé approxi-mant method to the investigation of some magne-tic properties of the Ising Model, Phys. Rev. 124768-74 (1961)
[6] J. Kim, Q. Pham, O. Philipsen, J. Scheunert, toappear
ResultsEstimation of τc and critical exponents for SU(2), SU(3). γSU(3)c ≈ 0.712 which is unknown universalityclass ∼ �rst order phase transition
0.135
0.14
0.145
0.15
0.155
0.16
0.165
0.17
0.175
0.18
2 4 6 8 10 12
τc
n
Pade result
Simulation result
0
0.5
1
1.5
2
SU2 SU3
γc
SU(2)SU(3)Ising
Unbiased estimations of critical points (κc, τc) at µ = 0 from our series
τ
κ
(κc, τc)
9 10 11 12 13 14order
0.10
0.12
0.14
0.16
0.18
c
order=10order=11order=12order=13
9 10 11 12 13 14order
0.0000
0.0025
0.0050
0.0075
0.0100
0.0125
0.0150
0.0175
0.0200c
order=10order=11order=12order=13
