Physics Stack Exchange is a question and answer site for active researchers, academics and students of physics. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

Let me use this lecture note as the reference.

  • I would like to know how in the above the expression (14) was obtained from expression (12).

In some sense it makes intuitive sense but I would like to know of the details of what happened in between those two equations. The point being that if there were no overall factor of "$\sqrt{N}$" in equation (12) then it would be a "textbook" case of doing the "method of steepest descent" in the asymptotic limit of "N".

  • I am wondering if in between there is an unwritten argument that in the "large N" limit one is absorbing the $\sqrt{N}$ into a re(un?)defined measure and then doing the steepest descent on only the exponential part of the integrand.

    I don't know how is the "method of steepest descent" supposed to be done on the entire integrand if the measure were not to be redefined.

  • But again if something like that is being done then why is there an approximation symbol in equation (14)?

After taking the thermodynamic limit and doing the steepest descent shouldn't the equation (14) become an equality with a sum over all the $\mu_s$ which solve equation (15)?

Though naively to me the expression (12) looks more amenable to a Dirac delta function interpretation since in the "large N" limit it seems to be looking like the standard representation of the Dirac delta function, $\frac{n}{\sqrt{\pi}} e^{-n^2x}$

  • I would like to know of some comments/explanation about this general philosophy/proof by which one wants to say that the "method of steepest descent" is exact in the "thermodynamic limit".
share|cite|improve this question

This is a straightforward application of the steepest descent method. So, we have the integral (taken from the lectures you cite)

$$Q_N=\sqrt{\frac{N\beta J}{2\pi}}\int_{-\infty}^\infty d\mu e^{[Nq(\beta J,\beta b,\mu)]}$$


$$q(\beta J,\beta b,\mu)=\ln\{2\cosh[\beta(J\mu+b)]\}-\frac{\beta J\mu^2}{2}.$$

Now, you will have, by applying steepest descent method to the integral,

$$\frac{\partial q}{\partial\mu}=\beta J\tanh[\beta(J\mu+b)]-\beta J\mu=0.$$

Let us call $\mu_s$ the solution of this equation and expand the argument of the exponential around this value. You will get

$$q(\beta J,\beta b,\mu)=q(\beta J,\beta b,\mu_s)-\frac{J\beta}{2}[1-J\beta(1-\mu_s^2)](\mu-\mu_s)^2.$$

This shows that

$$Q_N\approx e^{Nq(\beta J,\beta b,\mu_s)}\sqrt{\frac{N\beta J}{2\pi}}\int_{-\infty}^\infty d\mu e^{-\frac{NJ\beta}{2}[1-J\beta(1-\mu_s^2)](\mu-\mu_s)^2}$$


$$Q_N\approx e^{Nq(\beta J,\beta b,\mu_s)}\sqrt{\frac{1}{1-J\beta(1-\mu_s^2)}}$$

Note that the $N$ factor is completely removed after integration and remains just into the argument of the exponential. Eq. (16-18) follow straightforwardly.

share|cite|improve this answer
Thanks for the reply but I guess my question was not framed clearly. I CAN see what you have done BUT this is NOT what I would have thought the "method of steepest descent" to do! You seem to only do a Taylor expansion of the exponent. Point being that the only kind of steepest descent analysis that I have seen works on integrals of the form $F(s) = \int f(z)e^{sg(z)} dz$ in the large $s$ limit. This is not in that standard form. May be there is some other version of "method of steepest descent" that you have in mind and it would be great if you can give a reference for that. – user6818 Feb 13 '12 at 16:30
Also I want to know how this (or any valid steepest descent method) on this particular example will effectively manage to take the "large N" (thermodynamic) limit. – user6818 Feb 13 '12 at 16:32
@user6818 : this is standard textbook material. You should read, e.g., chapter VIII in this nice book. – Yvan Velenik Feb 13 '12 at 17:39
But THIS IS the steepest descent (it would be better to call it Laplace) method and the one that was applied in that lectures. Firstly, I find the extremum and then I expand around it and the integral becomes a Gaussian one and so computable. Just check – Jon Feb 13 '12 at 17:39
@user6818: The problem is that here you have no clear at all the mathematics. It is an everyday fact that the integral of an exponential like this takes its most important contribution where the argument of the exponential has an extremum. This is an asymptotic approximation to the integral and holds also for the oscillating case (stationary phase). I urge you to read some good book about asymptotic approximations like… . In this particular case it is just the integral that must be evaluated. – Jon Feb 14 '12 at 7:28

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.