
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

39237548
70954
10.1038/s41598-024-70954-x
Article
Twists through turbidity: propagation of light carrying orbital angular momentum through a complex scattering medium
Khanom Fatima 1
Mohamed Nawal 1
Lopushenko Ivan 2
Sdobnov Anton 2
Doronin Alexander 3
Bykov Alexander 2
Rafailov Edik 1
Meglinski Igor i.meglinski@aston.ac.uk

1
1 https://ror.org/05j0ve876 grid.7273.1 0000 0004 0376 4727 College of Engineering and Physical Sciences, Aston University, Birmingham, B4 7ET UK
2 https://ror.org/03yj89h83 grid.10858.34 0000 0001 0941 4873 Optoelectronics and Measurement Techniques, University of Oulu, P.O. Box 4500, 90014 Oulu, Finland
3 https://ror.org/0040r6f76 grid.267827.e 0000 0001 2292 3111 School of Engineering and Computer Science, Victoria University of Wellington, Wellington, 6140 New Zealand
5 9 2024
5 9 2024
2024
14 2066211 4 2024
22 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
We explore the propagation of structured vortex laser beams-shaped light carrying orbital angular momentum (OAM)-through complex multiple scattering medium. These structured vortex beams consist of a spin component, determined by the polarization of electromagnetic fields, and an orbital component, arising from their spatial structure. Although both spin and orbital angular momenta are conserved when shaped light propagates through a homogeneous, low-scattering medium, we investigate the conservation of these angular momenta during the propagation of Laguerre–Gaussian (LG) beams with varying topological charges through a turbid multiple scattering environment. Our findings demonstrate that the OAM of the LG beam is preserved, exhibiting a distinct phase shift indicative of the ‘twist of light’ through the turbid medium. This preservation of OAM within such environments is confirmed by in-house developed Monte Carlo simulations, showing strong agreement with experimental studies. Our results suggest exciting prospects for leveraging OAM in sensing applications, opening avenues for groundbreaking fundamental research and practical applications in optical communications and remote sensing.

Keywords

Vortex beams
Orbital angular momentum
Light propagation
Turbid scattering medium
Twist of light
Subject terms

Optical physics
Optical physics
501100002341 Academy of Finland (Suomen Akatemia) issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Navigating through complex media presents significant challenges in optical imaging and information transmission, primarily due to the erratic scattering of light1. Nano-scale spatial and/or temporal variations in the refractive index within complex random, heterogeneous, emulsions, powders, porous, polycrystalline materials and fractal media, such as biological tissue, multi-mode fibers, and atmospheric phenomena like clouds and fog lead to optical scattering, which distorts the wavefront of incident light2. The endeavor to elucidate, comprehend, and ultimately forecast the propagation of waves presents a formidable set of challenges to researchers within this domain. Since the inception of this field in the early twentieth century, it became apparent that the introduction of novel concepts and significant approximations was imperative. Noteworthy among these conceptual advancements are the effective medium theory3,4 and the radiative transport theory5,6, which remain pivotal to contemporary investigations. A principal challenge associated with these approximative frameworks arises from the often comparable scales of inhomogeneities within the complex medium to the wavelength of interest. These theories are primarily applicable in scenarios where the inhomogeneities’ length scales significantly exceed the wavelength. Consequently, numerous relevant cases emerge where the effective medium theory, radiative transport theory, or both may not adequately apply, leading to considerable discrepancies in their predictive capabilities.

In recent years, research efforts have concentrated on mitigating the impact of scattering to increase the penetration depths of light in turbid media7. Significant progress has been made in the understanding and application of Optical Reciprocity (OR) principles8. This advancement has facilitated the development of innovative techniques for focusing light deeply within tissue-like scattering media9,10 and unveiled the OR relationship between incident probing light and back-scattered light in turbid tissue-like medium.

Wavefront shaping emerges as a promising technique, offering precise control over the phase and amplitude of light waves. By manipulating light wavefronts, it becomes possible to shape the spatial distribution of light intensity effectively11. Recent advancements in optical wavefront shaping and phase recording techniques have significantly enhanced our capability to overcome the challenges posed by random light scattering in turbid media12. Traditional methodologies, including the diffusion equation, frequently failed to account for the complex interference patterns of scattered light, thereby exacerbating the challenge. Nevertheless, the advent of innovative approaches, such as optical wavefront shaping, precise phase recording, the employment of structured light13,14, and the application of complex optics, have effectively surmounted these historical limitations. These developments represent a paradigm shift in our approach to manipulating and controlling light in highly scattering environments, opening new avenues for research and application in fields ranging from biomedical imaging to optical communication1. These breakthroughs enable researchers to exert control over the coherent transport of light within highly scattering media, facilitating significant progress in various applications. Improved focusing and transmission of light through complex systems and the performance of intricate tasks within turbid environments, such as optical micro-manipulation, are now achievable. The capability to manipulate and control light waves in turbid environments marks a significant achievement in the field of optical engineering, heralding a wide range of practical applications and scientific breakthroughs. Techniques such as wavefront shaping, phase-conjugation, and the employment of spatial light modulators (SLMs) are instrumental in effectively directing the propagation of light15. The shaped light with a helical wavefront, possessing Orbital Angular Momentum (OAM)16 and also referred to as vortex beams (VBs)17, presents new unique opportunities in the field of biological imaging.

The OAM carried by vortex beams represents a novel set of carrier signals distinct from amplitude, phase, polarization, and frequency. When combined with traditional multiplexing methods, OAM can significantly enhance channel transmission capacity. Research has revealed that partially coherent vortex light fields, possessing both OAM and partial coherence, offer unique benefits in applications such as beam shaping, ghost imaging, optical communications, and information encryption13,18. Recent research underscores that OAM beams achieve superior penetration depths through complex, tissue-like scattering media compared to conventional Gaussian beams16,19–21. This attribute is highly advantageous for biomedical applications, where deeper penetration is crucial for achieving enhanced imaging resolution within deep tissue structures.

In current paper, we investigate the predictive capabilities of OAM shaped light associated with Laguerre-Gaussian (LG) beams22. This is achieved by analyzing the alterations in their helical wavefront as they propagate through a complex scattering medium, utilizing a specially developed Monte Carlo (MC) method presented below. The study encompasses various wave properties of light, including coherence, polarization, interference, and the phenomena of reflection and refraction at the boundary of the medium. These properties facilitate a comprehensive interpretation of light behavior within complex tissue-like scattering medium.

Results

For LG beams, we rely on the following representation valid in paraxial approximation23–25:1 LGpℓ(ρ,ϕ,z)=2p!π(|ℓ|+p)!w2(z)2ρw(z)|ℓ|Lp|ℓ|2ρ2w2(z)exp-ρ2w2(z)××exp[i(2p+|ℓ|+1)arctan(z/zR)]exp-ikρ2z2(z2+zR2)exp[-iℓφ]exp[-ikz].

Here, k=2π/λ, w(z)=w(0)1+(z/zR)2, zR=πw2(0)/λ, w(0) corresponds to the zero-order Gaussian beam waist and is adjusted to fit experimental image, {ρ,ϕ,z} is cylindrical coordinate system with beam propagating along z axis. Then phase evolution of the LG beam is defined as25:2 Ψ(ρ,ϕ,z)=argLGpℓ(ρ,ϕ,z)=-kρ2z2(z2+zR2)-ℓϕ-kz+G(z),

where G(z)=(2p+|ℓ|+1)arctan(z/zR) corresponds to the Gouy phase. The helical phase of LG beam, characterized by a phase front that spirals around the beam’s axis, results from the azimuthal phase term in the beam’s electric field expression. This term correlates with the beam’s topological charge (ℓ and/or radial index (p), defining the beam’s unique spatial structure.

The results of computational MC modeling shows that during transmission, the beam not only loses intensity due to phenomena like scattering, refraction, and absorption but also experiences beam spreading (Fig. 1). As one can see the speckle interference pattern, resulting from the superposition of helical wavefront components, exhibits spatial variations both in intensity and phase distributions. Such patterns lead to the distortion of the LG beam’s doughnut structure even in low-scattering (d/l∗=2.5) environments (see Fig. 1-(top)) and its complete breakdown in environments characterized by higher scattering (d/l∗=5).

Here, the complex turbid media, distinguished by either low or multiple light scattering, are defined by their optical depth (d/l∗), which quantitatively represents the degree of light attenuation and scattering intensity as it travels through the medium; where, d denotes the medium’s thickness or the depth of the beam penetration, and l∗ is the transport mean free path (l∗=1/μs, where μs is the scattering coefficient of the medium). It is generally accepted that d/l∗ of 5–6 or greater indicates a medium undergoing multiple or diffuse scattering, while values between approximately 2 to 5 suggest single or intermediate scattering, characterized as ‘snake-like photons’26.

In contrast, in medium with low scattering (d/l∗∼2), the phase of LG beam progresses with minimal distortion, preserving its original OAM state, petal-like helical phase structure (see Fig. 1-(bottom). Conversely, in higher scattering (d/l∗=5), the LG beam’s phase structure diffusively spreads, and its helical phase front is disrupted, leading to the formation of a complex phase speckle pattern (refer to Fig. 2).Fig. 1 Intensity (top) and phase (bottom) distributions of the LG05 beam along propagation from Spatial Light Modulator (SLM) into the complex scattering medium. Insets show, respectively, 2D spatial distributions of intensity (normalized) and phase at the surface of the medium (d/l∗=0) and within the medium at the depth d/l∗=2.5 and d/l∗=5.0.

The numerical OAM sorter employed to confirm presence of the topological charge in the scattered signal is adopted from27. It simulates optical elements that, firstly, transform the attenuated helical light into beam with a transverse phase gradient and, secondly, focus the OAM state present in the resulting beam to a specific lateral position defined by the value of topological charge ℓ. This method, originally introduced and experimentally confirmed to efficiently distinguish between OAM states with different charges, exhibits certain potential to detect OAM presence in biomedical diagnostic applications.Fig. 2 Intensity (a) and phase (b) distributions of the LG03 beam after propagation through a complex scattering medium with 1 mm thickness and specified scattering coefficients μs=2,4,6,10mm-1; absorption coefficient μa=0.01mm-1 and anisotropy factor g=0.8. Row (c) depicts intensity output of the OAM sorter algorithm adopted from27. Dotted red line on left figure indicates position of the intensity peak which corresponds to the Gaussian beam without OAM27. The scale bar in (a) represents 0.5 mm, which applies to all images.

Figure 2 demonstrates that the numerical OAM sorter27 can detect the OAM memory—the retention of the helical phase structure by the LG beam as it propagates through the scattering medium. More specifically we apply a method for efficiently OAM states of light by employing two static optical elements that transform the helical phase structure of OAM beams into a linear phase gradient27. This transformation is achieved through a Cartesian to log-polar coordinate mapping, followed by a phase correction to eliminate distortions. The transformed LG beams are then focused to distinct lateral positions corresponding to their OAM states. The horizontal position of the intensity blob(s) on the detector plane (see Fig. 2) is a direct indicator of the OAM state of the detected light. Each unique OAM state results in a distinct and predictable lateral position, allowing for efficient sorting and detection of multiple OAM states simultaneously.

The preservation of the OAM phase, known as OAM phase memory, is ascribed to the congruence in the alterations of the spiral trajectories encircling the beam axis, arising from rotational symmetry. The spiral trajectories undergo substantial modifications from their initial LG beam configurations due to the interaction with the scattering medium. This leads to disparate partial components of the LG beam’s helical wavefront undergoing varying extents of phase distortions. Owing to the inherent rotational symmetry, these alterations in the spiral trajectories exhibit uniformity in both magnitude and direction around the beam’s axis. Consequently, despite the significant phase distortions resulting from multiple scattering events (d/l∗∼10), the overall retention of the LG beam’s OAM is distinctly evident in the phase behavior of the speckle pattern. This pattern mirrors the modulation of the LG beam’s initial phase as configured by the SLM (see the Supplementary Video).

The results of this phenomenological model are well agreed with the experimental results28, presented in Fig. 3. The LG beam preserves its initial phase structure and spatial intensity profile, demonstrating its ability to maintain phase coherence and OAM characteristics in a medium with moderate (d/l∗=2) scattering properties (see Fig. 3a). Propagation through the multiple-scattering medium (d/l∗∼10) caused a diffusive spread and disruption of the LG beam’s helical phase front, resulting in a complex speckle pattern (see Fig. 3b). Despite significant phase distortions, the phase speckle pattern retained initial modulation, indicating partial OAM memory, especially within the axial annular region (see Fig. 3c,d, respectively).Fig. 3 Experimentally observed phase distributions for the LG beam propagating through (a) low-scattering (d/l∗=2) and (b) multiple-scattering (d/l∗=9.6) media. The central annular zone (highlighted by contours) corresponds to the LG03 beam as if it were passing through a non-scattering medium. Phase change at the selected area of the LG03 axial annular area for low-scattering (c), multiple-scattering (d) media and a comparison of relative phase change (e) according to the initial phase configuration (-3π/10⩽Ψ⩽3π/10) at the SLM (f).

Thus, despite the prevalence of strong diffuse scattering, the phase within the speckle pattern retains modulation reflective of the LG beam’s initial phase, illustrating the persistence of OAM’s memory effect amidst multiple scattering.

Discussion

This study embarked on an exploration of the propagation dynamics of LG vortex beams carrying OAM through turbid, complex scattering media. The results presented above illustrate the impact of scattering on the phase preservation of an LG beam as it propagates through different scattering media. In the low-scattering medium (d/l∗=2), the LG beam maintains its helical phase structure, with minimal disruption to its initial phase and intensity profile (see Fig. 3a). This demonstrates that the LG beam can retain its OAM characteristics and phase coherence in environments with moderate scattering properties. Conversely, in the multiple-scattering medium (d/l∗=9.6), the LG beam undergoes significant phase distortion, resulting in a complex speckle pattern (see Fig. 3b). Despite these distortions, the speckle pattern retains some modulation of the initial phase, indicating a degree of OAM memory. This partial preservation of the OAM phase is particularly noticeable within the axial annular region, while it diminishes in the central and outer regions beyond this annular zone (see Fig. 3).

Utilizing both experimental methods and in-house developed MC simulations, we investigated the conservation of OAM and the accompanying phase shifts, indicative of the unique ‘twist of light’ phenomena through these challenging environments. Our findings demonstrate a robust preservation of OAM despite the multiple scattering processes, underpinning the potential of OAM for enhanced sensing capabilities in optical communications and remote sensing. The employment of a numerical OAM sorter further validated the persistence of OAM in scattered signals, affirming the foundational concept of OAM phase memory within the scattering media.Fig. 4 Comparison between simulated and experimentally measured phase patterns for the LG03 beam propagated through a complex scattering medium μs=4mm-1, μa=0.05mm-1, g=0.8.

The investigation reveals that despite the inherent scattering challenges presented by complex media, OAM’s unique properties remain largely intact, exhibiting only minor distortions even in high-scattering environments. This resilience of OAM through turbid media opens new avenues for advanced sensing applications, offering a novel paradigm for optical communication systems that require high penetration depths and robust signal integrity. The successful implementation of MC simulations aligns closely with experimental observations (Fig. 4), providing a powerful tool for predicting light propagation behaviors in similarly complex environments. This study contributes to the fundamental understanding of light-matter interactions within turbid medium, as well as paves the way for the practical application of OAM in diverse fields such as biomedical imaging, atmospheric sensing, and secure communications. Future work will aim to refine the predictive capabilities of our models and explore the potential of multiplexed OAM states for even more complex communication and sensing tasks, thereby leveraging the full spectrum of opportunities afforded by the unique properties of structured light beams.

Methodology

Experimental setup

A Mach–Zehnderinterferometer29,30, typically used to characterize the relative phase shift between two laser beams, allows to observe the phase structure of the LG beam interfered with the reference plane wave of the same frequency31. In this study, the modified Mach–Zehnder-based interferometer32 is used to examine an evolution of OAM of the LG beams propagated through the medium (see Fig.5).Fig. 5 LD—laser diode. P—polarizer. FM—fiber mount. BC—beam colimator. SLM—spatial light modulator. M1, M2—mirrors. L1, L2, L3, L4—lenses. PH—pinhole. PBS—polarizing beam splitter. HWP—half wave plate. S—cuvette filled with the sample liquid. BS—beam splitter. NF—neutral filter. O—objective. CCD—camera. The detailed description of the optical setup is presented in the main article text.

In the experiment the coherent laser source (60 mW, OBIS, Coherent, USA) emitting at 633 nm is used as a light source. To clear the optical mode, laser beam was focused into a single-mode optical fiber (Thorlabs, USA) F. The output beam was collimated using beam collimator (Thorlabs, USA) BC to obtain Gaussian beam with 1.6 mm waist diameter. Also, to obtain the horizontal linear polarization the polarizer (Thorlabs, USA) P was been used after collimator. Further, the Gaussian beam has been split in sample and reference beams using polarizing beam splitter (Thorlabs, USA) PBS. The sample beam illuminated the phase only spatial light modulator (PLUTO-2-NIR-011, Holoeye, Germany) SLM operating in reflective regime. To produce LG beam with different moments the corresponding forked diffraction patterns were generated at SLM. The light diffracted from SLM has been directed using set of mirrors M to the lens L3 (f=45mm, Thorlabs, USA). This lens was used to focus the first order diffraction through pinhole (Thorlabs, USA). Further the LG beam was re-collimated using lens L4 (f=45mm, Thorlabs, USA).

Finally, sample beam was passed through the custom made turbid tissue-like media of low scattering: d/l∗=2 (thickness of d=1mm and μs=10mm-1) and multiple-scattering: d/l∗∼10 (thickness of d=8mm and μs=6mm-1). The anisotropy of scattering in both cases was g=0.8. The detailed methodology for the preparation of such phantoms is thoroughly described in the study by Sieryi et al., as referenced in33. The reference beam was passed through the half-wave plate (Thorlabs, USA) to control the polarization orientation of Gaussian beam. Further, reference beam was expanded by set of lenses L1 (f=30mm Thorlabs, USA) and L2 (f=70mm, Thorlabs, USA) and directed to the beam splitter (Thorlabs, USA) BS where expanded Gaussian beam interfered with LG beam. Finally, the interference pattern was registered using CMOS camera (DCC3240M, 1280×1024, Thorlabs, USA) CAM in combination with objective (10×, Nikon, Japan) O. The proposed setup allowed to obtain interference patterns in both on-axis34 and off-axis35 regimes. The corresponding simulated intensity profiles for Gaussian beam, expanded Gaussian beam, LG05 beam and corresponding interference pattern for on-axis regime are shown in “jet” colour palette at Fig.5. Black arrows shows the polarization direction. The corresponding simulated phases for LG05 beam at the SLM, before the cuvette with sample liquid and after the cuvette with sample liquid are show in “parula” colour palette at Fig. 5.

Simulation

MC method is vastly used as an effective tool for imitation of image transfer through turbid scattering media36, including complex structured ones like biological tissues37 where the technique is considered as ‘gold standard’38,39. Recently, a number of various MC methods have been developed to simulate propagation of polarized light within turbid tissue-like scattering media40–43. The MC approach, informed by the exact analytic Milne solution and utilizing the iterative solution to the Bethe–Salpeter equation44–46 with Jones vector formalism, effectively tracks the polarization of MC-photons within turbid tissue-like media and simulates coherent backscattering47,48. Advantages of the Bethe–Salpeter-based approach involve a direct relation to the analytic Milne solution49 and intuitive physical interpretation of multiple scattering process via ladder diagrams. The approach has subsequently progressed to include polarized light, extending the computational capabilities of biomedical optics diagnostic toolbox50. Fundamental ground for the polarized MC approaches was established by the Vector Radiative Transfer Equation (VRTE) which represents a system of equations for each Stokes parameter and can be rigorously derived from the Maxwell electromagnetic theory51–53. Recently, it has been shown on the fundamental level that VRTE and Bethe–Salpeter-based approaches are equivalent under certain conditions54. Modern implementations of the polarization-resolved MC50 aim to provide a comprehensive description of polarized light scattering with either Jones or Mueller formalism, depending on the representation of the polarization state55. Recent biomedical polarimetry advancements highlight the efficacy of circularly polarized light (CPL) in characterizing tissues with optical anisotropy56–58. Accurate simulation tools are essential for studying CPL-tissue interactions59,60.

The inclusion of shaped light carrying OAM in MC modeling has garnered attention for its potential to enhance the biomedical diagnosis toolkit by offering deeper insights into the structural complexities of turbid tissue-like media and overcoming existing limitations61–63. In MC modeling of OAM light, initial phase ψ0j is determined according to (2). Correspondingly, initial direction is determined utilizing Poynting vector direction (p→={pρ,pϕ,pz})23:3 pr=ωkrzzR2+z2|LGpℓ(r,ϕ,z)|2,pϕ=ωlr|LGpℓ(r,ϕ,z)|2,pz=ωk|LGpℓ(r,ϕ,z)|2,

where ω is the angular frequency of the light. Regarding the MC method, within the context of corresponding Cartesian coordinates, this is represented by a direction vector (s→={sx,sy,sz}):sx=kzxz2+zR2-ℓyx2+y2,sy=kzyz2+zR2+ℓxx2+y2,sz=k,

As previously above, the polarization tracing framework, rooted in the iterative solution to the Bethe–Salpeter equation and incorporating both Jones and Stokes–Mueller formalisms50, effectively models optical phenomena including reflection and refraction for linear, elliptical, and circular polarizations at the medium’s surface. Within the Bethe–Salpeter-based MC model44,45, a large amount (Ninc>109) of MC-photons with pre-defined statistical weight Wj,j=[1…Ninc] is launched from the source oriented under θi angle to the medium surface, propagates through turbid medium and statistics is collected from those Nph≤Ninc arrived on the detector. Turbid medium is defined by scattering coefficient μs, absorption coefficient μa, anisotropy of scattering g and refractive index n64. Each MC-photon defined at the OAM light source is characterized by the initial statistical weight W0j, Cartesian coordinates (x0j,y0j,0), propagation direction s→0j, initial polarization state and, most importantly, by the initial phase ψ0j. The topological charge of the LG beam determines both s→0j and ψ0j. To track the polarization state along the trajectory of the MC-photon, we introduce a real-valued vector P→, representing the direction of the linearly polarized electric field E→45,46,50.

After launch, all MC-photons undergo surface (z=0) interaction and are transmitted to the turbid medium layer with account for Snell’s law and appropriate Fresnel coefficients influencing MC-photon weights, directions and polarization. In turbid medium (z>0) each MC-photon trajectory is modeled as a sequence of the elementary simulations containing limited amount of scattering events Nscatt. This procedure has been thoroughly covered in previous works50,64. At each i’th scattering event, i=[1…Nscatt], the following computational steps are performed: random path length li=-lnξ/μs is computed (in this paper, we assume that μa≪μs and ξ∈(0,1] is a uniformly distributed random number), MC-photon is moved to the next position r→i=r→i-1+s→ili with weight attenuated according to the Beer–Lambert law (Wi=Wi-1e-μali), and next propagation direction s→i+1 is evaluated via inversion of the Henyey–Greenstein (HG) phase function65pHG(cosθ′)=14π1-g21+g2-2gcosθ′3/2,

where θ′ is the polar scattering angle in the MC-photon reference plane. Here, we have used position vector r→i=(xi,yi,zi) and unit direction for each scattering event s→i=[sX,sY,sZ]i=[sinθcosφ,sinθsinφ,cosθ]i, with θ,φ as azimuthal and polar angles that correspond to the global Cartesian coordinates. It should be noted that, basically, any phase function p can be used47,66. If analytical inversion of p is not possible, then table lookup method is involved to ensure fast computational speed. At each step we check if MC-photon path crosses the medium boundary and invoke surface refraction-transmission and detection procedures if this is the case. Evolution of each linearly polarized state P→x=Pxx,Pxy,Pxz, P→y=Pyx,Pyy,Pyz can be traced along MC-photon trajectory r→i, i=[1…Nscatt] via the following procedure which is obtained from the iterative solution to Bethe–Salpeter equation50,63:4 P→i=-s→i×s→i×P→i-1=I^-s→i⊗s→iP→i-1,

where I^ is the third-rank unit tensor and ⊗ indicates a direct product.

We repeat outlined computational steps for each scattering event until one of the following conditions is met: either Wi<10-4 (statistical weight becomes negligible as follows from Beer–Lambert law) or the amount of scattering events Nscatt becomes larger than 103. These limitations ensure proper trajectory tracing cut-off63. We continue launching MC-photons until the certain amount (no less than Nph=107) arrives on the detector. Detection procedure consists of the two checks: MC-photon coordinates should lie within the detector area (-rd≤xN≤rd,-rd≤yN≤rd,zN=0), and refracted direction s→N should meet the detector numerical aperture (NA) requirements. We would limit those directions by using acos(sN·sd)<NA, where sd=[sin(-θd),0,cos(-θd)] is a unit vector collinear to the detector axis. Both here and in the subsequent sections N is considered to be an index of the detection event.

Thus, each detected MC-photon has defined by the following total statistical weight influenced by its path within medium and depolarization:WNj=W0jPxx2+Pyx2+Pxy2+Pyy2NjΓRjexp-μa∑i=1Njli,

where 0<Nj<Nscatt is the index of detection event for j’th MC-photon, li is the path length between two neighbouring scattering events and ΓR is the Rayleigh factor45,49,63. Based on our Bathe-Salpeter based polarization tracking approach, we are also able to introduce polarized Wj‖ and depolarized Wj⊥ weights for each MC-photon:5 Wj‖=W0jPxx2+Pyx2NjΓRjexp-μa∑i=1Njli,

6 Wj⊥=W0jPxy2+Pyy2NjΓRjexp-μa∑i=1Njli.

By ensemble averaging within each detector pixel we are able to obtain both intensity and phase values on the detector:Ipx(x,y)=∑j=1NpxWNj(x,y)+2∑i=1Npx∑j>iNpxWNi(x,y)WNj(x,y)cos(ΨNi(x,y)-ΨNj(x,y)),ψ(x,y)=arctan∑j=1NpxWNj(x,y)sinΨNj(x,y)∑j=1NpxWNj(x,y)cosΨNj(x,y).

Here, we assume that Npx<Nph is the amount of photons that arrived at the specific pixel, ΨNj is the phase of the j-th detected photon, Ipx is the resulting intensity in this pixel and ψ(x,y) is the resulting phase value in this pixel.

The MC method, a foundational technique in photon transport simulation, has historically encountered the challenge of balancing detailed, accurate modeling with significant computational demands. It has been successfully used for the imitation of image transfer through the complex scattering media, taking into account major experimental parameter, including density, propagation distance, medium optical properties, beam waist, coherence, polarization, etc36. The MC is a robust and versatile tool for modeling light propagation through scattering media, valid across various scattering regimes, including Rayleigh (X/λ≪1), Mie (X/λ≈1), and geometrical optics (X/λ≫1) regimes; Here X is the smallest dimension of the main scatterer. Renowned for its exceptional precision in representing the intricate interactions between light and biological tissues, the MC method is widely regarded as the ‘gold standard’39.

Despite this, the MC method’s reliance on extensive computational resources has been a longstanding issue, with traditional simulations relying heavily on central processing units (CPUs). While the shift from CPU to GPU processing marked a revolutionary decrease in computation time—transforming processes that once spanned hours into mere seconds—this evolution brought with it an increased burden on power consumption and a heightened environmental impact, reflecting a growing concern in our climate-conscious era.

The above mentioned realization of MC implements, for the first time, the energy-effective MC algorithm for complex light propagation utilizing the recently introduced Apple M-family processors67. These cutting-edge chips feature a unified memory architecture, which became a revolutionary advancement in our needs of simulating complex light transfer. This development effectively eliminates the previous constraints on the number of photon packets or statistics that need to be stored and processed, dramatically broadening the scope for MC simulations of photon behaviour. Such capability enables the detailed study/statistical analyses of light propagation through turbid media without the computational and environmental costs historically associated with such efforts. In our evaluations, that the efficiency of these processors has been highlighted by their ability to conduct simulations using approximately 100 times less energy than their traditional GPU counterparts and achieve 300 to 10000 times speed increases than those possible with conventional MC simulations.

The M-family chips possess a sophisticated architecture uniquely optimized for parallel processing, a key factor in their superior performance. We utilize the Metal Compute framework developed by Apple, which plays a crucial role in harnessing these capabilities for MC simulations. Metal provides a low-overhead, hardware-accelerated graphics and compute Application Programming Interface (API) finely tuned for parallel data processing, offering a C-style programming language that is both powerful and efficient. Within this framework, compute shaders act as programmable kernel functions, specially designed to execute MC algorithms on the GPU, enabling the detailed and complex simulation of photon trajectories with unparalleled efficiency and speed.

Supplementary Information

Supplementary Information.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-70954-x.

Author contributions

Conceptualization: IM; Methodology: IM, AD, IL,AS, NM; Experiment: AS, AB, IL; Data acquisition and computer modeling: FK, NM, IL, AD; Funding acquisition: IM, AB, ER; Project administration: IM, AB, ER; Supervision: IM, AB, ER; Writing—original draft: FK, NM, IL,AS, AD, IM; Writing—review and editing: IM.

Funding

This article is based upon work from COST Action CA21159 - Understanding interaction light—biological surfaces: possibility for new electronic materials and devices (PhoBioS) and , supported by COST (European Cooperation in Science and Technology). Authors also acknowledge the support from the UKKi UK-Israel innovation researcher mobility, Academy of Finland (grant projects 325097, 351068) and the support of the EPSRC project EP/W002868/1.

Data availability

All data related to the experiments and computational modeling described in this article are archived on a lab computer at Aston University. All data are available from the corresponding author upon reasonable request.

Code availability

Our models are open-source and also available online as web applications, providing a variety of other MC simulations such as sampling volume, fluence rates, spectra, and colour. The codes used for the modeling of conversion of OAM for different types of LG beams and evaluation the twist of OAM in experimental studies are available from the corresponding author upon reasonable request.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Cao H Mosk AP Rotter S Shaping the propagation of light in complex media Nat. Phys. 2022 18 9 994 1007 10.1038/s41567-022-01677-x
Cao, H., Mosk, A. P. & Rotter, S. Shaping the propagation of light in complex media. Nat. Phys. 18(9), 994–1007 (2022).10.1038/s41567-022-01677-x
2. Mishchenko MI Travis LD Lacis AA Multiple Scattering of Light by Particles: Radiative Transfer and Coherent Backscattering 2006 Cambridge Cambridge University Press
Mishchenko, M. I., Travis, L. D. & Lacis, A. A. Multiple Scattering of Light by Particles: Radiative Transfer and Coherent Backscattering (Cambridge University Press, Cambridge, 2006).
3. Choy TC Effective Medium Theory 2016 Oxford Oxford University Press
Choy, T. C. Effective Medium Theory (Oxford University Press, Oxford, 2016).
4. Wenshan C Shalaev V Optical Metamaterials: Fundamentals and Application 2009 New York Springer
Wenshan, C. & Shalaev, V. Optical Metamaterials: Fundamentals and Application (Springer, New York, 2009).
5. Chandrasekhar S Radiative Transfer 1950 New York Oxford University Press
Chandrasekhar, S. Radiative Transfer (Oxford University Press, New York, 1950).
6. Ishimaru A Wave Propagation and Scattering in Random Media, vol I and II 1978 New York Academic
Ishimaru, A. Wave Propagation and Scattering in Random Media, vol I and II (Academic, New York, 1978).
7. Lu B Morgan SP Crowe JA Stockford IM Comparison of methods for reducing the effects of scattering in spectrophotometry Appl. Spectrosc. 2006 60 1157 66 10.1366/000370206778664725 17059668
Lu, B., Morgan, S. P., Crowe, J. A. & Stockford, I. M. Comparison of methods for reducing the effects of scattering in spectrophotometry. Appl. Spectrosc. 60, 1157–66 (2006).17059668 10.1366/000370206778664725
8. Lee H Yoon S Loohuis P Hong JH Kang S Choi W High-throughput volumetric adaptive optical imaging using compressed time-reversal matrix Light Sci. Appl. 2022 11 16 10.1038/s41377-021-00705-4 35027538
Lee, H. et al. High-throughput volumetric adaptive optical imaging using compressed time-reversal matrix. Light Sci. Appl. 11, 16 (2022).35027538 10.1038/s41377-021-00705-4
9. Sanjeev A Kapellner Y Shabairou N Gur E Sinvani M Zalevsky Z Non-invasive imaging through scattering medium by using a reverse response wavefront shaping technique Sci. Rep. 2020 10 6029 10.1038/s41598-020-62442-9 32238830
Sanjeev, A. et al. Non-invasive imaging through scattering medium by using a reverse response wavefront shaping technique. Sci. Rep. 10, 6029 (2020).32238830 10.1038/s41598-020-62442-9
10. Sanjeev A Trivedi V Zalevsky Z Optical reciprocity induced wavefront shaping for axial and lateral shifting of focus through a scattering medium Sci. Rep. 2022 11 6387 10.1038/s41598-022-10378-7
Sanjeev, A., Trivedi, V. & Zalevsky, Z. Optical reciprocity induced wavefront shaping for axial and lateral shifting of focus through a scattering medium. Sci. Rep. 11, 6387 (2022).10.1038/s41598-022-10378-7
11. Gigan S Katz O Aguiar HB Andresen ER Aubry A Bertolotti J Bossy E Bouchet D Brake J Brasselet S Bromberg Y Cao H Chaigne T Cheng Z Choi W Čižmár T Cui M Curtis VR Defienne H Hofer M Horisaki R Horstmeyer R Ji N LaViolette AK Mertz J Moser C Mosk AP Pégard NC Piestun R Popoff S Phillips DB Psaltis D Rahmani B Rigneault H Rotter S Tian L Vellekoop IM Waller L Wang L Weber T Xiao S Xu C Yamilov A Yang C Yılmaz H Roadmap on wavefront shaping and deep imaging in complex media J. Phys. Photon. 2022 4 4 042501 10.1088/2515-7647/ac76f9
...Gigan, S. et al. Roadmap on wavefront shaping and deep imaging in complex media. J. Phys. Photon. 4(4), 042501 (2022).10.1088/2515-7647/ac76f9
12. Yu Z Li H Zhong T Park JH Cheng S Woo CM Zhao Q Yao J Zhou Y Huang X Pang W Yoon H Shen Y Liu H Zheng Y Park Y Wang LV Lai P Wavefront shaping: A versatile tool to conquer multiple scattering in multidisciplinary fields Innovation (Cambridge) 2022 3 100292
Yu, Z. et al. Wavefront shaping: A versatile tool to conquer multiple scattering in multidisciplinary fields. Innovation (Cambridge) 3, 100292 (2022).
13. Shen Y Wang X Xie Z Min C Fu X Liu Q Gong M Yuan X Optical vortices 30 years on: OAM manipulation from topological charge to multiple singularities Light Sci. Appl. 2019 8 90 10.1038/s41377-019-0194-2 31645934
Shen, Y. et al. Optical vortices 30 years on: OAM manipulation from topological charge to multiple singularities. Light Sci. Appl. 8, 90 (2019).31645934 10.1038/s41377-019-0194-2
14. He C Shen Y Forbes A Towards higher-dimensional structured light Light Sci. Appl. 2022 11 205 10.1038/s41377-022-00897-3 35790711
He, C., Shen, Y. & Forbes, A. Towards higher-dimensional structured light. Light Sci. Appl. 11, 205 (2022).35790711 10.1038/s41377-022-00897-3
15. Cao H Čižmár T Turtaev S Tyc T Rotter S Controlling light propagation in multimode fibers for imaging, spectroscopy, and beyond Adv. Opt. Photon. 2023 15 2 524 612 10.1364/AOP.484298
Cao, H., Čižmár, T., Turtaev, S., Tyc, T. & Rotter, S. Controlling light propagation in multimode fibers for imaging, spectroscopy, and beyond. Adv. Opt. Photon. 15(2), 524–612 (2023).10.1364/AOP.484298
16. Wang WB Gozali R Shi L Lindwasser L Alfano RR Deep transmission of Laguerre–Gaussian vortex beams through turbid scattering media Opt. Lett. 2016 41 9 2069 2072 10.1364/OL.41.002069 27128076
Wang, W. B., Gozali, R., Shi, L., Lindwasser, L. & Alfano, R. R. Deep transmission of Laguerre–Gaussian vortex beams through turbid scattering media. Opt. Lett. 41(9), 2069–2072 (2016).27128076 10.1364/OL.41.002069
17. Angelsky OV Mokhun II Bekshaev AY Zenkova CY Zheng J Polarization singularities: Topological and dynamical aspects Front. Phys. 2023 11 1147788 10.3389/fphy.2023.1147788
Angelsky, O. V., Mokhun, I. I., Bekshaev, A. Y., Zenkova, C. Y. & Zheng, J. Polarization singularities: Topological and dynamical aspects. Front. Phys. 11, 1147788 (2023).10.3389/fphy.2023.1147788
18. Forbes A Advances in orbital angular momentum lasers J. Light. Technol. 2023 41 7 2079 2086 10.1109/JLT.2022.3220509
Forbes, A. Advances in orbital angular momentum lasers. J. Light. Technol. 41(7), 2079–2086 (2023).10.1109/JLT.2022.3220509
19. Mamani S Shintre S Li Z Rodríguez-Contreras A Shi L Alfano R OAM transmission of polarized multipole laser beams in rat cerebellum tissue Opt. Commun. 2023 532 129241 10.1016/j.optcom.2022.129241
Mamani, S. et al. OAM transmission of polarized multipole laser beams in rat cerebellum tissue. Opt. Commun. 532, 129241 (2023).10.1016/j.optcom.2022.129241
20. Biton N Kupferman J Arnon S OAM light propagation through tissue Sci. Rep. 2021 11 1 2407 10.1038/s41598-021-82033-6 33510283
Biton, N., Kupferman, J. & Arnon, S. OAM light propagation through tissue. Sci. Rep. 11(1), 2407 (2021).33510283 10.1038/s41598-021-82033-6
21. Wang, W. B. et al. Optical vortex beam transmission with different OAM in scattering beads and brain tissue media. In Complex Light and Optical Forces X, Vol. 9764 114–119 (2016).
22. Wang Y Bai L Huang C Xie J Zhang D Guo L Orbital angular momentum of Laguerre–Gaussian beams with non-zero radial index at limited aperture size Results Phys. 2023 48 106436 10.1016/j.rinp.2023.106436
Wang, Y. et al. Orbital angular momentum of Laguerre–Gaussian beams with non-zero radial index at limited aperture size. Results Phys. 48, 106436 (2023).10.1016/j.rinp.2023.106436
23. Allen L Padgett MJ Babiker M The Orbital Angular Momentum of Light 1999 North-Holland Elsevier 291 372
Allen, L., Padgett, M. J. & Babiker, M. The Orbital Angular Momentum of Light Vol. 39, 291–372 (Elsevier, North-Holland, 1999).
24. Rosales-Guzmán C Ndagano B Forbes A A review of complex vector light fields and their applications J. Opt. 2018 20 12 123001 10.1088/2040-8986/aaeb7d
Rosales-Guzmán, C., Ndagano, B. & Forbes, A. A review of complex vector light fields and their applications. J. Opt. 20(12), 123001 (2018).10.1088/2040-8986/aaeb7d
25. Berry MV McDonald KT Exact and geometrical optics energy trajectories in twisted beams J. Opt. A Pure Appl. Opt. 2008 10 3 035005 10.1088/1464-4258/10/3/035005
Berry, M. V. & McDonald, K. T. Exact and geometrical optics energy trajectories in twisted beams. J. Opt. A Pure Appl. Opt. 10(3), 035005 (2008).10.1088/1464-4258/10/3/035005
26. Yoo KM Liu F Alfano RR When does the diffusion approximation fail to describe photon transport in random media? Phys. Rev. Lett. 1990 64 2647 2650 10.1103/PhysRevLett.64.2647 10041774
Yoo, K. M., Liu, F. & Alfano, R. R. When does the diffusion approximation fail to describe photon transport in random media?. Phys. Rev. Lett. 64, 2647–2650 (1990).10041774 10.1103/PhysRevLett.64.2647
27. Berkhout GCG Lavery MPJ Courtial J Beijersbergen MW Padgett MJ Efficient sorting of orbital angular momentum states of light Phys. Rev. Lett. 2010 105 153601 10.1103/PhysRevLett.105.153601 21230900
Berkhout, G. C. G., Lavery, M. P. J., Courtial, J., Beijersbergen, M. W. & Padgett, M. J. Efficient sorting of orbital angular momentum states of light. Phys. Rev. Lett. 105, 153601 (2010).21230900 10.1103/PhysRevLett.105.153601
28. Meglinski, I., Sdobnov, A., Lopushenko, I. & Bykov, A. Phase memory of orbital angular momentum in multiple scattering environment. Laser Sci. Appl. (arXiv preprint arXiv:2312.08928) (2024–in press)
29. Mach L Ein neuer Interferenzrefraktor Z. Instrum. 1892 12 89 93
Mach, L. Ein neuer Interferenzrefraktor. Z. Instrum. 12, 89–93 (1892).
30. Zehnder L Ein neuer Interferenzrefraktor Z. Instrum. 1891 11 275 285
Zehnder, L. Ein neuer Interferenzrefraktor. Z. Instrum. 11, 275–285 (1891).
31. Padgett M Arlt J Simpson N Allen L An experiment to observe the intensity and phase structure of Laguerre–Gaussian laser modes Am. J. Phys. 1996 64 1 77 82 10.1119/1.18283
Padgett, M., Arlt, J., Simpson, N. & Allen, L. An experiment to observe the intensity and phase structure of Laguerre–Gaussian laser modes. Am. J. Phys. 64(1), 77–82 (1996).10.1119/1.18283
32. Kumar P Nishchal NK Modified Mach–Zehnder interferometer for determining the high-order topological charge of Laguerre-Gaussian vortex beams J. Opt. Soc. Am. A 2019 36 8 1447 1455 10.1364/JOSAA.36.001447
Kumar, P. & Nishchal, N. K. Modified Mach–Zehnder interferometer for determining the high-order topological charge of Laguerre-Gaussian vortex beams. J. Opt. Soc. Am. A 36(8), 1447–1455 (2019).10.1364/JOSAA.36.001447
33. Sieryi O Popov A Kalchenko V Bykov A Meglinski I Tissue-mimicking phantoms for biomedical applications Proc. SPIE 2020 11363 1136312
Sieryi, O., Popov, A., Kalchenko, V., Bykov, A. & Meglinski, I. Tissue-mimicking phantoms for biomedical applications. Proc. SPIE 11363, 1136312 (2020).
34. Cui S Xu B Luo S Xu H Cai Z Luo Z Pu J Chávez-Cerda S Determining topological charge based on an improved fizeau interferometer Opt. Expr. 2019 27 9 12774 12779 10.1364/OE.27.012774
Cui, S. et al. Determining topological charge based on an improved fizeau interferometer. Opt. Expr. 27(9), 12774–12779 (2019).10.1364/OE.27.012774
35. Vayalamkuzhi P Bhattacharya S Transform-based phase retrieval techniques from a single off-axis interferogram Appl. Opt. 2021 60 19 5523 5533 10.1364/AO.422900 34263840
Vayalamkuzhi, P. et al. Transform-based phase retrieval techniques from a single off-axis interferogram. Appl. Opt. 60(19), 5523–5533 (2021).34263840 10.1364/AO.422900
36. Berrocal E Meglinski I Greenhalgh DA Linne MA Image transfer through the complex scattering turbid media Laser Phys. Lett. 2006 3 9 464 468 10.1002/lapl.200610035
Berrocal, E., Meglinski, I., Greenhalgh, D. A. & Linne, M. A. Image transfer through the complex scattering turbid media. Laser Phys. Lett. 3(9), 464–468 (2006).10.1002/lapl.200610035
37. Carles G Zammit P Harvey AR Holistic Monte–Carlo optical modelling of biological imaging Sci. Rep. 2019 9 15832 10.1038/s41598-019-51850-1 31676825
Carles, G., Zammit, P. & Harvey, A. R. Holistic Monte–Carlo optical modelling of biological imaging. Sci. Rep. 9, 15832 (2019).31676825 10.1038/s41598-019-51850-1
38. Zhu R Avsievich T Popov A Meglinski I Optical tweezers in the studies of red blood cells Cells 2020 9 3 545 10.3390/cells9030545 32111018
Zhu, R., Avsievich, T., Popov, A. & Meglinski, I. Optical tweezers in the studies of red blood cells. Cells 9(3), 545 (2020).32111018 10.3390/cells9030545
39. Periyasamy V Pramanik M Advances in Monte Carlo simulation for light propagation in tissue IEEE Rev. Biomed. Eng. 2017 10 122 135 10.1109/RBME.2017.2739801 28816674
Periyasamy, V. & Pramanik, M. Advances in Monte Carlo simulation for light propagation in tissue. IEEE Rev. Biomed. Eng. 10, 122–135 (2017).28816674 10.1109/RBME.2017.2739801
40. Xu M Electric field Monte Carlo simulation of polarized light propagation in turbid media Opt. Express 2004 12 6530 6539 10.1364/OPEX.12.006530 19488304
Xu, M. Electric field Monte Carlo simulation of polarized light propagation in turbid media. Opt. Express 12, 6530–6539 (2004).19488304 10.1364/OPEX.12.006530
41. Gangnus SV Matcher SJ Meglinski I Monte Carlo modeling of polarized light propagation in biological tissues Laser Phys. 2004 14 6 886 891
Gangnus, S. V., Matcher, S. J. & Meglinski, I. Monte Carlo modeling of polarized light propagation in biological tissues. Laser Phys. 14(6), 886–891 (2004).
42. Ramella-Roman JC Prahl SA Jacques SL Three Monte Carlo programs of polarized light transport into scattering media: Part I Opt. Express 2005 13 12 4420 4438 10.1364/OPEX.13.004420 19495358
Ramella-Roman, J. C., Prahl, S. A. & Jacques, S. L. Three Monte Carlo programs of polarized light transport into scattering media: Part I. Opt. Express 13(12), 4420–4438 (2005).19495358 10.1364/OPEX.13.004420
43. Ramella-Roman JC Prahl SA Jacques SL Three Monte Carlo programs of polarized light transport into scattering media: part II Opt. Express 2005 13 25 10392 10405 10.1364/OPEX.13.010392 19503254
Ramella-Roman, J. C., Prahl, S. A. & Jacques, S. L. Three Monte Carlo programs of polarized light transport into scattering media: part II. Opt. Express 13(25), 10392–10405 (2005).19503254 10.1364/OPEX.13.010392
44. Kuzmin VL Meglinski IV Coherent multiple scattering effects and Monte Carlo method JETP Lett. 2004 79 109 112 10.1134/1.1719124
Kuzmin, V. L. & Meglinski, I. V. Coherent multiple scattering effects and Monte Carlo method. JETP Lett. 79, 109–112 (2004).10.1134/1.1719124
45. Meglinski I Kuzmin VL Churmakov DY Greenhalgh DA Monte Carlo simulation of coherent effects in multiple scattering Proc. R. Soc. A 2005 463 43 53 10.1098/rspa.2004.1369
Meglinski, I., Kuzmin, V. L., Churmakov, D. Y. & Greenhalgh, D. A. Monte Carlo simulation of coherent effects in multiple scattering. Proc. R. Soc. A 463, 43–53 (2005).10.1098/rspa.2004.1369
46. Kuzmin VL Meglinski IV Coherent effects of multiple scattering for scalar and electromagnetic fields: Monte–Carlo simulation and Milne-like solutions Opt. Commun. 2007 273 2 307 310 10.1016/j.optcom.2007.01.025
Kuzmin, V. L. & Meglinski, I. V. Coherent effects of multiple scattering for scalar and electromagnetic fields: Monte–Carlo simulation and Milne-like solutions. Opt. Commun. 273(2), 307–310 (2007).10.1016/j.optcom.2007.01.025
47. Doronin A Radosevich AJ Backman V Meglinski I Two electric field Monte Carlo models of coherent backscattering of polarized light J. Opt. Soc. Am. A 2014 31 11 2394 2400 10.1364/JOSAA.31.002394
Doronin, A., Radosevich, A. J., Backman, V. & Meglinski, I. Two electric field Monte Carlo models of coherent backscattering of polarized light. J. Opt. Soc. Am. A 31(11), 2394–2400 (2014).10.1364/JOSAA.31.002394
48. Meglinski I Kuz’min VL Coherent backscattering of circularly polarized optical radiation from a disperse random medium Prog. Electromagn. Res. M 2011 21 1972 1977
Meglinski, I. & Kuz’min, V. L. Coherent backscattering of circularly polarized optical radiation from a disperse random medium. Prog. Electromagn. Res. M 21, 1972–1977 (2011).
49. Akkermans E Wolf PE Maynard R Maret G Theoretical study of the coherent backscattering of light by disordered media J. Phys. France 1988 49 77 98 10.1051/jphys:0198800490107700
Akkermans, E., Wolf, P. E., Maynard, R. & Maret, G. Theoretical study of the coherent backscattering of light by disordered media. J. Phys. France 49, 77–98 (1988).10.1051/jphys:0198800490107700
50. Lopushenko I Sieryi O Bykov A Meglinski I Exploring the evolution of circular polarized light backscattered from turbid tissue-like disperse medium utilizing generalized Monte Carlo modeling approach with a combined use of Jones and Stokes–Mueller formalisms J. Biomed. Opt. 2024 29 052913 38089555
Lopushenko, I., Sieryi, O., Bykov, A. & Meglinski, I. Exploring the evolution of circular polarized light backscattered from turbid tissue-like disperse medium utilizing generalized Monte Carlo modeling approach with a combined use of Jones and Stokes–Mueller formalisms. J. Biomed. Opt. 29, 052913 (2024).38089555
51. Mishchenko MI Vector radiative transfer equation for arbitrarily shaped and arbitrarily oriented particles: A microphysical derivation from statistical electromagnetics Appl. Opt. 2002 41 7114 7134 10.1364/AO.41.007114 12463259
Mishchenko, M. I. Vector radiative transfer equation for arbitrarily shaped and arbitrarily oriented particles: A microphysical derivation from statistical electromagnetics. Appl. Opt. 41, 7114–7134 (2002).12463259 10.1364/AO.41.007114
52. Raković MJ Kattawar GW Mehrűbeoğlu M Cameron BD Wang LV Rastegar S Coté GL Light backscattering polarization patterns from turbid media: Theory and experiment Appl. Opt. 1999 38 3399 3408 10.1364/AO.38.003399 18319938
Raković, M. J. et al. Light backscattering polarization patterns from turbid media: Theory and experiment. Appl. Opt. 38, 3399–3408 (1999).18319938 10.1364/AO.38.003399
53. Tynes HH Kattawar GW Zege EP Katsev IL Prikhach AS Chaikovskaya LI Monte Carlo and multicomponent approximation methods for vector radiative transfer by use of effective Mueller matrix calculations Appl. Opt. 2001 40 400 412 10.1364/AO.40.000400 18357013
Tynes, H. H. et al. Monte Carlo and multicomponent approximation methods for vector radiative transfer by use of effective Mueller matrix calculations. Appl. Opt. 40, 400–412 (2001).18357013 10.1364/AO.40.000400
54. Doicu A Mishchenko MI An overview of methods for deriving the radiative transfer theory from the Maxwell equations. II: Approach based on the Dyson and Bethe–Salpeter equations J. Quant. Spectrosc. Radiat. Transf. 2019 224 25 36 10.1016/j.jqsrt.2018.10.032 30713354
Doicu, A. & Mishchenko, M. I. An overview of methods for deriving the radiative transfer theory from the Maxwell equations. II: Approach based on the Dyson and Bethe–Salpeter equations. J. Quant. Spectrosc. Radiat. Transf. 224, 25–36 (2019).30713354 10.1016/j.jqsrt.2018.10.032
55. Günhan Akarçay, H., Hohmann, A., Kienle, A., Frenz, M., & Rička, J.: Monte Carlo modeling of polarized light propagation: Stokes vs. Jones. Part I. Appl. Opt. 53(31), 7576–7585 (2014).
56. Ivanov D Dremin V Bykov A Borisova E Genova T Popov A Ossikovski R Novikova T Meglinski I Colon cancer detection via Poincaré sphere representation and 2D polarimetric mapping of ex vivo tissue samples J. Biophoton. 2020 13 202000082 10.1002/jbio.202000082
Ivanov, D. et al. Colon cancer detection via Poincaré sphere representation and 2D polarimetric mapping of ex vivo tissue samples. J. Biophoton. 13, 202000082 (2020).10.1002/jbio.202000082
57. Borovkova MA Bykov AV Popov A Meglinski IV Role of scattering and birefringence in phase retardation revealed by locus of Stokes vector on Poincaré sphere J. Biomed. Opt. 2020 25 5 057001 10.1117/1.JBO.25.5.057001 32436372
Borovkova, M. A., Bykov, A. V., Popov, A. & Meglinski, I. V. Role of scattering and birefringence in phase retardation revealed by locus of Stokes vector on Poincaré sphere. J. Biomed. Opt. 25(5), 057001 (2020).32436372 10.1117/1.JBO.25.5.057001
58. Singh MD Vitkin IA Spatial helicity response metric to quantify particle size and turbidity of heterogeneous media through circular polarization imaging Sci. Rep. 2023 13 1 2231 10.1038/s41598-023-29444-9 36755076
Singh, M. D. & Vitkin, I. A. Spatial helicity response metric to quantify particle size and turbidity of heterogeneous media through circular polarization imaging. Sci. Rep. 13(1), 2231 (2023).36755076 10.1038/s41598-023-29444-9
59. Nishizawa N Kuchimaru T Depth estimation of tumor invasion in early gastric cancer using scattering of circularly polarized light: Monte carlo simulation study J. Biophoton. 2022 15 10 202200062 10.1002/jbio.202200062
Nishizawa, N. & Kuchimaru, T. Depth estimation of tumor invasion in early gastric cancer using scattering of circularly polarized light: Monte carlo simulation study. J. Biophoton. 15(10), 202200062 (2022).10.1002/jbio.202200062
60. Lopushenko I Bykov A Meglinski I Depolarization composition of backscattered circularly polarized light Phys. Rev. A 2023 108 041502 10.1103/PhysRevA.108.L041502
Lopushenko, I., Bykov, A. & Meglinski, I. Depolarization composition of backscattered circularly polarized light. Phys. Rev. A 108, 041502 (2023).10.1103/PhysRevA.108.L041502
61. Doronin A Milione G Meglinski I Alfano RR Propagation and scattering of vector light beam in turbid scattering medium Proc. SPIE 2014 8940 894006 10.1117/12.2038818
Doronin, A., Milione, G., Meglinski, I. & Alfano, R. R. Propagation and scattering of vector light beam in turbid scattering medium. Proc. SPIE 8940, 894006 (2014).10.1117/12.2038818
62. Doronin, A., Vera, N., Staforelli, J.P., Coelho, P. & Meglinski, I. Propagation of cylindrical vector laser beams in turbid tissue-like scattering media. Photonics 6(2), 56–67 (2019)
63. Doronin, A., Novikova, V. N. Tatiana, Staforelli, J. P. & Meglinski, I. Assessment of twisted light localization in turbid tissue-like scattering media using 3D geometrical exploration. In Proceedings of SPIE PC12373 1237308 (2023).
64. Meglinski I Doronin A Wang RK Tuchin VV Chapter 1 Advanced Biophotonics: Tissue Optical Sectioning 2013 Boca Raton CRC Press 1 72
Meglinski, I. & Doronin, A. Chapter 1. In Advanced Biophotonics: Tissue Optical Sectioning (eds Wang, R. K. & Tuchin, V. V.) 1–72 (CRC Press, Boca Raton, 2013).
65. Henyey LG Greenstein JL Diffuse radiation in the galaxy Astrophys. J. 1941 93 70 83 10.1086/144246
Henyey, L. G. & Greenstein, J. L. Diffuse radiation in the galaxy. Astrophys. J. 93, 70–83 (1941).10.1086/144246
66. Kuz’min VL Val’kov AY Zubkov LA Photon diffusion in random media and anisotropy of scattering in the Henyey–Greenstein and Rayleigh–Gans models J. Exp. Theor. Phys. 2019 128 396 406 10.1134/S1063776119020109
Kuz’min, V. L., Val’kov, A. Y. & Zubkov, L. A. Photon diffusion in random media and anisotropy of scattering in the Henyey–Greenstein and Rayleigh–Gans models. J. Exp. Theor. Phys. 128, 396–406 (2019).10.1134/S1063776119020109
67. Clennell A Nguyen V Yakovlev VS Doronin A Neu(t)ralMC: Energy-efficient open source Monte Carlo algorithm for assessing photon transport in turbid media Opt. Express 2023 31 19 30921 30931 10.1364/OE.496516 37710624
Clennell, A., Nguyen, V., Yakovlev, V. S. & Doronin, A. Neu(t)ralMC: Energy-efficient open source Monte Carlo algorithm for assessing photon transport in turbid media. Opt. Express 31(19), 30921–30931 (2023).37710624 10.1364/OE.496516
