1. Introduction

Optical beam modes are transverse eigensolutions of the paraxial wave equation (PWE) of electrodynamics [1]. The modes hold for beams that are well collimated such that the curvature in the propagation direction can be neglected. Three families of solutions have been found by separation of variables and take on the geometry of their respective coordinate system such as Cartesian, cylindrical and elliptical polar coordinates [1, 2 ]. The fundamental mode in all three families is the commonly encountered Gaussian beam produced by many laser systems. Of these solutions, optical vortices have attracted considerable attention from the scientific community because they possess an external degree of freedom known as optical orbital angular momentum (OAM) in addition to the internal spin angular momentum [3]. Applications of optical modes include: optical manipulation of Bose Einstein (BECs) condensates, such as tight confinement using Hermite Gaussian (HG) modes [4]; transfer of angular momentum with Laguerre Gaussian (LG) modes [5, 6 ]; optical spanners [7] and pumps [8]; imaging of cytoskeletal filaments (stimulated emission depletion) [9] multiplexing and cryptography [10, 11 ] and the search for exoplanets (optical vortex coronagraph) [12]. The success of many high-field experiments is possible by advances in the generation and analysis of beam modes in femtosecond fields [2, 13, 14 ]. However, investigations of highly nonlinear interactions of beam modes with matter have remained largely unexplored.

Recently a number of investigations involving the nonlinear interaction of optical vortices with matter appeared in the literature [15—17 ]. In experiments presented in [15], the authors crossed time-delayed chirped pulses in a lead tungsten PbWO₄ crystal and demonstrated for the first time the coherent transfer of optical angular momentum in the generation of cascaded Raman sidebands. In these experiments, topological charge measurements verified a predicted orbital angular momentum algebra (OAM-algebra). This OAM-algebra was later verified in experiments given in [16, 17 ]. In [16] the authors employed two optical parametric amplifiers to produce two spectrally separated beams that were crossed in a Raman active crystal, and in [17] the authors utilized a dichroic beam splitter and spectral filters. In each of these previous experiments, a spiral phase plate was used producing only the ${\text{LG}}_{\rho =0,\ell =1}$ transverse mode having $\ell =\pm 1$. This therefore limited the scope of the previous experiments.

In the present work, we extend upon our earlier work [15] by investigating the production of Raman sidebands with higher-order helical beams such as the Laguerre Gaussian (LG) and Ince Gaussian (IG) modes; the non-helical even Laguerre Gaussian (${\text{LG}}^e$) and odd Laguerre Gaussian (${\text{LG}}^o$), and Hermite Gaussian modes (HG). From our experiments, we verify that the OAM-algebra holds for the high order helical-LG beams as well as for the helical-IG modes. In the case of the non-helical modes, such as the even and odd LG modes or the HG modes, the OAM-algebra derived in [15] cannot be used; therefore, a new OAM-algebra is developed that describes the quantum coherent transfer of phase and amplitude information into the sideband orders.

Throughout this work all optical beam modes were produced using computer generated holograms (CGH) displayed on a liquid crystal spatial light modulator (SLM). The CGHs were produced by interfering a plane wave with the paraxial solutions of the desired beam mode. These solutions, which correspond to different coordinate geometries are:

2. Experimental setup

In our experiments, amplified radiation from a Ti:sapphire laser system was produced by chirped-pulse amplification (Coherent Legend) producing ~50 fs laser pulses at a repetition rate of 1 kHz. The central wavelength was 800 nm, and the energy per pulse was ~1 mJ. Output radiation from the laser was sent through a variable iris to control both mode quality [16] and beam power in the setup. Following the iris, the radiation was sent through a four-port Michelson interferometer Fig. 1(a) . The purpose of the four-port interferometer is to produce a beam encoded with the phase information of the desired optical mode and to produce a reference beam needed to investigate the phase content of the generated sidebands. In contrast to our previous setup [15], one of the mirrors in the stationary arm of the four-port was replaced by a parallel aligned liquid crystal on silicon spatial light modulator (Hamamatsu LCOS-SLM 104683). This device was addressed by a control unit and computer to display the CGHs. The inset in Fig. 1(a) shows two examples of such holograms. A routing mirror was used to steer radiation from the output of the four-port Michelson into the beam crossing setup shown Fig. 1(b). A lens (L) having a nominal focal length of 40 cm was used to focus the pump and Stokes beams into the Raman-active crystal. The separation distance between the beams at the position of the lens was ~2.25 cm and resulted in a full-angle of 3.2 degrees between the two beams. The nonlinear crystal was a 0.5 mm thick lead tungsten crystal (MTI Corporation). Images of the pump and Stokes beams, and the generated Stokes and anti-Stokes orders were captured using a charged-coupled device (Spiricon, SP503U) having a resolution of $640\times 480$pixels. Immediately before the crystal, the pump and Stokes beams had an energy per pulse of $~3\text{mJ}$. The pulse duration was extrapolated to be $~400fs$ using a Grenouille (UPM-8-20, 18fs—180fs). With a spot size of 2mm (CCD) at the lens and for the given parameters, the peak intensity was estimated to be about $~3\times {10}^{10}{\text{W/cm}}^{\text{2}}$, which assumes a diffraction-limited beam. For larger beam modes or for higher-order sidebands, the power needed to be increased to at most ~4mJ.

 

Fig. 1 Illustration of (a) the four-port Michelson and (b) the beam crossing setup with the following elements: beam splitter (BS), Mirror (M), Lens (L), translation stage (T), spatial light modulator (SLM) and lead tungsten crystal (PbWO4). The SLM, used to produce the desired beam mode, is positioned in the stationary arm of the four-port Michelson. Radiation from the variable arm was used as a reference.

Download Full Size | PDF

To produce Raman sidebands, we employed time-delayed chirped pulses in the beam crossing setup. Figure 2(a) illustrates how two delayed chirped-pulses in our experiments result in a frequency difference capable of being tuned to the Raman transitions in a nonlinear medium. For linear chirped pulses, the phase content has a quadratic time dependence, so that the instantaneous frequency is given by $\omega ={\omega }_0+bt$. The two up-chirped pulses shown in Fig. 2(a) are delayed by $t_d$ and result in a frequency difference of $\Delta \omega =b{t}_d$ for all times $t$. In our experiments, the first folding mirrors in the compressor of the laser was slightly detuned by about ~1mm. The sidebands were then found by adjusting the time delay between the two pulses using the translation stage shown in Fig. 2(b). In this way, the detuning $\Delta \omega =b{t}_d$ was found experimentally.

 

Fig. 2 Cascaded generation of sidebands: (a). Schematic representation of the time-delayed chirped pulse scheme employed to access the Raman transitions of the Raman-active crystal (PbWO4). Shown are the spectral density distribution versus time plots of two positively $b>0$, linearly chirped pulses delayed in time by $t_d$. At all times, the frequency difference between the pulses is $\Delta \omega =b{t}_d$. In order to access the Raman transitions, this difference must equal the Raman frequency ${\omega }_R=\omega t_d$. (b). Diagram of cascade generation of sideband frequencies.

Download Full Size | PDF

For completeness, the OAM-algebra derived by us in [15] for the helical LG beam of order $\left|\ell \right|=1$ is briefly presented. In all further discussions, the different Raman orders will be signified by SN for Stokes and ASN for anti-Stokes, where N denotes the order. The pump is denoted by P and the Stokes by S. In Fig. 2(b) the level diagram shows two real states and a number of virtual states. To find the frequency of the radiation in the AS1 order, the nonlinear polarization can be written as the product of plane wave solutions$e^{i{\omega }_{P}t}{e}^{-i{\omega }_{S}t}{e}^{i{\omega }_{P}t}$. The resulting electric field in the first anti-Stokes order then has a frequency of${\omega }_{AS1}=2{\omega }_P-{\omega }_S$. Multiplying through by $\hslash$ motivates a photon description in which the first anti-Stokes order can be described by the absorption of a pump photon followed by the stimulated emission of a Stokes photon with the subsequent absorption of another pump photon. The helical LG beams carry a phase factor of $e^{i\ell \theta }$ in which its argument $\ell \theta$ comes along with the frequency, and therefore the same algebra ${\ell }_{AS1}=2{\ell }_P-{\ell }_S$ follows. The amount of OAM in the AS2 order can be determined by the absorption of a pump photon followed by emission of a Stokes photon with the subsequent absorption of an AS1 photon, so that ${\ell }_{AS2}={\ell }_P-{\ell }_S+{\ell }_{AS1}=3{\ell }_P-2{\ell }_S$. The absorption of the AS1 photon in the generation of the AS2 radiation demonstrates why this is a cascade process.

Continuing with this photon description, the OAM-algebra for anti-Stokes orders is found to be${\ell }_n^{AS}=(n+1){\ell }_P-n{\ell }_S$. Using conservation of momentum ${\ell }_P+{\ell }_S={\ell }_n^{AS}+{\ell }_n^S$ the algebra for the Stokes orders is found to be ${\ell }_n^S=(n+1){\ell }_S-n{\ell }_P$, which can be deduced by reversing the direction of the pump and Stokes arrows in Fig. 2(b). The image [bottom inset in Fig. 2(b)] shows measured sidebands up to the 19th anti-Stokes order. Due to spatially-dependent nonlinearities and dispersion, the orders do not lie along a line. The distortion can be manipulated by varying the overlap between the two beams.

3. Sideband generation with helical LG and IG beams

To check the functionality of our modified setup and to compare with results from our previous setup, on-axis computer generated holograms of helical Laguerre Gaussian beams [insets in Fig. 1(a)] were displayed on the SLM positioned in the four-port Michelson. The produced optical beam modes were then sent into the beam crossing setup and allowed to interact in the Raman-active crystal. When the device is used to produce off-axis holograms, amplitude information can also be encoded in the computer generated holograms by redirecting radiation between orders [2, 14 ]. This phase-amplitude encoding produces high fidelity beam modes but at the expense of pulse energy. In our setup, reduction of pulse energy from the three beam splitters, the mode-improving iris, chirp, and the reduction due to the efficiency of operating the SLM in off-axis mode resulted in a focal intensity insufficient to reliably produce Raman sidebands. For this reason, the SLM was operated in on-axis mode, which converts a large proportion of its input radiation into the desired mode.

In our previous setup, when radiation with $\ell =\pm 1$ was passed through the setup with mirror M₅ in place, the topological charges in the pump and Stokes beams were found to have the same magnitude but were opposite in sign. In contrast, when mirror M₅ was removed from the setup, the topological charges had the same magnitude and sign. This property of the beam carries over to higher-order modes with $\ell >1$ and can be understood by a direct comparison with circularly polarized light. To show this, we decompose the Laguerre Gaussian beams with angular mode numbers $\ell$ and radial mode number $\rho =0$ in the Hermite Gaussian basis,

Figure 3 shows the beam profiles and interferograms of the generated sidebands using higher-order LG beams ($\ell =1,2,3$) with mirror M₅ in place (balanced arms). For balanced arms, the mode number of the pump and Stokes are equal ${\ell }_P={\ell }_S=\ell$. The OAM-algebra for this configuration reduces to ${\ell }_n^{AS}={\ell }_n^S=\ell$. The topological charges in the pump and Stokes beams can be verified from the interferograms in Rows 2, 4 columns 1—5, and row 2 columns 6—9. In these images it can be seen that the multifurcations and bifurcations of all orders are facing in the same direction, and this indicates that ${\ell }_P={\ell }_S=\ell$ verifying that the OAM-algebra is indeed satisfied.

 

Fig. 3 Generation of Raman sidebands with LG beams having the same helicity ${\ell }_P={\ell }_S$. Columns 1—5: S1, S, P, AS1 and AS2; rows 1 and 3 are sidebands generated with$\ell =1\text{and}\text{2}$respectively. Rows 2 and 4 are the interferograms of those beams in rows 1 and 3 with the reference beam from the 4-port Michelson. The measured topological charges are consistent with that predicted by the OAM-algebra. Columns 6—9: S1, S, P, and AS1; row 1 are sidebands produced with $\ell =3$ and row 2 shows the interferograms. Panels (a) and (b) show the intensity distribution and the interferogram of the AS1 order beam produced with $\ell =5$. As can be seen, the multifurcation has broken up into 5 single bifurcations [denoted by the red dots in panel (b)].

Download Full Size | PDF

With mirror M₅ removed from the setup, the arms are unbalanced and the relative signs of $\ell$ between the pump and Stokes and the Stokes and anti-Stokes orders are opposite ${\ell }_n^{AS}=-{\ell }_n^S=\ell$, so that the OAM-algebra reduces to ${\ell }_n^{AS}=(2n+1)\ell$ and ${\ell }_n^S=-(2n+1)\ell$. Figure 4 shows sidebands produced in this manner. In this case, interferograms of the pump and Stokes beams show multifurcation facing each other showing that indeed ${\ell }_n^{AS}=-{\ell }_n^S=\ell$. The OAM-algebra for each order was verified by fringe counting. In our experiments, it was noted that the beam profiles in many of the sidebands were distorted, and the interferograms showed that the multifurcations in the higher-order modes separated into single bifurcations. This instability has been previously investigated and is found to commonly occur for high-order LG modes [19]. Despite the distortions, verification of the OAM in the generated sidebands is possible by identifying and summing the number of bifurcations or counting the fringes around a cluster of multifurcations. All data was found to be consistent with the OAM-algebra derived in [15]. We found that the quality of the generated sidebands is sensitive to many experimental parameters. Such parameters include intensity, input beam size and chirp. When producing high order modes larger beam sizes result, and consequently the peak intensity drops. Further increase in intensity also results in distortion of the generated sidebands. Significant distortion and the presence of radiation from four-wave mixing were difficult to mitigate.

 

Fig. 4 Generation of Raman sidebands with LG beam of the same order but opposite helicity ${\ell }_P=-{\ell }_S$(unbalanced arms). Columns 1—5: S1, S, P, AS1 and AS2; rows 1 and 2, and rows 3 and 4 show sidebands generated with$\ell =1$ and $\ell =2$ respectively. Rows 2 and 4 are the interference of those beams in rows 1 and 3 with the reference beam for the Michelson. The measured topological charges are consistent with those predicted by the OAM-algebra ${\ell }_n^{AS}=(n+1){\ell }_P-n{\ell }_S$. Column (a) shows images of the AS1 beams generated with $\ell =$4, 5, 6 and 7 (rows 1—d respectively); and column (b) shows the interferograms. In all data, the topological charges were determined by counting multifurcation and fringes.

Download Full Size | PDF

A natural extension of the OAM-algebra can be made to the OAM-carrying helical Ince Gaussian beams Eq. (2). Like the helical-LG beams, the IG beams possess OAM. When the ellipticity parameter Eq. (1c) of the helical-IG beams is zero $\epsilon =0$, the helical-IG beams are equivalent to the helical-LG beams. This equivalence is a consequence of the relationship between cylindrical elliptic coordinates and cylindrical polar coordinates. When the ellipticity parameter is greater than zero$\epsilon >0$, the multifurcation of the beam splits into a string of bifurcations along a line joining the two foci of the coordinate. The total topological charge of the beam is then just equal to the sum of the bifurcations. For this reason, helical-IG beams of order $\ell >1$ and $\epsilon \ge 0$ were investigated.

Figure 5 shows generated sidebands of IG beams with $\ell =$2, 3, 4 and 5 with varying non-zero ellipticities. As with the helical-LG beams, we observed distortion of the generated sidebands, but by inspection of the interferograms, the transfer of OAM into the sidebands was found to follow the same OAM-algebra as that for the helical-LG beams. For the current experimental conditions, an interesting effect was observed when the ellipticity parameter was varied between $\epsilon =0.1-0.4$. When $\epsilon =0$, all sidebands demonstrated a distortion in such a way that the bifurcations were separated about a line at ~70 degrees from the horizontal. This is similar to the instability of the higher-order LG beams. By varying$\epsilon$, the beams were observed to take on a more symmetric shape (columns 2 and 7 in Fig. 5) and the interferograms showed a tendency for the individual bifurcations to group more towards the center of the beams. This effect may present a method to actively correct for distortions in high-order OAM-containing beams in the focus or when delivered to a target.

 

Fig. 5 Ellipticity-dependence of the generation Raman sidebands with Ince-Gaussian beams of the same order and helicity ${\ell }_P={\ell }_S$(balanced arms). All images are of the AS1 order. Columns 1—5: $\epsilon =(0.0,0.2,0.4,0.6,1.0)$ [columns 6—10: $\epsilon =(0.0,0.1,0.2,0.3,1.0)$] rows 1 and 3 (columns 1—5) are sidebands generated with $\ell =2,\text{and}\text{3}$. Rows 1 and 3, columns 6—10 are sidebands generated with $\ell =4,\text{and}\text{5}$. The interferograms appear under the images of the orders. The measured topological charges are determined by counting multifurcation and fringes and are consistent with those predicted by the OAM-algebra found for the LG beams. The general trend is for the multifurcations (bifurcations) to split up and move along a line as the ellipticity increases. An interesting consequence of the control of the bifurcation is that the initial instability of the higher order IG beams can be partially corrected for by changing the ellipticity parameter (beams in columns 2 and 7 have a more circular shape than in column 1 and 6).

Download Full Size | PDF

4. Even Laguerre Gaussian modes

As shown above and in our previous work [15], a straightforward derivation based on parametric four-wave mixing allowed us to determine the orbital angular momentum algebra of the generated sidebands for the helical-LG beams. We now turn our attention to the even and odd Laguerre Gaussian beams of Eq. (1) . The solutions of these beams are similar to those of the helical-LG beams with the exception that the total OAM sums to zero. Unlike the OAM-algebra for the helical-LG beams, the derivations for the even and odd beams are more involved. Our considerations in this section will be for the even LG beams; however the results directly carry over to the odd LG beams by adding $\pi /2$ to the arguments of the cosine functions. To gain understanding of how OAM from the even LG beams affect the generation process, we derive the OAM-algebra with the help of four-wave mixing relations and by using the energy level diagram in Fig. 2(b) as a guide. A consequence of this scheme is the neglecting of terms in the polarizability that do not satisfy the cascaded Raman process.

In nonlinear optics, the polarizability can generally be expanded as a power series in the electric field $P(t)={\epsilon }_0{\chi }^{(1)}E(t)+{\epsilon }_0{\chi }^{(2)}{E}^2(t)+{\epsilon }_0{\chi }^{(3)}{E}^3(t)+\cdots$, where the ${\chi }^{(n)}$ are the nonlinear susceptibilities [20]. The scalar wave equation for the generated radiation field is

To obtain an analytical expression, we note that the right hand side of Eq. (8) is periodic in $\theta$ and motivates a periodic solution of the form$E_{ASN}(\theta )={\displaystyle \sum a_n\cos (n\theta )}$. With this substitution in Eq. (8) and with$B=\left[\left(N+1\right){\omega }_P-N{\omega }_S\right]/{\epsilon }_0{c}^2$, Eq. (8) can then be written as

Experimental results of the generation of Raman sidebands with the pump and Stokes beams in the modes ${\text{LG}}_{0,2}^e$ and ${\text{LG}}_{0,3}^e$ are shown in Fig. 6 . Row 1 shows the S1, S, P, AS1 and AS2 beams (columns 1—5 respectively) for the ${\text{LG}}_{0,2}^e$ modes, and row 2 shows the same orders for the ${\text{LG}}_{0,3}^e$ beam. From Eq. (10), the OAM-algebra for the S1 and AS1 orders ($N=1$) suggests that the resulting beams will have contributions from $\ell$, $-3\ell$ and $3\ell$. For the ${\text{LG}}_{0,2}^e$ modes ($\ell =2$) in row 1, the expected contributions are from ${\text{LG}}_{0,2}^e$ and ${\text{LG}}_{0,6}^e$, and for the ${\text{LG}}_{0,3}^e$ modes in row 2 the expected contributions are from ${\text{LG}}_{0,3}^e$ and ${\text{LG}}_{0,9}^e$. For the $\ell =2$ modes in the top row, the ${\text{LG}}_{0,2}^e$ contribution is the same as the pump and Stokes and has 4nodes, while the contribution from the ${\text{LG}}_{0,6}^e$ mode has 12 nodes. This result is consistent with the measured S1 and AS1 orders in row 1. For the S1 and AS1 orders in row 2, we expect contributions from ${\text{LG}}_{0,3}^e$ which are the pump and the Stokes modes and from ${\text{LG}}_{0,9}^e$ which has 18 angular nodes. These results are in excellent agreement with the measured data. For the AS2 order, Eq. (10) predicts that this order is expected to have contribution from even LG beam of mode orders ${\text{LG}}_{0,3}^e$,${\text{LG}}_{0,9}^e$ and ${\text{LG}}_{0,15}^e$, which have angular nodes of 6, 18 and 30. Again these results are consistent with the measured data and suggest that Eq. (10) is capable of describing the angular momentum content of the generated sidebands.

 

Fig. 6 Generation of Raman sidebands with ^eLG beams. Columns a—e: S1, S, P, AS1 and AS2; row 1 are sidebands generated with $\ell =2$ and Row 2 are sidebands generated with $\ell =3$. For the even (odd) LG beams, the number of angular nodes is $2\ell$. From the OAM-algebra (see text, Eq. (10)), the AS1 and S1 orders are expected to have modal contributions from $\ell$ and $3\ell$. So for ${\text{LG}}_{0,2}^e$ in the first row, the S1 and AS1 order are expected to have similarities with ${\text{LG}}_{0,2}^e$ and ${\text{LG}}_{0,6}^e$modes having 4 and 12 nodes. In the bottom row $\ell =3$ the S1 and AS1 orders are expected to have contributions from the ${\text{LG}}_{0,3}^e$ and ${\text{LG}}_{0,9}^e$modes. These modes have 6 and 18 nodes respectively. For AS2 sideband, the OAM-algebra results in angular mode numbers of $\ell$, $3\ell$ and $5\ell$. For ${\text{LG}}_{0,2}^e$(top row) the AS2 order is expected to have contributions from ${\text{LG}}_{0,2}^e$, ${\text{LG}}_{0,6}^e$ and ${\text{LG}}_{0,10}^e$ (4, 12 and 20 angular nodes) and for ${\text{LG}}_{0,2}^e$(bottom row) the AS2 order is expected to have contributions from ${\text{LG}}_{0,3}^e$, ${\text{LG}}_{0,9}^e$ and ${\text{LG}}_{0,15}^e$(6, 18 and 30 angular nodes).

Download Full Size | PDF

Because the holographic generation of the pump and Stokes beams is an imperfect process, it may happen that impurities in the generated modes produce similar modal lobes in the generated Raman orders. Although inspection of the imaging suggests that this is not a likely conclusion, nevertheless we decided to investigate the generation process further by blocking a modal lobe of the Stokes beam. Figure 7 shows the result of this experiment. Row 1 shows the S, P, and AS1 orders (from left to right) with an even ${\text{LG}}_{0,2}^e$ beam. In this case, no modal lobes were blocked. The data in row 2 shows the measured result when a single modal lobe of the Stokes beam is blocked. In this case the generated AS1 order shows a nearly identical pattern as compared to the case with no blocking. This suggests that the lowest order modal contribution predicted by the OAM algebra is being generated in the nonlinear process. The result also indicates that the entire mode structure is coherently participating in the generation process, and that the highly nonlinear nature of the process is not causing each modal lobe to act independently. Finally, the Hermite Gaussian modes can be expanded in the even and odd LG basis by a transformation. For completeness, the expressions for the expansions are given in Appendix B. For the work presented here, it can be seen that the first row in Fig. 6 and the modes in Fig. 7 are both an ${\text{LG}}_{0,2}^o$, and a ${\text{HG}}_{1,1}$ rotated at a 45 degree angle and the expansion is given in Eq. (26).

 

Fig. 7 Generation of Raman sidebands with even ^eLG beams having a blocked lobe. Columns 1—3: S, P, and AS1. Row 1 is generated with $\ell =2$. Row 2 is generation with even ^eLG beams having$\ell =2$. In the last row, a modal lobe was blocked to investigate the generation of the lowest order contribution to the OAM-algebra which in this case is ${\text{LG}}_{0,2}^e$. The AS1 orders for both scenarios (blocked and unblocked) are nearly identical demonstrating the generation of the predicted modes and that the entire modes structures of the pump and Stokes are participating in the generation process.

Download Full Size | PDF

5. Discussion and summary

In the synthesis of few-cycle pulses of radiation from the combination of Raman sidebands, it is possible to produce a short pulse in a pure transversal mode when using the helical LG and IG beams. This is because when ${\ell }_S={\ell }_P$, the OAM-algebra shows that the OAM in each order is the same. The question as to whether a similar situation can occur for the even (odd) LG and the HG beams naturally arises. The extended algebra of Eq. (10) shows that all the sidebands will have contribution from many modes, but they will all have a contribution from the generating mode. It may therefore be possible to spatially filter out the mode impurities before combining the sidebands in the synthesis of few-cycle pulse

In conclusion, we have generated Raman sidebands using various beam modes and investigated the nonlinear interaction of these beams in a Raman active crystal. We have generalized the OAM-algebra to include the OAM content of higher-order LG and IG beams, even and odd LG beams, and HG beams. Measurements of interference patterns produced with the reference beam in a simultaneous Young double slit experiment provided quantitative confirmation of the derived results.