Steady-state properties of single-file systems with conversion

Generalized hydrodynamic analysis of transport through a finite open nanopore for two-component single-file systems

Single-file diffusion (SFD) in finite open nanopores is characterized by nonzero spatially varying tracer diffusion coefficients within a generalized hydrodynamic description. This contrasts with infinite SFD systems where tracer diffusivity vanishes. In standard tracer counterpermeation (TCP) analysis, two reservoirs, each containing a different species, are connected to opposite ends of a finite pore. We implement an extended TCP analysis to allow the two reservoirs to contain slightly different mixtures of the two species. Then, determination of diffusion fluxes through the pore allows extraction of diffusion coefficients for near-constant partial concentrations of the two species. This analysis is applied for a lattice-gas model describing two-component SFD through a finite linear pore represented by a one-dimensional array of cells. Two types of particles, A and B, can hop only to adjacent empty cells with generally different rates, h_{A} and h_{B}. Particles are noninteracting other than exclusion of multiple cell occupancy. Results reveal generalized hydrodynamic tracer diffusion coefficients which adopt small values inversely proportional to pore length in the pore center, but which are strongly enhanced near pore openings.

Catalytic conversion in nanoporous materials: Concentration oscillations and spatial correlations due to inhibited transport and intermolecular interactions

We show that steady-state catalytic conversion in nanoporous materials can occur in a quasi-counter-diffusion mode with the reactant (product) concentration strongly decaying (growing) into the pore, but also with oscillations in the total concentration. These oscillations reflect the response of the fluid to the transition from an extended to a confined environment near the pore opening. We focus on the regime of strongly inhibited transport in narrow pores corresponding to single-file diffusion. Here, limited penetration of the reactant into the pores and the associated low reaction yield is impacted by strong spatial correlations induced by both reaction (non-equilibrium correlations) and also by intermolecular interactions (thermodynamic correlations). We develop a generalized hydrodynamic formulation to effectively describe inhibited transport accounting for the effect of these correlations, and incorporate this description of transport into appropriate reaction-diffusion equations. These equations accurately describe both shorter-range concentration oscillations near the pore opening and the longer-range mesoscale variation of concentration profiles in the pore (and thus also describe reaction yield). Success of the analytic theory is validated by comparison with a precise kinetic Monte Carlo simulation of an appropriate molecular-level stochastic reaction-diffusion model. This work elucidates unconventional chemical kinetics in interacting confined systems.

Catalytic conversion reactions in nanoporous systems with concentration-dependent selectivity: Statistical mechanical modeling

Statistical mechanical modeling is developed to describe a catalytic conversion reaction A→B^{c} or B^{t} with concentration-dependent selectivity of the products, B^{c} or B^{t}, where reaction occurs inside catalytic particles traversed by narrow linear nanopores. The associated restricted diffusive transport, which in the extreme case is described by single-file diffusion, naturally induces strong concentration gradients. Furthermore, by comparing kinetic Monte Carlo simulation results with analytic treatments, selectivity is shown to be impacted by strong spatial correlations induced by restricted diffusivity in the presence of reaction and also by a subtle clustering of reactants, A.

Langevin and fOkker-Planck analyses of inhibited molecular passing processes controlling transport and reactivity in nanoporous materials

Inhibited passing of reactant and product molecules within the linear pores of nanoporous catalytic materials strongly reduces reactivity. The dependence of the passing propensity P on pore radius R is analyzed utilizing Langevin dynamics to account for solvent effects. We find that P ∼ (R-R(c))(σ), where passing is sterically blocked for R≤R(c), with σ below the transition state theory value. Deeper insight comes from analysis of the corresponding high-dimensional Fokker-Planck equation, which facilitates an effective small-P approximation, and dimensional reduction enabling utilization of conformal mapping ideas. We analyze passing for spherical molecules and also assess the effect of rotational degrees of freedom for elongated molecules.

Generalized hydrodynamic treatment of the interplay between restricted transport and catalytic reactions in nanoporous materials

Behavior of catalytic reactions in narrow pores is controlled by a delicate interplay between fluctuations in adsorption-desorption at pore openings, restricted diffusion, and reaction. This behavior is captured by a generalized hydrodynamic formulation of appropriate reaction-diffusion equations (RDE). These RDE incorporate an unconventional description of chemical diffusion in mixed-component quasi-single-file systems based on a refined picture of tracer diffusion for finite-length pores. The RDE elucidate the nonexponential decay of the steady-state reactant concentration into the pore and the non-mean-field scaling of the reactant penetration depth.

Catalytic conversion reactions mediated by single-file diffusion in linear nanopores: hydrodynamic versus stochastic behavior

We analyze the spatiotemporal behavior of species concentrations in a diffusion-mediated conversion reaction which occurs at catalytic sites within linear pores of nanometer diameter. Diffusion within the pores is subject to a strict single-file (no passing) constraint. Both transient and steady-state behavior is precisely characterized by kinetic Monte Carlo simulations of a spatially discrete lattice-gas model for this reaction-diffusion process considering various distributions of catalytic sites. Exact hierarchical master equations can also be developed for this model. Their analysis, after application of mean-field type truncation approximations, produces discrete reaction-diffusion type equations (mf-RDE). For slowly varying concentrations, we further develop coarse-grained continuum hydrodynamic reaction-diffusion equations (h-RDE) incorporating a precise treatment of single-file diffusion in this multispecies system. The h-RDE successfully describe nontrivial aspects of transient behavior, in contrast to the mf-RDE, and also correctly capture unreactive steady-state behavior in the pore interior. However, steady-state reactivity, which is localized near the pore ends when those regions are catalytic, is controlled by fluctuations not incorporated into the hydrodynamic treatment. The mf-RDE partly capture these fluctuation effects, but cannot describe scaling behavior of the reactivity.

Polymer length distributions for catalytic polymerization within mesoporous materials: non-Markovian behavior associated with partial extrusion

We analyze a model for polymerization at catalytic sites distributed within parallel linear pores of a mesoporous material. Polymerization occurs primarily by reaction of monomers diffusing into the pores with the ends of polymers near the pore openings. Monomers and polymers undergo single-file diffusion within the pores. Model behavior, including the polymer length distribution, is determined by kinetic Monte Carlo simulation of a suitable atomistic-level lattice model. While the polymers remain within the pore, their length distribution during growth can be described qualitatively by a Markovian rate equation treatment. However, once they become partially extruded, the distribution is shown to exhibit non-Markovian scaling behavior. This feature is attributed to the long-tail in the "return-time distribution" for the protruding end of the partially extruded polymer to return to the pore, such return being necessary for further reaction and growth. The detailed form of the scaled length distribution is elucidated by application of continuous-time random walk theory.

Statistical mechanical modeling of catalytic polymerization within surface-functionalized mesoporous materials

A discrete lattice model is developed to describe diffusion-mediated polymerization occurring within mesopores, where reaction is enhanced at catalytic sites distributed within the interior of the pores. Diffusive transport of monomers and polymers is one-dimensional, diffusion coefficients for the latter decreasing with polymer length. Kinetic Monte Carlo simulation is utilized to analyze model behavior focusing on a "clogging" regime, where the amount of polymer within the pores grows. We characterize the evolution of the overall and mean length of polymers, the mean number of polymers, as well as the polymer spatial and length distributions.

Infinitely fast diffusion in single-file systems

We have used dynamic Monte Carlo(DMC) methods and analytical techniques to analyze single-file systems for which diffusion is infinitely fast. We have simplified the master equation removing the fast reactions, and we have introduced a DMC algorithm for infinitely fast diffusion. The DMC method for fast diffusion give similar results as the standard DMC with high diffusion rates. We have investigated the influence of characteristic parameters, such as pipe length, adsorption, desorption, and conversion rate constants on the steady-state properties of single-file systems with a reaction, looking at cases when all the sites are reactive and when only some of them are reactive. We find that the effect of fast diffusion on single-file properties of the system is absent even when diffusion is infinitely fast. Diffusion is not important in these systems. Smaller systems are less reactive and the occupancy profiles for infinitely long systems show an exponential behavior.

Analytical results for the reactivity of a single-file system

We derive analytical expressions for the reactivity of a single-file system with fast diffusion and particles entering and leaving the system at one end. If the conversion reaction is fast, then the reactivity depends only very weakly on the system size, and the conversion is about 100%. If the reaction is slow, then the reactivity becomes proportional to the system size, the loading, and the reaction rate constant. If the system size increases the reactivity goes to the geometric mean of the reaction rate constant and the rate of particles entering and leaving the system. For large systems, the number of unconverted particles decreases exponentially with distance from the open end.