Classical novae are thermonuclear explosions that occur on the surfaces of white dwarf stars in interacting binary systems (Bode & Evans 2008). It has long been thought that the luminosity of classical novae is powered by continued nuclear burning on the surface of the white dwarf after the initial runaway (Gallaher & Starrfield 1978). However, recent observations of GeV $\gamma$-rays from classical novae have hinted that shocks internal to the nova ejecta may dominate the nova emission. Shocks have also been suggested to power the luminosity of events as diverse as stellar mergers (Metzger & Pejcha 2017), supernovae (Moriya et al. 2018), and tidal disruption events (Roth et al. 2016), but observational confirmation has been lacking. Here we report simultaneous space-based optical and $\gamma$-ray observations of the 2018 nova V906 Carinae (ASASSN-18fv), revealing a remarkable series of distinct correlated flares in both bands. The optical and $\gamma$-ray flares occur simultaneously, implying a common origin in shocks. During the flares, the nova luminosity doubles, implying that the bulk of the luminosity is shock-powered. Furthermore, we detect concurrent but weak X-ray emission from deeply embedded shocks, confirming that the shock power does not appear in the X-ray band and supporting its emergence at longer wavelengths. Our data, spanning the spectrum from radio to $\gamma$-ray, provide direct evidence that shocks can power substantial luminosity in classical novae and other optical transients.

We present one-dimensional hydrodynamical simulations including radiative losses, of internal shocks in the outflows from classical novae, to explore the role of shocks in powering multi-wavelength emission from radio to gamma-ray wavelengths. Observations support a picture in which the initial phases of some novae generate a slow, equatorially-focused outflow (directly from the outer Lagrange point, or from a circumbinary disk), which then transitions to, or is overtaken by, a faster more isotropic outflow from the white dwarf which collides and shocks the slower flow, powering gamma-ray and optical emission through reprocessing by the ejecta. However, the common occurence of multiple peaks in nova light curves suggests that the outflow's acceleration need not be monotonic, but instead can involve successive transitions between "fast" and "slow" modes. Such a time-fluctuating outflow velocity naturally can reproduce several observed properties of nova, such as correlated gamma-ray and optical flares, expansion of the photosphere coincident with (though lagging slightly) the peak flare luminosity, and complex time-evolution of spectral lines (including accelerating, decelerating, and merging velocity components). While the shocks are still deeply embedded during the gamma-ray emission, the onset of $\sim$ keV X-ray and $\sim 10$ GHz radio synchrotron emission is typically delayed until the forward shock of the outermost monolithic shell (created by merger of multiple internal shock-generated shells) reaches a sufficiently low column through the dense external medium generated by the earliest phase of the outburst.

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.

The tidal destruction of a star by a massive black hole, known as a tidal disruption event (TDE), is commonly modeled using the "frozen-in" approximation. Under this approximation, the star maintains exact hydrostatic balance prior to entering the tidal sphere (radius $r_{\rm t}$), after which point its internal pressure and self-gravity become instantaneously negligible and the debris undergoes ballistic free fall. We present a suite of hydrodynamical simulations of TDEs with high penetration factors $\beta \equiv r_{\rm t}/r_{\rm p} = 5-7$, where $r_{\rm p}$ is the pericenter of the stellar center of mass, calculated using a Voronoi-based moving-mesh technique. We show that basic assumptions of the frozen-in model, such as the neglect of self-gravity inside $r_{\rm t}$, are violated. Indeed, roughly equal fractions of the final energy spread accumulate exiting and entering the tidal sphere, though the frozen-in prediction is correct at the order-of-magnitude level. We also show that an $\mathcal{O}(1)$ fraction of the debris mass remains transversely confined by self-gravity even for large $\beta$ which has implications for the radio emission from the unbound debris and, potentially, for the circularization efficiency of the bound streams.

The discovery of GeV gamma-rays from classical novae has led to a reassessment of these garden-variety explosions, and highlighted their importance for understanding radiative shocks, particle acceleration, and dust formation in more exotic, distant transients. Recent collaboration between observers and theorists has revealed that shocks in novae are energetically important, and can even dominate their bolometric luminosity. Shocks may also explain long-standing mysteries in novae such as dust production, super-Eddington luminosities, and `flares' in optical light curves. Here, we highlight the multi-wavelength facilities of the next decade that will further test our nova shock model and fulfill the promise of novae as powerful astrophysical laboratories.

Radiative shocks, behind which gas cools faster than the dynamical time, play a key role in many astrophysical transients, including classical novae and young supernovae interacting with circumstellar material. The dense layer behind high Mach number $\mathcal{M} \gg 1$ radiative shocks is susceptible to thin-shell instabilities, creating a "corrugated" shock interface. We present two and three-dimensional hydrodynamical simulations of optically-thin radiative shocks to study their thermal radiation and acceleration of non-thermal relativistic ions. We employ a moving-mesh code and a specialized numerical technique to eliminate artificial heat conduction across grid cells. The fraction of the shock's luminosity $L_{\rm tot}$ radiated at X-ray temperatures $kT_{\rm sh} \approx (3/16)\mu m_p v_{\rm sh}^{2}$ expected from a one-dimensional analysis is suppressed by a factor $L(>T_{\rm sh}/3)/L_{\rm tot} \approx 4.5/\mathcal{M}^{4/3}$ for $\mathcal{M} \approx 4-36$. This suppression results in part from weak shocks driven into under-pressured cold filaments by hot shocked gas, which sap thermal energy from the latter faster than it is radiated. Combining particle-in-cell simulation results for diffusive shock acceleration with the inclination angle distribution across the shock (relative to an upstream magnetic field in the shock plane$-$the expected geometry for transient outflows), we predict the efficiency and energy spectrum of ion acceleration. Though negligible acceleration is predicted for adiabatic shocks, the corrugated shock front enables local regions to satisfy the quasi-parallel magnetic field geometry required for efficient acceleration, resulting in an average acceleration efficiency of $\epsilon_{\rm nth} \sim 0.005-0.02$ for $\mathcal{M} \approx 12-36$, in agreement with modeling of the gamma-ray nova ASASSN-16ma.

As part of our long-term campaign to understand how cold streams feed massive galaxies at high redshift, we study the Kelvin-Helmholtz instability (KHI) of a supersonic, cold, dense gas stream as it penetrates through a hot, dilute circumgalactic medium (CGM). A linear analysis (Paper I) showed that, for realistic conditions, KHI may produce nonlinear perturbations to the stream during infall. Therefore, we proceed here to study the nonlinear stage of KHI, still limited to a two-dimensional slab with no radiative cooling or gravity. Using analytic models and numerical simulations, we examine stream breakup, deceleration and heating via surface modes and body modes. The relevant parameters are the density contrast between stream and CGM ($\delta$), the Mach number of the stream velocity with respect to the CGM ($M_{\rm b}$) and the stream radius relative to the halo virial radius ($R_{\rm s}/R_{\rm v}$). We find that sufficiently thin streams disintegrate prior to reaching the central galaxy. The condition for breakup ranges from $R_{\rm s} < 0.03 R_{\rm v}$ for $(M_{\rm b} \sim 0.75, \delta \sim 10)$ to $R_{\rm s} < 0.003 R_{\rm v}$ for $(M_{\rm b} \sim 2.25, \delta \sim 100)$. However, due to the large stream inertia, KHI has only a small effect on the stream inflow rate and a small contribution to heating and subsequent Lyman-$\alpha$ cooling emission.

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.

The G2 object has recently passed its pericenter passage in our Galactic Center. While the $Br_\gamma$ emission shows clear signs of tidal interaction, the change in the observed luminosity is only of about a factor of 2, in contention with all previous predictions. We present high resolution simulations performed with the moving mesh code, RICH, together with simple analytical arguments that reproduce the observed $Br_\gamma$ emission. In our model, G2 is a gas cloud that undergoes tidal disruption in a dilute ambient medium. We find that during pericenter passage, the efficient cooling of the cloud results in a vertical collapse, compressing the cloud by a factor of $\sim5000$. By properly taking into account the ionization state of the gas, we find that the cloud is UV starved and are able to reproduce the observed $Br_\gamma$ luminosity. For densities larger than $\approx500\;\mathrm{cm}^{-3}$ at pericenter, the cloud fragments, due to cooling instabilities and the emitted radiation is inconsistent with observations. For lower densities, the cloud survives the pericenter passage intact and its emitted radiation matches the observed lightcurve. From the duration of $Br_\gamma$ emission which contains both redshifted and blueshifted components, we show that the cloud is not spherical but rather elongated with a size ratio of 4 at year 2001. The simulated cloud's elongation grows as it travels towards pericenter and is consistent with observations, due to viewing angles. The simulation is also consistent with having a spherical shape at apocenter.

Massive galaxies at high redshift are predicted to be fed from the cosmic web by narrow, dense, cold streams. These streams penetrate supersonically through the hot medium encompassed by a stable shock near the virial radius of the dark-matter halo. Our long-term goal is to explore the heating and dissipation rate of the streams and their fragmentation and possible breakup, in order to understand how galaxies are fed, and how this affects their star-formation rate and morphology. We present here the first step, where we analyze the linear Kelvin-Helmholtz instability (KHI) of a cold, dense slab or cylinder flowing through a hot, dilute medium in the transonic regime. The current analysis is limited to the adiabatic case with no gravity and assuming equal pressure in the stream and the medium. By analytically solving the linear dispersion relation, we find a transition from a dominance of the familiar rapidly growing surface modes in the subsonic regime to more slowly growing body modes in the supersonic regime. The system is parameterized by three parameters: the density contrast between the stream and the medium, the Mach number of stream velocity with respect to the medium, and the stream width with respect to the halo virial radius. We find that a realistic choice for these parameters places the streams near the mode transition, with the KHI exponential-growth time in the range 0.01-10 virial crossing times for a perturbation wavelength comparable to the stream width. We confirm our analytic predictions with idealized hydrodynamical simulations. Our linear-KHI estimates thus indicate that KHI may in principle be effective in the evolution of streams by the time they reach the galaxy. More definite conclusions await the extension of the analysis to the nonlinear regime and the inclusion of cooling, thermal conduction, the halo potential well, self-gravity and magnetic fields.

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 Asteroid Belt and the Kuiper Belt are relics from the formation of our solar system. Understanding the size and spin distribution of the two belts is crucial for a deeper understanding of the formation of our solar system and the dynamical process that govern it. In this paper, we investigate the effect of collisions on the evolution of the spin distribution of asteroids and KBO's. We find that the power law nature of the impactors' size distribution leads to a L\'evy distribution of the spin rates. This results in a power law tail of the spin distribution, in stark contrast to the usually quoted Maxwellian distribution. We show that for bodies larger than 10 km, collisions alone lead to spin rates peaking at 0.15-0.5 revolutions per day. Comparing that to the observed spin rates of large asteroids ($R>50$ km), we find that the spins of large asteroids, peaking at $\sim1-2$ revolutions per day, are dominated by a primordial component that reflects the formation mechanism of the asteroids. Similarly, the Kuiper Belt has undergone virtually no collisional spin evolution, assuming current day density. Collisions contribute a spin rate of $\sim0.01$ revolutions per day, thus the observed fast spin rates of KBOs are also primordial in nature.

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 new, deep near-infrared SINFONI @ VLT integral field spectroscopy of the gas cloud G2 in the Galactic Center, from late summer 2013 and spring 2014. G2 is visible in recombination line emission. The spatially resolved kinematic data track the ongoing tidal disruption. As expected for an observation near pericenter passage, roughly half of the gas in 2014 is found at the redshifted, pre-pericenter side of the orbit, while the other half is at the post-pericenter, blueshifted side. We also present an orbital solution for the gas cloud G1, which was discovered a decade ago in L'-band images when it was spatially almost coincident with Sgr A*. The orientation of the G1 orbit in the three angles is almost identical to the one of G2, but it has a lower eccentricity and smaller semi-major axis. We show that the observed astrometric positions and radial velocities of G1 are compatible with the G2 orbit, assuming that (i) G1 was originally on the G2 orbit preceding G2 by 13 years and (ii) a simple drag force acted on it during pericenter passage. Taken together with the previously described tail of G2, which we detect in recombination line emission and thermal broadband emission, we propose that G2 may be a bright knot in a much more extensive gas streamer. This matches purely gaseous models for G2, such as a stellar wind clump or the tidal debris from a partial disruption of a star.

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.

The Kepler mission has discovered a large number of planetary systems. We analyze the implications of the discovered single/multi-exoplanet systems from Kepler's data. As done in previous works, we test a simple model in which the intrinsic occurrence of planet is an independent process, and with equal probability around all planet producing stars. This leads to a Poisson distribution for the intrinsic number of planets around each host. However, the possibility of zero/low mutual inclination is taken into account, creating a correlation between detecting different planets in a given stellar system, leading to a non Poisson distribution for the number of transiting planets per system. Comparing the model's predictions with the observations made by Kepler, we find that the correlation produced by planarity is insufficient and a higher correlation is needed; either the formation of one planet in the system enhances the likelihood of other planets to form, and/or that some stars are considerably more fertile than others. Kepler's data presents evidences that both correlations might play a part, in particular a significant dependency in the radial distribution of planets in multi-planet systems is shown. Followup observations on Kepler planet's hosts can help pinpoint the physical nature of this correlation.

The rotation states of kilometer sized near earth asteroids are known to be affected by the YORP effect. In a related effect, Binary YORP (BYORP) the orbital properties of a binary asteroid evolves under a radiation effect mostly acting on a tidally locked secondary. The BYORP effect can alter the orbital elements in $\sim 10^{4-5}$ years for a $D_{p}=2\;km$ primary with a $D_{s}=0.4\; km$ secondary at $1\; AU$. It can either separate the binary components or cause them to collide. In this paper we devise a simple approach to calculate the YORP effect on asteroids and BYORP effect on binaries including $J_2$ effects due to primary oblateness and the sun. We apply this to asteroids with known shapes as well as a set of randomly generated bodies with various degrees of smoothness. We find a strong correlation between the strengths of an asteroid's YORP and BYORP effects. Therefore, a statistical knowledge on one, could be used to estimate the effect of the other. We show that the action of BYORP preferentially shrinks rather than expands the binary orbit and that YORP preferentially slows down asteroids. This conclusion holds for the two extremes of thermal conductivities studied in this work and assuming the asteroid reaches a stable point, but may break down for moderate thermal conductivity. The YORP and BYORP effects are shown to be smaller than what could be naively expected due to near cancellation of the effects on small scales. Taking this near cancellation into account, a simple order of magnitude estimate of the YORP and BYORP effects as function of the sizes and smoothness of the bodies is calculated. Finally, we provide a simple proof showing that there is no secular effect due to absorption of radiation in BYORP.