Stability analysis of TOI-6041 system. The phase space of the system is explored by varying the orbital period, P_i, and eccentricity, e_i, of each planet independently. For each initial condition (Table 4), the system is integrated over 10 kyr, and a stability criterion is derived with the frequency analysis of the mean longitude. The chaotic diffusion is measured by the variation in the frequencies. The color-scale corresponds to the logarithmic variation in the stability index used in Correia et al. (2010). The red zone corresponds to highly unstable orbits, while the dark blue region can be assumed to be stable on a billion-year timescale. Dashed white lines indicate the nominal orbital periods reported in Table 4.