2
Examples of binary mixtures Mixtures may be solid or liquid Brass (85% copper-15% zinc) in solid form is an example of a solid binary mixture. When raised to above the respective melting temperatures, a liquid binary mixture is produced. It is in this phase that the components can be mixed. Example of a liquid mixture at room temperature. Gasohol (10% ethanol-90% gasoline)

3
What is Phase Separation? The process of phase separation is the unmixing of a thermodynamically unstable solution. Essentially the mixture self segregates based on thermodynamic principles.

4
States of a mixture Homogeneous state In the homogeneous state, the respective concentration of the components is the same at all locations. Heterogeneous state In the heterogeneous state, there are local domains richer in one of the components

5
Phase Diagrams A phase diagram is a plot with variables: temperature (T) and concentration (  Composition, , is defined as  (r) =  a (r) -  b (r).  i (r) is the relative concentration of each component at location r State of mixture is described by its location on the phase diagram and defined by ( , T)

6
States of a mixture related to the phase diagram Which state the system is in is determined by its location on the phase diagram. If the state of the mixture is above the coexistence curve, it is homogeneous. If the state of the mixture is below the coexistence cure, the system is globally unstable in the homogeneous state.

8
Features of phase diagram Coexistence curve separate regions of homogeneous and heterogeneous states. Along this curve the mixed and unmixed states are in equilibrium. When quenched below the coexistence curve, the system is locally stable in the homogeneous state. Between the coexistence and spinodal curves, the homogeneous state is said to be metastable. Beneath the spinodal curve, the system is unstable in the homogeneous state.

18
Goals To observe that the elongated domains grow as a function of time: size(t) ~t  Compute , by calculating the Fourier transform (structure factor) as a function of time, simulating an optical study of the system. Compute , by using the correlation function to calculate time dependence of the correlation length. Calculate the critical exponent and compare to known value of  = 1/3

20
Mathematical Background The kinetics of the systems is described by the continuity equation: The diffusive current is: M represents the mobility of the particles Chemical potential, , is defined as the functional derivative of the Free Energy with respect to particle number:

21
Mathematical Background Free Energy is described by the functional: f(r, t) is the Ginzburg-Landau function: The free energy becomes:

22
Mathematical Background Using the functional derivative obtain: Minimizing the free energy and integrating by parts yields: Which is defined as the chemical potential:

23
Mathematical Background Combining the chemical potential with the definition of the diffusive current we obtain: When replaced in the continuity equation with use of the following partial derivatives: Yields the Cahn-Hilliard Equation:

24
Cahn-Hilliard Equation The Cahn-Hilliard equation, when integrated with respect to time, determines the time evolution of the composition of the system. Euler’s method will be used to solve this equation numerically giving us the composition as a function of position and time.  (r, t)

25
Length Scaling Of key interest is the size or length of the elongated domains as a function of time. The length is defined as: = average area of domains as a function of time A(t) is defined as the number of points making up the domain. = average perimeter of domains as a function of time P(t) is defined as the number of points on a domain that border a region made up of the opposite particle type. For a circular domain: A =  R 2 and P = 2  R, where R is the radius of the domain. The length parameter is then: L =R/2 or L(t) ~R(t)

26
Scaling law for droplets, non-rigorous approach Rewrite the continuity equation in terms of chemical potential: Using the Gibbs-Thompson relation for liquid droplets, the chemical potential is defined as: with  the surface tension, d is the dimensions of the system, and  is the change in order parameter at the droplet’s surface)

27
Scaling law for droplets, non-rigorous approach The continuity equation becomes: Which yields: After integrating the equation, dimensional analysis gives us a scaling law: R ~ t n, with n=1/3. Note that the critical exponent in the growth law is independent of the system’s dimensions.

28
Modeling a binary mixture undergoing spinodal decomposition (2-D study) 2-D square lattice is created with each lattice site receiving a random real number in the range (-1, 1). This is the initial value of  at each site,  (r, 0). Random number generator is weighted so that the overall particle concentrations are A = 52.5% and B = 47.5%. Goal is to find  (r, t)

29
Finding  (r, t)  (r, t) is found in the following manner: At each time step the Cahn-Hilliard equation is determined using the process described in the previous section. After being solved numerically, the Cahn-Hilliard equation produces a new value of  (r, t). These steps are repeated a set number of times for each time increment.

30
Qualitative description of particle dynamics Particle motion at each lattice site is determined by the interaction energy from 4 nearest neighbors. If adjacent sites are rich with the same particle type, energy is attractive. If adjacent sites do not share predominant particle type, energy is positive and repulsive.

31
Constraint on the system Continuity equation acts a restraint. Each particle type has a GLOBALLY conserved quantity. This conservation law is not enforced at each lattice site. The total energy of the system must be lowered for particle motion to be allowed. Energy may increase locally between neighbors if the total energy of the system is reduced.

32
Goals of the system (what the particles want to do) To reduce the total energy, the sum of the nearest neighbor energies, to a minimum value. This would occur if all the A-rich domains are in one continuous region, while all the B-rich domains are in a second continuous region. Examples include A-rich B-rich

33
Methods of observing the system Visual observation Images of the system Structure Factor Similar to conducting a light scattering experiment. Correlation length Using the Correlation function

34
Observing the system The program records a snapshot of the system after every time step. Each of these snapshots is a surface plot mapping of  (r) =  a (r) -  b (r). The time of each snapshot is determined by: Time = (snapshot #)(# of updates per scan)(delta time) The example below is the 10 th snapshot, with 5000 updates per scan, and dt = 0.01 units. Time = 10*5000*0.01=500 units.

36
Structure Factor The program calculates the 2-D Fourier Transform of each surface plot of  (r) =  a (r) -  b (r). With the radial wave number being:

37
Structure Factor K, is obtained from its respective wave number using the relation K = 2  k/D, where D is the size of the system. In our system, D is equal to 256. K max corresponds to the inverse of the most common length scale of the contiguous regions Of key interest is the time dependence of the maximum wave number/vector. k max (t) ~ t -  or K max (t) ~ t - 

39
Correlation length The correlation length is defined as the spatial range over which fluctuations in one region are correlated to fluctuations in another region. As domains increase in size, the number of adjacent lattice sites correlating with each other also will increase. The scaling law for the correlation length is: L(t) ~ t a The correlation length is also defined as the first zero of the correlation function :

41
Results Plots were made of the structure factor for time ranging up to 900 units

42
Time dependence of maximum wave number The power law found from the measurement resulted in a critical exponent of  = 0.3217. This compares favorably to the accepted value of  = 1/3.

43
The correlation function is shown for a range of times extending to 900 units

44
Time dependence of correlation length The power law found from the measurement resulted in a critical exponent of  = 0.2695. This is lower than the critical exponent for the structure factor plots by ~16%.

45
One point of interest is that by observing the snapshots up to a time value of 900 units, the domains are still predominantly elongated. The simulation was extended to 4950 units to see if the large circular domains became present. Check to see if the scaling law for the correlation length versus time changed over a longer time scale.

47
The correlation function is shown for a range of times extending to 4950 units

48
Time dependence of correlation length With the longer range of times, the critical exponent , was found to be 0.3091. This was in better agreement with the value for the critical exponent found from the structure factor data.

49
Conclusions Data from structure factor plots produced critical exponent of  = 0.3217, which was the closest calculate value of a to the accepted value of  = 1/3 Initially the correlation length data yielded  = 0.2695 When the time was extended by a factor of 5 the correlation length yielded  = 0.3091. Elongated domains remained present