I'm reading this article about coherent exciton transport in photosynthetic light harvesting and the role of quantized vibrations. Along the way, I came across a section where the article claimed the exciton-phonon Hamiltonian (that is, the Hamiltonian representing exciton-phonon coupling) can be represented as follows:
$$H_{\text{elec-vib}}=\sum_nL_n\sum_\lambda\kappa_{n\lambda}(a_{n\lambda}^\dagger+a_{n\lambda}),$$where the system operators $L_n$ are the projectors $L_n=-|n\rangle\langle n|$. The coupling constants $\kappa_{n\lambda}$... describe the strength of the coupling of electronic excitation of pigment $n$ to the vibrational mode $\lambda$ of this pigment.
The article also notes that $a_{n\lambda}^\dagger$ is the phonon creation operator for the $\lambda^{\text{th}}$ vibrational mode of the $n^{\text{th}}$ pigment. My question is: where does the expression $(a_{n\lambda}^\dagger+a_{n\lambda})$ come from? Why not just use the operator $a_{n\lambda}^\dagger a_{n\lambda}$ to "count" the number of phonons on a given mode and then multiply that number by their coupling strength to the excitons? The operator $(a_{n\lambda}^\dagger+a_{n\lambda})$ seems very weird and unintuitive to me, and I can't even figure out what its eigenvalues would be.