Yes it should be. Here is an example for a CCSD calculation of a custom (real) Hamiltonian: https://github.com/pyscf/pyscf/blob/master/examples/cc/40-ccsd_custom_hamiltonian.py

Depending on the way in which your Hamiltonian is complex, you'll have to choose an appropriate mean-field and CC method to use.

I mean Hamiltonian, Fock matrix with complex numbers. Provided examples works only with real Hamiltonian and eris

Does it break spin symmetry (like spin-orbit would) or no?