  The Bloch oscillations in the zigzag array of the optical waveguides are
considered. The multiple scattering formalism (MSF) is used for the numerical
simulation of the optical beam which propagates within the array. The effect of
the second-order coupling which depends on the geometrical parameters of the
array is investigated. The results obtained within the MSF are compared with
the calculation based on the phenomenological coupling modes model. The
calculations are performed for the waveguides fabricated in alkaline earth
boro-aluminosilicate glass sample, which are the most promising for the C-band
(1530-1565 nm).
