A transition state theory for calculating hopping times and diffusion in highly confined fluids

The effect of soft repulsive interactions on the diffusion of particles in quasi-one-dimensional channels: A hopping time approach

Fluids confined to quasi-one-dimensional channels exhibit a dynamic crossover from single file diffusion to normal diffusion as the channel becomes wide enough for particles to hop past each other. In the crossover regime, where hopping events are rare, the diffusion coefficient in the long time limit can be related to a hopping time that measures the average time it takes for a particle to escape the local cage formed by its neighbors. In this work, we show that a transition state theory (TST) that calculates the free energy barrier for two particles attempting to pass each other in the small system isobaric ensemble is able to quantitatively predict the hopping time in a system of two-dimensional soft repulsive disks [U(r_ij)=(σ/r_ij)^α] confined to a hard walled channel over a range of channel radii and degrees of particle softness measured in terms of 1/α. The free energy barrier exhibits a maximum at intermediate values of α that moves to smaller values of 1/α (harder particles) as the channel becomes narrower. However, the presence of the maximum is only observed in the hopping times for wide channels because the interaction potential dependence of the kinetic prefactor plays an increasingly important role for narrower channels. We also begin to explore how our TST approach can be used to optimize and control dynamics in confined quasi-one-dimensional fluids.

Two-tag correlations and nonequilibrium fluctuation-response relation in ageing single-file diffusion

Spatiotemporally correlated motions of interacting Brownian particles, confined in a narrow channel of infinite length, are studied in terms of statistical quantities involving two particles. A theoretical framework that allows analytical calculation of two-tag correlations is presented on the basis of the Dean-Kawasaki equation describing density fluctuations in colloidal systems. In the equilibrium case, the time-dependent Einstein relation holds between the two-tag displacement correlation and the response function corresponding to it, which is a manifestation of the fluctuation-dissipation theorem for the correlation of density fluctuations. While the standard procedure of closure approximation for nonlinear density fluctuations is known to be obstructed by inconsistency with the fluctuation-dissipation theorem, this difficulty is naturally avoided by switching from the standard Fourier representation of the density field to the label-based Fourier representation of the vacancy field. In the case of ageing dynamics started from equidistant lattice configuration, the time-dependent Einstein relation is violated, as the two-tag correlation depends on the waiting time for equilibration while the response function is not sensitive to it. Within linear approximation, however, there is a simple relation between the density (or vacancy) fluctuations and the corresponding response function, which is valid even if the system is out of equilibrium. This non-equilibrium fluctuation-response relation can be extended to the case of nonlinear fluctuations by means of closure approximation for the vacancy field.

Collective Motion of Repulsive Brownian Particles in Single-File Diffusion with and without Overtaking

Subdiffusion is commonly observed in liquids with high density or in restricted geometries, as the particles are constantly pushed back by their neighbors. Since this “cage effect” emerges from many-body dynamics involving spatiotemporally correlated motions, the slow diffusion should be understood not simply as a one-body problem but as a part of collective dynamics, described in terms of space–time correlations. Such collective dynamics are illustrated here by calculations of the two-particle displacement correlation in a system of repulsive Brownian particles confined in a (quasi-)one-dimensional channel, whose subdiffusive behavior is known as the single-file diffusion (SFD). The analytical calculation is formulated in terms of the Lagrangian correlation of density fluctuations. In addition, numerical solutions to the Langevin equation with large but finite interaction potential are studied to clarify the effect of overtaking. In the limiting case of the ideal SFD without overtaking, correlated motion with a diffusively growing length scale is observed. By allowing the particles to overtake each other, the short-range correlation is destroyed, but the long-range weak correlation remains almost intact. These results describe nested space–time structure of cages, whereby smaller cages are enclosed in larger cages with longer lifetimes.

Diffusion in quasi-one-dimensional channels: A small system n, p, T, transition state theory for hopping times

Particles confined to a single file, in a narrow quasi-one-dimensional channel, exhibit a dynamic crossover from single file diffusion to Fickian diffusion as the channel radius increases and the particles begin to pass each other. The long time diffusion coefficient for a system in the crossover regime can be described in terms of a hopping time, which measures the time it takes for a particle to escape the cage formed by its neighbours. In this paper, we develop a transition state theory approach to the calculation of the hopping time, using the small system isobaric-isothermal ensemble to rigorously account for the volume fluctuations associated with the size of the cage. We also describe a Monte Carlo simulation scheme that can be used to calculate the free energy barrier for particle hopping. The theory and simulation method correctly predict the hopping times for a two-dimensional confined ideal gas system and a system of confined hard discs over a range of channel radii, but the method breaks down for wide channels in the hard discs' case, underestimating the height of the hopping barrier due to the neglect of interactions between the small system and its surroundings.

Transition state theory for solvated reactions beyond recrossing-free dividing surfaces

The accuracy of rate constants calculated using transition state theory depends crucially on the correct identification of a recrossing-free dividing surface. We show here that it is possible to define such optimal dividing surface in systems with non-Markovian friction. However, a more direct approach to rate calculation is based on invariant manifolds and avoids the use of a dividing surface altogether, Using that method we obtain an explicit expression for the rate of crossing an anharmonic potential barrier. The excellent performance of our method is illustrated with an application to a realistic model for LiNC⇌LiCN isomerization.

Insights from Single-File Diffusion into Cooperativity in Higher Dimensions

Diffusion in colloidal suspensions can be very slow due to the cage effect, which confines each particle within a short radius on one hand, and involves large-scale cooperative motions on the other. In search of insight into this cooperativity, here the authors develop a formalism to calculate the displacement correlation in colloidal systems, mainly in the two-dimensional (2D) case. To clarify the idea for it, studies are reviewed on cooperativity among the particles in the one-dimensional (1D) case, i.e. the single-file diffusion (SFD). As an improvement over the celebrated formula by Alexander and Pincus on the mean-square displacement (MSD) in SFD, it is shown that the displacement correlation in SFD can be calculated from Lagrangian correlation of the particle interval in the one-dimensional case, and also that the formula can be extended to higher dimensions. The improved formula becomes exact for large systems. By combining the formula with a nonlinear theory for correlation, a correction to the asymptotic law for the MSD in SFD is obtained. In the 2D case, the linear theory gives description of vortical cooperative motion.

Equivalence of Brownian dynamics and dynamic Monte Carlo simulations in multicomponent colloidal suspensions

We propose a simple but powerful theoretical framework to quantitatively compare Brownian dynamics (BD) and dynamic Monte Carlo (DMC) simulations of multicomponent colloidal suspensions. By extending our previous study focusing on monodisperse systems of rodlike colloids, here we generalize the formalism described there to multicomponent colloidal mixtures and validate it by investigating the dynamics in isotropic and liquid crystalline phases containing spherical and rodlike particles. In order to investigate the dynamics of multicomponent colloidal systems by DMC simulations, it is key to determine the elementary time step of each species and establish a unique timescale. This is crucial to consistently study the dynamics of colloidal particles with different geometry. By analyzing the mean-square displacement, the orientation autocorrelation functions, and the self part of the van Hove correlation functions, we show that DMC simulation is a very convenient and reliable technique to describe the stochastic dynamics of any multicomponent colloidal system. Our theoretical formalism can be easily extended to any colloidal system containing size and/or shape polydisperse particles.

Relationship between single-file diffusion of mixed and pure gases in dipeptide nanochannels by high field diffusion NMR

Under single-file confinement, the relationship between diffusion rates of mixed and pure gases is studied experimentally for the first time and observed to differ from that for normal diffusion.