Next Article in Journal
Determining Single Photon Quantum States through Robust Waveguides on Chip
Previous Article in Journal
Design of Reflective Tunable Structural Color Metasurface Based on Guided-Mode Resonance Filter and Sb2S3
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Reproduction of Visible Absorbance Spectra of Highly Scattering Suspensions within an Integrating Sphere by Monte Carlo Simulation

Department of Physics, Faculty of Science, Tokyo University of Science, 1-3 Kagurazaka, Shinjuku-ku, Tokyo 162-8601, Japan
*
Author to whom correspondence should be addressed.
Photonics 2023, 10(7), 754; https://doi.org/10.3390/photonics10070754
Submission received: 27 April 2023 / Revised: 20 June 2023 / Accepted: 27 June 2023 / Published: 29 June 2023

Abstract

:
It is important to avoid the overestimation of absorption due to scattering when using absorption spectroscopy to measure scattering samples. We approached this issue by placing the sample inside an integrating sphere (IS) to collect the scattered light in all solid angles but encountered difficulty when determining the absorption coefficient from the absorbance because the light took various paths inside the IS and the sample. Therefore, by ray tracing inside the IS and the sample using Monte Carlo simulations (MC), we estimated the relationship between the absorption, scattering, anisotropy coefficients, and the measured absorbance. Scattering sample M, prepared by mixing polystyrene microspheres with trypan blue solution, and pure trypan blue solution for comparison were used as samples at various concentrations. MC reproduced the measurement results for the absorbance spectrum and its concentration dependence at 591 nm up to the measurement limit value. In addition, the saturated absorbance of sample M was lower than that of the trypan blue solution. This is because, from the distribution of distance d , light passed through the sample estimated by the MC, and more light with smaller d was detected due to scattering for higher concentration, resulting in a smaller increase in absorbance with the absorption coefficient.

1. Introduction

UV-visible absorption spectroscopy is a measurement method that investigates the electronic structure and composition of a material from the energy transitions due to optical absorption of the atoms or molecules that make it up. Because it can evaluate the physical properties of samples in a relatively simple, safe, and nondestructive manner, it is used in a wide range of fields, including materials science such as solar cell characterization [1], medical and biological fields such as protein concentration measurement [2], cell density measurement [3], and food testing [4]. However, since the size of the particles/clusters in these samples is often comparable with the wavelength of the measured light such that the scattering is not neglected [5,6,7], the transmitted light is attenuated not only by absorption but also by scattering. To accurately evaluate the amount of absorption of such scattering samples, it is necessary to detect not only transmitted light but also scattered light to estimate the contribution of scattering to the attenuation of transmitted light.
The instrument that makes this possible is the integrating sphere (IS), which is often used for absorption spectroscopy measurements of scattering samples. Generally, the sample is placed at the entrance of the IS as shown in Figure 1a. In this arrangement, forward scattered light as well as transmitted light is collected by the IS, enabling absorbance measurement without the influence of forward scattered light [8,9,10]. However, backscattered light does not enter the IS, so the measurement cannot be said to account for all scattered light [11].
Therefore, we placed a sample inside the IS (Figure 1b) to collect the scattered light at all solid angles with the IS, which enabled us to measure absorbance by taking into account the total scattering [12]. This arrangement has been used primarily in fluorescence quantum yield measurements [13,14,15,16,17]. Although it is also used in absorbance measurements, it is not yet commonly used [18,19]. One possible reason for this is that when the sample is placed inside the IS, it is difficult to estimate the absorption coefficient μ a from the measured absorbance A using the Lambert–Beer law ( A = μ a d log 10 e ) because the distance d light passing through the sample does not match the thickness h of the sample for the following reasons.
(1) Since the IS also collects reflected light from the sample surface, stray light that is detected without ever passing through the sample. If the sample is highly absorptive, the weak transmitted light will be buried by this reflected light [16,20,21,22,23];
(2) In the process, the IS may reach and pass through the sample again. In that case, the optical path length is longer than that of a single pass [16,24];
(3) Light leaks through apertures (entrance of irradiation and reference light to the sample) of the IS, which increases the apparent absorption [20];
(4) The IS is coated with highly reflective materials, such as BaSO4 and a special fluoropolymer, and their reflectance is greater than 90% and 95%, respectively, but not 100% in the visible region [25,26,27]. Therefore, light is slightly absorbed by the inner walls of the IS each time it reflects off these materials, which increases the apparent absorption.
Furthermore, although not directly related to placing the sample inside the IS, when measuring a scattering sample, (5) it must be considered that light does not travel straight through the sample in one direction due to scattering. In our previous work, we proposed an absorbance correction equation that reproduces the absorbance measured in the configuration shown in Figure 1b by taking into account (1) and (2) [12]. However, this correction formula does not fully take into account the effects of (3) and (5), and is not applicable to samples where multiple scattering is indicated. Therefore, we investigated a method to estimate the relationship between the absorption coefficient and measured absorbance from ray tracing using Monte Carlo simulations (MC).
In estimating the optical properties of a sample by combining IS and MC, one technique [28,29,30,31,32,33,34,35,36,37,38,39,40,41,42] that is often used is to measure the diffuse reflected light (backscattered light) and total transmitted light (forward scattered light and transmitted light) of a sample by IS and link the measured results with the absorption coefficient, scattering coefficient, and anisotropy coefficient by MC.
Diffuse reflected light and total transmitted light are mainly measured by the double integrating sphere (DIS) method and the single integrating sphere (SIS) method. In optical measurements with DIS [28,29,30,31,32,33,34], the sample is sandwiched between two ISs, and the diffuse reflectance is measured with the IS in front of the sample and the total transmittance is measured with the IS behind the sample. Since the reflection and transmission properties of a scattering sample can be obtained in a single measurement, optical measurements can be performed while modifying sample conditions, such as changes in temperature [28]. On the other hand, disadvantages of the method include the loss of side scattering light and an inability to measure reflection and transmission under the same conditions because the reflected light and transmitted light are measured with a different IS [38]. SIS is a technique for measuring total transmittance by placing a sample at the entrance of the IS; diffuse reflectance is measured by placing a sample at the diffuse reflectance measurement port (the opening opposite the entrance) of the IS [35,36,37,38,39,40,41,42]. Since only one IS is used, the reflected light and transmitted light can be measured under the same conditions. However, it is difficult to measure changes in the state of the sample, such as the temperature and time dependence of optical properties, because two measurements must be performed on a single sample.
In response to these, we attempted to estimate the relationship between the measured absorbance and the optical properties of the sample (absorption coefficient, scattering coefficient, and anisotropy coefficient) by MC [18,24] when the sample is placed inside an IS. With this method, it is possible to obtain the absorption properties of a non-fluorescent microparticle suspension in a single measurement with a single IS. In this paper, model samples with known optical properties were prepared at various concentrations, and the MC, with input information on the absorption coefficient, scattering coefficient, and scattering anisotropy of the samples, reproduced the concentration dependence of the absorbance measurement results from low absorbance to the measurement limit when the samples were placed inside the IS.

2. Materials and Methods

2.1. Sample Preparation

To evaluate the accuracy of the MC simulation, an absorption-scattering sample was prepared by mixing absorbers (non-scattering in the measurement wavelength region) and scatterers (non-absorbing in the measurement wavelength region) with known optical properties, where the measurement wavelength range was 300–800 nm. A trypan blue (Sigma-Aldrich, St. Louis, MO, USA) aqueous solution was used as the absorber and an aqueous suspension of polystyrene microspheres (Polysciences Inc., Warrington, PA, USA) with a diameter of 500 nm was used as the scatterer. A pseudo-scattering sample M was prepared by mixing 120 μM of trypan blue with polystyrene microspheres at a concentration of 0.0364 spheres/μm3. A set of sample Bs was prepared with nine levels of dilution of the mixed sample M as volume fractions of 2.5, 5.0, 12.5, 25.0, 37.5, 50.0, 62.5, 75.0%, and 100% (no dilution) with purified water.

2.2. Setup

The absorption spectra were measured by placing the samples in the arrangement shown in Figure 1b on the IS built into the spectrophotometer (SolidSpec-3700DUV, SHIMADZU CORPORATION, Tokyo, Japan). A top view of the setup is shown in Figure 2. The measurement was performed using the double-beam method, with the signal light hitting the sample through aperture S and the reference light entering the IS through aperture R. The light intensity was detected with a detector attached to the bottom of the IS. The detected intensity is the intensity of the signal light hitting the sample divided by the detected intensity of the reference light [9]. This reduced the intensity fluctuations of the light source [43]. The absorbance A was calculated from the detected intensity I 0 of the empty IS and the detected intensity I when the sample was placed inside the IS, using the following equation.
A = log 10 I I 0 .
A liquid sample of 600 μL was placed in a custom-made cylindrical cell made of quartz glass (Figure 2a), and the optical cell was placed into the IS as shown in Figure 1b and Figure 2b.

2.3. Optical Properties of Samples

To reproduce the absorbance measurement results by MC, information on the absorption coefficient of the absorber (trypan blue); scattering coefficient of the scatterer (polystyrene microspheres); and scattering anisotropy (the expected value of the cosine of the scattering angle) of sample M was needed. These optical properties were obtained by the following procedure, respectively. The absorption coefficient of the absorbent was obtained by the same procedure used in [12], where the absorbent, a trypan blue solution, was placed in a rectangular optical cell (made of quartz glass) with a 10 mm optical path length and measured in the arrangement shown in Figure 1a. The measurement results are shown in Figure 3. From the absorption spectrum (left), the concentration–absorption relationship at the absorption peak wavelength of 591 nm was extracted (right), and a linear fit was performed to obtain the relationship between the concentration and the absorption coefficient of trypan blue. The scattering and anisotropy coefficients of polystyrene microspheres were then calculated from the complex refractive index and diameter of the polystyrene microspheres based on the Mie scattering theory. For the calculation of the Mie scattering theory, we used the program “wiscombe-mie-code” written in Fortran, which was installed from [44]. The values of the real part of the complex refractive index of polystyrene microspheres were taken from [45]. The imaginary part was set to 0, as the absorption of visible light is negligible [46]. The input parameters used in the Mie scattering theory calculations and the calculated optical properties are shown in Table 1.

2.4. Monte Carlo Simulation (MC)

The behavior of the photon inside the IS and the sample was simulated in three dimensions by ray tracing with MC to estimate the measured absorbance A . MC was performed by an executable program (written in Fortran) that reproduced the actual measurement setup (Figure 2), including the size of the IS, the shape of the sample and its arrangement inside the IS, and the optical properties of each sample (absorption coefficient μ a , scattering coefficient μ s , anisotropy coefficient g obtained in the procedure in Section 2.3) as input parameters. The behavior of the photons inside the IS and the sample was simulated in three dimensions by MC. Following the flowchart shown in Figure 4, signal light (incident through the aperture S in Figure 2) and reference light (incident through the aperture R in Figure 2) simulations were performed for 50,000 photons. This number of photons was set to a sufficient number of simulations so that the standard deviation of the simulation results was sufficiently small (Appendix A). The photon is initially given a weight W = 1 , and each time there is absorption by the inner walls of the IS or by the sample, W is attenuated by an amount equivalent to the amount of absorption. For example, when traveling a distance p in an absorbing medium with absorption coefficient μ a , W is attenuated from W W e μ a p . The direction in which photons are reflected by the inner walls of the IS, whether they are reflected at the sample surface, and the direction of scattering by scatterers inside the sample were determined by random numbers based on the probability density distribution of each physical phenomenon [47]. The details of the simulation are as follows.

2.4.1. Photon Initialization

The initial position of the photon at the aperture of the IS was determined by a function based on the Gaussian distribution which assumed that the incident light had an intensity distribution that followed a Gaussian distribution. The spread of the incident light (dispersion of the Gaussian distribution) was set so that the saturated absorbance obtained from the measurement results matched the saturated absorbance obtained from the MC. This is because the spread of the incident light (dispersion of the Gaussian distribution) is reflected in the measurement limit of absorbance [12]. W was also initialized to W = 1 .

2.4.2. Inside the IS

A photon travels straight until it collides with the inner walls of the IS or the surface of the sample, and when the photon collides, an event is generated that corresponds to the point of impact. If a photon reaches the inner wall of the IS, it is reflected by the inner wall. The direction of reflection is determined randomly, assuming that the reflection has equal probability in all directions in the range of the elevation angle 0~π/2 and azimuth angle 0~2π. W is reduced by a small amount of absorption by the inner wall and is changed to W R i W . R i is the reflectance of the inner wall, which was determined by measuring the relative reflectance to the standard white plate with the diffuse reflectance measurement of the IS. When it reaches the surface of the sample, it is determined whether it is reflected or refracted on the sample surface by a random number according to the reflectance of that surface. The reflectance at the sample surface was derived from Fresnel’s equation. The reflectance was calculated from the angle of incidence on the sample under the following conditions: the incident light was unpolarized, the refractive index of the sample was 1.46 (quartz glass), which is the refractive index of the optical cell, and the refractive index outside the sample (inside the IS) was 1 (air). If it was transmitted, it was then processed inside the sample (Section 2.4.3). If it was reflected, the direction of travel was changed to the direction of specular reflection and sent straight inside the IS again.

2.4.3. Inside the Sample

The light that enters the sample travels straight for a distance s until it hits the scatterer and is scattered by the scatterer. s is determined by a random number based on the probability p d i s s = 1 e μ s s of traveling a distance s without ever hitting a scatterer. If it reaches the sample surface before reaching the scatterer, it travels to the boundary (distance p ) and determines whether it is reflected or refracted on the boundary by a random number using the same procedure as in Section 2.4.1. If the photon is refracted, it exits the sample and returns to the IS’s internal processing (Section 2.4.1). In the case of reflection, the direction of travel is changed to the direction of specular reflection, and the travel distance s to the scatterer is changed to s s p where s is subtracted from p , and the sample is moved straight through the sample again. If it reaches the scatterer, it is scattered and changes its direction of travel. The scattering direction is determined by a random number based on the scattering phase correlation number p s c a cos θ . p s c a cos θ is given by the Henyey–Greenstein function, an approximation of the Mie scattering theory [47], expressed by
p s c a cos θ = 1 g 2 2 1 + g 2 2 g cos θ .
Each time the photon travels inside the sample, it is absorbed from the absorbing medium, so the weight of the photon is changed to W W e μ a d for the travel distance d .

2.4.4. End of Photon Processing and Estimation of Absorbance

The three end conditions of the simulation for a single photon are: the photon reaches the detector and is detected, reaches the aperture, exits the IS and W is negligibly small ( W < W e n d = 1.0 × 10 6 ). If a photon is detected, the current W is added to the sum of the W of the detected photons, W t o t a l . The average value of W t o t a l , W t o t a l / N , corresponds to the ratio of the detected intensity to the incident light intensity ( W = 1 ), i.e., the measured transmission T ¯ . Because of the double-beam system, T ¯ is expressed as T ¯ = T ¯ s / T ¯ R , where T ¯ s is the incident light on the sample divided by T ¯ R , the reference light. Thus, the measured absorbance A is calculated by
A = log 10 T ¯ = log 10 T ¯ s T ¯ R .

3. Results

3.1. Absorbance Measured When a Sample Is Placed inside the IS

Figure 5 shows the absorbance measurement results for the setup in Figure 2. A broad absorption peak at 591 nm was observed for both the pure trypan blue solution (non-scattering sample), and the sample M (scattering sample) which was mixed with polystyrene microspheres as scatterers. In general, according to the Lambert–Beer law, absorbance and concentration should have a linear relationship, but when the sample is placed inside the IS, the relationship becomes nonlinear. In addition, sample M, which contained scatterers, and showed a smaller measurement limit than the pure trypan blue solution. Even in the case of pure trypan blue solution, the saturation absorbance is about Abs. = 1.5 which is less than the instrument limit (Abs. = 5).

3.2. Reproduction of Absorbance Measurement Results by MC

Next, the absorbance estimated by the MC was compared to the measured absorbance. Figure 6a and Figure 6b compare the absorbance estimated by MC with the measurements shown in Figure 5a and Figure 5b, respectively. Absorption spectrum estimates by MC were performed every 10 nm in the wavelength range 400~700 nm for concentrations of 9.0, 30.0, and 120.0 μM of trypan blue. Figure 6 shows that MC can reproduce the absorbance measurement results when the sample is placed inside the IS for both non-scattering and scattering samples.

4. Discussion

Sample M (the scattering sample), which was a trypan blue solution (non-scattering sample) mixed with non-absorbing scatterers, had the same absorption coefficient if the absorber, trypan blue, was at the same concentration. Therefore, it was expected that there would be no significant difference in the concentration dependence of the measured absorbance of the two when the absorbance measurement was performed considering total scattered light by placing a sample inside the IS. However, as is shown in Figure 6, there were differences in the shapes of the absorption peaks and saturation absorbance depending on the presence or absence of scattering. We considered this to be due to the fact that the distance that light travels inside the sample changed with the presence of scattering.
The absorbance is generally proportional to μ a d according to the Lambert–Beer law, which is the product of the absorption coefficient μ a and the distance d that light passes through the sample until it is detected. Therefore, the smaller the value of d , the smaller the amount of the increase in absorbance A with respect to the increase in absorption coefficient μ a . Based on the above, the reason why the increase in absorbance relative to the increase in concentration (increase in μ a ) becomes smaller as the concentration (absorption) of the sample increases is considered as follows: light with relatively small d is detected more often as the concentration of the sample increases. In order to verify this, we estimated the distribution of d of the detected light for signal light (incident light from the aperture S shown in Figure 2) by MC. Figure 7 is a histogram of d by the MC showing what percentage of the photons detected were classified into each optical path length range. The region L m i d is the range of d that can be taken when signal light travels straight through the sample, and it is classified into region L l o n g or L s h o r t if it is longer or shorter than L m i d , respectively. If the photon does not transmit through the sample due to reflection on its surface, the photon is classified as d = 0 .
For non-scattering samples, there are no photons classified in the region L s h o r t . This is because light does not exit on its way through the sample by scattering as the non-scattering sample has no scatterers. Therefore, light travels straight through the sample if it is not reflected at the sample’s surface. In the case of a high-absorption sample, since the transmitted light becomes weak, the light of L l o n g , L m i d is buried in the light of d = 0 . For light with d = 0 , A does not increase even if μ a increases and the increase in absorbance plateaus. This is also explained in the absorbance correction formula when the sample is placed inside the IS, as shown in the previous paper [12].
On the other hand, in the case of the scattering sample, there are photons contained in L s h o r t , and the ratio increases as the concentration (absorption) of the sample increases. It is suggested that the presence of photons classified as L s h o r t makes the d of the scattering samples smaller than that of the non-scattering samples in the high-concentration (high-absorption) region, resulting in lower saturation absorbance. The reason why the ratio of L s h o r t is higher at higher concentrations is that the amount of light exits the sample while traveling straight through the sample due to scattering increases as the concentration of the scatterer increases. At the same time, there is also an increase in the optical path length due to scattering, but when d is sufficiently long in the high-absorption region, the intensity is so weak that almost no light is absorbed.

5. Conclusions

Absorbance measurements were performed considering the effect of scattering by placing the sample inside the IS and using it to collect scattered light at all solid angles. Using this method, the distance d that light passes through the sample before it is detected does not match the thickness of the sample due to the repeated scattering of light on the inner walls of the IS, which takes various optical paths inside the IS, and in the case of a scattering sample, takes various paths inside the sample due to scattering. Therefore, it is difficult to quantitatively evaluate the absorption coefficient μ a from the Lambert–Beer law ( A = μ a d log 10 e ). Therefore, we estimated the relationship between the measured absorbance and absorption coefficient by simulating the optical path tracking inside the IS and inside the sample, which reflected the actual optical setup by MC.
For the scattering sample M with known optical properties and the non-scattering sample for comparison, absorbance measurements were performed with the sample inside the IS, respectively, and the measured absorbance estimated by the MC was compared with the measured results. A trypan blue solution was used for the non-scattering sample, and sample M was made by mixing it with polystyrene microspheres as scatterers. The MC results reproduced the experimental results well, reproducing the shape of the absorption spectrum and the behavior of absorbance with respect to sample concentrations at the peak absorption wavelength (591 nm) for both non-scattering and scattering samples.
In the absorbance measurements inside the IS, the scattering samples showed a smaller saturation absorbance than the non-scattering samples. We consider that this is because, in the case of scattering samples, d decreases with increasing sample concentrations, resulting in a smaller increase in absorbance with increasing sample concentrations, because the absorbance depends on the product of the concentration of the sample and d according to to the Lambert–Beer law. This was verified by estimating the d distribution of the detected light by MC. The simulation results showed that for the scattering sample, the distribution of d of the detected light is biased toward smaller values as the concentration (absorption) of the sample increases, compared to the non-scattering sample. This is thought to be because the higher the concentration, the more scatterers there are, increasing the probability that light will travel out of the sample by scattering on its way straight through the sample.

Author Contributions

Conceptualization, E.T. and A.M.; methodology, A.M., K.Y. and E.T.; software, A.M.; validation, A.M.; formal analysis, A.M.; investigation, A.M.; resources, E.T. and K.Y.; data curation, A.M.; writing—original draft preparation, A.M.; writing—review and editing, E.T. and K.Y.; visualization, A.M.; supervision, E.T.; project administration, E.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors thank Takayoshi Kobayashi and Keisuke Seto for helpful discussions. This work was supported by JST SPRING, Grant Number JPMJSP2151.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The number of photons used in the Monte Carlo simulations was set to a sufficient number of simulations so that the standard deviation of the simulation results was sufficiently small. Figure A1a,b shows the variation in the simulation results when the simulation shown on the right side of Figure 6a,b of this manuscript was performed 10 times. The vertical axis is the standard deviation, and the horizontal axis is the number of photons for one simulation. The change in standard deviation from 50,000 photons to 75,000 photons remained almost flat, and we concluded that increasing the number of photons would only increase the simulation time. Table A1 and Table A2 show the absorbance average (Ā) and standard deviation (σ) for each photon number in the simulation.
Figure A1. Standard deviation in the MC-simulated absorbance of non-scattering samples (a); and scattering samples (b) for each photon number. In the main text, simulation results with 50,000 photons were adopted (Figure 6-right).
Figure A1. Standard deviation in the MC-simulated absorbance of non-scattering samples (a); and scattering samples (b) for each photon number. In the main text, simulation results with 50,000 photons were adopted (Figure 6-right).
Photonics 10 00754 g0a1
Table A1. The absorbance average (Ā) and standard deviation (σ) for each photon number in the simulation of non-scattering samples (50,000 photons in Figure 6a-right).
Table A1. The absorbance average (Ā) and standard deviation (σ) for each photon number in the simulation of non-scattering samples (50,000 photons in Figure 6a-right).
Concentration of
Absorbar
The Number of Photons
1000500010,00025,00050,00075,000
ĀσĀσĀσĀσĀσĀσ
3 μM0.1600.0330.1580.0130.1550.0090.1520.0070.1530.0050.1520.005
9 μM0.3070.0240.3300.0070.3290.0080.3300.0090.3290.0060.3270.004
15 μM0.4760.0440.4860.0320.4940.0110.4930.0080.4900.0030.4920.005
30 μM0.8540.0300.8350.0250.8430.0180.8430.0140.8430.0060.8420.007
45 μM1.1250.0791.0860.0241.1160.0191.1080.0131.1140.0121.1120.009
60 μM1.3580.0821.3080.0551.2840.0431.2850.0251.2870.0231.2800.014
75 μM1.4650.1351.3660.0591.3770.0381.3810.0201.3720.0221.3760.014
90 μM1.4620.1261.4270.0851.4170.0451.4280.0411.4270.0261.4100.016
120 μM1.4390.1181.4610.0541.4210.0471.4390.0431.4400.0231.4380.025
Table A2. The absorbance average (Ā) and standard deviation (σ) for each photon number in the simulation of scattering samples (50,000 photons in Figure 6b-right).
Table A2. The absorbance average (Ā) and standard deviation (σ) for each photon number in the simulation of scattering samples (50,000 photons in Figure 6b-right).
Concentration of
Absorbar
The Number of Photons
1000500010,00025,00050,00075,000
ĀσĀσĀσĀσĀσĀσ
3 μM0.1460.0240.1520.0150.1430.0140.1460.0060.1470.0060.1410.004
9 μM0.3460.0420.3470.0070.3460.0080.3470.0060.3520.0030.3520.005
15 μM0.5330.0530.5260.0140.5300.0110.5250.0090.5230.0040.5250.005
30 μM0.7990.0630.7860.0340.7720.0170.7880.0130.7770.0090.7840.008
45 μM0.8740.0860.9010.0340.8840.0110.8680.0130.8800.0070.8750.008
60 μM0.9140.0730.9170.0340.9150.0200.9180.0100.9180.0090.9240.006
75 μM0.9620.1090.9410.0350.9340.0290.9430.0120.9400.0060.9470.009
90 μM0.9530.0750.9510.0510.9440.0170.9500.0160.9600.0090.9530.007
120 μM0.9010.0710.9710.0260.9520.0260.9650.0110.9630.0120.9610.009

References

  1. Adimule, V.; Suryavanshi, A.; Bc, Y.; Nandi, S.S. A Facile Synthesis of Poly(3-Octyl Thiophene):Ni0.4Sr0.6TiO3 Hybrid Nanocomposites for Solar Cell Applications. Macromol. Symp. 2020, 392, 2000001. [Google Scholar] [CrossRef]
  2. McKechnie, W.S.; Tugcu, N.; Kandula, S. Accurate and Rapid Protein Concentration Measurement of In-Process, High Concentration Protein Pools. Biotechnol. Prog. 2018, 34, 1234–1241. [Google Scholar] [CrossRef] [PubMed]
  3. Drieschner, T.; Ostertag, E.; Boldrini, B.; Lorenz, A.; Brecht, M.; Rebner, K. Direct Optical Detection of Cell Density and Viability of Mammalian Cells by Means of UV/VIS Spectroscopy. Anal. Bioanal. Chem. 2020, 412, 3359–3371. [Google Scholar] [CrossRef] [PubMed]
  4. Aleixandre-Tudo, J.L.; Buica, A.; Nieuwoudt, H.; Aleixandre, J.L.; du Toit, W. Spectrophotometric Analysis of Phenolic Compounds in Grapes and Wines. J. Agric. Food Chem. 2017, 65, 4009–4026. [Google Scholar] [CrossRef]
  5. Numata, Y.; Nair, S.V.; Nakagawa, K.; Ishino, H.; Kobayashi, T.; Tokunaga, E. Optical Size Effect of Organic Nanocrystals Studied by Absorption Spectroscopy within an Integrating Sphere. Chem. Phys. Lett. 2014, 601, 128–133. [Google Scholar] [CrossRef]
  6. Wrigglesworth, E.G.; Johnston, J.H. Mie Theory and the Dichroic Effect for Spherical Gold Nanoparticles: An Experimental Approach. Nanoscale Adv. 2021, 3, 3530–3536. [Google Scholar] [CrossRef]
  7. Ishino, H.; Nair, S.V.; Nakagawa, K.; Kobayashi, T.; Tokunaga, E. Effect of Light Scattering on the Transmission Spectra of Organic Nanocrystals. Appl. Phys. Lett. 2011, 99, 053304. [Google Scholar] [CrossRef] [Green Version]
  8. Dazzi, A.; Deniset-Besseau, A.; Lasch, P. Minimising Contributions from Scattering in Infrared Spectra by Means of an Integrating Sphere. Analyst 2013, 138, 4191–4201. [Google Scholar] [CrossRef] [Green Version]
  9. Gaigalas, A.K.; He, H.-J.; Wang, L. Measurement of Absorption and Scattering With an Integrating Sphere Detector: Application to Microalgae. J. Res. Natl. Inst. Stand. Technol. 2009, 114, 69–81. [Google Scholar] [CrossRef]
  10. Merzlyak, M.N.; Razi Naqvi, K. On Recording the True Absorption Spectrum and the Scattering Spectrum of a Turbid Sample: Application to Cell Suspensions of the Cyanobacterium Anabaena Variabilis. J. Photochem. Photobiol. B Biol. 2000, 58, 123–129. [Google Scholar] [CrossRef]
  11. Röttgers, R.; Gehnke, S. Measurement of Light Absorption by Aquatic Particles: Improvement of the Quantitative Filter Technique by Use of an Integrating Sphere Approach. Appl. Opt. 2012, 51, 1336–1351. [Google Scholar] [CrossRef] [Green Version]
  12. Mori, A.; Yamashita, K.; Tabata, Y.; Seto, K.; Tokunaga, E. Absorbance Spectroscopy of Light Scattering Samples Placed inside an Integrating Sphere for Wide Dynamic Range Absorbance Measurement. Rev. Sci. Instrum. 2021, 92, 123103. [Google Scholar] [CrossRef]
  13. Würth, C.; Pauli, J.; Lochmann, C.; Spieles, M.; Resch-Genger, U. Integrating Sphere Setup for the Traceable Measurement of Absolute Photoluminescence Quantum Yields in the Near Infrared. Anal. Chem. 2012, 84, 1345–1352. [Google Scholar] [CrossRef]
  14. Leyre, S.; Coutino-Gonzalez, E.; Joos, J.J.; Ryckaert, J.; Meuret, Y.; Poelman, D.; Smet, P.F.; Durinck, G.; Hofkens, J.; Deconinck, G.; et al. Absolute Determination of Photoluminescence Quantum Efficiency Using an Integrating Sphere Setup. Rev. Sci. Instrum. 2014, 85, 123115. [Google Scholar] [CrossRef]
  15. Quintanilla, M.; Hemmer, E.; Marques-Hueso, J.; Rohani, S.; Lucchini, G.; Wang, M.; Zamani, R.R.; Roddatis, V.; Speghini, A.; Richards, B.S.; et al. Cubic versus Hexagonal—Phase, Size and Morphology Effects on the Photoluminescence Quantum Yield of NaGdF 4 :Er 3+ /Yb 3+ Upconverting Nanoparticles. Nanoscale 2022, 14, 1492–1504. [Google Scholar] [CrossRef]
  16. MacDougall, S.K.W.; Ivaturi, A.; Marques-Hueso, J.; Richards, B.S. Measurement Procedure for Absolute Broadband Infrared Up-Conversion Photoluminescent Quantum Yields: Correcting for Absorption/Re-Emission. Rev. Sci. Instrum. 2014, 85, 063109. [Google Scholar] [CrossRef]
  17. Jones, C.M.S.; Panov, N.; Skripka, A.; Gibbons, J.; Hesse, F.; Bos, J.-W.G.; Wang, X.; Vetrone, F.; Chen, G.; Hemmer, E.; et al. Effect of Light Scattering on Upconversion Photoluminescence Quantum Yield in Microscale-to-Nanoscale Materials. Opt. Express 2020, 28, 22803–22818. [Google Scholar] [CrossRef]
  18. Villanueva, Y.; Veenstra, C.; Steenbergen, W. Measuring Absorption Coefficient of Scattering Liquids Using a Tube inside an Integrating Sphere. Appl. Opt. 2016, 55, 3030–3038. [Google Scholar] [CrossRef]
  19. Darby, B.L.; Auguié, B.; Meyer, M.; Pantoja, A.E.; Le Ru, E.C. Modified Optical Absorption of Molecules on Metallic Nanoparticles at Sub-Monolayer Coverage. Nat. Photonics 2016, 10, 40–45. [Google Scholar] [CrossRef] [Green Version]
  20. Tassan, S.; Ferrari, G.M. Variability of Light Absorption by Aquatic Particles in the Near-Infrared Spectral Region. Appl. Opt. 2003, 42, 4802–4810. [Google Scholar] [CrossRef]
  21. Chess, D.J.; Billings, E.; Covian, R.; Glancy, B.; French, S.; Taylor, J.; de Bari, H.; Murphy, E.; Balaban, R.S. Optical Spectroscopy in Turbid Media Using an Integrating Sphere: Mitochondrial Chromophore Analysis during Metabolic Transitions. Anal. Biochem. 2013, 439, 161–172. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Haardt, H.; Maske, H. Specific in Vivo Absorption Coefficient of Chlorophyll a at 675 Nm1. Limnol. Oceanogr. 1987, 32, 608–619. [Google Scholar] [CrossRef] [Green Version]
  23. Aduev, B.P.; Nurmukhametov, D.R.; Belokurov, G.M.; Zvekov, A.A.; Kalenskii, A.V.; Nikitin, A.P.; Liskov, I.Y. Integrating Sphere Study of the Optical Properties of Aluminum Nanoparticles in Tetranitropentaerytrite. Tech. Phys. 2014, 59, 1387–1392. [Google Scholar] [CrossRef]
  24. Nelson, N.B.; Prézelin, B.B. Calibration of an Integrating Sphere for Determining the Absorption Coefficient of Scattering Suspensions. Appl. Opt. 1993, 32, 6710–6717. [Google Scholar] [CrossRef]
  25. How to Choose Integrating Spheres and White Reference Plates. Available online: https://www.shimadzu.com/an/service-support/technical-support/technical-information/uv-vis/uv-ap/white/index.html (accessed on 16 June 2023).
  26. Gaballah, A.; Abdelmageed, A. Investigation of Energy Efficiency Index for Indoor LED Lighting Units. Semicond. Phys. Quantum Electron. Optoelectron. 2023, 26, 97–104. [Google Scholar] [CrossRef]
  27. Chiu, Y.-C.; Liu, W.-R.; Yeh, Y.-T.; Jang, S.-M.; Chen, T.-M. The Luminescence Properties and Ce3+ → Tb3+ Energy Transfer of Ca3Y2(Si3O9)2. ECS Trans. 2009, 25, 157. [Google Scholar] [CrossRef]
  28. Kara, E.; Çilesiz, İ.; Gülsoy, M. Monitoring System for Investigating the Effect of Temperature Change on Optical Properties. Lasers Med. Sci. 2018, 33, 1763–1768. [Google Scholar] [CrossRef]
  29. Yust, B.G.; Mimun, L.C.; Sardar, D.K. Optical Absorption and Scattering of Bovine Cornea, Lens, and Retina in the near-Infrared Region. Lasers Med. Sci. 2012, 27, 413–422. [Google Scholar] [CrossRef] [Green Version]
  30. Fukutomi, D.; Ishii, K.; Awazu, K. Highly Accurate Scattering Spectra of Strongly Absorbing Samples Obtained Using an Integrating Sphere System by Considering the Angular Distribution of Diffusely Reflected Light. Lasers Med. Sci. 2015, 30, 1335–1340. [Google Scholar] [CrossRef]
  31. Fukutomi, D.; Ishii, K.; Awazu, K. Determination of the Scattering Coefficient of Biological Tissue Considering the Wavelength and Absorption Dependence of the Anisotropy Factor. Opt. Rev. 2016, 23, 291–298. [Google Scholar] [CrossRef]
  32. Hülsbusch, M.; Blazek, V. Photon-Tissue Interaction Modelled by Monte Carlo Method for Optimizing Optoelectronic Sensor Concepts. In Studies in Skin Perfusion Dynamics: Photoplethysmography and its Applications in Medical Diagnostics; Blazek, V., Kumar, V.J., Leonhardt, S., Mukunda Rao, M., Eds.; Biological and Medical Physics Biomedical Engineering; Springer: Singapore, 2021; pp. 163–175. ISBN 9789811554490. [Google Scholar]
  33. Saccomandi, P.; Larocca, E.S.; Rendina, V.; Schena, E.; D’Ambrosio, R.; Crescenzi, A.; Di Matteo, F.M.; Silvestri, S. Estimation of Optical Properties of Neuroendocrine Pancreas Tumor with Double-Integrating-Sphere System and Inverse Monte Carlo Model. Lasers Med. Sci. 2016, 31, 1041–1050. [Google Scholar] [CrossRef]
  34. Liu, X.; Wu, Y. Monte-Carlo Optical Model Coupled with Inverse Adding-Doubling for Building Integrated Photovoltaic Smart Window Design and Characterisation. Sol. Energy Mater. Sol. Cells 2021, 223, 110972. [Google Scholar] [CrossRef]
  35. Morales-Cruzado, B.; Pérez-Gutiérrez, F.G.; de Lange, D.F.; Romero-Méndez, R. Study of the Effect Introduced by an Integrating Sphere on the Temporal Profile Characterization of Short Laser Pulses Propagating through a Turbid Medium. Appl. Opt. 2015, 54, 2383–2390. [Google Scholar] [CrossRef]
  36. Friebel, M.; Helfmann, J.; Netz, U.J.; Meinke, M.C. Influence of Oxygen Saturation on the Optical Scattering Properties of Human Red Blood Cells in the Spectral Range 250 to 2000 Nm. J. Biomed. Opt. 2009, 14, 034001. [Google Scholar] [CrossRef]
  37. Salomatina, E.V.; Jiang, B.; Novak, J.; Yaroslavsky, A.N. Optical Properties of Normal and Cancerous Human Skin in the Visible and Near-Infrared Spectral Range. J. Biomed. Opt. 2006, 11, 064026. [Google Scholar] [CrossRef] [Green Version]
  38. Simpson, C.R.; Kohl, M.; Essenpreis, M.; Cope, M. Near-Infrared Optical Properties Ofex Vivohuman Skin and Subcutaneous Tissues Measured Using the Monte Carlo Inversion Technique. Phys. Med. Biol. 1998, 43, 2465–2478. [Google Scholar] [CrossRef]
  39. Bergmann, F.; Foschum, F.; Zuber, R.; Kienle, A. Precise Determination of the Optical Properties of Turbid Media Using an Optimized Integrating Sphere and Advanced Monte Carlo Simulations. Part 2: Experiments. Appl. Opt. 2020, 59, 3216–3226. [Google Scholar] [CrossRef]
  40. Hamdy, O.; Mohammed, H.S. Investigating the Transmission Profiles of 808 Nm Laser through Different Regions of the Rat’s Head. Lasers Med. Sci. 2021, 36, 803–810. [Google Scholar] [CrossRef]
  41. Nilsson, A.M.K.; Berg, R.; Andersson-Engels, S. Measurements of the Optical Properties of Tissue in Conjunction with Photodynamic Therapy. Appl. Opt. 1995, 34, 4609–4619. [Google Scholar] [CrossRef] [Green Version]
  42. Vincely, V.D.; Vishwanath, K. Accuracy of Retrieving Optical Properties from Liquid Tissue Phantoms Using a Single Integrating Sphere. Appl. Opt. 2022, 61, 375–385. [Google Scholar] [CrossRef]
  43. C101-E123 UV Talk Letter Vol. 9. Available online: https://www.shimadzu.com/an/sites/shimadzu.com.an/files/pim/pim_document_file/journal/talk_letters/9633/jpa113001.pdf (accessed on 16 June 2023).
  44. Mie Scattering. Available online: https://omlc.org/software/mie/ (accessed on 16 April 2023).
  45. Sultanova, N.; Kasarova, S.; Nikolov, I. Dispersion Properties of Optical Polymers. Acta Phys. Pol. A 2009, 116, 585–587. [Google Scholar] [CrossRef]
  46. Shahin, A.; Bachir, W.; Sayem El-Daher, M. Polystyrene Microsphere Optical Properties by Kubelka–Munk and Diffusion Approximation with a Single Integrating Sphere System: A Comparative Study. J. Spectrosc. 2019, 2019, e3406319. [Google Scholar] [CrossRef]
  47. Welch, A.J.; Gemert, M.J.C. Optical-Thermal Response of Laser-Irradiated Tissue; Plenum Press: New York, NY, USA, 1995. [Google Scholar]
Figure 1. Setup for absorbance measurement using the IS: (a) the sample is placed in the aperture. This arrangement is often used in absorbance measurements using the IS; and (b) the sample is placed inside the IS. This paper discusses absorbance measurements in this arrangement.
Figure 1. Setup for absorbance measurement using the IS: (a) the sample is placed in the aperture. This arrangement is often used in absorbance measurements using the IS; and (b) the sample is placed inside the IS. This paper discusses absorbance measurements in this arrangement.
Photonics 10 00754 g001
Figure 2. (a) Setup of the optical system used in the measurement (top view of Figure 1a). Light is input through aperture S to measure the absorption of the sample, and reference light is input through aperture R to correct for fluctuations in the light source; and (b) custom-made optical cell designed to be installed inside the IS. A liquid sample was placed in this cell and inserted inside the IS through the port for the diffuse reflection measurement, and the back of the cell was pressed down with a white plate made of the same material as the inner wall of the IS.
Figure 2. (a) Setup of the optical system used in the measurement (top view of Figure 1a). Light is input through aperture S to measure the absorption of the sample, and reference light is input through aperture R to correct for fluctuations in the light source; and (b) custom-made optical cell designed to be installed inside the IS. A liquid sample was placed in this cell and inserted inside the IS through the port for the diffuse reflection measurement, and the back of the cell was pressed down with a white plate made of the same material as the inner wall of the IS.
Photonics 10 00754 g002
Figure 3. Concentration dependence of the absorption spectrum of trypan blue (left) and the absorption coefficient at the 591 nm absorption peak wavelength (right) are shown. As demonstrated in the figure on the right, the absorption coefficient is proportional to the concentration. Therefore, an approximate straight line (dashed line) can be drawn to the measured values (circles). The absorption coefficient can be estimated from the concentration using this approximate line.
Figure 3. Concentration dependence of the absorption spectrum of trypan blue (left) and the absorption coefficient at the 591 nm absorption peak wavelength (right) are shown. As demonstrated in the figure on the right, the absorption coefficient is proportional to the concentration. Therefore, an approximate straight line (dashed line) can be drawn to the measured values (circles). The absorption coefficient can be estimated from the concentration using this approximate line.
Photonics 10 00754 g003
Figure 4. Flowchart of the MC to estimate the absorbance measured in the setup shown in Figure 2, beginning with the initialization of one photon at ‘Photon initialization’ and continuing until the end of the flowchart at ‘end’. For each sample, this flowchart was executed for 50,000 photons, with one incident at the sample through aperture S and the other through the aperture R.
Figure 4. Flowchart of the MC to estimate the absorbance measured in the setup shown in Figure 2, beginning with the initialization of one photon at ‘Photon initialization’ and continuing until the end of the flowchart at ‘end’. For each sample, this flowchart was executed for 50,000 photons, with one incident at the sample through aperture S and the other through the aperture R.
Photonics 10 00754 g004
Figure 5. Results of absorbance measurements of non-scattering samples, Trypan blue solution (a), and scattering samples, sample M (generated by mixing polystyrene spheres in an aqueous solution of trypan blue) (b), performed in the setup shown in Figure 2. The left panel shows the absorption spectrum, and the right panel shows the relationship between concentration and absorbance at the absorption peak wavelength of 591 nm. The concentrations shown in the graphs refer to the concentration of the absorbent trypan blue in both (a,b).
Figure 5. Results of absorbance measurements of non-scattering samples, Trypan blue solution (a), and scattering samples, sample M (generated by mixing polystyrene spheres in an aqueous solution of trypan blue) (b), performed in the setup shown in Figure 2. The left panel shows the absorption spectrum, and the right panel shows the relationship between concentration and absorbance at the absorption peak wavelength of 591 nm. The concentrations shown in the graphs refer to the concentration of the absorbent trypan blue in both (a,b).
Photonics 10 00754 g005
Figure 6. Results of absorbance measurements of non-scattering samples, Trypan blue solution (a), and scattering samples, sample M (b), performed in the setup shown in Figure 2. The left panel shows the absorption spectrum, and the right panel shows the relationship between concentration and absorbance at the absorption peak wavelength of 591 nm. The concentrations shown in the graphs refer to the concentration of the absorbent trypan blue in both (a,b). Error bars are not shown because the standard deviation of the simulation ( 0.026) is small relative to the width of the plot points (about ±0.03). The standard deviations of the simulations are given in Appendix A.
Figure 6. Results of absorbance measurements of non-scattering samples, Trypan blue solution (a), and scattering samples, sample M (b), performed in the setup shown in Figure 2. The left panel shows the absorption spectrum, and the right panel shows the relationship between concentration and absorbance at the absorption peak wavelength of 591 nm. The concentrations shown in the graphs refer to the concentration of the absorbent trypan blue in both (a,b). Error bars are not shown because the standard deviation of the simulation ( 0.026) is small relative to the width of the plot points (about ±0.03). The standard deviations of the simulations are given in Appendix A.
Photonics 10 00754 g006
Figure 7. Histograms of the optical path length d of the detected light for: (a) the non-scattering sample (trypan blue solution); and (b) the scattering sample (sample M) measured with the setup shown in Figure 2. The region L m i d ( 0.58   cm < d     0.8   cm ) is the range of d that can be taken if the light incident through aperture S travels straight through the sample, and is included in the region L h i g h ( d > 0.8   cm ) if d is larger than it or L l o w ( 0.0   cm < d     0.58   cm ) if d is smaller than it. If the sample is not transmitted through the sample due to reflections on the sample surface, d = 0 cm. (c) illustrates the light paths in each level of the histogram.
Figure 7. Histograms of the optical path length d of the detected light for: (a) the non-scattering sample (trypan blue solution); and (b) the scattering sample (sample M) measured with the setup shown in Figure 2. The region L m i d ( 0.58   cm < d     0.8   cm ) is the range of d that can be taken if the light incident through aperture S travels straight through the sample, and is included in the region L h i g h ( d > 0.8   cm ) if d is larger than it or L l o w ( 0.0   cm < d     0.58   cm ) if d is smaller than it. If the sample is not transmitted through the sample due to reflections on the sample surface, d = 0 cm. (c) illustrates the light paths in each level of the histogram.
Photonics 10 00754 g007
Table 1. Optical properties of polystyrene microspheres (diameter: 500 nm) irradiated with 591 nm light are calculated by a numerical algorithm based on Mie scattering theory. The scattering efficiency, absorption efficiency, and anisotropy coefficient of polystyrene microspheres were calculated by a numerical calculation program based on Mie scattering theory using the wavelength of irradiated light, particle diameter and complex refractive index of polystyrene microspheres as input parameters.
Table 1. Optical properties of polystyrene microspheres (diameter: 500 nm) irradiated with 591 nm light are calculated by a numerical algorithm based on Mie scattering theory. The scattering efficiency, absorption efficiency, and anisotropy coefficient of polystyrene microspheres were calculated by a numerical calculation program based on Mie scattering theory using the wavelength of irradiated light, particle diameter and complex refractive index of polystyrene microspheres as input parameters.
Input ParametersOutput Parameters
Wavelength [nm]Diameter [nm]Refractive IndexScattering Efficiency
μ s
Absorption Efficiency
μ a
Anisotropy Efficiency
g
RealImaginary
5915001.5850.0000.83462880.00000000.8394426
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mori, A.; Yamashita, K.; Tokunaga, E. Reproduction of Visible Absorbance Spectra of Highly Scattering Suspensions within an Integrating Sphere by Monte Carlo Simulation. Photonics 2023, 10, 754. https://doi.org/10.3390/photonics10070754

AMA Style

Mori A, Yamashita K, Tokunaga E. Reproduction of Visible Absorbance Spectra of Highly Scattering Suspensions within an Integrating Sphere by Monte Carlo Simulation. Photonics. 2023; 10(7):754. https://doi.org/10.3390/photonics10070754

Chicago/Turabian Style

Mori, Ayaka, Kyohei Yamashita, and Eiji Tokunaga. 2023. "Reproduction of Visible Absorbance Spectra of Highly Scattering Suspensions within an Integrating Sphere by Monte Carlo Simulation" Photonics 10, no. 7: 754. https://doi.org/10.3390/photonics10070754

APA Style

Mori, A., Yamashita, K., & Tokunaga, E. (2023). Reproduction of Visible Absorbance Spectra of Highly Scattering Suspensions within an Integrating Sphere by Monte Carlo Simulation. Photonics, 10(7), 754. https://doi.org/10.3390/photonics10070754

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop