Atomistic molecular dynamics simulations have been performed on heptyl phosphonic acid (HPA) to understand the dynamic hydrogen bonding network in the liquid phase. HPA is a phosphonic-acid functionalized alkane (heptane) and a model compound for one of the promising polymers for high temperature (> 100 degrees C) fuel cell polymer electrolyte membranes. For the simulation, a force field for this molecule has been generated with the help of quantum chemical calculations and optimized by simplex algorithm. The force field has been validated against experimentally measured properties, for example, density and self-diffusion constant. From molecular dynamics simulations conducted at different temperatures, we have confirmed the hypothesis of dynamic hydrogen bond network formation in this material.