Survival probability and order statistics of diffusion on disordered media

We investigate the first passage time t_{j,N} to a given chemical or Euclidean distance of the first j of a set of N>>1 independent random walkers all initially placed on a site of a disordered medium. To solve this order-statistics problem we assume that, for short times, the survival probability (the probability that a single random walker is not absorbed by a hyperspherical surface during some time interval) decays for disordered media in the same way as for Euclidean and some class of deterministic fractal lattices. This conjecture is checked by simulation on the incipient percolation aggregate embedded in two dimensions. Arbitrary moments of t_{j,N} are expressed in terms of an asymptotic series in powers of 1/ln N which is formally identical to those found for Euclidean and (some class of) deterministic fractal lattices. The agreement of the asymptotic expressions with simulation results for the two-dimensional percolation aggregate is good when the boundary is defined in terms of the chemical distance. The agreement worsens slightly when the Euclidean distance is used.Comment: 8 pages including 9 figure

Order statistics for d-dimensional diffusion processes

We present results for the ordered sequence of first passage times of arrival of N random walkers at a boundary in Euclidean spaces of d dimensions

Pair correlation function of short-ranged square-well fluids

We have performed extensive Monte Carlo simulations in the canonical (NVT) ensemble of the pair correlation function for square-well fluids with well widths $\lambda-1$ ranging from 0.1 to 1.0, in units of the diameter $\sigma$ of the particles. For each one of these widths, several densities $\rho$ and temperatures $T$ in the ranges $0.1\leq\rho\sigma^3\leq 0.8$ and $T_c(\lambda)\lesssim T\lesssim 3T_c(\lambda)$, where $T_c(\lambda)$ is the critical temperature, have been considered. The simulation data are used to examine the performance of two analytical theories in predicting the structure of these fluids: the perturbation theory proposed by Tang and Lu [Y. Tang and B. C.-Y. Lu, J. Chem. Phys. {\bf 100}, 3079, 6665 (1994)] and the non-perturbative model proposed by two of us [S. B. Yuste and A. Santos, J. Chem. Phys. {\bf 101}, 2355 (1994)]. It is observed that both theories complement each other, as the latter theory works well for short ranges and/or moderate densities, while the former theory does for long ranges and high densities.Comment: 10 pages, 10 figure

Serum amyloid P aids complement-mediated immunity to Streptococcus pneumoniae

The physiological functions of the acute phase protein serum amyloid P (SAP) component are not well defined, although they are likely to be important, as no natural state of SAP deficiency has been reported. We have investigated the role of SAP for innate immunity to the important human pathogen Streptococcus pneumoniae. Using flow cytometry assays, we show that SAP binds to S. pneumoniae, increases classical pathway–dependent deposition of complement on the bacteria, and improves the efficiency of phagocytosis. As a consequence, in mouse models of infection, mice genetically engineered to be SAP-deficient had an impaired early inflammatory response to S. pneumoniae pneumonia and were unable to control bacterial replication, leading to the rapid development of fatal infection. Complement deposition, phagocytosis, and control of S. pneumoniae pneumonia were all improved by complementation with human SAP. These results demonstrate a novel and physiologically significant role for SAP for complement-mediated immunity against an important bacterial pathogen, and provide further evidence for the importance of the classical complement pathway for innate immunity

Coagulation reaction in low dimensions: Revisiting subdiffusive A+A reactions in one dimension

We present a theory for the coagulation reaction A+A -> A for particles moving subdiffusively in one dimension. Our theory is tested against numerical simulations of the concentration of $A$ particles as a function of time (``anomalous kinetics'') and of the interparticle distribution function as a function of interparticle distance and time. We find that the theory captures the correct behavior asymptotically and also at early times, and that it does so whether the particles are nearly diffusive or very subdiffusive. We find that, as in the normal diffusion problem, an interparticle gap responsible for the anomalous kinetics develops and grows with time. This corrects an earlier claim to the contrary on our part.Comment: The previous version was corrupted - some figures misplaced, some strange words that did not belong. Otherwise identica

Mean Field Model of Coagulation and Annihilation Reactions in a Medium of Quenched Traps: Subdiffusion

We present a mean field model for coagulation ($A+A\to A$) and annihilation ($A+A\to 0$) reactions on lattices of traps with a distribution of depths reflected in a distribution of mean escape times. The escape time from each trap is exponentially distributed about the mean for that trap, and the distribution of mean escape times is a power law. Even in the absence of reactions, the distribution of particles over sites changes with time as particles are caught in ever deeper traps, that is, the distribution exhibits aging. Our main goal is to explore whether the reactions lead to further (time dependent) changes in this distribution.Comment: 9 pages, 3 figure

How `sticky' are short-range square-well fluids?

The aim of this work is to investigate to what extent the structural properties of a short-range square-well (SW) fluid of range $\lambda$ at a given packing fraction and reduced temperature can be represented by those of a sticky-hard-sphere (SHS) fluid at the same packing fraction and an effective stickiness parameter $\tau$. Such an equivalence cannot hold for the radial distribution function since this function has a delta singularity at contact in the SHS case, while it has a jump discontinuity at $r=\lambda$ in the SW case. Therefore, the equivalence is explored with the cavity function $y(r)$. Optimization of the agreement between y_{\sw} and y_{\shs} to first order in density suggests the choice for $\tau$. We have performed Monte Carlo (MC) simulations of the SW fluid for $\lambda=1.05$, 1.02, and 1.01 at several densities and temperatures $T^*$ such that $\tau=0.13$, 0.2, and 0.5. The resulting cavity functions have been compared with MC data of SHS fluids obtained by Miller and Frenkel [J. Phys: Cond. Matter 16, S4901 (2004)]. Although, at given values of $\eta$ and $\tau$, some local discrepancies between y_{\sw} and y_{\shs} exist (especially for $\lambda=1.05$), the SW data converge smoothly toward the SHS values as $\lambda-1$ decreases. The approximate mapping y_{\sw}\to y_{\shs} is exploited to estimate the internal energy and structure factor of the SW fluid from those of the SHS fluid. Taking for y_{\shs} the solution of the Percus--Yevick equation as well as the rational-function approximation, the radial distribution function $g(r)$ of the SW fluid is theoretically estimated and a good agreement with our MC simulations is found. Finally, a similar study is carried out for short-range SW fluid mixtures.Comment: 14 pages, including 3 tables and 14 figures; v2: typo in Eq. (5.1) corrected, Fig. 14 redone, to be published in JC

Contact values of the radial distribution functions of additive hard-sphere mixtures in d dimensions: A new proposal

The contact values $g_{ij}(\sigma_{ij})$ of the radial distribution functions of a $d$-dimensional mixture of (additive) hard spheres are considered. A `universality' assumption is put forward, according to which $g_{ij}(\sigma_{ij})=G(\eta, z_{ij})$, where $G$ is a common function for all the mixtures of the same dimensionality, regardless of the number of components, $\eta$ is the packing fraction of the mixture, and $z_{ij}$ is a dimensionless parameter that depends on the size distribution and the diameters of spheres $i$ and $j$. For $d=3$, this universality assumption holds for the contact values of the Percus--Yevick approximation, the Scaled Particle Theory, and, consequently, the Boublik--Grundke--Henderson--Lee--Levesque approximation. Known exact consistency conditions are used to express $G(\eta,0)$, $G(\eta,1)$, and $G(\eta,2)$ in terms of the radial distribution at contact of the one-component system. Two specific proposals consistent with the above conditions (a quadratic form and a rational form) are made for the $z$-dependence of $G(\eta,z)$. For one-dimensional systems, the proposals for the contact values reduce to the exact result. Good agreement between the predictions of the proposals and available numerical results is found for $d=2$, 3, 4, and 5.Comment: 10 pages, 11 figures; Figure 1 changed; Figure 5 is new; New references added; accepted for publication in J. Chem. Phy