“…A persistent challenge for quantitative prediction of Franck–Condon integrals, particularly for systems of high dimension, has been the computational time and associated resources (e.g., memory, disk storage) required for determination of the FCFs. ,, This is why early efforts were limited to 1-D methods, , and determination of 2-, , 3-, and 4-dimensional FC overlap integrals, and only for small molecules, primarily diatomics. ,,− Only in the past decade have efforts been extended toward the simplification of the computational methodology for more efficient and less costly analytic and recursion approaches, including our recent demonstration involving systems up to 30-D. , What remains is greater accessibility to much larger dimension systems (>50-D), which we have tackled in this work.…”