Kiser, AyasGurdal, Yeliz2026-02-272026-02-2720250360-31991879-348710.1016/j.ijhydene.2025.01.285http://dx.doi.org/10.1016/j.ijhydene.2025.01.285https://hdl.handle.net/20.500.14669/4352The type and coordination of ligands surrounding water reduction photo-catalysts affect the H2 production reaction pathway. In this respect, we examined two successive reductions followed by two protonations mechanism of H2 evolution reaction catalyzed by Cobalt-based pentapyridyl molecular catalyst in an aqueous solution employing Ab-initio Molecular Dynamics (AIMD), Density Functional Theory (DFT), and Free Energy Perturbation Theory. Each intermediate step of the H2 evolution reaction mechanism was simulated to determine the allowable spin states and solvent response surrounding the reaction center of the catalyst. A single electron transfer to the catalyst, leading to a singlet spin state, induces a conformational change in the water molecule in the first solvation shell, facilitating the formation of a Co-proton bond. Subsequent to the second electron transfer, the Co center acquires the proton from this water molecule, while the remaining OH- ion diffuses swiftly through the solvent. The first protonation step was determined to continue adopting a quartet spin state. Due to the catalyst's cage-like structure around the Co center, following the second protonation, H2 forms; however, it fails to diffuse from the reaction center into the solvent, suggesting the importance of maintaining an open structure around the Co center when designing ligand structures for water- splitting catalysts. The first and the second reduction free energies were calculated as -2 eV and -1.3 eV, respectively.eninfo:eu-repo/semantics/closedAccessWater splittingHydrogen evolution reactionAb-initio Molecular DynamicsDensity Functional TheoryReduction free energyArticle873864106WOS:001424854600001