We review the physics underlying thermal infrared emission, focusing on the processes that generate polarization. We also discuss a number of features of the theory that are not widely discussed in the literature, including the Straubel theorem for absorbing media, the mixed Poynting vector, and the emission of nonisothermal media. Using a model for emission from nonisothermal media, we present a quantitative discussion of the conditions under which Kirchhoff’s law is valid. Using this model, we show how it is possible for many blackbody simulators and infrared camera shutters to produce emission very close to true blackbodies. |
1.IntroductionAlthough research in infrared polarization has been pursued at least since Millikan’s studies on emission polarization in 1895,1,2 infrared polarimetry has never reached widespread use. As a result, there are many polarimetric features specific to infrared wavelengths that are not widely known, despite their fascinating physics. This review aims to survey the physics underlying infrared polarimetry and to explore a few of the less common features of the theory. In the infrared polarization literature, writers often refer to the thermal emission of surfaces, so it is easy to slip into thinking of infrared light emission as a surface phenomenon rather than a volumetric one. A search through the literature encounters “surface emissivity” as a common phrase, and Planck’s blackbody law, after all, is generally expressed in terms of surfaces and surface temperatures, rather than volumes and volumetric temperature distributions. Nevertheless, emission is a volumetric property, much like absorption but played out in reverse. For opaque materials, in which absorption completely attenuates incident light after passing only microns deep inside the material, this volumetric behavior is not obvious. For semitransparent (diathermous) materials, however, understanding the volumetric behavior is essential to accurately modeling their emission properties.3,4 In Sec. 2, we model the thermal light emission that constantly occurs inside every material and show that the surface emission that we see is simply the portion of that internal light that reaches the surface, transmits, and refracts into the air.3,5 This model holds even for highly reflective materials such as metals, and we find that the apparent skin-like quality of emission derives from the fact that so many materials in the infrared have strong absorption, so only the layers close to the surface are able to participate in generating the emitted light that we see. In Sec. 6, we also show that emitted light is polarized by the same Fresnel reflection and transmission effects that generate polarization from smooth visible-light surfaces. Because emission polarization is actually antagonistic to reflection polarization, the combination of emitted and reflected light often causes most of the polarimetric signature to cancel, producing the characteristically low degree of polarization values that are typically seen in this regime.6 In outdoor environments, this cancellation can become small for imaging quite hot surfaces and blue sky (very cold) backgrounds. The advantage of going carefully through the steps of deriving emission polarization from first principles shows that Kirchhoff’s law, which is often presented as an empirical formula, can be derived from basic physics. It also allows us to quantify the range inside which Kirchhoff’s law is valid (Sec. 7) and to demonstrate that, although thermal emission may be a totally incoherent process, interference effects still play a role in determining emission properties when thin films (with thicknesses on the order of a wavelength) are present or when using high-resolution spatial/spectral detection. In addition to improving understanding of the theory, some motivations for pursuing this work are to resolve some questions that are not easily clarified by a casual read through the literature. These questions include the following:
In the discussion below, we survey the theoretical tools needed to resolve these questions, though the discussion of anodize layers, and many of the details on paint layers are postponed until Part II of the review. To the extent possible, we aim to provide quantitative answers to these questions and bring together into one place the various work scattered through the literature on these topics. Part II of the review uses these tools to guide analysis of infrared imaging polarization measurements. 2.Calculating Emitted Surface Radiance with the Emission Density,The first widely known attempt to apply the tools of radiance transfer to the problem of volumetric emission from solids is that of McMahon in 1950.3 This work was later generalized by Gardon in 1956 to be applied to non-normal angles of emission.4 Despite its beautiful physics, the resulting model is not well known in the infrared imaging community even though it is widely used in the heat-transfer literature. In the discussion below, we modify the conventional derivation by centering the calculations around radiance rather than optical power. This has the advantage of simplifying expressions while also modernizing the units, so that they are more familiar to the modern reader. The calculation of the emissivity of an object can be quite involved if presented all at once, so we start with a simple setup and add complexity to the model in stages. At the first stage of the calculation, we aim to calculate the normal-angle emittance of a material with known refractive index and absorption index , i.e., a complex refractive index of . The electric field of a plane wave propagating along is then written as or where is the wavelength in vacuum and is the intensity absorption coefficient.In the discussion below, most variables will have a wavelength dependence. Writing an explicit wavelength dependence, such as , for each variable quickly becomes tedious and makes the equations more difficult to read. As a result, we generally drop the wavelength dependence from intermediate results and include the explicit dependence when it helps to clarify context. In addition, rather than “emitting medium,” we generally refer to the region from which we are calculating the emission as “the solid.” This simplifies the nomenclature, but readers should note that it can refer to liquids and gases as well, without difficulty. We start by postulating that the medium is uniform in temperature and composition and has negligible scattering and that each point within the solid can be characterized as emitting isotropic thermal radiation. For highly opaque materials, radiation emitted from within the volume is quickly absorbed by the surrounding matter and does not reach the surface. For semitransparent materials, however, this volumetric emission can reach the surface, so calculating the light emitted by the entire solid requires integrating over the object’s volume. The setup for this calculation is illustrated in Fig. 1. Figure 1 shows the radiation emitted from an element of volume in the interior of the solid that reaches a small area element of the surface, as originally constructed by Gardon.4 From this geometry, Gardon develops an integral over and the differential solid angle , taking advantage of the conical symmetry of the problem to define as a differential sector of a conical annulus. In this way, any region of the solid lying inside the region bounded by the blue lines [with differential angles and , as indicated in Fig. 1(b)], produces radiation propagating within the same element of the solid angle. The development of the integral is shown in Sec. 10, the end result of which is a simple one-dimensional (1D) integral along the ray between and . The quantity is the material’s spectral emission density, at a given temperature , expressed in units of —spectral radiance per unit length. Thus, we can expect that integrating over a length will lead to the familiar spectral radiance given by the Planck blackbody equation, viz., for Planck constant , Boltzmann constant , and the refractive index of the medium . The “bb” subscript on indicates that the radiance corresponds to a blackbody.In the plane parallel medium shown in Fig. 1, the radiance is generated in , lying at a depth of below the surface. As it propagates, the emitted light is attenuated by absorption, so only the fraction arrives at the surface. Therefore, the element of radiance arriving at the surface from this layer at depth is By integrating this over the thickness of the slab, from to , we obtain the total spectral radiance arriving at the surface as We simplify this as where we define as the bulk internal transmission of the entire slab, and we have assumed that the emission density is uniform throughout the volume. For inhomogeneous materials, or for materials with a temperature that varies throughout the volume, the depth dependence of prevents it from being pulled out of the integral as a constant.So far we have ignored any secondary reflections inside the volume: the portion of light that is not transmitted through the upper surface is instead reflected downward, and this radiance is . After traversing the volume of the slab (and attenuating due to absorption), it reaches the far side of the slab, where it reflects again. Traversing the slab one more time brings it back to the upper surface, where again a portion of the radiation transmits through the surface. This second component is therefore where is the two internal reflections, is the attenuation along that doubled path, and is the transmission through the surface. This process continues to occur, giving another contribution of and other higher-order terms. Summing all of these together gives the total upwelling radiance reaching the inside of the surface asSo far, we have only considered upwelling emission from the layers inside the slab. Because any thermal emitted radiance is isotropic, there is an equal amount emitted in the downward direction. For this portion, light reaching the upper surface of the slab must first undergo one reflection and suffer attenuation from traversing the volume of the medium, so the first portion of radiance reaching the upper surface is Likewise summing all contributions from multiple reflections, we obtain the total downwelling radiance reaching the inside of the top surface asThe total radiance reaching the underside of the surface is the sum of the upwelling and downwelling terms If the surface transmittance is , the spectral radiance that passes through the surface isHaving reached this point, we do not yet know how this emitted radiance relates to the standard method of representing surface-emitted radiance, viz., If we consider a material that is highly opaque (large ), then we know that the bulk transmission of the material is approximately zero (). In this situation, Kirchhoff’s law is unambiguous, giving . In our model of emission from the solid Eq. (14), we use Kirchhoff’s law to express the surface transmittance as . Thus, if we equate our calculated emitted radiance Eq. (14) with the familiar form of the object radiance Eq. (15), we have At normal incidence, we conclude that This means that the emission density divided by the absorption coefficient is equal to the blackbody radiance for the same conditions of wavelength and temperature.3.Straubel Theorem for Radiance Scaling Across an InterfaceSo far we have obtained the expression for the emitted radiance in terms of the internal angle but not yet in terms of the external angle . (Subscript 1 refers to quantities inside the solid, and subscript 2 refers to quantities in the external medium—usually air.) For this, we need to apply the angle-dependent refraction effects at the surface, making use of a well-known theorem in radiometry governing the scaling of the solid angle across an interface [shown in (Fig. 2)], apparently first discovered by Straubel in 1902.7,8 To start, from Snell’s law we know that . Taking the derivative of both sides with respect to , this becomes which leads to Using the definition of the differential solid angle, , we also write so that Substituting Eq. (19) into this, we obtain which is known as Straubel’s theorem.If we ignore for a moment the reflection losses on passing through a smooth interface, then the power entering through an element of solid angle is equal to the power exiting the interface through : Substituting Eq. (22) into this and incorporate surface transmission loss , we have This result, in which the emitted spectral radiance is proportional to the cosine of the exitance angle , is commonly described as Lambert’s cosine law.As we will see in Sec. 4 below, since the derivation of Straubel’s theorem assumes dielectric-only media, the theorem cannot hold in general for absorbing-material interfaces, yet how the theorem should be expressed in the more general case of absorbing media has apparently not been discussed in the literature. The silence on the subject is likely because media with significant are generally opaque, so they have little transmitted light to measure, making this a case of limited utility for most situations. Likewise, media with small transmit sufficient light for measurement, but the small value of minimizes the effect on the refraction angle. For intermediate cases, we need to use an expression giving the dependence of the refraction angle on the solid’s absorption index [given in Eq. (26)]. 4.Generalized Snell’s LawWhen working with absorbing materials, Snell’s law is used directly with the complex index : where we have explicitly labeled each of the complex variables with a tilde. The meaning of a complex-valued angle may be unclear, but it is convenient for the purposes of calculating the Fresnel reflection coefficients . However, to get the correct result for emitted light from the solid refracting into an ambient lossless medium (such as air), we can get a better feel for how the complex variables affect the refraction angles by working entirely in terms of real-valued geometry. For an interface between an absorbing emitter (the solid) and a nonabsorbing medium (the ambient air), we have and . For this case, the expression relating the real-valued angle inside the solid to the transmission angle in the ambient is given as where the Snell’s law correction factor is given by Bell et al. as9 where and is limited to cases in which the ambient medium is nonabsorbing ().In Fig. 3, a set of curves show the refraction angles given by Eq. (26), drawn together with a second set of dashed curves that are obtained by taking the real parts of and and calculated using the standard Snell’s law. We can see that, for angles close to normal (small ), the standard Snell’s law actually provides a good approximation to the generalized angle calculation. We can also see this mathematically, by taking the approximation in Eq. (27), from which we obtain (where means “both of and ”). The approximation is not valid when is comparable to or larger than . If we want to know the incident angle from within the solid for a given emission angle , we use Eq. (26) in reverse. Although this produces a transcendental function, this representation of the Snell correction factor is easy to apply due to the smooth monotonic relationship between and [see Fig. 3(b)]. Unfortunately, however, Eq. (27) does not provide us with the complex form of the angle inside the solid, , which we need in Sec. 5. The easiest way to get the complex is to work out Snell’s law Eq. (25) backward, either using a nonlinear solver or by interpolation. 5.Polarization and the Complex Fresnel EquationsSo far, all of our expressions for the emitted radiance and emissivity of the solid have been written in terms of the Fresnel intensity coefficients for surface transmittance and reflectivity . The Fresnel coefficients are more often written in terms of complex amplitudes and , from which the intensity coefficients are derived as and a similar expression for the transmission. (The superscript indicates a complex conjugate.) The most general form of the complex-valued Fresnel reflection and transmission, written for p-polarized light (electric field vibration in the plane of incidence at the surface) and for s-polarized light (electric field vibration perpendicular to the plane of incidence), are written as16,17 and the corresponding surface transmission amplitude coefficients are,17 Taking the magnitude squared of and , respectively, we obtain the intensity reflectivity and transmittance as18–21 and The fraction at the rightside of and in Eq. (31) expresses the change in the beam cross-sectional area upon refraction.For unpolarized light, it is common to define an unpolarized average reflectivity given as Because thermally emitted light from thick media is typically incoherent, the phase relationship between s and p reflected light can be considered random, that is, emitted radiance is generally partially linearly polarized, but not elliptically polarized, for thick media.An example set of reflectivity curves is shown in Fig. 4, at , for the cases of light incident from air and reflecting from an absorbing medium. Thus, in this case, and . Note that, although the water reflectivity curve looks like a curve for a nonabsorbing interface, it does not have a true Brewster angle.22 Instead, it has a pseudo-Brewster angle (for ) at an external angle of incidence of , where . A difficulty with using the complex Fresnel equations in our situation is that the radiance transfer approach allows us to know inside the solid, but it does not tell us the corresponding value for . If we naively attempt to use in place of in Eqs. (28) and (29), we will obtain unphysical reflection and transmission coefficient values. Thus, if we invert Eq. (27) for the generalized Snell’s law, we can calculate the correct angle of emission from the solid. Instead, we must invert Eq. (25) to obtain the corresponding complex angle for use in the Fresnel equations. Another phenomenon that may be unfamiliar to many is that, for absorbing incident media, interference between the incident wave and reflected wave causes nodes in the intensity distribution near the material interface, and this alters the behavior at the interface.23,24 As a result, the energy conservation requires the expression where can be positive or negative. This extra term, which has been termed the “mixed Poynting vector,” becomes negative when a maximum node (intensity peak) occurs near the interface itself and positive when a minimum node (intensity valley) occurs near the interface, causing a small increase or decrease in the amount of transmitted light depending on whether the node at the interface is a valley or a peak. This extra term is written as25,26 which at normal incidence simplifies to27Researchers who discuss the mixed Poynting vector almost universally assert that is “small” and therefore can be neglected.28–31 However, if one proceeds to substitute values for example materials into Eq. (34), one quickly finds that can be surprisingly large in the thermal infrared domain for a wide class of materials. It is quite common for metals to have , and in this regime, it is not unusual to have (i.e., transmission ) at higher angles of incidence, as shown in Fig. 5. Thus, what is missing from the literature is a description of the correct context for calling “small.” Figure 5 shows example curves of for light transmitted across an interface from an absorbing medium (in this case, silica and alumina) into air. For silica, we see that is small (), but it seems careless to refer to even a 2% change as negligible, especially as it occurs at normal incidence and not just at large angles of incidence. For alumina, the larger ratio between the absorption and refractive indices means that plays a larger role. Although at normal incidence, it reaches and at an angle of incidence . This does not violate Eq. (33) because we also find that and there. From these results, we can conclude that the nonlinear interaction between the incident wave and reflective wave at the interior boundary of the solid causes an enhancement (or in other cases, diminishment) of the internally reflected wave intensity. Because the absorption coefficient of these materials is large, this enhancement will be small in comparison with the absorption experienced by the wave as it propagates back toward the interior of the solid. For all but the thinnest of films, this likely makes the effect difficult to observe using the standard experimental approach of taking the ratio between the incident light and the transmitted light. However, semitransparent (and thus extremely thin) metal coatings on mylar films may provide the right combination of circumstances to observe this effect. However, as we see in Eq. (36), the surface transmittance is proportional to the object’s emittance. This means that any impact that the mixed Poynting vector has on the surface transmittance should be observable as a change in emittance. Either this impact has not been observed before, due to the limited range of wavelengths and angles over which it generally occurs, or the impact is negligible, for reasons that are not considered in the existing literature. Either case provides tasks for future work. 6.EmissivityThe emissivity is defined as the emitted radiance of a body in proportion to the emitted radiance of an equivalent (same temperature and same angle) blackbody emitter, so with Eq. (24) our radiative model gives the spectral emissivity as Here, “m” is used to index the polarization state as either p or s. When the solid is optically thick, so that as , we can see that this equation simplifies to —the emissivity of the object is equal to the object surface transmittance. This is from the perspective of light generated from inside the solid, refracting through the surface. On the other hand, we can look from the perspective of light in air incident on the surface of the solid. For an opaque solid, all light transmitted through the surface is absorbed, so this is also equivalent to , which is the basic expression of the Kirchhoff law. Figure 6 shows curves obtained from Eq. (36) for water, (alumina), (silica), graphite, and copper.At normal angles, when the ambient medium is air (, ), the emissivity Eq. (36) simplifies to which is shown in Fig. 7. The curve for may appear to be nonphysical because a nonabsorbing medium cannot be an emitter, but this result can be taken as an asymptotic limit of a medium in which is brought closer to zero, whereas the medium itself becomes semi-infinite in extent. If the limit is taken in this way, can be brought arbitrarily close to zero, if not quite all of the way to zero.From Eq. (27), we can calculate the incident angle from the emission angle , and by inverting Eq. (25) we can calculate the complex angle that is needed for the general Fresnel equations Eq. (31). Because is a function of the material thickness , we find that, for small values of thickness, increases approximately linearly with , but that as becomes large, approaches zero, and the emissivity approaches an asymptotic value. Following Edalatpour, we can define a “critical thickness” at which the emissivity reaches 99% of its asymptotic value.33 Taking the logarithm of Eq. (7) for . In the thermal infrared (3 to ), typical values for in metals are 30 to 90,34 giving critical thicknesses of 24 to 360 nm (i.e., less than a micron). Many nonmetal absorbing materials have ,35 giving to (i.e., well under 1 mm).In the outdoors, the surfaces exposed to water tend to be hydrophilic, so water tends to completely wet the surface as a film rather than drawing together into beads. This includes most glasses and ceramics and the most common metals (aluminum, steel, copper, etc.).36 After being saturated with rainwater and allowed to drain, many surfaces tend to retain a thin water film of roughly thickness.37 Because water at and has absorption indices of and 0.051, the critical thicknesses at which the water completely obscures the emission and reflection of the underlying surface is about . This thickness is only achieved as a layer fully covering a surface when the water cannot fully drain, so we expect that most wet surfaces are not fully obscured by their overlying water films. 6.1.Emission from Films and CoatingsAlthough Eq. (36) allows us to calculate the emissivity at any angle of incidence, for a homogeneous material, many materials naturally form a surface oxide layer upon exposure to air. Pure aluminum, for example, produces a natural amorphous oxide layer that is typically 2 to 4 nm thick.38 Naturally, the question arises about how these surface layers modify the appearance of these materials in the infrared. For nontransparent coatings, the object will take on the optical properties of the coating material if the layer is very thick. Conversely, it will maintain the optical properties of the substrate if the coating is quite thin. To define “thick” and “thin” for this usage, we once again make use of the critical thickness of the coating material to make a good approximation. References 39 to 41 develop a model for going beyond the approximate answer that we suggest here, for the case of transparent surface layers. For many coating layers, the coherence length of blackbody spectra exceeds the coating thicknesses. As a result, even though the emission of the substrate is incoherent, interference effects produced by multiple reflections within the surface layer will modify the coating optical properties. Modeling this requires the use of a coherent-light model that tracks the phase of a light wave passing through the coating, closely resembling the approach typically found in thin film optical models at shorter wavelengths. The resulting model accurately reproduces the modulations seen when monitoring the emissivity of a film as it is being grown on top of a substrate (as in Fig. 8). A general approach that allows for the thickness of the film to be of the same order as the coherence length of the emitted light requires use of a partial coherence theory of thin film optics.42,43 Such a partial coherence theory allows one to go smoothly between the coherent and incoherent extremes. We can also note that, although thermal emission is incoherent (has a very short coherence length), the appropriate coherence length to use when modeling a measurement is also determined by the angular size of the detection aperture and the spectral bandwidth of detection, rather than solely by the size and spectrum of the source alone. For small solid angles and narrow spectral bandwidths, even thermal light becomes coherent. Although a model of a transparent coating on an absorbing substrate is quite useful, transparent coatings are somewhat rare among materials in the thermal infrared. Most coatings, such as the natural oxide coatings that form on many metals upon exposure to air, are partially transparent and partially absorbing/emitting. To model this more general case, it is necessary to augment the coating model to include not only interference of partial reflections but also absorption and incoherent emission within the film. This is the approach taken in Ref. 33, which models the effect of growing a thin metal film or an absorbing dielectric film on top of an absorbing substrate. As we should expect, the model shows a continuous transition from the optical properties of the substrate to those of the film as the film increases in thickness. The change is most rapid when the film is at its thinnest and asymptotically approaches its final value when the film thickness exceeds (see Fig. 8). As with the transparent film, we also find oscillations in the emissivity curve as a function of thickness or wavelength. 7.Refining Kirchhoff’s LawThe basic form of Kirchhoff’s law expresses the equivalence of emission and absorption: for polarization state . In most of the literature, the conservation of energy for the interaction of light with media is invoked to say that probabilities of reflection, transmission, and absorption sum to 1: . Using the equivalence property of emission and absorption, this leads to the form of Kirchhoff’s law most commonly used in thermal imaging: The superscript ∘ is used to indicate that these are properties of the object (apparent properties) and not just a single surface (surface properties). For opaque media, , in which case Eq. (39) provides a simple method for obtaining the emissivity from the reflectivity .In the literature, Kirchhoff’s law is more often written as , i.e., with the Fresnel coefficients and rather than the apparent object reflectivity and transmittance. Doing so, however, ignores the effect of internal reflections within the object and internal scattering, so it is strictly only true for either opaque or diaphanous media. In addition, there is also a common assumption that, because the interface is an infinitesimally thin region, no absorption is possible at the surface. If true, this would allow the surface transmittance to be calculated from the complement of the reflectivity when the emissivity is negligible, . However, for absorbing media in which the waves are incident from within an absorbing medium, we have seen that there can also exist a small amount of (positive or negative) “apparent surface absorption” . A weakness of using the Kirchhoff law for deriving from is that its validity can be unclear for inhomogeneous media and particularly for nonisothermal media. The direct approach of modeling the emissivity from the spectral emission density provides better insight into the physical behavior. Using the incoherent sums of partial waves within the solid, McMahon derived expressions for the normal-angle emissivity, reflectivity, and transmittance as3 where and are the apparent spectral reflectivity and spectral transmittance of the object, respectively, and and retain their meaning as the Fresnel intensity reflectivity and transmittance of the interface between the solid and the ambient medium, respectively. The wavelength and angle dependence of and were suppressed for clarity. Adding all three together, with some algebra we obtain the Kirchhoff law Thus, for a homogeneous and isothermal medium, we have been able to derive Kirchhoff’s law as a consequence of the Fresnel equations, for all angles and all polarization states. One difference in Eq. (40) from McMahon’s results is that the expression for uses in place of McMahon’s . These are usually asserted to be equivalent, but we have seen from Eq. (33) that, although they are typically close, they are not exactly equivalent due to the presence of the mixed Poynting vector term . Thus, it seems that the mixed Poynting vector may point to a deviation from Kirchhoff’s law, for situations in which is significant.7.1.Nonisothermal MediaTreatments of Kirchhoff’s law generally state that the law is valid when the emitting object is in local thermodynamic equilibrium (LTE) and may further add that this is not satisfied when an object is not isothermal. Because LTE is never exactly satisfied in practice, it is even more important to be able to say quantitatively to what degree a given situation permits the law to provide accurate results, that is, rather than evaluating a situation as valid or invalid for applying the law, we should instead express how closely we can expect the law to hold. First, we start with a heuristic analysis to gain insight into the shape of the problem before going on to developing the path integral for the nonisothermal case. We can start with an example of a painted metal plate that is exposed to the open air and is cooled by air convection currents and radiative losses. The paint coating has an absorption index of at wavelength . In this case, the skin depth of the coating is . Regardless of the status of Kirchhoff’s law, we can expect that most of the emission from the solid will take place within this narrow region just below the surface. If the paint layer is nonisothermal, then any thermal gradient will have to be significant within this skin depth. Because most paint coatings are 5 to thick, this narrow region is the outermost surface of the paint. If the temperature difference between the warm plate (the “substrate” of the coating) and the surface of the coating is 5°C and the coating is thick, then the temperature delta across the skin depth is only . This gives the maximum possible effect that the thermal gradient within this coating can have on the emissivity, making the surface appear as if it were 35 mK warmer than it really is. In practice, the emission by deeper layers is attenuated by absorption before reaching the surface, so the change in emission is even smaller than this and therefore below the noise level of most thermal infrared imagers. With this result in mind, we can say that painted metal surfaces used as blackbody simulators will likely have an apparent radiative temperature that differs by less than ∼35 mK from the temperature of the underlying metal substrate as long as the high-emissivity coating layer is kept to thicknesses near . To make this more quantitative, we develop a model for the emission of a nonisothermal solid. The spectral emission density becomes a function of due to the linear depth-dependence of temperature within the solid: for surface temperature and temperature gradient (in units of ), so Inserting this into Eq. (4) gives the element of radiance emitted from within a layer of thickness (see Fig. 1). Integrating this over the thickness of the solid, from to , we obtain the spectral radiance arriving at the surface as Although this is a complicated integral, the expression is close to linear when the gradient term is small in comparison with the surface temperature, . If we take the Taylor series expansion of this expression about the point , then we obtain for . Putting this expansion into Eq. (43), we obtain where is the radiance obtained for an isothermal solid at . When the solid is semi-infinite in extent (), then this simplifies to where is the nonisothermal modifier term. Although following a different derivation, this last result agrees with that given by Ref. 44. We can note that , so that materials with large absorption index tend to have values closely agreeing with unless the thermal gradient is extreme.Figure 9 shows example curves of for different paint coating thicknesses from 2 to and a temperature difference between the two sides of the coating (the upper surface exposed to air and the lower surface touching the substrate). For very thin layers of paint that still have a thermal delta of 4°C, we see that the nonisothermal effect on the emittance can reach as high as 10%. However, for more ordinary thermal deltas, the value of quickly goes below our ability to measure. Because there is a wide variety of base media and additives in paints, there is a corresponding wide range of refractive and absorption indices. Oppenheim gives a representative measurement of , whereas Richmond gives , and Tyo finds .6,14,45 These values fall within the range of our own experimental measurements, in which the high- paints generally contain a large number of pigment particles, making them appear black in the visible spectrum. The paint layers simulated in Fig. 9 assume intermediate values of . In addition to the requirement that the emitting object be isothermal, one sometimes encounters the statement that the object must also be homogeneous for Kirchhoff’s law to hold.41 However, it is clear that inhomogeneous objects do not by themselves cause a violation of the equilibrium condition for Kirchhoff’s law because Kirchhoff’s law has been applied successfully to multilayer films and to composite materials such as concrete. 7.2.Nonisothermal Behavior in Paints and GlassesFigure 9 shows that a 4°C temperature difference across any paint coating more than thick produces such a small change in emissivity that it will be difficult to measure. Although this simulation assumes a temperature delta of 4°C, the author’s experience with thin paint coatings in the lab indicate that deltas of 0.2°C to 0.5°C are more likely (using contact thermisters to measure substrate and coating surface temperature). As a result, for blackbody simulator paint coatings, we should be able to neglect the effect of deviating from isothermal sources. If the thermal gradient becomes steeper and the absorption coefficient becomes smaller (so that deeper regions participate in emission), the change in emissivity can become significant. For example, for a coating layer with , a temperature delta of 2°C across a layer of thickness (i.e., a gradient of ) produces a 1.5% change in emissivity at normal incidence, falling off to 1.2% near the grazing incidence. We can see from Eq. (48) that is proportional to the thermal gradient and inversely proportional to the absorption index . Thus, increases linearly with the thermal change across the emitter’s skin depth. If we scale by , then we can use the same curve for all values of , , and , as shown in Fig. 10. The proportionality also shows why metals present no difficulty for Kirchhoff’s law—the large value and low of most metals make negligible in almost all circumstances. Using the curves shown in Fig. 10, we can interpret the results in the following way. The blue curve shows that, if there is a 4°C temperature difference between the substrate and the surface of the coating layer, for the coating layer thickness of , then the deviation of emission from the isothermal case will be about 5%. One can likely achieve these kinds of conditions by heating the substrate to a high temperature (such as 100°C using the painted kitchen pot shown in Fig. 3 of Part II of this review) while cooling the surface with a fan blowing cool air. This illustrates the conditions under which the nonisothermal effect becomes important. In the case of silica () and alumina () glasses, the emission skin depth at is 2.0 and , respectively. For alumina, however, happens to be close to an absorption peak, so it is not necessarily representative of the material’s overall behavior in the thermal infrared. In the midwave infrared (MWIR) (3 to ), for example, the absorption coefficient for silica drops to —three orders of magnitude lower than the value at . Similarly, the absorption index for alumina drops to 0.0018 and for paint (if we take polymethyl methacrylate as representative of a plastic-based paint) to 0.00065. Thus, in the MWIR at , the skin depths of silica, alumina, and acrylic-base paint become 2.3, 0.5, and 1.2 mm, respectively. The lower absorption coefficients in the MWIR versus longwave infrared (LWIR) mean that nonisothermal effects will be more likely in the midwave regime. Although most natural materials, under normal circumstances, do no deviate significantly from Kirchhoff’s law, one category of natural materials that consistently exhibits significant departures from Kirchhoff’s law is layered fine particulate such as quartz sand powder.46 Reference 46 measured up to 6% deviations from Kirchhoff’s law for dry unpacked quartz dust. The air gaps between particles reduce the thermal conductivity of the material while preventing convective cooling; at the same time, this allows for sunlight to penetrate deeply and heat up the underlying particle layers, producing a large thermal gradient. Another material having similar behavior is lunar soil.47,48 One might expect fresh snow to have similar properties, due to having a similar physical structure, but experimental measurements have shown that snow exhibits a close correspondence to Kirchhoff’s law.49 One can of course envisage cases in which nonisothermal behavior has a predominant effect on emission. For example, most “transparent” materials in the infrared actually have a small amount of absorption present. We can imagine a thick slab of such a medium, in which the slab thickness is larger than the material’s thermal penetration depth, but much shorter than the optical absorption depth. Within such a thick slab, emission can take place over a large range of distances, allowing the various emitting layers to have large differences in temperature. 7.3.Kirchhoff’s Law for GasesWhen modeling gases in the thermal infrared, it is common to write Kirchhoff’s law simply as absorption equals emission: . Although one can interpret this as equivalent to Eq. (41) when the reflectivity is zero and the transmission is , equating absorption and emission is actually the most fundamental form of Kirchhoff’s law.50 As a result, it is not necessary to know the scattering (the diffuse reflection) by the gas cloud to calculate the emission. Instead, one needs to know only the gas temperature and its absorbance along the relevant line of sight. For a gas cloud of thickness along the line of sight , the absorbance is51 for gas concentration (i.e., number density) and absorption cross-section . Equating the emissivity with the absorbance , the emitted radiance of a gas at uniform temperature is52 When the gas is nonuniform in temperature, however, one must divide the gas cloud into a series of separate parcels, each with their own locally uniform temperature and their own local absorption value (for light passing from one side to the other side of the parcel along the line of sight).53–55 Although this substantially increases the complexity of the calculation, it is a straightforward process to model. In the unlikely case of a linear thermal gradient along the line of sight through the gas cloud, one reaches a result analogous to that of the linear gradient in a solid given in Eq. (46). What is different from solids is that most gases are weak absorbers/emitters, so that light can propagate through a large depth of gas before extinction. In addition, the lack of direct thermal conduction means that large differences in temperature between parcels of gas are possible. As a result, modeling a gas cloud as an isothermal body with uniform absorbance/emissivity via Kirchhoff’s law is often subject to larger errors than in the case of solids.8.Rough Surfaces and Scattering MediaSo far we have considered only interfaces with smooth boundaries. When the boundaries are rough, incident light from one angle is scattered into a range of transmission and reflection angles, requiring a modification to the emissivity model that we have used so far. Jordan et al. provide an approach for doing this by modeling the orientation angle of the interface as comprising a probability distribution over angles.18,56,57 In general, we can expect that, as the surface roughness becomes large, multiple reflections among surface facets increase, causing a reduction in reflectivity and a consequent increase in the emissivity.14 The effect on emission polarization is more complex. Following Chang et al., we can model the emitting surface roughness as consisting of triangular corrugations of a given angle (see Fig. 11). From the standpoint of the viewer, the emission angle is , but the Fresnel equations for emission use the local surface normal to calculate the polarization properties. Thus, if the corrugations are ordered and not randomly oriented, this allows for polarized emission even at , as shown in Fig. 11. For emission angles above the local tilt , the surface emission is affected by shadowing, in which part of the facet emits directly into the external medium, but emission from the remainder of the facet is reflected from the opposite face, inducing a change in polarization on reflection.58,59 Although this shadowing effect makes analytical analysis of the problem difficult, numerical simulations have shown that the corrugation tends to increase overall emissivity, but the shadowing effect tends to reduce the degree of polarization in the emitted light.18 Goldstein and Chenault show polarimetric measurements of reflection from various paints for a range of angles of incidence.60,61 Although the measurements were taken in the visible and near-infrared regions, they also contained very useful tables of the RMS surface roughness for each paint type, showing that the RMS roughness ranges from 50 nm for glossy paints to as much as for matte paints.62 Germer considers the effect of surface microroughness—surface scattering features that are small with respect to the light wavelength—and its effect on optical scattering and observed polarization.63,64 Although Germer’s simulations provide only for the case of visible light, the overall approach is equally applicable to infrared thermal emission and reflection and indicate that, because microroughness and subsurface defects have a somewhat different angular dependence of polarization, the two effects can be separated from one another with careful measurements. When a material is itself inhomogeneous, so that light propagating through the medium experiences a significant amount of scatter as it traverses the volume, the reflectivity and emissivity values calculated using a homogeneous material model are no longer reliable. One approach to take in this case is to form a Monte Carlo model of the rays propagating through the media, but this does not provide much physical insight to apply across various different media. When the scattering is large, however, we can take the diffuse light approximation, in which case we can apply the Kubelka–Munk approach to calculating the reflectivity.65 In this approach, light is assumed to be evenly distributed in angle, so the calculated reflectivity is angle independent. This allows the problem to be simplified into a one-dimensional differential equation. Making use of the Kirchhoff law, we are then able to infer the emissivity from . In the diffuse light approach, we imagine the light as being composed of two streams, an upwelling flux and a downwelling flux (shown as red and blue, respectively, in Fig. 12). Each stream is subjected to absorption and scattering effects, for which we use coefficients and . If we examine the behavior of the two light streams as they traverse a layer of differential thickness inside the medium, we see that upwelling flux suffers from absorption loss () and scattering loss () upon traversing the layer, but gains from backscattered downwelling light (). The downwelling flux likewise suffers from absorption and scattering but gains from backscattered upwelling light. As a result, the differential changes in fluxes across the layer, defined as are given by the pair of coupled differential equations The solutions for these equations are45 for and , are constants to be determined by the boundary conditions of the problem.Although the and coefficients are intended as experimentally determined parameters, researchers have established that in many cases we can approximate the absorption coefficient as and the backscattering coefficient as , where we model the medium as composed of discrete scatterers of scattering cross-section and number density .66,67 A more detailed calculation of and that takes into account specific material properties is presented in Refs. 68 and 69. The steps toward calculating and are given in detail in Ref. 45, so we skip to the solution: where and is the “external” reflectivity for light incident from the ambient medium onto the surface of the coating layer; is the “internal” reflectivity for light incident on the surface of the coating layer from inside the coating; and is the “substrate” reflectivity for light incident from the coating onto the substrate. The substrate itself is assumed to be fully opaque. Each of these three reflectivities , , and are hemispherical average reflectivities obtained from the Fresnel equations. In the solution given in Eq. (55), the material reflectivity is the diffuse reflectance, so it is independent of angle or polarization.For a direct model of scattering media that does not derive the emissivity from applying Kirchhoff’s law to the reflectivity, it becomes necessary to account for the emission density within the medium itself. Xu and Shen show that this adds a new source term to each of the equations in Eq. (53), requiring modification of the solution.70 9.ConclusionOne of the primary aims of this review was to gather together material illustrating how the elements in the theory of emission polarization can be derived from the Fresnel laws and volumetric blackbody radiation inside absorbing materials. We also considered the generalization of the Straubel theorem to absorbing materials, quantified the requirements for Kirchhoff’s law to hold, and discussed the potential role of the mixed Poynting vector on emission. The physics of thermal emission is likely one of the few places in which the mixed Poynting vector may play a nontrivial role in practical measurement situations. Along the way through our survey, we also provided some of the answers to the questions that originally motivated gathering together the theory in one place. (Other answers await development until Part II of this review.) For example, we found that a layer of water fully obscures its underlying substrate when 1 mm thick or more. For hydrophilic surfaces in which water spreads widely without collecting into beads, water only forms layers this thick when it cannot drain. Wet surfaces that can drain retain only very thin layers of water, in which case the emission/reflection from the surface will be predominantly from the substrate and not the water layer. We also found that high-emissivity paints can be considered to fully obscure their substrates when over about thick, though this depends a lot on the type and density of colorants used in the paint. Surprisingly, a fully obscuring paint at visible wavelengths is generally a good guide that the paint will also be fully obscuring (or nearly so) in the thermal infrared as well. There is, however, some concern about the use of paints for blackbody simulators. If the heater underneath the paint layer is pushed to a high temperature, so that the thermal gradient across the paint becomes large, this gradient can cause measurable deviations from Kirchhoff’s law. Finally, the McMahon–Gardon model of emission (Sec. 2) shows that the emission from a thin (partially transparent) plastic film is not simple. Most plastics have a complex absorption spectrum in the thermal infrared, particularly in the LWIR. Although some wavelengths are highly transparent, others are significantly absorbing, making their spectral dependence a major factor. In addition, any experimental measurements are prone to differences in composition and any coloring agents (such as carbon particles) or impurities. For example, clear plastic bags are more transparent in the infrared than black plastic bags because the latter contain a large concentration of absorbing carbon particles. 10.Appendix A: Gardon’s IntegralIn Sec. 2, we said that a 1D integral along the ray is sufficient to model the emissivity from a solid slab. We can confirm this with a rigorous treatment of the differential solid angle elements and the differential surface area element. With the setup shown in Fig. 1, the volume element is part of an annular ring of radius , with radial thickness and axial thickness . From the center of the ring, the elemental volume subtends an angle of . In cylindrical polar coordinates , we write the elemental volume as However, a more convenient representation is to rewrite this in terms of and rather than and . Using the triangle formed by , , and , The geometry shown in Fig. 1 gives so and Therefore, in terms of the more convenient coordinates , the volume element isAgain making use of the geometrical relationships shown in Fig. 1, we write the solid angle element , which is the solid angle subtended by surface area element from , as Similarly, for the solid angle element defined within the blue lines drawn in Fig. 1(b) representing the solid angle subtended by from , we have The element of volume emits power in the direction of : The power arriving at from is less than this because of absorption along the path. Including path absorption, we have The corresponding elemental radiance reaching the surface from inside the solid is given by the power per units area per unit solid angle : With this expression, we can solve for the radiance arriving at from within the solid by integrating over all layers:AcknowledgmentsThe author is grateful to Photron Inc. for allowing the use of their FLIR x6900sc InSb camera in the MWIR polarimetry experiments. The author would also like to thank Shuhei Shibata for help with conducting the FLIR A35 LWIR polarimetry experiments. A portion of this research was funded by a grant from Denso Corp. and by JSPS KAKENHI (Grant No. JP 20K04516). ReferencesR. A. Millikan,
“A study of the polarization of the light emitted by incandescent solid and liquid surfaces,”
Phys. Rev., 3 81
–99
(1895). PHRVAO 0031-899X Google Scholar
R. A. Millikan,
“A study of the polarization of the light emitted by incandescent solid and liquid surfaces. II,”
Phys. Rev., 3 177
–192
(1895). PHRVAO 0031-899X Google Scholar
H. O. McMahon,
“Thermal radiation from partially transparent reflecting bodies,”
J. Opt. Soc. Am., 40 376
–380
(1950). https://doi.org/10.1364/JOSA.40.000376 JOSAAH 0030-3941 Google Scholar
R. Gardon,
“The emissivity of transparent materials,”
J. Am. Ceram. Soc., 39 278
–285
(1956). https://doi.org/10.1111/j.1151-2916.1956.tb15833.x JACTAW 0002-7820 Google Scholar
O. Sandus,
“A review of emission polarization,”
Appl. Opt., 4 1634
–1642
(1965). https://doi.org/10.1364/AO.4.001634 APOPAI 0003-6935 Google Scholar
J. S. Tyo et al.,
“The effects of thermal equilibrium and contrast in LWIR polarimetric images,”
Opt. Express, 15 15161
–15167
(2007). https://doi.org/10.1364/OE.15.015161 OPEXFF 1094-4087 Google Scholar
R. Straubel,
“Über einen allgemeinen Satz der geometrischen Optik und einige Anwendungen [About a general theorem of geometrical optics and some applications],”
Verhandl. Deut. Physik. Ges., 4 328
–334
(1902). Google Scholar
A. Gershun,
“The light field,”
J. Math. Phys., 18 51
–151
(1939). https://doi.org/10.1002/sapm193918151 Google Scholar
R. J. Bell et al.,
“Generalized laws of refraction and reflection,”
J. Opt. Soc. Am., 59 187
–189
(1969). https://doi.org/10.1364/JOSA.59.000187 JOSAAH 0030-3941 Google Scholar
J. Kischkat et al.,
“Mid-infrared optical properties of thin films of aluminum oxide, titanium dioxide, silicon dioxide, aluminum nitride, and silicon nitride,”
Appl. Opt., 51 6789
–6798
(2012). https://doi.org/10.1364/AO.51.006789 APOPAI 0003-6935 Google Scholar
G. M. Hale and M. R. Querry,
“Optical constants of water in the 200-nm to wavelength region,”
Appl. Opt., 12 555
–563
(1973). https://doi.org/10.1364/AO.12.000555 APOPAI 0003-6935 Google Scholar
S. Popova, T. Tolstykh and V. Vorobev,
“Optical characteristics of amorphous quartz in the region,”
Opt. Spectrosc., 33 444
–445
(1972). OPSUA3 0030-400X Google Scholar
M. R. Querry,
“Optical constants,”
(1985). Google Scholar
U. P. Oppenheim and Y. Feiner,
“Polarization of the reflectivity of paints and other rough surfaces in the infrared,”
Appl. Opt., 34 1664
–1671
(1995). https://doi.org/10.1364/AO.34.001664 APOPAI 0003-6935 Google Scholar
A. D. Rakic,
“Algorithm for the determination of intrinsic optical constants of metal films: application to aluminum,”
Appl. Opt., 34 4755
–4767
(1995). https://doi.org/10.1364/AO.34.004755 APOPAI 0003-6935 Google Scholar
O. S. Heavens, Optical Properties of Thin Solid Films, Dover(1991). Google Scholar
D. H. Goldstein, Polarized Light, 2nd ed.Marcel Dekker(2003). Google Scholar
P. C. Y. Chang et al.,
“Importance of shadowing and multiple reflections in emission polarization,”
Waves Rand. Media, 12
(1), 1
–19
(2002). https://doi.org/10.1088/0959-7174/12/1/301 Google Scholar
S. J. Orfanidis,
“Electromagnetic waves and antennas,”
(2016) https://www.ece.rutgers.edu/~orfanidi/ewa/ Google Scholar
S. J. Byrnes,
“Multilayer optical calculations,”
(2016). https://arxiv.org/abs/1603.02720 Google Scholar
F. L. Pedrotti, L. M. Pedrotti and L. S. Pedrotti, Introduction to Optics, 3rd ed.Cambridge University Press(2018). Google Scholar
J. A. Shaw,
“Degree of linear polarization in spectral radiances from water-viewing infrared radiometers,”
Appl. Opt., 38 3157
–3165
(1999). https://doi.org/10.1364/AO.38.003157 APOPAI 0003-6935 Google Scholar
M. A. Dupertuis, B. Acklin and M. Proctor,
“Generalized energy balance and reciprocity relations for thin-film optics,”
J. Opt. Soc. Am. A, 11 1167
–1174
(1994). https://doi.org/10.1364/JOSAA.11.001167 JOAOD6 0740-3232 Google Scholar
G. P. Ortiz and W. L. Mochán,
“Nonadditivity of Poynting vector within opaque media,”
J. Opt. Soc. Am. A, 22 2827
–2837
(2005). https://doi.org/10.1364/JOSAA.22.002827 JOAOD6 0740-3232 Google Scholar
H. A. Macleod, Thin Film Optical Filters, 3rd ed.Institute of Physics, Philadelphia
(2001). Google Scholar
H. Weber,
“The Fresnel equations for lossy dielectrics and conservation of energy,”
J. Mod. Opt., 61 1219
–1224
(2014). https://doi.org/10.1080/09500340.2014.928375 JMOPEW 0950-0340 Google Scholar
J. F. Lodenquai,
“Electromagnetic wave propagation in media with complex refractive indices and transmission through absorbing films,”
Am. J. Phys., 59 248
–254
(1991). https://doi.org/10.1119/1.16571 AJPIAS 0002-9505 Google Scholar
F. R. Keßler,
“Optics with gradients of free carrier concentration,”
Festkörperprobleme, 26 277
–308 Springer(1986). Google Scholar
C. von Fragstein,
“History of the mixed Poynting vector,”
Abh. Braunschw. Wiss. Ges., 34 25
–29
(1987). ABWGAZ 0068-0737 Google Scholar
A. MacLeod,
“The mixed Poynting vector,”
Soc. Vac. Coat. Bull., 14
(Summer), 30
–34
(2014). Google Scholar
N. Sun et al.,
“Effects of non-synchronized variations of electric and magnetic properties on transmitted waves at lossy interface,”
J. Quantum Spectrosc. Radiat. Transf., 138 50
–59
(2014). https://doi.org/10.1016/j.jqsrt.2014.01.018 JQSRAE 0022-4073 Google Scholar
J. R. Howell, R. Siegel and M. P. Mengüç, Thermal Radiation Heat Transfer, 5th ed.CRC Press(2010). Google Scholar
S. Edalatpour and M. Francoeur,
“Size effect on the emissivity of thin films,”
J. Quant. Spectrosc. Radiat. Transf., 118 75
–85
(2013). https://doi.org/10.1016/j.jqsrt.2012.12.012 JQSRAE 0022-4073 Google Scholar
III R. O. Carter,
“The effect of metal optical constants on far- and mid-infrared reflection-absorption,”
Spectrochim. Acta A, 47 551
–561
(1991). https://doi.org/10.1016/0584-8539(91)80071-P Google Scholar
Handbook of Optical Constants of Solids, Academic Press(1998). Google Scholar
J. D. Bernardin et al.,
“Contact angle temperature dependence for water droplets on practical aluminum surfaces,”
Int. J. Heat Mass Transf., 40 1017
–1033
(1997). https://doi.org/10.1016/0017-9310(96)00184-6 Google Scholar
B. M. Gallaway et al.,
“Influence of water depths on friction properties of various pavement types,”
(1974). Google Scholar
N. F. Mott,
“Oxidation of metals and the formation of protective films,”
Nature, 145 996
–1000
(1940). https://doi.org/10.1038/145996a0 Google Scholar
R. Raman and H. K. Sehgal,
“Emissivity of substrate supported films,”
Energy Convers. Manage., 22 213
–217
(1982). https://doi.org/10.1016/0196-8904(82)90044-9 ECMADL 0196-8904 Google Scholar
G. Neuer and F. Guntert,
“In situ measurements of layer thickness during oxidation of titanium,”
Thermochim. Acta, 133 299
–304
(1988). https://doi.org/10.1016/0040-6031(88)87173-9 THACAS 0040-6031 Google Scholar
P. Pigeat, D. Rouxel and B. Weber,
“Calculation of thermal emissivity for thin films by a direct method,”
Phys. Rev. B, 57 9293
–9300
(1998). https://doi.org/10.1103/PhysRevB.57.9293 Google Scholar
K. Richter, G. Chen and C.-L. Tien,
“Partial coherence theory of multilayer thin-film optical properties,”
Opt. Eng., 32 1897
–1903
(1993). https://doi.org/10.1117/12.147152 Google Scholar
C. F. Anderson and Y. Bayazitoglu,
“Radiative properties of films using partial coherence theory,”
J. Thermophys. Heat Transf., 10 26
–32
(1996). https://doi.org/10.2514/3.748 JTHTEO 0887-8722 Google Scholar
A. H. Madjid,
“Spectral emissivity of nonisothermal semi-infinite slabs of transparent materials,”
J. App. Phys., 44 5423
–5426
(1973). https://doi.org/10.1063/1.1662168 Google Scholar
J. C. Richmond,
“Relation of emittance to other optical properties,”
J. Res. NIST, 67 217
–226
(1963). Google Scholar
A. R. Korb, J. W. Salisbury and D. M. D’Aria,
“Thermal-infrared remote sensing and Kirchhoff’s law 2. Field measurements,”
J. Geophys. Res. B, 104 15339
(1999). https://doi.org/10.1029/97JB03537 JJGBDU Google Scholar
A. E. Wechsler and P. E. Glaser,
“Pressure effects on postulated lunar materials,”
Icarus, 4 335
–352
(1965). https://doi.org/10.1016/0019-1035(65)90038-2 ICRSA5 0019-1035 Google Scholar
L. Millán, I. Thomas and N. Bowles,
“Lunar regolith thermal gradients and emission spectra: modeling and validation,”
J. Geophys. Res. E, 116 12003
(2011). https://doi.org/10.1029/2011JE003874 Google Scholar
M. Hori et al.,
“Modeling angular-dependent spectral emissivity of snow and ice in the thermal infrared atmospheric window,”
Appl. Opt., 52 7243
–7255
(2013). https://doi.org/10.1364/AO.52.007243 APOPAI 0003-6935 Google Scholar
M. A. Weinstein,
“On the validity of Kirchhoff’s law for a freely radiating body,”
Am. J. Phys., 28 123
–125
(1960). https://doi.org/10.1119/1.1935075 AJPIAS 0002-9505 Google Scholar
C. F. Bohren and E. E. Clothiaux, Fundamentals of Atmospheric Radiation, Wiley(2006). Google Scholar
S. E. Golowich and D. G. Manolakis,
“Performance limits of LWIR gaseous plume quantification,”
Proc. SPIE, 8048 80481F
(2011). https://doi.org/10.1117/12.884360 PSISDG 0277-786X Google Scholar
K. Stamnes et al.,
“Numerically stable algorithm for discrete-ordinate-method radiative transfer in multiple scattering and emitting layered media,”
Appl. Opt., 27 2502
–2509
(1988). https://doi.org/10.1364/AO.27.002502 APOPAI 0003-6935 Google Scholar
G. E. Thomas and K. Stamnes, Radiative Transfer in the Atmosphere and Ocean, Cambridge(1999). Google Scholar
N. Hagen,
“Survey of autonomous gas leak detection and quantification with snapshot infrared spectral imaging,”
J. Opt., 22 103001
(2020). https://doi.org/10.1088/2040-8986/abb1cf Google Scholar
D. L. Jordan, G. D. Lewis and E. Jakeman,
“Emission polarization of roughened glass and aluminum surfaces,”
Appl. Opt., 35 3583
–3590
(1996). https://doi.org/10.1364/AO.35.003583 APOPAI 0003-6935 Google Scholar
J. A. Shaw and C. Marston,
“Polarized infrared emissivity for a rough water surface,”
Opt. Express, 7 375
–380
(2000). https://doi.org/10.1364/OE.7.000375 OPEXFF 1094-4087 Google Scholar
B. G. Smith,
“Geometrical shadowing of a random rough surface,”
IEEE Trans. Antennas Propag., 15
(5), 668
–671
(1967). https://doi.org/10.1109/TAP.1967.1138991 IETPAK 0018-926X Google Scholar
J. A. Ogilvy,
“Wave scattering from rough surfaces,”
Rep. Prog. Phys., 50 1553
–1608
(1987). https://doi.org/10.1088/0034-4885/50/12/001 RPPHAG 0034-4885 Google Scholar
D. H. Goldstein and D. B. Chenault,
“Spectropolarimetric reflectometer,”
Opt. Eng., 41 1013
–1020
(2002). https://doi.org/10.1117/1.1467933 Google Scholar
D. H. Goldstein,
“Polarization measurements of automobile paints,”
Proc. SPIE, 6972 69720V
(2008). https://doi.org/10.1117/12.780054 PSISDG 0277-786X Google Scholar
D. H. Goldstein,
“Polarimetric characterization of federal standard paints,”
Proc. SPIE, 4133 112
–122
(2000). https://doi.org/10.1117/12.406618 PSISDG 0277-786X Google Scholar
T. A. Germer,
“Angular dependence and polarization of out-of-plane optical scattering from particulate contamination, subsurface defects, and surface microroughness,”
Appl. Opt., 36
(33), 8798
–8805
(1997). https://doi.org/10.1364/AO.36.008798 APOPAI 0003-6935 Google Scholar
T. A. Germer and C. C. Asmail,
“Polarization of light scattered by microrough surfaces and subsurface defects,”
J. Opt. Soc. Am. A, 16
(6), 1326
–1332
(1999). https://doi.org/10.1364/JOSAA.16.001326 JOAOD6 0740-3232 Google Scholar
P. Kubelka,
“New contributions to the optics of intensely light-scattering materials. Part I,”
J. Opt. Soc. Am., 38 448
–457
(1948). https://doi.org/10.1364/JOSA.38.000448 JOSAAH 0030-3941 Google Scholar
P. S. Mudgett and L. W. Richards,
“Multiple scattering calculations for technology,”
Appl. Opt., 10 1485
–1502
(1971). https://doi.org/10.1364/AO.10.001485 APOPAI 0003-6935 Google Scholar
J. H. Nobbs,
“Kubelka–Munk theory and the prediction of reflectance,”
Rev. Prog. Coloration Rel. Top., 15 66
–75
(1985). https://doi.org/10.1111/j.1478-4408.1985.tb03737.x Google Scholar
L. Liu et al.,
“Calculation of emittance of a coating layer with the Kubelka–Munk theory and the Mie-scattering model,”
J. Opt. Soc. Am. A, 22 2424
–2429
(2005). https://doi.org/10.1364/JOSAA.22.002424 JOAOD6 0740-3232 Google Scholar
M. P. Diebold, Application of Light Scattering to Coatings: A User’s Guide, Springer(2014). Google Scholar
W. Xu and S. C. Shen,
“Infrared radiation and reflection in an inhomogeneous coating layer on a substrate,”
Appl. Opt., 31 4488
–4496
(1992). https://doi.org/10.1364/AO.31.004488 APOPAI 0003-6935 Google Scholar
Biography |