I struggle to understand the concept of free-energy landscape.

It seems to me the concept makes perfect sense for energies, but not for (canonical) free energies.

In my actual, hopefully to be bettered, understanding, a free energy landscape is built by fixing certain internal coordinates and computing the free energy under such constraints.

One should then be able to get information on the system behavior (for example, locating the equilibrium minima).

Let me clarify my doubts with a simple example: a particle in a harmonic well. There are two coordinates of interest: position $x$ (with related potential energy via a constant $\alpha$) and a fictitious internal coordinate $s$ (just for explanatory purposes, let us say the size of the particle): further we assume there is an energy related to the size, also quadratic w.r.t to the latter, via a constant $\gamma$.

For simplicity I completely neglect kinetic energy. So, the partition function for such a system would read

$$ Z = \int_0^{\infty} \int_{-\infty}^{\infty} \exp{[-\beta({\gamma s^2 + \alpha x^2})]} \mathrm{d}x\mathrm{d}s = \frac{1}{2} \sqrt {\frac{\pi}{\beta \alpha}} \sqrt{\frac{\pi}{\beta \gamma}}$$

From which the free energy at equilibrium follows. At equilibrium, both mean $x$ and $s$ will be different from zero.

Now I could be tempted to describe the free energy landscape by fixing the size $s$, and calculate the constrained free energy.

The “constrained” partition function reads $$ \bar{Z} (s) = \exp^{-\beta \gamma s^2} \sqrt {\frac{\pi}{\beta \alpha}} $$

The free energy landscape with respect to $s$ is a decreasing function of $s$, so one would conclude the system will try to minimize the size $s$ (to minimise its free energy), which is not correct.

I must be making quite some confusion between variables that can be constrained and variables for which constraining goes exactly against the meaning of canonical ensemble. Could somebody please help me to understand? Thanks in advance

LATER EDIT As a matter of fact, the more I think about I start to think the overall picture could be saved. Indeed the probability of having a size $s$ is proportional to $ \exp{-\beta \gamma s^2}$ but it is possibly incorrect to state this will result in the system selecting the smallest possible size $s$: indeed is the expected value that should matters. Any further hint would be most welcome.

  • $\begingroup$ I think you picked a poor example. Free energy landscapes are good for describing conformations of macromolecules (e.g., protein folding). $\endgroup$ Feb 16, 2018 at 12:19
  • 1
    $\begingroup$ I though I would pick an example in the simplest setting, just to expose "the flesh" of the matter. $\endgroup$
    – Smerdjakov
    Feb 16, 2018 at 12:38
  • 1
    $\begingroup$ @Smerdjakov "At equilibrium, both mean $x$ and $s$ will be different from zero"... It seems mean $x$ will be zero from the symmetry of the model $\endgroup$ Feb 16, 2018 at 15:28
  • $\begingroup$ @AlekseyDruggist, thanks for your comment. As a matter of fact you mention an aspect I might well completely misunderstand. For an oscillator in contact with thermal bath: indeed the symmetry of the potential dictates that $x=0$ is the equilibrium. Yet the average energy will be different from $0$, so does it make sense to talk about the equilibrium being $0$? $\endgroup$
    – Smerdjakov
    Feb 20, 2018 at 16:40
  • $\begingroup$ @Smerdjakov $\langle x^2\rangle \neq \langle x\rangle^2$. With your model, the average position (and size) is 0, even though their fluctuation is not. To have non-zero equilibrium size, you should use $(s-S)^2$, with $S=\langle s\rangle$. $\endgroup$
    – Adam
    Feb 22, 2018 at 7:52


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.