Stacking configurations, free energy fluctuation in the AIMD simulations at the temperature of 300 K, considered FM and AFM states, schemetic of the nearest and second-nearest exchange coupling parameters, radar charts of MAE in two polarization states.