We consider a notion of conservation for the heat semigroup associated with a generalized Dirac Laplacian acting on sections of a vector bundle over a noncompact manifold with a (possibly noncompact) boundary under mixed boundary conditions. Assuming that the geometry of the underlying manifold is controlled in a suitable way and imposing uniform lower bounds on the zero-order piece (Weitzenbock potential) of the Dirac Laplacian, and on the endomorphism defining the mixed boundary condition, we show that the corresponding conservation principle holds. A key ingredient in the proof is a domination property for the heat semigroup which follows from an extension to this setting of a Feynman–Kac formula recently proved by the author de Lima (Pac J Math 292(1):177–201, 2018) in the context of differential forms. When applied to the Hodge Laplacian acting on differential forms satisfying absolute boundary conditions, this extends previous results by Vesentini (Ann Math Pura Appl 182:1–19, 2003) and Masamune (Atti Accad Naz Lincei Rend Lincei Mat Appl 18(4):351–358, 2007) in the boundaryless case. Along the way, we also prove a vanishing result for $$L^2$$ harmonic sections in the broader context of generalized (not necessarily Dirac) Laplacians. These results are further illustrated with applications to the Dirac Laplacian acting on spinors and to the Jacobi operator acting on sections of the normal bundle of a free boundary minimal immersion.