The pair-correlation functions for fluid ionic mixtures in arbitrary spatial dimensions are computed in hypernetted chain (HNC) approximation. In the primitive model, all ions are approximated as non-overlapping hyperspheres with Coulomb interactions
. Our spectral HNC solver is based on a Fourier-Bessel transform introduced by Talman [J. Comput. Phys., 29, 35 (1978)], with logarithmically spaced computational grids. Numeric efficiency for arbitrary spatial dimensions is a commonly exploited virtue of this transform method. Here, we highlight another advantage of logarithmic grids, consisting in efficient sampling of pair-correlation functions for highly asymmetric ionic mixtures. For three-dimensional fluids, ion size- and charge-ratios larger than one thousand can be treated, corresponding to hitherto computationally not accessed micrometer-sized colloidal spheres in 1-1 electrolyte. Effective colloidal charge numbers are extracted from our primitive model results. For moderately large ion size- and charge-asymmetries, we present Molecular Dynamics simulation results that agree well with the approximate HNC pair correlations.
For sedimenting colloidal hard spheres, the propagation and broadening of the crystal-fluid interface is studied by Brownian dynamics computer simulations of an initially homogeneous sample. Two different types of interface broadenings are observed:
the first occurs during growth and is correlated with the interface velocity, the second is concomitant with the splitting of the crystal-fluid interface into the crystal-amorphous and amorphous-liquid interfaces. The latter width is strongly peaked as a function of the gravitational driving strength with a huge amplitude relative to its equilibrium counterpart.
We have performed simulations to study how increasing humidity affects the structure of Nafion-like ionomers under conditions of low sulfonate concentration and low humidity. At the onset of membrane hydration, the clusters split into smaller parts.
These subsequently swell, but then maintain constant the number of sulfonates per cluster. We find that the distribution of water in low-sulfonate membranes depends strongly on the sulfonate concentration. For a relatively low sulfonate concentration, nearly all the side-chain terminal groups are within cluster formations, and the average water loading per cluster matches the water content of membrane. However, for a relatively higher sulfonate concentration the water-to-sulfonate ratio becomes non-uniform. The clusters become wetter, while the inter-cluster bridges become drier. We note the formation of unusual shells of water-rich material that surround the sulfonate clusters.
Coarse-grained molecular-dynamics simulations were used to study the morphological changes induced in a Nafion$^{tiny textregistered}$-like ionomer by the imposition of a strong electric field. We observe the formation of novel structures aligned a
long the direction of the applied field. The polar head groups of the ionomer side chains aggregate into clusters, which then form rod-like formations which assemble into a hexatic array aligned with the direction of the field. Occasionally these lines of sulfonates and protons form a helical structure. Upon removal of the electric field, the hexatic array of rod-like structures persists, and has a lower calculated free energy than the original isotropic morphology.
This simulation study investigates the dependence of the structure of dry Nafion$^{tinytextregistered}$-like ionomers on the electrostatic interactions between the components of the molecules. In order to speed equilibration, a procedure was adopted
which involved detaching the side chains from the backbone and cutting the backbone into segments, and then reassembling the macromolecule by means of a strong imposed attractive force between the cut ends of the backbone, and between the non-ionic ends of the side chains and the midpoints of the backbone segments. Parameters varied in this study include the dielectric constant, the free volume, side-chain length, and strength of head-group interactions. A series of coarse-grained mesoscale simulations shows the morphlogy to depend sensitively on the ratio of the strength of the dipole-dipole interactions between the side-chain acidic end groups to the strength of the other electrostatic components of the Hamiltonian. Examples of the two differing morphologies proposed by Gierke and by Gebel emerge from our simulations.
