We show that multichannel quantum defect theory (MQDT) can be applied successfully as an efficient computational method for cold molecular collisions in
Li+NH, which has a deep and strongly anisotropic interaction potential. In this strongly coupled system, closed-channel poles restrict the range over which the MQDT $\mathbit{Y}$ can be interpolated. We present an improved procedure to transform the MQDT reference functions so that the poles are removed from the energy range of interest. Effects due to very-long-range spin dipolar couplings are outside the scope of MQDT, but can be added perturbatively. This procedure makes it possible to calculate the elastic and inelastic cross sections, over the entire range of energies and fields needed to evaluate the feasibility of sympathetic cooling of NH by Li, using coupled-channel calculations at only five combinations of energy and field.
Multichannel quantum defect theory (MQDT) can provide an efficient alternative to full coupled-channel calculations for low-energy molecular collisions. However, the efficiency relies on interpolation of the $\mathbit{Y}$ matrix that encapsulates the short-range dynamics, and there are poles in $\mathbit{Y}$ that may prevent interpolation over the range of energies of interest for cold molecular collisions. We show how the phases of the MQDT reference functions may be chosen so as to remove such poles from the vicinity of a reference energy and dramatically increase the range of interpolation. For the test case of Mg+NH, the resulting optimized $\mathbit{Y}$ matrix may be interpolated smoothly over an energy range of several Kelvin and a magnetic field range of over 1000 gauss. Calculations at additional energies and fields can then be performed at a computational cost that is proportional to the number of channels
N and not to
N3.