Source author record

Dirk Gillespie

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

Biological Physics Biomolecules cond-mat.soft cond-mat.stat-mech physics.chem-ph physics.comp-ph

Catalog footprint

Research overview

3works

6topics

4close collaborators

Actions

Is this your work?

Open graph Browse papers

preprint2022arXiv

Individual ion species chemical potentials in the Mean Spherical Approximation

The Mean Spherical Approximation (MSA) is a commonly-used thermodynamic theory for computing the energetics of ions in the primitive model (i.e., charged hard-sphere ions in a background dielectric). For the excess chemical potential, however, the early MSA formulations (which were widely adopted) only included the terms needed to compute the mean excess chemical potential (or the mean activity coefficient). Other terms for the chemical potential $μ_i$ of individual species $i$ were not included because they sum to $0$ in the mean chemical potential. Here, we derive these terms to give a complete MSA formulation of the chemical potential. The result is a simple additive term for $μ_i$ that we show is a qualitative improvement over the previous MSA version. In addition, our derivation shows that the MSA's assumption of global charge neutrality is not strictly necessary, so that the MSA is also valid for systems close to neutrality.

preprint2013arXiv

Dynamic Monte Carlo simulation of coupled transport through a narrow multiply-occupied pore

Dynamic Monte Carlo simulations are used to study coupled transport (co-transport) through sub-nanometer-diameter pores. In this classic Hodgkin-Keynes mechanism, an ion species uses the large flux of an abundant ion species to move against its concentration gradient. The efficiency of co-transport is examined for various pore parameters so that synthetic nanopores can be engineered to maximize this effect. In general, the pore must be narrow enough that ions cannot pass each other and the charge of the pore large enough to attract many ions so that they exchange momentum. Co-transport efficiency increases as pore length increases, but even very short pores exhibit co-transport, in contradiction to the usual perception that long pores are necessary. The parameter ranges where co-transport occurs is consistent with current and near-future synthetic nanopore geometry parameters, suggesting that co-transport of ions may be a new application of nanopores.

preprint2009arXiv

An Efficient Algorithm for Classical Density Functional Theory in Three Dimensions: Ionic Solutions

Classical density functional theory (DFT) of fluids is a valuable tool to analyze inhomogeneous fluids. However, few numerical solution algorithms for three-dimensional systems exist. Here we present an efficient numerical scheme for fluids of charged, hard spheres that uses $\mathcal{O}(N\log N)$ operations and $\mathcal{O}(N)$ memory, where $N$ is the number of grid points. This system-size scaling is significant because of the very large $N$ required for three-dimensional systems. The algorithm uses fast Fourier transforms (FFT) to evaluate the convolutions of the DFT Euler-Lagrange equations and Picard (iterative substitution) iteration with line search to solve the equations. The pros and cons of this FFT/Picard technique are compared to those of alternative solution methods that use real-space integration of the convolutions instead of FFTs and Newton iteration instead of Picard. For the hard-sphere DFT we use Fundamental Measure Theory. For the electrostatic DFT we present two algorithms. One is for the \textquotedblleft bulk-fluid\textquotedblright functional of Rosenfeld [Y. Rosenfeld. \textit{J. Chem. Phys.} 98, 8126 (1993)] that uses $\mathcal{O}(N\log N)$ operations. The other is for the \textquotedblleft reference fluid density\textquotedblright (RFD) functional [D. Gillespie et al., J. Phys.: Condens. Matter 14, 12129 (2002)]. This functional is significantly more accurate than the bulk-fluid functional, but the RFD algorithm requires $\mathcal{O}(N^{2})$ operations.

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.