Overture: An advanced object-oriented software system for moving overlapping grid computations

While the development of high-level, easy-to-use, software libraries for numerical computations has been successful in some areas (e.g. linear system solvers, ODE solvers, grid generation), this has been an elusive goal for developers of partial differential equation (PDE) solvers. The advent of new high level languages such as C++ has begun to make this an achievable goal. This report discusses an object- oriented environment that we are developing for solving problems on overlapping (Chimera) grids. The goal of this effort is to support flexible PDE solvers on adaptive, moving, overlapping grids that cover a domain and overlap where they meet. Solutions values at the overlap are determined by interpolation. The overlapping grid approach is particularly efficient for rapidly generating high- quality grids for moving geometries since as the component grids move, only the list of interpolation points changes, and the component grids do not have to be regenerated. We use structured component grids so that efficient, fast finite-difference algorithms can be used. Oliger-Berger-Corella type mesh refinement is used to efficiently resolve fine features of the flow

Overture: an object-oriented software system for solving partial differential equations in serial and parallel environments

The OVERTURE Framework is an object-oriented environment for solving PDEs on serial and parallel architectures. It is a collection of C++ libraries that enables the use of finite difference and finite volume methods at a level that hides the details of the associated data structures, as well as the details of the parallel implementation. It is based on the A++/P++ array class library and is designed for solving problems on a structured grid or a collection of structured grids. In particular, it can use curvilinear grids, adaptive mesh refinement and the composite overlapping grid method to represent problems with complex moving geometry

‘The Brick’ is not a brick: a comprehensive study of the structure and dynamics of the central molecular zone cloud G0.253+0.016

This article has been accepted for publication in Monthly Notices of the Royal Astronomical Society © 2019 The Author(s). Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reserved.In this paper we provide a comprehensive description of the internal dynamics of G0.253+0.016 (a.k.a. ‘the Brick’); one of the most massive and dense molecular clouds in the Galaxy to lack signatures of widespread star formation. As a potential host to a future generation of high-mass stars, understanding largely quiescent molecular clouds like G0.253+0.016 is of critical importance. In this paper, we reanalyse Atacama Large Millimeter Array cycle 0 HNCO J = 4(0, 4) − 3(0, 3) data at 3 mm, using two new pieces of software that we make available to the community. First, SCOUSEPY, a Python implementation of the spectral line fitting algorithm SCOUSE. Secondly, ACORNS (Agglomerative Clustering for ORganising Nested Structures), a hierarchical n-dimensional clustering algorithm designed for use with discrete spectroscopic data. Together, these tools provide an unbiased measurement of the line-of-sight velocity dispersion in this cloud, σvlos,1D=4.4±2.1 km s−1, which is somewhat larger than predicted by velocity dispersion-size relations for the central molecular zone (CMZ). The dispersion of centroid velocities in the plane of the sky are comparable, yielding σvlos,1D/σvpos,1D∼1.2±0.3⁠. This isotropy may indicate that the line-of-sight extent of the cloud is approximately equivalent to that in the plane of the sky. Combining our kinematic decomposition with radiative transfer modelling, we conclude that G0.253+0.016 is not a single, coherent, and centrally condensed molecular cloud; ‘the Brick’ is not a brick. Instead, G0.253+0.016 is a dynamically complex and hierarchically structured molecular cloud whose morphology is consistent with the influence of the orbital dynamics and shear in the CMZ

Star formation in 'the Brick': ALMA reveals an active protocluster in the Galactic centre cloud G0.253+0.016

Interstellar matter and star formatio

Star clusters near and far; tracing star formation across cosmic time

© 2020 Springer-Verlag. The final publication is available at Springer via https://doi.org/10.1007/s11214-020-00690-x.Star clusters are fundamental units of stellar feedback and unique tracers of their host galactic properties. In this review, we will first focus on their constituents, i.e.\ detailed insight into their stellar populations and their surrounding ionised, warm, neutral, and molecular gas. We, then, move beyond the Local Group to review star cluster populations at various evolutionary stages, and in diverse galactic environmental conditions accessible in the local Universe. At high redshift, where conditions for cluster formation and evolution are more extreme, we are only able to observe the integrated light of a handful of objects that we believe will become globular clusters. We therefore discuss how numerical and analytical methods, informed by the observed properties of cluster populations in the local Universe, are used to develop sophisticated simulations potentially capable of disentangling the genetic map of galaxy formation and assembly that is carried by globular cluster populations.Peer reviewedFinal Accepted Versio

Search for gravitational-wave transients associated with magnetar bursts in advanced LIGO and advanced Virgo data from the third observing run

Gravitational waves are expected to be produced from neutron star oscillations associated with magnetar giant f lares and short bursts. We present the results of a search for short-duration (milliseconds to seconds) and longduration (∼100 s) transient gravitational waves from 13 magnetar short bursts observed during Advanced LIGO, Advanced Virgo, and KAGRA’s third observation run. These 13 bursts come from two magnetars, SGR1935 +2154 and SwiftJ1818.0−1607. We also include three other electromagnetic burst events detected by FermiGBM which were identified as likely coming from one or more magnetars, but they have no association with a known magnetar. No magnetar giant flares were detected during the analysis period. We find no evidence of gravitational waves associated with any of these 16 bursts. We place upper limits on the rms of the integrated incident gravitational-wave strain that reach 3.6 × 10−²³ Hz at 100 Hz for the short-duration search and 1.1 ×10−²² Hz at 450 Hz for the long-duration search. For a ringdown signal at 1590 Hz targeted by the short-duration search the limit is set to 2.3 × 10−²² Hz. Using the estimated distance to each magnetar, we derive upper limits upper limits on the emitted gravitational-wave energy of 1.5 × 1044 erg (1.0 × 1044 erg) for SGR 1935+2154 and 9.4 × 10^43 erg (1.3 × 1044 erg) for Swift J1818.0−1607, for the short-duration (long-duration) search. Assuming isotropic emission of electromagnetic radiation of the burst fluences, we constrain the ratio of gravitational-wave energy to electromagnetic energy for bursts from SGR 1935+2154 with the available fluence information. The lowest of these ratios is 4.5 × 103

A joint Fermi-GBM and Swift-BAT analysis of gravitational-wave candidates from the third gravitational-wave observing run

We present Fermi Gamma-ray Burst Monitor (Fermi-GBM) and Swift Burst Alert Telescope (Swift-BAT) searches for gamma-ray/X-ray counterparts to gravitational-wave (GW) candidate events identified during the third observing run of the Advanced LIGO and Advanced Virgo detectors. Using Fermi-GBM onboard triggers and subthreshold gamma-ray burst (GRB) candidates found in the Fermi-GBM ground analyses, the Targeted Search and the Untargeted Search, we investigate whether there are any coincident GRBs associated with the GWs. We also search the Swift-BAT rate data around the GW times to determine whether a GRB counterpart is present. No counterparts are found. Using both the Fermi-GBM Targeted Search and the Swift-BAT search, we calculate flux upper limits and present joint upper limits on the gamma-ray luminosity of each GW. Given these limits, we constrain theoretical models for the emission of gamma rays from binary black hole mergers

Constraints on the cosmic expansion history from GWTC–3

We use 47 gravitational wave sources from the Third LIGO–Virgo–Kamioka Gravitational Wave Detector Gravitational Wave Transient Catalog (GWTC–3) to estimate the Hubble parameter H(z), including its current value, the Hubble constant H0. Each gravitational wave (GW) signal provides the luminosity distance to the source, and we estimate the corresponding redshift using two methods: the redshifted masses and a galaxy catalog. Using the binary black hole (BBH) redshifted masses, we simultaneously infer the source mass distribution and H(z). The source mass distribution displays a peak around 34 M⊙, followed by a drop-off. Assuming this mass scale does not evolve with the redshift results in a H(z) measurement, yielding ${H}_{0}={68}_{-8}^{+12}\,\mathrm{km}\ \,\ {{\rm{s}}}^{-1}\,{\mathrm{Mpc}}^{-1}$ (68% credible interval) when combined with the H0 measurement from GW170817 and its electromagnetic counterpart. This represents an improvement of 17% with respect to the H0 estimate from GWTC–1. The second method associates each GW event with its probable host galaxy in the catalog GLADE+, statistically marginalizing over the redshifts of each event's potential hosts. Assuming a fixed BBH population, we estimate a value of ${H}_{0}={68}_{-6}^{+8}\,\mathrm{km}\ \,\ {{\rm{s}}}^{-1}\,{\mathrm{Mpc}}^{-1}$ with the galaxy catalog method, an improvement of 42% with respect to our GWTC–1 result and 20% with respect to recent H0 studies using GWTC–2 events. However, we show that this result is strongly impacted by assumptions about the BBH source mass distribution; the only event which is not strongly impacted by such assumptions (and is thus informative about H0) is the well-localized event GW190814

Open data from the third observing run of LIGO, Virgo, KAGRA, and GEO

The global network of gravitational-wave observatories now includes five detectors, namely LIGO Hanford, LIGO Livingston, Virgo, KAGRA, and GEO 600. These detectors collected data during their third observing run, O3, composed of three phases: O3a starting in 2019 April and lasting six months, O3b starting in 2019 November and lasting five months, and O3GK starting in 2020 April and lasting two weeks. In this paper we describe these data and various other science products that can be freely accessed through the Gravitational Wave Open Science Center at https://gwosc.org. The main data set, consisting of the gravitational-wave strain time series that contains the astrophysical signals, is released together with supporting data useful for their analysis and documentation, tutorials, as well as analysis software packages

Search for eccentric black hole coalescences during the third observing run of LIGO and Virgo

Despite the growing number of binary black hole coalescences confidently observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include the effects of eccentricity. Here, we present observational results for a waveform-independent search sensitive to eccentric black hole coalescences, covering the third observing run (O3) of the LIGO and Virgo detectors. We identified no new high-significance candidates beyond those that have already been identified with searches focusing on quasi-circular binaries. We determine the sensitivity of our search to high-mass (total source-frame mass M > 70 M⊙) binaries covering eccentricities up to 0.3 at 15 Hz emitted gravitational-wave frequency, and use this to compare model predictions to search results. Assuming all detections are indeed quasi-circular, for our fiducial population model, we place a conservative upper limit for the merger rate density of high-mass binaries with eccentricities 0 < e ≤ 0.3 at 16.9 Gpc−3 yr−1 at the 90% confidence level