We investigate numerically pseudo-first-order irreversible bimolecular reactions of the type A + B -> B between hard spheres undergoing event-driven Brownian dynamics. We study the encounter rate and the survival probability of A particles as functions of the packing fraction phi in the trapping (a single particle diffusing among static non-overlapping traps) and target (many traps diffusing in the presence of a single static target particle) settings, as well as in the case of diffusing traps and particles (full mobility). We show that, since inertial effects are accounted for in our simulation protocol, the standard Smoluchowski theory of coagulation of non-interacting colloids is recovered only at times greater than a characteristic time Delta t, marking the transition from the under-damped to the over-damped regime. We show that the survival probability S(t) decays exponentially during this first stage, with a rate 1/tau(0) proportional to phi. Furthermore, we work out a simple analytical expression that is able to capture to an excellent extent the numerical results for t < Delta t at low and intermediate densities. Moreover, we demonstrate that the time constant of the asymptotic exponential decay of S(t) for diffusing traps and particles is k(S)(-1), where k(S) = 4 pi(D-A +/- D-B)R rho is the Smoluchowski rate. Detailed analyses of the effective decay exponent beta = dinverted right perpendicularlog(-log S(t))inverted left perpendicular/d(log t) and of the steady-state encounter rate reveal that the full mobility and trapping problem are characterized by very similar kinetics, rather different from the target problem. Our results do not allow one to ascertain whether the prediction S(t) proportional to exp(-at(3/2))(a = const.) as t -> infinity for the trapping problem in 3D is indeed recovered. In fact, at high density, S(t) is dominated by short encounter times, which makes it exceedingly hard to record the events corresponding to the exploration of large, trap-free regions. As a consequence, at high densities the steady-state rate simply tends to 1/tau(0). Finally, we work out an analytical formula for the rate that shows a remarkable agreement with the numerics up phi = 0.4.

Irreversible bimolecular reactions with inertia: from the trapping to the target setting at finite densities / Francesco, Piazza; Giuseppe, Foffi; DE MICHELE, Cristiano. - In: JOURNAL OF PHYSICS. CONDENSED MATTER. - ISSN 0953-8984. - 25:24(2013). [10.1088/0953-8984/25/24/245101]

Irreversible bimolecular reactions with inertia: from the trapping to the target setting at finite densities

DE MICHELE, CRISTIANO
2013

Abstract

We investigate numerically pseudo-first-order irreversible bimolecular reactions of the type A + B -> B between hard spheres undergoing event-driven Brownian dynamics. We study the encounter rate and the survival probability of A particles as functions of the packing fraction phi in the trapping (a single particle diffusing among static non-overlapping traps) and target (many traps diffusing in the presence of a single static target particle) settings, as well as in the case of diffusing traps and particles (full mobility). We show that, since inertial effects are accounted for in our simulation protocol, the standard Smoluchowski theory of coagulation of non-interacting colloids is recovered only at times greater than a characteristic time Delta t, marking the transition from the under-damped to the over-damped regime. We show that the survival probability S(t) decays exponentially during this first stage, with a rate 1/tau(0) proportional to phi. Furthermore, we work out a simple analytical expression that is able to capture to an excellent extent the numerical results for t < Delta t at low and intermediate densities. Moreover, we demonstrate that the time constant of the asymptotic exponential decay of S(t) for diffusing traps and particles is k(S)(-1), where k(S) = 4 pi(D-A +/- D-B)R rho is the Smoluchowski rate. Detailed analyses of the effective decay exponent beta = dinverted right perpendicularlog(-log S(t))inverted left perpendicular/d(log t) and of the steady-state encounter rate reveal that the full mobility and trapping problem are characterized by very similar kinetics, rather different from the target problem. Our results do not allow one to ascertain whether the prediction S(t) proportional to exp(-at(3/2))(a = const.) as t -> infinity for the trapping problem in 3D is indeed recovered. In fact, at high density, S(t) is dominated by short encounter times, which makes it exceedingly hard to record the events corresponding to the exploration of large, trap-free regions. As a consequence, at high densities the steady-state rate simply tends to 1/tau(0). Finally, we work out an analytical formula for the rate that shows a remarkable agreement with the numerics up phi = 0.4.
2013
bimolecular reaftions; brownian dynamics; computer simulation; theory
01 Pubblicazione su rivista::01a Articolo in rivista
Irreversible bimolecular reactions with inertia: from the trapping to the target setting at finite densities / Francesco, Piazza; Giuseppe, Foffi; DE MICHELE, Cristiano. - In: JOURNAL OF PHYSICS. CONDENSED MATTER. - ISSN 0953-8984. - 25:24(2013). [10.1088/0953-8984/25/24/245101]
File allegati a questo prodotto
File Dimensione Formato  
Piazza_Irreversible_2013.pdf

solo gestori archivio

Tipologia: Versione editoriale (versione pubblicata con il layout dell'editore)
Licenza: Tutti i diritti riservati (All rights reserved)
Dimensione 1.02 MB
Formato Adobe PDF
1.02 MB Adobe PDF   Contatta l'autore

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11573/517818
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 7
  • ???jsp.display-item.citation.isi??? 7
social impact