Characterizing liquid turbid media by frequency-domain photon-migration spectroscopy

Recent Advances in Light Penetration Depth for Postharvest Quality Evaluation of Fruits and Vegetables

Light penetration depth, as a characteristic parameter reflecting light attenuation and transmission in biological tissues, has been applied in nondestructive detection of fruits and vegetables. Recently, with emergence of new optical detection technologies, researchers have begun to explore methods evaluating optical properties of double-layer or even multilayer fruit and vegetable tissues due to the differences between peel and pulp in the chemical composition and physical properties, which has gradually promoted studies on light penetration depth. A series of demonstrated research on light penetration depth could ensure the accuracy of the optical information obtained from each layer of tissue, which is beneficial to enhance detection accuracy for quality assessment of fruits and vegetables. Therefore, the aim of this review is to give detailed outlines about the theory and principle of light penetration depth based on several emerging optical detection technologies and to focus primarily on its applications in the field of quality evaluation of fruits and vegetables, its future applicability in fruits and vegetables and the challenges it may face in the future.

Time-domain NIRS system based on supercontinuum light source and multi-wavelength detection: validation for tissue oxygenation studies

We present and validate a multi-wavelength time-domain near-infrared spectroscopy (TD-NIRS) system that avoids switching wavelengths and instead exploits the full capability of a supercontinuum light source by emitting and acquiring signals for the whole chosen range of wavelengths. The system was designed for muscle and brain oxygenation monitoring in a clinical environment. A pulsed supercontinuum laser emits broadband light and each of two detection modules acquires the distributions of times of flight of photons (DTOFs) for 16 spectral channels (used width 12.5 nm / channel), providing a total of 32 DTOFs at up to 3 Hz. Two emitting fibers and two detection fiber bundles allow simultaneous measurements at two positions on the tissue or at two source-detector separations. Three established protocols (BIP, MEDPHOT, and nEUROPt) were used to quantitatively assess the system's performance, including linearity, coupling, accuracy, and depth sensitivity. Measurements were performed on 32 homogeneous phantoms and two inhomogeneous phantoms (solid and liquid). Furthermore, measurements on two blood-lipid phantoms with a varied amount of blood and Intralipid provide the strongest validation for accurate tissue oximetry. The retrieved hemoglobin concentrations and oxygen saturation match well with the reference values that were obtained using a commercially available NIRS system (OxiplexTS) and a blood gas analyzer (ABL90 FLEX), except a discrepancy occurs for the lowest amount of Intralipid. In-vivo measurements on the forearm of three healthy volunteers during arterial (250 mmHg) and venous (60 mmHg) cuff occlusions provide an example of tissue monitoring during the expected hemodynamic changes that follow previously well-described physiologies. All results, including quantitative parameters, can be compared to other systems that report similar tests. Overall, the presented TD-NIRS system has an exemplary performance evaluated with state-of-the-art performance assessment methods.

Depletion-Induced Flocculation of Concentrated Emulsions Probed by Photon Density Wave Spectroscopy

Stable, creaming-free oil in water emulsions with high volume fractions of oil (ϕ = 0.05-0.40, density matched to water) and polysorbate 80 as an emulsifier were characterized without dilution by Photon Density Wave spectroscopy measuring light absorption and scattering behavior, the latter serving as the basis for droplet size distribution analysis. The emulsion with ϕ = 0.10 was used to investigate flocculation processes induced by xanthan as a semi-flexible linear nonabsorbing polymer. Different time regimes in the development of the reduced scattering coefficient μ_s' could be identified. First, a rapid, temperature-dependent change in μ_s' during the depletion process was observed. Second, the further decrease of μ_s' follows a power law in analogy to a spinodal demixing behavior, as described by the Cahn-Hilliard theory.

Correction of an adding-doubling inversion algorithm for the measurement of the optical parameters of turbid media

We present broadband measurements of the optical properties of tissue-mimicking solid phantoms using a single integrating sphere to measure the hemispherical reflectance and transmittance under a direct illumination at the normal incident angle. These measurements are traceable to reflectance and transmittance scales. An inversion routine using the output of the adding-doubling algorithm restricted to the reflectance and transmittance under a direct illumination was developed to produce the optical parameters of the sample along with an uncertainty budget at each wavelength. The results for two types of phantoms are compared to measurements by time-resolved approaches. The results between our method and these independent measurements agree within the estimated measurement uncertainties.

Algorithm for rapid determination of optical scattering parameters

Preliminary experiments at the NIST Spectral Tri-function Automated Reference Reflectometer (STARR) facility have been conducted with the goal of providing the diffuse optical properties of a solid reference standard with optical properties similar to human skin. Here, we describe an algorithm for determining the best-fit parameters and the statistical uncertainty associated with the measurement. The objective function is determined from the profile log likelihood, including both experimental and Monte Carlo uncertainties. Initially, the log likelihood is determined over a large parameter search box using a relatively small number of Monte Carlo samples such as 2·104. The search area is iteratively reduced to include the 99.9999% confidence region, while doubling the number of samples at each iteration until the experimental uncertainty dominates over the Monte Carlo uncertainty. Typically this occurs by 1.28·106 samples. The log likelihood is then fit to determine a 95% confidence ellipse. The inverse problem requires the values of the log likelihood on many points. Our implementation uses importance sampling to calculate these points on a grid in an efficient manner. Ultimately, the time-to-solution is approximately six times the cost of a Monte Carlo simulation of the radiation transport problem for a single set of parameters with the largest number of photons required. The results are found to be 64 times faster than our implementation of Particle Swarm Optimization.

Process analytical approaches for the coil-to-globule transition of poly(N-isopropylacrylamide) in a concentrated aqueous suspension

The coil-to-globule transition of poly(N-isopropylacrylamide) (PNIPAM) microgel particles suspended in water has been investigated in situ as a function of heating and cooling rate with four optical process analytical technologies (PAT), sensitive to structural changes of the polymer. Photon Density Wave (PDW) spectroscopy, Focused Beam Reflectance Measurements (FBRM), turbidity measurements, and Particle Vision Microscope (PVM) measurements are found to be powerful tools for the monitoring of the temperature-dependent transition of such thermo-responsive polymers. These in-line technologies allow for monitoring of either the reduced scattering coefficient and the absorption coefficient, the chord length distribution, the reflected intensities, or the relative backscatter index via in-process imaging, respectively. Varying heating and cooling rates result in rate-dependent lower critical solution temperatures (LCST), with different impact of cooling and heating. Particularly, the data obtained by PDW spectroscopy can be used to estimate the thermodynamic transition temperature of PNIPAM for infinitesimal heating or cooling rates. In addition, an inverse hysteresis and a reversible building of micrometer-sized agglomerates are observed for the PNIPAM transition process.

Limitations of turbidity process probes and formazine as their calibration standard

Turbidity measurements are frequently implemented for the monitoring of heterogeneous chemical, physical, or biotechnological processes. However, for quantitative measurements, turbidity probes need calibration, as is requested and regulated by the ISO 7027:1999. Accordingly, a formazine suspension has to be produced. Despite this regulatory demand, no scientific publication on the stability and reproducibility of this polymerization process is available. In addition, no characterization of the optical properties of this calibration material with other optical methods had been achieved so far. Thus, in this contribution, process conditions such as temperature and concentration have been systematically investigated by turbidity probe measurements and Photon Density Wave (PDW) spectroscopy, revealing an influence on the temporal formazine formation onset. In contrast, different reaction temperatures do not lead to different scattering properties for the final formazine suspensions, but give an access to the activation energy for this condensation reaction. Based on PDW spectroscopy data, the synthesis of formazine is reproducible. However, very strong influences of the ambient conditions on the measurements of the turbidity probe have been observed, limiting its applicability. The restrictions of the turbidity probe with respect to scatterer concentration are examined on the basis of formazine and polystyrene suspensions. Compared to PDW spectroscopy data, signal saturation is observed at already low reduced scattering coefficients.

Finite element modeling of light propagation in turbid media under illumination of a continuous-wave beam

Spatially resolved spectroscopy provides a means for measuring the optical properties of biological tissues, based on analytical solutions to diffusion approximation for semi-infinite media under the normal illumination of an infinitely small light beam. The method is, however, prone to error in measurement because the actual boundary condition and light beam often deviate from that used in deriving the analytical solutions. It is therefore important to quantify the effect of different boundary conditions and light beams on spatially resolved diffuse reflectance in order to improve the measurement accuracy of the technique. This research was aimed at using finite element method (FEM) to model light propagation in turbid media, subjected to normal illumination by a continuous-wave beam of infinitely small or finite size. Three types of boundary conditions [i.e., partial current (PCBC), extrapolated (EBC), and zero (ZBC)] were evaluated and compared against Monte Carlo (MC) simulations, since MC could provide accurate fluence rate and diffuse reflectance. The effect of beam size was also investigated. Overall results showed that FEM provided results as accurate as those of the analytical method when an appropriate boundary condition was applied. ZBC did not give satisfactory results in most cases. FEM-PCBC yielded a better fluence rate at the boundary than did FEM-EBC, while they were almost identical in predicting diffuse reflectance. Results further showed that FEM coupled with EBC effectively simulated spatially resolved diffuse reflectance under the illumination of a finite size beam. A large beam introduced more error, especially within the region of illumination. Research also confirmed an earlier finding that a light beam of less than 1 mm diameter should be used for estimation of optical parameters. FEM is effective for modeling light propagation in biological tissues and can be used for improving the optical property measurement by the spatially resolved technique.

Double-integrating-sphere system at the National Institute of Standards and Technology in support of measurement standards for the determination of optical properties of tissue-mimicking phantoms

There is a need for a common reference point that will allow for the comparison of the optical properties of tissue-mimicking phantoms. After a brief review of the methods that have been used to measure the phantoms for a contextual backdrop to our approach, this paper reports on the establishment of a standardized double-integrating-sphere platform to measure absorption and reduced scattering coefficients of tissue-mimicking biomedical phantoms. The platform implements a user-friendly graphical user interface in which variations of experimental configurations and model-based analysis are implemented to compute the coefficients based on a modified inverse adding-doubling algorithm allowing a complete uncertainty evaluation. Repeatability and validation of the measurement results of solid phantoms are demonstrated for three samples of different thicknesses, d = 5.08 mm, 7.09 mm, and 9.92 mm, with an absolute error estimate of 4.0% to 5.0% for the absorption coefficient and 11% to 12% for the reduced scattering coefficient (k = 2). The results are in accordance with those provided by the manufacturer. Measurements with different polarization angles of the incident light are also presented, and the resulting optical properties were determined to be equivalent within the estimated uncertainties.

Flexible tool for simulating the bulk optical properties of polydisperse spherical particles in an absorbing host: experimental validation

In this study, a flexible tool to simulate the bulk optical properties of polydisperse spherical particles in an absorbing host medium is described. The generalized Mie solution for Maxwell's equations is consulted to simulate the optical properties for a spherical particle in an absorbing host, while polydispersity of the particle systems is supported by discretization of the provided particle size distributions. The number of intervals is optimized automatically in an efficient iterative procedure. The developed tool is validated by simulating the bulk optical properties for two aqueous nanoparticle systems and an oil-in-water emulsion in the visible and near-infrared wavelength range, taking into account the representative particle sizes and refractive indices. The simulated bulk optical properties matched closely (R2 ≥ 0.899) with those obtained by reference measurements.

Optimization of the hyperspectral imaging-based spatially-resolved system for measuring the optical properties of biological materials

This paper reports on the optimization and assessment of a hyperspectral imaging-based spatially-resolved system for determination of the optical properties of biological materials over the wavelengths of 500-1,000 nm. Twelve model samples covering a wide range of absorption and reduced scattering coefficients were created to validate the hyperspectral imaging system, and their true values of absorption and reduced scattering coefficients were determined and then cross-validated using three commonly used methods (i.e., transmittance, integrating sphere, and empirical equation). Light beam and source-detector distance were optimized through Monte Carlo simulations and experiments for the model samples. The optimal light beam should be of Gaussian type with the diameter of less than 1 mm, and the optimal minimum and maximum source-detector distance should be 1.5 mm and 10-20 mean free paths, respectively. The optimized hyperspectral imaging-based spatially-resolved system achieved good estimation of the optical parameters.

Temperature-dependent optical properties of Intralipid measured with frequency-domain photon-migration spectroscopy

We present the temperature dependence of absorption and reduced scattering coefficients of 1.8% Intralipid measured by frequency-domain photon-migration spectroscopy between 710 and 850 nm. These measurements were made in the physiologically relevant 30 to 40 degrees C temperature range. The temperature coefficients for absorption were consistent during heating and cooling and follow closely other reported results. The change in absorption coefficient at 740 nm suggests that a minimum temperature change of 4 degrees C is observable within the error limits. We found that the reduced scattering coefficient shows a hysteresis with temperature at 740 nm. The temperature coefficient for reduced scattering determined from heating cycle measurements agrees with theory and other measurements within the error limits.