Fast Calculation of the Lomb-Scargle Periodogram Using Graphics Processing Units

I introduce a new code for fast calculation of the Lomb-Scargle periodogram, that leverages the computing power of graphics processing units (GPUs). After establishing a background to the newly emergent field of GPU computing, I discuss the code design and narrate key parts of its source. Benchmarking calculations indicate no significant differences in accuracy compared to an equivalent CPU-based code. However, the differences in performance are pronounced; running on a low-end GPU, the code can match 8 CPU cores, and on a high-end GPU it is faster by a factor approaching thirty. Applications of the code include analysis of long photometric time series obtained by ongoing satellite missions and upcoming ground-based monitoring facilities; and Monte-Carlo simulation of periodogram statistical properties.Comment: Accepted by ApJ. Accompanying program source (updated since acceptance) can be downloaded from http://www.astro.wisc.edu/~townsend/resource/download/code/culsp.tar.g

3D Reconstruction of the Density Field: An SVD Approach to Weak Lensing Tomography

We present a new method for constructing three-dimensional mass maps from gravitational lensing shear data. We solve the lensing inversion problem using truncation of singular values (within the context of generalized least squares estimation) without a priori assumptions about the statistical nature of the signal. This singular value framework allows a quantitative comparison between different filtering methods: we evaluate our method beside the previously explored Wiener filter approaches. Our method yields near-optimal angular resolution of the lensing reconstruction and allows cluster sized halos to be de-blended robustly. It allows for mass reconstructions which are 2-3 orders-of-magnitude faster than the Wiener filter approach; in particular, we estimate that an all-sky reconstruction with arcminute resolution could be performed on a time-scale of hours. We find however that linear, non-parametric reconstructions have a fundamental limitation in the resolution achieved in the redshift direction.Comment: 11 pages, 6 figures. Accepted for publication in Ap

Transit surveys for Earths in the habitable zones of white dwarfs

To date the search for habitable Earth-like planets has primarily focused on nuclear burning stars. I propose that this search should be expanded to cool white dwarf stars that have expended their nuclear fuel. I define the continuously habitable zone of white dwarfs, and show that it extends from ~0.005 to 0.02 AU for white dwarfs with masses from 0.4 to 0.9 solar masses, temperatures less than 10,000 K, and habitable durations of at least 3 Gyr. As they are similar in size to Earth, white dwarfs may be deeply eclipsed by terrestrial planets that orbit edge-on, which can easily be detected with ground-based telescopes. If planets can migrate inward or reform near white dwarfs, I show that a global robotic telescope network could carry out a transit survey of nearby white dwarfs placing interesting constraints on the presence of habitable Earths. If planets were detected, I show that the survey would favor detection of planets similar to Earth: similar size, temperature, rotation period, and host star temperatures similar to the Sun. The Large Synoptic Survey Telescope could place even tighter constraints on the frequency of habitable Earths around white dwarfs. The confirmation and characterization of these planets might be carried out with large ground and space telescopes.Comment: 5 pages, 4 figures, published in ApJL; Figure 4 corrected, other minor change

Maximizing survey volume for large-area multi-epoch surveys with Voronoi tessellation

The survey volume of a proper motion-limited sample is typically much smaller than a magnitude-limited sample. This is because of the noisy astrometric measurements from detectors that are not dedicated for astrometric missions. In order to apply an empirical completeness correction, existing works limit the survey depth to the shallower parts of the sky that hamper the maximum potential of a survey. The number of epoch of measurement is a discrete quantity that cannot be interpolated across the projected plane of observation, so that the survey properties change in discrete steps across the sky. This work proposes a method to dissect the survey into small parts with Voronoi tessellation using candidate objects as generating points such that each part defines a ‘mini-survey’ that has its own properties. Coupling with a maximum volume density estimator, the new method is demonstrated to be unbiased and recovered ∼20 per cent more objects than the existing method in a mock catalogue of a white dwarf-only solar neighbourhood with Pan–STARRS 1-like characteristics. Towards the end of this work, we demonstrate one way to increase the tessellation resolution with artificial generating points, which would be useful for analysis of rare objects with small number counts

Large-k Limit of Multi-Point Propagators in the RG Formalism

Renormalized versions of cosmological perturbation theory have been very successful in recent years in describing the evolution of structure formation in the weakly non-linear regime. The concept of multi-point propagators has been introduced as a tool to quantify the relation between the initial matter distribution and the final one and to push the validity of the approaches to smaller scales. We generalize the n-point propagators that have been considered until now to include a new class of multi-point propagators that are relevant in the framework of the renormalization group formalism. The large-k results obtained for this general class of multi-point propagators match the results obtained earlier both in the case of Gaussian and non-Gaussian initial conditions. We discuss how the large-k results can be used to improve on the accuracy of the calculations of the power spectrum and bispectrum in the presence of initial non-Gaussianities.Comment: 30 page

Dynamical Capture Binary Neutron Star Mergers

We study dynamical capture binary neutron star mergers as may arise in dense stellar regions such as globular clusters. Using general-relativistic hydrodynamics, we find that these mergers can result in the prompt collapse to a black hole or in the formation of a hypermassive neutron star, depending not only on the neutron star equation of state but also on impact parameter. We also find that these mergers can produce accretion disks of up to a tenth of a solar mass and unbound ejected material of up to a few percent of a solar mass. We comment on the gravitational radiation and electromagnetic transients that these sources may produce.Comment: 6 pages, 4 figures; revised to match published versio

Expected Large Synoptic Survey Telescope (LSST) Yield of Eclipsing Binary Stars

In this paper we estimate the Large Synoptic Survey Telescope (LSST) yield of eclipsing binary stars, which will survey ~20,000 square degrees of the southern sky during the period of 10 years in 6 photometric passbands to r ~ 24.5. We generate a set of 10,000 eclipsing binary light curves sampled to the LSST time cadence across the whole sky, with added noise as a function of apparent magnitude. This set is passed to the Analysis of Variance (AoV) period finder to assess the recoverability rate for the periods, and the successfully phased light curves are passed to the artificial intelligence-based pipeline EBAI to assess the recoverability rate in terms of the eclipsing binaries' physical and geometric parameters. We find that, out of ~24 million eclipsing binaries observed by LSST with S/N>10 in mission life-time, ~28% or 6.7 million can be fully characterized by the pipeline. Of those, ~25% or 1.7 million will be double-lined binaries, a true treasure trove for stellar astrophysics.Comment: 19 pages, 7 figures. Accepted to AJ, to appear in issue 142:2 (Aug 2011

Strongly lensed SNe Ia in the era of LSST: observing cadence for lens discoveries and time-delay measurements

The upcoming Large Synoptic Survey Telescope (LSST) will detect many strongly lensed Type Ia supernovae (LSNe Ia) for time-delay cosmography. This will provide an independent and direct way for measuring the Hubble constant $H_0$, which is necessary to address the current $4.4 \sigma$ tension in $H_0$ between the local distance ladder and the early Universe measurements. We present a detailed analysis of different observing strategies for the LSST, and quantify their impact on time-delay measurement between multiple images of LSNe Ia. For this, we produced microlensed mock-LSST light curves for which we estimated the time delay between different images. We find that using only LSST data for time-delay cosmography is not ideal. Instead, we advocate using LSST as a discovery machine for LSNe Ia, enabling time delay measurements from follow-up observations from other instruments in order to increase the number of systems by a factor of 2 to 16 depending on the observing strategy. Furthermore, we find that LSST observing strategies, which provide a good sampling frequency (the mean inter-night gap is around two days) and high cumulative season length (ten seasons with a season length of around 170 days per season), are favored. Rolling cadences subdivide the survey and focus on different parts in different years; these observing strategies trade the number of seasons for better sampling frequency. In our investigation, this leads to half the number of systems in comparison to the best observing strategy. Therefore rolling cadences are disfavored because the gain from the increased sampling frequency cannot compensate for the shortened cumulative season length. We anticipate that the sample of lensed SNe Ia from our preferred LSST cadence strategies with rapid follow-up observations would yield an independent percent-level constraint on $H_0$.Comment: 25 pages, 22 figures; accepted for publication in A&

Cosmological parameter constraints from galaxy–galaxy lensing and galaxy clustering with the SDSS DR7

Recent studies have shown that the cross-correlation coefficient between galaxies and dark matter is very close to unity on scales outside a few virial radii of galaxy haloes, independent of the details of how galaxies populate dark matter haloes. This finding makes it possible to determine the dark matter clustering from measurements of galaxy–galaxy weak lensing and galaxy clustering. We present new cosmological parameter constraints based on large-scale measurements of spectroscopic galaxy samples from the Sloan Digital Sky Survey (SDSS) data release 7. We generalize the approach of Baldauf et al. to remove small-scale information (below 2 and 4 h^(−1) Mpc for lensing and clustering measurements, respectively), where the cross-correlation coefficient differs from unity. We derive constraints for three galaxy samples covering 7131 deg^2, containing 69 150, 62 150 and 35 088 galaxies with mean redshifts of 0.11, 0.28 and 0.40. We clearly detect scale-dependent galaxy bias for the more luminous galaxy samples, at a level consistent with theoretical expectations. When we vary both σ_8 and Ω_m (and marginalize over non-linear galaxy bias) in a flat Λ cold dark matter model, the best-constrained quantity is σ_8(Ω_m/0.25)^(0.57) = 0.80 ± 0.05 (1σ, stat. + sys.), where statistical and systematic errors (photometric redshift and shear calibration) have comparable contributions, and we have fixed n_s = 0.96 and h = 0.7. These strong constraints on the matter clustering suggest that this method is competitive with cosmic shear in current data, while having very complementary and in some ways less serious systematics. We therefore expect that this method will play a prominent role in future weak lensing surveys. When we combine these data with Wilkinson Microwave Anisotropy Probe 7-year (WMAP7) cosmic microwave background (CMB) data, constraints on σ_8, Ω_m, H_0, w_(de) and ∑m_ν become 30–80 per cent tighter than with CMB data alone, since our data break several parameter degeneracies