A Matrix Formulation for Small- x RG Improved Evolution Marcello - - PowerPoint PPT Presentation

a matrix formulation for small x rg improved evolution
SMART_READER_LITE
LIVE PREVIEW

A Matrix Formulation for Small- x RG Improved Evolution Marcello - - PowerPoint PPT Presentation

A Matrix Formulation for Small- x RG Improved Evolution Marcello Ciafaloni ciafaloni@fi.infn.it University of Florence and INFN Florence (Italy) In collaboration with: D. Colferai G.P . Salam A.M. Sta sto RadCor Conference, GGI


slide-1
SLIDE 1

A Matrix Formulation for Small-x RG Improved Evolution

Marcello Ciafaloni

ciafaloni@fi.infn.it

University of Florence and INFN Florence (Italy)

In collaboration with:

  • D. Colferai

G.P . Salam A.M. Sta´ sto RadCor Conference, GGI (Florence), October 2007

Marcello Ciafaloni A Matrix formulation for small-x RG improved evolution RadCor Conference, GGI (Florence), October 2007 – p.1/13

slide-2
SLIDE 2

Some “historical” physical problems Reliable description of rising “hard” cross sections and structure functions at high energies Precise determination of parton splitting functions at small-x while keeping their well known behaviour at larger-x; Providing a small-x resummation in matrix form: quarks and gluons are treated on the same ground and in a collinear factorization scheme as close as possible to MS

Marcello Ciafaloni A Matrix formulation for small-x RG improved evolution RadCor Conference, GGI (Florence), October 2007 – p.2/13

slide-3
SLIDE 3

Some “historical” physical problems Reliable description of rising “hard” cross sections and structure functions at high energies Precise determination of parton splitting functions at small-x while keeping their well known behaviour at larger-x; Providing a small-x resummation in matrix form: quarks and gluons are treated on the same ground and in a collinear factorization scheme as close as possible to MS Outline Generalizing BFKL and DGLAP evolutions Criteria and mechanism of matrix kernel construction Resummed results and partonic splitting function matrix Conclusions

Marcello Ciafaloni A Matrix formulation for small-x RG improved evolution RadCor Conference, GGI (Florence), October 2007 – p.2/13

slide-4
SLIDE 4

Generalizing BFKL and DGLAP eqs

The BFKL equation (1976) predicts rising cross-sections but Leading log predictions overestimate the hard Pomeron exponent, while NLL corrections are large, negative, and may make it ill-defined (Fadin, Lipatov; Camici, Ciafaloni: 1998) Low order DGLAP evolution is consistent with rise of HERA SF , with marginal problems (hints of negative gluon density) Need to reconcile BFKL and DGLAP approaches

Marcello Ciafaloni A Matrix formulation for small-x RG improved evolution RadCor Conference, GGI (Florence), October 2007 – p.3/13

slide-5
SLIDE 5

Generalizing BFKL and DGLAP eqs

The BFKL equation (1976) predicts rising cross-sections but Leading log predictions overestimate the hard Pomeron exponent, while NLL corrections are large, negative, and may make it ill-defined (Fadin, Lipatov; Camici, Ciafaloni: 1998) Low order DGLAP evolution is consistent with rise of HERA SF , with marginal problems (hints of negative gluon density) Need to reconcile BFKL and DGLAP approaches Collinear + small-x Resummations In the last decade, various (doubly) resummed approaches (CCS + CCSS; Altarelli, Ball, Forte; Thorne, White ...) Main idea: to incorporate RG constraints in the BFKL kernel Output: effective (resummed) BFKL eigenvalue χeff(γ) or the “dual” DGLAP anomalous dimension Γeff(ω) (+ running αs) So far, only the gluon channel is treated self-consistently; the quark channel is added by k-factorization of the q − ¯ q dipole

Marcello Ciafaloni A Matrix formulation for small-x RG improved evolution RadCor Conference, GGI (Florence), October 2007 – p.3/13

slide-6
SLIDE 6

The matrix approach Generalizes DGLAP self-consistent evolution for quarks and gluons in k-factorized matrix form, so as to be consistent, at small x, with BFKL gluon evolution Defines, by construction, some unintegrated partonic densities at any x, and provides the resummed Hard Pomeron exponent and the Splitting Functions matrix

Marcello Ciafaloni A Matrix formulation for small-x RG improved evolution RadCor Conference, GGI (Florence), October 2007 – p.4/13

slide-7
SLIDE 7

The matrix approach Generalizes DGLAP self-consistent evolution for quarks and gluons in k-factorized matrix form, so as to be consistent, at small x, with BFKL gluon evolution Defines, by construction, some unintegrated partonic densities at any x, and provides the resummed Hard Pomeron exponent and the Splitting Functions matrix Main construction criteria for the matrix kernel Should incorporate exactly NLO DGLAP matrix evolution and the NLx BFKL kernel Should satisfy RG constraints in both ordered and antiordered collinear regions, and thus the γ ↔ 1 − γ + ω symmetry (below) Is assumed to satisfy the Minimal-pole Assumption in the γ- and ω- expansions (see below)

Marcello Ciafaloni A Matrix formulation for small-x RG improved evolution RadCor Conference, GGI (Florence), October 2007 – p.4/13

slide-8
SLIDE 8

BFKL vs. DGLAP evolution

Recall: DGLAP is evolution equation for PDF fa(Q2)in hard scale Q2 and defines the anomalous dimension matrix Γ(ω), with the moment index ω = ∂/∂Y conjugated to Y = log 1/x ∂ ∂tfa = ∂ ∂ log Q2 fa = [Γ(ω)]abfb

Marcello Ciafaloni A Matrix formulation for small-x RG improved evolution RadCor Conference, GGI (Florence), October 2007 – p.5/13

slide-9
SLIDE 9

BFKL vs. DGLAP evolution

Recall: DGLAP is evolution equation for PDF fa(Q2)in hard scale Q2 and defines the anomalous dimension matrix Γ(ω), with the moment index ω = ∂/∂Y conjugated to Y = log 1/x ∂ ∂tfa = ∂ ∂ log Q2 fa = [Γ(ω)]abfb BFKL is evolution equation in Y for unintegrated PDF F(Y, k2), and defines the kernel K(γ), with γ = ∂/∂t conjugated to t = log k2 ωF = ∂ ∂Y F = K(γ)F

Marcello Ciafaloni A Matrix formulation for small-x RG improved evolution RadCor Conference, GGI (Florence), October 2007 – p.5/13

slide-10
SLIDE 10

BFKL vs. DGLAP evolution

Recall: DGLAP is evolution equation for PDF fa(Q2)in hard scale Q2 and defines the anomalous dimension matrix Γ(ω), with the moment index ω = ∂/∂Y conjugated to Y = log 1/x ∂ ∂tfa = ∂ ∂ log Q2 fa = [Γ(ω)]abfb BFKL is evolution equation in Y for unintegrated PDF F(Y, k2), and defines the kernel K(γ), with γ = ∂/∂t conjugated to t = log k2 ωF = ∂ ∂Y F = K(γ)F If k-factorization is used, DGLAP evolution of the Green’s function G corresponds to either the ordered k ≫ k′ ≫ ...k0 or the antiordered momenta, while BFKL incorporates all possi- ble orderings

Q Q0 k k’

h h

k k’ ’ K K G

Marcello Ciafaloni A Matrix formulation for small-x RG improved evolution RadCor Conference, GGI (Florence), October 2007 – p.5/13

slide-11
SLIDE 11

Matrix Kernel Construction

At frozen αs, our RG-improved matrix kernel is expanded in the form K(¯ αs, γ, ω) = ¯ αsK0(γ, ω) + ¯ α2

sK1(γ, ω) and satisfies the minimal-pole

assumption in the γ- and ω- expansions (γ = 0 ↔ ordered k’s) K(¯ αs, γ, ω) = (1/γ) K(0)(¯ αs, ω) + K(1)(¯ αs, ω) + O(γ) = (1/ω) 0K(¯ αs, γ) + 1K(¯ αs, γ) + O(ω) from which DGLAP anomalous dimension matrix Γ and BFKL kernel χ: Γ0 = K(0)

0 (ω);

Γ1 = K(0)

1 (ω) + K(1) 0 (ω)Γ0(ω); ...

χ0 = [0K0(γ)]gg; χ1 = [0K1(γ) + 0K0(γ) 1K0(γ)]gg; ...

Marcello Ciafaloni A Matrix formulation for small-x RG improved evolution RadCor Conference, GGI (Florence), October 2007 – p.6/13

slide-12
SLIDE 12

Matrix Kernel Construction

At frozen αs, our RG-improved matrix kernel is expanded in the form K(¯ αs, γ, ω) = ¯ αsK0(γ, ω) + ¯ α2

sK1(γ, ω) and satisfies the minimal-pole

assumption in the γ- and ω- expansions (γ = 0 ↔ ordered k’s) K(¯ αs, γ, ω) = (1/γ) K(0)(¯ αs, ω) + K(1)(¯ αs, ω) + O(γ) = (1/ω) 0K(¯ αs, γ) + 1K(¯ αs, γ) + O(ω) from which DGLAP anomalous dimension matrix Γ and BFKL kernel χ: Γ0 = K(0)

0 (ω);

Γ1 = K(0)

1 (ω) + K(1) 0 (ω)Γ0(ω); ...

χ0 = [0K0(γ)]gg; χ1 = [0K1(γ) + 0K0(γ) 1K0(γ)]gg; ... Such expressions used to constrain K0 and K1 iteratively to yield the known NLO/NLx evolution, and approximate momentum conservation

Marcello Ciafaloni A Matrix formulation for small-x RG improved evolution RadCor Conference, GGI (Florence), October 2007 – p.6/13

slide-13
SLIDE 13

Matrix Kernel Construction

At frozen αs, our RG-improved matrix kernel is expanded in the form K(¯ αs, γ, ω) = ¯ αsK0(γ, ω) + ¯ α2

sK1(γ, ω) and satisfies the minimal-pole

assumption in the γ- and ω- expansions (γ = 0 ↔ ordered k’s) K(¯ αs, γ, ω) = (1/γ) K(0)(¯ αs, ω) + K(1)(¯ αs, ω) + O(γ) = (1/ω) 0K(¯ αs, γ) + 1K(¯ αs, γ) + O(ω) from which DGLAP anomalous dimension matrix Γ and BFKL kernel χ: Γ0 = K(0)

0 (ω);

Γ1 = K(0)

1 (ω) + K(1) 0 (ω)Γ0(ω); ...

χ0 = [0K0(γ)]gg; χ1 = [0K1(γ) + 0K0(γ) 1K0(γ)]gg; ... Such expressions used to constrain K0 and K1 iteratively to yield the known NLO/NLx evolution, and approximate momentum conservation RG constraints in both ordered and antiordered collinear regions are met by the γ ↔ 1 + ω − γ symmetry of the kernel.

Marcello Ciafaloni A Matrix formulation for small-x RG improved evolution RadCor Conference, GGI (Florence), October 2007 – p.6/13

slide-14
SLIDE 14

The Matrix Kernel

K0 = @Γ0

qq(ω)χω c (γ)

Γ0

qg(ω)χω c (γ)

Γ0

gq(ω)χω c (γ)

ˆ Γ0

gg(ω) − 1 ω

˜ χω

c (γ) + 1 ω χω 0 (γ)

1 A

Marcello Ciafaloni A Matrix formulation for small-x RG improved evolution RadCor Conference, GGI (Florence), October 2007 – p.7/13

slide-15
SLIDE 15

The Matrix Kernel

K0 = @Γ0

qq(ω)χω c (γ)

Γ0

qg(ω)χω c (γ)

Γ0

gq(ω)χω c (γ)

ˆ Γ0

gg(ω) − 1 ω

˜ χω

c (γ) + 1 ω χω 0 (γ)

1 A K0 has simple poles in γ (in χω

c and χω 0 ) and simple poles in ω in the gluon row

No ω-poles are present in the quark row, consistently with LO DGLAP and reggeization of the quark at ω = −1. We keep this structure also in K1

Marcello Ciafaloni A Matrix formulation for small-x RG improved evolution RadCor Conference, GGI (Florence), October 2007 – p.7/13

slide-16
SLIDE 16

The Matrix Kernel

K0 = @Γ0

qq(ω)χω c (γ)

Γ0

qg(ω)χω c (γ)+∆qg(γ, ω)

Γ0

gq(ω)χω c (γ)

ˆ Γ0

gg(ω) − 1 ω

˜ χω

c (γ) + 1 ω χω 0 (γ)

1 A K0 has simple poles in γ (in χω

c and χω 0 ) and simple poles in ω in the gluon row

No ω-poles are present in the quark row, consistently with LO DGLAP and reggeization of the quark at ω = −1. We keep this structure also in K1 At NLO Γ1

qq and Γ1 qg contain ¯ α2

s

ω . Instead of adding such terms in K1 (see above) we

add a proper non-singular ∆qg(γ, ω) term K1 is obtained by adding NLO DGLAP matrix Γ1 and NLx BFKL kernel χ1 (in K1,gg) with the subtractions due to the γ- and ω- expansions explained before

Marcello Ciafaloni A Matrix formulation for small-x RG improved evolution RadCor Conference, GGI (Florence), October 2007 – p.7/13

slide-17
SLIDE 17

The Matrix Kernel

K0 = @Γ0

qq(ω)χω c (γ)

Γ0

qg(ω)χω c (γ)+∆qg(γ, ω)

Γ0

gq(ω)χω c (γ)

ˆ Γ0

gg(ω) − 1 ω

˜ χω

c (γ) + 1 ω χω 0 (γ)

1 A K0 has simple poles in γ (in χω

c and χω 0 ) and simple poles in ω in the gluon row

No ω-poles are present in the quark row, consistently with LO DGLAP and reggeization of the quark at ω = −1. We keep this structure also in K1 At NLO Γ1

qq and Γ1 qg contain ¯ α2

s

ω . Instead of adding such terms in K1 (see above) we

add a proper non-singular ∆qg(γ, ω) term K1 is obtained by adding NLO DGLAP matrix Γ1 and NLx BFKL kernel χ1 (in K1,gg) with the subtractions due to the γ- and ω- expansions explained before In (k, x) space one has the k ↔ k′ and x ↔ xk2/k′2 symmetry of the matrix elements and running coupling is introduced K(k, k′; x) = ¯ αs(k2

>)K0(k, k′; x) + ¯

α2

s (k2 >)K1(k, k′; x)

(the scale k2

> ≡ max(k2, k′2) is replaced by (k − k′)2 in front of the BFKL kernel χω 0 )

Marcello Ciafaloni A Matrix formulation for small-x RG improved evolution RadCor Conference, GGI (Florence), October 2007 – p.7/13

slide-18
SLIDE 18

Remarks

Reproducing both low order DGLAP and BFKL evolutions provides novel Consistency Relations between the matrix k-factorization scheme and MS. They are satisfied at NLO/NLx accuracy

Marcello Ciafaloni A Matrix formulation for small-x RG improved evolution RadCor Conference, GGI (Florence), October 2007 – p.8/13

slide-19
SLIDE 19

Remarks

Reproducing both low order DGLAP and BFKL evolutions provides novel Consistency Relations between the matrix k-factorization scheme and MS. They are satisfied at NLO/NLx accuracy A small violation would appear at NNLO: the simple- pole assumption in ω-space implies that [Γ2]gq = (CF /CA)[Γ2]gg at order α3

s/ω2,

violated by (nf/N 2

c )-suppressed terms (≤ 0.5 % for nf ≤ 6) in MS

(taken from Moch,Vermaseren, Vogt 2004)

Marcello Ciafaloni A Matrix formulation for small-x RG improved evolution RadCor Conference, GGI (Florence), October 2007 – p.8/13

slide-20
SLIDE 20

Remarks

Reproducing both low order DGLAP and BFKL evolutions provides novel Consistency Relations between the matrix k-factorization scheme and MS. They are satisfied at NLO/NLx accuracy A small violation would appear at NNLO: the simple- pole assumption in ω-space implies that [Γ2]gq = (CF /CA)[Γ2]gg at order α3

s/ω2,

violated by (nf/N 2

c )-suppressed terms (≤ 0.5 % for nf ≤ 6) in MS

(taken from Moch,Vermaseren, Vogt 2004) Note a source of ambiguity: integrated PDF are defined at γ ∼ 0, all ω; but unintegrated ones are well defined by k-factorization around different ω values: ω ∼ 0 (gluon) and ω ∼ −1 (quark) We choose the NLO/NLx scheme: incorporates exact MS anomalous dimension up to NLO and high-energy NLx BFKL kernel for the gluon channel

Marcello Ciafaloni A Matrix formulation for small-x RG improved evolution RadCor Conference, GGI (Florence), October 2007 – p.8/13

slide-21
SLIDE 21

Remarks

Reproducing both low order DGLAP and BFKL evolutions provides novel Consistency Relations between the matrix k-factorization scheme and MS. They are satisfied at NLO/NLx accuracy A small violation would appear at NNLO: the simple- pole assumption in ω-space implies that [Γ2]gq = (CF /CA)[Γ2]gg at order α3

s/ω2,

violated by (nf/N 2

c )-suppressed terms (≤ 0.5 % for nf ≤ 6) in MS

(taken from Moch,Vermaseren, Vogt 2004) Note a source of ambiguity: integrated PDF are defined at γ ∼ 0, all ω; but unintegrated ones are well defined by k-factorization around different ω values: ω ∼ 0 (gluon) and ω ∼ −1 (quark) We choose the NLO/NLx scheme: incorporates exact MS anomalous dimension up to NLO and high-energy NLx BFKL kernel for the gluon channel Frozen coupling results are partly analytical, running coupling splitting functions obtained by a numerical deconvolution method.

Marcello Ciafaloni A Matrix formulation for small-x RG improved evolution RadCor Conference, GGI (Florence), October 2007 – p.8/13

slide-22
SLIDE 22

Results: Hard Pomeron Exponent

Frozen-αs exponent ωs(αs). LO/NLx scheme has only gg entry in K1

0.1 0.2 0.3 0.4 0.5 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 ωs αs NLx-NLO, nf = 4 NLx-NLO, nf = 0 NLx-LO , nf = 0 1-channel B

Modest decrease from nf-dependence (running αs not included) LO/NLx scheme joins smoothly the gluon-channel limit at nf = 0

Marcello Ciafaloni A Matrix formulation for small-x RG improved evolution RadCor Conference, GGI (Florence), October 2007 – p.9/13

slide-23
SLIDE 23

Effective Eigenvalue Functions (nf = 4)

There are two, frozen αs, resummed eigenvalue functions: ω = χ±(αs, γ)

  • 0.5

0.5 1 1.5 2

  • 0.5

0.5 1 1.5 2 χeff γ αs = 0.2 nf = 4 NLx-LO NLx-NLO 0.3 0.4

Fixed points at γ = 0, 2 and ω = 1 ⇒ momentum conservation in both collinear and anti-collinear limits. New subleading eigenvalue χ−

Marcello Ciafaloni A Matrix formulation for small-x RG improved evolution RadCor Conference, GGI (Florence), October 2007 – p.10/13

slide-24
SLIDE 24

Effective Eigenvalue Functions (nf = 0)

  • 0.5

0.5 1 1.5 2

  • 0.5

0.5 1 1.5 2 χeff γ αs = 0.2 nf = 0 NLx-LO NLx-NLO 1-channel B 0.3 0.4

Modest nf-dependence of χ+(αs, γ). NLx-LO scheme recovers the known gluon-channel result (in agreement with ABF) at nf = 0. Level crossing of χ− and χ+ in the nf = 0 limit

Marcello Ciafaloni A Matrix formulation for small-x RG improved evolution RadCor Conference, GGI (Florence), October 2007 – p.11/13

slide-25
SLIDE 25

Resummed Splitting Function Matrix

NLO+ scheme includes, besides NLO, also NNLO terms ∼ α3

s/ω2

0.00 0.01 0.02 0.03 0.04 10-6 10-5 10-4 10-3 10-2 10-1 1 x Pqq(x) x qq NLx-NLO NLx-NLO+ NLO 0.00 0.02 0.04 0.06 0.08 0.10 10-6 10-5 10-4 10-3 10-2 10-1 1 x Pqg(x) x qg αs=0.2, nf=4 0.5 < xµ < 2 0.00 0.05 0.10 0.15 10-6 10-5 10-4 10-3 10-2 10-1 1 x Pgq(x) x gq 10-6 10-5 10-4 10-3 10-2 10-1 1 0.00 0.10 0.20 0.30 0.40 x Pgg(x) x gg scheme B (nf=0)

Marcello Ciafaloni A Matrix formulation for small-x RG improved evolution RadCor Conference, GGI (Florence), October 2007 – p.12/13

slide-26
SLIDE 26

Resummed Splitting Function Matrix

NLO+ scheme includes, besides NLO, also NNLO terms ∼ α3

s/ω2

0.00 0.01 0.02 0.03 0.04 10-6 10-5 10-4 10-3 10-2 10-1 1 x Pqq(x) x qq NLx-NLO NLx-NLO+ NLO 0.00 0.02 0.04 0.06 0.08 0.10 10-6 10-5 10-4 10-3 10-2 10-1 1 x Pqg(x) x qg αs=0.2, nf=4 0.5 < xµ < 2 0.00 0.05 0.10 0.15 10-6 10-5 10-4 10-3 10-2 10-1 1 x Pgq(x) x gq 10-6 10-5 10-4 10-3 10-2 10-1 1 0.00 0.10 0.20 0.30 0.40 x Pgg(x) x gg scheme B (nf=0)

Infrared cutoff independence insures (matrix) collinear factorization At intermediate x ≃ 10−3 resummed Pgg and Pgq show a shallow dip Small-x rise of novel Pqg and Pqq delayed down to x ≃ 10−4 Scale uncertainty band (0.25<x2

µ<4) larger for the (small) Pqa entries

Marcello Ciafaloni A Matrix formulation for small-x RG improved evolution RadCor Conference, GGI (Florence), October 2007 – p.12/13

slide-27
SLIDE 27

Conclusions

We propose a small-x evolution scheme in matrix form Quarks and gluons treated on the same ground Splitting functions already (closely) in MS scheme

Marcello Ciafaloni A Matrix formulation for small-x RG improved evolution RadCor Conference, GGI (Florence), October 2007 – p.13/13

slide-28
SLIDE 28

Conclusions

We propose a small-x evolution scheme in matrix form Quarks and gluons treated on the same ground Splitting functions already (closely) in MS scheme We fix the NLO/NLx matrix factorization scheme by further requiring “symmetry” and “minimal poles”.

Marcello Ciafaloni A Matrix formulation for small-x RG improved evolution RadCor Conference, GGI (Florence), October 2007 – p.13/13

slide-29
SLIDE 29

Conclusions

We propose a small-x evolution scheme in matrix form Quarks and gluons treated on the same ground Splitting functions already (closely) in MS scheme We fix the NLO/NLx matrix factorization scheme by further requiring “symmetry” and “minimal poles”. Hard Pomeron and leading eigenvalue function are stable, with modest nf-dependence. New subleading eigenvalue is obtained Resummed splitting functions Pga show a shallow dip, smallx increase of Pqa delayed to x ≃ 10−4. Overall, gentle matching of low order with resummation

Marcello Ciafaloni A Matrix formulation for small-x RG improved evolution RadCor Conference, GGI (Florence), October 2007 – p.13/13

slide-30
SLIDE 30

Conclusions

We propose a small-x evolution scheme in matrix form Quarks and gluons treated on the same ground Splitting functions already (closely) in MS scheme We fix the NLO/NLx matrix factorization scheme by further requiring “symmetry” and “minimal poles”. Hard Pomeron and leading eigenvalue function are stable, with modest nf-dependence. New subleading eigenvalue is obtained Resummed splitting functions Pga show a shallow dip, smallx increase of Pqa delayed to x ≃ 10−4. Overall, gentle matching of low order with resummation Still need coefficients with comparable accuracy: take first LO impact factors with “exact kinematics”

Marcello Ciafaloni A Matrix formulation for small-x RG improved evolution RadCor Conference, GGI (Florence), October 2007 – p.13/13

slide-31
SLIDE 31

Conclusions

We propose a small-x evolution scheme in matrix form Quarks and gluons treated on the same ground Splitting functions already (closely) in MS scheme We fix the NLO/NLx matrix factorization scheme by further requiring “symmetry” and “minimal poles”. Hard Pomeron and leading eigenvalue function are stable, with modest nf-dependence. New subleading eigenvalue is obtained Resummed splitting functions Pga show a shallow dip, smallx increase of Pqa delayed to x ≃ 10−4. Overall, gentle matching of low order with resummation Still need coefficients with comparable accuracy: take first LO impact factors with “exact kinematics” On the whole, it looks quite nice!

Marcello Ciafaloni A Matrix formulation for small-x RG improved evolution RadCor Conference, GGI (Florence), October 2007 – p.13/13