Ginzburg-Landau Polynomials and the Asymptotic Behavior of the Magnetization in the Neighborhood of a Tricritical Point

Size: px
Start display at page:

Download "Ginzburg-Landau Polynomials and the Asymptotic Behavior of the Magnetization in the Neighborhood of a Tricritical Point"

Transcription

1 Richard S. Ellis, Asymptotic Behavior of the Magnetization Ginzburg-Landau Polynomials and the Asymptotic Behavior of the Magnetization in the Neighborhood of a Tricritical Point Richard S. Ellis Department of Mathematics and Statistics University of Massachusetts Amherst, MA 0003 Talk for the Workshop Dynamics and Thermodynamics of Systems with Long-Range Interactions: Theory and Experiments Assisi (Perugia), Italy 4-8 July 2007 Collaborators: Jonathan Machta Department of Physics University of Massachusetts Amherst Peter Tak-Hun Otto Department of Mathematics Willamette University Salem, OR rsellis@math.umass.edu

2 Richard S. Ellis, Asymptotic Behavior of the Magnetization 2 Outline of the Talk BEG Spin Model Phase transitions: second order (continuous bifurcation), first order (discontinuous bifurcation), tricritical point Main result Find x and θ such that magnetization m(β n, K n ) x/n θα 0, where (β n, K n ) converges to a second-order point or the tricritical point and α regulates the speed of approach. We work this out for 6 sequences. For each sequence, x is the unique, positive, global minimum point of the Ginzburg-Landau polynomial associated with the sequence. Ginzburg-Landau Polynomials Definition in terms of a limit of scaled free-energy functionals 8 applications. The G-L polynomials rigorize the G-L phenomenology of critical phenomena in the region of phase space through which (β n, K n ) passes. 2. The G-L polynomials reveal features of the phase-transition structure. 3. The asymptotic behavior of m(β n, K n ) 0 is expressed in terms of the unique, positive, global minimum point of G. 4. Scaling limit for total spin with limiting density const exp[ G]. 5. Moderate deviation principle (MDP) for total spin with rate function G const. 6. A consequence of the MDP led us to the asymptotic formula for m(β n, K n ). 7. MDP picks out physically relevant values of α in the asymptotic formula for m(β n, K n ). 8. Through m(β n, K n ) x/n θα, the G-L polynomials make rigorous Riedel s scaling approach to tricritical phase transitions. Contrast behavior near a second-order point versus near the tricritical point. We study 2 sequences (β n, K n ) converging to a second-order point. We study 4 sequences (β n, K n ) converging to the tricritical point. The behavior of the Ginzburg-Landau polynomials reveals the much more complicated phase transition structure near the tricritical point. References Phase Transitions in BEG Model: RSE, P. Otto, H. Touchette, AAP (2005) Limit Theorems in BEG Model: M. Costeniuc, RSE, P. T.-H. Otto, JSP (2007) Asymptotics of the Magnetization in the BEG Model: RSE, J. Machta, P. T.-H. Otto Scaling Approach to Tricritical Phase Transitions: E. K. Riedel, PRL (972)

3 Richard S. Ellis, Asymptotic Behavior of the Magnetization 3 Phase Transitions in the BEG Model K K(β) K(β) K(β ) c.0 K (β) β c β Figure : Sets describing phase-transition structure of the BEG model: the second-order curve {(β, K(β)), 0 < β < β c }, the first-order curve {(β, K (β)), β > β c }, and the tricritical point (β c, K(β c )). The phase-coexistence region consists of all (β, K) above the second-order curve, above the tricritical point, on the first-order curve, and above the first-order curve. The extension of the second-order curve to β > β c is called the spinodal curve. Define M β,k = the set of equilibrium values of the magnetization. Continuous bifurcation (second-order phase transition) for 0 < β < β c = log 4 as K K(β). Discontinuous bifurcation (first-order phase transition) for β > β c as K K (β). Tricritical point at (β c,k(β c )) separates the two phase transitions. Unique phase for (β, K) under the two phase-transition curves: M β,k = {0}. Double phase for (β, K) above the two phase-transition curves: M β,k = {±m(β,k)}, m(β,k) > 0. Triple phase for (β, K) on the first-order curve: M β,k = {±0,m(β,K)}, m(β,k) > 0.

4 Richard S. Ellis, Asymptotic Behavior of the Magnetization 4 Main Result (β n,k n ) (β,k(β)), 0 < β β c, second-order point or the tricritical point. α > 0 regulates speed of approach. Study fine structure of phase transitions. m(β n,k n ) 0. Plausible since M β,k = 0 for (β,k) = (β,k(β)). Find x and θ > 0 such that m(β n,k n ) x/n θα (n θα m(β n,k n ) x). K a 6 4b.09 K(β ) c.08 K(β) 4d.07.3 β c Figure 2: Sequences and 2 converge to a second-order point, sequences 3 6 converge to the tricritical point..4.5 K (β) 4c β 6 sequences (β n,k n ) with different asymptotic behavior of m(β n,k n ). Sequence : along ray above tangent line (K (β)b k < 0) β n = β + b/n α, K n = K(β) + k/n α 2. Sequence 2: along curve = second-order curve to order p 3. Sequence 3: along ray above tangent line (seq. for β = β c ) 4. Sequences 4a 4d: along curve tangent to spinodal curve (b > 0, l l c ) β n = β + b/n α, K n = K(β c ) + K (β c )b/n α + lb 2 /2n 2α 5. Sequence 5: along curve tangent to second-order curve (seq. 4 with b < 0, l > K (β c )) 6. Sequence 6: along curve = second-order curve to order p

5 Richard S. Ellis, Asymptotic Behavior of the Magnetization 5. Sequence : along ray above tangent line (K (β)b k < 0) β n = β + b/n α, K n = K(β) + k/n α 2. Sequence 2: along curve = second-order curve to order p (b 0, lb p > K (p) (β)b p ) [higher order correction to sequence ] β n = β + b/n α, K n = K(β) + p j= K(j) (β)b j /j!n jα + lb p /p!n pα 3. Sequence 3: along ray above tangent line (seq. for β = β c ) 4. Sequences 4a 4d: along curve tangent to spinodal curve (b > 0, l l c ) [higher order correction to sequence 3 for b > 0] β n = β + b/n α, K n = K(β c ) + K (β c )b/n α + lb 2 /2n 2α 5. Sequence 5: along curve tangent to second-order curve (seq. 4 with b < 0, l > K (β c )) 6. Sequence 6: along curve = second-order curve to order p (seq. 2 for β = β c with b < 0, lb p > K (p) (β)b p ) [higher order correction to sequence 3 for b < 0] Theorem. Asymptotic behavior of m(β n, K n ) for the 6 sequences. (a) Similar sequences give rise to different asymptotic behavior of m(β n,k n ) 0: compare and 3, 2 and 6 (influence of first-order curve). sequence converges to theorem asymptotic behavior in paper of m(β n, K n ) second-order point Thm. 3. m(β n,k n ) x/n α/2 2 second-order point Thm. 3.2 m(β n,k n ) x/n pα/2 3 tricritical point Thm. 6. m(β n,k n ) x/n α/4 4a 4d tricritical point Thm. 6.2 m(β n,k n ) x/n α/2 5 tricritical point Thm. 6.3 m(β n,k n ) x/n α/2 6 tricritical point Thm. 6.4 m(β n,k n ) x/n (p )α/2 (b) For each sequence (β n,k n ), x equals the unique, positive, global minimum point of the Ginzburg-Landau polynomial associated with (β n,k n ).

6 Richard S. Ellis, Asymptotic Behavior of the Magnetization 6 Ginzburg-Landau Polynomials Theorem. Free-energy functional and equilibrium values of the magnetization. Let G β,k denote the free-energy functional (see p. 8). Then M β,k = the set of global minimum points of G β,k. Definition of Ginzburg-Landau Polynomial G. Let α > 0 regulate the speed (β n,k n ) (β,k(β)) for 0 < β β c. Seek u R, γ > 0, and a polynomial G of degree 4 or 6 such that for all x R lim n u G βn,k n (x/n γ ) = G(x). n 8 Applications of the Ginzburg-Landau Polynomials. The Ginzburg-Landau polynomials G rigorize the Ginzburg-Landau phenomenology of critical phenomena in the region of phase space through which (β n,k n ) passes. 2. The Ginzburg-Landau polynomials G enable us to discover features of the phase-transition structure. 3. Asymptotic behavior of m(β n,k n ) is expressed in terms of the unique, positive, global minimum point x of G: m(β n,k n ) x/n θα. 4. Scaling limit with limiting density const exp[ G]: {S n /n γ dx} = exp[ G(x)]dx/ R exp[ G(y)]dy. 5. Moderate deviation principle (MDP) with rate function G const: {S n /n γ dx} exp[ n u (G(x) Ḡ)]dx, Ḡ = inf y R G(y). 6. We were led to the asymptotic formula m(β n,k n ) x/n θα by a consequence of the MDP: {S n /n γ dx} = ( 2 δ x + 2 δ x) (dx). 7. MDP picks out physically relevant values of α in the asymptotic formula m(β n,k n ) x/n θα. 8. Through m(β n,k n ) x/n θα, the Ginzburg-Landau polynomials make rigorous Riedel s scaling approach to tricritical phase transitions.

7 Richard S. Ellis, Asymptotic Behavior of the Magnetization 7 BEG Spin Model Mean-field approximation to Blume Capel model (966) and Blume Emery Griffiths model (97): study He 3 -He 4 mixtures and other physical systems. n spins ω i Λ = {,0,} 2. Microstates: ω = (ω,ω 2,...,ω n ) Λ n 3. Total spin: S n (ω) = n j= ω j 4. Hamiltonian or energy function (K > 0): H n,k (ω) = n j= ω2 j K n 5. Prior measure: n j,k= ω jω k = n j= ω2 j nk(s n(ω)/n) 2 P n (ω) = n j= ρ(ω j) = /3 n, ρ = 3 (δ + δ 0 + δ ) 6. Canonical ensemble (β > 0, K > 0): 7. Alignment effect: P n,β,k (ω) = Z n (β,k) exp[ βh n,k(ω)] /3 n, Z n (β,k) = ω Λ n exp[ βh n,k (ω)] /3 n where ρ β (ω j ) = e βω2 j /( + 2e β ), lim P n,β,k (ω) = P n (ω) = /3 n, β 0 + n lim P n,β,k (ω) = K 0 + j= ρ β(ω j ), lim P n,β,k (ω) = 2 (δ ω (ω) + δ ω +(ω)), K where for all j, ω j = and ω+ j = 8. Phase transition: persistence of alignment effect as n

8 Richard S. Ellis, Asymptotic Behavior of the Magnetization 8 Second-Order Phase Transition for β < β c Theorem. Large deviation principle and M β,k. (a) For β > 0, K > 0, we have the LDP with rate function I(β,K): P n,β,k {S n /n dx} exp[ ni β,k (x)]dx. (b) If I β,k (x) > 0, then P n,β,k {S n /n x} 0. (c) Define M β,k = {x R : I β,k (x) = 0} and the free-energy functional ( ) + e G β,k (x) = βkx 2 β (e 2βKx + e 2βKx ) log + 2e β = βkx 2 log Λ e2βkω ρ β (dω ) Then M β,k = {x R : x a global minimum point of G β,k } = M Gβ,K. (a) (b) Figure 3: Graphs of G β,k for 0 < β β c = log 4. (a) 0 < K K(β), (b) K > K(β) Theorem. Second-order phase transition for 0 < β β c = log 4. (a) For 0 < K K(β) = (e β + 2)/4β, M β,k = {0} (b) For K > K(β), m(β,k) > 0 such that M β,k = {±m(β,k)}. (c) m(β,k) 0 as K (K(β)) +. (d) M β,k has a continuous bifurcation at K = K(β) [second-order phase transition].

9 Richard S. Ellis, Asymptotic Behavior of the Magnetization 9 First-Order Phase Transition for β > β c (a) (b) Figure 4: Graphs of G β,k for β > β c. (a) 0 < K < κ(β), (b) κ(β) < K < K (β) Figure 5: Graph of G β,k for β > β c and K = K (β) (a) (b) Figure 6: Graphs of G β,k for β > β c. (a) K (β) < K < K(β), (b) K K(β) Theorem. First-order phase transition for β > β c = log 4. (a) For 0 < K < K (β), M β,k = {0}. (b) For K = K (β), m(β,k (β)) > 0 such that M β,k (β) = {0, ±m(β,k (β))}. (c) For K > K (β), m(β,k) > 0 such that M β,k = {±m(β,k)}. (d) m(β,k) m(β,k (β)) > 0 as K (K (β)) +. (e) M β,k has a discontinuous bifurcation at K = K (β) [first-order phase transition].

10 Richard S. Ellis, Asymptotic Behavior of the Magnetization 0 Ginzburg-Landau Phenomenology and Polynomials, I Phenomenology. For 0 < β < β c, K near K(β), and x near 0, c 4 (β) > 0 G β,k (x) G (2) β,k (0)x2 /2 + G (4) β,k (0)x4 /4! G β,k (x) = β(k(β) K)x 2 + c 4 (β)x 4. M Gβ,K has same form as M Gβ,K : point for 0 < K K(β), 2 points for K > K(β). Theorem. Ginzburg-Landau polynomial for sequence. For 0 < β < β c, α > 0, and K (β)b k < 0, define β n = β + b/n α, K n = K(β) + k/n α. (a) For u = 2α and γ = α/2 n u G βn,k n (x/n γ ) G(x) = β(k (β)b k)x 2 + c 4 (β)x 4, c 4 (β) > 0. (b) M G = {± x}, where x = β(k K (β)b)/2c 4 (β)). Proof. (a) Two-term Taylor expansion and asymptotics give Since n u G βn,k n (x/n γ ) n 2γ +uβ(k(β n) K n )x 2 + n 4γ +uc 4(β)x 4. K(β n ) K n = K(β + b/n α ) (K(β) + k/n α ) (K (β)b k)/n α, we have n u G βn,k n (x/n γ ) n 2γ+α +uβ(k (β)b k)x 2 + n 4γ +uc 4(β)x 4. Choose 2γ + α + u = 0 = 4γ + u or γ = α/2, u = 4γ = 2α, then n. (b) Solve G ( x) = 0.

11 Richard S. Ellis, Asymptotic Behavior of the Magnetization Ginzburg-Landau Polynomial Mirrors Phase-Transition Structure, I Theorem. Ginzburg-Landau polynomial for sequence. For 0 < β < β c, α > 0, and K (β)b k < 0, define β n = β + b/n α, K n = K(β) + k/n α. (a) For u = 2α and γ = α/2 n u G βn,k n (x/n γ ) G(x) = β(k (β)b k)x 2 + c 4 (β)x 4, c 4 (β) > 0. (b) M G = {± x}, where x = β(k K (β)b)/2c 4 (β)). (a) (b) Figure 7: Graphs of G β,k for 0 < β β c = log 4 and graphs of G. (a) Graph (a) shows the graph of G β,k for 0 < K K(β) and the graph of G for K (β)b k 0. (β n,k n ) (β,k(β)) from the single-phase region below the second-order curve. We have M Gβn,Kn = {0} = M G. (b) Graph (b) shows the graph of G β,k for K > K(β) and the graph of G for K (β)b k < 0. (β n,k n ) (β,k(β)) from the phase-coexistence region above the second-order curve. We have M Gβn,Kn = {±m(β n,k n )}, M G = {± x}, M Gβn,Kn = 2 = M G.

12 Richard S. Ellis, Asymptotic Behavior of the Magnetization 2 Asymptotics of m(β n, K n ) for Sequence Theorem. Ginzburg-Landau polynomial for sequence. For 0 < β < β c, α > 0, and K (β)b k < 0, define β n = β + b/n α, K n = K(β) + k/n α. (a) For u = 2α and γ = α/2 n u G βn,k n (x/n γ ) G(x) = β(k (β)b k)x 2 + c 4 (β)x 4, c 4 (β) > 0. (b) M G = {± x}, where x = β(k K (β)b)/2c 4 (β)). In the next theorem we give the asymptotics of m(β n,k n ) 0 for sequence. Theorem. Asymptotics of m(β n, K n ) for sequence. For 0 < β < β c, α > 0, and K (β)b k < 0, define Then β n = β + b/n α, K n = K(β) + k/n α. m(β n,k n ) x/n α/2. Proof. We write m n = n γ m(β n,k n ) and G n (x) = n u G βn,k n (x/n γ ). Since m(β n,k n ) is a global minimum point of G βn,k n, m n is a global minimum point of G n. We assert without proof that there exists x 0 such that m n x. Since G n ( m n ) G n (y) for all y and G n (x) G(x) for all x, G( x) = lim G n ( m n ) lim G n (y) = G(y). n n Thus x is a nonnegative global minimum point of G. It follows that x = x, the unique, nonnegative, global minimum point of G. Since γ = α/2, we obtain the desired conclusion: m n = n γ m(β n,k n ) = n α/2 m(β n,k n ) x; i.e.,m(β n,k n ) x/n α/2.

13 Richard S. Ellis, Asymptotic Behavior of the Magnetization 3 Ginzburg-Landau Phenomenology and Polynomials, II Phenomenology. For β > β c, β near β c, K near K(β c ), and x near 0, there exists c 4 > 0 and c 6 > 0 such that G β,k (x) G (2) β,k (0)x2 /2 + G (4) β,k (0)x4 /4! + G (6) β,k (0)x6 /6! G β,k (x) = β c (K(β) K)x 2 + c 4 (4 e β )x 4 + c 6 x 6. M G β,k has same form as M Gβ,K : point for 0 < K < K(β), 3 points for K = K (β), 2 points for K > K(β). K a 6 4b.09 K(β c).08 K(β) 4d.07 K (β) 4c.3 β c.4.5 β Theorem. Ginzburg-Landau polynomial for sequences 4a 4d. For β > β c, b > 0, α > 0, and l l c, define β n = β + b/n α, K n = K(β c ) + K (β)b/n α + lb 2 /2n 2 α. (a) For u = 3α and γ = α/2 n u G βn,k n (x/n γ ) G l (x) = 2 β c(k (β c ) l)b 2 x 2 4c 4 bx 4 + c 6 x 6. (b) For l = l c, M Glc = {0, ± x(l c )}. (c) For l > l c, M Gl = {± x(l)}. Proof. (a) Three-term Taylor expansion and asymptotics give n u G βn,k n (x/n γ ) n 2γ +uβ c(k(β n ) K n )x 2 + n 4γ +uc 4(4 e β n )x 4 + n 6γ c 6x 6. Since K(β n ) K n = K(β c + b/n α ) K n (K (β c ) l)b 2 /2n 2α, 4 e β n 4b/n α, n u G βn,k n (x/n γ ) n β 2γ+2α +u 2 c(k (β c ) l)x 2 n 4γ+α +u4c 4bx 4 + n 6γ +uc 6x 6. Choose 2γ + 2α + u = 0 = 4γ + α + u = 6γ + u or γ = α/2, u = 6γ = 3α, then n. (b) (c) Solve G l ( x) = 0.

14 Richard S. Ellis, Asymptotic Behavior of the Magnetization 4 Ginzburg-Landau Polynomial Mirrors Phase-Transition Structure, II Theorem. Ginzburg-Landau polynomial G l for sequences 4a 4d. For β > β c, b > 0, α > 0, and l l c = K (β c ) 5/4β c, define β n = β + b/n α, K n = K(β c ) + K (β)/n α + lb 2 /2n 2 α. (a) For u = 3α and γ = α/2 n u G βn,k n (x/n γ ) G l (x) = 2 β c(k (β c ) l)b 2 x 2 4c 4 bx 4 + c 6 x 6. (b) For l = l c, M Glc = {0, ± x(l c )}. (c) For l > l c, M Gl = {± x(l)}. Since β n β c = b/n α, (β n,k n ) (β c,k(β c )) along the curve (β, K(β)), where for β > β c K(β) = K(β c ) + K (β c )(β β c ) + l(β β c ) 2 /2. K a 6 4b.09 K(β c).08 K(β) 4d.07 K (β) 4c.3 β c.4.5 β (a) Sequence 4a corresponds to l > K (β c ). (b) Sequence 4b corresponds to l = K (β c ). (c) Sequence 4c corresponds to l = l c = K (β c ) 5/4β c ; G lc has 3 global minimum points. (d) Sequence 4d corresponds to l c < l < K (β c ).

15 Richard S. Ellis, Asymptotic Behavior of the Magnetization 5 (a) (b) Figure 8: Graphs of G β,k for β > β c, 0 < K < K (β) and graphs of G l for l < l c. Figure 9: Graph of G β,k for β > β c, K = K (β) and graph of G l for l = l c. (a) (b) Figure 0: Graphs of G β,k for β > β c, K > K (β) and graphs of G l for l > l c.. Figure 9. For l = l c, M Glc = {0, ± x(l c )}. We argue that (β n,k n ) (β c,k(β c )) along (β, K(β)), which coincides with first-order curve to order 2. This implies K (β c) = l c. M Glc = {0, ± x(l c )} mirrors M = {0, ±m(β Gβn,Kn n,k n )} for (β n,k n ) on first-order curve. 2. Figure 8. For l < l c, (β n,k n ) (β c,k(β c )) from the single-phase region under first-order curve. We have M Gβn,Kn = {0} = M G l. 3. Figure 0. For l > l c, (β n,k n ) (β c,k(β c )) from the phase-coexistence region above first-order curve. We have M Gβn,Kn = {±m(β n,k n )}, M Gl = {± x(l)}, M Gβn,Kn = 2 = M G l. 4. M Gl undergoes a first-order phase transition at l = l c, mirroring the first-order phase transition in M Gβn,Kn at K = K (β n ) for β n > β c.

16 Richard S. Ellis, Asymptotic Behavior of the Magnetization 6 Asymptotics of m(β n, K n ) for Sequences 4a 4d Theorem. Ginzburg-Landau polynomial for sequences 4a 4d. For β > β c, b > 0, α > 0, and l l c, define β n = β + b/n α, K n = K(β c ) + K (β)b/n α + lb 2 /2n 2 α. (a) For u = 3α and γ = α/2 n u G βn,k n (x/n γ ) G l (x) = 2 β c(k (β c ) l)b 2 x 2 4c 4 bx 4 + c 6 x 6. (b) For l = l c, M Glc = {0, ± x(l c )}. (c) For l > l c, M Gl = {± x(l)}. In the next theorem we give the asymptotics of m(β n,k n ) 0 for sequences 4a 4d. These are the only sequences for which higher order corrections do not affect the asymptotic behavior of m(β n,k n ). Theorem. Asymptotics of m(β n, K n ) for sequence 4a 4d. For β > β c, b > 0, α > 0, and l l c, let (β n,k n ) be the sequence in the preceding theorem. Then m(β n,k n ) x(l)/n α/2. Proof. For l > l c, the proof is the same as the proof for sequence (see p. 0). The proof for l = l c is much harder. We write m n = n γ m(β n, K n ) and G n (x) = n u G βn,k n (x/n γ ). Since m(β n, K n ) is a global minimum point of G βn,k n, m n is a global minimum point of G n. We assert without proof that there exists x 0 such that m n x. Since G n ( m n ) G n (y) for all y and G n (x) G lc (x) for all x, G lc ( x) = lim G n ( m n ) lim G n (y) = G lc (y). n n Thus x is a nonnegative global minimum point of G lc. In the proof for sequence, we were able to assert that x = x, the unique, nonnegative, global minimum point of G. However, G lc has two nonnegative global minimum points, 0 and x(l c ). If x = 0, then we do not obtain the desired asymptotic formula. In order to prove that x = x(l c ), one must use results from complex analysis, including Rouché s theorem. After one shows that x = x(l c ), the proof is completed as before. Since γ = α/2, the desired conclusion follows: m n = n γ m(β n, K n ) = n α/2 m(β n, K n ) x(l c ); i.e., m(β n, K n ) x(l c )/n α/2.

17 Richard S. Ellis, Asymptotic Behavior of the Magnetization 7 Ginzburg-Landau Polynomial and Scaling Limit Theorem. Ginzburg-Landau polynomial for sequence. For 0 < β < β c, α > 0, and K (β)b k < 0, define β n = β + b/n α, K n = K(β) + k/n α. (a) For u = 2α and γ = α/2 n u G βn,k n (x/n γ ) G(x) = β(k (β)b k)x 2 + c 4 (β)x 4, c 4 (β) > 0. (b) M G = {± x}, where x = β(k K (β)b)/2c 4 (β)). Theorem. Scaling limit for S n /n 3/4 for sequence. For 0 < β < β c, α > 0, and K (β)b k < 0, let (β n,k n ) be the sequence in the preceding theorem. Then {S n /n 3/4 dx} = exp[ G(x)]dx/ R exp[ G(y)]dy; i.e., for any bounded, continuous function f lim n Λ n f(s n/n 3/4 ) d = R f(x) exp[ G(x)]dx/ R exp[ G(y)]dy. Proof. In part (a) of the first theorem on this page, set u = 0. Then for γ = /4 ng βn,k n (x/n γ ) = ng βn,k n (x/n /4 ) G(x) = β(k (β)b k)x 2 + c 4 (β)x 4, c 4 (β). Omitting several technical steps, we have for any bounded, continuous function f lim n Λ n f(s n/n 3/4 )d = lim f(x) exp[ ng β n R n,k n (x/n /4 )] dx/ exp[ ng β R n,k n (y/n /4 )] dy = R f(x)exp[ G(x)] dx/ exp[ G(y)] dy. R This completes the proof. The first equality is a consequence of the Hubbard-Stratonovitch transformation.

18 Richard S. Ellis, Asymptotic Behavior of the Magnetization 8 Ginzburg-Landau Polynomial and Moderate Deviation Principle Theorem. Ginzburg-Landau polynomial for sequence. For 0 < β < β c, α > 0, and K (β)b k < 0, define β n = β + b/n α, K n = K(β) + k/n α. (a) For u = 2α = 4γ and γ = α/2 n u G βn,k n (x/n γ ) G(x) = β(k (β)b k)x 2 + c 4 (β)x 4, c 4 (β) > 0. (b) M G = {± x}, where x = β(k K (β)b)/2c 4 (β)). The weak limit in part (b) of the next theorem led us to the asymptotic formula for m(β n,k n ) in the theorem on page 2. Theorem. Moderate deviation principle for S n /n γ for sequence. For 0 < β < β c, γ (0,/4), α = 2γ (0,/2), and K (β)b k < 0, let (β n,k n ) be the sequence in the preceding theorem. (a) Define Γ(x) = G(x) min y R G(y). We have the MDP {S n /n γ dx} exp[ n 4γ Γ(x)]dx, (b) We have the weak limit {S n /n γ dx} = ( 2 δ x + 2 δ x) (dx). Definition of MDP and Laplace Principle. If A is any set in R, then Γ(A) denotes the infimum of Γ over A. The MDP states for every closed set F in R lim sup n n log P 4γ n,β n,k n {S n /n γ F } Γ(F) and that for every open set U in R lim inf n n log P 4γ n,β n,k n {S n /n γ U} Γ(U). The MDP is equivalent to the Laplace principle, which states that for any bounded, continuous function ψ lim n n 4γ Λ n exp[ n 4γ ψ(s n /n γ )]d = sup{ψ(x) Γ(x)}. x R

19 Richard S. Ellis, Asymptotic Behavior of the Magnetization 9 Proof of the Laplace principle. We prove that for any bounded, continuous function ψ {ψ(x) Γ(x)}. n 4γ lim n Λ n exp[ n 4γ ψ(s n /n γ )] d = sup x R According to part (a) of the theorem on page 8, for u = 2α = 4γ and γ = α/2 n u G βn,k n (x/n γ ) = n 4γ G βn,k n (x/n γ ) G(x) = β(k (β)b k)x 2 + c 4 (β)x 4, c 4 (β) > 0. For the Laplace principle and thus the MDP, we require u = 4γ > 0, which restricts γ and α = 2γ to the respective open intervals (0, /4) and (0, /2). Omitting several technical steps, we have for any bounded, continuous function ψ Λ n exp[ n 4γ ψ(s n /n γ )] d R exp[n 4γ {ψ(x) n 4γ G βn,k n (x/n γ )}] dx/ R exp[n 4γ {ψ(x) G(x)}] dx/ R exp[n 4γ { G(x)}] dx exp[n 4γ sup x R {ψ(x) G(x)}] / exp[ n 4γ inf x R {G(x)}] = exp[n 4γ sup x R {ψ(x) Γ(x)}]. R exp[n 4γ { n 4γ G βn,k n (x/n γ )}] dx The next to last line is a consequence of Laplace s method for the asymptotics of integrals over R. The limit giving the Laplace principle follows from this display. Proof of weak limit in part (b). We have G(± x) = min y R G(y) and Γ(± x) = 0 < Γ(y) for all y ± x. Let F be any closed set in R not containing ± x. Then Γ(F) > 0, and the moderate deviation upper bound implies that for all sufficiently large n {S n /n γ F } exp[ n 4γ Γ(F)] 0. A separate argument shows that the weak limit of the sequence {S n /n γ dx} exists. The last display implies that this weak limit equals ( 2 δ x + 2 δ x)(dx). Weak limit in part (b) led us to asymptotic formula for m(β n, K n ). According to the theorem on page 2, m(β n, K n ) x/n α/2 for any α > 0. To derive this asymptotic formula, we start with the heuristic formula {S n /n dx} ( 2 δ m(β n,k n) + 2 δ m(β n,k n))(dx), which expresses the fact that ±m(β n, K n ) are the equilibrium values of S n /n. Replace dx with dx/n γ and obtain {S n /n γ dx} ( 2 δ n γ m(β n,k n) + 2 δ n γ m(β n,k n))(dx). Compare this with the weak limit in part (b), which holds for α (0, /2): ( {S n /n γ dx} = ) 2 δ x + δ 2 x (dx). This leads to the desired asymptotic formula for α (0, /2): n γ m(β n, K n ) x(l) or m(β n, K n ) x/n α/2.

20 Richard S. Ellis, Asymptotic Behavior of the Magnetization 20 Ginzburg-Landau Polynomial and Moderate Deviation Principle Theorem. Ginzburg-Landau polynomial for sequence. For 0 < β < β c, α > 0, and K (β)b k < 0, define β n = β + b/n α, K n = K(β) + k/n α. (a) For u = 2α = 4γ and γ = α/2 n u G βn,k n (x/n γ ) G(x) = β(k (β)b k)x 2 + c 4 (β)x 4, c 4 (β) > 0. (b) M G = {± x}, where x = β(k K (β)b)/2c 4 (β)). The weak limit in part (b) of the next theorem led us to the asymptotic formula for m(β n,k n ) in the theorem on page 2. Theorem. Moderate deviation principle for S n /n γ for sequence. For 0 < β < β c, γ (0,/4), α = 2γ (0,/2), and K (β)b k < 0, let (β n,k n ) be the sequence in the preceding theorem: β n = β + b/n α, K n = K(β) + k/n α. (a) Define Γ(x) = G(x) inf y R G(y). We have the MDP {S n /n γ dx} exp[ n 4γ Γ(x)]dx, (b) We have the weak limit {S n /n γ dx} = ( 2 δ x + 2 δ x) (dx). Definition of MDP. If A is any set in R, then Γ(A) denotes the infimum of Γ over A. The MDP states for every closed set F in R lim sup n and that for every open set U in R lim inf n n 4γ log P n,β n,k n {S n /n γ F } Γ(F) n 4γ log P n,β n,k n {S n /n γ U} Γ(U).

21 Richard S. Ellis, Asymptotic Behavior of the Magnetization 2 Theorem. Moderate deviation principle for S n /n γ for sequence. For 0 < β < β c, γ (0,/4), α = 2γ (0,/2), and K (β)b k < 0, define β n = β + b/n α, K n = K(β) + k/n α and the Ginzburg-Landau polynomial G(x) = β(k (β)b k)x 2 + c 4 (β)x 4. (a) Define Γ(x) = G(x) inf y R G(y). We have the MDP {S n /n γ dx} exp[ n 4γ Γ(x)]dx, (b) We have the weak limit {S n /n γ dx} = ( 2 δ x + 2 δ x) (dx). Proof of weak limit in part (b). We have G(± x) = min y R G(y) and Γ(± x) = 0 < Γ(y) for all y ± x. Let F be any closed set in R not containing ± x. Then Γ(F) > 0, and for all sufficiently large n {S n /n γ F } exp[ n 4γ Γ(F)] 0. A separate argument shows that the weak limit of the sequence {S n /n γ dx} exists. The last display implies that this weak limit equals ( 2 δ x + 2 δ x)(dx). Weak limit in part (b) led us to asymptotic formula for m(β n, K n ). Theorem on page 2: for any α > 0,m(β n,k n ) x/n α/2. Since ±m(β n,k n ) = equilibrium values of S n /n, we have heuristic formula {S n /n dx} ( 2 δ m(β n,k n ) + 2 δ m(β n,k n ))(dx), Replace dx with dx/n γ and obtain {S n /n γ dx} ( 2 δ n γ m(β n,k n ) + 2 δ n γ m(β n,k n ))(dx). Compare with weak limit in part (b), which holds for α (0, /2): {S n /n γ dx} = ( 2 δ x + 2 δ x) (dx). This leads to n γ m(β n,k n ) x or m(β n,k n ) x/n α/2 for α (0,/2). Note different ranges of α. Here α (0,/2) versus any α > 0 in the theorem on page 2. The MDP is telling us something physically important.

22 Richard S. Ellis, Asymptotic Behavior of the Magnetization 22 MDP Picks Out Physically Relevant Values of α For 0 < β < β c, γ (0,/4), α = 2γ (0,/2), K (β)b k < 0, define β n = β + b/n α and K n = K(β) + k/n α. Case. α > /2. See figure (a). (β n,k n ) (β,k(β)) quickly, and we are in the critical regime, where finite-size scaling matters. m(β n,k n ), the positive global minimum point of G βn,k n, is much smaller than the fluctuations of S n /n: m(β n,k n ) x/n α/2 E n,βn,k n S n /n const/n /4. m(β n,k n ) has no physical meaning. S n /n has large fluctuations around 0. (a) (b) Gβ n, K n P n, β n, K( S dx) n n/n - m( β n, K n ) m( β n, K n ) - m( β n, K n ) m( β n, K n ) Figure : G βn,k n and {S n /n dx}. (a) α /2, (b) 0 < α < /2 Case 2. 0 < α < /2 (values of α for MDP). See figure (b). (β n,k n ) (β,k(β)) slowly, and we are in the two-phase regime, where the system is effectively infinite. G βn,k n has two, deep, global minimum points at ±m(β n,k n ), and the probability distribution of S n /n is sharply peaked at ±m(β n,k n ): m n x n α/2 E n,βn,k n { Sn /n m n Sn /n > δm n } const, δ (0,). n /2 α/2 E n,βn,k n S n /n m(β n,k n ) x/n α/2. m(β n,k n ) has physical meaning because it is the expected value of a measurement of S n /n.

23 Richard S. Ellis, Asymptotic Behavior of the Magnetization 23 Scaling Approach to Tricritical Phase Transitions K a 6 4b.09 K(β c).08 K(β) 4d.07 K (β) 4c.3 β c.4.5 β sequence converges to theorem asymptotic behavior in paper of m(β n, K n ) second-order point Thm. 3. m(β n, K n ) x/n α/2 2 second-order point Thm. 3.2 m(β n, K n ) x/n pα/2 3 tricritical point Thm. 6. m(β n, K n ) x/n α/4 4a 4d tricritical point Thm. 6.2 m(β n, K n ) x/n α/2 5 tricritical point Thm. 6.3 m(β n, K n ) x/n α/2 6 tricritical point Thm. 6.4 m(β n, K n ) x/n (p )α/2 Connection with scaling theory Sequence. /n α = distance from (β n,k n ) to second-order curve, and critical exponent for m(β n,k n ) is /2. Thus m(β n,k n ) 0 like /n α/2. This is the case p = 2 of sequence 2. Sequence 2. /n pα = distance from (β n,k n ) to second-order curve, and critical exponent for m(β n,k n ) is /2. Thus m(β n,k n ) 0 like /n pα/2. Sequence 3. /n α = distance from (β n,k n ) to first-order curve, and critical exponent for m(β n,k n ) is /4. Thus m(β n,k n ) 0 like /n α/4. Sequences 4a 4d. This can be seen via Riedel s scaling approach to tricritical phase transitions for region I in his paper. Sequence 5. This is the case p = 2 of sequence 6. Sequence 6. This can be seen via Riedel s scaling approach to tricritical phase transitions for region II in his paper.

24 Richard S. Ellis, Asymptotic Behavior of the Magnetization 24 Summary of Talk Theorem. Free-energy functional and equilibrium values of the magnetization. Let G β,k denote the free-energy functional (see p. 8). Then M β,k = the set of global minimum points of G β,k. Definition of Ginzburg-Landau Polynomial G. Let α > 0 regulate the speed (β n,k n ) (β,k(β)) for 0 < β β c. Seek u R, γ > 0, and a polynomial G of degree 4 or 6 such that for all x R lim n u G βn,k n (x/n γ ) = G(x). n 8 Applications of the Ginzburg-Landau Polynomials. The Ginzburg-Landau polynomials G rigorize the Ginzburg-Landau phenomenology of critical phenomena in the region of phase space through which (β n,k n ) passes. 2. The Ginzburg-Landau polynomials G enable us to discover features of the phase-transition structure. 3. Asymptotic behavior of m(β n,k n ) is expressed in terms of the unique, positive, global minimum point x of G: m(β n,k n ) x/n θα. 4. Scaling limit with limiting density const exp[ G]: {S n /n γ dx} = exp[ G(x)]dx/ R exp[ G(y)]dy. 5. Moderate deviation principle (MDP) with rate function G const: {S n /n γ dx} exp[ n u (G(x) Ḡ)]dx, Ḡ = inf y R G(y). 6. We were led to the asymptotic formula m(β n,k n ) x/n θα by a consequence of the MDP: {S n /n γ dx} = ( 2 δ x + 2 δ x) (dx). 7. MDP picks out physically relevant values of α in the asymptotic formula m(β n,k n ) x/n θα. 8. Through m(β n,k n ) x/n θα, the Ginzburg-Landau polynomials make rigorous Riedel s scaling approach to tricritical phase transitions.

Asymptotic Behavior of the Magnetization Near Critical and Tricritical Points via Ginzburg-Landau Polynomials

Asymptotic Behavior of the Magnetization Near Critical and Tricritical Points via Ginzburg-Landau Polynomials Asymptotic Behavior of the Magnetization Near Critical and Tricritical Points via Ginzburg-Landau Polynomials Richard S. Ellis rsellis@math.umass.edu Jonathan Machta 2 machta@physics.umass.edu Peter Tak-Hun

More information

Ginzburg-Landau Polynomials and the Asymptotic Behavior of the Magnetization Near Critical and Tricritical Points

Ginzburg-Landau Polynomials and the Asymptotic Behavior of the Magnetization Near Critical and Tricritical Points Ginzburg-Landau Polynomials and the Asymptotic Behavior of the Magnetization Near Critical and Tricritical Points Richard S. Ellis rsellis@math.umass.edu Jonathan Machta 2 machta@physics.umass.edu Peter

More information

ASYMPTOTIC BEHAVIOR OF THE FINITE-SIZE MAGNETIZATION AS A FUNCTION OF THE SPEED OF APPROACH TO CRITICALITY

ASYMPTOTIC BEHAVIOR OF THE FINITE-SIZE MAGNETIZATION AS A FUNCTION OF THE SPEED OF APPROACH TO CRITICALITY The Annals of Applied Probability 200, Vol. 20, No. 6, 28 26 DOI: 0.24/0-AAP679 Institute of Mathematical Statistics, 200 ASYMPTOTIC BEHAVIO OF THE FINITE-SIZE MAGNETIZATION AS A FUNCTION OF THE SPEED

More information

Path Coupling and Aggregate Path Coupling

Path Coupling and Aggregate Path Coupling Path Coupling and Aggregate Path Coupling 0 Path Coupling and Aggregate Path Coupling Yevgeniy Kovchegov Oregon State University (joint work with Peter T. Otto from Willamette University) Path Coupling

More information

ANALYSIS OF PHASE TRANSITIONS IN THE MEAN-FIELD BLUME EMERY GRIFFITHS MODEL

ANALYSIS OF PHASE TRANSITIONS IN THE MEAN-FIELD BLUME EMERY GRIFFITHS MODEL The Annals of Applied Probability 2005, Vol. 15, No. 3, 2203 2254 DOI 10.1214/105051605000000421 Institute of Mathematical Statistics, 2005 ANALYSIS OF PHASE TRANSITIONS IN THE MEAN-FIELD BLUME EMERY GRIFFITHS

More information

ENSEMBLE EQUIVALENCE AND PHASE TRANSITIONS FOR GENERAL MODELS IN STATISTICAL MECHANICS AND FOR THE CURIE-WEISS-POTTS MODEL

ENSEMBLE EQUIVALENCE AND PHASE TRANSITIONS FOR GENERAL MODELS IN STATISTICAL MECHANICS AND FOR THE CURIE-WEISS-POTTS MODEL ENSEMBLE EQUIVALENCE AND PHASE TRANSITIONS FOR GENERAL MODELS IN STATISTICAL MECHANICS AND FOR THE CURIE-WEISS-POTTS MODEL A Dissertation Presented by MARIUS F. COSTENIUC Submitted to the Graduate School

More information

The Generalized Canonical Ensemble and Its Universal Equivalence with the Microcanonical Ensemble

The Generalized Canonical Ensemble and Its Universal Equivalence with the Microcanonical Ensemble Journal of Statistical Physics, Vol. 119, Nos. 5/6, June 2005 ( 2005) DOI: 10.1007/s10955-005-4407-0 The Generalized Canonical Ensemble and Its Universal Equivalence with the Microcanonical Ensemble Marius

More information

Global Optimization, Generalized Canonical Ensembles, and Universal Ensemble Equivalence

Global Optimization, Generalized Canonical Ensembles, and Universal Ensemble Equivalence Richard S. Ellis, Generalized Canonical Ensembles 1 Global Optimization, Generalized Canonical Ensembles, and Universal Ensemble Equivalence Richard S. Ellis Department of Mathematics and Statistics University

More information

Global Optimization, the Gaussian Ensemble, and Universal Ensemble Equivalence

Global Optimization, the Gaussian Ensemble, and Universal Ensemble Equivalence Richard S. Ellis, Gaussian Ensemble and Universal Ensemble Equivalence 1 Global Optimization, the Gaussian Ensemble, and Universal Ensemble Equivalence Richard S. Ellis Department of Mathematics and Statistics

More information

Lecture 32: Taylor Series and McLaurin series We saw last day that some functions are equal to a power series on part of their domain.

Lecture 32: Taylor Series and McLaurin series We saw last day that some functions are equal to a power series on part of their domain. Lecture 32: Taylor Series and McLaurin series We saw last day that some functions are equal to a power series on part of their domain. For example f(x) = 1 1 x = 1 + x + x2 + x 3 + = ln(1 + x) = x x2 2

More information

Chapter 11. Taylor Series. Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27

Chapter 11. Taylor Series. Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27 Chapter 11 Taylor Series Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27 First-Order Approximation We want to approximate function f by some simple function. Best possible approximation

More information

Microscopic Derivation of Ginzburg Landau Theory. Mathematics and Quantum Physics

Microscopic Derivation of Ginzburg Landau Theory. Mathematics and Quantum Physics Microscopic Derivation of Ginzburg Landau heory Robert Seiringer IS Austria Joint work with Rupert Frank, Christian Hainzl, and Jan Philip Solovej J. Amer. Math. Soc. 25 (2012), no. 3, 667 713 Mathematics

More information

Ferromagnets and superconductors. Kay Kirkpatrick, UIUC

Ferromagnets and superconductors. Kay Kirkpatrick, UIUC Ferromagnets and superconductors Kay Kirkpatrick, UIUC Ferromagnet and superconductor models: Phase transitions and asymptotics Kay Kirkpatrick, Urbana-Champaign October 2012 Ferromagnet and superconductor

More information

Laplace s Equation. Chapter Mean Value Formulas

Laplace s Equation. Chapter Mean Value Formulas Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic

More information

1. Prove the following properties satisfied by the gamma function: 4 n n!

1. Prove the following properties satisfied by the gamma function: 4 n n! Math 205A: Complex Analysis, Winter 208 Homework Problem Set #6 February 20, 208. Prove the following properties satisfied by the gamma function: (a) Values at half-integers: Γ ( n + 2 (b) The duplication

More information

Real Analysis Problems

Real Analysis Problems Real Analysis Problems Cristian E. Gutiérrez September 14, 29 1 1 CONTINUITY 1 Continuity Problem 1.1 Let r n be the sequence of rational numbers and Prove that f(x) = 1. f is continuous on the irrationals.

More information

P653 HW4. Due Sept. 22, Problem 1. Mean Field Dynamics

P653 HW4. Due Sept. 22, Problem 1. Mean Field Dynamics P653 HW4 Due Sept. 22, 2005 Problem 1. Mean Field Dynamics Here we wish to explore a crude, yet generic, model for the time dynamics of the order parameter near a critical point. For concreteness, lets

More information

H = 1 2 τψ gψ4, (0.0.1)

H = 1 2 τψ gψ4, (0.0.1) 32.1 Landau theory We have derived the Ginzburg-Landau Hamiltonian in Lecture 18. If we apply the mean field theory, that is, if we ignore fluctuations, then H = 1 2 τψ2 + 1 4 gψ4, (0.0.1) should be interpreted

More information

General Power Series

General Power Series General Power Series James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University March 29, 2018 Outline Power Series Consequences With all these preliminaries

More information

Lecture 4: Fourier Transforms.

Lecture 4: Fourier Transforms. 1 Definition. Lecture 4: Fourier Transforms. We now come to Fourier transforms, which we give in the form of a definition. First we define the spaces L 1 () and L 2 (). Definition 1.1 The space L 1 ()

More information

X n D X lim n F n (x) = F (x) for all x C F. lim n F n(u) = F (u) for all u C F. (2)

X n D X lim n F n (x) = F (x) for all x C F. lim n F n(u) = F (u) for all u C F. (2) 14:17 11/16/2 TOPIC. Convergence in distribution and related notions. This section studies the notion of the so-called convergence in distribution of real random variables. This is the kind of convergence

More information

9.1 System in contact with a heat reservoir

9.1 System in contact with a heat reservoir Chapter 9 Canonical ensemble 9. System in contact with a heat reservoir We consider a small system A characterized by E, V and N in thermal interaction with a heat reservoir A 2 characterized by E 2, V

More information

MADHAVA MATHEMATICS COMPETITION, December 2015 Solutions and Scheme of Marking

MADHAVA MATHEMATICS COMPETITION, December 2015 Solutions and Scheme of Marking MADHAVA MATHEMATICS COMPETITION, December 05 Solutions and Scheme of Marking NB: Part I carries 0 marks, Part II carries 30 marks and Part III carries 50 marks Part I NB Each question in Part I carries

More information

AP Calculus Testbank (Chapter 9) (Mr. Surowski)

AP Calculus Testbank (Chapter 9) (Mr. Surowski) AP Calculus Testbank (Chapter 9) (Mr. Surowski) Part I. Multiple-Choice Questions n 1 1. The series will converge, provided that n 1+p + n + 1 (A) p > 1 (B) p > 2 (C) p >.5 (D) p 0 2. The series

More information

The Gamma Function. July 9, Louisiana State University SMILE REU. The Gamma Function. N. Cannady, T. Ngo, A. Williamson.

The Gamma Function. July 9, Louisiana State University SMILE REU. The Gamma Function. N. Cannady, T. Ngo, A. Williamson. The The Louisiana State University SMILE REU July 9, 2010 The Developed as the unique extension of the factorial to non-integral values. The Developed as the unique extension of the factorial to non-integral

More information

2 The Curie Weiss Model

2 The Curie Weiss Model 2 The Curie Weiss Model In statistical mechanics, a mean-field approximation is often used to approximate a model by a simpler one, whose global behavior can be studied with the help of explicit computations.

More information

the convolution of f and g) given by

the convolution of f and g) given by 09:53 /5/2000 TOPIC Characteristic functions, cont d This lecture develops an inversion formula for recovering the density of a smooth random variable X from its characteristic function, and uses that

More information

This ODE arises in many physical systems that we shall investigate. + ( + 1)u = 0. (λ + s)x λ + s + ( + 1) a λ. (s + 1)(s + 2) a 0

This ODE arises in many physical systems that we shall investigate. + ( + 1)u = 0. (λ + s)x λ + s + ( + 1) a λ. (s + 1)(s + 2) a 0 Legendre equation This ODE arises in many physical systems that we shall investigate We choose We then have Substitution gives ( x 2 ) d 2 u du 2x 2 dx dx + ( + )u u x s a λ x λ a du dx λ a λ (λ + s)x

More information

Assignment 3. Tyler Shendruk February 26, 2010

Assignment 3. Tyler Shendruk February 26, 2010 Assignment 3 Tyler Shendruk February 6, 00 Kadar Ch. 4 Problem 7 N diatomic molecules are stuck on metal. The particles have three states:. in plane and aligned with the ˆx axis. in plane but aligned with

More information

Economics 204. The Transversality Theorem is a particularly convenient formulation of Sard s Theorem for our purposes: with r 1+max{0,n m}

Economics 204. The Transversality Theorem is a particularly convenient formulation of Sard s Theorem for our purposes: with r 1+max{0,n m} Economics 204 Lecture 13 Wednesday, August 12, 2009 Section 5.5 (Cont.) Transversality Theorem The Transversality Theorem is a particularly convenient formulation of Sard s Theorem for our purposes: Theorem

More information

Reflected Brownian Motion

Reflected Brownian Motion Chapter 6 Reflected Brownian Motion Often we encounter Diffusions in regions with boundary. If the process can reach the boundary from the interior in finite time with positive probability we need to decide

More information

AB.Q103.NOTES: Chapter 2.4, 3.1, 3.2 LESSON 1. Discovering the derivative at x = a: Slopes of secants and tangents to a curve

AB.Q103.NOTES: Chapter 2.4, 3.1, 3.2 LESSON 1. Discovering the derivative at x = a: Slopes of secants and tangents to a curve AB.Q103.NOTES: Chapter 2.4, 3.1, 3.2 LESSON 1 Discovering the derivative at x = a: Slopes of secants and tangents to a curve 1 1. Instantaneous rate of change versus average rate of change Equation of

More information

Legendre-Fenchel transforms in a nutshell

Legendre-Fenchel transforms in a nutshell 1 2 3 Legendre-Fenchel transforms in a nutshell Hugo Touchette School of Mathematical Sciences, Queen Mary, University of London, London E1 4NS, UK Started: July 11, 2005; last compiled: August 14, 2007

More information

4 Invariant Statistical Decision Problems

4 Invariant Statistical Decision Problems 4 Invariant Statistical Decision Problems 4.1 Invariant decision problems Let G be a group of measurable transformations from the sample space X into itself. The group operation is composition. Note that

More information

Oblique derivative problems for elliptic and parabolic equations, Lecture II

Oblique derivative problems for elliptic and parabolic equations, Lecture II of the for elliptic and parabolic equations, Lecture II Iowa State University July 22, 2011 of the 1 2 of the of the As a preliminary step in our further we now look at a special situation for elliptic.

More information

Thermodynamic versus statistical nonequivalence of ensembles for the mean-field Blume-Emery- Griffiths model

Thermodynamic versus statistical nonequivalence of ensembles for the mean-field Blume-Emery- Griffiths model University of Massachusetts Amherst ScholarWorks@UMass Amherst Mathematics and Statistics Department Faculty Publication Series Mathematics and Statistics 2004 Thermodynamic versus statistical nonequivalence

More information

The Laplace Transform. Background: Improper Integrals

The Laplace Transform. Background: Improper Integrals The Laplace Transform Background: Improper Integrals Recall: Definite Integral: a, b real numbers, a b; f continuous on [a, b] b a f(x) dx 1 Improper integrals: Type I Infinite interval of integration

More information

A Ginzburg-Landau Type Problem for Nematics with Highly Anisotropic Elastic Term

A Ginzburg-Landau Type Problem for Nematics with Highly Anisotropic Elastic Term A Ginzburg-Landau Type Problem for Nematics with Highly Anisotropic Elastic Term Peter Sternberg In collaboration with Dmitry Golovaty (Akron) and Raghav Venkatraman (Indiana) Department of Mathematics

More information

Ferromagnets and the classical Heisenberg model. Kay Kirkpatrick, UIUC

Ferromagnets and the classical Heisenberg model. Kay Kirkpatrick, UIUC Ferromagnets and the classical Heisenberg model Kay Kirkpatrick, UIUC Ferromagnets and the classical Heisenberg model: asymptotics for a mean-field phase transition Kay Kirkpatrick, Urbana-Champaign June

More information

First-Passage Statistics of Extreme Values

First-Passage Statistics of Extreme Values First-Passage Statistics of Extreme Values Eli Ben-Naim Los Alamos National Laboratory with: Paul Krapivsky (Boston University) Nathan Lemons (Los Alamos) Pearson Miller (Yale, MIT) Talk, publications

More information

Mesoscale Simulation Methods. Ronojoy Adhikari The Institute of Mathematical Sciences Chennai

Mesoscale Simulation Methods. Ronojoy Adhikari The Institute of Mathematical Sciences Chennai Mesoscale Simulation Methods Ronojoy Adhikari The Institute of Mathematical Sciences Chennai Outline What is mesoscale? Mesoscale statics and dynamics through coarse-graining. Coarse-grained equations

More information

Thermodynamic limit for a system of interacting fermions in a random medium. Pieces one-dimensional model

Thermodynamic limit for a system of interacting fermions in a random medium. Pieces one-dimensional model Thermodynamic limit for a system of interacting fermions in a random medium. Pieces one-dimensional model Nikolaj Veniaminov (in collaboration with Frédéric Klopp) CEREMADE, University of Paris IX Dauphine

More information

Read carefully the instructions on the answer book and make sure that the particulars required are entered on each answer book.

Read carefully the instructions on the answer book and make sure that the particulars required are entered on each answer book. THE UNIVERSITY OF WARWICK FOURTH YEAR EXAMINATION: APRIL 2007 INTRODUCTION TO STATISTICAL MECHANICS Time Allowed: 3 hours Read carefully the instructions on the answer book and make sure that the particulars

More information

Nonlinear Integral Equation Formulation of Orthogonal Polynomials

Nonlinear Integral Equation Formulation of Orthogonal Polynomials Nonlinear Integral Equation Formulation of Orthogonal Polynomials Eli Ben-Naim Theory Division, Los Alamos National Laboratory with: Carl Bender (Washington University, St. Louis) C.M. Bender and E. Ben-Naim,

More information

Analysis Qualifying Exam

Analysis Qualifying Exam Analysis Qualifying Exam Spring 2017 Problem 1: Let f be differentiable on R. Suppose that there exists M > 0 such that f(k) M for each integer k, and f (x) M for all x R. Show that f is bounded, i.e.,

More information

Sequences and Series of Functions

Sequences and Series of Functions Chapter 13 Sequences and Series of Functions These notes are based on the notes A Teacher s Guide to Calculus by Dr. Louis Talman. The treatment of power series that we find in most of today s elementary

More information

In class midterm Exam - Answer key

In class midterm Exam - Answer key Fall 2013 In class midterm Exam - Answer key ARE211 Problem 1 (20 points). Metrics: Let B be the set of all sequences x = (x 1,x 2,...). Define d(x,y) = sup{ x i y i : i = 1,2,...}. a) Prove that d is

More information

The polar coordinates

The polar coordinates The polar coordinates 1 2 3 4 Graphing in polar coordinates 5 6 7 8 Area and length in polar coordinates 9 10 11 Partial deravitive 12 13 14 15 16 17 18 19 20 Double Integral 21 22 23 24 25 26 27 Triple

More information

Collective Effects. Equilibrium and Nonequilibrium Physics

Collective Effects. Equilibrium and Nonequilibrium Physics Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 3, 3 March 2006 Collective Effects in Equilibrium and Nonequilibrium Physics Website: http://cncs.bnu.edu.cn/mccross/course/ Caltech

More information

F (x) = P [X x[. DF1 F is nondecreasing. DF2 F is right-continuous

F (x) = P [X x[. DF1 F is nondecreasing. DF2 F is right-continuous 7: /4/ TOPIC Distribution functions their inverses This section develops properties of probability distribution functions their inverses Two main topics are the so-called probability integral transformation

More information

Long Run Average Cost Problem. E x β j c(x j ; u j ). j=0. 1 x c(x j ; u j ). n n E

Long Run Average Cost Problem. E x β j c(x j ; u j ). j=0. 1 x c(x j ; u j ). n n E Chapter 6. Long Run Average Cost Problem Let {X n } be an MCP with state space S and control set U(x) forx S. To fix notation, let J β (x; {u n }) be the cost associated with the infinite horizon β-discount

More information

Economics 204 Summer/Fall 2011 Lecture 5 Friday July 29, 2011

Economics 204 Summer/Fall 2011 Lecture 5 Friday July 29, 2011 Economics 204 Summer/Fall 2011 Lecture 5 Friday July 29, 2011 Section 2.6 (cont.) Properties of Real Functions Here we first study properties of functions from R to R, making use of the additional structure

More information

e x = 1 + x + x2 2! + x3 If the function f(x) can be written as a power series on an interval I, then the power series is of the form

e x = 1 + x + x2 2! + x3 If the function f(x) can be written as a power series on an interval I, then the power series is of the form Taylor Series Given a function f(x), we would like to be able to find a power series that represents the function. For example, in the last section we noted that we can represent e x by the power series

More information

Copyright 2010 Pearson Education, Inc. Publishing as Prentice Hall.

Copyright 2010 Pearson Education, Inc. Publishing as Prentice Hall. .1 Limits of Sequences. CHAPTER.1.0. a) True. If converges, then there is an M > 0 such that M. Choose by Archimedes an N N such that N > M/ε. Then n N implies /n M/n M/N < ε. b) False. = n does not converge,

More information

Generalized canonical ensembles and ensemble equivalence

Generalized canonical ensembles and ensemble equivalence University of Massachusetts Amherst ScholarWorks@UMass Amherst Mathematics and Statistics Department Faculty Publication Series Mathematics and Statistics 2006 Generalized canonical ensembles and ensemble

More information

Nonlinear Systems and Control Lecture # 19 Perturbed Systems & Input-to-State Stability

Nonlinear Systems and Control Lecture # 19 Perturbed Systems & Input-to-State Stability p. 1/1 Nonlinear Systems and Control Lecture # 19 Perturbed Systems & Input-to-State Stability p. 2/1 Perturbed Systems: Nonvanishing Perturbation Nominal System: Perturbed System: ẋ = f(x), f(0) = 0 ẋ

More information

The Ising model Summary of L12

The Ising model Summary of L12 The Ising model Summary of L2 Aim: Study connections between macroscopic phenomena and the underlying microscopic world for a ferromagnet. How: Study the simplest possible model of a ferromagnet containing

More information

Math Homework 2

Math Homework 2 Math 73 Homework Due: September 8, 6 Suppose that f is holomorphic in a region Ω, ie an open connected set Prove that in any of the following cases (a) R(f) is constant; (b) I(f) is constant; (c) f is

More information

MA 162 FINAL EXAM PRACTICE PROBLEMS Spring Find the angle between the vectors v = 2i + 2j + k and w = 2i + 2j k. C.

MA 162 FINAL EXAM PRACTICE PROBLEMS Spring Find the angle between the vectors v = 2i + 2j + k and w = 2i + 2j k. C. MA 6 FINAL EXAM PRACTICE PROBLEMS Spring. Find the angle between the vectors v = i + j + k and w = i + j k. cos 8 cos 5 cos D. cos 7 E. cos. Find a such that u = i j + ak and v = i + j + k are perpendicular.

More information

INTRODUCTION TO STATISTICAL MECHANICS Exercises

INTRODUCTION TO STATISTICAL MECHANICS Exercises École d hiver de Probabilités Semaine 1 : Mécanique statistique de l équilibre CIRM - Marseille 4-8 Février 2012 INTRODUCTION TO STATISTICAL MECHANICS Exercises Course : Yvan Velenik, Exercises : Loren

More information

Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( )

Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( ) Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio (2014-2015) Etienne Tanré - Olivier Faugeras INRIA - Team Tosca October 22nd, 2014 E. Tanré (INRIA - Team Tosca) Mathematical

More information

Lecture Characterization of Infinitely Divisible Distributions

Lecture Characterization of Infinitely Divisible Distributions Lecture 10 1 Characterization of Infinitely Divisible Distributions We have shown that a distribution µ is infinitely divisible if and only if it is the weak limit of S n := X n,1 + + X n,n for a uniformly

More information

Real Analysis Comprehensive Exam Fall A(k, ε) is of Lebesgue measure zero.

Real Analysis Comprehensive Exam Fall A(k, ε) is of Lebesgue measure zero. Real Analysis Comprehensive Exam Fall 2002 by XYC Good luck! [1] For ε>0andk>0, denote by A(k, ε) thesetofx Rsuch that x p q 1 for any integers p, q with q 0. k q 2+ε Show that R \ k=1 A(k, ε) is of Lebesgue

More information

5 Operations on Multiple Random Variables

5 Operations on Multiple Random Variables EE360 Random Signal analysis Chapter 5: Operations on Multiple Random Variables 5 Operations on Multiple Random Variables Expected value of a function of r.v. s Two r.v. s: ḡ = E[g(X, Y )] = g(x, y)f X,Y

More information

(x 3)(x + 5) = (x 3)(x 1) = x + 5. sin 2 x e ax bx 1 = 1 2. lim

(x 3)(x + 5) = (x 3)(x 1) = x + 5. sin 2 x e ax bx 1 = 1 2. lim SMT Calculus Test Solutions February, x + x 5 Compute x x x + Answer: Solution: Note that x + x 5 x x + x )x + 5) = x )x ) = x + 5 x x + 5 Then x x = + 5 = Compute all real values of b such that, for fx)

More information

Ginzburg-Landau Theory of Phase Transitions

Ginzburg-Landau Theory of Phase Transitions Subedi 1 Alaska Subedi Prof. Siopsis Physics 611 Dec 5, 008 Ginzburg-Landau Theory of Phase Transitions 1 Phase Transitions A phase transition is said to happen when a system changes its phase. The physical

More information

Estimates for the resolvent and spectral gaps for non self-adjoint operators. University Bordeaux

Estimates for the resolvent and spectral gaps for non self-adjoint operators. University Bordeaux for the resolvent and spectral gaps for non self-adjoint operators 1 / 29 Estimates for the resolvent and spectral gaps for non self-adjoint operators Vesselin Petkov University Bordeaux Mathematics Days

More information

1 The fundamental equation of equilibrium statistical mechanics. 3 General overview on the method of ensembles 10

1 The fundamental equation of equilibrium statistical mechanics. 3 General overview on the method of ensembles 10 Contents EQUILIBRIUM STATISTICAL MECHANICS 1 The fundamental equation of equilibrium statistical mechanics 2 2 Conjugate representations 6 3 General overview on the method of ensembles 10 4 The relation

More information

Gärtner-Ellis Theorem and applications.

Gärtner-Ellis Theorem and applications. Gärtner-Ellis Theorem and applications. Elena Kosygina July 25, 208 In this lecture we turn to the non-i.i.d. case and discuss Gärtner-Ellis theorem. As an application, we study Curie-Weiss model with

More information

Physics 212: Statistical mechanics II Lecture XI

Physics 212: Statistical mechanics II Lecture XI Physics 212: Statistical mechanics II Lecture XI The main result of the last lecture was a calculation of the averaged magnetization in mean-field theory in Fourier space when the spin at the origin is

More information

Econ 508B: Lecture 5

Econ 508B: Lecture 5 Econ 508B: Lecture 5 Expectation, MGF and CGF Hongyi Liu Washington University in St. Louis July 31, 2017 Hongyi Liu (Washington University in St. Louis) Math Camp 2017 Stats July 31, 2017 1 / 23 Outline

More information

PUTNAM PROBLEMS SEQUENCES, SERIES AND RECURRENCES. Notes

PUTNAM PROBLEMS SEQUENCES, SERIES AND RECURRENCES. Notes PUTNAM PROBLEMS SEQUENCES, SERIES AND RECURRENCES Notes. x n+ = ax n has the general solution x n = x a n. 2. x n+ = x n + b has the general solution x n = x + (n )b. 3. x n+ = ax n + b (with a ) can be

More information

Single-scaling-field approach for an isolated polymer chain

Single-scaling-field approach for an isolated polymer chain J. Phys. A: Math. Gen. 14 (1981) L55-L61. Printed in Great Britain LETTER TO THE EDITOR Single-scaling-field approach for an isolated polymer chain S Redner and P J Reynoldst Center for Polymer Studies2

More information

LEGENDRE POLYNOMIALS AND APPLICATIONS. We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.

LEGENDRE POLYNOMIALS AND APPLICATIONS. We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates. LEGENDRE POLYNOMIALS AND APPLICATIONS We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.. Legendre equation: series solutions The Legendre equation is

More information

International Competition in Mathematics for Universtiy Students in Plovdiv, Bulgaria 1994

International Competition in Mathematics for Universtiy Students in Plovdiv, Bulgaria 1994 International Competition in Mathematics for Universtiy Students in Plovdiv, Bulgaria 1994 1 PROBLEMS AND SOLUTIONS First day July 29, 1994 Problem 1. 13 points a Let A be a n n, n 2, symmetric, invertible

More information

The Γ-Limit of the Ginzburg-Landau energy in a Thin Domain with a Large Magnetic Field

The Γ-Limit of the Ginzburg-Landau energy in a Thin Domain with a Large Magnetic Field p. 1/2 The Γ-Limit of the Ginzburg-Landau energy in a Thin Domain with a Large Magnetic Field Tien-Tsan Shieh tshieh@math.arizona.edu Department of Mathematics University of Arizona Tucson, AZ 85719 p.

More information

Gases and the Virial Expansion

Gases and the Virial Expansion Gases and the irial Expansion February 7, 3 First task is to examine what ensemble theory tells us about simple systems via the thermodynamic connection Calculate thermodynamic quantities: average energy,

More information

Collective Effects. Equilibrium and Nonequilibrium Physics

Collective Effects. Equilibrium and Nonequilibrium Physics 1 Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 2, 24 March 2006 1 Collective Effects in Equilibrium and Nonequilibrium Physics Website: http://cncs.bnu.edu.cn/mccross/course/ Caltech

More information

Lecture 5: Expectation

Lecture 5: Expectation Lecture 5: Expectation 1. Expectations for random variables 1.1 Expectations for simple random variables 1.2 Expectations for bounded random variables 1.3 Expectations for general random variables 1.4

More information

2. Function spaces and approximation

2. Function spaces and approximation 2.1 2. Function spaces and approximation 2.1. The space of test functions. Notation and prerequisites are collected in Appendix A. Let Ω be an open subset of R n. The space C0 (Ω), consisting of the C

More information

Real Analysis Math 131AH Rudin, Chapter #1. Dominique Abdi

Real Analysis Math 131AH Rudin, Chapter #1. Dominique Abdi Real Analysis Math 3AH Rudin, Chapter # Dominique Abdi.. If r is rational (r 0) and x is irrational, prove that r + x and rx are irrational. Solution. Assume the contrary, that r+x and rx are rational.

More information

Chapter 11 - Sequences and Series

Chapter 11 - Sequences and Series Calculus and Analytic Geometry II Chapter - Sequences and Series. Sequences Definition. A sequence is a list of numbers written in a definite order, We call a n the general term of the sequence. {a, a

More information

Phase Transitions and Critical Behavior:

Phase Transitions and Critical Behavior: II Phase Transitions and Critical Behavior: A. Phenomenology (ibid., Chapter 10) B. mean field theory (ibid., Chapter 11) C. Failure of MFT D. Phenomenology Again (ibid., Chapter 12) // Windsor Lectures

More information

Renyi divergence and asymptotic theory of minimal K-point random graphs Alfred Hero Dept. of EECS, The University of Michigan, Summer 1999

Renyi divergence and asymptotic theory of minimal K-point random graphs Alfred Hero Dept. of EECS, The University of Michigan, Summer 1999 Renyi divergence and asymptotic theory of minimal K-point random graphs Alfred Hero Dept. of EECS, The University of Michigan, Summer 999 Outline Rényi Entropy and Rényi Divergence ffl Euclidean k-minimal

More information

A List of Problems in Real Analysis

A List of Problems in Real Analysis A List of Problems in Real Analysis W.Yessen & T.Ma December 3, 218 This document was first created by Will Yessen, who was a graduate student at UCI. Timmy Ma, who was also a graduate student at UCI,

More information

Monte Carlo study of the Baxter-Wu model

Monte Carlo study of the Baxter-Wu model Monte Carlo study of the Baxter-Wu model Nir Schreiber and Dr. Joan Adler Monte Carlo study of the Baxter-Wu model p.1/40 Outline Theory of phase transitions, Monte Carlo simulations and finite size scaling

More information

1. Thermodynamics 1.1. A macroscopic view of matter

1. Thermodynamics 1.1. A macroscopic view of matter 1. Thermodynamics 1.1. A macroscopic view of matter Intensive: independent of the amount of substance, e.g. temperature,pressure. Extensive: depends on the amount of substance, e.g. internal energy, enthalpy.

More information

The Laplace Transform. Background: Improper Integrals

The Laplace Transform. Background: Improper Integrals The Laplace Transform Background: Improper Integrals Recall: Definite Integral: a, b real numbers, a b; f continuous on [a, b] b a f(x) dx 1 Improper integrals: Type I Infinite interval of integration

More information

LECTURE 10: REVIEW OF POWER SERIES. 1. Motivation

LECTURE 10: REVIEW OF POWER SERIES. 1. Motivation LECTURE 10: REVIEW OF POWER SERIES By definition, a power series centered at x 0 is a series of the form where a 0, a 1,... and x 0 are constants. For convenience, we shall mostly be concerned with the

More information

Math Final Exam Review

Math Final Exam Review Math - Final Exam Review. Find dx x + 6x +. Name: Solution: We complete the square to see if this function has a nice form. Note we have: x + 6x + (x + + dx x + 6x + dx (x + + Note that this looks a lot

More information

Math 564 Homework 1. Solutions.

Math 564 Homework 1. Solutions. Math 564 Homework 1. Solutions. Problem 1. Prove Proposition 0.2.2. A guide to this problem: start with the open set S = (a, b), for example. First assume that a >, and show that the number a has the properties

More information

Invaded cluster dynamics for frustrated models

Invaded cluster dynamics for frustrated models PHYSICAL REVIEW E VOLUME 57, NUMBER 1 JANUARY 1998 Invaded cluster dynamics for frustrated models Giancarlo Franzese, 1, * Vittorio Cataudella, 1, * and Antonio Coniglio 1,2, * 1 INFM, Unità di Napoli,

More information

Behaviour of Lipschitz functions on negligible sets. Non-differentiability in R. Outline

Behaviour of Lipschitz functions on negligible sets. Non-differentiability in R. Outline Behaviour of Lipschitz functions on negligible sets G. Alberti 1 M. Csörnyei 2 D. Preiss 3 1 Università di Pisa 2 University College London 3 University of Warwick Lars Ahlfors Centennial Celebration Helsinki,

More information

Definitions & Theorems

Definitions & Theorems Definitions & Theorems Math 147, Fall 2009 December 19, 2010 Contents 1 Logic 2 1.1 Sets.................................................. 2 1.2 The Peano axioms..........................................

More information

1 Lecture 25: Extreme values

1 Lecture 25: Extreme values 1 Lecture 25: Extreme values 1.1 Outline Absolute maximum and minimum. Existence on closed, bounded intervals. Local extrema, critical points, Fermat s theorem Extreme values on a closed interval Rolle

More information

Principles of Equilibrium Statistical Mechanics

Principles of Equilibrium Statistical Mechanics Debashish Chowdhury, Dietrich Stauffer Principles of Equilibrium Statistical Mechanics WILEY-VCH Weinheim New York Chichester Brisbane Singapore Toronto Table of Contents Part I: THERMOSTATICS 1 1 BASIC

More information

MTH4101 CALCULUS II REVISION NOTES. 1. COMPLEX NUMBERS (Thomas Appendix 7 + lecture notes) ax 2 + bx + c = 0. x = b ± b 2 4ac 2a. i = 1.

MTH4101 CALCULUS II REVISION NOTES. 1. COMPLEX NUMBERS (Thomas Appendix 7 + lecture notes) ax 2 + bx + c = 0. x = b ± b 2 4ac 2a. i = 1. MTH4101 CALCULUS II REVISION NOTES 1. COMPLEX NUMBERS (Thomas Appendix 7 + lecture notes) 1.1 Introduction Types of numbers (natural, integers, rationals, reals) The need to solve quadratic equations:

More information

Magnetism at finite temperature: molecular field, phase transitions

Magnetism at finite temperature: molecular field, phase transitions Magnetism at finite temperature: molecular field, phase transitions -The Heisenberg model in molecular field approximation: ferro, antiferromagnetism. Ordering temperature; thermodynamics - Mean field

More information

A mathematical model for a copolymer in an emulsion

A mathematical model for a copolymer in an emulsion J Math Chem (2010) 48:83 94 DOI 10.1007/s10910-009-9564-y ORIGINAL PAPER A mathematical model for a copolymer in an emulsion F. den Hollander N. Pétrélis Received: 3 June 2007 / Accepted: 22 April 2009

More information

Clusters and Percolation

Clusters and Percolation Chapter 6 Clusters and Percolation c 2012 by W. Klein, Harvey Gould, and Jan Tobochnik 5 November 2012 6.1 Introduction In this chapter we continue our investigation of nucleation near the spinodal. We

More information