Showing 1–31 of 31 results for author: Chernoff, D F

Waveform Modelling for the Laser Interferometer Space Antenna

Authors: LISA Consortium Waveform Working Group, Niayesh Afshordi, Sarp Akçay, Pau Amaro Seoane, Andrea Antonelli, Josu C. Aurrekoetxea, Leor Barack, Enrico Barausse, Robert Benkel, Laura Bernard, Sebastiano Bernuzzi, Emanuele Berti, Matteo Bonetti, Béatrice Bonga, Gabriele Bozzola, Richard Brito, Alessandra Buonanno, Alejandro Cárdenas-Avendaño, Marc Casals, David F. Chernoff, Alvin J. K. Chua, Katy Clough, Marta Colleoni, Mekhi Dhesi, Adrien Druart , et al. (121 additional authors not shown)

Abstract: LISA, the Laser Interferometer Space Antenna, will usher in a new era in gravitational-wave astronomy. As the first anticipated space-based gravitational-wave detector, it will expand our view to the millihertz gravitational-wave sky, where a spectacular variety of interesting new sources abound: from millions of ultra-compact binaries in our Galaxy, to mergers of massive black holes at cosmologic… ▽ More LISA, the Laser Interferometer Space Antenna, will usher in a new era in gravitational-wave astronomy. As the first anticipated space-based gravitational-wave detector, it will expand our view to the millihertz gravitational-wave sky, where a spectacular variety of interesting new sources abound: from millions of ultra-compact binaries in our Galaxy, to mergers of massive black holes at cosmological distances; from the beginnings of inspirals that will venture into the ground-based detectors' view to the death spiral of compact objects into massive black holes, and many sources in between. Central to realising LISA's discovery potential are waveform models, the theoretical and phenomenological predictions of the pattern of gravitational waves that these sources emit. This white paper is presented on behalf of the Waveform Working Group for the LISA Consortium. It provides a review of the current state of waveform models for LISA sources, and describes the significant challenges that must yet be overcome. △ Less

Submitted 20 December, 2023; v1 submitted 2 November, 2023; originally announced November 2023.

Comments: 239 pages, 11 figures, white paper from the LISA Consortium Waveform Working Group, invited for submission to Living Reviews in Relativity, updated with comments from community

Modeling the Beam of Gravitational Radiation from a Cosmic String Loop

Authors: Namitha Suresh, David F. Chernoff

Abstract: We investigate exact and approximate techniques to calculate the emission of gravitational radiation from cosmic string loops in order to generate beam models covering the entire celestial sphere for a wide range of modes $m$. One approach entails summing over contributions of stationary and nearly-stationary points of individual, factorized, left- and right-moving modes. This "multipoint method"… ▽ More We investigate exact and approximate techniques to calculate the emission of gravitational radiation from cosmic string loops in order to generate beam models covering the entire celestial sphere for a wide range of modes $m$. One approach entails summing over contributions of stationary and nearly-stationary points of individual, factorized, left- and right-moving modes. This "multipoint method" generalizes traditional methods that rely on expansions around exact stationary points of mode products. A second complementary approach extends the method of steepest descent to generate an asymptotic description of the beam as $m \to \infty$. We present example calculations of the emission of power from cusp-containing loops and compare the results with those obtained by numerically exact techniques as well as by previous approaches. The multipoint method achieves its best results at an intermediate range of modes, improving over previous methods in terms of accuracy. It handles the emission from all regions of the loop not just those near cusps. We demonstrate this capability by making a detailed study of the "pseudocusp" phenomenon. △ Less

Submitted 26 June, 2024; v1 submitted 1 October, 2023; originally announced October 2023.

Comments: 34 pages, 36 figures. Revised discussion section

Journal ref: Phys. Rev. D 109, 123540 (2024)

Periodic Fast Radio Bursts from Axion Emission by Cosmic Superstrings

Authors: David F. Chernoff, Shing Yan Li, S. -H. Henry Tye

Abstract: We propose that the periodic fast radio bursts of FRB 180916.J0158+65 are sourced by axion emission (mass $m_{a} \sim 10^{-14}$ eV) from cosmic superstrings. Some of the emitted axions are converted to photons by magnetic fields as they travel along the line of sight to Earth. An impulsive burst of axion emission generates a photon signal typically lasting for milliseconds and varying with frequen… ▽ More We propose that the periodic fast radio bursts of FRB 180916.J0158+65 are sourced by axion emission (mass $m_{a} \sim 10^{-14}$ eV) from cosmic superstrings. Some of the emitted axions are converted to photons by magnetic fields as they travel along the line of sight to Earth. An impulsive burst of axion emission generates a photon signal typically lasting for milliseconds and varying with frequency in the observed manner. We find a range of parameters in our cosmic string network model consistent with the properties of FRB 180916.J0158+65. We suggest followup gravitational wave observations to test our model. △ Less

Submitted 16 March, 2020; originally announced March 2020.

Comments: 15 pages, 1 figure

Prospects of Cosmic Superstring Detection through Microlensing of Extragalactic Point-Like Sources

Authors: David F. Chernoff, Ariel Goobar, Janina J. Renk

Abstract: The existence of cosmic superstrings may be probed by astronomical time domain surveys. When crossing the line of sight to point-like sources, strings produce a distinctive microlensing signature. We consider two avenues to hunt for a relic population of superstring loops: frequent monitoring of (1) stars in Andromeda, lensed by loops in the haloes of the Milky-Way and Andromeda and (2) supernovae… ▽ More The existence of cosmic superstrings may be probed by astronomical time domain surveys. When crossing the line of sight to point-like sources, strings produce a distinctive microlensing signature. We consider two avenues to hunt for a relic population of superstring loops: frequent monitoring of (1) stars in Andromeda, lensed by loops in the haloes of the Milky-Way and Andromeda and (2) supernovae at cosmological distances, lensed by loops in the intergalactic medium. We assess the potential of such experiments to detect and/or constrain strings with a range of tensions, $10^{-15} \lesssim G μ/c^2 \lesssim 10^{-6}$. The practical sensitivity is tied to cadence of observations which we explore in detail. We forecast that high-cadence monitoring of $\sim 10^5$ stars on the far side of Andromeda over a year-long period will detect microlensing events if $Gμ/c^2 \sim 10^{-13}$, while $\sim 10^6$ stars will detect events if $10^{-13.5} <Gμ/c^2 < 10^{-11.5}$; the upper and lower bounds of the accessible tension range continue to expand as the number of stars rises. We also analyse the ability to reject models in the absence of fluctuations. While challenging, these studies are within reach of forthcoming time-domain surveys. Supernova observations can hypothetically constrain models with $10^{-12} < Gμ/c^2 < 10^{-6}$ without any optimisation of the survey cadence. However, the event rate forecast suggests it will be difficult to reject models of interest. As a demonstration, we use observations from the Pantheon Type Ia supernova cosmology data-set to place modest constraints on the number density of cosmic superstrings in a poorly tested region of the parameter space. △ Less

Submitted 11 October, 2019; v1 submitted 9 May, 2019; originally announced May 2019.

Comments: 17 pages, 7 figures + appendix. Main results summarised in Fig. 5; notation clarification & minor revisions w.r.t. v1 -- accepted for publication by MNRAS

Models and Simulations for the Photometric LSST Astronomical Time Series Classification Challenge (PLAsTiCC)

Authors: R. Kessler, G. Narayan, A. Avelino, E. Bachelet, R. Biswas, P. J. Brown, D. F. Chernoff, A. J. Connolly, M. Dai, S. Daniel, R. Di Stefano, M. R. Drout, L. Galbany, S. González-Gaitán, M. L. Graham, R. Hložek, E. E. O. Ishida, J. Guillochon, S. W. Jha, D. O. Jones, K. S. Mandel, D. Muthukrishna, A. O'Grady, C. M. Peters, J. R. Pierel , et al. (4 additional authors not shown)

Abstract: We describe the simulated data sample for the "Photometric LSST Astronomical Time Series Classification Challenge" (PLAsTiCC), a publicly available challenge to classify transient and variable events that will be observed by the Large Synoptic Survey Telescope (LSST), a new facility expected to start in the early 2020s. The challenge was hosted by Kaggle, ran from 2018 September 28 to 2018 Decembe… ▽ More We describe the simulated data sample for the "Photometric LSST Astronomical Time Series Classification Challenge" (PLAsTiCC), a publicly available challenge to classify transient and variable events that will be observed by the Large Synoptic Survey Telescope (LSST), a new facility expected to start in the early 2020s. The challenge was hosted by Kaggle, ran from 2018 September 28 to 2018 December 17, and included 1,094 teams competing for prizes. Here we provide details of the 18 transient and variable source models, which were not revealed until after the challenge, and release the model libraries at https://doi.org/10.5281/zenodo.2612896. We describe the LSST Operations Simulator used to predict realistic observing conditions, and we describe the publicly available SNANA simulation code used to transform the models into observed fluxes and uncertainties in the LSST passbands (ugrizy). Although PLAsTiCC has finished, the publicly available models and simulation tools are being used within the astronomy community to further improve classification, and to study contamination in photometrically identified samples of type Ia supernova used to measure properties of dark energy. Our simulation framework will continue serving as a platform to improve the PLAsTiCC models, and to develop new models. △ Less

Submitted 10 July, 2019; v1 submitted 27 March, 2019; originally announced March 2019.

Detection of Low Tension Cosmic Superstrings

Authors: David F. Chernoff, S. -H. Henry Tye

Abstract: Cosmic superstrings of string theory differ from conventional cosmic strings of field theory. We review how the physical and cosmological properties of the macroscopic string loops influence experimental searches for these relics from the epoch of inflation. The universe's average density of cosmic superstrings can easily exceed that of conventional cosmic strings having the same tension by two or… ▽ More Cosmic superstrings of string theory differ from conventional cosmic strings of field theory. We review how the physical and cosmological properties of the macroscopic string loops influence experimental searches for these relics from the epoch of inflation. The universe's average density of cosmic superstrings can easily exceed that of conventional cosmic strings having the same tension by two or more orders of magnitude. The cosmological behavior of the remnant superstring loops is qualitatively distinct because the string tension is exponentially smaller than the string scale in flux compactifications in string theory. Low tension superstring loops live longer, experience less recoil (rocket effect from the emission of gravitational radiation) and tend to cluster like dark matter in galaxies. Clustering enhances the string loop density with respect to the cosmological average in collapsed structures in the universe. The enhancement at the Sun's position is $\sim 10^5$. We develop a model encapsulating the leading order string theory effects, the current understanding of the string network loop production and the influence of cosmological structure formation suitable for forecasting the detection of superstring loops via optical microlensing, gravitational wave bursts and fast radio bursts. We evaluate the detection rate of bursts from cusps and kinks by LIGO- and LISA-like experiments. Clustering dominates rates for $G μ< 10^{-11.9}$ (LIGO cusp), $G μ<10^{-11.2}$ (LISA cusp), $G μ< 10^{-10.6}$ (LISA kink); we forecast experimentally accessible gravitational wave bursts for $G μ>10^{-14.2}$ (LIGO cusp), $G μ>10^{-15}$ (LISA cusp) and $G μ>10^{- 14.1}$ (LISA kink). △ Less

Submitted 5 February, 2018; v1 submitted 13 December, 2017; originally announced December 2017.

Comments: 81 pages, 15 figures; refs. added

Pulsar Timing Perturbations from Galactic Gravitational Wave Bursts with Memory

Authors: Dustin R. Madison, David F. Chernoff, James M. Cordes

Abstract: Pulsar timing arrays (PTAs) are used to search for long-wavelength gravitational waves (GWs) by monitoring a set of spin-stable millisecond pulsars. Most theoretical analyses assume that the relevant GW sources are much more distant from Earth than the pulsars comprising the array. Unlike ground- or solar system-based GW detectors, PTAs might well contain embedded GW sources. We derive the PTA res… ▽ More Pulsar timing arrays (PTAs) are used to search for long-wavelength gravitational waves (GWs) by monitoring a set of spin-stable millisecond pulsars. Most theoretical analyses assume that the relevant GW sources are much more distant from Earth than the pulsars comprising the array. Unlike ground- or solar system-based GW detectors, PTAs might well contain embedded GW sources. We derive the PTA response from sources at any distance, with a specific focus on GW bursts with memory (BWMs). We consider supernovae and compact binary mergers as potential Galactic BWM sources and evaluate the signature for an array with pulsars in globular clusters or in the Galactic center. Understanding the response of PTAs to nearby sources of BWM is a step towards investigating other more complex Galactic sources. △ Less

Submitted 13 October, 2017; originally announced October 2017.

Comments: 14 pages, 7 figures, submitted to PRD

Journal ref: Phys. Rev. D 96, 123016 (2017)

Inflation, String Theory and Cosmology

Authors: David F. Chernoff, S. -H. Henry Tye

Abstract: At its very beginning, the universe is believed to have grown exponentially in size via the mechanism of inflation. The almost scale-invariant density perturbation spectrum predicted by inflation is strongly supported by cosmological observations, in particular the cosmic microwave background radiation. However, the universe's precise inflationary scenario remains a profound problem for cosmology… ▽ More At its very beginning, the universe is believed to have grown exponentially in size via the mechanism of inflation. The almost scale-invariant density perturbation spectrum predicted by inflation is strongly supported by cosmological observations, in particular the cosmic microwave background radiation. However, the universe's precise inflationary scenario remains a profound problem for cosmology and for fundamental physics. String theory, the most-studied theory as the final physical theory of nature, should provide an answer to this question. Some of the proposals on how inflation is realized in string theory are reviewed. Since everything is made of strings, some string loops of cosmological sizes are likely to survive in the hot big bang that followed inflation. They appear as cosmic strings, which can have intricate properties. Because of the warped geometry in flux compactification of the extra spatial dimensions in string theory, some of the cosmic strings may have tensions substantially below the Planck or string scale. Such strings cluster in a manner similar to dark matter leading to hugely enhanced densities. As a result, numerous fossil remnants of the low tension cosmic strings may exist within the galaxy. They can be revealed through the optical lensing of background stars in the near future and studied in detail through gravitational wave emission. We anticipate that these cosmic strings will permit us to address central questions about the properties of string theory as well as the birth of our universe. △ Less

Submitted 20 December, 2014; v1 submitted 1 December, 2014; originally announced December 2014.

Comments: 61 pages, 10 figures Added references and summary section; expanded introduction

Dynamically formed black hole+millisecond pulsar binaries in globular clusters

Authors: Drew Clausen, Steinn Sigurdsson, David F. Chernoff

Abstract: The discovery of a binary comprising a black hole (BH) and a millisecond pulsar (MSP) would yield insights into stellar evolution and facilitate exquisitely sensitive tests of general relativity. Globular clusters (GCs) are known to harbor large MSP populations and recent studies suggest that GCs may also retain a substantial population of stellar mass BHs. We modeled the formation of BH+MSP binar… ▽ More The discovery of a binary comprising a black hole (BH) and a millisecond pulsar (MSP) would yield insights into stellar evolution and facilitate exquisitely sensitive tests of general relativity. Globular clusters (GCs) are known to harbor large MSP populations and recent studies suggest that GCs may also retain a substantial population of stellar mass BHs. We modeled the formation of BH+MSP binaries in GCs through exchange interactions between binary and single stars. We found that in dense, massive clusters most of the dynamically formed BH+MSP binaries will have orbital periods of 2 to 10 days, regardless of the mass of the BH, the number of BHs retained by the cluster, and the nature of the GC's binary population. The size of the BH+MSP population is sensitive to several uncertain parameters, including the BH mass function, the BH retention fraction, and the binary fraction in GCs. Based on our models, we estimate that there are $0.6\pm0.2$ dynamically formed BH+MSP binaries in the Milky Way GC system, and place an upper limit on the size of this population of $\sim 10$. Interestingly, we find that BH+MSP binaries will be rare even if GCs retain large BH populations. △ Less

Submitted 13 June, 2014; v1 submitted 29 April, 2014; originally announced April 2014.

Comments: 15 pages, 5 figures, 2 tables, accepted for publication in MNRAS, updated to match published version

Cosmic String Loop Microlensing

Authors: Jolyon K. Bloomfield, David F. Chernoff

Abstract: Cosmic superstring loops within the galaxy microlens background point sources lying close to the observer-string line of sight. For suitable alignments, multiple paths coexist and the (achromatic) flux enhancement is a factor of two. We explore this unique type of lensing by numerically solving for geodesics that extend from source to observer as they pass near an oscillating string. We characteri… ▽ More Cosmic superstring loops within the galaxy microlens background point sources lying close to the observer-string line of sight. For suitable alignments, multiple paths coexist and the (achromatic) flux enhancement is a factor of two. We explore this unique type of lensing by numerically solving for geodesics that extend from source to observer as they pass near an oscillating string. We characterize the duration of the flux doubling and the scale of the image splitting. We probe and confirm the existence of a variety of fundamental effects predicted from previous analyses of the static infinite straight string: the deficit angle, the Kaiser-Stebbins effect, and the scale of the impact parameter required to produce microlensing. Our quantitative results for dynamical loops vary by O(1) factors with respect to estimates based on infinite straight strings for a given impact parameter. A number of new features are identified in the computed microlensing solutions. Our results suggest that optical microlensing can offer a new and potentially powerful methodology for searches for superstring loop relics of the inflationary era. △ Less

Submitted 27 November, 2013; originally announced November 2013.

Comments: 20 pages, 19 figures

Journal ref: Phys. Rev. D 89, 124003 (2014)

Modelling non-linear evolution using Lagrangian Perturbation Theory (LPT) re-expansions

Authors: Sharvari Nadkarni-Ghosh, David F. Chernoff

Abstract: We present a new method to calculate formation of cosmological structure in the Newtonian limit. The method is based on Lagrangian perturbation theory plus two key theoretical extensions. One advance involves identifying and fixing a previously ignored gauge-like degree of freedom relating quantities calculated in LPT to those measured by a preferred Friedmann-Robertson-Walker (FRW) observer. Hand… ▽ More We present a new method to calculate formation of cosmological structure in the Newtonian limit. The method is based on Lagrangian perturbation theory plus two key theoretical extensions. One advance involves identifying and fixing a previously ignored gauge-like degree of freedom relating quantities calculated in LPT to those measured by a preferred Friedmann-Robertson-Walker (FRW) observer. Handling this connection between calculational and observer frames is physically essential and ensures a momentum conserving description. The second extension is to systematically re-expand the equations of motion to increase LPT's radius of convergence to the maximum future time prior to orbit crossing. The paper implements a complete algorithm and performs extensive "proof of principle" tests of the new method, including direct comparison to known solutions, evaluation of conserved quantities and formal convergence studies. All are satisfactory. We show convergence is exponential in grid size and Lagrangian order and polynomial in step size. There are three {\it powerful advantages} afforded by the new technique: (1) it employs a smooth representation of all fields and the results are not limited by particle induced shot-noise errors, (2) it permits the numerical error to be controlled by changing Lagrangian order and/or number of steps allowing, in principle, arbitrarily small errors to be achieved prior to orbit crossing and (3) it handles generic cold initial data (any periodic density and velocity fields, including those with initial rotational components). Together, these properties make the new technique well-suited to handle quasi-linear scales where analytic methods and/or numerical simulations fail to provide suitably accurate answers. △ Less

Submitted 28 February, 2013; v1 submitted 25 November, 2012; originally announced November 2012.

Comments: 31 pages, 8 figrues, 2 tables. Accepted in MNRAS

Journal ref: MNRAS, Volume 431, Issue 1, p.799-823, May 2013

Black Hole-Neutron Star Mergers in Globular Clusters

Authors: Drew Clausen, Steinn Sigurdsson, David F. Chernoff

Abstract: We model the formation of black hole-neutron star (BH-NS) binaries via dynamical interactions in globular clusters. We find that in dense, massive clusters, 16-61% of the BH-NS binaries formed by interactions with existing BH binaries will undergo mergers driven by the emission of gravitational radiation. If the BHs are retained by the cluster after merging with a NS, the BHs acquire subsequent NS… ▽ More We model the formation of black hole-neutron star (BH-NS) binaries via dynamical interactions in globular clusters. We find that in dense, massive clusters, 16-61% of the BH-NS binaries formed by interactions with existing BH binaries will undergo mergers driven by the emission of gravitational radiation. If the BHs are retained by the cluster after merging with a NS, the BHs acquire subsequent NS companions and undergo several mergers. Thus, the merger rate depends critically upon whether or not the BH is retained by the cluster after the merger. Results from numerical relativity suggest that kick imparted to a ~7 M_sun BH after it merges with a NS will greatly exceed the cluster's escape velocity. In this case, the models suggest that the majority of BH-NS mergers in globular clusters occur within 4 Gyrs of the cluster's formation and would be unobservable by Advanced LIGO. For more massive BHs, on the other hand, the post merger kick is suppressed and the BH is retained. Models with 35 M_sun BHs predict Advanced LIGO detection rates in the range 0.04 - 0.7 per year. On the pessimistic end of this range, BH-NS mergers resulting from binary-single star interactions in globular clusters could account for an interesting fraction of all BH-NS mergers. On the optimistic end, this channel may dominate the rate of detectable BH-NS mergers. △ Less

Submitted 30 October, 2012; originally announced October 2012.

Comments: 13 pages, 7 figures, 3 tabels, accepted for publication in MNRAS

Bayesian Methods for Analysis and Adaptive Scheduling of Exoplanet Observations

Authors: Thomas J. Loredo, James O. Berger, David F. Chernoff, Merlise A. Clyde, Bin Liu

Abstract: We describe work in progress by a collaboration of astronomers and statisticians developing a suite of Bayesian data analysis tools for extrasolar planet (exoplanet) detection, planetary orbit estimation, and adaptive scheduling of observations. Our work addresses analysis of stellar reflex motion data, where a planet is detected by observing the "wobble" of its host star as it responds to the gra… ▽ More We describe work in progress by a collaboration of astronomers and statisticians developing a suite of Bayesian data analysis tools for extrasolar planet (exoplanet) detection, planetary orbit estimation, and adaptive scheduling of observations. Our work addresses analysis of stellar reflex motion data, where a planet is detected by observing the "wobble" of its host star as it responds to the gravitational tug of the orbiting planet. Newtonian mechanics specifies an analytical model for the resulting time series, but it is strongly nonlinear, yielding complex, multimodal likelihood functions; it is even more complex when multiple planets are present. The parameter spaces range in size from few-dimensional to dozens of dimensions, depending on the number of planets in the system, and the type of motion measured (line-of-sight velocity, or position on the sky). Since orbits are periodic, Bayesian generalizations of periodogram methods facilitate the analysis. This relies on the model being linearly separable, enabling partial analytical marginalization, reducing the dimension of the parameter space. Subsequent analysis uses adaptive Markov chain Monte Carlo methods and adaptive importance sampling to perform the integrals required for both inference (planet detection and orbit measurement), and information-maximizing sequential design (for adaptive scheduling of observations). We present an overview of our current techniques and highlight directions being explored by ongoing research. △ Less

Submitted 10 May, 2018; v1 submitted 29 July, 2011; originally announced August 2011.

Comments: 29 pages, 11 figures. An abridged version is accepted for publication in Statistical Methodology for a special issue on astrostatistics, with selected (refereed) papers presented at the Astronomical Data Analysis Conference (ADA VI) held in Monastir, Tunisia, in May 2010. Update corrects equation (3)

Journal ref: Statistical Methodology 9 (2012) 101-114

Extending the domain of validity of the Lagrangian approximation

Authors: Sharvari Nadkarni-Ghosh, David F. Chernoff

Abstract: We investigate convergence of Lagrangian Perturbation Theory (LPT) by analyzing the model problem of a spherical homogeneous top-hat in an Einstein-deSitter background cosmology. We derive the formal structure of the LPT series expansion, working to arbitrary order in the initial perturbation amplitude. The factors that regulate LPT convergence are identified by studying the exact, analytic soluti… ▽ More We investigate convergence of Lagrangian Perturbation Theory (LPT) by analyzing the model problem of a spherical homogeneous top-hat in an Einstein-deSitter background cosmology. We derive the formal structure of the LPT series expansion, working to arbitrary order in the initial perturbation amplitude. The factors that regulate LPT convergence are identified by studying the exact, analytic solution expanded according to this formal structure. The key methodology is to complexify the exact solution, demonstrate that it is analytic and apply well-known convergence criteria for power series expansions of analytic functions.This analysis fully explains the previously reported observation that LPT fails to predict the evolution of an underdense, open region beyond a certain time. It also implies the existence of other examples, including overdense, closed regions, for which LPT predictions should also fail. We show that this is indeed the case by numerically computing the LPT expansion in these problematic cases. The formal limitations to the validity of LPT expansion are considerably more complicated than simply the first occurrence of orbit crossings as is often assumed. Evolution to a future time generically requires re-expanding the solution in overlapping domains that ultimately link the initial and final times, each domain subject to its own convergence criterion. We demonstrate that it is possible to handle all the problematic cases by taking multiple steps (LPT re-expansion). We characterize how the leading order numerical error for a solution generated by LPT re- expansion varies with the choice of Lagrangian order and of time step size. Convergence occurs when the Lagrangian order increases and/or the time step size decreases in a simple, well-defined manner. We develop a recipe for time step control for LPT re-expansion based on these results. △ Less

Submitted 14 September, 2010; v1 submitted 7 May, 2010; originally announced May 2010.

Comments: 40 pages, 24 figures. Abstract abbreviated. Accepted for publication in MNRAS. Updated version with minor comments from referee incorporated

Journal ref: MNRAS, Jan 2011, Volume 410, Issue 3, pp. 1454-1488

Pseudospectral Calculation of the Wavefunction of Helium and the Negative Hydrogen Ion

Authors: Paul E. Grabowski, David F. Chernoff

Abstract: We study the numerical solution of the non-relativistic Schrödinger equation for two-electron atoms in ground and excited S-states using pseudospectral (PS) methods of calculation. The calculation achieves convergence rates for the energy, Cauchy error in the wavefunction, and variance in local energy that are exponentially fast for all practical purposes. The method requires three separate subd… ▽ More We study the numerical solution of the non-relativistic Schrödinger equation for two-electron atoms in ground and excited S-states using pseudospectral (PS) methods of calculation. The calculation achieves convergence rates for the energy, Cauchy error in the wavefunction, and variance in local energy that are exponentially fast for all practical purposes. The method requires three separate subdomains to handle the wavefunction's cusp-like behavior near the two-particle coalescences. The use of three subdomains is essential to maintaining exponential convergence. A comparison of several different treatments of the cusps and the semi-infinite domain suggest that the simplest prescription is sufficient. For many purposes it proves unnecessary to handle the logarithmic behavior near the three-particle coalescence in a special way. The PS method has many virtues: no explicit assumptions need be made about the asymptotic behavior of the wavefunction near cusps or at large distances, the local energy is exactly equal to the calculated global energy at all collocation points, local errors go down everywhere with increasing resolution, the effective basis using Chebyshev polynomials is complete and simple, and the method is easily extensible to other bound states. This study serves as a proof-of-principle of the method for more general two- and possibly three-electron applications. △ Less

Submitted 26 March, 2010; v1 submitted 9 September, 2009; originally announced September 2009.

Comments: 23 pages, 20 figures, 2 tables, Final refereed version - Some references added, some stylistic changes, added paragraph to matrix methods section, added last sentence to abstract.

Journal ref: Phys.Rev.A81:032508,2010

Clustering of Superstring Loops

Authors: David F. Chernoff

Abstract: The fragmentation of horizon crossing superstrings creates loops that accrete and bind to collapsing matter perturbations if the string tension is less than a critical value and if the loops are large-scale. For typical model parameters (velocity dispersion of newly formed loops, loop length distribution, etc.) the galactic loop density is greatly enhanced as compared to the universe's mean loop… ▽ More The fragmentation of horizon crossing superstrings creates loops that accrete and bind to collapsing matter perturbations if the string tension is less than a critical value and if the loops are large-scale. For typical model parameters (velocity dispersion of newly formed loops, loop length distribution, etc.) the galactic loop density is greatly enhanced as compared to the universe's mean loop density. The limiting enhancement of loop (energy) density is 0.25-0.4 times the enhancement of cold dark matter (ratio of galaxy to universe's dark matter densities) and is fully achieved for dimensionless string tensions less than 10^-13 at galactic scales less than 100 kpc. The enhancement is substantial for any tension less than the critical value which is of order 10^-9. The fact that the string loop enhancement roughly tracks that of cold dark matter is a robust result for small tension and large-scale fragmentation. Experiments sensitive to the local Galactic population of loops, especially microlensing, should enjoy characteristic detection rates greatly boosted over homogeneous estimates. △ Less

Submitted 27 August, 2009; originally announced August 2009.

Cosmic String Detection via Microlensing of Stars

Authors: David F. Chernoff, S. -H. Henry Tye

Abstract: Cosmic superstrings are produced towards the end of the brane inflation. If the string tension is low enough, loops tend to be relatively long-lived. The resultant string network is expected to contain many loops which are smaller than typical Galactic scales. Cosmic expansion damps the center of mass motion of the loops which then cluster like cold dark matter. Loops will lens stars within the… ▽ More Cosmic superstrings are produced towards the end of the brane inflation. If the string tension is low enough, loops tend to be relatively long-lived. The resultant string network is expected to contain many loops which are smaller than typical Galactic scales. Cosmic expansion damps the center of mass motion of the loops which then cluster like cold dark matter. Loops will lens stars within the galaxy and local group. We explore microlensing of stars as a tool to detect and to characterize some of the fundamental string and string network properties, including the dimensionless string tension $G μ/c^2$ and the density of string loops within the Galaxy. As $G μ\to 0$ the intrinsic microlensing rate diverges as $1/\sqrt{G μ}$ but experimental detection will be limited by shortness of the lensing timescale and/or smallness of the bending angle which each vary $\propto G μ$. We find that detection is feasible for a range of tensions. As an illustration, the planned optical astrometric survey mission, Gaia, should be able to detect numerous microlensing events for string networks with tensions $10^{-10} \simless G μ\simless 10^{-6}$. A null result for optical microlensing implies $G μ\simless 10^{-10}$. If lensing of a given source is observed it will repeat because the internal motions of the loop are relativistic but the center of mass motion may be much smaller, of order the halo velocity. This distinctive hallmark $\sim 1000$ repetitions, suggests a useful method for confirmation of a potential lensing detection. △ Less

Submitted 9 September, 2007; v1 submitted 7 September, 2007; originally announced September 2007.

Comments: 31 pages, 11 figures

On the Origin of the most massive stars around R136

Authors: B. R. Brandl, S. Portegies Zwart, A. J. F. Moffat, D. F. Chernoff

Abstract: We discuss the signature of a peculiar constellation of very massive stars at a projected distance of 2-3 pc around R136a. We discuss various scenarios for its possible origin, such as independent clusters, triggered star formation, supernovae, and ejections via dynamical interactions. If the latter scenario were the correct one this would have significant implications on the way to probe the co… ▽ More We discuss the signature of a peculiar constellation of very massive stars at a projected distance of 2-3 pc around R136a. We discuss various scenarios for its possible origin, such as independent clusters, triggered star formation, supernovae, and ejections via dynamical interactions. If the latter scenario were the correct one this would have significant implications on the way to probe the conditions in dense stellar cores, and on the evolution of massive clusters in general. △ Less

Submitted 26 January, 2005; originally announced January 2005.

Comments: To appear in the proceedings of "Massive Stars in Interacting Binaries", eds. A.J.F. Moffat & N. St-Louis

The Velocity Distribution of Isolated Radio Pulsars

Authors: Z. Arzoumanian, D. F. Chernoff, J. M. Cordes

Abstract: (Abridged) We infer the velocity distribution of radio pulsars by modelling their birth, evolution, and detection in large-scale 0.4 GHz pulsar surveys, and by comparing model distributions of measurable pulsar properties with survey data using a likelihood function. We test models that characterize a population's birth rate, luminosity, shutoff of radio emission, birth locations, and birth velo… ▽ More (Abridged) We infer the velocity distribution of radio pulsars by modelling their birth, evolution, and detection in large-scale 0.4 GHz pulsar surveys, and by comparing model distributions of measurable pulsar properties with survey data using a likelihood function. We test models that characterize a population's birth rate, luminosity, shutoff of radio emission, birth locations, and birth velocities. We infer that the radio beam luminosity (i) is comparable to the energy flux of relativistic particles in models for spin-driven magnetospheres, signifying that radio emission losses reach nearly 100% for the oldest pulsars; and (ii) scales approximately as sqrt(Edot) which, in magnetosphere models, is proportional to the voltage drop available for acceleration of particles. We find that a two-component velocity distribution with characteristic velocities of 90 km/s and 500 km/s is greatly preferred to any one-component distribution. We explore some consequences of the preferred birth velocity distribution: (i) roughly 50% of pulsars in the solar neighborhood will escape the Galaxy, while ~15% have velocities greater than 1000 km/s; (ii) an important low-velocity population exists that increases the fraction of neutron stars retained by globular clusters and is consistent with the number of old objects that accrete from the interstellar medium; (iii) under standard assumptions for supernova remnant expansion and pulsar spin-down, ~10% of pulsars younger than 20 kyr will appear to lie outside of their host remnants. Finally, we comment on the ramifications of our birth velocity distribution for binary survival and the population of inspiraling binary neutron stars relevant to some GRB models and potential sources for LIGO. △ Less

Submitted 29 November, 2001; v1 submitted 8 June, 2001; originally announced June 2001.

Comments: 16 pages, 4 figures, ApJ (in press). Revisions to reflect refereed version include a new sub(sub)section on effects of the 'standard candle" assumption, and a revised Figure 2. Essential results unchanged

Journal ref: Astrophys.J. 568 (2002) 289-301

Pulsar Jets: Implications for Neutron Star Kicks and Initial Spins

Authors: Dong Lai, David F. Chernoff, James M. Cordes

Abstract: We study implications for the apparent alignment of the spin axes, proper-motions, and polarization vectors of the Crab and Vela pulsars. The spin axes are deduced from recent Chandra X-ray Observatory images that reveal jets and nebular structure having definite symmetry axes. The alignments indicate these pulsars were born either in isolation or with negligible velocity contributions from bina… ▽ More We study implications for the apparent alignment of the spin axes, proper-motions, and polarization vectors of the Crab and Vela pulsars. The spin axes are deduced from recent Chandra X-ray Observatory images that reveal jets and nebular structure having definite symmetry axes. The alignments indicate these pulsars were born either in isolation or with negligible velocity contributions from binary motions. We examine the effects of rotation and the conditions under which spin-kick alignment is produced for various models of neutron star kicks. If the kick is generated when the neutron star first forms by asymmetric mass ejection or/and neutrino emission, then the alignment requires that the protoneutron star possesses an original spin with period $P_s$ much less than the kick timescale, thus spin-averaging the kick forces. The kick timescale ranges from 100 ms to 10 s depending on whether the kick is hydrodynamically driven or neutrino-magnetic field driven. For hydrodynamical models, spin-kick alignment further requires the rotation period of an asymmetry pattern at the radius near shock breakout (>100 km) to be much less than ~100 ms; this is difficult to satisfy unless rotation plays a dynamically important role in the core collapse and explosion ($P_s\lo 1$ ms). Aligned kick and spin vectors are inherent to the slow process of asymmetric electromagnetic radiation from an off-centered magnetic dipole. We reassess the viability of this effect, correcting a factor of 4 error in Harrison and Tademaru's calculation that increases the size of the effect. To produce a kick velocity of order a few hundred km/s requires that the neutron star be born with an initial spin close to 1 ms and that spindown due to r-mode driven gravitational radiation be inefficient compared to standard magnetic braking. △ Less

Submitted 1 December, 2000; v1 submitted 18 July, 2000; originally announced July 2000.

Comments: Small changes/additions; final version to be published in ApJ, Vol.549 (March 10, 2001)

Effects of magnetic fields on radiatively overstable shock waves

Authors: Paul A. Kimoto, David F. Chernoff

Abstract: We discuss high-resolution simulations of one-dimensional, plane-parallel shock waves with mean speeds between 150 and 240 km/s propagating into gas with Alfven velocities up to 40 km/s and outline the conditions under which these radiative shocks experience an oscillatory instability in the cooling length, shock velocity, and position of the shock front. We investigate two forms of postshock co… ▽ More We discuss high-resolution simulations of one-dimensional, plane-parallel shock waves with mean speeds between 150 and 240 km/s propagating into gas with Alfven velocities up to 40 km/s and outline the conditions under which these radiative shocks experience an oscillatory instability in the cooling length, shock velocity, and position of the shock front. We investigate two forms of postshock cooling: a truncated single power law and a more realistic piecewise power law. The degree of nonlinearity of the instability depends strongly on the cooling power law and the Alfven Mach number: for power-law indices α< 0 typical magnetic field strengths may be insufficient either to stabilize the fundamental oscillatory mode or to prevent the oscillations from reaching nonlinear amplitudes. △ Less

Submitted 31 July, 1997; originally announced July 1997.

Comments: 11 text pages, LaTeX/AASTeX (aaspp4); 5 figures; accepted by ApJ

Neutron Star Population Dynamics.II: 3D Space Velocities of Young Pulsars

Authors: J. M. Cordes, David F. Chernoff

Abstract: We use astrometric, distance and spindown data on pulsars to: (1) estimate three-dimensional velocity components, birth distances from the galactic plane, and ages of individual objects; (2) determine the distribution of space velocities and the scale height of pulsar progenitors; (3) test spindown laws for pulsars; (4) test for correlations between space velocities and other pulsar parameters;… ▽ More We use astrometric, distance and spindown data on pulsars to: (1) estimate three-dimensional velocity components, birth distances from the galactic plane, and ages of individual objects; (2) determine the distribution of space velocities and the scale height of pulsar progenitors; (3) test spindown laws for pulsars; (4) test for correlations between space velocities and other pulsar parameters; and (5) place empirical requirements on mechanisms than can produce high velocity neutron stars. Our approach incorporates measurement errors, uncertainties in distances, deceleration in the Galactic potential, and differential galactic rotation. We find that the scale height of the progenitors is approximately 0.13 kpc, that the 3D velocities are distributed in two components with characteristic speeds of 175(+20,-30) km/s and 700(+200,-150) km/s representing 83% and 17% of the population respectively. These results are insensitive to the explicit relation of chronological and spindown ages. We infer that the most probable chronological ages are typically smaller than conventional spindown ages by factors as large as two. We assess mechanisms for producing high-velocity neutron stars in view of the derived velocity distribution function. △ Less

Submitted 28 July, 1997; originally announced July 1997.

Comments: 54 pages, 13 Postscript figures, uses AAS LaTeX v. 4.0

Journal ref: Astrophys.J. 505 (1998) 315-338

Molecule Formation In and On Grains.I: Physical Regimes

Authors: Amir Levinson, David F. Chernoff, Edwin E. Salpeter

Abstract: We study molecular hydrogen formation in and on solids. We construct a model with surface sites and bulk sites capable of describing (1) the motion and exchange of H and H_2 between surface and bulk, (2) the recombination of H and dissociation of H_2 in and on the solid, and (3) the injection of H from the gas phase and the loss of H and H_2 from the solid. The basic physical processes include t… ▽ More We study molecular hydrogen formation in and on solids. We construct a model with surface sites and bulk sites capable of describing (1) the motion and exchange of H and H_2 between surface and bulk, (2) the recombination of H and dissociation of H_2 in and on the solid, and (3) the injection of H from the gas phase and the loss of H and H_2 from the solid. The basic physical processes include thermally activated reactions, collisionally induced reactions, and tunneling reactions. Our main application is to the astronomical problem of H_2 formation on grains in space but the model has more general applicability. We investigate the steady-state H and H_2 concentrations in and on the solid when gas phase atoms or ions stick to the surface or penetrate the body of the grain. The model identifies ranges of physical parameters for which the solid becomes saturated (surface and/or bulk) with particles (H and/or H_2) and facilitates the calculation of the efficiency of molecule formation (the fraction of the gas phase atoms that leave as molecules). These solutions are highly degenerate in the sense that they depend only on a small number of dimensionless parameters. We find that a variety of recombination pathways operate under a broad range of conditions. As an example we study H_2 formation in and on carbon grains. △ Less

Submitted 25 June, 1997; originally announced June 1997.

Comments: 28 pages, 7 postscript figures, uses revtex

Neutron Star Population Dynamics.I: Millisecond Pulsars

Authors: J. M. Cordes, David F. Chernoff

Abstract: We study the field millisecond pulsar (MSP) population to infer its intrinsic distribution in spin period and luminosity and to determine its spatial distribution within the Galaxy. Our likelihood analysis on data from extant surveys (22 pulsars with periods <20 ms) accounts for important selection effects. We infer a minimum period cutoff P(min) > 0.65 ms (99% confidence), a period distribution… ▽ More We study the field millisecond pulsar (MSP) population to infer its intrinsic distribution in spin period and luminosity and to determine its spatial distribution within the Galaxy. Our likelihood analysis on data from extant surveys (22 pulsars with periods <20 ms) accounts for important selection effects. We infer a minimum period cutoff P(min) > 0.65 ms (99% confidence), a period distribution proportional to P^{-2.0 +- 0.33} and a pseudo-luminosity distribution proportional to L_p^{-2.0 +- 0.2} (where L_p = flux density times distance^2, for L_p >= 1.1 mJy kpc^2). We find a vertical scale height 0.65{+0.16,-0.12} kpc. We use our results to estimate the total number and birthrate of MSPs in the disk of the Galaxy. We limit the density contribution of a diffuse halo-like component to <1% of the midplane value. The MSP velocity dispersion is smaller that that of young, long-period pulsars by about a factor of 5. Our best estimate of the 1D velocity kick that is unique to MSP evolution is approximately 40 km s^-1. We discuss the evolutionary relationship of MSPs and low-mass X-ray binaries and prospects for future searches for MSPs. △ Less

Submitted 16 June, 1997; originally announced June 1997.

Comments: 53 pages, 12 Postscript figures, LaTeX/AASTeX (aasp4), to appear Ap. J. 482, 1997 June 20

Radiative instabilities in simulations of spherically symmetric supernova blast waves

Authors: Paul A. Kimoto, David F. Chernoff

Abstract: High-resolution simulations of the cooling regions of spherically symmetric supernova remnants demonstrate a strong radiative instability. This instability, whose presence is dependent on the shock velocity, causes large-amplitude fluctuations in the shock velocity. The fluctuations begin almost immediately after the radiative phase begins (upon shell formation) if the shock velocity lies in the… ▽ More High-resolution simulations of the cooling regions of spherically symmetric supernova remnants demonstrate a strong radiative instability. This instability, whose presence is dependent on the shock velocity, causes large-amplitude fluctuations in the shock velocity. The fluctuations begin almost immediately after the radiative phase begins (upon shell formation) if the shock velocity lies in the unstable range; they last until the shock slows to speeds less than approximately 130 km/s. We find that shock-velocity fluctuations from the reverberations of waves within the remnant are small compared to those due to the instability. Further, we find (in plane-parallel simulations) that advected inhomogeneities from the external medium do not interfere with the qualitative nature of the instability-driven fluctuations. Large-amplitude inhomogeneities may alter the phases of shock-velocity fluctuations, but do not substantially reduce their amplitudes. △ Less

Submitted 30 May, 1997; originally announced May 1997.

Comments: 18 pages text, LaTeX/AASTeX (aaspp4); 10 figures; accepted by ApJ

An adaptive hierarchical particle-mesh code with isolated boundary conditions

Authors: Sergio Gelato, David F. Chernoff, Ira Wasserman

Abstract: This article describes a new, fully adaptive Particle-Multiple-Mesh numerical simulation code developed primarily for simulations of small regions (such as a group of galaxies) in a cosmological context. It integrates the equations of motion of a set of particles subject to their mutual gravitational interaction and to an arbitrary external field. The interactions are computed using a hierarchy… ▽ More This article describes a new, fully adaptive Particle-Multiple-Mesh numerical simulation code developed primarily for simulations of small regions (such as a group of galaxies) in a cosmological context. It integrates the equations of motion of a set of particles subject to their mutual gravitational interaction and to an arbitrary external field. The interactions are computed using a hierarchy of nested grids constructed anew at each integration step to enhance the spatial resolution in high-density regions of interest. Significant effort has gone into supporting isolated boundary conditions at the top grid level. This makes our method also applicable to non-cosmological problems, at the cost of some complications which we discuss. We point out the implications of some differences between our approach and those of other authors of similar codes, in particular with respect to the handling of the interface between regions of different spatial resolution. We present a selection of tests performed to verify the correctness and performance of our implementation. The conclusion suggests possible further improvements in the areas of independent time steps and particle softening lengths. △ Less

Submitted 21 November, 1996; v1 submitted 29 July, 1996; originally announced July 1996.

Comments: 35 pages, LaTeX, uses aaspp4.sty, 8 figures included ApJ, in press. Only modest changes from earlier posted version

Report number: SISSA 120/96/A, CRSR 1116

Journal ref: Ap.J. 480 (1997) 115

Truncation of the Binary Distribution Function in Globular Cluster Formation

Authors: E. Vesperini, D. F. Chernoff

Abstract: We investigate a population of primordial binaries during the initial stage of evolution of a star cluster. For our calculations we assume that equal mass stars form rapidly in a tidally truncated gas cloud, that $\sim10$\% of the stars are in binaries and that the resulting star cluster undergoes an epoch of violent relaxation. We study the collisional interaction of the binaries and single sta… ▽ More We investigate a population of primordial binaries during the initial stage of evolution of a star cluster. For our calculations we assume that equal mass stars form rapidly in a tidally truncated gas cloud, that $\sim10$\% of the stars are in binaries and that the resulting star cluster undergoes an epoch of violent relaxation. We study the collisional interaction of the binaries and single stars, in particular, the ionization of the binaries and the energy exchange between binaries and single stars. We find that for large $N$ systems ($N > 10^3$) even the most violent beginning leaves the binary distribution function largely intact. Hence, the binding energy originally tied up in the cloud's protostellar pairs is preserved during the relaxation process and the binaries are available to interact at later times within the virialized cluster. △ Less

Submitted 25 October, 1995; originally announced October 1995.

Comments: 37 pages, LaTeX, 12 compressed postscript figures, accepted by The Astrophysical Journal

Convergence properties of finite-difference hydrodynamics schemes in the presence of shocks

Authors: Paul A. Kimoto, David F. Chernoff

Abstract: We investigate asymptotic convergence in the~$Δx \!\rightarrow\! 0$ limit as a tool for determining whether numerical computations involving shocks are accurate. We use one-dimensional operator-split finite-difference schemes for hydrodynamics with a von Neumann artificial viscosity. An internal-energy scheme converges to demonstrably wrong solutions. We associate this failure with the presence… ▽ More We investigate asymptotic convergence in the~$Δx \!\rightarrow\! 0$ limit as a tool for determining whether numerical computations involving shocks are accurate. We use one-dimensional operator-split finite-difference schemes for hydrodynamics with a von Neumann artificial viscosity. An internal-energy scheme converges to demonstrably wrong solutions. We associate this failure with the presence of discontinuities in the limiting solution. Our extension of the Lax-Wendroff theorem guarantees that certain conservative, operator-split schemes converge to the correct continuum solution. For such a total-energy scheme applied to the formation of a single shock, convergence of a Cauchy error approaches the expected rate slowly. We relate this slowness to the effect of varying diffusion, due to varying linear artificial-viscous length, on small-amplitude waves. In an appendix we discuss the scaling of shock-transition regions with viscous lengths, and exhibit several difficulties for attempts to make extrapolations. △ Less

Submitted 6 November, 1993; originally announced November 1993.

Comments: 29 pages + 21 figures (compressed, uuencoded PostScript, appended), LaTeX 2.09 (AAS preprint substyle v3.0)

Polarized Scattering in the Vicinty of Galaxies

Authors: B. W. Murphy, D. F. Chernoff

Abstract: Some bright cD galaxies in cluster cooling flows have Thomson optical depths exceeding 0.01. A few percent of their luminosity is scattered and appears as diffuse polarized emission. We calculate the scattering process for different geometric combinations of luminosity sources and scattering media. We apply our results to galaxies, with and without active nuclei, immersed in cooling flows. We mo… ▽ More Some bright cD galaxies in cluster cooling flows have Thomson optical depths exceeding 0.01. A few percent of their luminosity is scattered and appears as diffuse polarized emission. We calculate the scattering process for different geometric combinations of luminosity sources and scattering media. We apply our results to galaxies, with and without active nuclei, immersed in cooling flows. We model observations of NGC 1275 and M87 (without active nuclei) in the presence of sky and galactic background fluxes which hinder the measurement of the scattered light at optical wavelengths. Current instruments are unable to detect the scattered light from such objects. However, when a galaxy has an active nucleus of roughly the same luminosity as the remainder of the galaxy in V, both the total and polarized scattered intensity should observable on large scales (5--30kpc), meaning intensity levels greater than 1% of the background level. For typical AGN and galaxy spectral distributions, the scattering is most easily detected at short (U) wavelengths. We point out that a number of such cases will occur. We show that the radiation pattern from the central nuclear region can be mapped using the scattering. We also show that the scattered light can be used to measure inhomogeneities in the cooling flow. △ Less

Submitted 11 June, 1993; originally announced June 1993.

Comments: 29 pages of TEX, 14 figs, CRSR-1046, in ApJ Nov 20, 1993

Gravitational Radiation, Inspiraling Binaries, and Cosmology

Authors: David F. Chernoff, Lee Samuel Finn

Abstract: We show how to measure cosmological parameters using observations of inspiraling binary neutron star or black hole systems in one or more gravitational wave detectors. To illustrate, we focus on the case of fixed mass binary systems observed in a single Laser Interferometer Gravitational-wave Observatory (LIGO)-like detector. Using realistic detector noise estimates, we characterize the rate of… ▽ More We show how to measure cosmological parameters using observations of inspiraling binary neutron star or black hole systems in one or more gravitational wave detectors. To illustrate, we focus on the case of fixed mass binary systems observed in a single Laser Interferometer Gravitational-wave Observatory (LIGO)-like detector. Using realistic detector noise estimates, we characterize the rate of detections as a function of a threshold signal-to-noise ratio $ρ_0$, the Hubble constant $H_0$, and the binary ``chirp'' mass. For $ρ_0 = 8$, $H_0 = 100$ km/s/Mpc, and $1.4 \msun$ neutron star binaries, the sample has a median redshift of $0.22$. Under the same assumptions but independent of $H_0$, a conservative rate density of coalescing binaries ($8\times10^{-8}\,{\rm yr}^{-1}\,{\rm Mpc}^{-3}$) implies LIGO will observe $\sim 50\,{\rm yr}^{-1}$ binary inspiral events. The precision with which $H_0$ and the deceleration parameter $q_0$ may be determined depends on the number of observed inspirals. For fixed mass binary systems, $\sim 100$ observations with $ρ_0 = 10$ in the LIGO detector will give $H_0$ to 10\% in an Einstein-DeSitter cosmology, and 3000 will give $q_0$ to 20\%. For the conservative rate density of coalescing binaries, 100 detections with $ρ_0 = 10$ will require about 4~yrs. △ Less

Submitted 15 April, 1993; originally announced April 1993.

Comments: 13 pages (plus 1 figure), LaTeX/aasms, NU-GR-4

Journal ref: Astrophys.J. 411 (1993) L5-L8

Observing binary inspiral in gravitational radiation: One interferometer

Authors: Lee Samuel Finn, David F. Chernoff

Abstract: We investigate the sensitivity of individual LIGO/VIRGO-like interferometers and the precision with which they can determine the characteristics of an inspiralling binary system. Since the two interferometers of the LIGO detector share nearly the same orientation, their joint sensitivity is similar to that of a single, more sensitive interferometer. We express our results for a single interferom… ▽ More We investigate the sensitivity of individual LIGO/VIRGO-like interferometers and the precision with which they can determine the characteristics of an inspiralling binary system. Since the two interferometers of the LIGO detector share nearly the same orientation, their joint sensitivity is similar to that of a single, more sensitive interferometer. We express our results for a single interferometer of both initial and advanced LIGO design, and also for the LIGO detector in the limit that its two interferometers share exactly the same orientation. We approximate the evolution of a binary system as driven exclusively by leading order quadrupole gravitational radiation. To assess the sensitivity, we calculate the rate at which sources are expected to be observed, the range to which they are observable, and the precision with which characteristic quantities describing the observed binary system can be determined. Assuming a conservative rate density for coalescing neutron star binary systems we expect that the advanced LIGO detector will observe approximately 69~yr${}^{-1}$ with an amplitude SNR greater than 8. Of these, approximately 7~yr${}^{-1}$ will be from binaries at distances greater than 950~Mpc. We explore the sensitivity of these results to a tunable parameter in the interferometer design (the recycling frequency). The optimum choice of the parameter is dependent on the goal of the observations, e.g., maximizing the rate of detections or maximizing the precision of measurement. We determine the optimum parameter values for these two cases. △ Less

Submitted 6 January, 1993; originally announced January 1993.

Comments: 40 pages (plus 7 figures), LaTeX/REVTEX3.0, NU-GR-3

Journal ref: Phys.Rev. D47 (1993) 2198-2219