Showing 1–50 of 224 results for author: Brown, D A

The VISTA Variables in the Vía Láctea eXtended (VVVX) ESO public survey: Completion of the observations and legacy

Authors: R. K. Saito, M. Hempel, J. Alonso-García, P. W. Lucas, D. Minniti, S. Alonso, L. Baravalle, J. Borissova, C. Caceres, A. N. Chené, N. J. G. Cross, F. Duplancic, E. R. Garro, M. Gómez, V. D. Ivanov, R. Kurtev, A. Luna, D. Majaess, M. G. Navarro, J. B. Pullen, M. Rejkuba, J. L. Sanders, L. C. Smith, P. H. C. Albino, M. V. Alonso , et al. (121 additional authors not shown)

Abstract: The ESO public survey VISTA Variables in the Vía Láctea (VVV) surveyed the inner Galactic bulge and the adjacent southern Galactic disk from $2009-2015$. Upon its conclusion, the complementary VVV eXtended (VVVX) survey has expanded both the temporal as well as spatial coverage of the original VVV area, widening it from $562$ to $1700$ sq. deg., as well as providing additional epochs in… ▽ More The ESO public survey VISTA Variables in the Vía Láctea (VVV) surveyed the inner Galactic bulge and the adjacent southern Galactic disk from $2009-2015$. Upon its conclusion, the complementary VVV eXtended (VVVX) survey has expanded both the temporal as well as spatial coverage of the original VVV area, widening it from $562$ to $1700$ sq. deg., as well as providing additional epochs in $JHK_{\rm s}$ filters from $2016-2023$. With the completion of VVVX observations during the first semester of 2023, we present here the observing strategy, a description of data quality and access, and the legacy of VVVX. VVVX took $\sim 2000$ hours, covering about 4% of the sky in the bulge and southern disk. VVVX covered most of the gaps left between the VVV and the VISTA Hemisphere Survey (VHS) areas and extended the VVV time baseline in the obscured regions affected by high extinction and hence hidden from optical observations. VVVX provides a deep $JHK_{\rm s}$ catalogue of $\gtrsim 1.5\times10^9$ point sources, as well as a $K_{\rm s}$ band catalogue of $\sim 10^7$ variable sources. Within the existing VVV area, we produced a $5D$ map of the surveyed region by combining positions, distances, and proper motions of well-understood distance indicators such as red clump stars, RR Lyrae, and Cepheid variables. In March 2023 we successfully finished the VVVX survey observations that started in 2016, an accomplishment for ESO Paranal Observatory upon 4200 hours of observations for VVV+VVVX. The VVV+VVVX catalogues complement those from the Gaia mission at low Galactic latitudes and provide spectroscopic targets for the forthcoming ESO high-multiplex spectrographs MOONS and 4MOST. △ Less

Submitted 24 June, 2024; originally announced June 2024.

Comments: 17 pages, 11 figures (+ appendix). Accepted for publication in Astronomy and Astrophysics in section 14: Catalogs and data

Towards accelerated nuclear-physics parameter estimation from binary neutron star mergers: Emulators for the Tolman-Oppenheimer-Volkoff equations

Authors: Brendan T. Reed, Rahul Somasundaram, Soumi De, Cassandra L. Armstrong, Pablo Giuliani, Collin Capano, Duncan A. Brown, Ingo Tews

Abstract: Gravitational-wave observations of binary neutron-star (BNS) mergers have the potential to revolutionize our understanding of the nuclear equation of state (EOS) and the fundamental interactions that determine its properties. However, Bayesian parameter estimation frameworks do not typically sample over microscopic nuclear-physics parameters that determine the EOS. One of the major hurdles in doin… ▽ More Gravitational-wave observations of binary neutron-star (BNS) mergers have the potential to revolutionize our understanding of the nuclear equation of state (EOS) and the fundamental interactions that determine its properties. However, Bayesian parameter estimation frameworks do not typically sample over microscopic nuclear-physics parameters that determine the EOS. One of the major hurdles in doing so is the computational cost involved in solving the neutron-star structure equations, known as the Tolman-Oppenheimer-Volkoff (TOV) equations. In this paper, we explore approaches to emulating solutions for the TOV equations: Multilayer Perceptrons (MLP), Gaussian Processes (GP), and a data-driven variant of the reduced basis method (RBM). We implement these emulators for three different parameterizations of the nuclear EOS, each with a different degree of complexity represented by the number of model parameters. We find that our MLP-based emulators are generally more accurate than the other two algorithms whereas the RBM results in the largest speedup with respect to the full, high-fidelity TOV solver. We employ these emulators for a simple parameter inference using a potentially loud BNS observation, and show that the posteriors predicted by our emulators are in excellent agreement with those obtained from the full TOV solver. △ Less

Submitted 30 May, 2024; originally announced May 2024.

Comments: 13 pages and 9 figures. Comments Welcome

Report number: LA-UR-24-25009

Measuring Neutron Star Radius with second and third generation Gravitational Wave Detector Networks

Authors: Ananya Bandopadhyay, Keisi Kacanja, Rahul Somasundaram, Alexander H. Nitz, Duncan A. Brown

Abstract: The next generation of ground-based interferometric gravitational wave detectors will observe mergers of black holes and neutron stars throughout cosmic time. A large number of the binary neutron star merger events will be observed with extreme high fidelity, and will provide stringent constraints on the equation of state of nuclear matter. In this paper, we investigate the systematic improvement… ▽ More The next generation of ground-based interferometric gravitational wave detectors will observe mergers of black holes and neutron stars throughout cosmic time. A large number of the binary neutron star merger events will be observed with extreme high fidelity, and will provide stringent constraints on the equation of state of nuclear matter. In this paper, we investigate the systematic improvement in the measurability of the equation of state with increase in detector sensitivity by combining constraints obtained on the radius of a $1.4 \, \mathrm{M}_{\odot}$ neutron star from a simulated source population. Since the measurability of the equation of state depends on its stiffness, we consider a range of realistic equations of state that span the current observational constraints. We show that a single 40km Cosmic Explorer detector can pin down the neutron star radius for a soft, medium and stiff equation of state to an accuracy of 10m within a decade, whereas the current generation of ground-based detectors like the Advanced LIGO-Virgo network would take $\mathcal{O}(10^5)$ years to do so for a soft equation of state. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Comments: 14 pages, 3 figures, 1 table, supplemental materials at https://github.com/sugwg/bns-eos-nggw

Report number: LA-UR-24-21031

A Fully Bayesian Approach for Comprehensive Mapping of Magnitude and Phase Brain Activation in Complex-Valued fMRI Data

Authors: Zhengxin Wang, Daniel B. Rowe, Xinyi Li, D. Andrew Brown

Abstract: Functional magnetic resonance imaging (fMRI) plays a crucial role in neuroimaging, enabling the exploration of brain activity through complex-valued signals. These signals, composed of magnitude and phase, offer a rich source of information for understanding brain functions. Traditional fMRI analyses have largely focused on magnitude information, often overlooking the potential insights offered by… ▽ More Functional magnetic resonance imaging (fMRI) plays a crucial role in neuroimaging, enabling the exploration of brain activity through complex-valued signals. These signals, composed of magnitude and phase, offer a rich source of information for understanding brain functions. Traditional fMRI analyses have largely focused on magnitude information, often overlooking the potential insights offered by phase data. In this paper, we propose a novel fully Bayesian model designed for analyzing single-subject complex-valued fMRI (cv-fMRI) data. Our model, which we refer to as the CV-M&P model, is distinctive in its comprehensive utilization of both magnitude and phase information in fMRI signals, allowing for independent prediction of different types of activation maps. We incorporate Gaussian Markov random fields (GMRFs) to capture spatial correlations within the data, and employ image partitioning and parallel computation to enhance computational efficiency. Our model is rigorously tested through simulation studies, and then applied to a real dataset from a unilateral finger-tapping experiment. The results demonstrate the model's effectiveness in accurately identifying brain regions activated in response to specific tasks, distinguishing between magnitude and phase activation. △ Less

Submitted 11 January, 2024; originally announced January 2024.

Efficient Fully Bayesian Approach to Brain Activity Mapping with Complex-Valued fMRI Data

Authors: Zhengxin Wang, Daniel B. Rowe, Xinyi Li, D. Andrew Brown

Abstract: Functional magnetic resonance imaging (fMRI) enables indirect detection of brain activity changes via the blood-oxygen-level-dependent (BOLD) signal. Conventional analysis methods mainly rely on the real-valued magnitude of these signals. In contrast, research suggests that analyzing both real and imaginary components of the complex-valued fMRI (cv-fMRI) signal provides a more holistic approach th… ▽ More Functional magnetic resonance imaging (fMRI) enables indirect detection of brain activity changes via the blood-oxygen-level-dependent (BOLD) signal. Conventional analysis methods mainly rely on the real-valued magnitude of these signals. In contrast, research suggests that analyzing both real and imaginary components of the complex-valued fMRI (cv-fMRI) signal provides a more holistic approach that can increase power to detect neuronal activation. We propose a fully Bayesian model for brain activity mapping with cv-fMRI data. Our model accommodates temporal and spatial dynamics. Additionally, we propose a computationally efficient sampling algorithm, which enhances processing speed through image partitioning. Our approach is shown to be computationally efficient via image partitioning and parallel computation while being competitive with state-of-the-art methods. We support these claims with both simulated numerical studies and an application to real cv-fMRI data obtained from a finger-tapping experiment. △ Less

Submitted 27 October, 2023; originally announced October 2023.

Cosmic Explorer: A Submission to the NSF MPSAC ngGW Subcommittee

Authors: Matthew Evans, Alessandra Corsi, Chaitanya Afle, Alena Ananyeva, K. G. Arun, Stefan Ballmer, Ananya Bandopadhyay, Lisa Barsotti, Masha Baryakhtar, Edo Berger, Emanuele Berti, Sylvia Biscoveanu, Ssohrab Borhanian, Floor Broekgaarden, Duncan A. Brown, Craig Cahillane, Lorna Campbell, Hsin-Yu Chen, Kathryne J. Daniel, Arnab Dhani, Jennifer C. Driggers, Anamaria Effler, Robert Eisenstein, Stephen Fairhurst, Jon Feicht , et al. (51 additional authors not shown)

Abstract: Gravitational-wave astronomy has revolutionized humanity's view of the universe, a revolution driven by observations that no other field can make. This white paper describes an observatory that builds on decades of investment by the National Science Foundation and that will drive discovery for decades to come: Cosmic Explorer. Major discoveries in astronomy are driven by three related improvements… ▽ More Gravitational-wave astronomy has revolutionized humanity's view of the universe, a revolution driven by observations that no other field can make. This white paper describes an observatory that builds on decades of investment by the National Science Foundation and that will drive discovery for decades to come: Cosmic Explorer. Major discoveries in astronomy are driven by three related improvements: better sensitivity, higher precision, and opening new observational windows. Cosmic Explorer promises all three and will deliver an order-of-magnitude greater sensitivity than LIGO. Cosmic Explorer will push the gravitational-wave frontier to almost the edge of the observable universe using technologies that have been proven by LIGO during its development. With the unprecedented sensitivity that only a new facility can deliver, Cosmic Explorer will make discoveries that cannot yet be anticipated, especially since gravitational waves are both synergistic with electromagnetic observations and can reach into regions of the universe that electromagnetic observations cannot explore. With Cosmic Explorer, scientists can use the universe as a laboratory to test the laws of physics and study the nature of matter. Cosmic Explorer allows the United States to continue its leading role in gravitational-wave science and the international network of next-generation observatories. With its extraordinary discovery potential, Cosmic Explorer will deliver revolutionary observations across astronomy, physics, and cosmology including: Black Holes and Neutron Stars Throughout Cosmic Time, Multi-Messenger Astrophysics and Dynamics of Dense Matter, New Probes of Extreme Astrophysics, Fundamental Physics and Precision Cosmology, Dark Matter and the Early Universe. △ Less

Submitted 23 June, 2023; originally announced June 2023.

Measuring the properties of $f-$mode oscillations of a protoneutron star by third generation gravitational-wave detectors

Authors: Chaitanya Afle, Suman Kumar Kundu, Jenna Cammerino, Eric R Coughlin, Duncan A. Brown, David Vartanyan, Adam Burrows

Abstract: Core-collapse supernovae are among the astrophysical sources of gravitational waves that could be detected by third-generation gravitational-wave detectors. Here, we analyze the gravitational-wave strain signals from two- and three-dimensional simulations of core-collapse supernovae generated using the code F{\sc{ornax}}. A subset of the two-dimensional simulations has non-zero core rotation at th… ▽ More Core-collapse supernovae are among the astrophysical sources of gravitational waves that could be detected by third-generation gravitational-wave detectors. Here, we analyze the gravitational-wave strain signals from two- and three-dimensional simulations of core-collapse supernovae generated using the code F{\sc{ornax}}. A subset of the two-dimensional simulations has non-zero core rotation at the core bounce. A dominant source of time changing quadrupole moment is the $l=2$ fundamental mode ($f-$ mode) oscillation of the proto-neutron star. From the time-frequency spectrogram of the gravitational-wave strain we see that, starting $\sim 400$ ms after the core bounce, most of the power lies within a narrow track that represents the frequency evolution of the $f-$mode oscillations. The $f-$mode frequencies obtained from linear perturbation analysis of the angle-averaged profile of the protoneutron star corroborate what we observe in the spectrograms of the gravitational-wave signal. We explore the measurability of the $f-$mode frequency evolution of protoneutron star for a supernova signal observed in the third-generation gravitational-wave detectors. Measurement of the frequency evolution can reveal information about the masses, radii, and densities of the proto-neutron stars. We find that if the third generation detectors observe a supernova within 10 kpc, we can measure these frequencies to within $\sim$90\% accuracy. We can also measure the energy emitted in the fundamental $f-$mode using the spectrogram data of the strain signal. We find that the energy in the $f-$mode can be measured to within 20\% error for signals observed by Cosmic Explorer using simulations with successful explosion, assuming source distances within 10 kpc. △ Less

Submitted 9 April, 2023; originally announced April 2023.

Comments: 17 pages, 11 figures, 2 tables

Reproducing the results for NICER observation of PSR J0030+0451

Authors: Chaitanya Afle, Patrick R. Miles, Silvina Caino-Lores, Collin D. Capano, Ingo Tews, Karan Vahi, Ewa Deelman, Michela Taufer, Duncan A. Brown

Abstract: NASA's Neutron Star Interior Composition Explorer (NICER) observed X-ray emission from the pulsar PSR J0030+0451 in 2018. Riley et al. reported Bayesian parameter measurements of the mass and the star's radius using pulse-profile modeling of the X-ray data. This paper reproduces their result using the open-source software X-PSI and publicly available data within expected statistical errors. We not… ▽ More NASA's Neutron Star Interior Composition Explorer (NICER) observed X-ray emission from the pulsar PSR J0030+0451 in 2018. Riley et al. reported Bayesian parameter measurements of the mass and the star's radius using pulse-profile modeling of the X-ray data. This paper reproduces their result using the open-source software X-PSI and publicly available data within expected statistical errors. We note the challenges we faced in reproducing the results and demonstrate that the analysis can be reproduced and reused in future works by changing the prior distribution for the radius and the sampler configuration. We find no significant change in the measurement of the mass and radius, demonstrating that the original result is robust to these changes. Finally, we provide a containerized working environment that facilitates third-party reproduction of the measurements of mass and radius of PSR J0030+0451 using the NICER observations. △ Less

Submitted 31 January, 2024; v1 submitted 3 April, 2023; originally announced April 2023.

Comments: 13 pages, 4 figures, 2 tables. Final version accepted for publication in Computing in Science & Engineering

Report number: LA-UR-22-29359

The Present and Future of QCD

Authors: P. Achenbach, D. Adhikari, A. Afanasev, F. Afzal, C. A. Aidala, A. Al-bataineh, D. K. Almaalol, M. Amaryan, D. Androić, W. R. Armstrong, M. Arratia, J. Arrington, A. Asaturyan, E. C. Aschenauer, H. Atac, H. Avakian, T. Averett, C. Ayerbe Gayoso, X. Bai, K. N. Barish, N. Barnea, G. Basar, M. Battaglieri, A. A. Baty, I. Bautista , et al. (378 additional authors not shown)

Abstract: This White Paper presents the community inputs and scientific conclusions from the Hot and Cold QCD Town Meeting that took place September 23-25, 2022 at MIT, as part of the Nuclear Science Advisory Committee (NSAC) 2023 Long Range Planning process. A total of 424 physicists registered for the meeting. The meeting highlighted progress in Quantum Chromodynamics (QCD) nuclear physics since the 2015… ▽ More This White Paper presents the community inputs and scientific conclusions from the Hot and Cold QCD Town Meeting that took place September 23-25, 2022 at MIT, as part of the Nuclear Science Advisory Committee (NSAC) 2023 Long Range Planning process. A total of 424 physicists registered for the meeting. The meeting highlighted progress in Quantum Chromodynamics (QCD) nuclear physics since the 2015 LRP (LRP15) and identified key questions and plausible paths to obtaining answers to those questions, defining priorities for our research over the coming decade. In defining the priority of outstanding physics opportunities for the future, both prospects for the short (~ 5 years) and longer term (5-10 years and beyond) are identified together with the facilities, personnel and other resources needed to maximize the discovery potential and maintain United States leadership in QCD physics worldwide. This White Paper is organized as follows: In the Executive Summary, we detail the Recommendations and Initiatives that were presented and discussed at the Town Meeting, and their supporting rationales. Section 2 highlights major progress and accomplishments of the past seven years. It is followed, in Section 3, by an overview of the physics opportunities for the immediate future, and in relation with the next QCD frontier: the EIC. Section 4 provides an overview of the physics motivations and goals associated with the EIC. Section 5 is devoted to the workforce development and support of diversity, equity and inclusion. This is followed by a dedicated section on computing in Section 6. Section 7 describes the national need for nuclear data science and the relevance to QCD research. △ Less

Submitted 4 March, 2023; originally announced March 2023.

Comments: QCD Town Meeting White Paper, as submitted to 2023 NSAC LRP committee on Feb. 28, 2023

Journal ref: Nucl.Phys.A 1047 (2024) 122874

Dense Nuclear Matter Equation of State from Heavy-Ion Collisions

Authors: Agnieszka Sorensen, Kshitij Agarwal, Kyle W. Brown, Zbigniew Chajęcki, Paweł Danielewicz, Christian Drischler, Stefano Gandolfi, Jeremy W. Holt, Matthias Kaminski, Che-Ming Ko, Rohit Kumar, Bao-An Li, William G. Lynch, Alan B. McIntosh, William G. Newton, Scott Pratt, Oleh Savchuk, Maria Stefaniak, Ingo Tews, ManYee Betty Tsang, Ramona Vogt, Hermann Wolter, Hanna Zbroszczyk, Navid Abbasi, Jörg Aichelin , et al. (111 additional authors not shown)

Abstract: The nuclear equation of state (EOS) is at the center of numerous theoretical and experimental efforts in nuclear physics. With advances in microscopic theories for nuclear interactions, the availability of experiments probing nuclear matter under conditions not reached before, endeavors to develop sophisticated and reliable transport simulations to interpret these experiments, and the advent of mu… ▽ More The nuclear equation of state (EOS) is at the center of numerous theoretical and experimental efforts in nuclear physics. With advances in microscopic theories for nuclear interactions, the availability of experiments probing nuclear matter under conditions not reached before, endeavors to develop sophisticated and reliable transport simulations to interpret these experiments, and the advent of multi-messenger astronomy, the next decade will bring new opportunities for determining the nuclear matter EOS, elucidating its dependence on density, temperature, and isospin asymmetry. Among controlled terrestrial experiments, collisions of heavy nuclei at intermediate beam energies (from a few tens of MeV/nucleon to about 25 GeV/nucleon in the fixed-target frame) probe the widest ranges of baryon density and temperature, enabling studies of nuclear matter from a few tenths to about 5 times the nuclear saturation density and for temperatures from a few to well above a hundred MeV, respectively. Collisions of neutron-rich isotopes further bring the opportunity to probe effects due to the isospin asymmetry. However, capitalizing on the enormous scientific effort aimed at uncovering the dense nuclear matter EOS, both at RHIC and at FRIB as well as at other international facilities, depends on the continued development of state-of-the-art hadronic transport simulations. This white paper highlights the essential role that heavy-ion collision experiments and hadronic transport simulations play in understanding strong interactions in dense nuclear matter, with an emphasis on how these efforts can be used together with microscopic approaches and neutron star studies to uncover the nuclear EOS. △ Less

Submitted 25 January, 2024; v1 submitted 30 January, 2023; originally announced January 2023.

Comments: White paper prepared for the 2023 Long Range Plan. v3: Updated version as published in Progress in Particle and Nuclear Physics. Note: the published version does not include the executive summary; in the updated arXiv version, the executive summary is included as an appendix. v4: Corrected list of authors

Report number: INT-PUB-23-001, LA-UR-23-20514, LLNL-TR-844629

Journal ref: Prog.Part.Nucl.Phys. 134 (2024) 104080

Detectability of Sub-Solar Mass Neutron Stars Through a Template Bank Search

Authors: Ananya Bandopadhyay, Brendan Reed, Surendra Padamata, Erick Leon, C. J. Horowitz, Duncan A. Brown, David Radice, F. J. Fattoyev, J. Piekarewicz

Abstract: We study the detectability of gravitational-wave signals from sub-solar mass binary neutron star systems by the current generation of ground-based gravitational-wave detectors. We find that finite size effects from large tidal deformabilities of the neutron stars and lower merger frequencies can significantly impact the sensitivity of the detectors to these sources. By simulating a matched-filter… ▽ More We study the detectability of gravitational-wave signals from sub-solar mass binary neutron star systems by the current generation of ground-based gravitational-wave detectors. We find that finite size effects from large tidal deformabilities of the neutron stars and lower merger frequencies can significantly impact the sensitivity of the detectors to these sources. By simulating a matched-filter based search using injected binary neutron star signals with tidal deformabilities derived from physically motivated equations of state, we calculate the reduction in sensitivity of the detectors. We conclude that the loss in sensitive volume can be as high as $78.4 \%$ for an equal mass binary system of chirp mass $0.17 \, \textrm{M}_{\odot}$, in a search conducted using binary black hole template banks. We use this loss in sensitive volume, in combination with the results from the search for sub-solar mass binaries conducted on data collected by the LIGO-Virgo observatories during their first three observing runs, to obtain a conservative upper limit on the merger rate of sub-solar mass binary neutron stars. Since the discovery of a low-mass neutron star would provide new insight into formation mechanisms of neutron stars and further constrain the equation of state of dense nuclear matter, our result merits a dedicated search for sub-solar mass binary neutron star signals. △ Less

Submitted 9 May, 2023; v1 submitted 7 December, 2022; originally announced December 2022.

Comments: 12 pages, 7 figures, supplemental materials at https://github.com/sugwg/sub-solar-ns-detectability

Journal ref: Phys. Rev. D 107, 103012 (2023)

Prospects for a precise equation of state measurement from Advanced LIGO and Cosmic Explorer

Authors: Daniel Finstad, Laurel V. White, Duncan A. Brown

Abstract: Gravitational-wave observations of neutron star mergers can probe the nuclear equation of state by measuring the imprint of the neutron star's tidal deformability on the signal. We investigate the ability of future gravitational-wave observations to produce a precise measurement of the equation of state from binary neutron star inspirals. Since measurability of the tidal effect depends on the equa… ▽ More Gravitational-wave observations of neutron star mergers can probe the nuclear equation of state by measuring the imprint of the neutron star's tidal deformability on the signal. We investigate the ability of future gravitational-wave observations to produce a precise measurement of the equation of state from binary neutron star inspirals. Since measurability of the tidal effect depends on the equation of state, we explore several equations of state that span current observational constraints. We generate a population of binary neutron stars as seen by a simulated Advanced LIGO-Virgo network, as well as by a planned Cosmic Explorer observatory. We perform Bayesian inference to measure the parameters of each signal, and we combine measurements across each population to determine $R_{1.4}$, the radius of a $1.4M_{\odot}$ neutron star. We find that with 321 signals the LIGO-Virgo network is able to measure $R_{1.4}$ to better than 2% precision for all equations of state we consider, however we find that achieving this precision could take decades of observation, depending on the equation of state and the merger rate. On the other hand we find that with one year of observation, Cosmic Explorer will measure $R_{1.4}$ to better than 0.6% precision. In both cases we find that systematic biases, such as from an incorrect mass prior, can significantly impact measurement accuracy and efforts will be required to mitigate these effects. △ Less

Submitted 2 November, 2022; originally announced November 2022.

Comments: 11 pages, 6 figures

A novel Machine-Learning method for spin classification of neutron resonances

Authors: G. P. A. Nobre, D. A. Brown, S. J. Hollick, S. Scoville, P. Rodríguez

Abstract: The performance of nuclear reactors and other nuclear systems depends on a precise understanding of the neutron interaction cross sections for materials used in these systems. These cross sections exhibit resonant structure whose shape is determined in part by the angular momentum quantum numbers of the resonances. The correct assignment of the quantum numbers of neutron resonances is, therefore,… ▽ More The performance of nuclear reactors and other nuclear systems depends on a precise understanding of the neutron interaction cross sections for materials used in these systems. These cross sections exhibit resonant structure whose shape is determined in part by the angular momentum quantum numbers of the resonances. The correct assignment of the quantum numbers of neutron resonances is, therefore, paramount. In this project, we apply machine learning to automate the quantum number assignments using only the resonances' energies and widths and not relying on detailed transmission or capture measurements. The classifier used for quantum number assignment is trained using stochastically generated resonance sequences whose distributions mimic those of real data. We explore the use of several physics-motivated features for training our classifier. These features amount to out-of-distribution tests of a given resonance's widths and resonance-pair spacings. We pay special attention to situations where either capture widths cannot be trusted for classification purposes or where there is insufficient information to classify resonances by the total spin $J$. We demonstrate the efficacy of our classification approach using simulated and actual $^{52}$Cr resonance data. △ Less

Submitted 16 February, 2023; v1 submitted 28 September, 2022; originally announced September 2022.

Comments: 22 pages, 5 tables, 13 figures. Version resubmitted to PRC after comments from referee (including bug fix) and accepted for publication

Journal ref: Phys. Rev. C 107, 034612, 2023

Advancing the Landscape of Multimessenger Science in the Next Decade

Authors: Kristi Engel, Tiffany Lewis, Marco Stein Muzio, Tonia M. Venters, Markus Ahlers, Andrea Albert, Alice Allen, Hugo Alberto Ayala Solares, Samalka Anandagoda, Thomas Andersen, Sarah Antier, David Alvarez-Castillo, Olaf Bar, Dmitri Beznosko, Łukasz Bibrzyck, Adam Brazier, Chad Brisbois, Robert Brose, Duncan A. Brown, Mattia Bulla, J. Michael Burgess, Eric Burns, Cecilia Chirenti, Stefano Ciprini, Roger Clay , et al. (69 additional authors not shown)

Abstract: The last decade has brought about a profound transformation in multimessenger science. Ten years ago, facilities had been built or were under construction that would eventually discover the nature of objects in our universe could be detected through multiple messengers. Nonetheless, multimessenger science was hardly more than a dream. The rewards for our foresight were finally realized through Ice… ▽ More The last decade has brought about a profound transformation in multimessenger science. Ten years ago, facilities had been built or were under construction that would eventually discover the nature of objects in our universe could be detected through multiple messengers. Nonetheless, multimessenger science was hardly more than a dream. The rewards for our foresight were finally realized through IceCube's discovery of the diffuse astrophysical neutrino flux, the first observation of gravitational waves by LIGO, and the first joint detections in gravitational waves and photons and in neutrinos and photons. Today we live in the dawn of the multimessenger era. The successes of the multimessenger campaigns of the last decade have pushed multimessenger science to the forefront of priority science areas in both the particle physics and the astrophysics communities. Multimessenger science provides new methods of testing fundamental theories about the nature of matter and energy, particularly in conditions that are not reproducible on Earth. This white paper will present the science and facilities that will provide opportunities for the particle physics community renew its commitment and maintain its leadership in multimessenger science. △ Less

Submitted 18 March, 2022; originally announced March 2022.

Comments: 174 pages, 12 figures. Contribution to Snowmass 2021. Solicited white paper from CF07. Comments and endorsers welcome. Still accepting contributions (contact editors)

Snowmass2021 Cosmic Frontier White Paper: Future Gravitational-Wave Detector Facilities

Authors: Stefan W. Ballmer, Rana Adhikari, Leonardo Badurina, Duncan A. Brown, Swapan Chattopadhyay, Matthew Evans, Peter Fritschel, Evan Hall, Jason M. Hogan, Karan Jani, Tim Kovachy, Kevin Kuns, Ariel Schwartzman, Daniel Sigg, Bram Slagmolen, Salvatore Vitale, Christopher Wipf

Abstract: The next generation of gravitational-wave observatories can explore a wide range of fundamental physics phenomena throughout the history of the universe. These phenomena include access to the universe's binary black hole population throughout cosmic time, to the universe's expansion history independent of the cosmic distance ladders, to stochastic gravitational-waves from early-universe phase tran… ▽ More The next generation of gravitational-wave observatories can explore a wide range of fundamental physics phenomena throughout the history of the universe. These phenomena include access to the universe's binary black hole population throughout cosmic time, to the universe's expansion history independent of the cosmic distance ladders, to stochastic gravitational-waves from early-universe phase transitions, to warped space-time in the strong-field and high-velocity limit, to the equation of state of nuclear matter at neutron star and post-merger densities, and to dark matter candidates through their interaction in extreme astrophysical environments or their interaction with the detector itself. We present the gravitational-wave detector concepts than can drive the future of gravitational-wave astrophysics. We summarize the status of the necessary technology, and the research needed to be able to build these observatories in the 2030s. △ Less

Submitted 29 March, 2022; v1 submitted 15 March, 2022; originally announced March 2022.

Comments: 31 pages, 5 figures, contribution to Snowmass 2021

Snowmass2021 Cosmic Frontier White Paper: Observational Facilities to Study Dark Matter

Authors: Sukanya Chakrabarti, Alex Drlica-Wagner, Ting S. Li, Neelima Sehgal, Joshua D. Simon, Simon Birrer, Duncan A. Brown, Rebecca Bernstein, Alberto D. Bolatto, Philip Chang, Kyle Dawson, Paul Demorest, Daniel Grin, David L. Kaplan, Joseph Lazio, Jennifer Marshall, Eric J. Murphy, Scott Ransom, Brant E. Robertson, Rajeev Singh, Anže Slosar, Tommaso Treu, Yu-Dai Tsai, Benjamin F. Williams

Abstract: We present an overview of future observational facilities that will significantly enhance our understanding of the fundamental nature of dark matter. These facilities span a range of observational techniques including optical/near-infrared imaging and spectroscopy, measurements of the cosmic microwave background, pulsar timing, 21-cm observations of neutral hydrogen at high redshift, and the measu… ▽ More We present an overview of future observational facilities that will significantly enhance our understanding of the fundamental nature of dark matter. These facilities span a range of observational techniques including optical/near-infrared imaging and spectroscopy, measurements of the cosmic microwave background, pulsar timing, 21-cm observations of neutral hydrogen at high redshift, and the measurement of gravitational waves. Such facilities are a critical component of a multi-pronged experimental program to uncover the nature of dark matter, while often providing complementary measurements of dark energy, neutrino physics, and inflation. △ Less

Submitted 11 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021. 31 pages, 6 figures, 2 tables

A Kernel-Based Approach for Modelling Gaussian Processes with Functional Information

Authors: John Nicholson, Peter Kiessler, D. Andrew Brown

Abstract: Gaussian processes are among the most useful tools in modeling continuous processes in machine learning and statistics. If the value of a process is known at a finite collection of points, one may use Gaussian processes to construct a surface which interpolates these values to be used for prediction and uncertainty quantification in other locations. However, it is not always the case that the avai… ▽ More Gaussian processes are among the most useful tools in modeling continuous processes in machine learning and statistics. If the value of a process is known at a finite collection of points, one may use Gaussian processes to construct a surface which interpolates these values to be used for prediction and uncertainty quantification in other locations. However, it is not always the case that the available information is in the form of a finite collection of points. For example, boundary value problems contain information on the boundary of a domain, which is an uncountable collection of points that cannot be incorporated into typical Gaussian process techniques. In this paper we construct a Gaussian process model which utilizes reproducing kernel Hilbert spaces to unify the typical finite case with the case of having uncountable information by exploiting the equivalence of conditional expectation and orthogonal projections. We discuss this construction in statistical models, including numerical considerations and a proof of concept. △ Less

Submitted 26 January, 2022; originally announced January 2022.

Comments: 21 pages, 4 figures

Next Generation Observatories -- Report from the Dawn VI Workshop; October 5-7 2021

Authors: D. H. Shoemaker, Stefan Ballmer, Matteo Barsuglia, E. Berger, Emanuele Berti, Duncan A. Brown, Poonam Chandra, Matthew Evans, Ke Fang, Wen-fai Fong, Andreas Freise, Peter Fritschel, Jenny Greene, C. J. Horowitz, Jeff Kissel, Brian Lantz, Paul D. Lasky, Harald Lueck, M. Coleman Miller, Alexander H. Nitz, David Ottaway, Hiranya V. Peiris, Michele Punturo, D. H. Reitze, Gary H. Sanders , et al. (11 additional authors not shown)

Abstract: The workshop Dawn VI: Next Generation Observatories took place online over three days, 5-7 October, 2021. More than 200 physicists and astronomers attended to contribute to, and learn from, a discussion of next-generation ground-based gravitational-wave detectors. The program was centered on the next generation of ground-based gravitational-wave observatories and their synergy with the greater lan… ▽ More The workshop Dawn VI: Next Generation Observatories took place online over three days, 5-7 October, 2021. More than 200 physicists and astronomers attended to contribute to, and learn from, a discussion of next-generation ground-based gravitational-wave detectors. The program was centered on the next generation of ground-based gravitational-wave observatories and their synergy with the greater landscape of scientific observatories of the 2030s. Cosmic Explorer (CE), a concept developed with US National Science Foundation support, was a particular focus; Einstein Telescope (ET), the European next generation concept, is an important complement and partner in forming a network. The concluding summary of the meeting expressed the sentiment that the observational science accessible to CE and ET, also in combination with data from other non-GW observatories, will stimulate a very broad community of analysts and yield insights which are exciting given the access to GWs from the entire universe. The need, and desire, for closer collaboration between ET and CE was expressed; a three-detector network is optimal for delivering much of the science. The science opportunities afforded by CE and ET are broad and compelling, impacting a wide range of disciplines in physics and high energy astrophysics. There was a consensus that CE is a concept that can deliver the promised science. A strong endorsement of Cosmic Explorer, as described in the CE Horizon Study, is a primary outcome of DAWN VI. △ Less

Submitted 20 February, 2022; v1 submitted 23 December, 2021; originally announced December 2021.

Comments: Proceedings of the Dawn VI workshop

Bayesian Regularization for Functional Graphical Models with Applications to Neuroimaging

Authors: Jiajing Niu, Boyoung Hur, John Absher, D. Andrew Brown

Abstract: Graphical models, used to express conditional dependence between random variables observed at various nodes, are used extensively in many fields such as genetics, neuroscience, and social network analysis. While most current statistical methods for estimating graphical models focus on scalar data, there is interest in estimating analogous dependence structures when the data observed at each node a… ▽ More Graphical models, used to express conditional dependence between random variables observed at various nodes, are used extensively in many fields such as genetics, neuroscience, and social network analysis. While most current statistical methods for estimating graphical models focus on scalar data, there is interest in estimating analogous dependence structures when the data observed at each node are functional, such as signals or images. In this paper, we propose a fully Bayesian regularization scheme for estimating functional graphical models. We first consider a direct Bayesian analog of the functional graphical lasso proposed by Qiao et al. (2019). We then propose a regularization strategy via the graphical horseshoe. We compare these approaches via simulation study and apply our proposed functional graphical horseshoe to two motivating applications, electroencephalography data for comparing brain activation between an alcoholic group and controls, as well as changes in structural connectivity in the presence of traumatic brain injury (TBI). Our results yield insight into how the brain attempts to compensate for disconnected networks after injury. △ Less

Submitted 21 October, 2022; v1 submitted 11 October, 2021; originally announced October 2021.

Comments: 29 pages, 3 figures

A Horizon Study for Cosmic Explorer: Science, Observatories, and Community

Authors: Matthew Evans, Rana X Adhikari, Chaitanya Afle, Stefan W. Ballmer, Sylvia Biscoveanu, Ssohrab Borhanian, Duncan A. Brown, Yanbei Chen, Robert Eisenstein, Alexandra Gruson, Anuradha Gupta, Evan D. Hall, Rachael Huxford, Brittany Kamai, Rahul Kashyap, Jeff S. Kissel, Kevin Kuns, Philippe Landry, Amber Lenon, Geoffrey Lovelace, Lee McCuller, Ken K. Y. Ng, Alexander H. Nitz, Jocelyn Read, B. S. Sathyaprakash , et al. (7 additional authors not shown)

Abstract: This Horizon Study describes a next-generation ground-based gravitational-wave observatory: Cosmic Explorer. With ten times the sensitivity of Advanced LIGO, Cosmic Explorer will push gravitational-wave astronomy towards the edge of the observable universe ($z \sim 100$). The goals of this Horizon Study are to describe and evaluate design concepts for Cosmic Explorer; to plan for the United States… ▽ More This Horizon Study describes a next-generation ground-based gravitational-wave observatory: Cosmic Explorer. With ten times the sensitivity of Advanced LIGO, Cosmic Explorer will push gravitational-wave astronomy towards the edge of the observable universe ($z \sim 100$). The goals of this Horizon Study are to describe and evaluate design concepts for Cosmic Explorer; to plan for the United States' leadership in gravitational-wave astronomy; and to envisage the role of Cosmic Explorer in the international effort to build a "Third-Generation" (3G) observatory network that will make discoveries transformative across astronomy, physics, and cosmology. △ Less

Submitted 6 October, 2021; v1 submitted 20 September, 2021; originally announced September 2021.

Comments: 173 pages

Report number: CE-P2100003-v7

Evaluation of $^{86}$Kr Cross Sections For Use in Fusion Diagnostics

Authors: M. Vorabbi, G. P. A. Nobre, D. A. Brown, A. M. Lewis, E. Rubino, S. Mughabghab

Abstract: The National Ignition Facility at Lawrence Livermore National Laboratory uses $^{86}$Kr as a diagnostic tool to measure the neutron flux produced by fusion reactions. As krypton is chemically inert, it can be implanted directly into the fuel capsule, and the reaction products can be measured to determine the flux of fusion neutrons. $^{86}$Kr cross sections also provide model constraints for the… ▽ More The National Ignition Facility at Lawrence Livermore National Laboratory uses $^{86}$Kr as a diagnostic tool to measure the neutron flux produced by fusion reactions. As krypton is chemically inert, it can be implanted directly into the fuel capsule, and the reaction products can be measured to determine the flux of fusion neutrons. $^{86}$Kr cross sections also provide model constraints for the $^{85}$Kr branching point in the s-process and the neutron flux in stars. In this work, experimental data on the neutron production, radiative capture, inelastic scattering, and total cross sections of $^{86}$Kr were used in conjunction with the fast region nuclear reaction code EMPIRE and a new resonance-region evaluation to produce a new evaluation of neutron-induced reactions on $^{86}$Kr. For the EMPIRE calculations, we fitted the optical model potential up to 12 MeV to simultaneously reproduce the experimental data for the total cross section and the main inelastic gamma transition from the $2^+$ state to the $0^+$ ground state. For energies above 12 MeV, due to large fluctuations and uncertainties in the total cross section data, we preferred to adopt the Koning-Delaroche global spherical optical model potential. With these models and corrections to the structure of $^{86}$Kr, the evaluated cross sections matched the experimental data. The new evaluation has been submitted for incorporation in the next release of the ENDF/B nuclear reaction library. △ Less

Submitted 16 September, 2021; originally announced September 2021.

Eccentric Binary Neutron Star Search Prospects for Cosmic Explorer

Authors: Amber K. Lenon, Duncan A. Brown, Alexander H. Nitz

Abstract: We determine the ability of Cosmic Explorer, a proposed third-generation gravitational-wave observatory, to detect eccentric binary neutron stars and to measure their eccentricity. We find that for a matched-filter search, template banks constructed using binaries in quasi-circular orbits are effectual for eccentric neutron star binaries with $e_{7} \leq 0.004$ ($e_{7} \leq 0.003$) for CE1 (CE2),… ▽ More We determine the ability of Cosmic Explorer, a proposed third-generation gravitational-wave observatory, to detect eccentric binary neutron stars and to measure their eccentricity. We find that for a matched-filter search, template banks constructed using binaries in quasi-circular orbits are effectual for eccentric neutron star binaries with $e_{7} \leq 0.004$ ($e_{7} \leq 0.003$) for CE1 (CE2), where $e_7$ is the binary's eccentricity at a gravitational-wave frequency of 7~Hz. We show that stochastic template placement can be used to construct a matched-filter search for binaries with larger eccentricities and construct an effectual template bank for binaries with $e_{7} \leq 0.05$. We show that the computational cost of both the search for binaries in quasi-circular orbits and eccentric orbits is not significantly larger for Cosmic Explorer than for Advanced LIGO and is accessible with present-day computational resources. We investigate Cosmic Explorer's ability to distinguish between circular and eccentric binaries. We estimate that for a binary with a signal-to-noise ratio of 8 (800), Cosmic Explorer can distinguish between a circular binary and a binary with eccentricity $e_7 \gtrsim 10^{-2}$ ($10^{-3}$) at 90% confidence. △ Less

Submitted 17 September, 2021; v1 submitted 25 March, 2021; originally announced March 2021.

Comments: 8 pages, 7 figures, supplemental materials at https://github.com/gwastro/cosmic-explorer-bns-eccentricity, revised to match version accepted by PRD

Journal ref: Phys. Rev. D 104, 063011 (2021)

Newly Evaluated Neutron Reaction Data on Chromium Isotopes

Authors: G. P. A. Nobre, M. T. Pigni, D. A. Brown, R. Capote, A. Trkov, K. H. Guber, R. Arcilla, J. Gutierrez, A. Cuadra, G. Arbanas, B. Kos, D. Bernard, P. Leconte

Abstract: Neutron reaction data for the set of major chromium isotopes were reevaluated from the thermal energy range up to 20 MeV. In the low energy region, updates to the thermal values together with an improved $R$-matrix analysis of the resonance parameters characterizing the cluster of large $s$-wave resonances for $^{50,53}$Cr isotopes were performed. In the intermediate and high energy range up to 20… ▽ More Neutron reaction data for the set of major chromium isotopes were reevaluated from the thermal energy range up to 20 MeV. In the low energy region, updates to the thermal values together with an improved $R$-matrix analysis of the resonance parameters characterizing the cluster of large $s$-wave resonances for $^{50,53}$Cr isotopes were performed. In the intermediate and high energy range up to 20 MeV, the evaluation methodology used statistical nuclear reaction models implemented in the EMPIRE code within the Hauser-Feshbach framework to evaluate the reaction cross sections and angular distributions. Exceptionally, experimental data were used to evaluate relevant cross sections above the resonance region up to 5 MeV in the major $^{52}$Cr isotope. Evaluations were benchmarked with Monte Carlo simulations of a small suite of critical assemblies highly sensitive to Chromium data, and with the Oktavian shielding benchmark to judge deep penetration performance with a 14-MeV D-T neutron source. A significant improvement in performance is demonstrated compared to existing evaluations. △ Less

Submitted 23 February, 2021; originally announced February 2021.

Comments: Submitted to Nuclear Data Sheets. 41 pages, 44 Figures, 13 Tables

Journal ref: Nuclear Data Sheets Volume 173, March-April 2021, Pages 1-41

Reproducing GW150914: the first observation of gravitational waves from a binary black hole merger

Authors: Duncan A. Brown, Karan Vahi, Michela Taufer, Von Welch, Ewa Deelman

Abstract: In 2016, LIGO and Virgo announced the first observation of gravitational waves from a binary black hole merger, known as GW150914. To establish the confidence of this detection, large-scale scientific workflows were used to measure the event's statistical significance. They used code written by the LIGO/Virgo and were executed on the LIGO Data Grid. The codes are publicly available, but there has… ▽ More In 2016, LIGO and Virgo announced the first observation of gravitational waves from a binary black hole merger, known as GW150914. To establish the confidence of this detection, large-scale scientific workflows were used to measure the event's statistical significance. They used code written by the LIGO/Virgo and were executed on the LIGO Data Grid. The codes are publicly available, but there has not yet been an attempt to directly reproduce the results, although several analyses have replicated the analysis, confirming the detection. We attempt to reproduce the result presented in the GW150914 discovery paper using publicly available code on the Open Science Grid. We show that we can reproduce the main result but we cannot exactly reproduce the LIGO analysis as the original data set used is not public. We discuss the challenges we encountered and make recommendations for scientists who wish to make their work reproducible. △ Less

Submitted 2 March, 2021; v1 submitted 14 October, 2020; originally announced October 2020.

Comments: 11 pages, 3 figures, 2 tables, 1 listing. Revised version accepted by IEEE Computing in Science and Engineering. Code available from https://doi.org/10.5281/zenodo.4085984 and data available from http://dx.doi.org/10.21227/c634-qh33

Inferring physical properties of stellar collapse by third-generation gravitational-wave detectors

Authors: Chaitanya Afle, Duncan A Brown

Abstract: Galactic core-collapse supernovae are among the possible sources of gravitational waves. We investigate the ability of gravitational-wave observatories to extract the properties of the collapsing progenitor from the gravitational waves radiated. We use simulations of supernovae that explore a variety of progenitor core rotation rates and nuclear equations of state and examine the ability of curren… ▽ More Galactic core-collapse supernovae are among the possible sources of gravitational waves. We investigate the ability of gravitational-wave observatories to extract the properties of the collapsing progenitor from the gravitational waves radiated. We use simulations of supernovae that explore a variety of progenitor core rotation rates and nuclear equations of state and examine the ability of current and future observatories to determine these properties using gravitational-wave parameter estimation. We use principal component analysis of the simulation catalog to determine the dominant features of the waveforms and create a map between the measured properties of the waveform and the physical properties of the progenitor. We use Bayesian parameter inference and the parameter map to calculate posterior probabilities for the physical properties given a gravitational-wave observation. We estimate the ratio of the progenitor's core rotational kinetic energy to potential energy ($β$) and the post bounce oscillation frequency. For a supernovae at the distance of the galactic center (8.1 kpc) with $β= 0.02$ our method can estimate $β$ with a $90\%$ credible interval of $0.004$ for Advanced LIGO, improving to $0.0008$ for Cosmic Explorer. We demonstrate that if the core is rotating sufficiently rapidly for a signal observed by Cosmic Explorer, our method can also extract the post bounce oscillation frequency of the protoneutron star to a precision of within $5$~Hz ($90\%$ credible interval) allowing us to constrain the nuclear equation of state. For a supernovae at the distance of the Magellanic Clouds (48.5 kpc) Cosmic Explorer's ability to measure these parameters decreases slightly to $0.003$ for rotation and $11$~Hz for the postbounce oscillation frequency ($90\%$ credible interval). Sources in Magellanic Clouds will be too distant for Advanced LIGO to measure these properties. △ Less

Submitted 27 January, 2021; v1 submitted 1 October, 2020; originally announced October 2020.

Comments: 10 pages, 9 figures. Updated to accepted version

Journal ref: Phys. Rev. D 103, 023005 (2021)

Fast Parameter Estimation of Binary Mergers for Multimessenger Followup

Authors: Daniel Finstad, Duncan A. Brown

Abstract: Significant human and observational resources have been dedicated to electromagnetic followup of gravitational-wave events detected by Advanced LIGO and Virgo. As the sensitivity of LIGO and Virgo improves, the rate of sources detected will increase. Margalit & Metzger (2019; arXiv:1904.11995) have suggested that it may be necessary to prioritize observations of future events. Optimal prioritizati… ▽ More Significant human and observational resources have been dedicated to electromagnetic followup of gravitational-wave events detected by Advanced LIGO and Virgo. As the sensitivity of LIGO and Virgo improves, the rate of sources detected will increase. Margalit & Metzger (2019; arXiv:1904.11995) have suggested that it may be necessary to prioritize observations of future events. Optimal prioritization requires a rapid measurement of a gravitational-wave event's masses and spins, as these can determine the nature of any electromagnetic emission. We extend the relative binning method of Zackay et al. (2018; arXiv:1806.08792) to a coherent detector-network statistic. We show that the method can be seeded from the output of a matched-filter search and used in a Bayesian parameter measurement framework to produce marginalized posterior probability densities for the source's parameters within 20 minutes of detection on 32 CPU cores. We demonstrate that this algorithm produces unbiased estimates of the parameters with the same accuracy as running parameter estimation using the standard gravitational-wave likelihood. We encourage the adoption of this method in future LIGO-Virgo observing runs to allow fast dissemination of the parameters of detected events so that the observing community can make best use of its resources. △ Less

Submitted 29 September, 2020; originally announced September 2020.

Comments: 11 pages, 5 figures

DQSEGDB: A time-interval database for storing gravitational wave observatory metadata

Authors: Ryan P. Fisher, Gary Hemming, Marie-Anne Bizouard, Duncan A. Brown, Peter F. Couvares, Florent Robinet, Didier Verkindt

Abstract: The Data Quality Segment Database (DQSEGDB) software is a database service, backend API, frontend graphical web interface, and client package used by the Laser Interferometer Gravitational-Wave Observatory (LIGO), Virgo, GEO600 and the Kamioka Gravitational wave detector for storing and accessing metadata describing the status of their detectors. The DQSEGDB has been used in the analysis of all pu… ▽ More The Data Quality Segment Database (DQSEGDB) software is a database service, backend API, frontend graphical web interface, and client package used by the Laser Interferometer Gravitational-Wave Observatory (LIGO), Virgo, GEO600 and the Kamioka Gravitational wave detector for storing and accessing metadata describing the status of their detectors. The DQSEGDB has been used in the analysis of all published detections of gravitational waves in the advanced detector era. The DQSEGDB currently stores roughly 600 million metadata entries and responds to roughly 600,000 queries per day with an average response time of 0.223 ms. △ Less

Submitted 25 August, 2020; originally announced August 2020.

Comments: 11 pages, 2 tables, 1 figure, submitted to SoftwareX

Measuring the eccentricity of GW170817 and GW190425

Authors: Amber K. Lenon, Alexander H. Nitz, Duncan A. Brown

Abstract: Two binary neutron star mergers, GW170817 and GW190425, have been detected by Advanced LIGO and Virgo. These signals were detected by matched-filter searches that assume the star's orbit has circularized by the time their gravitational-wave emission is observable. This suggests that their eccentricity is low, but full parameter estimation of their eccentricity has not yet been performed. We use gr… ▽ More Two binary neutron star mergers, GW170817 and GW190425, have been detected by Advanced LIGO and Virgo. These signals were detected by matched-filter searches that assume the star's orbit has circularized by the time their gravitational-wave emission is observable. This suggests that their eccentricity is low, but full parameter estimation of their eccentricity has not yet been performed. We use gravitational-wave observations to measure the eccentricity of GW170817 and GW190425. We find that the eccentricity at a gravitational-wave frequency of 10 Hz is $e \leq 0.024$ and $e \leq 0.048$ for GW170817 and GW190425, respectively (90% confidence). This is consistent with the binaries being formed in the field, as such systems are expected to have circularized to $e \leq 10^{-4}$ by the time they reach the LIGO-Virgo band. Our constraint is a factor of two smaller that an estimate based on GW170817 being detected by searches that neglect eccentricity. However, we caution that we find significant prior dependence in our limits, suggesting that there is limited information in the signals. We note that other techniques used to constrain binary neutron star eccentricity without full parameter estimation may miss degeneracies in the waveform, and that for future signals it will be important to perform full parameter estimation with accurate waveform templates. △ Less

Submitted 20 July, 2020; v1 submitted 28 May, 2020; originally announced May 2020.

Comments: 6 pages, 3 figures, 1 table, supplemental materials at https://github.com/gwastro/bns-eccentric-pe, revised to match version accepted by MNRAS

Search for Eccentric Binary Neutron Star Mergers in the first and second observing runs of Advanced LIGO

Authors: Alexander H. Nitz, Amber Lenon, Duncan A. Brown

Abstract: We present a search for gravitational waves from merging binary neutron stars which have non-negligible eccentricity as they enter the LIGO observing band. We use the public Advanced LIGO data which covers the period from 2015 through 2017 and contains $\sim164$ days of LIGO-Hanford and LIGO-Livingston coincident observing time. The search was conducted using matched-filtering using the PyCBC tool… ▽ More We present a search for gravitational waves from merging binary neutron stars which have non-negligible eccentricity as they enter the LIGO observing band. We use the public Advanced LIGO data which covers the period from 2015 through 2017 and contains $\sim164$ days of LIGO-Hanford and LIGO-Livingston coincident observing time. The search was conducted using matched-filtering using the PyCBC toolkit. We find no significant binary neutron star candidates beyond GW170817, which has previously been reported by searches for binaries in circular orbits. We place a 90% upper limit of $\sim1700$ mergers $\textrm{Gpc}^{-3} \textrm{Yr}^{-1}$ for eccentricities $\lesssim 0.43$ at a dominant-mode gravitational-wave frequency of 10 Hz. The absence of a detection with these data is consistent with theoretical predictions of eccentric binary neutron star merger rates. Using our measured rate we estimate the sensitive volume of future gravitational-wave detectors and compare this to theoretical rate predictions. We find that, in the absence of a prior detection, the rate limits set by six months of Cosmic Explorer observations would constrain all current plausible models of eccentric binary neutron star formation. △ Less

Submitted 19 December, 2019; v1 submitted 11 December, 2019; originally announced December 2019.

Comments: 7 pages, 3 figures, 1 table, supplemental materials at https://github.com/gwastro/eccentric-bns-search, revised to match version accepted by ApJ

Journal ref: 2020 ApJ 890 1

Detectability of the subdominant mode in a binary black hole ringdown

Authors: Swetha Bhagwat, Miriam Cabero, Collin D. Capano, Badri Krishnan, Duncan A. Brown

Abstract: The ringdown is the late part of the post-merger signature emitted during the coalescence of two black holes and comprises of a superposition of quasi-normal-modes. Within general relativity, because of the no-hair theorems, the frequencies and damping times of these modes are entirely determined by the mass and angular momentum of the final Kerr black hole. A detection of multiple ringdown modes… ▽ More The ringdown is the late part of the post-merger signature emitted during the coalescence of two black holes and comprises of a superposition of quasi-normal-modes. Within general relativity, because of the no-hair theorems, the frequencies and damping times of these modes are entirely determined by the mass and angular momentum of the final Kerr black hole. A detection of multiple ringdown modes would potentially allow us to test the no-hair theorem from observational data. The parameters which determine whether sub-dominant ringdown modes can be detected are primarily the overall signal-to-noise ratio present in the ringdown signal, and on the amplitude of the subdominant mode with respect to the dominant mode. In this paper, we use Bayesian inference to determine the detectability of a subdominant mode in a set of simulated analytical ringdown signals. Focusing on the design sensitivity of the Advanced LIGO detectors, we systematically vary the signal-to-noise ratio of the ringdown signal, and the mode amplitude ratio in order to determine what kind of signals are promising for performing black hole spectroscopy. △ Less

Submitted 30 October, 2019; v1 submitted 29 October, 2019; originally announced October 2019.

Comments: 7 figures, 2 tables

Journal ref: Phys. Rev. D 102, 024023 (2020)

Constraining level densities through quantitative correlations with cross-section data

Authors: G. P. A. Nobre, D. A. Brown, M. W. Herman, A. Golas

Abstract: The adopted level densities (LD) for the nuclei produced through different reaction mechanisms significantly impact the calculation of cross sections for the many reaction channels. Common LD models make simplified assumptions regarding the overall behavior of the total LD and the intrinsic spin and parity distributions of the excited states. However, very few experimental constraints are taken in… ▽ More The adopted level densities (LD) for the nuclei produced through different reaction mechanisms significantly impact the calculation of cross sections for the many reaction channels. Common LD models make simplified assumptions regarding the overall behavior of the total LD and the intrinsic spin and parity distributions of the excited states. However, very few experimental constraints are taken into account: LD at neutron separation energy coming from average resonance spacings, whenever they have been previously measured, and the sometimes subjective extrapolation of discrete levels. These, however, constrain the LD only for very specific spins, parities and excitation energies. This work aims to establish additional experimental constraints on LD through quantitative correlations between cross sections and LD. This allows for the fitting and determination of detailed structures in LD. For this we use the microscopic Hartree-Fock-Bogoliubov (HFB) LD to associate variations predicted by the model with the structure observed in double-differential spectra at low outgoing neutron energy, which is dominated by the LD input. We also use \nuc{56}{Fe} ($n,p$) as an example cross sections are extremely sensitive to LD. For comparison purposes we also perform calculations with the GC model. With this approach we are able to perform fits of the LD based on actual experimental data, constraining the model and ensuring its consistency. This approach can be particularly useful in extrapolating the LD to nuclei for which high-excited discrete levels and/or resonance spacings are unknown. It also predicts inelastic gamma cross sections that can significantly differ from more standard phenomenological LD. △ Less

Submitted 28 February, 2020; v1 submitted 20 September, 2019; originally announced September 2019.

Comments: 16 pages, 12 Figures. Accepted by PRC

Report number: BNL-213738-2020-JAAM, LA-UR-19-29483

Journal ref: Phys. Rev. C 101, 034608 (2020)

Stringent constraints on neutron-star radii from multimessenger observations and nuclear theory

Authors: Collin D. Capano, Ingo Tews, Stephanie M. Brown, Ben Margalit, Soumi De, Sumit Kumar, Duncan A. Brown, Badri Krishnan, Sanjay Reddy

Abstract: The properties of neutron stars are determined by the nature of the matter that they contain. These properties can be constrained by measurements of the star's size. We obtain stringent constraints on neutron-star radii by combining multimessenger observations of the binary neutron-star merger GW170817 with nuclear theory that best accounts for density-dependent uncertainties in the equation of st… ▽ More The properties of neutron stars are determined by the nature of the matter that they contain. These properties can be constrained by measurements of the star's size. We obtain stringent constraints on neutron-star radii by combining multimessenger observations of the binary neutron-star merger GW170817 with nuclear theory that best accounts for density-dependent uncertainties in the equation of state. We construct equations of state constrained by chiral effective field theory and marginalize over these using the gravitational-wave observations. Combining this with the electromagnetic observations of the merger remnant that imply the presence of a short-lived hyper-massive neutron star, we find that the radius of a $1.4\,\rm{M}_\odot$ neutron star is $R_{1.4\,\mathrm{M}_\odot} = 11.0^{+0.9}_{-0.6}~{\rm km}$ (90% credible interval). Using this constraint, we show that neutron stars are unlikely to be disrupted in neutron-star black-hole mergers; subsequently, such events will not produce observable electromagnetic emission. △ Less

Submitted 24 March, 2020; v1 submitted 27 August, 2019; originally announced August 2019.

Comments: 39 pages, 4 figures

Report number: INT-PUB-19-037, LA-UR-19-28442

Journal ref: Nature Astronomy, 2020

Coupling material and mechanical design processes via computer model calibration

Authors: Carl Ehrett, D. Andrew Brown, Evan Chodora, Christopher Kitchens, Sez Atamturktur

Abstract: Computer model calibration typically operates by choosing parameter values in a computer model so that the model output faithfully predicts reality. By using performance targets in place of observed data, we show that calibration techniques can be repurposed to wed engineering and material design, two processes that are traditionally carried out separately. This allows materials to be designed wit… ▽ More Computer model calibration typically operates by choosing parameter values in a computer model so that the model output faithfully predicts reality. By using performance targets in place of observed data, we show that calibration techniques can be repurposed to wed engineering and material design, two processes that are traditionally carried out separately. This allows materials to be designed with specific engineering targets in mind while quantifying the associated sources of uncertainty. We demonstrate our proposed approach by "calibrating" material design settings to performance targets for a wind turbine blade. △ Less

Submitted 13 January, 2020; v1 submitted 22 July, 2019; originally announced July 2019.

Comments: 20 pages, 7 figures. Supplementary material and computer code available from the authors

Cosmic Explorer: The U.S. Contribution to Gravitational-Wave Astronomy beyond LIGO

Authors: David Reitze, Rana X Adhikari, Stefan Ballmer, Barry Barish, Lisa Barsotti, GariLynn Billingsley, Duncan A. Brown, Yanbei Chen, Dennis Coyne, Robert Eisenstein, Matthew Evans, Peter Fritschel, Evan D. Hall, Albert Lazzarini, Geoffrey Lovelace, Jocelyn Read, B. S. Sathyaprakash, David Shoemaker, Joshua Smith, Calum Torrie, Salvatore Vitale, Rainer Weiss, Christopher Wipf, Michael Zucker

Abstract: This white paper describes the research and development needed over the next decade to realize "Cosmic Explorer," the U.S. node of a future third-generation detector network that will be capable of observing and characterizing compact gravitational-wave sources to cosmological redshifts. This white paper describes the research and development needed over the next decade to realize "Cosmic Explorer," the U.S. node of a future third-generation detector network that will be capable of observing and characterizing compact gravitational-wave sources to cosmological redshifts. △ Less

Submitted 10 July, 2019; originally announced July 2019.

Comments: Astro2020 Decadal Survey ground-based technology development white paper

Report number: LIGO-P1900316

Journal ref: 2019 BAAS 51(7) 035

Detection Prospects of Core-Collapse Supernovae with Supernova-Optimized Third-Generation Gravitational-wave Detectors

Authors: Varun Srivastava, Stefan Ballmer, Duncan A. Brown, Chaitanya Afle, Adam Burrows, David Radice, David Vartanyan

Abstract: We optimize the third-generation gravitational-wave detector to maximize the range to detect core-collapse supernovae. Based on three-dimensional simulations for core-collapse and the corresponding gravitational-wave waveform emitted, the corresponding detection range for these waveforms is limited to within our galaxy even in the era of third-generation detectors. The corresponding event rate is… ▽ More We optimize the third-generation gravitational-wave detector to maximize the range to detect core-collapse supernovae. Based on three-dimensional simulations for core-collapse and the corresponding gravitational-wave waveform emitted, the corresponding detection range for these waveforms is limited to within our galaxy even in the era of third-generation detectors. The corresponding event rate is two per century. We find from the waveforms that to detect core-collapse supernovae with an event rate of one per year, the gravitational-wave detectors need a strain sensitivity of 3$\times10^{-27}~$Hz$^{-1/2}$ in a frequency range from 100~Hz to 1500~Hz. We also explore detector configurations technologically beyond the scope of third-generation detectors. We find with these improvements, the event rate for gravitational-wave observations from CCSN is still low, but is improved to one in twenty years. △ Less

Submitted 31 May, 2019; originally announced June 2019.

Journal ref: Phys. Rev. D 100, 043026 (2019)

Constraining level densities using spectral data

Authors: G. P. A. Nobre, D. A. Brown, M. W. Herman

Abstract: Several models of level densities exist and they often make simplified assumptions regarding the overall behavior of the total level densities (LD) and the intrinsic spin and parity distributions of the excited states. Normally, such LD models are constrained only by the measured $D_0$, i.e. the density of levels at the neutron separation energy of the compound nucleus (target plus neutron), and t… ▽ More Several models of level densities exist and they often make simplified assumptions regarding the overall behavior of the total level densities (LD) and the intrinsic spin and parity distributions of the excited states. Normally, such LD models are constrained only by the measured $D_0$, i.e. the density of levels at the neutron separation energy of the compound nucleus (target plus neutron), and the sometimes subjective extrapolation of discrete levels. In this work we use microscopic Hartree-Fock-Bogoliubov (HFB) level densities, which intrinsically provide more realistic spin and parity distributions, and associate variations predicted by the HFB model with the observed double-differential cross sections at low outgoing neutron energy, region that is dominated by the LD input. With this approach we are able to perform fits of the LD based on actual experimental data, constraining the model and ensuring its consistency. This approach can be particularly useful in extrapolating the LD to nuclei for which high-excited discrete levels and/or values of $D_0$ are unknown. It also predicts inelastic gamma (n,n$^{\prime}γ$) cross sections that in some cases can differ significantly from more standard LD models such as Gilbert-Cameron. △ Less

Submitted 22 May, 2019; originally announced May 2019.

Comments: 6 pages, 4 figures. Submitted to the proceedings of the CNR*18 conference

Report number: BNL-211548-2019-JAAM

Journal ref: In: Escher J. et al. (eds) Compound-Nuclear Reactions. Springer Proceedings in Physics, vol 254. Springer, Cham., pp 133-138, 2020

Custom Execution Environments with Containers in Pegasus-enabled Scientific Workflows

Authors: Karan Vahi, Mats Rynge, George Papadimitriou, Duncan A. Brown, Rajiv Mayani, Rafael Ferreira da Silva, Ewa Deelman, Anirban Mandal, Eric Lyons, Michael Zink

Abstract: Science reproducibility is a cornerstone feature in scientific workflows. In most cases, this has been implemented as a way to exactly reproduce the computational steps taken to reach the final results. While these steps are often completely described, including the input parameters, datasets, and codes, the environment in which these steps are executed is only described at a higher level with end… ▽ More Science reproducibility is a cornerstone feature in scientific workflows. In most cases, this has been implemented as a way to exactly reproduce the computational steps taken to reach the final results. While these steps are often completely described, including the input parameters, datasets, and codes, the environment in which these steps are executed is only described at a higher level with endpoints and operating system name and versions. Though this may be sufficient for reproducibility in the short term, systems evolve and are replaced over time, breaking the underlying workflow reproducibility. A natural solution to this problem is containers, as they are well defined, have a lifetime independent of the underlying system, and can be user-controlled so that they can provide custom environments if needed. This paper highlights some unique challenges that may arise when using containers in distributed scientific workflows. Further, this paper explores how the Pegasus Workflow Management System implements container support to address such challenges. △ Less

Submitted 20 May, 2019; originally announced May 2019.

Comments: 10 pages, 7 figures, submitted to eScience 2019

Posterior samples of the parameters of binary black holes from Advanced LIGO, Virgo's second observing run

Authors: Soumi De, Christopher M. Biwer, Collin D. Capano, Alexander H. Nitz, Duncan A. Brown

Abstract: This paper presents a parameter estimation analysis of the seven binary black hole mergers---GW170104, GW170608, GW170729, GW170809, GW170814, GW170818, and GW170823---detected during the second observing run of the Advanced LIGO and Virgo observatories using the gravitational-wave open data. We describe the methodology for parameter estimation of compact binaries using gravitational-wave data, an… ▽ More This paper presents a parameter estimation analysis of the seven binary black hole mergers---GW170104, GW170608, GW170729, GW170809, GW170814, GW170818, and GW170823---detected during the second observing run of the Advanced LIGO and Virgo observatories using the gravitational-wave open data. We describe the methodology for parameter estimation of compact binaries using gravitational-wave data, and we present the posterior distributions of the inferred astrophysical parameters. We release our samples of the posterior probability density function with tutorials on using and replicating our results presented in this paper. △ Less

Submitted 3 June, 2019; v1 submitted 22 November, 2018; originally announced November 2018.

Comments: 15 pages, 4 figures

Journal ref: Scientific Data 6, 81 (2019)

Investigating the noise residuals around the gravitational wave event GW150914

Authors: Alex B. Nielsen, Alexander H. Nitz, Collin D. Capano, Duncan A. Brown

Abstract: We use the Pearson cross-correlation statistic proposed by Liu and Jackson, and employed by Creswell et al., to look for statistically significant correlations between the LIGO Hanford and Livingston detectors at the time of the binary black hole merger GW150914. We compute this statistic for the calibrated strain data released by LIGO, using both the residuals provided by LIGO and using our own s… ▽ More We use the Pearson cross-correlation statistic proposed by Liu and Jackson, and employed by Creswell et al., to look for statistically significant correlations between the LIGO Hanford and Livingston detectors at the time of the binary black hole merger GW150914. We compute this statistic for the calibrated strain data released by LIGO, using both the residuals provided by LIGO and using our own subtraction of a maximum-likelihood waveform that is constructed to model binary black hole mergers in general relativity. To assign a significance to the values obtained, we calculate the cross-correlation of both simulated Gaussian noise and data from the LIGO detectors at times during which no detection of gravitational waves has been claimed. We find that after subtracting the maximum likelihood waveform there are no statistically significant correlations between the residuals of the two detectors at the time of GW150914. △ Less

Submitted 12 February, 2019; v1 submitted 9 November, 2018; originally announced November 2018.

Comments: 14 pages, 7 figures. Minor text and figure changes in final v3. Notebooks for generating the results are available at https://github.com/gwastro/gw150914_investigation

1-OGC: The first open gravitational-wave catalog of binary mergers from analysis of public Advanced LIGO data

Authors: Alexander H. Nitz, Collin Capano, Alex B. Nielsen, Steven Reyes, Rebecca White, Duncan A. Brown, Badri Krishnan

Abstract: We present the first Open Gravitational-wave Catalog (1-OGC), obtained by using the public data from Advanced LIGO's first observing run to search for compact-object binary mergers. Our analysis is based on new methods that improve the separation between signals and noise in matched-filter searches for gravitational waves from the merger of compact objects. The three most significant signals in ou… ▽ More We present the first Open Gravitational-wave Catalog (1-OGC), obtained by using the public data from Advanced LIGO's first observing run to search for compact-object binary mergers. Our analysis is based on new methods that improve the separation between signals and noise in matched-filter searches for gravitational waves from the merger of compact objects. The three most significant signals in our catalog correspond to the binary black hole mergers GW150914, GW151226, and LVT151012. We assume a common population of binary black holes for these three signals by defining a region of parameter space that is consistent with these events. Under this assumption, we find that LVT151012 has a 97.6\% probability of being astrophysical in origin. No other significant binary black hole candidates are found, nor did we observe any significant binary neutron star or neutron star--black hole candidates. We make available our complete catalog of events, including the sub-threshold population of candidates. △ Less

Submitted 27 February, 2019; v1 submitted 5 November, 2018; originally announced November 2018.

Comments: 9 pages, 3 figures, 2 tables. Updated to match published version. Supplementary data at http://www.github.com/gwastro/1-ogc

Efficient Marginalization-based MCMC Methods for Hierarchical Bayesian Inverse Problems

Authors: Arvind K. Saibaba, Johnathan Bardsley, D. Andrew Brown, Alen Alexanderian

Abstract: Hierarchical models in Bayesian inverse problems are characterized by an assumed prior probability distribution for the unknown state and measurement error precision, and hyper-priors for the prior parameters. Combining these probability models using Bayes' law often yields a posterior distribution that cannot be sampled from directly, even for a linear model with Gaussian measurement error and Ga… ▽ More Hierarchical models in Bayesian inverse problems are characterized by an assumed prior probability distribution for the unknown state and measurement error precision, and hyper-priors for the prior parameters. Combining these probability models using Bayes' law often yields a posterior distribution that cannot be sampled from directly, even for a linear model with Gaussian measurement error and Gaussian prior. Gibbs sampling can be used to sample from the posterior, but problems arise when the dimension of the state is large. This is because the Gaussian sample required for each iteration can be prohibitively expensive to compute, and because the statistical efficiency of the Markov chain degrades as the dimension of the state increases. The latter problem can be mitigated using marginalization-based techniques, but these can be computationally prohibitive as well. In this paper, we combine the low-rank techniques of Brown, Saibaba, and Vallelian (2018) with the marginalization approach of Rue and Held (2005). We consider two variants of this approach: delayed acceptance and pseudo-marginalization. We provide a detailed analysis of the acceptance rates and computational costs associated with our proposed algorithms, and compare their performances on two numerical test cases---image deblurring and inverse heat equation. △ Less

Submitted 6 June, 2019; v1 submitted 2 November, 2018; originally announced November 2018.

Comments: 27 pages, 8 figures, 5 tables

Constraints on nonlinear tides due to $p$-$g$ mode coupling from the neutron-star merger GW170817

Authors: Steven Reyes, Duncan A. Brown

Abstract: It has been suggested by Weinberg et al. (2013) that an instability due to the nonlinear coupling of a neutron star's tide to its $p$- and $g$-modes could affect the gravitational-wave phase evolution of a neutron-star binary. Weinberg (2016) suggests that this instability can turn on as the gravitational-waves pass through the sensitive band of ground-based detectors, although the size of the eff… ▽ More It has been suggested by Weinberg et al. (2013) that an instability due to the nonlinear coupling of a neutron star's tide to its $p$- and $g$-modes could affect the gravitational-wave phase evolution of a neutron-star binary. Weinberg (2016) suggests that this instability can turn on as the gravitational-waves pass through the sensitive band of ground-based detectors, although the size of the effect is not known. The discovery of the binary neutron star merger GW170817 provides an opportunity to look for evidence of nonlinear tides from $p$-$g$ mode coupling. We compute Bayesian evidences that compare waveform models that include the $p$-$g$ mode coupling to models that do not. Using the waveform model and priors of Essick et al. (2016) we find that the observation of GW170817 is consistent with $p$-$g$ mode coupling, in agreement with Abbott, B. P., et al. (2019a). We investigate the properties of the model and find that this consistency is due to degeneracy in a large region of the parameter space between the the model that includes nonlinear tides from $p$-$g$ mode coupling and the standard post-Newtonian model that does not. We investigate the consistency of the GW170817 signal with regions of the parameter space where the effect of nonlinear tides is not degenerate with the standard model. Regions of the nonlinear tide parameter space that have a fitting factor of less than 99$\%$ (98.5$\%$) are disfavored by a Bayes factor of 15 (25). We conclude that the consistency of the GW170817 signal with the model of Essick et al. (2016) is due to parameter degeneracy and that regions where nonlinear tides produce a measurable effect are strongly disfavored. △ Less

Submitted 7 January, 2020; v1 submitted 21 August, 2018; originally announced August 2018.

Comments: 15 pages, 3 figures, submitted to ApJ, posteriors and evidences available from https://github.com/sugwg/gw170817-pg-modes

PyCBC Inference: A Python-based parameter estimation toolkit for compact binary coalescence signals

Authors: C. M. Biwer, Collin D. Capano, Soumi De, Miriam Cabero, Duncan A. Brown, Alexander H. Nitz, V. Raymond

Abstract: We introduce new modules in the open-source PyCBC gravitational- wave astronomy toolkit that implement Bayesian inference for compact-object binary mergers. We review the Bayesian inference methods implemented and describe the structure of the modules. We demonstrate that the PyCBC Inference modules produce unbiased estimates of the parameters of a simulated population of binary black hole mergers… ▽ More We introduce new modules in the open-source PyCBC gravitational- wave astronomy toolkit that implement Bayesian inference for compact-object binary mergers. We review the Bayesian inference methods implemented and describe the structure of the modules. We demonstrate that the PyCBC Inference modules produce unbiased estimates of the parameters of a simulated population of binary black hole mergers. We show that the posterior parameter distributions obtained used our new code agree well with the published estimates for binary black holes in the first LIGO-Virgo observing run. △ Less

Submitted 26 July, 2018; originally announced July 2018.

Impact of alternative transmission coefficient parameterizations on Hauser-Feshbach theory

Authors: David A. Brown, Gustavo P. A. Nobre, Michal W. Herman

Abstract: We investigate different formulations of the transmission coefficient $T_c$, including the form implied by Moldauer's ``sum rule for resonance reactions'' [P.A. Moldauer, Phys. Rev. Lett. 19, 1047 (1967)], the SPRT method [G. Noguere, et al. EPJ Web Conf. 146, 02036 (2017)] and the Moldauer-Simonius form [M. Simonius, Phys. Lett. 52B, 279 (1974); P.A. Moldauer, Phys. Rev. 157, 907 (1967)]. Within… ▽ More We investigate different formulations of the transmission coefficient $T_c$, including the form implied by Moldauer's ``sum rule for resonance reactions'' [P.A. Moldauer, Phys. Rev. Lett. 19, 1047 (1967)], the SPRT method [G. Noguere, et al. EPJ Web Conf. 146, 02036 (2017)] and the Moldauer-Simonius form [M. Simonius, Phys. Lett. 52B, 279 (1974); P.A. Moldauer, Phys. Rev. 157, 907 (1967)]. Within these different formulations, we compute the neutron transmission coefficients in the resolved and unresolved resonance regions, allowing a direct comparison with the transmission coefficients computed using an optical model potential. For nuclei for which there are no measured resonances, these approaches allow one to predict the average neutron resonance parameters directly from the optical model and level densities. Some of the approaches are valid in both the strong and weak coupling limits (i.e., any value of the average width and mean level spacing). Finally, both the Moldauer-Simonius and Moldauer's Sum Rule forms approaches suggest that superradiance, that is, the quantum chaotic enhancement of certain channels, may be a common phenomena in nuclear collisions. Our results suggest why superradiance has been previously overlooked. We apply our approach to neutron reactions on the closed shell $^{90}$Zr nucleus and the mid-shell $^{197}$Au nucleus. △ Less

Submitted 13 July, 2018; originally announced July 2018.

Comments: 12 pages, 6 figures, to be submitted to Phys. Rev. C

Report number: BNL-207837-2018-JAAM

Journal ref: Phys. Rev. C 98, 024616 (2018)

Tidal Deformabilities and Radii of Neutron Stars from the Observation of GW170817

Authors: Soumi De, Daniel Finstad, James M. Lattimer, Duncan A. Brown, Edo Berger, Christopher M. Biwer

Abstract: We use gravitational-wave observations of the binary neutron star merger GW170817 to explore the tidal deformabilities and radii of neutron stars. We perform Bayesian parameter estimation with the source location and distance informed by electromagnetic observations. We also assume that the two stars have the same equation of state; we demonstrate that for stars with masses comparable to the compo… ▽ More We use gravitational-wave observations of the binary neutron star merger GW170817 to explore the tidal deformabilities and radii of neutron stars. We perform Bayesian parameter estimation with the source location and distance informed by electromagnetic observations. We also assume that the two stars have the same equation of state; we demonstrate that for stars with masses comparable to the component masses of GW170817, this is effectively implemented by assuming that the stars' dimensionless tidal deformabilities are determined by the binary's mass ratio $q$ by $Λ_1/Λ_2 = q^6$. We investigate different choices of prior on the component masses of the neutron stars. We find that the tidal deformability and 90$\%$ credible interval is $\tildeΛ=222^{+420}_{-138}$ for a uniform component mass prior, $\tildeΛ=245^{+453}_{-151}$ for a component mass prior informed by radio observations of Galactic double neutron stars, and $\tildeΛ=233^{+448}_{-144}$ for a component mass prior informed by radio pulsars. We find a robust measurement of the common areal radius of the neutron stars across all mass priors of $8.9 \le \hat{R} \le 13.2$ km, with a mean value of $\langle \hat{R} \rangle = 10.8$ km. Our results are the first measurement of tidal deformability with a physical constraint on the star's equation of state and place the first lower bounds on the deformability and areal radii of neutron stars using gravitational waves. △ Less

Submitted 30 October, 2018; v1 submitted 23 April, 2018; originally announced April 2018.

Comments: Main paper: 6 pages, 3 figures; Supplemental Material: 5 pages, 6 figures, Erratum 2 pages, 1 figure

Journal ref: Phys. Rev. Lett. 121, 091102 (2018)

Measuring the viewing angle of GW170817 with electromagnetic and gravitational waves

Authors: Daniel Finstad, Soumi De, Duncan A. Brown, Edo Berger, Christopher M. Biwer

Abstract: The joint detection of gravitational waves (GWs) and electromagnetic (EM) radiation from the binary neutron star merger GW170817 ushered in a new era of multi-messenger astronomy. Joint GW-EM observations can be used to measure the parameters of the binary with better precision than either observation alone. Here, we use joint GW-EM observations to measure the viewing angle of GW170817, the angle… ▽ More The joint detection of gravitational waves (GWs) and electromagnetic (EM) radiation from the binary neutron star merger GW170817 ushered in a new era of multi-messenger astronomy. Joint GW-EM observations can be used to measure the parameters of the binary with better precision than either observation alone. Here, we use joint GW-EM observations to measure the viewing angle of GW170817, the angle between the binary's angular momentum and the line of sight. We combine a direct measurement of the distance to the host galaxy of GW170817 (NGC 4993) of $40.7\pm 2.36$ Mpc with the Laser Interferometer Gravitational-wave Observatory (LIGO)/Virgo GW data and find that the viewing angle is $32^{+10}_{-13}\,\pm 1.7$ degrees (90% confidence, statistical, and systematic errors). We place a conservative lower limit on the viewing angle of $\ge 13^\circ$, which is robust to the choice of prior. This measurement provides a constraint on models of the prompt $γ$-ray and radio/X-ray afterglow emission associated with the merger; for example, it is consistent with the off-axis viewing angle inferred for a structured jet model. We provide for the first time the full posterior samples from Bayesian parameter estimation of LIGO/Virgo data to enable further analysis by the community. △ Less

Submitted 8 June, 2018; v1 submitted 11 April, 2018; originally announced April 2018.

Comments: 7 pages, 3 figures

Journal ref: D. Finstad et al. 2018 ApJL 860 L2

A Large Scale Spatio-temporal Binomial Regression Model for Estimating Seroprevalence Trends

Authors: Stella Watson Self, Christopher McMahan, D. Andrew Brown, Robert Lund, Jenna Gettings, Michael Yabsley

Abstract: This paper develops a large-scale Bayesian spatio-temporal binomial regression model for the purpose of investigating regional trends in antibody prevalence to Borrelia burgdorferi, the causative agent of Lyme disease. The proposed model uses Gaussian predictive processes to estimate the spatially varying trends and a conditional autoregressive model to account for spatio-temporal dependence. Care… ▽ More This paper develops a large-scale Bayesian spatio-temporal binomial regression model for the purpose of investigating regional trends in antibody prevalence to Borrelia burgdorferi, the causative agent of Lyme disease. The proposed model uses Gaussian predictive processes to estimate the spatially varying trends and a conditional autoregressive model to account for spatio-temporal dependence. Careful consideration is made to develop a novel framework that is scalable to large spatio-temporal data. The proposed model is used to analyze approximately 16 million Borrelia burgdorferi test results collected on dogs located throughout the conterminous United States over a sixty month period. This analysis identifies several regions of increasing canine risk. Specifically, this analysis reveals evidence that Lyme disease is getting worse in some endemic regions and that it could potentially be spreading to other non-endemic areas. Further, given the zoonotic nature of this vector-borne disease, this analysis could potentially reveal areas of increasing human risk. △ Less

Submitted 29 March, 2018; originally announced March 2018.

Comments: 19 pages without figures. All figures are available as ancillary files

A Search for Tensor, Vector, and Scalar Polarizations in the Stochastic Gravitational-Wave Background

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Afrough, B. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, G. Allen, A. Allocca , et al. (1075 additional authors not shown)

Abstract: The detection of gravitational waves with Advanced LIGO and Advanced Virgo has enabled novel tests of general relativity, including direct study of the polarization of gravitational waves. While general relativity allows for only two tensor gravitational-wave polarizations, general metric theories can additionally predict two vector and two scalar polarizations. The polarization of gravitational w… ▽ More The detection of gravitational waves with Advanced LIGO and Advanced Virgo has enabled novel tests of general relativity, including direct study of the polarization of gravitational waves. While general relativity allows for only two tensor gravitational-wave polarizations, general metric theories can additionally predict two vector and two scalar polarizations. The polarization of gravitational waves is encoded in the spectral shape of the stochastic gravitational-wave background, formed by the superposition of cosmological and individually-unresolved astrophysical sources. Using data recorded by Advanced LIGO during its first observing run, we search for a stochastic background of generically-polarized gravitational waves. We find no evidence for a background of any polarization, and place the first direct bounds on the contributions of vector and scalar polarizations to the stochastic background. Under log-uniform priors for the energy in each polarization, we limit the energy-densities of tensor, vector, and scalar modes at 95% credibility to $Ω^T_0 < 5.6 \times 10^{-8}$, $Ω^V_0 < 6.4\times 10^{-8}$, and $Ω^S_0 < 1.1\times 10^{-7}$ at a reference frequency $f_0 = 25$ Hz. △ Less

Submitted 2 October, 2019; v1 submitted 27 February, 2018; originally announced February 2018.

Comments: Minor updates to match published version

Journal ref: Phys. Rev. Lett. 120, 201102 (2018)

Full Band All-sky Search for Periodic Gravitational Waves in the O1 LIGO Data

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Afrough, B. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, B. Allen, G. Allen, A. Allocca, P. A. Altin , et al. (1077 additional authors not shown)

Abstract: We report on a new all-sky search for periodic gravitational waves in the frequency band 475-2000 Hz and with a frequency time derivative in the range of [-1.0e-8, +1e-9] Hz/s. Potential signals could be produced by a nearby spinning and slightly non-axisymmetric isolated neutron star in our galaxy. This search uses the data from Advanced LIGO's first observational run O1. No gravitational wave… ▽ More We report on a new all-sky search for periodic gravitational waves in the frequency band 475-2000 Hz and with a frequency time derivative in the range of [-1.0e-8, +1e-9] Hz/s. Potential signals could be produced by a nearby spinning and slightly non-axisymmetric isolated neutron star in our galaxy. This search uses the data from Advanced LIGO's first observational run O1. No gravitational wave signals were observed, and upper limits were placed on their strengths. For completeness, results from the separately published low frequency search 20-475 Hz are included as well. Our lowest upper limit on worst-case (linearly polarized) strain amplitude h_0 is 4e-25 near 170 Hz, while at the high end of our frequency range we achieve a worst-case upper limit of 1.3e-24. For a circularly polarized source (most favorable orientation), the smallest upper limit obtained is ~1.5e-25. △ Less

Submitted 14 February, 2018; originally announced February 2018.

Journal ref: Phys. Rev. D 97, 102003 (2018)

Constraints on cosmic strings using data from the first Advanced LIGO observing run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Afrough, B. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, G. Allen, A. Allocca , et al. (1020 additional authors not shown)

Abstract: Cosmic strings are topological defects which can be formed in GUT-scale phase transitions in the early universe. They are also predicted to form in the context of string theory. The main mechanism for a network of Nambu-Goto cosmic strings to lose energy is through the production of loops and the subsequent emission of gravitational waves, thus offering an experimental signature for the existence… ▽ More Cosmic strings are topological defects which can be formed in GUT-scale phase transitions in the early universe. They are also predicted to form in the context of string theory. The main mechanism for a network of Nambu-Goto cosmic strings to lose energy is through the production of loops and the subsequent emission of gravitational waves, thus offering an experimental signature for the existence of cosmic strings. Here we report on the analysis conducted to specifically search for gravitational-wave bursts from cosmic string loops in the data of Advanced LIGO 2015-2016 observing run (O1). No evidence of such signals was found in the data, and as a result we set upper limits on the cosmic string parameters for three recent loop distribution models. In this paper, we initially derive constraints on the string tension $Gμ$ and the intercommutation probability, using not only the burst analysis performed on the O1 data set, but also results from the previously published LIGO stochastic O1 analysis, pulsar timing arrays, cosmic microwave background and Big-Bang nucleosynthesis experiments. We show that these data sets are complementary in that they probe gravitational waves produced by cosmic string loops during very different epochs. Finally, we show that the data sets exclude large parts of the parameter space of the three loop distribution models we consider. △ Less

Submitted 2 May, 2018; v1 submitted 4 December, 2017; originally announced December 2017.

Comments: Physical Review D, in-press

Journal ref: Phys. Rev. D 97, 102002 (2018)