Monte Carlo simulation of light transport in turbid medium with embedded object-spherical, cylindrical, ellipsoidal, or cuboidal objects embedded within multilayered tissues

Monte Carlo modeling of light transport in multilayered tissue (MCML) is modified to incorporate objects of various shapes (sphere, ellipsoid, cylinder, or cuboid) with a refractive-index mismatched boundary. These geometries would be useful for modeling lymph nodes, tumors, blood vessels, capillaries, bones, the head, and other body parts. Mesh-based Monte Carlo (MMC) has also been used to compare the results from the MCML with embedded objects (MCML-EO). Our simulation assumes a realistic tissue model and can also handle the transmission/reflection at the object-tissue boundary due to the mismatch of the refractive index. Simulation of MCML-EO takes a few seconds, whereas MMC takes nearly an hour for the same geometry and optical properties. Contour plots of fluence distribution from MCML-EO and MMC correlate well. This study assists one to decide on the tool to use for modeling light propagation in biological tissue with objects of regular shapes embedded in it. For irregular inhomogeneity in the model (tissue), MMC has to be used. If the embedded objects (inhomogeneity) are of regular geometry (shapes), then MCML-EO is a better option, as simulations like Raman scattering, fluorescent imaging, and optical coherence tomography are currently possible only with MCML. (C) 2014 Society of Photo-Optical Instrumentation Engineers (SPIE)

Approximation of Internal Refractive Index Variation Improves Image Guided Diffuse Optical Tomography of Breast

Effective usage of image guidance by incorporating the refractive index (RI) variation in computational modeling of light propagation in tissue is investigated to assess its impact on optical-property estimation. With the aid of realistic patient breast three-dimensional models, the variation in RI for different regions of tissue under investigation is shown to influence the estimation of optical properties in image-guided diffuse optical tomography (IG-DOT) using numerical simulations. It is also shown that by assuming identical RI for all regions of tissue would lead to erroneous estimation of optical properties. The a priori knowledge of the RI for the segmented regions of tissue in IG-DOT, which is difficult to obtain for the in vivo cases, leads to more accurate estimates of optical properties. Even inclusion of approximated RI values, obtained from the literature, for the regions of tissue resulted in better estimates of optical properties, with values comparable to that of having the correct knowledge of RI for different regions of tissue.

Mesh Simplification Based on Edge Collapsing Could Improve Computational Efficiency in Near Infrared Optical Tomographic Imaging

The diffusion equation-based modeling of near infrared light propagation in tissue is achieved by using finite-element mesh for imaging real-tissue types, such as breast and brain. The finite-element mesh size (number of nodes) dictates the parameter space in the optical tomographic imaging. Most commonly used finite-element meshing algorithms do not provide the flexibility of distinct nodal spacing in different regions of imaging domain to take the sensitivity of the problem into consideration. This study aims to present a computationally efficient mesh simplification method that can be used as a preprocessing step to iterative image reconstruction, where the finite-element mesh is simplified by using an edge collapsing algorithm to reduce the parameter space at regions where the sensitivity of the problem is relatively low. It is shown, using simulations and experimental phantom data for simple meshes/domains, that a significant reduction in parameter space could be achieved without compromising on the reconstructed image quality. The maximum errors observed by using the simplified meshes were less than 0.27% in the forward problem and 5% for inverse problem.

Analytical solutions for diffuse fluorescence spectroscopy/imaging in biological tissues. Part I: zero and extrapolated boundary conditions

The mathematical model for diffuse fluorescence spectroscopy/imaging is represented by coupled partial differential equations (PDEs), which describe the excitation and emission light propagation in soft biological tissues. The generic closed-form solutions for these coupled PDEs are derived in this work for the case of regular geometries using the Green's function approach using both zero and extrapolated boundary conditions. The specific solutions along with the typical data types, such as integrated intensity and the mean time of flight, for various regular geometries were also derived for both time-and frequency-domain cases. (C) 2013 Optical Society of America

Effect of ion temperature on propagation of ion-acoustic solitary waves of small amplitudes in collisionless plasma

By using a perturbation technique, the Korteweg-de Vries equation is derived for a mixture of warm-ion fluid and hot, isothermal electrons. Stationary solutions are obtained for this equation and are compared with the corresponding solutions for a mixture consisting of cold-ion fluid and hot, isothermal electrons.

Propagation of ion-acoustic waves in a multi-component plasma

The authors derive the Korteweg-de Vries equation in a multicomponent plasma that includes any number of positive and negative ions. The solitary wave solutions are also found explicitly for the case of isothermal and non-isothermal electrons.

Effect of ionic temperature on propagation of ion-acoustic solitary waves in a collisionless plasma of warm-ion fluid and non-isothermal electrons

Using a perturbation technique, we derive Modified Korteweg—de Vries (MKdV) equations for a mixture of warm-ion fluid (γ i = 3) and hot and non-isothermal electrons (γ e> 1), (i) when deviations from isothermality are finite, and (ii) when deviations from isothermality are small. We obtain stationary solutions for these equations, and compare them with the corresponding solutions for a mixture of warm-ion fluid (γ i = 3) and hot, isothermal electrons (γ i = 1).

Alfven surface waves along a cylindrical plasma column surrounded by neutral gas

The dispersion characteristics of Alfven surface waves along a cylindrical plasma column insulated by a neutral gas are discussed. There is no qualitative change in the characteristic curves below the critical magnetic field, given by vA approximately=s, as compared to the propagation of surface waves along the plasma-plasma interface. For magnetic fields above this critical value, there exists a cut-off wave number kc, which depends upon the azimuthal wave number, the radius of the cylinder, the strength of the magnetic field above the critical value and the gas pressure, such that surface waves do not exist for kplasma column.

Study of turbid media with light: Recovery of mechanical and optical properties from boundary measurement of intensity autocorrelation of light

We discuss the inverse problem associated with the propagation of the field autocorrelation of light through a highly scattering object like tissue. In the first part of the work, we reconstructed the optical absorption coefficient mu(u) and particle diffusion coefficient D-B from simulated measurements which are integrals of a quantity computed from the measured intensity and intensity autocorrelation g(2)(tau) at the boundary. In the second part we recover the mean square displacement (MSD) distribution of particles in an inhomogeneous object from the sampled g(2)(tau) measure on the boundary. From the MSD, we compute the storage and loss moduli distributions in the object. We have devised computationally easy methods to construct the sensitivity matrices which are used in the iterative reconstruction algorithms for recovering these parameters from the measurements. The results of the reconstruction of mu(a), D-B, MSD and the viscoelastic parameters, which are presented, show reasonable good position and quantitative accuracy.

Low-frequency degenerate surface modes of an electron-hole plasma

The computations of Flahive and Quinn1 of the dispersion curves of low frequency degenerate surface (DS) modes propagating along the magnetic field in an electron-hole plasma are extended to higher values of the wavenumber. We find that beyond a certain value of the wavenumber the DS mode re-enters the allowed region of surface wave propagation and tends to an asymptotic frequency ωR (<ωLH). These low frequency resonances of an electron-hole plasma are discussed with reference to the experimental observations.

Instability of Coupled Longitudinal and Transverse Electromagnetic Modes of Wave‐Propagation in Bounded Streaming Plasmas

The instability of coupled longitudinal and transverse electromagnetic modes associated with long wavelengths is studied in bounded streaming plasmas. The main conclusions are as follows: (i) For long waves for which O (k 2)=0, in the absence of relative streaming motion of electrons and ions and aωp/c<0.66, the whole spectrum of harmonic waves is excited due to finite temperature and boundary effects consisting of two subseries. One of these subseries can be identified with Tonks-Dattner resonance oscillations for the electrons, and arises primarily due to the electrons with frequencies greater than the electrostatic plasma frequency corresponding to the electron density in the midplane in the undisturbed state. The other series arises primarily due to ion motion. When aωp/c>0.66, in addition to the above spectrum of harmonic waves, the system admits an infinite number of growing and decaying waves. The instability associated with these modes is found to arise due to the interaction of the waves inside the plasma with the external electromagnetic field. (ii) For modes with comparatively shorter wavelengths for which O (k3)=0, the coupling due to finite temperature sets in, and it is found that the two series of harmonic waves obtained in (i) deriving energy from the transverse modes also become unstable. Thus, for these wavelengths the system admits three sets of growing and decaying modes, first two for all values of aωp/c and the third for (aωp/c) > 0.66. (iii) The presence of streaming velocities introduces various other coupling mechanisms, and we find that even for the wavelengths for which O (k2)=0, we get three sets of growing and decaying waves. The numerical values for the growth rates show that the streaming velocities enhance the growth rates of instability significantly.

Propagation of an isothermal shock in stellar envelopes

We have studied in this paper the propagation of an isothermal shock in the radiative envelopes of the Bosman-Crespin model for a hot star and Boury’s model for a giant star. A spherically symmetric disturbance is supposed to be originated at or outside the surface of the convective core. We have used Whitham’s rule to study the variation in the shock strength and the shock velocity after modifying it for inclusion of pressure, energy and flux of radiation. We find the shock increases in strength as it propagates through the envelopes of decreasing density, pressure and temperature. The velocity of the shock decreases for very weak initial shock strengths, for intermediate initial shock strength it first decreases and then increases, while for large initial shock strength, it always increases. This aspect of the problem throws some light on the stability of the models under consideration.

Electromagnetic wave propagation at the interface between two polar semiconductors

The interface between two polar semiconductors can support three types of phonon-plasmon-polariton modes propagating in three well-defined frequency windows ??1?[min(?1,?3),?R1], ??2?[max(?2,?4),?R2], and ??3?[min(?2,?4),?R3]. The limiting frequencies ?1,2,3,4 are defined by ?1(?)=0, ?2(?)=0, and ?R1,2,3 by ?1(?)+?2(?)=0, where ?i(?) are dielectric functions of the two media with i=1,2. The dispersion, decay distances, and polarization of the three modes are discussed. The variation of the limiting frequencies with the interface plasma parameter ???p22/?p12 reveals an interesting feature in the dispersion characteristics of these modes. For the interfaces for which the bulk coupled phonon-plasmon frequencies of medium 1 are greater than the LO frequency or are less than the TO frequency of medium 2, there exist two values of ?=?1 and ?2(

The modulation of radiation in an electron-positron plasma

The modulational instability of a large-amplitude, linearly polarized electromagnetic wave propagating in an electron-positron plasma is considered, including the combined effect of relativistic mass variation of the plasma particles, harmonic generation, and the non-resonant, finite-frequency electrostatic density perturbations, all caused by the large-amplitude radiation field. The radiation from many strong sources, such as AGN and pulsars, has been observed to vary over a host of time-scales. It is possible that the extremely rapid variations in the non-thermal continuum of AGN, as well as in the non-thermal radio radiation from pulsars, can be accounted for by the modulational instabilities to which radiation may be subjected during its propagation out of the emission region.

Ultrasound-modulated optical tomography: recovery of amplitude of vibration in the insonified region from boundary measurement of light correlation

We address a certain inverse problem in ultrasound-modulated optical tomography: the recovery of the amplitude of vibration of scatterers [p(r)] in the ultrasound focal volume in a diffusive object from boundary measurement of the modulation depth (M) of the amplitude autocorrelation of light [phi(r, tau)] traversing through it. Since M is dependent on the stiffness of the material, this is the precursor to elasticity imaging. The propagation of phi(r, tau) is described by a diffusion equation from which we have derived a nonlinear perturbation equation connecting p(r) and refractive index modulation [Delta n(r)] in the region of interest to M measured on the boundary. The nonlinear perturbation equation and its approximate linear counterpart are solved for the recovery of p(r). The numerical results reveal regions of different stiffness, proving that the present method recovers p(r) with reasonable quantitative accuracy and spatial resolution. (C) 2011 Optical Society of America