In most of the discussion presented so far in this book, the quantum character of atoms and electrons has been ignored. The Ising spin models have been an exception, but since the Ising Hamiltonian is diagonal (in the absence of a transverse magnetic field), all energy eigenvalues are known and the Monte Carlo sampling can be carried out just as in the case of classical statistical mechanics. Furthermore, the physical properties are in accord with the third law of thermodynamics for Ising-type Hamiltonians (e.g. entropy S and specific heat vanish for temperature T → 0, etc.) in contrast to the other truly classical models dealt with in previous chapters (e.g. classical Heisenberg spin models, classical fluids and solids, etc.) which have many unphysical low temperature properties.
In most of the discussion presented so far in this book, the quantum character of atoms and electrons has been ignored. The Ising spin models have been an exception, but since the Ising Hamiltonian is diagonal (in the absence of a transverse magnetic field), all energy eigenvalues are known and the Monte Carlo sampling can be carried out just as in the case of classical statistical mechanics. Furthermore, the physical properties are in accord with the third law of thermodynamics for Ising-type Hamiltonians (e.g. entropy S and specific heat vanish for temperature T → 0, etc.) in contrast to the other truly classical models dealt with in previous chapters (e.g. classical Heisenberg spin models, classical fluids and solids, etc.) which have many unphysical low temperature properties.
INTRODUCTION In most of the discussion presented so far in this book, the quantum character of atoms and electrons has been ignored. The Ising spin models have been an exception, but since the Ising Hamiltonian is diagonal (in the absence of a transverse magnetic field!), all energy eigenvalues are known and the Monte Carlo sampling can be carried out just as in the case of classical statistical mechanics. Furthermore, the physical properties are in accord with the third law of thermodynamics for Ising-type Hamiltonians (e.g. entropy S and specific heat vanish for temperature T → 0, etc.) in contrast to the other truly classical models dealt with in previous chapters (e.g. classical Heisenberg spin models, classical fluids and solids, etc.) which have many unphysical low temperature properties. A case in point is a classical solid for which the specific heat follows the Dulong–Petit law, C = 3 Nk B , as T → 0, and the entropy has unphysical behavior since S → – ∞. Also, thermal expansion coefficients tend to non-vanishing constants for T → 0 while the third law implies that they must be zero. While the position and momentum of a particle can be specified precisely in classical mechanics, and hence the groundstate of a solid is a perfectly rigid crystal lattice (motionless particles localized at the lattice points), in reality the Heisenberg uncertainty principle forbids such a perfect rigid crystal, even at T → 0; due to zero point motions which ‘smear out’ the particles over some region around these lattice points.
INTRODUCTION In most of the discussion presented so far in this book, the quantum character of atoms and electrons has been ignored. The Ising spin models have been an exception, but since the Ising Hamiltonian is diagonal (in the absence of a transverse magnetic field!), all energy eigenvalues are known and the Monte Carlo sampling can be carried out just as in the case of classical statistical mechanics. Furthermore, the physical properties are in accord with the third law of thermodynamics for Ising-type Hamiltonians (e.g. entropy S and specific heat vanish for temperature T → 0, etc.) in contrast to the other truly classical models dealt with in previous chapters (e.g. classical Heisenberg spin models, classical fluids and solids, etc.) which have many unphysical low temperature properties. A case in point is a classical solid for which the specific heat follows the Dulong–Petit law, C = 3 Nk B , as T → 0, and the entropy has unphysical behavior since S → −∞. Also, thermal expansion coefficients tend to non-vanishing constants for T → 0 while the third law implies that they must be zero. While the position and momentum of a particle can be specified precisely in classical mechanics, and hence the groundstate of a solid is a perfectly rigid crystal lattice (motionless particles localized at the lattice points), in reality the Heisenberg uncertainty principle forbids such a perfect rigid crystal, even at T → 0; due to zero point motions which ‘smear out’ the particles over some region around these lattice points.
Introduction In most of the discussion presented so far in this book, the quantum character of atoms and electrons has been ignored. The Ising spin models have been an exception, but since the Ising Hamiltonian is diagonal (in the absence of a transverse magnetic field), all energy eigenvalues are known and the Monte Carlo sampling can be carried out just as in the case of classical statistical mechanics. Furthermore, the physical properties are in accord with the third law of thermodynamics for Ising-type Hamiltonians (e.g. entropy S and specific heat vanish for temperature T → 0, etc.) in contrast to the other truly classical models dealt with in previous chapters (e.g. classical Heisenberg spin models, classical fluids and solids, etc.) which have many unphysical low temperature properties. A case in point is a classical solid for which the specific heat follows the Dulong–Petit law, C = 3 Nk B , as T → 0, and the entropy has unphysical behavior since S → –∞. Also, thermal expansion coefficients tend to non-vanishing constants for T → 0 while the third law implies that they must be zero. While the position and momentum of a particle can be specified precisely in classical mechanics, and hence the groundstate of a solid is a perfectly rigid crystal lattice (motionless particles localized at the lattice points), in reality the Heisenberg uncertainty principle forbids such a perfect rigid crystal, even at T → 0, due to zero point motions which ‘smear out’ the particles over some region around these lattice points. This delocalization of quantum-mechanical particles increases as the atomic mass is reduced; therefore, these quantum effects are most pronounced for light atoms like hydrogen in metals, or liquid helium. Spectacular phenomena like superfluidity are a consequence of the quantum nature of the particles and have no classical counterpart at all. Even for heavier atoms, which do not show superfluidity because the fluid–solid transition intervenes before a transition from normal fluid to superfluid could occur, there are genuine effects of quantum nature. Examples include the isotope effects (remember that in classical statistical mechanics the kinetic energy part of the Boltzmann factor cancels out from all averages, and thus in thermal equilibrium no property depends explicitly on the mass of the particles).
As we have seen on several occasions, a great many of the physical phenomena of interest concern collective behaviour or strongly correlated particles. While it is sometimes possible to obtain useful insight using approximation methods, such as perturbation theory, reliable controlled approximations do not exist. In such situations, numerical simulations have become an indispensable tool. We shall not attempt in this chapter to discuss computer programming nor shall we discuss many of the powerful methods, both algorithmic and in data analysis, that have been developed since numerical simulations have become an important research tool. Instead, we shall assume that the reader already has enough familiarity with computer programming to apply the notions which we discuss here. Our goal is to explain enough of the basics of the Monte Carlo method (classical and quantum) to allow the reader to apply it readily to interesting equilibrium statistical physics problems of the type discussed in this book. The problems in this chapter have been tailored with this in mind. The programs needed are relatively simple and the physics illustrates many of the phenomena already seen: phase transitions with and without spontaneous symmetry breaking, critical exponents, diverging susceptibilities, scaling, vortices, superfluids etc. Consequently, some of the problems are somewhat long and may be better considered as mini-projects. Markov chains, convergence and detailed balance Even relatively innocent looking models are impossible to solve exactly. The exact solution of the two-dimensional Ising model is known but not in the presence of an external magnetic field. In three dimensions, even in the absence of a magnetic field, the exact solution is not known. In fact, relatively few models can be solved exactly; see Reference [14] for an excellent review of the subject.
We investigate the thermodynamic properties of the frustrated bilayer quantum Heisenberg antiferromagnet at low temperatures in the vicinity of the saturation magnetic field. The low-energy degrees of freedom of the spin model are mapped onto a hard-square gas on a square lattice. We use exact diagonalization data for finite spin systems to check the validity of such a description. Using a classical Monte Carlo method we give a quantitative description of the thermodynamics of the spin model at low temperatures around the saturation field. The main peculiarity of the considered two-dimensional Heisenberg antiferromagnet is related to a phase transition of the hard-square model on the square lattice, which belongs to the two-dimensional Ising model universality class. It manifests itself in a logarithmic (low-) temperature singularity of the specific heat of the spin system observed for magnetic fields just below the saturation field.
Thermodynamic properties of S=2 antiferromagnetic Heisenberg chains are studied not only under the periodic boundary condition but also under the open one employing a quantum Monte Carlo method. Temperature and size dependences of the energy, the specific heat, and the magnetic susceptibility are calculated and edge effects on them are investigated in detail. The specific heat shows a well-pronounced Schottky anomaly but the maximum is located at a temperature much larger than the Haldane gap of the system. As temperature goes to zero, the magnetic susceptibility vanishes for the periodic chains, while it diverges for the open chains. The edge contribution in the open-chain susceptibility is attributed to the two S=1 effective spins localized at the chain ends at low temperatures, while to an S=2 free spin at high temperatures. This is a visualization of a quantum-classical crossover and an evidence that the present model is a Haldane system. High-temperature behaviors of the thermodynamic quantities are also discussed with the help of a series-expansion method. \textcopyright{} 1996 The American Physical Society.