Abstract
The application of the adiabatic geometric phase (AGP) to nonlinear frequency conversion may help to develop new types of all-optical devices, which leads to all-optical modulation of the phase front of one wave by the intensity of other waves. In this paper, we develop the canonical Hamilton equation and a corresponding geometric representation for two schemes of four-wave mixing (FWM) processes (ω1 + ω2 = ω3 + ω4 and ω1 + ω2 + ω3 = ω4), which can precisely describe and calculate the AGP controlled by the quasi-phase matching technique. The AGPs of the idler (ω1) and signal (ω4) waves for these two schemes of FWM are studied systematically when the two pump waves (ω2 and ω3) are in either the undepleted or in the depleted pump cases, respectively. The analysis reveals that the proposed methods for calculating the AGP are universal in both cases. We expect that the analysis of AGP in FWM processes can be applied to all-optically shaping or encoding of ultrafast light pulse.
© 2021 Optical Society of America under the terms of the OSA Open Access Publishing Agreement
1. Introduction
When an eigenstate slowly evolves around a closed trajectory in the parameter space by varying its Hamiltonian, an adiabatic geometric phase (AGP) factor is acquired in addition to the dynamical phase, which is given by the time integral of the eigenenergy [1]. Moreover, the AGP can be interpreted as the flux of a virtual field through the surface enclosed by the closed circuit in the parameter space [2]. Over the past couple of decades, important applications of the AGP have been demonstrated in quantum physics [3], for example, in high-precision quantum measurement [4], quantum information processing [5], quantum computing [6], Bose-Einstein condensates [7,8] and condensed-matter physics [9]. AGPs also have extensive applications in optics, in which they are used for beam shaping [10,11], beam shearing [12], guiding [13], holography [14] and encoding both classical [15] and quantum [16] information.
Recently, the geometric phase for $\chi ^{(2)}$ nonlinear frequency conversion was studied [17–19]. Some manifestations of the geometric phase for light were obtained by changing the polarization or the photon direction in nonlinear metasurfaces [20–24]. A different case is when the AGP is accumulated in a bulk material by nonlinear modulation pattern along the nonlinear crystal through the quasi-phase matching (QPM) technique. The analysis in literatures [18,19] provides an effective geometric representation of the latter case and a method to describing and calculating the AGP in a precise way. This method is universal for the undepleted pump (the intensity of the pump wave is much stronger than the intensities of the weak signal and idler) and depleted pump (the intensity of the pump wave is comparable to the intensities of the signal and idler) cases. Moreover, nonreciprocal transmission and wavefront shaping through AGPs were realized experimentally by circular rotation of the QPM parameters using sum-frequency generation (SFG) in the $\chi ^{(2)}$ process [25]. These methods provide new ways to manipulate AGP by nonlinear optics.
Adiabatic conversion in nonlinear process is not limited only to three wave mixing processes. Recently, adiabatic four wave mixing was theoretically studied, through either a longitudinally varying silicon waveguide structure or a linearly tapered fluoride step-index fiber [26]. It was assumed that two of the four waves serve as undepleted pump waves, so that the adiabatic exchange of energy occurs between the two remaining waves. Very recently, such an undepleted adiabatic passage in FWM was experimentally observed in a photonic crystal fiber with a tapered core diameter [27]. Here we would like to address the question of how to accumulate geometric phase in a four wave mixing process and how it can be calculated. Unlike Refs. [26,27], our analysis is not limited to the case of undepleted pump waves. In order to analyze the systems dynamics, we develop a canonical Hamiltonian equation and a corresponding geometrical representation of all the system states on a closed surface in three dimensions. This surface is the Bloch sphere in the undepleted case, but for the depleted case the shape becomes elongated and one of the poles attains a sharp conical shape if the system reaches the full depletion condition.
The remainder of this paper is structured as follows. The theory and model for the QPM in the FWM process are illustrated in Sec. II; two schemes of the FWM process ($\omega _{1}+\omega _{2}=\omega _{3}+\omega _{4}$ and $\omega _{1}+\omega _{2}+\omega _{3}=\omega _{4}$) are considered in this section. The AGPs in schemes I and II for the undepleted, and then depleted pump cases are discussed in Sec. III and IV, respectively. The main results of this paper are summarized in Sec. V.
2. Model and theory
Generally, the adiabatic process for nonlinear wave conversion requires to control the wave-vector mismatch parameter. Reference [26] introduces a small periodic change in the waveguide width with a longitudinal chirp to the period to achieve this parameter, as illustrated in Fig. 1(a), where the local period $\Lambda (z)$ adds a quasi momentum $K_g(z)=2\pi /\Lambda (z)$ to the wave-vector difference between the guided modes of the interacting waves. This type of corrugated waveguide can be readily fabricated using planar integrated optics technologies. By introducing a duty cycle and its relative position within the periodic cycle (i.e., the phase), three QPM parameters, namely, the modulation wave vector, modulation phase, and duty cycle, are introduced in this setting. In this paper, two schemes of the FWM process are considered; their energy diagrams are displayed in Figs. 1(b,c). Following the definition in Ref. [26], we also employ waves $\omega _{1}$ and $\omega _{4}$ to denote the idler and signal, respectively, and we use waves $\omega _{2}$ and $\omega _{3}$ to denote the two pumps, respectively. Moreover, we discuss not only the AGP in the undepleted case, i.e., $I_{2,3}\gg I_{1,4}$ (where $I_{j}$ is the intensity of wave $\omega _{j}$), but also the AGP in the depleted case, i.e., $I_{2,3}\sim I_{1,4}$.
2.1 Dimensionless equation for the FWM process of scheme I
The energy diagram for the FWM process of scheme I is shown in Fig. 1(b), which satisfies $\omega _{1}+\omega _{2}=\omega _{3}+\omega _{4}$. The coupled amplitude equations for this process with QPM [28] are
If an adiabatic modulation is added to $d(z)$, as mentioned above, periodically changing the width of the waveguide [26] [see Fig. 1(a)], $d(z)$ can be expanded by a Fourier series with slowly varying components:
The following definitions are applied:
2.2 Dimensionless equation for the FWM process of scheme II
The energy diagram for the FWM process of scheme II is illustrated in Fig. 1(c), which satisfies $\omega _{4}=\omega _{1}+\omega _{2}+\omega _{3}$. Therefore, the coupled amplitude equations for this scheme are
Equations (19) can also be presented via a canonical Hamiltonian structure, $dq_{j}/d\tau =2i\left (\partial H/\partial q^{\ast }_{j}\right )$, where
2.3 Geometrical representation and calculation of the AGP
To depict the geometric motion of Eqs. (10) and (19) in the state space, we introduce a set of coordinates $(X,Y,Z)$ with the surface $\varphi =0$ for these two schemes, which adopt the similar rule of three-wave mixing in Refs. [19,29,31] (using the real and the imaginary parts of product of the waves in the interacting part of the Hamiltonian as the coordinates $X$ and $Y$ , respectively, and using the intensity of one wave as the coordinate $Z$):
For the QPM control parameters, we can still use the design of the circular rotation of $[\Gamma (\tau ), \Xi (\tau ),\phi _{d}(\tau )]$ in Refs. [18,19]:
Here, the symbols “$+$” and “$-$” correspond to the clockwise and counterclockwise rotation, respectively, of the QPM vector, expressed as $\mathbf {Q}=(\Xi \cos \phi _{d}, \Xi \sin \phi _{d}, \Delta \Gamma )$. $\Theta$ is the angle of the normal vector $\hat {\mathbf {\Theta }}$ of the plane containing the circular trajectory of vector $\mathbf {Q}$ with respect to the vertical axis [see Fig. 1(d)]. Hence, this circular rotation is fully characterized by the angle $\Theta$. $\Omega =2\pi /T$ is the angular velocity, which denotes the speed of the rotation (here, $T=\eta L$ is the scaled length of the waveguide). For convenience, we fix $\Gamma _0=1$ in the numerical simulation. Moreover, to ensure the rotation satisfies the adiabatic condition, we select $T=1000$, which is a sufficient scale for $\Omega$ to perform an adiabatic rotation.Following the framework of three-wave mixing [19], the AGP in the QPM process for wave $\omega _j$ can be generally expressed by
where $\Phi ^{\pm }_{j}$ is the total phase shift of the wave $\omega _j$ for the clockwise and counter-clockwise rotation, respectively, $\mathcal {D}_{j}$ is the dynamical phase, which equals to the circular integral of the eigenenergy along the rotation, $\gamma _{j}$ is the acquired geometric phase. Hence, the AGP can be calculated by two equivalent methods. The first one is based on the half of the difference of the total phase between the clockwise and the counter-clockwise rotation of the vector $\mathbf {Q}$. The expression of this method is written as where $\Phi ^{\pm }_{j}$ are the total phase shift of wave $\omega _{j}$ for the clockwise and counter-clockwise circular rotation, respectively. The second one equals to the difference between the total phase accumulated along the circular trajectory and that along its vertical projection. The expression of this method is written as3. Numerical results for the AGP in the FWM process of scheme I
3.1 Undepleted limit with $I_{1}\ll I_{2}$
An undepleted case occurs at
at which the intensities of the idler and signal are much smaller than those of the two pumps. If we assume $I_4=0$, the dynamics starts from the bottom of the Bloch sphere. We first discuss the simplest case, $I_{2}=I_{3}$, at which the two pumps have a the same intensity. According to the conservative condition in Eq. (8), the limit of this condition is $\tilde {q}^{\ast }_{3}\tilde {q}_{2}=\tilde {q}^{\ast }_{2}\tilde {q}_{3}=|\tilde {q}_{2}|^{2}=|\tilde {q}_{3}|^{2}\approx 1/2$, and the magnitudes of $|\tilde {q}_{1}|^{2}$ and $|\tilde {q}_{4}|^{2}$ can be neglected. Hence, Eqs. (10) are simplified to the spin-$1/2$ Bloch equation asA comparison between $\gamma _{l}$ and the numerical $\gamma _{1}$ for different values of $\Theta$ is shown in Fig. 2(a), which demonstrates that the numerical results are in accordance with the theoretical prediction. Hence, the two methods for calculating the AGP in three-wave mixing remain valid for the case of FWM. Typical examples for 3 types of trajectories on the Bloch sphere, namely, clockwise, counterclockwise, and round-trip motion, are displayed in Fig. 2(b). The AGP accumulation calculated via Eq. (33) is compared with that calculated via Eq. (34) in Fig. 2(c), which demonstrates that these two methods remain equivalent to each other for the case of FWM. Hence, we will use Eq. (34) to calculate the AGP accumulations in the following discussion.
However, if we assume $I_{2}\neq I_{3}$, a factor, which is the modulus of $\tilde {q}^{\ast }_{2}\tilde {q}_{3}$, is multiplied in front of $g(\tau )$. This factor transforms the Bloch sphere into an ellipse and changes trajectory of the state vector. Figure 2(d) shows a typical example of the comparison between the Bloch spheres as well as the trajectories of the state vector at the case of $I_{2}=I_{3}$ and $I_{2}\neq I_{3}$. These comparison indicates that the elliptical deformation the Bloch sphere and the trajectory of the state vector create a phase shift, i.e., $\Delta \gamma _{1}=\gamma _{1}-\gamma _{l}$ on the prediction of $\gamma _{l}$. Figure 2(f) shows that this additional phase shift at the case of $I_{2}\approx I_{3}/3$ as a function of $\Theta$. It shows that the maximum additional phase shift is appears at $\Theta \approx \pi /3$.
3.2 Partial depletion with $I_{1}<I_{2}$
Depletion of the system in scheme I occurs when $I_{1}$ and $I_{2}$ become comparable. Under this circumstance, the AGP at $\omega _1$ plays a dominant role when $I_{1}<I_{2}$. If the initial condition switches to $I_{1}>I_{2}$, the AGP at $\omega _{2}$ plays a dominant role. For convenience, we will select the condition $I_{1}<I_{2}$ for consideration. Here, the full depletion occurs at $I_{1}=I_{2}$, which we will discuss in the next subsection. A typical example for the trajectories of the state vector on the Bloch sphere, i.e., a geometric representation, in this case are displayed in Fig. 3(a), and the corresponding accumulations of the geometric phases for $\tilde {q}_{1,2,3}$ are displayed in in Fig. 3(b). These panels indicate that the Bloch sphere remains a smooth and symmetric surface and the magnitude of the geometric phase shift is mainly accumulated on $\gamma _{1}$.
Figure 4 displays the phase shifts of the different waves $\gamma _{1,2,3}$ versus $\Theta$ with different ratios of $(I_{1},I_{2},I_{3})$. Figure 4(a1) shows $\gamma _{1}(\Theta )$ upon fixing $I_{2}$ and different values of $I_{1}$. At $\Theta =0$ and $\pi /2$, $\gamma _{1}=\gamma _{l}$. In the region of $0<\Theta <\pi /2$, $\gamma _{1}$ deviates from $\gamma _{l}$ if the system becomes depleted. Moreover, at a typical angle of $\Theta$, namely, $\Theta _{\mathrm {c}}$, the AGP becomes independent of the depletion. Hence, in the region of $0<\Theta <\pi /2$, $\gamma _{1}$ satisfies
Table 3 shows that $\Delta$ increases as $I_{1}$ increases, which indicates that the deviation between $\gamma _{1}$ and $\gamma _{l}=-\pi (1-\cos \Theta )$ (the phase in the undepleted case with two equally pumps) increases as the depletion increases.
Figure 4(a2) shows $\gamma _{1}(\Theta )$ with fixed values of $I_{1}$ and different proportions of $(I_{2},I_{3})$. The results indicate that the ratios of $(I_{2},I_{3})$ also influence the deviation rate between $\gamma _{1}$ and $\gamma _{l}$. Finally, as mentioned above, Fig. 4(b1,b2,c1,c2) shows that the AGPs can also accumulate for pump waves $\tilde {q}_{2,3}$, but their absolute values, i.e., $|\gamma _{2,3}|$, are much smaller than $|\gamma _{1}|$.
3.3 Full depletion with $I_{1}=I_{2}$
At the limit of $I_{1}=I_{2}=I_{\mathrm {P}}$, the system is fully depleted. At this limit, $\gamma _{1}=\gamma _{2}$, and the Bloch sphere becomes elongated to an asymmetric surface, where one of the poles is transferred into a cone. We find a threshold angle of $\Theta$, $\Theta _{\mathrm {th}}$, such that the AGP can be successfully accumulated only if $\Theta <\Theta _{\mathrm {th}}$. Typical examples of the geometric representation and of the AGP accumulation of $\gamma _{1}$ with $\Theta <\Theta _{\mathrm {th}}$ and $\Theta >\Theta _{\mathrm {th}}$ are shown in Figs. 5(a1,a2) and (b1,b2), respectively. When $\Theta >\Theta _{\mathrm {th}}$, the cusp effect is induced, which is similar to the behavior that was observed in three-wave mixing [15]. Under this circumstance, the AGP fails to accumulate. The numerical results in Fig. 5(c) indicates that $\Theta _{\mathrm {th}}$ is a function of $I_{\mathrm {P}}$.
So far we considered the initial condition $I_4=0$, but if instead we start with $I_1=0$, this would mean that rather than starting from the bottom of the Bloch sphere, we start from the top of the sphere. The accumulated geometric phase at $\omega _4$ will have an opposite sign with respect to what was calculated so far for $\omega _1$. We note that a similar phenomenon of sign inversion (between the signal wave, vs the idler was) was observed in a circular rotation in three wave mixing process [18].
4. Numerical results for the AGP in the FWM process of scheme II
4.1 Undepleted limit
Similarly, an undepleted case occurs at the initial condition of $I_{1,4}\ll I_{2,3}$. Here, we consider first the undepleted case with two equally strong pumps. According to the normalization in Eq. (20), the limits of this condition are $I_{2}=I_{3}\approx 1/2$ and $I_{1}=I_{4}\approx 0$. Under this circumstance, $q_{2,3}(\tau )\approx \sqrt {I_{2,3}}e^{i\phi _{2,3}(\tau )}$ and $|q_{1,4}|^{2}$ can be neglected. Hence, Eqs. (19) can be simplified as
According to the above descriptions, in the case of $\Delta _{s}=\phi _{P}=0$, Eq. (42) driven by Eqs. (29)–(31) can give rise to the AGP as $\gamma _{l}=\mp \pi (1-\cos \Theta )$ for idler and signal waves (i.e., $\tilde {q}_{1}$ and $\tilde {q}_{4}$) if we adopt $I_{4}=0$ and $I_{1}=0$, respectively.
In comparison to to these the two theoretical results, $\gamma _{l}=\mp (1-\cos \Theta )$, we also select the following two initial conditions for the numerical calculation:
Obviously, these two initial conditions satisfy the undepleted condition. In the numerical simulation, we use Eq. (24) to define all the Kerr coefficients in Eqs. (19). Typical examples of the geometric representation and the accumulations of the geometric phase for the dynamics being initialed from these two conditions are displayed in Figs. 6(a1,a2) and (b1,b2), respectively. They show that the geometric phases calculated via Eqs. (33) and (34) are perfectly in accordance with each other, which demonstrates the validity of the calculation. Hence, we plot the results of the AGPs for the idler and signal (i.e., $\gamma _{1,4}$) versus $\Theta$ for these two conditions in Fig. 6(c). Unlike the AGP in the FWM process of scheme I, which can perfectly match the theoretical prediction at the undepleted limit with equal pump waves, the AGPs in the FWM process of scheme II show little deviation from the theoretical prediction. If we eliminate all the Kerr coefficients, i.e., if we set $s_{jk}=0$ ($j,k=1,2,3,4$), such deviations can be removed and $\gamma _{1,4}$ can perfectly match the theoretical prediction, $\gamma _{l}=\mp \pi (1+\cos \Theta )$, again. It demonstrates that these deviations originate from the nonzero values of $\Delta _{s}$ and $\phi _{P}$, which are owing to the Stark shift..4.2 FWM with depletion
4.2.1 From the bottom of the Bloch sphere ($I_{4}=0$): sum-frequency generation
First, we will consider the accumulation of the AGP of the idler rotated from the bottom of the Bloch sphere. This process is termed SFG because the initial signal is zero (i.e., $I_{4}=0$). and the three waves at $\omega _{1,2,3}$ sum up to generate the higher frequency wave at $\omega _{4}$. When the magnitudes of the intensities of the other three inputs, i.e., $I_{1,2,3}$, are comparable to each other, the system becomes depleted. For simplification, we only study the following two cases:
In the first case [cf. Equation (47)], we can focus mainly on $\gamma _1$ because the AGP accumulates primarily at the idler, and the contributions from the two pumps (i.e., $\gamma _{2}$ and $\gamma _{3}$) can be neglected. In the second case [cf. Equation (48)], $\gamma _{1}=\gamma _{2}$, while $\gamma _{3}$ can be neglected. For these two cases, the limits of full depletion both appear at $I_{1}=I_{2}=I_{3}=1/3$, and the domain under consideration is $0<I_{1}\leq 1/3$. If $I_{1}>1/3$, the AGP accumulates mainly on one of the pump wave, which owns a smaller initial intensity.In these two cases, the depletion condition makes the curves of $\gamma _{1}(\Theta )$ deviate from their counterparts at the undepleted limit. An example for the geometric representation for the case in Eq. (47) is displayed in Fig. 7(a1). It shows that the Bloch sphere is elongated to an asymmetric egg-shape surface when the system becomes depleted.
Figure 7(a2) displays some examples for $\gamma _{1}(\Theta )$ with different values of $I_{1}$ for the case of $I_{1}<I_{2}=I_{3}$ and shows that the degree of deviation increases as the level of depletion increases. Note that the level of depletion depends on the value of $I_{1}$ (here, $0\leq I_{1}\leq 1/3$, where $I_{1}\rightarrow 0$ is the undepleted limit and $I_{1}=1/3$ is the fully depleted limit). For the case of $I_{1}<I_{2}=I_{3}$, when $I_{1}>0.2$, a threshold value of $\Theta$ appears, i.e., $\Theta _{\mathrm {th}}$. When $\Theta >\Theta _{\mathrm {th}}$, the cusp effect is induced, as the AGP fails to accumulate under these circumstances. A typical example of a geometric representation for such a cusp effect is displayed in Fig. 7(b1). It shows that the spiral lines appear after the state vector passes the cusp, which indicates that the state vector cannot follow the motion of the vector of the QPM parameters after they pass through the cusp of the sphere. Figure 7(b) displays $\Theta _{\mathrm {th}}$ versus $I_{1}$ in the domain of $I_{1}\in (0.2,1/3]$. In contrast, for the case of $I_{1}=I_{2}<I_{3}$ [c.f. Equation (48)], $\Theta _{\mathrm {th}}$ exists for the entire domain of $I_{1}$ (i.e., $I_{1}\in (0,1/3]$). The cusp effect is induced if $\Theta >\Theta _{\mathrm {th}}$. A typical example of the geometric representation for the cusp effect in this case is shown in Fig. 7(c1) and $\Theta _{\mathrm {th}}$ versus $I_{1}$ is displayed in Fig. 7(c2).
4.2.2 From the top of the Bloch sphere ($I_{1}=0$): difference-frequency generation
In this subsection, we consider the motion of the state vector from the top of the Bloch sphere. In this case, we need to set one value of $I_{1,2,3}$ to zero (here, for convenience, we set the initial intensity of the idler to zero, i.e., $I_{1}=0$). Under this circumstance, the system is in the process of difference-frequency generation, since a low frequency signal at $\omega _1=\omega _4-(\omega _3+\omega _2)$ is generated. The depletion for this process is induced when the initial intensity of the signal becomes comparable to the intensities of the pumps. Similar to the case of three-wave mixing, once the initial intensity is set at the top of the Bloch sphere, the AGP accumulates only at the signal wave.
Similarly, we also consider two cases in this subsection:
This study finds that the AGP of the signal can be acquired in the domain of $0<I_{4}\leq 1/3$. Figures 8(a) and (b) display the curves for $\gamma _{4}(\Theta )$ with different values of $I_{4}$ for these two cases. Similar to the case study in SFG, the depletion creates a deviation of the AGP from the undepleted limit. Here, we find that this depletion causes $\gamma _{4}$ to acquire a larger phase magnitude than the undepleted counterpart when $\Theta <\sim 0.44\pi$. Typical examples of the 3 types of trajectories on the Bloch sphere are shown in Fig. 9(a1), and the calculations to this $\gamma _{4}$ by Eqs. (33) and (34) are shown in Fig. 9(a2). These two panels demonstrate the validation of the related calculations.However, when $I_{4}>1/3$, the state vector $\mathbf {W}$ cannot immediately follow the motion of the QPM vector $\mathbf {Q}$, therefore an angle is created between them. This angle induces a precession of $\mathbf {W}$ around $\mathbf {Q}$, which can be approximately described by
Obviously, if $\mathbf {W}$ can immediately follow the motion of $\mathbf {Q}$, the angle $\sim 0$, which yields $\mathbf {W}\times \mathbf {Q}\sim 0$, and the precession of $\mathbf {W}$ around $\mathbf {Q}$ is suppressed. This is what we observed at the case of $I_{4}<1/3$. On the other hand, vector $\mathbf {Q}$ itself performs a circular rotation (or round-trip motion) and drags the vector $\mathbf {W}$ follow with it. Hence, when $I_{4}>1/3$, the combination of these two types of motion on $\mathbf {W}$ (rotating around $\mathbf {Q}$ and following the motion of $\mathbf {Q}$) draws a spiral trajectory on the Bloch sphere. A typical example of this phenomenon is displayed in Fig. 9(b1). Moreover, this spiral motion creates strong oscillations in the accumulation of $\gamma _{4}$, as illustrated in Fig. 9(b2). In this panel, it shows that the amplitude of the oscillation of $\gamma _{4}$ gradually increases along $\tau$, which indicates that the angle between $\mathbf {W}$ and $\mathbf {Q}$ is also gradually increase with $\tau$. Finally, if $I_{4}$ continues to increase, the magnitude of the angle between $\mathbf {W}$ and $\mathbf {Q}$ is also increased. The spiral motion of $\mathbf {W}$ becomes very strong and eventually destroys all 3 types of motion on the Bloch sphere. Hence, the accumulation of $\gamma _{4}$ is also completely destroyed under this circumstance [see an example in Figs. 9(c1,c2)].5. Conclusion
In this paper, we study the accumulation of the adiabatic geometric phase through quasi-phase matching (QPM) for two schemes of a nonlinear four-wave mixing process. A geometric representation and a canonical Hamiltonian equation are built to describe the dynamics of these two processes. The parameters of the Hamiltonian are presented by three QPM parameters, namely, the period, phase, and duty cycle of the nonlinear modulation pattern. Furthermore, the geometric representation maps all the states of the system onto a close surface. The circular rotation of these QPM parameters is characterized by an angle, $\Theta$. Under adiabatic conditions, the state vector follows the rotation of the QPM vector and draws a circular trajectory on the state surface. The AGPs of these waves are equal to the half of the difference of the total phase between the clockwise and the counter-clockwise rotation of the vector $\mathbf {Q}$. An identical result is obtained by the difference between the total phase accumulated along the circular trajectory and that accumulated along its vertical projection. We note that, in this respect, the calculation of the geometric phase follows the same rules as the case of three-wave mixing [19].
We studied two different schemes of four wave mixing, and the main results are summarized in Tables 1 and 2. Under the undepleted condition, where the two pump waves are much stronger than the idler and signal waves, the system of scheme I ($\omega _{1}+\omega _{2}=\omega _{3}+\omega _{4}$) is completely in accordance with the spin-$1/2$ dynamics, and the AGPs for the idler and signal can be predicted analytically. However, for the system of scheme II ($\omega _{1}+\omega _{2}+\omega _{3}=\omega _{4}$), the existence of the Stark shift and the phase shift of the pump waves induce a deviation in the AGP from the theoretical result of the spin-$1/2$ system. In the depleted case, where the idler and signal waves become comparable to the pump waves, the magnitudes of the AGPs for the idler and signal deviate from their counterparts at undepleted limits. In addition, the deviation ratio increases as the level of depletion increases. In the case of full depletion, the threshold angle of $\Theta$ is found, beyond which the system becomes unstable and the accumulation of the AGP fails. The ability to control the geometric phase can find interesting applications in photonic devices. One possible application is that of an all-optical phased array device for fast steering and focusing, where the array elements are a set of parallel waveguides, and the phase of of each waveguide can be all-optically controlled by varying the pump intensities.
Funding
Israel Science Foundation (1415/17); National Natural Science Foundation of China (11874112, 11974146); Natural Science Foundation of Guangdong Province (2017B030306009).
Acknowledgments
This work was supported by the Israel Science Foundation (ISF), grant no. 1415/17, by the National Natural Science Foundation of China through Grant No. 11874112, 11974146, by the Natural Science Foundation of Guangdong Province (China) through Grant No. 2017B030306009.
Disclosures
The authors declare no conflicts of interest.
References
1. M. V. Berry, “Quantal phase factors accompanying adiabatic changes,” Proc. R. Soc. Lond. A 392(1802), 45–57 (1984). [CrossRef]
2. J. Liu, S. Li, L. Fu, and D. Ye, “Changes in atmospheric consituents and in radiative forcing,” in Nonlinear Adiabatic Evolution of Quantum Systems (Springer Nature Singapore Pte Ltd, 2018).
3. B. Wu, J. Liu, and Q. Niu, “Geometric phase for adiabatic evolutions of general quantum states,” Phys. Rev. Lett. 94(14), 140402 (2005). [CrossRef]
4. E. Martín-Martínez, A. Dragan, R. B. Mann, and I. Fuentes, “Berry phase quantum thermometer,” New J. Phys. 15(5), 053036 (2013). [CrossRef]
5. P. J. Leek, J. M. Fink, A. Blais, R. Bianchetti, M. Göppl, J. M. Gambetta, D. I. Schuster, L. Frunzio, R. J. Schoelkopf, and A. Wallraff, “Observation of berry’s phase in a solid-state qubit,” Science 318(5858), 1889–1892 (2007). [CrossRef]
6. J. Xu, S. Li, T. Chen, and Z. Xue, “Nonadiabatic geometric quantum computation with optimal control on superconducting circuits,” Front. Phys. 15(4), 41503 (2020). [CrossRef]
7. S. Li, L. Fu, and J. Liu, “Adiabatic geometric phase for a bose-einstein condensate coupled to a cavity,” Phys. Rev. A 84(5), 053610 (2011). [CrossRef]
8. S. Li, J. Liu, and L. Fu, “Berry phase and hannay angle of an interacting boson system,” Phys. Rev. A 83(4), 042107 (2011). [CrossRef]
9. D. Xiao, M. Chang, and Q. Niu, “Berry phase effects on electronic properties,” Rev. Mod. Phys. 82(3), 1959–2007 (2010). [CrossRef]
10. Z. Bomzon, V. Kleiner, and E. Hasman, “Pancharatnam-berry phase in space-variant polarization-state manipulations with subwavelength gratings,” Opt. Lett. 26(18), 1424 (2001). [CrossRef]
11. L. Marrucci, C. Manzo, and D. Paparo, “Pancharatnam-berry phase optical elements for wave front shaping in the visible domain: switchable helical mode generation,” Appl. Phys. Lett. 88(22), 221102 (2006). [CrossRef]
12. A. Alemán-Castaneda, B. Piccirillo, E. Santamato, L. Marrucci, and M. A. Alonso, “Shearing interferometry via geometric phase,” Optica 6(4), 396 (2019). [CrossRef]
13. S. Slussarenko, A. Alberucci, C. P. Jisha, B. Piccirillo, E. Santamato, G. Assanto, and L. Marrucci, “Guiding light via geometric phases,” Nat. Photonics 10(9), 571–575 (2016). [CrossRef]
14. T. Zhan, J. Xiong, Y.-H. Lee, R. Chen, and S.-T. Wu, “Fabrication of pancharatnam-berry phase optical elements with highly stable polarization holography,” Opt. Express 27(3), 2632 (2019). [CrossRef]
15. G. Milione, T. A. Nguyen, J. Leach, D. A. Nolan, and R. R. Alfano, “Using the nonseparability of vector beams to encode information for optical communication,” Opt. Lett. 40(21), 4887 (2015). [CrossRef]
16. T. Stav, A. Faerman, E. Maguid, D. Oren, V. Kleiner, E. Hasman, and M. Segev, “Quantum entanglement of the spin and orbital angular momentum of photons using metamaterials,” Science 361(6407), 1101–1104 (2018). [CrossRef]
17. K. Wang, Y. Shi, A. S. Solntsev, S. Fan, A. A. Sukhorukov, and D. N. Neshev, “Non-reciprocal geometric phase in nonlinear frequency conversion,” Opt. Lett. 42(10), 1990 (2017). [CrossRef]
18. A. Karnieli and A. Arie, “Fully controllable adiabatic geometric phase in nonlinear optics,” Opt. Express 26(4), 4920 (2018). [CrossRef]
19. Y. Li, O. Yesharim, I. Hurvitz, A. Karnieli, S. Fu, G. Porat, and A. Arie, “Adiabatic geometric phase in fully nonlinear three-wave mixing,” Phys. Rev. A 101(3), 033807 (2020). [CrossRef]
20. M. Tymchenko, J. S. Gomez-Diaz, J. Lee, N. Nookala, M. A. Belkin, and A. Alù, “Gradient nonlinear pancharatnam-berry metasurfaces,” Phys. Rev. Lett. 115(20), 207403 (2015). [CrossRef]
21. G. Li, S. Chen, N. Pholchai, B. Reineke, P. Wing, H. Wong, E. Y. B. Pun, K. W. Cheah, T. Zentgraf, and S. Zhang, “Continuous control of the nonlinearity phase for harmonic generations,” Nat. Mater. 14(6), 607–612 (2015). [CrossRef]
22. G. Li, G. Sartorello, S. Chen, L. H. Nicholls, K. F. Li, T. Zentgraf, S. Zhang, and A. V. Zayats, “Spin and geometric phase control four-wave mixing from metasurfaces,” Laser Photonics Rev. 12(6), 1800034 (2018). [CrossRef]
23. R. C. Devlin, A. Ambrosio, D. Wintz, S. L. Oscurato, A. Y. Zhu, M. Khorasaninejad, J. Oh, P. Maddalena, and F. Capasso, “Spin-to-orbital angular momentum conversion in dielectric metasurfaces,” Opt. Express 25(1), 377 (2017). [CrossRef]
24. H. Liang, Q. Lin, X. Xie, Q. Sun, Y. Wang, L. Zhou, L. Liu, X. Yu, J. Zhou, T. F. Krauss, and J. Li, “Ultrahigh numerical aperture metalens at visible wavelengths,” Nano Lett. 18(7), 4460–4466 (2018). [CrossRef]
25. A. Karnieli, S. Trajtenberg-Mills, G. D. Domenico, and A. Arie, “Experimental observation of the geometric phase in nonlinear frequency conversion,” Optica 6(11), 1401 (2019). [CrossRef]
26. E. Bahar, X. Ding, A. Dahan, H. Suchowski, and J. Moses, “Adiabatic four-wave mixing frequency conversion,” Opt. Express 26(20), 25582 (2018). [CrossRef]
27. X. Ding, D. Heberle, K. Harrington, N. Flemens, W. Chang, T. A. Birks, and J. Moses, “Observation of rapid adiabatic passage in optical four-wave mixing,” Phys. Rev. Lett. 124(15), 153902 (2020). [CrossRef]
28. G. Agrawal, “Nonlinear fiber optics,” in Optics and Photonics (Academic Press, 2012).
29. G. G. Luther, M. S. Alber, J. E. Marsden, and J. M. Robbins, “Geometric analysis of optical frequency conversion and its control in quadratic nonlinear media,” J. Opt. Soc. Am. B 17(6), 932 (2000). [CrossRef]
30. G. Porat and A. Arie, “Efficient, broadband, and robust frequency conversion by fully nonlinear adiabatic three-wave mixing,” J. Opt. Soc. Am. B 30(5), 1342 (2013). [CrossRef]
31. C. R. Phillips, C. Langrock, D. Chang, Y. W. Lin, L. Gallmann, and M. M. Fejer, “Apodization of chirped quasi-phasematching devices,” J. Opt. Soc. Am. B 30(6), 1551 (2013). [CrossRef]
32. G. Porat, Y. Silberberg, A. Arie, and H. Suchowski, “Two photon frequency conversion,” Opt. Express 20(4), 3613 (2012). [CrossRef]