Search

Predicting the fate of an interacting system in the limit where the electronic bandwidth is quenched is often highly non-trivial. The complex interplay between interactions and quantum fluctuations driven by the band geometry can drive competition between various ground states, such as charge density wave order and superconductivity. In this work, we study an electronic model of topologically-trivial flat bands with a continuously tunable Fubini-Study metric in the presence of on-site attraction and nearest-neighbor repulsion, using numerically exact quantum Monte Carlo simulations. By varying the electron filling and the minimal spatial extent of the localized flat-band Wannier wavefunctions, we obtain a number of intertwined orders. These include a phase with coexisting charge density wave (CDW) order and superconductivity, i.e., a supersolid. In spite of the non-perturbative nature of the problem, we identify an analytically tractable limit associated with a 'small' spatial extent of the Wannier functions and derive a low-energy effective Hamiltonian that can well describe our numerical results. We also provide unambiguous evidence for the violation of any putative lower bound on the zero-temperature superfluid stiffness in geometrically non-trivial flat bands., 5 pages, 3 figures; Supplement with 4 pages, 2 figures; (v2) Includes additional discussion and references

We report the experimental observation of dielectric relaxation by quantum critical magnons. Complex capacitance measurements reveal a dissipative feature with a temperature-dependent amplitude due to low-energy lattice excitations and an activation behavior of the relaxation time. The activation energy softens close to a field-tuned magnetic quantum critical point at $H=H_c$ and follows single-magnon energy for $H>H_c$, showing its magnetic origin. Our study demonstrates the electrical activity of coupled low-energy spin and lattice excitations, an example of quantum multiferroic behavior., 6 pages, 4 figures (SM 7 pages, 6 figures)

We argue that spin and valley-polarized metallic phases recently observed in graphene bilayers and trilayers support chiral edge modes that allow spin waves to propagate ballistically along system boundaries without backscattering. The chiral edge behavior originates from the interplay between the momentum-space Berry curvature in Dirac bands and the geometric phase of a spin texture in position space. The edge modes are weakly confined to the edge, featuring dispersion which is robust and insensitive to the detailed profile of magnetization at the edge. This unique character of edge modes reduces their overlap with edge disorder and enhances the mode lifetime. The mode propagation direction reverses upon reversing valley polarization, an effect that provides a clear testable signature of geometric interactions in isospin-polarized Dirac bands., 8pgs, 2fgs

We show that bilayer graphene in the presence of a 2D superlattice potential provides a highly tunable setup that can realize a variety of flat band phenomena. We focus on two regimes: (i) topological flat bands with non-zero Chern numbers, C, including bands with higher Chern numbers |C| > 1; and (ii) an unprecedented phase consisting of a stack of nearly perfect flat bands with C = 0. For realistic values of the potential and superlattice periodicity, this stack can span nearly 100 meV, encompassing nearly all of the low-energy spectrum. We further show that in the topological regime, the topological flat band has a favorable band geometry for realizing a fractional Chern insulator (FCI) and use exact diagonalization to show that the FCI is in fact the ground state at 1/3 filling. Our results provide a realistic guide for future experiments to realize a new platform for flat band phenomena., Comment: 14 pages including supplementary materials. Published version

Excitonic Bose-Einstein condensation (EBEC) has drawn increasing attention recently with the emergence of 2D materials. A general criterion for EBEC, as expected in an excitonic insulator (EI) state, is to have negative exciton formation energies in a semiconductor. Here, using exact diagonalization of multi-exciton Hamiltonian modelled in a diatomic Kagome lattice, we demonstrate that the negative exciton formation energies are only a prerequisite but insufficient condition for realizing an EI. By a comparative study between the cases of both a conduction and valence flat bands (FBs) versus that of a parabolic conduction band, we further show that the presence and increased FB contribution to exciton formation provide an attractive avenue to stabilize the EBEC, as confirmed by calculations and analyses of multi-exciton energies, wave functions and reduced density matrices. Our results warrant a similar many-exciton analysis for other known/new candidates of EIs, and demonstrate the FBs of opposite parity as a unique platform for studying exciton physics, paving the way to material realization of spinor BEC and spin-superfluidity.

We show that the Klein bottle entropy [Phys. Rev. Lett. 119, 261603 (2017)] for conformal field theories (CFTs) perturbed by a relevant operator is a universal function of the dimensionless coupling constant. The universal scaling of the Klein bottle entropy near criticality provides an efficient approach to extract the scaling dimension of lattice operators via data collapse. As paradigmatic examples, we validate the universal scaling of the Klein bottle entropy for Ising and Z3 parafermion CFTs with various perturbations using numerical simulation with continuous matrix product operator approach., 5+5 pages, 3 figures, published version

We study the low-energy limit of Wilson lines (charged impurities) in conformal gauge theories in 2+1 and 3+1 dimensions. As a function of the representation of the Wilson line, certain defect operators can become marginal, leading to interesting renormalization group flows and for sufficiently large representations to complete or partial screening by charged fields. This result is universal: in large enough representations, Wilson lines are screened by the charged matter fields. We observe that the onset of the screening instability is associated with fixed-point mergers. We study this phenomenon in a variety of applications. In some cases, the screening of the Wilson lines takes place by dimensional transmutation and the generation of an exponentially large scale. We identify the space of infrared conformal Wilson lines in weakly coupled gauge theories in 3+1 dimensions and determine the screening cloud due to bosons or fermions. We also study QED in 2+1 dimensions in the large $N_f$ limit and identify the nontrivial conformal Wilson lines. We briefly discuss 't Hooft lines in 3+1-dimensional gauge theories and find that they are screened in weakly coupled gauge theories with simply connected gauge groups. In non-Abelian gauge theories with S-duality, this together with our analysis of the Wilson lines gives a compelling picture for the screening of the line operators as a function of the coupling., Comment: 5 pages, 2 figures, v2 references updated, v3 minor changes, journal version

The topological superconducting state is a highly sought-after quantum state hosting topological order and Majorana excitations. In this work, we explore the mechanism to realize the topological superconductivity (TSC) in the doped Mott insulators with time-reversal symmetry (TRS). Through large-scale density matrix renormalization group study of an extended triangular-lattice $t$-$J$ model on the six- and eight-leg cylinders, we identify a $d+id$-wave chiral TSC with spontaneous TRS breaking, which is characterized by a Chern number $C=2$ and quasi-long-range superconducting order. We map out the quantum phase diagram by tuning the next-nearest-neighbor (NNN) electron hopping and spin interaction. In the weaker NNN-coupling regime, we identify a pseudogaplike phase with a charge stripe order coexisting with fluctuating superconductivity, which can be tuned into $d$-wave superconductivity by increasing the doping level and system width. The TSC emerges in the intermediate-coupling regime, which has a transition to a $d$-wave superconducting phase with larger NNN couplings. The emergence of the TSC is driven by geometrical frustrations and hole dynamics which suppress spin correlation and charge order, leading to a topological quantum phase transition., 6+12 pages, 5+11 figures

Using {\it ab initio} approaches for extended Hubbard interactions coupled to phonons, we reveal that the intersite Coulomb interaction plays important roles in determining various distinctive phases of the paradigmatic charge ordered materials of Ba$_{1-x}$K$_x A$O$_3$ ($A=$ Bi and Sb). We demonstrated that all their salient doping dependent experiment features such as breathing instabilities, anomalous phonon dispersions, and transition between charge-density wave and superconducting states can be accounted very well if self-consistently obtained nearest neighbor Hubbard interaction are included, thus establishing a minimal criterion for reliable descriptions of spontaneous charge orders in solids., Comment: 4 pages, 2 additional pages for references and 4 pages supplementary materials, title and abstract are modified

In gauge theory, it is commonly stated that time-reversal symmetry only exists at $\theta=0$ or $\pi$ for a $2\pi$-periodic $\theta$-angle. In this paper, we point out that in both the free Maxwell theory and massive QED, there is a non-invertible time-reversal symmetry at every rational $\theta$-angle, i.e., $\theta= \pi p/N$. The non-invertible time-reversal symmetry is implemented by a conserved, anti-linear operator without an inverse. It is a composition of the naive time-reversal transformation and a fractional quantum Hall state. We also find similar non-invertible time-reversal symmetries in non-Abelian gauge theories, including the $\mathcal{N}=4$ $SU(2)$ super Yang-Mills theory along the locus $|\tau|=1$ on the conformal manifold., Comment: 30 pages

Strongly correlated layered 2D systems are of central importance in condensed matter physics, but their numerical study is very challenging. Motivated by the enormous successes of tensor networks for 1D and 2D systems, we develop an efficient tensor network approach based on infinite projected entangled-pair states (iPEPS) for layered 2D systems. Starting from an anisotropic 3D iPEPS ansatz, we propose a contraction scheme in which the weakly-interacting layers are effectively decoupled away from the center of the layers, such that they can be efficiently contracted using 2D contraction methods while keeping the center of the layers connected in order to capture the most relevant interlayer correlations. We present benchmark data for the anisotropic 3D Heisenberg model on a cubic lattice, which shows close agreement with quantum Monte Carlo (QMC) and full 3D contraction results. Finally, we study the dimer to N\'eel phase transition in the Shastry-Sutherland model with interlayer coupling, a frustrated spin model which is out of reach of QMC due to the negative sign problem., Comment: 7 pages, 5 figures, 6 pages of supplemental material

Motivated by the observation of even denominator fractional quantum Hall effect in the $n=3$ Landau level of monolayer graphene [Y. Kim $\textit{et al.}$, Nature Physics $\textbf{15}$, 154 (2019)], we consider a Bardeen-Cooper-Schrieffer variational state for composite fermions and find that the composite-fermion Fermi sea in this Landau level is unstable to an $f$-wave pairing. Analogous calculation suggests the possibility of a $p$-wave pairing of composite fermions at half filling in the $n=2$ graphene Landau level, whereas no pairing instability is found at half filling in the $n=0$ and $1$ graphene Landau levels. The relevance of these results to experiments is discussed., Comment: 13 pages, 7 figures

Fractional quantum Hall states (FQHSs) at even-denominator Landau level filling factors ($\nu$) are of prime interest as they are predicted to host exotic, topological states of matter. We report here the observation of a FQHS at $\nu=1/2$ in a two-dimensional electron system of exceptionally high quality, confined to a wide AlAs quantum well, where the electrons can occupy multiple conduction-band valleys with an anisotropic effective mass. The anisotropy and multi-valley degree of freedom offer an unprecedented tunability of the $\nu=1/2$ FQHS as we can control both the valley occupancy via the application of in-plane strain, and the ratio between the strengths of the short- and long-range Coulomb interaction by tilting the sample in the magnetic field to change the electron charge distribution. Thanks to this tunability, we observe phase transitions from a compressible Fermi liquid to an incompressible FQHS and then to an insulating phase as a function of tilt angle. We find that this evolution and the energy gap of the $\nu=1/2$ FQHS depend strongly on valley occupancy.

Skyrmions are of interest both from a fundamental and technological point of view, due to their potential to act as information carriers. But one challenge concerns their manipulation, especially at high temperature where thermal fluctuations eventually disintegrate them. Here we study the competition between skyrmions and a chiral spin liquid, using the latter as an entropic buffer to impose a quasi-vacuum of skyrmions. As a result, the temperature becomes a knob to tune the skyrmion density from a dense liquid to a diluted gas, protecting the integrity of each skyrmion from paramagnetic disintegration. With this additional knob in hand, we find at high field a topological spin glass made of zero- and one-dimensional topological defects (resp. skyrmions and bimerons)., 6 pages, 4 figures. Accepted in Physical Review Letters

Many-body entanglement in condensed matter systems can be diagnosed from equilibrium response functions through the use of entanglement witnesses and operator-specific quantum bounds. Here, we investigate the applicability of this approach for detecting entangled states in quantum systems driven out of equilibrium. We use a multipartite entanglement witness, the quantum Fisher information, to study the dynamics of a paradigmatic fermion chain undergoing a time-dependent change of the Coulomb interaction. Our results show that the quantum Fisher information is able to witness distinct signatures of multipartite entanglement both near and far from equilibrium that are robust against decoherence. We discuss implications of these findings for probing entanglement in light-driven quantum materials with time-resolved optical and x-ray scattering methods., Comment: 8 pages, 4 figures

We extend the concept of dual unitary quantum gates to quantum lattice models in $2 + 1$ dimensions, by introducing and studying ternary unitary four-particle gates, which are unitary in time and both spatial dimensions. When used as building blocks of lattice models with periodic boundary conditions in time and space (corresponding to infinite temperature states), dynamical correlation functions exhibit a light-ray structure. We also generalize solvable MPS to two spatial dimensions with cylindrical boundary conditions, by showing that the analogous solvable PEPS can be identified with matrix product unitaries. In the resulting tensor network for evaluating equal-time correlation functions, the bulk ternary unitary gates cancel out. We delineate and implement a numerical algorithm for computing such correlations by contracting the remaining tensors., Comment: 16 pages, 6 figures

We study two-dimensional materials where electrons are coupled to the vacuum electromagnetic field of a cavity. We show that, at the onset of the superradiant phase transition towards a macroscopic photon occupation of the cavity, the critical electromagnetic fluctuations, consisting of photons strongly overdamped by their interaction with electrons, can in turn lead to the absence of electronic quasiparticles. Since transverse photons couple to the electronic current, the appearance of non-Fermi-Liquid behavior strongly depends on the lattice. In particular, we find that in a square lattice the phase space for electron-photon scattering is reduced in such a way to preserve the quasiparticles, while in a honeycomb lattice the latter are removed due to a non-analytical frequency dependence of the damping $\propto |\omega|^{2/3}$. Standard cavity probes could allow to measure the characteristic frequency spectrum of the overdamped critical electromagnetic modes responsible for the non-Fermi-liquid behavior., Comment: main text: 7 pages, 2 figures; supplemental: 8 pages

We study theoretically a moir\'e superlattice geometry consisting of mirror-symmetric twisted trilayer graphene surrounded by identical transition metal dichalcogenide layers. We show that this setup allows to switch on/off and control the spin-orbit splitting of the Fermi surfaces via application of a perpendicular displacement field $D_0$, and explore two manifestations of this control: first, we compute the evolution of superconducting pairing with $D_0$; this features a complex admixture of singlet and triplet pairing and, depending on the pairing state in the parent trilayer system, phase transitions between competing superconducting phases. Second, we reveal that, with application of $D_0$, the spin-orbit-induced spin textures exhibit vortices which lead to "M\"obius fermi surfaces'' in the interior of the Brillouin zone: diabatic electron trajectories, which are predicted to dominate quantum oscillation experiments, require encircling the $\Gamma$ point twice, making their M\"obius nature directly observable. We further show that the superconducting order parameter inherits the unconventional, M\"obius spin textures. Our findings suggest that this system provides a promising experimental avenue for studying systematically the impact of spin-orbit coupling on the multitude of topological and correlated phases in near-magic-angle twisted trilayer graphene., Comment: 17 pages, 8 figures, 1 table

We construct an exactly solvable lattice model for a deconfined quantum critical point (DQCP) in (1+1) dimensions. This DQCP occurs in an unusual setting, namely at the edge of a (2+1) dimensional bosonic symmetry protected topological phase (SPT) with $\mathbb{Z}_2\times\mathbb{Z}_2$ symmetry. The DQCP describes a transition between two gapped edges that break different $\mathbb{Z}_2$ subgroups of the full $\mathbb{Z}_2\times\mathbb{Z}_2$ symmetry. Our construction is based on an exact mapping between the SPT edge theory and a $\mathbb{Z}_4$ spin chain. This mapping reveals that DQCPs in this system are directly related to ordinary $\mathbb{Z}_4$ symmetry breaking critical points., 4.25 pages with 3 figures + 3 page supplemental material with 1 figure. References, discussion, and intro updated

While the ground state phase diagram of the correlated flat-band systems have been intensively investigated, the dynamic and thermodynamic properties of such lattice models are less explored, but it is the latter which is most relevant to the experimental probes (transport, quantum capacitance and spectroscopy) of the quantum moir\'e materials such as twisted bilayer graphene and transition metal dichalcogenides. Here we show, by means of momentum-space quantum Monte Carlo and exact diagonalization, there exists a unique thermodynamic characteristic for the correlated flat-band models with interaction-driven quantum anomalous Hall (QAH) ground state, namely, the transition from the QAH insulator to the metallic state takes place at a much lower temperature compared with the zero-temperature single-particle gap generated by the long-range Coulomb interaction. Such low transition temperature comes from the proliferation of excitonic particle-hole excitations, which "quantum teleport" the electrons across the gap between different topological bands to restore the broken time-reversal symmetry and give rise to a pronounced enhancement in the charge compressibility. Future experiments, to verify such generic thermodynamic characteristics, are proposed., Comment: 5 pages, 3 figures with supplemental material

Quantum oscillation measurements in the new unconventional superconductor YFe$_2$Ge$_2$ resolve all four Fermi surface pockets expected from band structure calculations, which predict an electron pocket in the Brillouin zone corner and three hole pockets enveloping the centers of the top and bottom of the Brillouin zone. The carrier masses are uniformly renormalized by about a factor of five and broadly account for the enhanced heat capacity Sommerfeld coefficient $\simeq 100$ mJ/molK$^2$. Our data highlight the key role of the electron pocket, which despite its small volume accounts for about half the total density of states, and point towards a predominantly local mechanism underlying the mass renormalization in YFe$_2$Ge$_2$., Comment: 5 pages, 5 figures

In the presence of crystalline symmetries, topological phases of matter acquire a host of invariants leading to non-trivial quantized responses. Here we study a particular invariant, the discrete shift $\mathscr{S}$, for the square lattice Hofstadter model of free fermions. $\mathscr{S}$ is associated with a $\mathbb{Z}_M$ classification in the presence of $M$-fold rotational symmetry and charge conservation. $\mathscr{S}$ gives quantized contributions to (i) the fractional charge bound to a lattice disclination, and (ii) the angular momentum of the ground state with an additional, symmetrically inserted magnetic flux. $\mathscr{S}$ forms its own `Hofstadter butterfly', which we numerically compute, refining the usual phase diagram of the Hofstadter model. We propose an empirical formula for $\mathscr{S}$ in terms of density and flux per plaquette for the Hofstadter bands, and we derive a number of general constraints. We show that bands with the same Chern number may have different values of $\mathscr{S}$, although odd and even Chern number bands always have half-integer and integer values of $\mathscr{S}$ respectively., Comment: 6+14 pages. Minor edits to the main text and the appendix

In this work, we establish a new theoretical paradigm for vortex Majorana physics in the recently discovered topological iron-based superconductors (tFeSCs). While tFeSCs are widely accepted as an exemplar of topological insulators (TIs) with intrinsic $s$-wave superconductivity, our theory implies that such common belief could be oversimplified. Our main finding is that the normal-state bulk Dirac nodes, usually ignored in TI-based vortex Majorana theories for tFeSCs, will play a key role of determining the vortex state topology. In particular, the interplay between TI and Dirac nodal bands will lead to multiple competing topological phases for a superconducting vortex line in tFeSCs, including an unprecedented hybrid topological vortex state that carries both Majorana bound states and a gapless dispersion. Remarkably, this exotic hybrid vortex phase generally exists in the vortex phase diagram for our minimal model for tFeSCs and is directly relevant to tFeSC candidates such as LiFeAs.When the four-fold rotation symmetry is broken by vortex-line tilting or curving, the hybrid vortex gets topologically trivialized and becomes Majorana-free, which could explain the puzzle of ubiquitous trivial vortices observed in LiFeAs. Origin of Majorana signal in other tFeSC candidates such as FeTe$_x$Se$_{1-x}$ and CaKFe$_4$As$_4$ are also interpreted within our theory framework. Our theory sheds new light on theoretically understanding and experimentally engineering Majorana physics in high-temperature iron-based systems., Comment: 7+19 pages, 4+9 figures

We develop a nonperturbative theory for hole dynamics in antiferromagnetic spin lattices, as described by the $t$-$J$ model. This is achieved by generalizing the selfconsistent Born approximation to nonequilibrium systems, making it possible to calculate the full time-dependent many-body wave function. Our approach reveals three distinct dynamical regimes, ultimately leading to the formation of magnetic polarons. Following the initial ballistic stage of the hole dynamics, coherent formation of string excitations gives rise to characteristic oscillations in the hole density. Their damping eventually leaves behind magnetic polarons that undergo ballistic motion with a greatly reduced velocity. The developed theory provides a rigorous framework for understanding nonequilibrium physics of defects in quantum magnets and quantitatively explains recent observations from cold-atom quantum simulations in the strong coupling regime., 17 pages, 7 figures

We study near-neighbour and dipolar Ising models on a lattice of corner-sharing octahedra. In an extended parameter range of both models, frustration between antiferromagnetism and a spin-ice-like three-in-three-out rule stabilises a Coulomb phase with correlated dipolar and quadrupolar spin textures, both yielding distinctive neutron-scattering signatures. Strong further-neighbour perturbations cause the two components to order independently, resulting in unusual multi-$k$ orders. We propose experimental realisations of our model in rare-earth antiperovskites., 6+10 pages, 5+5 figures with supplementary text; ancillary mcif files illustrate the ordered phases

We study pole skipping in holographic CFTs dual to diffeomorphism invariant theories containing an arbitrary number of bosonic fields in the large $N$ limit. Defining a weight to organize the bulk equations of motion, a set of general pole-skipping conditions are derived. In particular, the frequencies simply follow from general covariance and weight matching. In the presence of higher spin fields, we find that the imaginary frequency for the highest-weight pole-skipping point equals the higher-spin Lyapunov exponent which lies outside of the chaos bound. Without higher spin fields, we show that the energy density Green's function has its highest-weight pole skipping happening at a location related to the OTOC for arbitrary higher-derivative gravity, with a Lyapunov exponent saturating the chaos bound and a butterfly velocity matching that extracted from a shockwave calculation. We also suggest an explanation for this matching at the metric level by obtaining the on-shell shockwave solution from a regularized limit of the metric perturbation at the skipped pole., Comment: 7 pages; v2: comments added

We study the electrodynamics of spin triplet superconductors including dipolar interactions, which give rise to an interplay between the collective spin dynamics of the condensate and orbital Meissner screening currents. Within this theory, we identify a class of spin waves that originate from the coupled dynamics of the spin-symmetry breaking triplet order parameter and the electromagnetic field. In particular, we study magnetostatic spin wave modes that are localized to the sample surface. We show that these surface modes can be excited and detected using experimental techniques such as microwave spin wave resonance spectroscopy or nitrogen-vacancy magnetometry, and propose that the detection of these modes offers a means for the identification of spin triplet superconductivity., Comment: 13 pages, 8 figures

Competition among exchange interactions is able to induce novel spin correlations on a bipartite lattice without geometrical frustration. A prototype example is the spiral spin liquid, which is a correlated paramagnetic state characterized by sub-dimensional degenerate propagation vectors. Here, using spectral graph theory, we show that spiral spin liquids on a bipartite lattice can be approximated by a further-neighbor model on the corresponding line-graph lattice that is non-bipartite, thus broadening the space of candidate materials that may support the spiral spin liquid phases. As illustrations, we examine neutron scattering experiments performed on two spinel compounds, ZnCr$_2$Se$_4$ and CuInCr$_4$Se$_8$, to demonstrate the feasibility of this new approach and expose its possible limitations in experimental realizations., Comment: 20 pages, 18 figures

There is a number of contradictory findings with regard to whether the theory describing easy-plane quantum antiferromagnets undergoes a second-order phase transition. The traditional Landau-Ginzburg-Wilson approach suggests a first-order phase transition, as there are two different competing order parameters. On the other hand, it is known that the theory has the property of self-duality which has been connected to the existence of a deconfined quantum critical point (DQCP). The latter regime suggests that order parameters are not the elementary building blocks of the theory, but rather consist of fractionalized particles that are confined in both phases of the transition and only appear - deconfine - at the critical point. Nevertheless, many numerical Monte Carlo simulations disagree with the claim of a DQCP in the system, indicating instead a first-order phase transition. Here we establish from exact lattice duality transformations and renormalization group analysis that the easy-plane CP1 antiferromagnet does feature a DQCP. We uncover the criticality starting from a regime analogous to the zero temperature limit of a certain classical statistical mechanics system which we therefore dub "frozen". At criticality our bosonic theory is dual to a fermionic one with two massless Dirac fermions, which thus undergoes a second-order phase transition as well., Comment: v2: published version includes Supplemental Material containing details of calculations and 2 figures; an error in original calculations was fixed; references added

The two-dimensional Wigner crystal (WC) occurs in the strongly interacting regime ($r_s \gg 1$) of the two-dimensional electron gas (2DEG). The magnetism of a pure WC is determined by tunneling processes that induce multi-spin ring-exchange interactions, resulting in fully polarized ferromagnetism for large enough $r_s$. Recently, Hossain et al. [PNAS 117 (51) 32244-32250] reported the occurrence of a fully polarized ferromagnetic insulator at $r_s \gtrsim 35$ in an AlAs quantum well, but at temperatures orders of magnitude larger than the predicted exchange energies for the pure WC. Here, we analyze the large $r_s$ dynamics of an interstitial defect in the WC, and show that it produces local ferromagnetism with much higher energy scales. Three hopping processes are dominant, which favor a large, fully polarized ferromagnetic polaron. Based on the above results, we speculate concerning the phenomenology of the magnetism near the metal-insulator transition of the 2DEG., 4 pages + references & Supplementary (5 pages), 2+4 figures, 1 table. Minor revision based on comments

We use vortex photon fields with orbital and spin angular momentum to probe chiral fluctuations within liquid crystals. In the regime of iridescence with a well-defined pitch length of chirality, we find low energy Raman scattering that can be decomposed into helical and chiral components depending on the scattering vector and the topological charge of the incident photon field. Based on the observation of an anomalous dispersion we attribute quasi-elastic scattering to a transfer of angular momenta to roton-like quasiparticles. The latter are due to a competition of short-range repulsive and long-range dipolar interactions. Our approach using a transfer of orbital angular momentum opens up an avenue for the advanced characterization of chiral and optically active devices and materials., Comment: 5 pages, 3 figures, Physical Rev. Lett., in print (2022)

A variety of novel quantum phases have been achieved in twist bilayer graphene (tBLG) and other moire superlattices recently, including correlated insulators, superconductivity, magnetism, and topological states. These phenomena are very sensitive to the moire superlattices, which can hardly be changed rapidly or intensely. Here, we report the experimental realization of a high-order moire pattern (a high-order interference pattern) in graphene-monolayer xenon heterostructure (G/mXe), with moire period in-situ tuned from few nanometers to infinity by changing the lattice constant of Xe through different annealing temperatures and pressures. We use angle-resolved photoemission spectroscopy to directly observe that replicas of graphene Dirac cone emerge and move close to each other in momentum-space as moire pattern continuously expands in real-space. When the moire period approaches infinity, the replicas finally overlap with each other and an energy gap is observed at the Dirac point induced by intervalley coupling, which is a manifestation of Kekule distortion. We construct a continuum moire Hamiltonian, which can explain the experimental results well. The form of moire Hamiltonian in G/mXe is similar to that in tBLG, and moire band with narrow bandwidth is predicted in G/mXe. However, the moire Hamiltonian couples Dirac fermions from different valleys in G/mXe, instead of ones from different layers in tBLG. Our work demonstrates a novel platform to study the continuous evolution of moire pattern and its modulation effect on electronic structure, and provides an unprecedented approach for tailoring Dirac fermions with tunable intervalley coupling., Comment: 17 pages, 4 figures, supplementary materials available from the authors, submitted Feb. 2022

Publisher Copyright: © 2022 American Physical Society. We present an analysis of positron lifetimes in solids with unprecedented depth. Instead of modeling correlation effects with density functionals, we study positron-electron wave functions with long-range correlations included. This gives new insight in understanding positron annihilation in metals, insulators, and semiconductors. By using a new quantum Monte Carlo approach for computation of positron lifetimes, an improved accuracy compared to previous computations is obtained for a representative set of materials when compared with experiment. Thus, we present a method without free parameters as a useful alternative to the already existing methods for modeling positrons in solids.

We reveal the crucial effect of strong spin-charge coupling on high-harmonic generation (HHG) in Mott insulators. In a system with antiferromagnetic correlations, the HHG signal is drastically enhanced with decreasing temperature, even though the gap increases and the production of charge carriers is suppressed. This anomalous behavior, which has also been observed in recent HHG experiments on Ca$_2$RuO$_4$, originates from a cooperative effect between the spin-charge coupling and the thermal ensemble, and the strongly temperature-dependent coherence between charge carriers. We argue that the peculiar temperature dependence of HHG is a generic feature of Mott insulators, which can be controlled via the Coulomb interaction and dimensionality of the system. Our results demonstrate that correlations between different degrees of freedom, which are a characteristic feature of strongly correlated solids, have significant and nontrivial effects on nonlinear optical responses., Comment: 6 pages, 4 figures + 9 pages, 12 figures

Persistent oscillatory dynamics in non-equilibrium many-body systems is a tantalizing manifestation of ergodicity breakdown that continues to attract much attention. Recent works have focused on two classes of such systems: discrete time crystals and quantum many-body scars (QMBS). While both systems host oscillatory dynamics, its origin is expected to be fundamentally different: discrete time crystal is a phase of matter which spontaneously breaks the $\mathbb{Z}_2$ symmetry of the external periodic drive, while QMBS span a subspace of non-thermalizing eigenstates forming an su(2) algebra representation. Here we ask a basic question: is there a physical system that allows to tune between these two dynamical phenomena? In contrast to much previous work, we investigate the possibility of a \emph{continuous} time crystal (CTC) in undriven, energy-conserving systems exhibiting prethermalization. We introduce a long-range XYZ spin model and show that it encompasses both a CTC phase as well as QMBS. We map out the dynamical phase diagram using numerical simulations based on exact diagonalization and time-dependent variational principle in the thermodynamic limit. We identify a regime where QMBS and CTC order co-exist, and we discuss experimental protocols that reveal their similarities as well as key differences., 7 pages, 3 figures. (6 pages, 5 figures in Supplementary Material)

Spontaneous symmetry breaking (SSB) is a property of Hamiltonian equilibrium states which, in the thermodynamic limit, retain a finite average value of an order parameter even after a field coupled to it is adiabatically turned off. In the case of quantum spin models with continuous symmetry, we show that this adiabatic process is also accompanied by the suppression of the fluctuations of the symmetry generator -- namely, the collective spin component along an axis of symmetry. In systems of $S=1/2$ spins or qubits, the combination of the suppression of fluctuations along one direction and of the persistence of transverse magnetization leads to spin squeezing -- a much sought-after property of quantum states, both for the purpose of entanglement detection as well as for metrological uses. Focusing on the case of XXZ models spontaneously breaking a U(1) (or even SU(2)) symmetry, we show that the adiabatically prepared states have nearly minimal spin uncertainty; that the minimum phase uncertainty that one can achieve with these states scales as $N^{-3/4}$ with the number of spins $N$; and that this scaling is attained after an adiabatic preparation time scaling linearly with $N$. Our findings open the door to the adiabatic preparation of strongly spin-squeezed states in a large variety of quantum many-body devices including e.g. optical lattice clocks., 5+3 pages, 3+6 figures

Magnetic materials show rich optical responses related to the magnetic order. These phenomena reflect the nature of their excitations, providing a powerful probe for the magnetic states and a way to control them. In recent years, such studies were extended to the optical control of spin current using nonlinear optical response similar to the photogalvanic effect. However, neither a candidate material nor a general formula for calculating the photogalvanic spin current is known so far. In this work, we develop a general theory for the photogalvanic spin current through a magnetic resonance process. Using the nonlinear response formalism, we find the nonlinear conductivity consists of two contributions that involve one and two magnon bands; the latter is a contribution unknown to date. We argue that the two-band process produces a large photogalvanic spin current in the antiferromagnetic phase of bilayer CrI$_3$ and CrBr$_3$, whose resonance frequency can be tuned between GHz-THz range by an external magnetic field. Our findings open a route to the studies on the photogalvanic effect of spin angular momentum in realistic setups., Comment: 4 pages, 2 figures

We introduce a model of non-unitary quantum dynamics that exhibits infinitely long-lived discrete spatiotemporal order robust against any unitary or dissipative perturbation. Ergodicity is evaded by combining a sequence of projective measurements with a local feedback rule that is inspired by Toom's `North-East-Center' classical cellular automaton. The measurements in question only partially collapse the wavefunction of the system, allowing some quantum coherence to persist. We demonstrate our claims using numerical simulations of a Clifford circuit in two spatial dimensions which allows access to large system sizes, and also present results for more generic dynamics on modest system sizes. We also devise explicit experimental protocols realising this dynamics using one- and two-qubit gates that are available on present-day quantum computing platforms., 6+7 pages, 2+6 figures

It is generally accepted that the dynamical mean field theory gives a good solution of the Holstein model, but only in dimensions greater than two. Here, we show that this theory, which becomes exact in the weak coupling and in the atomic limit, provides an excellent, numerically cheap, approximate solution for the spectral function of the Holstein model in the whole range of parameters, even in one dimension. To establish this, we make a detailed comparison with the spectral functions that we obtain using the newly developed momentum-space numerically exact hierarchical equations of motion method, which yields electronic correlation functions directly in real time. We crosscheck these conclusions with our path integral quantum Monte Carlo and exact diagonalization results, as well as with the available numerically exact results from the literature., 6 pages, 4 figures + supp.mat. 25 pages, 28 figures

It has recently been established that cluster-like states -- states that are in the same symmetry-protected topological phase as the cluster state -- provide a family of resource states that can be utilized for Measurement-Based Quantum Computation. In this work, we ask whether it is possible to prepare cluster-like states in finite time without breaking the symmetry protecting the resource state. Such a symmetry-preserving protocol would benefit from topological protection to errors in the preparation. We answer this question in the positive by providing a Hamiltonian in one higher dimension whose finite-time evolution is a unitary that acts trivially in the bulk, but pumps the desired cluster state to the boundary. Examples are given for both the 1D cluster state protected by a global symmetry, and various 2D cluster states protected by subsystem symmetries. We show that even if unwanted symmetric perturbations are present in the driving Hamiltonian, projective measurements in the bulk along with feedforward correction is sufficient to recover a cluster-like state., 6+6 pages. 5+4 figures. Published version: a feedforward protocol replaces the need for postselection along with other minor improvements and references added

Majorana zero modes are a much sought-after consequence of one-dimensional topological superconductivity. Here, we show that, in turn, zero modes accompanying dynamical instanton events strongly enhance -- in some cases even enable -- superconductivity. We find that the dynamics of a one-dimensional topological triplet superconductor is governed by a $\theta$ term in the action. For isotropic triplets, this term enables algebraic charge-$2e$ superconductivity, which is destroyed by fluctuations in nontopological superconductors. For anisotropic triplets, zero modes suppress quantum phase slips and stabilize superconductivity over a large region of the phase diagram. We present predictions of correlation functions and thermodynamics for states of topologically enhanced superconductivity., Comment: 15 pages, 3 figures; Accepted and revised version: the structure has been revised, the discussion of the anisotropic case has mostly been relegated to the appendix. A discussion of the influence of dynamical gauge fields, as well as some numerical results on dynamic zero modes have been added. The exponent of the isotropic single-fermion correlator has been corrected

It is of great current interest to establish toy models of ergodicity breaking transitions in quantum many-body systems. Here, we study a model that is expected to exhibit an ergodic to nonergodic transition in the thermodynamic limit upon tuning the coupling between an ergodic quantum dot and distant particles with spin-1/2. The model is effectively zero dimensional; however, a variant of the model was proposed by De Roeck and Huveneers to describe the avalanche mechanism of ergodicity breaking transition in one-dimensional disordered spin chains. We show that exact numerical results based on the spectral form factor calculation accurately agree with theoretical predictions, and hence unambiguously confirm existence of the ergodicity breaking transition in this model. We benchmark specific properties that represent hallmarks of the ergodicity breaking transition in finite systems.

We apply modular flow -- entanglement generated dynamics -- to characterize quantum orders of ground state wavefunctions. In particular, we study the linear response of the entanglement entropy of a simply connected region with respect to modular flow. First, we apply it to (1+1)D conformal field theories and demonstrate its relationship to the chiral central charge -- or equivalently the perturbative gravitational anomaly -- which is shown to vanish. Next, we apply it to (2+1)D gapped ground states where it reduces to a recently proposed formula by Kim et. al. that is conjectured to compute the edge chiral central charge. Modular flow provides an intuitive picture for this conjecture based on bulk-edge correspondence. We also provide numerics on free fermion models that lend support to our picture., 5 pages, 4 figures

We use resonant inelastic x-ray scattering (RIXS) to probe the propagation of plasmons in the electron-doped cuprate superconductor Sr$_{0.9}$La$_{0.1}$CuO$_2$ (SLCO). We detect a plasmon gap of $\sim$~120 meV at the two-dimensional Brillouin zone center, indicating that low-energy plasmons in SLCO are not strictly acoustic. The plasmon dispersion, including the gap, is accurately captured by layered $t$-$J$-$V$ model calculations. A similar analysis performed on recent RIXS data from other cuprates suggests that the plasmon gap is generic and its size is related to the magnitude of the interlayer hopping $t_z$. Our work signifies the three-dimensionality of the charge dynamics in layered cuprates and provides a new method to determine $t_z$., Comment: 17 pages, 10 figures, includes Supplemental Material. Accepted for publication in Physical Review Letters

Photoinduced phase transitions and charge dynamics in the interacting Dirac-electron system with a charge-ordered ground state are theoretically studied by taking an organic salt $\alpha$-(BEDT-TTF)$_2$I$_3$. By analysing the extended Hubbard model for this compound using a combined method of numerical simulations based on the time-dependent Schr\"odinger equation and the Floquet theory, we observe successive dynamical phase transitions from the charge-ordered insulator to a gapless Dirac semimetal and, eventually, to a Chern insulator phase under irradiation with circularly polarized light. These phase transitions occur as a consequence of two major effects of circularly polarized light, i.e., closing of the charge gap through melting the charge order and opening of the topological gap by breaking the time reversal symmetry at the Dirac points. We demonstrate that these photoinduced phenomena are governed by charge dynamics of driven correlated Dirac electrons., Comment: 6 pages, 4 figures

The interplay between the Fermi sea anisotropy, electron-electron interaction, and localization phenomena can give rise to exotic many-body phases. An exciting example is an anisotropic two-dimensional (2D) Wigner solid (WS), where electrons form an ordered array with an anisotropic lattice structure. Such a state has eluded experiments up to now as its realization is extremely demanding: First, a WS entails very low densities where the Coulomb interaction dominates over the kinetic (Fermi) energy. Attaining such low densities while keeping the disorder low is very challenging. Second, the low-density requirement has to be fulfilled in a material that hosts an anisotropic Fermi sea. Here, we report transport measurements in a clean (low-disorder) 2D electron system with anisotropic effective mass and Fermi sea. The data reveal that at extremely low electron densities, when the r_{s} parameter, the ratio of the Coulomb to the Fermi energy, exceeds ≃38, the current-voltage characteristics become strongly nonlinear at small dc biases. Several key features of the nonlinear characteristics, including their anisotropic voltage thresholds, are consistent with the formation of a disordered, anisotropic WS pinned by the ubiquitous disorder potential.

The Kitaev model is a fascinating example of an exactly solvable model displaying a spin-liquid ground state in two dimensions. However, deviations from the original Kitaev model are expected to appear in real materials. In this Letter, we investigate the fate of Kitaev's spin liquid in the presence of disorder -- bond defects or vacancies -- for an extended version of the model. Considering static flux backgrounds, we observe a power-law divergence in the low-energy limit of the density of states with a nonuniversal exponent. We link this power-law distribution of energy scales to weakly coupled droplets inside the bulk, in an uncanny similarity to the Griffiths phase often present in the vicinity of disordered quantum phase transitions. If time-reversal symmetry is broken, we find that power-law singularities are tied to the destruction of the topological phase of the Kitaev model in the presence of bond disorder alone. However, there is a transition from this topologically trivial phase with power-law singularities to a topologically nontrivial one for weak to moderate site dilution. Therefore, diluted Kitaev materials are potential candidates to host Kitaev's chiral spin-liquid phase., Comment: Main text: 7 pages, 2 figures. Supplemental material: 8 pages, 12 figures. Published version

We show how the index of the fermion operator from the Euclidean action can be used to uncover the existence of gapless modes living on defects (such as edges and vortices) in topological insulators and superconductors. The 1-loop Feynman diagram that computes the index reveals an analog of the Quantum Hall current flowing on and off the defect -- even in systems without conserved currents or chiral anomalies -- and makes explicit the interplay between topology in momentum and coordinate space. We provide several explicit examples.

The correlated insulator (CI) states and the recently discovered density wave (DW) states in magic-angle twisted bilayer graphene (TBG) have stimulated intense research interest. However, up to date, the nature of these "featureless" correlated states with zero Chern numbers are still elusive, and are lack of characteristic experimental signature. Thus, an experimental probe to identify the characters of these featureless CI and DW states are urgently needed. In this work, we theoretically study the correlated insulators and density-wave states at different integer and fractional fillings of the flat bands in magic-angle TBG based on extended unrestricted Hartree-Fock calculations including the Coulomb screening effects from the remote bands. We further investigate the nonlinear optical response of the various correlated states, and find that the nonlinear optical conductivities can be used to identify the nature of these CI and DW states at most of the fillings. Therefore, we propose that nonlinear optical response can serve as a promising experimental probe to unveil the nature of the CI and DW states observed in magic-angle TBG., 6+19 pages

The entanglement entropy (EE) encodes key properties of quantum many-body systems. It is usually calculated for subregions of finite volume (or area in 2d). In this work, we study the EE of skeletal regions that have \textit{no} volume, such as a line in 2d. We show that skeletal entanglement displays new behavior compared to its bulk counterpart, and leads to distinct universal quantities. We provide non-perturbative bounds for the skeletal area-law coefficient of a large family of quantum states. We then explore skeletal scaling for the toric code, conformal bosons and Dirac fermions, Lifshitz critical points, and Fermi liquids. We discover signatures including skeletal topological EE, novel corner terms, and strict area-law scaling for metals. These findings suggest that skeletal entropy serves as a measure for the range of entanglement. We discuss the possibility of a continuum description involving the fusion of defect operators. Finally, we outline open questions relating to other systems, and measures such as the logarithmic negativity., 5+8 pages, 3+2 figures; v2: discussion and presentation of our results improved; v3: match published version