Let L=-div(A) be a second-order divergence form elliptic operator, where A is an accretive n×n matrix with bounded measurable complex coefficients in Rn. In this paper, we mainly establish the Lp boundedness for the commutators generated by bϵIα(BMO) and the square function related to fractional differentiation for second-order ellipt...