Abstract. Employing separate cluster ansatz in time-independent and time-dependent wave-operators, coupled-cluster (CC) response theory is generalized to multireference (MR) expansion spaces. For state energies, this corresponds to the MR secular problem with an arbitrary similarity-transformed effective Hamiltonian, H=1H. The effective Hamiltonian can be generated via size-extensive CC methods. Thus the states in MR linear response theory (MRLRT) maintain the usual CC core-extensive properties. We have used the Gelfand unitary group basis of the spin-adapted configurations to construct the matrix of H in the MR excitation space. As a preliminary application, the CC singles and doubles effective Hamiltonian is applied to excitation and photoionization energies of the CH+ and N2 molecules, and is compared with experimental results and results from other numerical procedures including conventional CC linear response theory (CC-LRT), MR and full configuration interaction (MRCI and FCI) methods. The numerical results indicate that MRLRT reproduces valence and external excited states quantitatively, combining the best features of CC-LRT and MRCI.