In this paper, we present a novel approach based on Bayesian Expectation Maximization-Maximization (BEMM) to address this challenge. Unlike traditional optimization methods, which may struggle with high-dimensional and nonlinear optimization problems, BEMM offers a robust framework that combines the benefits of Bayesian inference with the flexibility of expectation maximization and maximization techniques. By iteratively updating parameter estimates based on observed data and maximizing the likelihood of the model, BEMM effectively navigates the solution space to converge on accurate estimates of the unspecified variables in Proton Exchange Membrane Fuel Cells (PEMFCs)models. Through extensive experimentation and comparison with other metaheuristic techniques, including Arithmetic Optimization Algorithm (AOA), Gravitational Search Algorithm (GSA), Flower Pollination Algorithm (FPA), and Biogeography-Based Optimization (BBO), we demonstrate the superior performance of our BEMM approach. Our results show that BEMM outperforms these alternative methods in terms of precision, convergence speed, and stability across various scenarios involving different numbers of unspecified variables. The implications of our findings are significant for both researchers and practitioners in the field of PEMFC modeling and optimization. By providing a more efficient and reliable method for estimating model parameters, our approach can facilitate more accurate predictions of PEMFC performance, leading to better-informed decision-making in the design, operation, and optimization of PEMFC systems. Furthermore, the robustness and versatility of BEMM make it well-suited for a wide range of optimization problems beyond PEMFC modeling, highlighting its potential impact across various domains of engineering and scientific research. In the sensitivity analysis, as the population size increases from 10 to 40, there is a significant improvement in solution quality by approximately 100 %. However, beyond a population size of 40, the marginal gains diminish, with only marginal improvements of less than 1 % observed despite further increases in population size up to 200.