“…Moreover, the weighted L p (R n )-boundedness of these operators was studied in [37]. Recently, Ly [27] proved that V L −1 and ∇ 2 L −1 are bounded from the Hardy space H p L (R n ), associated with L, to L p (R n ) when p ∈ (0, 1], and ∇ 2 L −1 is also bounded from H p L (R n ) to the classical Hardy space H p (R n ) when p ∈ ( n n+1 , 1], via the boundedness of ∇ 2 L −1 on L p (R n ) with some p ∈ (1, ∞) and some Sobolev type estimates for the heat kernel of L. Moreover, the boundedness of ∇ 2 L −1 and V L −1 on the Musielak-Orlicz-Hardy space H ϕ, L (R n ) was independently studied in [7] by different method. Very recently, the boundedness of V L −1 , V 1/2 (∇ − i a)L −1 and (∇ − i a) 2 L −1 from the Musielak-Orlicz-Hardy space H ϕ, L (R n ), associated with the magnetic Schrödinger operator L, to the Musielak-Orlicz space L ϕ (R n ) was established in [8], where L := −(∇ − i a) • (∇ − i a) + V with a := (a 1 , a 2 , .…”