Siegel der Universität

Universität zu Köln
line
Mathematisch-Naturwissenschaftliche Fakultät
Fachgruppe Physik

I. Physikalisches Institut

Fourier gratings

Printer-friendly versionPDF version
developments

Developments

Workgroup for SubMillimeter and Terahertz Receivers


Fourier gratings

Contents

 

Introduction

Owing to their high diffraction efficiency, phase gratings are widely used in optics to shape and split monochromatic beams. With the advent of multibeam heterodyne receivers in submillimeter astronomy, they have also become popular in this part of the spectrum for local-oscillator beam multiplexing. Splitting the beam from one oscillator into a given number of equal beams by use of a phase grating is a very elegant way to distribute the local-oscillator power efficiently to an array of mixers.

Contents 

Stepped grating Structures

 

  • Standard in common optics ==> easy to develope.
  • Also used in heterodyne array receivers (CHAMP of the MPIfR).
  • Good, but not optimum performance
  • As reflection grating (they perform better than transmission gratings) no 2-dimensional beam patterns possible, because the 2-D structures are very difficult to manufacture.

These grating structures have been the starting point in our development of a new type of grating for the use as LO-multiplexer.

Structure of a Dammann grating ==> Difraction pattern of the Dammann grating
Structure of a 4 level stepped grating ==> Difraction pattern of the 4 level grating

Fig. 1: Stepped grating structure with the resulting beampattern (measurement of the far field intensity distribution). 2 Layer grating (top row) and an optimized 4 layer structure (second row).

Contents

Theory of Fourier gratings (FG)

A phase grating has a transmission function like

Transmission function

using the phase modulation at the position (x,y). The structure of a one dimensional symmetric Fourier grating is given by a finite Fourier series expansion :

Structure of a Fourier grating

where the unit cell extends from -Delta x/2 to +Delta x/2. The transmission function is accordingly written as

Transmission function of a FG

The diffraction pattern of the general Fourier grating is given by the Fourier transform of t(x). This results in a multiple convolution of the diffraction fields of the individual Fourier components:

Diffraction Pattern of a FG

where Jq(a) denotes the Bessel-function of the first kind of order q. Although this is a simple equation we can not analytically calculate the coefficients an. We thus developed a numerical optimization method to calculate them. A direct numerical FFT is usually faster than evaluating T(k), and is therefore used in our calculation process.

 

Structure of a 2 beam FG Structure of a 5 beam FG Structure of a 8 beam FG

Fig. 2: Structures of 13 coefficients Fourier gratings producing 2, 5 and 8 beams (from left to right). The theoretical efficiencies are 81%, 92% and 96%.

 

A comparison of the intensity distribution in the diffraction pattern using a standard Dammann grating (from left to right in Fig. 3), a multilevel grating, a 13 coefficient Fourier grating and and the optimum Fourier grating (right hand panel) shows the increase of the efficiency by using our type of grating.

Comparison Dammann grating with FG
Fig. 3

The efficiency of Fourier gratings producing a given number of equal beams (Fig. 4). The dotted line marks the maximum efficiency achievable with N=5 Fourier coefficients, the solid line gives the efficiency for N-->infinity.

Efficiencies for various FG's
Fig. 4

Number of Fourier coefficients is calculated from the tradeoff between:

  • minimum radius of curvature of the structure (limiting the size of the endmill)
  • efficiency of the grating (grows with the number of coefficients)

 

==>results in approximatly 13 Fourier coefficients (Fig. 5).

Number of Fourier coefficients
Fig. 5

 

Bandwidth of the gratings:

Detuning the frequency results either in position errors of the beams and in a power inbalance between them. Our SIS-system can tolerate a frequency shift of approximatley +-5%. For higher bandwidth we have to exchange the grating. Because the reference planes are machined together with the grating structure, exchanging the gratings is simple and reproduceble.

 

Contents

1-dimensional Fourier gratings

For the calculation of one-dimensional structures we use the following 2-step method:

  • Random search of a good set of coefficients an (e.g. 4 N sets for N Fourier-components). The calculated diffraction patterns of these parameter sets are evaluated by comparing the RMS deviation from the desired pattern.
  • Linear optimization of the best (e.g. 20) parameter sets to get the best possible grating structure for N Fourier-components.

The grating structures can be superimposed in different orthogonal directions to get a 2-dimensional beampattern. Of course, this method only produces rectangular arrengements.

 

Measured intensity distribution Structure of the 8 beam FG (superimposed 2 and 4 beam FG)

Fig. 6: Measured intensity distribution of a 490GHz beam, split by an 8 beam grating (left hand side). Contour levels are from 5% to 95% of the maximum peak intensity in steps of 10%. The grating consists of a 2 beam grating, superimposed with a 4 beam grating. The resulting surface topology of the unit cell using 13 Fourier coefficients for each direction is shown on the right hand side. Contour levels are in steps of .5 radians. The 1-dimensonal structures are plotted along the axis of the graph.

Contents

2-dimensional Fourier gratings

 

Calculation of two-dimensional structures:

Because of the higher number of coefficients N2 a random search is inefficient. We thus established the following method:

  • Modeling the grating as a beam combiner by summing up the fields of the desired diffraction orders as plane waves on the grating plane produces a grating with varying amplitude and phase. Minimizing the amplitude variation by optimizing the phases of the incident waves and taking only the phase-variation gives a good first guess for the optimization process.
  • Optimization of the best (e.g. 10) grating structures to get the optimum grating structure.

 

Measured intensity distribution Structure of the 8 beam FG (superimposed 2 and 4 beam FG)

Fig. 7: Measured intensity distribution of a 490GHz beam, split by the grating for STAR (left hand side). Contour levels are from 5% to 95% of the maximum peak intensity in steps of 10%. The theoretical (in parenthesis) and the measured efficiency is: (83.0%) 84.0%. And the surface topology of the unit cell using 13 Fourier coefficients for each direction (right hand side). Contour levels are from -1.5 to 1.5 radians in steps of .25 radians.

 

To show the power of the design process, we calculated some letters as Fourier gratings using also 13 coefficients per lateral dimension. Fig. 8 shows the unit cells of the letters (top row) and their calculated diffraction pattern. These grating structures would be very hard to manufacture as stepped gratings, but easily as Fourier gratings.

Letters as FG's
Fig. 8

Contents

Collimating Fourier gratings (CFG)

Using planar gratings, the LO source had to be reimaged to form a plane wave at the location of the grating. The new collimating Fourier grating is matched to the wave front curvature of the divergent beam and does not require an additional collimating element. Because of the smooth grating structure it is possible to produce the collimating Fourier gratings on a commercial numerically controlled milling machine.

LO-unit with a standard FG LO-unit with a CFG

Fig. 9: Comparison of a LO-unit using a standard flat Fourier grating with a collimating mirror and (left hand side) with LO-unit using the collimating Fourier grating.

Measured intensity distribution Structure of the unit cell

Fig. 10 : Measured intensity distribution of a 350GHz beam, split by the grating for the Desert-Star project of SORAL (left hand side). Contour levels are from 5% to 95% of the maximum peak intensity in steps of 10%. The surface topology of the unit cell using 8 complex Fourier coefficients for each direction is shown at the right hand side. Contour levels here are from 0 to 6 radians in steps of .5 radians. The unit cell ist streched in one direction by sqrt(3) to match the hexagonal arrangement.

 

The collimating Fourier gratings can as well be manufactured easily on our in house CNC milling machine. Fig. 11 shows a photograph of the grating for the Desert-Star project. The 4 corners are manufactured as referenc planes for easy repositioning of the grating.

Photograph of a CFG
Fig. 11

Contents

Contact

For Further informations contact: Dr. Urs Graf