The aim of this article is threefold. Firstly, we study Hardy spaces, hpL(Ω), associated with an operator L which is either the Dirichlet Laplacian ∆D or the Neumann Laplacian ∆N on a bounded Lipschitz domain Ω in Rn, for 0 < p ≤ 1. We obtain equivalent characterizations of these function spaces in terms of maximal functions and atomic decompositions. Secondly, we establish regularity results for the Green operators, regarded as the inverse of the Dirichlet and Neumann Laplacians, in the context of Hardy spaces associated with these operators on a bounded semiconvex domain Ω in Rn. Thirdly, we study relations between the Hardy spaces associated with operators and the standard Hardy spaces hr(Ω) and h p z(Ω), then establish regularity of the Green operators for the Dirichlet problem on a bounded semiconvex domain Ω in Rn, and for the Neumann problem on a bounded convex domain Ω in Rn, in the context of the standard Hardy spaces hr(Ω) and hz(Ω). This gives a new solution to the conjecture made by D.-C. Chang, S. Krantz and E.M. Stein regarding the regularity of Green operators for the Dirichlet and Neumann problems on hr(Ω) and h p z(Ω), respectively, for all n n+1 < p ≤ 1.