Lattice formalism

The sub-lattice model is now the predominant model used in most CALPHAD calculations, whether it be to model an interstitial solid solution, an intermetallic compound such as 7-TiAl or an ionic solution. Numerous early papers, often centred around Fe-X-C systems, showed how the Hillert-Staffansson sub-lattice formalism (Hillert and Staffansson 1970) could be applied (see for example Lundberg et al. (1977) on Fe-Cr-C (Fig. 10.8) and Chatfteld and Hillert (1977) on Fe-Mo-C (Fig. 10.9)). Later work on systems such as Cr-Fe (Andersson and Sundman 1987) (Fig. 10.10) showed how a more generalised sub-lattice treatment developed by Sundman and Agren (1981) could be applied to multi-sub-lattice phases such as a. 

It should be noted in conclusion of this Section that preliminary results obtained by means of the discrete lattice formalism are presented in [85], This study demonstrates clearly the cooperative nature of the aggregation of two kinds of the Frenkel defects, vacancies and interstitials. 

The two processes that contribute to the photoemission current are the direct emission from the adsorbate orbital into a plane wave final state /) and the indirect emission from the adatom orbital via backscattering from the substrate lattice. Formally both effects are taken into account by replacing the matrix element of Eq. (14) by... 

But the concentrated phase is in equilibrium with a diluted phase which is not consistently described by a Flory-Huggins model it is difficult to imagine that a diluted solution of the two types of polyelectrolytes with an opposite global charge and a rather stiff conformation can be treated by the lattice formalism. 

Ewald s formalism reduces the infinite lattice sum to a serial complexity of in the number of particles n, which has been reduced to n logn in more recent formulations. A review of variants on Ewald summation methods which includes a more complete derivation of the basic method is in [3]. 

Jaquet and Miller [1985] have studied the transfer of hydrogen atom between neighbouring equilibrium positions on the (100) face of W by using a model two-dimensional chemosorption PES [McGreery and Wolken 1975]. In that calculation, performed for fairly high temperatures (T> rj the flux-flux formalism along with the vibrationally adiabatic approximation (section 3.6) were used. It has been noted that the increase of the coupling to the lattice vibrations and decrease of the frequency of the latter increase the transition probability. 

Here Tq are coordinates in a reference volume Vq and r = potential energy of Ar crystals has been computed [288] as well as lattice constants, thermal expansion coefficients, and isotope effects in other Lennard-Jones solids. In Fig. 4 we show the kinetic and potential energy of an Ar crystal in the canonical ensemble versus temperature for different values of P we note that in the classical hmit (P = 1) the low temperature specific heat does not decrease to zero however, with increasing P values the quantum limit is approached. In Fig. 5 the isotope effect on the lattice constant (at / = 0) in a Lennard-Jones system with parameters suitable for Ne atoms is presented, and a comparison with experimental data is made. Please note that in a classical system no isotope effect can be observed, x "" and the deviations between simulations and experiments are mainly caused by non-optimized potential parameters. 

It should be realized that unlike the study of equilibrium thermodynamics for which a model is often mapped onto Ising system, elementary mechanism of atomic motion plays a deterministic role in the kinetic study. In an actual alloy system, diffusion of an atomic species is mainly driven by vacancy mechanism. The incorporation of the vacancy mechanism into PPM formalism, however, is not readily achieved, since the abundant freedom of microscopic path of atomic movement demands intractable number of variational parameters. The present study is, therefore, limited to a simple spin kinetics, known as Glauber dynamics [14] for which flipping events at fixed lattice points drive the phase transition. Hence, the present study for a spin system is regarded as a precursor to an alloy kinetics. The limitation of the model is critically examined and pointed out in the subsequent sections. 

We shall extensively employ the notation of graph theory it provides a powerful and elegant formalism for the description of both the structure of the discrete lattices on which the CA live, and the complete dynamics (he. the global state transitions) induced by those structures. Graph theory also allows the correspondence between CA configurations and the words of a regular language to be made in a very natural fashion. 

Fig. 2.9 The state transition graph Gc, computed for a Tdim lattice consisting of iV = 4 points (with periodic boundary conditions), and totalistic rule T2 . The vertices labeled ti represent transient configurations those labeled cd represent cyclic states, and give rise to the formal cycle cum decomposition C[ j =[3, lj-f[2,2].

In this section we briefly outline a general mean-field theory approach to arbitrary PCA and then apply the formalism to a particular class of one-parameter rules. We then compare the theoretical predictions to numerical simulations on lattices of dimension 1 < d < 4. 

In order to write down the microscopic equations of motion more formally, we consider a size N x N 4-neighbor lattice with periodic boundary conditions. At each site (i, j) there are four cells, each of which is associated with one of the four neighbors of site (i,j). Each cell at time t can be in one of two states defined by a Boolean variable where d = 1,..., 4 labels, respectively, the east, north,... 

More formally, let us introduce the food -value of lattice site i,j) at time t, (= green or yellow) an occupancy variable, aij t), which is equal to zero unless the site i,j) contains at least one vant, in which case Qij t) = 1 and the vant vector = i,j,d), specifying the state of the vant at time t, where... 

The point of this terse introduction is that cellular automata represent not just a formalism for describing a certain particular class of behaviors (lattice gas simulations of fluid dynamics, models of chemical reactions and diffusion processes, etc.), but a much more general template for original and heretofore untapped ways of looking at a large class of unsolved or only poorly understood fundamental problems. 

All of these approaches, however, have at least this one basic point in common the only reason that a lattice structure is introduced is to formally prevent the appearance of ultraviolet divergences and thus simplify the otherwise difficult calculations that must be made on the continuum. In other words, it is the... 

Having thus established at least a formal equivalency between a discretized field theory on a lattice and CA, Svozil invokes the so-called no-go theorem to show that field theory cannot be discretized in this simple fashion. The no-go theorem (see [karstSl] and [nielSl]) states essentially that under a set of only mild assumptions it is impossible to formulate a local, unitary, charge conserving lattice held theory without effectively doubling the size of the predicted fermion population (i.e. species doubling see discussion box). 

Jourjine [jour85] generalizes Euclidean lattice field theory on a d-dimensional lattice to a cell complex. He uses homology theory to replace points by cells of various dimensions and fields by functions on cells, the cochains, in hopes of developing a formalism that treats space-time as a dynamical variable and describes the change in the dimension of space-time as a phase transition (see figure 12.19). 

The lowest coordination number of tantalum or niobium permitted by crystal chemistry formalism is 6, which corresponds to an octahedral configuration. X Me ratios that equal 3, 2 or 1 can, therefore, be obtained by corresponding substitutions in the cationic sub-lattice. A condition for such substitution is no doubt steric similarity between the second cation and the tantalum or niobium ion so as to enable its replacement in the octahedral polyhedron. In such cases, the structure of the compound consists of oxyfluoride octahedrons that are linked by their vertexes, sides or faces, according to the compound type, MeX3, MeX2 or MeX respectively. Table 37 lists compounds that have a coordination-type structure [259-261]. 

The thermodynamic transition between different forms as the above described is formally discontinuous. The difference between polymorphs is shown in general also by a different metrical description of the corresponding lattices. 

The requirement I > 2 can be understood from the symmetry considerations. The case of no restoring force, 1=1, corresponds to a domain translation. Within our picture, this mode corresponds to the tunneling transition itself. The translation of the defects center of mass violates momentum conservation and thus must be accompanied by absorbing a phonon. Such resonant processes couple linearly to the lattice strain and contribute the most to the phonon absorption at the low temperatures, dominated by one-phonon processes. On the other hand, I = 0 corresponds to a uniform dilation of the shell. This mode is formally related to the domain growth at T>Tg and is described by the theory in Xia and Wolynes [ 1 ]. It is thus possible, in principle, to interpret our formalism as a multipole expansion of the interaction of the domain with the rest of the sample. Harmonics with I > 2 correspond to pure shape modulations of the membrane. 

---