Research Paper Notes on Phase Space Integral

Research Paper Notes on Phase Space Integral

本文档除了包括推导, 疑惑 外,做 读书重点 的记录

文献列表

 * Phys. Rev. Lett. 54, 510 (1985) Single-Particle Momentum Distribution at High Energies and Concept of Partition Temperature. by T.T. Chou, Chen Ning Yang, E. Yen


 * Phys. Rev. D 46, 160 (1992) Semi-inclusive rapidity distribution and a critical analysis of the concept of partition temperature. by Yogiro Hama, Michael Plumer


 * Phys. Rev. D 87, 114007 (2013) Tsallis fits to pt spectra and multiple hard scattering in pp collisions at the LHC. by Cheuk-Yin Wong, Grzegorz Wilk

Phys. Rev. Lett. 54, 510 (1985) Single-Particle Momentum Distribution at High Energies and Concept of Partition Temperature.
这里的一个重要的思想是系综的微正则分布导致单粒子的统计分布.

Phys. Rev. D 46, 160 (1992) Semi-inclusive rapidity distribution and a critical analysis of the concept of partition temperature.
本工作就是上面杨振宁一文思想的具体数学推导.具体可以参考Yogiro的详细推导手稿

本文的物理思想是,从微正则系综和粒子数守恒出发能够一般的得到正则系综从而单粒子分布具有温度的概念. 这里通过一个相对论重离子碰撞中粒子分布的具体计算,得到解析的结果,并和实验数据分析对照确认推导的合理性.

(2.4)

文中指出,比值(2.4)不依赖与$$\beta$$.因为如果依赖的话,那么本文的结果就相对依赖于假设(2.11). 而结果的一般性主要体现在本文的结果不依赖于(2.11)的具体形式.

具体考虑分子分母的形式.比如分子是粒子占有数$$n_k$$的加权平均值.每个权重是一个指数函数的积分,其中指数上涉及$$\beta$$的项为
 * $$\prod_l \left(\exp(-\beta\sqrt{p_{Tl}^2+m^2}\cosh y)\right)^{n_l}

=\prod_l \left(\exp(-\beta E_l)\right)^{n_l} =\exp\left(-\sum_l \beta {n_l}E_l\right) =\exp(- \beta W)$$ 其中最后一步是因为能量守恒的$$\delta$$函数.容易看到这个与$$\beta$$有关的因子对分子和分母是一样的.所以最后结果严格的与$$\beta$$无关.但是正如Yogiro指出的,做了近似后的结果是可能与$$\beta$$有关的,但是不应该敏感的依赖于$$\beta$$.

(2.7)

指数上有因子$$-(s \cosh y-t \sinh y)$$并对y积分,为了这个积分收敛,指数上的实部需要为复数,故$$s$$的实部要大于$$t$$的实部,即$$\epsilon_0 > \epsilon_1$$.

注意到在表达式$$A,B$$中对$$s$$的积分,指数上$$s$$前面的系数是负定的,而原始的积分是在虚轴上的.比较约当引理的证明,原始的积分是在实轴上,而指数上被积变量前的系数是纯虚数$$ia$$,其中$$a$$为正数或者复数决定了积分可以从实轴的上方或者下方绕行.按证明约当引理的思路,我们知道$$s$$在虚轴上的积分可以在虚轴的右侧(实轴正方向)的大圆上绕行.

对$$t$$的积分情况更复杂,一般的证明并不是很显然,但是我们延续上面的思路,不难想象如果对$$t$$的在大圆上绕行的积分始终同时伴随对$$s$$的大圆上绕行的积分,而且保证在大圆上的两个积分路径上指数的实部始终小于零,那么在大圆上的绕行就是可行的.

(2.9-10)

这里不考虑横动量守恒,从而对$$u_T$$的积分看成$$du_T \rightarrow \delta(u_T)du_T$$.

(2.11)

Yogiro指出,这里的第二项只是为了有效的斩断横动量,因为对于jet而言,粒子(相对于jet方向垂直)的横动量非常小,可以视为零.

(3.1)

按(2.22)对$$p_T$$部分积分,考虑零质量$$y=\eta$$,$$p_T=m_T$$,$$m\delta=0$$,由(2.22):


 * $$\frac{dn}{d\eta}=\frac{dn}{dy}

\propto \delta^2\int 2\pi p_T dp_T \exp(-\delta p_T) \exp(-p_T \cosh y/T_p) =\delta^2\int 2\pi p_T dp_T \exp(-\delta p_T) \exp(-p_T \cosh y/T_p) = \delta^2 2\pi \int p_Tdp_T \exp[-(\delta+\cosh y/T_p)p_T] = \delta^2 \frac{\partial}{\partial a} \int dx \exp(-a x) = \delta^2 \frac{1}{a^2} =\frac{1}{(1+\cosh y/\delta T_p)^2}$$

其中定义$$a=\delta+\cosh y/T_p$$.

(3.2)

这里因子$$2$$来源于两阶导数,$$T_p$$来自$$\delta^2$$没有消去的部分,由于积分中$$E=p_T \cosh\eta$$,似乎分子分母上都是除以$$\cosh^2\eta$$.

Table I

按CYY一文的计算采用的$$\langle p_T(0)\rangle$$和$$W$$,(2.23)得到的结果与CYY的拟合得到的结果很吻合.

这里最后三列,分别是通过实验数据得到的$$\langle p_T(0)\rangle$$和通过对$$dn/d\eta$$积分到入射能量决定的最大快度而得到的与CYY不同的jet的总能量$$W$$,来计算$$T_p$$,得到很不一样的结果.

(A4)

首先注意到$$E_k=E_{l=k}$$. 考虑$$a_l=q_l \exp[-(E_l s \cdots)+i p_{Tl}\cdot u_T+iv]$$,分子相当于计算
 * $$\Pi_l \sum_{n_l} n_k \frac{a_l^{n_l}}{n_l !}

=a_k \frac{\partial}{\partial a_k}\Pi_l \sum_nl \frac{a_l^{n_l}}{n_l !} =a_k \frac{\partial}{\partial a_k}\Pi_l \exp(a_l) =a_k \Pi_l \exp(a_l) =a_k \exp(\sum_l a_l)$$

(A5)

这里对$$v$$的积分限为$$(0,2\pi)$$因为$$n$$和$$n_l$$都是整数,所以展开对象是周期函数.换言之,这里是傅里叶级数而非傅里叶积分. 利用(A5)的提示,得到:
 * $$\int_0^{2\pi} dv e^{-i(n-1)v}\exp[F e^{iv}]

=\sum_l \frac{F^l}{l!}\int_0^{2\pi} dv e^{-i(n-1-l)v} =\sum_l \frac{F^l}{l!}\delta(n-1-l) =\frac{F^{n-1}}{(n-1)!}$$

(B6)

这里Ei对应指数积分.

Fig.1

这里有一个过程要利用复共轭来把复数的反超三角函数的实部和虚部分开.具体可以参考这个计算过程和Yogiro的详细推导手稿.

这里我们给出一个对Yogiro推导手稿(见上面的链接)的总结.首先,推导的主要思路包括(1)按(2.15)的变换,最好的讨论方法是把所有的复数都写成分量的形式,这样在$$s,t$$一头由于复数的实部固定,有两个变量,在$$\eta\xi$$一头暂时有4个分量.手稿首先尝试从$$s,t$$一头出发讨论$$\eta\xi$$的范围,然后反过来从$$\eta\xi$$出发验证是否的确完全覆盖了$$s,t$$参数空间(2)随着讨论,会要自然放宽(2.6)中对$$s$$和$$t$$的积分路径,按复变函数积分理论,我们的确被允许这样做.比如从$$\eta\xi$$一头出发,固定$$\eta$$讨论$$\xi$$变化导致的在$$s,t$$空间的曲线,手稿通过去$$\epsilon_1=0$$来消去一个变量.但是这是不难验证,这时$$\epsilon_0$$不再是常数了.

我们这里对手稿最后的在$$S,T$$空间图做出说明.这个图的每根实线对应一组固定$$a,b$$,实线上的点对应不同的$$\alpha$$值,而$$\gamma$$由$$\epsilon_1=0$$不是自由变量.只需要考虑$$a>0$$的情况,在$$b>0$$和$$b<0$$时,曲线分别处于$$T$$轴的左右两边,曲线分别有四个角的无穷大的渐进行为.但是考虑覆盖全空间,还需要考虑这些曲线是否覆盖在原点附近的$$S,T$$参数空间.这些曲线与$$S$$轴有截距,这个截距在$$b$$趋于零的时候趋于零.这些曲线在$$|b|>a$$时与$$S=T$$相交,交点在$$b$$趋于零时趋于零.如果直接取$$b=0$$得到一根在$$T$$轴上从$$a$$点向上的割线.这个割线的渐进行为通过考察$$b$$趋于零时的渐进行为得到.考虑曲线与$$T=a$$的交点,在$$b$$很小时,对应$$\alpha$$趋于正负零从而$$S$$趋于正负零.所以曲线在$$b$$不但在$$S$$轴上的截距趋近于原点,而且无限逼近$$T=a,S=0$$点.但是注意到即便在$$b$$趋于零时,当$$\alpha$$趋于正负无限大时曲线的$$S$$值仍然趋于无限大,如图所示其渐近线同样是$$S=\pm T$$,而不是包裹上述割线.这个渐进行为使得$$S,T$$参数空间存在如图虚线上(下)方的空洞.就此,手稿得到结论,需要通过把$$a$$推到无穷大才能正常覆盖全部$$T,S$$空间.

(C10)

文中计算了沿着实轴方向的两阶导数,但是鞍点积分是沿着虚轴方向的.因为$$I_\alpha$$和$$K_\alpha$$分别是指数增加和指数衰减,函数沿着实轴方向导数为零的点为最小值.接着文中用沿着实轴的两阶导数来估算鞍点积分的结果. 这是因为对解析函数,导数也是解析函数,所以极值必然是鞍点,而在鞍点附近对小的复数展开满足$$\Phi(\xi)=\Phi(\xi_0)+\frac{1}{2}\Phi''(\xi-\xi_0)^2$$.

这里总结一下积分的办法,首先是直接积分,其次考虑特殊函数(利用特殊函数的积分表示),再次有留数法积分,鞍点积分法,换元法,最后可以用mathematica积分.

Phys. Rev. D 87, 114007 (2013) Tsallis fits to pt spectra and multiple hard scattering in pp collisions at the LHC.
这里的问题是是否可以用jet的分布和粒子发射的微正则分布导出末态重子的分布形式.