Molecular Dynamics Analysis of the Scattering Phenomena of Oxygen Molecules on an Ionomer Surface in Catalyst Layer of Fuel Cell
Authors: Mizuki Keisuke , Mabuchi Takuya , Kinefuchi Ikuya , Tokumasu Takashi 
Published: February 2024

As energy demand increases and global warming progresses, expectations for fuel cells, which generate electricity through chemical reactions between hydrogen and oxygen, are rising. There are various types of fuel cells, which are classified according to the electrolyte. Polymer electrolyte fuel cells (PEFCs) are expected to be used in stationary power sources for home use and fuel cell vehicles because of their low operating temperature, short start-up time, and ease of miniaturization. The cost has been an issue in the widespread use of fuel cells. The amount of platinum used in the catalyst layer (CL) must be reduced to reduce costs. However, it increases the current density per unit area of platinum. There are three factors that contribute to voltage drop: ohmic polarization, activation polarization, and diffusion polarization. At high current densities, the voltage drop due to diffusion polarization is dominant. The main cause of diffusion polarization is the transport resistance of oxygen molecules in the CL, and improving the oxygen transport properties in the CL will lead to lower cost PEFCs. In this study, we analyzed the transport properties of oxygen molecules in the CLs of PEFCs using the Monte Carlo (MC) method. The effect of oxygen scattering on the overall oxygen transport properties was investigated by considering the behavior on the surface of the ionomer film (surface diffusion) into the MC method, which simulates the overall transport in the PEFC CL. The objective of this study is to reproduce accurate oxygen transport by comparing the results of MC and molecular dynamics (MD) methods. First, we used the MC method to analyze oxygen scattering phenomena. In this study, we used the 3D structure data of CL in the MC simulations. The three-dimensional structure of the sample was reconstructed from sequential slice images. Molecules were assumed to reflect specularly on the boundaries of the simulation domain. Oxygen molecules were randomly placed in the simulation system. The time step was set at 5 ps and the simulation time was set at 4 μs. In order to examine the effect of surface diffusion on the transport properties of oxygen molecules, the simulation was performed using various conditions of surface diffusion. We calculated the effective diffusion coefficient of oxygen in the CL as the transport property of oxygen molecules in the MC method. As the surface diffusion coefficient increases, the effective diffusion coefficient of oxygen increases. When the surface diffusion coefficient is small, the effective diffusion coefficient of oxygen increases due to a decrease in the time constant, which determines the probability distribution of the surface residence time. However, when the surface diffusion coefficient is large, a different trend emerges. As the time constant increases, the effective diffusion coefficient has a peak value and then decreases. Next, to verify the model of surface diffusion in the MC method, we analyzed the scattering phenomenon of oxygen molecules in a simulation system that simulates the pores in the CL using the MD method. We modeled the nano pores in the CL as a slit between the ionomer walls in the MD method. Nafion with the equivalent weight (EW) of 1146 was used as polymer model. The number of Nafion chains is set at 4, and the water content λ is set at 7. Three-dimensional periodic boundary conditions were applied to the simulation domain. After the system was equilibrated, 1 ns of NVT simulation was performed for production. The temperature was set at 297 K. The behavior of oxygen molecules impacting on the ionomer wall was analyzed. We analyzed the behavior of oxygen molecules colliding with the ionomer membrane using the MD method. The average residence time of surface diffusion on the ionomer thin film was of the order of 10-10s. When the calculation results were compared with the MC results, it was found that the average residence time obtained by the MD method was the value where the surface diffusion coefficient affects much on the effective diffusion coefficient. Comparing the results of the MD method with the surface diffusion model used in the MC method, there are some discrepancies between the results and the model. For example, the MC method used a model in which oxygen molecules reflect diffusively on the ionomer surface. However, in the MD method, it was found that oxygen molecules tend to reflect in the direction of travel (Fig. 1). Therefore, a more accurate model of the MC method needs to be developed in the future.

Figure 1