High-energy resummation in $\Lambda_c$ baryon production

We present a study on inclusive emissions of a double $\Lambda_c$ or of a $\Lambda_c$ plus a light-flavored jet system as probe channels in the semi-hard regime of QCD. Our formalism relies on the so-called hybrid high-energy/collinear factorization, where the standard collinear description is supplemented by the $t$-channel resummation \`a la BFKL of energy logarithms up to the next-to-leading accuracy. We make use of the JETHAD modular interface, suited to the analysis of different semi-hard reactions, employing the novel KKSS19 parameterization for the description of parton fragmentation into $\Lambda_c$ baryons. We provide predictions for rapidity distributions and azimuthal correlations, that can be studied at current and forthcoming LHC configurations, hunting for possible stabilizing effects of the high-energy series.


Hors d'oeuvre
Fixed-order perturbative QCD is a fundamental and successful approach for the description of strong-interaction processes ruled by one or more hard kinematic scales, i.e. scales much larger than Λ QCD . In some peculiar kinematic regimes, however, fixed-order perturbative calculations are no longer adequate to build sensible predictions and all-order resummations of some definite classes of contributions to the perturbative series become mandatory. One such kinematic regime is the semi-hard one [1], characterized by a center-of-mass energy, √ s, much larger than the hard scale(s) of the process and, therefore, of the QCD mass scale. In this regime, large energy logarithms enter the perturbative series with a power which increases with the perturbative order, thus systematically compensating the smallness of the coupling and calling for all-order resummation. The approximation where only leading-order energy logarithms are resummed is called LLA; if also next-to-leading logarithms are resummed, we have the next-to-leading approximation or NLA. The LHC, with its world-record beam energies, is the best possible stage for semi-hard processes and a unique opportunity for testing theoretical approaches based on energy resummations in pertubative QCD.
A procedure for the systematic inclusion of large energy logarithms to all orders in perturbative QCD, both in the LLA and in the NLA, is offered by the Balitsky-Fadin-Kuraev-Lipatov (BFKL) [2][3][4][5]. In this framework, the cross sections of hadronic processes take a peculiar factorized form, given by the convolution of two impact factors, related to the fragmentation of each colliding particle to an identified final-state object, and a process-independent Green's function. The evolution of the BFKL Green's function is built out of an integral equation, whose kernel is known at the next-to-leading order (NLO) both for forward scattering (i.e. for t = 0 and color singlet in the t-channel) [6,7] and for any fixed, not growing with s, momentum transfer t and any possible two-gluon color state in the t-channel [8][9][10][11][12].
In this paper we propose the inclusive production of Λ c baryons in semi-hard regimes as a further probe of the high-energy resummation in kinematic ranges typical of current and upcoming experimental studies at the LHC. In particular, we focus on final states featuring the emission of a forward Λ c particle accompanied by a backward Λ c (panel a) of Fig. 1) or by a jet (panel b) of Fig. 1). The final-state inclusiveness is warranted by the emission of undetected hard partons strongly ordered in rapidity, while energy scales at play are large enough to make a variable-flavor number-scheme (VFNS) description valid (see, e.g., Refs. [96,97]). Λ c production at the LHC in central-rapidity ALICE and in forward-rapidity LHCb regimes was considered in Ref. [98] within k T -factorization.
Similar final-state configurations have been considered with NLA accuracy in [59], with Λ hyperons playing the role of the Λ c . The motivations brought in Ref. [59] in support of the Λ hyperon case plainly extend to the Λ c . What makes the latter case interesting on its own is the following consideration. It is well known that NLA BFKL suffers from severe instabilities, due to being the NLO corrections in the BFKL kernel and in impact factors large in size and opposite in sign with respect to the LO. In Ref. [59] and in our previous NLA BFKL studies the procedure adopted to tame this instability was the use of the MOM scheme for the strong coupling renormalization with Brodsky-Lepage-Mackenzie (BLM) optimization [99][100][101][102] of the renormalization scale fixing. This allowed to obtain stable results for cross sections and azimuthal correlations or ratios between them, at the price of fixing the renormalization scale at values much larger than the kinematic ones, typical of the considered process. The interesting feature of the processes suggested in this paper is that the BLM optimization scales turn out to be much smaller than in previous studies, due to a subtle interplay between proton PDFs and Λ c fragmentation functions (FFs), discussed in details below. The obvious interpretation of this phenomenon is that Λ c production, due to the particular shape of the corresponding FFs, leads to an improved stability of the BFKL series, making the suggested processes worthy of attention in LHC experimental analyses. Figure 1: Diagrammatic representation of a) the double Λ c and of b) the Λ c + jet production in hybrid high-energy/collinear factorization. Red blobs stand for proton collinear PDFs, whereas bordeaux rectangles denote baryon collinear FFs and green arrows refer to the jet selection function. The BFKL gluon Green's function, represented by the yellow central blob, is connected to impact factors through Reggeon (zigzag) lines. Diagrams were realized by making use of the JaxoDraw 2.0 interface [103].
The paper is organized as follows: in Section 2 we introduce our theoretical setup, in Section 3 we show and discuss results of our numerical analysis, in Section 4 we draw our conclusions and propose future prospects.
2 Hybrid factorization for Λ c production 2.1 NLA-resummed cross section The two hadronic reactions object of our investigation are (see Fig. 1): where a Λ c particle (we are inclusive on the baryon charge) is accompanied by another Λ c or by a light-flavored jet 1 . Produced objects have transverse momenta, | p 1,2 | Λ QCD , large enough 1 Since we are working in the VFNS, where all parton quarks are treated as light ones, the jet can be generated by any quark, up to the active flavor number, n f , or by a gluon.
to justify the use of pQCD. At the same time, their large distance in rapidity, ∆Y ≡ y 1 − y 2 , ensures diffractive final states. An undetected gluon activity, denoted as X, accompanies manyparticle emissions. We will consider values of ∆Y larger than zero, so that the 1-labeled Λ c will be always a forward particle, while the 2-labeled object will be always a backward one. We take as Sudakov basis the one generated by momenta of incoming protons, P a,b , thus having P 2 a,b = 0 and (P a ·P b ) = s/2, with s the hadronic center-of-mass energy squared. In this way, one can decompose the final-state transverse momenta on that basis: where the subscript O 1 refers to the first emitted object (always a baryon) and O 2 refers to the second one (another baryon or a jet). In the center-of-mass frame, the following relations between longitudinal fractions, rapidities, and transverse momenta of emitted particles hold A pure fixed-order treatment relies on the standard collinear factorization, where the onshell partonic cross section is convoluted with PDFs and FFs. In the double Λ c channel (panel a) of Fig. 1 For the sake of readability, the explicit dependence of PDFs, FFs, and the partonic cross section on the factorization scale has been dropped everywhere. At variance with the standard collinear approach, we build our hybrid setup by starting from the high-energy factorization which naturally emerges inside the BFKL formalism, and then we improve our description by plugging collinear ingredients in. It is convenient to write the cross section as a Fourier series of the azimuthal-angle coefficients, C n≥0 , this leading to the following general expression where ϕ 1,2 are the azimuthal angles of the tagged objects and ϕ ≡ ϕ 1 − ϕ 2 − π. The BFKL approach provides us with a consistent definition of NLA-resummed azimuthal coefficients, C n = C NLA n , whose definition in the MS renormalization scheme reads (details on the derivation can be found in Ref. [39]) with γ E the Euler-Mascheroni constant and ψ(z) ≡ Γ (z)/Γ(z) the logarithmic derivative of the Gamma function, whereasχ(n, ν) comes from the NLO correction to the BFKL kernel and was calculated in Ref. [104] (see also Ref. [105]). Then, the c 1,2 functions stand for the LO impact factors, calculated in the Mellin space, describing the emissions of forward/backward objects. The LO impact factor for the production of a Λ c baryon reads Similarly, the light jet is portrayed by the corresponding LO impact factor Furthermore, the f (ν) function is defined as The remaining quantities in Eq. (8),ĉ 1,2 , are the NLO corrections to the hadron impact factor, whose analytic expression was calculated in Ref. [22], and the jet impact-factor one. For this one we will employ a simple version (see Ref. [20]), suited to numerical studies, which encodes a jet selection function calculated in the "small-cone" approximation (SCA) [106,107]. Eqs. (8), (10), and (11) clearly show how our hybrid factorization is realized. The cross section is high-energy factorized in terms of gluon Green's function and impact factors, and the latter ones embody collinear PDFs and FFs. We remark that the description of Λ c particles in terms of light-hadron impact factors (Eq. (10)) is adequate, provided that energy scales are much larger than the Λ c mass. This condition is guaranteed by the transverse-momentum ranges of our interest (see Section 3.1).
By truncating the NLA azimuthal coefficients in Eq. (8) to the O(α 3 s ), we obtain a fixedorder formula that acts an effective high-energy DGLAP counterpart of our BKFL-resummed expression. Thus, we keep the leading-power asymptotic pattern of a pure NLO DGLAP description, eliminating at the same time those terms which are suppressed by inverse powers ofŝ. Our DGLAP formula can be cast in the form where an expansion up to terms proportional to α s (µ R ) replaces the BFKL exponentiated kernel. Starting from Eqs. (8) and (13), it is possible to obtain corresponding expressions in the MOM scheme by performing the finite renormalization where 2.3439 and ξ a gauge parameter, fixed at zero in the following.

BLM scale optimization
According to the BLM method, the optimal renormalization-scale value, labeled as µ BLM R , is the value of µ R that cancels the non-conformal, β 0 -dependent part of a given observable. A suitable procedure, recently set up [46], allows us to remove all non-conformal contributions that appear both in the NLA BFKL kernel and in the NLO non-universal impact factors of high-energy distributions. Its application makes µ BLM R dependent on the energy of the process and thus on ∆Y .
As a preliminary step, we need to perform a finite renormalization from the MS scheme to the MOM one (see Eq. (14)). Thus, the condition for the BLM scale setting for a given azimuthal coefficient, C n , is given as a solution of the integral equation with dΦ(y 1,2 , | p 1,2 |, ∆Y ) the final-state differential phase space (see Section 3.1), and For the sake of convenience, we introduce the ratio between the BLM scale and the natural one suggested by the kinematic of the process, namely µ N ≡ √ m 1⊥ m 2⊥ , so that C BLM µ ≡ µ BLM R /µ N , and we look for values of C µ R that solve Eq. (16). Here, m i⊥ stands for the i-th particle transverse mass. Therefore, one always has m 1⊥ = | p 1 | 2 + m 2 Λc , with m Λc = 2.286 GeV. Then, m 2⊥ = | p 2 | 2 + m 2 Λc in the double Λ c production, while m 2⊥ coincides with the jet p T in the other channel. We set µ F = µ R everywhere, as assumed by most of the existent PDF parameterizations.
Finally, the BLM scale value is inserted into expressions for the integrated coefficients, and we get the following NLA BFKL formula in the MOM renormalization scheme wherec 1,2 (n, ν, | p 1,2 |, x 1,2 ) stand for the NLO impact-factor corrections after removing the β 0dependent terms, which can be universally expressed via the LO impact factors, c 1,2 . One hasc Analogously, we expand and truncate to the O(α 3 s ) the BFKL kernel in Eq. (19), this getting a BLM-MOM expression for the high-energy DGLAP limit (21) The corresponding formulae of Eqs. (19) and (21) in the MS scheme are obtained by making the substitutions (note that the value of the renormalization scale is left unchanged) 3 Numerical analysis

Rapidity and azimuthal distributions
Key ingredients for building our distributions are the azimuthal coefficients integrated over rapidity and transverse momenta of the two tagged object, their rapidity separation being kept fixed: The C n and C n coefficients can refer to the corresponding NLA BFKL calculations (see Eq. (8)) or the ones taken in the high-energy DGLAP limit (see Eq. (13)).
In order to match realistic LHC configurations, we allow the rapidity of Λ c baryons to be in the range from −2.0 and 2.0, while their transverse momentum goes from 10 GeV to p max H 21.5 GeV. These cuts are borrowed from typical analyses of the Λ b particle at CMS [108], whereas the p max H value is constrained by the energy-scale lower cutoff on the FF sets (see Section 3.2). As for the jet, we consider standard CMS configurations [109], namely |y J | < 4.7 and 35 GeV < p J < 60 GeV. A major benefit of allowing for jet detection also by the CASTOR ultra-backward detector (−6.6 < y J < −5.2) [110,111] is the possibility to test our observables on larger values of rapidity intervals, say ∆Y 9. However, in our previous work on Λ hyperons [59] we have highlighted how, in this kinematics regime, large values of partons' longidudinal fractions effectively restrict the weight of the undetected gluon radiation. This leads to the appearance of large Sudakov-type double logarithms (threshold double logarithms) in the perturbative series, that have to be resummed to all orders. Since this resummation mechanism has not been yet embodied in our approach, we postpone the investigation of Λ c emissions in CASTOR-jet configurations to future, dedicated studies.
The integrated coefficients defined in Eq. (23) permit us to study the ϕ-summed cross section, C 0 , and the azimuthal-correlation ratios, R nm = C m /C m , as functions of ∆Y . The R n0 ratios have an immediate physical interpretation, being the moments cos nϕ , while the ones without zero indices correspond to ratios of cosines, cos nϕ / cos mϕ [112,113]. In our study we fix the center-of-mass energy at √ s = 13 TeV.

JETHAD settings
We performed our phenomenological studies by making use of the JETHAD modular interface under development at our Group. Numerical computations of distributions as well as BLM scales were done via the Fortran 2008 work package implemented in JETHAD, whereas an automatized Python 3.0 analyzer was developed for elaboration of results. Collinear PDFs were calculated via the MMHT14 NLO PDF set [114] as provided by the LHAPDFv6.2.1 interpolator [115], and a two-loop running coupling with α s (M Z ) = 0.11707 and a dynamic-flavor threshold was chosen. We described the parton fragmentation into Λ c baryons in terms of the novel KKSS19 NLO FF set [116] (see also Refs. [117,118]), whose native implementation was directly linked to JETHAD. This parameterization mainly relies on a descriptionà la Bowler [119] for c and b quark/antiquark flavors. Technical details on the fitting procedure are presented in Section IV of Ref. [116]. Here, we limit ourselves to saying that the use of a given VFNS PDF or FF set is valid in our approach, provided that energy scales at work are much larger than thresholds for DGLAP evolution of heavier quarks. Since our scales are always higher than 10 GeV (see Section 3.1), while KKSS19 thresholds for c and b quarks are respectively 1.5 GeV and 5 GeV, this requirement is fulfilled. Lighter-hadron emissions (Λ hyperons, pions, kaons, and protons) in Figs. 2 and 3 were described in terms AKK08 NLO FFs [120], which are the closest in technology to KKSS19.
Nested integrations over phase space, ν variable, and longitudinal fractions inside impact factors were mainly evaluated via an adaptative-quadrature strategy provided by the JETHAD integrators. In all cases, their global uncertainty was kept below 1%.
We performed a dedicated study on the sensitivity of our predictions on scale variation. More in particular, we gauged the effect of concurrently varying µ R and µ F around their natural values and their BLM optimal ones, in the range 1/2 to two. The C µ parameter entering the figures in Section 3.3 stands for the ratio C µ = µ R,F /µ N . Error bands in all our plots show the combined uncertainty coming from numerical integration and scale variation, this latter being the dominant one. All calculations were done in the MS scheme.
An extension of our analysis that includes all systematic uncertainties together with a comparison of results obtained with different FF sets, as the DMS20 one [121], is postponed to a future work.

Discussion
As a preliminary analysis, we compare results for the BLM-scale parameter, C BLM µ , as a function of the rapidity interval, ∆Y , for different species of emitted hadrons (Λ, π, K and p). Notably, BLM scales for Λ c emissions are much lower (although still larger than natural ones) than the ones obtained when lighter-hadron species are detected. The effect is much more evident in the double Λ c channel (left panels of Fig. 2) with respect to the Λ plus jet one (right panels of Fig. 2), this making us speculate that the production Λ c baryons could act as a stabilizer of the BFKL series under higher-order corrections. Indeed, since applying the BLM method effectively translates into a rise of energy-scale values that reduces the weight of next-to-leading contributions, lower values of C BLM µ indicate that the high-energy series was already (partially) stable, before applying BLM. A straightforward way to corroborate this statement consists in comparing predictions for observables of our interest both at natural scales and at BLM ones.
Upper panels of Fig. 3 show the ∆Y -dependence of the ϕ-summed cross section, C 0 , in the double Λ c channel, together with corresponding predictions for detection of Λ hyperons. Here, two competing effects come into play. On the one hand, partonic cross sections grow with energy, as predicted by BFKL. On the other hand, collinear densities quench predictions when ∆Y increases. The net result is a downtrend with ∆Y of C 0 distributions. We note that NLA bands are almost nested (except for very large values of ∆Y ) inside LLA ones and they are generally narrower in the Λ c case. This is a clear effect of a (partially) reached stability of the high-energy series, for both hadron emissions. However, while predictions for hyperons lose almost one order of magnitude when passing from natural scales to the expanded BLM ones (from left to right panel), results for Λ c baryons are much more stable, the NLA band becoming even wider in the BLM case. The stability is partially lost when Λ c particle is accompanied by a jet, as shown in lower panels of Fig. 3. Here, LLA and NLA bands are almost disjoined at natural scales (left panel), while in the BLM case (right panel) they come closer to each other for hyperon plus jet, and almost entirely contained for Λ c plus jet. This study on C 0 clearly highlights how Λ c emissions allow for a stabilization of the resummed series, that cannot be obtained with lighter hadrons 2 , nor with the associated production of a light-flavored jet. The detection of Λ c particles makes our distributions from one to two order of magnitude lower than the corresponding Λ hyperon cases. This helps to dampen, from the experimental point of view, minimum-bias contaminations, in a more effective way when Λ c are tagged. At the same time, statistics is favorable, since values of C 0 for Λ c emissions are always larger than 10 −2 nb. We checked that the different behavior of our predictions between Λ c baryons and other hadrons is not artificially generated by the different lower cutoff for the momentum fraction in the native FF grids, 10 −4 for KKSS19 versus 5 × 10 −2 for AKK08. According to our numerical tests, both the C BLM µ parameters and the C 0 distributions are left almost unchanged when the lower cutoff for KKSS19 is raised up to the AKK08 one. This feature was expected, since the major contribution to cross section is given by FF longitudinal fractions larger than 10 −1 (see discussion Section in Ref. [52]).
Predictions for R nm azimuthal ratios in the double Λ c channel at natural and at BLM scales are presented in Figs. 4 and 5, respectively. The downtrend of all these ratios when ∆Y grows is a well-know signal of the onset of high-energy dynamics. Larger rapidity distances heighten the weight of undetected gluons, thus leading to a decorrelation pattern in the azimuthal plane, which is more pronounced in pure LLA series. At variance with the C 0 case, here we observe a reversed situation, where NLA BLM-optimized results exhibit a narrower uncertainty bands and are closer to the corresponding LLA ones, with respect to what happens at natural scales. This dichotomy is much more emphasized in the Λ c plus jet channel. In particular, instabilities rising at natural scales are strong enough (although being milder than the ones observed in the Mueller-Navelet dijet channel) to prevent any realistic analysis. Therefore, we show only the behavior of R nm ratios after BLM optimization (Fig. 6), whose patterns are in line with corresponding predictions for the double Λ c production. The R 21 ratio exhibits a fair stability under NLA corrections for both scale choices and both final-state channels (see lower left panels of Figs. 4-6). In the presented plots, the value of the R 10 moment exceeds one in the small-∆Y region. This unphysical effect is fairly explained by the fact that contributions which are powersuppressed in energy and are not included in our BFKL treatment start to become relevant in those kinematic ranges, thus worsening the accuracy of our predictions.
We complete our analysis by comparing NLA predictions for Λ c plus jet final-state configurations with the corresponding ones calculated in our high-energy DGLAP limit. The ∆Y dependence of R 10 and R 20 azimuthal correlations with BLM scale optimization is presented in Fig. 7. As expected, the distance between BFKL and DGLAP, already marked at low ∆Y , becomes sharper and sharper when the rapidity interval grows. This phenomenon, already observed in the dijet [47] as well as in the hadron plus jet channel [59,122], is easily explained. At variance with the BFKL case, the limited number of inclusive gluon emissions due to the truncation of the perturbative series make the two final-state objects almost fully correlated (namely quasi back-to-back produced) independently of the value of ∆Y . The genuinely asymmetric kinematic configuration provided by the Λ c plus jet reaction suppresses the Born contribution, thus magnifying the distance between BFKL and DGLAP.
The overall outcome is that the inclusive detection of Λ c of baryons in semi-hard reactions allows for a stabilization of the high-energy resummation under higher-order corrections. A similar effect has been already observed in processes involving the production of massive bosons, such as the recently proposed Higgs plus jet channel [75]. However, while in that case the large energy scales provided by the Higgs transverse mass act as "natural" stabilizers for the BFKL series, here it comes as an intrinsic feature of the Λ c production. Moreover, at variance with the Higgs plus jet case, where the formal description relies on a partial NLA treatment, here the stabilizing effect is manifest in the full NLA BFKL.

More on Λ c fragmentation
In this Section we investigate the connection between Λ c FFs and the stability of cross sections under energy-scale variation. In upper panels of Fig. 8 we compare the µ F -dependence of KKSS19 Λ c FFs (left) with AKK08 Λ ones (right) for a value of the hadron momentum fraction typical of our analysis, namely z = 0.5. As expected, heavy flavors (c-and b-quarks) heavily dominate in Λ c fragmentation, whereas s-quark prevails in Λ emission, here the lighter-quark species and the gluon giving however a more appreciable contribution. Notably, KKSS19 FFs smoothly increase with µ F until they reach a constant value (apart from the b-quark, which decreases and then becomes constant). Conversely, AKK08 functions soften when µ F raises. This difference turns out to be relevant when FFs are convoluted with PDFs in our LO impact factors (Eq. (10)). Since the dominant contribution to PDFs in the kinematic sector of our interest, where the longitudinal-momentum fraction x ranges approximately from 10 −4 to 10 −2 , is given by the gluon (see lower panels of Fig. 8), the behavior of the gluon FF also becomes relevant. Indeed, the employment of large scales, such as the BLM ones, gives rise to two competing effects. On the one side, larger µ R values translate in a numerically smaller running coupling, both in the exponentiated kernel and the impact factors. On the other side, larger µ F values heighten the gluon-PDF contribution. When this last is convoluted with a smooth-behaved, non-decreasing gluon FF, such as the Λ c , the two features offset each other, thus generating the stability of cross sections under scale variation discussed in Section 3.3. Vice versa, the downtrend with µ F of the Λ-hyperon gluon FF spoils the balance between the two effects, thus preventing cross sections from reaching stability. This statement is supported by predictions shown in Fig. 9. Here, we test C 0 for double production of Λ c baryons (upper left panel) or Λ hyperons (upper right panel) without applying the BLM method, under a progressive variation of energy-scales in a wider range that includes the typical BLM ones, 1 < C µ < 30. We compare these results with corresponding ones for the inclusive emission, in the same kinematic domain, of two toy hadrons (lower panel) described in term of the following, flavor and µ F -independent model of FFs with N 1.5 × 10 −5 and λ = 0.2. The z-dependence of the model does not play a crucial role in our test. We clearly note that C 0 exhibits a fair stability under progressive scale variation both in the Λ c and in the toy-hadron channels, while its sensitivity spans over almost one order of magnitude in the hyperon case. This corroborates our assumption that smoothly-behaved, non-decreasing with µ F FFs (µ F -constant, in the toy case) stabilize cross sections. Furthermore, the flavor independence of toy FFs removes any possible modulation on parton densities, thus confirming that the gluon channel drives the growth with µ F of the convolution between PDFs and FFs. The main objection against our statement could be that, at NLA, the large contribution of c-and b-quark FFs for Λ c is not anymore dampened by the smallness of intrinsic heavy-quark PDFs. Indeed, the term proportional to the C gq non-diagonal coefficient that appears in the hadron NLO impact-factor correction (see Eq. (4.58) of Ref. [22]) is enhanced by large heavy-flavor FFs, whose sum is in turn multiplied by the gluon PDF. Thus, the production channel that opens up at NLA can in principle compete with the pure LO one and its effect could spoil the description presented above. We numerically checked, however, that in kinematic ranges typical of our investigation the C gg diagonal coefficient strongly prevails over the non-diagonal ones, its regular part being larger 50 to 10 4 times than C gq . Therefore, gluon dynamics still dominates at NLA and this confirms our statement that the peculiar behavior of the Λ c gluon FF is responsible for the stability of our distributions under higher-order corrections.

Summary and prospects
We studied the inclusive hadroproduction of a forward Λ c baryon in association with a backward object (another Λ c or a light jet) as a new probe channel of the semi-hard regime of QCD. Results for rapidity-distance distributions and azimuthal-angle correlations, calculated within a hybrid factorization that combines collinear ingredients inside a full NLA BFKL treatment, offer a promising statistics and exhibit standard features of the high-energy dynamics. We provided with a two-fold analysis on the sensitivity of energy scales, based on the variation of µ R (and µ F ) around their natural values as well as around the ones prescribed by the BLM optimization scheme. We discovered that the tag of Λ c particles allows us to dampen the instabilities of the BFKL series, this resulting in a partial stabilization of resummed distributions under higherorder logarithmic corrections. This effect is more pronounced in the double Λ c production, while further studies on the Λ c plus jet channel will gauge the dependence of our observables on intrinsic features of the description of jet emissions, such as the selection algorithm. We plan to extend our program on semi-hard phenomenology by hunting for stabilizing effects via the inclusive production of heavier particles, such as Λ b baryons, heavy-light mesons and quarkonium states (see Refs. [123,124] for a selection of phenomenological prospects in wider kinematic ranges, and Ref. [125] for applications at the EIC). proton(p 1 ) + proton(p 2 ) → H(p H1 , y H1 ) + X + H(p H2 , y H2 )    Figure 4: Behavior of azimuthal-correlation moments, R nm ≡ C n /C m , as functions of ∆Y , in the double Λ c channel, at natural scales, and for √ s = 13 TeV. Error bands provide with the combined uncertainty coming from scale variation and numerical integrations. Text boxes inside panels show transverse-momentum and rapidity ranges.   Figure 5: Behavior of azimuthal-correlation moments, R nm ≡ C n /C m , as functions of ∆Y , in the double Λ c channel, at BLM scales, and for √ s = 13 TeV. Error bands provide with the combined uncertainty coming from scale variation and numerical integrations. Text boxes inside panels show transverse-momentum and rapidity ranges.   proton(P 1 ) + proton(P 2 ) → Λ c (p Λ , y Λ ) + X + jet(p J , y J ) LLA NLA Figure 6: Behavior of azimuthal-correlation moments, R nm ≡ C n /C m , as functions of ∆Y , in the Λ c plus jet channel, at BLM scales, and for √ s = 13 TeV. Error bands provide with the combined uncertainty coming from scale variation and numerical integrations. Text boxes inside panels show transverse-momentum and rapidity ranges.    Figure 9: Behavior of the ϕ-summed cross section, C 0 , as a function of ∆Y , in the dihadron production channel, and for √ s = 13 TeV. A study on progressive energy-scale variation in the range 1 < C µ < 30 is done for Λ c emissions (upper left), for Λ detections (upper right), and for the production of a generic hadron species (lower) described by toy-FF parameterizations.