Diffusion and Trapping on a one-dimensional lattice

The properties of a particle diffusing on a one-dimensional lattice where at each site a random barrier and a random trap act simultaneously on the particle are investigated by numerical and analytical techniques. The combined effect of disorder and traps yields a decreasing survival probability with broad distribution (log-normal). Exact enumerations, effective-medium approximation and spectral analysis are employed. This one-dimensional model shows rather rich behaviours which were previously believed to exist only in higher dimensionality. The possibility of a trapping-dominated super universal class is suggested.Comment: 20 pages, Revtex 3.0, 13 figures in compressed format using uufiles command, to appear in Phys. Rev. E, for an hard copy or problems e-mail to: [email protected]

Diffusion with critically correlated traps and the slow relaxation of the longest wavelength mode

We study diffusion on a substrate with permanent traps distributed with critical positional correlation, modeled by their placement on the perimeters of a critical percolation cluster. We perform a numerical analysis of the vibrational density of states and the largest eigenvalue of the equivalent scalar elasticity problem using the method of Arnoldi and Saad. We show that the critical trap correlation increases the exponent appearing in the stretched exponential behavior of the low frequency density of states by approximately a factor of two as compared to the case of no correlations. A finite size scaling hypothesis of the largest eigenvalue is proposed and its relation to the density of states is given. The numerical analysis of this scaling postulate leads to the estimation of the stretch exponent in good agreement with the density of states result.Comment: 15 pages, LaTeX (RevTeX

A pseudo-spectral method for the Kardar-Parisi-Zhang equation

We discuss a numerical scheme to solve the continuum Kardar-Parisi-Zhang equation in generic spatial dimensions. It is based on a momentum-space discretization of the continuum equation and on a pseudo-spectral approximation of the non-linear term. The method is tested in (1+1)- and (2+1)- dimensions, where it is shown to reproduce the current most reliable estimates of the critical exponents based on Restricted Solid-on-Solid simulations. In particular it allows the computations of various correlation and structure functions with high degree of numerical accuracy. Some deficiencies which are common to all previously used finite-difference schemes are pointed out and the usefulness of the present approach in this respect is discussed.Comment: 12 pages, 13 .eps figures, revetx4. A few equations have been corrected. Erratum sent to Phys. Rev.

A numerical study of a binary Yukawa model in regimes characteristic of globular proteins in solutions

The main goal of this paper is to assess the limits of validity, in the regime of low concentration and strong Coulomb coupling (high molecular charges), for a simple perturbative approximation to the radial distribution functions (RDF), based upon a low-density expansion of the potential of mean force and proposed to describe protein-protein interactions in a recent Small-Angle-Scattering (SAS) experimental study. A highly simplified Yukawa (screened Coulomb) model of monomers and dimers of a charged globular protein ($\beta$-lactoglobulin) in solution is considered. We test the accuracy of the RDF approximation, as a necessary complementary part of the previous experimental investigation, by comparison with the fluid structure predicted by approximate integral equations and exact Monte Carlo (MC) simulations. In the MC calculations, an Ewald construction for Yukawa potentials has been used to take into account the long-range part of the interactions in the weakly screened cases. Our results confirm that the perturbative first-order approximation is valid for this system even at strong Coulomb coupling, provided that the screening is not too weak (i.e., for Debye length smaller than monomer radius). A comparison of the MC results with integral equation calculations shows that both the hypernetted-chain (HNC) and the Percus-Yevick (PY) closures have a satisfactory behavior under these regimes, with the HNC being superior throughout. The relevance of our findings for interpreting SAS results is also discussed.Comment: Physical Review E, in press (2005

Prevalence of Mycobacterium avium subsp. paratuberculosis in milk and dairy cattle in Southern Italy; preliminary results

Paratuberculosis affects all ruminants worldwide. Mycobacterium avium subsp. paratuberculosis could have a role in human diseases like Crohn\u2019s. Some extra EU countries request importation of MAP-free products. Italy has not yet actualized a control program and the diffusion of the infection is still unknown in Southern Italy. The aim of this study was to evaluate the prevalence of the infection in five regions of Southern Italy. Bulk tank milk and in-line milk filters were sampled in 780 dairy cattle herds and respectively analyzed by ELISA and real time PCR. One hundred and fifty-five out of 780 herds (19.9%) were found positive by ELISA and/or real time PCR. Individual milk samples were then collected from all the producing animals of positive herds and from a selection of negative herds. The estimated prevalence varies from region to region between 2.8% and 5.5%. Our results indicate that the disease is widespread in the five regions. The observed prevalence could be underestimated

Microscopic formulation of the Zimm-Bragg model for the helix-coil transition

A microscopic spin model is proposed for the phenomenological Zimm-Bragg model for the helix-coil transition in biopolymers. This model is shown to provide the same thermophysical properties of the original Zimm-Bragg model and it allows a very convenient framework to compute statistical quantities. Physical origins of this spin model are made transparent by an exact mapping into a one-dimensional Ising model with an external field. However, the dependence on temperature of the reduced external field turns out to differ from the standard one-dimensional Ising model and hence it gives rise to different thermophysical properties, despite the exact mapping connecting them. We discuss how this point has been frequently overlooked in the recent literature.Comment: 11 pages, 2 figure

Effect of Polydispersity and Anisotropy in Colloidal and Protein Solutions: an Integral Equation Approach

Application of integral equation theory to complex fluids is reviewed, with particular emphasis to the effects of polydispersity and anisotropy on their structural and thermodynamic properties. Both analytical and numerical solutions of integral equations are discussed within the context of a set of minimal potential models that have been widely used in the literature. While other popular theoretical tools, such as numerical simulations and density functional theory, are superior for quantitative and accurate predictions, we argue that integral equation theory still provides, as in simple fluids, an invaluable technique that is able to capture the main essential features of a complex system, at a much lower computational cost. In addition, it can provide a detailed description of the angular dependence in arbitrary frame, unlike numerical simulations where this information is frequently hampered by insufficient statistics. Applications to colloidal mixtures, globular proteins and patchy colloids are discussed, within a unified framework.Comment: 17 pages, 7 figures, to appear in Interdiscip. Sci. Comput. Life Sci. (2011), special issue dedicated to Prof. Lesser Blu

Interstellar dust in the BOOMERanG maps

Interstellar dust (ISD) emission is present in the mm-wave maps obtained by the BOOMERanG experiment at intermediate and high Galactic latitudes. We find that, while being sub-dominant at the lower frequencies (90,150, 240 GHz), thermal emission from ISD is dominant at 410 GHz, and is well correlated with the IRAS map at 100 µm. We find also that the angular power spectrum of ISD fluctuations at 410 GHz is a power law, and its level is negligible with respect to the angular power spectrum of the Cosmic Microwave Background (CMB) at 90 and 150 GHz

Behaviour of Escherichia coli O157:H7 during the manufacture and ripening of an Italian traditional raw goat milk cheese

Formagelle di capra is a raw goat cheese produced from whole chilled goat milk; traditional technology involving unpasteurised milk and indigenous lactic starter cultures is employed for its production in Italy. The purpose of this study was to assess the behavior of Escherichia coli O157:H7 during the manufacturing and ripening of this raw goat milk cheese. Raw milk was experimentally inoculated with E. coli O157:H7 in a laboratory scale plant and the count was monitored during production and 30 days of ripening required for this cheese. Results showed that E. coli O157:H7 count increased to more than 1.5 Log cfu g-1 during cheese production and remained constant until the end of ripening. The evidence that E. coli O157:H7 is able to survive during the manufacturing and ripening process suggests that the 30-day ripening period alone is insufficient to eliminate levels of viable E. coli O157:H7 in Formaggelle di capra cheese and that the presence of low numbers of E. coli O157:H7 in milk destined for the production of raw goat milk cheeses could represent a potential source of infection for humans and a threat for consumer

Self-assembly mechanism in colloids: perspectives from Statistical Physics

Motivated by recent experimental findings in chemical synthesis of colloidal particles, we draw an analogy between self-assembly processes occurring in biological systems (e.g. protein folding) and a new exciting possibility in the field of material science. We consider a self-assembly process whose elementary building blocks are decorated patchy colloids of various types, that spontaneously drive the system toward a unique and predetermined targeted macroscopic structure. To this aim, we discuss a simple theoretical model -- the Kern-Frenkel model -- describing a fluid of colloidal spherical particles with a pre-defined number and distribution of solvophobic and solvophilic regions on their surface. The solvophobic and solvophilic regions are described via a short-range square-well and a hard-sphere potentials, respectively. Integral equation and perturbation theories are presented to discuss structural and thermodynamical properties, with particular emphasis on the computation of the fluid-fluid (or gas-liquid) transition in the temperature-density plane. The model allows the description of both one and two attractive caps, as a function of the fraction of covered attractive surface, thus interpolating between a square-well and a hard-sphere fluid, upon changing the coverage. By comparison with Monte Carlo simulations, we assess the pros and the cons of both integral equation and perturbation theories in the present context of patchy colloids, where the computational effort for numerical simulations is rather demanding.Comment: 14 pages, 7 figures, Special issue for the SigmaPhi2011 conferenc