Source author record

Douglas W. McKay

Douglas W. McKay 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

25works
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

25 published item(s)

preprint2016arXiv

Light Sterile Neutrinos, Lepton Number Violating Interactions and the LSND Anomaly

We develop the consequences of introducing a purely leptonic, lepton number violating non-standard interaction (NSI) and standard model neutrino mixing with a fourth, sterile neutrino in the analysis of short-baseline, neutrino experiments. We focus on the muon decay at rest (DAR) result from the Liquid Scintillation Neutrino Experiment (LSND) and the Karlsruhe and Rutherford Medium Energy Neutrino Experiment (KARMEN). We make a comprehensive analysis of lepton number violating, NSI effective operators and find nine that affect muon decay relevant to LSND results. Two of these preserve the standard model (SM) value 3/4 for the Michel rho and delta parameters and, overall, show favorable agreement with precision data and the electron anti-neutrino signal from LSND data. We display theoretical models that lead to these two effective operators. In the model we choose to apply to DAR data, both electron anti-neutrino appearance from muon anti-neutrino oscillation and electron anti-neutrino survival after production from NSI decay of the positive muon contribute to the expected signal. This is a unique feature of our scheme. We find a range of parameters where both experiments can be accommodated consistently with recent global, sterile neutrino fits to short baseline data. We comment on implications of the models for new physics searches at colliders and comment on further implications of the lepton number violating interactions plus sterile neutrino-standard model neutrino mixing.

preprint2014arXiv

Optimizing the Determination of the Neutrino Mixing Angle $θ_{13}$ from Reactor Data

The technical breakthroughs of multiple detectors developed by Daya Bay and RENO collaborations have gotten great attention. Yet the optimal determination of neutrino mixing parameters from reactor data depends on the statistical method and demands equal attention. We find that a straightforward method using a minimal parameters will generally outperform a multi-parameter method by delivering more reliable values with sharper resolution. We review standard confidence levels and statistical penalties for models using extra parameters, and apply those rules to our analysis. We find that the methods used in recent work of the Daya Bay and RENO collaborations have several undesirable properties. The existing work also uses non-standard measures of significance which we are unable to explain. A central element of the current methods consists of variationally fitting many more parameters than data points. As a result the experimental resolution of $\sin ^{2}(2θ_{13})$ is degraded. The results also become extremely sensitive to certain model parameters that can be adjusted arbitrarily. The number of parameters to include in evaluating significance is an important issue that has generally been overlooked. The measures of significance applied previously would be consistent if and only if all parameters but one were considered to have no physical relevance for the experiment's hypothesis test. Simpler, more transparent methods can improve the determination of the mixing angle $θ_{13}$ from reactor data, and exploit the advantages from superb hardware technique of the experiments. We anticipate that future experimental analysis will fully exploit those advantages.

preprint2014arXiv

Short Baseline Reactor Antineutrino-Electron Scattering Experiments and Non-Standard Neutrino Interactions at Source and Detector

We investigate non-standard interaction effects in antineutrino-electron scattering experiments with baselines short enough to ignore standard oscillation phenomena. The setup is free of ambiguities from the interference between new physics and oscillation effects and is sensitive to both semileptonic new physics at the source and purely leptonic new physics in the weak interaction scattering at the detector. We draw on the TEXONO experiment as the model system, extending its analysis of non-standard interaction effects at the detector to include the generally allowed non-standard interaction phase at the detector and both non-universal and flavor changing new physics at the reactor source. We confirm that the current data allows for new physics constraints at the detector of the same order as those currently published, but we find that constraints on the source new physics are at least an order of magnitude weaker. The new physics phase effects are at the 5% level, noticeable in the 90% C.L. contour plots but not significantly affecting the conclusions. Based on projected increase in sensitivity with an upgraded TEXONO experiment, we estimate the improvement of sensitivity to both source and detector non-standard interactions. We find that the bounds on source parameters improve by an order of magnitude, but do not reach parameter space beyond current limits. On the other hand, the detector new physics sensitivity would push current limits by factors 5 to 10 smaller.

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.

preprint2013arXiv

Sensitivity of medium-baseline reactor neutrino mass-hierarchy experiments to nonstandard interactions

We explore the impact of nonstandard interactions at source and detector on the interpretation of reactor electron neutrino disappearance experiments with short- and medium-baseline designs. We use the constraints from the recent results from short-baseline experiments and generalize current estimates of medium-baseline event rates to include charged current interactions at source and detector with standard Lorentz structure but with nonstandard flavor structure. We find that the average spectrum of observed events at a baseline of 50 km, in the middle of the currently favored region, provides a probe of new interactions. We show that an improvement in sensitivity to nonstandard interactions is possible if combined with improved precision of input mixing parameters in independent experiments, despite ambiguity in interpretation of medium-baseline data. We show that nonstandard interactions can enhance or suppress the sensitivity of experiments to the mass hierarchy, depending on the phases of the parameters and the CP-violating phase in the standard three-neutrino mixing picture.

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.

preprint2011arXiv

Comment on "Ultrahigh-Energy Neutrino-Nucleon Deep-Inelastic Scattering and the Froissart Bound": Phys. Rev. Lett. 106, 231802 (2011)

The authors of a recent paper, "Ultrahigh-Energy Neutrino-Nucleon Deep-Inelastic Scattering and the Froissart Bound", A. Illarianov, B. Kniehl and A. Kotikov, Phys. Rev. Lett. 106, 231802 (2011), derive an approximate formula for the UHE limit of $σ_{νN}(s)$ in a class of models that includes our own and assert that they are led "to the important observation that $σ_{BBT}^{νN} \propto ln^3s$, which manifestly violates the Froissart bound [2] in contrast to what is stated in Refs. [6-8]", the latter reference being to our work and the $σ_{BBT}^{νN}$ to the cross sections we reported there. We here correct their erroneous implication that $σ_{νN}(s) should satisfy the Froissart bound and their mistaken assertion that we state that $σ_{BBT}^{νN}$ satisfies it.

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).

preprint2010arXiv

New physics, the cosmic ray spectrum knee, and $pp$ cross section measurements

We explore the possibility that a new physics interaction can provide an explanation for the knee just above $10^6$ GeV in the cosmic ray spectrum. We model the new physics modifications to the total proton-proton cross section with an incoherent term that allows for missing energy above the scale of new physics. We add the constraint that the new physics must also be consistent with published $pp$ cross section measurements, using cosmic ray observations, an order of magnitude and more above the knee. We find that the rise in cross section required at energies above the knee is radical. The increase in cross section suggests that it may be more appropriate to treat the scattering process in the black disc limit at such high energies. In this case there may be no clean separation between the standard model and new physics contributions to the total cross section. We model the missing energy in this limit and find a good fit to the Tibet III cosmic ray flux data. We comment on testing the new physics proposal for the cosmic ray knee at the Large Hadron Collider.

preprint2010arXiv

Ultrahigh energy neutrino scattering: an update

We update our estimates of charged and neutral current neutrino total cross sections on isoscalar nucleons at ultrahigh energies using a global (x, Q^2) fit, motivated by the Froissart bound, to the F_2 (electron-proton) structure function utilizing the most recent analysis of the complete ZEUS and H1 data sets from HERA I. Using the large Q^2, small Bjorken-x limits of the "wee" parton model, we connect the ultrahigh energy neutrino cross sections directly to the large Q^2, small-x extrapolation of our new fit, which we assume saturates the Froissart bound. We compare both to our previous work, which utilized only the smaller ZEUS data set, as well as to recent results of a calculation using the ZEUS-S based global perturbative QCD parton distributions using the combined HERA I results as input. Our new results substantiate our previous conclusions, again predicting significantly smaller cross sections than those predicted by extrapolating pQCD calculations to neutrino energies above 10^9 GeV.

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.

preprint2006arXiv

Bounds on Low Scale Gravity from RICE data and Cosmogenic Neutrino Flux Models

We explore limits on low scale gravity models set by results from the Radio Ice Cherenkov Experiment's (RICE) ongoing search for cosmic ray neutrinos in the cosmogenic, or GZK, energy range. The bound on M_D, the fundamental scale of gravity, depends upon cosmogenic flux model, black hole formation and decay treatments, inclusion of graviton mediated elastic neutrino processes, and the number of large extra dimensions, d. Assuming proton-based cosmogenic flux models that cover a broad range of flux possibilities, we find bounds in the interval 0.9 TeV <M_D< 10 TeV. Heavy nucleus-based models generally lead to smaller fluxes and correspondingly weaker bounds. Values d = 5, 6 and 7, for which laboratory and astrophysical bounds on LSG models are less restrictive, lead to essentially the same limits on M_D.

preprint2002arXiv

Coherent Radio Pulses From GEANT Generated Electromagnetic Showers In Ice

Radio Cherenkov radiation is arguably the most efficient mechanism for detecting showers from ultra-high energy particles of 1 PeV and above. Showers occuring in Antarctic ice should be detectable at distances up to 1 km. We report on electromagnetic shower development in ice using a GEANT Monte Carlo simulation. We have studied energy deposition by shower particles and determined shower parameters for several different media, finding agreement with published results where available. We also report on radio pulse emission from the charged particles in the shower, focusing on coherent emission at the Cherenkov angle. Previous work has focused on frequencies in the 100 MHz to 1 GHz range. Surprisingly, we find that the coherence regime extends up to tens of Ghz. This may have substantial impact on future radio-based neutrino detection experiments as well as any test beam experiment which seeks to measure coherent Cherenkov radiation from an electromagnetic shower. Our study is particularly important for the RICE experiment at the South Pole.

preprint1998arXiv

Limits on Exotic Quarks in the SU(3)xU(1) Extension of the Standard Model from SUSY Search Data

We study the $p\bar p$ production and decay of exotic quarks that are predicted by the SU(3)$\times$U(1) extension of the Standard Model. We show that recent experimental searches for SUSY particles at the Tevatron limit the mass of these quarks to be above 250 GeV. Run II will extend the reach to 320 GeV. This is one example of SUSY search signatures that apply directly to another, quite different, new physics model.

preprint1997arXiv

Single Top from Technipion Production

We investigate the contribution of technicolor mechanisms to the production of single top quarks at hadron colliders. Technipions with mass larger than the top mass will decay predominantly to a top quark plus a bottom antiquark. We investigate two promising sub-processes: technipion plus W-boson via gluon-gluon fusion and technipion plus quark production via quark gluon interaction. The decay chain of technipion to top plus bottom quarks and then top to W plus bottom yields final states for the two subprocesses with, respectively, two W's and two bottom quarks and one W, two bottom quarks and a light quark. We calculate the total cross sections and the $p_{T}$ distributions for these technipion production mechanisms at Tevatron and LHC energies for a range of technipion masses, starting at 200 GeV. We study the backgrounds to our processes and the kinematic cuts that enhance the signal to background ratio and we report event rate estimates for the upgraded Tevatron and the LHC. Only the LHC has the potential to observe these processes.

preprint1996arXiv

Study of Quark Propagator Solutions to the Dyson--Schwinger Equation in a Confining Model

We solve the Dyson--Schwinger equation for the quark propagator in a model with singular infrared behavior for the gluon propagator. We require that the solutions, easily found in configuration space, be tempered distributions and thus have Fourier transforms. This severely limits the boundary conditions that the solutions may satisify. The sign of the dimensionful parameter that characterizes the model gluon propagator can be either positive or negative. If the sign is negative, we find a unique solution. It is singular at the origin in momentum space, falls off like $1/p^2$ as $p^2\rightarrow +/-\infty$, and it is truly nonperturbative in that it is singular in the limit that the gluon--quark interaction approaches zero. If the sign of the gluon propagator coefficient is positive, we find solutions that are, in a sense that we exhibit, unconstrained linear combinations of advanced and retarded propagators. These solutions are singular at the origin in momentum space, fall off like $1/p^2$ asympotically, exhibit ``resonant--like" behavior at the position of the bare mass of the quark when the mass is large compared to the dimensionful interaction parameter in the gluon propagator model, and smoothly approach a linear combination of free--quark, advanced and retarded two--point functions in the limit that the interaction approaches zero. In this sense, these solutions behave in an increasingly ``particle--like" manner as the quark becomes heavy. The Feynman propagator and the Wightman function are not tempered distributions and therefore are not acceptable solutions to the Schwinger--Dyson equation in our model. On this basis we advance several arguments to show that the Fourier--transformable solutions we find are consistent with quark confinement, even though they have singularities on the

preprint1995arXiv

A Minimal Four-Family Supergravity Model

We investigate the phenomenology of minimal four-family MSSM supergravity theories containing an additional generation of heavy fermions along with their superpartners. We demand: gauge coupling constant unification at high energy scales; perturbative values for all Yukawa couplings for energy scales $\leq \mgut$; radiative electroweak symmetry breaking via renormalization group evolution; and a neutral LSP. The perturbative constraints imply a light fourth-family quark and lepton spectrum, and $\tanb\lsim 3$. The lightest CP-even Higgs mass is increased. Fourth-family Yukawa coupling contributions to the evolution of scalar masses lead to unexpected mass hierarchies; \eg\ the $\staupone$ is generally the lightest slepton and the lightest squark is the $\wt\bpr_1$. A significant lower bound is placed on the gluino mass by the requirement that the $\staupone$ not be the LSP. Sleptons of the first two families are much more massive compared to the LSP and other neutralinos and charginos than in the three-family models, and could easily lie beyond the reach of a $\sqrt s=500\gev$ $\epem$ collider. Relations between slepton masses and gaugino masses are shown to be very sensitive to the presence of a fourth generation. The most important near-future experimental probes of the four-family models are reviewed. A scenario with $\mt\sim\mw$ and $t\rta {\wt t_1}\cnone$ is shown to be inconsistent with universal soft-SUSY-breaking boundary conditions. Full four-family evolution of $α_s$ is shown to lead to a significant enhancement in inclusive jet and di-jet spectra at Tevatron energies when all sparticle masses are near their lower bounds.

preprint1995arXiv

On Radio Detection of Ultra-High Energy Neutrinos in Antarctic Ice

Interactions of ultrahigh energy neutrinos of cosmological origin in large volumes of dense, radio-transparent media can be detected via coherent Cherenkov emission from accompanying electromagnetic showers. Antarctic ice meets the requirements for an efficient detection medium for a radio frequency neutrino telescope. We carefully estimate the sensitivity of realistic antennas embedded deep in the ice to 100 MHz - 1 GHz signals generated by predicted neutrino fluxes from active galactic nuclei. Our main conclusion is that a {\it single radio receiver} can probe a $\sim 1$ ${\rm km}^3$ volume for events with primary energy near 2 PeV and that the total number of events registered would be roughly 200 to 400 ${\rm year}^{-1}$ in our most conservative estimate. An array of such receivers would increase sensitivity dramatically. A radio neutrino telescope could directly observe and test our understanding of the most powerful particle accelerators in the universe, simultaneously testing the standard theory of particle physics at unprecedented energies.

preprint1995arXiv

Possible Higgs Boson Effects on the Running of Third and Fourth Generation Quark Masses and Mixings.

The Schwinger-Dyson equation for the quark self-energy is solved for the case of the third and fourth quark generations. The exchanges of standard model gluons and Higgs bosons are taken into account. It is found that Higgs boson exchange dominates quark self-energy in the ultraviolet region for sufficiently large input quark masses, causing the running quark propagator mass to increase with energy-scale. No running of the quark mixing angles is found for input quark masses up to and including 500 GeV.

preprint1994arXiv

Gauge-Coupling Unification and the Minimal SUSY Model: A Fourth Generation Below the Top?

We explore the possibility of a fourth generation in the gauge-coupling-unified, minimal supersymmetric (MSSM) framework. We find that a sequential fourth generation (with a heavy neutrino $ν'$) can still fit, surviving all present experimental constraints, provided $λ_b(M_U)=λ_τ(M_U)$ Yukawa unification is relaxed. For the theory to remain perturbative up to $M_U$, the new leptonic generation must lie within reach of LEP-II and the new $b',t'$ must have masses within the reach of the Tevatron. For example, for $\mt>150\gev$ we find $m_{ν'},m_{τ'}< 86\gev$, $m_{t'}<178$, and $m_{b'}<156\gev$. Experiments at Fermilab are already sensitive to the latter mass regions; we comment on direct $b'$ searches and on the $m_{t'}\simeq \mt$ case in light of new CDF data. Discovery may involve novel decay signatures; however, CDF and LEP-II will confirm or exclude an MSSM fourth generation in the near future.

preprint1993arXiv

Isospin Multiplet Structure in Ultra--Heavy Fermion Bound States

The coupled Bethe--Salpeter bound state equations for a $Q\bar Q$ system, where $Q=(U,D)$ is a degenerate, fourth generation, super--heavy quark doublet, are solved in several ladder approximation models. The exchanges of gluon, Higgs and Goldstone modes in the standard model are calculated in the ultra--heavy quark limit where weak $γ, W^\pm$ and $Z^0$ contributions are negligible. A natural $I=0$ and $I=1$ multiplet pattern is found, with large splittings occuring between the different weak iso--spin states when $M_Q$, the quark masses, are larger than values in the range $0.4 TeV<M_Q<0.8 TeV$, depending on which model is used. Consideration of ultra--heavy quark lifetime constraints and $U-D$ mass splitting constraints are reviewed to establish the plausibility of lifetime and mass degeneracy requirements assumed for this paper.