Source author record

Daniel E. Holz

Daniel E. Holz appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

gr-qc astro-ph.CO astro-ph.HE astro-ph astro-ph.IM hep-ph astro-ph.GA astro-ph.SR hep-th math.ST Statistics Theory

Catalog footprint

What is connected

57works

11topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2026arXiv

False positive control in time series coincidence detection

We study the problem of coincidence detection in time series data, where we aim to determine whether the appearance of simultaneous or near-simultaneous events in two time series is indicative of some shared underlying signal or synchronicity, or might simply be due to random chance. This problem arises across many applications, such as astrophysics (e.g., detecting astrophysical events such as gravitational waves, with two or more detectors) and neuroscience (e.g., detecting synchronous firing patterns between two or more neurons). In this work, we consider methods based on time-shifting, where the timeline of one data stream is randomly shifted relative to another, to mimic the types of coincidences that could occur by random chance. Our theoretical results establish rigorous finite-sample guarantees controlling the probability of false positives, under weak assumptions that allow for dependence within the time series data, providing reassurance that time-shifting methods are a reliable tool for inference in this setting. Empirical results with simulated and real data validate the strong performance of time-shifting methods in dependent-data settings.

preprint2022arXiv

Avoiding a Cluster Catastrophe: Retention Efficiency and the Binary Black Hole Mass Spectrum

The population of binary black hole mergers identified through gravitational waves has uncovered unexpected features in the intrinsic properties of black holes in the universe. One particularly surprising and exciting result is the possible existence of black holes in the pair-instability mass gap, $\sim50-120~M_\odot$. Dense stellar environments can populate this region of mass space through hierarchical mergers, with the retention efficiency of black hole merger products strongly dependent on the escape velocity of the host environment. We use simple toy models to represent hierarchical merger scenarios in various dynamical environments. We find that hierarchical mergers in environments with high escape velocities $(\gtrsim300~\mathrm{km\,s}^{-1})$ are efficiently retained. If such environments dominate the binary black hole merger rate, this would lead to an abundance of high-mass mergers that is potentially incompatible with the empirical mass spectrum from the current catalog of binary black hole mergers. Models that efficiently generate hierarchical mergers, and contribute significantly to the observed population, must therefore be tuned to avoid a "cluster catastrophe" of overproducing binary black hole mergers within and above the pair-instability mass gap.

preprint2022arXiv

Bridging the Gap: Categorizing Gravitational-Wave Events at the Transition Between Neutron Stars and Black Holes

We search for features in the mass distribution of detected compact binary coalescences which signify the transition between neutron stars and black holes. We analyze all gravitational wave detections by LIGO-Virgo-KAGRA made through the end of the first half of the third observing run, and find clear evidence for two different populations of compact objects based solely on gravitational wave data. We confidently (99.3%) find a steepening relative to a single power law describing neutron stars and low-mass black holes below $2.4^{+0.5}_{-0.5}\text{ M}_{\odot}$, which is consistent with many predictions for the maximum neutron star mass. We find suggestions of the purported lower mass gap between the most massive neutron stars and the least massive black holes, but are unable to conclusively resolve it with current data. If it exists, we find the lower mass gap's edges to lie at $2.2^{+0.7}_{-0.5}\text{ M}_{\odot}$ and $6.0^{+2.4}_{-1.4}\text{ M}_{\odot}$. We re-examine events that have been deemed "exceptional" by the LIGO-Virgo-KAGRA collaborations in the context of these features. We analyze GW190814 self-consistently in the context of the full population of compact binaries, finding support for its secondary to be either a neutron star or a lower mass gap object, consistent with previous claims. Our models are the first to accommodate this event, which is an outlier with respect to the binary black hole population. We find that GW200105 and GW200115 probe the edges of, and may have components within, the lower mass gap. As future data improve global population models, the classification of these events will also improve.

preprint2022arXiv

Cosmology Intertwined: A Review of the Particle Physics, Astrophysics, and Cosmology Associated with the Cosmological Tensions and Anomalies

In this paper we will list a few important goals that need to be addressed in the next decade, also taking into account the current discordances between the different cosmological probes, such as the disagreement in the value of the Hubble constant $H_0$, the $σ_8$--$S_8$ tension, and other less statistically significant anomalies. While these discordances can still be in part the result of systematic errors, their persistence after several years of accurate analysis strongly hints at cracks in the standard cosmological scenario and the necessity for new physics or generalisations beyond the standard model. In this paper, we focus on the $5.0\,σ$ tension between the {\it Planck} CMB estimate of the Hubble constant $H_0$ and the SH0ES collaboration measurements. After showing the $H_0$ evaluations made from different teams using different methods and geometric calibrations, we list a few interesting new physics models that could alleviate this tension and discuss how the next decade's experiments will be crucial. Moreover, we focus on the tension of the {\it Planck} CMB data with weak lensing measurements and redshift surveys, about the value of the matter energy density $Ω_m$, and the amplitude or rate of the growth of structure ($σ_8,fσ_8$). We list a few interesting models proposed for alleviating this tension, and we discuss the importance of trying to fit a full array of data with a single model and not just one parameter at a time. Additionally, we present a wide range of other less discussed anomalies at a statistical significance level lower than the $H_0$--$S_8$ tensions which may also constitute hints towards new physics, and we discuss possible generic theoretical approaches that can collectively explain the non-standard nature of these signals.[Abridged]

preprint2022arXiv

Snowmass2021 Cosmic Frontier White Paper: Fundamental Physics and Beyond the Standard Model

Gravitational wave detectors are formidable tools to explore strong-field gravity, especially black holes and neutron stars. These compact objects are extraordinarily efficient at producing electromagnetic and gravitational radiation. As such, they are ideal laboratories for fundamental physics and have an immense discovery potential. The detection of black hole binaries by third-generation Earth-based detectors, space-based detectors and pulsar timing arrays will provide exquisite tests of general relativity. Loud "golden" events and extreme mass-ratio inspirals can strengthen the observational evidence for horizons by mapping the exterior spacetime geometry, inform us on possible near-horizon modifications, and perhaps reveal a breakdown of Einstein's gravity. Measurements of the black-hole spin distribution and continuous gravitational-wave searches can turn black holes into efficient detectors of ultralight bosons across ten or more orders of magnitude in mass. A precise monitoring of the phase of inspiralling binaries can constrain the existence of additional propagating fields and characterize the environment in which the binaries live, bounding the local dark matter density and properties. Gravitational waves from compact binaries will probe general relativity and fundamental physics in previously inaccessible regimes, and allow us to address fundamental issues in our current understanding of the cosmos.

preprint2022arXiv

Spectral sirens: cosmology from the full mass distribution of compact binaries

We explore the use of the mass spectrum of neutron stars and black holes in gravitational-wave compact binary sources as a cosmological probe. These standard siren sources provide direct measurements of luminosity distance. In addition, features in the mass distribution, such as mass gaps or peaks, will redshift, and thus provide independent constraints on their redshift distribution. We argue that the entire mass spectrum should be utilized to provide cosmological constraints. For example, we find that the mass spectrum of LIGO--Virgo--KAGRA events introduces at least five independent mass ``features'': the upper and lower edges of the pair instability supernova (PISN) gap, the upper and lower edges of the neutron star--black hole gap, and the minimum neutron star mass. We find that although the PISN gap dominates the cosmological inference with current detectors (2G), as shown in previous work, it is the lower mass gap that will provide the most powerful constraints in the era of Cosmic Explorer and Einstein Telescope (3G). By using the full mass distribution, we demonstrate that degeneracies between mass evolution and cosmological evolution can be broken, unless an astrophysical conspiracy shifts all features of the full mass distribution simultaneously following the (non-trivial) Hubble diagram evolution. We find that this self-calibrating ``spectral siren'' method has the potential to provide precision constraints of both cosmology and the evolution of the mass distribution, with 2G achieving better than $10\%$ precision on $H(z)$ at $z\lesssim1$ within a year, and 3G reaching $\lesssim1\%$ at $z\gtrsim2$ within one month.

preprint2022arXiv

Target of Opportunity Observations of Gravitational Wave Events with Vera C. Rubin Observatory

The discovery of the electromagnetic counterpart to the binary neutron star merger GW170817 has opened the era of gravitational-wave multi-messenger astronomy. Rapid identification of the optical/infrared kilonova enabled a precise localization of the source, which paved the way to deep multi-wavelength follow-up and its myriad of related science results. Fully exploiting this new territory of exploration requires the acquisition of electromagnetic data from samples of neutron star mergers and other gravitational wave sources. After GW170817, the frontier is now to map the diversity of kilonova properties and provide more stringent constraints on the Hubble constant, and enable new tests of fundamental physics. The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) can play a key role in this field in the 2020s, when an improved network of gravitational-wave detectors is expected to reach a sensitivity that will enable the discovery of a high rate of merger events involving neutron stars (about tens per year) out to distances of several hundred Mpc. We design comprehensive target-of-opportunity observing strategies for follow-up of gravitational-wave triggers that will make the Rubin Observatory the premier instrument for discovery and early characterization of neutron star and other compact object mergers, and yet unknown classes of gravitational wave events.

preprint2021arXiv

Distance measures in gravitational-wave astrophysics and cosmology

We present quantities which characterize the sensitivity of gravitational-wave observatories to sources at cosmological distances. In particular, we introduce and generalize the horizon, range, response, and reach distances. These quantities incorporate a number of important effects, including cosmologically well-defined distances and volumes, cosmological redshift, cosmological time dilation, and rate density evolution. In addition, these quantities incorporate unique aspects of gravitational wave detectors, such as the variable sky sensitivity of the detectors and the scaling of the sensitivity with inverse distance. An online calculator (https://users.rcc.uchicago.edu/~dholz/gwc/) and python notebook (https://github.com/hsinyuc/distancetool) to determine GW distances are available. We provide answers to the question: "How far can gravitational-wave detectors hear?"

preprint2021arXiv

Please repeat: Strong lensing of gravitational waves as a probe of compact binary and galaxy populations

Strong gravitational lensing of gravitational wave sources offers a novel probe of both the lens galaxy and the binary source population. In particular, the strong lensing event rate and the time delay distribution of multiply-imaged gravitational-wave binary coalescence events can be used to constrain the mass distribution of the lenses as well as the intrinsic properties of the source population. We calculate the strong lensing event rate for a range of second (2G) and third generation (3G) detectors, including Advanced LIGO/Virgo, A+, Einstein Telescope (ET), and Cosmic Explorer (CE). For 3G detectors, we find that {$\sim0.1\%$} of observed events are expected to be strongly lensed. We predict detections of {$\sim 1$} lensing pair per year with A+, and {$\sim 50$} pairs {per year} with ET/CE. These rates are highly sensitive to the characteristic galaxy velocity dispersion, $σ_*$, implying that observations of the rates will be a sensitive probe of lens properties. We explore using the time delay distribution between multiply-imaged gravitational-wave sources to constrain properties of the lenses. We find that 3G detectors would constrain $σ_*$ to {$\sim21\%$ after 5 years}. Finally, we show that the presence or absence of strong lensing {within the detected population} provides useful insights into the source redshift and mass distribution out to redshifts beyond the peak of the star formation rate, which can be used to constrain formation channels and their relation to the star formation rate and delay time distributions for these systems.

preprint2021arXiv

Probing Extremal Gravitational-Wave Events with Coarse-Grained Likelihoods

As catalogs of gravitational-wave transients grow, new records are set for the most extreme systems observed to date. The most massive observed black holes probe the physics of pair instability supernovae while providing clues about the environments in which binary black hole systems are assembled. The least massive black holes, meanwhile, allow us to investigate the purported neutron star-black hole mass gap, and binaries with unusually asymmetric mass ratios or large spins inform our understanding of binary and stellar evolution. Existing outlier tests generally implement leave-one-out analyses, but these do not account for the fact that the event being left out was by definition an extreme member of the population. This results in a bias in the evaluation of outliers. We correct for this bias by introducing a coarse-graining framework to investigate whether these extremal events are true outliers or whether they are consistent with the rest of the observed population. Our method enables us to study extremal events while testing for population model misspecification. We show that this ameliorates biases present in the leave-one-out analyses commonly used within the gravitational-wave community. Applying our method to results from the second LIGO--Virgo transient catalog, we find qualitative agreement with the conclusions of Abbott et al, ApJL 913 L7 (2021). GW190814 is an outlier because of its small secondary mass. We find that neither GW190412 nor GW190521 are outliers.

preprint2021arXiv

The Gravity Collective: A Search for the Electromagnetic Counterpart to the Neutron Star-Black Hole Merger GW190814

We present optical follow-up imaging obtained with the Katzman Automatic Imaging Telescope, Las Cumbres Observatory Global Telescope Network, Nickel Telescope, Swope Telescope, and Thacher Telescope of the LIGO/Virgo gravitational wave (GW) signal from the neutron star-black hole (NSBH) merger GW190814. We searched the GW190814 localization region (19 deg$^{2}$ for the 90th percentile best localization), covering a total of 51 deg$^{2}$ and 94.6% of the two-dimensional localization region. Analyzing the properties of 189 transients that we consider as candidate counterparts to the NSBH merger, including their localizations, discovery times from merger, optical spectra, likely host-galaxy redshifts, and photometric evolution, we conclude that none of these objects are likely to be associated with GW190814. Based on this finding, we consider the likely optical properties of an electromagnetic counterpart to GW190814, including possible kilonovae and short gamma-ray burst afterglows. Using the joint limits from our follow-up imaging, we conclude that a counterpart with an $r$-band decline rate of 0.68 mag day$^{-1}$, similar to the kilonova AT 2017gfo, could peak at an absolute magnitude of at most $-17.8$ mag (50% confidence). Our data are not constraining for ''red'' kilonovae and rule out ''blue'' kilonovae with $M>0.5 M_{\odot}$ (30% confidence). We strongly rule out all known types of short gamma-ray burst afterglows with viewing angles $<$17$^{\circ}$ assuming an initial jet opening angle of $\sim$$5.2^{\circ}$ and explosion energies and circumburst densities similar to afterglows explored in the literature. Finally, we explore the possibility that GW190814 merged in the disk of an active galactic nucleus, of which we find four in the localization region, but we do not find any candidate counterparts among these sources.

preprint2020arXiv

Black Hole Coagulation: Modeling Hierarchical Mergers in Black Hole Populations

Data from the LIGO and Virgo detectors has confirmed that stellar-mass black holes can merge within a Hubble time, leaving behind massive remnant black holes. In some astrophysical environments such as globular clusters and AGN disks, it may be possible for these remnants to take part in further compact-object mergers, producing a population of hierarchically formed black holes. In this work, we present a parameterized framework for describing the population of binary black hole mergers, while self-consistently accounting for hierarchical mergers. The framework casts black holes as particles in a box which can collide based on an effective cross-section, but allows inputs from more detailed astrophysical simulations. Our approach is relevant to any population which is comprised of second or higher generation black holes, such as primordial black holes or dense cluster environments. We describe some possible inputs to this generic model and their effects on the black hole merger populations, and use the model to perform Bayesian inference on the catalog of black holes from LIGO and Virgo's first two observing runs. We find that models with a high rate of hierarchical mergers are disfavored, consistent with previous population analyses. Future gravitational-wave events will further constrain the inputs to this generic hierarchical merger model, enabling a deeper look into the formation environments of binary black holes.

preprint2020arXiv

Cosmological Inference using Gravitational Wave Standard Sirens: A Mock Data Challenge

The observation of binary neutron star merger GW170817, along with its optical counterpart, provided the first constraint on the Hubble constant $H_0$ using gravitational wave standard sirens. When no counterpart is identified, a galaxy catalog can be used to provide the necessary redshift information. However, the true host might not be contained in a catalog which is not complete out to the limit of gravitational-wave detectability. These electromagnetic and gravitational-wave selection effects must be accounted for. We describe and implement a method to estimate $H_0$ using both the counterpart and the galaxy catalog standard siren methods. We perform a series of mock data analyses using binary neutron star mergers to confirm our ability to recover an unbiased estimate of $H_0$. Our simulations used a simplified universe with no redshift uncertainties or galaxy clustering, but with different magnitude-limited catalogs and assumed host galaxy properties, to test our treatment of both selection effects. We explore how the incompleteness of catalogs affects the final measurement of $H_0$, as well as the effect of weighting each galaxy's likelihood of being a host by its luminosity. In our most realistic simulation, where the simulated catalog is about three times denser than the density of galaxies in the local universe, we find that a 4.4\% measurement precision can be reached using galaxy catalogs with 50\% completeness and $\sim 250$ binary neutron star detections with sensitivity similar to that of Advanced LIGO's second observing run.

preprint2020arXiv

Does Matter Matter? Using the mass distribution to distinguish neutron stars and black holes

Gravitational-wave detectors have opened a new window through which we can observe black holes (BHs) and neutron stars (NSs). Analyzing the 11 detections from LIGO/Virgo's first gravitational-wave catalog, GWTC-1, we investigate whether the power-law fit to the BH mass spectrum can also accommodate the binary neutron star (BNS) event GW170817, or whether we require an additional feature, such as a mass gap, in between the NS and BH populations. We find that with respect to the power-law fit to binary black hole (BBH) masses, GW170817 is an outlier at the 0.13\% level, suggesting a distinction between NS and BH masses. A single power-law fit across the entire mass range is in mild tension with: (a) the detection of one source in the BNS mass range ($\sim 1$--$2.5 \,M_\odot$), (b) the absence of detections in the "mass-gap" range ($\sim 2.5$--$5 \,M_\odot$), and (c) the detection of 10 sources in the BBH mass range ($\gtrsim 5 \,M_\odot$). Instead, the data favor models with a feature between NS and BH masses, including a mass gap (Bayes factor of 4.6) and a break in the power law, with a steeper slope at NS masses compared to BH masses (91\% credibility). We estimate the merger rates of compact binaries based on our fit to the global mass distribution, finding $\mathcal{R}_\mathrm{BNS} = 871^{+3015}_{-805} \ \mathrm{Gpc}^{-3} \ \mathrm{yr}^{-1}$ and $\mathcal{R}_\mathrm{BBH} = 47.5^{+57.9}_{-28.8} \ \mathrm{Gpc}^{-3} \ \mathrm{yr}^{-1}$. We conclude that, even in the absence of any prior knowledge of the difference between NSs and BHs, the gravitational-wave data alone already suggest two distinct populations of compact objects.

preprint2020arXiv

Picky Partners: The Pairing of Component Masses in Binary Black Hole Mergers

We examine the relationship between individual black hole (BH) masses in merging binary black hole (BBH) systems. Analyzing the ten BBH detections from LIGO/Virgo's first two observing runs, we find that the masses of the component BHs comprising each binary are unlikely to be randomly drawn from the same underlying distribution. Instead, the two BHs of a given binary prefer to be of comparable mass. We show that it is $\sim 5$ times more likely that the component BHs in a given binary are always equal (to within 5\%) than that they are randomly paired. If we assume that the probability of a merger between two BHs scales with the mass ratio $q$ as $q^β$, so that $β=0$ corresponds to random pairings, we find $β>0$ is favored at credibility 0.987. By modeling the mass distribution, we find that the median mass ratio is $q_{50\%} = 0.91^{+0.05}_{-0.17}$ at 90\% credibility. While the pairing between BHs depends on their mass ratio, we find no evidence that it depends on the total mass of the system: it is $\sim 6$ times more likely that the pairing depends purely on the mass ratio than on the total mass. We predict that 99\% of BBHs detected by LIGO/Virgo will have mass ratios $q > 0.5$. We conclude that merging black holes do not form random pairings; instead they are selective about their partners, preferring to mate with black holes of a similar mass. The details of these selective pairings provide insight into the underlying formation channels of merging binaries.

preprint2019arXiv

Counting on Short Gamma-Ray Bursts: Gravitational-Wave Constraints of Jet Geometry

The detection of GW170817 in gravitational waves and gamma rays revealed that short gamma-ray bursts are associated with the merger of neutron-stars. Gamma rays are thought to result from the formation of collimated jets, but the details of this process continue to elude us. One fundamental observable is the emission profile of the jet as a function of viewing angle. We present two methods to measure the effective angular width, $θ_B$, of short gamma-ray burst (sGRB) jets using gravitational wave and gamma-ray data, assuming all sGRBs have the same angular dependence for their luminosities. The first is a counting experiment, where we combine the known detection thresholds of the LIGO/Virgo and Fermi Gamma Ray Burst Monitor detectors to infer parameters of systems that are detected in gravitational waves. This method requires minimal knowledge about each event, beyond whether or not they were detected in gamma-rays. The second method uses additional information from the gravitational-wave and electromagnetic data to estimate parameters of the source, and thereby improve constraints on jet properties. Applying our methods to GW170817, we find only weak constraints on the sGRB luminosity profile, with statistical uncertainty dominating differences between models. We also analyze simulated events from future observing runs, and find that with 5 and 100 BNS detections, the counting method constrains the relative uncertainty in $θ_B$ to within 51% and 12%, respectively. Incorporating gravitational-wave parameter estimation would further tighten these constraints to 43% and 9.6%. In the limit of many detections, incorporating parameter estimation achieves only marginal improvements; we conclude that the majority of the information about jet structure comes from the relative sensitivities of gravitational-wave and gamma-ray detectors as encoded in simple counting experiments.

preprint2019arXiv

Nonparametric Inference of Neutron Star Composition, Equation of State, and Maximum Mass with GW170817

The detection of GW170817 in gravitational waves provides unprecedented constraints on the equation of state (EOS) of the ultra-dense matter within the cores of neutron stars (NSs). We extend the nonparametric analysis first introduced in Landry & Essick (2019), and confirm that GW170817 favors soft EOSs. We infer macroscopic observables for a canonical 1.4 $M_{\odot}$ NS, including the tidal deformability $Λ_{1.4} = 211^{+312}_{-137}$ ($491^{+216}_{-181}$) and radius $R_{1.4}= 10.86^{+2.04}_{-1.42}$ ($12.51^{+1.00}_{-0.88}$) km, as well as the maximum mass for nonrotating NSs, $M_{max} = 2.064^{+0.260}_{-1.34}$ ($2.017^{0.238}_{-0.087}$) $M_\odot$, with nonparametric priors loosely (tightly) constrained to resemble candidate EOSs from the literature. Furthermore, we find weak evidence that GW170817 involved at least one NS based on gravitational-wave data alone ($B^{NS}_{BBH}= 3.3 \pm 1.4$), consistent with the observation of electromagnetic counterparts. We also investigate GW170817's implications for the maximum spin frequency of millisecond pulsars, and find that the fastest known pulsar is spinning at more than 50% of its breakup frequency at 90% confidence. We additionally find modest evidence in favor of quark matter within NSs, and GW170817 favors the presence of at least one disconnected hybrid star branch in the mass--radius relation over a single stable branch by a factor of 2. Assuming there are multiple stable branches, we find a suggestive posterior preference for a sharp softening around nuclear density followed by stiffening around twice nuclear density, consistent with a strong first-order phase transition. While the statistical evidence in favor of new physics within NS cores remains tenuous with GW170817 alone, these tantalizing hints reemphasize the promise of gravitational waves for constraining the supranuclear EOS.

preprint2019arXiv

The Most Massive Binary Black Hole Detections and the Identification of Population Outliers

Advanced LIGO and Virgo detected ten binary black holes (BBHs) in their first two observing runs (O1 and O2). Analysis of these events found strong evidence for a dearth of BBHs with component masses greater than $\sim45 \ M_\odot$, as would be expected from a pair-instability mass gap. Meanwhile, a standalone analysis of the merger GW170729 found its primary mass $m_1 = {51.2^{+16.2}_{-11.0} \ M_\odot}$, with the majority of its posterior support at $m_1 > 45 \ M_\odot$. Although this appears to be in contradiction with the existence of a limit at $\sim45\ M_\odot$, we argue that individual events cannot be evaluated without reference to the entire population. When GW170729 is analyzed jointly with the rest of the detections, as part of a full hierarchical population analysis, its inferred primary mass tightens considerably, to $m_1 = {38.9^{+7.3}_{-4.5} \ M_\odot}$. For a large sample of events in the presence of noise, apparent outliers in the detected distribution are inevitable, even if the underlying population forbids outliers. We discuss methods of distinguishing between statistical fluctuations and population outliers using posterior predictive tests. Applying these tests to the primary mass distribution in O1 and O2, we find that the ten detections are consistent with even the simplest power-law plus maximum-mass model considered by the LVC. This supports the claim that GW170729 is not a population outlier. We also provide non-parametric constraints on the rate of high-mass mergers and conservatively bound the rate of mergers with $m_1 > 45 \ M_\odot$ at $2.8^{+5.4}_{-2.0}\%$ of the total merger rate. After 100 detections like those of O1 and O2 from a population with a maximum primary mass of $45 \, M_\odot$, it would be common for the most massive system to have an observed maximum-likelihood mass $m_1 \gtrsim 70 \, M_\odot$.

preprint2016arXiv

Compact Binary Merger Rates: Comparison with LIGO/Virgo Upper Limits

We compare evolutionary predictions of double compact object merger rate densities with initial and forthcoming LIGO/Virgo upper limits. We find that: (i) Due to the cosmological reach of advanced detectors, current conversion methods of population synthesis predictions into merger rate densities are insufficient. (ii) Our optimistic models are a factor of 18 below the initial LIGO/Virgo upper limits for BH-BH systems, indicating that a modest increase in observational sensitivity (by a factor of 2.5) may bring the first detections or first gravitational wave constraints on binary evolution. (iii) Stellar-origin massive BH-BH mergers should dominate event rates in advanced LIGO/Virgo and can be detected out to redshift z=2 with templates including inspiral, merger, and ringdown. Normal stars (<150 Msun) can produce such mergers with total redshifted mass up to 400 Msun. (iv) High black hole natal kicks can severely limit the formation of massive BH-BH systems (both in isolated binary and in dynamical dense cluster evolution), and thus would eliminate detection of these systems even at full advanced LIGO/Virgo sensitivity. We find that low and high black hole natal kicks are allowed by current observational electromagnetic constraints. (v) The majority of our models yield detections of all types of mergers (NS-NS, BH-NS, BH-BH) with advanced detectors. Numerous massive BH-BH merger detections will indicate small (if any) natal kicks for massive BHs.

preprint2016arXiv

Finding the One: Identifying the Host Galaxies of Gravitational-Wave Sources

We explore the localization of compact binary coalescences with ground-based gravitational-wave detector networks. We simulate tens of thousands of binary events, and present the distributions of localization sky areas and localization volumes for a range of sources and network configurations. We show that generically there exists a tail of particularly well-localized events, with 2D and 3D localizations of $<10\,\mbox{deg}^2$ and $<1000\,\mbox{Mpc}^3$ achievable, respectively, starting in LIGO/Virgo's third observing run. Incorporating estimates for the galaxy density and the binary event rates, we argue that future gravitational-wave detector networks will localize a small number of binary systems per year to a sufficiently small volume that the unique host galaxy might be identified. For these golden events, which are generally the closest and loudest ones, the gravitational-wave detector networks will point (in 3D; the length of the finger matters) directly at the source. This will allow for studies of the properties of the host galaxies of compact binary mergers, which may be an important component in exploring the formation channels of these sources. In addition, since the host will provide an independent measurement of the redshift, this will allow the use of the event as a standard siren to measure cosmology. Furthermore, identification of a small number of host galaxies can enable deep follow-up searches for associated electromagnetic transients.

preprint2016arXiv

The first gravitational-wave source from the isolated evolution of two 40-100 Msun stars

The merger of two massive 30 Msun black holes has been detected in gravitational waves (1,GW150914). This discovery validates recent predictions (2-4) that massive binary black holes would constitute the first detection. However, previous calculations have not sampled the relevant binary black hole progenitors---massive, low-metallicity binary stars---with sufficient accuracy and input physics to enable robust predictions to better than several orders of magnitude (5-10). Here we report a suite of high-precision numerical simulations of binary black hole formation via the evolution of isolated binary stars, providing a framework to interpret GW150914 and predict the properties of subsequent binary black hole gravitational-wave events. Our models imply that these events form in an environment where the metallicity is less than 10 percent of solar; have initial masses of 40-100 Msun; and interact through mass transfer and a common envelope phase. Their progenitors likely form either at 2 Gyr, or somewhat less likely, at 11 Gyr after the Big Bang. Most binary black holes form without supernova explosions, and their spins are nearly unchanged since birth, but do not have to be parallel. The classical field formation of binary black holes proposed in this study, with low natal kicks and restricted common envelope evolution, produces 40 times more binary black holes than dynamical formation channels involving globular clusters (11) and is comparable to the rate from homogeneous evolution channels (12-15). Our calculations predict detections of about 1,000 black hole mergers per year with total mass of 20-80 Msun once second generation ground-based gravitational wave observatories reach full sensitivity.

preprint2015arXiv

Neutron stars versus black holes: probing the mass gap with LIGO/Virgo

Inspirals and mergers of black hole (BHs) and/or neutron star (NSs) binaries are expected to be abundant sources for ground-based gravitational-wave (GW) detectors. We assess the capabilities of Advanced LIGO and Virgo to measure component masses using inspiral waveform models including spin-precession effects using a large ensemble of GW sources {\bf randomly oriented and distributed uniformly in volume. For 1000 sources this yields signal-to-noise ratios between 7 and 200}. We make quantitative predictions for how well LIGO and Virgo will distinguish between BHs and NSs and appraise the prospect of using LIGO/Virgo observations to definitively confirm, or reject, the existence of a putative "mass gap" between NSs ($m\leq3\ M_\odot$) and BHs ($m\geq 5\ M_\odot$). We find sources with the smaller mass component satisfying $m_2 \lesssim1.5\ M_\odot$ to be unambiguously identified as containing at least one NS, while systems with $m_2\gtrsim6\ M_\odot$ will be confirmed binary BHs. Binary BHs with $m_2<5\ M_\odot$ (i.e., in the gap) cannot generically be distinguished from NSBH binaries. High-mass NSs ($2<m<3$ $M_\odot$) are often consistent with low-mass BH ($m<5\ M_\odot$), posing a challenge for determining the maximum NS mass from LIGO/Virgo observations alone. Individual sources will seldom be measured well enough to confirm objects in the mass gap and statistical inferences drawn from the detected population will be strongly dependent on the underlying distribution. If nature happens to provide a mass distribution with the populations relatively cleanly separated in chirp mass space, as some population synthesis models suggest, then NSs and BHs are more easily distinguishable.

preprint2014arXiv

Comparison of LIGO/Virgo upper limits with predicted compact binary merger rates

We compare the current LIGO/Virgo upper limits on double compact object volumetric merger rates with our theoretical predictions. Our optimistic models are a factor of 3 below the existing upper limits for massive BH-BH systems with total mass 50-70 Msun, suggesting that a small increase in observational sensitivity may bring the first detections. The LIGO/Virgo gravitational wave observatories are currently being upgraded to advanced design sensitivity. If a sizeable population of BH-BH binaries is detected, the maximum total binary mass of this population will discriminate between two general families of common envelope models. If no binaries are detected, the new upper limits will provide astrophysically useful information about the environment and physical processes (e.g., metallicity of host galaxies or BH natal kicks) crucial to the formation of binaries containing black holes. For NS-NS systems, our predicted rates are 3 orders of magnitude below the current upper limits; even if advanced instruments reach their design sensitivities (factors of 10 times in distance, and 1,000 times in volumetric rate) the detection of NS-NS systems is not assured. However, we note that although our predicted NS-NS merger rates are consistent with estimates derived from Galactic NS-NS binaries and short GRBs, they are on the low side of these empirical estimates.

preprint2014arXiv

Dark Sky Simulations: Early Data Release

The Dark Sky Simulations are an ongoing series of cosmological N-body simulations designed to provide a quantitative and accessible model of the evolution of the large-scale Universe. Such models are essential for many aspects of the study of dark matter and dark energy, since we lack a sufficiently accurate analytic model of non-linear gravitational clustering. In July 2014, we made available to the general community our early data release, consisting of over 55 Terabytes of simulation data products, including our largest simulation to date, which used $1.07 \times 10^{12}~(10240^3)$ particles in a volume $8h^{-1}\mathrm{Gpc}$ across. Our simulations were performed with 2HOT, a purely tree-based adaptive N-body method, running on 200,000 processors of the Titan supercomputer, with data analysis enabled by yt. We provide an overview of the derived halo catalogs, mass function, power spectra and light cone data. We show self-consistency in the mass function and mass power spectrum at the 1% level over a range of more than 1000 in particle mass. We also present a novel method to distribute and access very large datasets, based on an abstraction of the World Wide Web (WWW) as a file system, remote memory-mapped file access semantics, and a space-filling curve index. This method has been implemented for our data release, and provides a means to not only query stored results such as halo catalogs, but also to design and deploy new analysis techniques on large distributed datasets.

preprint2014arXiv

The Formation and Gravitational-Wave Detection of Massive Stellar Black-Hole Binaries

If binaries consisting of two 100 Msun black holes exist they would serve as extraordinarily powerful gravitational-wave sources, detectable to redshifts of z=2 with the advanced LIGO/Virgo ground-based detectors. Large uncertainties about the evolution of massive stars preclude definitive rate predictions for mergers of these massive black holes. We show that rates as high as hundreds of detections per year, or as low as no detections whatsoever, are both possible. It was thought that the only way to produce these massive binaries was via dynamical interactions in dense stellar systems. This view has been challenged by the recent discovery of several stars with mass above 150 Msun in the R136 region of the Large Magellanic Cloud. Current models predict that when stars of this mass leave the main sequence, their expansion is insufficient to allow common envelope evolution to efficiently reduce the orbital separation. The resulting black-hole--black-hole binary remains too wide to be able to coalesce within a Hubble time. If this assessment is correct, isolated very massive binaries do not evolve to be gravitational-wave sources. However, other formation channels exist. For example, the high multiplicity of massive stars, and their common formation in relatively dense stellar associations, opens up dynamical channels for massive black hole mergers (e.g., via Kozai cycles or repeated binary-single interactions). We identify key physical factors that shape the population of very massive black-hole--black-hole binaries. Advanced gravitational-wave detectors will provide important constraints on the formation and evolution of very massive stars.

preprint2014arXiv

The Loudest Gravitational Wave Events

As first emphasized by Bernard Schutz, there exists a universal distribution of signal-to-noise ratios for gravitational wave detection. Because gravitational waves (GWs) are almost impossible to obscure via dust absorption or other astrophysical processes, the strength of the detected signal is dictated solely by the emission strength and the distance to the source. Assuming that the space density of an arbitrary population of GW sources does not evolve, we show explicitly that the distribution of detected signal-to-noise (SNR) values depends solely on the detection threshold; it is independent of the detector network (interferometer or pulsar timing array), the individual detector noise curves (initial or Advanced LIGO), the nature of the GW sources (compact binary coalescence, supernova, or some other discrete source), and the distributions of source variables (only non-spinning neutron stars of mass exactly $1.4\,M_\odot$ or a complicated distribution of masses and spins). We derive the SNR distribution for each individual detector within a network as a function of the relative detector orientations and sensitivities. While most detections will have SNR near the detection threshold, there will be a tail of events to higher SNR. We derive the SNR distribution of the loudest (highest SNR) events in any given sample of detections. We find that the median SNR of the loudest out of the first four events should have an $\mbox{SNR}=22$ (for a threshold of 12, appropriate for the Advanced LIGO/Virgo network), increasing to a median value for the loudest SNR of 47 for 40 detections. We expect these loudest events to provide particularly powerful constraints on their source parameters, and they will play an important role in extracting astrophysics from gravitational wave sources. These distributions also offer an important internal calibration of the response of the GW detector networks.

preprint2013arXiv

Determining the Hubble constant from gravitational wave observations of merging compact binaries

Recent observations have accumulated compelling evidence that some short gamma-ray bursts (SGRBs) are associated with the mergers of neutron star (NS) binaries. This would indicate that the SGRB event is associated with a gravitational-wave (GW) signal corresponding to the final inspiral of the compact binary. In addition, the radioactive decay of elements produced in NS binary mergers may result in transients visible in the optical and infrared with peak luminosities on hours-days timescales. Simultaneous observations of the inspiral GWs and signatures in the electromagnetic band may allow us to directly and independently determine both the luminosity distance and redshift to a binary. These standard sirens (the GW analog of standard candles) have the potential to provide an accurate measurement of the low-redshift Hubble flow. In addition, these systems are absolutely calibrated by general relativity, and therefore do not experience the same set of astrophysical systematics found in traditional standard candles, nor do the measurements rely on a distance ladder. We show that 15 observable GW and EM events should allow the Hubble constant to be measured with 5% precision using a network of detectors that includes advanced LIGO and Virgo. Measuring 30 beamed GW-SGRB events could constrain H_0 to better than 1%. When comparing to standard Gaussian likelihood analysis, we find that each event's non-Gaussian posterior in H_0 helps reduce the overall measurement errors in H_0 for an ensemble of NS binary mergers.

preprint2013arXiv

Double Compact Objects II: Cosmological Merger Rates

The development of advanced gravitational wave (GW) observatories, such as Advanced LIGO and Advanced Virgo, provides impetus to refine theoretical predictions for what these instruments might detect. In particular, with the range increasing by an order of magnitude, the search for GW sources is extending beyond the "local" Universe and out to cosmological distances. Double compact objects (neutron star-neutron star (NS-NS), black hole-neutron star (BH-NS) and black hole-black hole (BH-BH) systems) are considered to be the most promising gravitational wave sources. In addition, NS-NS and/or BH-NS systems are thought to be the progenitors of gamma ray bursts (GRBs), and may also be associated with kilonovae. In this paper we present the merger event rates of these objects as a function of cosmological redshift. We provide the results for four cases, each one investigating a different important evolution parameter of binary stars. Each case is also presented for two metallicity evolution scenarios. We find that (i) in most cases NS-NS systems dominate the merger rates in the local Universe, while BH-BH mergers dominate at high redshift; (ii) BH-NS mergers are less frequent than other sources per unit volume, for all time; and (iii) natal kicks may alter the observable properties of populations in a significant way, allowing the underlying models of binary evolution and compact object formation to be easily distinguished. This is the second paper in a series of three. The third paper will focus on calculating the detection rates of mergers by gravitational wave telescopes.

preprint2013arXiv

Finding the First Cosmic Explosions I: Pair-Instability Supernovae

The first stars are the key to the formation of primitive galaxies, early cosmological reionization and chemical enrichment, and the origin of supermassive black holes. Unfortunately, in spite of their extreme luminosities, individual Population III stars will likely remain beyond the reach of direct observation for decades to come. However, their properties could be revealed by their supernova explosions, which may soon be detected by a new generation of NIR observatories such as JWST and WFIRST. We present light curves and spectra for Pop III pair-instability supernovae calculated with the Los Alamos radiation hydrodynamics code RAGE. Our numerical simulations account for the interaction of the blast with realistic circumstellar envelopes, the opacity of the envelope, and Lyman absorption by the neutral IGM at high redshift, all of which are crucial to computing the NIR signatures of the first cosmic explosions. We find that JWST will detect pair-instability supernovae out to z > 30, WFIRST will detect them in all-sky surveys out to z ~ 15 - 20 and LSST and Pan-STARRS will find them at z ~ 7 - 8. The discovery of these ancient explosions will probe the first stellar populations and reveal the existence of primitive galaxies that might not otherwise have been detected.

preprint2013arXiv

Finding the First Cosmic Explosions II: Core-Collapse Supernovae

Understanding the properties of Pop III stars is prerequisite to elucidating the nature of primeval galaxies, the chemical enrichment and reionization of the early IGM, and the origin of supermassive black holes. While the primordial IMF remains unknown, recent evidence from numerical simulations and stellar archaeology suggests that some Pop III stars may have had lower masses than previously thought, 15 - 50 \Ms in addition to 50 - 500 \Ms. The detection of Pop III supernovae by JWST, WFIRST or the TMT could directly probe the primordial IMF for the first time. We present numerical simulations of 15 - 40 \Ms Pop III core-collapse SNe done with the Los Alamos radiation hydrodynamics code RAGE. We find that they will be visible in the earliest galaxies out to z ~ 10 - 15, tracing their star formation rates and in some cases revealing their positions on the sky. Since the central engines of Pop III and solar-metallicity core-collapse SNe are quite similar, future detection of any Type II supernovae by next-generation NIR instruments will in general be limited to this epoch.

preprint2013arXiv

Illuminating the Primeval Universe with Type IIn Supernovae

The detection of Pop III supernovae could directly probe the primordial IMF for the first time, unveiling the properties of the first galaxies, early chemical enrichment and reionization, and the seeds of supermassive black holes. Growing evidence that some Pop III stars were less massive than 100 solar masses may complicate prospects for their detection, because even though they would have been more plentiful they would have died as core-collapse supernovae, with far less luminosity than pair-instability explosions. This picture greatly improves if the SN shock collides with a dense circumstellar shell ejected during a prior violent LBV type eruption. Such collisions can turn even dim SNe into extremely bright ones whose luminosities can rival those of pair-instability SNe. We present simulations of Pop III Type IIn SN light curves and spectra performed with the Los Alamos RAGE and SPECTRUM codes. Taking into account Lyman-alpha absorption in the early universe and cosmological redshifting, we find that 40 solar mass Pop III Type IIn SNe will be visible out to z ~ 20 with JWST and out to z ~ 7 with WFIRST. Thus, even low mass Pop III SNe can be used to probe the primeval universe.

preprint2013arXiv

Supermassive Seeds for Supermassive Black Holes

Recent observations of quasars powered by supermassive black holes (SMBHs) out to z > 7 constrain both the initial seed masses and the growth of the most massive black holes (BHs) in the early universe. Here we elucidate the implications of the radiative feedback from early generations of stars and from BH accretion for popular models for the formation and growth of seed BHs. We show that by properly accounting for (1) the limited role of mergers in growing seed BHs as inferred from cosmological simulations of early star formation and radiative feedback, (2) the sub-Eddington accretion rates of BHs expected at the earliest times, and (3) the large radiative efficiencies (e_rad) of the most massive BHs inferred from observations of active galactic nuclei at high redshift (e_rad > 0.1), we are led to the conclusion that the initial BH seeds may have been as massive as > 10^5 M_Sun. This presents a strong challenge to the Population III seed model, which calls for seed masses of ~ 100 M_Sun and, even with constant Eddington-limited accretion, requires e_rad < 0.09 to explain the highest-z SMBHs in today's standard LambdaCDM cosmological model. It is, however, consistent with the prediction of the direct collapse scenario of SMBH seed formation, in which a supermassive primordial star forms in a region of the universe with a high molecule-dissociating background radiation field, and collapses directly into a 10^4--10^6 M_Sun seed BH. These results corroborate recent cosmological simulations and observational campaigns which suggest that these massive BHs were the seeds of a large fraction of the SMBHs residing in the centers of galaxies today.

preprint2012arXiv

Evidence for Type Ia Supernova Diversity from Ultraviolet Observations with the Hubble Space Telescope

We present ultraviolet (UV) spectroscopy and photometry of four Type Ia supernovae (SNe 2004dt, 2004ef, 2005M, and 2005cf) obtained with the UV prism of the Advanced Camera for Surveys on the Hubble Space Telescope. This dataset provides unique spectral time series down to 2000 Angstrom. Significant diversity is seen in the near maximum-light spectra (~ 2000--3500 Angstrom) for this small sample. The corresponding photometric data, together with archival data from Swift Ultraviolet/Optical Telescope observations, provide further evidence of increased dispersion in the UV emission with respect to the optical. The peak luminosities measured in uvw1/F250W are found to correlate with the B-band light-curve shape parameter dm15(B), but with much larger scatter relative to the correlation in the broad-band B band (e.g., ~0.4 mag versus ~0.2 mag for those with 0.8 < dm15 < 1.7 mag). SN 2004dt is found as an outlier of this correlation (at > 3 sigma), being brighter than normal SNe Ia such as SN 2005cf by ~0.9 mag and ~2.0 mag in the uvw1/F250W and uvm2/F220W filters, respectively. We show that different progenitor metallicity or line-expansion velocities alone cannot explain such a large discrepancy. Viewing-angle effects, such as due to an asymmetric explosion, may have a significant influence on the flux emitted in the UV region. Detailed modeling is needed to disentangle and quantify the above effects.

preprint2012arXiv

GRB beaming and gravitational-wave observations

Using the observed rate of short-duration gamma-ray bursts (GRBs) it is possible to make predictions for the detectable rate of compact binary coalescences in gravitational-wave detectors. These estimates rely crucially on the growing consensus that short gamma-ray bursts are associated with the merger of two neutron stars or a neutron star and a black hole, but otherwise make no assumptions beyond the observed rate of short GRBs. In particular, our results do not assume coincident gravitational wave and electromagnetic observations. We show that the non-detection of mergers in the existing LIGO/Virgo data constrains the progenitor masses and beaming angles of gamma-ray bursts. For future detectors, we find that the first detection of a NS-NS binary coalescence associated with the progenitors of short GRBs is likely to happen within the first 16 months of observation, even in the case of a modest network of observatories (e.g., only LIGO-Hanford and LIGO-Livingston) operating at modest sensitivities (e.g., advanced LIGO design sensitivity, but without signal recycling mirrors), and assuming a conservative distribution of beaming angles (e.g. all GRBs beamed at θ=30 deg). Less conservative assumptions reduce the waiting time until first detection to weeks to months. Alternatively, the compact binary coalescence model of short GRBs can be ruled out if a binary is not seen within the first two years of operation of a LIGO-Hanford, LIGO-Livingston, and Virgo network at advanced design sensitivity. We also demonstrate that the rate of GRB triggered sources is less than the rate of untriggered events if θ<30 deg, independent of the noise curve, network configuration, and observed GRB rate. Thus the first detection in GWs of a binary GRB progenitor is unlikely to be associated with a GRB.

preprint2012arXiv

Seeing the First Supernovae at the Edge of the Universe with JWST

The first stars ended the cosmic Dark Ages and created the first heavy elements necessary for the formation of planets and life. The properties of these stars remain uncertain, and it may be decades before individual Pop III stars are directly observed. Their masses, however, can be inferred from their supernova explosions, which may soon be found in both deep-field surveys by JWST and in all-sky surveys by WFIRST. We have performed radiation hydrodynamical simulations of the near infrared signals of Pop III pair-instability supernovae in realistic circumstellar environments with Lyman absorption by the neutral intergalactic medium. We find that JWST and WFIRST will detect these explosions out to z ~ 30 and 20, respectively, unveiling the first generation of stars in the universe.

preprint2011arXiv

Localizing compact binary inspirals on the sky using ground-based gravitational wave interferometers

The inspirals and mergers of compact binaries are among the most promising events for ground-based gravitational-wave (GW) observatories. The detection of electromagnetic (EM) signals from these sources would provide complementary information to the GW signal. It is therefore important to determine the ability of gravitational-wave detectors to localize compact binaries on the sky, so that they can be matched to their EM counterparts. We use Markov Chain Monte Carlo techniques to study sky localization using networks of ground-based interferometers. Using a coherent-network analysis, we find that the Laser Interferometer Gravitational Wave Observatory (LIGO)-Virgo network can localize 50% of their ~8 sigma detected neutron star binaries to better than 50 sq.deg. with 95% confidence region. The addition of the Large Scale Cryogenic Gravitational Wave Telescope (LCGT) and LIGO-Australia improves this to 12 sq.deg.. Using a more conservative coincident detection threshold, we find that 50% of detected neutron star binaries are localized to 13 sq.deg. using the LIGO-Virgo network, and to 3 sq.deg. using the LIGO-Virgo-LCGT-LIGO-Australia network. Our findings suggest that the coordination of GW observatories and EM facilities offers great promise.

preprint2011arXiv

Outflowing Galactic Winds in Post-starburst and AGN Host Galaxies at 0.2<z<0.8

We present Keck/LRIS-B spectra for a sample of ten AEGIS X-ray AGN host galaxies and thirteen post-starburst galaxies from SDSS and DEEP2 at 0.2<z<0.8 in order to investigate the presence, properties, and influence of outflowing galactic winds at intermediate redshifts. We focus on galaxies that either host a low-luminosity AGN or have recently had their star formation quenched to test whether these galaxies have winds of sufficient velocity to potentially clear gas from the galaxy. We find, using absorption features of Fe II, Mg II, and Mg I, that six of the ten (60%) X-ray AGN host galaxies and four of the thirteen (31%) post-starburst galaxies have outflowing galactic winds, with typical velocities of ~200 km/s. We additionally find that most of the galaxies in our sample show line emission, possibly from the wind, in either Fe II* or Mg II. A total of 100% of our X-ray AGN host sample (including four red sequence galaxies) and 77% of our post-starburst sample has either blueshifted absorption or line emission. Several K+A galaxies have small amounts of cool gas absorption at the systemic velocity, indicating that not all of the cool gas has been expelled. We conclude that while outflowing galactic winds are common in both X-ray low-luminosity AGN host galaxies and post-starburst galaxies at intermediate redshifts, the winds are likely driven by supernovae (as opposed to AGN) and do not appear to have sufficiently high velocities to quench star formation in these galaxies.

preprint2010arXiv

Exploring short gamma-ray bursts as gravitational-wave standard sirens

Recent observations support the hypothesis that a large fraction of "short-hard" gamma-ray bursts (SHBs) are associated with compact binary inspiral. Since gravitational-wave (GW) measurements of well-localized inspiraling binaries can measure absolute source distances, simultaneous observation of a binary's GWs and SHB would allow us to independently determine both its luminosity distance and redshift. Such a "standard siren" (the GW analog of a standard candle) would provide an excellent probe of the relatively nearby universe's expansion, complementing other standard candles. In this paper, we examine binary measurement using a Markov Chain Monte Carlo technique to build the probability distributions describing measured parameters. We assume that each SHB observation gives both sky position and the time of coalescence, and we take both binary neutron stars and black hole-neutron star coalescences as plausible SHB progenitors. We examine how well parameters particularly distance) can be measured from GW observations of SHBs by a range of ground-based detector networks. We find that earlier estimates overstate how well distances can be measured, even at fairly large signal-to-noise ratio. The fundamental limitation to determining distance proves to be a degeneracy between distance and source inclination. Overcoming this limitation requires that we either break this degeneracy, or measure enough sources to broadly sample the inclination distribution. (Abridged)

preprint2010arXiv

Reducing the weak lensing noise for the gravitational wave Hubble diagram using the non-Gaussianity of the magnification distribution

Gravitational wave sources are a promising cosmological standard candle because their intrinsic luminosities are determined by fundamental physics (and are insensitive to dust extinction). They are, however, affected by weak lensing magnification due to the gravitational lensing from structures along the line of sight. This lensing is a source of uncertainty in the distance determination, even in the limit of perfect standard candle measurements. It is commonly believed that the uncertainty in the distance to an ensemble of gravitational wave sources is limited by the standard deviation of the lensing magnification distribution divided by the square root of the number of sources. Here we show that by exploiting the non-Gaussian nature of the lensing magnification distribution, we can improve this distance determination, typically by a factor of 2--3; we provide a fitting formula for the effective distance accuracy as a function of redshift for sources where the lensing noise dominates.

preprint2010arXiv

The most massive objects in the Universe

We calculate the most massive object in the Universe, finding it to be a cluster of galaxies with total mass M_200=3.8e15 Msun at z=0.22, with the 1 sigma marginalized regions being 3.3e15 Msun<M<4.4e15 Msun and 0.12<z<0.36. We restrict ourselves to self-gravitating bound objects, and base our results on halo mass functions derived from N-body simulations. Since we consider the very highest mass objects, the number of candidates is expected to be small, and therefore each candidate can be extensively observed and characterized. If objects are found with excessively large masses, or insufficient objects are found near the maximum expected mass, this would be a strong indication of the failure of LambdaCDM. The expected range of the highest masses is very sensitive to redshift, providing an additional evolutionary probe of LambdaCDM. We find that the three most massive clusters in the recent SPT 178 deg^2 catalog match predictions, while XMMU J2235.3--2557 is roughly 3 sigma inconsistent with LambdaCDM. We discuss Abell 2163 and Abell 370 as candidates for the most massive cluster in the Universe, although uncertainties in their masses preclude definitive comparisons with theory. Our findings motivate further observations of the highest mass end of the mass function. Future surveys will explore larger volumes, and the most massive object in the Universe may be identified within the next decade. The mass distribution of the largest objects in the Universe is a potentially powerful test of LambdaCDM, probing non-Gaussianity and the behavior of gravity on large scales.

preprint2010arXiv

Weak lensing and dark energy: the impact of dark energy on nonlinear dark matter clustering

We examine the influence of percent-level dark energy corrections to the nonlinear matter power spectrum on constraints of the dark energy equation of state from future weak lensing probes. We explicitly show that a poor approximation (off by > 10%) to the nonlinear corrections causes a > 1 sigma bias on the determination of the dark energy equation of state. Future weak lensing surveys must therefore incorporate dark energy modifications to the nonlinear matter power spectrum accurate to the percent-level, to avoid introducing significant bias in their measurements. For the WMAP5 cosmology, the more accurate power spectrum is more sensitive to dark energy properties, resulting in a factor of two improvement in dark energy equation of state constraints. We explore the complementary constraints on dark energy from future weak lensing and supernova surveys. A space-based, JDEM-like survey measures the equation of state in five independent redshift bins to ~10%, while this improves to ~5% for a wide-field ground-based survey like LSST. These constraints are contingent upon our ability to control weak lensing systematic uncertainties to the sub-percent level.

preprint2009arXiv

Astro2010 Decadal Survey Whitepaper: Coordinated Science in the Gravitational and Electromagnetic Skies

It is widely expected that the coming decade will witness the first direct detection of gravitational waves (GWs). The ground-based LIGO and Virgo GW observatories are being upgraded to advanced sensitivity, and are expected to observe a significant binary merger rate. The launch of The Laser Interferometer Space Antenna (LISA) would extend the GW window to low frequencies, opening new vistas on dynamical processes involving massive (M >~ 10^5 M_Sun) black holes. GW events are likely to be accompanied by electromagnetic (EM) counterparts and, since information carried electromagnetically is complementary to that carried gravitationally, a great deal can be learned about an event and its environment if it becomes possible to measure both forms of radiation in concert. Measurements of this kind will mark the dawn of trans-spectral astrophysics, bridging two distinct spectral bands of information. The aim of this whitepaper is to articulate future directions in both theory and observation that are likely to impact broad astrophysical inquiries of general interest. What will EM observations reflect on the nature and diversity of GW sources? Can GW sources be exploited as complementary probes of cosmology? What cross-facility coordination will expand the science returns of gravitational and electromagnetic observations?

preprint2009arXiv

No Evidence for Dark Energy Dynamics from a Global Analysis of Cosmological Data

We use a variant of principal component analysis to investigate the possible temporal evolution of the dark energy equation of state, w(z). We constrain w(z) in multiple redshift bins, utilizing the most recent data from Type Ia supernovae, the cosmic microwave background, baryon acoustic oscillations, the integrated Sachs-Wolfe effect, galaxy clustering, and weak lensing data. Unlike other recent analyses, we find no significant evidence for evolving dark energy; the data remains completely consistent with a cosmological constant. We also study the extent to which the time-evolution of the equation of state would be constrained by a combination of current- and future-generation surveys, such as Planck and the Joint Dark Energy Mission.

preprint2009arXiv

On The Origin Of The Highest Redshift Gamma-Ray Bursts

GRB 080913 and GRB 090423 are the most distant gamma-ray bursts (GRBs) known to-date, with spectroscopically determined redshifts of z=6.7 and z=8.1, respectively. The detection of bursts at this early epoch of the Universe significantly constrains the nature of GRBs and their progenitors. We perform population synthesis studies of the formation and evolution of early stars, and calculate the resulting formation rates of short and long-duration GRBs at high redshift. The peak of the GRB rate from Population II stars occurs at z=7 for a model with efficient/fast mixing of metals, while it is found at z=3 for an inefficient/slow metallicity evolution model. We show that in the redshift range 6<z<10 essentially all GRBs originate from Population II stars, regardless of metallicity evolution model. These stars (having small, but non-zero metallicity) are the most likely progenitors for both long GRBs (collapsars) and short GRBs (NS-NS or BH-NS mergers) at this epoch. Although the predicted intrinsic rates of long and short GRBs are similar at these high redshifts, observational selection effects lead to higher (factor of 10) observed rates for long GRBs. We conclude that the two recently observed high-z GRB events are most likely long GRBs originating from Population II collapsars.

preprint2009arXiv

Ultra-high precision cosmology from gravitational waves

We show that the Big Bang Observer (BBO), a proposed space-based gravitational-wave (GW) detector, would provide ultra-precise measurements of cosmological parameters. By detecting ~300,000 compact-star binaries, and utilizing them as standard sirens, BBO would determine the Hubble constant to 0.1%, and the dark energy parameters w_0 and w_a to ~0.01 and 0.1,resp. BBO's dark-energy figure-of-merit would be approximately an order of magnitude better than all other proposed dark energy missions. To date, BBO has been designed with the primary goal of searching for gravitational waves from inflation. To observe this inflationary background, BBO would first have to detect and subtract out ~300,000 merging compact-star binaries, out to z~5. It is precisely this foreground which would enable high-precision cosmology. BBO would determine the luminosity distance to each binary to ~percent accuracy. BBO's angular resolution would be sufficient to uniquely identify the host galaxy for most binaries; a coordinated optical/infrared observing campaign could obtain the redshifts. Combining the GW-derived distances and EM-derived redshifts for such a large sample of objects leads to extraordinarily tight constraints on cosmological parameters. Such ``standard siren'' measurements of cosmology avoid many of the systematic errors associated with other techniques. We also show that BBO would be an exceptionally powerful gravitational lensing mission, and we briefly discuss other astronomical uses of BBO.

preprint2008arXiv

Implications of Two Type Ia Supernova Populations for Cosmological Measurements

Recent work suggests that Type Ia supernovae (SNe) are composed of two distinct populations: prompt and delayed. By explicitly incorporating properties of host galaxies, it may be possible to target and eliminate systematic differences between these two putative populations. However, any resulting {\em post}-calibration shift in luminosity between the components will cause a redshift-dependent systematic shift in the Hubble diagram. Utilizing an existing sample of 192 SNe Ia, we find that the average luminosity difference between prompt and delayed SNe is constrained to be $(4.5 \pm 8.9)%$. If the absolute difference between the two populations is 0.025 mag, and this is ignored when fitting for cosmological parameters, then the dark energy equation of state (EOS) determined from a sample of 2300 SNe Ia is biased at $\sim1σ$. By incorporating the possibility of a two-population systematic, this bias can be eliminated. However, assuming no prior on the strength of the two-population effect, the uncertainty in the best-fit EOS is increased by a factor of 2.5, when compared to the equivalent sample with no underlying two-population systematic. To avoid introducing a bias in the EOS parameters, or significantly degrading the measurement accuracy, it is necessary to control the post-calibration luminosity difference between prompt and delayed SN populations to better than 0.025 mag.

preprint2008arXiv

Measuring dark energy spatial inhomogeneity with supernova data

The gravitational lensing distortion of distant sources by the large-scale distribution of matter in the Universe has been extensively studied. In contrast, very little is known about the effects due to the large-scale distribution of dark energy. We discuss the use of Type Ia supernovae as probes of the spatial inhomogeneity and anisotropy of dark energy. We show that a shallow, almost all-sky survey can limit rms dark energy fluctuations at the horizon scale down to a fractional energy density of ~10^-4

preprint2005arXiv

Precision Determination of the Mass Function of Dark Matter Halos

The predicted mass function of dark matter halos is essential in connecting observed galaxy cluster counts and models of galaxy clustering to the properties of the primordial density field. We determine the mass function in the concordance $Λ$CDM cosmology, as well as its uncertainty, using sixteen $1024^3$-particle nested-volume dark-matter simulations, spanning a mass range of over five orders of magnitude. Using the nested volumes and single-halo tests, we find and correct for a systematic error in the friends-of-friends halo-finding algorithm. We find a fitting form and full error covariance for the mass function that successfully describes the simulations' mass function and is well-behaved outside the simulations' resolutions. Estimated forecasts of uncertainty in cosmological parameters from future cluster count surveys have negligible contribution from remaining statistical uncertainties in the central cosmology multiplicity function. There exists a potentially non-negligible cosmological dependence (non-universality) of the halo multiplicity function.

preprint2005arXiv

Using gravitational-wave standard sirens

Gravitational waves (GWs) from supermassive binary black hole (BBH) inspirals are potentially powerful standard sirens (the GW analog to standard candles) (Schutz 1986, 2002). Because these systems are well-modeled, the space-based GW observatory LISA will be able to measure the luminosity distance (but not the redshift) to some distant massive BBH systems with 1-10% accuracy. This accuracy is largely limited by pointing error: GW sources generally are poorly localized on the sky. Localizing the binary independently (e.g., through association with an electromagnetic counterpart) greatly reduces this positional error. An electromagnetic counterpart may also allow determination of the event's redshift. In this case, BBH coalescence would constitute an extremely precise (better than 1%) standard candle visible to high redshift. In practice, gravitational lensing degrades this precision, though the candle remains precise enough to provide useful information about the distance-redshift relation. Even if very rare, these GW standard sirens would complement, and increase confidence in, other standard candles.

preprint2004arXiv

Consequences of gravitational radiation recoil

Coalescing binary black holes experience an impulsive kick due to anisotropic emission of gravitational waves. We discuss the dynamical consequences of the recoil accompanying massive black hole mergers. Recoil velocities are sufficient to eject most coalescing black holes from dwarf galaxies and globular clusters, which may explain the apparent absence of massive black holes in these systems. Ejection from giant elliptical galaxies would be rare, but coalescing black holes are displaced from the center and fall back on a time scale of order the half-mass crossing time. Displacement of the black holes transfers energy to the stars in the nucleus and can convert a steep density cusp into a core. Radiation recoil calls into question models that grow supermassive black holes from hierarchical mergers of stellar-mass precursors.

preprint2004arXiv

Gravitational Waves from Stellar Collapse: Correlations to Explosion Asymmetries

The collapse of massive stars not only produces observable outbursts across the entire electromagnetic spectrum but, for Galactic (or near-Galactic) supernovae, detectable signals for ground-based neutrino and gravitational wave detectors. Gravitational waves and neutrinos provide the only means to study the actual engine behind the optical outbursts: the collapsed stellar core. While the neutrinos are most sensitive to details of the equation of state, gravitational waves provide a means to study the mass asymmetries in this central core. We present gravitational wave signals from a series of 3-dimensional core-collapse simulations with asymmetries derived from initial perturbations caused by pre-collapse convection, core rotation, and low-mode convection in the explosion engine itself. A Galactic supernovae will allow us to differentiate these different sources of asymmetry. Combining this signal with other observations of the supernova, from neutrinos to gamma-rays to the compact remnant, dramatically increases the predictive power of the gravitational wave signal. We conclude with a discussion of the gravitational wave signal arising from collapsars, the leading engine for long-duration gamma-ray bursts.

preprint2004arXiv

Safety in numbers: Gravitational Lensing Degradation of the Luminosity Distance-Redshift Relation

Observation of the expansion history of the Universe allows exploration of the physical properties and energy density of the Universe's various constituents. Standardizable candles such as Type Ia supernovae remain one of the most promising and robust tools in this endeavor, by allowing for a direct measure of the luminosity distance-redshift curve, and thereby producing detailed studies of the dark energy responsible for the Universe's currently accelerating expansion. As such observations are pushed to higher redshifts, the observed flux is increasingly affected by gravitational lensing magnification due to intervening structure along the line-of-sight. We simulate and analyze the non-Gaussian probability distribution function of de/amplification due to lensing of standard candles, quantify the effect of a convolution over many independent sources (which acts to restore the intrinsic average (unlensed) luminosity due to flux conservation), and compute the additional uncertainty due to lensing on derived cosmological parameters. For example, the ``degradation factor'' due to lensing is a factor of three reduction in the effective number of usable supernovae at z=1.5 (for sources with intrinsic flux dispersion of 10%). We also derive a useful expression for the effective increased dispersion in standard candles due to lensing, as a function of redshift.

preprint2001arXiv

Constant Crunch Coordinates for Black Hole Simulations

We reinvestigate the utility of time-independent constant mean curvature foliations for the numerical simulation of a single spherically-symmetric black hole. Each spacelike hypersurface of such a foliation is endowed with the same constant value of the trace of the extrinsic curvature tensor, $K$. Of the three families of $K$-constant surfaces possible (classified according to their asymptotic behaviors), we single out a sub-family of singularity-avoiding surfaces that may be particularly useful, and provide an analytic expression for the closest approach such surfaces make to the singularity. We then utilize a non-zero shift to yield families of $K$-constant surfaces which (1) avoid the black hole singularity, and thus the need to excise the singularity, (2) are asymptotically null, aiding in gravity wave extraction, (3) cover the physically relevant part of the spacetime, (4) are well behaved (regular) across the horizon, and (5) are static under evolution, and therefore have no ``grid stretching/sucking'' pathologies. Preliminary numerical runs demonstrate that we can stably evolve a single spherically-symmetric static black hole using this foliation. We wish to emphasize that this coordinatization produces $K$-constant surfaces for a single black hole spacetime that are regular, static and stable throughout their evolution.

preprint2001arXiv

Gravitational Wave Emission From Core-Collapse of Massive Stars

We derive estimates for the characteristics of gravitational radiation from stellar collapse, using recent models of the core-collapse of Chandrasekhar mass white dwarfs (accretion induced collapse), core-collapse supernovae and collapsars, and the collapse of very massive stars (~> 300 Msun). We study gravitational-wave emission mechanisms using several estimation techniques, including two-dimensional numerical computation of quadrupole wave emission, estimates of bar-mode strength, estimates of r-mode emission, and estimates of waves from black hole ringing. We also review the rate at which the relevant collapses are believed to occur, which has a major impact on their relevance as astrophysical sources. Although the latest supernova progenitor simulations produce cores rotating much slower than those used in the past, we find that bar-mode and r-mode instabilities from core-collapse supernovae remain among the leading candidate sources for LIGO-II. Accretion induced collapse (AIC) of a white dwarf could produce gravitational-wave signals similar to those from core-collapse. In the models that we examine, such collapses are not unstable to bar modes; we note that models recently examined by Liu and Lindblom, which have slightly more angular momentum, are certainly unstable to bar formation. Because AIC events are probably 1,000 times less common than core-collapse supernovae, the typical AIC event will be much further away, and thus the observed waves will be much weaker. In the most optimistic circumstances, we find it may be possible to detect gravitational waves from the collapse of 300 Msun Population III stars.

preprint1998arXiv

A New Method for Determining Cumulative Gravitational Lensing Effects in Inhomogeneous Universes

We present a new approach to calculating the statistical distributions for magnification, shear, and rotation of images of cosmological sources due to gravitational lensing. In this approach one specifies an underlying Robertson-Walker cosmological model together with relevant information on the clumping of matter on scales much smaller than the Hubble radius. The geodesic deviation equation is then integrated backwards in time until the desired redshift is reached, using a Monte Carlo procedure wherein each photon beam in effect ``creates its own universe'' as it propagates. The approach is somewhat similar to that used in ``Swiss cheese'' models, but the ``cheese'' has been completely eliminated, the matter distribution in the ``voids'' need not be spherically symmetric, and the total mass in each void need equal the corresponding Robertson-Walker mass only on average. We show that when the matter in the universe consists of point masses, the statistical distributions of the lensing images are essentially independent of both the mass spectrum and clustering properties of the point masses, provided that the clustering is spherical. We apply our results (i) to argue that the positive correlation recently found between quasar luminosity and the number of absorption line systems is not likely to be due to lensing, and (ii) to determine the amount of ``noise'' and possible bias produced by lensing in measurements of $q_0$ using distant supernovae.

preprint1995arXiv

Photon Statistics Limits for Earth-Based Parallax Measurements of MACHO Events

We analyze the limitations imposed by photon counting statistics on extracting useful information about MACHOs from Earth-based parallax observations of microlensing events. We find that if one or more large (say $2.5\,\rm m$) telescopes are dedicated to observing a MACHO event for several nights near maximum amplification, then it is possible, in principle, to measure the velocity of the MACHO well enough to distinguish between disk and halo populations for events with $ω{A_m}\!^2 \gta 1\,\rm day^{-1}$, where $ω^{-1}$ denotes the timescale of the event and $A_m$ denotes its maximum amplification. Thus, if it turns out to be possible to reduce all other sources of error to below that of photon statistics, it may be feasible to do useful Earth-based parallax observations for high amplification events.

preprint1994arXiv

Spin Dynamics of the LAGEOS Satellite in Support of a Measurement of the Earth's Gravitomagnetism

LAGEOS is an accurately-tracked, dense spherical satellite covered with 426 retroreflectors. The tracking accuracy is such as to yield a medium term (years to decades) inertial reference frame determined via relatively inexpensive observations. This frame is used as an adjunct to the more difficult and data intensive VLBI absolute frame measurements. There is a substantial secular precession of the satellite's line of nodes consistent with the classical, Newtonian precession due to the non-sphericity of the earth. Ciufolini has suggested the launch of an identical satellite (LAGEOS-3) into an orbit supplementary to that of LAGEOS-1: LAGEOS-3 would then experience an equal and opposite classical precession to that of LAGEOS-1. Besides providing a more accurate real-time measurement of the earth's length of day and polar wobble, this paired-satellite experiment would provide the first direct measurement of the general relativistic frame-dragging effect. Of the five dominant error sources in this experiment, the largest one involves surface forces on the satellite, and their consequent impact on the orbital nodal precession. The surface forces are a function of the spin dynamics of the satellite. Consequently, we undertake here a theoretical effort to model the spin ndynamics of LAGEOS. In this paper we present our preliminary results.

Daniel E. Holz

What is connected

Connect this record

See the researcher in context

Building this map preview

57 published item(s)

False positive control in time series coincidence detection

Avoiding a Cluster Catastrophe: Retention Efficiency and the Binary Black Hole Mass Spectrum

Bridging the Gap: Categorizing Gravitational-Wave Events at the Transition Between Neutron Stars and Black Holes

Cosmology Intertwined: A Review of the Particle Physics, Astrophysics, and Cosmology Associated with the Cosmological Tensions and Anomalies

Snowmass2021 Cosmic Frontier White Paper: Fundamental Physics and Beyond the Standard Model

Spectral sirens: cosmology from the full mass distribution of compact binaries

Target of Opportunity Observations of Gravitational Wave Events with Vera C. Rubin Observatory

Distance measures in gravitational-wave astrophysics and cosmology

Please repeat: Strong lensing of gravitational waves as a probe of compact binary and galaxy populations

Probing Extremal Gravitational-Wave Events with Coarse-Grained Likelihoods

The Gravity Collective: A Search for the Electromagnetic Counterpart to the Neutron Star-Black Hole Merger GW190814

Black Hole Coagulation: Modeling Hierarchical Mergers in Black Hole Populations

Cosmological Inference using Gravitational Wave Standard Sirens: A Mock Data Challenge

Does Matter Matter? Using the mass distribution to distinguish neutron stars and black holes

Picky Partners: The Pairing of Component Masses in Binary Black Hole Mergers

Counting on Short Gamma-Ray Bursts: Gravitational-Wave Constraints of Jet Geometry

Nonparametric Inference of Neutron Star Composition, Equation of State, and Maximum Mass with GW170817

The Most Massive Binary Black Hole Detections and the Identification of Population Outliers

Compact Binary Merger Rates: Comparison with LIGO/Virgo Upper Limits

Finding the One: Identifying the Host Galaxies of Gravitational-Wave Sources

The first gravitational-wave source from the isolated evolution of two 40-100 Msun stars

Neutron stars versus black holes: probing the mass gap with LIGO/Virgo

Comparison of LIGO/Virgo upper limits with predicted compact binary merger rates

Dark Sky Simulations: Early Data Release

The Formation and Gravitational-Wave Detection of Massive Stellar Black-Hole Binaries

The Loudest Gravitational Wave Events

Determining the Hubble constant from gravitational wave observations of merging compact binaries

Double Compact Objects II: Cosmological Merger Rates

Finding the First Cosmic Explosions I: Pair-Instability Supernovae

Finding the First Cosmic Explosions II: Core-Collapse Supernovae

Illuminating the Primeval Universe with Type IIn Supernovae

Supermassive Seeds for Supermassive Black Holes

Evidence for Type Ia Supernova Diversity from Ultraviolet Observations with the Hubble Space Telescope

GRB beaming and gravitational-wave observations

Seeing the First Supernovae at the Edge of the Universe with JWST

Localizing compact binary inspirals on the sky using ground-based gravitational wave interferometers

Outflowing Galactic Winds in Post-starburst and AGN Host Galaxies at 0.2<z<0.8

Exploring short gamma-ray bursts as gravitational-wave standard sirens

Reducing the weak lensing noise for the gravitational wave Hubble diagram using the non-Gaussianity of the magnification distribution

The most massive objects in the Universe

Weak lensing and dark energy: the impact of dark energy on nonlinear dark matter clustering

Astro2010 Decadal Survey Whitepaper: Coordinated Science in the Gravitational and Electromagnetic Skies

No Evidence for Dark Energy Dynamics from a Global Analysis of Cosmological Data

On The Origin Of The Highest Redshift Gamma-Ray Bursts

Ultra-high precision cosmology from gravitational waves

Implications of Two Type Ia Supernova Populations for Cosmological Measurements

Measuring dark energy spatial inhomogeneity with supernova data

Precision Determination of the Mass Function of Dark Matter Halos

Using gravitational-wave standard sirens

Consequences of gravitational radiation recoil

Gravitational Waves from Stellar Collapse: Correlations to Explosion Asymmetries

Safety in numbers: Gravitational Lensing Degradation of the Luminosity Distance-Redshift Relation

Constant Crunch Coordinates for Black Hole Simulations

Gravitational Wave Emission From Core-Collapse of Massive Stars

A New Method for Determining Cumulative Gravitational Lensing Effects in Inhomogeneous Universes

Photon Statistics Limits for Earth-Based Parallax Measurements of MACHO Events

Spin Dynamics of the LAGEOS Satellite in Support of a Measurement of the Earth's Gravitomagnetism