We study the violent phase of the merger of massive binary white dwarf systems. Our aim is to characterize the conditions for explosive burning to occur, and identify a possible explosion mechanism of Type Ia supernovae. The primary components of our model systems are carbon-oxygen (C/O) white dwarfs, while the secondaries are made either of C/O or of pure helium. We account for tidal effects in the initial conditions in a self-consistent way, and consider initially well-separated systems with slow inspiral rates. We study the merger evolution using an adaptive mesh refinement, reactive, Eulerian code in three dimensions, assuming symmetry across the orbital plane. We use a co-rotating reference frame to minimize the effects of numerical diffusion, and solve for self-gravity using a multi-grid approach. We find a novel detonation mechanism in C/O mergers with massive primaries. Here the detonation occurs in the primary's core and relies on the combined action of tidal heating, accretion heating, and self-heating due to nuclear burning. The exploding structure is compositionally stratified, with a reverse shock formed at the surface of the dense ejecta. The existence of such a shock has not been reported elsewhere. The explosion energy ($1.6\times 10^{51}$ erg) and $^{56}$Ni mass (0.86 M$_\odot$) are consistent with a SN Ia at the bright end of the luminosity distribution, with an approximated decline rate of $\Delta m_{15}(B)\approx 0.99$. Our study does not support double-detonation scenarios in the case of a system with a 0.6 M$_\odot$ helium secondary and a 0.9 M$_\odot$ primary. Although the accreted helium detonates, it fails to ignite carbon at the base of the boundary layer or in the primary's core.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

We present simulations of the magnetized interstellar medium (ISM) in models of massive star forming (40 Msun / yr) disk galaxies with high gas surface densities (~100 Msun / pc^2) similar to observed star forming high-redshift disks. We assume that type II supernovae deposit 10 per cent of their energy into the ISM as cosmic rays and neglect the additional deposition of thermal energy or momentum. With a typical Galactic diffusion coefficient for CRs (3e28 cm^2 / s) we demonstrate that this process alone can trigger the local formation of a strong low density galactic wind maintaining vertically open field lines. Driven by the additional pressure gradient of the relativistic fluid the wind speed can exceed 1000 km/s, much higher than the escape velocity of the galaxy. The global mass loading, i.e. the ratio of the gas mass leaving the galactic disk in a wind to the star formation rate becomes of order unity once the system has settled into an equilibrium. We conclude that relativistic particles accelerated in supernova remnants alone provide a natural and efficient mechanism to trigger winds similar to observed mass-loaded galactic winds in high-redshift galaxies. These winds also help explaining the low efficiencies for the conversion of gas into stars in galaxies as well as the early enrichment of the intergalactic medium with metals. This mechanism can be at least of similar importance than the traditionally considered momentum feedback from massive stars and thermal and kinetic feedback from supernova explosions.

We investigate streaming instability using two-fluid approximation (neutral gas and dust) in a quasi-global, unstratified protoplanetary disc, with the help of PIERNIK code. We compare amplification rate of the eigen-mode in numerical simulations, with the corresponding growth resulting from the linear stability analysis of full system of Euler's equation including aerodynamic drag. Following Youdin & Goodman (2005) we show that (1) rapid dust clumping occurs due to the difference in azimuthal velocities of gas and dust, coupled by the drag force, (2) initial density perturbations are amplified by several orders of magnitude. We demonstrate that the multi-fluid extension of the simple and efficient Relaxing TVD scheme, implemented in PIERNIK, leads to results, which are compatible with those obtained with other methods.

(abridged) We study the hydrodynamic evolution of a non-spherical core-collapse supernova in two spatial dimensions. We find that our model displays a strong tendency to expand toward the pole. We demonstrate that this expansion is a physical property of the low-mode, SASI instability. The SASI leaves behind a large lateral velocity gradient in the post shock layer which affects the evolution for minutes and hours later. This results in a prolate deformation of the ejecta and a fast advection of Ni-rich material from moderate latitudes to the polar regions. This effect might actually be responsible for the global asymmetry of the nickel lines in SN 1987A. The simulations demonstrate that significant radial and lateral motions in the post-shock region, produced by convective overturn and the SASI during the early explosion phase, contribute to the evolution for minutes and hours after shock revival. They lead to both later clump formation, and a significant prolate deformation of the ejecta which are observed even as late as one week after the explosion. As pointed out recently by Kjaer et al., such an ejecta morphology is in good agreement with the observational data of SN 1987A. Systematic future studies are needed to investigate how the SASI-induced late-time lateral expansion depends on the dominant mode of the SASI, and to which extent it is affected by the dimensionality of the simulations. The impact on and importance of the SASI for the distribution of iron group nuclei and the morphology of the young SNR argues for future three-dimensional explosion and post-explosion studies on singularity-free grids that cover the entire sphere. Given the results of our 2D resolution study, present 3D simulations must be regarded as underresolved, and their conclusions must be verified by a proper numerical convergence analysis in three dimensions.

We investigate numerically the orbital evolution of massive extrasolar planets within central cavities of their parent protoplanetary discs. Assuming that they arrive at the inner edge of the disc due to type II migration, we show that they spiral further in. We find that in magnetospheric cavities more massive planets stop migrating at a larger distance from the edge of the disc. This effect may qualitatively explain the correlation between masses and orbital periods found for massive planets with P shorter than 5 days.

We review recent results of SPH simulations of gravitational instability in gaseous protoplanetary disks,emphasizing the role of thermodynamics in both isolated and binary systems. Contradictory results appeared in the literature regarding disk fragmentation at tens of AU from the central star are likely due to the different treatment of radiation physics as well as reflecting different initial conditions. Further progress on the subject requires extensive comparisons between different codes with the requirement that the same initial conditions are adopted. It is discussed how the local conditions of the disks undergoing fragmentation at $R < 25$ AU in recent SPH simulations are in rough agreement with the prediction of analytical models, with small differences being likely related to the inability of analytical models to account for the dynamics and thermodynamics of three-dimensional spiral shocks. We report that radically different adaptive hydrodynamical codes, SPH and adaptive mesh refinement (AMR), yield very similar results on disk fragmentation at comparable resolution in the simple case of an isothermal equation of state. A high number of refinements in AMR codes is necessary but not sufficient to correctly follow fragmentation, rather an initial resolution of the grid high enough to capture the wavelength of the strongest spiral modes when they are still barely nonlinear is essential. These tests represent a useful benchmark and a starting point for a forthcoming code comparison with realistic radiation physics.

We investigate both analytically and numerically the motion of massless particles orbiting primary star in a close circular binary system with particular focus on the gas drag effects. These are the first calculations with particles ranging in size from 1m to 10km, which account for the presence of a tidally perturbed gaseous disk. We have found numerically that the radial mass transport by the tidal waves plays a crucial role in the orbital evolution of particles. Numerical results are confirmed analytical calculations that do not assume anything about origin of the radial flow in the disk. We demonstrate that the migration rate of a particle in a disk out of radial equilibrium is enhanced due to the enhanced mass flux of gas colliding with the particle and the migration is always directed inward regardless of the sign of the radial gas velocity. Within the framework of the perturbation theory we derive general, formulae for short-term variations of the particle semi-major axis, eccentricity and inclination in such disk. The formulae account for departures from axial symmetry by introducing effective components of the gas velocity. They agree with numerical results within several percent. We have also found in numerical simulations that the tidal waves introduce coherence in periastron longitude and eccentricity for particles on neighbouring orbits. The degree of the coherence depends on the particle size and on the distance from the primary star, being most prominent for particles with 10m radius. The results are important mainly in the context of planetesimal formation and, to a lesser degree, during the early planetesimal accretion stage.

We have carried out a hydrodynamical code comparison study of interacting multiphase fluids. The two commonly used techniques of grid and smoothed particle hydrodynamics (SPH) show striking differences in their ability to model processes that are fundamentally important across many areas of astrophysics. Whilst Eulerian grid based methods are able to resolve and treat important dynamical instabilities, such as Kelvin-Helmholtz or Rayleigh-Taylor, these processes are poorly or not at all resolved by existing SPH techniques. We show that the reason for this is that SPH, at least in its standard implementation, introduces spurious pressure forces on particles in regions where there are steep density gradients. This results in a boundary gap of the size of the SPH smoothing kernel over which information is not transferred.

We perform numerical simulations of a disc-planet system using various grid-based and smoothed particle hydrodynamics (SPH) codes. The tests are run for a simple setup where Jupiter and Neptune mass planets on a circular orbit open a gap in a protoplanetary disc during a few hundred orbital periods. We compare the surface density contours, potential vorticity and smoothed radial profiles at several times. The disc mass and gravitational torque time evolution are analyzed with high temporal resolution. There is overall consistency between the codes. The density profiles agree within about 5% for the Eulerian simulations while the SPH results predict the correct shape of the gap although have less resolution in the low density regions and weaker planetary wakes. The disc masses after 200 orbital periods agree within 10%. The spread is larger in the tidal torques acting on the planet which agree within a factor 2 at the end of the simulation. In the Neptune case the dispersion in the torques is greater than for Jupiter, possibly owing to the contribution from the not completely cleared region close to the planet.

We study axisymmetric models of layered protoplanetary discs taking radiative transfer effects into account, and allowing for a residual viscosity in the dead zone. We also explore the effect of different viscosity prescriptions. In addition to the ring instability reported in the first paper of the series we find an oscillatory instability of the dead zone, accompanied by variations of the accretion rate onto the central star. We provide a simplified analytical description explaining the mechanism of the oscillations. Finally, we find that the residual viscosity enables stationary accretion in large regions of layered discs. Based on results obtained with the help of a simple 1-D hydrocode we identify these regions, and discuss conditions in which layered discs can give rise to FU~Orionis phenomena.

We present three-dimensional models of accretion disks in U Gem - like systems and calculate their Doppler tomograms. The tomograms are based on two different assumptions concerning the origin of line emission from the disk. The assumption of lines originating due to irradiation of the surface layer of the disk by the central source leads to a better agreement with observations. We argue that fully three-dimensional modelling is necessary to properly interpret the observed tomograms.

The slow and dense wind from a symbiotic red giant can be significantly deflected toward the orbital plane by the gravitational pull of the companion star. In such an environment, the ionizing radiation from the companion creates a highly asymmetric HII region. We present three-dimensional models of HII regions in symbiotic S-type systems, for which we calculate radio maps and radio spectra. We show that the standard assumption of spherically symmetric RG wind results in wrong shapes, sizes and spectra of ionized regions, which in turn affects the observational estimates of orbital separation and mass loss rate. A sample of radio maps and radio spectra of our models is presented and the results are discussed in relation to observational data.