“…The 3D intensity distribution can be calculated using the approach proposed in ref. [24]. In the case of a squared regular grating, the light field diffracted behind the grating can be presented in the form Here, and , , k is the wavevector, w is the beam spot radius, is the reciprocal lattice vector, and the subscripts are integers.…”