Source author record

Joseph D. Romano

Joseph D. Romano appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

23 published item(s)

preprint2026arXiv

The NANOGrav 15 yr Data Set: Piecewise Power-Law Reconstruction of the Gravitational-Wave Background

The NANOGrav 15-year (NG15) data set provides evidence for a gravitational-wave background (GWB) signal at nanohertz frequencies, which is expected to originate either from a cosmic population of inspiraling supermassive black-hole binaries or new particle physics in the early Universe. A firm identification of the source of the NG15 signal requires an accurate reconstruction of its frequency spectrum. In this paper, we provide such a spectral characterization of the NG15 signal based on a piecewise power-law (PPL) ansatz that strikes a balance between existing alternatives in the literature. Our PPL reconstruction is more flexible than the standard constant-power-law model, which describes the GWB spectrum in terms of only two parameters: an amplitude A and a spectral index gamma. Concurrently, it better approximates physically realistic GWB spectra -- especially those of cosmological origin -- than the free spectral model, since the latter allows for arbitrary variations in the GWB amplitude from one frequency bin to the next. Our PPL reconstruction of the NG15 signal relies on individual PPL models with a fixed number of internal nodes (i.e., constant power law, broken power law, doubly broken power law, etc.) that are ultimately combined in a Bayesian model average. The data products resulting from our analysis provide the basis for fast refits of spectral GWB models.

preprint2022arXiv

Comparison of maximum likelihood mapping methods for gravitational-wave backgrounds

Detection of a stochastic background of gravitational waves is likely to occur in the next few years. Beyond searches for the isotropic component of SGWBs, there have been various mapping methods proposed to target anisotropic backgrounds. Some of these methods have been applied to data taken by the Laser Interferometer Gravitational-wave Observatories (LIGO) and Virgo. Specifically, these directional searches have focused on mapping the intensity of the signal on the sky via maximum likelihood solutions. We compare this intensity mapping approach to a previously proposed, but never employed, amplitude-phase mapping method to understand whether this latter approach may be employed in future searches. We build up our understanding of the differences between these two approaches by analysing simple toy models of time-stream data, and run mock-data mapping tests for the two methods. We find that the amplitude-phase method is only applicable to the case of a background which is phase-coherent on large scales or, at the very least, has an intrinsic coherence scale that is larger than that of the detector. Otherwise, the amplitude-phase mapping method leads to a loss of overall information, with respect to both phase and amplitude. Since we do not expect these phase-coherent properties to hold for any of the gravitational-wave background signals we hope to detect in the near future, we conclude that intensity mapping is the preferred method for such backgrounds.

preprint2020arXiv

Estimating the angular power spectrum of the gravitational-wave background in the presence of shot noise

There has been much recent interest in studying anisotropies in the astrophysical gravitational-wave (GW) background, as these could provide us with interesting new information about galaxy clustering and large-scale structure. However, this information is obscured by shot noise, caused by the finite number of GW sources that contribute to the background at any given time. We develop a new method for estimating the angular spectrum of anisotropies, based on the principle of combining statistically-independent data segments. We show that this gives an unbiased estimate of the true, astrophysical spectrum, removing the offset due to shot noise power, and that in the limit of many data segments, it is the most efficient (i.e. lowest-variance) estimator possible.

preprint2015arXiv

A Mock Data and Science Challenge for Detecting an Astrophysical Stochastic Gravitational-Wave Background with Advanced LIGO and Advanced Virgo

The purpose of this mock data and science challenge is to prepare the data analysis and science interpretation for the second generation of gravitational-wave experiments Advanced LIGO-Virgo in the search for a stochastic gravitational-wave background signal of astrophysical origin. Here we present a series of signal and data challenges, with increasing complexity, whose aim is to test the ability of current data analysis pipelines at detecting an astrophysically produced gravitational-wave background, test parameter estimation methods and interpret the results. We introduce the production of these mock data sets that includes a realistic observing scenario data set where we account for different sensitivities of the advanced detectors as they are continuously upgraded toward their design sensitivity. After analysing these with the standard isotropic cross-correlation pipeline we find that we are able to recover the injected gravitational-wave background energy density to within $2σ$ for all of the data sets and present the results from the parameter estimation. The results from this mock data and science challenge show that advanced LIGO and Virgo will be ready and able to make a detection of an astrophysical gravitational-wave background within a few years of operations of the advanced detectors, given a high enough rate of compact binary coalescing events.

preprint2015arXiv

Mapping gravitational-wave backgrounds in modified theories of gravity using pulsar timing arrays

We extend our previous work on applying CMB techniques to the mapping of gravitational-wave backgrounds to backgrounds which have non-GR polarisations. Our analysis and results are presented in the context of pulsar-timing array observations, but the overarching methods are general, and can be easily applied to LIGO or eLISA observations using appropriately modified response functions. Analytic expressions for the pulsar-timing response to gravitational waves with non-GR polarisation are given for each mode of a spin-weighted spherical-harmonic decomposition of the background, which permit the signal to be mapped across the sky to any desired resolution. We also derive the pulsar-timing overlap reduction functions for the various non-GR polarisations, finding analytic forms for anisotropic backgrounds with scalar-transverse ("breathing") and vector-longitudinal polarisations, and a semi-analytic form for scalar-longitudinal backgrounds. Our results indicate that pulsar-timing observations will be completely insensitive to scalar-transverse mode anisotropies in the polarisation amplitude beyond dipole, and anisotropies in the power beyond quadrupole. Analogously to our previous findings that pulsar-timing observations lack sensitivity to tensor-curl modes for a transverse-traceless tensor background, we also find insensitivity to vector-curl modes for a vector-longitudinal background.

preprint2015arXiv

Phase-coherent mapping of gravitational-wave backgrounds using ground-based laser interferometers

We extend the formalisms developed in Gair et al. and Cornish and van Haasteren to create maps of gravitational-wave backgrounds using a network of ground-based laser interferometers. We show that in contrast to pulsar timing arrays, which are insensitive to half of the gravitational-wave sky (the curl modes), a network of ground-based interferometers is sensitive to both the gradient and curl components of the background. The spatial separation of a network of interferometers, or of a single interferometer at different times during its rotational and orbital motion around the Sun, allows for recovery of both components. We derive expressions for the response functions of a laser interferometer in the small-antenna limit, and use these expressions to calculate the overlap reduction function for a pair of interferometers. We also construct maximum-likelihood estimates of the + and x-polarization modes of the gravitational-wave sky in terms of the response matrix for a network of ground-based interferometers, evaluated at discrete times during Earth's rotational and orbital motion around the Sun. We demonstrate the feasibility of this approach for some simple simulated backgrounds (a single point source and spatially-extended distributions having only grad or curl components), calculating maximum-likelihood sky maps and uncertainty maps based on the (pseudo)inverse of the response matrix. The distinction between this approach and standard methods for mapping gravitational-wave power is also discussed.

preprint2015arXiv

Understanding the gravitational-wave Hellings and Downs curve for pulsar timing arrays in terms of sound and electromagnetic waves

Searches for stochastic gravitational-wave backgrounds using pulsar timing arrays look for correlations in the timing residuals induced by the background across the pulsars in the array. The correlation signature of an isotropic, unpolarized gravitational-wave background predicted by general relativity follows the so-called Hellings and Downs curve, which is a relatively simple function of the angle between a pair of Earth-pulsar baselines. In this paper, we give a pedagogical discussion of the Hellings and Downs curve for pulsar timing arrays, considering simpler analogous scenarios involving sound and electromagnetic waves. We calculate Hellings-and-Downs-type functions for these two scenarios and develop a framework suitable for doing more general correlation calculations.

preprint2015arXiv

When is a gravitational-wave signal stochastic?

We discuss the detection of gravitational-wave backgrounds in the context of Bayesian inference and suggest a practical definition of what it means for a signal to be considered stochastic---namely, that the Bayesian evidence favors a stochastic signal model over a deterministic signal model. A signal can further be classified as Gaussian-stochastic if a Gaussian signal model is favored. In our analysis we use Bayesian model selection to choose between several signal and noise models for simulated data consisting of uncorrelated Gaussian detector noise plus a superposition of sinusoidal signals from an astrophysical population of gravitational-wave sources. For simplicity, we consider co-located and co-aligned detectors with white detector noise, but the method can be extended to more realistic detector configurations and power spectra. The general trend we observe is that a deterministic model is favored for small source numbers, a non-Gaussian stochastic model is preferred for intermediate source numbers, and a Gaussian stochastic model is preferred for large source numbers. However, there is very large variation between individual signal realizations, leading to fuzzy boundaries between the three regimes. We find that a hybrid, trans-dimensional model comprised of a deterministic signal model for individual bright sources and a Gaussian-stochastic signal model for the remaining confusion background outperforms all other models in most instances.

preprint2014arXiv

Mapping gravitational-wave backgrounds using methods from CMB analysis: Application to pulsar timing arrays

We describe an alternative approach to the analysis of gravitational-wave backgrounds, based on the formalism used to characterise the polarisation of the cosmic microwave background. In contrast to standard analyses, this approach makes no assumptions about the nature of the background and so has the potential to reveal much more about the physical processes that generated it. An arbitrary background can be decomposed into modes whose angular dependence on the sky is given by gradients and curls of spherical harmonics. We derive the pulsar timing overlap reduction functions for the individual modes, which are given by simple combinations of spherical harmonics evaluated at the pulsar locations. We show how these can be used to recover the components of an arbitrary background, giving explicit results for both isotropic and anisotropic uncorrelated backgrounds. We also find that the response of a pulsar timing array to curl modes is identically zero, so half of the gravitational-wave sky will never be observed using pulsar timing, no matter how many pulsars are included in the array. An isotropic, unpolarised and uncorrelated background can be accurately represented using only three modes, and so a search of this type will be only slightly more complicated than the standard cross-correlation search using the Hellings and Downs overlap reduction function. However, by measuring the components of individual modes of the background and checking for consistency with isotropy, this approach has the potential to reveal much more information. Each individual mode on its own describes a background that is correlated between different points on the sky. A measurement of the components that indicates the presence of correlations in the background on large angular scales would suggest startling new physics.

preprint2014arXiv

Time-domain Implementation of the Optimal Cross-Correlation Statistic for Stochastic Gravitational-Wave Background Searches in Pulsar Timing Data

Supermassive black hole binaries, cosmic strings, relic gravitational waves from inflation, and first order phase transitions in the early universe are expected to contribute to a stochastic background of gravitational waves in the 10^(-9) Hz-10^(-7) Hz frequency band. Pulsar timing arrays (PTAs) exploit the high precision timing of radio pulsars to detect signals at such frequencies. Here we present a time-domain implementation of the optimal cross-correlation statistic for stochastic background searches in PTA data. Due to the irregular sampling typical of PTA data as well as the use of a timing model to predict the times-of-arrival of radio pulses, time-domain methods are better suited for gravitational wave data analysis of such data. We present a derivation of the optimal cross-correlation statistic starting from the likelihood function, a method to produce simulated stochastic background signals, and a rigorous derivation of the scaling laws for the signal-to-noise ratio of the cross-correlation statistic in the two relevant PTA regimes: the weak signal limit where instrumental noise dominates over the gravitational wave signal at all frequencies, and a second regime where the gravitational wave signal dominates at the lowest frequencies.

preprint2013arXiv

A Rømer time-delay determination of the gravitational-wave propagation speed

In 1676 Olaus Rømer presented the first observational evidence for a finite light velocity $\cem$. He formed his estimate by attributing the periodically varying discrepancy between the observed and expected occultation times of the Galilean satellite Io by its planetary host Jupiter to the time it takes light to cross Earth's orbital diameter. Given a stable celestial clock that can be observed in gravitational waves the same principle can be used to measure the propagation speed $\cgw$ of gravitational radiation. Space-based "LISA"-like detectors will, and terrestrial LIGO-like detectors may, observe such clocks and thus be capable of directly measuring the propagation velocity of gravitational waves. In the case of space-based detectors the clocks will be galactic close white dwarf binary systems; in the case of terrestrial detectors, the most likely candidate clock is the periodic gravitational radiation from a rapidly rotating non-axisymmetric neutron star. Here we evaluate the accuracy that may be expected of such a Rømer-type measurement of $\cgw$ by foreseeable future space-based and terrestrial detectors. For space-based, LISA-like detectors, periodic sources are plentiful: by the end of the first year of scientific operations a LISA-like detector will have measured $\cgw$ to better than a part in a thousand. Periodic sources may not be accessible in terrestrial detectors available to us in the foreseeable future; however, if such a source is detected then with a year of observations we could measure $\cgw$ to better than a part in a million.

preprint2013arXiv

The stochastic background: scaling laws and time to detection for pulsar timing arrays

We derive scaling laws for the signal-to-noise ratio of the optimal cross-correlation statistic, and show that the large power-law increase of the signal-to-noise ratio as a function of the the observation time $T$ that is usually assumed holds only at early times. After enough time has elapsed, pulsar timing arrays enter a new regime where the signal to noise only scales as $\sqrt{T}$. In addition, in this regime the quality of the pulsar timing data and the cadence become relatively un-important. This occurs because the lowest frequencies of the pulsar timing residuals become gravitational-wave dominated. Pulsar timing arrays enter this regime more quickly than one might naively suspect. For T=10 yr observations and typical stochastic background amplitudes, pulsars with residual RMSs of less than about $1\,μ$s are already in that regime. The best strategy to increase the detectability of the background in this regime is to increase the number of pulsars in the array. We also perform realistic simulations of the NANOGrav pulsar timing array, which through an aggressive pulsar survey campaign adds new millisecond pulsars regularly to its array, and show that a detection is possible within a decade, and could occur as early as 2016.

preprint2013arXiv

Towards a unified treatment of gravitational-wave data analysis

We present a unified description of gravitational-wave data analysis that unites the template-based analysis used to detect deterministic signals from well-modeled sources, such as binary-black-hole mergers, with the cross-correlation analysis used to detect stochastic gravitational-wave backgrounds. We also discuss the connection between template-based analyses and those that target poorly-modeled bursts of gravitational waves, and suggest a new approach for detecting burst signals.

preprint2012arXiv

In an expanding universe, what doesn't expand?

The expansion of the universe is often viewed as a uniform stretching of space that would affect compact objects, atoms and stars, as well as the separation of galaxies. One usually hears that bound systems do not take part in the general expansion, but a much more subtle question is whether bound systems expand partially. In this paper, a very definitive answer is given for a very simple system: a classical "atom" bound by electrical attraction. With a mathemical description appropriate for undergraduate physics majors, we show that this bound system either completely follows the cosmological expansion, or -- after initial transients -- completely ignores it. This "all or nothing" behavior can be understood with techniques of junior-level mechanics. Lastly, the simple description is shown to be a justifiable approximation of the relativistically correct formulation of the problem.

preprint2012arXiv

Treatment of Calibration Uncertainty in Multi-Baseline Cross-Correlation Searches for Gravitational Waves

Uncertainty in the calibration of gravitational-wave (GW) detector data leads to systematic errors which must be accounted for in setting limits on the strength of GW signals. When cross-correlation measurements are made using data from a pair of instruments, as in searches for a stochastic GW background, the calibration uncertainties of the individual instruments can be combined into an uncertainty associated with the pair. With the advent of multi-baseline GW observation (e.g., networks consisting of multiple detectors such as the LIGO observatories and Virgo), a more sophisticated treatment is called for. We describe how the correlations between calibration factors associated with different pairs can be taken into account by marginalizing over the uncertainty associated with each instrument.

preprint2009arXiv

Probing the anisotropies of a stochastic gravitational-wave background using a network of ground-based laser interferometers

We present a maximum-likelihood analysis for estimating the angular distribution of power in an anisotropic stochastic gravitational-wave background using ground-based laser interferometers. The standard isotropic and gravitational-wave radiometer searches (optimal for point sources) are recovered as special limiting cases. The angular distribution can be decomposed with respect to any set of basis functions on the sky, and the single-baseline, cross-correlation analysis is easily extended to a network of three or more detectors-that is, to multiple baselines. A spherical harmonic decomposition, which provides maximum-likelihood estimates of the multipole moments of the gravitational-wave sky, is described in detail. We also discuss: (i) the covariance matrix of the estimators and its relationship to the detector response of a network of interferometers, (ii) a singular-value decomposition method for regularizing the deconvolution of the detector response from the measured sky map, (iii) the expected increase in sensitivity obtained by including multiple baselines, and (iv) the numerical results of this method when applied to simulated data consisting of both point-like and di#use sources. Comparisions between this general method and the standard isotropic and radiometer searches are given throughout, to make contact with the existing literature on stochastic background searches.

preprint1999arXiv

Detecting an association between Gamma Ray and Gravitational Wave Bursts

If $γ$-ray bursts (GRBs) are accompanied by gravitational wave bursts (GWBs) the correlated output of two gravitational wave detectors evaluated in the moments just prior to a GRB will differ from that evaluated at times not associated with a GRB. We can test for this difference independently of any model of the GWB signal waveform. If we invoke a model for the GRB source population and GWB radiation spectral density we can find a confidence interval or upper limit on the root-mean-square GWB signal amplitude in the detector waveband. To illustrate we adopt a simple, physically motivated model and estimate that initial LIGO detector observations coincident with 1000 GRBs could lead us to exclude, with 95% confidence, associated GWBs with $h_{RMS} \gtrsim 1.7 \times 10^{-22}$. This result does not require the detector noise be Gaussian or that any inter-detector correlated noise be measured or measurable; it does not require advanced or a priori knowledge of the source waveform; and the limits obtained on the wave-strength improve with the number of observed GRBs.

preprint1999arXiv

Quasistationary binary inspiral. I. Einstein equations for the two Killing vector spacetime

The geometry of two infinitely long lines of mass moving in a fixed circular orbit is considered as a toy model for the inspiral of a binary system of compact objects due to gravitational radiation. The two Killing fields in the toy model are used, according to a formalism introduced by Geroch, to describe the geometry entirely in terms of a set of tensor fields on the two-manifold of Killing vector orbits. Geroch's derivation of the Einstein equations in this formalism is streamlined and generalized. The explicit Einstein equations for the toy model spacetime are derived in terms of the degrees of freedom which remain after a particular choice of gauge.

preprint1996arXiv

Dirac Constraint Quantization of a Dilatonic Model of Gravitational Collapse

We present an anomaly-free Dirac constraint quantization of the string-inspired dilatonic gravity (the CGHS model) in an open 2-dimensional spacetime. We show that the quantum theory has the same degrees of freedom as the classical theory; namely, all the modes of the scalar field on an auxiliary flat background, supplemented by a single additional variable corresponding to the primordial component of the black hole mass. The functional Heisenberg equations of motion for these dynamical variables and their canonical conjugates are linear, and they have exactly the same form as the corresponding classical equations. A canonical transformation brings us back to the physical geometry and induces its quantization.

preprint1995arXiv

Gravitational Constraints which Generate a Lie Algebra

The coupling of gravity to dust helps to discover simple quadratic combinations of the gravitational super-Hamiltonian and supermomentum whose Poisson brackets strongly vanish. This leads to a new form of vacuum constraints which generate a true Lie algebra. We show that the coupling of gravity to a massless scalar field leads to yet another set of constraints with the same property, albeit not as simple as that based on the coupling to dust.

preprint1995arXiv

Spherically Symmetric Scalar Field Collapse: An Example of the Spacetime Problem of Time

A canonical formalism for spherical symmetry, originally developed by Kuchař to describe vacuum Schwarzschild black holes, is extended to include a spherically symmetric, massless, scalar field source. By introducing the ADM mass as a canonical coordinate on phase space, one finds that the super-Hamiltonian and supermomentum constraints for the coupled system simplify considerably. Yet, despite this simplification, it is difficult to find a functional time formalism for the theory. First, the configuration variable that played the role of time for the vacuum theory is no longer a spacetime scalar once spherically symmetric matter is coupled to gravity. Second, although it is possible to perform a canonical transformation to a new set of variables in terms of which the super-Hamiltonian and supermomentum constraints can be solved, the new time variable also fails to be a spacetime scalar. As such, our solutions suffer from the so-called {\it spacetime problem of time}. A candidate for a time variable that {\it is} a spacetime scalar is presented. Problems with turning this variable into a canonical coordinate on phase space are discussed.

preprint1992arXiv

Degenerate Extensions of General Relativity

General relativity has previously been extended to incorporate degenerate metrics using Ashtekar's hamiltonian formulation of the theory. In this letter, we show that a natural alternative choice for the form of the hamiltonian constraints leads to a theory which agrees with GR for non-degenerate metrics, but differs in the degenerate sector from Ashtekar's original degenerate extension. The Poisson bracket algebra of the alternative constraints closes in the non-degenerate sector, with structure functions that involve the {\it inverse} of the spatial triad. Thus, the algebra does {\it not} close in the degenerate sector. We find that it must be supplemented by an infinite number ofsecondary constraints, which are shown to be first class (although their explicit form is not worked out in detail). All of the constraints taken together are implied by, but do not imply, Ashtekar's original form of constraints. Thus, the alternative constraints give rise to a different degenerate extension of GR. In the corresponding quantum theory, the single loop and intersecting loop holonomy states found in the connection representation satisfy {\it all} of the constraints. These states are therefore exact (formal) solutions to this alternative degenerate extension of quantum gravity, even though they are {\it not} solutions to the usual vector constraint.