Simulation study of field-induced morphological changes in a proton-conducting ionomer

Atomistic simulations of the structure of highly crosslinked sulfonated poly(styrene-co-divinylbenzene) ion exchange resins

The microscopic structures of highly crosslinked sulfonated poly(styrene-co-divinylbenzene) resins have been modeled by generating atomistic microstructures using stochastic-like algorithms, which are subsequently relaxed using molecular dynamics. Two different generation algorithms have been tested. The relaxation of the microstructures generated by the first algorithm, which is based on a homogeneous construction of the resin, leads to a significant overestimation of the experimental density as well as to an unsatisfactory description of the porosity. In contrast, the generation approach that combines algorithms for the heterogeneous growing and branching of the chains enables the formation of crosslinks with different topologies. In particular, the intrinsic heterogeneity observed in these resins is efficiently reproduced when the topological loops, which are defined by two or more crosslinks closing a cycle, are present in their microscopic description. Thus, the apparent density, porosity and pore volume estimated using microstructures with these topological loops, called super-crosslinks, are in very good agreement with the experimental results. Although the backbone dihedral angle distribution of the generated and relaxed models is not influenced by the topology, the number and type of crosslinks affect the medium- and long-range atomic disposition of the backbone atoms and the distribution of sulfonic groups. An analysis of the distribution of the local density indicates that super-crosslinks are responsible for the heterogeneous homogenization observed during the MD relaxation. Finally the π-π stacking interactions have been analyzed. Results indicate that those in which the two rings adopt a T-shaped disposition are considerably more abundant as compared to those with the co-facially oriented rings, independently of the resin topology.

Self-assembly in Nafion membranes upon hydration: water mobility and adsorption isotherms

By means of dissipative particle dynamics (DPD) and Monte Carlo (MC) simulations, we explored geometrical, transport, and sorption properties of hydrated Nafion-type polyelectrolyte membranes. Composed of a perfluorinated backbone with sulfonate side chains, Nafion self-assembles upon hydration and segregates into interpenetrating hydrophilic and hydrophobic subphases. This segregated morphology determines the transport properties of Nafion membranes that are widely used as compartment separators in fuel cells and other electrochemical devices, as well as permselective diffusion barriers in protective fabrics. We introduced a coarse-grained model of Nafion, which accounts explicitly for polymer rigidity and electrostatic interactions between anionic side chains and hydrated metal cations. In a series of DPD simulations with increasing content of water, a classical percolation transition from a system of isolated water clusters to a 3D network of hydrophilic channels was observed. The hydrophilic subphase connectivity and water diffusion were studied by constructing digitized replicas of self-assembled morphologies and performing random walk simulations. A non-monotonic dependence of the tracer diffusivity on the water content was found. This unexpected behavior was explained by the formation of large and mostly isolated water domains detected at high water content and high equivalent polymer weight. Using MC simulations, we calculated the chemical potential of water in the hydrated polymer and constructed the water sorption isotherms, which extended to the oversaturated conditions. We determined that the maximum diffusivity and the onset of formation of large water domains corresponded to the saturation conditions at 100% humidity. The oversaturated membrane morphologies generated in the canonical ensemble DPD simulations correspond to the metastable and unstable states of Nafion membrane that are not realized in the experiments.

Ion effects at electrode/solid polymer electrolyte membrane interfaces

The differential capacity and the potential distribution at electrode/solid polymer electrolyte membrane/solution interfaces are calculated through an analytical approach. The model considers coions' and counterions' permeation through the membrane from the solvent phase and the ions' partitioning equilibrium at the SPEM/solution interface. The latter effects are included by incorporating the Donnan equilibrium, the steric hindrance, the solvation energy change when ions move from water to membrane pores and ion electrostatic interactions. It is shown that capacitance maxima in capacitance-potential curves may appear because of the acid-base dissociation process inside the membrane and the change in the ions' total interaction energy with the applied potential. For low dielectric constants inside membrane pores, εp, sharp peaks can be obtained. These peaks broaden, decrease in magnitude and shift to positive potentials once εp is increased. Finally, model predictions are discussed in light of recent experimental data obtained on Nafion® covered Pt(111) electrodes, providing a theoretical framework for the qualitative electroanalysis of these systems.

Transport in supported polyelectrolyte brushes

Functional polymers have a wide variety of applications ranging from energy storage to drug delivery. For energy storage applications, desirable material properties include low cost, high charge storage and/or mobility, and low rates of degradation. Isotropic thin films have been used for many of these types of applications, but research suggests that different structures such as polymer brushes can improve charge transport by an order of magnitude. Supported polymer brush structures produced by "grafting-from" polymerization methods offer a framework for a controlled study of these materials on the molecular scale. Using these materials, researchers can study the basis of hindered diffusion because they contain a relatively homogeneous polyelectrolyte membrane. In addition, researchers can use fluorescent molecular probes with different charges to examine steric and Coulombic contributions to transport near and within polymer brushes. In this Account, we discuss recent progress in using fluorescence correlation spectroscopy, single-molecule polarization-resolved spectroscopy, and a novel three-dimensional orientational technique to understand the transport of charged dye probes interacting with the strong polyanionic brush, poly(styrene sulfonate). Our preliminary experiments demonstrate that a cationic dye, Rhodamine 6G, probes the brush as a counterion, and diffusion is therefore dominated by Coulombic forces, which results in a 10,000-fold decrease in the diffusion coefficient in comparison with free diffusion. We also support our experimental results with molecular dynamics simulations. Further experiments show that, up to 50% of the time, Rhodamine 6G translates within the brush without significant rotational diffusion, which indicates a strong deviation from Fickian transport mechanisms (in which translational and rotational diffusion are related directly through parameters such as chemical potential, size, solution viscosity, and thermal properties). To understand this oriented transport, we discuss the development of an experimental technique that allows us to quantify the three-dimensional orientation on the time scale of intrabrush transport. This method allowed us to identify a unique orientational transport direction for Rhodamine 6G within the poly(styrene sulfonate) brush and to report preliminary evidence for orientational dye "hopping".

Simulation study of poled low-water ionomers with different architectures

The role of the ionomer architecture in the formation of ordered structures in poled membranes is investigated by molecular dynamics computer simulations. It is shown that the length of the sidechain L(s) controls both the areal density of cylindrical aggregates N(c) and the diameter of these cylinders in the poled membrane. The backbone segment length L(b) tunes the average diameter D(s) of cylindrical clusters and the average number of sulfonates N(s) in each cluster. A simple empirical formula is noted for the dependence of the number density of induced rod-like aggregates on the sidechain length L(s) within the parameter range considered in this study.

A simulation study of field-induced proton-conduction pathways in dry ionomers

The morphological changes that can be induced in a dry ionomer by application of a strong electric field have been studied by means of computer simulation. The internal energy of the membrane at first slowly decreases with increasing field, but then rapidly increases after a certain threshold field is reached. This effect is interpreted as the reorganization of interacting head group dipoles in response to the external perturbation. The resulting morphology contains continuous channels of hydrophilic material capable of facilitating proton conduction. Upon removal of the poling field, the system does not return to its original morphology, but retains the anisotropic structure of the poled material. The poled structure appears to be thermodynamically stable, as confirmed by calculations of the Helmholtz energy of the original and poled samples.