  This paper introduces an ultra-weak space-time DPG method for the heat
equation. We prove well-posedness of the variational formulation with broken
test functions and verify quasi-optimality of a practical DPG scheme. Numerical
experiments visualize beneficial properties of an adaptive and parabolically
scaled mesh-refinement driven by the built-in error control of the DPG method.
