Stress tensor distribution around static quarks in hot medium
Ryosuke Yanagihara (Osaka University) For FlowQCD collaboration : Takumi Iritani, Masakiyo Kitazawa, Masayuki Asakawa, Tetsuo Hatsuda
FLQCD 2019 @ YITP (2019/04/18)
Stress tensor distribution around static quarks in hot medium - - PowerPoint PPT Presentation
Stress tensor distribution around static quarks in hot medium Ryosuke Yanagihara (Osaka University) For FlowQCD collaboration : Takumi Iritani, Masakiyo Kitazawa, Masayuki Asakawa, Tetsuo Hatsuda FLQCD 2019 @ YITP (2019/04/18) Confined vs.
Stress tensor distribution around static quarks in hot medium
Ryosuke Yanagihara (Osaka University) For FlowQCD collaboration : Takumi Iritani, Masakiyo Kitazawa, Masayuki Asakawa, Tetsuo Hatsuda
FLQCD 2019 @ YITP (2019/04/18)Confined vs. Deconfined
π
FLQCD 2019 @ YITP (2019/04/18)Critical temperature π
π
Confined Deconfined
gluon quark
1
ζΈ©εΊ¦ π
FLQCD 2019 @ YITP (2019/04/18)ιιγθΎΌγηΈ
π
π Burkert et al.,
Nature 557 (2018) 396. Shanahan et al., PRL 122 (2019) no7, 072003. Exp. Th. Kumano et al., PRD 97 (2018) 014020.
Confined vs. Deconfined
Confined
gluon quark
Pressure distribution inside Hadrons
1
Pressure distribution inside hadrons vs. Our study
π ΰ΄€ π
FLQCD 2019 @ YITP (2019/04/18)2
Burkert et al., Nature 557 (2018) 396. Shanahan et al., PRL 122 (2019) no7, 072003.Exp. Th.
Kumano et al., PRD 97 (2018) 014020.Pressure distribution inside Hadrons Our study
Flux tube
QED QCD
οΌ Flux tube, squeezed one-dimensionally οΌ Confinement potential οΌ Electric field spreads all over the space οΌ Coulomb potential
π ΰ΄€ π
FLQCD 2019 @ YITP (2019/04/18)3
Flux tube
QED QCD
οΌ Flux tube, squeezed one-dimensionally οΌ Confinement potential οΌ Electric field spreads all over the space οΌ Coulomb potential
π ΰ΄€ π Maxwell stress
Local interaction
FLQCD 2019 @ YITP (2019/04/18)3
Cardoso et al., PRD86 (2013) 054501. Cea et al., PRD88 (2012) 054504.
Action density Color electric field
π¦ π¨ π§ mid π§π¨ πππ¨ ππ ππ
A lot of previous studies
FLQCD 2019 @ YITP (2019/04/18)4
Cardoso et al., PRD86 (2013) 054501. Cea et al., PRD88 (2012) 054504.
Action density Color electric field
π¦ π¨ π§ mid π§π¨ πππ¨ ππ ππ
A lot of previous studies
FLQCD 2019 @ YITP (2019/04/18)More direct physical quantity : Stress tensor !!
4
Energy momentum tensor (EMT)
οΌ Stress is force per unit area π
π = πππππ ; πππ = βπππ
π
ππ =
π00 π01 π02 π03 π
10
π
11
π
12
π
13
π20 π30 π21 π31 π22 π23 π32 π33
Energy density
Momentum density Pressure
Stress tensor
Landau and Lifshitz
FLQCD 2019 @ YITP (2019/04/18)5
οΌ Stress is force per unit area π
π = πππππ ; πππ = βπππ
π
ππ =
π00 π01 π02 π03 π
10
π
11
π
12
π
13
π20 π30 π21 π31 π22 π23 π32 π33
Energy density
Momentum density Pressure
Landau and Lifshitz
FLQCD 2019 @ YITP (2019/04/18)rubber
Energy momentum tensor (EMT) Stress tensor
5
Maxwell stress
π
ππ = π0 πΉππΉ π β πππ2 πΉ2 + 1 π0 πΆππΆ
π β πππ2 πΆ2 πΉ
οΌ Stress tensor
π
ππππ (π) = ππππ (π)(π, π = 1,2,3 ; π = 1,2,3) οΌ Perpendicular plane: ππ < 0 οΌ Parallel plane: ππ > 0
Length of arrows= ππ Τ¦ π Τ¦ π
FLQCD 2019 @ YITP (2019/04/18)6
Measurement on the lattice
β Prepare π ΰ΄€ π on the lattice β‘Measure EMT around π ΰ΄€ π
To do
7
Measurement on the lattice
β Prepare π ΰ΄€ π on the lattice β‘Measure EMT around π ΰ΄€ π
To do
7
π π = β lim
πββ1 π logβ¨π π, π β©
Ground state potential
Wilson Loop
Confinement potential β¨π π, π β© = π·0exp βπ π π + π·1exp βπ
1 π π + β―οΌ quenched SU(3) Yang-Mills οΌ πΎ = 6.600 (π = 0.038 fm)
EMT defined via gradient flow
π ππ π’, π¦ = 1 π½π π’ πππ π’, π¦ + πππ 4π½πΉ(π’) πΉ π’, π¦ β πΉ π’, π¦ + π(π’) Suzuki (2013)Gradient flow
ππΆ π π’, π¦ ππ’ = βπ0 2 ππ[πΆ] ππΆ π(π’, π¦)Flow eq.
L α· uscher (2010)πΆ
π : smeared fieldIritani et al. (2018) Entropy density vs. temperature οΌ 2-loop coefficient is now available !
Harlander et al. (2018)Measurement on the lattice
β Prepare π ΰ΄€ π on the lattice β‘Measure EMT around π ΰ΄€ π
To do
7
Set up
FLQCD 2019 @ YITP (2019/04/18)8
οΌ Quenched SU(3) Yang-Mills gauge theory οΌ Wilson gauge action οΌ Clover operator οΌ Continuum limit οΌ APE smearing for spatial links οΌ Multihit improvement in temporal links οΌ Simulation using BlueGene/Q @ KEK πΈ Lattice spacing Lattice size # of statistics 6.304 0.057 fm 484 140 6.465 0.046 fm 484 440 6.513 0.043 fm 484 600 6.600 0.038 fm 484 1500 6.819 0.029 fm 644 1000 0.92 fm 0.69 fm 0.46 fm
π§ mid-plane π§π¨-plane
π
π¨ π¦ π π π ΰ΄€ π
FLQCD 2019 @ YITP (2019/04/18)A lattice study of stress distribution around π ΰ΄€ π in vacuum
FlowQCD, PLB 789 (2019) 210.
Stress distribution in terms of local interaction
9
Stress distribution around π ΰ΄€ π
π¦ π¨ π§ mid π§π¨ ποΌ π = 0.029 fm οΌ π’/π2 = 2.0 οΌ π = 0.69 fm οΌ Length of arrows = ππ
FLQCD 2019 @ YITP (2019/04/18)οΌ Gauge invariant οΌ Local interaction οΌ squeezed
FlowQCD, PLB 789 (2019) 210.
10
Stress distribution around π ΰ΄€ π : Cylindrical coordinates
π
ππ =
π
44
π
π¨π¨
π
π π
πππ Diagonalized EMT
(Cylindrical / Parity symmetry)
π π
Degeneracy (Maxwell Theory) π
44 = π π¨π¨ = π π π = |πππ|
ππ¨ ππ ππ
π ΰ΄€ π
11
Stress distribution around π ΰ΄€ π
FLQCD 2019 @ YITP (2019/04/18)ππ¨ ππ ππ
π¦ π¨ π§ mid π§π¨ ποΌ π
44 β π π¨π¨, π π π β π ππ (Degeneracy)οΌ π
44 β π π π (Separation)οΌ Οπ π
ππ β 0 (Trace anomalyβ 0)Properties in non-Abelian theory
( Note : after double limit )
FlowQCD, PLB 789 (2019) 210.
12
EMT and confinement potential
FLQCD 2019 @ YITP (2019/04/18)confinement potential π β βπΊstress β ΰΆ± mid π
π¨π¨ π ΰ΄€ π π2π¦From EMT π π = π + ππ + Ξ€ π π γγγ«ζ°εΌγε ₯εγγΎγγ πΊpot β β ππ(π) ππ
13
EMT and confinement potential
FLQCD 2019 @ YITP (2019/04/18)confinement potential π β βπΊstress β ΰΆ± mid π
π¨π¨ π ΰ΄€ π π2π¦From EMT π π = π + ππ + Ξ€ π π γγγ«ζ°εΌγε ₯εγγΎγγ πΊpot β β ππ(π) ππ
Good agreement !!
13
Toward analysis at nonzero temperature
π
FLQCD 2019 @ YITP (2019/04/18)Stress distribution around π ΰ΄€ π /π at nonzero temperature
Critical temperature π
π
14
Measurement on the lattice
β Prepare π ΰ΄€ π on the lattice β‘Measure EMT around π / ΰ΄€ π
To do
FLQCD 2019 @ YITP (2019/04/18)15
πβπΊ π /π = 1 3 Tr Ξ©β Τ¦ π¦ Ξ© ( Τ¦ π§)
Color singlet free energy (We use Coulomb gauge fixing)
Polyakov Loop
Free energy
οΌ quenched SU(3) Yang-Mills οΌ πΎ = 6.600 (π = 0.038 fm)
π¦4 Τ¦ π¦
π ΰ΄€ π Measurement on the lattice
β Prepare π ΰ΄€ π on the lattice β‘Measure EMT around π / ΰ΄€ π
To do
15
EMT defined via gradient flow
π ππ π’, π¦ = 1 π½π π’ πππ π’, π¦ + πππ 4π½πΉ(π’) πΉ π’, π¦ β πΉ π’, π¦ + π(π’) Suzuki (2013)Gradient flow
ππΆ π π’, π¦ ππ’ = βπ0 2 ππ[πΆ] ππΆ π(π’, π¦)Flow eq.
L α· uscher (2010)πΆ
π : smeared fieldIritani et al. (2018) Entropy density vs. temperature οΌ 2-loop coefficient is now available !
Harlander et al. (2018)Measurement on the lattice
β Prepare π ΰ΄€ π on the lattice β‘Measure EMT around π / ΰ΄€ π
To do
15
Set up (quarkβanti-quark, single quark)
οΌ Quenched SU(3) Yang-Mills gauge theory οΌ Wilson gauge action οΌ Clover operator οΌ Fixed π, π’ οΌ Multihit improvement in temporal links οΌ Simulation using OCTOPUS, Reedbush
FLQCD 2019 @ YITP (2019/04/18)πΈ Lattice spacing Spatial size Temporal size πΌ/πΌπ # of statistics 6.600 0.038 fm 483 12 1.44 640
π ΰ΄€ π
0.69 fm 0.46 fm
16
Maxwell stress (revisit)
π
ππ = π0 πΉππΉ π β πππ2 πΉ2 + 1 π0 πΆππΆ
π β πππ2 πΆ2 πΉ
οΌ Stress tensor
π
ππππ (π) = ππππ (π)(π, π = 1,2,3 ; π = 1,2,3) οΌ Perpendicular plane: ππ < 0 οΌ Parallel plane: ππ > 0
Length of arrows= ππ Τ¦ π Τ¦ π
FLQCD 2019 @ YITP (2019/04/18)6β
Preliminary
οΌ singlet οΌ π = 0.038 fm (fixed) οΌ π’/π2 = 2.0 (fixed) οΌ π = 0.69 fm οΌ Length of arrows = ππ
Stress distribution around π ΰ΄€ π
17
π Preliminary Critical temperature π
π
Stress distribution around π ΰ΄€ π
18
Stress distribution around π ΰ΄€ π : Cylindrical coordinates
π
ππ =
π
44
π
π¨π¨
π
π π
πππ Diagonalized EMT
(Cylindrical / Parity symmetry)
π π
Degeneracy (Maxwell Theory) π
44 = π π¨π¨ = π π π = |πππ|
ππ¨ ππ ππ
π ΰ΄€ π
19
ππ¨ ππ ππ
Preliminary
( Note : singlet )
Stress distribution around π ΰ΄€ π
20
ππ¨ ππ ππ Preliminary Preliminary
οΌ π
π π β π ππ (Degeneracy)οΌ π
44 β π π π (Separation)οΌ Οπ π
ππ β 0 (Trace anomalyβ 0)οΌ Damping β Debye mass οΌ ππΊ/ππ β Χ¬ π
π¨π¨ π2π¦Stress distribution around π ΰ΄€ π
20
Preliminary ππ¨ ππ ππ
Stress distribution around π ΰ΄€ π
οΌ π
π π β π ππ (Degeneracy)οΌ π
44 β π π π (Separation)οΌ Οπ π
ππ β 0 (Trace anomalyβ 0)οΌ Damping β Debye mass οΌ ππΊ/ππ β Χ¬ π
π¨π¨ π2π¦20
Preliminary
πΊstress β 0.30 GeV/fm πΊFree β 0.41 GeV/fm
ππ¨ ππ ππ
Stress distribution around π ΰ΄€ π
οΌ π
π π β π ππ (Degeneracy)οΌ π
44 β π π π (Separation)οΌ Οπ π
ππ β 0 (Trace anomalyβ 0)οΌ Damping β Debye mass οΌ ππΊ/ππ β Χ¬ π
π¨π¨ π2π¦20
Stress distribution around π : Spherical coordinates
π π π, π π, π
π
ππ =
π
44
π
π π
ππ’π’ ππ’π’
π π
Diagonalized EMT
(Spherical symmetry)
Degeneracy (Maxwell Theory) π
44 = π π π = |ππ’π’|
21
(π’: π, π)
Stress distribution around π
π π, π π, π
οΌ Overlap : π β² 2π’ οΌ Separation οΌ Around originβΌ π½π‘(π )/π 4 οΌ DampingβΌ πβππΈπ οΌ Conservation lawοΌ ππ π 2π
π π = π π π’π’Preliminary 22
Pressure distribution inside hadrons vs. Our study
FLQCD 2019 @ YITP (2019/04/18)23
Burkert et al., Nature 557 (2018) 396. Shanahan et al., PRL 122 (2019) no7, 072003.Exp. Th.
Kumano et al., PRD 97 (2018) 014020.Pressure distribution inside Hadrons Our study
24
Summary and Outlook
Summary οΌ We first measure stress distribution around π ΰ΄€ π /π at zero/nonzero temperature on the lattice Outlook οΌ π, π’ β 0 (double limit) οΌ Temperature dependence οΌ ApplicationοΌπ π γπ π π γexcited stateγhadron (full QCD)β¦
Single quark QuarkβAnti-Quark Preliminary Preliminary Preliminary
Back up
FLQCD 2019 @ YITP (2019/04/18)Flow time dependence (single quark system)
Interference b/w singlet and octet state
FLQCD 2019 @ YITP (2019/04/18)Preliminary Preliminary