Publications

Order by: publication date number of citations

• Pauli crystal melting in shaken optical traps

  Jiabing Xiang, Paolo Molignini, Miriam Büttner, Axel U. J. Lode

  SciPost Phys. 14, 003 (2023) · published 13 January 2023 | Toggle abstract · pdf

  Pauli crystals are ordered geometric structures that emerge in trapped noninteracting fermionic systems due to their underlying Pauli repulsion. The deformation of Pauli crystals - often called melting - has been recently observed in experiments, but the mechanism that leads to it remains unclear. We address this question by studying the melting dynamics of $N=6$ fermions as a function of periodic driving and experimental imperfections in the trap (anisotropy and anharmonicity) by employing a combination of numerical simulations and Floquet theory. Surprisingly, we reveal that the melting of Pauli crystals is not simply a direct consequence of an increase in system energy, but is instead related to the trap geometry and the population of the Floquet modes. We show that the melting is absent in traps without imperfections and triggered only by a sufficiently large shaking amplitude in traps with imperfections.

  1 citation

• Topological defects in a double-mirror quadrupole insulator displace diverging charge

  Isidora Araya Day, Anton R. Akhmerov, Dániel Varjas

  SciPost Phys. Core 5, 053 (2022) · published 15 December 2022 | Toggle abstract · pdf

  We show that topological defects in quadrupole insulators do not host quantized fractional charges, contrary to what their Wannier representation indicates. In particular, we test the charge quantization hypothesis based on the Wannier representation of a disclination and a parametric defect. Since disclinations necessarily strain the lattice and parametric defects require closed curves in parameter space, both defects break four-fold rotation symmetry, even away from their origin. The Wannier representation of the defects is thus determined by local reflection symmetries. Contrary to the hypothesis, we find that the local charge density decays as $\sim 1/r^2$ with distance, leading to a diverging defect charge. Because topological defects are incompatible with four-fold rotation symmetry, we conclude that defect charge quantization is protected by sublattice symmetry, and not higher order topology.

  1 citation

• Floquet engineering of non-equilibrium superradiance

  Lukas Broers, Ludwig Mathey

  SciPost Phys. 14, 018 (2023) · published 13 February 2023 | Toggle abstract · pdf

  We demonstrate the emergence of a non-equilibrium superradiant phase in the dissipative Rabi-Dicke model. This phase is characterized by a photonic steady state that oscillates with a frequency close to the cavity frequency, in contrast to the constant photonic steady state of the equilibrium superradiant phase in the Dicke model. We relate this superradiant phase to the population inversion of Floquet states by introducing a Schwinger representation of the driven two-level systems in the cavity. This inversion is depleted near Floquet energies that are resonant with the cavity frequency to sustain a coherent light-field. In particular, our model applies to solids within a two-band approximation, in which the electrons act as Schwinger fermions. We propose to use this Floquet-assisted superradiant phase to obtain controllable optical gain for a laser-like operation.

  1 citation

• Automatic transformation of irreducible representations for efficient contraction of tensors with cyclic group symmetry

  Yang Gao, Phillip Helms, Garnet Kin-Lic Chan, Edgar Solomonik

  SciPost Phys. Codebases 10 (2023) · published 24 February 2023 | Toggle abstract · pdf

  Tensor contractions are ubiquitous in computational chemistry and physics, where tensors generally represent states or operators and contractions express the algebra of these quantities. In this context, the states and operators often preserve physical conservation laws, which are manifested as group symmetries in the tensors. These group symmetries imply that each tensor has block sparsity and can be stored in a reduced form. For nontrivial contractions, the memory footprint and cost are lowered, respectively, by a linear and a quadratic factor in the number of symmetry sectors. State-of-the-art tensor contraction software libraries exploit this opportunity by iterating over blocks or using general block-sparse tensor representations. Both approaches entail overhead in performance and code complexity. With intuition aided by tensor diagrams, we present a technique, irreducible representation alignment, which enables efficient handling of Abelian group symmetries via only dense tensors, by using contraction-specific reduced forms. This technique yields a general algorithm for arbitrary group symmetric contractions, which we implement in Python and apply to a variety of representative contractions from quantum chemistry and tensor network methods. As a consequence of relying on only dense tensor contractions, we can easily make use of efficient batched matrix multiplication via Intel's MKL and distributed tensor contraction via the Cyclops library, achieving good efficiency and parallel scalability on up to 4096 Knights Landing cores of a supercomputer.

  1 citation

• Codebase release 1.3 for SymTensor

  Yang Gao, Phillip Helms, Garnet Kin-Lic Chan, Edgar Solomonik

  SciPost Phys. Codebases 10-r1.3 (2023) · published 24 February 2023 | Toggle abstract · src

  Tensor contractions are ubiquitous in computational chemistry and physics, where tensors generally represent states or operators and contractions express the algebra of these quantities. In this context, the states and operators often preserve physical conservation laws, which are manifested as group symmetries in the tensors. These group symmetries imply that each tensor has block sparsity and can be stored in a reduced form. For nontrivial contractions, the memory footprint and cost are lowered, respectively, by a linear and a quadratic factor in the number of symmetry sectors. State-of-the-art tensor contraction software libraries exploit this opportunity by iterating over blocks or using general block-sparse tensor representations. Both approaches entail overhead in performance and code complexity. With intuition aided by tensor diagrams, we present a technique, irreducible representation alignment, which enables efficient handling of Abelian group symmetries via only dense tensors, by using contraction-specific reduced forms. This technique yields a general algorithm for arbitrary group symmetric contractions, which we implement in Python and apply to a variety of representative contractions from quantum chemistry and tensor network methods. As a consequence of relying on only dense tensor contractions, we can easily make use of efficient batched matrix multiplication via Intel's MKL and distributed tensor contraction via the Cyclops library, achieving good efficiency and parallel scalability on up to 4096 Knights Landing cores of a supercomputer.

  1 citation

• Adiabatic and irreversible classical discrete time crystals

  Adrian Ernst, Anna M. E. B. Rossi, Thomas M. Fischer

  SciPost Phys. 13, 091 (2022) · published 11 October 2022 | Toggle abstract · pdf

  We simulate the dynamics of paramagnetic colloidal particles that are placed above a magnetic hexagonal pattern and exposed to an external field periodically changing its direction along a control loop. The conformation of three colloidal particles above one unit cell adiabatically responds with half the frequency of the external field creating a time crystal at arbitrary low frequency. The adiabatic time crystal occurs because of the non-trivial topology of the stationary manifold. When coupling colloidal particles in different unit cells, many body effects cause the formation of topologically isolated time crystals and dynamical phase transitions between different adiabatic reversible and non-adiabatic irreversible space-time-crystallographic arrangements.

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

  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)]

• Efficient scheme for time-dependent thermal pure quantum state: Application to the Kitaev model with armchair edges

  Hirokazu Taguchi, Akihisa Koga, Yuta Murakami

  SciPost Phys. Proc. 11, 003 (2023) · published 5 June 2023 | Toggle abstract · pdf

  We consider the time-dependent thermal pure quantum state method and introduce the efficient scheme to evaluate the change in physical quantities induced by the time-dependent perturbations, which has been proposed in our previous paper [H. Taguchi et al., Phys. Rev. B 105, 125137 (2022)]. Here, we treat the Kitaev model to consider the Majorana-mediated spin transport, as an example. We demonstrate how efficient our scheme is to evaluate spin oscillations induced by the magnetic field pulse.

• Real-space spectral simulation of quantum spin models: Application to generalized Kitaev models

  Francisco M. O. Brito, Aires Ferreira

  SciPost Phys. Core 7, 006 (2024) · published 21 February 2024 | Toggle abstract · pdf

  The proliferation of quantum fluctuations and long-range entanglement presents an outstanding challenge for the numerical simulation of interacting spin systems with exotic ground states. Here, we present a toolset of Chebyshev polynomial-based iterative methods that provides a unified framework to study the thermodynamical properties, critical behavior and dynamics of frustrated quantum spin models with controlled accuracy. Similar to previous applications of the Chebyshev spectral methods to condensed matter systems, the algorithmic complexity scales linearly with the Hilbert space dimension and the Chebyshev truncation order. Using this approach, we study two paradigmatic quantum spin models on the honeycomb lattice: the Kitaev-Heisenberg (K-H) and the Kitaev-Ising (K-I) models. We start by applying the Chebyshev toolset to compute nearest-neighbor spin correlations, specific heat and entropy of the K-H model on a 24-spin cluster. Our results are benchmarked against exact diagonalization and a popular iterative method based on thermal pure quantum states. The transitions between a variety of magnetic phases, namely ferromagnetic, Néel, zigzag and stripy antiferromagnetic and quantum spin liquid phases are obtained accurately and efficiently. We also determine the temperature dependence of the spin correlations, over more than three decades in temperature, by means of a finite temperature Chebyshev polynomial method introduced here. Finally, we report novel dynamical signatures of the quantum phase transitions in the K-I model. Our findings suggest that the efficiency, versatility and low-temperature stability of the Chebyshev framework developed here could pave the way for previously unattainable studies of quantum spin models in two dimensions.

• Quantum spin nematic liquid in the low-dimensional anisotropic magnets -S=1/2 delta spin chain with the anisotropic ferromagnetic interaction in magnetic field-

  Toru Sakai, Rito Furuchi, Hiroki Nakano, Kiyomi Okamoto

  SciPost Phys. Proc. 11, 001 (2023) · published 5 June 2023 | Toggle abstract · pdf

  The magnetization process of the $S=1/2$ delta chain with the anisotropic ferromagnetic interaction is investigated using the numerical diagonalization of finite-size clusters. It is found that the spin nematic liquid phase appears in higher magnetization region, as well as the SDW liquid one in lower region.