This paper introduced a new formulation of the dispersion interaction from the Longuet-Higgins type integral over the product of coupled-cluster polarization propagators of the monomers. Unlike previous formulations, instead of numerical integration over imaginary frequencies, we performed analytical integration and introduced a new type of dispersion excitation operator σ. The σ operator can be obtained from the solution of the linear equations, similarly to CC response theory. The general, nested commutator expansion, The new formulation of the dispersion energy might find applications in the future for highly accurate calculations of the interaction potentials. A big advantage over CCD+ST(CCD) method and TD-DFT approaches is its systematic improvability: the general equation (18) was derived for arbitrary order of T . However, similarly to time-independent response theory it is possible to introduce truncation schemes valid to desired level of W operator.