Contributor info: Prof. Kai Phillip Schmidt

J. A. Koziol, G. Morigi, K. P. Schmidt

SciPost Phys. 17, 111 (2024) · published 10 October 2024 | Toggle abstract · pdf

We analyse the ground-state quantum phase diagram of hardcore Bosons interacting with repulsive dipolar potentials. The bosons dynamics is described by the extended-Bose-Hubbard Hamiltonian on a two-dimensional lattice. The ground state results from the interplay between the lattice geometry and the long-range interactions, which we account for by means of a classical spin mean-field approach limited by the size of the considered unit cells. This extended classical spin mean-field theory accounts for the long-range density-density interaction without truncation. We consider three different lattice geometries: square, honeycomb, and triangular. In the limit of zero hopping the ground state is always a devil's staircase of solid (gapped) phases. Such crystalline phases with broken translational symmetry are robust with respect to finite hopping amplitudes. At intermediate hopping amplitudes, these gapped phases melt, giving rise to various lattice supersolid phases, which can have exotic features with multiple sublattice densities. At sufficiently large hoppings the ground state is a superfluid. The stability of phases predicted by our approach is gauged by comparison to the known quantum phase diagrams of the Bose-Hubbard model with nearest-neighbour interactions as well as quantum Monte Carlo simulations for the dipolar case on the square and triangular lattice. Our results are of immediate relevance for experimental realisations of self-organised crystalline ordering patterns in analogue quantum simulators, e.g., with ultracold dipolar atoms in an optical lattice.

L. Schamriß, M. R. Walther, K. P. Schmidt

SciPost Phys. 17, 094 (2024) · published 26 September 2024 | Toggle abstract · pdf

Directly evaluated enhanced perturbative continuous unitary transformations (deepCUTs) are used to calculate non-perturbatively extrapolated numerical data for the ground-state energy and the energy gap. The data coincide with numerical evaluations of the truncated perturbative series and provide robust extrapolations beyond the perturbative regime. We develop a general scheme to extract quantum-critical properties from the deepCUT data based on critical scaling and a strict correspondence between the truncation used for deepCUT and the length scale of correlations at the critical point. We apply our approach to transverse-field Ising models (TFIMs) as paradigmatic systems for quantum phase transitions of various universality classes depending on the lattice geometry and the choice of antiferromagnetic or ferromagnetic coupling. In particular, we focus on the quantum phase diagram of the bilayer antiferromagnetic TFIM on the triangular lattice with an Ising-type interlayer coupling. Without a field, the model is known to host a classically disordered ground state, and in the limit of decoupled layers it exhibits the 3d-XY 'order by disorder' transition of the corresponding single-layer model. Our starting point for the unknown parts of the phase diagram is a high-order perturbative calculation about the limit of isolated dimers where the model is in a gapped phase.

Calvin Krämer, Jan Alexander Koziol, Anja Langheld, Max Hörmann, Kai Phillip Schmidt

SciPost Phys. 17, 061 (2024) · published 23 August 2024 | Toggle abstract · pdf

We study the ferromagnetic transverse-field Ising model with quenched disorder at $T = 0$ in one and two dimensions by means of stochastic series expansion quantum Monte Carlo simulations using a rigorous zero-temperature scheme. Using a sample-replication method and averaged Binder ratios, we determine the critical shift and width exponents $\nu_\mathrm{s}$ and $\nu_\mathrm{w}$ as well as unbiased critical points by finite-size scaling. Further, scaling of the disorder-averaged magnetisation at the critical point is used to determine the order-parameter critical exponent $\beta$ and the critical exponent $\nu_{\mathrm{av}}$ of the average correlation length. The dynamic scaling in the Griffiths phase is investigated by measuring the local susceptibility in the disordered phase and the dynamic exponent $z'$ is extracted. By applying various finite-size scaling protocols, we provide an extensive and comprehensive comparison between the different approaches on equal footing. The emphasis on effective zero-temperature simulations resolves several inconsistencies in existing literature.

Viktor Kott, Matthias Mühlhauser, Jan Alexander Koziol, Kai Phillip Schmidt

SciPost Phys. 17, 053 (2024) · published 15 August 2024 | Toggle abstract · pdf

We investigate the quantum robustness of the topological order in the toric code on the honeycomb lattice in the presence of a uniform parallel field. For a field in $z$-direction, the low-energy physics is in the flux-free sector and can be mapped to the transverse-field Ising model on the honeycomb lattice. One finds a second-order quantum phase transition in the 3D Ising$^\star$ universality class for both signs of the field. The same is true for a postive field in $x$-direction where an analogue mapping in the charge-free sector yields a ferromagnetic transverse-field Ising model on the triangular lattice and the phase transition is still 3D Ising$^\star$. In contrast, for negative $x$-field, the charge-free sector is mapped to the highly frustrated antiferromagnetic transverse-field Ising model on the triangular lattice which is known to host a quantum phase transition in the 3D XY$^\star$ universality class. Further, the charge-free sector does not always contain the low-energy physics for negative $x$-fields and a first-order phase transition to the polarized phase in the charge-full sector takes place at larger negative field values. We quantify the location of this transition by comparing quantum Monte Carlo simulations and high-field series expansions. The full extension of the topological phase in the presence of $x$- and $z$-fields is determined by perturbative linked-cluster expansions using a full graph decomposition. Extrapolating the high-order series of the charge and the flux gap allows to estimate critical exponents of the gap closing. This analysis indicates that the topological order breaks down by critical lines of 3D Ising$^\star$ and 3D XY$^\star$ type with interesting potential multi-critical crossing points. All findings for the toric code on the honeycomb lattice can be transferred exactly to the toric code on the triangular lattice.

Andreas Schellenberger, Kai Phillip Schmidt

SciPost Phys. Core 7, 038 (2024) · published 2 July 2024 | Toggle abstract · pdf

We investigate classes of interacting quantum spin systems in a single-mode cavity with a Dicke coupling, as a paradigmatic example of strongly correlated light-matter systems. Coming from the limit of weak light-matter couplings and large number of matter entities, we map the relevant low-energy sector of a broad class of models in the non-superradiant phases onto the exactly solvable Dicke model. We apply the outcomes to the Dicke-Ising model as a paradigmatic example [Phys. Rev. Res. 2, 023131 (2020), Sci. Rep. 4, 4083 (2014)], in agreement with results obtained by mean-field theory [Sci. Rep. 4, 4083 (2014)]. We further accompany and verify our findings with finite-size calculations, using exact diagonalization and the series expansion method pcst$^{++}$ [Phys. Rev. A 108, 013323 (2023)]

Matthias Mühlhauser, Viktor Kott, Kai Phillip Schmidt

SciPost Phys. Core 7, 031 (2024) · published 30 May 2024 | Toggle abstract · pdf

In this work we describe how to systematically implement a full graph decomposition to set up a linked-cluster expansion for the topological phase of Kitaev's toric code in a field. This demands to include the non-local effects mediated by the mutual anyonic statistics of elementary charge and flux excitations. Technically, we describe how to consistently integrate such non-local effects into a hypergraph decomposition for single excitations. The approach is demonstrated for the ground-state energy and the elementary excitation energies of charges and fluxes in the perturbed topological phase.

Max Hörmann, Kai P. Schmidt

SciPost Phys. 15, 097 (2023) · published 15 September 2023 | Toggle abstract · pdf

We construct a projection-based cluster-additive transformation that block-diagonalizes wide classes of lattice Hamiltonians $\mathcal{H}=\mathcal{H}_0 +V$. Its cluster additivity is an essential ingredient to set up perturbative or non-perturbative linked-cluster expansions for degenerate excitation subspaces of $\mathcal{H}_0$. Our transformation generalizes the minimal transformation known amongst others under the names Takahashi's transformation, Schrieffer-Wolff transformation, des Cloiseaux effective Hamiltonian, canonical van Vleck effective Hamiltonian or two-block orthogonalization method. The effective cluster-additive Hamiltonian and the transformation for a given subspace of $\mathcal{H}$, that is adiabatically connected to the eigenspace of $\mathcal{H}_0$ with eigenvalue $e_0^n$, solely depends on the eigenspaces of $\mathcal{H}$ connected to $e_0^m$ with $e_0^m≤ e_0^n$. In contrast, other cluster-additive transformations like the multi-block orthogonalization method or perturbative continuous unitary transformations need a larger basis. This can be exploited to implement the transformation efficiently both perturbatively and non-perturbatively. As a benchmark, we perform perturbative and non-perturbative linked-cluster expansions in the low-field ordered phase of the transverse-field Ising model on the square lattice for single spin-flips and two spin-flip bound-states.

Patrick Adelhardt, Kai Phillip Schmidt

SciPost Phys. 15, 087 (2023) · published 11 September 2023 | Toggle abstract · pdf

We investigate the quantum-critical behavior between the rung-singlet phase with hidden string order and the Néel phase with broken SU(2)-symmetry in quantum spin ladders with algebraically decaying unfrustrated long-range Heisenberg interactions. To this end, we determine high-order series expansions of energies and observables in the thermodynamic limit about the isolated rung-dimer limit. This is achieved by extending the method of perturbative continuous unitary transformations (pCUT) to long-range Heisenberg interactions and to the calculation of generic observables. The quantum-critical breakdown of the rung-singlet phase then allows us to determine the critical phase transition line and the entire set of critical exponents as a function of the decay exponent of the long-range interaction. We demonstrate long-range mean-field behavior as well as a non-trivial regime of continuously varying critical exponents implying the absence of deconfined criticality contrary to a recent suggestion in the literature.

Jan A. Koziol, Antonia Duft, Giovanna Morigi, Kai P. Schmidt

SciPost Phys. 14, 136 (2023) · published 30 May 2023 | Toggle abstract · pdf

We propose a general approach to analyse diagonal ordering patterns in bosonic lattice models with algebraically decaying density-density interactions on arbitrary lattices. The key idea is a systematic search for the energetically best order on all unit cells of the lattice up to a given extent. Using resummed couplings we evaluate the energy of the ordering patterns in the thermodynamic limit using finite unit cells. We apply the proposed approach to the atomic limit of the extended Bose-Hubbard model on the triangular lattice at fillings $f=1/2$ and $f=1$. We investigate the ground-state properties of the antiferromagnetic long-range Ising model on the triangular lattice and determine a six-fold degenerate plain-stripe phase to be the ground state for finite decay exponents. We also probe the classical limit of the Fendley-Sengupta-Sachdev model describing Rydberg atom arrays. We focus on arrangements where the atoms are placed on the sites or links of the Kagome lattice. Our method provides a general framework to treat cristalline structures resulting from long-range interactions.

A. Langheld, J. A. Koziol, P. Adelhardt, S. C. Kapfer, K. P. Schmidt

SciPost Phys. 13, 088 (2022) · published 10 October 2022 | Toggle abstract · pdf

The hyperscaling relation and standard finite-size scaling (FSS) are known to break down above the upper critical dimension due to dangerous irrelevant variables. We establish a coherent formalism for FSS at quantum phase transitions above the upper critical dimension following the recently introduced Q-FSS formalism for thermal phase transitions. Contrary to long-standing belief, the correlation sector is affected by dangerous irrelevant variables. The presented formalism recovers a generalized hyperscaling relation and FSS form. Using this new FSS formalism, we determine the full set of critical exponents for the long-range transverse-field Ising chain in all criticality regimes ranging from the nearest-neighbor to the long-range mean field regime. For the same model, we also explicitly confirm the effect of dangerous irrelevant variables on the characteristic length scale.

M. Mühlhauser, K. P. Schmidt , J. Vidal, M. R. Walther

SciPost Phys. 12, 069 (2022) · published 22 February 2022 | Toggle abstract · pdf

We study the competition between two different topological orders in three dimensions by considering the X-cube model and the three-dimensional toric code. The corresponding Hamiltonian can be decomposed into two commuting parts, one of which displaying a self-dual spectrum. To determine the phase diagram, we compute the high-order series expansions of the ground-state energy in all limiting cases. Apart from the topological order related to the toric code and the fractonic order related to the X-cube model, we found two new phases which are adiabatically connected to classical limits with nontrivial sub-extensive degeneracies. All phase transitions are found to be first order.

David A. Reiss, Kai Phillip Schmidt

SciPost Phys. 6, 078 (2019) · published 27 June 2019 | Toggle abstract · pdf

We study the robustness of 3D intrinsic topogical order under external perturbations by investigating the paradigmatic microscopic model, the 3D toric code in an external magnetic field. Exact dualities as well as variational calculations reveal a ground-state phase diagram with first and second-order quantum phase transitions. The variational approach can be applied without further approximations only for certain field directions. In the general field case, an approximative scheme based on an expansion of the variational energy in orders of the variational parameters is developed. For the breakdown of the 3D intrinsic topological order, it is found that the (im-)mobility of the quasiparticle excitations is crucial in contrast to their fractional statistics.

Michael Powalski, Kai P. Schmidt, Götz S. Uhrig

SciPost Phys. 4, 001 (2018) · published 17 January 2018 | Toggle abstract · pdf

The Heisenberg model for S=1/2 describes the interacting spins of electrons localized on lattice sites due to strong repulsion. It is the simplest strong-coupling model in condensed matter physics with wide-spread applications. Its relevance has been boosted further by the discovery of curate high-temperature superconductors. In leading order, their undoped parent compounds realize the Heisenberg model on square-lattices. Much is known about the model, but mostly at small wave vectors, i.e., for long-range processes, where the physics is governed by spin waves (magnons), the Goldstone bosons of the long-range ordered antiferromagnetic phase. Much less, however, is known for short-range processes, i.e., at large wave vectors. Yet these processes are decisive for understanding high-temperature superconductivity. Recent reports suggest that one has to resort to qualitatively different fractional excitations, spinons. By contrast, we present a comprehensive picture in terms of dressed magnons with strong mutual attraction on short length scales. The resulting spectral signatures agree strikingly with experimental data.