Sustainable management of miombo woodlands in the Northern part of Mozambique (Niassa National Reserve - NNR).

Poster presented at Commiting Science to Global Development. Lisbon (Portugal). 29-30 Sep 2009

Towards a method for rigorous development of generic requirements patterns

We present work in progress on a method for the engineering, validation and verification of generic requirements using domain engineering and formal methods. The need to develop a generic requirement set for subsequent system instantiation is complicated by the addition of the high levels of verification demanded by safety-critical domains such as avionics. Our chosen application domain is the failure detection and management function for engine control systems: here generic requirements drive a software product line of target systems. A pilot formal specification and design exercise is undertaken on a small (twosensor) system element. This exercise has a number of aims: to support the domain analysis, to gain a view of appropriate design abstractions, for a B novice to gain experience in the B method and tools, and to evaluate the usability and utility of that method.We also present a prototype method for the production and verification of a generic requirement set in our UML-based formal notation, UML-B, and tooling developed in support. The formal verification both of the structural generic requirement set, and of a particular application, is achieved via translation to the formal specification language, B, using our U2B and ProB tools

Computation of thermodynamic and transport properties to predict thermophoretic effects in an argon-krypton mixture

Thermophoresis is the movement of molecules caused by a temperature gradient. Here we report the results of a study of thermophoresis using non-equilibrium molecular dynamics simulations of a confined argon-krypton fluid subject to two different temperatures at thermostated walls. The resulting temperature profile between the walls is used along with the Soret coefficient to predict the concentration profile that develops across the channel. We obtain the Soret coefficient by calculating the mutual diffusion and thermal diffusion coefficients. We report an appropriate method for calculating the transport coefficients for binary systems, using the Green-Kubo integrals and radial distribution functions obtained from equilibrium molecular dynamics simulations of the bulk fluid. Our method has the unique advantage of separating the mutual diffusion and thermal diffusion coef- ficients, and calculating the sign and magnitude of their individual contributions to thermophoresis in binary mixtures

Asymptotic analysis for the generalized langevin equation

Various qualitative properties of solutions to the generalized Langevin equation (GLE) in a periodic or a confining potential are studied in this paper. We consider a class of quasi-Markovian GLEs, similar to the model that was introduced in \cite{EPR99}. Geometric ergodicity, a homogenization theorem (invariance principle), short time asymptotics and the white noise limit are studied. Our proofs are based on a careful analysis of a hypoelliptic operator which is the generator of an auxiliary Markov process. Systematic use of the recently developed theory of hypocoercivity \cite{Vil04HPI} is made.Comment: 27 pages, no figures. Submitted to Nonlinearity

New Langevin and Gradient Thermostats for Rigid Body Dynamics

We introduce two new thermostats, one of Langevin type and one of gradient (Brownian) type, for rigid body dynamics. We formulate rotation using the quaternion representation of angular coordinates; both thermostats preserve the unit length of quaternions. The Langevin thermostat also ensures that the conjugate angular momenta stay within the tangent space of the quaternion coordinates, as required by the Hamiltonian dynamics of rigid bodies. We have constructed three geometric numerical integrators for the Langevin thermostat and one for the gradient thermostat. The numerical integrators reflect key properties of the thermostats themselves. Namely, they all preserve the unit length of quaternions, automatically, without the need of a projection onto the unit sphere. The Langevin integrators also ensure that the angular momenta remain within the tangent space of the quaternion coordinates. The Langevin integrators are quasi-symplectic and of weak order two. The numerical method for the gradient thermostat is of weak order one. Its construction exploits ideas of Lie-group type integrators for differential equations on manifolds. We numerically compare the discretization errors of the Langevin integrators, as well as the efficiency of the gradient integrator compared to the Langevin ones when used in the simulation of rigid TIP4P water model with smoothly truncated electrostatic interactions. We observe that the gradient integrator is computationally less efficient than the Langevin integrators. We also compare the relative accuracy of the Langevin integrators in evaluating various static quantities and give recommendations as to the choice of an appropriate integrator.Comment: 16 pages, 4 figure

A lesson on interrogations from detainees: Predicting self-reported confessions and cooperation

The ability to predict confessions and cooperation from the elements of an interrogation was examined. Incarcerated men (N = 100) completed a 50-item questionnaire about their most recent police interrogation, and regression analyses were performed on self-reported decisions to confess and cooperate. Results showed that the likelihood of an interrogation resulting in a confession was greatest when evidence strength and score on a humanitarian interviewing scale were high, and when the detainee had few previous convictions or did not seek legal advice. We also found that the level of cooperation was greatest when the humanitarian interviewing score was high, and when previous convictions were low. The implications of the findings for interrogation practices are discussed

Boundary lubrication properties of materials with expansive freezing

We have performed molecular dynamics simulations of solid-solid contacts lubricated by a model fluid displaying many of the properties of water, particularly its expansive freezing. Near the region where expansive freezing occurs, the lubricating film remains fluid, and the friction force decreases linearly as the shear velocity is reduced. No sign of stick-slip motion is observed even at the lowest velocities. We give a simple interpretation of these results, and suggest that in general good boundary lubrication properties will be found in the family of materials with expansive freezing.Comment: Version to appear in Phys. Rev. Let

Postcopulatory sexual selection

The female reproductive tract is where competition between the sperm of different males takes place, aided and abetted by the female herself. Intense postcopulatory sexual selection fosters inter-sexual conflict and drives rapid evolutionary change to generate a startling diversity of morphological, behavioural and physiological adaptations. We identify three main issues that should be resolved to advance our understanding of postcopulatory sexual selection. We need to determine the genetic basis of different male fertility traits and female traits that mediate sperm selection; identify the genes or genomic regions that control these traits; and establish the coevolutionary trajectory of sexes

Epithelial Immunization Induces Polyfunctional CD8+ T Cells and Optimal Mousepox Protection.

We assessed several routes of immunization with vaccinia virus (VACV) in protecting mice against ectromelia virus (ECTV). By a wide margin, skin scarification provided the greatest protection. Humoral immunity and resident-memory T cells notwithstanding, several approaches revealed that circulating, memory CD8(+) T cells primed via scarification were functionally superior and conferred enhanced virus control. Immunization via the epithelial route warrants further investigation, as it may also provide enhanced defense against other infectious agents

General framework for fluctuating dynamic density functional theory

We introduce a versatile bottom-up derivation of a formal theoretical framework to describe (passive) soft-matter systems out of equilibrium subject to fluctuations. We provide a unique connection between the constituent-particle dynamics of real systems and the time evolution equation of their measurable (coarse-grained) quantities, such as local density and velocity. The starting point is the full Hamiltonian description of a system of colloidal particles immersed in a fluid of identical bath particles. Then, we average out the bath via Zwanzig's projection-operator techniques and obtain the stochastic Langevin equations governing the colloidal-particle dynamics. Introducing the appropriate definition of the local number and momentum density fields yields a generalisation of the Dean-Kawasaki (DK) model, which resembles the stochastic Navier-Stokes (NS) description of a fluid. Nevertheless, the DK equation still contains all the microscopic information and, for that reason, does not represent the dynamical law of observable quantities. We address this controversial feature of the DK description by carrying out a nonequilibrium ensemble average. Adopting a natural decomposition into local-equilibrium and nonequilibrium contribution, where the former is related to a generalised version of the canonical distribution, we finally obtain the fluctuating-hydrodynamic equation governing the time-evolution of the mesoscopic density and momentum fields. Along the way, we outline the connection between the ad-hoc energy functional introduced in previous DK derivations and the free-energy functional from classical density-functional theory (DFT). The resultant equation has the structure of a dynamical DFT (DDFT) with an additional fluctuating force coming from the random interactions with the bath. We show that our fluctuating DDFT formalism corresponds to a particular version of the fluctuating NS equations, originally derived by Landau and Lifshitz. Our framework thus provides the formal apparatus for ab-initio derivations of fluctuating DDFT equations capable of describing the dynamics of soft-matter systems in and out of equilibrium. We believe that the derivation offered here represents the current state of the art in the field