nueva página del texto (beta)
Español (pdf)
Artículo en XML
Referencias del artículo
Traducción automática
Enviar artículo por email
Citado por SciELO 
Similares en
SciELO 
Permalink
PACS: 78.67 Bf; 73.20 Mf; 63.2.Pw; 78.67.-n; 42.25.Bs
1. Introducción
La nanoplasmónica se define como el estudio de las excitaciones en superficies metálicas localizadas en nanoestructuras. Es un campo de investigación que ha tomado fuerza desde la primera década del presente siglo 1, es de interés su existencia en coloides cuya componente sólida posee simetría esférica. El estudio de las propiedades ópticas de una partícula de forma esférica y estructuras esféricas más complejas comienza en el principio del siglo XX. En esta época, G. Mie 2 estudia la difracción de una onda plana electromagnética monocromática por una esfera y encuentra una solución exacta a este problema, también aplica la solución al estudio óptico de coloides formados por partículas metálicas. Al mismo tiempo P. Debye resuelve un problema semejante e innova introduciendo los llamados potenciales de Debye 3. Poco antes J.C. Maxwell-Garnett 4 explicó los colores de coloides metálicos.
Recientemente publicamos un trabajo asociado a la difracción de luz por esferas dieléctricas 5 y excitaciones producidas en estas. La teoría general de la difracción de luz en esferas es válida tanto para esferas dieléctricas como met álicas, por lo cual la teoría básica que aparece en ese trabajo será utilizada aquí, referenciándonos a las ecuaciones que allí se presentan. También gran número de referencias bibliográficas que allí se muestran son de nuestro interés. Dos magníficos libros en el tema de agregados metálicos son el de Kreibig and Vollmer 6 en 1995 y otro editado por Mishchenko et al. 7 en el 2000. Un amplio resumen de resultados para partículas metálicas se puede encontrar en el trabajo de Mulvaney 8 y en el crecimiento de nano-partículas de plata en los trabajos de Slistan-Grijalva et al., 9,10. En 2007 Derkachova y Kolwas 11 obtienen las frecuencias complejas de plasmones de Superficie (SP) en función del radio de la esfera metálica. Ancey et al., en 2009 calculan la sección eficaz de esparcimiento y de absorción respecto a la energía de excitación 12. Es conocido que las frecuencias asociadas a los modos SP dependen del tamaño, la forma, y composición de la partícula, así como del medio en que se encuentra. Además, no es necesario usar algún medio de acoplamiento para excitarlos 13.
Para el desarrollo de este estudio se utilizan los resultados del trabajo 5 que introducen los armónicos esféricos para resolver la ecuación de la onda esférica escalar con los cuales se resuelven las ecuaciones de Maxwell mediante un desarrollo multipolar donde se usan los vectores armónicos esféricos, comenzando con el desarrollo de una onda plana, para posteriormente estudiar el problema del esparcimiento de la esfera. Se utiliza la notación del libro Classical Electrodynamics por D. Jackson 14, En el capítulo 10 trata de forma breve el problema de Mie. Los vectores armónicos esféricos son la parte medular de la expansión multipolar de un campo electromagnético.
En la siguiente sección se presentan los resultados numéricos, criterio de convergencia, el campo cercano, la sección de extinción de una esfera metálica. En detalle se estudia el caso de esferas metálicas y el plasmón de superficie. En la última sección vertemos nuestras conclusiones.
2. Difracción de Ondas Electromagnéticas por Esferas Metálicas
Se dan las ecuaciones de partida para calcular el campo electromagnético de una onda plana difractada por una esfera, metálica de radio
Si se asume que la dependencia temporal de los campos es armónica, es decir,
y
La representación multipolar de los campos está dada por una serie infinita de los denominados esféricos armónicos vectoriales,
donde l = 1, 2, ,…, ꝏ, m = -l, …, -1, 0, +1, …, +l, ε, es la constante dieléctrica y
Donde
La sección eficaz de esparcimiento (SES) total de la esfera 15,16, es:
La sección eficaz de absorción es
La sección eficaz de extinción (SEE) es la suma de las secciones eficaces de esparcimiento y de la absorción
Existen modos que pueden ser plasmones de superficie para el caso de la esfera metálica se manifiestan en un máximo o máximos en la sección eficaz de esparcimiento. Cuando no hay absorción en la esfera, la sección eficaz de extinción es igual a la sección eficaz de esparcimiento.
3. Cálculos y resultados
Para obtener resultados numéricos de la sección eficaz de esparcimiento (4) o del campo eléctrico (2), las sumatorias infinitas se truncan hasta un valor L, tal que los datos calculados con L-1 y con L no difieran en 1%. Una regla gruesa es que valores pequeños del radio de la esfera requiere L pequeños, mientras que radios grandes requieren L grandes. Así que l, en (1) y (2), tomará valores desde 1 hasta L.
En esta sección presentamos resultados numéricos de la difracción de las ondas electromagnéticas por esferas metálicas con función dieléctrica tipo Drude. Los radios fueron variados según fuera procedente.
El libro de Bohren y Huffman 16 tiene un apéndice con un código en FORTRAN para el caso de una esfera con el que se puede calcular la sección eficaz de extinción, de esparcimiento, así como la intensidad en la aproximación de campo lejano, pero no se puede calcular el campo cercano, así que este código fue modificado para la incluir el cálculo del campo cercano. Este código modificado fue empleado para obtener los resultados del campo cercano y sección eficaz debido a la difracción por esferas dieléctricas 5. Es pertinente mencionar que actualmente existen programas orientados a realizar estos cálculos en el rango nanométrico tales como el MieLab o Scantlay, basados en la teoría de Mie.
3.1. Plasmones de Superficie en Esferas Metálicas
En esta sección se estudiará la difracción de luz por esferas metálicas, una referencia útil es 17. Para poder comparar con algunos resultados mostrados en la Ref. 12, se presenta la constante dieléctrica tipo-Drude, Ec. (2) de la misma referencia, que ellos usaron en su artículo,
Los parámetros empleados son la constante dieléctrica de fondo
Las frecuencias de los modos soportados por esferas metálicas con respuesta tipo Drude en el caso no-retardado están dadas por una sencilla expresión, ver Engelman and Ruppin en Ref. 18,
Estas frecuencias son aplicables para radios mucho menores que la longitud de onda. En el caso retardado las frecuencias de los modos son complejas, lo cual refleja el hecho de que los plasmones de superficie tienen una vida media finita debido a su decaimiento radiativo y se calculan haciendo cero el denominador de la Ec. (20) de la Ref. 5. Estas frecuencias fueron calculadas en la Ref. 11.
Si el radio de esfera es muy pequeño comparado con la longitud de onda (
Donde
Sustituyendo (10) en la Ec. (6) la SEE es
La derivación de (11) fue hecha por Genzel and Martin 19. De (10) podemos obtener el valor de la frecuencia para la cual (10) tiene un máximo
Primero pasaremos a estudiar las convergencias de las soluciones, para esto en la Fig. 1 se muestra la intensidad del campo cercano, el módulo cuadrado del campo eléctrico, calculado en el plano z-y para a) L = 1, b) L = 2, c) L = 3 y d) L = 4. La esfera metálica tiene un radio de 10 nm. La onda electromagnética incidente está linealmente polarizada y se propaga en la dirección z. La longitud de onda es de 633 nm, la cual corresponde a una frecuencia de
Figura 1. La intensidad de campo cercano proyectado en el plano z-y debido a una esfera metálica de radio 10 nm, para diferentes valores de L: a) L = 1, b) L = 2, c) L = 3 y d) L = 4. La flecha indica la dirección de la onda incidente, λ = 633 nm, y γ = 0.001. Se aprecian ligeras diferencias entre las gráficas calculadas con L = 1 y L = 2. Mientras que los datos obtenidos con L = 3 y L = 4 son iguales.
De las figuras se puede apreciar que la convergencia se alcanzó para L = 4, esto se verificó numéricamente de acuerdo al criterio establecido en la sección anterior.
En las imágenes se puede observar que la esfera presenta un comportamiento tipo dipolar con orientación en el eje-y, que coincide con la dirección de oscilación del campo eléctrico incidente. Este es un comportamiento típico de la difracción de luz por nano-esferas metálicas. También se observa que los dos máximos de intensidad del campo se localizan muy cerca de la superficie, típicamente del orden de a/10.
En la Fig. 2 se muestra la sección eficaz de
esparcimiento (SES) como función de la frecuencia dividida entre la frecuencia
de plasma para una esfera metálica de 10 nm de radio. La función dieléctrica de
la esfera metálica se supone tipo Drude. La curva de color rojo corresponde a γ
= 0, la de color negro a γ = 0.001 y la de color azul a la aproximación dada por
(10), calculada para este último valor de γ, en la Fig. 2. La SES despliega tres máximos, en el caso exacto. Para poder
resolver uno de los picos la resolución espectral fue de
Figure 2. La sección eficaz de esparcimiento total como función de la frecuencia adimensional ω/ ωp , para una esfera metálica de radio a = 10 nm rodeada de aire ε1 = n1 = 1.0, la curva de color rojo corresponde a γ = 0.00, la negra a γ = 0.01, la curva mostrada en color azul fue calculada usando la forma aproximada (10). Junto a esta última curva se muestra la expresión para la frecuencia de resonancia, el máximo y el ancho de las SES.
El valor asintótico de la frecuencia cuando l tiende a infinito en la fórmula de Ruppin es la frecuencia del plasmón de superficie en una interfaz plana meta-aire o sea
La SEE total cuando γ = 0 tiene un valor máximo más grande que cuando γ = 0.001, pero una anchura menor, esto se debe a la absorción del medio.
Una consecuencia de la excitación de un plasmón de superficie es el incremento del campo eléctrico muy cerca de la superficie. Para observar este efecto se estudia el campo cercano de la esfera metálica a la frecuencia de resonancia.
En la Fig. 3 se muestra los resultados para SEE, enfatizando la variación del radio de la esfera metálica: a) 10 - 30 nm, b) 40 -60 nm, c) 70- 90 nm, d) 100 -120 nm. Se observa la tendencia a obtener varias excitaciones de SP, conforme el radio aumenta.
Figura 3. Se muestra los resultados graficados para SEE, al variar el radio de la esfera metálica desde 10 nm a 120 nm: a) 10-30 nm, b) 40-60 nm, c) 70-90 nm, d) 100-120 nm. Las imágenes ilustran claramente la tendencia a una correlación entre el número de excitaciones de SP y el incremento del radio de la esfera bajo estudio.
En la Fig. 4 se muestra la intensidad del campo cercano calculada en el plano a) z-y y a) z - y para L = 4, a = 10 nm, γ = 0.01 y a la frecuencia de resonancia ω1,1 = 0.55918ωp . Se puede observar de nuevo el comportamiento dipolar anteriormente descrito. Los máximos de intensidad del campo eléctrico están muy cerca de la superficie, pero son aproximadamente 1000 veces mayor que los máximos mostrado en la Fig. 1, además la intensidad es notable a distancias del orden del radio, al menos a lo largo del eje-y.
Figura 4. La intensidad del campo cercano calculada en el plano a) z - y y b) z - x producida por una esfera metálica con valores de radio a = 10 nm, γ = 0.01, L = 4 y a la frecuencia de resonancia ω1,1 = 0.55918ωp , la cual corresponde a una longitud de onda γ = 221.7660 n.
En la Fig. 5 se muestra el módulo cuadrado de la componente φ del campo eléctrico como función del ángulo de esparcimiento, para a = 10 nm y γ = 0.01, la curva de color rojo corresponde al caso resonante y la negra al caso no resonante.
Figura 5. La intensidad del campo cercano calculada en el plano a) z - y y b) z - x producida por una esfera metálica con valores de radio a = 10 nm, γ = 0.01, L = 4 y a la frecuencia de resonancia ω1,1 = 0.55918ωp , la cual corresponde a una longitud de onda γ = 221.7660 n.
Observamos que la magnitud del campo cuando hay resonancia es mayor que cuando no la hay, la distribución angular es la misma. Sólo hay un mínimo alrededor de 90°.
La intensidad muestra un comportamiento cuasi-lambertiano sin ninguna otra característica importante.
En el caso retardado, varios autores como Ruppin 18, Martinos 17, Ancey et al.12 reportaron varias resonancias, así que para observarlas escogemos un radio mayor con el objetivo de reproducir los resultados de la Fig. 1 a 1a mostrada en 12.
En la Fig. 6 se muestra la sección eficaz de
esparcimiento total como función de
Figura 6. Sección eficaz de esparcimiento total como función de la frecuencia ωa/c. El radio de la esfera es a = 2πc/ωp . La curva solida de color negro se calculó con γ = 0.00 y la curva en color rojo con γ = 0.01. La SEE despliega 9 picos asociados con la excitación de plasmones de superficie esféricos.
En la Fig. 7 se incluye el símbolo del coeficiente de la onda parcial relacionado con cada resonancia. Para lograr asociar a cada resonancia una onda parcial se encontró el mínimo del denominador, que debe ser aproximadamente cero, asociado con cada una de ellas. Los modos excitados son del tipo TM solamente.
Figura 7. El módulo de la componente φ del campo eléctrico esparcido
En la Tabla I se muestra el valor de la frecuencia reducida obtenida de la SES y los valores dados por la fórmula de Ruppin. Los valores de las frecuencias, son aproximadamente iguales, lo cual establece que la aproximación no-retardada es válida para radio de esferas mucho menores que la longitud de onda.
Tabla I. Valores de la frecuencia reducida obtenida de la SEE y los valores dados por la fórmula de Ruppin.
En la Tabla II se muestra el valor de la frecuencia reducida obtenida de la SEE total y los valores dados por la fórmula de Ruppin.
Tabla II. Se muestran valores de la frecuencia reducida obtenida de la SES total y los valores dados por la fórmula de Ruppin.
En este caso no hay coincidencia alguna entre las frecuencias de los modos, obtenida de la SEE de la Fig. 6, con aquellas dadas por la fórmula cuasi-estática de Ruppin.
Ahora pasaremos a estudiar la luz esparcida y el campo cercano producidas por la esfera metálica a dos frecuencias distintas, una resonante ωa/c = 3.416 y la otra no resonante ωa/c = 2.180, ambas marcadas con flechas de color azul en la figura anterior.
En la Fig. 7 se muestra
La intensidad del campo lejano presenta máximos y mínimos relacionado con el radio de la esfera. A mayor tamaño mayor número de máximos. De estas curvas no es posible determinar si la luz ha excitado a uno de los modos del sistema.
La intensidad del campo lejano presenta máximos y mínimos relacionado con el radio de la esfera. A mayor tamaño mayor número de máximos. De estas curvas no es posible determinar si la luz ha excitado a uno de los modos del sistema.
En la Fig. 8 se muestra el campo cercano proyectado en el plano z - y para L = 9, radio
Figura 8. La intensidad del campo cercano proyectada en plano z - y por una esfera de radio a = 2πc=ωp , para γ = 0.00 y γ = 0.01 calculado a la frecuencia de resonancia ωa/c = 3.416.
El plasmón de superficie excitado despliega diez lóbulos que rodean a la esfera. El número de estos corresponde adecuadamente al orden multipolar asociado con la onda parcial b5,1.
La extensión del campo fuera de la superficie de la esfera metálica es aproximadamente de a/10, por otro lado se observa que el campo cercano es ligeramente más intenso cuando γ = 0, debido a que no hay absorción.
En las Figs. 8 aparecen dos máximos de intensidad que no hemos podido explicar su naturaleza, solo un análisis detallado de la onda parcial asociada podrá tal vez dar la respuesta al punto. Aunque se intentó encontrar una relación entre el número de oscilaciones o lóbulos dentro de esfera con los parámetros involucrados, como por ejemplo el radio, la longitud de onda, etc., no encontramos alguna relación.
En la Fig. 9 se muestra el campo cercano proyectado en el plano z - y con L = 9, radio
Al igual que en la Fig. 8, en la Fig. 9 se puede observar dos máximos de intensidad, por lo que se puede decir que es un fenómeno que no depende de la resonancia.
4. Conclusiones
Estudiamos teóricamente la excitación de modos electromagnéticos en esferas metálicas. Se asume que tiene un comportamiento tipo Drude 20,21, modelo representativo de los metales, que ha resistido la auscultación por más de un siglo, todavía muy utilizado para el estudio de propiedades ópticas de metales 21. La sección eficaz de esparcimiento total está en excelente acuerdo con el resultado mostrado en la Ref. 12.
La SEE tiene asociado mayor número de excitaciones PS según varia el radio de la esfera metálica.
Se excitó un plasmón de superficie y se observa que los máximos de intensidad del campo eléctrico están muy cerca de la superficie, y son aproximadamente 1000 veces mayores que los que los encontrados cuando no hay excitación. Este comportamiento se observa también en el campo lejano.
Cuando el parámetro γ aumenta, las resonancias en la sección eficaz de esparcimiento total disminuyen su amplitud, e incluso desaparecen y el campo cercano es menos intenso; lo anterior se debe a que la esfera metálica absorbe energía en tales condiciones. Los modos asociados a una esfera metálica son TM únicamente.
Las perspectivas a futuro es el estudio de una esfera cerca de una interfaz plana que divide a dos medios semi-infinitos.
El tema aquí tratado puede ser útil en programas de materias avanzadas de licenciaturas en Física e Ingenierías o posgrados en estos mismos campos.
