Search

The surface and bulk properties of the two-dimensional Q > 4 state Potts model in the vicinity of the first order bulk transition point have been studied by exact calculations and by density matrix renormalization group techniques. For the surface transition the complete analytical solution of the problem is presented in the $Q \to \infty$ limit, including the critical and tricritical exponents, magnetization profiles and scaling functions. According to the accurate numerical results the universality class of the surface transition is independent of the value of Q > 4. For the bulk transition we have numerically calculated the latent heat and the magnetization discontinuity and we have shown that the correlation lengths in the ordered and in the disordered phases are identical at the transition point., Comment: 11 pages, RevTeX, 6 PostScript figures included. Manuscript substantially extended, details on the analytical and numerical calculations added. To appear in Phys. Rev. B

We consider an Ising model confined in an $L \times \infty$ geometry with identical surface fields at the boundaries. According to the Kelvin equation the bulk coexistence field scales as 1/L for large L; thermodynamics and scaling arguments predict higher order corrections of the type $1/L^2$ and $1/L^{5/3}$ at partial and complete wetting respectively. Our numerical results, obtained by density-matrix renormalization techniques for systems of widths up to L = 144, are in agreement with a $1/L^2$ correction in the partial wet regime. However at complete wetting we find a large range of surface fields and temperatures with a correction to scaling of type $1/L^{4/3}$. We show that this term is generated by a {\it thin} wetting layer whose free energy is dominated by the contacts with the wall. For $L$ sufficiently large we expect a crossover to a $1/L^{5/3}$ correction as predicted by the theory for a {\it thick} wetting layer. This crossover is indeed found in a solid-on-solid model which provides a simplified description of the wetting layer and allows the study of much larger systems than those available in a density-matrix renormalization calculation., Comment: 9 pages, RevTeX, 7 PostScript figures included. Some misprints corrected. Final version as published

We use the density matrix renormalization group (DMRG) to perform a detailed study of the critical properties of the two dimensional Q state Potts model, including the magnetization and energy-density profiles, bulk and surface critical exponents and the Casimir amplitudes. We apply symmetry breaking boundary conditions to a $L \times \infty$ strip and diagonalize the corresponding transfer matrix for a series of moderately large systems ($L \le 64$) by the DMRG method. The numerically very accurate finite lattice results are then extrapolated by efficient sequence extrapolation techniques. The critical density profiles and the Casimir amplitudes are found to follow precisely the conformal predictions for Q=2 and 3. Similarly, the bulk and surface critical exponents of the models are in very good agreement with the conformal and exact values: their accuracy has reached or even exceeded the accuracy of other available numerical methods. For the Q=4 model both the profiles and the critical exponents show strong logarithmic corrections, which are also studied., Comment: 12 pages, RevTeX, 14 PostScript figures included

The properties of a fluid, or Ising magnet, confined in a $L \times \infty$ geometry with opposing surface fields at the walls are studied by density matrix renormalization techniques. In particular we focus on the effect of gravity on the system, which is modeled by a bulk field whose strength varies linearly with the distance from the walls. It is well known that in the absence of gravity phase coexistence is restricted to temperatures below the wetting temperature. We find that gravity restores phase coexistence up to the bulk critical temperature, in agreement with previous mean field results. A detailed study of the scaling to the critical point, as $L \to \infty$, is performed. The temperature shift scales as $1/L^{y_T}$, while the gravitational constant scales as $1/L^{1+y_H}$, with $y_T$ and $y_H$ the bulk thermal and magnetic exponents respectively. For weak surface fields and $L$ not too large, we also observe a regime where the gravitational constant scales as $1/L^{1+y_H - \Delta_1 y_T}$ ($\Delta_1$ is the surface gap exponent) with a crossover, for sufficiently large $L$, to a scaling of type $1/L^{1+y_H}$., Comment: 9 pages, RevTeX, 11 PostScript figures included. Minor corrections. Final version as published

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., Comment: 4 pages, RevTeX and 4 PostScript figures included. Final version as published. To appear in Phys. Rev. Lett

Networks of steps, seen in STM observations of vicinal surfaces on Au and Pt (110), are analyzed. A simple model is introduced for the calculation of the free energy of the networks as function of the slope parameters, valid at low step densities. It predicts that the networks are unstable, or at least metastable, against faceting and gives an equilibrium crystal shape with sharp edges either between the (110) facet and rounded regions or between two rounded regions. Experimental observations of the equilibrium shapes of Au or Pt crystals at sufficiently low temperatures, i.e. below the deconstruction temperature of the (110) facet, could check the validity of these predictions., Comment: 14 pages, RevTex, uses epsf.sty, the file rotate.sty is included. To appear in the Proceedings of the Traverse City workshop on The Dynamics of Crystal Surfaces (August 1996); this work is an extension of Phys. Rev. Lett. {\bf 76}, 4191 (1996)

The phase diagram of the staggered six vertex, or body centered solid on solid model, is investigated by transfer matrix and finite size scaling techniques. The phase diagram contains a critical region, bounded by a Kosterlitz-Thouless line, and a second order line describing a deconstruction transition. In part of the phase diagram the deconstruction line and the Kosterlitz-Thouless line approach each other without merging, while the deconstruction changes its critical behaviour from Ising-like to a different universality class. Our model has the same type of symmetries as some other two-dimensional models, such as the fully frustrated XY model, and may be important for understanding their phase behaviour. The thermal behaviour for weak staggering is intricate. It may be relevant for the description of surfaces of ionic crystals of CsCl structure., Comment: 13 pages, RevTex, 1 Postscript file with all figures, to be published in Phys. Rev. B