In situ estimation of optical properties of rat and monkey brains using femtosecond time-resolved measurements

Influence of Tumor Volume on the Fluence Rate Within Human Breast Model Using Continuous-Wave Diffuse Optical Imaging: A Simulation Study

Objective: This article investigates the effect of varying breast tumor size on the fluence rate distribution within a breast model during the diffuse optical imaging procedure. Background: Early detection of breast cancer is of significant importance owing to its wide spread among women worldwide. Mastectomy surgery became very common due to the late detection of breast cancers by the conventional diagnostic methods such as X-ray mammography and magnetic resonance imaging. On the contrary, optical imaging techniques provide a safe and more sensitive methodology, which is suitable for the early detection criteria. Methods: The implementation was performed based on simulating multiple detectors placed on the outer surface of a human breast model to compute the optical fluence rate after probing the breast (normal and different tumor sizes) with laser irradiation. Different laser wavelengths ranging from the red to near-infrared rays spectral range were examined to determine the optimum fluence rate that shows the highest capability to differentiate between normal and cancerous breasts. A three-dimensional breast model was created using the COMSOL multiphysics package where the optical fluence rate was estimated based on the finite-element solution of the diffusion equation. Results: To evaluate the efficiency of the suggested technique for identifying cancers and discriminate them from normal breast at various wavelengths (600-1000 nm) and several tumor sizes. Conclusions: The obtained results reveal different fluence rate distributions in the breast with different radius tumors, especially at 600 nm due to the significant differences in the scattering coefficient between malignancies and healthy tissue.

Optical Characterization of Biological Tissues Based on Fluorescence, Absorption, and Scattering Properties

Optical diagnostics methods are significantly appealing in biological applications since they are non-destructive, safe, and minimally invasive. Laser-induced fluorescence is a promising optical spectrochemical analytical technique widely employed for tissue classification through molecular analysis of the studied samples after excitation with appropriate short-wavelength laser light. On the other hand, diffuse optics techniques are used for tissue monitoring and differentiation based on their absorption and scattering characteristics in the red to the near-infrared spectra. Therefore, it is strongly foreseen to obtain promising results by combining these techniques. In the present work, tissues under different conditions (hydrated/dry skin and native/boiled adipose fat) were distinguished according to their fluorescence emission, absorption, and scattering properties. The selected tissues' optical absorption and scattering parameters were determined via Kubelka-Munk mathematical model according to the experimental tissue reflectance and transmittance measurements. Such measurements were obtained using an optical configuration of integrating sphere and spectrometer at different laser wavelengths (808, 830, and 980 nm). Moreover, the diffusion equation was solved for the fluence rate at the sample surface using the finite element method. Furthermore, the accuracy of the obtained spectroscopic measurements was evaluated using partial least squares regression statistical analysis with 0.87 and 0.89 R-squared values for skin and adipose fat, respectively.

Experimental imaging and Monte Carlo modeling of ultrafast pulse propagation in thin scattering slabs

SIGNIFICANCE

Most radiative transport problems in turbid media are typically associated with mm or cm scales, leading to typical time scales in the range of hundreds of ps or more. In certain cases, however, much thinner layers can also be relevant, which can dramatically alter the overall transport properties of a scattering medium. Studying scattering in these thin layers requires ultrafast detection techniques and adaptations to the common Monte Carlo (MC) approach.

AIM

We aim to discuss a few relevant aspects for the simulation of light transport in thin scattering membranes, and compare the obtained numerical results with experimental measurements based on an all-optical gating technique.

APPROACH

A thin membrane with controlled scattering properties based on polymer-dispersed TiO2 nanoparticles is fabricated for experimental validation. Transmittance measurements are compared against a custom open-source MC implementation including specific pulse profiles for tightly focused femtosecond laser pulses.

RESULTS

Experimental transmittance data of ultrafast pulses through a thin scattering sample are compared with MC simulations in the spatiotemporal domain to retrieve its scattering properties. The results show good agreement also at short distances and time scales.

CONCLUSIONS

When simulating light transport in scattering membranes with thicknesses in the orders of tens of micrometer, care has to be taken when describing the temporal, spatial, and divergence profiles of the source term, as well as the possible truncation of step length distributions, which could be introduced by simple strategies for the generation of random exponentially distributed variables.

Silent microscopy to explore a brain that hears butterflies' wings

A silent two-photon laser-scanning microscopy system, which eliminates mechanical vibrations in the audible range, has enabled the detection of auditory cortical neurons with responses at sound pressure levels as low as 5 dB in nonhuman primates.

Model equations of light scattering properties and a characteristic time of light propagation for polydisperse colloidal suspensions at different volume fractions

We developed model equations of light scattering properties and a characteristic time of light propagation for polydisperse colloidal suspensions at different volume fractions. By the model equations, we examined numerical results using the first-order (dependent) scattering theory (FST) and radiative transfer theory in 600-980 nm wavelength. The model equations efficiently treat the interference of electric fields scattered from colloidal particles by a single effective coefficient, providing fast computation. Meanwhile, the FST provides accurate but complicated treatment. We found the interference effects on the scattering properties and characteristic time depend linearly on wavelength. Dimensionless analysis showed a simple mechanism of the interference effects, independently of wavelength and source-detector distance.

Numerical Study of Near-Infrared Light Propagation in Aqueous Alumina Suspensions Using the Steady-State Radiative Transfer Equation and Dependent Scattering Theory

Understanding light propagation in liquid phantoms, such as colloidal suspensions, involves fundamental research of near-infrared optical imaging and spectroscopy for biological tissues. Our objective is to numerically investigate light propagation in the alumina colloidal suspensions with the mean alumina particle diameter of 55 nm at the volume fraction range 1–20%. We calculated the light scattering properties using the dependent scattering theory (DST) on a length scale comparable to the optical wavelength. We calculated the steady-state radiative transfer and photon diffusion equations (RTE and PDE) using the DST results based on the finite difference method in a length scale of the mean free path. The DST calculations showed that the scattering and reduced scattering coefficients become more prominent at a higher volume fraction. The anisotropy factor is almost zero at all the volume fractions, meaning the scattering is isotropic. The comparative study of the RTE with the PDE showed that the diffusion approximation holds at the internal region with all the volume fractions and the boundary region with the volume fraction higher than 1%. Our findings suggest the usefulness of the PDE as a light propagation model for the alumina suspensions rather than the RTE, which provides accurate but complicated computation.

Photon transport model for dense polydisperse colloidal suspensions using the radiative transfer equation combined with the dependent scattering theory

In near-infrared imaging and spectroscopy, high-fidelity modeling of photon transport for dense polydisperse colloidal suspensions is crucial. We developed photon transport models using the radiative transfer equation (RTE) with the dependent scattering theory (DST) at volume fractions up to 20%. The polydispersity and interference effects strongly influence results of the scattering properties and the RTE in cases of small mean diameter and large variance of the particle size distribution. We compared the RTE-results for the Henyey-Greenstein (conventional) function with those for the phase function using the DST. The RTE-results differ between both functions at low volume fractions for forward scattering media, suggesting the limitation of the conventional function.

Investigating the transmission profiles of 808 nm laser through different regions of the rat's head

Studying light penetration in biological tissues became a very important concern in various medical applications. It is an essential factor required to resolve the optical dose in many diagnostic and therapeutic procedures. The absorption and scattering properties of the inspected tissue control how deep the light will travel inside the tissue. However, these optical properties are highly dependent on the wavelength of the light source. In this work, the light transmission through different regions of the rat's head was investigated and the minimum laser power required to reach different parts of the head is also determined using 808-nm semiconductor laser diode. The power variation in different regions of the head is estimated using Monte Carlo simulation. Absorption and scattering coefficients of the head layers were calculated using integrating sphere measurements and Kubelka-Munk model. The absorption coefficient of the skin was 0.19 ± 0.071 mm^-1, 0.024 ± 0.11 mm^-1 for skull, and 0.35 ± 0.13 mm^-1 for the brain, while the scattering coefficients were 7.35 ± 1.09, 2.71 ± 0.37, and 13.04 ± 0.36 mm^-1 for skin, skull, and brain, respectively. The obtained results provide a relationship between laser incident power and the depth in the rat's head showing a higher optical transmission at the frontal part of the head than the middle or back regions due to the variations in the skull thickness. Therefore, the study revealed that the transmitted power of 808 nm laser at different incident locations on the head is nonlinear and variable due to different skull's thickness.

Short-wave Infrared Neural Stimulation Drives Graded Sciatic Nerve Activation Across A Continuum of Wavelengths

Infrared neural stimulation (INS) is an optical stimulation technique which uses coherent light to stimulate nerves and neurons and which shows increased spatial selectivity compared to electrical stimulation. This could improve deep brain, high channel count, or vagus nerve stimulation. In this study, we seek to understand the wavelength dependence of INS in the near-infrared optical window. Rat sciatic nerves were excised ex vivo and stimulated with wavelengths between 700 and 900 nm. Recorded compound nerve action potentials (CNAPs) showed that stimulation was maximized in the 700 nm window despite comparable laser power levels across wavelengths. Computational models demonstrated that wavelength-based activation dependencies were not a result of passive optical properties. This data demonstrates that INS is both wavelength and power level dependent, which inform stimulation systems to actively target neural microcircuits in humans.

Restored interlaced volumetric imaging increases image quality and scanning speed during intravital imaging in living mice

Dynamic intravital imaging is essential for revealing ongoing biological phenomena within living organisms and is influenced primarily by several factors: motion artifacts, optical properties and spatial resolution. Conventional imaging quality within a volume, however, is degraded by involuntary movements and trades off between the imaged volume, imaging speed and quality. To balance such trade-offs incurred by two-photon excitation microscopy during intravital imaging, we developed a unique combination of interlaced scanning and a simple image restoration algorithm based on biological signal sparsity and a graph Laplacian matrix. This method increases the scanning speed by a factor of four for a field size of 212 μm × 106 μm × 130 μm, and significantly improves the quality of four-dimensional dynamic volumetric data by preventing irregular artifacts due to the movement observed with conventional methods. Our data suggest this method is robust enough to be applied to multiple types of soft tissue.

Facilitated glass separation by asymmetric Bessel-like beams

The modification of transparent materials is enabled by focused ultrashort laser pulses. Single pass processing up to several millimeters can be achieved by the usage of elongated beam profiles. We studied the mechanical separability in dependence of the material thickness. As simulations show, asymmetric beam profiles can cause modifications with preferential direction reducing the necessary breaking force. Pump-probe microscopy is implemented to examine the laser-matter-interaction. We present a measured 3D-reconstruction of the transient interaction inside the material and elucidate the desired crack formation. We demonstrate beam shaping concepts to create a new, efficient and robust class of Bessel-like beams, which can be used to achieve a preferred crack direction. We verify the concept by modification and separation of silicate glasses.

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.