A solution formula and the R-boundedness for the generalized Stokes resolvent problem in an infinite layer with Neumann boundary condition

We consider the generalized Stokes resolvent problem in an infinite layer with Neumann boundary conditions. This problem arises from a free boundary problem describing the motion of incompressible viscous one-phase fluid flow without surface tension in an infinite layer bounded both from above and from below by free surfaces. We derive a new exact solution formula to the generalized Stokes resolvent problem and prove the $\mathcal{R}$-boundedness of the solution operator families with resolvent parameter $\lambda$ varying in a sector ${\mathrm{\Sigma }}_{\epsilon ,{\gamma }_0}$ for any ${\gamma }_0>0$ and $0<\epsilon <\pi \mathrm{/}2$, where $\Sigma_{\varepsilon,\gamma_0} ={ \lambda\in$\mathbb{C}$\setminus{0} \mid |\arg\lambda|\leq\pi-\varepsilon,
|\lambda|>\gamma_0 }$\mathrm{.}Asapplications,weobtainthemaximal$L_p$-$L_q$ regularity for the nonstationary Stokes problem and then establish the well-posedness locally in time of the nonlinear free boundary problem mentioned above in $L_p$-$L_q$ setting. We make full use of the solution formula to take ${\gamma }_0>0$ arbitrarily, while in general domains we only know the $\mathcal{R}$-boundedness for ${\gamma }_0\gg 1$ from the result by Shibata. As compared with the case of Neumann-Dirichlet boundary condition studied by Saito, analysis is even harder on account of higher singularity of the symbols in the solution formula.