Lateral light scattering in fibrous media

Lateral light scattering in fibrous media is investigated by computing the modulation transfer function (MTF) of 22 paper samples using a Monte Carlo model. The simulation tool uses phase functions from infinitely long homogenous cylinders and the directional inhomogeneity of paper is achieved by aligning the cylinders in the plane. The inverse frequency at half maximum of the MTF is compared to both measurements and previous simulations with isotropic and strongly forward single scattering phase functions. It is found that the conical scattering by cylinders enhances the lateral scattering and therefore predicts a larger extent of lateral light scattering than models using rotationally invariant single scattering phase functions. However, it does not fully reach the levels of lateral scattering observed in measurements. It is argued that the hollow lumen of a wood fiber or dependent scattering effects must be considered for a complete description of lateral light scattering in paper. © 2013 Optical Society of America OCIS codes: (010.5620) Radiative transfer; (100.2810) Halftone image reproduction; (290.4210) Multiple scattering; (290.7050) Turbid media. References and links 1. J. A. C. Yule and W. J. Nielsen, “The penetration of light into paper and its effect on halftone reproduction,” Proceedings of TAGA 3, 65–67 (1951). 2. T. Linder and T. Löfqvist, “Anisotropic light propagation in paper,” Nord. Pulp Pap. Res. J. 27, 500–506 (2012). 3. M. Neuman, L. G. Coppel, and P. Edstrom, “Point spreading in turbid media with anisotropic single scattering,” Opt. Express 19, 1915–1920 (2011). 4. J. Arney, C. Arney, M. Katsube, and P. Engeldrum, “An MTF analysis of papers,” J. Imaging Sci. Technol. 40, 19–25 (1996). 5. J. S. Arney, J. Chauvin, J. Nauman, and P. G. Anderson, “Kubelka-Munk theory and the MTF of paper,” J. Imaging Sci. Technol. 47, 339–345 (2003). 6. L. G. Coppel, M. Neuman, and P. Edström, “Lateral light scattering in paper MTF simulation and measurement,” Opt. Express 19, 25181–25187 (2011). 7. C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (John Wiley and Sons, 1983). 8. A. Kienle, F. K. Forster, R. Diebolder, and R. Hibst, “Light propagation in dentin: influence of microstructure on anisotropy,” Phys. Med. Biol. 48, N7–N14 (2003). 9. A. Kienle, F. K. Forster, and R. Hibst, “Anisotropy of light propagation in biological tissue,” Opt. Lett. 29, 2617–2619 (2004). 10. H. He, N. Zeng, R. Liao, T. Yun, W. Li, Y. He, and H. Ma, “Application of sphere-cylinder scattering model to skeletal muscle,” Opt. Express 18, 15104–15112 (2010). 11. A. Kienle, C. Wetzel, A. Bassi, D. Comelli, P. Taroni, and A. Pifferi, “Determination of the optical properties of anisotropic biological media using an isotropic diffusion model,” J. Biomed. Opt. 12 (2007). #185971 $15.00 USD Received 25 Feb 2013; revised 15 Mar 2013; accepted 16 Mar 2013; published 22 Mar 2013 (C) 2013 OSA 25 March 2013 / Vol. 21, No. 6 / OPTICS EXPRESS 7835 12. A. Kienle, C. D’Andrea, F. Foschum, P. Taroni, and A. Pifferi, “Light propagation in dry and wet softwood,” Opt. Express 16, 9895–9906 (2008). 13. B. Peng, T. Ding, and P. Wang, “Propagation of polarized light through textile material,” Appl. Opt. 51, 6325– 6334 (2012). 14. T. Yun, N. Zeng, W. Li, D. Li, X. Jiang, and H. Ma, “Monte Carlo simulation of polarized photon scattering in anisotropic media,” Opt. Express 17, 16590–16602 (2009). 15. C. Fellers and B. Norman, Pappersteknik, 3rd ed. (Department of Pulp and Paper Chemistry and Technology, Royal Institute of Technology, 1996). 16. M. Ukishima, H. Kaneko, T. Nakaguchi, N. Tsumura, M. Hauta-Kasari, J. Parkkinen, and Y. Miyake, “A Simple Method to Measure MTF of Paper and Its Application for Dot Gain Analysis,” IEICE Trans. Fundam. Electron. Commun. Comput. Sci. E92A, 3328–3335.


Introduction
The underlying physics of light scattering in turbid media is important in a wide variety of fields.Light scattering in paper has a fundamental role as the optical appearance is one of its main functions.Lateral light scattering has a key role in the optical dot gain in paper, the so called Yule-Nielsen effect [1], where printed dots due to light scattering appear larger.It has also been suggested that monitoring scattering in paper can potentially be used in measuring fiber formation [2].
Proper models for lateral light scattering in turbid media are requried to correctly relate to the optical properties in materials such as paper.It has been shown by for example Neuman et al. [3] that both single scattering anisotropy and the medium mean free path has a big impact on the point spread function.Modulation transfer function (MTF) analysis is a convienient way of evaluating lateral light scattering [4].Arney et al. [5] compared measured values of the inverse frequency at full width half maximum of the MTF, k p , with a model based on the Kubelka-Munk (KM) scattering coefficient S.Recently Coppel et al. [6] compared these results with Monte Carlo simulations based on the Henyey-Greenstein phase function with two different asymmetry factors.They found that neither isotropic scattering nor strong forward scattering correctly could predict the large lateral scattering observed in the measurements.Since paper is highly anisotropic, based upon layers of cellulose fibers, phase functions that depend on the absolute direction of light within the structure are called for.Applying such phase functions to address similar observations has been suggested for other turbid media.Anisotropic directional dependent scattering can for example be modelled by using phase functions derived from an analytical solution of electromagnetic scattering by an infinite cylinder [7].This was first proposed by Kienle et al. [8,9] who observed the anisotropic shapes of backscattered point spread functions (PSF) for different stochastic representations of the cylinder alignments.Similar Monte Carlo simulations have since then been used to model light scattering in anisotropic structures such as biological tissue [10,11], softwood [12] and textile [13].A detailed description of the implementation of such a Monte Carlo model can be found in Yun et al. [14].One benefit of using analytical solutions is that the phase function is derived from material parameters like particle size and refractive index.This means that the asymmetry factor g is no longer required in the same way as in for example the Henyey-Greenstein phase function.Another benefit is that the scattering efficiency Q s can be obtained, it can be used to define how scattering distance depend on the absolute direction of photons propagating in aligned fiber structures.
The purpose of the present work is to utilize the cylinder phase function through Monte Carlo simulations to test if it can be used to mimic scattering in paper.The structural anisotropy is considered by aligning the cylinders isotropically in the plane of the paper and with a small Gaussian distribution in the thickness direction.The simulated MTFs are then compared to both the measured [5] and simulated values [6] of k p to test whether they can predict the amount of lateral scattering in the real paper samples.

Method
The Monte Carlo model in this work utilizes phase functions from an analytical solution of scattering by an infinite cylinder [7].For perpendicular incidence (ζ = 90 • ) light scatters in a disc around the longitudinal axis of the cylinder and for oblique incidence light scatters in a cone around the longitudinal axis where the phase function specifies the distributed intensity.The half angle of the scattering cone equals the incident angle ζ governing the directional dependencies observed by for example Kienle et al. [8].μ s is related to the density of cylinders C a , cylinder diameter d and scattering efficiency Q s for any incident angle ζ as The cylinder density is defined as the total length of cylinders per volume unit.It can be valid to test the performance of the model in relation to C a as it is the concentration of particles, especially since the cylinder phase function only is valid for plane electromagnetic waves.
The scattering parameters are determined by matching the reflectance and transmittance factors to those of the previously published measurements for each of a set of 22 samples.The shape of the simulated MTFs are then compared to the Monte Carlo simulations using the Henyey-Greenstain phase function with isotropic single scattering (g = 0.0) and forward single scattering (g = 0.8) and the simulated values of k p are also compared to the measured values.

Material parameters and parameter estimation
Paper is essentially an entangled network of fibers, however the model used in this work assumes that the fibers are not in contact with one another.Surface effects are not considered as the scattering only take place at the air-cylinder interfaces in the model.The cylinders are aligned isotropically in the plane of the simulated microstructure and are assumed to be following a Gaussian distribution with standard deviation σ z = 5 • in the thickness direction.A more detailed description of the cylinder alignment has been reported in [2].The refractive index of cellulose (n = 1.55) is used for the cylinders and the surrounding medium is air (n = 1.0).The cylinder diameter is chosen to be 20 μm as typical wood fiber dimension [15].The wavelength of the light was chosen to be λ = 510 nm as the experiments in Arney et al. used a green filter.We optimized C a iteratively using a Gauss-Newton method in the 0 • /d geometry with the Monte Carlo simulation tool to match the measured total transmission and reflection for each of the paper samples with their given thicknesses.The obtained cylinder density C a , scattering coefficient in the thickness direction μ s (z) and absorption coefficient μ a are shown in Table 1.

Simulation of the edge response
The Monte Carlo simulation tool (available at http://fibermc.sourceforge.net/) is used to compute the point spread function (PSF) for each of the 22 sets of paper sheets.Both for a single paper sheet and for an opaque pad of paper sheets.The resolution was chosen to 10 μm and the incident trajectories of the photons was tilted with a 20 • angle toward the surface normal to match the measurements.A 2D convolution between the simulated point spread function and a intensity distribution i(x, y) was used to derive the edge response ESF(x, y) = PSF(x, y) * i(x, y). ( The intensity distribution is uniformly distributed and cut off by a sharp knife edge, i.e. a matrix consisting of 1's on one side and 0's on the other.A similar convolution was used by Ukishima et al. [16] who considered the intensity distribution of pencil light together with the point spread function.This approach is different compared to Coppel et al. who simulated the edge response directly by distributing each of the photons over an area.Since the statistical response of a PSF

MTF
The line spread function (LSF) is obtained by taking the derivative of the edge spread function (ESF, see Fig. 1(a)).We then obtain the MTF by taking the Fourier transform of the LSF.The inverse frequency at full width half maximum of the MTF, k p , is defined as Only measured values of k p are compared to the simulations since Arney et al.only reported k p .The full MTF would have been preferred over the metric k p as it only indicates the amount of lateral light scattering and lose a lot of the information held by the MTF [6].

Results
The simulated MTFs of samples 2 and 17 plotted for the three phase functions are shown in Fig 1(b-c).The shape of the simulated MTFs is narrower at low frequencies for the simulations with the cylinder phase function compared to the previous Monte Carlo simulations but seems to agree well with simulations with anisotropic scattering (g = 0.8) at higher frequencies.The Figure 3 shows the values of k p in relation to C a for a single sheet and an opaque pad of paper samples.At low concentrations, or low scattering coefficients, the cylinder model seems to agree better with the previous simulations than with the measurements.It appears that the best correlation to the measurements are found in the middle region where C a lie between 3• 10 9 and 5• 10 9 m −2 .For larger values of C a a decresing correlation can be observed.model clearly predicted a higher value of k p than models using rotationally invariant isotropic and strong forward scattering.However, the cylinder model still underestimates k p compared to the measured values.This means that the model predicting less lateral scattering than observed.

Discussion and conclusions
The conical scattering by a cylinder has a tendency to sustain photon propagation in the plane of the paper thus increasing the lateral scattering.We originally thought that low correlation between simulations and measurements for the samples with low optical thickness (Fig. 3) were due to photons only scattering once before leaving the medium.Increased lateral scattering can only be observed if the photons manage to scatter so that they start propagating in the plane.However, the trend remains for the low concentrations in an opaque pad (Fig. 3(b)) and therefore contradicts this.The decreasing correlation at the larger values of the scattering coefficient, i.e. for very dense papers, is not suprising.The phase functions for scattering by cylinders is only valid for plane waves which means that the scatterers need to be in the farfield of one another.A possible explanation for the large values of the cylinder density is that the fibers in paper often are of banded or elliptical shapes.It shall be emphasized that it is a rough approximation to model a wood fiber as a homogenous cylinder.Additional simulations were made with adjusted values of the material parameters that affect the phase function, like e.g. the diameter of the scattering cylinders d and wavelength λ .We observed that this did not have any significant effect on the shape of the PSF and the resulting MTFs.This indicates that the model cannot predict more lateral light scattering for parameters within the natural parameter range of wood fibers.Arney suggested that the hollow lumen of a wood fiber could increase the lateral light scattering through a light-piping effect.This seem to be a reasonable explanation to why the model still predicts less lateral scattering than measured.Another reasonable explanation is that the concentration of particles is high, causing dependent scattering effects giving rise to interactions which are different from the interdependent scattering theory used in this work.Part of the lateral scattering can, however, be explained by the conical scattering by cylindrical objects.Compared to the Monte Carlo simulations using the Henyey-Greenstain phase function it roughly closes half the gap between earlier MC simulations and the experimental results.The cylinder phase function gives a good idea of how the lateral light scattering in fibrous materials is generated and has potential to better model for example optical dot gain.

#Fig. 1 .Fig. 2 .
Fig. 1.Simulated ESF and corresponding LSF in (a) and simulated MTFs for single paper sheet sample 2 in (b) and sample 17 in (c).The main difference between the samples is that sample 17 has a about four times larger scattering coefficient compared to sample 2. The red dots indicate the measured value of k p presented by Arney et al.

Fig. 3 .
Fig. 3. Measured and simulated values of k p for a single paper sheet (a) and an opaque pad of paper sheets (b) plotted against the density of cylinders C a .

Table 1 .
Sample thickness t, cylinder density C a , scattering coefficient μ s (z), absorption coefficient μ a and simulated values of k p for a single sheet and an opaque pad (k ∞ p ).