Showing 1–11 of 11 results for author: Barrett, J W

Rapid stellar and binary population synthesis with COMPAS

Authors: Team COMPAS, :, Jeff Riley, Poojan Agrawal, Jim W. Barrett, Kristan N. K. Boyett, Floor S. Broekgaarden, Debatri Chattopadhyay, Sebastian M. Gaebel, Fabian Gittins, Ryosuke Hirai, George Howitt, Stephen Justham, Lokesh Khandelwal, Floris Kummer, Mike Y. M. Lau, Ilya Mandel, Selma E. de Mink, Coenraad Neijssel, Tim Riley, Lieke van Son, Simon Stevenson, Alejandro Vigna-Gomez, Serena Vinciguerra, Tom Wagg , et al. (1 additional authors not shown)

Abstract: Compact Object Mergers: Population Astrophysics and Statistics (COMPAS; https://compas.science) is a public rapid binary population synthesis code. COMPAS generates populations of isolated stellar binaries under a set of parametrized assumptions in order to allow comparisons against observational data sets, such as those coming from gravitational-wave observations of merging compact remnants. It i… ▽ More Compact Object Mergers: Population Astrophysics and Statistics (COMPAS; https://compas.science) is a public rapid binary population synthesis code. COMPAS generates populations of isolated stellar binaries under a set of parametrized assumptions in order to allow comparisons against observational data sets, such as those coming from gravitational-wave observations of merging compact remnants. It includes a number of tools for population processing in addition to the core binary evolution components. COMPAS is publicly available via the github repository https://github.com/TeamCOMPAS/COMPAS/, and is designed to allow for flexible modifications as evolutionary models improve. This paper describes the methodology and implementation of COMPAS. It is a living document which will be updated as new features are added to COMPAS; the current document describes COMPAS v02.21.00. △ Less

Submitted 28 December, 2021; v1 submitted 20 September, 2021; originally announced September 2021.

Comments: Code publicly available via https://compas.science . Minor updates to match version accepted to ApJS

The effect of the metallicity-specific star formation history on double compact object mergers

Authors: Coenraad J. Neijssel, Alejandro Vigna-Gómez, Simon Stevenson, Jim W. Barrett, Sebastian M. Gaebel, Floor Broekgaarden, Selma E. de Mink, Dorottya Szécsi, Serena Vinciguerra, Ilya Mandel

Abstract: We investigate the impact of uncertainty in the metallicity-specific star formation rate over cosmic time on predictions of the rates and masses of double compact object mergers observable through gravitational waves. We find that this uncertainty can change the predicted detectable merger rate by more than an order of magnitude, comparable to contributions from uncertain physical assumptions rega… ▽ More We investigate the impact of uncertainty in the metallicity-specific star formation rate over cosmic time on predictions of the rates and masses of double compact object mergers observable through gravitational waves. We find that this uncertainty can change the predicted detectable merger rate by more than an order of magnitude, comparable to contributions from uncertain physical assumptions regarding binary evolution, such as mass transfer efficiency or supernova kicks. We statistically compare the results produced by the COMPAS population synthesis suite against a catalog of gravitational-wave detections from the first two Advanced LIGO and Virgo observing runs. We find that the rate and chirp mass of observed binary black hole mergers can be well matched under our default evolutionary model with a star formation metallicity spread of $0.39$ dex around a mean metallicity $\left<Z\right>$ that scales with redshift $z$ as $\left<Z\right>=0.035 \times 10^{-0.23 z}$, assuming a star formation rate of $0.01 \times (1+z)^{2.77} / (1+((1+z)/2.9)^{4.7}) \, \rm{M}_\odot$ Mpc$^{-3}$ yr$^{-1}$. Intriguingly, this default model predicts that 80\% of the approximately one binary black hole merger per day that will be detectable at design sensitivity will have formed through isolated binary evolution with only dynamically stable mass transfer, i.e., without experiencing a common-envelope event. △ Less

Submitted 19 June, 2019; originally announced June 2019.

Comments: 22 pages, 17 figures. Data from the COMPAS simulation will be publicly available

STROOPWAFEL: Simulating rare outcomes from astrophysical populations, with application to gravitational-wave sources

Authors: Floor S. Broekgaarden, Stephen Justham, Selma E. de Mink, Jonathan Gair, Ilya Mandel, Simon Stevenson, Jim W. Barrett, Alejandro Vigna-Gómez, Coenraad J. Neijssel

Abstract: Gravitational-wave observations of double compact object (DCO) mergers are providing new insights into the physics of massive stars and the evolution of binary systems. Making the most of expected near-future observations for understanding stellar physics will rely on comparisons with binary population synthesis models. However, the vast majority of simulated binaries never produce DCOs, which mak… ▽ More Gravitational-wave observations of double compact object (DCO) mergers are providing new insights into the physics of massive stars and the evolution of binary systems. Making the most of expected near-future observations for understanding stellar physics will rely on comparisons with binary population synthesis models. However, the vast majority of simulated binaries never produce DCOs, which makes calculating such populations computationally inefficient. We present an importance sampling algorithm, STROOPWAFEL, that improves the computational efficiency of population studies of rare events, by focusing the simulation around regions of the initial parameter space found to produce outputs of interest. We implement the algorithm in the binary population synthesis code COMPAS, and compare the efficiency of our implementation to the standard method of Monte Carlo sampling from the birth probability distributions. STROOPWAFEL finds $\sim$25-200 times more DCO mergers than the standard sampling method with the same simulation size, and so speeds up simulations by up to two orders of magnitude. Finding more DCO mergers automatically maps the parameter space with far higher resolution than when using the traditional sampling. This increase in efficiency also leads to a decrease of a factor $\sim$3-10 in statistical sampling uncertainty for the predictions from the simulations. This is particularly notable for the distribution functions of observable quantities such as the black hole and neutron star chirp mass distribution, including in the tails of the distribution functions where predictions using standard sampling can be dominated by sampling noise. △ Less

Submitted 10 September, 2019; v1 submitted 2 May, 2019; originally announced May 2019.

Comments: Accepted. Data and scripts to reproduce main results is publicly available. The code for the STROOPWAFEL algorithm will be made publicly available. Early inquiries can be addressed to the lead author

On the formation history of Galactic double neutron stars

Authors: Alejandro Vigna-Gómez, Coenraad J. Neijssel, Simon Stevenson, Jim W. Barrett, Krzysztof Belczynski, Stephen Justham, Selma E. de Mink, Bernhard Müller, Philipp Podsiadlowski, Mathieu Renzo, Dorottya Szécsi, Ilya Mandel

Abstract: Double neutron stars (DNSs) have been observed as Galactic radio pulsars, and the recent discovery of gravitational waves from the DNS merger GW170817 adds to the known DNS population. We perform rapid population synthesis of massive binary stars and discuss model predictions, including formation rates, mass distributions, and delay time distributions. We vary assumptions and parameters of physica… ▽ More Double neutron stars (DNSs) have been observed as Galactic radio pulsars, and the recent discovery of gravitational waves from the DNS merger GW170817 adds to the known DNS population. We perform rapid population synthesis of massive binary stars and discuss model predictions, including formation rates, mass distributions, and delay time distributions. We vary assumptions and parameters of physical processes such as mass transfer stability criteria, supernova kick distributions, remnant mass distributions and common-envelope energetics. We compute the likelihood of observing the orbital period-eccentricity distribution of the Galactic DNS population under each of our population synthesis models, allowing us to quantitatively compare the models. We find that mass transfer from a stripped post-helium-burning secondary (case BB) onto a neutron star is most likely dynamically stable. We also find that a natal kick distribution composed of both low (Maxwellian $σ=30\rm~km~s^{-1}$) and high ($σ=265\rm~km~s^{-1}$) components is preferred over a single high-kick component. We find that the observed DNS mass distribution can place strong constraints on model assumptions. △ Less

Submitted 6 September, 2018; v1 submitted 21 May, 2018; originally announced May 2018.

Comments: 23 pages, 10 main figures, 2 tables, 3 figures in appendices; typos added, references added; Formation rates corrected; Table 1 has an extra column

Optimizing spectroscopic follow-up strategies for supernova photometric classification with active learning

Authors: E. E. O. Ishida, R. Beck, S. Gonzalez-Gaitan, R. S. de Souza, A. Krone-Martins, J. W. Barrett, N. Kennamer, R. Vilalta, J. M. Burgess, B. Quint, A. Z. Vitorelli, A. Mahabal, E. Gangler

Abstract: We report a framework for spectroscopic follow-up design for optimizing supernova photometric classification. The strategy accounts for the unavoidable mismatch between spectroscopic and photometric samples, and can be used even in the beginning of a new survey -- without any initial training set. The framework falls under the umbrella of active learning (AL), a class of algorithms that aims to mi… ▽ More We report a framework for spectroscopic follow-up design for optimizing supernova photometric classification. The strategy accounts for the unavoidable mismatch between spectroscopic and photometric samples, and can be used even in the beginning of a new survey -- without any initial training set. The framework falls under the umbrella of active learning (AL), a class of algorithms that aims to minimize labelling costs by identifying a few, carefully chosen, objects which have high potential in improving the classifier predictions. As a proof of concept, we use the simulated data released after the Supernova Photometric Classification Challenge (SNPCC) and a random forest classifier. Our results show that, using only 12\% the number of training objects in the SNPCC spectroscopic sample, this approach is able to double purity results. Moreover, in order to take into account multiple spectroscopic observations in the same night, we propose a semi-supervised batch-mode AL algorithm which selects a set of $N=5$ most informative objects at each night. In comparison with the initial state using the traditional approach, our method achieves 2.3 times higher purity and comparable figure of merit results after only 180 days of observation, or 800 queries (73% of the SNPCC spectroscopic sample size). Such results were obtained using the same amount of spectroscopic time necessary to observe the original SNPCC spectroscopic sample, showing that this type of strategy is feasible with current available spectroscopic resources. The code used in this work is available in the COINtoolbox: https://github.com/COINtoolbox/ActSNClass . △ Less

Submitted 3 January, 2019; v1 submitted 10 April, 2018; originally announced April 2018.

Comments: 18 pages, 15 figures - replace to match journal version

Journal ref: MNRAS, Volume 483, Issue 1, 11 February 2019, Pages 2-18

Aldebaran b's temperate past uncovered in planet search data

Authors: Will M. Farr, Benjamin J. S. Pope, Guy R. Davies, Thomas S. H. North, Timothy R. White, Jim W. Barrett, Andrea Miglio, Mikkel N. Lund, Victoria Antoci, Mads Fredslund Andersen, Frank Grundahl, Daniel Huber

Abstract: The nearby red giant Aldebaran is known to host a gas giant planetary companion from decades of ground-based spectroscopic radial velocity measurements. Using Gaussian Process-based Continuous Auto-Regressive Moving Average (CARMA) models, we show that these historic data also contain evidence of acoustic oscillations in the star itself, and verify this result with further dedicated ground-based s… ▽ More The nearby red giant Aldebaran is known to host a gas giant planetary companion from decades of ground-based spectroscopic radial velocity measurements. Using Gaussian Process-based Continuous Auto-Regressive Moving Average (CARMA) models, we show that these historic data also contain evidence of acoustic oscillations in the star itself, and verify this result with further dedicated ground-based spectroscopy and space-based photometry with the Kepler Space Telescope. From the frequency of these oscillations we determine the mass of Aldebaran to be $1.16 \pm 0.07 \, M_\odot$, and note that this implies its planet will have been subject to insolation comparable to the Earth for some of the star's main sequence lifetime. Our approach to sparse, irregularly sampled time series astronomical observations has the potential to unlock asteroseismic measurements for thousands of stars in archival data, and push to lower-mass planets around red giant stars. △ Less

Submitted 27 February, 2018; originally announced February 2018.

Comments: 24 pages, 7 figures (including appendices); submitted to ApJL; paper text, figures, data, and code at https://github.com/farr/Aldebaran

Accuracy of inference on the physics of binary evolution from gravitational-wave observations

Authors: Jim W. Barrett, Sebastian M. Gaebel, Coenraad J. Neijssel, Alejandro Vigna-Gómez, Simon Stevenson, Christopher P. L. Berry, Will M. Farr, Ilya Mandel

Abstract: The properties of the population of merging binary black holes encode some of the uncertain physics of the evolution of massive stars in binaries. The binary black hole merger rate and chirp mass distribution are being measured by ground-based gravitational-wave detectors. We consider isolated binary evolution and explore how accurately the physical model can be constrained with such observations… ▽ More The properties of the population of merging binary black holes encode some of the uncertain physics of the evolution of massive stars in binaries. The binary black hole merger rate and chirp mass distribution are being measured by ground-based gravitational-wave detectors. We consider isolated binary evolution and explore how accurately the physical model can be constrained with such observations by applying the Fisher information matrix to the merging black hole population simulated with the rapid binary population synthesis code COMPAS. We investigate variations in four COMPAS parameters: common envelope efficiency, kick velocity dispersion, and mass loss rates during the luminous blue variable and Wolf--Rayet stellar evolutionary phases. We find that 1000 observations would constrain these model parameters to a fractional accuracy of a few percent. Given the empirically determined binary black hole merger rate, we can expect gravitational-wave observations alone to place strong constraints on the physics of stellar and binary evolution within a few years. △ Less

Submitted 16 May, 2018; v1 submitted 16 November, 2017; originally announced November 2017.

Comments: 12 pages, 9 figures; version accepted by Monthly Notices of the Royal Astronomical Society

Journal ref: MNRAS; 477(4): 4685-4695; 2018

Exploring the Parameter Space of Compact Binary Population Synthesis

Authors: Jim W. Barrett, Ilya Mandel, Coenraad J. Neijssel, Simon Stevenson, Alejandro Vigna-Gomez

Abstract: As we enter the era of gravitational wave astronomy, we are beginning to collect observations which will enable us to explore aspects of astrophysics of massive stellar binaries which were previously beyond reach. In this paper we describe COMPAS (Compact Object Mergers: Population Astrophysics and Statistics), a new platform to allow us to deepen our understanding of isolated binary evolution and… ▽ More As we enter the era of gravitational wave astronomy, we are beginning to collect observations which will enable us to explore aspects of astrophysics of massive stellar binaries which were previously beyond reach. In this paper we describe COMPAS (Compact Object Mergers: Population Astrophysics and Statistics), a new platform to allow us to deepen our understanding of isolated binary evolution and the formation of gravitational-wave sources. We describe the computational challenges associated with their exploration, and present preliminary results on overcoming them using Gaussian process regression as a simulation emulation technique. △ Less

Submitted 12 April, 2017; originally announced April 2017.

Comments: Accepted for Proceedings of IAU Symposium 325

Formation of the first three gravitational-wave observations through isolated binary evolution

Authors: Simon Stevenson, Alejandro Vigna-Gómez, Ilya Mandel, Jim W. Barrett, Coenraad J. Neijssel, David Perkins, Selma E. de Mink

Abstract: During its first 4 months of taking data, Advanced LIGO has detected gravitational waves from two binary black hole mergers, GW150914 and GW151226, along with the statistically less significant binary black hole merger candidate LVT151012. We use our rapid binary population synthesis code COMPAS to show that all three events can be explained by a single evolutionary channel -- classical isolated b… ▽ More During its first 4 months of taking data, Advanced LIGO has detected gravitational waves from two binary black hole mergers, GW150914 and GW151226, along with the statistically less significant binary black hole merger candidate LVT151012. We use our rapid binary population synthesis code COMPAS to show that all three events can be explained by a single evolutionary channel -- classical isolated binary evolution via mass transfer including a common envelope phase. We show all three events could have formed in low-metallicity environments (Z = 0.001) from progenitor binaries with typical total masses $\gtrsim 160 M_\odot$, $\gtrsim 60 M_\odot$ and $\gtrsim 90 M_\odot$, for GW150914, GW151226, and LVT151012, respectively. △ Less

Submitted 10 April, 2017; v1 submitted 5 April, 2017; originally announced April 2017.

Comments: Published in Nature Communications

Journal ref: Nature Communications 8, Article number: 14906 (2017)

Dual frequency 230/690 GHz interferometry at the Submillimeter Array

Authors: Todd R. Hunter, John W. Barrett, Raymond Blundell, Robert D. Christensen, Robert S. Kimberk, Steven P. Leiker, Daniel P. Marrone, Scott N. Paine, D. Cosmo Papa, Nimesh Patel, Patricia Riddle, Michael J. Smith, T. K. Sridharan, C. Y. Edward Tong, Ken H. Young, Jun-Hui Zhao

Abstract: The Submillimeter Array (SMA), a collaboration between the Smithsonian Astrophysical Observatory and the Academica Sinica Institute for Astronomy and Astrophysics of Taiwan, is an eight-element radio-interferometer designed to operate throughout the major atmospheric windows from about 180 to 900 GHz. In an effort to mitigate the effects of atmospheric instabilities which limit the phase coheren… ▽ More The Submillimeter Array (SMA), a collaboration between the Smithsonian Astrophysical Observatory and the Academica Sinica Institute for Astronomy and Astrophysics of Taiwan, is an eight-element radio-interferometer designed to operate throughout the major atmospheric windows from about 180 to 900 GHz. In an effort to mitigate the effects of atmospheric instabilities which limit the phase coherence of the array especially in the higher frequency bands, the array was designed to allow simultaneous operation of a low frequency receiver (<350 GHz) with a high frequency receiver (>330 GHz). The overlap region of 330-350 GHz was included to facilitate dual polarization measurements in the frequency range considered to offer the highest sensitivity for continuum observations with the array. So far, the array is equipped with working SIS receivers covering the frequency ranges 176-256 GHz, 260-350 GHz, and 600-700 GHz, and single frequency operation has been routine in the lower two frequency bands for the past year. More recently, with the completion of IF hardware required to make full use of the SMA cross-correlator, dual receiver operation became possible. We have since made a number of Galactic and extra-galactic astronomical observations in dual-band mode with the hopes of using the 230 GHz receiver as a phase reference to enable improved interferometry in the 650 GHz band. We will present the current antenna and receiver performance, some of the first interferometric images in the 650 GHz receiver band, and our initial attempts at phase transfer. △ Less

Submitted 15 September, 2005; originally announced September 2005.

Comments: 6 pages containing text and 11 embedded figures, written in the IEEE Transactions style Latex template, based on an Oral Presentation at the 16th International Symposium on Space Terahertz Technology held at Chalmers University of Technology in Goteborg, Sweden on May 2-4, 2005

The Molecular Interstellar Medium in Ultraluminous Infrared Galaxies

Authors: P. M. Solomon, D. Downes, S. J. E. Radford, J. W. Barrett

Abstract: We present CO observations of a large sample of ultraluminous IR galaxies out to z = 0.3. Most of the galaxies are interacting, but not completed mergers. All but one have high CO(1-0) luminosities, log(Lco [K-km/s-pc^2]) = 9.92 +/- 0.12. The dispersion in Lco is only 30%, less than that in the FIR luminosity. The integrated CO intensity correlates Strongly with the 100 micron flux density, as e… ▽ More We present CO observations of a large sample of ultraluminous IR galaxies out to z = 0.3. Most of the galaxies are interacting, but not completed mergers. All but one have high CO(1-0) luminosities, log(Lco [K-km/s-pc^2]) = 9.92 +/- 0.12. The dispersion in Lco is only 30%, less than that in the FIR luminosity. The integrated CO intensity correlates Strongly with the 100 micron flux density, as expected for a black body model in which the mid and far IR radiation are optically thick. We use this model to derive sizes of the FIR and CO emitting regions and the enclosed dynamical masses. Both the IR and CO emission originate in regions a few hundred parsecs in radius. The median value of Lfir/Lco = 160 Lsun/(K-km/s-pc^2), within a factor of two of the black body limit for the observed FIR temperatures. The entire ISM is a scaled up version of a normal galactic disk with densities a factor of 100 higher, making even the intercloud medium a molecular region. Using three different techniques of H2 mass estimation, we conclude that the ratio of gas mass to Lco is about a factor of four lower than for Galactic molecular clouds, but that the gas mass is a large fraction of the dynamical mass. Our analysis of CO emission reduces the H2 mass from previous estimates of 2-5e10 Msun to 0.4-1.5e10 Msun, which is in the range found for molecular gas rich spiral galaxies. A collision involving a molecular gas rich spiral could lead to an ultraluminous galaxy powered by central starbursts triggered by the compression of infalling preexisting GMC's. △ Less

Submitted 21 October, 1996; originally announced October 1996.

Comments: 34 pages LaTeX with aasms.sty, 14 Postscript figures, submitted to ApJ Higher quality versions of Figs 2a-f and 7a-c available by anonymous FTP from ftp://sbast1.ess.sunysb.edu/solomon/ .

Journal ref: Astrophys.J.478:144,1997