ESA and NASA are moving forward with plans to launch LISA around 2034. With data from the Illustris cosmological simulation, we provide analysis of LISA detection rates accompanied by characterization of the merging massive black hole population. Massive black holes of total mass $\sim10^5-10^{10} M_\odot$ are the focus of this study. We evolve Illustris massive black hole mergers, which form at separations on the order of the simulation resolution ($\sim$kpc scales), through coalescence with two different treatments for the binary massive black hole evolutionary process. The coalescence times of the population, as well as physical properties of the black holes, form a statistical basis for each evolutionary treatment. From these bases, we Monte Carlo synthesize many realizations of the merging massive black hole population to build mock LISA detection catalogs. We analyze how our massive black hole binary evolutionary models affect detection rates and the associated parameter distributions measured by LISA. With our models, we find massive black hole binary detection rates with LISA of $\sim0.5-1$ yr$^{-1}$ for massive black holes with masses greater than $10^5M_\odot$. This should be treated as a lower limit primarily because our massive black hole sample does not include masses below $10^5M_\odot$, which may significantly add to the observed rate. We suggest reasons why we predict lower detection rates compared to much of the literature.

Most studies of mass transfer in binary systems assume circular orbits at the onset of Roche lobe overflow. However, there are theoretical and observational indications that mass transfer could occur in eccentric orbits. In particular, eccentricity could be produced via sudden mass loss and velocity kicks during supernova explosions, or Lidov-Kozai (LK) oscillations in hierarchical triple systems, or, more generally, secular evolution in multiple-star systems. However, current analytic models of eccentric mass transfer are faced with the problem that they are only well defined in the limit of very high eccentricities, and break down for less eccentric and circular orbits. This provides a major obstacle to implementing such models in binary and higher-order population synthesis codes, which are useful tools for studying the long-term evolution of a large number of systems. Here, we present a new analytic model to describe the secular orbital evolution of binaries undergoing conservative mass transfer. The main improvement of our model is that the mass transfer rate is a smoothly varying function of orbital phase, rather than a delta function centered at periapsis. Consequently, our model is in principle valid for any eccentricity, thereby overcoming the main limitation of previous works. We implement our model in an easy-to-use and publicly available code that can be used as a basis for implementations of our model into population synthesis codes. We investigate the implications of our model in a number of applications with circular and eccentric binaries, and triples undergoing LK oscillations.

Nuclear star clusters around massive black holes are expected to be abundant in stellar mass black holes and black hole binaries. These binaries form a hierarchical triple system with the massive black hole at the center. Gravitational perturbations from the massive black hole can cause high eccentricity excitation. During this process, the eccentricity may approach unity, and the pericenter distance may become sufficiently small that gravitational wave emission drives the binary to merge. In this paper, we consider a simple proof of concept and explore the effect of the eccentric Kozai-Lidov mechanism for unequal mass binaries. We perform a set of Monte Carlo simulations on BH-BH binaries in galactic nuclei with quadrupole and octupole-level secular perturbations, general relativistic precession, and gravitational wave emission. For a nominal number of steady-state BH-BH binaries, our model gives a total merger rate $\sim 1 - 3$$Gpc^{-3} yr^{-1}$, depending on the assumed density profile. Thus, our model potentially competes with other dynamical mechanisms, such as the dynamical formations and mergers of BH binaries in globular clusters or dense nuclear clusters without a massive black hole. We provide predictions for the distributions of these LIGO sources in galactic nuclei.

Many giant exoplanets are found near their Roche limit and in mildly eccentric orbits. In this study we examine the fate of such planets through Roche-lobe overflow as a function of the physical properties of the binary components, including the eccentricity and the asynchronicity of the rotating planet. We use a direct three-body integrator to compute the trajectories of the lost mass in the ballistic limit and investigate the possible outcomes. We find three different outcomes for the mass transferred through the Lagrangian point $L_{1}$: (i) self-accretion by the planet, (ii) direct impact on the stellar surface, (iii) disk formation around the star. We explore the parameter space of the three different regimes and find that at low eccentricities, $e\lesssim 0.2$, mass overflow leads to disk formation for most systems, while for higher eccentricities or retrograde orbits self-accretion is the only possible outcome. We conclude that the assumption often made in previous work that when a planet overflows its Roche lobe it is quickly disrupted and accreted by the star is not always valid.

The supermassive black holes originally in the nuclei of two merging galaxies will form a binary in the remnant core. The early evolution of the massive binary is driven by dynamical friction before the binary becomes "hard" and eventually reaches coalescence through gravitational wave emission. { We consider the dynamical friction evolution of massive binaries consisting of a secondary hole orbiting inside a stellar cusp dominated by a more massive central black hole.} In our treatment we include the frictional force from stars moving faster than the inspiralling object which is neglected in the standard Chandrasekhar's treatment. We show that the binary eccentricity increases if the stellar cusp density profile rises less steeply than $\rho\propto r^{-2}$. In cusps shallower than $\rho\propto r^{-1}$ the frictional timescale can become very long due to the deficit of stars moving slower than the massive body. Although including the fast stars increases the decay rate, low mass-ratio binaries ($q\lesssim 10^{-3}$) in sufficiently massive galaxies have decay timescales longer than one Hubble time. During such minor mergers the secondary hole stalls on an eccentric orbit at a distance of order one tenth the influence radius of the primary hole (i.e., $\approx 10-100\rm pc$ for massive ellipticals). We calculate the expected number of stalled satellites as a function of the host galaxy mass, and show that the brightest cluster galaxies should have $\gtrsim 1$ of such satellites orbiting within their cores. Our results could provide an explanation to a number of observations, which include multiple nuclei in core ellipticals, off-center AGNs and eccentric nuclear disks.

We have endeavoured to understand the formation and evolution of the black hole (BH) X-ray binary LMC X-3. We estimate the properties of the system at 4 evolutionary stages: 1) at the Zero Age Main Sequence (ZAMS), 2) just prior to the supernova (SN) explosion of the primary, 3) just after the SN, and 4) at the moment of RLO onset.We use a hybrid approach, combining detailed stellar structure and binary evolution calculations with approximate population synthesis models. This allows us to estimate potential natal kicks and the evolution of the BH spin. In the whole analysis we incorporate as model constraints the most up-to-date observational information, encompassing the binary's orbital properties, the companion star mass, effective temperature, surface gravity and radius, as well as the black hole's mass and spin. We find that LMC X-3 began as a ZAMS system with the mass of the primary star in the range $M_{\rm{1,ZAMS}}$ = 22-31 $\rm{M_{\odot}}$ and a secondary star of $M_{\rm{2,ZAMS}} = 5.0-8.3M_{\odot}$, in a wide ($P_{ZAMS} \gtrsim 2.000\, \rm days$) and eccentric ($e_{\rm{ZAMS}} \gtrsim 0.23$) orbit. Just prior to the SN, the primary has a mass of $M_{\mathrm{1,preSN}} = 11.1-18.0\,\rm M_{\odot}$, with the secondary star largely unaffected. The orbital period decreases to $0.6-1.7\, \rm days$, and is still eccentric $0 \leq e_{\rm{preSN}} \leq 0.44$. We find that a symmetric SN explosion with no or small natal kicks (a few tens of $\rm km\, s^{-1}$) imparted on the BH cannot be formally excluded, however, large natal kicks in excess of $\gtrsim 120 \,\rm km\, s^{-1}$ increase the estimated formation rate by an order of magnitude. Following the SN, the system has a BH $M_{\mathrm{BH,postSN}} = 6.4-8.2\,\rm M_{\odot}$ and is put into an eccentric orbit. At the RLO onset the orbit is circularised and it has an orbital period of $P_{\rm{RLO}} = 0.8-1.4\,\rm days$.

Observations reveal that mass-transferring binary systems may have non-zero orbital eccentricities. The time-evolution of the orbital semi-major axis and eccentricity of mass-transferring eccentric binary systems is an important part of binary evolution theory and has been widely studied. However, various different approaches and assumptions on the subject have made the literature difficult to comprehend and comparisons between different orbital element time-evolution equations not easy to make. Consequently, no self-consistent treatment of this phase has been ever included in binary population synthesis codes. In this paper, we present a general formalism to derive the time-evolution equations of the binary orbital elements, treating mass-loss and mass-transfer as perturbations to the general two-body problem. We present the self-consistent form of the perturbing acceleration and the phase-dependent time-evolution equations for the orbital elements under different mass-loss/transfer processes. First, we study the cases of isotropic and anisotropic wind mass-loss. Then, we proceed with the non-isotropic ejection and accretion in conservative as well as non-conservative manner for both point masses and extended bodies. Comparison of the derived equations with similar work in the literature is made and explanation of the existing discrepancies is provided.

Finite eccentricities in mass-transferring eccentric binary systems can be explained by taking into account mass-loss and mass-transfer processes that often occur in these systems. These processes can be treated as perturbations to the general two-body problem. The time-evolution equations for the semi-major axis and the eccentricity derived from perturbative methods are in general phase-dependent. The osculating semi-major axis and eccentricity change over the orbital timescale and they are not easy to implement in binary evolution codes like MESA. However, the secular orbital element evolution equations can be simplified averaging over the rapidly varying true anomalies. In this paper, we derive the secular time-evolution equations for the semi-major axis and the eccentricity for various mass-loss/transfer processes using either the adiabatic approximation or the assumption of delta-function mass-loss/transfer at periastron. We begin with the cases of isotropic and anisotropic wind mass-loss. We continue with conservative and non-conservative non-isotropic mass ejection/accretion (including RLOF) for both point-masses and extended bodies. We conclude with the case of phase-dependent mass accretion. Comparison of the derived equations with similar work in the literature is included and explanation of the existing discrepancies is provided.

We use a general relativistic approach to investigate the effects of weak cosmological magnetic fields on linear rotational perturbations during the radiation and dust epochs of the universe. This includes ordinary kinematic vorticity, as well as vortex-like inhomogeneities in the density distribution of the cosmic medium. Our study confirms that magnetic fields can source both types of perturbations and shows that their presence also helps cosmic rotation to survive longer. In agreement with previous Newtonian studies, we find that during the dust era vorticity decays slower than in non-magnetised cosmologies. Prior to equipartition the magnetic effect is more pronounced, since it helps the aforementioned rotational distortions to maintain constant magnitude. Overall, magnetised universes seem to provide a much better environment for the survival of cosmic vorticity than their magnetic-free counterparts.

We study the role of viscosity and the effects of a magnetic field on a rotating, self-gravitating fluid, using Newtonian theory and adopting the ideal magnetohydrodynamic approximation. Our results confirm that viscosity can generate vorticity in inhomogeneous environments, while the magnetic tension can produce vorticity even in the absence of fluid pressure and density gradients. Linearizing our equations around an Einstein-de Sitter cosmology, we find that viscosity adds to the diluting effect of the universal expansion. Typically, however, the dissipative viscous effects are confined to relatively small scales. We also identify the characteristic length bellow which the viscous dissipation is strong and beyond which viscosity is essentially negligible. In contrast, magnetism seems to favor cosmic rotation. The magnetic presence is found to slow down the standard decay-rate of linear vortices, thus leading to universes with more residual rotation than generally anticipated.