Holomorphic functions in several complex variables were generalized to regular functions in several quaternionic variables, and further to monogenic functions of several vector variables, which are annihilated by several Dirac operators on k copies of the Euclidean space \({\mathbb {R}}^n\) . As the Dolbeault complex in complex analysis, the Dirac complex resolving several Dirac operators plays the fundamental role to investigate monogenic functions. Although the spaces in the Dirac complex are complicated irreducible modules of \(\textrm{GL}(k)\times \textrm{Spin}(n),\) we give a simple characterization of the first four spaces, which allows us to write down first three operators in the Dirac complex explicitly and to show this part to be an elliptic complex. Then the PDE method can be applied to obtain solutions to the non-homogeneous several Dirac equations under the compatibility condition, which implies the Hartogs’ phenomenon for monogenic functions. Moreover, we find the boundary version of several Dirac operators and introduce the notion of a tangentially monogenic function, corresponding to tangential Cauchy-Riemann operator and CR functions in several complex variables, and establish the Hartogs-Bochner extension for tangentially monogenic functions on the boundary of a domain.
