Search

OBJECTIVES: To provide estimates for how different treatment pathways for the management of severe aortic stenosis (AS) may affect National Health Service (NHS) England waiting list duration and associated mortality. DESIGN: We constructed a mathematical model of the excess waiting list and found the closed-form analytic solution to that model. From published data, we calculated estimates for how the strategies listed under Interventions may affect the time to clear the backlog of patients waiting for treatment and the associated waiting list mortality. SETTING: The NHS in England. PARTICIPANTS: Estimated patients with AS in England. INTERVENTIONS: (1) Increasing the capacity for the treatment of severe AS, (2) converting proportions of cases from surgery to transcatheter aortic valve implantation and (3) a combination of these two. RESULTS: In a capacitated system, clearing the backlog by returning to pre-COVID-19 capacity is not possible. A conversion rate of 50% would clear the backlog within 666 (533-848) days with 1419 (597-2189) deaths while waiting during this time. A 20% capacity increase would require 535 (434-666) days, with an associated mortality of 1172 (466-1859). A combination of converting 40% cases and increasing capacity by 20% would clear the backlog within a year (343 (281-410) days) with 784 (292-1324) deaths while awaiting treatment. CONCLUSION: A strategy change to the management of severe AS is required to reduce the NHS backlog and waiting list deaths during the post-COVID-19 'recovery' period. However, plausible adaptations will still incur a substantial wait to treatment and many hundreds dying while waiting., BG is supported by the NIHR Bristol Biomedical Research Centre at the University of Bristol and University Hospitals Bristol and Weston NHS Foundation Trust. JHFR is part-supported by the NIHR Cambridge Biomedical Research Centre, the British Heart Foundation, HEFCE, the EPSRC Cambridge Centre for Mathematics of Information in Healthcare and the Wellcome Trust.

Background From January to May 2021 the alpha variant (B.1.1.7) of SARS-CoV-2 was the most commonly detected variant in the UK. Following this, the Delta variant (B.1.617.2) then became the predominant variant. The UK COVID-19 vaccination programme started on 8th December 2020. Prior to the Delta variant, most vaccine effectiveness studies focused on the alpha variant. We therefore aimed to estimate the effectiveness of the BNT162b2 (Pfizer-BioNTech) and the ChAdOx1 nCoV-19 (Oxford-AstraZeneca) vaccines in preventing symptomatic and asymptomatic infection with respect to the Delta variant in a UK setting. Methods We used anonymised public health record data linked to infection data (PCR) using the Combined Intelligence for Population Health Action resource. We then constructed an SIR epidemic model to explain SARS-CoV-2 infection data across the Cheshire and Merseyside region of the UK. Vaccines were assumed to be effective after 21 days for 1 dose and 14 days for 2 doses. Results We determined that the effectiveness of the Oxford-AstraZeneca vaccine in reducing susceptibility to infection is 39% (95% credible interval [34, 43]) and 64% (95% credible interval [61, 67]) for a single dose and a double dose respectively. For the Pfizer-BioNTech vaccine, the effectiveness is 20% (95% credible interval [10, 28]) and 84% (95% credible interval [82, 86]) for a single-dose and a double dose respectively. Conclusion Vaccine effectiveness for reducing susceptibility to SARS-CoV-2 infection shows noticeable improvement after receiving two doses of either vaccine. Findings also suggest that a full course of the Pfizer-BioNTech provides the optimal protection against infection with the Delta variant. This reinforces the need to complete the full course programme to maximise individual protection and reduce transmission.

Our knowledge of viral host ranges remains limited. Completing this picture by identifying unknown hosts of known viruses is an important research aim that can help identify and mitigate zoonotic and animal-disease risks, such as spill-over from animal reservoirs into human populations. To address this knowledge-gap we apply a divide-and-conquer approach which separates viral, mammalian and network features into three unique perspectives, each predicting associations independently to enhance predictive power. Our approach predicts over 20,000 unknown associations between known viruses and susceptible mammalian species, suggesting that current knowledge underestimates the number of associations in wild and semi-domesticated mammals by a factor of 4.3, and the average potential mammalian host-range of viruses by a factor of 3.2. In particular, our results highlight a significant knowledge gap in the wild reservoirs of important zoonotic and domesticated mammals’ viruses: specifically, lyssaviruses, bornaviruses and rotaviruses., A more comprehensive map of viral host ranges can help identify and mitigate zoonotic and animal-disease risks. A divide-and-conquer approach which separates viral, mammalian and network features predicts over 20,000 unknown associations between known viruses and susceptible mammalian species.

Deterministic approximations to stochastic Susceptible–Infectious–Susceptible models typically predict a stable endemic steady-state when above threshold. This can be hard to relate to the underlying stochastic dynamics, which has no endemic steady-state but can exhibit approximately stable behaviour. Here, we relate the approximate models to the stochastic dynamics via the definition of the quasi-stationary distribution (QSD), which captures this approximately stable behaviour. We develop a system of ordinary differential equations that approximate the number of infected individuals in the QSD for arbitrary contact networks and parameter values. When the epidemic level is high, these QSD approximations coincide with the existing approximation methods. However, as we approach the epidemic threshold, the models deviate, with these models following the QSD and the existing methods approaching the all susceptible state. Through consistently approximating the QSD, the proposed methods provide a more robust link to the stochastic models.

Our knowledge of viral host ranges remains limited. Completing this picture by identifying unknown hosts of known viruses is an important research aim that can help identify zoonotic and animal-disease risks. Furthermore, such understanding can be used to mitigate against viral spill-over from animal reservoirs into human population. To address this knowledge-gap we apply a divide-and-conquer approach which separates viral, mammalian and network features into three unique perspectives, each predicting associations independently to enhance predictive power. Our approach predicts over 20,000 unknown associations between known viruses and mammalian hosts, suggesting that current knowledge underestimates the number of associations in wild and semi-domesticated mammals by a factor of 4.3, and the average mammalian host-range of viruses by a factor of 3.2. In particular, our results highlight a significant knowledge gap in the wild reservoirs of important zoonotic and domesticated mammals’ viruses: specifically, lyssaviruses, bornaviruses and rotaviruses.

Diseases that spread to humans from animals, zoonoses, pose major threats to human health. Identifying animal reservoirs of zoonoses and predicting future outbreaks are increasingly important to human health and well-being and economic stability, particularly where research and resources are limited. Here, we integrate complex networks and machine learning approaches to develop a new approach to identifying reservoirs. An exhaustive dataset of mammal–pathogen interactions was transformed into networks where hosts are linked via their shared pathogens. We present a methodology for identifying important and influential hosts in these networks. Ensemble models linking network characteristics with phylogeny and life-history traits are then employed to predict those key hosts and quantify the roles they undertake in pathogen transmission. Our models reveal drivers explaining host importance and demonstrate how these drivers vary by pathogen taxa. Host importance is further integrated into ensemble models to predict reservoirs of zoonoses of various pathogen taxa and quantify the extent of pathogen sharing between humans and mammals. We establish predictors of reservoirs of zoonoses, showcasing host influence to be a key factor in determining these reservoirs. Finally, we provide new insight into the determinants of zoonosis-sharing, and contrast these determinants across major pathogen taxa.

As ecosystems evolve, species can become extinct due to fluctuations in the environment. This leads to the evolutionary adaption known as bet-hedging, where species hedge against these fluctuations to reduce their likelihood of extinction. Environmental variation can be either within or between generations. Previous work has shown that selection for bet-hedging against within-generational variation should not occur in large populations. However, this work has been limited by assumptions of well-mixed populations, whereas real populations usually have some degree of structure. Using the framework of evolutionary graph theory, we show that through adding competition structure to the population, within-generational variation can have a significant impact on the evolutionary process for any population size. This complements research using subdivided populations, which suggests that within-generational variation is important when local population sizes are small. Together, these conclusions provide evidence to support observations by some ecologists that are contrary to the widely held view that only between-generational environmental variation has an impact on natural selection. This provides theoretical justification for further empirical study into this largely unexplored area. Supplementary Information The online version supplementary material available at 10.1007/s00285-021-01597-z.

Highlights • 87% of English and Welsh nodes are reachable from Scotland via the live fish movement network. • 72% of Scottish nodes are reachable from England or Wales via the live fish movement network. • 7.2% of all live fish movements cross the England-Scotland border. • Targeted surveillance on a handful of sites is effective in identifying and controlling outbreaks. • The combination of different mechanisms of transmission increases the chance of large epidemics., We analyse the network structure of the British salmonid aquaculture industry from the perspective of infectious disease control. We combine for the first time live fish transport (or movement) data covering England and Wales with data covering Scotland and include network layers representing potential transmission by rivers, sea water and local transmission via human or animal vectors in the immediate vicinity of each farm or fishery site. We find that 7.2% of all live fish transports cross the England-Scotland border and network analysis shows that 87% of English and Welsh nodes and 72% of Scottish nodes are reachable from cross-border connections via live fish transports alone. Consequently, from a disease-control perspective, the contact structures of England and Wales and of Scotland should not be considered in isolation. We also show that large epidemics require the live fish movement network and so control strategies targeting movements can be very effective. While there is relatively low risk of widespread epidemics on the live fish transport network alone, the potential risk is substantially amplified by the combined interaction of multiple network layers.

A biologically motivated individual-based framework for evolution in network-structured populations is developed that can accommodate eco-evolutionary dynamics. This framework is used to construct a network birth and death model. The evolutionary graph theory model, which considers evolutionary dynamics only, is derived as a special case, highlighting additional assumptions that diverge from real biological processes. This is achieved by introducing a negative ecological feedback loop that suppresses ecological dynamics by forcing births and deaths to be coupled. We also investigate how fitness, a measure of reproductive success used in evolutionary graph theory, is related to the life-history of individuals in terms of their birth and death rates. In simple networks, these ecologically motivated dynamics are used to provide new insight into the spread of adaptive mutations, both with and without clonal interference. For example, the star network, which is known to be an amplifier of selection in evolutionary graph theory, can inhibit the spread of adaptive mutations when individuals can die naturally.

We show that eigenvector centrality exhibits localization phenomena on networks that can be easily partitioned by the removal of a vertex cut set, the most extreme example being networks with a cut vertex. Three distinct types of localization are identified in these structures. One is related to the well-established hub node localization phenomenon and the other two are introduced and characterized here. We gain insights into these problems by deriving the relationship between eigenvector centrality and Katz centrality. This leads to an interpretation of the principal eigenvector as an approximation to more robust centrality measures which exist in the full span of an eigenbasis of the adjacency matrix.

Population structure can have a significant effect on evolution. For some systems with sufficient symmetry, analytic results can be derived within the mathematical framework of evolutionary graph theory which relate to the outcome of the evolutionary process. However, for more complicated heterogeneous structures, computationally intensive methods are required such as individual-based stochastic simulations. By adapting methods from statistical physics, including moment closure techniques, we first show how to derive existing homogenised pair approximation models and the exact neutral drift model. We then develop node-level approximations to stochastic evolutionary processes on arbitrarily complex structured populations represented by finite graphs, which can capture the different dynamics for individual nodes in the population. Using these approximations, we evaluate the fixation probability of invading mutants for given initial conditions, where the dynamics follow standard evolutionary processes such as the invasion process. Comparisons with the output of stochastic simulations reveal the effectiveness of our approximations in describing the stochastic processes and in predicting the probability of fixation of mutants on a wide range of graphs. Construction of these models facilitates a systematic analysis and is valuable for a greater understanding of the influence of population structure on evolutionary processes.

The duration of the infectious period is a crucial determinant of the ability of an infectious disease to spread. We consider an epidemic model that is network based and non-Markovian, containing classic Kermack-McKendrick, pairwise, message passing, and spatial models as special cases. For this model, we prove a monotonic relationship between the variability of the infectious period (with fixed mean) and the probability that the infection will reach any given subset of the population by any given time. For certain families of distributions, this result implies that epidemic severity is decreasing with respect to the variance of the infectious period. The striking importance of this relationship is demonstrated numerically. We then prove, with a fixed basic reproductive ratio (R_{0}), a monotonic relationship between the variability of the posterior transmission probability (which is a function of the infectious period) and the probability that the infection will reach any given subset of the population by any given time. Thus again, even when R_{0} is fixed, variability of the infectious period tends to dampen the epidemic. Numerical results illustrate this but indicate the relationship is weaker. We then show how our results apply to message passing, pairwise, and Kermack-McKendrick epidemic models, even when they are not exactly consistent with the stochastic dynamics. For Poissonian contact processes, and arbitrarily distributed infectious periods, we demonstrate how systems of delay differential equations and ordinary differential equations can provide upper and lower bounds, respectively, for the probability that any given individual has been infected by any given time.

Background: Hypoxia inducible factor-α (HIF-α) is the main transcription factor activated in low oxygen conditions. Results: Single cell imaging reveals pulses in nuclear levels of HIF-α. Conclusion: The transient nature of the HIF-α nuclear accumulation is required to avoid cell death. Significance: The duration of HIF-α response depends on cellular oxygenation, and can encode information and dictate cell fate., Intracellular signaling involving hypoxia-inducible factor (HIF) controls the adaptive responses to hypoxia. There is a growing body of evidence demonstrating that intracellular signals encode temporal information. Thus, the dynamics of protein levels, as well as protein quantity and/or localization, impacts on cell fate. We hypothesized that such temporal encoding has a role in HIF signaling and cell fate decisions triggered by hypoxic conditions. Using live cell imaging in a controlled oxygen environment, we observed transient 3-h pulses of HIF-1α and -2α expression under continuous hypoxia. We postulated that the well described prolyl hydroxylase (PHD) oxygen sensors and HIF negative feedback regulators could be the origin of the pulsatile HIF dynamics. We used iterative mathematical modeling and experimental analysis to scrutinize which parameter of the PHD feedback could control HIF timing and we probed for the functional redundancy between the three main PHD proteins. We identified PHD2 as the main PHD responsible for HIF peak duration. We then demonstrated that this has important consequences, because the transient nature of the HIF pulse prevents cell death by avoiding transcription of p53-dependent pro-apoptotic genes. We have further shown the importance of considering HIF dynamics for coupling mathematical models by using a described HIF-p53 mathematical model. Our results indicate that the tight control of HIF transient dynamics has important functional consequences on the cross-talk with key signaling pathways controlling cell survival, which is likely to impact on HIF targeting strategies for hypoxia-associated diseases such as tumor progression and ischemia.

Methods for efficiently controlling dynamics propagated on networks are usually based on identifying the most influential nodes. Knowledge of these nodes can be used for the targeted control of dynamics such as epidemics, or for modifying biochemical pathways relating to diseases. Similarly they are valuable for identifying points of failure to increase network resilience in, for example, social support networks and logistics networks. Many measures, often termed ‘centrality’, have been constructed to achieve these aims. Here we consider Katz centrality and provide a new interpretation as a steady-state solution to continuous-time dynamics. This enables us to implement a sensitivity analysis which is similar to metabolic control analysis used in the analysis of biochemical pathways. The results yield a centrality which quantifies, for each node, the net impact of its absence from the network. It also has the desirable property of requiring a node with a high centrality to play a central role in propagating the dynamics of the system by having the capacity to both receive flux from others and then to pass it on. This new perspective on Katz centrality is important for a more comprehensive analysis of directed networks.

We present results for the static inter-quark potential, lightest glueballs, light hadron spectrum and topological susceptibility using a non-perturbatively improved action on a $16^3\times 32$ lattice at a set of values of the bare gauge coupling and bare dynamical quark mass chosen to keep the lattice size fixed in physical units ($\sim 1.7$ fm). By comparing these measurements with a matched quenched ensemble, we study the effects due to two degenerate flavours of dynamical quarks. With the greater control over residual lattice spacing effects which these methods afford, we find some evidence of charge screening and some minor effects on the light hadron spectrum over the range of quark masses studied ($M_{PS}/M_{V}\ge0.58$). More substantial differences between quenched and unquenched simulations are observed in measurements of topological quantities., Comment: 53 pages, LaTeX/RevTeX, 16 eps figures; corrected clover action expression and various typos, no results changed

The relationship between system-level and subsystem-level master equations is investigated and then utilised for a systematic and potentially automated derivation of the hierarchy of moment equations in a susceptible-infectious-removed (SIR) epidemic model. In the context of epidemics on contact networks we use this to show that the approximate nature of some deterministic models such as mean-field and pair-approximation models can be partly understood by the identification of implicit anomalous terms. These terms describe unbiological processes which can be systematically removed up to and including the nth order by nth order moment closure approximations. These terms lead to a detailed understanding of the correlations in network-based epidemic models and contribute to understanding the connection between individual-level epidemic processes and population-level models. The connection with metapopulation models is also discussed. Our analysis is predominantly made at the individual level where the first and second order moment closure models correspond to what we term the individual-based and pair-based deterministic models, respectively. Matlab code is included as supplementary material for solving these models on transmission networks of arbitrary complexity.

Epidemics are frequently simulated on redundantly wired contact networks, which have many more links between sites than are minimally required to connect all. Consequently, the modelled pathogen can travel numerous alternative routes, complicating effective containment strategies. These networks have moreover been found to exhibit ‘scale-free’ properties and percolation, suggesting resilience to damage. However, realistic H5N1 avian influenza transmission probabilities and containment strategies, here modelled on the British poultry industry network, show that infection dynamics can additionally express characteristic scales. These system-preferred scales constitute small areas within an observed power law distribution that exhibit a lesser slope than the power law itself, indicating a slightly increased relative likelihood. These characteristic scales are here produced by a network-pervading intranet of so-called hotspot sites that propagate large epidemics below the percolation threshold. This intranet is, however, extremely vulnerable; targeted inoculation of a mere 3–6% (depending on incorporated biosecurity measures) of the British poultry industry network prevents large and moderate H5N1 outbreaks completely, offering an order of magnitude improvement over previously advocated strategies affecting the most highly connected ‘hub’ sites. In other words, hotspots and hubs are separate functional entities that do not necessarily coincide, and hotspots can make more effective inoculation targets. Given the ubiquity and relevance of networks (epidemics, Internet, power grids, protein interaction), recognition of this spreading regime elsewhere would suggest a similar disproportionate sensitivity to such surgical interventions.

We prove that, for Poisson transmission and recovery processes, the classic susceptible→infected→recovered (SIR) epidemic model of Kermack and McKendrick provides, for any given time t>0, a strict lower bound on the expected number of susceptibles and a strict upper bound on the expected number of recoveries in the general stochastic SIR epidemic. The proof is based on the recent message passing representation of SIR epidemics applied to a complete graph.

We first generalise ideas discussed by Kiss et al. (2015) to prove a theorem for generating exact closures (here expressing joint probabilities in terms of their constituent marginal probabilities) for susceptible-infectious-removed (SIR) dynamics on arbitrary graphs (networks). For Poisson transmission and removal processes, this enables us to obtain a systematic reduction in the number of differential equations needed for an exact 'moment closure' representation of the underlying stochastic model. We define 'transmission blocks' as a possible extension of the block concept in graph theory and show that the order at which the exact moment closure representation is curtailed is the size of the largest transmission block. More generally, approximate closures of the hierarchy of moment equations for these dynamics are typically defined for the first and second order yielding mean-field and pairwise models respectively. It is frequently implied that, in principle, closed models can be written down at arbitrary order if only we had the time and patience to do this. However, for epidemic dynamics on networks, these higher-order models have not been defined explicitly. Here we unambiguously define hierarchies of approximate closed models that can utilise subsystem states of any order, and show how well-known models are special cases of these hierarchies.

In cases where there are limited resources for the eradication of an epidemic, or where we seek to minimise possible adverse impacts of interventions, it is essential to optimise the efficacy of control measures. We introduce a new approach, Epidemic Control Analysis (ECA), to design effective targeted intervention strategies to mitigate and control the propagation of infections across heterogeneous contact networks. We exemplify this methodology in the context of a newly developed individual-level deterministic Susceptible-Infectious-Susceptible (SIS) epidemiological model (we also briefly consider applications to Susceptible-Infectious-Removed (SIR) dynamics). This provides a flexible way to systematically determine the impact of interventions on endemic infections in the population. Individuals are ranked based on their influence on the level of infectivity. The highest-ranked individuals are prioritised for targeted intervention. Many previous intervention strategies have determined prioritisation based mainly on the position of individuals in the network, described by various local and global network centrality measures, and their chance of being infectious. Comparisons of the predictions of the proposed strategy with those of widely used targeted intervention programmes on various model and real-world networks reveal its efficiency and accuracy. It is demonstrated that targeting central individuals or individuals that have high infection probability is not always the best strategy. The importance of individuals is not determined by network structure alone, but can be highly dependent on the infection dynamics. This interplay between network structure and infection dynamics is effectively captured by ECA.

The message passing approach of Karrer and Newman [Phys. Rev. E 82, 016101 (2010)] is an exact and practicable representation of susceptible-infected-recovered dynamics on finite trees. Here we show that, assuming Poisson contact processes, a pair-based moment closure representation [Sharkey, J. Math. Biol. 57, 311 (2008)] can be derived from their equations. We extend the applicability of both representations and discuss their relative merits. On arbitrary time-independent networks, as was shown for the message passing formalism, the pair-based moment closure equations also provide a rigorous lower bound on the expected number of susceptibles at all times., Comment: 7 pages, 2 figures

Understanding models which represent the invasion of network-based systems by infectious agents can give important insights into many real-world situations, including the prevention and control of infectious diseases and computer viruses. Here we consider Markovian susceptible-infectious-susceptible (SIS) dynamics on finite strongly connected networks, applicable to several sexually transmitted diseases and computer viruses. In this context, a theoretical definition of endemic prevalence is easily obtained via the quasi-stationary distribution (QSD). By representing the model as a percolation process and utilising the property of duality, we also provide a theoretical definition of invasion probability. We then show that, for undirected networks, the probability of invasion from any given individual is equal to the (probabilistic) endemic prevalence, following successful invasion, at the individual (we also provide a relationship for the directed case). The total (fractional) endemic prevalence in the population is thus equal to the average invasion probability (across all individuals). Consequently, for such systems, the regions or individuals already supporting a high level of infection are likely to be the source of a successful invasion by another infectious agent. This could be used to inform targeted interventions when there is a threat from an emerging infectious disease., Comment: 16 pages, 5 figures. Supplementary data available with published version at http://dx.doi.org/10.1371/journal.pone.0069028

The susceptibility of the English and Welsh fish farming and fisheries industry to emergent diseases is assessed using a stochastic simulation model. The model dynamics operate on a network comprising directed transport and river contacts, as well as undirected local and fomite transmissions. The directed connections cause outward transmission risk to be geographically more confined than inward risk. We consider reactive, proactive, and hybrid methods of control which correspond to a mixture of policy and the ease of disease detection. An explicit investigation of the impact of laboratory capacity is made. General quantified guidelines are derived to mitigate future epidemics.

SulfoSYS (Sulfolobus Systems Biology) focuses on the study of the CCM (central carbohydrate metabolism) of Sulfolobus solfataricus and its regulation under temperature variation at the systems level. In Archaea, carbohydrates are metabolized by modifications of the classical pathways known from Bacteria or Eukarya, e.g. the unusual branched ED (Entner-Doudoroff) pathway, which is utilized for glucose degradation in S. solfataricus. This archaeal model organism of choice is a thermoacidophilic crenarchaeon that optimally grows at 80°C (60-92°C) and pH 2-4. In general, life at high temperature requires very efficient adaptation to temperature changes, which is most difficult to deal with for organisms, and it is unclear how biological networks can withstand and respond to such changes. This integrative project combines genomic, transcriptomic, proteomic and metabolomic, as well as kinetic and biochemical information. The final goal of SulfoSYS is the construction of a silicon cell model for this part of the living cell that will enable computation of the CCM network. In the present paper, we report on one of the first archaeal systems biology projects. © The Authors Journal compilation © 2009 Biochemical Society.

Highly pathogenic avian influenza and in particular the H5N1 strain has resulted in the culling of millions of birds and continues to pose a threat to poultry industries worldwide. The recent outbreak of H5N1 in the UK highlights the need for detailed assessment of the consequences of an incursion and of the efficacy of control strategies. Here, we present results from a model of H5N1 propagation within the British poultry industry. We find that although the majority of randomly seeded incursions do not spread beyond the initial infected premises, there is significant potential for widespread infection. The efficacy of the European Union strategy for disease control is evaluated and our simulations emphasize the pivotal role of duck farms in spreading H5N1.

The process of infection during an epidemic can be envisaged as being transmitted via a network of routes represented by a contact network. Most differential equation models of epidemics are mean-field models. These contain none of the underlying spatial structure of the contact network. By extending the mean-field models to pair-level, some of the spatial structure can be contained in the model. Some networks of transmission such as river or transportation networks are clearly asymmetric, whereas others such as airborne infection can be regarded as symmetric. Pair-level models have been developed to describe symmetric contact networks. Here we report on work to develop a pair-level model that is also applicable to asymmetric contact networks. The procedure for closing the model at the level of pairs is discussed in detail. The model is compared against stochastic simulations of epidemics on asymmetric contact networks and against the predictions of the symmetric model on the same networks.

As existing hydrocarbon reservoirs are developed and exploited, the need to improve techniques for identifying new reservoirs and describing existing reservoirs grows. Indeed, it has been reported that unless a step change in resolution is achieved, production from the North Sea fields will start to decline mid decade. This paper explores the potential for tracking techniques to achieve that step change.

We report on very strong evidence of the occurrence of power terms in αs(p), the QCD running coupling constant in the MOM scheme, by analyzing non-perturbative measurements from the lattice three-gluon vertex between 2.0 and 10.0 GeV at zero flavor. While putting forward the caveat that this definition of the coupling is a gauge dependent one, the general relevance of such an occurrence is discussed. We fit ΛMSbar(nf = 0) = 237±3 + 0−10 MeV in perfect agreement with the result obtained by the ALPHA group with a totally different method. The power correction to αs(p) is fitted to (0.63±0.03 + 0.0− 0.13) GeV2/p2.

In many fields of science including population dynamics, the vast state spaces inhabited by all but the very simplest of systems can preclude a deterministic analysis. Here, a class of approximate deterministic models is introduced into the field of epidemiology that reduces this state space to one that is numerically feasible. However, these reduced state space master equations do not in general form a closed set. To resolve this, the equations are approximated using closure approximations. This process results in a method for constructing deterministic differential equation models with a potentially large scope of application including dynamic directed contact networks and heterogeneous systems using time dependent parameters. The method is exemplified in the case of an SIR (susceptible-infectious-removed) epidemiological model and is numerically evaluated on a range of networks from spatially local to random. In the context of epidemics propagated on contact networks, this work assists in clarifying the link between stochastic simulation and traditional population level deterministic models.