In this talk, we will discuss existence, regularity, and qualitative properties of solutions to the Hamiltonian elliptic system $$ -\Delta u = |v|^{q-1} v\ \ \ \text{in} \ \Omega,\quad -\Delta v = |u|^{p-1} u\ \ \ \text{in} \ \Omega,\quad \partial_\nu u=\partial_\nu v=0\ \ \ \text{on} \ \partial\Omega,$$with $\Omega\subset \mathbb R^N$ bounded, both in the sublinear $pq< 1$ and superlinear $pq>1$ problems, in the subcritical regime. In balls and annuli, we show that least energy solutions are not radial functions, but only partially symmetric (namely foliated Schwarz symmetric). A key element in the proof is a new $L^t$-norm-preserving transformation, which combines a suitable flipping with a decreasing rearrangement. This combination allows us to treat annular domains, sign-changing functions, and Neumann problems, which are nonstandard settings to use rearrangements and symmetrizations. Our theorems also apply to the scalar associated model, where our approach provides new results as well as alternative proofs of known facts.