SciPost Submission Page
Colour unwound  disentangling colours for azimuthal asymmetries in DrellYan scattering
by Daniël Boer, Tom van Daal, Jonathan R. Gaunt, Tomas Kasemets, Piet J. Mulders
This Submission thread is now published as
Submission summary
Authors (as registered SciPost users):  Jonathan Gaunt · Tom van Daal 
Submission information  

Preprint Link:  http://arxiv.org/abs/1709.04935v3 (pdf) 
Date accepted:  20171214 
Date submitted:  20171104 01:00 
Submitted by:  van Daal, Tom 
Submitted to:  SciPost Physics 
Ontological classification  

Academic field:  Physics 
Specialties: 

Approach:  Theoretical 
Abstract
It has been suggested that a colourentanglement effect exists in the DrellYan cross section for the 'double Todd' contributions at low transverse momentum $Q_T$, rendering the colour structure different from that predicted by the usual factorisation formula [1]. These Todd contributions can come from the BoerMulders or Sivers transverse momentum dependent distribution functions. The different colour structure should be visible already at the lowest possible order that gives a contribution to the double BoerMulders (dBM) or double Sivers (dS) effect, that is at the level of two gluon exchanges. To discriminate between the different predictions, we compute the leadingpower contribution to the low$Q_T$ dBM cross section at the twogluon exchange order in the context of a spectator model. The computation is performed using a method of regions analysis with Collins subtraction terms implemented. The results conform with the predictions of the factorisation formula. In the cancellation of the colour entanglement, diagrams containing the threegluon vertex are essential. Furthermore, the Glauber region turns out to play an important role  in fact, it is possible to assign the full contribution to the dBM cross section at the given order to the region in which the two gluons have Glauber scaling. A similar disentanglement of colour is found for the dS effect.
Author comments upon resubmission
List of changes
1. The origin of the entanglement in [1] stems from the form of the link structure assumed (in eq. (2) and fig. 2 of [1]), which is based on rules given in a paper by Bomhof, Mulders, and Pijlman (Eur.Phys.J. C47 (2006) 147162), which were not derived for the case of dBM or dS. For these cases, the rules would lead to a colour entanglement as in [1], just as for instance for quarkquarkgluon correlations (with nonzero gluon momentum) in the situation of double twistthree contributions, although that situation may have other problems with factorisation. It is precisely the fact that one is dealing with zeromomentum gluons that leads to a larger set of diagrams than the set accounted for in [1], as demonstrated in the present paper. We have now pointed this out in the Introduction.
2. What the CSS proof shows is the cancellation of finalstate poles that prevent deformations out of the Glauber region into the soft/collinear regions. After the finalstate pole cancellation, effects from the Glauber region of momentum do still remain, which may then be absorbed into the soft and collinear functions (i.e. the TMDs). This point is emphasised in section 3. Moreover, in various places in the text, we have emphasised that even though the contribution at the order considered can be viewed as coming entirely from the Glauber region, it can also be absorbed entirely into the standard BM TMDs  see the sentence beginning with “As a byproduct...” in the Introduction, a similar sentence in the Conclusions, and the detailed discussion at the end of section 4.1. The apparent contradiction thus comes from the (often used) incorrect statement that Glauber \textit{contributions} cancel.
3. We have added the reference and reformulated the sentence, also incorporating point 5. It now reads: “To make matters worse, a recent analysis suggested that also in the DY process ‘colourentangled’ contributions can arise [1], i.e. contributions that, at best, come in a factorised form with a colour factor different from that predicted by the factorisation theorem. The affected contributions involve two Todd TMDs, such as the BoerMulders (BM) function [13] and the Sivers function [14,15].”
4. We have added references here.
5. See point 3.
6. We have reformulated this sentence to read: “Factorisation of DY scattering into PDFs and a perturbatively calculable hard factor was established by Collins, Soper, and Sterman (CSS) during the eighties in [33–35] and fully completed in 2011 also for the TMD case in [4].”
7. We have added the following sentence after eq. (6): “Eq. (5) is not in fact the full definition of the TMD – one has to accompany the bilocal matrix element by a soft factor that removes rapidity divergences and avoids double counting between the TMDs (see [4] and section 3). We do not consider this soft factor further here, however, as it will not appear in our model calculation in section 4.”
8. We were indeed a bit sloppy here; the definition of the Wilson line is now more explicit, see eqs. (7)–(9).
9. We have reformulated this sentence to read: “Although the TMD factorisation proof in [4] was intended to apply to all polarised cases, no explicit proton or quark polarisation was in fact considered. The colourentanglement effect in [1] would signal a loophole in the proof for the double Todd contributions that involve polarisation.”
10. We have modified the paragraph to state that the pinch surfaces only appear in the massless limit.
11. We fully agree that it makes more sense to bring the two lists together, and have done so. We have also added references.
12. We have added the ultrasoft scaling along with the soft scaling to the list of regions, noting that these are typically treated together in the CSS methodology and citing [2] and references therein.
13. In our paper we now write explicitly that “$\lambda$ is a small parameter which should in practice be of order $\Lambda/Q$”, hoping that no reader will assume it is an integration parameter.
14. Since the paper is entirely focussed on the TMD case, the discussion of the collinear case has simply been removed now.
15. We have changed the GrammerYennie approximation to the one of [4] – see point 16. We mention one key place in which the GrammerYennie approximation fails for the Glauber region (see the discussion at the end of p.9). At the top of p.10, we have explained why the sign of $i\epsilon$ in the denominators has been chosen as it has. We hope that this is sufficiently clear to the reader.
16. We have now altered the discussion to use the GrammerYennie approximation of [4]. This brief discussion of the GrammerYennie approximation is only there to illustrate how softcollinear and collinearhard attachments can be stripped away by Ward identities, but that the same procedure does not work for Glauber gluons. Eq. (17) is not used anywhere else in the paper, and in our opinion, a rather detailed technical discussion of the relative merits of the prescriptions in [2] and [4] would detract from the flow of the paper.
17. We have inserted ‘dBM’ here.
18. We have added a sentence on p.12 (at the end of the second paragraph) pointing out that if $\ell_1$ or $\ell_2$ is soft or ultrasoft in figures (i)–(iii), the contribution to the dBM cross section is power suppressed. The region where $\ell_1$ is ultrasoft and $\ell_2$ is Glauber can contribute at leading power in diagram (v), but this contribution also cancels  we have added some discussion of this in appendix B.
19. We have added a paragraph in section 4.1 noting that this paper involves the first application of the Collins subtraction method in which the Glauber region is treated in a distinct way. There is some work along similar lines that makes use of other subtraction schemes – this work is cited in the added paragraph.
Published as SciPost Phys. 3, 040 (2017)
Reports on this Submission
Anonymous Report 1 on 2017119 (Invited Report)
 Cite as: Anonymous, Report on arXiv:1709.04935v3, delivered 20171109, doi: 10.21468/SciPost.Report.274
Report
Generally, the authors have dealt with my concerns properly in the revised version. However, there are still a few areas where I still find problems.
1. The response to item 1, includes the following: "... the (often used) incorrect statement that the Glauber contributions cancel". I don't see that this is an incorrect statement, at least for a certain valid use of the terminology, as I will try to explain below. But there is clearly a different usage where one can say that the Glauber contributions do not cancel. The view with the cancellation is essential to deriving factorization, while the view with the noncancellation is useful for calculations such as are made in Sec. 4.
Now, a key point in the LibbySterman analysis is that, in determining the infrared regions to be considered in a graph, it is only necessary to consider those regions where the integration contour is trapped. If one deforms the contour away from regions without a pinch, contributions to a process only arise from regions corresponding to pinches. In individual cut graphs for the DrellYan process there are pinches in the Glauber region. One can try moving the contour away from pinches involving Glauber momenta, but only at the price of extra contributions from finalstate poles that are crossed. These contributions can have Glauber momenta, but they cancel to leading power after a sum over cuts. What is left is an integral over the deformed contour, which does not go through Glauber momenta. With this view, which is correctly summarized in the paper, I think it is correct to say that the Glauber contributions cancel (as always with the proviso "to leading power").
A different point of view is also possible, which leads to the calculational methods devised in the paper with its explicit use of contributions from the Glauber region, and associated subtractions.
I think a resolution can be motivated by the situation in the BrodskyHwangSchmidt (BHS) calculation of a onegluon contribution to the Sivers asymmetry in SIDIS. To get the asymmetry an imaginary part in an amplitude is needed, and BHS calculate this by setting the appropriate intermediate state onshell. In the calculation, the gluon is in the Glauber region, and one can therefore say that the asymmetry is given by the Glauber contribution. (Essentially, the same idea appears in calculation in Sec. 4.2 in the present paper for the BM function.)
But in an analysis to derive factorization, one can deform the integration out of the Glauber region into the collinear region for the gluon, and thereby derive factorization. The imaginary part now comes from the fact that the integral over the longitudinal momentum fraction of the gluon goes through complex values to avoid a pole. In this view, there is no Glauber contribution, only a collinear contribution.
Therefore, unless I have misunderstood something, the above argument leads to the need for the paper to contain further explanation of these issues. I am pretty sure that these different views are known and understood by the authors. But making them and their relations as fully explicit as possible would enable readers to clearly understood them as well, and would considerably assist future work.
2. Concerning the justification that the soft and ultrasoft regions do not need to be considered for the dBM cross section: The authors have added a statement making the assertion that these regions are power suppressed. I think it would be very helpful also to add a summary of how this is known, at least if this can be done briefly. A couple of sentences may suffice. In addition, I now see that some similar remarks are needed for the calculation of the BM function in Sec. 4.2. These changes would improve the reproducibility of the results. I suspect that readers who are not intimately familiar with the issues involved will probably have great difficulty verifying this result in the absence of suitable hints. They may well know that soft and/or ultrasoft contributions are generally leading (and even divergent) in the graphs considered, but may not find it obvious that the situation changes for the BM/dBM case.
3. A related issue, particularly prominent for the calculation of the BM function, is about the need for a collinear term in (31). (A generalization is also relevant for the calculation of the dBM cross section.) The final contribution to the BM function at this order is only from the Glauber region. If there is no collinear contribution, why bother with the collinear term in (31)? After all, there are no (leadingpower) soft and ultrasoft contributions, so no corresponding terms are used in (31).
I think the reason is that the power suppression for the soft and ultrasoft regions in the BM function happens graphbygraph, whereas the collinear contribution is unsuppressed for individual graphs. The collinear contribution only cancels after adding the hermitian conjugate graph, so in a graphbygraph analysis it must be retained. In addition, the use of the collinear term is needed to allow a matching by contour deformation with the factorization point of view with its attribution of the BM function to the collinear region with complex momentum fraction. Could the authors confirm whether or not my understanding is correct? If so, it would be useful to add something to Sec. 4.2 (and/or elsewhere) to explain near the beginning why they use formulas containing collinear terms even though the final answer is only from Glauber contributions.
4. I am still rather confused by the explanation of how the derivations in Ref. [1] came to be in error. But since that issue concerns a different paper, I will not pursue it further.
5. On p. 6, it is stated that the TMD factorization proof in [4] did not explicitly include proton or quark polarization effects. It is indeed true that the chapter that treats the DY case seems to restrict its mention of these effects to inserting the words "polarization effects" in factorization formulas. However, the chapter should be read in the context of previous chapters, where polarization effects are explicitly treated. The relevant parts include (a) Sec. 11.9 for the Ward identities, especially Sec. 11.9.5 for the generalization beyond DIS (since Ch. 11 is about DIS); (b) Chs. 12 and 13. There the methods explicitly include polarization; See, for example, (13.22), where projectors on Dirac spinor space are applied that enable all leading power polarization effects to be included. The application of Ward identities for collinear gluon attachments to hard subgraphs and for soft gluon attachments to collinear subgraphs is such that they factorize independently of polarization. (Earlier chapters contain the details.)
It is surely intended to be taken for granted that the same procedures are to be applied to DrellYan. The extra step for DY is the demonstration of cancellation of final state terms, in Sec. 14.4. There is no need to explicitly mention polarization, since all necessary polarization dependent terms are inside the various shaded subgraphs.
Of course, it is always possible there is a mistake; the treatment of the final state poles is very nontrivial, and the dBM and dS processes are particularly sensitive to any problems. Certainly it would be nice to improve the proof. But in identifying potential sources of breakdown or improvement of the proof it is important to correctly identify what was done and not done. If after considering my remarks, the authors are still concerned that issues concerning polarization are missed, it would help for them to specify more precisely what their concern is, if that is possible.
6. In addition, I have noticed the following minor issue: In the second paragraph of Sec. 2, the references for proofs of factorization in [3335] include the Bodwin paper, but the only authors mentioned are CSS. Either include Bodwin in the list of people or remove also CSS.
Requested changes
See report.
Author: Tom van Daal on 20171121 [id 194]
(in reply to Report 1 on 20171109)We would like to thank the referee for his/her very careful reading of the manuscript and providing valuable feedback. Below we comment on each one of the points raised and provide suggestions for changes where needed.
Currently we have below eq. (38) in section 4.2:
“Inasmuch as the function $\chi^j$ is real, only the imaginary part of the $\ell_1^+$ integral contributes to $h_1^\perp$ as its real part is canceled by the Hermitian conjugate term. This imaginary part comes from the region where $\ell_1^+$ is sensitive to the $i\epsilon$ term in the denominator, which is the case when $\ell_1^+ \to 0$ – i.e. when $\ell_1$ has Glauber scaling.”
We would like to add to this the following sentence containing a reference to the BHS papers:
“Note that similar arguments were used in \cite{Brodsky:2002cx,Brodsky:2002rv} to obtain singlespin asymmetries.”
“Note that none of these regions involve a soft or ultrasoft scaling for either $\ell_1$ or $\ell_2$ – if either $\ell_1$ or $\ell_2$ is soft or ultrasoft, then the contribution to the dBM cross section from the graph is power suppressed.”
We agree with the referee that a brief expansion on this statement would be helpful. Hence, we propose to add the following:
“In the case in which $\ell_1$ or $\ell_2$ is soft, the graphs become power suppressed as too many quark lines are brought off shell to virtualities of order $\Lambda Q$ by the soft momentum. The same power suppression would also hold for these graphs in the unpolarised case. By contrast, the power suppression of the graphs when $\ell_1$ and/or $\ell_2$ is ultrasoft is specific to the spindependent case – here the suppression occurs in the numerator traces of eq. (22).”
In section 4.2, just above eq. (31), we currently write:
“To calculate the BM function, we first identify two nontrivial momentum regions for the gluon momentum $\ell_1$ that give a leadingpower contribution, namely $G_1$ and $C_1$. Summing over these regions gives, according to eqs. (19) and (20),”
We propose to change this into the following:
“To calculate the BM function, we first identify two nontrivial momentum regions for the gluon momentum $\ell_1$ that give a leadingpower contribution, namely $G_1$ and $C_1$. The $S$ and $U$ regions give powersuppressed contributions for the same reasons as discussed earlier in section 4.1. Summing over the leading regions gives, according to eqs. (19) and (20),”
Given this, we suggest to include the following after the added paragraph for point (2):
“Note that ultimately we will see that the $C_1G$, $GC_2$, and $C_1C_2$ regions also vanish at leading power. However, this happens in a highly nontrivial way only after the sum over graphs and possible finalstate cuts, and only when the appropriate subtraction terms for smaller regions are included. Furthermore, this is related to the rapidity regulator that we use (see discussion below). Thus, we consider these regions explicitly here, detailing how and why this cancellation happens.”
Okay.
We realise that we were not completely accurate here and agree with the referee. Currently, we state on p.6:
“Although the TMD factorisation proof in [4] was intended to apply to all polarised cases, no explicit proton or quark polarisation was in fact considered. The colourentanglement effect in [1] would signal a loophole in the proof for the double Todd contributions that involve polarisation.”
We propose to change this into the following:
“The colourentanglement effect in [1] would signal a loophole in the TMD factorisation proof of [4] for double Todd contributions that involve polarisation.”
Additionally, in the first sentence of the first complete paragraph on p.17, we should omit the word “unpolarised”.
Furthermore, on p.26 we currently have:
“We can conclude that the dBM cross section factorises in the same way as the double unpolarised contribution that was considered in the original CSS factorisation proof.”
We propose to change this sentence into the following:
“We can conclude that the dBM cross section precisely factorises as already anticipated by the CSS works; no loophole in their original proof for this double Todd contribution is found.”
“Factorisation of DY scattering into PDFs and a perturbatively calculable hard factor was established by Collins, Soper, and Sterman (CSS) during the eighties in [33–35] and fully completed in 2011 also for the TMD case in [4].”
We propose to replace this by the following:
“Factorisation of DY scattering into PDFs and a perturbatively calculable hard factor was established by Collins, Soper, and Sterman (CSS) during the eighties in [33–34], with important work in this direction also being done by Bodwin [35]. The factorisation proof for the TMD case largely proceeds along the same lines and is covered in [4].”
Furthermore, we currently have at the end of section 3:
“The final step of the factorisation proof is the partitioning of the soft subgraph between the two collinear subgraphs – once this is done one has the final factorised form with two TMDs and a hard function, as in eq. (1).”
We suggest to replace this by the following:
“The final step of the factorisation proof is the partitioning of the soft subgraph between the two collinear subgraphs, for which recently an allorder proof was provided in \cite{Vladimirov:2017ksc}. The result of this procedure is a factorised form with two TMDs and a hard function, as in eq. (1).”