The orientational ordering transition is investigated in the quantum generalization of the anisotropic-planar-rotor model in the low temperature regime.The phase diagram of the model is first analyzed within the mean-field approximation. This predicts at T = 0 a phase transition from the ordered to the disordered state when the strength of quantum fluctuations, characterized by the rotational constant Θ, exceeds a critical value Θ M F c . As a function of temperature, mean-field theory predicts a range of values of Θ where the system develops long-range order upon cooling, but enters again into a disordered state at sufficiently low temperatures (reentrance). The model is further studied by means of path integral Monte Carlo simulations in combination with finite-size scaling techniques, concentrating on the region of parameter space where reentrance is predicted to occur. The phase diagram determined from the simulations does not seem to exhibit reentrant behavior; at intermediate temperatures a pronounced increase of short-range order is observed rather 1 than a genuine long-range order. 05.70.Fh, 64.60.Cn, 68.35.Rh Typeset using REVT E X 2
I. MOTIVATIONPhysisorbates are experimental realizations of quasi two-dimensional systems that display an extremely rich phase behavior due to the competition between intermolecular and molecule-surface interactions, as documented, e.g., in Refs. [1][2][3][4][5]. Correspondingly there is a wealth of phase transitions between the various ordered phases as a function of temperature and coverage. Since many of these transitions occur at fairly low temperatures, quantum effects might play an important or even crucial role [3,4], as recently demonstrated for the ordering of hydrogen isotopes on graphite [6]. Molecular systems are particularly interesting as they possess orientational degrees of freedom which can order in addition to the positions [5].In case of linear molecules, the anisotropic-planar-rotor (APR) model [7,8] was devised to describe the herringbone (quadrupolar orientational two-sublattice) ordering transition [5], e.g., in commensurate N 2 monolayers on graphite. The classical two-dimensional APR model consists of planar rotators pinned with their center of rotation on a triangular lattice and interacting via nearest-neighbor quadrupolar interactions only; a three-dimensional version has also been proposed and investigated in various approximations [9].Over the years the APR Hamiltonian acquired the status of a statistical mechanical model on its own right, see Refs. [10,5]. Many of these activities arose because the order of the APR phase transition turned out to be extremely challenging to determine [11], finally favoring a first order phase transition that is "weak" and fluctuation driven [12], see the detailed review in Ref. [5]. The plain APR model was generalized to include vacancies or impurities [13][14][15], as well as multipole interactions other than those of quadrupolar symmetry [14,16]. Simplified Hamiltonians were obtained after discreti...