Magnetic fields are ubiquitous in the Universe. The Sun's magnetic field drives the solar wind and causes solar flares and other energetic surface phenomena that profoundly affect space weather here on Earth. The first magnetic field in a star other than the Sun was detected in 1947 in the peculiar A-type star 78 Vir. It is now known that the magnetic fields of the Sun and other low-mass stars (<1.5 solar masses) are generated in-situ by a dynamo process in their turbulent, convective envelopes. Unlike such stars, intermediate-mass and high-mass stars (>1.5 solar masses; referred to as "massive" stars here) have relatively quiet, radiative envelopes where a solar-like dynamo cannot operate. However, about 10% of them, including 78 Vir, have strong, large-scale surface magnetic fields whose origin has remained a major mystery. The massive star $\tau$ Sco is a prominent member of this group and appears to be surprisingly young compared to other presumably coeval members of the Upper Scorpius association. Here, we present the first 3D magneto-hydrodynamical simulations of the coalescence of two massive main-sequence stars and 1D stellar evolution computations of the subsequent evolution of the merger product that can explain $\tau$ Sco's magnetic field, apparent youth and other observed characteristics. We argue that field amplification in stellar mergers is a general mechanism to form strongly-magnetised massive stars. These stars are promising progenitors of those neutron stars that host the strongest magnetic fields in the Universe, so-called magnetars, and that may give rise to some of the enigmatic fast radio bursts. Strong magnetic fields affect the explosions of core-collapse supernovae and, moreover, those magnetic stars that have rapidly-rotating cores at the end of their lives might provide the right conditions to power long-duration gamma-ray bursts and super-luminous supernovae.

In spring 2015, the Leibniz Supercomputing Centre (Leibniz-Rechenzentrum, LRZ), installed their new Peta-Scale System SuperMUC Phase2. Selected users were invited for a 28 day extreme scale-out block operation during which they were allowed to use the full system for their applications. The following projects participated in the extreme scale-out workshop: BQCD (Quantum Physics), SeisSol (Geophysics, Seismics), GPI-2/GASPI (Toolkit for HPC), Seven-League Hydro (Astrophysics), ILBDC (Lattice Boltzmann CFD), Iphigenie (Molecular Dynamic), FLASH (Astrophysics), GADGET (Cosmological Dynamics), PSC (Plasma Physics), waLBerla (Lattice Boltzmann CFD), Musubi (Lattice Boltzmann CFD), Vertex3D (Stellar Astrophysics), CIAO (Combustion CFD), and LS1-Mardyn (Material Science). The projects were allowed to use the machine exclusively during the 28 day period, which corresponds to a total of 63.4 million core-hours, of which 43.8 million core-hours were used by the applications, resulting in a utilization of 69%. The top 3 users were using 15.2, 6.4, and 4.7 million core-hours, respectively.

We present results based on observations of SN 2015H which belongs to the small group of objects similar to SN 2002cx, otherwise known as type Iax supernovae. The availability of deep pre-explosion imaging allowed us to place tight constraints on the explosion epoch. Our observational campaign began approximately one day post-explosion, and extended over a period of about 150 days post maximum light, making it one of the best observed objects of this class to date. We find a peak magnitude of M$_r$ = -17.27 $\pm$ 0.07, and a ($\Delta m_{15})_r$ = 0.69 $\pm$ 0.04. Comparing our observations to synthetic spectra generated from simulations of deflagrations of Chandrasekhar mass carbon-oxygen white dwarfs, we find reasonable agreement with models of weak deflagrations that result in the ejection of ~0.2 M$_{\odot}$ of material containing ~0.07 M$_{\odot}$ of 56Ni. The model light curve however, evolves more rapidly than observations, suggesting that a higher ejecta mass is to be favoured. Nevertheless, empirical modelling of the pseudo-bolometric light curve suggests that $\lesssim$0.6 M_sun of material was ejected, implying that the white dwarf is not completely disrupted, and that a bound remnant is a likely outcome.

There has been much debate about the origin of the diffuse $\gamma$--ray background in the MeV range. At lower energies, AGNs and Seyfert galaxies can explain the background, but not above $\simeq$0.3 MeV. Beyond $\sim$10 MeV blazars appear to account for the flux observed. That leaves an unexplained gap for which different candidates have been proposed, including annihilations of WIMPS. One candidate are Type Ia supernovae (SNe Ia). Early studies concluded that they were able to account for the $\gamma$--ray background in the gap, while later work attributed a significantly lower contribution to them. All those estimates were based on SN Ia explosion models which did not reflect the full 3D hydrodynamics of SNe Ia explosions. In addition, new measurements obtained since 2010 have provided new, direct estimates of high-z SNe Ia rates beyond $z\sim$2. We take into account these new advances to see the predicted contribution to the gamma--ray background. We use here a wide variety of explosion models and a plethora of new measurements of SNe Ia rates. SNe Ia still fall short of the observed background. Only for a fit, which would imply $\sim$150\% systematic error in detecting SNe Ia events, do the theoretical predictions approach the observed fluxes. This fit is, however, at odds at the highest redshifts with recent SN Ia rates estimates. Other astrophysical sources such as FSRQs do match the observed flux levels in the MeV regime, while SNe Ia make up to 30--50\% of the observed flux.

The measurement of gamma-ray lines from the decay chain of 56Ni provides unique information about the explosion in supernovae. The 56Ni freshly-produced in the supernova powers the optical light curve, as it emits gamma-rays upon its radioactive decay to 56Co and then 56Fe. Gamma-ray lines from 56Co decay are expected to become directly visible through the overlying white dwarf material several weeks after the explosion, as they progressively penetrate the overlying material of the supernova envelope, diluted as it expands. The lines are expected to be Doppler-shifted or broadened from the kinematics of the 56Ni ejecta. With the SPI spectrometer on INTEGRAL and using an improved instrumental background method, we detect the two main lines from 56Co decay at 847 and 1238 keV from SN2014J at 3.3 Mpc, significantly Doppler-broadened, and at intensities (3.65+/-1.21) 10^-4 and (2.27+/-0.69) 10^-4 ph cm^-2s^-1, respectively, at brightness maximum. We measure their rise towards a maximum after about 60-100 days and decline thereafter. The intensity ratio of the two lines is found consistent with expectations from 56Co decay (0.62+/-0.28 at brightness maximum, expected is 0.68). We find that the broad lines seen in the late, gamma-ray transparent phase are not representative for the early gamma-ray emission, and rather notice the emission spectrum to be complex and irregular until the supernova is fully transparent to gamma-rays, with progressive uncovering of the bulk of 56Ni. We infer that the explosion morphology is not spherically symmetric, both in the distribution of 56Ni and of the unburnt material which occults the 56Co emission. Comparing light curves from different plausible models, the resulting 56Ni mass is determined as 0.49+/-0.09 M_sol.

We perform a model independent reconstruction of the cosmic expansion rate based on type Ia supernova data. Using the Union 2.1 data set, we show that the Hubble parameter behaviour allowed by the data without making any hypothesis about cosmological model or underlying gravity theory is consistent with a flat LCDM universe having H_0 = 70.43 +- 0.33 and Omega_m=0.297 +- 0.020, weakly dependent on the choice of initial scatter matrix. This is in closer agreement with the recently released Planck results (H_0 = 67.3 +- 1.2, Omega_m = 0.314 +- 0.020) than other standard analyses based on type Ia supernova data. We argue this might be an indication that, in order to tackle subtle deviations from the standard cosmological model present in type Ia supernova data, it is mandatory to go beyond parametrized approaches.

In a companion paper, Seitenzahl et al. (2013) presented a set of three-dimensional delayed detonation models for thermonuclear explosions of near-Chandrasekhar mass white dwarfs (WDs). Here, we present multi-dimensional radiative transfer simulations that provide synthetic light curves and spectra for those models. The model sequence explores both changes in the strength of the deflagration phase (controlled by the ignition configuration) and the WD central density. In agreement with previous studies, we find that the strength of the deflagration significantly affects the explosion and the observables. Variations in the central density also have an influence on both brightness and colour, but overall it is a secondary parameter in our set of models. In many respects, the models yield a good match to normal Type Ia supernovae (SNe Ia): peak brightness, rise/decline time scales and synthetic spectra are all in reasonable agreement. There are, however, several differences. In particular, the models are too red around maximum light, manifest spectral line velocities that are a little too high and yield I-band light curves that do not match observations. Although some of these discrepancies may simply relate to approximations made in the modelling, some pose real challenges to the models. If viewed as a complete sequence, our models do not reproduce the observed light-curve width-luminosity relation (WLR) of SNe Ia: all our models show similar B-band decline rates, irrespective of peak brightness. This suggests that simple variations in the strength of the deflagration phase in Chandrasekhar-mass deflagration-to-detonation models do not readily explain the observed diversity of normal SNe Ia. This may imply that some other parameter within the Chandrasekhar-mass paradigm is key to the WLR, or that a substantial fraction of normal SNe Ia arise from an alternative explosion scenario.

Motivated by the fact that calibrated light curves of Type Ia supernovae (SNe Ia) have become a major tool to determine the expansion history of the Universe, considerable attention has been given to, both, observations and models of these events over the past 15 years. Here, we summarize new observational constraints, address recent progress in modeling Type Ia supernovae by means of three-dimensional hydrodynamic simulations, and discuss several of the still open questions. It will be be shown that the new models have considerable predictive power which allows us to study observable properties such as light curves and spectra without adjustable non-physical parameters. This is a necessary requisite to improve our understanding of the explosion mechanism and to settle the question of the applicability of SNe Ia as distance indicators for cosmology. We explore the capabilities of the models by comparing them with observations and we show how such models can be applied to study the origin of the diversity of SNe Ia.

Nucleosynthetic yield predictions for multi-dimensional simulations of thermonuclear supernovae generally rely on the tracer particle method to obtain isotopic information of the ejected material for a given supernova simulation. We investigate how many tracer particles are required to determine converged integrated total nucleosynthetic yields. For this purpose, we conduct a resolution study in the number of tracer particles for different hydrodynamical explosion models at fixed spatial resolution. We perform hydrodynamic simulations on a co-expanding Eulerian grid in two dimensions assuming rotational symmetry for both pure deflagration and delayed detonation Type Ia supernova explosions. Within a given explosion model, we vary the number of tracer particles to determine the minimum needed for the method to give a robust prediction of the integrated yields of the most abundant nuclides. For the first time, we relax the usual assumption of constant tracer particle mass and introduce a radially vary- ing distribution of tracer particle masses. We find that the nucleosynthetic yields of the most abundant species (mass fraction > 10E-5) are reasonably well predicted for a tracer number as small as 32 per axis and direction - more or less independent of the explosion model. We conclude that the number of tracer particles that were used in extant published works appear to have been sufficient as far as integrated yields are concerned for the most copiously produced nuclides. Additionally we find that a suitably chosen tracer mass distribution can improve convergence for nuclei produced in the outer layer of the supernova where the constant tracer mass prescription suffers from poor spatial resolution.

Motivated by recent models involving off-centre ignition of type Ia supernova explosions, we undertake three-dimensional time-dependent radiation transport simulations to investigate the range of bolometric light curve properties that could be observed from supernovae in which there is a lop-sided distribution of the products from nuclear burning. We consider both a grid of artificial toy models which illustrate the conceivable range of effects and a recent three-dimensional hydrodynamical explosion model. We find that observationally significant viewing angle effects are likely to arise in such supernovae and that these may have important ramifications for the interpretation of the observed diversity of type Ia supernova and the systematic uncertainties which relate to their use as standard candles in contemporary cosmology.

The recent discovery of a very bright type Ia supernova, SNLS-03D3bb ($\equiv$ SN 2003fg), in the Supernova Legacy Survey (SNLS) has raised the question of whether super-Chandrasekhar-mass white dwarf stars are needed to explain such explosions. In principle, such a progenitor could form by the mergers of two rather massive white dwarfs. Binary systems of two white dwarfs in close orbit where their total mass significantly exceeds the Chandrasekhar mass, have not yet been found. Therefore SNLS-03D3bb could establish the first clear case of a double-degenerate progenitor of a (peculiar) type Ia supernovae. Moreover, if this interpretation is correct, it casts some doubt on the universality of the calibration relations used to make SNe Ia distance indicators for cosmology. Here we present a critical discussion of the expected observational finger prints of super-Chandrasekhar-mass explosions; in important respects, these are not consistent with the observations of SNLS-03D3bb. We demonstrate that the lop-sided explosion of a Chandrasekhar-mass white dwarf could provide a better explanation.