L boundedness of commutator operator associated with schrödinger operators on heisenberg group

Let L = −ΔHⁿ + V be a Schrödinger operator on Heisenberg group Hⁿ, where ΔHⁿ is the sublaplacian and the nonnegative potential V belongs to the reverse Hölder class B_Q/2, where Q is the homogeneous dimension of Hⁿ  . Let T₁=(−ΔHn+V)⁻¹V,T₂=(−ΔHn+V)^−1/V₂1/2$T_1={(-\Delta H^n+V)}^{-1}V,T_2={(-\Delta H^n+V)}^{-1/2_{V}1/2}$, and T₃=(−Δ_Hn+V)^−1/2_∇Hn$T_3={(-{\Delta }_{H^n}+V)}^{-1/2}{\nabla }_{H^n}$, then we verify that b, T_i], i = 1,2,3 are bounded on some L^p(Hⁿ), where b ∈ BMO(Hⁿ). Note that the kernel of T_i, i = 1,2,3 has no smoothness.