We have used the kinetic Monte Carlo simulation technique to analyze magnetic relaxation. We have used the same algorithm in the present work, described in detail in the references [1]}, [2]}, [3]}. Therefore, we do not reiterate it to avoid repetitions. In this procedure, we first saturate all the magnetic moments along the \(y\) -direction of the system by applying a huge magnetic field of strength \(\mu _oH_{\mathrm {o}} =20\)  T. Next, we divide the total simulation time into 2000 equal steps and switch off the external field \(\mu _oH_{\mathrm {o}}\)  at \(t=0\)  s. We then study the time evolution of magnetization of the underlying system using the kMC simulation. Finally, we fit the so-obtained magnetization-decay curve with \(M(t) = M^{}_s \exp (-t/\tau ^{}_N)\)  to extract effective Néel relaxation time \(\tau ^{}_N\)  of the underlying system.
