Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Tailored diffractions of asymmetric columns and symmetric rows in two-dimensional multi-element phase gratings

Open Access Open Access

Abstract

Two-dimensional multi-element phase gratings can be engineered to show an even symmetry along one direction while an odd symmetry along the other direction in terms of offset refractive indices in each unit cell. The interplay of such even and odd symmetries has been explored to tailor diffraction columns and rows on demand by making offset refractive indices to satisfy specific requirements and hence attain different types of destructive interference. The resultant tailoring effects include the directional column elimination, the grouped column elimination, and the directional column selection as well as the natural row absence, the grouped row elimination, and the central row selection.

© 2023 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

Diffraction gratings, owing to engineered beam splitting functions, exhibit important applications in many frontier science and technology fields including physics, astronomy, chemistry, and biology [14]. The underlying physics of far-field light diffraction is known in particular as Fourier transform optics [5], which has been employed to accurately predict and control diverse diffraction behaviors indispensable in realistic applications. In other words, it provides a theoretical basis for designing etched phase gratings with neither loss nor gain, like holographic gratings, blazed gratings, and Dammann gratings [610], to achieve specific diffraction effects characterized by different intensity distributions along different directions. This is relevant to binary optics [1113] allowing for making micron-level relief phase elements based on computer-aided designs and has been extended to metasurfaces composed of subwavelength building blocks to gain a more flexible and complete control over the flow of light [1416].

Diffraction patterns of traditional phase gratings are generally symmetric, but asymmetric diffractions are becoming more and more essential in light-wave communications, micro-nano lithography, and other application scenarios though basically difficult to attain. This prompts us to recall the Friedel’s law [17] arguing that light diffraction from a real potential is expected to be symmetric when (a) scattering is weak and (b) there is no absorption in the corresponding optical crystals or structures [1820]. That indicates, we have two possible methods to realize asymmetric diffractions by breaking different necessary conditions of the Friedel’s law. The familiar one is to utilize parity-time (PT) symmetric or more general non-Hermitian gratings [2123] where absorption and/or gain is introduced to display an out-of-phase spatial modulation with respect to refractive index, hence yielding a complex optical potential with misaligned real and imaginary parts. The method of non-Hermitian modulation, when incorporated into electromagnetically induced gratings (EIGs) [2426], owns the advantage of coherent generation and dynamic manipulation, including both one-dimensional (1D) [27,24,2830] and two-dimensional (2D) [3133] examples. A main limitation of this method lies in that it is not viable to tailor specific diffraction orders on demand though diffraction asymmetry may be controlled on the whole, e.g. to quench diffracted beams in one side (1D case) or some quadrants (2D case).

To date, few works have been done to attain asymmetric diffractions, especially in the optical regime, by modulating the spatial distributions of only refractive indices in the presence of strong scattering [3436]. In a recent work, we have shown that 1D multi-element phase gratings could exhibit asymmetric diffractions as offset refractive indices of different elements in a unit cell are not in the spatially even symmetry and largest asymmetries occur always in the case of spatially odd-symmetric offset refractive indices [37]. It is more important that appropriate engineerings of odd-symmetric offset refractive indices allow us to attain a few interesting effects including the directional eliminations of a single diffraction order, all even or odd diffraction orders, and all but one or two diffraction orders. Then a question arises, what will happen and what are new if we consider 2D multi-element phase gratings?

In the spirit of extended binary optics, here we examine appropriate ways for engineering 2D multi-element phase gratings with offset refractive indices in each unit cell being odd symmetric along both diagonal directions $y=\pm x$, i.e. odd and even symmetric along the $y$ and $x$ directions, respectively. First, we show that it is viable to eliminate a few equidistant diffraction columns via paired destructive interference, eliminate all even or odd diffraction columns via dual paired destructive interference, and select a few equidistant diffraction columns via successive destructive interference as offset refractive indices satisfy specific requirements along the $y$ direction. Second, we show that these diffraction effects on tailored columns are accompanied by the natural absence of paired symmetric diffraction rows and it is viable to further attain the elimination of all even or odd diffraction rows and the definite selection of the central diffraction row by making offset refractive indices to satisfy some extra requirements along the $x$ direction. The joint tailorings of diffraction columns and rows allow for a flexible and versatile control over all diffracted beams and could even leave only two diffracted beams at desired output angles for instance.

2. Model and method

We start by considering a 2D phase grating of same period $a$ along both $x$ and $y$ directions, which will result in the far-field diffraction of a light beam incident along the $z$ direction into a few discrete orders marked by $i\in [-I,I]$ along the $x$ direction and $j\in [-J,J]$ along the $y$ direction as shown in Fig. 1(a). Each unit cell of this grating is designed to have $2m_{g}\times 2n_{g}$ elements with identical areas $\delta x\times \delta y=\frac {a}{2m_{g}}\times \frac {a}{2n_{g}}$ but different refractive indices $n_{m,n}=n_{0}+\delta n_{m,n}$ for $m\in \pm [1,m_{g}]$ and $n\in \pm [1,n_{g}]$, with $n_{0}$ being the mean value while $\delta n_{m,n}$ the offset values, as shown in Fig. 1(b). Taking $\mathcal {L}$ as the common thickness of all elements and considering a probe field with wavelength $\lambda _{p}$ and wavenumber $k_{p}=2\pi /\lambda _p$, for the $(m,n)$ element, we can write down the individual transmission function

$$T_{m,n}=e^{\imath n_{m,n}k_{p}\mathcal{L}}=e^{\imath\delta n_{m,n}k_{p}\mathcal{L}},$$
where the common and hence unimportant phase factor $e^{\imath n_{0}k_{p}\mathcal {L}}$ has been removed and $\imath$ stands for the imaginary unit. Then, each unit cell of the 2D multi-element phase grating can be described by the total transmission function
$$T(x,y)=\sum_{m={-}m_{g}}^{m_{g}}\sum_{n={-}n_{g}}^{n_{g}} T_{m,n}\cdot rect\left(\frac{x-x_{m}}{\delta x}\right)\cdot rect\left(\frac{y-y_{n}}{\delta y}\right),$$
with $x_{m}=(2m\mp 1)a/4m_{g}$ for $m\in \pm [1,m_{g}]$ while $y_{n}=(2n\mp 1)a/4n_{g}$ for $n\in \pm [1,n_{g}]$. That is, $(x_{m},y_{n})$ has been adopted to denote the center of the $(m,n)$ element. It is also worth noting that the terms referring to $m=0$ or $n=0$ are not present in the summation above as well as in the following one.

 figure: Fig. 1.

Fig. 1. $(a)$ Schematic of a phase grating of period $a$ along both $x$ and $y$ directions, which can diffract a light beam of amplitude $E_{in}$ incident along the $z$ direction into a few light beams of amplitudes $E_{i,j}$ deviating with angles $(\theta _{i},\theta _{j})$ from the $z$ direction. $(b)$ Fine structure in a period consisting of $2m_{g}\times 2n_{g}$ elements with a common thickness $\mathcal {L}$, identical widths $\delta x=a/2m_{g}$ and $\delta y=a/2n_{g}$, while different refractive indices $n_{m,n}=n_{0}+\delta n_{m,n}$ restricted by $\delta n_{m,n}=\delta n_{-m,n}$ and $\delta n_{m,n}=-\delta n_{m,-n}$.

Download Full Size | PDF

Further considering the translational invariance, we now make a Fourier transform to expand the total transmission function into $T(x,y)=\int _{-k_{p}}^{k_{p}}\int _{-k_{p}}^{k_{p}}E(\theta _{x},\theta _{y}){e^{\imath k_{x}x}}{e^{\imath k_{y}y}}dk_{x}dk_{y}$. Here $\theta _{x}$ ($\theta _{y}$) corresponds to an angle of diffracted beams deviating from the $z$ direction in the $xz$ ($yz$) plane, while $k_{x}=k_{p}\sin \theta _{x}$ ($k_{y}=k_{p}\sin \theta _{y}$) denotes a projection of probe wavenumber $k_{p}$ from the direction determined by angle $\theta _{x}$ ($\theta _{y}$) to the $x$ ($y$) direction. Then, an inverse procedure of this Fourier transform yields the following dimensionless field amplitude

$$E(\theta_{x},\theta_{y})=\frac{4\sin(k_{x}\delta x/2)\sin(k_{y}\delta y/2)}{a^2k_{x}k_{y}}\sum_{m={-}m_{g}}^{m_g}\sum_{n={-}n_g}^{n_g}e^{\imath(\delta n_{m,n}k_{p}\mathcal{L}-k_{x}x_{m}-k_{y}y_{n})},$$
for a normally incident beam diffracted into the direction defined by angles $\theta _{x}$ and $\theta _{y}$. It is thus straightforward to write down the diffraction intensity contributed by $M_{x}\times M_{y}$ irradiated periods in the $xy$ plane
$$I(\theta_{x},\theta_{y})=\vert E(\theta_{x},\theta_{y})\vert^2\left[\frac{\sin(\pi M_{x}R_x\sin\theta_{x})}{M_{x}\sin(\pi R_x\sin\theta_{x})}\frac{\sin(\pi M_{y}R_y\sin\theta_{y})}{M_{y}\sin(\pi R_y\sin\theta_{y})}\right]^2,$$
where beam radius $W_{b}$ has been taken into account in terms of ratios $M_{x}=M_{y}=W_{b}/a$. Diffraction peaks are known to occur at angles $\theta _{i}$ and $\theta _{j}$ determined by $k_{i}=k_{p}\sin \theta _{i}=i2\pi /a$ and $k_{j}=k_{p}\sin \theta _{j}=j2\pi /a$ for $i,j\in \{0,\pm 1,\pm 2,\ldots \}$, which can be further translated into $i=R_{x}\sin \theta _{i}$ and $j=R_{y}\sin \theta _{j}$ with $R_x=R_y=a/\lambda _p$. Based on this consideration, we would pay special attention to the $(i,j)$-order diffraction by examining $E_{i,j}\equiv E(\theta _{i},\theta _{j})$ in the following analytical discussions, though numerical calculations will be done still for $I(\theta _{x},\theta _{y})$ in Eq. (4).

In what follows, we restrict our discussions to the specific case where $\delta n_{\pm m,n}=-\delta n_{\pm m,-n}$ and $\delta n_{-m,\pm n}=\delta n_{m,\pm n}$ are simultaneously satisfied. That means, our 2D phase grating exhibits an even symmetry along the $x$ direction and meanwhile an odd symmetry along the $y$ direction, which then yields an odd symmetry along both diagonal directions defined by $y=\pm x$. In this case, Eq. (3) can be recast into a more compact form

$$E_{i,j}=\frac{16\sin(k_{i}\delta x/2)\sin(k_{j}\delta y/2)}{a^{2}k_{i}k_{j}}\sum_{m=1}^{m_g} \cos(k_{i}x_{m})\sum_{n=1}^{n_g}\cos(\delta n_{m,n}k_{p}\mathcal{L}-k_{j}y_{n}),$$
where $x_{m}=-x_{-m}$ and $y_{n}=-y_{-n}$ have been considered with $\cos (k_{i}x_{m})\cos (\delta n_{m,n}k_{p}\mathcal {L}-k_{j}y_{n})$ standing for a joint contribution of the four $(\pm m,\pm n)$ elements. It is thus clear that we can deal with the diffraction columns and rows along the $x$ and $y$ directions, respectively, to explore some intriguing effects though there exists an interplay between them bridged by the longitudinal traveling phase $\delta n_{m,n}k_{p}\mathcal {L}$. This will be done in the next two sections via analytical discussions and numerical calculations in the weak modulation limit of $|\delta n_{m,n}|\ll n_{0}$ so as to avoid multiple reflections between the phase grating and a surrounding medium of refractive index $n=n_{0}$.

3. Tailored asymmetric diffraction columns

In this section, we consider the case where the odd symmetry along the $y$ direction, namely $\delta n_{\pm m,n}=-\delta n_{\pm m,-n}$, is modulated while the even symmetry along the $x$ direction, namely $\delta n_{m,\pm n}=\delta n_{-m,\pm n}$, remains fixed.

Directional column elimination.-First, we note from Eq. (5) that the diffraction amplitude $E_{i,j}$ will become vanishing for a certain column order $j\in [-J,J]$, independent of the row order $i\in [-I,I]$, if we require $\delta n_{m,n}k_{p}\mathcal {L}=\delta n_{m,n}^{j}k_{p}\mathcal {L}\equiv k_{j}y_{n}+(l_{m,n}^j-1/2)\pi$, with $l_{m,n}^j$ being an integer. A simple rearrangement of this requirement yields

$$\delta n_{m,n}^j=\frac{[j(2n-1)+(2l_{m,n}^j-1)n_{g}]\lambda_{p}}{4n_{g}\mathcal{L}},$$
where $l_{m,n}^j$ can be chosen at will for $m\in [1,m_{g}]$ and $n\in [1,n_{g}]$ to realize the directional elimination of a diffraction column of order $j\in [-J,J]$. Sufficient conditions for realizing a directional column elimination are summarized in the left column of Table 1, according to which two examples on the integer matrices $\hat {l}^{1}_{m,n}$ and $\hat {l}^{-2}_{m,n}$ have been chosen to show that the $1$st and $-2$nd diffraction columns can be eliminated on demand in Fig. 2($a_1$) and Fig. 2($a_2$), respectively. Note, however, that the $-5$th and $4$th diffraction columns also disappear in Fig. 2($a_1$) and Fig. 2($a_2$), respectively, which can be easily understood by considering $k_{j\pm 2n_{g}}y_{n}-k_{j}y_{n}=\pm (2n-1)\pi$ in Eq. (5) so that $\cos (\delta n_{m,n}k_{p}\mathcal {L}-k_{j}y_{n})=0$ must result in $\cos (\delta n_{m,n}k_{p}\mathcal {L}-k_{j\pm 2n_{g}}y_{n})=0$. Moreover, we note from Fig. 2($a_1$) and Fig. 2($a_2$) that other surviving diffraction columns are clearly asymmetric as manifested by $E_{i,j'}\ne E_{i,-j'}$ with $j'\ne j$.

 figure: Fig. 2.

Fig. 2. Intensity $I$ against $\sin \theta _x$ and $\sin \theta _y$ with $E_{i,-5}=E_{i,1}=0$ $(a_1)$ while $E_{i,-2}=E_{i,4}=0$ ($a_2$) as we set integer matrices $\hat {l}^{1}_{m,n}=\hat {l}^{-2}_{m,n}=\{-1,0,2;-2,0,1;0,1,2\}$ to attain offset refractive indices $\delta n_{m,n}^{1}$ in ($b_1$) while $\delta n_{m,n}^{-2}$ in ($b_2$). Other parameters used in calculations are $\lambda _{p}=0.8$ $\mu$m, $m_{g}=3, n_{g}=3$, $R=30$, $M_x=M_y=10$, $\mathcal {L}=10\lambda _{p}$, and $n_{0}=2.0$.

Download Full Size | PDF

Tables Icon

Table 1. Conditions for three types of tailored asymmetric diffraction columns.

As in the 1D case [37], the directional column elimination arises from the destructive interference between paired elements in each unit cell while the asymmetric diffraction columns are resulted from a violation of the Friedel’s law. The paired destructive interference can be understood as follows: (i) different parts of a light beam scattered by the two $(m,\pm n)$ elements, respectively, suffer the opposite deflected phases $\pm \beta _{n}^{j}=\pm k_{j}y_{n}$ and the opposite forward phases $\pm \alpha _{m, n}^{j}=\pm \delta n_{m,n}^{j}k_{p}\mathcal {L}$ but share a common deflected phase $\gamma _{m}^{i}=k_{i}x_{m}$, which then results in a combined amplitude proportional to $e^{\imath \gamma _{m}^{i}}\cos {(\alpha _{m,n}^{j}-\beta _{n}^{j})}$; (ii) similarly, the two $(-m,\pm n)$ elements will contribute a combined amplitude proportional to $e^{-\imath \gamma _{m}^{i}}\cos {(\alpha _{m,n}^{j}-\beta _{n}^{j})}$ so that we have $\mathcal {E}^{i,j}_{m,n}\propto \cos {(\gamma _{m}^{i})}\cos {(\alpha _{m,n}^{j}-\beta _{n}^{j})}$ for the four $(\pm m,\pm n)$ elements as shown in Eq. (5); (iii) hence, $E_{i,j}$ for a certain column order $j\in [-J,J]$ will be vanishing when perfect destructive interference occurs between the $(m,n)$ and $(m,-n)$ elements as well as the $(-m,n)$ and $(-m,-n)$ elements in the case of $\alpha _{m,n}^{j}-\beta _{n}^{j}=(l_{m,n}^j-1/2)\pi$ as indicated by Eq. (6). On the other hand, we know that a violation of the Friedel’s law requires non-negligible higher-order scattering processes [37], which can be seen by making a power series expansion upon $\delta n_{m,n}k_{p}\mathcal {L}$ in Eq. (5) and explains why asymmetric diffraction columns are observed in Fig. 2.

Grouped column elimination.-Second, on the basis of $E_{i,j}=0$ for a certain value of integer $j$, we find from Eq. (5) and Eq. (6) that the diffraction column of order $j\pm 2k$ with $k\in \{1,2,3,\ldots \}$ exhibits the amplitude

$$E_{i,j\pm2k}=\frac{16\sin(k_{i}\delta x/2)\sin(k_{j\pm 2k}\delta y/2)}{a^2k_{i}k_{j\pm 2k}}\sum_{m=1}^{m_g}\cos(k_i x_m)\sum_{n=1}^{n_g}\sin{(l_{m,n}^{j}\pi\mp k\delta_{n})},$$
with $\delta _{n}=(2n-1)\pi /n_{g}$ being a $n$-dependent phase shift. Noting that $\delta _{n}+\delta _{n^{\ast }}=2\pi$ with $n^{\ast }=n_{g}+1-n$ for two conjugate columns, it is easy to find that $E_{i,j\pm 2k}$ will also become vanishing if we further require
$$\frac{\sin{(l_{m,n}^{j}\pi\mp k\delta_{n})}+\sin{(l_{m,n^{{\ast}}}^{j}\pi\mp k\delta_{n^{{\ast}}})}}{2}=\sin{\left(\frac{l_{m,n}^{j}+l_{m,n^{{\ast}}}^{j}}{2}\pi\right)}\cos{\left(\frac{l_{m,n}^{j}-l_{m,n^{{\ast}}}^{j}}{2}\pi \mp \frac{2n-1}{n_{g}}k\pi\right)}=0,$$
no matter $n_{g}$ is an even or odd integer. It is thus clear that all even or odd diffraction columns will disappear provided (i) $j$ is an even or odd integer and (ii) $l_{m,n}^{j}+l_{m,n^{\ast }}^{j}$ is chosen as an even integer, independent of the value of integer $k$. Sufficient conditions for realizing a grouped column elimination are summarized in the middle column of Table 1, according to which two examples on the integer matrices $\hat {l}^{1}_{m,n}$ and $\hat {l}^{-2}_{m,n}$ have been chosen to show that the odd and even diffraction columns can be eliminated on demand in Fig. 3($a_1$) and Fig. 3($a_2$), respectively. That is, on the basis of $E_{i,1}=0$ in Fig. 3($a_1$) and $E_{i,-2}=0$ in Fig. 3($a_2$), other odd and even diffraction columns have been eliminated by further requiring $l_{m,n}^{1}+l_{m,n^{\ast }}^{1}=2k$ and $l_{m,n}^{-2}+l_{m,n^{\ast }}^{-2}=2k$, respectively, with $k\in \{1,2,3,\ldots \}$.

 figure: Fig. 3.

Fig. 3. Intensity $I$ against $\sin \theta _x$ and $\sin \theta _y$ with all odd columns eliminated on the basis of $E_{i,1}=0$ $(a_1)$ while all even columns eliminated on the basis of $E_{i,-2}=0$ $(a_2)$. Parameters are the same as in Fig. 2 except $\hat {l}^{1}_{m,n}=\hat {l}^{-2}_{m,n}=\{-1,0,1;-2,0,0;0,1,2\}$ restricted by $l^{1,-2}_{m,n}+l^{1,-2}_{m,n^\ast }\in \{0,\pm 2,\pm 4,\ldots \}$ are reset to attain $\delta n_{m,n}^{1}$ in ($b_1$) while $\delta n_{m,n}^{-2}$ in ($b_2$).

Download Full Size | PDF

The grouped column elimination arises instead from the destructive interference between dual paired elements and can be understood as follows: (i) the $j$th diffraction column has been eliminated with $E_{i,j}=\sum _{m=1}^{m_{g}}\sum _{n=1}^{n_{g}}\mathcal {E}^{i,j}_{m,n}=0$ by taking $\alpha ^{j}_{m,n}-\beta ^{j}_{n}=(l_{m,n}^{j}-1/2)\pi$ in $\mathcal {E}^{i,j}_{m,n}\propto \cos {(\alpha _{m,n}^{j}-\beta _{n}^{j})}$ for the four $(\pm m,\pm n)$ elements; (ii) as to the diffraction column of order $j\pm 2k$, we have $\mathcal {E}^{i,j\pm 2k}_{m,n}\propto \sin (l_{m,n}^{j}\pi \mp k\delta _{n})$ for the $(\pm m,\pm n)$ elements and $\mathcal {E}^{i,j\pm 2k}_{m,n^{\ast }}\propto \sin (l_{m,n^{\ast }}^{j}\pi \mp k\delta _{n^{\ast }})$ for the $(\pm m,\pm n^{\ast })$ elements by considering $\alpha ^{j\pm 2k}_{m,n}=\alpha ^{j}_{m,n}$ and $\beta ^{j\pm 2k}_{m,n}=\beta ^{j}_{n}\pm k\delta _{n}$; (iii) taking $\delta _{n}+\delta _{n^{\ast }}=2\pi$ into account, it is not difficult to find that the $(\pm m,\pm n)$ and $(\pm m,\pm n^{\ast })$ elements jointly contribute a beam superposition of amplitude $\mathcal {E}^{i,j\pm 2k}_{m,n}+\mathcal {E}^{i,j\pm 2k}_{m,n^{\ast }}\propto \sin [(l_{m,n}^{j}+l_{m,n^{\ast }}^{j})\pi /2]$; (iv) hence, the diffraction column of order $j\pm 2k$ will disappear with $E_{i,j\pm 2k}=0$ due to perfect destructive interference between the ‘dual’ paired $(m,\pm n)$ and $(m,\pm n^{\ast })$ elements as well as $(-m,\pm n)$ and $(-m,\pm n^{\ast })$ elements provided $l^j_{m,n}+l^j_{m,n^{\ast }}$ is an even integer.

Directional column selection.-Finally, on the basis of $E_{i,j}=0$ for a certain value of integer $j$, we find from Eq. (5) and Eq. (6) that the diffraction column of order $j\pm k$ with $k\in \{1,2,3,\ldots \}$ exhibits the amplitude

$$E_{i,j\pm k}=\frac{4\sin(i\pi/2m_g)\sin[(j\pm k)\pi/2n_{g}]}{i(j\pm k)\pi^2}\sum_{m=1}^{m_g}\cos(k_i x_m)\sum_{n=1}^{n_g}\sin{\left(l_{m,n}^{j}\pi\mp \frac{k\delta_{n}}{2}\right)}.$$

This diffraction amplitude, upon the multiplication of a non-zero term $\cos (k\delta _{1}/2)$ in the case of $k\ne \widetilde {n}_{g}\in \{1,3,5,\ldots \}n_{g}$, will become vanishing if we further require

$$\sum_{n=1}^{n_g}\cos\left(\frac{k\delta_{1}}{2}\right)\sin{\left(l_{m,n}^{j}\pi\mp \frac{k\delta_{n}}{2}\right)} =\frac{1}{2}\sum_{n=1}^{n_g}\left[\sin{\left(l_{m,n}^{j}\pi\mp \frac{k(n-1)}{n_{g}}\pi\right)}+\sin{\left(l_{m,n}^{j}\pi\mp \frac{kn}{n_{g}}\pi\right)}\right]=0.$$

By defining $l_{m,n_{g}+1}^{j}\equiv l_{m,1}^{j}\pm k$, we can translate this requirement into

$$\sum_{n=1}^{n_g}\sin{\left(\frac{l_{m,n}^{j}+l_{m,n+1}^{j}}{2}\pi\mp \frac{kn}{n_{g}}\pi\right)}\cos{\left(\frac{l_{m,n+1}^{j}-l_{m,n}^{j}}{2}\pi\right)}=0,$$
after a sequential summation of the $2n_{g}$ sine functions. It is thus clear that we can attain $E_{i,j\pm k}=0$ for $k\ne \widetilde {n}_{g}$ if $l_{m,n+1}^{j}-l_{m,n}^{j}$ is chosen as an odd integer to ensure $\cos {[(l_{m,n+1}^{j}-l_{m,n}^{j})\pi /2]}=0$. Sufficient conditions for realizing a directional column selection are summarized in the right column of Table 1, according to which two examples on the integer matrices $\hat {l}^{1}_{m,n}$ and $\hat {l}^{-2}_{m,n}$ have been chosen to show that only the diffraction columns of order $j=1\pm \widetilde {n}_{g}$ and $j=-2\pm \widetilde {n}_{g}$ survive in Fig. 4($a_1$) and Fig. 4($a_2$), respectively. That is, the $-2$nd and $4$th diffraction columns survive in Fig. 4($a_1$) on the basis of $E_{i,1}=0$ by requiring $l_{m,n+1}^{1}-l_{m,n}^{1}=2k-1$ while the $-5$th and $1$st diffraction columns survive in Fig. 4($a_2$) on the basis of $E_{i,-2}=0$ by requiring $l_{m,n+1}^{-2}-l_{m,n}^{-2}=2k-1$.

 figure: Fig. 4.

Fig. 4. Intensity $I$ against $\sin \theta _x$ and $\sin \theta _y$ with two even columns selected on the basis of $E_{i,1}=0$ $(a_1)$ while two odd columns selected on the basis of $E_{i,-2}=0$ $(a_2)$. Parameters are the same as in Fig. 2 except $\hat {l}^{1}_{m,n}=\hat {l}^{-2}_{m,n}=\{-1,0,1;-2,-1,0;0,1,2\}$ restricted by $l^{1,-2}_{m,n+1}-l^{1,-2}_{m,n}\in \{\pm 1,\pm 3,\ldots \}$ are reset to attain $\delta n_{m,n}^{1}$ in ($b_1$) while $\delta n_{m,n}^{-2}$ in ($b_2$).

Download Full Size | PDF

The directional column selection arises however from the destructive interference between successive paired elements and can be understood as follows: (i) the $j$th diffraction column has been eliminated with $E_{i,j}=\sum _{m=1}^{m_{g}}\sum _{n=1}^{n_{g}}\mathcal {E}^{i,j}_{m,n}=0$ by taking $\alpha ^{j}_{m,n}-\beta ^{j}_{n}=(l_{m,n}^{j}-1/2)\pi$ in $\mathcal {E}^{i,j}_{m,n}\propto \cos {(\alpha _{m,n}^{j}-\beta _{n}^{j})}$ for the four $(\pm m,\pm n)$ elements; (ii) as to the diffraction column of order $j\pm k$, we have $\mathcal {E}^{i,j\pm k}_{m,n}\propto \sin (l_{m,n}^{j}\pi \mp k\delta _{n}/2)$ for the $(\pm m,\pm n)$ elements, which can be rewritten as a sum of inward $\mathcal {E}^{i,j\pm k,in}_{m,n}\propto \sin [l_{m,n}^{j}\pi \mp k(n-1)\pi /n_{g}]$ and outward $\mathcal {E}^{i,j\pm k,out}_{m,n}\propto \sin (l_{m,n}^{j}\pi \mp kn\pi /n_{g})$ amplitudes after multiplied by a non-zero $\cos (k\delta _{1}/2)$ for $k\ne \widetilde {n}_{g}$; (iii) it is not difficult to find that the sum of $\mathcal {E}^{i,j\pm k,out}_{m,n}$ contributed by the $(\pm m,\pm n)$ elements and $\mathcal {E}^{i,j\pm k,in}_{m,n+1}$ contributed by the $(\pm m,\pm n\pm 1)$ elements is proportional to $\cos [(l_{m,n+1}^{j}-l_{m,n}^{j})\pi /2]$; (iv) hence, all diffraction columns except those of orders $j\pm \widetilde {n}_{g}$ will disappear due to destructive interference between the successive $(\pm m,\pm n)$ and $(\pm m,\pm n\pm 1)$ elements provided $l^j_{m,n+1}-l^j_{m,n}$ is an odd integer.

4. Tailored symmetric diffraction rows

In this section, we examine whether it is viable to further tailor the diffraction rows while reserving the phenomena of directional column elimination, grouped column elimination, and directional column selection. This will be done by exerting extra requirements on integers $l_{m,n}^{j}$ with respect to subscript $m$ in the case of $\delta n_{\pm m,n}=-\delta n_{\pm m,-n}$ (odd symmetry along the $y$ direction) and $\delta n_{m,\pm n}=\delta n_{-m,\pm n}$ (even symmetry along the $x$ direction).

Natural row absence.-First, we note that Eq. (9) turns into

$$E_{q m_g,j\pm k}=\frac{4\sin[(j\pm k)\pi/2n_g]}{qm_{g}(j\pm k)\pi^2} \sum_{m=1}^{m_g}({-}1)^{qm}\sin(q\pi/2)\cos(q\pi/2)\sum_{n=1}^{n_g}\sin{\left(l_{m,n}^{j}\pi\mp \frac{k\delta_{n}}{2}\right)},$$
for the diffraction rows of orders $q{m}_g$ with $q\in \{\pm 1,\pm 2,\ldots \}$. It is clear that these diffraction rows naturally disappear due to $\sin (q\pi /2)\cos (q\pi /2)=0$ without extra requirements beyond $l^{j}_{m,n}$ being an odd or even integer for the directional column elimination, $l^{j}_{m,n}+l^{j}_{m,n^{\ast }}$ being an even integer for the grouped column elimination, and $l^{j}_{m,n+1}-l^{j}_{m,n}$ being an odd integer for the directional column selection. The absence of $\pm 3$th diffraction rows is intrinsic to the odd and even symmetries of $\delta n_{m,n}$ along the $y$ and $x$ directions, respectively, and has been marked by the yellow dashed lines in Fig. 2 - Fig. 4. This natural row absence arises from the destructive interference between paired elements in each column of a unit cell and can be understood as follows: (i) in $E_{qm_{g},j\pm k}=\sum _{m=1}^{m_{g}}\sum _{n=1}^{n_{g}}\mathcal {E}^{qm_{g},j\pm k}_{m,n}$ with $j\pm k$ covering all diffraction columns, we have $\mathcal {E}^{qm_{g},j\pm k}_{m,n}\propto \sin (q\pi /2)\cos (q\pi /2)\sin (l_{m,n}^{j}\pi \mp k\delta _{n}/2)$ as a contribution of the $(\pm m,\pm n)$ elements; (ii) hence, the diffraction rows of orders $qm_{g}$ disappear due to destructive interference between the $(m,n)$ and $(-m,n)$ elements as well as the $(m,-n)$ and $(-m,-n)$ elements, independent of $j\pm k$.

Grouped row elimination.-Second, on the basis of $E_{i,j}=0$ for a certain value of integer $j$, we rewrite Eq. (9) into $E_{i,j\pm k}=\sum _{m=1}^{m_{g}}\sum _{n=1}^{n_{g}}\mathcal {E}^{i,j\pm k}_{m,n}$ where $\mathcal {E}^{i,j\pm k}_{m,n}\propto \cos (k_i x_m)\sin (l_{m,n}^{j}\pi \mp k\delta _{n}/2)$ is a contribution from the four $(\pm m,\pm n)$ elements. Similarly, we have $\mathcal {E}^{i,j\pm k}_{m^{\ast },n}\propto \cos (k_i x_{m^{\ast }})\sin (l_{m^{\ast },n}^{j}\pi \mp k\delta _{n}/2)$ contributed by other four $(\pm m^{\ast },\pm n)$ conjugate elements with $m^{\ast }=m_{g}+1-m$. Then, it is easy to prove $\cos (k_i x_{m^{\ast }})=\cos (i\pi )\cos (k_i x_{m})$ and attain

$$\mathcal{E}^{i,j\pm k}_{m,n}+\mathcal{E}^{i,j\pm k}_{m^{{\ast}},n}\propto2\cos(k_i x_m) \sin\left[\frac{(l_{m,n}^{j}+l_{m^\ast,n}^{j})\pi}{2}\mp \frac{k \delta_n}{2}\right]\cos\left[\frac{(l_{m,n}^{j}-l_{m^\ast,n}^{j})\pi}{2}\right],$$
if $i$ is an even integer while
$$\mathcal{E}^{i,j\pm k}_{m,n}+\mathcal{E}^{i,j\pm k}_{m^{{\ast}},n}\propto2\cos(k_i x_m) \sin\left[\frac{(l_{m,n}^{j}+l_{m^\ast,n}^{j})\pi}{2}\mp \frac{k \delta_n}{2}\right]\sin\left[\frac{(l_{m,n}^{j}-l_{m^\ast,n}^{j})\pi}{2}\right],$$
if $i$ is an odd integer. That means, all even (odd) rows will disappear if we further require $l_{m,n}^j-l_{m^\ast,n}^{j}$ to be an odd (even) integer, in addition to $l_{m,n}^j$ being an odd or even integer for the directional column elimination, $l_{m,n}^j+l_{m,n^\ast }^{j}$ being an even integer for the grouped column elimination, and $l_{m,n+1}^j-l_{m,n}^{j}$ being an odd integer for the directional column selection, which has been listed in the left (right) column in Table 2. An example on the grouped elimination of all odd rows on the basis of $E_{i,2}=0$ is given in Fig. 5, while an example on the grouped elimination of all even rows on the basis of $E_{i,1}=0$ is given in Fig. 6, where directional column elimination, grouped column elimination, and directional column selection can be found, respectively, in panels ($a_{1}$), ($a_{2}$), and ($a_{3}$).

 figure: Fig. 5.

Fig. 5. Intensity $I$ against $\sin \theta _x$ and $\sin \theta _y$ with all odd rows eliminated on the basis of $E_{i,2}=0$ in the presence of directional column elimination $(a_1)$; grouped column elimination $(a_2)$; directional column selection $(a_3)$. Parameters are the same as in Fig. 2 except corresponding integer matrices restricted by $l^{2}_{m,n}-l^{2}_{m^\ast,n}\in \{0,\pm 2,\pm 4,\ldots \}$ are shown in $(b_1)$, $(b_2)$, and $(b_3)$ while corresponding offset refractive indices are shown in $(c_1)$, $(c_2)$, and $(c_3)$.

Download Full Size | PDF

 figure: Fig. 6.

Fig. 6. Intensity $I$ against $\sin \theta _x$ and $\sin \theta _y$ with all even rows eliminated on the basis of $E_{i,1}=0$ in the presence of directional column elimination $(a_1)$; grouped column elimination $(a_2)$; directional column selection $(a_3)$. Parameters are the same as in Fig. 2 except $m_g=n_g=4$ and corresponding integer matrices restricted by $l^{1}_{m,n}-l^{1}_{m^\ast,n}\in \{\pm 1,\pm 3,\ldots \}$ are shown in $(b_1)$, $(b_2)$, and $(b_3)$ while corresponding offset refractive indices are shown in $(c_1)$, $(c_2)$, and $(c_3)$.

Download Full Size | PDF

Tables Icon

Table 2. Conditions for three types of tailored symmetric diffraction rows.

The grouped row elimination arises from the destructive interference between dual paired elements in each column of a unit cell and can be understood as follows: (i) in $E_{i,j\pm k}=\sum _{m=1}^{m_{g}}\sum _{n=1}^{n_{g}}\mathcal {E}^{i,j\pm k}_{m,n}$ with $j\pm k$ covering all diffraction columns, we have $\mathcal {E}^{i,j\pm k}_{m,n}\propto \cos (k_{i}x_{m})\sin (l_{m,n}^{j}\pi \mp k\delta _{n}/2)$ and $\mathcal {E}^{i,j\pm k}_{m^\ast,n}\propto \cos (k_{i}x_{m^\ast })\sin (l_{m,n}^{j}\pi \mp k\delta _{n}/2)$ contributed by the $(\pm m,\pm n)$ and $(\pm m^\ast,\pm n)$ elements, respectively; (ii) taking $\cos (k_i x_{m^{\ast }})=\cos (i\pi )\cos (k_i x_{m})$ into account, it is not difficult to find that the $(\pm m,\pm n)$ and $(\pm m^\ast,\pm n)$ elements jointly contribute a beam superposition of amplitude $\mathcal {E}^{i,j\pm k}_{m,n}+\mathcal {E}^{i,j\pm k}_{m^{\ast },n}\propto \cos [(l_{m,n}^{j}-l_{m^\ast,n}^{j})\pi /2]$ for $i\in \{0,\pm 2,\pm 4,\ldots \}$ while $\mathcal {E}^{i,j\pm k}_{m,n}+\mathcal {E}^{i,j\pm k}_{m^{\ast },n}\propto \sin [(l_{m,n}^{j}-l_{m^\ast,n}^{j})\pi /2]$ for $i\in \{\pm 1,\pm 3,\ldots \}$; (iii) hence, all odd (even) diffraction rows will disappear due to perfect destructive interference between the dual paired conjugate $(\pm m, n)$ and $(\pm m^\ast,n)$ elements as well as the dual paired conjugate $(\pm m, -n)$ and $(\pm m^\ast, -n)$ elements provided $l^j_{m,n}-l^j_{m^{\ast },n}$ is an even (odd) integer, independent of $j\pm k$.

Central row selection.-Finally, we note that Eq. (9) will turn into

$$E_{i,j\pm k}=\frac{4\sin(i\pi/2m_g)\sin[(j\pm k)\pi/2n_{g}]}{i(j\pm k)\pi^2}\sum_{m=1}^{m_g}\cos(k_i x_m)\sum_{n=1}^{n_g}\sin{\left(l_{m_{g},n}^{j}\pi\mp \frac{k\delta_{n}}{2}\right)},$$
due to $\sum _{n=1}^{n_g}\sin {(l_{1,n}^{j}\pi \mp k\delta _{n}/2)}=\sum _{n=1}^{n_g}\sin {(l_{2,n}^{j}\pi \mp k\delta _{n}/2)}=\cdots =\sum _{n=1}^{n_g}\sin {(l_{m_{g},n}^{j}\pi \mp k\delta _{n}/2)}$ if we require $l_{m+1,n}^j-l_{m,n}^j$ to be an even integer. That means, there exists no interplay again between the two summations in $E_{i,j\pm k}$ over $m$ and $n$, respectively, so that they can be dealt with independently. On this basis, it is not difficult to prove that we have both $\sum _{m=1}^{m_g}\cos (k_i x_m)\ne 0$ and $\sin (i\pi /2m_{g})\ne 0$ only in the case of $i=0$, hence all diffraction rows except the central one will disappear in the case of $l_{m+1,n}^j-l_{m,n}^j\in \{0,\pm 2,\pm 4,\ldots \}$, which has been listed in the middle column of Table 2. An example on this central row selection is given in Fig. 7 where directional column elimination, grouped column elimination, and directional column selection can be found, respectively, in panels ($a_{1}$), ($a_{2}$), and ($a_{3}$), leading to a few or just two asymmetric diffracted beams at desired output angles along the $y$ direction.

 figure: Fig. 7.

Fig. 7. Intensity $I$ against $\sin \theta _x$ and $\sin \theta _y$ with the central row selected on the basis of $E_{i,1}=0$ in the presence of directional column elimination $(a_1)$; grouped column elimination $(a_2)$; directional column selection $(a_3)$. Parameters are the same as in Fig. 2 except $m_g=n_g=4$ and corresponding integer matrices restricted by $l^{1}_{m+1,n}-l^{1}_{m,n}\in \{0,\pm 2,\pm 4,\ldots \}$ are shown in $(b_1)$, $(b_2)$, and $(b_3)$ while corresponding offset refractive indices are shown in $(c_1)$, $(c_2)$, and $(c_3)$.

Download Full Size | PDF

The central row selection arises from the collective destructive interference between all $m_{g}$ elements in each column as the total contribution of all $n_{g}$ elements is identical for each row. This can be understood as follows: (i) we note from Eq. (15) that the summation $\sum _{m=1}^{m_g}\cos (k_i x_m)$ is independent of $n$ and the summation $\sum _{n=1}^{n_g}\sin {(l_{m_{g},n}^{j}\pi \mp k\delta _{n}/2)}$ is independent of $m$ in the case of $l^{j}_{m+1,n}-l^{j}_{m,n}\in \{0,\pm 2,\pm 4,\ldots \}$; (ii) the separable summation $\sum _{m=1}^{m_g}\cos (k_i x_m)$ is definitely zero for $i\ne 0$ because the angles $\{k_{i}x_{1},k_{i}x_{2},\ldots,k_{i}x_{m_{g}}\}$ are symmetrically distributed on a circle so that their cosine functions must cancel each other out; (iii) hence, all diffraction rows except the central one will disappear due to collective destructive interference as a joint contribution of all $m_g$ elements in each column.

5. Conclusions

In summary, a 2D multi-element phase grating has been designed to exhibit diverse odd-symmetric structures along two diagonal $y=\pm x$ directions in terms of offset refractive indices in each unit cell. It is of interest that this grating with neither loss nor gain is characterized by asymmetric diffraction columns along the $y$ direction while symmetric diffraction rows along the $x$ direction. By engineering offset refractive indices along the $y$ direction to satisfy specific requirements $\delta _{m,n}k_p\mathcal {L}-k_jy_n= (l_{m,n}^{j} -1/2)\pi$, we can realize on demand the elimination of a few equidistant diffraction columns, the elimination of all even or odd diffraction columns, or the selection of a few equidistant diffraction columns, all accompanied by the natural absence of symmetric row pairs. Similar engineerings of offset refractive indices along the $x$ direction allow us to further eliminate all even or odd diffraction rows, or definitely select the central diffraction row, on the basis of tailored diffraction columns mentioned above. These diffraction effects of tailored columns and rows can be explained in terms of paired, dual paired, successive, or collective destructive interference between different elements in each column or row of a unit cell, and may find applications in devising photonic devices requiring peculiar light filtering functions like eliminating or selecting one or more specific diffracted beams, especially when extended to some materials with tunable refractive indices [38,39]. Our results can also be generalized to realize diagonal row selection with other 2D multi-element phase gratings exhibiting, e.g., odd-symmetric offset refractive indices along a certain diagonal direction. A combination of our analytical method and numerical approaches based on iterative Fourier-transform algorithms (IFTA) [40] will be more efficient in realizing various beam filtering and shaping functions.

Funding

National Natural Science Foundation of China (12074061); National Key Research and Development Program of China (2021YFE0193500).

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

1. G. Quaranta, G. Basset, O. Martin, and B. Gallinet, “Recent advances in resonant waveguide gratings,” Laser Photonics Rev. 12(9), 1800017 (2018). [CrossRef]  

2. N. Bonod and J. Neauport, “Diffraction gratings: from principles to applications in high-intensity lasers,” Adv. Opt. Photonics 8(1), 156–199 (2016). [CrossRef]  

3. F. Esposito, L. Sansone, A. Srivastava, F. Baldini, S. Campopiano, F. Chiavaioli, M. Giordano, A. Giannetti, and A. Iadicicco, “Long period grating in double cladding fiber coated with graphene oxide as high-performance optical platform for biosensing,” Biosens. Bioelectron. 172, 112747 (2021). [CrossRef]  

4. C. Francesco, G. Carlos, J. Pedro, and B. Francesco, “Towards a uniform metrological assessment of grating-based optical fiber sensors: From refractometers to biosensors,” Biosensors 7(4), 23 (2017). [CrossRef]  

5. Z. Yang, T. Albrow-Owen, W. Cai, and T. Hasan, “Miniaturization of optical spectrometers,” Science 371(6528), eabe0722 (2021). [CrossRef]  

6. Y. L. Li, N. N. Li, D. Wang, F. Chu, S. D. Lee, Y. W. Zheng, and Q. H. Wang, “Tunable liquid crystal grating based holographic 3d display system with wide viewing angle and large size,” Light: Sci. Appl. 11(1), 188–1584 (2022). [CrossRef]  

7. L. Zhao, H. Tian, H. Liu, W. Zhang, F. Zhao, X. Song, and J. Shao, “Bio-inspired soft-rigid hybrid smart artificial muscle based on liquid crystal elastomer and helical metal wire,” Small 19(17), 2370113 (2023). [CrossRef]  

8. C. Ziqi, Y. Qian, S. Keita, L. Peng, H. Yupeng, H. Yao, W. Xibin, Z. Tianfeng, and M. Masayoshi, “High-precision and high-efficiency fabrication of blazed grating by ultrasonic-assisted ultraprecision planing,” J. Mater. Process. Technol. 311, 117802 (2023). [CrossRef]  

9. J. Yu, C. Zhou, W. Jia, J. Ma, and A. Hu, “Distorted dammann grating,” Opt. Lett. 38(4), 474–476 (2013). [CrossRef]  

10. R. V. Vinu, M. K. Sharma, R. K. Singh, and P. Senthilkumaran, “Generation of spatial coherence comb using dammann grating,” Opt. Lett. 39(8), 2407–2410 (2014). [CrossRef]  

11. M. P. Cagigal, A. Fuentes, P. J. Valle, V. F. Canales, and M. A. Cagigas, “Optical-component-only adaptive optics,” Opt. Lett. 46(14), 3452–3455 (2021). [CrossRef]  

12. K. Sun, H. Jiang, D. A. Bykov, V. Van, U. Levy, Y. Cai, and Z. Han, “1d quasi-bound states in the continuum with large operation bandwidth in the omega similar to k space for nonlinear optical applications,” Photonics Res. 10(7), 1575–1581 (2022). [CrossRef]  

13. S. Marchesini and A. Sakdinawat, “Shaping coherent x-rays with binary optics,” Opt. Express 27(2), 907–917 (2019). [CrossRef]  

14. M. Khorasaninejad and F. Capasso, “Metalenses: Versatile multifunctional photonic components,” Science 358(6367), aam8100 (2017). [CrossRef]  

15. A. Kristensen, J. K. W. Yang, S. I. Bozhevolnyi, S. Link, P. Nordlander, N. J. Halas, and N. A. Mortensen, “Plasmonic colour generation,” Nat. Rev. Mater. 2(1), 16088 (2017). [CrossRef]  

16. D. Lin, P. Fan, E. Hasman, and M. L. Brongersma, “Dielectric gradient metasurface optical elements,” Science 345(6194), 298–302 (2014). [CrossRef]  

17. G. Friedel, “Sur les symétries cristallines que peut révéler la diffraction des rayons röntgen,” C. R. Acad. Sci. (Paris) 157, 1533–1536 (1913).

18. S. Miyake and R. Uyeda, “Friedel’s law in the dynamical theory,” Acta Crystallogr. 8(6), 335–342 (1955). [CrossRef]  

19. A. Ambuj, E. Walla, S. Andaloro, S. Nomoto, R. Vyas, and S. Singh, “Symmetry in the diffraction of beams carrying orbital angular momentum,” Phys. Rev. A 99(1), 013846 (2019). [CrossRef]  

20. R. F. Berger, S. Lee, J. Johnson, B. Nebgen, F. Sha, and J. Xu, “The mystery of perpendicular fivefold axes and the fourth dimension in intermetallic structures,” Chem. Eur. J. 14(13), 3908–3930 (2008). [CrossRef]  

21. C. Chen, Y. Liu, L. Zhao, X. Hu, and Y. Fu, “Asymmetric nonlinear-mode-conversion in an optical waveguide with pt symmetry,” Front. Phys. 17(5), 52504 (2022). [CrossRef]  

22. J. T. Cole, K. G. Makris, Z. H. Musslimani, D. N. Christodoulides, and S. Rotter, “Twofold pt symmetry in doubly exponential optical lattices,” Phys. Rev. A 93(1), 013803 (2016). [CrossRef]  

23. M. V. Berry, “Lop-sided diffraction by absorbing crystals,” J. Phys. A: Math. Gen. 31(15), 3493–3502 (1998). [CrossRef]  

24. T. Shui, W. X. Yang, S. Liu, L. Li, and Z. Zhu, “Asymmetric diffraction by atomic gratings with optical pt symmetry in the raman-nath regime,” Phys. Rev. A 97(3), 033819 (2018). [CrossRef]  

25. D. Ma, D. Yu, X. D. Zhao, and J. Qian, “Unidirectional and controllable higher-order diffraction by a rydberg electromagnetically induced grating,” Phys. Rev. A 99(3), 033826 (2019). [CrossRef]  

26. C. Hang, W. Li, and G. Huang, “Nonlinear light diffraction by electromagnetically induced gratings with pt symmetry in a rydberg atomic gas,” Phys. Rev. A 100(4), 043807 (2019). [CrossRef]  

27. J. Gao, C. Hang, and G. Huang, “Linear and nonlinear bragg diffraction by electromagnetically induced gratings with pt symmetry and their active control in a rydberg atomic gas,” Phys. Rev. A 105(6), 063511 (2022). [CrossRef]  

28. X. Y. Zhu, Y. L. Xu, Y. Zou, X. C. Sun, C. He, M.-H. Lu, X.-P. Liu, and Y.-F. Chen, “Asymmetric diffraction based on a passive parity-time grating,” Appl. Phys. Lett. 109(11), 111101 (2016). [CrossRef]  

29. F. Gao, Y.-M. Liu, X.-D. Tian, C.-L. Cui, and J.-H. Wu, “Intrinsic link of asymmetric reflection and diffraction in non-hermitian gratings,” Opt. Express 26(26), 33818–33829 (2018). [CrossRef]  

30. Y. M. Liu, F. Gao, C. H. Fan, and J. H. Wu, “Asymmetric light diffraction of an atomic grating with pt symmetry,” Opt. Lett. 42(21), 4283–4286 (2017). [CrossRef]  

31. Y. M. Liu, F. Gao, J. H. Wu, M. Artoni, and G. C. La Rocca, “Lopsided diffractions of distinct symmetries in two-dimensional non-hermitian optical gratings,” Phys. Rev. A 100(4), 043801 (2019). [CrossRef]  

32. A. Hafiz and S. Qamar, “Electromagnetically induced holographic imaging using monolayer graphene,” Opt. Express 28(2), 1970–1986 (2020). [CrossRef]  

33. T. Shui, W. X. Yang, L. Li, and X. Wang, “Lop-sided raman-nath diffraction in pt-antisymmetric atomic lattices,” Opt. Lett. 44(8), 2089–2092 (2019). [CrossRef]  

34. Y. S. Hsieh, Y. C. Ho, S. Y. Lee, C. W. Lu, C. P. Jiang, C. C. Chuang, C. Y. Wang, and C.-W. Sun, “Subgingival calculus imaging based on swept-source optical coherence tomography,” J. Biomed. Opt. 16(7), 071409 (2011). [CrossRef]  

35. R. Birabassov, A. Yesayan, and T. Galstyan, “Nonreciprocal diffraction by spatial modulation of absorption and refraction,” Opt. Lett. 24(23), 1669–1671 (1999). [CrossRef]  

36. Y. Li, C. Shen, Y. Xie, J. Li, W. Wang, S. A. Cummer, and Y. Jing, “Tunable asymmetric transmission via lossy acoustic metasurfaces,” Phys. Rev. Lett. 119(3), 035501 (2017). [CrossRef]  

37. S. Hua, Y. M. Liu, G. E. Lio, X. J. Zhang, J. H. Wu, M. Artoni, and G. C. La Rocca, “Tailored diffraction asymmetries from spatially odd-symmetric phase gratings,” Phys. Rev. Res. 4(2), 023113 (2022). [CrossRef]  

38. S. Lisinski, J. Krause, D. Schaniel, L. Ratke, and T. Woike, “Variable refractive indices in mixed oxide xerogels,” Scr. Mater. 58(7), 553–555 (2008). [CrossRef]  

39. D. M. Wu, M. L. Solomon, G. V. Naik, A. Garcia-Etxarri, M. Lawrence, A. Salleo, and J. A. Dionne, “Chemically responsive elastomers exhibiting unity-order refractive index modulation,” Adv. Mater. 30(7), 1703912 (2018). [CrossRef]  

40. O. Ripoll, V. Kettunen, and H. P. Herzig, “Review of iterative Fourier-transform algorithms for beam shaping applications,” Opt. Eng. 43(11), 2549–2556 (2004). [CrossRef]  

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (7)

Fig. 1.
Fig. 1. $(a)$ Schematic of a phase grating of period $a$ along both $x$ and $y$ directions, which can diffract a light beam of amplitude $E_{in}$ incident along the $z$ direction into a few light beams of amplitudes $E_{i,j}$ deviating with angles $(\theta _{i},\theta _{j})$ from the $z$ direction. $(b)$ Fine structure in a period consisting of $2m_{g}\times 2n_{g}$ elements with a common thickness $\mathcal {L}$ , identical widths $\delta x=a/2m_{g}$ and $\delta y=a/2n_{g}$ , while different refractive indices $n_{m,n}=n_{0}+\delta n_{m,n}$ restricted by $\delta n_{m,n}=\delta n_{-m,n}$ and $\delta n_{m,n}=-\delta n_{m,-n}$ .
Fig. 2.
Fig. 2. Intensity $I$ against $\sin \theta _x$ and $\sin \theta _y$ with $E_{i,-5}=E_{i,1}=0$ $(a_1)$ while $E_{i,-2}=E_{i,4}=0$ ( $a_2$ ) as we set integer matrices $\hat {l}^{1}_{m,n}=\hat {l}^{-2}_{m,n}=\{-1,0,2;-2,0,1;0,1,2\}$ to attain offset refractive indices $\delta n_{m,n}^{1}$ in ( $b_1$ ) while $\delta n_{m,n}^{-2}$ in ( $b_2$ ). Other parameters used in calculations are $\lambda _{p}=0.8$ $\mu$ m, $m_{g}=3, n_{g}=3$ , $R=30$ , $M_x=M_y=10$ , $\mathcal {L}=10\lambda _{p}$ , and $n_{0}=2.0$ .
Fig. 3.
Fig. 3. Intensity $I$ against $\sin \theta _x$ and $\sin \theta _y$ with all odd columns eliminated on the basis of $E_{i,1}=0$ $(a_1)$ while all even columns eliminated on the basis of $E_{i,-2}=0$ $(a_2)$ . Parameters are the same as in Fig. 2 except $\hat {l}^{1}_{m,n}=\hat {l}^{-2}_{m,n}=\{-1,0,1;-2,0,0;0,1,2\}$ restricted by $l^{1,-2}_{m,n}+l^{1,-2}_{m,n^\ast }\in \{0,\pm 2,\pm 4,\ldots \}$ are reset to attain $\delta n_{m,n}^{1}$ in ( $b_1$ ) while $\delta n_{m,n}^{-2}$ in ( $b_2$ ).
Fig. 4.
Fig. 4. Intensity $I$ against $\sin \theta _x$ and $\sin \theta _y$ with two even columns selected on the basis of $E_{i,1}=0$ $(a_1)$ while two odd columns selected on the basis of $E_{i,-2}=0$ $(a_2)$ . Parameters are the same as in Fig. 2 except $\hat {l}^{1}_{m,n}=\hat {l}^{-2}_{m,n}=\{-1,0,1;-2,-1,0;0,1,2\}$ restricted by $l^{1,-2}_{m,n+1}-l^{1,-2}_{m,n}\in \{\pm 1,\pm 3,\ldots \}$ are reset to attain $\delta n_{m,n}^{1}$ in ( $b_1$ ) while $\delta n_{m,n}^{-2}$ in ( $b_2$ ).
Fig. 5.
Fig. 5. Intensity $I$ against $\sin \theta _x$ and $\sin \theta _y$ with all odd rows eliminated on the basis of $E_{i,2}=0$ in the presence of directional column elimination $(a_1)$ ; grouped column elimination $(a_2)$ ; directional column selection $(a_3)$ . Parameters are the same as in Fig. 2 except corresponding integer matrices restricted by $l^{2}_{m,n}-l^{2}_{m^\ast,n}\in \{0,\pm 2,\pm 4,\ldots \}$ are shown in $(b_1)$ , $(b_2)$ , and $(b_3)$ while corresponding offset refractive indices are shown in $(c_1)$ , $(c_2)$ , and $(c_3)$ .
Fig. 6.
Fig. 6. Intensity $I$ against $\sin \theta _x$ and $\sin \theta _y$ with all even rows eliminated on the basis of $E_{i,1}=0$ in the presence of directional column elimination $(a_1)$ ; grouped column elimination $(a_2)$ ; directional column selection $(a_3)$ . Parameters are the same as in Fig. 2 except $m_g=n_g=4$ and corresponding integer matrices restricted by $l^{1}_{m,n}-l^{1}_{m^\ast,n}\in \{\pm 1,\pm 3,\ldots \}$ are shown in $(b_1)$ , $(b_2)$ , and $(b_3)$ while corresponding offset refractive indices are shown in $(c_1)$ , $(c_2)$ , and $(c_3)$ .
Fig. 7.
Fig. 7. Intensity $I$ against $\sin \theta _x$ and $\sin \theta _y$ with the central row selected on the basis of $E_{i,1}=0$ in the presence of directional column elimination $(a_1)$ ; grouped column elimination $(a_2)$ ; directional column selection $(a_3)$ . Parameters are the same as in Fig. 2 except $m_g=n_g=4$ and corresponding integer matrices restricted by $l^{1}_{m+1,n}-l^{1}_{m,n}\in \{0,\pm 2,\pm 4,\ldots \}$ are shown in $(b_1)$ , $(b_2)$ , and $(b_3)$ while corresponding offset refractive indices are shown in $(c_1)$ , $(c_2)$ , and $(c_3)$ .

Tables (2)

Tables Icon

Table 1. Conditions for three types of tailored asymmetric diffraction columns.

Tables Icon

Table 2. Conditions for three types of tailored symmetric diffraction rows.

Equations (15)

Equations on this page are rendered with MathJax. Learn more.

T m , n = e ı n m , n k p L = e ı δ n m , n k p L ,
T ( x , y ) = m = m g m g n = n g n g T m , n r e c t ( x x m δ x ) r e c t ( y y n δ y ) ,
E ( θ x , θ y ) = 4 sin ( k x δ x / 2 ) sin ( k y δ y / 2 ) a 2 k x k y m = m g m g n = n g n g e ı ( δ n m , n k p L k x x m k y y n ) ,
I ( θ x , θ y ) = | E ( θ x , θ y ) | 2 [ sin ( π M x R x sin θ x ) M x sin ( π R x sin θ x ) sin ( π M y R y sin θ y ) M y sin ( π R y sin θ y ) ] 2 ,
E i , j = 16 sin ( k i δ x / 2 ) sin ( k j δ y / 2 ) a 2 k i k j m = 1 m g cos ( k i x m ) n = 1 n g cos ( δ n m , n k p L k j y n ) ,
δ n m , n j = [ j ( 2 n 1 ) + ( 2 l m , n j 1 ) n g ] λ p 4 n g L ,
E i , j ± 2 k = 16 sin ( k i δ x / 2 ) sin ( k j ± 2 k δ y / 2 ) a 2 k i k j ± 2 k m = 1 m g cos ( k i x m ) n = 1 n g sin ( l m , n j π k δ n ) ,
sin ( l m , n j π k δ n ) + sin ( l m , n j π k δ n ) 2 = sin ( l m , n j + l m , n j 2 π ) cos ( l m , n j l m , n j 2 π 2 n 1 n g k π ) = 0 ,
E i , j ± k = 4 sin ( i π / 2 m g ) sin [ ( j ± k ) π / 2 n g ] i ( j ± k ) π 2 m = 1 m g cos ( k i x m ) n = 1 n g sin ( l m , n j π k δ n 2 ) .
n = 1 n g cos ( k δ 1 2 ) sin ( l m , n j π k δ n 2 ) = 1 2 n = 1 n g [ sin ( l m , n j π k ( n 1 ) n g π ) + sin ( l m , n j π k n n g π ) ] = 0.
n = 1 n g sin ( l m , n j + l m , n + 1 j 2 π k n n g π ) cos ( l m , n + 1 j l m , n j 2 π ) = 0 ,
E q m g , j ± k = 4 sin [ ( j ± k ) π / 2 n g ] q m g ( j ± k ) π 2 m = 1 m g ( 1 ) q m sin ( q π / 2 ) cos ( q π / 2 ) n = 1 n g sin ( l m , n j π k δ n 2 ) ,
E m , n i , j ± k + E m , n i , j ± k 2 cos ( k i x m ) sin [ ( l m , n j + l m , n j ) π 2 k δ n 2 ] cos [ ( l m , n j l m , n j ) π 2 ] ,
E m , n i , j ± k + E m , n i , j ± k 2 cos ( k i x m ) sin [ ( l m , n j + l m , n j ) π 2 k δ n 2 ] sin [ ( l m , n j l m , n j ) π 2 ] ,
E i , j ± k = 4 sin ( i π / 2 m g ) sin [ ( j ± k ) π / 2 n g ] i ( j ± k ) π 2 m = 1 m g cos ( k i x m ) n = 1 n g sin ( l m g , n j π k δ n 2 ) ,
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.