Un enfoque DFT-1/2 preciso para estados de defecto poco profundos: cálculo eficiente de las energías de enlace de donores en silicio

Por qué importan pequeños ajustes en el silicio

Cada chip y cada célula solar depende de impurezas cuidadosamente añadidas, o “dopantes”, que controlan cuán fácilmente fluye la electricidad. Para algunas tecnologías de vanguardia —desde transistores ultrarrápidos hasta qubits basados en donores— necesitamos saber con precisión cuán fuertemente queda ligado un electrón extra a un átomo dopante dentro de un cristal semiconductor. Este artículo presenta una forma más rápida y práctica de calcular esa energía de enlace con alta precisión, especialmente para dopantes en silicio, la columna vertebral de la electrónica moderna.

Átomos que donan electrones extra

En el silicio puro, los átomos comparten electrones en enlaces ordenados y repetitivos, y el material conduce poco a temperatura ambiente. Añadir una pequeña cantidad de un elemento del grupo V como fósforo, arsénico, antimonio o bismuto aporta un electrón extra por dopante. Ese electrón adicional no queda libre; en cambio, se sitúa en una nube similar a la de un átomo de hidrógeno, débilmente ligada al dopante y al silicio circundante. La fuerza de ese enlace —la energía de enlace del donante— determina qué tan fácilmente puede liberarse el electrón para conducir corriente o participar en operaciones cuánticas. Medir estas energías en el laboratorio está bien establecido, pero predecirlas de forma fiable desde cálculos de primer principio ha resultado difícil y costoso.

Por qué los cálculos estándar se quedan cortos

Los modelos computacionales basados en la teoría del funcional de la densidad (DFT) son la herramienta de trabajo del diseño de materiales, pero tienden a subestimar cuán localizados están los electrones y a colocar incorrectamente los bordes de las bandas de energía en semiconductores. Para donores poco profundos, cuyas nubes electrónicas se extienden sobre decenas de átomos, esto significa que DFT suele predecir energías de enlace mucho menores de las reales. Métodos más avanzados, como los funcionales híbridos y los cálculos GW, pueden corregir estos problemas, pero a un coste computacional enorme, especialmente cuando se necesitan cajas de simulación grandes con miles de átomos para capturar el estado extendido del donante. Enfoques previos “en tándem” tuvieron que mezclar distintos niveles de teoría en diferentes tamaños de celda y luego ensamblar los resultados, haciendo el flujo de trabajo complejo y dependiente del sistema.

Una corrección simple con gran beneficio

Los autores se basan en una técnica llamada DFT‑1/2, que añade una corrección aproximada de autoenergía directamente en DFT estándar. En términos prácticos, modifican ligeramente el potencial efectivo de átomos seleccionados al conceptualmente eliminar medio electrón de orbitales atómicos concretos. Primero aplican esta corrección al silicio en bloque para que la brecha de banda calculada concuerde mucho mejor con el experimento y proporcione una referencia fiable para la banda de conducción. A continuación, examinan el carácter electrónico del estado del donante y encuentran que, para todos los dopantes del grupo V, está dominado por el orbital s del átomo dopante. Luego aplican una corrección de medio electrón adaptada a ese orbital y afinan un único radio de corte para maximizar la separación entre el nivel del donante y el estado de conducción vacío más cercano. Es importante que esta corrección optimizada siga siendo válida cuando se amplía la caja de simulación, por lo que puede reutilizarse en superceldas que contienen hasta miles de átomos.

Qué tan bien funciona el método

Con esta corrección en dos pasos —primero al silicio anfitrión y luego al propio donante— el método produce energías de enlace de donantes que coinciden estrechamente con los valores experimentales. Para arsénico en silicio, la energía predicha difiere del experimento en solo 0,3 milielectrónvoltios, un acuerdo prácticamente perfecto y comparable con cálculos híbridos mucho más costosos. Para antimonio y fósforo, los errores son de aproximadamente 5 y 8 milielectrónvoltios, respectivamente, una mejora notable frente al DFT sin corregir. Para el bismuto, un dopante muy pesado, los autores también incluyen acoplamiento espín‑órbita, un efecto relativista que modifica ligeramente los niveles de energía. Esto reduce la energía de enlace calculada a un valor dentro de unos 5 milielectrónvoltios del experimento y pone de manifiesto física que métodos anteriores, más exigentes, habían pasado por alto. Para mostrar que el enfoque no está limitado al silicio, aplican con éxito el mismo flujo de trabajo a un donante de hidrógeno en óxido de zinc, reproduciendo de nuevo las energías de enlace medidas con una discrepancia de solo unos pocos milielectrónvoltios.

Una herramienta práctica para diseñar chips del futuro

Para no especialistas, el mensaje clave es que los autores han creado una receta que conserva el bajo coste y la simplicidad de los cálculos DFT estándar mientras alcanza la precisión de métodos mucho más pesados. Corrigiendo sistemáticamente tanto la estructura de bandas global del material anfitrión como el entorno local alrededor del dopante, su protocolo DFT‑1/2 ofrece energías de enlace de donantes fiables en celdas de simulación muy grandes. Esto lo convierte en una herramienta potente y general para estudiar los dopantes que controlan la electrónica cotidiana y los dispositivos cuánticos emergentes, ayudando a los ingenieros a diseñar materiales en los que las impurezas individuales se comporten exactamente como se pretende.

Cita: Claes, J., Partoens, B., Lamoen, D. et al. An accurate DFT-1/2 approach for shallow defect states: efficient calculation of donor binding energies in silicon. npj Comput Mater 12, 153 (2026). https://doi.org/10.1038/s41524-026-02003-2

Palabras clave: donores poco profundos, dopantes de silicio, teoría del funcional de la densidad, materiales cuánticos, defectos en semiconductores