Contributor info: Dr Roman Rausch

Robert Peters, Roman Rausch

SciPost Phys. 14, 166 (2023) · published 21 June 2023 | Toggle abstract · pdf

Using variational matrix product states, we analyze the finite temperature behavior of a half-filled periodic Anderson model in one dimension, a prototypical model of a Kondo insulator. We present an extensive analysis of single-particle Green's functions, two-particle Green's functions, and transport functions creating a broad picture of the low-temperature properties. We confirm the existence of energetically low-lying spin excitations in this model and study their energy-momentum dispersion and temperature dependence. We demonstrate that charge-charge correlations at the Fermi energy exhibit a different temperature dependence than spin-spin correlations. While energetically low-lying spin excitations emerge approximately at the Kondo temperature, which exponentially depends on the interaction strength, charge correlations vanish already at high temperatures. Furthermore, we analyze the charge and thermal conductivity at finite temperatures by calculating the time-dependent current-current correlation functions. While both charge and thermal conductivity can be fitted for all interaction strengths by gapped systems with a renormalized band gap, the gap in the system describing the thermal conductivity is generally smaller than the system describing the charge conductivity. Thus, two-particle correlations affect the charge and heat conductivities in a different way resulting in a temperature region where the charge conductivity of this one-dimensional Kondo insulator is already decreasing while the heat conductivity is still increasing.

Roman Rausch, Matthias Peschke, Cassian Plorin, Jürgen Schnack, Christoph Karrasch

SciPost Phys. 14, 052 (2023) · published 28 March 2023 | Toggle abstract · pdf

The ferrimagnetic phase of the sawtooth chain with mixed ferromagnetic nearest-neighbour interactions $J$ and antiferromagnetic next-nearest-neighbour interactions $J'$ (within the isotropic Heisenberg model) was previously characterized as a phase with commensurate order. In this paper, we demonstrate that the system in fact exhibits an incommensurate quantum spin spiral. Even though the ground state is translationally invariant in terms of the local spin expectations $\left<{\vec{S}_i}\right>$, the spiral can be detected via the connected spin-spin correlations $\left<{\vec{S}_i\cdot\vec{S}_j}\right>-\left<{\vec{S}_i}\right>\cdot\left<{\vec{S}_j}\right>$ between the apical spins. It has a long wavelength that grows with $J'$ and that soon exceeds finite-system sizes typically employed in numerical simulations. A faithful treatment thus requires the use of state-of-the-art simulations for large, periodic systems. In this work, we are able to accurately treat up to $L=400$ sites (200 unit cells) with periodic boundary conditions using the density-matrix renormaliztion group (DMRG). Exploiting the SU(2) symmetry allows us to directly compute the lowest-energy state for a given total spin. Our results are corroborated by variational uniform matrix product state (VUMPS) calculations, which work directly in the thermodynamic limit at the cost of a lower accuracy.

Roman Rausch, Matthias Peschke, Cassian Plorin, Christoph Karrasch

SciPost Phys. 12, 143 (2022) · published 2 May 2022 | Toggle abstract · pdf

We compute ground-state properties of the isotropic, antiferromagnetic Heisenberg model on the sodalite cage geometry. This is a 60-spin spherical molecule with 24 vertex-sharing tetrahedra which can be regarded as a molecular analogue of a capped kagome lattice and which has been synthesized with high-spin rare-earth atoms. Here, we focus on the $S=1/2$ case where quantum effects are strongest. We employ the SU(2)-symmetric density-matrix renormalization group (DMRG). We find a threefold degenerate ground state that breaks the spatial symmetry and that splits up the molecule into three large parts which are almost decoupled from each other. This stands in sharp contrast to the behaviour of most known spherical molecules. On a methodological level, the disconnection leads to ``glassy dynamics'' within the DMRG that cannot be targeted via standard techniques. In the presence of finite magnetic fields, we find broad magnetization plateaus at 4/5, 3/5, and 1/5 of the saturation, which one can understand in terms of localized magnons, singlets, and doublets which are again nearly decoupled from each other. At the saturation field, the zero-point entropy is $S=\ln(182)\approx 5.2$ in units of the Boltzmann constant.

Roman Rausch, Cassian Plorin, Matthias Peschke

SciPost Phys. 10, 087 (2021) · published 23 April 2021 | Toggle abstract · pdf

We solve the quantum-mechanical antiferromagnetic Heisenberg model with spins positioned on vertices of the truncated icosahedron using the density-matrix renormalization group (DMRG). This describes magnetic properties of the undoped C$_{60}$ fullerene at half filling in the limit of strong on-site interaction $U$. We calculate the ground state and correlation functions for all possible distances, the lowest singlet and triplet excited states, as well as thermodynamic properties, namely the specific heat and spin susceptibility. We find that unlike smaller C$_{20}$ or C$_{32}$ that are solvable by exact diagonalization, the lowest excited state is a triplet rather than a singlet, indicating a reduced frustration due to the presence of many hexagon faces and the separation of the pentagonal faces, similar to what is found for the truncated tetrahedron. This implies that frustration may be tuneable within the fullerenes by changing their size. The spin-spin correlations are much stronger along the hexagon bonds and exponentially decrease with distance, so that the molecule is large enough not to be correlated across its whole extent. The specific heat shows a high-temperature peak and a low-temperature shoulder reminiscent of the kagomé lattice, while the spin susceptibility shows a single broad peak and is very close to the one of C$_{20}$.