We study a series of N oscillators, each coupled to its nearest neighbors, and linearly to a phonon field through the oscillator's number operator. We show that the Hamiltonian of a pair of adjacent oscillators, or a dimer, within the series of oscillators can be transformed into a form in which they are collectively coupled to the phonon field as a composite unit. In the weak coupling and rotating-wave approximation, the system behaves effectively as the trilinear boson model in the one excitation subspace of the dimer subsystem. The reduced dynamics of the one excitation subspace of the dimer subsystem coupled weakly to a phonon bath is similar to that of a two-level system, with a metastable state against the vacuum. The decay constant of the subsystem is proportional to the dephasing rate of the individual oscillator in a phonon bath, attenuated by a factor that depends on site asymmetry, intersite coupling, and the resonance frequency between the transformed oscillator modes, or excitons. As a result of the collective effect, the excitation relaxation lifetime is prolonged over the dephasing lifetime of an individual oscillator coupled to the same bath.
We study the reduced dynamics of a pair of nondegenerate oscillators coupled collectively to a thermal bath. The model is related to the trilinear boson model where the idler mode is promoted to a field. Due to nonlinear coupling, the Markovian master equation for the pair of oscillators admits non-Gaussian equilibrium states, where the modes distribute according to the Bose-Einstein statistics. These states are metastable before the nonlinear coupling is taken over by linear coupling between the individual oscillators and the field. The Gibbs state for the individual modes lies in the subspace with infinite occupation quantum number. We present the time evolution of a few states to illustrate the behaviors of the system.
We consider the reduced dynamics of a molecular chain weakly coupled to a phonon bath. With a small and constant inhomogeneity in the coupling, the excitation relaxation rates are obtained in closed form. They are dominated by transitions between exciton modes lying next to each other in the energy spectrum. The rates are quadratic in the number of sites in a long chain. Consequently, the evolution of site occupation numbers exhibits longer coherence lifetime for short chains only. When external source and sink are added, the rate equations of exciton occupation numbers are similar to those obtained earlier by Fröhlich to explain energy storage and energy transfer in biological systems. There is a clear separation of timescale into a faster one pertaining to internal influence of the chain and phonon bath, and a slower one determined by external influence, such as the pumping rate of the source, the absorption rate of the sink, and the rate of radiation loss. The energy transfer efficiency at steady state depends strongly on these external parameters and is robust against a change in the internal parameters, such as temperature and inhomogeneity. Excitations are predicted to concentrate to the lowest energy mode when the source power is sufficiently high. In the site basis, this implies that when sustained by a high power source, a sink positioned at the center of the chain is more efficient in trapping energy than a sink placed at its end. Analytic expressions of energy transfer efficiency are obtained in the high power and low-power source limit. Parameters of a photosynthetic system are used as examples to illustrate the results.