Mathematical model of corneal surface smoothing after laser refractive surgery

Purpose To construct a quantitative model of corneal surface smoothing after laser ablation for refractive correction. Design Experimental study, interventional case series, and meta-analysis of literature. Methods A theory of epithelial smoothing in response to corneal contour change is derived from differential equations that describe epithelial migration, growth, and loss. Computer simulations calculate the effects on postoperative epithelial thickness, topography, refraction, and spherical aberration. Model parameter is matched with laser in situ keratomileusis (LASIK) outcome in literature and in a retrospective study of primary spherical myopic (77 eyes) and hyperopic (19 eyes) corrections. Surgically induced refractive change was the main outcome measure. Results Simulated epithelial remodeling after myopic ablation produces central epithelial thickening, reduction in achieved correction, and induction of oblate spherical aberration. Simulation of hyperopic ablation shows peripheral epithelial thickening, a larger reduction in correction, and induction of prolate spherical aberration. Simulation using a minus cylinder laser ablation pattern shows decreased astigmatism correction and increased hyperopic shift. In the LASIK series, linear regression of achieved correction vs ablation setting in hyperopic and minus cylinder corrections shows slopes of 0.97, 0.71, and 0.74, respectively. These clinical results match model predictions when the smoothing constant is set at 0.32, 0.63, and 0.55 mm, respectively. Conclusions Epithelial thickness modulations after ablation can be modeled mathematically to explain clinically observed regression and induction of aberration. The cornea appears to smooth over ablated features smaller than approximately 0.5 mm. The model provides an approach for designing ablation patterns that precompensate for the smoothing to improve final outcome.