Abstract
An approximate beam propagation method is proposed as an intuitive simulation of the optics of Pancharatnam–Berry phase and polarization volume hologram devices. Using this method, the connection between and polarization properties of these two types of devices are made clear.
© 2023 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement
1. INTRODUCTION
Polarization volume hologram (PVH) gratings and Pancharatnam–Berry phase devices (PPD) have wide application in optical systems. For example, a pancake optics design using a PVH has been proposed as one of the potential solutions for a lightweight headset in a virtual reality (VR) system [1]. PPD devices, proposed in [2,3], have been used in optical beam steering devices [4] and metamaterial [5–8].
These two types of devices, when considered separately, have much different optical effects and intuitive explanations as shown in Fig. 1. In the figure, as in the remaining discussion, it is assumed that the thickness of each device is close to $\lambda /({2}\Delta {n})$. The optical response of PPD devices is to have a polarization dependent steering of light between two angles symmetric about the input angle of light, which has been made clear using the Jones calculus. PVH devices are similar in their optical properties to Bragg reflectors except that light of orthogonal polarization states is either deflected to the Bragg angle, or not deflected. PPD devices have a relatively high efficiency for a wide range of incident light angles, while PVH devices have a narrower range of incident angles that yield high efficiency optical deflection.
For the given structure, Px is positive, and the grating vector is along the ${+}x$ direction, as shown in Fig. 1 with the black arrow. The grating vector angle $\phi$ makes ${+}{90}^\circ$ with the ${+}z$ axis. The angle $\tau$ describes the director rotation about the $z$ axis due to Px. ${\theta _0}$ is the incident angle in the grating. When the incident angle has projections on the ${+}x$ and ${+}z$ axes, the incident angle is positive (in Fig. 1, ${\theta _0}$ is positive). The polarization state of light represented by the red arrow is right-hand circular polarization (RCP), while the green arrow is left-hand circular polarization (LCP). For PPD, with RCP incidence the zeroth and ${-}{1}$st diffraction orders are observed. With LCP incidence, the zeroth and ${+}{1}$st diffraction orders are observed. For PVH, with RCP incidence only the zeroth diffraction order is observed. With LCP incidence, the zeroth and ${+}{1}$st diffraction orders are observed. The intensity of the two diffraction orders depends on the device thickness, while the highest diffract efficiency (${-}{1}$st or ${+}{1}$st order) is reached for the case of the device thickness of $\lambda /({2}\Delta {n})$.
Both of these effects can be quantitatively modeled using the rigorous coupled wave approach (RCWA) or by using the finite difference time domain (FDTD) method. The two methods are rigorous without approximations and give the same accurate result [9], but these computational approaches do not provide an intuitive analytical connection between the PPD and the PVH type devices, or a predictive analytical expression to explain their properties.
For isotropic materials, Gaylord and Moharam [10] considered a similar question in exploring the connection between Raman–Nath and Bragg type devices. They used an approximate coupled wave analysis to provide an intuitive derivation of the “$\rho$” parameter, which is used for distinguishing the Raman–Nath grating and Bragg grating effects in [11].
In this paper, the approach of Gaylord is extended to birefringent materials using a beam propagation method to clearly show the analytical connection between the optical properties of PPD and PVH devices.
2. APPROXIMATE COUPLED WAVE APPROACH FOR BIREFRINGENCE MATERIAL
In this work, the approach of Gaylord and Moharam for isotropic material [10] and the “$\rho$” parameter in [11] is extended to birefringence materials using a modified version of the dielectric tensor derivation of Sasaki et al. [12].
Following [10], to find the diffraction efficiency, the starting point is the Helmholtz equation:
${E},\;{ k}$, and $\varepsilon$ are the electric field, wave vector, and dielectric tensor. To use Eq. (1) to show the evolution of the amplitude and polarization states of light as it propagates in a periodic birefringent structure, light is considered to be decomposed into plane waves with defined directions corresponding to the allowed diffraction orders under the Floquet condition. Further, the periodic birefringent layer is considered to be discretized into layers that are parallel to the plane of the film. We can think the film consists of ${ N}$ layers, each with a very small thickness $d$ ($d = {D}/{N}$).
Two assumptions are used to simplify the calculation:
- 1. The $k$ vector directions of the $s$ and $p$ modes for each layer are assumed to be at the same angle.
- 2. Boundary effects are ignored.
Equations are here derived to find the amplitude and polarization state of each Floquet order in a particular layer, given the amplitude and polarization state of each of the orders in the layer before it.
The calculation steps shown in Fig. 2 that are detailed in Supplement 1, give the equation that relates the values of the complex wave amplitude for the $l$th diffraction order, ${ S}_l^{\rm{layer}}$, in a layer to those in the next layer (${\rm layer} + {1}$). Specifically, the complex wave amplitude for the $l$th diffraction order ${S}_l^{\rm layer + 1}$, in the light propagation direction, is
The details of this derivation are in Supplement 1, and the definitions of the variables in Eq. (2) are defined below.
In Eq. (2), the complex amplitude of ${{S}_l}$ in a layer is given by
${J}$ is the Jones vector describing the polarization state of light in a frame with the $z$ axis along the $k$ vector direction. ${\theta _l}$ is the diffraction angle for the $l$th diffraction order. $R({{\theta _l}})$ is the rotation matrix from the considered ray’s $k$ vector frame to the lab frame:
$\Lambda$ is the grating pitch. $\lambda$ is the vacuum wavelength of the incident beam. The value of $m$ is equal to an integer when the Bragg condition is met.
${{A}_l}$ is
Also, the components of the dielectric tensor that appear in Eq. (2) are ${\varepsilon_a,\;{{\rm M}_ +}}$, and ${\rm M} $-, where ${\varepsilon _a}$ is the part of the dielectric tensor that is not a function of the periodicity of the material (thus not a function of $\tau$).
For slanted transmissive polarization volume hologram (T-PVH) [13], it is
$\Delta {{N}^2}$ is defined as
${{\rm M}_ +}$ and ${\rm M}$- are matrices:
At the Bragg condition, the incident, deflected, and grating vector forms a vector triangle [14]. Thus, the Bragg diffraction order is ${+}{1}$st order when the incident angle is positive ($m = {1}$). The Bragg diffraction order is ${-}{1}$st order when the incident angle is negative ($m = - {1}$).
For the case of positive incident angle, and grating vector of 90°, writing out Eq. (2) for the ${-}{1}$st, zeroth, and ${+}{1}$st orders,
These equations can be used to find the output amplitude and polarization state of light exiting a film of many layers. It turns out that this simple method provides quite accurate results compared with RCWA, at least for the case in which the grating $k$ vector is in the plane of the layers ($\phi = {90}^\circ$). A comparison between the proposed beam propagation method and RCWA can be found in Supplement 1. Given the verification of this method to predict the amplitude and phase of light exiting a film of a periodic birefringent media, we can use the simple form of these equations to analyze the evolution of light going through the film as a function the parameters they contain.
3. EFFECT OF THE INCIDENT LIGHT POLARIZATION STATE ON THE RESPONSE OF PPD AND PVH DEVICES
In Eq. (2), on the right-hand side in the bracket, there are two terms containing the complex amplitude for the considered diffraction order in the previous layer (${S}_{l}^{{\rm layer}}$), and terms containing the complex amplitude for the adjacent diffraction orders (${S}_{l + 1}^{{\rm layer}}$ and ${S}_{l - 1}^{{\rm layer}}$).
There is an extra term [the term with $({{\varepsilon _a} - {n}_{{\rm ave}}^2})$ in the brackets of Eqs. (2) and (3)] as compared the equation for an isotropic material in [10]. When the grating vector is close to 90°, and light is normally incident (electric field has a small component along the $z$ direction), this term is very close to 0.
With the above assumptions, Eq. (2) becomes Eq. (4), and Eq. (3) becomes Eq. (5):
For slanted T-PVH, when the grating vector angle $\phi = {90}^\circ$, ${\rm M} +$ and ${\rm M} -$ are simplified to
When the incident polarization state is LCP, and the incident angle is close to 0° (normal incidence), Eq. (5) yields
We can see from here that for the ${+}{1}$st order beam, excited from the zeroth order LCP incident beam, is RCP.
For ${{S}_{- 1}}$,
We can see here that there is no energy flowing into the ${-}{1}$st order.
The same calculation steps can be applied for RCP incidence, in which case the ${-}{1}$st order is activated instead of the ${+}{1}$st order.
Note that when the incident angle is negative($m = \; - {1}$), the same calculation can be applied to find whether the same incident polarization state will activate the same diffraction order as when the incident angle is positive.
To summarize the polarization response for a PPD/PVH structure with positive Px,
- 1. LCP incidence can only activate the ${+}{1}$st order.
- 2. RCP incidence can only activate the ${-}{1}$st order.
4. DISTINCTION BETWEEN PPD AND PVH DEVICES
In [11], the “$\rho$” parameter is derived for an isotropic material that is the ratio of the two terms for the complex amplitude ${{S}_{l}}$. In this section, the $\rho$ parameter is extracted for a birefringent material from Eq. (2) applying a similar method as for the isotropic material. In this section a positive angle of incidence is assumed.
In Eq. (5), we can see that the change in the zeroth and ${+}{1}$st orders is related to their coupling with adjacent orders. But the change in the ${-}{1}$st order has a different dependence. It is related to its coupling with the adjacent zeroth order, and to itself.
It will be shown below that, for the ${-}{1}$st order, the ratio of the self-coupling term and the coupling to zeroth order terms is significant to differentiating the PPD and PVH optical effects. The ratio is
This ratio will be designated as the “$\rho$” parameter for birefringence material,
The effect of the value of the “$\rho$” parameter can be clearly seen by looking at the propagation of light in the birefringent periodic material though its first two layers [the first two iterations of Eq. (5)].
Light incident on the first layer only has energy in the zeroth order. With RCP incident light, using Eq. (5), the complex amplitudes of the zeroth, ${+}{1}$st, and ${-}{1}$st orders exiting the first layer are
As shown in Section 3, with the assumption that the incident angle is close to 0, energy does not flow to the ${+}{1}$st order with RCP incidence. We therefore only need to calculate the zeroth and ${-}{1}$st orders when exiting the second layer.
Exiting the second layer, for the zeroth order,
The phase of the zeroth order remains unchanged after propagating through the first two layers.
Exiting the second layer, for the ${-}{1}$st order,
The term coupling from the zeroth order in the second layer has the same phase as the ${-}{1}$st order in the first layer. However, the self-coupling term (boxed) has phase of zero that is ${-}\pi /{2}$ different from the ${-}{1}$st order in the first layer.
If the “$\rho$” parameter is small, the self-coupling term is negligible, which results in the phase of the ${-}{1}$st order being constant as it propagates from layer to layer through the material. However, if the “$\rho$” parameter is large, the self-coupling term is not negligible. In this case the phase of the ${-}{1}$st order changes when it propagates from layer to layer, and its phase is changing with respect to the other orders.
In [9], the second term (red boxed) is called a “dephasing” term. This effect can be illustrated by plotting the phase of each of the considered three orders as a function of film thickness as shown in Fig. 3.
When the “$\rho$” parameter is very small, the ${-}{1}$st order has a constant phase [Fig. 3(a)]. As the “$\rho$” parameter keeps increasing, the variation of the phase of the ${-}{1}$st order across the thickness of the device starts to be obvious [Fig. 3(c)]. When the “$\rho$” parameter is large, we can see that the phase of the ${-}{1}$st order is rapidly changing relative to the other two orders [Fig. 3(d)].
From Eq. (5), we can see that the energy is only transferred between adjacent orders. Further, if the ${-}{1}$st order does not have a constant phase with respect to the zeroth order, then as this diffraction order travels through the periodic media the energy transferred to the ${-}{1}$st order from the zeroth order is averaged out to zero. It can also be seen that relative phase between the ${+}{1}$st order and the zeroth order is approximately constant, so energy is transferred between these two orders independent of the value of the “$\rho$” parameter.
The physical explanation of this rapidly varying phase of the ${-}{1}$st order with respect to the zeroth order for the case of a large “$\rho$” value can be roughly seen by considering the simplified ray trace diagrams in Fig. 4.
For PPD, because the grating period is large, there is not expected to be a large phase difference for the ${+}{1}$st (long dashes) and ${-}{1}$st (dotted line) order beams. For PVH, with a much smaller grating period, the zeroth order deflected beam (solid line) and the ${+}{1}$st order beam traverse regions of similar index of refraction and are expected to have a similar phase variation as they propagate through the media. However, the ${-}{1}$st order sees a greater variation in the index of refraction, and therefore is expected to have a changing phase relationship to the zeroth order beam.
Note that when the incident angle is negative, $m = - {1}$, the ${+}{1}$st order has a dephasing term (instead of the ${-}{1}$st order for the positive incident angle). Following the same calculation in this section, the ${+}{1}$st order will not have a constant phase with a large “$\rho$” parameter because of the self-coupling effect illustrated.
5. EFFECT THE “$\rho$” PARAMETER AND THE TRANSITION FROM PPD TO PVH BEHAVIOR
In this section, we evaluate the effect of the “$\rho$” parameter on the PVH to PPD transition when RCP light incident on the ${+}{\rm Px}$ structure (as in Fig. 1).
In Fig. 5(a), where the value of “$\rho$” is low, it is seen, where the value of $\Delta {n}*{ D}/\;\lambda = {1/2}$, that nearly all intensity is transferred from the incident zeroth order to the deflected ${-}{1}$st order as expected for a PPD device. However, as the value of “$\rho$” increases [Figs. 5(b)–5(d)], less power is transferred from the incident zeroth order to the ${-}{1}$st order, showing the transition from PPD type optical response to PVH type optical response.
Also, for the PPD device, the oscillation of power between the allowed diffracted order and the zeroth order as a function of thickness is made clear.
6. SUMMARY AND DISCUSSION
In this work, the approximate coupled wave approach is extended from isotropic to birefringent materials. Compared with the RCWA, the approximate approach has its limitations (close to Bragg incidence and small vertical grating vector), but is more intuitively understandable.
Using the obtained equation from the beam propagation method of the solution to the wave equation, an explanation is provided for the difference of the PPD and PVH gratings. First, the incident polarization state only has the possibility to excite either the ${+}{1}$st or ${-}{1}$st order depending on its polarization state being LCP or RCP (for an assumed ${+}{\rm Px}$ structure as in Fig. 1). Second, for the case in which the incident angle and the grating $k$ vector have the same sign, the distinction between the optical effect of a PPD versus a PVH is related to the changing relative phase between the ${-}{1}$st order and the zeroth order that results in no net energy being transferred from the zeroth order to the ${-}{1}$st order if the “$\rho$” parameter is large. In the case in which the incident angle and the grating $k$ vector are of the opposite sign, the distinction between the optical effect of a PPD versus a PVH is related to the changing relative phase between the ${+}{1}$st order and the zeroth order that results in no net energy being transferred from the zeroth order to the ${+}{1}$st order if the “$\rho$” parameter is large. These results are summarized in Fig. 6 for the case in which the grating $k$ vector and incident angle are both positive.
In Fig. 6, as shown in Section 3, energy transfer to the order indicated by the black arrow is prevented because the order is not activated by the incident polarization state. The difference between PPD and PVH devices was shown in Section 4. For the PVH device, the energy transfer to the ${-}{1}$st order is prevented because of the rapid changing phase relationship between this order and the zeroth order incident beam. The polarization explanation as shown in Section 3 and the phase explanation as shown in Section 4 combine to give us a complete answer to why Fig. 1 is observed in reality.
7. LIMITATION OF THE WAVE PROPAGATION SIMULATION
Unlike FDTD and RCWA, the proposed beam propagation method is very transparent in showing the intensity and phase change of each diffraction order when propagating through an anisotropic media. However, the drawback is that the proposed beam propagation method is not as accurate as FDTD and RCWA because of the necessary approximations that we used to arrive at a simple expression of complex amplitude.
The first approximation is that the $k$ vector directions of the $s$ and $p$ modes for each layer are assumed to be at the same angle.
This limits the simulation to be less accurate for large birefringence. In future work, a more accurate beam propagation method can be worked on by studying the $s$ and $p$ modes separately.
The second approximation is the boundary effects are ignored. This means that the reflections at the interfaces are not considered.
The third approximation is that the term with $({{\varepsilon _a} - {n_{{\rm ave}}^2}}){S}_l^{{\rm layer}}$ in the brackets of Eqs. (2) and (3) is assumed to be 0.
When the $\phi = {90}^\circ$, as shown in Supplement 1, the expression simplifies to
This shows that the third approximation yields a simulation result that is less accurate for large birefringence, large incident angle (nonzero ${Sz}$), and when the grating vector is different from 90°.
Funding
Meta.
Disclosures
The authors declare no conflicts of interest.
Data availability
Data underlying the results presented in this paper are available in Ref. [15].
Supplemental document
See Supplement 1 for supporting content.
REFERENCES
1. A. Maimone and J. Wang, “Holographic optics for thin and lightweight virtual reality,” ACM Trans. Graph. 39, 67 (2020). [CrossRef]
2. S. Pancharatnam, “Generalized theory of interference and its applications,” Proc. Indian Acad. Sci. A 44, 398 (1956). [CrossRef]
3. M. V. Berry, “Quantal phase factors accompanying adiabatic changes,” Proc. R. Soc. Lond. A 392, 45 (1984). [CrossRef]
4. C. Yousefzadeh, A. Van Rynbach, and P. J. Bos, “Tunable liquid crystal beam steering device based on Pancharatnam phase,” Proc. SPIE 11303, 113030S (2020). [CrossRef]
5. Y. Guo, M. Pu, F. Zhang, M. Xu, X. Li, X. Ma, and X. Luo, “Classical and generalized geometric phase in electromagnetic metasurfaces,” Photon. Insights 1, R03 (2022). [CrossRef]
6. Q. Xu, X. Su, X. Zhang, L. Dong, L. Liu, Y. Shi, Q. Wang, M. Kang, A. Alù, S. Zhang, J. Han, and W. Zhang, “Mechanically reprogrammable Pancharatnam–Berry metasurface for microwaves,” Adv. Photon. 4, 016002 (2022). [CrossRef]
7. X. Xie, M. Pu, J. Jin, M. Xu, Y. Guo, X. Li, P. Gao, X. Ma, and X. Luo, “Generalized Pancharatnam-Berry phase in rotationally symmetric meta-atoms,” Phys. Rev. Lett. 126, 183902 (2021). [CrossRef]
8. X. Zhang, S. Yang, W. Yue, Q. Xu, C. Tian, X. Zhang, E. Plum, S. Zhang, J. Han, and W. Zhang, “Direct polarization measurement using a multiplexed Pancharatnam–Berry metahologram,” Optica 6, 1190–1198 (2019). [CrossRef]
9. X. Xiang and M. J. Escuti, “Numerical modeling of polarization gratings by rigorous coupled wave analysis,” Proc. SPIE 9769, 976918 (2016). [CrossRef]
10. T. Gaylord and M. Moharam, “Planar dielectric grating diffraction theories,” Appl. Phys. B 28, 1–14 (1982). [CrossRef]
11. M. G. Moharam and L. Young, “Criterion for Bragg and Raman-Nath diffraction regimes,” Appl. Opt. 17, 1757–1759 (1978). [CrossRef]
12. T. Sasaki, K. Miura, O. Hanaizumi, A. Emoto, and H. Ono, “Coupled-wave analysis of vector holograms: effects of modulation depth of anisotropic phase retardation,” Appl. Opt. 49, 5205–5211 (2010). [CrossRef]
13. X. Feng, L. Lu, O. Yaroshchuk, and P. Bos, “Closer look at transmissive polarization volume holograms: geometry, physics, and experimental validation,” Appl. Opt. 60, 580–592 (2021). [CrossRef]
14. H. Kogelnik, “Coupled wave theory for thick hologram gratings,” Bell Syst. Tech. J. 48, 2909–2947 (1969). [CrossRef]
15. F. Xi and P. Bos, “Intuitive understanding of the connection between Pancharatnam-Berry optical beam deflectors and polarization volume holograms: supplemental document,” figshare, (2023), https://doi.org/10.6084/m9.figshare.22015922.