Theoretical treatment of channel mixing in excited [Formula Presented] and [Formula Presented] ultracold molecules: Determination of predissociation lifetimes with coordinate mapping

The treatment of the dynamics of ultracold molecules requires new theoretical tools. Previous work of the present authors [J. Chem. Phys. 110, 9865 (1999)] for calculation of vibrational levels by a Fourier grid representation with use of adaptative coordinates is generalized here to the treatment of the bound-continuum interaction in a two-channel problem. Two numerical methods are presented: a time-dependent method using a Chebyshev propagator to compute the correlation function and a time-independent method with diagonalization of a Hamiltonian that includes an absorbing optical potential. In both cases the adaptative coordinate is defined by a numerical rather than an analytical procedure. Lifetimes are reported for the predissociated levels of the [Formula Presented] and [Formula Presented] spectra, where [Formula Presented] The two numerical methods give similar results. The lifetimes increase with the vibrational quantum number proportionally to the classical vibration period estimated from the Le Roy–Bernstein law for an asymptotic [Formula Presented] potential, and the energy variation can be fitted to an analytical formula. The results are shown to be very sensitive to the molecular parameters, potentials, and couplings. The measured width of 8.5 GHz reported by Cline et al. [Phys. Rev. Lett. 73, 632 (1994)] for one predissociated level of [Formula Presented] is reproduced. A strong isotopic effect is found for the rubidium dimer, the lifetimes of [Formula Presented] and [Formula Presented] levels differing by a factor of 3. Finally, we present a third approach, in the framework of a generalized two-channel quantum-defect theory, where lifetimes are determined by extrapolation of parameters fitted to Lu-Fano plots of computed bound levels below the [Formula Presented] dissociation limit. Excellent agreement is obtained with the numerical results, suggesting the possibility of fitting to experimental spectra.