US20070201136A1  Thin Film Interference Filter and Bootstrap Method for Interference Filter Thin Film Deposition Process Control  Google Patents
Thin Film Interference Filter and Bootstrap Method for Interference Filter Thin Film Deposition Process Control Download PDFInfo
 Publication number
 US20070201136A1 US20070201136A1 US11/684,779 US68477907A US2007201136A1 US 20070201136 A1 US20070201136 A1 US 20070201136A1 US 68477907 A US68477907 A US 68477907A US 2007201136 A1 US2007201136 A1 US 2007201136A1
 Authority
 US
 United States
 Prior art keywords
 reflectance
 δ
 method
 cos
 wavelength
 Prior art date
 Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
 Abandoned
Links
 239000010409 thin films Substances 0 title claims abstract description 25
 238000004886 process control Methods 0 title claims description 5
 239000010408 films Substances 0 claims abstract description 112
 239000010410 layers Substances 0 claims abstract description 84
 230000001747 exhibited Effects 0 claims abstract description 7
 230000000875 corresponding Effects 0 claims abstract description 3
 230000003287 optical Effects 0 claims description 38
 230000003466 anticipated Effects 0 claims description 8
 238000000151 deposition Methods 0 claims description 8
 238000005259 measurements Methods 0 claims description 8
 238000004364 calculation methods Methods 0 description 30
 239000002609 media Substances 0 description 14
 238000001228 spectrum Methods 0 description 11
 239000000463 materials Substances 0 description 10
 238000000034 methods Methods 0 description 5
 238000004458 analytical methods Methods 0 description 2
 238000005365 production Methods 0 description 2
 239000000654 additives Substances 0 description 1
 239000000562 conjugates Substances 0 description 1
 238000010291 electrical methods Methods 0 description 1
 238000004519 manufacturing process Methods 0 description 1
 238000010297 mechanical methods and processes Methods 0 description 1
 238000005498 polishing Methods 0 description 1
Images
Classifications

 G—PHYSICS
 G01—MEASURING; TESTING
 G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
 G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using infrared, visible or ultraviolet light
 G01N21/84—Systems specially adapted for particular applications
 G01N21/8422—Investigating thin films, e.g. matrix isolation method

 G—PHYSICS
 G01—MEASURING; TESTING
 G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
 G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using infrared, visible or ultraviolet light
 G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
 G01N21/55—Specular reflectivity

 G—PHYSICS
 G01—MEASURING; TESTING
 G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
 G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using infrared, visible or ultraviolet light
 G01N21/84—Systems specially adapted for particular applications
 G01N21/8422—Investigating thin films, e.g. matrix isolation method
 G01N2021/8438—Mutilayers

 G—PHYSICS
 G02—OPTICS
 G02B—OPTICAL ELEMENTS, SYSTEMS, OR APPARATUS
 G02B5/00—Optical elements other than lenses
 G02B5/20—Filters
 G02B5/28—Interference filters
Abstract
A thin film interference filter system includes a plurality of stacked films having a determined reflectance; a modeled monitor curve; and a topmost layer configured to exhibit a wavelength corresponding to one of the determined reflectance or the modeled monitor curve. The topmost layer is placed on the plurality of stacked films and can be a lowindex film such as silica or a high index film such as niobia.
Description
 This application claims benefit of U.S. Provisional Patent Application, Ser. No. 60/609,406, entitled “Bootstrap Method for Interference Filter Thin Film Deposition Process Control,” filed Sep. 13, 2004.
 Research described in this application was sponsored by the United States Air Force Research Laboratory grant number F226150026059.
 The present invention relates to thin film optical devices. More particularly, the invention relates to complex interference filters.
 Highly accurate optical interference filters can be manufactured using thin film deposition processes. These optical interference filters are used for multivariate optical computing, multiplebandpass, and the like and can exhibit complex optical spectra defined over a range of wavelengths. These filters are typically constructed by depositing alternating layers of transparent materials where one layer possesses a much larger refractive index relative to the other layer. Theoretically, the proper choice of composition, thickness and quantity of layers could result in a device with any desired transmission spectrum.
 Among the simplest devices is the single cavity bandpass filter; i.e., the thinfilm form of an etalon. This device consists of three sets of layers. The first stack is a dielectric mirror, a next thicker layer forms a spacer, and a second stack forms another dielectric mirror. The mirror stacks are typically fabricated by depositing alternating transparent materials that have an optical thickness that is one quarter of the optical wavelength of light.
 To achieve theoretical optical performance, each layer must possess a precise and specific physical thickness and refractive index. Any nonuniformity in the deposition of the layers can affect the spectral placement and transmission or reflection characteristics of the device. A design that requires very tight manufacturing tolerances over large substrate areas could result in the costly rejection of many devices. Given these manufacturing limits, it would be desirable to analyze the devices after construction and alter the devices that do not meet a predetermined optical transmission or reflection specification by some electrical or mechanical means. For example, if the peak transmission wavelength of a manufactured optical bandpass cavity filter was slightly out of tolerance, it would be desirable to have a mechanism or process for shifting the peak back to the desired spectral location. It is also desirable that the optical filters have precise rejection bands and passbands that are electrically or mechanically selectable.
 Mechanical methods of achieving a variable transmission spectrum device are well known. This includes changing a prism or grating angle, or altering the optical spacing between mirrors of an etalon. To overcome the performance, size and cost disadvantages of using mechanical schemes, many have conceived of electrical methods for varying a transmission spectrum. For example, U.S. Pat. No. 5,150,236, issued Sep. 22, 1992 to Patel, discloses a tunable liquid crystal etalon filter. The liquid crystal fills the space between dielectric mirrors. Electrodes on the mirrors are used to apply an electric field, which changes the orientation of the liquid crystal that changes the optical length for tuning. The change in the optical length corresponds to a change in the location of the passband. U.S. Pat. No. 5,103,340, issued Apr. 7, 1992 to Dono et al., discloses piezoelectric elements placed outside the optical path that are used to change the spacing between cascaded cavity filters. Furthermore, U.S. Pat. No. 5,799,231, issued Aug. 25, 1998 to Gates et al., discloses a variable index distributed mirror. This is a dielectric mirror with half of the layers having a variable refractive index that is matched to other layers. Changing the applied field increases the index difference that increases the reflectance. The mathematics that describes the transmission characteristics of multilayer films composed of electrooptic and dielectric materials are well known.
 Another electrically actuated thin film optical filter uses a series of crossed polarizers and liquid crystalline layers that allow electrical controls to vary the amount of polarization rotation in the liquid by applying an electric field in such a way that some wavelengths are selectively transmitted. However, these electrically actuated thin film optical filters have the characteristic that the light must be polarized and that the frequencies of light not passed are absorbed, not reflected. Another electrically actuated thin film optical device is the tunable liquid crystal etalon optical filter. The tunable liquid crystal etalon optical filter uses a liquid crystal between two dielectric mirrors.
 The common cavity filter, such as the etalon optical filter, is an optical filter with one or more spacer layers that are deposited in the stack and define the wavelength of the rejection and pass bands. The optical thickness of the film defines the placement of the passband. U.S. Pat. No. 5,710,655, issued Jan. 20, 1998 to Rumbaugh et al., discloses a cavity thickness compensated etalon filter.
 In the tunable liquid crystal etalon optical filter, an electric field is applied to the liquid crystal that changes the optical length between the two mirrors so as to change the passband of the etalon. Still another tunable optical filter device tunes the passband by using piezoelectric elements to mechanically change the physical spacing between mirrors of an etalon filter.
 Bulk dielectrics are made by subtractive methods like polishing from a larger piece; whereas thin film layer are made by additive methods like vapor or liquid phase deposition. A bulk optical dielectric, e.g., greater than ten microns, disposed between metal or dielectric mirrors suffers from excessive manufacturing tolerances and costs. Moreover, the bulk material provides unpredictable, imprecise, irregular, or otherwise undesirable passbands. These electrical and mechanical optical filters disadvantageously do not provide precise rejection bands and passbands that are repeatably manufactured.
 In an attempt to avoid some of the foregoing problems, modeling of interference filters can be conducted during online fabrication with insitu optical spectroscopy of the filter during deposition. The current state of the art for online correction of the deposition involves fitting the observed spectra to a multilayer model composed of “ideal” films based on a model for each film. The resulting model spectra are approximations of the actual spectra. To use reflectance as an example: the measured reflectance of a stack of films can be approximately matched to a theoretical reflectance spectrum by modeling. Layers remaining to be deposited can then be adjusted to compensate for errors in the film stack already deposited, provided the film stack has been accurately modeled. However, films vary in ways that cannot be readily modelled using any fixed or simple physical model. Heterogeneities in the films that cannot be predicted or compensated by this method cause the observed spectra to deviate more and more from the model. This makes continued automatic deposition very difficult; complex film stacks are therefore very operatorintensive and have a high failure rate. To improve efficiency in fabrication, laboratories that fabricate these stacks strive to make their films as perfectly as possible so the models are as accurate as possible.
 As outlined above, many thin films are usually designed in a stack to produce complex spectra and small variations in deposition conditions make it difficult to accurately model in situ film spectra for feedback control of a continuous deposition process because it is practically impossible to obtain full knowledge of the detailed structure of the stack from reflectance, transmittance, ellipsometry, mass balance or other methods. Thus, a thin film interference filter is needed that is less difficult to manufacture, which will address varying refractive indices of thin films and varying absorptions with deposition parameters.
 In general, the present invention is directed to a layered, thin film interference filter and related bootstrap methods. A bootstrap method according to one aspect of the invention permits a user to focus on a single layer of a film stack as the layer is deposited to obtain an estimate of the properties of the stack. Although the single layer model is a guideline and not a basis for compensating errors, only the mostrecentlydeposited layer—and not the alreadydeposited film stack—need be modeled according to an aspect of the present invention. Thus, the user can neglect deviations of the stack from ideality for all other layers. The singlelayer model can then be fit exactly to the observed spectra of the film stack at each stage of deposition to allow accurate updating of the remaining film stack for continued deposition.
 The present invention works with any type of films, whether absorbing or nonabsorbing, and regardless of whether the control of the deposition conditions are stateoftheart or not. The methods of the invention are relatively straightforward, and a resultant thin film interference filter is economical to produce and use.
 According to a particular aspect of the invention, a method using experimental measurements to determine reflectance phase and complex reflectance for arbitrary thin film stacks includes the steps of determining reflectance of a stack of a plurality of films before depositing a topmost layer; considering a modeled monitor curve for a wavelength of a highindex layer; and discarding a plurality of monitor curves without maxima in their reflectance during the topmost layer deposition. In this aspect of the invention, the topmost layer can be a niobia layer.
 The exemplary method can also include the steps of determining an anticipated standard deviation in φ_{k }for a plurality of monitor wavelengths in the niobia layer and discarding any with σ greater than 0.9 degrees. Another step in this aspect is computing expected error in δ for wavelengths with σ less than 0.9 degrees at a target thickness of the niobia layer. When no wavelengths have an error less than 0.9 degrees, a further step according to this method is to proceed with a full model deposition of the niobia layer. When there are no maxima in the reflectance of each of the plurality of monitor curves, another step according to the exemplary method is to use only the modeled monitor curve during the topmost layer deposition.
 The method according to this aspect of the invention can also include the step of computing a value of δ for all wavelengths based on a value calculated for the monitor wavelength.
 The method according to this aspect of the invention can further include the step of computing two possible values of phase angle for each wavelength other than the monitor wavelength.
 Additional steps according to the exemplary method include using information extracted from the model for r_{k }at each wavelength and the computed best value of δ, and computing an estimated standard deviation of phase at all wavelengths except the monitor.
 The method may further include the steps of using the computed phase closest to the model phase for r_{k }at each wavelength, measured R_{f }and R_{k }values and the computed best value of δ, and computing the estimated standard deviation of phase at all wavelengths for which the magnitude of r_{k }was estimated other than the monitor.
 Further steps according to the exemplary method include the steps of determining if a phase error estimate is less than about 1.3 degrees and averaging calculated and modeled reflectance and phase values to obtain a new value for use in subsequent modeling at that wavelength.
 In yet another aspect according to the exemplary method, the topmost layer can be a silica film. Accordingly, the method can include the step of replacing the magnitude of the amplitude reflectance at each wavelength with √{square root over (R_{k})} whenever measuring a latest depositing silica film having an intensity reflectance greater than 9%. The method can also include the step of determining if a phase error estimate is less than about 1.3 degrees when the magnitude of the amplitude reflectance at each wavelength has been replaced with √{square root over (R_{k})} and averaging calculated and modeled reflectance and phase values to obtain a new value for use in subsequent modeling at that wavelength.
 In yet another aspect of the invention, a method for correcting thin film stack calculations for accurate deposition of complex optical filters can include the steps of determining phase angle φ_{k }at a monitor wavelength from r′_{k} and Rk using a first equation expressed as:
$\mathrm{cos}\text{\hspace{1em}}\left(\pm \text{\hspace{1em}}{\varphi}_{k}\right)=\frac{{\uf603{r}_{k}^{\prime}\uf604}^{2}\left(1+{R}_{k}{\uf603{r}_{2}\uf604}^{2}\right){\uf603{r}_{2}\uf604}^{2}{R}_{k}}{2\text{\hspace{1em}}\uf603{r}_{2}\uf604\sqrt{{R}_{k}}\left(1{\uf603{r}_{k}^{\prime}\uf604}^{2}\right)};$  and estimating r′_{k }using a second equation
${r}_{k}^{\prime}=\frac{{r}_{k}{r}_{2}}{1{r}_{2}{r}_{k}};$
and obtaining a value for phase at a monitor wavelength.  Still another aspect of the invention includes a method for automated deposition of complex optical interference filters including the steps of determining from a measurement of intensity reflectance at a topmost interface a phase angle φ at an interface k according to an equation expressed as:
$\mathrm{cos}\text{\hspace{1em}}\left({\varphi}_{k}\right)=\left(\frac{A\left(1+{r}_{2}^{2}\right)\mathrm{sin}\text{\hspace{1em}}\left(\delta \right)\pm B\text{\hspace{1em}}\mathrm{cos}\left(\delta \right)}{C}\right)$ $\begin{array}{c}A={R}_{f}+{r}_{2}^{4}\left({R}_{f}{R}_{k}\right){R}_{k}+\\ 2\text{\hspace{1em}}{r}_{\text{\hspace{1em}}2}^{\text{\hspace{1em}}2}\left(\left(1{R}_{\text{\hspace{1em}}f}\right)\left(1+{R}_{\text{\hspace{1em}}k}\right)\mathrm{cos}\text{\hspace{1em}}\left(2\text{\hspace{1em}}\delta \right)\left(1{R}_{\text{\hspace{1em}}f}{R}_{\text{\hspace{1em}}k}\right)\right)\end{array}$ $B=\sqrt{D\left(1+{r}_{2}^{12}\right)+F\left({r}_{2}^{2}+{r}_{2}^{10}\right)+G\left({r}_{2}^{4}+{r}_{2}^{8}\right)+{\mathrm{Hr}}_{2}^{6}}$ $C=\mathrm{sin}\text{\hspace{1em}}\left(\delta \right)\left(4\text{\hspace{1em}}{r}_{2}\left(1{R}_{f}\right){R}_{k}^{1/2}\left(2\text{\hspace{1em}}{r}_{2}^{2}\mathrm{cos}\text{\hspace{1em}}\left(2\text{\hspace{1em}}\delta \right)1{r}_{2}^{4}\right)\right)$ $D={\left({R}_{f}{R}_{k}\right)}^{2}$ $F=2\left(\begin{array}{c}{R}_{k}\left(2+{R}_{k}\right)+{R}_{f}^{2}\left(1+2\text{\hspace{1em}}{R}_{k}\right)+2\text{\hspace{1em}}{R}_{f}\left(15\text{\hspace{1em}}{R}_{k}+{R}_{k}^{2}\right)\\ 2\left(1{R}_{f}\right)\left(1{R}_{k}\right)\left({R}_{f}+{R}_{k}\right)\mathrm{cos}\text{\hspace{1em}}\left(2\text{\hspace{1em}}\delta \right)\end{array}\right)$ $\begin{array}{c}G=64\text{\hspace{1em}}{R}_{f}5\text{\hspace{1em}}{R}_{f}^{2}4\text{\hspace{1em}}{R}_{k}+38\text{\hspace{1em}}{R}_{f}{R}_{k}4\text{\hspace{1em}}{R}_{f}^{2}{R}_{k}5\text{\hspace{1em}}{R}_{k}^{2}\\ 4\text{\hspace{1em}}{R}_{f}{R}_{k}^{2}6\text{\hspace{1em}}{R}_{f}^{2}{R}_{k}^{2}+8\left(1{R}_{f}^{2}\right)\left(1{R}_{k}^{2}\right)\mathrm{cos}\text{\hspace{1em}}\left(2\text{\hspace{1em}}\delta \right)\\ 2{\left(1{R}_{f}\right)}^{2}{\left(1{R}_{k}\right)}^{2}\mathrm{cos}\text{\hspace{1em}}\left(4\text{\hspace{1em}}\delta \right)\end{array}$ $H=4\left(\begin{array}{c}\begin{array}{c}3+2\text{\hspace{1em}}{R}_{f}^{2}10\text{\hspace{1em}}{R}_{f}{R}_{k}+2\text{\hspace{1em}}{R}_{k}^{2}+3\text{\hspace{1em}}{R}_{f}^{2}{R}_{k}^{2}\\ 2\left(1{R}_{\text{\hspace{1em}}f}\right)\left(1{R}_{\text{\hspace{1em}}k}\right)\left(2+{R}_{f}+{R}_{k}+2\text{\hspace{1em}}{R}_{f}{R}_{k}\right)\mathrm{cos}\text{\hspace{1em}}\left(2\text{\hspace{1em}}\delta \right)+\end{array}\\ {\left(1{R}_{f}\right)}^{2}{\left(1{R}_{k}\right)}^{2}\mathrm{cos}\text{\hspace{1em}}\left(4\text{\hspace{1em}}\delta \right)\end{array}\right)\text{}$  According to this exemplary method, a process control for a deposition system is bootstrapped by detaching the deposition system from all but the topmost interface.
 The method can further include the step of validating two resultant solutions according to the expression:
${R}_{f}=\frac{2\text{\hspace{1em}}{r}_{2}^{2}+{R}_{k}\left(1+{r}_{2}^{4}\right)+2\text{\hspace{1em}}{r}_{2}Q}{1+{r}_{2}^{4}+2\text{\hspace{1em}}{r}_{2}^{2}{R}_{k}+2\text{\hspace{1em}}{r}_{2}Q}$ $\begin{array}{c}Q={r}_{2}^{2}{R}_{k}^{1/2}\mathrm{cos}\text{\hspace{1em}}\left(2\text{\hspace{1em}}\delta +{\varphi}_{k}\right)+{R}_{k}^{1/2}\mathrm{cos}\text{\hspace{1em}}\left(2\text{\hspace{1em}}\delta {\varphi}_{k}\right)\\ {r}_{2}\left(1+{R}_{k}\right)\mathrm{cos}\text{\hspace{1em}}\left(2\text{\hspace{1em}}\delta \right)\left(1+{r}_{2}^{2}\right){R}_{k}^{1/2}\mathrm{cos}\text{\hspace{1em}}\left({\varphi}_{k}\right)\end{array}$  The method can also include the step of averaging calculated and modeled reflectance and phase values to obtain a new value to be used in all future modeling at a given wavelength.
 According to another aspect of the invention, a thin film interference filter system includes a plurality of stacked films having a determined reflectance; a modeled monitor curve; and a topmost layer configured to exhibit a wavelength corresponding to one of the determined reflectance or the modeled monitor curve, the topmost layer being disposed on the plurality of stacked films. The topmost layer according to this aspect can be a lowindex film such as silica or a high index film such as niobia.
 Other aspects and advantages of the invention will be apparent from the following description and the attached drawings, or can be learned through practice of the invention.
 A full and enabling disclosure of the present invention, including the best mode thereof to one of ordinary skill in the art, is set forth more particularly in the remainder of the specification, including reference to the accompanying figures in which:

FIG. 1 is a schematic, cross sectional view of a stack of films according to an aspect of the invention; 
FIG. 2 is a schematic, cross sectional view of a stack of films similar toFIG. 1 showing an internal interface in the stack of films according to another aspect of the invention; 
FIG. 3 is a schematic, cross sectional view of a stack of films according to another aspect of the invention similar toFIG. 1 but showing amplitude reflectance of a buried interface; 
FIG. 4 is a representation of regions of permissible values of Rk and Rmax, particularly showing lowest values of σ_{φ}k for silica according to an aspect of the invention; 
FIG. 5 is similar toFIG. 4 but for niobia according to another aspect of the invention; 
FIG. 6 is a histogram for silica as inFIG. 4 ; 
FIG. 7 is a histogram for niobia as inFIG. 5 ; 
FIG. 8 is a perspective plot showing an error logarithm versus Rmax and δ in accordance with an aspect of the invention; and 
FIG. 9 is a histogram ofFIG. 8 data.  Detailed reference will now be made to the drawings in which examples embodying the present invention are shown. Repeat use of reference characters in the drawings and detailed description is intended to represent like or analogous features or elements of the present invention.
 The drawings and detailed description provide a full and detailed written description of the invention and the manner and process of making and using it, so as to enable one skilled in the pertinent art to make and use it. The drawings and detailed description also provide the best mode of carrying out the invention. However, the examples set forth herein are provided by way of explanation of the invention and are not meant as limitations of the invention. The present invention thus includes modifications and variations of the following examples as come within the scope of the appended claims and their equivalents.
 Turning now to the figures,
FIG. 1 shows a thin film interference filter 10, which broadly includes a substrate 12 upon which a stack of films 14AX is deposited (where x represents a theoretically infinite number of film layers). As shown, the last (alternatively, final, top or topmost) deposited film is designated by the alphanumeral 14A while previously deposited or lower level films are designated 14BX. An incoming ray 18 is shown inFIG. 1 being reflected at an interface 16 (also referred to herein as top or top surface and when mathematically referenced as k). The reflected ray is designated by the number 20. For simplicity, any contribution from multiple incoherent reflections in the substrate 12 is ignored in the following discussion and only reflections with respect to the film stack 14AX are described.  Typically, reflectance of the top surface 16 is obtained using a matrix calculation that is in turn built from the characteristic matrices of each of the preceding films 14AX. As shown in
FIG. 1 , a computed value of an electric field amplitude reflectance r_{k }is obtained by optimizing the thickness of the topmost film 14A to provide the best fit over the fall spectrum consistent with an understanding of the existing film stack 14BX. The calculated value of r_{k }can be considered an estimate of the actual value of r_{k }exhibited by the film stack 14AX. Since r_{k }is a complex value, it cannot be measured directly.  Normal Matrxtype Calculations:
 For standard calculations, the complex reflectance of a stack of films is computed using the admittances of the incident medium (often air), and the first interface of the stack. For s and p polarization, this reflectance is:
$\begin{array}{cc}{r}_{s}^{\left({f}_{1}\right)}=\frac{{\eta}_{s}^{\left(\mathrm{inc}\right)}{\eta}_{s}^{\left({f}_{1}\right)}}{{\eta}_{s}^{\left(\mathrm{inc}\right)}+{\eta}_{s}^{\left({f}_{1}\right)}}\text{}{r}_{p}^{\left({f}_{1}\right)}=\frac{{\eta}_{p}^{\left(\mathrm{inc}\right)}{\eta}_{p}^{\left({f}_{1}\right)}}{{\eta}_{p}^{\left(\mathrm{inc}\right)}+{\eta}_{p}^{\left({f}_{1}\right)}}& 6.\end{array}$
where the subscript s or p indicates for s or ppolarized light, η is a complex admittance, the superscript “inc” indicates the incident medium and (f_{1}) indicates the first interface the light strikes when coming from the incident medium.  For s and p polarized light, the admittances of the incident media are written:
$\begin{array}{cc}{\eta}_{s}^{\left(\mathrm{inc}\right)}=\sqrt{{\varepsilon}_{\mathrm{inc}}{\varepsilon}_{\mathrm{inc}}{\mathrm{sin}\text{\hspace{1em}}}^{2}\left({\theta}_{\mathrm{inc}}\right)}\text{}{\eta}_{p}^{\left(\mathrm{inc}\right)}=\frac{{\varepsilon}_{\mathrm{inc}}}{\sqrt{\text{\hspace{1em}}{\varepsilon}_{\mathrm{inc}}{\varepsilon}_{\mathrm{inc}}\text{\hspace{1em}}{\mathrm{sin}}^{2}\text{\hspace{1em}}(\text{\hspace{1em}}{\theta}_{\mathrm{inc}})}}& 7.\end{array}$
where ∈ indicates a complex dielectric constant, and θ_{inc }is the angle of incidence in the medium of incidence.  The problematic part of the calculation is how to express the admittance of the initial interface. The matrix calculation proceeds by relating the admittance of the initial interface to that of the second interface, the admittance of the second to the third, etcetera, through a series of 2×2 matrices, until the calculation is related to the final interface. At the final interface, the admittance (ratio of magnetic to electric fields) is equal to the admittance of the exit medium, which is simple to compute because there is only a single ray (the transmitted ray), rather than rays propagating in two different directions.
 For a single layer stack the admittance of the initial interface is related to the exit medium admittance according to the following equations:
$\begin{array}{cc}\begin{array}{c}\left(\begin{array}{c}{E}_{x}^{\left({f}_{1}\right)}\\ {H}_{y}^{\left({f}_{1}\right)}\end{array}\right)=\left(\begin{array}{cc}{s}_{11}& {s}_{12}\\ {s}_{21}& {s}_{22}\end{array}\right)\left(\begin{array}{c}{E}_{x}^{\left({f}_{2}\right)}\\ {H}_{y}^{\left({f}_{2}\right)}\end{array}\right)\\ =\left(\begin{array}{cc}{s}_{11}& {s}_{12}\\ {s}_{21}& {s}_{22}\end{array}\right)\left(\begin{array}{c}{E}_{x}^{\left(3\right)}\\ {H}_{y}^{\left(3\right)}\end{array}\right)\therefore \left(\begin{array}{c}{E}_{x}^{\left({f}_{1}\right)}/{E}_{x}^{\left(3\right)}\\ {H}_{y}^{\left({f}_{1}\right)}/{E}_{x}^{\left(3\right)}\end{array}\right)\\ =\left(\begin{array}{cc}{s}_{11}& {s}_{12}\\ {s}_{21}& {s}_{22}\end{array}\right)\left(\begin{array}{c}1\\ {\eta}_{s}^{\left(3\right)}\end{array}\right)\therefore {\eta}_{s}^{\left({f}_{1}\right)}\\ =\frac{{H}_{y}^{\left({f}_{1}\right)}}{{E}_{x}^{\left({f}_{1}\right)}}=\frac{{s}_{11}+{s}_{12}\left({\eta}_{s}^{\left(3\right)}\right)}{{s}_{21}+{s}_{22}\left({\eta}_{s}^{\left(3\right)}\right)}\end{array}& 8.\end{array}$  In these equations, the superscript (3) indicates the exit medium. The admittance of the exit medium for spolarized light is given by η_{s} ^{(3)}=√{square root over (∈_{3}−∈_{inc }sin^{2}(θ_{inc}))}. The negative sign in front of η_{s} ^{(3) }in the second line results from defining light as propagating in the negative z direction. For spolarized light, the magnetic and electric fields are of opposite signs in this case. For ppolarized light, they are of the same sign. The 2×2 matrix for the single film is described below.
 For ppolarized light, the admittance of the initial interface is arrived at in the same way:
$\begin{array}{cc}\begin{array}{c}\left(\begin{array}{c}{E}_{y}^{\left({f}_{1}\right)}\\ {H}_{x}^{\left({f}_{1}\right)}\end{array}\right)=\left(\begin{array}{cc}{p}_{11}& {p}_{12}\\ {p}_{21}& {p}_{22}\end{array}\right)\left(\begin{array}{c}{E}_{y}^{\left({f}_{2}\right)}\\ {H}_{x}^{\left({f}_{2}\right)}\end{array}\right)\\ =\left(\begin{array}{cc}{p}_{11}& {p}_{12}\\ {p}_{21}& {p}_{22}\end{array}\right)\left(\begin{array}{c}{E}_{y}^{\left(3\right)}\\ {H}_{x}^{\left(3\right)}\end{array}\right)\therefore \left(\begin{array}{c}{E}_{y}^{\left({f}_{1}\right)}/{E}_{y}^{\left(3\right)}\\ {H}_{x}^{\left({f}_{1}\right)}/{E}_{y}^{\left(3\right)}\end{array}\right)\\ =\left(\begin{array}{cc}{p}_{11}& {p}_{12}\\ {p}_{21}& {p}_{22}\end{array}\right)\left(\begin{array}{c}1\\ {\eta}_{p}^{\left(3\right)}\end{array}\right)\therefore {\eta}_{p}^{\left({f}_{1}\right)}\\ =\frac{{H}_{x}^{\left({f}_{1}\right)}}{{E}_{y}^{\left({f}_{1}\right)}}=\frac{{p}_{11}+{p}_{12}{\eta}_{p}^{\left(3\right)}}{{p}_{21}+{p}_{22}{\eta}_{p}^{\left(3\right)}}\end{array}& 9.\end{array}$  For ppolarized light, the admittance of the exit medium is written as η_{p} ^{(3)}=∈_{3}/√{square root over (∈_{3}−∈_{inc }sin^{2}(θ_{inc}))}. The 2×2 matrices for s and p polarizations are defined by:
$\begin{array}{cc}{s}_{12}={s}_{22}={p}_{11}={p}_{22}=\mathrm{cos}\text{\hspace{1em}}\left({\delta}_{\mathrm{film}}\right)\text{}{s}_{12}=\frac{i}{{\eta}_{s}^{\left(\mathrm{film}\right)}}\mathrm{sin}\text{\hspace{1em}}\left({\delta}_{\mathrm{film}}\right)\text{}{s}_{21}=i\text{\hspace{1em}}{\eta}_{s}^{\left(\mathrm{film}\right)}\mathrm{sin}\text{\hspace{1em}}\left({\delta}_{\mathrm{film}}\right)\text{}{p}_{12}=\frac{i}{{\eta}_{p}^{\left(\mathrm{film}\right)}}\mathrm{sin}\text{\hspace{1em}}\left({\delta}_{\mathrm{film}}\right)\text{}{p}_{21}=i\text{\hspace{1em}}{\eta}_{p}^{\left(\mathrm{film}\right)}\mathrm{sin}\text{\hspace{1em}}\left({\delta}_{\mathrm{film}}\right)& 10.\end{array}$  In these equations, the admittances of the film are written in the same form as the admittances of the exit medium given above, but with ∈film replacing ∈3. The value δfilm is the phase thickness of the film, given by
$\begin{array}{cc}{\delta}_{\mathrm{film}}=\frac{2\text{\hspace{1em}}\pi \text{\hspace{1em}}{d}_{\mathrm{film}}\sqrt{{\varepsilon}_{\mathrm{film}}{\varepsilon}_{\mathrm{inc}}{\mathrm{sin}}^{2}\left({\theta}_{\mathrm{inc}}\right)}}{{\lambda}_{0}}& 11.\end{array}$
where dfilm is the physical thickness of the film and λ0 is the freespace wavelength of the incident light.  If there are multiple films, the matrix for a stack of films is obtained from
$\begin{array}{cc}S=\left(\begin{array}{cc}{s}_{11}& {s}_{12}\\ {s}_{21}& {s}_{22}\end{array}\right)=\prod _{\mathrm{inc}}^{\mathrm{exit}}{S}_{\mathrm{films}}& 12.\end{array}$
where the product is over the 2×2 matrices of each individual film from the entrance to the exit. The final product matrix is used as though it described a single equivalent layer.  Notably in the preceding calculation, the matrices describing the film are used as “transfer” matrices. This permits propagation of the calculation of the admittance of the initial interface down through a stack of films. The downward propagation is stopped at the substrate because, once there are no longer rays propagating in both directions, a simple form (the admittance of the exit medium) can be written. Thus, the matrices allow an impossible calculation to be related to a simplified calculation via a 2×2 matrix.
 A critical piece of understanding results from the following discussion. Referring to Equation 8 above for s polarization, −η_{s} ^{(3)}, is the admittance of the final interface (thus the basis for preserving the negative sign). In comparison, p polarization, η_{p} ^{(3)}, in Equation 6 is the admittance of the final interface. Thus, both the s and p calculation of admittance for the initial interface can be written in a generic form:
$\begin{array}{cc}{\eta}^{\left({f}_{1}\right)}=\frac{{H}^{\left({f}_{1}\right)}}{{E}^{\left({f}_{1}\right)}}=\frac{{m}_{11}+{m}_{12}{\eta}^{\left({f}_{\omega}\right)}}{{m}_{21}+{m}_{22}{\eta}^{\left({f}_{\omega}\right)}}& 13.\end{array}$
where the matrix elements are for either the s or p matrices, and η^{(fω) }is the admittance of the final interface. The final interface is always chosen because a simple expression for its admittance can be written in terms of the admittance of the exit medium.  A Bootstrap Method according to an aspect of the invention depends on finding experimental values for the complex reflectance at a given interface in a film stack. Thus, an initial matter of using the amplitude reflectance of a film surface alone to complete the matrix calculation for the films above the film surface in question will be described.
 Solving the Top of the Stack
 Turning now to a problem illustrated in
FIG. 2 , reflectance at some interface 160 (mathematically, k′) below a film stack 115AX is assumed. Earlier layers 114AX are shown “grayed out” and with diagonal lines to indicate a vague idea of what those layers 114AX are. The reflectance of the “known” layer 114A differs in value from the reflectance described above (thus, a “prime” symbol on k at the interface 116 indicates the different value). However, it is still desirable to be able to compute the reflectance spectrum of the stack of films 115AX layered on top of the “known” layer 114A.  It is possible to compute the spectrum of a film stack when the amplitude reflectance at one interface at the bottom is known. To appreciate how this calculation is done, it is useful to review the known method for calculating reflectance and understand what the assumptions are.
 If the amplitude reflectivity for an interface k (which can be any interface, including the final interface) is known, an equivalent expression to Equation 13 can be obtained in terms of the admittance of that interface in lieu of carrying the calculation all the way down to the substrate. The admittance for the k^{th }interface can be written as follows:
$\begin{array}{cc}\mathrm{for}\text{\hspace{1em}}s\text{:}\left\{\begin{array}{c}\left(\begin{array}{c}{E}_{x}^{\left({f}_{k}\right)}\\ {H}_{y}^{\left({f}_{k}\right)}\end{array}\right)=\left(\begin{array}{c}\left({r}_{k}^{\prime}+1\right){E}_{x}^{\left({i}_{2}\right)}\\ \left({r}_{k}^{\prime}1\right){\eta}_{s}^{\left(2\right)}{\hat{E}}_{x}^{\left({i}_{2}\right)}\end{array}\right)\\ {\eta}_{s}^{\left({f}_{k}\right)}={\eta}_{s}^{\left(2\right)}\frac{\left({r}_{k}^{\prime}1\right)}{\left({r}_{k}^{\prime}+1\right)}\end{array}\right\}\text{}\mathrm{for}\text{\hspace{1em}}p\text{:}\left\{\begin{array}{c}\left(\begin{array}{c}{E}_{y}^{\left({f}_{k}\right)}\\ {H}_{x}^{\left({f}_{k}\right)}\end{array}\right)=\left(\begin{array}{c}\left(1+{r}_{k}^{\prime}\right){E}_{y}^{\left({i}_{2}\right)}\\ \left(1{r}_{k}^{\prime}\right){\eta}_{p}^{\left(2\right)}{E}_{y}^{\left({i}_{2}\right)}\end{array}\right)\\ {\eta}_{p}^{\left({f}_{k}\right)}={\eta}_{p}^{\left(2\right)}\frac{\left(1+{r}_{k}^{\prime}\right)}{\left(1{r}_{k}^{\prime}\right)}\end{array}\right\}& 14.\end{array}$  In this expression for the admittance of the known interface, “2” indicates the admittance of the film deposited directly on the known interface. This can then used as the starting point for computing the reflectance of a stack of films above the known interface. The admittance of the top interface can be written as:
$\begin{array}{cc}{\eta}_{s}^{\left({f}_{1}\right)}=\frac{{H}_{y}^{\left({f}_{1}\right)}}{{\hat{E}}_{x}^{\left({f}_{1}\right)}}=\frac{{s}_{21}^{\prime}+{s}_{22}^{\prime}{\eta}_{s}^{\left({f}_{k}\right)}}{{s}_{11}^{\prime}+{s}_{12}^{\prime}{\eta}_{s}^{\left({f}_{k}\right)}}=\frac{{s}_{21}^{\prime}+{s}_{22}^{\prime}{\eta}_{s}^{\left(2\right)}\left(\frac{{r}_{k}^{\prime}1}{{r}_{k}^{\prime}+1}\right)}{{s}_{12}^{\prime}+{s}_{12}^{\prime}{\eta}_{s}^{\left(2\right)}\left(\frac{{r}_{k}^{\prime}1}{{r}_{k}^{\prime}+1}\right)}\text{}{\eta}_{p}^{\left({f}_{1}\right)}=\frac{{H}_{x}^{\left({f}_{1}\right)}}{{\hat{E}}_{y}^{\left({f}_{1}\right)}}=\frac{{p}_{21}^{\prime}+{p}_{22}^{\prime}{\eta}_{p}^{\left({f}_{k}\right)}}{{p}_{11}^{\prime}+{p}_{12}^{\prime}{\eta}_{p}^{\left({f}_{k}\right)}}=\frac{{p}_{21}^{\prime}+{p}_{22}^{\prime}{\eta}_{p}^{\left(2\right)}\left(\frac{1{r}_{k}^{\prime}}{1+{r}_{k}^{\prime}}\right)}{{p}_{11}^{\prime}+{p}_{12}^{\prime}{\eta}_{p}^{\left(2\right)}\left(\frac{1{r}_{k}^{\prime}}{1+{r}_{k}^{\prime}}\right)}& 15.\end{array}$  In addition to the change in definition for the terminal interface of the calculation, another difference is the matrix elements come from a modified 2×2 matrix for the film stack. The modified matrix is computed as:
$\begin{array}{cc}{S}^{\prime}=\left(\begin{array}{cc}{s}_{11}^{\prime}& {s}_{12}^{\prime}\\ {s}_{21}^{\prime}& {s}_{22}^{\prime}\end{array}\right)=\prod _{\mathrm{inc}}^{k+}{S}_{\mathrm{films}}\text{}{P}^{\prime}=\left(\begin{array}{cc}{p}_{11}^{\prime}& {p}_{12}^{\prime}\\ {p}_{21}^{\prime}& {p}_{22}^{\prime}\end{array}\right)=\prod _{\mathrm{inc}}^{k+}{P}_{\mathrm{films}}& 16.\end{array}$
where the product is taken in the order of incident light penetrating the stack as before, but the calculation ends with the film deposited directly onto the known interface. The symbol k+ in Equation 16 is used to indicate that the product terminates with a layer 115A directly above the known interface 114A. The change is a significant one, in that the 2×2 matrices for any of the films below the known interface 114A no longer have to be computed.  Returning to Equation 6, the following can be expressed:
$\begin{array}{cc}{r}_{s}^{\left({f}_{1}\right)}=\frac{\begin{array}{c}{\eta}_{s}^{\left(\mathrm{inc}\right)}{s}_{11}^{\prime}\left(1+{r}_{k}^{\prime}\right){\eta}_{s}^{\left(\mathrm{inc}\right)}{\eta}_{s}^{\left(2\right)}{s}_{12}^{\prime}\left(1{r}_{k}^{\prime}\right)+\\ {s}_{21}^{\prime}\left(1+{r}_{k}^{\prime}\right)+{\eta}_{s}^{\left(2\right)}{s}_{22}^{\prime}\left(1{r}_{k}^{\prime}\right)\end{array}}{\begin{array}{c}{\eta}_{s}^{\left(\mathrm{inc}\right)}{s}_{11}^{\prime}\left(1+{r}_{k}^{\prime}\right){\eta}_{s}^{\left(\mathrm{inc}\right)}{\eta}_{s}^{\left(2\right)}{s}_{12}^{\prime}\left(1{r}_{k}^{\prime}\right)\\ {s}_{21}^{\prime}\left(1+{r}_{k}^{\prime}\right)+{\eta}_{s}^{\left(2\right)}{s}_{22}^{\prime}\left(1{r}_{k}^{\prime}\right)\end{array}}\text{}{r}_{p}^{\left({f}_{1}\right)}=\frac{\begin{array}{c}{\eta}_{p}^{\left(\mathrm{inc}\right)}{p}_{11}^{\prime}\left(1+{r}_{k}^{\prime}\right)+{\eta}_{p}^{\left(\mathrm{inc}\right)}{\eta}_{p}^{\left(2\right)}{p}_{12}^{\prime}\left(1{r}_{k}^{\prime}\right)\\ {p}_{21}^{\prime}\left(1+{r}_{k}^{\prime}\right){\eta}_{p}^{\left(2\right)}{p}_{22}^{\prime}\left(1{r}_{k}^{\prime}\right)\end{array}}{\begin{array}{c}{\eta}_{p}^{\left(\mathrm{inc}\right)}{p}_{11}^{\prime}\left(1+{r}_{k}^{\prime}\right)+{\eta}_{p}^{\left(\mathrm{inc}\right)}{\eta}_{p}^{\left(2\right)}{p}_{12}^{\prime}\left(1{r}_{k}^{\prime}\right)+\\ {p}_{21}^{\prime}\left(1+{r}_{k}^{\prime}\right)+{\eta}_{p}^{\left(2\right)}{p}_{22}^{\prime}\left(1{r}_{k}^{\prime}\right)\end{array}}& 17.\end{array}$  Note that these equations feature a complex quantity called r′_{k}, which, as mentioned above, is not the same as r_{k}, the amplitude reflectance of the top of the film stack before the topmost layer was added. The two things are related to one another, however, as evident in the following discussion.
 Amplitude Reflectance of a Buried Interface
 As discussed above, the magnitude of the reflectance of the interface 116 (mathematically, k) in air can be learned by measuring its intensity reflectance but not the phase of the reflectance in the complex plane. With reference to
FIG. 3 , when the interface 116 (k) is covered by another material 115, the known reflectance is changed. If how the reflectance changes cannot be computed in a simple way, then having learned anything about that reflectance is of no use. Fortunately, there is a straightforward way to relate this to the amplitude reflectance of the interface in air.  In Optical Properties of Thin Solid Films (Dover Publications, Inc., Mineola, USA, 1991), O. S. Heavens gives an expression for the amplitude reflectance of a film in terms of the reflectance of the two interfaces of the film:
$\begin{array}{cc}{r}_{\mathrm{film}}=\frac{{r}_{\mathrm{top}}+{r}_{\mathrm{bot}}{e}^{\mathrm{i2\delta}}}{1+{r}_{\mathrm{top}}{r}_{\mathrm{bot}}{e}^{\mathrm{i2\delta}}}& 18.\end{array}$
where δ is the optical phase change$\delta =\frac{2\pi \text{\hspace{1em}}d\sqrt{{\varepsilon}_{2}{\varepsilon}_{\mathrm{inc}}{\mathrm{sin}}^{2}\left({\theta}_{\mathrm{inc}}\right)}}{{\lambda}_{0}},$
directly proportional to the physical thickness of the film. If rtop is replaced with r2 (the Fresnel coefficient for reflectance off an infinite slab of the film material with dielectric constant ∈2), and rbot is allowed to be r′_{k}, the reflectance of the multilayer stack when the entrance medium is an infinite slab of film material, then the reflectance of the film's top interface can be written as:$\begin{array}{cc}r\left({f}_{1}\right)=\frac{{r}_{2}+{r}_{k}^{\prime}{e}^{{\mathrm{i2\delta}}_{2}}}{1+{r}_{2}{r}_{k}^{\prime}{e}^{{\mathrm{i2\delta}}_{2}}}& 19.\end{array}$  When the thickness of the film goes to zero, the exponential equals 1 and the film reflectance must be identical to r_{k}. This allows r′_{k }to be solved in terms of r_{k }as:
$\begin{array}{cc}{r}_{k}^{\prime}=\frac{{r}_{k}{r}_{2}}{1{r}_{2}{r}_{k}}& 20.\end{array}$  This provides an estimate of r′_{k }that is partially independent of the preceding film stack, since r_{2 }does not depend on it at all and r_{k }has been modified, keeping only the phase determined by the film stack calculation.
 The Bootstrap Method.
 In light of the foregoing introduction, a Bootstrap method for film deposition and refinement is described in the following sections; more particularly, steps to perform Bootstrap refinement of the optical model of a thin film stack are provided as follows.
 Step 1. Determine the reflectance of an existing film stack prior to the deposition of a new layer.
 The reflectance of a film stack provides some information regarding the complex amplitude reflectance that can be used to refine the model of the reflectance, and that is totally independent of any modeling. If one does not measure reflectance directly, it can be obtained by noting that transmission plus reflectance for an absorptionfree thin film stack is unity.
 The relationship between the amplitude and intensity reflectance is that the intensity reflectance is the absolute square of the amplitude reflectance. Considering the amplitude reflectance for a moment, it will be clear that it can be expressed in standard Cartesian coordinates on a complex plane, or in complex polar coordinates:
$\begin{array}{cc}{r}_{k}=a+i\text{\hspace{1em}}b=\uf603{r}_{k}\uf604{e}^{{\mathrm{i\varphi}}_{k}}\text{}\uf603{r}_{k}\uf604=\sqrt{{a}^{2}+{b}^{2}}\text{}{\varphi}_{k}={\mathrm{tan}}^{1}\left(\frac{b}{a}\right).& 1\end{array}$  If the amplitude reflectance is expressed in polar coordinates, it is the magnitude of the amplitude reflectance that is provided by a measure of intensity reflectance, r_{k}=√{square root over (R_{k})}.
 It is possible at this point to replace the magnitude of the amplitude reflectivity in Equation 1, r_{k}, with the square root of the intensity reflectance. This can be done whenever the anticipated error in future reflectance values due to errors in this initial measurement of R_{k }is expected to be small. How to obtain this relationship is shown in the following.
 Calculating the Worstcase Future Reflectance
 The standard deviation of the magnitude of the amplitude reflectance is given by Equation 3:
$\begin{array}{cc}\uf603{r}_{k}\uf604=\sqrt{{R}_{k}}\text{}{\varepsilon}_{\uf603{r}_{k}\uf604}=\frac{1}{\sqrt{{R}_{k}}}{\varepsilon}_{{R}_{k}}\text{}{\sigma}_{\uf603{r}_{k}\uf604}=\frac{{\sigma}_{{R}_{k}}}{\sqrt{{R}_{k}}}.& 3\end{array}$
In the worstcase scenario, Rk is a minimum (the amplitude reflection is on the real axis nearest the origin). This would make the error in magnitude relatively larger. Further, the next layer could result in this vector being advanced by π, crossing the real axis at the furthest point from the origin, producing a reflectance maximum. Again, in the worst case scenario, the maximum reflectance that could be generated as a result of the observed Rk is:$\begin{array}{cc}{R}_{\mathrm{max}}\left(\mathrm{max}\right)={\left(\frac{\left(1+{r}_{2}^{2}\right){R}_{k}^{1/2}2{r}_{2}}{1+{r}_{2}^{2}2{r}_{2}{R}_{k}^{1/2}}\right)}^{2}.& 4\end{array}$  To assure that replacing the reflectance amplitude with a measured value does not affect future reflectance measurements by more than the standard deviation of the reflectance measurement, the worstcase scenario must be known. This is:
$\begin{array}{cc}{\sigma}_{{R}_{\mathrm{max}}}\left(\mathrm{max}\right)=\frac{{\left({r}_{2}^{2}1\right)}^{2}\left(1+{r}_{2}^{2}2{r}_{2}{R}_{k}^{1/2}\right)}{{\left(1+{r}_{2}^{2}2{r}_{2}{R}_{k}^{1/2}\right)}^{3}}{\sigma}_{{R}_{k}}.& 5\end{array}$
Solving Equation 5 for a factor of 1σ is not easy, and the result is very complicated. However, a numerical solution is straightforward. For silica (r2=−0.2), this value is about Rk=0.19 or 19% reflectance. For niobia (r2=−0.4), the value is about 9% reflectance. In other words, when about to deposit a silica layer, values of r_{k} should not be adjusted when the measured reflectance is less than 19%. When about to deposit a niobia layer, values should not be replaced when the measured reflectance is less than 9%. Instead, assume the modeled reflectance is more accurate in these cases, although it is not necessary to do this often. In the following sections, reasons are discussed to perform bootstrapping only on highindex layers, so by extension, this step is recommended only when a lowindex layer is completed. 
 Step 2: Replace the magnitude of the amplitude reflectance at each wavelength with √{square root over (R_{k})} whenever measuring a freshly completed silica film with an intensity reflectance greater than 9% or a lowindex film with an intensity reflectance greater than the limiting value of the highindex material.
 While the intensity measurement provides useful information (most of the time) about the magnitude of reflectance, it unfortunately provides no information about the phase angle in the complex plane, φ. Much of the remainder of the present description relates to how to obtain these phase angles in at least some circumstances.
 Monitor Curves Can Give NonRedundant Calculations of Phase
 In most cases, the magnitude of r_{k }at the base of a niobia layer can be obtained from the measured reflectance of the film stack terminating in a fresh silica layer. The phase of the amplitude reflectance is more difficult to ascertain, but there are two general approaches. The first is to consider what values of phase are consistent with the final value of reflectance after the next layer is added. To use this information, the optical thickness of the next layer must be known. This is sometimes a redundant calculation since the estimation of optical thickness is usually based on an understanding of the initial reflectance. This is, in fact, a weakness of the usual matrix modeling approach—the calculation is somewhat redundant.
 Without additional information, redundant calculation would normally be the only option. However, monitor curves are usually recorded during deposition, and those curves contain all the information necessary to compute the phase, φ_{k}, without the need for redundant calculations. This involves the use of reflectance maxima in the monitor curves.

 Step 3. Consider the modeled monitor curve for each wavelength of a niobia (highindex) layer. Discard any monitor curves without maxima in their reflectance during the niobia layer deposition. If none meet this criterion, deposit the layer using a pure model approach.
 Based on Equation 14 above, the maxima and minima of a monitor curve can be shown to depend solely on r′_{k}, the magnitude of the buried interface's reflectance, and not at all on its phase. The maximum and minimum reflectance during the monitor curve are given by Equation 21:
$\begin{array}{cc}{R}_{\mathrm{max}}={\left(\frac{\uf603{r}_{k}^{\prime}\uf604+\uf603{r}_{2}\uf604}{1+\uf603{r}_{2}\uf604\uf603{r}_{k}^{\prime}\uf604}\right)}^{2}\text{}{R}_{\mathrm{min}}={\left(\frac{\uf603{r}_{k}^{\prime}\uf604\uf603{r}_{2}\uf604}{1\uf603{r}_{2}\uf604\uf603{r}_{k}^{\prime}\uf604}\right)}^{2}.& 21\end{array}$  Therefore, R_{max }or R_{min }can be used in the monitor curve to convey the magnitude of the buried reflectance. This magnitude can be related to the reflectance as follows.
$\begin{array}{cc}\uf603{r}_{k}^{\prime}\uf604=\frac{\sqrt{{R}_{\mathrm{max}}}\uf603{r}_{2}\uf604}{1\uf603{r}_{2}\uf604\sqrt{{R}_{\mathrm{max}}}}\text{}\mathrm{or}\text{}\uf603{r}_{k}^{\prime}\uf604=\frac{\sqrt{{R}_{\mathrm{min}}}+\uf603{r}_{2}\uf604}{1+\uf603{r}_{2}\uf604\sqrt{{R}_{\mathrm{min}}}};\frac{\uf603{r}_{2}\uf604\sqrt{{R}_{\mathrm{min}}}}{1\uf603{r}_{2}\uf604\sqrt{{R}_{\mathrm{min}}}}.& 22\end{array}$
Thus, from a monitor curve covering at least a quarter wave at the monitor curve wavelength, r′_{k} can be determined.  A caveat to using these equations is as follows. The R_{min }expression has two failings. First, there are two possible solutions for r′_{k} based on R_{min }depending on whether r′_{k} is less than or greater than r_{2}. If it is less than r_{2}, the righthand solution is appropriate. If it is greater than r_{2}, the lefthand solution is appropriate. The R_{max }expression also has two solutions in principle but can be discarded because it provides nonphysical results. The second problem with the equation from R_{min }is the issue of experimental error. The error expected in the estimation of r′_{k} is related to the error in measurement of R_{min }and R_{max }by Equation 23:
$\begin{array}{cc}{\varepsilon}_{\uf603{r}_{k}^{\prime}\uf604}=\frac{\left(1\uf603{r}_{2}^{2}\uf604\right)}{2\sqrt{{R}_{\mathrm{min}}}{\left(1+\uf603{r}_{2}\uf604\sqrt{{R}_{\mathrm{min}}}\right)}^{2}}{\varepsilon}_{{R}_{\mathrm{min}}}\approx \frac{{\varepsilon}_{{R}_{\mathrm{min}}}}{2\sqrt{{R}_{\mathrm{min}}}}\text{}{\varepsilon}_{\uf603{r}_{k}^{\prime}\uf604}=\frac{\left(1\uf603{r}_{2}^{2}\uf604\right)}{2\sqrt{{R}_{\mathrm{max}}}{\left(1\uf603{r}_{2}\uf604\sqrt{{R}_{\mathrm{max}}}\right)}^{2}}{\varepsilon}_{{R}_{\mathrm{max}}}\approx \frac{{\varepsilon}_{{R}_{\mathrm{max}}}}{2\sqrt{{R}_{\mathrm{max}}}}.& 23\end{array}$
In other words, the error goes up as the key reflectance diminishes. Since the minimum is, by definition, smaller than the maximum, the error expected in estimating r′_{k} goes up accordingly. Thus, for both reasons, the calculation of magnitude from a maximum reflectance (i.e., a minimum in the transmission monitor curve) is preferred.
Choosing the Best Monitor Wavelength  The phase angle φk at the monitor wavelength can be determined from r′_{k} and the Rk as:
$\begin{array}{cc}\mathrm{cos}\left(\pm {\varphi}_{k}\right)=\frac{{\uf603{r}_{k}^{\prime}\uf604}^{2}\left(1+{R}_{k}{\uf603{r}_{2}\uf604}^{2}\right){\uf603{r}_{2}\uf604}^{2}{R}_{k}}{2\uf603{r}_{2}\uf604\sqrt{{R}_{k}}\left(1{\uf603{r}_{k}^{\prime}\uf604}^{2}\right)}.& 25\end{array}$
This, of course, provides 2 solutions. Once the phase angle φk from Equation 25 is determined, r′_{k }can be computed using Equation 20 and a monitor curve can be generated if desired for comparison with the actual to help determine which solution is better. Afterwards, a value for phase at the monitor wavelength should have been obtained that is as correct as possible. It depends, of course, on accurately measuring the maximum reflectance value and Rk. Thus, not all wavelengths are created equal as potential monitor wavelengths. A fullspectrum monitor (acquiring many monitor wavelengths) is the best solution, but if only a single wavelength is available, then there is a systematic approach to choosing the best.  Equation 25 depends, ultimately, on only two measurements: the measurement of the initial reflectance and the measurement of the maximum reflectance. For those wavelengths that exhibit a maximum reflectance during deposition, these can be evaluated quantitatively as possible monitor wavelengths.
 It can be shown that the anticipated standard deviation of the phase calculation can be written as:
$\begin{array}{cc}{\sigma}_{{\varphi}_{k}}=\sqrt{A\text{\hspace{1em}}{\sigma}_{{R}_{k}}^{2}+B\text{\hspace{1em}}{\sigma}_{{R}_{\mathrm{max}}}^{2}}\text{}A=\frac{{\left(\left(1+{r}_{2}^{2}\right){R}_{k}+2{r}_{2}\left(1+{R}_{k}\right)\sqrt{{R}_{\mathrm{max}}}+\left(1+{r}_{2}^{2}\right){R}_{\mathrm{max}}\right)}^{2}}{4{R}_{k}^{2}\left(\begin{array}{c}4{r}_{2}^{2}{{R}_{k}\left(1{R}_{\mathrm{max}}\right)}^{2}\\ {\left(\begin{array}{c}2{r}_{2}\sqrt{{R}_{\mathrm{max}}}+\left(1+{r}_{2}^{2}\right){R}_{\mathrm{max}}\\ {R}_{k}\left(1+{r}_{2}^{2}+2{r}_{2}\sqrt{{R}_{\mathrm{max}}}\right)\end{array}\right)}^{2}\end{array}\right)}\text{}B=\frac{{\left(1{R}_{k}\right)}^{2}{\left({r}_{2}+\left(1+{r}_{2}^{2}\right)\sqrt{{R}_{\mathrm{max}}}+{r}_{2}{R}_{\mathrm{max}}\right)}^{2}}{{{R}_{\mathrm{max}}\left(1{R}_{\mathrm{max}}\right)}^{2}\left(\begin{array}{c}4{r}_{2}^{2}{{R}_{k}\left(1{R}_{\mathrm{max}}\right)}^{2}\\ {\left(\begin{array}{c}2{r}_{2}\sqrt{{R}_{\mathrm{max}}}+\left(1+{r}_{2}^{2}\right){R}_{\mathrm{max}}\\ {R}_{k}\left(1+{r}_{2}^{2}+2{r}_{2}\sqrt{{R}_{\mathrm{max}}}\right)\end{array}\right)}^{2}\end{array}\right)}.& 27\end{array}$  This equation is helpful selecting the best monitor wavelength for the purpose of determining the phase of the amplitude reflectivity at the monitor wavelength.

FIG. 4 shows a representation of the regions of permissible values of Rk and Rmax, with a color code for the lowest values of σ_{φ}k for silica (assuming r2=−0.2 and the standard deviation of the reflectance measurements is 0.003). The lowest value possible under these conditions is 0.0262 radians (1.5 degrees). The lower axis, Rmax, represents possible values of Rmax, while the left axis, Rk, gives possible values of Rk. Note that large regions of reflectance are not possible—there are many combinations of Rk and Rmax that cannot coexist. On the boundaries of those disallowed regions, the error in estimating the phase angle, φk, becomes infinite.  The same plot for niobia, assuming r_{2}=−0.4, is given in
FIG. 5 . The minimum value of error here under the same conditions is 0.00875 radians (0.5 degrees)—a much better phase calculation.  It is possible to develop a histogram of the number of combinations of allowed Rmax and Rk that provide a specific level of error in φk. This is accomplished first by considering the range of possible Rmax values when depositing a layer: it cannot, as
FIGS. 4 and 5 illustrate, be less than the reflectance of the thin film material being deposited itself (e.g., Rmax for niobia inFIG. 5 cannot be less than −0.4^{2}=0.2 as shown in a leftmost portion ofFIG. 5 in dark yellow). Possible values of Rk can be evaluated for each value of Rmax, and the values of phase precision those values provide can be determined using Equation 27. This is accomplished by dividing the phase thickness of the top layer into increments and computing the reflectance at each increment (this is done because the reflectances at the turning points are more likely than those in between the turning point values). In the end, a histogram of the resulting precision values can be formed and a determination made as to how likely each will appear for a given film material. This has been done for silica (assuming r_{2}=−0.2, and reflectance standard deviations of 0.003) and for niobia (assuming r_{2}=−0.4) as shown respectively inFIGS. 6 and 7 .  For silica, 30% of all observed combinations will have a phase error given by Equation 24 that is less than 2.4 degrees. No values less than about 1.5 degrees error in φk is possible for silica under these conditions. For niobia, the same fraction will have φk errors less than 0.9 degrees, as illustrated in the
FIG. 7 , also derived from a full numeric simulation.  Thus, the calculated phase error of possible monitor wavelengths will tend to be considerably better for niobia films than for silica. If a limit of 0.9 degrees phase error is placed on the monitors before this calculation is performed, only niobia will give possible monitor wavelengths, and 30% of all wavelengths (overall) will meet this criterion. On some layers, it is possible that no wavelengths will meet this criterion, while on others many may do so.

 Step 4. For the remaining possible monitor wavelengths in a niobia layer deposition, determine the anticipated standard deviation in φk. Discard any with cy greater than 0.9 degrees (0.016 radians). If none remain, proceed with a pure model deposition.
 For lowindex layers, bootstrapping is not recommended. For layers with larger magnitudes of r_{2 }(such as niobia), the precision of the bootstrap is almost always better, but there is no guarantee that a specific layer will include a set of R_{k }and R_{max }values anticipated to provide excellent precision in calculating the phase φk. If no wavelength with an anticipated reflectance maximum meets this criterion, modeling alone should be relied upon to deposit the layer until a suitable bootstrap layer is reached.
 An ingenious characteristic about this calculation is that it provides φk—and thus also φ′_{k}—that is consistent with the monitor curve and is independent of δ, the phase thickness of the film. Once a valid solution for rk is determined, the valid solution for r′_{k }can be obtained. Thus, the final transmission value at the monitor curve wavelength can be used to determine what value of δ is most accurate for the monitor wavelength.
 If a wavelength meets the criterion specified above, then consideration should be given as to whether the anticipated end of the layer will have a reflectance, Rf, suitable for estimating δ, the phase thickness of the layer with some precision. This calculation is performed with Equation 26, where δ depends on Rf and rk′.
$\begin{array}{cc}\delta =\frac{{\varphi}_{k}^{\prime}}{2}\pm \frac{1}{2}{\mathrm{cos}}^{1}\left(\frac{{R}_{f}+{R}_{f}{\uf603{r}_{2}\uf604}^{2}{\uf603{r}_{k}^{\prime}\uf604}^{2}{\uf603{r}_{2}\uf604}^{2}{\uf603{r}_{k}^{\prime}\uf604}^{2}}{2\uf603{r}_{2}\uf604\uf603{r}_{k}^{\prime}\uf604\left(1{R}_{f}\right)}\right).& 26\end{array}$  This equation provides fairly unique solutions for δ that can be used to correct the reflectance at all wavelengths. The possible solutions for δ that are obtained can be tested against the observed monitor curve to determine which is correct.
 For a given monitor curve with a given expectation of Rmax, and with a known value of r2, the sensitivity of δ to errors in Rf can be determined according to Equation 28:
$\begin{array}{cc}{\sigma}_{\delta}=\uf603\frac{\left(1\uf603{r}_{2}^{2}\uf604\right)\left(1{R}_{\mathrm{max}}\right)}{2\left({R}_{f}1\right)\sqrt{A}}\uf604{\sigma}_{{R}_{f}}\text{}A=\left({R}_{\mathrm{max}}{R}_{f}\right)\left(\begin{array}{c}2\uf603{r}_{2}^{2}\uf604\left({R}_{\mathrm{max}}{R}_{f}+{R}_{f}{R}_{\mathrm{max}}2\right)\\ \left(1+\uf603{r}_{2}^{4}\uf604\right)\left({R}_{\mathrm{max}}{R}_{f}\right)+\\ 4\left(\uf603{r}_{2}\uf604+\uf603{r}_{2}^{3}\uf604\right)\left(1{R}_{f}\right)\sqrt{{R}_{\mathrm{max}}}\end{array}\right).& 28\end{array}$  If the type of numerical analysis above using Equation 28 is repeated, a plot as shown in
FIG. 8 is obtained (shown as the logarithm of error vs. Rmax and δ because otherwise the scale would be difficult to see). This plot is made over the entire range of possible values of Rmax between r_{2} ^{2 }and 1 (on the receding axis) and δ angle between 0 and 2π (the front axis).  The data in
FIG. 8 can also be rendered as a histogram as shown inFIG. 9 . The histogram inFIG. 9 implies that the error in phase thickness is usually satisfactory compared to the error in φk. At a cutoff of 0.9 degrees error in δ, about 51% of the remaining monitor wavelengths for niobia, for instance, should be usable. Thus, given a good monitor wavelength for determining the phase angle φk, there is a good chance of having one that also provides a good precision in δ.  To conserve time, and when modeling is running fairly well, it is reasonable to only select monitor wavelengths for niobia (in a silica/niobia stack), and only when they meet these two criteria (standard deviation of φk<0.9 degrees, and standard deviation of δ<0.9 degrees at the end of the layer). When no monitor wavelengths meet these criteria, it is reasonable to proceed with a pure model matrix approach to depositing the next layer.

 Step 5. Compute the expected error in δ for the remaining wavelengths at the target thickness of the niobia layer. If no wavelengths have an error less than 0.9 degrees, proceed with a full model deposition of the layer. If some do meet this criterion, select the lowest error in this category.
Reflectance of a New Film and Determination of the Old Phase at Wavelengths Other than the Monitor Wavelength.
 Step 5. Compute the expected error in δ for the remaining wavelengths at the target thickness of the niobia layer. If no wavelengths have an error less than 0.9 degrees, proceed with a full model deposition of the layer. If some do meet this criterion, select the lowest error in this category.
 Why bother with determining the φk and δ from a monitor curve as precisely as possible? First, if a monitor wavelength for bootstrapping has been selected successfully, all connection to the previous dependence on the matrix calculation for the monitor wavelength may be avoided. For all other wavelengths, there is at least the opportunity to determine the magnitude of rk, but only the original modeled estimate of phase. The question arises: how to “repair” the phases of all other wavelengths? Since in a singlechannel monitor there are no monitor curves at those wavelengths, the phase at each wavelength cannot be directly obtained. (If there was spectrograph recording all wavelengths all of the time, such as with an FTIR system, the steps provided below would not be needed). However, to reach this point in bootstrapping, there must be a good value for δ for the monitor wavelength, plus good values of Rmax and Rk. With this, at least some of the other measurements can be “fixed” to accord with measurements already taken. This will restrict errors to those of a single layer at the worst for those wavelengths where the following calculation is possible.
 From δ for the monitor wavelength, the physical thickness of the layer consistent with the modeled refractive index of the film material can be estimated. This physical thickness and the modeled refractive index of the film can be used to estimate δ for all other wavelengths.

 Step 6. Compute the value of δ for all wavelengths based on the value calculated for the monitor wavelength.
 Returning to Equation 16 and replacing rk′ with the definition in Equation 20, the following is obtained:
$\begin{array}{cc}{r}_{f}=\frac{{r}_{2}\left(1{r}_{2}{r}_{k}\right)+\left({r}_{k}{r}_{2}\right)\left(\mathrm{cos}\left(2\delta \right)i\text{\hspace{1em}}\mathrm{sin}\left(2\delta \right)\right)}{1{r}_{2}{r}_{k}+{r}_{k}\left({r}_{k}{r}_{2}\right)\left(\mathrm{cos}\left(2\delta \right)i\text{\hspace{1em}}\mathrm{sin}\left(2\delta \right)\right)}.& 29\end{array}$  In Equation 29, the exponential has been replaced with a trigonometric expression using Euler's relation. If r_{k }is replaced with a+i b, a conventional expression for r_{f }can be obtained. The complex conjugate of r_{f }can be formed and the product taken of the two. This provides the intensity reflectance of the top interface in terms of values from the new film, plus a and b. The following can then replace a and b:
a=cos (φ_{k})√{square root over (R _{k})}
b=sin (φ_{k})√{square root over (R _{k})} (30.)
where use is made of the intensity reflectance measured in vacuum for interface k to represent the magnitude of the amplitude reflectance vector at the interface k in vacuum.  The resulting expression can be simplified as Equation 31:
$\begin{array}{cc}\text{\hspace{1em}}{R}_{f}=\frac{2{r}_{2}^{2}+{R}_{k}\left(1+{r}_{2}^{4}\right)+2{r}_{2}Q}{1+{r}_{2}^{4}+2{r}_{2}^{2}{R}_{k}+2{r}_{2}Q}\text{}Q={r}_{2}^{2}{R}_{k}^{1/2}\mathrm{cos}\left(2\delta +{\varphi}_{k}\right)+{R}_{k}^{1/2}\mathrm{cos}\left(2\delta {\varphi}_{k}\right){r}_{2}\left(1+{R}_{k}\right)\mathrm{cos}\left(2\delta \right)\left(1+{r}_{2}^{2}\right){R}_{k}^{1/2}\mathrm{cos}\left({\varphi}_{k}\right).& 31\end{array}$  In this expression, everything is known EXCEPT the phase angle φ at interface k, allowing it to be determined from the measurement of intensity reflectance at the subsequent interface.
 Dispensing with numerical solution methods, this expression can be solved for the Cosine of the angle:
$\begin{array}{cc}\text{\hspace{1em}}\mathrm{cos}\left({\varphi}_{k}\right)=\left(\frac{A\left(1+{r}_{2}^{2}\right)\mathrm{sin}\left(\delta \right)\pm B\text{\hspace{1em}}\mathrm{cos}\left(\delta \right)}{C}\right)\text{}A={R}_{f}+{r}_{2}^{4}\left({R}_{f}{R}_{k}\right){R}_{k}+2{r}_{2}^{2}\left(\begin{array}{c}\left(1{R}_{f}\right)\left(1+{R}_{k}\right)\mathrm{cos}\left(2\delta \right)\\ \left(1{R}_{f}{R}_{k}\right)\end{array}\right)\text{}\text{\hspace{1em}}B=\sqrt{D\left(1+{r}_{2}^{12}\right)+F\left({r}_{2}^{2}+{r}_{2}^{10}\right)+G\left({r}_{2}^{4}+{r}_{8}^{2}\right)+H\text{\hspace{1em}}{r}_{6}^{2}}\text{}\text{\hspace{1em}}C=\mathrm{sin}\left(\delta \right)\left(4{r}_{2}\left(1{R}_{f}\right){R}_{k}^{1/2}\left(2{r}_{2}^{2}\mathrm{cos}\left(2\delta \right)1{r}_{2}^{4}\right)\right)\text{}\text{\hspace{1em}}D={\left({R}_{f}{R}_{k}\right)}^{2}\text{}\text{\hspace{1em}}F=2\left(\begin{array}{c}{R}_{k}\left(2+{R}_{k}\right)+{R}_{f}^{2}\left(1+2{R}_{k}\right)+2{R}_{f}\left(15{R}_{k}+{R}_{k}^{2}\right)\\ 2\left(1{R}_{f}\right)\left(1{R}_{k}\right)\left({R}_{f}+{R}_{k}\right)\mathrm{cos}\left(2\delta \right)\end{array}\right)\text{}G=64{R}_{f}5{R}_{f}^{2}4{R}_{k}+38{R}_{f}{R}_{k}4{R}_{f}^{2}{R}_{k}5{R}_{k}^{2}4{R}_{f}{R}_{k}^{2}6{R}_{f}^{2}{R}_{k}^{2}+8\left(1{R}_{f}^{2}\right)\left(1{R}_{k}^{2}\right)\mathrm{cos}\left(2\delta \right)2{\left(1{R}_{f}\right)}^{2}{\left(1{R}_{k}\right)}^{2}\mathrm{cos}\left(4\delta \right)\text{}\text{\hspace{1em}}H=4\left(\begin{array}{c}3+2{R}_{f}^{2}10{R}_{f}{R}_{k}+2{R}_{k}^{2}+3{R}_{f}^{2}{R}_{k}^{2}\\ 2\left(1{R}_{f}\right)\left(1{R}_{k}\right)\left(2+{R}_{f}+{R}_{k}+2{R}_{f}{R}_{k}\right)\mathrm{cos}\left(2\delta \right)+\\ {\left(1{R}_{f}\right)}^{2}{\left(1{R}_{k}\right)}^{2}\mathrm{cos}\left(4\delta \right)\end{array}\right)& 32\end{array}$  In principle, Equation 32 can be used to solve for the phase angle. By doing so, the deposition system process control is effectively “bootstrapped” by detaching the system completely from everything that came before the last layer. Equation 32 provides four (4) solutions for the phase angle; two come from the +/− portion of the calculation; two more from the fact that cosine is an even function, so positive and negative angles both work equally well. However, only two of these solutions are consistent with the measured value of R_{f}. Thus, solutions should be checked via Equation 31 for validity. Only two solutions should be left after this process is complete.

 Step 7. Compute the two possible values of phase angle for each wavelength other than the monitor wavelength.
Estimating Error for NonMonitor Wavelengths
 Step 7. Compute the two possible values of phase angle for each wavelength other than the monitor wavelength.
 The phase angle is dependent on three reflectivity measurements (R_{max}, R_{f }and R_{k}), those being scrambled together in Equation 32. While this works well for a hypothetical system with no noise, a real spectrometer exhibits errors in measurement of the intensity transmittance.
 Based on work already done, the analysis of Equation 32 is fairly straightforward for extending phase information to other wavelengths. The following expression can be constructed from it:
$\begin{array}{cc}{\varepsilon}_{{\varphi}_{k}}=\frac{1}{\mathrm{sin}\left({\varphi}_{k}\right)}\left(\left(\frac{\partial \mathrm{cos}\left({\varphi}_{k}\right)}{\partial {R}_{k}}\right){\varepsilon}_{{R}_{k}}+\left(\frac{\partial \mathrm{cos}\left({\varphi}_{k}\right)}{\partial {R}_{f}}\right){\varepsilon}_{{R}_{f}}+\left(\frac{\partial \mathrm{cos}\left({\varphi}_{k}\right)}{\partial \delta}\right){\varepsilon}_{\delta}\right)\text{}\text{\hspace{1em}}{\sigma}_{{\varphi}_{k}}\approx \frac{{\sigma}_{R}}{\uf603\mathrm{sin}\left({\varphi}_{k}\right)\uf604}{\left({\left(\frac{\partial \mathrm{cos}\left({\varphi}_{k}\right)}{\partial {R}_{k}}\right)}^{2}+{\left(\frac{\partial \mathrm{cos}\left({\varphi}_{k}\right)}{\partial {R}_{f}}\right)}^{2}\right)}^{1/2}.& 33\end{array}$  In Equation 33 error in phase thickness has been omitted from the calculation since it is of minor concern at this stage. The phase thickness was settled previously for the sake of argument; therefore, making that approximation, computing the sine of the angle and the sensitivity of the cosine functions to errors in R_{k }and R_{f }remain. There is a subtle issue to be considered at this point. What angle should be focused on—the modeled angle or the computed angle using Equations 31/32? In principle, both values are available.
 Thus, the conservative answer to the foregoing question is “both”. Assume, for example, that the model has a complex reflectance rk at a given wavelength from which the phase angle φk and the predicted value of Rk can be obtained. Using a best estimate of the phase thickness at the wavelength, Rf can be computed using Equation 31. Next, compute the phase angle using Equation 32 after varying Rk and Rf each by a small amount—e.g., by 10^{−5}, 1/100^{th }of a percent transmission. Since the modeled phase is known exactly, it will be trivial to identify which results are the ones closest to the modeled phase; therefore, the sensitivities are computed as:
$\begin{array}{cc}\frac{\partial \mathrm{cos}\left({\varphi}_{k}\right)}{\partial {R}_{k}}\approx \frac{\mathrm{cos}\left({\varphi}_{k}\left({R}_{k,m},{R}_{f,m},\delta \right)\right)\mathrm{cos}\left({\varphi}_{k}\left({R}_{k,m}{10}^{5},{R}_{f,m},\delta \right)\right)}{{10}^{5}}& 34\\ \frac{\partial \mathrm{cos}\left({\varphi}_{k}\right)}{\partial {R}_{f}}\approx \frac{\mathrm{cos}\left({\varphi}_{k}\left({R}_{k,m},{R}_{f,m},\delta \right)\right)\mathrm{cos}\left({\varphi}_{k}\left({R}_{k,m},{R}_{f,m}{10}^{5},\delta \right)\right)}{{10}^{5}}.& \text{\hspace{1em}}\end{array}$  From the model, the sine of the phase angle that appears in Equation 33 is trivial to obtain.
 Thus, the standard deviation of the phase calculation can be estimated from the model. If the model were absolutely trustworthy, this would be sufficient; however, this is not the case.

 Step 8. Using information extracted from the model for r_{k }at each wavelength and the computed best value of δ, compute the estimated standard deviation of phase at all wavelengths except the monitor.
 Since the model is not completely trustworthy, the process is repeated using values taken from the calculated value of phase closest to the model value of phase. In other words, for the two possible values of φk, the one closest numerically to the model phase is selected. Now, keying on that value, compute the sensitivities according to Equation 34 using the measured values of R_{k}, R_{f }and the estimate of δ for the wavelength being tested.

 Step 9. Using the computed phase closest to the model phase for r_{k }at each wavelength, the measured R_{f }and R_{k }values and the computed best value of δ, compute the estimated standard deviation of phase at all wavelengths for which the magnitude of r_{k }was estimated (Step 2) other than the monitor.
 Now, compute the estimated error in phase as follows:
σ_{φk}(estimated)=(σ_{φk} ^{2}(model)+σ_{φk} ^{2}(calculated))^{1/2} 35.  What are the probable limits to this calculation? A bit of experimentation suggests that the value of standard deviation in phase angle could be between about that of the monitor measurement (on the low end) to nearly infinite. A sufficient approach is to replace the modeled values with calculated values only if a rather conservative cutoff is obtained. If the error estimates are less than about 1.3 degrees from Equation 35 (square root of two times 0.9), and if the reflectance was greater than 9 percent at the start of the layer (Step 2), then the model and calculated values are averaged together to obtain a new estimate. If either the phase error OR the reflectance criteria are not met, the estimate is returned alone to the model.

 Step 10. If the phase error estimate is less than 1.3 degrees AND the criterion of Step 2 is met, average the calculated and modeled reflectance and phase values to obtain a new value that will be used in all future modeling at that wavelength.
 A further refinement of this approach is to fit the phase values obtained above to a KramersKronig model of the phase to fill in values of the phase that have not been determined previously.
 While preferred embodiments of the invention have been shown and described, those of ordinary skill in the art will recognize that changes and modifications may be made to the foregoing examples without departing from the scope and spirit of the invention. Furthermore, those of ordinary skill in the art will appreciate that the foregoing description is by way of example only, and is not intended to limit the invention so further described in such appended claims. It is intended to claim all such changes and modifications as fall within the scope of the appended claims and their equivalents.
Claims (22)
1. A method using experimental measurements to determine reflectance phase and complex reflectance for arbitrary thin film stacks, the method comprising the steps of:
determining reflectance of a stack of a plurality of films before depositing a topmost layer;
considering a modeled monitor curve for a wavelength of a highindex layer; and
discarding a plurality of monitor curves without maxima in their reflectance during the topmost layer deposition.
2. The method as in claim 1 , wherein the topmost layer is a niobia layer.
3. The method as in claim 2 , further comprising the steps of determining an anticipated standard deviation in φk for a plurality of monitor wavelengths in the niobia layer and discarding any with σ greater than 0.9 degrees.
4. The method as in claim 3 , further comprising the step of computing expected error in δ for wavelengths with σ less than 0.9 degrees at a target thickness of the niobia layer.
5. The method as in claim 4 , further comprising the step proceeding with a full model deposition of the niobia layer when no wavelengths have an error less than 0.9 degrees.
6. The method as in claim 1 , further comprising the step of using only the modeled monitor curve during the topmost layer deposition when there is no maxima in the reflectance of each of the plurality of monitor curves.
7. The method as in claim 1 , further comprising the step of computing a value of δ for all wavelengths based on a value calculated for the monitor wavelength.
8. The method as in claim 1 , further comprising the step of computing two possible values of phase angle for each wavelength other than the monitor wavelength.
9. The method as in claim 1 , further comprising the steps of using information extracted from the model for r_{k }at each wavelength and the computed best value of δ, and computing an estimated standard deviation of phase at all wavelengths except the monitor.
10. The method as in claim 9 , further comprising the steps of using the computed phase closest to the model phase for r_{k }at each wavelength, measured R_{f }and R_{k }values and the computed best value of δ, and computing the estimated standard deviation of phase at all wavelengths for which the magnitude of r_{k }was estimated other than the monitor.
11. The method as in claim 9 , further comprising the steps of determining if a phase error estimate is less than about 1.3 degrees and averaging calculated and modeled reflectance and phase values to obtain a new value for use in subsequent modeling at that wavelength.
12. The method as in claim 1 , wherein the topmost layer is a silica film.
13. The method as in claim 12 , further comprising the step of replacing the magnitude of the amplitude reflectance at each wavelength with √{square root over (R_{k})} whenever measuring a latest depositing silica film having an intensity reflectance greater than 9%.
14. The method as in claim 13 , further comprising the step of determining if a phase error estimate is less than about 1.3 degrees when the magnitude of the amplitude reflectance at each wavelength has been replaced with √{square root over (R_{k})} and averaging calculated and modeled reflectance and phase values to obtain a new value for use in subsequent modeling at that wavelength.
15. A method for correcting thin film stack calculations for accurate deposition of complex optical filters, the method comprising the steps of:
$\mathrm{cos}\left(\pm {\varphi}_{k}\right)=\frac{{\uf603{r}_{k}^{\prime}\uf604}^{2}\left(1+{R}_{k}{\uf603{r}_{2}\uf604}^{2}\right){\uf603{r}_{2}\uf604}^{2}{R}_{k}}{2\uf603{r}_{2}\uf604\sqrt{{R}_{k}}\left(1{\uf603{r}_{k}^{\prime}\uf604}^{2}\right)};$ ${r}_{k}^{\prime}=\frac{{r}_{k}{r}_{2}}{1{r}_{2}{r}_{k}};$
determining phase angle φk at a monitor wavelength from r′_{k} and Rk using a first equation expressed as:
and
estimating r′_{k }using a second equation
and
obtaining a value for phase at a monitor wavelength.
16. A method for automated deposition of complex optical interference filters, the method comprising the steps of:
$\text{\hspace{1em}}\mathrm{cos}\left({\varphi}_{k}\right)=\left(\frac{A\left(1+{r}_{2}^{2}\right)\mathrm{sin}\left(\delta \right)\pm B\text{\hspace{1em}}\mathrm{cos}\left(\delta \right)}{C}\right)$ $\text{\hspace{1em}}A={R}_{f}+{r}_{2}^{4}\left({R}_{f}{R}_{k}\right){R}_{k}+2{r}_{2}^{2}\left(\begin{array}{c}\left(1{R}_{f}\right)\left(1+{R}_{k}\right)\mathrm{cos}\left(2\delta \right)\\ \left(1{R}_{f}{R}_{k}\right)\end{array}\right)$ $\text{\hspace{1em}}B=\sqrt{D\left(1+{r}_{2}^{12}\right)+F\left({r}_{2}^{2}+{r}_{2}^{10}\right)+G\left({r}_{2}^{4}+{r}_{8}^{2}\right)+H\text{\hspace{1em}}{r}_{6}^{2}}$ $\text{\hspace{1em}}C=\mathrm{sin}\left(\delta \right)\left(4{r}_{2}\left(1{R}_{f}\right){R}_{k}^{1/2}\left(2{r}_{2}^{2}\mathrm{cos}\left(2\delta \right)1{r}_{2}^{4}\right)\right)$ $\text{\hspace{1em}}D={\left({R}_{f}{R}_{k}\right)}^{2}$ $\text{\hspace{1em}}F=2\left(\begin{array}{c}{R}_{k}\left(2+{R}_{k}\right)+{R}_{f}^{2}\left(1+2{R}_{k}\right)+2{R}_{f}\left(15{R}_{k}+{R}_{k}^{2}\right)\\ 2\left(1{R}_{f}\right)\left(1{R}_{k}\right)\left({R}_{f}+{R}_{k}\right)\mathrm{cos}\left(2\delta \right)\end{array}\right)$ $G=64{R}_{f}5{R}_{f}^{2}4{R}_{k}+38{R}_{f}{R}_{k}4{R}_{f}^{2}{R}_{k}5{R}_{k}^{2}4{R}_{f}{R}_{k}^{2}6{R}_{f}^{2}{R}_{k}^{2}+8\left(1{R}_{f}^{2}\right)\left(1{R}_{k}^{2}\right)\mathrm{cos}\left(2\delta \right)2{\left(1{R}_{f}\right)}^{2}{\left(1{R}_{k}\right)}^{2}\mathrm{cos}\left(4\delta \right)$ $\text{\hspace{1em}}H=4\left(\begin{array}{c}3+2{R}_{f}^{2}10{R}_{f}{R}_{k}+2{R}_{k}^{2}+3{R}_{f}^{2}{R}_{k}^{2}\\ 2\left(1{R}_{f}\right)\left(1{R}_{k}\right)\left(2+{R}_{f}+{R}_{k}+2{R}_{f}{R}_{k}\right)\mathrm{cos}\left(2\delta \right)+\\ {\left(1{R}_{f}\right)}^{2}{\left(1{R}_{k}\right)}^{2}\mathrm{cos}\left(4\delta \right)\end{array}\right)$
determining from a measurement of intensity reflectance at a topmost interface a phase angle φ at an interface k according to an equation expressed as:
17. The method as in claim 16 , wherein a process control for a deposition system is bootstrapped by detaching the deposition system from all but the topmost interface.
18. The method as in claim 16 , further comprising the step of validating two resultant solutions according to the expression:
19. The method as in claim 16 , further comprising the step of averaging calculated and modeled reflectance and phase values to obtain a new value to be used in all future modeling at a given wavelength.
20. A thin film interference filter system, comprising:
a plurality of stacked films having a determined reflectance;
a modeled monitor curve; and
a topmost layer configured to exhibit a wavelength corresponding to one of the determined reflectance or the modeled monitor curve, the topmost layer being disposed on the plurality of stacked films.
21. The thin film interference filter system as in claim 20 , wherein the topmost layer is a silica film.
22. The thin film interference filter system as in claim 20 , wherein the topmost layer is a niobia film.
Priority Applications (3)
Application Number  Priority Date  Filing Date  Title 

WOUS2005/032420  20040913  
PCT/US2005/032420 WO2006031733A2 (en)  20040913  20050913  Thin film interference filter and bootstrap method for interference filter thin film deposition process control 
US11/684,779 US20070201136A1 (en)  20040913  20070312  Thin Film Interference Filter and Bootstrap Method for Interference Filter Thin Film Deposition Process Control 
Applications Claiming Priority (2)
Application Number  Priority Date  Filing Date  Title 

US11/684,779 US20070201136A1 (en)  20040913  20070312  Thin Film Interference Filter and Bootstrap Method for Interference Filter Thin Film Deposition Process Control 
US12/819,560 US8184371B2 (en)  20050913  20100621  Thin film interference filter and bootstrap method for interference filter thin film deposition process control 
Related Child Applications (1)
Application Number  Title  Priority Date  Filing Date 

US12/819,560 Division US8184371B2 (en)  20040913  20100621  Thin film interference filter and bootstrap method for interference filter thin film deposition process control 
Publications (1)
Publication Number  Publication Date 

US20070201136A1 true US20070201136A1 (en)  20070830 
Family
ID=38443717
Family Applications (2)
Application Number  Title  Priority Date  Filing Date 

US11/684,779 Abandoned US20070201136A1 (en)  20040913  20070312  Thin Film Interference Filter and Bootstrap Method for Interference Filter Thin Film Deposition Process Control 
US12/819,560 Active 20270409 US8184371B2 (en)  20040913  20100621  Thin film interference filter and bootstrap method for interference filter thin film deposition process control 
Family Applications After (1)
Application Number  Title  Priority Date  Filing Date 

US12/819,560 Active 20270409 US8184371B2 (en)  20040913  20100621  Thin film interference filter and bootstrap method for interference filter thin film deposition process control 
Country Status (1)
Country  Link 

US (2)  US20070201136A1 (en) 
Cited By (21)
Publication number  Priority date  Publication date  Assignee  Title 

US20070282647A1 (en) *  20060530  20071206  Freese Robert P  Methods of assessing and designing an application specific measurement system 
US20100302539A1 (en) *  20070330  20101202  Myrick Michael L  Novel multianalyte optical computing system 
US20110090342A1 (en) *  20091015  20110421  University Of South Carolina  Multimode imaging in the thermal infrared for chemical contrast enhancement 
US8049881B2 (en)  20051128  20111101  Halliburton Energy Services, Inc.  Optical analysis system and methods for operating multivariate optical elements in a normal incidence orientation 
US8154726B2 (en)  20051128  20120410  Halliburton Energy Services, Inc.  Optical analysis system and method for real time multivariate optical computing 
US8184371B2 (en)  20050913  20120522  Halliburton Energy Services, Inc.  Thin film interference filter and bootstrap method for interference filter thin film deposition process control 
US8184295B2 (en)  20070330  20120522  Halliburton Energy Services, Inc.  Tablet analysis and measurement system 
US8208147B2 (en)  20051128  20120626  Halliburton Energy Services, Inc.  Method of highspeed monitoring based on the use of multivariate optical elements 
US8212213B2 (en)  20080407  20120703  Halliburton Energy Services, Inc.  Chemicallyselective detector and methods relating thereto 
US8212216B2 (en)  20070330  20120703  Halliburton Energy Services, Inc.  Inline process measurement systems and methods 
WO2012108885A1 (en) *  20110211  20120816  Halliburton Energy Services Inc.  Method for fabrication of a multivariate optical element 
US8283633B2 (en)  20071130  20121009  Halliburton Energy Services, Inc.  Tuning D* with modified thermal detectors 
US8345251B2 (en)  20031231  20130101  Halliburton Energy Services, Inc.  Thinlayer porous optical sensors for gases and other fluids 
US8345234B2 (en)  20051128  20130101  Halliburton Energy Services, Inc.  Self calibration methods for optical analysis system 
US9170154B2 (en)  20060626  20151027  Halliburton Energy Services, Inc.  Data validation and classification in optical analysis systems 
US9182282B2 (en)  20061102  20151110  Halliburton Energy Services, Inc.  Multianalyte optical computing system 
EP2917486A4 (en) *  20131218  20160706  Halliburton Energy Services Inc  Systems and methods to reduce delamination in integrated computational elements used downhole 
US20170262674A1 (en) *  20160311  20170914  Smart Vision Lights  Machine vision systems incorporating polarized lectromagnetic radiation emitters 
US9863870B2 (en)  20110211  20180109  Halliburton Energy Services, Inc.  Method and apparatus to use multiple spectroscopic envelopes to determine components with greater accuracy and dynamic range 
US9885147B2 (en)  20150424  20180206  University Of South Carolina  Reproducible sample preparation method for quantitative stain detection 
US10041866B2 (en)  20150424  20180807  University Of South Carolina  Reproducible sample preparation method for quantitative stain detection 
Families Citing this family (32)
Publication number  Priority date  Publication date  Assignee  Title 

US8661874B2 (en) *  20110428  20140304  Honeywell International Inc.  Photoacoustic detector with background signal correction 
US9222892B2 (en)  20110805  20151229  Halliburton Energy Services, Inc.  Systems and methods for monitoring the quality of a fluid 
US9261461B2 (en)  20110805  20160216  Halliburton Energy Services, Inc.  Systems and methods for monitoring oil/gas separation processes 
US9206386B2 (en)  20110805  20151208  Halliburton Energy Services, Inc.  Systems and methods for analyzing microbiological substances 
US8908165B2 (en)  20110805  20141209  Halliburton Energy Services, Inc.  Systems and methods for monitoring oil/gas separation processes 
US9182355B2 (en)  20110805  20151110  Halliburton Energy Services, Inc.  Systems and methods for monitoring a flow path 
US8912477B2 (en)  20120426  20141216  Halliburton Energy Services, Inc.  Methods and devices for optically determining a characteristic of a substance 
US8823939B2 (en)  20120426  20140902  Halliburton Energy Services, Inc.  Methods and devices for optically determining a characteristic of a substance 
US9019501B2 (en)  20120426  20150428  Halliburton Energy Services, Inc.  Methods and devices for optically determining a characteristic of a substance 
US9702811B2 (en)  20120426  20170711  Halliburton Energy Services, Inc.  Methods and devices for optically determining a characteristic of a substance using integrated computational elements 
US9658149B2 (en)  20120426  20170523  Halliburton Energy Services, Inc.  Devices having one or more integrated computational elements and methods for determining a characteristic of a sample by computationally combining signals produced therewith 
US9013702B2 (en)  20120426  20150421  Halliburton Energy Services, Inc.  Imaging systems for optical computing devices 
US8780352B2 (en)  20120426  20140715  Halliburton Energy Services, Inc.  Methods and devices for optically determining a characteristic of a substance 
US9013698B2 (en)  20120426  20150421  Halliburton Energy Services, Inc.  Imaging systems for optical computing devices 
US8941046B2 (en)  20120426  20150127  Halliburton Energy Services, Inc.  Methods and devices for optically determining a characteristic of a substance 
US8879053B2 (en)  20120426  20141104  Halliburton Energy Services, Inc.  Devices having an integrated computational element and a proximal interferent monitor and methods for determining a characteristic of a sample therewith 
US9383307B2 (en) *  20120426  20160705  Halliburton Energy Services, Inc.  Methods and devices for optically determining a characteristic of a substance 
US9080943B2 (en)  20120426  20150714  Halliburton Energy Services, Inc.  Methods and devices for optically determining a characteristic of a substance 
US9170208B2 (en)  20120831  20151027  Halliburton Energy Services, Inc.  Handheld characteristic analyzer and methods of using the same 
US9103716B2 (en)  20120831  20150811  Halliburton Energy Services, Inc.  Handheld characteristic analyzer and methods of using the same 
US8619256B1 (en)  20120914  20131231  Halliburton Energy Services, Inc.  Systems and methods for monitoring the properties of a fluid cement composition in a flow path 
US9176052B2 (en)  20120914  20151103  Halliburton Energy Services, Inc.  Systems and methods for inspecting and monitoring a pipeline 
US9222896B2 (en)  20120914  20151229  Halliburton Energy Services, Inc.  Systems and methods for inspecting and monitoring a pipeline 
US9086383B2 (en)  20120914  20150721  Halliburton Energy Services, Inc.  Systems and methods for monitoring chemical processes 
US8765061B2 (en)  20120914  20140701  Halliburton Energy Services, Inc.  Systems and methods for inspecting and monitoring a pipeline 
US9228940B2 (en)  20120914  20160105  Halliburton Energy Services, Inc.  Systems, methods, and apparatuses for in situ monitoring of cement fluid compositions and setting processes thereof 
CA2897779A1 (en) *  20130211  20140814  Halliburton Energy Services, Inc.  Fluid analysis system with integrated computation element formed using atomic layer deposition 
BR112015029784A2 (en)  20130709  20170725  Halliburton Energy Services Inc  system, measuring tool and method 
US9733183B2 (en)  20140220  20170815  Halliburton Energy Services, Inc.  Designs for integrated computational elements 
US9708908B2 (en)  20140613  20170718  Halliburton Energy Services, Inc.  Integrated computational element with multiple frequency selective surfaces 
WO2016032437A1 (en)  20140826  20160303  Halliburton Energy Services, Inc.  Systems and methods for in situ monitoring of cement slurry locations and setting processes thereof 
US10344197B2 (en)  20140826  20190709  Halliburton Energy Services, Inc.  Systems and methods for analyzing the characteristics and compositions of cement additives 
Citations (5)
Publication number  Priority date  Publication date  Assignee  Title 

US5103340A (en) *  19890221  19920407  International Business Machines Corporation  Multiplecavity optical filter using change of cavity length 
US5150236A (en) *  19900831  19920922  Bell Communications Research, Inc.  Tunable liquid crystal etalon filter 
US5710655A (en) *  19930721  19980120  Apeldyn Corporation  Cavity thickness compensated etalon filter 
US5799231A (en) *  19960725  19980825  International Business Machines Corporation  Variable index distributed mirror 
US6573999B1 (en) *  20000714  20030603  Nanometrics Incorporated  Film thickness measurements using light absorption 
Family Cites Families (146)
Publication number  Priority date  Publication date  Assignee  Title 

JPS5126720B1 (en) *  19700403  19760807  
US3761724A (en)  19720706  19730925  Resalab Inc  Double beam hydrocarbon gas detector 
JPS57142546A (en)  19810228  19820903  Shimadzu Corp  Infrared multiple reflection type oil concentration measuring apparatus 
US4687337A (en) *  19810902  19870818  The United States Of America As Represented By The Secretary Of The Air Force  Atmospheric Aerosol extinctiometer 
US4499378A (en) *  19820309  19850212  Horiba, Ltd.  Infrared radiation gas analyzer 
US4607914A (en) *  19831012  19860826  Environmental Research Institute Of Michigan  Optical system design techniques using holographic optical element 
JPS60135757A (en)  19831223  19850719  Hochiki Corp  Gas sensor 
US4595832A (en) *  19840330  19860617  The United States Of America As Represented By The United States Department Of Energy  Thermal sensor with an improved coating 
US5259381A (en)  19860818  19931109  PhysioControl Corporation  Apparatus for the automatic calibration of signals employed in oximetry 
GB8712582D0 (en)  19870528  19870701  Neotronics Ltd  Acidic gas sensors 
JPH01167646A (en) *  19871223  19890703  Tekune Yoko:Kk  Electric capacity type hygrometer 
US5622868A (en) *  19890427  19970422  Microbiological Research Authority Camr (Centre For Applied Microbiology & Research)  Analytical apparatus utilizing a colorimetric or other optically detectable effect 
US4981332A (en) *  19890629  19910101  Hughes Aircraft Company  Dispersioncompensated windshield hologram virtual image display 
US5090807A (en) *  19900110  19920225  Environmental Research Institute Of Michigan  Real time optical predetection processing of multispectral image data 
US5137364A (en) *  19910131  19920811  Mccarthy Cornelius J  Optical spectral analysis apparatus 
US5223715A (en) *  19910920  19930629  Amoco Corporation  Process for spectrophotometric analysis 
US5406082A (en) *  19920424  19950411  Thiokol Corporation  Surface inspection and characterization system and process 
EP0600334B1 (en) *  19921201  20020206  Ulrich Walter Glaus  Method for the determination of materials and/or their properties and apparatus therefor 
US5946088A (en) *  19940503  19990831  Pfizer Inc.  Apparatus for mixing and detecting online homogeneity 
US5504332A (en) *  19940826  19960402  Merck & Co., Inc.  Method and system for determining the homogeneity of tablets 
GB2302731B (en)  19950616  19971112  Sun Electric Uk Ltd  Method and apparatus for gas monitoring 
US5905571A (en) *  19950830  19990518  Sandia Corporation  Optical apparatus for forming correlation spectrometers and optical processors 
US5760399A (en) *  19951002  19980602  Foss Nirsystems, Inc.  Measurement of transmission spectra of pharmaceutical tablets 
WO1997015819A1 (en)  19951025  19970501  University Of Washington  Surface plasmon resonance light pipe sensor 
US5641962A (en) *  19951205  19970624  Exxon Research And Engineering Company  Non linear multivariate infrared analysis method (LAW362) 
US7028899B2 (en) *  19990607  20060418  Metrologic Instruments, Inc.  Method of specklenoise pattern reduction and apparatus therefore based on reducing the temporalcoherence of the planar laser illumination beam before it illuminates the target object by applying temporal phase modulation techniques during the transmission of the plib towards the target 
US5747806A (en) *  19960202  19980505  Instrumentation Metrics, Inc  Method and apparatus for multispectral analysis in noninvasive nir spectroscopy 
US5734098A (en) *  19960325  19980331  Nalco/Exxon Energy Chemicals, L.P.  Method to monitor and control chemical treatment of petroleum, petrochemical and processes with online quartz crystal microbalance sensors 
DE59703919D1 (en)  19960816  20010802  Zeptosens Ag Witterswil  Optical detection device 
JP2002506358A (en)  19960906  20020226  メルク アンド カンパニー，インコーポレーテッド  Customerspecific packing line 
US5781289A (en) *  19961105  19980714  Sabsabi; Mohamad  Method and apparatus for rapid in situ analysis of preselected components of homogeneous solid compositions, especially pharmaceutical compositions 
AU5547798A (en)  19961115  19980610  Optosens Optische Spektroskopie Und Sensortechnik Gmbh  Method and device for combined absorption and reflectance spectroscopy 
EP0855592B1 (en)  19970125  20060712  Siemens Schweiz AG  Optoacoustic gas sensor 
US6040914A (en) *  19970610  20000321  New Focus, Inc.  Simple, low cost, laser absorption sensor system 
AU8164898A (en) *  19970627  19990119  Baker Hughes Incorporated  Drilling system with sensors for determining properties of drilling fluid downhole 
US6198531B1 (en) *  19970711  20010306  University Of South Carolina  Optical computational system 
US6982431B2 (en) *  19980831  20060103  Molecular Devices Corporation  Sample analysis systems 
US5941821A (en) *  19971125  19990824  Trw Inc.  Method and apparatus for noninvasive measurement of blood glucose by photoacoustics 
US7267948B2 (en) *  19971126  20070911  UtBattelle, Llc  SERS diagnostic platforms, methods and systems microarrays, biosensors and biochips 
US5939717A (en) *  19980129  19990817  Schlumberger Technology Corporation  Methods and apparatus for determining gasoil ratio in a geological formation through the use of spectroscopy 
US6424407B1 (en) *  19980309  20020723  Otm Technologies Ltd.  Optical translation measurement 
EP1086366B1 (en)  19980612  20100414  Radiometer Medical ApS  A method in quality control of a spectrophotometer 
JP3020154B2 (en) *  19980612  20000315  東京大学長  Method for producing porous diamond body 
US6137108A (en)  19980617  20001024  Foss Nirsystems Incorporated  Instrument and method for spectroscopic analysis by reflectance and transmittance 
US5946089A (en) *  19980713  19990831  Jordan Valley Applied Radiation Ltd.  Plasma spectrometer with shutter assembly 
SE9802690D0 (en)  19980807  19980807  Astra Ab  Mixing Apparatus 
US6306658B1 (en)  19980813  20011023  Symyx Technologies  Parallel reactor with internal sensing 
EP1105710B1 (en) *  19980814  20111228  OnSite Analysis, Inc.  Method of normalizing the spectrum intensity measured by optical emission spectrometry 
US6430513B1 (en) *  19980904  20020806  Perkinelmer Instruments Llc  Monitoring constituents of an animal organ using statistical correlation 
US6630663B2 (en)  19981021  20031007  Raytheon Company  Miniature ion mobility spectrometer 
US6690464B1 (en) *  19990219  20040210  Spectral Dimensions, Inc.  Highvolume online spectroscopic composition testing of manufactured pharmaceutical dosage units 
US6341257B1 (en) *  19990304  20020122  Sandia Corporation  Hybrid least squares multivariate spectral analysis methods 
US7123844B2 (en)  19990406  20061017  Myrick Michael L  Optical computational system 
US6529276B1 (en) *  19990406  20030304  University Of South Carolina  Optical computational system 
US6304854B1 (en)  19990909  20011016  Dunhill Holdings, Corp.  System and method for providing a comparable branded product based on a current branded product for noncomparison shopped products 
US6347131B1 (en) *  19991102  20020212  Cintex Limited  Noncontact volume measurement 
US6507401B1 (en) *  19991202  20030114  Aps Technology, Inc.  Apparatus and method for analyzing fluids 
US6748334B1 (en) *  19991206  20040608  Jorge E. Perez  Specialty gas analysis system 
FR2803081B1 (en) *  19991228  20021206  Trixell Sas  Method for temperature compensation of an image sensor 
US6887703B2 (en)  20000216  20050503  Arturus Bioscience, Inc.  Transfer film for laser microcapture 
WO2001060507A1 (en)  20000217  20010823  Astrazeneca Uk Limited  Mixing apparatus and method 
US6420708B2 (en) *  20000310  20020716  Wilks Enterprise, Inc.  Spectroscopy analyzer using a detector array 
US6760359B2 (en) *  20000428  20040706  Photodigm, Inc.  Gratingoutcoupled surfaceemitting lasers with flared gain regions 
DE10021378A1 (en) *  20000502  20011108  Leica Microsystems Wetzlar Gmbh  Optical measuring arrangement with an ellipsometer 
DE60134844D1 (en) *  20000510  20080828  Dsm Ip Assets Bv  Method for identifying objects by means of an optical spectrometer and a transport system 
FR2809816B1 (en) *  20000530  20030418  Gaz De France  Method and device for detecting gas leaks 
US6437326B1 (en) *  20000627  20020820  Schlumberger Technology Corporation  Permanent optical sensor downhole fluid analysis systems 
WO2002016976A2 (en) *  20000821  20020228  3M Innovative Properties Company  Loss enhanced reflective optical filters 
US6595427B1 (en)  20000831  20030722  Polaroid Corporation  Method and apparatus for encoding and decoding information in a nonvisible manner 
US7138156B1 (en)  20000926  20061121  Myrick Michael L  Filter design algorithm for multivariate optical computing 
US6476384B1 (en)  20001010  20021105  Schlumberger Technology Corporation  Methods and apparatus for downhole fluids analysis 
JP3674504B2 (en) *  20001211  20050720  ウシオ電機株式会社  Spectral reflectance measuring apparatus and spectral reflectance measuring method 
EP1229356A3 (en) *  20010131  20040121  Planar Systems, Inc.  Methods and apparatus for the production of optical filters 
US6894772B2 (en) *  20010212  20050517  Analytical Spectral Devices  System and method for grouping reflectance data 
US6853447B2 (en) *  20010212  20050208  Analytical Spectral Devices, Inc.  System and method for the collection of spectral image data 
US6765212B2 (en) *  20010212  20040720  Analytical Spectral Devices, Inc.  System and method for combining reflectance data 
US6667802B2 (en)  20010212  20031223  Analytical Spectral Devices, Inc.  System and method for selfreferencing calibration 
US7023553B2 (en) *  20010907  20060404  Wallac Oy  Intelligent instrumentation with changeable optical components 
GB2379976B (en)  20010920  20050202  Ndc Infrared Eng Ltd  Optical sampling window 
US6891606B2 (en) *  20011011  20050510  Baker Hughes Incorporated  Realtime online sensing and control of mineral scale deposition from formation fluids 
US6870629B1 (en) *  20011029  20050322  Precision Photonics Corporation  Optical frequency sweep control and readout by using a phase lock 
CN1659424A (en) *  20020306  20050824  埃斯柏克特瑞克斯公司  Method and apparatus for radiation encoding and analysis 
US6771369B2 (en) *  20020312  20040803  Analytical Spectral Devices, Inc.  System and method for pharmacy validation and inspection 
US6741348B2 (en)  20020429  20040525  The Curators Of The University Of Missouri  Ultrasensitive spectrophotometer 
US20030223068A1 (en)  20020604  20031204  Baker Hughes Incorporated  Method and apparatus for a high resolution downhole spectrometer 
US6839137B2 (en) *  20020719  20050104  Exxonmobil Research And Engineering Company  Asphaltene aggregation in petroleum oil mixtures determined by small angle light scattering 
EP1534368A1 (en)  20020725  20050601  Glaxo Group Limited  Medicament dispenser 
JP2004149909A (en) *  20021101  20040527  Advance:Kk  Processes for forming silver mirror film and coated film comprising the same 
EP1558931B1 (en) *  20021104  20070613  Transform Pharmaceuticals, Inc.  Analysis of pharmaceutical solubility and stability 
CA2520261A1 (en) *  20030224  20050106  Cdex, Inc.  System and methods for detection and identification of chemical substances 
US7671973B2 (en) *  20031219  20100302  Koninklijke Philips Electronics N.V.  Optical analysis system using multivariate optical elements 
CN100443867C (en) *  20021219  20081217  皇家飞利浦电子股份有限公司  Optical analysis system 
US6897951B2 (en) *  20030214  20050524  Raman Systems, Inc.  Probe assemblies for Raman spectroscopy 
US20040227086A1 (en)  20030307  20041118  Pfizer Inc  Method and device for NIR reflection spectroscopy 
US7173239B2 (en) *  20030314  20070206  Baker Hughes Incorporated  Method and apparatus for downhole quantification of methane using near infrared spectroscopy 
US20050032235A1 (en) *  20030618  20050210  Srinivas Tummala  Method of monitoring the blending of a mixture 
US6952267B2 (en)  20030707  20051004  Cymer, Inc.  Method and apparatus for measuring bandwidth of a laser output 
JP2007515653A (en) *  20031222  20070614  コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ  Optical analysis system, blood analysis system, and method for determining the amplitude of a principal component 
JP2007523321A (en)  20031231  20070816  ユニヴァーシティー オブ サウスカロライナ  Thin layer porous optical sensor for gases and other fluids 
JP4904165B2 (en) *  20040115  20120328  デポジッション サイエンス インク  Method and apparatus for monitoring optical properties of thin films during a deposition process 
WO2005111580A1 (en)  20040507  20051124  Optiscan Biomedical Corporation  Sample element with fringingreduction capabilities 
US20050288906A1 (en)  20040629  20051229  Drennen James K Iii  Spectroscopic pharmacy verification and inspection system 
US7271883B2 (en)  20040713  20070918  Newell Benjamin E  Refractive index sensor using internally reflective light beams 
EP1626278A3 (en)  20040803  20060621  OnChip Cellomics Consortium  Cellomics system 
US7411729B2 (en)  20040812  20080812  Olympus Corporation  Optical filter, method of manufacturing optical filter, optical system, and imaging apparatus 
US20080309930A1 (en)  20040826  20081218  Koninklijke Philips Electronics N.V.  Calibration for Spectroscopic Analysis 
US7218822B2 (en) *  20040903  20070515  Chemimage Corporation  Method and apparatus for fiberscope 
US20070201136A1 (en)  20040913  20070830  University Of South Carolina  Thin Film Interference Filter and Bootstrap Method for Interference Filter Thin Film Deposition Process Control 
JP2008512730A (en)  20040913  20080424  ユニヴァーシティー オブ サウスカロライナ  Thin film interference filter and bootstrap method for controlling deposition process of interference filter 
US8240189B2 (en)  20041004  20120814  Halliburton Energy Services, Inc.  Thermal selectivity multivariate optical computing 
US7993276B2 (en) *  20041015  20110809  Pulse Tracer, Inc.  Motion cancellation of optical input signals for physiological pulse measurement 
EP1805499A4 (en)  20041021  20100721  Optiscan Biomedical Corp  Method and apparatus for determining an analyte concentration in a sample having interferent 
US20060093523A1 (en) *  20041029  20060504  Hyperteq, Lp  System, method and apparatus for mudgas extraction, detection and analysis thereof 
WO2006063094A1 (en)  20041209  20060615  Caleb Brett Usa Inc.  In situ optical computation fluid analysis system and method 
CN101080617A (en) *  20041215  20071128  皇家飞利浦电子股份有限公司  Calibration of optical analysis making use of multivariate optical elements 
US7302148B2 (en) *  20050113  20071127  Komag, Inc.  Test head for optically inspecting workpieces 
JP2006214758A (en)  20050201  20060817  Denso Corp  Infrared detector 
US7330746B2 (en)  20050607  20080212  Chem Image Corporation  Noninvasive biochemical analysis 
US7279678B2 (en) *  20050815  20071009  Schlumber Technology Corporation  Method and apparatus for composition analysis in a logging environment 
US7834999B2 (en)  20051128  20101116  University Of South Carolina  Optical analysis system and optical train 
EP1969326A1 (en)  20051128  20080917  Ometric Corporation  Optical analysis system and method for real time multivariate optical computing 
WO2007062202A1 (en)  20051128  20070531  University Of South Carolina  Novel multivariate optical elements for optical analysis system 
EP1955032A1 (en)  20051128  20080813  University of South Carolina  Process for selecting spectral elements and components for optical analysis systems 
WO2007064579A1 (en) *  20051128  20070607  University Of South Carolina  Optical analysis system and elements to isolate spectral region 
WO2007061436A1 (en)  20051128  20070531  University Of South Carolina  Self calibration methods for optical analysis system 
US8208147B2 (en)  20051128  20120626  Halliburton Energy Services, Inc.  Method of highspeed monitoring based on the use of multivariate optical elements 
US20070166245A1 (en) *  20051128  20070719  Leonard Mackles  Propellant free foamable toothpaste composition 
US7458258B2 (en) *  20051216  20081202  Schlumberger Technology Corporation  Methods and apparatus for oil composition determination 
US7623233B2 (en) *  20060310  20091124  Ometric Corporation  Optical analysis systems and methods for dynamic, highspeed detection and realtime multivariate optical computing 
US8027855B2 (en)  20060530  20110927  Halliburton Energy Services Inc.  Methods of assessing and designing an application specific measurement system 
US20070294094A1 (en)  20060620  20071220  Ometric Corporation  Methods of customizing, licensing and sustaining a technology option to meet a customer requirement 
EP2033196A2 (en)  20060626  20090311  University of South Carolina  Data validation and classification in optical analysis systems 
US7652767B2 (en) *  20061019  20100126  Sporian Microsystems, Inc.  Optical sensor with chemically reactive surface 
WO2008057905A2 (en)  20061102  20080515  Ometric Corporation  Selfcontained multivariate optical computing and analysis system 
WO2008057912A2 (en)  20061102  20080515  University Of South Carolina  Multianalyte optical computing system 
WO2008057913A2 (en)  20061102  20080515  University Of South Carolina  Improved signal processing for optical computing system 
US8213012B2 (en)  20061102  20120703  Halliburton Energy Services, Inc.  Stability for optical computing system 
US20080111064A1 (en) *  20061110  20080515  Schlumberger Technology Corporation  Downhole measurement of substances in earth formations 
WO2008106391A1 (en) *  20070228  20080904  University Of South Carolina  Design of multivariate optical elements for nonlinear calibration 
US8184295B2 (en)  20070330  20120522  Halliburton Energy Services, Inc.  Tablet analysis and measurement system 
WO2008121715A1 (en) *  20070330  20081009  Ometric Corporation  Inline process measurement systems and methods 
US8213006B2 (en)  20070330  20120703  Halliburton Energy Services, Inc.  Multianalyte optical computing system 
WO2009082418A2 (en)  20071012  20090702  RealTime Analyzers, Inc.  Method and apparatus for determining properties of fuels 
US8283633B2 (en) *  20071130  20121009  Halliburton Energy Services, Inc.  Tuning D* with modified thermal detectors 
US8212213B2 (en)  20080407  20120703  Halliburton Energy Services, Inc.  Chemicallyselective detector and methods relating thereto 
US7853104B2 (en)  20090105  20101214  Yokogawa Electric Corporation  Bidirectional optical module and optical time domain reflectometer 

2007
 20070312 US US11/684,779 patent/US20070201136A1/en not_active Abandoned

2010
 20100621 US US12/819,560 patent/US8184371B2/en active Active
Patent Citations (5)
Publication number  Priority date  Publication date  Assignee  Title 

US5103340A (en) *  19890221  19920407  International Business Machines Corporation  Multiplecavity optical filter using change of cavity length 
US5150236A (en) *  19900831  19920922  Bell Communications Research, Inc.  Tunable liquid crystal etalon filter 
US5710655A (en) *  19930721  19980120  Apeldyn Corporation  Cavity thickness compensated etalon filter 
US5799231A (en) *  19960725  19980825  International Business Machines Corporation  Variable index distributed mirror 
US6573999B1 (en) *  20000714  20030603  Nanometrics Incorporated  Film thickness measurements using light absorption 
Cited By (30)
Publication number  Priority date  Publication date  Assignee  Title 

US8345251B2 (en)  20031231  20130101  Halliburton Energy Services, Inc.  Thinlayer porous optical sensors for gases and other fluids 
US8184371B2 (en)  20050913  20120522  Halliburton Energy Services, Inc.  Thin film interference filter and bootstrap method for interference filter thin film deposition process control 
US8345234B2 (en)  20051128  20130101  Halliburton Energy Services, Inc.  Self calibration methods for optical analysis system 
US8049881B2 (en)  20051128  20111101  Halliburton Energy Services, Inc.  Optical analysis system and methods for operating multivariate optical elements in a normal incidence orientation 
US8154726B2 (en)  20051128  20120410  Halliburton Energy Services, Inc.  Optical analysis system and method for real time multivariate optical computing 
US8208147B2 (en)  20051128  20120626  Halliburton Energy Services, Inc.  Method of highspeed monitoring based on the use of multivariate optical elements 
US8358418B2 (en)  20051128  20130122  Halliburton Energy Services, Inc.  Optical analysis system for dynamic realtime detection and measurement 
US8027855B2 (en) *  20060530  20110927  Halliburton Energy Services Inc.  Methods of assessing and designing an application specific measurement system 
US20070282647A1 (en) *  20060530  20071206  Freese Robert P  Methods of assessing and designing an application specific measurement system 
US9170154B2 (en)  20060626  20151027  Halliburton Energy Services, Inc.  Data validation and classification in optical analysis systems 
US9182282B2 (en)  20061102  20151110  Halliburton Energy Services, Inc.  Multianalyte optical computing system 
US8213006B2 (en)  20070330  20120703  Halliburton Energy Services, Inc.  Multianalyte optical computing system 
US20100302539A1 (en) *  20070330  20101202  Myrick Michael L  Novel multianalyte optical computing system 
US8184295B2 (en)  20070330  20120522  Halliburton Energy Services, Inc.  Tablet analysis and measurement system 
US8212216B2 (en)  20070330  20120703  Halliburton Energy Services, Inc.  Inline process measurement systems and methods 
US8283633B2 (en)  20071130  20121009  Halliburton Energy Services, Inc.  Tuning D* with modified thermal detectors 
US8212213B2 (en)  20080407  20120703  Halliburton Energy Services, Inc.  Chemicallyselective detector and methods relating thereto 
US20110090342A1 (en) *  20091015  20110421  University Of South Carolina  Multimode imaging in the thermal infrared for chemical contrast enhancement 
US8823802B2 (en)  20091015  20140902  University Of South Carolina  Multimode imaging in the thermal infrared for chemical contrast enhancement 
GB2501641A (en) *  20110211  20131030  Halliburton Energy Serv Inc  Method for fabrication of a multivariate optical element 
CN103354915A (en) *  20110211  20131016  哈利伯顿能源服务公司  Method for fabrication of a multivariate optical element 
WO2012108885A1 (en) *  20110211  20120816  Halliburton Energy Services Inc.  Method for fabrication of a multivariate optical element 
US10094955B2 (en)  20110211  20181009  Halliburton Energy Services, Inc.  Method for fabrication of a multivariate optical element 
US9863870B2 (en)  20110211  20180109  Halliburton Energy Services, Inc.  Method and apparatus to use multiple spectroscopic envelopes to determine components with greater accuracy and dynamic range 
GB2501641B (en) *  20110211  20180221  Halliburton Energy Services Inc  Method for fabrication of a multivariate optical element 
EP2917486A4 (en) *  20131218  20160706  Halliburton Energy Services Inc  Systems and methods to reduce delamination in integrated computational elements used downhole 
US9885147B2 (en)  20150424  20180206  University Of South Carolina  Reproducible sample preparation method for quantitative stain detection 
US10041866B2 (en)  20150424  20180807  University Of South Carolina  Reproducible sample preparation method for quantitative stain detection 
US10067069B2 (en) *  20160311  20180904  Smart Vision Lights  Machine vision systems incorporating polarized electromagnetic radiation emitters 
US20170262674A1 (en) *  20160311  20170914  Smart Vision Lights  Machine vision systems incorporating polarized lectromagnetic radiation emitters 
Also Published As
Publication number  Publication date 

US8184371B2 (en)  20120522 
US20100305741A1 (en)  20101202 
Similar Documents
Publication  Publication Date  Title 

Ulrich et al.  Measurement of thin film parameters with a prism coupler  
JP4679365B2 (en)  Processing end point detection method using broadband reflected light measurement method  
US6979578B2 (en)  Process endpoint detection method using broadband reflectometry  
Woollam et al.  Overview of variableangle spectroscopic ellipsometry (VASE): I. Basic theory and typical applications  
Nestell et al.  Derivation of optical constants of metals from thinfilm measurements at oblique incidence  
EP0577399B1 (en)  Apparatus and method for performing thin film layer thickness metrology on a thin film layer having shape deformations and local slope variations  
US5365340A (en)  Apparatus and method for measuring the thickness of thin films  
US7009680B2 (en)  Narrow band tunable filter with integrated detector  
Macleod  Turning value monitoring of narrowband alldielectric thinfilm optical filters  
US7046442B2 (en)  Wire grid polarizer  
US7203001B2 (en)  Optical retarders and related devices and systems  
CN100398694C (en)  Method and device for controlling thickness of optical film, insulation multilayer film and making device  
Tompkins et al.  Spectroscopic Ellipsometry: Practical Application to Thin Film Characterization  
Vriens et al.  Optical constants of absorbing thin solid films on a substrate  
US5354575A (en)  Ellipsometric approach to antireflection coatings of semiconductor laser amplifiers  
KR100508007B1 (en)  Process for forming a thin film and apparatus therefor  
US5333049A (en)  Apparatus and method for interferometrically measuring the thickness of thin films using full aperture irradiation  
Connell et al.  Amorphous germanium III. Optical properties  
US7196793B2 (en)  Method for analyzing thinfilm layer structure using spectroscopic ellipsometer  
Rovira et al.  Analysis of specular and textured SnO 2: F films by high speed fourparameter Stokes vector spectroscopy  
Liddiard  Application of interferometric enhancement to selfabsorbing thin film thermal IR detectors  
Konstantinov et al.  Analysis of errors in thinfilm optical parameters derived from spectrophotometric measurements at normal light incidence  
US4707611A (en)  Incremental monitoring of thin films  
US20040263983A1 (en)  Antireflective coatings and structures  
DE60314706T2 (en)  Wire grid polarizer 
Legal Events
Date  Code  Title  Description 

AS  Assignment 
Owner name: UNIVERSITY OF SOUTH CAROLINA, SOUTH CAROLINA Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:MYRICK, MICHAEL L;REEL/FRAME:018995/0283 Effective date: 20070228 

STCB  Information on status: application discontinuation 
Free format text: ABANDONED  FAILURE TO RESPOND TO AN OFFICE ACTION 

AS  Assignment 
Owner name: HALLIBURTON ENERGY SERVICES, INC., TEXAS Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:UNIVERSITY OF SOUTH CAROLINA;REEL/FRAME:025790/0442 Effective date: 20110121 