The new technique of nuclear resonance vibrational spectroscopy (NRVS) has increased the range and quality of dynamical data from Fe-containing molecules that when combined with Raman and infrared spectroscopies impose stricter constraints on normal mode simulations, especially at lower frequencies. Going beyond the usual single molecule approximation, a classical normal-mode analysis that includes intermolecular coupling and the full crystalline symmetry is found to produce a better fit with fewer free parameters for the heme compound iron octaethylporphyrin (FeOEP), using NRVS data from polycrystalline material. Off-diagonal force constants were completely unnecessary, indicating that their role in previous single molecule fits was just to emulate intermolecular coupling. Sound velocities deduced from the calculated phonon dispersion curves are compared to NRVS measurements to further constrain the intermolecular force constants. The NRVS data by themselves are insufficient to rigorously determine all unknown force constants for molecules of this size, but the improved crystal model fit indicates the necessity of including intermolecular interactions for normal-mode analyses.