Surface detonation in type Ia supernova explosions?

We explore the evolution of thermonuclear supernova explosions when the progenitor white dwarf star ignites asymmetrically off-center. Several numerical simulations are carried out in two and three dimensions to test the consequences of different initial flame configurations such as spherical bubbles displaced from the center, more complex deformed configurations, and teardrop-shaped ignitions. The burning bubbles float towards the surface while releasing energy due to the nuclear reactions. If the energy release is too small to gravitationally unbind the star, the ash sweeps around it, once the burning bubble approaches the surface. Collisions in the fuel on the opposite side increase its temperature and density and may -- in some cases -- initiate a detonation wave which will then propagate inward burning the core of the star and leading to a strong explosion. However, for initial setups in two dimensions that seem realistic from pre-ignition evolution, as well as for all three-dimensional simulations the collimation of the surface material is found to be too weak to trigger a detonation.Comment: 5 pages, 3 figures, in: Proceedings of the SciDAC 2006 Meeting, Denver June 25-26 2006, also available at http://herald.iop.org/jpcs46/m51/gbr//link/40

Spin 0 and spin 1/2 particles in a constant scalar-curvature background

We study the Klein-Gordon and Dirac equations in the presence of a background metric ds^2 = -dt^2 + dx^2 + e^{-2gx}(dy^2 + dz^2) in a semi-infinite lab (x>0). This metric has a constant scalar curvature R=6g^2 and is produced by a perfect fluid with equation of state p=-\rho /3. The eigenfunctions of spin-0 and spin-1/2 particles are obtained exactly, and the quantized energy eigenvalues are compared. It is shown that both of these particles must have nonzero transverse momentum in this background. We show that there is a minimum energy E^2_{min}=m^2c^4 + g^2c^2\hbar^2$ for bosons E_{KG} > E_{min}, while the fermions have no specific ground state E_{Dirac}>mc^2.Comment: LateX, 10 page

Optimisation of patch distribution strategies for AMR applications

As core counts increase in the world's most powerful supercomputers, applications are becoming limited not only by computational power, but also by data availability. In the race to exascale, efficient and effective communication policies are key to achieving optimal application performance. Applications using adaptive mesh refinement (AMR) trade off communication for computational load balancing, to enable the focused computation of specific areas of interest. This class of application is particularly susceptible to the communication performance of the underlying architectures, and are inherently difficult to scale efficiently. In this paper we present a study of the effect of patch distribution strategies on the scalability of an AMR code. We demonstrate the significance of patch placement on communication overheads, and by balancing the computation and communication costs of patches, we develop a scheme to optimise performance of a specific, industry-strength, benchmark application

Emergent Mesoscale Phenomena in Magnetized Accretion Disc Turbulence

We study how the structure and variability of magnetohydrodynamic (MHD) turbulence in accretion discs converge with domain size. Our results are based on a series of vertically stratified local simulations, computed using the Athena code, that have fixed spatial resolution, but varying radial and azimuthal extent (from \Delta R = 0.5H to 16H, where H is the vertical scale height). We show that elementary local diagnostics of the turbulence, including the Shakura-Sunyaev {\alpha} parameter, the ratio of Maxwell stress to magnetic energy, and the ratio of magnetic to fluid stresses, converge to within the precision of our measurements for spatial domains of radial size Lx \geq 2H. We obtain {\alpha} = 0.02-0.03, consistent with recent results. Very small domains (Lx = 0.5H) return anomalous results, independent of spatial resolution. The convergence with domain size is only valid for a limited set of diagnostics: larger spatial domains admit the emergence of dynamically important mesoscale structures. In our largest simulations, the Maxwell stress shows a significant large scale non-local component, while the density develops long-lived axisymmetric perturbations (zonal flows) at the 20% level. Most strikingly, the variability of the disc in fixed-sized patches decreases strongly as the simulation volume increases. We find generally good agreement between our largest local simulations and global simulations with comparable spatial resolution. There is no direct evidence that the presence of curvature terms or radial gradients in global calculations materially affect the turbulence, except to perhaps introduce an outer radial scale for mesoscale structures. The demonstrated importance of mean magnetic fields, seen in both large local and global simulations implies that the growth and saturation of these fields is likely of critical importance for the evolution of accretion discs. (abridged)Comment: 18 pages, 20 figures, accepted to MNRA

Simulating Turbulence Using the Astrophysical Discontinuous Galerkin Code TENET

In astrophysics, the two main methods traditionally in use for solving the Euler equations of ideal fluid dynamics are smoothed particle hydrodynamics and finite volume discretization on a stationary mesh. However, the goal to efficiently make use of future exascale machines with their ever higher degree of parallel concurrency motivates the search for more efficient and more accurate techniques for computing hydrodynamics. Discontinuous Galerkin (DG) methods represent a promising class of methods in this regard, as they can be straightforwardly extended to arbitrarily high order while requiring only small stencils. Especially for applications involving comparatively smooth problems, higher-order approaches promise significant gains in computational speed for reaching a desired target accuracy. Here, we introduce our new astrophysical DG code TENET designed for applications in cosmology, and discuss our first results for 3D simulations of subsonic turbulence. We show that our new DG implementation provides accurate results for subsonic turbulence, at considerably reduced computational cost compared with traditional finite volume methods. In particular, we find that DG needs about 1.8 times fewer degrees of freedom to achieve the same accuracy and at the same time is more than 1.5 times faster, confirming its substantial promise for astrophysical applications.Comment: 21 pages, 7 figures, to appear in Proceedings of the SPPEXA symposium, Lecture Notes in Computational Science and Engineering (LNCSE), Springe

An Adaptive Semi-Implicit Scheme for Simulations of Unsteady Viscous Compressible Flows

A numerical scheme for simulation of unsteady, viscous, compressible flows is considered. The scheme employs an explicit discretization of the inviscid terms of the Navier-Stokes equations and an implicit discretization of the viscous terms. The discretization is second order accurate in both space and time. Under appropriate assumptions, the implicit system of equations can be decoupled into two linear systems of reduced rank. These are solved efficiently using a Gauss-Seidel method with multigrid convergence acceleration. When coupled with a solution-adaptive mesh refinement technique, the hybrid explicit-implicit scheme provides an effective methodology for accurate simulations of unsteady viscous flows. The methodology is demonstrated for both body-fitted structured grids and for rectangular (Cartesian) grids

Numerical Methods for the Stochastic Landau-Lifshitz Navier-Stokes Equations

The Landau-Lifshitz Navier-Stokes (LLNS) equations incorporate thermal fluctuations into macroscopic hydrodynamics by using stochastic fluxes. This paper examines explicit Eulerian discretizations of the full LLNS equations. Several CFD approaches are considered (including MacCormack's two-step Lax-Wendroff scheme and the Piecewise Parabolic Method) and are found to give good results (about 10% error) for the variances of momentum and energy fluctuations. However, neither of these schemes accurately reproduces the density fluctuations. We introduce a conservative centered scheme with a third-order Runge-Kutta temporal integrator that does accurately produce density fluctuations. A variety of numerical tests, including the random walk of a standing shock wave, are considered and results from the stochastic LLNS PDE solver are compared with theory, when available, and with molecular simulations using a Direct Simulation Monte Carlo (DSMC) algorithm

Inhibition of Bromodomain and Extraterminal Domain (BET) Proteins by JQ1 Unravels a Novel Epigenetic Modulation to Control Lipid Homeostasis

The homeostatic control of lipid metabolism is essential for many fundamental physiological processes. A deep understanding of its regulatory mechanisms is pivotal to unravel prospective physiopathological factors and to identify novel molecular targets that could be employed to design promising therapies in the management of lipid disorders. Here, we investigated the role of bromodomain and extraterminal domain (BET) proteins in the regulation of lipid metabolism. To reach this aim, we used a loss-of-function approach by treating HepG2 cells with JQ1, a powerful and selective BET inhibitor. The main results demonstrated that BET inhibition by JQ1 efficiently decreases intracellular lipid content, determining a significant modulation of proteins involved in lipid biosynthesis, uptake and intracellular trafficking. Importantly, the capability of BET inhibition to slow down cell proliferation is dependent on the modulation of cholesterol metabolism. Taken together, these data highlight a novel epigenetic mechanism involved in the regulation of lipid homeostasis

Directional Charge Separation in Isolated Organic Semiconductor Crystalline Nanowires

One of the fundamental design paradigms in organic photovoltaic device engineering is based on the idea that charge separation is an extrinsically driven process requiring an interface for exciton fission. This idea has driven an enormous materials science engineering effort focused on construction of domain sizes commensurate with a nominal exciton diffusion length of order 10 nm. Here, we show that polarized optical excitation of isolated pristine crystalline nanowires of a small molecule n-type organic semiconductor, 7,8,15,16-tetraazaterrylene, generates a significant population of charge-separated polaron pairs along the π-stacking direction. Charge separation was signalled by pronounced power-law photoluminescence decay polarized along the same axis. In the transverse direction, we observed exponential decay associated with excitons localized on individual monomers. We propose that this effect derives from an intrinsic directional charge-transfer interaction that can ultimately be programmed by molecular packing geometry