Showing 1–50 of 137 results for author: Schmidt, K P

Quantum melting of long-range ordered quantum antiferromagnets investigated by momentum-space continuous similarity transformations

Authors: Dag-Björn Hering, Matthias R. Walther, Kai P. Schmidt, Götz S. Uhrig

Abstract: We apply continuous similarity transformations (CSTs) to study the zero-temperature breakdown of long-range ordered quantum antiferromagnets. The CST flow equations are truncated in momentum space by the scaling dimension so that all operators with scaling dimension up to two are taken into account. We determine the quantum phase transition out of the Néel-ordered phase in the unfrustrated square… ▽ More We apply continuous similarity transformations (CSTs) to study the zero-temperature breakdown of long-range ordered quantum antiferromagnets. The CST flow equations are truncated in momentum space by the scaling dimension so that all operators with scaling dimension up to two are taken into account. We determine the quantum phase transition out of the Néel-ordered phase in the unfrustrated square lattice Heisenberg bilayer as well as the quantum melting of the Néel-ordered and columnar phase in the highly frustrated $J_1$-$J_2$ model on the square lattice. In all cases the CST is set up to isolate the ground state so that the stability of the flow equations, the ground-state energy, and the sublattice magnetization are used to explore the long-range ordered phases. We extract quantum-critical points which agree well with values in the literature. Further, we estimate the associated critical exponents $α$ and $β$ which turns out to be a challenging task for the CST approach. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: 16 pages, 15 figures

Quantitative description of long-range order in the anisotropic spin-1/2 Heisenberg antiferromagnet on the square lattice

Authors: Nils Caci, Dag-Björn Hering, Matthias R. Walther, Kai P. Schmidt, Stefan Wessel, Götz S. Uhrig

Abstract: The quantitative description of long-range order remains a challenge in quantum many-body physics. We provide zero-temperature results from two complementary methods for the ground-state energy per site, the sublattice magnetization, the spin gap, and the transverse spin correlation length for the spin-1/2 anisotropic quantum Heisenberg antiferromagnet on the square lattice. On the one hand, we us… ▽ More The quantitative description of long-range order remains a challenge in quantum many-body physics. We provide zero-temperature results from two complementary methods for the ground-state energy per site, the sublattice magnetization, the spin gap, and the transverse spin correlation length for the spin-1/2 anisotropic quantum Heisenberg antiferromagnet on the square lattice. On the one hand, we use exact, large-scale quantum Monte Carlo (QMC) simulations. On the other hand, we use the semi-analytic approach based on continuous similarity transformations in terms of elementary magnon excitations. Our findings confirm the applicability and quantitative validity of both approaches along the full parameter axis from the Ising point to the symmetry-restoring phase transition at the Heisenberg point and further provide quantitative reference results in the thermodynamic limit. In addition, we analytically derive the relation between the dispersion and the correlation length at zero temperature in arbitrary dimension, and discuss improved second-moment QMC estimators. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Comments: 12 pages, 5 figures

Quantum-critical properties of the one- and two-dimensional random transverse-field Ising model from large-scale quantum Monte Carlo simulations

Authors: C. Krämer, J. A. Koziol, A. Langheld, M. Hörmann, K. P. Schmidt

Abstract: 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 $ν_\mathrm{s}$ and $ν_\mathrm{w}$ as well as unbiased cr… ▽ More 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 $ν_\mathrm{s}$ and $ν_\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 $β$ and the critical exponent $ν_{\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. △ Less

Submitted 5 June, 2024; v1 submitted 8 March, 2024; originally announced March 2024.

Comments: 46 pages, 28 figures

Monte Carlo based techniques for quantum magnets with long-range interactions

Authors: P. Adelhardt, J. A. Koziol, A. Langheld, K. P. Schmidt

Abstract: Long-range interactions are relevant for a large variety of quantum systems in quantum optics and condensed matter physics. In particular, the control of quantum-optical platforms promises to gain deep insights in quantum-critical properties induced by the long-range nature of interactions. From a theoretical perspective, long-range interactions are notoriously complicated to treat. Here, we give… ▽ More Long-range interactions are relevant for a large variety of quantum systems in quantum optics and condensed matter physics. In particular, the control of quantum-optical platforms promises to gain deep insights in quantum-critical properties induced by the long-range nature of interactions. From a theoretical perspective, long-range interactions are notoriously complicated to treat. Here, we give an overview of recent advancements to investigate quantum magnets with long-range interactions focusing on two techniques based on Monte Carlo integration. First, the method of perturbative continuous unitary transformations where classical Monte Carlo integration is applied within the embedding scheme of white graphs. This linked-cluster expansion allows to extract high-order series expansions of energies and observables in the thermodynamic limit. Second, stochastic series expansion quantum Monte Carlo which enables calculations on large finite systems. Finite-size scaling can then be used to determine physical properties of the infinite system. In recent years, both techniques have been applied successfully to one- and two-dimensional quantum magnets involving long-range Ising, XY, and Heisenberg interactions on various bipartite and non-bipartite lattices. Here, we summarise the obtained quantum-critical properties including critical exponents for all these systems in a coherent way. Further, we review how long-range interactions are used to study quantum phase transitions above the upper critical dimension and the scaling techniques to extract these quantum critical properties from the numerical calculations. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Comments: 141 pages, 38 figures

Extracting quantum-critical properties from directly evaluated enhanced perturbative continuous unitary transformations

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

Abstract: 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 coincides with the perturbative series up to the order with respect to which the deepCUT is truncated. We develop a general scheme to extract quantum-critical properties from the deepCUT… ▽ More 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 coincides with the perturbative series up to the order with respect to which the deepCUT is truncated. 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. △ Less

Submitted 29 February, 2024; originally announced February 2024.

Comments: 31 pages

Quantum robustness of the toric code in a parallel field on the honeycomb and triangular lattice

Authors: V. Kott, M. Mühlhauser, J. A. Koziol, K. P. Schmidt

Abstract: 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 cl… ▽ More 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. △ Less

Submitted 10 June, 2024; v1 submitted 23 February, 2024; originally announced February 2024.

Comments: 26 pages, 8 figures

(Almost) Everything is a Dicke model -- Mapping non-superradiant correlated light-matter systems to the exactly solvable Dicke model

Authors: Andreas Schellenberger, Kai Phillip Schmidt

Abstract: 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… ▽ More 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, in agreement with results obtained by mean-field theory. We further accompany and verify our findings with finite-size calculations, using exact diagonalization and the series expansion method pcst++. △ Less

Submitted 2 May, 2024; v1 submitted 23 February, 2024; originally announced February 2024.

Comments: 26 pages, 7 figures; Added proposed changes by the editors during SciPost review process

Order-by-disorder and long-range interactions in the antiferromagnetic transverse-field Ising model on the triangular lattice -- A perturbative point of view

Authors: J. A. Koziol, M. Mühlhauser, K. P. Schmidt

Abstract: We study the low-field ground-state (GS) properties of the antiferromagnetic transverse-field Ising model with long-range interactions (afLRTFIM) on the triangular lattice. We use the method of perturbative continuous unitary transformations (pCUT) to derive an effective model for the degenerate GS space of the antiferromagnetic nearest-neighbour (NN) Ising model on a finite system, by treating th… ▽ More We study the low-field ground-state (GS) properties of the antiferromagnetic transverse-field Ising model with long-range interactions (afLRTFIM) on the triangular lattice. We use the method of perturbative continuous unitary transformations (pCUT) to derive an effective model for the degenerate GS space of the antiferromagnetic nearest-neighbour (NN) Ising model on a finite system, by treating the transverse-field (TF) and the long-range interactions (LRI) as a perturbation. We determine a level-crossing between the plain stripe phase at small TF and the clock-ordered phase at intermediate TF at $h\cong0.129$ for $α=6$, $N=36$ spins in order three perturbation theory. We discuss the qualitative layout of the quantum phase diagram of the afLRTFIM on the triangular lattice. △ Less

Submitted 16 February, 2024; originally announced February 2024.

Comments: 3 pages, 1 figure

Incorporating non-local anyonic statistics into a graph decomposition

Authors: M. Mühlhauser, V. Kott, K. P. Schmidt

Abstract: 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 hy… ▽ More 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. △ Less

Submitted 12 January, 2024; originally announced January 2024.

Comments: 29 pages

Order-by-disorder in the antiferromagnetic $J_1$-$J_2$-$J_3$ transverse-field Ising model on the ruby lattice

Authors: A. Duft, J. A. Koziol, P. Adelhardt, M. Mühlhauser, K. P. Schmidt

Abstract: We investigate the quantum phase diagram of the $J_1$-$J_2$-$J_3$ antiferromagnetic transverse-field Ising model on the ruby lattice. In the low-field limit we derive an effective quantum dimer model, analyzing how the extensive ground-state degeneracy at zero field is lifted by an order-by-disorder scenario. We support our analysis by studying the gap-closing of the high-field phase using series… ▽ More We investigate the quantum phase diagram of the $J_1$-$J_2$-$J_3$ antiferromagnetic transverse-field Ising model on the ruby lattice. In the low-field limit we derive an effective quantum dimer model, analyzing how the extensive ground-state degeneracy at zero field is lifted by an order-by-disorder scenario. We support our analysis by studying the gap-closing of the high-field phase using series expansions. For $J_2>J_3$, we find a columnar phase at low fields, followed by a clock-ordered phase stabilized by resonating plaquettes at intermediate field values, and an emergent 3d-XY quantum phase transition to the polarized high-field phase. For $J_3>J_2$, an order-by-disorder mechanism stabilizes a distinct $k=(0,0)$ order and a quantum phase transition in the 3d-Ising universality class is observed. Further, we discuss the possible implementation of the columnar- and clock-ordered phase in existing Rydberg atom quantum simulators. When taking into account the full algebraically decaying long-range interactions on the ruby lattice, we find that long-range interactions favor the same ground state as the quantum fluctuations induced by a transverse field, which could make the ruby lattice a promising candidate for the realization of a clock-ordered phase. △ Less

Submitted 13 May, 2024; v1 submitted 20 December, 2023; originally announced December 2023.

Comments: 23 pages, 16 figures

Quantum phases of hardcore bosons with repulsive dipolar density-density interactions on two-dimensional lattices

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

Abstract: 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 limi… ▽ More 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. △ Less

Submitted 10 July, 2024; v1 submitted 17 November, 2023; originally announced November 2023.

Comments: 33 pages, 10 figures

Engineering Photon-mediated Long-Range Spin Interactions in Mott Insulators

Authors: Paul Fadler, Jiajun Li, Kai Phillip Schmidt, Martin Eckstein

Abstract: We investigate the potential to induce long-range spin interactions in a Mott insulator via the quantum electromagnetic field of a cavity. The coupling between light and spins is inherently non-linear, and occurs via multi-photon processes like Raman scattering and two-photon absorption/emission with electronically excited intermediate states. Based on this, two pathways are elucidated: (i) In the… ▽ More We investigate the potential to induce long-range spin interactions in a Mott insulator via the quantum electromagnetic field of a cavity. The coupling between light and spins is inherently non-linear, and occurs via multi-photon processes like Raman scattering and two-photon absorption/emission with electronically excited intermediate states. Based on this, two pathways are elucidated: (i) In the absence of external driving, long-range interactions are mediated by the exchange of at least two virtual cavity photons. We show that these vacuum-mediated interactions can surpass local Heisenberg interactions in mesoscopic setups such as sufficiently small split-ring resonators. (ii) In a laser-driven cavity, interactions can be tailored through a hybrid scheme involving both external laser photons and cavity photons. This offers a versatile pathway for Floquet engineering of long-range interactions in macroscopic systems. In general, the derivation of these interactions requires careful consideration: Notably, we demonstrate that a simple phenomenological approach, based on a spin-photon Hamiltonian that captures Raman and two-photon processes with effective matrix elements, can be used only if the cavity is resonantly driven. Outside of these narrow resonant regimes as well as for the undriven case, a fourth-order series expansion within the underlying electronic model is necessary, which we perform to obtain long-range four-spin interactions in the half-filled Hubbard model. △ Less

Submitted 2 November, 2023; originally announced November 2023.

Comments: 25 pages, 9 figures

Journal ref: Phys. Rev. B 109, 085149 (2024)

Hybrid quantum-classical algorithm for the transverse-field Ising model in the thermodynamic limit

Authors: Sumeet, M. Hörmann, K. P. Schmidt

Abstract: We describe a hybrid quantum-classical approach to treat quantum many-body systems in the thermodynamic limit. This is done by combining numerical linked-cluster expansions (NLCE) with the variational quantum eigensolver (VQE). Here, the VQE algorithm is used as a cluster solver within the NLCE. We test our hybrid quantum-classical algorithm (NLCE$+$VQE) for the ferromagnetic transverse-field Isin… ▽ More We describe a hybrid quantum-classical approach to treat quantum many-body systems in the thermodynamic limit. This is done by combining numerical linked-cluster expansions (NLCE) with the variational quantum eigensolver (VQE). Here, the VQE algorithm is used as a cluster solver within the NLCE. We test our hybrid quantum-classical algorithm (NLCE$+$VQE) for the ferromagnetic transverse-field Ising model on the one-dimensional chain and the two-dimensional square lattice. The calculation of ground-state energies on each open cluster demands a modified Hamiltonian variational ansatz for the VQE. One major finding is convergence of NLCE$+$VQE to the conventional NLCE result in the thermodynamic limit when at least $N/2$ layers are used in the VQE ansatz for each cluster with $N$ sites. Our approach demonstrates the fruitful connection of techniques known from correlated quantum many-body systems with hybrid algorithms explored on existing quantum-computing devices. △ Less

Submitted 11 October, 2023; originally announced October 2023.

Comments: 14 pages, 7 figures

Cooperative effects in dense cold atomic gases including magnetic dipole interactions

Authors: N. S. Bassler, I. Varma, M. Proske, P. Windpassinger, K. P. Schmidt, C. Genes

Abstract: We theoretically investigate cooperative effects in cold atomic gases exhibiting both electric and magnetic dipole-dipole interactions, such as occurring for example in clouds of dysprosium atoms. We distinguish between the quantum degenerate case, where we take a many body physics approach and the quantum non-degenerate case, where we use the formalism of open system dynamics. For quantum non-deg… ▽ More We theoretically investigate cooperative effects in cold atomic gases exhibiting both electric and magnetic dipole-dipole interactions, such as occurring for example in clouds of dysprosium atoms. We distinguish between the quantum degenerate case, where we take a many body physics approach and the quantum non-degenerate case, where we use the formalism of open system dynamics. For quantum non-degenerate gases, we illustrate the emergence of tailorable spin models in the high-excitation limit. In the low-excitation limit, we provide analytical and numerical results detailing the effect of magnetic interactions on the directionality of scattered light and characterize sub- and superradiant effects. For quantum degenerate gases, we study the interplay between sub- and superradiance effects and the fermionic or bosonic quantum statistics nature of the ensemble. △ Less

Submitted 20 June, 2023; originally announced June 2023.

Journal ref: Phys. Rev. Research 6, 023147 (2024)

Metasurface-based hybrid optical cavities for chiral sensing

Authors: Nico S. Bassler, Andrea Aiello, Kai P. Schmidt, Claudiu Genes, Michael Reitz

Abstract: Quantum metasurfaces, i.e., two-dimensional subwavelength arrays of quantum emitters, can be employed as mirrors towards the design of hybrid cavities, where the optical response is given by the interplay of a cavity-confined field and the surface modes supported by the arrays. We show that, under external magnetic field control, stacked layers of quantum metasurfaces can serve as helicity-preserv… ▽ More Quantum metasurfaces, i.e., two-dimensional subwavelength arrays of quantum emitters, can be employed as mirrors towards the design of hybrid cavities, where the optical response is given by the interplay of a cavity-confined field and the surface modes supported by the arrays. We show that, under external magnetic field control, stacked layers of quantum metasurfaces can serve as helicity-preserving cavities. These structures exhibit ultranarrow resonances and can enhance the intensity of the incoming field by orders of magnitude, while simultaneously preserving the handedness of the field circulating inside the resonator, as opposed to conventional cavities. The rapid phase shift in the cavity transmission around the resonance can be exploited for the sensitive detection of chiral scatterers passing through the cavity. We discuss possible applications of these resonators as sensors for the discrimination of chiral molecules. △ Less

Submitted 1 June, 2023; originally announced June 2023.

Comments: 6+11 pages, 2+3 figures

Projective cluster-additive transformation for quantum lattice models

Authors: M. Hörmann, K. P. Schmidt

Abstract: 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 amon… ▽ More 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\leq e_0^n$. In contrast, other cluster-additive transformations like the multi-block orthognalization 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. △ Less

Submitted 8 March, 2023; originally announced March 2023.

Comments: 28 pages, 5 figures

Journal ref: SciPost Phys. 15, 097 (2023)

Absence of fractal quantum criticality in the quantum Newman-Moore model

Authors: R. Wiedmann, L. Lenke., M. Mühlhauser, K. P. Schmidt

Abstract: The quantum phase transition between the low-field fracton phase with type-II fracton excitations and the high-field polarized phase is investigated in the two-dimensional self-dual quantum Newman-Moore model. We apply perturbative and numerical linked-cluster expansions to calculate the ground-state energy per site in the thermodynamic limit revealing a level crossing at the self-dual point. In a… ▽ More The quantum phase transition between the low-field fracton phase with type-II fracton excitations and the high-field polarized phase is investigated in the two-dimensional self-dual quantum Newman-Moore model. We apply perturbative and numerical linked-cluster expansions to calculate the ground-state energy per site in the thermodynamic limit revealing a level crossing at the self-dual point. In addition, high-order series expansions of the relevant low-energy gaps are determined using perturbative continuous unitary transformations indicating no gap-closing. Our results therefore predict a first-order phase transition between the low-field fracton and the high-field polarized phase at the self-dual point. △ Less

Submitted 3 February, 2023; originally announced February 2023.

Comments: 8 pages, 4 figures

Series expansions in closed and open quantum many-body systems with multiple quasiparticle types

Authors: L. Lenke, A. Schellenberger, K. P. Schmidt

Abstract: The established approach of perturbative continuous unitary transformations (pCUTs) constructs effective quantum many-body Hamiltonians as perturbative series that conserve the number of one quasiparticle type. We extend the pCUT method to similarity transformations - dubbed $\mathrm{pcst}^{\texttt{++}}$ - allowing for multiple quasiparticle types with complex-valued energies. This enlarges the fi… ▽ More The established approach of perturbative continuous unitary transformations (pCUTs) constructs effective quantum many-body Hamiltonians as perturbative series that conserve the number of one quasiparticle type. We extend the pCUT method to similarity transformations - dubbed $\mathrm{pcst}^{\texttt{++}}$ - allowing for multiple quasiparticle types with complex-valued energies. This enlarges the field of application to closed and open quantum many-body systems with unperturbed operators corresponding to arbitrary superimposed ladder spectra. To this end, a generalized counting operator is combined with the quasiparticle generator for open quantum systems recently introduced by Schmiedinghoff and Uhrig (arXiv:2203.15532). The $\mathrm{pcst}^{\texttt{++}}$ then yields model-independent block-diagonal effective Hamiltonians and Lindbladians allowing a linked-cluster expansion in the thermodynamic limit similar to the conventional pCUT method. We illustrate the application of the $\mathrm{pcst}^{\texttt{++}}$ method by discussing representative closed, open, and non-Hermitian quantum systems. △ Less

Submitted 4 October, 2023; v1 submitted 2 February, 2023; originally announced February 2023.

Comments: 12 pages, 4 figures

Journal ref: Phys. Rev. A 108, 013323 (2023)

Systematic Analysis of Crystalline Phases in Bosonic Lattice Models with Algebraically Decaying Density-Density Interactions

Authors: J. A. Koziol, A. Duft, G. Morigi, K. P. Schmidt

Abstract: 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 us… ▽ More 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. \changed{Our method provides a general framework to treat cristalline structures resulting from long-range interactions. △ Less

Submitted 4 March, 2023; v1 submitted 5 December, 2022; originally announced December 2022.

Comments: 35 pages, 11 figures

Journal ref: SciPost Phys. 14, 136 (2023)

Continuous similarity transformation for critical phenomena: easy-axis antiferromagnetic XXZ model

Authors: Matthias R. Walther, Dag-Björn Hering, Götz S. Uhrig, Kai P. Schmidt

Abstract: We apply continuous similarity transformations (CSTs) to the easy-axis antiferromagnetic XXZ-model on the square lattice. The CST flow equations are truncated in momentum space by the scaling dimension $d$ so that all contributions with $d\le 2$ are taken into account. The resulting quartic magnon-conserving effective Hamiltonian is analyzed in the zero-, one-, and two-magnon sector. In this way,… ▽ More We apply continuous similarity transformations (CSTs) to the easy-axis antiferromagnetic XXZ-model on the square lattice. The CST flow equations are truncated in momentum space by the scaling dimension $d$ so that all contributions with $d\le 2$ are taken into account. The resulting quartic magnon-conserving effective Hamiltonian is analyzed in the zero-, one-, and two-magnon sector. In this way, a quantitative description of the ground-state energy, the one-magnon dispersion and its gap as well as of two-magnon bound states is gained for anisotropies ranging from the gapped Ising model to the gapless Heisenberg model. We discuss the critical properties of the gap closing as well as the evolution of the one-magnon roton mininum. The excitation energies of two-magnon bound states are calculated and their decay into the two-magnon continuum is determined via the inverse participation ratio. △ Less

Submitted 10 November, 2022; originally announced November 2022.

Comments: 12 pages, 8 figures

Journal ref: Phys. Rev. Res. 5, 013132 (2023)

Linear optical elements based on cooperative subwavelength emitter arrays

Authors: Nico S. Bassler, Michael Reitz, Kai P. Schmidt, Claudiu Genes

Abstract: We describe applications of two-dimensional subwavelength quantum emitter arrays as efficient optical elements in the linear regime. For normally incident light, the cooperative optical response, stemming from emitter-emitter dipole exchanges, allows the control of the array's transmission, its resonance frequency, and bandwidth. Operations on fully polarized incident light, such as generic linear… ▽ More We describe applications of two-dimensional subwavelength quantum emitter arrays as efficient optical elements in the linear regime. For normally incident light, the cooperative optical response, stemming from emitter-emitter dipole exchanges, allows the control of the array's transmission, its resonance frequency, and bandwidth. Operations on fully polarized incident light, such as generic linear and circular polarizers as well as phase retarders can be engineered and described in terms of Jones matrices. Our analytical approach and accompanying numerical simulations identify optimal regimes for such operations and reveal the importance of adjusting the array geometry and of the careful tuning of the external magnetic fields amplitude and direction. △ Less

Submitted 7 September, 2022; originally announced September 2022.

Comments: 12 pages main text + 4 Appendix, 8 figures

Continuously varying critical exponents in long-range quantum spin ladders

Authors: P. Adelhardt, K. P. Schmidt

Abstract: 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 on quantum spin ladders with algebraically decaying unfrustrated long-range Heisenberg interactions. Combining perturbative continuous unitary transformations (pCUT) with a white-graph expansion and Monte Carlo simulations yields high-order series exp… ▽ More 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 on quantum spin ladders with algebraically decaying unfrustrated long-range Heisenberg interactions. Combining perturbative continuous unitary transformations (pCUT) with a white-graph expansion and Monte Carlo simulations yields high-order series expansions of energies and observables in the thermodynamic limit about the isolated rung-dimer limit. The breakdown of the rung-singlet phase allows to determine the critical line and the entire set of critical exponents as a function of the decay exponent of the long-range interaction. A non-trivial regime of continuously varying critical exponents as well as long-range mean-field behavior is demonstrated reminiscent of the long-range transverse-field Ising model. △ Less

Submitted 2 September, 2022; originally announced September 2022.

Comments: 13 pages, 4 figures

Dynamic structure factor of the antiferromagnetic Kitaev model in large magnetic fields

Authors: A. Schellenberger, M. Hörmann, K. P. Schmidt

Abstract: We investigate the dynamic structure factor of the antiferromagnetic Kitaev honeycomb model in a magnetic field by applying perturbative continuous unitary transformations about the high-field limit. One- and two-quasiparticle properties of the dressed elementary spin-flip excitations of the high-field polarized phase are calculated which account for most of the spectral weight in the dynamic stru… ▽ More We investigate the dynamic structure factor of the antiferromagnetic Kitaev honeycomb model in a magnetic field by applying perturbative continuous unitary transformations about the high-field limit. One- and two-quasiparticle properties of the dressed elementary spin-flip excitations of the high-field polarized phase are calculated which account for most of the spectral weight in the dynamic structure factor. We discuss the evolution of spectral features in these quasiparticle sectors in terms of one-quasiparticle dispersions, two-quasiparticle continua, the formation of antibound states, and quasiparticle decay. In particular, a comparably strong spectral feature above the upper edge of the upmost two-quasiparticle continuum represents three antibound states which form due to nearest-neighbor density-density interactions. △ Less

Submitted 6 December, 2022; v1 submitted 25 March, 2022; originally announced March 2022.

Comments: 14 pages, 10 figures

Journal ref: Physical Review B 106, 104403 (2022)

Scaling at quantum phase transitions above the upper critical dimension

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

Abstract: 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… ▽ More 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. △ Less

Submitted 8 July, 2022; v1 submitted 15 March, 2022; originally announced March 2022.

Comments: 29 pages, 5 figures

Journal ref: SciPost Phys. 13, 088 (2022)

Linked Cluster Expansions via Hypergraph Decompositions

Authors: M. Mühlhauser, K. P. Schmidt

Abstract: We propose a hypergraph expansion which facilitates the direct treatment of quantum spin models with many-site interactions via perturbative linked cluster expansions. The main idea is to generate all relevant subclusters and sort them into equivalence classes essentially governed by hypergraph isomorphism. Concretely, a reduced König representation of the hypergraphs is used to make the equivalen… ▽ More We propose a hypergraph expansion which facilitates the direct treatment of quantum spin models with many-site interactions via perturbative linked cluster expansions. The main idea is to generate all relevant subclusters and sort them into equivalence classes essentially governed by hypergraph isomorphism. Concretely, a reduced König representation of the hypergraphs is used to make the equivalence relation accessible by graph isomorphism. During this procedure we determine the embedding factor for each equivalence class, which is used in the final resummation in order to obtain the final result. As an instructive example we calculate the ground-state energy and a particular excitation gap of the plaquette Ising model in a transverse field on the three-dimensional cubic lattice. △ Less

Submitted 7 February, 2022; originally announced February 2022.

Comments: 12 pages, 7 figures

Quantum phase transitions in the $K$-layer Ising toric code

Authors: L. Schamriss, L. Lenke, M. Mühlhauser, K. P. Schmidt

Abstract: We investigate the quantum phase diagram of the $K$-layer Ising toric code corresponding to $K$ layers of two-dimensional toric codes coupled by Ising interactions. While for small Ising interactions the system displays $\mathbb{Z}_2^K$ topological order originating from the toric codes in each layer, the system shows $\mathbb{Z}_2$ topological order in the high-Ising limit. The latter is demonstr… ▽ More We investigate the quantum phase diagram of the $K$-layer Ising toric code corresponding to $K$ layers of two-dimensional toric codes coupled by Ising interactions. While for small Ising interactions the system displays $\mathbb{Z}_2^K$ topological order originating from the toric codes in each layer, the system shows $\mathbb{Z}_2$ topological order in the high-Ising limit. The latter is demonstrated for general $K$ by deriving an effective low-energy model in $K^{\rm th}$-order degenerate perturbation theory, which is given as an effective anisotropic single-layer toric code in terms of collective pseudo-spins 1/2 refering to the two ground states of isolated Ising chain segments. For the specific cases $K=3$ and $K=4$ we apply high-order series expansions to determine the gap series in the low- and high-Ising limit. Extrapolation of the elementary energy gaps gives convincing evidence that the ground-state phase diagram consists of a single quantum critical point in the 3d Ising* universality class for both $K$ separating both types of topological order, which is consistent with former findings for the bilayer Ising toric code. △ Less

Submitted 25 January, 2022; originally announced January 2022.

Comments: 16 pages, 6 figures

Kitaev honeycomb antiferromagnet in a field: quantum phase diagram for general spin

Authors: Saeed S. Jahromi, Max Hörmann, Patrick Adelhardt, Sebastian Fey, Roman Orus, Kai Phillip Schmidt

Abstract: We combine tensor-network approaches and high-order linked-cluster expansions to investigate the quantum phase diagram of the antiferromagnetic Kitaev's honeycomb model in a magnetic field for general spin values. For the pure Kitaev model, tensor network calculations confirm the absence of fluxes and spin-spin correlations beyond nearest neighbor in the ground state, but signal a breaking of the… ▽ More We combine tensor-network approaches and high-order linked-cluster expansions to investigate the quantum phase diagram of the antiferromagnetic Kitaev's honeycomb model in a magnetic field for general spin values. For the pure Kitaev model, tensor network calculations confirm the absence of fluxes and spin-spin correlations beyond nearest neighbor in the ground state, but signal a breaking of the discrete orientational symmetry for $S\in\{1,3/2,2\}$ inline with the semiclassical limit. An intermediate region between Kitaev phases and the high-field polarized phase is demonstrated for all considered spin values. The analysis of the high-field gap and the associated spectral weight of the polarized phase for general spin S is consistent with an unconventional quantum critical breakdown of the high-field polarized phase in accordance with the presence of exotic physics at intermediate Kitaev couplings. △ Less

Submitted 11 November, 2021; originally announced November 2021.

Comments: 4 pages, 4 figures and supplementary material with 7 pages and 3 figure

High-order series expansion of non-Hermitian quantum spin models

Authors: L. Lenke, M. Mühlhauser, K. P. Schmidt

Abstract: We investigate the low-energy physics of non-Hermitian quantum spin models with $PT$-symmetry. To this end we consider the one-dimensional Ising chain and the two-dimensional toric code in a non-Hermitian staggered field. For both systems dual descriptions in terms of non-Hermitian staggered Ising interactions in a conventional transverse field exist. We perform high-order series expansions about… ▽ More We investigate the low-energy physics of non-Hermitian quantum spin models with $PT$-symmetry. To this end we consider the one-dimensional Ising chain and the two-dimensional toric code in a non-Hermitian staggered field. For both systems dual descriptions in terms of non-Hermitian staggered Ising interactions in a conventional transverse field exist. We perform high-order series expansions about the high- and low-field limit for both systems to determine the ground-state energy per site and the one-particle gap. The one-dimensional non-Hermitian Ising chain is known to be exactly solvable. Its ground-state phase diagram consists of second-order quantum phase transitions, which can be characterized by logarithmic singularities of the second derivative of the ground-state energy and, in the symmetry-broken phase, the gap closing of the low-field gap. In contrast, the gap closing from the high-field phase is not accessible perturbatively due to the complex energy and the occurrence of exceptional lines in the high-field gap expression. For the two-dimensional toric code in a non-Hermitian staggered field we study the quantum robustness of the topologically ordered phase by the gap closing of the low-field gap. We find that the well-known second-order quantum phase transition of the toric code in a uniform field extends into a large portion of the non-Hermitian parameter space. However, the series expansions become unreliable for a dominant anti-Hermitian field. Interestingly, the analysis of the high-field gap reveals the potential presence of an intermediate region. △ Less

Submitted 12 August, 2021; originally announced August 2021.

Comments: 15 pages, 9 figures

Competing topological orders in three dimensions: X-cube versus toric code

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

Abstract: 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 displays 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.… ▽ More 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 displays 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. △ Less

Submitted 18 March, 2022; v1 submitted 10 June, 2021; originally announced June 2021.

Comments: 14 pages, 5 figures

Journal ref: SciPost Physics 12, 069 (2022)

Quantum-critical properties of the long-range transverse-field Ising model from quantum Monte Carlo simulations

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

Abstract: The quantum-critical properties of the transverse-field Ising model with algebraically decaying interactions are investigated by means of stochastic series expansion quantum Monte Carlo, on both the one-dimensional linear chain and the two-dimensional square lattice. We extract the critical exponents $ν$ and $β$ as a function of the decay exponent of the long-range interactions. For ferromagnetic… ▽ More The quantum-critical properties of the transverse-field Ising model with algebraically decaying interactions are investigated by means of stochastic series expansion quantum Monte Carlo, on both the one-dimensional linear chain and the two-dimensional square lattice. We extract the critical exponents $ν$ and $β$ as a function of the decay exponent of the long-range interactions. For ferromagnetic Ising interactions, we resolve the limiting regimes known from field theory, ranging from the nearest-neighbor Ising to the long-range Gaussian universality classes, as well as the intermediate regime with continuously varying critical exponents. In the long-range Gaussian regime, we treat the effect of dangerous irrelevant variables on finite-size scaling forms. For antiferromagnetic and therefore competing Ising interactions, the stochastic series expansion algorithm displays growing auto-correlation times leading to a reduced performance. Nevertheless, our results are consistent with the nearest-neighbor Ising universality for all investigated interaction ranges both on the linear chain and the square lattice. △ Less

Submitted 17 March, 2021; originally announced March 2021.

Comments: 11 pages, 4 figures

Journal ref: Phys. Rev. B 103, 245135 (2021)

Evidence for pressure induced unconventional quantum criticality in the coupled spin ladder antiferromagnet C$_9$H$_{18}$N$_2$CuBr$_4$

Authors: Tao Hong, Tao Ying, Qing Huang, Sachith E. Dissanayake, Yiming Qiu, Mark M. Turnbull, Andrey A. Podlesnyak, Yan Wu, Huibo Cao, Yaohua Liu, Izuru Umehara, Jun Gouchi, Yoshiya Uwatoko, Masaaki Matsuda, David A. Tennant, Gia-Wei Chern, Kai P. Schmidt, Stefan Wessel

Abstract: Quantum phase transitions in quantum matter occur at zero temperature between distinct ground states by tuning a nonthermal control parameter. Often, they can be accurately described within the Landau theory of phase transitions, similarly to conventional thermal phase transitions. However, this picture can break down under certain circumstances. Here, we present a comprehensive study of the effec… ▽ More Quantum phase transitions in quantum matter occur at zero temperature between distinct ground states by tuning a nonthermal control parameter. Often, they can be accurately described within the Landau theory of phase transitions, similarly to conventional thermal phase transitions. However, this picture can break down under certain circumstances. Here, we present a comprehensive study of the effect of hydrostatic pressure on the magnetic structure and spin dynamics of the spin-1/2 ladder compound C$_9$H$_{18}$N$_2$CuBr$_4$. Single-crystal heat capacity and neutron diffraction measurements reveal that the N$\rm \acute{e}$el-ordered phase breaks down beyond a critical pressure of $P_{\rm c}$$\sim$1.0 GPa through a continuous quantum phase transition. Estimates of the critical exponents suggest that this transition may fall outside the traditional Landau paradigm. The inelastic neutron scattering spectra at 1.3 GPa are characterized by two well-separated gapped modes, including one continuum-like and another resolution-limited excitation in distinct scattering channels, which further indicates an exotic quantum-disordered phase above $P_{\rm c}$. △ Less

Submitted 19 May, 2022; v1 submitted 9 November, 2020; originally announced November 2020.

Comments: The revised version includes 10 pages and 5 figures in the main text. To appear in Nature Communications. Any comments are welcome

Journal ref: Nat. Commun. 13, 3073 (2022)

Quantum critical phase transition between two topologically-ordered phases in the Ising toric code bilayer

Authors: R. Wiedmann, L. Lenke, M. R. Walther, M. Mühlhauser, K. P. Schmidt

Abstract: We demonstrate that two toric code layers on the square lattice coupled by an Ising interaction display two distinct phases with intrinsic topological order. The second-order quantum phase transition between the weakly-coupled $\mathbb{Z}_2\times\mathbb{Z}_2$ and the strongly-coupled $\mathbb{Z}_2$ topological order can be described by the condensation of bosonic quasiparticles from both sides and… ▽ More We demonstrate that two toric code layers on the square lattice coupled by an Ising interaction display two distinct phases with intrinsic topological order. The second-order quantum phase transition between the weakly-coupled $\mathbb{Z}_2\times\mathbb{Z}_2$ and the strongly-coupled $\mathbb{Z}_2$ topological order can be described by the condensation of bosonic quasiparticles from both sides and belongs to the 3d Ising$^*$ universality class. This can be shown by an exact duality transformation to the transverse-field Ising model on the square lattice, which builds on the existence of an extensive number of local $\mathbb{Z}_2$ conserved parities. These conserved quantities correspond to the product of two adjacent star operators on different layers. Notably, we show that the low-energy effective model derived about the limit of large Ising coupling is given by an effective single-layer toric code in terms of the conserved quantities of the Ising toric code bilayer. The two topological phases are further characterized by the topological entanglement entropy which serves as a non-local order parameter. △ Less

Submitted 12 October, 2020; originally announced October 2020.

Comments: 6 pages, 3 figures

Journal ref: Phys. Rev. B 102, 214422 (2020)

Effects of domain walls in bilayer graphene in an external magnetic field

Authors: N. S. Bassler, K. P. Schmidt

Abstract: We investigate bilayer graphene systems with layer switching domain walls separating the two energetically equivalent Bernal stackings in the presence of an external magnetic field. To this end we calculate quantum transport and local densities of three microscopic models for a single domain wall: a hard wall, a defect due to shear, and a defect due to tension. The quantum transport calculations a… ▽ More We investigate bilayer graphene systems with layer switching domain walls separating the two energetically equivalent Bernal stackings in the presence of an external magnetic field. To this end we calculate quantum transport and local densities of three microscopic models for a single domain wall: a hard wall, a defect due to shear, and a defect due to tension. The quantum transport calculations are performed with a recursive Green's function method. Technically, we discuss an explicit algorithm for the separation of a system into subsystems for the recursion and we present an optimization of the well known iteration scheme for lead self-energies for sparse chain couplings. We find strong physical differences for the three different types of domain walls in the integer quantum Hall regime. For a domain wall due to shearing of the upper graphene layer there is a plateau formation in the magnetoconductance for sufficiently wide defect regions. For wide domain walls due to tension in the upper graphene layer there is only an approximate plateau formation with fluctuations of the order of the elementary conuctance quantum $σ_0$. A direct transition between stacking regions like for the hard wall domain wall shows no plateau formation and is therefore not a good model for either of the previously mentioned extended domain walls. △ Less

Submitted 14 September, 2020; originally announced September 2020.

Comments: 17 pages, 20 figures

Journal ref: Phys. Rev. B 103, 155422 (2021)

Quantum criticality and excitations of a long-range anisotropic $XY$-chain in a transverse field

Authors: P. Adelhardt, J. A. Koziol, A. Schellenberger, K. P. Schmidt

Abstract: The critical breakdown of a one-dimensional quantum magnet with long-range interactions is studied by investigating the high-field polarized phase of the anisotropic XY model in a transverse field for the ferro- and antiferromagnetic case. While for the limiting case of the isotropic long-range XY model we can extract the elementary one quasi-particle dispersion analytically and calculate two quas… ▽ More The critical breakdown of a one-dimensional quantum magnet with long-range interactions is studied by investigating the high-field polarized phase of the anisotropic XY model in a transverse field for the ferro- and antiferromagnetic case. While for the limiting case of the isotropic long-range XY model we can extract the elementary one quasi-particle dispersion analytically and calculate two quasi-particle excitation energies quantitatively in a numerical fashion, for the long-range Ising limit as well as in the intermediate regime we use perturbative continuous unitary transformations on white graphs in combination with classical Monte Carlo simulations for the graph embedding to extract high-order series expansions in the thermodynamic limit. This enables us to determine the quantum-critical breakdown of the high-field polarized phase by analyzing the gap-closing including associated critical exponents and multiplicative logarithmic corrections. In addition, for the ferromagnetic isotropic XY model we determined the critical exponents $z$ and $ν$ analytically by a bosonic quantum-field theory. △ Less

Submitted 17 November, 2020; v1 submitted 31 July, 2020; originally announced July 2020.

Comments: 11 pages, 4 figures

Journal ref: Phys. Rev. B 102, 174424 (2020)

Dynamic structure factor of Heisenberg bilayer dimer phases in the presence of quenched disorder and frustration

Authors: M. Hörmann, K. P. Schmidt

Abstract: We investigate the influence of quenched disorder on the dynamic structure factor of Heisenberg bilayers on the square, triangular, and kagome lattice in the quantum paramagnetic phase. Perturbative continuous unitary transformations and white graphs are employed to calculate the one-triplon contribution up to high orders in perturbation about the dimer limit for bimodal and continuous disorder. F… ▽ More We investigate the influence of quenched disorder on the dynamic structure factor of Heisenberg bilayers on the square, triangular, and kagome lattice in the quantum paramagnetic phase. Perturbative continuous unitary transformations and white graphs are employed to calculate the one-triplon contribution up to high orders in perturbation about the dimer limit for bimodal and continuous disorder. For the square lattice we find that the lifetime of the gap mode is increased by stronger quantum correlations while stronger disorder effects are observed for the triangular lattice due to geometric frustration. For intra-dimer disorder, in-band energy gaps are observed for both lattices which can be understood in terms of a level repulsion on dimers with low and high intra-dimer exchange that are close in energy at the momentum where the in-band gap opens. For the highly frustrated kagome lattice disorder even allows to decrease the gap energy. In addition, the localization length of the low-energy flat band is increased up to order $7$ in perturbation theory. The interplay of quenched disorder, geometric frustration, and strong correlations leads therefore to rich structures in the dynamical structure factor of two-dimensional quantum magnets. △ Less

Submitted 1 April, 2020; originally announced April 2020.

Comments: 20 pages, 21 figures

Journal ref: Phys. Rev. B 102, 094427 (2020)

The Ising model in a light-induced quantized transverse field

Authors: J. Rohn, M. Hörmann, Claudiu Genes, K. P. Schmidt

Abstract: We investigate the influence of light-matter interactions on correlated quantum matter by studying the paradigmatic Ising model subject to a quantum Rabi coupling. This type of coupling to a confined, spatially delocalized bosonic light mode, such as provided by an optical resonator, resembles a quantized transverse magnetic field of tunable strength. As a consequence, the symmetry-broken magnetic… ▽ More We investigate the influence of light-matter interactions on correlated quantum matter by studying the paradigmatic Ising model subject to a quantum Rabi coupling. This type of coupling to a confined, spatially delocalized bosonic light mode, such as provided by an optical resonator, resembles a quantized transverse magnetic field of tunable strength. As a consequence, the symmetry-broken magnetic state breaks down for strong enough light-mater interactions to a paramagnetic state. The non-local character of the bosonic mode can change the quantum phase transition in a drastic manner, which we analyze quantitatively for the simplest case of a chain geometry (Dicke-Ising chain). The results show a direct transition between a magnetically ordered phase with zero photon density and a magnetically polarized phase with lasing behaviour of the light. Our predictions are equally valid for the dual quantized Ising chain in a conventional transverse magnetic field. △ Less

Submitted 17 March, 2020; v1 submitted 12 March, 2020; originally announced March 2020.

Comments: 17 pages, 6 figures

Journal ref: Phys. Rev. Research 2, 023131 (2020)

Exact plaquette singlet phases in an orthogonal-plaquette model

Authors: C. Boos, K. P. Schmidt

Abstract: We introduce a quantum spin-1/2 model hosting two exact plaquette singlet ground states in extended parameter regimes. There is an exact phase transition between both phases, at which the system has an extensive ground-state degeneracy. Further, the model exhibits an extensive number of conserved quantities, which allow the prediction of additional phases. We exploit this feature and explore the p… ▽ More We introduce a quantum spin-1/2 model hosting two exact plaquette singlet ground states in extended parameter regimes. There is an exact phase transition between both phases, at which the system has an extensive ground-state degeneracy. Further, the model exhibits an extensive number of conserved quantities, which allow the prediction of additional phases. We exploit this feature and explore the phase diagram in detail for two specific parameter regimes. The model is based solely on Heisenberg interactions and seems sufficiently simple to be realized in quantum materials. A general scheme to determine exact singlet product states is briefly discussed. △ Less

Submitted 2 May, 2020; v1 submitted 19 December, 2019; originally announced December 2019.

Comments: 6 pages, 3 figures

Journal ref: Phys. Rev. B 101, 144428 (2020)

Quantum robustness of fracton phases

Authors: M. Mühlhauser, M. R. Walther, D. A. Reiss, K. P. Schmidt

Abstract: The quantum robustness of fracton phases is investigated by studying the influence of quantum fluctuations on the X-Cube model and Haah's code, which realize a type-I and type-II fracton phase, respectively. To this end a finite uniform magnetic field is applied to induce quantum fluctuations in the fracton phase resulting in zero-temperature phase transitions between fracton phases and polarized… ▽ More The quantum robustness of fracton phases is investigated by studying the influence of quantum fluctuations on the X-Cube model and Haah's code, which realize a type-I and type-II fracton phase, respectively. To this end a finite uniform magnetic field is applied to induce quantum fluctuations in the fracton phase resulting in zero-temperature phase transitions between fracton phases and polarized phases. Using high-order series expansions and a variational approach, all phase transitions are classified as strongly first order, which turns out to be a consequence of the (partial) immobility of fracton excitations. Indeed, single fractons as well as few-fracton composites can hardly lower their excitation energy by delocalization due to the intriguing properties of fracton phases, as demonstrated in this work explicitly in terms of fracton quasi-particles. △ Less

Submitted 14 January, 2020; v1 submitted 29 November, 2019; originally announced November 2019.

Comments: 22 pages, 30 figures

Journal ref: Phys. Rev. B 101, 054426 (2020)

Fine-Grained Tensor Network Methods

Authors: Philipp Schmoll, Saeed S. Jahromi, Max Hörmann, Matthias Mühlhauser, K. P. Schmidt, Román Orús

Abstract: We develop a strategy for tensor network algorithms that allows to deal very efficiently with lattices of high connectivity. The basic idea is to fine-grain the physical degrees of freedom, i.e., decompose them into more fundamental units which, after a suitable coarse-graining, provide the original ones. Thanks to this procedure, the original lattice with high connectivity is transformed by an is… ▽ More We develop a strategy for tensor network algorithms that allows to deal very efficiently with lattices of high connectivity. The basic idea is to fine-grain the physical degrees of freedom, i.e., decompose them into more fundamental units which, after a suitable coarse-graining, provide the original ones. Thanks to this procedure, the original lattice with high connectivity is transformed by an isometry into a simpler structure, which is easier to simulate via usual tensor network methods. In particular this enables the use of standard schemes to contract infinite 2d tensor networks - such as Corner Transfer Matrix Renormalization schemes - which are more involved on complex lattice structures. We prove the validity of our approach by numerically computing the ground-state properties of the ferromagnetic spin-1 transverse-field Ising model on the 2d triangular and 3d stacked triangular lattice, as well as of the hard-core and soft-core Bose-Hubbard models on the triangular lattice. Our results are benchmarked against those obtained with other techniques, such as perturbative continuous unitary transformations and graph projected entangled pair states, showing excellent agreement and also improved performance in several regimes. △ Less

Submitted 27 May, 2020; v1 submitted 12 November, 2019; originally announced November 2019.

Comments: 5 pages, 6 figures. Supplementary material with 8 pages, 6 figures. Revised version

Journal ref: Phys. Rev. Lett. 124, 200603 (2020)

Quantum criticality of the transverse-field Ising model with long-range interactions on triangular-lattice cylinders

Authors: J. Koziol, S. Fey, S. C. Kapfer, K. P. Schmidt

Abstract: To gain a better understanding of the interplay between frustrated long-range interactions and zero-temperature quantum fluctuations, we investigate the ground-state phase diagram of the transverse-field Ising model with algebraically-decaying long-range Ising interactions on quasi one-dimensional infinite-cylinder triangular lattices. Technically, we apply various approaches including low- and hi… ▽ More To gain a better understanding of the interplay between frustrated long-range interactions and zero-temperature quantum fluctuations, we investigate the ground-state phase diagram of the transverse-field Ising model with algebraically-decaying long-range Ising interactions on quasi one-dimensional infinite-cylinder triangular lattices. Technically, we apply various approaches including low- and high-field series expansions. For the classical long-range Ising model, we investigate cylinders with an arbitrary even circumference. We show the occurrence of gapped stripe-ordered phases emerging out of the infinitely-degenerate nearest-neighbor Ising ground-state space on the two-dimensional triangular lattice. Further, while cylinders with circumferences $6$, $10$, $14$ et cetera are always in the same stripe phase for any decay exponent of the long-range Ising interaction, the family of cylinders with circumferences $4$, $8$, $12$ et cetera displays a phase transition between two different types of stripe structures. For the full long-range transverse-field Ising model, we concentrate on cylinders with circumference four and six. The ground-state phase diagram consists of several quantum phases in both cases including an $x$-polarized phase, stripe-ordered phases, and clock-ordered phases which emerge from an order-by-disorder scenario already present in the nearest-neighbor model. In addition, the generic presence of a potential intermediate gapless phase with algebraic correlations and associated Kosterlitz-Thouless transitions is discussed for both cylinders. △ Less

Submitted 24 July, 2019; originally announced July 2019.

Comments: 12 pages, 10 figures

Journal ref: Phys. Rev. B 100, 144411 (2019)

Competition between intermediate plaquette phases in SrCu$_2$(BO$_3$)$_2$ under pressure

Authors: C. Boos, S. P. G. Crone, I. A. Niesen, P. Corboz, K. P. Schmidt, F. Mila

Abstract: Building on the growing evidence based on NMR, magnetization, neutron scattering, ESR, and specific heat that, under pressure, SrCu$_2$(BO$_3$)$_2$ has an intermediate phase between the dimer and the Néel phase, we study the competition between two candidate phases in the context of a minimal model that includes two types of intra- and inter-dimer interactions without enlarging the unit cell. We s… ▽ More Building on the growing evidence based on NMR, magnetization, neutron scattering, ESR, and specific heat that, under pressure, SrCu$_2$(BO$_3$)$_2$ has an intermediate phase between the dimer and the Néel phase, we study the competition between two candidate phases in the context of a minimal model that includes two types of intra- and inter-dimer interactions without enlarging the unit cell. We show that the empty plaquette phase of the Shastry-Sutherland model is quickly replaced by a quasi-1D full plaquette phase when intra- and/or inter-dimer couplings take different values, and that this full plaquette phase is in much better agreement with available experimental data than the empty plaquette one. △ Less

Submitted 18 June, 2019; v1 submitted 19 March, 2019; originally announced March 2019.

Comments: 19 pages

Journal ref: Phys. Rev. B 100, 140413 (2019)

Quantum robustness and phase transitions of the 3D Toric Code in a field

Authors: D. A. Reiss, K. P. Schmidt

Abstract: 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… ▽ More 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. △ Less

Submitted 27 May, 2019; v1 submitted 11 February, 2019; originally announced February 2019.

Comments: 45 pages, 14 figures, minor changes

Journal ref: SciPost Phys. 6, 078 (2019)

The Higgs Mode of Planar Coupled Spin-Ladders and its Observation in C$_9$H$_{18}$N$_2$CuBr$_4$

Authors: Tao Ying, Kai P. Schmidt, Stefan Wessel

Abstract: Polarized inelastic neutron scattering experiments recently identified the amplitude (Higgs) mode in C$_9$H$_{18}$N$_2$CuBr$_4$, a two-dimensional near-quantum-critical spin-1/2 two-leg ladder compound, which exhibits a weak easy-axis exchange anisotropy. Here, we theoretically examine the dynamic spin structure factor of such planar coupled spin-ladder systems using large-scale quantum Monte Carl… ▽ More Polarized inelastic neutron scattering experiments recently identified the amplitude (Higgs) mode in C$_9$H$_{18}$N$_2$CuBr$_4$, a two-dimensional near-quantum-critical spin-1/2 two-leg ladder compound, which exhibits a weak easy-axis exchange anisotropy. Here, we theoretically examine the dynamic spin structure factor of such planar coupled spin-ladder systems using large-scale quantum Monte Carlo simulations. This allows us to provide a quantitative account of the experimental neutron scattering data within a consistent quantum spin model. Moreover, we trance the details of the continuous evolution of the amplitude mode from a two-particle bound state of coupled ladders in the classical Ising limit all the way to the quantum spin-1/2 Heisenberg limit with fully restored SU(2) symmetry, where it gets overdamped by the two-magnon continuum in neutron scattering. △ Less

Submitted 14 December, 2018; originally announced December 2018.

Comments: 9 pages, 10 figures

Journal ref: Phys. Rev. Lett. 122, 127201 (2019)

Bilayer Kitaev models: Phase diagrams and novel phases

Authors: Urban F. P. Seifert, Julian Gritsch, Erik Wagner, Darshan G. Joshi, Wolfram Brenig, Matthias Vojta, Kai P. Schmidt

Abstract: Kitaev's honeycomb-lattice spin-$1/2$ model has become a paradigmatic example for $\mathbb{Z}_2$ quantum spin liquids, both gapped and gapless. Here we study the fate of these spin-liquid phases in differently stacked bilayer versions of the Kitaev model. Increasing the ratio between the inter-layer Heisenberg coupling $J_\perp$ and the intra-layer Kitaev couplings $K^{x,y,z}$ destroys the topolog… ▽ More Kitaev's honeycomb-lattice spin-$1/2$ model has become a paradigmatic example for $\mathbb{Z}_2$ quantum spin liquids, both gapped and gapless. Here we study the fate of these spin-liquid phases in differently stacked bilayer versions of the Kitaev model. Increasing the ratio between the inter-layer Heisenberg coupling $J_\perp$ and the intra-layer Kitaev couplings $K^{x,y,z}$ destroys the topological spin liquid in favor of a paramagnetic dimer phase. We study phase diagrams as a function of $J_\perp/K$ and Kitaev coupling anisotropies using Majorana-fermion mean-field theory, and we employ different expansion techniques in the limits of small and large $J_\perp/K$. For strongly anisotropic Kitaev couplings, we derive effective models for the different layer stackings which we use to discuss the quantum phase transition out of the Kitaev phase. We find that the phase diagrams depend sensitively on the nature of the stacking and anisotropy strength. While in some stackings and at strong anisotropies we find a single transition between the Kitaev and dimer phases, other stackings are more involved: Most importantly, we prove the existence of two novel macro-spin phases which can be understood in terms of Ising chains which can be either coupled ferromagnetically or remain degenerate, thus realizing a classical spin liquid. In addition, our results suggest the existence of a flux phase with spontaneous inter-layer coherence. △ Less

Submitted 1 October, 2018; v1 submitted 5 June, 2018; originally announced June 2018.

Comments: 23 pages, 18 figures; expanded summary of results and minor improvements

Journal ref: Phys. Rev. B 98, 155101 (2018)

Dynamic structure factor of disordered quantum spin ladders

Authors: M. Hörmann, P. Wunderlich, K. P. Schmidt

Abstract: We investigate the impact of quenched disorder on the dynamical correlation functions of two-leg quantum spin ladders. Perturbative continuous unitary transformations with the help of white graphs and bond-operator mean-field theory are used to calculate the one- and two-triplon contribution of the zero-temperature dynamical structure factor. Disorder results in huge effects on quasi-particles as… ▽ More We investigate the impact of quenched disorder on the dynamical correlation functions of two-leg quantum spin ladders. Perturbative continuous unitary transformations with the help of white graphs and bond-operator mean-field theory are used to calculate the one- and two-triplon contribution of the zero-temperature dynamical structure factor. Disorder results in huge effects on quasi-particles as well as composite bound states due to localization. This leads to intriguing quantum structures in dynamical correlation functions well observable in spectroscopic experiments. △ Less

Submitted 5 June, 2018; originally announced June 2018.

Journal ref: Phys. Rev. Lett. 121, 167201 (2018)

Tuning the $S = 1/2$ square-lattice antiferromagnet Sr$_2$Cu(Te$_{1-x}$W$_x$)O$_6$ from Néel order to quantum disorder to columnar order

Authors: O. Mustonen, S. Vasala, K. P. Schmidt, E. Sadrollahi, H. C. Walker, I. Terasaki, F. J. Litterst, E. Baggio-Saitovitch, M. Karppinen

Abstract: The spin-1/2 square-lattice Heisenberg model is predicted to have a quantum disordered ground state when magnetic frustration is maximized by competing nearest-neighbor $J_1$ and next-nearest-neighbor $J_2$ interactions ($J_2/J_1 \approx 0.5$). The double perovskites Sr$_2$CuTeO$_6$ and Sr$_2$CuWO$_6$ are isostructural spin-1/2 square-lattice antiferromagnets with Néel ($J_1$ dominates) and column… ▽ More The spin-1/2 square-lattice Heisenberg model is predicted to have a quantum disordered ground state when magnetic frustration is maximized by competing nearest-neighbor $J_1$ and next-nearest-neighbor $J_2$ interactions ($J_2/J_1 \approx 0.5$). The double perovskites Sr$_2$CuTeO$_6$ and Sr$_2$CuWO$_6$ are isostructural spin-1/2 square-lattice antiferromagnets with Néel ($J_1$ dominates) and columnar ($J_2$ dominates) magnetic order, respectively. Here we characterize the full isostructural solid solution series Sr$_2$Cu(Te$_{1-x}$W$_x$)O$_6$ ($0 \leq x \leq 1$) tunable from Néel order to quantum disorder to columnar order. A spin-liquid-like ground state was previously observed for the $x$ = 0.5 phase, but we show that the magnetic order is suppressed below 1.5 K in a much wider region of $x \approx$ 0.1-0.6. This coincides with significant $T$-linear terms in the low-temperature specific heat. However, density functional theory calculations predict most of the materials are not in the highly frustrated $J_2/J_1 \approx 0.5$ region square-lattice Heisenberg model. Thus, a combination of both magnetic frustration and quenched disorder is the likely origin of the spin-liquid-like state in $x$ = 0.5. △ Less

Submitted 26 July, 2018; v1 submitted 23 April, 2018; originally announced April 2018.

Comments: 20+5 pages, 6+4 figures. Accepted for publication in PRB

Journal ref: Phys. Rev. B 98, 064411 (2018)

Quantum criticality of two-dimensional quantum magnets with long-range interactions

Authors: S. Fey, Sebastian C. Kapfer, K. P. Schmidt

Abstract: We study the critical breakdown of two-dimensional quantum magnets in the presence of algebraically decaying long-range interactions by investigating the transverse-field Ising model on the square and triangular lattice. This is achieved technically by combining perturbative continuous unitary transformations with classical Monte Carlo simulations to extract high-order series for the one-particle… ▽ More We study the critical breakdown of two-dimensional quantum magnets in the presence of algebraically decaying long-range interactions by investigating the transverse-field Ising model on the square and triangular lattice. This is achieved technically by combining perturbative continuous unitary transformations with classical Monte Carlo simulations to extract high-order series for the one-particle excitations in the high-field quantum paramagnet. We find that the unfrustrated systems change from mean-field to nearest-neighbor universality with continuously varying critical exponents, while the system remains in the universality class of the nearest-neighbor model in the frustrated cases independent of the long-range nature of the interaction. △ Less

Submitted 19 February, 2018; originally announced February 2018.

Comments: 13 pages, 7 figures

Journal ref: Phys. Rev. Lett. 122, 017203 (2019)

Spin-liquid-like state in a spin-1/2 square-lattice antiferromagnet perovskite induced by $d^{10}$-$d^0$ cation mixing

Authors: O. Mustonen, S. Vasala, E. Sadrollahi, K. P. Schmidt, C. Baines, H. C. Walker, I. Terasaki, F. J. Litterst, E. Baggio-Saitovitch, M. Karppinen

Abstract: A quantum spin liquid state has long been predicted to arise in spin-1/2 Heisenberg square-lattice antiferromagnets at the boundary region between Néel (nearest-neighbor interaction dominates) and columnar (next-nearest-neighbor dominates) antiferromagnetic order. However, there are no known compounds in this region. Here we use $d^{10}$-$d^0$ cation mixing to tune the magnetic interactions on the… ▽ More A quantum spin liquid state has long been predicted to arise in spin-1/2 Heisenberg square-lattice antiferromagnets at the boundary region between Néel (nearest-neighbor interaction dominates) and columnar (next-nearest-neighbor dominates) antiferromagnetic order. However, there are no known compounds in this region. Here we use $d^{10}$-$d^0$ cation mixing to tune the magnetic interactions on the square lattice while simultaneously introducing disorder. We find spin-liquid-like behavior in the double perovskite Sr$_2$Cu(Te$_{0.5}$W$_{0.5}$)O$_6$, where the isostructural end phases Sr$_2$CuTeO$_6$ and Sr$_2$CuWO$_6$ are Néel and columnar type antiferromagnets, respectively. We show that magnetism in Sr$_2$Cu(Te$_{0.5}$W$_{0.5}$)O$_6$ is entirely dynamic down to 19 mK. Additionally, we observe at low temperatures for Sr$_2$Cu(Te$_{0.5}$W$_{0.5}$)O$_6$, similar to several spin liquid candidates, a plateau in muon spin relaxation rate and a strong $T$-linear dependence in specific heat. Our observations for Sr$_2$Cu(Te$_{0.5}$W$_{0.5}$)O$_6$ highlight the role of disorder in addition to magnetic frustration in spin liquid physics. △ Less

Submitted 13 February, 2018; originally announced February 2018.

Comments: 19+4 pages, 4+3 figures. Accepted for publication in Nature Communications

Journal ref: Nat. Commun. 9 (2018) 1085

Time-reversal symmetry breaking Abelian chiral spin liquid in Mott phases of three-component fermions on the triangular lattice

Authors: C. Boos, C. J. Ganahl, M. Lajkó, P. Nataf, A. M. Läuchli, K. Penc, K. P. Schmidt, F. Mila

Abstract: We provide numerical evidence in favor of spontaneous chiral symmetry breaking and the concomitant appearance of an Abelian chiral spin liquid for three-component fermions on the triangular lattice described by an SU(3) symmetric Hubbard model with hopping amplitude $-t$ ($t>0$) and on-site interaction $U$. This chiral phase is stabilized in the Mott phase with one particle per site in the presenc… ▽ More We provide numerical evidence in favor of spontaneous chiral symmetry breaking and the concomitant appearance of an Abelian chiral spin liquid for three-component fermions on the triangular lattice described by an SU(3) symmetric Hubbard model with hopping amplitude $-t$ ($t>0$) and on-site interaction $U$. This chiral phase is stabilized in the Mott phase with one particle per site in the presence of a uniform $π$-flux per plaquette, and in the Mott phase with two particles per site without any flux. Our approach relies on effective spin models derived in the strong-coupling limit in powers of $t/U$ for general SU$(N)$ and arbitrary uniform charge flux per plaquette, which are subsequently studied using exact diagonalizations and variational Monte Carlo simulations for $N=3$, as well as exact diagonalizations of the SU($3$) Hubbard model on small clusters. Up to third order in $t/U$, and for the time-reversal symmetric cases (flux $0$ or $π$), the low-energy description is given by the $J$-$K$ model with Heisenberg coupling $J$ and real ring exchange $K$. The phase diagram in the full $J$-$K$ parameter range contains, apart from three already known, magnetically long-range ordered phases, two previously unreported phases: i) a lattice nematic phase breaking the lattice rotation symmetry and ii) a spontaneous time-reversal and parity symmetry breaking Abelian chiral spin liquid. For the Hubbard model, an investigation that includes higher-order itinerancy effects supports the presence of a phase transition inside the insulating region, occurring at $(t/U)_{\rm c}\approx 0.07$ [$(U/t)_{\rm c} \approx 13$] between the three-sublattice magnetically ordered phase at small $t/U$ and this Abelian chiral spin liquid. △ Less

Submitted 2 May, 2020; v1 submitted 9 February, 2018; originally announced February 2018.

Comments: 21 pages, 23 figures

Journal ref: Phys. Rev. Research 2, 023098 (2020)

Quantum phase transitions to topological Haldane phases in spin-one chains studied by linked-cluster expansions

Authors: P. Adelhardt, J. Gritsch, M. Hille, D. A. Reiss, K. P. Schmidt

Abstract: We use linked-cluster expansions to analyze the quantum phase transitions between symmetry unbroken trivial and topological Haldane phases in two different spin-one chains. The first model is the spin-one Heisenberg chain in the presence of a single-ion anisotropy while the second one is the dimerized spin-one Heisenberg chain. For both models we determine the ground-state energy and the one-parti… ▽ More We use linked-cluster expansions to analyze the quantum phase transitions between symmetry unbroken trivial and topological Haldane phases in two different spin-one chains. The first model is the spin-one Heisenberg chain in the presence of a single-ion anisotropy while the second one is the dimerized spin-one Heisenberg chain. For both models we determine the ground-state energy and the one-particle gap inside the non-topological phase as a high-order series using perturbative continuous unitary transformations. Extrapolations of the gap series are applied to locate the quantum critical point and to extract the associated critical exponent. We find that this approach works unsatisfactory for the anisotropic chain, since the quality of the extrapolation appears insufficient due to the large correlation length exponent. In contrast, extrapolation schemes display very good convergence for the gap closing in the case of the dimerized spin-one Heisenberg chain. △ Less

Submitted 9 October, 2017; originally announced October 2017.

Comments: 7 pages, 4 figures

Journal ref: Phys. Rev. B 96, 235123 (2017)