118430259
Fast Algorithms and Efficient Statistics: Density Estimation in Large Astronomical Datasets
In this paper, we outline the use of Mixture Models in density estimation of large astronomical databases. This method of density estimation has been known in Statistics for some time but has not been implemented because of the large computational cost. Herein, we detail an implementation of the Mixture Model density estimation based on multi-resolutional KD-trees which makes this statistical technique into a computationally tractable problem. We provide the theoretical and experimental background for using a mixture model of Gaussians based on the Expectation Maximization (EM) Algorithm. Applying these analyses to simulated data sets we show that the EM algorithm - using the AIC penalized likelihood to score the fit - out-performs the best kernel density estimate of the distribution while requiring no ``fine--tuning'' of the input algorithm parameters. We find that EM can accurately recover the underlying density distribution from point processes thus providing an efficient adaptive smoothing method for astronomical source catalogs. To demonstrate the general application of this statistic to astrophysical problems we consider two cases of density estimation: the clustering of galaxies in redshift space and the clustering of stars in color space. From these data we show that EM provides an adaptive smoothing of the distribution of galaxies in redshift space (describing accurately both the small and large-scale features within the data) and a means of identifying outliers in multi-dimensional color-color space (e.g. for the identification of high redshift QSOs). Automated tools such as those based on the EM algorithm will be needed in the analysis of the next generation of astronomical catalogs (2MASS, FIRST, PLANCK, SDSS) and ultimately in in the development of the National Virtual Observatory.
118985428
Three-Point Correlation Functions of SDSS Galaxies in Redshift Space: Morphology, Color, and Luminosity Dependence
We present measurements of the redshift-space three-point correlation function of galaxies in the Sloan Digital Sky Survey (SDSS). For the first time, we analyze the dependence of this statistic on galaxy morphology, color, and luminosity. In order to control systematics due to selection effects, we used r-band, volume-limited samples of galaxies, constructed from magnitude-limited SDSS data (14.5 <r< 17.5), and further divided the samples into two morphological types (early and late) or two color populations (red and blue). The three-point correlation function of SDSS galaxies follows the hierarchical relation well, and the reduced three-point amplitudes in redshift-space are almost scale-independent (Qz =0 .5 ∼ 1.0). In addition, their dependence on the morphology, color, and luminosity is not statistically significant. Given the robust morphological, color, and luminosity dependences of the two-point correlation function, this implies that galaxy biasing is complex on weakly non-linear to non-linear scales. We show that a simple deterministic linear relation with the underlying mass could not explain our measurements on these scales.
16921302
Toward a resolution of the discrepancy between different estimators of star formation rate
Different wavelength regimes and methods for estimating the space density of the star formation rate (SFR) result in discrepant values. While it is recognized that ultraviolet (UV) and Hα emission-line data must be corrected for the effects of extinction, the magnitude of the required correction is uncertain. Even when these corrections are made there remains a significant discrepancy between SFRs derived from UV and Hα measurements compared with those derived from far-infrared (FIR) and radio luminosities. Since the FIR-radio–derived SFRs are not affected by extinction and simple corrections to reconcile the UV and Hα measurement with these do not fully account for the discrepancies, a more sophisticated correction may be required. Recent results suggest that at least part of the solution may be a form of extinction that increases with increasing SFR (or luminosity, given the common assumption that SFR is proportional to luminosity). We present an analysis of the effects of a dust reddening dependent on star formation rate applied to estimators of SFR. We show (1) that the discrepancies between Hα and FIR-radio SFR estimates may be explained by such an effect and we present an iterative method for applying the correction and (2) that UV-based estimates of SFR are harder to reconcile with FIR-radio estimates using this method, although the extent of the remaining discrepancy is less than for a non-SFR–dependent correction. Particularly at high redshift, our understanding of extinction at UV wavelengths may require a still more complex explanation.
119214510
INTERPOLATING MASKED WEAK-LENSING SIGNAL WITH KARHUNEN-LOÈVE ANALYSIS
We explore the utility of Karhunen-Lo?ve (KL) analysis in solving practical problems in the analysis of gravitational shear surveys. Shear catalogs from large-field weak-lensing surveys will be subject to many systematic limitations, notably incomplete coverage and pixel-level masking due to foreground sources. We develop a method to use two-dimensional KL eigenmodes of shear to interpolate noisy shear measurements across masked regions. We explore the results of this method with simulated shear catalogs, using statistics of high-convergence regions in the resulting map. We find that the KL procedure not only minimizes the bias due to masked regions in the field, it also reduces spurious peak counts from shape noise by a factor of ~3 in the cosmologically sensitive regime. This indicates that KL reconstructions of masked shear are not only useful for creating robust convergence maps from masked shear catalogs, but also offer promise of improved parameter constraints within studies of shear peak statistics.
8278022
First Results from the SPICES Survey
We present first results from SPICES, the Spectroscopic, Photometric, Infrared-Chosen Extragalactic Survey. SPICES is comprised of four ~30 square arcminute high Galactic latitude fields with deep BRIzJK imaging reaching depths of ~25th magnitude (AB) in the optical and ~23rd magnitude (AB) in the near-infrared. To date we have 626 spectroscopic redshifts for infrared-selected SPICES sources with K 1. One of the SPICES fields has been the target of a deep 190 ksec Chandra exposure; we discuss initial results from analysis of that data set. Finally, we briefly discuss a successful campaign to identify high-redshift sources in the SPICES fields.
119191264
Classification of Stellar Spectra with LLE
We investigate the use of dimensionality reduction techniques for the classification of stellar spectra selected from the SDSS. Using local linear embedding (LLE), a technique that preserves the local (and possibly non-linear) structure within high dimensional data sets, we show that the majority of stellar spectra can be represented as a one dimensional sequence within a three dimensional space. The position along this sequence is highly correlated with spectral temperature. Deviations from this “stellar locus” are indicative of spectra with strong emission lines (including misclassified galaxies) or broad absorption lines (e.g. Carbon stars). Based on this analysis, we propose a hierarchical classification scheme using LLE that progressively identifies and classifies stellar spectra in a manner that requires no feature extraction and that can reproduce the classic MK classifications to an accuracy of one type.
34851507
Constraints on Cosmic Neutrino Fluxes from the ANITA Experiment
P. F. Dowkontt, 4 M. A. DuVernois,5 P. A. Evenson, 6 D. Goldstein, 1 P. W. Gorham, 9 C. L. Hebert, 9 M. H. Israel,4 J. G. Learned, 9 K. M. Liewer,10 J. T. Link,9 S. Matsuno, 9 P. Miočinović,9 J. Nam, 1 C. J. Naudet, 10 R. Nichol,2 K. Palladino, 2 M. Rosen, 9 D. Saltzberg, 7 D. Seckel, 6 A. Silvestri,1 B. T. Stokes, 9 G. S. Varner, 9 and F. Wu1 1Department of Physics and Astronomy, University of California at Irvine, Irvine, California 2Department of Physics, Ohio State University, Columbus, Ohio 3Department of Physics and Astronomy, University of Kansas, Lawrence, Kansas 4Department of Physics, Washington University in St. Louis, St. Louis, Missouri 5School of Physics and Astronomy, University of Minnesota, Minneapolis, Minnesota 6Bartol Research Institute, University of Delaware, Newark, Delaware 7Department of Physics and Astronomy, University of California at Los Angeles, Los Angeles, California 8Department of Astronomy and Astrophysics, Pennsylvania State University, University Park, Pennsylvania 9Department of Physics and Astronomy, University of Hawaii at Manoa, Honolulu, Hawaii 10Jet Propulsion Laboratory, Pasadena, California
45943105
Tuning into UHE Neutrinos in Antarctica - The ANITA Experiment
ice will produce a radio pulse by the means of the Askaryan effect. The large radio transparency of ice allows for such a pulse to be recorded by a cluster of balloon-borne antennas. The details of the ANITA instrument, now in a construction phase, and the science we hope to achieve is discussed. In order to prepare for the main mission, we have flown ANITA-lite during the 2003/04 austral season. ANITA-lite consisted of two quad-ridge horn antennas and a prototype RF (radio frequency) triggering and recording system. Here we present the results of an impulsive RF background survey of Antarctica, as well as proof-of-principle gain, tracking, and timing calibrations conducted by observing solar radio emissions and calibration radio-pulses. A preliminary UHE neutrino flux limit based on ANITA-lite data is also presented.
120506645
Results from the ANITA experiment (Analysis A)
The ANITA (ANtarctic Impulsive Transient Antenna) is a balloon-borne neutrino telescope which consists of an array of 32 broadband horn antennas. It successfully completed a 35 day flight over Antarctica during the 2006-2007 austral summer. The primary goal of ANITA is to search for astrophysical neutrinos with energies E > 1019 eV by detecting radio Cherenkov signals from neutrino-induced showers in the Antarctic ice. There are two independent data analysis groups in the ANITA experiment based on data blinding and analysis methods. We present results from one of the analysis groups with an event blinding approach.
123575242
Introduction to the SalSA, a saltdome shower array as a GZK neutrino observatory
The observed spectrum of ultra-high energy cosmic rays virtually guarantees the presence of ultra-high energy neutrinos due to their interaction with the cosmic microwave background. Every one of these neutrinos will point back to its source and, unlike cosmic rays, will arrive at the Earth unattenuated, from sources perhaps as distant as z=20. The neutrino telescopes currently under construction, should discover a handful of these events, probably too few for detailed study. In this talk I will describe how an array of VHF and UHF antennas embedded in a large salt dome, SalSA (Saltdome Shower Array) promises to yield a teraton detector (> 500 km3-sr) for contained neutrino events with energies above 1017 eV. Our simulations show that such a detector may observe several hundreds of these neutrinos over its lifetime. Our simulations also show how such interactions will provide high energy physicists with an energy frontier for weak interactions an order-of-magnitude larger than that of the LHC. The flavor ID capalities of SALSA, combined with the extreme L/E of these neutrinos, will provide a window on neutrino oscillations and decay times eight orders of magnitude higher than laboratory experiments. In addition to the latest simulation results, we describe progress on detectors and site selection.
119288671
The Effective Number Density of Galaxies for Weak Lensing Measurements in the LSST Project
Future weak lensing surveys potentially hold the highest statistical power for constraining cosmological parameters compared to other cosmological probes. The statistical power of a weak lensing survey is determined by the sky coverage, the inverse of the noise in shear measurements, and the galaxy number density. The combination of the latter two factors is often expressed in terms of $n_{\rm eff}$ -- the "effective number density of galaxies used for weak lensing measurements". In this work, we estimate $n_{\rm eff}$ for the Large Synoptic Survey Telescope (LSST) project, the most powerful ground-based lensing survey planned for the next two decades. We investigate how the following factors affect the resulting $n_{\rm eff}$ of the survey with detailed simulations: (1) survey time, (2) shear measurement algorithm, (3) algorithm for combining multiple exposures, (4) inclusion of data from multiple filter bands, (5) redshift distribution of the galaxies, and (6) masking and blending. For the first time, we quantify in a general weak lensing analysis pipeline the sensitivity of $n_{\rm eff}$ to the above factors. 
We find that with current weak lensing algorithms, expected distributions of observing parameters, and all lensing data ($r$- and $i$-band, covering 18,000 degree$^{2}$ of sky) for LSST, $n_{\rm eff} \approx37$ arcmin$^{-2}$ before considering blending and masking, $n_{\rm eff} \approx31$ arcmin$^{-2}$ when rejecting seriously blended galaxies and $n_{\rm eff} \approx26$ arcmin$^{-2}$ when considering an additional 15% loss of galaxies due to masking. With future improvements in weak lensing algorithms, these values could be expected to increase by up to 20%. Throughout the paper, we also stress the ways in which $n_{\rm eff}$ depends on our ability to understand and control systematic effects in the measurements.