saprEMo: a simplified algorithm for predicting detections of electromagnetic transients in surveys

The multi-wavelength detection of GW170817 has inaugurated multi-messenger astronomy. The next step consists in interpreting observations coming from population of gravitational wave sources. We introduce saprEMo, a tool aimed at predicting the number of electromagnetic signals characterised by a specific light curve and spectrum, expected in a particular sky survey. By looking at past surveys, saprEMo allows us to constrain models of electromagnetic emission or event rates. Applying saprEMo to proposed astronomical missions/observing campaigns provides a perspective on their scientific impact and tests the effect of adopting different observational strategies. For our first case study, we adopt a model of spindown-powered X-ray emission predicted for a binary neutron star merger producing a long-lived neutron star. We apply saprEMo on data collected by XMM-Newton and Chandra and during $10^4$ s of observations with the mission concept THESEUS. We demonstrate that our emission model and binary neutron star merger rate imply the presence of some signals in the XMM-Newton catalogs. We also show that the new class of X-ray transients found by Bauer et al. in the Chandra Deep Field-South is marginally consistent with the expected rate. Finally, by studying the mission concept THESEUS, we demonstrate the substantial impact of a much larger field of view in searches of X-ray transients

Exploring the Parameter Space of Compact Binary Population Synthesis

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

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

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

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

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

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

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

Maximizing efficiency of rumen microbial protein production.

Rumen microbes produce cellular protein inefficiently partly because they do not direct all ATP toward growth. They direct some ATP toward maintenance functions, as long-recognized, but they also direct ATP toward reserve carbohydrate synthesis and energy spilling (futile cycles that dissipate heat). Rumen microbes expend ATP by vacillating between (1) accumulation of reserve carbohydrate after feeding (during carbohydrate excess) and (2) mobilization of that carbohydrate thereafter (during carbohydrate limitation). Protozoa account for most accumulation of reserve carbohydrate, and in competition experiments, protozoa accumulated nearly 35-fold more reserve carbohydrate than bacteria. Some pure cultures of bacteria spill energy, but only recently have mixed rumen communities been recognized as capable of the same. When these communities were dosed glucose in vitro, energy spilling could account for nearly 40% of heat production. We suspect that cycling of glycogen (a major reserve carbohydrate) is a major mechanism of spilling; such cycling has already been observed in single-species cultures of protozoa and bacteria. Interconversions of short-chain fatty acids (SCFA) may also expend ATP and depress efficiency of microbial protein production. These interconversions may involve extensive cycling of intermediates, such as cycling of acetate during butyrate production in certain butyrivibrios. We speculate this cycling may expend ATP directly or indirectly. By further quantifying the impact of reserve carbohydrate accumulation, energy spilling, and SCFA interconversions on growth efficiency, we can improve prediction of microbial protein production and guide efforts to improve efficiency of microbial protein production in the rumen

The formation of compact object binaries through isolated binary evolution

Observations indicate that most stars are in binary or higher multiplicity systems (Preibisch et al., 1999; Sana et al., 2012, 2013; Duchêne and Kraus, 2013; Chini et al., 2013; Sota et al., 2014; Kobulnicky et al., 2014; Dunstall et al., 2015; Moe and Di Stefano, 2017; Sana, 2017). In a binary two stars orbit each other, bound by their mutual gravitational pull. If the orbital separation is short enough the stars interact, drastically altering their evolution. In the past decades, the means to perform complex calculations have drastically improved, giving us the chance to explore the physics of binary evolution in greater detail. This helped explain several observed properties, amongst others, why some stars are more luminous than expected or have peculiar surface abundances. Nonetheless, large uncertainties persist in the field of binary star physics. In 2015, the gravitational waves from two colliding black holes were detected for the first time (Abbott et al., 2016b). For decades it has been hypothesised that two massive stars in an isolated binary could interact without external influences and form a binary black hole system tight enough to merge in a Hubble time (van den Heuvel and De Loore, 1973; Tutukov and Yungelson, 1973). In this dissertation I follow in the footsteps of many other studies and assume the observed gravitational-wave events come from isolated binary evolution, even though other formation channels for the mergers of neutron stars and black holes are also possible. The aim is to study what constraints, if any, properties of gravitational-wave events can place on the evolution of massive stars in binaries. The general approach in of this dissertation is to evolve a population of stars under various model assumptions and estimate the rates and properties of gravitational wave mergers for each model. The predicted distributions enable a quantitative or qualitative assessment of the impact of current uncertainties in binary-star physics on estimates of the rates and masses of gravitational-wave events. Evaluating the effect of uncertainties is crucial to determine whether comparisons between synthetic populations of gravitational-wave sources and observations can place meaningful constraints on binary-star physics. If the uncertainties are large and model-dependent features are not predicted, then the detections of gravitational-wave mergers may only provide marginal constraints. In this dissertation I assess the impact of the following model assumptions. In chapter 4 I investigate the uncertainties in the rate and initial chemical composition with which stars form and the impact of these uncertainties on the predictions of the merger rate of neutron stars and black holes. In chapter 5 I vary the response of stars to mass loss and explore how it alters the interactions that lead to the formation of binary black holes. In chapter 6 I examine whether Cygnus X-1 may evolve into a binary black hole system and if the observed mass of the black hole in Cygnus X-1 provides constraints on the wind mass-loss rates of stars. Chapter 1 and chapter 2 provide introductory material for the reader on stellar evolution and binary interactions. Chapter 3 summarises the theoretical model used in this dissertation to evolve a population of stars. Chapter 7 provides a summary and personal view on the conclusions of this dissertation

The role of magnesium and calcium ions in the glucose dehydrogenase activity of <i>Klebsiella pneumoniae</i> NCTC 418

Magnesium-limited chemostat cultures of Klebsiella pneumoniae NCTC 418 with 20 μM CaCl2 in the medium showed a low rate of gluconate plus 2-ketogluconate production relative to potassium- or phosphate-limited cultures. However, when the medium concentration of CaCl2 was increased to 1 mM, the glucose dehydrogenase (GDH) activities also increased and became similar to those observed in potassium- or phosphate limited cultures. It is concluded that this is due to Mg2+ and Ca2+ ions being involved in the binding of pyrroloquinoline quinone (PQQ) to the GDH apoenzyme. There seems to be an absolute requirement of divalent cations for proper enzyme functioning and in this respect Ca2+ ions could replace Mg2+ ions. The high GDH activity which has been found in cells grown under Mg2−-limited conditions in the presence of higher concentrations of Ca2+ ions, is compatible with the earlier proposal that GDH functions as an auxiliary energy generating system involved in the maintenance of high transmembrane ion gradients.Centro de Investigación y Desarrollo en Fermentaciones Industriale