It is known that in a box potential, when we set $V = 0$ inside and $V = \infty$ on the boundaries, the solution to the equation $$ - \frac{\hbar}{2m} \bigg( \frac{\partial^2}{\partial x^2} + \frac{ \partial^2}{\partial y^2} + \frac{\partial^2}{\partial z^2} \bigg) \psi(x,y,z) = E \psi(x,y,z) $$ is given by $$ \psi_{n_x, n_y, n_z} (x,y,z) = C \sin( \frac{n_x \pi x}{L}) \sin (\frac{n_y \pi y}{L}) \sin(\frac{n_z \pi z}{L}). $$ Here $L$ is the dimension of the box. The energies are $E_n = \frac{\hbar^2 \pi^2 }{2mL^2} n^2 $ with $n^2 = n_x^2 + n_y^2 + n_z^2$. Now, I was wondering, what if the potential inside the box is not zero, but we let it be $V = V_0$, a constant? What is the solution of the Schrödinger equation then? And what are the energy eigenvalues?
Does this amount to adding a phase factor to the solution (the product of sines)?