“…Motivated by this, the study for the real-variable theory of various function spaces on R n or domains in R n , especially, the Hardy-type spaces, associated with different differential operators, has inspired great interests in recent years; see, for example, [3,4,31,32,34,46,49,52,73,78] for the case of Hardy spaces, [11,61,72] for the case of weighted Hardy spaces, [1,88,89,90] for the case of variable exponent Hardy spaces and [2,50,51,80,82,83,85] for the case of (Musielak-)Orlicz Hardy spaces. Recall that the Musielak-Orlicz space was originated by Nakano [67] and developed by Musielak and Orlicz [64,65], which is a natural generalization of many important spaces such as (weighted) Lebesgue spaces, variable Lebesgue spaces and Orlicz spaces and not only has its own interest, but is also very useful in partial differential equations [6,7,44,40], in calculus of variations [27], in image restoration [43,54] and in fluid dynamics [77,62]. The Musielak-Orlicz Hardy space on R n has proved useful in harmonic analysis (see, for example, [56,18,57,79]) and, especially, naturally appears in the endpoint estimate for both the div-curl lemma and the commutator of Cald...…”