skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Efficient Stochastic Template Bank Using Inner Product Inequalities
Abstract Gravitational wave searches are crucial for studying compact sources such as neutron stars and black holes. Many sensitive modeled searches use matched filtering to compare gravitational strain data to a set of waveform models known as template banks. We introduce a new stochastic placement method for constructing template banks, offering efficiency and flexibility to handle arbitrary parameter spaces, including orbital eccentricity, tidal deformability, and other extrinsic parameters. This method can be computationally limited by the ability to compare proposal templates with the accepted templates in the bank. To alleviate this computational load, we introduce the use of inner product inequalities to reduce the number of required comparisons. We also introduce a novel application of Gaussian Kernel Density Estimation to enhance waveform coverage in sparser regions. Our approach has been employed to search for eccentric binary neutron stars, low-mass neutron stars, primordial black holes, and supermassive black hole binaries. We demonstrate that our method produces self-consistent banks that recover the required minimum fraction of signals. For common parameter spaces, our method shows comparable computational performance and similar template bank sizes to geometric placement methods and stochastic methods, while easily extending to higher-dimensional problems. The time to run a search exceeds the time to generate the bank by a factor of O ( 10 5 ) for dedicated template banks, such as geometric, mass-only stochastic, and aligned spin cases, O ( 10 4 ) for eccentric and O ( 10 3 ) for the tidal deformable bank. With the advent of efficient template bank generation, the primary area for improvement is developing more efficient search methodologies.  more » « less
Award ID(s):
2309240
PAR ID:
10580402
Author(s) / Creator(s):
; ; ; ; ; ; ;
Publisher / Repository:
The Astrophysical Journal
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
975
Issue:
2
ISSN:
0004-637X
Page Range / eLocation ID:
212
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The Simulating eXtreme Spacetimes Collaboration’s codeSpECcan now routinely simulate binary black hole mergers undergoing 25 orbits, with the longest simulations undergoing nearly 180 orbits. While this sounds impressive, the mismatch between the highest resolutions for this long simulation is O ( 10 1 ) . Meanwhile, the mismatch between resolutions for the more typical simulations tends to be O ( 10 4 ) , despite the resolutions being similar to the long simulations’. In this note, we explain why mismatch alone gives an incomplete picture of code—and waveform—quality, especially in the context of providing waveform templates for LISA and 3G detectors, which require templates with O ( 10 3 ) O ( 10 5 ) orbits. We argue that to ready the GW community for the sensitivity of future detectors, numerical relativity groups must be aware of this caveat, and also run future simulations with at least three resolutions to properly assess waveform accuracy. 
    more » « less
  2. Abstract A test of lepton flavor universality in B ± K ± μ + μ and B ± K ± e + e decays, as well as a measurement of differential and integrated branching fractions of a nonresonant B ± K ± μ + μ decay are presented. The analysis is made possible by a dedicated data set of proton-proton collisions at s = 13 TeV recorded in 2018, by the CMS experiment at the LHC, using a special high-rate data stream designed for collecting about 10 billion unbiased b hadron decays. The ratio of the branching fractions B ( B ± K ± μ + μ ) to B ( B ± K ± e + e ) is determined from the measured double ratio R ( K ) of these decays to the respective branching fractions of the B ± J / ψ K ± with J / ψ μ + μ and e + e decays, which allow for significant cancellation of systematic uncertainties. The ratio R ( K ) is measured in the range 1.1 < q 2 < 6.0 GeV 2 , whereqis the invariant mass of the lepton pair, and is found to be R ( K ) = 0.78 0.23 + 0.47 , in agreement with the standard model expectation R ( K ) 1 . This measurement is limited by the statistical precision of the electron channel. The integrated branching fraction in the sameq2range, B ( B ± K ± μ + μ ) = ( 12.42 ± 0.68 ) × 10 8 , is consistent with the present world-average value and has a comparable precision. 
    more » « less
  3. Abstract JWST is revealing a remarkable new population of high-redshift (z ≳ 4), low-luminosity active galactic nuclei in deep surveys and detecting the host galaxy's stellar light in the most luminous and massive quasars atz ∼ 6 for the first time. Recent findings claim that supermassive black holes (SMBHs) in these systems are significantly more massive than predicted by the local black hole (BH) mass–stellar mass ( M BH - M ) relation and that this is not due to sample selection effects. Through detailed statistical modeling, we demonstrate that the coupled effects of selection biases (i.e., finite detection limit and requirements for detecting broad lines) and measurement uncertainties can largely explain the reported offset and flattening in the observed M BH - M relation toward the upper envelope of the local relation, even for those at M BH < 1 0 8 M . We further investigate the possible evolution of the M BH - M relation atz ≳ 4 with careful treatment of observational biases and consideration of the degeneracy between intrinsic evolution and dispersion in this relation. The bias-corrected intrinsic M BH - M relation in the low-mass regime ( M 1 0 10 M ) suggests a large population of low-mass BHs ( M BH 1 0 5 M ), possibly originating from lighter seeds, may remain undetected or unidentified. These results underscore the importance of forward modeling observational biases to better understand BH seeding and SMBH–galaxy coevolution mechanisms in the early universe, even with the deepest JWST surveys. 
    more » « less
  4. Abstract Analytic continuation from (3, 1) signature Minkowski to (2, 2) signature Klein space has emerged as a useful tool for the understanding of scattering amplitudes and flat space holography. Under this continuation, past and future null infinity merge into a single boundary ( J ) which is the product of a null line with a (1, 1) signature torus. The Minkowskian S -matrix continues to a Kleinian S -vector which in turn may be represented by a Poincaré-invariant vacuum state | C in the Hilbert space built on J . | C contains all information about S in a novel, repackaged form. We give an explicit construction of | C in a Lorentz/conformal basis for a free massless scalar. J separates into two halves J ± which are the asymptotic null boundaries of the regions timelike and spacelike separated from the origin. | C is shown to be a maximally entangled state in the product of the J ± Hilbert spaces. 
    more » « less
  5. Abstract We present a new suite of numerical simulations of the star-forming interstellar medium (ISM) in galactic disks using the TIGRESS-NCR framework. Distinctive aspects of our simulation suite are (1) sophisticated and comprehensive numerical treatments of essential physical processes including magnetohydrodynamics, self-gravity, and galactic differential rotation, as well as photochemistry, cooling, and heating coupled with direct ray-tracing UV radiation transfer and resolved supernova feedback and (2) wide parameter coverage including the variation in metallicity over Z Z / Z 0.1 - 3 , gas surface density Σgas∼ 5–150Mpc−2, and stellar surface density Σstar∼ 1–50Mpc−2. The range of emergent star formation rate surface density is ΣSFR∼ 10−4–0.5Mkpc−2yr−1, and ISM total midplane pressure isPtot/kB= 103–106cm−3K, withPtotequal to the ISM weight W . For given Σgasand Σstar, we find Σ SFR Z 0.3 . We provide an interpretation based on the pressure-regulated feedback-modulated (PRFM) star formation theory. The total midplane pressure consists of thermal, turbulent, and magnetic stresses. We characterize feedback modulation in terms of the yield ϒ, defined as the ratio of each stress to ΣSFR. The thermal feedback yield varies sensitively with both weight and metallicity as ϒ th W 0.46 Z 0.53 , while the combined turbulent and magnetic feedback yield shows weaker dependence ϒ turb + mag W 0.22 Z 0.18 . The reduction in ΣSFRat low metallicity is due mainly to enhanced thermal feedback yield, resulting from reduced attenuation of UV radiation. With the metallicity-dependent calibrations we provide, PRFM theory can be used for a new subgrid star formation prescription in cosmological simulations where the ISM is unresolved. 
    more » « less