Source author record

Loyal Durand

Loyal Durand appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

21works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

21 published item(s)

preprint2022arXiv

Fractional operators and multi-integral representations for associated Legendre functions

In a recent paper, Cohl and Costas-Santos derived a number of interesting multi-derivative and multi-integral relations for associated Legendre and Ferrers functions in which the orders of those functions are changed in integral steps. These are of potential use in a number of physical problems. We show here how their results can be derived simply from more general relations involving non-integer changes in the order obtained using the fractional group operator methods developed earlier for SO(2,1), E(2,1) and its conformal extension, and SO(3). We also present general integral relations for fractional changes of the degrees of the functions, and related multi-derivative and multi-integral representations.

preprint2016arXiv

Comment on "More on Heisenberg's model for high energy nucleon-nucleon scattering"

We comment on the treatment of asymptotic black-disk scattering in a recent paper of Nastase and Sonnenschein, Phys.\ Rev.\ D\ {\bf 92}, 015028 (2015), on scattering in an updated version of the Heisenberg model which gives $pp$ and $\bar{p}p$ cross sections which increase at very high energies as $\ln^2s$. We show that the total cross section they define does not correspond to that measured in experiments, with the result that their limit for the ratio $σ_{\rm elas}/σ_{\rm tot}$ is too small by a factor 2. The correct ratio for black-disk scattering, $σ_{\rm elas}/σ_{\rm tot} \rightarrow 1/2$ for $s\rightarrow\infty$, is strongly supported by experiment.

preprint2016arXiv

The slope, curvature, and higher parameters in $pp$ and $\bar{p}p$ scattering, and the extrapolation of measurements of $dσ(s,t)/dt$ to $t=0$

We study the effects of curvature in the expansion of the logarithm of the differential elastic scattering cross section near $t=0$ as $dσ(s,t)/dt=dσ(s,0)/dt\,\times\exp(Bt+Ct^2+Dt^3\cdots)$ in an eikonal model for $pp$ and $\bar{p}p$ scattering, and use the results to discuss the extrapolation of measured differential cross sections and the slope parameters $B$ to $t=-q^2=0$. We find that the curvature effects represented by the parameters $C$ and $D$, while small, lead to significant changes in the forward slope parameter relative to that determined in a purely exponential fit, and to smaller but still significant changes in the forward elastic scattering and total cross sections. Curvature effects should therefore be considered in future analyses or reanalyses of the elastic scattering data.

preprint2015arXiv

Comprehensive fits to high energy data for $σ$, $ρ$, and $B$ and the asymptotic black-disk limit

We demonstrate that the entirety of the data on proton--proton and antiproton--proton forward scattering between 6 GeV and 57 TeV center-of-mass energy is sufficient to show that $σ_{\rm elas}/σ_{\rm tot} \rightarrow 1/2$, and that $8πB/σ_{\rm tot}\rightarrow 1$ at very high energies, where $B$ the forward slope parameter for the differential elastic scattering cross sections. The relations demonstrate convincingly that the asymptotic $pp$ and $\bar{p}p$ scattering amplitudes approach those of scattering from a black disk. This result obviously has implications for any new physics that modifies the forward scattering amplitudes.

preprint2015arXiv

Eikonal fit to $pp$ and $\bar{p}p$ scattering and the edge in the scattering amplitude

We make a detailed eikonal fit to current data on the total and elastic scattering cross sections, the ratios $ρ$ of the real to the imaginary parts of the forward elastic scattering amplitudes, and the logarithmic slopes $B$ of the differential cross sections $dσ/dt$ at $t=0$, for proton-proton and antiproton-proton scattering at center-of-mass energies $W$ from 5 GeV to 57 TeV. The fit allows us to investigate the structure of the eikonal amplitudes in detail, including the impact-parameter structure of the energy-independent edge in the scattering amplitude shown to exist by Block {\em et al.} \cite{edge}. We show that the edge region has an essentially fixed shape with a peak at approximately the "black disk" radius $R_{\rm tot}=\sqrt{σ_{\rm tot}/2π}$ of the scattering amplitude, a constant width $t_{\rm edge}\approx 1$ fm, and migrates to larger impact parameters with increasing energy proportionally to $R_{\rm tot}$. We comment on possible physical mechanisms which could lead to the edge. We show that the eikonal results for the cross sections and $ρ$ values are described to high accuracy by analytic expressions of the forms used in earlier analyses by Block and Halzen, and extend the result to the elastic-scattering slope parameter $B$. These expressions provide simple extrapolations of the results to much higher energies. Finally, we calculate the survival probabilities for large rapidity gaps in the scattering.

preprint2015arXiv

Evidence for a Constant `Edge' in Proton-Proton Scattering at Very High Energies

Accurate fits to $pp$ and $\bar pp$ cross section data up to Tevatron energies, incorporating the constraints imposed by analyticity and unitarity, successfully predict the results of recent LHC and cosmic ray measurements, and suggest that the cross sections approach a black disc limit asymptotically. The approach to the limit is, however, very slow. We present a simple geometric picture which explains these features in a natural way. A black disc of logarithmically growing radius is supplemented by a soft `edge' whose properties are invariant with energy. The constancy of the edge results in the prediction that the quantity $(σ^{TOT}-2σ^{El})/\surdσ^{TOT}$ approaches a constant at high energy. Using the existing fits, this prediction appears to be verified. The value of the limiting constant allows an estimate of the thickness of the edge, which turns out to be on the order of $1\,{\rm fm}$. One thus arrives at a picture where the proton-proton scattering at lower energies is dominated by what becomes the edge, while at higher energies it is dominated by the disc. The crossover between the two regimes is only at $\surd s\geq $ 10 TeV, accounting for the slow approach to asymptotic behavior. Some questions as to the nature of the edge are discussed.

preprint2014arXiv

Connection of the virtual $γ^*p$ cross section of $ep$ deep inelastic scattering to real $γp$ scattering, and the implications for $νN$ and $ep$ total cross sections

We show that it is possible to fit all of the HERA DIS (deep inelastic scattering) data on $F_2^{γp}$ at small values of Bjorken $x$, including the data at {\em very low} $Q^2$, using a new model for $F_2^{γp}$ which both includes an asymptotic (high energy) part that satisfies a saturated Froissart bound behavior, with a vector-dominance like mass factor in the parameterization, and extends smoothly to $Q^2=0$. We require that the corresponding part of the virtual $γ^* p$ cross section match the known asymptotic part of the real $γp$ cross section at $Q^2=0$, a cross section which is determined by strong interactions and asymptotically satisfies a saturated Froissart bound of the form $α+β\ln s+γ\ln^2s$. Using this model for the asymptotic part of $F_2^{γp}$ plus a known valence contribution, we fit the asymptotic high energy part of the HERA data with $x\le 0.1$ and $W\ge 25$ GeV; the fit is excellent. We find that the mass parameter in the fit lies in the region of the light vector mesons, somewhat above the $ρ$ meson mass, and is compatible with vector dominance. We use this fit to obtain accurate results for the high energy $ep$ and isoscalar $νN$ total cross sections. Both cross sections obey an analytic expression of the type $a +b \ln E +c \ln^2 E +d \ln^3 E$ at large energies $E$ of the incident particle, reflecting the fact that the underlying strong interaction parts of the $γ^*p$, $Z^*N$ and $W^*N$ cross sections satisfy the saturated Froissart bound. Since approximately 50% of the $νN$ center of mass (cms) energy is found in $W$---the cms energy of the strongly interacting intermediate vector boson-nucleon system---a study of ultra-high-energy neutrino-nucleon cross sections would allow us, for the first time, to explore {\em strong interactions at incredibly high energies}.

preprint2013arXiv

Implications of a Froissart bound saturation of $γ^*$-$p$ deep inelastic scattering. Part I. Quark distributions at ultra small $x$

We argue that the deep inelastic structure function $F_2^{γp}(x, Q^2)$, regarded as a cross section for virtual $γ^*p$ scattering, is hadronic in nature. This implies that its growth is limited by the Froissart bound at high hadronic energies, giving a $\ln^2 (1/x)$ bound on $F_2^{γp}$ as Bjorken $x\rightarrow 0$. The same bound holds for the individual quark distributions. In earlier work, we obtained a very accurate global fit to the combined HERA data on $F_2^{γp}$ using a fit function which respects the Froissart bound at small $x$, and is equivalent in its $x$ dependence to the function used successfully to describe all high energy hadronic cross sections, including $γp$ scattering. We extrapolate that fit by a factor of $\lesssim$3 beyond the HERA region in the natural variable $\ln(1/x)$ to the values of $x$ down to $x=10^{-14}$ and use the results to derive the quark distributions needed for the reliable calculation of neutrino cross sections at energies up to $E_ν=10^{17}$ GeV. These distributions do not satisfy the Feynman "wee parton" assumption, that they all converge toward a common distribution $xq(x,Q^2)$ at small $x$ and large $Q^2$. This was used in some past calculations to express the dominant neutrino structure function $F_2^{ν(\barν)}$ directly in terms of $F_2^{γp}$. We show that the correct distributions nevertheless give results for $F_2^{ν(\barν)}$ which differ only slightly from those obtained assuming that the wee parton limit holds. In two Appendices, we develop simple analytic results for the effects of QCD evolution and operator-product corrections on the distribution functions at small $x$, and show that these effects amount mainly to shifting the values of $\ln(1/x)$ in the initial distributions.

preprint2013arXiv

Implications of a Froissart bound saturation of $γ^*$-$p$ deep inelastic scattering. Part II. Ultra-high energy neutrino interactions

In Part I (in this journal) we argued that the structure function $F_2^{γp}(x,Q^2)$ in deep inelastic $ep$ scattering, regarded as a cross section for virtual $γ^*p$ scattering, has a saturated Froissart-bounded form behaving as $\ln^2 (1/x)$ at small $x$. This form provides an excellent fit to the low $x$ HERA data, including the very low $Q^2$ regions, and can be extrapolated reliably to small $x$ using the natural variable $\ln(1/x)$. We used our fit to derive quark distributions for values of $x$ down to $x=10^{-14}$. We use those distributions here to evaluate ultra-high energy (UHE) cross sections for neutrino scattering on an isoscalar nucleon, $N=(n+p)/2$, up to laboratory neutrino energies $E_ν\sim 10^{16}$-$10^{17}$ GeV where there are now limits on neutrino fluxes. We estimate that these cross sections are accurate to $\sim$2% at the highest energies considered, with the major uncertainty coming from the errors in the parameters that were needed to fit $F_2^{γp}(x,Q^2)$. We compare our results to recently published neutrino cross sections derived from NLO parton distribution functions, which become much larger at high energies because of the use of power-law extrapolations of quark distributions to small $x$. We argue that our calculation of the UHE $νN$ cross sections is the best one can make based the existing experimental deep inelastic scattering data. Further, we show that the strong interaction Froissart bound of $\ln^2 (1/x)$ on $F_2^{γp}$ translates to an exact bound of $\ln^3E_ν$ for leading-order-weak $νN$ scattering. The energy dependence of $νN$ total cross section measurements consequently has important implications for hadronic interactions at enormous cms (center-of-mass) energies not otherwise accessible.

preprint2011arXiv

A new numerical method for inverse Laplace transforms used to obtain gluon distributions from the proton structure function

We recently derived a very accurate and fast new algorithm for numerically inverting the Laplace transforms needed to obtain gluon distributions from the proton structure function $F_2^{γp}(x,Q^2)$. We numerically inverted the function $g(s)$, $s$ being the variable in Laplace space, to $G(v)$, where $v$ is the variable in ordinary space. We have since discovered that the algorithm does not work if $g(s)\rightarrow 0$ less rapidly than $1/s$ as $s\rightarrow\infty$, e.g., as $1/s^β$ for $0<β<1$. In this note, we derive a new numerical algorithm for such cases, which holds for all positive and non-integer negative values of $β$. The new algorithm is {\em exact} if the original function $G(v)$ is given by the product of a power $v^{β-1}$ and a polynomial in $v$. We test the algorithm numerically for very small positive $β$, $β=10^{-6}$ obtaining numerical results that imitate the Dirac delta function $δ(v)$. We also devolve the published MSTW2008LO gluon distribution at virtuality $Q^2=5$ GeV$^2$ down to the lower virtuality $Q^2=1.69$ GeV$^2$. For devolution, $ β$ is negative, giving rise to inverse Laplace transforms that are distributions and not proper functions. This requires us to introduce the concept of Hadamard Finite Part integrals, which we discuss in detail.

preprint2011arXiv

An analytic solution to LO coupled DGLAP evolution equations: a new pQCD tool

We have analytically solved the LO pQCD singlet DGLAP equations using Laplace transform techniques. Newly-developed highly accurate numerical inverse Laplace transform algorithms allow us to write fully decoupled solutions for the singlet structure function F_s(x,Q^2)and G(x,Q^2) as F_s(x,Q^2)={\cal F}_s(F_{s0}(x), G_0(x)) and G(x,Q^2)={\cal G}(F_{s0}(x), G_0(x)). Here {\cal F}_s and \cal G are known functions of the initial boundary conditions F_{s0}(x) = F_s(x,Q_0^2) and G_{0}(x) = G(x,Q_0^2), i.e., the chosen starting functions at the virtuality Q_0^2. For both G and F_s, we are able to either devolve or evolve each separately and rapidly, with very high numerical accuracy, a computational fractional precision of O(10^{-9}). Armed with this powerful new tool in the pQCD arsenal, we compare our numerical results from the above equations with the published MSTW2008 and CTEQ6L LO gluon and singlet F_s distributions, starting from their initial values at Q_0^2=1 GeV^2 and 1.69 GeV^2, respectively, using their choices of α_s(Q^2). This allows an important independent check on the accuracies of their evolution codes and therefore the computational accuracies of their published parton distributions. Our method completely decouples the two LO distributions, at the same time guaranteeing that both G and F_s satisfy the singlet coupled DGLAP equations. It also allows one to easily obtain the effects of the starting functions on the evolved gluon and singlet structure functions, as functions of both Q^2 and Q_0^2, being equally accurate in devolution as in evolution. Further, it can also be used for non-singlet distributions, thus giving LO analytic solutions for individual quark and gluon distributions at a given x and Q^2, rather than the numerical solutions of the coupled integral-differential equations on a large, but fixed, two-dimensional grid that are currently available.

preprint2011arXiv

Applications of the leading-order Dokshitzer-Gribov-Lipatov-Altarelli-Parisi evolution equations to the combined HERA data on deep inelastic scattering

We recently derived explicit solutions of the leading-order Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations for the $Q^2$ evolution of the singlet structure function $F_s(x,Q^2)$ and the gluon distribution $G(x,Q^2)$ using very efficient Laplace transform techniques. We apply our results here to a study of the HERA data on deep inelastic $ep$ scattering as recently combined by the H1 and ZEUS groups. We use initial distributions $F_2^{γp}(x,Q_0^2)$ and $G(x,Q_0^2)$ fixed by a global fit to the HERA data. From $F_2^{γp}(x,Q_0^2)$ we obtain the singlet quark distribution $F_s(x,Q_0^2)$---using small non-singlet quark distributions taken from either the CTEQ6L or the MSTW2008LO analyses---evolve to arbitrary $Q^2$, and then convert the results to individual quark distributions. Finally, we show directly from a study of systematic trends in a comparison of the evolved $F_2^{γp}(x,Q^2)$ with the HERA data, that the assumption of leading-order DGLAP evolution is inconsistent with those data.

preprint2010arXiv

Decoupling the coupled DGLAP evolution equations: an analytic solution to pQCD

Using Laplace transform techniques, along with newly-developed accurate numerical inverse Laplace transform algorithms, we decouple the solutions for the singlet structure function $F_s(x,Q^2)$ and $G(x,Q^2)$ of the two leading-order coupled singlet DGLAP equations, allowing us to write fully decoupled solutions: F_s(x,Q^2)={\cal F}_s(F_{s0}(x), G_0(x)), G(x,Q^2)={\cal G}(F_{s0}(x), G_0(x)). Here ${\cal F}_s$ and $\cal G$ are known functions---found using the DGLAP splitting functions---of the functions $F_{s0}(x) \equiv F_s(x,Q_0^2)$ and $G_{0}(x) \equiv G(x,Q_0^2)$, the chosen starting functions at the virtuality $Q_0^2$. As a proof of method, we compare our numerical results from the above equations with the published MSTW LO gluon and singlet $F_s$ distributions, starting from their initial values at $Q_0^2=1 GeV^2$. Our method completely decouples the two LO distributions, at the same time guaranteeing that both distributions satisfy the singlet coupled DGLAP equations. It furnishes us with a new tool for readily obtaining the effects of the starting functions (independently) on the gluon and singlet structure functions, as functions of both $Q^2$ and $Q_0^2$. In addition, it can also be used for non-singlet distributions, thus allowing one to solve analytically for individual quark and gluon distributions values at a given $x$ and $Q^2$, with typical numerical accuracies of about 1 part in $10^5$, rather than having to evolve numerically coupled integral-differential equations on a two-dimensional grid in $x, Q^2$, as is currently done.

preprint2010arXiv

Decoupling the NLO coupled DGLAP evolution equations: an analytic solution to pQCD

Using repeated Laplace transform techniques, along with newly-developed accurate numerical inverse Laplace transform algorithms, we transform the coupled, integral-differential NLO singlet DGLAP equations first into coupled differential equations, then into coupled algebraic equations, which we can solve iteratively. After Laplace inverting the algebraic solution analytically, we numerically invert the solutions of the decoupled differential equations. Finally, we arrive at the decoupled NLO evolved solutions F_s(x,Q^2)=calF_s(F_{s0}(x),G_0(x)) and G(x,Q^2)=calG(F_{s0}(x),G_0(x)), where calF_s and calG are known functions - determined using the DGLAP splitting functions up to NLO in the strong coupling constant alpha_s(Q^2). The functions F_{s0}(x)=F_s(x,Q_0^2) and G_0(x)=G(x,Q_0^2) are the starting functions for the evolution at Q_0^2. This approach furnishes us with a new tool for readily obtaining, independently, the effects of the starting functions on either the evolved gluon or singlet structure functions, as a function of both Q^2 and Q_0^2. It is not necessary to evolve coupled integral-differential equations numerically on a two-dimensional grid, as is currently done. The same approach can be used for NLO non-singlet distributions where it is simpler, only requiring one Laplace transform. We make successful NLO numerical comparisons to two non-singlet distributions, using NLO quark distributions published by the MSTW collaboration, over a large range of x and Q^2. Our method is readily generalized to higher orders in the strong coupling constant alpha_s(Q^2).

preprint2009arXiv

Analytic treatment of leading-order parton evolution equations: theory and tests

We recently derived an explicit expression for the gluon distribution function G(x, Q^2) = xg(x, Q^2) in terms of the proton structure function F_2^{γp} (x, Q^2) in leading-order (LO) QCD by solving the the LO DGLAP equation for the Q^2 evolution of F_2^{γp} (x, Q^2) analytically, using a differential-equation method. We re-derive and extend the results here using a Laplace-transform technique, and show that the singlet quark structure function F_S(x,Q^2) can be determined directly in terms of G from the DGLAP gluon evolution equation. To illustrate the method and check the consistency of existing LO quark and gluon distributions, we used the published values of the LO quark distributions from the CTEQ5L and MRST2001LO analyses to form F_2^{γp} (x, Q^2), and then solved analytically for G(x,Q^2). We find that the analytic and fitted gluon distributions from MRST2001LO agree well with each other for all x and Q^2, while those from CTEQ5L differ significantly from each other for large x values, x>~0.03 - 0.05 at all Q^2. We conclude that the published CTEQ5L distributions are incompatible in this region. Using a non-singlet evolution equation, we obtain a sensitive test of quark distributions which holds in both LO and NLO perturbative QCD. We find in either case that the CTEQ5 quark distributions satisfy the tests numerically for small x, but fail the tests for x>~0.03 - 0.05 - their use could potentially lead to significant shifts in predictions of quantities sensitive to large x. We encountered no problems with the MRST2001LO distributions or later CTEQ distributions. We suggest caution in the use of the CTEQ5 distributions.

preprint2008arXiv

Analytic derivation of the leading-order gluon distribution function G(x,Q^2) = xg(x,Q^2) from the proton structure function F_2^p(x,Q^2)

We derive a second-order linear differential equation for the leading order gluon distribution function G(x,Q^2) = xg(x,Q^2) which determines G(x,Q^2) directly from the proton structure function F_2^p(x,Q^2). This equation is derived from the leading order DGLAP evolution equation for F_2^p(x,Q^2), and does not require knowledge of either the individual quark distributions or the gluon evolution equation. Given an analytic expression that successfully reproduces the known experimental data for F_2^p(x,Q^2) in a domain x_min<=x<=x_max, Q_min^2<=Q^2<=Q_max^2 of the Bjorken variable x and the virtuality Q^2 in deep inelastic scattering, G(x,Q^2) is uniquely determined in the same domain. We give the general solution and illustrate the method using the recently proposed Froissart bound type parametrization of F_2^p(x,Q^2) of E. L. Berger, M. M. Block and C-I. Tan, PRL 98, 242001, (2007). Existing leading-order gluon distributions based on power-law description of individual parton distributions agree roughly with the new distributions for x>~10^-3 as they should, but are much larger for x<~10^-3.

preprint2005arXiv

Electromagnetic corrections to baryon masses

We analyze the electromagnetic contributions to the octet and decuplet baryon masses using the heavy baryon approximation in chiral effective field theory and methods we developed in earlier analyses of the baryon masses and magnetic moments. Our methods connect simply to Morpurgo's general parametrization of the electromagnetic contributions and to semirelativistic quark models. Our calculations are carried out including the one-loop mesonic corrections to the basic electromagnetic interactions, so to two loops overall. We find that to this order in the chiral loop expansion there are no three-body contributions. The Coleman-Glashow relation and other sum rules derived in quark models with only two-body terms therefore continue to hold, and violations involve at least three-loop processes and can be expected to be quite small. We present the complete formal results and some estimates of the matrix elements here. Numerical calculations will be presented separately.

preprint2002arXiv

Fractional operators and special functions. II. Legendre functions

Most of the special functions of mathematical physics are connected with the representation of Lie groups. The action of elements $D$ of the associated Lie algebras as linear differential operators gives relations among the functions in a class, for example, their differential recurrence relations. In this paper, we apply the fractional generalizations $D^μ$ of these operators developed in an earlier paper in the context of Lie theory to the group SO(2,1) and its conformal extension. The fractional relations give a variety of interesting relations for the associated Legendre functions. We show that the two-variable fractional operator relations lead directly to integral relations among the Legendre functions and to one- and two-variable integral representations for those functions. Some of the relations reduce to known fractional integrals for the Legendre functions when reduced to one variable. The results enlarge the understanding of many properties of the associated Legendre functions on the basis of the underlying group structure.

preprint2001arXiv

Effective field theory and the quark model, II. Structure of loop corrections

We analyze the structure of meson loop corrections to the O($m_s$) expressions for the baryon masses and magnetic moments in heavy-baryon chiral perturbation theory (HBChPT), and show in detail how the bulk of the corrections are absorbed into redefinitions of the unknown parameters of the O($m_s$) chiral expansion. To effect this analysis, we use the three-flavor-index representation of the effective baryon fields and HBChPT developed in a preceding paper, and a decomposition of the corrections in terms of effective one-, two-, and three-quark operators. The results show why the loop corrections have so little effect on fits to the masses and moments, and do not seriously disrupt the Gell-Mann--Okubo relations for the masses and the Okubo relation for the moments even though individual loops can be quite large. We also examine the momentum structure of the residual loop corrections, and comment on limits on their validity in HBChPT. The structural analysis can be generalized straightforwardly to other problems in HBChPT using the three-flavor-index representation of the effective baryon fields, and provides a fairly simple way to determine what parts of the loop corrections are actually significant in a given setting.

preprint1994arXiv

Two-loop ${\rm O}\left(G_F^2M_H^4\right)$ corrections to the fermionic decay rates of the Higgs boson

We calculate the dominant ${\rm O}\left(G_F^2M_H^4\right)$ two-loop electroweak corrections to the fermi\-onic decay widths of a heavy Higgs boson in the Standard Model. Use of the Goldstone-boson equivalence theorem reduces the problem to one involving only the physical Higgs boson $H$ and the Goldstone bosons $w^\pm$ and $z$ of the unbroken theory. The two-loop corrections are opposite in sign to the one-loop electroweak corrections, exceed the one-loop corrections in magnitude for $M_H>1114\ {\rm GeV}$, and increase in relative magnitude as $M_H^2$ for larger values of $M_H$. We conclude that the perturbation expansion in powers of $G_FM_H^2$ breaks down for $M_H\approx 1100\ {\rm GeV}$. We discuss briefly the QCD and the complete one-loop electroweak corrections to $H\rightarrow b\bar{b}, \,t\bar{t}$, and comment on the validity of the equivalence theorem. Finally we note how a very heavy Higgs boson could be described in a phenomenological manner.

preprint1993arXiv

Two-loop renormalization constants and high energy $2\rightarrow 2$ scattering amplitudes in the Higgs sector of the standard model

We calculate the complete matrix of two-body scattering amplitudes for the scattering of longitudinally polarized gauge bosons $W_L^\pm$, $Z_L$ and Higgs bosons to two loops in the high-energy, heavy-Higgs limit $\sqrt{s}\gg M_H\gg M_W$. Use of the Goldstone boson equivalence theorem reduces the problem to one involving only the scalar fields $w^\pm$, $z$ (the Goldstone bosons of the original theory) and the Higgs boson. Renormalization of the scattering amplitudes requires the calculation of the self-energy functions $Π_i^0(M_i^2)$, the renormalization constants $Z_i$, and the bare quartic Higgs coupling $λ_0$ to two loops. The results will be useful in other calculations. To facilitate the calculations, we introduce a powerful new technique for evaluating integrals over Feynman parameters in dimensional regularization which is based on a Barnes' type representation of the binomial expansion. We also collect some useful integrals which extend the tables given by Devoto and Duke.