Nonlinear diffusion model for annealed proton-exchanged waveguides in zirconium-doped lithium niobate