重夸克偶素作为一个非常特殊的物理系统，在探究强互作用相关的微扰和非微扰方面都有着重要的意义。近年来更是得到了广泛的研究和探讨。我们在非相对论量子色动力学框架下应用围绕计算研究了一些重夸克偶素的产生过程。　　我们考察了次领头阶水平上HERA上光生J/Ψ介子的极化。结果显示J/Ψ粒子产生的横动量分布和能量分布与实验测量值符合得不好。由于极化参数的能量z分布对重整化和因子化标度的不同取值依赖很大，所以其理论结果在次领头阶无法给出精确的预言。在Pt＞3GeV的区域，J/Ψ粒子的Pt分布的不确定性比较小，但是甚至在这个区域最终的数值结果还是没有能描述实验数据。　　我们还研究了色单态框架下J/Ψ(Υ)在z玻色子衰变中的单举产生。我们计算了Z→ Quarkonium+QQ这一过程的次领头阶修正，这是在数图水平上贡献最大的道。计算过程中z玻色子和费米子耦合顶角的矢量部分和轴矢部分是分开处理的。结果表明，矢量部分和轴矢部分有着相同的K因子(1.13)，参数μ=2mc，mc=1.5GeV。当使用BLM方案确定重整化标度时，K因子降到0.918，参数μBLM=2.28GeV，mc=1.5GeV。通过计入次重要的过程的贡献，主要是光子和胶子碎裂过程，Z→ J/Ψprompt+X的分支比可以达到(7.3～10.0)×10-5。这里不确定性来自于重整化标度和Charm夸克质量。这结果约是实验测量中心值的一半。进一步的，J/Ψ粒子的能量分布也不能与实验值一致。因此在次领头阶，色单态模型模型下的z→ Quarkonium+Q(Q)无法给出实验值需要的理论预言。这里应该有来自八重态的贡献。我们定义量Rc(c)=Γ(Z→J/ψc(c)X)/Γ(Z→J/ψX)，并在色单态模型下得到Rcc=0.84，在色八重态模型下有Rcc=0.49。Rcc值的测量将有助于澄清色八重态的贡献。　　这里我们研究了重夸克碎裂到重夸克偶素的碎裂函数。其小z区域的非物理的碎裂函数值用三种情况替换并研究了影响。数值结果表明，这三种情况间的差异随着碎裂标度逐渐增大的演化而变小。而且这种差异对唯象的计算研究几乎没有影响。使用NLO的重夸克到重夸克偶素的碎裂函数，我们研究了Υ伴随一个Bottom夸克喷注在强子对撞机上的产生。并和固定阶微扰论的领头阶结果相对比。在大Pt区域，碎裂机制的结果是领头阶过程p+p→Υ+b(b)+X的结果的1.5倍。