A Simple Model of Epidemics with Pathogen Mutation

We study how the interplay between the memory immune response and pathogen mutation affects epidemic dynamics in two related models. The first explicitly models pathogen mutation and individual memory immune responses, with contacted individuals becoming infected only if they are exposed to strains that are significantly different from other strains in their memory repertoire. The second model is a reduction of the first to a system of difference equations. In this case, individuals spend a fixed amount of time in a generalized immune class. In both models, we observe four fundamentally different types of behavior, depending on parameters: (1) pathogen extinction due to lack of contact between individuals, (2) endemic infection (3) periodic epidemic outbreaks, and (4) one or more outbreaks followed by extinction of the epidemic due to extremely low minima in the oscillations. We analyze both models to determine the location of each transition. Our main result is that pathogens in highly connected populations must mutate rapidly in order to remain viable.Comment: 9 pages, 11 figure

Real and imaginary chemical potential in 2-color QCD

In this paper we study the finite temperature SU(2) gauge theory with staggered fermions for non-zero imaginary and real chemical potential. The method of analytical continuation of Monte Carlo results from imaginary to real chemical potential is tested by comparison with simulations performed {\em directly} for real chemical potential. We discuss the applicability of the method in the different regions of the phase diagram in the temperature -- imaginary chemical potential plane.Comment: 15 pages, 7 figures; a few comments added; version published on Phys. Rev.

Gauge equivalence in QCD: the Weyl and Coulomb gauges

The Weyl-gauge ($A_0^a=0)$ QCD Hamiltonian is unitarily transformed to a representation in which it is expressed entirely in terms of gauge-invariant quark and gluon fields. In a subspace of gauge-invariant states we have constructed that implement the non-Abelian Gauss's law, this unitarily transformed Weyl-gauge Hamiltonian can be further transformed and, under appropriate circumstances, can be identified with the QCD Hamiltonian in the Coulomb gauge. We demonstrate an isomorphism that materially facilitates the application of this Hamiltonian to a variety of physical processes, including the evaluation of $S$-matrix elements. This isomorphism relates the gauge-invariant representation of the Hamiltonian and the required set of gauge-invariant states to a Hamiltonian of the same functional form but dependent on ordinary unconstrained Weyl-gauge fields operating within a space of ``standard'' perturbative states. The fact that the gauge-invariant chromoelectric field is not hermitian has important implications for the functional form of the Hamiltonian finally obtained. When this nonhermiticity is taken into account, the ``extra'' vertices in Christ and Lee's Coulomb-gauge Hamiltonian are natural outgrowths of the formalism. When this nonhermiticity is neglected, the Hamiltonian used in the earlier work of Gribov and others results.Comment: 25 page

Mesonic Wavefunctions in the three-dimensional Gross-Neveu model

We present results from a numerical study of bound state wavefunctions in the (2+1)-dimensional Gross-Neveu model with staggered lattice fermions at both zero and nonzero temperature. Mesonic channels with varying quantum numbers are identified and analysed. In the strongly coupled chirally broken phase at T=0 the wavefunctions expose effects due to varying the interaction strength more effectively than straightforward spectroscopy. In the weakly coupled chirally restored phase information on fermion - antifermion scattering is recovered. In the hot chirally restored phase we find evidence for a screened interaction. The T=0 chirally symmetric phase is most readily distinguished from the symmetric phase at high T via the fermion dispersion relation.Comment: 18 page

The spread of epidemic disease on networks

The study of social networks, and in particular the spread of disease on networks, has attracted considerable recent attention in the physics community. In this paper, we show that a large class of standard epidemiological models, the so-called susceptible/infective/removed (SIR) models can be solved exactly on a wide variety of networks. In addition to the standard but unrealistic case of fixed infectiveness time and fixed and uncorrelated probability of transmission between all pairs of individuals, we solve cases in which times and probabilities are non-uniform and correlated. We also consider one simple case of an epidemic in a structured population, that of a sexually transmitted disease in a population divided into men and women. We confirm the correctness of our exact solutions with numerical simulations of SIR epidemics on networks.Comment: 12 pages, 3 figure

Complex Ashtekar variables and reality conditions for Holst's action

From the Holst action in terms of complex valued Ashtekar variables additional reality conditions mimicking the linear simplicity constraints of spin foam gravity are found. In quantum theory with the results of You and Rovelli we are able to implement these constraints weakly, that is in the sense of Gupta and Bleuler. The resulting kinematical Hilbert space matches the original one of loop quantum gravity, that is for real valued Ashtekar connection. Our result perfectly fit with recent developments of Rovelli and Speziale concerning Lorentz covariance within spin-form gravity.Comment: 24 pages, 2 picture

Application of the Maximum Entropy Method to the (2+1)d Four-Fermion Model

We investigate spectral functions extracted using the Maximum Entropy Method from correlators measured in lattice simulations of the (2+1)-dimensional four-fermion model. This model is particularly interesting because it has both a chirally broken phase with a rich spectrum of mesonic bound states and a symmetric phase where there are only resonances. In the broken phase we study the elementary fermion, pion, sigma and massive pseudoscalar meson; our results confirm the Goldstone nature of the pi and permit an estimate of the meson binding energy. We have, however, seen no signal of sigma -> pi pi decay as the chiral limit is approached. In the symmetric phase we observe a resonance of non-zero width in qualitative agreement with analytic expectations; in addition the ultra-violet behaviour of the spectral functions is consistent with the large non-perturbative anomalous dimension for fermion composite operators expected in this model.Comment: 25 pages, 13 figure

Wilson function transforms related to Racah coefficients

The irreducible $*$-representations of the Lie algebra $su(1,1)$ consist of discrete series representations, principal unitary series and complementary series. We calculate Racah coefficients for tensor product representations that consist of at least two discrete series representations. We use the explicit expressions for the Clebsch-Gordan coefficients as hypergeometric functions to find explicit expressions for the Racah coefficients. The Racah coefficients are Wilson polynomials and Wilson functions. This leads to natural interpretations of the Wilson function transforms. As an application several sum and integral identities are obtained involving Wilson polynomials and Wilson functions. We also compute Racah coefficients for U_q(\su(1,1)), which turn out to be Askey-Wilson functions and Askey-Wilson polynomials.Comment: 48 page

Optimisation of reverse osmosis based wastewater treatment system for the removal of chlorophenol using genetic algorithms

YesReverse osmosis (RO) has found extensive applications in industry as an efficient separation process in comparison with thermal process. In this study, a one-dimensional distributed model based on a wastewater treatment spiral-wound RO system is developed to simulate the transport phenomena of solute and water through the membrane and describe the variation of operating parameters along the x-axis of membrane. The distributed model is tested against experimental data available in the literature derived from a chlorophenol rejection system implemented on a pilot-scale cross-flow RO filtration system with an individual spiral-wound membrane at different operating conditions. The proposed model is then used to carry out an optimisation study using a genetic algorithm (GA). The GA is developed to solve a formulated optimisation problem involving two objective functions of RO wastewater system performance. The model code is written in MATLAB, and the optimisation problem is solved using an optimisation platform written in C++. The objective function is to maximize the solute rejection at different cases of feed concentration and minimize the operating pressure to improve economic aspects. The operating feed flow rate, pressure and temperature are considered as decision variables. The optimisation problem is subjected to a number of upper and lower limits of decision variables, as recommended by the module’s manufacturer, and the constraint of the pressure loss along the membrane length to be within the allowable value. The algorithm developed has yielded a low optimisation execution time and resulted in improved unit performance based on a set of optimal operating conditions

Dynamical model and nonextensive statistical mechanics of a market index on large time windows

The shape and tails of partial distribution functions (PDF) for a financial signal, i.e. the S&P500 and the turbulent nature of the markets are linked through a model encompassing Tsallis nonextensive statistics and leading to evolution equations of the Langevin and Fokker-Planck type. A model originally proposed to describe the intermittent behavior of turbulent flows describes the behavior of normalized log-returns for such a financial market index, for small and large time windows, both for small and large log-returns. These turbulent market volatility (of normalized log-returns) distributions can be sufficiently well fitted with a $\chi^2$-distribution. The transition between the small time scale model of nonextensive, intermittent process and the large scale Gaussian extensive homogeneous fluctuation picture is found to be at $ca.$ a 200 day time lag. The intermittency exponent ($\kappa$) in the framework of the Kolmogorov log-normal model is found to be related to the scaling exponent of the PDF moments, -thereby giving weight to the model. The large value of $\kappa$ points to a large number of cascades in the turbulent process. The first Kramers-Moyal coefficient in the Fokker-Planck equation is almost equal to zero, indicating ''no restoring force''. A comparison is made between normalized log-returns and mere price increments.Comment: 40 pages, 14 figures; accepted for publication in Phys Rev