Abstract A flag domain in ℝ 3 {\mathbb{R}^{3}} is a subset of ℝ 3 {\mathbb{R}^{3}} of the form { ( x , y , t ) : y < A ( x ) } {\{(x,y,t):y<A(x)\}} , where A : ℝ → ℝ {A\colon\mathbb{R}\to\mathbb{R}} is a Lipschitz function. We solve the Dirichlet and Neumann problems for the sub-elliptic Kohn–Laplacian △ ♭ = X 2 + Y 2 {\bigtriangleup^{\flat}=X^{2}+Y^{2}} in flag domains Ω ⊂ ℝ 3 {\Omega\subset\mathbb{R}^{3}} , with L 2 {L^{2}} -boundary values. We also obtain improved regularity for solutions to the Dirichlet problem if the boundary values have first order L 2 {L^{2}} -Sobolev regularity. Our solutions are obtained as sub-elliptic single and double layer potentials, which are best viewed as integral operators on the first Heisenberg group. We develop the theory of these operators on flag domains, and their boundaries.