Linearly polarized gluon density in the rescaling model

The low-x behavior of the linearly polarized gluon density h_{g}(x,k_t^2, Q^2) in nuclei is studied in the rescaling model at small transverse momentum k_t.


Introduction
Determination of parton (quark and gluon) distribution functions (PDFs) in a proton and nuclei is a rather important task for modern high energy physics.In particular, detailed knowledge on the gluon densities is necessary for experiments planned at the Large Hadron Collider (LHC) and future colliders, such as Electron-Ion Collider (EIC), Future Circular hadronelectron Collider (FCC-he), Electron-Ion Collider in China (EicC) and Nuclotron-based Ion Collider fAcility (NICA) [1,2,3,4,5].For unpolarized cases, there are a distribution of unpolarized gluons, denoted as f g (x, Q2 ), and a distribution of linearly polarized gluons h g (x, k 2 t , Q 2 ), which corresponds to interference between ±1 gluon helicity states 1 .Compared to f g (x, Q 2 ), function h g (x, k 2 t , Q 2 ) is currently poorly known 2 in comparison with f g (x, Q 2 ) and depends on the gluon transverse momentum k t (so called Transverse Momentum Dependent, or TMD gluon density).A theoretical upper bound for h g (x, k 2 t , Q 2 ) was obtained [7,8].Previously, we have derived an analytical expressions for linearly polarized gluon density in a proton and investigated its behavior at low x [9].Our analysis was based on the small-x asymptotics for sea quark and gluon densities calculated in the generalized double asymptotic scale (DAS) approach [10,11] (see also [12]) and was done with leading order (LO) accuracy.
In the present note we extend the consideration [9] for nuclei.
As it is known, the study of deep inelastic scattering (DIS) of leptons on nuclei reveals the appearance of a significant nuclear effect, which excludes the naive idea of the nucleus as a system of quasi-free nucleons (see, for example, [13,14] for review).This effect was first discovered by the European Muon Collaboration in the domain of valence quark dominance, so it was called the EMC effect [15].Currently, there are two main approaches to study the nuclear PDFs (nPDFs).In the first, which is currently more common, nPDFs are extracted from global fits (see the recent review [16] and references therein) to nuclear data using empirical parametrization of their normalizations and numerical solution of Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations [17] to describe the corresponding QCD evolution (dependence on scale Q 2 ).The second strategy is based on some nPDF models (see [18] for more information).Here we follow the rescaling model [19] based on the assumption [20] that the effective size of gluon and quark confinement in the nucleus is greater than in the free nucleon.Within the perturbative QCD, it was pointed out [19,20] that this confinement rescaling predicts that nPDFs and PDFs can be connected by simply scaling the argument Q 2 (see also a review [21]).Thus, one can say that the rescaling model demonstrates the features inherent in both approaches: there are certain relationships between conventional and nuclear PDFs that arise as a result of shifting the values of the kinematic variable Q 2 and, at the same time, both densities obey DGLAP equations.
Initially, the rescaling model was proposed for the domain of valence quarks dominance, 0.2 ≤ x ≤ 0.8, where x is the Bjorken variable.Recently it was extended to a small x [22,23], where certain shadowing and antishadowing effects3 were found for the sea quark and gluon densities.Our main goal is to apply the rescaling model to linearly polarized gluon density h g (x, k 2 t , Q 2 ) and show its nuclear modification for small x values.

Approach
At LO of perturbation theory the linearly polarized gluon density in a proton h g (x, k 2 t , Q 2 ) at low k t values has the following form [28] (all PDFs are multiplied by x): where is related with the conventional strong coupling constant α s (Q 2 ) and β 0 = 11 − 2f /3 is the first coefficient of QCD β-function in the MS-scheme and f is the number of active quarks.Considering low x asymptotics for gluon density f g (x, Q 2 ), we have [11]: where A g and A q are magnitudes of gluon and (sea) quark densities at some initial scale Q 2 0 , I 0 is modified Bessel function.Here also are the singular and regular parts of the anomalous dimensions.
The linearly polarized gluon density in a proton h g (x, Q 2 ), which is important mostly at low x reqion, has the following form [9] where I i is modified Bessel function.

Rescaling model
In the rescaling model [19], the valence part of quark densities are changed in the case of a A nucleus at intermediate and large values of variable x (0.2 ≤ x ≤ 0.8) as follows where the new scale i.e. the kernel modification of the main variable s depends on the Q 2 -independent parameter δ A V having small values (see Tables 2 and 3 in [22]).

Rescaling model at low x
In [22], the PDF asymptotics shown in (3) were applied to the small x region of the EMC effect, using the simple fact that the sea quark and gluons densities increase with increasing Q 2 .Thus, in the case of nuclei, the PDF evolution scale is less than Q 2 , and this can directly reproduce the shadowing effect observed with global fits.Since there are two components for each gluon density, see ( 3) and ( 6), we have two free parameters Q 2 A,± , which can be determined from the analysis of experimental data for the EMC effect at low x values.
Note that it is usually convenient to study the following ratio x 7 Li, AKL 7 Li, ACKL 16 O, AKL 16 O, ACKL 16 O, EPPS21 16 O, nNNPDF3.0 The blue and yellow solid lines are for R AD h (x, Q 2 ) obtained in the present paper for 7 Li and 16 O, respectively.The purple dotted and red dashed lines are for R AD g (x, Q 2 ) obtained in [30].The blue dash-dotted and green lines and bands are borrowed from [16].
for both conventional f g (x, Q 2 ) and linearly polarized h g (x, Q 2 ) gluon densities (here t = f, h).Taking advantage of the fact that the nuclear effect in the deuteron is very small4 : t D g (x, Q 2 ) ≈ t g (x, Q 2 ), we can assume that The expressions for f ± g (x, Q 2 ) and h ± g (x, Q 2 ) are given in ( 3) and ( 6), respectively, and the corresponding values of s AD ± turned out to be where the results for δ AD ± can be found [22].

Results
Our numerical results obtained for R AD g (x, Q 2 ) and R AD h (x, Q 2 ) are shown in Fig. 1.As it was discussed [30], our predictions for R AD g (x, Q 2 ) are very close to the ones [22] obtained at x ≤ 10 −2 , since the parameters δ AD ± are taken from this article.Moreover, in the low x our calculations are also close to the recent results [16] obtained by fitting experimental data.We see that predicted R AD h (x, Q 2 ) are quite similar to the R AD g (x, Q 2 ), but they are less affected by nuclear effects.Indeed, in both cases: in the shadowing area (x ≤ 0.05) and in the antishadowing area (0.05 ≤ x ≤ 0.1), the values of R AD h (x, Q 2 ) less differ from 1 than the corresponding values of R AD g (x, Q 2 ).So that, the derived expressions could be useful for subsequent phenomenological applications.
In our forthcoming studies, we plan to extend our analysis for x ≥ 0.1 using recently obtained PDF set [32].It is a combination of exact analytical solutions of DGLAP equations for small and large values of x and, therefore, could be applicable over the entire range of x.We also plan to apply other modification models for nPDFs [13,14,16,18,33].Studies in a wide x range could be supplemented with IMParton framework [34].In addition, we will study nuclear modifications of transverse momentum dependent PDFs (or TMD gluon and quark densities) [32,35], which are close relatives of h g (x, k 2 t , Q 2 ).Now the these quantities become very popular (see [36] and references and discussions therein).
Researches described in Section 2 were supported by the Russian Science Foundation under grant 22-22-00387.Studies described and performed in Sections 3 and 4 were supported by the Russian Science Foundation under grant 22-22-00119.