Molecular Simulation of Nafion Swelling in the Cathode Catalyst Layer of Polymer Electrolyte Fuel Cells
Authors: Koda Naoki , Sato Kohei , Otic Clint John Cortes , Kinefuchi Ikuya 
Published: February 2024

As the world moves towards carbon neutrality, polymer electrolyte fuel cells (PEFCs) have emerged as a vital technology in recent years. The cathode catalyst layer of PEFCs is characterized by its porous structure, which promotes the production of water molecules from oxygen and protons through chemical reactions. The ionomer thin films, predominantly composed of Nafion, are critical in regulating the transport and mobility of molecules and ions in the cathode catalyst layer. Sulfonic groups in the side chains of Nafion enable accumulation of water molecules, causing the ionomer thin films to swell.
Previous studies have examined the swelling of Nafion under wet conditions[1][2]. However, precise swelling mechanisms remain unclear. While privious research has analyzed the relationship between relative humidity and the amount of water within bulk Nafion samples[3], the swelling behavior in films with thickness of several nanometers differs significantly, necessitating further investigation[4]. This study therefore investigates the swelling mechanisms and the properties of ionomer thin films used in the cathode catalyst layer of PEFCs. This study aims to establish guidelines for designing the molecular structures of ionomers by elucidating the mechanisms underlying ionomer swelling in the cathode catalyst layer. Grand canonical Monte Carlo (GCMC) simulations and molecular dynamics (MD) simulations were employed to reproduce swelling in Nafion thin films and examine the impact of Nafion's structure on water content and swelling. These molecular simulation techniques enable the elucidation of the influence exerted by polymer structure on water content and swelling at the molecular scale.
In this study, we firstly created a simulation model of a Nafion system using LAMMPS[5], the DREIDING Force Field[6], and Mulliken charges[7] in which parameter settings were referred from literature values[8]. The Nafion system contains 10 Nafion molecules with a degree of polymerization of 10, water molecules, and 100 oxonium ions for electrical neutrality. The density of 1.9 g/cm3 achieved during simulations closely agrees with literature values[8].
To reproduce Nafion swelling under humid conditions, GCMC simulations and MD simulations were employed. The rationale behind utilizing GCMC simulations lies in the fact that they enable the insertion of water molecules into the Nafion system under constant chemical potential conditions, thereby inducing Nafion swelling. Subsequent to the GCMC simulations, MD simulations were performed under the NPT ensemble, which facilitated the system to reach equilibrium. This alternating cycle of GCMC and MD simulations is repeated to realize the final swelling state. Preliminary results indicate lower water content values than those experimentally measured under similar conditions. Possible causes include inadequate force field parameters that lead to polymer morphology different from the actual Nafion structure.
Future work will explore the causes of these discrepancies and examine the effects of Nafion's structure on water content and swelling under various conditions. This will involve testing alternative force field parameters and initial structures to better represent the molecular-level structure of ionomer thin films on carbon support particles in the cathode catalyst. Additionally, we will investigate how different relative humidity conditions, temperatures, and ionomer compositions influence swelling behavior, which could provide insights into optimizing PEFCs performance under various operating environments.
[1] G. C. Abuin et al., J. Membr. Sci.
428, 507 (2013).
[2] J. Catalano et al., Int. J. Hydrog. Energy
37, 6308 (2012).
[3] T. A. Zawadzinski et al., J. Electrochem. Soc.
140, 1981 (1993).
[4] M. A. Modestino et al., Macromolecules
46, 867 (2013).
[5] A. P. Thompson et al., Comp. Phys. Comm.
271, 10817 (2022).
[6] S. L. Mayo et al., J. Phys. Chem.
94, 8897 (1990).
[7] S. S. Jang et al., J. Phys. Chem.
B
108, 3149 (2004).
[8] Y. Kurihara et al., J. Power Sources
414, 263 (2019).
Fig. 1 An ionomer model composed of Nafion molecules, oxonium ions, and water molecules. The red and white spheres represent oxygen and hydrogen atoms making up water molecules. The size of the water molecules is emphasized for clarity.

Figure 1