Tracing boron diffusion into a textured silicon solar cell using electron beam induced current in a scanning transmission electron microscope

Resumen Técnico: Rastreo de la Difusión de Boro en Células Solares de Silicio Texturizadas mediante STEM-EBIC

Planteamiento del Problema
La caracterización de la difusión de dopantes en dispositivos optoelectrónicos modernos, particularmente en aquellos con estructuras superficiales fuertemente texturizadas como las estructuras tipo pirámide en células solares de silicio, presenta desafíos significativos. Las técnicas macroscópicas como la capacitancia-voltaje electroquímica, SIMS o la espectrometría de masas por tiempo de vuelo de plasma carecen de la resolución lateral requerida para estructuras heterogéneas. Por el contrario, los métodos de alta resolución como la tomografía de sonda atómica o nano SIMS detectan átomos de dopantes independientemente de su estado de ionización, fallando en capturar las especies eléctricamente activas críticas para el rendimiento del dispositivo. Si bien la Microscopía de Resistencia de Dispersión de Escaneo (SSRM) y la Microscopía de Capacitancia de Escaneo ofrecen una alta resolución, sufren de dificultades de calibración debido al desgaste de la punta. La Corriente Inducida por Haz de Electrones (EBIC) ofrece alta resolución espacial y sensibilidad a los dopantes activos, pero se ve complicada por los fuertes efectos de recombinación superficial, lo que dificulta la interpretación cuantitativa en geometrías complejas. Además, un modelo cuantitativo específicamente aplicable a la geometría de Microscopía Electrónica de Transmisión de Barrido (STEM) —que utiliza muestras delgadas y transparentes a los electrones— se ha establecido recientemente.

Metodología
Los autores emplearon un marco de trabajo autoconstante que combina mediciones experimentales de STEM-EBIC con simulaciones de método de elementos finitos (FEM) de la difusión de boro y el transporte de portadores de carga.

• Preparación de la Muestra: Se preparó una célula solar de silicio texturizada con un emisor de boro difundido. Se utilizó un haz de iones focalizados (FIB) para extraer una lámina transparente a los electrones de la zona cubierta por aluminio entre dos pirámides. Se realizaron cortes verticales durante la preparación para evitar cortocircuitos por material redepositado.
• Caracterización Experimental:
  • STEM-EBIC: Los experimentos se llevaron a cabo en un E-TEM FEI Titan 80-300 E a 300 kV. Un haz de electrones (42 pA) escaneó la muestra y se midió la corriente de cortocircuito para generar mapas 2D de EBIC.
  • SSRM: Las mediciones de SSRM transversales se realizaron en el borde de una muestra seccionada y pulida para proporcionar una estimación aproximada de la posición de la unión química.
  • Referencia: Se procesó simultáneamente una muestra de referencia plana para determinar la resistencia de hoja mediante mediciones de cuatro puntas (4PP).
• Marco de Simulación:
  • Simulación de Difusión: Utilizando Sentaurus Process, los autores modelaron el proceso de difusión de boro (deposición, conducción, post-oxidación, enfriamiento). Consideraron el crecimiento de una capa de vidrio de boro y la segregación de boro en el silicio y el óxido. El modelo fue calibrado utilizando la resistencia de hoja planar y la posición de la unión química derivada de los datos experimentales. Se aplicó un factor de escala de difusión mejorada por oxígeno (OED) para coincidir con las observaciones experimentales.
  • Simulación de EBIC: Se implementó un modelo de deriva-difusión 2D (omitiendo la dirección de propagación del haz) en COMSOL Multiphysics. El modelo resolvió las ecuaciones de van-Roosbroeck (ecuaciones de Poisson y de continuidad) para el potencial electrostático y las concentraciones de electrones y huecos.
  • Modelado de Recombinación: Para dar cuenta de la recombinación superficial en la delgada lámina, se empleó un modelo de recombinación efectiva. Este modelo proyecta el sistema a lo largo de la dirección del haz de electrones, tratando la recombinación superficial y volumétrica mediante tiempos de vida efectivos ($\tau_h$ y $\tau_e$) dentro de un marco de Shockley-Read-Hall. Este enfoque permite la comparación de mapas experimentales 2D completos con las simulaciones.

Resultados Clave

• Observaciones Experimentales: El mapa STEM-EBIC reveló una señal notablemente por debajo del contacto de aluminio, lo que refleja la textura superficial pero apareciendo aplanada. El análisis de los perfiles verticales arrojó longitudes de difusión efectivas de 63 nm (lado p) y 103 nm (lado n), valores significativamente menores que las expectativas de volumen debido a la recombinación superficial. Las mediciones de SSRM indicaron una profundidad de la unión química de aproximadamente 500–550 nm por debajo del contacto de aluminio.
• Acuerdo de la Simulación: La distribución de boro simulada, calibrada para coincidir con la resistencia de hoja experimental y la posición de la unión, mostró un perfil retrógrado con una concentración máxima de ~0.12 µm por debajo de la superficie.
• Comparación de EBIC: Hubo un buen acuerdo entre los mapas 2D experimentales y simulados de EBIC respecto a la forma de la señal y la posición del máximo de EBIC. Las longitudes de decaimiento de las colas de la señal en las regiones neutras también coincidieron, validando el modelo de recombinación efectiva.
• Discrepancias: Se observó un desplazamiento de aproximadamente 80 nm entre la posición del máximo campo eléctrico (unión química) y el máximo de EBIC. Además, se notaron desviaciones dentro de la región de carga espacial, particularmente en el lado n. Los autores atribuyen estas discrepancias a posibles cargas superficiales (por ejemplo, de la preparación por FIB), términos de recombinación simplificados o heterogeneidad 3D no capturada por el modelo 2D.

Contribuciones Clave

1. Primer Modelado Numérico de Mapas 2D de STEM-EBIC: El artículo demuestra, por primera vez, que los mapas completos de STEM-EBIC 2D de sistemas heterogéneos y texturizados pueden modelarse numéricamente. Esto se logra mediante el uso de un modelo de recombinación efectiva que proyecta la geometría de la muestra delgada 3D en un marco 2D.
2. Marco de Trabajo Autoconstante: Los autores establecieron un flujo de trabajo donde las simulaciones de difusión proporcionan la entrada de distribución de dopantes para las simulaciones de EBIC, que luego se comparan directamente con los mapas experimentales. Esto permite rastrear distribuciones de dopantes debajo de texturas superficiales complejas.
3. Validación de Modelos Efectivos: El estudio confirma que los modelos efectivos que consideran la recombinación volumétrica y superficial pueden describir con éxito las señales de EBIC en regiones de semiconductores neutros en muestras delgadas, permitiendo comparaciones cuantitativas entre experimento y simulación.

Significancia y Afirmaciones
El artículo afirma que este enfoque abre "posibilidades sin precedentes" para observar la dinámica de portadores de exceso, campos eléctricos débiles y distribuciones de dopantes en dispositivos funcionales a nanoescala. Los autores enfatizan que, si bien el modelo actual muestra un excelente acuerdo en las regiones neutras, las discrepancias en la región de carga espacial sugieren que los parámetros de recombinación efectiva pueden depender del campo eléctrico.

Los autores sugieren modestamente que, una vez establecido un modelo para la dependencia del campo eléctrico, la técnica STEM-EBIC podría servir como un método altamente sensible para medir campos eléctricos débiles, superando potencialmente los límites de sensibilidad de técnicas complementarias como la holografía electrónica o el 4D-STEM, particularmente en sustratos de bajo dopaje. El trabajo destaca el valor de obtener retroalimentación sobre las distribuciones de dopantes en geometrías heterogéneas a la escala nanométrica, cerrando la brecha entre la imagen de alta resolución y el análisis cuantitativo de propiedades electrónicas.