Hi. I am also facing the problem when building bilayer graphene. I tried as what you said, to start by graphite in database, and then change the coordinates of the second layer very close to the first layer. But, it seems that this does not work when calculating the Band strcutrue. Could you please give me hand?