We propose in this White Paper a concept for a space experiment using cold atoms to search for ultra-light dark matter, and to detect gravitational waves in the frequency range between the most sensitive ranges of LISA and the terrestrial LIGO/Virgo/KAGRA/INDIGO experiments. This interdisciplinary experiment, called Atomic Experiment for Dark Matter and Gravity Exploration (AEDGE), will also complement other planned searches for dark matter, and exploit synergies with other gravitational wave detectors. We give examples of the extended range of sensitivity to ultra-light dark matter offered by AEDGE, and how its gravitational-wave measurements could explore the assembly of super-massive black holes, first-order phase transitions in the early universe and cosmic strings. AEDGE will be based upon technologies now being developed for terrestrial experiments using cold atoms, and will benefit from the space experience obtained with, e.g., LISA and cold atom experiments in microgravity. This paper is based on a submission (v1) in response to the Call for White Papers for the Voyage 2050 long-term plan in the ESA Science Programme. ESA limited the number of White Paper authors to 30. However, in this version (v2) we have welcomed as supporting authors participants in the Workshop on Atomic Experiments for Dark Matter and Gravity Exploration held at CERN: ({\tt https://indico.cern.ch/event/830432/}), as well as other interested scientists, and have incorporated additional material.

Since the very beginning of astronomy the location of objects on the sky has been a fundamental observational quantity that has been taken for granted. While precise two dimensional positional information is easy to obtain for observations in the electromagnetic spectrum, the positional accuracy of current and near future gravitational wave detectors is limited to between tens and hundreds of square degrees, which makes it extremely challenging to identify the host galaxies of gravitational wave events or to confidently detect any electromagnetic counterparts. Gravitational wave observations provide information on source properties and distances that is complementary to the information in any associated electromagnetic emission and that is very hard to obtain in any other way. Observing systems with multiple messengers thus has scientific potential much greater than the sum of its parts. A gravitational wave detector with higher angular resolution would significantly increase the prospects for finding the hosts of gravitational wave sources and triggering a multi-messenger follow-up campaign. An observatory with arcminute precision or better could be realised within the Voyage 2050 programme by creating a large baseline interferometer array in space and would have transformative scientific potential. Precise positional information of standard sirens would enable precision measurements of cosmological parameters and offer new insights on structure formation; a high angular resolution gravitational wave observatory would allow the detection of a stochastic background and resolution of the anisotropies within it; it would also allow the study of accretion processes around black holes; and it would have tremendous potential for tests of modified gravity and the discovery of physics beyond the Standard Model.

Detailed observation by the High-Altitude Water Cerenkov Observatory of the regions within about 20 pc of the nearby Geminga and PSR B0656+14 pulsars have suggested the need for the diffusion coefficient of 10-100 TeV leptons yielding the detected gamma-rays to be about 4x10^27 cm^2 s^-1, two orders of magnitude less than the value usually required to explain galactic cosmic ray propagation at these energies. This result casts doubt on the ability of pulsars to explain the positron excess above the value determined by secondary cosmic ray production, although use of two zone diffusion models mitigates the problem. It is the purpose of this paper to give a simple demonstration of the possibility that the gradient in the pulsar accelerated electrons/positrons can give rise to the Alfven instability whereby the particles are self-confined. An estimated value of the Halo diffusion coefficient of 3x10^27 cm^2 s^-1 at 47 TeV is found.

We discuss two extraordinary increases of cosmic ray intensity observed by Voyager 1 in the last 1.1 AU before heliopause crossing, Aug 2012 at 121.7 AU. The two increases are roughly similar in amplitude and result in a total increase of cosmic ray nuclei around 1 GV of over 50 percent and of 0.01 GV electrons of a factor about 2. During the first increase, the changes in the magnetic, B, field are small. After the first increase, the B field variations and cosmic ray changes become large and during the second increase the B field variations and the cosmic ray changes are correlated to within a day. The intensity variations of H and He nuclei during these time intervals are measured from 0.1 to over 1 GV. The total GCR increse over the two events resemble those expected from a simple force-field 'like' solar modulation model with a modulation potential of about 80MV. This is nearly one third of the total modulation potential of about 250 GV required to produce the modulation of these particles observed at the earth at the 2009 sunspot minimum and adds a new aspect to ideas about heliospheric modulation.

Voyager 1 has explored the solar wind-interstellar medium interaction region between the terminal shock and heliopause following the intensity distribution of galactic cosmic ray protons above 200 MeV energy. Before this component reached the galactic level at 121.7 AU, 4 episodes of rapid intensity change occured similar to the Forbush Decreases found near the sun, rather than the expected result of models related to those describing Long Term Modulation in the inner solar system. Because the mean solar wind flow is both expected and observed to be perpendicular to the radial direction close to the heliopause, explanation is given in terms of transient radial flows related to possible heliopause boundary flapping. It is necessary that radial flows are at the sound speed found for conditions downstream of the teminal shock and that the relevant perpendicular cosmic ray diffusion is controlled by 'slab' field fluctuations accounting for 20 percent or less of the total power in field variance. However, additional radial drift motion related to possible north to south gradients in the magnetic field may allow the inclusion of some diffusion according to 2-D turbulence theory. The required field gradients may arise due to variation in the field carried by the solar plasma deflected away from the solar equatorial plane. Modulation amounting to a total 30 percent drop in galactic intensity requires explanation by a combination of several transient episodes.

Voyager 1 has explored the solar wind-interstellar medium interaction region between the Terminal Shock and the Heliopause, following the intensity distribution of the shock accelerated anomalous component of cosmic rays in the MeV energy range. The sudden disappearance of this component at 121.7 AU from the sun is discussed in terms of three models for the transition into the interstellar plasma flow. Particles trapped flowing parallel to the boundary may penetrate up to one Larmour radius beyond. If the boundary is stationary, Voyager 1 directly samples this distance. The boundary could flap, depending on Heliosheath pressure changes and Voyager 1 then samples the extent of this motion. Finally, a turbulent boundary layer is considered in which the MeV particle distribution falls off with distance, thus measuring diffusion within the layer.

We studied the application of statistical reconstruction algorithms, namely maximum likelihood and least squares methods, to the problem of event reconstruction in a dual phase liquid xenon detector. An iterative method was developed for in-situ reconstruction of the PMT light response functions from calibration data taken with an uncollimated gamma-ray source. Using the techniques described, the performance of the ZEPLIN-III dark matter detector was studied for 122 keV gamma-rays. For the inner part of the detector (R<100 mm), spatial resolutions of 13 mm and 1.6 mm FWHM were measured in the horizontal plane for primary and secondary scintillation, respectively. An energy resolution of 8.1% FWHM was achieved at that energy. The possibility of using this technique for improving performance and reducing cost of scintillation cameras for medical applications is currently under study.

We report experimental upper limits on WIMP-nucleon elastic scattering cross sections from the second science run of ZEPLIN-III at the Boulby Underground Laboratory. A raw fiducial exposure of 1,344 kg.days was accrued over 319 days of continuous operation between June 2010 and May 2011. A total of eight events was observed in the signal acceptance region in the nuclear recoil energy range 7-29 keV, which is compatible with background expectations. This allows the exclusion of the scalar cross-section above 4.8E-8 pb near 50 GeV/c^2 WIMP mass with 90% confidence. Combined with data from the first run, this result improves to 3.9E-8 pb. The corresponding WIMP-neutron spin-dependent cross-section limit is 8.0E-3 pb. The ZEPLIN programme reaches thus its conclusion at Boulby, having deployed and exploited successfully three liquid xenon experiments of increasing reach.

Plastic scintillators are widely used in industry, medicine and scientific research, including nuclear and particle physics. Although one of their most common applications is in neutron detection, experimental data on their response to low-energy nuclear recoils are scarce. Here, the relative scintillation efficiency for neutron-induced nuclear recoils in a polystyrene-based plastic scintillator (UPS-923A) is presented, exploring recoil energies between 125 keV and 850 keV. Monte Carlo simulations, incorporating light collection efficiency and energy resolution effects, are used to generate neutron scattering spectra which are matched to observed distributions of scintillation signals to parameterise the energy-dependent quenching factor. At energies above 300 keV the dependence is reasonably described using the semi-empirical formulation of Birks and a kB factor of (0.014+/-0.002) g/MeVcm^2 has been determined. Below that energy the measured quenching factor falls more steeply than predicted by the Birks formalism.

ZE3RA is the software package responsible for processing the raw data from the ZEPLIN-III dark matter experiment and its reduction into a set of parameters used in all subsequent analyses. The detector is a liquid xenon time projection chamber with scintillation and electroluminescence signals read out by an array of 31 photomultipliers. The dual range 62-channel data stream is optimised for the detection of scintillation pulses down to a single photoelectron and of ionisation signals as small as those produced by single electrons. We discuss in particular several strategies related to data filtering, pulse finding and pulse clustering which are tuned to recover the best electron/nuclear recoil discrimination near the detection threshold, where most dark matter elastic scattering signatures are expected. The software was designed assuming only minimal knowledge of the physics underlying the detection principle, allowing an unbiased analysis of the experimental results and easy extension to other detectors with similar requirements.

The ZEPLIN-III experiment is operating in its second phase at the Boulby Underground Laboratory in search of dark matter WIMPs. The major upgrades to the instrument over its first science run include lower background photomultiplier tubes and installation of a plastic scintillator veto system. Performance results from the veto detector using calibration and science data in its first six months of operation in coincidence with ZEPLIN-III are presented. With fully automated operation and calibration, the veto system has maintained high stability and achieves near unity live time relative to ZEPLIN-III. Calibrations with a neutron source demonstrate a rejection of 60% of neutron-induced nuclear recoils in ZEPLIN-III that might otherwise be misidentified as WIMPs. This tagging efficiency reduces the expected untagged nuclear recoil background from neutrons during science data taking to a very low rate of ~0.2 events per year in the WIMP acceptance region. Additionally, the veto detector provides rejection of 28% of gamma-ray induced background events, allowing the sampling of the dominant source of background in ZEPLIN-III - multiple scatter gamma-rays with rare topologies. Since WIMPs will not be tagged by the veto detector, and tags due to gamma-rays and neutrons are separable, this population of multiple scatter events may be characterised without biasing the analysis of candidate WIMP signals in the data.

WMAP's detection of the Sunyaev-Zel'dovich effect at a much reduced level among several large samples of rich clusters is interpreted in terms of conventional physics. It has been suggested that the central soft X-ray and EUV excess found in some clusters cannot be of thermal origin, due to problems with rapid gas cooling and the persistent non-detection of the O VII line, but may arise from inverse-Compton scattering between intracluster relativistic electrons and the cosmic microwave background. In fact, recent XMM observations of the soft X-rays from Coma and Abell 3112 are equally well fitted by a power law or a thermal virialized gas. Therefore the missing Sunyaev-Zel'dovich flux could partly be due to an overestimate of the central density of virialized electrons which scatter the CMB. Synchrotron radiation in an intracluster magnetic field of strength of a few $\mu$G is responsible for significant additional electron energy loss. Equipartition between relativistic particle and magnetic field energy densities is a realistic possibility. GHz radiation data from a Coma cluster halo yields information on the high energy steepening of the cluster relativistic electron spectrum. Cluster microwave emission in the WMAP passbands by higher energy cosmic ray electrons and gamma ray emission from an accompanying cosmic ray proton flux are also considered. The energetic electrons could originate from AGN jet injection, then distributed cluster-wide by Alfven wave sweeping, with accompanying {\it in situ} Fermi acceleration.

The design, optimisation and construction of an anti-coincidence veto detector to complement the ZEPLIN-III direct dark matter search instrument is described. One tonne of plastic scintillator is arranged into 52 bars individually read out by photomultipliers and coupled to a gadolinium-loaded passive polypropylene shield. Particular attention has been paid to radiological content. The overall aim has been to achieve a veto detector of low threshold and high efficiency without the creation of additional background in ZEPLIN-III, all at a reasonable cost. Extensive experimental measurements of the components have been made, including radioactivity levels and performance characteristics. These have been used to inform a complete end-to-end Monte Carlo simulation that has then been used to calculate the expected performance of the new instrument, both operating alone and as an anti-coincidence detector for ZEPLIN-III. The veto device will be capable of rejecting over 65% of coincident nuclear recoil events from neutron background in the energy range of interest in ZEPLIN-III. This will reduce the background in ZEPLIN-III from ~0.4 to ~0.14 events per year in the WIMP acceptance region, a significant factor in the event of a non-zero observation. Furthermore, in addition to providing valuable diagnostic capabilities, the veto is capable of tagging over 15% for gamma-ray rejection, all whilst contributing no significant additional background. In conjunction with the replacement of the internal ZEPLIN-III photomultiplier array, the new veto is expected to improve significantly the sensitivity of the ZEPLIN-III instrument to dark matter, allowing spin independent WIMP-nucleon cross sections below 1E-8 pb to be probed.

The ZEPLIN-III experiment in the Palmer Underground Laboratory at Boulby uses a 12kg two-phase xenon time projection chamber to search for the weakly interacting massive particles (WIMPs) that may account for the dark matter of our Galaxy. The detector measures both scintillation and ionisation produced by radiation interacting in the liquid to differentiate between the nuclear recoils expected from WIMPs and the electron recoil background signals down to ~10keV nuclear recoil energy. An analysis of 847kg.days of data acquired between February 27th 2008 and May 20th 2008 has excluded a WIMP-nucleon elastic scattering spin-independent cross-section above 8.1x10(-8)pb at 55GeV/c2 with a 90% confidence limit. It has also demonstrated that the two-phase xenon technique is capable of better discrimination between electron and nuclear recoils at low-energy than previously achieved by other xenon-based experiments.

Diffusive shock acceleration is invoked to explain non-thermal particle acceleration in Supernova Remnants, Active Galactic Nuclei (AGN) Jets, Gamma ray Bursts (GRBs) and various large scale cosmic structures. The importance of achieving the highest observed particle energies by such a mechanism in a given astrophysical situation is a recurring theme. In this work, shock acceleration in relativistic shocks is discussed, mostly focusing on a numerical study concerning proton acceleration efficiency by subluminal and superluminal shocks, emphasising on the dependence of the scattering model, bulk Lorentz factor and the angle between the magnetic field and the shock flow. We developed a diffuse cosmic ray model based on the study of different shock boost factors, which shows that spectra from AGN fit current observations of ultra high energy cosmic rays, above 5.7 x 10^10 GeV, much better than GRBs, indicating that AGN are the primary candidates to explain the UHECR flux. Recent Fermi observations of GRB090816c indicate very flat spectra which are expected within our model predictions and support evidence that GRB particle spectra can be flat, when the shock Lorentz factor is of order ~1000.

The flux of Ultra High Energy Cosmic Rays (UHECRs) at $E>10^{18.5}$ eV is believed to arise in plasma shock environments in extragalactic sources. In this paper, we present a systematic study of particle acceleration by relativistic shocks, in particular concerning the dependence on bulk Lorentz factor and the angle between the magnetic field and the shockflow. For the first time, simulation results of super- and subluminal shocks with boost factors up to $\Gamma=1000$ are investigated and compared systematically. While superluminal shocks are shown to be inefficient at the highest energies ($E>10^{18.5}$ eV), subluminal shocks may provide particles up to $10^{21}$ eV, limited only by the Hillas-criterion. For the subluminal case, we find that mildly relativistic shocks, thought to occur in jets of Active Galactic Nuclei (AGN, $\Gamma\sim 10-30$) yield energy spectra of $dN/dE\sim E^{-2}$. Highly relativistic shocks expected in Gamma Ray Bursts (GRBs, $100<\Gamma<1000$), on the other hand, have spectra as flat as $E^{-1.5}$. The model results are compared to the measured flux of Cosmic Rays at the highest energies and it is shown that, while AGN spectra are well-suited, GRB spectra are too flat to explain the observed flux. The first evidence of a correlation between the Cosmic Ray flux above $5.7\cdot 10^{10}$ GeV and the distribution of AGN by Auger are explained by the model. Neutrino production is expected in GRBs, either in mildly or highly relativistic shocks and although these sources are excluded as the principle origin of UHECRs, superluminal shocks in particular may be observable via neutrino and photon fluxes, rather than as protons.

The understanding of the particle spectra resulting from acceleration in relativistic shocks as they occur in extragalactic sources, is essential for the interpretation of the cosmic ray spectrum above the ankle ($E_p>3\cdot 10^{18}$ eV). It is believed that extragalactic sources like Active Galactic Nuclei and Gamma Ray Bursts can produce particle spectra up to $E_p\sim 10^{21}$ eV. In this contribution, subluminal shocks are investigated with respect to different shock boost factors $\Gamma$ and the inclination angle between the shock normal and the magnetic field $\psi$. A correlation between the boost factor and the spectral behavior of the emitted particles is found. The results are compared to Active Galactic Nuclei and Gamma Ray Burst diffuse cosmic ray contribution and the observed cosmic ray spectrum at the highest energies.

Results are presented from the first underground data run of ZEPLIN-II, a 31 kg two phase xenon detector developed to observe nuclear recoils from hypothetical weakly interacting massive dark matter particles. Discrimination between nuclear recoils and background electron recoils is afforded by recording both the scintillation and ionisation signals generated within the liquid xenon, with the ratio of these signals being different for the two classes of event. This ratio is calibrated for different incident species using an AmBe neutron source and Co-60 gamma-ray sources. From our first 31 live days of running ZEPLIN-II, the total exposure following the application of fiducial and stability cuts was 225 kgxdays. A background population of radon progeny events was observed in this run, arising from radon emission in the gas purification getters, due to radon daughter ion decays on the surfaces of the walls of the chamber. An acceptance window, defined by the neutron calibration data, of 50% nuclear recoil acceptance between 5 keVee and 20 keVee, had an observed count of 29 events, with a summed expectation of 28.6+/-4.3 gamma-ray and radon progeny induced background events. These figures provide a 90% c.l. upper limit to the number of nuclear recoils of 10.4 events in this acceptance window, which converts to a WIMP-nucleon spin-independent cross-section with a minimum of 6.6x10^-7 pb following the inclusion of an energy dependent, calibrated, efficiency. A second run is currently underway in which the radon progeny will be eliminated, thereby removing the background population, with a projected sensitivity of 2x10^-7 pb for similar exposures as the first run.

We present details of the technical design and manufacture of the ZEPLIN-III dark matter experiment. ZEPLIN-III is a two-phase xenon detector which measures both the scintillation light and the ionisation charge generated in the liquid by interacting particles and radiation. The instrument design is driven by both the physics requirements and by the technology requirements surrounding the use of liquid xenon. These include considerations of key performance parameters, such as the efficiency of scintillation light collection, restrictions placed on the use of materials to control the inherent radioactivity levels, attainment of high vacuum levels and chemical contamination control. The successful solution has involved a number of novel design and manufacturing features which will be of specific use to future generations of direct dark matter search experiments as they struggle with similar and progressively more demanding requirements.

We present results from a GEANT4-based Monte Carlo tool for end-to-end simulations of the ZEPLIN-III dark matter experiment. ZEPLIN-III is a two-phase detector which measures both the scintillation light and the ionisation charge generated in liquid xenon by interacting particles and radiation. The software models the instrument response to radioactive backgrounds and calibration sources, including the generation, ray-tracing and detection of the primary and secondary scintillations in liquid and gaseous xenon, and subsequent processing by data acquisition electronics. A flexible user interface allows easy modification of detector parameters at run time. Realistic datasets can be produced to help with data analysis, an example of which is the position reconstruction algorithm developed from simulated data. We present a range of simulation results confirming the original design sensitivity of a few times $10^{-8}$ pb to the WIMP-nucleon cross-section.

We perform Monte Carlo simulations of diffusive shock acceleration at highly relativistic oblique shock waves. High upstream flow Lorentz gamma factors are used, which are relevant to models of ultra relativistic particle shock acceleration in Active Galactic Nuclei (AGN) central engines and relativistic jets and Gamma Ray Burst (GRB) fireballs. We investigate numerically the acceleration properties -in the ultra relativistic flow regime of $\Gamma \sim 10-10^{3}$- such as angular distribution, acceleration time constant, particle energy gain versus number of crossings and spectral shapes. We perform calculations for sub-luminal and super-luminal shocks, using two different approaches respectively. The $\Gamma^{2}$ energization for the first crossing cycle and the significantly large energy gain for subsequent crossings as well as the high 'speed up' factors found, are important in supporting the Vietri and Waxman models on GRB ultra-high energy cosmic ray, neutrino, and gamma-ray output.

Monte-Carlo computations for highly relativistic parallel shock particle acceleration are presented for upstream flow gamma factors, $\Gamma=(1-V_{1}^{2}/c^{2})^{-0.5}$ with values between 5 and $10^{3}$. The results show that the spectral shape at the shock depends on whether or not the particle scattering is small angle with $\delta \theta < \Gamma^{-1}$ or large angle, which is possible if $\lambda > 2r_{g} \Gamma^{2}$ where $\lambda$ is the scattering mean free path along the field line and $r_{g}$ the gyroradius, these quantities being measured in the plasma flow frame. The large angle scattering case exhibits distinctive structure superimposed on the basic power-law spectrum, largely absent in the pitch angle case. Also, both cases yield an acceleration rate faster than estimated by the conventional, non-relativistic expression, $t_{acc}=[c/(V_{1}-V_{2})] [\lambda_{1}/V_{1}+\lambda_{2}/V_{2}]$ where '1' and '2' refer to upstream and downstream and $\lambda$ is the mean free path. A $\Gamma^{2}$ energy enhancement factor in the first shock crossing cycle and a significant energy multiplication in the subsequent shock cycles are also observed. The results may be important for our understanding of the production of very high energy cosmic rays and the high energy neutrino and gamma-ray output from Gamma Ray Bursts (GRB) and Active Galactic Nuclei (AGN).

The status of dark matter searches with inorganic scintillator detectors at Boulby mine is reviewed and the results of tests with a CsI(Tl) crystal are presented. The objectives of the latter experiment were to study anomalous fast events previously observed and to identify ways to remove this background. Clear indications were found that these events were due to surface contamination of crystals by alphas, probably from radon decay. A new array of unencapsulated NaI(Tl) crystals immersed either in liquid paraffin or pure nitrogen atmosphere is under construction at Boulby. Such an approach allows complete control of the surface of the crystals and the ability to remove any surface contamination. First data from the unencapsulated NaI(Tl) do not show the presence of anomalous fast events.

In this paper we present an investigation of numerical Monte Carlo simulations of the diffusive shock acceleration in the test particle limit. Very high gamma flow astrophysical plasmas, have been used, from $\gamma_{up}$ $\sim50$ up to $\gamma_{up}$ $\sim1000$, which could be relevant to the suggested models of AGNs Jets and their Central Engines as well as the ultra-relativistic shock particle acceleration in Gamma Ray Burst (GRB) fireballs. Particularly the energy gain per shock crossing and the time constant for the above high relativistic processes is numerically calculated. We explicitly find a considerable $\gamma^{2}$ energy boosting in the first shock cycle, and in all subsequent shock cycles the particle energy is multiplied by a large factor. Also a noted acceleration speed-up for the same acceleration process has been observed. Both of those results are connected with theoretical suggestions that a significant enhancement of the acceleration process is possible due to relativistic effects not present at lower plasma flow speeds. The acceleration speed-up and the mean energy-gain per shock cycle found could efficiently justify the theory of the origin of Ultra High Energy Cosmic Rays (UHECR) from the sites of GRBs.

Numerical modelling is performed for extreme relativistic parallel shocks with upstream Lorentz factor $\Gamma=50$. Assuming the scattering is either large angle or over pitch angles $\sim \Gamma^{-1}$, spectral flattening and shock aaceleration speed-up is found. The energy gain for the first shock cycle is $\sim \Gamma^{2}$. The likely output from relativistic shocks due to the infall from the accretion disk to the AGN black hole is computed. Neutrinos from proton-gamma interactions may be detectable with planned neutrino telescopes.

Numerical Monte Carlo simulations of the diffusive shock acceleration in the test particle limit are investigated. We simulate high relativistic flow astrophysical plasmas for upstream $\gamma$ $\sim5$ and up to $\gamma$ $\sim1000$. These gamma values are relevant to the diffusive shock acceleration models of Active Galactic Nuclei (AGN) and Gamma Ray Bursts (GRBs). The spectral shape for the above high relativistic processes will be calculated. This work presents distinctive spectral shapes with plateau-like characteristics due to the relativistic effects and the spectral indexes are found to be as flat as the flow becomes relativistic; the latter is in a very good agreement with previous theoretical and numerical works concerning relativistic plasma flows.

We give an overview of the current status and scientific goals of VERITAS, a proposed hexagonal array of seven 10 m aperture imaging Cherenkov telescopes. The selected site is Montosa Canyon (1390 m a.s.l.) at the Whipple Observatory, Arizona. Each telescope, of 12 m focal length, will initially be equipped with a 499 element photomultiplier camera covering a 3.5 degree field of view. A central station will initiate the readout of 500 MHz FADCs upon receipt of multiple telescope triggers. The minimum detectable flux sensitivity will be 0.5% of the Crab Nebula flux at 200 GeV. Detailed simulations of the array's performance are presented elsewhere at this meeting. VERITAS will operate primarily as a gamma-ray observatory in the 50 GeV to 50 TeV range for the study of active galaxies, supernova remnants, pulsars and gamma-ray bursts.

Clustering properties of two samples of X-ray selected AGNs are compared with IR/Optical galactic samples. These are the Einstein Observatory Extended Medium Sensitivity Survey: EMSS; and a sample of 654 ROSAT detected radio quasars. The cosmological dipole moment is calculated for these samples, yielding the values of $b \Omega_o^{-0.6}$. We find this factor to be $0.68 \pm 0.22$ for the EMSS sample, and $1.02 \pm 0.014$ for the quasars sample. Comparing these values with that of the IRAS 1.2 Jy sample $(1.95 \pm 0.096)$, suggests that these AGNs are less clustered than the IRAS galaxies. For each sample there is a depth at which the amplitude of the dipole growth curve saturates. This convergence radius $(R_{conv})$ is found to be $\sim 300-400$ h$^{-1}$ Mpc for the EMSS sample and $\sim 450-500$ h$^{-1}$ Mpc for the ROSAT detected quasars sample, which suggests that X-ray AGNs have a much deeper contribution to the Local Group (LG) motion compared to the galaxy and cluster samples. Using the properties of the angular and spatial correlation functions for these samples, we find: $b_{IRAS}/b_{quasars} = 3.61 \pm 2.7$; $b_{IRAS}/b_{EMSS}= 8.46 \pm 6.4$ and $b_{quasars} / b_{EMSS} = 2.35 \pm 0.2$. Using Chi-squared minimisation, we have fitted the angular correlation functions using the standard power law model of the form: $W(\theta)= A \theta^{-\delta}$, and the spatial correlation functions using the power law model of the form: $\xi(r) =(r/r_{0})^{-\gamma}$. For the angular correlation functions, we find $A \sim0.6$ for the AGN auto-correlations and $A \sim0.2-0.3$ for the AGN-IRAS galaxy cross-correlations. Also $\delta$ is found to be $~ 0.8-1.1$. In the spatial correlation analysis, we find $r_{o}\sim3.5-5.1$ $h^{-1}Mpc$ and $\gamma\sim1.4-1.7$ .

We present the results of the deepest optically identified X-ray survey yet made. The X-ray survey was made with the ROSAT PSPC and reaches a flux limit of 1.6x10^-15 erg cm^-2 s^-1 (0.5--2.0 keV). Above a flux limit of 2x10^-15 erg cm^-2 s^-1 we define a complete sample of 70 sources of which 59 are identified. Some (5) other sources have tentative identifications and in a further 4 the X-ray error-boxes are blank to R=23 mag. At the brighter flux levels (>= 10^-14 erg cm^-2 s^-1) we confirm the results of previous less deep X-ray surveys with 84% of the sources begin QSOs. At fainter fluxes, however, the survey is dominated by a population of galaxies with narrow optical emission lines (NELGs). In addition, a number of groups and clusters of galaxies are found at intermediate fluxes. Most of these are poor systems of low X-ray luminosity and are generally found at redshifts of > 0.3. Their numbers are consistent with a zero evolutionary scenario, in contrast to the situation for high luminosity clusters at the same redshift. We discuss the significance of these results to the determination of the cosmic soft X-ray background (XRB) and show that at 2x10^-15 erg cm^-2 s^-1, we have resolved more than 50% of the background. We also briefly consider the probable importance of NELG objects to the residual background and look at some of the properties of these unusual objects.