Next Article in Journal
Luminosity Measurements at the LHC at CERN Using Medipix, Timepix and Timepix3 Devices
Previous Article in Journal
Cosmological Formation of (2 + 1)-Dimensional Soliton Structures in Models Possessing Potentials with Local Peaks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Critical Angle Refractometry for Lossy Media with a Priori Known Extinction Coefficient

by
Spyridon Koutsoumpos
,
Panagiotis Giannios
and
Konstantinos Moutzouris
*
Laboratory of Electronic Devices and Materials, Department of Electrical & Electronic Engineering, University of West Attica, 12244 Egaleo, Greece
*
Author to whom correspondence should be addressed.
Current address: Barcelona Institute of Science and Technology, Institute for Research in Biomedicine (IRB), 08028 Barcelona, Spain.
Submission received: 4 June 2021 / Revised: 14 July 2021 / Accepted: 22 July 2021 / Published: 3 August 2021
(This article belongs to the Section Applied Physics)

Abstract

:
Critical angle refractometry is an established technique for determining the refractive index of liquids and solids. For transparent samples, the critical angle refractometry precision is limited by incidence angle resolution. For lossy samples, the precision is also affected by reflectance measurement error. In the present study, it is demonstarted that reflectance error can be practically eliminated, provided that the sample’s extinction coefficient is a priori known with sufficient accuracy (typically, better than 5 % ) through an independent measurement. Then, critical angle refractometry can be as precise with lossy media as with transparent ones.

1. Introduction

Critical angle refractometry is the standard for determining the refractive index of transparent media [1,2,3,4,5,6]. The method relies on measurements of the reflectance, R ( θ ) , of the interface, formed by a transparent front medium (which is commonly a prism) of known refractive index, n o , and the sample of unknown index, n, for a range of incidence angles, θ , that includes the critical angle, θ c , of total internal reflection (TIR) [7,8,9].
As it is shown in Figure 1, the critical angle is located at a sharp discontinuity of the θ -derivative, R ( θ ) ; the refractive index of the sample is obtained from the TIR condition:
n / n o = sin θ c ,
with the corresponding relative error:
δ n / n = δ θ c / tan θ c .
In determining the TIR critical angle, the error, δ θ c , is mainly regulated by the angular resolution of the experimental setup, which is typically between 20 μ rad and 100 μ rad for state-of-the-art refractometers [10].
Critical angle refractometry has also been used with lossy media where the refractive index turns into the complex number, n = n r + i n i , the imaginary part of which (imaginary index) incorporates loss effects, namely, absorption and/or scattering. The reflectance profile, R ( θ ) , becomes smoother, and the critical angle, θ c , is located at the maximum of R ( θ ) , which now marks the gradual transition from the attenuated TIR (a-TIR) to the partial reflection regime; see Figure 1. The underlying assumption is that Equation (1) is still valid (at least approximately) and that it can be used to obtain an estimate of the real index of the lossy sample [11,12]. This assumption introduces systematic errors that have been the subject of extensive discussion; see, e.g., [13,14,15] Let us note that fitting the reflectance profile, R ( θ ) , to Fresnel equations is another way to compute the complex optical constants [16,17,18]. However, bearing its own strengths and weaknesses, data regression is not a point method and, thus, cannot be taken as a “critical angle” refractometry approach in itself.
Recently, the universal condition of a-TIR,
n r / n o = g 1 ( θ c , R c ) and n i / n o = g 2 ( θ c , R c ) ,
was proposed, where R c is the critical reflectance at θ c , and functions g 1 and g 2 are derived from the Fresnel equations [19]. In the case of the s-polarisation, which is the main foucus of the present study, these equations are given in Section 2. Using Equation (3), the complex refractive index can be determined by θ c and R c .
Then, the relative error of the real index, δ n r / n r , depends on both δ θ c and δ R c :
δ n r n r = 1 g 1 · g 1 θ c · δ θ c 2 + g 1 R c · δ R c 2 .
In determining the critical reflectance, the error δ R c is mainly due to intensity fluctuations of the laser sources and, therefore, to be of the order of 10 2 , which is at least two orders of magnitude larger than δ θ c . As a result, critical angle refractometry is less precise with lossy media than it is with transparent ones.
In this paper, it is demonstrated that the critical reflectance can be numerically computed with higher accuracy than it can be experimentally obtained, provided that the extinction coefficient of the sample is a priori known with sufficient accuracy (<5%, typically). Then, the real index relative error for lossy samples (being δ θ c and δ R c dependent; see Equation (4)) can become equal to the relative error for transparent samples (which is only δ θ c dependent; see Equation (2)). This is a significant advancement in the field of optical characterisation of absorbing and/or scattering media, which include, among others, various forms of biological matter, non-transparent liquids, colloids and food products.

2. Background Theory and Initial Observations

The functions g 1 and g 2 [19] are compactly expressed in terms of the sample’s complex dielectric constant, ϵ r + i · ϵ i = ( n r + i · n i ) 2 , and the prism’s dielectric constant, ϵ o = n p 2 :
g 1 = ϵ r 2 + ϵ i 2 2 ϵ o + ϵ r 2 ϵ o 1 / 2 and g 2 = ϵ r 2 + ϵ i 2 2 ϵ o ϵ r 2 ϵ o 1 / 2 ,
where
ϵ r ϵ o = α + t 2 1 + t 2 and ϵ i ϵ o = γ 2 α 2 1 + t 2 ,
with
α = ( 1 + γ ) 2 2 ρ 2 γ , γ = 2 t ( 3 ρ 2 2 ρ 2 ) t + ρ ( 9 ρ 2 12 ρ 8 ) t 2 4 ,
and
t = tan θ c , ρ = 1 + R c 1 R c .
As it is described in [19], the preceding formalism, which is valid for s-polarisation, calculates output values ( n r , n i ) from refractometric input pairs ( θ c , R c ) by successive algebraic substitutions into Equations (5)–(8). Let us remark that Equations (5) and (6) are also valid for p-polarised light, in which case however, there exist no explicit solutions for α and γ such as those provided by Equation (7). Instead, these parameters can be calculated as solutions to two algebraic equations, as is analytically shown in [19]. In that sense, the forthcoming analysis may be replicated for p-polarisation. Despite some more complication in the calculations, the results, reported here for s-polarisation, are qualitatively similar to those of p-polarisation.
Figure 2 gives a graphical insight by the isoangular curve which consists of pairs ( n r / n o , n i / n o ) , computed by keeping θ c constant ( θ c = π / 4 in this example), while letting R c to vary. There is no loss of generality, associated with the specific choice of θ c , because larger (smaller) critical angles simply shift the curve to the right (left). An indicative measurement (with θ c = π / 4 , R c varying) defines a single point in the ( n r / n o , n i / n o ) -plane, whose coordinates uniquely determine the complex refractive index.
In Figure 2, one can observe quite wide range, spanned by the imaginary index as the critical reflectance varies, especially compared to the relatively modest shift for the real index, which is nevertheless non-negligible. The vertical dashed line in Figure 2 marks the transparency limit n r / n o = sin ( π / 4 ) 0.707 , which is the estimate of the TIR-condition (Equation (1)) for each point of the isoangular curve. Unsurprisingly, this estimate practically coincides with the exact result for small n i / n o , e.g., n i / n o 10 4 , a range that includes, for example, most liquids at wavelengths far from their ultraviolet or infrared resonances, but excludes most forms of biological matter, even in the therapeutic optical window. The most interesting feature in Figure 2 is the turning point of the isoangular curve, which vividly demonstrates that the systematic error, introduced by use of Equation (1) with lossy media, does not increase monotonically with attenuation but decreases above the turning point, and even becomes zero when the isoangular curve crosses the transparency limit marked by the dashed vertical line. This behaviour, which now emerges as a natural consequence of the universal a-TIR condition, had been previously observed [13,20] and labeled as inexplicable by other authors [14].

3. Method’s Application with a Priori Known Extinction Coefficient

3.1. Main Concept

The sample’s extinction coefficient, μ , i.e., the sum of the absorption coefficient, μ a , and the scattering coefficient, μ s , can be a priori known via an independent experimental method such as absorption spectroscopy or collimated transmittance spectroscopy. This type of independent extinction coefficient measurement has been reported earlier for pure water [21], various emulsions [22], bioliquids [23], oils [24] and semiconductors [25]; typically, the relative error, σ μ = δ μ / μ , in the measurement of the extinction coefficient is of 0.1 % to 10 % [26,27,28,29,30].
The formalism, described in Section 2, can be used to calculate the real part of the refractive index (real index) from the value of the critical angle alone. The direct measurement of critical reflectance is no longer needed, since R c can be determined from the value of the extinction coefficient, which is related to the imaginary index [31]:
μ = 4 π n i λ ,
so that the following equation holds true ( λ is the light wavelength):
g 2 ( θ c , R c ) = n i n o = μ λ 4 π n o .
Solving Equation (10) for R c (the only unknown variable) is a simple task for iterative computing software such as Mathematica’s FindRoot or Mathcad’s Find. As long as the initial guess root is kept within its physically meaningful range, 0 < R c < 1 , the solution is always unique and unambiguous. Substituting the measured value of θ c and the numerically retrieved value of R c in g 1 ( θ c , R c ) yields the real index of the lossy sample.
The procedure described has a straightforward graphical interpretation. Let us consider the measurement of a critical angle, for example, θ c = π / 4 , so that one can refer to Figure 2. An a priori known value of the extinction coefficient (equivalently, n i ) defines a unique horizontal line that crosses the isoangular curve at a single point; locating its position is equivalent to numerically computing the value of R c . The abscissa of that point is the real index of the lossy medium.

3.2. Critical Reflectance Error

The formalism of Section 2 enables the calculation of n r from the values of θ c and R c . The method, introduced in Section 3.1, retrieves the numerical value of R c from the extinction coefficient (when the extinction coefficient is already known) instead of its direct measurement. To better appreciate the advantage of the proposed method, one has to remember that the experimental error of a direct reflectance measurement is typically greater than 10 2 . Breaking this threshold does not seem possible with standard unstabilized laser sources of moderate cost, such as those used in common refractometers.
Let us now show that the propagated error, δ R c , of the numerically retrieved R c can be kept below the 10 2 threshold. To this end, one accounts for the fact that R c is numerically computed from input values of θ c and μ , or equivalently, θ c and n i , which are assumed to be independent variables, so that the covariance between them can be taken to be zero. Therefore, propagated error δ R c conveys corresponding experimental errors δ θ c and δ n i :
δ R c = R c θ c · δ θ c 2 + R c n i · δ n i 2 .
Critical reflectance, R c , is not expressed explicitly as a function of θ c and n i . Instead, Equation (10) can be restated in the implicit form: F ( θ c , R c , n i ) = F ( g 2 ( θ c , R c ) , n i ) = g 2 ( θ c , R c ) n i / n o = 0 . Then, the derivatives in Equation (11) can be obtained using the implicit function theorem:
R c θ c = F θ c / F R c = g 2 θ c / g 2 R c ,
and
R c n i = F n i / F R c = 1 n o / g 2 R c .
Moreover, one can introduce the relative error, σ μ , which is related to δ n i and δ μ :
σ μ = δ μ μ = δ n i n i = δ n i g 2 · n o .
Substituting Equations (12)–(14) into Equation (11) yields:
δ R c = | g 2 R c | 1 g 2 θ c · δ θ c 2 + ( σ μ · g 2 ) 2 .
Values of g 2 and its two partial derivatives that appear in Equation (15) can be easily computed for any pair ( θ c , R c ) . Therefore, error δ R c can be estimated as a function of the variables θ c , R c , δ θ c , and σ μ . Numerical investigation reveals that δ R c is practically independent of δ θ c , at least when this parameter is kept within its reasonable range, 20 μ rad δ θ c 100 μ rad. This observation reflects the easily verifiable fact that
g 2 θ c · δ θ c σ μ · g 2 ,
which reduces Equation (15) to the approximate form:
δ R c | g 2 R c | 1 · σ μ · g 2 ,
and narrows δ R c down to a function of only three variables, namely, of θ c , R c , and σ μ .
Figure 3 shows the results of calculations of δ R c versus σ μ for indicative values θ c = π / 4 and R c = 0.5 : the solid line, produced via the exact Equation (15). In agreement with the approximate Equation (17), δ R c is not affected by δ θ c and increases linearly with σ μ , which is taken within the range 0.1 % σ μ 10 % . One can observe that δ R c 10 4 (or 10 3 ) when σ μ 0.1 % ( 1 % ), indicating that critical reflectance error can be two (one) orders of magnitude smaller than the 10 2 threshold, when R c is numerically retrieved from the a priori known extinction coefficient.

3.3. Real Index Error

The relative error, δ n r / n r , can be computed as a function of the variables θ c , R c , δ θ c , and σ μ by substituting Equation (15) into Equation (4).
Figure 3 shows δ n r / n r versus σ μ for θ c = π / 4 and R c = 0.5 : the dashed lines correspond to δ θ c = 100 , 50 and 20 μ rad, as indicated. The error lines become horizontal for small σ μ at a constant value δ n r / n r δ θ c / tan θ c (in this example, tan θ c = 1 ). In the horizontal region, the relative real index error is practically equal to that of a transparent sample (cf. Equation (2)). This observation alone proves that refractometry with lossy media can be as precise as the refractometry with transparent ones, provided that the extinction coefficient is a priori known with sufficient accuracy. The horizontal regions terminate at a maximum permissible σ μ , above which δ n r / n r continues to increase beyond the transparency limit. In the example considered, this transition occurs in the vicinity of σ μ 1 % when δ θ c = 20   μ rad, at σ μ 2 % when δ θ c = 50   μ rad, and at σ μ 4 % when δ θ c = 100   μ rad.
Note that the upper and lower dashed lines in Figure 3 flatten out at δ n r / n r 10 4 (when δ θ c = 100   μ rad) and δ n r / n r 2 × 10 5 (when δ θ c = 20   μ rad). These levels match the specifications of modern refractometers operating with transparent samples, which typically prescribe an error δ n / n between ∼ 10 4 (standard precision) and ∼ 2 × 10 5 (ultimate precision) [32].
As a final exercise, let us consider the refractometric measurement of a lossy sample with n r = 1.5 by use of a reference front medium with n o = 2 , so that n r / n o = 0.75 . With these assumptions, as given, it is straightforward to calculate the magnitude of σ μ that is required to reach (i) the ultimate precision δ n r / n r 2 × 10 5 , assuming that the refractometer measures critical angles with an ultrahigh accuracy δ θ c = 20   μ rad, and (ii) a standard precision δ n r / n r 10 4 , with a more tolerant δ θ c = 100   μ rad. The calculations are shown in Figure 4 by solid and dashed lines, respectively.
Therein, one observes a general trend: the maximum permissible σ μ decreases with increasing n i . This trend clearly indicates that the need for accurate measurements of the extinction coefficient is more acute when optical loss is growing. There is one exception to this, which manifests itself as a characteristic sharp peak at larger n i . This peak reflects the vertical slope of the isoangular curve at the turning point, see Figure 2. Quite naturally, σ μ is found increasing again at the transparent limit, where the slope of the isoangular curve is again nearly vertical. Likewise, σ μ decreases rapidly for larger n i after the peak because the slope of isoangular curve above there is zero
To better understand the interpretation made, let us point that a vertical slope in Figure 2 suggests that small enough change in the imaginary index causes zero shift to the real index. Similarly, a horizontal slope suggests that small change in the imaginary index causes a maximum shift in the real index. Hence, the real index error minimizes (and the maximum permissible σ μ maximizes) as soon as the slope turns vertical, and vice versa.
Meantime, the most important observation in Figure 4 is that the real index of the sample can be measured with an extreme precision, δ n r / n r 2 × 10 5 , as long as n i < 10 1 and σ μ 0.1 % (solid line). Note that although being strict, this situation is indeed realistic. For example, the extinction coefficient of water has been measured with a relative error even less than 0.1 % near its infrared absorption peak (at λ = 1410 nm), where n i = 2.3675 × 10 4  [21]; see discussion in Section 3.4 below. This is also possible with a much more lax σ μ 1 % for all n i < 10 2 . An adequate for many applications error δ n r / n r 10 4 is reached with σ μ 4 % for all n i < 4 × 10 2 (dashed line). These observations illustrate the functionality of the proposed method with all but the most extremely attenuating media.

3.4. Comments on Implementation Issues

The proposed method involves the tandem use of two instruments. The first determines the reflectance profile at the sample’s interface with a prism; prism coupling refractometers operate routinely in various laboratories and they are also commercially available; see, e.g. Metricons’s 2010/M model. The second measures the sample’s extinction coefficient. Collimated transmittance setups (typically, homemade) suit ideally this purpose; see, e.g., [21,26,27]. The corresponding operating principle is straightforward. The transmitted portion of a collimated light beam, travelling through sample sections of variable length, is monitored; the extinction coefficient is then deduced by fitting experimental results to the known Lambert-Beer law. This technique is commonly used for fluids where liquid cells of different sizes facilitate tuning of the path length. The technique is also applicable for solid samples, provided that slices of variable thickness can be cut.
This tandem measurement is actually known. In [21], the (real) refractive index and the extinction coefficient of several types of liquids, including water, were measured by the combined use of critical angle refractometetry and collimated transmittance spectroscopy. Unaware of the attenuated-TIR condition (Equation (3)) at that time, the standard TIR law of transparent media (Equation (1)) was exploited in [21] to determine the real index. Let us see, through a didactic example, how the method, proposed in this paper, adjusts those results.
Let us start by selecting the reported data for water at a wavelength λ = 1410 nm, where the extinction coefficient, μ = 21.1 cm 1 , was measured [21] translating via Equation (9) into the imaginary index, n i = 2.3675 × 10 4 . The relative error in this measurement was as small as σ μ = 0.05 % . Using an SF13 glass prism with n o = 1.7068 , the real index was approximated by the standard TIR condition, leading to the value n r = 1.3183 , which, given the index of the prism, corresponds to the critical angle θ c = 50.568 .
Using the required input pair ( θ c = 50.568 , μ = 21.1 cm 1 ), the method, proposed in the current paper, yields, first, the critical reflectance R c = 91.99 % and, then, the real index of water, corrected to a new value, n r = 1.31842 . The correction from the previous estimate ( n r = 1.3183 ) is relatively modest, while non negligible for applications requiring high accuracy. The modest correction effect is expected (cf. Figure 2), since the imaginary index in our example just exceeds 10 4 ; it would have been much stronger, had the sample been more attenuating.
Using Equation (15), the prescribed σ μ = 0.05 % value, combined with an assumed error in the measurement of critical angle, 20 μ rad δ θ c 100   μ rad, leads to a propagated error in the numerical estimation of the critical reflectance within δ R c 0.002 % . Such an accuracy is surely orders of magnitude higher than the ∼ 1 % uncertainty with which critical reflectance can be directly measured. This finding manifests the main advantage of the method here described. Then, within the same δ θ c range, this method yields the real index with an error of 2 × 10 5 δ n r 1 × 10 4 , the latter calculated as it is described in Section 3.3. It is worth noticing the compliance of this result with Figure 3.

4. Conclusions

In this paper, the universal attenuated total internal reflection (a-TIR) condition is exploited in order to accurately determine the real part of the refractive index (real index), nr, of lossy media from the critical incidence angle, θ c , at s-polarisation given the extinction coefficient, μ . This is accomplished in a two-step procedure. First, the critical reflectance, R c , is numerically retrieved from the pair ( θ c , μ ) . Then, the real index is calculated from θ c and R c .
Numerical investigation reveals that R c can be recovered (in the first step) more accurately than it can be directly measured. As a consequence the determination of the real index (in the second step) becomes more precise. The relative error, δ n r / n r , can be reduced down to 2 × 10 5 when the respective error in the independent measurement of μ is in the range of 0.1 % to 1 % .
The results, obtained here, demonstrate that refractometry of lossy media can be as precise as with transparent media, a development that is of interest to various applications in biomedical optics, material characterisation, analytical chemistry and quality control.

Author Contributions

Conceptualization, K.M.; methodology, S.K., K.M.; validation, P.G., S.K.; investigation, S.K., K.M.; writing—original draft preparation, S.K., P.G.; writing—review and editing, K.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data supporting this article are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Castrejón-Pita, J.R.; Morales, A.; Castrejón-García, R. Critical angle laser refractometer. Rev. Sci. Instrum. 2006, 77, 035101. [Google Scholar] [CrossRef]
  2. Moutzouris, K.; Hloupis, G.; Stavrakas, I.; Triantis, D.; Chou, M.H. Temperature-dependent visible to near-infrared optical properties of 8 mol% Mg-doped lithium tantalate. Opt. Mater. Express 2011, 1, 458. [Google Scholar] [CrossRef]
  3. Moutzouris, K.; Stavrakas, I.; Triantis, D.; Enculescu, M. Temperature-dependent refractive index of potassium acid phthalate (KAP) in the visible and near-infrared. Opt. Mater. 2011, 33, 812–816. [Google Scholar] [CrossRef]
  4. Moutzouris, K.; Papamichael, M.; Betsis, S.C.; Stavrakas, I.; Hloupis, G.; Triantis, D. Refractive, dispersive and thermo-optic properties of twelve organic solvents in the visible and near-infrared. Appl. Phys. B 2013, 116, 617–622. [Google Scholar] [CrossRef]
  5. Dongare, M.; Buchade, P.; Shaligram, A. Refractive index based optical Brix measurement technique with equilateral angle prism for sugar and Allied Industries. Optik 2015, 126, 2383–2385. [Google Scholar] [CrossRef]
  6. Chiappe, C.; Margari, P.; Mezzetta, A.; Pomelli, C.S.; Koutsoumpos, S.; Papamichael, M.; Giannios, P.; Moutzouris, K. Temperature effects on the viscosity and the wavelength-dependent refractive index of imidazolium-based ionic liquids with a phosphorus-containing anion. Phys. Chem. Chem. Phys. 2017, 19, 8201–8209. [Google Scholar] [CrossRef] [PubMed]
  7. Chen, J.Y.; Xie, Z.H.; Li, W.N.; Lin, S.B.; Zhang, L.L.; Liu, C.X. Construction and investigation of a planar waveguide in photo-thermal-refractive glass by proton implantation. Optik 2020, 207, 164461. [Google Scholar] [CrossRef]
  8. Revathi, V.; Rajendran, V. Investigation about nonlinear optics and antibacterial activity of pyrrolidine-2-carboxylic acid cadmium chloride hydrate single crystal. Optik 2018, 154, 234–241. [Google Scholar] [CrossRef]
  9. Liu, C.X.; Shen, X.L.; Guo, H.T.; Li, W.N.; Wei, W. Proton-implanted optical waveguides fabricated in Er3+-doped phosphate glasses. Optik 2017, 131, 132–137. [Google Scholar] [CrossRef]
  10. Koutsoumpos, S.; Giannios, P.; Moutzouris, K. Extended derivative method of critical-angle refractometry for attenuating media: Error analysis. Meas. Sci. Technol. 2021, 32, 105007. [Google Scholar] [CrossRef]
  11. Zeng, H.; Wang, J.; Ye, Q.; Deng, Z.; Mei, J.; Zhou, W.; Zhang, C.; Tian, J. Study on the refractive index matching effect of ultrasound on optical clearing of bio-tissues based on the derivative total reflection method. Biomed. Opt. Express 2014, 5, 3482. [Google Scholar] [CrossRef] [Green Version]
  12. Sun, T.Q.; Ye, Q.; Wang, X.W.; Wang, J.; Deng, Z.C.; Mei, J.C.; Zhou, W.Y.; Zhang, C.P.; Tian, J.G. Scanning focused refractive-index microscopy. Sci. Rep. 2014, 4, 5647. [Google Scholar] [CrossRef] [Green Version]
  13. Meeten, G.H.; North, A.N.; Willmouth, F.M. Errors in critical-angle measurement of refractive index of optically absorbing materials. J. Phys. E Sci. Instrum. 1984, 17, 642–643. [Google Scholar] [CrossRef]
  14. Goyal, K.G.; Dong, M.L.; Kane, D.G.; Makkar, S.S.; Worth, B.W.; Bali, L.M.; Bali, S. Note: Refractive index sensing of turbid media by differentiation of the reflectance profile: Does error-correction work? Rev. Sci. Instrum. 2012, 83, 086107. [Google Scholar] [CrossRef] [Green Version]
  15. Morales-Luna, G.; García-Valenzuela, A. Viability and fundamental limits of critical-angle refractometry of turbid colloids. Meas. Sci. Technol. 2017, 28, 125203. [Google Scholar] [CrossRef] [Green Version]
  16. Giannios, P.; Toutouzas, K.G.; Matiatou, M.; Stasinos, K.; Konstadoulakis, M.M.; Zografos, G.C.; Moutzouris, K. Visible to near-infrared refractive properties of freshly-excised human-liver tissues: Marking hepatic malignancies. Sci. Rep. 2016, 6, 27910. [Google Scholar] [CrossRef] [Green Version]
  17. Giannios, P.; Koutsoumpos, S.; Toutouzas, K.G.; Matiatou, M.; Zografos, G.C.; Moutzouris, K. Complex refractive index of normal and malignant human colorectal tissue in the visible and near-infrared. J. Biophotonics 2016, 10, 303–310. [Google Scholar] [CrossRef] [PubMed]
  18. Räty, J.; Pääkkönen, P.; Peiponen, K.E. Assessment of wavelength dependent complex refractive index of strongly light absorbing liquids. Opt. Express 2012, 20, 2835. [Google Scholar] [CrossRef] [PubMed]
  19. Koutsoumpos, S.; Giannios, P.; Stavrakas, I.; Moutzouris, K. The derivative method of critical-angle refractometry for attenuating media. J. Opt. 2020, 22, 075601. [Google Scholar] [CrossRef]
  20. Meeten, G.H. Refractive index errors in the critical-angle and the Brewster-angle methods applied to absorbing and heterogeneous materials. Meas. Sci. Technol. 1997, 8, 728–733. [Google Scholar] [CrossRef]
  21. Kedenburg, S.; Vieweg, M.; Gissibl, T.; Giessen, H. Linear refractive index and absorption measurements of nonlinear optical liquids in the visible and near-infrared spectral region. Opt. Mater. Express 2012, 2, 1588. [Google Scholar] [CrossRef]
  22. Ninni, P.D.; Martelli, F.; Zaccanti, G. Intralipid: Towards a diffusive reference standard for optical tissue phantoms. Phys. Med. Biol. 2010, 56, N21–N28. [Google Scholar] [CrossRef] [PubMed]
  23. Zhao, Y.; Qiu, L.; Sun, Y.; Huang, C.; Li, T. Optimal hemoglobin extinction coefficient data set for near-infrared spectroscopy. Biomed. Opt. Express 2017, 8, 5151. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Zhao, D.; Zhang, G.; Zhang, X.; Li, D. Optical properties of paraffin at temperature range from 40 to 80 C. Optik 2018, 157, 184–189. [Google Scholar] [CrossRef] [Green Version]
  25. Forouhi, A.R.; Bloomer, I. Optical properties of crystalline semiconductors and dielectrics. Phys. Rev. B 1988, 38, 1865–1874. [Google Scholar] [CrossRef]
  26. Li, X.; Zhao, J.M.; Wang, C.C.; Liu, L.H. Improved transmission method for measuring the optical extinction coefficient of micro/nano particle suspensions. Appl. Opt. 2016, 55, 8171. [Google Scholar] [CrossRef] [Green Version]
  27. Ninni, P.D.; Martelli, F.; Zaccanti, G. Effect of dependent scattering on the optical properties of Intralipid tissue phantoms. Biomed. Opt. Express 2011, 2, 2265. [Google Scholar] [CrossRef] [Green Version]
  28. Feder, I.; Duadi, H.; Fixler, D. Single wavelength measurements of absorption coefficients based on iso-pathlength point. Biomed. Opt. Express 2020, 11, 5760. [Google Scholar] [CrossRef]
  29. Marchesini, R.; Bertoni, A.; Andreola, S.; Melloni, E.; Sichirollo, A.E. Extinction and absorption coefficients and scattering phase functions of human tissues in vitro. Appl. Opt. 1989, 28, 2318. [Google Scholar] [CrossRef] [PubMed]
  30. Ogusu, K.; Suzuki, K.; Nishio, H. Simple and accurate measurement of the absorption coefficient of an absorbing plate by use of the Brewster angle. Opt. Lett. 2006, 31, 909. [Google Scholar] [CrossRef]
  31. Wu, S.; Lian, J.; Song, P.; Gao, S.; Wang, X.; Ma, Z.; Wang, Y.; Guan, W. Optical properties of Sr3NbGa3Si2O14 crystal. Opt. Int. J. Light Electron Opt. 2013, 124, 686–688. [Google Scholar] [CrossRef]
  32. Zhou, W.; Zhou, Y.; Albert, J. A true fiber optic refractometer. Laser Photonics Rev. 2017, 11, 1600157. [Google Scholar] [CrossRef]
Figure 1. (a) Fresnel reflectance, R ( θ ) , profile, as a function of the incidence angle, θ , at an interface between a transparent sample and a transparent prism, such that the ratio of unknown refractive index to the known index is n / n o = 0.75 . Long-dashed arrow indicates the location of the corresponding critical angle, θ c . (b) Fresnel reflectance profile at an interface between a lossy sample and a transparent prism, such that n r / n o = 0.75 and n i / n o = 0.005 , where n r and n i are the real and inginary parts of the refractive index, respectively. Short-dashed arrows indicate the location of the corresponding critical angle and critical reflectance, R c . (c) The critical angle rests at the peak of the reflectance derivative, which is also depicted (in arbitrary units).
Figure 1. (a) Fresnel reflectance, R ( θ ) , profile, as a function of the incidence angle, θ , at an interface between a transparent sample and a transparent prism, such that the ratio of unknown refractive index to the known index is n / n o = 0.75 . Long-dashed arrow indicates the location of the corresponding critical angle, θ c . (b) Fresnel reflectance profile at an interface between a lossy sample and a transparent prism, such that n r / n o = 0.75 and n i / n o = 0.005 , where n r and n i are the real and inginary parts of the refractive index, respectively. Short-dashed arrows indicate the location of the corresponding critical angle and critical reflectance, R c . (c) The critical angle rests at the peak of the reflectance derivative, which is also depicted (in arbitrary units).
Physics 03 00036 g001
Figure 2. Isoangular curves in the ( n r / n o , n i / n o ) -plane for the s-polarisation, θ c = π / 4 and variable R c . As one moves along the direction of the arrow, R c increases in steps of 0.005 , within the range [ 0.385 , 0.990 ] . The vertical dashed line marks the transparency limit, n r / n o = sin ( π / 4 ) .
Figure 2. Isoangular curves in the ( n r / n o , n i / n o ) -plane for the s-polarisation, θ c = π / 4 and variable R c . As one moves along the direction of the arrow, R c increases in steps of 0.005 , within the range [ 0.385 , 0.990 ] . The vertical dashed line marks the transparency limit, n r / n o = sin ( π / 4 ) .
Physics 03 00036 g002
Figure 3. Solid line: error in the numerical determination of the critical reflectance, δ R c , as a function of the exstinction coefficient relative error, σ μ . Calculations are based on Equation (15) for θ c = π / 4 and R c = 0.5 . Dashed lines: relative error, δ n r / n r , as a function of σ μ for the same pair of θ c and R c values; δ θ c varies from 100 via 50 to 20 μ rad, as indicated. Calculations result from substituting Equation (15) into Equation (4).
Figure 3. Solid line: error in the numerical determination of the critical reflectance, δ R c , as a function of the exstinction coefficient relative error, σ μ . Calculations are based on Equation (15) for θ c = π / 4 and R c = 0.5 . Dashed lines: relative error, δ n r / n r , as a function of σ μ for the same pair of θ c and R c values; δ θ c varies from 100 via 50 to 20 μ rad, as indicated. Calculations result from substituting Equation (15) into Equation (4).
Physics 03 00036 g003
Figure 4. Maximum permissible relative error, σ μ , as a function of n i for ultimate refractometric precision δ n r / n r = 2 × 10 5 , when δ θ c = 20   μ rad (solid line), and standard refractometric precision δ n r / n r = 10 4 , when δ θ c = 100   μ rad (dashed line). A sample with n r = 1.5 and a prism with n o = 2 are considered.
Figure 4. Maximum permissible relative error, σ μ , as a function of n i for ultimate refractometric precision δ n r / n r = 2 × 10 5 , when δ θ c = 20   μ rad (solid line), and standard refractometric precision δ n r / n r = 10 4 , when δ θ c = 100   μ rad (dashed line). A sample with n r = 1.5 and a prism with n o = 2 are considered.
Physics 03 00036 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Koutsoumpos, S.; Giannios, P.; Moutzouris, K. Critical Angle Refractometry for Lossy Media with a Priori Known Extinction Coefficient. Physics 2021, 3, 569-578. https://0-doi-org.brum.beds.ac.uk/10.3390/physics3030036

AMA Style

Koutsoumpos S, Giannios P, Moutzouris K. Critical Angle Refractometry for Lossy Media with a Priori Known Extinction Coefficient. Physics. 2021; 3(3):569-578. https://0-doi-org.brum.beds.ac.uk/10.3390/physics3030036

Chicago/Turabian Style

Koutsoumpos, Spyridon, Panagiotis Giannios, and Konstantinos Moutzouris. 2021. "Critical Angle Refractometry for Lossy Media with a Priori Known Extinction Coefficient" Physics 3, no. 3: 569-578. https://0-doi-org.brum.beds.ac.uk/10.3390/physics3030036

Article Metrics

Back to TopTop