Here, we report extended ensemble molecular dynamics simulations of ammonia-cellulose I complex crystal models to evaluate the diffusion behavior of the guest ammonia molecules and the potential of mean force (PMF), namely, the free energy change along the chosen reaction coordinate, for migration of an ammonia molecule in the crystal models. Accelerated molecular dynamics simulations confirmed that ammonia molecules almost exclusively diffused through the hydrophilic channel even when the crystal framework was retained. Adaptive steered molecular dynamics simulations detected distinct PMF peaks with heights of approximately 7 kcal/mol as the ammonia molecule passed through the cellulose-chain layers. Introducing hybrid quantum mechanical and molecular mechanics theory to the adaptive steered molecular dynamics simulation effectively lowered the heights of the PMF peaks to approximately 5 kcal/mol, accompanied by a slight decrease in the baseline. Removal of the ammonia molecules in the neighboring channels resulted in a continuous increase in the baseline for the migration of an ammonia molecule in the hydrophilic channel. When the halves of the crystal model were separated to widen the hydrophilic channel to 0.2 nm, the PMF profiles exhibited an unexpected increase. This resulted from water structuring in the expanded hydrophilic channel, which disappeared with further expansion of the hydrophilic channel to 0.3 nm.