Simulation Of Light Propagation Research Articles

Field propagation study in ultra-small core photonic crystal fiber beyond the fundamental mode cut-off

We present the study of the behavior of the electromagnetic field propagating in the ultra-small core photonic crystal fiber (PCF) beyond its fundamental mode cut-off. We first use plane waves expansion (PWE) method to find necessary parameters and then involve finite-difference time-domain (FDTD) method accelerated with graphic processing unit (GPU) for long-range light propagation simulation. The confinement losses beyond $\lambda = 1200 nm$ achieves tens of dB/cm which can be lowered significantly by increasing the number of cladding rings.

Complete head cerebral sensitivity mapping for diffuse correlation spectroscopy using subject-specific magnetic resonance imaging models.

We characterize cerebral sensitivity across the entire adult human head for diffuse correlation spectroscopy, an optical technique increasingly used for bedside cerebral perfusion monitoring. Sixteen subject-specific magnetic resonance imaging-derived head models were used to identify high sensitivity regions by running Monte Carlo light propagation simulations at over eight hundred uniformly distributed locations on the head. Significant spatial variations in cerebral sensitivity, consistent across subjects, were found. We also identified correlates of such differences suitable for real-time assessment. These variations can be largely attributed to changes in extracerebral thickness and should be taken into account to optimize probe placement in experimental settings.

Rationale for the Choice of the Ellipsoidal Reflector Parameters for Biomedical Photometers

Biomedical photometersʼ information-measuring systems with ellipsoidal reflectors have acceptable results in determining of biological tissues optical properties in the visible and near-infrared spectral range. These photometers make it possible to study the optical radiation propagation in turbid media for direct and inverse problems of light-scattering optics. The purpose of this work is to study the influence of the ellipsoidal reflectors design parameters on the results of biomedical photometry when simulating the optical radiation propagation in a system of biological tissue and reflectors in transmitted and reflected light.The paper substantiates the choice of the ellipsoidal reflectors’ focal parameter for efficient registration of forward and backscattered light. The methodology of the process is illustrated by the results of a model experiment using the Monte Carlo simulation for samples of human brain white and gray matter at the visible range of 405 nm, 532 nm, and 650 nm. The total transmittance, diffuse reflectance, and absorption graphs depending on the sample thickness were obtained. Based on the introduced concepts of the ellipsoidal reflector efficiency index and its efficiency factor, the expediency of choosing the ellipsoidal reflectors focal parameter is analyzed to ensure the registration of the maximum amount of scattered light. The graphs of efficiency index in reflected and transmitted light for different thickness samples of white and gray matter and efficiency factors depending on the sample thickness were obtained.The influence of the reflectors ellipticity on the illuminance of various zones of photometric images using the example of an absorbing biological medium – pig liver tissue – at wavelength of 405 nm with a Monte Carlo simulation was analyzed.The optical properties of biological media (scattering and absorption coefficients, scattering anisotropy factor, refractive index) and the samples’ geometric dimensions, particularly the thickness, are predetermined when choosing the ellipsoidal reflectors parameters for registration of the scattered light. Coordinates of the output of photons and their statistical weight obtained in the Monte Carlo simulation of light propagation in biological tissue have a physical effect on a characteristic scattering spot formation in the receiving plane of a biomedical photometer with ellipsoidal reflectors.

Tissue biomolecular and microstructure profiles in optical colorectal cancer delineation

Colorectal cancer (CRC) is the 3rd most common and the 2nd most deadly type of cancer worldwide. Understanding the biochemical and microstructural aspects of carcinogenesis is a critical step towards developing new technologies for accurate CRC detection. To date, optical detection through analyzing tissue chromophore concentrations and scattering parameters has been mostly limited to chromophores in the visible region and analytical light diffusion models. In this study, tissue parameters were extracted by fitting diffuse reflectance spectra (DRS) within the range 350–1900 nm based on reflectance values from a look-up table built using Monte Carlo simulations of light propagation in tissues. This analysis was combined with machine learning models to estimate parameter thresholds leading to best differentiation between mucosa and tumor tissues based on almost 3000 DRS recorded from fresh ex vivo tissue samples from 47 subjects. DRS spectra were measured with a probe for superficial tissue and another for slightly deeper tissue layers. By using the classification and regression tree algorithm, the most important parameters for CRC detection were the total lipid content (f lipid), the reduced scattering amplitude (α′), and the Mie scattering power (b Mie). Successful classification with an area under the receiver operating characteristic curve higher than 90% was achieved. To the best of our knowledge, this is the first study to evaluate the potential tissue biomolecule concentrations and scattering properties in superficial and deeper tissue layers for CRC detection in the luminal wall. This may have important clinical applications for the rapid diagnosis of colorectal neoplasia.

GiBUU based neutrino interaction simulations in KM3NeT

The simulation of the neutrino interaction is a crucial step in the simulation chainof a neutrino experiment. The different processes taking part in the neutrino scatteringon a nucleus require several approximations in order to make the simulation possibleand to realize reasonable computation times.This can be realised in different ways, e.g. by parametrised models for the different scattering processes and energy regimes as it is implemented in GENIE.The GiBUU neutrino generator utilises the Boltzmann-Uehling-Uhlenbeck equation tosimulate the particle flow after the neutrino interactions, the so-called final stateinteractions.The detector-specific results in form of the visible energy in the detector afterthe light propagation simulation and the KM3NeT event reconstruction are presented.In addition to that, the comparison to the GENIE based simulation environment in KM3NeT (gSeaGen) is drawn.

Determination of Optical Properties and Photodynamic Threshold of Lung Tissue for Treatment Planning of In Vivo Lung Perfusion Assisted Photodynamic Therapy

BackgroundIsolated lung metastases in sarcoma and colorectal cancer patients are inadequately treated with current standard therapies. In Vivo Lung Perfusion, a novel platform, could overcome limitations to photodynamic therapy treatment volumes by using low cellular perfusate, removing blood, theoretically allowing greater light penetration. To develop personalized photodynamic therapy protocols requires in silico light propagation simulations based on optical properties and maximal permissible photodynamic threshold dose of lung tissue. This study presents quantification of optical properties for two perfusates and the photodynamic threshold for 5-ALA and Chlorin e6. MethodsPorcine and human lungs were placed on Ex Vivo Lung Perfusion, and perfused with acellular solution or blood. Isotropic diffusers were placed within bronchi and on lung surface for light transmission measurements, from which absorption and light scattering properties were calculated at multiple wavelengths. Separately, pigs were injected with 5-ALA or Chlorin e6, and lung tissue was irradiated at increasing doses. Resultant lesion sizes were measured by CT and histology to quantify the photodynamic threshold. ResultsLow cellular perfusate reduced the tissue absorption coefficient significantly, increasing penetration depth of light by 3.3 mm and treatment volumes 3-fold. The photodynamic threshold for lung exposed to 5-ALA was consistent with other malignancies. Chlorin e6 levels were undetectable in lung tissue and did not demonstrate photodynamic-induced necrosis. ConclusionsLight penetration with low cellular perfusate is significantly greater and could enable treatments for diffuse disease. This data aids photodynamic treatment planning and will guide clinical translation of photodynamic therapy protocols in the lung, especially during lung perfusion.

Comparative Evaluation of Photodynamic Therapeutic Probes for Peripheral Lung Cancer Based on Light Propagation Simulation

本研究では，計算機シミュレーションに基づく非臨床試験により，末梢型肺癌に対する光線力学治療用プローブを評価する．モンテカルロ法に基づく光伝搬シミュレーションにより，臨床利用が検討されている直射プローブ，側射プローブにて光照射した際の組織内光分布を算出し，それぞれの照射体積を比較した．直射プローブと比べ，側射プローブが腫瘍組織に対してより効率的にレーザー光照射が可能であることを定量的に示した．

Analysis of the applicability of the classical probabilistic parameters of the Monte Carlo algorithm for problems of light transport in turbid biological media with continuous absorption and discrete scattering

Simulation of light propagation by the statistical Monte Carlo (MC) method is widely used in many fields, especially in astrophysics, atmospheric optics, ocean optics, and nuclear medicine. In the optics of biological tissues, the MC method is used to simulate the luminous flux, which is formed during various medical therapeutic or diagnostic procedures inside a biological tissue and on its surface. In such calculations, the MC method is commonly considered as a reference one, which ensures an arbitrarily high accuracy with an increase in the number of ‘photons’. However, it can be shown that this is not always the case. In this paper, in the methodological aspect, the idealised one-dimensional problems of the transport theory for a turbid medium with continuous absorption and scattering and a turbid medium with discrete scatterers inside a continuously absorbing medium are considered. Their exact analytical solutions are presented and compared with the results of statistical modelling by the MC method. It is found that the use of classical probabilistic parameters for a medium with continuous absorption and scattering in the MC algorithm leads to a systematic method error in determining the values of radiation fluxes for biological media with discrete scattering, up to 10% for fluxes at the boundary in some cases. The causes of the error are discussed and it is shown how to modify the probabilistic parameters of the MC algorithm to eliminate it.

Accurate Simulation of Light Propagation in Complex Skin Tissues Using an Improved Tetrahedron-Based Monte Carlo Method

Understanding light transportation in skin tissues can help improve clinical efficacy in the laser treatment of dermatosis, such as port-wine stains (PWS). Patient-specific cross-bridge PWS vessels are structurally complicated and considerably influence laser energy deposition due to shading effects. The shading effect of PWS vessels is investigated using a tetrahedron-based Monte Carlo (MC) method with extended boundary condition (TMCE). In TMCE, body-fitted tetrahedra are generated in different tissues, and the precision of photon–surface interaction can be considerably improved via mesh refinement. Such improvement is difficult to achieve with the widely used voxel-based MC method. To fit the real physical boundary, the extended boundary condition is adapted by extending photon propagation to the semi-infinite tissue layers while restricting the statistics of photon propagation in the computational domain. Results indicate that the shading parameters, such as the cross angle, vessel distance, and geometric shadow (GS), of cross-bridge blood vessel pairs determine the peak characteristics of photon deposition in deep vessels by affecting the relative deposition of collimated and diffused light. Collimated light is shaded, attenuated, and partially transformed into diffused light due to the increase in vessel distance and GS of vessel pairs, resulting in difficulty in treating deep and shallow vessels with one laser pulse. The TMCE method can be used for the individualized and precise forecasting of laser energy deposition based on the morphology and embedding characteristics of vascular lesions.

Parallel, multi-purpose Monte Carlo code for simulation of light propagation in segmented tissues

We have developed a massive-parallel, multi-purpose Monte Carlo (MC) code for simulation of light propagation in complex structures modelled with e.g. magnetic resonance images. The code is designed to execute parallel threads on a Linux-based cluster of computers equipped with multiple graphical processing units (GPU) utilizing NVidia CUDA technology. We show steps one can take to implements such code itself. Furthermore, we provide methodology of building a MRI-based head model and populating it with realistic optical properties at excitation and fluorescence emission wavelengths during inflow of fluorescence agent (indocyanine green – ICG). The proposed code provides following original features: (i) Simulation of fluorescence light propagation in media with spatial distribution of multiple different fluorophores characterized by concentration, quantum yield and fluorescence lifetime. (ii) The fluorescence light tracking does not need extra photon tracking and works in parallel with tracking photons at the excitation wavelength, introducing execution time overhead of 0.02% only. (iii) Calculation of high-resolution spatial distributions of sensitivity factors mapping voxels absorption change to parameters measured on a model surface: statistical parameters (moments) of the distribution of time of flight of photons. (iv) Random number generators states are preserved between runs to greatly improve calculation time of the sensitivity factors. (v) Simulation of fluorescence lifetime wide field imaging in diffusively scattering media. The code is designed to handle big models, divided into e.g. 200 million voxels or more, which other methods struggle to handle. Simulations on a human head model with 0.3 mm voxel size require 2 GB of a GPU memory only. This is supported by developed non-standard floating-point data storage format. The proposed code is cross-validated with field-leading MC and finite element methods on the same hardware environment for varying model sizes and temporal resolutions. Tests revealed competitive execution time and high temporal resolution of boundary data for high spatial resolution of MRI-based head model.

One-Shot Fabrication of Polymeric Hollow Microneedles by Standard Photolithography.

Microneedles (MNs) are an emerging technology in pharmaceutics and biomedicine, and are ready to be commercialized in the world market. However, solid microneedles only allow small doses and time-limited administration rates. Moreover, some well-known and already approved drugs need to be re-formulated when supplied by MNs. Instead, hollow microneedles (HMNs) allow for rapid, painless self-administrable microinjection of drugs in their standard formulation. Furthermore, body fluids can be easily extracted for analysis by a reverse use of HMNs, thus making them perfect for sensing issues and theranostics applications. The fabrication of HMNs usually requires several many-step processes, increasing the costs and consequently decreasing the commercial interest. Photolithography is a well-known fabrication technique in microelectronics and microfluidics that fabricates MNs. In this paper, authors show a proof of concept of a patented, easy and one-shot fabrication of two kinds of HMNs: (1) Symmetric HMNs with a “volcano” shape, made by using a photolithographic mask with an array of transparent symmetric rings; and (2) asymmetric HMNs with an oblique aperture, like standard hypodermic steel needles, made by using an array of transparent asymmetric rings, defined by two circles, which centers are slightly mismatched. Simulation of light propagation, fabrication process, and preliminary results on ink microinjection are presented.

A novel design of all optical half-subtractor using a square lattice photonic crystals

In this research, an all-optical half-subtractor is designed and simulated using two-dimensional photonic crystals. First, a photonic crystal structure is created using Si rods in the air context to obtain the optical half-subtractor. Afterward, using point and line defects, two waveguides are created for the input and two waveguides are created for the outputs. A high logical value and a low logical value are defined based on the optical power in each port. The FDTD method is used in the simulation of light propagation in the structure. The simulation results show that the designed half-subtractor has high optical power values for logic “1” and low values for logic “0”. The small size of the designed structure is among the advantages of this structure. Moreover, given that this half-subtractor is devoid of ring resonators, it can be used in high-speed integrated optical circuits. Another advantage of the proposed half-subtractor is that the optical powers in the outputs are similar in the high and low optical states.

Probe Design Optimization for Time-Domain NIROT "Pioneer" System for Imaging the Oxygenation of the Preterm Brain.

In preterm infants, there is a risk of life-lasting impairments due to hemorrhagic/ischemic lesions. Our time-domain (TD) near-infrared optical tomography (NIROT) system "Pioneer" aims at detecting both disorders with high spatial resolution. Successfully tested on phantoms, "Pioneer" entered the phase of improvements and enhancements. The current probe (A-probe) was adapted for an optoacoustics instrument. A new probe (B-probe) optimized for TD measurements is required. Our aim is to determine the optimal arrangement of light sources in the B-probe to increase the sensitivity and the resolution of Pioneer and to improve the ability of the system to detect both ischemia and hemorrhage. To do this, we simulated TD-NIROT signals in NIRFAST, a MATLAB-based package used to model near-infrared light propagation through tissue. We used 16×16 detector array, with ~2.2 mm distance between the detectors. Light sources were arranged around the field of view (FoV). We performed forward simulations of light propagation through a "homogeneous case" (HC) tissue (μ's=5.6 cm-1, μa=0.07 cm-1). Next, we simulated light propagation through "inhomogeneous case" -tissue' (IC) tissue by adding ischemia (μa=μa·2.5 cm-1) or hemorrhage (μa=μa·50 cm-1) to HT as a spherical inclusion of 5 mm radius at different depths in the FoV center and identified the source location that provides the higher contrast on the FoV: maxi∈I (FoVContrastSOURCE). It was found that sources located closer to the FoV center generate greater contrast for late photons. This study suggests the light sources in B-probe should be closer to the FoV center. The higher sensitivity is expected to lead to a higher image quality.

Face Illumination Transfer and Swapping via Dense Landmark and Semantic Parsing

The image-based virtual illumination technology directly changes the illumination effect of objects on the image. It does not require complex light propagation simulation calculations, it is an image-based rendering technology. The image mainly relies on the imaging of the visual sensor, and uses the virtual illumination technology of the image to illuminate the image obtained by the visual sensor. This is a new cross-cutting direction in computer vision, virtual reality and other fields. Face Illumination Swapping via Dense Landmark and Semantic Parsing is a major branch. Keeping the geometrical features of the target images and relighting the entire image instead of the face area are problems to be solved in the research. This paper based on the three-dimensional model to analyze the illumination information of the face images and re-render the illumination of the target face, and finally achieves the illumination swap between the two face images. We designed and implemented a 3DDFA-based face image illumination transfer method. First, 3DDFA is used to reconstruct the target face image. Estimate the surface normal and albedo. Then align and fill the surface normal and face parsing to illuminate the face image for light rendering and illumination transfer of the face images. Finally, the illumination analysis and re-rendering of face images based on 3DDFA are expanded to achieve the swap of illumination between face images. Experimental results show that this method can generate good effects of face image illumination transmission and swap while keeping the geometric features of the target images.

Light propagation within N95 filtered face respirators: A simulation study for UVC decontamination.

This study presents numerical simulations of UVC light propagation through seven different filtered face respirators (FFR) to determine their suitability for Ultraviolet germicidal inactivation (UVGI). UV propagation was modeled using the FullMonte program for two external light illuminations. The optical properties of the dominant three layers were determined using the inverse adding doubling method. The resulting fluence rate volume histograms and the lowest fluence rate recorded in the modeled volume, sometimes in the nW cm-2 , provide feedback on a respirator's suitability for UVGI and the required exposure time for a given light source. While UVGI can present an economical approach to extend an FFR's useable lifetime, it requires careful optimization of the illumination setup and selection of appropriate respirators.

Simulation of light propagation in citrus fruit using monte carlo multi-layered (MCML) method

The interaction of visible and near-infrared light with citrus tissues is of great importance for the design of new optical sensors for fruit quality inspection. In this study, the light propagation through citrus fruit was simulated using the Monte Carlo multi-layered (MCML) technique based on the bulk optical properties (BOP) of citrus tissues. Comparison of the diffuse reflectance profiles from MCML and from the laser scatter images confirmed the accuracy of the simulations. Simulation results indicated that most of the light was diffusely reflected or absorbed in the flavedo. Only a small fraction around 800 nm was transmitted through the three layers. Flavedo and albedo were found to contribute mostly to the diffuse reflectance spectra at 800 nm and their contributions respectively decreased and increased with source-detector distance. Although the contribution of juice vesicles increased with distance, the maximal diffuse reflectance intensity contributed by juice vesicles was located within 0.56 cm. Based on these results, the combination of the power intensity and the source-detector distance can be optimized to obtain more information about the layer of interest, which lays the foundations for designing new quality detection sensors for citrus fruit.

TLD1433-Mediated Photodynamic Therapy with an Optical Surface Applicator in the Treatment of Lung Cancer Cells In Vitro.

Intra-operative photodynamic therapy (IO-PDT) in combination with surgery for the treatment of non-small cell lung cancer and malignant pleural mesothelioma has shown promise in improving overall survival in patients. Here, we developed a PDT platform consisting of a ruthenium-based photosensitizer (TLD1433) activated by an optical surface applicator (OSA) for the management of residual disease. Human lung adenocarcinoma (A549) cell viability was assessed after treatment with TLD1433-mediated PDT illuminated with either 532- or 630-nm light with a micro-lens laser fiber. This TLD1433-mediated PDT induced an EC50 of 1.98 μM (J/cm2) and 4807 μM (J/cm2) for green and red light, respectively. Cells were then treated with 10 µM TLD1433 in a 96-well plate with the OSA using two 2-cm radial diffusers, each transmitted 532 nm light at 50 mW/cm for 278 s. Monte Carlo simulations of the surface light propagation from the OSA computed light fluence (J/cm2) and irradiance (mW/cm2) distribution. In regions where 100% loss in cell viability was measured, the simulations suggest that >20 J/cm2 of 532 nm was delivered. Our studies indicate that TLD1433-mediated PDT with the OSA and light simulations have the potential to become a platform for treatment planning for IO-PDT.

Minimal required PDT light dosimetry for nonmuscle invasive bladder cancer.

.Significance: Photodynamic therapy (PDT) could become a treatment option for nonmuscle invasive bladder cancer when the current high morbidity rate associated with red light PDT and variable PDT dose can be overcome through a combination of intravesical instillation of the photosensitizer and the use of green light creating a steep PDT dose gradient.Aim: To determine how a high PDT selectivity can be maintained throughout the bladder wall considering other efficacy determining parameters, in particular, the average optical properties of the mucosal layer governing the fluence rate multiplication factor, as well as the bladder shape and the position of the emitter in relationship to the bladder wall.Approach: We present three irradiance monitoring systems and evaluate their ability to enable selective bladder PDT considering previously determined photodynamic threshold values for the bladder cancer, mucosa and urothelium in a preclinical model, and the photosensitizer’s specific uptake ratio. Monte Carlo-based light propagation simulations performed for six human bladders at the time of therapy for a range of tissue optical properties. The performance of one irradiance sensing device in a clinical phase 1B trial is presented to underline the impact of irradiance monitoring, and it is compared to the Monte Carlo-derived dose surface histogram.Results: Monte Carlo simulations showed that irradiance monitoring systems need to comprise at least three sensors. Light scattering inside the bladder void needs to be minimized to prevent increased heterogeneity of the irradiance. The dose surface histograms vary significantly depending on the bladder shape and bladder volume but are less dependent on tissue optical properties.Conclusions: We demonstrate the need for adequate irradiance monitoring independent of a photosensitizer’s specific uptake ratio.

A complete optical model for liquid-scintillator detectors

Liquid scintillator (LS) is widely used in various neutrino oscillation experiments, in particular, the reactor neutrino experiments. The complex absorption and re-emission processes of optical photons are known to be an important source of the non-linear and non-uniform response of LS detectors. Precise simulation of light propagation in LS is highly desirable to model the detector response and reduce the systematic errors. In this paper, we develop a novel optical model which can completely deal with the competitive photon absorption and subsequent re-emission processes of the LS components. It allows to directly plug in the laboratory measurements of the LS components to model any LS composition. Extensive measurements have been performed to obtain the essential optical parameters for this model. We validate the model with a bench-top experiment featuring a small LS volume. Furthermore, we demonstrate that for any given detector geometry, this model provides the capability of optimizing the LS recipe to maximize the light collection. It is valuable for designing future LS-based detectors and improving the agreement between Monte Carlo and data for current neutrino experiments.

Simulation of Mach-Zehnder Interferometer plastic optical fibre for temperature and refractive index measurement

This work presents simulation of light propagation inside tapered structure-plastic optical fibre (POF) subjected to various surrounding temperatures and refractive index (RI). The tapered structure served as Mach-Zehnder interferometer (MZI) which split light into core modes and clad modes that recombined at the end of the taper. Simulation was done by using Beam Propagation Method (BPM) solved using Finite Difference Method (FDM). The light launched in to the MZI was varied from 611nm to 655nm. In order to give temperature effect, the temperature was varied from 30°C to 80°C with increment of 10°C. The results showed that the free spectral range (FSR) of the MZI are shifted as temperature increased. The sensitivity for the peak and dip of the FSR are -0.0158 nm/°C and -0.0153 nm/°C, respectively. For RI sensitivity investigation, the surrounding RI was varied from 1.333 to 1.349. The sensitivity of peak FSR and dip FSR to RI change are -98.638 nm/RIU and 105.4 nm/RIU. Therefore, the MZI can be used to measure RI and temperature simultaneously.