First-order nonadiabatic coupling matrix elements using coupled cluster methods. I. Theory