Source author record

Oskar Vafek

Oskar Vafek 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

29works
6topics
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

29 published item(s)

preprint2022arXiv

Trions in Twisted Bilayer Graphene

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.

preprint2021arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2016arXiv

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

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.

preprint2015arXiv

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

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.

preprint2015arXiv

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

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.

preprint2015arXiv

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

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.

preprint2014arXiv

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

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.

preprint2014arXiv

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

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.

preprint2014arXiv

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

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.

preprint2014arXiv

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

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.

preprint2014arXiv

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

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.

preprint2013arXiv

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

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.

preprint2013arXiv

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

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.

preprint2013arXiv

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

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.

preprint2013arXiv

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

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.

preprint2013arXiv

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

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$.

preprint2012arXiv

Electronic multicriticality in bilayer graphene

We map out the possible ordered states in bilayer graphene at the neutrality point by extending the previous renormalization group treatment of many-body instabilities to finite temperature, trigonal warping and externally applied perpendicular electric field. We were able to analytically determine all outcomes of the RG flow equations for the 9 four-fermion coupling constants. While the full phase diagram exhibits a rich structure, we confirm that when forward scattering dominates, the only ordering tendency with divergent susceptibility at finite temperature is the nematic. At finite temperature this result is stable with respect to small back and layer imbalance scattering; further increasing their strength leads to the layer antiferromagnet. We also determine conditions for other ordered states to appear and compare our results to the special cases of attractive and repulsive Hubbard models where exact results are available.

preprint2012arXiv

Fermions on bilayer graphene: symmetry breaking for B=0 and nu=0

We extend previous analyses of fermions on a honeycomb bilayer lattice via weak-coupling renormalization group (RG) methods with extremely short-range and extremely long-range interactions to the case of finite-range interactions. In particular, we consider different types of interactions including screened Coulomb interactions, much like those produced by a point charge placed either above a single infinite conducting plate or exactly halfway between two parallel infinite conducting plates. Our considerations are motivated by the fact that, in some recent experiments on bilayer graphene there is a single gate while in others there are two gates, which can function as the conducting planes and which, we argue, can lead to distinct broken symmetry phases. We map out the phases that the system enters as a function of the range of the interaction. We discover that the system enters an antiferromagnetic phase for short ranges of the interaction and a nematic phase at long ranges, in agreement with previous work. While the antiferromagnetic phase results in a gap in the spectrum, the nematic phase is gapless, splitting the quadratic degeneracy points into two Dirac cones each. We also consider the effects of an applied magnetic field on the system in the antiferromagnetic phase via variational mean field theory. At low fields, we find that the antiferromagnetic order parameter, Delta(B)-Delta(0) \sim B^2. At higher fields, when omega_c > 2*Delta_0, we find that Delta(B)=omega_c/[ln(omega_c/Delta(0))+C], where C=0.67 and omega_c=eB/m^*c. We also determine the energy gap for creating electron-hole excitations in the system, and, at high fields, we find it to be a*omega_c+2*Delta(B), where a is a non-universal, interaction-dependent, constant.

preprint2012arXiv

Geometric phases of d-wave vortices in a model of lattice fermions

We study the local and topological features of Berry phases associated with the adiabatic transport of vortices in a d-wave superconductor of lattice fermions. At half filling, where the local Berry curvature must vanish due to symmetries, the phase associated with the exchange of two vortices is found to vanish as well, implying that vortices behave as bosons. Away from half filling, and in the limit where the magnetic length is large compared to the lattice constant, the local Berry curvature gives rise to an intricate flux pattern within the large magnetic unit cell. This renders the Berry phase associated with an exchange of two vortices highly path dependent. However, it is shown that "statistical" fluxes attached to the vortex positions are still absent. Despite the complicated profile of the Berry curvature away from half filling, we show that the average flux density associated with this curvature is tied to the average particle density. This is familiar from dual theories of bosonic systems, even though in the present case, the underlying particles are fermions.

preprint2012arXiv

Topology and symmetry breaking in ABC trilayer graphene

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.

preprint2011arXiv

Quantum Hall effect in a singly and doubly connected 3D topological insulator

The surface states of topological insulators, which behave as charged massless Dirac fermions, are studied in the presence of a quantizing uniform magnetic field. Using the method of D.H. Lee[1], analytical formula satisfied by the energy spectrum is found for a singly and doubly connected geometry. This is in turn used to argue that the way to measure the quantized Hall conductivity is to perform the Laughlin's flux ramping experiment and measure the charge transferred from the inner to the outer surface, analogous to the experiment in Ref.[2]. Unlike the Hall bar setup used currently, this has the advantage of being free of the contamination from the delocalized continuum of the surface edge states. In the presence of the Zeeman coupling, and/or interaction driven Quantum Hall ferromagnetism, which translate into the Dirac mass term, the quantized charge Hall conductivity sigma_{xy}=n e^2/h, with n=0,\pm 1,\pm 3,\pm 5... Backgating of one of the surfaces leads to additional Landau level splitting and in this case n can be any integer.

preprint2011arXiv

Spin-orbit coupling induced enhancement of superconductivity in a two-dimensional repulsive gas of fermions

We study a model of a two-dimensional repulsive Fermi gas with Rashba spin-orbit coupling $α_R$, and investigate the superconducting instability using renormalization group approach. We find that in general superconductivity is enhanced as the dimensionless ratio 1/2mα_R^2/E_F$ increases, resulting in unconventional superconducting states which break time reversal symmetry.

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

Interacting fermions on the honeycomb bilayer: from weak to strong coupling

Many-body instabilities of the half-filled honeycomb bilayer are studied using weak coupling renormalization group as well as strong coupling expansion. For spinless fermions and assuming parabolic degeneracy, there are 4-independent four-fermion contact couplings. While the dominant instability depends on the microscopic values of the couplings, the broken symmetry state is typically a gapped insulator with either broken inversion symmetry or broken time reversal symmetry, with a quantized anomalous Hall effect. Under certain conditions, the dominant instability may appear in the particle-particle (pairing) channel. For some non-generic fine-tuned initial conditions, weak coupling RG trajectories flow into the non-interacting fixed point, although generally we find runaway flows which we associate with ordering tendencies. Additionally, a tight binding model with nearest neighbor hopping and nearest neighbor repulsion is studied in weak and strong couplings and in each regime a gapped phase with inversion symmetry breaking is found. In the strong coupling limit, the ground state wavefunction is constructed for vanishing in-plane hopping but finite inter-plane hopping, which explicitly displays the broken inversion symmetry and a finite difference between the number of particles on the two layers. Finally, we discuss the spin-1/2 case and use Fierz identities to show that the number of independent 4-fermion contact couplings is 9. The corresponding RG equations in the spin-1/2 case are also presented, and used to show that, just as in strong coupling, the most dominant weak coupling instability of the repulsive Hubbard model (at half-filling) is an anti-ferromagnet.

preprint2010arXiv

Relaxation of Nuclear Magnetic Moments and Site-Selective NMR in d-Wave Superconductors

A new mechanism for relaxing the nuclear magnetic moments, in which a pair of spin-polarized BCS quasiparticles is emitted or absorbed, and which dominates at low temperature, is identified in type-II d-wave superconductors in an external magnetic field above Hc1. The results of the theory are compared with the NMR experiments on YBCO in high magnetic fields and found to agree without invoking antiferromagnetic order in the vortex core.

preprint2009arXiv

Many-body instability of Coulomb interacting bilayer graphene: RG approach

Low-energy electronic structure of (unbiased) bilayer graphene is made of two Fermi points with quadratic dispersions, if trigonal-warping and other high order contributions are ignored. We show that as a result of this qualitative difference from single-layer graphene, short-range (or screened Coulomb) interactions are marginally relevant. We use renormalization group to study their effects on low-energy properties of the system, and show that the two quadratic Fermi points spontaneously split into four Dirac points, at zero temperature. This results in a nematic state that spontaneously breaks the six-fold lattice rotation symmetry (combined with layer permutation) down to a two-fold one, with a finite transition temperature. Critical properties of the transition and effects of trigonal warping are also discussed.

preprint2007arXiv

Anomalous scaling and gapless fermions of d-wave superconductors in magnetic field

The problem of d-wave quasiparticles in a weakly disordered Abrikosov vortex lattice is studied. Starting with a periodic lattice, the topological structure of the magnetic crystal momenta of gapless fermions is found for the particle-hole symmetric case. If in addition the site centered inversion symmetry is present, both the location and the number of the gapless fermions can be determined using an index theorem. In the case of spatially aperiodic vortex array, Simon and Lee scaling is found to be violated due to a quantum anomaly. The electronic density of states is found to scale with the root-mean-square vortex displacement as $\sqrt{H}f(\bu_{rms}^2/ξ^2)$, while thermal conductivity is H-independent, but {\it different} from the H=0 case.