We develop stellar population models to predict the number of binaries with a single luminous member is Gaia and Hipparcos. Our models yield dozens of detections of black hole -luminous companion binaries (BHLC) and hundreds to thousands of neutron star -luminous companion binaries (NSLC) with Gaia. Interestingly, our models also yield a single detection of BHLC binary with Hipparcos, and a few NSLC binaries. We also show how the statistical distribution of detected binaries with a single luminous companion can be used to constrain the formation process of neutron stars and black holes.

We present here RICH, a state-of-the-art two-dimensional hydrodynamic code based on Godunov's method, on an unstructured moving mesh (the acronym stands for Racah Institute Computational Hydrodynamics). This code is largely based on the code AREPO. It differs from AREPO in the interpolation and time-advancement schemeS as well as a novel parallelization scheme based on Voronoi tessellation. Using our code, we study the pros and cons of a moving mesh (in comparison to a static mesh). We also compare its accuracy to other codes. Specifically, we show that our implementation of external sources and time-advancement scheme is more accurate and robust than is AREPO when the mesh is allowed to move. We performed a parameter study of the cell rounding mechanism (Lloyd iterations) and its effects. We find that in most cases a moving mesh gives better results than a static mesh, but it is not universally true. In the case where matter moves in one wayand a sound wave is traveling in the other way (such that relative to the grid the wave is not moving) a static mesh gives better results than a moving mesh. We perform an analytic analysis for finite difference schemes that reveals that a Lagrangian simulation is better than aEulerian simulation in the case of a highly supersonic flow. Moreover, we show that Voronoi-based moving mesh schemes suffer from an error, which is resolution independent, due to inconsistencies between the flux calculation and the change in the area of a cell. Our code is publicly available as open source and designed in an object-oriented, user-friendly way that facilitates incorporation of new algorithms and physical processes.

When a star gets too close to a supermassive black hole, it is torn apart by the tidal forces. Roughly half of the stellar mass becomes unbound and flies away at tremendous velocities -around 10 4 km/s. In this work we explore the idea that the shock produced by the interaction of the unbound debris with the ambient medium gives rise to the synchrotron radio emission observed in several TDEs. We use a moving mesh numerical simulation to study the evolution of the unbound debris and the bow shock around it. We find that as the periapse distance of the star decreases, the outflow becomes faster and wider. A tidal disruption event whose periapse distance is a factor of 7 smaller than the tidal radius can account for the radio emission observed in ASASSN-14li. This model also allows us to obtain a more accurate estimate for the gas density around the centre of the host galaxy of ASASSN-14li.

A linear stability analysis is performed for two counterstreaming proton beams in the presence of a thermal electron background. Growth rates and polarization properties of unstable modes are calculated for various density ratios of the proton beams. It is found that in most cases, two unstable modes grow simultaneously: an electromagnetic filamentary mode that propagates perpendicular to the beam and an electrostatic mode that propagates parallel to the beam. The growth rates of the two modes are comparable, so that one expects that the instability would result in the development of a filamentary structure with a superimposed electrostatic pattern.

Sub-Chandrasekhar CO white dwarfs accreting helium have been considered as candidates for Type Ia supernova(SNIa) progenitors since the early 1980s (helium shell mass > 0.1M ). These models, once detonated did not fit the observed spectra and light curve of typical SNIa observations. New theoretical work examined detonations on much less massive (< 0.05M ) envelopes. They find stable detonations that lead to light curves, spectra and abundances that compare relatively well with the observational data. The exact mechanism leading to the ignition of helium detonation is a key issue, since it is a mandatory first step for the whole scenario. As the flow of the accreted envelope is unstable to convection long before any hydrodynamic phenomena develops, a multidimensional approach is needed in order to study the ignition process. The complex convective reactive flow is challenging to any hydrodynamical solver. To the best of our knowledge, all previous 2D studies ignited the detonation artificially. We present here, for the first time, fully consistent results from two hydrodynamical 2D solvers that adopt two independent accurate schemes. For both solvers an effort was made to overcome the problematics raised by the finite resolution and numerical diffusion by the advective terms. Our best models lead to the ignition of a detonation in a convective cell. Our results are robust and the agreement between the two different numerical approaches is very good.