Boundary value problems for the Euclidean Hodge-Laplacian in three dimensions −Δ HL ∶= curl curl − grad div lead to variational formulations set in subspaces of H(curl, Ω) ∩ H(div, Ω), Ω ⊂ R 3 a bounded Lipschitz domain. Via a representation formula and Calderón identities, we derive corresponding first-kind boundary integral equations set in trace spaces of H 1 (Ω), H(curl, Ω), and H(div, Ω). They give rise to saddle-point variational formulations and feature kernels whose dimensions are linked to fundamental topological invariants of Ω. Kernels of the same dimensions also arise for the linear systems generated by low-order conforming Galerkin (BE) discretization. On their complements, we can prove stability of the discretized problems, nevertheless. We prove that discretization does not affect the dimensions of the kernels and also illustrate this fact by numerical tests. KEYWORDS boundary element method (BEM), Calderón indentities, first-kind boundary integral equations, harmonic vector fields, Hodge-Laplacian, representation formula, saddle-point problems MSC CLASSIFICATION 31A10; 45A05; 65N38