Consider this differential operator $$ \mathcal{H}(\phi(\mathbf{x})) = -\triangle + V(\mathbf{x})H_\epsilon (\phi(\mathbf{x})) $$ where $\mathbf{x} \in \mathbb{R}^2$, $\phi : \mathbb{R}^2 \rightarrow \mathbb{R}^2$, and $H_\epsilon$ is a relaxed version of the Heaviside step function.

I am interested in studying spectra of $\mathcal{H}$ as $\phi$ varies but am not sure about the definitions. Suppose $(v,\lambda)$ is an eigenpair of $\mathcal{H}$, I want to know $\nabla v$ and $\nabla \lambda$.

For example, I compute $\nabla \lambda$ as $$ \mathcal{H}v = \lambda v $$ $$ \implies \nabla (\mathcal{H}v) = \nabla (\lambda v) $$ $$ \implies \nabla \mathcal{H} v + \mathcal{H} \nabla v = \nabla \lambda v + \lambda \nabla v $$ $$ \implies \nabla \lambda = v^T \nabla \mathcal{H} v $$ $$ \implies \nabla \lambda = v^T V(\mathbf{x}) \delta_\epsilon(\phi(\mathbf{x})) v $$

which is weird since $\nabla \lambda$ is a number and I don't know that I computed $\nabla \mathcal{H}$ properly. It's the same way Terry does the computation in the single variable case (cf http://terrytao.wordpress.com/2008/10/28/when-are-eigenvalues-stable/ ).

Any advice? I hope to use $\phi$ in a level set method for domain decomposition.

Edit: Restrict everything to the discrete case (ie $\mathcal{H}$ is a matrix, like in the link above) and I still don't know how to do this. I guess the concise way to state my question is "Given a matrix $\mathcal{H}(\phi)$ does anyone know a formula for the first variation of the eigenvalues or eigenvectors of $\mathcal{H}$?"