A semantic Grid for molecular science

Proceedings of the 2003 UK e-Science All Hands Meeting, 31st August - 3rd September, Nottingham UKThe properties of molecules have very well defined semantics and allow the creation of a semantic GRID. Markup languages (CML - Chemical Markup Language) and dictionary-based ontologies have been designed to support a wide range of applications, including chemical supply, publication and the safety of compounds. Many properties can be computed by Quantum Mechanical (QM) programs and we have developed a "black-box" system based on XML wrappers for all components. This is installed on a Condor system on which we have computed properties for 250, 000 compounds. The results of this will be available in an OpenData/OpenSource peer-to-peer (P2P) system (WorldWide Molecular Matrix - WWMM)

Computing the common zeros of two bivariate functions via Bézout resultants

The common zeros of two bivariate functions can be computed by finding the common zeros of their polynomial interpolants expressed in a tensor Chebyshev basis. From here we develop a bivariate rootfinding algorithm based on the hidden variable resultant method and Bézout matrices with polynomial entries. Using techniques including domain subdivision, Bézoutian regularization, and local refinement we are able to reliably and accurately compute the simple common zeros of two smooth functions with polynomial interpolants of very high degree (≥ 1000). We analyze the resultant method and its conditioning by noting that the Bézout matrices are matrix polynomials. Two implementations are available: one on the Matlab Central File Exchange and another in the roots command in Chebfun2 that is adapted to suit Chebfun’s methodology

LES-based Study of the Roughness Effects on the Wake of a Circular Cylinder from Subcritical to Transcritical Reynolds Numbers

This paper investigates the effects of surface roughness on the flow past a circular cylinder at subcritical to transcritical Reynolds numbers. Large eddy simulations of the flow for sand grain roughness of size k/D = 0.02 are performed (D is the cylinder diameter). Results show that surface roughness triggers the transition to turbulence in the boundary layer at all Reynolds numbers, thus leading to an early separation caused by the increased momentum deficit, especially at transcritical Reynolds numbers. Even at subcritical Reynolds numbers, boundary layer instabilities are triggered in the roughness sublayer and eventually lead to the transition to turbulence. The early separation at transcritical Reynolds numbers leads to a wake topology similar to that of the subcritical regime, resulting in an increased drag coefficient and lower Strouhal number. Turbulent statistics in the wake are also affected by roughness; the Reynolds stresses are larger due to the increased turbulent kinetic energy production in the boundary layer and separated shear layers close to the cylinder shoulders.We acknowledge “Red Española de Surpercomputación” (RES) for awarding us access to the MareNostrum III machine based in Barcelona, Spain (Ref. FI-2015-2-0026 and FI-2015-3-0011). We also acknowledge PRACE for awarding us access to Fermi and Marconi Supercomputers at Cineca, Italy (Ref. 2015133120). Oriol Lehmkuhl acknowledges a PDJ 2014 Grant by AGAUR (Generalitat de Catalunya). Ugo Piomelli acknowledges the support of the Natural Sciences and Engineering Research Council (NSERC) of Canada under the Discovery Grant Programme (Grant No. RGPIN-2016-04391). Ricard Borrell acknowledges a Juan de la Cierva postdoctoral grant (IJCI-2014-21034). Ivette Rodriguez, Oriol Lehmkuhl, Ricard Borrell and Assensi Oliva acknowledge Ministerio de Economía y Competitividad, Secretaría de Estado de Investigación, Desarrollo e Innovación, Spain (ref. ENE2014-60577-R).Peer ReviewedPostprint (author's final draft

Coalescent-based genome analyses resolve the early branches of the euarchontoglires

Despite numerous large-scale phylogenomic studies, certain parts of the mammalian tree are extraordinarily difficult to resolve. We used the coding regions from 19 completely sequenced genomes to study the relationships within the super-clade Euarchontoglires (Primates, Rodentia, Lagomorpha, Dermoptera and Scandentia) because the placement of Scandentia within this clade is controversial. The difficulty in resolving this issue is due to the short time spans between the early divergences of Euarchontoglires, which may cause incongruent gene trees. The conflict in the data can be depicted by network analyses and the contentious relationships are best reconstructed by coalescent-based analyses. This method is expected to be superior to analyses of concatenated data in reconstructing a species tree from numerous gene trees. The total concatenated dataset used to study the relationships in this group comprises 5,875 protein-coding genes (9,799,170 nucleotides) from all orders except Dermoptera (flying lemurs). Reconstruction of the species tree from 1,006 gene trees using coalescent models placed Scandentia as sister group to the primates, which is in agreement with maximum likelihood analyses of concatenated nucleotide sequence data. Additionally, both analytical approaches favoured the Tarsier to be sister taxon to Anthropoidea, thus belonging to the Haplorrhine clade. When divergence times are short such as in radiations over periods of a few million years, even genome scale analyses struggle to resolve phylogenetic relationships. On these short branches processes such as incomplete lineage sorting and possibly hybridization occur and make it preferable to base phylogenomic analyses on coalescent methods

Dynamics of direct inter-pack encounters in endangered African wild dogs

Aggressive encounters may have important life history consequences due to the potential for injury and death, disease transmission, dispersal opportunities or exclusion from key areas of the home range. Despite this, little is known of their detailed dynamics, mainly due to the difficulties of directly observing encounters in detail. Here, we describe detailed spatial dynamics of inter-pack encounters in African wild dogs (Lycaon pictus), using data from custom-built high-resolution GPS collars in 11 free-ranging packs. On average, each pack encountered another pack approximately every 7 weeks and met each neighbour twice each year. Surprisingly, intruders were more likely to win encounters (winning 78.6% of encounters by remaining closer to the site in the short term). However, intruders did tend to move farther than residents toward their own range core in the short-term (1 h) post-encounter, and if this were used to indicate losing an encounter, then the majority (73.3%) of encounters were won by residents. Surprisingly, relative pack size had little effect on encounter outcome, and injuries were rare (<15% of encounters). These results highlight the difficulty of remotely scoring encounters involving mobile participants away from static defendable food resources. Although inter-pack range overlap was reduced following an encounter, encounter outcome did not seem to drive this, as both packs shifted their ranges post-encounter. Our results indicate that inter-pack encounters may be lower risk than previously suggested and do not appear to influence long-term movement and ranging

Environmental Costs of Government-Sponsored Agrarian Settlements in Brazilian Amazonia

Brazil has presided over the most comprehensive agrarian reform frontier colonization program on Earth, in which ~1.2 million settlers have been translocated by successive governments since the 1970's, mostly into forested hinterlands of Brazilian Amazonia. These settlements encompass 5.3% of this ~5 million km2 region, but have contributed with 13.5% of all land conversion into agropastoral land uses. The Brazilian Federal Agrarian Agency (INCRA) has repeatedly claimed that deforestation in these areas largely predates the sanctioned arrival of new settlers. Here, we quantify rates of natural vegetation conversion across 1911 agrarian settlements allocated to 568 Amazonian counties and compare fire incidence and deforestation rates before and after the official occupation of settlements by migrant farmers. The timing and spatial distribution of deforestation and fires in our analysis provides irrefutable chronological and spatially explicit evidence of agropastoral conversion both inside and immediately outside agrarian settlements over the last decade. Deforestation rates are strongly related to local human population density and road access to regional markets. Agrarian settlements consistently accelerated rates of deforestation and fires, compared to neighboring areas outside settlements, but within the same counties. Relocated smallholders allocated to forest areas undoubtedly operate as pivotal agents of deforestation, and most of the forest clearance occurs in the aftermath of government-induced migration

Charting the landscape of N=4 flux compactifications

We analyse the vacuum structure of isotropic Z_2 x Z_2 flux compactifications, allowing for a single set of sources. Combining algebraic geometry with supergravity techniques, we are able to classify all vacua for both type IIA and IIB backgrounds with arbitrary gauge and geometric fluxes. Surprisingly, geometric IIA compactifications lead to a unique theory with four different vacua. In this case we also perform the general analysis allowing for sources compatible with minimal supersymmetry. Moreover, some relevant examples of type IIB non-geometric compactifications are studied. The computation of the full N=4 mass spectrum reveals the presence of a number of non-supersymmetric and nevertheless stable AdS_4 vacua. In addition we find a novel dS_4 solution based on a non-semisimple gauging.Comment: Minor corrections and references added. Version published in JHE

Simulations of extensional flow in microrheometric devices

We present a detailed numerical study of the flow of a Newtonian fluid through microrheometric devices featuring a sudden contraction–expansion. This flow configuration is typically used to generate extensional deformations and high strain rates. The excess pressure drop resulting from the converging and diverging flow is an important dynamic measure to quantify if the device is intended to be used as a microfluidic extensional rheometer. To explore this idea, we examine the effect of the contraction length, aspect ratio and Reynolds number on the flow kinematics and resulting pressure field. Analysis of the computed velocity and pressure fields show that, for typical experimental conditions used in microfluidic devices, the steady flow is highly three-dimensional with open spiraling vortical structures in the stagnant corner regions. The numerical simulations of the local kinematics and global pressure drop are in good agreement with experimental results. The device aspect ratio is shown to have a strong impact on the flow and consequently on the excess pressure drop, which is quantified in terms of the dimensionless Couette and Bagley correction factors. We suggest an approach for calculating the Bagley correction which may be especially appropriate for planar microchannels

Impact of a hospice rapid response service on preferred place of death, and costs

Background: Many people with a terminal illness would prefer to die at home. A new palliative rapid response service (RRS) provided by a large hospice provider in South East England was evaluated (2010) to provide evidence of impact on achieving preferred place of death and costs. The RRS was delivered by a team of trained health care assistants and available 24/7. The purpose of this study was to (i) compare the characteristics of RRS users and non-users, (ii) explore differences in the proportions of users and non-users dying in the place of their choice, (iii) monitor the whole system service utilisation of users and non-users, and compare costs. Methods: All hospice patients who died with a preferred place of death recorded during an 18 month period were included. Data (demographic, preferences for place of death) were obtained from hospice records. Dying in preferred place was modelled using stepwise logistic regression analysis. Service use data (period between referral to hospice and death) were obtained from general practitioners, community providers, hospitals, social services, hospice, and costs calculated using validated national tariffs. Results: Of 688 patients referred to the hospice when the RRS was operational, 247 (35.9 %) used it. Higher proportions of RRS users than non-users lived in their own homes with a co-resident carer (40.3 % vs. 23.7 %); more non-users lived alone or in residential care (58.8 % vs. 76.3 %). Chances of dying in the preferred place were enhanced 2.1 times by being a RRS user, compared to a non-user, and 1.5 times by having a co-resident carer, compared to living at home alone or in a care home. Total service costs did not differ between users and non-users, except when referred to hospice very close to death (users had higher costs). Conclusions: Use of the RRS was associated with increased likelihood of dying in the preferred place. The RRS is cost neutral