Correlations in the low-temperature phase of the two-dimensional XY model

Monte Carlo simulations of the two-dimensional XY model are performed in a square geometry with fixed boundary conditions. Using a conformal mapping, it is very easy to deduce the exponent eta(sigma)(T) of the order parameter correlation function at any temperature in the critical phase of the model. The temperature behaviour of eta(sigma)(T) is obtained numerically with a good accuracy up to the Kosterlitz-Thouless transition temperature. At very low temperatures, a good agreement is found with Berezinskii's harmonic approximation. Surprisingly, we show some evidence that there are no logarithmic corrections to the behaviour of the order parameter density profile (with symmetry-breaking surface fields) at the Kosterlitz-Thouless transition temperature.