A multi-target detection system XAX, comprising concentric 10 ton targets of 136Xe and 129/131Xe, together with a geometrically similar or larger target of liquid Ar, is described. Each is configured as a two-phase scintillation/ionization TPC detector, enhanced by a full 4pi array of ultra-low radioactivity Quartz Photon Intensifying Detectors (QUPIDs) replacing the conventional photomultipliers for detection of scintillation light. It is shown that background levels in XAX can be reduced to the level required for dark matter particle (WIMP) mass measurement at a 10^-10 pb WIMP-nucleon cross section, with single-event sensitivity below 10^-11 pb. The use of multiple target elements allows for confirmation of the A^2 dependence of a coherent cross section, and the different Xe isotopes provide information on the spin-dependence of the dark matter interaction. The event rates observed by Xe and Ar would modulate annually with opposite phases from each other for WIMP mass >~100 GeV/c^2. The large target mass of 136Xe and high degree of background reduction allow neutrinoless double beta decay to be observed with lifetimes of 10^27-10^28 years, corresponding to the Majorana neutrino mass range 0.01-0.1 eV, the most likely range from observed neutrino mass differences. The use of a 136Xe-depleted 129/131Xe target will also allow measurement of the pp solar neutrino spectrum to a precision of 1-2%.

I present here a review of past and present multi-disciplinary research of the Pittsburgh Computational AstroStatistics (PiCA) group. This group is dedicated to developing fast and efficient statistical algorithms for analysing huge astronomical data sources. I begin with a short review of multi-resolutional kd-trees which are the building blocks for many of our algorithms. For example, quick range queries and fast n-point correlation functions. I will present new results from the use of Mixture Models (Connolly et al. 2000) in density estimation of multi-color data from the Sloan Digital Sky Survey (SDSS). Specifically, the selection of quasars and the automated identification of X-ray sources. I will also present a brief overview of the False Discovery Rate (FDR) procedure (Miller et al. 2001a) and show how it has been used in the detection of ``Baryon Wiggles'' in the local galaxy power spectrum and source identification in radio data. Finally, I will look forward to new research on an automated Bayes Network anomaly detector and the possible use of the Locally Linear Embedding algorithm (LLE; Roweis & Saul 2000) for spectral classification of SDSS spectra.