We report a multiscale modeling study for charged cylindrical nanopores using three modeling levels that include (1) an all-atom explicit-water model studied with molecular dynamics, and reduced models with implicit water containing (2) hard-sphere ions studied with the Local Equilibrium Monte Carlo simulation method (computing ionic correlations accurately), and (3) point ions studied with Poisson-Nernst-Planck theory (mean-field approximation). We show that reduced models are able to reproduce device functions (rectification and selectivity) for a wide variety of charge patterns, that is, reduced models are useful in understanding the mesoscale physics of the device (i.e., how the current is produced). We also analyze the relationship of the reduced implicit-water models with the explicit-water model and show that diffusion coefficients in the reduced models can be used as adjustable parameters with which the results of the explicit- and implicit-water models can be related. We find that the values of the diffusion coefficients are sensitive to the net charge of the pore but are relatively transferable to different voltages and charge patterns with the same total charge.

The purpose of this study is to provide data for the primitive model of the planar electrical double layer, where ions are modeled as charged hard spheres, the solvent as an implicit dielectric background (with dielectric constant ϵ = 78.5), and the electrode as a smooth, uniformly charged, hard wall. We use canonical and grand canonical Monte Carlo simulations to compute the concentration profiles, from which the electric field and electrostatic potential profiles are obtained by solving Poisson’s equation. We report data for an extended range of parameters including 1:1, 2:1, and 3:1 electrolytes at concentrations c = 0.0001 − 1 M near electrodes carrying surface charges up to σ = ±0.5 Cm⁻². The anions are monovalent with a fixed diameter d₋ = 3 Å, while the charge and diameter of cations are varied in the range z₊ = 1, 2, 3 and d₊ = 1.5, 3, 6, and 9 Å (the temperature is 298.15 K). We provide all the raw data in the supplementary material.

We present a novel molecular dynamics-based simulation technique for investigating gas transport through membranes. In our simulations, the main control parameters are the partial pressure for the components on the input side of the membrane and the total pressure on the output side. The essential point of our scheme is that this pressure control should be realised by adjusting the particle numbers in the input and output side control cells indirectly. Although this perturbation is applied sufficiently far from the membrane, the bulk-phase properties are well controlled in a simulation cell of common size. Numerical results are given for silicalite-1 membrane with permeating CH₄, CO₂, H₂ and N₂ gases as well as with binary mixtures of CO₂ with the other three components. To describe interactions between particles, we used the simple shifted and cut Lennard–Jones potential with parameters available in the literature. It is expected that the proposed technique can be applied to several other types of membranes and transported fluids in order to support the development of a deeper understanding of separation processes.

A systematic simulation study was performed to investigate the interlayer interactions in a 1:1 layered phyllosilicate clay, kaolinite. Atomistic simulations with classical realistic force fields (INTERFACE and ClayFF) were used to examine the influence of the individual non-bonded interactions on the interlayer binding in the kaolinite model system. By switching off selected pairwise interactions in the applied force fields (leaving the intralayer interactions intact), it was confirmed that the tetrahedral–octahedral-type pairwise interactions held the kaolinite plates together and that interlayer hydrogen bonding, modeled by Coulombic forces, played a dominant role in this. Furthermore, it was observed that the number of hydrogen bonds formed had a significant influence on the basal spacing, and thus there was a striking change in the layer–layer interaction strength when there were only two kaolinite plates in the system, rather than several plates, as in real kaolinite particles. Contrary to expectations, the dispersion forces of the studied force fields alone were found to be strong enough to hold the kaolinite plates together.