The functionality of weakly bound organic materials, either in Nanoelectronics or in Materials Science, is known to be strongly affected by their morphology. Theoretical predictions of the underlying structure-property relationships are frequently based on calculations performed on isolated dimers, but the optimized structure of the latter may significantly differ from experimental data even when dispersion-corrected methods are used for it. Here, we address this problem on two organic crystals, namely coronene and 5,6,11,12-tetrachlorotetracene, concluding that it is caused by the absence of the surrounding monomers present in the crystal, and that it can be efficiently cured when the dimer is embedded into a general Quantum Mechanics/Molecular Mechanics (QM/MM) geometry optimization scheme. We also investigate how the size of the MM region affects the results. These findings may be helpful for the simulation of the morphology of active materials in crystalline or glassy samples.