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

2021 ◽  
Vol 51 (5) ◽  
pp. 408-414
Author(s):  
A P Tarasov ◽  
S Persheyev ◽  
D A Rogatkin
Keyword(s):  
Monte Carlo ◽  
Light Transport ◽  
Monte Carlo Algorithm ◽  
Biological Media ◽  
Continuous Absorption

scholarly journals A GAMOS plug-in for GEANT4 based Monte Carlo simulation of radiation-induced light transport in biological media

2013 ◽  
Vol 4 (5) ◽  
pp. 741 ◽  
Author(s):  
Adam K. Glaser ◽  
Stephen C. Kanick ◽  
Rongxiao Zhang ◽  
Pedro Arce ◽  
Brian W. Pogue
Keyword(s):  
Monte Carlo Simulation ◽  
Monte Carlo ◽  
Light Transport ◽  
Biological Media ◽  
Radiation Induced

Parallel Simulation of Steady State Light Transport Through Solid Media Using Logistic Map

2020 ◽  
Vol 17 (4) ◽  
pp. 1606-1609
Author(s):  
K. Sathish ◽  
Siddharth Singh ◽  
Anirudh Agarwal ◽  
Pranab Kumar Shukla
Keyword(s):  
Monte Carlo ◽  
Photodynamic Therapy ◽  
Gpu Computing ◽  
Logistic Map ◽  
Parallel Simulation ◽  
Light Transport ◽  
Monte Carlo Algorithm ◽  
Photon Transport ◽  
Solid Media ◽  
Wide Range

The Monte Carlo algorithm has been extensively used for photon transport simulations in medical imaging to assists doctors in Photodynamic Therapy for treatment of wide range of medical conditions including varieties of cancer by eliciting phototoxicity in cells. Previously this was done using static 2 Dimensional models on traditional CPUs. With the advent of GPU Computing further work was done to extend this model by separating the PRNG.

scholarly journals MCMLpar and MCSLinv: A Parallel Version of MCML and an Inverse Monte Carlo Algorithm to Calculate Optical Scattering Parameters

2017 ◽  
Vol 122 ◽  
Author(s):  
Richelle H. Streater ◽  
Anne-Michelle R. Lieberson ◽  
Adam L. Pintar ◽  
Zachary H. Levine
Keyword(s):  
Monte Carlo ◽  
Quadratic Function ◽  
Search Space ◽  
Monte Carlo Sampling ◽  
Optical Scattering ◽  
Light Transport ◽  
Monte Carlo Algorithm ◽  
Optimum Parameters ◽  
Inverse Monte Carlo ◽  
Log Likelihood

The MCML program for Monte Carlo modeling of light transport in multi-layered tissues has been widely used in the past 20 years or so. Here, we have re-implemented MCML for solving the inverse problem. Our formulation features optimizing the profile log likelihood which takes into account uncertainties due to both experimental and Monte Carlo sampling. We limit the search space for the optimum parameters with relatively few Monte Carlo trials and then iteratively double the number of Monte Carlo trials until the search space stabilizes. At this point, the log likelihood can be fit with a quadratic function to find the optimum. The time-to-solution is only a few minutes in typical cases because we use importance sampling to determine the log likelihood on a grid of parameters at each iteration. Also, our implementation uses OpenMP and SPRNG to generate Monte Carlo trials in parallel.

scholarly journals A hybrid mesh and voxel based Monte Carlo algorithm for accurate and efficient photon transport modeling in complex bio-tissues

2020 ◽  
Author(s):  
Shijie Yan ◽  
Qianqian Fang
Keyword(s):  
Monte Carlo ◽  
Graphics Processing Units ◽  
Light Transport ◽  
Monte Carlo Algorithm ◽  
Brain Atlas ◽  
Simulation Accuracy ◽  
Mc Simulation ◽  
Voxel Data ◽  
Graphics Processing ◽  
Improved Accuracy

AbstractOver the past decade, an increasing body of evidence has suggested that threedimensional (3-D) Monte Carlo (MC) light transport simulations are affected by the inherent limitations and errors of voxel-based domain boundaries. In this work, we specifically address this challenge using a hybrid MC algorithm, namely split-voxel MC or SVMC, that combines both mesh and voxel domain information to greatly improve MC simulation accuracy while remaining highly flexible and efficient in parallel hardware, such as graphics processing units (GPU). We achieve this by applying a marching-cubes algorithm to a pre-segmented domain to extract and encode sub-voxel information of curved surfaces, which is then used to inform ray-tracing computation within boundary voxels. This preservation of curved boundaries in a voxel data structure demonstrates significantly improved accuracy in several benchmarks, including a human brain atlas. The accuracy of the SVMC algorithm is comparable to that of mesh-based MC (MMC), but runs 2x-6x faster and requires only a lightweight preprocessing step. The proposed algorithm has been implemented in our open-source software and is freely available at http://mcx.space.

scholarly journals Light transport modeling in highly complex tissues using implicit mesh-based Monte Carlo algorithm

2020 ◽  
Author(s):  
Yaoshen Yuan ◽  
Shijie Yan ◽  
Qianqian Fang
Keyword(s):  
Monte Carlo ◽  
Light Propagation ◽  
Light Transport ◽  
Monte Carlo Algorithm ◽  
Memory Usage ◽  
Microvascular Networks ◽  
Multi Scale ◽  
Tissue Structures ◽  
Vessel Networks ◽  
Shape Complexity

AbstractThe mesh-based Monte Carlo (MMC) technique has grown tremendously since its initial publication nearly a decade ago. It is now recognized as one of the most accurate Monte Carlo (MC) methods, providing accurate reference solutions for the development of novel biophotonics techniques. In this work, we aim to further advance MMC to address a major challenge in biophotonics modeling, i.e. light transport within highly complex tissues, such as dense microvascular networks, porous media and multi-scale tissue structures. Although the current MMC framework is capable of simulating light propagation in such media given its generality, the run-time and memory usage grow rapidly with increasing media complexity and size. This greatly limits our capability to explore complex and multi-scale tissue structures. Here, we propose a highly efficient implicit mesh-based Monte Carlo (iMMC) method that incorporates both mesh- and shape-based tissue representations to create highly complex yet memory efficient light transport simulations. We demonstrate that iMMC is capable of providing accurate solutions for dense vessel networks and porous tissues while reducing memory usage by greater than a hundred- or even thousand-fold. In a sample network of microvasculature, the reduced shape complexity results in nearly 3x speed acceleration. The proposed algorithm is now available in our open-source MMC software at http://mcx.space/#mmc.

scholarly journals A New Discrete Implicit Monte Carlo Scheme for Simulating Radiative Transfer Problems

2022 ◽  
Vol 258 (1) ◽  
pp. 14
Author(s):  
Elad Steinberg ◽  
Shay I. Heizler
Keyword(s):  
Monte Carlo ◽  
Radiative Transfer ◽  
Monte Carlo Algorithm ◽  
Basic Algorithm ◽  
Time Step ◽  
Radiative Shock ◽  
Continuous Absorption ◽  
Discrete Nature ◽  
The One ◽  
Transfer Problems

Abstract We present a new algorithm for radiative transfer—based on a statistical Monte Carlo approach—that does not suffer from teleportation effects, on the one hand, and yields smooth results, on the other hand. Implicit Monte Carlo (IMC) techniques for modeling radiative transfer have existed from the 1970s. When they are used for optically thick problems, however, the basic algorithm suffers from “teleportation” errors, where the photons propagate faster than the exact physical behavior, due to the absorption-blackbody emission processes. One possible solution is to use semianalog Monte Carlo, in its new implicit form (ISMC), which uses two kinds of particles, photons and discrete material particles. This algorithm yields excellent teleportation-free results, but it also produces noisier solutions (relative to classic IMC), due to its discrete nature. Here, we derive a new Monte Carlo algorithm, Discrete Implicit Monte Carlo (DIMC), which also uses the idea of two kinds of discrete particles, and thus does not suffer from teleportation errors. DIMC implements the IMC discretization and creates new radiation photons for each time step, unlike ISMC. Using the continuous absorption technique, DIMC yields smooth results like classic IMC. One of the main elements of the algorithm is the avoidance of the explosion of the particle population, by using particle merging. We test the new algorithm on 1D and 2D cylindrical problems, and show that it yields smooth, teleportation-free results. We finish by demonstrating the power of the new algorithm on a classic radiative hydrodynamic problem—an opaque radiative shock wave. This demonstrates the power of the new algorithm for astrophysical scenarios.

Sign in / Sign up

Export Citation Format

Share Document