• Laser & Optoelectronics Progress
  • Vol. 58, Issue 6, 600003 (2021)
Lü Chenyang and Zhan Renjun*
Author Affiliations
  • College of Equipment Management and Support, Engineering University of People’s Armed Police, Xi''an, Shaanxi 710086, China
  • show less
    DOI: 10.3788/LOP202158.0600003 Cite this Article Set citation alerts
    Lü Chenyang, Zhan Renjun. Theoretical Models of Light Distribution in Biological Tissues Irradiated by Laser[J]. Laser & Optoelectronics Progress, 2021, 58(6): 600003 Copy Citation Text show less
    Dynamic analysis model of photothermal effect of biological tissues
    Fig. 1. Dynamic analysis model of photothermal effect of biological tissues
    Hierarchical structure of photon transport theory
    Fig. 2. Hierarchical structure of photon transport theory
    ParameterPhysical meaningParameterPhysical meaning
    JAttenuated irradiancesDirection before scattering
    r=(r, z)Spatial locationsDirection after scattering
    dω'Micro stereoangle in direction of sρRadial position vector
    θAngle between s’ and sSOptical path length
    pScattering phase functiongAnisotropic factor
    JcCoherent irradianceI0Incident light intensity
    JdDiffuse irradianceITransmitted light intensity
    QPhoton sourcedOptical depth
    PLaser powerφLuminous flux density
    ω0Spot radiusφ0Incident luminous flux density
    DDiffusion coefficientGPhoton time dispersion
    τLaser pulse durationtCalculating time
    μaAbsorption coefficientRSpecular reflectance
    μtTotal attenuation coefficientμsScattering coefficient
    LPhysical thickness of samplecSpeed of light in medium
    μsReduced scattering coefficientμt’Reduced attenuation coefficient
    kscScattered photon distributionhSpatial angular distribution function
    kTotal photon distribution
    Table 1. Meaning of physical quantities in equations in this paper
    MethodAcceleration relative to standard MCRelative error in simulated optical measurementAdvantageDisadvantage
    SMC200Less than 4%No approximation is made, and it is accurate and fastCurrently only be applicable to hierarchical organization models
    PMC1300Can be less than 4% , depending on magnitude of perturbationIt is applicable to tissuewith complex structuresSensitive to disturbances in scattering
    HMC300Around 5%It has a larger applicablerange than PMCCalculation is relatively complex and specific area must be homogeneous
    Variancereduction300Around 5%There are a variety ofchoices available.Limitation varies withspecific technique
    Table 2. Comparison of characteristics of different accelerated MC methods
    Quantitative analysis methodAdvantageLimitation
    First-order scatteringWhen the diffuse anti-illumination rate is very small; the analytical solution can be obtainedLimit to plane incident waves;not take multiple scattering into account; not be applicable for high scattering bands
    Monte Carlo simulationConsider the light loss of the sample edge that does not match the interface, any phase function of the medium, and the size and arbitrary angular distribution of the incident beamBecause of its statistical reasons, it requires a lot of calculations and a long time to ensure good statistical convergence
    Beam broadening modeMost powerful solution for describing light distribution in approximate high forward scattering mediaEquivalent to MC method in accuracy
    Diffuse approximationHigh calculation speed and high precision;analytical solutions can be given in some specific casesNot be suitable for studying the light transmission and distribution near the light source and the marginal area of biological tissue
    Beam spreading functionConsider the influence of higher-order photons propagating through paths of different lengths;allow to calculate the temporal dispersion of light intensity; much faster than MC methodNot consider the wave characteristics of light
    Inverse adding-doublingIterative solution it can be obtained quickly by computer; it can be used for more than one type of tissue blockCan only be solved under the conditions of limited tissue thickness and uniform optical parameter distribution; the solution process is independent of time
    Kubelka-MunkSimple mathematical formOnly be used for one-dimensional diffuse radiation with scattering far exceeding absorption;only be used be used be suitable for processing one-dimensional structural geometry; theoretical model is not perfect
    Table 3. Classification and characteristic comparison of quantitative analysis methods