We use cookies to distinguish you from other users and to provide you with a better experience on our websites. Close this message to accept cookies or find out how to manage your cookie settings.
To save content items to your account,
please confirm that you agree to abide by our usage policies.
If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your account.
Find out more about saving content to .
To save content items to your Kindle, first ensure [email protected]
is added to your Approved Personal Document E-mail List under your Personal Document Settings
on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part
of your Kindle email address below.
Find out more about saving to your Kindle.
Note you can select to save to either the @free.kindle.com or @kindle.com variations.
‘@free.kindle.com’ emails are free but can only be saved to your device when it is connected to wi-fi.
‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.
In the discussion in the previous chapter we neglected the fact that polymers are almost always immersed in a solvent. A good solvent is defined as a solvent in which it is energetically more favourable for the monomers of the polymer to be surrounded by molecules of the solvent than by other monomers. As a consequence, one can imagine that there exists round each monomer a region (the excluded volume) in which the chance of finding another monomer is very small. This will lead to a more open, more expanded structure for the polymer than if the excluded volume effects were absent.
The most popular model to describe this effect is the self avoiding walk. Here one considers only the subset of random walks which never visit the same site again. An example is given in figure 2.1. When one compares this figure with that of a random walk, the excluded volume effect is obvious.
Thus, the equilibrium properties of a polymer with excluded volume effects are studied by making averages over the set of all N-step self avoiding walks (SAW) (we will encounter a ‘continuum version’ of the SAW model in chapter 4). All energy effects are taken into account by limiting the set of allowed configurations to the self avoiding ones. For the moment all SAWs therefore have the same energy and thus when we calculate averages, we weight all configurations equally. Note that the self avoidance constraint doesn't come from the fact that no two monomers can be in the same place, as is often stated.
In the previous chapter we learned the importance of the O(n)-model for the study of polymers. In this chapter we will see how in two dimensions the critical behaviour of the O(n)-model has been determined exactly. The critical exponents of this model were first conjectured from renormalisation group arguments by Cardy and Hamber. These conjectures were then confirmed by Nienhuis using an approximate mapping onto the ‘Coulomb gas’. Since the Coulomb gas can be renormalised exactly, this lent support to the belief that the exponents of Cardy and Hamber were indeed exact. In 1986, Baxter succeeded in solving the O(n)-model exactly on the hexagonal lattice. His results, which were obtained in the thermodynamic limit, were later extended to finite systems by Batchelor and Blöte. In more recent years an O(n)-model on the square lattice has received a lot of attention since it has a very rich critical behaviour.
We conclude this chapter with a discussion of the SAW on fractal lattices.
The Coulomb gas approach to the SAW in d = 2
In this section we discuss how the O(n)-model (2.17) can be related to the Coulomb gas. This relation holds for the O(n)-model on the hexagonal lattice. It will give exact results for exponents ν and γ which because of universality should also hold on other lattices. Furthermore an exact value for μ follows from this mapping. When we talk about an exact solution we must note that the result is derived by nonrigorous means but that nevertheless it is generally believed that the result is the exact one. All numerical calculations furthermore support these conjectured values.
In this chapter we study self avoiding surfaces on a lattice. These surfaces are not immediately relevant for the study of polymers, although they could be of interest in the study of β-sheet polymers, which are important building blocks in proteins. Surfaces are used as models in the study of membranes or interfaces. Moreover, as we will see below, they allow the generalisation of vesicles to d = 3. The reason why we discuss surfaces in this book is mainly to show how the methods introduced in the statistical mechanics of polymers can be used in the study of other, but related, problems.
Several kinds of surfaces have been introduced in the literature. A distinction has to be made between surfaces which are models of polymerised membranes and those which describe liquid membranes. In the first case, the number of nearest neighbours of a given monomer is fixed. An interesting model is that of so called tethered surfaces introduced by Kantor, Kardar and Nelson. For liquid surfaces, on the other hand, the number of neighbours is not fixed. In this chapter, we will limit ourselves to a study of a lattice model of liquid surfaces, the ‘plaquette’ surfaces.
The critical behaviour of these surfaces is closely related to that of branched polymers. This is one of the many relations between surfaces and polymers.
Let us begin by defining the objects which we will study in this chapter and which will be referred to as plaquette surfaces or as self avoiding surfaces (SAS). An example of such a surface is shown in figure 11.1. The surface is built out of plaquettes of the cubic lattice.
Polymers are long chain molecules consisting of a large number of units (the monomers), which are held together by chemical bonds. These units may all be the same (in which case we speak of homopolymers) or may be different (heteropolymers).
Chemists spend most of their time developing polymers with specific chemical or physical properties. Such properties are often determined by the characteristics of the monomers and their mutual binding. In other words, they are determined on a local scale. In contrast, physicists work in the spirit of Richard Feynman and “have a habit of taking the simplest example of any phenomenon and calling it ‘physics’, leaving the more complicated examples to become the concern of other fields.” This attitude is taken to the extreme in the statistical mechanics of polymers, where one is interested mainly in universal properties, i.e. those properties that depend only on the fact that the polymer is a long linear molecule, and are determined by ‘large scale quantities’ such as the quality of the solvent in which the polymer is immersed, the temperature, the presence of surfaces (on which a polymer can adsorb) and so on.
Having this in mind, we can introduce a description of polymers in terms of random and self avoiding walks. When we look at the polymer on a microscopic scale we remember from our chemistry courses that one of the binding angles between successive monomers is essentially fixed (like the well known 105° angle between the two H–O bonds in a molecule of water), leaving one rotational degree of freedom (figure 1.1) for the chemical bond.
In d = 3 neither conformal invariance nor the Coulomb gas technique is very helpful in determining the critical behaviour of self avoiding walks. The transfer matrix can only reach up to small widths W, series can only get up to rather small N, and so on. So we have to look for different methods. There are essentially three of these. The first is probably the most obvious one; we can perform experiments on real polymers. In this book, we will only mention the results of these. Secondly, a powerful numerical method which so far has not been discussed is the Monte Carlo technique. It can of course be applied more easily and more accurately in d = 2, but in d = 3 it hats less competition from other methods. That's why we will discuss it in section 4.2 mainly from the point of view of learning about polymers in d = 3. We begin the chapter with a brief discussion of the third method, which is the RG approach to the critical behaviour of polymers.
Direct renormalisation of the Edwards model
We already encountered the exact RG methods for SAWs on fractal lattices in the previous chapter. But such a real space approach can only work very approximately on Euclidean lattices. The most precise RG calculations for polymers are performed with continuum techniques. A first method uses the O(n)-model, calculates the exponents of that model using techniques such as the ∈-expansion, and then sends n → 0 in the final equations. Here we will not get into these calculations; they have been very well described elsewhere.
Our discussion of lattice-gas models now takes a qualitative turn. We continue to study fluid mixtures as in the previous chapter, but now they will exhibit some surprising behavior—they won't like to mix!
This change in direction also steers us towards the heart of this book: models for complex hydrodynamics. The particular kind of complexity we introduce in this chapter relates to interfaces in immiscible fluids such as one might find in a mixture of oil and water. We are all familiar with the kind of bubbly complexity that that can entail. So it seems all the more remarkable that only a revised set of collision rules are needed to simulate it with lattice gases. Indeed, the models of immiscible fluids that we shall introduce are so close to the models of the previous chapters that we call them immiscible lattice gases.
This chapter, an introduction to immiscible lattice-gas mixtures, is limited to a discussion of two-dimensional models. In the next chapter, we introduce a lattice-Boltzmann method that is the “Boltzmann equivalent” of the immiscible lattice gas. That then sets the stage for our discussion of three-dimensional immiscible lattice gases in Chapter 11.
Color-dependent collisions
In the miscible lattice gases of the previous chapter, the collision rules were independent of color. The diffusive behavior derived instead from the redistribution of color after generic colorblind collisions were performed. Aside from some diffusion, the color simply went with the flow.
In this chapter we give a full derivation of the Navier-Stokes equation for the lattice gas. The first step, the Boltzmann approximation, is an approximation of the exact Liouville dynamics. The Boltzmann approximation should not be confused with the lattice-Boltzmann method of Chapter 6 and Chapter 7, but results that we obtain here for the Boltzmann approximation and the Navier-Stokes equations are also useful for the Boltzmann method. One of these results is the H-theorem for lattice gases. In this chapter we also reopen the tricky issue of the spurious invariants of lattice-gas or Boltzmann dynamics. We specifically discuss non-uniform global linear invariants for which, unlike the general nonlinear invariants, some theoretical results are known. Among them we find the staggered momentum invariants. We discuss their effect on hydrodynamics, which leads us to corrections to the Euler equation of Chapter 2.
General Boolean dynamics
It is useful to express the Boolean dynamics as a sequence of Boolean calculations, as we did in Section 2.6. In this section we shall denote by s or n local configurations. We further define a field of “rate bits”, which are random, Boolean variables, defined independently on each site and denoted by ass′(x, t). They are equal to one with probability 〈ass′〉 = A(s, s′). Thus if the pre-collision configuration is n the post-collision configuration is that single s′ for which ans′ = 1.
Our objectives for this chapter are twofold. First, we review some elementary aspects of fluid mechanics. We include in that discussion a classical derivation of the Navier-Stokes equations from the conservation of mass and momentum in a continuum fluid. We then discuss the analogous conservation relations in a lattice gas. Finally, we briefly describe the derivation of hydrodynamic equations for the lattice gas, but defer our first detailed discussion of this subject to the following chapter.
Molecular dynamics versus continuum mechanics
The study of fluids typically proceeds in either of two ways. Either one begins at the microscopic scale of molecular interactions, or one assumes that at a particular macroscopic scale a fluid may be described as a smoothly varying continuum. The latter approach allows us to write conservation equations in the form of partial-differential equations. Before we do so, however, it is worthwhile to recall the basis of such a point of view.
The macroscopic description of fluids corresponds to our everyday experience of flows. Figure 2.1 shows that a flow may have several characteristic length scales li. These lengths scales may be related either to geometric properties of the flow such as channel width or the diameter of obstacles or to intrinsic properties such as the size of vortical structures. The smallest of these length scales will be called Lhydro.
We consider a simple fluid such as air or water. This fluid is described by a number of thermodynamical fields, such as pressure, density, etc., as well as a velocity field u. A variable is called specific when it gives a quantity per unit mass. For instance, let E be the internal energy of a finite volume V of fluid. Let M = ρV be the mass of this volume of fluid. Then e = E/M is the specific internal energy. Table C.1 lists all the thermodynamic variables used in this appendix.
In addition to thermodynamic variables there are variables describing the external actions on the fluid. The effect of gravity, for instance may be represented by the acceleration f = g. The heating rate per unit mass q represents sources of heat, for instance from radiation. There may be heat and momentum exchanges inside the fluid, by heat conduction or viscous forces. These are not taken into account in the non-dissipative description.
For a one-component gas, there are only two independent thermodynamic variables. As a special choice, one may choose ρ, T as independent variables and express all quantities such as p, e, h, etc., in terms of ρ and T, but other choices are possible as well. These relations are linked through equations of state. For instance p = p(ρ, T) is the standard form for an equation of state.