Generally, one has to truncate the L-VCI expansion. Here, we use an L-VCISD expansion, which includes the ground state as well as singly and doubly excited states, which is the minimum necessary for modeling 2D IR spectroscopy. This approximation is usually sufficient when considering weakly coupled localized modes. By diagonalizing the L-VCI matrix,