Band Structures of Bilayer Graphene Superlattices

We formulate a low energy effective Hamiltonian to study superlattices in bilayer graphene (BLG) using a minimal model which supports quadratic band touching points. We show that a one dimensional (1D) periodic modulation of the chemical potential or the electric field perpendicular to the layers leads to the generation of zero-energy anisotropic massless Dirac fermions and finite energy Dirac points with tunable velocities. The electric field superlattice maps onto a coupled chain model comprised of ``topological'' edge modes. 2D superlattice modulations are shown to lead to gaps on the mini-Brillouin zone boundary but do not, for certain symmetries, gap out the quadratic band touching point. Such potential variations, induced by impurities and rippling in biased BLG, could lead to subgap modes which are argued to be relevant to understanding transport measurements.

Laterally structured ripple and square phases with one and two dimensional thickness modulations in a model bilayer system

Molecular dynamics simulations of bilayers in a surfactant/co-surfactant/water system with explicit solvent molecules show formation of topologically distinct gel phases depending upon the bilayer composition. At low temperatures, the bilayers transform from the tilted gel phase, L beta', to the one dimensional (1D) rippled, P beta' phase as the surfactant concentration is increased. More interestingly, we observe a two dimensional (2D) square phase at higher surfactant concentration which, upon heating, transforms to the gel L beta' phase. The thickness modulations in the 1D rippled and square phases are asymmetric in two surfactant leaflets and the bilayer thickness varies by a factor of similar to 2 between maximum and minimum. The 1D ripple consists of a thinner interdigitated region of smaller extent alternating with a thicker non-interdigitated region. The 2D ripple phase is made up of two superimposed square lattices of maximum and minimum thicknesses with molecules of high tilt forming a square lattice translated from the lattice formed with the thickness minima. Using Voronoi diagrams we analyze the intricate interplay between the area-per-head-group, height modulations and chain tilt for the different ripple symmetries. Our simulations indicate that composition plays an important role in controlling the formation of low temperature gel phase symmetries and rippling accommodates the increased area-per-head-group of the surfactant molecules.