The experimental results of the diffusion coefficients for Rf and oxygen are found in
Table 1.
Rf diffusion coefficients of 4 × 10
−7 cm
2/s were obtained at 25°C and 6 × 10
−7 cm
2/s at 35°C for a 0.1% Rf solution in distilled water. These values are very similar to those found in the literature. Spoerl et al.
12 have used a value of 6.5 × 10
−7 cm
2/s, on the basis of previously measured diffusion coefficient for sodium fluorescein. Because Rf and sodium fluorescein have similar polarity and molecular weight, it was assumed that the diffusivity of Rf would be similar. Araie and Maurice
14 have reported diffusivity of fluorescein in rabbit corneas of 6 × 10
−7 cm
2/s at 37°C in an earlier study, and Nagataki et al.
15 reported a value of (1.2 ± 0.2) × 10
−6 cm
2/s at 19°C (in this case, the value represents diffusion along [parallel to] the stroma and not across the stroma, as in our case). After a 20-minute instillation with Rf, the flap concentration at 100 μm was more than 80% of the Rf concentration applied at the corneal surface, and it was assumed that the Rf had reached a steady-state value, thus simplifying the theoretical calculations.
Illumination with UV-A irradiation caused a rapid (within several seconds) depletion of oxygen in Rf-containing cornea, and turning the UV light off led to replenishment of the oxygen to its original level within 3 to 4 minutes (
Fig. 4). In these experiments, the observed period needed to reach a steady-state flow of oxygen through a corneal flap of known thickness matched expected oxygen diffusivity in the cornea. Oxygen diffusion coefficients of 4 × 10
−6 cm
2/s at 25°C and 6 × 10
−6 cm
2/s at 35°C were obtained. Diffusion of oxygen in the cornea has been studied by numerous authors. The diffusivity values vary because of their strong dependence on corneal hydration and temperature, which are difficult to control. A recent model for human corneal stroma reported by Larrea and Buchler
16 uses an oxygen diffusion coefficient of 28 × 10
−6 cm
2/s (35°C), although earlier studies
17–19 all have used an in vivo value close to 7 × 10
−6 cm
2/s. The oxygen diffusion coefficient for rabbit corneal stroma has been reported as (5 ± 1.8) × 10
−6 cm
2/s (25°C) by Takahashi et al.,
20 and Harvitt and Bonnano
21 have used a value close to 6 × 10
−6 cm
2/s in a later study.
Figure 5 shows the time-dependent depletion of dissolved oxygen below the 100-μm thick corneal flap with Rf cross-linking for a 3 mW/cm
2 UV-A irradiance at 25°C and 35°C both theoretically and experimentally; the oxygen concentration (mg/L) fell to zero at 15 and 10 seconds, respectively. The theoretically modeled results were obtained by using parameters found in
Table 1. The theoretical and experimental results correlated particularly well. (The mean square error [MSE] and the coefficient of determination [
R 2] were 0.14 and 0.98, respectively, for 25°C; 0.33 and 0.90, respectively, for 35°C).
Figure 6 shows the time-dependent depletion of dissolved oxygen below the 100-μm thick corneal flap with Rf cross-linking for a 3 mW/cm
2 and 30 mW/cm
2 UV-A irradiance at 25°C both theoretically and experimentally; the oxygen concentration (mg/L) fell to zero at 15 seconds and 5 seconds, respectively. The response time of the oxygen sensor was 2 seconds per data point, which limited the ability to correlate to the theoretical model. Again, the theoretically modeled results were obtained by using parameters found in
Table 1. The theoretical and experimental results correlated particularly well, even with the response limits of the detector (for 30 mW/cm
2, MSE = 0.64,
R 2 = 0.89).
Figure 7 shows the time-dependent depletion and gradual replenishment of dissolved oxygen below a 100-μm thick corneal flap with Rf cross-linking for a 3 mW/cm
2 UV-A irradiance at 25°C. The oxygen concentration (mg/L) fell to zero at 15 seconds and gradually started to increase after approximately 10 minutes, getting back to approximately one-tenth its starting value after 30 minutes.
At 30 mW/cm
2, a photodegradation reaction of Rf in this model was disregarded, as it seemed that it did not matter when comparing the model's result with the experimental data for the first 10 seconds of the irradiation (
Fig. 6). Two major products of Rf photodegradation, lumiflavin and lumichrome, are reported to be less efficient than Rf in reaction with oxygen producing superoxide radical anions.
22 Additionally, lumichrome has a 60% lower yield in producing singlet oxygen than riboflavin
22 and therefore, the overall effect on oxygen depletion would be weaker in the authors' model, which maintained a steady Rf level and did not include generation of lumiflavin or lumichrome.