Porovnáváme čtyři různé typy pohybových rovnic pro matici hustoty systému molekulárních excitonů interagujících s termodynamickou lázní. Diskutujeme vhodnost těchto metod k simulaci přenosu energie v systému a identifikujeme vlastnosti dvojdimenzionálních koherentní spekter, které umožňují odlišit časově nelokální efekty od lokálních.