Estimating scattering of pure water from density fluctuation of the refractive index

The use of density derivative of the refractive index from the classic Lorentz-Lorenz equation or its variations performed poorly in estimating the scattering by water, leading to the alternative use of pressure derivative instead, which however has been scarcely measured due to its extremely low sensitivity. Recently, density derivative has been deduced directly from theoretical models. Three characterizations of density derivative of the refractive index were evaluated and scattering of water thus calculated converge with each other within 3.5% and agree with the measurement by Morel (Cahiers Oceanographiques, 20, 157, 1968) within 2% (with depolarization ratio = 0.039), all improving over the earlier estimates based on either density or pressure derivatives. Taking into account of uncertainty associated with the depolarization ratio, the prediction based on the model by Proutiere et al. (J. Phy. Chem., 96, 3485, 1992) still agrees with the measurement within the experimental errors (2%). 2009 Optical Society of America OCIS codes: (290.5840) Scattering, molecules; (010.4450) Oceanic optics References and links 1. A. Morel, "Optical properties of pure water and pure sea water," in Optical Aspects of Oceanography, N. G. Jerlov, and E. S. Nielsen, eds., (Academic Press, New York, 1974), pp. 1-24. 2. K. S. Shifrin, Physical Optics of Ocean Water (American Institute of Physics, New York, 1988). 3. I. L. Fabelinskii, Molecular Scattering of Light (Plenum Press, New York, 1968). 4. A. Morel, "Etude Experimentale de la diffusion de la lumiere par l'eau, les solutions de chlorure de sodium et l'eau de mer optiquement pures," Journal de Chimie Physique 10, 1359-1366 (1966). 5. A. Morel, "Note au sujet des constantes de diffusion de la lumiere pour l'eau et l'eau de mer optiquement pures," Cahiers Oceanographiques 20, 157-162 (1968). 6. M. Kerker, The Scattering of Light and Other Electromagnetic Radiation (Academic Press, New York, 1969). 7. H. Buiteveld, J. H. M. Hakvoort, and M. Donze, "The optical properties of pure water," Proc. SPIE 2258, 174-183 (1994). 8. M. S. Twardowski, H. Claustre, S. A. Freeman, D. Stramski, and Y. Huot, "Optical backscattering properties of the "clearest" natural waters," Biogeosciences 4, 1041-1058 (2007). 9. A. Morel, B. Gentili, H. Claustre, M. Babin, A. Bricaud, J. Ras, and F. Tieche, "Optical properties of the "clearest" natural waters," Limnol. Oceanogr. 52, 217-229 (2007). 10. E. W. Washburn, International Critical Tables of Numerical Data, Physics, Chemistry and Technology (McGraw Book Co., New York, 1930). 11. A. M. Evtyushenkov,and Y. F. Kiyachenko, "Determination of the dependence of liquid refractive index on pressure and temperature," Opt. Spectrosc. 52, 95-98 (1982). 12. R. M. Waxler, C. E. Weir, and H. W. Schamp, "Effect of pressure and temperature upon the optical dispersion of Benzene, carbon tetrachloride and water," J. Res. Natl. Bur. Stand. 5, 489-498 (1964). 13. R. W. Austin and G. Halikas, "The index of refraction of seawater," (Scripps Institute of Oceanography, La Jolla, 1974), p. 121. 14. Z. Niedrich, "Dispersion interactions and the refractive index of liquids," Physica 128B, 69-75 (1985). 15. H. Eisenberg, "Equation for the Refractive Index of Water," J. Chem. Phys. 43, 3887-3892 (1965). 16. A. Proutiere, E. Megnassan, and H. Hucteau, "Refractive index and density variations in pure liquids: A new theoretical relation," J. Phys. Chem. 96, 3485-3489 (1992). 17. A. Bot, "Comment on "Refractive index variations in pure liquids. A new theorectical relation," J. Phys. Chem. 97, 2804 (1993). #104801 $15.00 USD Received 4 Dec 2008; revised 16 Jan 2009; accepted 21 Jan 2009; published 27 Jan 2009 (C) 2009 OSA 2 February 2009 / Vol. 17, No. 3 / OPTICS EXPRESS 1671 18. L. W. Tilton and J. K. Taylor, "Refractive index and dispersion of distilled water for visible radiation at temperature 0 to 60 C," J. Res. Natl. Bur. Stand. 20, 419 (1938). 19. M. Jonasz and G. R. Fournier, Light Scattering by Particles in Water: Theoretical and Experimental Foundations (Academic Press, New York, 2007). 20. R. S. Farinato and R. L. Rowell, "New values of the light scattering depolarization and anisotropy of water," J. Chem. Phys. 65, 593-595 (1976). 21. G. S. Kell, "Isothermal compressibility of liquid water at 1 atm," J. Chem. Eng. Data 15, 119-122 (1970). 22. X. Quan and E. S. Fry, "Empirical equation for the index of refraction of seawater," Appl. Opt. 34, 34773480 (1995). 23. P. D. T. Huibers, "Models for the wavelength dependence of the index of refraction of water," Appl. Opt. 36, 3785-3787 (1997). 24. P. E. Ciddor, "Refractive index of air: new equation for the visible and near infrared," Appl. Opt. 35, 15661573 (1996). 25. G. Cohen and H. Eisenberg, "Light scattering of water, deuterium oxide, and other pure liquids," J. Chem. Phys. 43, 3881-3887 (1965).


Introduction
According to Einstein-Smoluchowski theory, scattering of light in a particle-free medium is due to fluctuations in the dielectric constant (ε), caused by the random motion of molecules.Following this theory, the volume scattering function at 90°, β(90), for a pure liquid is modeled as [1,2]: where 2 ε ∆ represents the fluctuations of dielectric constant within a volume of V ∆ , which is small as compared to the wavelength λ and yet large enough for the laws of statistical thermodynamics to apply.The Cabannes factor 6 6 6 7 f δ δ , where δ is the depolarization ratio due to the fluctuation in the orientation of anisotropic molecules, such as H 2 O.For pure liquid, ε ∆ (hence n ∆ , where n is the refractive index and ε = n 2 ) can be expanded as a function of two independent thermodynamic variables, density ρ and temperature T [3], The second term in Eq. ( 2), the fluctuation of refractive index due to temperature, is typically < 1% of the first term, the fluctuation due to density [2,3].Omitting the second term in Eq.
(2) and using the thermodynamic statistics, (90) ( ) ( ) 2 where k is the Boltzmann constant, β T is the isothermal compressibility, and T is absolute temperature.
Scattering by pure water has been measured experimentally by Morel [4,5] at five wavelengths of 366, 405, 436, 546, and 578 nm with a relative experimental error of 2% (his results are listed in row 1, Table 1).Comparing to the earlier experimental determinations, the measurements by Morel gave the smallest values (Table 1 of [1]).For the theoretical calculation of molecular scattering by pure water, various representations of 2 n ρ ρ ∂ ∂ have been tried but none of them agreed with the experiment [1,2].Kerker [6] compared the scattering measurements for a variety of liquids with the corresponding estimates calculated ∂ ∂ , either derived theoretically or empirically, and concluded that none of them gave satisfactory agreement for all liquids.As an example, the estimates of β(90) using the density derivative that are deduced from two classical equations are shown in Table 1 (row 2 for the Lorentz-Lorenz equation and row 3 for the Laplace equation), each with an average difference of 32% and -18%, respectively.Table 1.β(90) measured by Morel (row 1) and its comparison with various theoretical calculations with δ = 0.039 at 20°C.In the last column, the numbers in the parenthesis are the results if δ = 0.051 or 0.09 were used in the calculations.
The relative experimental error of the measurements [4].
Because of these discrepancy, the density derivative was replaced with the pressure derivative, i.e., , which can be measured relatively easier.Alternatively, Equation ( 4) has been used by Morel [1], Shifrin [2], and Buiteveld et al. [7] to estimate scattering by pure water.Twardowski et al [8] reviewed these early studies and recommended the use of the estimate by Buiteveld et al [7], which was based on the most recent experimental results and showed a better agreement with the measurements by Morel [1].Their estimate was also used by Morel et al. [9] in their recent study of Southern Pacific Ocean.The results by Buiteveld et al, estimated as having a relative error of 6%, are shown in row 4 of Table 1.The percentage difference between Morel's measurement and their estimates is ~3.2%.Note, Buiteveld et al used a value of 0.051 for δ, which would give a percentage difference of ~5.9%, and in Table 1 we used a value of 0.039, which has led to a better agreement.
Typical values of ( ) T n P ∂ ∂ for water are ~1.5×10-10 Pa -1 varying about ±10% with the temperature.Given that the typical experimental precision in measuring n is on the order of 10 -5 , P ∆ needs to span at least ~ 1 Atm, which we believe is beyond the microscopic fluctuation range for pressure.Among the few experiments determining this pressure derivative, P ∆ ranged from 1.5 Atm [10,11] to ~100 Atm [12].Austin and Halikas [13] commented that among all the measurements, n(P) is of the least quality as compared to those of n(T), n(λ), and n(S) (S is salinity), and it is well known that a derivative such as n P ∂ ∂ is very sensitive to the errors in the function of n.
In the meantime, several theoretical relationships between refractive index and density for liquid water have been developed and verified with experimental measurements [14][15][16].Given the inherent uncertainties associated with determining ( ) T n P ∂ ∂ experimentally, it is of great interest to re-evaluating the scattering of pure water directly from its underlying physics: density fluctuation of the refractive index.

Density fluctuation of the refractive index
From the spherical cavity model, Proutiere et al. [16] where ε 0 is the dielectric constant for vacuum, N the number of molecules per unit volume, α the molecular polarizability, v the volume of the molecular cavity, and the bar indicates a mean value within V ∆ .The authors referred Eq. ( 5) as the generalized formulation for the Lorentz-Lorenz equation.Assuming that can be obtained, where , M is the molar mass and N A the Avogadro number.Proutiere et al. [16] pointed out that since the classic equations (such as Lorentz-Lorenz or Laplace) are only approximations, they may fail to give accurate values for the derivatives, such as 2 n ρ ρ ∂ ∂ .By imposing the same approximations but only after differentiating n 2 in Eq. ( 5), they derived, The comparison with the measurements of 75 different liquids including water showed that Eq. ( 7) has an average error of 4.1% [16].Note, Proutiere et al. [16] only derived the isobaric term, ∂ ∂ [17].We will refer Eq. ( 7) as PMH model.Based on Van der Waals's equation of state and assuming a liquid be represented by uniformly distributed, statistically permanent pairs of spherical molecules, Niedrich [14] deduced the general formulae for both We will refer Eq. ( 8) as Niedrich model.By examining the data of Tilton and Taylor of the refractive of index of pure water [18], which, with its 7 significant decimal digits, remains by far the most accurate measurement, Eisenberg [15] found that the right hand side of Eq. ( 6), which was supposed to be independent of T and P, actually decreases continuously with the increase in temperature of water.This finding further confirmed that the Lorentz-Lorenz equation cannot be used directly in calculating water scattering.He proposed an improved Lorentz-Lorenz equation, where A, B, and C are coefficients dependent on the wavelength only.Equation ( 9) is valid for water between 0° and 60°C at any wavelength in the visible and agrees with the Tilton and Taylor's measurements [18] to the 7 th decimal digit.Both coefficients B and C, though derived empirically, have thermodynamic meaning and are related to the thermal expansion coefficient at constant pressure.From Eq. ( 9), we have, Except for a correction parameter B, Eq. ( 10) is the same as that derived directly from the Lorentz-Lorenz equation (row 2, Table 1).From the same relationship as Eq. ( 9), Jonasz and Fournier [Eq.2.68 in 19], however, derived a different representation of , which is in the form similar to the Laplace equation (row 3, Table 1).Equation ( 9) is based on the measurements taken under one atmospheric pressure, but the test with the by Waxler et al [12] showed that it applies to the pressure up to 1100 bar.We will refer to Eq. ( 10) as Eisenberg model.

Results and discussion
In their respective theoretical estimates of β(90), Morel [1] and Shifrin [2] use a value of 0.09 for the depolarization ratio, δ, and Buiteveld et al. [7] adopted a value of 0.051.Recently, Jonasz and Fournier [19] recommended a value of 0.039 determined by Farinato and Rowell [20] and showed that a better agreement with the measurements by Morel [4] could be achieved if the lower value of δ was used in Buiteveld et al. (0.039 vs. 0.051).Farinato and Rowell [20] attributed the lowered value of δ to the correction of the detector's finite acceptance angle effect and to the reduction of stray light contamination.By progressively reducing the stray light from using no filter, to a medium band filter (bandwidth = 22.5 nm) and finally to a narrow band filter (bandwidth = 0.46 nm), their experiment led to a series of determinations of δ with decreasing values of 0.051, 0.045, and 0.039, respectively.They also reported that similar trends of improvement had been observed for other molecules, such as gaseous methane and liquid benzene.
We would assume that 0.039 is the most accurate for δ by far because the narrow band filter would remove most (if not completely) of the stray light.However, there is still uncertainty in this parameter; and we will use the other two values for δ to show its effect.Everything else the same, using 0.039 for δ instead of 0.051 or 0.09 would translate into a reduction in the values of β(90) by a factor of 0.97 or 0.89, respectively.
For isothermal compressibility, β T , we used the equation reported in Kell [21], where β T in bar -1 , T c in Degree Celsius, and the values for a 0-5 and b 0-1 are 50.88630,0.7171582, 0.7819867×10 -3 , 31.62214×10 -6 , -0.1323594×10 -6 , 0.6345750×10 -9 , 1.0, 21.65928×10 -3 , respectively.Equation ( 11) is valid for T c from 0 to 110°C at one atmospheric pressure with a standard error of 0.002×10 -6 bar -1 .Note, Eq. ( 11) is effectively the same as the one used in Buiteveld et al. [7], who generated a quadratic fit to the data produced by Eq. (11).For n, we used the empirical equation by Quan and Fry [22], who fitted the experimental data selected by Austin and Halikas [13] with an average error of 1.5×10 -5 .Their equation for pure water is: The coefficients have the following values: n 0 = 1.31405, n 4 = -2.02×10 - , n 5 = 15.868,n 7 = -0.00423,n 8 = -4382, n 9 = 1.1455×10 6 .Their model, originally developed from the visible region, fits the available data well over an extended range covering the UV to the near-IR (200 -1100 nm) [23].The index of refraction of water in Austin and Halikas [13] was defined relative to the air, and so is Eq. ( 12).The n in Eqs. ( 7), ( 8) and ( 10) is the index in vacuum, which can be derived by multiplying Eq. ( 12) with the refractive index of air [24], where ν is the wavenumber (reciprocal of the vacuum wavelength) in µm -1 , and the other coefficients have the following values: k 0 = 238.0185,k 1 = 5792105, k 2 = 57.362,and k 3 = 167917 µm -2 .Equation ( 13) is for standard air at 15°C, 1 Atm, and 0% humidity.The refractive index of air does change with the temperature, pressure and water content.The variations, however, would affect β(90) only on the 8 th decimal digit, and therefore can be safely ignored for the purpose of this study.We calculated β(90) using Eq. ( 3) with 2 ( ) T n ρ ρ ∂ ∂ estimated using the three models, Eqs. ( 7), (8), and ( 10), respectively and with δ = 0.039.The results are shown in Table 1 (rows 5, 6, and 7 respectively).Apparently, in terms of comparison with Morel's measurements, the latest development in characterizing density fluctuation of the refractive index has led to an improved estimate of scattering by water molecules.The average differences with the measurements (-0.67%, 0.72%, and 1.89%, for Eqs. ( 7), (8), and (10), respectively) are all within the relative experimental error (2%) and much smaller than those for the earlier estimates using the same equation (Eq.3).They also perform better than the model by Buiteveld et al., whose results were based on a different equation (Eq.4).In particular, estimates based on PMH model are about 0.5-1.5% lower than Morel's experimental values, which might be expected if there were extremely low levels of unavoidable contamination in Morel's sample, which in turn would cause additional scattering.
We also calculated β(90) with δ = 0.051 and 0.09 and their comparisons with the measurements are shown in the last column (inside the parentheses) of Table 1.Except for the estimates based on the Laplace equation (row 3), δ = 0.09 leads to significantly larger differences between various models and the measurements, which might suggest that this value be too big for the depolarization ratio of water molecules.With δ = 0.051, although all the three new models performed better than the earlier models, only the estimates based on PMH model agree with Morel's measurements (1.97%) within the experimental error (2%).Based on Table 1 alone, however, we still cannot completely reject 0.051 as a possible value of δ.Regardless of the values used for δ, the two theoretical models (PMH and Niedrich) performed better than the empirical one (Eisenberg).
Spectral β(90) from 350 to 700 nm calculated by Buiteveld et al. [7], Morel [Table 4 of Ref. 1], Niedrich model, and Eisenberg model relative to that by PMH model are plotted in Fig. 1 (a) for δ = 0.039 and (b) for δ = 0.051.Also shown in Fig. 1 are the measured β(90) by Morel [4,5].It can be seen from Fig. 1, estimates based on the three models agree to each other within 3.5%, with better agreement in the shorter wavelengths.It can also be seen from Fig. 1(a), where δ = 0.039 was used, that the values of the measurements by Morel fall between the estimates based on the two theoretical models, PMH and Niedrich.
As we have mentioned above, PMH model and the Lorentz-Lorenz equation share the same theoretical basis with difference in the stages where the approximations are applied.Eisenberg model is basically the same as the Lorentz-Lorenz equation but with an empirical adjustment based on the measurements.Niedrich model was derived from the Van de Waals equation of state.Yet, Table 1 and Fig. 1 show that they all have converged on estimating the molecular scattering by water, indicating a closure is achieved among different models in characterizing the density fluctuation of water, which in turn leads to a better closure between the theory and the observation.
Theoretically, Eqs. ( 3) and ( 4 3) can be expressed as a simple function of n, whose values for water have been measured with a relatively high precision (typically 10 -5 but up to 10 -7 as in [18]).Based on Eqs.(7), (8), and (10), the relative errors in 2 ( ) T n ρ ρ ∂ ∂ due to the uncertainty in n are at most 0.001%.Therefore the total error is dominated by the modeling uncertainty of the density derivative.The convergence of the estimates of water scattering based on the three models as shown in Table 1 and Fig. 1 suggested that the uncertainty associated with the modeling of the density derivative has been significantly reduced.
On the other hand, the accuracy of ( ) T n P ∂ ∂ is limited by the experiment.With no analytical form existing, ( ) T n P ∂ ∂ can only be approximated as n P ∆ ∆ .Among the few experiments, P ∆ were set at least ~ 1-2 Atm in order for n ∆ to be detectable (typical value for ( ) T n P ∂ ∂ is only ~ 1.5×10 -10 Pa -1 ).Even though linearity has been assumed, the function n(P) does behave in a nonlinear manner over a large pressure range [13].In principal, we still do not know ( ) T n P ∂ ∂ at the microscopic scales, under which the Einstein-Smoluchowski theory applies.
Among the three models of the density derivative, we recommend the PMH model because 1) the prediction based on it agrees closest to the measurement within the experimental error, even considering uncertainty in δ, ranging from -0.67% for δ = 0.039 to 1.97% for δ = 0.051; and 2) its underlying theoretical basis can also lead to the classic Lorentz-Lorenz equation.The MATLAB code to compute these coefficients will be provided upon request or can be downloaded from ftp://ftp.umac.org/zhang/betaw_ZH2009.m The molecular scattering of water predicted with PMH model together with Eq. ( 3) show that the spectral dependence of scattering as in 0 ( ) s λ λ has a slope s of value -4.28 when the anchor wavelength λ 0 is at 450 nm.Fitting to Morel's experimental data (row 1 of Table 1) rendered an s of -4.32 with λ 0 at 436 nm, the same value was also reported by Morel [1] based on his theoretical calculation.Fitting to Buiteveld et al. [7] values gave an s of -4.14, the same as that estimated by Twardowski et al. [8].Shifrin [2] reported a value of -4.17.Based on our results, both estimates by Buiteveld et al. and Shifrin seemed a little bit low (in absolute values).The minimum scattering occurs at ~26°C, which is close to the minimum of 22°C determined both theoretically and experimentally by Cohen and Eisenberg [25].Note, Buiteveld et al. reported a maximum at 15°C, which is probably due to the relatively large uncertainty in the pressure derivative of the refractive index they used.The overall temperature variation is relatively small, about 3.7% for T from 0 to 26°C.Scattering decreases slowly as the pressure increases, by ~1.3% for an increase in P of 100 bar.

Conclusions
With the recent theoretical development in characterizing density derivative of the refractive index for liquid, we re-estimated the scattering by pure water using three different models (PMH, Niedrich, and Eisenberg).These theoretical estimates agree to each other within 3.5% and with the experiment [5] within 2% (δ = 0.039).Use of the PMH model would predict values about 0.5-1.5% lower than the measurements, which is expected because even extremely low level of unavoidable dust contamination in the sample would have caused additional scattering.Taking into account of uncertainty associated with the depolarization ratio, the prediction based on the PMH model still agrees with the measurement within 2%.Because of extremely low sensitivity in n P ∆ ∆ as an approximation of the pressure derivative of the refractive index, the use of density derivative offers numerical advantage.
their theory, it can be shown that the isothermal term is equal to the isobaric term, i.e.,

#
104801 -$15.00USD Received 4 Dec 2008; revised 16 Jan 2009; accepted 21 Jan 2009; published 27 Jan 2009 (C) 2009 OSA ) are the same; therefore the uncertainty associated with the density fluctuation term determines which one to use in practice.The density derivative (

Fig 1 .
Fig 1.Comparison of pure water scattering β(90) from the theories and the experiment.The percentage difference relative to the values of calculation using Eq. 3 and PMH model (Eq.7) are plotted.The value used for the depolarization ratio δ was 0.039 in (a) and 0.051 in (b).Values for Eisenberg and Niedrich are calculated using Eq. 3, and the Eisenberg (Eq.10) and Niedrich (Eq.8) models, respectively.Other values are from Buiteveld et al. [7], Morel [4, 5] for experiment, and Morel [table 4 of 1] for his theoretical calculation.Note Morel and Buiteveld et al. used δ = 0.09 and 0.051 in their respective calculations, here we scaled their values correspondingly.