We report the observation of radar echoes from the ionization trails of high-energy particle cascades. These data were taken at the SLAC National Accelerator Laboratory, where the full electron beam ($\sim$10$^9$ e$^-$ at $\sim$10 GeV/e$^-$) was directed into a plastic target to simulate an ultra high-energy neutrino interaction. This target was interrogated with radio waves, and coherent radio reflections from the cascades were detected, with properties consistent with theoretical expectations. This is the first definitive observation of radar echoes from high-energy particle cascades, which may lead to a viable neutrino detection technology for energies $\gtrsim 10^{16}$ eV.

Testbeam experiment 576 (T576) at the SLAC National Accelerator Laboratory sought to make the first measurement of coherent radio reflections from the ionization produced in the wake of a high-energy particle shower. The >10 GeV electron beam at SLAC End Station A was directed into a large high-density polyethylene target to produce a shower analogous to that produced by an EeV neutrino interaction in ice. Continuous wave (CW) radio was transmitted into the target, and receiving antennas monitored for reflection of the transmitted signal from the ionization left in the wake of the shower. We detail the experiment and report on preliminary hints of a signal consistent with a radio reflection at a significance of 2.36$\sigma$. We recommend another test-beam measurement in order to verify the signal.

We compare the statistics of parity even and odd multipoles of the cosmic microwave background (CMB) sky from PLANCK full mission temperature measurements. An excess power in odd multipoles compared to even multipoles has previously been found on large angular scales. Motivated by this apparent parity asymmetry, we evaluate directional statistics associated with even compared to odd multipoles, along with their significances. Primary tools are the \emph{Power Tensor} and \emph{Alignment Tensor} statistics. We limit our analysis to the first sixty multipoles i.e., $l=[2,61]$. We find no evidence for statistically unusual alignments of even parity multipoles. More than one independent statistic finds evidence for alignments of anisotropy axes of odd multipoles, with a significance equivalent to $\sim 2 \sigma$ or more. The robustness of alignment axes is tested by making galactic cuts and varying the multipole range. Very interestingly, the region spanned by the (a)symmetry axes is found to broadly contain other parity (a)symmetry axes previously observed in the literature.

Type 1a supernova magnitudes are used to fit cosmological parameters under the assumption the model will fit the observed redshift dependence. We test this assumption with the Union 2.1 compilation of 580 sources. Several independent tests find the existing model fails to account for a significant correlation of supernova color and redshift. The correlation of magnitude residuals relative to the $\Lambda CDM$ model and $color \times redshift$ has a significance equivalent to 13 standard deviations, as evaluated by randomly shuffling the data. Extending the existing $B-V$ color correction to a relation linear in redshift improves the goodness of fit $\chi^{2}$ by more than 50 units, an equivalent 7-$\sigma$ significance, while adding only one parameter. The $color-redshift$ correlation is quite robust, cannot be attributed to outliers, and passes several tests of consistency. We review previous hints of redshift dependence in color parameters found in bin-by-bin fits interpreted as parameter bias. We show that neither the bias nor the change $\Delta \chi^{2}$ of our study can be explained by those effects. The previously known relation that bluer supernovae have larger absolute luminosity tends to empirically flatten out with increasing redshift. The best-fit cosmological dark energy density parameter is revised from $ \Omega_{\Lambda} =0.71 \pm 0.02$ to $ \Omega_{\Lambda} = 0.74 \pm 0.02$ assuming a flat universe. One possible physical interpretation is that supernovae or their environments evolve significantly with increasing redshift.

The RICE experiment seeks observation of ultra-high energy (UHE; >100 PeV) neutrinos interacting in Antarctic ice, by measurement of the radiowavelength Cherenkov radiation resulting from the collision of a neutrino with an ice molecule. RICE was initiated in 1999 as a first-generation prototype for an eventual, large-scale in-ice UHE neutrino detector. Herein, we present updated limits on the diffuse UHE neutrino flux, based on twelve years of data taken between 1999 and 2010. We find no convincing neutrino candidates, resulting in 95% confidence-level model-dependent limits on the flux which improve on our previous (2006) results by a factor of two. We also report on ice properties' studies, including observation of radiowave reflections from layers internal to the ice sheet. We observe correlations of these reflections with the local ice flow direction, and also qualitatively verify the expected reduction in radio frequency ice transparency with increasing depth/temperature. A new experimental initiative (the Askaryan Radio Array, ARA) commenced at South Pole in 2010-11 will improve the RICE sensitivity (per unit time) by approximately two orders of magnitude.

The "cosmological principle" was set up early without realizing its implications for the horizon problem, and almost entirely without support from observational data. Consistent signals of anisotropy have been found in data on electromagnetic propagation, polarizations of QSOs and $CMB$ temperature maps. The axis of Virgo is found again and again in signals breaking isotropy, from independent observables in independent energy regimes. There are no satisfactory explanations of these effects in conventional astrophysics. Axion-photon mixing and propagation in axion condensates are capable of encompassing the data.

We analyze the effect of foregrounds on the observed alignment of CMBR quadrupole and octopole. The alignment between these multipoles is studied by using a symmetry based approach which assigns a principal eigenvector (PEV) or an axis with each multipole. We determine the significance of alignment between these multipoles by using the Internal Linear Combination (ILC) 5 and 7 year map s and also the maps obtained by using the Internal Power Spectrum Estimation (IPSE) procedure. The effect of foreground cleaning is studied in detail within the framework of the IPSE method both analytically and numerically. By using simulated CMBR data, we study how the PEVs of the pure simulated CMB map differ from those of the final cleaned map. We find that, in general, the shift in the PEVs is relatively small and in random directions. Due to the random nature of the shift we conclude that it can only lead to misalignment rather than alignment of multipoles. We also directly estimate the significance of alignment by using simulated cleaned maps. We find that the results in this case are identical to those obtained by simple analytic estimate or by using simulated pure CMB maps.

Many experiments seek dark matter by detecting relatively low energy nuclear recoils. Yet since events from ordinary physics with energies in the 1-100 KeV range are commonplace, all claims of signals or their absence hinge on exhaustive calibrations and background rejection. We document many curious and consistent discrepancies between the backgrounds which neutrons can produce versus the picture of neutrons and claims of neutron calibration found in dark matter literature. Much of the actual physics of neutrons is either under-recognized or under-reported, opening up new interpretations of current data. All signals seen so far, including those presented tentatively such as CoGENT, or the bold claims and time dependence of DAMA/LIBRA, appear to be consistent with neutron-induced backgrounds. At the same time it is the burden of proof of experimental groups to support their claims no possible background could matter, not ours. The existing hypotheses about backgrounds stated by experiments, accepted at face value and as published, make possible a variety of neutron-induced events to be registered as dark matter signals.

We find new and universal relations for the properties of dark matter particles consistent with standard relic abundances. Analysis is based on first characterizing the $s$-channel resonant annihilation process in great detail, keeping track of all velocity-dependence, the presence of multiple scales and treating each physical regime above, below, and close to thresholds separately. The resonant regime as well as extension to include non-resonant processes are then reduced to analytic formulas and inequalities that describe the full range of multi-dimensional numerical work. These results eliminate the need to recompute relic abundance model by model, and reduce calculations to verifying certain scale and parameter combinations are consistent. Remarkably simple formulas describe the relation between the total width of an $s$-channel intermediate particle, the masses and the couplings involved. Eliminating the width in terms of the mass produces new consistency relations between dark matter masses and the intermediate masses. The formulas are general enough to test directly whether new particles can be identified as dark matter. Resonance mass and total width are quantities directly observable at accelerators such as the LHC, and will be sufficient to establish whether new discoveries are consistent with the cosmological bounds on dark matter.

We find model-independent upper limits on rates of dark matter annihilation in galactic halos. The Born approximation generally fails, while exotic threshold enhancements akin to "Sommerfeld factors" also turn out to be baseless The most efficient annihilation mechanism involves perturbatively small decay widths that have largely been ignored. Widths that are very small compared to TeV mass scales suffice to effectively saturate unitarity bounds on annihilation rates. Bound state formation in weakly coupled theories produces small effects due to wave function normalizations. Unitarity shows the Sommerfeld factor cannot produce large changes in cross sections, and serves to identify where those approximations break down.

The space-time correlations of streams of photons can provide fundamentally new channels of information about the Universe. Today's astronomical observations essentially measure certain amplitude coherence functions produced by a source. The spatial correlations of wave fields has traditionally been exploited in Michelson-style amplitude interferometry. However the technology of the past was largely incapable of fine timing resolution and recording multiple beams. When time and space correlations are combined it is possible to achieve spectacular measurements that are impossible by any other means. Stellar intensity interferometry is ripe for development and is one of the few unexploited mechanisms to obtain potentially revolutionary new information in astronomy. As we discuss below, the modern use of stellar intensity interferometry can yield unprecedented measures of stellar diameters, binary stars, distance measures including Cepheids, rapidly rotating stars, pulsating stars, and short-time scale fluctuations that have never been measured before.

Recently a symmetry-based method to test for statistical isotropy of the cosmic microwave background was developed. We apply the method to template-cleaned 3-year and 5-year WMAP-$DA$ maps. We examine a wide range of angular multipoles from $2 < l < 300$. The analysis detects statistically signicant signals of anisotropy inconsistent with an isotropic CMB in some of the foreground cleaned maps. We are unable to resolve whether the anomalies have a cosmological, local astrophysical or instrumental origin. Assuming the anisotropy arises due to residual foreground contamination, we estimate the residual foreground power in the maps. For the W band maps, we also find a highly improbable degree of isotropy we cannot explain. We speculate that excess isotropy may be caused by faulty modeling of detector noise.

We report an upper limit on the flux of relativistic monopoles based on the non-observation of in-ice showers by the Radio Ice Cherenkov Experiment (RICE) at the South Pole. We obtain a 95% C.L. limit of order 10^{-18}/(cm^2-s-sr) for intermediate mass monopoles of 10^7<gamma<10^{12} at the anticipated energy E=10^{16} GeV. This bound is over an order of magnitude stronger than all previously published experimental limits for this range of boost parameters gamma, and exceeds two orders of magnitude improvement over most of the range. We review the physics of radio detection, describe a Monte Carlo simulation including continuous and stochastic energy losses, and compare to previous experimental limits.

We model the effect of photon and ultra-light pseudoscalar mixing on the propagation of electromagnetic radiation through the extragalactic medium. The medium is modelled as a large number of magnetic domains, uncorrelated with one another. We obtain an analytic expression for the different Stokes parameters in the limit of small mixing angle. The different Stokes parameters are found to increase linearly with the number of domains. We also verify this result by direct numerical simulations. We use this formalism to estimate the effect of pseudoscalar-photon mixing on the Cosmic Microwave Background (CMB) polarization. We impose limits on the model parameters by the CMB observations. We find that the currently allowed parameter range admits a CMB circular polarization up to order $10^{-7}$.

We introduce new symmetry-based methods to test for isotropy in cosmic microwave background radiation. Each angular multipole is factored into unique products of power eigenvectors, related multipoles and singular values that provide 2 new rotationally invariant measures mode by mode. The power entropy and directional entropy are new tests of randomness that are independent of the usual CMB power. Simulated galactic plane contamination is readily identified, and the new procedures mesh perfectly with linear transformations employed for windowed-sky analysis. The ILC -WMAP data maps show 7 axes well aligned with one another and the direction Virgo. Parameter free statistics find 12 independent cases of extraordinary axial alignment, low power entropy, or both having 5% probability or lower in an isotropic distribution. Isotropy of the ILC maps is ruled out to confidence levels of better than 99.9%, whether or not coincidences with other puzzles coming from the Virgo axis are included. Our work shows that anisotropy is not confined to the low l region, but extends over a much larger l range.

Radiation propagating over cosmological distances can probe light weakly interacting pseudoscalar (or scalar) particles. The existence of a spin-0 field changes the dynamical symmetries of electrodynamics. It predicts spontaneous generation of polarization of electromagnetic waves due to mode mixing in the presence of background magnetic field.We illustrate this by calculations of propagation in a uniform medium, as well as in a slowly varying background medium, and finally with resonant mixing. Highly complicated correlations between different Stokes parameters are predicted depending on the parameter regimes. The polarization of propagating waves shows interesting and complex dependence on frequency, the distance of propagation, coupling constants, and parameters of the background medium such as the plasma density and the magnetic field strength. For the first time we study the resonant mixing of electromagnetic waves with the scalar field, which occurs when the background plasma frequency becomes equal to the mass of the scalar field at some point along the path. Dynamical effects are found to be considerably enhanced in this case. We also formulate the condition under which the adiabatic approximation can be used consistently, and find caveats about comparing different frequency regimes.

We examine the Type 1a supernova data in order to determine if it shows any signal of large scale anisotropy. The anisotropy is modelled by an extended G\"{o}del metric, which incorporates expansion along with rotation. The model is smoothly connected to the usual FRW type, while expressing anisotropic metric effects depending on certain parameters. We find no significant signal of anisotropy in the data. We obtain bounds on an anisotropic redshift versus magnitude relationship, and accompanying parameters of the G\"{o}del-Obukhov metric.

We revisit the classification of polarization observables of the cosmic microwave background. There exists a unified approach to the $3 \times 3$ density matrix by which intensity, linear and circular polarization are treated on an equal footing. The representations of the rotation group on the sphere contain certain right-left setup symmetries which have not been taken into account before. Left-right symmetries revise the construction of invariants and certain predictions based on symmetries. There are two true rotational invariants mode by mode in polarization data. Independent of models of cosmological perturbations, we emphasize methods to construct invariant distributions in order to test predictions of isotropy, axial or parity symmetries, and generate covariant and invariant statistics order by order.

We show that distances of objects at cosmological distances can be measured directly using interferometry. Our approach to interferometric parallax comes from analysis of 4-point amplitude and intensity correlations that can be generated from pairs of well-separated detectors. The baseline required to measure cosmological distances of Gigaparsec order are within the reach of the next generation of space-borne detectors. The semi-classical interpretation of intensity correlations uses a notion of a single photon taking two paths simultaneously. Semi-classically a single photon can simultaneously enter four detectors separated by an astronomical unit, developing correlations feasible to measure with current technology.

Type 1a supernova magnitudes conventionally include an additive parameter called the extinction coefficient. We find that the extinction coefficients of a popular ``gold'' set are well correlated with the deviation of magnitudes from Hubble diagrams. If the effect is due to bias, extinctions have been overestimated, which makes supernovas appear more dim. The statistical significance of the extinction-acceleration correlation has a random chance probability of less than one in a million. The hypothesis that extinction coefficients should be corrected empirically provides greatly improved fits to both accelerating and non-accelerating models, with the independent feature of eliminating any significant correlation of residuals.

We reconsider analysis of data on the cosmic microwave background on the largest angular scales. Temperature multipoles of any order factor naturally into a direct product of axial quantities and cosets. Striking coincidences exist among the axes associated with the dipole, quadrupole, and octupole CMB moments. These axes also coincide well with two other axes independently determined from polarizations at radio and optical frequencies propagating on cosmological scales. The five coincident axes indicate physical correlation and anisotropic properties of the cosmic medium not predicted by the conventional Big Bang scenario. We consider various mechanisms, including foreground corrections, as candidates for the observed correlations. We also consider whether the propagation anomalies may be a signal of ``dark energy'' in the form of a condensed background field. Perhaps {\it light propagation} will prove to be an effective way to look for the effects of {\it dark energy}.

Radiation traversing the observable universe provides powerful ways to probe anisotropy of electromagnetic propagation. A controversial recent study claimed a signal of dipole character. Here we test a new and independent data set of 361 points under the null proposal of {\it statistical independence} of linear polarization alignments relative to galaxy axes, versus angular positions. The null hypothesis is tested via maximum likelihood analysis of best fits among numerous independent types of factored distributions. We also examine single-number correlations which are parameter free, invariant under coordinate transformations, and distributed very robustly. The statistics are shown explicitly not to depend on the uneven distribution of sources on the sky. We find that the null proposal is not supported at the level of less than 5% to less than 0.1% by several independent statistics. The signal of correlation violates parity, that is, symmetry under spatial inversion, and requires a statistic which transforms properly. The data indicate an axis of correlation, on the basis of likelihood determined to be $[{\rm R.A.}=(0^{\rm h},9^{\rm m}) \pm (1^{\rm h},0^{\rm m})$, ${\rm Decl.} = -1^o\pm 15^o]$.

The distribution of arrival directions of ultra-high energy cosmic rays may yield clues to their mysterious origin. We introduce a method of {\it invariant statistics} to analyze cosmic ray data which eliminates coordinate-dependent artifacts. When combined with maximum likelihood analysis, the method is capable of quantifying deviations of the distribution from isotropy with high reliability. We test our method against published AGASA events with energies above $4 \times 10^{19}$ eV. Angular cuts from observational limitations are taken into account. A model based on the Fisher distribution reveals the rotation of the Earth with the axis $\hat n$ along the direction ($5^h 53.36^m, 85.78^{\circ}$) in $(RA,DEC)$ coordinates, which is within $5^{\circ}$ of the equatorial north pole. Global anisotropy of the data, if any, hinges on finer understanding of detector acceptance than what is available from the published literature.

We reevaluate our published calculations of electromagnetic showers generated by GEANT 3.21 and the radio frequency pulses they produce in ice. We are prompted by a recent report showing that GEANT 3.21-modeled showers are sensitive to internal settings in the electron tracking subroutine. We report the shower and pulse characteristics obtained with different settings of GEANT 3.21 and with GEANT 4. The default setting of electron tracking in GEANT 3.21 we used in previous work speeds up the shower simulation at the cost of information near the end of the tracks. We find that settings tracking electron and positron to lower energy yield a more accurate calculation, a more intense shower, and proportionately stronger radio pulses at low frequencies. At high frequencies the relation between shower tracking algorithm and pulse spectrum is more complex. We obtain radial distributions of shower particles and phase distributions of pulses from 100 GeV showers that are consistent with our published results.

Cosmic rays in the regime of 10^19-10^20 eV should not exist, according to the GZK bound. Conflicts between data and the bound have led to suggestions that laws of physics must be modified. We find that physical processes neglected in the early history terminate the LPM effect, important in calibrating the energy of electromagnetic and hadronic showers in the GZK energy regime. The processes neglected are the direct losses of electrons and positrons, overlooked by LPM and greatly augmented by ultra-high energy circumstances. We present an exact formula for direct energy losses of e^+e^- pairs in electromagnetic showers. Numerical integration in current hadronic models shows that actual energy losses are always substantially larger than LPM estimates, but the energy at which crossover occurs depends on models.

Cubic kilometer neutrino telescopes are capable of probing fundamental questions of ultra-high energy neutrino interactions. There is currently great interest in neutrino interactions caused by low-scale, extra dimension models. Above 1 PeV the cross section in low scale gravity models rises well above the total Standard Model cross section. We assess the observability of this effect in the 1 PeV - 100 PeV energy range of kilometer-scale detectors with several new points of emphasis that hinge on enhanced neutral current cross sections. A major point is the importance of ``feed-down'' regeneration of upward neutrino flux, driven by new-physics neutral current interactions in the flux evolution equations. Feed-down is far from negligible, and it is essential to include its effect. We then find that the angular distribution of events has high discriminating value in separating models. In particular the ``up-to-down'' ratio between upward and downward-moving neutrino fluxes is a practical diagnostic tool which can discriminate between models in the near future. The slope of the angular distribution, in the region of maximum detected flux, is also substantially different in low-scale gravity and the Standard Model. These observables are only weakly dependent on astrophysical flux uncertainties. We conclude that angular distributions can reveal a breakdown of the Standard Model and probe the new physics beyond, as soon as data become available.

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.

In models based on large extra dimensions where massive spin 2 exchange can dominate at high energies, the neutrino-proton cross section can rise to typical hadronic values at energies above 10^20 eV. The neutrino then becomes a candidate for the primary that initiates the highest energy cosmic ray showers. We investigate characteristics of neutrino-induced showers compared to proton-induced showers. The comparison includes study of starting depth, profile with depth, lateral particle distribution at ground and muon lateral distribution at ground level. We find that for cross sections above 20 mb there are regions of parameter space where the two types of showers are essentially indistinguishable. We conclude that the neutrino candidate hypothesis cannot be ruled out on the basis of shower characteristics.

Angular correlations of ultra-high energy cosmic rays with cosmologically distant sources may provide clues to these mysterious events. We compare cosmic ray tracks with energies above $ 10^{20}eV$ to a compilation of radio-loud compact QSO positions. The statistical method emphasizes invariant quantities and a test of {\it statistical independence} of track and source distributions. Statistical independence is ruled out by several independent statistics at confidence levels of less that $10^{-3}$ ($99.9% $.), indicating that track directions and QSO source positions are correlated at a highly significant level.

Cosmic ray events above $10^{20}$ eV are on the verge of confronting fundamental particle physics. The neutrino is the only candidate primary among established particles capable of crossing 100 Mpc intergalactic distances unimpeded. The magnitude of $\nu N$ cross sections indicated by events, plus consistency with the Standard Model at low-energy, point to new physics of massive spin-2 exchange. In models based on extra dimensions, we find that the $\nu N$ cross section rises to typical hadronic values of between 1 and 100 mb at energies above $10^{20}$ eV. Our calculations take into account constraints of unitarity. We conclude that air-showers observed with energies above $10^{19}$ eV are consistent with neutrino primaries and extra-dimension models. An {\it upper bound} of 1-10 TeV on the mass scale at which graviton exchange becomes strong in current Kaluza-Klein models follows.

Radio Cherenkov emission underlines detection of high energy particles via a signal growing like the particle-energy-squared. Cosmic ray-induced electromagnetic showers are a primary application. While many studies have treated the phenomenon approximately, none have attempted to incorporate all the physical scales involved in problems with time- or spatially- evolving charges. We find it is possible to decompose the calculated fields into the product of a form factor, characterizing a moving charge distribution, multiplying a general integral which depends on the charge evolution. In circumstances of interest for cosmic ray physics, the resulting expressions can be evaluated explicitely in terms of a few parameters obtainable from shower codes. The classic issues of Frauhofer and Fresnel zones play a crucial role in the coherence.

We report on the feasibility of using an isotropic flux of cosmic neutrinos in the energy range of 10 to 10,000 TeV to study the interior geophysical structure of the Earth. The angular distribution of events in a $\sim {\rm km}^3$-scale neutrino telescope can be inverted to yield the Earth's nucleon density distribution. The approach is independent of previous seismic methods. The energy spectrum of the neutrino primaries can also be determined from consistency with the angular distribution. Depending on neutrino flux, a reasonable density determination can be obtained with a few hundred receivers operating for a few years.

We study the feasibility of using an isotropic flux of cosmic neutrinos in the energy range of 10 to 10000 TeV to study the interior structure of Earth. The angular distribution of events in a $\sim {\rm km}^3$-scale neutrino telescope can be inverted to yield information on the Earth's mass distribution that is independent of other methods. The energy spectrum of the neutrino primaries is also determined from consistency with the angular distribution. It is possible to make a model independent determination of the density profile of Earth's interior, separate from the absolute normalization of the incident cosmic neutrinos.

Polarizations of electromagnetic waves from distant galaxies are known to be correlated with the source orientations. These quantities have been used to search for signals of cosmological birefringence. We review and classify transformation properties of the polarization and source orientation observables. The classifications give a firm foundation to certain practices which have sprung up informally in the literature. Transformations under parity play a central role, showing that parity violation in emission or in the subsequent propagation is an observable phenomenon. We also discuss statistical measures, correlations and distributions which transform properly and which can be used for systematic data analysis.

RICE is designed to detect ultrahigh energy (>100 TeV) neutrinos from astrophysical sources. It will consist of an array of compact radio (100 to 1000 MHz) receivers buried in ice at the South Pole. During the 1995-96 and 1996-97 austral summers, several receivers transmitters were deployed in bore holes drilled for the AMANDA project, at depths of 141 to 260 m. This was the first in situ test of radio receivers in deep ice for neutrino astronomy.

We discuss a new mechanism which could cause a rotation of polarization of electromagnetic waves due to magnetic fields on cosmological scales. The effect is due to the geometrical phase of Pancharatnam and Berry, and causes a corkscrew twisting of the plane of polarization. The new effect represents an additional tool that allows possible intergalactic and cosmological magnetic fields to be studied using radio propagation.

We review evidence for a new phenomenon in the propagation of radio waves across the Universe, an anisotropic rotation of the plane of polarization not accounted for by conventional physics.

Recently, J. P. Leahy commented (astro-ph/9704285) on the indication of anisotropy in the propagation of radio waves over cosmological distances, published by Nodland and Ralston (B. Nodland and J. P. Ralston, Phys. Rev. Lett. 78, 3043 (1997); astro-ph/9704196). Unfortunately, Leahy ignores the need for a statistical comparison, egregiously selects data, compares incompatible measures, uses a variable incapable of resolving birefringence, and makes mistakes in the calculations he presents.

Eisenstein and Bunn have widely circulated a Comment (astro-ph/9704247) suggesting a non-standard null hypothesis in a test for cosmological birefringence. The Comment misrepresents the procedure used by Nodland and Ralston (B. Nodland and J. P. Ralston, Phys. Rev. Lett. 78, 3043 (1997); astro-ph/9704196) and lacks statistical basis; no calculations were reported, but sweeping conclusions were drawn from eyeballing a single scatter plot. The results of the suggested procedure range from an underestimate of the statistical significance of well-correlated data, to a failure of detecting a perfect correlation in the limit of strong cuts. We verified these faults by performing the actual calculation which EB suggested, but neglected to carry out themselves. Furthermore, the calculation showed that the original correlation remains statistically significant, with the lego-plot of 1/P versus s showing a persistent bump, and the anisotropy direction s remaining in the same direction as previously reported.

We present evidence for a new phenomenon in the propagation of electromagnetic radiation across the Universe, a corkscrew rotation of the plane of polarization not accounted for by conventional physics.

We report a systematic rotation of the plane of polarization of electromagnetic radiation propagating over cosmological distances. The effect is extracted independently from Faraday rotation, and found to be correlated with the angular positions and distances to the sources. Monte Carlo analysis yields probabilistic P-values of order 10^(-3) for this to occur as a fluctuation. A fit yields a birefringence scale of order 10^(25) meters. Dependence on redshift z rules out a local effect. Barring hidden systematic bias in the data, the correlation indicates a new cosmological effect.

The ultra-high energy neutrino nucleon cross section grows at a surprising rate with energy due to QCD effects in the target. Recent electroproduction data allows an update of earlier predictions. We compare the results of our own calculations with those of other groups, and critically review the foundation and reliability of the calculations. The question of the `new physics' potential of neutrino telescopes sensitive to the total cross section in multi-PeV energy domain is considered. We point out a loophole in the arguments which might be an important consideration in extrapolating the cross section to extremely high energies.

Interactions of ultrahigh energy neutrinos of cosmological origin in large volumes of radio-transparent South Polar ice can be detected via coherent Cherenkov emission from accompanying electromagnetic showers. A pilot experiment employing buried radio receivers has been successfully deployed at the South Pole and data are now being collected. The physics of coherent radio emission together with the properties of radio-pulse propagation in Antarctic ice clearly distinguishes the radio method from phototube detection. In the context of the proposed km^3-scale neutrino telescope, these two detection modes provide complementary information on UHE neutrino interactions.

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.

New structure function data at small Bjorken $x$ from HERA are used along with next-to-leading order QCD evolution to predict a cross section for charged-current interactions of ultrahigh energy neutrinos with nucleons. This new result is over twice the size of previous estimates and has important implications for cosmic ray experiments now underway as well as for KM3 arrays (cubic kilometer-scale neutrino telescopes) now in the planning stages.