Evolution systems for paraxial wave equations of Schrödinger-type with non-smooth coefficients