One can expect that using correlated monomer states, which already include local strong correlation effects, as a basis to represent the aggregate state (Eq. ( )) should mitigate this problem. Note that, in principle, a full account of the correlation effects including the nonlocal static ones is possible within Frenkel formalism with the two-particle (or in general n-particle) basis in the spirit of Eq. ( ) and systematic inclusion of couplings between all possible states. However, in this case, one again approaches the almost exact supermolecular calculation (excluding intermonomer charge transfer transitions) and loses the main advantage of the approach -its efficiency. Nevertheless, the systematic increase of the number of included interstate couplings is a way to improve the description; this strategy is not feasible for linear-response TDDFT, since many of the coupled states will have more than one-electron difference between them, which renders them inaccessible for a linear-response method.