arxiv: v1 [math.pr] 16 Jun 2010

Size: px
Start display at page:

Download "arxiv: v1 [math.pr] 16 Jun 2010"

Transcription

1 The GS inequality for the Potts Model arxiv: v1 [math.pr] 16 Jun 2010 Sérgio de Carvalho Bezerra affiliation: Universidade Federal de Pernambuco address: Av. Fernando Simões Barbosa Boa Viagem-Recife-PE-Brazil, CEP November 8, 2017 Abstract In this note we are concerned about the generalization of the GS inequality for the Potts model. We also obtain by a different method the proof of the GS inequality for the Ising model. We take advantage of a polynomial expansion and we derive that some of the variables of the polynomial can be separated. Keywords: Potts model, GS inequality and local magnetization. AMS:60K35, 82B99 1 Introduction In Statistical Mechanics the Ising model is one of the most important model. One possible extension of the Ising model is the Potts one. An excellent review is [1]. Recent works about the Potts model are [2],[3] and [4]. The Potts model can be described as follows: Consider N particles and for each particle i we associate a number i Σ {1,2,...,r} called the spin of the particle i. r is a positive integer bigger than one. The vector ( 1,..., N Σ N receives the name of the system configuration. In the nineteenth century people started using one function to describe a physical system: the energy function. In the Potts Model the energy function is: E( N J i,jδ( i, j B iδ(1, i (1 1 i<j N wherethenumbers(j i,j 1 i<j N,(B i 1 i N will beconsidernon-negative real numbers and the function δ(x,y 1 if x y and zero otherwise. The J i,j represents the exchange interaction between the particles i and j. As J i,j 0 we have only ferromagnetic interactions. The B i is the external field in the direction of the particle i. The probabilistic model is defined by the triple (Σ N,P(Σ N,P where i1 P( ee( Z N (2 1

2 thetermofnormalization Z N Σ N e E( isnamedpartitionfunction. We obtain the Ising model if we take r 2. Let J (J 1,2,...,J N 1,N and B (B 1,...,B N, we define the local magnetization as: m i(j,b < δ(1, i > Σ N δ(1, ip(. (3 It was proved in [5] which is called GS inequality that Theorem 1.1 Let r be equal to two. Then, for all i,j,k {1,...,N} m i(j,b B j B k 0. (4 Alternative proofs can be found in [6] or [7]. Naturally, we can ask for what does it happen with r bigger than two? This works gives an answer to this question. In fact, we prove: Theorem 1.2 (GS inequality for the potts modelgiven non-negative real numbers J and B the local magnetization has the following properties i if r 2 then ii if r 3 then 2 m i(j,b B j B k 0 (5 2 m i(j,b B j B k 0. (6 Item i is the theorem 1.1 which we obtain by a different method. It was obtained in [8] the proof of the concavity of the magnetization for a model of spins that are real-valued random variables. In the work [9] they proved the GS inequality for families of random variables which arise in certain ferromagnetic models of statistical mechanics and quantum field theory. In the next section, we obtain an expression for the second derivative of the local magnetization. In the following section, we determine that this derivative can be seen as a polynomial. After, we verify that this polynomial can have almost all the variables separated. In the last section, we evaluate some coefficient and prove the GS inequality for the Potts model. 2 The second derivative of the local magnetization The first step is to evaluate the second derivative of the local magnetization. 2

3 Proposition 2.1 Let r be a positive integer bigger than one. J and B vectors with entries non-negative. Then, we have 2 m i(j,b B j B k < δ(1, iδ(1, jδ(1, k > < δ(1, iδ(1, k >< δ(1, j > < δ(1, iδ(1, j >< δ(1, k > < δ(1, k δ(1, j >< δ(1, i > 2 < δ(1, i >< δ(1, k >< δ(1, j >. Proof if we take the first derivative we have: m i(j,b B k Σ N δ(1, iδ(1, k P( Now, we take the second derivative: 2 m i(j,b B j B k Σ N ˆ Σ N δ(1, iδ(1,ˆ k P(P(ˆ δ(1, iδ(1, k δ(1, jp( Σ N δ(1, iδ(1, k δ(1,ˆ jp(p(ˆ Σ N ˆ Σ N δ(1, iδ(1,ˆ k δ(1, jp(p(ˆ Σ N ˆ Σ N δ(1, iδ(1,ˆ k δ(1,ˆ jp(p(ˆ Σ N ˆ Σ N 2 δ(1, iδ(1,ˆ k δ(1,ˇ jp(pˆp(ˇ ˇ Σ N Σ N ˆ Σ N We pick up the idea introduced in [5] and we consider a ghost spin 0. Consequently, we consider the energy as E( 1 i<j N J i,jδ( i, j N B iδ( 0, i i1 0 i<j N J i,jδ( i, j where J 0,i B i and if we set up i 1,j 2,k 3 then the second derivative becomes: 2 m 1(J,B < δ( 0, 1δ( 0, 2δ( 0, 3 > B 2 B 3 < δ( 0, 1δ( 0, 3 >< δ( 0, 2 > < δ( 0, 1δ( 0, 2 >< δ( 0, 3 > < δ( 0, 3δ( 0, 2 >< δ( 0, 1 > 2 < δ( 0, 1 >< δ( 0, 2 >< δ( 0, 3 >. We will study the signal of this last expression. By one remark in [5] this signal implies the signal of the second derivative of the Proposition

4 3 The recursive property of the second derivative signal We define ( exp( 0 i<j N Thus, we pay attention for the following term: I(J,B: ZN 3 2 m 1(J,B B 2 B J i,jδ( i, j (7 Given a positive integer s let A s be the set of matrix A (a i,j with index 1 i ( 2 and 1 j 3 where ai,j {0,1} if the index i is bigger or equal to ( ( 2 s 1 and less or equal than 2.Further, a i,j 0 otherwise. We consider that the pairs of particles are in a specific order O {(0,1,(0,2,...,(N 1,N} {o 1,...o ( }. The row p in A is associated to the pair of particles o p.let s be s( exp( J i,jδ( i, j. For each matrix A (i,j the first ( 2 s terms in O A s we associate the term I A I 1,A I 2,A I 3,A I 4,A I 5,A with I 1,A ( s i j if o r (i,j,a r, i j if o r (i,j,a r,3 1 s s i j if o r (i,j,a r,

5 I 2,A ( I 3,A ( I 4,A ( I 5,A 2 ( i j if o r (i,j,a r, i j if o r (i,j,a r,3 1 i j if o r (i,j,a r, i j if o r (i,j,a r,3 1 i j if o r (i,j,a r, i j if o r (i,j,a r,3 1 s s s s s s 0 1 i j if o r (i,j,a r,1 1 s. 0 3 i j if o r (i,j,a r,3 1 s i j if o r (i,j,a r, i j if o r (i,j,a r, i j if o r (i,j,a r, i j if o r (i,j,a r,2 1 s s s s 5

6 It means that we have many constraints of the type i j which are defined by the element A. As For each row of A we have a pair o p (i,j we consider its weight defined by n(p: a(p,1 a(p,2 a(p,3 and also we designate by J p: J i,j. Now, we can announce the following result which can be interpreted as I(J,B is a polynomial in the variables X p: (e Jp 1 with,..., ( 2. Proposition 3.1 Let J and B be non-negative real numbers and s be a positive integer that belongs to {1,..., ( 2 }. Then I(J,B I n(( AX 2 s1...x n(( 2 (8 ( 2 s1 ( 2 A A s A particular case in which we are interested is when s ( 2 we have I(J,B I AX n(1 2 (9 A A ( X n(( ( 2 we remark that in this case ( 2 1. Proof The proof is by induction on s. Let s be equal to one. We observe that N 1 N N 1 N N 1 N e J ( 2 N 1 N J 1e ( 2 (e J ( 2 1 N 1 N. (10 In the second step we have just used the fact that N 1 N N 1 N. We use the notation [i 1,...,i n] if i1... in. Then, If we replace expression (10 for each one of the elements in the equation (7 we obtain: I(J,B ( ( (e J ( 2 1 (e J ( 2 1 6

7 ( [0,1,2,3] ( (e J ( 2 1 (e J ( 2 1 [0,1,2,3] ( [0,1,2] (e J ( 2 1 [0,1,2] ( [0,3] (e J ( 2 1 [0,3] ( (e J ( 2 1 ( [0,1,3] (e J ( 2 1 [0,1,3] ( [0,2] (e J ( 2 1 [0,2] ( (e J ( 2 1 ( [0,2,3] (e J ( 2 1 [0,2,3] 7

8 ( [0,1] 2 ( [0,1] ( [0,2] (e J ( 2 1 [0,1] (e J ( 2 1,[0,1] (e J ( 2 1 [0,2] ( [0,3] (e J ( 2 1 [0,3]. Thus, the result follows. Now, we consider that we proved the formula for a positive integer s. That means, I(J,B 2 s1 2. (11 A A s I AX n(( ( 2 s1 Let p ( N 2 s 1 with op (i,j. Then, s s1e Jpδ( i, j...x n(( ( 2 (12 we have that condition on the k s s condition on the k s i j condition on the k s i j condition on the k s i j condition on the k s s1e Jp s1 s1(e Jp 1 s1. (13 8

9 In the last step, we used the fact that s1 condition on the k s i j condition on the k s condition on the k s i j s1 s1. If we use expression (13 in the defintion of the term I A and we replace the result in the formula (11 we prove the proposition. 4 The separation formula At this moment, we obtain the separation formula of the variable X i s. Let O 1 O 2 O 3 where O 2 {(1,2,(1,3,(2,3}and the set O 3 be equal to {(0,1,(0,2,(0,3}. The pair of particles of O 2 will take a particular importance, we denote them by p 1 (1,2,p 2 (1,3,p 3 (2,3. Proposition 4.1 (Separation formula I(J,B ( 2 o p / O 1 ( 2 (1r 1 X p 3 (1(12r 1 X p (2r 1 r 2 X 2 p r 2 X 3 p 2 if o p / O 2 a(p,1 a(p,2 a(p,3 0 I AX n(p 1 p 1 X n(p 2 p 2 X n(p 3 p 3. Proof We start by the expression (9. Given A we define  as the matrix with entries equal to the entries of A in all rows except in row p. At row p we have â(p,1 â(p,2 â(p,3 0. We consider the case that p / O 1. Let n(p be the weight of the row p of A. We observe the following relations, if n(p 0 then I A IÂ. In the case that n( then I A r 1 IÂ. For n(p 2 we have I A r 2 IÂ. Finally, when n(p 3 we get I A r 3 IÂ. We do the calculations for n(. let p be the pair (i,j then n( means that appears the condition i j in one term of I 1,A,I 2,A,I 3,A,I 4,A and I 5,A defined in (8. Furthermore, this is the only additional constraint if we compare the conditions implied by A and Â. These constraints define a partition of the 9

10 elements {0,1,...,N} where two elements l,ˆl belong to the same subset if l ˆl. When we sum up some constraints of the type l ˆl we obtain r S where S is the number of subsets of the partition defined by the constraints. In the case that we are considering, we have one additional constraint if we compare A with Â. Thus, the partition implied by the constraints of A has the number of subsets of the partition implied by  less one. Further, as (i,j / O 1 the constraint i j is a new constraint. It is not an existent condition as in I 1,A. Thus, I A r 1 IÂ. We observe that there exists eight elements belong to A ( that we associate the same Â. They differ one from each other 2 by the possible values of a(p,1,a(p,2 and a(p,3. We have one element with n(p 0, three elements with n(, three elements with n(p 2 and one element with n(p 3. ence, the formula I(J,B 2 n(p 0 2 n( 2 n(p 2 2 n(p 3 I AX n(1 1...Xp 0...X n(( 1 2 ( 2 I AX n(1 1...Xp 1...X n(( 2 ( 2 I AX n(1 1...X 2 p...x n(( 2 ( 2 I AX n(1 1...Xp 3...X n(( 2 ( 2 10

11 can be simplified I(J,B 2 n(p 0 3r 1 2 n(p 0 3r 2 I AX n( ˆX p...x n(( 2 ( 2 I AX n( ˆX p...x n(( 2 n(p 0 I AX n( ˆX p...x n(( 2 n(p 0 2 ( 2 I AX n( ˆX p...x n(( 2 ( 2 2 X ( p 3 2 X p X 2 p where ˆX meansthatthistermdoesnotappear. Withmoresimplifications, we obtain I(J,B 2 n(p 0 I AX n( ˆX p...x n(( 2 ( 2 3r 2 Xp 2 r 3 Xp 3 I AX n( ˆX p...x n(( ( 2 2 n(p 0 2 We can repeat the procedure for all p / O 3 which produces (13r 1 X p (1r 1 X p 3. I(J,B ( 2 o p / O 1 (1r 1 X p 3 2 if o p / O 3 a(p,1 a(p,2 a(p,3 0 ( 2 I A X n(p p. (14 Now, we pay attention to the terms that belong to O 1. Let p 4 be (0,1. For a matrix with entries a(p,1 a(p,2 a(p,3 0 for 2 all p / O 3 we associate the matrix  with the same entries except for row 11

12 p 4 we have â(p,1 â(p,2 â(p,3 0. We can show that if n(p 4 0 then we have I A IÂ. Also, we obtain 2 n(p 4 1 n(p 0 p O 3 ( 2 I AX p4 (12r 1 X p4 For n(p 4 2 we have p p 4 2 n(p 4 0 n(p 0 p O 3 X n(p p ( 2 I A p p 4 X n(p p. (15 2 n(p 4 2 n(p 0 p O 3 ( 2 I AXp 2 4 (2r 1 r 2 X 2 p 4 p p 4 2 n(p 4 0 n(p 0 p O 3 X n(p p ( 2 I A p p 4 X n(p p. (16 And if n(p 4 3 we get I A r 2 IÂ. Again, we have eight matrix A which are associated for the same Â. We consider the three matrix with n(p 4 1, i.e., the matrix A 1 with a(p 4,1 1 the matrix A 2 with a(p 4,2 1 and A 3 with a(p 4,3 1. We have that I 1,A1 r 1 I 1, Â,I2,A 1 r 1 I 2, Â,I3,A 1 r 1 I 3, Â,I4,A 1 r 1 I 4, Â and I 5,A1 I 5, Â. In the last term does not appear the factor r 1 because the condition 0 1 appears independently of the matrix A 1. We also obtain that I 1,A2 r 1 I 1, Â,I2,A 2 I 2, Â,I3,A 2 I 3, Â,I4,A 2 r 1 I 4, Â and I 5,A2 r 1 I 5, Â. There are two terms that do not appear r 1 it happens because again the condition 0 1 occurs independently of the matrix A 2. For A 3 we have I 1,A3 I 1, Â,I2,A 3 r 1 I 2, Â,I3,A 3 r 1 I 3, Â,I4,A 3 I 4, Â and I5,A 3 r 1 I 5, Â. Again the terms that do not appear r 1 are consequence of the constraint 0 1 happens independently of the matrix A 3. Besides, equations (15 is a consequence of these facts. Now we look at the three matrices with n(p 2 let A 4 be the matrix with a(p 4,1 a(p 4,2 1 A 5 be the matrix with a(p 4,1 a(p 4,3 1 and A 5 the one with a(p 4,2 a(p 4,3 1. We obtain that I 1,A4 r 2 I 1, Â,I2,A 4 r 1 I 2, Â,I3,A 4 r 1 I 3, Â,I4,A 4 r 2 I 4, Â and I 5,A4 r 1 I 5, Â. The terms r 1 appear when one of the conditions 12

13 0 1 already occurs independently of A 4 and the term r 2 appears when both conditions 0 1 are new. For the matrix A 5 we obtain that I 1,A5 r 1 I 1, Â,I2,A 5 r 2 I 2, Â,I3,A 5 r 2 I 3, Â,I4,A 5 r 1 I 4, Â and I 5,A5 r 1 I 5, Â with the same reasoning. And if we see the matrix A6 we get I 1,A6 r 1 I 1, Â,I2,A 6 r 1 I 2, Â,I3,A 6 r 1 I 3, Â,I4,A 6 r 1 I 4, Â and I 5,A6 r 2 I 5, Â. Thus, we prove equation (16. As a matter of facts we obtain from equation (14 that I(J,B ( 2 o p / O 1 (1r 1 X p 3 2 n(p 4 0 if o p / O 3 n(p 0 ( 2 I A p p 4 X n(p p 2 n(p 4 1 if o p / O 3 n(p 0 ( 2 I AX p4 p p 4 X n(p p 2 n(p 4 2 if o p / O 3 n(p 0 ( 2 I AXp 2 4 p p 4 X n(p p 2 n(p 4 3 if o p / O 3 n(p 0 ( 2 I AXp 3 4 p p 4 X n(p p 13

14 By the remarks above and equation (15 and (16 we conclude that I(J,B ( 2 o p / O 1 (1r 1 X p 3 (1(12r 1 X p4 (2r 1 r 2 X 2 p 4 r 2 X 3 p 4 2 n(p 4 0 if o p / O 3 n(p 0 ( 2 I A p p 4 X n(p p. (17 We can repeat this last procedure for the other two elements of O 3 which ends the proof. 5 Some coefficients and the main theorem Now, We will evaluate the constants I A for 2 with n(p 0 if o p / O 2. These matrices can have entries different from zero only in three rows, we will represent these matrices by matrices 3 3 where the first row is associated to the pair (1,2 the second one to the pair (1,3 and the third one to (2,3. We also define the coefficients α(x,y,z 2 n(p 1 x,n(p 2 y n(p 3 z, for the other p n(p 0 Using the fact that if 1 2 and 1 3 implies that 2 3 and the symmetries, we get α(3,3,3 α(3,0,3 α(0,3,3 α(3,3,0 I ( r n 1 r n 1 r n 2 r n 1 r n 2 r n 2 r n 1 r n 2 r n 2 r n 1 r n 2 r n 2 2r n 2 r n 2 r n 2 r 3n 6 (r 2 3r 2. We observe by the last remark, we have α(3,2,3 α(2,3,3 α(3,3,2 α(3,1,3 α(1,3,3 α(3,3,1 I ( I ( r 3n 6 (r 2 3r 2 I ( I A 14

15 and α(3,2,2 α(2,3,2 α(2,2,3 I ( I ( I ( I ( I ( I ( I ( I ( I ( r 3n 5 ( 2r 1r 3n 6 (r 2 3r 2r 3n 6 (r 2 3r 2 r 3n 6 (r 2 3r 2r 3n 4 (r 1r 3n 6 (r 2 3r 2 r 3n 6 (r 2 3r 2r 3n 6 (r 2 3r 2r 3n 5 (r 2 3r 2 r 3n 6 (2r 3 15r 12. We also obtain α(3,2,1 α(2,3,1 α(3,1,2 α(1,3,2 α(1,2,3 α(2,1,3 I ( I ( I ( I ( I ( I ( I ( I ( I ( r 3n 5 ( 2r 1r 3n 5 ( 2r 1r 3n 6 (r 2 3r 2r 3n 4 (r 1 r 3n 6 (r 2 3r 2r 3n 4 (r 1r 3n 6 (r 2 3r 2 r 3n 5 (r 2 3r 2r 3n 5 (r 2 3r 2 r 3n 6 (4r 3 9r 2 3r 6. For the element α(3,2,0 α(2,3,0 α(3,0,2 α(0,3,2 α(0,2,3 α(2,0,3 I ( I ( I ( r 3n 5 ( 2r 1r 3n 4 (r 1r 3n 5 (r 2 3r 2 r 3n 5 (2r 2 6r 3. Now, the element α(3,1,1 α(1,3,1 α(1,1,3 I ( I ( I ( I ( I ( I ( I ( I ( I ( r 3n 5 ( 2r 1r 3n 4 (r 1r 3n 5 ( 2r 1r 3n 4 ( 2r2 r 3n 5 (r 2 3r 2r 3n 4 (r 1r 3n 5 (r 2 3r 2r 3n 3 (r 1 r 3n 5 (r 3 r 2 10r 6 and α(3,1,0 α(0,3,1 α(0,1,3 α(3,0,1 α(1,0,3 α(1,3,0 I ( I ( I ( r 3n 4 ( 2r 2r 3n 3 (r 1 r 3n 4 (r 2 3r 2 15

16 also α(3,0,0 α(0,3,0 α(0,0,3 I ( One of the longest terms is and α(2,2,2 I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( r 3n 4 ( 3r 3r 3n 5 ( 2r 2r 3n 5 ( 2r 2r 3n 5 ( 2r 2 r 3n 4 (r 1r 3n 6 (r 2 3r 2r 3n 5 ( 2r2r 3n 6 (r 2 3r 2 r 3n 5 (r 2 3r 2r 3n 5 ( 2r 2r 3n 4 (r 1r 3n 6 (r 2 3r 2 r 3n 4 (r 1r 3n 5 (r 3 3r 2r 3n 4 (r 1r 3n 6 (r 2 3r 2 r 3n 4 (r 1r 3n 5 (r 2 3r 2r 3n 5 ( 2r2 r 3n 6 (r 2 3r 2r 3n 5 (r 2 3r 2r 3n 6 (r 2 3r 2 r 3n 4 (r 1r 3n 5 (r 2 3r 2r 3n 5 (r 2 3r 2 r 3n 5 (r 2 3r 2r 3n 4 (r 2 3r 2 r 3n 6 (2r 4 6r 3 28r 2 8r 12 α(2,2,1 α(2,1,2 α(1,2,2 I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( r 3n 4 ( 3r 3r 3n 4 ( 3r 3r 3n 5 ( 2r20 I ( I ( I ( I ( I ( r 3n 5 ( 2r 2r 3n 4 (r 1r 3n 5 ( 2r 2r 3n 4 ( 2r 2 r 3n 5 (r 2 3r 20r 3n 5 ( 2r 2r 3n 4 (r 1 r 3n 5 (r 3 3r 2r 3n 4 (r 1r 3n 5 (r 3 3r 2r 3n 4 (r 1 r 3n 5 (r 2 3r 2r 3n 3 (r 1r 3n 5 ( 2r 2r 3n 4 ( 2r 2 16

17 r 3n 5 (r 2 3r 2r 3n 4 (r 1r 3n 5 (r 2 3r 2r 3n 3 (r 1 r 3n 5 (r 2 3r 2r 3n 4 (r 2 3r 2r 3n 4 (r 2 3r 2 r 3n 5 (5r 3 7r 2 22r 24. Now, the element α(2,2,0 α(2,0,2 α(0,2,2 I ( and I ( I ( I ( I ( I ( I ( I ( I ( r 3n 4 ( 3r 30r 3n 4 ( 2r 20r 3n 5 (r 3 3r 2 r 3n 3 (r 1r 3n 4 ( 2r 2r 3n 3 (r 1r 3n 4 (r 2 3r 2 r 3n 5 (4r 3 12r 2 6r 2 α(2,1,1 α(1,1,2 α(1,2,1 I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( r 3n 3 ( r1r 3n 4 ( 3r 30r 3n 4 ( 3r 3r 3n 3 ( 3r 3 r 3n 4 ( 2r 20r 3n 4 ( 2r 2r 3n 3 (r 1r 3n 3 (r 1 0r 3n 5 (r 3 3r 20r 3n 4 ( 2r 2r 3n 3 (r 1 r 3n 5 (r 3 3r 2r 3n 3 (r 1r 3n 3 (r 2 10r 3n 4 ( 2r 2 r 3n 3 (r 1r 3n 4 ( 2r 2r 3n 3 ( 2r 2r 3n 4 (r 2 3r 2 r 3n 3 (r 1r 3n 4 (r 2 3r 2r 3n 2 (r 1 r 3n 5 (2r 4 3r 3 23r 2 14r 4. We also obtain α(2,1,0 α(0,1,2 α(0,2,1 α(2,0,1 α(1,0,2 α(1,2,0 I ( I ( I ( I ( I ( I ( I ( I ( I ( r 3n 3 ( r1r 3n 3 ( 3r 30 r 3n 3 (r 10r 3n 3 (r 2 10r 3n 3 ( 2r 2r 3n 2 (r 1 r 3n 3 (2r 2 6r 4. 17

18 ere, we get α(2,0,0 α(0,0,2 α(2,0,0 I ( r 3n 2 ( r1r 3n 2 (r 10 0 and we can also get I ( I ( α(1,1,1 I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( I ( 0r 3n 3 ( r 1r 3n 3 (r 1r 3n 3 ( r 1r 3n 3 ( 3r 3 0r 3n 3 (r 10r 3n 3 (r 2 1r 3n 3 ( r1 r 3n 3 ( 3r 30r 3n 3 ( 3r 3r 3n 2 ( 3r 3r 3n 3 ( 2r 2 0r 3n 3 ( 2r2r 3n 2 (r 1r 3n 3 (r 10 r 3n 3 (r 2 10r 3n 3 ( 2r 2r 3n 2 (r 1r 3n 3 (r 2 1 r 3n 2 (r 1r 3n 2 (r 2 1 r 3n 3 (r 3 3r 2 16r 12. For the element α(1,1,0 α(1,0,1 α(0,1,1 I ( I ( I ( I ( I ( I ( I ( I ( I ( 0r 3n 2 ( r1r 3n 2 (r 1r 3n 2 ( r 1r 3n 2 ( 3r 3 0r 3n 2 (r 10r 3n 2 (r 2 1 r 3n 2 (r 2 3r 2 and α(1,0,0 α(0,0,1 α(0,1,0 I ( I ( 0r 3n 2 ( r 1r 3n 1 ( r1r 3n 1 (r 1 0. Finally, we have α(0,0,0 I ( 0 At this moment we can prove the theorem 1.2 : I ( 18

19 Proof (Theorem GS inequality for the Potts model The study of the signal of the quantity 2 m i (J,B B j B k is related to the signal of I(J,B (see equation (7. By the separation formula the signal of I(J, B is defined by the signal of the polynomial 2 if o p / O 2 a(p,1 a(p,2 a(p,3 0 n(p 1,n(p 2,n(p 3 I AX n(p 1 p 1 X n(p 2 p 2 X n(p 3 p 3 α(n(p 1,n(p 2,n(p 3X n(p 1 p 1 X n(p 2 p 2 X n(p 3 p 3 (18 But we can easily verify that the coefficients α(n(p 1,n(p 2,n(p 3 are all of them non positive if r 2 and are all of them non-negative if r 3. References [1] WuF., The Potts model. Reviewof modern Physics. Volume 54, Number 1, [2] Chang S. and Shrock R., Some exact results on the Potts model partition function in a magnetic field. J Phys. A: Math. Theor. Volume 42,Number 38, [3] Botirov G. and Rozikov U. Potts model with competing interactions on the Cayley tree: The contour method. Teoret. Mat. Fiz.,Volume 153, Number 1, 2007 pp [4] Chang S. and Shrock R., Structural properties of Potts model partition functions and chromatic polynomials for lattice strips. Physica A: Statistical Mechanics and its Applications. Volume 296, Number 1-2,2001, pp [5] Robert B. Griffths,C.A. urst,s. Sherman, Concavity of Magnetization of an Ising Ferromagnetic in a Positive External Field. Journal of Mathematical Physics. Volume 11, Number 3,1970, pp [6] Ellis R and Monroe J., A simple proof of the GS and further inequalities. Commun. math. Phys.. Volume 41,1975, pp [7] Lebowitz J., GS and other inequalities. Commun. math. Phys.. Volume 35, 1974, pp [8] Ellis R., Concavity of magnetization for a class of even ferromagnets. Bulletim of the American Mathematical Society. Volume 81, Number 5,1975, pp [9] Ellis R., Monroe J. and Newman C., The GS and other Correlation Inequalities for a class of even ferromagnets. Commun. math. Phys.. Volume 46, 1976, pp

PHYS Handout 6

PHYS Handout 6 PHYS 060 Handout 6 Handout Contents Golden Equations for Lectures 8 to Answers to examples on Handout 5 Tricks of the Quantum trade (revision hints) Golden Equations (Lectures 8 to ) ψ Â φ ψ (x)âφ(x)dxn

More information

Ising Model with Competing Interactions on Cayley Tree of Order 4: An Analytic Solution

Ising Model with Competing Interactions on Cayley Tree of Order 4: An Analytic Solution Journal of Physics: Conference Series OPEN ACCESS Ising Model with Competing Interactions on Cayley Tree of Order 4: An Analytic Solution To cite this article: Rukiah bin Ali et al 2013 J. Phys.: Conf.

More information

6.241 Dynamic Systems and Control

6.241 Dynamic Systems and Control 6.241 Dynamic Systems and Control Lecture 5: Matrix Perturbations Readings: DDV, Chapter 5 Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology February 16, 2011 E. Frazzoli

More information

Codes over Subfields. Chapter Basics

Codes over Subfields. Chapter Basics Chapter 7 Codes over Subfields In Chapter 6 we looked at various general methods for constructing new codes from old codes. Here we concentrate on two more specialized techniques that result from writing

More information

Characterizations of Strongly Regular Graphs: Part II: Bose-Mesner algebras of graphs. Sung Y. Song Iowa State University

Characterizations of Strongly Regular Graphs: Part II: Bose-Mesner algebras of graphs. Sung Y. Song Iowa State University Characterizations of Strongly Regular Graphs: Part II: Bose-Mesner algebras of graphs Sung Y. Song Iowa State University sysong@iastate.edu Notation K: one of the fields R or C X: a nonempty finite set

More information

j=1 [We will show that the triangle inequality holds for each p-norm in Chapter 3 Section 6.] The 1-norm is A F = tr(a H A).

j=1 [We will show that the triangle inequality holds for each p-norm in Chapter 3 Section 6.] The 1-norm is A F = tr(a H A). Math 344 Lecture #19 3.5 Normed Linear Spaces Definition 3.5.1. A seminorm on a vector space V over F is a map : V R that for all x, y V and for all α F satisfies (i) x 0 (positivity), (ii) αx = α x (scale

More information

ACO Comprehensive Exam March 17 and 18, Computability, Complexity and Algorithms

ACO Comprehensive Exam March 17 and 18, Computability, Complexity and Algorithms 1. Computability, Complexity and Algorithms (a) Let G(V, E) be an undirected unweighted graph. Let C V be a vertex cover of G. Argue that V \ C is an independent set of G. (b) Minimum cardinality vertex

More information

An Application of Catalan Numbers on Cayley Tree of Order 2: Single Polygon Counting

An Application of Catalan Numbers on Cayley Tree of Order 2: Single Polygon Counting BULLETIN of the Malaysian Mathematical Sciences Society http://math.usm.my/bulletin Bull. Malays. Math. Sci. Soc. (2) 31(2) (2008), 175 183 An Application of Catalan Numbers on Cayley Tree of Order 2:

More information

From Satisfiability to Linear Algebra

From Satisfiability to Linear Algebra From Satisfiability to Linear Algebra Fangzhen Lin Department of Computer Science Hong Kong University of Science and Technology Clear Water Bay, Kowloon, Hong Kong Technical Report August 2013 1 Introduction

More information

a 11 a 12 a 11 a 12 a 13 a 21 a 22 a 23 . a 31 a 32 a 33 a 12 a 21 a 23 a 31 a = = = = 12

a 11 a 12 a 11 a 12 a 13 a 21 a 22 a 23 . a 31 a 32 a 33 a 12 a 21 a 23 a 31 a = = = = 12 24 8 Matrices Determinant of 2 2 matrix Given a 2 2 matrix [ ] a a A = 2 a 2 a 22 the real number a a 22 a 2 a 2 is determinant and denoted by det(a) = a a 2 a 2 a 22 Example 8 Find determinant of 2 2

More information

arxiv:math/ v1 [math.pr] 10 Jun 2004

arxiv:math/ v1 [math.pr] 10 Jun 2004 Slab Percolation and Phase Transitions for the Ising Model arxiv:math/0406215v1 [math.pr] 10 Jun 2004 Emilio De Santis, Rossella Micieli Dipartimento di Matematica Guido Castelnuovo, Università di Roma

More information

arxiv: v1 [math.na] 22 Jun 2010

arxiv: v1 [math.na] 22 Jun 2010 On computing Bézier curves by Pascal matrix methods arxiv:1006.4327v1 [math.na] 22 Jun 2010 Licio Hernanes Bezerra a, Leonardo Koller Sacht b a Universidade Federal de Santa Catarina, Departamento de Matemática,

More information

Rank Analysis of Cubic Multivariate Cryptosystems

Rank Analysis of Cubic Multivariate Cryptosystems Rank Analysis of Cubic Multivariate Cryptosystems John Baena 1 Daniel Cabarcas 1 Daniel Escudero 2 Karan Khathuria 3 Javier Verbel 1 April 10, 2018 1 Universidad Nacional de Colombia, Colombia 2 Aarhus

More information

Comment on Conjectures on exact solution of three-dimensional (3D) simple orthorhombic Ising lattices

Comment on Conjectures on exact solution of three-dimensional (3D) simple orthorhombic Ising lattices arxiv:0811.1802v2 [cond-mat.stat-mech] 22 Nov 2008 Comment on Conectures on exact solution of three-dimensional (3D) simple orthorhombic Ising lattices Jacques H.H. Perk 145 Physical Sciences, Oklahoma

More information

1 T 1 = where 1 is the all-ones vector. For the upper bound, let v 1 be the eigenvector corresponding. u:(u,v) E v 1(u)

1 T 1 = where 1 is the all-ones vector. For the upper bound, let v 1 be the eigenvector corresponding. u:(u,v) E v 1(u) CME 305: Discrete Mathematics and Algorithms Instructor: Reza Zadeh (rezab@stanford.edu) Final Review Session 03/20/17 1. Let G = (V, E) be an unweighted, undirected graph. Let λ 1 be the maximum eigenvalue

More information

Legendre s Equation. PHYS Southern Illinois University. October 18, 2016

Legendre s Equation. PHYS Southern Illinois University. October 18, 2016 Legendre s Equation PHYS 500 - Southern Illinois University October 18, 2016 PHYS 500 - Southern Illinois University Legendre s Equation October 18, 2016 1 / 11 Legendre s Equation Recall We are trying

More information

POSITIVE AND NEGATIVE CORRELATIONS FOR CONDITIONAL ISING DISTRIBUTIONS

POSITIVE AND NEGATIVE CORRELATIONS FOR CONDITIONAL ISING DISTRIBUTIONS POSITIVE AND NEGATIVE CORRELATIONS FOR CONDITIONAL ISING DISTRIBUTIONS CAMILLO CAMMAROTA Abstract. In the Ising model at zero external field with ferromagnetic first neighbors interaction the Gibbs measure

More information

Travelling Salesman Problem

Travelling Salesman Problem Travelling Salesman Problem Fabio Furini November 10th, 2014 Travelling Salesman Problem 1 Outline 1 Traveling Salesman Problem Separation Travelling Salesman Problem 2 (Asymmetric) Traveling Salesman

More information

Math 343 Midterm I Fall 2006 sections 002 and 003 Instructor: Scott Glasgow

Math 343 Midterm I Fall 2006 sections 002 and 003 Instructor: Scott Glasgow Math 343 Midterm I Fall 006 sections 00 003 Instructor: Scott Glasgow 1 Assuming A B are invertible matrices of the same size, prove that ( ) 1 1 AB B A 1 = (11) B A 1 1 is the inverse of AB if only if

More information

Math 502 Fall 2005 Solutions to Homework 3

Math 502 Fall 2005 Solutions to Homework 3 Math 502 Fall 2005 Solutions to Homework 3 (1) As shown in class, the relative distance between adjacent binary floating points numbers is 2 1 t, where t is the number of digits in the mantissa. Since

More information

A necessary and sufficient condition for the existence of a spanning tree with specified vertices having large degrees

A necessary and sufficient condition for the existence of a spanning tree with specified vertices having large degrees A necessary and sufficient condition for the existence of a spanning tree with specified vertices having large degrees Yoshimi Egawa Department of Mathematical Information Science, Tokyo University of

More information

Dynamic Programming. Shuang Zhao. Microsoft Research Asia September 5, Dynamic Programming. Shuang Zhao. Outline. Introduction.

Dynamic Programming. Shuang Zhao. Microsoft Research Asia September 5, Dynamic Programming. Shuang Zhao. Outline. Introduction. Microsoft Research Asia September 5, 2005 1 2 3 4 Section I What is? Definition is a technique for efficiently recurrence computing by storing partial results. In this slides, I will NOT use too many formal

More information

Math , Fall 2012: HW 5 Solutions

Math , Fall 2012: HW 5 Solutions Math 230.0, Fall 202: HW 5 Solutions Due Thursday, October 4th, 202. Problem (p.58 #2). Let X and Y be the numbers obtained in two draws at random from a box containing four tickets labeled, 2, 3, 4. Display

More information

A tree-decomposed transfer matrix for computing exact partition functions for arbitrary graphs

A tree-decomposed transfer matrix for computing exact partition functions for arbitrary graphs A tree-decomposed transfer matrix for computing exact partition functions for arbitrary graphs Andrea Bedini 1 Jesper L. Jacobsen 2 1 MASCOS, The University of Melbourne, Melbourne 2 LPTENS, École Normale

More information

MATH 2030: MATRICES ,, a m1 a m2 a mn If the columns of A are the vectors a 1, a 2,...,a n ; A is represented as A 1. .

MATH 2030: MATRICES ,, a m1 a m2 a mn If the columns of A are the vectors a 1, a 2,...,a n ; A is represented as A 1. . MATH 030: MATRICES Matrix Operations We have seen how matrices and the operations on them originated from our study of linear equations In this chapter we study matrices explicitely Definition 01 A matrix

More information

arxiv: v3 [math.ra] 22 Aug 2014

arxiv: v3 [math.ra] 22 Aug 2014 arxiv:1407.0331v3 [math.ra] 22 Aug 2014 Positivity of Partitioned Hermitian Matrices with Unitarily Invariant Norms Abstract Chi-Kwong Li a, Fuzhen Zhang b a Department of Mathematics, College of William

More information

Pyramids and monomial blowing-ups

Pyramids and monomial blowing-ups Pyramids and monomial blowing-ups arxiv:math/0409446v1 [math.ac] 23 Sep 2004 M.J. Soto soto@us.es Departamento de Algebra Universidad de Sevilla October 4, 2018 José L. Vicente jlvc@us.es Departamento

More information

HOMEWORK 9 solutions

HOMEWORK 9 solutions Math 4377/6308 Advanced Linear Algebra I Dr. Vaughn Climenhaga, PGH 651A Fall 2013 HOMEWORK 9 solutions Due 4pm Wednesday, November 13. You will be graded not only on the correctness of your answers but

More information

Out-colourings of Digraphs

Out-colourings of Digraphs Out-colourings of Digraphs N. Alon J. Bang-Jensen S. Bessy July 13, 2017 Abstract We study vertex colourings of digraphs so that no out-neighbourhood is monochromatic and call such a colouring an out-colouring.

More information

Some Correlation Inequalities for Ising Antiferromagnets

Some Correlation Inequalities for Ising Antiferromagnets Some Correlation Inequalities for Ising Antiferromagnets David Klein Wei-Shih Yang Department of Mathematics Department of Mathematics California State University University of Colorado Northridge, California

More information

AB ExamSolutions Texas A&M High School Math Contest November 8, 2014

AB ExamSolutions Texas A&M High School Math Contest November 8, 2014 AB ExamSolutions Texas A&M High School Math Contest November 8, 2014 1. What is the largest power of 2 that divides 2 2013 + 10 2013? ANSWER: 2 2014 Solution: 2 2013 + 10 2013 = 2 2013 (1 + 5 2013 ). Since

More information

Andrew/CS ID: Midterm Solutions, Fall 2006

Andrew/CS ID: Midterm Solutions, Fall 2006 Name: Andrew/CS ID: 15-780 Midterm Solutions, Fall 2006 November 15, 2006 Place your name and your andrew/cs email address on the front page. The exam is open-book, open-notes, no electronics other than

More information

Integer Linear Programs

Integer Linear Programs Lecture 2: Review, Linear Programming Relaxations Today we will talk about expressing combinatorial problems as mathematical programs, specifically Integer Linear Programs (ILPs). We then see what happens

More information

Methods in Computer Vision: Introduction to Matrix Lie Groups

Methods in Computer Vision: Introduction to Matrix Lie Groups Methods in Computer Vision: Introduction to Matrix Lie Groups Oren Freifeld Computer Science, Ben-Gurion University June 14, 2017 June 14, 2017 1 / 46 Definition and Basic Properties Definition (Matrix

More information

Lecture 5: Probabilistic tools and Applications II

Lecture 5: Probabilistic tools and Applications II T-79.7003: Graphs and Networks Fall 2013 Lecture 5: Probabilistic tools and Applications II Lecturer: Charalampos E. Tsourakakis Oct. 11, 2013 5.1 Overview In the first part of today s lecture we will

More information

Math 489AB Exercises for Chapter 2 Fall Section 2.3

Math 489AB Exercises for Chapter 2 Fall Section 2.3 Math 489AB Exercises for Chapter 2 Fall 2008 Section 2.3 2.3.3. Let A M n (R). Then the eigenvalues of A are the roots of the characteristic polynomial p A (t). Since A is real, p A (t) is a polynomial

More information

Notes on the decomposition result of Karlin et al. [2] for the hierarchy of Lasserre by M. Laurent, December 13, 2012

Notes on the decomposition result of Karlin et al. [2] for the hierarchy of Lasserre by M. Laurent, December 13, 2012 Notes on the decomposition result of Karlin et al. [2] for the hierarchy of Lasserre by M. Laurent, December 13, 2012 We present the decomposition result of Karlin et al. [2] for the hierarchy of Lasserre

More information

MATH 323 Linear Algebra Lecture 6: Matrix algebra (continued). Determinants.

MATH 323 Linear Algebra Lecture 6: Matrix algebra (continued). Determinants. MATH 323 Linear Algebra Lecture 6: Matrix algebra (continued). Determinants. Elementary matrices Theorem 1 Any elementary row operation σ on matrices with n rows can be simulated as left multiplication

More information

1. Which one of the following points is a singular point of. f(x) = (x 1) 2/3? f(x) = 3x 3 4x 2 5x + 6? (C)

1. Which one of the following points is a singular point of. f(x) = (x 1) 2/3? f(x) = 3x 3 4x 2 5x + 6? (C) Math 1120 Calculus Test 3 November 4, 1 Name In the first 10 problems, each part counts 5 points (total 50 points) and the final three problems count 20 points each Multiple choice section Circle the correct

More information

CSC Discrete Math I, Spring Relations

CSC Discrete Math I, Spring Relations CSC 125 - Discrete Math I, Spring 2017 Relations Binary Relations Definition: A binary relation R from a set A to a set B is a subset of A B Note that a relation is more general than a function Example:

More information

On the calculation of inner products of Schur functions

On the calculation of inner products of Schur functions 172 J.A. Castilho Alcarás and V.K.B. Kota On the calculation of inner products of Schur functions J.A. Castilho Alcarás Instituto de Física Teórica, Universidade Estadual Paulista, UNESP, 01140-070, São

More information

MODEL ANSWERS TO THE FIRST QUIZ. 1. (18pts) (i) Give the definition of a m n matrix. A m n matrix with entries in a field F is a function

MODEL ANSWERS TO THE FIRST QUIZ. 1. (18pts) (i) Give the definition of a m n matrix. A m n matrix with entries in a field F is a function MODEL ANSWERS TO THE FIRST QUIZ 1. (18pts) (i) Give the definition of a m n matrix. A m n matrix with entries in a field F is a function A: I J F, where I is the set of integers between 1 and m and J is

More information

Math 203, Solution Set 4.

Math 203, Solution Set 4. Math 203, Solution Set 4. Problem 1. Let V be a finite dimensional vector space and let ω Λ 2 V be such that ω ω = 0. Show that ω = v w for some vectors v, w V. Answer: It is clear that if ω = v w then

More information

22m:033 Notes: 3.1 Introduction to Determinants

22m:033 Notes: 3.1 Introduction to Determinants 22m:033 Notes: 3. Introduction to Determinants Dennis Roseman University of Iowa Iowa City, IA http://www.math.uiowa.edu/ roseman October 27, 2009 When does a 2 2 matrix have an inverse? ( ) a a If A =

More information

2 Eigenvectors and Eigenvalues in abstract spaces.

2 Eigenvectors and Eigenvalues in abstract spaces. MA322 Sathaye Notes on Eigenvalues Spring 27 Introduction In these notes, we start with the definition of eigenvectors in abstract vector spaces and follow with the more common definition of eigenvectors

More information

arxiv:cond-mat/ v2 [cond-mat.stat-mech] 26 May 1998

arxiv:cond-mat/ v2 [cond-mat.stat-mech] 26 May 1998 Self-dual property of the Potts model in one dimension arxiv:cond-mat/9805301v2 [cond-mat.stat-mech] 26 May 1998 F. Y. Wu Department of Physics Northeastern University, Boston, Massachusetts 02115 Abstract

More information

MATRIX ALGEBRA AND SYSTEMS OF EQUATIONS. + + x 1 x 2. x n 8 (4) 3 4 2

MATRIX ALGEBRA AND SYSTEMS OF EQUATIONS. + + x 1 x 2. x n 8 (4) 3 4 2 MATRIX ALGEBRA AND SYSTEMS OF EQUATIONS SYSTEMS OF EQUATIONS AND MATRICES Representation of a linear system The general system of m equations in n unknowns can be written a x + a 2 x 2 + + a n x n b a

More information

Tutte Polynomials of Bracelets. Norman Biggs

Tutte Polynomials of Bracelets. Norman Biggs Tutte Polynomials of Bracelets Norman Biggs Department of Mathematics London School of Economics Houghton Street London WC2A 2AE U.K. n.l.biggs@lse.ac.uk January 2009 CDAM Research Reports Series LSE-CDAM

More information

arxiv: v2 [nlin.si] 27 Jun 2015

arxiv: v2 [nlin.si] 27 Jun 2015 On auto and hetero Bäcklund transformations for the Hénon-Heiles systems A. V. Tsiganov St.Petersburg State University, St.Petersburg, Russia e mail: andrey.tsiganov@gmail.com arxiv:1501.06695v2 [nlin.si]

More information

The Shortest Vector Problem (Lattice Reduction Algorithms)

The Shortest Vector Problem (Lattice Reduction Algorithms) The Shortest Vector Problem (Lattice Reduction Algorithms) Approximation Algorithms by V. Vazirani, Chapter 27 - Problem statement, general discussion - Lattices: brief introduction - The Gauss algorithm

More information

Total positivity in Markov structures

Total positivity in Markov structures 1 based on joint work with Shaun Fallat, Kayvan Sadeghi, Caroline Uhler, Nanny Wermuth, and Piotr Zwiernik (arxiv:1510.01290) Faculty of Science Total positivity in Markov structures Steffen Lauritzen

More information

arxiv: v1 [math.na] 30 Jun 2011

arxiv: v1 [math.na] 30 Jun 2011 arxiv:11066263v1 [mathna] 30 Jun 2011 Another proof of Pell identities by using the determinant of tridiagonal matrix Meral Yaşar &Durmuş Bozkurt Department of Mathematics, Nigde University and Department

More information

From Potts to Tutte and back again... A graph theoretical view of statistical mechanics

From Potts to Tutte and back again... A graph theoretical view of statistical mechanics From Potts to Tutte and back again... A graph theoretical view of statistical mechanics Jo Ellis-Monaghan e-mail: jellis-monaghan@smcvt.edu website: http://academics.smcvt.edu/jellis-monaghan 10/27/05

More information

1 Last time: least-squares problems

1 Last time: least-squares problems MATH Linear algebra (Fall 07) Lecture Last time: least-squares problems Definition. If A is an m n matrix and b R m, then a least-squares solution to the linear system Ax = b is a vector x R n such that

More information

Math/CS 466/666: Homework Solutions for Chapter 3

Math/CS 466/666: Homework Solutions for Chapter 3 Math/CS 466/666: Homework Solutions for Chapter 3 31 Can all matrices A R n n be factored A LU? Why or why not? Consider the matrix A ] 0 1 1 0 Claim that this matrix can not be factored A LU For contradiction,

More information

Lecture 12 (Tue, Mar 5) Gaussian elimination and LU factorization (II)

Lecture 12 (Tue, Mar 5) Gaussian elimination and LU factorization (II) Math 59 Lecture 2 (Tue Mar 5) Gaussian elimination and LU factorization (II) 2 Gaussian elimination - LU factorization For a general n n matrix A the Gaussian elimination produces an LU factorization if

More information

Week 15-16: Combinatorial Design

Week 15-16: Combinatorial Design Week 15-16: Combinatorial Design May 8, 2017 A combinatorial design, or simply a design, is an arrangement of the objects of a set into subsets satisfying certain prescribed properties. The area of combinatorial

More information

Lee Yang zeros and the Ising model on the Sierpinski gasket

Lee Yang zeros and the Ising model on the Sierpinski gasket J. Phys. A: Math. Gen. 32 (999) 57 527. Printed in the UK PII: S35-447(99)2539- Lee Yang zeros and the Ising model on the Sierpinski gasket Raffaella Burioni, Davide Cassi and Luca Donetti Istituto Nazionale

More information

Statistical Mechanics

Statistical Mechanics Franz Schwabl Statistical Mechanics Translated by William Brewer Second Edition With 202 Figures, 26 Tables, and 195 Problems 4u Springer Table of Contents 1. Basic Principles 1 1.1 Introduction 1 1.2

More information

Lecture 2: Lattices and Bases

Lecture 2: Lattices and Bases CSE 206A: Lattice Algorithms and Applications Spring 2007 Lecture 2: Lattices and Bases Lecturer: Daniele Micciancio Scribe: Daniele Micciancio Motivated by the many applications described in the first

More information

Maximum union-free subfamilies

Maximum union-free subfamilies Maximum union-free subfamilies Jacob Fox Choongbum Lee Benny Sudakov Abstract An old problem of Moser asks: how large of a union-free subfamily does every family of m sets have? A family of sets is called

More information

Chapter 2: Linear Independence and Bases

Chapter 2: Linear Independence and Bases MATH20300: Linear Algebra 2 (2016 Chapter 2: Linear Independence and Bases 1 Linear Combinations and Spans Example 11 Consider the vector v (1, 1 R 2 What is the smallest subspace of (the real vector space

More information

Linear Systems and Matrices

Linear Systems and Matrices Department of Mathematics The Chinese University of Hong Kong 1 System of m linear equations in n unknowns (linear system) a 11 x 1 + a 12 x 2 + + a 1n x n = b 1 a 21 x 1 + a 22 x 2 + + a 2n x n = b 2.......

More information

What is A + B? What is A B? What is AB? What is BA? What is A 2? and B = QUESTION 2. What is the reduced row echelon matrix of A =

What is A + B? What is A B? What is AB? What is BA? What is A 2? and B = QUESTION 2. What is the reduced row echelon matrix of A = STUDENT S COMPANIONS IN BASIC MATH: THE ELEVENTH Matrix Reloaded by Block Buster Presumably you know the first part of matrix story, including its basic operations (addition and multiplication) and row

More information

Brazilian Journal of Physics ISSN: Sociedade Brasileira de Física Brasil

Brazilian Journal of Physics ISSN: Sociedade Brasileira de Física Brasil Brazilian Journal of Physics ISSN: 0103-9733 luizno.bjp@gmail.com Sociedade Brasileira de Física Brasil Braga, Gastão A.; Sá Barreto, F. C.; Santos, Jander P. Upper Bounds on the Critical Temperature of

More information

Decay of correlations in 2d quantum systems

Decay of correlations in 2d quantum systems Decay of correlations in 2d quantum systems Costanza Benassi University of Warwick Quantissima in the Serenissima II, 25th August 2017 Costanza Benassi (University of Warwick) Decay of correlations in

More information

CSC 2414 Lattices in Computer Science September 27, Lecture 4. An Efficient Algorithm for Integer Programming in constant dimensions

CSC 2414 Lattices in Computer Science September 27, Lecture 4. An Efficient Algorithm for Integer Programming in constant dimensions CSC 2414 Lattices in Computer Science September 27, 2011 Lecture 4 Lecturer: Vinod Vaikuntanathan Scribe: Wesley George Topics covered this lecture: SV P CV P Approximating CVP: Babai s Nearest Plane Algorithm

More information

MATH STUDENT BOOK. 12th Grade Unit 9

MATH STUDENT BOOK. 12th Grade Unit 9 MATH STUDENT BOOK 12th Grade Unit 9 Unit 9 COUNTING PRINCIPLES MATH 1209 COUNTING PRINCIPLES INTRODUCTION 1. PROBABILITY DEFINITIONS, SAMPLE SPACES, AND PROBABILITY ADDITION OF PROBABILITIES 11 MULTIPLICATION

More information

Summary: A Random Walks View of Spectral Segmentation, by Marina Meila (University of Washington) and Jianbo Shi (Carnegie Mellon University)

Summary: A Random Walks View of Spectral Segmentation, by Marina Meila (University of Washington) and Jianbo Shi (Carnegie Mellon University) Summary: A Random Walks View of Spectral Segmentation, by Marina Meila (University of Washington) and Jianbo Shi (Carnegie Mellon University) The authors explain how the NCut algorithm for graph bisection

More information

Math 315: Linear Algebra Solutions to Assignment 7

Math 315: Linear Algebra Solutions to Assignment 7 Math 5: Linear Algebra s to Assignment 7 # Find the eigenvalues of the following matrices. (a.) 4 0 0 0 (b.) 0 0 9 5 4. (a.) The characteristic polynomial det(λi A) = (λ )(λ )(λ ), so the eigenvalues are

More information

Matrices: 2.1 Operations with Matrices

Matrices: 2.1 Operations with Matrices Goals In this chapter and section we study matrix operations: Define matrix addition Define multiplication of matrix by a scalar, to be called scalar multiplication. Define multiplication of two matrices,

More information

arxiv: v6 [math.nt] 16 May 2017

arxiv: v6 [math.nt] 16 May 2017 arxiv:5.038v6 [math.nt] 6 May 207 An alternative proof of the Dirichlet prime number theorem Haifeng Xu May 7, 207 Abstract Dirichlet s theorem on arithmetic progressions called as Dirichlet prime number

More information

1.8 Dual Spaces (non-examinable)

1.8 Dual Spaces (non-examinable) 2 Theorem 1715 is just a restatement in terms of linear morphisms of a fact that you might have come across before: every m n matrix can be row-reduced to reduced echelon form using row operations Moreover,

More information

Sparsity of Matrix Canonical Forms. Xingzhi Zhan East China Normal University

Sparsity of Matrix Canonical Forms. Xingzhi Zhan East China Normal University Sparsity of Matrix Canonical Forms Xingzhi Zhan zhan@math.ecnu.edu.cn East China Normal University I. Extremal sparsity of the companion matrix of a polynomial Joint work with Chao Ma The companion matrix

More information

Determinants Chapter 3 of Lay

Determinants Chapter 3 of Lay Determinants Chapter of Lay Dr. Doreen De Leon Math 152, Fall 201 1 Introduction to Determinants Section.1 of Lay Given a square matrix A = [a ij, the determinant of A is denoted by det A or a 11 a 1j

More information

A probabilistic proof of Perron s theorem arxiv: v1 [math.pr] 16 Jan 2018

A probabilistic proof of Perron s theorem arxiv: v1 [math.pr] 16 Jan 2018 A probabilistic proof of Perron s theorem arxiv:80.05252v [math.pr] 6 Jan 208 Raphaël Cerf DMA, École Normale Supérieure January 7, 208 Abstract Joseba Dalmau CMAP, Ecole Polytechnique We present an alternative

More information

RANDOM MATRIX THEORY AND TOEPLITZ DETERMINANTS

RANDOM MATRIX THEORY AND TOEPLITZ DETERMINANTS RANDOM MATRIX THEORY AND TOEPLITZ DETERMINANTS David García-García May 13, 2016 Faculdade de Ciências da Universidade de Lisboa OVERVIEW Random Matrix Theory Introduction Matrix ensembles A sample computation:

More information

A theoretical investigation for low-dimensional molecular-based magnetic materials

A theoretical investigation for low-dimensional molecular-based magnetic materials J. Phys.: Condens. Matter 10 (1998) 3003 3017. Printed in the UK PII: S0953-8984(98)87508-5 A theoretical investigation for low-dimensional molecular-based magnetic materials T Kaneyoshi and Y Nakamura

More information

Bounds for pairs in partitions of graphs

Bounds for pairs in partitions of graphs Bounds for pairs in partitions of graphs Jie Ma Xingxing Yu School of Mathematics Georgia Institute of Technology Atlanta, GA 30332-0160, USA Abstract In this paper we study the following problem of Bollobás

More information

Approximate inference, Sampling & Variational inference Fall Cours 9 November 25

Approximate inference, Sampling & Variational inference Fall Cours 9 November 25 Approimate inference, Sampling & Variational inference Fall 2015 Cours 9 November 25 Enseignant: Guillaume Obozinski Scribe: Basile Clément, Nathan de Lara 9.1 Approimate inference with MCMC 9.1.1 Gibbs

More information

MATH 323 Linear Algebra Lecture 12: Basis of a vector space (continued). Rank and nullity of a matrix.

MATH 323 Linear Algebra Lecture 12: Basis of a vector space (continued). Rank and nullity of a matrix. MATH 323 Linear Algebra Lecture 12: Basis of a vector space (continued). Rank and nullity of a matrix. Basis Definition. Let V be a vector space. A linearly independent spanning set for V is called a basis.

More information

1 Last time: inverses

1 Last time: inverses MATH Linear algebra (Fall 8) Lecture 8 Last time: inverses The following all mean the same thing for a function f : X Y : f is invertible f is one-to-one and onto 3 For each b Y there is exactly one a

More information

CSL361 Problem set 4: Basic linear algebra

CSL361 Problem set 4: Basic linear algebra CSL361 Problem set 4: Basic linear algebra February 21, 2017 [Note:] If the numerical matrix computations turn out to be tedious, you may use the function rref in Matlab. 1 Row-reduced echelon matrices

More information

Angular Momentum set II

Angular Momentum set II Angular Momentum set II PH - QM II Sem, 7-8 Problem : Using the commutation relations for the angular momentum operators, prove the Jacobi identity Problem : [ˆL x, [ˆL y, ˆL z ]] + [ˆL y, [ˆL z, ˆL x

More information

THE UNIVERSITY OF HONG KONG DEPARTMENT OF STATISTICS AND ACTUARIAL SCIENCE

THE UNIVERSITY OF HONG KONG DEPARTMENT OF STATISTICS AND ACTUARIAL SCIENCE THE UNIVERSITY OF HONG KONG DEPARTMENT OF STATISTICS AND ACTUARIAL SCIENCE STAT1301 PROBABILITY AND STATISTICS I EXAMPLE CLASS 2 Review Definition of conditional probability For any two events A and B,

More information

Linear Algebra M1 - FIB. Contents: 5. Matrices, systems of linear equations and determinants 6. Vector space 7. Linear maps 8.

Linear Algebra M1 - FIB. Contents: 5. Matrices, systems of linear equations and determinants 6. Vector space 7. Linear maps 8. Linear Algebra M1 - FIB Contents: 5 Matrices, systems of linear equations and determinants 6 Vector space 7 Linear maps 8 Diagonalization Anna de Mier Montserrat Maureso Dept Matemàtica Aplicada II Translation:

More information

Lecture 2. MATH3220 Operations Research and Logistics Jan. 8, Pan Li The Chinese University of Hong Kong. Integer Programming Formulations

Lecture 2. MATH3220 Operations Research and Logistics Jan. 8, Pan Li The Chinese University of Hong Kong. Integer Programming Formulations Lecture 2 MATH3220 Operations Research and Logistics Jan. 8, 2015 Pan Li The Chinese University of Hong Kong 2.1 Agenda 1 2 3 2.2 : a linear program plus the additional constraints that some or all of

More information

Wavelet Analysis. Willy Hereman. Department of Mathematical and Computer Sciences Colorado School of Mines Golden, CO Sandia Laboratories

Wavelet Analysis. Willy Hereman. Department of Mathematical and Computer Sciences Colorado School of Mines Golden, CO Sandia Laboratories Wavelet Analysis Willy Hereman Department of Mathematical and Computer Sciences Colorado School of Mines Golden, CO 8040-887 Sandia Laboratories December 0, 998 Coordinate-Coordinate Formulations CC and

More information

1. Introduction and notation

1. Introduction and notation Pré-Publicações do Departamento de Matemática Universidade de Coimbra Preprint Number 08 54 DYNAMICS AND INTERPRETATION OF SOME INTEGRABLE SYSTEMS VIA MULTIPLE ORTHOGONAL POLYNOMIALS D BARRIOS ROLANÍA,

More information

Question: Given an n x n matrix A, how do we find its eigenvalues? Idea: Suppose c is an eigenvalue of A, then what is the determinant of A-cI?

Question: Given an n x n matrix A, how do we find its eigenvalues? Idea: Suppose c is an eigenvalue of A, then what is the determinant of A-cI? Section 5. The Characteristic Polynomial Question: Given an n x n matrix A, how do we find its eigenvalues? Idea: Suppose c is an eigenvalue of A, then what is the determinant of A-cI? Property The eigenvalues

More information

arxiv: v1 [math.gt] 15 Mar 2017

arxiv: v1 [math.gt] 15 Mar 2017 ENUMERATION ON RAPH MOSAICS KYUNPYO HON AND SEUNSAN OH arxiv:1703.04868v1 [math.t 15 Mar 2017 Abstract. Since the Jones polynomial was discovered, the connection between knot theory and quantum physics

More information

Generic maximum nullity of a graph

Generic maximum nullity of a graph Generic maximum nullity of a grap Leslie Hogben Bryan Sader Marc 5, 2008 Abstract For a grap G of order n, te maximum nullity of G is defined to be te largest possible nullity over all real symmetric n

More information

x 1 + x 2 2 x 1 x 2 1 x 2 2 min 3x 1 + 2x 2

x 1 + x 2 2 x 1 x 2 1 x 2 2 min 3x 1 + 2x 2 Lecture 1 LPs: Algebraic View 1.1 Introduction to Linear Programming Linear programs began to get a lot of attention in 1940 s, when people were interested in minimizing costs of various systems while

More information

Probabilistic Systems Analysis Spring 2018 Lecture 6. Random Variables: Probability Mass Function and Expectation

Probabilistic Systems Analysis Spring 2018 Lecture 6. Random Variables: Probability Mass Function and Expectation EE 178 Probabilistic Systems Analysis Spring 2018 Lecture 6 Random Variables: Probability Mass Function and Expectation Probability Mass Function When we introduce the basic probability model in Note 1,

More information

Rapid Introduction to Machine Learning/ Deep Learning

Rapid Introduction to Machine Learning/ Deep Learning Rapid Introduction to Machine Learning/ Deep Learning Hyeong In Choi Seoul National University 1/24 Lecture 5b Markov random field (MRF) November 13, 2015 2/24 Table of contents 1 1. Objectives of Lecture

More information

Tutte polynomials of bracelets

Tutte polynomials of bracelets J Algebr Comb (2010) 32: 389 398 DOI 10.1007/s10801-010-0220-1 Tutte polynomials of bracelets Norman Biggs Received: 15 May 2009 / Accepted: 15 February 2010 / Published online: 2 March 2010 Springer Science+Business

More information

CS5314 Randomized Algorithms. Lecture 18: Probabilistic Method (De-randomization, Sample-and-Modify)

CS5314 Randomized Algorithms. Lecture 18: Probabilistic Method (De-randomization, Sample-and-Modify) CS5314 Randomized Algorithms Lecture 18: Probabilistic Method (De-randomization, Sample-and-Modify) 1 Introduce two topics: De-randomize by conditional expectation provides a deterministic way to construct

More information

Quadratic forms. Here. Thus symmetric matrices are diagonalizable, and the diagonalization can be performed by means of an orthogonal matrix.

Quadratic forms. Here. Thus symmetric matrices are diagonalizable, and the diagonalization can be performed by means of an orthogonal matrix. Quadratic forms 1. Symmetric matrices An n n matrix (a ij ) n ij=1 with entries on R is called symmetric if A T, that is, if a ij = a ji for all 1 i, j n. We denote by S n (R) the set of all n n symmetric

More information

Instructions Please answer the five problems on your own paper. These are essay questions: you should write in complete sentences.

Instructions Please answer the five problems on your own paper. These are essay questions: you should write in complete sentences. Instructions Please answer the five problems on your own paper. These are essay questions: you should write in complete sentences.. Recall that P 3 denotes the vector space of polynomials of degree less

More information

A reverse Sidorenko inequality Independent sets, colorings, and graph homomorphisms

A reverse Sidorenko inequality Independent sets, colorings, and graph homomorphisms A reverse Sidorenko inequality Independent sets, colorings, and graph homomorphisms Yufei Zhao (MIT) Joint work with Ashwin Sah (MIT) Mehtaab Sawhney (MIT) David Stoner (Harvard) Question Fix d. Which

More information