A Monte Carlo model for the detection efficiency of cylindrical geometry used in fission isomer studies

This paper describes a Monte Carlo calculation of the detection efficiency of cylindrical geometry, using solid state nuclear track detectors. A method is presented for determining the half-life of a fission isomer, which gives results in good agreement with direct electronic timing measurements.     

Monte Carlo simulations and laser Doppler flow measurements with high penetration depth in biological tissuelike head phantoms

Laser Doppler flow measurements on biological tissuelike phantoms have shown that penetration depths of 30 mm could be obtained, thus exceeding the penetration depth of commercial instruments for the measurement of skin perfusion by more than an order of magnitude. Monte Carlo simulations were performed and compared with measurement results obtained on a headlike tissue model to quantify the influence of perfusion of the scalp on the cortex perfusion results. We found Doppler frequency spectra to be independent of the mean scattering angle and could be fitted with a sum of Gaussian functions, using a simple analytical model.     

Development of a bioengineered tissue model and its application in the investigation of the depth selectivity of polarization gating

Understanding the propagation of polarized light in tissue is crucial for a number of biomedical optics applications. Here we report the development of a bioengineered connective tissue model fabricated by the combination of scaffolding and cross-linking techniques to study light transport in biological tissue. It demonstrates great similarity to real connective tissue in its optical properties as well as microarchitecture. Moreover, the optical properties of the model can be reproducibly controlled. As an example, we report the utilization of this model to study the effect of epithelium and the underlying connective tissue on the depth selectivity of polarization gating.     

Crossed source-detector geometry for a novel spray diagnostic: Monte Carlo simulation and analytical results

Sprays and other industrially relevant turbid media can be quantitatively characterized by light scattering. However, current optical diagnostic techniques generate errors in the intermediate scattering regime where the average number of light scattering is too great for the single scattering to be assumed, but too few for the diffusion approximation to be applied. Within this transitional single-to-multiple scattering regime, we consider a novel crossed source-detector geometry that allows the intensity of single scattering to be measured separately from the higher scattering orders. We verify Monte Carlo calculations that include the imperfections of the experiment against analytical results. We show quantitatively the influence of the detector numerical aperture and the angle between the source and the detector on the relative intensity of the scattering orders in the intermediate single-to-multiple scattering regime. Monte Carlo and analytical calculations of double light-scattering intensity are made with small particles that exhibit isotropic scattering. The agreement between Monte Carlo and analytical techniques validates use of the Monte Carlo approach in the intermediate scattering regime. Monte Carlo calculations are then performed for typical parameters of sprays and aerosols with anisotropic (Mie) scattering in the intermediate single-to-multiple scattering regime.     

Controlled Monte Carlo method for light propagation in tissue of semi-infinite geometry

The controlled Monte Carlo method is generalized to model photon migration in turbid media of arbitrary geometries. Its implementation for the reflection geometry is exemplified in this paper. The most probable diffuse direction of light is used as the local attractive vector that serves as the basis of biased sampling of scattering angles. Consequently, path-length resolved photon trajectories can be generated with a significantly improved efficiency. We report a more than 29 times reduction in simulation time for early arriving photons in a typical configuration.     

Multiple scattering from Chebyshev particles: Monte Carlo simulations for backscattering in lidar geometry

Lidar measurements are often interpreted on the basis of two fundamental assumptions: absence of multiple scattering and sphericity of the particles that make up the diffusing medium. There are situations in which neither holds true. We focus our interest on multiply-scattered returns from homogeneous layers of monodisperse, randomly oriented, axisymmetric nonspherical particles. T(2) Chebyshev particles have been chosen and their single-scattering properties have been reviewed. A Monte Carlo procedure has been employed to calculate the backscattered signal for several fields of view. Comparisons with the case of scattering from equivalent (equal-volume) spheres have been carried out (narrow polydispersions have been used to smooth the phase functions' oscillations). Our numerical effort highlights a considerable variability in the intensity of the multiply-scattered signal, which is a consequence of the strong dependence of the backscattering cross section on deformation of the particles. Even more striking effects have been noted for depolarization; peculiar behavior was observed at moderate optical depths when particles characterized by a large backscattering depolarization ratio were employed in our simulations. The sensitivity of depolarization to even small departures from sphericity, in spite of random orientation of the particles, has been confirmed. The results obtained with the Monte Carlo codes have been successfully checked with an analytical formula for double scattering.     

Comparison of radiance and polarization values observed in the Mediterranean Sea and simulated in a Monte Carlo model

Measurements of the radiance and degree of polarization made in 1971 in the Mediterranean Sea are presented along with the simulation of all observed quantities by a Monte Carlo technique. It is shown that our independent scattering treatment utilizing a Stokes vector formalism to describe the polarization state of the light field produces remarkably good agreement with those values measured in situ.     

Forward-adjoint fluorescence model: Monte Carlo integration and experimental validation

The adjoint form of the photon transport equation is applied to a generalized fluorescence detection problem, and its accuracy is empirically tested. This approach can be interpreted as mathematically reversing the temporal flow of fluorescent photons; that is, they are tracked from the detector back to potential sites of origin in the scattering medium. The result is a distribution of potential fluorescing sites that, when properly normalized, gives a probability field of the relative importance of the photon starting position and direction to the resulting signal. This adjoint solution can be combined with the temporally forward-derived distribution of absorbed excitation photons to evaluate the fluorescence excitation detection scheme. This bypasses the normal, temporal derivation wherein the fluorescence transport solution is dependent on the result of the excitation transport solution.     

Half-ball lens couples a beveled fiber probe for depth-resolved spectroscopy: Monte Carlo simulations

In this study, we propose a beveled fiber-optic probe coupled with a half-ball lens for improving the depth-resolved fluorescence measurements of epithelial tissue. The Monte Carlo (MC) simulation results show that for a given excitation-collection fiber separation, the probe design with a bevel-angled collection fiber is more sensitive to detect fluorescence photons emitted from the shallow layer of tissue, whereas the flat-tip collection fiber is in favor of probing fluorescence photons originating from deeper tissue areas. This compact half-ball lens-beveled fiber probe design has the potential to facilitate the depth-resolved fluorescence detection of epithelial tissue.     

Transillumination optical tomography of tissue-engineered blood vessels: a Monte Carlo simulation

A Monte Carlo technique has been developed to simulate the transillumination laser computed tomography of tissue-engineered blood vessels. The blood vessel was modeled as a single cylinder layer mounted on a tubular mandrel. Sequences of images were acquired while rotating the mandrel. The tomographic image was reconstructed by applying a standard Radon transform. Angular discrimination was applied to simulate a spatial filter, which was used to reject multiply scattered photons. The simulation results indicated that the scattering effect can be overcome with angular discrimination because of the thin tissue thickness. However, any refractive-index mismatch among the tissue, the surrounding media, and the mandrel could produce significant distortions in the reconstructed image.     

Comparison of simplified Monte Carlo simulation and diffusion approximation for the fluorescence signal from phantoms with typical mouse tissue optical properties

A simplified approach is proposed to simulate the fluorescence signal from a fluorophore submerged inside a turbid medium using the Monte Carlo method. Based on the reversibility of photon propagation, the fluorescence signal can be obtained from a single Monte Carlo simulation of the excitation light. This is computationally less expensive and also allows for the direct use of well-validated nonfluorescence photon migration Monte Carlo codes. Fluorescence signals from a mouse tissuelike phantom were computed using both the simplified Monte Carlo simulation and the diffusion approximation. The relative difference of signal intensity was found to be at most 30% for a fluorophore placed in the medium at various depths and horizontally midway between a source-detector pair separated by 3 mm. The difference in time characteristics of the signal is also examined.     

Markov-chain Monte Carlo approach to the design of multilayer thin-film optical coatings

We reconsider the problem of locating the globally optimal solution of a multilayer-optical-coating design problem, within some predetermined space of parameters, with the aim of obtaining a robust technique that requires a minimum of user intervention. The approach we adopt centers on exploring the space of the parameters of interest by using a Markov-chain Monte Carlo sampling algorithm. This technique enables one to locate the global optimum automatically with high confidence and without the need for a good starting design. It also allows the trivial inclusion of prior constraints on the variables and provides a natural means for investigating the robustness of the optimal solution.     

The influence of sample geometry on the magnetoresistance of Ni-Fe films

The magnetoresistance of Ni-Fe (81–91) films has been measured as a function of film geometry (thickness to width ratio N). It has been found that under the condition of large magnetization dispersion the behaviour of films with a high value of N differs distinctly from that of films with N ≈ 0. Moreover, the easy axis coercive force H_c increases markedly when the width of the film is of the order of 20 μm and less.     

Determination of optical probe interrogation field of near-infrared reflectance: phantom and Monte Carlo study

An optical probe used to localize human brain tissues in vivo has been reported previously. It was able to sense the underlying tissue structure with an optical interrogation field, termed as "look ahead distance" (LAD). A new side-firing probe has been designed with its optical window along its side. We have defined the optical interrogation field of the new side probe as "look aside distance" (LASD). The purpose of this study is to understand the dependence of the LAD and LASD on the optical properties of tissue, the light source intensity, and the integration time of the detector, using experimental and computational methods. The results show that a decrease in light intensity does decrease the LAD and LASD and that an increase in integration time of detection may not necessarily improve the depths of LAD and LASD. Furthermore, Monte Carlo simulation results suggest that the LAD/LASD decreases with an increase in reduced scattering coefficient to a point, after which the LAD/LASD remains constant. We expect that an optical interrogation field of a tip or side probe is approximately 1-2 mm in white matter and 2-3.5 mm in gray matter. These conclusions will help us optimally manipulate the parameter settings during surgery and determine the spatial resolution of the probe.     

Monte Carlo-based inverse model for calculating tissue optical properties. Part I: Theory and validation on synthetic phantoms

A flexible and fast Monte Carlo-based model of diffuse reflectance has been developed for the extraction of the absorption and scattering properties of turbid media, such as human tissues. This method is valid for a wide range of optical properties and is easily adaptable to existing probe geometries, provided a single phantom calibration measurement is made. A condensed Monte Carlo method was used to speed up the forward simulations. This model was validated by use of two sets of liquid-tissue phantoms containing Nigrosin or hemoglobin as absorbers and polystyrene spheres as scatterers. The phantoms had a wide range of absorption (0-20 cm(-1)) and reduced scattering coefficients (7-33 cm(-1)). Mie theory and a spectrophotometer were used to determine the absorption and reduced scattering coefficients of the phantoms. The diffuse reflectance spectra of the phantoms were measured over a wavelength range of 350-850 nm. It was found that optical properties could be extracted from the experimentally measured diffuse reflectance spectra with an average error of 3% or less for phantoms containing hemoglobin and 12% or less for phantoms containing Nigrosin.     

Monte Carlo-based inverse model for calculating tissue optical properties. Part II: Application to breast cancer diagnosis

The Monte Carlo-based inverse model of diffuse reflectance described in part I of this pair of companion papers was applied to the diffuse reflectance spectra of a set of 17 malignant and 24 normal-benign ex vivo human breast tissue samples. This model allows extraction of physically meaningful tissue parameters, which include the concentration of absorbers and the size and density of scatterers present in tissue. It was assumed that intrinsic absorption could be attributed to oxygenated and deoxygenated hemoglobin and beta-carotene, that scattering could be modeled by spheres of a uniform size distribution, and that the refractive indices of the spheres and the surrounding medium are known. The tissue diffuse reflectance spectra were evaluated over a wavelength range of 400-600 nm. The extracted parameters that showed the statistically most significant differences between malignant and nonmalignant breast tissues were hemoglobin saturation and the mean reduced scattering coefficient. Malignant tissues showed decreased hemoglobin saturation and an increased mean reduced scattering coefficient compared with nonmalignant tissues. A support vector machine classification algorithm was then used to classify a sample as malignant or nonmalignant based on these two extracted parameters and produced a cross-validated sensitivity and specificity of 82% and 92%, respectively.     

Monte Carlo model for studying the effects of melanin concentrations on retina light absorption

We developed a Monte Carlo model to calculate light absorption in human and mouse retinas. The retina was modeled as a five-layer spherical structure. The effects of melanin concentrations in the retinal pigment epithelium (RPE) and choroid layer were studied. Variations of blood content in choroid were also considered in the simulation. Our simulation results indicated that light absorption in neural retina was at least 20% higher in albino subjects than in pigmented subjects under both photobleaching and dark-adapted conditions. It can be four times higher at optical wavelengths corresponding to minimal hemoglobin absorption. The elevated absorption at neural retina was attributed to the light backscattered from the choroid and sclera layers. This simulation model may provide useful information in studying light-induced retina damage.     

Twilight polarization and optical depth of stratospheric aerosols over Beijing after the Pinatubo volcanic eruption

After the volcanic eruption of Mt. Pinatubo the degree of polarization of skylight during twilight over Beijing was monitored with a polarimeter aimed at the local zenith. We analyze the effect of changes in the scattering coefficient of atmospheric aerosols for the case of multiple scattering on skylight polarization at the zenith and then discuss the evolution of skylight polarization over Beijing during the posteruption period. As a reference and for comparison we also discuss the evolution of the aerosol optical depth retrieved from the combination of skylight polarization and backscattering ratio measured by the polarimeter and a lidar for the period beginning with the eruption of Mt. Pinatubo through the end of 1993. The contributions of atmospheric aerosols at different altitudes to the ground-observed twilight polarization depend on the solar zenith angle. For larger solar zenith angles, the skylight polarization is mostly sensitive to aerosol variations in the upper layer that range from 15 to 30 km. The twilight polarization at the zenith from June 1991 to mid-1994 shows different features for three periods: (1) From October 1991 to February 1992, volcanic dust traveled to mid-latitudes, and the degree of polarization decreased substantially. (2) From February 1992 to November 1993, volcanic dust was dispersed the minimum degree of polarization at the solar zenith angle of 93.5 degrees disappeared and the maximum increased. In addition, polarization for solar zenith angles less than 90 degrees also increased. (3) From November 1993 to May 1994, most of the volcanic dust had fallen off, the atmosphere was restored to the background state, and the skylight polarization approached the preeruption condition.