Patlak plot: Difference between revisions

From formulasearchengine
Jump to navigation Jump to search
No edit summary
 
en>Monkbot
 
Line 1: Line 1:
I'm Yoshiko Oquendo. Delaware is the only place I've been residing in. The occupation I've been occupying for many years is a bookkeeper but I've currently utilized for another one. To maintain birds is 1 of the issues he enjoys most.<br><br>Have a look at my weblog: [http://109.230.254.127/index.php?mod=users&action=view&id=29690 109.230.254.127]
A '''vertex model''' is a type of [[statistical mechanics]] [[Scientific modelling|model]] in which the [[Partition function (statistical mechanics)|Boltzmann weight]]s are associated with a [[vertex (graph theory)|vertex]] in the model (representing an [[atom]] or particle).<ref>R.J. Baxter, ''Exactly solved models in statistical mechanics'', London, Academic Press, 1982</ref><ref>V. Chari and A.N. Pressley, ''A Guide to Quantum Groups'' Cambridge University Press, 1994</ref> This contrasts with a nearest-neighbour model, such as the [[Ising model]], in which the energy, and thus the Boltzmann weight of a statistical microstate is attributed to the bonds connecting two neighbouring particles. The energy associated with a vertex in the lattice of particles is thus dependent on the state of the bonds which connect it to adjacent vertices. It turns out that every solution of the [[Yang-Baxter equation]] with spectral parameters in a tensor product of [[vector space]]s <math> V\otimes V </math> yields an exactly-solvable vertex model.
 
[[File:lattice2.jpg|450px|A 2-dimensional vertex model|thumb|right]]
Although the model can be applied to various [[geometry|geometries]] in any number of dimensions, with any number of possible states for a given bond, the most fundamental examples occur for two dimensional lattices, the simplest being a [[square lattice]] where each bond has two possible states. In this model, every particle is connected to four other particles, and each of the four bonds adjacent to the particle has two possible states, indicated by the direction of an arrow on the bond. In this model, each vertex can adopt <math>2^4</math> possible configurations. The [[energy]] for a given vertex can be given by <math>\varepsilon_{ij}^{k\ell}</math>,
[[File:vertex.jpg|100px|right|thumb|A vertex in the square lattice vertex model]] with a state of the lattice is an assignment of a state of each bond, with the total energy of the state being the sum of the vertex energies. As the energy is often divergent for an infinite lattice, the model is studied for a finite lattice as the lattice approaches infinite size. [[periodic function|Periodic]] or domain wall<ref>V.E. Korepin et al., ''Quantum inverse scattering method and correlation functions'', New York, Press Syndicate of the University of Cambridge, 1993</ref> [[boundary condition]]s may be imposed on the model.
 
==Discussion==
For a given state, the Boltzmann weight can be written in terms of the product of the Boltzmann weights of the corresponding vertices
 
:<math>\exp(-\beta \varepsilon(\mbox{state})) = \prod_\mbox{vertices} \exp(-\beta \varepsilon_{ij}^{k\ell})</math>
 
where the Boltzmann weights for the vertices are written
 
:<math>R_{ij}^{k\ell} = \exp(-\beta \varepsilon_{ij}^{k\ell})</math>.
 
The [[probability]] of the system being in any given state at a particular time, and hence the properties of the system are determined by the [[partition function (statistical mechanics)|partition function]], for which an analytic solution is desired.
 
:<math>\mathbb{Z} = \sum_\mbox{states} \exp(-\beta \varepsilon(\mbox{state})) </math>
 
where β=''1/kT'', ''T'' is [[temperature]] and ''k'' is [[Boltzmann's constant]]. The probability that the system is in any given [[microstate]] is given by
 
:<math>\frac{\exp(-\beta \varepsilon(\mbox{state}))}{\mathbb{Z}}</math>
 
so that the average value of the energy of the system is given by
 
:<math>
\langle \varepsilon \rangle = \frac{\sum_\mbox{states} \varepsilon \exp(-\beta \varepsilon)}{\sum_\mbox{states} \exp(-\beta \varepsilon)}
= kT^2 \frac{\partial}{\partial T} \ln \mathbb{Z}
</math>
 
In order to evaluate the partition function, firstly examine the states of a row of vertices.
 
[[File:vertexrow.jpg|400px|right|thumb|A row of vertices in the square lattice vertex model]]
 
The external edges are free variables, with summation over the internal bonds. Hence, form the row partition function
 
:<math>T_{i_1 k_1 \dots k_N}^{i'_1 \ell_1 \dots l_N} = \sum_{r_1,\dots,r_{N-1}}  R_{i_1 k_1}^{r_1 \ell_1} R_{r_1 k_2}^{r_2 \ell_2} \cdots R_{r_{N-1} k_N}^{i'_1 \ell_N}
</math>
 
This can be reformulated in terms of an auxiliary ''n''-dimensional vector space ''V'', with a [[basis]]{{disambiguation needed|date=April 2012}} <math>\{v_1, \ldots, v_n\}</math>, and <math>R \in End(V \otimes V)</math> as
 
:<math>R(v_i \otimes v_j) = \sum_{k,\ell} R_{ij}^{k\ell} v_k \otimes v_\ell </math>
 
and <math>T \in End(V \otimes V^{\otimes N})</math> as
 
:<math>T(v_{i_1} \otimes v_{k_1} \otimes \cdots \otimes v_{k_N}) = \sum_{i'_1,\ell_1, \dots \ell_N} T_{i_1 k_1 \dots k_N}^{i'_1 \ell_1 \dots \ell_N} v_{i'_1} \otimes v_{\ell_1} \otimes \cdots \otimes v_{\ell_N}</math>
 
thereby implying that ''T'' can be written as
:<math>T = R_{01}R_{02}\cdots R_{0N}</math>
 
where the indices indicate the factors of the [[tensor product]] <math> V \otimes V^{\otimes N}</math> on which ''R'' operates. Summing over the states of the bonds in the first row with the periodic boundary conditions <math>i_1 = i'_1</math>, gives
 
:<math>(\operatorname{trace}_{V}(T))_{k_1 \dots k_N }^{\ell_1 \dots \ell_N}</math>
 
where <math>\tau = \operatorname{trace}_{V}(T)</math> is the row-transfer matrix.
[[File:twovertexrow.jpg|400px|centre|thumb|Two rows of vertices in the square lattice vertex model]]
 
By summing the contributions over two rows, the result is
:<math>(\operatorname{trace}_{V}(T))_{k_1 \dots k_N }^{\ell_1 \dots \ell_N} (\operatorname{trace}_{V}(T))_{j_1 \dots j_N}^{k_1 \dots k_N} </math>
 
which upon summation over the vertical bonds connecting the first two rows gives:<math>((\operatorname{trace}_{V}(T))^2)_{j_1 \dots j_N }^{\ell_1 \dots \ell_N} </math>
 
for ''M'' rows, this gives
:<math>((\operatorname{trace}_{V}(T))^M)_{\ell'_1 \dots \ell'_N }^{\ell_1 \dots \ell_N} </math>
 
and then applying the periodic boundary conditions to the vertical columns, the partition function can be expressed in terms of the transfer matrix <math>tau</math> as
 
:<math>\mathbb{Z}= \operatorname{trace}_{V^{\otimes N}}(\tau^M)
\sim \lambda_{max}^M </math>
 
where <math>\lambda_{max}</math> is the largest [[eigenvalue]] of <math>\tau</math>. The approximation follows from the fact that the eigenvalues of <math>\tau^M</math> are the eigenvalues of <math>\tau</math> to the power of ''M'', and as <math>M \rightarrow \infty</math>, the power of the largest eigenvalue becomes much larger than the others. As the [[trace (linear algebra)|trace]] is the sum of the eigenvalues, the problem of calculating <math>\mathbb{Z}</math> reduces to the problem of finding the maximum eigenvalue of <math>\tau</math>. This in it itself is another field of study. However, a standard approach to the problem of finding the largest eigenvalue of <math>\tau</math> is to find a large family of operators which commute with <math>\tau</math>. This implies that the [[eigenspace]]s are common, and restricts the possible space of solutions. Such a family of commuting operators is usually found by means of the [[Yang-Baxter equation]], which thus relates statistical mechanics to the study of [[quantum group]]s.
 
== Integrability ==
'''Definition''': A vertex model is ''integrable'' if, <math>\forall \mu, \nu, \exists \lambda</math> such that
 
:<math> R_{12}(\lambda)R_{13}(\mu)R_{23}(\nu) = R_{23}(\nu)R_{13}(\mu)R_{12}(\lambda)</math>
 
This is a parameterized version of the Yang-Baxter equation, corresponding to the possible dependence of the vertex energies,and hence the Boltzmann weights ''R'' on external parameters, such as temperature, external fields, etc.
 
The integrability condition implies the following relation.
 
'''Proposition''': For an integrable vertex model, with <math>\lambda, \mu</math> and <math>\nu</math> defined as above, then
 
:<math>R(\lambda)(1 \otimes T(\mu))(T(\nu) \otimes 1) = (T(\nu) \otimes 1)(1 \otimes T(\mu))R(\lambda) </math>
 
as [[endomorphism]]s of <math>V \otimes V \otimes V^{\otimes N}</math>, where <math>R(\lambda)</math> acts on the first two vectors of the tensor product.
 
It follows by multiplying both sides of the above equation on the right by <math> R(\lambda)^{-1}</math> and using the cyclic property of the trace operator that the following corollary holds.
 
'''Corollary''': For an integrable vertex model for which <math>R(\lambda)</math> is invertible <math>\forall \lambda</math>, the transfer matrix <math>\tau(\mu)</math> commutes with <math>\tau(\nu), \ \forall \mu, \nu</math>.
 
This illustrates the role of the Yang-Baxter equation in the solution of solvable lattice models. Since the transfer matrices <math>\tau</math> commute for all <math>\lambda, \nu</math>, the eigenvectors of <math>\tau</math> are common, and hence independent of the parameterization. It is a recurring theme which appears in many other types of statistical mechanical models to look for these commuting transfer matrices.
 
From the definition of ''R'' above, it follows that for every solution of the Yang-Baxter equation in the tensor product of two ''n''-dimensional vector spaces, there is a corresponding 2-dimensional solvable vertex model where each of the bonds can be in the possible states <math>\{1,\ldots,n\}</math>, where ''R'' is an endomorphism in the space spanned by <math>\{|a \rangle \otimes |b \rangle\}, 1 \leq a,b \leq n </math>. This motivates the classification of all the finite-dimensional irreducible [[representation theory|representations]] of a given [[Quantum algebra]] in order to find solvable models coreesponding to it.
 
== Notable vertex models ==
* [[Six-vertex model]]
**[[File:Sixvertex.jpg|centre|200px|Configurations of the Six vertex model with Boltzmann weights]]
* [[Eight-vertex model]]
**[[File:Eightvertex.jpg|centre|200px|Configurations of the Eight vertex model with Boltzmann weights]]
*[http://projecteuclid.org/DPubS?service=UI&version=1.0&verb=Display&handle=euclid.cmp/1103909051 Nineteen-vertex model] (Izergin-Korepin model) <ref>[http://www.springerlink.com/content/v866837136883838/?p=8b77f6f1d1bf4e5aaed3e580361fd581&pi=0 A. G. Izergin and V. E. Korepin, The inverse scattering method approach to the quantum Shabat-Mikhailov model. ''Communications in Mathematical Physics'', '''79''', 303 (1981)]</ref>
 
== References ==
{{reflist}}
 
[[Category:Statistical mechanics]]
[[Category:Lattice models]]

Latest revision as of 06:57, 28 January 2014

A vertex model is a type of statistical mechanics model in which the Boltzmann weights are associated with a vertex in the model (representing an atom or particle).[1][2] This contrasts with a nearest-neighbour model, such as the Ising model, in which the energy, and thus the Boltzmann weight of a statistical microstate is attributed to the bonds connecting two neighbouring particles. The energy associated with a vertex in the lattice of particles is thus dependent on the state of the bonds which connect it to adjacent vertices. It turns out that every solution of the Yang-Baxter equation with spectral parameters in a tensor product of vector spaces VV yields an exactly-solvable vertex model.

A 2-dimensional vertex model

Although the model can be applied to various geometries in any number of dimensions, with any number of possible states for a given bond, the most fundamental examples occur for two dimensional lattices, the simplest being a square lattice where each bond has two possible states. In this model, every particle is connected to four other particles, and each of the four bonds adjacent to the particle has two possible states, indicated by the direction of an arrow on the bond. In this model, each vertex can adopt 24 possible configurations. The energy for a given vertex can be given by εijk,

A vertex in the square lattice vertex model

with a state of the lattice is an assignment of a state of each bond, with the total energy of the state being the sum of the vertex energies. As the energy is often divergent for an infinite lattice, the model is studied for a finite lattice as the lattice approaches infinite size. Periodic or domain wall[3] boundary conditions may be imposed on the model.

Discussion

For a given state, the Boltzmann weight can be written in terms of the product of the Boltzmann weights of the corresponding vertices

exp(βε(state))=verticesexp(βεijk)

where the Boltzmann weights for the vertices are written

Rijk=exp(βεijk).

The probability of the system being in any given state at a particular time, and hence the properties of the system are determined by the partition function, for which an analytic solution is desired.

=statesexp(βε(state))

where β=1/kT, T is temperature and k is Boltzmann's constant. The probability that the system is in any given microstate is given by

exp(βε(state))

so that the average value of the energy of the system is given by

ε=statesεexp(βε)statesexp(βε)=kT2Tln

In order to evaluate the partition function, firstly examine the states of a row of vertices.

A row of vertices in the square lattice vertex model

The external edges are free variables, with summation over the internal bonds. Hence, form the row partition function

Ti1k1kNi'11lN=r1,,rN1Ri1k1r11Rr1k2r22RrN1kNi'1N

This can be reformulated in terms of an auxiliary n-dimensional vector space V, with a basisTemplate:Disambiguation needed {v1,,vn}, and REnd(VV) as

R(vivj)=k,Rijkvkv

and TEnd(VVN) as

T(vi1vk1vkN)=i'1,1,NTi1k1kNi'11Nvi'1v1vN

thereby implying that T can be written as

T=R01R02R0N

where the indices indicate the factors of the tensor product VVN on which R operates. Summing over the states of the bonds in the first row with the periodic boundary conditions i1=i'1, gives

(traceV(T))k1kN1N

where τ=traceV(T) is the row-transfer matrix.

Two rows of vertices in the square lattice vertex model

By summing the contributions over two rows, the result is

(traceV(T))k1kN1N(traceV(T))j1jNk1kN

which upon summation over the vertical bonds connecting the first two rows gives:((traceV(T))2)j1jN1N

for M rows, this gives

((traceV(T))M)'1'N1N

and then applying the periodic boundary conditions to the vertical columns, the partition function can be expressed in terms of the transfer matrix tau as

=traceVN(τM)λmaxM

where λmax is the largest eigenvalue of τ. The approximation follows from the fact that the eigenvalues of τM are the eigenvalues of τ to the power of M, and as M, the power of the largest eigenvalue becomes much larger than the others. As the trace is the sum of the eigenvalues, the problem of calculating reduces to the problem of finding the maximum eigenvalue of τ. This in it itself is another field of study. However, a standard approach to the problem of finding the largest eigenvalue of τ is to find a large family of operators which commute with τ. This implies that the eigenspaces are common, and restricts the possible space of solutions. Such a family of commuting operators is usually found by means of the Yang-Baxter equation, which thus relates statistical mechanics to the study of quantum groups.

Integrability

Definition: A vertex model is integrable if, μ,ν,λ such that

R12(λ)R13(μ)R23(ν)=R23(ν)R13(μ)R12(λ)

This is a parameterized version of the Yang-Baxter equation, corresponding to the possible dependence of the vertex energies,and hence the Boltzmann weights R on external parameters, such as temperature, external fields, etc.

The integrability condition implies the following relation.

Proposition: For an integrable vertex model, with λ,μ and ν defined as above, then

R(λ)(1T(μ))(T(ν)1)=(T(ν)1)(1T(μ))R(λ)

as endomorphisms of VVVN, where R(λ) acts on the first two vectors of the tensor product.

It follows by multiplying both sides of the above equation on the right by R(λ)1 and using the cyclic property of the trace operator that the following corollary holds.

Corollary: For an integrable vertex model for which R(λ) is invertible λ, the transfer matrix τ(μ) commutes with τ(ν),μ,ν.

This illustrates the role of the Yang-Baxter equation in the solution of solvable lattice models. Since the transfer matrices τ commute for all λ,ν, the eigenvectors of τ are common, and hence independent of the parameterization. It is a recurring theme which appears in many other types of statistical mechanical models to look for these commuting transfer matrices.

From the definition of R above, it follows that for every solution of the Yang-Baxter equation in the tensor product of two n-dimensional vector spaces, there is a corresponding 2-dimensional solvable vertex model where each of the bonds can be in the possible states {1,,n}, where R is an endomorphism in the space spanned by {|a|b},1a,bn. This motivates the classification of all the finite-dimensional irreducible representations of a given Quantum algebra in order to find solvable models coreesponding to it.

Notable vertex models

References

43 year old Petroleum Engineer Harry from Deep River, usually spends time with hobbies and interests like renting movies, property developers in singapore new condominium and vehicle racing. Constantly enjoys going to destinations like Camino Real de Tierra Adentro.

  1. R.J. Baxter, Exactly solved models in statistical mechanics, London, Academic Press, 1982
  2. V. Chari and A.N. Pressley, A Guide to Quantum Groups Cambridge University Press, 1994
  3. V.E. Korepin et al., Quantum inverse scattering method and correlation functions, New York, Press Syndicate of the University of Cambridge, 1993
  4. A. G. Izergin and V. E. Korepin, The inverse scattering method approach to the quantum Shabat-Mikhailov model. Communications in Mathematical Physics, 79, 303 (1981)