Salve, sono alle prese con lo studio riguardante transizioni fra stati elettronici molecolari in potenziali di Morse cercando di calcolare i Fattori di Franck-Condon in via analitica. Premettendo che non sono del tutto sicuro del risultato a cui sono arrivato volevo chiedervi aiuto per la risoluzione di tale integrale in modo da poter avere una formula facilmente programmabile.
\(\displaystyle \int_0^\infty e^{-(x/2)(1+e^{\beta'd})} x^{(1/2)(a+a')} L^{a}_{n}(x) L^{a'}_{n'}(xe^{\beta'd)} dx\)
dove \(\displaystyle x= C e^{-\beta(r-r_{e})} \) e gli \(\displaystyle L^{a}_{n}(x) \)sono i polinomi di Laguerre generalizzati. Sono arrivato a questo integrale partendo dalle autofunzioni del potenziale di Morse e andando a calcolare l'integrale di sovrapposizione fra due di queste aventi parametri diversi \(\displaystyle a, a', \beta, \beta', n, n', x, x', r_{e}, r_{e}', C, C'\) e sapendo che \(\displaystyle r_{e}- r_{e}'=d\) da cui, effettuando la sostituzione, sono arrivato al risultato di sopra. Ho trovato testi poco chiari che risolvono sfruttando funzione gamma di Eulero e la equazione ipergeometrica confluente (entrambe facilmente programmabili) ma seguono percorsi molto diversi da quello da me intrapreso. Qualsiasi consiglio sulla risoluzione o su strade più percorribili è ben accetto. Grazie in anticipo.