Chaste: an open source C++ library for computational physiology and biology

Chaste - Cancer, Heart And Soft Tissue Environment - is an open source C++ library for the computational simulation of mathematical models developed for physiology and biology. Code development has been driven by two initial applications: cardiac electrophysiology and cancer development. A large number of cardiac electrophysiology studies have been enabled and performed, including high performance computational investigations of defibrillation on realistic human cardiac geometries. New models for the initiation and growth of tumours have been developed. In particular, cell-based simulations have provided novel insight into the role of stem cells in the colorectal crypt. Chaste is constantly evolving and is now being applied to a far wider range of problems. The code provides modules for handling common scientific computing components, such as meshes and solvers for ordinary and partial differential equations (ODEs/PDEs). Re-use of these components avoids the need for researchers to "re-invent the wheel" with each new project, accelerating the rate of progress in new applications. Chaste is developed using industrially-derived techniques, in particular test-driven development, to ensure code quality, re-use and reliability. In this article we provide examples that illustrate the types of problems Chaste can be used to solve, which can be run on a desktop computer. We highlight some scientific studies that have used or are using Chaste, and the insights they have provided. The source code, both for specific releases and the development version, is available to download under an open source Berkeley Software Distribution (BSD) licence at http://www.cs.ox.ac.uk/chaste, together with details of a mailing list and links to documentation and tutorials

Anisotropic Structure of the Order Parameter in FeSe0.45Te0.55 Revealed by Angle Resolved Specific Heat

The symmetry and structure of the superconducting gap in the Fe-based superconductors are the central issue for understanding these novel materials. So far the experimental data and theoretical models have been highly controversial. Some experiments favor two or more constant or nearly-constant gaps, others indicate strong anisotropy and yet others suggest gap zeros ("nodes"). Theoretical models also vary, suggesting that the absence or presence of the nodes depends quantitatively on the model parameters. An opinion that has gained substantial currency is that the gap structure, unlike all other known superconductors, including cuprates, may be different in different compounds within the same family. A unique method for addressing this issue, one of the very few methods that are bulk and angle-resolved, calls for measuring the electronic specific heat in a rotating magnetic field, as a function of field orientation with respect to the crystallographic axes. In this Communication we present the first such measurement for an Fe-based high-Tc superconductor (FeBSC). We observed a fourfold oscillation of the specific heat as a function of the in-plane magnetic field direction, which allowed us to identify the locations of the gap minima (or nodes) on the Fermi surface. Our results are consistent with the expectations of an extended s-wave model with a significant gap anisotropy on the electron pockets and the gap minima along the \Gamma M (or Fe-Fe bond) direction.Comment: 32 pages, 7 figure

Gene set of nuclear-encoded mitochondrial regulators is enriched for common inherited variation in obesity

There are hints of an altered mitochondrial function in obesity. Nuclear-encoded genes are relevant for mitochondrial function (3 gene sets of known relevant pathways: (1) 16 nuclear regulators of mitochondrial genes, (2) 91 genes for oxidative phosphorylation and (3) 966 nuclear-encoded mitochondrial genes). Gene set enrichment analysis (GSEA) showed no association with type 2 diabetes mellitus in these gene sets. Here we performed a GSEA for the same gene sets for obesity. Genome wide association study (GWAS) data from a case-control approach on 453 extremely obese children and adolescents and 435 lean adult controls were used for GSEA. For independent confirmation, we analyzed 705 obesity GWAS trios (extremely obese child and both biological parents) and a population-based GWAS sample (KORA F4, n = 1,743). A meta-analysis was performed on all three samples. In each sample, the distribution of significance levels between the respective gene set and those of all genes was compared using the leading-edge-fraction-comparison test (cut-offs between the 50(th) and 95(th) percentile of the set of all gene-wise corrected p-values) as implemented in the MAGENTA software. In the case-control sample, significant enrichment of associations with obesity was observed above the 50(th) percentile for the set of the 16 nuclear regulators of mitochondrial genes (p(GSEA,50) = 0.0103). This finding was not confirmed in the trios (p(GSEA,50) = 0.5991), but in KORA (p(GSEA,50) = 0.0398). The meta-analysis again indicated a trend for enrichment (p(MAGENTA,50) = 0.1052, p(MAGENTA,75) = 0.0251). The GSEA revealed that weak association signals for obesity might be enriched in the gene set of 16 nuclear regulators of mitochondrial genes

The catalytic subunit of the system L1 amino acid transporter (S<i>lc7a5</i>) facilitates nutrient signalling in mouse skeletal muscle

The System L1-type amino acid transporter mediates transport of large neutral amino acids (LNAA) in many mammalian cell-types. LNAA such as leucine are required for full activation of the mTOR-S6K signalling pathway promoting protein synthesis and cell growth. The SLC7A5 (LAT1) catalytic subunit of high-affinity System L1 functions as a glycoprotein-associated heterodimer with the multifunctional protein SLC3A2 (CD98). We generated a floxed Slc7a5 mouse strain which, when crossed with mice expressing Cre driven by a global promoter, produced Slc7a5 heterozygous knockout (Slc7a5+/-) animals with no overt phenotype, although homozygous global knockout of Slc7a5 was embryonically lethal. Muscle-specific (MCK Cre-mediated) Slc7a5 knockout (MS-Slc7a5-KO) mice were used to study the role of intracellular LNAA delivery by the SLC7A5 transporter for mTOR-S6K pathway activation in skeletal muscle. Activation of muscle mTOR-S6K (Thr389 phosphorylation) in vivo by intraperitoneal leucine injection was blunted in homozygous MS-Slc7a5-KO mice relative to wild-type animals. Dietary intake and growth rate were similar for MS-Slc7a5-KO mice and wild-type littermates fed for 10 weeks (to age 120 days) with diets containing 10%, 20% or 30% of protein. In MS-Slc7a5-KO mice, Leu and Ile concentrations in gastrocnemius muscle were reduced by ∼40% as dietary protein content was reduced from 30 to 10%. These changes were associated with >50% decrease in S6K Thr389 phosphorylation in muscles from MS-Slc7a5-KO mice, indicating reduced mTOR-S6K pathway activation, despite no significant differences in lean tissue mass between groups on the same diet. MS-Slc7a5-KO mice on 30% protein diet exhibited mild insulin resistance (e.g. reduced glucose clearance, larger gonadal adipose depots) relative to control animals. Thus, SLC7A5 modulates LNAA-dependent muscle mTOR-S6K signalling in mice, although it appears non-essential (or is sufficiently compensated by e.g. SLC7A8 (LAT2)) for maintenance of normal muscle mass

Fire decline in dry tropical ecosystems enhances decadal land carbon sink

The terrestrial carbon sink has significantly increased in the past decades, but the underlying mechanisms are still unclear. The current synthesis of process-based estimates of land and ocean sinks requires an additional sink of 0.6 PgC yr⁻¹ in the last decade to explain the observed airborne fraction. A concurrent global fire decline was observed in association with tropical agriculture expansion and landscape fragmentation. Here we show that a decline of 0.2 ± 0.1 PgC yr⁻¹ in fire emissions during 2008–2014 relative to 2001–2007 also induced an additional carbon sink enhancement of 0.4 ± 0.2 PgC yr⁻¹ attributable to carbon cycle feedbacks, amounting to a combined sink increase comparable to the 0.6 PgC yr⁻¹ budget imbalance. Our results suggest that the indirect effects of fire, in addition to the direct emissions, is an overlooked mechanism for explaining decadal-scale changes in the land carbon sink and highlight the importance of fire management in climate mitigation

Magnetism and its microscopic origin in iron-based high-temperature superconductors

High-temperature superconductivity in the iron-based materials emerges from, or sometimes coexists with, their metallic or insulating parent compound states. This is surprising since these undoped states display dramatically different antiferromagnetic (AF) spin arrangements and N$\rm \acute{e}$el temperatures. Although there is general consensus that magnetic interactions are important for superconductivity, much is still unknown concerning the microscopic origin of the magnetic states. In this review, progress in this area is summarized, focusing on recent experimental and theoretical results and discussing their microscopic implications. It is concluded that the parent compounds are in a state that is more complex than implied by a simple Fermi surface nesting scenario, and a dual description including both itinerant and localized degrees of freedom is needed to properly describe these fascinating materials.Comment: 14 pages, 4 figures, Review article, accepted for publication in Nature Physic

Frozen and Invariant Quantum Discord under Local Dephasing Noise

In this chapter, we intend to explore and review some remarkable dynamical properties of quantum discord under various different open quantum system models. Specifically, our discussion will include several concepts connected to the phenomena of time invariant and frozen quantum discord. Furthermore, we will elaborate on the relation of these two phenomena to the non-Markovian features of the open system dynamics and to the usage of dynamical decoupling protocols.Comment: 29 pages, 8 figure

Observation of a ppb mass threshoud enhancement in \psi^\prime\to\pi^+\pi^-J/\psi(J/\psi\to\gamma p\bar{p}) decay

The decay channel $\psi^\prime\to\pi^+\pi^-J/\psi(J/\psi\to\gamma p\bar{p})$ is studied using a sample of $1.06\times 10^8$ $\psi^\prime$ events collected by the BESIII experiment at BEPCII. A strong enhancement at threshold is observed in the $p\bar{p}$ invariant mass spectrum. The enhancement can be fit with an $S$-wave Breit-Wigner resonance function with a resulting peak mass of $M=1861^{+6}_{-13} {\rm (stat)}^{+7}_{-26} {\rm (syst)} {\rm MeV/}c^2$ and a narrow width that is $\Gamma<38 {\rm MeV/}c^2$ at the 90% confidence level. These results are consistent with published BESII results. These mass and width values do not match with those of any known meson resonance.Comment: 5 pages, 3 figures, submitted to Chinese Physics

Lattice Boltzmann simulations of soft matter systems

This article concerns numerical simulations of the dynamics of particles immersed in a continuum solvent. As prototypical systems, we consider colloidal dispersions of spherical particles and solutions of uncharged polymers. After a brief explanation of the concept of hydrodynamic interactions, we give a general overview over the various simulation methods that have been developed to cope with the resulting computational problems. We then focus on the approach we have developed, which couples a system of particles to a lattice Boltzmann model representing the solvent degrees of freedom. The standard D3Q19 lattice Boltzmann model is derived and explained in depth, followed by a detailed discussion of complementary methods for the coupling of solvent and solute. Colloidal dispersions are best described in terms of extended particles with appropriate boundary conditions at the surfaces, while particles with internal degrees of freedom are easier to simulate as an arrangement of mass points with frictional coupling to the solvent. In both cases, particular care has been taken to simulate thermal fluctuations in a consistent way. The usefulness of this methodology is illustrated by studies from our own research, where the dynamics of colloidal and polymeric systems has been investigated in both equilibrium and nonequilibrium situations.Comment: Review article, submitted to Advances in Polymer Science. 16 figures, 76 page