benchmarks.wiki / Public workspace

CritPt / Challenge_23_main / In perturbative quantum chromodynamics (QCD), one always encounters divergent…

Problem

Answer published by the source. Consult the official source to check your work against its answer.

problem description

# Problem setup: In perturbative quantum chromodynamics (QCD), one always encounters divergent Feynman diagrams in loop corrections. In the Large Momentum Effective Theory (LaMET) framework for calculating parton physics, the perturbative matching coefficient is calculated from Feynman diagrams that involves one external quark or gluon with momentum pp, which is set to be on-shell and massless, i.e., p2=0p^2=0 and \begin{align} I &= \int {d^d k\over (2\pi)^d} {N(k,p)\over D(k,p)} \,, \end{align} where dd is the dimension of the Minkowski space time with the metric gμν=diag{1,1,1,1}g^{\mu\nu}=\mbox{diag}\{1,-1,-1,-1\}. There are two types of divergences, ultraviolet (UV) and infrared (IR). Before one carries out the loop integration over kν=(k0,kx,ky,kz)k^\nu=(k^0, k^x, k^y, k^z), it is necessary to first identify the UV and IR divergences of each diagram through the superficial degree of divergence (SDD). *UV divergences.* In the limit of kνk^\nu \to \infty for each ν=0,1,2,,d\nu=0,1,2,\ldots,d, the SDD is dtd-t if N(k,p)/D(k,p)1/ktN(k,p)/D(k,p)\rightarrow 1/k^t, where tt is an integer. *IR divergences.* There are two types of IR divergences, soft and collinear divergences. - Soft divergence: In the limit of kν0k^\nu \to 0 for each ν=0,1,2,,d\nu=0,1,2,\ldots,d, N(k,p)/D(k,p)1/ktN(k,p)/D(k,p)\rightarrow 1/k^t; - Collinear divergence: In the limit of kνλpνk^\nu \to \lambda p^\nu, with 0<λ<10<\lambda< 1, N(k,p)/D(k,p)1/(p2)t/2N(k,p)/D(k,p)\rightarrow 1/(p^2)^{t/2}. One way to regulate these divergent integrals is dimensional regularization, where the space-time dimension is set to be d=42ϵd=4-2\epsilon with a continuous parameter ϵ\epsilon. It works under the assumption that the regulated results can be analytically continued in ϵ\epsilon. The UV divergences are regulated with $\epsilon_{\rm UV}>0$, and the IR divergences are regulated with $\epsilon_{\rm IR}<0$. # Main problem: The so-called sail-diagram contribution to the one-loop matrix element of the quasi-PDF in a quark state with momentum pν=(pz,0,0,pz)p^\nu=(p^z,0,0,p^z) is \begin{align} \tilde q_{\rm sail}(x,p^z,\epsilon,\mu) &=(\mu^2)^\epsilon {\exp(\epsilon\gamma_E)\over (4\pi)^\epsilon}\int_{-\infty}^\infty {dz \over 2\pi}e^{ixp^z z} \int {d^d k\over (2\pi)^d} {k^0+k^z\over k^2 (p-k)^2 (p^z-k^z)}e^{-ik^z z}. \end{align} Carry out the above integral in three intervals: x<0x < 0, 0<x<10<x<1 and x>1x>1. Then expand the final result in ϵ\epsilon to O(ϵ0)O(\epsilon^0) in the $\overline{\rm MS}$ scheme. The final results of $\tilde q_{\rm sail}(x,p^z,\epsilon,\mu)$ should be expressed as functions in xx, pzp^z, $\epsilon_{\rm UV}$ or $\epsilon_{\rm IR}$, and μ\mu.
Plain-text mathematical notation (without MathML)
# Problem setup:
In perturbative quantum chromodynamics (QCD), one always encounters divergent Feynman diagrams in loop corrections. In the Large Momentum Effective Theory (LaMET) framework for calculating parton physics, the perturbative matching coefficient is calculated from Feynman diagrams that involves one external quark or gluon with momentum p, which is set to be on-shell and massless, i.e., p²=0 and
\begin{align}
    I &= \int {d^d k\over (2\pi)^d} {N(k,p)\over D(k,p)} \,,
\end{align}
where d is the dimension of the Minkowski space time with the metric g^(μν)=diag{1,−1,−1,−1}.

There are two types of divergences, ultraviolet (UV) and infrared (IR). Before one carries out the loop integration over k^(ν)=(k⁰,k^(x),k^(y),k^(z)), it is necessary to first identify the UV and IR divergences of each diagram through the superficial degree of divergence (SDD).

*UV divergences.* In the limit of k^(ν)→∞ for each ν=0,1,2,…,d, the SDD is d−t if N(k,p)/D(k,p)→1/k^(t), where t is an integer.

*IR divergences.* There are two types of IR divergences, soft and collinear divergences.
- Soft divergence: In the limit of k^(ν)→0 for each ν=0,1,2,…,d, N(k,p)/D(k,p)→1/k^(t);
- Collinear divergence: In the limit of k^(ν)→λp^(ν), with 0<λ<1,  N(k,p)/D(k,p)→1/(p²)^(t/2).

One way to regulate these divergent integrals is dimensional regularization, where the space-time dimension is set to be d=4−2ϵ with a continuous parameter ϵ. It works under the assumption that the regulated results can be analytically continued in ϵ. The UV divergences are regulated with $\epsilon_{\rm UV}>0$, and the IR divergences are regulated with $\epsilon_{\rm IR}<0$.

# Main problem:

The so-called sail-diagram contribution to the one-loop matrix element of the quasi-PDF in a quark state with momentum p^(ν)=(p^(z),0,0,p^(z)) is

\begin{align}
    \tilde q_{\rm sail}(x,p^z,\epsilon,\mu) &=(\mu^2)^\epsilon {\exp(\epsilon\gamma_E)\over (4\pi)^\epsilon}\int_{-\infty}^\infty {dz \over 2\pi}e^{ixp^z z} \int {d^d k\over (2\pi)^d} {k^0+k^z\over k^2 (p-k)^2 (p^z-k^z)}e^{-ik^z z}.
\end{align}
Carry out the above integral in three intervals: x<0, 0<x<1 and x>1. Then expand the final result in ϵ to O(ϵ⁰) in the $\overline{\rm MS}$ scheme. The final results of $\tilde q_{\rm sail}(x,p^z,\epsilon,\mu)$ should be expressed as functions in x, p^(z), $\epsilon_{\rm UV}$ or $\epsilon_{\rm IR}$, and μ.
Original LaTeX notation
# Problem setup:
In perturbative quantum chromodynamics (QCD), one always encounters divergent Feynman diagrams in loop corrections. In the Large Momentum Effective Theory (LaMET) framework for calculating parton physics, the perturbative matching coefficient is calculated from Feynman diagrams that involves one external quark or gluon with momentum $p$, which is set to be on-shell and massless, i.e., $p^2=0$ and
\begin{align}
    I &= \int {d^d k\over (2\pi)^d} {N(k,p)\over D(k,p)} \,,
\end{align}
where $d$ is the dimension of the Minkowski space time with the metric $g^{\mu\nu}=\mbox{diag}\{1,-1,-1,-1\}$.

There are two types of divergences, ultraviolet (UV) and infrared (IR). Before one carries out the loop integration over $k^\nu=(k^0, k^x, k^y, k^z)$, it is necessary to first identify the UV and IR divergences of each diagram through the superficial degree of divergence (SDD).

*UV divergences.* In the limit of $k^\nu \to \infty$ for each $\nu=0,1,2,\ldots,d$, the SDD is $d-t$ if $N(k,p)/D(k,p)\rightarrow 1/k^t$, where $t$ is an integer.

*IR divergences.* There are two types of IR divergences, soft and collinear divergences.
- Soft divergence: In the limit of $k^\nu \to 0$ for each $\nu=0,1,2,\ldots,d$, $N(k,p)/D(k,p)\rightarrow 1/k^t$;
- Collinear divergence: In the limit of $k^\nu \to \lambda p^\nu$, with $0<\lambda< 1$,  $N(k,p)/D(k,p)\rightarrow 1/(p^2)^{t/2}$.

One way to regulate these divergent integrals is dimensional regularization, where the space-time dimension is set to be $d=4-2\epsilon$ with a continuous parameter $\epsilon$. It works under the assumption that the regulated results can be analytically continued in $\epsilon$. The UV divergences are regulated with $\epsilon_{\rm UV}>0$, and the IR divergences are regulated with $\epsilon_{\rm IR}<0$.

# Main problem:

The so-called sail-diagram contribution to the one-loop matrix element of the quasi-PDF in a quark state with momentum $p^\nu=(p^z,0,0,p^z)$ is

\begin{align}
    \tilde q_{\rm sail}(x,p^z,\epsilon,\mu) &=(\mu^2)^\epsilon {\exp(\epsilon\gamma_E)\over (4\pi)^\epsilon}\int_{-\infty}^\infty {dz \over 2\pi}e^{ixp^z z} \int {d^d k\over (2\pi)^d} {k^0+k^z\over k^2 (p-k)^2 (p^z-k^z)}e^{-ik^z z}.
\end{align}
Carry out the above integral in three intervals: $x < 0$, $0<x<1$ and $x>1$. Then expand the final result in $\epsilon$ to $O(\epsilon^0)$ in the $\overline{\rm MS}$ scheme. The final results of $\tilde q_{\rm sail}(x,p^z,\epsilon,\mu)$ should be expressed as functions in $x$, $p^z$, $\epsilon_{\rm UV}$ or $\epsilon_{\rm IR}$, and $\mu$.

code template

Code

import sympy as sp

x, p_z, epsilon_UV, epsilon_IR, mu = sp.symbols('x p_z epsilon_UV epsilon_IR mu')

def answer(x, p_z, epsilon_UV, epsilon_IR, mu):
    r"""
    Return the expressions of $\tilde q_{\rm sail}(x,p^z,\epsilon,\mu)$ in Sympy format.

    Inputs
    ----------
    x: sympy.Symbol, longitudinal momentum fraction $x$
    p_z: sympy.Symbol, longitudinal momentum $p^z$
    epsilon_UV: sympy.Symbol, dimensional–regularization parameter for UV divergences, $\epsilon_{\rm UV}$
    epsilon_IR: sympy.Symbol, dimensional–regularization parameter for IR divergences, $\epsilon_{\rm IR}$
    mu: sympy.Symbol, renormalization scale $\mu$

    Outputs
    ----------
    expr_lt0:  sympy.Expr,  $\tilde q_{\rm sail}(x,p^z,\epsilon,\mu)$ for $x < 0$, to $O(\epsilon_{\rm UV}^0)$
    expr_mid:  sympy.Expr,  $\tilde q_{\rm sail}(x,p^z,\epsilon,\mu)$ for $0<x<1$, to $O(\epsilon_{\rm UV}^0)$
    expr_gt1:  sympy.Expr,  $\tilde q_{\rm sail}(x,p^z,\epsilon,\mu)$ for $x > 1$, to $O(\epsilon_{\rm UV}^0)$
    """

    # ------------------ FILL IN YOUR RESULTS BELOW ------------------
    expr_lt0 = ...
    expr_mid = ...
    expr_gt1 = ...
    # ---------------------------------------------------------------

    return expr_lt0, expr_mid, expr_gt1

Discussion

Discussion

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

Official source

initial import