• We present a new analysis of $A_N$ in $p^\uparrow p\to \pi\, X$ within the collinear twist-3 factorization formalism. We incorporate recently derived Lorentz invariance relations into our calculation and focus on input from the kinematical twist-3 functions, which are weighted integrals of transverse momentum dependent (TMD) functions. In particular, we use the latest extractions of the Sivers and Collins functions with TMD evolution to compute certain terms in $A_N$. Consequently, we are able to constrain the remaining contributions from the lesser known dynamical twist-3 correlators.
  • We study jet mass distributions measured in the single inclusive jet production in proton-proton collisions $pp\to \text{jet}+X$ at the LHC. We consider both standard ungroomed jets as well as soft drop groomed jets. Within the Soft Collinear Effective Theory (SCET), we establish QCD factorization theorems for both cases and we study their relation. The developed framework allows for the joint resummation of several classes of logarithmic corrections to all orders in the strong coupling constant. For the ungroomed case we resum logarithms in the jet radius parameter and in the small jet mass. For the groomed case we resum in addition the logarithms in the soft threshold parameter $z_{\text{cut}}$ which is introduced by the soft drop grooming algorithm. In this way we are able to reliably determine the absolute normalization of the groomed jet mass distribution in proton-proton collisions. All logarithmic corrections are resummed to the next-to-leading logarithmic accuracy. We present numerical results and compare with the available data from the LHC. For both the groomed and ungroomed jet mass distributions we find very good agreement after including non-perturbative corrections.
  • We study jet angularity measurements for single-inclusive jet production at the LHC. Jet angularities depend on a continuous parameter $a$ allowing for a smooth interpolation between different traditional jet shape observables. We establish a factorization theorem within Soft Collinear Effective Theory (SCET) where we consistently take into account in- and out-of-jet radiation by making use of semi-inclusive jet functions. For comparison, we elaborate on the differences to jet angularities measured on an exclusive jet sample. All the necessary ingredients for the resummation at next-to-leading logarithmic (NLL) accuracy are presented within the effective field theory framework. We expect semi-inclusive jet angularity measurements to be feasible at the LHC and we present theoretical predictions for the relevant kinematic range. In addition, we investigate the potential impact of jet angularities for quark-gluon discrimination.
  • We study the production and polarization of $J/\psi$ mesons within a jet in proton-proton collisions at the LHC. We define the $J/\psi$-jet fragmentation function as a ratio of differential jet cross sections with and without the reconstructed $J/\psi$ in the jet. We demonstrate that this is a very useful observable to help explore the $J/\psi$ production mechanism, and to differentiate between different NRQCD global fits based on inclusive $J/\psi$ cross sections. Furthermore, we propose to measure the polarization of $J/\psi$ mesons inside the jet, which can provide even more stringent constraints for the heavy quarkonium production mechanism.
  • We present a framework that describes the energy distribution of subjets of radius $r$ within a jet of radius $R$. We consider both an inclusive sample of subjets as well as subjets centered around a predetermined axis, from which the jet shape can be obtained. For $r \ll R$ we factorize the physics at angular scales $r$ and $R$ to resum the logarithms of $r/R$. For central subjets, we consider both the standard jet axis and the winner-take-all axis, which involve double and single logarithms of $r/R$, respectively. All relevant one-loop matching coefficients are given, and an inconsistency in some previous results for cone jets is resolved. Our results for the standard jet shape differ from previous calculations at next-to-leading logarithmic order, because we account for the recoil of the standard jet axis due to soft radiation. Numerical results are presented for an inclusive subjet sample for $pp \to {\rm jet}+X$ at next-to-leading order plus leading logarithmic order.
  • We investigate the Collins azimuthal asymmetry of hadrons produced inside jets in transversely polarized proton-proton collisions. Recently, the quark transversity distributions and the Collins fragmentation functions have been extracted within global analyses from data of the processes semi-inclusive deep inelastic scattering and electron-positron annihilation. We calculate the Collins azimuthal asymmetry for charged pions inside jets using these extractions for RHIC kinematics at center-of-mass energies of 200 and 500 GeV. We compare our results with recent data from the STAR Collaboration at RHIC and find good agreement, which confirms the universality of the Collins fragmentation functions. In addition, we further explore the impact of transverse momentum dependent evolution effects.
  • We study the transverse momentum distribution of hadrons within jets, where the transverse momentum is defined with respect to the standard jet axis. We consider the case where the jet substructure measurement is performed for an inclusive jet sample $pp\to\text{jet}+X$. We demonstrate that this observable provides new opportunities to study transverse momentum dependent fragmentation functions (TMDFFs) which are currently poorly constrained from data, especially for gluons. The factorization of the cross section is obtained within Soft Collinear Effective Theory (SCET), and we show that the relevant TMDFFs are the same as for the more traditional processes semi-inclusive deep inelastic scattering (SIDIS) and electron-positron annihilation. Different than in SIDIS, the observable for the in-jet fragmentation does not depend on TMD parton distribution functions which allows for a cleaner and more direct probe of TMDFFs. We present numerical results and compare to available data from the LHC.
  • Vector boson-tagged jet production in collisions of heavy nuclei opens new opportunities to study parton shower formation and propagation in strongly interacting matter. It has been argued to provide a golden channel that can constrain the energy loss of jets in the quark-gluon plasma created in heavy ion reactions. We present theoretical results for isolated photon-tagged and $Z^0$ boson-tagged jet production in Pb+Pb collisions with $\sqrt{s_{NN}} = 5.02$ TeV at the LHC. Specifically, we evaluate the transverse momentum imbalance ${\rm x_{JV}}$ distribution and nuclear modification factor ${\rm I_{AA}}$ of tagged jets and compare our theoretical calculations to recent experimental measurements by ATLAS and CMS collaborations. Our analysis, which includes both collisional and radiative energy losses, sheds light on their relative importance versus the strength of jet-medium interactions and helps quantify the amount of out-of-cone radiation of predominantly prompt quark-initiated jets.
  • Future experiments at the Jefferson Lab 12 GeV upgrade, in particular, the Solenoidal Large Intensity Device (SoLID), aim at a very precise data set in the region where the partonic structure of the nucleon is dominated by the valence quarks. One of the main goals is to constrain the quark transversity distributions. We apply recent theoretical advances of the global QCD extraction of the transversity distributions to study the impact of future experimental data from the SoLID experiments. Especially, we develop a simple strategy based on the Hessian matrix analysis that allows one to estimate the uncertainties of the transversity quark distributions and their tensor charges extracted from SoLID data simulation. We find that the SoLID measurements with the proton and the effective neutron targets can improve the precision of the u- and d-quark transversity distributions up to one order of magnitude in the range 0.05 < x < 0.6.
  • We develop a new formalism to describe the inclusive production of small radius jets in heavy-ion collisions, which is consistent with jet calculations in the simpler proton-proton system. Only at next-to-leading order (NLO) and beyond, the jet radius parameter $R$ and the jet algorithm dependence of the jet cross section can be studied and a meaningful comparison to experimental measurements is possible. We are able to consistently achieve NLO accuracy by making use of the recently developed semi-inclusive jet functions within Soft Collinear Effective Theory (SCET). In addition, single logarithms of the jet size parameter $\alpha_s^n\ln^n R$ are resummed to next-to-leading logarithmic (NLL$_R$) accuracy. The medium modified semi-inclusive jet functions are obtained within the framework of SCET with Glauber gluons that describe the interaction of jets with the medium. We present numerical results for the suppression of inclusive jet cross sections in heavy ion collisions at the LHC and the formalism developed here can be extended directly to corresponding jet substructure observables.
  • Within the framework of higher-twist collinear factorization, transverse momentum broadening for the final hadrons in semi-inclusive deeply inelastic $e+A$ collisions is studied at the next-to-leading order (NLO) in perturbative QCD. Through explicit calculations of real and virtual corrections at twist-4, the transverse-momentum-weighted differential cross section due to double scattering is shown to factorize at NLO and can be expressed as a convolution of twist-4 nuclear parton correlation functions, the usual twist-2 fragmentation functions and hard parts which are finite and free of any divergences. A QCD evolution equation is also derived for the renormalized twist-4 quark-gluon correlation function which can be applied to future phenomenological studies of transverse momentum broadening and jet quenching at NLO.
  • We review recent progress in the theoretical description of hard probes in heavy-ion collisions within the framework of Soft Collinear Effective Theory (SCET). Firstly, we consider the inclusive production of heavy flavor mesons and jets with a small radius parameter in proton-proton collisions. Secondly, we describe how in-medium effects can be incorporated consistently with the proton-proton baseline calculations. The in-medium interactions are described by Glauber gluon interactions. We present results for the nuclear modification factor $R_{AA}$ which is most commonly used to study the quenching of hadron or jet production yields in heavy-ion collisions and we compare to recent data from the LHC.
  • We review a new approach to open heavy flavor production in heavy ion collisions based on Soft Collinear Effective Theory (SCET). We include both finite heavy quark masses in the SCET Lagrangian and Glauber gluons that describe the interaction of collinear partons with the hot and dense QCD medium. From the new effective field theory, we derive massive in-medium splitting kernels and propose a new framework for including in-medium interactions consistent with next-to-leading order calculations in QCD. We present numerical results for the suppression of $D$- and $B$-mesons and compare to results obtained within the traditional approach to parton energy loss. We find good agreement when comparing to existing data from the LHC at $\sqrt{s_{\mathrm{NN}}}=5.02$ TeV and 2.76 TeV.
  • We propose a new method to evaluate jet substructure observables in inclusive jet measurements, based upon semi-inclusive jet functions in the framework of Soft Collinear Effective Theory (SCET). As a first example, we consider the jet fragmentation function, where a hadron $h$ is identified inside a fully reconstructed jet. We introduce a new semi-inclusive fragmenting jet function ${\mathcal G}^h_i(z= \omega_J/\omega,z_h=\omega_h/\omega_J,\omega_J, R,\mu)$, which depends on the jet radius $R$ and the large light-cone momenta of the parton `$i$' initiating the jet ($\omega$), the jet ($\omega_J$), and the hadron $h$ ($\omega_h$). The jet fragmentation function can then be expressed as a semi-inclusive observable, in the spirit of actual experimental measurements, rather than as an exclusive one. We demonstrate the consistency of the effective field theory treatment and standard perturbative QCD calculations of this observable at next-to-leading order (NLO). The renormalization group (RG) equation for the semi-inclusive fragmenting jet function ${\mathcal G}_i^h(z,z_h, \omega_J, R,\mu)$ are also derived and shown to follow exactly the usual timelike DGLAP evolution equations for fragmentation functions. The newly obtained RG equations can be used to perform the resummation of single logarithms of the jet radius parameter $R$ up to next-to-leading logarithmic (NLL$_R$) accuracy. In combination with the fixed NLO calculation, we obtain NLO+NLL$_R$ results for the hadron distribution inside the jet. We present numerical results for $pp\to(\mathrm{jet}\,h)X$ in the new framework, and find excellent agreement with existing LHC experimental data.
  • We introduce a new kind of jet function: the semi-inclusive jet function $J_i(z, \omega_J, \mu)$, which describes how a parton $i$ is transformed into a jet with a jet radius $R$ and energy fraction $z = \omega_J/\omega$, with $\omega_J$ and $\omega$ being the large light-cone momentum component of the jet and the corresponding parton $i$ that initiates the jet, respectively. Within the framework of Soft Collinear Effective Theory (SCET) we calculate both $J_q(z, \omega_J, \mu)$ and $J_g(z, \omega_J, \mu)$ to the next-to-leading order (NLO) for cone and anti-k$_{\rm T}$ algorithms. We demonstrate that the renormalization group (RG) equations for $J_i(z, \omega_J, \mu)$ follow exactly the usual DGLAP evolution, which can be used to perform the $\ln R$ resummation for {\it inclusive} jet cross sections with a small jet radius $R$. We clarify the difference between our RG equations for $J_i(z, \omega_J, \mu)$ and those for the so-called unmeasured jet functions $J_i(\omega_J, \mu)$, widely used in SCET for {\it exclusive} jet production. Finally, we present applications of the new semi-inclusive jet functions to inclusive jet production in $e^+e^-$ and $pp$ collisions. We demonstrate that single inclusive jet production in these collisions shares the same short-distance hard functions as single inclusive hadron production, with only the fragmentation functions $D_i^h(z, \mu)$ replaced by $J_i(z, \omega_J, \mu)$. This can facilitate more efficient higher-order analytical computations of jet cross sections. We further match our $\ln R$ resummation at both LL$_{R}$ and NLL$_{R}$ to fixed NLO results and present the phenomenological implications for single inclusive jet production at the LHC.
  • We develop a version of Soft Collinear Effective Theory (SCET) which includes finite quark masses, as well as Glauber gluons that describe the interaction of collinear partons with QCD matter. In the framework of this new effective field theory, labeled SCET$_{\mathrm{M,G}}$, we derive the massive splitting functions in the vacuum and the QCD medium for the processes $Q\to Qg$, $Q\to gQ$ and $g\to Q\bar Q$. The numerical effects due to finite quark masses are sizable and our results are consistent with the traditional approach to parton energy loss in the soft gluon emission limit. In addition, we present a new framework for including the medium-induced full splitting functions consistent with next-to-leading order calculations in QCD for inclusive hadron production. Finally, we show numerical results for the suppression of $D$- and $B$-mesons in heavy ion collisions at $\sqrt{s_{\mathrm{NN}}}=5.02$ TeV and 2.76 TeV and compare to available data from the LHC.
  • We review the definition of semi-inclusive jet functions within Soft Collinear Effective Theory (SCET) and their application to inclusive jet cross sections. As an example, we consider both the inclusive production of jets and the jet fragmentation function in proton-proton collisions. The semi-inclusive jet functions satisfy renormalization group (RG) equations which take the form of standard timelike DGLAP evolution equations, analogous to collinear fragmentation functions. By solving these RG equations, the resummation of potentially large single logarithms $(\alpha_s \ln R)^n$ can be achieved. We present numerical results at NLO+NLL$_R$ accuracy and compare to existing data from the LHC.
  • Jim Alexander, Marco Battaglieri, Bertrand Echenard, Rouven Essig, Matthew Graham, Eder Izaguirre, John Jaros, Gordan Krnjaic, Jeremy Mardon, David Morrissey, Tim Nelson, Maxim Perelstein, Matt Pyle, Adam Ritz, Philip Schuster, Brian Shuve, Natalia Toro, Richard G Van De Water, Daniel Akerib, Haipeng An, Konrad Aniol, Isaac J. Arnquist, David M. Asner, Henning O. Back, Keith Baker, Nathan Baltzell, Dipanwita Banerjee, Brian Batell, Daniel Bauer, James Beacham, Jay Benesch, James Bjorken, Nikita Blinov, Celine Boehm, Mariangela Bondí, Walter Bonivento, Fabio Bossi, Stanley J. Brodsky, Ran Budnik, Stephen Bueltmann, Masroor H. Bukhari, Raymond Bunker, Massimo Carpinelli, Concetta Cartaro, David Cassel, Gianluca Cavoto, Andrea Celentano, Animesh Chaterjee, Saptarshi Chaudhuri, Gabriele Chiodini, Hsiao-Mei Sherry Cho, Eric D. Church, D. A. Cooke, Jodi Cooley, Robert Cooper, Ross Corliss, Paolo Crivelli, Francesca Curciarello, Annalisa D'Angelo, Hooman Davoudiasl, Marzio De Napoli, Raffaella De Vita, Achim Denig, Patrick deNiverville, Abhay Deshpande, Ranjan Dharmapalan, Bogdan Dobrescu, Sergey Donskov, Raphael Dupre, Juan Estrada, Stuart Fegan, Torben Ferber, Clive Field, Enectali Figueroa-Feliciano, Alessandra Filippi, Bartosz Fornal, Arne Freyberger, Alexander Friedland, Iftach Galon, Susan Gardner, Francois-Xavier Girod, Sergei Gninenko, Andrey Golutvin, Stefania Gori, Christoph Grab, Enrico Graziani, Keith Griffioen, Andrew Haas, Keisuke Harigaya, Christopher Hearty, Scott Hertel, JoAnne Hewett, Andrew Hime, David Hitlin, Yonit Hochberg, Roy J. Holt, Maurik Holtrop, Eric W. Hoppe, Todd W. Hossbach, Lauren Hsu, Phil Ilten, Joe Incandela, Gianluca Inguglia, Kent Irwin, Igal Jaegle, Robert P. Johnson, Yonatan Kahn, Grzegorz Kalicy, Zhong-Bo Kang, Vardan Khachatryan, Venelin Kozhuharov, N. V. Krasnikov, Valery Kubarovsky, Eric Kuflik, Noah Kurinsky, Ranjan Laha, Gaia Lanfranchi, Dale Li, Tongyan Lin, Mariangela Lisanti, Kun Liu, Ming Liu, Ben Loer, Dinesh Loomba, Valery E. Lyubovitskij, Aaron Manalaysay, Giuseppe Mandaglio, Jeremiah Mans, W. J. Marciano, Thomas Markiewicz, Luca Marsicano, Takashi Maruyama, Victor A. Matveev, David McKeen, Bryan McKinnon, Dan McKinsey, Harald Merkel, Jeremy Mock, Maria Elena Monzani, Omar Moreno, Corina Nantais, Sebouh Paul, Michael Peskin, Vladimir Poliakov, Antonio D Polosa, Maxim Pospelov, Igor Rachek, Balint Radics, Mauro Raggi, Nunzio Randazzo, Blair Ratcliff, Alessandro Rizzo, Thomas Rizzo, Alan Robinson, Andre Rubbia, David Rubin, Dylan Rueter, Tarek Saab, Elena Santopinto, Richard Schnee, Jessie Shelton, Gabriele Simi, Ani Simonyan, Valeria Sipala, Oren Slone, Elton Smith, Daniel Snowden-Ifft, Matthew Solt, Peter Sorensen, Yotam Soreq, Stefania Spagnolo, James Spencer, Stepan Stepanyan, Jan Strube, Michael Sullivan, Arun S. Tadepalli, Tim Tait, Mauro Taiuti, Philip Tanedo, Rex Tayloe, Jesse Thaler, Nhan V. Tran, Sean Tulin, Christopher G. Tully, Sho Uemura, Maurizio Ungaro, Paolo Valente, Holly Vance, Jerry Vavra, Tomer Volansky, Belina von Krosigk, Andrew Whitbeck, Mike Williams, Peter Wittich, Bogdan Wojtsekhowski, Wei Xue, Jong Min Yoon, Hai-Bo Yu, Jaehoon Yu, Tien-Tien Yu, Yue Zhang, Yue Zhao, Yiming Zhong, Kathryn Zurek
    Aug. 30, 2016 hep-ph, hep-ex, nucl-ex, astro-ph.CO
    This report, based on the Dark Sectors workshop at SLAC in April 2016, summarizes the scientific importance of searches for dark sector dark matter and forces at masses beneath the weak-scale, the status of this broad international field, the important milestones motivating future exploration, and promising experimental opportunities to reach these milestones over the next 5-10 years.
  • We present the nuclear transverse momentum broadening for Drell-Yan lepton pair production in p+A collisions at next-to-leading order (NLO) in powers of strong coupling constant $\alpha_s$. We verify that the transverse momentum weighted differential cross section in NLO perturbative QCD (pQCD) at twist-4 can be factorized into the convolution of parton distribution function of an active parton in the projectile proton, a twist-4 multiparton correlation function of the target nucleus, and perturbatively calculable hard coefficient function. We identify a QCD evolution equation for such a twist-4 nuclear gluon-quark correlation function, and verify its universality -- its process independence. This finding demonstrates the prediction power of pQCD factorization approach for studying parton multiple scattering in nuclear medium.
  • The jet fragmentation function describes the longitudinal momentum distribution of hadrons inside a reconstructed jet. We study the jet fragmentation function in proton-proton collisions in the framework of soft-collinear effective theory (SCET). We find that, up to power corrections, the jet fragmentation function can be expressed as the ratio of the fragmenting jet function and the unmeasured jet function. Using renormalization group techniques, we are able to resum large logarithms of jet radii R in the perturbative expansion of the cross section. We use our theoretical formalism to describe the jet fragmentation functions for light hadron and heavy meson production measured at the Large Hadron Collider (LHC). Our calculations agree very well with the experimental data for the light hadron production. On the other hand, although our calculations for the heavy meson production inside jets are consistent with the PYTHIA simulation, they fail to describe the LHC data. We find that the jet fragmentation function for heavy meson production is very sensitive to the gluon-to-heavy-meson fragmentation function.
  • Recent advances in soft-collinear effective theory with Glauber gluons have led to the development of a new method that gives a unified description of inclusive hadron production in reactions with nucleons and heavy nuclei. We show how this approach, based on the generalization of the DGLAP evolution equations to include final-state medium-induced parton showers, can be combined with initial-state effects for applications to jet quenching phenomenology. We demonstrate that the traditional parton energy loss calculations can be regarded as a special soft-gluon emission limit of the general QCD evolution framework. We present phenomenological comparison of the SCET$_{\rm G}$-based results on the suppression of inclusive charged hadron and neutral pion production in $\sqrt{s_{NN}}=2.76$ TeV lead-lead collisions at the Large Hadron Collider to experimental data. We also show theoretical predictions for the upcoming $\sqrt{s_{NN}} \simeq 5.1$ TeV Pb+Pb run at the LHC.
  • Recent measurements of the centrality and rapidity dependence of single inclusive jet production in p+Pb collisions at the LHC have revealed large and non-trivial nuclear modification of the production cross section for this process. In this paper, we explore to what extent such nuclear modification can be understood by the framework of standard cold nuclear matter effects, in particular initial-state cold nuclear matter energy loss. We demonstrate quantitatively that theoretical calculations which include medium-induced radiative corrections can describe rather reasonably the attenuation of the jet production yields in the large transverse momentum region in d+Au collisions at RHIC and p+Pb collisions at the LHC for central to semi-central collisions. We further show that the observed scaling behavior of the nuclear modification factor as a function of the total jet energy $p_T\, {\cosh} (y)$ for various rapidity intervals has a natural explanation in the picture of cold nuclear matter energy loss. On the other hand, the observed enhancement in peripheral collisions is not described in this picture and could have a different origin.
  • Tagged jet measurements in high energy hadronic and nuclear reactions provide constraints on the energy and parton flavor origin of the parton shower that recoils against the tagging particle. Such additional insight can be especially beneficial in illuminating the mechanisms of heavy flavor production in proton-proton collisions at the LHC and their modification in the heavy ion environment, which are not fully understood. With this motivation, we present theoretical results for isolated-photon-tagged and B-meson-tagged b-jet production at center-of-mass energy 5.1 TeV for comparison to the upcoming lead-lead data. We find that photon-tagged b-jets exhibit smaller momentum imbalance shift in nuclear matter, and correspondingly smaller energy loss, than photon-tagged light flavor jets. Our results show that B-meson tagging is most effective in ensuring that the dominant fraction of recoiling jets originate from prompt b-quarks. Interestingly, in this channel the large suppression of the cross section is not accompanied by a significant momentum imbalance shift.
  • We study the cross section for vector boson ($W^{\pm}/Z^0/\gamma^*$) production in polarized nucleon-nucleon collisions for low transverse momentum of the observed vector boson. For the case where one measures the transverse momentum and azimuthal angle of the vector bosons, we present the cross sections and the associated spin asymmetries in terms of transverse momentum dependent parton distribution functions (TMDs) at tree level within the TMD factorization formalism. To assess the feasibility of experimental measurements, we estimate the spin asymmetries for $W^{\pm}/Z^0$ boson production in polarized proton-proton collisions at the Relativistic Heavy Ion Collider (RHIC) by using current knowledge of the relevant TMDs. We find that some of these asymmetries can be sizable if the suppression effect from TMD evolution is not too strong. The $W$ program at RHIC can, thus, test and constrain spin theory by providing unique information on the universality properties of TMDs, TMD evolution, and the nucleon structure. For example, the single transverse spin asymmetries could be used to probe the well-known Sivers function $f_{1T}^{\perp q}$, as well as the transversal helicity distribution $g_{1T}^{q}$ via the parity-violating nature of $W$ production.
  • We study the Sivers asymmetry in semi-inclusive hadron production in deep inelastic scattering. We concentrate on the contribution from the photon-gluon fusion channel at ${\cal O}(\alpha_{em}^2 \alpha_s)$, where three-gluon correlation functions play a major role within the twist-3 collinear factorization formalism. We establish the correspondence between such a formalism with three gluon correlation functions and the usual transverse momentum dependent (TMD) factorization formalism at moderate hadron transverse momenta. We derive the coefficient functions used in the usual TMD evolution formalism related to the quark Sivers function expansion in terms of the three-gluon correlation functions. We further perform the next-to-leading order calculation for the transverse-momentum-weighted spin-dependent differential cross section, and identify the off-diagonal contribution from the three-gluon correlation functions to the QCD collinear evolution of the twist-3 Qiu-Sterman function.