Scaling of Phase Diagram and Critical Point Parameters in Liquid-Vapour Phase Transition of Metallic Fluids

Abstract

The first objective of this paper is to investigate the scaling behavior of liquid-vapor phase transition in FCC and BCCmetals starting from the zero-temperature four-parameter formula for cohesive energy. The effective potentials between the atoms in the solid are determined while using lattice inversion techniques as a function of scaling variables in the four-parameter formula. These potentials are split into repulsive and attractive parts, as per the Weeks–Chandler–Anderson prescription, and used in the coupling-parameter expansion for solving the Ornstein–Zernike equation that was supplemented with an accurate closure. Thermodynamic quantities obtained via the correlation functions are used in order to obtain critical point parameters and liquid-vapor phase diagrams. Their dependence on the scaling variables in the cohesive energy formula are also determined. An equally important second objective of the paper is to revisit coupling parameter expansion for solving the Ornstein–Zernike equation. The Newton–Armijo non-linear solver and Krylov-space based linear solvers are employed in this regard. These methods generate a robust algorithm that can be used to span the entire fluid region, except very low temperatures. The accuracy of the method is established by comparing the phase diagrams with those that were obtained via computer simulation. The avoidance of the ’no-solution-region’ of the Ornstein-Zernike equation in coupling-parameter expansion is also discussed. Details of the method and complete algorithm provided here would make this technique more accessible to researchers investigating the thermodynamic properties of one component fluids.