I have created a cif file containing two different organic 2D layers with large vacuum in the z-axis. Previously I've calculated the the band gap of each layer separately using different cif files. But now I need to calculate the band gap of each of these two layers, while they are on top of each other (to see the change). Now if I do scf calculation for theses two layers in a single cif, it will show me the HOMO-LUMO of the entire system and not each layer.
So I wonder how I should approach that, so that it will give me HOMO-LUMO of each layer separately. Thanks!