Thermodynamic and electronic properties are obtained for a lattice-gas model fluid with self-consistent, partial, occupation of its sites; the self consistency consists in obtaining ionic configurations from grand-canonical Monte Carlo simulations based on fits to the exact, electronic, tight-binding energies of isothermal ensembles of those same ionic configurations. The energy of an ion is found to be a concave-up function of its local coordination. Liquid-vapor coexistence densities and the electrical conductivity, which shows a metal-nonmetal transition, have been obtained.
1.IntroductionThe statistical mechanics of simple insulating liquids is a well developed subject, with different approaches being used to obtain the thermodynamics and the correlation structure, from the additive pairwise interaction potential between atoms or molecules [1]. A quantummechanical study of the atomic or molecular structures provides the interatomic potential, but, in all other aspects, the interaction is decoupled from the classical statistical-physics problem of obtaining the positions and correlations of the atoms. One of the more striking deviations from this simple liquid behavior is provided by liquid metals, in which the conduction electrons are fully delocalized and the system has to be treated as a mixture of ions ( with classical statistics) and electrons (with Fermi-Dirac statistics). The study of dense liquid metals, near the triple-point temperature, has been based on a double perturbative expansion [2] around a reference simple fluid for the ions, and around the jellium model for the conduction electrons. The vapor, at coexistence with the liquid metal, has a qualitatively different electronic structure, with the valence electrons localized in neutral atoms or clusters. At low temperature, the vapor has extremely low density and is almost trivial. In the neighborhood of the critical point, and of the metal-nonmetal transition region, the interrelation between electronic delocalization and ionic structure becomes crucial, and the approaches valid at low temperature fail qualitatively.Recent experimental data on the critical region of the alkali fluids [3] provide strong motivation for a more extensive theoretical study of these systems [4,5]. Our main objective here is to set a minimal model, including the main relevant features of the coupling between the electronic and the ionic structure, to analyze the statisticalphysics problem posed by these systems. We have reported [6] preliminary calculations using this approach. In our model, we forsake the 'state of the art' description of liquid metals in condensed matter physics, and our model neglects aspects of the problem which are certainly relevant for some properties of these systems. However, we show that many qualitative effects of the coupling between electronic and ionic structures are already present in a very simple model. This approach allows the study of the critical region through Monte Carlo simulation with system of large size, compared with ...