$N$-body simulations of non-resonant tightly-packed planetary systems have found that their survival time (i.e. time to first close encounter) grows exponentially with their interplanetary spacing and planetary masses. Although this result has important consequences for the assembly of planetary systems by giants collisions and their long-term evolution, this underlying exponential dependence is not understood from first principles, and previous attempts based on orbital diffusion have only yielded power-law scalings. We propose a different picture, where large deviations of the system from its initial conditions is due to few slowly developing high-order resonances. Thus, we show that the survival time of the system $T$ can be estimated using a heuristic motivated by Nekhoroshev's theorem, and obtain a formula for systems away from overlapping two-body mean-motion resonances as: $T/P=c_1 \frac{a}{\Delta a} \exp \left(c_2 \frac{\Delta a}{a} /\mu^{1/4}\right)$, where $P$ is the average Keplerian period, $a$ is the average semi major axis, $\Delta a\ll a$ is the difference between the semi major axes of neighbouring planets, $\mu$ is the planet to star mass ratio, and $c_1$ and $c_2$ are dimensionless constants. We show that this formula is in good agreement with numerical N-body experiments for $c_1=5 \cdot 10^{-4}$ and $c_2=8$.

We present a mathematical model for the propagation of the shock waves that occur during planetary collisions. Such collisions are thought to occur during the formation of terrestrial planets, and they have the potential to erode the planet's atmosphere. We show that under certain assumptions, this evolution of the shock wave can be determined using the method of self similar solutions. This self similar solution is of type II, which means that it only applies to a finite region behind the shock front. This region is bounded by the shock front and the sonic point. Energy and matter continuously flow through the sonic point, so that energy in the self similar region is not conserved, as is the case for type I solutions. Instead, the evolution of the shock wave is determined by boundary conditions at the shock front and at the sonic point. We show how the evolution can be determined for different equations of state, allowing these results to be readily used to calculate the atmospheric mass loss from planetary cores made of different materials.

We study the hydrodynamic evolution of an explosion close to the stellar surface, and give predictions for the radiation from such an event. We show that such an event will give rise to a multi-wavelength transient. We apply this model to describe a precursor burst to the peculiar supernova iPTF14hls, which occurred in 1954, sixty year before the supernova. We propose that the new generation of optical surveys might detect similar transients, and they can be used to identify supernova progenitors well before the explosion.

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.

Using moving mesh hydrodynamic simulations, we determine the shock propagation and resulting ground velocities for a planet hit by a high velocity impactor. We use our results to determine the atmospheric mass loss caused by the resulting ground motion due to the impact shock wave. We find that there are two distinct shock propagation regimes: In the limit in which the impactor is significantly smaller than the target ($R_i<< R_t$), the solutions are self-similar and the shock velocity at a fixed point on the target scale as $m_i^{2/3}$, where $m_i$ is the mass of the impactor. In addition, the ground velocities follow a universal profile given by $v_g/v_i=(14.2x^2-25.3x+11.3)/(x^2-2.5x+1.9) +2\ln{R_i/R_t}$, where $x=\sin\left(\theta/2\right)$, $\theta$ is the latitude on the target measured from the impact site, and $v_g$ and $v_i$ are the ground velocity and impact velocity, respectively. In contrast, in the limit in which the impactor is comparable to the size of the target ($R_i \sim R_t$), we find that shock velocities decline with the mass of the impactor significantly more weakly than $m_i^{2/3}$. We use the resulting surface velocity profiles to calculate the atmospheric mass loss for a large range of impactor masses and impact velocities and apply them to the Kepler-36 system and the Moon forming impact. Finally, we present and generalise our results in terms of the $v_g/v_i$ and the impactor to target size ratio ($R_i/R_t$) such that they can easily be applied to other collision scenarios.

Fast radio bursts (FRBs) are, as the name implies, short and intense pulses of radiation at wavelengths of roughly one metre. FRBs have extremely high brightness temperatures, which points to a coherent source of radiation. The energy of a single burst ranges from $10^{36}$ to $10^{39}$ erg. At the high end of the energy range, FRBs have enough energy to unbind an earth-sized planet, and even at the low end, there is enough energy to vaporise and unbind the atmosphere and the oceans. We therefore propose that FRBs are signatures of an artificial terraformer, capable of eradicating life on another planet, or even destroy the planet entirely. The necessary energy can be harvested from Wolf-Rayet stars with a Dyson sphere ($\sim 10^{38}$ erg s$^{-1}$) , and the radiation can be readily produced by astrophysical masers. We refer to this mechanism as Volatile Amplification of a Destructive Emission of Radiation (VADER). We use the observational information to constrain the properties of the apparatus. We speculate that the non-repeating FRBs are low-energy pulses used to exterminate life on a single planet, but leaving it otherwise intact, and that the stronger repeating FRB is part of an effort to destroy multiple objects in the same solar system, perhaps as a preventative measure against panspermia. In this picture, the persistent synchrotron source associated with the first repeating FRB arises from the energy harvesting process. Finally we propose that Oumuamua might have resulted from a destruction of a planet in this manner.

The distances to fast radio bursts (FRBs) are crucial for understanding their underlying engine, and for their use as cosmological probes. In this paper, we provide three statistical estimates of the distance to ASKAP FRBs. First, we show that the number of events of similar luminosity in ASKAP does not scale as distance cubed, as one would expect, when directly using the observed dispersion measure (DM) to infer distance. Second, by comparing the average DMs of FRBs observed with different instruments, we estimated the average redshift of ASKAP FRBs to be $z\sim 0.01$ using CHIME and ASKAP, and $z\lesssim0.07$ using Parkes and ASKAP. Both values are much smaller than the upper limit $z\sim0.3$ estimated directly from the DM. Third, we cross-correlate the locations of the ASKAP FRBs with existing large-area redshift surveys, and see a 3$\sigma$ correlation with the 2MASS Redshift Survey and a 5$\sigma$ correlation with the HI Parkes All Sky Survey at $z\sim0.007$. This corresponds well with the redshift of the most likely host galaxy of ASKAP FRB 171020, which is at $z=0.00867$. These arguments combined suggest an extremely nearby origin of ASKAP FRBs and a local environment with accumulated electrons that contribute a DM of several hundred pc/cm$^3$, which should be accounted for in theoretical models.

Type II supernova progenitors are expected to emit copious amounts of mass in a dense stellar wind prior to the explosion. When the progenitor is a member of a binary, the orbital motion modulates the density of this wind. When the progenitor explodes, the high-velocity ejecta collides with the modulated wind, which in turn produces a modulated radio signal. In this work we derive general analytic relations between the parameters of the radio signal modulations and binary parameter in the limit of large member mass ratio. We use these relations to infer the semi major axis of SN1979c and a lower bound for the mass of the companion. We further constrain the analytic estimates by numerical simulations using the AMUSE framework. In these calculations we simulate the progenitor binary system including the wind and the gravitational effect of a companion star. The simulation output is compared to the observed radio signal in supernova SN1979C. We find that it must have been a binary with an orbital period of about 2000 year. If the exploding star evolved from a $\sim 18 M_{\odot}$ zero-age main-sequence at solar metalicity, we derive a companion mass of $5$ to $12 M_{\odot}$ in an orbit with an eccentricity lower than about 0.8.

Tidal disruption events of stars by supermassive black holes have so far been discovered months to years after the fact. In this paper we explore the short, faint and hard burst of radiation is emitted at maximum compression, as a result of shock breakout. The detection of this burst can be used to capture tidal disruption events in real time. We verify that shock breakout from main sequence stars produces radiation in the X-ray range, but find that it is difficult to detect using all sky X-ray surveying telescopes. In the case of shock breakout from red giants, most of the radiation is emitted in the UV and visible range, which is significantly easier to detect. A similar burst of UV/optical radiation will also be emitted by stars puffed by tidal heating from a previous passage close to the central black hole. This radiation can be detected by surveys like ZTF and LSST. We calculate detection rates for different types of galactic nuclei. For the case of a very full loss cone we predict a detection rate of once per month with LSST, whereas for the case of a very empty loss cone we predict a rate of once per year with LSST. Evidence from a recent tidal disruption event, ASASSN-14li, seems to favour a very full loss cone, in which case LSST is expected to detect one such event every month.

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 develop an analytic, steady-state model for the gas environment in quiescent galactic nuclei. We assume that the mass is constantly supplied by a spherically symmetric distribution of wind emitting stars, and that gravity is solely due to a central supermassive black hole. We show that at some finite radius, where the Keplerian velocity is comparable to the wind velocity, the bulk velocity vanishes. Matter generated below that radius will be accreted onto the black hole, while matter outside it will escape the system. Under certain conditions, the flow may become supersonic at both domains. We obtain radial profiles of the hydrodynamic variables and verify them using a time-dependent hydrodynamic simulation. We delineate the conditions under which radiative cooling can be neglected, and predict the luminosity and spectrum of the free-free X-ray emission from such a system. We discuss applications of our solution to our own Galactic Centre and other quiescent galactic nuclei.

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.1 M_\odot $). 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.05 M_\odot $) 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.

It has been proposed that an intermediate mass black hole can produce compact molecular clouds with large velocity dispersions. This model is called into question due to the discovery of several such short lived clouds, which suggests they are common, and hence it is unlikely that their formation involves an exotic object such as an intermediate mass black hole. In this paper we propose an alternative model, where supernovae produce such clouds. We apply our model to a such a cloud called CO-0.40-0.22. A previous work pointed out the fact that the kinetic energy of this cloud is comparable to that from a supernova, but disqualified the supernova scenario due to lack of thermal emission. We explain the lack of thermal emission as a consequence of the large density of the medium. In such a medium, the supernova shock is radiation dominated, and hence the initial temperatures are lower by many orders of magnitude compared with a typical supernova. Finally, we discuss some of the consequences of having a significant population of intermediate mass black holes in the galactic centre.

Supernovae near the galactic center evolve differently from regular galactic supernovae. This is mainly due to the environment into which the supernova remnants propagate. Instead of a static, uniform density medium, SNRs near the galactic center propagate into a wind swept environment with a velocity away from the galactic center, and a graded density profile. This causes these SNRs to be non - spherical, and to evolve faster than their galactic counterparts. We develop an analytic theory for the evolution of explosions within a stellar wind, and verify it using a hydrodynamic code. We show that such explosions can evolve in one of three possible morphologies. Using these results we discuss the association between the two SNRs (SGR East and SGR A's bipolar radio/X-ray Lobes) and the two neutron stars (the cannonball and SGR J1745-2900) near the galactic center. We show that, given the morphologies of the SNR and positions of the neutron stars, the only possible association is between SGR A's bipolar radio/X-ray Lobes and SGR J1745-2900. The compact object created in the explosion of SGR East remains undetected, and the SNR of the supernova that created the cannonball has already disappeared.

We investigate the relativistic breakout of a shock wave from the surface of a star. In this process, each fluid shell is endowed with some kinetic and thermal energy by the shock, and then continues to accelerate adiabatically by converting thermal energy into kinetic energy. This problem has been previously studied for a mildly relativistic breakout, where the acceleration ends close to the surface of the star. The current work focuses on the case where the acceleration ends at distances much greater than the radius of the star. We derive an analytic description for the hydrodynamic evolution of the ejecta in this regime, and validate it using a numerical simulation. We also provide predictions for the expected light curves and spectra from such an explosion. The relevance to astrophysical explosions is discussed, and it is shown that such events require more energy than is currently believed to result from astrophysical explosions.

Current Voronoi based moving mesh hydro codes suffer from "grid noise". We identify the cause of this noise as the volume inconsistency error, where the volume that is transferred between cells is inconsistent with the hydrodynamical calculations. As a result, the codes do not achieve second order convergence. In this paper we describe how a simple fix allows Voronoi based moving mesh codes to attain second order convergence. The fix is based on the understanding that the volume exchanged between cells should be consistent with the hydrodynamical calculations. We benchmark our fix with three test problems and show that it can significantly improve the computational accuracy. We also examine the effect of initial mesh initialization and present an improved model for the Green-Gauss based gradient estimator.

The problem of a cold gas flowing past a stationary object is considered. It is shown that at large distances from the obstacle the shock front forms a parabolic solid of revolution. The interior of the shock front is obtained by solution of the hydrodynamic equations in parabolic coordinates. The results are verified with a hydrodynamic simulation. The drag force and expected spectra are calculated for such shock, both in case of an optically thin and thick media. Finally, relations to astrophysical bow shocks and other analytic works on oblique shocks are discussed.

The use of numerical simulations in science is ever increasing and with it the computational size. In many cases single processors are no longer adequate and simulations are run on multiple core machines or supercomputers. One of the key issues when running a simulation on multiple CPUs is maintaining a proper load balance throughout the run and minimizing communications between CPUs. We propose a novel method of utilizing a Voronoi diagram to achieve a nearly perfect load balance without the need of any global redistributions of data. As a show case, we implement our method in RICH, a 2D moving mesh hydrodynamical code, but it can be extended trivially to other codes in 2D or 3D. Our tests show that this method is indeed efficient and can be used in a large variety of existing hydrodynamical codes as well as other applications.

We present here RICH, a state of the art 2D 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 scheme 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 AREPO's, when the mesh is allowed to move. We performed a parameter study of the cell rounding mechanism (Llyod iterations) and it 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. Moreover, we show that Voronoi based moving mesh schemes suffer from an error, that is resolution independent, due to inconsistencies between the flux calculation and 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.

We quantify the atmospheric mass loss during planet formation by examining the contributions to atmospheric loss from both giant impacts and planetesimal accretion. Giant impacts cause global motion of the ground. Using analytic self-similar solutions and full numerical integrations we find (for isothermal atmospheres with adiabatic index ($\gamma=5/3$) that the local atmospheric mass loss fraction for ground velocities $v_g < 0.25 v_{esc}$ is given by $\chi_{loss}=(1.71 v_g/v_{esc})^{4.9}$, where $v_{esc}$ is the escape velocity from the target. Yet, the global atmospheric mass loss is a weaker function of the impactor velocity $v_{Imp}$ and mass $m_{Imp}$ and given by $X_{loss} ~ 0.4x+1.4x^2-0.8x^3$ (isothermal atmosphere) and $X_{loss} ~ 0.4x+1.8x^2-1.2x^3$ (adiabatic atmosphere), where $x=(v_{Imp}m/v_{esc}M)$. Atmospheric mass loss due to planetesimal impacts proceeds in two different regimes: 1) Large enough impactors $m > \sqrt{2} \rho_0 (\pi h R)^{3/2}$ (25~km for the current Earth), are able to eject all the atmosphere above the tangent plane of the impact site, which is $h/2R$ of the whole atmosphere, where $h$, $R$ and $\rho_0$ are the atmospheric scale height, radius of the target, and its atmospheric density at the ground. 2) Smaller impactors, but above $m>4 \pi \rho_0 h^3$ (1~km for the current Earth) are only able to eject a fraction of the atmospheric mass above the tangent plane. We find that the most efficient impactors (per unit impactor mass) for atmospheric loss are planetesimals just above that lower limit and that the current atmosphere of the Earth could have resulted from an equilibrium between atmospheric erosion and volatile delivery to the atmosphere from planetesimals. We conclude that planetesimal impacts are likely to have played a major role in atmospheric mass loss over the formation history of the terrestrial planets. (Abridged)

In this paper we extend the Sedov - Taylor - Von Neumann model for a strong explosion to account for small angular and radial variations in the density. We assume that the density profile is given by $\rho\left(r,\theta,\phi\right)=kr^{-\omega}\left(1+\varepsilon\left(\frac{r}{r_{0}}\right)^{q}Y_{lm}\left(\theta,\phi\right)\right)$, where $\varepsilon\ll1$ and $\omega\le\frac{7-\gamma}{\gamma+1}$. In order to verify our results we compare them to analytical approximations and full hydrodynamic simulations. We demonstrate how this method can be used to describe arbitrary (not just self similar) angular perturbations. This work complements our previous analysis on radial, spherically symmetric perturbations, and allows one to calculate the response of an explosion to arbitrary perturbations in the upstream density. Together, they settle an age old controversy about the inner boundary conditions

We present new solutions to the strong explosion problem in a non power law density profi{}le. The unperturbed self similar solutions developed by Sedov, Taylor and Von Neumann describe strong Newtonian shocks propagating into a cold gas with a density profile falling off as $r^{-\omega}$, where $\omega\le\frac{7-\gamma}{\gamma+1}$ (filled type I solutions), and $\gamma$ is the adiabatic index of the gas. The perturbations we consider are spherically symmetric and log periodic with respect to the radius. While the unperturbed solutions are continuously self similar, the log periodicity of the density perturbations leads to a discrete self similarity of the perturbations, i.e., the solution repeats itself up to a scaling at discrete time intervals. We discuss these solutions and verify them against numerical integrations of the time dependent hydrodynamic equations. This is an extension of a previous investigation on type II solutions and helps clarifying boundary conditions for perturbations to type I self similar solutions.

Self similarity allows for analytic or semi-analytic solutions to many hydrodynamics problems. Most of these solutions are one dimensional. Using linear perturbation theory, expanded around such a one-dimensional solution, we find self-similar hydrodynamic solutions that are two- or three-dimensional. Since the deviation from a one-dimensional solution is small, we call these slightly two-dimensional and slightly three-dimensional self-similar solutions, respectively. As an example, we treat strong spherical explosions of the second type. A strong explosion propagates into an ideal gas with negligible temperature and density profile of the form rho(r,theta,phi)=r^{-omega}[1+sigma*F(theta,phi)], where omega>3 and sigma << 1. Analytical solutions are obtained by expanding the arbitrary function F(theta,phi) in spherical harmonics. We compare our results with two dimensional numerical simulations, and find good agreement.