Coupled lattice methods

In addition to the direct solution of PDEs corresponding to reaction-diffusion equations, in recent years attention has begun to be focused on the use of coupled lattice methods. In this approach, diffusion is not treated explicitly, but, rather, a lattice of elements in which the kinetic processes occur are coupled together in a variety of ways. The simulation of excitable media by cellular automata techniques has grown in popularity because they offer much greater computational efficiency for the two- and three-dimensional configurations required to study complex wave activity such as spirals and scroll waves. 

Computational techniques are centrally important at every stage of investigation of nonlinear dynamical systems. We have reviewed the main theoretical and computational tools used in studying these problems among these are bifurcation and stability analysis, numerical techniques for the solution of ordinary differential equations and partial differential equations, continuation methods, coupled lattice and cellular automata methods for the simulation of spatiotemporal phenomena, geometric representations of phase space attractors, and the numerical analysis of experimental data through the reconstruction of phase portraits, including the calculation of correlation dimensions and Lyapunov exponents from the data. 

These considerations suggest that lattice methods are somewhat more flexible and versatile for soft-matter simulations. On the other hand, the coupling between solvent and immersed particles is less straightforward than for a pure particle system. The coupling between solid particles and a lattice-based fluid model will be discussed in detail in Sect. 4. 

With the availabihty of computers, the transfer matrix method [14] emerged as an alternative and powerful technique for the study of cooperative phenomena of adsorbates resulting from interactions [15-17]. Quantities are calculated exactly on a semi-infinite lattice. Coupled with finite-size scaling towards the infinite lattice, the technique has proved popular for the determination of phase diagrams and critical-point properties of adsorbates [18-23] and magnetic spin systems [24—26], and further references therein. Application to other aspects of adsorbates, e.g., the calculation of desorption rates and heats of adsorption, has been more recent [27-30]. Sufficient accuracy can usually be obtained for the latter without scaling and essentially exact results are possible. In the following, we summarize the elementary but important aspects of the method to emphasize the ease of application. Further details can be found in the above references. 

To lessen experimental time, the null-point method may be employed by locating the pulse spacing, tnun, for which no magnetization is observed after the 180°-1-90° pulse-sequence. The relaxation rate is then obtained directly by using the relationship / , = 0.69/t n. In this way, a considerable diminution of measuring time is achieved, which is especially desirable in measurements of very low relaxation-rates, or for samples that are not very stable. In addition, estimates of relaxation rates for overlapping resonances can often be achieved. However, as the recovery curves for coupled spin-systems are, more often than not, nonexponential, observation of the null point may violate the initial-slope approximation. Hence, this method is best reserved for preliminary experiments that serve to establish the time scale for spin-lattice relaxation, and for qualitative conclusions. 

The essence of the ASAP method is based on the different spin-lattice relaxation behaviour of these two types of protons. In the case of a HMBC experiment, the acceptor protons are those directly bound to 13C with large spin-spin couplings they are the spins that give rise to the final spectrum. By contrast, the donor protons have negligible couplings to 13C and are therefore essentially unaffected by this polarization transfer sequence, which simply returns them to the z axis. In their method, Kupce and Freeman have proposed to replace the usual relaxation delay with a short cross-polarization (HOHAHA) interval. This offers a... 

This process involves the suspension of the biocatalyst in a monomer solution which is polymerized, and the enzymes are entrapped within the polymer lattice during the crosslinking process. This method differs from the covalent binding that the enzyme itself does not bind to the gel matrix. Due to the size of the biomolecule it will not diffuse out of the polymer network but small substrate or product molecules can transfer across or within it to ensure the continuous transformation. For sensing purposes, the polymer matrix can be formed directly on the surface of the fiber, or polymerized onto a transparent support (for instance, glass) that is then coupled to the fiber. The most popular matrices include polyacrylamide (Figure 5), silicone rubber, poly(vinyl alcohol), starch and polyurethane. 

Furthermore, the method of orientation selection can only be applied to systems with an electron spin-spin cross relaxation time Tx much larger than the electron spin-lattice relaxation time Tle77. In this case, energy exchange between the spin packets of the polycrystalline EPR spectrum by spin-spin interaction cannot take place. If on the other hand Tx < Tle, the spin packets are coupled by cross relaxation, and a powder-like ENDOR signal will be observed77. Since T 1 is normally the dominant relaxation rate in transition metal complexes, the orientation selection technique could widely be applied in polycrystalline and frozen solution samples of such systems (Sect. 6). 

It is often beneficial to define a coordinate Rti that describes the center of mass of the top layer. There are three common ways to set up the top layer. (1) The positions of top layer atoms r are confined to (lattice) sites rW)o, which are connected rigidly to the top layer. (2) The top layer atoms are coupled elastically to sites rra 0 fixed relative to the top layer, e.g., with springs of stiffness k. (3) An effective potential, such as a Steele potential Vs34 is applied between embedded atoms and the top layer. Specific advantages and disadvantages are associated with each method. Approach (1) may be the one that is most easily coded, (2) allows one to thermostat the outermost layer in an effective manner, whereas (3) is probably cheapest in terms of CPU time. 

While nonlinear in g2 (the coefficients 7,.r and b are lengthy integral expressions), the partial differential equation is linear in the derivatives. It can thus be solved by the method of characteristics, with the boundary conditions given by the coupling at /x = 0, as obtained from finite-T lattice QCD. 

---