Source author record

Igor F. Herbut

Igor F. Herbut appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

47works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

47 published item(s)

preprint2022arXiv

Nematic Quantum Criticality in Dirac Systems

We investigate nematic quantum phase transitions in two different Dirac fermion models. The models feature twofold and fourfold, respectively, lattice rotational symmetries that are spontaneously broken in the ordered phase. Using negative-sign-free quantum Monte Carlo simulations and an $ε$-expansion renormalization group analysis, we show that both models exhibit continuous phase transitions. In contrast to generic Gross-Neveu dynamical mass generation, the quantum critical regime is characterized by large velocity anisotropies, with fixed-point values being approached very slowly. Hence both experimental and numerical investigations will not be representative of the infrared fixed point, but of a crossover regime characterized by drifting exponents.

preprint2022arXiv

p-wave superconductivity in Luttinger semimetals

We consider the three-dimensional spin-orbit-coupled Luttinger semimetal of "spin" $3/2$ particles in presence of weak attractive interaction in the $l=1$ (p-wave) channel, and determine the low-temperature phase diagram for both particle- and hole-dopings. The phase diagram depends crucially on the sign of the chemical potential, with two different states (with total angular momentum $j=0$ and $j=3$) competing on the hole-doped side, and three (one $j=1$ and two different $j=2$) states on the particle-doped side. The ground state condensates of Cooper pairs with the total angular momentum $j=1,2,3$ are selected by the quartic, and even sextic terms in the Ginzburg-Landau free energy. Interestingly, we find that all the p-wave ground states that appear in the phase diagram, while displaying different patterns of reduction of the rotational symmetry, preserve the time reversal symmetry. The resulting quasiparticle spectrum is either fully gapped or with point nodes, with nodal lines being absent.

preprint2020arXiv

d-wave superconductivity and Bogoliubov-Fermi surfaces in Rarita-Schwinger-Weyl semimetals

We uncover the properties of complex tensor (d-wave) superconducting order in three-dimensional Rarita-Schwinger-Weyl semimetals that host pseudospin-3/2 fermions at a fourfold linear band crossing point. Although the general theory of d-wave order was originally developed for materials displaying quadratic band touching, it directly applies to the case of semimetals with linear dispersion, several candidate compounds of which have been discovered experimentally very recently. The spin-3/2 nature of the fermions allows for the formation of spin-2 Cooper pairs which may be described by a complex second-rank tensor order parameter. In the case of linear dispersion, for the chemical potential at the Fermi point and at strong coupling, the energetically preferred superconducting state is the uniaxial nematic state, which preserves time-reversal symmetry and provides a full (anisotropic) gap for quasiparticle excitations. In contrast, at a finite chemical potential, we find that the usual weak-coupling instability is towards the "cyclic state", well known from the studies of multicomponent Bose-Einstein condensates, which breaks time reversal symmetry maximally, has vanishing average value of angular momentum, and features 16 small Bogoliubov--Fermi surfaces. The Rarita-Schwinger-Weyl semimetals provide therefore the first example of weakly coupled, three-dimensional, isotropic d-wave superconductors where the ground state is uniquely determined by the quartic expansion of the mean-field free energy, and not afflicted by the accidental degeneracy first noticed by Mermin over 40 years ago. We discuss the appearance and stability of the Bogoliubov-Fermi surfaces in absence of inversion symmetry in the electronic Hamiltonian, as in the case at hand.

preprint2020arXiv

Hydrodynamic transport in the Luttinger-Abrikosov-Beneslavskii non-Fermi liquid

We determine the shear viscosity and the dc electrical conductivity of interacting three-dimensional Luttinger semimetals, which have a quadratic band touching point in the energy spectrum, in the hydrodynamic regime. It is well-known that when the chemical potential is right at the band touching point the long-range Coulomb interaction induces the Luttinger-Abrikosov-Beneslavskii (LAB) phase at T = 0, which is an interacting, scale-invariant, non-Fermi-liquid state of electrons. Upon combining the renormalization-group (RG) analysis near the upper critical spatial dimension of four with the Boltzmann kinetic equation, we determine the universal ratio of viscosity over entropy, and the electrical dc conductivity of the system at the interacting LAB fixed point of the RG flow. The projection of the Coulomb interaction on the eigenstates of the system is found to play an important quantitative role for the scattering amplitude in the collision integral, and the so-called Auger-processes make a large numerical contribution to the inverse scattering time in the transport quantities. The obtained leading order result suggests that the universal ratio of the viscosity over entropy, when extrapolated to the physical three spatial dimensions, is above, but could be rather close to, the Kovtun-Son-Starinets lower bound.

preprint2016arXiv

Chiral symmetry breaking in three-dimensional quantum electrodynamics as fixed point annihilation

Spontaneous chiral symmetry breaking in three dimensional ($d=3$) quantum electrodynamics is understood as annihilation of an infrared-stable fixed point that describes the large-N conformal phase by another unstable fixed point at a critical number of fermions $N=N_c$. We discuss the root of universality of $N_c$ in this picture, together with some features of the phase boundary in the $(d,N)$ plane. In particular, it is shown that as $d\rightarrow 4$, $N_c\rightarrow 0$ with a constant slope, our best estimate of which suggests that $N_c = 2.89$ in $d=3$.

preprint2016arXiv

Emergent Lorentz symmetry near fermionic quantum critical points in two and three dimensions

We study the renormalization group flow of the velocities in the field theory describing the coupling of the massless quasi-relativistic fermions to the bosons through the Yukawa coupling, as well as with both bosons and fermions coupled to a fluctuating $U(1)$ gauge field in two and three spatial dimensions. Different versions of this theory describe quantum critical behavior of interacting Dirac fermions in various condensed-matter systems. We perform an analysis using one-loop $ε-$expansion about three spatial dimensions, which is the upper critical dimension in the problem. In two dimensions, we find that velocities of both charged fermions and bosons ultimately flow to the velocity of light, independently of the initial conditions, the number of fermionic and bosonic flavors, and the value of the couplings at the critical point. In three dimensions, due to the analyticity of the gauge field propagator, both the $U(1)$ charge and the velocity of light flow, which leads to a richer behavior than in two dimensions. We show that all three velocities ultimately flow to a common terminal velocity, which is non-universal and different from the original velocity of light. Therefore, emergence of the Lorentz symmetry in the ultimate infrared regime seems to be a rather universal feature of this class of theories in both two and three dimensions.

preprint2016arXiv

Excitonic instability of three-dimensional gapless semiconductors: Large-N theory

Three-dimensional gapless semiconductors with quadratic band touching, such as HgTe, $α$-Sn, or Pr$_2$Ir$_2$O$_7$ are believed to display a non-Fermi-liquid ground state due to long-range electron-electron interaction. We argue that this state is inherently unstable towards spontaneous formation of a (topological) excitonic insulator. The instability can be parameterized by a critical fermion number $N_c$. For $N < N_c$ the rotational symmetry is spontaneously broken, the system develops a gap in the spectrum, and features a finite nematic order parameter. To leading order in the 1/N expansion and in the static approximation, the analogy with the problem of dynamical mass generation in (2+1)-dimensional quantum electrodynamics yields $N_c = 16/[3π(π-2)]$. Taking the important dynamical screening effects into account, we find that $N_c \geq 2.6(2)$ and therefore safely above the physical value of $N = 1$. Some experimental consequences of the nematic ground state are discussed.

preprint2016arXiv

Gauge-field-assisted Kekulé quantum criticality

We study the quantum phase transition of $U(1)$ - charged Dirac fermions Yukawa-coupled to the Kekulé valence bond solid order parameter with $Z_3$ symmetry of the honeycomb lattice. The symmetry allows for the presence of the term in the action which is cubic in the Kekulé order parameter, and which is expected to prevent the quantum phase transition in question from being continuous. The Gross-Neveu-Yukawa theory for the transition is investigated using a perturbative renormalization group and within the $ε$ expansion close to four space-time dimensions. For a vanishing $U(1)$ charge we show that quantum fluctuations may render the phase transition continuous only sufficiently far away from 3+1 dimensions, where the validity of the conclusions based on the leading order $ε$ expansion appear questionable. In the presence of a fluctuating gauge field, on the other hand, we find quantum critical behavior even at weak coupling to appear close to 3+1 dimensions, that is, within the domain of validity of the perturbation theory. We also determine the renormalization group scaling of the cubic coupling at higher-loop orders and for a large number of Dirac fermions for vanishing charge.

preprint2016arXiv

Half vortex and fractional electrical charge in two dimensions

Despite fermion doubling, a two-dimensional quasi-relativistic spin-1/2 system can still lead to true fractionalization of electrical charge, when a massive ordered phase supports a "half-vortex". Such topological defect is possible when the order parameter in form of Dirac mass is described by two $U(1)$ angles, and each of them winds by an angle $π$ around a point. We demonstrate that such a mass configuration in an eight-dimensional Dirac Hamiltonian exhibits only a single bound zero mode, and therefore binds the charge of $e/2$. In graphene, for example, such an ordered phase is provided by the easy-plane spin-triplet Kekule valence bond solid. We argue that an application of an in-plane magnetic field can cause an excitonic instability toward such ordered phase, even for weak repulsion, when the on-site, nearest-neighbor and second neighbor components of it are of comparable strengths.

preprint2016arXiv

Note on "Quantum superconducting criticality in graphene and topological insulators"

We correct our previous conclusion regarding the fate of a charged quantum critical point across the superconducting transition for two dimensional massless Dirac fermion. Within the leading order $ε$ expansion, we now find that the requisite number of four-component Dirac fermion flavors ($N_f$) for the continuous phase transition through a charged critical point is $N_f>18.2699$. For $N_f\geq1/2$, the critical number of bosonic flavors for this transition is significantly reduced as compared to the value determined in the absence of the Dirac fermions in the theory.

preprint2016arXiv

Quadratic band touching with long range interactions in and out of equilibrium

Motivated by recent advances in cold atomic systems, we study the equilibrium and quench properties of two dimensional fermions with quadratic band touching at the Fermi level, in the presence of infinitely long range interactions. Unlike when only short range interactions are present, both nematic and quantum anomalous Hall (QAH) states state appear at weak interactions, separated by a narrow coexistence region, whose boundaries mark second and third order quantum phase transitions. After an interaction quench, the QAH order exhibits three distinct regions: persistent or damped oscillations and exponential decay to zero. In contrast, the nematic order always reaches a non-zero stationary value through power law damped oscillations, due to the interplay of the symmetry of the interaction and the specific topology of the quadratic band touching.

preprint2016arXiv

Superconducting quantum criticality in three-dimensional Luttinger semimetals

We study a simple model of three-dimensional fermions close to a quadratic band touching point, built from the celebrated Luttinger single-particle Hamiltonian and an attractive contact interaction between the particles. Such a system displays a novel quantum critical point between the semimetallic and an s-wave superconducting phase at which the low-energy "Luttinger fermions" are inextricably coupled to the order parameter fluctuations. The quantum critical point is perturbatively accessible near four spatial dimensions, where it features nontrivial scaling with dynamical exponent 1<z<2 and emergent rotational and particle-hole symmetries. Some features of the criticality, such as oscillatory corrections to scaling and its enhanced symmetry, are discussed.

preprint2015arXiv

Competition of density waves and quantum multicritical behavior in Dirac materials from functional renormalization

We study the competition of spin- and charge-density waves and their quantum multicritical behavior for the semimetal-insulator transitions of low-dimensional Dirac fermions. Employing the effective Gross-Neveu-Yukawa theory with two order parameters as a model for graphene and a growing number of other two-dimensional Dirac materials allows us to describe the physics near the multicritical point at which the semimetallic and the spin- and charge-density-wave phases meet. With the help of a functional renormalization group approach, we are able to reveal a complex structure of fixed points, the stability properties of which decisively depend on the number of Dirac fermions $N_f$. We give estimates for the critical exponents and observe crucial quantitative corrections as compared to the previous first-order $ε$ expansion. For small $N_f$, the universal behavior near the multicritical point is determined by the chiral Heisenberg universality class supplemented by a decoupled, purely bosonic, Ising sector. At large $N_f$, a novel fixed point with nontrivial couplings between all sectors becomes stable. At intermediate $N_f$, including the graphene case ($N_f = 2$) no stable and physically admissible fixed point exists. Graphene's phase diagram in the vicinity of the intersection between the semimetal, antiferromagnetic and staggered density phases should consequently be governed by a triple point exhibiting first-order transitions.

preprint2015arXiv

Critical $O(2)$ and $O(3)$ $ϕ^4$ theories near six dimensions

We consider $O(N)$-symmetric bosonic $ϕ^4$ field theories above four dimensions, and propose a new reformulation in terms of an irreducible tensorial field with a cubic and Yukawa terms. The $ϕ^4$ field theory so rewritten exhibits real and nontrivial IR-stable fixed points near and below six dimension, for low values of $N$ such as $N=2$ and $N=3$. The so-defined UV completions of the $O(2)$ and $O(3)$ models hence constitute precious examples of asymptotically safe quantum field theories. The possibility of an extension of our results to five dimensions is discussed.

preprint2015arXiv

Fermionic quantum criticality in honeycomb and $π$-flux Hubbard models: Finite-size scaling of renormalization-group-invariant observables from quantum Monte Carlo

We numerically investigate the critical behavior of the Hubbard model on the honeycomb and the $π$-flux lattice, which exhibits a direct transition from a Dirac semimetal to an antiferromagnetically ordered Mott insulator. We use projective auxiliary-field quantum Monte Carlo simulations and a careful finite-size scaling analysis that exploits approximately improved renormalization-group-invariant observables. This approach, which is successfully verified for the three-dimensional XY transition of the Kane-Mele-Hubbard model, allows us to extract estimates for the critical couplings and the critical exponents. The results confirm that the critical behavior for the semimetal to Mott insulator transition in the Hubbard model belongs to the Gross-Neveu-Heisenberg universality class on both lattices.

preprint2015arXiv

Mott multicriticality of Dirac electrons in graphene

We study the multicritical behavior for the semimetal-insulator transitions on graphene's honeycomb lattice using the Gross-Neveu-Yukawa effective theory with two order parameters: the SO(3) (Heisenberg) order parameter describes the antiferromagnetic transition, and the $\mathbb{Z}_2$ (Ising) order parameter describes the transition to a staggered density state. Their coupling induces multicritical behavior which determines the structure of the phase diagram close to the multicritical point. Depending on the number of fermion flavors $N_f$ and working in the perturbative regime in vicinity of three (spatial) dimensions, we observe first order or continuous phase transitions at the multicritical point. For the graphene case of $N_f=2$ and within our low order approximation, the phase diagram displays a tetracritical structure.

preprint2015arXiv

Nematic quantum criticality in three-dimensional Fermi system with quadratic band touching

We construct and discuss the field theory for tensorial nematic order parameter coupled to gapless four-component fermions at the quadratic band touching point in three (spatial) dimensions. Within a properly formulated epsilon-expansion this theory is found to have a quantum critical point, which describes the (presumably continuous) transition from the semimetal into a (nematic) Mott insulator. The latter phase breaks the rotational, but not the time-reversal symmetry, and may be relevant to materials such as gray tin or mercury telluride at low temperatures. The critical point represents a simple quantum analogue of the familiar classical isotropic-to-nematic transition in liquid crystals. The properties and the consequences of this quantum critical point are discussed. Its existence supports the scenario of the "fixed-point collision", according to which three-dimensional Fermi systems with quadratic band touching and long-range Coulomb interactions are unstable towards the gapped nematic ground state at low temperatures.

preprint2014arXiv

Interaction-induced anomalous quantum Hall state on the honeycomb lattice

We examine the existence of the interaction-generated quantum anomalous Hall phase on the honeycomb lattice. For the spinless model at half filling, the existence of a quantum anomalous Hall phase (Chern insulator phase) has been predicted using mean-field methods. However, recent exact diagonalization studies for small clusters with periodic boundary condition have not found a clear sign of an interaction-driven Chern insulator phase. We use exact diagonalization method to study properties of small clusters with open boundary condition and, contrary to previous studies, we find clear signatures of the topological phase transition for finite size clusters. We also examine applicability of the entangled-plaquette state (correlator-product state) ansatz to describe the ground states of the system. Within this approach the lattice is covered with plaquettes and the ground state wave-function is written in terms of the plaquette coefficients. Configurational weights can then be optimized using a variational Monte Carlo algorithm. Using the entangled-plaquette state ansatz we study the ground state properties of the system for larger system sizes and show that the results agree with the exact diagonalization results for small clusters. This confirms validity of the entangled-plaquette state ansatz to describe the ground states of the system and provides further confirmation of the existence of the quantum anomalous Hall phase in the thermodynamic limit, as predicted by the mean-field theory calculations.

preprint2014arXiv

Nematic, topological and Berry phases when a flat and a parabolic band touch

A (single flavor) quadratic band crossing in two dimensions is known to have a generic instability towards a quantum anomalous Hall (QAH) ground state for infinitesimal repulsive interactions. Here we introduce a generalization of a quadratic band crossing which is protected only by rotational symmetry. By focusing on the representative case of a parabolic and flat band touching, which also allows for a straightforward lattice realization, the interaction induced nematic phase is found generally to compete successfully with the QAH insulator, and to become the dominant instability in certain parts of the phase diagram already at weak coupling. The full phase diagram of the model, together with its topological properties, is mapped out using a perturbative renormalization group, strong coupling analysis, the mean-field theory. Interestingly, the Berry flux varies continuously in the single flavour limit with various control parameters.

preprint2014arXiv

Topological Mott insulator in three-dimensional systems with quadratic band touching

We argue that a three dimensional electronic system with the Fermi level at the quadratic band touching point such as HgTe could be unstable with respect to the spontaneous formation of the (topological) Mott insulator at arbitrary weak long range Coulomb interaction. The mechanism of the instability can be understood as the collision of Abrikosov's non-Fermi liquid fixed point with another, quantum critical, fixed point, which approaches it in the coupling space as the system's dimensionality $d\rightarrow d_{low} +$, with the "lower critical dimension" $2< d_{low}<4$. Arguments for the existence of the quantum critical point based on considerations in the large-$N$ limit in $d=3$, as well as close to $d=2$, are given. In the one-loop calculation we find that $d_{low} = 3.26$, and thus above, but not far from three dimensions. This translates into a temperature/energy window $(T_c, T_*)$ over which the NFL scaling should still be observable, before the Mott transition finally takes place at the critical temperature $T_c \sim T_* \exp(-z C/(d_{low} -d)^{1/2})$. We estimate $C=π/1.1$, dynamical critical exponent $z\approx 1.8$, and the temperature scale $k_B T_* \approx (4 m/ m_{el} ε^2)13.6 eV$, with $m$ as the band mass and $ε$ as the dielectric constant.

preprint2014arXiv

Zero-modes and global antiferromagnetism in strained graphene

A novel magnetic ground state is reported for the Hubbard Hamiltonian in strained graphene. When the chemical potential lies close to the Dirac point, the ground state exhibits locally both the Néel and ferromagnetic orders, even for weak Hubbard interaction. Whereas the Néel order parameter remains of the same sign in the entire system, the magnetization at the boundary takes the opposite sign from the bulk. The total magnetization this way vanishes, and the magnetic ground state is globally only an antiferromagnet. This peculiar ordering stems from the nature of the strain-induced single particle zero-energy states, which have support on one sublattice of the honeycomb lattice in the bulk, and on the other sublattice near the boundary of a finite system. We support our claim with the self-consistent numerical calculation of the order parameters, as well as by the Monte Carlo simulations of the Hubbard model in both uniformly and non-uniformly strained honeycomb lattice. The present result is contrasted with the magnetic ground state of the same Hubbard model in the presence of a true magnetic field (and for vanishing Zeeman coupling), which is exclusively Néel ordered, with zero local magnetization everywhere in the system.

preprint2013arXiv

Coupling, merging, and splitting Dirac points by electron-electron interaction

The manipulation and movement of Dirac points in the Brillouin zone by the electron-electron interaction is considered within leading order perturbation theory. At the merging point, an infinitesimal interaction is shown to cause opening of the gap or splitting of the Dirac points, depending on the inter- or intrasublattice nature of the merging and the sign of the interaction. The topology of the spectrum can therefore be efficiently changed by simply tuning the interaction between particles, as opposed to the usual careful band structure engineering. This is illustrated around the merging transition of one, two, and three dimensional Dirac-Weyl fermions. A simple Weyl-like Hamiltonian that describes the quadratic band-crossing in three dimensions is also proposed, and its stability under interactions is addressed.

preprint2013arXiv

Density of States Scaling at the Semimetal to Metal Transition in Three Dimensional Topological Insulators

The quantum phase transition between the three dimensional Dirac semimetal and the diffusive metal can be induced by increasing disorder. Taking the system of disordered $\mathbb{Z}_2$ topological insulator as an important example, we compute the single particle density of states by the kernel polynomial method. We focus on three regions: the Dirac semimetal at the phase boundary between two topologically distinct phases, the tricritical point of the two topological insulator phases and the diffusive metal, and the diffusive metal lying at strong disorder. The density of states obeys a novel single parameter scaling, collapsing onto two branches of a universal scaling function, which correspond to the Dirac semimetal and the diffusive metal. The diverging length scale critical exponent $ν$ and the dynamical critical exponent $z$ are estimated, and found to differ significantly from those for the conventional Anderson transition. Critical behavior of experimentally observable quantities near and at the tricritical point is also discussed.

preprint2013arXiv

Pinning the order: the nature of quantum criticality in the Hubbard model on honeycomb lattice

In numerical simulations, spontaneously broken symmetry is often detected by computing two-point correlation functions of the appropriate local order parameter. This approach, however, computes the square of the local order parameter, and so when it is {\it small}, very large system sizes at high precisions are required to obtain reliable results. Alternatively, one can pin the order by introducing a local symmetry breaking field, and then measure the induced local order parameter infinitely far from the pinning center. The method is tested here at length for the Hubbard model on honeycomb lattice, within the realm of the projective auxiliary field quantum Monte Carlo algorithm. With our enhanced resolution we find a direct and continuous quantum phase transition between the semi-metallic and the insulating antiferromagnetic states with increase of the interaction. The single particle gap in units of the Hubbard $U$ tracks the staggered magnetization. An excellent data collapse is obtained by finite size scaling, with the values of the critical exponents in accord with the Gross-Neveu universality class of the transition.

preprint2013arXiv

Quantum superconducting criticality in graphene and topological insulators

The field theory of the semimetal-superconductor quantum phase transition for graphene and surface states of topological insulators is presented. The Lagrangian possesses the global U(1) symmetry, with the self-interacting complex bosonic order-parameter and the massless Dirac fermions coupled through a Yukawa term. The same theory also governs the quantum critical behavior of graphene near the transition towards the bond-density-wave (Kekule) insulator. The local U(1) gauged version of the theory which describes the quantum semimetal-superconductor transition in the ultimate critical regime is also considered. Due to the Yukawa coupling the transitions are found to be always continuous, both with and without the fluctuating gauge field. The critical behavior is addressed within the dimensional regularization near four space-time dimensions, and the calculation of various universal quantities, including critical exponents and the universal mass-ratio, is reported.

preprint2013arXiv

Topological insulators in strained graphene at weak interaction

The nature of the electronic ground states in strained undoped graphene at weak interaction between electrons is discussed. After providing a lattice realization of the strain-induced axial magnetic field we numerically find the self-consistent solution for the time reversal symmetry breaking quantum anomalous Hall order-parameter, at weak second-nearest-neighbor repulsion between spinless fermions. The anomalous Hall state is obtained in both uniform and nonuniform axial magnetic fields, with the spatial profile of the order-parameter resembling that of the axial field itself. When the electron spin is included, the time reversal symmetric anomalous spin Hall state becomes slightly preferred energetically at half filling, but the additional anomalous Hall component should develop at a finite doping.

preprint2013arXiv

Universal conductivity of graphene in the ultra-relativistic regime

We calculate the optical (cutoff >> frequency >> temperature) conductivity in clean graphene in the ultimate low-energy regime, when retardation effects of the electromagnetic interaction become important and when the full Lorentz symmetry emerges. In contrast to what happens with the short range or with the Coulomb long-range instantaneous interactions, the optical conductivity is now no longer equal to its non interacting value, but acquires universal corrections in powers of the fine structure constant. The coefficient of the first order correction is computed, and found to be of order one. We also present the result for the conductivity in the large-N limit, with $N$ as the number of Dirac fermions species, to the order $1/N^2$.

preprint2012arXiv

Conserved charges of order-parameter textures in Dirac systems

A simple expression for the induced fermion current in the presence of a texture in mass-order-parameters in two-dimensional condensed-matter Dirac systems is derived using the representation theory of Clifford algebras. In particular, it is shown that every texture in three mutually anticommuting order parameters, in graphene for example, implies an induced density of a properly defined conserved charge. The sufficient condition for the general charge to be the familiar electrical charge is that the remaining two anticommuting order parameters allowed by the particle-hole symmetry are the two phase components of some superconducting order. This allows eight different types of electrically charged textures in graphene or in the $π$-flux Hamiltonian on the square lattice. Generalized charge of mass-textures on the surfaces of thin films of topological insulators, or in spinless Dirac fermions hopping on the honeycomb lattice is also discussed.

preprint2012arXiv

Isospin of topological defects in Dirac systems

We study the Dirac quasiparticles in $d$-dimensional lattice systems of electrons in the presence of domain walls ($d=1$), vortices ($d=2$), or hedgehogs ($d=3$) of superconducting and/or insulating, order parameters, which appear as mass terms in the Dirac equation. Such topological defects have been known to carry non-trivial quantum numbers such as charge and spin. Here we discuss their additional internal degree of freedom: irrespectively of the dimensionality of space and the nature of orders that support the defect, an extra mass-order-parameter is found to emerge in their core. Six linearly independent local orders, which close two mutually commuting three-dimensional Clifford algebras are proven to be in general possible. We show how the particle-hole symmetry restricts the defects to always carry the quantum numbers of a single effective isospin-1/2, quite independently of the values of their electric charge or true spin. Examples of this new degree of freedom in graphene and on surfaces of topological insulators are discussed.

preprint2012arXiv

Zero modes and charged Skyrmions in graphene bilayer

We show that the electric charge of the Skyrmion in the vector order parameters that characterize the quantum anomalous spin Hall state and the layer-antiferromagnet in a graphene bilayer is four and zero, respectively. The result is based on the demonstration that a vortex configuration in two broken symmetry states in bilayer graphene with the quadratic band crossing has the number of zero modes doubled relative to the single layer. The doubling can be understood as a result of Kramers' theorem implied by the "pseudo time reversal" symmetry of the vortex Hamiltonian. Disordering the quantum anomalous spin Hall state by Skyrmion condensation should produce a superconductor of an elementary charge 4e.

preprint2011arXiv

Inhomogeneous magnetic catalysis on graphene's honeycomb lattice

We investigate the ordering instability of interacting (and for simplicity, spinless) fermions on graphene's honeycomb lattice by numerically computing the Hartree self-consistent solution for the charge-density-wave order parameter in presence of both uniform and non-uniform magnetic fields. For a uniform field the overall behavior of the order parameter is found to be in accord with the continuum theory. In the inhomogeneous case, the spatial profile of the order parameter resembles qualitatively the form of the magnetic field itself, at least when the interaction is not overly strong. We find that right at the zero-field critical point of the infinite system the local order parameter scales as the square-root of the local strength of the magnetic field, apparently independently of the assumed field's profile. The finite size effects on various parameters of interest, such as the critical interaction and the universal amplitude ratio of the interaction-induced gap to the Landau level energy at criticality are also addressed.

preprint2011arXiv

Spectrum of the Dirac Hamiltonian with the mass-hedgehog in arbitrary dimension

It is shown that the square of the Dirac Hamiltonian with the isotropic mass-hedgehog potential in d dimensions is the number operator of fictitious bosons and fermions over d quantum states. This result allows one to obtain the complete spectrum and degeneracies of the Dirac Hamiltonian with the hedgehog mass configuration in any dimension. The result pertains to low-energy states in the core of a general superconducting or insulating vortex in graphene in two dimensions, and in the superconducting vortex at the topological - trivial insulator interface in three dimensions, for example. The spectrum in d=2 is also understood in terms of the underlying accidental SU(2) symmetry and the supersymmetry of the Hamiltonian.

preprint2011arXiv

Supersymmetric Runge-Lenz-Pauli vector for Dirac vortex in topological insulators and graphene

The Dirac mass-vortex at the surface of a topological insulator or in graphene is considered. Within the linear approximation for the vortex amplitude's radial dependence, the spectrum is a series of degenerate bound states, which can be classified by a set of accidental SU(2) and supersymmetry generators (I. F. Herbut and C.-K. Lu, Phys. Rev. B 83 125412 (2011)). Here we discuss further the properties and manifestations of the supersymmetry of the vortex Hamiltonian, and point out some interesting analogies to the Runge-Lenz-Pauli vector in the non-relativistic hydrogen atom. Symmetry breaking effects due to a finite chemical potential, and the Zeeman field are also analyzed. We find that a residual accidental degeneracy remains only in the special case of equal magnitudes of both terms, whereas otherwise it becomes removed entirely.

preprint2011arXiv

Time reversal, fermion doubling, and the masses of lattice Dirac fermions in three dimensions

Motivated by recent examples of three-dimensional lattice Hamiltonians with massless Dirac fermions in their (bulk) spectrum, I revisit the problem of fermion doubling on bipartite lattices. The number of components of the Dirac fermion in a time-reversal and parity invariant d-dimensional lattice system is determined by the minimal representation of the Clifford algebra of $d+1$ Hermitian Dirac matrices that allows a construction of the time-reversal operator with the square of unity, and it equals $2^d$ for $d=2,3$. Possible mass-terms for (spinless) Dirac fermions are listed and discussed. In three dimensions there are altogether eight independent masses, out of which four are even, and four are odd under time reversal. A specific violation of time-reversal symmetry that leads to (minimal) four-component massless Dirac fermion in three dimensions at low energies is constructed.

preprint2010arXiv

Antilinear spectral symmetry and the vortex zero-modes in topological insulators and graphene

We construct the general extension of the four-dimensional Jackiw-Rossi-Dirac Hamiltonian that preserves the antilinear reflection symmetry between the positive and negative energy eigenstates. Among other systems, the resulting Hamiltonian describes the s-wave superconducting vortex at the surface of the topological insulator, at a finite chemical potential, and in the presence of both Zeeman and orbital couplings to the external magnetic field. Here we find that the bound zero-mode exists only when the Zeeman term is below a critical value. Other physical realizations pertaining to graphene are considered, and some novel zero-energy wave functions are analytically computed.

preprint2010arXiv

Conductivity of interacting massless Dirac particles in graphene: Collisionless regime

We provide detailed calculation of the a.c. conductivity in the case of 1/r-Coulomb interacting massless Dirac particles in graphene in the collisionless limit when ω>> T. The analysis of the electron self-energy, current vertex function and polarization function, which enter into the calculation of physical quantities including the a.c. conductivity, is carried out by checking the Ward-Takahashi identities associated with the electrical charge conservation and making sure that they are satisfied at each step. We adopt a variant of the dimensional regularization of Veltman and t'Hooft by taking the spatial dimension D=2-ε, for ε> 0. The procedure adopted here yields a result for the conductivity correction which, while explicitly preserving charge conservation laws, is nevertheless different from the results reported previously in literature.

preprint2010arXiv

Pairing symmetry and vortex zero-mode for superconducting Dirac fermions

We study the vortex zero-energy bound states in presence of pairing among the low-energy Dirac fermions on the surface of a topological insulator. The pairing symmetries considered include the $s$-wave, $p$-wave, and, in particular, the mixed-parity symmetry, which arises in absence of the inversion symmetry on the surface. The zero-mode is analyzed within the generalized Jackiw-Rossi-Dirac Hamiltonian that contains a momentum-dependent mass-term, and includes the effects of the electromagnetic gauge field and the Zeeman coupling as well. At a finite chemical potential, as long as the spectrum without the vortex is fully gapped, the presence of a single Fermi surface with a definite helicity always leads to one Majorana zero-mode, in which both electron's spin projections participate. In particular, the critical effects of the Zeeman coupling on the zero-mode are discussed.

preprint2010arXiv

Topological insulator in the core of the superconducting vortex in graphene

The core of the vortex in a general superconducting order parameter in graphene is argued to be ordered, with the possible order parameters forming the algebra U(1) X Cl(3), where Cl(3) is the three dimensional Clifford algebra. A sufficiently strong Zeeman coupling of the magnetic field of the vortex to the electron spin breaks the degeneracy in the core in favor of the anomalous quantum Hall state. I consider a variety of superconducting condensates on the honeycomb lattice and demonstrate the surprising universality of this result. A way to experimentally determine the outcome of the possible competition between different types of orders in the core is proposed.

preprint2010arXiv

Unconventional superconductivity on honeycomb lattice: the theory of Kekule order parameter

A spatially non-uniform superconducting phase is proposed as the electronic variational ground state for the attractive interactions between nearest neighbors on graphene's honeycomb lattice, close to and right at the filling one half. The state spontaneously breaks the translational invariance of the lattice into the Kekule pattern of bond order parameters, and it is gapped, spin triplet, and odd under the sublattice exchange. With the increase of attractive interactions we first find the transition from the semimetallic phase into the p-Kekule superconductor, defined as being odd under the exchange of Dirac points, with the additional discontinuous superconductor-superconductor transition into the even s-Kekule state, deep within the superconducting phase. Topological excitations of the Kekule superconductor and its competition with other superconducting states on the honeycomb lattice are discussed.

preprint2010arXiv

Zero-energy states of massive Dirac equation in magnetic fields

The Dirac equation with a U(1) vortex in the mass-term is solved in the presence of magnetic-like fields at zero energy. By drawing an analogy to classical mechanics, it is shown that the four-component Dirac equation in arbitrary magnetic field always yields one zero-energy state. In the time-reversal preserving, pseudo-magnetic field, however, the number of zero-energy states may depend on the field's profile and sign. Some explicit examples are worked out. Possible implications of these results for the charge of the vortex and for the behavior of graphene in magnetic field are discussed.

preprint2009arXiv

Coulomb interaction at the metal-insulator critical point in graphene

We compute the renormalization group flow of the long-ranged electron-electron interaction at the Gross-Neveu quantum critical point between the semimetal and the excitonic insulator in graphene, perturbatively in the small parameter $ε=d-1$, with $d$ as the spatial dimension. The $O(ε)$ correction to the usual beta-function makes the long-range interaction only more irrelevant at the critical than at the Gaussian fixed point. A weak long-range tail of the Coulomb interaction is found to be marginally irrelevant also in arbitrary dimension when the number of Dirac fermions is large. Its ultimate irrelevancy notwithstanding, it is shown that the metal-insulator transition may still be induced by increasing only the long-range tail of the Coulomb interaction.

preprint2007arXiv

Excitons in QED3 and spin response in a phase-fluctuating d-wave superconductor

We study the particle-hole exciton mode in the QED3 theory of a phase-fluctuating d-wave superconductor in ladder approximation. We derive a Schroedinger-like equation for the exciton bound state and determine the conditions for its existence. We find the dispersion of this mode below the particle-hole continuum and compare our results with the resonance observed in neutron scattering measurements in cuprates.

preprint2004arXiv

Quantum fluctuations, pseudogap, and the T=0 superfluid density in strongly correlated d-wave superconductors

I study the effect of Coulomb interaction on superconducting order in a d-wave lattice superconductor at T=0 by considering the superconducting saddle point in the two dimensional t-J-U model with an on-site repulsion U. The theory of low energy phase fluctuations around this saddle point is derived in terms of the effective hard-core bosons (representing the density of spin-up electrons and the order parameter phase), interacting with the fluctuating density of spin-down electrons. Whereas the saddle point value of the gap is found to monotonically increase towards half filling, the phase stiffness at T=0 has a maximum, and then decreases with further underdoping. Right at half filling the stiffness vanishes for large U. This argues that the pseudogap phenomenon of the type observed in cuprates is in principle possible without a development of any competing order, purely as a result of growing correlations in the superconducting state. The effects of finite temperature and disorder are discussed qualitatively.

preprint2003arXiv

Absence of U(1) spin liquids in two dimensions

Many popular models of fractionalized spin liquids contain neutral fermionic spinon excitations on a Fermi surface, carrying unit charges under a compact U(1) gauge force. We argue that instanton effects generically render such states unstable to confinement in two spatial dimensions, so that all elementary excitations are gauge neutral, and there is no spinon Fermi surface. Similar results are expected to apply to SU(2) spin liquids. However, fractionalized states can appear when the gauge symmetry is broken down to a discrete subgroup by the Higgs mechanism. Our argument generalizes earlier results on confinement in the pure gauge theory, and on the instability of the U(1) staggered flux and algebraic spin liquids with a Dirac spectrum for fermionic spinons.

preprint2003arXiv

Permanent confinement in the compact QED3 with fermionic matter

We argue that the compact three dimensional electrodynamics with massless fermions is always confining, in spite of the bare interaction between the magnetic monopoles being rendered logarithmic by fermions. The effect is caused by screening by other dipoles, which transforms the logarithmic back into the Coulomb interaction at large distances. Possible implications for the chiral symmetry breaking for fermions are discussed, and the global phase diagram of the theory is proposed.

preprint2003arXiv

Self-consistent theory of compact QED3 with relativistic fermions

We consider three dimensional quantum electrodynamics (cQED3) with massless relativistic fermions coupled to a compact gauge field using a combined perturbative variational approach. Coupling to matter renders the bare interaction between magnetic monopoles logarithmic at large distances, suggesting the possibility of a confinement-deconfinement transition of the Berezinskii-Kosterlitz-Thouless type in the theory. Our self-consistent calculation suggests, however, that screening effects always destabilise the confined phase, in agreement with the previous renormalisation group study of the same model.

preprint2002arXiv

Fine structure of chiral symmetry breaking in the $QED_3$ theory of underdoped high-$T_c$ superconductors

We study the effects of the leading irrelevant perturbations on chiral symmetry breaking in the effective $QED_3$ theory of d-wave superconductor-insulator transition in underdoped cuprates. For weak symmetry breaking, the effect of a perturbation on energies of various insulating states can be classified according to its engineering dimension in the maximally symmetric theory. Considering the velocity anisotropy, repulsive interactions, and higher order derivatives we show that the insulating state with the lowest energy is the spin density wave.