Metal-mesh optical filter
Metal-mesh optical filters are optical filters made from stacks of metal meshes and dielectric. They are used as part of an optical path to filter the incoming light to allow frequencies of interest to pass while reflecting other frequencies of light.
Metal-mesh filters have many applications for use in the far infrared (FIR)<ref name=melo08> Arline M. Melo; Mariano A. Kornberg; Pierre Kaufmann; Maria H. Piazzetta; et al. (Nov 2008). "Metal mesh resonant filters for terahertz frequencies". Applied Optics. 47 (32): 6064–9. Bibcode:2008ApOpt..47.6064M. doi:10.1364/AO.47.006064. PMID 19002231. </ref> and submillimeter regions of the electromagnetic spectrum. These filters have been used in FIR and submillimeter astronomical instruments for over 4 decades,<ref name=ade06> Ade, Peter A. R.; Pisano, Giampaolo; Tucker, Carole; Weaver, Samuel (Jul 2006). Zmuidzinas, Jonas; Holland, Wayne S; Withington, Stafford; Duncan, William D (eds.). "A Review of Metal Mesh Filters" (PDF). Millimeter and Submillimeter Detectors and Instrumentation for Astronomy III. Proceedings of the SPIE. Millimeter and Submillimeter Detectors and Instrumentation for Astronomy III. 6275: 62750U. Bibcode:2006SPIE.6275E..0UA. doi:10.1117/12.673162. S2CID 7647444. </ref> in which they serve two main purposes: band-pass or low-pass filters are cooled and used to lower the noise equivalent power of cryogenic bolometers (detectors) by blocking excess thermal radiation outside of the frequency band of observation,<ref name=porterfield94> D. W. Porterfield; J. L. Hesler; R. Densing; E. R. Mueller; et al. (Sep 1994). "Resonant metal-mesh bandpass filters for the far infrared". Applied Optics. 33 (25): 6046–52. Bibcode:1994ApOpt..33.6046P. doi:10.1364/AO.33.006046. PMID 20936018. </ref> and band-pass filters can be used to define the observation band of the detectors. Metal-mesh filters can also be designed for use at 45° to split an incoming optical signal into several observation paths, or for use as a polarizing half-wave plate.<ref name=pisano08> Giampaolo Pisano; Giorgio Savini; Peter A. R. Ade; and Vic Haynes (2008). "Metal-mesh achromatic half-wave plate for use at submillimeter wavelengths". Applied Optics. 47 (33): 6251–6256. Bibcode:2008ApOpt..47.6251P. doi:10.1364/AO.47.006251. PMID 19023391. </ref>
Transmission line theory can be applied to metallic meshes to understand how they work and the overall light transmission properties of groups of metallic meshes grouped together.<ref name=ulrich67> R. Ulrich (Mar 1967). "Far infrared properties of metallic mesh and its complementary structure". Infrared Physics. 7 (1): 37–50. Bibcode:1967InfPh...7...37U. doi:10.1016/0020-0891(67)90028-0. </ref> Modeling the properties of these metallic meshes allows for reliable manufacture of filters with the desired transmission properties.
Theory
In 1967 Ulrich showed that the optical transmission properties of a metallic mesh can be modeled by considering the mesh to be a simple circuit element on a free space transmission line. To develop the theory of metallic meshes, he focused on the properties of two types of mesh structure: a metallic grid with square openings; and a grid of metallic squares supported on a thin dielectric substrate. Using the transmission line method, he then modeled the behavior of each of these meshes as either lumped inductance (square openings) or a lumped capacitance (free-standing squares). These two types of meshes are commonly referred to as inductive or capacitive meshes.<ref name="ade06"/><ref name = ulrich67 />
The theory developed by Ulrich to explain light transmission by metallic meshes makes a few assumptions and idealizations, which will be used here as well in explaining the theory. This theory is valid for thin meshes, i.e. <math>t<<a</math>, but the following equations assume that the grid is infinitely thin, the metallic parts are perfectly conducting, and the supporting dielectric film in the capacitive grids has no effect. The electromagnetic theory can then be applied to develop a model of an oscillating circuit on a transmission line model that explains the transmission properties of these meshes quite well as long as the wavelength of light is larger than the size of the metallic element (<math>\lambda > g</math>).<ref name = ulrich67 />
Electromagnetic theory
Electromagnetic theory of light can be used to describe how light incident on both capacitive and inductive metallic meshes will behave in transmission, reflection, and absorption.
Transmission and reflection
If an incident plane wave of electromagnetic radiation hits a metallic grid of either type perpendicular to its path it will scatter, and the only propagating parts will be the zeroth order reflected wave and the zeroth order transmitted wave.<ref name = ulrich67 /> The frequency of both of these electric fields will be equal, and the ratio of their amplitudes is <math>\Gamma(\omega)</math>, where <math>\Gamma</math> is the reflection coefficient, and <math>\omega = g/\lambda</math> is the normalized frequency. If we assume that the incident wave had unit amplitude, we can add the incident wave to the transmitted scattered wave to get the total amplitude of the transmitted wave, <math>\tau(\omega)</math>:
<math>\tau(\omega) = \left[1 + \Gamma(\omega)\right]</math>.
Since we are neglecting losses, the amplitude squared of the reflected and transmitted waves must equal unity:
<math>\left|\Gamma(w)\right|^2 + \left|\tau(\omega)\right|^2 = 1</math>.
Given these two relations, the phase of the reflection coefficient, <math>\phi_{\Gamma}(\omega)</math>, and the phase of the transmission coefficient <math>\phi_{\tau}(\omega)</math>can be simply related to the transmitted power, <math>\left|\tau(\omega)\right|^2</math>, which can be directly measured in experiments with metallic meshes.
<math>\sin^2 \phi_{\Gamma} = 1 - \left|\tau(\omega)\right|^2</math>
<math>\sin^2 \phi_{\tau} = \left|\tau(\omega)\right|^2</math>
Solving these equations lets us find the amplitude of the scattered wave in terms of the phases of the reflected and transmitted waves:
<math>\left|\Gamma(\omega)\right|^2 = \sin^2 \phi_{\Gamma}(\omega) = 1 - \sin^2 \phi_{\tau}(\omega)</math>.
The result of drawing <math>\Gamma(\omega)</math> vs. <math>\omega</math> in the complex plane is a unit half circle centered on the point <math>\left[Re(-1/2), Im(0)\right]</math> which is in the upper half-plane <math>\left(Im(\Gamma(\omega)) > 0\right)</math> for inductive grids and in the lower half-plane <math>\left(Im(\Gamma(\omega)) < 0\right)</math> for capacitive grids. At all frequencies <math>\omega</math> the transmitted and reflected waves are out of phase <math>\left(\phi_\tau(\omega) \ne \phi_\Gamma(\omega)\right)</math>.<ref name = ulrich67 />
Until now, the theory has been general—whether the mesh was inductive or capacitive has not been specified. Since <math>\tau(\omega)</math> and <math>\Gamma(\omega)</math> are independent of polarization, we can apply Babinet's principle to the capacitive and inductive grids. Concisely, Babinet's principle states that if we swap the metallic parts of a grid for the gaps, (i.e., make a complementary mesh), then the sum of the transmitted wave from the original structure and the structure's complement must equal the original incident wave.<ref name=born99> Max Born and Emil Wolf (1999). Principles of Optics (7th ed.). Cambridge University Press. ISBN 978-0-521-78449-8. </ref> Therefore, if we have complementary capacitive and inductive grids,
<math>\left[\tau_{ind} + \tau_{cap}\right] = 1</math>.
Given the relations between the reflected and transmitted waves found earlier, this means that the transmitted wave in an inductive grid is equal to the negative of the reflected wave in a capacitive grid and vice versa, and also that the transmitted powers for capacitive and inductive grids sum to unity for a unit incident wave.
<math>\tau_{ind}\left(\omega\right) = -\Gamma_{cap}(\omega)</math>
<math>\tau_{cap}\left(\omega\right) = -\Gamma_{ind}(\omega)</math>
<math>\left|\tau_{cap}(w)\right|^2 + \left|\tau_{ind}(\omega)\right|^2 = 1</math>.<ref name = ulrich67 />
Solving for the exact form of <math>\tau_{cap}\left(\omega\right)</math> or <math>\tau_{ind}\left(\omega\right)</math> requires solving Maxwell's equations on the grids, which for the general case can only be solved numerically. However, in an inductive grid the metal is continuous, and hence DC currents can exist. Considering the limiting case of <math> \omega \rightarrow 0</math>, the inductive grid must reflect the entire incident wave<ref name = ulrich67 /> because of the boundary conditions for the electric field at the surface of a conductor.<ref name=griffiths89> David J. Griffiths (1989). Introduction to Electrodynamics (2nd ed.). Prentice Hall. </ref> The relations derived above therefore show that a capacitive mesh will transmit the entire incident wave in this case.
<math>\tau_{ind}\left(\omega \rightarrow 0\right) = 0</math>
<math>\tau_{cap}\left(\omega \rightarrow 0\right) = 1</math>
Because the grids are complements of each other, these equations show that a capacitive mesh is a low-pass filter and an inductive mesh is a high-pass filter.<ref name = ulrich67 />
Absorption
Up until now, the theory has only been considering the ideal case where the grids are infinitely thin and perfectly conducting. In principle grids with finite dimensions could also absorb some of the incident radiation either through ohmic losses or losses in the dielectric supporting material.
Assuming that the skin depth of the metal being used in the grids is much smaller than the thickness of the grid, the real part of the surface impedance of the metal is <math>\rho = 1/\delta \sigma</math> where <math>\sigma</math> is the conductivity of the metal and <math>\delta</math> is the skin depth of the metal. With a reflected wave <math>\Gamma(\omega)</math>, the change in the magnetic field amplitude across the grid is <math>2\Gamma(\omega)</math> because of surface currents on both sides of the grid. The average surface currents on both sides of the grid are <math>\bar{J} = \Gamma(\omega)* c/4\pi</math>.<ref name = ulrich67 />
Given the average surface current and the surface impedance, we could calculate the power dissipated as <math>P_D = 2\rho \bar{J}^2</math>. However, because the actual extent of the metal in the grids is different between the capacitive and inductive grids and a flat sheet of metal, we need to introduce a factor <math>\eta</math> which is the ratio of the area the grid to that of a flat sheet. For capacitive grids, <math>\eta = g/2a</math> and for inductive grids <math>\eta = 1/(1-2a/g)</math>. This modifies the power dissipated to be <math>P_D = 2\rho\eta\bar{J}^2</math>. Using the definition of skin depth, the unitless absorptivity, <math>A = P_d/P_o</math> where <math>P_o</math> is the incident power, of the grid is
<math>A = \left|\Gamma\right|^2 2 \rho \eta = \left|\Gamma\right|^2 \eta \left(\frac{c}{\lambda \sigma}\right)^{1/2}</math>.<ref name = ulrich67 />
For microwave and infrared radiation incident on copper, this unitless absorptivity comes out to be <math>10^{-4}</math> to <math>10^{-2}</math>, which means that the initial assumption that absorption could be ignored in this ideal model was a good one. The dielectric losses can likewise be ignored.<ref name = ulrich67 />
Comparison to measurements
For single layer metallic grids, the simple theory Ulrich laid out works quite well. The functions <math>\left|\tau_{cap}(\omega)\right|^2</math> and <math>\left|\tau_{ind}(\omega)\right|^2</math> can be determined by measuring the transmission through the filter, and the phases <math>\phi_{cap}\left(\omega\right)</math> and <math>\phi_{ind}\left(\omega\right)</math> can be measured by setting two identical grids variable distances apart and measuring the interference maximum of <math>\Gamma(\omega)\phi(\omega)</math> as a function of separation. Measurements of very thin nearly ideal grids show the expected behavior and have very low absorptive loss.<ref name = ulrich67 />
In order to build filters out of metallic meshes with the desired properties, it is necessary to stack many metallic meshes together, and while the simple electromagnetic theory laid out above works well for one grid, it becomes more complicated when more than one element is introduced. However, these filters can be modeled as elements in a transmission line, which has easily calculable transmission properties.<ref name = ade06 /><ref name = ulrich67 />
Transmission line model
A transmission line model of metallic meshes is easy to work with, flexible, and is readily adapted for use in electronic modeling software. It not only handles the case of a single metallic grid, but is easily extended to many stacked grids.
Theoretical model
Under the conditions of normal incidence and <math>\omega<1</math> the electric field across a metallic grid is continuous, but the magnetic field is not,<ref name=born99 /> so a transmission line with an admittance <math>2Y(\omega)</math> between the two lines can be used to model the transmission and reflection from a metallic filter. If, for example, three identical grids were stacked, then there would be three admittance shunts in parallel across the transmission line. Using simple transmission line theory, the reflection coefficient <math>\Gamma(\omega)</math> and transmission coefficient <math>\tau(\omega)</math> are calculated to be
<math>\Gamma(\omega) = \frac{- Y(\omega)}{1 + Y(\omega)}</math>
<math>\tau(\omega) = \frac{1}{1 + Y(\omega)}</math>
which of course satisfy the original relation between the transmission and reflection coefficients:
<math>\tau(\omega) = \left[1 + \Gamma(\omega)\right]</math>.
In a lossless circuit, the admittance becomes a purely imaginary susceptance, <math>Y\left(\omega\right) = iB(\omega)</math> where <math> B\left(\omega\right)</math> is a real function of <math>\omega</math>. Because of the complementary nature of the grids, we also know that <math>B_{ind}\left(\omega\right) B_{cap}(\omega) = -1</math>.<ref name=ulrich67 />
To calculate the behavior of an ideal metallic grid, only <math>B\left(\omega\right)</math> needs to be found. The standard approach is not to characterize the equivalent circuit by <math>B\left(\omega\right)</math>, but instead to parameterize it with values of <math>L</math>, <math>C</math>, and <math>R</math> which duplicate the transmission properties of the filters. At low frequencies, a reasonable model is to replace the shunt in the transmission line with a capacitor of value <math>2C</math> for capacitive meshes and an inductor of value <math>L/2</math> for inductive meshes, where for complementary grids <math>L_{ind} = C_{cap}</math>. However, at high frequencies this model fails to reflect the behavior of real metallic meshes correctly. The measured transmissions as <math>\omega \rightarrow 1</math> are
<math>\tau_{cap} \left(\omega \rightarrow 1 \right) = 0</math>
<math>\tau_{ind} \left(\omega \rightarrow 1 \right) = 1</math>.<ref name=ulrich67 />
The behavior of the transmission in the two limiting cases can be replicated with the transmission line model by adding an extra element. In addition, losses can be taken into account by adding another resistance <math>R</math>. At resonance <math>\left(\omega = \omega_o\right)</math>, the impedance of capacitors and inductors are <math>Z_o = i\omega L = 1/i \omega C</math>. Typically, <math>Z_o</math> and <math>\omega_o</math> have to be measured based on transmission properties of the grids, and both depend on the parameter <math>a/g</math>. The <math>R</math> included in the 2-element equivalent circuit is consistent with the earlier calculation of absorptivity, which gives <math>R = \eta/2 \left(\frac{c}{\lambda \sigma}\right)^{1/2}</math>. The following table summarizes all the parameters to go from equivalent circuit parameters to expected reflection and transmission coefficients.<ref name=ulrich67 />
| Capacitive circuit | Inductive Circuit | |
|---|---|---|
| Normalized Impedance <math>Z_o\left(\omega_o\right)</math> | <math>Z_o = i\omega L = 1/\left(i \omega C\right)</math> | |
| Generalized frequency <math>\Omega\left(\omega\right)</math> | <math>\Omega\left(\omega\right) = \left(\omega / \omega_o \right) - \left( \omega_o / \omega \right) = \left( \lambda_o / \lambda \right) - \left( \lambda / \lambda_o \right)</math> | |
| Normalized admittance <math>Y\left(\omega\right)</math> | <math>\frac{1}{1+iZ_o \Omega}</math> | <math>\frac{1}{1-iZ_o / \Omega}</math> |
| \Gamma\left(\omega\right)\right|^2</math> | <math>\frac{1}{(1+R)^2+Z_o^2\Omega^2}</math> | <math>\frac{1}{(1+R)^2+Z_o^2/\Omega^2}</math> |
| \tau\left(\omega\right)\right|^2</math> | <math>\frac{R^2 + Z_o^2\Omega^2}{(1+R)^2+Z_o^2\Omega^2}</math> | <math>\frac{R^2 + Z_o^2/\Omega^2}{(1+R)^2+Z_o^2/\Omega^2}</math> |
| Reflected phase <math>\phi_{\Gamma}\left(\omega\right)</math> | <math>\pi - \arctan{\left(\frac{Z_o\Omega}{(1+R)}\right)}</math> | <math>\pi + \arctan{\left(\frac{Z_o}{(1+R)\Omega}\right)}</math> |
| Transmitted phase <math>\phi_{\tau}\left(\omega\right)</math> | <math>\arctan{\left(\frac{Z_o\Omega}{R(1+R)+Z_o^2\Omega^2}\right)}</math> | <math>-\arctan{\left(\frac{Z_o/\Omega}{R(1+R)+Z_o^2/\Omega^2}\right)}</math> |
| Absorptivity <math>A\left(\omega\right)</math> | <math>2R\left|\Gamma\right|^2</math> | |
The real power in this model is it allows prediction of the transmission properties of many metallic grids stacked together with spacers to form interference filters. Stacks of capacitive grids make a low-pass filter with a sharp frequency cutoff above which transmission is almost zero. Likewise, stacks of inductive grids make a high-pass filter with a sharp frequency cutoff below which transmission is almost zero. Stacked inductive and capacitive meshes can be used to make band-pass filters.<ref name=ade06 />
Comparison to measurements
The transmission line model gives the expected first-order transmission of the stacked metal mesh filters; however, it cannot be used to model transmission of light that is incident at an angle, loss in the supporting dielectric materials, or the transmission properties when <math>\lambda < g</math> due to diffraction. To model those effects, scientists have used a cascade scattering matrix approach to model dielectric loss, and other modeling tools such as High Frequency Structure Simulator and Floquet mode analysis.<ref name = ade06 />
Manufacture
The manufacture of metal-mesh filters starts with photolithography of copper on a substrate, which allows fine control over the parameters <math>a</math>, <math>g</math>, and <math>t</math>. The metallic grids are made of thin copper film on top of a dielectric substrate such as mylar or polypropylene. The copper is <math>\approx .4\mu m</math> thick, and the dielectric ranges from <math>.9\mu m</math> to <math> 1.5\mu m</math>.<ref name = ade06 />
There are two ways to create a multi-layer metal-mesh filter. The first is to suspend the separate layers in supporting rings with a small gap which is either filled with air or under vacuum between the layers. However, these filters are mechanically delicate. The other way to build a multi-layer filer is to stack sheets of dielectric between the layers of metallic mesh and hot press the whole stack together. This results in a filter that is one solid piece. Hot pressed filters are mechanically robust and when impedance matched to vacuum show a pass-band fringe due to Fabry-Perot interference in the underlying dielectric material.<ref name = ade06 />
Use in experiments
These filters have been used in FIR and submillimeter astronomical instruments for over 4 decades,<ref name = ade06 /> in which they serve two main purposes: band-pass or low-pass filters are cooled and used to lower the noise equivalent power of cryogenic bolometers by blocking excess thermal radiation outside of the frequency band of observation,<ref name = porterfield94 /> and band-pass filters can be used to define the observation band of the detectors. Metal-mesh filters can also be designed for use at 45° to split an incoming optical signal into several observation paths, or for use as a polarizing half-wave plate.<ref name=pisano08 />
-
Low-pass filter at 4 Kelvin used for blocking excess thermal radiation installed in the South Pole Telescope receiver.
-
Low-pass metal-mesh filters at 250mK used to define the upper edge of the detector observation band in the South Pole Telescope receiver.
