'Half of the total angular spread.
Examples
>>> from sympy.physics.optics import BeamParameter
>>> p = BeamParameter(530e-9, 1, w=1e-3)
>>> p.divergence
0.00053/pi'

@property
def divergence(self):

    return ((self.wavelen / pi) / self.w_0)
