Monte carlo study of the phase diagrams of binary alloys with face centered cubic lattice structure

Universal critical behavior of curvature-dependent interfacial tension

From the analysis of Monte Carlo simulations of a binary Lennard-Jones mixture in the coexistence region, we provide evidence that the curvature dependence of the interfacial tension can be described by a simple theoretical function σ(R)ξ(2)=C(1)/[1+C(2)(ξ/R)(2)], where ξ is the correlation length and R is the droplet radius. The universal constants C(1) and C(2) are estimated. In the model, a Tolman length is strictly absent, but, since its critical behavior is believed to be much weaker than ξ, we argue that it only provides a correction to scaling and does not affect the leading critical behavior, which should be described by the above function for any system in the Ising universality class. The large value of C(2)≃32 implies that conventional nucleation theory becomes inaccurate even for a significantly large droplet radius.

Thermodynamic properties of a symmetrical binary mixture in the coexistence region

A three-dimensional symmetric binary fluid is studied, as a function of temperature, in the two-phase (liquid-liquid) coexistence region via Monte Carlo simulations. Particular focus has been in the understanding of curvature-dependent interfacial tension, which is observed to vary as σ(R) = σ(∞)/[1+2(ℓ/R)(2)], implying that a Tolman length is zero in the limit R → ∞. The length ℓ is found to have a critical divergence the same as the correlation length, but its amplitude is significantly larger (ℓ ~/= 4ξ). Our findings hence imply that the barrier against homogeneous nucleation is significantly reduced (in comparison with the classical nucleation theory) in the critical region. We also report results for the critical behavior of the flat interfacial tension σ(∞) and the concentration susceptibility, as well as the amplitude ratios involving these thermodynamic quantities. Noting that the interatomic potential in our model is described by the Lennard-Jones form that decays faster that 1/r(3), all of our results for critical phenomena are expectedly consistent with the Ising universality class of three spatial dimensions.

On the Feasibility of AB Initio Calculations of Ordering Alloy Phase Diagrams

Many stable or metastable intermetallic phases useful to the alloy designer have crystal structures which are ordered superstructures of a parent disordered phase. A highly reliable statistical mechanical method (CVM) has now been developed for calculating such superstructure phase equilibria derived from say, the fcc parent lattice. To obtain phase diagrams, one needs certain physical parameters, such as effective pair interaction ratios. It is possible, in principle, to extract these parameters from band structure calculations in the coherent potential approximation (CPA), particularly from recently developed cluster-CPA techniques. If sufficient accuracy can be achieved, truly first-principles phase diagram calculations may soon become feasible.

Static and dynamic critical behavior of a symmetrical binary fluid: A computer simulation

A symmetrical binary, A+B Lennard-Jones mixture is studied by a combination of semi-grand-canonical Monte Carlo (SGMC) and molecular dynamics (MD) methods near a liquid-liquid critical temperature T(c). Choosing equal chemical potentials for the two species, the SGMC switches identities (A-->B-->A) to generate well-equilibrated configurations of the system on the coexistence curve for T<T(c) and at the critical concentration, x(c)=12, for T>T(c). A finite-size scaling analysis of the concentration susceptibility above T(c) and of the order parameter below T(c) is performed, varying the number of particles from N=400 to 12 800. The data are fully compatible with the expected critical exponents of the three-dimensional Ising universality class. The equilibrium configurations from the SGMC runs are used as initial states for microcanonical MD runs, from which transport coefficients are extracted. Self-diffusion coefficients are obtained from the Einstein relation, while the interdiffusion coefficient and the shear viscosity are estimated from Green-Kubo expressions. As expected, the self-diffusion constant does not display a detectable critical anomaly. With appropriate finite-size scaling analysis, we show that the simulation data for the shear viscosity and the mutual diffusion constant are quite consistent both with the theoretically predicted behavior, including the critical exponents and amplitudes, and with the most accurate experimental evidence.