Skip to main content
Chemistry LibreTexts

14.3: The Dynamic Friction Kernel

  • Page ID
    5308
  • \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    \( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

    ( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\id}{\mathrm{id}}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\kernel}{\mathrm{null}\,}\)

    \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\)

    \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\)

    \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

    \( \newcommand{\vectorA}[1]{\vec{#1}}      % arrow\)

    \( \newcommand{\vectorAt}[1]{\vec{\text{#1}}}      % arrow\)

    \( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vectorC}[1]{\textbf{#1}} \)

    \( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

    \( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

    \( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

    \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    \(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)

    The convolution integral term

    \[ \int_0^t\;d\tau\;\dot{q}(\tau)\zeta(t-\tau) \nonumber \]

    is called the memory integral because it depends, in general, on the entire history of the evolution of \( {q} \). Physically it expresses the fact that the bath requires a finite time to respond to any fluctuation in the motion of the system (\( {q} \)). This, in turn, affects how the bath acts back on the system. Thus, the force that the bath exerts on the system presently depends on what the system coordinate \( {q} \) did in the past. However, we have seen previously the regression of fluctuations (their decay to 0) over time. Thus, we expect that what the system did very far in the past will no longer the force it feels presently, i.e., that the lower limit of the memory integral (which is rigorously 0) could be replaced by \( {t-t_{\rm mem}} \), where \( {t _{mem} }\) is the maximum time over which memory of what the system coordinate did in the past is important. This can be interpreted as a indicating a certain decay time for the friction kernel \( \zeta(t) \). In fact, \(\zeta(t) \) often does decay to 0 in a relatively short time. Often this decay takes the form of a rapid initial decay followed by a slow final decay, as shown in the figure below:

    Consider the extreme case that the bath is capable of responding infinitely quickly to changes in the system coordinate \( {q} \). This would be the case, for example, if there were a large mass disparity between the system and the bath ( \( {m \gg m_{\alpha} } \)). Then, the bath retains no memory of what the system did in the past, and we could take \(\zeta(t) \) to be a \(\delta\)-function in time:

    \[ \zeta(t) = 2\zeta_0\delta(t) \nonumber \]

    Then

    \[ \int_0^t\;d\tau\;\dot{q}(\tau)\zeta(t-\tau) =\int_0^t\;d\tau \dot {q} (t - \tau ) \zeta (\tau) = 2 \zeta _0 \int _0^t\;d\tau\;\delta(\tau)\dot{q}(t-\tau) = \zeta_0\dot{q}(t) \nonumber \]

    and the GLE becomes

    \[ m\ddot{q} = -{\partial \phi \over \partial q} - \zeta_0\dot{q} + R(t) \nonumber \]

    This simpler equation of motion is known as the Langevin equation and it is clearly a special case of the more generalized equation of motion. It is often invoked to describe Brownian motion where clearly such a mass disparity is present. The constant \(\zeta _0\) is known as the static friction and is given by

    \[ \zeta_0 = \int_0^{\infty}\;dt\;\zeta(t) \nonumber \]

    In fact, this is a general relation for determining the static friction constant.

    The other extreme is a very sluggish bath that responds slowly to changes in the system coordinate. In this case, we may take \(\zeta (t) \) to be a constant \(\zeta \equiv \zeta(0) \), at least, for times short compared to the response time of the bath. Then, the memory integral becomes

    \[ \int_0^t\;d\tau\;\dot{q}(\tau)\zeta(t-\tau) \approx \zeta(q(t)-q(0)) \nonumber \]

    and the GLE becomes

    \[ m\ddot{q} = -{\partial \over \partial q}\left(\phi(q) +{1 \over 2}\zeta (q-q_0)^2\right) + R(t) \nonumber \]

    where the friction term now manifests itself as an extra harmonic term added to the potential. Such a term has the effect of trapping the system in certain regions of configuration space, an effect known as dynamic caging. An example of this is a dilute mixture of small, light particles in a bath of heavy, large particles. The light particles can get trapped in regions of space where many bath particles are in a sort of spatial 'cage.'' Only the rare fluctuations in the bath that open up larger holes in configuration space allow the light particles to escape the cage, occasionally, after which, they often get trapped again in a new cage for a similar time interval.


    This page titled 14.3: The Dynamic Friction Kernel is shared under a CC BY-NC-SA 4.0 license and was authored, remixed, and/or curated by Mark E. Tuckerman.