Linearization of a PDE

I have been struggling with some linearization argument of the following paper: “M. Weinstein: Modulational stability of ground states of NLS”. In order to give a bit of context to my question, let us consider the NLS equation $ $ 2i\phi_t+\Delta\phi+\vert\phi\vert^{2\sigma}\phi=0, \quad 0<\sigma<\tfrac{2}{n-2}.$ $ This equation has very interesting “localized” solutions of the form: $ $ \phi(t,x)=u(x)e^{it/2},$ $ where $ u(t,x)$ solves $ \Delta u-u+\vert u\vert^{2\sigma}u=0$ . Besides, the latter equation has an even more interesting real, positive and radial $ H^1(\mathbb{R}^n)$ solution called “Ground state” and denoted by $ R(x)$ .

Now let me try to explain my question. Consider the perturbed Initial Valued Problem (IVP): $ $ 2i\phi_t^\varepsilon+\Delta \phi^\varepsilon+\vert \phi^\varepsilon\vert^{2\sigma}\phi^\varepsilon=\varepsilon F(\vert \phi^\varepsilon\vert)\phi^\varepsilon, \quad \phi^\varepsilon(t=0,x)=R(x)+\varepsilon S(x)$ $ We will seek solutions of the previous equation of the form $ $ \phi^\varepsilon(t,x)=(R(x)+\varepsilon w_1+\varepsilon^2 w_2+…)e^{it/2}.$ $ According to Weinstein if you reeplace this function into the perturbed equation and linearize you will get the following IVP for the linearized perturbation $ w$ : $ $ 2iw_t+\Delta w-w+(\sigma+1)R^{2\sigma}w+\sigma R^{2\sigma}\overline{w}=F(R)R, \quad w(0,x)=0.$ $ Now my problem is: I do not really understand how to obtain this linearization, can someone explain a little bit how to do it? Or recommend some references to learn about it. I tried replacing $ \phi^\varepsilon$ (truncated after $ \varepsilon w$ ) and then I took the derivative with respect to $ \varepsilon$ and evaluate at $ \varepsilon=0$ , but I cannot recover the equation claimed by Weinstein, so I think that I’m not undersitanding how to linearize.

Note2: The parameter $ n$ denotes the dimension.