Source author record

John T. Whelan

John T. Whelan 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

15works
9topics
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

15 published item(s)

preprint2021arXiv

Bradley-Terry Modeling with Multiple Game Outcomes with Applications to College Hockey

The Bradley-Terry model has previously been used in both Bayesian and frequentist interpretations to evaluate the strengths of sports teams based on win-loss game results. It has also been extended to handle additional possible results such as ties. We implement a generalization which includes multiple possible outcomes such as wins or losses in regulation, overtime, or shootouts. A natural application is to ice hockey competitions such as international matches, European professional leagues, and NCAA hockey, all of which use a zero-sum point system which values overtime and shootout wins as 1/3 of a win, and overtime and shootout losses as 1/3 of a win. We incorporate this into the probability model, and evaluate the posterior distributions for the associated strength parameters using techniques such as Gaussian expansion about maximum a posteriori estimates, and Hamiltonian Monte Carlo.

preprint2021arXiv

Template Lattices for a Cross-Correlation Search for Gravitational Waves from Scorpius X-1

We describe the application of the lattice covering problem to the placement of templates in a search for continuous gravitational waves from the low-mass X-Ray binary Scorpius X-1. Efficient placement of templates to cover the parameter space at a given maximum mismatch is an application of the sphere covering problem, for which an implementation is available in the LatticeTiling software library. In the case of Sco X-1, potential correlations, in both the prior uncertainty and the mismatch metric, between the orbital period and orbital phase, lead to complications in the efficient construction of the lattice. We define a shearing coordinate transformation which simultaneously minimizes both of these sources of correlation, and allows us to take advantage of the small prior orbital period uncertainty. The resulting lattices have a factor of about 3 fewer templates than the corresponding parameter space grids constructed by the prior straightforward method, allowing a more sensitive search at the same computing cost and maximum mismatch.

preprint2015arXiv

Constraining star formation rates in cool-core brightest cluster galaxies

We used broad-band imaging data for 10 cool-core brightest cluster galaxies (BCGs) and conducted a Bayesian analysis using stellar population synthesis to determine the likely properties of the constituent stellar populations. Determination of ongoing star formation rates (SFRs), in particular, has a direct impact on our understanding of the cooling of the intracluster medium (ICM), star formation and AGN-regulated feedback. Our model consists of an old stellar population and a series of young stellar components. We calculated marginalized posterior probability distributions for various model parameters and obtained 68% plausible intervals from them. The 68% plausible interval on the SFRs is broad, owing to a wide range of models that are capable of fitting the data, which also explains the wide dispersion in the star formation rates available in the literature. The ranges of possible SFRs are robust and highlight the strength in such a Bayesian analysis. The SFRs are correlated with the X-ray mass deposition rates (the former are factors of 4 to 50 lower than the latter), implying a picture where the cooling of the ICM is a contributing factor to star formation in cool-core BCGs. We find that 9 out of 10 BCGs have been experiencing starbursts since 6 Gyr ago. While four out of 9 BCGs seem to require continuous SFRs, 5 out of 9 seem to require periodic star formation on intervals ranging from 20 Myr to 200 Myr. This time scale is similar to the cooling-time of the ICM in the central (< 5 kpc) regions.

preprint2015arXiv

Model-Based Cross-Correlation Search for Gravitational Waves from Scorpius X-1

We consider the cross-correlation search for periodic GWs and its potential application to the LMXB Sco X-1. This method coherently combines data from different detectors at the same time, as well as different times from the same or different detectors. By adjusting the maximum time offset between a pair of data segments to be coherently combined, one can tune the method to trade off sensitivity and computing costs. In particular, the detectable signal amplitude scales as the inverse fourth root of this coherence time. The improvement in amplitude sensitivity for a search with a coherence time of 1hr, compared with a directed stochastic background search with 0.25Hz wide bins is about a factor of 5.4. We show that a search of 1yr of data from Advanced LIGO and Advanced Virgo with a coherence time of 1hr would be able to detect GWs from Sco X-1 at the level predicted by torque balance over a range of signal frequencies from 30-300Hz; if the coherence time could be increased to 10hr, the range would be 20-500Hz. In addition, we consider several technical aspects of the cross-correlation method: We quantify the effects of spectral leakage and show that nearly rectangular windows still lead to the most sensitive search. We produce an explicit parameter-space metric for the cross-correlation search in general and as applied to a neutron star in a circular binary system. We consider the effects of using a signal template averaged over unknown amplitude parameters: the search is sensitive to a combination of the intrinsic signal amplitude and the inclination of the neutron star rotation axis, and the peak of the expected detection statistic is systematically offset from the true signal parameters. Finally, we describe the potential loss of SNR due to unmodelled effects such as signal phase acceleration within the Fourier transform timescale and gradual evolution of the spin frequency.

preprint2014arXiv

New Coordinates for the Amplitude Parameter Space of Continuous Gravitational Waves

The parameter space for continuous gravitational waves (GWs) can be divided into amplitude parameters (signal amplitude, inclination and polarization angles describing the orientation of the source, and an initial phase) and phase-evolution parameters. The division is useful in part because the Jaranowski-Krolak-Schutz (JKS) coordinates on the four-dimensional amplitude parameter space allow the GW signal to be written as a linear combination of four template waveforms with the JKS coordinates as coefficients. We define a new set of coordinates on the amplitude parameter space, with the same properties, which is more closely connected to the physical amplitude parameters. These naturally divide into two pairs of Cartesian-like coordinates on two-dimensional subspaces, one corresponding to left- and the other to right-circular polarization. We thus refer to these as CPF (circular polarization factored) coordinates. The corresponding two sets of polar coordinates (known as CPF-polar) can be related in a simple way to the physical parameters. We illustrate some simplifying applications for these various coordinate systems, such as: a calculation of Jacobians between various coordinate systems; an illustration of the signal coordinate singularities associated with left- and right-circular polarization, which correspond to the origins of the two two-dimensional subspaces; and an elucidation of the form of the log-likelihood ratio between hypotheses of Gaussian noise with and without a continuous GW signal. These are used to illustrate some of the prospects for approximate evaluation of a Bayesian detection statistic defined by marginalization over the physical parameter space. Additionally, in the presence of simplifying assumptions about the observing geometry, we are able to explicitly evaluate the integral for the Bayesian detection statistic, and compare it to the approximate results.

preprint2013arXiv

Improving the sensitivity of a search for coalescing binary black holes with non-precessing spins in gravitational wave data

We demonstrate the implementation of a sensitive search pipeline for gravitational waves from coalescing binary black holes whose components have spins aligned with the orbital angular momentum. We study the pipeline recovery of simulated gravitational wave signals from aligned-spin binary black holes added to real detector noise, comparing the pipeline performance with aligned-spin filter templates to the same pipeline with non-spinning filter templates. Our results exploit a three-parameter phenomenological waveform family that models the full inspiral-merger-ringdown coalescence and treats the effect of aligned spins with a single effective spin parameter χ. We construct template banks from these waveforms by a stochastic placement method and use these banks as filters in the recently-developed gstlal search pipeline. We measure the observable volume of the analysis pipeline for binary black hole signals with total mass in [15,25] solar masses and χin [0, 0.85]. We find an increase in observable volume of up to 45% for systems with 0.2 <= χ<= 0.85 with almost no loss of sensitivity to signals with 0 <= χ<= 0.2. We demonstrate this analysis on 25.9 days of data obtained from the Hanford and Livingston detectors in LIGO's fifth observation run.

preprint2012arXiv

Herschel observations of extended atomic gas in the core of the Perseus cluster

We present Herschel observations of the core of the Perseus cluster of galaxies. The brightest cluster galaxy, NGC 1275, is surrounded by a network of filaments previously imaged extensively in Hα and CO. In this work, we report detections of FIR lines with Herschel. All but one of the lines are spatially extended, with the [CII] line emission extending up to 25 kpc from the core. There is spatial and kinematical correlation among [CII], Hα and CO, which gives us confidence to model the different components of the gas with a common heating model. With the help of FIR continuum Herschel measurements, together with a suite of coeval radio, submm and infrared data, we performed a SED fitting of NGC 1275 using a model that contains contributions from dust emission as well as synchrotron AGN emission. The data indicate a low dust emissivity index, beta ~ 1, a total dust mass close to 10^7 solar mass, a cold dust component with temperature 38 \pm 2 K and a warm dust component with temperature of 116 \pm 9 K. The FIR-derived star formation rate (SFR) is 24 \pm 1 solar mass per yr, in close agreement with the FUV-derived SFR. We investigated in detail the source of the Herschel FIR and Hα emissions emerging from a core region 4 kpc in radius. Based on simulations conducted using the radiative transfer code, CLOUDY, a heating model comprising old and young stellar populations is sufficient to explain these observations. We have also detected [CII] in three well-studied regions of the filaments. We find a [OI]/[CII] ratio about 1 dex smaller than predicted by the otherwise functional Ferland (2009) model. The line ratio suggests that the lines are optically thick, as is typical of galactic PDRs, and implies that there is a large reservoir of cold atomic gas. [abridged]

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.

preprint2011arXiv

Designing a cross-correlation search for continuous-wave gravitational radiation from a neutron star in the supernova remnant SNR 1987A

A strategy is devised for a semi-coherent cross-correlation search for a young neutron star in the supernova remnant SNR 1987A, using science data from the Initial LIGO and/or Virgo detectors. An astrophysical model for the gravitational wave phase is introduced which describes the star's spin down in terms of its magnetic field strength $B$ and ellipticity $ε$, instead of its frequency derivatives. The model accurately tracks the gravitational wave phase from a rapidly decelerating neutron star under the restrictive but computationally unavoidable assumption of constant braking index, an issue which has hindered previous searches for such young objects. The theoretical sensitivity is calculated and compared to the indirect, age-based wave strain upper limit. The age-based limit lies above the detection threshold in the frequency band 75\,Hz $\lesssim ν\lesssim 450$\,Hz. The semi-coherent phase metric is also calculated and used to estimate the optimal search template spacing for the search. The range of search parameters that can be covered given our computational resources ($\sim 10^9$ templates) is also estimated. For Initial LIGO sensitivity, in the frequency band between 50\,Hz and 500\,Hz, in the absence of a detected signal, we should be able to set limits of $B \gtrsim 10^{11}$\,G and $ε\lesssim 10^{-4}$.

preprint2010arXiv

Searching for Galactic White Dwarf Binaries in Mock LISA Data using an F-Statistic Template Bank

We describe an F-statistic search for continuous gravitational waves from galactic white-dwarf binaries in simulated LISA Data. Our search method employs a hierarchical template-grid based exploration of the parameter space. In the first stage, candidate sources are identified in searches using different simulated laser signal combinations (known as TDI variables). Since each source generates a primary maximum near its true "Doppler parameters" (intrinsic frequency and sky position) as well as numerous secondary maxima of the F-statistic in Doppler parameter space, a search for multiple sources needs to distinguish between true signals and secondary maxima associated with other, "louder" signals. Our method does this by applying a coincidence test to reject candidates which are not found at nearby parameter space positions in searches using each of the three TDI variables. For signals surviving the coincidence test, we perform a fully coherent search over a refined parameter grid to provide an accurate parameter estimation for the final candidates. Suitably tuned, the pipeline is able to extract 1989 true signals with only 5 false alarms. The use of the rigid adiabatic approximation allows recovery of signal parameters with errors comparable to statistical expectations, although there is still some systematic excess with respect to statistical errors expected from Gaussian noise. An experimental iterative pipeline with seven rounds of signal subtraction and re-analysis of the residuals allows us to increase the number of signals recovered to a total of 3419 with 29 false alarms.

preprint2010arXiv

The Mock LISA Data Challenges: from Challenge 3 to Challenge 4

The Mock LISA Data Challenges are a program to demonstrate LISA data-analysis capabilities and to encourage their development. Each round of challenges consists of one or more datasets containing simulated instrument noise and gravitational waves from sources of undisclosed parameters. Participants analyze the datasets and report best-fit solutions for the source parameters. Here we present the results of the third challenge, issued in Apr 2008, which demonstrated the positive recovery of signals from chirping Galactic binaries, from spinning supermassive--black-hole binaries (with optimal SNRs between ~ 10 and 2000), from simultaneous extreme-mass-ratio inspirals (SNRs of 10-50), from cosmic-string-cusp bursts (SNRs of 10-100), and from a relatively loud isotropic background with Omega_gw(f) ~ 10^-11, slightly below the LISA instrument noise.

preprint2009arXiv

Status of NINJA: the Numerical INJection Analysis project

The 2008 NRDA conference introduced the Numerical INJection Analysis project (NINJA), a new collaborative effort between the numerical relativity community and the data analysis community. NINJA focuses on modeling and searching for gravitational wave signatures from the coalescence of binary system of compact objects. We review the scope of this collaboration and the components of the first NINJA project, where numerical relativity groups shared waveforms and data analysis teams applied various techniques to detect them when embedded in colored Gaussian noise.

preprint2009arXiv

Testing gravitational-wave searches with numerical relativity waveforms: Results from the first Numerical INJection Analysis (NINJA) project

The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave data analysis communities. The purpose of NINJA is to study the sensitivity of existing gravitational-wave search algorithms using numerically generated waveforms and to foster closer collaboration between the numerical relativity and data analysis communities. We describe the results of the first NINJA analysis which focused on gravitational waveforms from binary black hole coalescence. Ten numerical relativity groups contributed numerical data which were used to generate a set of gravitational-wave signals. These signals were injected into a simulated data set, designed to mimic the response of the Initial LIGO and Virgo gravitational-wave detectors. Nine groups analysed this data using search and parameter-estimation pipelines. Matched filter algorithms, un-modelled-burst searches and Bayesian parameter-estimation and model-selection algorithms were applied to the data. We report the efficiency of these search methods in detecting the numerical waveforms and measuring their parameters. We describe preliminary comparisons between the different search methods and suggest improvements for future NINJA analyses.

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.

preprint1995arXiv

Generalized quantum mechanics of nonabelian gauge theories

Hartle's generalized quantum mechanics in the sum-over-histories formalism is used to describe a nonabelian gauge theory. Predictions are made for certain alternatives, with particular attention given to coarse-grainings involving the constraint. In this way, the theory is compared to other quantum-mechanical descriptions of gauge theories in which the constraints are imposed by hand. The vanishing of the momentum space constraint is seen to hold, both through a simple formal argument and via a more careful description of the Lorentzian path integral as defined on a spacetime lattice. (Incidentally, the treatment of the time slicing in the path integral may be of general technical interest.) The configuration space realization of the constraint is shown to behave in a more complicated fashion. For some coarse grainings, we recover the known result from an abelian theory, that coarse grainings by values of the constraint either predict its vanishing or fail to decohere. However, sets of alternatives defined in terms of a more complicated quantity in the abelian case are exhibited where definite predictions can be made which disagree with the assumption that the constraints vanish. Finally, the configuration space sum-over-histories theory is exhibited in a manifestly Lorentz-invariant formulation.