Development and clinical verification of numerical simulation for laser in situ keratomileusis

To develop and validate numerical models of the laser in situ keratomileusis (LASIK) procedure through considering its effect on corneal biomechanical behavior. 3D finite element models of the human eye were developed to simulate LASIK. The models' predictions of post-operative corneal elevation, corneal refractive power with vector decomposition (M_-c-pos, J_0-c-pos, J_45-c-pos) and refractive error correction (M_-rec, J_0-rec, J_45-rec) were compared against clinical data obtained for 28 eyes of 28 patients. A parallel exercise was conducted to estimate the post-operative corneal shape using a shape subtraction method (SSM) - which does not consider the effects of LASIK on corneal mechanical behavior - and the results are compared with the finite element method (FEM). A significant decrease in elevation differences between FEM predictions and clinical data was found compared with the differences between SSM results and clinical data (p = 0.000). In addition, there were no significant differences in post-operative equivalent sperical corneal refractive power between FEM results and corresponding clinical data (M_-c-pos: p = 0.501), while SSM showed significant differences with clinical data (M_-c-pos: p = 0.000). Further, FEM achieved a predicted value of M_-c-pos within ± 1.00D accuracy in 100% of cases, compared with 57% achieved by the SSM. M_-rec predicted by FEM was not significantly different from clinical results (p = 0.085), while SSM overestimated it (p = 0.000). The match between LASIK numerical model predictions with clinical measurements improved significantly when the procedure's effect on corneal biomechanical behavior was considered. This outcome has important implications on efforts to develop planning tools for refractive surgery.