GOAL: To examine the effect of Optical Properties, Solver and Source Configuration on Diffuse Reflectance.

I. Impact of Optical Properties on Spatially-Resolved Reflectance

1. On the "Forward/Analysis" Panel select: (a) Forward Model: Standard Diffusion (Analytic - Point Source) and (b) Solution Domain: Steady State (select R(ρ)).
2. Select start and stop locations to 0.5 and 9.5 mm, respectively with 19 points (every 0.5 mm).
3. In Optical Properties: enter μ_a = 0.01mm^-1, μ'_s=1mm^-1.
4. Press "Plot Reflectance" button.
5. Press "Hold On".
6. Fix μ_s'=1mm^-1. Repeat the above steps for μ_a = 0.1 and 1.0 mm^-1.
7. Please note the trend of decreasing reflectance with increasing absorption.
8. Now toggle the plots with a logarithmic y-axis spacing. Note the linear behavior at larger detector location.

   Question: Can you relate this to the underlying analytic approximation?

9. Press "Clear All" and toggle back to "Linear" y-axis spacing.
10. Select start and stop locations to 0.5 and 9.5 mm, respectively with 19 points (every 0.5 mm).
11. In Optical Properties: enter μ_a = 0.01mm^-1, μ'_s = 1 mm^-1.
12. Press "Plot Reflectance" button.
13. Press "Hold On".
14. Fix μ'_s=1mm^-1. Repeat the above steps for μ'_s = 0.5 and 1.5 mm^-1.
15. Now toggle the plots with a logarithmic y-axis spacing.

    Questions: Note the trend of increasing reflectance with increasing scattering close to the source but the opposite far from the source. Is this expected? Why or why not?

II. Compare SDA and MC predictions for Spatially-Resolved Reflectance

1. On the "Forward/Analysis" Panel select: (a) Forward Model: Standard Diffusion (Analytic - Point Source) and (b) Solution Domain: Steady State (select R(ρ)).
2. Select start and stop locations to 0.5 and 9.5 mm, respectively with 19 points (every 0.5 mm).
3. In Optical Properties: enter μ_a = 0.01mm^-1, μ'_s = 1 mm^-1.
4. Plot Reflectance.
5. Press "Hold On".
6. Now select: Forward Model: Scaled Monte Carlo (g=0.8).
7. Plot Reflectance.
8. Repeat the steps for μ_a = 0.1 and 1 mm^-1.

Questions:

1. How do the SDA and MC models compare close to ρ = 0?
2. Now switch to a logarithmic y-axis spacing. How do the models compare far from the source?

Press "Clear All" and return to linear axis spacing.

III. Impact of Optical Properties on Temporally-Resolved Reflectance

1. On the "Forward/Analysis" Panel select: (a) Forward Model: Scaled Monte Carlo (g=0.8) and (b) Solution Domain: Time-Domain, Collimated Beam Source (select R(ρ,t)).
2. For the Independent axis choose t and set ρ = 10 mm.
3. Choose "Start" = 0 ns and "Stop" = 0.5 ns with 101 time points (1 point every 5 ps).
4. In Optical Properties: enter μ_a = 0.01mm^-1, μ'_s = 1mm^-1.
5. Plot time-resolved reflectance, R(t).
6. Repeat the above for μ_a = 0.1 and 1 mm^-1.

   Question: Note also the difference in the amplitude and shape of these plots. What do you believe is responsible for this? Hint: It may help to use both linear and log y-axis spacing.

7. Press "Clear All" and return and In Optical Properties: enter μ_a = 0.01mm^-1, μ'_s = 1mm^-1.
8. Plot time-resolved reflectance, R(t).
9. Repeat the above for μ'_s = 0.5 and 1.5 mm^-1.

Questions:

1. Note that no photons are detected before a finite time in the time-resolved reflectance signal. Can you independently calculate the minimal delay time? Note: It may help to visualize this delay using logarithmic x-axis spacing.
2. Note that the peak reflectance values are different and not located at the same time point. Can you speculate as to the origin of these features? Hint: It may help to use both linear and log y-axis spacing.
3. You are designing a time-resolved optical imaging system to detect early formation of a fibrous tumor. What is an approximate time resolution and source detector separation necessary to differentiate normal breast with μ'_s=0.6 mm^-1 from a fibroid tumor with μ'_s=1.2 mm^-1 with such a system?