This article describes the operation of the near-infrared wavefront sensing based Adaptive Optics (AO) system CIAO. The Coud\'e Infrared Adaptive Optics (CIAO) system is a central auxiliary component of the Very Large Telescope (VLT) interferometer (VLTI). It enables in particular the observations of the Galactic Center (GC) using the GRAVITY instrument. GRAVITY is a highly specialized beam combiner, a device that coherently combines the light of the four 8-m telescopes and finally records interferometric measurements in the K-band on 6 baselines simultaneously. CIAO compensates for phase disturbances caused by atmospheric turbulence, which all four 8 m Unit Telescopes (UT) experience during observation. Each of the four CIAO units generates an almost diffraction-limited image quality at its UT, which ensures that maximum flux of the observed stellar object enters the fibers of the GRAVITY beam combiner. We present CIAO performance data obtained in the first 3 years of operation as a function of weather conditions. We describe how CIAO is configured, calibrated, and finally used for observations with GRAVITY. In addition, we focus on the outstanding features of the near-infrared sensitive Saphira detector, which is used for the first time on Paranal, and show how it works as a wavefront sensor detector.

Ultra-hot giant exoplanets receive thousands of times Earth's insolation. Their high-temperature atmospheres (>2,000 K) are ideal laboratories for studying extreme planetary climates and chemistry. Daysides are predicted to be cloud-free, dominated by atomic species and substantially hotter than nightsides. Atoms are expected to recombine into molecules over the nightside, resulting in different day-night chemistry. While metallic elements and a large temperature contrast have been observed, no chemical gradient has been measured across the surface of such an exoplanet. Different atmospheric chemistry between the day-to-night ("evening") and night-to-day ("morning") terminators could, however, be revealed as an asymmetric absorption signature during transit. Here, we report the detection of an asymmetric atmospheric signature in the ultra-hot exoplanet WASP-76b. We spectrally and temporally resolve this signature thanks to the combination of high-dispersion spectroscopy with a large photon-collecting area. The absorption signal, attributed to neutral iron, is blueshifted by -11+/-0.7 km s-1 on the trailing limb, which can be explained by a combination of planetary rotation and wind blowing from the hot dayside. In contrast, no signal arises from the nightside close to the morning terminator, showing that atomic iron is not absorbing starlight there. Iron must thus condense during its journey across the nightside.

We present predictions for the gravitational-wave (GW) emission of three-dimensional supernova (SN) simulations performed for a 15 solar-mass progenitor with the Prometheus-Vertex code using energy-dependent, three-flavor neutrino transport. The progenitor adopted from stellar evolution calculations including magnetic fields had a fairly low specific angular momentum (j_Fe <~ 10^{15} cm^2/s) in the iron core (central angular velocity ~0.2 rad/s), which we compared to simulations without rotation and with artificially enhanced rotation (j_Fe <~ 2*10^{16} cm^2/s; central angular velocity ~0.5 rad/s). Our results confirm that the time-domain GW signals of SNe are stochastic, but possess deterministic components with characteristic patterns at low frequencies (<~200 Hz), caused by mass motions due to the standing accretion shock instability (SASI), and at high frequencies, associated with gravity-mode oscillations in the surface layer of the proto-neutron star (PNS). Non-radial mass motions in the post-shock layer as well as PNS convection are important triggers of GW emission, whose amplitude scales with the power of the hydrodynamic flows. There is no monotonic increase of the GW amplitude with rotation, but a clear correlation with the strength of SASI activity. Our slowly rotating model is a fainter GW emitter than the non-rotating model because of weaker SASI activity and damped convection in the post-shock layer and PNS. In contrast, the faster rotating model exhibits a powerful SASI spiral mode during its transition to explosion, producing the highest GW amplitudes with a distinctive drift of the low-frequency emission peak from ~80-100 Hz to ~40-50 Hz. This migration signifies shock expansion, whereas non-exploding models are discriminated by the opposite trend.

Cassiopeia A, one of the most intriguing galactic supernova remnants, has been a target of many observational efforts including most recent observations by ALMA, Hubble, Herschel, Spitzer, NuSTAR, Integral, and other observatories. We use recent gamma-ray lines observations of the radioactive products of Cas A supernova explosive nucleosynthesis as well as spectral energy densities derived for Cas A at infrared wavelengths to speculate about the possibility of radioactive 44Ti being locked into large dust grains. This suggestion is also supported by the possible observation of a pre-supernova outburst about 80 years before the actual Cas A supernova explosion in 1671 AD by Italian astronomer G.D. Cassini. The plausibility of such a scenario is discussed also with reference to recent supernovae, and to the contribution of core-collapse supernovae to the overall dust production in the Galaxy.

This work presents an interferometric study of the massive-binary fraction in the Orion Trapezium Cluster with the recently comissioned GRAVITY instrument. We observe a total of 16 stars of mainly OB spectral type. We find three previously unknown companions for $\theta ^1$ Ori B, $\theta ^2$ Ori B, and $\theta ^2$ Ori C. We determine a separation for the previously suspected companion of NU Ori. We confirm four companions for $\theta ^1$ Ori A, $\theta ^1$ Ori C, $\theta ^1$ Ori D, and $\theta ^2$ Ori A, all with substantially improved astrometry and photometric mass estimates. We refine the orbit of the eccentric high-mass binary $\theta ^1$ Ori C and we are able to derive a new orbit for $\theta ^1$ Ori D. We find a system mass of 21.7 $M_{\odot}$ and a period of $53$ days. Together with other previously detected companions seen in spectroscopy or direct imaging, eleven of the 16 high-mass stars are multiple systems. We obtain a total number of 22 companions with separations up to 600 AU. The companion fraction of the early B and O stars in our sample is about 2, significantly higher than in earlier studies of mostly OB associations. The separation distribution hints towards a bimodality. Such a bimodality has been previously found in A stars, but rarely in OB binaries, which up to this point have been assumed to be mostly compact with a tail of wider companions. We also do not find a substantial population of equal-mass binaries. The observed distribution of mass ratios declines steeply with mass, and like the direct star counts, indicates that our companions follow a standard power law initial mass function. Again, this is in contrast to earlier findings of flat mass ratio distributions in OB associations. We exclude collision as a dominant formation mechanism but find no clear preference for core accretion or competitive accretion.

Eta Car is one of the most intriguing luminous blue variables in the Galaxy. Observations and models at different wavelengths suggest a central binary with a 5.54 yr period residing in its core. 2D and 3D radiative transfer and hydrodynamic simulations predict a primary with a dense and slow stellar wind that interacts with the faster and lower density wind of the secondary. The wind-wind collision scenario suggests that the secondary's wind penetrates the primary's wind creating a low-density cavity in it, with dense walls where the two winds interact. We aim to trace the inner ~5-50 au structure of Eta Car's wind-wind interaction, as seen through BrG and, for the first time, through the He I 2s-2p line. We have used spectro-interferometric observations with GRAVITY at the VLTI. Our modeling of the continuum allows us to estimate its FWHM angular size close to 2 mas and an elongation ratio of 1.06 +/- 0.05 over a PA = 130 +/- 20 deg. Our CMFGEN modeling helped us to confirm that the role of the secondary should be taken into account to properly reproduce the observed BrG and He I lines. Chromatic images across BrG reveal a southeast arc-like feature, possibly associated to the hot post-shocked winds flowing along the cavity wall. The images of He I 2s-2p served to constrain the 20 mas structure of the line-emitting region. The observed morphology of He I suggests that the secondary is responsible for the ionized material that produces the line profile. Both the BrG and the He I 2s-2p maps are consistent with previous hydrodynamical models of the colliding wind scenario. Future dedicated simulations together with an extensive interferometric campaign are necessary to refine our constraints on the wind and stellar parameters of the binary, which finally will help us predict the evolutionary path of Eta Car.

The highly elliptical, 16-year-period orbit of the star S2 around the massive black hole candidate Sgr A* is a sensitive probe of the gravitational field in the Galactic centre. Near pericentre at 120 AU, ~1400 Schwarzschild radii, the star has an orbital speed of ~7650 km/s, such that the first-order effects of Special and General Relativity have now become detectable with current capabilities. Over the past 26 years, we have monitored the radial velocity and motion on the sky of S2, mainly with the SINFONI and NACO adaptive optics instruments on the ESO Very Large Telescope, and since 2016 and leading up to the pericentre approach in May 2018, with the four-telescope interferometric beam-combiner instrument GRAVITY. From data up to and including pericentre, we robustly detect the combined gravitational redshift and relativistic transverse Doppler effect for S2 of z ~ 200 km/s / c with different statistical analysis methods. When parameterising the post-Newtonian contribution from these effects by a factor f, with f = 0 and f = 1 corresponding to the Newtonian and general relativistic limits, respectively, we find from posterior fitting with different weighting schemes f = 0.90 +/- 0.09 (stat) +\- 0.15 (sys). The S2 data are inconsistent with pure Newtonian dynamics.

The birth properties of neutron stars yield important information on the still debated physical processes that trigger the explosion and on intrinsic neutron-star physics. These properties include the high space velocities of young neutron stars with average values of several 100 km/s, whose underlying "kick" mechanism is not finally clarified. There are two competing possibilities that could accelerate NSs during their birth: anisotropic ejection of either stellar debris or neutrinos. We here present new evidence from X-ray measurements that chemical elements between silicon and calcium in six young gaseous supernova remnants are preferentially expelled opposite to the direction of neutron star motion. There is no correlation between the kick velocities and magnetic field strengths of these neutron stars. Our results support a hydrodynamic origin of neutron-star kicks connected to asymmetric explosive mass ejection, and they conflict with neutron-star acceleration scenarios that invoke anisotropic neutrino emission caused by particle and nuclear physics in combination with very strong neutron-star magnetic fields.

Geometrically thick tori with constant specific angular momentum have been widely used in the last decades to construct numerical models of accretion flows onto black holes. Such discs are prone to a global non-axisymmetric hydrodynamic instability, known as Papaloizou-Pringle instability (PPI), which can redistribute angular momentum and also lead to an emission of gravitational waves. It is, however, not clear yet how the development of the PPI is affected by the presence of a magnetic field and by the concurrent development of the magnetorotational instability (MRI). We present a numerical analysis using three-dimensional GRMHD simulations of the interplay between the PPI and the MRI considering, for the first time, an analytical magnetized equilibrium solution as initial condition. In the purely hydrodynamic case, the PPI selects as expected the large-scale $m=1$ azimuthal mode as the fastest growing and non-linearly dominant mode. However, when the torus is threaded by a weak toroidal magnetic field, the development of the MRI leads to the suppression of large-scale modes and redistributes power across smaller scales. If the system starts with a significantly excited $m=1$ mode, the PPI can be dominant in a transient phase, before being ultimately quenched by the MRI. Such dynamics may well be important in compact star mergers and tidal disruption events.

A large fraction of core-collapse supernovae (CCSNe), 30-50%, are expected to originate from the low-mass end of progenitors with $M_{\rm ZAMS}~= 8-12~M_\odot$. However, degeneracy effects make stellar evolution modelling of such stars challenging, and few predictions for their supernova light curves and spectra have been presented. Here we calculate synthetic nebular spectra of a 9 $M_\odot$ Fe CCSN model exploded with the neutrino mechanism. The model predicts emission lines with FWHM$\sim$1000 km/s, including signatures from each deep layer in the metal core. We compare this model to observations of the three subluminous IIP SNe with published nebular spectra; SN 1997D, SN 2005cs, and SN 2008bk. The prediction of both line profiles and luminosities are in good agreement with SN 1997D and SN 2008bk. The close fit of a model with no tuning parameters provides strong evidence for an association of these objects with low-mass Fe CCSNe. For SN 2005cs, the interpretation is less clear, as the observational coverage ended before key diagnostic lines from the core had emerged. We perform a parameterised study of the amount of explosively made stable nickel, and find that none of these three SNe show the high $^{58}$Ni/$^{56}$Ni ratio predicted by current models of electron capture SNe (ECSNe) and ECSN-like explosions. Combined with clear detection of lines from O and He shell material, these SNe rather originate from Fe core progenitors. We argue that the outcome of self-consistent explosion simulations of low-mass stars, which gives fits to many key observables, strongly suggests that the class of subluminous Type IIP SNe is the observational counterpart of the lowest mass CCSNe.

We discuss torsional oscillations of highly magnetised neutron stars (magnetars) using two-dimensional, magneto-elastic-hydrodynamical simulations. Our model is able to explain both the low- and high-frequency quasi-periodic oscillations (QPOs) observed in magnetars. The analysis of these oscillations provides constraints on the breakout magnetic-field strength, on the fundamental QPO frequency, and on the frequency of a particularly excited overtone. More importantly, we show how to use this information to generically constraint properties of high-density matter in neutron stars, employing Bayesian analysis. In spite of current uncertainties and computational approximations, our model-dependent Bayesian posterior estimates for SGR 1806-20 yield a magnetic-field strength $\bar B\sim 2.1^{+1.3}_{-1.0}\times10^{15}\,$G and a crust thickness of $\Delta r = 1.6^{+0.7}_{-0.6}$ km, which are both in remarkable agreement with observational and theoretical expectations, respectively (1-$\sigma$ error bars are indicated). Our posteriors also favour the presence of a superfluid phase in the core, a relatively low stellar compactness, $M/R<0.19$, indicating a relatively stiff equation of state and/or low mass neutron star, and high shear speeds at the base of the crust, $c_s>1.4\times10^8\,$cm/s. Although the procedure laid out here still has large uncertainties, these constraints could become tighter when additional observations become available.

Magneto-elastic oscillations of neutron stars are believed to explain observed quasi-periodic oscillations (QPOs) in the decaying tail of the giant flares of highly magnetized neutron stars (magnetars). Strong efforts of the theoretical modelling from different groups have increased our understanding of this phenomenon significantly. Here, we discuss some constraints on the matter in neutron stars that arise if the interpretation of the observations in terms of superfluid, magneto-elastic oscillations is correct. To explain the observed modulation of the light curve of the giant flare, we describe a model that allows the QPOs to couple to the stellar exterior through the magnetic field. In this magnetosphere, the shaking magnetic field induces currents that provide scattering targets for resonant cyclotron scattering of photons, which is calculated with a Monte-Carlo approach and coupled to a code that calculates the momentum distribution of the charge carriers as a one-dimensional accelerator problem. We show first results of a simplified, but self-consistent momentum distribution, i.e. a waterbag distribution, and of the corresponding spectra.

To investigate the inner regions of protoplanetary disks, we performed near-infrared interferometric observations of the classical TTauri binary system S CrA. We present the first VLTI-GRAVITY high spectral resolution ($R\sim$4000) observations of a classical TTauri binary, S CrA (composed of S CrA N and S CrA S and separated by $\sim$1.4"), combining the four 8-m telescopes in dual-field mode. Our observations in the near-infrared K-band continuum reveal a disk around each binary component, with similar half-flux radii of about 0.1 au at d$\sim$130 pc, inclinations ($i=$28$\pm$3$^o$\ and $i=$22$\pm$6$^o$), and position angles (PA=0$^o\pm$6$^o$ and PA=-2$^o\pm$12$^o$), suggesting that they formed from the fragmentation of a common disk. The S CrA N spectrum shows bright HeI and Br$\gamma$ line emission exhibiting inverse P-Cygni profiles, typically associated with infalling gas. The continuum-compensated Br$\gamma$ line visibilities of S CrA N show the presence of a compact Br$\gamma$ emitting region the radius of which is about $\sim$0.06 au, which is twice as big as the truncation radius. This component is mostly tracing a wind. Moreover, a slight radius change between the blue- and red-shifted Br$\gamma$ line components is marginally detected. The presence of an inverse P-Cygni profile in the HeI and Br$\gamma$ lines, along with the tentative detection of a slightly larger size of the blue-shifted Br$\gamma$ line component, hint at the simultaneous presence of a wind and magnetospheric accretion in S CrA N.

We discuss the relevance of the magnetorotational instability (MRI) in core-collapse supernovae (CCSNe). Our recent numerical studies show that in CCSNe, the MRI is terminated by parasitic instabilities of the Kelvin-Helmholtz type. To determine whether the MRI can amplify initially weak magnetic fields to dynamically relevant strengths in CCSNe, we performed three-dimensional simulations of a region close to the surface of a differentially rotating proto-neutron star in non-ideal magnetohydrodynamics with two different numerical codes. We find that under the conditions prevailing in proto-neutron stars, the MRI can amplify the magnetic field by (only) one order of magnitude. This severely limits the role of MRI channel modes as an agent amplifying the magnetic field in proto-neutron stars starting from small seed fields.

We observe the HMXB BP Cru using interferometry in the near-infrared K band with VLTI/GRAVITY. Continuum visibilities are at most partially resolved, consistent with the predicted size of the hypergiant. Differential visibility amplitude ($\Delta |V| \sim 5\%$) and phase ($\Delta \phi \sim 2 \deg$) signatures are observed across the HeI $2.059 \mu$m and Br$\gamma$ lines, the latter seen strongly in emission, unusual for the donor star's spectral type. For a baseline $B \sim 100$m, the differential phase RMS $\sim 0.2 \deg$ corresponds to an astrometric precision of $\sim 2 \mu$as. A model-independent analysis in the marginally resolved limit of interferometry reveals asymmetric and extended emission with a strong wavelength dependence. We propose geometric models based on an extended and distorted wind and/or a high density gas stream, which has long been predicted to be present in this system. The observations show that optical interferometry is now able to resolve HMXBs at the spatial scale at which accretion takes place, and therefore probe the effects of the gravitational and radiation fields of the compact object on its environment.

We present the first optical observation at sub-milliarcsecond (mas) scale of the microquasar SS 433 obtained with the GRAVITY instrument on the VLT interferometer. The 3.5 hour exposure reveals a rich K-band spectrum dominated by hydrogen Br$\gamma\ $ and \ion{He}{i} lines, as well as (red-shifted) emission lines coming from the jets. The K-band continuum emitting region is dominated by a marginally resolved point source ($<$ 1 mas) embedded inside a diffuse background accounting for 10\% of the total flux. The jet line positions agree well with the ones expected from the jet kinematic model, an interpretation also supported by the consistent sign (i.e. negative/positive for the receding/approaching jet component) of the phase shifts observed in the lines. The significant visibility drop across the jet lines, together with the small and nearly identical phases for all baselines, point toward a jet that is offset by less than 0.5 mas from the continuum source and resolved in the direction of propagation, with a typical size of 2 mas. The jet position angle of $\sim$80$^{\circ}$ is consistent with the expected one at the observation date. Jet emission so close to the central binary system would suggest that line locking, if relevant to explain the amplitude and stability of the 0.26c jet velocity, operates on elements heavier than hydrogen. The Br$\gamma\ $ profile is broad and double peaked. It is better resolved than the continuum and the change of the phase signal sign across the line on all baselines suggests an East-West oriented geometry alike the jet direction and supporting a (polar) disk wind origin.

GRAVITY is a new instrument to coherently combine the light of the European Southern Observatory Very Large Telescope Interferometer to form a telescope with an equivalent 130 m diameter angular resolution and a collecting area of 200 m$^2$. The instrument comprises fiber fed integrated optics beam combination, high resolution spectroscopy, built-in beam analysis and control, near-infrared wavefront sensing, phase-tracking, dual beam operation and laser metrology [...]. This article gives an overview of GRAVITY and reports on the performance and the first astronomical observations during commissioning in 2015/16. We demonstrate phase tracking on stars as faint as m$_K$ ~ 10 mag, phase-referenced interferometry of objects fainter than m$_K$ ~ 15 mag with a limiting magnitude of m$_K$ ~ 17 mag, minute long coherent integrations, a visibility accuracy of better than 0.25 %, and spectro-differential phase and closure phase accuracy better than 0.5{\deg}, corresponding to a differential astrometric precision of better than 10 microarcseconds ({\mu}as). The dual-beam astrometry, measuring the phase difference of two objects with laser metrology, is still under commissioning. First observations show residuals as low as 50 {\mu}as when following objects over several months. We illustrate the instrument performance with the observations of archetypical objects for the different instrument modes. Examples include the Galactic Center supermassive black hole and its fast orbiting star S2 for phase referenced dual beam observations and infrared wavefront sensing, the High Mass X-Ray Binary BP Cru and the Active Galactic Nucleus of PDS 456 for few {\mu}as spectro-differential astrometry, the T Tauri star S CrA for a spectro-differential visibility analysis, {\xi} Tel and 24 Cap for high accuracy visibility observations, and {\eta} Car for interferometric imaging with GRAVITY.

We propose a simple ansatz for estimating the value of the numerical resistivity and the numerical viscosity of any Eulerian MHD code. We test this ansatz with the help of simulations of the propagation of (magneto)sonic waves, Alfven waves, and the tearing mode instability using the MHD code Aenus. By comparing the simu- lation results with analytical solutions of the resistive-viscous MHD equations and an empirical ansatz for the growth rate of tearing modes we measure the numerical viscosity and resistivity of Aenus. The comparison shows that the fast-magnetosonic speed and wavelength are the characteristic velocity and length, respectively, of the aforementioned (relatively simple) systems. We also determine the dependance of the numerical viscosity and resistivity on the time integration method, the spatial reconstruction scheme and (to a lesser extent) the Riemann solver employed in the simulations. From the measured results we infer the numerical resolution (as a function of the spatial reconstruction method) required to properly resolve the growth and saturation level of the magnetic field amplified by the magnetorotational instability in the post-collapsed core of massive stars. Our results show that it is to the best advantage to resort to ultra-high order methods (e.g., 9th-order Monotonicity Preserving method) to tackle this problem properly, in particular in three dimensional simulations.

We present a new fourth-order, finite-volume hydrodynamics code named Apsara. The code employs a high-order, finite-volume method for mapped coordinates with extensions for nonlinear hyperbolic conservation laws. Apsara can handle arbitrary structured curvilinear meshes in three spatial dimensions. The code has successfully passed several hydrodynamic test problems, including the advection of a Gaussian density profile and a nonlinear vortex and the propagation of linear acoustic waves. For these test problems, Apsara produces fourth-order accurate results in case of smooth grid mappings. The order of accuracy is reduced to first-order when using the nonsmooth circular grid mapping. When applying the high-order method to simulations of low-Mach number flows, for example, the Gresho vortex and the Taylor-Green vortex, we discover that Apsara delivers superior results to codes based on the dimensionally split, piecewise parabolic method (PPM) widely used in astrophysics. Hence, Apsara is a suitable tool for simulating highly subsonic flows in astrophysics. In the first astrophysical application, we perform implicit large eddy simulations (ILES) of anisotropic turbulence in the context of core collapse supernova (CCSN) and obtain results similar to those previously reported.

Stochastic differential equations describe well many physical, biological and sociological systems, despite the simplification often made in their derivation. Here the usage of simple stochastic differential equations to characterize and classify complex dynamical systems is proposed within a Bayesian framework. To this end, we develop a dynamic system classifier (DSC). The DSC first abstracts training data of a system in terms of time dependent coefficients of the descriptive stochastic differential equation. Thereby the DSC identifies unique correlation structures within the training data. For definiteness we restrict the presentation of DSC to oscillation processes with a time dependent frequency {\omega}(t) and damping factor {\gamma}(t). Although real systems might be more complex, this simple oscillator captures many characteristic features. The {\omega} and {\gamma} timelines represent the abstract system characterization and permit the construction of efficient signal classifiers. Numerical experiments show that such classifiers perform well even in the low signal-to-noise regime.

Recent numerical simulations suggest that the model by Zahn (1992, A&A, 265, 115) for the turbulent mixing of chemical elements due to differential rotation in stellar radiative zones is valid. We investigate the robustness of this result with respect to the numerical configuration and Reynolds number of the flow. We compare results from simulations performed with two different numerical codes, including one that uses the shearing-box formalism. We also extensively study the dependence of the turbulent diffusion coefficient on the turbulent Reynolds number. The two numerical codes used in this study give consistent results. The turbulent diffusion coefficient is independent of the size of the numerical domain if at least three large turbulent structures fit in the box. Generally, the turbulent diffusion coefficient depends on the turbulent Reynolds number. However, our simulations suggest that an asymptotic regime is obtained when the turbulent Reynolds number is larger than $10^3$. Shear mixing in the regime of small P\'eclet numbers can be investigated numerically both with shearing-box simulations and simulations using explicit forcing. Our results suggest that Zahn's model is valid at large turbulent Reynolds numbers.

Whether the magnetorotational instability (MRI) can amplify initially weak magnetic fields to dynamically relevant strengths in core collapse supernovae is still a matter of active scientific debate. Recent numerical studies have shown that the first phase of MRI growth dominated by channel flows is terminated by parasitic instabilities of the Kelvin-Helmholtz type that disrupt MRI channel flows and quench further magnetic field growth. However, it remains to be prop- erly assessed by what factor the initial magnetic field can be amplified and how it depends on the initial field strength and the amplitude of the perturbations. Different termination criteria leading to different estimates of the amplification factor were proposed within the parasitic model. To determine the amplification factor and test which criterion is a better predictor of the MRI termination, we perform three-dimensional shearing-disc and shearing-box simula- tions of a region close to the surface of a differentially rotating proto-neutron star in non-ideal MHD with two different numerical codes. We find that independently of the initial magnetic field strength, the MRI channel modes can amplify the magnetic field by, at most, a factor of 100. Under the conditions found in proto-neutron stars a more realistic value for the mag- netic field amplification is of the order of 10. This severely limits the role of the MRI channel modes as an agent amplifying the magnetic field in proto-neutron stars starting from small seed fields. A further amplification should therefore rely on other physical processes, such as for example an MRI-driven turbulent dynamo.

We study the effect of superfluidity on torsional oscillations of highly magnetised neutron stars (magnetars) with a microphysical equation of state by means of two-dimensional, magnetohydrodynamical- elastic simulations. The superfluid properties of the neutrons in the neutron star core are treated in a parametric way in which we effectively decouple part of the core matter from the oscillations. Our simulations confirm the existence of two groups of oscillations, namely continuum oscillations that are confined to the neutron star core and are of Alfv\'enic character, and global oscillations with constant phase and that are of mixed magneto-elastic type. The latter might explain the quasi-periodic oscillations observed in magnetar giant flares, since they do not suffer from the additional damping mechanism due to phase mixing, contrary to what happens for continuum oscillations. However, we cannot prove rigorously that the coherent oscillations with constant phase are normal modes. Moreover, we find no crustal shear modes for the magnetic field strengths typical for magnetars.We provide fits to our numerical simulations that give the oscillation frequencies as functions of magnetic field strength and proton fraction in the core.

We study the influence of numerical methods and grid resolution on the termination of the magnetorotational instability (MRI) by means of parasitic instabilities in three-dimensional shearing-disc simulations reproducing typical conditions found in core-collapse supernovae. Whether or not the MRI is able to amplify weak magnetic fields in this context strongly depends, among other factors, on the amplitude at which its growth terminates. The qualitative results of our study do not depend on the numerical scheme. In all our models, MRI termination is caused by Kelvin-Helmholtz instabilities, consistent with theoretical predictions. Quantitatively, however, there are differences, but numerical convergence can be achieved even at relatively low grid resolutions if high-order reconstruction methods are used.

We investigate a method to construct parametrized progenitor models for core-collapse supernova simulations. Different from all modern core-collapse supernova studies, which rely on progenitor models from stellar evolution calculations, we follow the methodology of Baron & Cooperstein (1990) to construct initial models. Choosing parametrized spatial distributions of entropy and electron fraction as a function of mass coordinate and solving the equation of hydrostatic equilibrium, we obtain the initial density structures of our progenitor models. First, we calculate structures with parameters fitting broadly the evolutionary model s11.2 of Woosley et al. (2002). We then demonstrate the reliability of our method by performing general relativistic hydrodynamic simulations in spherical symmetry with the isotropic diffusion source approximation to solve the neutrino transport. Our comprehensive parameter study shows that initial models with a small central entropy ($\lesssim 0.4\,k_B$ nucleon$^{-1}$) can explode even in spherically symmetric simulations. Models with a large entropy ($\gtrsim 6\,k_B$ nucleon$^{-1}$) in the Si/O layer have a rather large explosion energy ($\sim 4\times 10^{50}$ erg) at the end of the simulations, which is still rapidly increasing.

The magnetorotational instability (MRI) can be a powerful mechanism amplifying the magnetic field in core-collapse supernovae. Whether initially weak magnetic fields can be amplified by this instability to dynamically relevant strengths is still a matter of debate. One of the main uncertainties concerns the process that terminates the growth of the instability. Parasitic instabilities of both Kelvin-Helmholtz and tearing-mode type have been suggested to play a crucial role in this process, disrupting MRI channel flows and quenching magnetic field amplification. We perform two-dimensional and three-dimensional sheering-disc simulations of a differentially rotating protoneutron star layer in non-ideal magnetohydrodynamics with unprecedented high numerical accuracy, finding that Kelvin-Helmholtz parasitic modes dominate tearing modes in the regime of large hydrodynamic and magnetic Reynolds numbers, as encountered close to the surface of protoneutron stars. They also determine the maximum magnetic field stress achievable during the exponential growth of the MRI. Our results are consistent with the theory of parasitic instabilities based on a local stability analysis. To simulate the Kelvin-Helmholtz instabilities properly, a very high numerical resolution is necessary. Using ninth-order spatial reconstruction schemes, we find that at least eight grid zones per MRI channel are necessary to simulate the growth phase of the MRI and reach an accuracy of ~10 per cent in the growth rate, while more than ~60 zones per channel are required to achieve convergent results for the value of the magnetic stress at MRI termination.

Observations in all electromagnetic bands show that many supernova remnants (SNRs) have a very aspherical shape. This can be the result of asymmetries in the supernova explosion or a clumpy circumstellar medium. We study the generation of inhomogeneities and the mixing of elements arising from these two sources in multidimensional hydrodynamic simulations of the propagation of a supernova blast wave into a cloudy environment. We model a specific SNR, Vela Jr (RX J0852.0-4622). By comparing our results with recent observations, we can constrain the properties of the explosion. We find that a very energetic explosion of several 10^{51} erg occurring roughly about 800 years ago is consistent with the shape and emission of the SNR, as well as a supernova with an energy closer to the canonical value of 10^{51} erg a few thousand years ago.

We couple internal torsional, magneto-elastic oscillations of highly magnetized neutron stars (magnetars) to their magnetospheres. The corresponding axisymmetric perturbations of the external magnetic field configuration evolve as a sequence of linear, force-free equilibria that are completely determined by the background magnetic field configuration and by the perturbations of the magnetic field at the surface. The perturbations are obtained from simulations of magneto-elastic oscillations in the interior of the magnetar. While such oscillations can excite travelling Alfv\'en waves in the exterior of the star only in a very limited region close to the poles, they still modulate the near magnetosphere by inducing a time-dependent twist between the foot-points of closed magnetic field lines that exit the star at a polar angle $\gtrsim 0.19\,$rad. Moreover, we find that for a dipole-like background magnetic field configuration the magnetic field modulations in the magnetosphere, driven by internal oscillations, can only be symmetric with respect to the equator. This is in agreement with our previous findings, where we interpreted the observed quasi-periodic oscillations in the X-ray tail of magnetar bursts as driven by the family of internal magneto-elastic oscillations with symmetric magnetic field perturbations.

In this paper we report on what we believe is the first successful implementation of relativistic hydrodynamics, coupled to dynamical spacetimes, in spherical polar coordinates without symmetry assumptions. We employ a high-resolution shock-capturing scheme, which requires that the equations be cast in flux-conservative form. One example of such a form is the :Valencia" formulation, which has been adopted in numerous applications, in particular in Cartesian coordinates. Here we generalize this formulation to allow for a reference-metric approach, which provides a natural framework for calculations in curvilinear coordinates. In spherical polar coordinates, for example, it allows for an analytical treatment of the singular r and sin(\theta) terms that appear in the equations. We experiment with different versions of our generalized Valencia formulation in numerical implementations of relativistic hydrodynamics for both fixed and dynamical spacetimes. We consider a number of different tests -- non-rotating and rotating relativistic stars, as well as gravitational collapse to a black hole -- to demonstrate that our formulation provides a promising approach to performing fully relativistic astrophysics simulations in spherical polar coordinates.

Our numerical simulations show that axisymmetric, torsional, magneto-elastic oscillations of magnetars with a superfluid core can explain the whole range of observed quasi-periodic oscillations (QPOs) in the giant flares of soft gamma-ray repeaters. There exist constant phase, magneto-elastic QPOs at both low (f<150 Hz) and high frequencies (f>500 Hz), in full agreement with observations. The range of magnetic field strengths required to match the observed QPO frequencies agrees with that from spin-down estimates. These results strongly suggest that neutrons in magnetar cores are superfluid.

We study magneto-elastic oscillations of highly magnetized neutron stars (magnetars) which have been proposed as an explanation for the quasi-periodic oscillations (QPOs) appearing in the decaying tail of the giant flares of soft gamma-ray repeaters (SGRs). We extend previous studies by investigating various magnetic field configurations, computing the Alfv\'en spectrum in each case and performing magneto-elastic simulations for a selected number of models. By identifying the observed frequencies of 28 Hz (SGR 1900+14) and 30 Hz (SGR 1806-20) with the fundamental Alfv\'en QPOs, we estimate the required surface magnetic field strength. For the magnetic field configurations investigated (dipole-like poloidal, mixed toroidal-poloidal with a dipole-like poloidal component and a toroidal field confined to the region of field lines closing inside the star, and for poloidal fields with an additional quadrupole-like component) the estimated dipole spin-down magnetic fields are between 8x10^14 G and 4x10^15 G, in broad agreement with spin-down estimates for the SGR sources producing giant flares. A number of these models exhibit a rich Alfv\'en continuum revealing new turning points which can produce QPOs. This allows one to explain most of the observed QPO frequencies as associated with magneto-elastic QPOs. In particular, we construct a possible configuration with two turning points in the spectrum which can explain all observed QPOs of SGR 1900+14. Finally, we find that magnetic field configurations which are entirely confined in the crust (if the core is assumed to be a type I superconductor) are not favoured, due to difficulties in explaining the lowest observed QPO frequencies (f<30 Hz).

In the absence of symmetry assumptions most numerical relativity simulations adopt Cartesian coordinates. While Cartesian coordinates have some desirable properties, spherical polar coordinates appear better suited for certain applications, including gravitational collapse and supernova simulations. Development of numerical relativity codes in spherical polar coordinates has been hampered by the need to handle the coordinate singularities at the origin and on the axis, for example by careful regularization of the appropriate variables. Assuming spherical symmetry and adopting a covariant version of the BSSN equations, Montero and Cordero-Carri\'on recently demonstrated that such a regularization is not necessary when a partially implicit Runge-Kutta (PIRK) method is used for the time evolution of the gravitational fields. Here we report on an implementation of the BSSN equations in spherical polar coordinates without any symmetry assumptions. Using a PIRK method we obtain stable simulations in three spatial dimensions without the need to regularize the origin or the axis. We perform and discuss a number of tests to assess the stability, accuracy and convergence of the code, namely weak gravitational waves, "hydro-without-hydro" evolutions of spherical and rotating relativistic stars in equilibrium, and single black holes.

We study supernova neutrino flux variations in the IceCube detector, using 3D models based on a simplified neutrino transport scheme. The hemispherically integrated neutrino emission shows significantly smaller variations compared with our previous study of 2D models, largely because of the reduced SASI activity in this set of 3D models which we interpret as a pessimistic extreme. For the studied cases, intrinsic flux variations up to about 100 Hz frequencies could still be detected in a supernova closer than about 2 kpc.

By means of two dimensional, general-relativistic, magneto-hydrodynamical simulations we investigate the oscillations of magnetized neutron star models (magnetars) including the description of an extended solid crust. The aim of this study is to understand the origin of the QPOs observed in the giant flares of SGRs. We confirm the existence of three different regimes: (a) a weak magnetic field regime B<5 x 10^13 G, where crustal shear modes dominate the evolution; (b) a regime of intermediate magnetic fields 5 x 10^13 G<B< 10^15 G, where Alfv\'en QPOs are mainly confined to the core of the neutron star and the crustal shear modes are damped very efficiently; and (c) a strong field regime B>10^15 G, where magneto-elastic oscillations reach the surface and approach the behavior of purely Alfv\'en QPOs. When the Alfv\'en QPOs are confined to the core of the neutron star, we find qualitatively similar QPOs as in the absence of a crust. The lower QPOs associated with the closed field lines of the dipolar magnetic field configuration are reproduced as in our previous simulations without crust, while the upper QPOs connected to the open field lines are displaced from the polar axis. Additionally, we observe a family of edge QPOs. Our results do not leave much room for a crustal-mode interpretation of observed QPOs in SGR giant flares, but can accommodate an interpretation of these observations as originating from Alfv\'en-like, global, turning-point QPOs in models with dipolar magnetic field strengths in the narrow range of 5 x 10^15 G < B < 1.4 x 10^16 G. This range is somewhat larger than estimates for magnetic field strengths in known magnetars. The discrepancy may be resolved in models including a more complicated magnetic field structure or with models taking superfluidity of the neutrons and superconductivity of the protons in the core into account.

Time-dependent and direction-dependent neutrino and gravitational-wave (GW) signatures are presented for a set of 3D hydrodynamic models of parametrized, neutrino-driven supernova explosions of non-rotating 15 and 20 solar mass stars. We employ an approximate treatment of neutrino transport. Due to the excision of the high-density core of the proto-neutron star and the use of an axis-free overset grid, the models can be followed from the post-bounce accretion phase for more than one second without imposing any symmetry restrictions. GW and neutrino emission exhibit the generic time-dependent features known from 2D models. Non-radial hydrodynamic mass motions in the accretion layer and their interaction with the outer layers of the proto-neutron star together with anisotropic neutrino emission give rise to a GW signal with an amplitude of ~5-20 cm and frequencies 100--500 Hz. The GW emission from mass motions reaches a maximum before the explosion sets in. Afterwards the GW signal exhibits a low-frequency modulation, in some cases describing a quasi-monotonic growth, associated with the non-spherical expansion of the explosion shock wave and the large-scale anisotropy of the escaping neutrino flow. Variations of the mass-quadrupole moment due to convective activity inside the nascent neutron star contribute a high-frequency component to the GW signal during the post-explosion phase. The GW signals exhibit strong variability between the two polarizations, different explosion simulations and different observer directions, and does not possess any template character. The neutrino emission properties show fluctuations over the neutron star surface on spatial and temporal scales that reflect the different types of non-spherical mass motions. The modulation amplitudes of the measurable neutrino luminosities and mean energies are significantly smaller than predicted by 2D simulations.

During most stages of stellar evolution the nuclear burning of lighter to heavier elements results in a radial composition profile which is stabilizing against buoyant acceleration, with light material residing above heavier material. However, under some circumstances, such as off-center ignition, the composition profile resulting from nuclear burning can be destabilizing, and characterized by an outwardly increasing mean molecular weight. The potential for instabilities under these circumstances, and the consequences that they may have on stellar structural evolution, remain largely unexplored. In this paper we study the development and evolution of instabilities associated with unstable composition gradients in regions which are initially stable according to linear Schwarzschild and Ledoux criteria. In particular, we explore the mixing taking place under various conditions with multi-dimensional hydrodynamic convection models based on stellar evolutionary calculations of the core helium flash in a 1.25 \Msun star, the core carbon flash in a 9.3\,\Msun star, and of oxygen shell burning in a star with a mass of 23\,\Msun. The results of our simulations reveal a mixing process associated with regions having outwardly increasing mean molecular weight that reside below convection zones. The mixing is not due to overshooting from the convection zone, nor is it due directly to thermohaline mixing which operates on a timescale several orders of magnitude larger than the simulated flows. Instead, the mixing appears to be due to the presence of a wave field induced in the stable layers residing beneath the convection zone which enhances the mixing rate by many orders of magnitude and allows a thermohaline type mixing process to operate on a dynamical, rather than thermal, timescale. We discuss our results in terms of related laboratory phenomena and associated theoretical developments.

We extend a general-relativistic ideal magneto-hydrodynamical code to include the effects of elasticity. Using this numerical tool we analyse the magneto-elastic oscillations of highly magnetised neutron stars (magnetars). In simulations without magnetic field we are able to recover the purely crustal shear oscillations within an accuracy of about a few per cent. For dipole magnetic fields between 5 x 10^13 and 10^15 G the Alfv\'en oscillations become modified substantially by the presence of the crust. Those quasi-periodic oscillations (QPOs) split into three families: Lower QPOs near the equator, Edge QPOs related to the last open field line and Upper QPOs at larger distance from the equator. Edge QPOs are called so because they are related to an edge in the corresponding Alfv\'en continuum. The Upper QPOs are of the same kind, while the Lower QPOs are turning-point QPOs, related to a turning point in the continuous spectrum.

In a realistic model of magneto-elastic oscillations in magnetars, we find that crustal shear oscillations, often invoked as an explanation of quasi-periodic oscillations (QPOs) seen after giant flares in soft gamma-ray repeaters (SGRs), are damped by resonant absorption on timescales of at most 0.2s, for a lower limit on the dipole magnetic field strength of 5 x 10^13 G. At higher magnetic field strengths (typical in magnetars) the damping timescale is even shorter, as anticipated by earlier toy-models. We have investigated a range of equations of state and masses and if magnetars are dominated by a dipole magnetic field, our findings exclude torsional shear oscillations of the crust from explaining the observed low-frequency QPOs. In contrast, we find that the Alfv\'en QPO model is a viable explanation of observed QPOs, if the dipole magnetic field strength exceeds a minimum strength of about several times 10^14 G to 10^15 G. Then, Alfv\'en QPOs are no longer confined to the fluid core, but completely dominate in the crust region and have a maximum amplitude at the surface of the star.

We present the first stationary, axisymmetric neutron star models with meridional circulation in general relativity. For that purpose, we developed GRNS, a new code based on a fixed point iteration. We find a two-dimensional set of meridional circulation modes, which differ by the number of vortices in the stream lines of the neutron star fluid. For expected maximal meridional circulation velocities of about 1000 km/s, the vortices cause surface deformations of about a percent. The deformations depend on the shape of the vortices close to the surface and increase with the meridional circulation velocity. We also computed models of rotating neutron stars with meridional circulation, where neither the surface rotates nor does the rotation velocity exceed the circulation velocity.

Context. Global MHD simulations show Kelvin-Helmholtz (KH) instabilities at the contact surface of two merging neutron stars. That region has been identified as the site of efficient amplification of magnetic fields. However, these global simulations, due to numerical limitations, were unable to determine the saturation level of the field strength, and thus the possible back-reaction of the magnetic field onto the flow. Aims. We investigate the amplification of initially weak fields in KH unstable shear flows, and the back-reaction of the field onto the flow. Methods. We use a high-resolution ideal MHD code to perform 2D and 3D local simulations of shear flows. Results. In 2D, the magnetic field is amplified in less than 0.01ms until it reaches locally equipartition with the kinetic energy. Subsequently, it saturates due to resistive instabilities that disrupt the KH vortex and decelerate the shear flow on a secular time scale. We determine scaling laws of the field amplification with the initial field strength and the grid resolution. In 3D, this hydromagnetic mechanism may be dominated by purely hydrodynamic instabilities limiting the amplification. We find maximum magnetic fields of 10^16 G locally, and r.m.s. maxima within the box of 10^15 G. However, such strong fields exist only for a short period. In the saturated state, the magnetic field is mainly oriented parallel to the shear flow for strong initial fields, while weaker initial fields tend to lead to a more balanced distribution of the field energy. In all models the flow shows small-scale features. The magnetic field is at most in equipartition with the decaying shear flow. (abridged)

A type of overlapping grid in spherical coordinates called the Yin-Yang grid is successfully implemented into a 3D version of the explicit Eulerian grid-based code PROMETHEUS including self-gravity. The modified code successfully passed several standard hydrodynamic tests producing results which are in very good agreement with analytic solutions. Moreover, the solutions obtained with the Yin-Yang grid exhibit no peculiar behaviour at the boundary between the two grid patches. The code has also been successfully used to model astrophysically relevant situations, namely equilibrium polytropes, a Taylor-Sedov explosion, and Rayleigh-Taylor instabilities. According to our results, the usage of the Yin-Yang grid greatly enhances the suitability and efficiency of 3D explicit Eulerian codes based on spherical polar coordinates for astrophysical flows.

We assess the importance of the magneto-rotational instability in core-collapse supernovae by an analysis of the growth rates of unstable modes in typical post-collapse systems and by numerical simulations of simplified models. The interplay of differential rotation and thermal stratification defines different instability regimes which we confirm in our simulations. We investigate the termination of the growth of the MRI by parasitic instabilities, establish scaling laws characterising the termination amplitude, and study the long-term evolution of the saturated turbulent state.

Possible effects of magnetic fields in core collapse supernovae rely on an efficient amplification of the weak pre-collapse fields. It has been suggested that the magneto-rotational instability (MRI) leads to rapid field growth. Although MRI studies exist for accretion discs, the application of their results to core collapse supernovae is inhibited as the physics of supernova cores is substantially different from that of accretion discs. We address the problem of growth and saturation of the MRI by means of semi-global simulations, which combine elements of global and local simulations by taking the presence of global background gradients into account and using a local computational grid. We analyze the dispersion relation of the MRI to identify different regimes of the instability. This analysis is complemented by simulations, where we consider a local computational box rotating at sub-Keplerian velocity, and where we allow for a radial entropy gradient. We identify six regimes of the MRI depending on the ratio of the entropy and angular velocity gradient. Our numerical models confirm the instability criteria and growth rates for all relevant regimes. The MRI grows exponentially within milliseconds the flow and magnetic field geometries being dominated by channel flows. The MRI growth ceases once the channels are disrupted by resistive instabilities (due to finite numerical conductivity), and MHD turbulence sets in. From an analysis of the growth rates of the resistive instabilities, we deduce scaling laws for the termination amplitude of the MRI which agree well with our numerical models. We determine the dependence of the development of coherent flow structures in the saturated state on the aspect ratio of the simulation boxes. [abridged]

We present a new numerical code which solves the general relativistic magneto-hydrodynamics (GRMHD) equations coupled to the Einstein equations for the evolution of a dynamical spacetime within the conformally-flat approximation. This code has been developed with the main objective of studying astrophysical scenarios in which both, high magnetic fields and strong gravitational fields appear, such as the magneto-rotational collapse of stellar cores, the collapsar model of GRBs, and the evolution of neutron stars. The code is based on an existing and thoroughly tested purely hydrodynamics code and on its extension to accommodate weakly magnetized fluids (passive magnetic field approximation). The numerical code we present here is based on high-resolution shock-capturing schemes to solve the GRMHD equations together with the flux constraint transport method to ensure the solenoidal condition of the magnetic field. Since the astrophysical applications envisaged do not deviate much from spherical symmetry, the conformal flatness condition approximation is used for the formulation of the Einstein equations. In addition, the code can handle several equations of state, from simple analytical expressions to microphysical tabulated ones. In this paper we present stringent tests of our new GRMHD numerical code, which show its ability to handle all aspects appearing in the astrophysical scenarios for which the code is intended, namely relativistic shocks, highly magnetized fluids, and equilibrium configurations of magnetized neutron stars. As an application, magneto-rotational core collapse simulations of a realistic progenitor are presented, comparing the results with our previous finding in the passive magnetic field approximation.

We performed cosmological, magneto-hydrodynamical simulations to follow the evolution of magnetic fields in galaxy clusters, exploring the possibility that the origin of the magnetic seed fields are galactic outflows during the star-burst phase of galactic evolution. To do this we coupled a semi-analytical model for magnetized galactic winds as suggested by \citet{2006MNRAS.370..319B} to our cosmological simulation. We find that the strength and structure of magnetic fields observed in galaxy clusters are well reproduced for a wide range of model parameters for the magnetized, galactic winds and do only weakly depend on the exact magnetic structure within the assumed galactic outflows. Although the evolution of a primordial magnetic seed field shows no significant differences to that of galaxy clusters fields from previous studies, we find that the magnetic field pollution in the diffuse medium within filaments is below the level predicted by scenarios with pure primordial magnetic seed field. We therefore conclude that magnetized galactic outflows and their subsequent evolution within the intra-cluster medium can fully account for the observed magnetic fields in galaxy clusters. Our findings also suggest that measuring cosmological magnetic fields in low-density environments such as filaments is much more useful than observing cluster magnetic fields to infer their possible origin.

We show by two-dimensional and first three-dimensional simulations of neutrino-driven supernova explosions that low (l=1,2) modes can dominate the flow pattern in the convective postshock region on timescales of hundreds of milliseconds after core bounce. This can lead to large global anisotropy of the supernova explosion and pulsar kicks in excess of 500 km/s.