The MD calculations of the dispersed disordered model for Form I showed several changes in the hydrogen bonding motif during the simulation and this supports the case for this model over the domain disordered one. Attempts to further support this claim by comparing the calculated IR spectra from these models to experiment are frustrating. The predictions for both models are very similar above 150 cm -1 and all compare well with experiment, suggesting that the mid-IR region appears to be insensitive to the disorder and even the S-H stretch region is difficult to interpret as the broad experimental peak widths mask the effect of the different hydrogen bonding environments.