Source author record

S. Wessel

S. Wessel appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2021arXiv

Quantum Monte Carlo simulations in the trimer basis: first-order transitions and thermal critical points in frustrated trilayer magnets

The phase diagrams of highly frustrated quantum spin systems can exhibit first-order quantum phase transitions and thermal critical points even in the absence of any long-ranged magnetic order. However, all unbiased numerical techniques for investigating frustrated quantum magnets face significant challenges, and for generic quantum Monte Carlo methods the challenge is the sign problem. Here we report on a general quantum Monte Carlo approach with a loop-update scheme that operates in any basis, and we show that, with an appropriate choice of basis, it allows us to study a frustrated model of coupled spin-1/2 trimers: simulations of the trilayer Heisenberg antiferromagnet in the spin-trimer basis are sign-problem-free when the intertrimer couplings are fully frustrated. This model features a first-order quantum phase transition, from which a line of first-order transitions emerges at finite temperatures and terminates in a thermal critical point. The trimer unit cell hosts an internal degree of freedom that can be controlled to induce an extensive entropy jump at the quantum transition, which alters the shape of the first-order line. We explore the consequences for the thermal properties in the vicinity of the critical point, which include profound changes in the lines of maxima defined by the specific heat. Our findings reveal trimer quantum magnets as fundamental systems capturing in full the complex thermal physics of the strongly frustrated regime.

preprint2016arXiv

Thermodynamic properties of highly frustrated quantum spin ladders: influence of many-particle bound states

Quantum antiferromagnets have proven to be some of the cleanest realizations available for theoretical, numerical, and experimental studies of quantum fluctuation effects. At finite temperatures, however, the additional effects of thermal fluctuations in the restricted phase space of a low-dimensional system have received much less attention, particularly the situation in frustrated quantum magnets, where the excitations may be complex collective (bound or even fractionalized) modes. We investigate this problem by studying the thermodynamic properties of the frustrated two-leg S=1/2 spin ladder, with particular emphasis on the fully frustrated case. We present numerical results for the magnetic specific heat and susceptibility, obtained from exact diagonalization and quantum Monte Carlo studies, which we show can be rendered free of the sign problem even in a strongly frustrated system and which allow us to reach unprecedented sizes of L=200 ladder rungs. We find that frustration effects cause an unconventional evolution of the thermodynamic response across the full parameter regime of the model. However, close to the first-order transition they cause a highly anomalous reduction in temperature scales with no concomitant changes in the gap; the specific heat shows a very narrow peak at very low energies and the susceptibility rises abruptly at extremely low temperatures. Unusually, the two quantities have different gaps over an extended region of the parameter space. We demonstrate that these results reflect the presence of large numbers of multi-particle bound-state excitations, whose energies fall below the one-triplon gap in the transition region.

preprint2015arXiv

Dynamical structure factors and excitation modes of the bilayer Heisenberg model

Using quantum Monte Carlo simulations along with higher-order spin-wave theory, bond-operator and strong-coupling expansions, we analyse the dynamical spin structure factor of the spin-half Heisenberg model on the square-lattice bilayer. We identify distinct contributions from the low-energy Goldstone modes in the magnetically ordered phase and the gapped triplon modes in the quantum disordered phase. In the antisymmetric (with respect to layer inversion) channel, the dynamical spin structure factor exhibits a continuous evolution of spectral features across the quantum phase transition, connecting the two types of modes. Instead, in the symmetric channel we find a depletion of the spectral weight when moving from the ordered to the disordered phase. While the dynamical spin structure factor does not exhibit a well-defined distinct contribution from the amplitude (or Higgs) mode in the ordered phase, we identify an only marginally-damped amplitude mode in the dynamical singlet structure factor, obtained from interlayer bond correlations, in the vicinity of the quantum critical point. These findings provide quantitative information in direct relation to possible neutron or light scattering experiments in a fundamental two-dimensional quantum-critical spin system.

preprint2015arXiv

Finite-size effects in Luther-Emery phases of Holstein and Hubbard models

The one-dimensional Holstein model and its generalizations have been studied extensively to understand the effects of electron-phonon interaction. The half-filled case is of particular interest, as it describes a transition from a metallic phase with a spin gap due to attractive backscattering to a Peierls insulator with charge-density-wave (CDW) order. Our quantum Monte Carlo results support the existence of a metallic phase with dominant power-law charge correlations, as described by the Luther-Emery fixed point. We demonstrate that for Holstein and also for purely fermionic models the spin gap significantly complicates finite-size numerical studies, and explains inconsistent previous results for Luttinger parameters and phase boundaries. On the other hand, no such complications arise in spinless models. The correct low-energy theory of the spinful Holstein model is argued to be that of singlet bipolarons with a repulsive, mutual interaction. This picture naturally explains the existence of a metallic phase, but also implies that gapless Luttinger liquid theory is not applicable.

preprint2012arXiv

Quantum phase transitions in the Kane-Mele-Hubbard model

We study the two-dimensional Kane-Mele-Hubbard model at half filling by means of quantum Monte Carlo simulations. We present a refined phase boundary for the quantum spin liquid. The topological insulator at finite Hubbard interaction strength is adiabatically connected to the groundstate of the Kane-Mele model. In the presence of spin-orbit coupling, magnetic order at large Hubbard U is restricted to the transverse direction. The transition from the topological band insulator to the antiferromagnetic Mott insulator is in the universality class of the three-dimensional XY model. The numerical data suggest that the spin liquid to topological insulator and spin liquid to Mott insulator transitions are both continuous.

preprint2011arXiv

The ALPS project release 2.0: Open source software for strongly correlated systems

We present release 2.0 of the ALPS (Algorithms and Libraries for Physics Simulations) project, an open source software project to develop libraries and application programs for the simulation of strongly correlated quantum lattice models such as quantum magnets, lattice bosons, and strongly correlated fermion systems. The code development is centered on common XML and HDF5 data formats, libraries to simplify and speed up code development, common evaluation and plotting tools, and simulation programs. The programs enable non-experts to start carrying out serial or parallel numerical simulations by providing basic implementations of the important algorithms for quantum lattice models: classical and quantum Monte Carlo (QMC) using non-local updates, extended ensemble simulations, exact and full diagonalization (ED), the density matrix renormalization group (DMRG) both in a static version and a dynamic time-evolving block decimation (TEBD) code, and quantum Monte Carlo solvers for dynamical mean field theory (DMFT). The ALPS libraries provide a powerful framework for programers to develop their own applications, which, for instance, greatly simplify the steps of porting a serial code onto a parallel, distributed memory machine. Major changes in release 2.0 include the use of HDF5 for binary data, evaluation tools in Python, support for the Windows operating system, the use of CMake as build system and binary installation packages for Mac OS X and Windows, and integration with the VisTrails workflow provenance tool. The software is available from our web server at http://alps.comp-phys.org/.

preprint2011arXiv

Two-particle decay and quantum criticality in dimerized antiferromagnets

In certain Mott-insulating dimerized antiferromagnets, triplet excitations of the paramagnetic phase can decay into the two-particle continuum. When such a magnet undergoes a quantum phase transition into a magnetically ordered state, this coupling becomes part of the critical theory provided that the lattice ordering wavevector is zero. One microscopic example is the staggered-dimer antiferromagnet on the square lattice, for which deviations from O(3) universality have been reported in numerical studies. Using both symmetry arguments and microscopic calculations, we show that a non-trivial cubic term arises in the relevant order-parameter quantum field theory, and assess its consequences using a combination of analytical and numerical methods. We also present finite-temperature quantum Monte Carlo data for the staggered-dimer antiferromagnet which complement recently published results. The data can be consistently interpreted in terms of critical exponents identical to that of the standard O(3) universality class, but with anomalously large corrections to scaling. We argue that the two-particle decay of critical triplons, although irrelevant in two spatial dimensions, is responsible for the leading corrections to scaling due to its small scaling dimension.

preprint2010arXiv

Quantum spin-liquid emerging in two-dimensional correlated Dirac fermions

At sufficiently low temperatures, condensed-matter systems tend to develop order. An exception are quantum spin-liquids, where fluctuations prevent a transition to an ordered state down to the lowest temperatures. While such states are possibly realized in two-dimensional organic compounds, they have remained elusive in experimentally relevant microscopic two-dimensional models. Here, we show by means of large-scale quantum Monte Carlo simulations of correlated fermions on the honeycomb lattice, a structure realized in graphene, that a quantum spin-liquid emerges between the state described by massless Dirac fermions and an antiferromagnetically ordered Mott insulator. This unexpected quantum-disordered state is found to be a short-range resonating valence bond liquid, akin to the one proposed for high temperature superconductors. Therefore, the possibility of unconventional superconductivity through doping arises. We foresee its realization with ultra-cold atoms or with honeycomb lattices made with group IV elements.

preprint2009arXiv

Classical and quantum anisotropic Heisenberg antiferromagnets

We study classical and quantum Heisenberg antiferromagnets with exchange anisotropy of XXZ-type and crystal field single-ion terms of quadratic and cubic form in a field. The magnets display a variety of phases, including the spin-flop (or, in the quantum case, spin-liquid) and biconical (corresponding, in the quantum lattice gas description, to supersolid) phases. Applying ground-state considerations, Monte Carlo and density matrix renormalization group methods, the impact of quantum effects and lattice dimension is analysed. Interesting critical and multicritical behaviour may occur at quantum and thermal phase transitions.

preprint2009arXiv

Magnetocaloric effect in quantum spin-s chains

We compute the entropy of antiferromagnetic quantum spin-s chains in an external magnetic field using exact diagonalization and Quantum Monte Carlo simulations. The magnetocaloric effect, i.e., temperature variations during adiabatic field changes, can be derived from the isentropes. First, we focus on the example of the spin-s=1 chain and show that one can cool by closing the Haldane gap with a magnetic field. We then move to quantum spin-s chains and demonstrate linear scaling with $s$ close to the saturation field. In passing, we propose a new method to compute many low-lying excited states using the Lanczos recursion.

preprint2008arXiv

Three-body Interactions in one Dimension

We determine the phase-diagram of a one-dimensional system of hard-core lattice bosons interacting via repulsive three-body interactions by analytic methods and extensive quantum Monte-Carlo simulations. Such three-body interactions can be derived from a microscopic theory for polar molecules trapped in an optical lattice. Depending on the strength of the interactions and the particle density, we find superfluid and solid phases, the latter appearing at an unconventional filling of the lattice and displaying a coexistence of charge-density-wave and bond orders.

preprint2007arXiv

Classical and quantum two-dimensional anisotropic Heisenberg antiferromagnets

The classical and the quantum, spin $S=1/2, versions of the uniaxially anisotropic Heisenberg antiferromagnet on a square lattice in a field parallel to the easy axis are studied using Monte Carlo techniques. For the classical version, attention is drawn to biconical structures and fluctuations at low temperatures in the transition region between the antiferromagnetic and spin-flop phases. For the quantum version, the previously proposed scenario of a first-order transition between the antiferromagnetic and spin-flop phases with a critical endpoint and a tricritical point is scrutinized.

preprint2006arXiv

Quenching of the Haldane gap in LiVSi2O6 and related compounds

We report results of susceptibility chi and 7Li NMR measurements on LiVSi2O6. The temperature dependence of the magnetic susceptibility chi(T) exhibits a broad maximum, typical for low-dimensional magnetic systems. Quantitatively it is in agreement with the expectation for an S=1 spin chain, represented by the structural arrangement of V ions. The NMR results indicate antiferromagnetic ordering below T_N=24 K. The intra- and interchain coupling J and J_p for LiVSi2O6, and also for its sister compounds LiVGe2O6, NaVSi2O6 and NaVGe2O6, are obtained via a modified random phase approximation which takes into account results of quantum Monte Carlo calculations. While J_p is almost constant across the series, J varies by a factor of 5, decreasing with increasing lattice constant along the chain direction. The comparison between experimental and theoretical susceptibility data suggests the presence of an easy-axis magnetic anisotropy, which explains the formation of an energy gap in the magnetic excitation spectrum below T_N, indicated by the variation of the NMR spin-lattice relaxation rate at T << T_N.

preprint2004arXiv

Bosons in optical lattices - from the Mott transition to the Tonks-Girardeau gas

We present results from quantum Monte Carlo simulations of trapped bosons in optical lattices, focusing on the crossover from a gas of softcore bosons to a Tonks-Girardeau gas in a one-dimensional optical lattice. We find that depending on the quantity being measured, the behavior found in the Tonks-Girardeau regime is observed already at relatively small values of the interaction strength. A finite critical value for entering the Tonks-Girardeau regime does not exist. Furthermore, we discuss the computational efficiency of two quantum Monte Carlo methods to simulate large scale trapped bosonic systems: directed loops in stochastic series expansions and the worm algorithm.