The Landau–de Gennes Q-tensor theory provides a fundamental framework for modeling nematic liquid crystals. Numerical simulations of the corresponding Q-tensor gradient flows are challenging due to the presence of stiff nonlinear terms and the need to preserve essential physical constraints, in particular the energy dissipation law and the maximum bound principle. In this talk, I will present a class of first- and second-order stabilized exponential time differencing schemes for the Q-tensor flow. By introducing a linear stabilizing parameter that is consistent with the Lipschitz constant of the nonlinear bulk potential, we develop fully discrete linear schemes in which the tensor components are completely decoupled. We rigorously prove that the proposed schemes unconditionally preserve the discrete MBP and the original energy dissipation law. Moreover, optimal-order error estimates are established for both schemes. Numerical experiments will be presented to validate the theoretical results, demonstrating that the proposed schemes achieve high-order temporal accuracy while strictly preserving the physical bounds and energy dissipation property, even for relatively large time step sizes.