Understanding the behavior of light in granular media is necessary for determining the sample size, shape, and weight when probing using fiber optic setups. This is required for a correct estimate of the active pharmaceutical ingredient content in a pharmaceutical blend via near-infrared spectroscopy. Several strategies to describe the behavior of light in granular and turbid media exist. A common approach is the Monte-Carlo simulation of individual photons and their description using mean free path lengths for scattering and absorption. In this work, we chose a complementary method by approximating these parameters via real physical counterparts, i.e., the particle size, shape, and density and the resulting chord lengths. Additionally, the wavelength dependence of refractive indices is incorporated. The obtained results were compared with those obtained in an experimental setup that included the SAM-Spec Felin probe head by Indatech for detecting spatially resolved spectra of samples. Our method facilitates the interpretation of the acquired experimental results by contrasting the optical response, the physical particle attributes, and the simulation results.
U.S. Food and Drug Administration. Guidance for Industry: PAT – A Framework for Innovative Pharmaceutical Development, Manufacturing and Quality Assurance, Silver Spring, MD: U.S. Food and Drug Administration, 2004.
2.
BakeevK.A.Process Analytical Technology: Spectroscopic Tools and Implementation Strategies for the Chemical and Pharmaceutical Industries, Oxford, UK: Wiley, 2010.
3.
O’DonnellC.P.FaganC.CullenP.J.Process Analytical Technology for the Food Industry, New York, NY: Springer Science and Business Media, 2014.
4.
ChewW.SharrattP.“Trends in Process Analytical Technology”. Anal. Methods. 2010. 2(10): 1412.
5.
SimonL.L.PatakiH.MarosiG.MeemkenF.et al.“Assessment of Recent Process Analytical Technology (PAT) Trends: A Multiauthor Review”. Org. Process Res. Dev. 2015. 19(1): 3–62.
6.
JamrógiewiczM.“Application of the Near-infrared Spectroscopy in the Pharmaceutical Technology”. J Pharm. Biomed. Anal. 2012. 66: 1–10.
O. Scheibelhofer, D.M. Koller, P. Kerschhaggl, J.G. Khinast. “Continuous Powder Flow Monitoring via Near-infrared Hyperspectral Imaging”. In: IEEE 2012 International Instrumentation and Measurement Technology Conference Proceedings. Graz, May 13–16, 2012. Pp. 748–753. doi: 10.1109/I2MTC.2012.6229431.
9.
KesslerR.W.“Perspectives in Process Analysis”. J. Chemom. 2013. 27(11): 369–378.
10.
ScheibelhoferO.BalakN.KollerD.M.KhinastJ.G.“Spatially Resolved Monitoring of Powder Mixing Processes via Multiple NIR-Probes”. Powder Technol. 2013. 243: 161–170.
11.
WestadF.SchmidtA.KermitM.“Incorporating Chemical Band-Assignment in Near Infrared Spectroscopy Regression Models”. J. Near Infrared Spectrosc. 2008. 16(3): 265–273.
12.
CiurczakE.W.IgneB.Pharmaceutical and Medical Applications of Near-Infrared Spectroscopy, Boca Raton, FL: CRC Press, 2015.
13.
PasikatanM.C.SteeleJ.L.SpillmanC.K.HaqueE.“Near Infrared Reflectance Spectroscopy for Online Particle Size Analysis of Powders and Ground Materials”. J. Near Infrared Spectrosc. 2001. 9(1): 153.
14.
OtsukaM.MouriY.MatsudaY.“Chemometric Evaluation of Pharmaceutical Properties of Antipyrine Granules by Near-infrared Spectroscopy”. AAPS PharmSciTech. 2003. 4(3): E47.
15.
ScheibelhoferO.GrabnerB.Bondi JrR.W.IgneB.et al.“Designed Blending for Near Infrared Calibration”. J. Pharm. Sci. 2015. 104(7): 2312–2322.
16.
DahmD.J.DahmK.D.“The Physics of Near Infrared Scattering”. In: WilliamsP.NorrisK. (eds). Near Infrared Technology in the Agricultural and Food Industry, St. Paul, MN: American Association of Cereal Chemists, 2004. Chap. 1, Pp. 1–17.
17.
KesslerW.OelkrugD.KesslerR.“Using Scattering and Absorption Spectra as MCR-hard Model Constraints for Diffuse Reflectance Measurements of Tablets”. Anal. Chim. Acta. 2009. 642(1–2): 127–134.
18.
ShiZ.AndersonC.A.“Pharmaceutical Applications of Separation of Absorption and Scattering in Near-Infrared Spectroscopy (NIRS)”. J. Pharm. Sci. 2010. 99(12): 4766–4783.
KokhanovskyA.A.Light Scattering Reviews - Single and Multiple Light Scattering, Chichester, UK: Praxis Publishing, 2006.
21.
PoludniowskiG.G.EvansP.M.“Optical Photon Transport in Powdered-Phosphor Scintillators. Part 1. Multiple-scattering and Validity of the Boltzmann Transport Equation”. Med. Phys. 2013. 40(4): 041904.
22.
HohmannF.VoitJ.Schäferet al.“Comparison of Monte Carlo Simulations with Exact Maxwell Solutions for Polarized Light Scattering by Multiple Absorbing Spheres”. J. Phys. Conf. Ser. 2012. 369: 012007.
23.
BohrenC.F.HuffmanD.R.Absorption and Scattering of Light by Small Particles, New York, NY: Wiley-VCH Verlag GmbH and Co. KG, 2004.
24.
XuY.L.GustafsonB.Å.S.“A Generalized Multiparticle Mie-solution: Further Experimental Verification”. J. Quant. Spectrosc. Radiat. Transfer. 2001. 70(4–6): 395–419.
25.
JacquesS.L.“Optical Properties of Biological Tissues: A Review”. Phys. Med. Biol. 2013. 58(11): R37–61.
26.
KervellaM.D’AbzacF.-X.HacheF.HespelL.et al.“Picosecond Time Scale Modification of Forward Scattered Light Induced by Absorption Inside Particles”. Opt Express. 2012. 20(1): 32–41.
27.
RogersJ.D.RadosevichA.J.YiJ.BackmanV.“Modeling Light Scattering in Tissue as Continuous Random Media Using a Versatile Refractive Index Correlation Function”. IEEE J Sel Topics in Quantum Electron. 2014. 20(2): 7000514.
28.
WangL.JacquesS.L.ZhengL.“MCML–Monte Carlo Modeling of Light Transport in Multi-Layered Tissues”. Comput. Meth. Programs Biomed. 1995. 47: 131–146.
29.
D.K. McAllister, A. Lastra, W. Heidrich. “Efficient Rendering of Spatial Bi-directional Reflectance Distribution Functions”. In: Proceedings of the ACM SIGGRAPH/EUROGRAPHICS Conference on Graphics Hardware. Saarbrücken: Eurographics Association, 2002. Pp. 79–89.
30.
GovaertsY.M.VerstraeteM.M.“Raytran: A Monte Carlo Ray-tracing Model to Compute Light Scattering in Three-Dimensional Heterogeneous Media”. IEEE Trans. Geosci. Remote Sens. 1998. 36(2): 493–505.
31.
RitschelT.DachsbacherC.GroschT.KautzJ.“The State of the Art in Interactive Global Illumination”. Comput. Graph. Forum. 2012. 31(1): 160–188.
32.
GoedeckeG.H.“Radiative Transfer in Closely Packed Media”. J. Opt. Soc. Am. 1977. 67(10): 1339–1348.
33.
KubelkaP.“New Contributions to the Optics of Intensely Light-Scattering Materials Part I”. J. Opt. Soc. Am. 1948. 38(5): 448.
34.
ZhangY.McGeorgeG.“Quantitative Analysis of Pharmaceutical Bilayer Tablets Using Transmission Raman Spectroscopy”. J. Pharm. Innovation. 2015. 10(3): 269–280.
35.
BerntssonO.DanielssonL.-G.FolestadS.“Estimation of Effective Sample Size when Analysing Powders with Diffuse Reflectance Near-Infrared Spectrometry”. Anal. Chim. Acta. 1998. 364(1–3): 243–251.
36.
RoggoY.ChalusP.MaurerL.Lema-MartinezC.et al.“A Review of Near Infrared Spectroscopy and Chemometrics in Pharmaceutical Technologies”. J. Pharm. Biomed. Anal. 2007. 44(3): 683–700.
37.
KanickS.C.RobinsonD.J.SterenborgH.J.C.M.AmelinkA.“Monte Carlo Analysis of Single Fiber Reflectance Spectroscopy: Photon Path Length and Sampling Depth”. Phys. Med. Biol. 2009. 54(22): 6991–7008.
38.
PapaioannouT.PreyerN.W.FangQ.BrightwellA.et al.“Effects of Fiber-optic Probe Design and Probe-to-Target Distance on Diffuse Reflectance Measurements of Turbid Media: An Experimental and Computational Study at 337 nm”. Appl. Opt. 2004. 43(14): 2846–2860.
39.
ArhaliassA.SnabreP.“Multiple Light Scattering in Random Systems: Analysis of the Backscattering Spot Image”. Eur. Phys. J. Appl. Phys. 1999. 6(1): 81–87.
40.
ShiZ.AndersonC.A.“Application of Monte Carlo Simulation-Based Photon Migration for Enhanced Understanding of Near-Infrared (NIR) Diffuse Reflectance. Part I: Depth of Penetration in Pharmaceutical Materials”. J. Pharm. Sci. 2010. 99(5): 2399–2412.
41.
AbrahamssonC.JohanssonJ.Andersson-EngelsS.SvanbergS.et al.“Time-Resolved NIR Spectroscopy for Quantitative Analysis of Intact Pharmaceutical Tablets”. Anal. Chem. 2005. 77(4): 1055–1059.
42.
PilorgetC.VincendonM.PouletF.“A Radiative Transfer Model to Simulate Light Scattering in a Compact Granular Medium using a Monte-Carlo Approach: Validation and First Applications”. J. Geophys. Res.: Planets. 2013. 118(12): 2488–2501.
43.
OkamotoT.KumagawaT.MotodaM.IgarashiT.et al.“Monte Carlo Simulation of Light Reflection from Cosmetic Powder Particles Near the Human Skin Surface”. J. Biomed. Opt. 2013. 18(6): 61232.
44.
TranN.T.CampbellC.G.ShiF.G.“Study of Particle Size Effects on an Optical Fiber Sensor Response Examined with Monte Carlo Simulation”. Appl. Opt. 2006. 45(29): 7557–7566.
45.
SadjadiZ.MiriM.“Diffusion of Light in Two-dimensional Granular Materials”. Am. Inst. Phys. Conf. Proc. 2013. 1542: 109–112.
46.
SvenssonT.VynckK.GrisiM.SavoR.et al.“Holey Random Walks: Optics of Heterogeneous Turbid Composites”. Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 2013. 87(2): 022120.
47.
SvenssonT.VynckK.AdolfssonE.FarinaA.et al.“Light Diffusion in Quenched Disorder: Role of Step Correlations”. Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 2014. 89(2): 022141.
48.
KostinskiA.B.“On the Extinction of Radiation by a Homogeneous but Spatially Correlated Random Medium”. J. Opt. Soc. Am. A. 2001. 198(8): 1929–1933.
49.
OlsonG.L.MillerD.S.LarsenE.W.MorelJ.E.“Chord Length Distributions in Binary Stochastic Media in Two and Three Dimensions”. J. Quant. Spectrosc. Radiat. Transfer. 2006. 101(2): 269–283.
50.
DlugachJ.M.MishchenkoM.I.LiuL.MackowskiD.W.“Numerically Exact Computer Simulations of Light Scattering by Densely Packed, Random Particulate Media”. J. Quant. Spectrosc. Radiat. Transfer. 2011. 112(13): 2068–2078.
51.
TishkovetsV.P.PetrovaE.V.MishchenkoM.I.“Scattering of Electromagnetic Waves by Ensembles of Particles and Discrete Random Media”. J. Quant. Spectrosc. Radiat. Transfer. 2011. 112(13): 2095–2127.
KienleM.S.Patterson. “Determination of the Optical Properties of Turbid Media from a Single Monte Carlo Simulation”. Phys. Med. Biol. 1996. 41(10): 2221–2227.
54.
WenX.ZhongX.YuT.ZhuD.“A Monte Carlo Based Lookup Table for Spectrum Analysis of Turbid Media in the Reflectance Probe Regime”. Quant. Electron. 2014. 44(7): 641–645.
55.
YeongC.L.Y.TorquatoS.“Reconstructing Random Media”. Phys. Rev. E. 1998. 57(1): 495–506.
56.
VianaM.JouanninP.PontierC.ChuliaD.“About Pycnometric Density Measurements”. Talanta. 2002. 57(3): 583–93.
57.
F. Chauchard, S. Roussel. Sonde spectroscopique et procédé pour détecter une inhomogénéité. EP 2435816 B1. WO2010136732 A1. Filed 2010. Issued 2010.
58.
DehghaniH.LeblondF.PogueB.W.ChauchardF.“Application of Spectral Derivative Data in Visible and Near-Infrared Spectroscopy”. Phys. Med. Biol. 2010. 55(12): 3381–3399.
59.
F. Chauchard, J. Revel, M. Zeaiter. “SAM-Spec® Spectroscopy used to Measure Simultaneously API Concentration and Particle Size Distribution Example of Application on Drug Layering Process”. In: Chimiométrie 2009 Conference. Paris, Nov 30–Dec 1, 2009. Pp.
60.
BoiretM.ChauchardF.“Use of Near-Infrared Spectroscopy and Multipoint Measurements for Quality Control of Pharmaceutical Drug Products”. Anal. Bioanal. Chem. 2017. 409(3): 683–691.
61.
ShiZ.AndersonC.A.“Scattering Orthogonalization of Near-Infrared Spectra for Analysis of Pharmaceutical Tablets”. Anal. Chem. 2009. 81(4): 1389–1396.
62.
KattawarG.“A Three-Parameter Analytic Phase Function for Multiple Scattering Calculations”. J. Quant. Spectrosc. Radiat. Transfer. 1975. 15: 839–849.
63.
G. Fournier, J.L. Forand. “Analytic Phase Function for Ocean Water”. In: J.S. Jaffe, editor. Ocean Optics XII. Proc. SPIE. 1994. 2258: pp. 194–201. doi :10.1117/12.190063.
I. Sassi, M. Khemiri. “Scattering of Light from Rough Surfaces, the Limits of Validity of Geometric Optics Approximation Method”. In: Progress in Electromagnetics Research Symposium Proceedings. Taipei, 2013. Pp. 623–625.
66.
WatanabeA.“A Trial Production of a Table of the Optical Crystallographic Characteristics of Crystalline Drugs Including Crystal Habits (Study of Crystalline Drugs by Means of a Polarizing Microscope. XIX)”. Yakugaku Zasshi. 2002. 122(8): 595–606.
67.
SultanovaN.KasarovaS.NikolovI.“Dispersion Properties of Optical Polymers”. Acta Phys. Pol. A. 2009. 116(4): 585–587.
MazzoloA.RoesslingerB.“Monte-Carlo Simulation of the Chord Length Distribution Function across Convex Bodies, Non-Convex Bodies and Random Media”. Monte Carlo Methods and Applications. 2004. 10(3–4): 443–454.
70.
RufA.WorlitschekJ.MazzottiM.“Modeling and Experimental Analysis of PSD Measurements through FBRM”. Part. Part. Syst. Char. 2000. 17(4): 167–179.
71.
HansenS.“Estimation of Chord Length Distributions from Small-Angle Scattering Using Indirect Fourier Transformation”. J. Appl. Crystallogr. 2003. 36(5): 1190–1196.
72.
P.D. Bradley. The Development of a Novel Silicon Microdosimeter for High LET Radiation Therapy. [Ph.D. Thesis]. Wollongong, Australia: University of Wollongong, 2000.
73.
StoyanD.WagnerA.HermannH.ElsnerA.“Statistical Characterization of the Pore Space of Random Systems of Hard Spheres”. J. Non-Cryst. Solids. 2011. 357(6): 1508–1515.
74.
ReboulN.VincensE.CambouB.“A Statistical Analysis of Void Size Distribution in a Simulated Narrowly Graded Packing of Spheres”. Granular Matter. 2008. 10(6): 457–468.
75.
A. Mazzolo. “Chord Length Distributions in Dense Stochastic Media”. In: First International Conference on Physics and Technology of Reactors and Applications. Marrakech, 2007.
76.
JiW.MartinW.R.“Determination of Chord Length Distributions in Stochastic Media Composed of Dispersed Microspheres”. Trans. Am. Nucl. Soc. 2007. 96: 467–469.
77.
ShkuratovY.G.GrynkoY.S.“Light Scattering by Media Composed of Semitransparent Particles of Different Shapes in Ray Optics Approximation: Consequences for Spectroscopy, Photometry, and Polarimetry of Planetary Regoliths”. Icarus. 2005. 173(1): 16–28.
78.
ScheibelhoferO.SacherS.KhinastJ.G.“Near Infrared Light Spectroscopy: Applicability and New Developments”. European Pharmaceutical Review. 2014, 1: 51–55.
79.
FooteP.D.“Illumination from a Radiating Disk”. Bulletin of the Bureau of Standards. 1915. 12: 583–586.
80.
WeberM.J.Handbook of Laser Science and Technology–Supplement 2: Optical Materials, Boca Raton, FL: CRC Press, 1995.
81.
KogelnikH.“On the Propagation of Gaussian Beams of Light Through Lenslike Media Including those with a Loss or Gain Variation”. Appl. Opt. 1965. 4(12): 1562.
82.
LinL.ZhangM.“Evaluation of Scaling Monte Carlo Methods for Backscattering Properties of Turbid Media with Gaussian Incidence”. Opt. Commun. 2015. 336: 134–139.
K.F. Eckerman, J.C. Ryman, A.C. Taner, G.D. Kerr. “Traversal of Cells by Radiation and Absorbed Fraction Estimates for Electrons and Alpha Particles”. In: A.T. Schlafke-Stelson, E.E. Watson, editor. Proceedings of the Fourth International Radiopharmaceutical Dosimetry Symposium. Oak Ridge, TN: Oak Ridge Associated Universities, 1986. Pp. 26–36.
85.
SjöstrandN.“What is the Average Chord Length?”. Ann. Nucl. Energy. 2002. 29(13): 1607–1608.
86.
J.-P. Schäfer. Implementierung und Anwendung analytischer und numerischer Verfahren zur Lösung der Maxwellgleichungen für die Untersuchung der Lichtausbreitung in biologischem Gewebe. [Ph.D. Thesis]. Ulm, Germany: University of Ulm, 2011.
87.
SchäferJ.LeeS.C.KienleA.“Calculation of the Near Fields for the Scattering of Electromagnetic Waves by Multiple Infinite Cylinders at Perpendicular Incidence”. J. Quant. Spectrosc. Radiat. Transfer. 2012. 113(16): 2113–2123.
88.
GombásÁ.AntalI.Szabó-RévészP.MartonS.et al.“Quantitative Determination of Crystallinity of Alpha-Lactose Monohydrate by Near Infrared Spectroscopy (NIRS)”. Int. J. Pharm. 2003. 256(1–2): 25–32.
89.
AenuguH.P.Kumar, SrisudharsonD.S.ParthibanN.et al.“Near Infra Red Spectroscopy - An Overview”. Ing. J. ChemTech Res. 2011. 3(2): 825–836.
90.
SekkalM.LegrandP.“A Vibrational Molecular Force Field of Model Compounds with Biological Interest - III. Harmonic Dynamics of α- and β-D-galactose in the Crystalline State”. Spectrochim. Acta. 1992. 48(7): 959–973.
91.
GaffourH.M.BouterfasM.BekhtiN.DerrarS.N.et al.“Harmonic Dynamics of α-D-Lactose in the Crystalline State”. J. Mol. Imaging Dyn. 2011. 01(01): 1000102.
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
