I have a graphene bilayer structure whose capacitance I wish to find out using QuantumATK. I have made it into a device configuration with 6 layers and L-R electrodes I have followed the steps as given in the tutorial

https://docs.quantumatk.com/tutorials/atomic-scale_capacitance/atomic-scale_capacitance.html for the calculation of atomic scale capacitance. And I am facing some problems and some queries that I wished to clarify.

1. When the capacitance is being calculated using the induced charge method (Mulliken charge), I am getting a very high value of the plate distance d. When viewed from the structure, the actual value of plate distance d = 10.26 Angstrom whereas from the induced charge analysis, I am getting the value of plate distance as 160.84 Angstrom. I am aware that the value of distance is calculated from the capacitance value. But even then why am I getting such a value of capacitance? Attached along with are theinput .py file along with log file (mulliken_analysis.py and mulliken_analysis.log) and the induced charge plot (Device_Graphene_mulliken.png).

2. When I am calculating the value of capacitance using the electrostatic energy analysis, I am getting a plot of Induced density at 1 volt (Device_Graphene_induced.png file) which has two maximas (dips). I wanted to know whether the distance between the dips is same as that of the plate distance d?