GPU-based Monte Carlo simulation for light propagation in complex heterogeneous tissues

Dictionary Learning Method Based on K-Sparse Approximation and Orthogonal Procrustes Analysis for Reconstruction in Bioluminescence Tomography

Bioluminescence tomography (BLT) is one kind of noninvasive optical molecular imaging technology, widely used to study molecular activities and disease progression inside live animals. By combining the optical propagation model and inversion algorithm, BLT enables three-dimensional imaging and quantitative analysis of light sources within organisms. However, challenges like light scattering and absorption in tissues, and the complexity of biological structures, significantly impact the accuracy of BLT reconstructions. Here, we propose a dictionary learning method based on K-sparse approximation and Orthogonal Procrustes analysis (KSAOPA). KSAOPA uses an iterative alternating optimization strategy, enhancing solution sparsity with k-coefficients Lipschitzian mappings for sparsity(K-LIMAPS) in the sparse coding stage, and reducing errors with Orthogonal Procrustes analysis in the dictionary update stage, leading to stable and precise reconstructions. We assessed the method performance through simulations and in vivo experiments, which showed that KSAOPA excels in localization accuracy, morphological recovery, and in vivo applicability compared to other methods.

Logarithmic total variation regularization via preconditioned conjugate gradient method for sparse reconstruction of bioluminescence tomography

BACKGROUND AND OBJECTIVE

Bioluminescence Tomography (BLT) is a powerful optical molecular imaging technique that enables the noninvasive investigation of dynamic biological phenomena. It aims to reconstruct the three-dimensional spatial distribution of bioluminescent sources from optical measurements collected on the surface of the imaged object. However, BLT reconstruction is a challenging ill-posed problem due to the scattering effect of light and the limitations in detecting surface photons, which makes it difficult for existing methods to achieve satisfactory reconstruction results. In this study, we propose a novel method for sparse reconstruction of BLT based on a preconditioned conjugate gradient with logarithmic total variation regularization (PCG-logTV).

METHOD

This PCG-logTV method incorporates the sparsity of overlapping groups and enhances the sparse structure of these groups using logarithmic functions, which can preserve edge features and achieve more stable reconstruction results in BLT. To accelerate the convergence of the algorithm solution, we use the preconditioned conjugate gradient iteration method on the objective function and obtain the reconstruction results. We demonstrate the performance of our proposed method through numerical simulations and in vivo experiment.

RESULTS AND CONCLUSIONS

The results show that the PCG-logTV method obtains the most accurate reconstruction results, and the minimum position error (LE) is 0.254mm, which is 26%, 31% and 34% of the FISTA (0.961), IVTCG (0.81) and L1-TV (0.739) methods, and the root mean square error (RMSE) and relative intensity error (RIE) are the smallest, indicating that it is closest to the real light source. In addition, compared with the other three methods, the PCG-logTV method also has the highest DICE similarity coefficient, which is 0.928, which means that this method can effectively reconstruct the three-dimensional spatial distribution of bioluminescent light sources, has higher resolution and robustness, and is beneficial to the preclinical and clinical studies of BLT.

Regularized reconstruction based on joint smoothly clipped absolute deviation regularization and graph manifold learning for fluorescence molecular tomography

Objective.Fluorescence molecular tomography (FMT) is an optical imaging modality that provides high sensitivity and low cost, which can offer the three-dimensional distribution of biomarkers by detecting the fluorescently labeled probe noninvasively. In the field of preclinical cancer diagnosis and treatment, FMT has gained significant traction. Nonetheless, the current FMT reconstruction results suffer from unsatisfactory morphology and location accuracy of the fluorescence distribution, primarily due to the light scattering effect and the ill-posed nature of the inverse problem.Approach.To address these challenges, a regularized reconstruction method based on joint smoothly clipped absolute deviation regularization and graph manifold learning (SCAD-GML) for FMT is presented in this paper. The SCAD-GML approach combines the sparsity of the fluorescent sources with the latent manifold structure of fluorescent source distribution to achieve more accurate and sparse reconstruction results. To obtain the reconstruction results efficiently, the non-convex gradient descent iterative method is employed to solve the established objective function. To assess the performance of the proposed SCAD-GML method, a comprehensive evaluation is conducted through numerical simulation experiments as well asin vivoexperiments.Main results.The results demonstrate that the SCAD-GML method outperforms other methods in terms of both location and shape recovery of fluorescence biomarkers distribution.Siginificance.These findings indicate that the SCAD-GML method has the potential to advance the application of FMT inin vivobiological research.

A Monte Carlo simulation for Moving Light Source in Intracavity PDT

We developed a simulation method for modeling the light fluence delivery in intracavity Photodynamic Therapy (icav-PDT) for pleural lung cancer using a moving light source. Due to the large surface area of the pleural lung cavity, the light source needs to be moved to deliver a uniform dose around the entire cavity. While multiple fixed detectors are used for dosimetry at a few locations, an accurate simulation of light fluence and fluence rate is still needed for the rest of the cavity. We extended an existing Monte Carlo (MC) based light propagation solver to support moving light sources by densely sampling the continuous light source trajectory and assigning the proper number of photon packages launched along the way. The performance of Simphotek GPU CUDA-based implementation of the method - PEDSy-MC - has been demonstrated on a life-size lung-shaped phantom, custom printed for testing icav-PDT navigation system at the Perlman School of Medicine (PSM) - calculations completed under a minute (for some cases) and within minutes have been achieved. We demonstrate results within a 5% error of the analytic solution for multiple detectors in the phantom. PEDSy-MC is accompanied by a dose-cavity visualization tool that allows real-time inspection of dose values of the treated cavity in 2D and 3D, which will be expanded to ongoing clinical trials at PSM. PSM has developed a technology to measure 8-detectors in a pleural cavity phantom using Photofrin-mediated PDT that has been used during validation.

ABPO-TVSCAD: alternating Bregman proximity operators approach based on TVSCAD regularization for bioluminescence tomography

Objective.Bioluminescence tomography (BLT) is a promising non-invasive optical medical imaging technique, which can visualize and quantitatively analyze the distribution of tumor cells in living tissues. However, due to the influence of photon scattering effect and ill-conditioned inverse problem, the reconstruction result is unsatisfactory. The purpose of this study is to improve the reconstruction performance of BLT.Approach.An alternating Bregman proximity operators (ABPO) method based on TVSCAD regularization is proposed for BLT reconstruction. TVSCAD combines the anisotropic total variation (TV) regularization constraints and the non-convex smoothly clipped absolute deviation (SCAD) penalty constraints, to make a trade-off between the sparsity and edge preservation of the source. ABPO approach is used to solve the TVSCAD model (ABPO-TVSCAD for short). In addition, to accelerate the convergence speed of the ABPO, we adapt the strategy of shrinking the permission source region, which further improves the performance of ABPO-TVSCAD.Main results.The results of numerical simulations andin vivoxenograft mouse experiment show that our proposed method achieved superior accuracy in spatial localization and morphological reconstruction of bioluminescent source.Significance.ABPO-TVSCAD is an effective and robust reconstruction method for BLT, and we hope that this method can promote the development of optical molecular tomography.

Adaptive permissible region based random Kaczmarz reconstruction method for localization of carotid atherosclerotic plaques in fluorescence molecular tomography

Objective. In this study, we propose the adaptive permissible region based random Kaczmarz method as an improved reconstruction method to recover small carotid atherosclerotic plaque targets in rodents with high resolution in fluorescence molecular tomography (FMT). Approach. We introduce the random Kaczmarz method as an advanced minimization method to solve the FMT inverse problem. To satisfy the special condition of this method, we proposed an adaptive permissible region strategy based on traditional permissible region methods to flexibly compress the dimension of the solution space. Main results. Monte Carlo simulations, phantom experiments, and in vivo experiments demonstrate that the proposed method can recover the small carotid atherosclerotic plaque targets with high resolution and accuracy, and can achieve lower root mean squared error and distance error (DE) than other traditional methods. For targets with 1.5 mm diameter and 0.5 mm separation, the DE indicators can be improved by up to 40%. Moreover, the proposed method can be utilized for in vivo locating atherosclerotic plaques with high accuracy and robustness. Significance. We applied the random Kaczmarz method to solve the inverse problem in FMT and improve the reconstruction result via this advanced minimization method. We verified that the FMT technology has a great potential to locate and quantify atherosclerotic plaques with higher accuracy, and can be expanded to more preclinical research.

Rapid Quantification of Tissue Perfusion Properties with a Two-Stage Look-Up Table

Tissue perfusion properties reveal crucial information pertinent to clinical diagnosis and treatment. Multispectral spatial frequency domain imaging (SFDI) is an emerging imaging technique that has been widely used to quantify tissue perfusion properties. However, slow processing speed limits its usefulness in real-time imaging applications. In this study, we present a two-stage look-up table (LUT) approach that accurately and rapidly quantifies optical (absorption and reduced scattering maps) and perfusion (total hemoglobin and oxygen saturation maps) properties using stage-1 and stage-2 LUTs, respectively, based on reflectance images at 660 and 850 nm. The two-stage LUT can be implemented on both CPU and GPU computing platforms. Quantifying tissue perfusion properties using the simulated diffuse reflectance images, we achieved a quantification speed of 266, 174, and 74 frames per second for three image sizes 512 × 512, 1024 × 1024, and 2048 × 2048 pixels, respectively. Quantification of tissue perfusion properties was highly accurate with only 3.5% and 2.5% error for total hemoglobin and oxygen saturation quantification, respectively. The two-stage LUT has the potential to be integrated with dual-sensor cameras to enable real-time quantification of tissue hemodynamics.

Two-step verification method for Monte Carlo codes in biomedical optics applications

SIGNIFICANCE

Code verification is an unavoidable step prior to using a Monte Carlo (MC) code. Indeed, in biomedical optics, a widespread verification procedure for MC codes is still missing. Analytical benchmarks that can be easily used for the verification of different MC routines offer an important resource.

AIM

We aim to provide a two-step verification procedure for MC codes enabling the two main tasks of an MC simulator: (1) the generation of photons' trajectories and (2) the intersections of trajectories with boundaries separating the regions with different optical properties. The proposed method is purely based on elementary analytical benchmarks, therefore, the correctness of an MC code can be assessed with a one-sample t-test.

APPROACH

The two-step verification is based on the following two analytical benchmarks: (1) the exact analytical formulas for the statistical moments of the spatial coordinates where the scattering events occur in an infinite medium and (2) the exact invariant solutions of the radiative transfer equation for radiance, fluence rate, and mean path length in media subjected to a Lambertian illumination.

RESULTS

We carried out a wide set of comparisons between MC results and the two analytical benchmarks for a wide range of optical properties (from non-scattering to highly scattering media, with different types of scattering functions) in an infinite non-absorbing medium (step 1) and in a non-absorbing slab (step 2). The deviations between MC results and exact analytical values are usually within two standard errors (i.e., t-tests not rejected at a 5% level of significance). The comparisons show that the accuracy of the verification increases with the number of simulated trajectories so that, in principle, an arbitrary accuracy can be obtained.

CONCLUSIONS

Given the simplicity of the verification method proposed, we envision that it can be widely used in the field of biomedical optics.

MCDataset: a public reference dataset of Monte Carlo simulated quantities for multilayered and voxelated tissues computed by massively parallel PyXOpto Python package

SIGNIFICANCE

Current open-source Monte Carlo (MC) method implementations for light propagation modeling are many times tedious to build and require third-party licensed software that can often discourage prospective researchers in the biomedical optics community from fully utilizing the light propagation tools. Furthermore, the same drawback also limits rigorous cross-validation of physical quantities estimated by various MC codes.

AIM

Proposal of an open-source tool for light propagation modeling and an easily accessible dataset to encourage fruitful communications amongst researchers and pave the way to a more consistent comparison between the available implementations of the MC method.

APPROACH

The PyXOpto implementation of the MC method for multilayered and voxelated tissues based on the Python programming language and PyOpenCL extension enables massively parallel computation on numerous OpenCL-enabled devices. The proposed implementation is used to compute a large dataset of reflectance, transmittance, energy deposition, and sampling volume for various source, detector, and tissue configurations.

RESULTS

The proposed PyXOpto agrees well with the original MC implementation. However, further validation reveals a noticeable bias introduced by the random number generator used in the original MC implementation.

CONCLUSIONS

Establishing a common dataset is highly important for the validation of existing and development of MC codes for light propagation in turbid media.

Scalable and accessible personalized photodynamic therapy optimization with FullMonte and PDT-SPACE

SIGNIFICANCE

Open-source software packages have been extensively used in the past three decades in medical imaging and diagnostics, aiming to study the feasibility of the application ex vivo. Unfortunately, most of the existing open-source tools require some software engineering background to install the prerequisite libraries, choose a suitable computational platform, and combine several software tools to address different applications.

AIM

To facilitate the use of open-source software in medical applications, enabling computational studies of treatment outcomes prior to the complex in-vivo setting.

APPROACH

FullMonteWeb, an open-source, user-friendly web-based software with a graphical user interface for interstitial photodynamic therapy (iPDT) modeling, visualization, and optimization, is introduced. The software can perform Monte Carlo simulations of light propagation in biological tissues, along with iPDT plan optimization. FullMonteWeb installs and runs the required software and libraries on Amazon Web Services (AWS), allowing scalable computing without complex set up.

RESULTS

FullMonteWeb allows simulation of large and small problems on the most appropriate compute hardware, enabling cost improvements of 10  ×   versus always running on a single platform. Case studies in optical property estimation and diffuser placement optimization highlight FullMonteWeb's versatility.

CONCLUSION

The FullMonte open source suite enables easier and more cost-effective in-silico studies for iPDT.

Cuffless Blood Pressure Estimation Based on Monte Carlo Simulation Using Photoplethysmography Signals

Blood pressure measurements are one of the most routinely performed medical tests globally. Blood pressure is an important metric since it provides information that can be used to diagnose several vascular diseases. Conventional blood pressure measurement systems use cuff-based devices to measure the blood pressure, which may be uncomfortable and sometimes burdensome to the subjects. Therefore, in this study, we propose a cuffless blood pressure estimation model based on Monte Carlo simulation (MCS). We propose a heterogeneous finger model for the MCS at wavelengths of 905 nm and 940 nm. After recording the photon intensities from the MCS over a certain range of blood pressure values, the actual photoplethysmography (PPG) signals were used to estimate blood pressure. We used both publicly available and self-made datasets to evaluate the performance of the proposed model. In case of the publicly available dataset for transmission-type MCS, the mean absolute errors are 3.32 ± 6.03 mmHg for systolic blood pressure (SBP), 2.02 ± 2.64 mmHg for diastolic blood pressure (DBP), and 1.76 ± 2.8 mmHg for mean arterial pressure (MAP). The self-made dataset is used for both transmission- and reflection-type MCSs; its mean absolute errors are 2.54 ± 4.24 mmHg for SBP, 1.49 ± 2.82 mmHg for DBP, and 1.51 ± 2.41 mmHg for MAP in the transmission-type case as well as 3.35 ± 5.06 mmHg for SBP, 2.07 ± 2.83 mmHg for DBP, and 2.12 ± 2.83 mmHg for MAP in the reflection-type case. The estimated results of the SBP and DBP satisfy the requirements of the Association for the Advancement of Medical Instrumentation (AAMI) standards and are within Grade A according to the British Hypertension Society (BHS) standards. These results show that the proposed model is efficient for estimating blood pressures using fingertip PPG signals.

Tutorial on the Use of Deep Learning in Diffuse Optical Tomography

Diffuse optical tomography using deep learning is an emerging technology that has found impressive medical diagnostic applications. However, creating an optical imaging system that uses visible and near-infrared (NIR) light is not straightforward due to photon absorption and multi-scattering by tissues. The high distortion levels caused due to these effects make the image reconstruction incredibly challenging. To overcome these challenges, various techniques have been proposed in the past, with varying success. One of the most successful techniques is the application of deep learning algorithms in diffuse optical tomography. This article discusses the current state-of-the-art diffuse optical tomography systems and comprehensively reviews the deep learning algorithms used in image reconstruction. This article attempts to provide researchers with the necessary background and tools to implement deep learning methods to solve diffuse optical tomography.

MCX Cloud-a modern, scalable, high-performance and in-browser Monte Carlo simulation platform with cloud computing

SIGNIFICANCE

Despite the ample progress made toward faster and more accurate Monte Carlo (MC) simulation tools over the past decade, the limited usability and accessibility of these advanced modeling tools remain key barriers to widespread use among the broad user community.

AIM

An open-source, high-performance, web-based MC simulator that builds upon modern cloud computing architectures is highly desirable to deliver state-of-the-art MC simulations and hardware acceleration to general users without the need for special hardware installation and optimization.

APPROACH

We have developed a configuration-free, in-browser 3D MC simulation platform-Monte Carlo eXtreme (MCX) Cloud-built upon an array of robust and modern technologies, including a Docker Swarm-based cloud-computing backend and a web-based graphical user interface (GUI) that supports in-browser 3D visualization, asynchronous data communication, and automatic data validation via JavaScript Object Notation (JSON) schemas.

RESULTS

The front-end of the MCX Cloud platform offers an intuitive simulation design, fast 3D data rendering, and convenient simulation sharing. The Docker Swarm container orchestration backend is highly scalable and can support high-demand GPU MC simulations using MCX over a dynamically expandable virtual cluster.

CONCLUSION

MCX Cloud makes fast, scalable, and feature-rich MC simulations readily available to all biophotonics researchers without overhead. It is fully open-source and can be freely accessed at http://mcx.space/cloud.

Modeling optical design parameters for fine stimulation in sciatic nerve of optogenetic mice

Optogenetics presents an alternative method for interfacing with the nervous system over the gold-standard of electrical stimulation. While electrical stimulation requires electrodes to be surgically embedded in tissue for in vivo studies, optical stimulation offers a less-invasive approach that may yield more specific, localized stimulation. The advent of optogenetic laboratory animals-whose motor neurons can be activated when illuminated with blue light-enables research into refining optical stimulation of the mammalian nervous system where subsets of nerve fibers within a nerve may be stimulated without embedding any device directly into the nerve itself. However, optical stimulation has a major drawback in that light is readily scattered and absorbed in tissue thereby limiting the depth with which a single emission source can penetrate. We hypothesize that the use of multiple, focused light emissions deployed around the circumference of a nerve can overcome these light-scattering limitations. To understand the physical parameters necessary to produce pinpointed light stimulation within a single nerve, we employed a simplified Monte Carlo simulation to estimate the size of nerves where this technique may be successful, as well as the necessary optical lens design for emitters to be used during future in vivo studies. By modeling multiple focused beams, we find that only fascicles within a nerve diameter less than 1 mm are fully accessible to focused optical stimulation; a minimum of 4 light sources is required to generate a photon intensity at a point in a nerve over the initial contact along its surface. To elicit the same effect in larger nerves, focusing lenses would require a numerical aperture [Formula: see text]. These simulations inform on the design of instrumentation capable of stimulating disparate motor neurons in mouse sciatic nerve to control hindlimb movement.

Rapid quantification of tissue perfusion properties with a two-stage look-up table: a simulation study.. [DOI: 10.1101/2021.11.04.467306] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Tissue perfusion properties reveal crucial information pertinent to clinical diagnosis and treatment. Multispectral spatial frequency domain imaging (SFDI) is an emerging imaging technique that has been widely used to quantify tissue perfusion properties. However, slow processing speed limits its usefulness in real-time imaging applications. In this study, we present a two-stage look-up table (LUT) approach that accurately and rapidly quantifies optical (absorption and reduced scattering maps) and perfusion (total hemoglobin and oxygen saturation maps) properties using stage-1 and stage-2 LUTs, respectively, based on reflectance images at 660nm and 850nm. The two-stage LUT can be implemented on both CPU and GPU computing platforms. Quantifying tissue perfusion properties using the simulated diffuse reflectance images, we achieved a quantification speed of 266, 174, and 74 frames per second for three image sizes 512×512, 1024×1024, and 2048×2048 pixels, respectively. Quantification of tissue perfusion properties was highly accurate with only 3.5% and 2.5% error for total hemoglobin and oxygen saturation quantification, respectively. The two-stage LUT has the potential to be adopted in existing SFDI applications to enable real-time imaging capability of tissue hemodynamics.

Noninvasive In Vivo Estimation of Blood-Glucose Concentration by Monte Carlo Simulation

Continuous monitoring of blood-glucose concentrations is essential for both diabetic and nondiabetic patients to plan a healthy lifestyle. Noninvasive in vivo blood-glucose measurements help reduce the pain of piercing human fingertips to collect blood. To facilitate noninvasive measurements, this work proposes a Monte Carlo photon simulation-based model to estimate blood-glucose concentration via photoplethysmography (PPG) on the fingertip. A heterogeneous finger model was exposed to light at 660 nm and 940 nm in the reflectance mode of PPG via Monte Carlo photon propagation. The bio-optical properties of the finger model were also deduced to design the photon simulation model for the finger layers. The intensities of the detected photons after simulation with the model were used to estimate the blood-glucose concentrations using a supervised machine-learning model, XGBoost. The XGBoost model was trained with synthetic data obtained from the Monte Carlo simulations and tested with both synthetic and real data (n = 35). For testing with synthetic data, the Pearson correlation coefficient (Pearson’s r) of the model was found to be 0.91, and the coefficient of determination (R²) was found to be 0.83. On the other hand, for tests with real data, the Pearson’s r of the model was 0.85, and R² was 0.68. Error grid analysis and Bland–Altman analysis were also performed to confirm the accuracy. The results presented herein provide the necessary steps for noninvasive in vivo blood-glucose concentration estimation.

Spatial-Frequency Domain Imaging: An Emerging Depth-Varying and Wide-Field Technique for Optical Property Measurement of Biological Tissues

Measurement of optical properties is critical for understanding light-tissue interaction, properly interpreting measurement data, and gaining better knowledge of tissue physicochemical properties. However, conventional optical measuring techniques are limited in point measurement, which partly hinders the applications on characterizing spatial distribution and inhomogeneity of optical properties of biological tissues. Spatial-frequency domain imaging (SFDI), as an emerging non-contact, depth-varying and wide-field optical imaging technique, is capable of measuring the optical properties in a wide field-of-view on a pixel-by-pixel basis. This review first describes the typical SFDI system and the principle for estimating optical properties using the SFDI technique. Then, the applications of SFDI in the fields of biomedicine, as well as food and agriculture, are reviewed, including burn assessment, skin tissue evaluation, tumor tissue detection, brain tissue monitoring, and quality evaluation of agro-products. Finally, a discussion on the challenges and future perspectives of SFDI for optical property estimation is presented.

Adaptive Grouping Block Sparse Bayesian Learning Method for Accurate and Robust Reconstruction in Bioluminescence Tomography

OBJECTIVE

Bioluminescence tomography (BLT) is a promising modality that is designed to provide non-invasive quantitative three-dimensional information regarding the tumor distribution in living animals. However, BLT suffers from inferior reconstructions due to its ill-posedness. This study aims to improve the reconstruction performance of BLT.

METHODS

We propose an adaptive grouping block sparse Bayesian learning (AGBSBL) method, which incorporates the sparsity prior, correlation of neighboring mesh nodes, and anatomical structure prior to balance the sparsity and morphology in BLT. Specifically, an adaptive grouping prior model is proposed to adjust the grouping according to the intensity of the mesh nodes during the optimization process.

RESULTS

Numerical simulations and in vivo experiments demonstrate that AGBSBL yields a high position and morphology recovery accuracy, stability, and practicality.

CONCLUSION

The proposed method is a robust and effective reconstruction algorithm for BLT. Moreover, the proposed adaptive grouping strategy can further increase the practicality of BLT in biomedical applications.

Fourier transform infrared absorption (FTIR) on dry stratum corneum, corneocyte-lipid interfaces: experimental and vibrational spectroscopy calculations

Many questions concerning the biophysical and physiological properties of skin are still open. Skin aging, permeability, dermal absorption, hydration, and drug transdermal delivery, are few examples of processes with unveiled underlying mechanisms. In this work, it is presented a comparison between Fourier transform infrared absorption (FTIR) of dry stratum corneum and stratum corneum under lipase action supported by first-principles density functional vibrational calculations. The molecular structure of stratum corneum was modeled by an archetype of its hygroscopic proteic portion inside the corneocytes, the natural moisturizing factor, coupled to glycerol molecules which represent the lipid fraction of stratum corneum. Vibrational spectra were calculated and compared to experimental data obtained on the animal model of stratum corneum. The experimental results indicated prominent spectral differences between dry and lipase-treated stratum corneum. Principal components analysis and hyerarchical clustering indicated that 1200, 1650, and 1695 cm^-1 bands are the most influential on the discrimination. It is noticed that bands in the fingerprint region (800-1800 cm^-1) were correctly assigned. Moreover, the calculations revealed the existence of two coupled vibration between the hydroxyl group of lipid and methylene (1120 and 1160 cm^-1), which are of special interest since they probe the lipid-amino acid coupling. The model was also able to predict the shear modulus of dry stratum corneum in excellent agreement with the reported values from the literature. Other physical/chemical properties could be calculated exploring the chemical accuracy and molecular resolution of this model. Research in dermatology, cosmetology, and biomedical engineering in the specific topics of drug delivery and/or mechanical properties of skin are examples of fields that would potentially take advantage of this approach.

GPU-Accelerated Flexible Molecular Docking

Virtual screening is a key enabler of computational drug discovery and requires accurate and efficient structure-based molecular docking. In this work, we develop algorithms and software building blocks for molecular docking that can take advantage of graphics processing units (GPUs). Specifically, we focus on MedusaDock, a flexible protein-small molecule docking approach and platform. We accelerate the performance of the coarse docking phase of MedusaDock, as this step constitutes nearly 70% of total running time in typical use-cases. We perform a comprehensive evaluation of the quality and performance with single-GPU and multi-GPU acceleration using a data set of 3875 protein-ligand complexes. The algorithmic ideas, data structure design choices, and performance optimization techniques shed light on GPU acceleration of other structure-based molecular docking software tools.

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

The 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.

Hybrid mesh and voxel based Monte Carlo algorithm for accurate and efficient photon transport modeling in complex bio-tissues

Over the past decade, an increasing body of evidence has suggested that three-dimensional (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.

Parallel Monte Carlo simulation algorithm for the spectral reflectance and transmittance of the wind-generated bubble layers in the upper ocean using CUDA

The parallel Monte Carlo software CUDAMCML used in the bio-optics field was developed by Erik Alerstam et al. (J. Biomed. Opt., 13, 060504, 2008) based on the Compute Unified Device Architecture (CUDA) and can simulate light transport in multilayered media. In the present study, CUDAMCML is extended to form the new program CUDAMCML-OCEAN using the average sampling method. This new program can handle multiple types of particle seawater containing elements such as colored dissolved organic matter (CDOM) and bubbles. The accuracy and speedup of the new program are analyzed. The results show that when the parameters are set appropriately, the speedup of CUDAMCML-OCEAN is more than 200 times compared with serial code. And the accuracies of the spectral reflectance and transmittance all reached a satisfactory level for different wind speeds and chlorophyll concentrations.

Efficient inversion strategies for estimating optical properties with Monte Carlo radiative transport models

SIGNIFICANCE

Indirect imaging problems in biomedical optics generally require repeated evaluation of forward models of radiative transport, for which Monte Carlo is accurate yet computationally costly. We develop an approach to reduce this bottleneck, which has significant implications for quantitative tomographic imaging in a variety of medical and industrial applications.

AIM

Our aim is to enable computationally efficient image reconstruction in (hybrid) diffuse optical modalities using stochastic forward models.

APPROACH

Using Monte Carlo, we compute a fully stochastic gradient of an objective function for a given imaging problem. Leveraging techniques from the machine learning community, we then adaptively control the accuracy of this gradient throughout the iterative inversion scheme to substantially reduce computational resources at each step.

RESULTS

For example problems of quantitative photoacoustic tomography and ultrasound-modulated optical tomography, we demonstrate that solutions are attainable using a total computational expense that is comparable to (or less than) that which is required for a single high-accuracy forward run of the same Monte Carlo model.

CONCLUSIONS

This approach demonstrates significant computational savings when approaching the full nonlinear inverse problem of optical property estimation using stochastic methods.

Sparse-graph manifold learning method for bioluminescence tomography

In preclinical researches, bioluminescence tomography (BLT) has widely been used for tumor imaging and monitoring, imaged-guided tumor therapy, and so forth. For these biological applications, both tumor spatial location and morphology analysis are the leading problems needed to be taken into account. However, most existing BLT reconstruction methods were proposed for some specific applications with a focus on sparse representation or morphology recovery, respectively. How to design a versatile algorithm that can simultaneously deal with both aspects remains an impending problem. In this study, a Sparse-Graph Manifold Learning (SGML) method was proposed to balance the source sparseness and morphology, by integrating non-convex sparsity constraint and dynamic Laplacian graph model. Furthermore, based on the nonconvex optimization theory and some iterative approximation, we proposed a novel iteratively reweighted soft thresholding algorithm (IRSTA) to solve the SGML model. Numerical simulations and in vivo experiments result demonstrated that the proposed SGML method performed much superior to the comparative methods in spatial localization and tumor morphology recovery for various source settings. It is believed that the SGML method can be applied to the related optical tomography and facilitate the development of optical molecular tomography.

3D Monte Carlo simulation of light distribution in mouse brain in quantitative photoacoustic computed tomography

Background

Photoacoustic computed tomography (PACT) detects light-induced ultrasound (US) waves to reconstruct the optical absorption contrast of the biological tissues. Due to its relatively deep penetration (several centimeters in soft tissue), high spatial resolution, and inherent functional sensitivity, PACT has great potential for imaging mouse brains with endogenous and exogenous contrasts, which is of immense interest to the neuroscience community. However, conventional PACT either assumes homogenous optical fluence within the brain or uses a simplified attenuation model for optical fluence estimation. Both approaches underestimate the complexity of the fluence heterogeneity and can result in poor quantitative imaging accuracy.

Methods

To optimize the quantitative performance of PACT, we explore for the first time 3D Monte Carlo (MC) simulation to study the optical fluence distribution in a complete mouse brain model. We apply the MCX MC simulation package on a digital mouse (Digimouse) brain atlas that has complete anatomy information. To evaluate the impact of the brain vasculature on light delivery, we also incorporate the whole-brain vasculature in the Digimouse atlas. k-wave toolbox was used to investigate the effect of inhomogeneous illumination on the reconstructed images and chromophore concentration estimation.

Results

The simulation results clearly show that the optical fluence in the mouse brain is heterogeneous at the global level and can decrease by a factor of five with increasing depth. Moreover, the strong absorption and scattering of the brain vasculature also induce the fluence disturbance at the local level.

Conclusions

Both global and local fluence heterogeneity contributes to the reduced quantitative accuracy of the reconstructed PACT images of mouse brain. Correcting the optical fluence distribution can improve the quantitative accuracy of PACT.

Modeling voxel-based Monte Carlo light transport with curved and oblique boundary surfaces

SIGNIFICANCE

Monte Carlo (MC) light transport simulations are most often performed in regularly spaced three-dimensional voxels, a type of data representation that naturally struggles to represent boundary surfaces with curvature and oblique angles. Not accounting properly for such boundaries with an index of refractivity, mismatches can lead to important inaccuracies, not only in the calculated angles of reflection and transmission but also in the amount of light that transmits through or reflects from these mismatched boundary surfaces.

AIM

A new MC light transport algorithm is introduced to deal with curvature and oblique angles of incidence when simulated photons encounter mismatched boundary surfaces.

APPROACH

The core of the proposed algorithm applies the efficient preprocessing step of calculating a gradient map of the mismatched boundaries, a smoothing step on this calculated 3D vector field to remove surface roughness due to discretization and an interpolation scheme to improve the handling of curvature.

RESULTS

Through simulations of light hitting the side of a sphere and going through a lens, the agreement of this approach with analytical solutions is shown to be strong.

CONCLUSIONS

The MC method introduced here has the advantage of requiring only slight implementation changes from the current state-of-the-art to accurately simulate mismatched boundaries and readily exploit the acceleration of general-purpose graphics processing units. A code implementation, mcxyzn, is made available and maintained at https://omlc.org/software/mc/mcxyzn/.

Improved Block Sparse Bayesian Learning Method Using K-Nearest Neighbor Strategy for Accurate Tumor Morphology Reconstruction in Bioluminescence Tomography

OBJECTIVE

Bioluminescence tomography (BLT) is a non-invasive technique designed to enable three-dimensional (3D) visualization and quantification of viable tumor cells in living organisms. However, despite the excellent sensitivity and specificity of bioluminescence imaging (BLI), BLT is limited by the photon scattering effect and ill-posed inverse problem. If the complete structural information of a light source is considered when solving the inverse problem, reconstruction accuracy will be improved.

METHODS

This article proposed a block sparse Bayesian learning method based on K-nearest neighbor strategy (KNN-BSBL), which incorporated several types of a priori information including sparsity, spatial correlations among neighboring points, and anatomical information to balance over-sparsity and morphology preservation in BLT. Furthermore, we considered the Gaussian weighted distance prior in a light source and proposed a KNN-GBSBL method to further improve the performance of KNN-BSBL.

RESULTS

The results of numerical simulations and in vivo glioma-bearing mouse experiments demonstrated that KNN-BSBL and KNN-GBSBL achieved superior accuracy for tumor spatial positioning and morphology reconstruction.

CONCLUSION

The proposed method KNN-BSBL incorporated several types of a priori information is an efficient and robust reconstruction method for BLT.

Graphics processing unit-accelerated mesh-based Monte Carlo photon transport simulations

The mesh-based Monte Carlo (MMC) algorithm is increasingly used as the gold-standard for developing new biophotonics modeling techniques in 3-D complex tissues, including both diffusion-based and various Monte Carlo (MC)-based methods. Compared to multilayered and voxel-based MCs, MMC can utilize tetrahedral meshes to gain improved anatomical accuracy but also results in higher computational and memory demands. Previous attempts of accelerating MMC using graphics processing units (GPUs) have yielded limited performance improvement and are not publicly available. We report a highly efficient MMC-MMCL-using the OpenCL heterogeneous computing framework and demonstrate a speedup ratio up to 420× compared to state-of-the-art single-threaded CPU simulations. The MMCL simulator supports almost all advanced features found in our widely disseminated MMC software, such as support for a dozen of complex source forms, wide-field detectors, boundary reflection, photon replay, and storing a rich set of detected photon information. Furthermore, this tool supports a wide range of GPUs/CPUs across vendors and is freely available with full source codes and benchmark suites at http://mcx.space/#mmc.

FullMonteCUDA: a fast, flexible, and accurate GPU-accelerated Monte Carlo simulator for light propagation in turbid media

Optimizing light delivery for photodynamic therapy, quantifying tissue optical properties or reconstructing 3D distributions of sources in bioluminescence imaging and absorbers in diffuse optical imaging all involve solving an inverse problem. This can require thousands of forward light propagation simulations to determine the parameters to optimize treatment, image tissue or quantify tissue optical properties, which is time-consuming and computationally expensive. Addressing this problem requires a light propagation simulator that produces results quickly given modelling parameters. In previous work, we developed FullMonteSW: currently the fastest, tetrahedral-mesh, Monte Carlo light propagation simulator written in software. Additional software optimizations showed diminishing performance improvements, so we investigated hardware acceleration methods. This work focuses on FullMonteCUDA: a GPU-accelerated version of FullMonteSW which targets NVIDIA GPUs. FullMonteCUDA has been validated across several benchmark models and, through various GPU-specific optimizations, achieves a 288-936x speedup over the single-threaded, non-vectorized version of FullMonteSW and a 4-13x speedup over the highly optimized, hand-vectorized and multi-threaded version. The increase in performance allows inverse problems to be solved more efficiently and effectively.

Dual-grid mesh-based Monte Carlo algorithm for efficient photon transport simulations in complex three-dimensional media

The mesh-based Monte Carlo (MMC) method is an efficient algorithm to model light propagation inside tissues with complex boundaries, but choosing appropriate mesh density can be challenging. A fine mesh improves the spatial resolution of the output but requires more computation. We propose an improved MMC-dual-grid mesh-based Monte Carlo (DMMC)-to accelerate photon simulations using a coarsely tessellated tetrahedral mesh for ray-tracing computation and an independent voxelated grid for output data storage. The decoupling between ray-tracing and data storage grids allows us to simultaneously achieve faster simulations and improved output spatial accuracy. Furthermore, we developed an optimized ray-tracing technique to eliminate unnecessary ray-tetrahedron intersection tests in optically thick mesh elements. We validate the proposed algorithms using a complex heterogeneous domain and compare the solutions with those from MMC and voxel-based Monte Carlo. We found that DMMC with an unrefined constrained Delaunay tessellation of the boundary nodes yielded the highest speedup, ranging from 1.3  ×   to 2.9  ×   for various scattering settings, with nearly no loss in accuracy. In addition, the optimized ray-tracing technique offers excellent acceleration in high-scattering media, reducing the ray-tetrahedron test count by over 100-fold. Our DMMC software can be downloaded at http://mcx.space/mmc.

Graphics processing units-accelerated adaptive nonlocal means filter for denoising three-dimensional Monte Carlo photon transport simulations

The Monte Carlo (MC) method is widely recognized as the gold standard for modeling light propagation inside turbid media. Due to the stochastic nature of this method, MC simulations suffer from inherent stochastic noise. Launching large numbers of photons can reduce noise but results in significantly greater computation times, even with graphics processing units (GPU)-based acceleration. We develop a GPU-accelerated adaptive nonlocal means (ANLM) filter to denoise MC simulation outputs. This filter can effectively suppress the spatially varying stochastic noise present in low-photon MC simulations and improve the image signal-to-noise ratio (SNR) by over 5 dB. This is equivalent to the SNR improvement of running nearly 3.5  ×   more photons. We validate this denoising approach using both homogeneous and heterogeneous domains at various photon counts. The ability to preserve rapid optical fluence changes is also demonstrated using domains with inclusions. We demonstrate that this GPU-ANLM filter can shorten simulation runtimes in most photon counts and domain settings even combined with our highly accelerated GPU MC simulations. We also compare this GPU-ANLM filter with the CPU version and report a threefold to fourfold speedup. The developed GPU-ANLM filter not only can enhance three-dimensional MC photon simulation results but also be a valuable tool for noise reduction in other volumetric images such as MRI and CT scans.

Perspective: Prospects of non-invasive sensing of the human brain with diffuse optical imaging

Since the initial demonstration of near-infrared spectroscopy (NIRS) for noninvasive measurements of brain perfusion and metabolism in the 1970s, and its application to functional brain studies (fNIRS) in the 1990s, the field of noninvasive optical studies of the brain has been continuously growing. Technological developments, data analysis advances, and novel areas of application keep advancing the field. In this article, we provide a view of the state of the field of cerebral NIRS, starting with a brief historical introduction and a description of the information content of the NIRS signal. We argue that NIRS and fNIRS studies should always report data of both oxy- and deoxyhemoglobin concentrations in brain tissue, as they complement each other to provide more complete functional and physiological information, and may help identify different types of confounds. One significant challenge is the assessment of absolute tissue properties, be them optical or physiological, so that relative measurements account for the vast majority of NIRS and fNIRS applications. However, even relative measurements of hemodynamics or metabolic changes face the major problem of a potential contamination from extracerebral tissue layers. Accounting for extracerebral contributions to fNIRS signals is one of the most critical barriers in the field. We present some of the approaches that were proposed to tackle this challenge in the study of cerebral hemodynamics and functional connectivity. Finally, we critically compare fNIRS and functional magnetic resonance imaging (fMRI) by relating their measurements in terms of signal and noise, and by commenting on their complementarity.

Optical and thermal simulations for the design of optodes for minimally invasive optogenetics stimulation or photomodulation of deep and large cortical areas in non-human primate brain

The use of optogenetics or photobiomodulation in non-human primate (NHP) requires the ability to noninvasively stimulate large and deep cortical brain tissues volumes. In this context, the optical and geometrical parameters of optodes are critical. Methods and general guidelines to optimize these parameters have to be defined.

OBJECTIVE

We propose the design of an optode for safe and efficient optical stimulation of a large volume of NHP cortex, down to 3-5 mm depths without inserting fibers into the cortex.

APPROACH

Monte Carlo simulations of optical and thermal transport have been carried out using the Geant4 application for tomographic emission (GATE) platform. Parameters such as the fiber diameter, numerical aperture, number of fibers and their geometrical arrangement have been studied. Optimal hardware parameters are proposed to obtain homogeneous fluence above the fluence threshold for opsin activation without detrimental thermal effects.

MAIN RESULTS

The simulations show that a large fiber diameter and a large numerical aperture are preferable since they allow limiting power concentration and hence the resulting thermal increases at the brain surface. To obtain a volume of 200-500 mm³ of brain tissues receiving a fluence above the opsin activation threshold for optogenetics or below a phototocixity threshold for photobiomodulation, a 4 fibers configuration is proposed. The optimal distance between the fibers was found to be 4 mm. A practical implementation of the optode has been performed and the corresponding fluence and thermal maps have been simulated.

SIGNIFICANCE

The present study defines a method to optimize the design of optode and the choice of stimulation parameters for optogenetics and more generally light delivery to deep and large volumes of tissues in NHP brain with a controlled irradiance dosimetry. The general guidelines are the use of silica fibers with a large numerical aperture and a large diameter. The combination of several fibers is required if large volumes need to be stimulated while avoiding thermal effects.

Recent methodology advances in fluorescence molecular tomography

Molecular imaging (MI) is a novel imaging discipline that has been continuously developed in recent years. It combines biochemistry, multimodal imaging, biomathematics, bioinformatics, cell & molecular physiology, biophysics, and pharmacology, and it provides a new technology platform for the early diagnosis and quantitative analysis of diseases, treatment monitoring and evaluation, and the development of comprehensive physiology. Fluorescence Molecular Tomography (FMT) is a type of optical imaging modality in MI that captures the three-dimensional distribution of fluorescence within a biological tissue generated by a specific molecule of fluorescent material within a biological tissue. Compared with other optical molecular imaging methods, FMT has the characteristics of high sensitivity, low cost, and safety and reliability. It has become the research frontier and research hotspot of optical molecular imaging technology. This paper took an overview of the recent methodology advances in FMT, mainly focused on the photon propagation model of FMT based on the radiative transfer equation (RTE), and the reconstruction problem solution consist of forward problem and inverse problem. We introduce the detailed technologies utilized in reconstruction of FMT. Finally, the challenges in FMT were discussed. This survey aims at summarizing current research hotspots in methodology of FMT, from which future research may benefit.

Parallelized Monte Carlo software to efficiently simulate the light propagation in arbitrarily shaped objects and aligned scattering media

A GPU-based Monte Carlo software (MCtet) was developed to calculate the light propagation in arbitrarily shaped objects, like a human tooth, represented by a tetrahedral mesh. A unique feature of MCtet is a concept to realize different kinds of light-sources illuminating the complex-shaped surface of an object, for which no preprocessing step is needed. With this concept, it is also possible to consider photons leaving a turbid media and reentering again in case of a concave object. The correct implementation was shown by comparison with five other Monte Carlo software packages. A hundredfold acceleration compared with central processing units-based programs was found. MCtet can simulate anisotropic light propagation, e.g., by accounting for scattering at cylindrical structures. The important influence of the anisotropic light propagation, caused, e.g., by the tubules in human dentin, is shown for the transmission spectrum through a tooth. It was found that the sensitivity to a change in the oxygen saturation inside the pulp for transmission spectra is much larger if the tubules are considered. Another "light guiding" effect based on a combination of a low scattering and a high refractive index in enamel is described.

Scalable and massively parallel Monte Carlo photon transport simulations for heterogeneous computing platforms

We present a highly scalable Monte Carlo (MC) three-dimensional photon transport simulation platform designed for heterogeneous computing systems. Through the development of a massively parallel MC algorithm using the Open Computing Language framework, this research extends our existing graphics processing unit (GPU)-accelerated MC technique to a highly scalable vendor-independent heterogeneous computing environment, achieving significantly improved performance and software portability. A number of parallel computing techniques are investigated to achieve portable performance over a wide range of computing hardware. Furthermore, multiple thread-level and device-level load-balancing strategies are developed to obtain efficient simulations using multiple central processing units and GPUs.

Toward real-time diffuse optical tomography: accelerating light propagation modeling employing parallel computing on GPU and CPU

Parameter recovery in diffuse optical tomography is a computationally expensive algorithm, especially when used for large and complex volumes, as in the case of human brain functional imaging. The modeling of light propagation, also known as the forward problem, is the computational bottleneck of the recovery algorithm, whereby the lack of a real-time solution is impeding practical and clinical applications. The objective of this work is the acceleration of the forward model, within a diffusion approximation-based finite-element modeling framework, employing parallelization to expedite the calculation of light propagation in realistic adult head models. The proposed methodology is applicable for modeling both continuous wave and frequency-domain systems with the results demonstrating a 10-fold speed increase when GPU architectures are available, while maintaining high accuracy. It is shown that, for a very high-resolution finite-element model of the adult human head with ∼600,000 nodes, consisting of heterogeneous layers, light propagation can be calculated at ∼0.25  s/excitation source.

Advances in Monte Carlo Simulation for Light Propagation in Tissue

Monte Carlo (MC) simulation for light propagation in tissue is the gold standard for studying the light propagation in biological tissue and has been used for years. Interaction of photons with a medium is simulated based on its optical properties. New simulation geometries, tissue-light interaction methods, and recording techniques recently have been designed. Applications, such as whole mouse body simulations for fluorescence imaging, eye modeling for blood vessel imaging, skin modeling for terahertz imaging, and human head modeling for sinus imaging, have emerged. Here, we review the technical advances and recent applications of MC simulation.

Finite element simulation of light transfer in turbid media under structured illumination

The spatial-frequency domain (SFD) imaging technique allows us to estimate the optical properties of biological tissues in a wide field of view. The technique is, however, prone to error in measurement because the two crucial assumptions used for deriving the analytical solution to the diffusion approximation cannot be met perfectly in practical applications. This research mainly focused on modeling light transfer in turbid media under the normal incidence of structured illumination using the finite element method (FEM). Finite element simulations were performed for 50 simulation samples with different combinations of optical absorption and scattering coefficients for varying spatial frequencies, and the results were then compared with the analytical method and Monte Carlo simulation. Relationships between diffuse reflectance and dimensionless absorption and dimensionless scattering coefficients were investigated. The results indicated that the FEM provided reasonable results for diffuse reflectance, compared with the analytical method. Both the FEM and the analytical method overestimated the reflectance for μ_tr/f_x values of greater than 2 and underestimated the reflectance for μ_tr/f_x values of smaller than 2. Larger values of μs'/μ_a yielded better diffuse reflectance estimations than did those of smaller than 10. The reflectance increased nonlinearly with the dimensionless scattering, whereas the reflectance decreased linearly with the dimensionless absorption. It was also observed that diffuse reflectance was relatively stable and insensitive to μs' when the dimensionless scattering was larger than 50. Overall results demonstrate that the FEM is effective for modeling light transfer in turbid media and can be used to explore the effects of crucial parameters for the SFD imaging technique.

Fluorescence Diffusion in the Presence of Optically Clear Tissues in a Mouse Head Model

Diffuse Optical Tomography commonly neglects or assumes as insignificant the presence of optically clear regions in biological tissues, estimating their contribution as a small perturbation to light transport. The inaccuracy introduced by this practice is examined in detail in the context of a complete, based on realistic geometry, virtual fluorescence Diffuse Optical Tomography experiment where a mouse head is imaged in the presence of cerebral spinal fluid. Despite the small thickness of such layer, we point out that an error is introduced when neglecting it from the model with possibly reduction in the accuracy of the reconstruction and localization of the fluorescence distribution within the brain. The results obtained in the extensive study presented here suggest that fluorescence diffuse neuroimaging studies can be improved in terms of quantitative and qualitative reconstruction by accurately taking into account optically transparent regions especially in the cases where the reconstruction is aided by the prior knowledge of the structural geometry of the specimen. Thus, this has only recently become an affordable choice, thanks to novel computation paradigms that allow to run Monte Carlo photon propagation on a simple graphic card, hence speeding up the process a thousand folds compared to CPU-based solutions.

Experimental assessment of the safety and potential efficacy of high irradiance photostimulation of brain tissues

Optogenetics is widely used in fundamental neuroscience. Its potential clinical translation for brain neuromodulation requires a careful assessment of the safety and efficacy of repeated, sustained optical stimulation of large volumes of brain tissues. This study was performed in rats and not in non-human primates for ethical reasons. We studied the spatial distribution of light, potential damage, and non-physiological effects in vivo, in anesthetized rat brains, on large brain volumes, following repeated high irradiance photo-stimulation. We generated 2D irradiance and temperature increase surface maps based on recordings taken during optical stimulation using irradiance and temporal parameters representative of common optogenetics experiments. Irradiances of 100 to 600 mW/mm2 with 5 ms pulses at 20, 40, and 60 Hz were applied during 90 s. In vivo electrophysiological recordings and post-mortem histological analyses showed that high power light stimulation had no obvious phototoxic effects and did not trigger non-physiological functional activation. This study demonstrates the ability to illuminate cortical layers to a depth of several millimeters using pulsed red light without detrimental thermal damages.

GPU acceleration of Monte Carlo simulations for polarized photon scattering in anisotropic turbid media

In earlier studies, we developed scattering models and the corresponding CPU-based Monte Carlo simulation programs to study the behavior of polarized photons as they propagate through complex biological tissues. Studying the simulation results in high degrees of freedom that created a demand for massive simulation tasks. In this paper, we report a parallel implementation of the simulation program based on the compute unified device architecture running on a graphics processing unit (GPU). Different schemes for sphere-only simulations and sphere-cylinder mixture simulations were developed. Diverse optimizing methods were employed to achieve the best acceleration. The final-version GPU program is hundreds of times faster than the CPU version. Dependence of the performance on input parameters and precision were also studied. It is shown that using single precision in the GPU simulations results in very limited losses in accuracy. Consumer-level graphics cards, even those in laptop computers, are more cost-effective than scientific graphics cards for single-precision computation.