Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: single - transverse spin asymmetries ( ssas ) in high - energy collisions with one transversely polarized hadron are important phenomena that have been observed for more than three decades in various physical processes @xcite . in these processes , the observed final - state hadrons show an asymmetric distribution in a plane perpendicular to the beam direction when the transversely polarized hadron scatters off an unpolarized hadron ( or a virtual photon ) . the ssa is defined as @xmath2 , the ratio of the difference and the sum of ( differential ) cross sections when the hadron s spin vector , @xmath3 , is flipped . recent experimental measurements of ssas in polarized semi - inclusive lepton - nucleon deep inelastic scattering ( sidis ) @xcite and in hadronic collisions @xcite have renewed the theoretical interest in ssas and in understanding their roles in hadron structure and quantum chromodynamics ( qcd ) . although it was realized a long time ago @xcite that perturbative qcd can be used to study the effects of transverse spin , the size of the observed asymmetries came as a surprise and has posed a challenge for researchers in this field @xcite .    within a model calculation @xcite , brodsky , hwang , and schmidt showed that the final state interaction in deep inelastic scattering ( dis ) can generate a phase required for a nonzero ssa in sidis . it was later realized that this final state interaction can be factorized into the gauge link of the gauge invariant transverse momentum dependent ( tmd ) quark distributions @xcite . the nonvanishing ssa obtained in ref . @xcite is a consequence of the existence of a naively time - reversal - odd tmd quark distribution , the so - called sivers function @xcite . when applying the same calculation to the drell - yan production of lepton pairs at hadron colliders , the final state interaction in sidis becomes an initial state interaction in drell - yan , and the phase changes sign , which leads to a prediction of a sign change in the ssas between these two processes @xcite . this nontrivial `` universality '' property associated with the tmd parton distributions is the consequence of gauge interactions in qcd @xcite , and of the qcd factorization theorems for these two processes @xcite . experimental tests of this prediction will be crucial for our understanding of the origin of ssas in qcd  @xcite .    in ref . @xcite , it was proposed to study the sivers functions by means of a ssa in azimuthal - angular correlations of two jets produced nearly back - to - back at hadron colliders . measurements of this ssa for dijet production have begun at rhic @xcite , complementing the measurements in sidis . unlike the sidis or drell - yan process , dijet production at hadron colliders involves both initial and final state interactions that may produce the phase needed for a nonvanishing ssa . consequently , the sign and the size of the asymmetry will depend on the relative strength of these interactions . following the previous works on sidis and drell - yan , the authors of ref . @xcite developed a systematic approach to describe the role of initial- and final - state interactions in generating ssas in hadronic collisions , and they found that summing all initial / final state interactions into the gauge link of the tmd parton distributions leads to a very complicated functional form of the gauge link . in particular , the tmd parton distributions studied in dijet correlations in hadronic scattering will have no connection to those in the sidis and drell - yan processes , because their definitions are completely different @xcite . one then has to question the universality of the tmd parton distributions , and the predictive power of perturbative qcd calculations , which relies on comparing physical observables with the same factorized long - distance physics while having different perturbatively calculable short - distance dynamics .    in ref . @xcite , we briefly reported a new result for the ssa in dijet production in hadronic collisions in the twist-3 efremov - teryaev - qiu - sterman ( etqs ) approach @xcite . we considered the spin - dependent cross section , @xmath4 , for the process @xmath5 with the jet momenta @xmath6 and @xmath7 . when both @xmath8 and @xmath9 , the transverse components of the momenta @xmath10 and @xmath11 , respectively , are much larger than @xmath12 , a nonvanishing @xmath13 can be generated by the etqs mechanism in the collinear factorization approach . in this framework , the ssas are attributed to the spin - dependent twist - three quark - gluon correlation functions , which correspond to a quantum interference between different partonic scattering amplitudes . since the incoming partons are approximated to be collinear to the corresponding incoming hadrons in this approach , the momentum imbalance of the two jets is generated by producing a three - parton final - state . @xcite , we calculated the contribution from initial - state gluon radiation to @xmath13 in the kinematic region where @xmath14 . we presented the final result for the leading contributions to @xmath13 in the expansion of the partonic scattering in @xmath15 involving a hard @xmath16 subprocess . in this paper , we will provide the detailed derivations of this result , and we will also present the full contributions from all other partonic subprocesses at the same order . we find that although both initial - state and final - state interaction lead to a nonvanishing ssa , the final - state interactions give the dominant contributions to @xmath13 . we therefore expect that the ssa in dijet production will have the same sign as the sivers asymmetry in sidis . we find that at leading order in the @xmath15 expansion , the perturbatively calculated partonic parts can be further factorized into a single - scale ( @xmath8 ) hard part and perturbatively generated tmd parton distributions with transverse momenta @xmath17 . we also find that our perturbatively calculated result is equal to the leading - order term in the @xmath18 expansion of the following generalized tmd factorization formula @xcite : @xmath19_c\ , \delta^{(2)}(\vec{k}_{1\perp}+\vec{k}_{2\perp}+ \vec{\lambda}_\perp-\vec{q}_\perp ) \ , , \nonumber\end{aligned}\ ] ] where @xmath20 runs over all parton flavors , @xmath21 and @xmath22 are partonic hard and soft factors , respectively , and the @xmath23_c$ ] represents a trace in color space between the hard and soft factors due to the color flow into the jets @xcite . the hard factor in eq . ( [ e4 ] ) only depends on the single hard scale @xmath8 in terms of partonic mandelstam variables of the reaction @xmath24 : @xmath25 with @xmath26 , @xmath27 and @xmath28 and @xmath29 the rapidities of the two jets . in eq . ( [ e4 ] ) , @xmath30 and @xmath31 denote the transverse - spin dependent tmd quark distributions ( known as the sivers function ) and the unpolarized tmd quark distribution , respectively ; these tmd parton distributions were chosen to follow their definitions in the semi - inclusive dis process . for example , for a polarized proton with momentum @xmath32 with @xmath33 and transverse spin vector @xmath34 , the tmd distribution for quark flavor @xmath35 can be defined through the decomposition of the following matrix element @xcite , @xmath36 \ , \label{e3}\end{aligned}\ ] ] where @xmath37 is a hadronic mass scale introduced to keep the tmd distributions @xmath38 and @xmath39 at the same mass dimension , and the gauge link @xmath40 is defined in a covariant gauge as @xmath41 with the path link extended to @xmath42 . @xmath43 is a vector conjugate to the momentum vector @xmath10 . since we will work in a covariant gauge throughout this paper , the vector @xmath43 could be chosen to be a light - cone vector with @xmath44 and @xmath45 . an off - light - cone vector ( @xmath46 ) will have to be used when high order corrections are taken into account  @xcite . if we work in a singular gauge , like the light - cone gauge , an additional gauge link at the spatial infinity ( @xmath47 ) will have to be included in order to ensure the gauge invariance of the above definitions @xcite . we have chosen tmd parton distributions defined in sidis because of the dominance of final - state interactions . choosing tmd parton distributions defined according to the drell - yan process would change the sign of the partonic hard factors , but not affect the overall sign of the physical cross section . we note that the various factors in eq . ( [ e3 ] ) , apart from being functions of transverse momentum , also depend on the renormalization and factorization scales , and especially on the gluon rapidity cut - off when higher - order corrections are taken into account  . the latter dependence is governed by the collins - soper evolution equation  @xcite , which leads to a resummation of large logarithms of the form @xmath48 $ ] in the perturbative series  @xcite .    in order to evaluate the ssa for dijet production , we also calculate in this paper the contributions from initial - state gluon radiation to the spin - averaged dijet cross section , @xmath49 , in the same kinematic region where @xmath50 . using the collinear factorization approach , we find that like the spin - dependent case , the leading contribution in the @xmath15 expansion of the perturbatively calculated partonic scatterings can be further factorized into a hard part at a single - scale @xmath8 and perturbatively generated unpolarized tmd parton distributions with transverse momenta @xmath17 . our leading contribution in the @xmath15 expansion is equal to the leading - order term of the following tmd factorization formula @xcite , @xmath51_c \ , \delta^{(2)}(\vec{k}_{a\perp}+\vec{k}_{b\perp }              + \lambda_\perp-\vec{q}_\perp ) \ ; , \label{e2}\end{aligned}\ ] ] where the superscript `` @xmath52 '' indicates the scattering of an unpolarized beam off an unpolarized target / beam . for consistency , we again express the factorization formula in terms of tmd parton distributions defined in the sidis process . actually , the unpolarized tmd parton distributions defined in sidis and the drell - yan process are identical , because the unpolarized parton distributions are invariant under the naive - time - reversal transformation . a key feature of any qcd factorization is that the perturbatively calculated short - distance hard factors should not be sensitive to details of the _ factorized _ long distance physics . in this paper , we will derive all short - distance hard factors in eqs . ( [ e4]),([e2 ] ) also by using the brodsky - hwang - schmidt model for ssas @xcite . we find that the hard - scattering factors derived in this model are in fact the same as those derived in the collinear factorization approach , despite of the clear difference in the treatment of the nucleon . both the spin - averaged and the spin - dependent cross sections calculated in this model are therefore consistent with the generalized factorization formulas in eqs . ( [ e4]),([e2 ] ) .    to fully investigate the above factorization formalism , we would have to consider the contributions from gluon interactions and radiations in all possible regions of the phase space , and to all orders in perturbation theory , which is not what we are trying to do in this paper . instead , as a step to test the above factorization formalism , we will study the contribution at the first non - trivial order with a single gluon radiated nearly parallel to one of the incident nucleons . from this study , we will be able to show how to factorize this gluon into the tmd parton distributions of the incoming nucleons , and to verify the definition of the tmd parton distributions used in the generalized factorization formulas . our first order calculation at @xmath53 clearly shows the factorization of a hard part at @xmath54 and that the radiation of a collinear gluon can be absorbed into the relevant tmd parton distributions . however , our work does not address the factorization between the various tmd parton distributions at the scale @xmath55 . our calculation should be regarded merely as a `` necessary condition '' for such a factorization to hold . since the predictive power of perturbative - qcd calculations involving measured hadrons relies on factorization , a full proof or disproof of this factorization remains an important challenge @xcite . the rest of this paper is organized as follows . in sec . ii , we calculate the partonic hard factors in a simple model - inspired approach . using the quark - diquark model for the proton and the brodsky - hwang - schmidt model for ssas , we extract the hard factors by the factorization of the tmd parton distributions . we find that the hard factors for the unpolarized scattering in eq . ( [ e4 ] ) are identical to the partonic differential cross sections @xmath56 , and the ones for the transverse - spin dependent cross section are different due to the initial- and final - state color interactions . in sec . iii , we calculate the ssa in the collinear factorization approach , and show that the collinear gluons parallel to the polarized proton can be factorized into the sivers function , while those parallel to the unpolarized proton can be factorized into the unpolarized tmd parton distribution . we also demonstrate that the hard factors calculated in both approaches are the same and independent of the details of the factorized long - distance physics . in sec . iv , we extend our formulas to the @xmath9-weighted ssa , and compare with previous results in the literature . finally , we summarize our paper in section v. in this section , we will calculate the hard factors at @xmath57 in the generalized factorization formulas in eqs . ( [ e4]),([e2 ] ) , based on a model in which an energetic parton scatters off a nucleon that is made of a quark and a scalar diquark @xcite . we will follow the approach of brodsky et al . @xcite and incorporate the correct color interactions between the gauge boson and the scalar diquark in the proton wave function . we extract the hard factor by factorizing the model - sensitive long - distance physics from the differential cross section into the parton distributions . because the model - dependence was factorized into the parton distributions , the extracted hard factors resulting from our calculations do not depend on how the nucleon couples to the partons , and therefore , should be model - independent . " represents the color - index for the gauge boson , and @xmath58 the color - indices for the quark and/or the di - quark . _ , width=340 ]    -0.4 cm      in this subsection , we derive the hard factors @xmath59 for the unpolarized differential cross section in eq . ( [ e2 ] ) . we will first give a simple example of how the hard factors may be extracted in the context of the quark - diquark model of @xcite . we will also give a more general approach based on power counting techniques to calculate the hard factors directly from partonic diagrams .    to specify the quark - diquark model of the proton , we give in fig . [ f2 ] graphical rules for the vertex coupling the quark - diquark to the proton , and the vertex coupling the diquark to a gauge boson . in order to correctly take into account the color degree of freedom of the strong interactions , we included proper color factors for these vertices . using this quark - diquark model to extract the hard factors not only verifies the fact that the hard factors are independent of the dynamics at the scale of the proton , but also provides a relatively simpler demonstration of how to factorize the cross sections into the tmd quark distributions and corresponding hard factors at the leading order in the @xmath15 expansion .    in this particular model , we calculate the cross section for an energetic parton of momentum @xmath60 scattering off a parton of momentum @xmath61 from the ( modeled ) nucleon , producing two jets of momenta @xmath62 and @xmath63 plus an unobserved particle with momentum @xmath64 in the final state , as shown in fig . [ f3 ] , where the feynman diagram represents the contribution from a quark - quark scattering channel : @xmath16 with the quark @xmath11 from the nucleon . we extract the hard factors by comparing the calculated cross section with the factorized formulas in either eq . ( [ e4 ] ) or ( [ e2 ] ) . the diquark in the final state has a transverse momentum @xmath65 to balance the small momentum imbalance of the two jets , @xmath66 . channel , for the unpolarized cross section . _ , width=207 ]    -0.4 cm    according to the factorization formulas in eqs . ( [ e4]),([e2 ] ) , the transverse momentum imbalance @xmath9 of the two jets is a consequence of adding the transverse momenta of the two active partons and from the soft factor . since the incoming quark in this model calculation has only a longitudinal momentum component , its transverse as well as longitudinal momentum distributions should be given by delta functions , @xmath67 similarly , at this order , the soft factor is also a delta function of the transverse momentum , i.e. , @xmath68 . therefore , the only contribution to the transverse momentum imbalance of the two jets comes from the transverse momentum of the quark emerging from the nucleon , and the generalized factorization formulas can be reduced to @xmath69 for the unpolarized cross section , and @xmath70 for the single transverse - spin dependent cross section .    in order to extract the hard factor from the differential cross sections in eqs . ( [ e7]),([e8 ] ) , we will need to calculate the tmd quark ( unpolarized and sivers ) distributions within the same model . they are given by the diagrams shown in fig . [ f4 ] , and are available @xcite . we can summarize their results as follows : @xmath71 \ , \label{e9}\\ q_t^{\rm sidis}(x , k_\perp)&=&c_f\frac{g_n^2\alpha_s}{16\pi^3}\frac{(1-x)(m_q+xm_p ) } { \lambda^2(x , k_\perp^2)}\frac{m_p}{k_\perp^2 } \ln\frac{\lambda^2(x , k_\perp^2)}{\lambda^2(x,0 ) } \ , \label{e10}\end{aligned}\ ] ] where @xmath72 and @xmath73 , and @xmath74 , @xmath75 and @xmath76 are masses for the proton , the quark and the gauge boson , respectively . in the above formulas , we have included the color - factors arising from the vertices described in fig . [ f2 ] .        -0.4 cm    we now present details of the calculation of the contribution of the quark - quark scattering channel @xmath77 to the unpolarized scattering cross section . we also provide full results for all other channels , which can be derived similarly . the feynman diagram for @xmath77 has been shown in fig . [ f3 ] . in the limit @xmath78 , the two final - state jets have approximately the same transverse momentum , i.e. , @xmath79 . we define the following kinematic variables : @xmath80 where @xmath81 . from these , we immediately find the incident partons momentum fractions @xmath82 and @xmath83 , and the scalar diquark s momentum can be written as @xmath84 in the evaluation of the cross sections , we use the power counting analysis @xcite , keep only the leading power contributions and neglect all higher order corrections in @xmath85 , where @xmath86 represents any lower mass scale like @xmath9 , @xmath74 , @xmath75 and @xmath76 . in this limit , the differential cross section for fig . [ f3 ] can be written as @xmath87 is the scattering amplitude for the diagram . expanding the above delta function for the phase space integral of @xmath64 , we find that in the leading power contributions @xcite , @xmath88 where the `` plus '' distribution follows the usual definition @xcite . this delta function will help to simplify our calculations for contributions from different kinematic regions of @xmath64 . for example , if @xmath64 is parallel to @xmath89 , which means @xmath90 , we only have a contribution from the first term of the above expansion , which also implies @xmath91 in this limit . furthermore , the quark propagator in fig . [ f3 ] reads @xmath92 where @xmath93 is defined above after eq . ( [ e10 ] ) , and is order of @xmath94 . combining this propagator with the delta function expansion , we find that indeed only the first term in the delta function contributes . neglecting all higher order terms in @xmath85 , we obtain the amplitude squared for fig . [ f3 ] as @xmath95 \ , \ ] ] where @xmath96 is the color - factor , and the hadronic mandelstam variables @xmath97 , @xmath98 , and @xmath99 are defined as @xmath81 , @xmath100 , @xmath101 . in the limit @xmath91 , these variables can be related to the partonic mandelstam variables as : @xmath102 , @xmath103 , @xmath104 at @xmath105 , and we will have @xmath106 . substituting the above results into the differential cross section , we obtain , @xmath107 \frac{2\pi\alpha_s^2}{\hat s^2}\ , \frac{\hat s^2+\hat u^2}{\hat t^2}\ , \delta(x_b-1 ) \ .\ ] ] comparing this result with eq .  ( [ e7 ] ) , and identifying the tmd parton distribution @xmath108 defined in eq . ( [ e9 ] ) , we extract the partonic hard factor for the differential cross section from the diagram in fig .  3 as@xmath109 where @xmath110 are the partonic mandelstam variables defined before . this hard factor is just the well - known differential cross section @xmath111 for the partonic process @xmath16 . the hard factors for all other partonic channels can be calculated similarly , and they are all identical to the partonic differential cross section @xmath112 , and collected in appendix a.    instead of calculating the differential cross section directly , we now show how to use the power counting techniques @xcite to analyze the contribution from all partonic channels and to demonstrate the factorization of the hard factors . we then apply this approach to calculate the contributions to the spin dependent differential cross section . we again first consider the simple example of the @xmath16 channel . the jet transverse momentum @xmath8 is the only hard scale , while all other momentum scales , such as the imbalance @xmath9 of the two jets and all other mass scales @xmath74 , @xmath75 and @xmath76 , are relatively soft and much smaller than the hard scale . we expand the scattering amplitude squared in the small parameter @xmath15 . using power counting techniques , we identify the leading power contributions and match them to the momentum regions of the unobserved ( or integrated ) parton of momentum @xmath64 . in our current model , where a parton of momentum @xmath113 scatters off a nucleon of momentum @xmath89 , the leading power contribution is from the region when @xmath64 is parallel to @xmath89 . if @xmath64 is parallel to @xmath113 , the contribution is power suppressed , because the denominator of the quark propagator @xmath114 will be of order @xmath115 , instead of @xmath116 as is the case when @xmath64 is parallel to @xmath89 .     into the tmd quark distribution of the proton and the hard factor represented by the partonic scattering feynman diagram.__,width=491 ]    -0.4 cm    when the momentum @xmath64 is parallel to @xmath89 , the quark propagator in fig . [ f3 ] has a virtuality @xmath117 , which is much smaller than the hard scale @xmath118 . compared to the short distance physics in the partonic scattering taking place at the time scale @xmath119 , the incident quark of momentum @xmath120 is long - lived , and the contribution from this subprocess can be separated into two parts : the parton distribution part relevant at the soft scale @xmath121 and the hard partonic part at the scale @xmath122 , representing the long distance and short distance physics , respectively . we demonstrate this factorization in fig . [ f5 ] , where the contribution to the differential cross section from the process in fig . [ f3 ] is factorized into a parton distribution multiplied by a hard factor . we can further separate the hard factor into a partonic scattering amplitude squared stripped of any color , and the color - factor @xmath123 . in this way we can take into account the color decomposition of the subprocess in the case of unpolarized scattering . this is trivial for this particularly simple partonic channel , but becomes a convenience for the more complicated ones . after subtraction of the parton distribution from the differential cross section , the hard factor in fig . [ f5 ] can then be written as @xmath124 where @xmath125 represents a common factor and is introduced to simplify the notation , and @xmath126 . this of course agrees with the above result in eq . ( [ hqqp ] ) . most importantly , it demonstrates that the hard factor is model independent . it only depends on the partonic scattering cross section and the color - factor associated with the partonic diagram . it does not depend on the model we used to describe how the nucleon couples to the quark and diquark .          -0.4 cm    a non - vanishing single transverse - spin asymmetry requires initial / final state interactions generating a phase . because all initial and final partonic states relevant for dijet production are colored , both initial and final state interactions have to be taken into account for a complete result . in fig . [ f6 ] , we show all possible diagrams with gluon exchange between the diquark and the hard scattering part which may contribute to the ssa for the dijet - correlation , again for the case of an underlying @xmath16 process . the initial and final state interaction diagrams in fig .  [ f6](a - c ) contribute to a soft gluon pole , which corresponds to the quark sivers function contribution we are considering in this paper . however , the gluon interaction with the internal gluon propagator shown in fig . [ f6](d ) does not contribute to a soft gluon pole , but to a soft fermion pole @xcite . this soft fermion pole contribution is not related to the quark sivers function and will not be discussed in this paper . therefore , in the following calculations , we will only consider the contributions from the diagrams ( a)-(c ) . we will follow the same factorization approach that was used in last subsection for calculating the hard factors for the unpolarized cross section . we use the power counting method to factorize the additional gluon interaction between the hard scattering and the diquark into the sivers parton distribution for every partonic channel . as discussed above , the leading power contribution is from the region of phase space where the outgoing diquark momentum @xmath64 is parallel to the polarized proton @xmath89 : @xmath127 , @xmath128 and @xmath129 . to obtain a leading power contribution , the gluon momentum @xmath130 also needs to be parallel to @xmath89 . this is so because the diquark propagator in fig . [ f6 ] reads @xmath131 which will be power suppressed unless @xmath132 . this result can also be derived from the existence of a `` pinch '' singularity in the @xmath133 integral over the quark propagator of momentum @xmath134 and the diquark propagator of momentum @xmath135 @xcite . since the momentum @xmath130 of the gluon is parallel to @xmath89 , its polarization will also be proportional to @xmath89 , because the gauge - boson - scalar vertex is proportional to @xmath136 , as shown in fig . therefore , the gluon attaching to the initial- or final - state quark is longitudinally polarized , i.e. , its polarization is along its momentum . because of this , we can further decouple the gluon interaction with the external particles by using the eikonal approximation . for example , for the diagram with final - state interaction in fig . [ f6](b ) , the gluon interaction part can be reduced to @xmath137 where the gamma matrix @xmath138 appears because of the interaction with a longitudinally polarized gluon , and where @xmath35 is the color - index for this gluon . in the above derivations , we have only kept the leading power contributions .     for ( 1,2,4,5 ) where @xmath35 is the color index for the gluon with momentum @xmath130 and @xmath58 are indices for the quark lines ; @xmath139 for ( 3,6 ) where @xmath140 are indices for incoming and outgoing gluons.__,width=529 ]    -0.4 cm    similarly , the initial state interaction in fig .  [ f6](a ) and the other final state interaction in fig .  [ f6](c ) can be simplified following the same eikonal approximation . the eikonal approximation will give rise to eikonal propagators and vertices . in fig . [ f7 ] , we summarize these eikonalized expressions for all the initial and final state interactions , including the attachments of the longitudinally polarized gluon to a quark or anti - quark line , or to a gluon line as well . the eikonal approximation for quark- and antiquark - gluon interactions with incoming or outgoing momentum can be derived following the above example . for the three - gluon interactions , we can derive the eikonal vertex and propagator by choosing physical polarizations for the external gluons . with this choice , we find the vertex and gluon propagator can be reduced to , for example for fig .  [ f7](6 ) : @xmath141\frac{i(-g^{\nu\nu'})}{(p_b+l)^2+i\epsilon}g^{\rho -}\nonumber\\&&\approx   \frac{(2p_b+l)^- g_{\mu\nu}}{(p_b+l)^2+i\epsilon}(-if_{abc})=\frac{g}{-l^+-i\epsilon } g_{\mu\nu}(-if_{abc } ) \ , \end{aligned}\ ] ] where we have used the fact that the incoming gluon is transversely polarized , so that the polarization tensor can be chosen as @xmath142 similarly , we can simplify the longitudinally polarized gluon attachment to the outgoing gluon line as shown in fig . [ f7](3 ) . these eikonal propagators and the associated vertices can be absorbed into the parton distribution for the polarized proton , after factorizing the diagram . the only difference between this factorization and the one for the unpolarized scattering is the color - flow , because of the additional gluon exchange between the diquark and the hard part . the color - flow is sensitive to the difference between initial- and final - state interactions and depends on the specific partonic subprocess . therefore , the generalized factorization would immediately fail if the eikonal propagators and associated vertices , which are absorbed into the tmd distributions , were sensitive to the process dependent color flows . however , at the leading order , one can always attribute the color - flow to the partonic hard part and keep the corresponding tmd distribution insensitive to it . therefore , as compared to the unpolarized case ( see fig .  [ f5 ] ) , for the single - spin - dependent cross section we shall obtain a more complicated factorization . for the case of the three diagrams of fig . [ f6 ] the factorization is shown in fig . the contribution from each diagram can be factorized into the sivers function and the hard partonic scattering amplitude squared , multiplied by a color - factor that takes into account the color - flow of the initial / final state interaction . however , in order to verify the generalized factorization beyond the leading order , one would have to show that additional initial - state as well as final - state gluon interactions could also be factorized into a sum of terms , each of which is factorized into a product of higher order tmd distributions and corresponding color factors and partonic hard parts . for example , to verify the factorization at the next - to - leading order with two eikonal gluon interactions , one would have to show that the initial- and final - state interactions can be factorized into the same color factors and the lowest order hard parts as in fig . 8 , multiplied by the corresponding second order tmd distribution . we will leave the detailed study of the generalized factorization at higher orders to a future publication .    . _ , width=529 ]    -0.4 cm    at least to the leading order , the extracted color factors @xmath143 with @xmath144 in fig . [ f8 ] are insensitive to the long - distance physics . that is , they do not depend on how the quarks and the gluon couple to the nucleon . this is because the factorization of the color - flow between the hard partonic part and the tmd part results in a color configuration connecting an incoming and an outgoing quark ( with color indices @xmath145 ) and a gluon ( with color @xmath35 ) . the only possible color structure to describe this connection is @xmath146 , due to @xmath147 invariance and the wigner - eckart theorem . therefore , one can reduce all the color matrices from the hard part to a simple color matrix @xmath146 , and the coefficient will be exactly the color - factor in the factorization formula . in this way , the extracted color - factor can not depend on the bottom part of the diagram describing the coupling of the quarks and the gluon to the nucleon , but will only depend on the color - flow in the hard partonic diagrams with an additional gluon insertion . for example , the color factors appearing in fig . [ f8 ] can be derived in the following way , @xmath148 where the color - matrices on the left side of the equations come from the partonic diagrams with one gluon attachment to the initial- or final - state quark line . the extracted color - factor @xmath149 is for the initial - state interaction , and @xmath150 and @xmath151 are for the final - state interactions on the lines with the momentum @xmath62 and @xmath63 , respectively .    in the factorization shown in fig . [ f8 ] , the leading order sivers functions for all three diagrams are the same . this is because we normalize all eikonal vertices to the one used in the calculations of the sivers function in sidis ( shown in fig . the difference from the initial or final state interaction effects is summarized into the relevant color - factors . for other partonic channels , it is in some cases necessary to introduce an extra sign when the eikonal propagator contributes with an opposite phase compared to that in the sivers function in sidis . for example , the eikonal propagators in the diagrams ( 2,5,6 ) of fig . [ f7 ] contribute an opposite sign , whereas those of ( 1,3,4 ) contribute the same sign as in sidis . because all the initial and final state interactions in fig . [ f8 ] contribute the same sign , there is no sign change for their color - factors . another important point is that the partonic scattering amplitude squared is the same for all the three diagrams in fig . it can be calculated directly from the relevant feynman rules for the partonic diagrams and is also the same as the one for unpolarized scattering given above . after factorizing out the sivers function from the cross section contributions for the diagrams in fig . [ f8 ] , the color factors for the three diagrams can then be summed to give the final result for the hard factor for the @xmath16 subprocess : @xmath152h_{qq'\to qq'}= \frac{\alpha_s^2\pi}{\hat s^2}\frac{n_c^2 - 5}{4n_c^2}\frac{2(\hat s^2+\hat u^2)}{\hat t^2 } \ , \label{hqq}\ ] ] where @xmath153 has been defined above after eq . ( [ hqu ] ) . the above derivations can be extended to all other partonic channels involving the quark sivers function . they can also be extended to the case of sub - processes initiated by the gluon sivers function for which , however , more color - structures ( like @xmath154 and @xmath155 for the three gluon coupling to the nucleon state ) will emerge @xcite . in this paper , we will focus on the processes with a quark sivers function , which include all quark - quark and quark - antiquark scattering processes , such as @xmath16 calculated above , @xmath156 , @xmath157 , @xmath158 and @xmath159 , as well as the quark - gluon process @xmath160 and quark - antiquark annihilation into a gluon pair , @xmath161 . following the above derivations for the @xmath16 process , the hard factors for the unpolarized and single - spin - dependent cross sections for these channels can be summarized by the following master formulas : @xmath162 where @xmath163 labels an individual feynman diagram ( meant as the square of an amplitude or an interference between two amplitudes , see below ) with @xmath164 the associated expression for the hard factor . @xmath165 is the color - factor for the unpolarized cross section , @xmath166 for the initial state interaction for the single - spin dependent cross section , and @xmath167 and @xmath168 for the final state interactions when the gluon attachment is to the lines of momentum @xmath62 and @xmath63 respectively .    from the above analysis , we can easily see that the hard factors in eqs . ( [ e29]),([e30 ] ) are model - independent . the model - dependence is removed once we factorize out the unpolarized quark distribution or the sivers function from the differential cross sections . from the above master formulas , the hard factors depend on the partonic scattering amplitude squared , and the color - factors , derived from partonic diagrams with or without the additional gluon attachment . both the color and partonic factors , and therefore the hard factors , are independent of the model for the proton used in our calculation .        -0.4 cm    .the color - factors and hard cross sections of individual diagrams of fig . [ f9 ] for the quark - quark scattering channels . [ cols="<,^,^,^,^,^,^",options="header " , ]     in summary , in the intermediate transverse momentum region @xmath169 , we will only have leading - power contributions to the derivative terms from the expansion of the delta function for the on - shell condition of the radiated gluon @xmath64 , and when @xmath64 is parallel to @xmath89 . we have shown these diagrams in fig . this part of the contribution is relatively easy to derive , as we mentioned above , by multiplying a factor @xmath170 to the unpolarized cross sections . the only additional involvement is the color factor . we list the color - factors for all diagrams of fig . [ f15 ] in table iv . each of the color factors has a factor @xmath171 . this can be easily seen for the color factors for diagrams @xmath172 , @xmath173 and @xmath174 , for which we have the following simplifications , @xmath175 where @xmath149 , @xmath150 and @xmath151 have been defined above in eqs . ( [ e25]),([e26]),([e27 ] ) . combining these color factors with the squared amplitude calculated in the last subsection for the unpolarized cross section , we find that the contributions to the derivative terms can be factorized into the hard factor given in eq . ( [ hqq ] ) times the derivative terms of the perturbatively generated sivers function for the polarized nucleon in eq . ( [ e38 ] ) . we illustrate this factorization in fig . for example , the derivative contribution from diagram ( @xmath172 ) is factorized into a soft pole diagram contribution to the derivative term in the sivers function where the radiated gluon attaches to the quark line shown in fig . [ f17 ] , multiplied by the color - factor @xmath149 and the partonic amplitude squared for @xmath16 . @xmath173 and @xmath174 will be factorized into the same diagram contribution to the quark sivers function , but with different color - factors @xmath150 and @xmath151 , respectively . the same conclusion holds for all other diagrams . for example , the sum of diagrams ( @xmath176-@xmath177 ) is factorized into a soft pole diagram contribution to the sivers function where the radiated gluon attaches to the gauge link , multiplied by the color factor @xmath149 and the partonic amplitude squared for @xmath16 . is parallel to @xmath89 . this part can be factorized into a soft pole contribution to the perturbatively generated sivers function in polarized hadron @xmath89 . this factorization applies both to the derivative and the non - derivative contributions from the soft poles . the circle with the attached gluon in the feynman diagrams represents all possible gluon attachments . likewise , the gluon attachment to the circle in the sivers function diagram includes the diagrams with gluon attachments to the gauge link and to the quark line.__,width=453 ]    -0.4 cm    finally , we present the total result for the derivative contribution by the @xmath178 process to the single - spin asymmetry for the dijet - correlation , @xmath179 \nonumber\\ & \times & \frac{\alpha_s}{2\pi^2}\ , \left[-\frac{1}{2n_c}\right ] \int dx\ , \xi \left[x\frac{\partial}{\partial x}t_f(x , x)\right ] \left[1+\xi^2\right ] \ , , \label{x - ssa - d}\end{aligned}\ ] ] where the superscript `` d '' indicates the derivative term , and where the single - scale partonic hard part @xmath180 has been defined in eq . ( [ hqq ] ) . the result in eq . ( [ x - ssa - d ] ) indeed reproduces the derivative terms in the first order tmd factorization formula in eq . ( [ e39 ] ) . in addition to the derivative contributions , the soft poles also contribute to the non - derivative terms . one way of calculating the non - derivative terms is to follow the @xmath181-flow in the scattering amplitude squared , summing up all contributions . keeping in mind that @xmath182 , we have to keep the full dependence on @xmath183 and @xmath184 in the calculations . it is easy to keep track of the explicit @xmath185 dependence coming from momentum conservation . in addition , we also need to include the `` indirect '' @xmath185 dependence resulting from kinematic constrains : the on - shell condition of the radiated gluon and the double poles in the final state interactions . for example , the on - shell delta function for the momentum @xmath64 of the radiated gluon requires that the longitudinal momentum component depends on @xmath181 , which will lead to a @xmath181 dependence for the longitudinal momentum fractions of @xmath186 and @xmath187 as well .    as mentioned above , if the polarized gluon attaches to the left side of the cut , the radiated gluon will have momentum @xmath188 . since @xmath187 has transverse momentum @xmath183 , the transverse component of @xmath189 will be equal to @xmath190 . furthermore , the minus component of @xmath189 is fixed by momentum conservation to @xmath191 , where the contribution by @xmath192 has been neglected , because it is of order @xmath193 . therefore , @xmath194 does not depend linearly on @xmath181 . the plus component of @xmath189 , on the other hand , will have linear @xmath181 dependence because of the on - shell condition . we therefore parameterize @xmath189 as @xmath195 where we have indicated the term linear in @xmath183 in the plus component of @xmath189 . as mentioned above , the minus component @xmath196 does not depend linearly on @xmath181 and can be determined to be @xmath197 . by using the on - shell condition @xmath198 , we can then obtain the @xmath181-dependent term in the plus component of @xmath189 , and the final parameterization result for @xmath189 is @xmath199 it is easy to check that @xmath198 up to terms quadratic in @xmath183 . from this result , we can also determine the @xmath183-dependence of the plus components of @xmath186 and @xmath187 . for example , @xmath187 can be directly calculated from momentum conservation , and we find @xmath200p_a^\mu+k_{2\perp}^\mu \ , \label{e76}\ ] ] where the minus component of @xmath187 has again been neglected . we next obtain the parameterization for @xmath186 by @xmath201 , where @xmath202 is determined by the pole and will be different for different poles . for the initial state interaction , the pole @xmath203 gives @xmath204 , and the longitudinal momentum fraction does not depend on the linear term in @xmath181 . however , for the final state interaction with gluon attachment to the line with momentum @xmath62 , the pole will be given by @xmath205 , and the plus component of @xmath202 will depend on @xmath181 through @xmath206 . similarly , the final state interaction with @xmath63 will lead to @xmath207 . in this way , the parameterizations for @xmath186 and @xmath202 will vary for different poles , whereas that for @xmath187 will always remain the same . this is because the parameterization of @xmath187 comes from the on - shell condition for the radiated gluon , which is not sensitive to the location of the poles . when the polarized gluon attaches to the right of the cut lines , on the other hand , it is @xmath186 that is determined by the on - shell condition for the radiated gluon , @xmath208 . following the same method as above , we find @xmath209p_a^\mu+k_{1\perp}^\mu \ .\end{aligned}\ ] ] furthermore , in this case @xmath210 , with @xmath202 determined from the position of the poles ; @xmath187 will be different for the different poles .    with the above parameterizations for all the external particles momenta , the calculations for the linear expansion in @xmath181 are straightforward . the final result for the soft pole contributions to the non - derivative terms when @xmath64 is nearly parallel to @xmath89 have the following form : @xmath211 \nonumber\\ & \times & \frac{\alpha_s}{2\pi^2}\ , \left[-\frac{1}{2n_c}\right ] \int dx\ , \xi \ , t_f(x , x ) \left[\frac{2\xi^3 - 3\xi^2 - 1}{1-\xi}\right]\ , , \label{x - ssa - nd - sa}\end{aligned}\ ] ] where @xmath180 is given in eq . ( [ hqq ] ) . the separate factorization for initial and final state interactions diagrams in fig . [ f17 ] , which we used to describe the factorization of the derivative contributions , also applies to the non - derivative ones . for example , the contributions from diagrams ( @xmath172-@xmath177 ) of fig . [ f15 ] can be factorized into the sivers function shown in fig . [ f17 ] , multiplied by @xmath149 and the partonic scattering function , including both the derivative and non - derivative contributions . the same happens for the final state interaction diagrams ( @xmath173-@xmath212 ) and ( @xmath174-@xmath213 ) . there are also hard pole contributions when @xmath64 is parallel to @xmath89 . we have shown these diagrams in fig .  [ f18 ] . as discussed above , the hard pole diagrams only contribute to the non - derivative terms . to calculate their contributions , we follow the above method as for the soft pole contributions . because the hard poles appear in the internal propagator , we have to take _ all _ possible gluon attachments into account . for example , diagrams @xmath172-@xmath214 of fig . [ f18 ] display all contributions with the same hard pole ( indicated by the short bars in the propagators ) , with three different attachments of the polarized gluon . another important feature for the hard pole contributions is that , after taking the pole , the polarized gluon will have nonzero longitudinal momentum fraction , i.e. , @xmath215 as we mentioned above . this is in fact the reason for referring to these poles as `` hard '' poles . for example , the hard poles in the diagrams @xmath172-@xmath214 have the momentum @xmath216 with @xmath217 , and their pole arises at the following value for @xmath218 : @xmath219 where the last equality comes from the fact that @xmath64 is parallel to @xmath89 and is dominated by its plus component as shown in eq . ( [ e66 ] ) . we have also neglected all higher power corrections in @xmath220 in the above derivations . from this result , we find the associated twist - three quark - gluon correlation function will take the form @xmath221 , that is , its two arguments are not at the same . this is different from the soft pole contributions where @xmath222 and @xmath223 are equal in the twist-3 quark - gluon correlation matrix element . is parallel to @xmath89 . this part can be factorized into a hard pole contribution to the sivers function in the polarized hadron @xmath224 . the little circle in the diagrams represents all possible gluon attachments , as required by gauge - invariance . for example , the first one in the left panel represents all three diagrams @xmath172-@xmath214 of fig . [ f18].__,width=453 ]    -0.4 cm    in our calculations , we further choose physical polarization for the radiated gluon , with the following polarization tensor : @xmath225 with this choice , we only have contributions from the three - gluon vertex diagrams @xmath176 , @xmath226 and @xmath227 and , apart from the color factor , each of them contributes a term @xmath228 . their color factors are straightforward to decompose : @xmath229 we then directly find that the contributions of these diagrams can be factorized into the hard pole contributions to the sivers function , multiplied by @xmath149 ( or @xmath150 and @xmath151 ) and the partonic hard - scattering function . we display this factorization graphically in fig . the total contribution is given by @xmath230 \nonumber\\ & \times & \frac{\alpha_s}{2\pi^2}\ , \left[-\left(\frac{1}{2n_c } + c_f\right)\right ] \int dx\ , \xi \ , t_f(x , x-\hat x_g ) \left[\frac{1+\xi}{1-\xi}\right ] \ , \label{x - ssa - nd - ha}\end{aligned}\ ] ] where @xmath231 and @xmath180 is given in eq . ( [ hqq ] ) . indeed , the above result reproduces the contribution to eq . ( [ e39 ] ) that arises from the hard pole part of the quark sivers function .    in summary , by combining all the soft - pole and hard - pole contributions when the radiated gluon is nearly parallel to the polarized nucleon , we can reproduce the contribution from the perturbatively generated sivers function to the ssa in dijet - correlations in eq . ( [ e39 ] ) , derived from the tmd factorization formula in the regime @xmath169 . with these explicit calculations , we have demonstrated that at this perturbative order the contribution to the ssa when @xmath64 is parallel to @xmath89 can be factorized into the quark sivers function , multiplied by the initial / final state interaction color - factors ( @xmath149 , @xmath150 , and @xmath151 ) , and the partonic scattering amplitude squared . the calculation for the contributions when @xmath64 is parallel to @xmath113 follows the same procedure as above . as we have shown earlier , there is no contribution to the derivative terms when @xmath64 is parallel to @xmath113 , and so we only need to consider the non - derivative contributions . we will also have both soft and hard poles contributions . in fig . [ f20 ] , we have shown the soft pole diagrams , and in figs . [ f21],[f22 ] the hard pole ones . again we have to keep track of the @xmath181-flow in the scattering amplitudes . first , we need to parameterize the momenta of all the external particles including their linear dependence on @xmath181 . for example , when the polarized gluon attaches to the left of the cut line in the diagrams of figs .  [ f20]-[f22 ] , the momentum @xmath189 of the radiated gluon and the quark momentum @xmath187 can be parameterized as @xmath232 and @xmath233 with @xmath202 determined by the pole . note that these parameterizations are different from those in eqs . ( [ e75]),([e76 ] ) . this is because , when @xmath64 is parallel to @xmath113 , its minus component is dominant and is parameterized by @xmath234 . again , this minus component does not depend linearly on @xmath181 . the plus component of @xmath189 does have linear dependence , which can be calculated from the on - shell condition @xmath198 . the linear @xmath181-dependence for the plus component of @xmath187 can be calculated as well .    likewise , when the gluon attaches to the right of the cut , we will have @xmath235 and @xmath210 with @xmath202 determined from the position of the poles . again , different poles will give different results for @xmath202 , following the above discussions for the soft pole contributions . is parallel to @xmath113 . this contribution is factorized into the tmd quark distribution in the unpolarized nucleon @xmath236 . the little circle in the diagrams represents all possible gluon attachments . the bars indicate the poles . the second and third diagrams in the left panel represent the soft pole diagrams @xmath173-@xmath212 and @xmath174-@xmath213 of fig . [ f20 ] respectively . the first one contains both soft poles ( represented by the bar on the circle ( diagrams @xmath172-@xmath177 of fig . [ f20 ] ) ) and hard poles ( indicated by the bar on the straight line ( all diagrams in fig .  [ f21])).__,width=453 ]    -0.4 cm    by using the above parameterizations for the momenta of the external particles , we can calculate the contributions from the soft pole diagrams in fig . [ f20 ] , and the hard pole diagrams in figs . [ f21 ] and [ f22 ] . we further find that the contribution from the hard pole diagrams in fig . [ f22 ] is power suppressed by @xmath15 . so , for the hard pole contributions , we only need to consider the diagrams shown in fig . [ f21 ] , which all have the same hard pole indicated by the bars on the propagators . as discussed in the last subsection , the hard pole contributions may in general lead to the polarized gluon having nonzero longitudinal momentum fraction , @xmath215 . however , when the radiated gluon s momentum @xmath64 is nearly parallel to @xmath113 , the hard poles lead to @xmath237 at leading power , as was the case for the soft pole contributions . this is very different from what we found for the hard pole contributions in the last subsection in the case @xmath64 parallel to @xmath89 . the momentum giving the hard pole , @xmath238 with @xmath217 , is the same as in the last subsection . however , taking the pole , we now find @xmath239 where the last equation holds in the limit @xmath240 . the reason for this is that , when @xmath64 becomes parallel to @xmath113 , it will be dominated by its minus component and can be parameterized as @xmath241 . substituting this into the above equation  ( [ xdef ] ) , we find that the numerator is of order @xmath242 , whereas the denominator is of order @xmath118 . so , in the leading power approximation , @xmath243 for the hard pole contributions when @xmath64 is nearly parallel to @xmath113 . the contributions from the soft poles in fig . [ f20 ] and the hard poles in fig . [ f21 ] is added to obtain the final results . we emphasize again that our power - counting based analysis is only valid for the leading power terms in the @xmath15 expansion .    the final result for the contributions from figs .  [ f20 ] and [ f21 ] , plus their mirror diagrams , is @xmath244 \nonumber\\ & \times & \frac{\alpha_s}{2\pi^2}\ , c_f \ , \int dx'\ , \xi ' \ , q'(x ' ) \left[\frac{1+\xi^{\prime 2}}{1-\xi'}\right ] \ , \label{x - ssa - nd - b}\end{aligned}\ ] ] which is again proportional to the same hard factor @xmath245 , and reproduces the last term in the tmd factorization formula eq . ( [ e39 ] ) . it is important to note that a part of the hard pole contributions from fig . [ f21 ] cancels out the soft pole contributions from diagrams @xmath172-@xmath177 of fig . the remainder of the hard pole contribution is factorized into the quark distribution multiplied by @xmath149 and the partonic scattering function . the soft pole diagrams @xmath173-@xmath212 and @xmath174-@xmath213 are separately factorized into the same quark distribution for the unpolarized nucleon multiplied by @xmath150 and @xmath151 , respectively , and the partonic scattering function . we illustrate this factorization in fig . [ f23 ] .    after summing up all contributions from eqs . ( [ x - ssa - d ] ) , ( [ x - ssa - nd - sa ] ) , ( [ x - ssa - nd - ha ] ) , and ( [ x - ssa - nd - b ] ) , we obtain the total leading - power contribution to @xmath13 from the @xmath246 partonic subprocess in the @xmath15 expansion . it reproduces the factorized formula in eq . ( [ e39 ] ) , which is the leading term in the corresponding tmd factorization formula , eq . ( [ e4 ] ) .    in summary , we have calculated the single transverse - spin asymmetry in dijet - correlation in the twist - three approach when the radiated gluon is parallel to either the polarized nucleon or the unpolarized nucleon . by using the power counting technique , we have shown that the contribution to the ssa can be factorized into the perturbatively generated quark sivers function when the radiated gluon is parallel to the polarized nucleon , and into the unpolarized quark distribution when it is parallel to the unpolarized nucleon . we have demonstrated that the result reproduces the leading order terms in the corresponding tmd factorization formalism in the regime @xmath247 . this is a nontrivial result , especially for the single spin asymmetry , because the calculations for the quark sivers function at this order have to take into account the perturbative expansion of the gauge link at next - to - leading order ( @xmath248 ) . the single spin dependent differential cross section term in eq . ( [ e4 ] ) can be further simplified by taking a moment in @xmath9 . such a moment was also considered in @xcite . we can rewrite @xmath249 where @xmath250 is the so - called bi - sector angle of the two jets : @xmath251 with @xmath252 and @xmath253 the azimuthal angles of the two jets relative to the polarization vector @xmath34 . the angle @xmath254 measures how far the two jets are away from the back - to - back configuration . all these azimuthal angles are defined in a frame in which the polarized proton is moving in the @xmath255 direction .    since @xmath256 , the @xmath9-moment of the asymmetry is also related to the @xmath257-moment . one interesting moment is the following : @xmath258 where @xmath259 is the twist-3 matrix element of the quark - gluon correlation function defined in eq . ( [ tf ] ) @xcite . it is also related to the @xmath260-moment of the tmd quark sivers function , see eq . ( [ e34 ] ) . an advantage of taking this moment of the asymmetry is that the transverse momentum integrals of the various factors in the factorization formula decouple from each other , without further assumptions for the @xmath260-dependence of the tmd distributions . in fact , the ssa in the @xmath9 moment of the dijet momentum imbalance is effectively a physical quantity with only a _ single _ large observed scale , @xmath8 . for such an observable , the conventional collinear factorization is more appropriate . in this case , all partons transverse momentum dependence is integrated into the usual collinear parton distributions . the hard factors we obtained above for all the partonic channels agree with those given in @xcite , although a very different approach was adopted there . in @xcite , in order to investigate the @xmath9-moment of the single spin asymmetry , the appropriate gauge links were derived and expanded to first order in @xmath261 . since both approaches have included the initial / final state interactions in their formalisms , they should agree with each other at this order . in this paper , we have studied the asymmetric production of two jets in hadronic collisions . using the collinear factorization approach in perturbative qcd , we calculated both the spin - averaged and the spin - dependent differential cross section for the hadronic production of dijets with momentum imbalance @xmath9 , in a kinematic region where @xmath262 . at the first nonvanishing order , the momentum imbalance is generated by radiating a gluon with transverse momentum equal to the imbalance . in the limit when the imbalance @xmath9 is much less than the averaged jet momentum @xmath8 , we derived the leading - power contribution to the cross sections in the @xmath15 expansion when the radiated gluon is nearly parallel to one of the incoming hadrons . we found that the perturbatively calculated leading contributions to both the spin - averaged and the spin - dependent cross sections can be factorized into a partonic hard part which is a function of @xmath8 , and into perturbatively generated tmd parton distributions at @xmath248 . our results derived in the collinear factorization approach in the limit @xmath263 reproduce the same factorized expressions in eqs . ( [ e37]),([e39 ] ) that were derived as the leading order terms of the generalized tmd factorization formulas in eqs . ( [ e4]),([e2 ] ) when the observed @xmath9 is solely due to the perturbatively generated tmd parton distribution of one of the incoming hadrons . the consistency between the results derived in the collinear factorization approach and those derived by expanding the generalized tmd factorization formulas in eqs .  ( [ e4]),([e2 ] ) to the same order naturally leads us to ask if the generalized tmd factorization formulas in eqs . ( [ e4]),([e2 ] ) are actually valid for the hadronic dijet production to all orders in perturbative qcd when @xmath264 , like the tmd factorization formulas for drell - yan and sidis . _ if this were true _ , the consistency that we demonstrated in this paper in the overlap region where @xmath50 would effectively provide a unified picture and treatment for the hadronic dijet production over the full kinematic range of jet momenta : using the tmd factorization formalism when the dijet momentum imbalance is small and the collinear factorization approach when it is large . this is exactly what was achieved in ref . @xcite for the drell - yan and sidis processes . however , as we mentioned at the beginning of sec . iii , showing the consistency between two factorized formalisms at the leading order when @xmath265 is not sufficient to prove the tmd factorization formula over its full kinematic region . when we started with the collinear factorization approach to calculate the momentum imbalance of the dijet cross section , we already took as a fact that the long - distance interaction between the hadrons / jets at their mass scales , @xmath266 , are factorized for the given @xmath267 . that is , our results derived in the limit @xmath265 are valid when @xmath268 ; but , our derivation does not address the issues concerning the factorization of the interactions between the hadrons / jets when @xmath269 or the mass scale of the hadrons / jets involved . the consistency that we demonstrated here is certainly a necessary condition that needs to be satisfied should the tmd factorization formulas be valid .    for the tmd factorization to be true , we need to establish two key facts : the partonic hard part at @xmath54 must be insensitive to any physics at the scale @xmath78 , and the tmd parton distributions at the scale @xmath9 should be process independent and universal , or at least , `` quasi''-universal up to a sign change @xcite . in this paper , we also calculated in sec .  ii the partonic hard parts at @xmath57 by using the brodsky - hwang - schmidt model with the proper color factor for the interaction vertices . we derived the partonic hard parts from the lowest order scattering diagrams and showed that the dynamics of partonic scattering at scale @xmath57 is independent of the model of nucleon that we used in our derivation . since the ssa can be generated by either initial- or final - state interactions , the color flow of the partonic scattering in the dijet momentum imbalance is very different from that of drell - yan and sidis . since all poles from the initial- and final - state one gluon interactions contribute to the ssa , and the color factor for each individual scattering diagram is insensitive to the long - distance details of the nucleon because there is only one unique color structure for the active quark - gluon combination at this order , we were able to factor the process - dependent as well as the diagram - specific color factors along with the partonic hard parts , and leave the physics at scale @xmath270 in the tmd parton distributions defined in sidis . by doing that , we found that the resulting full partonic hard parts are the same as those derived in the collinear factorization approach . however , for a full tmd factorization , we will have to show that the poles from multiple gluon interactions can be exponentiated into the gauge link that defines the tmd parton distributions , while leaving the same partonic hard parts with the same color factors . due to the mixture of initial- and final - state interactions and the complications of color flows , this task is very non - trivial and is beyond the scope of this paper .    to better understand the dijet momentum imbalance in hadronic collisions , a number of extensions can be performed based on our results . first , in this paper , we have only studied the contribution of one gluon radiation collinear to one of the incident hadrons . this should be extended to other kinematic regions of the gluon . for example , when the gluon is radiated parallel to the final state jets or hadrons , it is important to see how the corresponding contribution can be accommodated by the jet definition , or be factorized into the fragmentation functions for the hadron , especially in the spin - dependent case . another important kinematic regime is when the radiated gluon becomes soft . it is very important to see if the resulting contribution can be factorized into a soft factor , as predicted by the tmd factorization formalism . in addition , in this paper we have only discussed the ssa for dijet correlations , which has sivers - type contributions . for the related di - hadron correlations , one expects also the collins mechanism to be important , and it will be interesting to explore it and the factorization properties of the corresponding spin - dependent cross section . the method and procedure used in this paper should provide useful guidance in all these related studies .    we finally note that if the tmd factorization indeed fails , a careful comparison between our results and experimental data as the momentum imbalance @xmath15 decreases could furnish an ideal test of the breaking of factorization and provide new opportunities to explore qcd dynamics . we are grateful to cedran bomhof , john collins , andreas metz , and piet mulders for useful discussions . we thank xiangdong ji for his general remarks and his reminding us the wigner - eckart theorem for the color - factor decomposition . we especially thank george sterman for continued discussions on the factorization issues in hadronic reactions . is supported in part by the u. s. department of energy under grant no . de - fg02 - 87er-40371 . w.v . are finally grateful to riken , brookhaven national laboratory and the u.s . department of energy ( contract number de - ac02 - 98ch10886 ) for providing the facilities essential for the completion of their work . thanks the high energy theory group at argonne national laboratory for its hospitality during the writing of this work .      in this appendix , we list the hard factors for the dijet - correlation at hadron colliders , for the unpolarized and the single - transverse spin dependent cases . for unpolarized scattering , the hard factors are well - known in the literature . they may also be obtained from the results shown in tables i - iii , using eq . ( [ e29 ] ) . for completeness , we list them here : @xmath271-\frac{n_c^2 - 1}{4n_c^3}\frac{4\hat s^2}{\hat t\hat u}\right\ } \ , \nonumber\\ h_{q\bar q\to q\bar q}^{uu}(\hat s,\hat t,\hat u)&=&\frac{\alpha_s^2\pi}{\hat s^2}\left\{\frac{n_c^2 - 1}{4n_c^2}\left[\frac{2(\hat s^2+\hat u^2)}{\hat t^2}+\frac{2(\hat u^2+\hat t^2)}{\hat s^2}\right]-\frac{n_c^2 - 1}{4n_c^3}\frac{4\hat u^2}{\hat s\hat t}\right\ } \ , \nonumber\\ h_{qg\to qg}^{uu}(\hat s,\hat t,\hat u)&=&\frac{\alpha_s^2\pi}{\hat s^2}\left\{\frac{1}{2}\frac{2(\hat s^2+\hat u^2)}{\hat t^2}-\frac{c_f}{2n_c}\frac{2(\hat s^2+\hat u^2)}{\hat s\hat t}\right\ } \ , \nonumber\\ h_{q\bar q\to gg}^{uu}(\hat s,\hat t,\hat u)&=&\frac{\alpha_s^2\pi}{\hat s^2}\left\{\frac{c_f^2}{n_c}\frac{2(\hat t^2+\hat u^2)}{\hat t\hat u}-c_f\frac{2(\hat t^2+\hat u^2)}{\hat s^2}\right\ } \ , \nonumber\end{aligned}\ ] ] @xmath272 for the single - transverse spin dependent case , using eq.([e30 ] ) , we obtain @xmath273+\frac{n_c^2 + 3}{4n_c^3}\frac{4\hat s^2}{\hat t\hat u}\right\ } \ , \nonumber\\ h_{q\bar q\to q\bar q}^{\rm sivers}(\hat s,\hat t,\hat u)&=&\frac{\alpha_s^2\pi}{\hat s^2}\left\{-\frac{n_c^2 - 3}{4n_c^2}\frac{2(\hat s^2+\hat u^2)}{\hat t^2}+\frac{n_c^2 + 1}{4n_c^2}\frac{2(\hat u^2+\hat t^2)}{\hat s^2}-\frac{n_c^2 + 1}{4n_c^3}\frac{4\hat u^2}{\hat s\hat t}\right\ } \ , \nonumber\\ h_{qg\to qg}^{\rm sivers}(\hat s,\hat t,\hat u)&=&\frac{\alpha_s^2\pi}{\hat s^2}\left\{-\frac{n_c^2}{4(n_c^2 - 1)}\frac{2(\hat s^2+\hat u^2)}{\hat t^2}\left[\frac{\hat s}{\hat u}-\frac{\hat u}{\hat s}\right]-\frac{1}{2(n_c^2 - 1)}\frac{2(\hat s^2+\hat u^2)}{\hat t^2}\right.\nonumber\\ & & ~~~~~~~\left.-\frac{1}{4n_c^2(n_c^2 - 1)}\frac{2(\hat s^2+\hat u^2)}{\hat s\hat u}\right\ } \ , \nonumber\\ h_{q\bar q\to gg}^{\rm sivers}(\hat s,\hat t,\hat u)&=&\frac{\alpha_s^2\pi}{\hat s^2}\left\{-\frac{1}{2n_c}\frac{2(\hat t^2+\hat u^2)}{\hat s^2}+\frac{n_c}{4}\frac{2(\hat t^2+\hat u^2)}{\hat s^2}\left[\frac{\hat t}{\hat u}+\frac{\hat u}{\hat t}\right]\right.\nonumber\\ & & ~~~~~~~\left.-\frac{2n_c^2 + 1}{4n_c^3}\frac{2(\hat t^2+\hat u^2)}{\hat t\hat u}\right\ } \ , \nonumber\\ h_{qg\to q\gamma}^{\rm sivers}(\hat s,\hat t,\hat u)&=&\frac{\alpha_s\alpha_e e_q^2\pi}{\hat s^2}\left(-\frac{n_c^2 + 1}{2n_c(n_c^2 - 1)}\right)\frac{2(\hat s^2+\hat u^2)}{-\hat s\hat u } \ , \nonumber\\ h_{q\bar q\to g\gamma}^{\rm sivers}(\hat s,\hat t,\hat u)&=&\frac{\alpha_s\alpha_e e_q^2\pi}{\hat s^2}\left(\frac{n_c^2 + 1}{2n_c^2}\right)\frac{2(\hat t^2+\hat u^2)}{\hat t\hat u } \ .\end{aligned}\ ] ] the above results are for @xmath98-channel scattering . in dijet - production , both @xmath98 and @xmath99 channel processes will contribute . the hard factors for the @xmath99-channel scattering can be obtained from the results above by exchanging @xmath274 . g.  bunce _ et al . _ , lett .   * 36 * , 1113 ( 1976 ) . a.  airapetian _ et al . _ [ hermes collaboration ] , phys . lett .   * 84 * , 4047 ( 2000 ) ; phys . lett .   * 94 * , 012002 ( 2005 ) . a.  bravar [ spin muon collaboration ] , nucl . a * 666 * , 314 ( 2000 ) ; h.  avakian [ clas collaboration ] , talk presented at the rbrc workshop `` single - spin asymmetries '' , brookhaven national laboratory , upton , new york , june 1 - 3 , 2005 , rbrc proceedings volume 75 ( bnl-74717 - 2005 ) ; v.  y.  alexakhin _ et al . _ [ compass collaboration ] , phys . lett .   * 94 * , 202002 ( 2005 ) ; e.  s.  ageev _ et al . _ [ compass collaboration ] , arxiv : hep - ex/0610068 ; f.  bradamante [ compass collaboration ] , talk presented at the `` 17th international spin physics symposium ( spin 2006 ) '' , kyoto , japan , october 2 - 7 , 2006 . j.  adams _ et al . _ [ star collaboration ] , phys .  rev .  lett .   * 92 * , 171801 ( 2004 ) ; l.  nogach , talk presented at the `` 17th international spin physics symposium ( spin 2006 ) '' , kyoto , japan , october 2 - 7 , 2006 , arxiv : hep - ex/0612030 . s.  s.  adler [ phenix collaboration ] , phys . * 95 * , 202001 ( 2005 ) ; m.  chiu , talk presented at the `` 17th international spin physics symposium ( spin 2006 ) '' , kyoto , japan , october 2 - 7 , 2006 ; c.  a.  aidala , ph.d . thesis , columbia u. , arxiv : hep - ex/0601009 . f.  videbaek [ brahms collaboration ] , aip conf . * 792 * , 993 ( 2005 ) ; j.  h.  lee [ brahms collaboration ] , talk presented at the `` 17th international spin physics symposium ( spin 2006 ) '' , kyoto , japan , october 2 - 7 , 2006 . g.  l.  kane , j.  pumplin and w.  repko , phys . lett .   * 41 * , 1689 ( 1978 ) . for reviews , see : m.  anselmino , a.  efremov and e.  leader , phys . * 261 * , 1 ( 1995 ) [ erratum - ibid .   * 281 * , 399 ( 1997 ) ] ; z.  t.  liang and c.  boros , int . j.  mod . a * 15 * , 927 ( 2000 ) ; v.  barone , a.  drago and p.  g.  ratcliffe , phys . * 359 * , 1 ( 2002 ) . s.  j.  brodsky , d.  s.  hwang and i.  schmidt , phys . b * 530 * , 99 ( 2002 ) ; nucl . b * 642 * , 344 ( 2002 ) . j.  c.  collins , phys . b * 536 * , 43 ( 2002 ) . d.  w.  sivers , phys . d * 41 * , 83 ( 1990 ) ; phys . d * 43 * , 261 ( 1991 ) . x.  ji and f.  yuan , phys . b * 543 * , 66 ( 2002 ) ; a.  v.  belitsky , x.  ji and f.  yuan , nucl .  phys . b * 656 * , 165 ( 2003 ) . d.  boer , p.  j.  mulders and f.  pijlman , nucl .  phys . b * 667 * , 201 ( 2003 ) . j.  c.  collins and d.  e.  soper , nucl . b * 193 * , 381 ( 1981 ) [ erratum - ibid . b * 213 * , 545 ( 1983 ) ] ; nucl . b * 197 * , 446 ( 1982 ) . j.  c.  collins , d.  e.  soper and g.  sterman , nucl . b * 250 * , 199 ( 1985 ) . d.  boer and w.  vogelsang , phys .  rev . d * 69 * , 094025 ( 2004 ) . b.  i.  abelev , _ et al . _ , [ star collaboration ] arxiv:0705.4629 [ hep - ex ] ; j.  balewski , talk presented at the spin 2006 symposium , kyoto , japan , october 2 - 7 , 2006 , arxiv : hep - ex/0612036 . c.  j.  bomhof , p.  j.  mulders and f.  pijlman , phys . b * 596 * , 277 ( 2004 ) ; eur . j.  c * 47 * , 147 ( 2006 ) ; a.  bacchetta , c.  j.  bomhof , p.  j.  mulders and f.  pijlman , phys . d * 72 * , 034030 ( 2005 ) ; c.  j.  bomhof and p.  j.  mulders , jhep * 0702 * , 029 ( 2007 ) . j.  w.  qiu , w.  vogelsang and f.  yuan , arxiv:0704.1153 [ hep - ph ] , phys . b in press . a.  v.  efremov and o.  v.  teryaev , sov . j.  nucl . * 36 * , 140 ( 1982 ) [ yad . 36 * , 242 ( 1982 ) ] ; a.  v.  efremov and o.  v.  teryaev , phys . b * 150 * , 383 ( 1985 ) . j.  qiu and g.  sterman , phys . lett .   * 67 * , 2264 ( 1991 ) ; nucl . b * 378 * , 52 ( 1992 ) ; phys . d * 59 * , 014004 ( 1999 ) . j.  botts and g.  sterman , nucl . b * 325 * , 62 ( 1989 ) . n.  kidonakis and g.  sterman , nucl . b * 505 * , 321 ( 1997 ) ; nucl . b * 525 * , 299 ( 1998 ) ; n.  kidonakis , g.  oderda and g.  sterman , nucl . b * 531 * , 365 ( 1998 ) . x.  ji , j.  w.  qiu , w.  vogelsang and f.  yuan , phys .  rev .  lett . * 97 * , 082002 ( 2006 ) ; phys . d * 73 * , 094017 ( 2006 ) ; phys . b * 638 * , 178 ( 2006 ) . w.  vogelsang and f.  yuan , phys . d * 72 * , 054028 ( 2005 ) . j.  collins and j.  w.  qiu , arxiv:0705.2141 [ hep - ph ] . g.  altarelli and g.  parisi , nucl . b * 126 * , 298 ( 1977 ) . g.  t.  bodwin , phys . d * 31 * , 2616 ( 1985 ) [ erratum - ibid . d * 34 * , 3932 ( 1986 ) ] . x.  ji , phys . b * 289 * , 137 ( 1992 ) . a.  bacchetta , c.  bomhof , u.  dalesio , p.  j.  mulders and f.  murgia , arxiv : hep - ph/0703153 .    c.  j.  bomhof , p. j.  mulders , w.  vogelsang and f.  yuan , phys .  rev . d * 75 * , 074019 ( 2007 ) . j.  w.  qiu and g.  sterman , aip conf . proc .   * 223 * , 249 ( 1991 ) ; nucl . b * 353 * , 137 ( 1991 ) . c.  kouvaris , j.  w.  qiu , w.  vogelsang and f.  yuan , phys . d * 74 * , 114013 ( 2006 ) . h.  eguchi , y.  koike and k.  tanaka , nucl . b * 752 * , 1 ( 2006 ) ; nucl . b * 763 * , 198 ( 2007 ) ; y.  koike and k.  tanaka , phys . b * 646 * , 232 ( 2007 ) ; arxiv : hep - ph/0703169 . g.  ratcliffe and o.  v.  teryaev , arxiv : hep - ph/0703293 .
Abstract:

Answer: we study the single transverse - spin asymmetry for dijet production in hadronic collisions in both the collinear qcd factorization approach and the brodsky - hwang - schmidt model . we show that a nonvanishing asymmetry is generated by both initial - state and final - state interactions , and that the final - state interactions dominate . we find that in the leading kinematic region where the transverse momentum imbalance of the two jets , @xmath0 , is much less than the momentum of either jet , the contribution from the lowest non - trivial perturbative order to both the spin - averaged and the spin - dependent dijet cross sections can be factorized into a hard part that is a function only of the averaged jet momentum @xmath1 , and perturbatively generated transverse momentum dependent ( tmd ) parton distributions . we show that the spin asymmetry at this non - trivial perturbative order can be described by the tmd parton distributions defined in either semi - inclusive dis or the drell - yan process . we derive the same hard parts from both the collinear factorization approach and in the context of the brodsky - hwang - schmidt model , verifying that they are not sensitive to details of the factorized long distance physics .