This is a microcanonical ensemble approach to a simple model of a two dimensional polymer bundle. The professor of the course I took does a lot of research in the area of polymer physics and so set a few problems pertaining to them. They aren’t easily found in textbooks or online either (this website notwithstanding) but are nevertheless quite interesting in themselves. We also approached this same problem from the canonical ensemble which I’ll upload soon.

Question

A polymer is represented by a chain of \(N\) segments of length \(a\). \(N_x^+\) segments are oriented in the positive \(x\)-direction \(N_x^-\) are oriented in the negative \(x\)-direction. \(N_y^+\) and \(N_y^-\) segments are oriented in the positive and negative \(y\)-direction, respectively.

  1. Assign an energy penalty \(\varepsilon\) for every polymer segment oriented in \(y\)-direction. This penalty reflects the confinement in a polymer bundle.
  2. For prescribed extensions \(L_x\) and \(L_y\) and energy \(E\), calculate the entropy \(S(E,N,L_x,L_y)\) of the chain.
  3. Set \(L_y=0\) in the following. Derive the mechanical equation of state for the force \(F_x(E,N,L_x)\) in the \(x\)-direction.
  4. Derive the thermal equation of state for the temperature \(T(E,N,L_x)\).
  5. Eliminate \(E/\varepsilon\) from the two equations of state and show that \begin{equation} L_x(T,N,F_x)=aN\frac{\sinh(\frac{F_x a}{k_B T})}{\cosh(\frac{F_x a}{k_B T}) + \exp(\frac{-\varepsilon}{k_B T})} \label{eq:lx} \end{equation}
  6. Expand \(L_x\) for small force and extract the entropic spring constant.

Solution

1. Energy Penalty

To include an energy penalty for segments oriented in the \(y\)-direction, we simply write \(E=\varepsilon [N_y^++N_y^-]\).

2. Entropy

To calculate the entropy of the chain, we first need to deduce the number of states. In the one-dimensional case this is simply the result for the one-dimensional random walk: $$\Omega=\frac{N!}{N^+! N^-!}$$

Extending this to two dimensions, we obtain

\begin{equation} \Omega=\frac{N!}{N_x^+! N_x^-! N_y^+! N_y^-!} \label{eq:omega} \end{equation}

We find the entropy using the well-known formula \(S=k_B \ln \Omega\) and applying Stirling’s approximation: $$S = k_B \left[ N\ln N - N_x^+ \ln N_x^+ - N_x^- \ln N_x^- - N_y^+ \ln N_y^+ - N_y^- \ln N_y^- \right]$$

If we define \(\rho_{x/y}^{\pm}=\frac{N_{x/y}^{\pm}}{N}\), we can write it a bit more succinctly:

\begin{equation} S = -k_B N \left[ \rho_x^+ \ln \rho_x^+ + \rho_x^- \ln \rho_x^- + \rho_y^+ \ln \rho_y^+ + \rho_y^- \ln \rho_y^- \right] \label{eq:srho} \end{equation}

We set \(L_y=0\). Now let us briefly take pause to define some variables for convenience:

For convenience, \(\eta := \frac{E}{N \varepsilon}\), \(x:=\frac{L_x}{L_0}\), and \(L_0:=Na\), where \(L_0\) is the “stretched out” (i.e. \(N=N_x\)) length of the polymer.

We also note that the effective number of steps in the positive \(x\)-direction is the effective length over the segment length: \(N_x^+ -N_x^-=\frac{L_x}{a}\) and the same for \(y\): \(N_y^+ -N_y^-=\frac{L_y}{a}\).

Combining these variables, it is possible to write \[\frac{N_x^{\pm}}{N}=\frac 12 \big[1-\eta \pm x\big]\]

3. Mechanical Equation of State

So, to the question. The equation for the force in the \(x\)-direction \(F_x\) is given by \(\frac{-F_x}{T}=\frac{\partial S}{\partial L_x}\), so

\begin{equation} \frac{-F_x}{T}=\frac{\partial S}{\partial L_x}=\frac{\partial S}{\partial x}\frac{\partial x}{\partial L_x}=\frac{1 }{L_0}\frac{\partial S}{\partial x} \label{eq:fx} \end{equation}

In order to find the derivative \(\frac{\partial S}{\partial x}\) we need to rewrite the entropy with the new constraint of \(L_y=0\).

\begin{equation} \frac{S}{-k_B N} = \frac{1}{2} \left(1-\eta + x \right)\ln \left(1-\eta + x \right) + \frac{1}{2} \left(1- \eta – x \right)\ln \left(1-\eta – x \right) + 2 \left( \frac{\eta}{2} \right) \ln \left( \frac{\eta}{2} \right) \label{eq:news} \end{equation}

Taking the derivative for \(x\) then, we find $$\frac{\partial S}{\partial x}=-\frac{k_B N}{2} \ln \left[\frac{1-\eta + x}{1-\eta - x}\right]$$

The mechanical equation of state for the force \(F_x\), found by rearranging \eqref{eq:fx} and using the result just obtained, is then:

\begin{equation} F_x=\frac{k_B T}{2a} \ln \left[ \frac{1-\eta + \frac{L_x}{L_0}}{1-\eta - \frac{L_x}{L_0}} \right] \label{eq:fx2} \end{equation}

4. Thermal Equation of State

To obtain the thermal equation of state, we use the formula \( \frac 1T = \frac{\partial S}{\partial E} \)

Using \eqref{eq:news}:

\[\frac{2\varepsilon}{k_B T} = \ln \left[ \frac 12 (1-\eta + x) \right] + \ln \left[ \frac 12 (1-\eta - x) \right] -2 \ln \frac \eta 2\]

which we can write much more neatly as:

\begin{equation} \frac{2\varepsilon}{k_B T} = \ln \left( \frac{1-\eta + x}{\eta} \right) + \ln \left( \frac{1-\eta , – , x}{\eta} \right) \label{eq:dsde} \end{equation}

5. Calculating the Extension

The equation we are to derive is quite an long process as the equations get somewhat bulky. In fact, the same result can be obtained much simpler (just a few lines) in the canonical ensemble approach to this problem. The steps here are only an outline as a guide.

Start by solving \eqref{eq:fx2} for \(\eta\):

\begin{equation} \eta=1-x·\coth\frac{F_x a}{k_B T} \label{eq:etacoth} \end{equation}

If we now sub. this into equation \eqref{eq:dsde}:

\[\frac{2 \varepsilon}{k_B T}=\ln \left[ \frac{\frac{L_x}{Na}\coth\left(\frac{F_x a}{k_B T} \right) +\frac{L_x}{Na}}{1-\frac{L_x}{Na}\coth\left(\frac{F_x a}{k_B T} \right)} \right] + \ln \left[ \frac{\frac{L_x}{Na}\coth\left(\frac{F_x a}{k_B T} \right) -\frac{L_x}{Na}}{1-\frac{L_x}{Na}\coth\left(\frac{F_x a}{k_B T} \right)} \right]\]

Simplifying and taking the exponent of both sides:

\[e^{\frac{2 \varepsilon}{k_B T}}=\frac{\frac{L_x^2}{N^2 a^2} \coth^2 \left( \frac{F_x a}{k_B T} \right) – \frac{L_x^2}{N^2 a^2} }{ \left[ 1-\frac{L_x}{N a} \coth \left( \frac{F_x a}{k_B T} \right) \right]^2}\]

Using that \(\coth^2(x)-1=\frac{1}{\sinh^2(x)}\):

\[e^{\frac{\varepsilon}{k_B T}}=\frac{\frac{L_x}{N a} \sqrt{\coth^2 \left(\frac{F_x a}{k_B T}\right) -1}}{1-\frac{L_x}{N a} \coth \left( \frac{F_x a}{k_B T} \right)}\]

Finally, we solve for \(L_x\):

\begin{equation} \frac{L_x}{Na}=\frac{\sinh(\frac{F_x a}{k_B T})}{\cosh(\frac{F_x a}{k_B T}) + \exp(\frac{-\varepsilon}{k_B T})} \label{eq:finally} \end{equation}

which is the equation \eqref{eq:lx} given.

6. Extracting the Entropic Spring Constant

We need to take the limit of \eqref{eq:lx} for small force. For small \(x\): \(\sinh(x) \to x\) and \(\cosh(x) \to 1\).

Therefore \(L_x\) becomes

\[L_x \approx aN \frac{\left(\frac{F_x a}{k_B T}\right)}{1+\exp \left(\frac{-\varepsilon}{k_B T}\right)}\]

If we re-arrange for \(F_x\):

\[F_x=\frac{k_B T}{a^2N}\left[1+\exp \left(\frac{-\varepsilon}{k_B T}\right) \right] L_x\]

And compare it with Hooke’s Law \(F=-kx\):

\[k_{\text{entropy}}=\frac{k_B T}{a^2N}\left[1+\exp \left(\frac{-\varepsilon}{k_B T}\right) \right].\]