# The harmonic bath Hamiltonian

In the theory of chemical reactions, it is often possible to isolate a small number or even a single degree of freedom in the system that can be used to characterize the reaction. This degree of freedom is coupled to other degrees of freedom (for example, reactions often take place in solution). Isomerization or dissociation of a diatomic molecule in solution is an excellent example of this type of system. The degree of freedom of paramount interest is the distance between the two atoms of the molecule - this is the degree of freedom whose detailed dynamics we would like to elucidate. The dynamics of the ``bath'' or environment to which is couples is less interesting, but still must be accounted for in some manner. A model that has maintained a certain level of both popularity and success is the so called ``harmonic bath'' model, in which the environment to which the special degree(s) of freedom couple is replaced by an effective set of harmonic oscillators. We will examine this model for the case of a single degree of freedom of interest, which we will designate \(\underline {q} \). For the case of the isomerizing or dissociating diatomic, \(\underline {q} \) could be the coordinate \(r - \langle r \rangle \), where \(r \) is the distance between the atoms. The particular definition of \(\underline {q} \) ensures that \(\langle q \rangle = 0 \). The degree of freedom \(\underline {q} \) is assumed to couple to the bath linearly, giving a Hamiltonian of the form

\[ H = {p^2 \over 2m} + \phi(q) + \sum_{\alpha}\left[{p_{\alpha}^2 \over 2m_{\alpha} } + {1 \over 2} m_{\alpha} \omega _{\alpha}^2 \left (x_{\alpha} + {g_{\alpha}\over m_{\alpha}\omega_{\alpha}^2}q\right)^2\right] \]

where the index \(\alpha \) runs over all the bath degrees of freedom, \(\underline {\omega _{\alpha} } \) are the harmonic bath frequencies, \(\underline {m_{\alpha } } \) are the harmonic bath masses, and \(\underline {g_{\alpha} } \) are the coupling constants between the bath and the coordinate \(\underline {q} \). \(\underline {p} \) is a momentum conjugate to \(\underline {q} \), and \(m\) is the mass associated with this degree of freedom (e.g., the reduced mass \(\underline {\mu} \) in the case of a diatomic). The coordinate \(\underline {q} \) is assumed to be subject to a potential \(\phi (q) \) as well (e.g., an internal bond potential). The form of the coupling between the system ( \(\underline {q} \) ) and the bath (\(\underline {x_{\alpha}} \)) is known as *bilinear*.

Below, using a completely classical treatment of this Hamiltonian, we will derive an equation for the detailed dynamics of \(\underline {q} \) alone. This equation is known as the generalized Langevin equation (GLE).