In the era of advanced electromagnetic and gravitational wave detectors, it has become increasingly important to effectively combine and study the impact of stellar evolution on binaries and dynamical systems of stars. Systematic studies dedicated to exploring uncertain parameters in stellar evolution are required to account for the recent observations of the stellar populations. We present a new approach to the commonly used Single-Star Evolution (SSE) fitting formulae, one that is more adaptable: Method of Interpolation for Single Star Evolution (METISSE). It makes use of interpolation between sets of pre-computed stellar tracks to approximate evolution parameters for a population of stars. We have used METISSE with detailed stellar tracks computed by the Modules for Experiments in Stellar Astrophysics (MESA), Bonn Evolutionary Code (BEC) and Cambridge STARS code. METISSE better reproduces stellar tracks computed using the STARS code compared to SSE, and is on average three times faster. Using stellar tracks computed with MESA and BEC, we apply METISSE to explore the differences in the remnant masses, the maximum radial expansion, and the main sequence lifetime of massive stars. We find that different physical ingredients used in the evolution of stars, such as the treatment of radiation dominated envelopes, can impact their evolutionary outcome. For stars in the mass range 9-100 M$_\odot$, the predictions of remnant masses can vary by up to 20 M$_\odot$, while the maximum radial expansion achieved by a star can differ by an order of magnitude between different stellar models.

Massive and very massive stars can play important roles in stellar populations by ejecting strong stellar winds and exploding in energetic phenomena. It is therefore of high importance that their behaviour is properly accounted for in synthetic model populations. Here we present nine grids of stellar evolutionary model sequences, together with finely resolved interpolated sequences and synthetic populations, of stars with 9-500 Msun and with metallicities ranging from Solar down to 1/250 Solar. The stellar models were computed with the 'Bonn' evolutionary code (covering core-hydrogen- and core-helium-burning phases, both complete). Post-processing for publication has been done using optimized methods, developed by our team, for massive and very massive stars. Interpolation and population synthesis were also performed on the models by our newly developed routine synStars. Eight of the grids represent slowly rotating massive stars with normal/classical evolution, while one grid represents fast rotating, chemically-homogeneously evolving models. Apart from the common stellar parameters such as mass, radius, surface temperature, luminosity and mass loss rate, we present stellar wind properties such as estimated wind velocity and kinetic energy of the wind. Additionally, we provide complete chemical yields of 34 isotopes, and estimates for the masses of the compact object remnants. The 'Bonn' Optimized Stellar Tracks (BoOST) project is published as simple tables - including stellar models, interpolated tracks and synthetic populations - thus ideal for further scientific applications. For example, star-formation studies could be done with BoOST to cover broad metallicity ranges, and so could be simulations of high-redshift galaxies. Additionally, gravitational-wave event rate predictions could be refined using BoOST by accounting for very massive stars at low-metallicity.

A population of binary black hole mergers has now been observed in gravitational waves by Advanced LIGO and Virgo. The masses of these black holes appear to show evidence for a pile-up between $30$--$45$ $M_\odot$ and a cut-off above $\sim 45$ $M_\odot$. One possible explanation for such a pile-up and subsequent cut-off are pulsational pair-instability supernovae (PPISNe) and pair-instability supernovae (PISNe) in massive stars. We investigate the plausibility of this explanation in the context of isolated massive binaries. We study a population of massive binaries using the rapid population synthesis software COMPAS, incorporating models for PPISNe and PISNe. Our models predict a maximum black hole mass of $40$ $M_\odot$. We expect $\sim 10$\% of all binary black hole mergers at redshift z = 0 will include at least one component that went through a PPISN (with mass $30$--$40$ $M_\odot$), constituting $\sim 20$--$50$\% of binary black hole mergers observed during the first two observing runs of Advanced LIGO and Virgo. Empirical models based on fitting the gravitational-wave mass measurements to a combination of a power law and a Gaussian find a fraction too large to be associated with PPISNe in our models. The rates of PPISNe and PISNe track the low metallicity star formation rate, increasing out to redshift $z = 2$. These predictions may be tested both with future gravitational-wave observations and with observations of superluminous supernovae.

Cosmic History has witnessed the lives and deaths of multiple generations of massive stars, all of them invigorating their host galaxies with ionizing photons, kinetic energy, fresh material and stellar-mass black holes. Ubiquitous engines as they are, Astrophysics needs a good understanding of their formation, evolution, properties and yields throughout the history of the Universe, and with decreasing metal content mimicking the environment at the earliest epochs. Ultimately, a physical model that could be extrapolated to zero metallicity would enable tackling long-standing questions such as "What did the First, very massive stars of the Universe look like?" or "What was their role in the re-ionization of the Universe?". Yet, most our knowledge of metal-poor massive stars is drawn from one single point in metallicity. Massive stars in the Small Magellanic Cloud (SMC, $\sim 1/5 Z_{\odot}$) currently serve as templates for low-metallicity objects in the early Universe, even though significant differences with respect to massive stars with poorer metal content have been reported. This White Paper summarizes the current knowledge on extremely (sub-SMC) metal poor massive stars, highlighting the most outstanding open questions and the need to supersede the SMC as standard. A new paradigm can be built from nearby extremely metal-poor galaxies that make a new metallicity ladder, but massive stars in these galaxies are out of reach to current observational facilities. Such task would require an L-size mission, consisting of a 10m-class space telescope operating in the optical and the ultraviolet ranges. Alternatively, we propose that ESA unites efforts with NASA to make the LUVOIR mission concept a reality, thus continuing the successful partnership that made Hubble Space Telescope one of the greatest observatories of all time.

We investigate the impact of uncertainty in the metallicity-specific star formation rate over cosmic time on predictions of the rates and masses of double compact object mergers observable through gravitational waves. We find that this uncertainty can change the predicted detectable merger rate by more than an order of magnitude, comparable to contributions from uncertain physical assumptions regarding binary evolution, such as mass transfer efficiency or supernova kicks. We statistically compare the results produced by the COMPAS population synthesis suite against a catalog of gravitational-wave detections from the first two Advanced LIGO and Virgo observing runs. We find that the rate and chirp mass of observed binary black hole mergers can be well matched under our default evolutionary model with a star formation metallicity spread of $0.39$ dex around a mean metallicity $\left<Z\right>$ that scales with redshift $z$ as $\left<Z\right>=0.035 \times 10^{-0.23 z}$, assuming a star formation rate of $0.01 \times (1+z)^{2.77} / (1+((1+z)/2.9)^{4.7}) \, \rm{M}_\odot$ Mpc$^{-3}$ yr$^{-1}$. Intriguingly, this default model predicts that 80\% of the approximately one binary black hole merger per day that will be detectable at design sensitivity will have formed through isolated binary evolution with only dynamically stable mass transfer, i.e., without experiencing a common-envelope event.

Context. Metal-poor massive stars are supposed to be progenitors of certain supernovae, gamma-ray bursts and compact object mergers, potentially contributing to the early epochs of the Universe with their strong ionizing radiation. However, they remain mainly theoretical as individual spectroscopic observations of such objects have rarely been carried out below the metallicity of the SMC. Aims. This work aims at exploring what our state-of-the-art theories of stellar evolution combined with those of stellar atmospheres predict about a certain type of metal-poor (0.02 Z$_{\odot}$) hot massive stars, the chemically homogeneously evolving ones, called TWUIN stars. Methods. Synthetic spectra corresponding to a broad range in masses (20-130 M$_{\odot}$) and covering several evolutionary phases from the zero-age main-sequence up to the core helium-burning stage were computed. Results. We find that TWUIN stars show almost no emission lines during most of their {core hydrogen-burning} lifetimes. Most metal lines are completely absent, including nitrogen. During their core helium-burning stage, lines switch to emission and even some metal lines (oxygen and carbon, but still almost no nitrogen) show up. Mass loss and clumping play a significant role in line-formation in later evolutionary phases, particularly during core helium-burning. Most of our spectra are classified as an early O type giant or supergiant, and we find Wolf-Rayet stars of type WO in the core helium-burning phase. Conclusions. An extremely hot, early O type star observed in a low-metallicity galaxy could be the outcome of chemically homogeneous evolution $-$ and therefore the progenitor of a long-duration gamma-ray burst or a type Ic supernova. TWUIN stars may play an important role in reionizing the Universe due to their being hot without showing prominent emission lines during the majority of their lifetimes.

Multiple stellar populations are observed in almost all globular-clusters, but the origin of this phenomenon is still debated. We investigate the role cool supergiants may have played. To do this, we combine two investigative methods: state-of-the-art massive stellar evolution and calculations of the hydrodynamic structure of the cluster-gas. This approach allows us to study how star-formation in young massive clusters depends on the energy- and mass-input of the first-generation of stars, while predicting the chemical composition of the second-generation. We find that the presence of massive (9-500 M$_{\odot}$) metal-poor supergiants in the young cluster leads to a star-formation episode within the first 4 Myr of the cluster's lifetime, that is, before the first core-collapse supernovae explode or the gas is expelled. The stellar winds accumulate in the cluster center, forming the second-generation there. Its composition is predicted to show variations in O & Na abundances, consistently with observations. The abundance of helium is, similarly to other scenarios involving massive stars, higher than what is referred from observations. Supposing dynamical removal of stars from the outskirts of the cluster, or applying a top-heavy initial-mass-function, we can predict a number ratio of the second-generation as high as 20-80%. The effect of metallicity is shown to be important, as the most luminous supergiants are only predicted at low-metallicity, thus limiting $-$ but not excluding $-$ the extent of a polluted second-generation at high-metallicity. These massive stars becoming black-holes suggests globular-clusters hosting gravitational-wave progenitors. Our scenario predicts a correlation between the mass of the cluster and the extent of the multiple population phenomenon.

Double neutron stars (DNSs) have been observed as Galactic radio pulsars, and the recent discovery of gravitational waves from the DNS merger GW170817 adds to the known DNS population. We perform rapid population synthesis of massive binary stars and discuss model predictions, including formation rates, mass distributions, and delay time distributions. We vary assumptions and parameters of physical processes such as mass transfer stability criteria, supernova kick distributions, remnant mass distributions and common-envelope energetics. We compute the likelihood of observing the orbital period-eccentricity distribution of the Galactic DNS population under each of our population synthesis models, allowing us to quantitatively compare the models. We find that mass transfer from a stripped post-helium-burning secondary (case BB) onto a neutron star is most likely dynamically stable. We also find that a natal kick distribution composed of both low (Maxwellian $\sigma=30\rm~km~s^{-1}$) and high ($\sigma=265\rm~km~s^{-1}$) components is preferred over a single high-kick component. We find that the observed DNS mass distribution can place strong constraints on model assumptions.

The recent discovery of the electromagnetic counterpart of the gravitational wave source GW170817, has demonstrated the huge informative power of multi-messenger observations. During the next decade the nascent field of multi-messenger astronomy will mature significantly. Around 2030, third generation gravitational wave detectors will be roughly ten times more sensitive than the current ones. At the same time, neutrino detectors currently upgrading to multi km^3 telescopes, will include a 10 km^3 facility in the Southern hemisphere that is expected to be operational around 2030. In this review, we describe the most promising high frequency gravitational wave and neutrino sources that will be detected in the next two decades. In this context, we show the important role of the Transient High Energy Sky and Early Universe Surveyor (THESEUS), a mission concept proposed to ESA by a large international collaboration in response to the call for the Cosmic Vision Programme M5 missions. THESEUS aims at providing a substantial advancement in early Universe science as well as playing a fundamental role in multi-messenger and time-domain astrophysics, operating in strong synergy with future gravitational wave and neutrino detectors as well as major ground- and space-based telescopes. This review is an extension of the THESEUS white paper (Amati et al. 2017), also in light of the discovery of GW170817/GRB170817A that was announced on October 16th, 2017.

Anomalous surface abundances are observed in a fraction of the low-mass stars of Galactic globular clusters, that may originate from hot-hydrogen-burning products ejected by a previous generation of massive stars. We present and investigate a scenario in which the second generation of polluted low-mass stars can form in shells around cool supergiant stars within a young globular cluster. Simulations of low-Z massive stars (M$_{\rm i}\sim$ 150$-$600 M$_{\odot}$) show that both core-hydrogen-burning cool supergiants and hot ionizing stellar sources are expected to be present simulaneously in young globular clusters. Under these conditions, photoionization-confined shells form around the supergiants. We find that the shell is gravitationally unstable on a timescale that is shorter than the lifetime of the supergiant, and the Bonnor-Ebert mass of the overdense regions is low enough to allow star formation. Since the low-mass stellar generation formed in this shell is made up of the material lost from the supergiant, its composition necessarily reflects the composition of the supergiant wind. We show that the wind contains hot-hydrogen-burning products, and that the shell-stars therefore have very similar abundance anomalies that are observed in the second generation stars of globular clusters. Considering the mass-budget required for the second generation star-formation, we offer two solutions. Either a top-heavy initial mass function is needed with an index of $-$1.71..$-$2.07. Alternatively, we suggest the shell-stars to have a truncated mass distribution, and solve the mass budget problem by justifiably accounting for only a fraction of the first generation. Even without forming a photoionizaton-confined shell, the cool supergiant stars predicted at low-Z could contribute to the pollution of the interstellar medium of the cluster from which the second generation was born.

This review describes the most common theories behind long-duration gamma-ray burst progenitors. I discuss two astrophysical scenarios: the collapsar and the magnetar models. According to their requirements, the progenitor should be an envelope-free massive star with a fast rotating, collapsing iron core. Such an object, called a TWUIN star, may be produced by chemically homogeneous evolution either from a massive single star or a massive binary system. Various outcomes of this evolutionary path (e.g. supernova explosions and gravitational wave production) are also mentioned, and directions for future research are suggested. In the era of multi-messenger astronomy, my hope is to present a timely overview on how stellar astrophysicists are searching for progenitor models of long-duration gamma-ray bursts, and what they have found so far.

The origin of gamma-ray bursts (GRBs) is still a fascinating field of research nowadays. While we have been collecting more and more observationally constrained properties of GRB-physics, new theoretical results on the progenitor evolution (be it stellar or compact object) have also emerged. I review some of the most promising progenitor theories for forming a short-duration GRB. A special emphasis is put on the hypothetical case of forming a short-duration GRB through the double black hole merger scenario -- in which case we may expect to observe a gravitational wave emission too. The chemically homogeneous channel for forming a black hole binary is discussed, and the stellar progenitors (so called TWUIN stars) are introduced. The birth place of these short-duration GRBs with a gravitational wave counterpart may be low-metallicity, starforming dwarf galaxies.

Recently it has been found that models of massive stars reach the Eddington limit in their interior, which leads to dilute extended envelopes. We perform a comparative study of the envelope properties of massive stars at different metallicities, with the aim to establish the impact of the stellar metallicity on the effect of envelope inflation. We analyse published grids of core-hydrogen burning massive star models computed with metallicities appropriate for massive stars in the Milky Way, the LMC and the SMC, the very metal poor dwarf galaxy I Zwicky 18, and for metal-free chemical composition. Stellar models of all the investigated metallicities reach and exceed the Eddington limit in their interior, aided by the opacity peaks of iron, helium and hydrogen, and consequently develop inflated envelopes. Envelope inflation leads to a redward bending of the zero-age main sequence and a broadening of the main sequence band in the upper part of the Hertzsprung-Russell diagram. We derive the limiting L/M-values as function of the stellar surface temperature above which inflation occurs, and find them to be larger for lower metallicity. While Galactic models show inflation above ~29 Msun, the corresponding mass limit for Population III stars is ~150 Msun. While the masses of the inflated envelopes are generally small, we find that they can reach 1-100 Msun in models with effective temperatures below ~8000 K, with higher masses reached by models of lower metallicity. Envelope inflation is expected to occur in sufficiently massive stars at all metallicities, and is expected to lead to rapidly growing pulsations, high macroturbulent velocities, and might well be related to the unexplained variability observed in Luminous Blue Variables like S Doradus and Eta Carina.

The Fermi collaboration identified a possible electromagnetic counterpart of the gravitational wave event of September 14, 2015. Our goal is to provide an unsupervised data analysis algorithm to identify similar events in Fermi's Gamma-ray Burst Monitor CTTE data stream. We are looking for signals that are typically weak. Therefore, they can only be found by a careful analysis of count rates of all detectors and energy channels simultaneously. Our Automatized Detector Weight Optimization (ADWO) method consists of a search for the signal, and a test of its significance. We developed ADWO, a virtual detector analysis tool for multi-channel multi-detector signals, and performed successful searches for short transients in the data-streams. We have identified GRB150522B, as well as possible electromagnetic candidates of the transients GW150914 and LVT151012. ADWO is an independently developed, unsupervised data analysis tool that only relies on the raw data of the Fermi satellite. It can therefore provide a strong, independent test to any electromagnetic signal accompanying future gravitational wave observations.

Massive rotating single stars with an initial metal composition appropriate for the dwarf galaxy I Zw 18 ([Fe/H]=$-$1.7) are modelled during hydrogen burning for initial masses of 9-300 M$_{\odot}$ and rotational velocities of 0-900 km s$^{-1}$. Internal mixing processes in these models were calibrated based on an observed sample of OB-type stars in the Magellanic Clouds. Even moderately fast rotators, which may be abundant at this metallicity, are found to undergo efficient mixing induced by rotation resulting in quasi chemically-homogeneous evolution. These homogeneously-evolving models reach effective temperatures of up to 90 kK during core hydrogen burning. This, together with their moderate mass-loss rates, make them Transparent Wind Ultraviolet INtense stars (TWUIN star), and their expected numbers might explain the observed HeII ionizing photon flux in I Zw 18 and other low-metallicity HeII galaxies. Our slowly rotating stars above $\sim$80 M$_{\odot}$ evolve into late B- to M-type supergiants during core hydrogen burning, with visual magnitudes up to 19$^{\mathrm{m}}$ at the distance of I Zw 18. Both types of stars, TWUIN stars and luminous late-type supergiants, are only predicted at low metallicity. Massive star evolution at low metallicity is shown to differ qualitatively from that in metal-rich environments. Our grid can be used to interpret observations of local star-forming dwarf galaxies and high-redshift galaxies, as well as the metal-poor components of our Milky Way and its globular clusters.

We present a method for determining the background of Fermi GBM GRBs using the satellite positional information and a physical model. Since the polynomial fitting method typically used for GRBs is generally only indicative of the background over relatively short timescales, this method is particularly useful in the cases of long GRBs or those which have Autonomous Repoint Request (ARR) and a background with much variability on short timescales. We give a Direction Dependent Background Fitting (DDBF) method for separating the motion effects from the real data and calculate the duration (T90 and T50, as well as confidence intervals) of the nine example bursts, from which two resulted an ARR. We also summarize the features of our method and compare it qualitatively with the official GBM Catalogue. Our background filtering method uses a model based on the physical information of the satellite position. Therefore, it has many advantages compared to previous methods. It can fit long background intervals, remove all the features caused by the rocking behaviour of the satellite, and search for long emissions or not-triggered events. Furthermore, many part of the fitting have now been automatised, and the method have been shown to work for both Sky Survey mode and ARR mode data. Future work will provide a burst catalogue with DDBF.

The Fermi Gamma-ray Burst Monitor (GBM) detects gamma-rays in the energy range 8 keV - 40 MeV. We developed a new background fitting process of these data, based on the motion of the satellite. Here we summarize this method, called Direction Dependent Background Fitting (DDBF), regarding the GBM triggered catalog. We also give some preliminary results and compare the duration parameters with the 2-years Fermi Burst Catalog.

The Swift satellite made a real break through with measuring simultaneously the gamma X-ray and optical data of GRBs, effectively. Although, the satellite measures the gamma, X-ray and optical properties almost in the same time a significant fractions of GRBs remain undetected in the optical domain. In a large number of cases only an upper bound is obtained. Survival analysis is a tool for studying samples where a part of the cases has only an upper (lower) limit. The obtained survival function may depend on some other variables. The Cox regression is a way to study these dependencies. We studied the dependence of the optical brightness (obtained by the UVOT) on the gamma and X-ray properties, measured by the BAT and XRT on board of the Swift satellite. We showed that the gamma peak flux has the greatest impact on the afterglow's optical brightness while the gamma photon index and the X-ray flux do not. This effect probably originates in the energetics of the jet launched from the central engine of the GRB which triggers the afterglow.