Showing 1–50 of 83 results for author: Vafek, O

Strongly interacting Hofstadter states in magic-angle twisted bilayer graphene

Authors: Minhao He, Xiaoyu Wang, Jiaqi Cai, Jonah Herzog-Arbeitman, Takashi Taniguchi, Kenji Watanabe, Ady Stern, B. Andrei Bernevig, Matthew Yankowitz, Oskar Vafek, Xiaodong Xu

Abstract: Magic-angle twisted bilayer graphene (MATBG) hosts a multitude of strongly correlated states at partial fillings of its flat bands. In a magnetic field, these flat bands further evolve into a unique Hofstadter spectrum renormalized by strong Coulomb interactions. Here, we study the interacting Hofstadter states spontaneously formed within the topological magnetic subbands of an ultraclean MATBG de… ▽ More Magic-angle twisted bilayer graphene (MATBG) hosts a multitude of strongly correlated states at partial fillings of its flat bands. In a magnetic field, these flat bands further evolve into a unique Hofstadter spectrum renormalized by strong Coulomb interactions. Here, we study the interacting Hofstadter states spontaneously formed within the topological magnetic subbands of an ultraclean MATBG device, notably including symmetry-broken Chern insulator (SBCI) states and fractional quantum Hall (FQH) states. The observed SBCI states form a cascade with their Chern numbers mimicking the main sequence correlated Chern insulators. The FQH states in MATBG form in Jain sequence; however, they disappear at high magnetic field, distinct from conventional FQH states which strengthen with increasing magnetic field. We reveal a unique magnetic field-driven phase transition from composite fermion phases to a dissipative Fermi liquid. Our theoretical analysis of the magnetic subbands hosting FQH states predicts non uniform quantum geometric properties far from the lowest Landau level. This points towards a more natural interpretation of these FQH states as in-field fractional Chern insulators of the magnetic subbands. △ Less

Submitted 2 August, 2024; originally announced August 2024.

Interacting phase diagram of twisted bilayer MoTe$_2$ in magnetic field

Authors: Minxuan Wang, Xiaoyu Wang, Oskar Vafek

Abstract: We study electron-electron interaction induced states of twisted bilayer MoTe$_2$ in an out-of-plane magnetic field $B\hat{\bf z}$ near one hole per moiré unit cell filling. The 3D phase diagram showing the evolution of competing phases with $B$, interaction strength and an out-of-plane electric field is presented at electron fillings that follow the Diophantine equation along Chern number… ▽ More We study electron-electron interaction induced states of twisted bilayer MoTe$_2$ in an out-of-plane magnetic field $B\hat{\bf z}$ near one hole per moiré unit cell filling. The 3D phase diagram showing the evolution of competing phases with $B$, interaction strength and an out-of-plane electric field is presented at electron fillings that follow the Diophantine equation along Chern number $-\text{sign}\left(B\right)$ line, that is pointing away from the charge neutral filling, where we find prominent Chern insulators consistent with the experiments. We also explain the experimental absence of prominent Chern insulators along the Chern number $+\text{sign}\left(B\right)$ line. △ Less

Submitted 3 June, 2024; v1 submitted 23 May, 2024; originally announced May 2024.

Comments: 6+9 pages, 4+7 figures typo in Fig 2g corrected, a reference added, minor edits

Heavy Fermions as an Efficient Representation of Atomistic Strain and Relaxation in Twisted Bilayer Graphene

Authors: Jonah Herzog-Arbeitman, Jiabin Yu, Dumitru Călugăru, Haoyu Hu, Nicolas Regnault, Oskar Vafek, Jian Kang, B. Andrei Bernevig

Abstract: Although the strongly interacting flat bands in twisted bilayer graphene (TBG) have been approached using the minimal Bistritzer-MacDonald (BM) Hamiltonian, there is mounting evidence that strain and lattice relaxation are essential in correctly determining the order of the correlated insulator groundstates. These effects can be incorporated in an enhanced continuum model by introducing additional… ▽ More Although the strongly interacting flat bands in twisted bilayer graphene (TBG) have been approached using the minimal Bistritzer-MacDonald (BM) Hamiltonian, there is mounting evidence that strain and lattice relaxation are essential in correctly determining the order of the correlated insulator groundstates. These effects can be incorporated in an enhanced continuum model by introducing additional terms computed from the relaxation profile. To develop an analytical and physical understanding of these effects, we include strain and relaxation in the topological heavy fermion (HF) model of TBG. We find that strain and relaxation are very well captured in first order perturbation theory by projection onto the fully symmetric HF Hilbert space, and remarkably do not alter the interacting terms in the periodic Anderson model. Their effects are fully incorporated in the single-particle HF Hamiltonian, and can be reproduced in a minimal model with only 4 symmetry-breaking terms. Our results demonstrate that the heavy fermion framework of TBG is an efficient and robust representation of the perturbations encountered in experiment. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Topological Heavy Fermion Principle For Flat (Narrow) Bands With Concentrated Quantum Geometry

Authors: Jonah Herzog-Arbeitman, Jiabin Yu, Dumitru Călugăru, Haoyu Hu, Nicolas Regnault, Chaoxing Liu, Sankar Das Sarma, Oskar Vafek, Piers Coleman, Alexei Tsvelik, Zhi-da Song, B. Andrei Bernevig

Abstract: We propose a general principle for the low-energy theory of narrow bands with concentrated Berry curvature and Fubini-Study metric in the form of a map to Anderson-"+" models composed of heavy fermions hybridizing and interacting with semi-metallic modes. This map resolves the obstruction preventing topological bands from being realized in a local Hamiltonian acting on the low-energy degrees of fr… ▽ More We propose a general principle for the low-energy theory of narrow bands with concentrated Berry curvature and Fubini-Study metric in the form of a map to Anderson-"+" models composed of heavy fermions hybridizing and interacting with semi-metallic modes. This map resolves the obstruction preventing topological bands from being realized in a local Hamiltonian acting on the low-energy degrees of freedom. The concentrated quantum geometry is reproduced through band inversion with a dispersive semi-metal, leaving a nearly flat, trivial band which becomes the heavy fermion. This representation is natural when the narrow band is not energetically isolated on the scale of the interaction and an enlarged Hilbert space is inescapable, but also provides analytical insight into the projected-interaction limit. First exemplified in twisted bilayer graphene (TBG), we extend it to (1) the twisted checkerboard, which we find has a chiral symmetric stable anomaly that forbids a lattice realization at all energies, and (2) the Lieb lattice with gapless flat bands, where we show the heavy fermions can be obtained by minimizing a Euclidean instanton action to saturate its BPS bound. The heavy fermion approach is widely applicable and physically transparent: heavy electrons carry the strong correlations and dispersive electrons carry the topology. This simple picture unifies the dichotomous phenomena observed in TBG and points to connections between moiré and stoichiometric materials. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Moiré Fractional Chern Insulators II: First-principles Calculations and Continuum Models of Rhombohedral Graphene Superlattices

Authors: Jonah Herzog-Arbeitman, Yuzhi Wang, Jiaxuan Liu, Pok Man Tam, Ziyue Qi, Yujin Jia, Dmitri K. Efetov, Oskar Vafek, Nicolas Regnault, Hongming Weng, Quansheng Wu, B. Andrei Bernevig, Jiabin Yu

Abstract: The experimental discovery of fractional Chern insulators (FCIs) in rhombohedral pentalayer graphene twisted on hexagonal boron nitride (hBN) has preceded theoretical prediction. Supported by large-scale first principles relaxation calculations at the experimental twist angle of $0.77^\circ$, we obtain an accurate continuum model of $n=3,4,5,6,7$ layer rhombohedral graphene-hBN moiré systems. Focu… ▽ More The experimental discovery of fractional Chern insulators (FCIs) in rhombohedral pentalayer graphene twisted on hexagonal boron nitride (hBN) has preceded theoretical prediction. Supported by large-scale first principles relaxation calculations at the experimental twist angle of $0.77^\circ$, we obtain an accurate continuum model of $n=3,4,5,6,7$ layer rhombohedral graphene-hBN moiré systems. Focusing on the pentalayer case, we analytically explain the robust $|C|=0,5$ Chern numbers seen in the low-energy single-particle bands and their flattening with displacement field, making use of a minimal two-flavor continuum Hamiltonian derived from the full model. We then predict nonzero valley Chern numbers at the $νにゅー= -4,0$ insulators observed in experiment. Our analysis makes clear the importance of displacement field and the moiré potential in producing localized "heavy fermion" charge density in the top valence band, in addition to the nearly free conduction band. Lastly, we study doubly aligned devices as additional platforms for moiré FCIs with higher Chern number bands. △ Less

Submitted 21 November, 2023; originally announced November 2023.

Comments: Second paper in the moiré FCI series

Theory of correlated Chern insulators in twisted bilayer graphene

Authors: Xiaoyu Wang, Oskar Vafek

Abstract: Magic-angle twisted bilayer graphene is the best studied physical platform featuring moire potential induced narrow bands with non-trivial topology and strong electronic correlations. Despite their significance, the Chern insulating states observed at a finite magnetic field -- and extrapolating to a band filling, $s$, at zero field -- remain poorly understood. Unraveling their nature is among the… ▽ More Magic-angle twisted bilayer graphene is the best studied physical platform featuring moire potential induced narrow bands with non-trivial topology and strong electronic correlations. Despite their significance, the Chern insulating states observed at a finite magnetic field -- and extrapolating to a band filling, $s$, at zero field -- remain poorly understood. Unraveling their nature is among the most important open problems in the province of moiré materials. Here we present the first comprehensive study of interacting electrons in finite magnetic field while varying the electron density, twist angle and heterostrain. Within a panoply of correlated Chern phases emerging at a range of twist angles, we uncover a unified description for the ubiquitous sequence of states with the Chern number $t$ for $(s,t)=\pm (0,4), \pm(1,3),\pm(2,2)$ and $\pm(3,1)$. We also find correlated Chern insulators at unconventional sequences with $s+t\neq \pm 4$, as well as with fractional $s$, and elucidate their nature. △ Less

Submitted 27 March, 2024; v1 submitted 24 October, 2023; originally announced October 2023.

Comments: major revisions compared to the v2, including refined B-SCHF algorithm, and connections to B=0 physics of TBG

Journal ref: Phys. Rev. X 14, 021042 (2024)

Fractional Chern Insulators vs. Non-Magnetic States in Twisted Bilayer MoTe$_2$

Authors: Jiabin Yu, Jonah Herzog-Arbeitman, Minxuan Wang, Oskar Vafek, B. Andrei Bernevig, Nicolas Regnault

Abstract: Fractionally filled Chern bands with strong interactions may give rise to fractional Chern insulator (FCI) states, the zero-field analogue of the fractional quantum Hall effect. Recent experiments have demonstrated the existence of FCIs in twisted bilayer MoTe$_2$ without external magnetic fields -- most robust at $νにゅー=-2/3$ -- as well as Chern insulators (CIs) at $νにゅー=-1$. Although the appearance of… ▽ More Fractionally filled Chern bands with strong interactions may give rise to fractional Chern insulator (FCI) states, the zero-field analogue of the fractional quantum Hall effect. Recent experiments have demonstrated the existence of FCIs in twisted bilayer MoTe$_2$ without external magnetic fields -- most robust at $νにゅー=-2/3$ -- as well as Chern insulators (CIs) at $νにゅー=-1$. Although the appearance of both of these states is theoretically natural in an interacting topological system, experiments repeatedly observe nonmagnetic states (lacking FCIs) at $νにゅー=-1/3$ and $-4/3$, a puzzling result which has not been fully theoretically explained. In this work, we perform Hartree-Fock and exact diagonalization calculations to test whether the standard MoTe$_2$ moiré model with the (greatly varying) parameter values available in the literature can reproduce the non-magnetic states at $νにゅー=-1/3$ and $-4/3$ in unison with the FCI at $νにゅー=-2/3$ and CI state at $νにゅー= -1$. We focus on the experimentally relevant twist angles and, crucially, include remote bands. We find that the parameters proposed in [Wang et al. (2023)] can nearly capture the experimental phenomena at $νにゅー=-1/3,-2/3,-1,-4/3$ simultaneously, though the predicted ground states at $νにゅー=-1/3$ are still mostly fully-spin-polarized and a larger dielectric constant $εいぷしろん>10$ than is typical of hexagonal boron nitride (h-BN) substrate $εいぷしろん\sim 6$ is required. Our results show the importance of remote bands in identifying the competing magnetic orders and lay the groundwork for further study of the realistic phase diagram. △ Less

Submitted 25 September, 2023; originally announced September 2023.

Comments: 16+32 pages, 10+28 figures, 2+1 tables

Topological heavy fermions in magnetic field

Authors: Keshav Singh, Aaron Chew, Jonah Herzog-Arbeitman, B. Andrei Bernevig, Oskar Vafek

Abstract: The recently introduced topological heavy fermion model (THFM) provides a means for interpreting the low-energy electronic degrees of freedom of the magic angle twisted bilayer graphene as hybridization amidst highly dispersing topological conduction and weakly dispersing localized heavy fermions. In order to understand the Landau quantization of the ensuing electronic spectrum, a generalization o… ▽ More The recently introduced topological heavy fermion model (THFM) provides a means for interpreting the low-energy electronic degrees of freedom of the magic angle twisted bilayer graphene as hybridization amidst highly dispersing topological conduction and weakly dispersing localized heavy fermions. In order to understand the Landau quantization of the ensuing electronic spectrum, a generalization of THFM to include the magnetic field B is desired, but currently missing. Here we provide a systematic derivation of the THFM in B and solve the resulting model to obtain the interacting Hofstadter spectra for single particle charged excitations. While naive minimal substitution within THFM fails to correctly account for the total number of magnetic subbands within the narrow band i.e. its total Chern number, our method -- based on projecting the light and heavy fermions onto the irreducible representations of the magnetic translation group -- reproduces the correct total Chern number. Analytical results presented here offer an intuitive understanding of the nature of the (strongly interacting) Hofstadter bands. △ Less

Submitted 22 June, 2024; v1 submitted 14 May, 2023; originally announced May 2023.

Comments: 16+65 pages, 5+21 figures. The manuscript layout is restructured and the references are modified. The derivation of effective models is added to the supplementary notes

Journal ref: Nature Communications volume 15, 5257 (2024)

Revisiting Bloch electrons in magnetic field: Hofstadter physics via hybrid Wannier states

Authors: Xiaoyu Wang, Oskar Vafek

Abstract: We revisit the Hofstadter butterfly for a subset of topologically trivial Bloch bands arising from a continuum free electron Hamiltonian in a periodic lattice potential. We employ the recently developed procedure -- which was previously used to analyze the case of topologically non-trivial bands [\href{https://journals.aps.org/prb/abstract/10.1103/PhysRevB.106.L121111}{Phys. Rev. B \textbf{106}, L… ▽ More We revisit the Hofstadter butterfly for a subset of topologically trivial Bloch bands arising from a continuum free electron Hamiltonian in a periodic lattice potential. We employ the recently developed procedure -- which was previously used to analyze the case of topologically non-trivial bands [\href{https://journals.aps.org/prb/abstract/10.1103/PhysRevB.106.L121111}{Phys. Rev. B \textbf{106}, L121111 (2022)}] -- to construct the finite field Hilbert space from the zero-field hybrid Wannier basis states. Such states are Bloch extended along one direction and exponentially localized along the other. The method is illustrated for square and triangular lattice potentials and is shown to reproduce all the main features of the Hofstadter spectrum obtained from a numerically exact Landau level expansion method. In the regime when magnetic length is much longer than the spatial extent of the hybrid Wannier state in the localized direction we recover the well known Harper equation. Because the method applies to both topologically trivial and non-trivial bands, it provides an alternative and efficient approach to moiré materials in magnetic field. △ Less

Submitted 7 December, 2023; v1 submitted 28 March, 2023; originally announced March 2023.

Comments: added an appendix E comparing Peierls substitution method to hWS method; added journal ref

Journal ref: Phys. Rev. B 108, 245109 (2023)

Electronic anisotropy in magic-angle twisted trilayer graphene

Authors: Naiyuan J. Zhang, Yibang Wang, Kenji Watanabe, Takashi Taniguchi, Oskar Vafek, J. I. A. Li

Abstract: Due to its potential connection with nematicity, electronic anisotropy has been the subject of intense research effort on a wide variety of material platforms. The emergence of spatial anisotropy not only offers a characterization of material properties of metallic phases, which cannot be accessed via conventional transport techniques, but it also provides a unique window into the interplay betwee… ▽ More Due to its potential connection with nematicity, electronic anisotropy has been the subject of intense research effort on a wide variety of material platforms. The emergence of spatial anisotropy not only offers a characterization of material properties of metallic phases, which cannot be accessed via conventional transport techniques, but it also provides a unique window into the interplay between Coulomb interaction and broken symmetry underlying the electronic order. In this work, we utilize a new scheme of angle-resolved transport measurement (ARTM) to characterize electron anisotropy in magic-angle twisted trilayer graphene. By analyzing the dependence of spatial anisotropy on moiré band filling, temperature and twist angle, we establish the first experimental link between electron anisotropy and the cascade phenomenon, where Coulomb interaction drives a number of isospin transitions near commensurate band fillings. Furthermore, we report the coexistence between electron anisotropy and a novel electronic order that breaks both parity and time reversal symmetry. Combined, the link between electron anisotropy, cascade phenomenon and PT-symmetry breaking sheds new light onto the nature of electronic order in magic-angle graphene moiré systems. △ Less

Submitted 2 November, 2022; originally announced November 2022.

Comments: Main text 7 pages, 5 figures. Total 17 pages, 15 figures. arXiv admin note: text overlap with arXiv:2209.12964

Phase diagram of twisted bilayer graphene at filling factor $νにゅー=\pm3$

Authors: Fang Xie, Jian Kang, B. Andrei Bernevig, Oskar Vafek, Nicolas Regnault

Abstract: We study the correlated insulating phases of twisted bilayer graphene (TBG) in the absence of lattice strain at integer filling $νにゅー=\pm3$. Using the self-consistent Hartree-Fock method on a particle-hole symmetric model and allowing translation symmetry breaking terms, we obtain the phase diagram with respect to the ratio of $AA$ interlayer hopping $(w_0)$ and $AB$ interlayer hopping $(w_1)$. When… ▽ More We study the correlated insulating phases of twisted bilayer graphene (TBG) in the absence of lattice strain at integer filling $νにゅー=\pm3$. Using the self-consistent Hartree-Fock method on a particle-hole symmetric model and allowing translation symmetry breaking terms, we obtain the phase diagram with respect to the ratio of $AA$ interlayer hopping $(w_0)$ and $AB$ interlayer hopping $(w_1)$. When the interlayer hopping ratio is close to the chiral limit ($w_0/w_1 \lesssim 0.5$), a quantum anomalous Hall state with Chern number $νにゅー_c = \pm 1$ can be observed consistent with previous studies. Around the realistic value $w_0/w_1 \approx 0.8$, we find a spin and valley polarized, translation symmetry breaking, state with $C_{2z}T$ symmetry, a charge gap and a doubling of the moiré unit cell, dubbed the $C_{2z}T$ stripe phase. The real space total charge distribution of this $C_{2z}T$ stripe phase in the flat band limit does not have modulation between different moiré unit cells, although the charge density in each layer is modulated, and the translation symmetry is strongly broken. Other symmetries, including $C_{2z}$, $C_{2x}$ and particle-hole symmetry $P$, and the topology of the $C_{2z}T$ stripe phase are also discussed in detail. We observed braiding and annihilation of the Dirac nodes by continuously turning on the order parameter to its fully self-consistent value, and provide a detailed explanation of the mechanism for the charge gap opening despite preserving $C_{2z}T$ and valley $U(1)$ symmetries. In the transition region between the quantum anomalous Hall phase and the $C_{2z}T$ stripe phase, we find an additional competing state with comparable energy corresponding to a phase with a tripling of the moiré unit cell. △ Less

Submitted 27 February, 2023; v1 submitted 28 September, 2022; originally announced September 2022.

Comments: 39 pages, 33 figures

Journal ref: Phys. Rev. B 107, 075156 (2023)

Anisotropic resistivity tensor from disk geometry magneto-conductance

Authors: Oskar Vafek

Abstract: Magneto-transport measurements on two dimensional van der Waals heterostructures have recently shown signatures of uniaxial anisotropy. Such measurements are almost exclusively performed in a Hall bar geometry which makes it difficult to extract the full resistivity tensor. The goal of this paper is to theoretically analyze anisotropic magneto-conductance in a homogeneous disk geometry and to prov… ▽ More Magneto-transport measurements on two dimensional van der Waals heterostructures have recently shown signatures of uniaxial anisotropy. Such measurements are almost exclusively performed in a Hall bar geometry which makes it difficult to extract the full resistivity tensor. The goal of this paper is to theoretically analyze anisotropic magneto-conductance in a homogeneous disk geometry and to provide a closed form expression for the electrical potential anywhere on the disk if the current source and drain are located somewhere on the circumference. This expression can then be used to experimentally extract the full conductivity tensor, and by a simple inversion, the full resistivity tensor. △ Less

Submitted 27 November, 2023; v1 submitted 16 September, 2022; originally announced September 2022.

Comments: 7 pages, 2 figures (clarifications added)

Journal ref: Phys. Rev. Applied 20, 064008 (2023)

Unusual magnetotransport in twisted bilayer graphene from strain-induced open Fermi surfaces

Authors: Xiaoyu Wang, Joe Finney, Aaron L. Sharpe, Linsey K. Rodenbach, Connie L. Hsueh, Kenji Watanabe, Takashi Taniguchi, M. A. Kastner, Oskar Vafek, David Goldhaber-Gordon

Abstract: Anisotropic hopping in a toy Hofstadter model was recently invoked to explain a rich and surprising Landau spectrum measured in twisted bilayer graphene away from the magic angle. Suspecting that such anisotropy could arise from unintended uniaxial strain, we extend the Bistritzer-MacDonald model to include uniaxial heterostrain. We find that such strain strongly influences band structure, shiftin… ▽ More Anisotropic hopping in a toy Hofstadter model was recently invoked to explain a rich and surprising Landau spectrum measured in twisted bilayer graphene away from the magic angle. Suspecting that such anisotropy could arise from unintended uniaxial strain, we extend the Bistritzer-MacDonald model to include uniaxial heterostrain. We find that such strain strongly influences band structure, shifting the three otherwise-degenerate van Hove points to different energies. Coupled to a Boltzmann magnetotransport calculation, this reproduces previously-unexplained non-saturating $B^2$ magnetoresistance over broad ranges of density near filling $νにゅー=\pm 2$, and predicts subtler features that had not been noticed in the experimental data. In contrast to these distinctive signatures in longitudinal resistivity, the Hall coefficient is barely influenced by strain, to the extent that it still shows a single sign change on each side of the charge neutrality point -- surprisingly, this sign change no longer occurs at a van Hove point. The theory also predicts a marked rotation of the electrical transport principal axes as a function of filling even for fixed strain and for rigid bands. More careful examination of interaction-induced nematic order versus strain effects in twisted bilayer graphene could thus be in order. △ Less

Submitted 16 September, 2022; originally announced September 2022.

Comments: 11 pages main text (4 figures) + 8 pages supplementary material (11 figures)

Journal ref: PNAS 120 (34) e2307151120 (2023)

Triangular lattice Hubbard model physics at intermediate temperatures

Authors: Kyungmin Lee, Prakash Sharma, Oskar Vafek, Hitesh J. Changlani

Abstract: Moire systems offer an exciting playground to study many-body effects of strongly correlated electrons in regimes that are not easily accessible in conventional material settings. Motivated by a recent experiment on $\text{WSe}_2/\text{WS}_2$ moire bilayers [Y. Tang et al., Nature 579, 353-358 (2020)], which realizes a triangular superlattice with a small hopping (of approximately 10 Kelvin), with… ▽ More Moire systems offer an exciting playground to study many-body effects of strongly correlated electrons in regimes that are not easily accessible in conventional material settings. Motivated by a recent experiment on $\text{WSe}_2/\text{WS}_2$ moire bilayers [Y. Tang et al., Nature 579, 353-358 (2020)], which realizes a triangular superlattice with a small hopping (of approximately 10 Kelvin), with tunable density of holes, we explore the Hubbard model on the triangular lattice for intermediate temperatures $t \lesssim T < U$. Employing finite temperature Lanczos calculations, and closely following the fitting protocols used in the experiment, we recover the observed trends in the reported Curie-Weiss temperature $Θしーた$ with filling, using the reported interaction strength $U/t=20$. We focus on the large increase of $|Θしーた|$ on decreasing the density below half filling and the sign change of $Θしーた$ at higher fillings, which signals the onset of ferromagnetism. The increase in $|Θしーた|$ is also seen in the $t$-$J$ model (the low energy limit of the Hubbard model) in the intermediate temperature range, which we clarify is opposite to the trend in its high temperature limit. Differences between the low, intermediate and high temperature regimes are discussed. Our numerical calculations also capture the crossover between short-range antiferromagnetic to ferromagnetic order in the intermediate temperature regime, a result broadly consistent with the experimental findings. We find that this behavior is a finite-temperature remnant of the underlying zero temperature phase transition, which we explore with ground state density matrix renormalization group calculations. We provide evidence of ferromagnetism characterized by weak (but robust) correlations that explain the small $Θしーた$ seen in the experiment. △ Less

Submitted 5 June, 2023; v1 submitted 1 September, 2022; originally announced September 2022.

Comments: 11 pages, 12 figures. v2 closely matches the published version. Sections expanded and supplemental material converted to appendices. The first two authors contributed equally

Journal ref: Phys. Rev. B 107, 235105 (2023)

Twisted-bilayer FeSe and the Fe-based superlattices

Authors: P. Myles Eugenio, Oskar Vafek

Abstract: We derive BM-like continuum models for the bands of superlattice heterostructures formed out of Fe-chalcogenide monolayers: (${\bf\text I}$) a single monolayer experiencing an external periodic potential, and (${\bf\text II}$) twisted bilayers with long-range moire tunneling. A symmetry derivation for the inter-layer moire tunnelling is provided for both the $Γがんま$ and $M$ high-symmetry points. In th… ▽ More We derive BM-like continuum models for the bands of superlattice heterostructures formed out of Fe-chalcogenide monolayers: (${\bf\text I}$) a single monolayer experiencing an external periodic potential, and (${\bf\text II}$) twisted bilayers with long-range moire tunneling. A symmetry derivation for the inter-layer moire tunnelling is provided for both the $Γがんま$ and $M$ high-symmetry points. In this paper, we focus on moire bands formed from hole-band maxima centered on $Γがんま$, and show the possibility of moire bands with $C=0$ or $\pm 1$ topological quantum numbers without breaking time-reversal symmetry. In the $C=0$ region for $θしーた\rightarrow 0$ (and similarly in the limit of large superlattice period for ${\bf\text I}$), the system becomes a square lattice of 2D harmonic oscillators. We fit our model to FeSe and argue that it is a viable platform for the simulation of the square Hubbard model with tunable interaction strength. △ Less

Submitted 22 May, 2023; v1 submitted 23 August, 2022; originally announced August 2022.

Comments: 35 pages; 18 figs; Changed name from "Minimal continuum models in Fe-based superlattices"

Journal ref: SciPost Phys. 15, 081 (2023)

Pseudo-magnetic fields, particle-hole asymmetry, and microscopic effective continuum Hamitonians of twisted bilayer graphene

Authors: Jian Kang, Oskar Vafek

Abstract: Using the method developed in the companion paper, we construct the effective continuum theories for two different microscopic tight binding models of the twisted bilayer graphene at the twist angle of $1.05^\circ$, one Slater-Koster based and the other ab-initio Wannier based. The energy spectra obtained from the continuum theory -- either for rigid twist or including lattice relaxation -- are fo… ▽ More Using the method developed in the companion paper, we construct the effective continuum theories for two different microscopic tight binding models of the twisted bilayer graphene at the twist angle of $1.05^\circ$, one Slater-Koster based and the other ab-initio Wannier based. The energy spectra obtained from the continuum theory -- either for rigid twist or including lattice relaxation -- are found to be in nearly perfect agreement with the spectra from the tight binding models when the gradient expansion is carried out to second order, demonstrating the validity of the method. We also analyze the properties of the Bloch states of the resulting narrow bands, finding non-negligible particle-hole symmetry breaking near the $Γがんま$ point in our continuum theory constructed for the ab-initio based microscopic model due to a term in the continuum theory that was previously overlooked. This reveals the difference with all existing continuum models where the particle-hole symmetry of the narrow band Hilbert space is nearly perfect. △ Less

Submitted 19 February, 2023; v1 submitted 11 August, 2022; originally announced August 2022.

Comments: 20 pages, 13 figures, and 18 pages appendix. Typos are corrected

Journal ref: Phys. Rev. B 107, 075408 (2023)

Continuum effective Hamiltonian for graphene bilayers for an arbitrary smooth lattice deformation from microscopic theories

Authors: Oskar Vafek, Jian Kang

Abstract: We provide a systematic real space derivation of the continuum Hamiltonian for a graphene bilayer starting from a microscopic lattice theory, allowing for an arbitrary inhomogeneous smooth lattice deformation, including a twist. Two different microscopic models are analyzed: first, a Slater-Koster like model and second, ab-initio derived model. We envision that our effective Hamiltonian can be use… ▽ More We provide a systematic real space derivation of the continuum Hamiltonian for a graphene bilayer starting from a microscopic lattice theory, allowing for an arbitrary inhomogeneous smooth lattice deformation, including a twist. Two different microscopic models are analyzed: first, a Slater-Koster like model and second, ab-initio derived model. We envision that our effective Hamiltonian can be used in conjunction with an experimentally determined atomic lattice deformation in twisted bilayer graphene in a specific device to predict and compare the electronic spectra with scanning tunneling spectroscopy measurements. As a byproduct, our approach provides electron-phonon couplings in the continuum Hamiltonian from microscopic models for any bilayer stacking. In the companion paper we analyze in detail the continuum models for relaxed atomic configurations of magic angle twisted bilayer graphene. △ Less

Submitted 17 February, 2023; v1 submitted 11 August, 2022; originally announced August 2022.

Comments: 9 pages, 1 figure, typos are corrected

Journal ref: Phys. Rev. B 107, 075123 (2023)

Trions in Twisted Bilayer Graphene

Authors: Frank Schindler, Oskar Vafek, B. Andrei Bernevig

Abstract: The strong coupling phase diagram of magic angle twisted bilayer graphene (TBG) predicts a series of exact one particle charge $\pm 1$ gapped excitations on top of the integer-filled ferromagnetic ground-states. Finite-size exact diagonalization studies showed that these are the lowest charge $\pm 1$ excitations in the system (for $10$nm screening length), with the exception of charge $+1$ at fill… ▽ More The strong coupling phase diagram of magic angle twisted bilayer graphene (TBG) predicts a series of exact one particle charge $\pm 1$ gapped excitations on top of the integer-filled ferromagnetic ground-states. Finite-size exact diagonalization studies showed that these are the lowest charge $\pm 1$ excitations in the system (for $10$nm screening length), with the exception of charge $+1$ at filling $-1$ in the chiral limit. In the current paper we show that this "trion bound state", a $3$-particle, charge $1$ excitation of the insulating ferromagnetic ground-state of the projected Hamiltonian of TBG is the lowest charge $+1$ overall excitation at $νにゅー=-1$, and also for some large ($\approx 20$nm) screening lengths at $νにゅー=-2$ in the chiral limit and with very small binding energy. At other fillings, we show that trion bound states do exist, but only for momentum ranges that do not cover the entire moiré Brillouin zone. The trion bound states (at different momenta) exist for finite parameter range $w_0/w_1$ but they all disappear in the continuum far below the realistic values of $w_0/w_1= 0.8$. We find the conditions for the existence of the trion bound state, a good variational wavefunction for it, and investigate its behavior for different screening lengths, at all integer fillings, on both the electron and hole sides. △ Less

Submitted 22 April, 2022; v1 submitted 23 December, 2021; originally announced December 2021.

Comments: 30 pages, 19 figures

Journal ref: Phys. Rev. B 105, 155135 (2022)

Narrow bands in magnetic field and strong-coupling Hofstadter spectra

Authors: Xiaoyu Wang, Oskar Vafek

Abstract: We develop a new, efficient, and general method to determine the Hofstadter spectrum of isolated narrow bands. The method works for topological as well as for trivial narrow bands by projecting the zero $\mathbf{B}$-field hybrid Wannier states -- which are localized in one direction and Bloch extended in another direction -- onto a representation of the magnetic translation group in the Landau gau… ▽ More We develop a new, efficient, and general method to determine the Hofstadter spectrum of isolated narrow bands. The method works for topological as well as for trivial narrow bands by projecting the zero $\mathbf{B}$-field hybrid Wannier states -- which are localized in one direction and Bloch extended in another direction -- onto a representation of the magnetic translation group in the Landau gauge. We then apply this method to find, for the first time, the Hofstadter spectrum for the exact single particle charged excitations in the strong coupling limit of the magic angle twisted bilayer graphene at the charge neutrality point and at $|νにゅー|=2$ down to low magnetic fields when the flux through the moiré unit cell is only $\sim 1/25$ of the electronic flux quantum i.e. $\sim 1$T at the first magic angle. The resulting spectra provide a means to investigate Landau quantization of the quasiparticles even if their dispersion is interaction induced. △ Less

Submitted 11 January, 2022; v1 submitted 15 December, 2021; originally announced December 2021.

Comments: Version 2 with 17 pages and 10 figures. Strong coupling excitation spectra (heavy hole and light electron) at nu=2, and supplementary material added

Spectroscopy of Twisted Bilayer Graphene Correlated Insulators

Authors: Dumitru Călugăru, Nicolas Regnault, Myungchul Oh, Kevin P. Nuckolls, Dillon Wong, Ryan L. Lee, Ali Yazdani, Oskar Vafek, B. Andrei Bernevig

Abstract: We analytically compute the scanning tunneling microscopy (STM) signatures of integer-filled correlated ground states of the magic angle twisted bilayer graphene (TBG) narrow bands. After experimentally validating the strong-coupling approach at $\pm 4$ electrons/moiré unit cell, we consider the spatial features of the STM signal for 14 different many-body correlated states and assess the possibil… ▽ More We analytically compute the scanning tunneling microscopy (STM) signatures of integer-filled correlated ground states of the magic angle twisted bilayer graphene (TBG) narrow bands. After experimentally validating the strong-coupling approach at $\pm 4$ electrons/moiré unit cell, we consider the spatial features of the STM signal for 14 different many-body correlated states and assess the possibility of Kekulé distortion (KD) emerging at the graphene lattice scale. Remarkably, we find that coupling the two opposite graphene valleys in the intervalley-coherent (IVC) TBG insulators does not always result in KD. As an example, we show that the Kramers IVC state and its nonchiral $\mathrm{U} \left( 4 \right)$ rotations do not exhibit any KD, while the time-reversal-symmetric IVC state does. Our results, obtained over a large range of energies and model parameters, show that the STM signal and Chern number of a state can be used to uniquely determine the nature of the TBG ground state. △ Less

Submitted 15 September, 2022; v1 submitted 28 October, 2021; originally announced October 2021.

Comments: 8+80 pages, 2+65 figures. New version matches the published version

Journal ref: Phys. Rev. Lett. 129, 117602, 8 September 2022

Lattice model for the Coulomb interacting chiral limit of the magic angle twisted bilayer graphene: symmetries, obstructions and excitations

Authors: Oskar Vafek, Jian Kang

Abstract: We revisit the localized Wannier state description of the twisted bilayer graphene, focusing on the chiral limit. We provide a simple method for constructing such 2D exponentially localized -- yet valley polarized -- Wannier states, centered on the sites of the honeycomb lattice, paying particular attention to maintaining all the unobstructed symmetries. This includes the unitary particle-hole sym… ▽ More We revisit the localized Wannier state description of the twisted bilayer graphene, focusing on the chiral limit. We provide a simple method for constructing such 2D exponentially localized -- yet valley polarized -- Wannier states, centered on the sites of the honeycomb lattice, paying particular attention to maintaining all the unobstructed symmetries. This includes the unitary particle-hole symmetry, and the combination of $C_2\mathcal{T}$ and the chiral particle-hole symmetry. The $C_2\mathcal{T}$ symmetry alone remains topologically obstructed and is not represented in a simple site-to-site fashion. We also analyze the gap and the dispersion of single particle and single hole excitations above a strong coupling ground state at integer fillings, which we find to be dominated by the on-site and the nearest neighbor terms of a triangular lattice hopping model, with a minimum at the center of the moire Brillouin zone. Finally, we use the insight gained from this real space description to understand the dependence of the gap and the effective mass on the range of the screened Coulomb interaction. △ Less

Submitted 10 June, 2021; originally announced June 2021.

Comments: 19 pages; (14p main text + 1.2p references + 4p supplementary) 7 figures

Journal ref: Phys. Rev. B 104, 075143 (2021)

Cascades between light and heavy fermions in the normal state of magic angle twisted bilayer graphene

Authors: Jian Kang, B. Andrei Bernevig, Oskar Vafek

Abstract: We present a framework for understanding the recently observed cascade transitions and the Landau level degeneracies at every integer filling of twisted bilayer graphene. The Coulomb interaction projected onto narrow bands causes the charged excitations at an integer filling to disperse, forming new bands. If the excitation moves the filling away from the charge neutrality point, then it has a ban… ▽ More We present a framework for understanding the recently observed cascade transitions and the Landau level degeneracies at every integer filling of twisted bilayer graphene. The Coulomb interaction projected onto narrow bands causes the charged excitations at an integer filling to disperse, forming new bands. If the excitation moves the filling away from the charge neutrality point, then it has a band minimum at the moire Brillouin zone center with a small mass that compares well with the experiment; if towards the charge neutrality point, then it has a much larger mass and a higher degeneracy. At a non-zero density away from an integer filling the excitations interact. The system on the small mass side has a large bandwidth and forms a Fermi liquid. On the large mass side the bandwidth is narrow, the compressibility is negative and the Fermi liquid is likely unstable. This explains the observed sawtooth features in compressibility, the Landau fans pointing away from charge neutrality as well as their degeneracies. By providing a description of the charge itineracy in the normal state this framework sets the stage for superconductivity at lower temperatures. △ Less

Submitted 2 April, 2021; originally announced April 2021.

Comments: 18 pages; (4.5p main text + 2.5p references + 11p supplementary) 2+4 figures

Realization of Topological Mott Insulator in a Twisted Bilayer Graphene Lattice Model

Authors: Bin-Bin Chen, Yuan Da Liao, Ziyu Chen, Oskar Vafek, Jian Kang, Wei Li, Zi Yang Meng

Abstract: Magic-angle twisted bilayer graphene has recently become a thriving material platform realizing correlated electron phenomena taking place within its topological flat bands. Several numerical and analytical methods have been applied to understand the correlated phases therein, revealing some similarity with the quantum Hall physics. In this work, we provide a Mott-Hubbard perspective for the TBG s… ▽ More Magic-angle twisted bilayer graphene has recently become a thriving material platform realizing correlated electron phenomena taking place within its topological flat bands. Several numerical and analytical methods have been applied to understand the correlated phases therein, revealing some similarity with the quantum Hall physics. In this work, we provide a Mott-Hubbard perspective for the TBG system. Employing the large-scale density matrix renormalization group on the lattice model containing the projected Coulomb interactions only, we identify a first-order quantum phase transition between the insulating stripe phase and the quantum anomalous Hall state with the Chern number of $\pm 1$. Our results not only shed light on the mechanism of the quantum anomalous Hall state discovered at three-quarters filling, but also provide an example of the topological Mott insulator, i.e., the quantum anomalous Hall state in the strong coupling limit. △ Less

Submitted 16 September, 2021; v1 submitted 15 November, 2020; originally announced November 2020.

Comments: 6+7 pages, 3+7 figures

Journal ref: Nature Communications 12, 5480 (2021)

Towards the hidden symmetry in Coulomb interacting twisted bilayer graphene: renormalization group approach

Authors: Oskar Vafek, Jian Kang

Abstract: We develop a two stage renormalization group which connects the continuum Hamiltonian for twisted bilayer graphene at length scales shorter than the moire superlattice period to the Hamiltonian for the active narrow bands only which is valid at distances much longer than the moire period. In the first stage, the Coulomb interaction renormalizes the Fermi velocity and the interlayer tunnelings in s… ▽ More We develop a two stage renormalization group which connects the continuum Hamiltonian for twisted bilayer graphene at length scales shorter than the moire superlattice period to the Hamiltonian for the active narrow bands only which is valid at distances much longer than the moire period. In the first stage, the Coulomb interaction renormalizes the Fermi velocity and the interlayer tunnelings in such a way as to suppress the ratio of the same sublattice to opposite sublatice tunneling, hence approaching the so-called chiral limit. In the second stage, the interlayer tunneling is treated non-perturbatively. Via a progressive numerical elimination of remote bands the relative strength of the one-particle-like dispersion and the interactions within the active narrow band Hamiltonian is determined, thus quantifying the residual correlations and justifying the strong coupling approach in the final step. We also calculate exactly the exciton energy spectrum from the Coulomb interactions projected onto the renormalized narrow bands. The resulting softening of the collective modes marks the propinquity of the enlarged ("hidden") $U(4)\times U(4)$ symmetry in the magic angle twisted bilayer graphene. △ Less

Submitted 11 November, 2020; v1 submitted 20 September, 2020; originally announced September 2020.

Comments: 11 pages; (4p main text + 2.5p references + 4p supplementary) 3+3 figures; added a new calculation of the RG flow of the particle-hole asymmetry term and an outline of the numerical recipe for stage 2 RG

Journal ref: PRL 125, 257602 (2020)

Tuning electron correlation in magic-angle twisted bilayer graphene using Coulomb screening

Authors: Xiaoxue Liu, Zhi Wang, K. Watanabe, T. Taniguchi, Oskar Vafek, J. I. A. Li

Abstract: The ability to control the strength of interaction is essential for studying quantum phenomena emerging from a system of correlated fermions. For example, the isotope effect illustrates the effect of electron-phonon coupling on superconductivity, providing an important experimental support for the BCS theory. In this work, we report a new device geometry where the magic-angle twisted bilayer graph… ▽ More The ability to control the strength of interaction is essential for studying quantum phenomena emerging from a system of correlated fermions. For example, the isotope effect illustrates the effect of electron-phonon coupling on superconductivity, providing an important experimental support for the BCS theory. In this work, we report a new device geometry where the magic-angle twisted bilayer graphene (tBLG) is placed in close proximity to a Bernal bilayer graphene (BLG) separated by a 3 nm thick barrier. Using charge screening from the Bernal bilayer, the strength of electron-electron Coulomb interaction within the twisted bilayer can be continuously tuned. Transport measurements show that tuning Coulomb screening has opposite effect on the insulating and superconducting states: as Coulomb interaction is weakened by screening, the insulating states become less robust, whereas the stability of superconductivity is enhanced. Out results demonstrate the ability to directly probe the role of Coulomb interaction in magic-angle twisted bilayer graphene. Most importantly, the effect of Coulomb screening points toward electron-phonon coupling as the dominant mechanism for Cooper pair formation, and therefore superconductivity, in magic-angle twisted bilayer graphene. △ Less

Submitted 8 July, 2020; v1 submitted 24 March, 2020; originally announced March 2020.

Comments: 6 pages, 4 figures plus supplementary information. Updated with results from new measurements performed with low-pass filtering

Non-Abelian Dirac node braiding and near-degeneracy of correlated phases at odd integer filling in magic angle twisted bilayer graphene

Authors: Jian Kang, Oskar Vafek

Abstract: We use the DMRG to study the correlated electron states favored by the Coulomb interaction projected onto the narrow bands of twisted bilayer graphene within a spinless one-valley model. The Hilbert space of the narrow bands is constructed from a pair of hybrid Wannier states with opposite Chern numbers. Depending on the parameters in the BM model, the DMRG in this basis determines the ground stat… ▽ More We use the DMRG to study the correlated electron states favored by the Coulomb interaction projected onto the narrow bands of twisted bilayer graphene within a spinless one-valley model. The Hilbert space of the narrow bands is constructed from a pair of hybrid Wannier states with opposite Chern numbers. Depending on the parameters in the BM model, the DMRG in this basis determines the ground state at one particle per unit cell to be either QAH state or a state with no Hall effect which is nearly a product state. Based on this form, we then apply the variational method to study their competition, thus identifying three states: the QAH, a gapless $C_2T$ symmetric nematic, and a gapped $C_2T$ symmetric stripe. All three states are nearly degenerate at the realistic parameters of the BM model. The single particle spectrum of the nematic contains either a quadratic node or two close Dirac nodes near $Γがんま$. Motivated by the Landau level degeneracy found in this state, we propose it to be the state observed at the charge neutrality point once spin and valley degeneracies are restored. The optimal period for the $C_2T$ stripe state is found to be $2$ unit cells. In addition, using the fact that the topological charge of the nodes in the $C_2T$ nematic phase is no longer described simply by their winding numbers once the translation symmetry is broken, but rather by certain elements of a non-Abelian group that was recently pointed out, we identify the mechanism of the gap opening within the $C_2T$ stripe state. Although the nodes at the Fermi energy are locally stable, they can be annihilated after braiding with other nodes connecting them to adjacent (folded) bands. Therefore, if the translation symmetry is broken, the gap at one particle per unit cell can open even if the system preserves the $C_2T$ and valley $U(1)$ symmetries, and the gap to remote bands remains open. △ Less

Submitted 13 July, 2020; v1 submitted 24 February, 2020; originally announced February 2020.

Comments: 33 pages and 20 figures

Journal ref: Phys. Rev. B 102, 035161 (2020)

A diagnosis of explicit symmetry breaking in the tight-binding constructions for symmetry-protected topological systems

Authors: Xiaoyu Wang, Oskar Vafek

Abstract: It has been well established that for symmetry protected topological systems, the non-trivial topology prevents a real space representation using exponentially localized Wannier wavefunctions (WFs) in all directions, unless the protecting symmetry is sacrificed as an on-site transformation. This makes it challenging to determine the symmetry of various physical observables represented using such W… ▽ More It has been well established that for symmetry protected topological systems, the non-trivial topology prevents a real space representation using exponentially localized Wannier wavefunctions (WFs) in all directions, unless the protecting symmetry is sacrificed as an on-site transformation. This makes it challenging to determine the symmetry of various physical observables represented using such WFs. In this work, we propose a practical method for overcoming such challenges using the Kane-Mele model as a concrete example. We present a systematic procedure for diagnosing the symmetry of any observables, as well as a method for constructing symmetric operators up to arbitrary truncation accuracy. △ Less

Submitted 20 August, 2020; v1 submitted 5 February, 2020; originally announced February 2020.

Comments: 6 pages + 5 figures

Journal ref: Phys. Rev. B 102, 075142 (2020)

Strong coupling phases of partially filled twisted bilayer graphene narrow bands

Authors: Jian Kang, Oskar Vafek

Abstract: We identify states favored by Coulomb interactions projected onto the Wannier basis of the four narrow bands of the "magic angle" twisted bilayer graphene. At the filling of two electrons/holes per moire unit cell, such interactions favor an insulating SU(4) ferromagnet. The kinetic terms select the ground state in which the two valleys with opposite spins are equally mixed, with vanishing magneti… ▽ More We identify states favored by Coulomb interactions projected onto the Wannier basis of the four narrow bands of the "magic angle" twisted bilayer graphene. At the filling of two electrons/holes per moire unit cell, such interactions favor an insulating SU(4) ferromagnet. The kinetic terms select the ground state in which the two valleys with opposite spins are equally mixed, with vanishing magnetic moment per particle. We also find extended excited states, the gap to which decreases in magnetic field. An insulating stripe ferromagnetic phase is favored at one electron/hole per unit cell. △ Less

Submitted 19 October, 2018; originally announced October 2018.

Comments: 5+10 pages

Journal ref: Phys. Rev. Lett. 122, 246401 (2019)

Symmetry, maximally localized Wannier states, and low energy model for the twisted bilayer graphene narrow bands

Authors: Jian Kang, Oskar Vafek

Abstract: We build symmetry adapted maximally localized Wannier states, and construct the low energy tight binding model for the four narrow bands of the twisted bilayer graphene. We do so when the twist angle is commensurate, near the `magic' value, and the narrow bands are separated from the rest of the bands by energy gaps. On each layer and sublattice, every Wannier state has three peaks near the triang… ▽ More We build symmetry adapted maximally localized Wannier states, and construct the low energy tight binding model for the four narrow bands of the twisted bilayer graphene. We do so when the twist angle is commensurate, near the `magic' value, and the narrow bands are separated from the rest of the bands by energy gaps. On each layer and sublattice, every Wannier state has three peaks near the triangular Moire lattice sites. However, each Wannier state is localized and centered around a site of the honeycomb lattice that is dual to the triangular Moire lattice. Space group and the time reversal symmetries are realized locally. The corresponding tight binding model provides a starting point for studying the correlated many-body phases. △ Less

Submitted 19 May, 2018; v1 submitted 13 May, 2018; originally announced May 2018.

Comments: 9 pages; typos corrected, discussion of the relation to Refs.[5] and [6] modified

Journal ref: Phys. Rev. X 8, 031088 (2018)

Quantum Anomalous Hall Insulator Stabilized By Competing Interactions

Authors: Shouvik Sur, Shou-Shu Gong, Kun Yang, Oskar Vafek

Abstract: We study the quantum phases driven by interaction in a semimetal with a quadratic band touching at the Fermi level. By combining the density matrix renormalization group (DMRG), analytical power expanded Gibbs potential method, and the weak coupling renormalization group, we study a spinless fermion system on a checkerboard lattice at half-filling, which has a quadratic band touching in the absenc… ▽ More We study the quantum phases driven by interaction in a semimetal with a quadratic band touching at the Fermi level. By combining the density matrix renormalization group (DMRG), analytical power expanded Gibbs potential method, and the weak coupling renormalization group, we study a spinless fermion system on a checkerboard lattice at half-filling, which has a quadratic band touching in the absence of interaction. In the presence of strong nearest-neighbor ($V_1$) and next-nearest-neighbor ($V_2$) interactions, we identify a site nematic insulator phase, a stripe insulator phase, and a phase separation region, in agreement with the phase diagram obtained analytically in the strong coupling limit (i.e. in the absence of fermion hopping). In the intermediate interaction regime, we establish a quantum anomalous Hall phase in the DMRG as evidenced by the spontaneous time-reversal symmetry breaking and the appearance of a quantized Chern number $C = 1$. For weak interaction, we utilize the power expanded Gibbs potential method that treats $V_1$ and $V_2$ on equal footing, as well as the weak coupling renormalization group. Our analytical results reveal that not only the repulsive $V_1$ interaction, but also the $V_2$ interaction (both repulsive and attractive), can drive the quantum anomalous Hall phase. We also determine the phase boundary in the $V_1$-$V_2$ plane that separates the semimetal from the quantum anomalous Hall state. Finally, we show that the nematic semimetal, which was proposed for $|V_2| \gg V_1$ at weak coupling in a previous study, is absent, and the quantum anomalous Hall state is the only weak coupling instability of the spinless quadratic band touching semimetal. △ Less

Submitted 28 September, 2018; v1 submitted 3 May, 2018; originally announced May 2018.

Comments: 25 pages, 19 figures

Journal ref: Phys. Rev. B 98, 125144 (2018)

Classification of Symmetry Derived Pairing at M Point in FeSe

Authors: P. Myles Eugenio, Oskar Vafek

Abstract: Using the constraints imposed by the crystalline symmetry of FeSe and the experimentally observed phenomenology, we analyze the possible pairing symmetry of the superconducting order parameter focusing on intercalated and monolayer FeSe compounds. Such analysis leads to three possible pairing symmetry states -- s-wave, d-wave, and helical p-wave. Despite the differences in the pairing symmetry, ea… ▽ More Using the constraints imposed by the crystalline symmetry of FeSe and the experimentally observed phenomenology, we analyze the possible pairing symmetry of the superconducting order parameter focusing on intercalated and monolayer FeSe compounds. Such analysis leads to three possible pairing symmetry states -- s-wave, d-wave, and helical p-wave. Despite the differences in the pairing symmetry, each of these states is fully gapped with gap minimum centered above the normal state Fermi surface, in agreement with photoemission data of Y. Zhang et al. The analysis provides additional insights into the possible pairing mechanism for each of these states, highlighting the detrimental role of the renormalized repulsive intra-orbital Hubbard $U$ and inter-orbital $U'$ , and the beneficial role of the pair hopping $J'$ and the Hunds $J$ terms, as well as the spin-orbit coupling in the effective low-energy Hamiltonian. △ Less

Submitted 2 July, 2018; v1 submitted 21 April, 2018; originally announced April 2018.

Comments: 17 pages, 13 figures. Minor corrections in order to match PRB publication, this includes: removed erroneous citation, added citation, minor corrections to wordings, and added paragraph. Further, the author corrected the pairing symmetry of the inter-band spin-triplet (w/ d-vector out of plane) pairs in Table III. This was a labeling mistake and has no impact on the rest of the paper

Journal ref: Phys. Rev. B 98, 014503 (2018)

Hund interaction, spin-orbit coupling and the mechanism of superconductivity in strongly hole-doped iron pnictides

Authors: Oskar Vafek, Andrey V. Chubukov

Abstract: We present a novel mechanism of s-wave pairing in Fe-based superconductors. The mechanism involves holes near dxz/dyz pockets only and is applicable primarily to strongly hole doped materials. We argue that as long as the renormalized Hund's coupling J exceeds the renormalized inter-orbital Hubbard repulsion U', any finite spin-orbit coupling gives rise to s-wave superconductivity. This holds even… ▽ More We present a novel mechanism of s-wave pairing in Fe-based superconductors. The mechanism involves holes near dxz/dyz pockets only and is applicable primarily to strongly hole doped materials. We argue that as long as the renormalized Hund's coupling J exceeds the renormalized inter-orbital Hubbard repulsion U', any finite spin-orbit coupling gives rise to s-wave superconductivity. This holds even at weak coupling and regardless of the strength of the intra-orbital Hubbard repulsion U. The transition temperature grows as the hole density decreases. The pairing gaps are four-fold symmetric, but anisotropic, with the possibility of eight accidental nodes along the larger pocket. The resulting state is consistent with the experiments on KFe2As2. △ Less

Submitted 17 November, 2016; originally announced November 2016.

Comments: 5 pages (main text) + 5 pages (supplementary material), 3 figures

Journal ref: Phys. Rev. Lett. 118, 087003 (2017)

Entanglement of Exact Excited Eigenstates of the Hubbard Model in Arbitrary Dimension

Authors: Oskar Vafek, Nicolas Regnault, B. Andrei Bernevig

Abstract: We compute exactly the von Neumann entanglement entropy of the eta-pairing states - a large set of exact excited eigenstates of the Hubbard Hamiltonian. For the singlet eta-pairing states the entropy scales with the logarithm of the spatial dimension of the (smaller) partition. For the eta-pairing states with finite spin magnetization density, the leading term can scale as the volume or as the are… ▽ More We compute exactly the von Neumann entanglement entropy of the eta-pairing states - a large set of exact excited eigenstates of the Hubbard Hamiltonian. For the singlet eta-pairing states the entropy scales with the logarithm of the spatial dimension of the (smaller) partition. For the eta-pairing states with finite spin magnetization density, the leading term can scale as the volume or as the area-times-log, depending on the momentum space occupation of the Fermions with flipped spins. We also compute the corrections to the leading scaling. In order to study the eigenstate thermalization hypothesis (ETH), we also compute the entanglement Renyi entropies of such states and compare them with the corresponding entropies of thermal density matrix in various ensembles. Such states, which we find violate strong ETH, may provide a useful platform for a detailed study of the time-dependence of the onset of thermalization due to perturbations which violate the total pseudospin conservation. △ Less

Submitted 3 December, 2017; v1 submitted 23 August, 2016; originally announced August 2016.

Comments: 5 pages (main text) + 10 pages (suppl. info); 1 figure and 1 table; References added and minor revisions

Journal ref: SciPost Phys. 3, 043 (2017)

Displacement and annihilation of Dirac gap-nodes in d-wave iron-based superconductors

Authors: Andrey V. Chubukov, Oskar Vafek, Rafael M. Fernandes

Abstract: Several experimental and theoretical arguments have been made in favor of a $d-$wave symmetry for the superconducting state in some Fe-based materials. It is a common belief that a $d-$wave gap in the Fe-based superconductors must have nodes on the Fermi surfaces centered at the $Γがんま$ point of the Brillouin zone. Here we show that, while this is the case for a single Fermi surface made out of a sing… ▽ More Several experimental and theoretical arguments have been made in favor of a $d-$wave symmetry for the superconducting state in some Fe-based materials. It is a common belief that a $d-$wave gap in the Fe-based superconductors must have nodes on the Fermi surfaces centered at the $Γがんま$ point of the Brillouin zone. Here we show that, while this is the case for a single Fermi surface made out of a single orbital, the situation is more complex if there is an even number of Fermi surfaces made out of different orbitals. In particular, we show that for the two $Γがんま$-centered hole Fermi surfaces made out of $d_{xz}$ and $d_{yz}$ orbitals, the nodal points still exist near $T_{c}$ along the symmetry-imposed directions, but are are displaced to momenta between the two Fermi surfaces. If the two hole pockets are close enough, pairs of nodal points can merge and annihilate at some $T<T_{c}$, making the $d-$wave state completely nodeless. These results imply that photoemission evidence for a nodeless gap on the $d_{xz}/d_{yz}$ Fermi surfaces of KFe$_{2}$As$_{2}$ does not rule out $d-$wave gap symmetry in this material, while a nodeless gap observed on the $d_{xy}$ pocket in K$_{x}$Fe$_{2-y}$Se$_{2}$ is truly inconsistent with the $d-$wave gap symmetry. △ Less

Submitted 7 September, 2016; v1 submitted 20 August, 2016; originally announced August 2016.

Comments: 7 pages, 4 figures, submitted to PRB, references added

Journal ref: Phys. Rev. B 94, 174518 (2016)

Intrinsic thermal Hall conductivity in the mixed state of d-wave superconductors: from wavepacket dynamics to scaling

Authors: Oskar Vafek

Abstract: Recent numerical calculation of the intrinsic thermal Hall conductivity of nodal d-wave superconductors in the mixed state revealed a rapid increase of this quantity above an onset temperature. Interestingly, this defines a measurable energy scale in an otherwise gapless state. Using the mathematics of magnetic coherent states, in this paper such energy scale is related to a dynamical process asso… ▽ More Recent numerical calculation of the intrinsic thermal Hall conductivity of nodal d-wave superconductors in the mixed state revealed a rapid increase of this quantity above an onset temperature. Interestingly, this defines a measurable energy scale in an otherwise gapless state. Using the mathematics of magnetic coherent states, in this paper such energy scale is related to a dynamical process associated with the Andreev scattering of an electron wavepacket moving along the constant energy contours in the momentum space. This energy scale is then used to obtain an improved scaling collapse of numerically calculated thermal Hall conductivity in a tight-binding model as a function of temperature, magnetic field and the d-wave pairing amplitude at various band fillings. The results indicate that the mentioned onset temperature is associated with the ability of the quasiparticle wavepacket to complete its semiclassical orbit before it is appreciably scattered by the superconducting condensate. △ Less

Submitted 9 August, 2015; originally announced August 2015.

Comments: 11 pages, 2 figures

Journal ref: Phys. Rev. B 92, 174508 (2015)

Majorana bands, Berry curvature, and thermal Hall conductivity in the vortex state of a chiral p-wave superconductor

Authors: James M. Murray, Oskar Vafek

Abstract: Majorana quasiparticles localized in vortex cores of a chiral p-wave superconductor hybridize with one another to form bands in a vortex lattice. We begin by solving a fully microscopic theory describing all quasiparticle bands in a chiral p-wave superconductor in magnetic field, then use this solution to build localized Wannier wavefunctions corresponding to Majorana quasiparticles. A low-energy… ▽ More Majorana quasiparticles localized in vortex cores of a chiral p-wave superconductor hybridize with one another to form bands in a vortex lattice. We begin by solving a fully microscopic theory describing all quasiparticle bands in a chiral p-wave superconductor in magnetic field, then use this solution to build localized Wannier wavefunctions corresponding to Majorana quasiparticles. A low-energy tight-binding theory describing the intervortex hopping of these is then derived, and its topological properties---which depend crucially on the signs of the imaginary intervortex hopping parameters---are studied. We show that the energy gap between the Majorana bands may be either topologically trivial or nontrivial, depending on whether the Chern number contributions from the Majorana bands and those from the background superconducting condensate add constructively or destructively. This topology directly affects the temperature-dependent thermal Hall conductivity, which we also calculate. △ Less

Submitted 26 June, 2015; v1 submitted 22 June, 2015; originally announced June 2015.

Comments: 8 pages, 6 figures. Added references and added discussion regarding edge states

Journal ref: Phys. Rev. B 92, 134520 (2015)

Incommensurate spin density wave at a ferromagnetic quantum critical point in a three-dimensional parabolic semimetal

Authors: James M. Murray, Oskar Vafek, Leon Balents

Abstract: We explore the ferromagnetic quantum critical point in a three-dimensional semimetallic system with upward- and downward-dispersing bands touching at the Fermi level. Evaluating the static spin susceptibility to leading order in the coupling between the fermions and the fluctuating ferromagnetic order parameter, we find that the ferromagnetic quantum critical point is masked by an incommensurate,… ▽ More We explore the ferromagnetic quantum critical point in a three-dimensional semimetallic system with upward- and downward-dispersing bands touching at the Fermi level. Evaluating the static spin susceptibility to leading order in the coupling between the fermions and the fluctuating ferromagnetic order parameter, we find that the ferromagnetic quantum critical point is masked by an incommensurate, longitudinal spin density wave phase. We first analyze an idealized model which, despite having strong spin-orbit coupling, still possesses O(3) rotational symmetry generated by the total angular momentum operator. In this case, the direction of the incommensurate spin density wave propagation can point anywhere, while the magnetic moment is aligned along the direction of propagation. Including symmetry-allowed anisotropies in the fermion dispersion and the coupling to the order parameter field, however, the ordering wavevector instead breaks a discrete symmetry and aligns along either the [111] or [100] direction, depending on the signs and magnitudes of these two types of anisotropy. △ Less

Submitted 28 May, 2015; originally announced May 2015.

Comments: 12 pages

Journal ref: Phys. Rev. B 92, 035137 (2015)

Distinguishing spin-orbit coupling and nematic order in the electronic spectrum of iron-based superconductors

Authors: Rafael M. Fernandes, Oskar Vafek

Abstract: The low-energy electronic states of the iron-based superconductors are strongly affected by both spin-orbit coupling and, when present, by the nematic order. These two effects have different physical origins, yet they can lead to similar gap features in the electronic spectrum. Here we show how to disentangle them experimentally in the iron superconductors with one Fe plane per unit cell. Although… ▽ More The low-energy electronic states of the iron-based superconductors are strongly affected by both spin-orbit coupling and, when present, by the nematic order. These two effects have different physical origins, yet they can lead to similar gap features in the electronic spectrum. Here we show how to disentangle them experimentally in the iron superconductors with one Fe plane per unit cell. Although the splitting of the low energy doublet at the Brillouin zone center ($Γがんま$-point) can be due to either the spin-orbit coupling or the nematic order, or both, the degeneracy of each of the doublet states at the zone corner ($M$-point) is protected by the space group symmetry even when spin-orbit coupling is taken into account. Therefore, any splitting at $M$ must be due to lowering of the crystal symmetry, such as due to the nematic order. We further analyze a microscopic tight-binding model with two different contributions to the nematic order: $d_{xz}/d_{yz}$ onsite energy anisotropy and the $d_{xy}$ hopping anisotropy. We find that a precise determination of the former, which has been widely used to characterize the nematic phase, requires a simultaneous measurement of the splittings of the $Γがんま$-point doublet and at the two low-energy $M$-point doublets. We also discuss the impact of twin domains and show how our results shed new light on ARPES measurements in the normal state of these materials. △ Less

Submitted 16 November, 2014; originally announced November 2014.

Comments: 7 pages, 4 figures

Journal ref: Phys. Rev. B 90, 214514 (2014)

Emergence of superconductivity in a doped single-valley quadratic band crossing system of spin-1/2 fermions

Authors: Kelly Ann Pawlak, James M. Murray, Oskar Vafek

Abstract: For two-dimensional single-valley quadratic band crossing systems with weak repulsive electron-electron interactions, we show that upon introducing a chemical potential, particle-hole order is suppressed and superconductivity becomes the leading instability. In contrast to the two-valley case realized in bilayer graphene, the single-valley quadratic band touching is protected by crystal symmetries… ▽ More For two-dimensional single-valley quadratic band crossing systems with weak repulsive electron-electron interactions, we show that upon introducing a chemical potential, particle-hole order is suppressed and superconductivity becomes the leading instability. In contrast to the two-valley case realized in bilayer graphene, the single-valley quadratic band touching is protected by crystal symmetries, and the different symmetries and number of fermion flavors can lead to distinct phase instabilities. Our results are obtained using a weak-coupling Wilsonian renormalization group procedure on a low-energy effective Hamiltonian relevant for describing electrons on checkerboard or kagomé lattices. In 4-fold symmetric systems we find that $d$-wave and $s$-wave superconductivity are realized for short-ranged (Hubbard) and longer-ranged (forward scattering), respectively. In the 6-fold symmetric case, we find either $s$-wave superconductivity or no superconducting instability. △ Less

Submitted 13 November, 2014; v1 submitted 13 November, 2014; originally announced November 2014.

Comments: 7+3 Pages, 9 Figures

Journal ref: Phys. Rev. B 91, 134509 (2015)

Berry phases and the intrinsic thermal Hall effect in high temperature cuprate superconductors

Authors: Vladimir Cvetkovic, Oskar Vafek

Abstract: The Bogoliubov quasiparticles move in a practically uniform magnetic field in the vortex state of high temperature cuprate superconductors. Do the quasiparticles experience a Lorentz force when set in motion by an externally applied heat current ${\bf j}_Q$, bending their trajectories and causing the temperature gradient perpendicular to ${\bf j}_Q$ and the applied field ${\bf H}$, or is the therm… ▽ More The Bogoliubov quasiparticles move in a practically uniform magnetic field in the vortex state of high temperature cuprate superconductors. Do the quasiparticles experience a Lorentz force when set in motion by an externally applied heat current ${\bf j}_Q$, bending their trajectories and causing the temperature gradient perpendicular to ${\bf j}_Q$ and the applied field ${\bf H}$, or is the thermal Hall effect a consequence of Berry phases as in an intrinsic anomalous Hall effect of a semiconductor/metal with spin-orbit coupling? Here we show that it is the latter, and for the first time, calculate the temperature, ${\bf H}$-field and the $d$-wave pairing gap $Δでるた$ dependence of the intrinsic thermal Hall conductivity, $κかっぱ_{xy}$. We find that the intrinsic contribution to $κかっぱ_{xy}$ displays a rapid onset with increasing temperature, which compares favourably with existing experiments at high ${\bf H}$-fields on the highest purity samples. This finding may help to settle a much-debated question of the bulk value of the pairing strength in cuprate superconductors in magnetic field. △ Less

Submitted 25 July, 2014; v1 submitted 22 July, 2014; originally announced July 2014.

Comments: 5 pages, 3 figures; Supplementary: 11 pages, 3 figures; v2: Figure 1 updated

Journal ref: Nature Communications 6, 6518 (2015)

Renormalization group study of interaction-driven quantum anomalous Hall and quantum spin Hall phases in quadratic band crossing systems

Authors: James M. Murray, Oskar Vafek

Abstract: It is shown that topological insulating phases driven by interactions can be realized without the need for spin-orbit coupling or large intersite interaction in a two-dimensional system of spin-1/2 fermions with a single pair of parabolically touching bands. Using a weak-coupling, Wilsonian renormalization group procedure, we show that a quantum anomalous Hall phase is realized for Hubbard interac… ▽ More It is shown that topological insulating phases driven by interactions can be realized without the need for spin-orbit coupling or large intersite interaction in a two-dimensional system of spin-1/2 fermions with a single pair of parabolically touching bands. Using a weak-coupling, Wilsonian renormalization group procedure, we show that a quantum anomalous Hall phase is realized for Hubbard interaction, while a quantum spin Hall phase is favored for longer-ranged interactions. △ Less

Submitted 14 May, 2014; v1 submitted 19 March, 2014; originally announced March 2014.

Comments: 4+3 pages. Published version

Journal ref: Phys. Rev. B 89, 201110 (2014)

Thermodynamic signature of a magnetic-field-driven phase transition within the superconducting state of an underdoped high-temperature superconductor

Authors: J. B. Kemper, O. Vafek, J. B. Betts, F. F. Balakirev, W. N. Hardy, Ruixing Liang, D. A. Bonn, G. S. Boebinger

Abstract: More than a quarter century after the discovery of the high temperature superconductor (HTS) YBa$_2$Cu$_3$O$_{6+δでるた}$ (YBCO) studies continue to uncover complexity in its phase diagram. In addition to HTS and the pseudogap there is growing evidence for multiple phases with boundaries which are functions of temperature ($T$), doping (p) and magnetic field. Here we report the low temperature electroni… ▽ More More than a quarter century after the discovery of the high temperature superconductor (HTS) YBa$_2$Cu$_3$O$_{6+δでるた}$ (YBCO) studies continue to uncover complexity in its phase diagram. In addition to HTS and the pseudogap there is growing evidence for multiple phases with boundaries which are functions of temperature ($T$), doping (p) and magnetic field. Here we report the low temperature electronic specific heat (C$_{elec}$) of YBCO6.47 (p=0.08) up to a magnetic field (H) of 34.5 teslas (T), a poorly understood region of the underdoped H-$T$-p phase space. We observe two regimes: below a characteristic magnetic field H'$\approx$10 T, C$_{elec}/T$ obeys an expected H$^{1/2}$ behavior, however, near H' there is a sharp inflection followed by a linear-in-H behavior. H' rests deep within the superconducting phase and the linear-in-H behavior is observed in the zero resistance regime. In the limit of zero temperature, C$_{elec}/T$ is proportional to the zero-energy electronic density of states. Thus this inflection is evidence of a magnetic-field-driven quantum phase transition. △ Less

Submitted 9 September, 2015; v1 submitted 14 March, 2014; originally announced March 2014.

Comments: 10 pages, 5 figures

Excitonic and superconducting orders from repulsive interaction on the doped honeycomb bilayer

Authors: James M. Murray, Oskar Vafek

Abstract: Using a weak-coupling renormalization group formalism, we study competing ordered phases for repulsively interacting fermions on the bilayer honeycomb lattice away from half-filling, which is realized experimentally as doped bilayer graphene. As electrons are added to the system, excitonic order is suppressed, and unconventional superconductivity appears generically in its place. In general it is… ▽ More Using a weak-coupling renormalization group formalism, we study competing ordered phases for repulsively interacting fermions on the bilayer honeycomb lattice away from half-filling, which is realized experimentally as doped bilayer graphene. As electrons are added to the system, excitonic order is suppressed, and unconventional superconductivity appears generically in its place. In general it is found that the maximum critical temperature for superconductivity appears directly adjacent to the dome of particle-hole order, illustrating the importance of fluctuations in these channels for the formation of unconventional superconductivity. We obtain the phase diagram showing characteristic ordering temperatures for both short- and long-ranged interactions, and show that the most likely superconducting instabilities occur in $d$-wave, $f$-wave, and pair density wave channels. The nature of and competition between these phases are further analyzed using both free energy expansion and self-consistent mean-field theory. The effects of finite temperature and trigonal warping due to further-neighbor hopping are studied, and implications for experiments on bilayer graphene are discussed. △ Less

Submitted 23 May, 2014; v1 submitted 24 December, 2013; originally announced December 2013.

Comments: 22 pages. Added refs and minor changes. Published version

Journal ref: Phys. Rev. B 89, 205119 (2014)

Cyclotron resonance of single valley Dirac fermions in gapless HgTe quantum well

Authors: J. Ludwig, Yu. B. Vasilyev, N. N. Mikhailov, J. M. Poumirol, Z. Jiang, O. Vafek, D. Smirnov

Abstract: We report on Landau level spectroscopy studies of two HgTe quantum wells (QWs) near or at the critical well thickness, where the band gap vanishes. In magnetic fields up to $B$=16T, oriented perpendicular to the QW plane, we observe a $\sqrt{B}$ dependence for the energy of the dominant cyclotron resonance (CR) transition characteristic of two-dimensional Dirac fermions. The dominant CR line exhib… ▽ More We report on Landau level spectroscopy studies of two HgTe quantum wells (QWs) near or at the critical well thickness, where the band gap vanishes. In magnetic fields up to $B$=16T, oriented perpendicular to the QW plane, we observe a $\sqrt{B}$ dependence for the energy of the dominant cyclotron resonance (CR) transition characteristic of two-dimensional Dirac fermions. The dominant CR line exhibits either a single or double absorption lineshape for the gapless or gapped QW. Using an effective Dirac model, we deduce the band velocity of single valley Dirac fermions in gapless HgTe quantum wells, $v_F=6.4 \times10^5$ m/s, and interpret the double absorption of the gapped QW as resulting from the addition of a small relativistic mass. △ Less

Submitted 18 October, 2013; originally announced October 2013.

Superconductivity on the brink of spin-charge order in doped honeycomb bilayer

Authors: Oskar Vafek, James M. Murray, Vladimir Cvetkovic

Abstract: Using a controlled weak-coupling renormalization group approach, we establish the mechanism of unconventional superconductivity in the vicinity of spin or charge ordered excitonic states for the case of electrons on the Bernal stacked bilayer honeycomb lattice. With one electron per site this system exhibits nearly parabolically touching conduction and valence bands. Such a state is unstable towar… ▽ More Using a controlled weak-coupling renormalization group approach, we establish the mechanism of unconventional superconductivity in the vicinity of spin or charge ordered excitonic states for the case of electrons on the Bernal stacked bilayer honeycomb lattice. With one electron per site this system exhibits nearly parabolically touching conduction and valence bands. Such a state is unstable towards a spontaneous symmetry breaking, and repulsive interactions favor excitonic order, such as a charge nematic and/or a layer antiferromagnet. We find that upon adding charge carriers to the system, the excitonic order is suppressed, and unconventional superconductivity appears in its place, before it is replaced by a Fermi liquid. We focus on firmly establishing this phenomenon using the RG formalism within an idealized model with parabolic touching. △ Less

Submitted 19 September, 2013; v1 submitted 12 September, 2013; originally announced September 2013.

Comments: 5+3 pages, 4 figures; typos corrected

Journal ref: Phys. Rev. Lett. 112, 147002 (2014)

Unconventional superconductivity in a two-dimensional repulsive gas of fermions with spin-orbit coupling

Authors: Luyang Wang, Oskar Vafek

Abstract: We investigate the superconducting instability of a two-dimensional repulsive fermion gas with Rashba spin-orbit coupling $\al_R$. Using renormalization group approach, we find the superconducting transition temperature as a function of the dimensionless ratio $Θしーた={1}{2}m\al_R^2/E_F$ where $E_F=0$ when the smaller Fermi surface shrinks to a (Dirac) point. The general trend is that superconductivity… ▽ More We investigate the superconducting instability of a two-dimensional repulsive fermion gas with Rashba spin-orbit coupling $\al_R$. Using renormalization group approach, we find the superconducting transition temperature as a function of the dimensionless ratio $Θしーた={1}{2}m\al_R^2/E_F$ where $E_F=0$ when the smaller Fermi surface shrinks to a (Dirac) point. The general trend is that superconductivity is enhanced as $Θしーた$ increases, but in an intermediate regime $Θしーた\sim0.1$, a dome-like behavior appears. At a very small value of $Θしーた$, the angular momentum channel $j_z$ in which superconductivity occurs is quite high. With increasing $Θしーた$, $j_z$ decreases with a step of 2 down to $j_z=6$, after which we find the sequence $j_z=6, 4, 6, 2$, the last value of which continues to $Θしーた\rightarrow\infty$. In an extended range of $Θしーた$, the superconducting gap predominantly resides on the large Fermi surface, while Josephson coupling induces a much smaller gap on the small Fermi surface. Below the superconducting transition temperature, we apply mean field theory to derive the self-consistent equations and find the condensation energies. The state with the lowest condensation energy is an unconventional superconducting state which breaks time reversal symmetry, and in which singlet and triplet pairings are mixed. In general, these states are topologically nontrivial, and the Chern number of the state with total angular momentum $j_z$ is $C=2j_z$. △ Less

Submitted 30 July, 2013; originally announced July 2013.

Comments: 19 pages, 7 figures

Journal ref: Physica C Volume 497, Pages 6-18 (2014)

Dirac Fermions in Solids - from High Tc cuprates and Graphene to Topological Insulators and Weyl Semimetals

Authors: Oskar Vafek, Ashvin Vishwanath

Abstract: Understanding Dirac-like Fermions has become an imperative in modern condensed matter sciences: all across its research frontier, from graphene to high T$_c$ superconductors to the topological insulators and beyond, various electronic systems exhibit properties which can be well described by the Dirac equation. Such physics is no longer the exclusive domain of quantum field theories and other esot… ▽ More Understanding Dirac-like Fermions has become an imperative in modern condensed matter sciences: all across its research frontier, from graphene to high T$_c$ superconductors to the topological insulators and beyond, various electronic systems exhibit properties which can be well described by the Dirac equation. Such physics is no longer the exclusive domain of quantum field theories and other esoteric mathematical musings; instead, real physics of real systems is governed by such equations, and important materials science and practical implications hinge on our understanding of Dirac particles in two and three dimensions. While the physics that gives rise to the massless Dirac Fermions in each of the above mentioned materials is different, the low energy properties are governed by the same Dirac kinematics. The aim of this article is to review a selected cross-section of this vast field by highlighting the generalities, and contrasting the specifics, of several physical systems. △ Less

Submitted 10 June, 2013; originally announced June 2013.

Comments: 46 pages, 7 figures. Manuscript submitted to Annual Reviews of Condensed Matter Physics

Journal ref: Annual Review of Condensed Matter Physics Vol. 5: 83-112 (2014)

Space group symmetry, spin-orbit coupling and the low energy effective Hamiltonian for iron based superconductors

Authors: Vladimir Cvetkovic, Oskar Vafek

Abstract: We construct the symmetry adapted low energy effective Hamiltonian for the electronic states in the vicinity of the Fermi level in iron based superconductors. We use Luttinger's method of invariants, expanding about Gamma and M points in the Brillouin zone corresponding to two iron unit cell, and then matching the coefficients of the expansion to the 5- and 8-band models. We then use the method of… ▽ More We construct the symmetry adapted low energy effective Hamiltonian for the electronic states in the vicinity of the Fermi level in iron based superconductors. We use Luttinger's method of invariants, expanding about Gamma and M points in the Brillouin zone corresponding to two iron unit cell, and then matching the coefficients of the expansion to the 5- and 8-band models. We then use the method of invariants to study the effects of the spin-density wave order parameters on the electronic spectrum, with and without spin-orbit coupling included. Among the results of this analysis is the finding that the nodal spin-density wave is unstable once spin-orbit coupling is included. Similar analysis is performed for the A_{1g} spin singlet superconducting state. Without spin-orbit coupling there is one pairing invariant near the Gamma point, but two near the M point. This leads to an isotropic spectral gap at the hole Fermi surface near Gamma, but anisotropic near M. The relative values of these three parameters determine whether the superconducting state is s_{++}, s_{+-}, or nodal. Inclusion of spin-orbit coupling leads to additional mixing of spin triplet pairing, with one additional pairing parameter near Gamma and one near M. This leads to an anisotropic spectral gap near both hole and electron Fermi surfaces, the latter no longer cross, but rather split. △ Less

Submitted 12 April, 2013; originally announced April 2013.

Comments: 38 pages, 18 figures, 10 tables

Journal ref: Phys. Rev. B 88, 134510 (2013)

Quantum oscillations of the specific heat in d-wave superconductors with loop current order

Authors: Luyang Wang, Oskar Vafek

Abstract: We report numerical results of quantum oscillations of the specific heat in the vortex state of a $d_{x^2-y^2}$-wave superconductor in the presence of loop current order, which gives rise to Fermi pockets coexisting with nodal $d_{x^2-y^2}$-wave superconductivity. Within a lattice tight-binding model, we find that in an intermediate temperature range, the oscillations seem to approximately follow… ▽ More We report numerical results of quantum oscillations of the specific heat in the vortex state of a $d_{x^2-y^2}$-wave superconductor in the presence of loop current order, which gives rise to Fermi pockets coexisting with nodal $d_{x^2-y^2}$-wave superconductivity. Within a lattice tight-binding model, we find that in an intermediate temperature range, the oscillations seem to approximately follow Onsager relation with an effective charge comparable to the electric charge. However, the quasiparticle spectrum does not resemble Landau levels. In order to understand the origin of the oscillations, we also perform Franz-Tesanovic transformation in the presence of the loop order and find that in addition to scalar and Berry potentials, one component of the gauge invariant superfluid velocity couples to the low lying Dirac particles as a component of a vector potential. The magnetic field associated with this vector potential vanishes on average but is highly non-uniform in the magnetic unit cell. We also compare the results with the model without the loop order but with Zeeman-like coupling which also induces Fermi pockets in the superconducting state. △ Less

Submitted 9 July, 2013; v1 submitted 23 January, 2013; originally announced January 2013.

Comments: 13 pages, 15 figures

Journal ref: Physical Review B 88, 024506 (2013)

Topology and symmetry breaking in ABC trilayer graphene

Authors: Vladimir Cvetkovic, Oskar Vafek

Abstract: The effects of topology and electron-electron interactions on the phase diagram of ABC stacked trilayer graphene (TLG) at the neutrality point are investigated within a weak coupling renormalization group approach. We find that the leading instability of TLG with only a forward scattering density-density interaction is towards a mirror-breaking gapless state. Addition of a small, but finite back s… ▽ More The effects of topology and electron-electron interactions on the phase diagram of ABC stacked trilayer graphene (TLG) at the neutrality point are investigated within a weak coupling renormalization group approach. We find that the leading instability of TLG with only a forward scattering density-density interaction is towards a mirror-breaking gapless state. Addition of a small, but finite back scattering favors gapped phases, allowing us to make connections to the existing experiments on TLG. We identify a fundamental symmetry difference between TLG and bilayer graphene (BLG) which is responsible for disfavoring nematic states in TLG under the same conditions that favor nematic states in BLG. △ Less

Submitted 17 October, 2012; originally announced October 2012.

Comments: 5 pages, 3 figures; Supplementary Material included (5 pages)