On a shape derivative formula for the Robin $p$-Laplace eigenvalue

We obtain shape derivative formulae for the first eigenvalue of the Robin $p$-Laplace operator. This result is used to study the variation of the first eigenvalue with respect to perturbations of the domain. In particular, we prove that for large values of the boundary parameter, the first eigenvalue is monotonic with respect to domain inclusion for smooth domains.