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

A deep learning approach to halo merger tree construction

A key ingredient for semi-analytic models of galaxy formation is the mass assembly history of haloes, encoded in a tree structure. The most commonly used method to construct halo merger histories is based on the outcomes of high-resolution, computationally intensive N-body simulations. We show that machine learning (ML) techniques, in particular Generative Adversarial Networks (GANs), are a promising new tool to tackle this problem with a modest computational cost and retaining the best features of merger trees from simulations. We train our GAN model with a limited sample of merger trees from the Evolution and Assembly of GaLaxies and their Environments (EAGLE) simulation suite, constructed using two halo finders-tree builder algorithms: SUBFIND-D-TREES and ROCKSTAR-ConsistentTrees. Our GAN model successfully learns to generate well-constructed merger tree structures with high temporal resolution, and to reproduce the statistical features of the sample of merger trees used for training, when considering up to three variables in the training process. These inputs, whose representations are also learned by our GAN model, are mass of the halo progenitors and the final descendant, progenitor type (main halo or satellite), and distance of a progenitor to that in the main branch. The inclusion of the latter two inputs greatly improves the final learned representation of the halo mass growth history, especially for SUBFIND-like ML trees. When comparing equally sized samples of ML merger trees with those of the EAGLE simulation, we find better agreement for SUBFIND-like ML trees. Finally, our GAN-based framework can be utilized to construct merger histories of low-and intermediate-mass haloes, the most abundant in cosmological simulations.Fil: Robles, Sandra. Universidad Autónoma de Madrid; España. Kings College London (kcl); . University of Melbourne; AustraliaFil: Gómez, Jonathan S. Universidad Católica de Chile; Chile. Universidad Autónoma de Madrid; España. Pontificia Universidad Católica de Chile; ChileFil: Ramírez Rivera, Adín. University of Oslo; NoruegaFil: Padilla, Nelson David. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Instituto de Astronomía Teórica y Experimental. Universidad Nacional de Córdoba. Observatorio Astronómico de Córdoba. Instituto de Astronomía Teórica y Experimental; ArgentinaFil: Dujovne, Diego. Universidad Diego Portales; Chil

Development of Antibacterial Steel Surfaces Through Laser Texturing

[Abstract] The aim of the present study was to develop novel antibacterial touch surfaces through the laser texturing optimization of stainless steel. A wide range of laser fluence (2.11 J/cm2–5.64 J/cm2) and scanning interval (10 µm–30 μm) parameters were explored. The impact of surfaces with different patterns, wettability, and oxidation states on the antimicrobial behavior of Escherichia coli K-12 and biofilm hyper-producing Acinetobacter baumannii MAR002 was assessed. Modification of laser input enacted topographical changes with high scanning intervals leading to ordered surface patterns, while increasing the laser fluence to 5.64 J/cm2 created larger and less ordered plateaus and valleys. Texturing also drove a transition from a hydrophilic starting surface with a contact angle of 80.67° ± 3.35° to hydrophobic (138°–148°). Antimicrobial analysis and bioluminescence assays of E. coli, alongside biofilm forming test through A. baumannii MAR002 indicated the ability of laser texturing to produce effective bactericidal touch surfaces. No simple correlation was found between wettability and bacterial behavior, revealing that proliferation is dependent on roughness, oxidation, and wettability. For the conditions selected in this study, a laser fluence of 5.64 J/cm2 and a scanning interval of 10 µm showcased the lowest amount of recovered bacteria after 30 min.This research was supported by Project Nos. 592 p-01216A and IJCI-2016-29524 (awarded to A.P.G.), funded by the Spanish Society of Infectious Diseases and Clinical Microbiology (SEIMC) and MINECO, respectively. This manuscript is part of Process Design to Prevent Prosthetic Infections (Grant No. EP/P02341X/1)Sociedad Española de Enfermedades Infecciosas y Microbiología Clínica; 592 p-01216AReino Unido. Engineering and Physical Sciences Research Council; EP/P02341X/

Halo merger tree comparison: impact on galaxy formation models

We examine the effect of using different halo finders and merger tree building algorithms on galaxy properties predicted using the GALFORM semi-analytical model run on a high resolution, large volume dark matter simulation. The halo finders/tree builders HBT, ROCKSTAR, SUBFIND, and VELOCI RAPTOR differ in their definitions of halo mass, on whether only spatial or phase-space information is used, and in how they distinguish satellite and main haloes; all of these features have some impact on the model galaxies, even after the trees are post-processed and homogenized by GALFORM. The stellar mass function is insensitive to the halo and merger tree finder adopted. However, we find that the number of central and satellite galaxies in GALFORM does depend slightly on the halo finder/tree builder. The number of galaxies without resolved subhaloes depends strongly on the tree builder, with VELOCIRAPTOR, a phase-space finder, showing the largest population of such galaxies. The distributions of stellar masses, cold and hot gas masses, and star formation rates agree well between different halo finders/tree builders. However, because VELOCIRAPTOR has more early progenitor haloes, with these trees GALFORM produces slightly higher star formation rate densities at high redshift, smaller galaxy sizes, and larger stellar masses for the spheroid component. Since in all cases these differences are small we conclude that, when all of the trees are processed so that the main progenitor mass increases monotonically, the predicted GALFORM galaxy populations are stable and consistent for these four halo finders/tree builders

How birds cope physiologically and behaviourally with extreme climatic events

As global climate change progresses, the occurrence of potentially disruptiveclimatic events such as storms are increasing in frequency, duration and inten-sity resulting in higher mortality and reduced reproductive success. Whatconstitutes an extreme climatic event? First we point out that extreme climaticevents in biological contexts can occur in any environment. Focusing on fieldand laboratory data on wild birds we propose a mechanistic approach to defin-ing and investigating what extreme climatic events are and how animals copewith them at physiological and behavioural levels. The life cycle of birds ismade up of life-history stages such as migration, breeding and moult thatevolved to match a range of environmental conditions an individual mightexpect during the year. When environmental conditions deteriorate anddeviate from the expected range then the individual must trigger copingmechanisms (emergency life-history stage) that will disrupt the temporal pro-gression of life-history stages, but enhance survival. Using the framework ofallostasis, we argue that an extreme climatic event in biological contexts canbe defined as when the cumulative resources available to an individual areexceeded by the sum of its energetic costs—a state called allostatic overload.This allostatic overload triggers the emergency life-history stage that tempor-arily allows the individual to cease regular activities in an attempt to surviveextreme conditions. We propose that glucocorticoid hormones play a majorrole in orchestrating coping mechanisms and are critical for enduring extremeclimatic events.This article is part of the themed issue ‘Behavioural, ecological andevolutionary responses to extreme climatic events’

Genetic analysis of 17 Y-STRs in a Mestizo population from the Central Valley of Mexico

This study aims to portray the complex diversity of the Mexican Mestizo population, which represents 98.8% of the entire population of Mexico. We compiled extended haplotype data of the Y chromosome from populations in the Central Valley of Mexico (CVM), which were compared to other Mestizo and parental (Amerindian, European and African) populations. A complex ancestral relationship was found in the CVM population, suggesting cosmopolitan origins. Nevertheless, the most preeminent lineages point towards a European ancestry, where the R1b was the most frequent. In addition, important frequencies of Amerindian linages were also found in the Mestizo sample studied. Interestingly, the Amerindian ancestry showed a remarkable substructure, which was represented by the two main founding lineages: QL54 (x M3) and M3. However, even within each lineage a high diversity was found despite the small number of samples bearers of these lineages. Further, we detected important genetic differences between the CVM populations and the Mexican Mestizo populations from the north and south. This result points to the fact that Mestizo populations present different ancestral proportions, which are related to the demographic events that gave origin to each population. Finally, we provide additional forensic statistical parameters that are useful in the interpretation of genetic analysis where autosomal loci are limited. Our findings illustrate the complex genetic background of the Mexican Mestizo population and reinforce the need to encompass more geographic regions to generate more robust data for forensic applications

Immunosuppressive Minimization Strategies in Kidney Transplantation

The long-term graft survival in renal transplantation results is still controversial, the toxicity and adverse reactions of the immunosuppressive drugs are implicated, as well as cellular and humoral antigen-specific immune mechanisms; therefore, different strategies for adapting immunosuppression are used to reduce the complications associated with the use of these drugs. Calcineurin inhibitors (CNI) require an adequate dose-dependent concentration leading to the appearance of drug-related adverse reactions. The variability in the required dose of CNI leads to minimization strategies that do not result in a higher acute rejection (AR) incidence when compared to other immunosuppressive agents. Early steroid withdrawal is another strategy, although with an increase in AR, but without an impact on the function and survival of the renal graft. The reduction of mycophenolate mofetil to 1.5 g/day seems to be a therapeutic option, decreasing the infectious, hematological and gastrointestinal adverse reactions. Finally, alemtuzumab, bortezomib, belatacept and cellular therapies are in the search for the new treatments, whose premise is the induction of donor-specific nonresponse in the context of operational tolerance or mixed chimerism. The use of adapted and adequate immunosuppression has led to variable results and some are very encouraging; however, they must be validated with experimental studies

Type II Critical Collapse of a Self-Gravitating Nonlinear σ\sigma-Model

We report on the existence and phenomenology of type II critical collapse within the one-parameter family of SU(2) $\sigma$-models coupled to gravity. Numerical investigations in spherical symmetry show discretely self-similar (DSS) behavior at the threshold of black hole formation for values of the dimensionless coupling constant \ccbeta ranging from 0.2 to 100; at 0.18 we see small deviations from DSS. While the echoing period $\Delta$ of the critical solution rises sharply towards the lower limit of this range, the characteristic mass scaling has a critical exponent $\gamma$ which is almost independent of \ccbeta, asymptoting to $0.1185 \pm 0.0005$ at large \ccbeta. We also find critical scaling of the scalar curvature for near-critical initial data. Our numerical results are based on an outgoing-null-cone formulation of the Einstein-matter equations, specialized to spherical symmetry. Our numerically computed initial-data critical parameters $p^*$ show 2nd order convergence with the grid resolution, and after compensating for this variation in $p^*$, our individual evolutions are uniformly 2nd order convergent even very close to criticality.Comment: 23 pages, includes 10 postscript figure files, uses REVTeX, epsf, psfrag, and AMS math fonts (amstex + amssymb); to appear in PRD15. Summary of revisions from v2: fix wrong formula in figure 6 caption and y-axis label, also minor wording changes and update publication status of refs 5-