Molecular Modeling of Proton Transport in the Short-Side-Chain Perfluorosulfonic Acid Ionomer

Morphological Effect of Side Chain Length in Sulfonated Poly(arylene ether sulfone)s Polymer Electrolyte Membranes via Molecular Dynamics Simulation

With the recognition of the multiple advantages of sulfonated hydrocarbon-based polymers that possess high chemical and mechanical stability with significant low cost, we employed molecular dynamics simulation to explore the morphological effects of side chain length in sulfonated polystyrene grafted poly(arylene ether sulfone)s (SPAES) proton exchange membranes. The calculated diffusion coefficients of hydronium ions (H₃O⁺) are in range of 0.61-1.15 × 10^-7 cm²/s, smaller than that of water molecules, due to the electrical attraction between the oppositely charged sulfonate group and H₃O⁺. The investigation into the radial distribution functions suggests that phase segregation in the SPAES membrane is more probable with longer side chains. As the hydration level of the membranes in this study is relatively low (λ = 3), longer side chains correspond to more water molecules in the amorphous cell, which provides better solvent effects for the distribution of sulfonated side chains. The coordination number of water molecules and hydronium ions around the sulfonate group increases from 1.67 to 2.40 and from 2.45 to 5.66, respectively, with the increase in the side chain length. A significant proportion of the hydronium ions appear to be in bridging configurations coordinated by multiple sulfonate groups. The microscopic conformation of the SPAES membrane is basically unaffected by temperature during the evaluated temperature range. Thus, it can be revealed that the side chain length plays a key role in the configuration of the polymer chain and would contribute to the formation of the microphase separation morphology, which profits proton transport in the hydrophilic domains.

Evaluation of proton transport and solvation effect in hydrated Nafion membrane with degradation

In proton exchange membrane fuel cells (PEMFCs), free radicals easily attack ionomers, resulting in membrane degradation.

Effects of Hydration and Temperature on the Microstructure and Transport Properties of Nafion Polyelectrolyte Membrane: A Molecular Dynamics Simulation

To investigate the effects of temperature and hydration on the microstructure of polymer electrolyte membrane and the transport of water molecules and hydronium ions, molecular dynamics simulations are performed on Nafion 117 for a series of water contents at different temperatures. The interactions among the sulfonate groups, hydronium ions, and water molecules are studied according to the analysis of radial distribution functions and coordination numbers. The sizes and connectivity of water clusters are also discussed, and it is found that the hydration level plays a key role in the phase separation of the membrane. However, the effect of the temperature is slight. When the water content increases from 3.5 to 16, the size of water clusters in the membrane increases, and the clusters connect to each other to form continuous channels for diffusion of water molecules and hydronium ions. The diffusion coefficients are estimated by studying the mean square displacements. The results show that the diffusion of water molecules and hydronium ions are both enhanced by the increase of the temperature and hydration level. Furthermore, the diffusion coefficient of water molecules is always much larger than that of hydronium ions. However, the ratio of the diffusion coefficient of water molecules to that of hydronium ions decreases with the increase of water content.

Morphological effect of side chain on H3O+ transfer inside polymer electrolyte membranes across polymeric chain via molecular dynamics simulation

Performance and durability of polymer electrolyte membrane are critical to fuel cell quality. As fuel cell vehicles become increasingly popular, membrane fundamentals must be understood in detail. Here, this study used molecular dynamic simulations to explore the morphological effects of perfluorosulfonic acid (PFSA)-based membranes on ionic conductivity. In particular, I developed an intuitive quantitative approach focusing principally on hydronium adsorbing to, and desorbing from, negatively charged sulfonate groups, while conventional ionic conductivity calculations featured the use of mean square displacements that included natural atomic vibrations. The results revealed that shorter side-chains caused more hydroniums to enter the conductive state, associated with higher ion conductivity. In addition, the hydronium path tracking showed that shorter side-chains allowed hydroniums to move among host groups, facilitating chain adsorption, in agreement with a mechanism suggested in earlier studies.

Simple and Reliable Method for Studying the Adsorption Behavior of Aquivion Ionomers on Carbon Black Surfaces

A better understanding of the interactions of carbon black and perfluorinated sulfonic acid (PFSA) ionomer helps to improve the effectiveness of polymer electrolyte membrane fuel cells. We present a simple and fast method for quantitative PFSA ionomer analysis based on suspension density measurements. After validation of the reliability of our method by thermogravimetric analysis and isothermal titration calorimetry (ITC), we investigate the adsorption equilibrium of short-side-chain PFSA ionomers of different equivalent weights (EW) and polarities on carbon black. The measured adsorption isotherms exhibit a plateau in the ionomer surface concentration for ionomer equilibrium concentrations ≤2 g/L. In this concentration range, the adsorption isotherms are described by the Langmuir model, whereby the surface concentrations in the plateau region are between 0.041 and 0.070 g/g. The plateau value of the ionomer surface concentration increases with EW and therefore with decreasing number of side chains with terminal sulfonic acid group per ionomer molecule, while the amount of adsorbed sulfonic acid groups remains constant for all investigated ionomers, resulting in similar ζ-potentials and sedimentation stability of the suspensions. The free energies of adsorption Δ G calculated from the association constants of the adsorption isotherms agree well with Δ G values obtained by isothermal titration calorimetry (ITC) and thus validate the adsorption isotherm measurement method. From the values of adsorption enthalpy Δ H ((-7.3 ± 0.8) kJ/mol) and entropy Δ S (ca. 100 J/(mol K)), which were extracted from ITC, we conclude that the ionomer adsorption on carbon black is a spontaneous physisorption process.

Ab initio and density functional theory (DFT) studies on triflic acid with water and protonated water clusters

The structure, stability and infrared spectral signatures of triflic acid (TA) with water clusters (W_n) and protonated water clusters (TAH⁺W_n, n = 1 - 6) were computed using DFT and MP2 methods. Our calculations show that a minimum of three water molecules are necessary to stabilize the dissociated zwitterionic form of TA. It can be seen from the results that there is no significant movement of protons in smaller (n = 1 and 2) and linear (n = 1 - 6) types of water clusters. Further, the geometries of TAW_n clusters first form a neutral pair (NP) to contact ion pair (CIP), then form a solvent separated ion pair (SSIP) in a water hexamer. These findings reveal that proton transfer may take place through NP to CIP and then CIP to SSIP. The calculated binding energies (BEs) of ion pair clusters is always higher than that of NP clusters (i.e., more stable than the NP). Existing excess proton linear chain clusters transfer a proton to adjacent water molecules via a Grotthuss mechanism, whereas the same isomers in the branched motifs do not conduct protons. Examination of geometrical parameters and infrared frequencies reveals hydronium ion (H₃O⁺ also called Eigen cation) formation in both TAW_n and protonated TAW_n clusters. The stability of Eigen water clusters is three times higher than that of other non-Eigen water clusters. Our study shows clearly that formation of ion pairs in TAW_n and TAH⁺W_n clusters greatly favors proton transfer to neighboring water molecules and also enhances the stability of these complexes.

Ab initio molecular dynamics simulations of aqueous triflic acid confined in carbon nanotubes

Ab initio molecular dynamics simulations were performed to investigate the effects of nanoscale confinement on the structural and dynamical properties of aqueous triflic acid (CF3SO3H). Single-walled carbon nanotubes (CNTs) with diameters ranging from ∼11 to 14 Å were used as confinement vessels, and the inner surface of the CNT were either left bare or fluorinated to probe the influence of the confined environment on structural and dynamical properties of the water and triflic acidic. The systems were simulated at hydration levels of n = 1-3 H2O/CF3SO3H. Proton dissociation expectedly increased with increasing hydration. Along with the level of hydration, hydrogen bond connectivity between the triflic acid molecules, both directly and via a single water molecule, played a role on proton dissociation. Direct hydrogen bonding between the CF3SO3H molecules, most commonly found in the larger bare CNT, also promoted interactions between water molecules allowing for greater separation of the dissociated protons from the CF3SO3(-) as the hydration level was increased. However, this also resulted in a decrease in the overall proportion of dissociated protons. The confinement dimensions altered both the hydrogen bond network and the distribution of water molecules where the H2O in the fluorinated CNTs tended to form small clusters with less proton dissociation at n = 1 and 2 but the highest at n = 3. In the absence of nearby hydrogen bond accepting sites from H2O or triflic acid SO3H groups, the water molecules formed weak hydrogen bonds with the fluorine atoms. In the bare CNT systems, these involved the CF3 groups of triflic acid and were more frequently observed when direct hydrogen bonding between CF3SO3H hindered potential hydrogen bonding sites. In the fluorinated tubes, interactions with the covalently bound fluorine atoms of the CNT wall dominated which appear to stabilize the hydrogen bond network. Increasing the hydration level increased the frequency of the OH···F (CNT) hydrogen bonding which was highly pronounced in the smaller fluorinated CNT indicating an influence on the confinement dimensions on these interactions.

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.

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.

Proton transport in triflic acid pentahydrate studied via ab initio path integral molecular dynamics

Trifluoromethanesulfonic acid hydrates provide a well-defined system to study proton dissociation and transport in perfluorosulfonic acid membranes, typically used as the electrolyte in hydrogen fuel cells, in the limit of minimal water. The triflic acid pentahydrate crystal (CF(3)SO(3)H·5H(2)O) is sufficiently aqueous that it contains an extended three-dimensional water network. Despite it being extended, however, long-range proton transport along the network is structurally unfavorable and would require considerable rearrangement. Nevertheless, the triflic acid pentahydrate crystal system can provide a clear picture of the preferred locations of local protonic defects in the water network, which provides insights about related structures in the disordered, low-hydration environment of perfluorosulfonic acid membranes. Ab initio molecular dynamics simulations reveal that the proton defect is most likely to transfer to the closest water that has the expected presolvation and only contains water in its first solvation shell. Unlike the tetrahydrate of triflic acid (CF(3)SO(3)H·4H(2)O), there is no evidence of the proton preferentially transferring to a water molecule bridging two of the sulfonate groups. However, this could be an artifact of the crystal structure since the only such water molecule is separated from the proton by long O-O distances. Hydrogen bonding criteria, using the two-dimensional potential of mean force, are extracted. Radial distribution functions, free energy profiles, radii of gyration, and the root-mean-square displacement computed from ab initio path integral molecular dynamics simulations reveal that quantum effects do significantly extend the size of the protonic defect and increase the frequency of proton transfer events by nearly 15%. The calculated IR spectra confirm that the dominant protonic defect mostly exists as an Eigen cation but contains some Zundel ion characteristics. Chain lengths and ring sizes determined from the hydrogen bond network, counted using graph theory techniques, are only moderately sensitive to quantum effects. Deliberately introducing a structural defect into the native crystal yields a protonic defect with one hydrogen bond to a sulfonate group that was found to be metastable for at least 10 ps.

Mesoscale modeling of hydrated morphologies of 3M perfluorosulfonic acid-based fuel cell electrolytes

Dissipative particle dynamics (DPD) simulations have been carried out to study the hydrated morphology of 3M perfluorosulfonic acid (PFSA) fuel cell membranes as a function of the equivalent weight (EW), molecular weight (MW), and hydration level. The 3M PFSA ionomers were modeled using typical EWs of 578, 640, and 790 g/mol, and molecular weights were varied from about 45,000 to 90,000 g/mol in order to be close to the experimental range. The morphology changes corresponding to the EW, MW, and hydration level were comparatively investigated by inspecting the water distributions, followed by quantitative analysis by radial distribution functions and Bragg spacing according to the periodicity of water domains. Compared to the morphologies of short-side-chain PFSA membrane (Wu, D.; Paddison, S. J.; Elliott, J. A. Macromolecules 2009, 42, 3358-3367), the longer side chain in 3M PFSA membrane provides more flexibility for the sulfonate-terminated side chains and generally results in the stronger aggregation of water clusters. This results in lower water uptake for higher EW, corresponding to a lower ion-exchange capacity (IEC), which is attributed experimentally to a higher crystallinity of the fluorocarbon phase, although our simulations were not able to observe the crystallites directly.

Ab initio molecular dynamics simulations investigating proton transfer in perfluorosulfonic acid functionalized carbon nanotubes

Proton dissociation and transfer were examined with ab initio molecular dynamics (AIMD) simulations of carbon nanotubes (CNT) functionalized with perfluorosulfonic acid (-CF(2)SO(3)H) groups with 1-3 H(2)O/SO(3)H. The CNT systems were constructed both with and without fluorine atoms covalently bound to the walls to elucidate the effects of the presence of a strongly hydrophobic environment, the fluorine, on proton dissociation, hydration, and stabilization. The simulations revealed that the dissociated proton was preferentially stabilized as a hydrated hydronium cation (i.e., Eigen like) in the fluorinated CNTs but as a Zundel (H(5)O(2)(+)) cation in the nonfluorinated CNTs. This feature is attributed to the fluorine atoms forming hydrogen bonds with the water molecules coordinated to the central hydronium ion.

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

A simulation study was made of the effects of strong electric fields on the morphology of a Nafion-like ionomer at various levels of hydration. The results of united-atom molecular-dynamics computations showed a self-organization of the side chain terminal groups into cylindrical clusters. The walls of these clusters contain the sulfonate dipoles, while the interior holds the majority of the water molecules. These cylindrical structures then align to form an hexatic array aligned along the direction of the applied electric field. The hexatic morphology persists after the removal of the field. A calculation by means of the Kirkwood coupling parameter method shows the Helmholtz free energy of the hexatic morphology of the poled membrane to be lower than that of the initial isotropic material, even in the absence of the applied field.

Proton transport in triflic acid hydrates studied via path integral car-parrinello molecular dynamics

The mono-, di-, and tetrahydrates of trifluoromethanesulfonic acid, which contain characteristic H(3)O(+), H(5)O(2)(+), and H(9)O(4)(+) structures, provide model systems for understanding proton transport in materials with high perfluorosulfonic acid density such as perfluorosulfonic acid membranes commonly employed in hydrogen fuel cells. Ab initio molecular dynamics simulations indicate that protons in these solids are predisposed to transfer to the water most strongly bound to sulfonate groups via a Grotthuss-type mechanism, but quickly return to the most solvated defect structure either due to the lack of a nearby species to stabilize the new defect or a preference for the proton to be maximally hydrated. Path integral molecular dynamics of the mono- and dihydrate reveal significant quantum effects that facilitate proton transfer to the "presolvated" water or SO(3)(-) in the first solvation shell and increase the Zundel character of all the defects. These trends are quantified in free energy profiles for each bonding environment. Hydrogen bonding criteria for HOH-OH(2) and HOH-O(3)S are extracted from the two-dimensional potential of mean force. The quantum radial distribution function, radius of gyration, and root-mean-square displacement position correlation function show that the protonic charge is distributed over two or more water molecules. Metastable structural defects with one excess proton shared between two sulfonate groups and another Zundel or Eigen type cation defect are found for the mono- and dihydrate but not for the tetrahydrate crystal. Results for the tetrahydrate native crystal exhibit minor differences at 210 and 250 K. IR spectra are calculated for all native and stable defect structures. Graph theory techniques are used to characterize the chain lengths and ring sizes in the hydrogen bond network. Low conductivities when limited water is present may be attributable to trapping of protons between SO(3)(-) groups and the increased probability that protons transfer to waters bridging two different sulfonate groups.

Determining proton diffusion in polymer films by lifetimes of luminescent complexes measured in the frequency domain

Polymer-supported luminescent metal complexes represent an important class of oxygen, pH, and ion sensors. The diffusion properties of the analyte into the sensing film are important for rational sensor and support design and development. We describe a technique using lifetime measurements in the frequency domain for determining the diffusion coefficient of hydrochloric acid through various polymeric pH sensor films. Two types of polymers are doped with [Ru(4,7-diphenyl-1,10-phenanthroline)2(4,4'-dicarboxy-2,2'-bipyridine)]Cl2. We monitor the phase shift of luminescence (from which we calculate the apparent lifetime, tau(app)) versus time after applying a step increase in the aqueous HCl concentration at the surfaces of the film. We model the decrease in tau(app) as a function of time using the diffusion coefficient of HCl in the polymer as the only adjustable parameter. The model accurately predicts the lifetime versus time curves, and the resulting diffusion coefficients are highly dependent on the polymer. Relative to bulk water, diffusion of protons within very hydrophilic hydrated D4 polymer (a polyethylene oxide cross-linked siloxane ring polymer) films is hindered approximately 4-fold, while within a more hydrophobic sol gel it is hindered by over 1 order of magnitude. The methodology is adaptable for measuring diffusion coefficients of a variety of analytes in different sensor films as long as the bound and unbound forms luminescence and the excited states have different lifetimes.

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.