Abstract:Modulating ion transport through nanoporous membranes is critical to many important chemical and biological separation processes. The corresponding transport timescales, however, are often too long to capture accurately using conventional molecular dynamics (\textsc{Md}). Recently, path sampling techniques such as forward-flux sampling (\textsc{Ffs}) have emerged as attractive alternatives for efficiently and accurately estimating arbitrarily long ionic passage times. Here, we use non-equilibrium \textsc{Md} and \textsc{Ffs} to explore how the kinetics and mechanism of pressure-driven chloride transport through a nanoporous graphitic membrane are affected by its lateral dimensions. We not only find ionic passage times and free energy barriers to decrease dramatically upon increasing the membrane surface area, but also observe an abrupt and discontinuous change in the locus of the transition state. These strong finite size effects arise due to the cumulative effect of the periodic images of the leading ion entering the pore on the distribution of the induced excess charge at the membrane surface in the feed. By assuming that the feed is an ideal conductor, we analytically derive a finite size correction term that can be computed from the information obtained from a single simulation and successfully use it to obtain corrected free energy profiles with no dependence on system size. We then estimate ionic passage times in the thermodynamic limit by assuming an Eyring-type dependence of rates on barriers with a size-independent prefactor. This approach constitutes a universal framework for removing finite size artifacts in molecular simulations of ion transport through nanoporous membranes and biological channel proteins.

Usefulness of Higher Order of System-Size Correction for Diffusion Coefficients of Macromolecules: A Molecular Dynamics Study

Correlation of Infinite-Dilution Diffusion Coefficients in Supercritical Fluids

SIMULATION STUDIES ON THE DIFFUSION COEFFICIENT OF PROTEIN-LIKE CHAINS

Application of Molecular Dynamics Simulations in Molecular Property Prediction II: Diffusion Coefficient.

System Size Dependence of the Diffusion Coefficients in MD Simulations: A Simple Correction Formula for Pure Dense Fluids

Optimal estimates of diffusion coefficients from molecular dynamics simulations

On the system size dependence of the diffusion coefficients in MD simulations: A simple correction formula for pure dense fluids

Accurate Estimation of Diffusion Coefficients and their Uncertainties from Computer Simulation

Molecular Dynamics Simulation of Self-Diffusion Coefficient and Its Relation with Temperature Using Simple Lennard-Jones Potential

Analytical And Numerical Approximation of Effective Diffusivities in The Cytoplasm of Biological Cells

Enhanced solvation force extrapolation for speeding up molecular dynamics simulations of complex biochemical liquids

Ion sizes and finite-size corrections for ionic-solvation free energies

Ideal Conductor Model: An analytical finite-size correction for non-equilibrium molecular dynamics simulations of ion transport through nanoporous membranes

Computing Accurate True Self-Diffusion Coefficients and Shear Viscosities Using the OrthoBoXY-Approach

A simple and effective solution to the constrained QM/MM simulations

Momentum removal to obtain the position‐dependent diffusion constant in constrained molecular dynamics simulation

OrthoBoXY: A Simple Way to Compute True Self-Diffusion Coefficients from MD Simulations with Periodic Boundary Conditions Without Prior Knowledge of the Viscosity

Analysis of the Distribution and Influencing Factors of Diffusion Coefficient Model Parameters Based on Molecular Dynamics Simulations

MD2D: a python module for accurate determination of diffusion coefficient from molecular dynamics

Adaptive Molecular Resolution Approach in Hamiltonian Form: An Asymptotic Analysis

Statistical variances of diffusional properties from ab initio molecular dynamics simulations