Publications

Order by: publication date number of citations

• Two $T$-linear scattering-rate regimes in the triangular lattice Hubbard model

  Jérôme Fournier, Pierre-Olivier Downey, Charles-David Hébert, Maxime Charlebois, André-Marie Tremblay

  SciPost Phys. 17, 072 (2024) · published 4 September 2024 | Toggle abstract · pdf

  In recent years, the $T$-linear scattering rate found at low temperatures, defining the strange metal phase of cuprates, has been a subject of interest. Since a wide range of materials have a scattering rate that obeys the equation $ \hbar / \tau ≈ k_B T$, the idea of a universal Planckian limit on the scattering rate has been proposed. However, there is no consensus on proposed theories yet. In this work, we present our results for the $T$-linear scattering rate in the triangular lattice Hubbard model obtained using the dynamical cluster approximation. We find two regions with $T$-linear scattering rate in the $T$-$p$ phase diagram: one emerges from the pseudogap to correlated Fermi liquid phase transition at low doping, whereas the other is solely caused by large interaction strength at large doping.

  2 citations

• Comparative microscopic study of entropies and their production

  Philipp Strasberg, Joseph Schindler

  SciPost Phys. 17, 143 (2024) · published 25 November 2024 | Toggle abstract · pdf

  We study the time evolution of eleven microscopic entropy definitions (of Boltzmann-surface, Gibbs-volume, canonical, coarse-grained-observational, entanglement and diagonal type) and three microscopic temperature definitions (based on Boltzmann, Gibbs or canonical entropy). This is done for the archetypal nonequilibrium setup of two systems exchanging energy, modeled here with random matrix theory, based on numerical integration of the Schrödinger equation. We consider three types of pure initial states (local energy eigenstates, decorrelated and entangled microcanonical states) and three classes of systems: (A) two normal systems, (B) a normal and a negative temperature system and (C) a normal and a negative heat capacity system. We find: (1) All types of initial states give rise to the same macroscopic dynamics. (2) Entanglement and diagonal entropy sensitively depend on the microstate, in contrast to all other entropies. (3) For class B and C, Gibbs-volume entropies can violate the second law and the associated temperature becomes meaningless. (4) For class C, Boltzmann-surface entropies can violate the second law and the associated temperature becomes meaningless. (5) Canonical entropy has a tendency to remain almost constant. (6) For a Haar random initial state, entanglement or diagonal entropy behave similar or identical to coarse-grained-observational entropy.

  2 citations

• $\mathbb{Z}_4$ transitions in quantum loop models on a zig-zag ladder

  Bowy M. La Rivière, Natalia Chepiga

  SciPost Phys. 17, 144 (2024) · published 26 November 2024 | Toggle abstract · pdf

  We study the nature of quantum phase transitions out of $\mathbb{Z}_4$ ordered phases in quantum loop models on a zig-zag ladder. We report very rich critical behavior that includes a pair of Ising transitions, a multi-critical Ashkin-Teller point and a remarkably extended interval of a chiral transition. Although plaquette states turn out to be essential to realize chiral transitions, we demonstrate that critical regimes can be manipulated by deforming the model as to increase the presence of leg-dimerized states. This can be done to the point where the chiral transition turns into first order, we argue that this is associated with the emergence of a critical end point.

  2 citations

• Tunable second harmonic generation in 2D materials: Comparison of different strategies

  Simone Grillo, Elena Cannuccia, Maurizia Palummo, Olivia Pulci, Claudio Attaccalite

  SciPost Phys. Core 7, 081 (2024) · published 4 December 2024 | Toggle abstract · pdf

  Nonlinear optical frequency conversion, where optical fields interact with a nonlinear medium to generate new frequencies, is a key phenomenon in modern photonic systems. However, a major challenge with these techniques lies in the difficulty of tuning the nonlinear electrical susceptibilities that drive such effects in a given material. As a result, dynamic control of optical nonlinearities has remained largely confined to research laboratories, limiting its practical use as a spectroscopic tool. In this work, we aim to advance the development of devices with tunable nonlinear responses by exploring two potential mechanisms for electrically manipulating second-order optical nonlinearity in two-dimensional materials. Specifically, we consider two configurations: in the first, the material does not inherently exhibit second-harmonic generation (SHG), but this response is induced by an external field; in the second, an external field induces doping in a material that already exhibits SHG, altering the intensity of the nonlinear signal. In this work, we have studied these two configurations using a real-time ab-initio approach under an out-of-plane external field and including the effects of doping-induced variations in the screened electron-electron interaction. We then discuss the limitations of current computational methods and compare our results with experimental measurements.

  2 citations

• Understanding and utilizing the inner bonds of process tensors

  Moritz Cygorek, Erik M. Gauger

  SciPost Phys. 18, 024 (2025) · published 20 January 2025 | Toggle abstract · pdf

  Process tensor matrix product operators (PT-MPOs) enable numerically exact simulations for an unprecedentedly broad range of open quantum systems. By representing environment influences in MPO form, they can be efficiently compressed using established algorithms. The dimensions of inner bonds of the compressed PT-MPO may be viewed as an indicator of the complexity of the environment. Here, we show that the inner bonds themselves, not only their dimensions, have a concrete physical meaning: They represent the subspace of the full environment Liouville space which hosts environment excitations that may influence the subsequent open quantum systems dynamics the most. This connection can be expressed in terms of lossy linear transformations, whose pseudoinverses facilitate the extraction of environment observables. We demonstrate this by extracting the environment spin of a central spin problem, the current through a quantum system coupled to two leads, the number of photons emitted from quantum emitters into a structured environment, and the distribution of the total absorbed energy in a driven non-Markovian quantum system into system, environment, and interaction energy terms. Numerical tests further indicate that different PT-MPO algorithms compress environments to similar subspaces. Thus, the physical interpretation of inner bonds of PT-MPOs both provides a conceptional understanding and it enables new practical applications.

  2 citations

• Insulator phases of Bose-Fermi mixtures induced by intraspecies next-neighbor interactions

  Felipe Gómez-Lozada, Roberto Franco, Jereson Silva-Valencia

  SciPost Phys. Core 8, 007 (2025) · published 21 January 2025 | Toggle abstract · pdf

  We study a one-dimensional mixture of two-color fermions and scalar bosons at the hardcore limit, focusing on the effect that the intraspecies next-neighbor interactions have on the zero-temperature ground state of the system for different fillings of each carrier. Exploring the problem's parameters, we observed that the nearest-neighbor interaction could favor or harm the well-known mixed Mott and spin-selective Mott insulators. We also found the emergence of three unusual insulating states with charge density wave (CDW) structures in which the orders of the carriers are out of phase between each other. For instance, the immiscible CDW appears only at half-filling bosonic density, whereas the mixed CDW state is characterized by equal densities of bosons and fermions. Finally, the spin-selective CDW couples the bosons and only one kind of fermions. Appropriate order parameters were proposed for each phase to obtain the critical parameters for the corresponding superfluid-insulator transition. Our results can inspire or contribute to understanding experiments in cold-atom setups with long-range interactions or recent reports involving quasiparticles in semiconductor heterostructures.

  2 citations

• DanceQ: High-performance library for number conserving bases

  Robin Schäfer, David J. Luitz

  SciPost Phys. Codebases 48 (2025) · published 5 February 2025 | Toggle abstract · pdf

  The complexity of quantum many-body problems scales exponentially with the size of the system, rendering any finite-size scaling analysis a formidable challenge. This is particularly true for methods based on the full representation of the wave function, where one simply accepts the enormous Hilbert space dimensions and performs linear algebra operations, e.g., for finding the ground state of the Hamiltonian. If the system satisfies an underlying symmetry where an operator with a degenerate spectrum commutes with the Hamiltonian, it can be block-diagonalized, thus reducing the complexity at the expense of additional bookkeeping. At the most basic level required for Krylov space techniques (like the Lanczos algorithm), it is necessary to implement a matrix-vector product of a block of the Hamiltonian with arbitrary block-wavefunctions, potentially without holding the Hamiltonian block in memory. An efficient implementation of this operation requires the calculation of the position of an arbitrary basis vector in the canonical ordering of the basis of the block. We present here an elegant and powerful, multi-dimensional approach to this problem for the $U(1)$ symmetry appearing in problems with particle number conservation. Our divide-and-conquer algorithm uses multiple subsystems and hence generalizes previous approaches to make them scalable. In addition to the theoretical presentation of our algorithm, we provide DanceQ, a flexible and modern – header only – C++20 implementation to manipulate, enumerate, and map to its index any basis state in a given particle number sector as open source software under https://DanceQ.gitlab.io/danceq.

  2 citations

• Codebase release 1.0 for DanceQ

  Robin Schäfer, David J. Luitz

  SciPost Phys. Codebases 48-r1.0 (2025) · published 5 February 2025 | Toggle abstract · src

  The complexity of quantum many-body problems scales exponentially with the size of the system, rendering any finite-size scaling analysis a formidable challenge. This is particularly true for methods based on the full representation of the wave function, where one simply accepts the enormous Hilbert space dimensions and performs linear algebra operations, e.g., for finding the ground state of the Hamiltonian. If the system satisfies an underlying symmetry where an operator with a degenerate spectrum commutes with the Hamiltonian, it can be block-diagonalized, thus reducing the complexity at the expense of additional bookkeeping. At the most basic level required for Krylov space techniques (like the Lanczos algorithm), it is necessary to implement a matrix-vector product of a block of the Hamiltonian with arbitrary block-wavefunctions, potentially without holding the Hamiltonian block in memory. An efficient implementation of this operation requires the calculation of the position of an arbitrary basis vector in the canonical ordering of the basis of the block. We present here an elegant and powerful, multi-dimensional approach to this problem for the $U(1)$ symmetry appearing in problems with particle number conservation. Our divide-and-conquer algorithm uses multiple subsystems and hence generalizes previous approaches to make them scalable. In addition to the theoretical presentation of our algorithm, we provide DanceQ, a flexible and modern – header only – C++20 implementation to manipulate, enumerate, and map to its index any basis state in a given particle number sector as open source software under https://DanceQ.gitlab.io/danceq.

  2 citations

• Carlo.jl: A general framework for Monte Carlo simulations in Julia

  Lukas Weber

  SciPost Phys. Codebases 49 (2025) · published 7 February 2025 | Toggle abstract · pdf

  Carlo.jl is a Monte Carlo simulation framework written in Julia. It provides MPI-parallel scheduling (including parallel tempering), organized storage of input, checkpoint, and output files, as well as statistical postprocessing. With a minimalist design, it aims to aid the development of high-quality Monte Carlo codes, especially for demanding applications in condensed matter and statistical physics. This hands-on user guide shows how to implement a simple code with Carlo.jl, how to run an existing stochastic series expansion code based on Carlo.jl, and provides benchmarks to show its efficacy.

  2 citations

• Boosting determinant quantum Monte Carlo with submatrix updates: Unveiling the phase diagram of the 3D Hubbard model

  Fanjie Sun, Xiao Yan Xu

  SciPost Phys. 18, 055 (2025) · published 17 February 2025 | Toggle abstract · pdf

  Determinant Quantum Monte Carlo (DQMC) provides numerically exact solutions for strongly correlated fermionic systems but faces significant computational challenges with increasing system size. While submatrix updates were originally developed for Hirsch-Fye QMC with onsite interactions at finite temperatures, their comprehensive application in DQMC has remained unexplored despite noted algorithmic similarities. We present the first comprehensive application of submatrix updates in DQMC, significantly extending beyond the original scope by enabling simulations with extended interactions and at zero temperature. Building upon conventional fast updates and delay updates, our generalized implementation achieves an order-of-magnitude improvement in computational efficiency, enabling simulations of the half-filled Hubbard model on lattices up to 8,000 sites - a scale previously challenging with standard DQMC implementations. This enhanced computational capability allows us to accurately determine the finite-temperature phase diagram of the 3D Hubbard model at half-filling. Our findings not only shed light on the phase transitions within these complex systems but also pave the way for more effective simulations of strongly correlated electrons, potentially guiding experimental efforts in cold atom simulations of the 3D Hubbard model.

  2 citations