Effect of Gravity and Confinement on Phase Equilibria: A Density Matrix Renormalization Approach

Abstract

The phase diagram of the 2D Ising model confined between two infinite walls and subject to opposing surface fields and to a bulk "gravitational" field is calculated by means of density matrix renormalization methods. In absence of gravity two phase coexistence is restricted to temperatures below the wetting temperature. We find that gravity restores the two phase coexistence up to the bulk critical temperature, in agreement with previous mean-field predictions. We calculate the exponents governing the finite size scaling in the temperature and in the gravitational field directions. The former is the exponent which describes the shift of the critical temperature in capillary condensation. The latter agrees, for large surface fields, with a scaling assumption of Van Leeuwen and Sengers. Magnetization profiles in the two phase and in the single phase region are calculated. The profiles in the single phase region, where an interface is present, agree well with magnetization profiles calculated from a simple solid-on-solid interface hamiltonian.