One way to test our hypothesis that both intrinsic proton hopping and intrinsic OH rotation combine to produce long-range PT is to independently vary the values of the proton hopping barrier, E hop , and the OH rotation barrier, E rot , and then observe how the value of D s responds to these changes. It is impossible to use DFT or DP simulations for this. Hence, we developed a lattice Monte Carlo (LMC) model of PT on graphanol, which has E hop and E rot as parameters. Details of the LMC simulations are given in the Methodology section and details of the algorithm are presented in the Supporting Information. We systematically varied E hop and E rot from 52 meV to 518 meV and calculated D s to perform a sensitivity analysis. These calculations were performed at a constant temperature of T = 400 K. The results are tabulated in Figure . We found that both E hop and E rot impact the value of D s , with E hop having a larger influence on D s than E rot . This can be seen, for example, by noting that for E hop = 52 meV and E rot = 207 meV we get D s = 2.8 Å2 /ps, whereas for E hop = 207 meV and E rot = 52 meV we get D s = 0.6 Å2 /ps, which is a factor of 4.7 times smaller. Hence, the impact of hopping and OH rotation is asymmetric. We next computed E A for two different barrier combinations: E hop = 52, E rot = 259 meV and E hop = 259, E rot = 52 meV, by calculating D s as a function of temperature and fitting the data to the Arrhenius equation. These results are plotted in Figure . The overall E A appears to be more sensitive to E hop than E rot , as one might expect based on the D s sensitivity analysis.