Boundary-term shift of canonical field momenta

If the spatially integrated Lagrangian changes by $dB[\Phi ]/dt$, where $B$ depends on fields and their spatial derivatives but not velocities, the canonical momentum components change by ${\pi }_a\mapsto {\pi }_a+\delta B/\delta {\Phi }_a$. This follows by varying $dB/dt$ and using spatial integration by parts. The canonical one-form changes by the exact variation $\delta B$, so the symplectic form and canonical commutation relations are preserved. For the two Feynman gauge Maxwell densities, $B=\int A_0{\partial }_i{A}_i{d}^{3}x$ gives shifts $\Delta {\pi }^0={\partial }_i{A}_i$ and $\Delta {\pi }^i=-{\partial }_i{A}_0$. Boundary flux must vanish for this spatial functional calculation.