Publications

Order by: publication date number of citations

• The asymmetric Fermi surface of Bi2201

  Steef Smit, Kourosh L. Shirkoohi, Saumya Mukherjee, Sergio Barquero Pierantoni, Lewis Bawden, Erik van Heumen, Arnaud P. N. Tchiomo, Jans Henke, Jasper van Wezel, Yingkai Huang, Takeshi Kondo, Tsunehiro Takeuchi, Timur K. Kim, Cephise Cacho, Marta Zonno, Sergey Gorovikov, Stephen B. Dugdale, Jorge I. Facio, Mariia Roslova, Laura Folkers, Anna Isaeva, Nigel E. Hussey, Mark S. Golden

  SciPost Phys. 18, 191 (2025) · published 17 June 2025 | Toggle abstract · pdf

  High-resolution angle-resolved photoemission spectroscopy (ARPES) performed on the single-layered cuprate (Pb$_{1-y}$,Bi$_y$)$_2$Sr$_{2-x}$La$_x$CuO$_{6+\delta}$ (Bi2201) reveals a 6-10% difference in the nodal $k_F$ vectors along the $\Gamma$Y and $\Gamma$X directions. This asymmetry is notably larger than the 2% orthorhombic distortion in the CuO$_2$ plane lattice constants determined using X-ray crystallography from the same samples. First principles calculations indicate that crystal-field splitting of the bands lies at the root of the $k_{F}$ asymmetry. Concomitantly, the nodal Fermi velocities for the $\Gamma$Y quadrant exceed those for $\Gamma$X by 4%. Momentum distribution curve widths for the two nodal dispersions are also anisotropic, showing identical energy dependencies, bar a scaling factor of $\sim$ 1.17$± 0.05$ between $\Gamma$Y and $\Gamma$X. Consequently, the imaginary part of the self-energy is found to be 10-20% greater along $\Gamma$Y than $\Gamma$X. These results emphasize the need to account for Fermi surface asymmetry in the analysis of ARPES data on Bi-based cuprate high temperature superconductors such as Bi2201. To illustrate this point, an orthorhombic tight-binding model (with twofold in-plane symmetry) was used to fit ARPES Fermi surface maps spanning all four quadrants of the Brillouin zone, and the ARPES-derived hole-doping (Luttinger count) was extracted. Comparison of the Luttinger count with one assuming four-fold in-plane symmetry strongly suggests the marked spread in previously-reported Fermi surface areas from ARPES on Bi2201 results from the differences in $k_F$ along $\Gamma$Y and $\Gamma$X. Using this analysis, a new, linear relationship emerges between the hole-doping derived from ARPES ($p_{\text{ARPES}}$) and that derived using the Presland ($p_{\text{Presland}}$) relation such that $p_{\text{ARPES}} = p_{\text{Presland}}+0.11$. The implications for this difference between the ARPES- and Presland-derived estimates for $p$ are discussed and possible future directions to elucidate the origin of this discrepancy are presented.

• Simulation of the 1d XY model on a quantum computer

  Marc Farreras, Alba Cervera-Lierta

  SciPost Phys. Lect. Notes 95 (2025) · published 11 June 2025 | Toggle abstract · pdf

  The field of quantum computing has grown fast in recent years, both in theoretical advancements and the practical construction of quantum computers. These computers were initially proposed, among other reasons, to efficiently simulate and comprehend the complexities of quantum physics. In this paper, we present a comprehensive scheme for the exact simulation of the 1-D XY model on a quantum computer. We successfully diagonalize the proposed Hamiltonian, enabling access to the complete energy spectrum. Furthermore, we propose a novel approach to design a quantum circuit to perform exact time evolution. Among all the possibilities this opens, we compute the ground and excited state energies for the symmetric XY model with spin chains of $n=4$ and $n=8$ spins. Further, we calculate the expected value of transverse magnetization for the ground state in the transverse Ising model. Both studies allow the observation of a quantum phase transition from an antiferromagnetic to a paramagnetic state. Additionally, we have simulated the time evolution of the state all spins up in the transverse Ising model. The scalability and high performance of our algorithm make it an ideal candidate for benchmarking purposes, while also laying the foundation for simulating other integrable models on quantum computers.

  1 citation

• Geometric constructions of generalized dual-unitary circuits from biunitarity

  Michael A. Rampp, Suhail A. Rather, Pieter W. Claeys

  SciPost Phys. 18, 182 (2025) · published 10 June 2025 | Toggle abstract · pdf

  We present a general framework for constructing solvable lattice models of chaotic many-body quantum dynamics with multiple unitary directions using biunitary connections. We show that a network of biunitary connections on the kagome lattice naturally defines a multi-unitary circuit, where three 'arrows of time' directly reflect the lattice symmetry. These models unify various constructions of hierarchical dual-unitary and triunitary gates and present new families of models with solvable correlations and entanglement dynamics. Using multilayer constructions of biunitary connections, we additionally introduce multilayer circuits with monoclinic symmetry and higher level hierarchical dual-unitary solvability and discuss their (non-)ergodicity. Our work demonstrates how different classes of solvable models can be understood as arising from different geometric structures in spacetime.

  1 citation

• Effect of the Coulomb repulsion and oxygen level on charge distribution and superconductivity in the Emery model for cuprates superconductors

  Louis-Bernard St-Cyr, David Sénéchal

  SciPost Phys. Core 8, 043 (2025) · published 5 June 2025 | Toggle abstract · pdf

  The Emery model (aka the three-band Hubbard model) offers a simplified description of the copper-oxide planes that form the building blocks of high-temperature superconductors. By contrast with the even simpler one-band Hubbard model, it differentiates between copper and oxygen orbitals and thus between oxygen occupation ($n_p$) and copper occupation ($n_d$). Here we demonstrate, using cluster dynamical mean field theory, how the two occupations are related to the on-site Coulomb repulsion $U$ on the copper orbital and to the energy difference $\epsilon_p$ between oxygen and copper orbitals. Since the occupations ($n_p$ and $n_d$) have been estimated from NMR for a few materials (LCO, YBCO and NCCO), this allows us to estimate the value of $U-\epsilon_p$ for these materials, within this model. We compute the density of states for these and the effect of $(U,\epsilon_p)$ on the $n_d$-$n_p$ curve, superconductivity, and antiferromagnetism.

  1 citation

• Stabilizer disentangling of conformal field theories

  Martina Frau, Poetri Sonya Tarabunga, Mario Collura, Emanuele Tirrito, Marcello Dalmonte

  SciPost Phys. 18, 165 (2025) · published 22 May 2025 | Toggle abstract · pdf

  Understanding how entanglement can be reduced through simple operations is crucial for both classical and quantum algorithms. We investigate the entanglement properties of lattice models hosting conformal field theories cooled via local Clifford operations, a procedure we refer to as stabilizer disentangling. We uncover two distinct regimes: a constant gain regime, where disentangling is volume-independent, and a log-gain regime, where disentanglement increases with volume, characterized by a reduced effective central charge. In both cases, disentangling efficiency correlates with the target state magic, with larger magic leading to more effective cooling. The dichotomy between the two cases stems from mutual stabilizer Rényi entropy, which influences the entanglement cooling process. We provide an analytical understanding of such effect in the context of cluster Ising models, that feature disentangling global Clifford operations. Our findings indicate that matrix product states possess subclasses based on the relationship between entanglement and magic, and clarifying the potential of new classes of variational states embedding Clifford dynamics within matrix product states.

  5 citations

• Three-dimensional quantum Hall effect in topological amorphous metals

  Jiong-Hao Wang, Yong Xu

  SciPost Phys. 18, 146 (2025) · published 2 May 2025 | Toggle abstract · pdf

  Weyl semimetals have been theoretically predicted to become topological metals with anomalous Hall conductivity in amorphous systems. However, measuring the anomalous Hall conductivity in realistic materials, particularly those with multiple pairs of Weyl points, is a significant challenge. If a system respects time-reversal symmetry, then the anomalous Hall conductivity even vanishes. As such, it remains an open question how to probe the Weyl band like topology in amorphous materials. Here, we theoretically demonstrate that, under magnetic fields, a topological metal slab in amorphous systems exhibits three-dimensional quantum Hall effect, even in time-reversal invariant systems, thereby providing a feasible approach to exploring Weyl band like topology in amorphous materials. We unveil the topological origin of the quantized Hall conductance by calculating the Bott index. The index is carried by broadened Landau levels with bulk states spatially localized except at critical transition energies. The topological property also results in edge states localized at distinct hinges on two opposite surfaces.

  1 citation

• Single-particle spectra and magnetic susceptibility in the Emery model: A dynamical mean-field perspective

  Yi-Ting Tseng, Mário O. Malcolms, Henri Menke, Marcel Klett, Thomas Schäfer, Philipp Hansmann

  SciPost Phys. 18, 145 (2025) · published 1 May 2025 | Toggle abstract · pdf

  We investigate dynamical mean-field calculations of the three-band Emery model at the one- and two-particle level for material-realistic parameters of high-$T_c$ superconductors. Our study shows that even within dynamical mean-field theory, which accounts solely for temporal fluctuations, the intrinsic multi-orbital nature of the Emery model introduces effective non-local correlations. These correlations lead to a non-Curie-like temperature dependence of the magnetic susceptibility, consistent with nuclear magnetic resonance experiments in the pseudogap regime. By analyzing the temperature dependence of the uniform static spin susceptibility obtained by single-site and cluster dynamical mean-field theory, we find indications of emerging oxygen-copper singlet fluctuations, explicitly captured by the model. Despite correctly describing the hallmark of the pseudogap at the two-particle level, such as the drop in the Knight shift of nuclear magnetic resonance, dynamical mean-field theory fails to capture the spectral properties of the pseudogap.

• Self-congruent point in critical matrix product states: An effective field theory for finite-entanglement scaling

  Jan T. Schneider, Atsushi Ueda, Yifan Liu, Andreas M. Läuchli, Masaki Oshikawa, Luca Tagliacozzo

  SciPost Phys. 18, 142 (2025) · published 29 April 2025 | Toggle abstract · pdf

  We set up an effective field theory formulation for the renormalization flow of matrix product states (MPS) with finite bond dimension, focusing on systems exhibiting finite-entanglement scaling close to a conformally invariant critical fixed point. We show that the finite MPS bond dimension $\chi$ is equivalent to introducing a perturbation by a relevant operator to the fixed-point Hamiltonian. The fingerprint of this mechanism is encoded in the $\chi$-independent universal transfer matrix's gap ratios, which are distinct from those predicted by the unperturbed Conformal Field Theory (CFT). This phenomenon defines a renormalization group self-congruent point, where the relevant coupling constant ceases to flow due to a balance of two effects; When increasing $\chi$, the infrared scale, set by the correlation length $\xi(\chi)$, increases, while the strength of the perturbation at the lattice scale decreases. The presence of a self-congruent point does not alter the validity of the finite-entanglement scaling hypothesis, since the self-congruent point is located at a finite distance from the critical fixed point, well inside the scaling regime of the CFT. We corroborate this framework with numerical evidence from the exact solution of the Ising model and density matrix renormalization group (DMRG) simulations of an effective lattice model.

  1 citation

• Interaction quenches in Bose gases studied with a time-dependent hypernetted-chain Euler-Lagrange method

  Mathias Gartner, David Miesbauer, Michael Kobler, Julia Freund, Giuseppe Carleo, Robert E. Zillich

  SciPost Phys. 18, 123 (2025) · published 9 April 2025 | Toggle abstract · pdf

  We present a new variational method to study the dynamics of a closed bosonic many-body system, the time-dependent hypernetted-chain Euler-Lagrange method, tHNC. Based on the Jastrow ansatz, it accounts for quantum fluctuations in a non-perturbative way. The tHNC method scales well with the number of dimensions, as demonstrated by our results on one-, two-, and three-dimensional systems. We apply the tHNC method to interaction quenches, i.e. sudden changes of the interaction strength, in homogeneous Bose gases. When the quench is strong enough that the final state has roton excitations (as found and predicted for dipolar and Rydberg-dressed Bose-Einstein condensates, respectively), the pair distribution function exhibits stable oscillations. For validation, we compare tHNC results with time-dependent variational Monte Carlo results in one and two dimensions.

• Clock factorized quantum Monte Carlo method for long-range interacting systems

  Zhijie Fan, Chao Zhang, Youjin Deng

  SciPost Phys. Core 8, 036 (2025) · published 8 April 2025 | Toggle abstract · pdf

  Simulating long-range interacting systems is a challenging task due to its computational complexity that the computational effort for each local update is of order $\mathcal{O}(N)$, where $N$ is the size of the system. In this work, we introduce the clock factorized quantum Monte Carlo method, an efficient technique for simulating long-range interacting quantum systems. The method is developed by generalizing the clock Monte Carlo method for classical systems [Phys. Rev. E 99 010105 (2019)] to the path-integral representation of long-range interacting quantum systems, with some specific treatments for quantum cases and a few significant technical improvements in general. We first explain how the clock factorized quantum Monte Carlo method is implemented to reduce the computational overhead from $\mathcal{O}(N)$ to $\mathcal{O}(1)$. In particular, the core ingredients, including the concepts of bound probabilities and bound rejection events, the recursive sampling procedure, and the fast algorithms for sampling an extensive set of discrete and small probabilities, are elaborated. Next, we show how the clock factorized quantum Monte Carlo method can be flexibly implemented in various update strategies, like the Metropolis and worm-type algorithms. Finally, we demonstrate the high efficiency of the clock factorized quantum Monte Carlo algorithms using examples of three typical long-range interacting quantum systems, including the transverse field Ising model with long-range $z$-$z$ interaction, the extended Bose-Hubbard model with long-range density-density interactions, and the XXZ Heisenberg model with long-range spin interactions. We expect that the clock factorized quantum Monte Carlo method would find broad applications in statistical and condensed-matter physics.

  2 citations