That said, our preferred choice for complexes not involving ions is XSAPT + MBD/def2-ma-TZVP, with monomer wavefunctions obtained using the LRC-ω GDD PBE functional, tuned in a monomer-specific way, and with CM5 embedding charges. The latter perform well for ions, and we recommend XSAPT + ai D3 as a starting point for ionic complexes. These methods are available in the Q-Chem code, with improved multithreading capabilities in recent releases.