The thermal resistivity and its scaling function in quasi-two-dimensional (2D) He4 systems are studied by Monte Carlo and spin-dynamics simulations. We use the classical 3D XY model on L×L×H lattices with L H, applying open boundary conditions along the H direction and periodic boundary conditions along the L directions. A hybrid Monte Carlo algorithm is adopted to efficiently deal with the critical slowing down and to produce initial states for time integration. The fourth-order Suzuki-Trotter decomposition method of exponential operators is used to solve numerically the coupled equations of motion for each spin. The thermal conductivity is calculated by a dynamic current-current correlation function. Our results show that (i) the simulational data collapse onto a single curve for several values of H and temperature, thus supporting the concept of finite-size scaling theory and (ii) the calculated scaling function agrees well with the available experimental results for slabs using two free fitting parameters.
|Original language||English (US)|
|Journal||Physical Review B - Condensed Matter and Materials Physics|
|State||Published - May 22 2006|
ASJC Scopus subject areas
- Electronic, Optical and Magnetic Materials
- Condensed Matter Physics