On-the-fly kinetic Monte Carlo (KMC) simulations are performed to investigate slow relaxation of non-equilibrium systems. Point defects induced by 25 keV cascades in α-Fe are shown to lead to a characteristic time-evolution, described by the replenish and relax mechanism. Then, we produce an atomistically-based assessment of models proposed to explain the slow structural relaxation by focusing on the aggregation of 50 vacancies and 25 self-interstital atoms (SIA) in 10-latticeparameter α-Fe boxes, two processes that are closely related to cascade annealing and exhibit similar time signature. Four atomistic effects explain the timescales involved in the evolution: defect concentration heterogeneities, concentration-enhanced mobility, cluster-size dependent bond energies and defect-induced pressure. These findings suggest that the two main classes of models to explain slow structural relaxation, the Eyring model and the Gibbs model, both play a role to limit the rate of relaxation of these simple point-defect systems.Many off-equilibrium physical systems exhibit a slow structural relaxation toward their ground state. Examples include glasses [1,2], colloids [3], concrete [4] and amorphous solids [5,6]. The degree of relaxation of these systems, e.g. their potential energy, is a nearly linear function of the logarithm of time. For simplicity, the term logarithmic relaxation is used to describe this behavior.Most models describe such aging as a sequence of activated processes that permit relaxation. As the system relaxes, the energy barriers of these processes increase , which delays aging by growing orders of magnitude in time. The literature contains a large number of such propositions. For example, in the Eyring model [7], barriers are linked to the energy recovered during relaxation through their coupling to shear strain. This idea, where the degree of relaxation has an effect on the height of the barriers, has been adapted and modified to include hierarchically constrained dynamics [8], stress relaxation in the Burridge-Knopoff model [9], models with a stressinduced barrier increase [6,10] and glassy polymer relaxation [1, 11] after nano-indentation. On the other hand, the Gibbs model [12] and its variants [13][14][15] argue that the system initially possesses a distribution of relaxation events with a near-constant density as a function of activation barrier, or rates described by a multiplicative stochastic process [2], which leads to logarithmic relaxation. A newly proposed model links logarithmic timeevolution to the system moving from one local state to another, where the waiting time of each state is defined by a power law and where all states evolve simultaneously [16]. Due to a lack of atomistic evidence, the validity of many of these proposed models remains ambiguous and the identity of the drivers to logarithmic relaxation remains elusive.Recently, we proposed a novel description of logarithmic relaxation, coined replenish and relax [17], which was based on nanocalorimetric measurements combined with ...