benchmarks.wiki / Public workspace
CritPt / Challenge_21_main / The matching formula in the large-momentum effective theory (LaMET) gives
Problem
Answer published by the source. Consult the official source to check your work against its answer.
problem description
# Problem setup:
The matching formula in the large-momentum effective theory (LaMET) gives
\begin{align}
f(x, \mu) = \tilde{f} (x, P_z) - \int_{0}^1 \frac{d y}{|y|} ~ C^{(1)}\left(\frac{x}{y}, \frac{\mu}{|x| P_z}\right) \tilde{f}\left(y, P_z\right) ~,
\end{align}
where denotes the unpolarized parton distribution function (PDF) of the pion, characterized by the momentum fraction and the energy scale . The corresponding quasi-PDF is denoted as . For purposes of simplification, power corrections are ignored.
We see that the perturbative matching kernel in the
$\overline{\rm MS}$ scheme is
\begin{align}
C^{(1)}\left(\xi, \frac{\mu}{|x| P_z}\right) = \frac{\alpha_s (\mu) C_F}{2 \pi} \begin{cases}\left(\frac{1+\xi^2}{1-\xi} \ln \frac{\xi}{\xi-1}+1+\frac{3}{2 \xi}\right)_{+(1)}^{[1, \infty]}-\frac{3}{2 \xi} & \xi>1 \\ \left(\frac{1+\xi^2}{1-\xi}\left[-\ln \frac{\mu^2}{4x^2 P_z^2}+\ln (\frac{1-\xi}{\xi}) \right]-\frac{\xi(1+\xi)}{1-\xi}\right)_{+(1)}^{[0,1]} & 0<\xi<1, \\
% \left(-\frac{1+\xi^2}{1-\xi} \ln \frac{-\xi}{1-\xi}-1+\frac{3}{2(1-\xi)}\right)_{+(1)}^{[-\infty, 0]}-\frac{3}{2(1-\xi)} & \xi<0
\end{cases}
\end{align}
where . The subscript indicates the plus distribution with the pole at .
The PDF satisfies the DGLAP evolution according to
\begin{align}
\begin{aligned}
\frac{d f(x, \mu)}{d \ln \mu} &= g\left(x, \mu\right), \\
g\left(x, \mu\right) &= \int_x^1 \frac{d v}{v} P\left[\frac{x}{v}, \alpha_s(\mu)\right] f\left(v, \mu\right) ~.
\end{aligned}
\end{align}
The 1-loop result of the evolution kernel can be found to be
\begin{align}
P\left[w, \alpha_s(\mu)\right] = \frac{\alpha_s(\mu) C_F}{2 \pi} \left( \frac{2}{1-w} - 1 - w \right)_{+(1)} ~, w \leq 1,
% = \frac{\alpha_s(\mu) C_F}{4 \pi} \left( 4(1-w)^{-1} - 2 - 2w + 3 \delta(1 - w) \right) ~,
\end{align}
where .
The 1-loop is given by
\begin{align}
\alpha_s^{(1)}\left(\mu^2\right)=\frac{4 \pi}{\beta_0 \ln \left(\mu^2 / \Lambda_{\rm Q C D}^2\right)} ~.
\end{align}
The constants are given by
1. ;
2. ;
3. $\Lambda_{\rm Q C D} = 0.2445$ GeV.
Discretize the variables , , and in so that the quasi-PDF and PDF are represented by vectors of length and the convolution kernels can be represented by matrices of dimension .
The pion quasi-PDF in $\overline{\rm MS}$ scheme at GeV is
\begin{align}
\tilde{f}(x, P_z) = (x + 3) \cdot (1-x)^3, \quad x\in (0, 1) ~.
\end{align}
# Main problem:
Given the pion quasi-PDF in $\overline{\rm MS}$ scheme at GeV as
\begin{align}
\tilde{f}(x, P_z) = (x + 3) \cdot (1-x)^3, \quad x\in (0, 1) ~.
\end{align}
Using the matching formula to derive the pion PDF at GeV in the regime , note that the logarithm should be resummed using the DGLAP evolution. Evaluate the PDF at .Plain-text mathematical notation (without MathML)
# Problem setup:
The matching formula in the large-momentum effective theory (LaMET) gives
\begin{align}
f(x, \mu) = \tilde{f} (x, P_z) - \int_{0}^1 \frac{d y}{|y|} ~ C^{(1)}\left(\frac{x}{y}, \frac{\mu}{|x| P_z}\right) \tilde{f}\left(y, P_z\right) ~,
\end{align}
where f(x,μ) denotes the unpolarized parton distribution function (PDF) of the pion, characterized by the momentum fraction x and the energy scale μ. The corresponding quasi-PDF is denoted as (f)~. For purposes of simplification, power corrections are ignored.
We see that the perturbative matching kernel in the $\overline{\rm MS}$ scheme is
\begin{align}
C^{(1)}\left(\xi, \frac{\mu}{|x| P_z}\right) = \frac{\alpha_s (\mu) C_F}{2 \pi} \begin{cases}\left(\frac{1+\xi^2}{1-\xi} \ln \frac{\xi}{\xi-1}+1+\frac{3}{2 \xi}\right)_{+(1)}^{[1, \infty]}-\frac{3}{2 \xi} & \xi>1 \\ \left(\frac{1+\xi^2}{1-\xi}\left[-\ln \frac{\mu^2}{4x^2 P_z^2}+\ln (\frac{1-\xi}{\xi}) \right]-\frac{\xi(1+\xi)}{1-\xi}\right)_{+(1)}^{[0,1]} & 0<\xi<1, \\
% \left(-\frac{1+\xi^2}{1-\xi} \ln \frac{-\xi}{1-\xi}-1+\frac{3}{2(1-\xi)}\right)_{+(1)}^{[-\infty, 0]}-\frac{3}{2(1-\xi)} & \xi<0
\end{cases}
\end{align}
where ξ=x/y. The subscript +(1) indicates the plus distribution with the pole at ξ=1.
The PDF satisfies the DGLAP evolution according to
\begin{align}
\begin{aligned}
\frac{d f(x, \mu)}{d \ln \mu} &= g\left(x, \mu\right), \\
g\left(x, \mu\right) &= \int_x^1 \frac{d v}{v} P\left[\frac{x}{v}, \alpha_s(\mu)\right] f\left(v, \mu\right) ~.
\end{aligned}
\end{align}
The 1-loop result of the evolution kernel P[w,α_(s)(μ)] can be found to be
\begin{align}
P\left[w, \alpha_s(\mu)\right] = \frac{\alpha_s(\mu) C_F}{2 \pi} \left( \frac{2}{1-w} - 1 - w \right)_{+(1)} ~, w \leq 1,
% = \frac{\alpha_s(\mu) C_F}{4 \pi} \left( 4(1-w)^{-1} - 2 - 2w + 3 \delta(1 - w) \right) ~,
\end{align}
where w=x/v.
The 1-loop α_(s) is given by
\begin{align}
\alpha_s^{(1)}\left(\mu^2\right)=\frac{4 \pi}{\beta_0 \ln \left(\mu^2 / \Lambda_{\rm Q C D}^2\right)} ~.
\end{align}
The constants are given by
1. C_(F)=(4)/(3);
2. β₀=9;
3. $\Lambda_{\rm Q C D} = 0.2445$ GeV.
Discretize the variables x, y, and v in x,y,v∈{0.002,0.004,0.006…,0.994,0.996,0.998,1} so that the quasi-PDF (f)~(x,P_(z)) and PDF f(x,μ) are represented by vectors of length 500 and the convolution kernels can be represented by matrices of dimension 500×500.
The pion quasi-PDF in $\overline{\rm MS}$ scheme at P_(z)=2 GeV is
\begin{align}
\tilde{f}(x, P_z) = (x + 3) \cdot (1-x)^3, \quad x\in (0, 1) ~.
\end{align}
# Main problem:
Given the pion quasi-PDF in $\overline{\rm MS}$ scheme at P_(z)=2 GeV as
\begin{align}
\tilde{f}(x, P_z) = (x + 3) \cdot (1-x)^3, \quad x\in (0, 1) ~.
\end{align}
Using the matching formula to derive the pion PDF f(x,μ) at μ=2 GeV in the regime x∈(0,1), note that the logarithm should be resummed using the DGLAP evolution. Evaluate the PDF at x∈{0.4,0.5,0.6}.Original LaTeX notation
# Problem setup:
The matching formula in the large-momentum effective theory (LaMET) gives
\begin{align}
f(x, \mu) = \tilde{f} (x, P_z) - \int_{0}^1 \frac{d y}{|y|} ~ C^{(1)}\left(\frac{x}{y}, \frac{\mu}{|x| P_z}\right) \tilde{f}\left(y, P_z\right) ~,
\end{align}
where $f(x, \mu)$ denotes the unpolarized parton distribution function (PDF) of the pion, characterized by the momentum fraction $x$ and the energy scale $\mu$. The corresponding quasi-PDF is denoted as $\tilde{f}$. For purposes of simplification, power corrections are ignored.
We see that the perturbative matching kernel in the $\overline{\rm MS}$ scheme is
\begin{align}
C^{(1)}\left(\xi, \frac{\mu}{|x| P_z}\right) = \frac{\alpha_s (\mu) C_F}{2 \pi} \begin{cases}\left(\frac{1+\xi^2}{1-\xi} \ln \frac{\xi}{\xi-1}+1+\frac{3}{2 \xi}\right)_{+(1)}^{[1, \infty]}-\frac{3}{2 \xi} & \xi>1 \\ \left(\frac{1+\xi^2}{1-\xi}\left[-\ln \frac{\mu^2}{4x^2 P_z^2}+\ln (\frac{1-\xi}{\xi}) \right]-\frac{\xi(1+\xi)}{1-\xi}\right)_{+(1)}^{[0,1]} & 0<\xi<1, \\
% \left(-\frac{1+\xi^2}{1-\xi} \ln \frac{-\xi}{1-\xi}-1+\frac{3}{2(1-\xi)}\right)_{+(1)}^{[-\infty, 0]}-\frac{3}{2(1-\xi)} & \xi<0
\end{cases}
\end{align}
where $\xi = x / y$. The subscript $+(1)$ indicates the plus distribution with the pole at $\xi = 1$.
The PDF satisfies the DGLAP evolution according to
\begin{align}
\begin{aligned}
\frac{d f(x, \mu)}{d \ln \mu} &= g\left(x, \mu\right), \\
g\left(x, \mu\right) &= \int_x^1 \frac{d v}{v} P\left[\frac{x}{v}, \alpha_s(\mu)\right] f\left(v, \mu\right) ~.
\end{aligned}
\end{align}
The 1-loop result of the evolution kernel $P\left[w, \alpha_s(\mu)\right]$ can be found to be
\begin{align}
P\left[w, \alpha_s(\mu)\right] = \frac{\alpha_s(\mu) C_F}{2 \pi} \left( \frac{2}{1-w} - 1 - w \right)_{+(1)} ~, w \leq 1,
% = \frac{\alpha_s(\mu) C_F}{4 \pi} \left( 4(1-w)^{-1} - 2 - 2w + 3 \delta(1 - w) \right) ~,
\end{align}
where $w = x / v$.
The 1-loop $\alpha_s$ is given by
\begin{align}
\alpha_s^{(1)}\left(\mu^2\right)=\frac{4 \pi}{\beta_0 \ln \left(\mu^2 / \Lambda_{\rm Q C D}^2\right)} ~.
\end{align}
The constants are given by
1. $C_F = \frac{4}{3}$;
2. $\beta_0 = 9$;
3. $\Lambda_{\rm Q C D} = 0.2445$ GeV.
Discretize the variables $x$, $y$, and $v$ in $x, y, v \in \{0.002, 0.004, 0.006 \dots, 0.994, 0.996, 0.998, 1\}$ so that the quasi-PDF $\tilde{f}(x, P_z)$ and PDF $f(x, \mu)$ are represented by vectors of length $500$ and the convolution kernels can be represented by matrices of dimension $500 \times 500$.
The pion quasi-PDF in $\overline{\rm MS}$ scheme at $P_z = 2$ GeV is
\begin{align}
\tilde{f}(x, P_z) = (x + 3) \cdot (1-x)^3, \quad x\in (0, 1) ~.
\end{align}
# Main problem:
Given the pion quasi-PDF in $\overline{\rm MS}$ scheme at $P_z = 2$ GeV as
\begin{align}
\tilde{f}(x, P_z) = (x + 3) \cdot (1-x)^3, \quad x\in (0, 1) ~.
\end{align}
Using the matching formula to derive the pion PDF $f(x, \mu)$ at $\mu = 2$ GeV in the regime $x\in (0, 1)$, note that the logarithm should be resummed using the DGLAP evolution. Evaluate the PDF at $x \in \{0.4, 0.5, 0.6\}$.code template
Code
def answer():
r"""
Return the values of the PDF at $x \in \{0.4, 0.5, 0.6\}$
Inputs
----------
None
Outputs
----------
f_0p4: float, the pion PDF $f(x, \mu)$ at $\mu=2$ GeV and $x=0.4$
f_0p5: float, the pion PDF $f(x, \mu)$ at $\mu=2$ GeV and $x=0.5$
f_0p6: float, the pion PDF $f(x, \mu)$ at $\mu=2$ GeV and $x=0.6$
"""
# ------------------ FILL IN YOUR RESULTS BELOW ------------------
f_0p4 = ...
f_0p5 = ...
f_0p6 = ...
# ---------------------------------------------------------------
return f_0p4, f_0p5, f_0p6Discussion
No discussion posts on this page yet. State an approach you tried, the evidence it uses, and a specific question another participant could help resolve. Use the posting template.
See answer Answer published by the source
Artifacts
Code, notes and reproducible work shared by participants. Files are served from a separate origin.
No artifacts on this page yet. Share reproducible code or notes in a contribution. State an approach you tried, the evidence it uses, and a specific question another participant could help resolve. Use the posting template.
Source and history
initial import