Lecture 13 Bragg-Williams Theory

Similar documents
Time Domain Method of Moments

Lecture 15: Phase Transitions. Phase transitions

8.333: Statistical Mechanics I Problem Set # 4 Due: 11/13/13 Non-interacting particles

Lightning electromagnetic environment in the presence of a tall grounded strike object

Lecture 15 (Nov. 1, 2017)

Chapter 9. The excitation process

Hankel Optimal Model Order Reduction 1

MOLECULAR ORBITAL THEORY- PART I

Strauss PDEs 2e: Section Exercise 3 Page 1 of 13. u tt c 2 u xx = cos x. ( 2 t c 2 2 x)u = cos x. v = ( t c x )u

Determination of the reaction order

Chapter 2 Linear Elastic Fracture Mechanics

Analysis of discretization in the direct simulation Monte Carlo

Chemistry (Physical chemistry) Lecture 10.

The gravitational phenomena without the curved spacetime

ELECTROMAGNETIC WAVES

Numerical Tests of Nucleation Theories for the Ising Models. Abstract

LECTURE 2 Geometrical Properties of Rod Cross Sections (Part 2) 1 Moments of Inertia Transformation with Parallel Transfer of Axes.

Math 220A - Fall 2002 Homework 8 Solutions

Lecture 10: Persistent Random Walks and the Telegrapher s Equation

Quantum critical phenomena in the Schrodinger formulation: mapping to classical lattices

CRITICAL EXPONENTS TAKING INTO ACCOUNT DYNAMIC SCALING FOR ADSORPTION ON SMALL-SIZE ONE-DIMENSIONAL CLUSTERS

Wave Propagation through Random Media

Advanced Computational Fluid Dynamics AA215A Lecture 4

Electromagnetic radiation of the travelling spin wave propagating in an antiferromagnetic plate. Exact solution.

max min z i i=1 x j k s.t. j=1 x j j:i T j

Physica A. Calculation of free energy of a three-dimensional Ising-like system in an external field with the use of the ρ 6 model

STATISTICAL MECHANICS & THERMODYNAMICS

Probabilistic Graphical Models

A NORMALIZED EQUATION OF AXIALLY LOADED PILES IN ELASTO-PLASTIC SOIL

On Certain Singular Integral Equations Arising in the Analysis of Wellbore Recharge in Anisotropic Formations

Lecture 3 - Lorentz Transformations

Phys 561 Classical Electrodynamics. Midterm

A simple expression for radial distribution functions of pure fluids and mixtures

Study of EM waves in Periodic Structures (mathematical details)

Relativity in Classical Physics

Mass Transfer (Stoffaustausch) Fall 2012

Chapter 2: Solution of First order ODE

CMSC 451: Lecture 9 Greedy Approximation: Set Cover Thursday, Sep 28, 2017

Solving Right Triangles Using Trigonometry Examples

University of Washington Department of Chemistry Chemistry 452/456 Summer Quarter 2011

The Concept of Mass as Interfering Photons, and the Originating Mechanism of Gravitation D.T. Froedge

Lecture Notes 4 MORE DYNAMICS OF NEWTONIAN COSMOLOGY

Department of Mechanical Engineering

Entropy, Dynamics and Instantaneous Normal Modes in a Random Energy Model. T. Keyes

ON DIRICHLET PROBLEM WITH SINGULAR INTEGRAL BOUNDARY CONDITION

Maximum Entropy and Exponential Families

Wavetech, LLC. Ultrafast Pulses and GVD. John O Hara Created: Dec. 6, 2013

SURFACE WAVES OF NON-RAYLEIGH TYPE

MODELING MATTER AT NANOSCALES. 4. Introduction to quantum treatments Eigenvectors and eigenvalues of a matrix

Magnetic Ordering in a Quasi Two-Dimensional Planar~Ferrdmagnet:CoCl 2 in Graphite. H. Zabel Ruhr-University Bochum, Bochum,FRG

V. Interacting Particles

RC DEEP BEAMS ANALYSIS CONSIDERING LOCALIZATION IN COMPRESSION

We will show that: that sends the element in π 1 (P {z 1, z 2, z 3, z 4 }) represented by l j to g j G.

ELECTROMAGNETIC NORMAL MODES AND DISPERSION FORCES.

EXACT TRAVELLING WAVE SOLUTIONS FOR THE GENERALIZED KURAMOTO-SIVASHINSKY EQUATION

Aharonov-Bohm effect. Dan Solomon.

Subject: Introduction to Component Matching and Off-Design Operation % % ( (1) R T % (

The Hanging Chain. John McCuan. January 19, 2006

Heat exchangers: Heat exchanger types:

Special Relativity Simply Debunked in Five Steps!

Current density and forces for a current loop moving parallel over a thin conducting sheet

Strategy selection for arbitrary initial condition

Нацiональна академiя наук України IНСТИТУТ ФIЗИКИ КОНДЕНСОВАНИХ СИСТЕМ ЛЬВIВ. УДК: PACS: b

arxiv:cond-mat/ v1 [cond-mat.str-el] 3 Aug 2006

A population of 50 flies is expected to double every week, leading to a function of the x

arxiv:physics/ v1 14 May 2002

Brazilian Journal of Physics, vol. 29, no. 3, September, Classical and Quantum Mechanics of a Charged Particle

Case I: 2 users In case of 2 users, the probability of error for user 1 was earlier derived to be 2 A1

ME 354 Tutorial, Week#5 Refrigeration Cycle A large refrigeration plant is to be maintained at -15 o C, and it requires refrigeration at a rate of

Casimir self-energy of a free electron

arxiv:gr-qc/ v2 6 Feb 2004

Physics of Relaxation. Outline

Physics 486. Classical Newton s laws Motion of bodies described in terms of initial conditions by specifying x(t), v(t).

ABSTRACT. and since the phase transition of asymmetric one-component fluid or mixture

Remark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the function V ( x ) to be positive definite.

Most results in this section are stated without proof.

EECS 120 Signals & Systems University of California, Berkeley: Fall 2005 Gastpar November 16, Solutions to Exam 2

4 Puck s action plane fracture criteria

Discrete Bessel functions and partial difference equations

Four-dimensional equation of motion for viscous compressible substance with regard to the acceleration field, pressure field and dissipation field

What are the locations of excess energy in open channels?

11.4 Molecular Orbital Description of the Hydrogen Molecule Electron Configurations of Homonuclear Diatomic Molecules

MTH 142 Solution Practice for Exam 2

Discontinuity of Shannon information entropy for two-electron atoms

Chapter 8 Hypothesis Testing

Advances in Radio Science

Diagrammatic Monte Carlo simulation of quantum impurity models weak and strong coupling approach

Thermodynamic Properties of Supercritical Fluids: Example of n-hexane

Diagrammatic Monte Carlo simulation of quantum impurity models weak and strong coupling approach

Where as discussed previously we interpret solutions to this partial differential equation in the weak sense: b

22.54 Neutron Interactions and Applications (Spring 2004) Chapter 6 (2/24/04) Energy Transfer Kernel F(E E')

+Ze. n = N/V = 6.02 x x (Z Z c ) m /A, (1.1) Avogadro s number

Spectroscopic Signatures of Novel Oxygen-Defect Complexes in Stoichiometrically Controlled CdSe

MODELLING THE POSTPEAK STRESS DISPLACEMENT RELATIONSHIP OF CONCRETE IN UNIAXIAL COMPRESSION

7 Max-Flow Problems. Business Computing and Operations Research 608

Integration of the Finite Toda Lattice with Complex-Valued Initial Data

A NETWORK SIMPLEX ALGORITHM FOR THE MINIMUM COST-BENEFIT NETWORK FLOW PROBLEM

Hidden Momentum in a Spinning Sphere

In this problem, we are given the following quantities: We want to find: Equations and basic calculations:

depend on the renormalization sheme) only near the ritial point at zero oupling while at intermediate and strong ouplings this statement is still ques

Transcription:

Leture 13 Bragg-Williams Theory As noted in Chapter 11, an alternative mean-field approah is to derive a free energy, F, in terms of our order parameter,m, and then minimize F with respet to m. We begin with a slightly different definition for m, m = ( ) where = + is the total number of sites in the lattie. (13.1) The entropy assoiated with a given m is just the logarithm of the total number of onfigurations with a given and, S = k B ln = k B ln! (13.) (1+m)! (1m)! Applying Stirling s approximation, the above equation redues to, (1 + m) (1 + m) (1 + m) S = k B [ ln { ln (1 m) (1 m) (1 m) + ln }] (1 + m) (1 + m) (1 + m) = k B [ ln { ln + ln (1 + m) (1 m) (1 m) (1 m) (1 m) + ln + ln }] [ ] (1 + m) (1 + m) (1 m) (1 m) = k B ln ln and hene, s(m) S = k B [ln 1 (1 + m) ln(1 + m) 1 (1 m) ln(1 m) ] (13.3)

LECTURE 13. BRAGG-WILLIAMS THEORY 87 whih is just the entropy of mixing. What is the mean energy, Ē? To alulate this average, we have to evaluate the Hamiltonian reognizing that the spins interat with eah other. In the mean field spirit, as we saw for the Weiss Moleular field theory, we replae σ i in the Hamiltonian with its position independent average, m: Ē = J ij m = 1 Jqm (13.4) where, again, q = d is the lattie oordination number or the number of nearest neighbour sites for a d-dimensional hyperubi lattie. Let us now write down the free energy for this system, we have: f(t, m) F (T, m) = E T S = 1 Jqm + 1 k BT [(1 + m) ln(1 + m) + (1 m) ln(1 m)] (13.5) k B T ln Expansion of s(m) near T ow, in the viinity of T, m is small (as we have noted before), so we may expand s(m), and hene f(m), as follows, s(m) = s(0) + 1 s m + 1 4 s m 4 +... (13.6)! m 4! m 4 where, and, s m 4 s m 4 ( ) k B = (1 + m)(1 m) = k B 1 m + m = k B = k (1 + m) 3 (1 m) 3 B Therefore, the entropy is given as, s(m) = k B [ln 1 m 1 ] 1 m4 +... and the free energy beomes, f(m) = k B T ln + ( 1 Jqm + k B T 1 ) m + k B T 1 1 m4 or (13.7) f(m) k B = 1 (T T )m + T m4 1 T ln (13.8)

LECTURE 13. BRAGG-WILLIAMS THEORY 88 It is lear from the above equation that for T > T, the minimum in the free energy ours for m = 0. For T < T, where does the minimum in f(m) our? f(m)/k m = (T T )m + m = 3(T T ) T 4T m3 1 = 0 or m = ± 1/ T 3 T 1 (13.9) whih is the same result as derived for the Weiss Moleular Field theory! To expliitly see the behaviour of the free energy, as a funtion of the order parameter m, in the viinity of T, refer to the following figure, h = 0 F T>T T=T T<T m 13.1 Spins in an external field Let us now onsider our system of spins in the presene of an external magneti field, h, suh that the free energy beomes f f hm. Referring to Eq. 13.5 we have, f(t, m) = 1 Jqm hm + 1 k BT [(1 + m) ln(1 + m) + (1 m) ln(1 m)] k B T ln Minimizing the above expression with respet to m we have, f m = Jqm h + 1 1 + m kt ln = 0 1 m (13.10)

LECTURE 13. BRAGG-WILLIAMS THEORY 89 Solving for h we find, h = Jqm + 1 1 + m kt ln 1 m = Jqm + 1 [ ] kt m + m3 3 + m5 5 +..., m < 1 = Jqm + kt tanh 1m, or [ h m = tanh kt + T ] T m (13.11) whih again is the same result that we obtained from the Weiss Moleular Field Theory model. 13. Critial Behaviour of Mean Field Theories Using the expression for the average energy, Eq. 13.4 with the mean field result, Eq. 13.9, we have, Ē Jq T T 3 T (13.1) From this expression, we may alulate the heat apaity as follows, C h T T = E = 3Jq 3k T T T T (13.13) Sine m = 0 for T > T, this implies the average energy and, hene, the heat apaity are zero as we approah T from the positive side. The qualitative behaviour of the heat apaity atually demonstrates a usp disontinuity, as seen below. This behaviour is in ontrast to more exat theories and experiments whih demonstrate a power law dependene of the heat apaity, C h = B ± T T α (13.14)

LECTURE 13. BRAGG-WILLIAMS THEORY 90 C v disontinuity at T T T