In this paper we establish the product Hardy spaces associated with the Bessel Schrödinger operator introduced by Muckenhoupt and Stein, and provide equivalent characterizations in terms of the Bessel Riesz transforms, non-tangential and radial maximal functions, and Littlewood-Paley theory, which are consistent with the classical product Hardy space theory developed by Chang and Fefferman. Moreover, in this specific setting, we also provide another characterization via the Telyakovskií transform, which further implies that the product Hardy space associated with this Bessel Schrödinger operator is isomorphic to the subspace of suitable "odd functions" in the standard Chang-Fefferman product Hardy space. Based on the characterizations of these product Hardy spaces, we study the boundedness of the iterated commutator of the Bessel Riesz transforms and functions in the product BMO space associated with Bessel Schrödinger operator. We show that this iterated commutator is bounded above, but does not have a lower bound. √ S λ or T t := e −tS λ .We now define the product Hardy space via the Littlewood-Paley area functions as follows.We now define another version of the Littlewood-Paley area function. Let ∇ t 1 , y 1 := (∂ t 1 , ∂ y 1 ), ∇ t 2 , y 2 := (∂ t 2 , ∂ y 2 ).
Then the Littlewood-Paley area function∇ t 1 , y 1 e −t 1 √ S λ ∇ t 2 , y 2 e −t 2 √ S λ (f )(y 1 , y 2 ) 2 dy 1 dy 2 dt 1 dt 2 1 2. (1.4)