Molecular Dynamics Study of the Methanol Effect on the Membrane Morphology of Perfluorosulfonic Ionomers

Electro-osmotic Drag and Thermodynamic Properties of Water in Hydrated Nafion Membranes from Molecular Dynamics

One of the important parameters in water management of proton exchange membranes is the electro-osmotic drag (EOD) coefficient of water. The value of the EOD coefficient is difficult to justify, and available literature data on this for Nafion membranes show scattering from in experiments and simulations. Here, we use a classical all-atom model to compute the EOD coefficient and thermodynamic properties of water from molecular dynamics simulations for temperatures between 330 and 420 K, and for different water contents between λ = 5 and λ = 20. λ is the ratio between the moles of water molecules to the moles of sulfonic acid sites. This classical model does not capture the Grotthuss mechanism; however, it is shown that it can predict the EOD coefficient within the range of experimental values for λ = 5 where the vehicular mechanism dominates proton transfer. For λ > 5, the Grotthuss mechanism becomes dominant. To obtain the EOD coefficient, average velocities of water and ions are computed by imposing different electric fields to the system. Our results show that the velocities of water and hydronium scale linearly with the electric field, resulting in a constant ratio of ca. 0.4 within the error bars. We find that the EOD coefficient of water linearly increases from 2 at λ = 5 to 8 at λ = 20 and the results are not sensitive to temperature. The EOD coefficient at λ = 5 is within the range of experimental values, confirming that the model can capture the vehicular transport of protons well. At λ = 20, due to the absence of proton hopping in the model, the EOD coefficient is overestimated by a factor of 3 compared to experimental values. To analyze the interactions between water and Nafion, the partial molar enthalpies and partial molar volumes of water are computed from molecular dynamics simulations. At different water uptakes, multiple linear regression is used on raw simulation data within a narrow composition range of water inside the Nafion membrane. The partial molar volumes and partial molar excess enthalpies of water asymptotically approach the molar volumes and molar excess enthalpies of pure water for water uptakes above 5. This confirms the model can capture the bulklike behavior of water in the Nafion at high water uptakes.

Self-assembly of the cationic surfactant n-hexadecyl-trimethylammonium chloride in methyltrimethoxysilane aqueous solution: classical and reactive molecular dynamics simulations

A flexible aerogel polymerized from methyltrimethoxysilane (MTMS) shows great promise as a high-performance insulator owing to its substantially low thermal conductivity and mechanical flexibility, attributed to its porous microstructure and organic-inorganic hybridization, respectively, which promote its industrial applications. Conventionally, the cationic surfactant n-hexadecyltrimethylammonium chloride (CTAC) is utilized to experimentally control the nanoscale microstructure and, consequently, the flexibility of the MTMS aerogel; however, the mechanism through which CTAC prevents MTMS aggregation in the solution is not yet fully understood. This study unravels the role of CTAC in preventing MTMS aggregation in aqueous solution using both classical and reactive molecular dynamics simulations. We found that CTAC molecules can form self-aggregates even when the polymerization of MTMS progresses and then the MTMS-derived oligomer turns to be hydrophobic in aqueous solution. In summary, the self-assemblies of CTAC disperse among the MTMS associations and effectively prevent MTMS clustering, and this is considered as the key mechanism underlying the formation of a flexible microstructure of the hybrid aerogel.

Molecular dynamics simulations of sodium alginate/sulfonated graphene oxide membranes properties

The influence of methanol as a solvent on the properties of sodium alginate/sulfonated graphene oxide (SA/SGO) membranes was explored in water-methanol mixed conditions with various methanol concentrations and temperatures through molecular dynamics simulations. The methanol uptake of the membrane showed an isolation phase determined from the simulation results. The distance between the sulfonic acid groups increased in higher methanol concentrations, as observed from S-S RDFs. Furthermore, the distance between the SA-chain RDFs and the solvent molecules was analysed to determine a) the affinity of water towards the sulfonic acid groups and b) the affinity of the aromatic backbone of the SA towards methanol molecules. A decrease in water molecule diffusion led to an increase in methanol diffusion and uptake. SA/SGO membranes exhibited a smaller diffusion coefficient than that for the Nafion membranes, as calculated from simulation results and compared to the experimental work. Additionally, the diffusion ability increased at higher temperatures for all permeants. The interaction information obtained is useful for DMFC applications.

Hydration and proton transfer in highly sulfonated poly(phenylene sulfone) ionomers: an ab initio study

The need to operate proton exchange membrane fuel cells under hot and dry conditions has driven the synthesis and testing of sulfonated poly(phenylene) sulfone (sPSO(2)) ionomers. The primary hydration and energetics associated with the transfer of protons in oligomeric fragments of two sPSO(2)ionomers were evaluated through first-principles electronic structures calculations. Our results indicate that the interaction between neighboring sulfonic acid groups affect both theconformation and stability of the fragments. The number of water molecules required to affect the transfer of a proton in the first hydration shell was observed to be a function of the hydrogen bonding in proximity of the sulfonic acid groups: three H(2)O for the meta- and four H(2)O for the ortho-conformations. Calculations of the rotational energy surfaces indicate that the aromatic backbones of sPSO(2) are much stiffer than the polytetrafluoroethylene (PTFE) backbones in perfluorosulfonic acid (PFSA) ionomers: the largest energy penalty for rotating phenylene rings (i.e., 15.5 kcal/mol for ortho-ortho-sPSO(2)) is nearly twice that computed for the rotation of a CF(2) unit in a PTFE backbone. The energetics for the transfer of various protons in proximity to one or two sulfonate groups (-SO(3)(-)) was also determined. The computed energy barrier for proton transfer when only one sulfonic acid group is present is approximately 1.9 kcal/mol, which is 2.1 kcal/mol lower than similar calculations for PFSA systems. When two sulfonic acid groups are bridged by water molecules, a symmetric bidirectional transfer occurs, which gives a substantially small energy barrier of only 0.7 kcal/mol.

Comparison of the hydration and diffusion of protons in perfluorosulfonic acid membranes with molecular dynamics simulations. scui@utk.edu

Classical molecular dynamics (MD) simulations were performed to determine the hydrated morphology and hydronium ion diffusion coefficients in two different perfluorosulfonic acid (PFSA) membranes as functions of water content. The structural and transport properties of 1143 equivalent weight (EW) Nafion, with its relatively long perfluoroether side chains, are compared to the short-side-chain (SSC) PFSA ionomer at an EW of 977. The separation of the side chains was kept uniform in both ionomers consisting of -(CF 2) 15- units in the backbone, and the degree of hydration was varied from 5 to 20 weight % water. The MD simulations indicated that the distribution of water clusters is more dispersed in the SSC ionomer, which leads to a more connected water-channel network at the low water contents. This suggests that the SSC ionomer may be more inclined to form sample-spanning aqueous domains through which transport of water and protons may occur. The diffusion coefficients for both hydronium ions and water molecules were calculated at hydration levels of 4.4, 6.4, 9.6, and 12.8 H 2O/SO 3H for each ionomer. When compared to experimental proton diffusion coefficients, this suggests that as the water content is increased the contribution of proton hopping to the overall proton diffusion increases.

Lattice Boltzmann simulations for proton transport in 2-D model channels of Nafion

Proton conductance in a 2-D channel with a slab-like structure was studied to verify that the lattice Boltzmann method (LBM) can be used as a simulation tool for proton conduction in a Nafion membrane, which is a mesoscopic system with a highly disordered porous structure. Diffusion resulting from a concentration gradient and migration by an electrostatic force were considered as the origins of proton transport. The electrostatic force acting on a proton was computed by solving the Poisson equation. The proton concentration in the membrane is expressed as a continuous function and the sulfonic charge is placed discretely. The space-averaged conductance of protons in a nonequilibrium stationary state was evaluated as a function of the structural parameters: namely, channel width and distribution of the sulfonic groups. The resulting space-averaged conductance deviates from the bulk values, depending particularly on the sulfonic group distribution. Details of the simulation scheme are described and the applicability of the present scheme to real membranes is discussed.

Molecular dynamics studies of the Nafion®, Dow® and Aciplex® fuel-cell polymer membrane systems

The Nafion, Dow and Aciplex systems--where the prime differences lies in the side-chain length--have been studied by molecular dynamics (MD) simulation under standard pressure and temperature conditions for two different levels of hydration: 5 and 15 water molecules per (H)SO(3) end-group. Structural features such as water clustering, water-channel dimensions and topology, and the dynamics of the hydronium ions and water molecules have all been analysed in relation to the dynamical properties of the polymer backbone and side-chains. It is generally found that mobility is promoted by a high water content, with the side-chains participating actively in the H(3)O(+)/H(2)O transport mechanism. Nafion, whose side-chain length is intermediate of the three polymers studied, is found to have the most mobile polymer side-chains at the higher level of hydration, suggesting that there could be an optimal side-chain length in these systems. There are also some indications that the water-channel network connectivity is optimal for high water-content Nafion system, and that this could explain why Nafion appears to exhibit the most favourable overall hydronium/water mobility.

Modelling of morphology and proton transport in PFSA membranes

Computational modelling studies of the structure of perfluorosulfonic acid (PFSA) ionomer membranes consistently exhibit a nanoscopic phase-separated morphology in which the ionic side chains and aqueous counterions segregate from the fluorocarbon backbone to form clusters or channels. Although these investigations do not unambiguously predict the size or shape of the clusters, and whether or not the channels percolate the matrix or if the connections between them are more transient, the sequence of co-monomers along the main chain appears strongly to influence the domain size of the ionic regions, with more blocky sequences giving rise to larger domain sizes. The fundamental insight that substantial rearrangement of the sulfonic acid terminated side chains and fluorocarbon backbone takes place during swelling or shrinkage is borne out by both molecular and mesoscale simulations of model PFSA polymers, along with ab initio electronic structure calculations of minimally hydrated oligomeric fragments. Molecular-level modelling of proton transport in PFSA membranes attests to the complexity of the underlying mechanisms and the need to examine the chemical and physical processes at several distinct time and length scales. These investigations have revealed that the conformation of the fluorocarbon backbone, flexibility of the sidechains, and degree of aggregation and association of the sulfonic acid groups under minimally hydrated conditions collectively control the dissociation of the protons and the formation of Zundel and Eigen cations. The former appear to be the dominant charge carriers when the limiting water content allows only for the formation of a contact ion pair with the tethered sulfonate anion. As the water content increases, solvent-separated Eigen ions begin to appear, indicating that the dominant mechanism for diffusion of protons occurs over a region approximately 4 A away from the sulfonate groups. Finally, both the vehicular and Grotthuss shuttling mechanisms contribute to the mobility of the protons but, surprisingly, they are not always correlated, resulting in a lower overall diffusion coefficient. In summary, as the preceding observations indicate, the state of computational modelling of PFSA membranes has progressed sufficiently over the last decade to enable its use as a powerful predictive tool with which to guide the process of designing novel membrane materials for fuel cell applications.

Alcohol and Proton Transport in Perfluorinated Ionomer Membranes for Fuel Cells

To clarify the transport mechanisms of alcohols and proton in perfluorosulfonated ionomer (PFSI) membranes for fuel cells, four membranes having different equivalent weight (EW) values were examined. Membranes were immersed in methanol, ethanol, and 2-propanol to prepare a total of 12 samples, and membrane swelling, mass (alcohol and proton) transports, and interactions between alcohols and proton were investigated systematically in the fully penetrated state. The membrane expansion fraction theta and alcohol content lambda increased with decreasing the EW value for all the samples. The self-diffusion coefficients (D's) of the alkyl group and of OH (including protons) were measured separately by the pulsed-gradient spin-echo (PGSE)-NMR method and the D's also increased with decreasing the EW value. These results implied that the alcohols penetrate into the hydrophilic regions of the PFSI membranes and diffuse through the space expanded by the alcohols. The ionic cluster regions formed by the alcohols resemble those induced by water in the water swollen membrane, where protons dissociated from sulfonic acid groups transport through the regions together with water molecules. The D values decreased with increasing the molecular weight of alcohols. This trend was supported by activation energies Ea estimated from the Arrhenius plots of D in the temperature range from 30 to -40 degrees C. The PGSE-NMR measurements also revealed that protons move faster than the alkyl groups in the membranes. The proton transport by the Grotthuss (hopping) mechanism was facilitated by the increase of the alcohol content and the decrease of the molecular weight. This result was also supported by the experimental results of proton conductivity kappa and mobility u(H(+)). Density functional theory (DFT) calculations of the interaction energy DeltaE(int) between proton and alcohol (including OH) showed that the /DeltaE(int)/ increases with increasing the molecular weight of alcohols, which is in a inverse relationship with the kappa and u(H(+)) values. The proton transport depends strongly on the DeltaE(int) in the membranes.