Polarization Converter and Circuit Elements for Use in Optical Waveguides

This invention relates to optical waveguide systems and particularly to such devices modifying the properties of guided optical waves in integrated optical systems and therefore useful as circuit elements therein.

As used herein, the terms light and optical refer to the electromagnetic energy spectrum including visible light, infrared light, and ultraviolet light as to which the invention is immediately and obviously applicable; and to any immediately adjacent spectrum to which it can be extended.

In the practical utilization of the microwave spectrum, a large class of circuit elements have been developed including isolators, gyrators, cavities, directional couplers and the like.

At the present time there is a growing body of knowledge and engineering concerning the propagation of light in thin films particularly with reference to the propagation of coherent optical beams in such films and utilization of the same for practical purposes. These developments have been generally known as integrated optics, or optical waveguide systems. In general, there is a need for devices for use in the light spectrum which will correspond to the functions served by the previously enumerated microwave circuit elements.

In general, it is an object of the present invention to provide optical waveguide circuit elements for use in optical waveguides and also for the polarization conversion of guided light therein and to also provide circuit elements based on such polarization conversion for performing essential modulation and control functions in optical waveguide systems.

A further object is to provide circuit elements of the above character for use in optical waveguide corresponding to the gyrator, the isolator, a switch, and to provide an optical readout device capable of optically retrieving information stored in a magnetic media.

In order to understand the present invention, it is helpful to review and amplify certain basic concepts and to develop a theoretical basis for predicting the results of structures constructed in accordance with this invention.

In the following theoretical discussion, use is made of the editorial we, which should be understood to be taken in the first person singular. In the following general and theoretical description certain of the FIGURES and drawings herein will be used to illustrate both the structures proposed and the accompanying analysis. In that connection, the following is a brief description of those drawings utilized in connection with this theoretical description, the description of the remainder being given thereafter.

FIG. 1A is a schematic diagram illustrating schemes used in conventional laser systems for modulating a laser beam. YIG and KDP are used here as examples of magneto-optic and electro-optic substances in which the elements of the dielectric permittivity tensor may be modulated, respectively, by a bias magnetic or electric field.

FIG. 1B is a schematic diagram of a circuit element for use in integrated optics to control the propagation of a laser beam and constructed in accordance with the present invention.

FIG. 2 shows graphs depicting three basic configurations for magneto-optic effect: (a) polar, (b) longitudinal, and (c) equatorial.

FIG. 3 is a schematic diagram showing reflection and mode-conversion at the film-substrate boundary of a circuit element contructed in accordance with the present invention. In the film (waveguide medium), the two normal modes are the TE and TM modes. The various field components are denoted by: E.sub.E, H.sub.E for the incident TE wave; E.sub.M, H.sub.M for the incident TM wave; E'.sub.E, H'.sub.E for the incident TM wave; E'.sub.E, H'.sub.E for the reflected and converted TE wave; E'.sub.M, H'.sub.M for the reflected and converted TM wave. The time and spatial dependence is indicated by the exponent term at the top. In the substrate the two normal modes are evanescent waves and have all three field components E.sub.x, E.sub.y and E.sub.z. Since E.sub.x and E.sub.z are related to E.sub.y, we express them in terms of E.sub.y. The quantities E.sub.1 and E.sub.2 represent the respective E.sub.y associated with the two modes with decay constants P.sub. 1 and P.sub.2.

FIG. 4 shows graphs depicting magnitude of the mode-conversion factor plotted as a function of .sqroot.K.sub.12 for several values of .beta.. The other parameters are chosen as n.sub.f =2.50 and n.sub.s =2.00 for the respective index in the film and in the substrate. The various values of .beta. correspond to incident-angle variation from .theta.=79.degree.(.beta.=2.45) to .theta.=55.5.degree.(.beta.=2.05). The critical angle is 53.degree..

FIG. 5 shows graphs of phase angles of the reflection coefficients plotted as a function of .sqroot.K.sub.12 for several values of .beta.. The parameters used in the calculation are the same as those in FIG. 4. The dashed curves are for TE to TE reflection whereas the solid curves are for TM to TM reflection.

FIG. 6 shows the measured (dashed curve) and calculated (solid curve) output of a reflection experiment performed on a KDP-CS.sub.2 boundary. The difference between the incident angle and the critical angle is used as the abscissa. Therefore, the positive side represents the totally reflecting case whereas the negative side represents the refracting case. The refractive-index data used are; .sqroot.K.sub.11 =.sqroot.K.sub.22 =1.5064 and .sqroot.K.sub.33 =1.4664 for KDP and n.sub.f =1.6 for CS.sub.2.

FIG. 7 is a diagram showing multiple reflection in an optical waveguide. The basic scattering matrix relates the incident electric fields (TE).sub.0 and (TM).sub.0 at z.sub.0.

FIG. 8 shows diagrams a-d showing all four possible combinations by which an incident TE wave can be converted into a TM wave after making three reflections at the lower boundary and two reflections at the upper boundary.

FIG. 9 shows elements A.sub.m and B.sub.m of the multiple scattering matrix plotted as a function of the number m of complete reflections for a hypothetical case with K.sub.12 =0.01. The other material parameters used in the computation are: n.sub.f =2.50, n.sub.s =2.00, and .theta..sub.crit =53.degree.. The three cases are for three different incident angles.

FIGS. 10 show the normalized amplitude of the converted field (or the value of B.sub.m) versus normalized waveguide length L/2W tan .theta. for different incident angles. FIG. 10a is for YIG substrate with TiO.sub.2 film (solid curves) or As.sub.2 S.sub.3 glass film (dashed curves). FIG. 10b is for LiNbO.sub.3 substrate used as an electro-optic material (solid curves) or as an anisotropic material (dashed curves).

FIG. 11 depicts a light-coupling apparatus (a) experimental apparatus of a polarization converter constructed in accordance with the present invention and (b) photographs showing the waveguide mode lines resulting from operation of the apparatus of FIG. 11a. The two lines close together are the mode lines for the two near-degenerate TE and TM modes. One line was partly blocked by a polarizer for identifying the polarization. The photograph on the right (left) was taken when the input beam had a TE (TM) polarization and the output beam with the TE (TM) polarization was partially blocked. The fuzziness of the lines was due to nonuniform thickness of the photoresist film which was not polished.

FIG. 12 are diagrams showing the nonreciprocal nature of a birefringent crystal. The word "nonreciprocal" is used in the sense that the sign of the off-diagonal elements in the permittivity tensor must be changed for waves propagating in opposite directions if the righthand coordinate system is to be preserved.

FIGS. 13-20 will be described hereinafter in connection with the structural realizations of the preferred embodiments.

Recent advances .sup.1-5 in coupling a laser beam into a thin film optical waveguide and achieving laser action in thin films with an active dye ingredient have stimulated interest in integrated optical systems.sup.6. In that connection, cross reference is made to my related application entitled THIN FILM LASERS, Ser. No. 296,178, and to the continuation-in-part application thereto entitled THIN FILM DEVICES AND LASERS, Ser. No. 331,675 now U.S. Pat. No. 3,868,589 filed concurrently herewith. The present invention is directed to the enabling of integrated optical systems to be put to practical use; particularly for performing modulation and control functions on a laser beam propagating in an optical waveguide system. Specifically, the effect of gryotropic and anisotropic substrate material on the propagation of laser beams in thin film optical waveguides will be analyzed.

FIG. 1a shows arrangements used in conventional optical systems to modulate the propagation property of a laser beam as the same is passed through a magneto-optic material 10 or an electro-optic medium 11. The propagation direction of the wave is k. The initial polarization direction D.sub.1 is rotated as the wave passes through the medium to a new direction D.sub.2.

In the present invention, an optical waveguide structure is developed as shown in FIG. 1b in which the polarization conversion is accomplished at the boundary or interface between the waveguide layer 14 and either the adjacent substrate 16 which supports it or a top layer 18 disposed above it. In the drawing, the substrate is shown as having polarization-dependent propagation properties such that the off-diagonal elements in the dielectric-permittivity tensor of the substrate material provide the coupling, and hence the mechanism for conversion, between the two polarization waves (TE-transverse electric and TM-transverse magnetic) in the optical waveguide.

There are two basic requirements for the successful operation of such an optical waveguide device: (1) sufficient penetration of the evanescent wave into the substrate and (2) cumulative polarization conversion upon successive reflections at the boundary. If the penetration depth l.sub.pen is larger than the distance l.sub.guide traveled by the wave down the waveguide between two successive polarization-converting reflections, then the effective length of interaction between a laser beam and the substrate medium can be longer than the actual length of the waveguide. Thus, by operating a waveguide mode-converter with its incident angle near the critical angle one can reduce the linear dimensions of the mode-converter to sizes compatible with the integrated-optics concept.

The second requirement concerns the relative phases of the incident and converted modes. Since mode conversion for a single reflection at the film-substrate boundary is quite small, cumulative polarization conversion is essential. Significant and efficient conversion from the TE to the TM mode, or vice versa, is possible only if certain phase conditions are met. In an analysis by Tien and Ulrich,.sup.9 it has been shown that for constructive build-up of light intensity in an optical waveguide, the wave transmittedd to the guide must be in phase with the wave already existing in the guide. A similar condition exists for mode conversion. Herein is derived the condition for total mode conversion and schemes to achieve it are analyzed.

Once we have shown the practicability of mode conversion, one can consider the feasibility of practical devices. Various schemes for the realization of such devices as the gyrator, the isolator, the optical switch, and the nondestructive read-out of information recorded on magnetic films are disclosed.

Wave propagation in gyrotropic and anisotropic media has been discussed and reviewed by many workers..sup.10-12 Here we summarize the pertinent results and make them directly applicable to the thin-film optical waveguide work. At optical frequencies, the relative permeability .mu./.mu..sub.0 may be taken to be one.sup.13,14 and the permittivity .epsilon. can be represented by a permittivity tensor: ##EQU1## In the following discussion, all the media are assumed lossless. Therefore, the matrix is Hermitian and has real diagonal elements. Assuming an exp(i.omega.t-ik.multidot.r) dependence for the fields, the electromagnetic wave equation is

where K is the relative permittivity tensor and k.sub.0 =.omega..sqroot..mu..sub.0 .epsilon..sub.0 is the free-space wave number.

For waveguiding in the film, an electromagnetic wave must be totally reflected at the film-substrate boundary. If we choose a coordinate system as shown in FIG. 2 or 3, then the wave in the substrate must be decaying in the negative x direction. Therefore, the electric field in the substrate can be expressed as

where p is the decay constant normalized to k.sub.0.

The phase constant k.sub.z can be obtained from Snell's law. Since the waves in the film and the substrate have the same phase in the propagation direction, they must have the same k.sub.z. If we let k.sub.z =.beta.k.sub.0, the wave number k in the substrate can be written in terms of .beta. and p as

The negative sign in front of p.sup.2 indicates that the wave is an evanescent wave.

Now we return to the dielectric tensor of Eq. (1). Since the matrix is Hermitian, it can be diagonalized if the substrate is properly oriented with its principal axes coinciding with the x, y, and z axes of the chosen coordinate system. However, off-diagonal elements will appear under the following circumstances: (1) in a magnetooptic medium where there exists a nonvanishing magnetization, (2) in an anisotropic medium of which the principal axes are oriented along directions away from the coordinate axes, (3) in an electrooptic medium under an applied bias electric field, and (4) in an optically active medium. The off-diagonal elements are imaginary numbers in (1) and (4), and are real numbers in (2) and (3). We consider three simple cases, each with a single pair of off-diagonal elements: (a) the longitudinal configuration, K.sub.12 .noteq.0, K.sub.23 =K.sub.31 =0; (b) the polar configuration, K.sub.23 .noteq.0, K.sub.12 =K.sub.31 =0; (c) the equatorial configuration, K.sub.31 .noteq.0, K.sub.12 =K.sub.23 =0. We have borrowed the terminology from that used in describing the magnetooptic Kerr effect.

The three configurations for a waveguide structure using magnetooptic effects are illustrated in FIG. 2. In a magnetooptic medium, the off-diagonal elements can be written as K.sub.ij =i A.sub.ijl M.sub.l where A.sub.ijl is the magnetooptic constant and M.sub.l is the magnetization oriented in the l direction which is perpendicular to the ij plane. For example, if M is in the z direction (FIG. 2b), the only nonzero off-diagonal elements are K.sub.12 and K.sub.12.sup.* and thus we have the longitudinal configuration.

In an electrooptic medium, the change of dielectric constant for the linear electrooptic (Pockel's) effect can be written as .DELTA.K.sub.ij /K.sub.ii K.sub.jj.=.SIGMA.r.sub.ijl E.sub.l where r.sub.ijl is the linear electrooptic coefficient. Take KDP as an example. With the principal axes coinciding with the coordinate axes and with the applied bias field E.sub.l pointing in the z direction, the only nonzero off-diagonal elements are K.sub.12 and K.sub.21 (=K.sub.12) and we have the longitudinal configuration.

Using Eq. (3) in Eq. (2) and eliminating k.sup.2 through Eq. (4), we obtain a characteristic equation second-order in p.sup.2 (Appendix 2). The two roots of this equation are:

for the longitudinal configuration and

for the polar configuration. The three parameters are

The equatorial configuration will not be further discussed because there is no coupling between the two eigen modes in this configuration.

Now we consider the problem of reflection and mode conversion at the film-substrate boundary (FIG. 3). In the film, which is assumed isotropic, the two normal polarization modes are the TE (transverse electric) and TM (transverse magnetic) waves. Take the TE mode as an example. Since the three components E.sub.x, E.sub.y, and E.sub.z in the substrate medium are related to one another, the TE mode in the film will excite not only E.sub.y but E.sub.x and E.sub.z as well. The existence of E.sub.x and E.sub.z in the substrate medium, in turn, generates TM wave in the film. Therefore, at a film-substrate boundary, part of the TE wave is reflected and part is converted into a TM wave. A similar situation applies to the TM wave. The situation can best be described by a scattering matrix, ##EQU2## where E.sub.E, E.sub.E ', E.sub.M, and E.sub.M ' are the complex amplitudes of the various fields defined in FIG. 3. In Eq. (6), r.sub.EE and r.sub.MM are the TE to TE and TM to TM reflection coefficients and r.sub.EM and r.sub.ME are the TM to TE and TE to TM mode conversion factors.

The reflection coefficients and the mode-conversion factors can be found by equating the tangential components of E and H at the film-substrate boundary (Appendix 1) and by establishing the relations between the field components in the substrate (Appendix 2). The expressions for the matrix elements are given below: ##EQU3## The various parameters are defined as ##EQU4## The subscripts L and P refer to quantities associated with the longitudinal and polar configuration, respectively. The quantities without either subscript apply to both configurations.

The quantities important to mode conversion are the magnitude .vertline.r.sub.EM .vertline.=.vertline.r.sub.ME .vertline. of the mode-conversion factor and the phase angles (.phi..sub.EE and .phi..sub.MM) of the reflection coefficients r.sub.EE and r.sub.MM. To gain insight into the problem of mode conversion, we consider a hypothetical case where the substrate is isotropic with K.sub.11 =K.sub.22 =K.sub.33 and oriented in the longitudinal configuration; but the values of the off-diagonal elements K.sub.12 and K.sub.12.sup.* are unspecified. For isotropic media, the parameters L.sub.1 and L.sub.2 become

as an example, we use a value of n.sub.f =2.50 for the refractive index in the film and a value of n.sub.s =.sqroot.K.sub.11 =2.00 for the refractive index in the substrate. The other two parameters which are needed to completely specify the problem are the values of the off-diagonal element K.sub.12 and the incident angle .theta. (or the parameter .beta.=n.sub.f sin .theta.). We take K.sub.12 and .theta. as adjustable parameters, and plot the values of .vertline.r.sub.EM .vertline.=.vertline.r.sub.ME .vertline., .phi..sub.EE, and .phi..sub.MM in FIGS. 4 and 5 as functions of .sqroot.K.sub.12 and .beta..

Two important features of the computational results are the dependence of .vertline.r.sub.EM .vertline.=.vertline.r.sub.ME .vertline. upon the incident angle (or the value of .beta.) and the difference in the phase angles .phi..sub.EE and .phi..sub.MM. The magnitude of the mode-conversion factor .vertline.r.sub.EM .vertline. increases with K.sub.12 as expected. However, for a given material (that is, a given K.sub.12), the value of .vertline.r.sub.EM .vertline. increases as .theta. decreases toward the critical angle .theta..sub.crit =sin.sup.-1 (n.sub.s /n.sub.f) or .beta..sub.crit =n.sub.s. As .theta. approaches the critical angle, P.sub.1,2 decreases. A further penetration of the evanescent wave into the substrate increases the interaction between the TE and TM waves, and thus enhances the mode conversion.

Let us compare a thin-film waveguide mode-converter with a conventional mode-converter. For ordinary Faraday rotation (FIG. 1a), the conversion per unit distance is ##EQU5## For the waveguide mode-converter (FIG. 1b), ##EQU6## Since W is of the order of wavelength .lambda., the conversion ratio is ##EQU7##

To estimate the ratio G.sub.W /G.sub.F, we use the point P in FIG. 4 as an example. For n.sub.f =2.5, .sqroot.K.sub.11 =2.00 and .beta.=2.05, .theta.=55.degree. and .theta..sub.crit =53.degree.. At P, K.sub.12 =0.01 and .vertline.r.sub.EM .vertline.=0.18. Using these values in Eq. (10c), we find G.sub.W =4G.sub.F. This calculation shows that it is conceivable to achieve a ratio larger than 1 for G.sub.W /G.sub.F. A similar reduction in size is also possible for electrooptic modulators of the waveguide type as compared to the conventional type. As long as .theta. is within about ten degrees of the critical angle, a waveguide polarization converter can be as efficient as a conventional one.

To verify the variation of .vertline.r.sub.EM .vertline. with .theta. near the critical angle, we have performed an experiment using KDP as the substrate immersed in CS.sub.2. A He-Ne laser beam with a TM polarization was incident onto the KDP crystal, and the component of the reflected wave with the TE polarization was measured by a photomultiplier through a crossed polarizer. The KDP crystal was biased by a low-frequency electric field with a 3.6 KV/cm peak value and oriented in the longitudinal configuration. The recorder trace of the measured r.m.s. output (the dashed curve) of the photomultiplier is plotted in FIG. 6 as a function of .theta..

Since the input polarizer has a finite extinction ratio of about 10.sup.-3, the incident beam had a small TE component E.sub.1 besides a large TM component E.sub.2. The output beam, therefore, consists of two TE components one coomponent r.sub.EE E.sub.1 due to pure reflection and the other r.sub.EM E.sub.2 due to polarization conversion. Since r.sub.EE are known functions of .theta., the output can be calculated on a relative scale if the ratio E.sub.1 /E.sub.2 =A exp(i.gamma.) is known. The solid curve in FIG. 6 was calculated by using A=0.05 and .gamma.=155.degree.. The measured and calculated curves are in general agreement.

The dip in the output is caused by a rapid change in the phase difference between r.sub.EE and r.sub.EM near the critical angle. Outside .vertline..theta.-.theta..sub.crit .vertline.=0.05.degree., the phase variation becomes less rapid, and the change in the curves reflects the variation of .vertline.r.sub.EM .vertline. with .theta.. For waveguide applications, we want a device to be operated at a few degrees away from .theta..sub.crit (Sec. VI) but still close to .theta..sub.crit. The critical angle for the KDP-CS.sub.2 system is 68.42.degree.. Calculated results show that the value of .vertline.r.sub.EM .vertline. decreases by a factor of 2.8 if .theta. changes from 69.3.degree. to 73.3.degree. and by another factor of 3.3 from 74.degree. to 82.degree..

For most device applications, we would want to use magnetooptic and electrooptic materials as substrates so that wave propagation in a mode converter can be controlled by a bias magnetic or electric field. However, such devices are expected to have a small mode-conversion factor (that is, a small ratio of r.sub.ME /r.sub.EE or r.sub.EM /r.sub.MM) because of a small off-diagonal element. For example, the magnitude of K.sub.12 in YIG is calculated to be 3.times.10.sup.-4 from the Faraday rotation.sup.15 of 240.degree./cm and the off-diagonal element K.sub.ij in LiNbO.sub.3.sup.16 is about 7.times.10.sup.-4 at a bias field of 10 KV/cm. Although magnetooptic materials.sup.17,18 such as EuSe and EuO have been found to have a Faraday rotation of about 10.sup.5 degrees/cm and electrooptic materials.sup.19,20 such as SBN and KTN have been found to have large electrooptic effects, the value of K.sub.12 is still less than 0.1. To enhance magnetooptic and electrooptic effects, we propose to use a thin-film waveguide structure so that the mode-conversion effects can become cumulative. For constructive interference, certain phase conditions must be met.

For the derivation of these phase conditions, we refer to FIG. 7. At the film-substrate boundary, the field components are related by Eq. (6). In the film, the fields suffer a phase delay (.phi..sub.dE for the TE wave and .phi..sub.dM for the TM wave) in going from the point z.sub.0 to z.sub.2. At the film-top layer boundary (at z.sub.1), the fields suffer another phase change (.phi..sub.E for the TE wave and .phi..sub.M for the TM wave) but no amplitude change. Therefore, the various field components at z=z.sub.0 at z.sub.2 are related by ##EQU8##

In Eq. (11), the phase changes in the film and at the upper boundary are jointly represented by the first scattering matrix with .phi..sub.E '=.phi..sub.E -.phi..sub.dE and .phi..sub.M '=.phi..sub.M -.phi..sub.dM.

For simplicity, we assume that we start with a (TE).sub.0 wave. Upon reflection at z.sub.2, the TM wave is given by (TM).sub.3 =r.sub.ME (TE).sub.2 +r.sub.MM (TM).sub.2, or

with r.sub.EE =.vertline.r.sub.EE .vertline. exp(i.phi..sub.EE) and r.sub.MM =.vertline.r.sub.MM .vertline. exp(.phi..sub.MM). The two components in Eq. (12) are in phase if their phase angles differ by an integral multiple of 2.pi.. Since we start out with a TE wave, it is implied that the TE mode is the propagating mode of the waveguide. Therefore, the necessary conditions for phase synchronism are

where q and q' are integers.

To further examine whether the conditions stated in Eq. (13) are sufficient for subsequent reflections, we show in FIG. 8 all four possible combinations after the third reflection. The letters E and M represent the character of the wave. Adding the four components, we have

using the conditions stated in Eq. (13), we can easily show that the first three terms in Eq. (14) are indeed in phase. The fourth term is in phase opposition and it accounts for the fact that the amplitude of the TE wave progressively decreases as a result of conversion into the TM wave.

To see whether or not the conditions stated in Eq. (13) can be simultaneously satisfied, we examine the magnitudes of the phase angles. Referring to FIG. 5, we see that the values of .phi..sub.EE and .phi..sub.MM are not affected by K.sub.12 at small values of K.sub.12. Therefore, in obtaining phase information, we may set the off-diagonal elements equal to zero. As a result, the phase condition becomes identical for longitudinal and polar configurations. Also, this simplifying step enables us to treat the problems of mode conversion and phase synchronism independently. The step is taken to gain information that will guide us in working out schemes to achieve the phase-synchronous condition. In the calculations to be presented in Sec. VI for actual mode-converting structures, the exact formulas are used. with K.sub.ij =for i.noteq.j in Eq. (5), we find ##EQU9## Since K.sub.12 and K.sub.23 are assumed zero, the two evanescent modes become pure TE and TM waves. The subscripts E and M are used to indicate the character of the wave. Using Eq. (15) in Eq. (8e), we find L.sub.M =0. Thus, only the terms with L.sub.M.sup.-1 in Eq. (7) are important, and we have r.sub.EE =J.sub.E */J.sub.E and r.sub.MM =G.sub.M */G.sub.M. From Eq. (8d,b), we obtain ##EQU10## The expressions for .phi..sub.E and .phi..sub.M are similar to Eq. (16) except that the values of K.sub.11, K.sub.22, and K.sub.33 are different. In other words, there are two dielectric tensors, one for the substrate and the other for the top layer (FIG. 1). Further, K.sub.ij .ident.0 if i.noteq.j for the top layer.

The simplest way to satisfy Eq. (13) is to make the total phase changes equal. Since the film is assumed isotropic, .phi..sub.dE =.phi..sub.dM and the phase condition reduces to

if both the substrate and the top layer are isotropic, n.sub.f.sup.2 >K.sub.11 =K.sub.33. Thus, from Eq. (16), we see that .phi..sub.MM >.phi..sub.EE..sup.21 This conslusion is confirmed by the computational results shown in FIG. 5. Similarly, .phi..sub.M >.phi..sub.E. Therefore, the synchronous phase condition can be satisfied only if we use anisotropic material for either the substrate or the top layer in FIG. 1.

For an anisotropic substrate (or top layer), the phase angle .phi..sub.MM (or .phi..sub.M) can be made smaller than .phi..sub.EE (or .phi..sub.E) by choosing K.sub.33 >n.sub.f.sup.2 and yet keeping K.sub.11 and K.sub.22 <.beta..sup.2. This flexibility makes it possible to operate the mode converter degenerately with .phi..sup.E =.phi..sup.M. As an example, we consider the case for which both the substrate and the top layer are made of uniaxial crystals. The crystals are oriented with the optic axes in the z direction, so K.sub.11 =K.sub.22 =n.sub.0.sup.2 and K.sub.33 =n.sub.e.sup.2. From Eqs. (15) and (16), we find p.sub.M =(n.sub.e /n.sub.0)p.sub.E and ##EQU11## The two phase angles .phi..sub.EE and .phi..sub.MM are equal if

Since n.sub.0 <n.sub.f as required for waveguiding, it is necessary to have n.sub.e >n.sub.f >n.sub.0 (a positive uniaxial crystal). An equation similar to Eq. (19) but with a different set of values for n.sub.0 and n.sub.e applies to the top layer. Note that if Eq. (19) holds for both the substrate and the top layer, then the phase condition stated in Eq. (17) is met at all incident angles.

Another possible arrangement is one in which only the top layer is anisotropic. As an example, we choose the following values for the various parameters: n.sub.f =2.35, .beta.=2.30, and n.sub.s =.sqroot.K.sub.11 =2.25 (substrate index). From Eqs. (15) and (16), we find p.sub.E =p.sub.M =0.477, .phi..sub.EE =88.8.degree., and .phi..sub.MM =94.2.degree.. For n.sub.0 =2.25 in the top layer, .phi..sub.E =.phi..sub.EE =88.8.degree., and we need a value of 2.70 for n.sub.e in order to satisfy Eq. (17). If both the substrate and the top layer were anisotropic with n.sub.0 =2.25, a value of 2.45 for n.sub.e would be needed.

An alternate scheme to have .phi..sub.M <.phi..sub.E is to make p.sub.M <p.sub.E by letting K.sub.11 >K.sub.22 and yet keeping K.sub.11 and K.sub.22 <.beta..sup.2. As an example, we assume the top layer to be anisotropic with its optic axis oriented in the y direction (K.sub.22 =n.sub.e.sup.2 and K.sub.11 =K.sub.33 =n.sub.0.sup.2). For the same set of values for n.sub.f, .beta., and n.sub.s, we can choose n.sub.e =2.25 and n.sub.0 =2.27 for the top layer to satisfy Eq. (17). The requirement on birefringence for this case can easily be met. Practical examples of mode converters using this scheme will be given in Section VI.

The elements of the scattering matrix have some interesting and useful properties. From Eq. (7), we see that the quantities F r.sub.EE and F r.sub.MM are complex conjugate to each other. Thus, we let .vertline.r.sub.EE .vertline.=.vertline.r.sub.MM .vertline.=r.sub.11. If F=.vertline.F.vertline. exp(-i.phi.), then .phi..sub.EE +.phi..sub.MM =2.phi..

Next we examine the mode-conversion factors. For the longitudinal gyrotropic and the polar anisotropic cases, K.sub.12 is purely imaginary in Eq. (7c,d) and K.sub.23 is real in Eq. (7e,f). Thus, the quantities F r.sub.EM and F r.sub.ME are real numbers, and F r.sub.EM =-F r.sub.ME. We can write r.sub.EM =-r.sub.ME =r.sub.12 exp(i.phi.) where r.sub.12 is a real number. Further, we rewrite Eq. (11) as ##EQU12## The various matrix elements are

If the mode-propagating conditions stated in Eq. (13) are satisfied, then r.sub.EE ' and r.sub.MM ' are real and

In other words, exp(i.alpha.')=exp(-i.alpha.). Thus we can rewrite [r.sub.ij ] of Eq. (20) as ##EQU13## After 2m successive reflections at the two boundaries, the amplitudes of the two modes are given by ##EQU14##

For the longitudinal anisotropic and the polar gyrotropic cases, K.sub.12 is real in Eq. (7c,d) and K.sub.23 is purely imaginary in Eq. (7e,f). Thus, the quantities Fr.sub.EM and Fr.sub.EM are equal and imaginary. For these two cases, the relevant scattering matrices have the following forms: ##EQU15## which is the counterpart to Eq. (22a) and ##EQU16## which is the counterpart to Eq. (23a).

An important property of the scattering matrix is that under the phase-matched condition the product r.sub.EM 'r.sub.ME '=-r.sub.12.sup.2 is a negative number whereas r.sub.EE '.sup.2 and r.sub.MM '.sup.2 are positive numbers. The opposite sign means physically that if we start with a TE wave, the TE created by the converted TM wave through the process TE.fwdarw.TM.fwdarw.TE is always in phase opposition to the original TE wave.* This fact ensures a continuous conversion of TE into TM mode and vice versa. Another property of [r.sub.ij '].sup.m is that after many reflections, the diagonal elements remain real and the phase angles of the off-diagonal elements remain unchanged. Since there is a continuous conversion between TE and TM modes, we might expect an ultimate situation in which the energy would be equally divided between the two modes. The mode conversion does not stop, however, at the equal energy point. If A.sub.m =B.sub.m =1/.sqroot.2 in Eq. (23a), then

The multiple scattering matrices of Eq. (23) can be obtained by using Sylvester's theorem.sup.22 (Appendix 3). Since all the media are assumed lossless, electromagnetic energy must be conserved so that r.sub.11.sup.2 +r.sub.12.sup.2 =1. With this relation, the matrix elements are found to be

for most magnetooptic and electrooptic media, .vertline.r.sub.EM .vertline.<<1 and .vertline.r.sub.EE .vertline..apprxeq.1. Thus, we approximate tan.sup.-1 (.vertline.r.sub.EM .vertline./.vertline.r.sub.EE .vertline.) by r.sub.EM. This approximation is used in obtaining Eq. (25).

Since Eq. (25) is based on .phi..sup.E =.phi..sup.M, the conditions as stated in Eqs. (13) and (17) are not merely for cumulative mode conversion but also for complete mode conversion. To verify the sinusoidal variation of mode conversion, we evaluate the elements in Eq. (23) with the aid of a computer. FIG. 9 shows the computational results of the matrix elements A.sub.m (dashed curve) and B.sub.m (solid curve) as a function of the number m of reflections at the film-substrate boundary. Since the distance between z.sub.2 and z.sub.0 is 2W tan .theta. in FIG. 7, the number is simply equal to m=z/(2W tan .theta.). The point P.sub.1 (circled) indicates the point of half conversion (energy-wise) whereas the point P.sub.2 (also circled) marks the point of full conversion from one mode to the other.

There are three curves for B.sub.m in FIG. 9. We used the following values: n.sub.f =2.50, .theta..sub.crit =53.degree. for the critical angle (or n.sub.s =2.00), and K.sub.12 =0.01 for all three curves. The three curves are computed under phase-matched condition and for different values of the incident angle .theta. with (I) .theta.=61.67.degree. and .beta.=2.2, (II) .theta.=55.11.degree. and .beta.=2.05, and (III) .theta.=73.78.degree. and .beta.=2.40. As we can see, the nmber m of reflections required to reach full conversion decreases from a value of 650 to 90 as we change .theta. from 73.78.degree. to 55.11.degree.. This is a result of increasing penetration of the evanescent wave into the substrate.

The computational result shown in FIG. 9 presupposes .phi..sup.E =.phi..sup.M. To show that the synchronous phase condition can be satisfied in real cases and that the construction of waveguide mode converters is feasible, we present in Table 1 the performance characteristics of several mode-converting structures using materials with known dielectric properties. The angle .theta..sub.m is the incident angle under the phase-matched condition. From Eqs. (5), (7), and (8), the values of .phi..sub.EE and .phi..sub.MM can be computed for a selected value of .theta. or .beta. (=n.sub.f sin .theta.). The computation is repeated until .phi..sup.E =.phi..sup.M. Once .theta..sub.m is found, the value of .vertline.r.sub.EM .vertline. can be calculated.

Complete conversion occurs at m.vertline.r.sub.EM .vertline.=.pi./2 when B.sub.m =1 in Eq. (25). Since there are l/(2W tan .theta.) mode-converting reflections in a distance l, the waveguide length necessary for complete mode conversion is ##EQU18## The film thickness W can be computed from Eq. (13). The mode number q=q' is so chosen that W is in the neighborhood of 0.5 micron. As we can see from Table 1, the value of l.sub.0 falls in a very convenient range. The dimension is small so it suits the concept of integrated optics, and yet large enough so the value of l.sub.0 can be accurately controlled.

Complete mode conversion is possible only when the phase-matching condition is satisfied. If a phase mismatch .delta..phi.=.phi..sup.M -.phi..sup.E .noteq.0 exists, only a partial conversion an be achieved. It can be shown (Appendix 3) that a general expression for B.sub.m is

where .phi..sub.2 =cos.sup.-1 [r.sub.11 cos (.delta..phi./2)]. Under the phase-matched condition .delta..phi.=0, .phi..sub.2 .apprxeq.sin .phi..sub.2 =r.sub.12 and Eq. (27) reduces to Eq. (25). For small .delta..phi., Eq. (27) can be approximated by ##EQU19## where G.sub.0 =.vertline.r.sub.EM .vertline..sup.2 and G=.vertline.r.sub.EM .vertline..sup.2 +.vertline.r.sub.EE .vertline..sup.2 (.delta..phi./2).sup.2. When .delta..phi.=2.vertline.r.sub.EM .vertline./.vertline.r.sub.EE .vertline., maximum achievable conversion (that is, maximum B.sub.m) is 1/2 and requires a distance l=l.sub.0 /.sqroot.2. Since .vertline.r.sub.EE .vertline..apprxeq.1 for most cases, a practical condition for near-complete mode conversion can be simply stated as

fig. 10 shows the magnitude of B.sub.m (that is, the converted field amplitude as a fraction of the initial field amplitude) versus the distance down the waveguide under .delta..phi.=0 and .delta..phi..noteq.0 for several mode-converting structures. As these structures have been discussed elsewhere,.sup.32 it suffices to point out that the behavior of mode conversion predicted by Eqs. (25) and (28) has been confirmed by computational results.

The synchronous phase condition .delta..phi.=0 can be satisfied only at a specific incident angle .theta..sub.m for a given W. If .theta. deviates from .theta..sub.m, there will be a cumulative phase change .delta..phi. between the converted fields from successive mode-converting reflections because .delta..phi..sup.E /.delta..theta. and .delta..phi..sup.M /.delta..theta. are different. The tolerance in incident angle misalignment can be represented by a bandwidth .DELTA..theta. defined by the half-power points. We note from Table 1 that the values for .DELTA..theta. are of the same order of magnitude as, or larger than, the tolerance.sup.2 for light coupling into an optical waveguide. Therefore, the alignment necessary for mode converting should not be more difficult than for light guiding.

The tolerance set on .DELTA..theta. automatically imposes a condition on the tolerable variation in film thickness because Eq. (13) which governs the selection of propagating TE and TM modes relates W to .theta.. The synchronous phase condition .phi..sup.E =.phi..sup.M is met only at a certain film thickness W.sub.m. For W different from W.sub.m, the incident angles .theta..sub.E and .theta..sub.M become different for the TE and TM waves in order that Eq. (13) can be satisfied. As a result, .phi..sup.E .noteq..phi..sup.M and the guided phase velocities for the TE and TM waves are different. Consider the process TE.fwdarw.TM conversion. The converted TM field has the phase constant k.sub.z.sup.E of the parent TE field, and yet the eigen TM mode of the guide has a phase constant k.sub.z.sup.M. As a result of k.sub.z.sup.E and k.sub.z.sup.M being different, we have a problem of phase coherence. The problem is similar to that encountered in parametric interaction. In analogy, we define a coherence length ##EQU20## such that .vertline.n.sub.f (.beta..sub.M -.beta..sub.E)k.sub.0 .vertline.l.sub.c =.pi.. In Eq. (30), D.sub.E =A+.differential..phi..sup.E /.differential..theta., D.sub.M =A+.differential..phi..sup.M /.differential..theta., A=2k.sub.0 n.sub.f W.sub.m sin .theta. and .delta.W=W-W.sub.m. Equation (30) is obtained by finding .delta.W/.delta..theta. from Eq. (13) and then relating .beta..sub.M -.beta..sub.E to .delta..theta.. For .delta.W.noteq.0, .delta..phi. is also nonzero and only a partial conversion is possible. The value of B.sub.m drops to 1/.sqroot.2 (half-power points) when l.sub.c =l.sub.0. We use this definition to calculate the bandwidth .DELTA.W, that is, the tolerance in film thickness (Table 1). With optical polishing to .lambda./50 being currently available, we feel that the construction of these structures, with the exception of structure 4, is feasible.

The values of .DELTA..theta. and .DELTA.W are influenced primarily by two factors: (1) the differential phase change (.differential..phi..sup.M /.differential..theta.)-(.differential..phi..sup.E /.differential..theta.) with respect to the incident angle and (2) the magnitude of the mode-conversion factor .vertline.r.sub.EM .vertline.. If .delta..phi.=.phi..sup.M -.phi..sup.E is less sensitive to .theta., the two eigen modes will remain degenerate over a wider range of .theta. and W. A substrate with a larger .vertline.r.sub.EM .vertline. requires less reflections at the film-substrate boundary to complete mode conversion. For m reflections, the cumulative phase shift is .delta..phi.. The converted field ceases to grow when m.delta..phi.>.pi. or .delta..phi.>.pi./m. Thus, for larger .vertline.r.sub.EM .vertline. (which means smaller m), we can allow a larger .delta..phi. per mode-converting reflection. This is obvious from Eq. (29). As the data in Table 1 show, substrate materials with larger off-diagonal elements should give better performance characteristics. Further, the longitudinal configuration has better performance characteristics than the polar configuration as a result of a larger .vertline.r.sub.EM .vertline..

Besides large off-diagonal elements, substrate materials should also have a small birefringence. If K.sub.11 is quite different from K.sub.22, then K.sub.22 can be close to .beta..sup.2 in Eq. (15) while K.sub.11 is still away from .beta..sup.2. A large difference between P.sub.1 and P.sub.2 has two effects: a large differential phase change .differential.(.delta..phi.)/.differential..theta. and a smaller .vertline.r.sub.EM .vertline. for the same .vertline.K.sub.12 .vertline.. Both these effects are not desirable. The poor performance of structure 4 is largely due to the large birefringence of LiNbO.sub.3. The off-diagonal elements being the same, a higher .vertline.r.sub.EM .vertline. would require a deeper and more equal penetration of the two eigen modes into the substrate.

The choice of the film material is determined by the index of refraction. As mentioned earlier (FIGS. 4 and 9), the value of .vertline.r.sub.EM .vertline. increases as .phi. approaches the critical angle. The phase angles .phi..sup.E and .phi..sup.M, however, change rapidly near .theta..sub.crit. if we choose an n.sub.f such that .theta..sub.m is close to .theta..sub.crit, we raise not only the value of .vertline.r.sub.EM .vertline. but also the value of .differential.(.delta..phi.)/.differential..theta.. We find that the performance of a mode converter improves as .theta..sub.m and .theta..sub.crit move toward each other. However, when .theta..sub.m approaches either the critical angle of the substrate or the top layer, we witness a trade-off between .vertline.r.sub.EM .vertline. and .differential.(.delta..phi.)/.differential..theta.. The crystal SBN has a relatively large off-diagonal elements, and yet structure 7 has relatively small tolerances. The reason is that .theta..sub.m is within one degree of the critical angle at the substrate-film boundary.

The importance of selecting film index can be seen by comparing .DELTA..theta. and .DELTA.W of structures 7 and 9. We can rewrite Eq. (29) as .DELTA..theta.<.vertline.r.sub.EM .vertline./(.differential..phi./.differential..theta.). The value of (r.sub.EM).sub.L is larger than (r.sub.EM).sub.P by a factor of about 6 under similar operating conditions. If the same film index n.sub.f =2.39 were used in the polar configuration, we would expect a value of about 0.4.degree. for .DELTA..theta.. In Eq. (30), since l.sub.0 is proportional to .vertline.r.sub.EM .vertline., the main factor in determining .DELTA.W is also .vertline.r.sub.EM .vertline./(.delta..phi./.delta..theta.). Therefore, a value of about 450 A for .DELTA.W should be possible.

From the several mode-converting structures shown in Table 1, it appears optimum to have the film and top-layer indices so chosen that .theta..sub.m is above 75.degree. and the difference .theta..sub.m-.theta..sub.crit is about 4.degree.. The small tolerances for structure 8 are largely due to the field-induced anisotropy in KTN which makes phase-matching difficult. By reducing the bias field by a factor of 3, we can almost restore isotropy in KTN, and still have a relatively large value of K.sub.23 =2.23.times.10.sup.-3. With this relaxation on the requirement of the top-layer anisotropy, it should be possible to improve the performance of structure 8.

Now, we comment on the selection of the scheme for phase compensation. As discussed in Sec. IV, the two possible schemes are: (1) a positive uniaxial crystal oriented with K.sub.11 =K.sub.22 <K.sub.33 and (2) a positive or negative uniaxial crystal oriented with K.sub.11 >K.sub.22. The phases .phi..sup.E and .phi..sup.M vary more equally in the case with K.sub.11 =K.sub.22 than in the case with K.sub.11 >K.sub.22 as it is obvious from Eq. (15). Therefore, scheme 1 is expected to yield much larger tolerances in .DELTA..theta. and .DELTA.W than scheme 2. Calculations based on scheme 1 for GdIG mode converter show that an improvement by a factor of 10 would be obtained if a hypothetical crystal with a relatively large birefringence were used. The requirement on birefringence can be relaxed somewhat by making P.sub.E and P.sub.M small. Since common positive uniaxial crystals do not have a large birefringence (with tellurium being a notable exception), we have used scheme 2 for phase compensation. Scheme 2 has a requirement on birefringence which most unaxial crystals can meet. However, scheme 1 may be useful in other wavelength regions where suitable positive uniaxial crystals should become available.

As a preliminary test of the theory, we have performed a light-coupling experiment (FIG. 11) on a waveguide made of isotropic photoresist film and anisotropic quartz substrate. A He-Ne laser beam was coupled into and out of the waveguide through an isotropic Pb-glass prism. The use of isotropic materials for the film and prism removes any ambiguity as to the origin of polarization conversion observed in the experiment. The quartz crystal was oriented in the longitudinal configuration with its c-axis lying in the xy plane and making an angle of 5.degree. with the normal of the film surface. Since the case is similar to structure 1 of Table 1 with a film index (n.sub.f =1.6) close to the indices of quartz, a large tolerance .DELTA.W is expected. This large tolerance should result in a near-degenerate operation of the TE and TM modes, ensuring a significant polarization conversion.

As can be seen from FIG. 11b, there was indeed a significant polarization conversion. The fact that the two TE and TM mode-lines were close together indicates that the two modes were nearly degenerate. In our experimental procedure, the waveguide was first decoupled from the prism by releasing the pressure on the substrate, and the polarization of the input beam was so adjusted that the polarization of the output beam was purely TE or TM. Then the waveguide was pressed back against the prism to obtain mode coupling, and the mode-lines were analyzed for their polarization by a polarizer. The simultaneous presence of the two TE and TM lines close together, when the input beam had pure TE or TM polarization, strongly suggests that significant polarization conversion had taken place. A similar experiment with an isotropic glass slide (n.sub.s =1.52) as the substrate was performed, and only mode-lines with a polarization identical to the input-beam polarization were observed.

In the foregoing analysis, use was made of the concepts of geometrical optics to explain this invention. Such an explanation is taken from a point of view using reflection from the boundaries of the waveguide and the adjacent anisotropic or gyrotropic material. The following discussion will generalize the foregoing by making use of a wave optics point of view and therefore will facilitate the understanding and application of the invention to a wider range of possible structures.

In wave optics, wave propagation in the waveguide can be described by two eigen modes herein referred to as mode A and mode B. Away from the phase-matching condition, the A and B modes have pure TE and TM polarization as a limit; for example, A mode is pure TE, and B mode is pure TM or vice versa. Near the phase-matching condition, the interaction between the TE and TM polarization becomes so strong that the A and B modes contain a substantial amount of the other polarization. When the phase-matching condition .phi..sup.E =.phi..sup.M is satisfied, the TE and TM contents in the A and B modes are nearly equal. Under the circumstances one can express the complex amplitudes of the two modes as

where either a.sub.2 /a.sub.1 =b.sub.2 /b.sub.1 =1 or a.sub.2 /a.sub.1 =b.sub.2 /b.sub.1 =i. The A and B modes have different phase velocities .omega./k.sub.z.sup.A and .omega./k.sub.z.sup.B. As the two modes propagate in the waveguide, the composition of the resultant wave changes with distance. For example, when (k.sub.z.sup.A -k.sub.z.sup.B) z=even multiple of .pi. the resultant wave A+B becomes pure TE due to mutual cancellation of the TM components over the properly chosen interaction length.

The above conclusion based on wave-optics analysis is perfectly general, irrespective of which medium (the substrate, the film or the top layer) is made of material capable of coupling the TE and TM polarization. In other words, wave-optics analysis predicts TE to TM or TM to TE conversion in optical waveguides where one of medium has off-diagonal permittivity tensor elements. The conversion from one polarization to the other is nearly complete if .phi..sup.E =.phi..sup.M.

As used herein and in the claims, anisotropic and gyrotropic should be taken as subclassifications of materials whose propagation characteristics for electromagnetic waves are polarization dependent. A more complete classification is as follows.

Anisotropic materials refer to those in which wave propagation is described by linearly polarized waves. Anisotropic materials are also characterized in that the off-diagonal elements of the associated permittivity tensor are real and/or the on-diagonal elements are unequal.

Anisotropic materials may be either inherently anisotropic or may be induced by application of external fields. Examples of inherently anisotropic materials include quartz and calcite. Examples of induced anisotropic materials include the electro-optic materials: Lithium niobate and LiTaO.sub.3 ; and the photo-elastic materials: silicon, GaAs and rocksalt.

Gyrotropic materials refer to those in which wave propagation is described by circularly polarized wave and are further characterized in that the off-diagonal elements of the permittivity are purely imaginary. Gyrotropic materials may be further classified as those which are inherently so, commonly known as optically active materials and those which may be induced to become so, commonly known as magneto-optic materials. Examples of inherently optically active materials include: quartz and lactic acid in liquid solution. Examples of induced gyrotropic materials or magneto-optic materials include the following: YIG (Yttrium iron garnet) and mixed magnetic garnet Eu Er.sub.2 Fe.sub.4 Ga O.sub.12.

It will be noticed that in the examples given for induced anisotropic and induced gyrotropic materials the anisotropic case is generally developed by the application of electric fields, or mechanical fields such as mechanical stress, while the gyrotropic case is generally induced by the application of magnetic fields.

The foregoing objects and features of the present invention will become apparent from the following description of particular devices when taken in conjunction with the remaining accompanying drawings and claims.

FIG. 1b shows a cross section of a mode converter constructed in accordance with the present invention.

FIGS. 2 through 12 have been previously described.

FIG. 13 is a schematic perspective view illustrating a mode converter using microwave biasing field constructed in accordance with the present invention.

FIG. 14 is a schematic perspective view illustrating the mode converter using a microwave and biasing field of different field orientation, and constructed in accordance with the present invention.

FIG. 15 is a schematic perspective view illustrating a mode converter constructed in accordance with the present invention and using an acoustic biasing field.

FIG. 16 is a schematic diagram of optical waveguide switches constructed in accordance with the present invention and employing a polarization converter as disclosed herein.

FIGS. 17 are (a) a schematic diagram of an optical waveguide gyrator constructed in accordance with the present invention, (b) the functional structure, (c) the input-output relation and (d) the circuit representation of a thin-film-waveguide gyrator. The magnetic and dielectric mode converters must form a complementary pair.

FIG. 18 is a schematic diagram of an optical waveguide isolator constructed in accordance with the present invention.

FIG. 19 is a cross sectional view of an optical readout and magnetic memory constructed in accordance with the present invention.

FIG. 20 is a cross sectional view of another form of optical readout and magnetic memory constructed in accordance with the present invention.

Referring again to FIG. 1b one can now generalize the nature of the polarization converter constructed in accordance with the present invention as consisting of a substrate layer 16 which serves as support for the converter and which possesses a predetermined permittivity tensor .epsilon..sub.s. Upon this substrate is disposed, in intimate contact, a thin film optical waveguide 14 which also possesses a predetermined permittivity tensor .epsilon..sub.f. Above this layer is disposed means 18 having a predetermined permittivity tensor .epsilon..sub.t. The latter means forming the top layer may consist of a material or may be, for example, air. In some circumstances, it is envisaged that this means may consist of the absence of any material as, for example, a vacuum. In the sense of the present application, a vacuum defines a permittivity tensor as well as an index of refraction and is capable of supporting the propagation of electromagnetic waves and accordingly can be useful in certain applications of the invention as, for example, in outer space. While the foregoing discussion defines each of the layers as having a certain permittivity tensor it will be understood that such a characterization is necessitated by the requirement that one of the layers be anisotropic or gyrotropic and therefore polarization dependent.

In the more common sense, each of the substrate and top layers must possess an index of refraction which is less than that of the thin film layer so that the necessary condition for the establishment of the thin film layer as a waveguide having reflective boundaries can be met. In a more general and precise sense this same condition can be stated as follows: that the permittivity tensor .epsilon..sub.f of the film has certain diagonal elements which control the propagation of the wave and that these diagonal elements are selected to be greater than the corresponding diagonal elements of the permittivity tensors .epsilon..sub.t and .epsilon..sub.s. In the context of the present application, this requirement defines the waveguide. One of the layers, either the waveguide itself or the substrate or the top layer, is made of a material having propagation characteristics which are polarization dependent to thereby cause coupling between the TE and TM polarizations of the wave propagated in the waveguide. This also establishes two eigen modes of propagation through the waveguide each of which is a combination of TE and TM polarization waves, as previously discussed. These eigen modes have different phase velocities and accordingly are found to add and subtract from each other at intervals. The waveguide film 14 is made with a thickness such that the phase-matched condition .phi..sup.E =.phi..sup.M is satisfied, and the two eigen modes therefore contain nearly equal TE and TM components. In addition, the material having the anisotropic or gyrotropic property is made with a length, herein called interaction length, such that the relative phase of the eigen modes changes over said length and causes mutual cancellation of either the TE or TM components from both of the modes. Accordingly, conversion of TE to TM or TM to TE waves results as they are propagated through an interaction length, l, of this material. As previously discussed, the material having the property of propagation characteristics which are polarization dependent can be anisotropic or gyrotropic, either induced or inherently so.

Once the principle of polarization conversion in optical waveguides is known it can be extended to cases which are similar but have not been discussed in the present mathematical analysis. One obvious extension is the case where the film, instead of the substrate, is made of material capable of coupling the two polarizations. For this case, the fields in the substrate can be described as pure TE and TM waves, and the fields in the film can be described as TE-rich (with E.sub.y being large and E.sub.x and E.sub.z being small) and TM-rich (with E.sub.x and E.sub.z being large and E.sub.y being small) waves. Therefore, this situation is analogous to the one already analyzed except that the nature of the fields in the film and the substrate is interchanged. For waveguide polarization-converters using anisotropic gyrotropic, magneto-optic, or electro-optic material for the film, the condition for complete polarization conversion is:

which is a generalization of Eq. (17).

Another interesting extension is shown in FIGS. 13 and 14 where the applied bias electric field is of high frequency (microwave) instead of d-c or zero frequency as assumed in the mathematical analysis. Thus, a mode converter such as in FIG. 13 including a substrate 20, waveguide layer 22 and top layer 24 is sandwiched between the conductive strips 26, 28 of a microwave strip line to which a microwave source 30 is connected.

The E field lines are shown diagramatically as extending between the strips, their closeness of spacing indicating the strength of the field. Since the strips are above and below the plane of the mode converter the E field lines pass normal to that plane. FIG. 14 shows another arrangement in which the conductive strips 26', 28' of the strip line are disposed on each side of the mode converter so that the E field is oriented in alignment with the plane of the waveguide. Both of the structures of FIGS. 13 and 14 make use of electro-optic materials in one of the layers 20,22,24.

For example, if a microwave bias field is used, then the condition for phase synchronism in the z direction becomes

instead of

for the d-c or zero frequency case, where K.sub.z is the wave vector of the microwave bias field. It is also possible to use photoelastic effect to produce polarization conversion when an acoustic or elastic wave is present. For this case, K.sub.z om Eq. (17b) is the wave vector of the acoustic or elastic wave. For nonzero K.sub.z the condition for complete conversion becomes:

the requirement is that .delta..phi. be very small but not zero. Once a mode-converting structure is made, the values of .phi..sub.EE +.phi..sub.E -.phi..sub.dE and .phi..sub.MM +.phi..sub.M -.phi..sub.dM are fixed and so are the values of k.sub.z.sup.E and k.sub.z.sup.M which are now unequal.

For complete polarization conversion, we choose the frequency of the microwave bias field such that Eq. (17b) is satisfied. These arrangements have two obvious advantages over the arrangement using a d-c bias field. First the mode converter can be operated at a high modulating frequency. Second, the tolerance in the film thickness is greatly relaxed because we now have the additional flexibility of using the modulating frequency or K.sub.z as an adjustable parameter in satisfying the condition for complete polarization conversion.

The same arrangements employing microwave strip lines can also be used for polarization converters made of magneto-optic materials. In this case, the microwave magnetic field induces off-diagonal elements in the permittivity tensor to cause coupling between the TE and TM waves. FIG. 15 shows a possible arrangement for a polarization converter using the photo-elastic effect.

Thus a mode converting structure is provided consisting of a substrate 32 on which an optical waveguide (34') is disposed. The top layer can be air. The substrate 32 is made of a photo-elastic material and extends laterally away from and on each side of the waveguide. Means are provided for establishing an acoustic wave field over a certain length of the converter and include a pair of acoustic transducers 36, 38 disposed in coupled relation to the substrate on opposite sides of the waveguide. The transducers are mounted at spaced locations with respect to the orientation of the waveguide so that the resultant acoustic wave vector K makes an angle to the wave vector K.sub.z in the waveguide and therefore defines an average interaction length l. An elastic wave is launched into the substrate through the use of an acoustic wave transducer. A standing wave pattern of stress field is set up between the two launchers particularly in the region of the interface between the waveguide and the substrate. It is desirable to have the acoustic wave vector making a slight angle with the direction of light wave propagation so that the intersection of the acoustic wave front with the light waveguide clearly defines a region of interaction between the stress wave and the light wave. The length of this interaction region determines the interaction length.

Thus, a thin-film optical waveguide as in FIG. 1b using electro-optic substrate results in a mode converter from which the output polarization TM or TE can be switched by the application of a bias electric field. Such a mode conversion device can also be used to combine two signals originally entering the waveguides as TE and TM modes into a single mode either TE or TM for signal processing as will be now discussed.

Referring to FIG. 16a there is shown a polarization converter 40 located in an optical waveguide which terminates in TE and TM channels 41, 42 which are mutually exclusive, that is to say which only transmit TE or TM waves exclusively and reject the opposite wave. Input light waves are launched through the device from the opposite side. Means are provided for supplying a bias electric field which can be controlled either by being reversed or by having a two valued output as indicated at 43a. A phase shifter 44 will usually have to be employed between the input and the converter in order to assure proper phase relationship between the TE and TM waves before they enter the converter. For example, from Eq. (23a), we have

an interesting case is one in which (TM).sub.0 =(TE).sub.0 exp. (-i.alpha.) through the adjustment of the phase shifter 44. If the converter is operated at the half-conversion point, then A.sub.m =1/.sqroot.2 and B=.+-.1/.sqroot.2 depending on the direction of the bias field. For B.sub.m =-1/.sqroot.2, (TM).sub.2m =(TE).sub.0. In other words, the incoming signal contains equal TE and TM signals. Each of these is split after passing through the converter at half conversion into four equal energy TE and TM components. Since the phase is controllable, these split signal components are combinable to cancel either polarization leaving TE or TM as the output wherever the applied biasing signal is present.

FIG. 16b shows a similar arrangement operated at full conversion, like parts being similarly designated with the addition of the prime ('). If an input TE wave is supplied, the output will be TE when the converter is off, and will be TM when the converter is on. Initially this provides for selection of wave types, and after passage through the branch channels 41', 42' results in the signal in either branch being alternately "on" or "off" depending on the state of the applied bias field.

The polarization converter thus conceived has an important property, that is, nonreciprocal propagation characteristics. There are four possible configurations of the polarization converter: (1) longitudinal (with off-diagonal elements in the transverse plane) gyrotropic, including magneto-optic, (2) polar (with off-diagonal elements in the longitudinal plane) gyrotropic, (3) longitudinal anisotropic, and (4) polar anisotropic. Wave propagation in a mode-converting waveguide can be described by a 2.times.2 scattering matrix. If the wave reverses its direction of propagation, the off-diagonal elements change their sign in configurations (1) and (4) but remain the same in configurations (2) and (3). Using the foregoing it is possible to obtain realization of optical gyrator and isolator. The gyrator and isolator were made practical devices first in the microwave region..sup.33 The use of the Faraday effect for signal isolation was later extended to the optical region .sup.34, 35 for conventional laser systems. Here it is proposed to use the nonreciprocal nature of waveguide mode-converters to accomplish these circuit functions for use in integrated optics.

For nonreciprocal device applications, the following observation is useful. The multiple scattering matrices presented in Sec. V can be classified into two groups: [r].sub.I and [r].sub.II'. ##EQU21## with [r].sub.I for the longitudinal magnetic and polar dielectric cases and [r].sub.II for the polar magnetic and longitudinal dielectric cases. The sign of B changes for both the longitudinal magnetic and dielectric cases and remains the same for both the polar magnetic and dielectric cases if the wave reverses its direction of propagation. This sign reversal or constancy can be seen from Eq. (7) as .beta. changes to -.beta.. The sign of B changes, of course, if the sign of K.sub.12 or K.sub.23 changes when the applied bias field reverses its direction.

In the magneto-optic and electro-optic cases, the sign of K.sub.12 or K.sub.23 is fixed once the direction of the bias magnetic or electric field is specified. To visualize the situation in an anisotropic medium, we refer to FIG. 10. The three principal crystal axes are so chosen that .epsilon..sub.a =.epsilon..sub.0 K.sub.1, .epsilon..sub.b =.epsilon..sub.0 K.sub.2, and .epsilon..sub.c =.epsilon..sub.0 K.sub.3, For simplicity, we consider the case in which the coordinate x and y axes are chosen as the axes bisecting the angles between a and b. It can be shown by matrix transformation that the relative permittivity tensor is given by [K].sub.a of FIG. 12 in the chosen coordinate system (x, y, z) with K.sub.11 =(K.sub.1 +K.sub.2)/2 and K.sub.12 =(K.sub.2 -K.sub.1)/2.

If we keep the same coordinate system for waves propagating in the +z and -z direction, the permittivity tensor remains the same but the sign of .beta. changes for the two waves. On the other hand, if we keep +.beta. for both waves, then the coordinate axes must change from (x,y,z) to (x', y', z') with z' in the -z direction. Note that waves propagating in the -z direction are still incident on the film-substrate boundary from the film side. Therefore, x'=x and we must make y'=-y in order to have a right-hand coordinate system. In the coordinate system (x', y', z'), the permittivity tensor is given by [K].sub.b of FIG. 12. This analysis shows that we must change either the sign of K.sub.12 or that of .beta. in Eq. (7). In either case, the sign of (r.sub.EM).sub.L changes if the wave reverses its direction of propagation. By the same reasoning, the sign of K.sub.23 remains unchanged for the polar configuration.

Now we refer to the arrangement of FIG. 17a in which a gyrotropic (or more specifically magnetic) mode converter 50, a phase shifter 52, and an anisotropic (or more specifically dielectric) mode converter 54 are arranged in tandem. Mutually exclusive waveguide pairs 55,56 and 57,58 are connected at each end of the unit. To be more specific, we further assume that the mode converters 52,54 are operated in the longitudinal configuration. Thus, the magnetic and dielectric mode converters can be represented by [r].sub.I and [r].sub.II, respectively. For the phase shifter, we let the scattering matrix be ##EQU22## For waves propagating in the +z direction, the resultant matrix is ##EQU23##

For waves propagating in the -z direction, we have two changes. First, the signs of B.sub.1 and B.sub.2 in [r].sub.I and [r].sub.II change as k.sub.z changes from .beta. to -.beta. for both longitudinal magnetic and dielectric cases. We denote the new matrix by [r].sub.I ' and [r].sub.II '. Second, the order of multiplication is reversed. Thus, we have ##EQU24##

If we have a complementary pair of magnetic and dielectric mode converters such that A.sub.1 =B.sub.1 =A.sub.2 =B.sub.2 -/.sqroot.2 and .alpha..sub.1 =.alpha..sub.2 =.alpha. and if a'=90.degree., then the two matrices reduce to ##STR1## Therefore, for waves propagating in the +z direction, the TE wave remains a TE wave and the TM wave remains a TM wave. For waves propagating in the -z direction, however, the incoming TE wave changes into an outgoing TM wave and the incoming TM wave changes into an outgoing TE wave. This situation is schematically illustrated in FIG. 17c. As we can see, such a device behaves like a gyrator.

We should point out that a complementary pair of mode converters together with a phase shifter can also be operated in the polar configuration. The choice between longitudinal and polar configuration depends on the relative ease in depositing the film and in applying the bias fields. It is also possible to form a complementary pair without the phase shifter. The possible combinations are (1) a longitudinal magnetic mode converter and a polar dielectric mode converter and (2) a polar magnetic mode converter and a longitudinal dielectric mode converter. Take case (1) as an example. For a dielectric polar mode converter, ##EQU25## which is the same for both directions of propagation. Thus, a combination of I and IV gives the following characteristics: ##EQU26## for waves propagating in the +z direction and ##EQU27## for waves propagating in the -z direction. Equation (37) is obtained under the complementary-pair condition that .alpha..sub.1 =.alpha..sub.2, A.sub.1 =A.sub.2 =1/.sqroot.2, and B.sub.1 =B.sub.2 =1/.sqroot.2.

By a complementary pair of polarization converters is meant that each one of the pair splits the energy of each component of the incoming wave into one-half TM and one-half TE waves. In addition, the relative sign of the TE and TM components passing through one of the pairs remains unchanged as a function of direction of passage of the wave; while the relative sign of the TE and TM components of the other of the pair reverses when the direction is reversed. Each of the complementary pairs listed have this property. If the phase mode is correct or if a phase shifter is employed between the pair the arrangement will first turn a TE wave, for example, entering the input into equal TE.sub.1 and TM.sub.2 waves and then split them again into TE.sub.12 and TM.sub.12 waves and into TE.sub.22 and TM.sub.22. By adjusting the relative phase of these components TE.sub.12 and TE.sub.22 add, while TM.sub.12 and TM.sub.22 subtract or cancel.

For a TE wave entering from the output or right side, the effect is similar except that the TE.sub.12 and TE.sub.22 waves cancel while the TM.sub.12 and TM.sub.22 waves add.

In summary, the four useful combinations of the various configurations are:

(A) A longitudinal gyrotropic mode converter +a phase shifter +a longitudinal anisotropic mode converter

(B) A polar gyrotropic mode converter +a phase shifter +a polar anisotropic mode converter

(C) a longitudinal gyrotropic mode converter +a polar anisotropic mode converter

(D) a longitudinal anisotropic mode converter +a polar gyrotropic mode converter

If the waveguide length is chosen properly, the polarization of the wave remains unchanged (say, TE.fwdarw.TE and TM.fwdarw.TM) for one direction (say, from right to left) of propagation but changes to the other polarization (say, TE.fwdarw.TM and TM.fwdarw.TE) for the other direction (say, from left to right) of propagation.

It is well known in microwave circuitry that one can make an isolator out of a gyrator. Absorbers 60, 61 are loaded into the TM waveguides 56, 58 at the two ends of the gyrator as shown in FIG. 18 so that only one polarization mode (say, TE) is allowed to propagate through to the right end whereas only the other polarization mode (say, TM) is allowed to propagate through to the left end. With this arrangement, a wave can propagate only in one direction (say, entering the left end as a TM wave, converted into a TE wave in the gyrator and then propagating through the right end as a TE wave) but not in the opposite direction.

The principle of using multiple reflections in a waveguide structure for enhanced polarization conversion as herein disclosed can also be applied to magneto-optic memory devices such as are used in computer devices. In such devices, particularly in the so-called magnetic bubble, digital information is stored magnetically in the form of magnetic domains. The magnetically stored information can be read in various ways. It is now proposed that an optical beam upon interacting with the magnetic domain as will be explained has its polarization direction rotated by a certain angle. The magnitude of the optical readout signal is found to be proportional to the angle of rotation and can be made sufficiently cumulative as to be useful. In ordinary magneto-optic readout devices the interaction between the laser beam and the magnetic beam is rather limited and as a result the amount of rotation is small and is difficult to bring to useful application. The waveguide structure and optical readout proposed in the present invention enhances the interaction by making the polarization rotation effect cumulative through multiple reflections of the laser beam inside the waveguide layer. Thus, referring to FIG. 19 there is shown a schematic cross-sectional view of one of possible schemes for utilizing the present invention for an optical readout for a magnetic memory. Thus there is shown a means forming an optical waveguide layer 70 having a predetermined index of refraction n.sub.f. Means is provided for forming upper and lower boundary layers 72,74 disposed in intimate contact above and below or on each side of the waveguide layer 70. These boundary layers respectively have indices of refraction n.sub.d, n.sub.m. Each of the indices of refraction are related by n.sub.f >n.sub.d, n.sub.m and this condition defines the waveguide layer in such a way that propagation of energy therein can take place.

Coupling means is provided for introducing a beam of light through the upper layer into the waveguide layer and for receiving said beam after it is propagated at a predetermined distance in the waveguide layer. Such means typically takes the form of a coupling prism 76 of the type that is known in integrated optics art. As shown, two laser beams 78 and 80, for example, enter the coupling prism and are refracted downwardly and coupled through the upper boundary layer into the waveguide layer and subsequently after multiple internal reflections it re-emerges as output beams 82,84 the purpose of which will be explained hereinafter. In the examples shown, the waveguide layer is contained or is made of magnetic material having two alternate states having inducible magnetic flux orientation.

Means are provided for inducing the magnetic material to shift from one state to another and so positioned on the opposite side from the coupler. Such means can vary depending upon the nature of the magnetic material as shown in FIG. 19. One can obtain selective domain magnetization reversal in the x direction utilizing an overlay structure conductive grid 85 which guides the movement of the magnetized domains, commonly called magnetic bubbles. Typical material useful for carrying out construction of the layer 70 of the structure of FIG. 19 is

which has a magnetic anisotropy resulting in magnetization and reversal states oriented in the x direction. This material is sandwiched between boundary layers 72, 74 made of Eu.sub.x Er.sub.3-x Ga.sub.5 O.sub.12 or Gd.sub.3 Ga.sub.5 O.sub.12, for example. In the case of FIG. 20 the magnetization lies in the plane of waveguide, i.e. the z direction and can be excited by utilizing a moveable magnetic head 86 having a gap 88 positioned directly below the domain of interest. Such results in domain magnetization reversal in the z direction. In this case, the magnetic material has low magnetic anisotropic resulting in the z-directional control. Example of such a waveguide material in FIG. 20 includes the ferrites, (MFe.sub.2 O.sub.4) or garnets (M.sub.3 Fe.sub.5 O.sub.12) where M is a rare earth, and is deposited on a substrate material with matching lattice constant.

In connection with the example of FIG. 19 the resultant so-called magnetic bubbles can be moved under the influence of the overlay grid the current through which controls the location of the bubble. Consideration will now be given to the interaction of the light beam as it is coupled through the prism and interacts for a given length within the waveguide layer.

For evanescent modes, the decay constant p is of the order of the optical wave number k. Therefore, a film of 10 microns thick would behave as if it were a semi-infinite medium. If we choose the indices of refraction n.sub.m, n.sub.f, n.sub.d and n.sub.p for the various media in such a way that n.sub.p, n.sub.f >n.sub.d, n.sub.m, then a light beam can be coupled into and out of the optical waveguide. As compared to the conventional read-out device using Kerr effect, a significant improvement in the mode conversion efficiency is expected through multiple total reflections in the optical waveguide.

If the digital information is written by a magnetic head, then the domain magnetization will be either in the +z or -z direction and we have the longitudinal case. If the digital information is stored in the form of magnetic bubbles, then the domain magnetization will be either in the +x or -x direction and we have the polar case. In either the longitudinal or the polar case, the signs of r.sub.EM and r.sub.ME change as the magnetization direction reverses. As a result, the relative phase of TE and TM modes changes, and the polarization of the outgoing wave is rotated by an angle either +.theta..sub.F or -.theta..sub.F with respect to the incoming wave.

There are three basic factors which we must consider in comparing the performance of optical read-out devices. The first factor is the interaction length. Although it is possible in theory to achieve a complete mode conversion, the amount of conversion is limited in practice by the size of a magnetic domain. Taking z=50 microns to be the domain size, W=0.5 micron for the film thickness, and .theta.=55.degree., we have m=36 reflections. Note that the curves in FIG. 10 for different .theta..phi. almost coincide at small values to m. As long as m.delta..phi.<.tau., converted fields from successive reflections are almost in phase and hence additive. Therefore for small m, we may use the results calculated under the phase-matched condition. From FIG. 9, we find B.sub.m =0.5 if K.sub.12. Thus, in GdIG.sup.36 where K.sub.12 =2.times.10.sup.-3, we would expect a value of 0.1 for B.sub.m or a rotation of 5.8.degree. of the polarization.

Previously, multi-layer structures .sup.37, 38 have been proposed to enhance the read-out signal. In the previous proposals, the magnetic film is a part of the propagation medium. Significant improvement in the read-out signal has been reported..sup.39 However, the experiment was performed on magnetic films with uniform magnetization instead of domain magnetization. It is conceivable, therefore, that the reported improvement could be significantly reduced if the magnetic film were operated under small-size domains. The ratio of the domain size to the total thickness of the layered structure is the key factor because it determines the maximum number of multiple reflections.

The other two considerations regarding a read-out device are the loss in the magnetic medium and the adjustment for maximum read-out signal. In the structure we propose, the magnetic film serves as the waveguide wall where the wave is in the evanescent mode. Therefore, a magnetic substance with moderate loss may still be usable as a substrate material in cases where the ratio of the Verdet constant to the loss constant is reasonably large. It is expected that we have a wider range in the choice of magnetic materials for the waveguide case than for the layered-structure case. Second, the wave at the waveguide boundary is totally reflected once .theta.>.theta..sub.crit. For optimum performance of the layered structure, matching of impedance between the different dielectric layers is important. On the other hand, for small m, the phase-matching condition in the waveguide case becomes less critical. It seems that the waveguide structure offers some definite advantages over the multi-layer structure.

Reference is made to the previous polarization converter analysis and particularly to that portion of the analysis indicating conversion of TE or TM to the respective TM or TE polarizations. It is noted that the off diagonal elements appear in the magneto-optic case where there exists a non-vanishing magnetization. In accordance, K.sub.ij =iA.sub.ij pMp. When the sign of the induced magnetization changes the associated permittivity tensor component also changes. This results in switching the partially converted wave direction. When added to the unconverted portion, the result is a combined vector which shifts in orientation a certain amount upon each reflection and in accordance with the present invention, these shifts are cumulative. By using a polarizer arranged to discriminate completely at the output against one of the orientations the cumulative shift in output orientation can be detected by zeroing the output of the polarizer for that orientation corresponding to one state of magnetization of the material and then looking for a finite output whenever the material is induced into the second state. Thus, the output will be two valued, either zero or finite depending upon the direction of magnetization within the domain of the magnetic material. This shift is cumulative and the resultant orientations of the various vectors are indicated by the angle .theta..sub.f in the FIG. 19.

The use of gyrotropic and anisotropic material for mode conversion and signal control in thin-film optical waveguide application has been disclosed. A complete conversion from the TE to TM mode, or vice versa, is possible if the phase-matching condition is satisfied for a simultaneous propagation of degenerate TE and TM modes. Various schemes to achieve complete mode conversion are proposed. The essential property of the scattering matrix is r.sub.EM 'r.sub.ME '=-.vertline.r.sub.12 .vertline..sup.2 and this condition for complete conversion is satisifed in all optical media with a Hermitian permittivity tensor. The orientation of the gyrotropic and anisotropic material is so chosen that the mode conversion device is operated in either the longitudinal or the polar configuration.

The mode converter presented herein is a continuous one from which we can have any desired ratio of the TE to TM amplitude. The other important properties of the mode converter can be summarized as follows. If the direction of wave propagation reverses, the mode conversion factor changes its sign for the longitudinal gyrotropic and anisotropic cases and remains the same for the polar gyrotropic and anisotropic cases. On the other hand, the mode conversion factor always changes its sine if the bias field or the crystal orientation reverses its direction. These features are very useful because with them one can properly divide and mix the two modes to obtain the desired device function.

The condition which is crucial to the realization of the various device functions is the phase-matching condition .phi..sup.E =.phi..sup.M. These proposed several schemes have been proposed to fulfill this condition. In these calculations, substrate and film materials familiar to workers in the field have been disclosed to show that the thin-film optical waveguide mode converters proposed here are practicable. Two important practical considerations are the tolerances in the adjustment of the incident angle and in the control of the film thickness. Optimum design can be achieved by choosing a proper refractive index for the film and a proper amount of birefringence for the substrate or the top layer.

The requirements imposed on the control of the incident angle and the film thickness can be met in practice if we use a substrate material having large off-diagonal elements, say 10.sup.-3 or larger. Ferroelectric crystals operated near the Curie temperature but in the paraelectric phase seem to be most suitable for use as the electrooptic substrate material. There are magnetic materials which have a relatively large Faraday rotation but, at the same time, an appreciable extinction coefficient (as a result of either high conductivity or moderate absorption). It is expected that a resonable loss constant may be tolerable in the substrate material because of the evanescent nature of the wave. Furthermore, one can use an intermediate layer between the substrate and the film to moderate a high loss, a large birefringence, or a high conversion in the substrate. It is useful, therefore, to extend the analysis to such structures.

For an arbitrary direction of propagation, an electromagnetic wave in a gyrotropic or anisotropic medium possesses all three field components E.sub.x, E.sub.y, and E.sub.z. We express the field components in terms of E.sub.y :

where .eta.'=.sqroot..epsilon..sub.f /.mu..sub.0, .epsilon..sub.f is the permittivity of the thin film, and .theta. is the incident angle (FIG. 3). The coefficients A, B, C, and D are evaluated in Appendix 2. Since there are two eigen modes [p.sub.1 and p.sub.2 in Eq. (5)] in the substrate, we use subscripts 1 and 2 to indicate the field quantities associated with each mode.

In the isotropic film, we have: (1) for the TE mode, E in the y direction and H in the x-y plane and (2) for the TM mode, H in the y direction and E in the x-z plane. Applying the boundary conditions to the field components (FIG. 3), in the order of E.sub.y, E.sub.z, H.sub.y, and H.sub.z, we obtain

The fields E.sub.E ', E.sub.M ', E.sub.1 , and E.sub.2 can be expressed in terms of E.sub.E and E.sub.M by inverting the matrix. Let ##STR2## The various matrix elements are

(a) The Longitudinal Configuration

The relations between the electric fields can be obtained from the wave equation. From Eq. (2), we have ##EQU28## Thus, we find ##EQU29## The step taken in Eq. (A6b) is based on an identity derived from Eq. (5a). This step greatly simplifies the calculations later on. The magnetic fields can be obtained from .gradient..times.E. Using the relation .beta.=n.sub.f sin .theta. and the identity ##EQU30## we find the various coefficients in Eq. (A1) to be ##EQU31##

The matrix elements given in Eq. (7) can be found by substituting Eq. (A8) in Eq. (A4) and by letting F=R/N where N=(K.sub.12 sin.theta.)/(n.sub.f cos.sup.2 .theta.). Furthermore, the following identity

is used in obtaining Eq. (7d) from Eq. (A4c). The step is taken so that the relation between r.sub.ME and r.sub.EM can be clearly seen. For completeness, the matrix elements for the evanescent waves are given below: ##EQU32##

(b) The Polar Configuration

The field components in the polar configuration are related through the following equation: ##EQU33## Using Eq. (A11) and the identity ##EQU34## we obtain the following coefficients ##EQU35##

The various matrix elements for the polar configuration can be obtained by following a procedure similar to the one used for the longitudinal configuration if proper modifications are made. First, the common factor is N=i K.sub.23.sup.* /(n.sub.f.sup.2 cos.sup.2 .theta.) in F=R/N. Second, the following identity

is used obtaining Eq. (7f) from Eq. (A4c). Finally, the quantity K.sub.23 should replace K.sub.12 in Eq. (A10).

Sylvester's theorem states that for a matrix u with distinct eigen values .lambda..sub.1, . . . , .lambda..sub.n, the matrix polynominal can be written as ##EQU36## where ##EQU37## and I is the identity matrix. From Eq. (22), the two eigen values are

where .phi..sub.1 =tan.sup.-1 (r.sub.12 /r.sub.11) and r.sub.11.sup.2 +r.sub.12.sup.2 =1 on account of energy conservation. Letting P(u)=u.sup.m, we find ##EQU38## Thus, Eq. (23a) becomes ##EQU39## A similar equation is also obtained for Eq. (23b).

Now we turn to the case where the synchronous condition of Eq. (17) is not satisfied. We let .phi..sup.Ml =q' 2.pi.+.delta..phi./2 and .phi..sup.E =q2.pi.-.delta..phi./2. Thus, r.sub.EE '=r.sub.11 exp(-i.delta..phi./2) and r.sub.MM '=r.sub.11 exp(i.delta..phi./2). We notice, however, that .alpha.+.alpha.'=(q+q')2.pi. is still valid. Therefore, the elements r.sub.EM ' and r.sub.ME ' remain unchanged. The two eigen values of the scattering matrix are

where r=r.sub.11 cos (.delta..phi./2) and r'=(1-r.sup.2).sup.1/2. Again, the condition r.sub.11.sup.2 +r.sub.12.sup.2 =1 is used in the calculation. The matrix elements of [r.sub.ij '].sup.m are thus found to be

which replaces Eq. (23). The positive sign is for A.sub.EE and B.sub.EM while the minus sign is for A.sub.MM and B.sub.ME. Since energy is conserved in the system, the following identities exist: .vertline.A.sub.EE .vertline..sup.2 +.vertline.B.sub.ME .vertline..sup.2 =1 and .vertline.A.sub.MM .vertline..sup.2 +.vertline.B.sub.EM .vertline..sup.2 =1. These relations can be obtained from Eq. (A22) by noting that sin.sup.2 .phi..sub.2 =1-r.sub.11.sup.2 cos.sup.2 (.delta..phi./2).

1. p. k. tien, R. Ulrich, and R. J. Martin, Appl. Phys. Letters, Vol. 14, p. 291 (1969).

2. M. L. Dakss, L. Kuhn, P. F. Heindrich, and B. A. Scott, Appl. Phys. Letters, Vol. 16, p. 523 (1970)

3. D. Hall, A. Yariv and E. Garmire, Appl. Phys. Letters, Vol. 17, p. 127 (1970).

4. J. H. Harris and R. Shubert, IEEE Trans. Microwave Theory Tech., Vol. MTT-17, p. 269 (1971).

5. H. Kogelnik and C. V. Shank, Appl. Phys. Letters, Vol. 18, p. 152 (1971).

6. S. E. Miller, Bell Sys. Tech. J., Vol. 48, p. 2059 (1969).

7. N. Bloembergen, P. S. Pershan, and L. R. Wilson, Phys. Rev., Vol. 120, p. 2014 (1960).

8. I. P. Kaminov, Phys. Rev. Letters, Vol. 6, p. 528 (1961).

9. P. K. Tien and R. Ulrich, J. Opt. Soc. Am., Vol. 60, p. 1325 (1970).

10. D. O. Smith, Opt. Acta., Vol. 12, p. 13 and p. 195 (1965).

11. M. J. Freiser, IEEE Trans. on Magnetics, Vol. MAG-4, p. 152 (1968).

12. G. N. Ramachandran and S. Ramaseshan, "Crystal Optics," Encyclopedia of Physics (Handbuch der Physik), Vol. XXV/1, Springer, Berlin (1961).

13. L. D. Landau and E. M. Lifschitz, Electrodynamics of Continuous Media, Addison-Wesley Publishing Co., Reading, Mass. (1960).

14. P. S. Pershan, J. Appl. Phys., Vol. 38, p. 1482 (1967).

15. R. C. Le Craw, D. L. Wood, J. F. Dillon, Jr., and J. P. Remeika, Appl. Phys. Letters, Vol. 7, p. 27 (1965).

16. G. D. Boyd, R. C. Miller, K. Nassau, W. L. Bond, and A. Savage, Appl. Phys. Letters, Vol. 5, p. 234 (1964).

17. J. C. Suits, B. E. Argyle, M. J. Freiser, J. Appl. Phys., Vol. 37, p. 1391 (1966).

18. K. Y. Ahn and T. R. McGuire, J. Appl. Phys., Vol. 39, p. 5061 (1968).

19. P. V. Lenzo, E. G. Spencer, and A. A. Ballman, Appl. Phys. Letters, Vol. 11, p. 23 (1967).

20. J. E. Geusic S. K. Kurtz, L. G. Van Uitert, and S. H. Wemple, Appl. Phys. Letters, Vol. 4, p. 141 (1964).

21. M. Born and E. Wolf, Principle of Optics, p. 46, Pergamon Press, London (1959).

22. R. A. Frazer, W. J. Duncan, A. R. Collar, Elementary Matrices, p. 78, Cambridge University Press, Cambridge (1952).

23. Corning Glass Works Literature.

24. American Institute of Physics Handbook, McGraw-Hill Book Co., (1957) p. 6-24.

25. A. Smakula, Einkristalle, Springer, Berlin (1962).

26. C. D. Mee, Contemp. Phys., Vol. 8, p. 385 (1967).

27. T. M. Bleniewski and S. J. Czyzak, J. Opt. Soc. Am., Vol. 53, p. 496 (1963).

28. D. H. Hensler, J. D. Cuthbert, R. J. Martin, and P. K. Tien, Appl. Opt., Vol. 10, p. 1037 (1971).

29. R. Wolfe, A. J. Kurtzig, and R. C. LeCraw, J. Appl. Phys., Vol. 41, p. 1218 (1970).

30. L. W. Tilton, E. K. Plyler, R. E. Stephens, J. Opt. Soc. Am., Vol. 40, p. 540 (1950).

31. D. A. Pinnow, L. G. Van Uitert, A. N. Warner, and W. A. Bonner, Appl. Phys. Letters, Vol. 15, p. 83 (1969).

32. S. Wang, J. D. Crow, and M. Shah, to be published in Appl. Phys. Letters (Sept. 15 issue).

33. C. L. Hogan, Rev. Mod. Phys., Vol. 25, p. 253 (1953).

34. L. J. Aplet and J. W. Carson, Appl. Optics, Vol. 3, p. 544 (1964).

35. J. E. Geusic and H. E. D. Scovil, Bell Sys. Tech. J., Vol. 41, p. 1371 (1962).

36. J. T. Chang, J. F. Dillon, Jr., and U. F. Gianola, J. Appl. Phys., Vol. 30, p. 1110 (1965).

37. J. Franz and A. Schauer, Optik, Vol. 18, p. 186 (1961).

38. D. O. Smith, J. Appl. Phys., Vol. 35, p. 772 (1964).

39. A. S. Hoffman, W. R. Hagemeier, and S. H. Cushner, J. Appl. Phys., Vol. 41, p. 1407 (1970).