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

Number of distinct sites visited by N random walkers on a Euclidean lattice

The evaluation of the average number S_N(t) of distinct sites visited up to time t by N independent random walkers all starting from the same origin on an Euclidean lattice is addressed. We find that, for the nontrivial time regime and for large N, S_N(t) \approx \hat S_N(t) (1-\Delta), where \hat S_N(t) is the volume of a hypersphere of radius (4Dt \ln N)^{1/2}, \Delta={1/2}\sum_{n=1}^\infty \ln^{-n} N \sum_{m=0}^n s_m^{(n)} \ln^{m} \ln N, d is the dimension of the lattice, and the coefficients s_m^{(n)} depend on the dimension and time. The first three terms of these series are calculated explicitly and the resulting expressions are compared with other approximations and with simulation results for dimensions 1, 2, and 3. Some implications of these results on the geometry of the set of visited sites are discussed.Comment: 15 pages (RevTex), 4 figures (eps); to appear in Phys. Rev.

Structure of ternary additive hard-sphere fluid mixtures

Monte Carlo simulations on the structural properties of ternary fluid mixtures of additive hard spheres are reported. The results are compared with those obtained from a recent analytical approximation [S. B. Yuste, A. Santos, and M. Lopez de Haro, J. Chem. Phys. 108, 3683 (1998)] to the radial distribution functions of hard-sphere mixtures and with the results derived from the solution of the Ornstein-Zernike integral equation with both the Martynov-Sarkisov and the Percus-Yevick closures. Very good agreement between the results of the first two approaches and simulation is observed, with a noticeable improvement over the Percus-Yevick predictions especially near contact.Comment: 11 pages, including 8 figures; A minor change; accepted for publication in PR

Some exact results for the trapping of subdiffusive particles in one dimension

We study a generalization of the standard trapping problem of random walk theory in which particles move subdiffusively on a one-dimensional lattice. We consider the cases in which the lattice is filled with a one-sided and a two-sided random distribution of static absorbing traps with concentration c. The survival probability Phi(t) that the random walker is not trapped by time t is obtained exactly in both versions of the problem through a fractional diffusion approach. Comparison with simulation results is madeComment: 15 pages, 2 figure

A model for the atomic-scale structure of a dense, nonequilibrium fluid: the homogeneous cooling state of granular fluids

It is shown that the equilibrium Generalized Mean Spherical Model of fluid structure may be extended to nonequilibrium states with equation of state information used in equilibrium replaced by an exact condition on the two-body distribution function. The model is applied to the homogeneous cooling state of granular fluids and upon comparison to molecular dynamics simulations is found to provide an accurate picture of the pair distribution function.Comment: 29 pages, 11 figures Revision corrects formatting of the figure

Adsorption of Reactive Particles on a Random Catalytic Chain: An Exact Solution

We study equilibrium properties of a catalytically-activated annihilation $A + A \to 0$ reaction taking place on a one-dimensional chain of length $N$ ($N \to \infty$) in which some segments (placed at random, with mean concentration $p$) possess special, catalytic properties. Annihilation reaction takes place, as soon as any two $A$ particles land onto two vacant sites at the extremities of the catalytic segment, or when any $A$ particle lands onto a vacant site on a catalytic segment while the site at the other extremity of this segment is already occupied by another $A$ particle. Non-catalytic segments are inert with respect to reaction and here two adsorbed $A$ particles harmlessly coexist. For both "annealed" and "quenched" disorder in placement of the catalytic segments, we calculate exactly the disorder-average pressure per site. Explicit asymptotic formulae for the particle mean density and the compressibility are also presented.Comment: AMSTeX, 27 pages + 4 figure

Reaction Front in an A+B -> C Reaction-Subdiffusion Process

We study the reaction front for the process A+B -> C in which the reagents move subdiffusively. Our theoretical description is based on a fractional reaction-subdiffusion equation in which both the motion and the reaction terms are affected by the subdiffusive character of the process. We design numerical simulations to check our theoretical results, describing the simulations in some detail because the rules necessarily differ in important respects from those used in diffusive processes. Comparisons between theory and simulations are on the whole favorable, with the most difficult quantities to capture being those that involve very small numbers of particles. In particular, we analyze the total number of product particles, the width of the depletion zone, the production profile of product and its width, as well as the reactant concentrations at the center of the reaction zone, all as a function of time. We also analyze the shape of the product profile as a function of time, in particular its unusual behavior at the center of the reaction zone

Alignment of the CMS silicon tracker during commissioning with cosmic rays

This is the Pre-print version of the Article. The official published version of the Paper can be accessed from the link below - Copyright @ 2010 IOPThe CMS silicon tracker, consisting of 1440 silicon pixel and 15 148 silicon strip detector modules, has been aligned using more than three million cosmic ray charged particles, with additional information from optical surveys. The positions of the modules were determined with respect to cosmic ray trajectories to an average precision of 3–4 microns RMS in the barrel and 3–14 microns RMS in the endcap in the most sensitive coordinate. The results have been validated by several studies, including laser beam cross-checks, track fit self-consistency, track residuals in overlapping module regions, and track parameter resolution, and are compared with predictions obtained from simulation. Correlated systematic effects have been investigated. The track parameter resolutions obtained with this alignment are close to the design performance.This work is supported by FMSR (Austria); FNRS and FWO (Belgium); CNPq, CAPES, FAPERJ, and FAPESP (Brazil); MES (Bulgaria); CERN; CAS, MoST, and NSFC (China); COLCIENCIAS (Colombia); MSES (Croatia); RPF (Cyprus); Academy of Sciences and NICPB (Estonia); Academy of Finland, ME, and HIP (Finland); CEA and CNRS/IN2P3 (France); BMBF, DFG, and HGF (Germany); GSRT (Greece); OTKA and NKTH (Hungary); DAE and DST (India); IPM (Iran); SFI (Ireland); INFN (Italy); NRF (Korea); LAS (Lithuania); CINVESTAV, CONACYT, SEP, and UASLP-FAI (Mexico); PAEC (Pakistan); SCSR (Poland); FCT (Portugal); JINR (Armenia, Belarus, Georgia, Ukraine, Uzbekistan); MST and MAE (Russia); MSTDS (Serbia); MICINN and CPAN (Spain); Swiss Funding Agencies (Switzerland); NSC (Taipei); TUBITAK and TAEK (Turkey); STFC (United Kingdom); DOE and NSF (USA)

Commissioning and performance of the CMS pixel tracker with cosmic ray muons

This is the Pre-print version of the Article. The official published verion of the Paper can be accessed from the link below - Copyright @ 2010 IOPThe pixel detector of the Compact Muon Solenoid experiment consists of three barrel layers and two disks for each endcap. The detector was installed in summer 2008, commissioned with charge injections, and operated in the 3.8 T magnetic field during cosmic ray data taking. This paper reports on the first running experience and presents results on the pixel tracker performance, which are found to be in line with the design specifications of this detector. The transverse impact parameter resolution measured in a sample of high momentum muons is 18 microns.This work is supported by FMSR (Austria); FNRS and FWO (Belgium); CNPq, CAPES, FAPERJ, and FAPESP (Brazil); MES (Bulgaria); CERN; CAS, MoST, and NSFC (China); COLCIENCIAS (Colombia); MSES (Croatia); RPF (Cyprus); Academy of Sciences and NICPB (Estonia); Academy of Finland, ME, and HIP (Finland); CEA and CNRS/IN2P3 (France); BMBF, DFG, and HGF (Germany); GSRT (Greece); OTKA and NKTH (Hungary); DAE and DST (India); IPM (Iran); SFI (Ireland); INFN (Italy); NRF (Korea); LAS (Lithuania); CINVESTAV, CONACYT, SEP, and UASLP-FAI (Mexico); PAEC (Pakistan); SCSR (Poland); FCT (Portugal); JINR (Armenia, Belarus, Georgia, Ukraine, Uzbekistan); MST and MAE (Russia); MSTDS (Serbia); MICINN and CPAN (Spain); Swiss Funding Agencies (Switzerland); NSC (Taipei); TUBITAK and TAEK (Turkey); STFC (United Kingdom); DOE and NSF (USA)

Performance of the CMS drift-tube chamber local trigger with cosmic rays

The performance of the Local Trigger based on the drift-tube system of the CMS experiment has been studied using muons from cosmic ray events collected during the commissioning of the detector in 2008. The properties of the system are extensively tested and compared with the simulation. The effect of the random arrival time of the cosmic rays on the trigger performance is reported, and the results are compared with the design expectations for proton-proton collisions and with previous measurements obtained with muon beams