README

kp_tblg

Generates k-dot-p Hamiltonians based on
ab initio models of twisted bilayer graphene.
A C++ (python compatible) and a MATLAB version included.

For reference of the k-dot-p model, see (and please cite) the following two works:

"Minimal model for low-energy electronic states of twisted bilayer graphene."
Stephen Carr, Shiang Fang, Ziyan Zhu, Efthimios Kaxiras
https://arxiv.org/abs/1901.03420

"Angle-dependent ab initio low-energy hamiltonians for a relaxed twisted bilayer graphene heterostructure."
Shiang Fang, Stephen Carr, Ziyan Zhu, Daniel Massatt, Efthimios Kaxiras
https://arxiv.org/abs/1908.00058

The 8-band and 5-band reduced models can be found in (check the README!):
./MATLAB_vers/wannierization/ 

For reference of the 8-band and 5-band models, see (and please cite):

"Derivation of Wannier orbitals and minimal-basis tight-binding hamiltonians for twisted bilayer graphene: a first-principles approach."
Stephen Carr, Shiang Fang, Hoi Chun Po, Ashvin Vishwanath, Efthimios Kaxiras
https://arxiv.org/abs/1907.06282

=====

CONTACT:
stcarr.nj(at)gmail.com
Please let me know of any bugs, or if you have a feature request.

=====

NOTE ON DATASET:

There is currently one dataset of k-dot-p terms.
The terms are organized in "shells" of equal momenta magnitude,
organized via hexagonal Miller (crystallographic) indices.
This includes 10 "shells" of coupling for both interlayer and
intralayer coupling, and 1 "shell" of k-dependent interlayer terms.

The dataset "data/full_relax_kp_01-06-2019.dat" contains terms
for theta between 5.08 and 0.416 (degrees).

We are preparing additional datasets for smaller angles, but these will
be relased seperately. They require many more "shells" for accuracy in
the low-energy electronic states.

=====

NOTE ON UNITS:

Our model uses energy units of eV and length units of Angstroms.

=====

INSTALLATION:

This code requires a c++ compiler and access to the header-only libraries:
Eigen (required)
pybind11 (optional)

A python-readable library can be generated if the option is turned 
on in the CMakeLists.txt file in this directory.
If making a python library, you have to pull the pybind11 repo by hand:

 git submodule update --init --recursive

To compile:

 mkdir build
 cd build
 cmake ..
 make

[python note: if pybind11 is finding python.3.x instead of python.2.x, try:
 
 cmake -DPYTHON_EXECUTABLE:FILEPATH=/usr/bin/python ..

instead of < cmake .. >, this will usually force it to use python (not python3)]

=====

USE:

The class Kp_tblg_construct can be used to create accurate k-dot-p 
Hamiltonians including relaxation of the layers.

(1) Initalize:		  Kp_tblg_construct kp = Kp_tblg_construct(); 
(2) Load data file: 	  kp.loadFiles("data/full_relax_kp_01-06-2019.dat");
(3) Set twisting angle:   kp.setTwist(1); // in degrees
(4) Set sparsity pattern: kp.prepare();
(5) Access Hamiltonian:	  Eigen::MatrixXcd H = kp.getH(k); // k must be an Eigen::Vector2d

Then use the hamiltonain H to e.g. compute bandstructure by diagonalization

[In python, the usage is identical, except replace Eigen datatypes with their numpy counterparts]

=====

TESTING (C++):

Execute the application app/kp_tblg_test:

./app/kp_tblg_test

This requires the data file (data/full_relax_kp-01-06-2019.dat) to be
in the same folder as where the executable is called from.
Note that the Eigen diagonalization routine is sometimes very slow, so this may take a while

=====

TESTING (PYTHON):
 
For testing the python code, try using:

pysrc/generate_H_example.py
pysrc/bandcalc.py

Both the data-file (data/full_relax_kp_01-06-2019.dat)
and the tblg_kpy.so library must be in your current 
working directory when these scripts are called.