Search

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

Jean-Paul Zahn’s formalism for vertical shear mixing is used in several stellar evolution codes, but the physics of the shear instability in stellar radiative zones is still not completely understood. Over the last few years, numerical simulations have provided new constraints on the shear instability, including the effect of thermal diffusion and chemical stratification. We present here new simulations that show the effect of viscosity on the vertical turbulent transport due to the shear instability.

Six binary-merger progenitors of Supernova 1987A (SN 1987A) with properties close to those of the blue supergiant Sanduleak -69 202 are exploded by neutrino heating and evolved until long after shock breakout in three dimensions (3D), and continued for light-curve calculations in spherical symmetry. Our results confirm previous findings for single-star progenitors: (1) 3D neutrino-driven explosions with SN 1987A-like energies synthesize Ni-56 masses consistent with the radioactive light-curve tail; (2) hydrodynamic models mix hydrogen inward to minimum velocities below 40 km/s compatible with spectral observations of SN 1987A; and (3) for given explosion energy the efficiency of outward radioactive Ni-56 mixing depends mainly on high growth factors of Rayleigh-Taylor instabilities at the (C+O)/He and He/H composition interfaces and a weak interaction of fast plumes with the reverse shock occurring below the He/H interface. All binary-merger models possess presupernova radii matching the photometric radius of Sanduleak -69 202 and a structure of the outer layers allowing them to reproduce the observed initial luminosity peak in the first about 7 days. Models that mix about 0.5 Msun of hydrogen into the He-shell and exhibit strong outward mixing of Ni-56 with maximum velocities exceeding the 3000 km/s observed for the bulk of ejected Ni-56 have light-curve shapes in good agreement with the dome of the SN 1987A light curve. A comparative analysis of the best representatives of our 3D neutrino-driven explosion models of SN 1987A based on single-star and binary-merger progenitors reveals that only one binary model fulfills all observational constraints, except one., 31 pages, 14 figures, 6 tables. Accepted for publication in ApJ

The extent of mixed regions around convective zones is one of the biggest uncertainties in stellar evolution. 1D overshooting descriptions introduce a free parameter ($f_{ov}$) that is in general not well constrained from observations. Especially in small central convective regions the value is highly uncertain due to its tight connection to the pressure scale height. Long-term multi-dimensional hydrodynamic simulations can be used to study the size of the overshooting region and the involved mixing processes. Here we show how one can calibrate an overshooting parameter by performing 2D Maestro simulations of Zero-Age-Main-Sequence stars ranging from $1.3$ to $3.5 M_\odot$. The simulations cover the convective cores of the stars and a large fraction of the surrounding radiative envelope. We follow the convective flow for at least 20 convective turnover times, while the longest simulation covers 430 turnover time scales. This allows us to study how the mixing as well as the convective boundary evolve with time, and how the resulting entrainment can be interpreted in terms of overshooting parameters. We find that increasing the overshooting parameter $f_{ov}$ beyond a certain value in the initial model of our simulations, changes the mixing behaviour completely. This result can be used to put limits on the overshooting parameter. We find $0.010 < f_{ov} < 0.017$ to be in good agreement with our simulations of a $3.5 M_\odot$ mass star. We also identify a diffusive mixing component due to internal gravity waves (IGW) that is active throughout the convectively stable layer, but likely overestimated in our simulations. Furthermore, applying our calibration method to simulations of less massive stars suggests a need for a mass-dependent overshooting description where the mixing in terms of the pressure scale height is reduced for small convective cores., Accepted for publication in A&A

Adaptation to prolonged flooding was investigated using cuttings of two tree species from the Central Amazon white-water floodplain (V�rzea). Morphological features and oxygen distribution patterns were correlated with metabolic changes under hypoxia, such as alterations in alcohol dehydrogenase (ADH) activity and adenylate energy charge (AEC) of root cells. Salix martiana (Leyb.) was able to react to hypoxic growth conditions with formation of adventitious roots rich in lysigenous aerenchyma, which facilitates root aeration by longitudinal oxygen transport and rhizosphere oxidation by radial oxygen loss (ROL). The oxygen concentration on the surface of adventitious roots of S. martiana reached 2-3 mg O2 L–1. The low resistance to gas exchange in Salix roots was reflected by low ADH activities, which ranged between 0.03-0.1 μmol NADH mg –1 min–1, and AEC values of 0.8-1 under hypoxic conditions. Adventitious roots were also formed by Tabernaemontana juruana ([Markgr.] Schumann ex. J.F. Macbride) during growth under low-oxygen conditions, although at a later stage. The gas-space continuum in roots of T. juruana was less pronounced, resulting in a 10-fold lower oxygen concentration in the root cortex under oxygen stress compared with adventitious roots of Salix. The lower oxygen content was reflected in 6-fold higher ADH activities and decreased AEC values. ROL occurred only at the non-suberized root tip, suggesting that the suberized hypodermis functions as a barrier against gas exchange between the root and the rhizosphere. These findings indicate that different strategies of adaptation to low oxygen levels are realized in the two species under investigation that occur naturally in the same ecosystem but inhabit different elevation sites.

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., Comment: 15 pages, 13 figures, submitted to MNRAS

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., Comment: 7 pages, 4 figures, proceedings of stars2017 and 2017smfns

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., Comment: 8 pages, 5 figures; accepted for publication in MNRAS

We present a new method for solving the three-dimensional gravitational potential of a density field on the Yin-Yang grid. Our algorithm is based on a multipole decomposition and completely symmetric with respect to the two Yin-Yang grid patches. It is particularly efficient on distributed-memory machines with a large number of compute tasks, because the amount of data being explicitly communicated is minimized. All operations are performed on the original grid without the need for interpolating data onto an auxiliary spherical mesh., 8 pages, 4 figures; two minor additions after refereeing; accepted by ApJ

Context. η Car is one of the most intriguing luminous blue variables in the Galaxy. Observations and models of the X-ray, ultraviolet, optical, and infrared emission suggest a central binary in a highly eccentric orbit 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. However, the morphology of the cavity and its physical properties are not yet fully constrained. Aims. We aim to trace the inner ∼5–50 au structure of η Car’s wind-wind interaction, as seen through Brγ and, for the first time, through the He I 2s-2p line. Methods. We have used spectro-interferometric observations with the K-band beam-combiner GRAVITY at the VLTI. The analyses of the data include (i) parametrical model-fitting to the interferometric observables, (ii) a CMFGEN model of the source’s spectrum, and (iii) interferometric image reconstruction. Results. Our geometrical modeling of the continuum data allows us to estimate its FWHM angular size close to 2 mas and an elongation ratio ϵ = 1.06 ± 0.05 over a PA = 130° ± 20°. Our CMFGEN modeling of the spectrum helped us to confirm that the role of the secondary should be taken into account to properly reproduce the observed Brγ and He I lines. Chromatic images across the Brγ line reveal a southeast arc-like feature, possibly associated to the hot post-shocked winds flowing along the cavity wall. The images of the He I 2s-2p line served to constrain the 20 mas (∼50 au) 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 Brγ 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 η Car.

With the same method as used previously, we investigate neutrino-driven explosions of a larger sample of blue supergiant models. The larger sample includes three new presupernova stars. The results are compared with light-curve observations of the peculiar type IIP SN 1987A. The explosions were modeled in 3D with the neutrino-hydrodynamics code PROMETHEUS-HOTB, and light-curve calculations were performed in spherical symmetry with the radiation-hydrodynamics code CRAB. Our results confirm the basic findings of the previous work: 3D neutrino-driven explosions with SN 1987A-like energies synthesize an amount of Ni-56 that is consistent with the radioactive tail of the light curve. Moreover, the models mix hydrogen inward to minimum velocities below 400 km/s as required by spectral observations. Hydrodynamic simulations with the new progenitor models, which possess smaller radii than the older ones, show much better agreement between calculated and observed light curves in the initial luminosity peak and during the first 20 days. A set of explosions with similar energies demonstrated that a high growth factor of Rayleigh-Taylor instabilities at the (C+O)/He composition interface combined with a weak interaction of fast Rayleigh-Taylor plumes, where the reverse shock occurs below the He/H interface, provides a sufficient condition for efficient outward mixing of Ni-56 into the hydrogen envelope. This condition is realized to the required extent only in one of the older stellar models, which yielded a maximum velocity of around 3000 km/s for the bulk of ejected Ni-56, but failed to reproduce the helium-core mass of 6 Msun inferred from the absolute luminosity of the presupernova star. We conclude that none of the single-star progenitor models proposed for SN 1987A to date satisfies all constraints set by observations. (Abridged), Comment: 16 pages, 9 figures, 3 tables; accepted for publication in Astron. Astrophys

The magneto-rotational instability (MRI) is considered to be a promising mechanism to amplify the magnetic field in fast rotating protoneutron stars. In contrast to accretion disks, radial buoyancy driven by entropy and lepton fraction gradients is expected to have a dynamical role as important as rotation and shear. We investigate the poorly known impact of buoyancy on the non-linear phase of the MRI, by means of three dimensional numerical simulations of a local model in the equatorial plane of a protoneutron star. The use of the Boussinesq approximation allows us to utilise a shearing box model with clean shearing periodic boundary conditions, while taking into account the buoyancy driven by radial entropy and composition gradients. We find significantly stronger turbulence and magnetic fields in buoyantly unstable flows. On the other hand, buoyancy has only a limited impact on the strength of turbulence and magnetic field amplification for buoyantly stable flows in the presence of a realistic thermal diffusion. The properties of the turbulence are, however, significantly affected in the latter case. In particular, the toroidal components of the magnetic field and of the velocity become even more dominant with respect to the poloidal ones. Furthermore, we observed in the regime of stable buoyancy the formation of long lived coherent structures such as channel flows and zonal flows. Overall, our results support the ability of the MRI to amplify the magnetic field significantly even in stably stratified regions of protoneutron stars., 22 pages, 15 figures, accepted for publication in MNRAS

Extremely strong magnetic fields of the order of $10^{15}\,{\rm G}$ are required to explain the properties of magnetars, the most magnetic neutron stars. Such a strong magnetic field is expected to play an important role for the dynamics of core-collapse supernovae, and in the presence of rapid rotation may power superluminous supernovae and hypernovae associated to long gamma-ray bursts. The origin of these strong magnetic fields remains, however, obscure and most likely requires an amplification over many orders of magnitude in the protoneutron star. One of the most promising agents is the magnetorotational instability (MRI), which can in principle amplify exponentially fast a weak initial magnetic field to a dynamically relevant strength. We describe our current understanding of the MRI in protoneutron stars and show recent results on its dependence on physical conditions specific to protoneutron stars such as neutrino radiation, strong buoyancy effects and large magnetic Prandtl number., Comment: 6 pages, 4 figures, proceedings of IAUS 331

International audience; We observe the high-mass X-ray binary (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 (${\rm{\Delta }}| V| \sim 5 \% $) and phase (${\rm{\Delta }}\phi \sim 2^\circ $) signatures are observed across the He i $2.059\,\mu {\rm{m}}$ and Brγ 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\buildrel{\circ}\over{.} 2$ corresponds to an astrometric precision of $\sim 2\,\mu \mathrm{as}$. We generalize expressions for image centroid displacements and variances in the marginally resolved limit of interferometry to spectrally resolved data, and use them to derive model-independent properties of the emission such as its asymmetry, extension, and 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 where accretion takes place, and therefore to probe the effects of the gravitational and radiation fields of the compact object on its environment.

Type II-plateau supernovae (SNe IIP) are the most numerous subclass of core-collapse SNe originating from massive stars. In the framework of the neutrino-driven explosion mechanism, we study the SN outburst properties for a red supergiant progenitor model and compare the corresponding light curves with observations of the ordinary Type IIP SN 1999em. Three-dimensional (3D) simulations of (parametrically triggered) neutrino-driven explosions are performed with the (explicit, finite-volume, Eulerian, multifluid hydrodynamics) code PROMETHEUS, using a presupernova model of a 15 Msun star as initial data. At approaching homologous expansion, the hydrodynamical and composition variables of the 3D models are mapped to a spherically symmetric configuration, and the simulations are continued with the (implicit, Lagrangian radiation-hydrodynamics) code CRAB to follow the blast-wave evolution during the SN outburst. Our 3D neutrino-driven explosion model with an explosion energy of about 0.5x10^51 erg produces Ni-56 in rough agreement with the amount deduced from fitting the radioactively powered light-curve tail of SN 1999em. The considered presupernova model, 3D explosion simulations, and light-curve calculations can explain the basic observational features of SN 1999em, except for those connected to the presupernova structure of the outer stellar layers. Our 3D simulations show that the distribution of Ni-rich matter in velocity space is asymmetric with a strong dipole component that is consistent with the observations of SN 1999em. The monotonic luminosity decline from the plateau to the radioactive tail in ordinary SNe IIP is a manifestation of the intense turbulent mixing at the He/H composition interface., Comment: 16 pages, 13 figures, 2 tables; added figure, discussions, and references; accepted for publication in ApJ

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/R1.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., Comment: 14 pages, 8 figures, 6 tables, submitted to MNRAS

We present a new numerical tool to calculate the emission of highly magnetized neutron stars (magnetars) and apply it to describe the quasi-periodic oscillations (QPOs) observed in magnetar giant flares. In previous work we have developed a model of magneto-elastic oscillations of magnetars that allows to reproduce the observed frequencies. These QPOs can couple to the star's exterior through the magnetic field and induce currents in the magnetosphere that provide scattering targets for resonant cyclotron scattering of the photons. The scattering is calculated with a Monte-Carlo approach and it is coupled to a code that calculates the momentum distribution of the charge carriers as an one-dimensional accelerator problem. As a first test of the method we calculate the modulation of the quiescent emission of the neutron star by the magneto-elastic QPOs for a prescribed momentum distribution of the charge carriers. The necessary amplitudes of the QPOs at the surface of the star to modulate the emission significantly are ≲1km. (© 2014 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)

The spatial and velocity distributions of nuclear species synthesized in the innermost regions of core-collapse supernovae can yield important clues about explosion asymmetries and the operation of the still disputed explosion mechanism. Recent observations of radioactive $^{44}$Ti with high-energy satellite telescopes (Nuclear Spectroscopic Telescope Array [NuSTAR], INTEGRAL) have measured gamma-ray line details, which provide direct evidence of large-scale explosion asymmetries in SN 1987A and in Cassiopeia A (Cas A) even by mapping of the spatial brightness distribution (NuSTAR). Here we discuss a 3D simulation of a neutrino-driven explosion, using a parameterized neutrino engine, whose $^{44}$Ti distribution is mostly concentrated in one hemisphere pointing opposite to the neutron star (NS) kick velocity. Both exhibit intriguing resemblance to the observed morphology of the Cas A remnant, although neither the progenitor nor the explosion was fine-tuned for a perfect match. Our results demonstrate that the asymmetries observed in this remnant can, in principle, be accounted for by a neutrino-driven explosion, and that the high $^{44}$Ti abundance in Cas A may be explained without invoking rapid rotation or a jet-driven explosion, because neutrino-driven explosions generically eject large amounts of high-entropy matter. The recoil acceleration of the NS is connected to mass ejection asymmetries and is opposite to the direction of the stronger explosion, fully compatible with the gravitational tugboat mechanism. Our results also imply that Cas A and SN 1987A could possess similarly "one-sided" Ti and Fe asymmetries, with the difference that Cas A is viewed from a direction with large inclination angle to the NS motion, whereas the NS in SN 1987A should have a dominant velocity component pointing toward us., 23 pages, 12 figures; revised version; discussion of new 3D analysis of 44Ti in Cas A and 3 new figures added; accepted by ApJ

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., Comment: 11 pages, 8 figures

During the common envelope (CE) phase, a giant star in a binary system overflows its Roche lobe and unstable mass transfer leads to a spiral-in of the companion, resulting in a close binary system or in a merger of the stellar cores. Dynamo processes during the CE phase have been proposed as a mechanism to generate magnetic fields that are important for forming magnetic white dwarfs (MWDs) and for shaping planetary nebulae. Here, we present the first magnetohydrodynamics simulations of the dynamical spiral-in during a CE phase. We find that magnetic fields are strongly amplified in the accretion stream around the $1M_\odot$ companion as it spirals into the envelope of a $2M_\odot$ RG. This leads to field strengths of 10 to 100 kG throughout the envelope after 120 d. The magnetic field amplification is consistent with being driven by the magnetorotational instability. The field strengths reached in our simulation make the magnetic field interesting for diagnostic purposes, but they are dynamically irrelevant. They are also too small to explain the formation of the highest fields found in MWDs, but may be relevant for luminous red novae, and detecting magnetic fields in these events would support the scenario as proposed here., 5 pages, 2 figures, 1 table. Accepted for publication in MNRAS Letters

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., Comment: 16 pages, 12 figures, accepted by MNRAS

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., 12 pages, 11 figures, 5 tables; accepted for publication in MNRAS

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., 19 pages, 6 figures, accepted for publication in MNRAS

We present gravitational wave (GW) signal predictions from four 3D multi-group neutrino hydrodynamics simulations of core-collapse supernovae of progenitors with 11.2 Msun, 20 Msun, and 27 Msun. GW emission in the pre-explosion phase strongly depends on whether the post-shock flow is dominated by the standing accretion shock instability (SASI) or convection and differs considerably from 2D models. SASI activity produces a strong signal component below 250 Hz through asymmetric mass motions in the gain layer and a non-resonant coupling to the proto-neutron star (PNS). Both convection- and SASI-dominated models show GW emission above 250 Hz, but with considerably lower amplitudes than in 2D. This is due to a different excitation mechanism for high-frequency l=2 motions in the PNS surface, which are predominantly excited by PNS convection in 3D. Resonant excitation of high-frequency surface g-modes in 3D by mass motions in the gain layer is suppressed compared to 2D because of smaller downflow velocities and a lack of high-frequency variability in the downflows. In the exploding 20 Msun model, shock revival results in enhanced low-frequency emission due to a change of the preferred scale of the convective eddies in the PNS convection zone. Estimates of the expected excess power in two frequency bands suggests that second-generation detectors will only be able to detect very nearby events, but that third-generation detectors could distinguish SASI- and convection-dominated models at distances of ~10 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 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.

An overview of grid-based numerical methods used in relativistic hydrodynamics (RHD) and magnetohydrodynamics (RMHD) is presented. Special emphasis is put on a comprehensive review of the application of high-resolution shock-capturing methods. Results of a set of demanding test bench simulations obtained with different numerical methods are compared in an attempt to assess the present capabilities and limits of the various numerical strategies. Applications to three astrophysical phenomena are briefly discussed to motivate the need for and to demonstrate the success of RHD and RMHD simulations in their understanding. The review further provides FORTRAN programs to compute the exact solution of the Riemann problem in RMHD, and to simulate 1D RMHD flows in Cartesian coordinates. Electronic Supplementary Material Supplementary material is available for this article at 10.1007/lrca-2015-3.

We perform a set of neutrino-driven core-collapse supernova (CCSN) simulations studying the hydrodynamical neutron star kick mechanism in three-dimensions. Our simulations produce neutron star (NS) kick velocities in a range between ~100-600 km/s resulting mainly from the anisotropic gravitational tug by the asymmetric mass distribution behind the supernova shock. This stochastic kick mechanism suggests that a NS kick velocity of more than 1000 km/s may as well be possible. An enhanced production of heavy elements in the direction roughly opposite to the NS recoil direction is also observed as a result of the asymmetric explosion. This large scale asymmetry might be detectable and can be used to constrain the NS kick mechanism.

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 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

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., Comment: 24 pages, 12 figures, accepted for publication in The Astrophysical Journal

The Commission home page is maintained by Claus Leitherer and contains general information on the Commission structure and activities, including links to stellar structure resources that were made available by the owners. The resources contain evolutionary tracks and isochrones from various groups, nuclear reaction, EOS, and opacity data as well as links to main astronomical journals. As a routine activity, the Organizing Committee has commented on and ranked proposals for several IAU sponsored meetings. Our Commission acted as one of the coordinating bodies of a Symposium held at the IAU XXVI General Assembly in Prague, August 2006, (IAU Symposium No. 239Convection in Astrophysics, and participated in the organization of the following Joint Discussions: JD05Calibrating the Top of the Stellar Mass-Luminosity Relation, JD06Neutron Stars and Black Holes in Star Clusters, JD08Solar and Stellar Activity Cycles, JD11Pre-Solar Grains as Astrophysical Tools; JD14Modelling Dense Stellar Systems; and JD17Highlights of Recent Progress in the Seismology of the Sun and Sun-like Stars.

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., 21 pages, 14 figures, accepted for publication in MNRAS

Die besonderen physikalischen und chemischen Eigenschaften des Wassers haben bei der Entstehung des Lebens zweifellos eine entscheidende Rolle gespielt. Auch heute noch ist Wasser als Losungsmittel eine unverzichtbare Voraussetzung fur alle biochemischen und physiologischen Ablaufe („corpora non agunt nisi soluta“). Eine ahnlich entscheidende Rolle spielt die Temperatur: Mit sinkender Temperatur verringert sich normalerweise die Geschwindigkeit (bio-)chemischer Reaktionen; und sinkt die Umgebungstemperatur bei normalen Druckverhaltnissen (∽ 1 bar) unter den Schmelzpunkt von Wasser (0 °C), geht dieses in seinen festen Aggregatzustand (Eis) uber. Dann ist kein Stofftransport mehr moglich und die Eiskristalle konnen uberlebenswichtige Strukturen, z. B. Membranen, irreparabel beschadigen. Der Temperaturbereich, innerhalb dessen Lebensvorgange moglich sind, erstreckt sich daher in der Regel von etwa 0–45 °C. Lebensraume mit niedrigeren bzw. hoheren Temperaturen konnen nur von Organismen erfolgreich besiedelt werden, die Wege gefunden haben, die todlichen Gefahren auserhalb dieses Temperaturbereichs zu vermeiden.

We have recently presented a new approach for numerical relativity simulations in spherical polar coordinates, both for vacuum and for relativistic hydrodynamics. Our approach is based on a reference-metric formulation of the BSSN equations, a factoring of all tensor components, as well as a partially implicit Runge-Kutta method, and does not rely on a regularization of the equations, nor does it make any assumptions about the symmetry across the origin. In order to demonstrate this feature we present here several off-centered simulations, including simulations of single black holes and neutron stars whose center is placed away from the origin of the coordinate system, as well as the asymmetric head-on collision of two black holes. We also revisit our implementation of relativistic hydrodynamics and demonstrate that a reference-metric formulation of hydrodynamics together with a factoring of all tensor components avoids problems related to the coordinate singularities at the origin and on the axes. As a particularly demanding test we present results for a shock wave propagating through the origin of the spherical polar coordinate system., Comment: 13 pages, 11 figures; matches version published in PRD

Simulation is an indispensable and expedient tool for understanding many astrophysical phenomena. Addressing a few challenging computational problems in astrophysics, it is illustrated that simulation can provide the link between observations and theoretical models of astrophysical phenomenon. For obvious reasons, the discussion is restricted to a few specific timely problems, the main bias being the personal involvement of the author into the simulation of these astrophysics problems. In particular, recent work on core collapse supernovae and their gravitational wave signal, on extragalactic jets, and on the modeling of gamma-ray burst sources is presented.

In situ polymerase chain reaction (isPCR) has been applied in many fields that require detection of a genomic marker in combination with its topographic localization in tissue. We describe here a novel approach that circumvents the major drawbacks of in situ PCR, ie, low sensitivity, leakage of DNA from cells, and inability to quantify the DNA input. Frozen sections of a lymph node from a human immunodeficiency virus (HIV)-1-infected patient were fixed on glass microscope slides, and the glass was scored into square fragments of 0.5-mm edge length using a diamond cutting device. Slides were then attached to adhesive, elastic plastic foil and finally broken, and the foil was extended to allow sorting of fragments into PCR microtiter plates. The material was tested for HIV-1 proviral DNA by a sensitive real-time PCR protocol. Subjacent sections were stained for follicular dendritic cells to identify follicles. The fragmentation process prevented leakage of amplified DNA to neighboring areas as often experienced with in situ PCR. Provirus was clearly associated with follicular areas, in which provirus-carrying cells represented an average of 0.8% of the total cell population (peak density, 3.1% of all follicular cells). The results of this method suggest that the high density of provirus-containing cells in follicles may be important for the persistence of proviral DNA in infected persons.

SummaryWe investigate explosive nuclear burning in core collapse supernovae by coupling a tracer particle method to one- and two-dimensional Eulerian hydrodynamic calculations. Adopting the most recent experimental and theoretical nuclear data, we compute the nucleosynthetic yields for 15 M⊙ stars with solar metallicity, by post-processing the temperature and density history of advected tracer particles. We compare our results to 1D calculations published in the literature.

We present the results of a two-dimensional relativistic hydrodynamic simulation of collisions of dense shells of matter moving within a uniform jet. The non-thermal synchrotron radiation produced by the relativistic electrons injected at shocks is computed following their temporal and spatial evolution. We test different parameterizations of the shock acceleration process and compute the corresponding X-ray light curves. A time lag between hard and soft X-ray radiation is found. The collision has an efficiency of few times 10-3 in converting kinetic energy into radiation.

In a series of time dependent numerical simulations we have performed a parameter study of magnetised relativistic jets. We have found that the impact of the magnetic field on the morphology of a jet depends strongly on the configuration of the field.

Trees inhabiting central Amazon floodplain forests are subjected to an annual flood-pulse lasting up to 10 months, leading to both oxygen shortage and accumulation of high levels of reduced iron. To understand the mechanisms underlying the adaptation to these conditions, cuttings from three tree species typical of varzea inundation forests (Salix martiana, Tabernaemontana juruana, and Laetia corymbulosa), were cultivated either aerobically or anaerobically under different iron regimes in greenhouse experiments. Although all species are considered to be non-deciduous, Laetia corymbulosa lost and formed new leaves continuously during the experimental period. Although relative growth rates (RGRs) of all species declined in response to hypoxic conditions, no marked changes in RGRs were apparent among different iron concentrations in the growth medium, ranging from 50 to 500 µM, supplied in ferrous form as FeSO4. Whereas roots exhibited color changes due to the formation of iron precipitates, no visual symptoms of iron toxicity were observed in the leaves. Iron concentration increased in all organs of all species with increasing iron concentrations in the medium, except for leaves of S. martiana and T. juruana, suggesting an effective restriction of iron influx into the leaf symplast. Although the leaf iron concentration was at the upper limit of the critical range at high external iron levels, it is suggested that internal active transport rather than intracellular detoxification mechanisms contribute to the tolerance to supra-optimal iron levels. Anatomical traits such as suberization of peripheral cell walls and the formation of aerenchyma appear to be of minor importance for Fe tolerance.

The formation of suberized and lignified barriers in the exodermis is suggested to be part of a suite of adaptations to flooded or waterlogged conditions, adjusting transport of solutes and gases in and out of roots. In this study, the composition of apoplasmic barriers in hypodermal cell walls and oxygen profiles in roots and the surrounding medium of four Amazon tree species that are subjected to long-term flooding at their habitat was analyzed. In hypodermal cell walls of the deciduous tree Crateva benthami, suberization is very weak and dominated by monoacids, 2-hydroxy acids, and ω-hydroxycarboxylic acids. This species does not show any morphological adaptations to flooding and overcomes the aquatic period in a dormant state. Hypodermal cells of Tabernaemontana juruana, a tree which is able to maintain its leaf system during the aquatic phase, are characterized by extensively suberized walls, incrusted mainly by the unsaturated C18ω-hydroxycarboxylic acid and the α,ω-dicarboxylic acid analogon, known as typical suberin markers. Two other evergreen species,Laetia corymbulosa and Salix martiana, contained 3- to 4-fold less aliphatic suberin in the exodermis, but more than 85% of the aromatic moiety of suberin are composed ofpara-hydroxybenzoic acid, suggesting a function of suberin in pathogen defense. No major differences in the lignin content among the species were observed. Determination of oxygen distribution in the roots and rhizosphere of the four species revealed that radial loss of oxygen can be effectively restricted by the formation of suberized barriers but not by lignification of exodermal cell walls.

Objective: Due to the heterogeneous biologic behavior of stage T1 bladder carcinomas, there is a need for new markers allowing to assess the prognosis more accurately. To our knowledge, there are no reports on studies investigating minichromosome maintenance protein 2 (MCM2) expression in bladder carcinomas. Thus, we investigated the prognostic value of MCM2 immunoreactivity in stage T1 bladder tumors. Methods: Fifty-four tumors were analyzed using Biochip microarrays. Also p53 and Ki67 antigen expression were examined. Immunohistochemical scores were compared with the clinical outcome. Results: During a median follow-up of 43 months, tumor recurrence was registered in 43 and progression to stage T2 in 19 patients. Kaplan–Meier curves demonstrated that high-level MCM2 expression was significantly associated with early tumor recurrence when using a cutoff of 60% ( p =0.0035 by log-rank test), and with early tumor progression when using a cutoff of 20% ( p =0.0454). There was no relationship ( p =0.604) between MCM2 and p53, but a tendentious relationship ( p =0.082) between MCM2 and Ki67 antigen expression. MCM2 ( p =0.006), Ki67 antigen ( p =0.035) and p53 expression ( p =0.049) as well as tumor grade ( p =0.026) and age ( p =0.025) were found significantly associated with recurrence-free survival by univariate Cox regression analysis, among which only Ki67 antigen expression ( p =0.015) and age ( p =0.019) proved to be of independent predictive value by multivariate analysis. Concerning tumor progression, MCM2 expression was identified as the only predictive parameter by log-rank test, but it was not of independent predictive value by multivariate analysis ( p =0.101). Conclusion: Our data suggest that MCM2 expression may bear some prognostic relevance in stage T1 bladder carcinomas.

Varzeas are species-rich forest communities of the Central Amazon floodplains, inhabited by highly adapted tree species that can withstand long flooding periods. The leaf shedding behaviour and morphological traits that may contribute to adaptation to low oxygen levels were studied at the Ilha de Marchantaria on the lower Solimoes-Amazonas river, Brazil, and in greenhouse experiments with cuttings of six tree species typical of the Amazon floodplain. Comparison of deciduousness in situ revealed that four of the species under investigation, Salix martiana, Tabernaemontana juruana, Laetia corymbulosa and Pouteria glomerata, are able to maintain their leaf system during the aquatic period. Adventitious roots were formed by S. martiana and T. juruana, but anatomical characteristics differed between the species. Whereas S. martiana developed lysigenous aerenchyma in its roots during aerobic and anaerobic growth, only small intercellular spaces of schizogenous origin were formed in the root cortex of T. juruana. Similar to the latter species, such spaces were constitutively formed in the deciduous species Crateva benthamiand Vitex cymosa. Suberin deposits were observed in tangential and radial cell walls of the hypodermis of roots from T. juruana, L. corymbulosa and P. glomerata. Suberin deposits were less pronounced in roots of S. martiana and absent in V. cymosaand C. benthami. The data show that different, almost contrasting, survival mechanisms have evolved in roots of plants with similar life forms in the same habitat. The results further suggest that the morphological traits of the investigated trees are causally linked with the in situ leaf shedding behaviour.

We present results from simulations of axisymmetric relativistic rotational core collapse. The main objective of our investigation is to compute the waveforms of gravitational radiation emitted in such events, extending previous Newtonian simulations to relativity. The general relativistic hydrodynamic equations are formulated in flux-conservative form and solved using a high-resolution shock-capturing scheme. The Einstein equations are solved assuming a conformally flat 3-metric and the quadrupole formula is used to extract waveforms of the gravitational radiation emitted during the collapse. A comparison of our results with those of Newtonian simulations shows that gravitational wave amplitudes agree within 30%. Surprisingly, in some cases, relativistic effects actually diminish the amplitude of the gravitational wave signal. We further find that the parameter range of models suffering multiple coherent bounces due to centrifugal forces is considerably smaller than in Newtonian simulations.

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., Comment: 36 pages, 35 figures, accepted for publication in ApJS

We have generalized a covariant and conformal version of the Z4 system of the Einstein equations by adopting a reference metric approach, that we denote as fCCZ4, well suited for curvilinear as well as Cartesian coordinates. We implement this formalism in spherical polar coordinates under the assumption of spherical symmetry using a partially-implicit Runge-Kutta (PIRK) method, without using any regularization scheme, and show that our code can evolve both vacuum and non-vacuum spacetimes without encountering instabilities. We have performed several tests and compared the Hamiltonian constraint violations of the fCCZ4 system, for different choices of certain free parameters, with these of BSSN. For an optimal choice of these parameters, and for neutron star spacetimes, the violations of the Hamiltonian constraint can be between 1 and 3 orders of magnitude smaller in the fCCZ4 system than in the BSSN formulation. For black hole spacetimes, on the other hand, any advantages of fCCZ4 over BSSN are less evident.

We present 3D simulations of core-collapse supernovae from blast-wave initiation by the neutrino-driven mechanism to shock breakout from the stellar surface, considering two 15 Msun red supergiants (RSG) and two blue supergiants (BSG) of 15 Msun and 20 Msun. We demonstrate that the metal-rich ejecta in homologous expansion still carry fingerprints of asymmetries at the beginning of the explosion, but the final metal distribution is massively affected by the detailed progenitor structure. The most extended and fastest metal fingers and clumps are correlated with the biggest and fastest-rising plumes of neutrino-heated matter, because these plumes most effectively seed the growth of Rayleigh-Taylor (RT) instabilities at the C+O/He and He/H composition-shell interfaces after the passage of the SN shock. The extent of radial mixing, global asymmetry of the metal-rich ejecta, RT-induced fragmentation of initial plumes to smaller-scale fingers, and maximal Ni and minimal H velocities do not only depend on the initial asphericity and explosion energy (which determine the shock and initial Ni velocities) but also on the density profiles and widths of C+O core and He shell and on the density gradient at the He/H transition, which lead to unsteady shock propagation and the formation of reverse shocks. Both RSG explosions retain a great global metal asymmetry with pronounced clumpiness and substructure, deep penetration of Ni fingers into the H-envelope (with maximum velocities of 4000-5000 km/s for an explosion energy around 1.5 bethe) and efficient inward H-mixing. While the 15 Msun BSG shares these properties (maximum Ni speeds up to ~3500 km/s), the 20 Msun BSG develops a much more roundish geometry without pronounced metal fingers (maximum Ni velocities only ~2200 km/s) because of reverse-shock deceleration and insufficient time for strong RT growth and fragmentation at the He/H interface., 21 pages, 15 figures; revised version with minor changes in Sect.1; accepted by Astron. Astrophys

We adopt a reference-metric approach to generalize a covariant and conformal version of the Z4 system of the Einstein equations. We refer to the resulting system as ``fully covariant and conformal", or fCCZ4 for short, since it is well suited for curvilinear as well as Cartesian coordinates. We implement this fCCZ4 formalism in spherical polar coordinates under the assumption of spherical symmetry using a partially-implicit Runge-Kutta (PIRK) method and show that our code can evolve both vacuum and non-vacuum spacetimes without encountering instabilities. Our method does not require regularization of the equations to handle coordinate singularities, nor does it depend on constraint-preserving outer boundary conditions, nor does it need any modifications of the equations for evolutions of black holes. We perform several tests and compare the performance of the fCCZ4 system, for different choices of certain free parameters, with that of BSSN. Confirming earlier results we find that, for an optimal choice of these parameters, and for neutron-star spacetimes, the violations of the Hamiltonian constraint can be between 1 and 3 orders of magnitude smaller in the fCCZ4 system than in the BSSN formulation. For black-hole spacetimes, on the other hand, any advantages of fCCZ4 over BSSN are less evident., 13 pages, 10 figures

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., Comment: 9 pages, 5 figures, 2 tables, Accepted by MNRAS

The magnetorotational instability (MRI) is a promising mechanism to amplify the magnetic field in fast rotating protoneutron stars. The diffusion of neutrinos trapped in the PNS induces a transport of momentum, which can be modelled as a viscosity on length-scales longer than the neutrino mean free path. This neutrino-viscosity can slow down the growth of MRI modes to such an extent that a minimum initial magnetic field strength of $\gtrsim 10^{12}\, {\rm G}$ is needed for the MRI to grow on a sufficiently short time-scale to potentially affect the explosion. It is uncertain whether the magnetic field of fast rotating progenitor cores is strong enough to yield such an initial magnetic field in PNS. At MRI wavelengths shorter than the neutrino mean free path, on the other hand, neutrino radiation does not act as a viscosity but rather induces a drag on the velocity with a damping rate independent of the wavelength. We perform a linear analysis of the MRI in this regime, and apply our analytical results to the PNS structure from a one-dimensional numerical simulation. We show that in the outer layers of the PNS, the MRI can grow from weak magnetic fields at wavelengths shorter than the neutrino mean free path, while deeper in the PNS MRI growth takes place in the viscous regime and requires a minimum magnetic field strength., Comment: Accepted for publication in MNRAS, one figure added (Fig. 11), 14 pages, 11 figures

The development of brown trout (Salmo trutta f. fario L.) in water of two differently polluted streams and in a control situation was monitored in order to get insights into the impact of anthropogenic chemical stressors on the reproductive success of this fish species indigenous to both streams. The test streams, situated in the south of Stuttgart, Germany, were the complexly polluted Körsch stream and the less polluted Krähenbach stream. Bypass systems connected to the streams and a laboratory control system were used for continuous exposure of early brown trout stages shortly after fertilisation up to the end of the embryonic development. Temperature and oxygen conditions were standardised in all test series in order to minimise unspecific effects. The examined endpoints were: (1) mortality, (2) developmental rate, (3) time course of hatching, (4) malformations, and (5) growth. A retarded development, reduced growth rates and higher mortality rates of Körsch stream water exposed embryos indicated an embryotoxic potential for the more polluted stream. High infection-related mortality rates of embryos suggested the presence of confounding factors also in the less polluted Krähenbach stream. In parallel to the exposure experiment, physicochemical and limnochemical parameters as well as concentrations of organic contaminants and heavy metals were monitored. Analytical data confirm the different degrees of pollution of both streams.