Proton transfer and the diffusion of H+ and OH- ions along water wires.

Hydrogen and hydroxide ion transport in narrow carbon nanotubes (CNTs) of diameter 8.1 Å and lengths up to 582 Å are investigated by molecular dynamics simulations using a dissociating water model. The diffusion coefficients of the free ions in an open chain are significantly larger than in periodically replicated wires that necessarily contain D or L end defects, and both are higher than they are in bulk water. The free hydroxide ion diffuses faster than the free hydronium ion in short CNTs, unlike diffusion in liquid water, and both coefficients increase and converge to nearly the same value with increasing tube length. The diffusion coefficients of the two ions increase further when the tubes are immersed in a water reservoir and they move easily out of the tube, suggesting an additional pathway for proton transport via OH(-) ions in biological channels.

Proton transfer and the mobilities of the H+ and OH- ions from studies of a dissociating model for water.

Hydrogen (H(+)) and hydroxide (OH(-)) ions in aqueous solution have anomalously large diffusion coefficients, and the mobility of the H(+) ion is nearly twice that of the OH(-) ion. We describe molecular dynamics simulations of a dissociating model for liquid water based on scaling the interatomic potential for water developed by Ojamäe-Shavitt-Singer from ab initio studies at the MP2 level. We use the scaled model to study proton transfer that occurs in the transport of hydrogen and hydroxide ions in acidic and basic solutions containing 215 water molecules. The model supports the Eigen-Zundel-Eigen mechanism of proton transfer in acidic solutions and the transient hyper-coordination of the hydroxide ion in weakly basic solutions at room temperature. The free energy barriers for proton transport are low indicating significant proton delocalization accompanying proton transfer in acidic and basic solutions. The reorientation dynamics of the hydroxide ion suggests changes in the proportions of hyper-coordinated species with temperature. The mobilities of the hydrogen and hydroxide ions and their temperature dependence between 0 and 50 °C are in excellent agreement with experiment and the reasons for the large difference in the mobilities of the two ions are discussed. The model and methods described provide a novel approach to studies of liquid water, proton transfer, and acid-base reactions in aqueous solutions, channels, and interfaces.

Mesoscopic description of solvent effects on polymer dynamics.

Solvent effects on polymer dynamics and structure are investigated using a mesoscopic solvent model that accounts for hydrodynamic interactions among the polymer beads. The simulation method combines molecular dynamics of the polymer chain, interacting with the solvent molecules through intermolecular forces, with mesoscopic multiparticle collision dynamics for the solvent molecules. Changes in the intermolecular forces between the polymer beads and mesoscopic solvent molecules are used to vary the solvent conditions from those for good to poor solvents. Polymer collapse and expansion dynamics following changes in solvent conditions are studied for homopolymer and block copolymer solutions. The frictional properties of polymers are also investigated.

Two-particle friction in a mesoscopic solvent.

The effects of hydrodynamic interactions on the friction tensors for two particles in solution are studied. The particles have linear dimensions on nanometer scales and are either simple spherical particles interacting with the solvent through repulsive Lennard-Jones forces or are composite cluster particles whose atomic components interact with the solvent through repulsive Lennard-Jones forces. The solvent dynamics is modeled at a mesoscopic level through multiparticle collisions that conserve mass, momentum, and energy. The dependence of the two-particle relative friction tensors on the interparticle separation indicates the importance of hydrodynamic interactions for these nanoparticles.

Structure of nonuniform polymer melts: density functional perturbation approach.

A density functional perturbation approximation for polyatomic molecules, which is based on the fundamental-measure theory for hard-core repulsion and the hybrid weighted-density approximation for chain connectivity, was proposed to clarify the structure of polymer melts at interfaces. It was applied to predict the local density distributions, adsorption isotherms, and surface excess of a freely jointed tangent hard-sphere chain in hard slit pores. Wertheim's first-order perturbation theory extended by J. Chem. Phys. 112, 2368 (2002)] was used to calculate the weight function and second-order direct correlation function due to the chain connectivity. The theoretical results are in excellent agreement with the computer simulations.

Friction and diffusion of a Brownian particle in a mesoscopic solvent.

The friction and diffusion coefficients of a massive Brownian particle in a mesoscopic solvent are computed from the force and the velocity autocorrelation functions. The mesoscopic solvent is described in terms of free streaming of the solvent molecules, interrupted at discrete time intervals by multiparticle collisions that conserve mass, momentum, and energy. The Brownian particle interacts with the solvent molecules through repulsive Lennard-Jones forces. The decays of the force and velocity autocorrelation functions are analyzed in the microcanonical ensemble as a function of the number N of solvent molecules and Brownian particle mass and diameter. The simulations are carried out for large system sizes and long times to assess the N dependence of the friction coefficient. The decay rates of these correlations are confirmed to vary as N(-1) in accord with earlier predictions. Hydrodynamic effects on the velocity autocorrelation function and diffusion coefficient are studied as a function of Brownian particle mass and diameter.