Summary: The recent advances in alternating direct implicit (ADI) methods promise important new capability for time domain plasma simulations, namely the elimination of numerical stability limits on the time step. But the utility of these methods in simulations with charge and current sources, such as in electromagnetic particle-in-cell (EMPIC) computations, has been uncertain, as the methods introduced so far do not have the property of divergence preservation. This property is related to charge conservation and self-consistency, and is critical for accurate and robust EMPIC simulation. This paper contains a complete study of these ADI methods in the presence of charge and current sources. It is shown that there are four significantly distinct cases, with four more related by duality. Of those, only one preserves divergence and, thus, is guaranteed to be stable in the presence of moving charged particles. Computational verification of this property is accomplished by implementation in existing 3D-EMPIC simulation software. Of the other three cases, two are verified unstable, as expected, and one remains stable, despite the lack of divergence preservation. This other stable algorithm is shown to be related to the divergence preserving case by a similarity transformation, effectively providing the complement of the divergence preserving field in the finite-difference energy quantity.