Line-driven stellar winds are ubiquitous among hot massive stars. In some cases they can become so strong, that the whole star is cloaked by an optically thick wind. The strong outflow gives rise to large emission lines, defining the class of so-called Wolf-Rayet (WR) stars. While being major players in the evolution of massive stars, the formation of heavy black holes, and the distribution of elements, the occurrence and nature of WR winds is still quite enigmatic.

A promising instrument towards a better theoretical understanding are stellar atmospheres allowing for a consistent inclusion of the hydrodynamics. By coupling stellar and wind parameters and the inclusion of a detailed non-LTE radiative transfer, they allow us to go beneath the observable layers and study the onset of WR-type winds. Establishing larger sets of models, we were able to make ground-breaking progress by identifying trends with mass and metallicity that deviate significantly from present empirical descriptions. Our modelling efforts reveal a complex picture for WR-type winds with strong, non-linear dependencies. Besides covering metallicity and mass, we further identify surface hydrogen as an important ingredient to retain WR-type mass loss at lower metallicity. Here, we present a summary of recent insights on the nature and onset of WR-type winds in massive stars including the consequences for stellar evolution, remaining open questions, and current efforts to overcome them.

Classical Wolf-Rayet (WR) stars mark an important stage in the late evolution of massive stars. As hydrogen-poor massive stars, these objects have lost their outer layers, while still losing further mass through strong winds indicated by their prominent emission line spectra. Wolf-Rayet stars have been detected in a variety of different galaxies. Their strong winds are a major ingredient of stellar evolution and population synthesis models. Yet, a coherent theoretical picture of their strong mass-loss is only starting to emerge. In particular, the occurrence of WR stars as a function of metallicity (Z) is still far from being understood.

To uncover the nature of the complex and dense winds of Wolf-Rayet stars, we employ a new generation of model atmospheres including a consistent solution of the wind hydrodynamics in an expanding non-LTE situation. With this technique, we can dissect the ingredients driving the wind and predict the resulting mass-loss for hydrogen-depleted massive stars. Our modelling efforts reveal a complex picture with strong, non-linear dependencies on the luminosity-to-mass ratio and Z with a steep, but not totally abrupt onset for WR-type winds in helium stars. With our findings, we provide a theoretical motivation for a population of helium stars at low Z, which cannot be detected via WR-type spectral features. Our study of massive He-star atmosphere models yields the very first mass-loss recipe derived from first principles in this regime. Implementing our first findings in stellar evolution models, we demonstrate how traditional approaches tend to overpredict WR-type mass loss in the young Universe.

To decode the information stored within a spectrum, detailed modelling of the physical state is required together with accurate radiative transfer solution schemes. In the analysis of stellar spectra, the numerical model often needs to account for high velocity outflows, multi-dimensional structures, and the effects of binary companions. Focusing now on binary systems, we present the BOSS-3D spectral synthesis code, which is capable of calculating synthetic line profiles for a variety of binary systems. Assuming the state of the circumstellar material to be known, the standard pz-geometry is extended by defining individual coordinate systems for each object. By embedding these coordinate systems within the observer’s frame, BOSS-3D automatically accounts for outflows or discs within both involved systems, and includes all Doppler shifts. Moreover, the code accounts for different length-scales, and thus could also be used to analyse transit-spectra of planetary atmospheres. As a first application of BOSS-3D, we model the phase-dependent line profiles for the enigmatic binary (or multiple) system LB-1.

Over the last decade, our abilities to observe the internal kinematics of star clusters have drastically increased. Where a few years back only small numbers of bright stars with radial velocity measurements were available, we can nowadays study the three-dimensional motions of large stellar samples, thanks to the combined datasets gathered by state-of-art spectrographs and astrometric satellites. In this work, I summarise the contribution of integral-field spectrographs, in particular MUSE, to this paradigm change. Using dedicated software tools, we were able to overcome a fundamental limitation of spectroscopy and advance to the crowded cluster centres. This allowed us to study the central kinematics in unprecedented detail and to start uncovering the populations of black holes that reside in massive star clusters.

We investigate the dissolution process of star clusters embedded in an external tidal field and harboring a subsystem of stellar-mass black hole. For this purpose we analyzed the MOCCA models of real star clusters contained in the Mocca Survey Database I. We showed that the presence of a stellar-mass black hole subsystem in tidally filling star cluster can lead to abrupt cluster dissolution connected with the loss of cluster dynamical equilibrium. Such cluster dissolution can be regarded as a third type of cluster dissolution mechanism. We additionally argue that such a mechanism should also work for tidally under-filling clusters with a top-heavy initial mass function.

Sizeable number of stellar-mass black holes (BHs) in globular clusters (GCs) can strongly influence the dynamical evolution and observational properties of their host cluster. Using results from a large set of numerical simulations, we identify the key ingredients needed to sustain a sizeable population of BHs in GCs up to a Hubble time. We find that while BH natal kick prescriptions are essential in determining the initial retention fraction of BHs in GCs, the long-term survival of BHs is determined by the size, initial central density and half-mass relaxation time of the GC. Simulated GC models that contain many BHs are characterized by relatively low central surface brightness, large half-light and core radii values. We also discuss novel ways to compare simulated results with available observational data to identify GCs that are most likely to contain many BHs.

With the hundreds of merging binary black holes (BHs) expected to be detected by LIGO, LISA, and other upcoming instruments, the modelling of astrophysical channels that lead to the formation of compact BH binaries has become of crucial importance. BHs of any size can form bound systems in every astrophysical environment, from the field to galactic nuclei. If a binary is too wide, it needs a catalysis process to harden and merge, as in the case a third objects orbiting the BH binary on a distant orbit. In this case, Kozai-Lidov cycles can pump up the binary eccentricity, thus driving it to a merger thanks to efficient energy dissipation at the pericenter. Some remarkable scenarios where the Kozai-Lidov mechanism operates are in triple and quadruple systems of stellar BHs, and in intermediate-mass BH-stellar BH binaries in orbit around a central supermassive BH in galactic nuclei.

We present the study of horizontal branch morphology of the cluster NGC 6656. A blueward shift in temperature of about ~5000 K (nM-jump) in the color-color plot is detected.To explain this feature, we study the presence of stellar-mass black hole by plotting Projected density profile (PDP) in the central HST region. The PDP in the inner region (r < 10″) can be nicely reproduced by the king+BH model. The blue ward shift in temperature can be due the presence of stellar mass black holes in the centre.

This is an introduction to Bayesian inference with a focus on hierarchical models and hyper-parameters. We write primarily for an audience of Bayesian novices, but we hope to provide useful insights for seasoned veterans as well. Examples are drawn from gravitational-wave astronomy, though we endeavour for the presentation to be understandable to a broader audience. We begin with a review of the fundamentals: likelihoods, priors, and posteriors. Next, we discuss Bayesian evidence, Bayes factors, odds ratios, and model selection. From there, we describe how posteriors are estimated using samplers such as Markov Chain Monte Carlo algorithms and nested sampling. Finally, we generalise the formalism to discuss hyper-parameters and hierarchical models. We include extensive appendices discussing the creation of credible intervals, Gaussian noise, explicit marginalisation, posterior predictive distributions, and selection effects.

We report our investigations on the host galaxy properties of stellar binary black holes (SBBHs) by implementing simple recipes for SBBH formation and merger into cosmological galaxy formation model. If the time delay between SBBH formation and merger ranges from Gyr to the Hubble time, SBBH mergers at redshift z < 0.3 occur preferentially in big galaxies with stellar mass M* > 2 × 1010M⊙ and metallicities Z peaking around ~0.6Z⊙. However, the host galaxy stellar mass distribution of heavy SBBH mergers (with total black hole mass >50M⊙) is bimodal with one peak at ~109M⊙ and the other peak at ~2 × 1010M⊙. The contribution fraction from metal-poor host galaxies (Z < 0.2Z⊙) to heavy mergers is much larger than that to less heavy mergers. If SBBHs were formed in the early universe, their mergers detected at z < 0.3 occur preferentially in even more massive galaxies with M* > 3 × 1010M⊙ and in galaxies with metallicities mostly >0.2Z⊙ and peaking at Z ~ 0.6Z⊙.

We perform hydrodynamic supernova (SN) simulations in spherical symmetry for progenitor models with solar metallicity across the stellar mass range from 9.0 to 120 M⊙ to explore the progenitor-explosion and progenitor-remnant connections based on the neutrino-driven mechanism. We use an approximative treatment of neutrino transport and replace the high-density interior of the neutron star (NS) by an inner boundary condition based on an analytic proto-NS core-cooling model, whose free parameters are chosen to reproduce the observables of SN 1987A and the Crab SN for theoretical models of their progenitor stars.

Judging the fate of a massive star, either a neutron star (NS) or a black hole (BH), solely by its structure prior to collapse has been ambiguous. Our work and previous attempts find a non-monotonic variation of successful and failed supernovae with zero-age main-sequence mass. We identify two parameters based on the “critical luminosity” concept for neutrino-driven explosions, which in combination allows for a clear separation of exploding and non-exploding cases.

Continuing our simulations beyond shock break-out, we are able to determine nucleosynthesis, light curves, explosion energies, and remnant masses. The resulting NS initial mass function has a mean gravitational mass near 1.4 M⊙. The average BH mass is about 9 M⊙ if only the helium core implodes, and 14 M⊙ if the entire pre-SN star collapses. Only ~10% of SNe come from stars over 20 M⊙, and some of these are Type Ib or Ic.

Recently, the Laser Interferometer Gravitational-Wave Observatory (LIGO) has detected the gravitational wave (GW) event, GW150914, as a result of the merger of a ~30 M⊙ black hole (BH) binary. So far, the merger of stellar-mass BHs has been thought to result from the evolution of binary stars. Here, we propose a novel path of the merger stemming from non-binary isolated stars. In our previous studies, we have found that multiple non-binary stellar-mass BHs whose separations are larger than 1000 AU can merge with each other under the gas-rich environments through the gas dynamical friction and three-body interaction. In this case, a considerable amount of gas can accrete onto BHs before the merger, that is, the initial mass of BHs can be lower than 30 M⊙. To explore this possibility, we perform post-Newtonian N-body simulations on mergers of accreting stellar-mass BHs. Based on our simulations, we find that the BH merger in GW150914 from smaller seed BHs is likely to occur in galactic nuclear regions or dense interstellar cloud cores. Furthermore, we roughly estimate event rates to be ~0.4 yr−1 in galactic nuclear regions and ~8 yr−1 in dense interstellar cloud cores.

We numerically investigate the thermally unstable accretion discs around spinning black holes with different spins. We adopted an additional evolutionary viscosity equation to replace the standard alpha-prescription based on the results of two MHD simulations. We find an interesting oscillation when accretion switches to slim disc mode. The oscillation arises from the sonic point of accretion flow and propagates outwards. We mimic the bolometric light-curve and find a series of harmonics on its power spectrum. The frequency ratio of those harmonics is a regular integer series. The lowest frequency of the harmonics is identical to the prediction of trapped p-mode in QPO theory.