# Patent

## Title
System and method for optical fiber

## Abstract

In one embodiment, an optical fiber includes a first layer having a first refractive index and a second layer surrounding the first layer, where the second layer has a second refractive index, an inner radius, and an outer radius. The optical fiber also includes a third layer surrounding the second layer, where the third layer has a third refractive index, where the first refractive index is less than the second refractive index, where the third refractive index is less than the second refractive index, and where a ratio of the outer radius to the inner radius is less than 1.5.

## Background
This application claims the benefit of U.S. Provisional Application Ser. No. 61/890,669 filed on Oct. 14, 2013 and entitled “High Contrast Optical Fiber and Methods of Making and Using Thereof,” which application is hereby incorporated herein by reference.


TECHNICAL FIELD
The present invention relates to a system and method for optical communications, and, in particular, to a system and method for an optical fiber.
BACKGROUND
The capacity of optical fiber links is increasing. Optical fibers are used in optical fiber communications, datacenters, metropolitan area networks (MAN), access networks, wireless networks, and other networks. Improvements in optical fibers include low-loss single-mode fibers (SMFs), erbium doped fiber amplifiers (EDFAs), wavelength division multiplexing (WDM), and increased spectral efficiency using coherent communications. In optical fiber communications, optical signals are multiplexed in time, wavelength polarization, and phase. Optical signals may also be multiplexed in space by using multiple core fibers and/or higher mode fibers. Multicore fibers have a high complexity in fabrication and may be problematic to couple light into. Few mode fibers (FMF) and multi-mode fibers (MMFs) may also be used to multiplex optical data in space.
SUMMARY
An embodiment optical fiber includes a first layer having a first refractive index and a second layer surrounding the first layer, where the second layer has a second refractive index, an inner radius, and an outer radius. The optical fiber also includes a third layer surrounding the second layer, where the third layer has a third refractive index, where the first refractive index is less than the second refractive index, where the third refractive index is less than the second refractive index, and where a ratio of the outer radius to the inner radius is less than 1.5.
Another embodiment optical fiber includes a first layer having a first refractive index and a second layer surrounding the first layer, where the second layer has a second refractive index. The optical fiber also includes a third layer surrounding the second layer, where the third layer has a third refractive index and a fourth layer surrounding the third layer, where the fourth layer has a fourth refractive index, where the first refractive index is less than the second refractive index, where the third refractive index is less than the second refractive index, and where the third refractive index is less than the fourth refractive index.
An embodiment method of transmitting light includes generating a plurality of orbital angular momentum (OAM) modes and coupling the plurality of OAM modes into an optical fiber. The optical fiber includes a first layer having a first refractive index and a second layer surrounding the first layer, where the second layer has a second refractive index, an inner radius, and an outer radius, and a third layer surrounding the second layer, where the third layer has a third refractive index, where the first refractive index is less than the second refractive index, where the third refractive index is less than the second refractive index, and where a ratio of the outer radius to the inner radius is less than 1.5. The method also includes propagating the plurality of OAM modes along the optical fiber to produce a plurality of propagated OAM modes.
The foregoing has outlined rather broadly the features of an embodiment of the present invention in order that the detailed description of the invention that follows may be better understood. Additional features and advantages of embodiments of the invention will be described hereinafter, which form the subject of the claims of the invention. It should be appreciated by those skilled in the art that the conception and specific embodiments disclosed may be readily utilized as a basis for modifying or designing other structures or processes for carrying out the same purposes of the present invention. It should also be realized by those skilled in the art that such equivalent constructions do not depart from the spirit and scope of the invention as set forth in the appended claims.



BRIEF DESCRIPTION OF THE DRAWINGS
For a more complete understanding of the present invention, and the advantages thereof, reference is now made to the following descriptions taken in conjunction with the accompanying drawing, in which:
FIG. 1 illustrates Hermite-Gaussian mode distributions;
FIG. 2 illustrates Laguerre-Gaussian mode distributions;
FIG. 3 illustrates a graph change in refractive index and mode intensity verses position for a vortex fiber;
FIG. 4 illustrates a graph of effective index versus wavelength for modes in a vortex fiber;
FIG. 5 illustrates an image of an air core fiber and a graph of refractive index verses position for an air core fiber;
FIG. 6 illustrates a graph of effective refractive index verses wavelength for an air core fiber;
FIG. 7 illustrates generation of a computer generated hologram (CGH);
FIG. 8 illustrates generation of an orbital angular momentum (OAM) mode from a Gaussian mode using a CGH;
FIG. 9 illustrates OAM modes generated from Gaussian modes by a CGH;
FIG. 10 illustrates generation of OAM modes from Gaussian modes using a CGH;
FIG. 11 illustrates an OAM mode generated from a Gaussian mode by an SPP;
FIG. 12 illustrates OAM mode converting OAM modes to Gaussian modes using a CGH;
FIG. 13 illustrates an embodiment optical fiber profile;
FIG. 14 illustrates a graph of effective index versus wavelength for modes in an embodiment optical fiber;
FIG. 15 illustrates the profile of another embodiment optical fiber;
FIG. 16 illustrates a graph of effective index versus wavelength for modes in an embodiment optical fiber;
FIG. 17 illustrates a graph of fiber profile for embodiment optical fibers;
FIG. 18 illustrates an embodiment optical fiber;
FIG. 19 illustrates a system for testing optical fibers;
FIG. 20 illustrates determination of numerical aperture for an optical fiber;
FIG. 21 illustrates determination of numerical aperture for an embodiment optical fiber;
FIG. 22 illustrates measured refractive indices for an embodiment optical fiber;
FIG. 23 illustrates a graph of refractive indices for embodiment optical fibers;
FIG. 24 illustrates a graph of effective index versus wavelength of modes in an embodiment optical fiber;
FIG. 25 illustrates a graph of group index versus wavelength for modes in an embodiment optical fiber;
FIG. 26 illustrates a graph of dispersion verses wavelength for modes in an embodiment optical fiber;
FIG. 27 illustrates a flowchart of an embodiment method of propagating light in an optical fiber; and
FIG. 28 illustrates a block diagram of an embodiment computer system.



Corresponding numerals and symbols in the different figures generally refer to corresponding parts unless otherwise indicated. The figures are drawn to clearly illustrate the relevant aspects of the embodiments and are not necessarily drawn to scale.
DETAILED DESCRIPTION OF ILLUSTRATIVE EMBODIMENTS
It should be understood at the outset that although an illustrative implementation of one or more embodiments are provided below, the disclosed systems and/or methods may be implemented using any number of techniques, whether currently known or in existence. The disclosure should in no way be limited to the illustrative implementations, drawings, and techniques illustrated below, including the exemplary designs and implementations illustrated and described herein, but may be modified within the scope of the appended claims along with their full scope of equivalents.
An embodiment optical fiber transmits multiple orbital angular momentum (OAM) modes with a separation between the modes. An embodiment optical fiber has a large hollow core, surrounded by a ring-core, which is a thin layer of high refractive index material. Around the high refractive index material is another slightly wider layer of material with a low refractive index, surrounded by silica cladding. The high contrast between the hollow core and high refractive index achieve a significant difference between the effective indices of OAM modes in the optical fiber.
Multimode fibers may transmit multiple linearly polarized (LP) modes. However, a system with LP modes uses heavy multiple-input multiple-output (MIMO) processing, because LP modes contain hybrid eigenmodes with slightly different propagation constants, which may cause mode walk-off during propagation.
Spin angular momentum (SAM) of light causes light to be polarized. The SAM is associated with the spin of photons, and has a value of σħ. Also, light may spin and carry OAM. OAM is a natural property of helically phased beams. Helically phased beams are described by electric and magnetic fields with an azimuthal phase in the form:

u(r,φ,z)=u0(r,z)ejlφ,

where u is the envelope of the field and l is an integer. Because the phase is not defined at r=0, there is no intensity at the center of the beam. Thus, OAM beams are known as optical vortex or cylindrical vector beams (CVBs). Under the paraxial approximation, total angular momentum is the sum of spin and orbital angular momentum (j=σ+l), implying that spin and orbital angular momentum are independent.

OAM light has a variety of applications, including the ability to spin microscopic objects, trapping and guiding atoms, in imaging, and in quantum optics. OAM states may carry information, because OAM states are orthogonal to each other. OAM states may be transmitted over back-to-back systems, without a propagation distance, to demonstrate multiplexing and demultiplexing. Other examples are over C-Band and combining OAM states with wavelength division multiplexing (WDM), polarization division multiplexing (PDM), and coherent detection. Also, OAM states may be used to perform quantum cryptography.
Optical fibers may generate OAM modes using mechanically stressed fibers, acousto-optic interaction, fiber couplers, few-mode fibers (FMFs), and highly non-linear photonic crystal fibers, using fibers a few centimeters long, not used as propagation media. In one example, OAM modes are propagated into multimode fiber.
Optical modes exist in guiding media, such as optical fibers, as well as in free-space laser beams. Transverse electric (TE) fields have no magnetic field in the direction of propagation and a nonzero electric field in the direction of propagation, while transverse magnetic (TM) fields have a nonzero magnetic field in the direction of propagation and no electric field in the direction of propagation. Transverse electro-magnetic (TEM) fields have no electric or magnetic fields in the direction of propagation, while hybrid fields have nonzero electrical and magnetic fields in the direction of propagation.
Many laser beams may be approximated by a Gaussian function in a Gaussian mode. A Gaussian beam is described by:




  
    
      
        E
        →
      
      ⁡
      
        (
        
          r
          ,
          z
        
        )
      
    
    =
    
      
        E
        0
      
      ⁢
      
        
          w
          0
        
        
          w
          ⁡
          
            (
            z
            )
          
        
      
      ⁢
      
        exp
        ⁡
        
          (
          
            -
            
              
                r
                2
              
              
                
                  w
                  ⁡
                  
                    (
                    z
                    )
                  
                
                2
              
            
          
          )
        
      
      ⁢
      
        exp
        ⁡
        
          (
          
            -
            
              j
              ⁡
              
                [
                
                  kz
                  -
                  
                    arctan
                    ⁢
                    
                      z
                      
                        z
                        R
                      
                    
                  
                  +
                  
                    
                      kr
                      2
                    
                    
                      2
                      ⁢
                      
                        
                      
                      ⁢
                      
                        R
                        ⁡
                        
                          (
                          z
                          )
                        
                      
                    
                  
                
                ]
              
            
          
          )
        
      
    
  
  ,




where |E0| is the peak amplitude, ω0 is the beam radius at the beam waist, k=2π/λ is the wavenumber, the beam radius is given by:





  
    
      ω
      ⁡
      
        (
        z
        )
      
    
    =
    
      
        ω
        0
      
      ⁢
      
        
          1
          +
          
            
              (
              
                z
                
                  z
                  R
                
              
              )
            
            2
          
        
      
    
  
  ,




zR=2πω/λ is the Rayleigh length, and R(z) is the radius of the wavefront.

Many laser cavities are rectangular, and thus produce rectangular shaped Hermite-Gaussian modes, which are identified as HGn,m, where n+1 and m+1 are the number of rows and columns of intensity spots. HG0,0 is the fundamental Gaussian mode. Hermite-Gaussian modes are described by the product of a Gaussian function and a Hermite polynomial. The intensity of the electric field is given by:




  
    
      
        
          E
          →
        
        
          n
          ,
          m
        
      
      ⁡
      
        (
        
          x
          ,
          y
          ,
          z
        
        )
      
    
    =
    
      
        E
        0
      
      ⁢
      
        
          w
          0
        
        
          w
          ⁡
          
            (
            z
            )
          
        
      
      ⁢
      
        
          H
          n
        
        ⁡
        
          (
          
            
              2
            
            ⁢
            
              x
              
                w
                ⁡
                
                  (
                  z
                  )
                
              
            
          
          )
        
      
      ⁢
      
        exp
        ⁡
        
          (
          
            -
            
              
                x
                2
              
              
                
                  w
                  ⁡
                  
                    (
                    z
                    )
                  
                
                2
              
            
          
          )
        
      
      ⁢
      
        
          H
          m
        
        ⁡
        
          (
          
            
              2
            
            ⁢
            
              y
              
                w
                ⁡
                
                  (
                  z
                  )
                
              
            
          
          )
        
      
      ⁢
      
        exp
        ⁡
        
          (
          
            -
            
              
                x
                2
              
              
                
                  w
                  ⁡
                  
                    (
                    z
                    )
                  
                
                2
              
            
          
          )
        
      
      ⁢
      
        exp
        ⁡
        
          (
          
            -
            
              j
              ⁡
              
                [
                
                  kz
                  -
                  
                    
                      (
                      
                        1
                        +
                        n
                        +
                        m
                      
                      )
                    
                    ⁢
                    arctan
                    ⁢
                    
                      z
                      
                        z
                        R
                      
                    
                  
                  +
                  
                    
                      k
                      ⁡
                      
                        (
                        
                          
                            x
                            2
                          
                          +
                          
                            y
                            2
                          
                        
                        )
                      
                    
                    
                      2
                      ⁢
                      
                        
                      
                      ⁢
                      
                        R
                        ⁡
                        
                          (
                          z
                          )
                        
                      
                    
                  
                
                ]
              
            
          
          )
        
      
    
  
  ,




where Hn(x) is the Hermite polynomial, with non-negative integer index n, defined by:





  
    
      H
      n
    
    ⁡
    
      (
      z
      )
    
  
  =
  
    
      
        (
        
          -
          1
        
        )
      
      n
    
    ⁢
    
      e
      
        x
        2
      
    
    ⁢
    
      
        d
        n
      
      
        d
        
          x
          n
        
      
    
    ⁢
    
      
        e
        
          -
          
            x
            2
          
        
      
      .
    
  




FIG. 1 illustrates Hermite-Gaussian mode, where mode 782 is HG0,0, mode 784 is HG0,1, mode 786 is HG0,2, mode 788 is HG0,3, mode 790 is HG1,0, mode 792 is HG1,1, mode 794 is HG1,2, and mode 796 is HG1,3.

Laguerre-Gaussian modes describe cylindrical laser modes. Laguerre-Gaussian modes are denoted by LGn,m, where l=|n−m| is the azimuthal index and p=min(n, m) is the radial index. The field diameter depends on l, while the number of intensity rings is given by p+1. The fundamental mode is LG0,0. The electric field is described by:




  
    
      
        
          E
          →
        
        
          n
          ,
          m
        
      
      ⁡
      
        (
        
          r
          ,
          φ
          ,
          z
        
        )
      
    
    =
    
      
        E
        0
      
      ⁢
      
        
          w
          0
        
        
          w
          ⁡
          
            (
            z
            )
          
        
      
      ⁢
      
        
          L
          
            min
            ⁡
            
              (
              
                n
                ,
                m
              
              )
            
          
          
            
            
              n
              -
              m
            
            
          
        
        ⁡
        
          (
          
            
              2
              ⁢
              
                
              
              ⁢
              
                r
                2
              
            
            
              
                w
                ⁡
                
                  (
                  r
                  )
                
              
              2
            
          
          )
        
      
      ⁢
      
        
          (
          
            -
            1
          
          )
        
        
          min
          ⁡
          
            (
            
              n
              ,
              m
            
            )
          
        
      
      ⁢
      
        
          (
          
            
              2
            
            ⁢
            
              r
              
                w
                ⁡
                
                  (
                  z
                  )
                
              
            
          
          )
        
        
          
          
            n
            -
            m
          
          
        
      
      ⁢
      
        exp
        ⁡
        
          (
          
            -
            
              
                r
                2
              
              
                
                  w
                  ⁡
                  
                    (
                    z
                    )
                  
                
                2
              
            
          
          )
        
      
      ⁢
      
        exp
        ⁡
        
          (
          
            -
            
              j
              ⁡
              
                [
                
                  
                    
                      (
                      
                        n
                        -
                        m
                      
                      )
                    
                    ⁢
                    φ
                  
                  +
                  
                    
                      (
                      
                        1
                        +
                        n
                        +
                        m
                      
                      )
                    
                    ⁢
                    arctan
                    ⁢
                    
                      z
                      
                        z
                        R
                      
                    
                  
                  +
                  
                    
                      kr
                      2
                    
                    
                      2
                      ⁢
                      
                        
                      
                      ⁢
                      
                        R
                        ⁡
                        
                          (
                          z
                          )
                        
                      
                    
                  
                
                ]
              
            
          
          )
        
      
    
  
  ,




where Lpl is the generalized Laguerre polynomial, defined as:





  
    
      L
      p
      l
    
    ⁡
    
      (
      x
      )
    
  
  =
  
    
      
        
          x
          
            -
            l
          
        
        ⁢
        
          e
          x
        
      
      
        p
        !
      
    
    ⁢
    
      
        d
        p
      
      
        d
        
          x
          p
        
      
    
    ⁢
    
      
        (
        
          
            e
            
              -
              x
            
          
          ⁢
          
            x
            
              p
              +
              l
            
          
        
        )
      
      .
    
  




LGn,m modes may be built from linear combinations of HGn,m modes, for example:




  
    
      
        LG
        
          0
          ,
          2
        
      
      =
      
        
          
            1
            2
          
          ⁢
          
            HG
            
              2
              ,
              0
            
          
        
        +
        
          
            j
            
              2
            
          
          ⁢
          
            HG
            
              1
              ,
              1
            
          
        
        -
        
          
            1
            2
          
          ⁢
          
            HG
            
              0
              ,
              2
            
          
        
      
    
    ,
    
      

    
    ⁢
    
      
        LG
        
          1
          ,
          1
        
      
      =
      
        
          
            1
            
              2
            
          
          ⁢
          
            HG
            
              2
              ,
              0
            
          
        
        +
        
          
            1
            
              2
            
          
          ⁢
          
            HG
            
              0
              ,
              2
            
          
        
      
    
    ,
    
      

    
    ⁢
    and
  




  
    
      LG
      
        2
        ,
        0
      
    
    =
    
      
        
          1
          2
        
        ⁢
        
          HG
          
            2
            ,
            0
          
        
      
      -
      
        
          j
          
            2
          
        
        ⁢
        
          HG
          
            1
            ,
            1
          
        
      
      -
      
        
          1
          2
        
        ⁢
        
          
            HG
            
              0
              ,
              2
            
          
          .
        
      
    
  



FIG. 2 illustrates some LG modes, including mode 802, which is LG0,0, mode 804, which is LG0,1, mode 806, which is LG0,2, mode 808, which is LG0,3, mode 810, which is LG1,0, mode 812, which is LG1,1, mode 814, which is LG1,2, and mode 816, which is LG1,3.

Vector modes are the eigenmodes of cylindrical optical fibers. The electrical and magnetic fields of vector modes have nonzero components in all directions on the Cartesian plane. Vector modes are represented by EH when the electrical field dominates in the direction of propagation and HE when the magnetic field dominates in the direction of propagation. Vector modes are identified by the parameters v and m, where v is a non-negative integer corresponding to the number of azimuthal symmetry axis in the field components and m is a positive integer corresponding to the number of concentric rings in the intensity of the field. Also, vector modes may be even or odd, corresponding to the azimuthal symmetry of the fields. Modes are identified by ξv,mψ where ξ is HE, EH, TE, or TM, and ψ is either even or odd. When v=0, vector modes are not hybrid, but are either TE or TM modes. Solving the characteristic equation of a waveguide finds the supported molds. The value of v and the possible effective indices may be determined. Modes are ordered based on their effective indices, starting with the highest effective index, to yield:

HEv,1>EHv,1>HEv,2>EHv,2>HEv,3> . . .

TE0,1>TE0,2>TE0,3> . . .

TM0,1>TM0,2>TM0,3> . . .

In many cases, a weakly guided estimation may be used, when the indices of the core (n1) and the cladding (n2) are nearly identical. Fields of LP modes only have one component, either in the x or y axis. LP modes may be referred to as scalar modes. LP modes are TEM modes, because they have zero electrical and magnetic field in the direction of propagation. LP modes may have the form LPl,m polψ, where ψ is even or odd, l is a non-negative integer corresponding to the number of azimuthal symmetry axes in the fields, m is a positive integer corresponding to the number of concentric rings in the intensity fields, and pol is x or y, indicating the axis of polarization of the electrical field. When l=0, the intensity of the fields is circular, and there is no even and odd degeneracies. For other LP modes there are four variants, even or odd, and x or y. LP modes are combinations of vector modes with similar effective indices given by:




  
    LP
    
      l
      ,
      m
    
  
  =
  
    {
    
      
        
          
            
              
                HE
                
                  1
                  ,
                  m
                
                even
              
              +
              
                HE
                
                  1
                  ,
                  m
                
                odd
              
            
          
          
            
              l
              =
              0
            
          
        
        
          
            
              
                HE
                
                  2
                  ,
                  m
                
                even
              
              +
              
                HE
                
                  2
                  ,
                  m
                
                odd
              
              +
              
                TM
                
                  0
                  ,
                  m
                
              
              +
              
                TE
                
                  0
                  ,
                  m
                
              
            
          
          
            
              l
              =
              1
            
          
        
        
          
            
              
                HE
                
                  
                    l
                    +
                    1
                  
                  ,
                  m
                
                even
              
              +
              
                HE
                
                  
                    l
                    +
                    1
                  
                  ,
                  m
                
                odd
              
              +
              
                EH
                
                  
                    l
                    -
                    1
                  
                  ,
                  m
                
                even
              
              +
              
                EH
                
                  
                    l
                    -
                    1
                  
                  ,
                  m
                
                odd
              
            
          
          
            
              l
              >
              1
            
          
        
      
      .
    
  




A vortex fiber may propagate OAM modes. Mode OAMp,l∓ is defined where + or − indicates whether the mode is right circularly polarized or left circularly polarized, respectively, the number of concentric rings of intensity in the field is p+1, the absolute value of l is the number of times the phase of the field goes from 0 to 2π over the azimuthal circle, and the sign of l identifies whether the rotation of the phase of the field is clockwise (right handed) or counterclockwise (left handed). In multimode fibers, the HE2,1, TM0,1, and TE0 modes have similar propagation constants, causing them to couple to the LP1,1 mode. The vortex fiber separates the effective indices of those modes, pulling them more than 10−4 apart, facilitating the HE2,1 mode, an eigenmode of OAM0,1 mode, to propagate without degeneracy. The vortex fiber may propagate 2 OAM and 2 LP0,1 modes using low density parity check (LDPC) codes and quadrature phase-shift keying (QPSK) modulation. Also, a vortex fiber may transmit WDM with 16 quadrature amplitude modulation (QAM). FIG. 3 illustrates a graph of refractive index profile (curve 252) verses position, and of mode intensity (curve 254). FIG. 4 illustrates a graph of effective index of LP1,1, eigenmodes, versus position. Curve 250 shows the effective index for TE0,1, curve 256 shows the effective index for HE2,1, and curve 258 shows the effective index for TM0,1. This fiber supports four modes, OAM0,0∓, which corresponds to circularly polarized LP0,1, and OAM0,1∓ modes. Mode separation is about 1.8·10−4. The center core supports a Gaussian LP0,1 mode, while the outer ring supports the OAM0,1 mode. This vortex fiber does not support higher order OAM modes.
An example air core fiber transmits 12 OAM modes over 2 m and 2 OAM modes over 2 km. The air core fiber supports a large number of modes, with a large average mode separation, and a few modes which have a larger separation. FIG. 5 illustrates air core fiber 306 with air core 310 and cladding 308, and graph 312 with the index of refraction verses radius for an air core fiber. FIG. 6 illustrates a graph of refractive index mode separation verses versus wavelength. Curve 314 shows the refractive index separation for l=9, curve 316 shows the refractive index separation for l=8, curve 318 show the refractive index separation for l=7, and curves 320 show the refractive index separation for lower order modes. Modes OAM0,∓9∓ and OAM0,∓8∓ have a mode separation greater than 1·10−4.
OAM modes may use less MIMO processing than LP modes for spatial multiplexing, because eigenmodes of OAM modes are two degeneracies of the same mode with the same propagation constant, while eigenmodes of LP modes are different vector modes with different propagation constants. Also, the intensity fields of OAM modes are closer to Gaussian fields than LP modes, suggesting better coupling efficiency between the laser and fiber modes. Additionally, OAM modes are polarization maintaining.
There are several methods of generating OAM modes. OAM beams have a phase singularity in their light field. The phase of the OAM beam varies from 0 to 2πl, where l is an integer, either clockwise or counterclockwise around the singularity. OAM beams may be generated by transforming the phase of a laser beam or a Gaussian beam. The phase transformation may be done through computer generated hologram (CGH) or spiral phase pattern (SPP).
Passive elements, such as CHG and SPP, are relatively inexpensive, but lack versatility, are precisely manufactured, are wavelength dependent, and have losses from the use of free-space optics. CGH may be printed on a transparency using a printer, etched on a glass plate using photolithography, or transferred to a spatial light modulator (SLM) through a computer to produce different order OAM beams. OAM beams may be generated using a hologram constructed on film or glass. FIG. 7 illustrates the generation of OAM beams using a CGH. A hologram is printed on reflective medium 602. The construction of the hologram is done by illuminating the film with a reference wave:

ejk×sin(θ), 

and with

ejlφ.

In another example, the hologram is made by programming a spatial light modulator. FIG. 8 illustrates a reference wave on hologram 612 producing mode 614, a donut shaped intensity signal on an OAM beam examined on a charge coupled device (CCD) camera. FIG. 9 illustrates the conversion of a Gaussian beam to an OAM beam using CGH. Gaussian beam spot 672 is converted to OAM beam spots 678 and 682 using CGH 674. The OAM beams are formed at the diffraction order and spatially dispersed off center to OAM beams 678 and 682, while non-diffracted light remains at beam 680. Due to the nature of the conversion process, the available power is reduced, because several diffraction orders are produced from the same light source.

In another example, OAM beams are generated by transferring an SPP via an SLM. The entire Gaussian beam may be converted into an OAM beam. FIG. 10 illustrates the generation of an OAM beam using a CGH. Laser 622 generates a light beam which propagates along single mode fiber (SMF) 624. The beam passes through lens 626 to become a free space beam, which reflects off SLM 628 and mirror 630. The beam is received by CCD camera 632 and analyzed on computer 634, where image 636 is observed. Using different areas of the same SLM, many OAM beams may be simultaneously generated. Also, the SLM is reconfigurable, but are expensive, and the coupling between optical fibers and free space optics may be a source of loss. FIG. 11 illustrates the generation of OAM beams using an SPP. Gaussian beam spot 662 is converted to OAM beam spot 666 using SPP 664.
Photonic integrated circuits (PICs) may be used. PICs facilitate coupling the OAM beams directly into the optical fiber without free space. In another example, a multicore fiber is used to generate and recover OAM modes. Alternatively, each channel is a separate ring of light into the fiber.
OAM beams may be converted to Gaussian beams, for example for a receiver designed to detect a Gaussian intensity shaped beam. When only one OAM mode is present, an SPP technique may be used. When multiple OAM modes are used, demultiplexing may simultaneously be performed. For example, a CGH technique creates spatial separation of the converted modes, and it may be used to convert and demultiplex simultaneously. FIG. 12 illustrates the conversion of an OAM to three modes. When l=−1, OAM mode 692 is converted to Gaussian modes 698, 700, and 702 by CGH 694. When l=1, OAM mode 692 is converted to Gaussian modes 706, 708, and 710 by CGH 694.
Modeling may precisely describe OAM modes in optical fibers. Fiber modes may be derived from Maxwell's equations, and may be solved for fibers with many steps or layers. For a non-magnetic dielectric optical fiber with time dependency and no variation on the z propagation axis:



  
    
      
        
          ∇
          2
        
        ⁢
        
          E
          →
        
      
      +
      
        ∇
        
          (
          
            
              E
              →
            
            ·
            
              
                ∇
                ɛ
              
              ɛ
            
          
          )
        
      
      +
      
        
          ω
          2
        
        ⁢
        
          μ
          0
        
        ⁢
        ɛ
        ⁢
        
          
        
        ⁢
        
          E
          →
        
      
    
    =
    0
  




  
    
      
        
          ∇
          2
        
        ⁢
        
          H
          →
        
      
      +
      
        
          (
          
            
              ∇
              ɛ
            
            ɛ
          
          )
        
        ×
        
          ∇
          
            ×
            
              H
              →
            
          
        
      
      +
      
        
          ω
          2
        
        ⁢
        
          μ
          0
        
        ⁢
        ɛ
        ⁢
        
          
        
        ⁢
        
          H
          →
        
      
    
    =
    0.
  



For radial position r and azimuthal position φ at time t:

{right arrow over (E)}={right arrow over (E)}(r,φ,z,t)={right arrow over (E)}(r,φ)exp{j(ωt−βz)},

and

{right arrow over (H)}={right arrow over (H)}(r,φ,z,t)={right arrow over (H)}(r,φ)exp{j(ωt−βz)},

are the electric and magnetic fields, respectively, ∇2 is the Laplacian, ∇ is the gradient, ∇× is the curl operator, ∈ is the dielectric permittivity of the material, μ0=1.257·10−6 H/m, the vacuum permeability, and ω is the angular frequency of the light. Also, the propagation constant of a given mode is defines as:

β=kneff,

where neff is the effective index of the mode, and:





  k
  =
  
    
      
        2
        ⁢
        π
      
      λ
    
    =
    
      
        ω
        c
      
      .
    
  




Assuming the fiber is a step index material, the equations may be simplified by removing the gradient, because ∇∈=0. This yields expressions of the electric field and magnetic field in the form of Bessel functions Jv, Yv, Kv, and Iv. The order of the Bessel function, v, is a non-negative integer. The fields may be separated into three components: ez and hz in the direction of propagation, er and hr radially in the fiber, and eφ and hφ perpendicular to the fiber. These components are linked such that:



  
    
      
        
          e
          r
        
        ⁡
        
          (
          
            r
            ,
            φ
          
          )
        
      
      =
      
        
          ±
          j
        
        ⁢
        
          
            
              k
              ⁢
              
                
              
              ⁢
              
                ρ
                2
              
            
            
              U
              2
            
          
          ⁡
          
            [
            
              
                
                  n
                  eff
                
                ⁢
                
                  
                    ∂
                    
                      e
                      z
                    
                  
                  
                    ∂
                    r
                  
                
              
              +
              
                
                  
                    
                      μ
                      0
                    
                    
                      ɛ
                      0
                    
                  
                
                ⁢
                
                  1
                  r
                
                ⁢
                
                  
                    ∂
                    
                      h
                      z
                    
                  
                  
                    ∂
                    φ
                  
                
              
            
            ]
          
        
      
    
    ,
    
      

    
    ⁢
    
      
        
          e
          φ
        
        ⁡
        
          (
          
            r
            ,
            φ
          
          )
        
      
      =
      
        
          ±
          j
        
        ⁢
        
          
            
              k
              ⁢
              
                
              
              ⁢
              
                ρ
                2
              
            
            
              U
              2
            
          
          ⁡
          
            [
            
              
                
                  
                    n
                    eff
                  
                  r
                
                ⁢
                
                  
                    ∂
                    
                      e
                      z
                    
                  
                  
                    ∂
                    φ
                  
                
              
              -
              
                
                  
                    
                      μ
                      0
                    
                    
                      ɛ
                      0
                    
                  
                
                ⁢
                
                  
                    ∂
                    
                      h
                      z
                    
                  
                  
                    ∂
                    r
                  
                
              
            
            ]
          
        
      
    
    ,
    
      

    
    ⁢
    
      
        
          h
          r
        
        ⁡
        
          (
          
            r
            ,
            φ
          
          )
        
      
      =
      
        
          ±
          j
        
        ⁢
        
          
            
              k
              ⁢
              
                
              
              ⁢
              
                ρ
                2
              
            
            
              U
              2
            
          
          ⁡
          
            [
            
              
                
                  n
                  eff
                
                ⁢
                
                  
                    ∂
                    
                      h
                      z
                    
                  
                  
                    ∂
                    r
                  
                
              
              -
              
                
                  
                    
                      ɛ
                      0
                    
                    
                      μ
                      0
                    
                  
                
                ⁢
                
                  
                    n
                    2
                  
                  r
                
                ⁢
                
                  
                    ∂
                    
                      e
                      z
                    
                  
                  
                    ∂
                    φ
                  
                
              
            
            ]
          
        
      
    
    ,
    
      

    
    ⁢
    and
  




  
    
      
        
          h
          φ
        
        ⁡
        
          (
          
            r
            ,
            φ
          
          )
        
      
      =
      
        
          ±
          j
        
        ⁢
        
          
            
              k
              ⁢
              
                
              
              ⁢
              
                ρ
                2
              
            
            
              U
              2
            
          
          ⁡
          
            [
            
              
                
                  
                    n
                    eff
                  
                  r
                
                ⁢
                
                  
                    ∂
                    
                      h
                      z
                    
                  
                  
                    ∂
                    φ
                  
                
              
              +
              
                
                  
                    
                      ɛ
                      0
                    
                    
                      μ
                      0
                    
                  
                
                ⁢
                
                  n
                  2
                
                ⁢
                
                  
                    ∂
                    
                      e
                      z
                    
                  
                  
                    ∂
                    r
                  
                
              
            
            ]
          
        
      
    
    ,
  



where ∈0=8.854·10−12 F/m is the vacuum permittivity and

U=kρ√{square root over (|neff2−n2|)}.

Also, ρ is the outer radius of the fiber layer at point (r, φ) (or the inner radius for the cladding), and n is the refractive index at point (r, φ). The sign of er, hr, eφ, and hφ depend on whether neff is smaller than n or greater than n.

Differential equations for ez and hz vary based on whether the refractive index n of the current layer is smaller or greater than the effective index neff. There are two sets of solutions with different Bessel functions, given by:




  
    
      e
      z
    
    ⁡
    
      (
      
        r
        ,
        φ
      
      )
    
  
  =
  
    {
    
      
        
          
            
              
                [
                
                  
                    A
                    ⁢
                    
                      
                        
                          J
                          v
                        
                        ⁡
                        
                          (
                          
                            Ur
                            /
                            ρ
                          
                          )
                        
                      
                      
                        
                          J
                          v
                        
                        ⁡
                        
                          (
                          U
                          )
                        
                      
                    
                  
                  +
                  
                    B
                    ⁢
                    
                      
                        
                          Y
                          v
                        
                        ⁡
                        
                          (
                          
                            Ur
                            /
                            ρ
                          
                          )
                        
                      
                      
                        
                          Y
                          v
                        
                        ⁡
                        
                          (
                          U
                          )
                        
                      
                    
                  
                
                ]
              
              ⁢
              
                
                  f
                  v
                
                ⁡
                
                  (
                  φ
                  )
                
              
            
          
          
            
              
                if
                ⁢
                
                  
                
                ⁢
                
                  n
                  eff
                
              
              <
              n
            
          
        
        
          
            
              
                [
                
                  
                    A
                    ⁢
                    
                      
                        
                          K
                          v
                        
                        ⁡
                        
                          (
                          
                            Ur
                            /
                            ρ
                          
                          )
                        
                      
                      
                        
                          K
                          v
                        
                        ⁡
                        
                          (
                          U
                          )
                        
                      
                    
                  
                  +
                  
                    B
                    ⁢
                    
                      
                        
                          I
                          v
                        
                        ⁡
                        
                          (
                          
                            Ur
                            /
                            ρ
                          
                          )
                        
                      
                      
                        
                          I
                          v
                        
                        ⁡
                        
                          (
                          U
                          )
                        
                      
                    
                  
                
                ]
              
              ⁢
              
                
                  f
                  v
                
                ⁡
                
                  (
                  φ
                  )
                
              
            
          
          
            
              
                if
                ⁢
                
                  
                
                ⁢
                
                  n
                  eff
                
              
              >
              n
            
          
        
      
      ,
      
        

      
      ⁢
      
        
          and
          ⁢
          
            

          
          ⁢
          
            
              h
              z
            
            ⁡
            
              (
              
                r
                ,
                φ
              
              )
            
          
        
        =
        
          {
          
            
              
                
                  
                    
                      [
                      
                        
                          
                            A
                            ′
                          
                          ⁢
                          
                            
                              
                                J
                                v
                              
                              ⁡
                              
                                (
                                
                                  Ur
                                  /
                                  ρ
                                
                                )
                              
                            
                            
                              
                                J
                                v
                              
                              ⁡
                              
                                (
                                U
                                )
                              
                            
                          
                        
                        +
                        
                          
                            B
                            ′
                          
                          ⁢
                          
                            
                              
                                Y
                                v
                              
                              ⁡
                              
                                (
                                
                                  Ur
                                  /
                                  ρ
                                
                                )
                              
                            
                            
                              
                                Y
                                v
                              
                              ⁡
                              
                                (
                                U
                                )
                              
                            
                          
                        
                      
                      ]
                    
                    ⁢
                    
                      
                        g
                        v
                      
                      ⁡
                      
                        (
                        φ
                        )
                      
                    
                  
                
                
                  
                    
                      if
                      ⁢
                      
                        
                      
                      ⁢
                      
                        n
                        eff
                      
                    
                    <
                    n
                  
                
              
              
                
                  
                    
                      [
                      
                        
                          
                            A
                            ′
                          
                          ⁢
                          
                            
                              
                                K
                                v
                              
                              ⁡
                              
                                (
                                
                                  Ur
                                  /
                                  ρ
                                
                                )
                              
                            
                            
                              
                                K
                                v
                              
                              ⁡
                              
                                (
                                U
                                )
                              
                            
                          
                        
                        +
                        
                          
                            B
                            ′
                          
                          ⁢
                          
                            
                              
                                I
                                v
                              
                              ⁡
                              
                                (
                                
                                  Ur
                                  /
                                  ρ
                                
                                )
                              
                            
                            
                              
                                I
                                v
                              
                              ⁡
                              
                                (
                                U
                                )
                              
                            
                          
                        
                      
                      ]
                    
                    ⁢
                    
                      
                        g
                        v
                      
                      ⁡
                      
                        (
                        φ
                        )
                      
                    
                  
                
                
                  
                    
                      if
                      ⁢
                      
                        
                      
                      ⁢
                      
                        n
                        eff
                      
                    
                    >
                    n
                  
                
              
            
            ,
          
        
      
    
  




where A, B, A′, and B′ are constants.





  
    
      f
      v
    
    ⁡
    
      (
      φ
      )
    
  
  =
  
    {
    
      
        
          
            
              cos
              ⁢
              
                
              
              ⁢
              
                (
                
                  v
                  ⁢
                  
                    
                  
                  ⁢
                  φ
                
                )
              
            
          
          
            
              even
              ⁢
              
                
              
              ⁢
              mode
            
          
        
        
          
            
              sin
              ⁢
              
                
              
              ⁢
              
                (
                
                  v
                  ⁢
                  
                    
                  
                  ⁢
                  φ
                
                )
              
            
          
          
            
              odd
              ⁢
              
                
              
              ⁢
              mode
            
          
        
      
      ,
      
        

      
      ⁢
      
        
          and
          ⁢
          
            

          
          ⁢
          
            
              g
              v
            
            ⁡
            
              (
              φ
              )
            
          
        
        =
        
          {
          
            
              
                
                  
                    
                      -
                      sin
                    
                    ⁢
                    
                      
                    
                    ⁢
                    
                      (
                      
                        v
                        ⁢
                        
                          
                        
                        ⁢
                        φ
                      
                      )
                    
                  
                
                
                  
                    even
                    ⁢
                    
                      
                    
                    ⁢
                    mode
                  
                
              
              
                
                  
                    cos
                    ⁢
                    
                      
                    
                    ⁢
                    
                      (
                      
                        v
                        ⁢
                        
                          
                        
                        ⁢
                        φ
                      
                      )
                    
                  
                
                
                  
                    odd
                    ⁢
                    
                      
                    
                    ⁢
                    mode
                  
                
              
            
            .
          
        
      
    
  




Even and odd modes are naming conventions based on the orientation of the electrical field. The general solution may be expressed as:

fv(φ)=cos(vφ+ψ),

gv(φ)=−sin(vφ+ψ),

where ψ is the phase. In the convention, ψ=0 for even modes and ψ=π/2 for odd nodes.

In the center of the fiber, when r approaches 0, solutions are only functions of Jv, or Iv as Yv and Kv diverge. When the effective index is higher than the index in the cladding, the solution in the cladding is a function of Kv because Iv diverges when r goes to infinity. When the effective index is lower than the cladding, the solution oscillates, and there is a leaky mode, because the Jv and Yv Bessel functions oscillate.
A step-index optical fiber with only a core and a cladding has the characteristic equation:



  
    
      
        
          [
          
            
              
                J
                v
                ′
              
              
                uJ
                v
              
            
            +
            
              
                K
                v
                ′
              
              
                wK
                v
              
            
          
          ]
        
        ⁡
        
          [
          
            
              
                J
                v
                ′
              
              
                uJ
                v
              
            
            +
            
              
                
                  n
                  2
                  2
                
                
                  n
                  1
                  2
                
              
              ⁢
              
                
                  K
                  v
                  ′
                
                
                  wK
                  v
                
              
            
          
          ]
        
      
      =
      
        
          
            
              [
              
                
                  v
                  ⁢
                  
                    
                  
                  ⁢
                  β
                
                
                  kn
                  1
                
              
              ]
            
            2
          
          ⁡
          
            [
            
              
                1
                
                  u
                  2
                
              
              +
              
                1
                
                  w
                  2
                
              
            
            ]
          
        
        2
      
    
    ,
    
      

    
    ⁢
    where
  




  
    
      
        u
        2
      
      +
      
        w
        2
      
    
    =
    
      
        V
        2
      
      =
      
        
          
            [
            
              k
              ⁢
              
                
              
              ⁢
              ρ
              ⁢
              
                
                  
                    n
                    1
                    2
                  
                  -
                  
                    n
                    2
                    2
                  
                
              
            
            ]
          
          2
        
        .
      
    
  



Also, n1 and n2 are the respective indices of the core and cladding, and Jv′ and Kv′ are the first derivatives of Jv and Kv with respect to μ and ω. The values of μ and ω which solve this equation form a mode with a corresponding neff given by:





  
    n
    eff
  
  =
  
    
      
        
          n
          1
          2
        
        -
        
          
            (
            
              u
              
                k
                ⁢
                
                  
                
                ⁢
                ρ
              
            
            )
          
          2
        
      
    
    =
    
      
        
          
            n
            2
            2
          
          +
          
            
              (
              
                w
                
                  k
                  ⁢
                  
                    
                  
                  ⁢
                  ρ
                
              
              )
            
            2
          
        
      
      .
    
  




For a given v parameter, m corresponds to the mth root of the characteristic equation, starting from the highest neff. When v the right hand portion of the equation is zero, and the modes are either TE0,m or TM0,m for solutions of:





  
    
      [
      
        
          
            J
            v
            ′
          
          
            uJ
            v
          
        
        +
        
          
            K
            v
            ′
          
          
            wK
            v
          
        
      
      ]
    
    =
    0
  
  ,
  
    
      and
      ⁢
      
        
      
      [
      
        
          
            
              ɛ
              1
            
            
              ɛ
              2
            
          
          ⁢
          
            
              J
              v
              ′
            
            uJv
          
        
        +
        
          
            K
            v
            ′
          
          
            wK
            v
          
        
      
      ]
    
    =
    0
  
  ,




respectively. TE modes are odd and TM modes are even. When v>0, there are two sets of solution, positive solutions and negative solutions of the quadratic equation, where positive or negative refers to the sign in from of the square root. EHv,m modes are associated with the positive solutions, and HEv,m modes are associated with the negative solutions.

When n1 is very close to n2, a weakly guiding approximation is used, and it is assumed that:




  
    
      n
      2
      2
    
    
      n
      1
      2
    
  
  ≈
  1




The characteristic equation becomes:





  
    
      
        
          J
          l
        
        ⁡
        
          (
          u
          )
        
      
      
        
          uJ
          
            l
            -
            1
          
        
        ⁡
        
          (
          u
          )
        
      
    
    =
    
      -
      
        
          
            K
            l
          
          ⁡
          
            (
            w
            )
          
        
        
          
            wK
            
              l
              -
              1
            
          
          ⁡
          
            (
            w
            )
          
        
      
    
  
  ,




which may be solved for LP modes.

Solving modes for multilayered fibers is slightly more complicated, because the ez and hz field components are solved for each refractive index interface between layers. For the innermost and outermost layers, there are only two independent variables. For the innermost layer, the solution is either a function of Jv or Iv. For the outermost layer, the cladding, the solution is a function of Kv only, assuming that neff is higher than n of the cladding. For the other layers, the four parameters A, B, A′, and B′ are independent, and solutions are linear combinations of Jv and Yv, or Iv and Kv. Thus, the number of independent parameters to be found is 2 for the center, 2 for the cladding, and 4(N−2) for the intermediate layers where N is the number of layers of the fiber, for a total of 4(N−1) parameters. For the N−2 intermediate nodes there are four equations to find four unknowns. These equations stem from the continuity of ez, hz, eφ, and hφ field components at the interfaces between two adjacent layers.
One technique for solving these equations is the layer-by-layer method. Beginning at the center of the fiber, it is assumed that constants A and A′ of the first layer are 1 and α, respectively. The field continuity condition is used to solve for A, B, A′, and B′ for the next layer, as a function of 1 and α. This is performed for each layer. There is a 4×4 linear system to solve for each layer. At the last layer, it is determined whether a solution for α exists which ensures the continuity of the fields, corresponding to verifying that the determinant of a 2×2 matrix is equal to zero. If the condition is not met, another value of neff is used.
Graded index fibers may also be solved numerically. The values of ez and hz are not expressed using Bessel functions, but as second order coupled differential equations, which may be transformed to four first-order coupled differential equations, which may be solved using the fourth-order Runge-Kutta method. This may be performed for each layer using the layer-by-layer method.
OAM modes are not eigenmodes of optical fibers. However, OAM modes may be built from linear combinations of vector modes. Thus, optical fibers may support OAM modes.
Field components are given by expressions of the form:

u(r,φ)=u0(r)ξ(φ),

where u is any component of the electric or magnetic field, and ξ(φ) is a sine or cosine function based on the component of the field and whether the mode is even or odd. From Euler's identity:

u(r,φ)even+ju(r,φ)odd=u0(r)e±jvφ,

where the imaginary number j corresponds to a π/2 phase difference between the even and odd fields. Thus, any combination of even and odd degeneracy of a given mode with a π/2 phase difference leads to a mode with a total angular momentum of v. This total angular momentum is the sum of SAM σ and OAM l. The combination of HE modes gives a circular polarization with the same direction of rotation as the orbital angular momentum, while the combination of EH modes gives a circular polarization with a direction of rotation opposed to the direction of the orbital angular momentum. That is:





  
    OAM
    
      
        m
        -
        1
      
      ,
      
        ±
        l
      
    
    
      ±
      σ
    
  
  =
  
    
      
        (
        
          
            x
            ^
          
          ±
          
            jσ
            ⁢
            
              y
              ^
            
          
        
        )
      
      ⁢
      
        
          F
          →
        
        ⁡
        
          (
          r
          )
        
      
      ⁢
      
        e
        
          ±
          jlφ
        
      
    
    =
    
      {
      
        
          
            
              
                
                  
                    HE
                    
                      
                        l
                        +
                        1
                      
                      ,
                      m
                    
                    even
                  
                  ±
                  
                    j
                    ⁢
                    
                      
                    
                    ⁢
                    
                      HE
                      
                        
                          l
                          +
                          1
                        
                        ,
                        m
                      
                      odd
                    
                  
                
                ,
              
            
            
              
                σ
                =
                
                  ±
                  1
                
              
            
          
          
            
              
                
                  EH
                  
                    
                      l
                      -
                      1
                    
                    ,
                    m
                  
                  even
                
                ±
                
                  j
                  ⁢
                  
                    
                  
                  ⁢
                  
                    EH
                    
                      
                        l
                        -
                        1
                      
                      ,
                      m
                    
                    odd
                  
                
              
            
            
              
                σ
                =
                
                  ∓
                  1
                
              
            
          
        
        ,
      
    
  




where {circumflex over (x)} and ŷ are unit vectors and {right arrow over (F)}(r) is the radial dependence of either the electrical or magnetic field.

OAM modes may be laser modes in free-space or fiber modes. OAM modes cannot be built from HE1,m, because l=0 and:

HE1,meven±jHE1,modd=OAMm−1,0±.

The resulting mode only has a circular polarization. The combination of TM and TE modes gives a total angular momentum of zero, but leads to an unstable vortex, because these modes have different propagation constants, and:

TM0,m±jTE0,m=OAMm−1,±1∓,

HEv,m modes are the basis of OAMm−1,v−1 modes, with a circular polarization in the same direction as the field rotation, while EHv,m modes are the basis of OAMm−1,v+1 modes with a circular polarization in the opposite direction than the field rotation. This is summarized by:

OAMp,±l±=HEl+1,p+1even±jHEl+1,p+1odd 

OAMp,±l∓=EHl−1,p+1even±jEHl−1,p+1odd,

where the sign of the superscript denotes the direction of the circular polarization and the sign of l denotes the direction of field rotation. OAM modes made from HE modes have their polarization and direction of rotation in the same direction, while OAM modes made from EH modes have their polarization and direction of rotation in the opposite direction. HEl+1,m and EHl−1,m are eigenmodes for LPl,m modes. The first OAM mode is built from HE2,1, and has a circular polarization in the same direction as the orbital angular momentum with:

HE2,meven±jHE2,modd=OAMm−1,±1±.

Higher OAM modes exist with four variants, two circular polarizations and two OAM directions. Table 1 illustrates the four degeneracies of OAM modes. The number of supported modes includes all degeneracies in polarization and direction of rotation of the phase front of the fields. The two OAMp,0 modes are the HE1,p+1 modes with left and right circular polarization. The two OAMp,1 modes have the same direction of circular polarization and rotation. Higher order OAM modes yield four independent channels, because they have polarization and rotation in the same or in opposite directions.









TABLE 1






Counterclockwise
Clockwise



Polarization
Polarization







Counterclockwise Phase Rotation
OAM−p, −|l|
OAM+p, −|l|


Clockwise Phase Rotation
OAM−p, +|l|
OAM+p, +|l|









These equations are solved for the propagation constant β of the guided modes in a fiber. The propagation constant β and the effective index neff are directly related. From β, the dispersion and propagation of the made are calculated, including non-linear effects, using a variant of the non-linear Schrödinger equation. This applies to vector modes, and thus also to OAM modes. When there is interaction between the modes during propagation, the model may also consider mode coupling.
Mode coupling largely stems from irregularities in the fiber, and may be caused by diameter variations, loss, and isolated particles. Fiber micro and macro bends are also a source of mode coupling. Strong coupling occurs between modes with similar propagation constants, while weak coupling occurs between modes with propagation constant differences greater than a threshold. There are statistical models for crosstalk in multicore fibers. Modal dispersion may be considered as a higher order generalized form of polarization mode dispersion. In WDM systems, there may be coupling between different wavelengths caused by non-linear effects. There are two kinds of intra-modal coupling: crosstalk and multi-path interference (MPI). Crosstalk is more generic, and refers to the amount of power which leaks from a given mode to other modes. Crosstalk may occur at the multiplexing and demultiplexing levels, as well as during transmission. Multi-path interference is more specific, and refers to power which goes back and forth from one mode to other modes. Other impairments may occur during transmission, such as losses, for example bending losses. A model may be used for bending loss in few mode fibers. Also, simulation on how twisted fiber interacts with the generation and propagation of OAM modes may be used.
A model is developed to describe OAM modes propagating in an optical fiber. The model includes coupling between the laser source and the optical fiber, propagation of OAM modes, and intra-modal coupling. Laser modes and fiber modes have different characteristics. When a free-space OAM beam enters a fiber, there is a transient regime where part of the energy is converted to the OAM mode and transmitted through the fiber, and part of the energy is radiated out of the fiber. It is also possible that part of the energy is converted into other fiber modes. The model estimates crosstalk and losses induced by the coupling between the laser beam and fiber, as a function of the incident beam geometry and fiber geometry.
For a given fiber profile at a given wavelength, the supported vector modes are calculated. The refractive index of a layer may be a fixed index or a material, where the refractive index is a function of the wavelength. Example materials include silica, air, silica doped with germanium (SiO2GeO2), or silica doped with fluorine (Si2F). The dopant concentration may be varied. For silica, the refractive index is calculated using the Sellmeier equation, while for doped silica, the refractive index is calculated using the Classius-Mossotti equation. Software also calculates the refractive index difference with a reference layer, and the numerical aperture relative to the next layer. The inner radius, outer radius, and/or width of the layer may be entered by a user. There may be graded profiles modeled using small steps. The fiber profile is plotted.
Fiber modes may be solved for a range of wavelengths, as well as a range of v values. When the mode is solved, the list of solved wavelengths is displayed, along with the name of the mode, the effective index of the mode, the distance to the nearest mode, and the dispersion parameter.
The whole algorithm is considered as the characteristic equation of the fiber. To find the zeros of this function, the range is divided into small refractive index intervals, and the function is evaluated for those indices, and changes in sign are examined. When an interval where the sign changes is found, the Brent algorithm is used to find the position of the zero in this interval. Additional testing is performed to determine that the found zero is a zero, not a discontinuity. Multiple calculations may be performed in parallel.
For each solved mode, the electrical and magnetic fields may be plotted. The plots represent a component of the field, either in cylindrical coordinates or Cartesian coordinates. The Poynting vector may also be plotted. Additionally, the phase and intensity of the OAM mode may be plotted.
To perform optimization, a wavelength is selected for optimization. Then, parameters, such as radius and refractive index, are swept over. After the sweep over is performed, the results are displayed with the mode separation between the nearest modes, the number of supported vector modes, the number of modes with a mode separation greater than 10−4, and the average mode separation. Then, it may be solved for the new fiber profile over a range of wavelengths.
From the calculated vector modes, the OAM modes may be calculated. Also, intrinsic parameters, such as effective index of the modes, propagation constants, mode confinements, dispersion parameter, non-linear parameter, etc. may be calculated. The intrinsic parameters are put in relation with extrinsic parameters, such as fiber bend radius and statistical parameters, to simulate mode coupling and losses.
To transmit OAM modes, it is desirable that eigenmodes not easily couple to each other to form LP modes. A mode index separation of 10−4 may be sufficient to prevent this coupling. It may also be desirable for the fiber profile to mimic the distribution of the mode intensity.
Silica cladding has a refractive index of 1.457 at 633 nm. The fiber numerical aperture (NA) may be less than 0.4, because, when the GeO2 concentration in the doped region is too high, the preform may break. A NA of, for example, around 0.3 may be used. The refractive index of the doped region may be about 0.03 to 0.04 higher than the refractive index of the cladding. Also, when the region is doped with F to obtain a refractive index lower than the cladding, the minimum difference in refractive index between the doped region and silica may be −0.005.
In designing a fiber, first a simple fiber profile is designed. Then complexity is added, step by step. Optimization is performed by slowly varying the diameter of the layers by small increments, and observing the number of supported modes, and average and maximum mode index separation.
It is desirable to have a fiber with multiple OAM modes with a large mode effective index separation. For example, the fiber may support at least 16 modes in the C-band (1530 nm-1565 nm) with a mode index separation of at least 10−4. In the design, there may be a margin to account for margins in layer diameters and refractive indices and the smoothness of layer edges. Also, when the effective index of the mode is close to the index of the cladding, the mode is sensitive to fiber bends and imperfections.
The number of supported modes may be adjusted by varying the diameter and the refractive index of the region of the fiber where the modes are located. A fiber with a higher refractive index contrast or a larger fiber layer supports more modes.
The mode separation is affected by the refractive index contrast between the fiber layers. The LP modes are present in the weakly guiding approximation, where the refractive indices of the core and cladding are usually different. Thus, a higher refractive index contrast prevents LP modes, leading to a greater mode index separation.
The width of the doped region in the fiber is adjusted to adjust the number of supported vector modes to have fewer modes with a wider spacing. An external layer material having a refractive index lower than the cladding is added to increase the contrast of the refractive indices and limit the number of higher order modes.
The number of supported modes also affects the mode index separation, because the higher the number of supported modes, the greater the chances of having modes with similar effective indices, even between vector modes belonging to different LP modes.
A vortex fiber is characterized by a central core and an outer ring. The modes with l=0 are Gaussian, which is good for coupling efficiency to a laser beam. FIG. 13 illustrates index profile 340 for a vortex like fiber. Table 2 illustrates the characteristics of the layers of the vortex like fiber, including the outer diameter φ, outer radius ρ, layer width Δρ, refractive index n, and refractive index relative to the cladding Δn.











TABLE 2





φ (μm)
ρ (μm)
Δρ (μm)
n
Δn



















3.60
1.80
1.80
1.466
+0.022


7.78
3.89
2.09
1.440
−0.004


10.38
5.19
1.30
1.476
+0.032


14.52
7.26
2.07
1.435
−0.009


125.00
62.5
55.24
1.444
0.000









FIG. 14 illustrates a graph of mode indices over the C-band. Curve 342 shows the effective refractive index for mode HE1,1, curve 344 shows the effective refractive index for mode TE0,1, curve 346 shows the effective refractive index for mode HE2,1, curve 348 shows the effective refractive index for mode TM0,1, curve 350 shows the effective refractive index for mode EH1,1, curve 352 shows the effective refractive index for mode HE3,1, and curve 354 shows the effective refractive index for mode HE1,2. Table 3 illustrates the effective refractive index and the difference between the effective index of the nearest mode for the modes in the vortex-like fiber. This fiber potentially supports 10 OAM modes with mode separation of at least 1.96·10−4.










TABLE 3






Mode
neff
Δneff








HE1, 1
1.4544
1.93 · 10−3



TE0, 1
1.4525
2.84 · 10−4



HE2, 1
1.4522
2.84 · 10−4



TM0, 1
1.4518
3.34 · 10−4



EH0, 1
1.4488
1.99 · 10−4



HE3, 1
1.4486
1.96 · 10−4



HE1, 2
1.4484
1.96 · 10−4









A hollow core fiber has a high index contrast, because the refractive index of air is 1. An embodiment air core ring fiber supports relatively few well separated modes. FIG. 15 illustrates profile 360 for a hollow core fiber. This fiber may support 40 OAM modes, up to OAM0,∓10∓ with a mode index separation of at least 1.07·104. Table 4 shows the outer later diameter, outer radius, layer width, refractive index, and refractive index relative to the cladding.











TABLE 4





φ (μm)
ρ (μm)
Δρ (μm)
n
Δn



















18.40
9.20
9.20
1.000
−0.457


23.20
11.60
2.40
1.498
+0.041


31.20
15.60
4.00
1.452
−0.005


125.0
62.50
46.90
1.457
0.000









FIG. 16 illustrates a graph of effective index verses wavelength, where line 404 is 1.55 μm and curve 360 is the refractive index of the cladding. Curve 362 shows the effective refractive index for mode EH10,1, curve 364 shows the effective refractive index for mode HE12,1, curve 366 shows the effective refractive index for mode EH9,1, curve 368 shows the effective index for mode HE11,1, curve 370 shows the effective index for mode EH8,1, curve 372 shows the effective index for mode HE10,1, curve 374 shows the effective index for mode EH7,1, curve 376 shows the effective index for mode HE9,1, curve 378 shows the effective index for mode EH6,1, curve 380 shows the effective index for mode HE8,1, curve 382 shows the effective index for mode EH5,1, curve 384 shows the effective index for mode HE7,1, curve 386 shows the effective index for mode EH4,1, curve 388 shows the effective index for mode HE6,1, curve 390 shows the effective index for mode EH3,1, curve 392 shows the effective index for mode HE5,1, curve 394 shows the effective index for mode EH2,1, curve 396 shows the effective index for mode HE4,1, curve 398 shows the effective index for mode EH1,1, curve 400 shows the effective index for mode TM0,1, curve 402 shows the effective index for mode HE3,1, curve 405 shows the effective refractive index for mode HE2,1, curve 406 shows the effective index for mode HE1,1, and curve 408 shows the effective index for mode TE0,1. The fundamental mode is TE0,1.
An embodiment air core ring fiber supports 16 OAM modes. Other embodiments may support 40 OAM modes with a minimum mode index separation of 1.07·10−4. The mode effective indices of TE0,1 and TM0,1 are so well separated from HE2,1 that the mode effective index of TE0,1 is closer to HE1,1, and the mode effective index of TM0,1 falls between HE3,1 and EH1,1 (LP2,1). The supported modes have m=1, because of the thinness of the ring. This simplifies the multiplexing and demultiplexing of the OAM modes, and prevents crosstalk between modes belonging to different families, facilitating good mode separation across the C-band. Also, the ring shape of the fiber profile mimics the intensity profile of the OAM beams and helps the efficiency in the coupling between the free space OAM beams and the fiber. Table 5 illustrates the modes, their effective refractive index, and the minimum refractive index spacing.










TABLE 5






Mode
neff
Δneff








TE0, 1
1.4662
1.07 · 10−4



HE1, 1
1.4661
1.07 · 10−4



HE2, 1
1.4657
3.50 · 10−4



HE3, 1
1.4651
5.87 · 10−4



TM0, 1
1.4645
2.56 · 10−4



EH1, 1
1.4643
1.26 · 10−4



HE4, 1
1.4641
1.26 · 10−4



EH2, 1
1.4635
6.13 · 10−4



HE5, 1
1.4628
4.71 · 10−4



EH3, 1
1.4623
4.71 · 10−4



HE6, 1
1.4612
3.76 · 10−4



EH4, 1
1.4608
3.76 · 10−4



HE7, 1
1.4591
3.07 · 10−4



EH5, 1
1.4588
3.07 · 10−4



HE8, 1
1.4567
2.51 · 10−4



EH6, 1
1.4565
2.51 · 10−4



HE9, 1
1.4540
2.02 · 10−4



EH7, 1
1.4538
2.02 · 10−4



HE10, 1
1.4509
1.57 · 10−4



EH8, 1
1.4507
1.57 · 10−4



HE11, 1
1.4474
1.12 · 10−4



EH9, 1
1.4473
1.12 · 10−4









An embodiment optical fiber including a large low refractive index (e.g. hollow) core, surrounded by a thin layer of material with a high refractive index, transmits multiple OAM modes. The high refractive index material is surrounded by another slightly wider layer of material with a low refractive index, all surrounded by silica cladding. The high contrast is introduced between the low refractive index of the core and the high refractive index of the highly doped region to achieve a large discrepancy between the effective indices of the HEt+1,m and EHl−1,m modes. When the inner core is not hollow, it may be silica doped with fluorine to reduce the refractive index. Also, the low refractive index layer may be silica doped with fluorine to reduce the refractive index. In another example, boron trioxide is used to reduce the refractive index. The high refractive index region may be silica doped with germanium to increase the refractive index. In another example, aluminum oxide is added to raise the refractive index. Phosphorous may also be added to the fiber. The spacing in the effective refractive index in the modes is achieved by setting the refractive index of the highly doped area to be high, to increase the contrast between the high refractive index layer and both the core and the low refractive index material.
Modes propagate through the high refractive index layer. The ring shape of this region matches the intensity profile of OAM modes. The number of supported modes is determined by the inner and outer diameter of the high refractive index material, as well as the refractive index of the high refractive index material. Because the effective indices of the modes reside between the refractive index of the cladding and the refractive index of the highly doped layer, increasing the number of supported modes leads to more closely spaced effective refractive indices. Fewer modes are achieved by making the highly doped layer narrower. Also, the narrowness of the high refractive index material removes modes with m>1. Fiber dimensions determine the achievable refractive indices and the number of supported nodes.
A custom air-core ring fiber is fabricated. A fiber preform is fabricated using modified chemical vapor deposition (MCVD) in two steps. First, a layer is formed using SiO2, P2O5, and F. Then, a ring layer is deposited using SiO2 and GeO2 to produce a step index profile. Finally, the fiber is partially collapsed to produce the glass preform. The tolerance for the radius of the air core is 9.2±0.2 μm, while the tolerance for the high refractive index layer is 2.4±0.2 μm. Also, the tolerance for the low refractive index layer is 4.0 μm or greater. The refractive index for the high refractive index layer is 0.036 to 0.041 higher than the refractive index of the cladding. Other embodiments have different tolerances. For example, when the refractive index of the ring layer is lower, the thickness is adjusted.
FIG. 17 illustrates the refractive index profile of a designed fiber in curve 510 and the refractive index profile of a fabricated fiber in curve 512. The non-zero air core refractive index measurement is an artifact from the measurement technique. The refractive index of oil which seeps into the air core, rather than the refractive index of air. The fabricated fiber has a lower refractive index than the designed fiber in the high refractive index region. This is because doping silica to increase its refractive index weakens the glass, making it problematic to achieve a very high refractive index.
FIG. 18 illustrates fiber 500, a fabricated air core ring fiber. Fiber 500 contains air core 506, high refractive index region 504, low refractive index region 502, and cladding 508. Air core 506 has a radius of 8.16 μm and a diameter of 16.32 μm. Low refractive index layer 502 has an inner radius of 10.10 μm and an outer radius of 15.39 μm.
Different types of fibers and lenses may be used for free space light coupling into air core ring fibers. Fibers which may be used include single mode fiber, multimode fiber (MMF), and orbital angular momentum fibers. FIG. 19 illustrates system 720 using a SMF to couple light into a test fiber. Laser 722, a continuously tunable laser couples light into SMF 724. Lens 726 converts the light to a free space beam. Power attenuator 728, which contains half wave plate (HP) 730 and polarization beam splitter (PBS) 732, maintains the linear polarization of the light. Mirrors 734 and 736 control the height of the beam from the optical table. The light beam is coupled into fiber 742, the fiber under test, by lens or microscopic objective (MO) 738 through stage 740. The power at the fiber output is monitored by detector 744, a power meter. To form circularly polarized light, a quarter wave plate is placed between the beam splitter and the mirror.
The input power from the laser is 6 dBm. Table 6 below illustrates the characterizations of the three fibers, while table 7 below shows power coupling measurements using a lens, and table 8 shows power coupling measurements using a microscope objective.











TABLE 1








Inner Core
Outer Core


Fiber
Numerical Aperture
Length
Diameter
Diameter





















SMF
0.14
1
m
N/A
8.2
μm


MMF
0.2
2.4
m
N/A
50
μm


OAM
0.32
1
m
15 μm
35
μm





















TABLE 2








Lens
Power at lens
Power at fiber


Fiber
Length
Polarization
NA
output
output






















SMF
1 m
Linear
0.3
400
μW
130
μW






−4
dBm
−8.86
dBm


OAM
1 m
Linear
0.3
400
μW
30
μW






−4
dBm
−15.23
dBm


OAM
1 m
Circular
0.3
400
μW
60
μW






−4
dBm
−12.22
dBm





















TABLE 3








Lens
Power at lens
Power at fiber


Fiber
Length
Polarization
NA
output
output























SMF
1
m
Linear
0.17 &
~320
μW
32
μW













0.65
−4.95
dBm
−14.95
dBm















MMF
2.4
m
Linear
0.17 &
~320
μW
240
μW













0.65
−4.95
dBm
−6.2
dBm















OAM
1
m
Circular
0.17 &
~350
μW
115
μW













0.65
−4.56
dBm
−9.4
dBm















OAM
1
m
Linear
0.17 &
~320
μW
105
μW













0.65
−4.95
dBm
−9.8
dBm









The numerical aperture (NA) in a simple fiber is given by:

NA=√{square root over (n12−n22)}=(n)sin θmax,

where n1 is the refractive index of the core, n2 is the refractive index of the cladding, n is the refractive index of the incident medium, and θmax is the maximum acceptance half angle. FIG. 20 illustrates a simple fiber with core 756, cladding 754 and 758, and maximum angle 752.

An OAM fiber, illustrated by FIG. 21, has many layers. Fiber 760 has low refractive index core 761, high refractive index regions 762 and 763, low refractive index regions 764 and 765, cladding 766 and 767, and maximum angles 768 and 769. For a supported mode, the intensity of the field is located in high refractive index regions 762 and 763. Thus, rays of light are reflected either at the interface of core 761 (θ12) or low refractive index regions 764 and 765 (023), because the refractive index contrast is higher. Thus, the maximum acceptance angle for the fiber is limited by the lowest acceptance angle. The NA is given by:

NA=√{square root over (n32−n22)}.

The characteristics of the OAM fiber are slightly different at its ends. FIG. 22 illustrates a plot of the measured refractive indices on the x and y axis at both ends of the fiber. Curve 778 shows the refractive index at the start of the fiber on the x axis, curve 776 shows the refractive index of the fiber at the start on the y axis, curve 774 shows the refractive index of the fiber at the end on the x axis, and curve 772 shows the refractive index of the fiber at the end on the y axis. The NA is calculated using the maximum value of the indices as n2 and the minimum value as n3. Table 9 shows the NAs. These are the upward bound for the NA. To calculate a lower bound for NA, values from FIG. 22 which are slightly lower than the lowest refractive index in the n2 region and slightly higher than the highest refractive index in the n3 region. Using n3=1.483 and n2=1.451 yields a NA of 0.306.











TABLE 9






Position
n3
n2
NA








Start x axis
1.48692
1.45074
0.326



Start y axis
1.48712
1.45006
0.330



End x axis
1.48421
1.44995
0.317



End y axis
1.48406
1.45040
0.314









An embodiment fiber supports 16 OAM modes, up to OAM0,±4±, with a minimum mode index separation of 1.12·10−4. FIG. 23 illustrates the measured refractive index on the x-axis (curve 520) and the measured refractive index on the y-axis (curve 522). Curve 524 shows the refractive index from a simulation. Table 10 shows the outer diameter, outer radius, layer width, refractive index, and refractive index difference relative to the cladding for the fiber.











TABLE 10





φ (μm)
ρ (μm)
Δρ (μm)
n
Δn



















16.2
8.10
8.10
1.000
−0.444


20.00
10.00
1.90
1.472
+0.028


30.00
15.00
5.00
1.139
−0.005


125.0
62.50
47.50
1.444
0.000









An embodiment fiber experiences good separation between the effective refractive indices. At 1550 nm, the closest modes are TE0,1 and HE1,1, which have a separation of 1.1·10−4. FIG. 24 illustrates a graph of the modes in the fiber. Line 530 shows 1.55 μm. Curve 532 shows n4, the refractive index of the cladding, curve 534 shows mode EH3,1, curve 536 shows mode HE5,1, curve 538 shows mode EH2,1, curve 540 shows mode HE4,1, curve 542 shows mode EH1,1, curve 544 shows mode TM0,1, curve 546 shows mode HE3,1, curve 548 shows mode HE2,1, curve 550 shows mode HE1,1, and curve 552 shows mode TE0,1. Table 11 shows the refractive index and refractive index separation for modes in the fiber. Because all the vector modes except the TE and TM modes may be the basis of two OAM modes, one with clockwise rotation and one with counterclockwise rotation, the fiber may support 16 OAM modes. Furthermore, the effective index separation is high, leading to little inter-modal coupling.










TABLE 11






Mode
neff
Δneff








TE0, 1
1.4490
1.12 · 10−4



HE1, 1
1.4489
1.12 · 10−4



HE2, 1
1.4485
3.99 · 10−4



HE3, 1
1.4477
3.83 · 10−4



TM0, 1
1.4474
3.38 · 10−4



EH1, 1
1.4470
3.38 · 10−4



HE4, 1
1.4466
4.62 · 10−4



EH2, 1
1.4461
4.89 · 10−4



HE5, 1
1.4449
3.52 · 10−4



EH3, 1
1.4446
3.52 · 10−4









FIG. 25 illustrates a graph of group index verses wavelength for modes in an embodiment air core ring fiber. Curve 560 shows the group index for mode TM0,1 and curve 562 shows the group index for TE0,1 mode. Curve 564 shows the group index for HE1,1, which are two OAM0,0 modes. Also, curve 568 shows the group index for mode HE2,1, which are two OAM0,1 modes. Curve 566 shows mode EH1,1 and curve 570 shows the group index for mode HE3,1, which are four OAM0,2 modes. Additionally, curve 572 shows the group index for mode EH2,1 and curve 574 shows the group index for HE4,1, which are four OAM0,3 modes. Curve 576 shows the group index for mode EH3,1 and curve 578 shows the group index for mode HE5,1, which are OAM0,4, four modes. The OAM modes with higher topological charge may propagate slower than OAM modes with lower topological charge. Also, among OAM modes with the same topological charge, the modes with polarization in the same direction as the rotation propagate slower than the modes with polarization in the opposite direction as the rotation.
FIG. 26 shows the simulated dispersion versus wavelength for the modes in an embodiment air core ring fiber. Curve 580 shows the dispersion for mode EH3,1, and curve 588 shows the dispersion for mode HE5,1, which are four OAM0,4 modes. Also, curve 582 shows the dispersion for mode EH2,1, and curve 590 shows the dispersion for mode HE4,1, four OAM0,3 modes. Additionally, curve 584 shows the dispersion for mode TM0,1 and curve 598 shows the dispersion for mode TE0,1. Curve 586 shows the dispersion for mode EH1,1, and curve 592 shows the dispersion for mode HE3,1, four OAM0,2 modes. Also, curve 594 shows the dispersion for mode HE2,1, two OAM0,1 modes, and curve 596 shows the dispersion for mode HE1,1, two OAM0,0 modes. The dispersion for all of the modes across the C band is negative, because of the profile of the fiber. It is problematic to design a low dispersion fiber, especially across a range of wavelengths. A dispersion compensation system may be used. However, dispersion is not problematic over short distance. Also, the dispersion reduces the risk of nonlinear effects.
The refractive index of an embodiment air core ring fiber is measured using an optical fiber analyzer.
Using an SLM, a laser beam is converted into an OAM beam, coupled into an embodiment air core ring fiber. Then, the signal at the end of the fiber is converted back to free-space, to a Gaussian beam using the SLM, and the output power is recorded. This is done one mode at a time.
The effective index difference between OAM modes is determined at different wavelengths in the C-Band. This may be determined using interferometry.
Then, crosstalk is measured, by generating a single OAM mode. The power of surrounding modes is measured at the receiver. The impact of bends and twists on crosstalk may also be examined.
FIG. 27 illustrates flowchart 820 for a method of transmitting OAM modes in an air core ring optical fiber. In one example, the air core fiber includes a first layer with a low refractive index, for example air, a high refractive index ring layer surrounding the low refractive index layer, a relatively low refractive index layer surrounding the high refractive index layer, and a cladding surrounding the relatively low refractive index layer. Initially, in step 822, the OAM light mode is generated. For example, a Gaussian beam may be generated by a laser, and then converted to an OAM mode in free space. OAM beams may be generated by transforming the phase of a laser beam or a Gaussian beam. The phase transformation may be done through CGH or SPP. Alternatively, PICs are used. In another example, a multicore fiber is used to generate and recover OAM modes. Alternatively, each channel is a separate ring of light into the fiber.
Then, in step 824, the OAM is coupled in the fiber. The OAM mode may be coupled into the fiber using a lens or microscope objective.
Next, in step 826, the OAM mode propagates along the fiber. Many OAM modes may propagate along the fiber with a significant separation between modes.
Finally, in step 828, the transmitted OAM mode is detected. A CCD detector may be used. Alternatively, photodiodes are used. In one example, demultiplexing is performed at the same time as detection. In another example, the OAM modes are demultiplexed and then detected in separate steps. For example, the CGH technique creates spatial separation of the converted modes, and it may be used to convert and demultiplex simultaneously. The OAM beams may be directly detected. Also, OAM beams may be converted to Gaussian beams, for example so a receiver designed for a Gaussian intensity shaper.
An embodiment method of transmitting light includes generating a plurality of OAM modes and coupling the plurality of OAM modes into an optical fiber. The optical fiber includes a first layer having a first refractive index, a second layer surrounding the first layer, where the second layer has a second refractive index, a third layer surrounding the second layer, where the third layer has a third refractive index, and a fourth layer surrounding the third layer, where the fourth layer has a fourth refractive index, where the first refractive index is less than the second refractive index, where the third refractive index is less than the second refractive index, and where the third refractive index is less than the fourth refractive index. The method also includes propagating the plurality of OAM modes along the optical fiber to produce a plurality of propagated OAM modes.
FIG. 28 illustrates a block diagram of processing system 270 that may be used for implementing the devices and methods disclosed herein. For example, a processing system may be used to determine the width and refractive index of rings in an optical fiber. Specific devices may utilize all of the components shown, or only a subset of the components, and levels of integration may vary from device to device. Furthermore, a device may contain multiple instances of a component, such as multiple processing units, processors, memories, transmitters, receivers, etc. The processing system may comprise a processing unit equipped with one or more input devices, such as a microphone, mouse, touchscreen, keypad, keyboard, and the like. Also, processing system 270 may be equipped with one or more output devices, such as a speaker, a printer, a display, and the like. The processing unit may include central processing unit (CPU) 274, memory 276, mass storage device 278, video adapter 280, and I/O interface 288 connected to a bus.
The bus may be one or more of any type of several bus architectures including a memory bus or memory controller, a peripheral bus, video bus, or the like. CPU 274 may comprise any type of electronic data processor. Memory 276 may comprise any type of non-transitory system memory such as static random access memory (SRAM), dynamic random access memory (DRAM), synchronous DRAM (SDRAM), read-only memory (ROM), a combination thereof, or the like. In an embodiment, the memory may include ROM for use at boot-up, and DRAM for program and data storage for use while executing programs.
Mass storage device 278 may comprise any type of non-transitory storage device configured to store data, programs, and other information and to make the data, programs, and other information accessible via the bus. Mass storage device 278 may comprise, for example, one or more of a solid state drive, hard disk drive, a magnetic disk drive, an optical disk drive, or the like.
Video adaptor 280 and I/O interface 288 provide interfaces to couple external input and output devices to the processing unit. As illustrated, examples of input and output devices include the display coupled to the video adapter and the mouse/keyboard/printer coupled to the I/O interface. Other devices may be coupled to the processing unit, and additional or fewer interface cards may be utilized. For example, a serial interface card (not pictured) may be used to provide a serial interface for a printer.
The processing unit also includes one or more network interface 284, which may comprise wired links, such as an Ethernet cable or the like, and/or wireless links to access nodes or different networks. Network interface 284 allows the processing unit to communicate with remote units via the networks. For example, the network interface may provide wireless communication via one or more transmitters/transmit antennas and one or more receivers/receive antennas. In an embodiment, the processing unit is coupled to a local-area network or a wide-area network for data processing and communications with remote devices, such as other processing units, the Internet, remote storage facilities, or the like.
While several embodiments have been provided in the present disclosure, it should be understood that the disclosed systems and methods might be embodied in many other specific forms without departing from the spirit or scope of the present disclosure. The present examples are to be considered as illustrative and not restrictive, and the intention is not to be limited to the details given herein. For example, the various elements or components may be combined or integrated in another system or certain features may be omitted, or not implemented.
In addition, techniques, systems, subsystems, and methods described and illustrated in the various embodiments as discrete or separate may be combined or integrated with other systems, modules, techniques, or methods without departing from the scope of the present disclosure. Other items shown or discussed as coupled or directly coupled or communicating with each other may be indirectly coupled or communicating through some interface, device, or intermediate component whether electrically, mechanically, or otherwise. Other examples of changes, substitutions, and alterations are ascertainable by one skilled in the art and could be made without departing from the spirit and scope disclosed herein.

## Claims

1. An optical fiber comprising:
a first layer having a first refractive index, wherein the first layer is air;
a second layer surrounding the first layer, wherein the second layer has a second refractive index, an inner radius, and an outer radius; and
a third layer surrounding the second layer, wherein the third layer has a third refractive index, wherein the first refractive index is less than the second refractive index, wherein the third refractive index is less than the second refractive index, and wherein a ratio of the outer radius to the inner radius is less than 1.5, wherein the optical fiber supports at least 16 orbital angular momentum (OAM) modes with a mode spacing between all of the at least 16 OAM modes of at least 1·10−4.

2. The optical fiber of claim 1, wherein the first layer has a diameter of between 15 μm and 20 μm.

3. The optical fiber of claim 1, wherein the second layer is germanium doped silica.

4. The optical fiber of claim 1, wherein the ratio of the outer radius to the inner radius is less than 1.3.

5. The optical fiber of claim 1, wherein a width of the second layer is between 2 μm and 3 μm.

6. An optical fiber comprising:
a first layer having a first refractive index, wherein the first layer is air;
a second layer surrounding the first layer, wherein the second layer has a second refractive index;
a third layer surrounding the second layer, wherein the third layer has a third refractive index; and
a fourth layer surrounding the third layer, wherein the fourth layer has a fourth refractive index, wherein the first refractive index is less than the second refractive index, wherein the third refractive index is less than the second refractive index, and wherein the third refractive index is less than the fourth refractive index.

7. The optical fiber of claim 6, wherein the third layer is fluorine doped silica.

8. The optical fiber of claim 6, wherein the fourth layer is silica.

9. The optical fiber of claim 6, wherein the optical fiber is configured to support at least 16 orbital angular momentum (OAM) modes with a mode spacing of at least 1·10−4.

10. A method of transmitting light, the method comprising:
generating a plurality of orbital angular momentum (OAM) modes;
coupling the plurality of OAM modes into an optical fiber, wherein the optical fiber comprises
a first layer having a first refractive index,
a second layer surrounding the first layer, wherein the second layer has a second refractive index, an inner radius, and an outer radius, and
a third layer surrounding the second layer, wherein the third layer has a third refractive index, wherein the first refractive index is less than the second refractive index, wherein the third refractive index is less than the second refractive index, and wherein a ratio of the outer radius to the inner radius is less than 1.5; and

propagating the plurality of OAM modes along the optical fiber to produce a plurality of propagated OAM modes, wherein propagating the plurality of OAM modes comprises propagating at least 16 OAM modes with a mode spacing between all of the at least 16 OAM modes of at least 1·10−4.

11. The method of claim 10, wherein generating the plurality of OAM modes comprises utilizing a computer generated hologram (CGH).

12. The method of claim 10, wherein generating the plurality of OAM modes comprises utilizing a spiral phase pattern (SPP).

13. The method of claim 10, further comprising demultiplexing the plurality of propagated OAM modes.

14. The method of claim 10, further comprising detecting the plurality of propagated OAM modes.

15. An optical fiber comprising:
a first layer having a first refractive index;
a second layer surrounding the first layer, wherein the second layer has a second refractive index;
a third layer surrounding the second layer, wherein the third layer has a third refractive index; and
a fourth layer surrounding the third layer, wherein the fourth layer has a fourth refractive index, wherein the first refractive index is less than the second refractive index, wherein the third refractive index is less than the second refractive index, and wherein the third refractive index is less than the fourth refractive index, wherein the optical fiber is configured to support at least 16 orbital angular momentum (OAM) modes with a mode spacing of at least 1·10−4.

16. The optical fiber of claim 15, wherein the third layer is fluorine doped silica.

17. The optical fiber of claim 15, wherein the fourth layer is silica.