Source author record

Andreas Schäfer

Andreas Schäfer appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

74works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

74 published item(s)

preprint2026arXiv

Total Gluon Helicity Contribution to the Proton Spin from Lattice QCD

We report a state-of-the-art lattice QCD calculation of the total gluon helicity contribution to the proton spin, $ΔG$. The calculation is done on ensembles with three different lattice spacings $a=\{0.08, 0.09, 0.11\}$ fm. By employing distillation and momentum smearing for proton external states, we extract the bare matrix elements of the topological current $K^μ$ using 5-HYP smeared Coulomb gauge fixing configurations. Furthermore, we apply a non-perturbative $\mathrm{RI/MOM}$ renormalization scheme augmented by the Cluster Decomposition Error Reduction (CDER) technique to determine the renormalization constants of $K^μ$. The results obtained from different components $K^{t,i}$ (with $i$ being the direction of proton momentum or polarization) are consistent with Lorentz covariance within uncertainties. After extrapolating to the continuum limit, $ΔG$ is found to be $ΔG = 0.231(17)^{\mathrm{sta.}}(44)^{\mathrm{sym.}}$ at the $\overline{\mathrm{MS}}$ scale $μ^2=10\ \mathrm{GeV}^2$, which constitutes approximately $46(9)\%$ of the proton spin.

preprint2022arXiv

BITES: Balanced Individual Treatment Effect for Survival data

Estimating the effects of interventions on patient outcome is one of the key aspects of personalized medicine. Their inference is often challenged by the fact that the training data comprises only the outcome for the administered treatment, and not for alternative treatments (the so-called counterfactual outcomes). Several methods were suggested for this scenario based on observational data, i.e.~data where the intervention was not applied randomly, for both continuous and binary outcome variables. However, patient outcome is often recorded in terms of time-to-event data, comprising right-censored event times if an event does not occur within the observation period. Albeit their enormous importance, time-to-event data is rarely used for treatment optimization. We suggest an approach named BITES (Balanced Individual Treatment Effect for Survival data), which combines a treatment-specific semi-parametric Cox loss with a treatment-balanced deep neural network; i.e.~we regularize differences between treated and non-treated patients using Integral Probability Metrics (IPM). We show in simulation studies that this approach outperforms the state of the art. Further, we demonstrate in an application to a cohort of breast cancer patients that hormone treatment can be optimized based on six routine parameters. We successfully validated this finding in an independent cohort. BITES is provided as an easy-to-use python implementation.

preprint2022arXiv

Confinement/deconfinement transition in the D0-brane matrix model -- A signature of M-theory?

We study the confinement/deconfinement transition in the D0-brane matrix model (often called the BFSS matrix model) and its one-parameter deformation (the BMN matrix model) numerically by lattice Monte Carlo simulations. Our results confirm general expectations from the dual string/M-theory picture for strong coupling. In particular, we observe the confined phase in the BFSS matrix model, which is a nontrivial consequence of the M-theory picture. We suggest that these models provide us with an ideal framework to study the Schwarzschild black hole, M-theory, and furthermore, the parameter region of the phase transition between type IIA superstring theory and M-theory. A detailed study of M-theory via lattice Monte Carlo simulations of the D0-brane matrix model might be doable with much smaller computational resources than previously expected.

preprint2022arXiv

Mellin moments of spin dependent and independent PDFs of the pion and rho meson

We compute the second moments of pion and rho parton distribution functions (PDFs) in lattice QCD with $N_f = 2+1$ flavors of improved Wilson fermions. We determine both singlet and non-singlet flavor combinations and, for the first time, take disconnected contributions fully into account. In the case of the rho, we also calculate the additional contribution arising from the $b_1$ structure function. The numerical analysis includes 26 ensembles, mainly generated by the CLS effort, with pion masses ranging from 420MeV down to 214MeV and with 5 different lattice spacings in the range of 0.1fm to 0.05fm. This enables us to take the continuum limit, as well as to resolve the quark mass dependencies reliably. Additionally we discuss the contaminations of rho correlation functions by two-pion states.

preprint2022arXiv

Nonperturbative Determination of Collins-Soper Kernel from Quasi Transverse-Momentum Dependent Wave Functions

In the framework of large-momentum effective theory at one-loop matching accuracy, we perform a lattice calculation of the Collins-Soper kernel which governs the rapidity evolution of transverse-momentum-dependent (TMD) distributions. We first obtain the quasi TMD wave functions at three different meson momenta on a lattice with valence clover quarks on a dynamical HISQ sea and lattice spacing $a=0.12$~fm from MILC, and renormalize the pertinent linear divergences using Wilson loops. Through one-loop matching to the light-cone wave functions, we determine the Collins-Soper kernel with transverse separation up to 0.6~fm. We study the systematic uncertainties from operator mixing and scale dependence, as well as the impact from higher power corrections. Our results potentially allow for a determination of the soft function and other transverse-momentum dependent quantities at one-loop accuracy.

preprint2022arXiv

T-odd proton-helicity asymmetry in semi-inclusive DIS in perturbative QCD

We compute the single-spin asymmetry $A_{UL}$ in semi-inclusive deep-inelastic scattering of unpolarized leptons and longitudinally polarized protons at large transverse momentum of the produced hadron. Our calculation is performed in collinear factorization at the lowest order of QCD perturbation theory. For photon exchange the asymmetry is T-odd and receives contributions from the interference of the tree level and one-loop absorptive amplitudes. We consider the behavior of the spin asymmetry at low transverse momentum where contact to the formalism based on transverse-momentum dependent distribution functions can be made. We also present some phenomenological results relevant for the COMPASS and HERMES experiments and the future Electron-Ion Collider.

preprint2021arXiv

A Hybrid Renormalization Scheme for Quasi Light-Front Correlations in Large-Momentum Effective Theory

In large-momentum effective theory (LaMET), calculating parton physics starts from calculating coordinate-space-$z$ correlation functions $\tilde h(z, a,P^z)$ in a hadron of momentum $P^z$ in lattice QCD. Such correlation functions involve both linear and logarithmic divergences in lattice spacing $a$, and thus need to be properly renormalized. We introduce a hybrid renormalization procedure to match these lattice correlations to those in the continuum $\overline{\rm MS}$ scheme, without introducing extra non-perturbative effects at large $z$. We analyze the effect of ${\cal O}(Λ_{\rm QCD})$ ambiguity in the Wilson line self-energy subtraction involved in this hybrid scheme. To obtain the momentum-space distributions, we recommend to extrapolate the lattice data to the asymptotic $z$-region using the generic properties of the coordinate space correlations at moderate and large $P^z$, respectively.

preprint2021arXiv

Double parton distributions in the pion from lattice QCD

We perform a lattice study of double parton distributions in the pion, using the relationship between their Mellin moments and pion matrix elements of two local currents. A good statistical signal is obtained for almost all relevant Wick contractions. We investigate correlations in the spatial distribution of two partons in the pion, as well as correlations involving the parton polarisation. The patterns we observe depend significantly on the quark mass. We investigate the assumption that double parton distributions approximately factorise into a convolution of single parton distributions.

preprint2021arXiv

Holographic Kolmogorov-Sinai entropy and the quantum Lyapunov spectrum

In classical chaotic systems the entropy, averaged over initial phase space distributions, follows an universal behavior. While approaching thermal equilibrium it passes through a stage where it grows linearly, while the growth rate, the Kolmogorov-Sinai entropy, is given by the sum over all positive Lyapunov exponents. A natural question is whether a similar relation is valid for quantum systems. We argue that the Maldacena-Shenker-Stanford bound on quantum Lyapunov exponents $λ$ implies that the upper bound on the growth rate of the entropy, averaged over states in Hilbert space that evolve towards a thermal state with temperature $T$ and entropy $S_{eq}$, should be given by $S_{eq} πT =\sum_{λ>0}2 πT$. Strongly coupled, large $N$ theories with black hole duals should saturate the bound. By studying a large number of isotropization processes of random, spatially homogeneous, far from equilibrium initial states in large $N$, $\mathcal{N}=4$ Super Yang Mills theory at strong coupling and computing the ensemble averaged growth rate of the dual black hole's apparent horizon area, we find both an analogous behavior as in classical chaotic systems and numerical evidence that the conjectured bound on averaged entropy growth is saturated granted that the Lyapunov exponents are degenerate $λ= \pm 2 πT$. This fits to the behavior of classical systems with plus/minus symmetric Lyapunov spectra, a symmetry which implies the validity of Liouville's theorem.

preprint2021arXiv

Partial Deconfinement at Strong Coupling on the Lattice

We provide evidence for partial deconfinement -- the deconfinement of a SU($M$) subgroup of the SU($N$) gauge group -- by using lattice Monte Carlo simulations. We take matrix models as concrete examples. By appropriately fixing the gauge, we observe that the $M\times M$ submatrices deconfine. This gives direct evidence for partial deconfinement at strong coupling. We discuss the applications to QCD and holography.

preprint2020arXiv

Nucleon axial structure from lattice QCD

We present a new analysis method that allows one to understand and model excited state contributions in observables that are dominated by a pion pole. We apply this method to extract axial and (induced) pseudoscalar nucleon isovector form factors, which satisfy the constraints due to the partial conservation of the axial current up to expected discretization effects. Effective field theory predicts that the leading contribution to the (induced) pseudoscalar form factor originates from an exchange of a virtual pion, and thus exhibits pion pole dominance. Using our new method, we can recover this behavior directly from lattice data. The numerical analysis is based on a large set of ensembles generated by the CLS effort, including physical pion masses, large volumes (with up to $96^3 \times 192$ sites and $L m_π= 6.4$), and lattice spacings down to $0.039 \, \text{fm}$, which allows us to take all the relevant limits. We find that some observables are much more sensitive to the choice of parametrization of the form factors than others. On the one hand, the $z$-expansion leads to significantly smaller values for the axial dipole mass than the dipole ansatz ($M_A^{\text{$z$-exp}}=1.02(10) \, \text{GeV}$ versus $M_A^{\text{dipole}} = 1.31(8) \, \text{GeV}$). On the other hand, we find that the result for the induced pseudoscalar coupling at the muon capture point is almost independent of the choice of parametrization ($g_P^{\star \ \text{$z$-exp}} = 8.68(45)$ and $g_P^{\star \ \text{dipole}} = 8.30(24)$), and is in good agreement with both, chiral perturbation theory predictions and experimental measurement via ordinary muon capture. We also determine the axial coupling constant $g_A$.

preprint2020arXiv

Phenomenological implications of asymmetric $AdS_5$ shockwave collision studies for heavy ion physics

This paper discusses possible phenomenological implications for p+A and A+A collisions of the results of recent numerical AdS/CFT calculations examining asymmetric collisions of planar shocks. In view of the extreme Lorentz contraction, we model highly relativistic heavy ion collisions (HICs) as a superposition of collisions between many near-independent transverse "pixels" with differing incident longitudinal momenta. It was found that also for asymmetric collisions the hydrodynamization time is in good approximation a proper time, just like for symmetric collisions, depending on the geometric mean of the longitudinally integrated energy densities of the incident projectiles. For realistic collisions with fluctuations in the initial energy densities, these results imply a substantial increase in the hydrodynamization time for highly asymmetric pixels. However, even in this case the local hydrodynamization time still is significantly smaller than perturbative results for the thermalization time.

preprint2020arXiv

Transverse momentum dependent factorization for lattice observables

Using soft collinear effective field theory, we derive the factorization theorem for the quasi-transverse-momentum-dependent (quasi-TMD) operator. We check the factorization theorem at one-loop level and compute the corresponding coefficient function and anomalous dimensions. The factorized expression is built from the physical TMD distribution, and a nonperturbative lattice related factor. We demonstrate that lattice related functions cancel in appropriately constructed ratios. These ratios could be used to explore various properties of TMD distributions, for instance, the nonperturbative evolution kernel. A discussion of such ratios and the related continuum properties of TMDs is presented.

preprint2019arXiv

Unpolarized isovector quark distribution function from Lattice QCD: A systematic analysis of renormalization and matching

We present a detailed Lattice QCD study of the unpolarized isovector quark Parton Distribution Function (PDF) using large-momentum effective theory framework. We choose a quasi-PDF defined by a spatial correlator which is free from mixing with other operators of the same dimension. In the lattice simulation, we use a Gaussian-momentum-smeared source at $M_π=356$ MeV and $P_z \in \{1.8,2.3\}$ GeV. To control the systematics associated with the excited states, we explore {five different source-sink separations}. The nonperturbative renormalization is conducted in a regularization-independent momentum subtraction scheme, and the matching between the renormalized quasi-PDF and $\bar{\rm MS}$ PDF is calculated based on perturbative QCD up to one-loop order. Systematic errors due to renormalization and perturbative matching are also analyzed in detail. Our results for lightcone PDF are in reasonable agreement with the latest phenomenological analysis.

preprint2016arXiv

$ρ$ and $K^*$ resonances on the lattice at nearly physical quark masses and $N_f=2$

Working with a pion mass $m_π\approx 150$ MeV, we study $ππ$ and $Kπ$ scattering using two flavours of non-perturbatively improved Wilson fermions at a lattice spacing $a\approx 0.071$ fm. Employing two lattice volumes with linear spatial extents of $N_s=48$ and $N_s=64$ points and moving frames, we extract the phase shifts for p-wave $ππ$ and $Kπ$ scattering near the $ρ$ and $K^*$ resonances.Comparing our results to those of previous lattice studies, that used pion masses ranging from about 200 MeV up to 470 MeV, we find that the coupling $g_{ρππ}$ appears to be remarkably constant as a function of $m_π$.

preprint2016arXiv

Direct determinations of the nucleon and pion $σ$ terms at nearly physical quark masses

We present a high statistics study of the pion and nucleon light and strange quark sigma terms using $N_f=2$ dynamical non-perturbatively improved clover fermions with a range of pion masses down to $m_π\sim 150$ MeV and several volumes, $Lm_π=3.4$ up to $6.7$, and lattice spacings, $a=0.06-0.08$ fm, enabling a study of finite volume and discretisation effects for $m_π\gtrsim 260$ MeV. Systematics are found to be reasonably under control. For the nucleon we obtain $σ_{πN}=35(6)$ MeV and $σ_s=35(12)$ MeV, or equivalently in terms of the quark fractions, $f_{T_u}=0.021(4)$, $f_{T_d}=0.016(4)$ and $f_{T_s}=0.037(13)$, where the errors include estimates of both the systematic and statistical uncertainties. These values, together with perturbative matching in the heavy quark limit, lead to $f_{T_c}=0.075(4)$, $f_{T_b}=0.072(2)$ and $f_{T_t}=0.070(1)$. In addition, through the use of the (inverse) Feynman-Hellmann theorem our results for $σ_{πN}$ are shown to be consistent with the nucleon masses determined in the analysis. For the pion we implement a method which greatly reduces excited state contamination to the scalar matrix elements from states travelling across the temporal boundary. This enables us to demonstrate the Gell-Mann-Oakes-Renner expectation $σ_π=m_π/2$ over our range of pion masses.

preprint2016arXiv

Light-cone distribution amplitudes of the baryon octet

We present results of the first ab initio lattice QCD calculation of the normalization constants and first moments of the leading twist distribution amplitudes of the full baryon octet, corresponding to the small transverse distance limit of the associated S-wave light-cone wave functions. The P-wave (higher twist) normalization constants are evaluated as well. The calculation is done using $N_f=2+1$ flavors of dynamical (clover) fermions on lattices of different volumes and pion masses down to 222 MeV. Significant SU(3) flavor symmetry violation effects in the shape of the distribution amplitudes are observed.

preprint2016arXiv

Novel quark smearing for hadrons with high momenta in lattice QCD

Hadrons in lattice QCD are usually created employing smeared interpolators. We introduce a new quark smearing that allows us to maintain small statistical errors and good overlaps of hadronic wavefunctions with the respective ground states, also at high spatial momenta. The method is successfully tested for the pion and the nucleon at a pion mass $m_π\approx 295$ MeV and momenta as high as 2.8 GeV. We compare the results obtained to dispersion relations and suggest further optimizations.

preprint2016arXiv

Nucleon generalized form factors from lattice QCD with nearly physical quark masses

We determine generalized form factors of the nucleon from lattice simulations with $N_f = 2$ mass-degenerate non-perturbatively improved Wilson-Sheikholeslami-Wohlert fermions down to a pion mass of 150 MeV. We also present the resulting isovector quark angular momentum. Possible excited-state contaminations are investigated with correlated simultaneous fits.

preprint2016arXiv

The decays $Λ_{b,c}\to N^*\, l\,ν$ in QCD

We present an exploratory study of the $Λ_{c,b}\to N^*$-form factors and the semileptonic decay width within the framework of light-cone sum rules. We use two different methods and two different interpolating currents for the $Λ_{c,b}$. As interpolating currents we choose an axial-vector like and a pseudoscalar like current. Our results show that the procedure of eliminating negative parity partners is not well suited for the case at hand and that the second approach using structures with highest powers of $p_+$ with an axial-vector like interpolating current gives the most reliable results. Our predictions are based on the models obtained in \cite{Anikin:2015ita,Braun:2014wpa}. The largest uncertainty comes from the uncertainty of the twist 4 parameters $η_{10},\,η_{11}$ and we take the spread between the two models in \cite{Anikin:2015ita} as a measure for this. We get \begin{eqnarray} Γ(Λ_b\to N^*(1535) l ν)&=&\left(0.0058^{+0.0010}_{-0.0009}\right)\cdot\left(\dfrac{V_{ub}}{3.5\cdot10^{-3}}\right)^2,\quad\mbox{LCSR(1)} \nonumber Γ(Λ_b\to N^*(1535) l ν)&=&\left(0.00070^{+0.00012}_{-0.00011}\right)\cdot\left(\dfrac{V_{ub}}{3.5\cdot10^{-3}}\right)^2,\quad\mbox{LCSR(2)} \nonumber Γ(Λ_c\to N^*(1535) l ν)&=&\left(0.0064^{+0.0012}_{-0.0011}\right)\cdot \left(\dfrac{V_{cd}}{0.225}\right)^2,\quad\mbox{LCSR(1)}\nonumber Γ(Λ_c\to N^*(1535) l ν)&=&\left(0.00077^{+0.00016}_{-0.00014}\right)\cdot \left(\dfrac{V_{cd}}{0.225}\right)^2,\quad\mbox{LCSR(2)}\nonumber \end{eqnarray} as predictions for the respective decay widths, where LCSR(1) and LCSR(2) refer to the two different models of the distribution amplitudes of the $N^*$. It is seen that even a rough measurement of these decays will greatly help to discriminate different models.

preprint2015arXiv

Center clusters in full QCD at finite temperature and background magnetic field

We study the center structure of full dynamical QCD at finite temperatures and nonzero values of the background magnetic field using continuum extrapolated lattice data. We concentrate on two particular observables characterizing center clusters: their fractality and the probability for percolation. For temperatures below and around the transition region, the fractal dimension is found to be significantly smaller than three, leading to a vanishing mean free path inside the cluster structure. This finding might be relevant for center symmetry-based models of heavy-ion collisions. In addition, the percolation probability is employed to define the transition temperature and to map out the QCD phase diagram in the magnetic field-temperature plane.

preprint2015arXiv

Model-independent calculation of $SU(3)_f$ violation in baryon octet light-cone distribution amplitudes

In this work we present a minimal parametrization of the light-cone distribution amplitudes of the baryon octet including higher twist contributions. Simultaneously we obtain the quark mass dependence of the amplitudes at leading one-loop accuracy by the use of three-flavor baryon chiral perturbation theory (BChPT), which automatically yields model-independent results for the leading $SU(3)$ flavor breaking effects. For that purpose we have constructed the nonlocal light-cone three-quark operators in terms of baryon octet and meson fields and have carried out a next-to-leading order BChPT calculation. We were able to find a minimal set of distribution amplitudes (DAs) that do not mix under chiral extrapolation towards the physical point and naturally embed the $Λ$ baryon. Additionally they are chosen in such a way that all DAs of a certain symmetry class have a similar quark mass dependence (independent of the twist of the corresponding amplitude), which allows for a compact presentation. The results are well-suited for the extrapolation of lattice data and for model building.

preprint2015arXiv

Nucleon isovector couplings from $N_f=2$ lattice QCD

We compute the axial, scalar, tensor and pseudoscalar isovector couplings of the nucleon as well as the induced tensor and pseudoscalar charges in lattice simulations with $N_f=2$ mass-degenerate non-perturbatively improved Wilson-Sheikholeslami-Wohlert fermions. The simulations are carried out down to a pion mass of 150 MeV and linear spatial lattice extents of up to 4.6 fm at three different lattice spacings ranging from approximately 0.08 fm to 0.06 fm. Possible excited state contamination is carefully investigated and finite volume effects are studied. The couplings, determined at these lattice spacings, are extrapolated to the physical pion mass. In this limit we find agreement with experimental results, where these exist, with the exception of the magnetic moment. A proper continuum limit could not be performed, due to our limited range of lattice constants, but no significant lattice spacing dependence is detected. Upper limits on discretization effects are estimated and these dominate the error budget.

preprint2015arXiv

One-Loop Matching for Generalized Parton Distributions

We present the one-loop matching condition for the unpolarized and polarized generalized quark distributions in the non-singlet case. The matching condition links the quasi distributions defined in terms of spacelike correlators at finite nucleon momentum to the light-cone distributions, and is useful for extracting the latter from the former in a lattice QCD calculation. Our results show that at one-loop and leading power accuracy the matching for the light-cone generalized quark distribution $H$ ($\tilde H$) is non-trivial, whereas no matching is required for $E$ ($\tilde E$). Therefore, $E$ ($\tilde E$) can be smoothly approached by its quasi counterpart in the large momentum limit. We also present the matching for the distribution amplitude of the pion.

preprint2015arXiv

Spin Decomposition of Electron in QED

We perform a systematic study on the spin decomposition of an electron in QED at one-loop order. It is found that the electron orbital angular momentum defined in Jaffe-Manohar and Ji spin sum rules agrees with each other, and the so-called potential angular momentum vanishes at this order. The calculations are performed in both dimensional regularization and Pauli-Villars regularization for the ultraviolet divergences, and they lead to consistent results. We further investigate the calculations in terms of light-front wave functions, and find a missing contribution from the instantaneous interaction in light-front quantization. This clarifies the confusing issues raised recently in the literature on the spin decomposition of an electron, and will help to consolidate the spin physics program for nucleons in QCD.

preprint2015arXiv

Threshold Resummation for Polarized High-$p_T$ Hadron Production at COMPASS

We study the cross section for the photoproduction process $γN\rightarrow h X$ where the incident photon and nucleon are longitudinally polarized and a hadron $h$ is observed at high transverse momentum. Specifically, we address the "direct" part of the cross section, for which the photon interacts in a pointlike way. For this contribution we perform an all-order resummation of logarithmic threshold corrections generated by soft or collinear gluon emission to next-to-leading logarithmic accuracy. We present phenomenological results relevant for the COMPASS experiment and compare to recent COMPASS data.

preprint2015arXiv

Transverse single spin asymmetry in direct photon production in polarized pA collisions

We study the transverse single spin asymmetry in direct photon production in pA collisions with incoming protons being transversely polarized. To facilitate the calculation, we formulate a hybrid approach in which the nucleus is treated in the Color Glass Condensate (CGC) framework while the collinear twist-3 formalism is applied on the proton side. It has been found that an additional term which arises from color entanglement shows up in the spin dependent differential cross section. The fact that this additional term is perturbatively calculable allows us to quantitatively study color entanglement effects.

preprint2014arXiv

A lattice study of the jet quenching parameter

We present a first-principle computation of the jet quenching parameter, which describes the momentum broadening of a high-energy parton moving through the deconfined state of QCD matter at high temperature. Following an idea originally proposed by Caron-Huot, we explain how one can evaluate the soft contribution to the collision kernel characterizing this real-time phenomenon, analyzing certain gauge-invariant operators in a dimensionally reduced effective theory (electrostatic QCD), which can be studied non-perturbatively via simulations on a Euclidean lattice. Our high-precision numerical computations at two different temperatures indicate that soft contributions to the jet quenching parameter are large. After discussing the systematic uncertainties involved, we present a quantitative estimate for the jet quenching parameter in the temperature range accessible at heavy-ion colliders, and compare it to results from phenomenological models as well as to strong-coupling computations based on the holographic correspondence.

preprint2014arXiv

Color entanglement for $γ$-jet production in polarized pA collisions

A more reliable treatment of transverse momentum dependent physics and in particular transverse single spin asymmetries is urgently required, e.g. for polarized pA physics including novel effects like color entanglement. We argue that the measurement of azimuthal angular correlations of photons and jets produced in pA collisions provides a direct access to the novel gluon distribution $G_4$ that enters into many such processes.

preprint2014arXiv

First moments of nucleon generalized parton distributions in chiral perturbation theory at full one-loop order

In this paper we present chiral extrapolation formulas for the generalized form factors connected to the first moments of chiral-even generalized parton distributions. For the description of QCD at low energies we employ the framework of two-flavor covariant baryon chiral perturbation theory at full one-loop accuracy. Our results are well suited for the extrapolation of lattice QCD data to physical quark masses.

preprint2014arXiv

Investigating jet quenching on the lattice

Due to the dynamical, real-time, nature of the phenomenon, the study of jet quenching via lattice QCD simulations is not straightforward. In this contribution, however, we show how one can extract information about the momentum broadening of a hard parton moving in the quark-gluon plasma, from lattice calculations. After discussing the basic idea (originally proposed by Caron-Huot), we present a recent study, in which we estimated the jet quenching parameter non-perturbatively, from the lattice evaluation of a particular set of gauge-invariant operators.

preprint2014arXiv

Jet quenching from the lattice

We present a lattice study of the momentum broadening experienced by a hard parton in the quark-gluon plasma. In particular, the contributions to this real-time phenomenon from soft modes are extracted from a set of gauge-invariant operators in a dimensionally reduced effective theory (electrostatic QCD), which can be simulated on a Euclidean lattice. At the temperatures accessible to present experiments, the soft contributions to the jet quenching parameter are found to be quite large. We compare our results to phenomenological models and to holographic computations.

preprint2014arXiv

Spectral functions of charmonium with 2+1 flavours of dynamical quarks

Finite temperature charmonium spectral functions in the pseudoscalar(PS) and vector(V) channels are studied in lattice QCD with 2+1 flavours of dynamical Wilson quarks, on fine isotropic lattices (with a lattice spacing of 0.057fm), with a non-physical pion mass of 545MeV. The highest temperature studied is approximately 1.4Tc. Up to this temperature no significant variation of the spectral function is seen in the PS channel. The V channel shows some temperature dependence, which seems to be consistent with a temperature dependent low frequency peak related to heavy quark transport, plus a temperature independent term at omega > 0. These results are in accord with previous calculations using the quenched approximation.

preprint2014arXiv

The moment $\langle x\rangle_{u-d}$ of the nucleon from $N_f=2$ lattice QCD down to nearly physical quark masses

We present an update of our analysis [1] which includes additional ensembles at different quark masses, lattice spacings and volumes, all with high statistics. We use $N_f=2$ mass-degenerate quark flavours, employing the non-perturbatively improved clover action. The lattice matrix elements are converted to the $\overline{\rm MS}$ scheme via renormalization factors determined non-perturbatively in the RI$^\prime$-MOM scheme. We have systematically investigated excited state contributions, in particular, at the smallest, near physical, pion mass. While our results~(with much increased precision) are consistent with Ref.~[1], comparing with previous determinations we find that excited state contributions can be significant if the quark smearing is not suitably optimized, in agreement with other recent studies. The difference with respect to the value for $\langle x\rangle_{u-d}$ extracted from experimental data is reduced but not resolved. Using lattice sizes in the range $L m_π\sim 3.4-6.7$, no significant finite volume effects were observed. Performing a controlled continuum limit that may remove the discrepancy will require simulations at lattice spacings $a< 0.06$ fm.

preprint2013arXiv

Azimuthal asymmetries for quark pair production in pA collisions

We study the azimuthal asymmetries for quark pair production in proton-nucleus collisions using a hybrid approach in which the nucleus is treated in the Color Glass Condensate (CGC) framework while the Lipatov approximation is applied on the proton side. Our treatment goes beyond the large N_c limit. We particularly focus on the so-called correlation limit where the imbalance of the transverse momentum of the quark pair is much smaller than the out-going individual quark transverse momenta. In this kinematic region, a matching between the hybrid approach and a factorization in terms of transverse momentum dependent parton distributions(TMDs) has been found. It is shown which of the various unpolarized and linearly polarized gluon TMD distributions contribute to $\cos 2 ϕ$ and $\cos 4 ϕ$ modulations of quark pair production.

preprint2013arXiv

Chiral behavior of vector meson self energies

We employ a chiral Lagrangian framework with three dynamical flavors to calculate the masses of the lowest-lying vector mesons to one loop accuracy, and use the resulting formulae to extrapolate recent QCDSF lattice data on vector meson masses and mass ratios to the physical point. Our representation for the vector meson self energies also enables us to discuss loop corrections to the $ωϕ$ mixing amplitude.

preprint2013arXiv

Excited light and strange hadrons from the lattice with two Chirally Improved quarks

Results for excited light and strange hadrons from the lattice with two flavors of Chirally Improved sea quarks are presented. We perform simulations at several values of the pion mass ranging from 250 to 600 MeV and extrapolate to the physical pion mass. The variational method is applied to extract excited energy levels but also to discuss the content of the states. Among others, we explore the flavor singlet/octet content of Lambda states. In general, our results agree well with experiment, in particular we confirm the Lambda(1405) and its dominant flavor singlet structure.

preprint2013arXiv

Jet quenching in a strongly interacting plasma - A lattice approach

The phenomenon of jet quenching, related to the momentum broadening of a high-energy parton, provides important experimental evidence for the production of a strongly coupled, deconfined medium in heavy-ion collisions. Its theoretical description has been addressed in a number of works, both perturbatively and non-perturbatively (using the gauge-gravity duality). In this contribution, following a proposal by Caron-Huot, we discuss a novel approach to this problem, enabling one to extract non-perturbative information on this real-time phenomenon from simulations on a Euclidean lattice.

preprint2013arXiv

Model Independent Evolution of Transverse Momentum Dependent Distribution Functions (TMDs) at NNLL

We discuss the evolution of the eight leading twist transverse momentum dependent parton distribution functions, which turns out to be universal and spin independent. By using the highest order perturbatively calculable ingredients at our disposal, we perform the resummation of the large logarithms that appear in the evolution kernel of transverse momentum distributions up to next-to-next-to-leading logarithms (NNLL), thus obtaining an expression for the kernel with highly reduced model dependence. Our results can also be obtained using the standard CSS approach when a particular choice of the $b^*$ prescription is used. In this sense, and while restricted to the perturbative domain of applicability, we consider our results as a "prediction" of the correct value of $b_{\rm max}$ which is very close to $1.5 {\rm GeV}^{-1}$. We explore under which kinematical conditions the effects of the non-perturbative region are negligible, and hence the evolution of transverse momentum distributions can be applied in a model independent way. The application of the kernel is illustrated by considering the unpolarized transverse momentum dependent parton distribution function and the Sivers function.

preprint2013arXiv

Momentum broadening of partons on the light cone from the lattice

The jet quenching parameter describes the momentum broadening of a high-energy parton moving through the quark-gluon plasma. Following an approach originally proposed by Caron-Huot, we discuss how one can extract information on the collision kernel associated with the parton momentum broadening, from the analysis of certain gauge-invariant operators in dimensionally reduced effective theories, and present numerical results from a lattice study.

preprint2013arXiv

Nucleon structure from stochastic estimators

Using stochastic estimators for connected meson and baryon three-point functions has successfully been tried in the past years. Compared to the standard sequential source method we trade the freedom to compute the current-to-sink propagator independently of the hadron sink for additional stochastic noise in our observables. In the case of the nucleon we can use this freedom to compute many different sink-momentum/polarization combinations, which grants access to more virtualities. We will present preliminary results on the scalar, electro-magnetic and axial form factors of the nucleon in $N_f=2+1$ lattice QCD and contrast the performance of the stochastic method to the sequential source method. We find the stochastic method to be competitive in terms of errors at fixed cost.

preprint2013arXiv

On the process dependent nuclear $k_\perp$ broadening effect

We study the process dependent nuclear $k_\perp$ broadening effect by employing the transverse momentum dependent(TMD) factorization approach in combination with the Mclerran-Venugopalan(MV) model. More specifically, we investigate how the parton transverse momentum distributions are affected by the process dependent gauge links in cold nuclear matter. In particular, our analysis also applies to the polarized cases including the nuclear quark Boer-Mulders function and the linearly polarized gluon distribution. Our main focus is on the nuclear TMDs at intermediate or large x.

preprint2013arXiv

Transverse single spin asymmetry in hadronic $η_{c,b}$ production

We study the transverse single spin asymmetry in $η_{c,b}$ production in polarized hadron collisions, employing the collinear twist-3 approach in combination with the color singlet model. Our main focus lies on the contribution from the twist-3 Efremov-Teryaev-Qiu-Sterman function. By extrapolating the derived spin dependent cross section to the small transverse momentum region, consistency between the collinear twist-3 approach and the transverse momentum dependent factorization approach is confirmed. As a byproduct of this work, we identify a term contributing to the scale evolution of tri-gluon correlations in the flavor-singlet case which was previously missed.

preprint2012arXiv

Exploring quark transverse momentum distributions with lattice QCD

We discuss in detail a method to study transverse momentum dependent parton distribution functions (TMDs) using lattice QCD. To develop the formalism and to obtain first numerical results, we directly implement a bi-local quark-quark operator connected by a straight Wilson line, allowing us to study T-even, "process-independent" TMDs. Beyond results for x-integrated TMDs and quark densities, we present a study of correlations in x and transverse momentum. Our calculations are based on domain wall valence quark propagators by the LHP collaboration calculated on top of gauge configurations provided by MILC with 2+1 flavors of asqtad-improved staggered sea quarks.

preprint2012arXiv

Higgs boson production in high energy proton-nucleus collisions

We study Higgs boson production from gluon-gluon fusion at mid-rapidity in high energy proton-nucleus collisions. For this process the presently still little known gluon distribution function $h_1^{\perp g}$ gives a numerically relevant contribution. We show by explicite calculation that using CGC (color glass condensate) model input the result obtained in the naive k_t factorization approach matches the result obtained in the TMD factorization framework for a dilute medium. We also verify the earlier finding that the k_t factorization formalism for Higgs production breaks down in a dense medium. In doing so we formulate a hybrid model which allows one to treat such reactions theoretically.

preprint2012arXiv

The first PDF moments for three dynamical flavors in Baryon Chiral Perturbation Theory

We present a calculation of generalized baryon form factors in the framework of three-flavor covariant baryon chiral perturbation theory at leading one-loop order. This is needed for lattice calculations of the first moments of generalized parton distribution functions. The formulae we derive can be used to guide the chiral extrapolation of such lattice results.

preprint2012arXiv

Thermodynamics of SU(N) Yang-Mills theories in 2+1 dimensions II - The deconfined phase

We present a non-perturbative study of the equation of state in the deconfined phase of Yang-Mills theories in D=2+1 dimensions. We introduce a holographic model, based on the improved holographic QCD model, from which we derive a non-trivial relation between the order of the deconfinement phase transition and the behavior of the trace of the energy-momentum tensor as a function of the temperature T. We compare the theoretical predictions of this holographic model with a new set of high-precision numerical results from lattice simulations of SU(N) theories with N=2, 3, 4, 5 and 6 colors. The latter reveal that, similarly to the D=3+1 case, the bulk equilibrium thermodynamic quantities (pressure, trace of the energy-momentum tensor, energy density and entropy density) exhibit nearly perfect proportionality to the number of gluons, and can be successfully compared with the holographic predictions in a broad range of temperatures. Finally, we also show that, again similarly to the D=3+1 case, the trace of the energy-momentum tensor appears to be proportional to T^2 in a wide temperature range, starting from approximately 1.2 T_c, where T_c denotes the critical deconfinement temperature.

preprint2012arXiv

Topology of dynamical lattice configurations including results from dynamical overlap fermions

We investigate how the topological charge density in lattice QCD simulations is affected by violations of chiral symmetry caused by the fermion action. To this end we compare lattice configurations generated with a number of different actions including first configurations generated with exact dynamical overlap quarks. We visualize the topological profiles after mild smearing. In the topological charge correlator we measure the size of the positive core, which is known to shrink to zero extension in the continuum limit. To leading order we find the core size to scale linearly with the lattice spacing with the same coefficient for all actions, even including quenched simulations. In the subleading term the different actions vary over a range of about 10 %. Our findings suggest that non-chiral lattice actions at current lattice spacings do not differ much for observables related to topology, both among themselves and compared to overlap fermions.

preprint2011arXiv

Azimuthal asymmetry in SIDIS off nuclei as probe for $\hat q$

The energy loss parameter $\hat q$ is one of the fundamental transport parameters of hadronic matter. Using the twist-4 collinear approach, we show that the $\cos ϕ$ azimuthal asymmetry in unpolarized semi-inclusive deeply inelastic scattering (SIDIS) off a large nucleus at intermediate transverse momentum is a sensitive observable for its determination. The effect is due to the suppression of the azimuthal asymmetry by final-state multiple scattering.

preprint2011arXiv

Chiral extrapolation of nucleon wave function normalization constants

Within the framework of two-flavor covariant baryon chiral perturbation theory we have expressed the Chernyak-Zhitnitsky, Ioffe and Dosch currents in terms of chiral fields to provide leading one-loop extrapolation formulae for the leading and next-to-leading twist normalization constants $f_N$, $λ_1$ and $λ_2$. Finite volume effects due to pion loops have been taken into account. The occurring low energy constants are fitted to data obtained from recent lattice QCD simulations in order to extract the values at the physical point.

preprint2011arXiv

Entropy Creation in Relativistic Heavy Ion Collisions

We review current ideas on entropy production during the different stages of a relativistic nuclear collision. This includes recent results on decoherence entropy and the entropy produced during the hydrodynamic phase by viscous effects. We start by a discussion of decoherence caused by gluon bremsstrahlung in the very first interactions of gluons from the colliding nuclei. We then present a general framework, based on the Husimi distribution function, for the calculation of entropy growth in quantum field theories, which is applicable to the early ("glasma") phase of the collision during which most of the entropy is generated. The entropy calculated from the Husimi distribution exhibits linear growth when the quantum field contains unstable modes and the growth rateis asymptotically equal to the Kolmogorov-Sinaï (KS) entropy. We outline how the approach can be used to investigate the problem of entropy production in a relativistic heavy-ion reaction from first principles. Finally we discuss some recent results on entropy production in the strong coupling limit, as obtained from AdS/CFT duality.

preprint2011arXiv

Holographic Thermalization

Using the AdS/CFT correspondence, we probe the scale-dependence of thermalization in strongly coupled field theories following a quench, via calculations of two-point functions, Wilson loops and entanglement entropy in d=2,3,4. In the saddlepoint approximation these probes are computed in AdS space in terms of invariant geometric objects - geodesics, minimal surfaces and minimal volumes. Our calculations for two-dimensional field theories are analytical. In our strongly coupled setting, all probes in all dimensions share certain universal features in their thermalization: (1) a slight delay in the onset of thermalization, (2) an apparent non-analyticity at the endpoint of thermalization, (3) top-down thermalization where the UV thermalizes first. For homogeneous initial conditions the entanglement entropy thermalizes slowest, and sets a timescale for equilibration that saturates a causality bound over the range of scales studied. The growth rate of entanglement entropy density is nearly volume-independent for small volumes, but slows for larger volumes.

preprint2011arXiv

Some results on excited hadrons in 2-flavor QCD

Results of hadron spectroscopy with two dynamical mass-degenerate chirally improved quarks are presented. Three ensembles with pion masses of 322(5), 470(4) and 525(7) MeV, lattices of size 16^3 \times 32, and lattice spacings close to 0.15 fm are investigated. We discuss the possible appearance of scattering states, considering masses and eigenvectors. Partially quenched results in the scalar channel suggest the presence of a 2-particle state, however, in most channels we cannot identify them. Where available, we compare to results from quenched simulations using the same action.

preprint2011arXiv

Threshold-Resummed Cross Section for the Drell-Yan Process in Pion-Nucleon Collisions at COMPASS

We present a study of the Drell-Yan process in pion-proton collisions including next-to-leading-logarithmic threshold-resummed contributions. We analyze rapidity-integrated as well as rapidity-differential cross sections in the kinematic regime relevant for the COMPASS fixed target experiment. We find that resummation leads to a significant enhancement of the cross section compared to fixed-order calculations in this regime. Particularly large corrections arise at large forward and backward rapidities of the lepton pair. We also study the scale dependence of the cross section and find it to be substantially reduced by threshold resummation.

preprint2011arXiv

Topology of dynamical lattice configurations including results from dynamical overlap fermions

We investigate how the topological charge density in lattice QCD simulations is affected by violations of chiral symmetry in different fermion actions. To this end we compare lattice configurations generated with a number of different actions including first configurations generated with exact overlap quarks. We visualize the topological profiles after mild smearing. In the topological charge correlator we measure the size of the positive core, which is known to vanish in the continuum limit. To leading order we find the core size to scale linearly with the lattice spacing with the same coefficient for all actions, even including quenched simulations. In the subleading term the different actions vary over a range of about 10 %. Our findings suggest that non-chiral lattice actions at current lattice spacings do not differ much for this specific observable related to topology, both among themselves and compared to overlap fermions.

preprint2010arXiv

A lattice Dirac operator for QCD with light dynamical quarks

In QCD chiral symmetry is explicitly broken by quark masses, the effect of which can be described reliably by chiral perturbation theory. Effects of explicit chiral symmetry breaking by the lattice regularisation of the Dirac operator, typically parametrised by the residual mass, should be negligible for almost all observables if the residual mass of the Dirac operator is much smaller than the quark mass. However, maintaining a small residual mass becomes increasingly expensive as the quark mass decreases towards the physical value and the continuum limit is approached. We investigate the feasibility of using a new approximately chiral Dirac operator with a small residual mass as an alternative to overlap and domain wall fermions for lattice simulations. Our Dirac operator is constructed from a Zolotarev rational approximation for the matrix sign function that is optimal for bulk modes of the Hermitian kernel Dirac operator but not for the low-lying parts of its spectrum. We test our operator on various 32^3\times64 lattices, comparing the residual mass and the performance of the Hybrid Monte Carlo algorithm at a similar lattice spacing and pion mass with a hyperbolic tangent operator as used by domain wall fermions. We find that our approximations have a significantly smaller residual mass than domain wall fermions at a similar computational cost, and still admit topological charge change.

preprint2010arXiv

Charge Fluctuations from the Chiral Magnetic Effect in Nuclear Collisions

We derive a nonlocal effective Lagrangian for the chiral magnetic effect. An electric field is generated by winding number fluctuations of the nonabelian gauge field in the presence of a strong magnetic field. We estimate the magnitude of charge asymmetry fluctuations with respect to the reaction plane induced by the chiral magnetic effect in relativistic heavy ion collisions to be less than $10^{-6}$, several orders of magnitude smaller than the signal observed in the STAR experiment.

preprint2010arXiv

Euler-Heisenberg Lagrangian to all orders in the magnetic field and the Chiral Magnetic Effect

In high energy heavy ion collisions as well as in astrophysical objects like magnetars extreme magnetic field strengths are reached. Thus, there exists a need to calculate divers QED processes to all orders in the magnetic field. We calculate the vacuum polarization graph in second order of the electric field and all orders of the magnetic field resulting in a generalization of the Euler-Heisenberg Lagrangian. We perform the calculation in the effective Lagrangian approach of J. Schwinger as well as using modified Feynman rules. We find that both approaches give the same results provided that the different finite renormalization terms are taken into account. Our results imply that any quantitative explanation of the recently proposed Chiral Magnetic Effect has to take 'Strong QED' effects into account, because these corrections are huge.

preprint2010arXiv

Filtered topological structure of the QCD vacuum: Effects of dynamical quarks

We systematically compare filtering methods used to extract topological structures on SU(3) lattice configurations. We show that there is a strong correlation of the topological charge densities obtained by APE and Stout smearing. To get rid of artifacts of these methods, we analyse structures that are also seen by Laplace filtering and indeed identify artifacts for strong smearing. The topological charge density in this combined analysis is more fragmented in the presence of dynamical quarks. A power law exponent that characterises the distribution of filtered topological clusters turns out to be not far off the values of an instanton gas model.

preprint2010arXiv

Meson and baryon spectrum for QCD with two light dynamical quarks

We present results of meson and baryon spectroscopy using the Chirally Improved Dirac operator on lattices of size 16**3 x 32 with two mass-degenerate light sea quarks. Three ensembles with pion masses of 322(5), 470(4) and 525(7) MeV and lattice spacings close to 0.15 fm are investigated. Results on ground and excited states for several channels are given, including spin two mesons and hadrons with strange valence quarks. The analysis of the states is done with the variational method, including two kinds of Gaussian sources and derivative sources. We obtain several ground states fairly precisely and find radial excitations in various channels. Excited baryon results seem to suffer from finite size effects, in particular at small pion masses. We discuss the possible appearance of scattering states in various channels, considering masses and eigenvectors. Partially quenched results in the scalar channel suggest the presence of a 2-particle state, however, in most channels we cannot identify them. Where available, we compare our results to results of quenched simulations using the same action.

preprint2010arXiv

Soft-Gluon Resummation and the Valence Parton Distribution Function of the Pion

We determine the valence parton distribution function of the pion by performing a new analysis of data for the Drell-Yan process $π^- N\to μ^+μ^-X$. Compared to previous analyses, we include next-to-leading-logarithmic threshold resummation effects in the calculation of the Drell-Yan cross section. As a result of these, we find a considerably softer valence distribution at high momentum fractions $x$ than obtained in previous next-to-leading order analyses, in line with expectations based on perturbative-QCD counting rules or Dyson-Schwinger equations.

preprint2009arXiv

Comparing the vacuum structure of quenched and dynamical configurations

We systematically compare filtering methods used to extract topological structures on SU(3) lattice configurations. We show that there is a strong correlation of the topological charge densities obtained by APE and Stout smearing. To get rid of artifacts of these methods, we analyze structures that are also seen by Laplace filtering. This combined analysis shows that the topological charge density is more fragmented in the presence of dynamical quarks.

preprint2008arXiv

Heavy-light hadrons and their excitations

We study the excitations of hadrons containing a single heavy quark. We present meson and baryon mass splittings and ratios of meson decay constants resulting from quenched and dynamical two-flavor configurations. Light quarks are simulated using the Chirally Improved (CI) lattice Dirac operator. The heavy quark is approximated by a static propagator, appropriate for the $b$ quark on our lattices ($1/a \sim 1-2$ GeV). We also include some preliminary calculations of the heavy-quark kinetic corrections to the states.

preprint2008arXiv

Wave functions of the nucleon and its parity partner from lattice QCD

We compute moments of distribution amplitudes using gauge configurations with two flavors of clover fermions from QCDSF/DIK and operators which are optimized with respect to their behavior under the lattice symmetries. The knowledge of these quantities helps in understanding the internal structure of hadrons and in the analysis of (semi-)exclusive processes. We present results for the nucleon distribution amplitude which suggest that the asymmetries (the deviations from the asymptotic form) are smaller than indicated by sum rule calculations. Using the same approach we were also able to calculate the same quantities for the $N^\star(1535)$, the parity partner of the nucleon. These results show a stronger deviation from the asymptotic form.

preprint2004arXiv

Excited hadrons from improved interpolating fields

The calculation of quark propagators for Ginsparg-Wilson-type Dirac operators is costly and thus limited to a few different sources. We present a new approach for determining spatially optimized operators for lattice spectroscopy of excited hadrons. Jacobi smeared quark sources with different widths are combined to construct hadron operators with different spatial wave functions. We study the Roper state and excited rho and pion mesons.