Densidad de estados para el grafeno

No estoy seguro si es demasiado tarde para esta respuesta, pero por si acaso.

La densidad de estados se define como$D(E)=\int_{1st BZ}\delta(E-E(\mathbf{k}) d\mathbf{k}$.

Para resultados analíticos, el siguiente artículo fue el primero en calcular el DOS en una estructura de panal con salto de vecino más cercano (cuya relación de dispersión es la que mencionas): http://journals.aps.org/pr/abstract/10.1103 /PhysRev.89.662 .

Para resultados numéricos, simplemente aplique la definición de$D(E)$. Explícitamente, debe crear dos matrices que contengan valores de$k_x$y$k_y$(digamos 100 para cada uno pasando de sus valores mínimos a sus valores máximos). Calcula sus energías asociadas. Cortar la 1.ª BZ en trozos del mismo tamaño$\mathbf{k}$regiones. Finalmente calcule el número de estados en cada área (suma 1 por cada$E(\mathbf{k})$dentro de un área determinada delimitada por los valores de$E(\mathbf{k})$evaluados en los límites de cada uno de los$\mathbf{k}$regiones en la 1ª BZ). Obviamente, las dimensiones de su$\mathbf{k}$área debe ser mayor que la diferencia entre consecutivos$k_x$y$k_y$. Finalmente normaliza tu resultado.