We consider the tight-binding model of graphene with slowly spatially varying hopping functions. We develop a low energy approximation as a derivative expansion in a Dirac spinor that is perturbative in the hopping function deformation. The leading description is the Dirac equation in flat 2+1-d spacetime with (strain-)gauge field. Prior work considered subleading corrections written as non-trivial frame and spin connection terms. We previously argued that such corrections cannot be considered consistently without taking all the terms at the same order of approximation, which due to the unconventional power counting originating from the large gauge field, involve also higher covariant derivative terms. Here we confirm this, explicitly computing subleading terms. To the order we explore, the theory is elegantly determined by the gauge field and frame, both given by the hopping functions, the torsion free spin connection of the frame, together with coefficients for the higher derivative terms derived from lattice invariants. For the first time we compute the metric that the Dirac field sees - the `electrometric' - to quadratic order in the deformation allowing us to describe the subleading corrections to the dispersion relation for inhomogeneous deformations originating from corrections to the frame. Focussing on in-plane inhomogeneous strain, we use a simple model to relate the hopping functions to the strain field, finding the electrometric becomes curved at this quadratic order. Thus this lattice model yields an effective analog gravity description as a curved space Dirac theory, with large magnetic field, and Lorentz violating higher covariant derivative terms. We check this by comparison to numerical diagonalization. From this we conjecture a form for the effective theory for monolayer graphene in terms of the strain tensor, consistent up to quadratic order in the deformation.

39 pages, 7 figures, downloadable Mathematica notebook, abridged abstract; v2. typos corrected, references added