Nonconvex Carnot Limit Shapes for Branching Random Walks on the Heisenberg Group

We study a supercritical branching random walk on the discrete Heisenberg group driven by the four standard horizontal generators. Horizontal position has scale $n$, whereas the central coordinate records discrete area and has scale $n^2$. On the Galton-Watson survival event, we prove almost sure Hausdorff convergence, under Carnot scaling, of both the rescaled particle positions and the collections of enhanced ancestral paths. The path limit consists of area lifts of horizontal paths $\omega$ satisfying $J_t(\omega )\le t\log \rho$ for every $t\in [0,1]$, where $\rho$ is the offspring mean. In contrast to the classical Euclidean case, the deterministic limit shape is nonconvex for every $\rho >1$. We identify its horizontal projection and interval vertical fibers, determine the saturated shape for $\rho \ge 4$, and characterize the origin fiber below saturation by Wulff duality. A hyperbolic representation gives an elliptic integral formula for the Wulff area. The maximal central coordinate ${\mathrm{\Gamma }}_{\rho }$ is nondecreasing and concave in $\rho$, equals $1\mathrm{/}4$ for $\rho \ge 4$, and satisfies ${\mathrm{\Gamma }}_{\rho }=(\log \rho )\mathrm{/}\pi -(\log \rho {)}^2\mathrm{/}(4\pi )+O((\log \rho {)}^3)$ as $\rho \downarrow 1$. For fixed $1<\rho <4$, we obtain $M_n^Z={\mathrm{\Gamma }}_{\rho }{n}^2+O(n^{3\mathrm{/}2}\sqrt{\log n})$ almost surely. Above saturation the first correction is a random linear term; at $\rho =4$, under a finite second moment assumption, the deficit ${\mathrm{\Gamma }}_4{n}^2-M_n^Z$ has order $n\log \log n$ in probability.