Common-envelope (CE) evolution in massive binary systems is thought to be one of the most promising channels for the formation of compact binary mergers. In the case of merging binary black holes (BBHs), the essential CE phase takes place at a stage when the first BH is already formed and the companion star expands as a supergiant. We study which BH binaries with supergiant companions will evolve through and potentially survive a CE phase. To this end, we compute envelope binding energies from detailed massive stellar models at different evolutionary stages and metallicities. We make multiple physically extreme choices of assumptions that favor easier CE ejection as well as account for recent advancements in mass transfer stability criteria. We find that even with the most optimistic assumptions, a successful CE ejection in BH (and also NS) binaries is only possible if the donor is a massive convective-envelope giant, a red supergiant (RSG). In other words, pre-CE progenitors of BBH mergers are BH binaries with RSG companions. We find that due to its influence on the radial expansion of massive giants, metallicity has an indirect but a very strong effect on the envelope structure and binding energies of RSGs. Our results suggest that merger rates from population synthesis models could be severely overestimated, especially at low metallicity. Additionally, the lack of observed RSGs with luminosities above log($L/L_{\odot}$) = 5.6-5.8, corresponding to stars with $M > 40 M_{\odot}$, puts into question the viability of the CE channel for the formation of the most massive BBH mergers. Either such RSGs elude detection due to very short lifetimes, or they do not exist and the CE channel can only produce BBH systems with masses $< 50 M_{\odot}$. Finally, we discuss an alternative CE scenario, in which a partial envelope ejection is followed by a phase of possibly long and stable mass transfer.

Metallicity is known to significantly affect the radial expansion of a massive star: the lower the metallicity, the more compact the star, especially during its post-MS evolution. We study this effect in the context of binary evolution. Using the stellar-evolution code MESA, we computed evolutionary tracks of stars at different metallicities, exploring variations of factors known to affect the radial expansion (e.g. semiconvection, overshooting, rotation). We find observational support for an evolution in which already at metallicity $0.2Z_{\odot}$ massive stars remain relatively compact during the Hertzprung-Gap (HG) phase and most of their expansion occurs during core-helium burning (CHeB). Consequently, we show that metallicity has a strong influence on the type of mass transfer evolution in binary systems. At solar metallicity, a case-B mass transfer is initiated shortly after the end of MS, and a giant donor is almost always a rapidly expanding HG star. At lower metallicity, the parameter space for mass transfer from a more evolved CHeB star increases dramatically. This means that envelope stripping and formation of helium stars in low-metallicity environments occurs later in the evolution of the donor, implying a much shorter duration of the Wolf-Rayet phase (even by an order of magnitude) and higher final core masses. This metallicity effect is independent of the impact of metallicity-dependent stellar winds. At very low metallicities, a significant fraction of massive stars in binaries engages in the first episode of mass transfer very late into their evolution, when they already have a well-developed CO core. The remaining lifetime ($< 10^4$ yr) is unlikely to be enough to strip the entire H-rich envelope. We also briefly discuss the extremely small parameter space for mass transfer from massive convective-envelope donors in the context of binary black hole merger formation.

All ten LIGO/Virgo binary black hole (BH-BH) coalescences reported from the O1/O2 runs have near zero effective spins. There are only three potential explanations of this fact. If the BH spin magnitudes are large then (i) either both BH spin vectors must be nearly in the orbital plane or (ii) the spin angular momenta of the BHs must be oppositely directed and similar in magnitude. Or, (iii) the BH spin magnitudes are small. We test the third hypothesis within the framework of the classical isolated binary evolution scenario of the BH-BH merger formation. We test three models of angular momentum transport in massive stars: a mildly efficient transport by meridional currents (as employed in the Geneva code), an efficient transport by the Tayler-Spruit magnetic dynamo (as implemented in the MESA code), and a very-efficient transport (as proposed by Fuller et al.) to calculate natal BH spins. We allow for binary evolution to increase the BH spins through accretion and account for the potential spin-up of stars through tidal interactions. Additionally, we update the calculations of the stellar-origin BH masses, include revisions to the history of star formation and to the chemical evolution across cosmic time. We find that we can match simultaneously the observed BH-BH merger rate density, BH masses, and effective spins. Models with efficient angular momentum transport are favored. The updated stellar-mass weighted gas-phase metallicity evolution now used in our models appears to be a key in better reproducing the LIGO/Virgo merger rate estimate. Mass losses during the pair-instability pulsation supernova phase are likely overestimated if the merger GW170729 hosts a BH more massive than 50 Msun. We also estimate rate of BH-NS mergers from recent LIGO/Virgo observations. Our updated models of BH-BH, BH-NS and NS-NS mergers are now publicly available at www.syntheticuniverse.org.

Gaia16aye was a binary microlensing event discovered in the direction towards the northern Galactic disc and was one of the first microlensing events detected and alerted to by the Gaia space mission. Its light curve exhibited five distinct brightening episodes, reaching up to I=12 mag, and it was covered in great detail with almost 25,000 data points gathered by a network of telescopes. We present the photometric and spectroscopic follow-up covering 500 days of the event evolution. We employed a full Keplerian binary orbit microlensing model combined with the motion of Earth and Gaia around the Sun to reproduce the complex light curve. The photometric data allowed us to solve the microlensing event entirely and to derive the complete and unique set of orbital parameters of the binary lensing system. We also report on the detection of the first-ever microlensing space-parallax between the Earth and Gaia located at L2. The properties of the binary system were derived from microlensing parameters, and we found that the system is composed of two main-sequence stars with masses 0.57$\pm$0.05 $M_\odot$ and 0.36$\pm$0.03 $M_\odot$ at 780 pc, with an orbital period of 2.88 years and an eccentricity of 0.30. We also predict the astrometric microlensing signal for this binary lens as it will be seen by Gaia as well as the radial velocity curve for the binary system. Events such as Gaia16aye indicate the potential for the microlensing method of probing the mass function of dark objects, including black holes, in directions other than that of the Galactic bulge. This case also emphasises the importance of long-term time-domain coordinated observations that can be made with a network of heterogeneous telescopes.

In large and complicated stellar systems like galaxies it is difficult to predict the number and characteristics of a black hole population. Such populations may be modelled as an aggregation of homogeneous (i.e. having uniform star formation history and the same initial chemical composition) stellar populations. Using realistic evolutionary models we predict the abundances and properties of black holes formed from binaries in these environments. We show that the black hole population will be dominated by single black holes originating from binary disruptions and stellar mergers. Furthermore, we discuss how black hole populations are influenced by such factors as initial parameters, metallicity, initial mass function, and natal kick models. As an example application of our results, we estimate that about 26 microlensing events to happen every year in the direction of the Galactic Bulge due to black holes in a survey like OGLE-IV. Our results may be used to perform in-depth studies related to realistic black hole populations, e.g. observational predictions for space survey missions like Gaia, or Einstein Probe. We prepared a publicly available database with the raw data from our simulations to be used for more in-depth studies.

We briefly describe and discuss the set-up of the project MOCCA Survey Database I. The database contains more than 2000 Monte Carlo models of evolution of real star cluster performed with the MOCCA code. Then, we very briefly discuss results of analysis of the database regarding the following projects: formation of intermediate mass black holes, abrupt cluster dissolution harboring black hole subsystems, retention fraction of black hole - black hole mergers, and tidal disruption events with intermediate mass black holes.

X-ray binaries with black hole (BH) accretors and massive star donors at short orbital periods of a few days can evolve into close binary BH systems (BBH) that merge within the Hubble time through stable mass transfer evolution. From observational point of view, upon the Roche-lobe overflow, such systems will most likely appear as ultra-luminous X-ray sources (ULXs). To study this connection, we compute the mass transfer phase in systems with BH accretors and massive star donors ($M > 15 \,M_\odot$) at various orbital separations and metallicities using the MESA stellar evolution code. In the case of core-hydrogen and core-helium burning donors (cases A and C of mass transfer) we find the typical duration of super-Eddington mass transfer of up to $10^6$ and $10^5 \, \rm yr$ , with rates of $10^{-6}$ and $10^{-5} \,M_\odot \, \rm yr^{-1}$ , respectively. Given that roughly $0.5$ ULXs are found per unit of star formation rate ($\,M_\odot \, \rm yr^{-1}$), and assuming that $10\%$ of all the observed ULXs form merging BBH, we estimate the rate of BBH mergers from stable mass transfer evolution to be at most $10 {\rm ~Gpc}^{-3} {\rm ~yr}^{-1}$.

The distributions of the initial main-sequence binary parameters are one of the key ingredients in obtaining evolutionary predictions for compact binary (BH-BH / BH-NS / NS-NS) merger rates. Until now, such calculations were done under the assumption that initial binary parameter distributions were independent. Here, we implement empirically derived inter-correlated distributions of initial binary parameters primary mass (M1), mass ratio (q), orbital period (P), and eccentricity (e). Unexpectedly, the introduction of inter-correlated initial binary parameters leads to only a small decrease in the predicted merger rates by a factor of 2 $-$ 3 relative to the previously used non-correlated initial distributions. The formation of compact object mergers in the isolated classical binary evolution favors initial binaries with stars of comparable masses (q = 0.5 $-$ 1) at intermediate orbital periods (log P (days) = 2 $-$ 4). New distributions slightly shift the mass ratios towards smaller values with respect to the previously used flat q distribution, which is the dominant effect decreasing the rates. New orbital periods only negligibly increase the number of progenitors. Additionally, we discuss the uncertainty of merger rate predictions associated with possible variations of the massive-star initial mass function (IMF). We argue that evolutionary calculations should be normalized to a star formation rate (SFR) that is obtained from the observed amount of UV light at wavelength 1500{\AA} (SFR indicator). In this case, contrary to recent reports, the uncertainty of the IMF does not affect the rates by more than a factor of 2. Any change to the IMF slope for massive stars requires a change of SFR in a way that counteracts the impact of IMF variations on the merger rates. In contrast, we suggest that the uncertainty in cosmic SFR at low metallicity can be a significant factor at play.

We investigate the impact of combining Gaia astrometry from space with precise, high cadence OGLE photometry from the ground. For the archival event OGLE3-ULENS-PAR-02, which is likely a black hole, we simulate a realistic astrometric time-series of Gaia measurements and combine it with the real photometric data collected by the OGLE project. We predict that at the end of the nominal 5 years of the Gaia mission, for the events brighter than $G\approx15.5$ mag at the baseline, caused by objects heavier than 10 $M_{\odot}$, it will be possible to unambiguously derive masses of the lenses, with accuracy between a few to 15 per cent. We find that fainter events ($G<17.5$) can still have their lens masses determined, provided that they are heavier than 30 $M_{\odot}$. We estimate that the rate of astrometric microlensing events caused by the stellar-origin black holes is $\approx 4 \times 10^{-7} \, \rm yr^{-1}$, which implies, that after 5 years of Gaia operation and $\approx 5 \times 10^6$ bright sources in Gaia, it will be possible to identify few such events in the Gaia final catalogues.

We revisit double neutron star (DNS) formation in the classical binary evolution scenario in light of the recent LIGO/Virgo DNS detection (GW170817). The observationally estimated Galactic DNS merger rate of $R_{\rm MW}=21^{+28}_{-14}$ Myr$^{-1}$, based on 3 Galactic DNS systems, fully supports our standard input physics model with $R_{\rm MW} =24$ Myr$^{-1}$. This estimate for the Galaxy translates in a non-trivial way (due to cosmological evolution of progenitor stars in chemically evolving Universe) into a local ($z\approx0$) DNS merger rate density of $R_{\rm local}=48$ Gpc$^{-3}$yr$^{-1}$, which {\em is not} consistent with the current LIGO/Virgo DNS merger rate estimate ($1540^{+3200}_{-1220}$ Gpc$^{-3}$yr$^{-1}$). Within our study of the parameter space we find solutions that allow for DNS merger rates as high as $R_{\rm local} \approx 600^{+600}_{-300}$ Gpc$^{-3}$yr$^{-1}$ which are thus consistent with the LIGO/Virgo estimate. However, our corresponding BH-BH merger rates for the models with high DNS merger rates exceed the current LIGO/Virgo estimate of local BH-BH merger rate ($12$-$213$ Gpc$^{-3}$yr$^{-1}$). Apart from being particularly sensitive to the common envelope treatment, DNS merger rates are rather robust against variations of several of the key factors probed in our study (e.g. mass transfer, angular momentum loss, natal kicks). This might suggest that either common envelope development/survival works differently for DNS ($\sim$ 10-20 Msun stars) than for BH-BH ($\sim$ 40-100 Msun stars) progenitors, or high BH natal kicks are needed to meet observational constraints for both types of binaries. Note that our conclusion is based on a limited number of (21) evolutionary models and is valid only within this particular DNS and BH-BH isolated binary formation scenario.

Very wide binaries (> 500 AU) are subject to numerous encounters with flying-by stars in the Galactic field and can be perturbated into highly eccentric orbits (e ~ 0.99). For such systems tidal interactions at close pericenter passages can lead to orbit circularization and possibly mass transfer, consequently producing X-Ray binaries without the need for common envelope. We test this scenario for the case of Black Hole Low-Mass X-Ray Binaries (BH LMXBs) by performing a population synthesis from primordial binaries with numerical treatment of random stellar encounters. We test various models for the threshold pericenter distance under which tidal forces cause circularization. We estimate that fly-by interactions can produce a current population of ~ 60$-$220 BH LMXBs in the Galactic field. The results are sensitive to the assumption on tidal circularization efficiency and zero to very small BH natal kicks of a few km/s are required. We show that the most likely donors are low-mass stars (< 1 Msun; at the onset of mass transfer) as observed in the population of known sources (~ 20). However, the low number of systems formed along this route is in tension with most recent observational estimate of the number of dormant BH LMXBs in the Galaxy 10$^4-$10$^8$ (Tetarenko et al. 2016). If indeed the numbers are so high, alternative formation channels of BHs with low-mass donors need to be identified.

Real-time analysis and classification of observational data collected within synoptic sky surveys is a huge challenge due to constant growth of data volumes. Machine learning techniques are often applied in order to perform this task automatically. The current bottleneck of transients detection in most surveys is the process of filtering numerous artifacts from candidate detection. We present a new method for automated artifact filtering based on hierarchical unsupervised classifier employing Self-Organizing Maps (SOMs). The system accepts 97 % of real transients and removes 97.5 % of artifacts when tested on the OGLE-IV Transient Detection System. The improvement of the artifacts filtering allowed for single-frame based rapid detections of transients within OGLE-IV, which now alerts on transient discoveries in less than 15 minutes from the image acquisition.

The current bottleneck of transient detection in most surveys is the problem of rejecting numerous artifacts from detected candidates. We present a triple-stage hierarchical machine learning system for automated artifact filtering in difference imaging, based on self-organizing maps. The classifier, when tested on the OGLE-IV Transient Detection System, accepts ~ 97 % of real transients while removing up to ~ 97.5 % of artifacts.