Proton transfer via tunneling in hydrogen bonds is a phenomenon that occurs in many chemical compounds and biological systems. It explains several molecular properties, such as the interconversion of tautomeric forms of nucleic acids-base pairs and the conformational dynamics of proteins. Because of the complexity of these systems, it is often challenging to formulate a mathematical model capable of explaining the mechanism of proton tunneling. This study proposes an approach based on path decomposition expansion, which allows the modelling of the proton transfer reaction in multidimensional space, simplifying the investigation of complex phenomena such as the domino proton transfer or competing proton transfer mechanisms.