Mott Transition in the Hubbard Model on Anisotropic Honeycomb Lattice with Implications for Strained Graphene: Gutzwiller Variational Study

The modification of interatomic distances due to high pressure leads to exotic phenomena, including metallicity, superconductivity and magnetism, observed in materials not showing such properties in normal conditions. In two-dimensional crystals, such as graphene, atomic bond lengths can be modified by more than 10 percent by applying in-plane strain, i.e., without generating high pressure in the bulk. In this work, we study the strain-induced Mott transition on a honeycomb lattice by using computationally inexpensive techniques, including the Gutzwiller Wave Function (GWF) and different variants of Gutzwiller Approximation (GA), obtaining the lower and upper bounds for the critical Hubbard repulsion (U) of electrons. For uniaxial strain in the armchair direction, the band gap is absent, and electron correlations play a dominant role. A significant reduction in the critical Hubbard U is predicted. Model considerations are mapped onto the tight-binding Hamiltonian for monolayer graphene by the auxiliary Su-Schrieffer-Heeger model for acoustic phonons, assuming zero stress in the direction perpendicular to the strain applied. Our results suggest that graphene, although staying in the semimetallic phase even for extremely high uniaxial strains, may show measurable signatures of electron correlations, such as the band narrowing and the reduction in double occupancies.