SciPost: SciPost Physics: issue detail

SciPost Physics

SciPost Physics Vol. 11 issue 6 (December 2021)

Previous issue | Vol. 11 issue 5

Vol. 12 issue 1 | Next issue

GHZ-like states in the Qubit-Qudit Rabi model

Yuan Shen, Giampiero Marchegiani, Gianluigi Catelani, Luigi Amico, Ai Qun Liu, Weijun Fan, Leong-Chuan Kwek

SciPost Phys. 11, 099 (2021) · published 1 December 2021 | · pdf

We study a Rabi type Hamiltonian system in which a qubit and a d-level quantum system (qudit) are coupled through a common resonator. In the weak and strong coupling limits the spectrum is analysed through suitable perturbative schemes. The analysis show that the presence of the multilevels of the qudit effectively enhance the qubit-qudit interaction. The ground state of the strongly coupled system is found to be of Greenberger-Horne-Zeilinger (GHZ) type. Therefore, despite the qubit-qudit strong coupling, the nature of the specific tripartite entanglement of the GHZ state suppresses the bipartite entanglement. We analyze the system dynamics under quenching and adiabatic switching of the qubit-resonator and qudit-resonator couplings. In the quench case, we found that the non-adiabatic generation of photons in the resonator is enhanced by the number of levels in the qudit. The adiabatic control represents a possible route for preparation of GHZ states. Our analysis provides relevant information for future studies on coherent state transfer in qubit-qudit systems.
Optical lattice with spin-dependent sub-wavelength barriers

Edvinas Gvozdiovas, Povilas Račkauskas, Gediminas Juzeliūnas

SciPost Phys. 11, 100 (2021) · published 2 December 2021 | · pdf

We analyze a tripod atom light coupling scheme characterized by two dark states playing the role of quasi-spin states. It is demonstrated that by properly configuring the coupling laser fields, one can create a lattice with spin-dependent sub-wavelength barriers. This allows to flexibly alter the atomic motion ranging from atomic dynamics in the effective brick-wall type lattice to free motion of atoms in one dark state and a tight binding lattice with a twice smaller periodicity for atoms in the other dark state. Between the two regimes, the spectrum undergoes significant changes controlled by the laser fields. The tripod lattice can be produced using current experimental techniques. The use of the tripod scheme to create a lattice of degenerate dark states opens new possibilities for spin ordering and symmetry breaking.
Emergent fractal phase in energy stratified random models

Anton G. Kutlin, Ivan M. Khaymovich

SciPost Phys. 11, 101 (2021) · published 6 December 2021 | · pdf

We study the effects of partial correlations in kinetic hopping terms of long-range disordered random matrix models on their localization properties. We consider a set of models interpolating between fully-localized Richardson's model and the celebrated Rosenzweig-Porter model (with implemented translation-invariant symmetry). In order to do this, we propose the energy-stratified spectral structure of the hopping term allowing one to decrease the range of correlations gradually. We show both analytically and numerically that any deviation from the completely correlated case leads to the emergent non-ergodic delocalization in the system unlike the predictions of localization of cooperative shielding. In order to describe the models with correlated kinetic terms, we develop the generalization of the Dyson Brownian motion and cavity approaches basing on stochastic matrix process with independent rank-one matrix increments and examine its applicability to the above set of models.
Rotating 5D black holes: Interactions and deformations near extremality

Alejandra Castro, Juan F. Pedraza, Chiara Toldo, Evita Verheijden

SciPost Phys. 11, 102 (2021) · published 7 December 2021 | · pdf

We study a two-dimensional theory of gravity coupled to matter that is relevant to describe holographic properties of black holes with two equal angular momenta in five dimensions (with or without cosmological constant). We focus on the near-horizon geometry of the near-extremal black hole, where the effective theory reduces to Jackiw-Teitelboim (JT) gravity coupled to a massive scalar field. We compute the corrections to correlation functions due to cubic interactions present in this theory. A novel feature is that these corrections do not have a definite sign: for AdS$_5$ black holes the sign depends on the mass of the extremal solution. We discuss possible interpretations of these corrections from a gravitational and holographic perspective. We also quantify the imprint of the JT sector on the UV region, i.e. how these degrees of freedom, characteristic for the near-horizon region, influence the asymptotically far region of the black hole. This gives an interesting insight on how to interpret the IR modes in the context of their UV completion, which depends on the environment that contains the black hole.
Dual applications of Chebyshev polynomials method: Efficiently finding thousands of central eigenvalues for many-spin systems

Haoyu Guan, Wenxian Zhang

SciPost Phys. 11, 103 (2021) · published 8 December 2021 | · pdf

Computation of a large group of interior eigenvalues at the middle spectrum is an important problem for quantum many-body systems, where the level statistics provides characteristic signatures of quantum chaos. We propose an exact numerical method, dual applications of Chebyshev polynomials (DACP), to simultaneously find thousands of central eigenvalues, where the level space decreases exponentially with the system size. To disentangle the near-degenerate problem, we employ twice the Chebyshev polynomials, to construct an exponential semicircle filter as a preconditioning step and to generate a large set of proper basis states in the desired subspace. Numerical calculations on Ising spin chain and spin glass shards confirm the correctness and efficiency of DACP. As numerical results demonstrate, DACP is $30$ times faster than the state-of-the-art shift-invert method for the Ising spin chain while $8$ times faster for the spin glass shards. In contrast to the shift-invert method, the computation time of DACP is only weakly influenced by the required number of eigenvalues, which renders it a powerful tool for large scale eigenvalues computations. Moreover, the consumed memory also remains a small constant (5.6 GB) for spin-$1/2$ systems consisting of up to $20$ spins, making it desirable for parallel computing.
On computing non-equilibrium dynamics following a quench

Neil J. Robinson, Albertus J. J. M. de Klerk, Jean-Sébastien Caux

SciPost Phys. 11, 104 (2021) · published 9 December 2021 | · pdf

Computing the non-equilibrium dynamics that follows a quantum quench is difficult, even in exactly solvable models. Results are often predicated on the ability to compute overlaps between the initial state and eigenstates of the Hamiltonian that governs time evolution. Except for a handful of known cases, it is generically not possible to find these overlaps analytically. Here we develop a numerical approach to preferentially generate the states with high overlaps for a quantum quench starting from the ground state or an excited state of an initial Hamiltonian. We use these preferentially generated states, in combination with a "high overlap states truncation scheme" and a modification of the numerical renormalization group, to compute non-equilibrium dynamics following a quench in the Lieb-Liniger model. The method is non-perturbative, works for reasonable numbers of particles, and applies to both continuum and lattice systems. It can also be easily extended to more complicated scenarios, including those with integrability breaking.
Generalized eigenproblem without fermion doubling for Dirac fermions on a lattice

M. J. Pacholski, G. Lemut, J. Tworzydło, C. W. J. Beenakker

SciPost Phys. 11, 105 (2021) · published 14 December 2021 | · pdf

The spatial discretization of the single-cone Dirac Hamiltonian on the surface of a topological insulator or superconductor needs a special "staggered" grid, to avoid the appearance of a spurious second cone in the Brillouin zone. We adapt the Stacey discretization from lattice gauge theory to produce a generalized eigenvalue problem, of the form ${\mathcal H}\psi=E {\mathcal P}\psi$, with Hermitian tight-binding operators ${\mathcal H}$, ${\mathcal P}$, a locally conserved particle current, and preserved chiral and symplectic symmetries. This permits the study of the spectral statistics of Dirac fermions in each of the four symmetry classes A, AII, AIII, and D.
Exact relaxation to Gibbs and non-equilibrium steady states in the quantum cellular automaton Rule 54

Katja Klobas, Bruno Bertini

SciPost Phys. 11, 106 (2021) · published 20 December 2021 | · pdf

We study the out-of-equilibrium dynamics of the quantum cellular automaton Rule 54 using a time-channel approach. We exhibit a family of (non-equilibrium) product states for which we are able to describe exactly the full relaxation dynamics. We use this to prove that finite subsystems relax to a one-parameter family of Gibbs states. We also consider inhomogeneous quenches. Specifically, we show that when the two halves of the system are prepared in two different solvable states, finite subsystems at finite distance from the centre eventually relax to the non-equilibrium steady state (NESS) predicted by generalised hydrodynamics. To the best of our knowledge, this is the first exact description of the relaxation to a NESS in an interacting system and, therefore, the first independent confirmation of generalised hydrodynamics for an inhomogeneous quench.
Entanglement dynamics in Rule 54: Exact results and quasiparticle picture

Katja Klobas, Bruno Bertini

SciPost Phys. 11, 107 (2021) · published 20 December 2021 | · pdf

We study the entanglement dynamics generated by quantum quenches in the quantum cellular automaton Rule 54. We consider the evolution from a recently introduced class of solvable initial states. States in this class relax (locally) to a one-parameter family of Gibbs states and the thermalisation dynamics of local observables can be characterised exactly by means of an evolution in space. Here we show that the latter approach also gives access to the entanglement dynamics and derive exact formulas describing the asymptotic linear growth of all Rényi entropies in the thermodynamic limit and their eventual saturation for finite subsystems. While in the case of von Neumann entropy we recover exactly the predictions of the quasiparticle picture, we find no physically meaningful quasiparticle description for other Rényi entropies. Our results apply to both homogeneous and inhomogeneous quenches.
Conserved momenta of ferromagnetic solitons through the prism of differential geometry

Xingjian Di, Oleg Tchernyshyov

SciPost Phys. 11, 108 (2021) · published 21 December 2021 | · pdf

The relation between symmetries and conservation laws for solitons in a ferromagnet is complicated by the presence of gyroscopic (precessional) forces, whose description in the Lagrangian framework involves a background gauge field. This makes canonical momenta gauge-dependent and requires a careful application of Noether's theorem. We show that Cartan's theory of differential forms is a natural language for this task. We use it to derive conserved momenta of the Belavin--Polyakov skyrmion, whose symmetries include translation, global spin rotation, and dilation.
Engineering spectral properties of non-interacting lattice Hamiltonians

Ali G. Moghaddam, Dmitry Chernyavsky, Corentin Morice, Jasper van Wezel, Jeroen van den Brink

SciPost Phys. 11, 109 (2021) · published 22 December 2021 | · pdf

We investigate the spectral properties of one-dimensional lattices with position-dependent hopping amplitudes and on-site potentials that are smooth bounded functions of the position. We find an exact integral form for the density of states (DOS) in the limit of an infinite number of sites, which we derive using a mixed Bloch-Wannier basis consisting of piecewise Wannier functions. Next, we provide an exact solution for the inverse problem of constructing the position-dependence of hopping in a lattice model yielding a given DOS. We confirm analytic results by comparing them to numerics obtained by exact diagonalization for various incarnations of position-dependent hoppings and on-site potentials. Finally, we generalize the DOS integral form to multi-orbital tight-binding models with longer-range hoppings and in higher dimensions.
Full counting statistics for interacting trapped fermions

Naftali R. Smith, Pierre Le Doussal, Satya N. Majumdar, Grégory Schehr

SciPost Phys. 11, 110 (2021) · published 22 December 2021 | · pdf

We study $N$ spinless fermions in their ground state confined by an external potential in one dimension with long range interactions of the general Calogero-Sutherland type. For some choices of the potential this system maps to standard random matrix ensembles for general values of the Dyson index $\beta$. In the fermion model $\beta$ controls the strength of the interaction, $\beta=2$ corresponding to the noninteracting case. We study the quantum fluctuations of the number of fermions ${\cal N}_{\cal D}$ in a domain $\cal{D}$ of macroscopic size in the bulk of the Fermi gas. We predict that for general $\beta$ the variance of ${\cal N}_{\cal D}$ grows as $A_{\beta} \log N + B_{\beta}$ for $N \gg 1$ and we obtain a formula for $A_\beta$ and $B_\beta$. This is based on an explicit calculation for $\beta\in\left\{ 1,2,4\right\} $ and on a conjecture that we formulate for general $\beta$. This conjecture further allows us to obtain a universal formula for the higher cumulants of ${\cal N}_{\cal D}$. Our results for the variance in the microscopic regime are found to be consistent with the predictions of the Luttinger liquid theory with parameter $K = 2/\beta$, and allow to go beyond. In addition we present families of interacting fermion models in one dimension which, in their ground states, can be mapped onto random matrix models. We obtain the mean fermion density for these models for general interaction parameter $\beta$. In some cases the fermion density exhibits interesting transitions, for example we obtain a noninteracting fermion formulation of the Gross-Witten-Wadia model.
A roadmap for bootstrapping critical gauge theories: decoupling operators of conformal field theories in $d>2$ dimensions

Yin-Chen He, Junchen Rong, Ning Su

SciPost Phys. 11, 111 (2021) · published 22 December 2021 | · pdf

We propose a roadmap for bootstrapping conformal field theories (CFTs) described by gauge theories in dimensions $d>2$. In particular, we provide a simple and workable answer to the question of how to detect the gauge group in the bootstrap calculation. Our recipe is based on the notion of \emph{decoupling operator}, which has a simple (gauge) group theoretical origin, and is reminiscent of the null operator of $2d$ Wess-Zumino-Witten CFTs in higher dimensions. Using the decoupling operator we can efficiently detect the rank (i.e. color number) of gauge groups, e.g., by imposing gap conditions in the CFT spectrum. We also discuss the physics of the equation of motion, which has interesting consequences in the CFT spectrum as well. As an application of our recipes, we study a prototypical critical gauge theory, namely the scalar QED which has a $U(1)$ gauge field interacting with critical bosons. We show that the scalar QED can be solved by conformal bootstrap, namely we have obtained its kinks and islands in both $d=3$ and $d=2+\epsilon$ dimensions.