An Advanced Leakage Scheme for Neutrino Treatment in Astrophysical Simulations

We present an Advanced Spectral Leakage (ASL) scheme to model neutrinos in the context of core-collapse supernovae (CCSNe) and compact binary mergers. Based on previous gray leakage schemes, the ASL scheme computes the neutrino cooling rates by interpolating local production and diffusion rates (relevant in optically thin and thick regimes, respectively) separately for discretized values of the neutrino energy. Neutrino trapped components are also modeled, based on equilibrium and timescale arguments. The better accuracy achieved by the spectral treatment allows a more reliable computation of neutrino heating rates in optically thin conditions. The scheme has been calibrated and tested against Boltzmann transport in the context of Newtonian spherically symmetric models of CCSNe. ASL shows a very good qualitative and a partial quantitative agreement for key quantities from collapse to a few hundreds of milliseconds after core bounce. We have proved the adaptability and flexibility of our ASL scheme, coupling it to an axisymmetric Eulerian and to a three-dimensional smoothed particle hydrodynamics code to simulate core collapse. Therefore, the neutrino treatment presented here is ideal for large parameter-space explorations, parametric studies, high-resolution tests, code developments, and long-term modeling of asymmetric configurations, where more detailed neutrino treatments are not available or are currently computationally too expensive

Shiga toxin-producing Escherichia coli O157 associated with human infections in Switzerland, 2000-2009

Shiga toxin-producing Escherichia coli (STEC), an important foodborne pathogen, can cause mild to severe bloody diarrhoea (BD), sometimes followed by life-threatening complications such as haemolytic uraemic syndrome (HUS). A total of 44 O157 strains isolated from different patients from 2000 through 2009 in Switzerland were further characterized and linked to medical history data. Non-bloody diarrhoea was experienced by 15·9%, BD by 61·4% of the patients, and 29·5% developed HUS. All strains belonged to MLST type 11, were positive for stx2 variants (stx2 and/or stx2c), eae and ehxA, and only two strains showed antibiotic resistance. Of the 44 strains, nine phage types (PTs) were detected the most frequent being PT32 (43·2%) and PT8 (18·2%). By PFGE, 39 different patterns were found. This high genetic diversity within the strains leads to the conclusion that STEC O157 infections in Switzerland most often occur as sporadic case

A New Lower Bound on the Maximum Number of Satisfied Clauses in Max-SAT and its Algorithmic Applications

A pair of unit clauses is called conflicting if it is of the form $(x)$, $(\bar{x})$. A CNF formula is unit-conflict free (UCF) if it contains no pair of conflicting unit clauses. Lieberherr and Specker (J. ACM 28, 1981) showed that for each UCF CNF formula with $m$ clauses we can simultaneously satisfy at least \pp m clauses, where \pp =(\sqrt{5}-1)/2. We improve the Lieberherr-Specker bound by showing that for each UCF CNF formula $F$ with $m$ clauses we can find, in polynomial time, a subformula $F'$ with $m'$ clauses such that we can simultaneously satisfy at least \pp m+(1-\pp)m'+(2-3\pp)n"/2 clauses (in $F$), where $n"$ is the number of variables in $F$ which are not in $F'$. We consider two parameterized versions of MAX-SAT, where the parameter is the number of satisfied clauses above the bounds $m/2$ and $m(\sqrt{5}-1)/2$. The former bound is tight for general formulas, and the later is tight for UCF formulas. Mahajan and Raman (J. Algorithms 31, 1999) showed that every instance of the first parameterized problem can be transformed, in polynomial time, into an equivalent one with at most $6k+3$ variables and $10k$ clauses. We improve this to $4k$ variables and $(2\sqrt{5}+4)k$ clauses. Mahajan and Raman conjectured that the second parameterized problem is fixed-parameter tractable (FPT). We show that the problem is indeed FPT by describing a polynomial-time algorithm that transforms any problem instance into an equivalent one with at most $(7+3\sqrt{5})k$ variables. Our results are obtained using our improvement of the Lieberherr-Specker bound above

A well-balanced scheme for the simulation tool-kit A-MaZe: implementation, tests, and first applications to stellar structure

This is the final version. Available on open access from EDP Sciences via the DOI in this recordCharacterizing stellar convection in multiple dimensions is a topic at the forefront of stellar astrophysics. Numerical simulations are an essential tool for this task. We present an extension of the existing numerical tool-kit A-MaZe that enables such simulations of stratified flows in a gravitational field. The finite-volume based, cell-centered, and time-explicit hydrodynamics solver of A-MaZe was extended such that the scheme is now well-balanced in both momentum and energy. The algorithm maintains an initially static balance between gravity and pressure to machine precision. Quasi-stationary convection in slab-geometry preserves gas energy (internal plus kinetic) on average despite strong local up- and down-drafts. By contrast, a more standard numerical scheme is demonstrated to result in substantial gains of energy within a short time on purely numerical grounds. The test is further used to point out the role of dimensionality, viscosity, and Rayleigh number for compressible convection. Applications to a young sun in 2D and 3D, covering a part of the inner radiative zone as well as the outer convective zone, demonstrate that the scheme meets its initial design goal. Comparison with results obtained for a physically identical setup with a time-implicit code show qualitative agreement.European Research Council (ERC)French National Program for High Energies PNH

Black hole partition functions and duality

The macroscopic entropy and the attractor equations for BPS black holes in four-dimensional N=2 supergravity theories follow from a variational principle for a certain `entropy function'. We present this function in the presence of R^2-interactions and non-holomorphic corrections. The variational principle identifies the entropy as a Legendre transform and this motivates the definition of various partition functions corresponding to different ensembles and a hierarchy of corresponding duality invariant inverse Laplace integral representations for the microscopic degeneracies. Whenever the microscopic degeneracies are known the partition functions can be evaluated directly. This is the case for N=4 heterotic CHL black holes, where we demonstrate that the partition functions are consistent with the results obtained on the macroscopic side for black holes that have a non-vanishing classical area. In this way we confirm the presence of a measure in the duality invariant inverse Laplace integrals. Most, but not all, of these results are obtained in the context of semiclassical approximations. For black holes whose area vanishes classically, there remain discrepancies at the semiclassical level and beyond, the nature of which is not fully understood at present.Comment: 36 pages, Late

The Impact of Fission on R-Process Calculations

We have performed r-process calculations in neutron star mergers (NSM) and jets of magnetohydrodynamically driven (MHD) supernovae. In these very neutron-rich environments the fission model of heavy nuclei has an impact on the shape of the final abundance distribution and the second r-process peak in particular. We have studied the effect of different fission fragment mass distribution models in calculations of low-Ye ejecta, ranging from a simple parametrization to extensive statistical treatments (ABLA07). The r-process path ends when it reaches an area in the nuclear chart where fission dominates over further neutron captures. The position of this point is determined by the fission barriers and the neutron separation energies of the nuclei involved. As these values both depend on the choice of the nuclear mass model, so does the r-process path. Here we present calculations using the FRDM (Finite Range Droplet Model) and the ETFSI (Extended Thomas Fermi with Strutinsky Integral) mass model with the related TF and ETFSI fission barrier predictions. Utilizing sophisticated fission fragment distribution leads to a highly improved abundance distribution.Peer reviewe

Core-collapse supernova explosions triggered by a quark-hadron phase transition during the early post-bounce phase

We explore explosions of massive stars, which are triggered via the quark-hadron phase transition during the early post bounce phase of core-collapse supernovae. We construct a quark equation of state, based on the bag model for strange quark matter. The transition between the hadronic and the quark phases is constructed applying Gibbs conditions. The resulting quark-hadron hybrid equations of state are used in core-collapse supernova simulations, based on general relativistic radiation hydrodynamics and three flavor Boltzmann neutrino transport in spherical symmetry. The formation of a mixed phase reduces the adiabatic index, which induces the gravitational collapse of the central protoneutron star. The collapse halts in the pure quark phase, where the adiabatic index increases. A strong accretion shock forms, which propagates towards the protoneutron star surface. Due to the density decrease of several orders of magnitude, the accretion shock turns into a dynamic shock with matter outflow. This moment defines the onset of the explosion in supernova models that allow for a quark-hadron phase transition, where otherwise no explosions could be obtained. The shock propagation across the neutrinospheres releases a burst of neutrinos. This serves as a strong observable identification for the structural reconfiguration of the stellar core. The ejected matter expands on a short timescale and remains neutron-rich. These conditions might be suitable for the production of heavy elements via the r-process. The neutron-rich material is followed by proton-rich neutrino-driven ejecta in the later cooling phase of the protoneutron star where the vp-process might occur.Comment: 29 pages, 24 figures, submitted to Ap

Gravitational waves from supernova matter

We have performed a set of 11 three-dimensional magnetohydrodynamical core collapse supernova simulations in order to investigate the dependencies of the gravitational wave signal on the progenitor's initial conditions. We study the effects of the initial central angular velocity and different variants of neutrino transport. Our models are started up from a 15 solar mass progenitor and incorporate an effective general relativistic gravitational potential and a finite temperature nuclear equation of state. Furthermore, the electron flavour neutrino transport is tracked by efficient algorithms for the radiative transfer of massless fermions. We find that non- and slowly rotating models show gravitational wave emission due to prompt- and lepton driven convection that reveals details about the hydrodynamical state of the fluid inside the protoneutron stars. Furthermore we show that protoneutron stars can become dynamically unstable to rotational instabilities at T/|W| values as low as ~2 % at core bounce. We point out that the inclusion of deleptonization during the postbounce phase is very important for the quantitative GW prediction, as it enhances the absolute values of the gravitational wave trains up to a factor of ten with respect to a lepton-conserving treatment.Comment: 10 pages, 6 figures, accepted, to be published in a Classical and Quantum Gravity special issue for MICRA200

Like Will to Like: Abundances of Closely Related Species Can Predict Susceptibility to Intestinal Colonization by Pathogenic and Commensal Bacteria

The intestinal ecosystem is formed by a complex, yet highly characteristic microbial community. The parameters defining whether this community permits invasion of a new bacterial species are unclear. In particular, inhibition of enteropathogen infection by the gut microbiota ( = colonization resistance) is poorly understood. To analyze the mechanisms of microbiota-mediated protection from Salmonella enterica induced enterocolitis, we used a mouse infection model and large scale high-throughput pyrosequencing. In contrast to conventional mice (CON), mice with a gut microbiota of low complexity (LCM) were highly susceptible to S. enterica induced colonization and enterocolitis. Colonization resistance was partially restored in LCM-animals by co-housing with conventional mice for 21 days (LCMcon21). 16S rRNA sequence analysis comparing LCM, LCMcon21 and CON gut microbiota revealed that gut microbiota complexity increased upon conventionalization and correlated with increased resistance to S. enterica infection. Comparative microbiota analysis of mice with varying degrees of colonization resistance allowed us to identify intestinal ecosystem characteristics associated with susceptibility to S. enterica infection. Moreover, this system enabled us to gain further insights into the general principles of gut ecosystem invasion by non-pathogenic, commensal bacteria. Mice harboring high commensal E. coli densities were more susceptible to S. enterica induced gut inflammation. Similarly, mice with high titers of Lactobacilli were more efficiently colonized by a commensal Lactobacillus reuteri RR strain after oral inoculation. Upon examination of 16S rRNA sequence data from 9 CON mice we found that closely related phylotypes generally display significantly correlated abundances (co-occurrence), more so than distantly related phylotypes. Thus, in essence, the presence of closely related species can increase the chance of invasion of newly incoming species into the gut ecosystem. We provide evidence that this principle might be of general validity for invasion of bacteria in preformed gut ecosystems. This might be of relevance for human enteropathogen infections as well as therapeutic use of probiotic commensal bacteria