Quarkonium dissociation by anisotropy in a strongly coupled plasma - - PowerPoint PPT Presentation

quarkonium dissociation by anisotropy in a strongly
SMART_READER_LITE
LIVE PREVIEW

Quarkonium dissociation by anisotropy in a strongly coupled plasma - - PowerPoint PPT Presentation

Quarkonium dissociation by anisotropy in a strongly coupled plasma Mariano Chernicoff DAMTP, University of Cambridge Based on 1208.2672 (and 1202.3696) in collaboration with D. Fernandez, D. Mateos and D. Trancanelli


slide-1
SLIDE 1

Quarkonium dissociation by anisotropy in a strongly coupled plasma

Mariano ¡Chernicoff ¡ DAMTP, ¡University ¡of ¡Cambridge ¡

Based on 1208.2672 (and 1202.3696) in collaboration with

  • D. Fernandez, D. Mateos and D. Trancanelli
slide-2
SLIDE 2

Plan for the talk

  • Motivation
  • Quarkonium physics
  • The AdS/CFT toolkit
  • Conclusions
slide-3
SLIDE 3

Motivation

For the creation of the quark gluon plasma 400 nucleons go in 8000 hadrons are produced Energy of CM 200 GeV 1-10 GeV per hadron

slide-4
SLIDE 4

Motivation

Collision Far from equilibrium Anisotropic hydrodynamics Isotropic hydrodynamics

τout τiso

P? 6= Pk P? = Pk

longitudinal ¡direc>on ¡ transverse ¡direc>on ¡

z ≡ x, y ≡

[Romantschke ¡et ¡al; ¡Mrowczynski ¡et ¡al.] ¡

τiso . 1fm The time evolution of the quark gluon plasma

slide-5
SLIDE 5

Motivation

Things we know about the quark gluon plasma After a very short period of time the system is in thermal equilibrium

∼ 200 − 400MeV

slide-6
SLIDE 6

Motivation

Things we know about the quark gluon plasma After a very short period of time the system is in thermal equilibrium

∼ 200 − 400MeV

It is a strongly coupled system ( Perturbative methods are inapplicable and lattice QCD is not well suited to study real-time dynamics. ideal fluid + small )

η/s

slide-7
SLIDE 7

Motivation

Things we know about the quark gluon plasma After a very short period of time the system is in thermal equilibrium

∼ 200 − 400MeV

It is a strongly coupled system ( Perturbative methods are inapplicable and lattice QCD is not well suited to study real-time dynamics. ideal fluid + small )

η/s

Some observables are sensitive to the presence of an anisotropy

  • Quarkonium physics ( )

[Dumitru ¡et ¡al.; ¡Philipsen ¡et ¡al.] ¡ ¡

J/Ψ

  • Momentum broadening

[Dumitru ¡et ¡al.; ¡Mehtar-­‑Tani; ¡Romantschke ¡et ¡al.] ¡

slide-8
SLIDE 8

Motivation

Quarkonium refers to charm-anticharm mesons ( , , , ... )

J/Ψ Ψ0 χc

and bottom-antibottom mesons ( , , ... )

Υ Υ0

1s ¡state ¡

About quarkonium in heavy ion collisions:

slide-9
SLIDE 9

Motivation

About quarkonium in heavy ion collisions: Quarkonium refers to charm-anticharm mesons ( , , , ... )

J/Ψ Ψ0 χc

and bottom-antibottom mesons ( , , ... )

Υ Υ0

dissociation temperature that is higher than the deconfinement temperature ( lattice predicts: ). mesons survive as bound states in a hot medium up to some

J/Ψ Tc Td Td(J/Ψ) ' 2Tc

slide-10
SLIDE 10

RHIC data: - suppression in nucleus-nucleus collisions when compared to proton-proton collisions. dissociation temperature that is higher than the deconfinement temperature ( lattice predicts: ).

Motivation

About quarkonium in heavy ion collisions: Quarkonium refers to charm-anticharm mesons ( , , , ... )

J/Ψ Ψ0 χc

and bottom-antibottom mesons ( , , ... )

Υ Υ0

mesons survive as bound states in a hot medium up to some

J/Ψ J/Ψ Tc Td Td(J/Ψ) ' 2Tc

mesons are screened in the quark gluon plasma

J/Ψ

slide-11
SLIDE 11

Motivation

They might be moving with significant transverse momentum through the hot medium, what is the effect of such “wind”? It is important to understand how they are screened by the QGP What is the effect of the anisotropy? We will use the AdS/CFT correspondence to address this questions

slide-12
SLIDE 12

=

T = 0

λ ≡ L4 l4

s

= g2

YMNc

AdS/CFT toolkit

u = 0

=

N = 4 SYM AdS5 × S5

Type IIB on

z xy u = ∞ ds2 = L2 u2 ⇥ − dt2 + d~ x2 + du2⇤

  • n Minkowski 3+1
slide-13
SLIDE 13

AdS/CFT toolkit

z xy u = 0 u = uh

=

T 6= 0

=

at finite temperature Schwarzschild AdS black hole

N = 4 SYM

ds2 = L2 u2 ⇥ − f(u)dt2 + d~ x2 + du2 f(u) ⇤ Note: we will refer to this metric as the isotropic metric f(u) = 1 − u4 u4

h

with

slide-14
SLIDE 14

External quark

=

Fundamental string

AdS/CFT toolkit

u = 0 u = uh

=

Quark ¡

T 6= 0

  • The string endpoint represents the quark, while the rest of the string

codifies the profile of the gluonic field

  • A fundamental string extending from the boundary at to the

horizon corresponds to an infinitely massive quark.

u = 0

slide-15
SLIDE 15

meson (bound state) = U-shaped string

AdS/CFT toolkit

u = 0 u = uh

=

T 6= 0

` `

slide-16
SLIDE 16

meson (bound state) moving at constant velocity

=

U-shaped string moving at constant velocity

AdS/CFT toolkit

u = 0 u = uh

=

T 6= 0

slide-17
SLIDE 17

meson (bound state) moving at constant velocity

=

U-shaped string moving at constant velocity

AdS/CFT toolkit

u = 0 u = uh

=

T 6= 0

But we are interested in studying an anisotropic strongly coupled plasma How can we do that using the AdS/CFT correspondence?

slide-18
SLIDE 18

The gauge theory that we will consider is a deformation of N = 4 SYM

AdS/CFT toolkit

S = SN =4 + Z ✓(~ x)TrF ∧ F with

T 6= 0

dimensions ¡of ¡energy ¡

✓(~ x) = 2⇡nD7z

slide-19
SLIDE 19

z xy u = 0 u = uh

=

T 6= 0

S = SN =4 + Z ✓(~ x)TrF ∧ F The gauge theory that we will consider is a deformation of .

N = 4 SYM

AdS/CFT toolkit

with

χ = az

τ = θ 2π + 4πi g2

YM

= χ + ie−φ ✓(~ x) = 2⇡nD7z The field dual to the -term is the axion

θ

The axion is magnetically sourced by D7-branes

slide-20
SLIDE 20

z xy u = 0 u = uh

=

T 6= 0

S = SN =4 + Z ✓(~ x)TrF ∧ F The gauge theory that we will consider is a deformation of N = 4 SYM

AdS/CFT toolkit

with t x y z u S5 Nc D3 nD7 D7 x x x x x x x x ✓(~ x) = 2⇡nD7z

χ = az

The field dual to the -term is the axion

θ

slide-21
SLIDE 21

z xy u = 0 u = uh

=

T 6= 0

S = SN =4 + Z ✓(~ x)TrF ∧ F The gauge theory that we will consider is a deformation of N = 4 SYM

AdS/CFT toolkit

with

χ = az

a = λnD7 4πNc ✓(~ x) = 2⇡nD7z and

Density ¡of ¡D7-­‑branes ¡ distributed ¡along ¡z ¡ is ¡a ¡measure ¡of ¡the ¡anisotropy ¡

The field dual to the -term is the axion

θ

slide-22
SLIDE 22

z xy u = 0 u = uh

=

T 6= 0

S = SN =4 + Z ✓(~ x)TrF ∧ F The gauge theory that we will consider is a deformation of N = 4 SYM

AdS/CFT toolkit

with

χ = az

✓(~ x) = 2⇡nD7z Finally, putting all these ingredients together, and solving the eom... The field dual to the -term is the axion

θ

slide-23
SLIDE 23

ds2 = L2 u2 h − F(u)B(u)dt2 + dx2 + dy2 + H(u)dz2 + du2 F(u) i

AdS/CFT toolkit

The anisotropic metric is

[Mateos ¡and ¡Trancanelli] ¡

0. 0.2 0.4 0.6 0.8 1.

  • 2
  • 1

1 2 3 4 5 6

  • 2
  • 1

1 2 3 4 5 6 0. 0.2 0.4 0.6 0.8 1.

  • 0.2

0.2 0.4 0.6 0.8 1. 1.2

  • 0.2

0.2 0.4 0.6 0.8 1. 1.2

H H F F B B φ φ u/uh u/uh

χ(z) = az φ ≡ φ(u)

and

slide-24
SLIDE 24

AdS/CFT toolkit

The anisotropic metric is

[Mateos ¡and ¡Trancanelli] ¡

RG flow between AdS (UV) and Lifshitz type (IR) ds2 = L2 u2 h − F(u)B(u)dt2 + dx2 + dy2 + H(u)dz2 + du2 F(u) i The entropy density interpolates between

T a s ⇠ T 3 T ⌧ a s ⇠ a1/3T 8/3

and Regular on and outside the horizon There is an analytical expression for the near-boundary behavior

  • f metric functions O(u6)
slide-25
SLIDE 25

Quarkonium physics

Preliminaries

such that for it is energetically favorable for the pair The screening length is define as the separation between a

Ls

to be bound (unbound).

q¯ q ` < Ls (` > Ls)

slide-26
SLIDE 26

Preliminaries

such that for it is energetically favorable for the pair The screening length is define as the separation between a

Ls

to be bound (unbound).

q¯ q

We will determine by comparing the action of the pair to

Ls q¯ q

the action of the unbound system; i.e.

Sunb

The screening length is the maximum value of for which is positive

∆S

(In the Euclidean version, this criterion corresponds to determining which configuration has the lowest free energy)

` < Ls (` > Ls) S(`)

∆S(`) = S(`) − Sunb

`

Quarkonium physics

slide-27
SLIDE 27
  • 1. Static case (to warm up)

u = 0 u = uh

θ Given the rotational symmetry in the xy-plane, the most general case is to consider the dipole in the xz-plane.

z x

Choosing the static gauge , , and the string embedding:

τ = t σ = u

Quarkonium physics

Z(u) = z(u) cos θ X(u) = x(u) sin θ

and

slide-28
SLIDE 28
  • 1. Static case

The action for the U-shaped string takes the form S = − L2 2πα0 2 Z dt Z umax du 1 u2 q B(1 + FH cos2 θz02 + F sin2 θx02)

Quarkonium physics

slide-29
SLIDE 29
  • 1. Static case

The action for the U-shaped string takes the form S = − L2 2πα0 2 Z dt Z umax du 1 u2 q B(1 + FH cos2 θz02 + F sin2 θx02) Two conserved momenta and associated to translation

Πz

invariance in the x, z direction. Then, the on-shell action can be written as

Πx

S = − L2 2πα0 2 Z dt Z umax du 1 u2 B √ FH p FBH − u4(Π2

z + HΠ2 x)

where the turning point is determined from the condition

umax

x0(umax) = z0(umax) → ∞

Quarkonium physics

slide-30
SLIDE 30
  • 1. Static case

The action for the U-shaped string takes the form S = − L2 2πα0 2 Z dt Z umax du 1 u2 q B(1 + FH cos2 θz02 + F sin2 θx02) Two conserved momenta and associated to translation

Πz

invariance in the x, z direction. Then, the on-shell action can be written as

Πx

S = − L2 2πα0 2 Z dt Z umax du 1 u2 B √ FH p FBH − u4(Π2

z + HΠ2 x)

which in terms of the momenta reads FBH − u4(Π2

z + HΠ2 x)|umax = 0

umax ≡ umax(a, T, Πi)

Quarkonium physics

slide-31
SLIDE 31
  • 1. Static case

From the boundary conditions we obtain the relation between the

Πz Πx

momenta , and the quark-antiquark separation `

Quarkonium physics

l 2 = Z umax duX0 = Z umax duZ0

slide-32
SLIDE 32
  • 1. Static case

From the boundary conditions we obtain the relation between the

Πz Πx

momenta , and the quark-antiquark separation Finally, to determine , we need to subtract from the U-shaped string action, that of the unbound pair (i.e. two straight strings)

Ls

Sunb = − L2 2πα0 2 Z dt Z uh du √ B u2 The UV divergences associated to integrating all the way to the boundary cancel out in the difference, and there are no IR divergences.

`

Quarkonium physics

l 2 = Z umax duX0 = Z umax duZ0

slide-33
SLIDE 33
  • 1. Static case

Ls ' 0.24 T

[Rey ¡et ¡al; ¡Brandhuber ¡et ¡al] ¡

Ls ' 0.24(π2N 2

c

2s )1/3 The isotropic result ( )

a = 0 H → 0, B → 1, F → f

  • r

at ¡constant ¡entropy ¡ density ¡

Quarkonium physics

slide-34
SLIDE 34
  • 1. Static case

Ls ' 0.24 T

[Rey ¡et ¡al; ¡Brandhuber ¡et ¡al] ¡

Ls ' 0.24(π2N 2

c

2s )1/3 The isotropic result ( )

a = 0 H → 0, B → 1, F → f

  • r

The anisotropic results:

5 10 15 20 25 30 0.5 0.6 0.7 0.8 0.9 1.0

Ls/Liso(T) a/T

5 10 15 20 25 30 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2

Ls/Liso(s) aN 2/3

c

/s1/3

θ = π/2 θ = 0 θ = π/2 θ = 0

Quarkonium physics

slide-35
SLIDE 35
  • 2. Dipole in a plasma wind

We will first consider the ultra-relativistic limit. There are at least two important reasons for doing so: It is relevant for the experiments. It can be understood analytically. For the isotropic case:

Ls(T, v) ∼ (1 − v2)1/4

[Liu ¡et ¡al.] ¡

Quarkonium physics

slide-36
SLIDE 36
  • 2. Dipole in a plasma wind

: direction of the wind (velocity) with respect to the anisotropic direction

θv θ, φ : parametrized the orientation

  • f the dipole

Our set up:

z x y θv θ φ q ¯ q

Dipole rest frame

Quarkonium physics

slide-37
SLIDE 37
  • 2. Dipole in a plasma wind

As in the static case, to determine the screening length we need to compare the actions of a bound and an unbound quark-antiquark pair.

Quarkonium physics

slide-38
SLIDE 38
  • 2. Dipole in a plasma wind

As in the static case, to determine the screening length we need to compare the actions of a bound and an unbound quark-antiquark pair. The unbound action is that of a string moving with constant velocity Let us summarize the important steps:

[MCh, ¡Fernandez, ¡Mateos ¡and ¡Trancanelli] ¡

Quarkonium physics

slide-39
SLIDE 39
  • 2. Dipole in a plasma wind

As in the static case, to determine the screening length we need to compare the actions of a bound and an unbound quark-antiquark pair. The unbound action is that of a string moving with constant velocity The position of the turning point is now Let us summarize the important steps:

[MCh, ¡Fernandez, ¡Mateos ¡and ¡Trancanelli] ¡

It is easy to check that for a fixed separation of the string endpoints, Then the dynamics of the string can be determined using the near boundary expansion of the metric. umax ≡ umax(a, T, Πi, v)

lim

v→1 umax → 0

Quarkonium physics

slide-40
SLIDE 40
  • 2. Dipole in a plasma wind

As in the static case, to determine the screening length we need to compare the actions of a bound and an unbound quark-antiquark pair. The unbound action is that of a string moving with constant velocity The position of the turning point is now Let us summarize the important steps:

[MCh, ¡Fernandez, ¡Mateos ¡and ¡Trancanelli] ¡

It is easy to check that for a fixed separation of the string endpoints, The analytical expression for the near boundary metric

[Mateos ¡et ¡al.] ¡

T - independent F, B, H ∼ 1 + Wi(a)u2 + Gi(a, T)u4 + O(u6) umax ≡ umax(a, T, Πi, v)

lim

v→1 umax → 0

Quarkonium physics

slide-41
SLIDE 41
  • 2. Dipole in a plasma wind

We want to compare the two actions in the ultra-relativistic limit and see how they scale with

(1 − v2)

Quarkonium physics

slide-42
SLIDE 42
  • 2. Dipole in a plasma wind

We want to compare the two actions in the ultra-relativistic limit and see how they scale with

(1 − v2)

After some algebra:

motion outside the transverse plane motion within the transverse plane

∆S(l, v) ∼ (1 − v2)−1/2 × (finite integral) ∆S(l, v) ∼ (1 − v2)−1/4 × (finite integral) I(a, θv, Πi, O(u6)) I(a, T, θv, Πi, O(u6))

⇒ T independent!

Quarkonium physics

slide-43
SLIDE 43
  • 2. Dipole in a plasma wind

And finally, using the boundary conditions, we obtain how the screening length scales in the ultra-relativistic limit: We want to compare the two actions in the ultra-relativistic limit and see how they scale with

(1 − v2)

After some algebra:

motion outside the transverse plane motion within the transverse plane

∆S(l, v) ∼ (1 − v2)−1/2 × (finite integral) ∆S(l, v) ∼ (1 − v2)−1/4 × (finite integral) Ls ⇠    (1 v2)1/2 ⇥ I(a, Πi, O(u6)) if θv 6= π/2 (1 v2)1/4 ⇥ J (a, T, Πi, O(u6)) if θv = π/2

Quarkonium physics

slide-44
SLIDE 44
  • 2. Dipole in a plasma wind

A glimpse of the numerical results: Wind along “z” and dipole along “x”

0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0

Ls/Liso(T)

5 10 15 20 25 30 0.0 0.2 0.4 0.6 0.8 1.0

Ls/Liso(T) a/T

v

a/T = 744 a/T = 12.2 v = 0.9995 v = 0.25

vanishes as (1 − v2)1/4 in the limit v → 1

Ls/Liso

Quarkonium physics

slide-45
SLIDE 45

5 10 15 20 25 30 0.6 0.7 0.8 0.9 1.0 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0

  • 2. Dipole in a plasma wind

A glimpse of the numerical results: Wind along “x” and dipole along “z”

Ls/Liso(T) Ls/Liso(T) a/T

v

a/T = 744 a/T = 12.2 v = 0.9995 v = 0.25

approaches a finite, non-zero value as v → 1

Ls/Liso

Quarkonium physics

slide-46
SLIDE 46

Another interesting limit to consider: a/T 1

Quarkonium physics

slide-47
SLIDE 47

The proper velocity along “z” of a point on the string at some

u

vproper(u) = vz s H(u) F(u)B(u) increases from to , more steeply as increases, a/T

H(u) u = 0 u = uh F(u)B(u) has the opposite behavior. ¡

Maximum value of beyond which becomes superluminal

umax vproper

Quarkonium physics

Another interesting limit to consider: a/T 1 Two observations:

slide-48
SLIDE 48

The proper velocity along “z” of a point on the string at some

u

vproper(u) = vz s H(u) F(u)B(u) increases from to , more steeply as increases, a/T

H(u) u = 0 u = uh F(u)B(u) has the opposite behavior. ¡

Maximum value of beyond which becomes superluminal

umax vproper

decreases as increases.

umax

a/T We can show that for ,

vz 6= 0

Quarkonium physics

Another interesting limit to consider: a/T 1 Two observations:

lim

a/T 1 umax → 0

More over, Use near boundary metric to study Ls

slide-49
SLIDE 49

For a fixed value of vz 6= 0

Ls ∼ a−1

It is straight forward to check that T independent as in the ultra-relativistic

case

Quarkonium physics

Another interesting limit to consider: a/T 1

slide-50
SLIDE 50

For a fixed value of vz 6= 0

Ls ∼ a−1

T = 0 a meson will dissociate for sufficiently large anisotropy ¡

adiss

It is straight forward to check that The limit can be understood as at fixed , or

a/T 1 a → ∞ T

as at fixed .

T → 0 a

Quarkonium physics

Another interesting limit to consider: a/T 1

slide-51
SLIDE 51

For a fixed value of vz 6= 0

Ls ∼ a−1

The limit can be understood as at fixed , or

a/T 1 a → ∞ T

as at fixed .

T → 0 a T = 0 a meson will dissociate for sufficiently large anisotropy ¡

adiss

For motion within the transverse plane ( )

vz = 0 Ls ∼ f(a, T)

No when

adiss T → 0

It is straight forward to check that

Quarkonium physics

Another interesting limit to consider: a/T 1

slide-52
SLIDE 52

So far we have studied , but clearly we could also think of

Ls(a, T, v)

and i.e. characterizes the dissociation of a qq-pair of fixed size in a plasma with a given degree of anisotropy . Analogously ¡for ¡ ¡ ¡ ¡ ¡ ¡ ¡ ¡ ¡ ¡ ¡ ¡ ¡.

a adiss

Using our results for the screening length, we can study the behavior and adiss(T, `, v)

Tdiss(a, `, v) adiss(T, `, v) Tdiss(a, `, v) Tdiss(a, `, v) `

Quarkonium physics

  • f
slide-53
SLIDE 53

2 4 6 8 10 12 0.00 0.05 0.10 0.15 0.20 0.25

0.00 0.05 0.10 0.15 0.20 0.25 2 4 6 8 10 12

lTdiss l a

diss

al Tl x x x x z z z z Numerical results: As explained before, even at zero temperature a meson of size will dissociate if the anisotropy is increased above and the proportionality constant is a decreasing function of the velocity

at ¡rest ¡ moving ¡along ¡the ¡z-­‑direc>on ¡( ¡ ¡ ¡ ¡ ¡ ¡ ¡ ¡ ¡ ¡ ¡ ¡ ¡ ¡ ¡ ¡ ¡ ¡) ¡ v = 0.45

adiss(T = 0, `) ∝ 1/`

Quarkonium physics

`

slide-54
SLIDE 54

Quarkonium physics

Numerical results for motion within the transverse plane:

0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0

Tdiss(v)/Tdiss(0) Tdiss(v)/Tdiss(0) v v extending ¡along ¡z ¡ extending ¡along ¡x ¡

The behavior is qualitatively analogous to that of the isotropic case.

Tdiss(v) ' Tdiss(0)(1 v2)1/4

[Liu ¡et ¡al.] ¡

slide-55
SLIDE 55

Quarkonium physics

0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0

Tdiss(v)/Tdiss(0) Tdiss(v)/Tdiss(0) v v extending ¡along ¡z ¡ extending ¡along ¡x ¡ a = 0 a = 0 a = 25

There is a limiting velocity even at zero temperature!

vlim < 1

Numerical results for motion outside the transverse plane:

slide-56
SLIDE 56

5 10 15 20 0.0 0.2 0.4 0.6 0.8 1.0

Limiting velocity for a fixed anisotropy and , meson oriented along the x-direction and moving along the z-direction vlim al

T = 0

Quarkonium physics

slide-57
SLIDE 57

5 10 15 20 0.0 0.2 0.4 0.6 0.8 1.0

Limiting velocity for a fixed anisotropy and , meson oriented along the x-direction and moving along the z-direction vlim al

T = 0

In the case of ultra-relativistic motion and :

a/T 1

adiss ⇠ 1 ` (1 v2

lim)1/2 if ✓v 6= ⇡/2

Tdiss ∼ 1 ` (1 − v2

lim)1/4 if ✓v = ⇡/2

Quarkonium physics

slide-58
SLIDE 58

Conclusions

  • We have completely characterized the screening length for quarkonium

mesons moving with arbitrary velocities and orientations.

  • Mesons dissociate above certain critical value of the anisotropy, even

at zero temperature.

  • There is a limiting velocity for mesons moving through the plasma, even

at zero temperature .

  • The gravity calculation involves only the coupling of the string to the

background metric, so any anisotropy that gives rise to a qualitatively similar metric will yield qualitatively similar results.

slide-59
SLIDE 59

Back up slide

F = 1 + 11 24a2u2 + ✓ F4 + 7 12a4 log u ◆ u4 + O(u6) , B = 1 − 11 24a2u2 + ✓ B4 − 7 12a4 log u ◆ u4 + O(u6) , H = 1 + 1 4a2u2 − ✓2 7B4 − 5 4032a4 − 1 6a4 log u ◆ u4 + O(u6) The near-boundary behavior of metric functions: B4(a, T) F4(a, T) and with