Методом молекулярної динаміки з використанням програмного забезпечення LAMMPS досліджено процеси релаксації та реконструкції поверхневих шарів Fe (001), (011), (111) до та після нанесення графену для двох температур — 300 і 400 К. Варіювався також кут спряження кристалічних ґратниць графену та Fe. Досліджено 12 модельних систем, для яких розраховувались: повна, кінетична та потенціальна енергії системи; зміна міжплощинних відстаней вздовж напряму нормалі до поверхні (релаксація) відносно незбуреного стану (об’єм); розподіли атомів Заліза та Вуглецю вздовж нормалі до поверхні; пошарові функції радіального розподілу атомів С та Fe. Позиції атомів С та Fe, а також характер утвореного поверхневого рельєфу для різних систем Fe/графен, візуалізовано за допомогою 2$D$-мап атомного розподілу.