Lorenz gauge fixing Is it always possible to define function $\psi$ satisfying the Lorenz gauge equation
$$
\partial_{\mu}\partial^{\mu} \psi + \partial_{\mu}A^{\mu} = 0? 
$$
 A: Yes. If you define $f=-\partial_\mu A^\mu$ then you can write the equation in the form $$ \partial_\mu\partial^\mu\psi = f$$
This is the Klein-Gordon equation with a nonzero source ($f$) and can be solved via Green's function methods. Once you have the Klein-Gordon propagator* $G(x)$ (this is derived in any e.g. quantum field theory textbook) appropriate to the boundary conditions the solution can be written as $$ \psi(x)=\int d^4 x' G(x-x') f(x')$$ since Green's functions by definition satisfy $$\partial_\mu\partial^\mu G(x-x')= \delta(x-x')$$ where we take all differentiations to be with respect to x.
*You need the propagator in the position space representation to write this down. It is usually more convenient to write it in momentum space; you can go back and forth using (inverse) Fourier transforms.
A: Yes, sure, it is always possible to find $\psi$ so that your equation will be obeyed – i.e. that the new $A_\mu$ will obey the Lorenz gauge – assuming that $A_\mu$ obeys the appropriate continuity conditions etc. 
There are many function like that (in flat infinite spacetime). For example, you may choose $\psi(x,y,z,t=0)$ arbitrarily and study the condition above at each point $(x,y,z)$ separately. Then the equation is just a very simple ordinary differential equation that depends on time which may be solved $dt$ after $dt$.
