Continuum Theory of Edge States of Topological Insulators: Variational Principle and Boundary Conditions

We develop a continuum theory to model low energy excitations of a generic four-band time reversal invariant electronic system with boundaries. We propose a variational energy functional for the wavefunctions which allows us derive natural boundary conditions valid for such systems. Our formulation is particularly suited to develop a continuum theory of the protected edge/surface excitations of topological insulators both in two and three dimensions. By a detailed comparison of our analytical formulation with tight binding calculations of ribbons of topological insulators modeled by the Bernevig-Hughes-Zhang (BHZ) hamiltonian, we show that the continuum theory with the natural boundary condition provides an appropriate description of the low energy physics. As a spin-off, we find that in a certain parameter regime, the gap that arises in topological insulator ribbons of finite width due to the hybridization of edges states from opposite edges, depends non-monotonically on the ribbon width and can nearly vanish at certain "magic widths".Comment: 8 pages, 5 figure

Fermionic Superfluid from a Bilayer Band Insulator in an Optical Lattice

We propose a model to realize a fermionic superfluid state in an optical lattice circumventing the cooling problem. Our proposal exploits the idea of tuning the interaction in a characteristically low entropy state, a band-insulator in an optical bilayer system, to obtain a superfluid. By performing a detailed analysis of the model including fluctuations and augmented by a variational quantum Monte Carlo calculations of the ground state, we show that the superfluid state obtained has high transition temperature of the order of the hopping energy. Our system is designed to suppress other competing orders such as a charge density wave. We suggest a laboratory realization of this model via an orthogonally shaken optical lattice bilayer.Comment: 5 pages, 7 figures, typos fixed, figures modifie

Synchronous and Asynchronous Mott Transitions in Topological Insulator Ribbons

We address how the nature of linearly dispersing edge states of two dimensional (2D) topological insulators evolves with increasing electron-electron correlation engendered by a Hubbard like on-site repulsion $U$ in finite ribbons of two models of topological band insulators. Using an inhomogeneous cluster slave rotor mean-field method developed here, we show that electronic correlations drive the topologically nontrivial phase into a Mott insulating phase via two different routes. In a synchronous transition, the entire ribbon attains a Mott insulating state at one critical $U$ that depends weakly on the width of the ribbon. In the second, asynchronous route, Mott localization first occurs on the edge layers at a smaller critical value of electronic interaction which then propagates into the bulk as $U$ is further increased until all layers of the ribbon become Mott localized. We show that the kind of Mott transition that takes place is determined by certain properties of the linearly dispersing edge states which characterize the topological resilience to Mott localization.Comment: 4+ pages, 5 figure

Sensory organ like response determines the magnetism of zigzag-edged honeycomb nanoribbons

We present an analytical theory for the magnetic phase diagram for zigzag edge terminated honeycomb nanoribbons described by a Hubbard model with an interaction parameter U . We show that the edge magnetic moment varies as ln U and uncover its dependence on the width W of the ribbon. The physics of this owes its origin to the sensory organ like response of the nanoribbons, demonstrating that considerations beyond the usual Stoner-Landau theory are necessary to understand the magnetism of these systems. A first order magnetic transition from an anti-parallel orientation of the moments on opposite edges to a parallel orientation occurs upon doping with holes or electrons. The critical doping for this transition is shown to depend inversely on the width of the ribbon. Using variational Monte-Carlo calculations, we show that magnetism is robust to fluctuations. Additionally, we show that the magnetic phase diagram is generic to zigzag edge terminated nanostructures such as nanodots. Furthermore, we perform first principles modeling to show how such magnetic transitions can be realized in substituted graphene nanoribbons.Comment: 5 pages, 5 figure

Coexistence of magnetism and superconductivity in a t-J bilayer

We investigate coexistence of antiferromagnetic and superconducting correlations in bilayered materials using a two-dimensional t-J model with couplings across the layers using variational Monte Carlo calculations. It is found that the underdoped regime supports a coexisting phase, beyond which the (d-wave) superconducting state becomes stable. Further, the effects of interplanar coupling parameters on the magnetic and superconducting correlations as a function of hole doping are studied in details. The magnetic correlations are found to diminish with increasing interplanar hopping away from half filling, while the exchange across the layers strengthens interplanar antiferromagnetic correlations both at and away from half filling. The superconducting correlations show more interesting features where larger interplanar hopping considerably reduces planar correlations at optimal doping, while an opposite behaviour, i.e. stabilisation of the superconducting state is realised in the overdoped regime, with the interplanar exchange all the while playing a dormant role.Comment: 8 pages, 9 figures, RevTex4, Submitted to Phys. Rev.

Helical edge states in multiple topological mass domains

The two-dimensional topological insulating phase has been experimentally discovered in HgTe quantum wells (QWs). The low-energy physics of two-dimensional topological insulators (TIs) is described by the Bernevig-Hughes-Zhang (BHZ) model, where the realization of a topological or a normal insulating phase depends on the Dirac mass being negative or positive, respectively. We solve the BHZ model for a mass domain configuration, analyzing the effects on the edge modes of a finite Dirac mass in the normal insulating region (soft-wall boundary condition). We show that at a boundary between a TI and a normal insulator (NI), the Dirac point of the edge states appearing at the interface strongly depends on the ratio between the Dirac masses in the two regions. We also consider the case of multiple boundaries such as NI/TI/NI, TI/NI/TI and NI/TI/NI/TI.Comment: 11 pages, 15 figure

Helical edge states in multiple topological mass domains

The two-dimensional topological insulating phase has been experimentally discovered in HgTe quantum wells (QWs). The low-energy physics of two-dimensional topological insulators (TIs) is described by the Bernevig-Hughes-Zhang (BHZ) model, where the realization of a topological or a normal insulating phase depends on the Dirac mass being negative or positive, respectively. We solve the BHZ model for a mass domain configuration, analyzing the effects on the edge modes of a finite Dirac mass in the normal insulating region (soft-wall boundary condition). We show that at a boundary between a TI and a normal insulator (NI), the Dirac point of the edge states appearing at the interface strongly depends on the ratio between the Dirac masses in the two regions. We also consider the case of multiple boundaries such as NI/TI/NI, TI/NI/TI and NI/TI/NI/TI.Comment: 11 pages, 15 figure