**MgAl**spinel is ∼

_{2}O_{4}**40%**greater in the range of 3.8 to 5.0

*μ*m than the value predicted by the computer code OPTIMATR®, which has been used for window and dome design for more than 20 years. As a result, spinel and

*-plane sapphire windows designed to support the same external pressure with the same probability of survival have approximately the same infrared absorptance in the range 3.8 to 5.0*

**a***μ*m.

*-Plane sapphire has greater absorptance than spinel in the range 3.8 to 5.0*

**c***μ*m. Spinel has two weak absorption bands near 1.8 and 3.0

*μ*m. At 1.064

*μ*m, the laser calorimetric absorption coefficient of spinel is 10 to 50 times greater than that of sapphire. New measurements of specific heat capacity, thermal expansion, thermal conductivity, elastic constants, and refractive index (including

**d**) of spinel are reported.

*n*/dT## 1.

## Introduction

Sapphire, spinel, and aluminum oxynitride (ALON) are the most durable, commercially available midwave infrared (3 to 5 *μ*m) window materials that are candidates for applications requiring environmental durability, such as sand and rain erosion resistance or the ability to operate at elevated temperature.^{1} All three materials transmit from about 5 *μ*m in the infrared to at least as far as 0.3 *μ*m in the ultraviolet. Wickersham and LaFever^{2} observed as early as 1960 that spinel has superior transmittance to that of sapphire near 5 *μ*m. Figure 1 compares infrared transmittance of these materials, with the region 4 to 5 *μ*m highlighted. The long wavelength cutoff between 5 and 6 *μ*m increases in the order $\mathrm{ALON}<c\text{-plane sapphire}<a\text{-plane sapphire}<\text{spinel}$.

For applications requiring maximum transmittance or minimum optical emittance in the 4 to 5 *μ*m region, it appears from Fig. 1 that spinel would be the best choice if it is the same thickness as sapphire and has sufficient mechanical strength. Sapphire is mechanically stronger than spinel,^{1} so a thinner sapphire window might be optically superior to a thicker spinel window and still meet mechanical requirements. Of the two commonly available crystal orientations of sapphire, $a$-plane material has less absorption than $c$-plane material.^{3} However, $a$-plane sapphire has a refractive index anisotropy of 0.001 in the plane of the window, whereas $c$-plane sapphire has an isotropic refractive index in the plane of the window. Optical anisotropy might be an issue for some applications.^{4} $a$-Plane sapphire is currently available in larger sizes than $c$-plane material, so the choice of crystal plane can be dictated by the required window size.

Careful transmittance measurements^{5}^{,}^{6} of these materials in the 1980s provided input to a computer code called OPTIMATR®, which was available in the 1990s,^{7}^{,}^{8} but is no longer commercially available. This computer code implements a multioscillator, multiphonon model for the complex refractive index of each material over a range of wavelengths and temperatures.^{9} The code was used to generate the optical constants $n$ and $k$ that are listed for ALON,^{10} spinel,^{11} and sapphire^{3} in Palik’s *Handbook of Optical Constants of Solids*. The code is validated in the region of intrinsic absorption by these materials, but is not meaningful in the transparent window region where absorption can be governed by impurities and defects.

We are interested in the relative merits of spinel and sapphire for large sensor windows with lateral dimensions $>30\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{cm}$ that need high transmittance in the 4 to 5 *μ*m region and must also withstand external pressure. We measured the absorption of some of the highest optical quality currently available spinel and sapphire near 5 *μ*m by laser calorimetry to compare with absorption coefficients measured in the 1980s. Spinel was taken from large area blanks that could be used to make large windows. This paper presents the results of laser calorimetry, as well as new measurements of spectral transmittance, specific heat capacity, refractive index, thermal expansion, thermal conductivity, Young’s modulus, and Poisson’s ratio of spinel. The result is a revision to the OPTIMATR absorption coefficient of spinel in the 4 to 5 *μ*m region. We discuss window design trade between optical absorption and mechanical reliability.

## 2.

## Experimental

Polycrystalline spinel with a grain size of several hundred microns was provided by TA&T (Annapolis, MD) in the form of two polished blocks stated to represent high-optical-quality material. This material is representative of the quality and grain size commercially available today for large area windows with high optical quality. The blocks were cut from $35\times 40\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{cm}$ plates made by hot pressing followed by hot isostatic pressing. One polished block with a finished thickness of 5.98 mm was cut from a blank with a thickness of 7.4 mm. Another polished block with a finished thickness of 12.76 mm was cut from a blank with a thickness of 14.7 mm. Thermal, mechanical, and optical properties reported in this paper come from the TA&T material. For careful measurements of infrared transmittance, additional samples of polished spinel were provided by TA&T and by MER Corp. (Tucson, AZ).

$a$-Plane ($\pm 2\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{deg}$) and $c$-plane ($\pm 2\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{deg}$) sapphire were provided by GT Advanced Technology (Salem, MA) (formerly Crystal Systems) from their highest-quality, strain-free, single-crystal material grown by the heat-exchanger method. Boules were grown from aluminum oxide melt stock with a stated purity of 99.996%. Disks for calorimetry were machined and polished to a 10-5 scratch-dig surface finish by Crystal Systems.

Specific heat capacity (${c}_{p}$) at constant pressure ($P$) is the heat ($H$) required to raise the temperature of 1 g of substance by 1 K. Sapphire is a reference material for this measurement (ASTM E1269).^{12} Duplicate samples of spinel were measured by differential scanning calorimetry at the Thermophysical Properties Research Laboratory (TPRL), West Lafayette, Indiana. Two small compartments were heated electrically at $20\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{K}/\mathrm{min}$ while temperature and electrical heat input were monitored. Samples were disks with a diameter of 5.8 mm and thickness of 1.5 mm with one smooth ground face for thermal contact. Spinel or a sapphire standard was placed in one compartment, while the reference compartment was empty. The differential power required to heat the sample or standard at the same rate as the blank reference pan was determined. From the mass of the standard or sample, the differential power required to heat each item, and the known specific heat capacity of the sapphire standard, the specific heat capacity of the sample was computed. Samples were measured in a PerkinElmer DSC 4 differential scanning calorimeter from 23 to 500°C under Ar and then in a Netzsch DSC 404 from 23 to 1000°C under Ar. Results from two samples, which had an average difference of 0.3 to 1% at each temperature, were averaged to obtain the specific heat capacity.

Thermal expansion of spinel was measured at the TRPL by ASTM standard testing procedure E228 with a dual push-rod dilatometer (Theta Dilatronics II). Test specimens were $51\times 6\times 6\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mm}$ parallelepipeds with flat, parallel, ground ends. The differential expansion between a test specimen and a reference material was measured at 10-deg intervals from 20 to 1000°C during warming and then from 1000 to 200°C during cooling. Expansion of the test specimen was computed from this differential expansion and the known expansion of the standard. Six standards from the National Institute of Standards and Technology (NIST) included materials with low, moderate, and large expansions. For calibration and checkout, one standard was measured against another standard.

Thermal diffusivity of spinel was measured at the TRPL by the laser flash method (ASTM E1461). In this method, the front face of a disk is subjected to a short laser burst and the rear face temperature rise is recorded and analyzed. Ceramic disk samples with parallel, fine ground faces had a diameter of 12.7 mm and thickness of 2.54 mm. The apparatus consists of a Korad K2 1.064 *μ*m Nd:YAG laser, a high-vacuum system including a bell jar with windows for viewing the sample, a tantalum or stainless steel tube heater surrounding a sample holding assembly, a thermocouple or infrared detector, and a digital data acquisition system capable of taking data at $\ge 40\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mu \mathrm{s}$. Both surfaces of spinel were coated with a thin layer of black stove paint to absorb the laser energy.

Young’s modulus and Poisson’s ratio were measured at the University of Dayton Research Institute (Dayton, OH) at 21°C by impulse excitation of vibration (ASTM C 1259) with duplicate $38\text{-}\mathrm{mm}\text{\hspace{0.17em}}\text{diameter}\times 2\text{-}\mathrm{mm}\text{\hspace{0.17em}}\text{thick}$ disks. Disk faces were optically polished and edges were ground. Measurements were made with a GrindoSonic Model MK3s (J. W. Lemmens, Leuven, Belgium) fitted with a microphone. Specimens supported on packing foam were excited by a steel ball bearing on a polymer rod. The flexural frequency was measured at the circumference of the disk after impact at the center. The torsional frequency was found by impacting the nodal circle whose diameter was 0.681 times the specimen diameter and measuring the vibration 90 deg away on the nodal circle. Testing was complete when five consecutive readings agreed within 1%. Young’s modulus ($E$) and Poisson’s ratio ($\nu $) were calculated using the EMOD software supplied with the GrindoSonic.

Transmittance from 1.5 to 8 *μ*m was measured with a nitrogen-purged PerkinElmer Spectrum GX Fourier transform infrared (FTIR) spectrophotometer at the Naval Air Warfare Center. The instrument uses a 1350 K wire coil source, a deuterated triglycine sulfate detector, and a KBr beamsplitter. Each spectrum is an average of 32 scans with $4\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{cm}}^{-1}$ resolution at a wavelength of 2.5 *μ*m. Photometric accuracy is $\pm 0.25\%$. The 100% transmittance was measured before each sample scan. A new background was collected if the 100% transmittance changed by more than $\pm 0.25\%$. Transmittance from 0.25 to 2.5 *μ*m was measured with a 150-mm diameter integrating sphere on a PerkinElmer Lambda 950 UV-Vis-NIR spectrophotometer with a photometric accuracy of $\pm 0.25\%\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{T}$. The FTIR values matched the integrating sphere values to better than $\pm 0.1\%$ in the wavelength region of 1.5 to 2.5 *μ*m where the two instruments overlap.

Total integrated infrared optical scatter in the forward and reverse hemispheres was measured at the Naval Air Warfare Center with a Coblentz sphere that captured light scattered between 2.5 and 70 deg from the direction of the incident laser beam.^{13} The light source was an optical parametric oscillator laser from JP Innovations (Monroe, WA) tunable from 3.46 to 4.41 *μ*m. The idler output is $\sim 6\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mW}$. A 4-mm thick Ge window at the laser output blocks both the pump and signal radiation. Total integrated forward and back scatter of spinel from TA&T and MER were each in the range of 0.1 to 0.4% at 3.5 *μ*m for samples with a thickness of 7 to 13 mm.

Refractive index of one TA&T spinel prism was measured at the Air Force Research Laboratory by the method of minimum deviation.^{14}^{,}^{15} The prism had a $10\times 10\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mm}$ base and rectangular $10\times 13\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mm}$ faces. The apex angle was $45.1075\pm 0.0014\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{deg}$, where the uncertainty is the standard deviation of 10 measurements made with an autocollimator. Five measurements of the angle of minimum deviation were made on a Moller-Wedel divided circle spectrometer at each wavelength and the average was used to calculate the refractive index. For thermo-optic coefficient measurements, the advantage of the method is that it allows for direct measurement of $\mathrm{d}n/\mathrm{d}T$ (change in refractive index $n$ with respect to temperature $T$) without having to subtract the contribution of the thermal expansion coefficient from the overall change in optical path, which is measured interferometrically.^{16}^{,}^{17} Refractive indices of standards were within 0.0001 of NIST measurements.^{18}^{,}^{19} The estimated error in the refractive index is $\sim 0.0001$ and the error in the thermo-optic coefficients, $\mathrm{d}n/\mathrm{d}T$, is $\pm 10\%$.

Laser calorimetry was conducted at the University of Central Florida College of Optics and Photonics.^{20} The experiment consists of thermally equilibrating a suspended test sample in the calorimeter, then exposing it to a laser while measuring the temperature rise of the specimen, and then continuing to monitor sample temperature after halting the laser exposure. Prior to calorimetry, specimens were cleaned with acetone and a cotton bud, washed with methanol, and wiped with lens paper. An insulated box enclosed the experiment to thermally isolate the sample at ambient temperature ($295\pm 1\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{K}$) in the air. Laser light entered and exited the box through holes surrounded by insulating material to reduce convective heat transfer to the atmosphere.

Calorimetry samples were optically polished disks with a diameter of 10 mm and nominal thicknesses of 1.0, 2.0, and 3.0 mm. One edge was ground flat over a length of 2 to 3 mm for gluing a monofilament nylon fishing line to suspend the specimen in the path of the laser beam. The fishing line was attached to a thermally isolated frame mounted with 3 degrees of translational and rotational movement. A precision temperature measurement system by GEC Instruments (model S4TC) allowed the simultaneous reading of four temperatures with precision better than 10 mK over a broad range. Three 0.2-mm diameter thermocouples were glued to the edge of each sample. Care was taken to prevent scattered laser light from striking the temperature sensors. The sensitivity of the calorimeter was $\sim 10\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mK}$.

Lasers emitting at 1.064 *μ*m (Nd:YAG, $\sim 3\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{W}$ output), 2.94 *μ*m (Er:YAG, $\sim 100\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mW}$), $\sim 3.41\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mu \mathrm{m}$ (cascade diode laser from Maxion Technologies, College Park, MD, $\sim 15\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mW}$), and $\sim 4.85\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mu \mathrm{m}$ (cascade diode laser from Maxion, $\sim 90\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mW}$) were used to irradiate the samples in the calorimeter via beam-directing mirrors. Laser power was measured as it entered the calorimeter using two calibrated power meters to check for consistency. If a laser beam was too large in diameter, or was divergent, a lens was used to focus the beam to $<3\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mm}$. Power was measured after the lens. The 4.85-*μ*m laser radiated over the range of 4.79 to 4.94 *μ*m, with half of the intensity on either side of 4.85 *μ*m.

Finite element analysis was carried out at the Naval Air Warfare Center with NX 7.5 finite element analysis software (Seimens PLM Software, Munich, Germany). A rectangular window was divided into 8-node cubic elements, with 51,300 elements in one plane and five planes through the thickness of the window. A static load of 13.79 kPa ($2\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{lb}/{\mathrm{in.}}^{2}$) was applied to the upper surface of the window. The analysis identifies in-plane stresses on the upper and lower surfaces of the window and computes the maximum principal stress in each surface element. For a simply supported window, nodes along both lower long edges and one of the lower short edges were fixed in the vertical ($z$) direction, but were free to rotate and to move in the $x$ and $y$ directions. The row of nodes on the other short, lower edge of the window was constrained in all six degrees of freedom. The model therefore has one edge restrained and three simply supported. Tensile stress adjacent to the constrained edge on the tensile face was neglected in computing the probability of survival of the window, so as to approximate a fully simply supported window. The one condition that was modeled was a window with dimensions of $48.26\times 68.58\times \phantom{\rule{0ex}{0ex}}1.27\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{cm}$, with the length of each finite element being 0.254 cm. For other pressures and window thicknesses, stresses were scaled in proportion to $p/{b}^{2}$, where $p$ is the uniform pressure on the upper surface of the window and $b$ is the thickness of the window. For example, if the applied pressure was doubled, stresses were doubled. If the window thickness was doubled, stresses were decreased by a factor of 4.

## 3.

## Thermal and Elastic Properties of Spinel

## 3.1.

### Specific Heat Capacity

The specific heat capacity of spinel was fit by the equation

## (1)

$${c}_{p}[\mathrm{J}/(\mathrm{g}\xb7\mathrm{K})]={a}_{0}+{a}_{1}T+{a}_{2}{T}^{2}+{a}_{3}{T}^{3}\phantom{\rule{0ex}{0ex}}+{a}_{4}{T}^{4}(T=296\text{\hspace{0.17em}}\text{to}\text{\hspace{0.17em}}1273\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{K}),$$^{21}For sapphire at 295 K, ${c}_{p}=0.768\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{J}/(\mathrm{g}\xb7\mathrm{K})$.

^{12}

## 3.2.

### Thermal Expansion

Thermal expansion of spinel was fit to the polynomial

## (2)

$$\mathrm{\Delta}L/{L}_{0}={b}_{0}+{b}_{1}T+{b}_{2}{T}^{2}+{b}_{3}{T}^{3}(T=293\text{\hspace{0.17em}}\text{to}\text{\hspace{0.17em}}1273\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{K}),$$^{22}whose data are fit by a cubic polynomial with ${b}_{0}=-2.0110\times {10}^{-3}$, ${b}_{1}=5.8103\times {10}^{-6}$, ${b}_{2}=2.0489\times \phantom{\rule{0ex}{0ex}}{10}^{-9}$, and ${b}_{3}=\u20130.1339\times {10}^{-12}$ for $T=293$ to 1900 K. The differential coefficient of thermal expansion (CTE) derived from Eq. (2) is

## (3)

$$\mathrm{CTE}=\frac{1}{L}\frac{\mathrm{d}L}{\mathrm{d}T}=\frac{{b}_{1}+2{b}_{2}T+3{b}_{3}{T}^{2}}{1+{b}_{0}+{b}_{1}T+{b}_{2}{T}^{2}+{b}_{3}{T}^{3}}\phantom{\rule[-0.0ex]{2em}{0.0ex}}\phantom{\rule{0ex}{0ex}}(T=293\text{\hspace{0.17em}}\text{to}\text{\hspace{0.17em}}1273\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{K}).$$## 3.3.

### Density

The density of spinel ($3.576\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{g}/\mathrm{mL}$) at 293 K was calculated from the cubic unit cell lattice dimension of 0.8085 nm and eight formula masses ($142.266\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{g}/\mathrm{mol}$) per unit cell.^{21}^{,}^{23} At elevated temperature $T$, the density ${\rho}_{T}$ decreases from thermal expansion.

## 3.4.

### Thermal Conductivity

Thermal diffusivity ($D$) is a measure of the rate at which temperature diffuses through a substance when there is a spatial gradient of temperature. Thermal conductivity, $k$, is computed from the measured diffusivity by the equation

where $\rho $ is the density in Eq. (4), ${c}_{p}$ is the specific heat capacity in Eq. (1), and $D$ is the thermal diffusivity. Our data in Fig. 2 are fit by the equation where $A=443174$, $B=1.8692$, $c=7.3970$, and $T=296$ to 1273 K, and the residual standard deviation of the curve fit is $0.15\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{W}/(\mathrm{m}\xb7\mathrm{K})$. Figure 2 also shows previously reported values of the thermal conductivity of spinel.^{22}

^{,}

^{24}

^{,}

^{25}

## 3.5.

### Elastic Moduli and Poisson’s Ratio

Young’s modulus ($E$) and Poisson’s ratio ($\nu $) measured from vibrational frequencies of test specimens were used to compute the shear modulus ($G$) and bulk modulus ($K$). Results are $E=272\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{GPa}$, $\nu =0.27$, $G=106\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{GPa}$, and $K=\phantom{\rule{0ex}{0ex}}197\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{GPa}$ for $T=294\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{K}$. Our value of Young’s modulus can be compared with values of 277,^{26} 275,^{27} 266 to 270,^{28} 268,^{29} and 258^{30}^{,}^{31} GPa in previous reports.

## 4.

## Laser Calorimetry

## 4.1.

### Treatment of Laser Calorimetry Data

Laser calorimetry is recommended for measuring the low absorptance of optical components and window materials.^{32}33.34.^{–}^{35} In this technique, a small disk is irradiated with a laser while the temperature of the disk is measured by thermocouples attached to the edge. After irradiation, the temperature falls back to ambient in an exponential manner. The graph of temperature versus time in Fig. 3 can be used to extract the absorptance ($a$) of the material, which is the fraction of incident radiant power (${P}_{0}$) absorbed.

Advantages of laser calorimetry include the simplicity of the apparatus and ease of absolute calibration. Surface and bulk absorptance can be distinguished by studying samples with different thicknesses. If the surfaces are identical, then a plot of measured absorptance versus sample thickness extrapolated to zero thickness reveals the surface absorptance. The slope of the graph gives the bulk absorptance.

If the laser calorimetry test specimen were not losing heat during laser irradiation, and if the thermal conductivity of the sample is high enough that its temperature ($T$) is uniform, the rate at which its temperature would rise is equal to the power absorbed divided by the thermal mass ($\text{mass}\times \text{specific heat capacity}$) of the disk.

where $t$ is time and $m$ is the mass of the sample. The numerator on the right side of Eq. (8) is the power absorbed by the sample. The denominator is the energy required to raise the sample temperature by 1 K. In the absence of heating, the rate at which the sample cools is proportional to the difference in temperature between the sample and its surroundings ($\mathrm{\Delta}T$). where $\gamma $ is an empirical heat loss coefficient accounting for radiative and conductive heat losses.At any time during the calorimetry experiment, the rate of change of temperature is the sum of Eqs. (8) and (9).

## (10)

$$\text{Rate of change of temperature in laser calorimetry}=\frac{\mathrm{d}T}{\mathrm{d}t}\phantom{\rule{0ex}{0ex}}=\frac{a{P}_{0}}{m{c}_{p}}-\gamma \mathrm{\Delta}T.$$Integration of Eq. (10) gives the temperature rise of the sample versus time.

## (11)

$$\text{Duration irradiation:}\text{\hspace{0.17em}}\mathrm{\Delta}T=\frac{a{P}_{0}}{\gamma m{c}_{p}}[1-{e}^{-\gamma (t-{t}_{\text{start}})}].$$## (12)

$$\text{After irradiation:}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{\Delta}T=\frac{a{P}_{0}}{\gamma m{c}_{p}}[{e}^{-\gamma (t-{t}_{\text{stop}})}-{e}^{-\gamma (t-{t}_{\text{start}})}],$$Internal transmittance is the fraction of radiant power that has entered a specimen, which reaches the opposite side after traversing a path length $b$ through the specimen. The relation between internal transmittance and absorption coefficient ($\alpha $) is

Internal transmittance is independent of surface losses by reflection, absorption, and scatter.For perpendicular incidence of a laser beam on a sample with parallel surfaces, negligible optical scatter, and negligible surface absorption, the absorption coefficient is related to absorptance by Eq. (14),^{35} which is derived in the appendix.

When the product of the absorption coefficient and the sample thickness is small ($\alpha b\ll 1$), a first-order expansion of Eq. (14) gives

If there is surface absorption, a term accounting for both surfaces can be added to Eq. (17).

In the present work, surface absorption was negligible relative to the precision of the measurement of bulk absorption. Equation (17) is typically accurate enough for laser calorimetry measurements of weakly absorbing specimens. However, Eq. (14) was required for sapphire and spinel at 4.85 *μ*m, where there is significant absorption.

## 4.2.

### Laser Calorimetry Results

Figures 4 and 5 show the absorptance of sapphire and spinel at 4.85 *μ*m and the absorptance of spinel at 3.41, 2.94, and 1.064 *μ*m. Error bars on the data points are estimated standard deviations from several measurements. Straight lines through the origin were drawn through the spinel data for 3.41, 2.94, and 1.064 *μ*m in Fig. 5. The slope is the absorption coefficient in Eq. (17).

The absorptance of spinel and sapphire at 4.85 *μ*m is too high to use Eq. (17), which produces values of absorption coefficient that are 3 to 9% low. Therefore, absorptance was fit to Eq. (14), with the single-surface Fresnel reflectance and transmittance calculated from the literature refractive index^{36}37.^{–}^{38} with Eqs. (15) and (16). The only adjustable parameter in the curve fit is absorptance. Curve fitting was done with Microsoft Excel Solver® as described by Harris.^{39} The standard uncertainty in absorptance from the curve fit was computed with Billo’s Excel macro Solvstat.^{40} The 95% confidence interval for absorptance was obtained by multiplying the standard uncertainty by Student’s $t$ for $n\text{-}p$ degrees of freedom, where $n=3$ is the number of data points on each curve and $p=1$ is the number of adjustable parameters.

Resulting absorption coefficients are listed in Table 1, where comparison is made with previously reported laser calorimetry measurements at several shorter wavelengths.^{41}42.43.^{–}^{44} The absorption coefficients of sapphire and spinel at 4.85 *μ*m are thought to represent intrinsic properties of the materials. Absorption coefficients at 1.064, 2.94, and 3.41 *μ*m are thought to arise from extrinsic factors such as OH impurities. It is noteworthy that the spinel we measured has $>50$ times as much absorption as sapphire at 1.064 *μ*m, and substantial absorption in the OH stretching region at 2.94 and 3.41 *μ*m.

## Table 1

Absorption coefficients measured by laser calorimetry.

Material | Absorption coefficient, α (cm−1) | |||
---|---|---|---|---|

Present work | ||||

1.064 μm | 2.94 μm | 3.41 μm | 4.85 μm | |

c-Plane sapphire | 0.0004±0.0001 | * | * | 0.58±0.17 |

a-Plane sapphire | 0.0003±0.0001 | * | * | 0.42±0.13 |

Spinel | 0.021±0.001 | 0.08±0.01 | 0.05±0.04 | 0.34±0.09 |

Literature data | ||||

1.3 μm | 2.7 μm | 3.8 μm | ||

c-Plane sapphire | 0.0008,a,b 0.00005c | 0.00127,a 0.003d | 0.024,b,d 0.0285a,b | |

Spinel | 0.00655a | 0.0514a | 0.0404a |

Note: Uncertainty is 95% confidence interval. The absorption coefficient for sapphire at 2.94 and 3.41 μm was too small to measure with the laser power available.

## a

Reference 41.

## b

Reference 43.

## c

Reference 44.

## d

Reference 42.

Figure 6 superimposes our results from laser calorimetry with previous optical transmittance measurements made with laser diodes at discrete wavelengths or by FTIR spectrometry.^{5}^{,}^{6}^{,}^{45} Our calorimetry results are represented by three large open symbols at 4.85 *μ*m. Dashed lines are from OPTIMATR. Solid circles are reported points for $c$-plane sapphire. Solid squares are reported points for $a$-plane sapphire. Solid triangles are reported points for spinel. Two superimposed open symbols near the origin at 3.8 *μ*m are laser calorimetry points for $c$-plane sapphire from the literature listed in Table 1.^{42}^{,}^{43}

We obtained the solid curve in Fig. 6 by measuring the transmittance of spinel with thicknesses of 0.598 and 1.276 cm. If surface absorption and surface scatter from the two specimens are the same, the bulk extinction coefficient ($\beta ={\alpha}_{\text{absorption}}+{\alpha}_{\text{bulk scatter}}$), arising from bulk absorption plus bulk scatter, is given by

## (19)

$$\beta ={\alpha}_{\text{absorption}}+{\alpha}_{\text{bulk scatter}}=\frac{\mathrm{ln}({T}_{1}/{T}_{2})}{{b}_{2}-{b}_{1}},$$*μ*m is 0.346% for a 0.598-cm thick specimen and 0.418% for a 1.276-cm thick specimen. (The total integrated scatter from the two samples of spinel did not vary with wavelength in the interval of 3.5 to 4.4

*μ*m.) If surface scatter of both specimens is equal, the difference in bulk scatter is $0.418-0.346=0.072\%$ for a thickness difference of $1.276-\phantom{\rule{0ex}{0ex}}0.598=0.678\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{cm}$. The bulk scatter coefficient is therefore ${\alpha}_{\text{bulk scatter}}=0.072\%/0.678\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{cm}=0.0011\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{cm}}^{-1}$. This negligible scatter coefficient was subtracted from the extinction coefficient $\beta $ in Eq. (19) to obtain the absorption coefficient ${\alpha}_{\text{absorption}}$ shown by the solid line in Fig. 6.

## 5.

## Revised Spectral Absorption Coefficient of Spinel

In Fig. 6, experimental data points for the absorption coefficient of spinel (solid triangles) and our measurement from two thicknesses of spinel [solid line computed with Eq. (19)] are superimposed, but do not agree with the dotted line for spinel computed with OPTIMATR. By contrast, experimental data points for $c$-plane sapphire (solid circles) and $a$-plane sapphire (solid squares) are in good agreement with the dashed curves computed by OPTIMATR. Our laser calorimetry measurement for $a$-plane sapphire at 4.85 *μ*m (open square) lies on the OPTIMATR curve. Our laser calorimetry measurement for $c$-plane sapphire (open circle) lies 12% below the OPTIMATR curve, but within the estimated uncertainty of the calorimetry measurement. For spinel, our laser calorimetry measurement (open triangle) lies between the dashed and solid curves for spinel, but within experimental uncertainty of both. In summary, experimental data points are in good agreement with the OPTIMATR curves for sapphire, but experimental points for spinel do not lie on the OPTIMATR curve in the region from 4.2 to 5.2 *μ*m. OPTIMATR underestimates the absorptance of spinel in the 4.2 to 5.2 *μ*m region.

We therefore made spectroscopic measurements of the absorption coefficient of additional samples of spinel to see if the solid curve in Fig. 6 is reproducible. Figure 7 shows that spectroscopic data from new samples of spinel from TA&T and from MER Corp. lie near the solid curve from Fig. 6 and not near the OPTIMATR curve. Table 2 lists the average absorption coefficient of spinel from measurements of three pairs of thick and thin specimens. We found other spinel samples that gave higher extinction coefficients than those shown in Table 2. We assume that these materials had higher optical scatter and deliberately used samples that gave the lowest extinction coefficients so as not to overestimate the absorption of spinel.

## Table 2

Mean absorption coefficient of spinel from transmittance of three pairs of thin and thick specimens.

Wavenumber (cm−1) | Absorption coefficient (α, cm−1) | Wavenumber (cm−1) | Absorption coefficient (α, cm−1) |
---|---|---|---|

2600 | 0.011 | 2225 | 0.149 |

2575 | 0.014 | 2200 | 0.174 |

2550 | 0.016 | 2175 | 0.202 |

2525 | 0.018 | 2150 | 0.233 |

2500 | 0.021 | 2125 | 0.268 |

2475 | 0.025 | 2100 | 0.307 |

2450 | 0.029 | 2075 | 0.352 |

2425 | 0.035 | 2050 | 0.402 |

2400 | 0.042 | 2025 | 0.460 |

2375 | 0.051 | 2000 | 0.527 |

2350 | 0.062 | 1975 | 0.603 |

2325 | 0.075 | 1950 | 0.692 |

2300 | 0.090 | 1925 | 0.795 |

2275 | 0.107 | 1900 | 0.915 |

2250 | 0.127 | 1875 | 1.052 |

Note: A polynomial fitting absorption coefficient as a function of wave number (ν˜, cm−1) from 2600 to 1875 cm−1 is α (cm−1)=5.51852179 × 10−18ν˜6 − 8.88520838 × 10−14ν˜5 + 5.85969325 ×10−10ν˜4−2.03475841×10−6ν˜3+3.93709258×10−3ν˜2−4.03652837ν˜+1717.85263. Full nine-digit coefficients are required to reproduce α with three decimal precision.

The pair of spinel samples that gave the lowest absorption coefficient in the 4 to 5 *μ*m region was used to calculate the extinction coefficient at shorter wavelengths. Figure 8 shows two weak absorption features near 1.8 and 2.8 *μ*m, which are evident in all samples of spinel that are thick enough to see the absorption. Note that the absorption coefficients from laser calorimetry for spinel in Table 1 are roughly twice as great as the values in Fig. 8. The values in Fig. 8 could therefore be an underestimate.

For comparison of the absorption of sapphire and spinel in Sec. 7 of this paper, we calculated the absorptance ($a$) of a 1-cm-thick slab of each material from the absorption coefficient by using Eq. (14), where $\tau =1-R$ is the single-surface transmittance. We then computed the average absorptance over the wave number range $\tilde{\nu}=2000\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{to}\text{\hspace{0.17em}}\phantom{\rule{0ex}{0ex}}2632\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{cm}}^{-1}$ (3.8 to 5.0 *μ*m) with Eq. (20).

## (20)

$$\text{Average absorptance}=\frac{1}{2632-2000}{\int}_{2000}^{2632}{a}_{\tilde{v}}\mathrm{d}\tilde{v}.$$Results are shown in Table 3, which also includes the average absorptance for the range 3.0 to 5.0 *μ*m using the absorption coefficient for spinel from Fig. 8. The average absorptance of spinel (0.134) in the 3.8 to 5.0 *μ*m region is 92% as great as that of $a$-plane sapphire (0.145) and 71% as great as that of $c$-plane sapphire (0.189). OPTIMATR predicts significantly lower average absorptance for spinel (0.094). The average absorptance of spinel (0.073) in the wider range of 3.0 to 5.0 *μ*m is slightly greater than that of $a$-plane sapphire (0.071) because spinel has absorption in the 3.0 to 3.8 *μ*m region (Fig. 8) and sapphire does not. OPTIMATR does not model any spinel absorption below 3.8 *μ*m.

## Table 3

Average absorptance of spinel and sapphire.

Data source | Wavelength range | Wave number range | Average absorptance for 1 cm thickness, Eq. (20) | ||
---|---|---|---|---|---|

c-Plane sapphire | a-Plane sapphire | Spinel | |||

OPTIMATR (295 K) | 3.8 to 5.0 μm | 2632 to 2000 | 0.189 | 0.145 | 0.094 |

Present work | 3.8 to 5.0 μm | 2632 to 2000 | a | a | 0.134 |

OPTIMATR (295 K) | 3.0 to 5.0 μm | 3333 to 2000 | 0.092 | 0.071 | 0.045 |

Present work | 3.0 to 5.0 μm | 3333 to 2000 | a | a | 0.073 |

## a

Spectral transmittance of sapphire was not measured in the current work. Laser calorimetry and measurements of sapphire conform to OPTIMATR.

## 6.

## Refractive Index of Spinel

Published^{2}^{,}^{46}^{,}^{47} and unpublished measurements of the refractive index ($n$) of spinel at room temperature from 0.35 to 5.5 *μ*m were summarized in 1995 by the Sellmeier equation. ^{38}

## (21)

$${n}^{2}-1=\frac{1.8938{\lambda}^{2}}{{\lambda}^{2}-{(0.09942)}^{2}}+\frac{3.0755{\lambda}^{2}}{{\lambda}^{2}-{(15.826)}^{2}},$$*μ*m. Spinel samples included natural minerals and synthetic material from the 1980s and earlier.

We measured the refractive index of a single prism of commercial polycrystalline spinel manufactured in 2010. Measurements near 295 K were made at 60 wavelengths between 0.4 and 5.0 *μ*m. An additional 93 measurements were made near 350, 400, and 450 K at selected wavelengths. The entire data set was fit to the following Sellmeier equation by the Levenberg-Marquart algorithm with the Excel macro Optimiz.^{48}

## (22)

$${n}^{2}-1=\frac{[A+B({T}^{2}-{T}_{0}^{2})]{\lambda}^{2}}{{\lambda}^{2}-{[C+D({T}^{2}-{T}_{0}^{2})]}^{2}}+\frac{E{\lambda}^{2}}{{\lambda}^{2}-{F}^{2}},$$*μ*m, $T$ is the temperature in the range of 295 to 450 K, ${T}_{0}=295.15\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{K}$, and the least-squares parameters are $A=1.893592$, $B=5.904\times {10}^{-8}$, $C=0.097911$, $D=4.437\times {10}^{-9}$, $E=\phantom{\rule{0ex}{0ex}}4.617875$, and $F=19.0165$. The root-mean-square error in fitting the refractive index is 0.00014. The root-mean-square difference between Eqs. (21) and (22) at 295 K is 0.0009, with the value of $n$ from Eq. (22) being higher than the value from Eq. (21). A 2 mol% increase in the ${\mathrm{Al}}_{2}{\mathrm{O}}_{3}:\mathrm{MgO}$ ratio would increase the refractive index by 0.0009.

The refractive index increases linearly in the range of 25 to 175°C, with a slope $\mathrm{d}n/\mathrm{d}T\approx 1.4\times {10}^{-5}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{K}}^{-1}$. For comparison, Lang and Duncan^{49} found $n=1.6837+1.30\phantom{\rule{0ex}{0ex}}(\pm 0.04)\times {10}^{-5}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{T}-3.6\times {10}^{-9}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{T}}^{2}$ ($\mathrm{T}=\xb0\mathrm{C}$) at 0.633 *μ*m in the range of 25 to 525°C. The average slope of this function from 25 to 175°C is $1.2\times {10}^{-5}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{K}}^{-1}$. Vadem et al.^{47} reported $\mathrm{d}n/\mathrm{d}T=0.90\times {10}^{-5}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{K}}^{-1}$ at 0.589 *μ*m for the range of 25 to 700°C.

## 7.

## Trading Optical Absorption and Mechanical Strength

The question that motivated our study of sapphire and spinel is which of these infrared window materials can provide the least midwave infrared absorption while meeting the mechanical requirements for a window? From Fig. 1, the long held conventional wisdom has been that increased infrared transmittance in the 4 to 5 *μ*m range is provided in the order $c\text{-plane sapphire}<a\text{-plane sapphire}<\text{spinel}$. However, careful measurements in Fig. 6 show that the absorption coefficients of spinel and $a$-plane sapphire are almost equal in the range of 3.8 to 4.8 *μ*m and diverge only at longer wavelengths. Sapphire is mechanically stronger than spinel, so a thin sapphire window could possibly be used instead of a thicker spinel window, resulting in lower absorption by the thin sapphire than by the thicker spinel. Our goal is to make quantitative comparisons to select the material that has the least infrared absorption, while meeting mechanical requirements.

One way to design a window plus its frame employs finite element analysis to calculate the stress distribution in the window under a given load and to use Weibull statistics to estimate the probability of survival of the window. If the desired probability of survival cannot be achieved with a practical thickness of material, a proof test on a finished window could be used to demonstrate that a particular window survives required loads. Another approach to window design is to measure slow crack growth rate of the window material (ASTM C1368) and then design the window to last for a chosen lifetime under operational stresses.^{29}^{,}^{50}

Consider a simply supported rectangular window subjected to a uniform pressure on one face. Maximum tensile stress (${\sigma}_{\mathrm{max}}$) on the tensile face of the window is proportional to the applied pressure ($p$) and decreases as the inverse square of window thickness ($b$).^{51}

To compare strengths of materials, we employ the characteristic strength (${s}_{c}$) defined by Klein as the effective mean strength of a $1\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{cm}}^{2}$ area under tension in a mechanical test.^{52}53.^{–}^{54} The strength of an optical ceramic is a highly variable, extrinsic property that depends on flaws such as scratches, subsurface damage, weak grain boundaries, large grains, inclusions, pore nests, and voids. The distribution of flaw sizes determines the distribution of strength, which is characterized by the Weibull modulus, $m$. The smaller the Weibull modulus, the broader is the distribution of strength.

For window design consideration, we choose the parameters ${s}_{c}$ and $m$ in columns 2 and 3 of Table 4, which approximately represent particular lots of real specimens. Three sets of Weibull parameters for spinel in Table 4 encompass the behavior of large-grain, commercial material for which we have data. Weibull parameters for sapphire are roughly representative of several widely varying reports.^{52}^{,}^{55}56.^{–}^{57} Sapphire and spinel could have widely different characteristics from those listed in the table. For example, experimental coupons of spinel with a grain size of 2 *μ*m had a characteristic strength of $\sim 425\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{MPa}$ and a Weibull modulus of 6.9.^{58} In our experience, Table 4 approximately represents the range of strength and Weibull modulus of high optical quality spinel that can be currently purchased in the form of large windows with lateral dimensions $>30\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{cm}$. The characteristic strength of sapphire in Table 4 might appear high in comparison to strengths listed in some manufacturer’s bulletins. The reason is that the characteristic strength ${s}_{c}$ applies to $1\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{cm}}^{2}$ area under load. When scaled for area with a Weibull modulus of 4, the strength of sapphire would look more familiar. For example, for an area of $20\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{cm}}^{2}$ instead of $1\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{cm}}^{2}$, the strength of $c$-plane sapphire is predicted to be $(900\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{MPa}){(1\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{cm}}^{2}/20\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{cm}}^{2})}^{1/4}=426\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{MPa}$

## Table 4

Applied pressure and window thickness that give Ps=0.999 for spinel and sapphire windows with different strength distributions.

Material | sc (MPa) | m | Applied pressure p (kPa) | Thickness b (cm) | Max stress in window (MPa) | Mean absorptance 3.8 to 5.0 μm | Relative absorptance 3.8 to 5.0 μm |
---|---|---|---|---|---|---|---|

Spinel 1 | 100 | 8 | 39.5 | 1.27 | 21.0 | 0.164 | 1.00 |

Spinel 2 | 150 | 6 | 39.5 | 1.38 | 17.7 | 0.175 | 1.07 |

Spinel 3 | 200 | 6 | 39.5 | 1.20 | 23.6 | 0.156 | 0.95 |

c-Sapphire | 900 | 4 | 39.5 | 1.02 | 32.4 | 0.192 | 1.17 |

a-Sapphire | 600 | 4 | 39.5 | 1.25 | 21.6 | 0.174 | 1.06 |

Note: Simply supported window (48.26×68.58×b cm thick) with pressure p applied to one surface.

For each element of area on the tensile surface of a window in a finite element model, the probability of survival (${P}_{s}$) is given by the Weibull distribution, which we choose to write as suggested by Klein.^{52}53.^{–}^{54}

## (24)

$${P}_{s}={e}^{-(A/{A}_{0}){\left[\mathrm{\Gamma}\right(1+\frac{1}{m}\left)\right]}^{m}{\left(\frac{\sigma}{{s}_{c}}\right)}^{m}},$$Finite element analysis was conducted for a simply supported rectangular window with dimensions $48.26\times 68.58\times b\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{cm}$ with a uniform pressure $p$ applied to one face. The thickness ($b$) of the window was varied to produce a probability of survival of 0.999 for each set of Weibull parameters. In Table 4, for spinel 1 with ${s}_{c}=100\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{MPa}$ and $m=8$, a thickness of $b=1.27\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{cm}$ and applied pressure $p=39.5\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{kPa}$ give a probability of survival of 0.999. The maximum stress at the center of the tensile face is 21.0 MPa. For $c$-plane sapphire with ${s}_{c}=900\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{MPa}$ and $m=4$, the same applied pressure of 39.5 kPa yields a probability of survival of 0.999 if the thickness is 1.02 cm. The stress at the center of the tensile face is 32.4 MPa. Three cases for spinel and two cases for sapphire are listed in Table 4.

From the thickness of the window required to produce a survival probability of 0.999, the mean absorptance for the wavelength range of 3.8 to 5.0 *μ*m was computed by solving Eq. (14) for absorptance at each wavenumber and then integrating Eq. (20). Absorptance is shown in the next-to-last column of Table 4 and relative absorptance is shown in the last column. We find that $a$-plane sapphire has equivalent absorptance to spinel 2 whose characteristic strength is 150 MPa and whose Weibull modulus is 6. Spinel 1 with a higher Weibull modulus and spinel 3 with a higher strength have slightly less absorptance than $a$-plane sapphire. The absorptance of $c$-plane sapphire is 9 to 23% greater than that of the three types of spinel.

The conclusion of this analysis is that spinel and $a$-plane sapphire have similar mean absorptance in the 3.8 to 5.0 *μ*m region. When designed to bear the same load, either material could be better than the other, depending on their Weibull parameters. This conclusion contradicts the common belief that a spinel window will have less midwave infrared absorptance than a sapphire window. Spinel does have less absorptance than $c$-plane sapphire, but does not have a clear advantage over $a$-plane sapphire.

## 8.

## Conclusion

Careful measurements of infrared absorption of high-quality, commercial, polycrystalline $\mathrm{Mg}{\mathrm{Al}}_{2}{\mathrm{O}}_{4}$ spinel show that the absorptance is 43% greater in the range of 3.8 to 5.0 *μ*m than the value predicted by the computer code OPTIMATR®, which has been used for window and dome design for more than 20 years. The average absorptance of spinel in the wavelength range of 3.8 to 5.0 *μ*m is 92% as great as that of $a$-plane sapphire and 71% as great as that of $c$-plane sapphire. Spinel also exhibits weak absorptance in the range of 1.4 to 3.4 *μ*m. As a result, the average absorptance of spinel in the range of 3.0 to 5.0 *μ*m is 103% as great as that of $a$-plane sapphire and 79% as great as that of $c$-plane sapphire. The spinel that we examined has $>50$ times as much optical absorption as sapphire at 1.064 *μ*m. Specific heat capacity, thermal expansion, thermal conductivity, elastic constants, and refractive index of spinel are also reported. We examined the relative infrared absorptance of sapphire and spinel windows designed to provide the same probability of survival under an applied load. The conclusion is that spinel and $a$-plane sapphire have similar mean absorptance in the 3.8 to 5.0 *μ*m region. The absorptance of $c$-plane sapphire is greater than that of spinel.

## Acknowledgments

We thank Larry Fehrenbacher of TA&T and Juan Sepulveda of MER Corp. for providing spinel and sharing data on spinel. We are grateful to GT Advanced Technologies (formerly Crystal Systems) for providing sapphire. Raytheon thermal conductivity data were kindly provided by Howard Poisl.

## Appendices

## Appendix:

### Derivation of Eq. (14) Relating Absorption Coefficient to Absorptance

Absorptance can occur in the bulk material and at both surfaces. Following the derivation of Orchard et al.,^{35} we assume that surface absorption is negligible compared with bulk absorption. The derived formula therefore places an upper limit on bulk absorption.

In Fig. 9, radiant power ${P}_{0}$ is incident on the first surface of the optical disk at the left. If the single-surface transmittance is $\tau $, and if there is negligible surface absorption and scatter, then the power that enters the first surface is $\tau {P}_{0}$. When the forward ray reaches the right-hand surface, some is transmitted to the right and some is reflected back to the left. When the reflected beam reaches the left-hand surface, some is transmitted to the left and some is reflected back to the right. There are an infinite number of progressively smaller internal reflections. The total radiant power traveling in the forward direction (to the right) inside the sample is designated ${P}_{F}$ in Fig. 9. This power is the sum of the power that entered the sample and the infinite number of reflections from the left-hand surface. Similarly, the total radiant power traveling in the reverse direction (to the left) inside the sample is designated ${P}_{R}$.

At the output face, the total transmitted power is ${P}_{T}$, so the total power incident on the output face from inside the sample must be ${P}_{T}/\tau $ because the transmittance of the output face is $\tau $. The total power propagating from left to right inside the sample is ${P}_{F}$. It is attenuated over a path length $b$ by a material with a bulk absorption coefficient $\alpha $. We say^{35} that power ${P}_{T}/\tau $ reaching the output face is a result of power ${P}_{F}$ being attenuated by a factor ${e}^{\u2013\alpha b}$.

## (25)

$${P}_{T}/\tau ={P}_{F}{e}^{-\alpha b}\phantom{\rule[-0.0ex]{1.0em}{0.0ex}}\mathrm{or}\phantom{\rule[-0.0ex]{1.0em}{0.0ex}}{P}_{F}=({P}_{T}/\tau ){e}^{\alpha b}.$$If the single-surface reflection coefficient is $R$, and total power ${P}_{T}/\tau $ reaches the output face, then the total reflected power propagating to the left inside the sample is

The total power absorbed inside the sample is a result of powers ${P}_{F}$ and ${P}_{R}$ each being attenuated over path length $b$.

## (27)

$${P}_{\text{absorbed}}=\underset{\text{radiant power in disk}}{\underbrace{({P}_{F}+{P}_{R})}}\times \underset{\text{fraction of power absorbed}}{\underbrace{(1-{e}^{-\alpha b})}}=a{P}_{0},$$## (28)

$$a{P}_{0}=({P}_{F}+{P}_{R})(1-{e}^{-\alpha b})=(\frac{{P}_{T}}{\tau}{e}^{\alpha b}+\frac{{P}_{T}}{\tau}R)(1-{e}^{-\alpha b})\phantom{\rule{0ex}{0ex}}=\frac{{P}_{T}}{\tau}({e}^{\alpha b}+R)(1-{e}^{-\alpha b}).$$The total transmitted power (${P}_{T}$ in Fig. 9) results from the incident power ${P}_{0}$ being transmitted through two surfaces and suffering an infinite number of partial reflections between those two surfaces. The following equation accounts for all losses:

## (29)

$${P}_{T}=\underset{\genfrac{}{}{0ex}{}{\text{loss from tranmission}}{\text{through two surfaces}}}{\underbrace{{P}_{0}{\tau}^{2}{e}^{-\alpha b}}}\underset{\genfrac{}{}{0ex}{}{\text{one term for each}}{\text{round-trip reflection}}}{\underbrace{(1+{R}^{2}{e}^{-2\alpha b}+{R}^{4}{e}^{-4\alpha b}+{R}^{6}{e}^{-6\alpha b}+\cdots )}}.$$The parentheses contain an infinite geometric series of the form $\sum _{n=0}^{n=\infty}{x}^{n}$, where $x={R}^{2}{e}^{-2\alpha b}$. For $|x|<1$, the sum converges to $\sum {x}^{n}=1/(1-x)$, so

## (30)

$$1+{R}^{2}{e}^{-2\alpha b}+{R}^{4}{e}^{-4\alpha b}+{R}^{6}{e}^{-6\alpha b}+\cdots =\frac{1}{1-{R}^{2}{e}^{-2\alpha b}}\phantom{\rule{0ex}{0ex}}=\frac{1}{(1-R{e}^{-\alpha b})(1+R{e}^{-\alpha b})}.$$Substituting the expression at the right of Eq. (30) for the infinite series into Eq. (29) gives

Finally, we substitute the expression for ${P}_{T}$ from Eq. (31) into Eq. (28).

## (32)

$$a{P}_{0}=\frac{{P}_{T}}{\tau}({e}^{\alpha b}+R)(1-{e}^{-\alpha b})\phantom{\rule{0ex}{0ex}}=\frac{{P}_{0}{\tau}^{2}{e}^{-\alpha b}}{\tau (1-R{e}^{-\alpha b})(1+R{e}^{-\alpha b})}({e}^{\alpha b}+R)(1-{e}^{\alpha b})\phantom{\rule{0ex}{0ex}}a{P}_{0}=\frac{{P}_{0}{\tau}^{2}{e}^{-\alpha b}({e}^{\alpha b}+R)(1-{e}^{-\alpha b})}{\tau (1-R{e}^{-\alpha b})(1+R{e}^{-\alpha b})}\phantom{\rule{0ex}{0ex}}=\frac{{P}_{0}\tau (1+R{e}^{-\alpha b})(1-{e}^{-\alpha b})}{(1-R{e}^{-\alpha b})(1+R{e}^{-\alpha b})}\phantom{\rule{0ex}{0ex}}a=\frac{\tau (1-{e}^{-\alpha b})}{1-R{e}^{-\alpha b}}.$$Solving Eq. (32) for ${e}^{-\alpha b}$ gives

which is Eq. (14).Equation (14) relates the absorption coefficient ($\alpha $) to the measured absorptance ($a$). The single-surface Fresnel reflectance ($R$) is given by Eq. (15) and the Fresnel transmittance ($\tau $) is given by Eq. (16). If we know the refractive index of the material at the laser wavelength, Eq. (14) allows us to calculate the absorption coefficient from the measured absorptance.

## References

## Biography

**Daniel C. Harris** is a senior scientist and esteemed fellow at the Naval Air Systems Command, China Lake, California, where he is responsible for infrared window materials research and development programs. He has degrees in chemistry from Massachusetts Institute of Technology in 1968 and California Institute of Technology in 1973. He is the author of *Materials for Infrared Windows and Domes* and a series of widely used textbooks in analytical chemistry.

**Robert Seaver** received a BS in physics from the University of California, Los Angeles, in 2005. Since 2006, he has worked at the Naval Air Warfare Center Weapons Division, China Lake. His efforts focus on infrared optical metrology, including measurement of reflectance, scatter, absorption, and laser damage threshold.

**Giorgio Turri** graduated in 2001 at the Universita’ degli Studi di Milano and obtained his PhD in 2001 from the Politecnico di Milano. He began his research in experimental atomic and molecular physics, working at the Elettra Synchrotron in Trieste, Italy, and at the Advanced Light Source of the Lawrence Berkeley National Laboratory in California. His current interests are in optics and photonics, including laser calorimetry and characterization of lasing elements for solid-state lasers. He focuses on spectroscopic properties of rare-earth-doped crystals and glasses.

**Michael Bass** received a BS in physics from Carnegie-Mellon University and a PhD in physics from the University of Michigan in 1964. He has been active in quantum electronics research since his thesis work on optical rectification. He was among the first to demonstrate homogeneous broadening in dye lasers. He was responsible for the use of the statistical approach to understand laser-induced damage leading to the recognition of electron avalanche breakdown as the intrinsic cause of this phenomenon. He is the coinventor of the YAlO_{3} laser host and the use of fiber optics and lasers for the treatment of internal bleeding and tumors. He pioneered the use of laser calorimetry for the first absolute measurements of two-photon absorption. His current research centers on the design of improved solid-state lasers, materials with potential for use in upconversion lasers and displays, gain-guided index anti-guided fiber lasers, beam combining using volume Bragg gratings, and novel ceramic crystal optical materials. He is a fellow of OSA and a life fellow of IEEE.

*μ*m wavelengths and comparison with sapphire," Optical Engineering 52(8), 087113 (26 August 2013). https://doi.org/10.1117/1.OE.52.8.087113 . Submission: