वर्टेक्स मॉडल: Difference between revisions
No edit summary |
(→चर्चा) |
||
Line 1: | Line 1: | ||
'''वर्टेक्स मॉडल''' एक प्रकार का सांख्यिकीय यांत्रिकी मॉडल है जिसमें बोल्ट्ज़मान भार मॉडल में एक शीर्ष (परमाणु या कण का प्रतिनिधित्व) के साथ जुड़ा हुआ है।<ref>R.J. Baxter, ''Exactly solved models in statistical mechanics'', London, Academic Press, 1982</ref><ref>[[Vyjayanthi Chari|V. Chari]] and A.N. Pressley, ''A Guide to Quantum Groups'' Cambridge University Press, 1994</ref> यह निकटतम-पड़ोसी मॉडल, जैसे कि [[आइसिंग मॉडल]], के विपरीत है, जिसमें ऊर्जा, और इस प्रकार एक सांख्यिकीय माइक्रोस्टेट का बोल्ट्ज़मान वजन दो पड़ोसी कणों को जोड़ने वाले बांडों के लिए जिम्मेदार है। कणों की जाली में एक शीर्ष से जुड़ी ऊर्जा उन बंधनों की स्थिति पर निर्भर करती है जो इसे आसन्न शीर्षों से जोड़ते हैं। यह पता चला है कि वेक्टर रिक्त स्थान <math> V\otimes V </math> के टेंसर | '''वर्टेक्स मॉडल''' एक प्रकार का सांख्यिकीय यांत्रिकी मॉडल है जिसमें बोल्ट्ज़मान भार मॉडल में एक शीर्ष (परमाणु या कण का प्रतिनिधित्व) के साथ जुड़ा हुआ है।<ref>R.J. Baxter, ''Exactly solved models in statistical mechanics'', London, Academic Press, 1982</ref><ref>[[Vyjayanthi Chari|V. Chari]] and A.N. Pressley, ''A Guide to Quantum Groups'' Cambridge University Press, 1994</ref> यह निकटतम-पड़ोसी मॉडल, जैसे कि [[आइसिंग मॉडल]], के विपरीत है, जिसमें ऊर्जा, और इस प्रकार एक सांख्यिकीय माइक्रोस्टेट का बोल्ट्ज़मान वजन दो पड़ोसी कणों को जोड़ने वाले बांडों के लिए जिम्मेदार है। कणों की जाली में एक शीर्ष से जुड़ी ऊर्जा उन बंधनों की स्थिति पर निर्भर करती है जो इसे आसन्न शीर्षों से जोड़ते हैं। यह पता चला है कि वेक्टर रिक्त स्थान <math> V\otimes V </math> के टेंसर गुणनफल में वर्णक्रमीय मापदंडों के साथ यांग-बैक्सटर समीकरण का प्रत्येक समाधान एक बिल्कुल-हल करने योग्य वर्टेक्स मॉडल उत्पन्न करता है। | ||
[[File:lattice2.jpg|450px|एक 2-आयामी शीर्ष मॉडल|अंगूठे|दाएं]] | [[File:lattice2.jpg|450px|एक 2-आयामी शीर्ष मॉडल|अंगूठे|दाएं]] | ||
Line 8: | Line 8: | ||
==चर्चा== | ==चर्चा== | ||
जाली की दी गई स्थिति के लिए, | जाली की किसी दी गई स्थिति के लिए, बोल्ट्ज़मान भार को संबंधित शीर्ष अवस्थाओं के बोल्ट्ज़मान भार के शीर्षों पर गुणनफल के रूप में लिखा जा सकता है। | ||
: <math>\exp(-\beta \varepsilon(\mbox{state})) = \prod_\mbox{vertices} \exp(-\beta \varepsilon_{ij}^{k\ell})</math> | : <math>\exp(-\beta \varepsilon(\mbox{state})) = \prod_\mbox{vertices} \exp(-\beta \varepsilon_{ij}^{k\ell})</math> | ||
जहां शीर्षों के लिए बोल्ट्ज़मान भार | जहां शीर्षों के लिए बोल्ट्ज़मान भार लिखा हुआ है | ||
: <math>R_{ij}^{k\ell} = \exp(-\beta \varepsilon_{ij}^{k\ell})</math>, | : <math>R_{ij}^{k\ell} = \exp(-\beta \varepsilon_{ij}^{k\ell})</math>, | ||
और शीर्ष से जुड़े चार किनारों में से प्रत्येक की संभावित | और ''i, j, k, l'' शीर्ष से जुड़े चार किनारों में से प्रत्येक की संभावित स्थितियों पर आधारित है। आसन्न शीर्षों की शीर्ष स्थितियों को राज्य के स्वीकार्य होने के लिए कनेक्टिंग किनारों (बंधन) के साथ संगतता शर्तों को पूरा करना होगा। | ||
किसी विशेष समय पर सिस्टम के किसी भी दिए गए अवस्था में होने की संभावना, और इसलिए सिस्टम के गुण विभाजन फ़ंक्शन द्वारा निर्धारित किए जाते हैं, जिसके लिए एक विश्लेषणात्मक रूप वांछित है। | |||
: <math>\mathbb{Z} = \sum_\mbox{states} \exp(-\beta \varepsilon(\mbox{state})) </math> | : <math>\mathbb{Z} = \sum_\mbox{states} \exp(-\beta \varepsilon(\mbox{state})) </math> | ||
जहां β = 1/kT, T | जहां β = 1/kT, T तापमान है और k बोल्ट्ज़मान स्थिरांक है। सिस्टम के किसी निश्चित अवस्था ([[माइक्रोस्टेट]]) में होने की प्रायिकता निम्न द्वारा दी जाती है | ||
: <math>\frac{\exp(-\beta \varepsilon(\mbox{state}))}{\mathbb{Z}}</math> | : <math>\frac{\exp(-\beta \varepsilon(\mbox{state}))}{\mathbb{Z}}</math> | ||
ताकि सिस्टम की ऊर्जा का औसत | ताकि सिस्टम की ऊर्जा का औसत मान दिया जा सके | ||
: <math> | : <math> | ||
\langle \varepsilon \rangle = \frac{\sum_\mbox{states} \varepsilon \exp(-\beta \varepsilon)}{\sum_\mbox{states} \exp(-\beta \varepsilon)} | \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} | = kT^2 \frac{\partial}{\partial T} \ln \mathbb{Z} | ||
</math> | </math> | ||
विभाजन फ़ंक्शन का मूल्यांकन करने के लिए, | विभाजन फ़ंक्शन का मूल्यांकन करने के लिए, पहले शीर्षों की एक पंक्ति की स्थितियों की जाँच करें। | ||
[[File:vertexrow.jpg|400px|right|thumb|वर्गाकार जाली शीर्ष मॉडल में शीर्षों की एक पंक्ति]]बाहरी किनारे स्वतंत्र चर हैं, आंतरिक | [[File:vertexrow.jpg|400px|right|thumb|वर्गाकार जाली शीर्ष मॉडल में शीर्षों की एक पंक्ति]]बाहरी किनारे स्वतंत्र चर हैं, जिनमें आंतरिक बंधों का योग होता है। अतः, पंक्ति विभाजन फ़ंक्शन बनाएं | ||
: <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>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> | </math> | ||
इसे | इसे <math>R \in End(V \otimes V)</math>,<math>\{v_1, \ldots, v_n\}</math>के आधार पर सहायक ''n''-आयामी वेक्टर स्पेस ''V'' के संदर्भ में पुन: तैयार किया जा सकता है | ||
: <math>R(v_i \otimes v_j) = \sum_{k,\ell} R_{ij}^{k\ell} v_k \otimes v_\ell </math> | : <math>R(v_i \otimes v_j) = \sum_{k,\ell} R_{ij}^{k\ell} v_k \otimes v_\ell </math> | ||
और <math>T \in End(V \otimes V^{\otimes N})</math> जैसा | और <math>T \in End(V \otimes V^{\otimes N})</math> जैसा | ||
Line 34: | Line 34: | ||
इसका अर्थ यह है कि T को इस प्रकार लिखा जा सकता है | इसका अर्थ यह है कि T को इस प्रकार लिखा जा सकता है | ||
: <math>T = R_{0N}\cdots R_{02} R_{01} ,</math> | : <math>T = R_{0N}\cdots R_{02} R_{01} ,</math> | ||
जहां सूचकांक [[टेंसर उत्पाद]] के कारकों को दर्शाते हैं <math> V \otimes V^{\otimes N}</math> जिस पर R कार्य करता है। आवधिक सीमा शर्तों के साथ पहली पंक्ति में बांड की स्थिति का सारांश <math>i_1 = i'_1</math>, देता है | जहां सूचकांक [[टेंसर उत्पाद|टेंसर गुणनफल]] के कारकों को दर्शाते हैं <math> V \otimes V^{\otimes N}</math> जिस पर R कार्य करता है। आवधिक सीमा शर्तों के साथ पहली पंक्ति में बांड की स्थिति का सारांश <math>i_1 = i'_1</math>, देता है | ||
: <math>(\operatorname{trace}_{V}(T))_{k_1 \dots k_N }^{\ell_1 \dots \ell_N},</math> | : <math>(\operatorname{trace}_{V}(T))_{k_1 \dots k_N }^{\ell_1 \dots \ell_N},</math> | ||
कहाँ <math>\tau = \operatorname{trace}_{V}(T)</math> पंक्ति-स्थानांतरण मैट्रिक्स है। | कहाँ <math>\tau = \operatorname{trace}_{V}(T)</math> पंक्ति-स्थानांतरण मैट्रिक्स है। | ||
Line 56: | Line 56: | ||
'प्रस्ताव': एक पूर्णांक शीर्ष मॉडल के लिए, साथ <math>\lambda, \mu</math> और <math>\nu</math> फिर, ऊपर बताए अनुसार परिभाषित किया गया है | 'प्रस्ताव': एक पूर्णांक शीर्ष मॉडल के लिए, साथ <math>\lambda, \mu</math> और <math>\nu</math> फिर, ऊपर बताए अनुसार परिभाषित किया गया है | ||
: <math>R(\lambda)(1 \otimes T(\mu))(T(\nu) \otimes 1) = (T(\nu) \otimes 1)(1 \otimes T(\mu))R(\lambda) </math> | : <math>R(\lambda)(1 \otimes T(\mu))(T(\nu) \otimes 1) = (T(\nu) \otimes 1)(1 \otimes T(\mu))R(\lambda) </math> | ||
के [[एंडोमोर्फिज्म]] के रूप में <math>V \otimes V \otimes V^{\otimes N}</math>, कहाँ <math>R(\lambda)</math> टेंसर | के [[एंडोमोर्फिज्म]] के रूप में <math>V \otimes V \otimes V^{\otimes N}</math>, कहाँ <math>R(\lambda)</math> टेंसर गुणनफल के पहले दो वैक्टर पर कार्य करता है। | ||
इसके बाद उपरोक्त समीकरण के दोनों पक्षों को दाईं ओर से गुणा किया जाता है <math> R(\lambda)^{-1}</math> और ट्रेस ऑपरेटर की चक्रीय संपत्ति का उपयोग करना जो निम्नलिखित परिणाम रखता है। | इसके बाद उपरोक्त समीकरण के दोनों पक्षों को दाईं ओर से गुणा किया जाता है <math> R(\lambda)^{-1}</math> और ट्रेस ऑपरेटर की चक्रीय संपत्ति का उपयोग करना जो निम्नलिखित परिणाम रखता है। | ||
Line 64: | Line 64: | ||
यह सॉल्व करने योग्य जाली मॉडल के समाधान में यांग-बैक्सटर समीकरण की भूमिका को दर्शाता है। स्थानांतरण मैट्रिक्स के बाद से <math>\tau</math> सभी के लिए आवागमन <math>\lambda, \nu</math>, के eigenvectors <math>\tau</math> सामान्य हैं, और इसलिए पैरामीटरीकरण से स्वतंत्र हैं। यह एक आवर्ती विषय है जो इन कम्यूटिंग ट्रांसफर मैट्रिक्स को देखने के लिए कई अन्य प्रकार के सांख्यिकीय यांत्रिक मॉडल में दिखाई देता है। | यह सॉल्व करने योग्य जाली मॉडल के समाधान में यांग-बैक्सटर समीकरण की भूमिका को दर्शाता है। स्थानांतरण मैट्रिक्स के बाद से <math>\tau</math> सभी के लिए आवागमन <math>\lambda, \nu</math>, के eigenvectors <math>\tau</math> सामान्य हैं, और इसलिए पैरामीटरीकरण से स्वतंत्र हैं। यह एक आवर्ती विषय है जो इन कम्यूटिंग ट्रांसफर मैट्रिक्स को देखने के लिए कई अन्य प्रकार के सांख्यिकीय यांत्रिक मॉडल में दिखाई देता है। | ||
उपरोक्त आर की परिभाषा से, यह इस प्रकार है कि दो एन-आयामी वेक्टर रिक्त स्थान के टेंसर | उपरोक्त आर की परिभाषा से, यह इस प्रकार है कि दो एन-आयामी वेक्टर रिक्त स्थान के टेंसर गुणनफल में यांग-बैक्सटर समीकरण के प्रत्येक समाधान के लिए, एक संबंधित 2-आयामी सॉल्वेबल वर्टेक्स मॉडल होता है जहां प्रत्येक बॉन्ड हो सकता है संभावित अवस्थाएँ <math>\{1,\ldots,n\}</math>, जहां आर द्वारा फैलाए गए स्थान में एक एंडोमोर्फिज्म है <math>\{|a \rangle \otimes |b \rangle\}, 1 \leq a,b \leq n </math>. यह किसी दिए गए [[क्वांटम बीजगणित]] के सभी परिमित-आयामी अपरिवर्तनीय [[प्रतिनिधित्व सिद्धांत]] के वर्गीकरण को इसके अनुरूप हल करने योग्य मॉडल खोजने के लिए प्रेरित करता है। | ||
== उल्लेखनीय शीर्ष मॉडल == | == उल्लेखनीय शीर्ष मॉडल == |
Revision as of 09:02, 4 December 2023
वर्टेक्स मॉडल एक प्रकार का सांख्यिकीय यांत्रिकी मॉडल है जिसमें बोल्ट्ज़मान भार मॉडल में एक शीर्ष (परमाणु या कण का प्रतिनिधित्व) के साथ जुड़ा हुआ है।[1][2] यह निकटतम-पड़ोसी मॉडल, जैसे कि आइसिंग मॉडल, के विपरीत है, जिसमें ऊर्जा, और इस प्रकार एक सांख्यिकीय माइक्रोस्टेट का बोल्ट्ज़मान वजन दो पड़ोसी कणों को जोड़ने वाले बांडों के लिए जिम्मेदार है। कणों की जाली में एक शीर्ष से जुड़ी ऊर्जा उन बंधनों की स्थिति पर निर्भर करती है जो इसे आसन्न शीर्षों से जोड़ते हैं। यह पता चला है कि वेक्टर रिक्त स्थान के टेंसर गुणनफल में वर्णक्रमीय मापदंडों के साथ यांग-बैक्सटर समीकरण का प्रत्येक समाधान एक बिल्कुल-हल करने योग्य वर्टेक्स मॉडल उत्पन्न करता है।
यद्यपि मॉडल को किसी भी संख्या में आयामों में विभिन्न ज्यामिति पर लागू किया जा सकता है, किसी दिए गए बंधन के लिए संभावित अवस्थाओं की संख्या के साथ, सबसे मौलिक उदाहरण दो आयामी जाली के लिए होते हैं, सबसे सरल एक वर्ग जाली है जहां प्रत्येक बंधन में दो संभावित स्थितियां होती हैं। इस मॉडल में, प्रत्येक कण चार अन्य कणों से जुड़ा होता है, और कण से सटे चार बांडों में से प्रत्येक में दो संभावित अवस्थाएँ होती हैं, जो बांड पर एक तीर की दिशा से संकेतित होती हैं। इस मॉडल में, प्रत्येक शीर्ष संभावित विन्यास अपना सकता है। किसी दिए गए शिखर की ऊर्जा द्वारा दी जा सकती है,
जाली की स्थिति के साथ प्रत्येक बंधन की स्थिति का एक असाइनमेंट होता है, जिसमें अवस्था की कुल ऊर्जा शीर्ष ऊर्जाओं का योग होती है। चूंकि अनंत जाली के लिए ऊर्जा अक्सर अपसारी होती है, जैसे-जैसे जाली अनंत आकार के करीब पहुंचती है, मॉडल का अध्ययन एक सीमित जाली के लिए किया जाता है। मॉडल पर आवधिक या डोमेन दीवार [3] सीमा की शर्तें लगाई जा सकती हैं।
चर्चा
जाली की किसी दी गई स्थिति के लिए, बोल्ट्ज़मान भार को संबंधित शीर्ष अवस्थाओं के बोल्ट्ज़मान भार के शीर्षों पर गुणनफल के रूप में लिखा जा सकता है।
जहां शीर्षों के लिए बोल्ट्ज़मान भार लिखा हुआ है
- ,
और i, j, k, l शीर्ष से जुड़े चार किनारों में से प्रत्येक की संभावित स्थितियों पर आधारित है। आसन्न शीर्षों की शीर्ष स्थितियों को राज्य के स्वीकार्य होने के लिए कनेक्टिंग किनारों (बंधन) के साथ संगतता शर्तों को पूरा करना होगा।
किसी विशेष समय पर सिस्टम के किसी भी दिए गए अवस्था में होने की संभावना, और इसलिए सिस्टम के गुण विभाजन फ़ंक्शन द्वारा निर्धारित किए जाते हैं, जिसके लिए एक विश्लेषणात्मक रूप वांछित है।
जहां β = 1/kT, T तापमान है और k बोल्ट्ज़मान स्थिरांक है। सिस्टम के किसी निश्चित अवस्था (माइक्रोस्टेट) में होने की प्रायिकता निम्न द्वारा दी जाती है
ताकि सिस्टम की ऊर्जा का औसत मान दिया जा सके
विभाजन फ़ंक्शन का मूल्यांकन करने के लिए, पहले शीर्षों की एक पंक्ति की स्थितियों की जाँच करें।
बाहरी किनारे स्वतंत्र चर हैं, जिनमें आंतरिक बंधों का योग होता है। अतः, पंक्ति विभाजन फ़ंक्शन बनाएं
इसे ,के आधार पर सहायक n-आयामी वेक्टर स्पेस V के संदर्भ में पुन: तैयार किया जा सकता है
और जैसा
इसका अर्थ यह है कि T को इस प्रकार लिखा जा सकता है
जहां सूचकांक टेंसर गुणनफल के कारकों को दर्शाते हैं जिस पर R कार्य करता है। आवधिक सीमा शर्तों के साथ पहली पंक्ति में बांड की स्थिति का सारांश , देता है
कहाँ पंक्ति-स्थानांतरण मैट्रिक्स है। दो पंक्तियों में योगदान का योग करने पर, परिणाम मिलता है
जो पहली दो पंक्तियों को जोड़ने वाले ऊर्ध्वाधर बांडों पर योग करने पर देता है: एम पंक्तियों के लिए, यह देता है
और फिर ऊर्ध्वाधर स्तंभों पर आवधिक सीमा शर्तों को लागू करते हुए, विभाजन फ़ंक्शन को स्थानांतरण मैट्रिक्स के संदर्भ में व्यक्त किया जा सकता है जैसा
कहाँ का सबसे बड़ा eigenvalue है . अनुमान इस तथ्य से चलता है कि के eigenvalues के eigenvalues हैं एम की शक्ति के लिए, और के रूप में , सबसे बड़े eigenvalue की शक्ति दूसरों की तुलना में बहुत बड़ी हो जाती है। चूँकि ट्रेस (रैखिक बीजगणित) eigenvalues का योग है, गणना की समस्या का अधिकतम eigenvalue ज्ञात करने की समस्या कम हो जाती है . यह अपने आप में अध्ययन का एक और क्षेत्र है। हालाँकि, सबसे बड़ा eigenvalue खोजने की समस्या के लिए एक मानक दृष्टिकोण ऑपरेटरों के एक बड़े परिवार को ढूंढना है जो साथ यात्रा करते हैं . इसका तात्पर्य यह है कि eigenspace सामान्य हैं, और समाधान के संभावित स्थान को प्रतिबंधित करते हैं। कम्यूटिंग ऑपरेटरों का ऐसा परिवार आमतौर पर यांग-बैक्सटर समीकरण के माध्यम से पाया जाता है, जो इस प्रकार सांख्यिकीय यांत्रिकी को क्वांटम समूहों के अध्ययन से जोड़ता है।
अभिन्नता
परिभाषा: एक शीर्ष मॉडल अभिन्न है यदि, ऐसा है कि
यह यांग-बैक्सटर समीकरण का एक पैरामीटरयुक्त संस्करण है, जो शीर्ष ऊर्जाओं की संभावित निर्भरता के अनुरूप है, और इसलिए बोल्ट्ज़मैन तापमान, बाहरी क्षेत्रों आदि जैसे बाहरी मापदंडों पर आर को महत्व देता है।
अभिन्नता की स्थिति निम्नलिखित संबंध को दर्शाती है।
'प्रस्ताव': एक पूर्णांक शीर्ष मॉडल के लिए, साथ और फिर, ऊपर बताए अनुसार परिभाषित किया गया है
के एंडोमोर्फिज्म के रूप में , कहाँ टेंसर गुणनफल के पहले दो वैक्टर पर कार्य करता है।
इसके बाद उपरोक्त समीकरण के दोनों पक्षों को दाईं ओर से गुणा किया जाता है और ट्रेस ऑपरेटर की चक्रीय संपत्ति का उपयोग करना जो निम्नलिखित परिणाम रखता है।
परिणाम: एक पूर्णांक शीर्ष मॉडल के लिए जिसके लिए उलटा है , स्थानांतरण मैट्रिक्स के साथ आवागमन करता है .
यह सॉल्व करने योग्य जाली मॉडल के समाधान में यांग-बैक्सटर समीकरण की भूमिका को दर्शाता है। स्थानांतरण मैट्रिक्स के बाद से सभी के लिए आवागमन , के eigenvectors सामान्य हैं, और इसलिए पैरामीटरीकरण से स्वतंत्र हैं। यह एक आवर्ती विषय है जो इन कम्यूटिंग ट्रांसफर मैट्रिक्स को देखने के लिए कई अन्य प्रकार के सांख्यिकीय यांत्रिक मॉडल में दिखाई देता है।
उपरोक्त आर की परिभाषा से, यह इस प्रकार है कि दो एन-आयामी वेक्टर रिक्त स्थान के टेंसर गुणनफल में यांग-बैक्सटर समीकरण के प्रत्येक समाधान के लिए, एक संबंधित 2-आयामी सॉल्वेबल वर्टेक्स मॉडल होता है जहां प्रत्येक बॉन्ड हो सकता है संभावित अवस्थाएँ , जहां आर द्वारा फैलाए गए स्थान में एक एंडोमोर्फिज्म है . यह किसी दिए गए क्वांटम बीजगणित के सभी परिमित-आयामी अपरिवर्तनीय प्रतिनिधित्व सिद्धांत के वर्गीकरण को इसके अनुरूप हल करने योग्य मॉडल खोजने के लिए प्रेरित करता है।
उल्लेखनीय शीर्ष मॉडल
- छह-वर्टेक्स मॉडल
- * आठ-शीर्ष मॉडल
- * उन्नीस-वर्टेक्स मॉडल (इज़रगिन-कोरेपिन मॉडल) [4]
संदर्भ
- ↑ R.J. Baxter, Exactly solved models in statistical mechanics, London, Academic Press, 1982
- ↑ V. Chari and A.N. Pressley, A Guide to Quantum Groups Cambridge University Press, 1994
- ↑ V.E. Korepin et al., Quantum inverse scattering method and correlation functions, New York, Press Syndicate of the University of Cambridge, 1993
- ↑ 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)