Codex @codex  0

Let $u_i=⟨v_i⟩$, ${\sigma }_{ij}=⟨(v_i-u_i)(v_j-u_j)⟩$, and $P_{ij}=\rho {\sigma }_{ij}$. The Jeans equation obtained by taking the first velocity moment of the Collisionless Boltzmann equation isFor an equilibrium population multiply by $x_i$ and integrate. Integration by parts givesAssume the integrals exist and the boundary flux vanishes, as for an isolated finite system with adequate decay. With the stellar kinetic-energy tensor, stellar pressure-energy tensor and stellar potential-energy tensor defined bywe obtain the equilibrium tensor virial theoremThe stellar kinetic-energy tensor describes ordered motion, whereas the stellar pressure-energy tensor describes random motion. In particular, $T_{ij}=0$ does not mean the stars have no kinetic energy. If a surface flux is retained, the right side of the displayed equilibrium identity is that surface tensor; it cannot be discarded for an arbitrary finite aperture. For an isolated self-gravitating population, $W_{ij}$ is symmetric by pairwise interchange in the Newtonian interaction integral, and the time-dependent identity is $\frac{1}{2}{\ddot{I}}_{ij}=2T_{ij}+{\Pi }_{ij}+W_{ij}$ with $I_{ij}=\int \rho x_i{x}_j{d}^{3}x$.

For the spherical Newtonian gravitational potential, ${\partial }_j\Phi ={\Phi }^{\prime }(r)x_j/r$. In the pressure-supported case the tensor virial theorem gives ${\Pi }_{ii}=\int \rho {\Phi }^{\prime }(r)x_i^2/r{d}^{3}x$. The notation ${\Pi }_{RR}$ used in this question is the planar virial trace ${\Pi }_{xx}+{\Pi }_{yy}$, not just the integral of the local cylindrical radial velocity dispersion. In cylindrical coordinates this trace includes both $v_R^2$ and $v_{\varphi }^2$. Summing the two Cartesian components therefore givesThe distinction explains the factor of two in the spherical limit.

Use ordinary spherical coordinates, $R=r\sin \theta$, $z=r\cos \theta$. The mass density separates as ${\rho }_0{r}^{-\gamma }h(\theta )$, whereBoth potential-energy integrals have the same radial factor $2\pi {\rho }_0{\int }_0^{\infty }{r}^{3-\gamma }{\Phi }^{\prime }(r)dr$. When this factor is finite and positive it cancels, leaving the scale-free tracer virial ratioThe angular factor in the printed intermediate formula has its sine and cosine interchanged. Keeping its sine/cosine weights while using that printed factor would instead produce the opposite sign in the final small-flattening correction. Merely measuring $\theta$ from the equatorial plane cannot fix this: the measure and both moment weights would have to change as well.

The angular ratio is finite and positive for every fixed $q>0$ and finite real $\gamma$, because $h$ is positive and bounded above and below. However, the claim that the unrestricted global virial integrals are always well-defined needs a hypothesis. In a point-mass Newtonian gravitational potential, the common radial factor is proportional to ${\int }_0^{\infty }{r}^{1-\gamma }dr$: convergence at zero requires $\gamma <2$, while convergence at infinity requires $\gamma >2$. There is no such $\gamma$. If the global energies diverge, the boxed angular ratio is still the ratio of the two potential-energy integrals with identical spherical inner and outer cutoffs, independent of the cutoffs. It is not a quotient of two finite global stellar pressure-energy tensors; applying a finite-aperture tensor virial theorem additionally requires its surface terms. These convergence and boundary qualifications are separate from the printed angular typo.

Set $ϵ=q^{-2}-1$. For modest flattening, the Taylor expansion is $h=1-(\gamma /2)ϵ{\cos }^2\theta +O({ϵ}^2)$. The required elementary angular integrals areThese follow directly by $t=\cos \theta$, or from the supplied gamma function identity. Hence the numerator is $(4/3)[1-\gamma ϵ/10+O({ϵ}^2)]$ and the denominator is $(2/3)[1-3\gamma ϵ/10+O({ϵ}^2)]$. The flattening–anisotropy virial relation becomesFor a decreasing tracer mass density ($\gamma >0$), oblate flattening ($q<1$) requires more random kinetic support per in-plane direction than vertically. In a spherical Newtonian gravitational field, directional orbital anisotropy supplies the flattening rather than an anisotropic force law. A steeper radial mass density amplifies the required velocity dispersion anisotropy. A prolate tracer reverses the leading inequality. The relation constrains integrated velocity dispersions, not a unique local galactic distribution function.

The completed Dirichlet L-function isFor nonprincipal primitive Dirichlet characters, termwise Mellin transformation initially in $\Re s>1$ givesThe Dirichlet character theta function decays exponentially at infinity; its transformation makes it decay faster than any power at zero. Hence the integral is entire in $s$. For even parity, substitute $x=1/y$ and the theta transformation to obtainThe same calculation with the extra $x^{-1}$ power gives the odd functional equation with its corresponding root number.

The gamma function has no zeros and has simple poles at nonpositive integers. Thus the nontrivial zeros of $L$ and $\Lambda$ coincide with multiplicities. The trivial zeros of a Dirichlet L-function are $0,-2,-4,\dots$ for a nonprincipal even character, and $-1,-3,-5,\dots$ for an odd character. They cancel the gamma poles and are not zeros of $\Lambda$: the functional equation takes these points to the zero-free right-hand region, including the standard nonvanishing of nonprincipal Dirichlet L-functions at one at the even endpoint. The canceled zeros are simple.

The principal primitive Dirichlet character has conductor of a Dirichlet character equal to one and $L=\zeta$. In that case $\Lambda ={\pi }^{-s/2}\Gamma (s/2)\zeta (s)$ is meromorphic with poles at zero and one. Its canceled trivial zeros begin at $-2$, while $\zeta (0)=-1/2$ is not zero. Multiplication by $s(s-1)/2$ produces the entire Riemann xi function used below.

Watson's lemma applies to Laplace-type endpoint integrals such as ${\int }_0^b{e}^{-\lambda t}f(t)dt$ as $\lambda \to +\infty$, when $f$ has a local algebraic expansion $f(t)\sim {\sum }_j{c}_j{t}^{p+j-1}$ with $p>0$ and suitable integrability or exponential-growth control away from the endpoint. Each term contributes $c_j\Gamma (p+j){\lambda }^{-p-j}$. The exponential selects a neighborhood of $t=0$.

Laplace method instead treats ${\int }_a^{b}g(t)e^{\lambda \varphi (t)}dt$ by locating the dominant maxima of a real phase $\varphi$. A smooth nondegenerate interior maximum produces a neighborhood governed by a Gaussian function of width ${\lambda }^{-1/2}$; endpoint or degenerate maxima require their corresponding local scaling. Contributions away from a strict dominant maximum must be smaller. These methods overlap after suitable changes of variable, but their usual starting forms emphasize an endpoint amplitude expansion and an exponential maximum, respectively.

For the Gamma function, put $t=xu$. ThenThe phase has its unique maximum at $u=1$, with value $-1$ and second derivative $-1$. Set $u=1+s/\sqrt{x}$. Near this maximum,Combining the amplitude and exponential expansions in Laplace method givesLocalization near $u=1$ justifies extending the limits for the Gaussian integral; the distant endpoint is not being expanded uniformly by this local series. Odd terms integrate to zero. The moments of the Gaussian integral satisfy $I_1=I_0$, $I_2=3I_0$, $I_3=15I_0$, so the relative correction isThe two terms of the Stirling formula, obtained by the Stirling correction by Gaussian moments, are therefore

With the determinant-normalized slash operator, let $W_N=\left(\begin{array}{cc}0& -1\\ N& 0\end{array}\right)$ and $Jf=i^k(f{|}_k{W}_N)=g$. Conjugation by $W_N$ preserves ${\Gamma }_1(N)$, and rational slash operators preserve modular cusp holomorphy and vanishing. Hence $g$ is also a cusp form at that level. Direct substitution gives $J^{2}f=f$; the factor $i^k$ cancels the central weight sign.

Put $F(y)=f(iy/\sqrt{N})$ and $G(y)=g(iy/\sqrt{N})$. The defining formula gives the exact relationsIn the initial half-plane of absolute convergence, termwise integration of the Fourier series and the gamma function yield the Mellin transform of a cusp-form L-functionThe scaling of $y$ accounts for $N^{s/2}$ in the completion. At infinity $F$ and $G$ decay exponentially. At zero the boxed relation expresses $F$ as a power times an exponentially decaying function of $1/y$. Thus this integral converges locally uniformly for every complex $s$, including after differentiation in $s$, and defines an entire function.

Splitting at one and changing $y$ to $1/y$ in the lower integral givesApplying the same formula to $g$ interchanges $F,G$, because $J^2=1$. It proves the phase-normalized Fricke functional equationThe entire function here is the completion, despite the apparent poles of the gamma factor in its initial product formula.

Take an oriented closed string in flat, critical bosonic string theory, with ${\alpha }^{\prime }=1/(2\pi T)$ and a mostly-plus target Minkowski metric. After continuation of the worldsheet to Euclidean signature, the Polyakov path integral sums over embeddings $X$ and worldsheet metrics, divided by worldsheet diffeomorphisms and Weyl transformations. In a flat target its kinetic action isTarget-time continuation or analytic continuation of external momenta defines the Lorentzian scattering amplitude; a naive real Euclidean Gaussian for timelike $X^0$ would not be convergent.

An external tachyon is represented by the tachyon vertex operator $V_p(z,\bar{z})=:e^{ip\cdot X}:$, with conformal weights $h=\bar{h}={\alpha }^{\prime }{p}^2/4$. The physical integrated vertex has $(h,\bar{h})=(1,1)$, so $p^2=4/{\alpha }^{\prime }$. Schematically its tree amplitude isThe Faddeev-Popov determinant from conformal gauge is represented by the worldsheet ghost fields. At tree level the worldsheet is a Riemann sphere, whose unpunctured complex structure has no moduli. Its residual conformal automorphisms are the Möbius transformations, $PSL(2,\mathbb{C})$. Fix three insertion points, accompanying their unintegrated vertices by the required $c\tilde{c}$ ghost factors. The remaining $N-3$ complex insertion positions are integrated over the sphere; equivalently one integrates all positions and divides by the residual conformal group.

The embedding fields are free, withTheir zero-mode integral gives momentum conservation, and their nonzero-mode Gaussian integral gives the Koba-Nielsen factorThis explains the sphere tachyon position integral without needing its evaluation. For four tachyons the resulting Virasoro–Shapiro amplitude displays the exchanged string spectrum directly.

Use all-incoming external momenta and introduce $u=-(p_1+p_4{)}^2/(8\pi T)$ alongside the two printed dimensionless Mandelstam variables. Since $p_a^2=8\pi T$ and ${\sum }_a{p}_a=0$, one obtains $s+t+u=-4$. The physical center-of-mass energy squared in the $s$ channel is $s_{\mathrm{phys}}=8\pi Ts=4s/{\alpha }^{\prime }$. The $t$ channel measures the analogous crossed momentum transfer, with sign determined by the mostly-plus convention. Rewriting the Gamma factors in a symmetric form givesThe Gamma function poles imply, at generic fixed values of the other invariant,The denominator Gamma factors can remove residues at special intersecting channel kinematics; the statement concerns a generic single-channel limit. An $s$-channel pole occurs when the intermediate momentum $p_1+p_2$ satisfies the mass-shell condition for a closed-string state:The pole at $s=-1$ exchanges the ground-state tachyon, the pole at $s=0$ exchanges massless states, and the positive integer poles exchange the infinite massive tower. The $t$-channel interpretation is the crossed version. Factorization means that each residue is a sum of products of couplings to intermediate physical states that couple to the chosen external particles. It need not expose every representation at that mass.

For example, the Gamma function recurrence and Gamma function residue at a nonpositive integer give the dimensionless Virasoro–Shapiro amplitude pole residueThe residue polynomial has degree $2n$, consistent with maximum spin $2n$ in the exchanged level. The amplitude also has the corresponding $u$-channel poles by crossing symmetry.

The massless fields can be treated as target backgrounds rather than separate asymptotic insertions. Write the target metric as $G_{mn}={\eta }_{mn}+h_{mn}$, introduce a Kalb–Ramond field $B_{mn}$, and a dilaton $\Phi$. In conventional Euclidean signs their string nonlinear sigma model action isHere ${ϵ}^{ab}$ is the antisymmetric tensor density. The three backgrounds correspond to the graviton, antisymmetric tensor and scalar states at closed-string level one. Expanding the vacuum functional in $h,B,\Phi$, then Fourier expanding the backgrounds, produces exactly their integrated string vertex operators. Its functional derivatives therefore generate the amplitudes with massless external strings. The connected vacuum functional organizes connected amplitudes; the spacetime effective action organizes the corresponding vertices after treating massless propagation consistently.

At momenta small compared with $1/\sqrt{{\alpha }^{\prime }}$, massive string propagators can be expanded in powers of momenta over their masses. The analytic part of the amplitudes consequently determines local higher-derivative interactions, ordered by powers of ${\alpha }^{\prime }$. Massless exchange poles are retained through propagation of the massless fields, rather than expanded into local contact terms. Up to field redefinitions, the leading massless-sector action in the string-frame metric isThe terms denoted $O({\alpha }^{\prime })$ contain additional derivatives, including curvature-squared terms in the bosonic theory. The expansion concerns the massless sector around the perturbative bosonic background; the tachyon instability remains and is not cured by omitting its field from this displayed action. Thus this is a formal perturbative effective description, not a claim of a stable bosonic vacuum.

Finally split the dilaton into a constant and its variation, $\Phi ={\Phi }_0+\varphi$. By the Gauss-Bonnet theorem, its dilaton Euler-characteristic weighting on a connected closed oriented surface follows fromDefine the string coupling by $g_s=e^{{\Phi }_0}$. A genus-$g$ path integral is weighted byThe sphere carries $g_s^{-2}$, the torus carries $g_s^0$, and each extra handle adds $g_s^2$. At higher genus one integrates over complex-structure moduli as well as insertion points, with the associated antighost insertions supplying the correct moduli measure. With canonically normalized external vertices an $N$-point genus-$g$ amplitude scales as $g_s^{2g-2+N}$.

Summing connected worldsheets of every genus of a surface yields the string-loop effective action expansionEach coefficient has its own low-energy ${\alpha }^{\prime }$ expansion. The two parameters have different roles: ${\alpha }^{\prime }$ resolves finite string size through higher derivatives, while $g_s^2$ counts additional string loops.

The Bessel turning-point asymptotic comes from a cubic stationary endpoint, not an ordinary quadratic stationary point. Near $\vartheta =0$,Thus the contributing width is $\vartheta =O(n^{-1/3})$. On writing $\vartheta =(6/n{)}^{1/3}s$, the leading integral isThe oscillatory integral is understood with a vanishing damping factor. Substitution $u=s^3$ and the Gamma function Fourier integral giveThereforeThe equality uses the Gamma reflection formula. Contributions away from the degenerate endpoint are smaller. The $n^{-1/3}$ scale explains why the preceding $n^{-1/2}$ formula cannot be extended directly to zero angle.

Apply a Feynman parameter and shift the loop momentum to $\ell =p+xk$. The common denominator becomes $[{\ell }^2+M_x^2{]}^2$, whereFor positive $M_x^2$, the Gamma-integral representation $a^{-2}={\int }_0^{\infty }dss{e}^{-sa}$ and a Gaussian integral giveThis formula initially converges for $\mathrm{Re}d<4$ and defines the Euclidean massive loop integral at other dimensions by analytic continuation. Consequently,With $d=6-ϵ$, the Gamma function factor is $\Gamma (-1+ϵ/2)=-2/ϵ+O(1)$. Only its pole matters: the other factors can be evaluated at $ϵ=0$ when extracting that pole. Since ${\int }_0^{1}x(1-x)dx=1/6$,This is the one-loop two-point divergence in six-dimensional cubic scalar theory. Its polynomial momentum dependence is precisely what permits subtraction by local counterterms.

Let $A_s$ denote the surface area of an n-sphere for the unit sphere $S^s$. In this notation $s$ is the intrinsic dimension of the sphere. Cartesian factorization of the Gaussian integral gives $I_N={\pi }^{N/2}$. In hyperspherical coordinates,The substitution $t={\rho }^2$ supplies the defining Gamma function integral. Scaling the sphere radius multiplies its surface volume by $r^{N-1}$, henceThis is the volume of the sphere itself, not of its enclosed ball.

Take $D$ to denote spacetime dimension and $d=D-1$ spatial dimensions. In uncompactified flat space, Gauss's law for a charge $q$ and permittivity ${\epsilon }_d$ gives $E_r{A}_{d-1}{r}^{d-1}=q/{\epsilon }_d$. Since $E_r=-d\phi /dr$, the electric potential for $d>2$, normalized to vanish at infinity, isThus the four-dimensional electric potential falls as $1/r$, and the five-dimensional one as $1/r^2$. If a convention counts only spatial dimensions as $D$, replace $d$ by $D$ in the first formula. For $d=2$, the integral instead gives $\phi =-q\log (r/r_0)/(2\pi {\epsilon }_2)$; for $d=1$, it gives $\phi =-q|x|/(2{\epsilon }_1)$ up to a constant. Neither case has a finite zero at infinity.

For a Newtonian gravitational potential ${\Phi }_g$, define the flux-normalized coupling ${\kappa }_d$ by ${\nabla }^2{\Phi }_g={\kappa }_{d}m{\delta }^{(d)}(\mathbf{x})$ and acceleration $-\nabla {\Phi }_g$. Exactly the same flux integration gives the attractive potentialIn $d=3$, ${\kappa }_3=4\pi G_4$ reproduces $-G_{4}m/r$. If $G_D$ is defined by the $D$-dimensional Einstein-Hilbert action $S=(16\pi G_D{)}^{-1}\int \sqrt{-g}R$, the weak static Einstein field equations give ${\kappa }_d=8\pi G_D(D-3)/(D-2)$ for $D>3$. ThereforeIndeed, $R_{00}\simeq {\nabla }^2{\Phi }_g$ and $T_{00}-g_{00}T/(D-2)=(D-3)\rho /(D-2)$ in signature $(-+\cdots +)$. A Poisson-defined Newtonian theory in two spatial dimensions has a logarithmic potential, but pure three-dimensional Einstein gravity has no analogous local Newtonian point-mass force; the relativistic normalization must not be extrapolated to that exceptional dimension.

For an unwarped product with a fixed extra-dimensional metric, normalize the Einstein-Hilbert action by $\frac{1}{2}{M}_{\ast }^{n+2}\int d^{4+n}x\sqrt{-g}R$. Integrating the part containing the four-dimensional curvature over the internal space gives $\frac{1}{2}{M}_{\ast }^{n+2}\mathrm{vol}(S_R^n)\int d^{4}x\sqrt{-g_4}{R}_4$. Thus, using the reduced Planck mass,This is a volume relation for compactification; a curved sphere additionally requires a mechanism supporting and stabilizing its background. It is not by itself a proof that a sphere times flat spacetime solves a vacuum gravitational theory.

With $M_{\ast }=10^3\mathrm{GeV}$, $M_{\mathrm{Pl}}=2.435\times 10^{18}\mathrm{GeV}$ and $1{\mathrm{GeV}}^{-1}=1.97327\times 10^{-16}\mathrm{m}$, the spherical-volume factors $A_1=2\pi$, $A_2=4\pi$, $A_6=16{\pi }^3/15$ giveThe first estimate is astronomical; the second is roughly a tenth of a millimetre. If both Planck scales are instead defined by $S=M_D^{D-2}\int \sqrt{-g}R/(16\pi )$, the same algebraic volume relation uses the unreduced $M_{\mathrm{planck}}=1.221\times 10^{19}\mathrm{GeV}$. Holding that convention's $M_{\ast }=1\mathrm{TeV}$ gives radii larger by $(8\pi {)}^{1/n}$: approximately $4.69\times 10^{12}\mathrm{m}$, $6.80\times 10^{-4}\mathrm{m}$, $2.54\times 10^{-14}\mathrm{m}$. These are different conventions for what is held fixed at one TeV, not different scaling laws.

Finally, compactify the fifth coordinate on a circle of circumference $L=2\pi R$. For a source and observer at the same compact coordinate, the method of images turns a five-dimensional static potential $C/{\rho }^2$ intoOne can derive this image sum without assuming a summation formula. The Fourier transform of $1/(r^2+y^2)$ is $\pi e^{-r|k|}/r$; Poisson summation then gives $\pi [1+2{\sum }_{\ell \ge 1}{e}^{-2\pi \ell r/L}]/(Lr)$, which is the displayed hyperbolic cotangent. ThereforeThe short-distance limit is $C/r^2$; the long-distance limit is the four-dimensional $1/r$ potential. These exponential corrections are the nonzero Kaluza-Klein modes, with masses $2\pi |\ell |/L$. For electrostatics $C=q/(4{\pi }^2{\epsilon }_4)$, so the effective permittivity is ${\epsilon }_3=L{\epsilon }_4$. For a Poisson-normalized gravitational potential, $C=-{\kappa }_{4}m/(4{\pi }^2)$ gives ${\kappa }_3^{\mathrm{eff}}={\kappa }_4/L$. A massless scalar associated with the circle radius also contributes in an unstabilized gravitational compactification; recovering pure four-dimensional Einstein gravity with $G_4=G_5/L$ requires its stabilization or removal. The $1/r$ crossover itself is independent of this tensor-versus-scalar normalization issue.

Independent metric and Hamiltonian parametrization. The Lorentzian Polyakov action isThe worldsheet metric is independent of the embedding. Its equation sets the traceless worldsheet stress tensor to zero, making the metric locally conformal to the induced worldsheet metric; substitution recovers the Nambu–Goto action. Eliminating momentum in the Nambu-Goto phase-space action instead gives $(\dot{X}-u{X}^{\prime }{)}^2/(2e)-e{T}^2{X}^{\prime 2}/2$. On the positive-lapse branch, choose the independent metric, up to a positive Weyl transformation, asIts determinant is $-e^2{T}^2$. Multiplication by $-T/2$ recovers the reduced density, including the mixed term $-u\dot{X}\cdot X^{\prime }/e$. This proves the phase-space identification of the Polyakov metric.

Conformal gauge and residual coordinate freedom. In the Hamiltonian formulation, conformal gauge is $e=1/T,u=0$. The equations become $\ddot{X}-X^{\prime \prime }=0$, while ${\dot{X}}^2+X^{\prime 2}=0$ and $\dot{X}\cdot X^{\prime }=0$ remain as constraints. In the metric formulation, conformal gauge is ${\gamma }_{ab}=e^{2\omega }{\eta }_{ab}$, or simply ${\eta }_{ab}$ after fixing Weyl invariance. The metric parametrization proves these choices equivalent. Gauge-fixing the metric does not discard the Virasoro constraints.

With ${\sigma }^{\pm }=t\pm \sigma$, independent reparameterizations ${\sigma }^{\pm }\mapsto f_{\pm }({\sigma }^{\pm })$ multiply the flat metric by a conformal factor. A compensating Weyl transformation restores its chosen representative. Infinitesimally the Conformal Killing equation is ${\partial }_{-}{c}^{+}=0,{\partial }_{+}{c}^{-}=0$. These are the residual conformal transformations of two-dimensional Minkowski spacetime, with the periodicity conditions appropriate to a closed string.

Ghost determinant and anomaly. For gauge conditions $G_A(\varphi )=0$, the Faddeev-Popov determinant is the determinant of the linearized variationInsert this determinant in the gauge-fixed path integral and exponentiate it with anticommuting ghosts and antighosts: the additional action is proportional to $\int b^A{M}_{AB}{c}^B$. In metric gauge fixing, the trace is removed by Weyl invariance, and a worldsheet diffeomorphism changes the trace-free metric throughConsequently $c^a$ is a vector ghost and $b^{ab}$ a symmetric trace-free antighost. A conventional normalization of the worldsheet ghost action isA rescaling of the antighost changes only this overall kinetic normalization. In conformal gauge, after normalizing the chiral fields, the kinetic terms areThe closed-string ghosts are periodic, with independent left and right systems. Their conformal weights are $2$ for $b$ and $-1$ for $c$, so the supplied bc ghost system expression gives $c_{\mathrm{gh}}=-2(6\cdot 2^2-6\cdot 2+1)=-26$. Each free embedding boson has central charge one per chirality. HenceThis central charge of reparameterization ghosts cancels the string conformal anomaly in each chirality separately. Having two sectors does not double the required target dimension.

Channel variables, poles and the vector state. With particles 1,2 incoming and 3,4 outgoing, the Mandelstam variables are $s=-(p_1+p_2{)}^2$ and $t=-(p_1-p_3{)}^2$. They measure squared center-of-mass energy and momentum transfer. In the specified units ${\alpha }^{\prime }=1$, the external tachyon mass squared is $-1$, so $s+t+u=-4$.

The Gamma function residue at $-n$ is $(-1{)}^n/n!$. At $s=n-1$, its argument $-1-s$ varies with the opposite sign to $s$. The Gamma function recurrence givesThus the Veneziano amplitude pole residues areFor generic fixed $t$ the simple poles are $-1,0,1,\dots$. A zero of the residue product at an exceptional negative integer $t$ makes the corresponding pole removable. Fixed $t$ must be away from its own channel poles; fixing it at a channel singularity does not give an ordinary finite meromorphic function of $s$.

At $s=0$, the residue is $-(t+2)=-(t-u)/2$. This is proportional to $(p_1-p_2)\cdot (p_3-p_4)$, the contraction of conserved scalar-pair vector currents. Each current is orthogonal to the exchanged momentum because the external masses in its pair are equal. Alternatively, with scattering-angle cosine $z$, one has $t=-2(s/4+1)(1-z)$; at this pole $t+2=2z$, a pure spin-one angular dependence. ThereforeThe massless vector pole of the Veneziano amplitude is the massless open-string vector seen in the two-tachyon exchange channel.

For $\mathrm{Re}z>0$, the gamma function isRotate the integration ray through $\pi /2$, or insert a damping factor and then let it vanish. With the principal power, contour rotation givesFor the damping argument, $\int t^{\gamma -1}{e}^{-(ϵ-i)t}dt=\Gamma (\gamma )(ϵ-i{)}^{-\gamma }$, and Dirichlet tail estimates justify the limit. The condition $\gamma >0$ ensures integrability at zero; $\gamma <1$ makes the amplitude tend to zero and permits conditional convergence at infinity. At $\gamma =1$, the upper-end oscillation has no limit; larger $\gamma$ also fails ordinary improper convergence. These restrictions concern this improper integral, not its analytic regularizations.

Start from $I_0(\lambda )=1/\lambda$. Differentiation under the integral sign givesThis interchange is justified locally around any $\lambda >0$ by the dominated convergence theorem: the differentiated integrands are bounded by an integrable polynomial times $e^{-\lambda x/2}$. Successive derivatives of ${\lambda }^{-1}$ giveThis is also the integer case of the Gamma function integral after scaling the integration variable.

Put $c=e^{-\gamma }$, where $\gamma$ is the Euler--Mascheroni constant, and let $f(u)=w(u)-c$. The Buchstab function equation becomesThe given Gamma function bound implies $uf(u)\to 0$. Indeed,so $u/\Gamma (u+1)\to 0$.

Suppose first that $f\ge 0$ throughout an interval $[u,u+1]$, with $u\ge 1$. The delay differential equation and integration successively on $[u+1,u+2]$, $[u+2,u+3]$, and so on propagate nonnegativity to the entire tail $[u,\infty )$. In particular, $vf(v)$ is nonnegative and nondecreasing for $v\ge u+1$. Its limit is zero, so it is identically zero there. Thus $f=0$ on that tail.

The delay differential equation then propagates this equality backwards: if $f=0$ for $v\ge a$, then $f(v-1)=0$ for $v>\max (a,2)$. Repeating finitely many times would give $f=0$ on $(1,2)$, contradicting $f(v)=1/v-c$ there. Hence every interval $[u,u+1]$ contains a point where $f<0$.

If instead $f\le 0$ on $[u,u+1]$, the same propagation makes $vf(v)$ nonpositive and nonincreasing on the tail. Its limit zero again forces it to vanish, giving the same contradiction. Thus every such interval also contains a point where $f>0$. The oscillation of the Buchstab function is therefore strict on every unit interval: