We present a novel explanation of the well known $M_{\bullet}-\sigma$ relation. In a triaxial potential binaries with chaotic orbits within a sphere that encompass $\sim100$ times the mass of the super-massive black-hole (SMBH) have a finite probability to be tidally disrupted by the SMBH. As a result one component loses energy and might itself break apart tidally and accreted onto the SMBH. More significantly, the other component, which gains energy, returns to the bulge and equilibrates its excess energy with the environment thereby changing the kinetic temperature, hence the velocity dispersion. We develop a mathematical model and find that its results are in agreement with the observed relation.

Evolution of close binaries often proceeds through the common envelope stage. The physics of the envelope ejection (CEE) is not yet understood, and several mechanisms were suggested to be involved. These could give rise to different timescales for the CEE mass-loss. In order to probe the CEE-timescales we study wide companions to post-CE binaries. Faster mass-loss timescales give rise to higher disruption rates of wide binaries and result in larger average separations. We make use of data from Gaia DR2 to search for ultra-wide companions (projected separations $10^3$-$2\times 10^5$ a.u. and $M_2 > 0.4$ M$_\odot$) to several types of post-CEE systems, including sdBs, white-dwarf post-common binaries, and cataclysmic variables. We find a (wide-orbit) multiplicity fraction of $1.4\pm 0.2$ per cent for sdBs to be compared with a multiplicity fraction of $5.0\pm 0.2$ per cent for late-B/A/F stars which are possible sdB progenitors. The distribution of projected separations of ultra-wide pairs to main sequence stars and sdBs differs significantly and is compatible with prompt mass loss (upper limit on common envelope ejection timescale of $10^2$ years). The smaller statistics of ultra-wide companions to cataclysmic variables and post-CEE binaries provide weaker constraints. Nevertheless, the survival rate of ultra-wide pairs to the cataclysmic variables suggest much longer, $\sim10^4$ years timescales for the CEE in these systems, possibly suggesting non-dynamical CEE in this regime.

We present bounds on mirror dark matter scenario derived by using the effect of mirror matter on the luminosity of the Sun. In the perturbative regime where the mirror matter concentration is small relative to the ordinary matter we estimate the heat transfer from ordinary matter to the mirror sector by simple analytic consideration. That amount of heat transfer is radiated via mirror photons and increases the required energy production in order to maintain the observed luminosity. We then present more detailed numerical calculations of the total amount of this energy transfer.

Several scenarios were suggested for the origins of gravitational-wave (GW) sources from mergers of stellar binary black holes (BBHs). Here we propose a novel origin through catalyzed formation of GW-sources from ultra-wide binaries in the field. Such binaries experience perturbations from random stellar fly-bys which excite their eccentricities. Once a wide-binary is driven to a sufficiently small peri-center approach, GW-emission becomes significant, and the binary inspirals and merges. We derive an analytic model and verify it with numerical calculation to compute the merger rate to be $\sim10{\rm \times f_{wide}\,{\rm Gpc}^{-3}yr^{-1}}$ ($f_{{\rm wide}}$ is the fraction of wide BH-binaries), which is comparable to the observationally inferred rate. The observational signatures from this channel include spin-orbit misalignment; preference for high mass-ratio BBH; preference for high velocity-dispersion host-galaxies; and a uniform delay-time distribution.

Common envelope (CE) is an important phase in the evolution of interacting evolved binary systems. The interaction of the binary components during the CE evolution (CEE) stage gives rise to orbital inspiral and the formation of a short-period binary or a merger, on the expense of extending and/or ejecting the envelope. CEE is not well understood, as hydrodynamical simulations show that only a fraction of the CE-mass is ejected during the dynamical inspiral, in contrast with observations of post-CE binaries. Different CE models suggest different timescales are involved in the CE-ejection, and hence a measurement of the CE-ejection timescale could provide direct constraints on the CEE-process. Here we propose a novel method for constraining the mass-loss timescale of the CE, using post-CE binaries which are part of wide-orbit triple systems. The orbit/existence of a third companion constrains the CE mass-loss timescale, since rapid CE mass-loss may disrupt the triple system, while slower CE mass-loss may change the orbit of the third companion without disrupting it. As first test-cases we examine two observed post-CE binaries in wide triples, Wolf-1130 and GD-319. We follow their evolution due to mass-loss using analytic and numerical tools, and consider different mass-loss functions. We calculate a wide grid of binary parameters and mass-loss timescales in order to determine the most probable mass-loss timescale leading to the observed properties of the systems. We find that mass-loss timescales of the order of $10^{3}-10^{5}{\rm yr}$ are the most likely to explain these systems. Such long timescales are in tension with most of the CE mass-loss models, which predict shorter, dynamical timescales, but are potentially consistent with the longer timescales expected from the dust-driven winds model for CE ejection.

The mergers of neutron stars (NSs) and white dwarfs (WDs) could give rise to explosive transients, potentially observable with current and future transient surveys. However, the expected properties and distribution of such events is not well understood. Here we characterize the rates of such events, their delay time distribution, their progenitors and the distribution of their properties. We use binary populations synthesis models and consider a wide range of initial conditions and physical processes. In particular we consider different common-envelope evolution models and different NS natal kick distributions. We provide detailed predictions arising from each of the models considered. We find that the majority of NS-WD mergers are born in systems in which mass-transfer played an important role, and the WD formed before the NS. For the majority of the mergers the WDs have a carbon-oxygen composition (60-80%) and most of the rest are with oxygen-neon WDs. The rates of NS-WD mergers are in the range of 3-15% of the type Ia supernovae (SNe) rate. Their delay time distribution is very similar to that of type Ia SNe, but slightly biased towards earlier times. They typically explode in young 0.1-1Gyr environments, but have a tail distribution extending to long, Gyrs-timescales. Models including significant kicks give rise to relatively wide offset distribution extending to hundreds of kpcs. The demographic and physical properties of NS-WD mergers suggest they are likely to be peculiar type Ic-like SNe, mostly exploding in late type galaxies. Their overall properties could be related to a class of rapidly evolving SNe recently observed, while they are less likely to be related to the class of Ca-rich SNe.

Binary black-holes (BHs) and binary neutron-stars (NSs) mergers had been recently detected through gravitational-wave (GW) emission, with the latter followed by post-merger electromagnetic counterparts, appearing seconds up to weeks after the merger. While post-merger electromagnetic counterparts had been anticipated theoretically, very little electromagnetic precursors to GW-sources had been proposed, and non observed yet. Here we show that a fraction of ${\rm a}\,{\rm few\times}10^{-4}-10^{-1}$ of LIGO GW-sources and short-GRBs, could be preceded by supernovae-explosions years up to decades before the merger. Each of the BH/NS-progenitors in GW-sources are thought to form following a supernova, likely accompanied by a natal velocity-kick to the newly born compact object. The evolution and natal-kicks determine the orbits of surviving binaries, and hence the delay-time between the birth of the compact-binary and its final merger through GW-emission. We use data from binary evolution population-synthesis models to show that the delay-time distribution has a non-negligible tail of ultra-short delay-times between 1-100 yrs, thereby giving rise to potentially observable supernovae precursors to GW-sources. Moreover, future LISA/DECIGO GW space-detectors will enable the detection of GW-inspirals in the pre-mergers stage weeks to decades before the final merger. These ultra-short delay-time sources could therefore produce a unique type of promptly appearing LISA/DECIGO-GW-sources accompanied by \emph{coincident} supernovae. The archival (and/or direct) detection of precursor (coincident) SNe with GW and/or short-GRBs will provide unprecedented characterizations of the merging-binaries, and their prior evolution through supernovae and natal kicks, otherwise inaccessible through other means.

We present a numerical simulation of the dynamical interaction between the proposed Planet Nine and a debris disk around the Sun for $4{\rm Gyr}$, accounting for the secular perturbation of the four giant planets in two scenarios: (a) an initially thin circular disk around the Sun (b) inclined and eccentric disk. We show, in both scenarios, that Planet Nine governs the dynamics in between $1000{\rm -5000{\rm AU}}$ and forms spherical structure in the inner part ($\sim1000{\rm AU}$) and inclined disk. This structure is the outcome of mean motion resonances and secular interaction with Planet Nine. We compare the morphology of this structure with the outcome from a fly-by encounter of a star with the debris disk and show distinct differences between the two cases. We predict that this structure serves as a source of comets and calculate the resulting comet production rate to be detectable.

A fundamental aspect of the three-body problem is its stability. Most stability studies have focused on the co-planar three-body problem, deriving analytic criteria for the dynamical stability of such pro/retrograde systems. Numerical studies of inclined systems phenomenologically mapped their stability regions, but neither complement it by theoretical framework, nor provided satisfactory fit for their dependence on mutual inclinations. Here we present a novel approach to study the stability of hierarchical three-body systems at arbitrary inclinations, which accounts not only for the instantaneous stability of such systems, but also for the secular stability and evolution through Lidov-Kozai cycles and evection. We generalize the Hill-stability criteria to arbitrarily inclined triple systems, explain the existence of quasi-stable regimes and characterize the inclination dependence of their stability. We complement the analytic treatment with an extensive numerical study, to test our analytic results. We find excellent correspondence up to high inclinations $(\sim120^{\circ}$), beyond which the agreement is marginal. At such high inclinations the stability radius is larger, the ratio between the outer and inner periods becomes comparable, and our secular averaging approach is no longer strictly valid. We therefore combine our analytic results with polynomial fits to the numerical results to obtain a generalized stability formula for triple systems at arbitrary inclinations. Besides providing a generalized secular-based physical explanation for the stability of non co-planar systems, our results have direct implications for any triple systems, and in particular binary planets and moon/satellite systems; we briefly discuss the latter as a test case for our models.

Based on the observed high velocity of pulsars it is thought that neutron stars (NSs) receive a significant velocity kick at birth. Such natal kicks are considered to play an important role in the the evolution of binary-NS systems. The kick given to the NS (together with the effect of mass loss due to the supernova explosion of the NS progenitor) may result in the binary disruption or lead to a significant change of the binary orbital properties. Here we explore in detail the dynamical aftermath of natal kicks in binary systems, determine their possible outcomes and characterize their relative frequency, making use of analytic arguments and detailed population synthesis models. In a fraction of the cases the kick may cast the NS in such a trajectory as to collide with the binary companion, or pass sufficiently close to it as to disrupt it (micro tidal disruption event; $\mu$TDE), or alternatively it could be tidally-captured into a close orbit, eventually forming an X-ray binary. We calculate the rates of direct post-kick physical collisions and the possible potential production of Thorne-Zytkow objects or long-GRBs through this process, estimate the rates X-ray binaries formation and determine the rates of $\mu$TDEs and faint supernovae from white dwarf disruptions by NSs. Finally we suggest that natal kicks can produce BH-NS binaries with very short gravitational-wave merger time, possibly giving rise to a new type of promptly appearing eLISA gravitational wave (GW) sources, as well as producing aLIGO binary-merger GW sources with a unique (likely type Ib/c) supernova electromagnetic counterpart which precedes the GW merger.

We present a potentially efficient dynamical formation scenario for Low Mass X-ray Binaries (LMXBs) in the field, focusing on black-hole (BH) LMXBs. In this formation channel LMXBs are formed from wide binaries $(>1000$ AU) with a BH component and a stellar companion. The wide binary is perturbed by fly-by's of field stars and its orbit random-walks and changes over time. This diffusion process can drive the binary into a sufficiently eccentric orbit such that the binary components tidally interact at peri-center and the binary evolves to become a short period binary, which eventually evolves into an LMXB. The formation rate of LMXBs through this channel mostly depends on the number of such BH wide binaries progenitors, which in turn depends on the velocity kicks imparted to BHs (or NSs) at birth. We consider several models for the formation and survival of such wide binaries, and calculate the LMXB formation rates for each model. We find that models where BHs form through direct collapse with no/little natal kicks can give rise to high formation rates comparable with those inferred from observations. This formation scenario had several observational signatures: (1) the number density of LMXBs generally follows the background stellar density and (2) the mass function of the BH stellar companion should be comparable to the mass function of the background stellar population, likely peaking at $0.4-0.6$ M$_{\odot}$. The latter aspect, in particular, is unique to this model compared with previously suggested LMXB formation models following common envelope binary stellar evolution. We note that NS LMXBs can similarly form from wide binaries, but their formation rate through this channel is likely significantly smaller due to their much higher natal kicks.

The dwarf planet Pluto is known to host an extended system of five co-planar satellites. Previous studies have explored the formation and evolution of the system in isolation, neglecting perturbative effects by the Sun. Here we show that secular evolution due to the Sun can strongly affect the evolution of outer satellites and rings in the system, if such exist. Although precession due to extended gravitational potential from the inner Pluto-Charon binary quench such secular evolution up to $a_{crit}\sim0.0035$ AU ($\sim0.09$ $R_{Hill}$ the Hill radius; including all of the currently known satellites), outer orbits can be significantly altered. In particular, we find that \emph{co-planar} rings and satellites should not exist beyond $a_{crit}$; rather, satellites and dust particles in these regions secularly evolve on timescales ranging between $10^{4}-10^{6}$ yrs, and quasi-periodically change their inclinations and eccentricities through secular evolution (Lidov-Kozai oscillations). Such oscillations can lead to high inclinations and eccentricities, constraining the range where such satellites (and dust particles) can exist without crossing the orbits of the inner satellites, or crossing the outer Hill stability range. Outer satellites, if such exist are therefore likely to be \emph{irregular} satellites, with orbits limited to be non-circular and/or highly inclined. These could be potentially detected and probed by the New-Horizon mission, possibly providing direct evidence for the secular evolution of the Pluto satellite system, and shedding new light on its origins.

Recent studies have shown that secular evolution of triple systems can play a major role in the evolution and interaction of their inner binaries. Very few studies explored the stellar evolution of triple systems, and in particular the mass loss phase of the evolving stellar components. Here we study the dynamical secular evolution of hierarchical triple systems undergoing mass loss. We use the secular evolution equations and include the effects of mass-loss and mass-transfer, as well as general relativistic effects. We present various evolutionary channels taking place in such evolving triples, and discuss both the effects of mass-loss and mass-transfer in the inner binary system, as well as the effects of mass-loss/transfer from an outer third companion. We discuss several distinct types/regimes of triple secular evolution, where the specific behavior of a triple system can sensitively depend on its hierarchy and the relative importance of classical and general relativistic effects. We show that the orbital changes due to mass-loss and/or mass-transfer processes can effectively transfer a triple system from one dynamical regime to another. In particular, mass loss/transfer can both induce and quench high amplitude (Lidov-Kozai) variations in the eccentricity and inclination of the inner binaries of evolving triples. They can also change the system dynamics from an orderly periodic behavior to a chaotic one, and vice versa.