MODE CONVERTER AND QUADRANT PHOTODIODE FOR SENSING OPTICAL CAVITY MODE MISMATCH

Abstract

A new technique for sensing optical cavity mode mismatch using a mode converter formed from a cylindrical lens mode converting telescope, radio frequency quadrant photodiodes (RFQPDs), and a heterodyne detection scheme. The telescope allows the conversion of the Laguerre-Gauss basis to the Hermite-Gauss (HG) basis, which can be measured with quadrant photodiodes. Conversion to the HG basis is performed optically, measurement of mode mismatched signals is performed with the RFQPDs, and a feedback error signal is obtained with heterodyne detection.

Claims

1. A system for sensing optical cavity mode mismatch, comprising: a mode converter having a pair of lenses with a focal length, wherein the pair of lenses are spaced apart by the focal length multiplied by the square root of two; and a photodiode having quadrants aligned with the mode converter.

2. The system of claim 1, wherein the focal length of the pair of lenses is a function of a waist size of an incoming beam.

3. The system of claim 2, further comprising a Gouy phase telescope positioned between the mode converter and the incoming beam.

4. The system of claim 3, wherein the incoming beam comprises an optical cavity reflection containing mismatched signals.

5. The system of claim 4, further comprising a mixer coupled to the photodiode for demodulating a signal output from the photodiode.

6. The system of claim 5, further comprising a set of feedback electronics coupled to an output of the mixer and controlling a mode matching actuator positioned between an electro-optic modulator and an optical cavity that will produce the optical cavity reflection.

7. The system of claim 6, further comprising an oscillator coupled between the electro-optic modulator and the mixer.

8. The system of claim 7, wherein the feedback electronics is configured to compare the quadrants of the photodiode to produce a mode mismatching error signal.

9. The system of claim 8, wherein the focal length of the pair of lenses is determined by f(w.sub.0)=(w.sub.0.sup.2)/(1+1/2).

10. A method of sensing optical cavity mode mismatch, comprising the steps of: passing an incoming beam through a mode converter having a pair of lenses with a focal length, wherein the pair of lenses are spaced apart by the focal length multiplied by the square root of two; and detecting the incoming beam after the mode converter with a photodiode having quadrants.

11. The method of claim 10, wherein the focal length of the pair of lenses is a function of a waist size of an incoming beam.

12. The method of claim 11, further comprising the step of passing the incoming beam through a Gouy phase telescope before the mode converter.

13. The method of claim 12, wherein the incoming beam comprises an optical cavity reflection containing mismatched signals.

14. The method of claim 13, further comprising the step of using a mixer coupled to the photodiode for demodulating the signal output from the photodiode.

15. The method of claim 14, further comprising the step of using a set of feedback electronics coupled to an output of the mixer to control a mode matching actuator positioned between an electro-optic modulator and an optical cavity that will produce the optical cavity reflection.

16. The method of claim 15, further comprising an oscillator coupled between the electro-optic modulator and the mixer.

17. The method of claim 16, wherein the feedback electronics compares the quadrants of the photodiode to produce a mode mismatching error signal.

18. The method of claim 17, wherein the focal length of the pair of lenses is determined by f(w.sub.0)=(w.sub.0.sup.2)/(1+1/2).

Description

BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWING(S)

[0006] The present invention will be more fully understood and appreciated by reading the following Detailed Description in conjunction with the accompanying drawings, in which:

[0007] FIG. 1 is a perspective view of a mode mismatch signal conversion unit according to the present invention;

[0008] FIG. 2 is a schematic of a mode mismatch signal conversion unit according to the present invention;

[0009] FIG. 3 is a series of graphs of the passing of the LG.sub.10 mode through the mode converting telescope using Fourier transform functions in MATLAB according to the present invention;

[0010] FIG. 4 is a schematic of a quadrant photodiode pattern that takes advantage of the intensity pattern of the mode converter of the present invention;

[0011] FIG. 5 is a series of graphs of the error signals extracted from the mode converter using the program FINESSE of the present invention;

[0012] FIG. 6 is a series of graphs of the optical power of higher order misalignment modes of the mode converter of the present invention;

[0013] FIG. 7 is a series of graphs of the conversion of the rotated HG.sub.01 mode to an LG mode using MATLAB according to the present invention;

[0014] FIG. 8 is a schematic of the sensing and correcting mode mismatch of a laser to an optical cavity according to the present invention;

[0015] FIG. 9 is a series of MATLAB images demonstrating conversion of the LG.sub.10 mode to the HG.sub.11 mode according to the present invention; and

[0016] FIG. 10 is a series of MATLAB images and scripts generating mode mismatch and applying Fourier transforming telescope functions according to the present invention;

[0017] FIG. 11 is a schematic of testing of a mode mismatch signal conversion unit according to the present invention;

[0018] FIG. 12 is a series of graphs illustrating the results of testing using the arrangement of FIG. 11.

DETAILED DESCRIPTION OF THE INVENTION

[0019] Referring to the figures, wherein like numerals refer to like parts throughout, there is seen in FIGS. 1 and 2, a mode converter 10 for converting a mode mismatch signal carrying LG.sub.10 mode to a HG.sub.11 mode and sensed with a quadrant photodiode 16. The LG.sub.10 passes through two cylindrical lenses 12 and 14 spaced by their focal length multiplied by 2. It is critical that the focal length be chosen as a function of the incoming beam waist size f(w.sub.o)=(w.sub.o.sup.2)/(1+1/2).

[0020] Referring to FIG. 8, a laser beam 20 is mode mismatched to an optical cavity 22. Sidebands are created by the Electro-Optic Modulator (EOM) 24, which are used for the heterodyne detection. The optical cavity reflection contains mismatch signals and, using a beam splitter 26 and mirror 28, is guided through a Gouy phase telescope 30 (for proper waist location placement) and finally through a mode converter 10 such as that seen in FIG. 1 or FIG. 2. The radio frequency quadrant photodiode signal from a quadrant photodiode 16 is demodulated by a mixer 36 and passed through feedback electronics 38 that control the mode matching actuator 40 (depicted as a displacement actuator, but could also be thermally deformable). The feedback electronics 38 also compare the four quadrants to extract the mode mismatching error signal.

[0021] In order to implement a mode converter into an actively controlled optical system, it is necessary to derive an error signal from the output that is linearly proportional to the waist position or size of the optical cavity. This is done by applying a mask to the output images as seen in FIG. 3, FIG. 4 and FIG. 5 which isolates the power on each section of a quadrant photodiode. MATLAB was used to pass the LG.sub.10 mode through the mode converting telescope function. The output result is a HG.sub.11 mode rotated 45 degrees in the transverse direction. FIG. 4 shows the photodiode pattern which takes advantage of the intensity pattern provided by mode converter. FIG. 5 shows a comparison of the error signals derived by the mode converter QPD vs. the bullseye photodiode and the error signals extracted from the mode converter show a linear response to shifts in the waist location.

[0022] Typical optical cavity alignment sensing requires the ability to measure the HG.sub.01 and HG.sub.10 modes with a quadrant photodiode. Mode mismatch manifests itself as concentrically symmetric LG modes. Since the whole beam will be affected by the mode converter, it is worthwhile to examine what happens to the alignment signals after they pass through. Passing a well aligned LG mode through the mode converter requires no specific rotation angle since the LG mode is radially symmetric. The cylindrical lenses will always produce a HG mode that is rotated 45 degrees from the cylindrical lens focusing axis in FIG. 3. As seen in FIG. 5, optical power of higher order misalignment modes are not affected by the mode converting telescope. However, note that there is an electric field sign flip.

[0023] Passing HG.sub.01 or HG.sub.10 alignment signals through the mode converter does require a specific rotation angle. The alignment signals must enter at 0 degrees with respect to the focusing axis to maintain their shape. If this condition is not respected then mixing of the modes will occur, as seen in FIG. 7.

[0024] The mode converter does not affect the misalignment signals in a notable way as long as the alignment modes enter at 0 degrees with respect to the focusing axis.

[0025] FIG. 6 shows how alignment modes perpendicular to the lens axis are completely unchanged, while alignment signals parallel to the lens axis have their electric field flipped. Since power is amplitude squared, the photodiodes will not measure anything different and hence the alignment signals pass through the mode converter are just as useful as before.

EXAMPLE

[0026] Computer simulation provided a quick method for testing our prediction before performing the experiment. A combination of FINESSE, which mainly uses ABCD matrix math at its core, and MATLAB was selected for the experiment. FINESSE was chosen because it had previously been used to generate mode mismatching signals from an optical cavity that were sensed with bullseye photodiodes and thus provides a good baseline for simulation comparison. MATLAB was used as a means to process our transverse electric field expression using Fourier Optics. FINESSE was used to produce the mode mismatched cavity in FIG. 9. As seen in FIGS. 9 and 10, an exemplary MATLAB study shows that the present invention was able to convert the LG.sub.10 mode to the HG.sub.11 mode. Fourier optics were used to clearly show the conversion. As seen in FIG. 9, MATLAB was used to pass the LG.sub.10 mode through the mode converting telescope. The output result is clearly a HG.sub.11 mode rotated 45 degrees in the transverse direction. Referring to FIG. 10, a MATLAB script was used to plot higher order HG modes and subtract them pixel by pixel. The |HG.sub.11| mode is plotted directly in the top left of FIG. 10. In the bottom center of FIG. 10, the sum of |HG.sub.10| and |HG.sub.01| may be seen. On the top right, |HG.sub.20||HG.sub.02| may be seen. The peaks are much closer together in the HG.sub.11 mode. However, the beam composition looks similar enough to HG.sub.11 that this could still be detected on a quadrant photodiode. Thus, it is possible to measure mode mismatching signals using well developed alignment sensors.

APPENDIX

[0027] The complex beam parameter of a Gaussian beam with Rayleigh range z.sub.R is defined as:

q=+i.sub.R.(1)

[0028] Beam size w and phase front radius of curvature R are then given by

[00001]$\begin{array}{cc}\frac{1}{q}=\frac{1}{R}-i.Math.\frac{}{.Math..Math.w^2},& \left(2\right)\end{array}$

where =2/k is the wave length of the light. It allows expressing the Gaussian beam in a simple form:

[00002]$\begin{array}{cc}\left(x,y,q\right)=A\left(x,y,q\right).Math.e^{-\mathrm{ikz}}& \left(3\right)\\ A\left(x,y,q\right)=\frac{A}{q}.Math.e^{-\mathrm{ik}.Math.\frac{z^2+y^2}{2.Math.q}}& \left(4\right)\end{array}$

where A is a complex constant (amplitude). It can be helpful to introduce the field amplitude on the optical axis, =A/q, which now evolves along the z-axis due to the Gouy phase evolution. Thus, for any given location on the optical axis z, the Gaussian beam is completely described by the two complex parameters and q. The main advantage of this formalism becomes apparent when using ray-transfer matrices M defined in geometric optics (e.g. Saleh, Teich) to represent the action of a full optical system. The two complex parameters after the system (q.sub.f, .sub.f), are given in terms of the initial parameters (q.sub.i, .sub.i) by

[00003]$\begin{array}{cc}M\left(\begin{array}{c}\frac{1}{{}_i}\\ \frac{1}{{}_i.Math.q_i}\end{array}\right)=\left(\begin{array}{c}\frac{1}{{}_f}\\ \frac{1}{{}_f.Math.q_f}\end{array}\right),& \left(5\right)\end{array}$

and the change of the Gouy phase through the system, , is given by

[00004]$\begin{array}{cc}{e}^{i.Math..Math..Math..Math.}=\sqrt{\frac{{}_f}{{}_f^{*}}.Math.\frac{{}_i^{*}}{{}_i}}.& \left(6\right)\end{array}$

[0029] This expression is consistent with the usual definition of Gouy phase for a Gaussian beam as =arctan z/z.sub.R, and can be proven by verifying it for a pure free-space propagation and a pure lens.

[0030] If astigmatism is introduced, either intentionally with cylindrical lenses or accidentally through imperfections, cylindrical symmetry around the beam axis will be lost. As long as we introduce this astigmatism along a pre-determined axis (say the x-axis), we can simply proceed by introducing separate q-parameters for the x- and y-axis, q.sub.x and q.sub.y. Since ray-transfer matrices are introduced with only 1 transverse axis, the propagation of q.sub.x and q.sub.y is done with ray-transfer matrices defined for the corresponding transverse axis. Thus we now have a separately-defined beam size w.sub.x, w.sub.y, phase front radius of curvature R.sub.x, R.sub.y, Rayleigh range z.sub.Rx, z.sub.Ry and Gouy phase x and y for each of the two transverse directions. The corresponding fundamental Gaussian beam is given by

[00005]$\begin{array}{cc}\left(x,y,q_x,q_y\right)=A\left(x,y,q_x,q_y\right).Math.e^{-\mathrm{ikz}}& \left(7\right)\\ A\left(x,y,q_x,q_y\right)=\frac{A}{\sqrt{q_x.Math.q_y}}.Math.e^{-\mathrm{ik}.Math.\frac{x^2}{2.Math.\mathrm{qx}}}.Math.e^{-\mathrm{ik}.Math.\frac{y^2}{2.Math.q_y}}& \left(8\right)\end{array}$

[0031] where A is again a complex amplitude. Next we introduce the Hermite-Gaussian basis set corresponding to the fundamental Gaussian beam. In the literature this is typically done only relative to a single q-parameter, but directly generalizes to the case with separate q.sub.x and q.sub.y parameters:

[00006]$\begin{array}{cc}{}_{n.Math..Math.m}\left(x,y,q_x,q_y\right)=A_{n.Math..Math.m}\left(x,y,q_x,q_y\right).Math.e^{-\mathrm{ikz}}& \left(9\right)\\ A_{n.Math..Math.m}\left(x,y,q_x,q_y\right)={\mathrm{NA}}_n\left(x,q_x\right).Math.A_m\left(y,q_y\right)& \left(10\right)\\ A_p\left(,q_{}\right)=e^{i.Math..Math.p.Math..Math.{}_{}}.Math.\sqrt{\frac{1}{2^p.Math.p!}.Math.{}_{}}.Math.H_p\left(\sqrt{2}.Math.\frac{}{w_{}}\right).Math.e^{-\mathrm{ik}.Math.\frac{{}^2}{2.Math.q_{}}}& \left(11\right)\\ {}_{}=\sqrt{\frac{2}{}}.Math.\frac{e^{i.Math..Math.{}_{}}}{w_{}}=\sqrt{\frac{2.Math.z_R}{}}.Math.\frac{i}{q_{}}& \left(12\right)\\ H_0\left(\right)=1,H_{p+1}\left(\right)=2.Math..Math..Math.H_p\left(\right)-\frac{d}{d.Math..Math.}.Math.H_p\left(\right)& \left(13\right)\end{array}$

[0032] Here, we redefined the overall amplitude N such that the total power P in a mode is simply given by P=|.sub.nm|.sup.2dxdy=|N|.sup.2. That equation (8) is of the same form can be seen using the identity iz.sub.R/q=e.sup.iw.sub.0/w. Furthermore, .sub. was defined in analog to the field amplitude introduced after equation (4), that is the field amplitude on the optical axis of the fundamental mode. It thus evolves, together with q.sub., according to equations (5) and (6). Note though that there is an extra Gouy phase term for the higher order modes that is explicitly excluded from the definition of .sub.. As a result, the overall Gouy phase evolution of .sub.nm(x, y, q.sub.x, q.sub.y) is proportional to e.sup.i(n+1/2).sup.z.sup.+i(m+1/2).sup.y.

[0033] As expected, these modes still solves the paraxial Helmholtz equation

[00007]$\begin{array}{cc}\left({}_T-2.Math.\mathrm{ik}.Math.\frac{}{z}\right).Math.A_{n.Math..Math.m}\left(x,y,q_x,q_y\right)=0& \left(14\right)\end{array}$

exactly.

Example 2

[0034] Referring to FIG. 11, a typical heterodyne detection approach was tested using BPDs compared to the proposed QPD plus a mode converter. The error signal produced by either the BPDs or QPDs is then fed back to thermal lens actuators which can adjust both alignment and mode matching. The top right and left plots of FIG. 12 were produced using a BPD sensor in FINESSE. The bottom right and left plots of FIG. 12 were produced using MATLAB. The MATLAB script passes beams through a Fourier optic representation of the mode converter. After conversion, a reference 00 beam is beat against a higher order mode. The result is then segmented into quadrants and combined A+DBC to produce the error signals plotted.

[0035] For the present invention, some design considerations are worth noting. First, after the mode converter the signals for pitch and yaw will be in different Gouy phases. For ideal sensing of all 6 degrees of freedom, 3 QPDs are needed at 0, 45 and 90 Gouy phase separation. However, in practice a setup with 2 QPDs 67.5 apart may be sufficient. The beam should be focused so that the beam waist is in the middle of the mode converter. The cylindrical lens separation must be f{square root over (2)}, where f is the focal length of each of the cylindrical lenses. This constrains the Rayleigh range of the reference beam and therefore its beam size w.sub.0={square root over ((1+1/{square root over (2)})f/)}.

[0036] The benefits of mode converter with quadrant photodiodes according to the present invention include the ability to measure mode mismatch using existing QPD. The mode converter also preserves cavity misalignment signals. For the present invention, no new electronics or channels are needed and QPDs are easier to align than BPDs. Compared to BPDs, QPDs are off-the-shelf and have a much better matched quadrant capacities on optical gains. Only one mode converter, followed by normal Gouy phase telescopes for each sensor is needed to sense all four alignment and two mode-match degrees of freedom. The only optical component critically sensitive to beam size is the mode converter itself. It is much easier to change than the segment size of (multiple) BPD.

[0037] As described above, the present invention may be a system, a method, and/or a computer program associated therewith and is described herein with reference to flowcharts and block diagrams of methods and systems. The flowchart and block diagrams illustrate the architecture, functionality, and operation of possible implementations of systems, methods, and computer programs of the present invention. It should be understood that each block of the flowcharts and block diagrams can be implemented by computer readable program instructions in software, firmware, or dedicated analog or digital circuits. These computer readable program instructions may be implemented on the processor of a general purpose computer, a special purpose computer, or other programmable data processing apparatus to produce a machine that implements a part or all of any of the blocks in the flowcharts and block diagrams. Each block in the flowchart or block diagrams may represent a module, segment, or portion of instructions, which comprises one or more executable instructions for implementing the specified logical functions. It should also be noted that each block of the block diagrams and flowchart illustrations, or combinations of blocks in the block diagrams and flowcharts, can be implemented by special purpose hardware-based systems that perform the specified functions or acts or carry out combinations of special purpose hardware and computer instructions.