“…Considering the unitary transformation U 1 = exp(−iθσ y /2) with θ = arctan(Ω s /∆ s ) that intends to diagonalizing H 0 in the dressed qubit basis, (i.e., the dressed basis is defined as the eigenstates of H 0 ) the following dressed system HamiltonianH dress = U † 1 H R U 1 , (i.e., the Hamiltonian in the dressed basis)H dress =[∆ − Ω d sin θ cos(∆ d t + φ)] Z Ω d cos θ cos(∆ d t + φ) X Ω d sin(∆ d t + φ)Y denotes the microwave-dressed qubit detuning, and {X = cos θσ x − sin θσ z , Y = σ y , Z = cos θσ z + sin θσ x } represent the Pauli operators defined on the dressed basis. According to the following unitary transformation[35,59,60]U 2 = exp −i Z 2 ∆t − Ω d sin θ ∆ d sin(∆ d t + φ) , (A5)and using the Jacobi-Anger relations, one can obtain the effective HamiltonianH J = U † 2 H dress U 2 + i∂ t (U † 2 )U 2 , i.e., H J = Ω d cos θ 4 (e i∆ d t+iφ + e −i∆ d t−iφ ) × e i∆t S + ∞ n=−∞ J n ( Ω d sin θ ∆ d )e −in(∆ d t+φ) + h.c. − Ω d 4 (e i∆ d t+iφ − e −i∆ d t−iφ ) × e i∆t S + ∞ n=−∞ J n ( Ω d sin θ ∆ d )e −in(∆ d t+φ) − h.c. ,(A6)where h.c. denotes the Hermitian conjugate, J n is the nth order Bessel function of the first kind, and S ± = (X ± iY )/2.Applying the RWA and dropping high-order Bessel functions, one can obtain the following effective Hamiltonian describing the usual single-qubit driven terms,H eff = Ω J 0 ( Ω d sin θ ∆ d ) e −iφ S + e i(∆−∆ d ) + h.c. .…”