Polytetrafluoroethylene (PTFE) is known as a strongly scattering and diffusely reflecting material. The diffuse reflection is mainly due to the volume scattering of light inside the material. When the light enters a volume-scattering material, it is scattered by microstructures and propagates in random directions before escaping the medium. Due to its characteristics of being a nearly diffuse reflector, PTFE has been used as a diffuse-reflectance standard in the spectral range from by the National Institute of Standards and Technology (NIST),1 as a calibration standard for onboard sensors on satellites in remote sensing,2, 3, 4, 5 as whiteness standards in colorimetry,6, 7 and as a coating layer in integrating spheres.8 In addition, since both PTFE and most biological tissues are strongly scattering materials, researchers have used PTFE as tissue phantoms to simulate layers of skin for the study of burn depth.9
Light scattering in a turbid medium such as PTFE and tissue has been studied extensively. The angle-resolved scattering of a turbid medium could not be predicted until the more general radiative transfer equation (RTE) was developed.10 The common techniques for solving RTE include Chandrasekhar’s and functions,10 the discrete-ordinates method,11 the Monte Carlo method,12, 13, 14, 15 and the adding-doubling method.16 To describe light propagation in a PTFE slab using the RTE, the scattering coefficient , absorption coefficient , and the scattering phase function need to be known. The determination of these parameters typically requires several measurements, among which are:17 1. directional-hemispherical reflectance ; 2. directional-hemispherical transmittance ; 3. absorptance ; 4. collimated light transmittance ; and 5. angular distribution of scattered light from a sample thin enough to ensure that only single scattering occurs. The double-integrating-sphere method has been widely used for measuring and .18, 19 With the additional measurement of , these parameters can be determined with an inverse method.20, 21
Using an approach similar to the double-integrating-sphere method, Huber, Heitz, and Bauerle22 reported the scattering coefficient, absorption coefficient, and the asymmetric parameter in the Henyey–Greenstein phase function of PTFE films with thicknesses from . The scattering parameters of PTFE reported in their study are very similar to those of biological tissues.17 For example, the scattering and absorption coefficients at were and , respectively. However, it appears that these authors did not distinguish the scattered light in the parallel direction from the collimated light transmission. For PTFE samples with thicknesses greater than , the scattered light would dominate the transmittance even in the direction parallel to the incident light; in other words, since is essentially zero, the collimated light does not contribute to the detector signal. Unless the PTFE films can be made thinner than approximately four times the mean free path of the incident light, it is not feasible to unambiguously determine the scattering coefficient because the collimated light would be much weaker than the forward scattered light. This difficulty also applies to the determination of scattering phase function by measuring the angular distribution of scattered light from a sample because the single scattering condition cannot be guaranteed. Furthermore, the measurement of absorptance of PTFE using an integrating sphere is problematic because the absorption coefficient of this material is so small that the uncertainty of integrating sphere measurements often makes the absorptance imperceptible.
In the present study, the bidirectional reflectance distribution function (BRDF) and bidirectional transmittance distribution function (BTDF) of five PTFE films, whose thicknesses vary from , were measured using a laser scatterometer at . Integrating the measured BRDF and BTDF over the corresponding hemisphere yields the and for normal incidence. The reduced scattering coefficient , where is the asymmetric parameter in the Henyey–Greenstein phase function, was determined by fitting (i.e., the ratio of and ) of the thin-film samples with those calculated from the adding-doubling method. The ranges of , , and were also estimated by analyzing the measurement results. The polarization-dependent BRDF and BTDF were measured at oblique incidence to demonstrate the effect of surface scattering. Moreover, a monochromator and an integrating sphere were employed to measure and in the wavelength region from .
The five sintered Zenith PTFE samples were purchased from SphereOptics (Concord, New Hampshire). The thicknesses of samples 1 to 5 were measured with a micrometer to be , , , , and , respectively. Table 1 lists the sample parameters as well as some of the measured optical properties. The densities of samples range from . These samples were cut into pieces from a sheet, except for sample 5, which came as a piece from the manufacturer.
Parameters of the PTFE samples and the directional-hemispherical reflectance and transmittance obtained by integrating the measured BRDF and BTDF at λ=635nm .
|Sample||Thickness d (mm)||Density ρ (g∕cm3)||Reflectance R||Transmittance T||R+T||R∕T|
A three-axis automated scatterometer (TAAS)23 with a laser diode was used for the measurement of the BRDF and BTDF. The BRDF is defined as the reflected radiance divided by the incident irradiance,is the reflected radiance, is the incident radiance, is the angle of incidence, and is the solid angle of the incident light. The measurement equation of BRDF is given as and are respectively the incident and reflected optical powers reaching the detector, is the observation angle (i.e., polar angle of the detector), and is the solid angle of the detector. An instrument calibration factor is introduced and is determined later. The BTDF is defined in a similar way by replacing with (i.e., transmitted radiance) in Eq. 1. A collimation lens was placed in front of the laser diode, resulting in a highly collimated light with a beam divergence less than . A lock-in amplifier was used to produce a modulated voltage signal for the laser diode controller; thus the effect of stray light can be minimized. The samples were mounted in a rotary stage so that the incidence angle could be changed. A detector controlled by another rotary stage measured the scattered light in an observation angle ranging from . In front of the detector, there was a precision aperture with a diameter of . The distance between the aperture and the sample holder was measured to be . Consequently, the detector solid angle was and the half cone angle was approximately . Since the detector blocked the incident light, the BRDF within of the retroreflection direction could not be measured.23
Figure 1 shows the schematic of light scattering from a semitransparent film. The observation angle is defined separately for BRDF and BTDF as depicted in the figure. For BTDF measurement when , both the transmitted collimated light (i.e., unscattered light) and the scattered light may contribute to the detector signal. The measured BTDF should exhibit a noticeable peak at when the collimated light transmittance is significant enough.
The wavelength dependence of and was characterized with a monochromator and an integrating sphere at wavelengths from .24 The light source was a tungsten-halogen lamp. The rotation of two gratings inside the monochromator and the change of filters at the inlet of the monochromator are capable of achieving a wavelength resolution of . The rms fluctuation of power from the monochromator was estimated to be 1%. The light exiting the monochromator is directed to the sample mounted either before the sphere entrance port for measuring or at the back of the sphere wall for measuring . A chopper is used to obtain a phase-locked optical signal. A silicon detector was located at the bottom port with a baffle to prevent the direct illumination of the detector from the scattered light. The detector output is amplified by a transimpedance preamplifier before being sent to a lock-in amplifier. The integrating sphere has an inner wall of in diameter. The entrance port is in diameter. The inner wall is coated with PTFE with a reflectance between 0.98 and 0.99 in the visible spectral range.
Theory and Calculations
With the assumptions that the medium does not emit light (i.e., cold medium), the wave-like interactions are negligible, the medium is homogeneous, and the polarization state is neglected, the light scattering in a PTFE film can be modeled using the RTE asis the (spectral) radiance, is the physical distance that the light travels, defines the scattering albedo, and represent the light propagation directions, is the solid angle, and is the scattering phase function. In many applications, the Henyey–Greenstein function is used as the scattering phase function and is given by20 is the angle between the propagation directions of the incident light and the scattered light (i.e., ), and is called the asymmetric parameter (ranging from to 1). The Henyey–Greenstein function yields isotropic scattering when . If , there are more forward scattered photons than backward scattered photons. If approaches 1, all photons are scattered in the direction parallel to the incident light. However, if , more photons are scattered backward than forward. If approaches , all photons will be scattered in the opposite direction of the incident light. The RTE governs the changes of radiance along an infinitesimal distance inside the medium. As a beam of light travels a distance in the medium, the radiance is attenuated due to scattering and absorption. Conversely, the radiance is enhanced by the incoming light scattered from other directions. Due to the complexity of the integrodifferential equation, no effective method is available that can analytically solve the RTE without some sort of approximation.
If only the hemispherical properties are of interest, the adding-doubling method is an effective technique that results in high accuracy for obtaining and of an absorbing and scattering medium. The adding-doubling method was intro-duced by van de Hulst25 to solve the RTE in a parallel slab composed of multiple layers. In the adding method, if the reflection and transmission functions of each individual layer are known, the reflection and transmission functions of the composite can be calculated. The adding method is termed as the doubling method when the layers are identical in thickness and in all the RTE parameters.20 The reflection and transmission functions for an arbitrarily thick slab can be obtained by repeatedly adding and doubling the layers until the desired thickness is reached. Subsequently, the directional-hemispherical reflectance and transmittance can be calculated by integrating the reflection and transmission functions. Detailed discussion of the adding-doubling method can be found from Refs. 16, 20. The adding-doubling method requires inputs of the scattering albedo, optical thickness with being the film thickness, and the asymmetric parameter to predict the directional-hemispherical properties of PTFE films.
A Monte Carlo algorithm was developed to model the BRDF and BTDF by considering volume scattering but neglecting the effect of surface scattering. Furthermore, the Monte Carlo simulation does not consider the dependent multiple scattering and coherent wave-like interactions. After a photon bundle is released from air to the PTFE medium, the propagation step between subsequent scattering events is calculated by , where is a random number between 0 and 1 and is generated based on a quasirandom sequence.26 The previous equation results in the mean free path (i.e., average propagation step) of . The angle between the direction vectors before and after the scattering is determined by the Henyey–Greenstein scattering phase function given in Eq. 4 following the procedure described by Wang, Jacques, and Zheng.27 After each scattering event, the energy of the photon bundle reduces to due to the absorption, where is the energy of the photon bundle before scattering.
In the Monte Carlo simulation, is calculated as the summation of the reflected photon bundle energy divided by the total energy of the incident photon bundle, which is the product of the number of photon bundles and the energy assigned to each incident photon bundle. Similarly, is calculated as the summation of the transmitted photon bundle energy divided by the total energy of the incident photon bundles. The calculated and using the Monte Carlo simulation with a photon bundle number of is in agreement with that using the adding-doubling method within a statistical fluctuation of 0.5%. To determine the BRDF and BTDF, a virtual detector with a solid angle of is placed in intervals in the corresponding hemisphere. For the calculation of BRDF and BTDF, photon bundles are used and each run takes approximately of CPU time with a Pentium 4 processor when the thickness of the sample is . The resulting fluctuation is within 1.5% of the BRDF and BTDF values. Additional photon bundles can be used to reduce the statistical fluctuation with longer calculation time.
Measured Bidirectional Reflectance Distribution Function and Bidirectional Transmittance Distribution Function
Figure 2 shows the measured BRDF for all five samples and the BTDF for samples 1 to 4 for normal incidence at the wavelength . The transmittance of sample 5 is less than 0.01, and its BTDF cannot be measured due to the limited signal-to-noise ratio. The measured BRDF and BTDF exhibit nearly diffuse characteristics. For a perfectly diffuse reflector, the BRDF should be a constant of . In a round-robin test of the PTFE material, Early 28 also reported the similar variation of the BRDF of thick PTFE slabs. As the sample thickness decreases, the corresponding BRDF value decreases because of the increased transmission. Interestingly, the hump in the BRDF near becomes more obvious as the thickness decreases. Due to the measurement uncertainty of the laser scatterometer, the measured BRDF does not exhibit perfect symmetry with respect to .
As expected, the BTDF values increase as the sample thickness decreases. The measurements reveal that the BTDF of samples 2, 3, and 4 are flatter than the corresponding BRDF even though the values are smaller. This can be explained by the multiple scattering of light inside the film. In order for the photon bundles to transmit through the PTFE films, they need to experience numerous scattering events, resulting in random propagation in the medium. Therefore, the BTDF of sample 4 with thickness of is nearly uniform within observation angles.
To obtain the hemispherical properties, the BRDF and BTDF data need to be carefully analyzed and fitted to perform integration. Additional measurements with closely spaced data intervals were performed. As an example, Fig. 3 shows the average of the measured BTDF (circles) and the fitted (solid line) BTDF of sample 1 at normal incidence. Eight measurements were conducted at different locations on the sample. The location on the sample was changed by rotating it by after each measurement. Because the laser beam was not aligned to the center of rotation of the sample, the location of the laser spot changed when the sample was rotated. Error bars stand for the random uncertainty at the 95% confidence level. The measured results beyond are not reliable because of the alignment difficulty. Since the measurement exhibits some fluctuations due to the measurement uncertainty and sample inhomogeneity, a polynomial equation was used to fit the BTDF data. For BTDF of sample 1, a quadratic function is sufficient. The standard error of estimate of the fitted equation is , which is much smaller than the standard deviation of the data points. Notice that an extrapolation of the fitted equation for from is performed before the integration; however, the relative error caused by this extrapolation is estimated to be less than 0.5%. In addition, the bidirectional properties are assumed to be symmetric with respect to .
The polynomial fitting functions are integrated to obtain from the BRDF and from the BTDF of each sample. When is taken as 1, the summation of and is always greater than unity by approximately 6%. While the measurements of TAAS agree well with specular samples (within 2% for a smooth Si), it is about 5% higher for a rough silicon surface and about 5% higher for the -thick PTFE than those obtained from a standard reference instrument at NIST.23, 28 This may be caused by the detector collection geometry and the uncertainty of the aperture area. To reduce the bias uncertainty, all TAAS measurement results are reduced by 6%. This may be considered as a self-calibration procedure. With , the reflectance of sample 5 obtained by integration of BRDF is 0.988, which agrees well with the reported values.8 The summation of and for each of the thinner samples is about unity, as shown in Table 1. When the sample location is fixed, the variation of repeated BRDF and BTDF measurements is within 2%. The combined uncertainty of TAAS is within 5% with a confidence level of 95%. For BTDF of sample 1, however, the data variation at different locations on the sample is greater than the instrument uncertainty due to sample inhomogeneity.
The effects of incidence angle on the BRDF are shown in Fig. 4 for all samples. The measurements were performed at incidence angle , 30, 50, and with a linear polarizer. The results are plotted in terms of to make the comparison easier. Even at oblique incidence, the diffuse scattering feature is preserved in a large observation angle range except a peak in the forward scattering direction. The peak is much higher for -polarization than for -polarization, suggesting that surface scattering becomes important at large angles of incidence. This can be understood by two facts: 1. the surface looks smoother at glazing angles; and 2. the reflectance at the interface between air and PTFE is higher for -polarization than for -polarization. The Brewster angle is about for refractive index . At the Brewster angle, the surface reflectance for -polarization approaches zero. When the data are examined more carefully, the BRDF for is higher than that for at large observation angles. Hence, volume scattering must also play a role in the enhanced forward scattering. It can also be seen that for samples 4 and 5, the peaks at and are much lower than for samples 1 to 3. Visual observation reveals that sample 1, 2, and 3 may contain some surface crystallization that has improved the smoothness and specularity of these samples. This may be the reason for the specular peaks in samples 1 to 3 for the -polarization at and . At , samples 4 and 5 also exhibit specular peaks for -polarization but with a smaller peak than that of the thinner samples.
The increased BRDF at oblique incidence results in a reduction in the BTDF as depicted in Fig. 5 , where the ordinate is in terms of . It should be noted that the scale of the ordinates is different for each sample. The data points at and were missing for due to blocking of the beam by the rotary stage. As discussed earlier, the BRDF decreases with increasing thickness, and also the polarization has little effect on the BRDF at normal incidence. It can be clearly seen that more reduction of BTDF existed for -polarization than for -polarization. On the other hand, because the transmitted light experiences multiple scattering events, they are redirected into random directions. Hence, the BTDF shape appears symmetric about , regardless of the incidence angle.
Estimation of and the Ranges of , , and
If the absorption coefficient is negligible, the directional-hemispherical properties of a volume scattering medium depend only on the film thickness and the reduced scattering coefficient .20 In the present study, is estimated by comparison of the calculated -to- ratio for the given thickness with those obtained by integrating the measured BRDF and BTDF over the corresponding hemisphere. In fact, the ratio is independent of the bias uncertainty in the scatterometer measurement. Note that the refractive index of PTFE is around 1.36, and the reflectance at the interface between air and PTFE is approximately 2.3% near normal incidence. The reflected light will be distributed hemispherically due to surface roughness and may exhibit a broadened specular peak. In the present study, the adding-doubling method is employed to solve the RTE for the PTFE films, without considering surface scattering due to refractive index mismatch and surface roughness.
To determine , the ratio is calculated as a function of for the thickness equal to the thicknesses of the four thinner samples. By assuming that all samples have the same , the standard error of estimate between the calculated and measured for all four samples is a function of only. The standard error of estimate between the measured and the calculated values reaches a minimum at . Considering the measurement uncertainty and sample-to-sample variation, the reduced scattering coefficient is estimated to be . Here, the obtained value is approximately one order of magnitude higher than the reported value in Ref. 22.
Figure 6 shows the measured ratio and the calculated values using the adding-doubling method with . The inset shows the results for smaller values to compare the measurement with the calculation clearly. While in general the reflectance and transmittance are complicated functions of the film thickness, the ratio exhibits a linear dependence on the film thickness for from . A linear regression of the adding-doubling calculations yields that . Using different values of , it was found that a constant of 0.59 remains the same, regardless of the value of .
For strongly scattering and nonabsorbing media, the linear dependence of on can be derived from other approximation models such as the Kubelka–Munk model (or two-flux model)29 and the three-flux model.30 Star, Marijnissen, and van Gemert31 reported the transformation between the scattering coefficient used in the Kubelka–Munk model and the scattering coefficient of RTE. By solving the governing equations of these models, it can be shown that the ratio of is in linear relationship with , similar to that calculated from the adding-doubling method. The linear dependence of on film thickness for strongly scattering and nonabsorbing mediums can be obtained from all these models, including the adding-doubling method, Monte Carlo simulation, Kubelka–Munk model, and the three-flux model. Although there are some variations in the final analytical expressions, depending on different model approximations, the linear relationship provides a very simple approach for the determination of the reduced scattering coefficient of PTFE films by measuring and of a sample with known thickness.
Table 2 lists the calculated and of the -thick PTFE slab using the adding-doubling method. The parameters are set to be and (i.e., ). As mentioned in the previous section, the reduced scattering coefficient is the dominant factor in determining the hemispherical properties. Hence, the calculation shown in Table 2 is also applicable for different values as long as is fixed. The absorption coefficient must be less than for of the -thick sample to be greater than 0.977. It should be noted that the measured values are usually greater than 0.98.8
The influence of the absorption coefficient on the directional-hemispherical reflectance and transmittance of the 10-mm -thick PTFE slab when σλ=1670cm−1 and g=0.9 .
As can be seen from Table 2, the adding-doubling calculation suggests that the -thick slab can have a nonzero transmittance as large as 0.01 when the absorption is neglected. To assess the transmittance value, the integrating sphere is used with the laser diode as the light source. The sample is placed at the entrance port. The transmittance of sample 5 ( thick) is estimated to be , which indicates that the absorption coefficient should indeed be less than (referring to Table 2). The uncertainty in the integrating sphere measurement for the -thick sample is large due to the low signal-to-noise ratio. The uncertainty associated with the integrating sphere measurements for most samples is greater than that in the scatterometer measurement. Therefore, the integration of BRDF and BTDF allows a more accurate determination of and at the laser wavelength. Because the obtained is at least four orders of magnitude smaller than the scattering coefficient, the absorption in PTFE films is neglected hereafter. It should be mentioned that the absorption coefficient of PTFE is much smaller than typical biological tissues.
If the collimated light transmittance is significant, the scattering coefficient can easily be determined by using Beer’s law . For highly scattering samples, however, the optical signal received by the detector when contains both the collimated light transmitted through the sample without scattering and the scattered light. The correct way to determine is to subtract the power received by the detector at for normal incidence to that at , where is a small angle from that allows the detector to stay away from direct light exposure. Unfortunately, even for the thinnest sample in this study, no peaks in the BTDF can be observed around . The BTDF of sample 1 ( thick) for eight different measurements at normal incidence are plotted in Fig. 7 for with a interval. The fitted curve and the error bounds with 95% confidence interval are also shown. The resulting BTDF is rather flat in this region and no peak can be seen at , suggesting that the transmitted collimated light is overwhelmed by the scattered light in the direction parallel to the incidence. In other words, the contribution of the collimated light transmission to the BTDF must be less than the uncertainty, which is about at the 95% confidence level. Therefore, the maximum direct transmittance (transmitted power divided by the incident laser power) should be . Setting yields , which may be regarded as the lower limit of the scattering coefficient. If the scattering coefficient were less than , a distinguishable peak, exceeding the upper bound of the uncertainty, would have appeared in the BTDF at .
Under the condition that the wave-like interactions and the dependent scattering are negligible, the mean free path should be greater than several wavelengths. Suppose that the minimum mean free path is , which is approximately four times the laser wavelength (i.e., ), then the upper limit of the scattering coefficient is estimated to be . Based on and , the asymmetric parameter should be between 0.86 and 0.96, suggesting that PTFE is a strongly forward scattering material. Compared with the parameters of biological tissues, the scattering coefficient of PTFE is approximately ten times larger.17 Much thinner samples, e.g., , would be necessary to precisely determine the scattering coefficient of PTFE.
Monte Carlo Simulation
The purpose of using Monte Carlo simulation is to see how volume scattering affects the BRDF and BTDF, without considering surface scattering and polarization effects. In the calculation, the parameters are set to be and , i.e., . The simulation results of BRDF and BTDF at normal incidence are shown in Fig. 8 . Compared with the measurements shown in Fig. 2, the Monte Carlo simulation can predict the general trend of the BRDF but with a rapid drop at large observation angles. The Monte Carlo results do not show the hump near , because the hump is associated with the surface scattering effect. Note that the humps are more prominent for samples 1 to 3 than for samples 4 and 5. It should be noted that the Monte Carlo model predicted the same and as the adding-doubling method did. For the BTDF, the Monte Carlo simulation also captures the essential features qualitatively; but the resulting BTDF exhibits a much larger angular variation compared with the experimental results. There are two reasons for this discrepancy. First, if surface scattering were included, there will be a total internal reflection for light emerged from PTFE to air at polar angles (inside PTFE) greater than the critical angle . This critical angle may largely affect the energy distribution at large observation angles. Second, the Henyey–Greenstein scattering phase function is an oversimplified one, and its applicability for predicting the bidirectional scattering properties may be questionable.
Figure 9 shows the cosine modified BRDF and BTDF for sample 2 at oblique incidence. Even without considering polarization dependence and surface scattering, the BRDFs exhibit asymmetric characteristics due to the enhanced forward scattering and reduced backward scattering at increase incidence angles, although the enhancement and reduction are overpredicted in the Monte Carlo simulation. Furthermore, the Monte Carlo simulation cannot capture the specular peaks for -polarization as observed experimentally. For BTDF, the Monte Carlo simulation predicts the symmetric feature with respect to , but overpredicts the reduction of BTDF at large incidence angles. Future research is needed to develop a more sophisticated Monte Carlo model to include surface reflection and refraction by considering the surface roughness and polarization dependence.
Spectral Reflectance and Transmittance Measurements
The directional-hemispherical reflectance and transmittance of the PTFE films were measured using the system of monochromator and integrating sphere as discussed previously. Care must be taken for the transmittance measurement. The reference signal is obtained by sending the light directly to the back of the sphere. On the other hand, when the sample covers the front entrance port, the transmitted light is diffuse. For the reflectance measurement, sample 5 is used as the reference and it is placed on the back port of the sphere facing the laser. The other samples are interchanged to obtain the ratio of the reflected signal. It is assumed that the reflectance of the -thick PTFE sample is 0.988, as determined by integration of the BRDF. The results of and obtained by integrating BRDF and BTDF are shown as square marks. The measured and using the diode laser at with the sphere are shown with the diamond marks. All the measured results agree very well. The relative uncertainty of the reflectance and transmittance measurements was estimated to be 10% with a confidence level of 95%.
Figures 10a and 10b show the measured spectral reflectance and transmittance in the wavelength region from of samples 1 to 4. As the thickness decreases, the transmittance increases but the reflectance decreases. As the wavelength increases, the reflectance decreases but the transmittance increases. This suggests that the scattering is stronger toward short wavelengths, as expected for small particle scattering. The reduced scattering coefficient can be calculated using the formula for each sample at each measurement wavelength. The average for the four samples and its standard deviation are plotted in Fig. 10c. It can be clearly seen that the reduced scattering coefficient increases as the wavelength is reduced. Furthermore, the value obtained previously with the scatterometer measurement of at is slightly lower than that obtained with the integrating sphere; but the agreement is within the expanded uncertainty.
A comprehensive experimental study is performed for a set of sintered PTFE films that are semitransparent and scatter light volumetrically. The BRDF and BTDF are obtained at the wavelength of with a laser scatterometer at various angles of incidence for both polarizations. Integrating the BRDF and BTDF in the hemisphere results in directional-hemispherical reflectance and transmittance, respectively. Analyzing the reflectance-to-transmittance ratio allows the determination of a reduced scattering coefficient . The scattering coefficient cannot be easily determined because the collimated light transmission is negligible, even for the thinnest sample studied here. An effort is made to estimate the range of to be between 1200 and . The absorption coefficient is estimated to be less than . Furthermore, the effect of wavelength is investigated measuring and in the region from .
Contrary to previously reported scattering coefficient and absorption coefficients of PTFE being similar to those of biological tissues, this study reveals that the scattering coefficient of PTFE is nearly ten times greater than that of typical tissues, while the absorption coefficient of PTFE is much less. The present study calls for careful distinction between directly transmitted light and scattered light toward the direction parallel to the incidence. This is important for future research of light scattering in biological media for disease diagnostics and laser therapy.
A Monte Carlo simulation is performed without considering surface scattering and polarization effects. Some of the features in the BRDF and BTDF cannot be described by the simple Monte Carlo model, suggesting the need to develop more comprehensive theoretical models that can include scattering by rough surfaces as well as volume scattering. The polarization-dependent BRDF and BTDF data presented in this work will be useful for model validation. The microstructure of PTFE depends on crystallization, particle cluster size, porosity, density, etc. Further study is also needed to investigate the effect of microstructures on the interaction of light with PTFE and other scattering materials.
The authors thank S. A. Prahl at Oregon Health and Science University for the communication on the adding-doubling software. Authors Lee and Zhang acknowledge the support of the National Science Foundation (CBET-0500113), and Zhang also acknowledges the support of the Department of Energy (DE-FG02-06ER46343) and the Optical Technology Division of NIST. Certain commerical equipment, instruments, or materials are identified in this work to foster understanding. Such identification does not imply recommendation or endorsement by the National Institute of Standards and Technology, nor does it imply that the materials or equipment identified are necessarily the best available for the purpose.