next up previous contents
Next: Scaling picture of the Up: Random Manifolds Previous: Elastic manifolds

Phase diagram of elastic manifolds

 

The first thing we would like to know about the model is its phase diagram, at least schematically. Suppose we consider the system in the absence of disorder, i.e. for V=0. Then the Hamiltonian is Gaussian, and we can calculate any expectation value easily. To do this, we go to momentum space

equation158

The Hamiltonian becomes

equation164

where we inserted a high-momentum cutoff tex2html_wrap_inline2804 , which represents a sort of inverse minimum bending length for the manifold. We can now use the equipartition theorem to calculate

eqnarray171

Note that the correlations are power-law in momentum (and in real space for d<2). Such power-law correlations indicate scale invariance. In fact, this Gaussian theory is a trivial example of a RG fixed point. We can see this explicitly by performing a simple ``momentum shell'' RG. The idea is to separate out modes with high momentum and then integrate those out of the theory, thereby deriving an effective theory for the low-momentum modes. Here this is trivial, because the Hamiltonian has decoupled into a sum over momenta. In equations, we first split up the field

equation183

where

eqnarray188

Then we can write

equation204

where b>1 is called the rescaling factor. Integrating out the tex2html_wrap_inline2814 modes is then trivial, and we end up with just

equation217

where I have relabeled the field tex2html_wrap_inline2816 . Now this isn't quite the same form as we had before, because the cutoff has changed from tex2html_wrap_inline2804 to tex2html_wrap_inline2820 . But we can restore it by rescaling

eqnarray227

It's worth noting that in real space, these rescalings look a bit different, owing to the integrals defining Fourier transforms:

eqnarray236

After these rescalings, tex2html_wrap_inline2822 is restored to the original form of tex2html_wrap_inline2824 . We have therefore shown that the theory is at a fixed point.

I said earlier that the FP describing a phase should be absolutely stable. That means if we perturb the Hamiltonian slightly, it should be restored to the original form as the RG proceeds. When this is the case, we say that the perturbation is irrelevant. If, on the other hand, the perturbation actually grows under the RG, it is relevant, and destabilizes the FP. In the pure case, the allowed perturbations of tex2html_wrap_inline2824 are severely restricted by the requirement of translational invariance tex2html_wrap_inline2828 . This requires all additional terms in H to have many gradients, e.g.

equation250

Using just the real-space rescaling portion of the RG, we see that

eqnarray254

The mode integration has no effect on tex2html_wrap_inline2832 , because it too decouples into a mode sum. It can be shown that, although the mode integration for tex2html_wrap_inline2834 is nontrivial, it is actually justified to neglect its contribution to the renormalized Hamiltonian. Then the rescaling (or ``power counting'') result above is correct. The powers of b ( tex2html_wrap_inline2838 and tex2html_wrap_inline2840 ) are known as RG eigenvalues. Because they are negative, both perturbations are indeed irrelevant. This is true for any allowed perturbation because of the extra gradients.

So we have identified a fixed point and shown its stability in the pure system. Let's call this the thermal fixed point, because the fluctuations of the manifold are entirely due to thermal noise. Does the associated thermal phase survive in the presence of disorder? To answer this question, we should see if the random potential V is a relevant or irrelevant perturbation.

Let's try and answer this question at the same level we just did for the higher gradient perturbations, i.e. not worrying too much about the nasty mode integration. We need to consider the term

equation261

Under rescaling, this becomes

  equation265

It is not immediately apparent what to do with this, since the new rescaled Hamiltonian isn't just a constant multiplying the old Hamiltonian. The reason is that the disorder is really specified by a whole distribution, instead of a single parameter. We need to understand how the entire distribution rescales under the RG.

It IS possible to do this generally, but this turns out not to be so useful just now. Instead, let us restrict to the physical case of interest. For a narrowly distributed random potential, it is natural to describe the distribution in terms of its moments, the most important of which is just the variance. By translational invariance, this is a function

equation273

We would also like to consider short-range correlated potentials, in which tex2html_wrap_inline2844 decays rapidly as a function of both arguments. In this case we can make an ``ultra-local'' expansion of tex2html_wrap_inline2844 in terms of delta functions and their derivatives. This is equivalent to making a Taylor series expansion in the conjugate variables for the Fourier transform of tex2html_wrap_inline2844 . The leading term is then

equation283

Equation 32 tells us that the renormalized potential satisfies

equation290

Using the delta-function form for the correlations of V, we find

equation297

since tex2html_wrap_inline2852 . If we iterate this rescaling procedure, we get a recurrence relation. It looks particularly simple if we then take the limit of infinitesimal rescaling, so that tex2html_wrap_inline2854 . Then we have a differential ``flow'' equation

equation305

The criterion that g is irrelevant is thus

eqnarray310

PROBLEM: Show that when this is satisfied, all other random perturbations (i.e. higher moments and derivatives of delta-functions) are also irrelevant.

   figure315
Figure 3: N-d ``phase diagram'' for random manifolds.

This specifies a curve in the N-d plane (see Fig. 3). On one side the thermal phase is locally stable. We see that this never occurs for domain walls, which have N=1. For directed polymers, however, d=1 and the requirement reduces to N>2. Interestingly, for N=2 and d=1, i.e. the directed polymer in D=2+1=3 dimensions, tex2html_wrap_inline2876 . When the RG eigenvalue vanishes, we say that the perturbation is marginal.

What happens in this case? Well, so far we have discussed a purely linear theory of the RG. This is usually adequate to understand the local stability near a FP, but fails when one of the RG eigenvalues vanishes at that point. In this case, we need to develop the quadratic correction to the RG flow for g. This correction comes entirely from the mode elimination step of the RG, which I have essentially neglected up to now. It is in fact much more work to calculate this correction, and I won't go into the details of how you do so. Instead, let's see how far we can get just by general arguments.

The next order term in the renormalization of g would arise at second order in perturbation theory. Since we are only integrating out short wavelength modes, we should expect it to be analytic. In particular, there is no reason for the coefficient of tex2html_wrap_inline2882 to vanish or be in any way singular near the line where tex2html_wrap_inline2884 vanishes. Thus we expect an RG equation of the form

  equation326

The only important result that arises from a detailed calculation is that c>0. In fact, we could essentially rule out c<0 by various physical requirements, but I will just assume it and proceed.

Eq. 40 immediately answers our question. We see that for tex2html_wrap_inline2876 , disorder still grows under the rescaling. We say that it is marginally relevant. This means that the thermal phase is unstable for the 2+1-dimensional directed polymer.

But let's not stop there. There is more to be learned from Eq. 40. Consider small tex2html_wrap_inline2884 . Then there is a regime tex2html_wrap_inline2896 , when both the linear and quadratic terms are comparable, but much smaller than the neglected ( tex2html_wrap_inline2898 ) corrections. We then find that for tex2html_wrap_inline2900 , there is a new fixed point, for tex2html_wrap_inline2902 , which is unstable. This is then a critical fixed point, separating the thermal phase from some non-thermal phase, which occurs for stronger disorder. As tex2html_wrap_inline2904 , this fixed point moves into the origin, and the flows become unstable everywhere for tex2html_wrap_inline2906 .

The implication is that a disorder-dominated phase exists for any d and N. For some range of d and N, if the disorder is weak enough, the thermal phase continues to exist, but inevitably undergoes a transition to the disorder-dominated phase as the randomness is increased. In fact, the existence of a phase transition has been rigorously proven for N>2 and d=1, using the same kind of methods Marc described for the Sherrington-Kirkpatrick model[5]. By rough analogy to the bulk phases of the magnet we discussed earlier, the thermal fixed point describes the high-temperature phase of the system, where entropy and not energy is the dominant source of fluctuations. The strong disorder phase is presumably energy dominated, and we might guess that it is an ordered phase described by a zero-temperature fixed point. This is certainly consistent with the finding that g flows to strong coupling in this phase.

To really address this scenario, we need to consider the nature of the T=0 state of the random manifold. The stability of the zero temperature fixed point is addressed like any other fixed point, by determining the RG eigenvalue of temperature. It is conventional to write

equation344

where tex2html_wrap_inline2924 is the eigenvalue of temperature. For the ordered phase to be stable, we need tex2html_wrap_inline2926 .


next up previous contents
Next: Scaling picture of the Up: Random Manifolds Previous: Elastic manifolds

Leon Balents
Thu May 30 08:21:44 PDT 1996